SMRT Sequencing of Long Tandem Nucleotide Repeats in SCA10 Reveals Unique Insight of Repeat Expansion Structure

Short Tandem Repeats in the era of next-generation sequencing: from historical loci to population databases

In this study, we explore the landscape of short tandem repeats (STRs) within the human genome through the lens of evolving technologies to detect genomic variations. STRs, which encompass approximately 3% of our genomic DNA, are crucial for understanding human genetic diversity, disease mechanisms, and evolutionary biology. The advent of high-throughput sequencing methods has revolutionized our ability to accurately map and analyze STRs, highlighting their significance in genetic disorders, forensic science, and population genetics. We review the current available methodologies for STR analysis, the challenges in interpreting STR variations across different populations, and the implications of STRs in medical genetics. Our findings underscore the urgent need for comprehensive STR databases that reflect the genetic diversity of global populations, facilitating the interpretation of STR data in clinical diagnostics, genetic research, and forensic applications. This work sets the stage for future studies aimed at harnessing STR variations to elucidate complex genetic traits and diseases, reinforcing the importance of integrating STRs into genetic research and clinical practice.

Sequence composition changes in short tandem repeats: heterogeneity, detection, mechanisms and clinical implications

Short tandem repeats (STRs) are a class of repetitive elements, composed of tandem arrays of 1-6 base pair sequence motifs, that comprise a substantial fraction of the human genome. STR expansions can cause a wide range of neurological and neuromuscular conditions, known as repeat expansion disorders, whose age of onset, severity, penetrance and/or clinical phenotype are influenced by the length of the repeats and their sequence composition. The presence of non-canonical motifs, depending on the type, frequency and position within the repeat tract, can alter clinical outcomes by modifying somatic and intergenerational repeat stability, gene expression and mutant transcript-mediated and/or protein-mediated toxicities. Here, we review the diverse structural conformations of repeat expansions, technological advances for the characterization of changes in sequence composition, their clinical correlations and the impact on disease mechanisms.

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.

A long-read sequencing strategy with overlapping linkers on adjacent fragments (OLAF-Seq) for targeted resequencing and enrichment

In this report, we present OLAF-Seq, a novel strategy to construct a long-read sequencing library such that adjacent fragments are linked with end-terminal duplications. We use the CRISPR-Cas9 nickase enzyme and a pool of multiple sgRNAs to perform non-random fragmentation of targeted long DNA molecules (> 300kb) into smaller library-sized fragments (about 20 kbp) in a manner so as to retain physical linkage information (up to 1000 bp) between adjacent fragments. DNA molecules targeted for fragmentation are preferentially ligated with adaptors for sequencing, so this method can enrich targeted regions while taking advantage of the long-read sequencing platforms. This enables the sequencing of target regions with significantly lower total coverage, and the genome sequence within linker regions provides information for assembly and phasing. We demonstrated the validity and efficacy of the method first using phage and then by sequencing a panel of 100 full-length cancer-related genes (including both exons and introns) in the human genome. When the designed linkers contained heterozygous genetic variants, long haplotypes could be established. This sequencing strategy can be readily applied in both PacBio and Oxford Nanopore platforms for both long and short genes with an easy protocol. This economically viable approach is useful for targeted enrichment of hundreds of target genomic regions and where long no-gap contigs need deep sequencing.

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

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

A Chinese SCA36 pedigree analysis of NOP56 expansion region based on long-read sequencing

Introduction: Spinocerebellar ataxias 36 (SCA36) is the neurodegenerative disease caused by the GGCCTG Hexanucleotide repeat expansions in NOP56, which is too long to sequence using short-read sequencing. Single molecule real time (SMRT) sequencing can sequence across disease-causing repeat expansion. We report the first long-read sequencing data across the expansion region in SCA36.

Methods: We collected and described the clinical manifestations and imaging features of Han Chinese pedigree with three generations of SCA36. Also, we focused on structural variation analysis for intron 1 of the NOP56 gene by SMRT sequencing in the assembled genome.

Results: The main clinical features of this pedigree are late-onset ataxia symptoms, with a presymptomatic presence of affective and sleep disorders. In addition, the results of SMRT sequencing showed the specific repeat expansion region and demonstrated that the region was not composed of single GGCCTG hexanucleotides and there were random interruptions.

Discussion: We extended the phenotypic spectrum of SCA36. We applied SMRT sequencing to reveal the correlation between genotype and phenotype of SCA36. Our findings indicated that long-read sequencing is well suited to characterize known repeat expansion.

Thymidine Kinase 2 and Mitochondrial Protein COX I in the Cerebellum of Patients with Spinocerebellar Ataxia Type 31 Caused by Penta-nucleotide Repeats (TTCCA)n

Spinocerebellar ataxia type 31 (SCA31), an autosomal-dominant neurodegenerative disorder characterized by progressive cerebellar ataxia with Purkinje cell degeneration, is caused by a heterozygous 2.5-3.8 kilobase penta-nucleotide repeat of (TTCCA)_n in intron 11 of the thymidine kinase 2 (TK2) gene. TK2 is an essential mitochondrial pyrimidine-deoxyribonucleoside kinase. Bi-allelic loss-of-function mutations of TK2 lead to mitochondrial DNA depletion syndrome (MDS) in humans through severe (~ 70%) reduction of mitochondrial electron-transport-chain activity, and tk2 knockout mice show Purkinje cell degeneration and ataxia through severe mitochondrial cytochrome-c oxidase subunit I (COX I) protein reduction. To clarify whether TK2 function is altered in SCA31, we investigated TK2 and COX I expression in human postmortem SCA31 cerebellum. We confirmed that canonical TK2 mRNA is transcribed from exons far upstream of the repeat site, and demonstrated that an extended version of TK2 mRNA ("TK2-EXT"), transcribed from exons spanning the repeat site, is expressed in human cerebellum. While canonical TK2 was conserved among vertebrates, TK2-EXT was specific to primates. Reverse transcription-PCR demonstrated that both TK2 mRNAs were preserved in SCA31 cerebella compared with control cerebella. The TK2 proteins, assessed with three different antibodies including our original polyclonal antibody against TK2-EXT, were detected as ~ 26 kilodalton proteins on western blot; their levels were similar in SCA31 and control cerebella. COX I protein level was preserved in SCA31 compared to nuclear DNA-encoded protein. We conclude that the expression and function of TK2 are preserved in SCA31, suggesting a mechanism distinct from that of MDS.

