Long-read trio sequencing of individuals with unsolved intellectual disability

The Third-Generation Sequencing Challenge: Novel Insights for the Omic Sciences

The understanding of the human genome has been greatly improved by the advent of next-generation sequencing technologies (NGS). Despite the undeniable advantages responsible for their widespread diffusion, these methods have some constraints, mainly related to short read length and the need for PCR amplification. As a consequence, long-read sequencers, called third-generation sequencing (TGS), have been developed, promising to overcome NGS. Starting from the first prototype, TGS has progressively ameliorated its chemistries by improving both read length and base-calling accuracy, as well as simultaneously reducing the costs/base. Based on these premises, TGS is showing its potential in many fields, including the analysis of difficult-to-sequence genomic regions, structural variations detection, RNA expression profiling, DNA methylation study, and metagenomic analyses. Protocol standardization and the development of easy-to-use pipelines for data analysis will enhance TGS use, also opening the way for their routine applications in diagnostic contexts.

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.

Narrowing the diagnostic gap: Genomes, episignatures, long-read sequencing, and health economic analyses in an exome-negative intellectual disability cohort

PURPOSE

Genome sequencing (GS)-specific diagnostic rates in prospective tightly ascertained exome sequencing (ES)-negative intellectual disability (ID) cohorts have not been reported extensively.

METHODS

ES, GS, epigenetic signatures, and long-read sequencing diagnoses were assessed in 74 trios with at least moderate ID.

RESULTS

The ES diagnostic yield was 42 of 74 (57%). GS diagnoses were made in 9 of 32 (28%) ES-unresolved families. Repeated ES with a contemporary pipeline on the GS-diagnosed families identified 8 of 9 single-nucleotide variations/copy-number variations undetected in older ES, confirming a GS-unique diagnostic rate of 1 in 32 (3%). Episignatures contributed diagnostic information in 9% with GS corroboration in 1 of 32 (3%) and diagnostic clues in 2 of 32 (6%). A genetic etiology for ID was detected in 51 of 74 (69%) families. Twelve candidate disease genes were identified. Contemporary ES followed by GS cost US$4976 (95% CI: $3704; $6969) per diagnosis and first-line GS at a cost of $7062 (95% CI: $6210; $8475) per diagnosis.

CONCLUSION

Performing GS only in ID trios would be cost equivalent to ES if GS were available at $2435, about a 60% reduction from current prices. This study demonstrates that first-line GS achieves higher diagnostic rate than contemporary ES but at a higher cost.

Genome sequencing as a generic diagnostic strategy for rare disease

BACKGROUND

To diagnose the full spectrum of hereditary and congenital diseases, genetic laboratories use many different workflows, ranging from karyotyping to exome sequencing. A single generic high-throughput workflow would greatly increase efficiency. We assessed whether genome sequencing (GS) can replace these existing workflows aimed at germline genetic diagnosis for rare disease.

METHODS

We performed short-read GS (NovaSeq™6000; 150 bp paired-end reads, 37 × mean coverage) on 1000 cases with 1271 known clinically relevant variants, identified across different workflows, representative of our tertiary diagnostic centers. Variants were categorized into small variants (single nucleotide variants and indels < 50 bp), large variants (copy number variants and short tandem repeats) and other variants (structural variants and aneuploidies). Variant calling format files were queried per variant, from which workflow-specific true positive rates (TPRs) for detection were determined. A TPR of ≥ 98% was considered the threshold for transition to GS. A GS-first scenario was generated for our laboratory, using diagnostic efficacy and predicted false negative as primary outcome measures. As input, we modeled the diagnostic path for all 24,570 individuals referred in 2022, combining the clinical referral, the transition of the underlying workflow(s) to GS, and the variant type(s) to be detected.

RESULTS

Overall, 95% (1206/1271) of variants were detected. Detection rates differed per variant category: small variants in 96% (826/860), large variants in 93% (341/366), and other variants in 87% (39/45). TPRs varied between workflows (79-100%), with 7/10 being replaceable by GS. Models for our laboratory indicate that a GS-first strategy would be feasible for 84.9% of clinical referrals (750/883), translating to 71% of all individuals (17,444/24,570) receiving GS as their primary test. An estimated false negative rate of 0.3% could be expected.

