Evaluation of the detection of GBA missense mutations and other variants using the Oxford Nanopore MinION

Oxford nanopore sequencing-based assay for BTD gene screening: Design, clinical validation, and variant frequency assessment in the Turkish population

Biotinidase deficiency (BTD) is an autosomal recessive disorder characterized by impaired recycling of the water-soluble vitamin biotin which leads to a spectrum of clinical manifestations ranging from mild to severe, including mainly neurological and cutaneous symptoms. Biotin supplementation is a cornerstone of treatment, but diagnosis often relies on measuring serum enzyme activity, which needs to be confirmed by genetic analysis. Thus, molecular methods become necessary in the differential diagnosis of BTD. Accordingly, countries with a high-incidence have implemented next-generation sequencing (NGS) techniques to newborn screening programs for BT. Nevertheless, NGS platforms, while well-established, present challenges in cost, labor, accessibility, and duration for newborn screening programs targeting BTD, therefore these limitations necessitate the exploration of alternative systems to ensure efficient and widespread screening. Here, third-generation sequencing platforms, notably Oxford Nanopore Technology (ONT), present promising solutions to the associated challenges. Hence, in the present study, we aimed to develop an ONT-based assay for the screening of BTD gene. After designing and optimizing primers for long-PCR using reference DNA, we assessed the performance of the ONT assay in BTD patients previously diagnosed by enzyme assay and confirmed using Illumina-based sequencing. The results demonstrate a strong correlation between the two methods, indicating the reliability of the ONT-based assay. Moreover, this first in-house single gene testing specifically tailored for BTD successfully detected previously known genetic variants with high sequencing depths, affirming the effectiveness of ONT-based sequencing in human genetics.

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.

Detection of mosaic and population-level structural variants with Sniffles2

Calling structural variations (SVs) is technically challenging, but using long reads remains the most accurate way to identify complex genomic alterations. Here we present Sniffles2, which improves over current methods by implementing a repeat aware clustering coupled with a fast consensus sequence and coverage-adaptive filtering. Sniffles2 is 11.8 times faster and 29% more accurate than state-of-the-art SV callers across different coverages (5-50×), sequencing technologies (ONT and HiFi) and SV types. Furthermore, Sniffles2 solves the problem of family-level to population-level SV calling to produce fully genotyped VCF files. Across 11 probands, we accurately identified causative SVs around MECP2, including highly complex alleles with three overlapping SVs. Sniffles2 also enables the detection of mosaic SVs in bulk long-read data. As a result, we identified multiple mosaic SVs in brain tissue from a patient with multiple system atrophy. The identified SV showed a remarkable diversity within the cingulate cortex, impacting both genes involved in neuron function and repetitive elements.

Carrier frequency estimation of pathogenic variants of autosomal recessive and X-linked recessive mendelian disorders using exome sequencing data in 1,642 Thais

BACKGROUND

People with autosomal recessive disorders often were born without awareness of the carrier status of their parents. The American College of Medical Genetics and Genomics (ACMG) recommends screening 113 genes known to cause autosomal recessive and X-linked conditions in couples seeking to learn about their risk of having children with these disorders to have an appropriate reproductive plan.

METHODS

We analyzed the exome sequencing data of 1,642 unrelated Thai individuals to identify the pathogenic variant (PV) frequencies in genes recommended by ACMG.

RESULTS

In the 113 ACMG-recommended genes, 165 PV and likely PVs in 60 genes of 559 exomes (34%, 559/1642) were identified. The carrier rate was increased to 39% when glucose-6-phosphate dehydrogenase (G6PD) was added. The carrier rate was still as high as 14.7% when thalassemia and hemoglobinopathies were excluded. In addition to thalassemia, hemoglobinopathies, and G6PD deficiency, carrier frequencies of > 1% were found for Gaucher disease, primary hyperoxaluria, Pendred syndrome, and Wilson disease. Nearly 2% of the couples were at risk of having offsprings with the tested autosomal recessive conditions.

CONCLUSIONS

Based on the study samples, the expanded carrier screening, which specifically targeted common autosomal recessive conditions in Thai individuals, will benefit clinical outcomes, regarding preconception/prenatal genetic carrier screening.

Accurate long-read sequencing identified GBA1 as major risk factor in the Luxembourgish Parkinson's study

Heterozygous variants in the glucocerebrosidase GBA1 gene are an increasingly recognized risk factor for Parkinson's disease (PD). Due to the GBAP1 pseudogene, which shares 96% sequence homology with the GBA1 coding region, accurate variant calling by array-based or short-read sequencing methods remains a major challenge in understanding the genetic landscape of GBA1-associated PD. We analyzed 660 patients with PD, 100 patients with Parkinsonism and 808 healthy controls from the Luxembourg Parkinson's study, sequenced using amplicon-based long-read DNA sequencing technology. We found that 12.1% (77/637) of PD patients carried GBA1 variants, with 10.5% (67/637) of them carrying known pathogenic variants (including severe, mild, risk variants). In comparison, 5% (34/675) of the healthy controls carried GBA1 variants, and among them, 4.3% (29/675) were identified as pathogenic variant carriers. We found four GBA1 variants in patients with atypical parkinsonism. Pathogenic GBA1 variants were 2.6-fold more frequently observed in PD patients compared to controls (OR = 2.6; CI = [1.6,4.1]). Three novel variants of unknown significance (VUS) were identified. Using a structure-based approach, we defined a potential risk prediction method for VUS. This study describes the full landscape of GBA1-related parkinsonism in Luxembourg, showing a high prevalence of GBA1 variants as the major genetic risk for PD. Although the long-read DNA sequencing technique used in our study may be limited in its effectiveness to detect potential structural variants, our approach provides an important advancement for highly accurate GBA1 variant calling, which is essential for providing access to emerging causative therapies for GBA1 carriers.

GBA1 in Parkinson's disease: variant detection and pathogenicity scoring matters

BACKGROUND

GBA1 variants are the strongest genetic risk factor for Parkinson's disease (PD). However, the pathogenicity of GBA1 variants concerning PD is still not fully understood. Additionally, the frequency of GBA1 variants varies widely across populations.

OBJECTIVES

To evaluate Oxford Nanopore sequencing as a strategy, to determine the frequency of GBA1 variants in Norwegian PD patients and controls, and to review the current literature on newly identified variants that add to pathogenicity determination.

METHODS

We included 462 Norwegian PD patients and 367 healthy controls. We sequenced the full-length GBA1 gene on the Oxford Nanopore GridION as an 8.9 kb amplicon. Six analysis pipelines were compared using two aligners (NGMLR, Minimap2) and three variant callers (BCFtools, Clair3, Pepper-Margin-Deepvariant). Confirmation of GBA1 variants was performed by Sanger sequencing and the pathogenicity of variants was evaluated.

