Link between short tandem repeats and translation initiation site selection

A Hypermutable Region in the DISP2 Gene Links to Natural Selection and Late-Onset Neurocognitive Disorders in Humans

(CCG) short tandem repeats (STRs) are predominantly enriched in genic regions, mutation hotspots for C to T truncating substitutions, and involved in various neurological and neurodevelopmental disorders. However, intact blocks of this class of STRs are widely overlooked with respect to their link with natural selection. The human neuron-specific gene, DISP2 (dispatched RND transporter family member 2), contains a (CCG) repeat in its 5' untranslated region. Here, we sequenced this STR in a sample of 448 Iranian individuals, consisting of late-onset neurocognitive disorder (NCD) (N = 203) and controls (N = 245). We found that the region spanning the (CCG) repeat was highly mutated, resulting in several flanking (CCG) residues. However, an 8-repeat of the (CCG) repeat was predominantly abundant (frequency = 0.92) across the two groups. While the overall distribution of genotypes was not different between the two groups (p > 0.05), we detected four genotypes in the NCD group only (2% of the NCD genotypes, Mid-p = 0.02), consisting of extreme short alleles, 5- and 6-repeats, that were not detected in the control group. The patients harboring those genotypes received the diagnoses of probable Alzheimer's disease and vascular dementia. We also found six genotypes in the control group only (2.5% of the control genotypes, Mid-p = 0.01) that consisted of the 8-repeat and extreme long alleles, 9- and 10-repeats, of which the 10-repeat was not detected in the NCD group. The (CCG) repeat specifically expanded in primates. In conclusion, we report an indication of natural selection at a novel hypermutable region in the human genome and divergent alleles and genotypes in late-onset NhCDs and controls. These findings reinforce the hypothesis that a collection of rare alleles and genotypes in a number of genes may unambiguously contribute to the cognition impairment component of late-onset NCDs.

The human SMAD9 (GCC) repeat links to natural selection and late-onset neurocognitive disorders

INTRODUCTION

Whereas (GCC)-repeats are overrepresented in genic regions, and mutation hotspots, they are largely unexplored with regard to their link with natural selection. Across numerous primate species and tissues, SMAD9 (SMAD Family Member 9) reaches highest level of expression in the human brain. This gene contains a (GCC)-repeat in the interval between + 1 and + 60 of the transcription start site, which is in the high-ranking (GCC)-repeats with respect to length.

METHODS

Here we sequenced this (GCC)-repeat in 396 Iranian individuals, consisting of late-onset neurocognitive disorder (NCD) (N = 181) and controls (N = 215).

RESULTS

We detected two predominantly abundant alleles of 7 and 9 repeats, forming 96.2% of the allele pool. The (GCC)7/(GCC)9 ratio was in the reverse order in the NCD group versus controls (p = 0.005), resulting from excess of (GCC)7 in the NCD group (p = 0.003) and (GCC)9 in the controls (p = 0.01). Five genotypes, predominantly consisting of (GCC)7 and lacking (GCC)9 were detected in the NCD group only (p = 0.008). The patients harboring those genotypes received the diagnoses of Alzheimer's disease (AD) and vascular dementia (VD). Five genotypes consisting of (GCC)9 and lacking (GCC)7 were detected in the control group only (p = 0.002). The group-specific genotypes formed approximately 4% of the genotype pool in the human samples studied.

CONCLUSION

We propose natural selection and a novel locus for late-onset AD and VD at the SMAD9 (GCC)-repeat in humans.

Novel islands of GGC and GCC repeats coincide with human evolution

BACKGROUND

Because of high mutation rate, overrepresentation in genic regions, and link with various neurological, neurodegenerative, and movement disorders, GGC and GCC short tandem repeats (STRs) are prone to natural selection. Among a number of lacking data, the 3-repeats of these STRs remain widely unexplored.

RESULTS

In a genome-wide search in human, here we mapped GGC and GCC STRs of ≥3-repeats, and found novel islands of up to 45 of those STRs, populating spans of 1 to 2 kb of genomic DNA. RGPD4 and NOC4L harbored the densest (GGC)3 (probability 3.09061E-71) and (GCC)3 (probability 1.72376E-61) islands, respectively, and were human-specific. We also found prime instances of directional incremented density of STRs at specific loci in human versus other species, including the FOXK2 and SKI GGC islands. The genes containing those islands significantly diverged in expression in human versus other species, and the proteins encoded by those genes interact closely in a physical interaction network, consequence of which may be human-specific characteristics such as higher order brain functions.

CONCLUSION

We report novel islands of GGC and GCC STRs of evolutionary relevance to human. The density, and in some instances, periodicity of these islands support them as a novel genomic entity, which need to be further explored in evolutionary, mechanistic, and functional platforms.

