Genome-wide sequencing as a first-tier screening test for short tandem repeat expansions

Abstract

Background Screening for short tandem repeat (STR) expansions in next-generation sequencing data can enable diagnosis, optimal clinical management/treatment, and accurate genetic counseling of patients with repeat expansion disorders. We aimed to develop an efficient computational workflow for reliable detection of STR expansions in next-generation sequencing data and demonstrate its clinical utility. Methods We characterized the performance of eight STR analysis methods (lobSTR, HipSTR, RepeatSeq, ExpansionHunter, TREDPARSE, GangSTR, STRetch, and exSTRa) on next-generation sequencing datasets of samples with known disease-causing full-mutation STR expansions and genomes simulated to harbor repeat expansions at selected loci and optimized their sensitivity. We then used a machine learning decision tree classifier to identify an optimal combination of methods for full-mutation detection. In Burrows-Wheeler Aligner (BWA)-aligned genomes, the ensemble approach of using ExpansionHunter, STRetch, and exSTRa performed the best (precision = 82%, recall = 100%, F1-score = 90%). We applied this pipeline to screen 301 families of children with suspected genetic disorders. Results We identified 10 individuals with full-mutations in the AR, ATXN1, ATXN8, DMPK, FXN, or HTT disease STR locus in the analyzed families. Additional candidates identified in our analysis include two probands with borderline ATXN2 expansions between the established repeat size range for reduced-penetrance and full-penetrance full-mutation and seven individuals with FMR1 CGG repeats in the intermediate/premutation repeat size range. In 67 probands with a prior negative clinical PCR test for the FMR1, FXN, or DMPK disease STR locus, or the spinocerebellar ataxia disease STR panel, our pipeline did not falsely identify aberrant expansion. We performed clinical PCR tests on seven (out of 10) full-mutation samples identified by our pipeline and confirmed the expansion status in all, showing absolute concordance between our bioinformatics and molecular findings. Conclusions We have successfully demonstrated the application of a well-optimized bioinformatics pipeline that promotes the utility of genome-wide sequencing as a first-tier screening test to detect expansions of known disease STRs. Interrogating clinical next-generation sequencing data for pathogenic STR expansions using our ensemble pipeline can improve diagnostic yield and enhance clinical outcomes for patients with repeat expansion disorders.

Bulked Segregant RNA Sequencing Revealed Difference Between Virulent and Avirulent Brown Planthoppers

“…The reads were mapped to the latest chromosome-level BPH reference genome (BioProject accession no. PRJNA591478 ) ( Ma et al, 2021 ) using Burrows–Wheeler Aligner (BWA, v. 0.7.5a) ( Rajan-Babu et al, 2021 ). The bwa index command and BWA–MEM were used to create reference genome index and map the reference genome, respectively.…”

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

“…Hence, PacBio CCS platforms allow the utilization of smaller fragment sizes, with a global quality higher than Illumina ( 109 ). To complicate things further, improved analysis of Illumina short read data is now possible thanks to new algorithms like the ExpansionHunter and GangSTR ( 68 , 114 , 115 ). To the best of our knowledge, a systematic comparison of these short-read analyses to LRS is not available.…”

Lessons learned: next-generation sequencing applied to undiagnosed genetic diseases

“…Reliable determination of RE length and composition remains a problem, and current RE detection often still relies on PCR or Southern blot techniques ( Tables 1 and 2 , Figure 1 , and ref. 44 ).…”