Detection of Single Nucleotide and Copy Number Variants in the Fabry Disease-associated GLA Gene Using Nanopore Sequencing

doi:10.21203/rs.3.rs-579267/v1

Download PDF

Research Article

Detection of Single Nucleotide and Copy Number Variants in the Fabry Disease-associated GLA Gene Using Nanopore Sequencing

https://doi.org/10.21203/rs.3.rs-579267/v1

This work is licensed under a CC BY 4.0 License

Journal Publication

published 16 Nov, 2021

Read the published version in Scientific Reports →

You are reading this latest preprint version

Introduction: More than one thousand 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.

Aims: We aimed to design and validate a method for sequencing the GLA gene using long read Oxford Nanopore sequencing technology.

Methods: Twelve Fabry patients were blindly analyzed, both by conventional Sanger sequence and by long read sequencing of a 13kb PCR amplicon. We used minimap2 to align the long read data and Nanopolish and Sniffles to call variants.

Results: 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.

Conclusions: 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.

Molecular Biology

Molecular Genetics

GLA

Oxford Nanopore sequencing technology

Fabry disease

PCR amplicon

Fabry disease (FD; OMIM: 301500) is an X-linked lysosomal storage disorder caused by deficiency of the alpha galactosidase A enzyme, resulting in accumulation of glycosphingolipids, particularly globotriaosylceramide (GL-3, Gb3, CTH) and globotriaosylsphingosine (Lyso-GL-3, lyso-Gb3) ¹.

These lipids progressively accumulate in all cell types and organs, resulting in the development of a multisystem disorder.

FD is monogenic and caused by loss of function variants in the GLA gene. This gene is located on the long arm of the X chromosome at the Xq22 position. Most cases are hereditary, and cases of spontaneous pathogenic variants are rare ^2–4. Over 900 different pathogenic variants have been described as the cause of the disease ⁵.

The GLA gene is approximately 12 kb long and spans seven exons. FD can be caused by several types of molecular variants in this gene: missense (57%), nonsense (11%), partial deletions (6%), insertions (6%), and defects in the processing of RNA, which lead to aberrant splicing (6%) ⁶. The correlation between genotype and phenotype is complex, since the same variant may determine different clinical manifestations ⁷.

Variants determine the clinical phenotype - classic or late-onset ^8–10. E.g., frame-shift variants such as deletions and duplications, non-sense and some missense and splicing variants lead to a zero or very low α-GAL activity and, consequently, the classic phenotype. Such patients are at high risk of developing a small-fiber neuropathy, progressive proteinuric kidney disease, chronic diarrhea and abdominal pains, fibrotic cardiac disease resulting in rhythm and conduction disturbances, progressive hypertrophic cardiomyopathy, and cerebrovascular stroke ¹¹. Other missense and splicing variants lead to a significant residual α-GAL activity and, consequently, the late-onset phenotype. Such patients suffer from cardiac or, more rarely, renal disease.

Sometimes, although the clinical suspicion for Fabry disease is high, the genetic testing is difficult due to several pitfalls: mosaicism, deep intronic variants, copy number variants (CNVs), etc… Typically, the genetic diagnosis of FD is performed by Sanger sequencing of all exon and exon boundaries in the GLA gene. However, this method does not detect CNVs, deep intronic variants and mosaics.

Previously, Bae et al reported about a late-onset male Fabry disease patient with somatic mosaicism of a classical GLA pathogenic variant ¹². Next generation sequencing enabled diagnosis since the male patient was mosaic for an SNV (58% of the reads detected the mutated variant although 100% variant frequency would be expected for an X-linked allele). This was missed by Sanger sequencing.

Nanopore-based sequencing has recently emerged as a powerful technology for nucleic acid sequencing in all fields of biology ^13,14. This technology was successfully commercialized by Oxford Nanopore Technologies (ONT) due to its ability to sequence ultra-long reads, perform real-time basecalling and analysis, provide base modification detection, all with short sample preparation time and low instrument costs ¹⁵. These advantages allowed the successful complete assembly of human chromosomes ¹⁶, tumor structural and epigenetic variations associated with cancer ¹⁷, detection of imprinted DNA methylation ¹⁸, identification of novel transcript isoforms ¹⁹ and fast detection of viral and bacterial pathogens ^20–23, including the SARS-CoV-2 virus ^24,25.

