Structural variant calling and clinical interpretation in 6224 unsolved rare disease exomes.

Structural variants (SVs), including large deletions, duplications, inversions, translocations, and more complex events have the potential to disrupt gene function resulting in rare disease. Nevertheless, current pipelines and clinical decision support systems for exome sequencing (ES) tend to focus on small alterations such as single nucleotide variants (SNVs) and insertions-deletions shorter than 50 base pairs (indels). Additionally, detection and interpretation of large copy-number variants (CNVs) are frequently performed. However, detection of other types of SVs in ES data is hampered by the difficulty of identifying breakpoints in off-target (intergenic or intronic) regions, which makes robust identification of SVs challenging. In this paper, we demonstrate the utility of SV calling in ES resulting in a diagnostic yield of 0.4% (23 out of 5825 probands) for a large cohort of unsolved patients collected by the Solve-RD consortium. Remarkably, 8 out of 23 pathogenic SV were not found by comprehensive read-depth-based CNV analysis, resulting in a 0.13% increased diagnostic value.

Unravelling undiagnosed rare disease cases by HiFi long-read genome sequencing.

Solve-RD is a pan-European rare disease (RD) research program that aims to identify disease-causing genetic variants in previously undiagnosed RD families. We utilised 10-fold coverage HiFi long-read sequencing (LRS) for detecting causative structural variants (SVs), single nucleotide variants (SNVs), insertion-deletions (InDels), and short tandem repeat (STR) expansions in extensively studied RD families without clear molecular diagnoses. Our cohort includes 293 individuals from 114 genetically undiagnosed RD families selected by European Rare Disease Network (ERN) experts. Of these, 21 families were affected by so-called 'unsolvable' syndromes for which genetic causes remain unknown, and 93 families with at least one individual affected by a rare neurological, neuromuscular, or epilepsy disorder without genetic diagnosis despite extensive prior testing. Clinical interpretation and orthogonal validation of variants in known disease genes yielded thirteen novel genetic diagnoses due to de novo and rare inherited SNVs, InDels, SVs, and STR expansions. In an additional four families, we identified a candidate disease-causing SV affecting several genes including an MCF2 / FGF13 fusion and PSMA3 deletion. However, no common genetic cause was identified in any of the 'unsolvable' syndromes. Taken together, we found (likely) disease-causing genetic variants in 13.0% of previously unsolved families and additional candidate disease-causing SVs in another 4.3% of these families. In conclusion, our results demonstrate the added value of HiFi long-read genome sequencing in undiagnosed rare diseases.

A GGC-repeat expansion in ZFHX3 encoding polyglycine causes spinocerebellar ataxia type 4 and impairs autophagy.

Despite linkage to chromosome 16q in 1996, the mutation causing spinocerebellar ataxia type 4 (SCA4), a late-onset sensory and cerebellar ataxia, remained unknown. Here, using long-read single-strand whole-genome sequencing (LR-GS), we identified a heterozygous GGC-repeat expansion in a large Utah pedigree encoding polyglycine (polyG) in zinc finger homeobox protein 3 (ZFHX3), also known as AT-binding transcription factor 1 (ATBF1). We queried 6,495 genome sequencing datasets and identified the repeat expansion in seven additional pedigrees. Ultrarare DNA variants near the repeat expansion indicate a common distant founder event in Sweden. Intranuclear ZFHX3-p62-ubiquitin aggregates were abundant in SCA4 basis pontis neurons. In fibroblasts and induced pluripotent stem cells, the GGC expansion led to increased ZFHX3 protein levels and abnormal autophagy, which were normalized with small interfering RNA-mediated ZFHX3 knockdown in both cell types. Improving autophagy points to a therapeutic avenue for this novel polyG disease. The coding GGC-repeat expansion in an extremely G+C-rich region was not detectable by short-read whole-exome sequencing, which demonstrates the power of LR-GS for variant discovery.

Genomes in clinical care.