Repeat Detector: versatile sizing of expanded tandem repeats and identification of interrupted alleles from targeted DNA sequencing

Targeted DNA sequencing approaches will improve how the size of short tandem repeats is measured for diagnostic tests and preclinical studies. The expansion of these sequences causes dozens of disorders, with longer tracts generally leading to a more severe disease. Interrupted alleles are sometimes present within repeats and can alter disease manifestation. Determining repeat size mosaicism and identifying interruptions in targeted sequencing datasets remains a major challenge. This is in part because standard alignment tools are ill-suited for repetitive and unstable sequences. To address this, we have developed Repeat Detector (RD), a deterministic profile weighting algorithm for counting repeats in targeted sequencing data. We tested RD using blood-derived DNA samples from Huntington's disease and Fuchs endothelial corneal dystrophy patients sequenced using either Illumina MiSeq or Pacific Biosciences single-molecule, real-time sequencing platforms. RD was highly accurate in determining repeat sizes of 609 blood-derived samples from Huntington's disease individuals and did not require prior knowledge of the flanking sequences. Furthermore, RD can be used to identify alleles with interruptions and provide a measure of repeat instability within an individual. RD is therefore highly versatile and may find applications in the diagnosis of expanded repeat disorders and in the development of novel therapies.

Long-read sequencing for molecular diagnostics in constitutional genetic disorders

Long-read sequencing (LRS) has been around for more than a decade, but widespread adoption of the technology has been slow due to the perceived high error rates and high sequencing cost. This is changing due to the recent advancements to produce highly accurate sequences and the reducing costs. LRS promises significant improvement over short read sequencing in four major areas: (1) better detection of structural variation (2) better resolution of highly repetitive or nonunique regions (3) accurate long-range haplotype phasing and (4) the detection of base modifications natively from the sequencing data. Several successful applications of LRS have demonstrated its ability to resolve molecular diagnoses where short-read sequencing fails to identify a cause. However, the argument for increased diagnostic yield from LRS remains to be validated. Larger cohort studies may be required to establish the realistic boundaries of LRS's clinical utility and analytical validity, as well as the development of standards for clinical applications. We discuss the limitations of the current standard of care, and contrast with the applications and advantages of two major LRS platforms, PacBio and Oxford Nanopore, for molecular diagnostics of constitutional disorders, and present a critical argument about the potential of LRS in diagnostic settings.

ATTCT and ATTCC repeat expansions in the ATXN10 gene affect disease penetrance of spinocerebellar ataxia type 10

Spinocerebellar ataxia type 10 (SCA10) is an autosomal-dominant disorder caused by an expanded pentanucleotide repeat in the ATXN10 gene. This repeat expansion, when fully penetrant, has a size of 850–4,500 repeats. It has been shown that the repeat composition can be a modifier of disease, e.g., seizures. Here, we describe a Mexican kindred in which we identified both pure (ATTCT)_n and mixed (ATTCT)_n-(ATTCC)_n expansions in the same family. We used amplification-free targeted sequencing and optical genome mapping to decipher the composition of these repeat expansions. We found a considerable degree of mosaicism of the repeat expansion. This mosaicism was confirmed in skin fibroblasts from individuals with ATXN10 expansions with RNAScope in situ hybridization. All affected family members with the mixed ATXN10 repeat expansion showed typical clinical signs of spinocerebellar ataxia and epilepsy. In contrast, individuals with the pure ATXN10 expansion present with Parkinson's disease or are unaffected, even in individuals more than 20 years older than the average age at onset for SCA10. Our findings suggest that the pure (ATTCT)_n expansion is non-pathogenic, while repeat interruptions, e.g., (ATTCC)_n, are necessary to cause SCA10. This mechanism has been recently described for several other repeat expansions including SCA31 (BEAN1), SCA37 (DAB1), and three loci for benign adult familial myoclonic epilepsy BAFME (SAMD12, TNRC6A, RAPGEF2). Therefore, long-read sequencing and optical genome mapping of the entire genomic structure of repeat expansions are critical for clinical practice and genetic counseling, as variations in the repeat can affect disease penetrance, symptoms, and disease trajectory.

The genetic and molecular features of the intronic pentanucleotide repeat expansion in spinocerebellar ataxia type 10

Spinocerebellar ataxia type 10 (SCA10) is characterized by progressive cerebellar neurodegeneration and, in many patients, epilepsy. This disease mainly occurs in individuals with Indigenous American or East Asian ancestry, with strong evidence supporting a founder effect. The mutation causing SCA10 is a large expansion in an ATTCT pentanucleotide repeat in intron 9 of the ATXN10 gene. The ATTCT repeat is highly unstable, expanding to 280–4,500 repeats in affected patients compared with the 9–32 repeats in normal individuals, one of the largest repeat expansions causing neurological disorders identified to date. However, the underlying molecular basis of how this huge repeat expansion evolves and contributes to the SCA10 phenotype remains largely unknown. Recent progress in next-generation DNA sequencing technologies has established that the SCA10 repeat sequence has a highly heterogeneous structure. Here we summarize what is known about the structure and origin of SCA10 repeats, discuss the potential contribution of variant repeats to the SCA10 disease phenotype, and explore how this information can be exploited for therapeutic benefit.