CONCLUSIONS

GS can capture clinically relevant germline variants in a 'GS-first strategy' for the majority of clinical indications in a genetics diagnostic lab.

Bioinformatics of germline variant discovery for rare disease diagnostics: current approaches and remaining challenges

Next-generation sequencing (NGS) has revolutionized the field of rare disease diagnostics. Whole exome and whole genome sequencing are now routinely used for diagnostic purposes; however, the overall diagnosis rate remains lower than expected. In this work, we review current approaches used for calling and interpretation of germline genetic variants in the human genome, and discuss the most important challenges that persist in the bioinformatic analysis of NGS data in medical genetics. We describe and attempt to quantitatively assess the remaining problems, such as the quality of the reference genome sequence, reproducible coverage biases, or variant calling accuracy in complex regions of the genome. We also discuss the prospects of switching to the complete human genome assembly or the human pan-genome and important caveats associated with such a switch. We touch on arguably the hardest problem of NGS data analysis for medical genomics, namely, the annotation of genetic variants and their subsequent interpretation. We highlight the most challenging aspects of annotation and prioritization of both coding and non-coding variants. Finally, we demonstrate the persistent prevalence of pathogenic variants in the coding genome, and outline research directions that may enhance the efficiency of NGS-based disease diagnostics.

Reclassifying variations of unknown significance in diseases affecting Saudi Arabia's population reveal new associations

Introduction: Physicians face diagnostic dilemmas upon reports indicating disease variants of unknown significance (VUS). The most puzzling cases are patients with rare diseases, where finding another matched genotype and phenotype to associate their results is challenging. This study aims to prove the value of updating patient files with new classifications, potentially leading to better assessment and prevention. Methodology: We recruited retrospective phenotypic and genotypic data from King Saud Medical City, Riyadh, Kingdom of Saudi Arabia. Between September 2020 and December 2021, 1,080 patients' genetic profiles were tested in a College of American Pathologists accredited laboratory. We excluded all confirmed pathogenic variants, likely pathogenic variants and copy number variations. Finally, we further reclassified 194 VUS using different local and global databases, employing in silico prediction to justify the phenotype-genotype association. Results: Of the 194 VUS, 90 remained VUS, and the other 104 were reclassified as follows: 16 pathogenic, 49 likely pathogenic, nine benign, and 30 likely benign. Moreover, most of these variants had never been observed in other local or international databases. Conclusion: Reclassifying the VUS adds value to understanding the causality of the phenotype if it has been reported in another family or population. The healthcare system should establish guidelines for re-evaluating VUS, and upgrading VUS should reflect on individual/family risks and management strategies.

Comparison and benchmark of structural variants detected from long read and long-read assembly

Structural variant (SV) detection is essential for genomic studies, and long-read sequencing technologies have advanced our capacity to detect SVs directly from read or de novo assembly, also known as read-based and assembly-based strategy. However, to date, no independent studies have compared and benchmarked the two strategies. Here, on the basis of SVs detected by 20 read-based and eight assembly-based detection pipelines from six datasets of HG002 genome, we investigated the factors that influence the two strategies and assessed their performance with well-curated SVs. We found that up to 80% of the SVs could be detected by both strategies among different long-read datasets, whereas variant type, size, and breakpoint detected by read-based strategy were greatly affected by aligners. For the high-confident insertions and deletions at non-tandem repeat regions, a remarkable subset of them (82% in assembly-based calls and 93% in read-based calls), accounting for around 4000 SVs, could be captured by both reads and assemblies. However, discordance between two strategies was largely caused by complex SVs and inversions, which resulted from inconsistent alignment of reads and assemblies at these loci. Finally, benchmarking with SVs at medically relevant genes, the recall of read-based strategy reached 77% on 5X coverage data, whereas assembly-based strategy required 20X coverage data to achieve similar performance. Therefore, integrating SVs from read and assembly is suggested for general-purpose detection because of inconsistently detected complex SVs and inversions, whereas assembly-based strategy is optional for applications with limited resources.

Comprehensive de novo mutation discovery with HiFi long-read sequencing

BACKGROUND