In the era of precision medicine, genome sequencing (GS) has become more affordable and the importance of genomics and multi-omics in clinical care is increasingly being recognized. However, how to scale and effectively implement GS on an institutional level remains a challenge for many. Here, we present Genome First and Ge-Med, two clinical implementation studies focused on identifying the key pillars and processes that are required to make routine GS and predictive genomics a reality in the clinical setting. We describe our experience and lessons learned for a variety of topics including test logistics, patient care processes, data reporting, and infrastructure. Our model of providing clinical care and comprehensive genomic analysis from a single source may be used by other centers with a similar structure to facilitate the implementation of omics-based personalized health concepts in medicine.

ZSCAN10 deficiency causes a neurodevelopmental disorder with characteristic oto-facial malformations.

Neurodevelopmental disorders are major indications for genetic referral and have been linked to more than 1,500 loci including genes encoding transcriptional regulators. The dysfunction of transcription factors often results in characteristic syndromic presentations, however, at least half of these patients lack a genetic diagnosis. The implementation of machine learning approaches has the potential to aid in the identification of new disease genes and delineate associated phenotypes. Next generation sequencing was performed in seven affected individuals with neurodevelopmental delay and dysmorphic features. Clinical characterization included reanalysis of available neuroimaging datasets and 2D portrait image analysis with GestaltMatcher. The functional consequences of ZSCAN10 loss were modelled in mouse embryonic stem cells (mESC), including a knock-out and a representative ZSCAN10 protein truncating variant. These models were characterized by gene expression and Western blot analyses, chromatin immunoprecipitation and quantitative PCR (ChIP-qPCR), and immunofluorescence staining. Zscan10 knockout mouse embryos were generated and phenotyped. We prioritized bi-allelic ZSCAN10 loss-of-function variants in seven affected individuals from five unrelated families as the underlying molecular cause. RNA-Seq analyses in Zscan10-/- mESCs indicated dysregulation of genes related to stem cell pluripotency. In addition, we established in mESCs the loss-of-function mechanism for a representative human ZSCAN10 protein truncating variant by showing alteration of its expression levels and subcellular localization, interfering with its binding to DNA enhancer targets. Deep phenotyping revealed global developmental delay, facial asymmetry, and malformations of the outer ear as consistent clinical features. Cerebral MRI showed dysplasia of the semicircular canals as an anatomical correlate of sensorineural hearing loss. Facial asymmetry was confirmed as a clinical feature by GestaltMatcher and was recapitulated in the Zscan10 mouse model along with inner and outer ear malformations. Our findings provide evidence of a novel syndromic neurodevelopmental disorder caused by bi-allelic loss-of-function variants in ZSCAN10.

Lessons Learned from Translating Genome Sequencing to Clinical Routine: Understanding the Accuracy of a Diagnostic Pipeline.

The potential of genome sequencing (GS), which allows detection of almost all types of genetic variation across nearly the entire genome of an individual, greatly expands the possibility for diagnosing genetic disorders. The opportunities provided with this single test are enticing to researchers and clinicians worldwide for human genetic research as well as clinical application. Multiple studies have highlighted the advantages of GS for genetic variant discovery, emphasizing its added value for routine clinical use. We have implemented GS as first-line genetic testing for patients with rare diseases. Here, we report on our experiences in establishing GS as a reliable diagnostic method for almost all types of genetic disorders, from validating diagnostic accuracy of sequencing pipelines to clinical implementation in routine practice.

Mobile element insertions in rare diseases: a comparative benchmark and reanalysis of 60,000 exome samples.

Mobile element insertions (MEIs) are a known cause of genetic disease but have been underexplored due to technical limitations of genetic testing methods. Various bioinformatic tools have been developed to identify MEIs in Next Generation Sequencing data. However, most tools have been developed specifically for genome sequencing (GS) data rather than exome sequencing (ES) data, which remains more widely used for routine diagnostic testing. In this study, we benchmarked six MEI detection tools (ERVcaller, MELT, Mobster, SCRAMble, TEMP2 and xTea) on ES data and on GS data from publicly available genomic samples (HG002, NA12878). For all the tools we evaluated sensitivity and precision of different filtering strategies. Results show that there were substantial differences in tool performance between ES and GS data. MELT performed best with ES data and its combination with SCRAMble increased substantially the detection rate of MEIs. By applying both tools to 10,890 ES samples from Solve-RD and 52,624 samples from Radboudumc we were able to diagnose 10 patients who had remained undiagnosed by conventional ES analysis until now. Our study shows that MELT and SCRAMble can be used reliably to identify clinically relevant MEIs in ES data. This may lead to an additional diagnosis for 1 in 3000 to 4000 patients in routine clinical ES.

