Paralog Studies Augment Gene Discovery: DDX and DHX Genes.

Members of a paralogous gene family in which variation in one gene is known to cause disease are eight times more likely to also be associated with human disease. Recent studies have elucidated DHX30 and DDX3X as genes for which pathogenic variant alleles are involved in neurodevelopmental disorders. We hypothesized that variants in paralogous genes encoding members of the DExD/H-box RNA helicase superfamily might also underlie developmental delay and/or intellectual disability (DD and/or ID) disease phenotypes. Here we describe 15 unrelated individuals who have DD and/or ID, central nervous system (CNS) dysfunction, vertebral anomalies, and dysmorphic features and were found to have probably damaging variants in DExD/H-box RNA helicase genes. In addition, these individuals exhibit a variety of other tissue and organ system involvement including ocular, outer ear, hearing, cardiac, and kidney tissues. Five individuals with homozygous (one), compound-heterozygous (two), or de novo (two) missense variants in DHX37 were identified by exome sequencing. We identified ten total individuals with missense variants in three other DDX/DHX paralogs: DHX16 (four individuals), DDX54 (three individuals), and DHX34 (three individuals). Most identified variants are rare, predicted to be damaging, and occur at conserved amino acid residues. Taken together, these 15 individuals implicate the DExD/H-box helicases in both dominantly and recessively inherited neurodevelopmental phenotypes and highlight the potential for more than one disease mechanism underlying these disorders.

A Novel Multisystem Proteinopathy Caused by a Missense ANXA11 Variant.

OBJECTIVE: Protein misfolding plays a central role not only in amyotrophic lateral sclerosis (ALS), but also in other conditions, such as frontotemporal dementia (FTD), inclusion body myopathy (hIBM) or Paget's disease of bone. The concept of multisystem proteinopathies (MSP) was created to account for those rare families that segregate at least 2 out of these 4 conditions in the same pedigree. The calcium-dependent phospholipid-binding protein annexin A11 was recently associated to ALS in European pedigrees. Herein, we describe in detail 3 Brazilian families presenting hIBM (isolated or in combination with ALS/FTD) caused by the novel p.D40Y change in the gene encoding annexin A11 (ANXA11). METHODS: We collected clinical, genetic, pathological and skeletal muscle imaging from 11 affected subjects. Neuroimaging was also obtained from 8 patients and 8 matched controls. RESULTS: Clinico-radiological phenotype of this novel hIBM reveals a slowly progressive predominant limb-girdle syndrome, but with frequent axial (ptosis/dropped head) and distal (medial gastrocnemius) involvement as well. Muscle pathology identified numerous rimmed vacuoles with positive annexin A11, TDP-43 and p62 inclusions, but no inflammation. Central nervous system was also involved: two patients had FTD, but diffusion tensor imaging uncovered multiple areas of cerebral white matter damage in the whole group (including the corticospinal tracts and frontal subcortical regions). INTERPRETATION: These findings expand the phenotypic spectrum related to ANXA11. This gene should be considered the cause of a novel multisystem proteinopathy (MSP type 6), rather than just ALS. ANN NEUROL 2021;90:239-252.

Complete Chromosome-Scale Genome Sequence Resource for Sporisorium panici-leucophaei, the Causal Agent of Sourgrass Smut Disease.

Here, we present the first complete chromosome-level genome assembly of the smut fungus strain Sporisorium panici-leucophaei SPL10A, the causal agent of the sourgrass (Digitaria insularis) smut disease. Combining Illumina paired-end and Nanopore long reads, we generated a final assembly composed of 23 chromosomes (22 nuclear and one mitochondrial) with 18,915,934 bp. Gene prediction accomplished using extrinsic evidence from the sugarcane smut fungus Sporisorium scitamineum originated a total of 6,402 protein-encoding genes. The secretome (388 proteins) and the effectorome repertoires (68 candidates) were also predicted, given their crucial roles in plant-pathogen interactions. The complete telomere-to-telomere chromosome sequences of this poorly studied fungus will provide a valuable resource for future comparative genomic studies among smuts to unravel their underlying pathogenicity mechanisms.[Formula: see text] Copyright © 2021 The Author(s). This is an open access article distributed under the CC BY-NC-ND 4.0 International license.

Biallelic UBE4A loss-of-function variants cause intellectual disability and global developmental delay.

