Rapid detection of expanded short tandem repeats in personal genomics using hybrid sequencing

Aicardi Syndrome Is a Genetically Heterogeneous Disorder

Aicardi Syndrome (AIC) is a rare neurodevelopmental disorder recognized by the classical triad of agenesis of the corpus callosum, chorioretinal lacunae and infantile epileptic spasms syndrome. The diagnostic criteria of AIC were revised in 2005 to include additional phenotypes that are frequently observed in this patient group. AIC has been traditionally considered as X-linked and male lethal because it almost exclusively affects females. Despite numerous genetic and genomic investigations on AIC, a unifying X-linked cause has not been identified. Here, we performed exome and genome sequencing of 10 females with AIC or suspected AIC based on current criteria. We identified a unique de novo variant, each in different genes: KMT2B, SLF1, SMARCB1, SZT2 and WNT8B, in five of these females. Notably, genomic analyses of coding and non-coding single nucleotide variants, short tandem repeats and structural variation highlighted a distinct lack of X-linked candidate genes. We assessed the likely pathogenicity of our candidate autosomal variants using the TOPflash assay for WNT8B and morpholino knockdown in zebrafish (Danio rerio) embryos for other candidates. We show expression of Wnt8b and Slf1 are restricted to clinically relevant cortical tissues during mouse development. Our findings suggest that AIC is genetically heterogeneous with implicated genes converging on molecular pathways central to cortical development.

Exploring microbial functional biodiversity at the protein family level-From metagenomic sequence reads to annotated protein clusters

Metagenomics has enabled accessing the genetic repertoire of natural microbial communities. Metagenome shotgun sequencing has become the method of choice for studying and classifying microorganisms from various environments. To this end, several methods have been developed to process and analyze the sequence data from raw reads to end-products such as predicted protein sequences or families. In this article, we provide a thorough review to simplify such processes and discuss the alternative methodologies that can be followed in order to explore biodiversity at the protein family level. We provide details for analysis tools and we comment on their scalability as well as their advantages and disadvantages. Finally, we report the available data repositories and recommend various approaches for protein family annotation related to phylogenetic distribution, structure prediction and metadata enrichment.

Detection of repeat expansions in large next generation DNA and RNA sequencing data without alignment

Bioinformatic methods for detecting short tandem repeat expansions in short-read sequencing have identified new repeat expansions in humans, but require alignment information to identify repetitive motif enrichment at genomic locations. We present superSTR, an ultrafast method that does not require alignment. superSTR is used to process whole-genome and whole-exome sequencing data, and perform the first STR analysis of the UK Biobank, efficiently screening and identifying known and potential disease-associated STRs in the exomes of 49,953 biobank participants. We demonstrate the first bioinformatic screening of RNA sequencing data to detect repeat expansions in humans and mouse models of ataxia and dystrophy.

GC-rich repeat expansions: associated disorders and mechanisms

Noncoding repeat expansions are a well-known cause of genetic disorders mainly affecting the central nervous system. Missed by most standard technologies used in routine diagnosis, pathogenic noncoding repeat expansions have to be searched for using specific techniques such as repeat-primed PCR or specific bioinformatics tools applied to genome data, such as ExpansionHunter. In this review, we focus on GC-rich repeat expansions, which represent at least one third of all noncoding repeat expansions described so far. GC-rich expansions are mainly located in regulatory regions (promoter, 5′ untranslated region, first intron) of genes and can lead to either a toxic gain-of-function mediated by RNA toxicity and/or repeat-associated non-AUG (RAN) translation, or a loss-of-function of the associated gene, depending on their size and their methylation status. We herein review the clinical and molecular characteristics of disorders associated with these difficult-to-detect expansions.

Neurodegenerative diseases associated with non-coding CGG tandem repeat expansions

Non-coding CGG repeat expansions cause multiple neurodegenerative disorders, including fragile X-associated tremor/ataxia syndrome, neuronal intranuclear inclusion disease, oculopharyngeal myopathy with leukodystrophy, and oculopharyngodistal myopathy. The underlying genetic causes of several of these diseases have been identified only in the past 2-3 years. These expansion disorders have substantial overlapping clinical, neuroimaging and histopathological features. The shared features suggest common mechanisms that could have implications for the development of therapies for this group of diseases - similar therapeutic strategies or drugs may be effective for various neurodegenerative disorders induced by non-coding CGG expansions. In this Review, we provide an overview of clinical and pathological features of these CGG repeat expansion diseases and consider the likely pathological mechanisms, including RNA toxicity, CGG repeat-associated non-AUG-initiated translation, protein aggregation and mitochondrial impairment. We then discuss future research needed to improve the identification and diagnosis of CGG repeat expansion diseases, to improve modelling of these diseases and to understand their pathogenesis. We also consider possible therapeutic strategies. Finally, we propose that CGG repeat expansion diseases may represent manifestations of a single underlying neuromyodegenerative syndrome in which different organs are affected to different extents depending on the gene location of the repeat expansion.

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.

Advancing genomic technologies and clinical awareness accelerates discovery of disease-associated tandem repeat sequences