Diagnostic genome sequencing improves diagnostic yield: a prospective single-centre study in 1000 patients with inherited eye diseases.

PURPOSE: Genome sequencing (GS) is expected to reduce the diagnostic gap in rare disease genetics. We aimed to evaluate a scalable framework for genome-based analyses 'beyond the exome' in regular care of patients with inherited retinal degeneration (IRD) or inherited optic neuropathy (ION). METHODS: PCR-free short-read GS was performed on 1000 consecutive probands with IRD/ION in routine diagnostics. Complementary whole-blood RNA-sequencing (RNA-seq) was done in a subset of 74 patients. An open-source bioinformatics analysis pipeline was optimised for structural variant (SV) calling and combined RNA/DNA variation interpretation. RESULTS: A definite genetic diagnosis was established in 57.4% of cases. For another 16.7%, variants of uncertain significance were identified in known IRD/ION genes, while the underlying genetic cause remained unresolved in 25.9%. SVs or alterations in non-coding genomic regions made up for 12.7% of the observed variants. The RNA-seq studies supported the classification of two unclear variants. CONCLUSION: GS is feasible in clinical practice and reliably identifies causal variants in a substantial proportion of individuals. GS extends the diagnostic yield to rare non-coding variants and enables precise determination of SVs. The added diagnostic value of RNA-seq is limited by low expression levels of the major IRD disease genes in blood.

Genetic and functional analysis of chymotrypsin-like protease (CTRL) in chronic pancreatitis.

BACKGROUND: Genetic predisposition is crucial in the pathogenesis of early-onset chronic pancreatitis (CP). So far, several genetic alterations have been identified as risk factors, predominantly in genes encoding digestive enzymes. However, many early-onset CP cases have no identified underlying cause. Chymotrypsins are a family of serine proteases that can cleave trypsinogen and lead to its degradation. Because genetic alterations in the chymotrypsins CTRC, CTRB1, and CTRB2 are associated with CP, we genetically and functionally investigated chymotrypsin-like protease (CTRL) as a potential risk factor. METHODS: We screened 1005 non-alcoholic CP patients and 1594 controls for CTRL variants by exome sequencing. We performed Western blots and activity assays to analyse secretion and proteolytic activity. We measured BiP mRNA expression to investigate the potential impact of identified alterations on endoplasmic reticulum (ER) stress. RESULTS: We identified 13 heterozygous non-synonymous CTRL variants: five exclusively in patients and three only in controls. Functionality was unchanged in 6/13 variants. Four alterations showed normal secretion but reduced (p.G20S, p.G56S, p.G61S) or abolished (p.S208F) activity. Another three variants (p.C201Y, p.G215R and p.C220G) were not secreted and already showed reduced or no activity intracellularly. However, intracellular retention did not lead to ER stress. CONCLUSION: We identified several CTRL variants, some showing potent effects on protease function and secretion. We observed these effects in variants found in patients and controls, and CTRL loss-of-function variants were not significantly more common in patients than controls. Therefore, CTRL is unlikely to play a relevant role in the development of CP.

Biallelic MAD2L1BP (p31comet) mutation is associated with mosaic aneuploidy and juvenile granulosa cell tumors.

MAD2L1BP-encoded p31comet mediates Trip13-dependent disassembly of Mad2- and Rev7-containing complexes and, through this antagonism, promotes timely spindle assembly checkpoint (SAC) silencing, faithful chromosome segregation, insulin signaling, and homology-directed repair (HDR) of DNA double-strand breaks. We identified a homozygous MAD2L1BP nonsense variant, R253*, in 2 siblings with microcephaly, epileptic encephalopathy, and juvenile granulosa cell tumors of ovary and testis. Patient-derived cells exhibited high-grade mosaic variegated aneuploidy, slowed-down proliferation, and instability of truncated p31comet mRNA and protein. Corresponding recombinant p31comet was defective in Trip13, Mad2, and Rev7 binding and unable to support SAC silencing or HDR. Furthermore, C-terminal truncation abrogated an identified interaction of p31comet with tp53. Another homozygous truncation, R227*, detected in an early-deceased patient with low-level aneuploidy, severe epileptic encephalopathy, and frequent blood glucose elevations, likely corresponds to complete loss of function, as in Mad2l1bp-/- mice. Thus, human mutations of p31comet are linked to aneuploidy and tumor predisposition.