PURPOSE: To identify novel genes associated with intellectual disability (ID) in four unrelated families. METHODS: Here, through exome sequencing and international collaboration, we report eight individuals from four unrelated families of diverse geographic origin with biallelic loss-of-function variants in UBE4A. RESULTS: Eight evaluated individuals presented with syndromic intellectual disability and global developmental delay. Other clinical features included hypotonia, short stature, seizures, and behavior disorder. Characteristic features were appreciated in some individuals but not all; in some cases, features became more apparent with age. We demonstrated that UBE4A loss-of-function variants reduced RNA expression and protein levels in clinical samples. Mice generated to mimic patient-specific Ube4a loss-of-function variant exhibited muscular and neurological/behavioral abnormalities, some of which are suggestive of the clinical abnormalities seen in the affected individuals. CONCLUSION: These data indicate that biallelic loss-of-function variants in UBE4A cause a novel intellectual disability syndrome, suggesting that UBE4A enzyme activity is required for normal development and neurological function.

High-Quality Draft Genome Sequence Resources of Eight Xylella fastidiosa Strains Isolated from Citrus, Coffee, Plum, and Hibiscus in South America.

Xylella fastidiosa subsp. pauca, once confined to South America and infecting mainly citrus and coffee plants, has been found to be associated with other hosts and in other geographic regions. We present high-quality draft genome sequences of X. fastidiosa subsp. pauca strains J1a12, B111, U24D, and XRB isolated from citrus plants in Brazil, strain Fb7 isolated from a citrus plant in Argentina and strains 3124, Pr8x, and Hib4 isolated, respectively, from coffee, plum, and hibiscus plants in Brazil. Sequencing was performed using Roche 454-GS FLX, MiSeq-Illumina or Pacific Biosciences platforms. These high-quality genome assemblies will be useful for further studies about the genomic diversity, evolution, and biology of X. fastidiosa.

A defect in the RNA-processing protein HNRPDL causes limb-girdle muscular dystrophy 1G (LGMD1G).

Limb-girdle muscular dystrophies (LGMD) are a heterogeneous group of genetically determined muscle disorders with a primary or predominant involvement of the pelvic or shoulder girdle musculature. More than 20 genes with autosomal recessive (LGMD2A to LGMD2Q) and autosomal dominant inheritance (LGMD1A to LGMD1H) have been mapped/identified to date. Mutations are known for six among the eight mapped autosomal dominant forms: LGMD1A (myotilin), LGMD1B (lamin A/C), LGMD1C (caveolin-3), LGMD1D (desmin), LGMD1E (DNAJB6), and more recently for LGMD1F (transportin-3). Our group previously mapped the LGMD1G gene at 4q21 in a Caucasian-Brazilian family. We now mapped a Uruguayan family with patients displaying a similar LGMD1G phenotype at the same locus. Whole genome sequencing identified, in both families, mutations in the HNRPDL gene. HNRPDL is a heterogeneous ribonucleoprotein family member, which participates in mRNA biogenesis and metabolism. Functional studies performed in S. cerevisiae showed that the loss of HRP1 (yeast orthologue) had pronounced effects on both protein levels and cell localizations, and yeast proteome revealed dramatic reorganization of proteins involved in RNA-processing pathways. In vivo analysis showed that hnrpdl is important for muscle development in zebrafish, causing a myopathic phenotype when knocked down. The present study presents a novel association between a muscular disorder and a RNA-related gene and reinforces the importance of RNA binding/processing proteins in muscle development and muscle disease. Understanding the role of these proteins in muscle might open new therapeutic approaches for muscular dystrophies.

Homozygous missense mutation in MED25 segregates with syndromic intellectual disability in a large consanguineous family.

BACKGROUND: Intellectual disability (ID) is a highly heterogeneous condition affecting 2% of the population worldwide. In a field study conducted in a highly inbred area of Northeastern Brazil, we investigated a consanguineous family in which seven adults presented syndromic ID. METHODS: Genome-Wide Human SNP Array 6.0 (Affymetrix) microarray was used to determine regions of homozygosity-by-descent and whole exome sequencing (WES) was performed in one affected individual using Extended Nextera Rapid-Capture Exome and Illumina HiSeq2500. RESULTS: We found two regions with an logarithm of the odds (LOD) score of 3.234: a region spanning 4.0 Mb in 19q13.32-q13.33 and a pericentromeric 20 Mb area in chromosome 2 (2p12-q11.2). WES disclosed in the critical region of chromosome 19 a homozygous variant (c.418C>T, p.Arg140Trp) in Mediator complex subunit 25 (MED25), predicted as deleterious by PolyPhen-2, Provean, Mutation Taster and Sorting Intolerant From Tolerant (SIFT). MED25 is a component of the Mediator complex, involved in regulation of transcription of nearly all RNA polymerase II-dependent genes. Deleterious mutations in MED12, MED17 and MED23 have already been associated with ID. CONCLUSIONS: These findings demonstrate that the combination of field investigation of families in highly inbred regions with modern technologies is an effective way for identifying new genes associated with ID.