RESULTS

We found 95.8% (115/120) true-positive GBA1 variant calls, while 4.2% (5/120) variant calls were false-positive, with the NGMLR/Minimap2-BCFtools pipeline performing best. In total, 13 rare GBA1 variants were detected: two were predicted to be (likely) pathogenic and eleven were of uncertain significance. The odds of carrying one of the two common GBA1 variants, p.L483P or p.N409S, in PD patients were estimated to be 4.11 times the odds of carrying one of these variants in controls (OR = 4.11 [1.39, 12.12]).

CONCLUSIONS

In conclusion, we have demonstrated that Oxford long-read Nanopore sequencing, along with the NGMLR/Minimap2-BCFtools pipeline is an effective tool to investigate GBA1 variants. Further studies on the pathogenicity of GBA1 variants are needed to assess their effect on PD.

Structural variation across 138,134 samples in the TOPMed consortium

Ever larger Structural Variant (SV) catalogs highlighting the diversity within and between populations help researchers better understand the links between SVs and disease. The identification of SVs from DNA sequence data is non-trivial and requires a balance between comprehensiveness and precision. Here we present a catalog of 355,667 SVs (59.34% novel) across autosomes and the X chromosome (50bp+) from 138,134 individuals in the diverse TOPMed consortium. We describe our methodologies for SV inference resulting in high variant quality and >90% allele concordance compared to long-read de-novo assemblies of well-characterized control samples. We demonstrate utility through significant associations between SVs and important various cardio-metabolic and hematologic traits. We have identified 690 SV hotspots and deserts and those that potentially impact the regulation of medically relevant genes. This catalog characterizes SVs across multiple populations and will serve as a valuable tool to understand the impact of SV on disease development and progression.

Structural variation across 138,134 samples in the TOPMed consortium

Ever larger Structural Variant (SV) catalogs highlighting the diversity within and between populations help researchers better understand the links between SVs and disease. The identification of SVs from DNA sequence data is non-trivial and requires a balance between comprehensiveness and precision. Here we present a catalog of 355,667 SVs (59.34% novel) across autosomes and the X chromosome (50bp+) from 138,134 individuals in the diverse TOPMed consortium. We describe our methodologies for SV inference resulting in high variant quality and >90% allele concordance compared to long-read de-novo assemblies of well-characterized control samples. We demonstrate utility through significant associations between SVs and important various cardio-metabolic and hemotologic traits. We have identified 690 SV hotspots and deserts and those that potentially impact the regulation of medically relevant genes. This catalog characterizes SVs across multiple populations and will serve as a valuable tool to understand the impact of SV on disease development and progression.

PMS2-associated Lynch syndrome: Past, present and future

Carriers of any pathogenic variant in one of the MMR genes (path_MMR carriers) were traditionally thought to be at comparable risk of developing a range of different malignancies, foremost colorectal cancer (CRC) and endometrial cancer. However, it is now widely accepted that their cancer risk and cancer spectrum range notably depending on which MMR gene is affected. Moreover, there is increasing evidence that the MMR gene affected also influences the molecular pathogenesis of Lynch syndrome CRC. Although substantial progress has been made over the past decade in understanding these differences, many questions remain unanswered, especially pertaining to path_PMS2 carriers. Recent findings show that, while the cancer risk is relatively low, PMS2-deficient CRCs tend to show more aggressive behaviour and have a worse prognosis than other MMR-deficient CRCs. This, together with lower intratumoral immune infiltration, suggests that PMS2-deficient CRCs might have more in common biologically with sporadic MMR-proficient CRCs than with other MMR-deficient CRCs. These findings could have important consequences for surveillance, chemoprevention and therapeutic strategies (e.g. vaccines). In this review we discuss the current knowledge, current (clinical) challenges and knowledge gaps that should be targeted by future studies.

Application of Long-Read Nanopore Sequencing to the Search for Mutations in Hypertrophic Cardiomyopathy

Increasing evidence suggests that both coding and non-coding regions of sarcomeric protein genes can contribute to hypertrophic cardiomyopathy (HCM). Here, we introduce an experimental workflow (tested on four patients) for complete sequencing of the most common HCM genes (MYBPC3, MYH7, TPM1, TNNT2, and TNNI3) via long-range PCR, Oxford Nanopore Technology (ONT) sequencing, and bioinformatic analysis. We applied Illumina and Sanger sequencing to validate the results, FastQC, Qualimap, and MultiQC for quality evaluations, MiniMap2 to align data, Clair3 to call and phase variants, and Annovar's tools and CADD to assess pathogenicity of variants. We could not amplify the region encompassing exons 6-12 of MYBPC3. A higher sequencing error rate was observed with ONT (6.86-6.92%) than with Illumina technology (1.14-1.35%), mostly for small indels. Pathogenic variant p.Gln1233Ter and benign polymorphism p.Arg326Gln in MYBPC3 in a heterozygous state were found in one patient. We demonstrated the ability of ONT to phase single-nucleotide variants, enabling direct haplotype determination for genes TNNT2 and TPM1. These findings highlight the importance of long-range PCR efficiency, as well as lower accuracy of variant calling by ONT than by Illumina technology; these differences should be clarified prior to clinical application of the ONT method.

Enhancing Molecular Testing for Effective Delivery of Actionable Gene Diagnostics

There is a deep need to navigate within our genomic data to find, understand and pave the way for disease-specific treatments, as the clinical diagnostic journey provides only limited guidance. The human genome is enclosed in every nucleated cell, and yet at the single-cell resolution many unanswered questions remain, as most of the sequencing techniques use a bulk approach. Therefore, heterogeneity, mosaicism and many complex structural variants remain partially uncovered. As a conceptual approach, nanopore-based sequencing holds the promise of being a single-molecule-based, long-read and high-resolution technique, with the ability of uncovering the nucleic acid sequence and methylation almost in real time. A key limiting factor of current clinical genetics is the deciphering of key disease-causing genomic sequences. As the technological revolution is expanding regarding genetic data, the interpretation of genotype-phenotype correlations should be made with fine caution, as more and more evidence points toward the presence of more than one pathogenic variant acting together as a result of intergenic interplay in the background of a certain phenotype observed in a patient. This is in conjunction with the observation that many inheritable disorders manifest in a phenotypic spectrum, even in an intra-familial way. In the present review, we summarized the relevant data on nanopore sequencing regarding clinical genomics as well as highlighted the importance and content of pre-test and post-test genetic counselling, yielding a complex approach to phenotype-driven molecular diagnosis. This should significantly lower the time-to-right diagnosis as well lower the time required to complete a currently incomplete genotype-phenotype axis, which will boost the chance of establishing a new actionable diagnosis followed by therapeutical approach.