Dyads of GGC and GCC form hotspot colonies that coincide with the evolution of human and other great apes

BACKGROUND

GGC and GCC short tandem repeats (STRs) are of various evolutionary, biological, and pathological implications. However, the fundamental two-repeats (dyads) of these STRs are widely unexplored.

RESULTS

On a genome-wide scale, we mapped (GGC)2 and (GCC)2 dyads in human, and found monumental colonies (distance between each dyad < 500 bp) of extraordinary density, and in some instances periodicity. The largest (GCC)2 and (GGC)2 colonies were intergenic, homogeneous, and human-specific, consisting of 219 (GCC)2 on chromosome 2 (probability < 1.545E-219) and 70 (GGC)2 on chromosome 9 (probability = 1.809E-148). We also found that several colonies were shared in other great apes, and directionally increased in density and complexity in human, such as a colony of 99 (GCC)2 on chromosome 20, that specifically expanded in great apes, and reached maximum complexity in human (probability 1.545E-220). Numerous other colonies of evolutionary relevance in human were detected in other largely overlooked regions of the genome, such as chromosome Y and pseudogenes. Several of the genes containing or nearest to those colonies were divergently expressed in human.

CONCLUSION

In conclusion, (GCC)2 and (GGC)2 form unprecedented genomic colonies that coincide with the evolution of human and other great apes. The extent of the genomic rearrangements leading to those colonies support overlooked recombination hotspots, shared across great apes. The identified colonies deserve to be studied in mechanistic, evolutionary, and functional platforms.

A GCC repeat in RAB26 undergoes natural selection in human and harbors divergent genotypes in late-onset Alzheimer's disease

Although mainly located in genic regions and being mutation hotspots, intact blocks of CG-rich trinucleotide short tandem repeats (STRs) are largely overlooked with respect to their link with natural selection. The human RAB26 (member RAS oncogene family) directs synaptic and secretory vesicles into preautophagosomal structures, inhibition of which specifically disrupts axonal transport of degradative organelles and leads to an axonal dystrophy, resembling Alzheimer's disease (AD). Human RAB26 contains a GCC repeat in the top 1st percent in respect of length. Here we sequenced this STR in 441 Iranian individuals, consisting of late-onset neurocognitive disorder (NCD) (N = 216) and controls (N = 225). In both groups, the 12-repeat allele and the 12/12 genotype were predominantly abundant. We found excess of homozygosity for non-12 alleles in the NCD group (Mid-P exact = 0.027). Furthermore, divergent genotypes were detected that were specific to the NCD group (2.8% of genotypes) (Mid-P exact = 0.006) or controls (3.1% of genotypes) (Mid-P exact = 0.004). The patients harboring divergent genotypes received the diagnosis of AD. Based on the predominant abundance of the 12-repeat and 12/12 genotype in both groups, excess of non-12 homozygosity in the NCD group, and divergent genotypes across the NCD and control groups, we propose natural selection at this locus and link with late-onset AD. Our findings strengthen the hypothesis that a collection of rare genotypes unambiguously contribute to the pathogenesis of late-onset NCDs, such as AD.

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.

Tandem repeats ubiquitously flank and contribute to translation initiation sites

Background

While the evolutionary divergence of cis-regulatory sequences impacts translation initiation sites (TISs), the implication of tandem repeats (TRs) in TIS selection remains largely elusive. Here, we employed the TIS homology concept to study a possible link between TRs of all core lengths and repeats with TISs.

Methods

Human, as reference sequence, and 83 other species were selected, and data was extracted on the entire protein-coding genes (n = 1,611,368) and transcripts (n = 2,730,515) annotated for those species from Ensembl 102. Following TIS identification, two different weighing vectors were employed to assign TIS homology, and the co-occurrence pattern of TISs with the upstream flanking TRs was studied in the selected species. The results were assessed in 10-fold cross-validation.

Results

On average, every TIS was flanked by 1.19 TRs of various categories within its 120 bp upstream sequence, per species. We detected statistically significant enrichment of non-homologous human TISs co-occurring with human-specific TRs. On the contrary, homologous human TISs co-occurred significantly with non-human-specific TRs. 2991 human genes had at least one transcript, TIS of which was flanked by a human-specific TR. Text mining of a number of the identified genes, such as CACNA1A, EIF5AL1, FOXK1, GABRB2, MYH2, SLC6A8, and TTN, yielded predominant expression and functions in the human brain and/or skeletal muscle.

Conclusion