Expansions of gene-specific DNA tandem repeats (TRs), first described in 1991 as a disease-causing mutation in humans, are now known to cause >60 phenotypes, not just disease, and not only in humans. TRs are a common form of genetic variation with biological consequences, observed, so far, in humans, dogs, plants, oysters, and yeast. Repeat diseases show atypical clinical features, genetic anticipation, and multiple and partially penetrant phenotypes among family members. Discovery of disease-causing repeat expansion loci accelerated through technological advances in DNA sequencing and computational analyses. Between 2019 and 2021, 17 new disease-causing TR expansions were reported, totaling 63 TR loci (>69 diseases), with a likelihood of more discoveries, and in more organisms. Recent and historical lessons reveal that properly assessed clinical presentations, coupled with genetic and biological awareness, can guide discovery of disease-causing unstable TRs. We highlight critical but underrecognized aspects of TR mutations. Repeat motifs may not be present in current reference genomes but will be in forthcoming gapless long-read references. Repeat motif size can be a single nucleotide to kilobases/unit. At a given locus, repeat motif sequence purity can vary with consequence. Pathogenic repeats can be "insertions" within nonpathogenic TRs. Expansions, contractions, and somatic length variations of TRs can have clinical/biological consequences. TR instabilities occur in humans and other organisms. TRs can be epigenetically modified and/or chromosomal fragile sites. We discuss the expanding field of disease-associated TR instabilities, highlighting prospects, clinical and genetic clues, tools, and challenges for further discoveries of disease-causing TR instabilities and understanding their biological and pathological impacts-a vista that is about to expand.

Genome sequencing of 320 Chinese children with epilepsy: a clinical and molecular study

The aim of this study is to evaluate the diagnostic value of genome sequencing in children with epilepsy, and to provide genome sequencing-based insights into the molecular genetic mechanisms of epilepsy to help establish accurate diagnoses, design appropriate treatments and assist in genetic counselling. We performed genome sequencing on 320 Chinese children with epilepsy, and interpreted single-nucleotide variants and copy number variants of all samples. The complete pedigree and clinical data of the probands were established and followed up. The clinical phenotypes, treatments, prognoses and genotypes of the patients were analysed. Age at seizure onset ranged from 1 day to 17 years, with a median of 4.3 years. Pathogenic/likely pathogenic variants were found in 117 of the 320 children (36.6%), of whom 93 (29.1%) had single-nucleotide variants, 22 (6.9%) had copy number variants and two had both single-nucleotide variants and copy number variants. Single-nucleotide variants were most frequently found in SCN1A (10/95, 10.5%), which is associated with Dravet syndrome, followed by PRRT2 (8/95, 8.4%), which is associated with benign familial infantile epilepsy, and TSC2 (7/95, 7.4%), which is associated with tuberous sclerosis. Among the copy number variants, there were three with a length <25 kilobases. The most common recurrent copy number variants were 17p13.3 deletions (5/24, 20.8%), 16p11.2 deletions (4/24, 16.7%), and 7q11.23 duplications (2/24, 8.3%), which are associated with epilepsy, developmental retardation and congenital abnormalities. Four particular 16p11.2 deletions and two 15q11.2 deletions were considered to be susceptibility factors contributing to neurodevelopmental disorders associated with epilepsy. The diagnostic yield was 75.0% in patients with seizure onset during the first postnatal month, and gradually decreased in patients with seizure onset at a later age. Forty-two patients (13.1%) were found to be specifically treatable for the underlying genetic cause identified by genome sequencing. Three of them received corresponding targeted therapies and demonstrated favourable prognoses. Genome sequencing provides complete genetic diagnosis, thus enabling individualized treatment and genetic counselling for the parents of the patients. Genome sequencing is expected to become the first choice of methods for genetic testing of patients with epilepsy.

Finding long tandem repeats in long noisy reads

Motivation

Long tandem repeat expansions of more than 1000 nt have been suggested to be associated with diseases, but remain largely unexplored in individual human genomes because read lengths have been too short. However, new long-read sequencing technologies can produce single reads of 10 000 nt or more that can span such repeat expansions, although these long reads have high error rates, of 10–20%, which complicates the detection of repetitive elements. Moreover, most traditional algorithms for finding tandem repeats are designed to find short tandem repeats (<1000 nt) and cannot effectively handle the high error rate of long reads in a reasonable amount of time.

Results

Here, we report an efficient algorithm for solving this problem that takes advantage of the length of the repeat. Namely, a long tandem repeat has hundreds or thousands of approximate copies of the repeated unit, so despite the error rate, many short k-mers will be error-free in many copies of the unit. We exploited this characteristic to develop a method for first estimating regions that could contain a tandem repeat, by analyzing the k-mer frequency distributions of fixed-size windows across the target read, followed by an algorithm that assembles the k-mers of a putative region into the consensus repeat unit by greedily traversing a de Bruijn graph. Experimental results indicated that the proposed algorithm largely outperformed Tandem Repeats Finder, a widely used program for finding tandem repeats, in terms of sensitivity.

Availability and implementation

https://github.com/morisUtokyo/mTR.

30 years of repeat expansion disorders: What have we learned and what are the remaining challenges?

Tandem repeats represent one of the most abundant class of variations in human genomes, which are polymorphic by nature and become highly unstable in a length-dependent manner. The expansion of repeat length across generations is a well-established process that results in human disorders mainly affecting the central nervous system. At least 50 disorders associated with expansion loci have been described to date, with half recognized only in the last ten years, as prior methodological difficulties limited their identification. These limitations still apply to the current widely used molecular diagnostic methods (exome or gene panels) and thus result in missed diagnosis detrimental to affected individuals and their families, especially for disorders that are very rare and/or clinically not recognizable. Most of these disorders have been identified through family-driven approaches and many others likely remain to be identified. The recent development of long-read technologies provides a unique opportunity to systematically investigate the contribution of tandem repeats and repeat expansions to the genetic architecture of human disorders. In this review, we summarize the current and most recent knowledge about the genetics of repeat expansion disorders and the diversity of their pathophysiological mechanisms and outline the perspectives of developing personalized treatments in the future.

Different types of disease-causing noncoding variants revealed by genomic and gene expression analyses in families with X-linked intellectual disability