Long-read sequencing (LRS) techniques have been very successful in identifying structural variants (SVs). However, the high error rate of LRS made the detection of small variants (substitutions and short indels < 20 bp) more challenging. The introduction of PacBio HiFi sequencing makes LRS also suited for detecting small variation. Here we evaluate the ability of HiFi reads to detect de novo mutations (DNMs) of all types, which are technically challenging variant types and a major cause of sporadic, severe, early-onset disease.

METHODS

We sequenced the genomes of eight parent-child trios using high coverage PacBio HiFi LRS (~ 30-fold coverage) and Illumina short-read sequencing (SRS) (~ 50-fold coverage). De novo substitutions, small indels, short tandem repeats (STRs) and SVs were called in both datasets and compared to each other to assess the accuracy of HiFi LRS. In addition, we determined the parent-of-origin of the small DNMs using phasing.

RESULTS

We identified a total of 672 and 859 de novo substitutions/indels, 28 and 126 de novo STRs, and 24 and 1 de novo SVs in LRS and SRS respectively. For the small variants, there was a 92 and 85% concordance between the platforms. For the STRs and SVs, the concordance was 3.6 and 0.8%, and 4 and 100% respectively. We successfully validated 27/54 LRS-unique small variants, of which 11 (41%) were confirmed as true de novo events. For the SRS-unique small variants, we validated 42/133 DNMs and 8 (19%) were confirmed as true de novo event. Validation of 18 LRS-unique de novo STR calls confirmed none of the repeat expansions as true DNM. Confirmation of the 23 LRS-unique SVs was possible for 19 candidate SVs of which 10 (52.6%) were true de novo events. Furthermore, we were able to assign 96% of DNMs to their parental allele with LRS data, as opposed to just 20% with SRS data.

CONCLUSIONS

HiFi LRS can now produce the most comprehensive variant dataset obtainable by a single technology in a single laboratory, allowing accurate calling of substitutions, indels, STRs and SVs. The accuracy even allows sensitive calling of DNMs on all variant levels, and also allows for phasing, which helps to distinguish true positive from false positive DNMs.

Twist exome capture allows for lower average sequence coverage in clinical exome sequencing

BACKGROUND

Exome and genome sequencing are the predominant techniques in the diagnosis and research of genetic disorders. Sufficient, uniform and reproducible/consistent sequence coverage is a main determinant for the sensitivity to detect single-nucleotide (SNVs) and copy number variants (CNVs). Here we compared the ability to obtain comprehensive exome coverage for recent exome capture kits and genome sequencing techniques.

RESULTS

We compared three different widely used enrichment kits (Agilent SureSelect Human All Exon V5, Agilent SureSelect Human All Exon V7 and Twist Bioscience) as well as short-read and long-read WGS. We show that the Twist exome capture significantly improves complete coverage and coverage uniformity across coding regions compared to other exome capture kits. Twist performance is comparable to that of both short- and long-read whole genome sequencing. Additionally, we show that even at a reduced average coverage of 70× there is only minimal loss in sensitivity for SNV and CNV detection.

CONCLUSION

We conclude that exome sequencing with Twist represents a significant improvement and could be performed at lower sequence coverage compared to other exome capture techniques.

High diagnostic potential of short and long read genome sequencing with transcriptome analysis in exome-negative developmental disorders

Exome sequencing (ES) has become the method of choice for diagnosing rare diseases, while the availability of short-read genome sequencing (SR-GS) in a medical setting is increasing. In addition, new sequencing technologies, such as long-read genome sequencing (LR-GS) and transcriptome sequencing, are being increasingly used. However, the contribution of these techniques compared to widely used ES is not well established, particularly in regards to the analysis of non-coding regions. In a pilot study of five probands affected by an undiagnosed neurodevelopmental disorder, we performed trio-based short-read GS and long-read GS as well as case-only peripheral blood transcriptome sequencing. We identified three new genetic diagnoses, none of which affected the coding regions. More specifically, LR-GS identified a balanced inversion in NSD1, highlighting a rare mechanism of Sotos syndrome. SR-GS identified a homozygous deep intronic variant of KLHL7 resulting in a neoexon inclusion, and a de novo mosaic intronic 22-bp deletion in KMT2D, leading to the diagnosis of Perching and Kabuki syndromes, respectively. All three variants had a significant effect on the transcriptome, which showed decreased gene expression, mono-allelic expression and splicing defects, respectively, further validating the effect of these variants. Overall, in undiagnosed patients, the combination of short and long read GS allowed the detection of cryptic variations not or barely detectable by ES, making it a highly sensitive method at the cost of more complex bioinformatics approaches. Transcriptome sequencing is a valuable complement for the functional validation of variations, particularly in the non-coding genome.