Detection of rare thalassemia mutations using long-read single-molecule real-time sequencing

Gap- polymerase chain reaction (PCR), reverse dot-blot assay (RDB), real-time PCR based multicolor melting curve analysis (MMCA assay), multiplex ligation-dependent probe amplification (MLPA) and Sanger sequencing are conventional methods to diagnose thalassemia but all of them have limitations. In this study, we applied single-molecule real-time (SMRT) sequencing following multiplex long-range PCR to uncover rare mutations in nine patients and their family members. The patients with different results between Gap-PCR and MMCA assay or with phenotype not matching genotype were included. Using SMRT sequencing, we first identified the carriers with ααα^anti3.7/HKαα, -α^762bpα/αα (chr16:172,648-173,409), αα^fusion/α^QSα (in a trans configuration), two cases with novel gene rearrangements and another case with a novel 341 bp insertion in α-globin gene cluster, respectively. One carrier with --^SEA/ααα^anti4.2, and two carriers with the coexistence of globin variant and an α-globin gene duplication were also found. Most importantly, we could determine two defects in α-globin gene cluster being a cis or trans configuration in a single test. Our results showed that SMRT has great advantages in detection of α-globin gene triplications, rare deletions and determination of a cis or trans configuration. SMRT is a comprehensive and one-step method for thalassemia screening and diagnosis, especially for detection of rare thalassemia mutations.

Mechanistic and Therapeutic Insights into Ataxic Disorders with Pentanucleotide Expansions

Pentanucleotide expansion diseases constitute a special class of neurodegeneration. The repeat expansions occur in non-coding regions, have likely arisen from Alu elements, and often result in autosomal dominant or recessive phenotypes with underlying cerebellar neuropathology. When transcribed (potentially bidirectionally), the expanded RNA forms complex secondary and tertiary structures that can give rise to RNA-mediated toxicity, including protein sequestration, pentapeptide synthesis, and mRNA dysregulation. Since several of these diseases have recently been discovered, our understanding of their pathological mechanisms is limited, and their therapeutic interventions underexplored. This review aims to highlight new in vitro and in vivo insights into these incurable diseases.

Third-Generation Cytogenetic Analysis: Diagnostic Application of Long-Read Sequencing

Copy number variants (CNVs) play important roles in the pathogenesis of several genetic syndromes. Traditional and molecular karyotyping are considered the first-tier diagnostic tests to detect macroscopic and cryptic deletions/duplications. However, their time-consuming and laborious experimental protocols protract diagnostic times from 3 to 15 days. Nanopore sequencing has the ability to reduce time to results for the detection of CNVs with the same resolution of current state-of-the-art diagnostic tests. Nanopore sequencing was compared to molecular karyotyping for the detection of pathogenic CNVs of seven patients with previously diagnosed causative CNVs of different sizes and cellular fractions. Larger chromosomal anomalies included trisomy 21 and mosaic tetrasomy 12p. Among smaller CNVs, two genomic imbalances of 1.3 Mb, a small deletion of 170 kb, and two mosaic deletions (1.2 Mb and 408 kb) were tested. DNA was sequenced and data generated during runs were analyzed in online mode. All pathogenic CNVs were identified with detection time inversely proportional to size and cellular fraction. Aneuploidies were called after only 30 minutes of sequencing, whereas 30 hours were needed to call small CNVs. These results demonstrate the clinical utility of our approach that allows the molecular diagnosis of genomic disorders within a 30-minute to 30-hour time frame and its easy implementation as a routinary diagnostic tool.

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.

Molecular Mechanisms in Pentanucleotide Repeat Diseases

The number of neurodegenerative diseases resulting from repeat expansion has increased extraordinarily in recent years. In several of these pathologies, the repeat can be transcribed in RNA from both DNA strands producing, at least, one toxic RNA repeat that causes neurodegeneration by a complex mechanism. Recently, seven diseases have been found caused by a novel intronic pentanucleotide repeat in distinct genes encoding proteins highly expressed in the cerebellum. These disorders are clinically heterogeneous being characterized by impaired motor function, resulting from ataxia or epilepsy. The role that apparently normal proteins from these mutant genes play in these pathologies is not known. However, recent advances in previously known spinocerebellar ataxias originated by abnormal non-coding pentanucleotide repeats point to a gain of a toxic function by the pathogenic repeat-containing RNA that abnormally forms nuclear foci with RNA-binding proteins. In cells, RNA foci have been shown to be formed by phase separation. Moreover, the field of repeat expansions has lately achieved an extraordinary progress with the discovery that RNA repeats, polyglutamine, and polyalanine proteins are crucial for the formation of nuclear membraneless organelles by phase separation, which is perturbed when they are expanded. This review will cover the amazing advances on repeat diseases.

Characterization of FMR1 Repeat Expansion and Intragenic Variants by Indirect Sequence Capture

