A 12-kb structural variation in progressive myoclonic epilepsy was newly identified by long-read whole-genome sequencing

Long-read sequencing and optical genome mapping identify causative gene disruptions in noncoding sequence in two patients with neurologic disease and known chromosome abnormalities

Despite advances in next generation sequencing (NGS), genetic diagnoses remain elusive for many patients with neurologic syndromes. Long-read sequencing (LRS) and optical genome mapping (OGM) technologies improve upon existing capabilities in the detection and interpretation of structural variation in repetitive DNA, on a single haplotype, while also providing enhanced breakpoint resolution. We performed LRS and OGM on two patients with known chromosomal rearrangements and inconclusive Sanger or NGS. The first patient, who had epilepsy and developmental delay, had a complex translocation between two chromosomes that included insertion and inversion events. The second patient, who had a movement disorder, had an inversion on a single chromosome disrupted by multiple smaller inversions and insertions. Sequence level resolution of the rearrangements identified pathogenic breaks in noncoding sequence in or near known disease-causing genes with relevant neurologic phenotypes (MBD5, NKX2-1). These specific variants have not been reported previously, but expected molecular consequences are consistent with previously reported cases. As the use of LRS and OGM technologies for clinical testing increases and data analyses become more standardized, these methods along with multiomic data to validate noncoding variation effects will improve diagnostic yield and increase the proportion of probands with detectable pathogenic variants for known genes implicated in neurogenetic disease.

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.

Long read sequencing on its way to the routine diagnostics of genetic diseases

The clinical application of technological progress in the identification of DNA alterations has always led to improvements of diagnostic yields in genetic medicine. At chromosome side, from cytogenetic techniques evaluating number and gross structural defects to genomic microarrays detecting cryptic copy number variants, and at molecular level, from Sanger method studying the nucleotide sequence of single genes to the high-throughput next-generation sequencing (NGS) technologies, resolution and sensitivity progressively increased expanding considerably the range of detectable DNA anomalies and alongside of Mendelian disorders with known genetic causes. However, particular genomic regions (i.e., repetitive and GC-rich sequences) are inefficiently analyzed by standard genetic tests, still relying on laborious, time-consuming and low-sensitive approaches (i.e., southern-blot for repeat expansion or long-PCR for genes with highly homologous pseudogenes), accounting for at least part of the patients with undiagnosed genetic disorders. Third generation sequencing, generating long reads with improved mappability, is more suitable for the detection of structural alterations and defects in hardly accessible genomic regions. Although recently implemented and not yet clinically available, long read sequencing (LRS) technologies have already shown their potential in genetic medicine research that might greatly impact on diagnostic yield and reporting times, through their translation to clinical settings. The main investigated LRS application concerns the identification of structural variants and repeat expansions, probably because techniques for their detection have not evolved as rapidly as those dedicated to single nucleotide variants (SNV) identification: gold standard analyses are karyotyping and microarrays for balanced and unbalanced chromosome rearrangements, respectively, and southern blot and repeat-primed PCR for the amplification and sizing of expanded alleles, impaired by limited resolution and sensitivity that have not been significantly improved by the advent of NGS. Nevertheless, more recently, with the increased accuracy provided by the latest product releases, LRS has been tested also for SNV detection, especially in genes with highly homologous pseudogenes and for haplotype reconstruction to assess the parental origin of alleles with de novo pathogenic variants. We provide a review of relevant recent scientific papers exploring LRS potential in the diagnosis of genetic diseases and its potential future applications in routine genetic testing.

Whole exome sequencing identifies variable expressivity of CLN6 variants in Progressive myoclonic epilepsy affected families

Progressive myoclonic epilepsies (PMEs) are a group of neurodegenerative disorders, predominantly affecting adolescents and, characterized by generalized worsening myoclonus epilepsies, ataxia, cognitive deficits, and dementia. To date, several genes, having implications in diverse phenotypic expressions associated with PMEs, have been identified. Genetic diagnosis is available for most of the adolescence-onset myoclonic epilepsies. This study aimed to elucidate the genetic basis of PMEs in three multiplex Pakistani families exhibiting clinically variable phenotypes. Causative variant(s) in the studied families, and mode of segregation were identified by Whole Exome Sequencing (WES) of the probands, followed by bi-directional Sanger sequencing for final validation. We identified homozygous recessive CLN6 missense variant c.768 C>G (p.Asp256Glu) in Family 1, and c.889 C>A (p.Pro297Thr) variant in Family 2. While in Family 3, we found a homozygous variant (c.316dup) that caused a frameshift mutation, leading to a premature stop codon in the CLN6 protein, resulting in a truncated protein (p.Arg106ProfsTer26). Though CLN6 is previously identified to underlie late infantile and adolescent onset neuronal ceroid lipofuscinosis, this study supports and expands the phenotypic spectrum of CLN6 mutations and signifies diagnositc potential CLN6 variants for PMEs. Diverse pathological effects of variant c .768 C>G were observed in Family 1, with same genotypes, suggesting clinical heterogeneity and/or variable expressivity that might be the implication of pleiotropic effects of the gene in these cases.

Targeted long-read sequencing identified a causal structural variant in X-linked nephrogenic diabetes insipidus

BACKGROUND

X-linked nephrogenic diabetes insipidus (NDI) is a rare genetic renal disease caused by pathogenic variants in the AVPR2 gene. Single nucleotide variants and small insertions/deletions in AVPR2 are reliably detected by routine clinical sequencing. Nevertheless, structural variants involving AVPR2 are challenging to identify accurately by conventional genetic testing. Here, we report a novel deletion of AVPR2 in a Czech family identified for the first time by targeted long-read sequencing (T-LRS).

METHODS

A male proband with X-linked NDI underwent clinical sequencing of the AVPR2 gene that failed and thus indicated possible whole-gene deletion. Therefore, PCR mapping and subsequent targeted long-read sequencing (T-LRS) using a Pacific Biosciences sequencer were applied to search for the suspected deletion. To validate the deletion breakpoints and prove variant segregation in the family with X-linked NDI, Sanger sequencing of the deletion junction was performed. Quantitative real-time PCR was further carried out to confirm the carrier status of heterozygous females.

RESULTS

By T-LRS, a novel 7.5 kb deletion of AVPR2 causing X-linked NDI in the proband was precisely identified. Sanger sequencing of the deletion junction confirmed the variant breakpoints and detected the deletion in the probands´ mother, maternal aunt, and maternal cousin with X-linked NDI. The carrier status in heterozygous females was further validated by quantitative real-time PCR.

CONCLUSIONS

Identifying the 7.5 kb deletion gave a precise molecular diagnosis for the proband, enabled genetic counselling and genetic testing for the family, and further expanded the spectrum of structural variants causing X-linked NDI. Our results also show that T-LRS has significant potential for accurately identifying putative structural variants.

Application of long read sequencing in rare diseases: The longer, the better?