The pioneering discovery research of X-linked intellectual disability (XLID) genes has benefitted thousands of individuals worldwide; however, approximately 30% of XLID families still remain unresolved. We postulated that noncoding variants that affect gene regulation or splicing may account for the lack of a genetic diagnosis in some cases. Detecting pathogenic, gene-regulatory variants with the same sensitivity and specificity as structural and coding variants is a major challenge for Mendelian disorders. Here, we describe three pedigrees with suggestive XLID where distinctive phenotypes associated with known genes guided the identification of three different noncoding variants. We used comprehensive structural, single-nucleotide, and repeat expansion analyses of genome sequencing. RNA-Seq from patient-derived cell lines, reverse-transcription polymerase chain reactions, Western blots, and reporter gene assays were used to confirm the functional effect of three fundamentally different classes of pathogenic noncoding variants: a retrotransposon insertion, a novel intronic splice donor, and a canonical splice variant of an untranslated exon. In one family, we excluded a rare coding variant in ARX, a known XLID gene, in favor of a regulatory noncoding variant in OFD1 that correlated with the clinical phenotype. Our results underscore the value of genomic research on unresolved XLID families to aid novel, pathogenic noncoding variant discovery.

Y-LineageTracker: a high-throughput analysis framework for Y-chromosomal next-generation sequencing data

BACKGROUND

Y-chromosome DNA (Y-DNA) has been used for tracing paternal lineages and offers a clear path from an individual to a known, or likely, direct paternal ancestor. The advance of next-generation sequencing (NGS) technologies increasingly improves the resolution of the non-recombining region of the Y-chromosome (NRY). However, a lack of suitable computer tools prevents the use of NGS data from the Y-DNA studies.

RESULTS

We developed Y-LineageTracker, a high-throughput analysis framework that not only utilizes state-of-the-art methodologies to automatically determine NRY haplogroups and identify microsatellite variants of Y-chromosome on a fine scale, but also optimizes comprehensive Y-DNA analysis methods for NGS data. Notably, Y-LineageTracker integrates the NRY haplogroup and Y-STR analysis modules with recognized strategies to robustly suggest an interpretation for paternal genetics and evolution. NRY haplogroup module mainly covers haplogroup classification, clustering analysis, phylogeny construction, and divergence time estimation of NRY haplogroups, and Y-STR module mainly includes Y-STR genotyping, statistical calculation, network analysis, and estimation of time to the most recent common ancestor (TMRCA) based on Y-STR haplotypes. Performance comparison indicated that Y-LineageTracker outperformed existing Y-DNA analysis tools for the high performance and satisfactory visualization effect.

CONCLUSIONS

Y-LineageTracker is an open-source and user-friendly command-line tool that provide multiple functions to efficiently analyze Y-DNA from NGS data at both Y-SNP and Y-STR level. Additionally, Y-LineageTracker supports various formats of input data and produces high-quality figures suitable for publication. Y-LineageTracker is coded with Python3 and supports Windows, Linux, and macOS platforms, and can be installed manually or via the Python Package Index (PyPI). The source code, examples, and manual of Y-LineageTracker are freely available at https://www.picb.ac.cn/PGG/resource.php or CodeOcean ( https://codeocean.com/capsule/7424381/tree ).

TRiCoLOR: tandem repeat profiling using whole-genome long-read sequencing data

BACKGROUND

Tandem repeat sequences are widespread in the human genome, and their expansions cause multiple repeat-mediated disorders. Genome-wide discovery approaches are needed to fully elucidate their roles in health and disease, but resolving tandem repeat variation accurately remains a challenging task. While traditional mapping-based approaches using short-read data have severe limitations in the size and type of tandem repeats they can resolve, recent third-generation sequencing technologies exhibit substantially higher sequencing error rates, which complicates repeat resolution.

RESULTS

We developed TRiCoLOR, a freely available tool for tandem repeat profiling using error-prone long reads from third-generation sequencing technologies. The method can identify repetitive regions in sequencing data without a prior knowledge of their motifs or locations and resolve repeat multiplicity and period size in a haplotype-specific manner. The tool includes methods to interactively visualize the identified repeats and to trace their Mendelian consistency in pedigrees.

CONCLUSIONS

TRiCoLOR demonstrates excellent performance and improved sensitivity and specificity compared with alternative tools on synthetic data. For real human whole-genome sequencing data, TRiCoLOR achieves high validation rates, suggesting its suitability to identify tandem repeat variation in personal genomes.

Advances in repeat expansion diseases and a new concept of repeat motif-phenotype correlation

Recently repeat expansions have been found in more than 10 diseases in the past two years. Because the same repeat motifs are found in similar disease (as exemplified by benign adult familial myoclonic epilepsy) or in diseases with overlapping phenotype (as exemplified by fragile X tremor/ataxia syndrome, neuronal intranuclear inclusion disease, oculopharyngeal myopathy with leukoencephalopathy, and oculopharyngodistal myopathy), we propose a new concept of 'repeat motif-phenotype correlation', which argue for toxic gain-of-function mechanism caused by expanded repeats, rather than altered functions of genes harboring expanded repeats. The concept is expected to help identify repeat expansions taking the similar or overlapping clinical presentations as the clues. Although repeat expansions have been identified predominantly in autosomal dominant diseases, recent progresses have demonstrated that they are also observed in autosomal recessive diseases. Furthermore, repeat expansions are not infrequently observed in patients without family histories, which urges us to pay attention to sporadic diseases. We should expand our views toward repeat expansion diseases to accelerate discovery of diseases caused by repeat expansions, better understanding the disease mechanisms, and development of therapeutic measures.

Unstable TTTTA/TTTCA expansions in MARCH6 are associated with Familial Adult Myoclonic Epilepsy type 3

Familial Adult Myoclonic Epilepsy (FAME) is a genetically heterogeneous disorder characterized by cortical tremor and seizures. Intronic TTTTA/TTTCA repeat expansions in SAMD12 (FAME1) are the main cause of FAME in Asia. Using genome sequencing and repeat-primed PCR, we identify another site of this repeat expansion, in MARCH6 (FAME3) in four European families. Analysis of single DNA molecules with nanopore sequencing and molecular combing show that expansions range from 3.3 to 14 kb on average. However, we observe considerable variability in expansion length and structure, supporting the existence of multiple expansion configurations in blood cells and fibroblasts of the same individual. Moreover, the largest expansions are associated with micro-rearrangements occurring near the expansion in 20% of cells. This study provides further evidence that FAME is caused by intronic TTTTA/TTTCA expansions in distinct genes and reveals that expansions exhibit an unexpectedly high somatic instability that can ultimately result in genomic rearrangements.