Overcoming the Pitfalls of Next-Generation Sequencing-Based Molecular Diagnosis of Shwachman-Diamond Syndrome

Shwachman-Diamond syndrome (SDS) is the second most common cause of exocrine pancreatic insufficiency, and 90% of patients carry mutations in the SBDS gene, the most common being the c.183_184delinsCT and c.258+2T>C variants. However, precise detection of these most contributory variants by conventional short-read next-generation sequencing data analysis is limited because of the SBDS/SBDSP1 highly homologous sequences. In this study, an efficient approach was established to infer the haplotype of SBDS based on the expectation-maximization algorithm. The workflow was retrospectively applied to detect the two most common SBDS variants in a Chinese SDS high-risk cohort, and a systematic comparison of variant detection results was performed between the workflow and conventional next-generation sequencing analysis based on Sanger sequencing validation. Among the Chinese SDS high-risk cohort (n = 47) and their available parents (n = 64), the established workflow improved the diagnostic rate for these two variants by 27.7% (95% CI, 15.6%-42.6%) compared with conventional analysis. For overall variant detection, the established workflow achieved 100% (95% CI, 92.5%-100%) concordance with Sanger sequencing, whereas conventional analysis showed only 65.8% accuracy; these results included 25.2% with missed variant calls, 7.2% with diagnosed but inaccurate variant calls, and 1.8% with false-positive calls. With its favorable result in both SDS patient diagnosis and carrier detection performance, the provided workflow showed its potential in clinical application for SDS molecular diagnosis.

Genotype Complements the Phenotype: Identification of the Pathogenicity of an LMNA Splice Variant by Nanopore Long-Read Sequencing in a Large DCM Family

Dilated cardiomyopathy (DCM) is a common cause of heart failure (HF) and is of familial origin in 20−40% of cases. Genetic testing by next-generation sequencing (NGS) has yielded a definite diagnosis in many cases; however, some remain elusive. In this study, we used a combination of NGS, human-induced pluripotent-stem-cell-derived cardiomyocytes (iPSC-CMs) and nanopore long-read sequencing to identify the causal variant in a multi-generational pedigree of DCM. A four-generation family with familial DCM was investigated. Next-generation sequencing (NGS) was performed on 22 family members. Skin biopsies from two affected family members were used to generate iPSCs, which were then differentiated into iPSC-CMs. Short-read RNA sequencing was used for the evaluation of the target gene expression, and long-read RNA nanopore sequencing was used to evaluate the relevance of the splice variants. The pedigree suggested a highly penetrant, autosomal dominant mode of inheritance. The phenotype of the family was suggestive of laminopathy, but previous genetic testing using both Sanger and panel sequencing only yielded conflicting evidence for LMNA p.R644C (rs142000963), which was not fully segregated. By re-sequencing four additional affected family members, further non-coding LMNA variants could be detected: rs149339264, rs199686967, rs201379016, and rs794728589. To explore the roles of these variants, iPSC-CMs were generated. RNA sequencing showed the LMNA expression levels to be significantly lower in the iPSC-CMs of the LMNA variant carriers. We demonstrated a dysregulated sarcomeric structure and altered calcium homeostasis in the iPSC-CMs of the LMNA variant carriers. Using targeted nanopore long-read sequencing, we revealed the biological significance of the variant c.356+1G>A, which generates a novel 5′ splice site in exon 1 of the cardiac isomer of LMNA, causing a nonsense mRNA product with almost complete RNA decay and haploinsufficiency. Using novel molecular analysis and nanopore technology, we demonstrated the pathogenesis of the rs794728589 (c.356+1G>A) splice variant in LMNA. This study highlights the importance of precise diagnostics in the clinical management and workup of cardiomyopathies.

Lack of full sequencing GBA1 studies for patients with Parkinson’s disease in Latin America

Full sequencing of the GBA1 gene in patients with Parkinson’s disease provides a wide screening of pathogenic variants, but less developed regions of the world, like Latin America, may have difficulties in performing full sequencing. We performed a systematic review with meta-analysis to explore the prevalence and the odds ratio of specific GBA1 variants in Parkinson’s disease in Latin America. We noted a lack of full sequencing GBA1 studies in Latin America.

Comprehensive short and long read sequencing analysis for the Gaucher and Parkinson's disease-associated GBA gene

GBA variants carriers are at increased risk of Parkinson’s disease (PD) and Lewy body dementia (LBD). The presence of pseudogene GBAP1 predisposes to structural variants, complicating genetic analysis. We present two methods to resolve recombinant alleles and other variants in GBA: Gauchian, a tool for short-read, whole-genome sequencing data analysis, and Oxford Nanopore sequencing after PCR enrichment. Both methods were concordant for 42 samples carrying a range of recombinants and GBAP1-related mutations, and Gauchian outperformed the GATK Best Practices pipeline. Applying Gauchian to sequencing of over 10,000 individuals shows that copy number variants (CNVs) spanning GBAP1 are relatively common in Africans. CNV frequencies in PD and LBD are similar to controls. Gains may coexist with other mutations in patients, and a modifying effect cannot be excluded. Gauchian detects more GBA variants in LBD than PD, especially severe ones. These findings highlight the importance of accurate GBA analysis in these patients.

Two methods fully resolve the GBA gene: Gauchian, a tool for short-read, whole-genome sequencing data analysis, and Oxford Nanopore sequencing after PCR enrichment. The approach improves our understanding of the relationship between GBA, Gaucher disease and Parkinson disease.

Proof of concept for multiplex amplicon sequencing for mutation identification using the MinION nanopore sequencer

Rapid, cost-effective identification of genetic variants in small candidate genomic regions remains a challenge, particularly for less well equipped or lower throughput laboratories. The application of Oxford Nanopore Technologies’ MinION sequencer has the potential to fulfil this requirement. We demonstrate a proof of concept for a multiplexing assay that pools PCR amplicons for MinION sequencing to enable sequencing of multiple templates from multiple individuals, which could be applied to gene-targeted diagnostics. A combined strategy of barcoding and sample pooling was developed for simultaneous multiplex MinION sequencing of 100 PCR amplicons. The amplicons are family-specific, spanning a total of 30 loci in DNA isolated from 82 human neurodevelopmental cases and family members. The target regions were chosen for further interrogation because a potentially disease-causative variant had been identified in affected individuals following Illumina exome sequencing. The pooled MinION sequences were deconvoluted by aligning to custom references using the minimap2 aligner software. Our multiplexing approach produced an interpretable and expected sequence from 29 of the 30 targeted genetic loci. The sequence variant which was not correctly resolved in the MinION sequence was adjacent to a five nucleotide homopolymer. It is already known that homopolymers present a resolution problem with the MinION approach. Interestingly despite equimolar quantities of PCR amplicon pooled for sequencing, significant variation in the depth of coverage (127×–19,626×; mean = 8321×, std err = 452.99) was observed. We observed independent relationships between depth of coverage and target length, and depth of coverage and GC content. These relationships demonstrate biases of the MinION sequencer for longer templates and those with lower GC content. We demonstrate an efficient approach for variant discovery or confirmation from short DNA templates using the MinION sequencing device. With less than 130 × depth of coverage required for accurate genotyping, the methodology described here allows for rapid highly multiplexed targeted sequencing of large numbers of samples in a minimally equipped laboratory with a potential cost as much 200 × less than that from Sanger sequencing.