Traditional methods for the analysis of repeat expansions, which underlie genetic disorders, such as fragile X syndrome (FXS), lack single-nucleotide resolution in repeat analysis and the ability to characterize causative variants outside the repeat array. These drawbacks can be overcome by long-read and short-read sequencing, respectively. However, the routine application of next-generation sequencing in the clinic requires target enrichment, and none of the available methods allows parallel analysis of long-DNA fragments using both sequencing technologies. In this study, we investigated the use of indirect sequence capture (Xdrop technology) coupled to Nanopore and Illumina sequencing to characterize FMR1, the gene responsible of FXS. We achieved the efficient enrichment (> 200×) of large target DNA fragments (~60-80 kbp) encompassing the entire FMR1 gene. The analysis of Xdrop-enriched samples by Nanopore long-read sequencing allowed the complete characterization of repeat lengths in samples with normal, pre-mutation, and full mutation status (> 1 kbp), and correctly identified repeat interruptions relevant for disease prognosis and transmission. Single-nucleotide variants (SNVs) and small insertions/deletions (indels) could be detected in the same samples by Illumina short-read sequencing, completing the mutational testing through the identification of pathogenic variants within the FMR1 gene, when no typical CGG repeat expansion is detected. The study successfully demonstrated the parallel analysis of repeat expansions and SNVs/indels in the FMR1 gene at single-nucleotide resolution by combining Xdrop enrichment with two next-generation sequencing approaches. With the appropriate optimization necessary for the clinical settings, the system could facilitate both the study of genotype-phenotype correlation in FXS and enable a more efficient diagnosis and genetic counseling for patients and their relatives.

Deciphering Neurodegenerative Diseases Using Long-Read Sequencing

Neurodegenerative diseases exhibit chronic progressive lesions in the central and peripheral nervous systems with unclear causes. The search for pathogenic mutations in human neurodegenerative diseases has benefited from massively parallel short-read sequencers. However, genomic regions, including repetitive elements, especially with high/low GC content, are far beyond the capability of conventional approaches. Recently, long-read single-molecule DNA sequencing technologies have emerged and enabled researchers to study genomes, transcriptomes, and metagenomes at unprecedented resolutions. The identification of novel mutations in unresolved neurodegenerative disorders, the characterization of causative repeat expansions, and the direct detection of epigenetic modifications on naive DNA by virtue of long-read sequencers will further expand our understanding of neurodegenerative diseases. In this article, we review and compare 2 prevailing long-read sequencing technologies, Pacific Biosciences and Oxford Nanopore Technologies, and discuss their applications in neurodegenerative diseases.

An update on the neurological short tandem repeat expansion disorders and the emergence of long-read sequencing diagnostics

BACKGROUND

Short tandem repeat (STR) expansion disorders are an important cause of human neurological disease. They have an established role in more than 40 different phenotypes including the myotonic dystrophies, Fragile X syndrome, Huntington's disease, the hereditary cerebellar ataxias, amyotrophic lateral sclerosis and frontotemporal dementia.

MAIN BODY

STR expansions are difficult to detect and may explain unsolved diseases, as highlighted by recent findings including: the discovery of a biallelic intronic 'AAGGG' repeat in RFC1 as the cause of cerebellar ataxia, neuropathy, and vestibular areflexia syndrome (CANVAS); and the finding of 'CGG' repeat expansions in NOTCH2NLC as the cause of neuronal intranuclear inclusion disease and a range of clinical phenotypes. However, established laboratory techniques for diagnosis of repeat expansions (repeat-primed PCR and Southern blot) are cumbersome, low-throughput and poorly suited to parallel analysis of multiple gene regions. While next generation sequencing (NGS) has been increasingly used, established short-read NGS platforms (e.g., Illumina) are unable to genotype large and/or complex repeat expansions. Long-read sequencing platforms recently developed by Oxford Nanopore Technology and Pacific Biosciences promise to overcome these limitations to deliver enhanced diagnosis of repeat expansion disorders in a rapid and cost-effective fashion.

CONCLUSION

We anticipate that long-read sequencing will rapidly transform the detection of short tandem repeat expansion disorders for both clinical diagnosis and gene discovery.

Founder Effects of Spinocerebellar Ataxias in the American Continents and the Caribbean

Spinocerebellar ataxias (SCAs) comprise a heterogeneous group of autosomal dominant disorders. The relative frequency of the different SCA subtypes varies broadly among different geographical and ethnic groups as result of genetic drifts. This review aims to provide an update regarding SCA founders in the American continents and the Caribbean as well as to discuss characteristics of these populations. Clusters of SCAs were detected in Eastern regions of Cuba for SCA2, in South Brazil for SCA3/MJD, and in Southeast regions of Mexico for SCA7. Prevalence rates were obtained and reached 154 (municipality of Báguano, Cuba), 166 (General Câmara, Brazil), and 423 (Tlaltetela, Mexico) patients/100,000 for SCA2, SCA3/MJD, and SCA7, respectively. In contrast, the scattered families with spinocerebellar ataxia type 10 (SCA10) reported all over North and South Americas have been associated to a common Native American ancestry that may have risen in East Asia and migrated to Americas 10,000 to 20,000 years ago. The comprehensive review showed that for each of these SCAs corresponded at least the development of one study group with a large production of scientific evidence often generalizable to all carriers of these conditions. Clusters of SCA populations in the American continents and the Caribbean provide unusual opportunity to gain insights into clinical and genetic characteristics of these disorders. Furthermore, the presence of large populations of patients living close to study centers can favor the development of meaningful clinical trials, which will impact on therapies and on quality of life of SCA carriers worldwide.

The Role of Noncoding Variants in Heritable Disease

The genetic basis of disease has largely focused on coding regions. However, it has become clear that a large proportion of the noncoding genome is functional and harbors genetic variants that contribute to disease etiology. Here, we review recent examples of inherited noncoding alterations that are responsible for Mendelian disorders or act to influence complex traits. We explore both rare and common genetic variants and discuss the wide range of mechanisms by which they affect gene regulation to promote disease. We also debate the challenges and progress associated with identifying and interpreting the functional and clinical significance of genetic variation in the context of the noncoding regulatory landscape.

TandemTools: mapping long reads and assessing/improving assembly quality in extra-long tandem repeats

MOTIVATION

Extra-long tandem repeats (ETRs) are widespread in eukaryotic genomes and play an important role in fundamental cellular processes, such as chromosome segregation. Although emerging long-read technologies have enabled ETR assemblies, the accuracy of such assemblies is difficult to evaluate since there are no tools for their quality assessment. Moreover, since the mapping of error-prone reads to ETRs remains an open problem, it is not clear how to polish draft ETR assemblies.