Rare diseases encompass a diverse group of genetic disorders that affect a small proportion of the population. Identifying the underlying genetic causes of these conditions presents significant challenges due to their genetic heterogeneity and complexity. Conventional short-read sequencing (SRS) techniques have been widely used in diagnosing and investigating of rare diseases, with limitations due to the nature of short-read lengths. In recent years, long read sequencing (LRS) technologies have emerged as a valuable tool in overcoming these limitations. This minireview provides a concise overview of the applications of LRS in rare disease research and diagnosis, including the identification of disease-causing tandem repeat expansions, structural variations, and comprehensive analysis of pathogenic variants with LRS.

Long-read sequencing revealing intragenic deletions in exome-negative spastic paraplegias

Hereditary spastic paraplegias (HSPs) are a heterogeneous group of neurodegenerative disorders characterized by progressive spasticity and weakness in the lower extremities. To date, a total of 88 types of SPG are known. To diagnose HSP, multiple technologies, including microarray, direct sequencing, multiplex ligation-dependent probe amplification, and short-read next-generation sequencing, are often chosen based on the frequency of HSP subtypes. Exome sequencing (ES) is commonly used. We used ES to analyze ten cases of HSP from eight families. We identified pathogenic variants in three cases (from three different families); however, we were unable to determine the cause of the other seven cases using ES. We therefore applied long-read sequencing to the seven undetermined HSP cases (from five families). We detected intragenic deletions within the SPAST gene in four families, and a deletion within PSEN1 in the remaining family. The size of the deletion ranged from 4.7 to 12.5 kb and involved 1-7 exons. All deletions were entirely included in one long read. We retrospectively performed an ES-based copy number variation analysis focusing on pathogenic deletions, but were not able to accurately detect these deletions. This study demonstrated the efficiency of long-read sequencing in detecting intragenic pathogenic deletions in ES-negative HSP patients.

Long-Read DNA Sequencing: Recent Advances and Remaining Challenges

DNA sequencing has revolutionized medicine over recent decades. However, analysis of large structural variation and repetitive DNA, a hallmark of human genomes, has been limited by short-read technology, with read lengths of 100-300 bp. Long-read sequencing (LRS) permits routine sequencing of human DNA fragments tens to hundreds of kilobase pairs in size, using both real-time sequencing by synthesis and nanopore-based direct electronic sequencing. LRS permits analysis of large structural variation and haplotypic phasing in human genomes and has enabled the discovery and characterization of rare pathogenic structural variants and repeat expansions. It has also recently enabled the assembly of a complete, gapless human genome that includes previously intractable regions, such as highly repetitive centromeres and homologous acrocentric short arms. With the addition of protocols for targeted enrichment, direct epigenetic DNA modification detection, and long-range chromatin profiling, LRS promises to launch a new era of understanding of genetic diversity and pathogenic mutations in human populations.

Biallelic structural variations within FGF12 detected by long-read sequencing in epilepsy

We discovered biallelic intragenic structural variations (SVs) in FGF12 by applying long-read whole genome sequencing to an exome-negative patient with developmental and epileptic encephalopathy (DEE). We also found another DEE patient carrying a biallelic (homozygous) single-nucleotide variant (SNV) in FGF12 that was detected by exome sequencing. FGF12 heterozygous recurrent missense variants with gain-of-function or heterozygous entire duplication of FGF12 are known causes of epilepsy, but biallelic SNVs/SVs have never been described. FGF12 encodes intracellular proteins interacting with the C-terminal domain of the alpha subunit of voltage-gated sodium channels 1.2, 1.5, and 1.6, promoting excitability by delaying fast inactivation of the channels. To validate the molecular pathomechanisms of these biallelic FGF12 SVs/SNV, highly sensitive gene expression analyses using lymphoblastoid cells from the patient with biallelic SVs, structural considerations, and Drosophila in vivo functional analysis of the SNV were performed, confirming loss-of-function. Our study highlights the importance of small SVs in Mendelian disorders, which may be overlooked by exome sequencing but can be detected efficiently by long-read whole genome sequencing, providing new insights into the pathomechanisms of human diseases.

Deciphering the exact breakpoints of structural variations using long sequencing reads with DeBreak

Long-read sequencing has demonstrated great potential for characterizing all types of structural variations (SVs). However, existing algorithms have insufficient sensitivity and precision. To address these limitations, we present DeBreak, a computational method for comprehensive and accurate SV discovery. Based on alignment results, DeBreak employs a density-based approach for clustering SV candidates together with a local de novo assembly approach for reconstructing long insertions. A partial order alignment algorithm ensures precise SV breakpoints with single base-pair resolution, and a k-means clustering method can report multi-allele SV events. DeBreak outperforms existing tools on both simulated and real long-read sequencing data from both PacBio and Nanopore platforms. An important application of DeBreak is analyzing cancer genomes for potentially tumor-driving SVs. DeBreak can also be used for supplementing whole-genome assembly-based SV discovery.

Advances in sequencing technologies for amyotrophic lateral sclerosis research

Amyotrophic lateral sclerosis (ALS) is caused by upper and lower motor neuron loss and has a fairly rapid disease progression, leading to fatality in an average of 2-5 years after symptom onset. Numerous genes have been implicated in this disease; however, many cases remain unexplained. Several technologies are being used to identify regions of interest and investigate candidate genes. Initial approaches to detect ALS genes include, among others, linkage analysis, Sanger sequencing, and genome-wide association studies. More recently, next-generation sequencing methods, such as whole-exome and whole-genome sequencing, have been introduced. While those methods have been particularly useful in discovering new ALS-linked genes, methodological advances are becoming increasingly important, especially given the complex genetics of ALS. Novel sequencing technologies, like long-read sequencing, are beginning to be used to uncover the contribution of repeat expansions and other types of structural variation, which may help explain missing heritability in ALS. In this review, we discuss how popular and/or upcoming methods are being used to discover ALS genes, highlighting emerging long-read sequencing platforms and their role in aiding our understanding of this challenging disease.

A companion to the preclinical common data elements for genomics, transcriptomics, and epigenomics data in rodent epilepsy models. A report of the TASK3-WG4 omics working group of the ILAE/AES joint translational TASK force

The International League Against Epilepsy/American Epilepsy Society (ILAE/AES) Joint Translational Task Force established the TASK3 working groups to create common data elements (CDEs) for various preclinical epilepsy research disciplines. The aim of the CDEs is to improve the standardization of experimental designs across a range of epilepsy research-related methods. Here, we have generated CDE tables with key parameters and case report forms (CRFs) containing the essential contents of the study protocols for genomics, transcriptomics, and epigenomics in rodent models of epilepsy, with a specific focus on adult rats and mice. We discuss the important elements that need to be considered for genomics, transcriptomics, and epigenomics methodologies, providing a rationale for the parameters that should be collected. This is the first in a two-part series of omics papers with the second installment to cover proteomics, lipidomics, and metabolomics in adult rodents.