Multicentric pilot study to standardize clinical whole exome sequencing (WES) for cancer patients.

A growing number of druggable targets and national initiatives for precision oncology necessitate broad genomic profiling for many cancer patients. Whole exome sequencing (WES) offers unbiased analysis of the entire coding sequence, segmentation-based detection of copy number alterations (CNAs), and accurate determination of complex biomarkers including tumor mutational burden (TMB), homologous recombination repair deficiency (HRD), and microsatellite instability (MSI). To assess the inter-institution variability of clinical WES, we performed a comparative pilot study between German Centers of Personalized Medicine (ZPMs) from five participating institutions. Tumor and matched normal DNA from 30 patients were analyzed using custom sequencing protocols and bioinformatic pipelines. Calling of somatic variants was highly concordant with a positive percentage agreement (PPA) between 91 and 95% and a positive predictive value (PPV) between 82 and 95% compared with a three-institution consensus and full agreement for 16 of 17 druggable targets. Explanations for deviations included low VAF or coverage, differing annotations, and different filter protocols. CNAs showed overall agreement in 76% for the genomic sequence with high wet-lab variability. Complex biomarkers correlated strongly between institutions (HRD: 0.79-1, TMB: 0.97-0.99) and all institutions agreed on microsatellite instability. This study will contribute to the development of quality control frameworks for comprehensive genomic profiling and sheds light onto parameters that require stringent standardization.

UV-radiation and MC1R germline mutations are risk factors for the development of conventional and spitzoid melanomas in children and adolescents.

BACKGROUND: Genomic characterisation has led to an improved understanding of adult melanoma. However, the aetiology of melanoma in children is still unclear and identifying the correct diagnosis and therapeutic strategies remains challenging. METHODS: Exome sequencing of matched tumour-normal pairs from 26 paediatric patients was performed to study the mutational spectrum of melanomas. The cohort was grouped into different categories: spitzoid melanoma (SM), conventional melanoma (CM), and other melanomas (OT). FINDINGS: In all patients with CM (n = 10) germline variants associated with melanoma were found in low to moderate melanoma risk genes: in 8 patients MC1R variants, in 2 patients variants in MITF, PTEN and BRCA2. Somatic BRAF mutations were detected in 60% of CMs, homozygous deletions of CDKN2A in 20%, TERTp mutations in 30%. In the SM group (n = 12), 5 patients carried at least one MC1R variant; somatic BRAF mutations were detected in 8.3%, fusions in 25% of the cases. No SM showed a homozygous CDKN2A deletion nor a TERTp mutation. In 81.8% of the CM/SM cases the UV damage signatures SBS7 and/or DBS1 were detected. The patient with melanoma arising in giant congenital nevus (CNM) demonstrated the characteristic NRAS Q61K mutation. INTERPRETATION: UV-radiation and MC1R germline variants are risk factors in the development of conventional and spitzoid paediatric melanomas. Paediatric CMs share genomic similarities with adult CMs while the SMs differ genetically from the CM group. Consistent genetic characterization of all paediatric melanomas will potentially lead to better subtype differentiation, treatment, and prevention in the future. FUNDING: Found in Acknowledgement.

Clinical trio genome sequencing facilitates the interpretation of variants in cancer predisposition genes in paediatric tumour patients.