Building the sugarcane genome for biotechnology and identifying evolutionary trends.

BACKGROUND: Sugarcane is the source of sugar in all tropical and subtropical countries and is becoming increasingly important for bio-based fuels. However, its large (10 Gb), polyploid, complex genome has hindered genome based breeding efforts. Here we release the largest and most diverse set of sugarcane genome sequences to date, as part of an on-going initiative to provide a sugarcane genomic information resource, with the ultimate goal of producing a gold standard genome. RESULTS: Three hundred and seventeen chiefly euchromatic BACs were sequenced. A reference set of one thousand four hundred manually-annotated protein-coding genes was generated. A small RNA collection and a RNA-seq library were used to explore expression patterns and the sRNA landscape. In the sucrose and starch metabolism pathway, 16 non-redundant enzyme-encoding genes were identified. One of the sucrose pathway genes, sucrose-6-phosphate phosphohydrolase, is duplicated in sugarcane and sorghum, but not in rice and maize. A diversity analysis of the s6pp duplication region revealed haplotype-structured sequence composition. Examination of hom(e)ologous loci indicate both sequence structural and sRNA landscape variation. A synteny analysis shows that the sugarcane genome has expanded relative to the sorghum genome, largely due to the presence of transposable elements and uncharacterized intergenic and intronic sequences. CONCLUSION: This release of sugarcane genomic sequences will advance our understanding of sugarcane genetics and contribute to the development of molecular tools for breeding purposes and gene discovery.

Towards accurate indel calling for oncopanel sequencing through an international pipeline competition at precisionFDA.

Accurately calling indels with next-generation sequencing (NGS) data is critical for clinical application. The precisionFDA team collaborated with the U.S. Food and Drug Administration's (FDA's) National Center for Toxicological Research (NCTR) and successfully completed the NCTR Indel Calling from Oncopanel Sequencing Data Challenge, to evaluate the performance of indel calling pipelines. Top performers were selected based on precision, recall, and F1-score. The performance of many other pipelines was close to the top performers, which produced a top cluster of performers. The performance was significantly higher in high confidence regions and coding regions, and significantly lower in low complexity regions. Oncopanel capture and other issues may have occurred that affected the recall rate. Indels with higher variant allele frequency (VAF) may generally be called with higher confidence. Many of the indel calling pipelines had good performance. Some of them performed generally well across all three oncopanels, while others were better for a specific oncopanel. The performance of indel calling can further be improved by restricting the calls within high confidence intervals (HCIs) and coding regions, and by excluding low complexity regions (LCR) regions. Certain VAF cut-offs could be applied according to the applications.

Generation of four induced pluripotent stem cells lines from PBMC of the DFNA58 family members: Two hearing-impaired duplication carriers (USPi006-A e USPi007-A) and two normal-hearing noncarriers (USPi004-A and USPi005-A).