lordFAST: sensitive and Fast Alignment Search Tool for LOng noisy Read sequencing Data

Motivation

Recent advances in genomics and precision medicine have been made possible through the application of high throughput sequencing (HTS) to large collections of human genomes. Although HTS technologies have proven their use in cataloging human genome variation, computational analysis of the data they generate is still far from being perfect. The main limitation of Illumina and other popular sequencing technologies is their short read length relative to the lengths of (common) genomic repeats. Newer (single molecule sequencing - SMS) technologies such as Pacific Biosciences and Oxford Nanopore are producing longer reads, making it theoretically possible to overcome the difficulties imposed by repeat regions. Unfortunately, because of their high sequencing error rate, reads generated by these technologies are very difficult to work with and cannot be used in many of the standard downstream analysis pipelines. Note that it is not only difficult to find the correct mapping locations of such reads in a reference genome, but also to establish their correct alignment so as to differentiate sequencing errors from real genomic variants. Furthermore, especially since newer SMS instruments provide higher throughput, mapping and alignment need to be performed much faster than before, maintaining high sensitivity.

Results

We introduce lordFAST, a novel long-read mapper that is specifically designed to align reads generated by PacBio and potentially other SMS technologies to a reference. lordFAST not only has higher sensitivity than the available alternatives, it is also among the fastest and has a very low memory footprint.

Availability and implementation

lordFAST is implemented in C++ and supports multi-threading. The source code of lordFAST is available at https://github.com/vpc-ccg/lordfast.

Supplementary information

Supplementary data are available at Bioinformatics online.

Characterizing microsatellite polymorphisms using assembly-based and mapping-based tools

Microsatellite polymorphism has always been a challenge for genome assembly and sequence alignment due to sequencing errors, short read lengths, and high incidence of polymerase slippage in microsatellite regions. Despite the information they carry being very valuable, microsatellite variations have not gained enough attention to be a routine step in genome sequence analysis pipelines. After the completion of the 1000 Genomes Project, which aimed to establish the most detailed genetic variation catalog for humans, the consortium released only two microsatellite prediction sets generated by two tools. Many other large research efforts have failed to shed light on microsatellite variations. We evaluated the performance of three different local assembly methods on three different experimental settings, focusing on genotype-based performance, coverage impact, and preprocessing including flanking regions. All these experiments supported our initial expectations on assembly. We also demonstrate that overlap-layout-consensus (OLC)-basedassembly methods show higher sensitivity to microsatellite variant calling when compared to a de Bruijn graph-based approach. We conclude that assembly with OLC is the better method for genotyping microsatellites. Our pipeline is available at https://github.com/gulfemd/STRAssembly.

Iso-Seq Allows Genome-Independent Transcriptome Profiling of Grape Berry Development

Transcriptomics has been widely applied to study grape berry development. With few exceptions, transcriptomic studies in grape are performed using the available genome sequence, PN40024, as reference. However, differences in gene content among grape accessions, which contribute to phenotypic differences among cultivars, suggest that a single reference genome does not represent the species' entire gene space. Though whole genome assembly and annotation can reveal the relatively unique or "private" gene space of any particular cultivar, transcriptome reconstruction is a more rapid, less costly, and less computationally intensive strategy to accomplish the same goal. In this study, we used single molecule-real time sequencing (SMRT) to sequence full-length cDNA (Iso-Seq) and reconstruct the transcriptome of Cabernet Sauvignon berries during berry ripening. In addition, short reads from ripening berries were used to error-correct low-expression isoforms and to profile isoform expression. By comparing the annotated gene space of Cabernet Sauvignon to other grape cultivars, we demonstrate that the transcriptome reference built with Iso-Seq data represents most of the expressed genes in the grape berries and includes 1,501 cultivar-specific genes. Iso-Seq produced transcriptome profiles similar to those obtained after mapping on a complete genome reference. Together, these results justify the application of Iso-Seq to identify cultivar-specific genes and build a comprehensive reference for transcriptional profiling that circumvents the necessity of a genome reference with its associated costs and computational weight.

De novo repeat interruptions are associated with reduced somatic instability and mild or absent clinical features in myotonic dystrophy type 1

Myotonic dystrophy type 1 (DM1) is a multisystem disorder, caused by expansion of a CTG trinucleotide repeat in the 3'-untranslated region of the DMPK gene. The repeat expansion is somatically unstable and tends to increase in length with time, contributing to disease progression. In some individuals, the repeat array is interrupted by variant repeats such as CCG and CGG, stabilising the expansion and often leading to milder symptoms. We have characterised three families, each including one person with variant repeats that had arisen de novo on paternal transmission of the repeat expansion. Two individuals were identified for screening due to an unusual result in the laboratory diagnostic test, and the third due to exceptionally mild symptoms. The presence of variant repeats in all three expanded alleles was confirmed by restriction digestion of small pool PCR products, and allele structures were determined by PacBio sequencing. Each was different, but all contained CCG repeats close to the 3'-end of the repeat expansion. All other family members had inherited pure CTG repeats. The variant repeat-containing alleles were more stable in the blood than pure alleles of similar length, which may in part account for the mild symptoms observed in all three individuals. This emphasises the importance of somatic instability as a disease mechanism in DM1. Further, since patients with variant repeats may have unusually mild symptoms, identification of these individuals has important implications for genetic counselling and for patient stratification in DM1 clinical trials.

Polyglutamine Repeats in Viruses