The prevalence of pathogenic and likely pathogenic (P/LP) variants in genes associated with cancer predisposition syndromes (CPS) is estimated to be 8-18% for paediatric cancer patients. In more than half of the carriers, the family history is unsuspicious for CPS. Therefore, broad genetic testing could identify germline predisposition in additional children with cancer resulting in important implications for themselves and their families. We thus evaluated clinical trio genome sequencing (TGS) in a cohort of 72 paediatric patients with solid cancers other than retinoblastoma or CNS-tumours. The most prevalent cancer types were sarcoma (n = 26), neuroblastoma (n = 15), and nephroblastoma (n = 10). Overall, P/LP variants in CPS genes were identified in 18.1% of patients (13/72) and P/LP variants in autosomal-dominant CPS genes in 9.7% (7/72). Genetic evaluation would have been recommended for the majority of patients with P/LP variants according to the Jongmans criteria. Four patients (5.6%, 4/72) carried P/LP variants in autosomal-dominant genes known to be associated with their tumour type. With the immediate information on variant inheritance, TGS facilitated the identification of a de novo P/LP in NF1, a gonadosomatic mosaic in WT1 and two pathogenic variants in one patient (DICER1 and PALB2). TGS allows a more detailed characterization of structural variants with base-pair resolution of breakpoints which can be relevant for the interpretation of copy number variants. Altogether, TGS allows comprehensive identification of children with a CPS and supports the individualised clinical management of index patients and high-risk relatives.

First case of adult onset neuronal intranuclear inclusion disease with both typical radiological signs and NOTCH2NLC repeat expansions in a Caucasian individual.

BACKGROUND AND PURPOSE: Adult onset neuronal intranuclear inclusion disease (NIID) is a rare neurodegenerative disorder with a heterogeneous clinical presentation that can mimic stroke and various forms of dementia. To date, it has been described almost exclusively in Asian individuals. METHODS: This case presentation includes magnetic resonance imaging (MRI) of the neurocranium, histology by skin biopsy, and long-read genome sequencing. RESULTS: A 75-year-old Caucasian female presented with paroxysmal encephalopathy twice within a 14-month period. Brain MRI revealed high-intensity signals at the cerebral corticomedullary junction (diffusion-weighted imaging) and the paravermal area (fluid-attenuated inversion recovery), a typical distribution observed in adult onset NIID. The diagnosis was corroborated by skin biopsy, which demonstrated eosinophilic intranuclear inclusion bodies, and confirmed by long-read genome sequencing, showing an expansion of the GGC repeat in exon 1 of NOTCH2NLC. CONCLUSIONS: Our case proves adult onset NOTCH2NLC-GGC-positive NIID with typical findings on MRI and histology in a Caucasian patient and underscores the need to consider this diagnosis in non-Asian individuals.

Whole genome analysis of two sympatric human Mansonella: Mansonella perstans and Mansonella sp "DEUX".

Introduction: Mansonella species are filarial parasites that infect humans worldwide. Although these infections are common, knowledge of the pathology and diversity of the causative species is limited. Furthermore, the lack of sequencing data for Mansonella species, shows that their research is neglected. Apart from Mansonella perstans, a potential new species called Mansonella sp "DEUX" has been identified in Gabon, which is prevalent at high frequencies. We aimed to further determine if Mansonella sp "DEUX" is a genotype of M. perstans, or if these are two sympatric species. Methods: We screened individuals in the area of Fougamou, Gabon for Mansonella mono-infections and generated de novo assemblies from the respective samples. For evolutionary analysis, a phylogenetic tree was reconstructed, and the differences and divergence times are presented. In addition, mitogenomes were generated and phylogenies based on 12S rDNA and cox1 were created. Results: We successfully generated whole genomes for M. perstans and Mansonella sp "DEUX". Phylogenetic analysis based on annotated protein sequences, support the hypothesis of two distinct species. The inferred evolutionary analysis suggested, that M. perstans and Mansonella sp "DEUX" separated around 778,000 years ago. Analysis based on mitochondrial marker genes support our hypothesis of two sympatric human Mansonella species. Discussion: The results presented indicate that Mansonella sp "DEUX" is a new Mansonella species. These findings reflect the neglect of this research topic. And the availability of whole genome data will allow further investigations of these species.

Creation of a structured molecular genomics report for Germany as a local adaption of HL7's Genomic Reporting Implementation Guide.