Genotyping of familial Mediterranean fever gene (MEFV)-Single nucleotide polymorphism-Comparison of Nanopore with conventional Sanger sequencing

Background

Through continuous innovation and improvement, Nanopore sequencing has become a powerful technology. Because of its fast processing time, low cost, and ability to generate long reads, this sequencing technique would be particularly suitable for clinical diagnostics. However, its raw data accuracy is inferior in contrast to other sequencing technologies. This constraint still results in limited use of Nanopore sequencing in the field of clinical diagnostics and requires further validation and IVD certification.

Methods

We evaluated the performance of latest Nanopore sequencing in combination with a dedicated data-analysis pipeline for single nucleotide polymorphism (SNP) genotyping of the familial Mediterranean fever gene (MEFV) by amplicon sequencing of 47 clinical samples. Mutations in MEFV are associated with Mediterranean fever, a hereditary periodic fever syndrome. Conventional Sanger sequencing, which is commonly applied in clinical genetic diagnostics, was used as a reference method.

Results

Nanopore sequencing enabled the sequencing of 10 target regions within MEFV with high read depth (median read depth 7565x) in all samples and identified a total of 435 SNPs in the whole sample collective, of which 29 were unique. Comparison of both sequencing workflows showed a near perfect agreement with no false negative calls. Precision, Recall, and F1-Score of the Nanopore sequencing workflow were > 0.99, respectively.

Conclusions

These results demonstrated the great potential of current Nanopore sequencing for application in clinical diagnostics, at least for SNP genotyping by amplicon sequencing. Other more complex applications, especially structural variant identification, require further in-depth clinical validation.

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.

Congenital Hyperinsulinism: Current Laboratory-Based Approaches to the Genetic Diagnosis of a Heterogeneous Disease

Congenital hyperinsulinism is characterised by the inappropriate release of insulin during hypoglycaemia. This potentially life-threatening disorder can occur in isolation, or present as a feature of syndromic disease. Establishing the underlying aetiology of the hyperinsulinism is critical for guiding medical management of this condition especially in children with diazoxide-unresponsive hyperinsulinism where the underlying genetics determines whether focal or diffuse pancreatic disease is present. Disease-causing single nucleotide variants affecting over 30 genes are known to cause persistent hyperinsulinism with mutations in the KATP channel genes (ABCC8 and KCNJ11) most commonly identified in children with severe persistent disease. Defects in methylation, changes in chromosome number, and large deletions and duplications disrupting multiple genes are also well described in congenital hyperinsulinism, further highlighting the genetic heterogeneity of this condition. Next-generation sequencing has revolutionised the approach to genetic testing for congenital hyperinsulinism with targeted gene panels, exome, and genome sequencing being highly sensitive methods for the analysis of multiple disease genes in a single reaction. It should though be recognised that limitations remain with next-generation sequencing with no single application able to detect all reported forms of genetic variation. This is an important consideration for hyperinsulinism genetic testing as comprehensive screening may require multiple investigations.

Hidden biases in germline structural variant detection

BACKGROUND

Genomic structural variations (SV) are important determinants of genotypic and phenotypic changes in many organisms. However, the detection of SV from next-generation sequencing data remains challenging.

RESULTS

In this study, DNA from a Chinese family quartet is sequenced at three different sequencing centers in triplicate. A total of 288 derivative data sets are generated utilizing different analysis pipelines and compared to identify sources of analytical variability. Mapping methods provide the major contribution to variability, followed by sequencing centers and replicates. Interestingly, SV supported by only one center or replicate often represent true positives with 47.02% and 45.44% overlapping the long-read SV call set, respectively. This is consistent with an overall higher false negative rate for SV calling in centers and replicates compared to mappers (15.72%). Finally, we observe that the SV calling variability also persists in a genotyping approach, indicating the impact of the underlying sequencing and preparation approaches.

CONCLUSIONS

This study provides the first detailed insights into the sources of variability in SV identification from next-generation sequencing and highlights remaining challenges in SV calling for large cohorts. We further give recommendations on how to reduce SV calling variability and the choice of alignment methodology.

Using DNA testing for the precise, definite, and low-cost diagnosis of sickle cell disease and other Haemoglobinopathies: findings from Tanzania

Background

Sickle cell disease (SCD) is an important cause of under-five mortality. Tanzania is the 5th country in the world with the highest births prevalence of SCD individuals. Significant advances in the neonatal diagnosis of SCD using rapid point-of-care testing have been made. However genetic confirmation is still required for positive cases, in uncertain cases, in multiply transfused patients, to resolve compound heterozygosity (Hb S/ β⁰ Thal or Hb S/ β⁺ thal) not uncommon in the coastal regions of East Africa and increasingly also for pre-marital counselling and potentially for future curative approaches such as gene therapy. The currently available DNA tests are prohibitively expensive. Here, we describe an easy-to-use, affordable and accurate β-globin sequencing approach that can be easily integrated within existing NBS for SCD and other haemoglobinopathies especially in Low- and Middle-income Countries.

Aim

To evaluate an affordable DNA technology for the diagnosis of Sickle cell disease and other haemoglobinopathies in a resource-limited setting.

Methods

Laboratory-based validation study was conducted by Muhimbili University of Health and Allied Sciences and the University of Oxford involving sequencing of the entire β -haemoglobin locus using the Oxford Nanopore MinION platform. A total number of 36 Dried blood spots and whole blood samples were subjected to conventional protein-based methods (isoelectric focusing, HPLC), and/or sequenced by the Sanger method as comparators.

Results

Sequencing results for SCD using the MinION were 100% concordant with those from the Sanger method. In addition, the long-read DNA sequencing method enabled the resolution of cases with unusual phenotypes which make up 1% of all children in Tanzania. The cost is £11/ sample for consumables, which is cheaper compared to other sequencing platforms.

