Accurate repeat expansion genotyping from whole-genome and real-time targeted nanopore sequencing with Mahi-mahi
Gibson, S. B.; Palmer, R. N.; Hollizeck, S.; Koole, W.; Jeffery, H. M.; Albring, M.; Buchan, J. G.; Cheetham, R. K.; Elferink, M. G.; Goffena, J.; Hamlin, E. J.; Jansen, R.; Lockhart, P. J.; Martin, R.; Miller, A. L.; van Outersterp, I.; Prall, T.; Rosinski, J.; Rowley, S. M.; Sadedin, S.; Sewell, D.; Todd, M.; Zabelkin, A.; Zalusky, M. P. G.; van der Zwaag, B.; Nurk, S.; Juul, S.; van Gassen, K. L. I.; Lunke, S.; Miller, D. E.; Rescheneder, P.
Show abstract
Tandem repeat expansions cause more than 60 Mendelian disorders but are difficult to detect and size accurately with short-read sequencing because of their length and repetitive structure. Long-read sequencing (LRS) resolves these loci directly, yet analytic tooling for nanopore sequencing data has lagged behind that available for other long-read platforms. Here we present Mahi-mahi, a tandem repeat genotyper designed for nanopore reads that reports allele-level repeat length, motif composition, interruptions, and read-level flags for systematic basecalling artifacts. Across six Genome in a Bottle samples, whole-genome and adaptive sampling (real-time targeted) data returned calls within one repeat unit of each other for 829 of 832 alleles, and Mahi-mahi calls were within one repeat unit of orthogonal HiFi data at 95.6% of compared alleles after excluding three loci at which the two tools define the repeat differently. In 245 samples from the 1000 Genomes Project, Mahi-mahi calls were within one repeat unit of assembly-derived calls for 98.6% of 10,927 alleles at 57 loci, and Mahi-mahi distinguished pathogenic from nonpathogenic RFC1 alleles of overlapping size by motif. Across 14 Coriell reference lines, Mahi-mahi identified the expected expansion in every line and agreed with orthogonal assays to within two repeat units at HTT, ATXN1, DMPK, and PABPN1. Finally, in 103 samples from three laboratories, Mahi-mahi detected an expanded allele in all samples in which prior testing had reported a pathogenic expansion, assigned the same category as the referring laboratory for 98% of 180 genotypes, and identified known ATXN1 interruption structure. Our results demonstrate that Mahi-mahi accurately identifies repeat expansions at disease-associated loci and could serve as a single assay to evaluate repeat expansions.
Matching journals
The top 3 journals account for 50% of the predicted probability mass.
Similar papers in this journal
- Genome-Wide Sequencing as a First-Tier Screening Test for Short Tandem Repeat Expansions 95%
- Evaluating Genome Sequencing Strategies: Trio, Singleton, and Standard Testing in Rare Disease Diagnosis 93%
- COBT: A gene-based rare variant burden test for case-only study designs using aggregated genotypes from public reference cohorts. 93%
Similar papers in this journal
- The Importance of Automation in Genetic Diagnosis: Lessons from Analyzing an Inherited Retinal Degeneration Cohort with the Mendelian Analysis Toolkit (MATK) 95%
- IGenomic answers for children: Dynamic analyses of >1000 pediatric rare disease genomes 94%
- One in seven pathogenic variants can be challenging to detect by NGS: An analysis of 450,000 patients with implications for clinical sensitivity and genetic test implementation 93%
Similar papers in this journal
- Long-read genome sequencing for the diagnosis of neurodevelopmental disorders 94%
- IMPROVE-DD: Integrating Multiple Phenotype Resources Optimises Variant Evaluation in genetically determined Developmental Disorders 91%
- Leveraging TOPMed Imputation Server and Constructing a Cohort-Specific Imputation Reference Panel to Enhance Genotype Imputation among Cystic Fibrosis Patients 91%
Similar papers in this journal
- Diagnostic Utility of Genome-wide DNA Methylation Analysis in Genetically Unsolved Developmental and Epileptic Encephalopathies and Refinement of a CHD2 Episignature 92%
- K-mer analysis of long-read alignment pileups for structural variant genotyping 92%
- A comprehensive tandem repeat catalog of the human genome 92%
"Similar papers" are the closest papers from that journal in the model's embedding space. They show what the match is built on, but the ranking comes mostly from a classifier over the whole training set, not from these examples alone.