A genomic view of short tandem repeats

The genomic landscape of short tandem repeats in cattle

Short tandem repeats (STRs) are abundant and have high mutation rates across cattle genomes; however, comprehensive exploration of cattle STRs is needed. Here, we constructed a comprehensive map of 467 553 polymorphic STRs (pSTRs) constructed from 423 cattle genomes representing 59 breeds worldwide. We observed that pSTRs in coding sequences and 5'UTRs (Untranslated Regions) were under strong selective constraints and exhibited a relatively low level of diversity. Furthermore, we found that these pSTRs underwent more contraction than expansion. Population analysis showed a strong positive correlation (R = 1) between pSTR diversity and single nucleotide polymorphic heterozygosity. We also investigated STR differences between taurine and indicine cattle and detected 2301 highly divergent STRs, which might relate to immune, endocrine and neurodevelopmental pathways. In summary, our large-scale study characterizes the spectrum of STRs in cattle, expands the scale of known cattle STR variation and provides novel insights into differences among various cattle subspecies.

Multi-omic quantitative trait loci link tandem repeat size variation to gene regulation in human brain

Tandem repeat (TR) size variation is implicated in ~50 neurological disorders, yet its impact on gene regulation in the human brain remains largely unknown. In the present study, we quantified the impact of TR size variation on brain gene regulation across distinct molecular phenotypes, based on 4,412 multi-omics samples from 1,597 donors, including 1,586 newly sequenced ones. We identified ~2.2 million TR molecular quantitative trait loci (TR-xQTLs), linking ~139,000 unique TRs to nearby molecular phenotypes, including many known disease-risk TRs, such as the G2C4 expansion in C9orf72 associated with amyotrophic lateral sclerosis. Fine-mapping revealed ~18,700 TRs as potential causal variants. Our in vitro experiments further confirmed the causal and independent regulatory effects of three TRs. Additional colocalization analysis indicated the potential causal role of TR variation in brain-related phenotypes, highlighted by a 3'-UTR TR in NUDT14 linked to cortical surface area and a TG repeat in PLEKHA1, associated with Alzheimer's disease.

Integrative genetic analysis: cornerstone of precision psychiatry

The role of genetic testing in the domain of neurodevelopmental and psychiatric disorders (NPDs) is gradually changing from providing etiological explanation for the presence of NPD phenotypes to also identifying young individuals at high risk of developing NPDs before their clinical manifestation. In clinical practice, the latter implies a shift towards the availability of individual genetic information predicting a certain liability to develop an NPD (e.g., autism, intellectual disability, psychosis etc.). The shift from mostly a posteriori explanation to increasingly a priori risk prediction is the by-product of the systematic implementation of whole exome or genome sequencing as part of routine diagnostic work-ups during the neonatal and prenatal periods. This rapid uptake of genetic testing early in development has far-reaching consequences for psychiatry: Whereas until recently individuals would come to medical attention because of signs of abnormal developmental and/or behavioral symptoms, increasingly, individuals are presented based on genetic liability for NPD outcomes before NPD symptoms emerge. This novel clinical scenario, while challenging, also creates opportunities for research on prevention interventions and precision medicine approaches. Here, we review why optimization of individual risk prediction is a key prerequisite for precision medicine in the sphere of NPDs, as well as the technological and statistical methods required to achieve this ambition.

Microsatellites explorer: A database of short tandem repeats across genomes

Short tandem repeats (STRs) are widespread, repetitive elements, with a number of biological functions and are among the most rapidly mutating regions in the genome. Their distribution varies significantly between taxonomic groups in the tree of life and are highly polymorphic within the human population. Advances in sequencing technologies coupled with decreasing costs have enabled the generation of an ever-growing number of complete genomes. Additionally, the arrival of accurate long reads has facilitated the generation of Telomere-to-Telomere (T2T) assemblies of complete genomes. Nevertheless, there is no comprehensive database that encompasses the STRs found per genome across different organisms and for different human genomes across diverse ancestries. Here we introduce Microsatellites Explorer, a database of STRs found in the genomes of 117,253 organisms across all major taxonomic groups, 15 T2T genome assemblies of different organisms, and 94 human haplotypes from the human pangenome. The database currently hosts 406,758,798 STR sequences, serving as a centralized user-friendly repository to perform searches, interactive visualizations, and download existing STR data for independent analysis. Microsatellites Explorer is implemented as a web-portal for browsing, analyzing and downloading STR data. Microsatellites Explorer is publicly available at https://www.microsatellitesexplorer.com.

ULTRA-effective labeling of tandem repeats in genomic sequence

In the age of long read sequencing, genomics researchers now have access to accurate repetitive DNA sequence (including satellites) that, due to the limitations of short read-sequencing, could previously be observed only as unmappable fragments. Tools that annotate repetitive sequence are now more important than ever, so that we can better understand newly uncovered repetitive sequences, and also so that we can mitigate errors in bioinformatic software caused by those repetitive sequences. To that end, we introduce the 1.0 release of our tool for identifying and annotating locally repetitive sequence, ULTRA Locates Tandemly Repetitive Areas (ULTRA). ULTRA is fast enough to use as part of an efficient annotation pipeline, produces state-of-the-art reliable coverage of repetitive regions containing many mutations, and provides interpretable statistics and labels for repetitive regions.

Availability and implementation

ULTRA is released under an open source license, and is available for download at https://github.com/TravisWheelerLab/ULTRA.

Repeat expansion disorders

An increasing number of repeat expansion disorders have been found to cause both rare and common neurological disease. This is exemplified in recent discoveries of novel repeat expansions underlying a significant proportion of several late-onset neurodegenerative disorders, such as CANVAS (cerebellar ataxia, neuropathy and vestibular areflexia syndrome) and spinocerebellar ataxia type 27B. Most of the 60 described repeat expansion disorders to date are associated with neurological disease, providing substantial challenges for diagnosis, but also opportunities for management in a clinical neurology setting. Commonalities in clinical presentation, overarching diagnostic features and similarities in the approach to genetic testing justify considering these disorders collectively based on their unifying causative mechanism. In this review, we discuss the characteristics and diagnostic challenges of repeat expansion disorders for the neurologist and provide examples to highlight their clinical heterogeneity. With the ready availability of clinical-grade whole-genome sequencing for molecular diagnosis, we discuss the current approaches to testing for repeat expansion disorders and application in clinical practice.

Mechanistic Characterization of De Novo Generation of Variable Number Tandem Repeats in Circular Plasmids during Site-Directed Mutagenesis and Optimization for Coding Gene Application

