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.

Accepting or declining preconception expanded carrier screening: An exploratory study with 407 couples.

Rapidly evolving genomic technologies have made genetic expanded carrier screening (ECS) possible for couples considering a pregnancy. The aim of ECS is to identify couples at risk of having a child affected with a severe disorder and to facilitate their reproductive decision-making process. The ECS test we offer at our center, called BeGECS (Belgian Genetic ECS), consists of 1268 autosomal recessive (AR) and X-linked pathogenic genes, including severe childhood-onset disorders. However, thus far data are scarce regarding the actual uptake of preconception ECS in a clinical setting. Therefore, our aim was to describe the characteristics of 407 couples to whom ECS was offered at the Center for Medical Genetics of the University Hospital Ghent (CMGG). In addition, we aimed to identify their reasons for accepting or declining BeGECS. Between October 2019 and January 2023, 407 preconception couples were offered BeGECS and were asked to fill in a questionnaire after their decision. Of the 407 couples participating in the survey, 270 (66%) decided to take the test and 137 (34%) declined. We observed that age, highest education level as well as indication for consultation were statistically different between the group that accepted to take the test and the group that declined (p = 0.037). In particular, age and education level were substantially higher in the group that accepted the test. Major reasons for taking BeGECS include prevention, wishing to obtain all information possible, helping preparing their future reproductive decision and increasing their sense of control by being informed. However, couples that do not chose to take BeGECS stated that too much information would make them anxious, that the result would not change their decision to have children, that they do not want to spend money on something that will not happen and that they do not worry about their family history. These findings show that the majority of preconception couples that were offered ECS, accepted the test.

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.

Systematic analysis of paralogous regions in 41,755 exomes uncovers clinically relevant variation.

The short lengths of short-read sequencing reads challenge the analysis of paralogous genomic regions in exome and genome sequencing data. Most genetic variants within these homologous regions therefore remain unidentified in standard analyses. Here, we present a method (Chameleolyser) that accurately identifies single nucleotide variants and small insertions/deletions (SNVs/Indels), copy number variants and ectopic gene conversion events in duplicated genomic regions using whole-exome sequencing data. Application to a cohort of 41,755 exome samples yields 20,432 rare homozygous deletions and 2,529,791 rare SNVs/Indels, of which we show that 338,084 are due to gene conversion events. None of the SNVs/Indels are detectable using regular analysis techniques. Validation by high-fidelity long-read sequencing in 20 samples confirms >88% of called variants. Focusing on variation in known disease genes leads to a direct molecular diagnosis in 25 previously undiagnosed patients. Our method can readily be applied to existing exome data.

G protein-coupled receptor kinase 6 (GRK6) regulates insulin processing and secretion via effects on proinsulin conversion to insulin.

Recent studies identified a missense mutation in the gene coding for G protein-coupled receptor kinase 6 (GRK6) that segregates with type 2 diabetes (T2D). To better understand how GRK6 might be involved in T2D, we used pharmacological inhibition and genetic knockdown in the mouse ß-cell line, MIN6, to determine whether GRK6 regulates insulin dynamics. We show inhibition of GRK5 and GRK6 increased insulin secretion but reduced insulin processing while GRK6 knockdown revealed these same processing defects with reduced levels of cellular insulin. GRK6 knockdown cells also had attenuated insulin secretion but enhanced proinsulin secretion consistent with decreased processing. In support of these findings, we demonstrate GRK6 rescue experiments in knockdown cells restored insulin secretion after glucose treatment. The altered insulin profile appears to be caused by changes in the proprotein convertases, the enzymes responsible for proinsulin to insulin conversion, as GRK6 knockdown resulted in significantly reduced convertase expression and activity. To identify how the GRK6-P384S mutation found in T2D patients might affect insulin processing, we performed biochemical and cell biological assays to study the properties of the mutant. We found that while GRK6-P384S was more active than WT GRK6, it displayed a cytosolic distribution in cells compared to the normal plasma membrane localization of GRK6. Additionally, GRK6 overexpression in MIN6 cells enhanced proinsulin processing, while GRK6-P384S expression had little effect. Taken together, our data show that GRK6 regulates insulin processing and secretion in a glucose-dependent manner and provide a foundation for understanding the contribution of GRK6 to T2D.

Exploring the Mutational Landscape of Isolated Congenital Heart Defects: An Exome Sequencing Study Using Cardiac DNA.