Biallelic CC2D2A variants, SNV and LINE-1 insertion simultaneously identified in siblings using long-read whole-genome sequencing and haplotype phasing

Joubert syndrome (JBTS) is characterized by a magnetic resonance imaging appearance called 'molar tooth sign', neonatal breathing dysregulation and hypotonia, and developmental delay. Whole-exome analysis based on short-read sequencing has often contributed to the identification of causative single-nucleotide variants in patients clinically diagnosed with JBTS. However, ~10% of them are still undiagnosed even though a single possible pathogenic variant has been identified. We report a successful identification of biallelic variants using long-read whole-genome sequencing and haplotype phasing analysis in a family with two Japanese siblings having morphological brain abnormalities. The affected siblings had a novel nonsynonymous variant (CC2D2A:NM_001080522.2:c.4454A>G:p.(Tyr1485Cys)) and an exonic insertion of Long INterspercsed Element-1 (LINE-1). The allelicity of these variants was clearly proven without the data of parents. Finally, our survey of in-house genome sequencing data indicates that there are rare carriers of CC2D2A related diseases, who harbour the exonic LINE-1 insertion in the CC2D2A gene.

The Genetics of Intellectual Disability

Intellectual disability (ID) has a prevalence of ~2-3% in the general population, having a large societal impact. The underlying cause of ID is largely of genetic origin; however, identifying this genetic cause has in the past often led to long diagnostic Odysseys. Over the past decades, improvements in genetic diagnostic technologies and strategies have led to these causes being more and more detectable: from cytogenetic analysis in 1959, we moved in the first decade of the 21st century from genomic microarrays with a diagnostic yield of ~20% to next-generation sequencing platforms with a yield of up to 60%. In this review, we discuss these various developments, as well as their associated challenges and implications for the field of ID, which highlight the revolutionizing shift in clinical practice from a phenotype-first into genotype-first approach.

Approaches to long-read sequencing in a clinical setting to improve diagnostic rate

Over the past decade, advances in genetic testing, particularly the advent of next-generation sequencing, have led to a paradigm shift in the diagnosis of molecular diseases and disorders. Despite our present collective ability to interrogate more than 90% of the human genome, portions of the genome have eluded us, resulting in stagnation of diagnostic yield with existing methodologies. Here we show how application of a new technology, long-read sequencing, has the potential to improve molecular diagnostic rates. Whole genome sequencing by long reads was able to cover 98% of next-generation sequencing dead zones, which are areas of the genome that are not interpretable by conventional industry-standard short-read sequencing. Through the ability of long-read sequencing to unambiguously call variants in these regions, we discovered an immunodeficiency due to a variant in IKBKG in a subject who had previously received a negative genome sequencing result. Additionally, we demonstrate the ability of long-read sequencing to detect small variants on par with short-read sequencing, its superior performance in identifying structural variants, and thirdly, its capacity to determine genomic methylation defects in native DNA. Though the latter technical abilities have been demonstrated, we demonstrate the clinical application of this technology to successfully identify multiple types of variants using a single test.

Developing CIRdb as a catalog of natural genetic variation in the Canary Islanders

The current inhabitants of the Canary Islands have a unique genetic makeup in the European diversity landscape due to the existence of African footprints from recent admixture events, especially of North African components (> 20%). The underrepresentation of non-Europeans in genetic studies and the sizable North African ancestry, which is nearly absent from all existing catalogs of worldwide genetic diversity, justify the need to develop CIRdb, a population-specific reference catalog of natural genetic variation in the Canary Islanders. Based on array genotyping of the selected unrelated donors and comparisons against available datasets from European, sub-Saharan, and North African populations, we illustrate the intermediate genetic differentiation of Canary Islanders between Europeans and North Africans and the existence of within-population differences that are likely driven by genetic isolation. Here we describe the overall design and the methods that are being implemented to further develop CIRdb. This resource will help to strengthen the implementation of Precision Medicine in this population by contributing to increase the diversity in genetic studies. Among others, this will translate into improved ability to fine map disease genes and simplify the identification of causal variants and estimate the prevalence of unattended Mendelian diseases.

