Genome-scale portrait and evolutionary significance of human-specific core promoter tri- and tetranucleotide short tandem repeats

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 huntergatherer 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 pentaand 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 triand 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 humanspecific 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.    The numbers before the brackets represent the start site of the STR in respect of the corresponding transcription start site. "Variant no" corresponds to the Ensembl isoform number  The numbers before the brackets represent the start site of the STR in respect of the corresponding transcription start site. "Variant no" corresponds to the Ensembl isoform number 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 interand 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   [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., ARH-GEF35, 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 genomescale 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. Availability of data and materials Please contact author for data requests.