Clinical and genetic characterization of a cohort of 97 CLN6 patients tested at a single center

Background

Ceroid lipofuscinoses neuronal 6 (CLN6) disease belongs to the neuronal ceroid lipofuscinoses (NCLs), complex and genetically heterogeneous disorders with wide geographical and phenotypic variation. The first clinical signs usually appear between 18 months and 8 years, but examples of later-onset have also been reported. Common manifestations include ataxia, seizures, vision impairment, and developmental regression. Because these are shared by other neurological diseases, identification of CLN6 genetic variants is imperative for early diagnosis.

Results

We present one of the largest cohorts to date of genetically diagnosed CLN6 patients screened at a single center. In total 97 subjects, originating from 20 countries were screened between 2010 and 2020. They comprised 86 late-infantile, eight juvenile, and three adult-onset cases (two patients with Kufs disease type A, and one with teenage progressive myoclonic epilepsy). The male to female ratio was 1.06: 1.00. The age at referral was between six months and 33 years. The time from disease onset to referral ranged from less than 1 month to 8.3 years. The clinical phenotype consisted of a combination of symptoms, as reported before. We characterized a total of 45 distinct variants defining 45 distinct genotypes. Twenty-four were novel variants, some with distinct geographic associations. Remarkably, c.257A > G (p.H86R) was present in five out of 23 unrelated Egyptian individuals but in no patients from other countries. The most common genotype was homozygosity for the c.794_796del in-frame deletion. It was present in about one-third of CLN6 patients (28 unrelated cases, and 2 familial cases), all with late-infantile onset. Variants with a high likelihood of causing loss of CLN6 function were found in 21% of cases and made up 33% of all distinct variants. Forty-four percent of variants were classified as pathogenic or likely pathogenic.

Conclusions

Our study significantly expands the number of published clinical cases and the mutational spectrum of disease-associated CLN6 variants, especially for the Middle Eastern and North African regions. We confirm previous observations regarding the most prevalent symptoms and recommend including CLN6 in the genetic diagnosis of patients presenting with early-onset abnormalities of the nervous system, musculoskeletal system, and eye.

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 .

Comprehensive evaluation of structural variant genotyping methods based on long-read sequencing data

Background

Structural variants (SVs) play a crucial role in gene regulation, trait association, and disease in humans. SV genotyping has been extensively applied in genomics research and clinical diagnosis. Although a growing number of SV genotyping methods for long reads have been developed, a comprehensive performance assessment of these methods has yet to be done.

Results

Based on one simulated and three real SV datasets, we performed an in-depth evaluation of five SV genotyping methods, including cuteSV, LRcaller, Sniffles, SVJedi, and VaPoR. The results show that for insertions and deletions, cuteSV and LRcaller have similar F1 scores (cuteSV, insertions: 0.69–0.90, deletions: 0.77–0.90 and LRcaller, insertions: 0.67–0.87, deletions: 0.74–0.91) and are superior to other methods. For duplications, inversions, and translocations, LRcaller yields the most accurate genotyping results (0.84, 0.68, and 0.47, respectively). When genotyping SVs located in tandem repeat region or with imprecise breakpoints, cuteSV (insertions and deletions) and LRcaller (duplications, inversions, and translocations) are better than other methods. In addition, we observed a decrease in F1 scores when the SV size increased. Finally, our analyses suggest that the F1 scores of these methods reach the point of diminishing returns at 20× depth of coverage.

Conclusions

We present an in-depth benchmark study of long-read SV genotyping methods. Our results highlight the advantages and disadvantages of each genotyping method, which provide practical guidance for optimal application selection and prospective directions for tool improvement.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-022-08548-y.

Uncovering Essential Tremor Genetics: The Promise of Long-Read Sequencing

Long-read sequencing (LRS) technologies have been recently introduced to overcome intrinsic limitations of widely-used next-generation sequencing (NGS) technologies, namely the sequencing limited to short-read fragments (150–300 base pairs). Since its introduction, LRS has permitted many successes in unraveling hidden mutational mechanisms. One area in clinical neurology in need of rethinking as it applies to genetic mechanisms is essential tremor (ET). This disorder, among the most common in neurology, is a syndrome often exhibiting an autosomal dominant pattern of inheritance whose large phenotypic spectrum suggest a multitude of genetic etiologies. Exome sequencing has revealed the genetic etiology only in rare ET families (FUS, SORT1, SCN4A, NOS3, KCNS2, HAPLN4/BRAL2, and USP46). We hypothesize that a reason for this shortcoming may be non-classical genetic mechanism(s) underpinning ET, among them trinucleotide, tetranucleotide, or pentanucleotide repeat disorders. In support of this hypothesis, trinucleotide (e.g., GGC repeats in NOTCH2NLC) and pentanucleotide repeat disorders (e.g., ATTTC repeats in STARD7) have been revealed as pathogenic in patients with a past history of what has come to be referred to as “ET plus,” bilateral hand tremor associated with epilepsy and/or leukoencephalopathy. A systematic review of LRS in neurodegenerative disorders showed that 10 of the 22 (45%) genetic etiologies ascertained by LRS include tremor in their phenotypic spectrum, suggesting that future clinical applications of LRS for tremor disorders may uncover genetic subtypes of familial ET that have eluded NGS, particularly those with associated leukoencephalopathy or family history of epilepsy. LRS provides a pathway for potentially uncovering novel genes and genetic mechanisms, helping narrow the large proportion of “idiopathic” ET.

Alu element in the RNA binding motif protein, X-linked 2 (RBMX2) gene found to be linked to bipolar disorder

Objective

We have used long-read single molecule, real-time (SMRT) sequencing to fully characterize a ~12Mb genomic region on chromosome Xq24-q27, significantly linked to bipolar disorder (BD) in an extended family from a genetic sub-isolate. This family segregates BD in at least four generations with 24 affected individuals.

Methods

We selected 16 family members for targeted sequencing. The selected individuals either carried the disease haplotype, were non-carriers of the disease haplotype, or served as married-in controls. We designed hybrid capture probes enriching for 5-9Kb fragments spanning the entire 12Mb region that were then sequenced to screen for candidate structural variants (SVs) that could explain the increased risk for BD in this extended family.

Results

Altogether, 201 variants were detected in the critically linked region. Although most of these represented common variants, three variants emerged that showed near-perfect segregation among all BD type I affected individuals. Two of the SVs were identified in or near genes belonging to the RNA Binding Motif Protein, X-Linked (RBMX) gene family—a 330bp Alu (subfamily AluYa5) deletion in intron 3 of the RBMX2 gene and an intergenic 27bp tandem repeat deletion between the RBMX and G protein-coupled receptor 101 (GPR101) genes. The third SV was a 50bp tandem repeat insertion in intron 1 of the Coagulation Factor IX (F9) gene.