RESULTS

To address these problems, we developed the TandemTools software that includes the TandemMapper tool for mapping reads to ETRs and the TandemQUAST tool for polishing ETR assemblies and their quality assessment. We demonstrate that TandemTools not only reveals errors in ETR assemblies but also improves the recently generated assemblies of human centromeres.

AVAILABILITY AND IMPLEMENTATION

https://github.com/ablab/TandemTools.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

How have advances in genetic technology modified movement disorder nosology?

The role of genetics and its technological development have been fundamental in advancing the field of movement disorders, opening the door to precision medicine. Starting from the revolutionary discovery of the locus of the Huntington's disease gene, we review the milestones of genetic discoveries in movement disorders and their impact on clinical practice and research efforts. Before the 1980s, early techniques did not allow the identification of genetic alteration in complex diseases. Further advances increasingly defined a large number of pathogenic genetic alterations. Moreover, these techniques allowed epigenomic, transcriptomic and microbiome analyses. In the 2020s, these new technologies are poised to displace phenotype-based classifications towards a nosology based on genetic/biological data. Advances in genetic technologies are engineering a reversal of the phenotype-to-genotype order of nosology development, replacing convergent clinicopathological disease models with the genotypic divergence required for future precision medicine applications.

Exome sequencing in genetic disease: recent advances and considerations

Over the past decade, exome sequencing (ES) has allowed significant advancements to the field of disease research. By targeting the protein-coding regions of the genome, ES combines the depth of knowledge on protein-altering variants with high-throughput data generation and ease of analysis. New discoveries continue to be made using ES, and medical science has benefitted both theoretically and clinically from its continued use. In this review, we describe recent advances and successes of ES in disease research. Through selected examples of recent publications, we explore how ES continues to be a valuable tool to find variants that might explain disease etiology or provide insight into the biology underlying the disease. We then discuss shortcomings of ES in terms of variant discoveries made by other sequencing technologies that would be missed because of the scope and techniques of ES. We conclude with a brief outlook on the future of ES, suggesting that although newer and more thorough sequencing methods will soon supplant ES, its results will continue to be useful for disease research.

Pulse-Field capillary electrophoresis of repeat-primed PCR amplicons for analysis of large repeats in Spinocerebellar Ataxia Type 10

Large expansions of microsatellite DNA cause several neurological diseases. In Spinocerebellar ataxia type 10 (SCA10), the repeat interruptions change disease phenotype; an (ATTCC)n or a (ATCCT)n/(ATCCC)n interruption within the (ATTCT)n repeat is associated with the robust phenotype of ataxia and epilepsy while mostly pure (ATTCT)n may have reduced penetrance. Large repeat expansions of SCA10, and many other microsatellite expansions, can exceed 10,000 base pairs (bp) in size. Conventional next generation sequencing (NGS) technologies are ineffective in determining internal sequence contents or size of these expanded repeats. Using repeat primed PCR (RP-PCR) in conjunction with a high-sensitivity pulsed-field capillary electrophoresis fragment analyzer (FEMTO-Pulse, Agilent, Santa Clara, CA) (RP-FEMTO hereafter), we successfully determined sequence content of large expansion repeats in genomic DNA of SCA10 patients and transformed yeast artificial chromosomes containing SCA10 repeats. This RP-FEMTO is a simple and economical methodology which could complement emerging NGS for very long sequence reads such as Single Molecule, Real-Time (SMRT) and nanopore sequencing technologies.

Critical assessment of bioinformatics methods for the characterization of pathological repeat expansions with single-molecule sequencing data

A number of studies have reported the successful application of single-molecule sequencing technologies to the determination of the size and sequence of pathological expanded microsatellite repeats over the last 5 years. However, different custom bioinformatics pipelines were employed in each study, preventing meaningful comparisons and somewhat limiting the reproducibility of the results. In this review, we provide a brief summary of state-of-the-art methods for the characterization of expanded repeats alleles, along with a detailed comparison of bioinformatics tools for the determination of repeat length and sequence, using both real and simulated data. Our reanalysis of publicly available human genome sequencing data suggests a modest, but statistically significant, increase of the error rate of single-molecule sequencing technologies at genomic regions containing short tandem repeats. However, we observe that all the methods herein tested, irrespective of the strategy used for the analysis of the data (either based on the alignment or assembly of the reads), show high levels of sensitivity in both the detection of expanded tandem repeats and the estimation of the expansion size, suggesting that approaches based on single-molecule sequencing technologies are highly effective for the detection and quantification of tandem repeat expansions and contractions.

Long-read sequencing for rare human genetic diseases

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

Spinocerebellar ataxias: prospects and challenges for therapy development

The spinocerebellar ataxias (SCAs) comprise more than 40 autosomal dominant neurodegenerative disorders that present principally with progressive ataxia. Within the past few years, studies of pathogenic mechanisms in the SCAs have led to the development of promising therapeutic strategies, especially for SCAs caused by polyglutamine-coding CAG repeats. Nucleotide-based gene-silencing approaches that target the first steps in the pathogenic cascade are one promising approach not only for polyglutamine SCAs but also for the many other SCAs caused by toxic mutant proteins or RNA. For these and other emerging therapeutic strategies, well-coordinated preparation is needed for fruitful clinical trials. To accomplish this goal, investigators from the United States and Europe are now collaborating to share data from their respective SCA cohorts. Increased knowledge of the natural history of SCAs, including of the premanifest and early symptomatic stages of disease, will improve the prospects for success in clinical trials of disease-modifying drugs. In addition, investigators are seeking validated clinical outcome measures that demonstrate responsiveness to changes in SCA populations. Findings suggest that MRI and magnetic resonance spectroscopy biomarkers will provide objective biological readouts of disease activity and progression, but more work is needed to establish disease-specific biomarkers that track target engagement in therapeutic trials. Together, these efforts suggest that the development of successful therapies for one or more SCAs is not far away.