The DFNA58 locus contains a genomic duplication involving three protein-coding genes (CNRIP1, PLEK, and PPP3R1's exon 1) and other uncharacterized lncRNA genes (LOC101927723, LOC107985892 and LOC102724389). To clarify the role of these genes in hearing and precisely determine their role in hearing loss, four iPSC lines were generated from two carriers and two noncarriers of the duplication.

Exploring the Schistosoma mansoni adult male transcriptome using RNA-seq.

Schistosoma mansoni is one of the agents of schistosomiasis, a chronic and debilitating disease. Here we present a transcriptome-wide characterization of adult S. mansoni males by high-throughput RNA-sequencing. We obtained 1,620,432 high-quality ESTs from a directional strand-specific cDNA library, resulting in a 26% higher coverage of genome bases than that of the public ESTs available at NCBI. With a 15×-deep coverage of transcribed genomic regions, our data were able to (i) confirm for the first time 990 predictions without previous evidence of transcription; (ii) correct gene predictions; (iii) discover 989 and 1196 RNA-seq contigs that map to intergenic and intronic genomic regions, respectively, where no gene had been predicted before. These contigs could represent new protein-coding genes or non-coding RNAs (ncRNAs). Interestingly, we identified 11 novel Micro-exon genes (MEGs). These data reveal new features of the S. mansoni transcriptional landscape and significantly advance our understanding of the parasite transcriptome.

New Insights into the Identity of the DFNA58 Gene.

Hearing loss is the most common sensory deficit, affecting 466 million people worldwide. The vast and diverse genes involved reflect the complexity of auditory physiology, which requires the use of animal models in order to gain a fuller understanding. Among the loci with a yet-to-be validated gene is the DFNA58, in which ~200 Kb genomic duplication, including three protein-coding genes (PLEK, CNRIP1, and PPP3R1's exon1), was found to segregate with autosomal dominant hearing loss. Through whole genome sequencing, the duplication was found to be in tandem and inserted in an intergenic region, without the disruption of the topological domains. Reanalysis of transcriptomes data studies (zebrafish and mouse), and RT-qPCR analysis of adult zebrafish target organs, in order to access their orthologues expression, highlighted promising results with Cnrip1a, corroborated by zebrafish in situ hybridization and immunofluorescence. Mouse data also suggested Cnrip1 as the best candidate for a relevant role in auditory physiology, and its importance in hearing seems to have remained conserved but the cell type exerting its function might have changed, from hair cells to spiral ganglion neurons.

Detection of germline variants in Brazilian breast cancer patients using multigene panel testing.

Genetic diversity of germline variants in breast cancer (BC) predisposition genes is unexplored in miscegenated populations, such those living in Latin America. We evaluated 1663 Brazilian BC patients, who underwent hereditary multigene panel testing (20-38 cancer susceptibility genes), to determine the spectrum and prevalence of pathogenic/likely pathogenic (P/LP) variants and variants of uncertain significance (VUS). Associations between P/LP variants and BC risk were estimated in a case-control analysis of BC patients and 18,919 Brazilian reference controls (RC). In total, 335 (20.1%) participants carried germline P/LP variants: 167 (10.0%) in BRCA1/2, 122 (7.3%) in BC actionable non-BRCA genes and 47 (2.8%) in candidate genes or other cancer predisposition genes. Overall, 354 distinctive P/LP variants were identified in 23 genes. The most commonly mutated genes were: BRCA1 (27.4%), BRCA2 (20.3%), TP53 (10.5%), monoallelic MUTYH (9.9%), ATM (8.8%), CHEK2 (6.2%) and PALB2 (5.1%). The Brazilian variant TP53 R337H (c.1010G>A, p.Arg337His), detected in 1.6% of BC patients and 0.1% of RC, was strongly associated with risk of BC, OR = 17.4 (95% CI: 9.4-32.1; p < 0.0001); monoallelic MUTYH variants c.1187G>A and c.536A>G, detected in 1.2% (0.9% RC) and 0.8% (0.4% RC) of the patients, respectively, were not associated with the odds of BC, the former with OR = 1.4 (95% CI: 0.8-2.4; p = 0.29) and the latter with OR = 1.9 (95% CI: 0.9-3.9; p = 0.09). The overall VUS rate was 46.1% for the entire patient population. Concluding, the use of multigene panel testing almost doubled the identification of germline P/LP variants in clinically actionable predisposition genes in BC patients. In Brazil, special attention should be given to TP53 P/LP variants.

Chromosomal microarray analyses from 5778 patients with neurodevelopmental disorders and congenital anomalies in Brazil.

Chromosomal microarray analysis (CMA) has been recommended and practiced routinely since 2010 both in the USA and Europe as the first-tier cytogenetic test for patients with unexplained neurodevelopmental delay/intellectual disability, autism spectrum disorders, and/or multiple congenital anomalies. However, in Brazil, the use of CMA is still limited, due to its high cost and complexity in integrating the results from both the private and public health systems. Although Brazil has one of the world's largest single-payer public healthcare systems, nearly all patients referred for CMA come from the private sector, resulting in only a small number of CMA studies in Brazilian cohorts. To date, this study is by far the largest Brazilian cohort (n = 5788) studied by CMA and is derived from a joint collaboration formed by the University of São Paulo and three private genetic diagnostic centers to investigate the genetic bases of neurodevelopmental disorders and congenital abnormalities. We identified 2,279 clinically relevant CNVs in 1886 patients, not including the 26 cases of UPD found. Among detected CNVs, the corresponding frequency of each category was 55.6% Pathogenic, 4.4% Likely Pathogenic and 40% VUS. The diagnostic yield, by taking into account Pathogenic, Likely Pathogenic and UPDs, was 19.7%. Since the rational for the classification is mostly based on Mendelian or highly penetrant variants, it was not surprising that a second event was detected in 26% of those cases of predisposition syndromes. Although it is common practice to investigate the inheritance of VUS in most laboratories around the world to determine the inheritance of the variant, our results indicate an extremely low cost-benefit of this approach, and strongly suggest that in cases of a limited budget, investigation of the parents of VUS carriers using CMA should not be prioritized.

Genomic epidemiology of the SARS-CoV-2 epidemic in Brazil.

The high numbers of COVID-19 cases and deaths in Brazil have made Latin America an epicentre of the pandemic. SARS-CoV-2 established sustained transmission in Brazil early in the pandemic, but important gaps remain in our understanding of virus transmission dynamics at a national scale. We use 17,135 near-complete genomes sampled from 27 Brazilian states and bordering country Paraguay. From March to November 2020, we detected co-circulation of multiple viral lineages that were linked to multiple importations (predominantly from Europe). After November 2020, we detected large, local transmission clusters within the country. In the absence of effective restriction measures, the epidemic progressed, and in January 2021 there was emergence and onward spread, both within and abroad, of variants of concern and variants under monitoring, including Gamma (P.1) and Zeta (P.2). We also characterized a genomic overview of the epidemic in Paraguay and detected evidence of importation of SARS-CoV-2 ancestor lineages and variants of concern from Brazil. Our findings show that genomic surveillance in Brazil enabled assessment of the real-time spread of emerging SARS-CoV-2 variants.

Genomic monitoring of the SARS-CoV-2 B1.1.7 (WHO VOC Alpha) in the Sao Paulo state, Brazil.

The SARS-CoV-2 alpha VOC (also known as lineage B.1.1.7) initially described in the autumn, 2020 in UK, rapidly became the dominant lineage across much of Europe. Despite multiple studies reporting molecular evidence suggestive of its circulation in Brazil, much is still unknown about its genomic diversity in the state of São Paulo, the main Brazilian economic and transportation hub. To get more insight regarding its transmission dynamics into the State we performed phylogenetic analysis on all alpha VOC strains obtained between February and August 2021 from the Sao Paulo state Network for Pandemic Alert of Emerging SARS-CoV-2 variants. The performed phylogenetic analysis showed that most of the alpha VOC genomes were interspersed with viral strains sampled from different Brazilian states and other countries suggesting that multiple independent Alpha VOC introductions from Brazil and overseas have occurred in the São Paulo State over time. Nevertheless, large monophyletic clusters were also observed especially from the Central-West part of the São Paulo State (the city of Bauru) and the metropolitan region of the São Paulo city. Our results highlight the Alpha VOC molecular epidemiology in the São Paulo state and reinforce the need for continued genomic surveillance strategies for the real-time monitoring of potential emerging SARS-CoV-2 variants during the ever-growing vaccination process.

Genomic epidemiology reveals the impact of national and international restrictions measures on the SARS-CoV-2 epidemic in Brazil.

Brazil has experienced some of the highest numbers of COVID-19 cases and deaths globally and from May 2021 made Latin America a pandemic epicenter. Although SARS-CoV-2 established sustained transmission in Brazil early in the pandemic, important gaps remain in our understanding of virus transmission dynamics at the national scale. Here, we describe the genomic epidemiology of SARS-CoV-2 using near-full genomes sampled from 27 Brazilian states and a bordering country - Paraguay. We show that the early stage of the pandemic in Brazil was characterised by the co-circulation of multiple viral lineages, linked to multiple importations predominantly from Europe, and subsequently characterized by large local transmission clusters. As the epidemic progressed under an absence of effective restriction measures, there was a local emergence and onward international spread of Variants of Concern (VOC) and Variants Under Monitoring (VUM), including Gamma (P.1) and Zeta (P.2). In addition, we provide a preliminary genomic overview of the epidemic in Paraguay, showing evidence of importation from Brazil. These data reinforce the usefulness and need for the implementation of widespread genomic surveillance in South America as a toolkit for pandemic monitoring that provides a means to follow the real-time spread of emerging SARS-CoV-2 variants with possible implications for public health and immunization strategies.

Correction: Lesbon et al. Nucleocapsid (N) Gene Mutations of SARS-CoV-2 Can Affect Real-Time RT-PCR Diagnostic and Impact False-Negative Results. Viruses 2021, 13, 2474.

The authors hereby request the inclusion of two authors (Olivia Teixeira and Maria Cristina Nonato) in the recently published article in Viruses entitled "Nucleocapsid (N) gene mutations of SARS-CoV-2 can affect real-time RT-PCR diagnostic and impact false-negative results" [...].