Conclusions

Among the three genetically linked SVs, additional evidence supported the Alu element deletion in RBMX2 as the leading candidate for contributing directly to the disease development of BD type I in this extended family.

Structural variants in the Chinese population and their impact on phenotypes, diseases and population adaptation

A complete characterization of genetic variation is a fundamental goal of human genome research. Long-read sequencing has improved the sensitivity of structural variant discovery. Here, we conduct the long-read sequencing-based structural variant analysis for 405 unrelated Chinese individuals, with 68 phenotypic and clinical measurements. We discover a landscape of 132,312 nonredundant structural variants, of which 45.2% are novel. The identified structural variants are of high-quality, with an estimated false discovery rate of 3.2%. The concatenated length of all the structural variants is approximately 13.2% of the human reference genome. We annotate 1,929 loss-of-function structural variants affecting the coding sequence of 1,681 genes. We discover rare deletions in HBA1/HBA2/HBB associated with anemia. Furthermore, we identify structural variants related to immunity which differentiate the northern and southern Chinese populations. Our study describes the landscape of structural variants in the Chinese population and their contribution to phenotypes and disease.

Clinical and molecular spectrum of P/Q type calcium channel Cav2.1 in epileptic patients

Background

Epilepsy is a neurological disorder characterized by the potential to induce seizure and accompanied by cognitive, psychological, and social consequences. CACNA1A gene is a voltage-gated P/Q-type Cav2.1 channel that is broadly expressed in the central nervous system, and the pathogenic variants within this gene may be associated with the epileptic phenotype. In the present study, we collected clinical and molecular data related to epileptic patients with CACNA1A pathogenic variants and investigated possible meaningful relationship between age at onset, neurodevelopmental disorders, type of seizures, brain imaging abnormalities, genotype, and protein domains.

Results

In our retrospective literature studies, from among 890 articles reviewed, a total of 90 individuals were related to epilepsy phenotype. Our findings showed that about 90 percent of patients have shown the first symptoms in childhood and teenage years and different types of neurodevelopmental disorders, such as intellectual disability, developmental arrest, and behavioral disorders, have been common findings for these patients. Further, a wide range of abnormalities have been observed in their brain imaging, and generalized seizures have been the most type of seizures in these patients. However, our data showed no specific genotype–phenotype correlation in epileptic patients with CACNA1A pathogenic alterations.

Conclusions

Our study focused on epileptic phenotype in patients with CACNA1A pathogenic variants and showed a wide range of clinical and molecular heterogeneity with no specific genotype–phenotype correlation. It seems that incomplete penetrance, de-novo variants, and modifier genes are obstacles in predicting the clinical outcome.

Supplementary Information

The online version contains supplementary material available at 10.1186/s13023-021-02101-y.

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.

p.Asn77Lys homozygous CLN6 mutation in two unrelated Japanese patients with Kufs disease, an adult onset neuronal ceroid lipofuscinosis

BACKGROUND

The neuronal ceroid lipofuscinosis (NCL) are a group of autosomal recessive neurodegenerative disorders that are characterized by the accumulation of ceroid lipofuscins. The NCLs are categorized into four classes based on the age of onset. Kufs disease is a rare adult-onset NCL caused by mutations in the CLN6 gene, which is rarely observed in the Japanese population.

CASE

We previously reported a case study on a patient with Kufs disease, whose parents had a consanguineous marriage. Later, we observed another unrelated patient with Kufs. Here we present the case and mutational gene report in patients with Kufs disease.

CONCLUSIONS

Gene analysis results of the first patient revealed a homozygous mutation c231C > G, p.Asn77Lys in exon 3 and a homozygous c.297 + 48 A > T mutation in intron 3 in the CLN6 gene. The Asn amino acid is perfectly conserved among species. In silico analysis showed that the mutation is predicted to be probably damaging. Moreover, the second patient with Kufs disease also had the same homozygous mutations. These data suggest that the missense mutation must be pathogenic. Furthermore, the patients had lived in the same district; therefore, they both potentially inherited the founder effect mutations.

Whole-genome sequencing as a first-tier diagnostic framework for rare genetic diseases

Rare diseases affect nearly 300 million people globally with most patients aged five or less. Traditional diagnostic approaches have provided much of the diagnosis; however, there are limitations. For instance, simply inadequate and untimely diagnosis adversely affects both the patient and their families. This review advocates the use of whole genome sequencing in clinical settings for diagnosis of rare genetic diseases by showcasing five case studies. These examples specifically describe the utilization of whole genome sequencing, which helped in providing relief to patients via correct diagnosis followed by use of precision medicine.

Deciphering complex genome rearrangements in C. elegans using short-read whole genome sequencing

Genomic rearrangements cause congenital disorders, cancer, and complex diseases in human. Yet, they are still understudied in rare diseases because their detection is challenging, despite the advent of whole genome sequencing (WGS) technologies. Short-read (srWGS) and long-read WGS approaches are regularly compared, and the latter is commonly recommended in studies focusing on genomic rearrangements. However, srWGS is currently the most economical, accurate, and widely supported technology. In Caenorhabditis elegans (C. elegans), such variants, induced by various mutagenesis processes, have been used for decades to balance large genomic regions by preventing chromosomal crossover events and allowing the maintenance of lethal mutations. Interestingly, those chromosomal rearrangements have rarely been characterized on a molecular level. To evaluate the ability of srWGS to detect various types of complex genomic rearrangements, we sequenced three balancer strains using short-read Illumina technology. As we experimentally validated the breakpoints uncovered by srWGS, we showed that, by combining several types of analyses, srWGS enables the detection of a reciprocal translocation (eT1), a free duplication (sDp3), a large deletion (sC4), and chromoanagenesis events. Thus, applying srWGS to decipher real complex genomic rearrangements in model organisms may help designing efficient bioinformatics pipelines with systematic detection of complex rearrangements in human genomes.

Optical genome mapping identifies a germline retrotransposon insertion in SMARCB1 in two siblings with atypical teratoid rhabdoid tumors

