- Primary research
- Open Access
- Published:
Genome-scale portrait and evolutionary significance of human-specific core promoter tri- and tetranucleotide short tandem repeats
Human Genomics volume 12, Article number: 17 (2018)
Abstract
Background
While there is an ongoing trend to identify single nucleotide substitutions (SNSs) that are linked to inter/intra-species differences and disease phenotypes, short tandem repeats (STRs)/microsatellites may be of equal (if not more) importance in the above processes. Genes that contain STRs in their promoters have higher expression divergence compared to genes with fixed or no STRs in the gene promoters. In line with the above, recent reports indicate a role of repetitive sequences in the rise of young transcription start sites (TSSs) in human evolution.
Results
Following a comparative genomics study of all human protein-coding genes annotated in the GeneCards database, here we provide a genome-scale portrait of human-specific short- and medium-size (≥ 3-repeats) tri- and tetranucleotide STRs and STR motifs in the critical core promoter region between − 120 and + 1 to the TSS and evidence of skewing of this compartment in reference to the STRs that are not human-specific (Levene’s test p < 0.001). Twenty-five percent and 26% enrichment of human-specific transcripts was detected in the tri and tetra human-specific compartments (mid-p < 0.00002 and mid-p < 0.002, respectively).
Conclusion
Our findings provide the first evidence of genome-scale skewing of STRs at a specific region of the human genome and a link between a number of these STRs and TSS selection/transcript specificity. The STRs and genes listed here may have a role in the evolution and development of characteristics and phenotypes that are unique to the human species.
Introduction
Speciation and evolution are, at least in part, due to the plasticity (expansion or contraction) of short tandem repeats (STRs)/microsatellites, which can function as “tuning knobs” in response to the environment or other genes [1,2,3]. In line with the above, certain STRs are directionally expanded in the human species or co-occur identically in related taxa such as primates [4,5,6,7,8]. Genes that contain STRs in their promoters have higher expression divergence compared to genes with fixed or no STRs in the gene promoters [9]. Recent reports indicate a role of repetitive sequences in the rise of young transcription start sites (TSSs) in human evolution [10,11,12].
Preliminary data on the sequencing of a number of “exceptionally long” STRs (≥ 6-repeats), which compose 1–2% of all human core promoter STRs [3], support critical evolutionary adaptive roles for a number of these STRs. Human specificity of the predominant allele of the RIT2 core promoter STR in the human species, the presence of the shortest allele of this STR (5-repeat) in hunter-gatherer humans (BUSHMAN KB1: rs113265205), the lack of this allele in the agricultural modern humans (Genome Aggregation database: gnomad.broadinstitute.org), and its co-occurrence with schizophrenia provide the first indication of STR allele selection in humans [13]. A link between the CYTH4 core promoter STR (the longest tetranucleotide STR identified in a human gene core promoter) with the Old World monkeys and Apes and evidence of extreme “disease-only” genotypes at this STR with schizophrenia [14] provide the first link between a primate-specific STR and higher-order brain functions in human. The “exceptionally long” CA-repeat in the core promoter of SCGB2B2 is another example of directional STR expansion in the Old World monkeys and Apes [5]. The PAXBP1 gene is an extreme example in which expansion of a core promoter CT-repeat occurs in the Old World monkeys and reaches maximum length and complexity in human; OMIM: 617621 [4].
As “exceptionally long” STRs may be subject to natural selection, short- and medium-size alleles (≥ 3-repeats) might have had similar fate. This is indicated by the predominance of specific short- and medium-size penta- and hexanucleotide STRs and their cognate transcription factors (TFs) in the critical core promoter interval [15]. Indeed, shortening of a number of STRs and their identical co-occurrence is linked to the evolution of primates [8]. In line with the above findings, repeats associated with younger human TSSs tend to be shorter than those in older TSSs [10]. In the study reported here, we present genome-scale data on two categories of STRs, i.e., tri- and tetranucleotide STRs, and their implication in human evolution.
Materials and methods
The interval between − 120 and + 1 to the TSS of all human protein-coding genes annotated in the GeneCards database (version 3.0) (www.genecards.org) was screened for tri- and tetranucleotide STRs of ≥ 3-repeats, based on the Ensembl database (versions 87-91) (asia.ensembl.org) and using the Microsatellite Repeats Finder at the following link: http://insilico.ehu.es/mini_tools/microsatellites/
The evolutionary status of the identified STRs was analyzed in 25 species (N), including primates (N = 5), non-primate mammals (N = 12), birds and reptiles (N = 5), amphibians (N = 1), and fish (N = 2), based on the Ensembl database.
Human specificity of transcripts was evaluated based on the multiple and pair-wise %identity scoring of the TSS-flanking 5′ untranslated region (UTR), using the sequence alignment program Clustal Omega (https://www.ebi.ac.uk/Tools/msa/clustalo), and the overall composition of the transcript and encoded protein (i.e., length of the transcript, number of exons and amino acids). The threshold of sequence identity was set at 50%, which was based on the comparison of two randomly selected and unrelated sequences in the human genome.
The p value for the skewing of the human-specific STR compartment was calculated using Levene’s equality of variances test.
The p values for transcript enrichment were calculated using the two by two table analysis;
the human-specific tri- and tetranucleotide STR groups were compared against corresponding randomly selected STRs from the non-human-specific STRs. The comparison was set based on the sample size of the human-specific STRs (n) and the sample size of the non-human compartments (1.5n).
Results
Overall prevalence of tri- and tetranucleotide STR motifs across human protein-coding core promoter sequences
In total, 56 and 82 STR motifs were detected for the tri- and tetranucleotide repeats, respectively (Figs. 1 and 2). The most prevalent tri- and tetranucleotide STR motifs across the human protein-coding gene core promoters were GGC and GGGC, respectively (Figs. 1 and 2). In the category of non-GC STRs, GGA and TCCC were the most prevalent tri- and tetranucleotides, respectively.
Skewing of the human-specific core promoter tri- and tetranucleotide STRs
A significant skewing of the tri- and tetranucleotide STR distribution was found in the human-specific tri- (Fig. 1) and tetranucleotide (Fig. 2) compartments (Levene’s p < 0.001). While the most prevalent tri- and tetranucleotide repeats in the non-human-specific category were the GGC- and GGGC-repeats, respectively, the most prevalent human-specific STRs were of the GCC and CTCC motifs, respectively. Disproportionate distribution of human-specific STRs was also detected in other STRs such as CCT, GAA, CTCC, GTTT, and GAAA.
The human-specific tri- and tetranucleotide STRs were of a wide range of motifs, e.g., the CCA motif in ADCY6, the TCCC motif in ARHGEF35, GCCC in DRD2, and GTTT in MCTP2 (Tables 1 and 2).
In a number of instances, not only the STR, but also the genes containing those STRs, were human-specific (e.g., ARHGEF35, AMY1C, and C1orf204). Furthermore, a number of the tri- and tetranucleotide STRs were found to be unique to the human species at the specified interval of − 120 to + 1 TSS. For example, in the tetranucleotide compartment, CACC, GACA, CCGG, GATA, TCTG, GGCT, and TTTA STRs were detected in human only.
Enrichment of human-specific transcripts at the human-specific STR compartment
Based on sequence comparison and the overall composition of the transcript and encoded protein, 25 and 26% of the transcripts in the tri and tetra human-specific compartments were found to be human-specific (mid-p < 0.00002 and mid-p < 0.002), respectively). The %identity score of multiple sequence alignment for the human-specific transcripts was 0 (exemplified in Fig. 3), and pair-wise analysis (exemplified in Fig. 4) resulted in %identity scores ranging from 37 to 48%. In the trinucleotide category, 14 genes, MPRIP, NPAS1, PAQR9, PRSS1, R3HDM2, TMEM99, ZSCAN30, C22orf24, ECSCR, AMY1C, DDX58, C1orf204, RGPD6, and LCE2B, contained human-specific transcripts. In the tetranucleotide category, five genes, DRD2, DUX4, TEAD4, ARL17B, and ARHGEF35, contained human-specific transcripts.
A number of the identified STRs were linked to noncanonical translation in the following genes, TEAD4, ECSCR, MPRIP, PAQR9, PRSS1, and ZSCAN30.
Discussion
There is an ever-growing literature on the biological and pathological implications of STRs at the inter- and intraspecies levels [16,17,18,19,20,21,22,23,24,25,26,27]. The STRs listed in the present study are genetic codes that are unique to humans and are likely to be responsible for the human-specific regulation of the relevant genes. The significant enrichment of human-specific transcripts at the human-specific STR compartment indicates a link to a mechanism for TSS selection and transcript specificity.
A number of the identified STRs such as GTTT have established repressor activity [6, 28, 29] and are differentially expanded in certain genes in the Old World monkeys and Apes [14]. Purine STRs such as GAAA repeats are also functional in gene expression regulation, and their link to certain diseases unique to humans were previously reported [30, 31]. While the CG-rich STRs (e.g., CCG, GGGC) are subject to DNA methylation and can repress gene expression activity [32], they can also form G4 quadruplex structures, which have significant functions in gene expression regulation [33]. Several other identified STRs can form G4 structures with high overlap fraction (e.g., AGGG/CCCT, GCCC/GGGC).
It is not possible to estimate the number of crucial events that have led to the emergence of the human species. However, only a few genetic changes are needed to spur the evolution of new species in general, exemplified by the highly restricted initial divergence in butterfly hybridization models [34]. Accelerated evolution of a number of the identified genes in the present study (e.g., DRD2) has a well-established role in the origin of Homo sapiens [35]. Remarkably, a human-specific 7-amino acid transcript of this gene is flanked by a human-specific GCCC-repeat. Human-specific transcripts are increasingly recognized of having a role in the pathogenesis of diseases unique to the human species, such as schizophrenia [25, 36].
In a number of instances, not only the STR and the transcript, but also the gene containing these STRs and transcripts, were unique to humans, e.g., AMY1C, which is indicated in the evolution of the human phenotype during the Pleistocene [37].
For a number of the identified genes, sparse literature is available on the relevant function and pathways (e.g., ARHGEF35, CXorf40A, C22orf24, TMEM99, and ARL17B).
In a number of the identified genes, the STRs were linked to noncanonical (non-AUG) translation. Although the significance of this compartment is unknown for the most part, recent emerging data indicates likely biological functions [38].
The plasticity of STRs confers them unique ability to respond to adaptive evolutionary processes in a more efficient way than the quaternary codes provided by the SNSs. This potential aspect of STRs is vastly unknown at present, and it is expected that identification of STRs that have evolved differentially in humans vs. other species may pave the way for better understanding of the evolutionary implication of these highly mutable motifs.
This study warrants expansion to other vitally important gene regulatory sequences such as the distal promoter, 5′UTR, and 3′UTR. It is also necessary to sequence these STRs in characteristics and diseases that are unique to the human species. The recent reports of mass STR analysis using CRISPR/Cas9 [39] make it particularly more feasible to investigate STRs in the context of human evolution.
Conclusion
Our findings provide the first evidence of genome-scale skewing of STRs at a specific region of the human genome, and support a link between STRs and TSS selection/transcript specificity. The genes and STRs listed here may have a role in the divergence of humans from other species through the development of characteristics and phenotypes that are unique to the human species.
Abbreviations
- SNS:
-
Single nucleotide substitution
- STR:
-
Short tandem repeat
- TF:
-
Transcription factor
- TSS:
-
Transcription start site
- UTR:
-
Untranslated region
References
Press MO, Carlson KD, Queitsch C. The overdue promise of short tandem repeat variation for heritability. Trends Genet. 2014;30(11):504–12.
King DG. Evolution of simple sequence repeats as mutable sites. Adv Exp Med Biol. 2012;769:10-25.
Ohadi M, Mohammadparast S, Darvish H. Evolutionary trend of exceptionally long human core promoter short tandem repeats. Gene. 2012;507(1):61–7.
Mohammadparast S, Bayat H, Biglarian A, Ohadi M. Exceptional expansion and conservation of a CT-repeat complex in the core promoter of PAXBP1 in primates. Am J Primatol. 2014;76(8):747–56.
Nikkhah M, Rezazadeh M, Khorshid HK, Biglarian A, Ohadi M. An exceptionally long CA-repeat in the core promoter of SCGB2B2 links with the evolution of apes and Old World monkeys. Gene. 2016;576(1):109–14.
Rezazadeh M, Gharesouran J, Mirabzadeh A, Khorshid HK, Biglarian A, Ohadi M. A primate-specific functional GTTT-repeat in the core promoter of CYTH4 is linked to bipolar disorder in human. Prog Neuro-Psychopharmacol Biol Psychiatry. 2015;56:161–7.
Namdar-Aligoodarzi P, Mohammadparast S, Zaker-Kandjani B, Kakroodi ST, Vesiehsari MJ, Ohadi M. Exceptionally long 5′ UTR short tandem repeats specifically linked to primates. Gene. 2015;569(1):88–94.
Ohadi M, Valipour E, Ghadimi-Haddadan S, Namdar-Aligoodarzi P, Bagheri A, Kowsari A, et al. Core promoter short tandem repeats as evolutionary switch codes for primate speciation. Am J Primatol. 2015;77(1):34–43.
Sonay TB, Carvalho T, Robinson MD, Greminger MP, Krützen M, Comas D, et al. Tandem repeat variation in human and great ape populations and its impact on gene expression divergence. Genome Res. 2015;25(11):1591–9.
Li C, Lenhard B, Luscombe NM. Integrated analysis sheds light on evolutionary trajectories of young transcription start sites in the human genome. bioRxiv. 2017. https://doi.org/10.1101/192757.
Kramer M, Sponholz C, Slaba M, Wissuwa B, Claus RA, Menzel U, et al. Alternative 5′untranslated regions are involved in expression regulation of human heme oxygenase-1. PLoS One. 2013;8(10):e77224.
Li Y, Seidel K, Marschall P, Klein M, Hope A, Schacherl J, et al. A polymorphic microsatellite repeat within the ECE-1c promoter is involved in transcriptional start site determination, human evolution, and Alzheimer’s disease. J Neurosci. 2012;32(47):16807–20.
Emamalizadeh B, Movafagh A, Darvish H, Kazeminasab S, Andarva M, Namdar-Aligoodarzi P, et al. The human RIT2 core promoter short tandem repeat predominant allele is species-specific in length: a selective advantage for human evolution? Mol Gen Genomics. 2017;292(3):611–7.
Khademi E, Alehabib E, Shandiz EE, Ahmadifard A, Andarva M, Jamshidi J, et al. Support for “disease-only” genotypes and excess of homozygosity at the CYTH4 primate-specific GTTT-repeat in schizophrenia. Genet Test Mol Biomarkers. 2017;21(8):485–90.
Bushehri A, Barez MM, Mansouri S, Biglarian A, Ohadi M. Genome-wide identification of human-and primate-specific core promoter short tandem repeats. Gene. 2016;587(1):83–90.
Valipour E, Kowsari A, Bayat H, Banan M, Kazeminasab S, Mohammadparast S, et al. Polymorphic core promoter GA-repeats alter gene expression of the early embryonic developmental genes. Gene. 2013;531(2):175–9.
Heidari A, Fam ZNS, Esmaeilzadeh-Gharehdaghi E, Banan M, Hosseinkhani S, Mohammadparast S, et al. Core promoter STRs: novel mechanism for inter-individual variation in gene expression in humans. Gene. 2012;492(1):195–8.
Bagshaw AT. Functional mechanisms of microsatellite DNA in eukaryotic genomes. Genome Biol Evol. 2017;9(9):2428–43.
Bagshaw AT, Horwood LJ, Fergusson DM, Gemmell NJ, Kennedy MA. Microsatellite polymorphisms associated with human behavioural and psychological phenotypes including a gene-environment interaction. BMC Med Genet. 2017;18(1):12.
Gymrek M, Willems T, Guilmatre A, Zeng H, Markus B, Georgiev S, et al. Abundant contribution of short tandem repeats to gene expression variation in humans. Nat Genet. 2016;48(1):22–9.
Hammock EA, Young LJ. Microsatellite instability generates diversity in brain and sociobehavioral traits. Science. 2005;308(5728):1630–4.
Carrat GR, Hu M, Nguyen-Tu M-S, Chabosseau P, Gaulton KJ, van de Bunt M, et al. Decreased STARD10 expression is associated with defective insulin secretion in humans and mice. Am J Hum Genet. 2017;100(2):238–56.
Abe H, Gemmell NJ. Evolutionary footprints of short tandem repeats in avian promoters. Sci Rep. 2016;6:19421.
Fondon JW, Hammock EA, Hannan AJ, King DG. Simple sequence repeats: genetic modulators of brain function and behavior. Trends Neurosci. 2008;31(7):328–34.
Alizadeh F, Bozorgmehr A, Tavakkoly-Bazzaz J, Ohadi M. Skewing of the genetic architecture at the ZMYM3 human-specific 5′ UTR short tandem repeat in schizophrenia. Mol Gen Genomics. 2018; https://doi.org/10.1007/s00438-018-1415-8.
Hannan AJ. Tandem repeats mediating genetic plasticity in health and disease. Nat Rev Genet. 2018; https://doi.org/10.1038/nrg.2017.115.
Watts PC, Kallio ER, Koskela E, Lonn E, Mappes T, Mokkonen M. Stabilizing selection on microsatellite allele length at arginine vasopressin 1a receptor and oxytocin receptor loci. Proc Biol Sci. 2017;284(1869):20171896. https://doi.org/10.1098/rspb.2017.1896.
Andrioli LPM, Vasisht V, Theodosopoulou E, Oberstein A, Small S. Anterior repression of a Drosophila stripe enhancer requires three position-specific mechanisms. Development. 2002;129(21):4931–40.
Andrioli LP, Oberstein AL, Corado MS, Yu D, Small S. Groucho-dependent repression by sloppy-paired 1 differentially positions anterior pair-rule stripes in the Drosophila embryo. Dev Biol. 2004;276(2):541–51.
Darvish H, Heidari A, Hosseinkhani S, Movafagh A, Khaligh A, Jamshidi J, et al. Biased homozygous haplotypes across the human caveolin 1 upstream purine complex in Parkinson’s disease. J Mol Neurosci. 2013;51(2):389–93.
Heidari A, Hosseinkhani S, Talebi S, Meshkani R, Esmaeilzadeh-Gharedaghi E, Banan M, et al. Haplotypes across the human caveolin 1 gene upstream purine complex significantly alter gene expression: implication in neurodegenerative disorders. Gene. 2012;505(1):186–9.
Quilez J, Guilmatre A, Garg P, Highnam G, Gymrek M, Erlich Y, et al. Polymorphic tandem repeats within gene promoters act as modifiers of gene expression and DNA methylation in humans. Nucleic Acids Res. 2016;44(8):3750–62.
Sawaya S, Bagshaw A, Buschiazzo E, Kumar P, Chowdhury S, Black MA, et al. Microsatellite tandem repeats are abundant in human promoters and are associated with regulatory elements. PLoS One. 2013;8(2):e54710.
Kronforst MR, Hansen ME, Crawford NG, Gallant JR, Zhang W, Kulathinal RJ, et al. Hybridization reveals the evolving genomic architecture of speciation. Cell Rep. 2013;5(3):666–77.
Dorus S, Vallender EJ, Evans PD, Anderson JR, Gilbert SL, Mahowald M, et al. Accelerated evolution of nervous system genes in the origin of Homo sapiens. Cell. 2004;119(7):1027–40.
Li M, Jaffe AE, Straub RE, Tao R, Shin JH, Wang Y, et al. A human-specific AS3MT isoform and BORCS7 are molecular risk factors in the 10q24. 32 schizophrenia-associated locus. Nat Med. 2016;22(6):649–56.
Hardy K, Brand-Miller J, Brown KD, Thomas MG, Copeland L. The importance of dietary carbohydrate in human evolution. Q Rev Biol. 2015;90(3):251–68.
Na CH, Barbhuiya MA, Kim MS, Verbruggen S, Eacker SM, Pletnikova O, Troncoso JC, Halushka MK, Menschaert G, Overall CM, Pandey A. Discovery of noncanonical translation initiation sites through mass spectrometric analysis of protein N termini. Genome Res. 2018;28(1):25–36.
Shin G, Grimes SM, Lee H, Lau BT, Xia LC, Ji HP. CRISPR–Cas9-targeted fragmentation and selective sequencing enable massively parallel microsatellite analysis. Nat Commun. 2017;8:14291.
Funding
This research was funded by the University of Social Welfare and Rehabilitation Sciences, Tehran, Iran.
Availability of data and materials
Please contact author for data requests.
Author information
Authors and Affiliations
Contributions
NN, FA, and TAA carried out the bioinformatics studies. AD and RS were the advisors of the project and helped in coordination. MO conceived the study, designed the project, supervised the analysis and wrote the manuscript. All authors read and approved the final manuscript.
Corresponding author
Ethics declarations
Ethics approval and consent to participate
Not applicable.
Consent for publication
Not applicable.
Competing interests
The authors declare that they have no competing interests.
Publisher’s Note
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Rights and permissions
Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.
About this article
Cite this article
Nazaripanah, N., Adelirad, F., Delbari, A. et al. Genome-scale portrait and evolutionary significance of human-specific core promoter tri- and tetranucleotide short tandem repeats. Hum Genomics 12, 17 (2018). https://doi.org/10.1186/s40246-018-0149-3
Received:
Accepted:
Published:
DOI: https://doi.org/10.1186/s40246-018-0149-3
Keywords
- Short tandem repeat
- Core promoter
- Human-specific
- Trinucleotide
- Tetranucleotide