Profiling the genome-wide landscape of tandem repeat expansions

A holistic perspective on herpes simplex virus (HSV) ecology and evolution

Herpes simplex viruses (HSV) cause chronic infection in humans that are characterized by periodic episodes of mucosal shedding and ulcerative disease. HSV causes millions of infections world-wide, with lifelong bouts of viral reactivation from latency in neuronal ganglia. Infected individuals experience different levels of disease severity and frequency of reactivation. There are two distantly related HSV species, with HSV-1 infections historically found most often in the oral niche and HSV-2 infections in the genital niche. Over the last two decades, HSV-1 has emerged as the leading cause of first-episode genital herpes in multiple countries. While HSV-1 has the highest level of genetic diversity among human alpha-herpesviruses, it is not yet known how quickly the HSV-1 viral population in a human host adapts over time, or if there are population bottlenecks associated with viral reactivation and/or transmission. It is also unknown how the ecological environments in which HSV infections occur influence their evolutionary trajectory, or that of co-occurring viruses and microbes. In this review, we explore how HSV accrues genetic diversity within each new infection, and yet maintains its ability to successfully infect most of the human population. A holistic examination of the ecological context of natural human infections can expand our awareness of how HSV adapts as it moves within and between human hosts, and reveal the complexity of these lifelong human-virus interactions. These insights may in turn suggest new areas of exploration for other chronic pathogens that successfully evolve and persist among their hosts.

Genetic architecture of autism spectrum disorder: Lessons from large-scale genomic studies

Autism spectrum disorder (ASD) is a neurodevelopmental disorder with a strong genetic component. Recently developed genomic technologies, including microarray and next-generation sequencing (NGS), have enabled researchers to genetic analyses aimed at identifying genetic variations associated with ASD and to elucidate the genetic architecture of the disorder. Large-scale microarray, exome sequencing analyses, and robust statistical methods have resulted in successful gene discovery and identification of high-confidence ASD genes from among de novo and inherited variants. Efforts have been made to understand the genetic architecture of ASD using whole-genome sequencing and genome-wide association studies aimed at identifying noncoding mutations and common variants associated with ASD. In addition, the development of systems biology approaches has resulted in the integration of genetic findings with functional genomic datasets, thereby providing a unique insight into the functional convergence of ASD risk genes and their neurobiology. In this review, we summarize the latest findings of ASD genetic studies involving large cohorts and discuss their implications in ASD neurobiology and in clinical practice.

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.

Genome-wide characterization of human minisatellite VNTRs: population-specific alleles and gene expression differences

Variable Number Tandem Repeats (VNTRs) are tandem repeat (TR) loci that vary in copy number across a population. Using our program, VNTRseek, we analyzed human whole genome sequencing datasets from 2770 individuals in order to detect minisatellite VNTRs, i.e., those with pattern sizes ≥7 bp. We detected 35 638 VNTR loci and classified 5676 as commonly polymorphic (i.e. with non-reference alleles occurring in >5% of the population). Commonly polymorphic VNTR loci were found to be enriched in genomic regions with regulatory function, i.e. transcription start sites and enhancers. Investigation of the commonly polymorphic VNTRs in the context of population ancestry revealed that 1096 loci contained population-specific alleles and that those could be used to classify individuals into super-populations with near-perfect accuracy. Search for quantitative trait loci (eQTLs), among the VNTRs proximal to genes, indicated that in 187 genes expression differences correlated with VNTR genotype. We validated our predictions in several ways, including experimentally, through the identification of predicted alleles in long reads, and by comparisons showing consistency between sequencing platforms. This study is the most comprehensive analysis of minisatellite VNTRs in the human population to date.

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.

Pervasive cis effects of variation in copy number of large tandem repeats on local DNA methylation and gene expression