In a subset of pediatric cancers, a germline cancer predisposition is highly suspected based on clinical and pathological findings, but genetic evidence is lacking, which hampers genetic counseling and predictive testing in the families involved. We describe a family with two siblings born from healthy parents who were both neonatally diagnosed with atypical teratoid rhabdoid tumor (ATRT). This rare and aggressive pediatric tumor is associated with biallelic inactivation of SMARCB1, and in 30% of the cases, a predisposing germline mutation is involved. Whereas the tumors of both siblings showed loss of expression of SMARCB1 and acquired homozygosity of the locus, whole exome and whole genome sequencing failed to identify germline or somatic SMARCB1 pathogenic mutations. We therefore hypothesized that the insertion of a pathogenic repeat‐rich structure might hamper its detection, and we performed optical genome mapping (OGM) as an alternative strategy to identify structural variation in this locus. Using this approach, an insertion of ~2.8 kb within intron 2 of SMARCB1 was detected. Long‐range PCR covering this region remained unsuccessful, but PacBio HiFi genome sequencing identified this insertion to be a SINE‐VNTR‐Alu, subfamily E (SVA‐E) retrotransposon element, which was present in a mosaic state in the mother. This SVA‐E insertion disrupts correct splicing of the gene, resulting in loss of a functional allele. This case demonstrates the power of OGM and long‐read sequencing to identify genomic variations in high‐risk cancer‐predisposing genes that are refractory to detection with standard techniques, thereby completing the clinical and molecular diagnosis of such complex cases and greatly improving counseling and surveillance of the families involved. © 2021 The Authors. The Journal of Pathology published by John Wiley & Sons, Ltd. on behalf of The Pathological Society of Great Britain and Ireland.

Identification and Characterization of Splicing Defects by Single-Molecule Real-Time Sequencing Technology (PacBio)

Although DNA-sequencing is the most effective procedure to achieve a molecular diagnosis in genetic diseases, complementary RNA analyses are often required.Reverse-Transcription polymerase chain reaction (RT-PCR) is still a valuable option when the clinical phenotype and/or available DNA-test results address the diagnosis toward a gene of interest or when the splicing effect of a single variant needs to be assessed.We use Single-Molecule Real-Time sequencing to detect and characterize splicing defects and single nucleotide variants in well-known disease genes (DMD, NF1, TTN). After proper optimization, the procedure could be used in the diagnostic setting, simplifying the workflow of cDNA analysis.

Comparative Analysis for the Performance of Long-Read-Based Structural Variation Detection Pipelines in Tandem Repeat Regions

There has been growing recognition of the vital links between structural variations (SVs) and diverse diseases. Research suggests that, with much longer DNA fragments and abundant contextual information, long-read technologies have advantages in SV detection even in complex repetitive regions. So far, several pipelines for calling SVs from long-read sequencing data have been proposed and used in human genome research. However, the performance of these pipelines is still lack of deep exploration and adequate comparison. In this study, we comprehensively evaluated the performance of three commonly used long-read SV detection pipelines, namely PBSV, Sniffles and PBHoney, especially the performance on detecting the SVs in tandem repeat regions (TRRs). Evaluated by using a robust benchmark for germline SV detection as the gold standard, we thoroughly estimated the precision, recall and F1 score of insertions and deletions detected by the pipelines. Our results revealed that all these pipelines clearly exhibited better performance outside TRRs than that in TRRs. The F1 scores of Sniffles in and outside TRRs were 0.60 and 0.76, respectively. The performance of PBSV was similar to that of Sniffles, and was generally higher than that of PBHoney. In conclusion, our findings can be benefit for choosing the appropriate pipelines in real practice and are good complementary to the application of long-read sequencing technologies in the research of rare diseases.

Complete sequencing of expanded SAMD12 repeats by long-read sequencing and Cas9-mediated enrichment

A pentanucleotide TTTCA repeat insertion into a polymorphic TTTTA repeat element in SAMD12 causes benign adult familial myoclonic epilepsy. Although the precise determination of the entire SAMD12 repeat sequence is important for molecular diagnosis and research, obtaining this sequence remains challenging when using conventional genomic/genetic methods, and even short-read and long-read next-generation sequencing technologies have been insufficient. Incomplete information regarding expanded repeat sequences may hamper our understanding of the pathogenic roles played by varying numbers of repeat units, genotype-phenotype correlations, and mutational mechanisms. Here, we report a new approach for the precise determination of the entire expanded repeat sequence and present a workflow designed to improve the diagnostic rates in various repeat expansion diseases. We examined 34 clinically diagnosed benign adult familial myoclonic epilepsy patients, from 29 families using repeat-primed PCR, Southern blot, and long-read sequencing with Cas9-mediated enrichment. Two cases with questionable results from repeat-primed PCR and/or Southern blot were confirmed as pathogenic using long-read sequencing with Cas9-mediated enrichment, resulting in the identification of pathogenic SAMD12 repeat expansions in 76% of examined families (22/29). Importantly, long-read sequencing with Cas9-mediated enrichment was able to provide detailed information regarding the sizes, configurations, and compositions of the expanded repeats. The inserted TTTCA repeat size and the proportion of TTTCA sequences among the overall repeat sequences were highly variable, and a novel repeat configuration was identified. A genotype-phenotype correlation study suggested that the insertion of even short (TTTCA)14 repeats contributed to the development of benign adult familial myoclonic epilepsy. However, the sizes of the overall TTTTA and TTTCA repeat units are also likely to be involved in the pathology of benign adult familial myoclonic epilepsy. Seven unsolved SAMD12-negative cases were investigated using whole-genome long-read sequencing, and infrequent, disease-associated, repeat expansions were identified in two cases. The strategic workflow resolved two questionable SAMD12-positive cases and two previously SAMD12-negative cases, increasing the diagnostic yield from 69% (20/29 families) to 83% (24/29 families). This study indicates the significant utility of long-read sequencing technologies to explore the pathogenic contributions made by various repeat units in complex repeat expansions and to improve the overall diagnostic rate.

The genetics of mitochondrial disease: dissecting mitochondrial pathology using multi-omic pipelines

Mitochondria play essential roles in numerous metabolic pathways including the synthesis of adenosine triphosphate through oxidative phosphorylation. Clinically, mitochondrial diseases occur when there is mitochondrial dysfunction – manifesting at any age and affecting any organ system; tissues with high energy requirements, such as muscle and the brain, are often affected. The clinical heterogeneity is parallel to the degree of genetic heterogeneity associated with mitochondrial dysfunction. Around 10% of human genes are predicted to have a mitochondrial function, and defects in over 300 genes are reported to cause mitochondrial disease. Some involve the mitochondrial genome (mtDNA), but the vast majority occur within the nuclear genome. Except for a few specific genetic defects, there remains no cure for mitochondrial diseases, which means that a genetic diagnosis is imperative for genetic counselling and the provision of reproductive options for at‐risk families. Next‐generation sequencing strategies, particularly exome and whole‐genome sequencing, have revolutionised mitochondrial diagnostics such that the traditional muscle biopsy has largely been replaced with a minimally‐invasive blood sample for an unbiased approach to genetic diagnosis. Where these genomic approaches have not identified a causative defect, or where there is insufficient support for pathogenicity, additional functional investigations are required. The application of supplementary ‘omics’ technologies, including transcriptomics, proteomics, and metabolomics, has the potential to greatly improve diagnostic strategies. This review aims to demonstrate that whilst a molecular diagnosis can be achieved for many cases through next‐generation sequencing of blood DNA, the use of patient tissues and an integrated, multidisciplinary multi‐omics approach is pivotal for the diagnosis of more challenging cases. Moreover, the analysis of clinically relevant tissues from affected individuals remains crucial for understanding the molecular mechanisms underlying mitochondrial pathology. © 2021 The Authors. The Journal of Pathology published by John Wiley & Sons, Ltd. on behalf of The Pathological Society of Great Britain and Ireland.