We conclude that TRs ubiquitously flank and contribute to TIS selection at the trans-species level. Future functional analyses, such as a combination of genome editing strategies and in vitro protein synthesis may be employed to further investigate the impact of TRs on TIS selection.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12863-022-01075-5.

Global abundance of short tandem repeats is non-random in rodents and primates

Background

While of predominant abundance across vertebrate genomes and significant biological implications, the relevance of short tandem repeats (STRs) (also known as microsatellites) to speciation remains largely elusive and attributed to random coincidence for the most part. Here we collected data on the whole-genome abundance of mono-, di-, and trinucleotide STRs in nine species, encompassing rodents and primates, including rat, mouse, olive baboon, gelada, macaque, gorilla, chimpanzee, bonobo, and human. The collected data were used to analyze hierarchical clustering of the STR abundances in the selected species.

Results

We found massive differential STR abundances between the rodent and primate orders. In addition, while numerous STRs had random abundance across the nine selected species, the global abundance conformed to three consistent < clusters>, as follows: <rat, mouse>, <gelada, macaque, olive baboon>, and <gorilla, chimpanzee, bonobo, human>, which coincided with the phylogenetic distances of the selected species (p < 4E-05). Exceptionally, in the trinucleotide STR compartment, human was significantly distant from all other species.

Conclusion

Based on hierarchical clustering, we propose that the global abundance of STRs is non-random in rodents and primates, and probably had a determining impact on the speciation of the two orders. We also propose the STRs and STR lengths, which predominantly conformed to the phylogeny of the selected species, exemplified by (t)10, (ct)6, and (taa4). Phylogenetic and experimental platforms are warranted to further examine the observed patterns and the biological mechanisms associated with those STRs.

Genomic instability and the link to infertility: A focus on microsatellites and genomic instability syndromes

Infertility is associated to multiple types of different genomic instabilities and is a genetic feature of genomic instability syndromes. While the mismatch repair machinery contributes to the maintenance of genome integrity, surprisingly its potential role in infertility is overlooked. Defects in mismatch repair mechanisms contribute to microsatellite instability and genomic instability syndromes, due to the inability to repair newly replicated DNA. This article reviews the literature to date to elucidate the contribution of microsatellite instability to genomic instability syndromes and infertility. The key findings presented reveal microsatellite instability is poorly researched in genomic instability syndromes and infertility.

A Structural Potential of Rare Trinucleotide Repeat Tracts in RNA

Among types of trinucleotide repeats, there is some disproportion in the frequency of their occurrence in the human exome. This research presents new data describing the folding and thermodynamic stability of short, tandem RNA repeats of 23 types, focusing on the rare, yet poorly analyzed ones. UV-melting experiments included the presence of PEG or potassium and magnesium ions to determine their effect on the stability of RNA repeats structures. Rare repeats predominantly stayed single-stranded but had the potential for base pairing with other partially complementary repeat tracts. A coexistence of suitably complementary repeat types in a single RNA creates opportunities for interaction in the context of the secondary structure of RNA. We searched the human transcriptome for model RNAs in which different, particularly rare trinucleotide repeats coexist and selected the GABRA4 and CHIC1 RNAs to study intramolecular interactions between the repeat tracts that they contain. In vitro secondary structure probing results showed that the UAA and UUG repeat tracts, present in GABRA4 3' UTR, form a double helix, which separates one of its structural domains. For the RNA CHIC1 ORF fragment containing four short AGG repeat tracts and the CGU tract, we proved the formation of quadruplexes that blocked reverse transcription.

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.

Predominant monomorphism of the RIT2 and GPM6B exceptionally long GA blocks in human and enriched divergent alleles in the disease compartment