Variable number tandem repeats (VNTRs) are composed of large tandemly repeated motifs, many of which are highly polymorphic in copy number. However, because of their large size and repetitive nature, they remain poorly studied. To investigate the regulatory potential of VNTRs, we used read-depth data from Illumina whole-genome sequencing to perform association analysis between copy number of ∼70,000 VNTRs (motif size ≥ 10 bp) with both gene expression (404 samples in 48 tissues) and DNA methylation (235 samples in peripheral blood), identifying thousands of VNTRs that are associated with local gene expression (eVNTRs) and DNA methylation levels (mVNTRs). Using an independent cohort, we validated 73%-80% of signals observed in the two discovery cohorts, while allelic analysis of VNTR length and CpG methylation in 30 Oxford Nanopore genomes gave additional support for mVNTR loci, thus providing robust evidence to support that these represent genuine associations. Further, conditional analysis indicated that many eVNTRs and mVNTRs act as QTLs independently of other local variation. We also observed strong enrichments of eVNTRs and mVNTRs for regulatory features such as enhancers and promoters. Using the Human Genome Diversity Panel, we define sets of VNTRs that show highly divergent copy numbers among human populations and show that these are enriched for regulatory effects and preferentially associate with genes that have been linked with human phenotypes through GWASs. Our study provides strong evidence supporting functional variation at thousands of VNTRs and defines candidate sets of VNTRs, copy number variation of which potentially plays a role in numerous human phenotypes.

TRTools: a toolkit for genome-wide analysis of tandem repeats

SUMMARY

A rich set of tools have recently been developed for performing genome-wide genotyping of tandem repeats (TRs). However, standardized tools for downstream analysis of these results are lacking. To facilitate TR analysis applications, we present TRTools, a Python library and suite of command line tools for filtering, merging and quality control of TR genotype files. TRTools utilizes an internal harmonization module, making it compatible with outputs from a wide range of TR genotypers.

AVAILABILITY AND IMPLEMENTATION

TRTools is freely available at https://github.com/gymreklab/TRTools. Detailed documentation is available at https://trtools.readthedocs.io.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

CANVAS: A New Genetic Entity in the Otorhinolaryngologist's Differential Diagnosis

OBJECTIVE

The biallelic inheritance of an expanded intronic pentamer (AAGGG)_exp in the gene encoding replication factor C subunit 1 (RFC1) has been found to be a cause of cerebellar ataxia, neuropathy, and vestibular areflexia syndrome (CANVAS). This study describes clinical and genetic features of our patients with clinical suspicion of the syndrome.

STUDY DESIGN

A retrospective descriptive study from an ataxia database comprising 500 patients.

SETTING

The study was performed at the Otorhinolaryngology Department of a hospital in the north of Spain.

METHODS

Specific genetic testing for CANVAS was performed in 13 patients with clinical suspicion of complete or incomplete syndrome. The clinical diagnosis was supported by quantitative vestibular hypofunction, cerebellar atrophy, and abnormal sensory nerve conduction testing.

RESULTS

Nine of 13 (69%) patients met clinical diagnostic criteria for definite CANVAS disease. The first manifestation of the syndrome was lower limb dysesthesia in 8 of 13 patients and gait imbalance in 5 of 13. Eleven of 13 (85%) patients were carriers of the biallelic (AAGGG)_exp in RFC1.

CONCLUSION

A genetic cause of CANVAS has recently been discovered. We propose genetic screening for biallelic expansions of the AAGGG pentamer of RFC1 in all patients with clinical suspicion of CANVAS, since accurate early diagnosis could improve the quality of life of these patients.

Systematic analysis of short tandem repeats in 38,095 exomes provides an additional diagnostic yield

PURPOSE

Expansions of a subset of short tandem repeats (STRs) have been implicated in approximately 30 different human genetic disorders. Despite extensive application of exome sequencing (ES) in routine diagnostic genetic testing, STRs are not routinely identified from these data.

METHODS

We assessed diagnostic utility of STR analysis in exome sequencing by applying ExpansionHunter to 2,867 exomes from movement disorder patients and 35,228 other clinical exomes.

RESULTS

We identified 38 movement disorder patients with a possible aberrant STR length. Validation by polymerase chain reaction (PCR) and/or repeat-primed PCR technologies confirmed the presence of aberrant expansion alleles for 13 (34%). For seven of these patients the genotype was compatible with the phenotypic description, resulting in a molecular diagnosis. We subsequently tested the remainder of our diagnostic ES cohort, including over 30 clinically and genetically heterogeneous disorders. Optimized manual curation yielded 167 samples with a likely aberrant STR length. Validations confirmed 93/167 (56%) aberrant expansion alleles, of which 48 were in the pathogenic range and 45 in the premutation range.