Panel-based genetic testing for inherited retinal disease screening 176 genes

Background

This case series reports the performance of a next‐generation sequencing (NGS) panel of 176 retinal genes (NGS 176) in patients with inherited retinal disease (IRD).

Methods

Subjects are patients who underwent genetic testing between 1 August 2016 and 1 January 2018 at Moorfields Eye Hospital, London, UK. Panel‐based genetic testing was performed unless a specific gene (e.g., RS1) or small group of genes (e.g., ABCA4, PRPH2) were suspected. If a novel variant was identified, a further comment on their predicted pathogenicity and evolutionary conservation was offered and segregation studies performed. The main outcome measure is the likelihood of obtaining a genetic diagnosis using NGS 176.

Results

488 patients were included. A molecular diagnosis was obtained for 59.4% of patients. Younger patients were more likely to receive a molecular diagnosis; with 92% of children under the age of 6 years receiving a conclusive result. There was a change in their initially assigned inheritance pattern in 8.4% of patients following genetic testing. Selected IRD diagnoses (e.g., achromatopsia, congenital stationary night blindness) were associated with high diagnostic yields.

Conclusion

This study confirms that NGS 176 is a useful first‐tier genetic test for most IRD patients. Age and initial clinical diagnosis were strongly associated with diagnostic yield.

Long-Read Sequencing Improves the Detection of Structural Variations Impacting Complex Non-Coding Elements of the Genome

The advent of long-read sequencing offers a new assessment method of detecting genomic structural variation (SV) in numerous rare genetic diseases. For autism spectrum disorders (ASD) cases where pathogenic variants fail to be found in the protein-coding genic regions along chromosomes, we proposed a scalable workflow to characterize the risk factor of SVs impacting non-coding elements of the genome. We applied whole-genome sequencing on an Emirati family having three children with ASD using long and short-read sequencing technology. A series of analytical pipelines were established to identify a set of SVs with high sensitivity and specificity. At 15-fold coverage, we observed that long-read sequencing technology (987 variants) detected a significantly higher number of SVs when compared to variants detected using short-read technology (509 variants) (p-value < 1.1020 × 10-57). Further comparison showed 97.9% of long-read sequencing variants were spanning within the 1-100 kb size range (p-value < 9.080 × 10-67) and impacting over 5000 genes. Moreover, long-read variants detected 604 non-coding RNAs (p-value < 9.02 × 10-9), comprising 58% microRNA, 31.9% lncRNA, and 9.1% snoRNA. Even at low coverage, long-read sequencing has shown to be a reliable technology in detecting SVs impacting complex elements of the genome.

Long-read whole-genome sequencing identified a partial MBD5 deletion in an exome-negative patient with neurodevelopmental disorder

Whole-exome sequencing (WES) can detect not only single-nucleotide variants in causal genes, but also pathogenic copy-number variations using several methods. However, there may be overlooked pathogenic variations in the out of target genome regions of WES analysis (e.g., promoters), leaving many patients undiagnosed. Whole-genome sequencing (WGS) can potentially analyze such regions. We applied long-read nanopore WGS and our recently developed analysis pipeline "dnarrange" to a patient who was undiagnosed by trio-based WES analysis, and identified a heterozygous 97-kb deletion partially involving 5'-untranslated exons of MBD5, which was outside the WES target regions. The phenotype of the patient, a 32-year-old male, was consistent with haploinsufficiency of MBD5. The transcript level of MBD5 in the patient's lymphoblastoid cells was reduced. We therefore concluded that the partial MBD5 deletion is the culprit for this patient. Furthermore, we found other rare structural variations (SVs) in this patient, i.e., a large inversion and a retrotransposon insertion, which were not seen in 33 controls. Although we considered that they are benign SVs, this finding suggests that our pipeline using long-read WGS is useful for investigating various types of potentially pathogenic SVs. In conclusion, we identified a 97-kb deletion, which causes haploinsufficiency of MBD5 in a patient with neurodevelopmental disorder, demonstrating that long-read WGS is a powerful technique to discover pathogenic SVs.

Pathogenic 12-kb copy-neutral inversion in syndromic intellectual disability identified by high-fidelity long-read sequencing

We report monozygotic twin girls with syndromic intellectual disability who underwent exome sequencing but with negative pathogenic variants. To search for variants that are unrecognized by exome sequencing, high-fidelity long-read genome sequencing (HiFi LR-GS) was applied. A 12-kb copy-neutral inversion was precisely identified by HiFi LR-GS after trio-based variant filtering. This inversion directly disrupted two genes, CPNE9 and BRPF1, the latter of which attracted our attention because pathogenic BRPF1 variants have been identified in autosomal dominant intellectual developmental disorder with dysmorphic facies and ptosis (IDDDFP), which later turned out to be clinically found in the twins. Trio-based HiFi LR-GS together with haplotype phasing revealed that the 12-kb inversion occurred de novo on the maternally transmitted chromosome. This study clearly indicates that submicroscopic copy-neutral inversions are important but often uncharacterized culprits in monogenic disorders and that long-read sequencing is highly advantageous for detecting such inversions involved in genetic diseases.

Assembling and Validating Bioinformatic Pipelines for Next-Generation Sequencing Clinical Assays

CONTEXT.—

Clinical next-generation sequencing (NGS) is being rapidly adopted, but analysis and interpretation of large data sets prompt new challenges for a clinical laboratory setting. Clinical NGS results rely heavily on the bioinformatics pipeline for identifying genetic variation in complex samples. The choice of bioinformatics algorithms, genome assembly, and genetic annotation databases are important for determining genetic alterations associated with disease. The analysis methods are often tuned to the assay to maximize accuracy. Once a pipeline has been developed, it must be validated to determine accuracy and reproducibility for samples similar to real-world cases. In silico proficiency testing or institutional data exchange will ensure consistency among clinical laboratories.

OBJECTIVE.—

To provide molecular pathologists a step-by-step guide to bioinformatics analysis and validation design in order to navigate the regulatory and validation standards of implementing a bioinformatic pipeline as a part of a new clinical NGS assay.

DATA SOURCES.—

This guide uses published studies on genomic analysis, bioinformatics methods, and methods comparison studies to inform the reader on what resources, including open source software tools and databases, are available for genetic variant detection and interpretation.

CONCLUSIONS.—

This review covers 4 key concepts: (1) bioinformatic analysis design for detecting genetic variation, (2) the resources for assessing genetic effects, (3) analysis validation assessment experiments and data sets, including a diverse set of samples to mimic real-world challenges that assess accuracy and reproducibility, and (4) if concordance between clinical laboratories will be improved by proficiency testing designed to test bioinformatic pipelines.