OBJECTIVE: The objective was to develop a dataset definition, information model, and FHIR® specification for key data elements contained in a German molecular genomics (MolGen) report to facilitate genomic and phenotype integration in electronic health records. MATERIALS AND METHODS: A dedicated expert group participating in the German Medical Informatics Initiative reviewed information contained in MolGen reports, determined the key elements, and formulated a dataset definition. HL7's Genomics Reporting Implementation Guide (IG) was adopted as a basis for the FHIR® specification which was subjected to a public ballot. In addition, elements in the MolGen dataset were mapped to the fields defined in ISO/TS 20428:2017 standard to evaluate compliance. RESULTS: A core dataset of 76 data elements, clustered into 6 categories was created to represent all key information of German MolGen reports. Based on this, a FHIR specification with 16 profiles, 14 derived from HL7®'s Genomics Reporting IG and 2 additional profiles (of the FamilyMemberHistory and RiskAssessment resources), was developed. Five example resource bundles show how our adaptation of an international standard can be used to model MolGen report data that was requested following oncological or rare disease indications. Furthermore, the map of the MolGen report data elements to the fields defined by the ISO/TC 20428:2017 standard, confirmed the presence of the majority of required fields. CONCLUSIONS: Our report serves as a template for other research initiatives attempting to create a standard format for unstructured genomic report data. Use of standard formats facilitates integration of genomic data into electronic health records for clinical decision support.

Detection of circulating cell-free HPV DNA of 13 HPV types for patients with cervical cancer as potential biomarker to monitor therapy response and to detect relapse.

BACKGROUND: HPV-related cervical cancer (CC) is the fourth most frequent cancer in women worldwide. Cell-free tumour DNA is a potent biomarker to detect treatment response, residual disease, and relapse. We investigated the potential use of cell-free circulating HPV-DNA (cfHPV-DNA) in plasma of patients with CC. METHODS: cfHPV-DNA levels were measured using a highly sensitive next-generation sequencing-based approach targeting a panel of 13 high-risk HPV types. RESULTS: Sequencing was performed in 69 blood samples collected from 35 patients, of which 26 were treatment-naive when the first liquid biopsy sample was retrieved. cfHPV-DNA was successfully detected in 22/26 (85%) cases. A significant correlation between tumour burden and cfHPV-DNA levels was observed: cfHPV-DNA was detectable in all treatment-naive patients with advanced-stage disease (17/17, FIGO IB3-IVB) and in 5/9 patients with early-stage disease (FIGO IA-IB2). Sequential samples revealed a decrease of cfHPV-DNA levels in 7 patients corresponding treatment response and an increase in a patient with relapse. CONCLUSIONS: In this proof-of-concept study we demonstrated the potential of cfHPV-DNA as a biomarker for therapy monitoring in patients with primary and recurrent CC. Our findings facilitate the development of a sensitive and precise, non-invasive, inexpensive, and easily accessible tool in CC diagnosis, therapy monitoring and follow-up.

Clinical outcome of biomarker-guided therapies in adult patients with tumors of the nervous system.

Background: The clinical utility of molecular profiling and targeted therapies for neuro-oncology patients outside of clinical trials is not established. We aimed at investigating feasibility and clinical utility of molecular profiling and targeted therapy in adult patients with advanced tumors in the nervous system within a prospective observational study. Methods: molecular tumor board (MTB)@ZPM (NCT03503149) is a prospective observational precision medicine study for patients with advanced tumors. After inclusion of patients, we performed comprehensive molecular profiling, formulated ranked biomarker-guided therapy recommendations based on consensus by the MTB, and collected prospective clinical outcome data. Results: Here, we present initial data of 661 adult patients with tumors of the nervous system enrolled by December 31, 2021. Of these, 408 patients were presented at the MTB. Molecular-instructed therapy recommendations could be made in 380/408 (93.1%) cases and were prioritized by evidence levels. Therapies were initiated in 86/380 (22.6%) cases until data cutoff. We observed a progression-free survival ratio >1.3 in 31.3% of patients. Conclusions: Our study supports the clinical utility of biomarker-guided therapies for neuro-oncology patients and indicates clinical benefit in a subset of patients. Our data might inform future clinical trials, translational studies, and even clinical care.