CONCLUSION

Our work provides guidance for the implementation of STR analysis in clinical ES. Our results show that systematic STR evaluation may increase diagnostic ES yield by 0.2%, and recommend making STR evaluation a routine part of ES interpretation in genetic testing laboratories.

Variable number tandem repeats mediate the expression of proximal genes

Variable number tandem repeats (VNTRs) account for significant genetic variation in many organisms. In humans, VNTRs have been implicated in both Mendelian and complex disorders, but are largely ignored by genomic pipelines due to the complexity of genotyping and the computational expense. We describe adVNTR-NN, a method that uses shallow neural networks to genotype a VNTR in 18 seconds on 55X whole genome data, while maintaining high accuracy. We use adVNTR-NN to genotype 10,264 VNTRs in 652 GTEx individuals. Associating VNTR length with gene expression in 46 tissues, we identify 163 "eVNTRs". Of the 22 eVNTRs in blood where independent data is available, 21 (95%) are replicated in terms of significance and direction of association. 49% of the eVNTR loci show a strong and likely causal impact on the expression of genes and 80% have maximum effect size at least 0.3. The impacted genes are involved in diseases including Alzheimer's, obesity and familial cancers, highlighting the importance of VNTRs for understanding the genetic basis of complex diseases.

Laboratory testing for fragile X, 2021 revision: a technical standard of the American College of Medical Genetics and Genomics (ACMG)