Long-read human genome sequencing and its applications

Over the past decade, long-read, single-molecule DNA sequencing technologies have emerged as powerful players in genomics. With the ability to generate reads tens to thousands of kilobases in length with an accuracy approaching that of short-read sequencing technologies, these platforms have proven their ability to resolve some of the most challenging regions of the human genome, detect previously inaccessible structural variants and generate some of the first telomere-to-telomere assemblies of whole chromosomes. Long-read sequencing technologies will soon permit the routine assembly of diploid genomes, which will revolutionize genomics by revealing the full spectrum of human genetic variation, resolving some of the missing heritability and leading to the discovery of novel mechanisms of disease.

Long-read sequencing for non-small-cell lung cancer genomes

Here, we report the application of a long-read sequencer, PromethION, for analyzing human cancer genomes. We first conducted whole-genome sequencing on lung cancer cell lines. We found that it is possible to genotype known cancerous mutations, such as point mutations. We also found that long-read sequencing is particularly useful for precisely identifying and characterizing structural aberrations, such as large deletions, gene fusions, and other chromosomal rearrangements. In addition, we identified several medium-sized structural aberrations consisting of complex combinations of local duplications, inversions, and microdeletions. These complex mutations occurred even in key cancer-related genes, such as STK11, NF1, SMARCA4, and PTEN. The biological relevance of those mutations was further revealed by epigenome, transcriptome, and protein analyses of the affected signaling pathways. Such structural aberrations were also found in clinical lung adenocarcinoma specimens. Those structural aberrations were unlikely to be reliably detected by conventional short-read sequencing. Therefore, long-read sequencing may contribute to understanding the molecular etiology of patients for whom causative cancerous mutations remain unknown and therapeutic strategies are elusive.

New Diagnostic Approaches for Undiagnosed Rare Genetic Diseases

Accurate diagnosis is the cornerstone of medicine; it is essential for informed care and promoting patient and family well-being. However, families with a rare genetic disease (RGD) often spend more than five years on a diagnostic odyssey of specialist visits and invasive testing that is lengthy, costly, and often futile, as 50% of patients do not receive a molecular diagnosis. The current diagnostic paradigm is not well designed for RGDs, especially for patients who remain undiagnosed after the initial set of investigations, and thus requires an expansion of approaches in the clinic. Leveraging opportunities to participate in research programs that utilize new technologies to understand RGDs is an important path forward for patients seeking a diagnosis. Given recent advancements in such technologies and international initiatives, the prospect of identifying a molecular diagnosis for all patients with RGDs has never been so attainable, but achieving this goal will require global cooperation at an unprecedented scale.

An elongated tract of polyQ in the carboxyl‑terminus of human α1A calcium channel induces cell apoptosis by nuclear translocation

An aberrant elongated tract of glutamine residues (polyQ) in proteins induces multiple diseases treated in the clinic. In our previous study of progressive myoclonic epilepsy (PME), using whole-exome sequencing, a mutant Cav2.1 protein with an aberrant elongated polyQ tract was identified in PME patients. To investigate the molecular mechanism and cell biology of this aberrant elongated polyQ tract, wild-type Cav2.1 with 13 polyQ repeats (Cav2.1 wt-Q13) and mutant-type Cav2.1 with 26 polyQ repeats (Cav2.1 mt-Q26) were prepared and introduced into human SH-SY5Y neuroblastoma cells. Using a WST-1 assay, it was revealed that Cav2.1 mt-Q26 markedly suppressed the proliferation of the SH-SY5Y cells, a result not observed for the Cav2.1 wt-Q13-transfected cells. It was also revealed that Cav2.1 mt and its truncated molecules suppressed cell proliferation by inducing apoptosis rather than arresting the cell cycle. Further investigations indicated a nuclear translocation phenomenon associated with the Cav2.1 mt molecules. Mechanistically, it was revealed that the Cav2.1 mt molecules activated the Bcl-2/Bax, caspase-3 and poly ADP-ribose polymerase (PARP) apoptotic pathways. The present study may provide new insights for interpreting the pathogenesis of PME and the relationship among polyQ, CACNA1A gene mutations and PME.

The third generation sequencing: the advanced approach to genetic diseases

Genomic sequencing technologies have revolutionized mutation detection of the genetic diseases in the past few years. In recent years, the third generation sequencing (TGS) has been gaining insight into more genetic diseases owing to the single molecular and real time sequencing technology. This paper reviews the genomic sequencing revolutionary history first and then focuses on the genetic diseases discovered through the TGS and the clinical effects of the TGS, which is followed by the discussion of the improvement in the bioinformatic analysis for the TGS and its limitations. In summary, the TGS has been enhancing the diagnostic accuracy of genetic diseases in molecular level as well as paving a new way for basic researches and therapies.

Analysis of gene variants in the GASH/Sal model of epilepsy

Epilepsy is a complex neurological disorder characterized by sudden and recurrent seizures, which are caused by various factors, including genetic abnormalities. Several animal models of epilepsy mimic the different symptoms of this disorder. In particular, the genetic audiogenic seizure hamster from Salamanca (GASH/Sal) animals exhibit sound-induced seizures similar to the generalized tonic seizures observed in epileptic patients. However, the genetic alterations underlying the audiogenic seizure susceptibility of the GASH/Sal model remain unknown. In addition, gene variations in the GASH/Sal might have a close resemblance with those described in humans with epilepsy, which is a prerequisite for any new preclinical studies that target genetic abnormalities. Here, we performed whole exome sequencing (WES) in GASH/Sal animals and their corresponding controls to identify and characterize the mutational landscape of the GASH/Sal strain. After filtering the results, moderate- and high-impact variants were validated by Sanger sequencing, assessing the possible impact of the mutations by “in silico” reconstruction of the encoded proteins and analyzing their corresponding biological pathways. Lastly, we quantified gene expression levels by RT-qPCR. In the GASH/Sal model, WES showed the presence of 342 variations, in which 21 were classified as high-impact mutations. After a full bioinformatics analysis to highlight the high quality and reliable variants, the presence of 3 high-impact and 15 moderate-impact variants were identified. Gene expression analysis of the high-impact variants of Asb14 (ankyrin repeat and SOCS Box Containing 14), Msh3 (MutS Homolog 3) and Arhgef38 (Rho Guanine Nucleotide Exchange Factor 38) genes showed a higher expression in the GASH/Sal than in control hamsters. In silico analysis of the functional consequences indicated that those mutations in the three encoded proteins would have severe functional alterations. By functional analysis of the variants, we detected 44 significantly enriched pathways, including the glutamatergic synapse pathway. The data show three high-impact mutations with a major impact on the function of the proteins encoded by these genes, although no mutation in these three genes has been associated with some type of epilepsy until now. Furthermore, GASH/Sal animals also showed gene variants associated with different types of epilepsy that has been extensively documented, as well as mutations in other genes that encode proteins with functions related to neuronal excitability, which could be implied in the phenotype of the GASH/Sal. Our findings provide valuable genetic and biological pathway data associated to the genetic burden of the audiogenic seizure susceptibility and reinforce the need to validate the role of each key mutation in the phenotype of the GASH/Sal model.