This review explores the presence and functions of polyglutamine (polyQ) in viral proteins. In mammals, mutations in polyQ segments (and CAG repeats at the nucleotide level) have been linked to neural disorders and ataxias. PolyQ regions in normal human proteins have documented functional roles, in transcription factors and, more recently, in regulating autophagy. Despite the high frequency of polyQ repeats in eukaryotic genomes, little attention has been given to the presence or possible role of polyQ sequences in virus genomes. A survey described here revealed that polyQ repeats occur rarely in RNA viruses, suggesting that they have detrimental effects on virus replication at the nucleotide or protein level. However, there have been sporadic reports of polyQ segments in potyviruses and in reptilian nidoviruses (among the largest RNA viruses known). Conserved polyQ segments are found in the regulatory control proteins of many DNA viruses. Variable length polyQ tracts are found in proteins that contribute to transmissibility (cowpox A-type inclusion protein (ATI)) and control of latency (herpes viruses). New longer-read sequencing methods, using original biological samples, should reveal more details on the presence and functional role of polyQ in viruses, as well as the nucleotide regions that encode them. Given the known toxic effects of polyQ repeats, the role of these segments in neurovirulent and tumorigenic viruses should be further explored.

Expansions of intronic TTTCA and TTTTA repeats in benign adult familial myoclonic epilepsy

Epilepsy is a common neurological disorder, and mutations in genes encoding ion channels or neurotransmitter receptors are frequent causes of monogenic forms of epilepsy. Here we show that abnormal expansions of TTTCA and TTTTA repeats in intron 4 of SAMD12 cause benign adult familial myoclonic epilepsy (BAFME). Single-molecule, real-time sequencing of BAC clones and nanopore sequencing of genomic DNA identified two repeat configurations in SAMD12. Intriguingly, in two families with a clinical diagnosis of BAFME in which no repeat expansions in SAMD12 were observed, we identified similar expansions of TTTCA and TTTTA repeats in introns of TNRC6A and RAPGEF2, indicating that expansions of the same repeat motifs are involved in the pathogenesis of BAFME regardless of the genes in which the expanded repeats are located. This discovery that expansions of noncoding repeats lead to neuronal dysfunction responsible for myoclonic tremor and epilepsy extends the understanding of diseases with such repeat expansion.

Methods for Assessing DNA Repair and Repeat Expansion in Huntington's Disease

Huntington's disease (HD) is caused by a CAG repeat expansion in the HTT gene. Repeat length can change over time, both in individual cells and between generations, and longer repeats may drive pathology. Cellular DNA repair systems have long been implicated in CAG repeat instability but recent genetic evidence from humans linking DNA repair variants to HD onset and progression has reignited interest in this area. The DNA damage response plays an essential role in maintaining genome stability, but may also license repeat expansions in the context of HD. In this chapter we summarize the methods developed to assay CAG repeat expansion/contraction in vitro and in cells, and review the DNA repair genes tested in mouse models of HD. While none of these systems is currently ideal, new technologies, such as long-read DNA sequencing, should improve the sensitivity of assays to assess the effects of DNA repair pathways in HD. Improved assays will be essential precursors to high-throughput testing of small molecules that can alter specific steps in DNA repair pathways and perhaps ameliorate expansion or enhance contraction of the HTT CAG repeat.

Profiling of Short-Tandem-Repeat Disease Alleles in 12,632 Human Whole Genomes

Short tandem repeats (STRs) are hyper-mutable sequences in the human genome. They are often used in forensics and population genetics and are also the underlying cause of many genetic diseases. There are challenges associated with accurately determining the length polymorphism of STR loci in the genome by next-generation sequencing (NGS). In particular, accurate detection of pathological STR expansion is limited by the sequence read length during whole-genome analysis. We developed TREDPARSE, a software package that incorporates various cues from read alignment and paired-end distance distribution, as well as a sequence stutter model, in a probabilistic framework to infer repeat sizes for genetic loci, and we used this software to infer repeat sizes for 30 known disease loci. Using simulated data, we show that TREDPARSE outperforms other available software. We sampled the full genome sequences of 12,632 individuals to an average read depth of approximately 30× to 40× with Illumina HiSeq X. We identified 138 individuals with risk alleles at 15 STR disease loci. We validated a representative subset of the samples (n = 19) by Sanger and by Oxford Nanopore sequencing. Additionally, we validated the STR calls against known allele sizes in a set of GeT-RM reference cell-line materials (n = 6). Several STR loci that are entirely guanine or cytosines (G or C) have insufficient read evidence for inference and therefore could not be assayed precisely by TREDPARSE. TREDPARSE extends the limit of STR size detection beyond the physical sequence read length. This extension is critical because many of the disease risk cutoffs are close to or beyond the short sequence read length of 100 to 150 bases.

Phenotype Characterization of HD Intermediate Alleles in PREDICT-HD

BACKGROUND

Huntington disease (HD) is a neurodegenerative disease caused by a CAG repeat expansion on chromosome 4. Pathology is associated with CAG repeat length. Prior studies examining people in the intermediate allele (IA) range found subtle differences in motor, cognitive, and behavioral domains compared to controls.

OBJECTIVE

The purpose of this study was to examine baseline and longitudinal differences in motor, cognitive, behavioral, functional, and imaging outcomes between persons with CAG repeats in three ranges: normal (≤26), intermediate (27-35), and reduced penetrance (36-39).

METHODS

We examined longitudinal data from 389 participants in three allele groups: 280 normal controls (NC), 21 intermediate allele [IA], and 88 reduced penetrance [RP]. We used linear mixed models to identify differences in baseline and longitudinal outcomes between groups. Three models were tested: 1) no baseline or longitudinal differences; 2) baseline differences but no longitudinal differences; and 3) baseline and longitudinal differences.

RESULTS