Conclusions

This is the first report of a comprehensive single DNA assay as a definitive diagnostic test for SCD and other haemoglobinopathies. The test is fast, precise, accurate and affordable.

Long-read single molecule real-time (SMRT) sequencing of GBA1 locus in Gaucher disease national cohort from Argentina reveals high frequency of complex allele underlying severe skeletal phenotypes: Collaborative study from the Argentine Group for Diagnosis and Treatment of Gaucher Disease

Gaucher disease is reckoned for extreme phenotypic diversity that does not show consistent genotype/phenotype correlations. In Argentina, a national collaborative group, Grupo Argentino de Diagnóstico y Tratamiento de la Enfermedad de Gaucher, GADTEG, have delineated uniformly severe type 1 Gaucher disease manifestations presenting in childhood with large burden of irreversible skeletal disease. Here using Long-Read Single Molecule Real-Time (SMRT) Sequencing of GBA1 locus, we show that RecNciI allele is highly prevalent and associates with severe skeletal manifestations in childhood.

Detection of single nucleotide and copy number variants in the Fabry disease-associated GLA gene using nanopore sequencing

More than 900 variants have been described in the GLA gene. Some intronic variants and copy number variants in GLA can cause Fabry disease but will not be detected by classical Sanger sequence. We aimed to design and validate a method for sequencing the GLA gene using long-read Oxford Nanopore sequencing technology. Twelve Fabry patients were blindly analyzed, both by conventional Sanger sequence and by long-read sequencing of a 13 kb PCR amplicon. We used minimap2 to align the long-read data and Nanopolish and Sniffles to call variants. All the variants detected by Sanger (including a deep intronic variant) were also detected by long-read sequencing. One patient had a deletion that was not detected by Sanger sequencing but was detected by the new technology. Our long-read sequencing-based method was able to detect missense variants and an exonic deletion, with the added advantage of intronic analysis. It can be used as an efficient and cost-effective tool for screening and diagnosing Fabry disease.

Nanopore sequencing technology, bioinformatics and applications

Rapid advances in nanopore technologies for sequencing single long DNA and RNA molecules have led to substantial improvements in accuracy, read length and throughput. These breakthroughs have required extensive development of experimental and bioinformatics methods to fully exploit nanopore long reads for investigations of genomes, transcriptomes, epigenomes and epitranscriptomes. Nanopore sequencing is being applied in genome assembly, full-length transcript detection and base modification detection and in more specialized areas, such as rapid clinical diagnoses and outbreak surveillance. Many opportunities remain for improving data quality and analytical approaches through the development of new nanopores, base-calling methods and experimental protocols tailored to particular applications.

Glucocerebrosidase mutations: A paradigm for neurodegeneration pathways

Biallelic (homozygous or compound heterozygous) glucocerebrosidase gene (GBA) mutations cause Gaucher disease, whereas heterozygous mutations are numerically the most important genetic risk factor for Parkinson disease (PD) and are associated with the development of other synucleinopathies, notably Dementia with Lewy Bodies. This phenomenon is not limited to GBA, with converging evidence highlighting further examples of autosomal recessive disease genes increasing neurodegeneration risk in heterozygous mutation carriers. Nevertheless, despite extensive research, the cellular mechanisms by which mutations in GBA, encoding lysosomal enzyme β-glucocerebrosidase (GCase), predispose to neurodegeneration remain incompletely understood. Alpha-synuclein (A-SYN) accumulation, autophagic lysosomal dysfunction, mitochondrial abnormalities, ER stress and neuroinflammation have been proposed as candidate pathogenic pathways in GBA-linked PD. The observation of GCase and A-SYN interactions in PD initiated the development and evaluation of GCase-targeted therapeutics in PD clinical trials.

Do Not Miss the (Genetic) Diagnosis of Gaucher Syndrome: A Narrative Review on Diagnostic Clues and Management in Severe Prenatal and Perinatal-Lethal Sporadic Cases

With a growing number of proved therapies and clinical trials for many lysosomal storage disorders (LSDs), a lot of hope for many patients and families exists. However, there are sometimes cases with poor prognosis, fatal outcomes when our efforts must be directed towards a prompt and correct genetic diagnosis, which offers the only possibility of providing the family with appropriate prevention and treatment. To address this issue, in this article, we present the clinical and genetic hallmarks of the lethal form of Gaucher disease (PLGD) and discuss the potential management. We hope that this will draw attention to its specific manifestations (such as collodion-baby phenotype, ichthyosis, arthrogryposis), which differ from best-known GD complications and ensure appropriate diagnostic assessment to provide families at risk with reliable counselling and treatment to avoid the medical complication of GD.

The remote assessment of parkinsonism supporting ongoing development of interventions in Gaucher disease

Mutations in GBA which are causative of Gaucher disease in their biallelic form, are the most common genetic risk factor for Parkinson's disease (PD). The diagnosis of PD relies upon clinically defined motor features which appear after irreversible neurodegeneration. Prodromal symptoms of PD may provide a means to predict latent pathology, years before the onset of motor features. Previous work has reported prodromal features of PD in GBA mutation carriers, however this has been insufficiently sensitive to identify those that will develop PD. The Remote Assessment of Parkinsonism Supporting Ongoing Development of Interventions in Gaucher Disease (RAPSODI GD) study assesses a large cohort of GBA mutation carriers, to aid development of procedures for earlier diagnosis of PD.

Towards population-scale long-read sequencing

Long-read sequencing technologies have now reached a level of accuracy and yield that allows their application to variant detection at a scale of tens to thousands of samples. Concomitant with the development of new computational tools, the first population-scale studies involving long-read sequencing have emerged over the past 2 years and, given the continuous advancement of the field, many more are likely to follow. In this Review, we survey recent developments in population-scale long-read sequencing, highlight potential challenges of a scaled-up approach and provide guidance regarding experimental design. We provide an overview of current long-read sequencing platforms, variant calling methodologies and approaches for de novo assemblies and reference-based mapping approaches. Furthermore, we summarize strategies for variant validation, genotyping and predicting functional impact and emphasize challenges remaining in achieving long-read sequencing at a population scale.

The Genetics of Parkinson's Disease and Implications for Clinical Practice

The genetic landscape of Parkinson’s disease (PD) is characterised by rare high penetrance pathogenic variants causing familial disease, genetic risk factor variants driving PD risk in a significant minority in PD cases and high frequency, low penetrance variants, which contribute a small increase of the risk of developing sporadic PD. This knowledge has the potential to have a major impact in the clinical care of people with PD. We summarise these genetic influences and discuss the implications for therapeutics and clinical trial design.