Single molecule real-time (SMRT) sequencing comes of age: applications and utilities for medical diagnostics

Short read massive parallel sequencing has emerged as a standard diagnostic tool in the medical setting. However, short read technologies have inherent limitations such as GC bias, difficulties mapping to repetitive elements, trouble discriminating paralogous sequences, and difficulties in phasing alleles. Long read single molecule sequencers resolve these obstacles. Moreover, they offer higher consensus accuracies and can detect epigenetic modifications from native DNA. The first commercially available long read single molecule platform was the RS system based on PacBio's single molecule real-time (SMRT) sequencing technology, which has since evolved into their RSII and Sequel systems. Here we capsulize how SMRT sequencing is revolutionizing constitutional, reproductive, cancer, microbial and viral genetic testing.

Long-Read Sequencing Emerging in Medical Genetics

The wide implementation of next-generation sequencing (NGS) technologies has revolutionized the field of medical genetics. However, the short read lengths of currently used sequencing approaches pose a limitation for the identification of structural variants, sequencing repetitive regions, phasing of alleles and distinguishing highly homologous genomic regions. These limitations may significantly contribute to the diagnostic gap in patients with genetic disorders who have undergone standard NGS, like whole exome or even genome sequencing. Now, the emerging long-read sequencing (LRS) technologies may offer improvements in the characterization of genetic variation and regions that are difficult to assess with the prevailing NGS approaches. LRS has so far mainly been used to investigate genetic disorders with previously known or strongly suspected disease loci. While these targeted approaches already show the potential of LRS, it remains to be seen whether LRS technologies can soon enable true whole genome sequencing routinely. Ultimately, this could allow the de novo assembly of individual whole genomes used as a generic test for genetic disorders. In this article, we summarize the current LRS-based research on human genetic disorders and discuss the potential of these technologies to facilitate the next major advancements in medical genetics.

Tandem-genotypes: robust detection of tandem repeat expansions from long DNA reads

Tandemly repeated DNA is highly mutable and causes at least 31 diseases, but it is hard to detect pathogenic repeat expansions genome-wide. Here, we report robust detection of human repeat expansions from careful alignments of long but error-prone (PacBio and nanopore) reads to a reference genome. Our method is robust to systematic sequencing errors, inexact repeats with fuzzy boundaries, and low sequencing coverage. By comparing to healthy controls, we prioritize pathogenic expansions within the top 10 out of 700,000 tandem repeats in whole genome sequencing data. This may help to elucidate the many genetic diseases whose causes remain unknown.

Mapping the landscape of tandem repeat variability by targeted long read single molecule sequencing in familial X-linked intellectual disability

BACKGROUND

The etiology of more than half of all patients with X-linked intellectual disability remains elusive, despite array-based comparative genomic hybridization, whole exome or genome sequencing. Since short read massive parallel sequencing approaches do not allow the detection of larger tandem repeat expansions, we hypothesized that such expansions could be a hidden cause of X-linked intellectual disability.

METHODS

We selectively captured over 1800 tandem repeats on the X chromosome and characterized them by long read single molecule sequencing in 3 families with idiopathic X-linked intellectual disability.

RESULTS

In male DNA samples, full tandem repeat length sequences were obtained for 88-93% of the targets and up to 99.6% of the repeats with a moderate guanine-cytosine content. Read length and analysis pipeline allow to detect cases of > 900 bp tandem repeat expansion. In one family, one repeat expansion co-occurs with down-regulation of the neighboring MIR222 gene. This gene has previously been implicated in intellectual disability and is apparently linked to FMR1 and NEFH overexpression associated with neurological disorders.

CONCLUSIONS

This study demonstrates the power of single molecule sequencing to measure tandem repeat lengths and detect expansions, and suggests that tandem repeat mutations may be a hidden cause of X-linked intellectual disability.

Efficient CRISPR/Cas9-mediated editing of trinucleotide repeat expansion in myotonic dystrophy patient-derived iPS and myogenic cells

CRISPR/Cas9 is an attractive platform to potentially correct dominant genetic diseases by gene editing with unprecedented precision. In the current proof-of-principle study, we explored the use of CRISPR/Cas9 for gene-editing in myotonic dystrophy type-1 (DM1), an autosomal-dominant muscle disorder, by excising the CTG-repeat expansion in the 3'-untranslated-region (UTR) of the human myotonic dystrophy protein kinase (DMPK) gene in DM1 patient-specific induced pluripotent stem cells (DM1-iPSC), DM1-iPSC-derived myogenic cells and DM1 patient-specific myoblasts. To eliminate the pathogenic gain-of-function mutant DMPK transcript, we designed a dual guide RNA based strategy that excises the CTG-repeat expansion with high efficiency, as confirmed by Southern blot and single molecule real-time (SMRT) sequencing. Correction efficiencies up to 90% could be attained in DM1-iPSC as confirmed at the clonal level, following ribonucleoprotein (RNP) transfection of CRISPR/Cas9 components without the need for selective enrichment. Expanded CTG repeat excision resulted in the disappearance of ribonuclear foci, a quintessential cellular phenotype of DM1, in the corrected DM1-iPSC, DM1-iPSC-derived myogenic cells and DM1 myoblasts. Consequently, the normal intracellular localization of the muscleblind-like splicing regulator 1 (MBNL1) was restored, resulting in the normalization of splicing pattern of SERCA1. This study validates the use of CRISPR/Cas9 for gene editing of repeat expansions.