Model 1 was the best fitting model for most outcome variables. Models 2 and 3 were best fitting for some of the variables. We found baseline and longitudinal trends of declining performance across increasing CAG repeat length groups, but no significant differences between the NC and IA groups.

CONCLUSION

We did not find evidence to support differences in the IA group compared to the NC group. These findings are limited by a small IA sample size.

Satellite DNA: An Evolving Topic

Satellite DNA represents one of the most fascinating parts of the repetitive fraction of the eukaryotic genome. Since the discovery of highly repetitive tandem DNA in the 1960s, a lot of literature has extensively covered various topics related to the structure, organization, function, and evolution of such sequences. Today, with the advent of genomic tools, the study of satellite DNA has regained a great interest. Thus, Next-Generation Sequencing (NGS), together with high-throughput in silico analysis of the information contained in NGS reads, has revolutionized the analysis of the repetitive fraction of the eukaryotic genomes. The whole of the historical and current approaches to the topic gives us a broad view of the function and evolution of satellite DNA and its role in chromosomal evolution. Currently, we have extensive information on the molecular, chromosomal, biological, and population factors that affect the evolutionary fate of satellite DNA, knowledge that gives rise to a series of hypotheses that get on well with each other about the origin, spreading, and evolution of satellite DNA. In this paper, I review these hypotheses from a methodological, conceptual, and historical perspective and frame them in the context of chromosomal organization and evolution.

CoLoRMap: Correcting Long Reads by Mapping short reads

MOTIVATION

Second generation sequencing technologies paved the way to an exceptional increase in the number of sequenced genomes, both prokaryotic and eukaryotic. However, short reads are difficult to assemble and often lead to highly fragmented assemblies. The recent developments in long reads sequencing methods offer a promising way to address this issue. However, so far long reads are characterized by a high error rate, and assembling from long reads require a high depth of coverage. This motivates the development of hybrid approaches that leverage the high quality of short reads to correct errors in long reads.

RESULTS

We introduce CoLoRMap, a hybrid method for correcting noisy long reads, such as the ones produced by PacBio sequencing technology, using high-quality Illumina paired-end reads mapped onto the long reads. Our algorithm is based on two novel ideas: using a classical shortest path algorithm to find a sequence of overlapping short reads that minimizes the edit score to a long read and extending corrected regions by local assembly of unmapped mates of mapped short reads. Our results on bacterial, fungal and insect data sets show that CoLoRMap compares well with existing hybrid correction methods.

AVAILABILITY AND IMPLEMENTATION

The source code of CoLoRMap is freely available for non-commercial use at https://github.com/sfu-compbio/colormap

CONTACT

ehaghshe@sfu.ca or cedric.chauve@sfu.ca

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

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.

A genomic view of short tandem repeats

Short tandem repeats (STRs) are some of the fastest mutating loci in the genome. Tools for accurately profiling STRs from high-throughput sequencing data have enabled genome-wide interrogation of more than a million STRs across hundreds of individuals. These catalogs have revealed that STRs are highly multiallelic and may contribute more de novo mutations than any other variant class. Recent studies have leveraged these catalogs to show that STRs play a widespread role in regulating gene expression and other molecular phenotypes. These analyses suggest that STRs are an underappreciated but rich reservoir of variation that likely make significant contributions to Mendelian diseases, complex traits, and cancer.

Long Single-Molecule Reads Can Resolve the Complexity of the Influenza Virus Composed of Rare, Closely Related Mutant Variants

As a result of a high rate of mutations and recombination events, an RNA-virus exists as a heterogeneous "swarm" of mutant variants. The long read length offered by single-molecule sequencing technologies allows each mutant variant to be sequenced in a single pass. However, high error rate limits the ability to reconstruct heterogeneous viral population composed of rare, related mutant variants. In this article, we present two single-nucleotide variants (2SNV), a method able to tolerate the high error rate of the single-molecule protocol and reconstruct mutant variants. 2SNV uses linkage between single-nucleotide variations to efficiently distinguish them from read errors. To benchmark the sensitivity of 2SNV, we performed a single-molecule sequencing experiment on a sample containing a titrated level of known viral mutant variants. Our method is able to accurately reconstruct clone with frequency of 0.2% and distinguish clones that differed in only two nucleotides distantly located on the genome. 2SNV outperforms existing methods for full-length viral mutant reconstruction.

Epidemiology and molecular mechanism of frontotemporal lobar degeneration/amyotrophic lateral sclerosis with repeat expansion mutation in C9orf72

GGGGCC hexanucleotide repeat expansions in C9orf72 were identified in 2011 as the genetic cause of frontotemporal lobar degeneration (FTLD)/amyotrophic lateral sclerosis (ALS) linked to chromosome 9. Since then, a number of studies have been conducted to delineate the molecular epidemiology of the repeat expansions and the molecular pathophysiology of the disease. The frequency of the repeat expansions considerably varied among countries. The frequency of the repeat expansions was high in European populations and populations of European descent and a substantial proportion of sporadic FTLD or ALS patients also have the mutations in these populations. On the other hand, the frequency was extremely low in Asia or Oceania except for limited regions including Kii Peninsula of Japan. A founder effect seems to strongly influence the regional differences in the frequency, but there is no definitive evidence that supports the notion that the repeat expansions arose in a single founder or multiple founders. As a disease-causing mechanism, several molecular mechanisms have been proposed, including conformational changes of DNA (G-quadruplex formation and hypermethylation) or RNA (G-quadruplex formation) molecules, altered transcriptional levels of C9orf72, sequestration of RNA-binding proteins, bidirectional transcription, formation of RNA foci, and neurotoxicity of dipeptide repeat proteins generated by repeat-associated non-ATG-initiated translation. Further investigations on the molecular mechanisms of neurodegeneration are expected to lead to the development of therapeutic interventions for this disease as well as for other diseases associated with non-coding repeat expansions.