Construction of a trio-based structural variation panel utilizing activated T lymphocytes and long-read sequencing technology

Long-read sequencing technology enable better characterization of structural variants (SVs). To adapt the technology to population-scale analyses, one critical issue is to obtain sufficient amount of high-molecular-weight genomic DNA. Here, we propose utilizing activated T lymphocytes, which can be established efficiently in a biobank to stably supply high-grade genomic DNA sufficiently. We conducted nanopore sequencing of 333 individuals constituting 111 trios with high-coverage long-read sequencing data (depth 22.2x, N50 of 25.8 kb) and identified 74,201 SVs. Our trio-based analysis revealed that more than 95% of the SVs were concordant with Mendelian inheritance. We also identified SVs associated with clinical phenotypes, all of which appear to be stably transmitted from parents to offspring. Our data provide a catalog of SVs in the general Japanese population, and the applied approach using the activated T-lymphocyte resource will contribute to biobank-based human genetic studies focusing on SVs at the population scale.

Long-read sequencing on activated T-cells from a sample of 333 Japanese individuals (representing 111 parent-offspring trios) provides a useful reference dataset of structural variation in the Japanese population.

How to proceed after "negative" exome: A review on genetic diagnostics, limitations, challenges, and emerging new multiomics techniques

Exome sequencing (ES) in the clinical setting of inborn metabolic diseases (IMDs) has created tremendous improvement in achieving an accurate and timely molecular diagnosis for a greater number of patients, but it still leaves the majority of patients without a diagnosis. In parallel, (personalized) treatment strategies are increasingly available, but this requires the availability of a molecular diagnosis. IMDs comprise an expanding field with the ongoing identification of novel disease genes and the recognition of multiple inheritance patterns, mosaicism, variable penetrance, and expressivity for known disease genes. The analysis of trio ES is preferred over singleton ES as information on the allelic origin (paternal, maternal, "de novo") reduces the number of variants that require interpretation. All ES data and interpretation strategies should be exploited including CNV and mitochondrial DNA analysis. The constant advancements in available techniques and knowledge necessitate the close exchange of clinicians and molecular geneticists about genotypes and phenotypes, as well as knowledge of the challenges and pitfalls of ES to initiate proper further diagnostic steps. Functional analyses (transcriptomics, proteomics, and metabolomics) can be applied to characterize and validate the impact of identified variants, or to guide the genomic search for a diagnosis in unsolved cases. Future diagnostic techniques (genome sequencing [GS], optical genome mapping, long-read sequencing, and epigenetic profiling) will further enhance the diagnostic yield. We provide an overview of the challenges and limitations inherent to ES followed by an outline of solutions and a clinical checklist, focused on establishing a diagnosis to eventually achieve (personalized) treatment.

Reanalysis of exome negative patients with rare disease: a pragmatic workflow for diagnostic applications

Background

Approximately two third of patients with a rare genetic disease remain undiagnosed after exome sequencing (ES). As part of our post-test counseling procedures, patients without a conclusive diagnosis are advised to recontact their referring clinician to discuss new diagnostic opportunities in due time. We performed a systematic study of genetically undiagnosed patients 5 years after their initial negative ES report to determine the efficiency of diverse reanalysis strategies.

Methods

We revisited a cohort of 150 pediatric neurology patients originally enrolled at Radboud University Medical Center, of whom 103 initially remained genetically undiagnosed. We monitored uptake of physician-initiated routine clinical and/or genetic re-evaluation (ad hoc re-evaluation) and performed systematic reanalysis, including ES-based resequencing, of all genetically undiagnosed patients (systematic re-evaluation).

Results