A Census of Tandemly Repeated Polymorphic Loci in Genic Regions Through the Comparative Integration of Human Genome Assemblies

Polymorphic Tandem Repeat (PTR) is a common form of polymorphism in the human genome. A PTR consists in a variation found in an individual (or in a population) of the number of repeating units of a Tandem Repeat (TR) locus of the genome with respect to the reference genome. Several phenotypic traits and diseases have been discovered to be strongly associated with or caused by specific PTR loci. PTR are further distinguished in two main classes: Short Tandem Repeats (STR) when the repeating unit has size up to 6 base pairs, and Variable Number Tandem Repeats (VNTR) for repeating units of size above 6 base pairs. As larger and larger populations are screened via high throughput sequencing projects, it becomes technically feasible and desirable to explore the association between PTR and a panoply of such traits and conditions. In order to facilitate these studies, we have devised a method for compiling catalogs of PTR from assembled genomes, and we have produced a catalog of PTR for genic regions (exons, introns, UTR and adjacent regions) of the human genome (GRCh38). We applied four different TR discovery software tools to uncover in the first phase 55,223,485 TR (after duplicate removal) in GRCh38, of which 373,173 were determined to be PTR in the second phase by comparison with five assembled human genomes. Of these, 263,266 are not included by state-of-the-art PTR catalogs. The new methodology is mainly based on a hierarchical and systematic application of alignment-based sequence comparisons to identify and measure the polymorphism of TR. While previous catalogs focus on the class of STR of small total size, we remove any size restrictions, aiming at the more general class of PTR, and we also target fuzzy TR by using specific detection tools. Similarly to other previous catalogs of human polymorphic loci, we focus our catalog toward applications in the discovery of disease-associated loci. Validation by cross-referencing with existing catalogs on common clinically-relevant loci shows good concordance. Overall, this proposed census of human PTR in genic regions is a shared resource (web accessible), complementary to existing catalogs, facilitating future genome-wide studies involving PTR.

Expression of Pyridoxal 5'-Phosphate-Independent Racemases Can Reduce 2-Aminoacrylate Stress in Salmonella enterica

The RidA protein (PF01042) from Salmonella enterica is a deaminase that quenches 2-aminoacrylate (2AA) and other reactive metabolites. In the absence of RidA, 2AA accumulates, damages cellular enzymes, and compromises the metabolic network. In vitro, RidA homologs from all domains of life deaminate 2AA, and RidA proteins from plants, bacteria, yeast, and humans complement the mutant phenotype of a ridA mutant strain of S. enterica In the present study, a methanogenic archaeon, Methanococcus maripaludis S2, was used to probe alternative mechanisms to restore metabolic balance. M. maripaludis MMP0739, which is annotated as an aspartate/glutamate racemase, complemented a ridA mutant strain and reduced the intracellular 2AA burden. The aspartate/glutamate racemase YgeA from Escherichia coli or S. enterica, when provided in trans, similarly restored wild-type growth to a ridA mutant. These results uncovered a new mechanism to ameliorate metabolic stress, and they suggest that direct quenching by RidA is not the only strategy to quench 2AA.IMPORTANCE 2-Aminoacrylate is an endogenously generated reactive metabolite that can damage cellular enzymes if not directly quenched by the conserved deaminase RidA. This study used an archaeon to identify a RidA-independent mechanism to prevent metabolic stress caused by 2AA. The data suggest that a gene product annotated as an aspartate/glutamate racemase (MMP0739) produces a metabolite that can quench 2AA, expanding our understanding of strategies available to quench reactive metabolites.

Parkinson's disease associated with pure ATXN10 repeat expansion

Large, non-coding pentanucleotide repeat expansions of ATTCT in intron 9 of the ATXN10 gene typically cause progressive spinocerebellar ataxia with or without seizures and present neuropathologically with Purkinje cell loss resulting in symmetrical cerebellar atrophy. These ATXN10 repeat expansions can be interrupted by sequence motifs which have been attributed to seizures and are likely to act as genetic modifiers. We identified a Mexican kindred with multiple affected family members with ATXN10 expansions. Four affected family members showed clinical features of spinocerebellar ataxia type 10 (SCA10). However, one affected individual presented with early-onset levodopa-responsive parkinsonism, and one family member carried a large repeat ATXN10 expansion, but was clinically unaffected. To characterize the ATXN10 repeat, we used a novel technology of single-molecule real-time (SMRT) sequencing and CRISPR/Cas9-based capture. We sequenced the entire span of ~5.3-7.0 kb repeat expansions. The Parkinson's patient carried an ATXN10 expansion with no repeat interruption motifs as well as an unaffected sister. In the siblings with typical SCA10, we found a repeat pattern of ATTCC repeat motifs that have not been associated with seizures previously. Our data suggest that the absence of repeat interruptions is likely a genetic modifier for the clinical presentation of l-Dopa responsive parkinsonism, whereas repeat interruption motifs contribute clinically to epilepsy. Repeat interruptions are important genetic modifiers of the clinical phenotype in SCA10. Advanced sequencing techniques now allow to better characterize the underlying genetic architecture for determining accurate phenotype-genotype correlations.

Interrogating the "unsequenceable" genomic trinucleotide repeat disorders by long-read sequencing

Microsatellite expansion, such as trinucleotide repeat expansion (TRE), is known to cause a number of genetic diseases. Sanger sequencing and next-generation short-read sequencing are unable to interrogate TRE reliably. We developed a novel algorithm called RepeatHMM to estimate repeat counts from long-read sequencing data. Evaluation on simulation data, real amplicon sequencing data on two repeat expansion disorders, and whole-genome sequencing data generated by PacBio and Oxford Nanopore technologies showed superior performance over competing approaches. We concluded that long-read sequencing coupled with RepeatHMM can estimate repeat counts on microsatellites and can interrogate the “unsequenceable” genomic trinucleotide repeat disorders.