Congenital heart defects (CHD) are the most common congenital anomalies in liveborn children. In contrast to syndromic CHD (SCHD), the genetic basis of isolated CHD (ICHD) is complex, and the underlying pathogenic mechanisms appear intricate and are incompletely understood. Next to rare Mendelian conditions, somatic mosaicism or a complex multifactorial genetic architecture are assumed for most ICHD. We performed exome sequencing (ES) in 73 parent-offspring ICHD trios using proband DNA extracted from cardiac tissue. We identified six germline de novo variants and 625 germline rare inherited variants with 'damaging' in silico predictions in cardiac-relevant genes expressed in the developing human heart. There were no CHD-relevant somatic variants. Transmission disequilibrium testing (TDT) and association testing (AT) yielded no statistically significant results, except for the AT of missense variants in cilia genes. Somatic mutations are not a common cause of ICHD. Rare de novo and inherited protein-damaging variants may contribute to ICHD, possibly as part of an oligogenic or polygenic disease model. TDT and AT failed to provide informative results, likely due to the lack of power, but provided a framework for future studies in larger cohorts. Overall, the diagnostic value of ES on cardiac tissue is limited in individual ICHD cases.

DeNovoCNN: a deep learning approach to de novo variant calling in next generation sequencing data.

De novo mutations (DNMs) are an important cause of genetic disorders. The accurate identification of DNMs from sequencing data is therefore fundamental to rare disease research and diagnostics. Unfortunately, identifying reliable DNMs remains a major challenge due to sequence errors, uneven coverage, and mapping artifacts. Here, we developed a deep convolutional neural network (CNN) DNM caller (DeNovoCNN), that encodes the alignment of sequence reads for a trio as 160$ \times$164 resolution images. DeNovoCNN was trained on DNMs of 5616 whole exome sequencing (WES) trios achieving total 96.74% recall and 96.55% precision on the test dataset. We find that DeNovoCNN has increased recall/sensitivity and precision compared to existing DNM calling approaches (GATK, DeNovoGear, DeepTrio, Samtools) based on the Genome in a Bottle reference dataset and independent WES and WGS trios. Validations of DNMs based on Sanger and PacBio HiFi sequencing confirm that DeNovoCNN outperforms existing methods. Most importantly, our results suggest that DeNovoCNN is likely robust against different exome sequencing and analyses approaches, thereby allowing the application on other datasets. DeNovoCNN is freely available as a Docker container and can be run on existing alignment (BAM/CRAM) and variant calling (VCF) files from WES and WGS without a need for variant recalling.

Hypergastrinemia, a clue leading to the identification of an atypical form of diabetes mellitus type 2.

BACKGROUND: A hitherto undescribed form of diabetes mellitus type 2 is reported in a Flemish family. In these patients, markedly elevated gastrin levels were observed, which could not be linked to gastrointestinal symptoms. MATERIALS AND METHODS: Gel permeation chromatography was performed for gastrin, insulin, and proinsulin. Proprotein convertase subtilisin/kexin type (PCSK1 and PCSK2)] were sequenced. Whole-exome sequencing was performed on the genomic DNA extracted from leukocytes of the proband of the family. RESULTS: Gel permeation chromatography revealed that the apparent hypergastrinemia was caused by the accumulation of biologically inactive progastrin. Besides, high serum concentrations of proinsulin and intact fibroblast growth factor 23 (FGF23) were also detected. Sequencing of PCSK1 and PCSK2 genes did not reveal any mutations in these genes. Whole exome sequencing revealed a c.1150C > T (p.Pro384Ser) mutation in G protein-coupled receptor kinase 6 (GRK6), which cosegregated with the disease. Expression of the mutant enzyme in mammalian cells revealed that it was mislocalized compared to the wild-type GRK6. CONCLUSIONS: In the affected patients, prohormone processing is impaired likely due to the altered function of mutant GRK6. Delayed pro-insulin processing causes hypoglycaemia episodes a couple of hours following meals. In addition, increased plasma concentrations of progastrin and intact FGF23 in the affected individuals can be explained by incomplete processing of the precursor hormones.

Shortcutting the diagnostic odyssey: the multidisciplinary Program for Undiagnosed Rare Diseases in adults (UD-PrOZA).