Amplification-based targeted sequencing can be performed using both short and long read sequencing as well as Sanger sequencing. However, as ONT sequencing is able to sequence full amplicons in a single read, it has a clear advantage in detection of structural and intronic variants, variant phasing and accurate mapping in cases of highly homologues genes (e.g. pseudogene, duplicated gene). For the lysosomal storage disorder Gaucher disease, ONT sequencing was shown to accurately analyze the GBA gene, previously considered challenging due to a nearby pseudogene ²⁶. This analysis also identified intronic variants that are missed by the classical Sanger sequencing-based GLA genotyping. Similar approaches were used to identify variations in FCGR3A gene ²⁷ and TP53 ²⁸.

The aim of our study was to determine the accuracy of ONT long amplicon sequencing for detection of variants in the GLA gene.

The clinical signs and symptoms are presented in Table 1. All patients were previously diagnosed by measuring the alpha galactosidase A enzyme activity and Sanger sequencing. The samples tested were blinded except for gender, meaning that aside from male sex, the genetic variant and any relationship among the samples was not shared with our lab.

Table 1

Clinical presentation of the cohort of Fabry patients. Acrop - acroparestesias, Hypo - hypohidriosis, Angio - angioceratoma, Cornea - cornea verticillata, Nephro – nephropathy, Cardio – cardiomyopathy, Stroke – history of stroke.
Sample	Therapy	Acrop	Hypo	Angio	Cornea	Nephro	Cardio	Stroke	phenotype	Note
1	Agalsidase α	+							classic
2	Agalsidase β	+	+	+	+				classic
3	Agalsidase β	+	+		+	+			classic
4	Agalsidase β	+	+	+	+				classic
5	Agalsidase β	+	+	+	+	+	+		classic
6	No therapy								benign
7	Agalsidase α	+	+	+	+	+			classic	brother of 8
8	Agalsidase α	+	+	+	+	+	+	+	classic	brother of 7
9	Agalsidase α	+	+		+	+	+		classic
10	Agalsidase β						+		later-onset
11	Agalsidase α	+	+	+	+	+	+		classic
12	Agalsidase α	+	+	+	+	+	+		classic

Nanopore testing of the patients:

GLA amplicon sequencing

In order to detect genomic variants in the GLA locus, we designed a PCR amplicon that produces a 13kb product including the entire gene, and 800 bp and 2000 bp up- and down-stream sequences, respectively. Sequencing the pooled PCR products of the 12 samples on one MinION flow cell yielded a median of 88,800 reads and 617 million bases per sample (Table S1). Read length distributions of the 12 samples show a peak around 13kb, demonstrating that most reads are of full amplicons (Fig. S1). Mapping the reads to the human reference genome showed a median of 44,000X coverage per sample around the GLA region.

SNV and indel detection

The nanopolish tool was then used to call single nucleotide variants (SNV) and short indels. The resulting variants were filtered based on quality score to eliminate false positives calls. The filtered variants were first searched in the ClinVar database for any known classification. The variants that did not match any entry in known databases were manually classified for their predicted effect on the protein using several prediction tools (see Methods). A summary of the variants with likely pathogenic effect on the protein translation is shown in Table 2. In 10 of the 12 samples analyzed we identified a pathogenic or likely pathogenic exonic SNV (6 samples) or short indel (4 samples). For another sample, we detected an intronic SNV (IVS53 + 405T > G) at a possible branch site that is likely to affect splicing. This variant was detected previously by Sanger sequencing of multiple GLA intronic amplicons. With nanopore sequencing, the same variant was detected by a single PCR assay.

Table 2