Ad hoc re-evaluation was initiated for 45 of 103 patients and yielded 18 diagnoses (including 1 non-genetic). Subsequent systematic re-evaluation identified another 14 diagnoses, increasing the diagnostic yield in our cohort from 31% (47/150) to 53% (79/150). New genetic diagnoses were established by reclassification of previously identified variants (10%, 3/31), reanalysis with enhanced bioinformatic pipelines (19%, 6/31), improved coverage after resequencing (29%, 9/31), and new disease-gene associations (42%, 13/31). Crucially, our systematic study also showed that 11 of the 14 further conclusive genetic diagnoses were made in patients without a genetic diagnosis that did not recontact their referring clinician.

Conclusions

We find that upon re-evaluation of undiagnosed patients, both reanalysis of existing ES data as well as resequencing strategies are needed to identify additional genetic diagnoses. Importantly, not all patients are routinely re-evaluated in clinical care, prolonging their diagnostic trajectory, unless systematic reanalysis is facilitated. We have translated our observations into considerations for systematic and ad hoc reanalysis in routine genetic care.

Supplementary Information

The online version contains supplementary material available at 10.1186/s13073-022-01069-z.

SvAnna: efficient and accurate pathogenicity prediction of coding and regulatory structural variants in long-read genome sequencing

Structural variants (SVs) are implicated in the etiology of Mendelian diseases but have been systematically underascertained owing to sequencing technology limitations. Long-read sequencing enables comprehensive detection of SVs, but approaches for prioritization of candidate SVs are needed. Structural variant Annotation and analysis (SvAnna) assesses all classes of SVs and their intersection with transcripts and regulatory sequences, relating predicted effects on gene function with clinical phenotype data. SvAnna places 87% of deleterious SVs in the top ten ranks. The interpretable prioritizations offered by SvAnna will facilitate the widespread adoption of long-read sequencing in diagnostic genomics. SvAnna is available at https://github.com/TheJacksonLaboratory/SvAnn a .

Familial long-read sequencing increases yield of de novo mutations

Studies of de novo mutation (DNM) have typically excluded some of the most repetitive and complex regions of the genome because these regions cannot be unambiguously mapped with short-read sequencing data. To better understand the genome-wide pattern of DNM, we generated long-read sequence data from an autism parent-child quad with an affected female where no pathogenic variant had been discovered in short-read Illumina sequence data. We deeply sequenced all four individuals by using three sequencing platforms (Illumina, Oxford Nanopore, and Pacific Biosciences) and three complementary technologies (Strand-seq, optical mapping, and 10X Genomics). Using long-read sequencing, we initially discovered and validated 171 DNMs across two children-a 20% increase in the number of de novo single-nucleotide variants (SNVs) and indels when compared to short-read callsets. The number of DNMs further increased by 5% when considering a more complete human reference (T2T-CHM13) because of the recovery of events in regions absent from GRCh38 (e.g., three DNMs in heterochromatic satellites). In total, we validated 195 de novo germline mutations and 23 potential post-zygotic mosaic mutations across both children; the overall true substitution rate based on this integrated callset is at least 1.41 × 10^-8 substitutions per nucleotide per generation. We also identified six de novo insertions and deletions in tandem repeats, two of which represent structural variants. We demonstrate that long-read sequencing and assembly, especially when combined with a more complete reference genome, increases the number of DNMs by >25% compared to previous studies, providing a more complete catalog of DNM compared to short-read data alone.

MRSD: A quantitative approach for assessing suitability of RNA-seq in the investigation of mis-splicing in Mendelian disease

Variable levels of gene expression between tissues complicates the use of RNA sequencing of patient biosamples to delineate the impact of genomic variants. Here, we describe a gene- and tissue-specific metric to inform the feasibility of RNA sequencing. This overcomes limitations of using expression values alone as a metric to predict RNA-sequencing utility. We have derived a metric, minimum required sequencing depth (MRSD), that estimates the depth of sequencing required from RNA sequencing to achieve user-specified sequencing coverage of a gene, transcript, or group of genes. We applied MRSD across four human biosamples: whole blood, lymphoblastoid cell lines (LCLs), skeletal muscle, and cultured fibroblasts. MRSD has high precision (90.1%-98.2%) and overcomes transcript region-specific sequencing biases. Applying MRSD scoring to established disease gene panels shows that fibroblasts, of these four biosamples, are the optimum source of RNA for 63.1% of gene panels. Using this approach, up to 67.8% of the variants of uncertain significance in ClinVar that are predicted to impact splicing could be assayed by RNA sequencing in at least one of the biosamples. We demonstrate the utility and benefits of MRSD as a metric to inform functional assessment of splicing aberrations, in particular in the context of Mendelian genetic disorders to improve diagnostic yield.