BACKGROUND: In order to facilitate the diagnostic process for adult patients suffering from a rare disease, the Undiagnosed Disease Program (UD-PrOZA) was founded in 2015 at the Ghent University Hospital in Belgium. In this study we report the five-year results of our multidisciplinary approach in rare disease diagnostics. METHODS: Patients referred by a healthcare provider, in which an underlying rare disease is likely, qualify for a UD-PrOZA evaluation. UD-PrOZA uses a multidisciplinary clinical approach combined with state-of-the-art genomic technologies in close collaboration with research facilities to diagnose patients. RESULTS: Between 2015 and 2020, 692 patients (94% adults) were referred of which 329 (48%) were accepted for evaluation. In 18% (60 of 329) of the cases a definite diagnosis was made. 88% (53 of 60) of the established diagnoses had a genetic origin. 65% (39 of 60) of the genetic diagnoses were made through whole exome sequencing (WES). The mean time interval between symptom-onset and diagnosis was 19 years. Key observations included novel genotype-phenotype correlations, new variants in known disease genes and the identification of three new disease genes. In 13% (7 of 53), identifying the molecular cause was associated with therapeutic recommendations and in 88% (53 of 60), gene specific genetic counseling was made possible. Actionable secondary findings were reported in 7% (12 of 177) of the patients in which WES was performed. CONCLUSION: UD-PrOZA offers an innovative interdisciplinary platform to diagnose rare diseases in adults with previously unexplained medical problems and to facilitate translational research.

A MT-TL1 variant identified by whole exome sequencing in an individual with intellectual disability, epilepsy, and spastic tetraparesis.

The genetic etiology of intellectual disability remains elusive in almost half of all affected individuals. Within the Solve-RD consortium, systematic re-analysis of whole exome sequencing (WES) data from unresolved cases with (syndromic) intellectual disability (n = 1,472 probands) was performed. This re-analysis included variant calling of mitochondrial DNA (mtDNA) variants, although mtDNA is not specifically targeted in WES. We identified a functionally relevant mtDNA variant in MT-TL1 (NC_012920.1:m.3291T > C; NC_012920.1:n.62T > C), at a heteroplasmy level of 22% in whole blood, in a 23-year-old male with severe intellectual disability, epilepsy, episodic headaches with emesis, spastic tetraparesis, brain abnormalities, and feeding difficulties. Targeted validation in blood and urine supported pathogenicity, with heteroplasmy levels of 23% and 58% in index, and 4% and 17% in mother, respectively. Interestingly, not all phenotypic features observed in the index have been previously linked to this MT-TL1 variant, suggesting either broadening of the m.3291T > C-associated phenotype, or presence of a co-occurring disorder. Hence, our case highlights the importance of underappreciated mtDNA variants identifiable from WES data, especially for cases with atypical mitochondrial phenotypes and their relatives in the maternal line.

Comprehensive validation of a diagnostic strategy for sequencing genes with one or multiple pseudogenes using pseudoxanthoma elasticum as a model.

Pseudogenes are frequently encountered noncoding sequences with a high sequence similarity to their protein-coding paralogue. For this reason, their presence is often considered troublesome in molecular diagnostics. In pseudoxanthoma elasticum (PXE), a disease predominantly caused by mutations in ATP-binding cassette family C member 6 (ABCC6), the presence of two pseudogenes complicates the analysis of sequence data. With whole-exome sequencing (WES) becoming the standard of care in molecular diagnostics, we wanted to evaluate whether this technique is as reliable as gene-specific targeted enrichment analysis for the analysis of ABCC6. We established a PCR-based targeted enrichment and next-generation sequencing testing approach and demonstrated that the ABCC6-specific enrichment combined with the applied mapping algorithm overcomes the complication of ABCC6 pseudogene aspecificities, contrary to WES. We propose a time- and cost-efficient diagnostic strategy for comprehensive and accurate molecular genetic testing of PXE, which is highly automatable.

A clinical scoring system for congenital contractural arachnodactyly.

PURPOSE: Congenital contractural arachnodactyly (CCA) is an autosomal dominant connective tissue disorder manifesting joint contractures, arachnodactyly, crumpled ears, and kyphoscoliosis as main features. Due to its rarity, rather aspecific clinical presentation, and overlap with other conditions including Marfan syndrome, the diagnosis is challenging, but important for prognosis and clinical management. CCA is caused by pathogenic variants in FBN2, encoding fibrillin-2, but locus heterogeneity has been suggested. We designed a clinical scoring system and diagnostic criteria to support the diagnostic process and guide molecular genetic testing. METHODS: In this retrospective study, we assessed 167 probands referred for FBN2 analysis and classified them into a FBN2-positive (n = 44) and FBN2-negative group (n = 123) following molecular analysis. We developed a 20-point weighted clinical scoring system based on the prevalence of ten main clinical characteristics of CCA in both groups. RESULTS: The total score was significantly different between the groups (P < 0.001) and was indicative for classifying patients into unlikely CCA (total score <7) and likely CCA (total score ≥7) groups. CONCLUSIONS: Our clinical score is helpful for clinical guidance for patients suspected to have CCA, and provides a quantitative tool for phenotyping in research settings.