Across human protein-coding genes, the human neuron-specific genes, RIT2 and GPM6B, contain the two longest GA short tandem repeats (STRs) of 11 and 9-repeats, respectively, the length ranges of which are functional, and result in gene expression alteration. Here we sequenced the RIT2 and GPM6B STRs in 600 human subjects, consisting of late-onset neurocognitive disorder (n = 200), multiple sclerosis (n = 200), and controls (n = 200). Furthermore, we selected two large human databases, including the general-population-based gnomAD ( https://gnomad.broadinstitute.org ) and a mainly disease-phenotype-archiving database, TOPMed ( https://www.nhlbiwgs.org ), to compare allele frequencies in the general populations vs. the disease compartment. The RIT2 and GPM6B GA-repeats were monomorphic in the human subjects studied, at lengths of 11 and 9-repeats, respectively, and were predominantly human-specific in formula. Exception included a 9/11 genotype of the RIT2 GA-STR in an isolate case of female multiple sclerosis. Exceedingly rare alleles of the two GA repeats were significantly enriched in TOPMed vs. the gnomAD. We report prime instances of predominant monomorphism for specific lengths of STRs in human, and possible enrichment of rare divergent alleles in the disease phenotype compartment. While STRs are most attended because of their high polymorphic nature, STR monomorphism is an underappreciated feature, which may have a link with natural selection and disease.

The Short Tandem Repeat of the DMT1 Gene as a Molecular Marker of Elite Long-Distance Runners

The DMT1 gene encodes divalent metal transporter 1, a membrane iron transport protein. Divalent metal transporter 1 influences cellular iron availability, which might further affect aerobic exercise capacity. Short tandem repeat (STR) polymorphisms have been used as genetic markers in the literature, yet the STR polymorphisms of the DMT1 gene have not been well studied. In this current study, we explored the polymorphisms of the DMT1 gene in a group of elite long-distance runners and controls, by using the PCR-RFLP (Restriction Fragment Length Polymorphism) and Gene scan technology. We found that the genotype frequency of the homozygous 258 bp STR polymorphism of the DMT1 gene (258 bp/258 bp) was significantly higher in the athlete group than in the controls (χ² = 14.01, p = 0.006) so does the allele frequency of the 258 bp STR polymorphism (χ² = 12.867, p = 0.008). These data suggested that the STR polymorphism of the DMT1 gene might be correlated with aerobic exercise capacity and the 258 bp homozygous (25 bp/258 bp) could be used as a molecular marker for the talent identification of elite long-distance runners.

Disease-only alleles at the extreme ends of the human ZMYM3 exceptionally long 5' UTR short tandem repeat in bipolar disorder: A pilot study

OBJECTIVE

The X-linked ZMYM3 gene (also known as ZNF261) contains the longest STR, (GA)32, identified in a human protein-coding gene 5'UTR (ENST00000373998.5: ZMYM3-207). This STR reaches maximum length in human, and is located in a complex string of four consecutive GA-STRs with a human-specific formula across the complex. A previous study in Iranian male schizophrenia (SCZ) patients revealed co-occurrence of the extreme short and long alleles of the STR with SCZ. Here we studied the allelic distribution of this STR in bipolar disorder (BD) type I. The interval encompassing the human ZMYM3 STR complex was PCR-amplified and sequenced in 546 male subjects, consisting of 157 BD patients and 389 controls.

RESULTS

We found three alleles at the extreme short (17-repeat) and long (38- and 43-repeat) ends of the allele distribution curve in the BD cases (4.4% of the BD alleles) that were not detected in the controls (Mid p < 0.0001). These alleles overlapped with the extreme disease-only alleles detected previously in the SCZ patients. Domain reconstruction of the GA-STR complex revealed significant structural alteration as a result of various sequence repeats and nucleotide compositions at the inter and intraspecies levels.

CONCLUSION

The ZMYM3 "exceptionally long" 5' UTR STR findings may alter our perspective of disease pathogenesis in psychiatric disorders, and set an example in which the low frequency alleles at the extreme short and long ends of the human STRs are, at least in part, a result of natural selection against these alleles and their unambiguous link to major human disorders.

Patterns of microsatellite distribution across eukaryotic genomes

Background

Microsatellites, or Simple Sequence Repeats (SSRs), are short tandem repeats of 1–6 nt motifs present in all genomes. Emerging evidence points to their role in cellular processes and gene regulation. Despite the huge resource of genomic information currently available, SSRs have been studied in a limited context and compared across relatively few species.

Results

We have identified ~ 685 million eukaryotic microsatellites and analyzed their genomic trends across 15 taxonomic subgroups from protists to mammals. The distribution of SSRs reveals taxon-specific variations in their exonic, intronic and intergenic densities. Our analysis reveals the differences among non-related species and novel patterns uniquely demarcating closely related species. We document several repeats common across subgroups as well as rare SSRs that are excluded almost throughout evolution. We further identify species-specific signatures in pathogens like Leishmania as well as in cereal crops, Drosophila, birds and primates. We also find that distinct SSRs preferentially exist as long repeating units in different subgroups; most unicellular organisms show no length preference for any SSR class, while many SSR motifs accumulate as long repeats in complex organisms, especially in mammals.

Conclusions

We present a comprehensive analysis of SSRs across taxa at an unprecedented scale. Our analysis indicates that the SSR composition of organisms with heterogeneous cell types is highly constrained, while simpler organisms such as protists, green algae and fungi show greater diversity in motif abundance, density and GC content. The microsatellite dataset generated in this work provides a large number of candidates for functional analysis and for studying their roles across the evolutionary landscape.

Electronic supplementary material

The online version of this article (10.1186/s12864-019-5516-5) contains supplementary material, which is available to authorized users.