Xq27.1 palindrome mediated interchromosomal insertion likely causes familial congenital bilateral laryngeal abductor paralysis (Plott syndrome)

Bilateral laryngeal abductor paralysis is a rare entity and the second most common cause of stridor in newborns. So far, no conclusive genetic or chromosomal aberration has been reported for X-linked isolated bilateral vocal cord paralysis, also referred to as Plott syndrome. Via whole genome sequencing (WGS), we identified a complex interchromosomal insertion in a large family with seven affected males. The 404 kb inserted fragment originates from chromosome 10q21.3, contains no genes and is inserted inversionally into the intergenic chromosomal region Xq27.1, 82 kb centromeric to the nearest gene SOX3. The patterns found at the breakpoint junctions resemble typical characteristics that arise in replication-based mechanisms with long-distance template switching. Non protein-coding insertions into the same genomic region have been described to result in different phenotypes, indicating that the phenotypic outcome likely depends on the introduction of regulatory elements. In conclusion, our data adds Plott syndrome as another entity, likely caused by the insertion of non-coding DNA into the intergenic chromosomal region Xq27.1. In this regard, we demonstrate the importance of WGS as a powerful diagnostic test in unsolved genetic diseases, as this genomic rearrangement has not been detected by current first-line diagnostic tests, i.e., exome sequencing and chromosomal microarray analysis.

Long-read technologies identify a hidden inverted duplication in a family with choroideremia

The lack of molecular diagnoses in rare genetic diseases can be explained by limitations of current standard genomic technologies. Upcoming long-read techniques have complementary strengths to overcome these limitations, with a particular strength in identifying structural variants. By using optical genome mapping and long-read sequencing, we aimed to identify the pathogenic variant in a large family with X-linked choroideremia. In this family, aberrant splicing of exon 12 of the choroideremia gene CHM was detected in 2003, but the underlying genomic defect remained elusive. Optical genome mapping and long-read sequencing approaches now revealed an intragenic 1,752 bp inverted duplication including exon 12 and surrounding regions, located downstream of the wild-type copy of exon 12. Both breakpoint junctions were confirmed with Sanger sequencing and segregate with the X-linked inheritance in the family. The breakpoint junctions displayed sequence microhomology suggestive for an erroneous replication mechanism as the origin of the structural variant. The inverted duplication is predicted to result in a hairpin formation of the pre-mRNA with the wild-type exon 12, leading to exon skipping in the mature mRNA. The identified inverted duplication is deemed the hidden pathogenic cause of disease in this family. Our study shows that optical genome mapping and long-read sequencing have significant potential for the identification of (hidden) structural variants in rare genetic diseases.

Diagnostic Analyses of Retinal Dystrophy Genes: Current Status and Perspective

Over the past decade, novel high-throughput DNA sequencing technologies have revolutionised both research and diagnostic testing for monogenic disorders. This applies particularly to genetically very heterogeneous disorders like retinal dystrophies (RDs). Next-generation sequencing (NGS) today is considered as reliable as Sanger sequencing, which had been the gold standard for decades. Today, comprehensive NGS-based diagnostic testing reveals the causative mutations in the majority of RD patients, with important implications for genetic counselling for recurrence risks and personalised medical management (from interdisciplinary surveillance to prophylactic measures and, albeit yet rare, [gene] therapy). While DNA sequencing is - in most cases - no longer the diagnostic bottleneck, one needs to be aware of interpretation pitfalls and dead ends. The advent of new (NGS) technologies will solve some of these issues. However, specialised medical geneticists who are familiar with the peculiarities of certain RD genes and closely interact with ophthalmologists will remain key to successful RD research and diagnostic testing for the benefit of the patients. This review sheds light on the current state of the field, its challenges and potential solutions.