Intronic Haplotypes in the GBA Gene Do Not Predict Age at Diagnosis of Parkinson's Disease

BACKGROUND

GBA mutations are a common risk factor for Parkinson's disease (PD). A recent study has suggested that GBA haplotypes, identified by intronic variants, can affect age at diagnosis of PD.

OBJECTIVES

In this study, we assess this hypothesis using long reads across a large cohort and the publicly available Accelerating Medicines Partnership-Parkinson's Disease (AMP-PD) cohort.

METHODS

We recruited a PD cohort through the Remote Assessment of Parkinsonism Supporting Ongoing Development of Interventions in Gaucher Disease study (RAPSODI) and sequenced GBA using Oxford Nanopore technology. Genetic and clinical data on the full AMP-PD cohort were obtained from the online portal of the consortium.

RESULTS

A total of 1417 participants were analyzed. There was no significant difference in age at PD diagnosis between the two main haplotypes of the GBA gene.

CONCLUSIONS

GBA haplotypes do not affect age at diagnosis of PD in the two independent cohorts studied. © 2021 The Authors. Movement Disorders published by Wiley Periodicals LLC on behalf of International Parkinson and Movement Disorder Society.

Nanopore Sequencing and Its Clinical Applications

Nanopore sequencing is a method for determining the order and modifications of DNA/RNA nucleotides by detecting the electric current variations when DNA/RNA oligonucleotides pass through the nanometer-sized hole (nanopore). Nanopore-based DNA analysis techniques have been commercialized by Oxford Nanopore Technologies, NabSys, and Sequenom, and widely used in scientific researches recently including human genomics, cancer, metagenomics, plant sciences, etc., moreover, it also has potential applications in the field of healthcare due to its fast turn-around time, portable and real-time data analysis. Those features make it a promising technology for the point-of-care testing (POCT) and its potential clinical applications are briefly discussed in this chapter.

A deep intronic variant is a common cause of OTC deficiency in individuals with previously negative genetic testing

Pathogenic variants in non-coding regions of genes encoding enzymes or transporters of the urea cycle can lead to urea cycle disorders (UCDs). However, not all commercially available testing platforms interrogate these regions. Here, we used a gene panel based on massively parallel sequencing (MPS) in 10 individuals with clinical or pedigree-based evidence of a proximal UCD but without a molecular confirmation of the diagnosis. We identified causal variant(s) in 5 of 10 individuals, including in 3 of 7 individuals in whom prior molecular testing was unrevealing. We show that a deep-intronic pathogenic variant in OTC, c.540+265G>A, is an important cause of ornithine transcarbamylase (OTC) deficiency.

PERHAPS: Paired-End short Reads-based HAPlotyping from next-generation Sequencing data