Inheritance patterns of ATCCT repeat interruptions in spinocerebellar ataxia type 10 (SCA10) expansions

Spinocerebellar ataxia type 10 (SCA10), an autosomal dominant cerebellar ataxia disorder, is caused by a non-coding ATTCT microsatellite repeat expansion in the ataxin 10 gene. In a subset of SCA10 families, the 5’-end of the repeat expansion contains a complex sequence of penta- and heptanucleotide interruption motifs which is followed by a pure tract of tandem ATCCT repeats of unknown length at its 3’-end. Intriguingly, expansions that carry these interruption motifs correlate with an epileptic seizure phenotype and are unstable despite the theory that interruptions are expected to stabilize expanded repeats. To examine the apparent contradiction of unstable, interruption-positive SCA10 expansion alleles and to determine whether the instability originates outside of the interrupted region, we sequenced approximately 1 kb of the 5’-end of SCA10 expansions using the ATCCT-PCR product in individuals across multiple generations from four SCA10 families. We found that the greatest instability within this region occurred in paternal transmissions of the allele in stretches of pure ATTCT motifs while the intervening interrupted sequences were stable. Overall, the ATCCT interruption changes by only one to three repeat units and therefore cannot account for the instability across the length of the disease allele. We conclude that the AT-rich interruptions locally stabilize the SCA10 expansion at the 5’-end but do not completely abolish instability across the entire span of the expansion. In addition, analysis of the interruption alleles across these families support a parsimonious single origin of the mutation with a shared distant ancestor.

RNA toxicity and foci formation in microsatellite expansion diseases

More than 30 incurable neurological and neuromuscular diseases are caused by simple microsatellite expansions consisted of 3-6 nucleotides. These repeats can occur in non-coding regions and often result in a dominantly inherited disease phenotype that is characteristic of a toxic RNA gain-of-function. The expanded RNA adopts unusual secondary structures, sequesters various RNA binding proteins to form insoluble nuclear foci, and causes cellular defects at a multisystem level. Nuclear foci are dynamic in size, shape and colocalization of RNA binding proteins in different expansion diseases and tissue types. This review sets to provide new insights into the disease mechanisms of RNA toxicity and foci modulation, in light of recent advancement on bi-directional transcription, antisense RNA, repeat-associated non-ATG translation and beyond.

DNA repair in the trinucleotide repeat disorders

BACKGROUND

Inherited diseases caused by unstable repeated DNA sequences are rare, but together represent a substantial cause of morbidity. Trinucleotide repeat disorders are severe, usually life-shortening, neurological disorders caused by nucleotide expansions, and most have no disease-modifying treatments. Longer repeat expansions are associated with genetic anticipation (ie, earlier disease onset in successive generations), although the differences in age at onset are not entirely accounted for by repeat length. Such phenotypic variation within disorders implies the existence of additional modifying factors in pathways that can potentially be modulated to treat disease.

RECENT DEVELOPMENTS

A genome-wide association study detected genetic modifiers of age at onset in Huntington's disease. Similar findings were seen in the spinocerebellar ataxias, indicating an association between DNA damage-response and repair pathways and the age at onset of disease. These studies also suggest that a common genetic mechanism modulates age at onset across polyglutamine diseases and could extend to other repeat expansion disorders. Genetic defects in DNA repair underlie other neurodegenerative disorders (eg, ataxia-telangiectasia), and DNA double-strand breaks are crucial to the modulation of early gene expression, which provides a mechanistic link between DNA repair and neurodegeneration. Mismatch and base-excision repair are important in the somatic expansion of repeated sequences in mouse models of trinucleotide repeat disorders, and somatic expansion of the expanded CAG tract in HTT correlates with age at onset of Huntington's disease and other trinucleotide repeat disorders. WHERE NEXT?: To understand the common genetic architecture of trinucleotide repeat disorders and any further genetic susceptibilities in individual disorders, genetic analysis with increased numbers of variants and sample sizes is needed, followed by sequencing approaches to define the phenotype-modifying variants. The findings must then be translated into cell biology analyses to elucidate the mechanisms through which the genetic variants operate. Genes that have roles in the DNA damage response could underpin a common DNA repeat-based mechanism and provide new therapeutic targets (and hence therapeutics) in multiple trinucleotide repeat disorders.

Genomic analyses of multidrug resistant Pseudomonas aeruginosa PA1 resequenced by single-molecule real-time sequencing

As a third-generation sequencing (TGS) method, single-molecule real-time (SMRT) technology provides long read length, and it is well suited for resequencing projects and de novo assembly. In the present study, Pseudomonas aeruginosa PA1 was characterized and resequenced using SMRT technology. PA1 was also subjected to genomic, comparative and pan-genomic analyses. The multidrug resistant strain PA1 possesses a 6,498,072 bp genome and a sequence type of ST-782. The genome of PA1 was also visualized, and the results revealed the details of general genome annotations, virulence factors, regulatory proteins (RPs), secretion system proteins, type II toxin–antitoxin (T–A) pairs and genomic islands. Whole genome comparison analysis suggested that PA1 exhibits similarity to other P. aeruginosa strains but differs in terms of horizontal gene transfer (HGT) regions, such as prophages and genomic islands. Phylogenetic analyses based on 16S rRNA sequences demonstrated that PA1 is closely related to PAO1, and P. aeruginosa strains can be divided into two main groups. The pan-genome of P. aeruginosa consists of a core genome of approximately 4,000 genes and an accessory genome of at least 6,600 genes. The present study presented a detailed, visualized and comparative analysis of the PA1 genome, to enhance our understanding of this notorious pathogen.