Genotyping results based on ONT amplicon sequencing, and predicted consequence of the genetic variant on *GLA*. * - two samples from the same family were tested; N/A = not applicable
ONT Barcode	GLA nucleotide variant (Accession: NM_000169)	GLA protein variant (Accession: NP_000160)	Comments
1	IVS53 + 405T > G	N/A-deep intronic	deep intronic variant; possible branch site; may affect splicing
2	c.1147_1149delTTC	p.Phe383del	pathogenic
3	c.744_745delTA*	p.Phe248LeufsX7	pathogenic; same haplotype as barcode05 suggests blood relationship
4	c.559_560delAT	p.Met187Valfs*6	pathogenic
5	c.744_745delTA*	p.Phe248LeufsX7	pathogenic; same haplotype as barcode03 suggests blood relationship
6	c.352C > T	p.Arg118Cys	Conflicting_interpretations_of_pathogenicity
7	c.370-2A > G*	N/A- splicing	pathogenic; same haplotype as barcode08 suggests blood relationship
8	c.370-2A > G*	N/A- splicing	pathogenic; same haplotype as barcode07 suggests blood relationship
9	c.704C > A	p.Ser235Tyr	pathogenic
10	c.337T > C	p. Phe113Leu	pathogenic
11	Exon 2 deletion)	N/A	pathogenic; 2914bp deletion removes GLA exon2 (chrX:100658307–100661221
12	c.581C > T	p.Thr194Ile	likely pathogenic

Structural variants

Although structural variants (SVs), such as insertions and deletions longer than 50 bp, are rare compared to SNVs, many of them result in a pathogenic effect on the encoded protein. Only 3 SVs (all pathogenic) in GLA are currently described in ClinVar, compared to 226 pathogenic SNVs and short indels. Surprisingly, in one sample in our study cohort, we detected a 2914bp deletion between introns 1 and 2, which completely removes exon 2 (Fig. S2). The read length distribution of the amplicon sequencing shows a peak around 10kb, the expected amplicon length of the deletion variant (Fig. S1, sample 11). As in the case of the intronic SNV, the precise breakpoints of this deletion could only be detected using whole gene amplicon sequencing as both of its boundaries are deep intronic.

Variant phasing

One of the advantages of using ONT in this workflow is the fact that > 50% of the reads include the full length gene, allowing haplotype phase determination of variants that originate from a single molecule. As the analysis was performed blindly to the personal background data of the patients, we could detect haplotypes that are shared between two samples (samples 3 and 5, and samples 7 and 8). We predicted that the samples from each of these two sets belonged to patients with blood relationships. At the end of the analysis, when the blinded file was decoded these findings were validated by the fact that the two sets of samples were actually siblings. Moreover, all 12 SNVs/CNVs were confirmed by prior Sanger sequence clinical testing.

Downsampling for designing higher multiplexed sequencing

In this study we sequenced in multiplex 12 samples on one MinION flowcell, yielding an average coverage of 45,000X, much higher than needed for high quality variant calling using ONT reads. In order to evaluate how many samples could be multiplexed in future analyses, we randomly downsampled the reads output to several points. Then, we repeated the analysis workflow on the downsampled reads and evaluated the detection rate of the variants found in the full dataset. Downsampling up to 500 reads per sample allowed detection of all variants from all 12 samples (Fig. 1). While for several samples using only 30 reads were sufficient for 100% true positive detection, for other samples using 250 reads or lower achieved only partial detection, and the rate of several false positive variants passing the quality filter increased. For future GLA genotyping using this pipeline we estimate that 1000 reads per sample will detect all true variants and efficiently discriminate them from low quality false negatives. Thus, for example, sequencing 96 multiplexed samples on a single ONT Flongle flow cell, is estimated to achieve the desired coverage while significantly reducing costs.

Deep intronic and large copy number variant detection using long amplicon sequencing

Two patients in the cohort were of special note: one with the deep intronic variant IV53 + 405T > G and the second one with the deletion of exon 2.

The patient with the variant IV53 + 405T > G suffered acroparesthsias and typical Fabry pain crises and abdominal cramping from childhood. Then, a suspicion of FD was raised, and his enzyme activity was found 8% of normal. His pathogenic variant remained unidentified for 10 years until Sanger sequencing was performed on multiple amplicons spanning the entire non-coding sequence of GLA. The patient with the deletion of exon 2, and his brother and mother were all clinically diagnosed, with Fabry disease. The patient had zero enzyme activity, but the variant was not found for several years by conventional Sanger sequencing. By extracting mRNA and sequencing the cDNA the deletion was found, as mentioned above, with a delay of several years. In contrast, both variants were identified with ease using the ONT long amplicon sequencing method.

All but two variants in Table 2 were classified as pathogenic by a combination of variant effect prediction tools (CADD, REVEL, SIFT, MutationTaser, Polyphen2) and 4 of them were already classified by ClinVar as pathogenic. One patient carried the variant R118C that has conflicting pathogenicity interpretation ^29,30. Indeed, he did not present any of the signs and symptoms that are shown in Table 1. The second variant T194I was classified as likely pathogenic. This patient presented all the signs and symptoms described in Table 1 except for stroke.

The purpose of our study is to validate the technique of Oxford Nanopore long amplicon sequencing of the GLA gene.

Fabry disease is a rare lysosomal genetic disorder and delay in diagnosis between the appearances of the first symptoms until the disease is recognized can span many years ³¹. While measuring alpha galactosidase enzyme activity is accurate in males, in women, there is an overlap between Fabry diseased and healthy females regarding levels of enzyme activity. Furthermore, definitive diagnosis is recommended and it implies finding the genetic variant causing Fabry. A genetic diagnosis is not only important for the patient but also for screening of the entire family. This can be performed only when the disease causing variant of the proband patient is detected.

The GLA protein has a 48767 Da mass and the encoding gene encompasses 7 exons. Approximately 1000 pathogenic variants have been described including, SNVs, small deletions, insertions, CNV and IVS variants.

While in most cases Sanger sequencing is straightforward, some cases can be challenging and remain without a molecular diagnosis for many years. Two examples are described in this paper, in which standard Sanger sequencing of the GLA coding sequence was unable to detect a deep intronic variant in one patient and a complete exon deletion in another.

For these two patients, although a clinical diagnosis was suspected, the final diagnosis based on the genetic variant was delayed by 10 years. While screening in the high risk males from both families could be performed by GLA enzyme activity measurements, lack of a molecular genetic diagnosis precluded genetic screening of the high risk females in each family.

Given that FD is a multi-system disorder, it is often that the coding sequence of GLA is sequenced in order to confirm FD diagnosis. Screening is also necessary for the entire GLA gene sequence because most pathogenic variants in GLA are “private” and therefore not expected to resurface in other patients with similar ancestry. Hence, cost-effective genetic screening methods are necessary in order to conduct FD testing at high scale.

In the present study, 12 samples were long read sequenced at high depth (45,000X). Using downsampling experiments, we demonstrate that such high coverage is unnecessary to enable accurate genotyping of the full gene sequence. Indeed, we propose that 1,000x sequencing depth will be sufficient to genotype any haploid male or diploid female sample across the entire coding and non-coding GLA sequence. Using long amplicon sequencing, this strategy would enable 96 samples to be sequenced at once on a single ONT Flongle flow cell. At this plexity, we expect the cost per sample of ONT to be much lower than that of Sanger sequencing ($10 ONT vs. $70 Sanger cost per sample). This cost savings would also be in addition to the aforementioned enhanced molecular coverage of CNVs and deep intronic variants.

In conclusion, we describe an accurate and cost-effective method for complete gene sequencing of the FD-associated GLA gene. Our new molecular assay is well adapted for high throughput Fabry pathogenic variant screening in the clinic. In addition, we expect our assay to reduce time to FD diagnosis in the long run, due to its enhanced sequencing coverage of the disease-associated coding and non-coding sequences.

This study was conducted in accordance with the principles of the Helsinki Declaration. The research received the approval by the local IRBs of Shaare Zedek Medical Center and Zurich University, Switzerland, All patients signed a written informed consent. All authors have read and approved the manuscript.

Twelve male Fabry patients (age range 31–60 years) were recruited. Patients data included medical history, cardiac, renal, and neurological evaluations. The presence of stroke or TIA (transient ischemic attack) was evaluated during annual examinations by asking the patient and/or using the medical records. Standard transthoracic 2D-echocardiography was routinely performed in all patients. Cardiomyopathy was defined as the presence of Fabry-typical electrocardiogram (ECG)-changes and/or signs of diastolic dysfunction and/or left ventricular hypertrophy on echocardiography or heart MRI. Kidney involvement was defined as either having protein/creatinine coefficient > 0.015 g/mmol and/or estimated glomerular filtration rate (eGFR) according to Chronic Kidney Disease Epidemiology Collaboration (CKD-EPI) formula of < 90 mL/min/1.73 m². Cerebral involvement was defined as the presence of stroke or TIA. For the present study, all clinical and laboratory results were obtained from the patients` medical records.

All variants have been classified as coding for the classic or late onset phenotype based on genotype and residual α-Gal A activity in males and are published in the International Fabry Disease Genotype/Phenotype Database (www.dbFGP.org) and in previous studies ^9,32 [5 6]. The phenotypic assignments of the variants are supported by the clinical manifestations in males, the age of symptoms onset and by in vitro expression assays as reported previously ^33,34.

DNA extraction and amplification. DNA was extracted from 200µl blood samples using FlexiGene DNA kit (Qiagen) according to manufacturer's instructions. The GLA genomic region was amplified by long PCR using LA Taq polymerase (TaKaRa) with the following primers: GLA-Fwd: TTTCTGTTGGTGCTGATATTGCTTGGGAGGGAATAAGCTAGAGCCATCA; GLA-Rev: ACTTGCCTGTCGCTCTATCTTCCTTTGTCAAGCACGCATTTGCCTAGAT. The 5’ ends of each primer included adapter sequences for subsequent priming with PCR Barcoding Kit (ONT; SQK-PBK004) barcoded sequencing primers. Two rounds of long PCR were performed. First, 1.25 units of TaKaRa LA Taq® DNA Polymerase (catalog number: RR002T) were used to amplify a 13kb amplicon (capturing the entire GLA genomic DNA sequence including promoter region) in a 25ul reaction with 100ng input DNA and 200nM each of the aforementioned primers in 1X LA PCR Buffer ll (Mg2 + plus) and 400uM dNTPs. Thermocycling was as follows: 94°C for 1 minute followed by 20 cycles of 98°C for 10 seconds and 68°C for 13 minutes, then 72°C final extension for 10 minutes. PCR products were purified with 0.45X Ampure XP beads (Beckman Coulter), eluted in low TE (10 mM Tris-HCl (pH 8.0), 0.1 mM EDTA), then subjected to another round of PCR using 200nM of barcoded LWB primer pairs from the ONT PCR Barcoding Kit. The second PCR was at 50ul total volume and also included 2.5 units of TaKaRa LA Taq® DNA Polymerase, 1X LA PCR Buffer ll (Mg2 + plus), and 400uM dNTPs. Thermocycling was as follows: 94°C for 1 minute followed by 30 cycles of 98°C for 10 seconds, and 55°C for 30 seconds, and 68°C for 13 minutes; then 72°C final extension for 10 minutes. The second PCR products were purified with 0.45X Ampure XP beads (Beckman Coulter) and eluted in Tris-NaCl (10 mM Tris-HCl (pH 8.0), 5 mM NaCl).

Library preparation and WGS. DNA concentration of final purified PCR products was determined by Qubit BR (Thermo Fischer) measurement and 12 samples were pooled at equimolar concentrations. Subsequently, 50 femtomol of pooled PCR products were loaded onto a MinION flow cell according to the manufacturer’s protocol (PCR Barcoding Kit SQK-PBK004,,Oxford Nanopore Technologies).and sequencing was performed using a MinION device and MinKnow software (MinION Release 19.12.5) for 48 hours.

Bioinformatics. Raw nanopore events (fast5) files were basecalled using command line Guppy (version 3.4.4) and basecalled reads were quality filtered using NanoFilt (version 2.6.0, parameters '-q 5 -l 100 --headcrop 40'). Filtered reads were aligned to the human reference genome (hg19) using minimap2 (version 2.17,³⁵). Short variants were called using nanopolish (version 0.11.3, Oxford Nanopore Technologies). Deletions and insertions larger than 50bp were called using sniffles (version 1.0.11, ³⁶). All variants were phased using Whatshap (version1.0, ³⁷) and annotated with annovar ³⁸ based on several prediction tools (CADD, REVEL, SIFT, MutationTaser, Polyphen2). SAMtools (version 1.9) was used for sorting, indexing, downsampling and calculating depth of bam files.

Authors contributions:

A.N. was in charge of the clinical data. O.M. – preformed the bioinformatics analysis. T.M. and D.A.Z. – planned and performed the experiments. G.A. – designed the project and was in charge of funding. A.N., O.M., D.A.Z. and G.A. – wrote the manuscript.

Conflict of interest

The authors declare that there is no conflict of interest.

Details of funding:

The authors confirm independence from the sponsors; the content of the article has not been influenced by the sponsors.

Details of ethics approval:

The research received IRB approval in both Medical centers from Shaare Zedek Medical Center, Jerusalem Israel and Zurich University, Switzerland

A patient consent statement:

Informed consent was obtained from all patients and is available upon request.

Data availability statement

The data that support the findings of this study are available on request from the corresponding author. The data are not publicly available due to privacy or ethical restrictions.

Brady, R. O. et al. Enzymatic defect in Fabry’s disease. Ceramidetrihexosidase deficiency. N. Engl. J. Med. 276, 1163–1167 (1967).
Ashton-Prolla, P. et al. Fabry disease: comparison of enzymatic, linkage, and mutation analysis for carrier detection in a family with a novel mutation (30delG). Am. J. Med. Genet. 84, 420–424 (1999).
Germain, D. P. Fabry disease. Orphanet J. Rare Dis. 5, 30 (2010).
Knol, I. E. et al. Different phenotypic expression in relatives with fabry disease caused by a W226X mutation. Am. J. Med. Genet. 82, 436–439 (1999).
Koulousios, K. et al. Fabry disease due to D313Y and novel GLA mutations. BMJ Open 7, e017098 (2017).
Shabbeer, J., Yasuda, M., Luca, E. & Desnick, R. J. Fabry disease: 45 novel mutations in the alpha-galactosidase A gene causing the classical phenotype. Mol. Genet. Metab. 76, 23–30 (2002).
Germain, D. P., Shabbeer, J., Cotigny, S. & Desnick, R. J. Fabry disease: twenty novel alpha-galactosidase A mutations and genotype-phenotype correlations in classical and variant phenotypes. Mol. Med. Camb. Mass 8, 306–312 (2002).
Arends, M. et al. Characterization of Classical and Nonclassical Fabry Disease: A Multicenter Study. J. Am. Soc. Nephrol. JASN 28, 1631–1641 (2017).
Nowak, A. et al. Genotype, phenotype and disease severity reflected by serum LysoGb3 levels in patients with Fabry disease. Mol. Genet. Metab. 123, 148–153 (2018).
Ortiz, A. et al. Fabry disease revisited: Management and treatment recommendations for adult patients. Mol. Genet. Metab. 123, 416–427 (2018).
Desnick, R. J. Fabry Disease (α-Galactosidase A Deficiency). in Brenner’s Encyclopedia of Genetics: Second Edition (ed. Valle, D.) 8–11 ({McGraw}-Hill, 2013). doi:10.1016/B978-0-12-374984-0.00516-7.
Bae, E. H. et al. A late-onset male Fabry disease patient with somatic mosaicism of a classical GLA mutation: a case report. Ann. Palliat. Med. (2020) doi:10.21037/apm-19-635.
Jain, M., Olsen, H. E., Paten, B. & Akeson, M. The Oxford Nanopore MinION: delivery of nanopore sequencing to the genomics community. Genome Biol. 17, 239 (2016).
Kono, N. & Arakawa, K. Nanopore sequencing: review of potential applications in functional genomics. Dev. Growth Differ. 61, 316–326 (2019).
Brown, C. G. & Clarke, J. Nanopore development at Oxford Nanopore. Nat. Biotechnol. 34, 810–811 (2016).
Jain, M. et al. Linear assembly of a human centromere on the Y chromosome. Nat. Biotechnol. 36, 321–323 (2018).
Euskirchen, P. et al. Same-day genomic and epigenomic diagnosis of brain tumors using real-time nanopore sequencing. Acta Neuropathol. (Berl.) 134, 691–703 (2017).
Gigante, S. et al. Using long-read sequencing to detect imprinted DNA methylation. Nucleic Acids Res. 47, e46–e46 (2019).
Clark, M. B. et al. Long-read sequencing reveals the complex splicing profile of the psychiatric risk gene CACNA1C in human brain. Mol. Psychiatry 25, 37–47 (2020).
Charalampous, T. et al. Nanopore metagenomics enables rapid clinical diagnosis of bacterial lower respiratory infection. Nat. Biotechnol. 37, 783–792 (2019).
Moon, J. et al. Rapid diagnosis of bacterial meningitis by nanopore 16S amplicon sequencing: a pilot study. Int. J. Med. Microbiol. 309, 151338 (2019).
Quick, J. et al. Real-time, portable genome sequencing for Ebola surveillance. Nature 530, 228–232 (2016).
Sakai, J. et al. An identification protocol for ESBL-producing Gram-negative bacteria bloodstream infections using a MinION nanopore sequencer. J. Med. Microbiol. 68, 1219–1226 (2019).
Bull, R. A. et al. Analytical validity of nanopore sequencing for rapid SARS-CoV-2 genome analysis. Nat. Commun. 11, 6272 (2020).
James, P. et al. LamPORE: rapid, accurate and highly scalable molecular screening for SARS-CoV-2 infection, based on nanopore sequencing. medRxiv 2020.08.07.20161737 (2020) doi:10.1101/2020.08.07.20161737.
Leija-Salazar, M. et al. Evaluation of the detection of GBA missense mutations and other variants using the Oxford Nanopore MinION. Mol. Genet. Genomic Med. 7, e564 (2019).
Mahaweni, N. M. et al. A comprehensive overview of FCGR3A gene variability by full-length gene sequencing including the identification of V158F polymorphism. Sci. Rep. 8, 15983 (2018).
Minervini, C. F. et al. TP53 gene mutation analysis in chronic lymphocytic leukemia by nanopore MinION sequencing. Diagn. Pathol. 11, 96 (2016).
Ferreira, S. et al. The alpha-galactosidase A p.Arg118Cys variant does not cause a Fabry disease phenotype: data from individual patients and family studies. Mol. Genet. Metab. 114, 248–258 (2015).
Talbot, A. & Nicholls, K. Elevated Lyso-Gb3 Suggests the R118C GLA Mutation Is a Pathological Fabry Variant. JIMD Rep. 45, 95–98 (2019).
Reisin, R., Perrin, A. & García-Pavía, P. Time delays in the diagnosis and treatment of Fabry disease. Int. J. Clin. Pract. 71, (2017).
Nowak, A., Mechtler, T., Kasper, D. C. & Desnick, R. J. Correlation of Lyso-Gb3 levels in dried blood spots and sera from patients with classic and Later-Onset Fabry disease. Mol. Genet. Metab. 121, 320–324 (2017).
Benjamin, E. R. et al. The validation of pharmacogenetics for the identification of Fabry patients to be treated with migalastat. Genet. Med. Off. J. Am. Coll. Med. Genet. 19, 430–438 (2017).
Yasuda, M., Shabbeer, J., Osawa, M. & Desnick, R. J. Fabry disease: novel alpha-galactosidase A 3’-terminal mutations result in multiple transcripts due to aberrant 3’-end formation. Am. J. Hum. Genet. 73, 162–173 (2003).
Li, H. Minimap2: pairwise alignment for nucleotide sequences. Bioinformatics 34, 3094–3100 (2018).
Sedlazeck, F. J. et al. Accurate detection of complex structural variations using single-molecule sequencing. Nat. Methods 15, 461–468 (2018).
Martin, M. et al. WhatsHap: fast and accurate read-based phasing. bioRxiv 085050 (2016) doi:10.1101/085050.
Wang, K., Li, M. & Hakonarson, H. ANNOVAR: functional annotation of genetic variants from high-throughput sequencing data. Nucleic Acids Res. 38, e164–e164 (2010).

Table 1: Clinical presentation of the cohort of Fabry patients. Acrop - acroparestesias, Hypo - hypohidriosis, Angio - angioceratoma, Cornea - cornea verticillata, Nephro – nephropathy, Cardio – cardiomyopathy, Stroke – history of stroke.

Sample	Therapy	Acrop	Hypo	Angio	Cornea	Nephro	Cardio	Stroke	phenotype	Note
1	Agalsidase α	+							classic
2	Agalsidase β	+	+	+	+				classic
3	Agalsidase β	+	+		+	+			classic
4	Agalsidase β	+	+	+	+				classic
5	Agalsidase β	+	+	+	+	+	+		classic
6	No therapy								benign
7	Agalsidase α	+	+	+	+	+			classic	brother of 8
8	Agalsidase α	+	+	+	+	+	+	+	classic	brother of 7
9	Agalsidase α	+	+		+	+	+		classic
10	Agalsidase β						+		later-onset
11	Agalsidase α	+	+	+	+	+	+		classic
12	Agalsidase α	+	+	+	+	+	+		classic

Table 2: Genotyping results based on ONT amplicon sequencing, and predicted consequence of the genetic variant on GLA. * - two samples from the same family were tested; N/A =not applicable

ONT Barcode	*GLA* nucleotide variant (Accession: NM_000169)	GLA protein variant (Accession: NP_000160)	Comments
1	IVS53+405T>G	N/A-deep intronic	deep intronic variant; possible branch site; may affect splicing
2	c.1147_1149delTTC	p.Phe383del	pathogenic
3	c.744_745delTA*	p.Phe248LeufsX7	pathogenic; same haplotype as barcode05 suggests blood relationship
4	c.559_560delAT	p.Met187Valfs*6	pathogenic
5	c.744_745delTA*	p.Phe248LeufsX7	pathogenic; same haplotype as barcode03 suggests blood relationship
6	c.352C>T	p.Arg118Cys	Conflicting_interpretations_of_pathogenicity
7	c.370-2A>G*	N/A- splicing	pathogenic; same haplotype as barcode08 suggests blood relationship
8	c.370-2A>G*	N/A- splicing	pathogenic; same haplotype as barcode07 suggests blood relationship
9	c.704C>A	p.Ser235Tyr	pathogenic
10	c.337T>C	p. Phe113Leu	pathogenic
11	Exon 2 deletion)	N/A	pathogenic; 2914bp deletion removes GLA exon2 (chrX:100658307-100661221
12	c.581C>T	p.Thr194Ile	likely pathogenic

No competing interests reported.

GLAONTsupportinginformation.pdf
Table S1. Sequencing yield statistics. Fig S1. Size distribution of 5000 randomly subsampled ONT reads from each sample, showing a peak at 13 kb, the full PCR product length. In sample 11 which possess a 2.9 kb deletion, the peak is at 10 kb. Fig S2. IGV session showing the 2914 bp deletion in sample 11.

Download PDF

Journal Publication

published 16 Nov, 2021

Read the published version in Scientific Reports →

Editorial decision: Major revision
20 Jul, 2021
Reviews received at journal
30 Jun, 2021
Reviewers agreed at journal
09 Jun, 2021
Reviewers invited by journal
09 Jun, 2021
Editor assigned by journal
02 Jun, 2021
Editor invited by journal
02 Jun, 2021
Submission checks completed at journal
02 Jun, 2021
First submitted to journal
01 Jun, 2021

You are reading this latest preprint version

Detection of Single Nucleotide and Copy Number Variants in the Fabry Disease-associated GLA Gene Using Nanopore Sequencing

Status:

Journal Publication

Version 1

Abstract

Figures

Introduction:

Results:

Discussion:

Material And Methods:

Declarations

References

Tables

Additional Declarations

Supplementary Files

Status:

Journal Publication

Version 1