Site-directed mutagenesis for creating point mutations, sometimes, gives rise to plasmids carrying variable number tandem repeats (VNTRs) locally, which are arbitrarily regarded as polymerase chain reaction (PCR) related artifacts. Here, the alternative end-joining mechanism is reported rather than PCR artifacts accounts largely for that VNTRs formation and expansion. During generating a point mutation on GPLD1 gene, an unexpected formation of VNTRs employing the 31 bp mutagenesis primers is observed as the repeat unit in the pcDNA3.1-GPLD1 plasmid. The 31 bp VNTRs are formed in 24.75% of the resulting clones with copy number varied from 2 to 13. All repeat units are aligned with the same orientation as GPLD1 gene. 43.54% of the repeat junctions harbor nucleotide mutations while the rest don't. Their demonstrated short primers spanning the 3' part of the mutagenesis primers are essential for initial creation of the 2-copy tandem repeats (TRs) in circular plasmids. The dimerization of mutagenesis primers by the alternative end-joining in a correct orientation is required for further expansion of the 2-copy TRs. Lastly, a half-double priming strategy is established, verified the findings and offered a simple method for VNTRs creation on coding genes in circular plasmids without junction mutations.

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

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

Endogenous retroviruses mediate transcriptional rewiring in response to oncogenic signaling in colorectal cancer

Cancer cells exhibit rewired transcriptional regulatory networks that promote tumor growth and survival. However, the mechanisms underlying the formation of these pathological networks remain poorly understood. Through a pan-cancer epigenomic analysis, we found that primate-specific endogenous retroviruses (ERVs) are a rich source of enhancers displaying cancer-specific activity. In colorectal cancer and other epithelial tumors, oncogenic MAPK/AP1 signaling drives the activation of enhancers derived from the primate-specific ERV family LTR10. Functional studies in colorectal cancer cells revealed that LTR10 elements regulate tumor-specific expression of multiple genes associated with tumorigenesis, such as ATG12 and XRCC4. Within the human population, individual LTR10 elements exhibit germline and somatic structural variation resulting from a highly mutable internal tandem repeat region, which affects AP1 binding activity. Our findings reveal that ERV-derived enhancers contribute to transcriptional dysregulation in response to oncogenic signaling and shape the evolution of cancer-specific regulatory networks.

Spontaneous single-nucleotide substitutions and microsatellite mutations have distinct distributions of fitness effects

The fitness effects of new mutations determine key properties of evolutionary processes. Beneficial mutations drive evolution, yet selection is also shaped by the frequency of small-effect deleterious mutations, whose combined effect can burden otherwise adaptive lineages and alter evolutionary trajectories and outcomes in clonally evolving organisms such as viruses, microbes, and tumors. The small effect sizes of these important mutations have made accurate measurements of their rates difficult. In microbes, assessing the effect of mutations on growth can be especially instructive, as this complex phenotype is closely linked to fitness in clonally evolving organisms. Here, we perform high-throughput time-lapse microscopy on cells from mutation-accumulation strains to precisely infer the distribution of mutational effects on growth rate in the budding yeast, Saccharomyces cerevisiae. We show that mutational effects on growth rate are overwhelmingly negative, highly skewed towards very small effect sizes, and frequent enough to suggest that deleterious hitchhikers may impose a significant burden on evolving lineages. By using lines that accumulated mutations in either wild-type or slippage repair-defective backgrounds, we further disentangle the effects of 2 common types of mutations, single-nucleotide substitutions and simple sequence repeat indels, and show that they have distinct effects on yeast growth rate. Although the average effect of a simple sequence repeat mutation is very small (approximately 0.3%), many do alter growth rate, implying that this class of frequent mutations has an important evolutionary impact.

Optical Genome Mapping for Applications in Repeat Expansion Disorders

Short tandem repeat (STR) expansions are associated with more than 60 genetic disorders. The size and stability of these expansions correlate with the severity and age of onset of the disease. Therefore, being able to accurately detect the absolute length of STRs is important. Current diagnostic assays include laborious lab experiments, including repeat-primed PCR and Southern blotting, that still cannot precisely determine the exact length of very long repeat expansions. Optical genome mapping (OGM) is a cost-effective and easy-to-use alternative to traditional cytogenetic techniques and allows the comprehensive detection of chromosomal aberrations and structural variants >500 bp in length, including insertions, deletions, duplications, inversions, translocations, and copy number variants. Here, we provide methodological guidance for preparing samples and performing OGM as well as running the analysis pipelines and using the specific repeat expansion workflows to determine the exact repeat length of repeat expansions expanded beyond 500 bp. Together these protocols provide all details needed to analyze the length and stability of any repeat expansion with an expected repeat size difference from the expected wild-type allele of >500 bp. © 2024 The Authors. Current Protocols published by Wiley Periodicals LLC. Basic Protocol 1: Genomic ultra-high-molecular-weight DNA isolation, labeling, and staining Basic Protocol 2: Data generation and genome mapping using the Bionano Saphyr® System Basic Protocol 3: Manual De Novo Assembly workflow Basic Protocol 4: Local guided assembly workflow Basic Protocol 5: EnFocus Fragile X workflow Basic Protocol 6: Molecule distance script workflow.

Spontaneous single-nucleotide substitutions and microsatellite mutations have distinct distributions of fitness effects

The fitness effects of new mutations determine key properties of evolutionary processes. Beneficial mutations drive evolution, yet selection is also shaped by the frequency of small-effect deleterious mutations, whose combined effect can burden otherwise adaptive lineages and alter evolutionary trajectories and outcomes in clonally evolving organisms such as viruses, microbes, and tumors. The small effect sizes of these important mutations have made accurate measurements of their rates difficult. In microbes, assessing the effect of mutations on growth can be especially instructive, as this complex phenotype is closely linked to fitness in clonally evolving organisms. Here, we perform high-throughput time-lapse microscopy on cells from mutation-accumulation strains to precisely infer the distribution of mutational effects on growth rate in the budding yeast, Saccharomyces cerevisiae. We show that mutational effects on growth rate are overwhelmingly negative, highly skewed towards very small effect sizes, and frequent enough to suggest that deleterious hitchhikers may impose a significant burden on evolving lineages. By using lines that accumulated mutations in either wild-type or slippage repair-defective backgrounds, we further disentangle the effects of two common types of mutations, single-nucleotide substitutions and simple sequence repeat indels, and show that they have distinct effects on yeast growth rate. Although the average effect of a simple sequence repeat mutation is very small (~0.3%), many do alter growth rate, implying that this class of frequent mutations has an important evolutionary impact.

STRchive: a dynamic resource detailing population-level and locus-specific insights at tandem repeat disease loci