Challenges, Solutions, and Quality Metrics of Personal Genome Assembly in Advancing Precision Medicine

Even though each of us shares more than 99% of the DNA sequences in our genome, there are millions of sequence codes or structure in small regions that differ between individuals, giving us different characteristics of appearance or responsiveness to medical treatments. Currently, genetic variants in diseased tissues, such as tumors, are uncovered by exploring the differences between the reference genome and the sequences detected in the diseased tissue. However, the public reference genome was derived with the DNA from multiple individuals. As a result of this, the reference genome is incomplete and may misrepresent the sequence variants of the general population. The more reliable solution is to compare sequences of diseased tissue with its own genome sequence derived from tissue in a normal state. As the price to sequence the human genome has dropped dramatically to around $1000, it shows a promising future of documenting the personal genome for every individual. However, de novo assembly of individual genomes at an affordable cost is still challenging. Thus, till now, only a few human genomes have been fully assembled. In this review, we introduce the history of human genome sequencing and the evolution of sequencing platforms, from Sanger sequencing to emerging "third generation sequencing" technologies. We present the currently available de novo assembly and post-assembly software packages for human genome assembly and their requirements for computational infrastructures. We recommend that a combined hybrid assembly with long and short reads would be a promising way to generate good quality human genome assemblies and specify parameters for the quality assessment of assembly outcomes. We provide a perspective view of the benefit of using personal genomes as references and suggestions for obtaining a quality personal genome. Finally, we discuss the usage of the personal genome in aiding vaccine design and development, monitoring host immune-response, tailoring drug therapy and detecting tumors. We believe the precision medicine would largely benefit from bioinformatics solutions, particularly for personal genome assembly.

Polymorphic tandem repeats within gene promoters act as modifiers of gene expression and DNA methylation in humans

Despite representing an important source of genetic variation, tandem repeats (TRs) remain poorly studied due to technical difficulties. We hypothesized that TRs can operate as expression (eQTLs) and methylation (mQTLs) quantitative trait loci. To test this we analyzed the effect of variation at 4849 promoter-associated TRs, genotyped in 120 individuals, on neighboring gene expression and DNA methylation. Polymorphic promoter TRs were associated with increased variance in local gene expression and DNA methylation, suggesting functional consequences related to TR variation. We identified >100 TRs associated with expression/methylation levels of adjacent genes. These potential eQTL/mQTL TRs were enriched for overlaps with transcription factor binding and DNaseI hypersensitivity sites, providing a rationale for their effects. Moreover, we showed that most TR variants are poorly tagged by nearby single nucleotide polymorphisms (SNPs) markers, indicating that many functional TR variants are not effectively assayed by SNP-based approaches. Our study assigns biological significance to TR variations in the human genome, and suggests that a significant fraction of TR variations exert functional effects via alterations of local gene expression or epigenetics. We conclude that targeted studies that focus on genotyping TR variants are required to fully ascertain functional variation in the genome.

PacBio Sequencing and Its Applications

Single-molecule, real-time sequencing developed by Pacific BioSciences offers longer read lengths than the second-generation sequencing (SGS) technologies, making it well-suited for unsolved problems in genome, transcriptome, and epigenetics research. The highly-contiguous de novo assemblies using PacBio sequencing can close gaps in current reference assemblies and characterize structural variation (SV) in personal genomes. With longer reads, we can sequence through extended repetitive regions and detect mutations, many of which are associated with diseases. Moreover, PacBio transcriptome sequencing is advantageous for the identification of gene isoforms and facilitates reliable discoveries of novel genes and novel isoforms of annotated genes, due to its ability to sequence full-length transcripts or fragments with significant lengths. Additionally, PacBio’s sequencing technique provides information that is useful for the direct detection of base modifications, such as methylation. In addition to using PacBio sequencing alone, many hybrid sequencing strategies have been developed to make use of more accurate short reads in conjunction with PacBio long reads. In general, hybrid sequencing strategies are more affordable and scalable especially for small-size laboratories than using PacBio Sequencing alone. The advent of PacBio sequencing has made available much information that could not be obtained via SGS alone.

Next-Generation Sequencing Approaches in Cancer: Where Have They Brought Us and Where Will They Take Us?

Next-generation sequencing (NGS) technologies and data have revolutionized cancer research and are increasingly being deployed to guide clinicians in treatment decision-making. NGS technologies have allowed us to take an “omics” approach to cancer in order to reveal genomic, transcriptomic, and epigenomic landscapes of individual malignancies. Integrative multi-platform analyses are increasingly used in large-scale projects that aim to fully characterize individual tumours as well as general cancer types and subtypes. In this review, we examine how NGS technologies in particular have contributed to “omics” approaches in cancer research, allowing for large-scale integrative analyses that consider hundreds of tumour samples. These types of studies have provided us with an unprecedented wealth of information, providing the background knowledge needed to make small-scale (including “N of 1”) studies informative and relevant. We also take a look at emerging opportunities provided by NGS and state-of-the-art third-generation sequencing technologies, particularly in the context of translational research. Cancer research and care are currently poised to experience significant progress catalyzed by accessible sequencing technologies that will benefit both clinical- and research-based efforts.

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

A large, non-coding ATTCT repeat expansion causes the neurodegenerative disorder, spinocerebellar ataxia type 10 (SCA10). In a subset of SCA10 patients, interruption motifs are present at the 5’ end of the expansion and strongly correlate with epileptic seizures. Thus, interruption motifs are a predictor of the epileptic phenotype and are hypothesized to act as a phenotypic modifier in SCA10. Yet, the exact internal sequence structure of SCA10 expansions remains unknown due to limitations in current technologies for sequencing across long extended tracts of tandem nucleotide repeats. We used the third generation sequencing technology, Single Molecule Real Time (SMRT) sequencing, to obtain full-length contiguous expansion sequences, ranging from 2.5 to 4.4 kb in length, from three SCA10 patients with different clinical presentations. We obtained sequence spanning the entire length of the expansion and identified the structure of known and novel interruption motifs within the SCA10 expansion. The exact interruption patterns in expanded SCA10 alleles will allow us to further investigate the potential contributions of these interrupting sequences to the pathogenic modification leading to the epilepsy phenotype in SCA10. Our results also demonstrate that SMRT sequencing is useful for deciphering long tandem repeats that pose as “gaps” in the human genome sequence.