Structural variation in the sequencing era

Identifying structural variation (SV) is essential for genome interpretation but has been historically difficult due to limitations inherent to available genome technologies. Detection methods that use ensemble algorithms and emerging sequencing technologies have enabled the discovery of thousands of SVs, uncovering information about their ubiquity, relationship to disease and possible effects on biological mechanisms. Given the variability in SV type and size, along with unique detection biases of emerging genomic platforms, multiplatform discovery is necessary to resolve the full spectrum of variation. Here, we review modern approaches for investigating SVs and proffer that, moving forwards, studies integrating biological information with detection will be necessary to comprehensively understand the impact of SV in the human genome.

Genome analysis and knowledge-driven variant interpretation with TGex

BACKGROUND

The clinical genetics revolution ushers in great opportunities, accompanied by significant challenges. The fundamental mission in clinical genetics is to analyze genomes, and to identify the most relevant genetic variations underlying a patient's phenotypes and symptoms. The adoption of Whole Genome Sequencing requires novel capacities for interpretation of non-coding variants.

RESULTS

We present TGex, the Translational Genomics expert, a novel genome variation analysis and interpretation platform, with remarkable exome analysis capacities and a pioneering approach of non-coding variants interpretation. TGex's main strength is combining state-of-the-art variant filtering with knowledge-driven analysis made possible by VarElect, our highly effective gene-phenotype interpretation tool. VarElect leverages the widely used GeneCards knowledgebase, which integrates information from > 150 automatically-mined data sources. Access to such a comprehensive data compendium also facilitates TGex's broad variant annotation, supporting evidence exploration, and decision making. TGex has an interactive, user-friendly, and easy adaptive interface, ACMG compliance, and an automated reporting system. Beyond comprehensive whole exome sequence capabilities, TGex encompasses innovative non-coding variants interpretation, towards the goal of maximal exploitation of whole genome sequence analyses in the clinical genetics practice. This is enabled by GeneCards' recently developed GeneHancer, a novel integrative and fully annotated database of human enhancers and promoters. Examining use-cases from a variety of TGex users world-wide, we demonstrate its high diagnostic yields (42% for single exome and 50% for trios in 1500 rare genetic disease cases) and critical actionable genetic findings. The platform's support for integration with EHR and LIMS through dedicated APIs facilitates automated retrieval of patient data for TGex's customizable reporting engine, establishing a rapid and cost-effective workflow for an entire range of clinical genetic testing, including rare disorders, cancer predisposition, tumor biopsies and health screening.

CONCLUSIONS

TGex is an innovative tool for the annotation, analysis and prioritization of coding and non-coding genomic variants. It provides access to an extensive knowledgebase of genomic annotations, with intuitive and flexible configuration options, allows quick adaptation, and addresses various workflow requirements. It thus simplifies and accelerates variant interpretation in clinical genetics workflows, with remarkable diagnostic yield, as exemplified in the described use cases. TGex is available at http://tgex.genecards.org/.

Longshot enables accurate variant calling in diploid genomes from single-molecule long read sequencing

Whole-genome sequencing using sequencing technologies such as Illumina enables the accurate detection of small-scale variants but provides limited information about haplotypes and variants in repetitive regions of the human genome. Single-molecule sequencing (SMS) technologies such as Pacific Biosciences and Oxford Nanopore generate long reads that can potentially address the limitations of short-read sequencing. However, the high error rate of SMS reads makes it challenging to detect small-scale variants in diploid genomes. We introduce a variant calling method, Longshot, which leverages the haplotype information present in SMS reads to accurately detect and phase single-nucleotide variants (SNVs) in diploid genomes. We demonstrate that Longshot achieves very high accuracy for SNV detection using whole-genome Pacific Biosciences data, outperforms existing variant calling methods, and enables variant detection in duplicated regions of the genome that cannot be mapped using short reads.

Single-molecule sequencing (SMS) such as Pacific Biosciences and Oxford Nanopore generate long reads with high error rate. Here, the authors develop Longshot, a computational method that detects and phases single nucleotide variants (SNV) in diploid genomes using SMS data.

Long-read sequencing for rare human genetic diseases

During the past decade, the search for pathogenic mutations in rare human genetic diseases has involved huge efforts to sequence coding regions, or the entire genome, using massively parallel short-read sequencers. However, the approximate current diagnostic rate is <50% using these approaches, and there remain many rare genetic diseases with unknown cause. There may be many reasons for this, but one plausible explanation is that the responsible mutations are in regions of the genome that are difficult to sequence using conventional technologies (e.g., tandem-repeat expansion or complex chromosomal structural aberrations). Despite the drawbacks of high cost and a shortage of standard analytical methods, several studies have analyzed pathogenic changes in the genome using long-read sequencers. The results of these studies provide hope that further application of long-read sequencers to identify the causative mutations in unsolved genetic diseases may expand our understanding of the human genome and diseases. Such approaches may also be applied to molecular diagnosis and therapeutic strategies for patients with genetic diseases in the future.

Rapid progression of a walking disability in a 5-year-old boy with a CLN6 mutation

INTRODUCTION

Neuronal ceroid lipofuscinoses (NCLs; CLN) are mainly autosomal recessive neurodegenerative disorders characterized by the accumulation of autofluorescent lipopigments in neuronal and other cells. Symptoms include visual disabilities, motor decline, and epilepsy. Causative genes are CLN1, CLN2, CLN3, CLN5, CLN6, CLN7, CLN8, CLN10, CLN11, CLN12, CLN13, and CLN14. We present the fourth Japanese case with a CLN6 mutation.

CASE PRESENTATION

At 3 years of age, our patient became clumsy and fell down easily. He developed focal seizures with impaired consciousness and was started on carbamazepine. He showed ataxic walking and dysarthria with increased deep tendon reflexes. Interictal electroencephalogram revealed slow waves in the left temporal and occipital areas. Brain magnetic resonance imaging showed cerebellar atrophy and ventriculomegaly. In optical coherence tomography (OCT), the inner layer of the retina was thick and highly reflective. Exome sequencing revealed a known homozygous mutation, C.794_976del, p. (Ser265del) in CLN6.

DISCUSSION

A total of 130 cases of NCL with CLN6 mutations have been reported globally, of which only four were from Japan including the current patient. The deletion of serine at position 265 has been reported in six cases. Ser265 is located in a region of short repeated sequences that is susceptible to mutation. Clinical trials of gene therapy using adeno-associated virus serotype 9 have started for NCL6, making early diagnosis crucial. OCT examination might be helpful in achieving a diagnosis.