Approximately 3% of the human genome consists of repetitive elements called tandem repeats (TRs), which include short tandem repeats (STRs) of 1-6bp motifs and variable number tandem repeats (VNTRs) of 7+bp motifs. TR variants contribute to several dozen mono- and polygenic diseases but remain understudied and "enigmatic," particularly relative to single nucleotide variants. It remains comparatively challenging to interpret the clinical significance of TR variants. Although existing resources provide portions of necessary data for interpretation at disease-associated loci, it is currently difficult or impossible to efficiently invoke the additional details critical to proper interpretation, such as motif pathogenicity, disease penetrance, and age of onset distributions. It is also often unclear how to apply population information to analyses. We present STRchive (S-T-archive, http://strchive.org/ ), a dynamic resource consolidating information on TR disease loci in humans from research literature, up-to-date clinical resources, and large-scale genomic databases, with the goal of streamlining TR variant interpretation at disease-associated loci. We apply STRchive -including pathogenic thresholds, motif classification, and clinical phenotypes-to a gnomAD cohort of ∼18.5k individuals genotyped at 60 disease-associated loci. Through detailed literature curation, we demonstrate that the majority of TR diseases affect children despite being thought of as adult diseases. Additionally, we show that pathogenic genotypes can be found within gnomAD which do not necessarily overlap with known disease prevalence, and leverage STRchive to interpret locus-specific findings therein. We apply a diagnostic blueprint empowered by STRchive to relevant clinical vignettes, highlighting possible pitfalls in TR variant interpretation. As a living resource, STRchive is maintained by experts, takes community contributions, and will evolve as understanding of TR diseases progresses.

Genome-wide meta-analysis of short-tandem repeats for Parkinson's disease risk using genotype imputation

Idiopathic Parkinson's disease is determined by a combination of genetic and environmental factors. Recently, the first genome-wide association study on short-tandem repeats in Parkinson's disease reported on eight suggestive short-tandem repeat-based risk loci (α = 5.3 × 10-6), of which four were novel, i.e. they had not been implicated in Parkinson's disease risk by genome-wide association analyses of single-nucleotide polymorphisms before. Here, we tested these eight candidate short-tandem repeats in a large, independent Parkinson's disease case-control dataset (n = 4757). Furthermore, we combined the results from both studies by meta-analysis resulting in the largest Parkinson's disease genome-wide association study of short-tandem repeats to date (n = 43 844). Lastly, we investigated whether leading short-tandem repeat risk variants exert functional effects on gene expression regulation based on methylation quantitative trait locus data in human 'post-mortem' brain (n = 142). None of the eight previously reported short-tandem repeats were significantly associated with Parkinson's disease in our independent dataset after multiple testing correction (α = 6.25 × 10-3). However, we observed modest support for short-tandem repeats near CCAR2 and NCOR1 in the updated meta-analyses of all available data. While the genome-wide meta-analysis did not reveal additional study-wide significant (α = 6.3 × 10-7) short-tandem repeat signals, we identified seven novel suggestive Parkinson's disease short-tandem repeat risk loci (α = 5.3 × 10-6). Of these, especially a short-tandem repeat near MEIOSIN showed consistent evidence for association across datasets. CCAR2, NCOR1 and one novel suggestive locus identified here (LINC01012) emerged from colocalization analyses showing evidence for a shared causal short-tandem repeat variant affecting both Parkinson's disease risk and cis DNA methylation in brain. Larger studies, ideally using short-tandem repeats called from whole-sequencing data, are needed to more fully investigate their role in Parkinson's disease.

A comparison of software for analysis of rare and common short tandem repeat (STR) variation using human genome sequences from clinical and population-based samples

Short tandem repeat (STR) variation is an often overlooked source of variation between genomes. STRs comprise about 3% of the human genome and are highly polymorphic. Some cause Mendelian disease, and others affect gene expression. Their contribution to common disease is not well-understood, but recent software tools designed to genotype STRs using short read sequencing data will help address this. Here, we compare software that genotypes common STRs and rarer STR expansions genome-wide, with the aim of applying them to population-scale genomes. By using the Genome-In-A-Bottle (GIAB) consortium and 1000 Genomes Project short-read sequencing data, we compare performance in terms of sequence length, depth, computing resources needed, genotyping accuracy and number of STRs genotyped. To ensure broad applicability of our findings, we also measure genotyping performance against a set of genomes from clinical samples with known STR expansions, and a set of STRs commonly used for forensic identification. We find that HipSTR, ExpansionHunter and GangSTR perform well in genotyping common STRs, including the CODIS 13 core STRs used for forensic analysis. GangSTR and ExpansionHunter outperform HipSTR for genotyping call rate and memory usage. ExpansionHunter denovo (EHdn), STRling and GangSTR outperformed STRetch for detecting expanded STRs, and EHdn and STRling used considerably less processor time compared to GangSTR. Analysis on shared genomic sequence data provided by the GIAB consortium allows future performance comparisons of new software approaches on a common set of data, facilitating comparisons and allowing researchers to choose the best software that fulfils their needs.

Association of tyrosine hydroxylase 01 (TH01) microsatellite and insulin gene (INS) variable number of tandem repeat (VNTR) with type 2 diabetes and fasting insulin secretion in Mexican population

PURPOSE

A variable number of tandem repeats (VNTR) in the insulin gene (INS) control region may be involved in type 2 diabetes (T2D). The TH01 microsatellite is near INS and may regulate it. We investigated whether the TH01 microsatellite and INS VNTR, assessed via the surrogate marker single nucleotide polymorphism rs689, are associated with T2D and serum insulin levels in a Mexican population.

METHODS

We analyzed a main case-control study (n = 1986) that used univariate and multivariate logistic regression models to calculate the risk conferred by TH01 and rs689 loci for T2D development; rs689 results were replicated in other case-control (n = 1188) and cross-sectional (n = 1914) studies.

RESULTS

TH01 alleles 6, 8, 9, and 9.3 and allele A of rs689 were independently associated with T2D, with differences between sex and age at diagnosis. TH01 alleles with ≥ 8 repeats conferred an increased risk for T2D in males compared with ≤ 7 repeats (odds ratio, ≥ 1.46; 95% confidence interval, 1.1-1.95). In females, larger alleles conferred a 1.5-fold higher risk for T2D when diagnosed ≥ 46 years but conferred protection when diagnosed ≤ 45 years. Similarly, rs689 allele A was associated with T2D in these groups. In males, larger TH01 alleles and the rs689 A allele were associated with a significant decrease in median fasting plasma insulin concentration with age in T2D cases; the reverse occurred in controls.

CONCLUSION

Larger TH01 alleles and rs689 A allele may potentiate insulin synthesis in males without T2D, a process disabled in those with T2D.

CAG repeat expansions create splicing acceptor sites and produce aberrant repeat-containing RNAs

Expansions of CAG trinucleotide repeats cause several rare neurodegenerative diseases. The disease-causing repeats are translated in multiple reading frames and without an identifiable initiation codon. The molecular mechanism of this repeat-associated non-AUG (RAN) translation is not known. We find that expanded CAG repeats create new splice acceptor sites. Splicing of proximal donors to the repeats produces unexpected repeat-containing transcripts. Upon splicing, depending on the sequences surrounding the donor, CAG repeats may become embedded in AUG-initiated open reading frames. Canonical AUG-initiated translation of these aberrant RNAs may account for proteins that have been attributed to RAN translation. Disruption of the relevant splice donors or the in-frame AUG initiation codons is sufficient to abrogate RAN translation. Our findings provide a molecular explanation for the abnormal translation products observed in CAG trinucleotide repeat expansion disorders and add to the repertoire of mechanisms by which repeat expansion mutations disrupt cellular functions.

CRISPR/Cas9-mediated deletion of a GA-repeat in human GPM6B leads to disruption of neural cell differentiation from NT2 cells

The human neuron-specific gene, GPM6B (Glycoprotein membrane 6B), is considered a key gene in neural cell functionality. This gene contains an exceptionally long and strictly monomorphic short tandem repeat (STR) of 9-repeats, (GA)9. STRs in regulatory regions, may impact on the expression of nearby genes. We used CRISPR-based tool to delete this GA-repeat in NT2 cells, and analyzed the consequence of this deletion on GPM6B expression. Subsequently, the edited cells were induced to differentiate into neural cells, using retinoic acid (RA) treatment. Deletion of the GA-repeat significantly decreased the expression of GPM6B at the RNA (p < 0.05) and protein (40%) levels. Compared to the control cells, the edited cells showed dramatic decrease of the astrocyte and neural cell markers, including GFAP (0.77-fold), TUBB3 (0.57-fold), and MAP2 (0.2-fold). Subsequent sorting of the edited cells showed an increased number of NES (p < 0.01), but a decreased number of GFAP (p < 0.001), TUBB3 (p < 0.05), and MAP2 (p < 0.01), compared to the control cells. In conclusion, CRISPR/Cas9-mediated deletion of a GA-repeat in human GPM6B, led to decreased expression of this gene, which in turn, disrupted differentiation of NT2 cells into neural cells.

Modification of Huntington's disease by short tandem repeats

Expansions of glutamine-coding CAG trinucleotide repeats cause a number of neurodegenerative diseases, including Huntington's disease and several of spinocerebellar ataxias. In general, age-at-onset of the polyglutamine diseases is inversely correlated with the size of the respective inherited expanded CAG repeat. Expanded CAG repeats are also somatically unstable in certain tissues, and age-at-onset of Huntington's disease corrected for individual HTT CAG repeat length (i.e. residual age-at-onset), is modified by repeat instability-related DNA maintenance/repair genes as demonstrated by recent genome-wide association studies. Modification of one polyglutamine disease (e.g. Huntington's disease) by the repeat length of another (e.g. ATXN3, CAG expansions in which cause spinocerebellar ataxia 3) has also been hypothesized. Consequently, we determined whether age-at-onset in Huntington's disease is modified by the CAG repeats of other polyglutamine disease genes. We found that the CAG measured repeat sizes of other polyglutamine disease genes that were polymorphic in Huntington's disease participants but did not influence Huntington's disease age-at-onset. Additional analysis focusing specifically on ATXN3 in a larger sample set (n = 1388) confirmed the lack of association between Huntington's disease residual age-at-onset and ATXN3 CAG repeat length. Additionally, neither our Huntington's disease onset modifier genome-wide association studies single nucleotide polymorphism data nor imputed short tandem repeat data supported the involvement of other polyglutamine disease genes in modifying Huntington's disease. By contrast, our genome-wide association studies based on imputed short tandem repeats revealed significant modification signals for other genomic regions. Together, our short tandem repeat genome-wide association studies show that modification of Huntington's disease is associated with short tandem repeats that do not involve other polyglutamine disease-causing genes, refining the landscape of Huntington's disease modification and highlighting the importance of rigorous data analysis, especially in genetic studies testing candidate modifiers.

Rediscovering tandem repeat variation in schizophrenia: challenges and opportunities

Tandem repeats (TRs) are prevalent throughout the genome, constituting at least 3% of the genome, and often highly polymorphic. The high mutation rate of TRs, which can be orders of magnitude higher than single-nucleotide polymorphisms and indels, indicates that they are likely to make significant contributions to phenotypic variation, yet their contribution to schizophrenia has been largely ignored by recent genome-wide association studies (GWAS). Tandem repeat expansions are already known causative factors for over 50 disorders, while common tandem repeat variation is increasingly being identified as significantly associated with complex disease and gene regulation. The current review summarizes key background concepts of tandem repeat variation as pertains to disease risk, elucidating their potential for schizophrenia association. An overview of next-generation sequencing-based methods that may be applied for TR genome-wide identification is provided, and some key methodological challenges in TR analyses are delineated.

Genetic modifiers of repeat expansion disorders

Repeat expansion disorders (REDs) are monogenic diseases caused by a sequence of repetitive DNA expanding above a pathogenic threshold. A common feature of the REDs is a strong genotype-phenotype correlation in which a major determinant of age at onset (AAO) and disease progression is the length of the inherited repeat tract. Over a disease-gene carrier's life, the length of the repeat can expand in somatic cells, through the process of somatic expansion which is hypothesised to drive disease progression. Despite being monogenic, individual REDs are phenotypically variable, and exploring what genetic modifying factors drive this phenotypic variability has illuminated key pathogenic mechanisms that are common to this group of diseases. Disease phenotypes are affected by the cognate gene in which the expansion is found, the location of the repeat sequence in coding or non-coding regions and by the presence of repeat sequence interruptions. Human genetic data, mouse models and in vitro models have implicated the disease-modifying effect of DNA repair pathways via the mechanisms of somatic mutation of the repeat tract. As such, developing an understanding of these pathways in the context of expanded repeats could lead to future disease-modifying therapies for REDs.

Expanding horizons of tandem repeats in biology and medicine: Why 'genomic dark matter' matters

Approximately half of the human genome includes repetitive sequences, and these DNA sequences (as well as their transcribed repetitive RNA and translated amino-acid repeat sequences) are known as the repeatome. Within this repeatome there are a couple of million tandem repeats, dispersed throughout the genome. These tandem repeats have been estimated to constitute ∼8% of the entire human genome. These tandem repeats can be located throughout exons, introns and intergenic regions, thus potentially affecting the structure and function of tandemly repetitive DNA, RNA and protein sequences. Over more than three decades, more than 60 monogenic human disorders have been found to be caused by tandem-repeat mutations. These monogenic tandem-repeat disorders include Huntington's disease, a variety of ataxias, amyotrophic lateral sclerosis and frontotemporal dementia, as well as many other neurodegenerative diseases. Furthermore, tandem-repeat disorders can include fragile X syndrome, related fragile X disorders, as well as other neurological and psychiatric disorders. However, these monogenic tandem-repeat disorders, which were discovered via their dominant or recessive modes of inheritance, may represent the 'tip of the iceberg' with respect to tandem-repeat contributions to human disorders. A previous proposal that tandem repeats may contribute to the 'missing heritability' of various common polygenic human disorders has recently been supported by a variety of new evidence. This includes genome-wide studies that associate tandem-repeat mutations with autism, schizophrenia, Parkinson's disease and various types of cancers. In this article, I will discuss how tandem-repeat mutations and polymorphisms could contribute to a wide range of common disorders, along with some of the many major challenges of tandem-repeat biology and medicine. Finally, I will discuss the potential of tandem repeats to be therapeutically targeted, so as to prevent and treat an expanding range of human disorders.

Structural variants and short tandem repeats impact gene expression and splicing in bovine testis tissue

Structural variants (SVs) and short tandem repeats (STRs) are significant sources of genetic variation. However, the impacts of these variants on gene regulation have not been investigated in cattle. Here, we genotyped and characterized 19,408 SVs and 374,821 STRs in 183 bovine genomes and investigated their impact on molecular phenotypes derived from testis transcriptomes. We found that 71% STRs were multiallelic. The vast majority (95%) of STRs and SVs were in intergenic and intronic regions. Only 37% SVs and 40% STRs were in high linkage disequilibrium (LD) (R2 > 0.8) with surrounding SNPs/insertions and deletions (Indels), indicating that SNP-based association testing and genomic prediction are blind to a nonnegligible portion of genetic variation. We showed that both SVs and STRs were more than 2-fold enriched among expression and splicing QTL (e/sQTL) relative to SNPs/Indels and were often associated with differential expression and splicing of multiple genes. Deletions and duplications had larger impacts on splicing and expression than any other type of SV. Exonic duplications predominantly increased gene expression either through alternative splicing or other mechanisms, whereas expression- and splicing-associated STRs primarily resided in intronic regions and exhibited bimodal effects on the molecular phenotypes investigated. Most e/sQTL resided within 100 kb of the affected genes or splicing junctions. We pinpoint candidate causal STRs and SVs associated with the expression of SLC13A4 and TTC7B and alternative splicing of a lncRNA and CAPP1. We provide a catalog of STRs and SVs for taurine cattle and show that these variants contribute substantially to gene expression and splicing variation.

Microsatellites used in forensics are in regions enriched for trait-associated variants

The 20 short tandem repeat (STR) loci of the combined DNA index system (CODIS) are the basis of the vast majority of forensic genetics in the United States. One argument for permissive rules about the collection of CODIS genotypes is that the CODIS loci are thought to contain little information about ancestry or traits. However, in the past 20 years, a growing field has identified hundreds of thousands of genotype-trait associations. Here, we conduct a survey of the landscape of such associations surrounding the CODIS loci as compared with non-CODIS STRs. Although this study cannot establish or quantify associations between CODIS genotypes and phenotypes, we find that the regions around the CODIS loci are enriched for both known pathogenic variants (> 90th percentile) and for trait-associated SNPs identified in genome-wide association studies (GWAS) (≥ 95th percentile in 10kb and 100kb flanking regions), compared with other random sets of autosomal tetranucleotide-repeat STRs.

CAG repeat expansions create splicing acceptor sites and produce aberrant repeat-containing RNAs

Expansions of CAG trinucleotide repeats cause several rare neurodegenerative diseases. The disease-causing repeats are translated in multiple reading frames, without an identifiable initiation codon. The molecular mechanism of this repeat-associated non-AUG (RAN) translation is not known. We find that expanded CAG repeats create new splice acceptor sites. Splicing of proximal donors to the repeats produces unexpected repeat-containing transcripts. Upon splicing, depending on the sequences surrounding the donor, CAG repeats may become embedded in AUG-initiated open reading frames. Canonical AUG-initiated translation of these aberrant RNAs accounts for proteins that are attributed to RAN translation. Disruption of the relevant splice donors or the in-frame AUG initiation codons is sufficient to abrogate RAN translation. Our findings provide a molecular explanation for the abnormal translation products observed in CAG trinucleotide repeat expansion disorders and add to the repertoire of mechanisms by which repeat expansion mutations disrupt cellular functions.

VarSCAT: A computational tool for sequence context annotations of genomic variants

The sequence contexts of genomic variants play important roles in understanding biological significances of variants and potential sequencing related variant calling issues. However, methods for assessing the diverse sequence contexts of genomic variants such as tandem repeats and unambiguous annotations have been limited. Herein, we describe the Variant Sequence Context Annotation Tool (VarSCAT) for annotating the sequence contexts of genomic variants, including breakpoint ambiguities, flanking bases of variants, wildtype/mutated DNA sequences, variant nomenclatures, distances between adjacent variants, tandem repeat regions, and custom annotation with user customizable options. Our analyses demonstrate that VarSCAT is more versatile and customizable than the currently available methods or strategies for annotating variants in short tandem repeat (STR) regions or insertions and deletions (indels) with breakpoint ambiguity. Variant sequence context annotations of high-confidence human variant sets with VarSCAT revealed that more than 75% of all human individual germline and clinically relevant indels have breakpoint ambiguities. Moreover, we illustrate that more than 80% of human individual germline small variants in STR regions are indels and that the sizes of these indels correlated with STR motif sizes. VarSCAT is available from https://github.com/elolab/VarSCAT.

Short Tandem Repeats of Human Genome Are Intrinsically Unstable in Cultured Cells in vivo

Short tandem repeats (STRs) are a class of abundant structural or functional elements in the human genome and exhibit a polymorphic nature of repeat length and genetic variation within human populations. Interestingly, STR expansions underlie about 60 neurological disorders. However, "stutter" artifacts or noises render it difficult to investigate the pathogenesis of STR expansions. Here, we systematically investigated STR instability in cultured human cells using GC-rich CAG and AT-rich ATTCT tandem repeats as examples. We found that triplicate bidirectional Sanger sequencing with PCR amplification under proper conditions can reliably assess STR length. In addition, we found that next-generation sequencing with paired-end reads bidirectionally covering STR regions can accurately and reliably assay STR length. Finally, we found that STRs are intrinsically unstable in cultured human cell populations and during single-cell cloning. Our data suggest a general method for accurately and reliably assessing STR length and have important implications in investigating pathogenesis of STR expansion diseases.

Characterization of genome-wide STR variation in 6487 human genomes

Short tandem repeats (STRs) are abundant and highly mutagenic in the human genome. Many STR loci have been associated with a range of human genetic disorders. However, most population-scale studies on STR variation in humans have focused on European ancestry cohorts or are limited by sequencing depth. Here, we depicted a comprehensive map of 366,013 polymorphic STRs (pSTRs) constructed from 6487 deeply sequenced genomes, comprising 3983 Chinese samples (~31.5x, NyuWa) and 2504 samples from the 1000 Genomes Project (~33.3x, 1KGP). We found that STR mutations were affected by motif length, chromosome context and epigenetic features. We identified 3273 and 1117 pSTRs whose repeat numbers were associated with gene expression and 3'UTR alternative polyadenylation, respectively. We also implemented population analysis, investigated population differentiated signatures, and genotyped 60 known disease-causing STRs. Overall, this study further extends the scale of STR variation in humans and propels our understanding of the semantics of STRs.

Interplay Between Polymorphic Short Tandem Repeats and Gene Expression Variation in Caenorhabditis elegans

Short tandem repeats (STRs) have orders of magnitude higher mutation rates than single nucleotide variants (SNVs) and have been proposed to accelerate evolution in many organisms. However, only few studies have addressed the impact of STR variation on phenotypic variation at both the organismal and molecular levels. Potential driving forces underlying the high mutation rates of STRs also remain largely unknown. Here, we leverage the recently generated expression and STR variation data among wild Caenorhabditis elegans strains to conduct a genome-wide analysis of how STRs affect gene expression variation. We identify thousands of expression STRs (eSTRs) showing regulatory effects and demonstrate that they explain missing heritability beyond SNV-based expression quantitative trait loci. We illustrate specific regulatory mechanisms such as how eSTRs affect splicing sites and alternative splicing efficiency. We also show that differential expression of antioxidant genes and oxidative stresses might affect STR mutations systematically using both wild strains and mutation accumulation lines. Overall, we reveal the interplay between STRs and gene expression variation by providing novel insights into regulatory mechanisms of STRs and highlighting that oxidative stress could lead to higher STR mutation rates.

Native functions of short tandem repeats

Over a third of the human genome is comprised of repetitive sequences, including more than a million short tandem repeats (STRs). While studies of the pathologic consequences of repeat expansions that cause syndromic human diseases are extensive, the potential native functions of STRs are often ignored. Here, we summarize a growing body of research into the normal biological functions for repetitive elements across the genome, with a particular focus on the roles of STRs in regulating gene expression. We propose reconceptualizing the pathogenic consequences of repeat expansions as aberrancies in normal gene regulation. From this altered viewpoint, we predict that future work will reveal broader roles for STRs in neuronal function and as risk alleles for more common human neurological diseases.

RPTRF: A rapid perfect tandem repeat finder tool for DNA sequences

The sequencing of eukaryotic genomes has shown that tandem repeats are abundant in their sequences. In addition to affecting some cellular processes, tandem repeats in the genome may be associated with specific diseases and have been the key to resolving criminal cases. Any tool developed for detecting tandem repeats must be accurate, fast, and useable in thousands of laboratories worldwide, including those with not very advanced computing capabilities. The proposed method, the Rapid Perfect Tandem Repeat Finder (RPTRF), minimizes the need for excess character comparison processing by indexing the input file and significantly helps to accelerate and prepare the output without artifacts by using an interval tree in the filtering section. The experiments demonstrated that the RPTRF is very fast in discovering all perfect tandem repeats of all categories of any genomic sequences. Although the detection of imperfect TRs is not the focus of the RPTRF, comparisons show that it even outperforms some other tools (in five selected gold standards) designed explicitly for this purpose. The implemented tool and how to use it are available on GitHub.

CRISPR-Cas9 Targeted Enrichment and Next-Generation Sequencing for Mutation Detection

Despite the rapid application of next-generation sequencing (NGS) technologies, target sequencing in regions of the genome is often required to diagnose many genetic diseases. Target enrichment can be an effective factor in reducing the cost of sequencing and the duration of sequencing. Recently, several clustered system regularly interspaced short palindromic repeats (CRISPR)-based methods (amplification-free sequencing) have been developed to target enrichment in combination with one of the NGS platforms. CRISPR-based target enrichment strategies act as an auxiliary tool to improve NGS analytical performance, thereby indirectly facilitating nucleic acid detection. The direct DNA cleavage approach by CRISPR-Cas at genome-specific sites enhances the possibility of separating native large fragments from disease-related genomic regions. The CRISPR-Cas can isolate the target region without any amplification; subsequently, long-read sequencing technologies were also implemented. These methods, as promising tools, have the ability to assess genetic and epigenetic composition for clinical application and treatment responses in cancer precision medicine. By modifying CRISPR-based enrichment protocols, it was possible to identify different types of mutations, including structural variants, short tandem repeats, fusion genes, and mobile elements. The Cas9 can specifically eliminate wild-type sequences, and it also enables the enrichment and detection of small amounts of tumor DNA fragments among the highly heterogeneous fragments of wild-type DNA.

The role of structural variations in Alzheimer's disease and other neurodegenerative diseases

Dozens of single nucleotide polymorphisms (SNPs) related to Alzheimer's disease (AD) have been discovered by large scale genome-wide association studies (GWASs). However, only a small portion of the genetic component of AD can be explained by SNPs observed from GWAS. Structural variation (SV) can be a major contributor to the missing heritability of AD; while SV in AD remains largely unexplored as the accurate detection of SVs from the widely used array-based and short-read technology are still far from perfect. Here, we briefly summarized the strengths and weaknesses of available SV detection methods. We reviewed the current landscape of SV analysis in AD and SVs that have been found associated with AD. Particularly, the importance of currently less explored SVs, including insertions, inversions, short tandem repeats, and transposable elements in neurodegenerative diseases were highlighted.

Rare tandem repeat expansions associate with genes involved in synaptic and neuronal signaling functions in schizophrenia

Tandem repeat expansions (TREs) are associated with over 60 monogenic disorders and have recently been implicated in complex disorders such as cancer and autism spectrum disorder. The role of TREs in schizophrenia is now emerging. In this study, we have performed a genome-wide investigation of TREs in schizophrenia. Using genome sequence data from 1154 Swedish schizophrenia cases and 934 ancestry-matched population controls, we have detected genome-wide rare (<0.1% population frequency) TREs that have motifs with a length of 2-20 base pairs. We find that the proportion of individuals carrying rare TREs is significantly higher in the schizophrenia group. There is a significantly higher burden of rare TREs in schizophrenia cases than in controls in genic regions, particularly in postsynaptic genes, in genes overlapping brain expression quantitative trait loci, and in brain-expressed genes that are differentially expressed between schizophrenia cases and controls. We demonstrate that TRE-associated genes are more constrained and primarily impact synaptic and neuronal signaling functions. These results have been replicated in an independent Canadian sample that consisted of 252 schizophrenia cases of European ancestry and 222 ancestry-matched controls. Our results support the involvement of rare TREs in schizophrenia etiology.

USAT: a bioinformatic toolkit to facilitate interpretation and comparative visualization of tandem repeat sequences

Background

Tandem repeats (TR), highly variable genomic variants, are widely used in individual identification, disease diagnostics, and evolutionary studies. The recent advances in sequencing technologies and bioinformatic tools facilitate calling TR haplotypes genome widely. Both length-based and sequence-based TR alleles are used in different applications. However, sequence-based TR alleles could provide the highest precision in characterizing TR haplotypes. The need to identify the differences at the single nucleotide level between or among TR haplotypes with an easy-use bioinformatic tool is essential.

Results

In this study, we developed a Universal STR Allele Toolkit (USAT) for TR haplotype analysis, which takes TR haplotype output from existing tools to perform allele size conversion, sequence comparison of haplotypes, figure plotting, comparison for allele distribution, and interactive visualization. An exemplary application of USAT for analysis of the CODIS core STR loci for DNA forensics with benchmarking human individuals demonstrated the capabilities of USAT. USAT has user-friendly graphic interfaces and runs fast in major computing operating systems with parallel computing enabled.

Conclusion

USAT is a user-friendly bioinformatics software for interpretation, visualization, and comparisons of TRs.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12859-022-05021-1.

Replication dependent and independent mechanisms of GAA repeat instability

Trinucleotide repeat instability is a driver of human disease. Large expansions of (GAA)n repeats in the first intron of the FXN gene are the cause Friedreich's ataxia (FRDA), a progressive degenerative disorder which cannot yet be prevented or treated. (GAA)n repeat instability arises during both replication-dependent processes, such as cell division and intergenerational transmission, as well as in terminally differentiated somatic tissues. Here, we provide a brief historical overview on the discovery of (GAA)n repeat expansions and their association to FRDA, followed by recent advances in the identification of triplex H-DNA formation and replication fork stalling. The main body of this review focuses on the last decade of progress in understanding the mechanism of (GAA)n repeat instability during DNA replication and/or DNA repair. We propose that the discovery of additional mechanisms of (GAA)n repeat instability can be achieved via both comparative approaches to other repeat expansion diseases and genome-wide association studies. Finally, we discuss the advances towards FRDA prevention or amelioration that specifically target (GAA)n repeat expansions.

Establishment of a co-analysis system for personal identification and body fluid identification: a preliminary report

Analysis of genetic markers can provide clues for case investigation. Short tandem repeat (STR) detection and analysis are widely used for both personal identification and parentage testing. However, DNA analysis currently cannot provide sufficient information for body fluid identification. Tissue or cell sources of samples can be identified by detecting body fluid-specific mRNA markers, which have been studied thoroughly. Integrating STR profiling and mRNA expression patterns can provide more information than conventional methods for investigations and the reconstruction of crime scenes; this can be achieved by DNA/RNA co-extraction technology, which is economical, efficient, and suitable for low-template samples. Here, we propose a co-analysis system based on the PowerPlex 16 kit. This system can simultaneously amplify 25 markers, including 15 STRs, one non-STR amelogenin, and nine mRNA markers (three blood-specific, two saliva-specific, two semen-specific, and two housekeeping gene markers). The specificity and sensitivity of the co-analysis system were determined and aged and degraded samples were used to validate the stability of the co-analysis system. Finally, different DNA/RNA ratios and various carriers were evaluated. The results showed that the DNA/RNA co-analysis system correctly identified different types of body fluid stains. The STR profiles obtained using the co-analysis system were identical to those obtained using the PP16 kit, which demonstrates that the mRNA primers used did not affect STR profiling. Complete STR and mRNA profiles could be obtained from 1/8 portions of buccal swabs, 1/16 portions of swabs of blood and semen samples, 0.1 cm² of blood samples, 0.25 cm² of semen samples, and 1.0 cm² saliva samples. Additionally, our findings indicate that complete STR and mRNA profiles can be obtained with this system from blood and semen samples when the DNA/RNA ratio is 1:1/32. This study suggests that the co-analysis system could be used for simultaneous personal identification and body fluid identification.

The number of mitochondrial DNA mutations as a genetic feature for hair cycle arrest (alopecia X) in Pomeranian dogs

BACKGROUND

Hair cycle arrest (HCA) is a noninflammatory alopecic disease affecting various dog breeds, particularly Pomeranian dogs. This disease is probably a hereditary disorder considering the strong breed predisposition. Despite efforts to identify the pathogenesis of this disease, an underlying specific cause is unknown.

OBJECTIVE

To identify candidate gene mutations for HCA in Pomeranian dogs.

ANIMALS

Four Pomeranian dogs diagnosed with HCA and four unaffected Pomeranian dogs.

MATERIALS AND METHODS

Whole blood was used for DNA extraction. Whole-genome sequencing (WGS) was performed, and variants were analysed using a Genome Analysis Toolkit (GATK) and SnpEff. All reads were aligned to the reference genome, Dog10K_Boxer_Tasha. Sanger sequencing was performed to define the complex mutations.

RESULTS

A total of 113 variants of mitochondrial DNA were found to be effective gene mutations in the eight dogs. The affected dogs showed significantly increased effective mutations (average 57 variants) compared with unaffected dogs (average eight variants; p < 0.05). There was no significant difference in the number of chromosomal DNA mutations between the two groups.

CONCLUSION AND CLINICAL IMPORTANCE

We suggest that an increased number of mitochondrial gene mutations are features for HCA in Pomeranian dogs.

STRipy: A graphical application for enhanced genotyping of pathogenic short tandem repeats in sequencing data

Expansions of short tandem repeats (STRs) have been implicated as the causal variant in over 50 diseases known to date. There are several tools which can genotype STRs from high-throughput sequencing (HTS) data. However, running these tools out of the box only allows around half of the known disease-causing loci to be genotyped. Furthermore, the genotypes estimated at these loci are often underestimated with maximum lengths limited to either the read or fragment length, which is less than the pathogenic cutoff for some diseases. Although analysis tools can be customized to genotype extra loci, this requires proficiency in bioinformatics to set up, limiting their widespread usage by other researchers and clinicians. To address these issues, we have developed a new software called STRipy, which is able to target all known disease-causing STRs from HTS data. We created an intuitive graphical interface for STRipy and significantly simplified the detection of STRs expansions. Moreover, we genotyped all disease loci for over two and half thousand samples to provide population-wide distributions to assist with interpretation of results. We believe the simplicity and breadth of STRipy will increase the genotyping of STRs in sequencing data resulting in further diagnoses of rare STR diseases.

Forensic Efficiency Estimation of a Homemade Six-Color Fluorescence Multiplex Panel and In-Depth Anatomy of the Population Genetic Architecture in Two Tibetan Groups

The genetic information of the Chinese Tibetan group has been a long-standing research hotspot among population geneticists and archaeologists. Herein, 309 unrelated individuals from two Tibetan groups living in Qinghai Province, China (CTQ), and Tibet Autonomous Region, China (CTT), were successfully genotyped using a new homemade six-color fluorescence multiplex panel, which contained 59 autosomal deletion/insertion polymorphisms (au-DIPs), two mini short tandem repeats (miniSTRs), two Y-chromosomal DIPs, and one Amelogenin. The cumulative probability of matching and combined power of exclusion values for this new panel in CTQ and CTT groups were 1.9253E-27 and 0.99999729, as well as 1.5061E-26 and 0.99999895, respectively. Subsequently, comprehensive population genetic analyses of Tibetan groups and reference populations were carried out based on the 59 au-DIPs. The multitudinous statistical analysis results supported that Tibetan groups have close genetic affinities with East Asian populations. These findings showed that this homemade system would be a powerful tool for forensic individual identification and paternity testing in Chinese Tibetan groups and give us an important insight for further perfecting the genetic landscape of Tibetan groups.

Linked-read sequencing for detecting short tandem repeat expansions

Detection of short tandem repeat (STR) expansions with standard short-read sequencing is challenging due to the difficulty in mapping multicopy repeat sequences. In this study, we explored how the long-range sequence information of barcode linked-read sequencing (BLRS) can be leveraged to improve repeat-read detection. We also devised a novel algorithm using BLRS barcodes for distance estimation and evaluated its application for STR genotyping. Both approaches were designed for genotyping large expansions (> 1 kb) that cannot be sized accurately by existing methods. Using simulated and experimental data of genomes with STR expansions from multiple BLRS platforms, we validated the utility of barcode and phasing information in attaining better STR genotypes compared to standard short-read sequencing. Although the coverage bias of extremely GC-rich STRs is an important limitation of BLRS, BLRS is an effective strategy for genotyping many other STR loci.

Identification and Monitoring of Nucleotide Repeat Expansions Using Southern Blotting in Drosophila Models of C9orf72 Motor Neuron Disease and Frontotemporal Dementia

Repeat expansion diseases, including fragile X syndrome, Huntington's disease, and C9orf72-related motor neuron disease and frontotemporal dementia, are a group of disorders associated with polymorphic expansions of tandem repeat nucleotide sequences. These expansions are highly repetitive and often hundreds to thousands of repeats in length, making accurate identification and determination of repeat length via PCR or sequencing challenging. Here we describe a protocol for monitoring repeat length in Drosophila models carrying 1,000 repeat C9orf72-related dipeptide repeat transgenes using Southern blotting. This protocol has been used regularly to check the length of these lines for over 100 generations with robust and repeatable results and can be implemented for monitoring any repeat expansion in Drosophila.

Profiling the Genome-Wide Landscape of Short Tandem Repeats by Long-Read Sequencing

Background: Short tandem repeats (STRs) are highly variable elements that play a pivotal role in multiple genetic diseases and the regulation of gene expression. Long-read sequencing (LRS) offers a potential solution to genome-wide STR analysis. However, characterizing STRs in human genomes using LRS on a large population scale has not been reported.

Methods: We conducted the large LRS-based STR analysis in 193 unrelated samples of the Chinese population and performed genome-wide profiling of STR variation in the human genome. The repeat dynamic index (RDI) was introduced to evaluate the variability of STR. We sourced the expression data from the Genotype-Tissue Expression to explore the tissue specificity of highly variable STRs related genes across tissues. Enrichment analyses were also conducted to identify potential functional roles of the high variable STRs.

Results: This study reports the large-scale analysis of human STR variation by LRS and offers a reference STR database based on the LRS dataset. We found that the disease-associated STRs (dSTRs) and STRs associated with the expression of nearby genes (eSTRs) were highly variable in the general population. Moreover, tissue-specific expression analysis showed that those highly variable STRs related genes presented the highest expression level in brain tissues, and enrichment pathways analysis found those STRs are involved in synaptic function-related pathways.

Conclusion: Our study profiled the genome-wide landscape of STR using LRS and highlighted the highly variable STRs in the human genome, which provide a valuable resource for studying the role of STRs in human disease and complex traits.

Mapping short tandem repeats for liver gene expression traits helps prioritize potential causal variants for complex traits in pigs

BACKGROUND

Short tandem repeats (STRs) were recently found to have significant impacts on gene expression and diseases in humans, but their roles on gene expression and complex traits in pigs remain unexplored. This study investigates the effects of STRs on gene expression in liver tissues based on the whole-genome sequences and RNA-Seq data of a discovery cohort of 260 F6 individuals and a validation population of 296 F7 individuals from a heterogeneous population generated from crosses among eight pig breeds.

RESULTS

We identified 5203 and 5868 significantly expression STRs (eSTRs, FDR < 1%) in the F6 and F7 populations, respectively, most of which could be reciprocally validated (π1 = 0.92). The eSTRs explained 27.5% of the cis-heritability of gene expression traits on average. We further identified 235 and 298 fine-mapped STRs through the Bayesian fine-mapping approach in the F6 and F7 pigs, respectively, which were significantly enriched in intron, ATAC peak, compartment A and H3K4me3 regions. We identified 20 fine-mapped STRs located in 100 kb windows upstream and downstream of published complex trait-associated SNPs, which colocalized with epigenetic markers such as H3K27ac and ATAC peaks. These included eSTR of the CLPB, PGLS, PSMD6 and DHDH genes, which are linked with genome-wide association study (GWAS) SNPs for blood-related traits, leg conformation, growth-related traits, and meat quality traits, respectively.

CONCLUSIONS

This study provides insights into the effects of STRs on gene expression traits. The identified eSTRs are valuable resources for prioritizing causal STRs for complex traits in pigs.

Population-level genome-wide STR discovery and validation for population structure and genetic diversity assessment of Plasmodium species

Short tandem repeats (STRs) are highly informative genetic markers that have been used extensively in population genetics analysis. They are an important source of genetic diversity and can also have functional impact. Despite the availability of bioinformatic methods that permit large-scale genome-wide genotyping of STRs from whole genome sequencing data, they have not previously been applied to sequencing data from large collections of malaria parasite field samples. Here, we have genotyped STRs using HipSTR in more than 3,000 Plasmodium falciparum and 174 Plasmodium vivax published whole-genome sequence data from samples collected across the globe. High levels of noise and variability in the resultant callset necessitated the development of a novel method for quality control of STR genotype calls. A set of high-quality STR loci (6,768 from P. falciparum and 3,496 from P. vivax) were used to study Plasmodium genetic diversity, population structures and genomic signatures of selection and these were compared to genome-wide single nucleotide polymorphism (SNP) genotyping data. In addition, the genome-wide information about genetic variation and other characteristics of STRs in P. falciparum and P. vivax have been available in an interactive web-based R Shiny application PlasmoSTR (https://github.com/bahlolab/PlasmoSTR).

Revisiting tandem repeats in psychiatric disorders from perspectives of genetics, physiology, and brain evolution

Genome-wide association studies (GWASs) have revealed substantial genetic components comprised of single nucleotide polymorphisms (SNPs) in the heritable risk of psychiatric disorders. However, genetic risk factors not covered by GWAS also play pivotal roles in these illnesses. Tandem repeats, which are likely functional but frequently overlooked by GWAS, may account for an important proportion in the "missing heritability" of psychiatric disorders. Despite difficulties in characterizing and quantifying tandem repeats in the genome, studies have been carried out in an attempt to describe impact of tandem repeats on gene regulation and human phenotypes. In this review, we have introduced recent research progress regarding the genomic distribution and regulatory mechanisms of tandem repeats. We have also summarized the current knowledge of the genetic architecture and biological underpinnings of psychiatric disorders brought by studies of tandem repeats. These findings suggest that tandem repeats, in candidate psychiatric risk genes or in different levels of linkage disequilibrium (LD) with psychiatric GWAS SNPs and haplotypes, may modulate biological phenotypes related to psychiatric disorders (e.g., cognitive function and brain physiology) through regulating alternative splicing, promoter activity, enhancer activity and so on. In addition, many tandem repeats undergo tight natural selection in the human lineage, and likely exert crucial roles in human brain evolution. Taken together, the putative roles of tandem repeats in the pathogenesis of psychiatric disorders is strongly implicated, and using examples from previous literatures, we wish to call for further attention to tandem repeats in the post-GWAS era of psychiatric disorders.

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.

Molecular mechanisms underlying nucleotide repeat expansion disorders

The human genome contains over one million short tandem repeats. Expansion of a subset of these repeat tracts underlies over fifty human disorders, including common genetic causes of amyotrophic lateral sclerosis (ALS) and frontotemporal dementia (C9orf72), polyglutamine-associated ataxias and Huntington disease, myotonic dystrophy, and intellectual disability disorders such as Fragile X syndrome. In this Review, we discuss the four major mechanisms by which expansion of short tandem repeats causes disease: loss of function through transcription repression, RNA-mediated gain of function through gelation and sequestration of RNA-binding proteins, gain of function of canonically translated repeat-harbouring proteins, and repeat-associated non-AUG translation of toxic repeat peptides. Somatic repeat instability amplifies these mechanisms and influences both disease age of onset and tissue specificity of pathogenic features. We focus on the crosstalk between these disease mechanisms, and argue that they often synergize to drive pathogenesis. We also discuss the emerging native functions of repeat elements and how their dynamics might contribute to disease at a larger scale than currently appreciated. Lastly, we propose that lynchpins tying these disease mechanisms and native functions together offer promising therapeutic targets with potential shared applications across this class of human disorders.