Whole genome analyses reveal no pathogenetic single nucleotide or structural differences between monozygotic twins discordant for amyotrophic lateral sclerosis

The contribution of genetic and environmental factors to the pathogenesis of sporadic amyotrophic lateral sclerosis (ALS) remains unclear. To investigate the genetic component of the disease, we performed whole genome sequencing on ALS discordant monozygotic twins. Illumina whole genome sequencing on white blood cell DNA of five ALS-discordant monozygotic twin pairs (10 samples in total) yielded ∼30x coverage per individual. All single nucleotide variants, indels, and structural variants (copy number variants, inversions and translocations) were called and evaluated for functional consequence, evolutionary conservation, population frequency and overlap with known ALS associated variants and genes. Results showed that no validated discordant coding or regulatory single nucleotide variants or indels were found, and nor were any genome-wide discordant structural variants detected. Concordant variants of particular interest were: 1) two rare, highly-conserved heterozygous non-synonymous variants in SYT9 and EWSR1, genes previously associated with ALS (out of 2044 rare heterozygous variants detected); 2) three rare homozygous missense variants; and 3) three novel copy number deletions that overlapped genes. In conclusion, no convincing coding or regulatory nucleotide or genome-wide structural differences were found between ALS discordant monozygotic twins. The results suggest that more work is needed to elucidate possible environmental, epigenetic, oligogenic and somatic genetic factors that could underlie susceptibility to sporadic ALS.

The DDBJ Japanese Genotype-phenotype Archive for genetic and phenotypic human data

The DNA Data Bank of Japan Center (DDBJ Center; http://www.ddbj.nig.ac.jp) maintains and provides public archival, retrieval and analytical services for biological information. Since October 2013, DDBJ Center has operated the Japanese Genotype-phenotype Archive (JGA) in collaboration with our partner institute, the National Bioscience Database Center (NBDC) of the Japan Science and Technology Agency. DDBJ Center provides the JGA database system which securely stores genotype and phenotype data collected from individuals whose consent agreements authorize data release only for specific research use. NBDC has established guidelines and policies for sharing human-derived data and reviews data submission and usage requests from researchers. In addition to the JGA project, DDBJ Center develops Semantic Web technologies for data integration and sharing in collaboration with the Database Center for Life Science. This paper describes the overview of the JGA project, updates to the DDBJ databases, and services for data retrieval, analysis and integration.

GFP-based fluorescence assay for CAG repeat instability in cultured human cells

Trinucleotide repeats can be highly unstable, mutating far more frequently than point mutations. Repeats typically mutate by addition or loss of units of the repeat. CAG repeat expansions in humans trigger neurological diseases that include myotonic dystrophy, Huntington disease, and several spinocerebellar ataxias. In human cells, diverse mechanisms promote CAG repeat instability, and in mice, the mechanisms of instability are varied and tissue-dependent. Dissection of mechanistic complexity and discovery of potential therapeutics necessitates quantitative and scalable screens for repeat mutation. We describe a GFP-based assay for screening modifiers of CAG repeat instability in human cells. The assay exploits an engineered intronic CAG repeat tract that interferes with expression of an inducible GFP minigene. Like the phenotypes of many trinucleotide repeat disorders, we find that GFP function is impaired by repeat expansion, in a length-dependent manner. The intensity of fluorescence varies inversely with repeat length, allowing estimates of repeat tract changes in live cells. We validate the assay using transcription through the repeat and engineered CAG-specific nucleases, which have previously been reported to induce CAG repeat instability. The assay is relatively fast and should be adaptable to large-scale screens of chemical and shRNA libraries.

Resolving complex tandem repeats with long reads

MOTIVATION

Resolving tandemly repeated genomic sequences is a necessary step in improving our understanding of the human genome. Short tandem repeats (TRs), or microsatellites, are often used as molecular markers in genetics, and clinically, variation in microsatellites can lead to genetic disorders like Huntington's diseases. Accurately resolving repeats, and in particular TRs, remains a challenging task in genome alignment, assembly and variation calling. Though tools have been developed for detecting microsatellites in short-read sequencing data, these are limited in the size and types of events they can resolve. Single-molecule sequencing technologies may potentially resolve a broader spectrum of TRs given their increased length, but require new approaches given their significantly higher raw error profiles. However, due to inherent error profiles of the single-molecule technologies, these reads presents a unique challenge in terms of accurately identifying and estimating the TRs.

RESULTS

Here we present PacmonSTR, a reference-based probabilistic approach, to identify the TR region and estimate the number of these TR elements in long DNA reads. We present a multistep approach that requires as input, a reference region and the reference TR element. Initially, the TR region is identified from the long DNA reads via a 3-stage modified Smith-Waterman approach and then, expected number of TR elements is calculated using a pair-Hidden Markov Models-based method. Finally, TR-based genotype selection (or clustering: homozygous/heterozygous) is performed with Gaussian mixture models, using the Akaike information criteria, and coverage expectations.

Returning to more finished genomes

Genomic data have become commonplace in most branches of the biological sciences and have fundamentally altered the way research is conducted. However, the predominance of short-read sequence data from second-generation sequencing technologies has commonly resulted in fragmented and partial genomic data characteristics. In this opinion, I will highlight how long, unbiased reads from single molecule, real-time (SMRT) sequencing now allow for a return to more contiguous and comprehensive views of genomes.