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.

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.

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.

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" [...].

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 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.

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.

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 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.

Nucleocapsid (N) Gene Mutations of SARS-CoV-2 Can Affect Real-Time RT-PCR Diagnostic and Impact False-Negative Results.

The current COVID-19 pandemic demands massive testing by Real-time RT-PCR (Reverse Transcription Polymerase Chain Reaction), which is considered the gold standard diagnostic test for the detection of the SARS-CoV-2 virus. However, the virus continues to evolve with mutations that lead to phenotypic alterations as higher transmissibility, pathogenicity or vaccine evasion. Another big issue are mutations in the annealing sites of primers and probes of RT-PCR diagnostic kits leading to false-negative results. Therefore, here we identify mutations in the N (Nucleocapsid) gene that affects the use of the GeneFinder COVID-19 Plus RealAmp Kit. We sequenced SARS-CoV-2 genomes from 17 positive samples with no N gene detection but with RDRP (RNA-dependent RNA polymerase) and E (Envelope) genes detection, and observed a set of three different mutations affecting the N detection: a deletion of 18 nucleotides (Del28877-28894), a substitution of GGG to AAC (28881-28883) and a frameshift mutation caused by deletion (Del28877-28878). The last one cause a deletion of six AAs (amino acids) located in the central intrinsic disorder region at protein level. We also found this mutation in 99 of the 14,346 sequenced samples by the Sao Paulo state Network for Pandemic Alert of Emerging SARS-CoV-2 variants, demonstrating the circulation of the mutation in Sao Paulo, Brazil. Continuous monitoring and characterization of mutations affecting the annealing sites of primers and probes by genomic surveillance programs are necessary to maintain the effectiveness of the diagnosis of COVID-19.

A genome sequence resource for the genus Passiflora, the genome of the wild diploid species Passiflora organensis.

The genus Passiflora comprises a large group of plants popularly known as passionfruit, much appreciated for their exotic flowers and edible fruits. The species (∼500) are morphologically variable (e.g., growth habit, size, and color of flowers) and are adapted to distinct tropical ecosystems. In this study, we generated the genome of the wild diploid species Passiflora organensis Gardner by adopting a hybrid assembly approach. Passiflora organensis has a small genome of 259 Mbp and a heterozygosity rate of 81%, consistent with its reproductive system. Most of the genome sequences could be integrated into its chromosomes with cytogenomic markers (satellite DNA) as references. The repeated sequences accounted for 58.55% of the total DNA analyzed, and the Tekay lineage was the prevalent retrotransposon. In total, 25,327 coding genes were predicted. Passiflora organensis retains 5,609 singletons and 15,671 gene families. We focused on the genes potentially involved in the locus determining self-incompatibility and the MADS-box gene family, allowing us to infer expansions and contractions within specific subfamilies. Finally, we recovered the organellar DNA. Structural rearrangements and two mitoviruses, besides relics of other mobile elements, were found in the chloroplast and mt-DNA molecules, respectively. This study presents the first draft genome assembly of a wild Passiflora species, providing a valuable sequence resource for genomic and evolutionary studies on the genus, and support for breeding cropped passionfruit species.

Analysis of an NGS retinopathy panel detects chromosome 1 uniparental isodisomy in a patient with RPE65-related leber congenital amaurosis.

Purpose: This study aims to demonstrate the possibility of detecting segmental uniparental isodisomy (iUPD) using a next-generation sequencing gene panel by reporting a Leber congenital amaurosis (LCA) case caused by a homozygous pathogenic variant in RPE65 (c.1022 T > C:p.Leu341Ser) inherited exclusively from the proband's mother.Methods: Samples from the trio (proband, mother, and father) were sequenced with a next-generation sequencing (NGS) retinopathy gene panel (224 genes) and the VCF file containing all variants was used in order to determine single nucleotide variant (SNV) counts from each sample across all chromosomes.Results: Trio analysis showed that of 81 Chr1 inherited variants 41 were exclusively maternal, including 21 homozygous. The other 40 variants were common to both parents. On remaining autosomal chromosomes (Chr2-22) 645 inherited variants were found, 147 of them were exclusively maternal and 132 exclusively paternal. Based on these NGS data, it was possible to note that the proband's chromosomes 1 are more similar to his mother's chromosome 1 than his father's, suggesting the pathogenic homozygous variant found in this patient was inherited exclusively from the mother due to uniparental maternal isodisomy.Conclusions: This study presents a secondary analysis pipeline to identify responsible variants for a phenotype and the correct inheritance pattern, which is a critical step to the proper and accurate genetic counseling of all family members. In addition, this approach could be used to determine iUPD in different Mendelian disorders if the sequencing panel identifies variants spread throughout the genome.

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.

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.

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.