Molecular genetic testing of the FMR1 gene is commonly performed in clinical laboratories. Pathogenic variants in the FMR1 gene are associated with fragile X syndrome, fragile X-associated tremor ataxia syndrome (FXTAS), and fragile X-associated primary ovarian insufficiency (FXPOI). This document provides updated information regarding FMR1 pathogenic variants, including prevalence, genotype-phenotype correlations, and variant nomenclature. Methodological considerations are provided for Southern blot analysis and polymerase chain reaction (PCR) amplification of FMR1, including triplet repeat-primed and methylation-specific PCR.The American College of Medical Genetics and Genomics (ACMG) Laboratory Quality Assurance Committee has the mission of maintaining high technical standards for the performance and interpretation of genetic tests. In part, this is accomplished by the publication of the document ACMG Technical Standards for Clinical Genetics Laboratories, which is now maintained online ( http://www.acmg.net ). This subcommittee also reviews the outcome of national proficiency testing in the genetics area and may choose to focus on specific diseases or methodologies in response to those results. Accordingly, the subcommittee selected fragile X syndrome to be the first topic in a series of supplemental sections, recognizing that it is one of the most frequently ordered genetic tests and that it has many alternative methods with different strengths and weaknesses. This document is the fourth update to the original standards and guidelines for fragile X testing that were published in 2001, with revisions in 2005 and 2013, respectively.This versionClarifies the clinical features associated with different FMRI variants (Section 2.3)Discusses important reporting considerations (Section 3.3.1.3)Provides updates on technology (Section 4.1).

Ethnically biased microsatellites contribute to differential gene expression and glutathione metabolism in Africans and Europeans

Approximately three percent of the human genome is occupied by microsatellites: a type of short tandem repeat (STR). Microsatellites have well established effects on (a) the genetic structure of diverse human populations and (b) expression of nearby genes. These lines of inquiry have uncovered 3,984 ethnically biased microsatellite loci (EBML) and 28,375 expression STRs (eSTRs), respectively. We hypothesize that a combination of EBML, eSTRs, and gene expression data (RNA-seq) can be used to show that microsatellites contribute to differential gene expression and phenotype in human populations. In fact, our previous study demonstrated a degree of mutual overlap between EBML and eSTRs but fell short of quantifying effects on gene expression. The present work aims to narrow the gap. First, we identify 313 overlapping EBML/eSTRs and recapitulate their mutual overlap. The 313 EBML/eSTRs are then characterized across ethnicity and tissue type. We use RNA-seq data to pursue validation of 49 regions that affect whole blood gene expression; 32 out of 54 affected genes are differentially expressed in Africans and Europeans. We quantify the relative contribution of these 32 genes to differential expression; fold change tends to be less than other differentially expressed genes. Repeat length correlates with expression for 15 of the 32 genes; two are conspicuously involved in glutathione metabolism. Finally, we repurpose a mathematical model of glutathione metabolism to investigate how a single polymorphic microsatellite affects phenotype. We conclude with a testable prediction that microsatellite polymorphisms affect GPX7 expression and oxidative stress in Africans and Europeans.

Abundancy of polymorphic CGG repeats in the human genome suggest a broad involvement in neurological disease

Expanded CGG-repeats have been linked to neurodevelopmental and neurodegenerative disorders, including the fragile X syndrome and fragile X-associated tremor/ataxia syndrome (FXTAS). We hypothesized that as of yet uncharacterised CGG-repeat expansions within the genome contribute to human disease. To catalogue the CGG-repeats, 544 human whole genomes were analyzed. In total, 6101 unique CGG-repeats were detected of which more than 93% were highly variable in repeat length. Repeats with a median size of 12 repeat units or more were always polymorphic but shorter repeats were often polymorphic, suggesting a potential intergenerational instability of the CGG region even for repeats units with a median length of four or less. 410 of the CGG repeats were associated with known neurodevelopmental disease genes or with strong candidate genes. Based on their frequency and genomic location, CGG repeats may thus be a currently overlooked cause of human disease.

Genome diversity in Ukraine

Background

The main goal of this collaborative effort is to provide genome-wide data for the previously underrepresented population in Eastern Europe, and to provide cross-validation of the data from genome sequences and genotypes of the same individuals acquired by different technologies. We collected 97 genome-grade DNA samples from consented individuals representing major regions of Ukraine that were consented for public data release. BGISEQ-500 sequence data and genotypes by an Illumina GWAS chip were cross-validated on multiple samples and additionally referenced to 1 sample that has been resequenced by Illumina NovaSeq6000 S4 at high coverage.

Results

The genome data have been searched for genomic variation represented in this population, and a number of variants have been reported: large structural variants, indels, copy number variations, single-nucletide polymorphisms, and microsatellites. To our knowledge, this study provides the largest to-date survey of genetic variation in Ukraine, creating a public reference resource aiming to provide data for medical research in a large understudied population.

Conclusions

Our results indicate that the genetic diversity of the Ukrainian population is uniquely shaped by evolutionary and demographic forces and cannot be ignored in future genetic and biomedical studies. These data will contribute a wealth of new information bringing forth a wealth of novel, endemic and medically related alleles.

Patterns of de novo tandem repeat mutations and their role in autism

Autism Spectrum Disorder (ASD) is an early onset developmental disorder characterized by deficits in communication and social interaction and restrictive or repetitive behaviors^1,2. Family studies demonstrate that ASD has a significant genetic basis with contributions both from inherited and de novo variants^3,4. It has been estimated that de novo mutations may contribute to 30% of all simplex cases, in which only a single child is affected per family⁵. Tandem repeats (TRs), defined here as 1-20bp sequences repeated consecutively, comprise one of the largest sources of de novo mutations in humans⁶. TR expansions are implicated in dozens of neurological and psychiatric disorders⁷. Yet, de novo TR mutations have not been characterized on a genome-wide scale, and their contribution to ASD remains unexplored. We develop novel bioinformatics methods for identifying and prioritizing de novo TR mutations from sequencing data and then perform a genome-wide characterization of de novo TR mutations in ASD-affected probands and unaffected siblings. Compared to recent work on TRs in ASD⁸, we explicitly infer mutation events and their precise changes in repeat number, and primarily focus on more prevalent stepwise copy number changes rather than large expansions. Our results demonstrate a significant genome-wide excess of TR mutations in ASD probands. Mutations in probands tend to be larger, enriched in fetal brain regulatory regions, and predicted to be more evolutionarily deleterious. Overall, our results highlight the importance of considering repeat variants in future studies of de novo mutations.

Genetic testing in dementia - utility and clinical strategies

Techniques for clinical genetic testing in dementia disorders have advanced rapidly but remain to be more widely implemented in practice. A positive genetic test offers a precise molecular diagnosis, can help members of an affected family to determine personal risk, provides a basis for reproductive choices and can offer options for clinical trials. The likelihood of identifying a specific genetic cause of dementia depends on the clinical condition, the age at onset and family history. Attempts to match phenotypes to single genes are mostly inadvisable owing to clinical overlap between the dementias, genetic heterogeneity, pleiotropy and concurrent mutations. Currently, the appropriate genetic test in most cases of dementia is a next-generation sequencing gene panel, though some conditions necessitate specific types of test such as repeat expansion testing. Whole-exome and whole-genome sequencing are becoming financially feasible but raise or exacerbate complex issues such as variants of uncertain significance, secondary findings and the potential for re-analysis in light of new information. However, the capacity for data analysis and counselling is already restricting the provision of genetic testing. Patients and their relatives need to be given reliable information to enable them to make informed choices about tests, treatments and data sharing; the ability of patients with dementia to make decisions must be considered when providing this information.

Genetics and Genomics Approaches for Diagnosis and Research Into Hereditary Anemias

The hereditary anemias are a relatively heterogeneous set of disorders that can show wide clinical and genetic heterogeneity, which often hampers correct clinical diagnosis. The classical diagnostic workflow for these conditions generally used to start with analysis of the family and personal histories, followed by biochemical and morphological evaluations, and ending with genetic testing. However, the diagnostic framework has changed more recently, and genetic testing is now a suitable approach for differential diagnosis of these patients. There are several approaches to this genetic testing, the choice of which depends on phenotyping, genetic heterogeneity, and gene size. For patients who show complete phenotyping, single-gene testing remains recommended. However, genetic analysis now includes next-generation sequencing, which is generally based on custom-designed targeting panels and whole-exome sequencing. The use of next-generation sequencing also allows the identification of new causative genes, and of polygenic conditions and genetic factors that modify disease severity of hereditary anemias. In the research field, whole-genome sequencing is useful for the identification of non-coding causative mutations, which might account for the disruption of transcriptional factor occupancy sites and cis-regulatory elements. Moreover, advances in high-throughput sequencing techniques have now resulted in the identification of genome-wide profiling of the chromatin structures known as the topologically associating domains. These represent a recurrent disease mechanism that exposes genes to inappropriate regulatory elements, causing errors in gene expression. This review focuses on the challenges of diagnosis and research into hereditary anemias, with indications of both the advantages and disadvantages. Finally, we consider the future perspectives for the use of next-generation sequencing technologies in this era of precision medicine.

Biallelic Intronic AAGGG Expansion of RFC1 is Related to Multiple System Atrophy

OBJECTIVE

A recessive biallelic repeat expansion, (AAGGG)_exp , in the RFC1 gene has been reported to be a frequent cause of late-onset ataxia. For cerebellar ataxia, neuropathy, and vestibular areflexia syndrome (CANVAS), the recessive biallelic (AAGGG)_exp genotype was present in ~92% of cases. This study aimed to examine whether the pentanucleotide repeat (PNR) was related to multiple system atrophy (MSA), which shares a spectrum of symptoms with CANVAS.

METHODS

In this study, we screened the pathogenic (AAGGG)_exp repeat and 5 other PNRs in 104 Chinese sporadic adult-onset ataxia of unknown aetiology (SAOA) patients, 282 MSA patients, and 203 unaffected individuals. Multiple molecular genetic tests were used, including long-range polymerase chain reaction (PCR), repeat-primed PCR (RP-PCR), Sanger sequencing, and Southern blot. Comprehensive clinical assessments were conducted, including neurological examination, neuroimaging, nerve electrophysiology, and examination of vestibular function.

RESULTS

We identified biallelic (AAGGG)_exp in 1 SAOA patient and 3 MSA patients. Additionally, 1 MSA patient had the (AAGGG)_exp /(AAAGG)_exp genotype with uncertain pathogenicity. We also described the carrier frequency for different PNRs in our cohorts. Furthermore, we summarized the distinct phenotypes of affected patients, suggesting that biallelic (AAGGG)_exp in RFC1 could be associated with MSA and should be screened routinely in the MSA diagnostic workflow.

INTERPRETATION

Our results expanded the clinical phenotypic spectrum of RFC1-related disorders and raised the possibility that MSA might share the same genetic background as CANVAS, which is crucial for re-evaluating the current CANVAS and MSA diagnostic criteria. ANN NEUROL 2020;88:1132-1143.

popSTR2 enables clinical and population-scale genotyping of microsatellites

Summary

popSTR2 is an update and augmentation of our previous work ‘popSTR: a population-based microsatellite genotyper’. To make genotyping sensitive to inter-sample differences, we supply a kernel to estimate sample-specific slippage rates. For clinical sequencing purposes, a panel of known pathogenic repeat expansions is provided along with a script that scans and flags for manual inspection markers indicative of a pathogenic expansion. Like its predecessor, popSTR2 allows for joint genotyping of samples at a population scale. We now provide a binning method that makes the microsatellite genotypes more amenable to analysis within standard association pipelines and can increase association power.

Availability and implementation

https://github.com/DecodeGenetics/popSTR.

Supplementary information

Supplementary data are available at Bioinformatics online.

Large scale in silico characterization of repeat expansion variation in human genomes

Significant progress has been made in elucidating single nucleotide polymorphism diversity in the human population. However, the majority of the variation space in the genome is structural and remains partially elusive. One form of structural variation is tandem repeats (TRs). Expansion of TRs are responsible for over 40 diseases, but we hypothesize these represent only a fraction of the pathogenic repeat expansions that exist. Here we characterize long or expanded TR variation in 1,115 human genomes as well as a replication cohort of 2,504 genomes, identified using ExpansionHunter Denovo. We found that individual genomes typically harbor several rare, large TRs, generally in non-coding regions of the genome. We noticed that these large TRs are enriched in their proximity to Alu elements. The vast majority of these large TRs seem to be expansions of smaller TRs that are already present in the reference genome. We are providing this TR profile as a resource for comparison to undiagnosed rare disease genomes in order to detect novel disease-causing repeat expansions.

Chemical Probing of Thymine in the TGG/CGG Triad to Explore the Deamination of 5-Methylcytosine in the CGG Repeat

The methylation of cytosine in the full mutation of the expanded CGG repeat and subsequent deamination to thymine could be a measure of repeat instability. We report the synthesis of NCD-Bpy, which binds to the TGG/CGG site in the repeat hairpin. NCD-Bpy forces the thymine in the TGG/CGG site to flip out from the π-stack, recruits osmium tetroxide in the vicinity of the flipped-out T, and oxidizes the T. The piperidine-induced cleavage band successfully determined the position of the T in the expanded CGG repeat.

ExpansionHunter: a sequence-graph-based tool to analyze variation in short tandem repeat regions

SUMMARY

We describe a novel computational method for genotyping repeats using sequence graphs. This method addresses the long-standing need to accurately genotype medically important loci containing repeats adjacent to other variants or imperfect DNA repeats such as polyalanine repeats. Here we introduce a new version of our repeat genotyping software, ExpansionHunter, that uses this method to perform targeted genotyping of a broad class of such loci.

AVAILABILITY AND IMPLEMENTATION

ExpansionHunter is implemented in C++ and is available under the Apache License Version 2.0. The source code, documentation, and Linux/macOS binaries are available at https://github.com/Illumina/ExpansionHunter/.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

ExpansionHunter Denovo: a computational method for locating known and novel repeat expansions in short-read sequencing data

Repeat expansions are responsible for over 40 monogenic disorders, and undoubtedly more pathogenic repeat expansions remain to be discovered. Existing methods for detecting repeat expansions in short-read sequencing data require predefined repeat catalogs. Recent discoveries emphasize the need for methods that do not require pre-specified candidate repeats. To address this need, we introduce ExpansionHunter Denovo, an efficient catalog-free method for genome-wide repeat expansion detection. Analysis of real and simulated data shows that our method can identify large expansions of 41 out of 44 pathogenic repeats, including nine recently reported non-reference repeat expansions not discoverable via existing methods.

The Genetics of Epilepsy

Epilepsy encompasses a group of heterogeneous brain diseases that affect more than 50 million people worldwide. Epilepsy may have discernible structural, infectious, metabolic, and immune etiologies; however, in most people with epilepsy, no obvious cause is identifiable. Based initially on family studies and later on advances in gene sequencing technologies and computational approaches, as well as the establishment of large collaborative initiatives, we now know that genetics plays a much greater role in epilepsy than was previously appreciated. Here, we review the progress in the field of epilepsy genetics and highlight molecular discoveries in the most important epilepsy groups, including those that have been long considered to have a nongenetic cause. We discuss where the field of epilepsy genetics is moving as it enters a new era in which the genetic architecture of common epilepsies is starting to be unraveled.

Accuracy of short tandem repeats genotyping tools in whole exome sequencing data

Background: Short tandem repeats are an important source of genetic variation. They are highly mutable and repeat expansions are associated dozens of human disorders, such as Huntington's disease and spinocerebellar ataxias. Technical advantages in sequencing technology have made it possible to analyse these repeats at large scale; however, accurate genotyping is still a challenging task. We compared four different short tandem repeats genotyping tools on whole exome sequencing data to determine their genotyping performance and limits, which will aid other researchers in choosing a suitable tool and parameters for analysis. Methods: The analysis was performed on the Simons Simplex Collection dataset, where we used a novel method of evaluation with accuracy determined by the rate of homozygous calls on the X chromosome of male samples. In total we analysed 433 samples and around a million genotypes for evaluating tools on whole exome sequencing data. Results: We determined a relatively good performance of all tools when genotyping repeats of 3-6 bp in length, which could be improved with coverage and quality score filtering. However, genotyping homopolymers was challenging for all tools and a high error rate was present across different thresholds of coverage and quality scores. Interestingly, dinucleotide repeats displayed a high error rate as well, which was found to be mainly caused by the AC/TG repeats. Overall, LobSTR was able to make the most calls and was also the fastest tool, while RepeatSeq and HipSTR exhibited the lowest heterozygous error rate at low coverage. Conclusions: All tools have different strengths and weaknesses and the choice may depend on the application. In this analysis we demonstrated the effect of using different filtering parameters and offered recommendations based on the trade-off between the best accuracy of genotyping and the highest number of calls.

Familial adult myoclonic epilepsy type 1 SAMD12 TTTCA repeat expansion arose 17,000 years ago and is present in Sri Lankan and Indian families

Familial adult myoclonic epilepsy 1 (FAME1), first recognised in Japanese families, was recently shown to be caused by a TTTCA repeat insertion in intron 4 of SAMD12 on chromosome 8. We performed whole genome sequencing on two families with FAME, one of Sri Lankan origin and the other of Indian origin, and identified a TTTCA repeat insertion in SAMD12 in both families. Haplotype analysis revealed that both families shared the same core ancestral haplotype reported in Japanese and Chinese families with FAME1. Mutation dating, based on the length of shared haplotypes, estimated the age of the ancestral haplotype to be ~670 generations, or 17,000 years old. Our data extend the geographic range of this repeat expansion to Southern Asia and potentially implicate an even broader regional distribution given the age of the variant. This finding suggests patients of Asian ancestry with suspected FAME should be screened for the SAMD12 TTTCA expansion.

ExpansionHunter: a sequence-graph-based tool to analyze variation in short tandem repeat regions

SUMMARY

We describe a novel computational method for genotyping repeats using sequence graphs. This method addresses the long-standing need to accurately genotype medically important loci containing repeats adjacent to other variants or imperfect DNA repeats such as polyalanine repeats. Here we introduce a new version of our repeat genotyping software, ExpansionHunter, that uses this method to perform targeted genotyping of a broad class of such loci.

AVAILABILITY AND IMPLEMENTATION

ExpansionHunter is implemented in C++ and is available under the Apache License Version 2.0. The source code, documentation, and Linux/macOS binaries are available at https://github.com/Illumina/ExpansionHunter/.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Sequencing Technologies and Analyses: Where Have We Been and Where Are We Going?

A wave of technologies transformed sequencing over a decade ago into the high-throughput era, demanding research in new computational methods to analyze these data. The applications of these sequencing technologies have continuously expanded since then. The RECOMB Satellite Workshop on Massively Parallel Sequencing (RECOMB-Seq) meeting, established in 2011, brings together leading researchers in computational genomics and genomic biology to discuss emerging frontiers in algorithm development for massively parallel sequencing data. The ninth edition of this workshop was held in Washington, DC, in George Washington University on May 3 and 4, 2019. There was an exploration of several traditional topics in sequence analysis, including genome assembly, sequence alignment, and data compression, and development of methods for new sequencing technologies, including linked reads and single-molecule long-read sequencing. Here we revisit these topics and discuss the current status and perspectives of sequencing technologies and analyses.