The identification of rare haplotypes may greatly expand our knowledge in the genetic architecture of both complex and monogenic traits. To this aim, we developed PERHAPS (Paired-End short Reads-based HAPlotyping from next-generation Sequencing data), a new and simple approach to directly call haplotypes from short-read, paired-end Next Generation Sequencing (NGS) data. To benchmark this method, we considered the APOE classic polymorphism (*1/*2/*3/*4), since it represents one of the best examples of functional polymorphism arising from the haplotype combination of two Single Nucleotide Polymorphisms (SNPs). We leveraged the big Whole Exome Sequencing (WES) and SNP-array data obtained from the multi-ethnic UK BioBank (UKBB, N=48,855). By applying PERHAPS, based on piecing together the paired-end reads according to their FASTQ-labels, we extracted the haplotype data, along with their frequencies and the individual diplotype. Concordance rates between WES directly called diplotypes and the ones generated through statistical pre-phasing and imputation of SNP-array data are extremely high (>99%), either when stratifying the sample by SNP-array genotyping batch or self-reported ethnic group. Hardy-Weinberg Equilibrium tests and the comparison of obtained haplotype frequencies with the ones available from the 1000 Genome Project further supported the reliability of PERHAPS. Notably, we were able to determine the existence of the rare APOE*1 haplotype in two unrelated African subjects from UKBB, supporting its presence at appreciable frequency (approximatively 0.5%) in the African Yoruba population. Despite acknowledging some technical shortcomings, PERHAPS represents a novel and simple approach that will partly overcome the limitations in direct haplotype calling from short read-based sequencing.

A Long-Read Sequencing Approach for Direct Haplotype Phasing in Clinical Settings

The reconstruction of individual haplotypes can facilitate the interpretation of disease risks; however, high costs and technical challenges still hinder their assessment in clinical settings. Second-generation sequencing is the gold standard for variant discovery but, due to the production of short reads covering small genomic regions, allows only indirect haplotyping based on statistical methods. In contrast, third-generation methods such as the nanopore sequencing platform developed by Oxford Nanopore Technologies (ONT) generate long reads that can be used for direct haplotyping, with fewer drawbacks. However, robust standards for variant phasing in ONT-based target resequencing efforts are not yet available. In this study, we presented a streamlined proof-of-concept workflow for variant calling and phasing based on ONT data in a clinically relevant 12-kb region of the APOE locus, a hotspot for variants and haplotypes associated with aging-related diseases and longevity. Starting with sequencing data from simple amplicons of the target locus, we demonstrated that ONT data allow for reliable single-nucleotide variant (SNV) calling and phasing from as little as 60 reads, although the recognition of indels is less efficient. Even so, we identified the best combination of ONT read sets (600) and software (BWA/Minimap2 and HapCUT2) that enables full haplotype reconstruction when both SNVs and indels have been identified previously using a highly-accurate sequencing platform. In conclusion, we established a rapid and inexpensive workflow for variant phasing based on ONT long reads. This allowed for the analysis of multiple samples in parallel and can easily be implemented in routine clinical practice, including diagnostic testing.

Electronic DNA Analysis of CSF Cell-free Tumor DNA to Quantify Multi-gene Molecular Response in Pediatric High-grade Glioma

PURPOSE

Pediatric high-grade glioma (pHGG) diagnosis portends poor prognosis and therapeutic monitoring remains difficult. Tumors release cell-free tumor DNA (cf-tDNA) into cerebrospinal fluid (CSF), allowing for potential detection of tumor-associated mutations by CSF sampling. We hypothesized that direct, electronic analysis of cf-tDNA with a handheld platform (Oxford Nanopore MinION) could quantify patient-specific CSF cf-tDNA variant allele fraction (VAF) with improved speed and limit of detection compared with established methods.

EXPERIMENTAL DESIGN

We performed ultra-short fragment (100-200 bp) PCR amplification of cf-tDNA for clinically actionable alterations in CSF and tumor samples from patients with pHGG (n = 12) alongside nontumor CSF (n = 6). PCR products underwent rapid amplicon-based sequencing by Oxford Nanopore Technology (Nanopore) with quantification of VAF. Additional comparison to next-generation sequencing (NGS) and droplet digital PCR (ddPCR) was performed.

RESULTS

Nanopore demonstrated 85% sensitivity and 100% specificity in CSF samples (n = 127 replicates) with 0.1 femtomole DNA limit of detection and 12-hour results, all of which compared favorably with NGS. Multiplexed analysis provided concurrent analysis of H3.3A (H3F3A) and H3C2 (HIST1H3B) mutations in a nonbiopsied patient and results were confirmed by ddPCR. Serial CSF cf-tDNA sequencing by Nanopore demonstrated correlation of radiological response on a clinical trial, with one patient showing dramatic multi-gene molecular response that predicted long-term clinical response.

CONCLUSIONS

Nanopore sequencing of ultra-short pHGG CSF cf-tDNA fragments is feasible, efficient, and sensitive with low-input samples thus overcoming many of the barriers restricting wider use of CSF cf-tDNA diagnosis and monitoring in this patient population.

Genetic modifiers of risk and age at onset in GBA associated Parkinson's disease and Lewy body dementia

Parkinson's disease is a genetically complex disorder. Multiple genes have been shown to contribute to the risk of Parkinson's disease, and currently 90 independent risk variants have been identified by genome-wide association studies. Thus far, a number of genes (including SNCA, LRRK2, and GBA) have been shown to contain variability across a spectrum of frequency and effect, from rare, highly penetrant variants to common risk alleles with small effect sizes. Variants in GBA, encoding the enzyme glucocerebrosidase, are associated with Lewy body diseases such as Parkinson's disease and Lewy body dementia. These variants, which reduce or abolish enzymatic activity, confer a spectrum of disease risk, from 1.4- to >10-fold. An outstanding question in the field is what other genetic factors that influence GBA-associated risk for disease, and whether these overlap with known Parkinson's disease risk variants. Using multiple, large case-control datasets, totalling 217 165 individuals (22 757 Parkinson's disease cases, 13 431 Parkinson's disease proxy cases, 622 Lewy body dementia cases and 180 355 controls), we identified 1691 Parkinson's disease cases, 81 Lewy body dementia cases, 711 proxy cases and 7624 controls with a GBA variant (p.E326K, p.T369M or p.N370S). We performed a genome-wide association study and analysed the most recent Parkinson's disease-associated genetic risk score to detect genetic influences on GBA risk and age at onset. We attempted to replicate our findings in two independent datasets, including the personal genetics company 23andMe, Inc. and whole-genome sequencing data. Our analysis showed that the overall Parkinson's disease genetic risk score modifies risk for disease and decreases age at onset in carriers of GBA variants. Notably, this effect was consistent across all tested GBA risk variants. Dissecting this signal demonstrated that variants in close proximity to SNCA and CTSB (encoding cathepsin B) are the most significant contributors. Risk variants in the CTSB locus were identified to decrease mRNA expression of CTSB. Additional analyses suggest a possible genetic interaction between GBA and CTSB and GBA p.N370S induced pluripotent cell-derived neurons were shown to have decreased cathepsin B expression compared to controls. These data provide a genetic basis for modification of GBA-associated Parkinson's disease risk and age at onset, although the total contribution of common genetics variants is not large. We further demonstrate that common variability at genes implicated in lysosomal function exerts the largest effect on GBA associated risk for disease. Further, these results have implications for selection of GBA carriers for therapeutic interventions.

PhaseME: Automatic rapid assessment of phasing quality and phasing improvement

BACKGROUND

The detection of which mutations are occurring on the same DNA molecule is essential to predict their consequences. This can be achieved by phasing the genomic variations. Nevertheless, state-of-the-art haplotype phasing is currently a black box in which the accuracy and quality of the reconstructed haplotypes are hard to assess.

FINDINGS

Here we present PhaseME, a versatile method to provide insights into and improvement of sample phasing results based on linkage data. We showcase the performance and the importance of PhaseME by comparing phasing information obtained from Pacific Biosciences including both continuous long reads and high-quality consensus reads, Oxford Nanopore Technologies, 10x Genomics, and Illumina sequencing technologies. We found that 10x Genomics and Oxford Nanopore phasing can be significantly improved while retaining a high N50 and completeness of phase blocks. PhaseME generates reports and summary plots to provide insights into phasing performance and correctness. We observed unique phasing issues for each of the sequencing technologies, highlighting the necessity of quality assessments. PhaseME is able to decrease the Hamming error rate significantly by 22.4% on average across all 5 technologies. Additionally, a significant improvement is obtained in the reduction of long switch errors. Especially for high-quality consensus reads, the improvement is 54.6% in return for only a 5% decrease in phase block N50 length.

CONCLUSIONS

PhaseME is a universal method to assess the phasing quality and accuracy and improves the quality of phasing using linkage information. The package is freely available at https://github.com/smajidian/phaseme.

Nanopore sequencing of the pharmacogene CYP2D6 allows simultaneous haplotyping and detection of duplications

Aim: Long read sequencing offers the promise of overcoming some of the challenges in accurate genotyping of complex genes, along with the advantage of straightforward variant phasing. We have established methods for sequencing and haplotyping of the whole CYP2D6 gene using nanopore sequencing. Materials and methods: 32 samples covering various haplotypes including gene duplication were sequenced on the GridION platform. Results: Haplotypes of 52 alleles matched accurately to known star (*) allele subvariants, with the remaining 12 being assigned as new alleles, or new subvariants of known alleles. Duplicated alleles could be detected by analyzing the allelic balance. Conclusion: Nanopore sequencing of CYP2D6 offers a high throughput method for accurate haplotyping, detection of new variants and determination of duplicated alleles.

Targeted nanopore sequencing with Cas9-guided adapter ligation

Despite recent improvements in sequencing methods, there remains a need for assays that provide high sequencing depth and comprehensive variant detection. Current methods1-4 are limited by the loss of native modifications, short read length, high input requirements, low yield or long protocols. In the present study, we describe nanopore Cas9-targeted sequencing (nCATS), an enrichment strategy that uses targeted cleavage of chromosomal DNA with Cas9 to ligate adapters for nanopore sequencing. We show that nCATS can simultaneously assess haplotype-resolved single-nucleotide variants, structural variations and CpG methylation. We apply nCATS to four cell lines, to a cell-line-derived xenograft, and to normal and paired tumor/normal primary human breast tissue. Median sequencing coverage was 675× using a MinION flow cell and 34× using the smaller Flongle flow cell. The nCATS sequencing requires only ~3 μg of genomic DNA and can target a large number of loci in a single reaction. The method will facilitate the use of long-read sequencing in research and in the clinic.

Nanopore Sequencing in Blood Diseases: A Wide Range of Opportunities

The molecular pathogenesis of hematological diseases is often driven by genetic and epigenetic alterations. Next-generation sequencing has considerably increased our genomic knowledge of these disorders becoming ever more widespread in clinical practice. In 2012 Oxford Nanopore Technologies (ONT) released the MinION, the first long-read nanopore-based sequencer, overcoming the main limits of short-reads sequences generation. In the last years, several nanopore sequencing approaches have been performed in various "-omic" sciences; this review focuses on the challenge to introduce ONT devices in the hematological field, showing advantages, disadvantages and future perspectives of this technology in the precision medicine era.

A novel mutation deep within intron 7 of the GBA gene causes Gaucher disease

BACKGROUND

Mutations in the GBA gene that encodes the lysosomal enzyme acid β-glucocerebrosidase cause Gaucher disease (GD), the most common lysosomal storage disorder. Most of the mutations are missense/nonsense, however, a few splicing mutations within or close to conserved consensus donor or acceptor splice sites have also been described. The aim of the study was to identify the mutation(s) in a Cypriot patient with type I GD.

METHODS

The genomic DNA of the proband was screened for nine common mutations using Polymerase chain reaction-restriction fragment length polymorphism (PCR-RFLP) analysis. All exons and exon-intron boundaries, and the 5'UTR and 3'UTR regions of the GBA gene, were investigated by Sanger sequencing. RNA analysis was performed using standard procedures, and the abnormal transcript was further cloned into pGEM-T-Easy plasmid vector and sequenced. The relevant intronic region was further sequenced by the Sanger method to identify the genetic variant.

RESULTS

A novel point mutation, g.12599C > A (c.999 + 242C > A), was detected deep in intron 7 of the GBA gene. This type of mutation has been previously described for other diseases but this is the first time, as far as we know, that it is described for GD. This mutation creates a new donor splice site leading to aberrant splicing and resulting in the insertion of the first 239nt of intron 7 as a pseudoexon in the mRNA, creating a premature stop codon.

CONCLUSION

This study expands the mutation spectrum of GD and highlights the importance of RNA sequencing for the molecular diagnosis of patients bearing mutations in nonexonic regions.

A new era of long-read sequencing for cancer genomics

Cancer is a disease largely caused by genomic aberrations. Utilizing many rapidly emerging sequencing technologies, researchers have studied cancer genomes to understand the molecular statuses of cancer cells and to reveal their vulnerabilities, such as driver mutations or gene expression. Long-read technologies enable us to identify and characterize novel types of cancerous mutations, including complicated structural variants in haplotype resolution. In this review, we introduce three representative platforms for long-read sequencing and research trends of cancer genomics with long-read data. Further, we describe that aberrant transcriptome and epigenome statuses, namely, fusion transcripts, as well as aberrant transcript isoforms and the phase information of DNA methylation, are able to be elucidated by long-read sequencers. Long-read sequencing may shed light on novel types of aberrations in cancer genomics that are being missed by conventional short-read sequencing analyses.

Nanopore sequencing of the glucocerebrosidase (GBA) gene in a New Zealand Parkinson's disease cohort

INTRODUCTION

Bi-allelic mutations in the gene for glucocerebrosidase (GBA) cause Gaucher disease, an autosomal recessive lysosomal storage disorder. Gaucher disease causing GBA mutations in the heterozygous state are also high risk factors for Parkinson's disease (PD). GBA analysis is challenging due to a related pseudogene and structural variations (SVs) that can occur at this locus. We have applied and refined a recently developed nanopore DNA sequencing method to analyze GBA variants in a clinically assessed New Zealand longitudinal cohort of PD.

METHOD

We examined amplicons encompassing the coding region of GBA (8.9 kb) from 229 PD cases and 50 healthy controls using the GridION nanopore sequencing platform, and Sanger validation.

RESULTS

We detected 23 variants in 21 PD cases (9.2% of patients). We detected modest PD risk variant p.N409S (rs76763715) in one case, p.E365K (rs2230288) in 12 cases, and p.T408 M (rs75548401) in seven cases, one of whom also had p.E365K. We additionally detected the possible risk variants p.R78C (rs146774384) in one case, p.D179H (rs147138516) in one case which occurred on the same haplotype as p.E365K, and one novel variant c.335C > T or p.(L335 = ), that potentially impacts splicing of GBA transcripts. Additionally, we found a higher prevalence of dementia among patients with GBA variants.

CONCLUSION

This work confirmed the utility of nanopore sequencing as a high-throughput method to identify known and novel GBA variants, and to assign precise haplotypes. Our observations may contribute to improved understanding of the effects of variants on disease pathogenesis, and to the development of more targeted treatments.

Structural variant calling: the long and the short of it

Recent research into structural variants (SVs) has established their importance to medicine and molecular biology, elucidating their role in various diseases, regulation of gene expression, ethnic diversity, and large-scale chromosome evolution-giving rise to the differences within populations and among species. Nevertheless, characterizing SVs and determining the optimal approach for a given experimental design remains a computational and scientific challenge. Multiple approaches have emerged to target various SV classes, zygosities, and size ranges. Here, we review these approaches with respect to their ability to infer SVs across the full spectrum of large, complex variations and present computational methods for each approach.

Genome sequencing and implications for rare disorders

The practice of genomic medicine stands to revolutionize our approach to medical care, and to realize this goal will require discovery of the relationship between rare variation at each of the ~ 20,000 protein-coding genes and their consequent impact on individual health and expression of Mendelian disease. The step-wise evolution of broad-based, genome-wide cytogenetic and molecular genomic testing approaches (karyotyping, chromosomal microarray [CMA], exome sequencing [ES]) has driven much of the rare disease discovery to this point, with genome sequencing representing the newest member of this team. Each step has brought increased sensitivity to interrogate individual genomic variation in an unbiased method that does not require clinical prediction of the locus or loci involved. Notably, each step has also brought unique limitations in variant detection, for example, the low sensitivity of ES for detection of triploidy, and of CMA for detection of copy neutral structural variants. The utility of genome sequencing (GS) as a clinical molecular diagnostic test, and the increased sensitivity afforded by addition of long-read sequencing or other -omics technologies such as RNAseq or metabolomics, are not yet fully explored, though recent work supports improved sensitivity of variant detection, at least in a subset of cases. The utility of GS will also rely upon further elucidation of the complexities of genetic and allelic heterogeneity, multilocus rare variation, and the impact of rare and common variation at a locus, as well as advances in functional annotation of identified variants. Much discovery remains to be done before the potential utility of GS is fully appreciated.