Myhre syndrome: A first familial recurrence and broadening of the phenotypic spectrum.

Myhre syndrome is a rare multisystem connective tissue disorder, characterized by short stature, facial dysmorphology, variable intellectual disability, skeletal abnormalities, arthropathy, cardiopathy, laryngotracheal anomalies, and stiff skin. So far, all molecularly confirmed cases harbored a de novo heterozygous gain-of-function mutation in SMAD4, encoding the SMAD4 transducer protein required for both transforming growth factor-beta and bone morphogenic proteins signaling. We report on four novel patients (one female proband and her two affected children, and one male proband) with Myhre syndrome harboring the recurrent c.1486C>T (p.Arg496Cys) mutation in SMAD4. The female proband presented with a congenital heart defect, vertebral anomalies, and facial dysmorphic features. She developed severe tracheal stenosis requiring a total laryngectomy. With assisted reproductive treatment, she gave birth to two affected children. The second proband presented with visual impairment following lensectomy in childhood, short stature, brachydactyly, stiff skin, and decreased peripheral sensitivity. Transmission electron microscopy (TEM) of the dermis shows irregular elastin cores with globular deposits and almost absent surrounding microfibrils and suggests age-related increased collagen deposition. We report on the first familial case of Myhre syndrome and illustrate the variable clinical spectrum of the disorder. Despite the primarily fibrotic nature of the disease, TEM analysis mainly indicates elastic fiber anomalies.

CRISPR/Cas9-mediated homology-directed repair by ssODNs in zebrafish induces complex mutational patterns resulting from genomic integration of repair-template fragments.

Targeted genome editing by CRISPR/Cas9 is extremely well fitted to generate gene disruptions, although precise sequence replacement by CRISPR/Cas9-mediated homology-directed repair (HDR) suffers from low efficiency, impeding its use for high-throughput knock-in disease modeling. In this study, we used next-generation sequencing (NGS) analysis to determine the efficiency and reliability of CRISPR/Cas9-mediated HDR using several types of single-stranded oligodeoxynucleotide (ssODN) repair templates for the introduction of disease-relevant point mutations in the zebrafish genome. Our results suggest that HDR rates are strongly determined by repair-template composition, with the most influential factor being homology-arm length. However, we found that repair using ssODNs does not only lead to precise sequence replacement but also induces integration of repair-template fragments at the Cas9 cut site. We observed that error-free repair occurs at a relatively constant rate of 1-4% when using different repair templates, which was sufficient for transmission of point mutations to the F1 generation. On the other hand, erroneous repair mainly accounts for the variability in repair rate between the different repair templates. To further improve error-free HDR rates, elucidating the mechanism behind this erroneous repair is essential. We show that the error-prone nature of ssODN-mediated repair, believed to act via synthesis-dependent strand annealing (SDSA), is most likely due to DNA synthesis errors. In conclusion, caution is warranted when using ssODNs for the generation of knock-in models or for therapeutic applications. We recommend the application of in-depth NGS analysis to examine both the efficiency and error-free nature of HDR events.This article has an associated First Person interview with the first author of the paper.

Publisher Correction: BATCH-GE: Batch analysis of Next-Generation Sequencing data for genome editing assessment.

A correction to this article has been published and is linked from the HTML and PDF versions of this paper. The error has not been fixed in the paper.

BATCH-GE: Analysis of NGS Data for Genome Editing Assessment.

Due to its simple nature, the clustered regularly interspaced short palindromic repeats (CRISPR)-Cas9 technique is massively used nowadays to modify genomic loci in a wide range of model systems. The possibility to interrogate gene function on a genome-wide scale is revolutionizing fundamental life sciences and will lead to new clinical breakthroughs. Its strength is even more pronounced when it is used in tandem with next-generation sequencing (NGS). The high throughput and low cost cause NGS to be the method of choice for exploring CRISPR-Cas9 experimental results. To analyze the NGS reads from genome editing experiments only few bioinformatics tools are available. BATCH-GE is a flexible and easy-to-use tool, which is especially useful for dealing with large amounts of data. It detects and reports indel mutations and other precise genome editing events and calculates the corresponding mutagenesis efficiencies for a large number of samples in parallel.