Genome-wide identification of human- and primate-specific core promoter short tandem repeats

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.

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.

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.

Transcriptional Interference Regulates the Evolutionary Development of Speech

The human capacity to speak is fundamental to our advanced intellectual, technological and social development. Yet so very little is known regarding the evolutionary genetics of speech or its relationship with the broader aspects of evolutionary development in primates. In this study, we describe a large family with evolutionary retrograde development of the larynx and wrist. The family presented with severe speech impairment and incremental retrograde elongations of the pisiform in the wrist that limited wrist rotation from 180° to 90° as in primitive primates. To our surprise, we found that a previously unknown primate-specific gene TOSPEAK had been disrupted in the family. TOSPEAK emerged de novo in an ancestor of extant primates across a 540 kb region of the genome with a pre-existing highly conserved long-range laryngeal enhancer for a neighbouring bone morphogenetic protein gene GDF6. We used transgenic mouse modelling to identify two additional GDF6 long-range enhancers within TOSPEAK that regulate GDF6 expression in the wrist. Disruption of TOSPEAK in the affected family blocked the transcription of TOSPEAK across the 3 GDF6 enhancers in association with a reduction in GDF6 expression and retrograde development of the larynx and wrist. Furthermore, we describe how TOSPEAK developed a human-specific promoter through the expansion of a penta-nucleotide direct repeat that first emerged de novo in the promoter of TOSPEAK in gibbon. This repeat subsequently expanded incrementally in higher hominids to form an overlapping series of Sp1/KLF transcription factor consensus binding sites in human that correlated with incremental increases in the promoter strength of TOSPEAK with human having the strongest promoter. Our research indicates a dual evolutionary role for the incremental increases in TOSPEAK transcriptional interference of GDF6 enhancers in the incremental evolutionary development of the wrist and larynx in hominids and the human capacity to speak and their retrogression with the reduction of TOSPEAK transcription in the affected family.

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.

Suppression of lncRNA NLRP3 inhibits NLRP3-triggered inflammatory responses in early acute lung injury

Acute lung injury (ALI) is a common lung pathology that is accompanied by alveolar macrophage (AM) activation and inflammatory response. This study investigated the role of the long non-coding RNA NONRATT004344 (hereafter named lncRNA NLRP3) in regulating the Nod-like receptor protein 3 (NLRP3)-triggered inflammatory response in early ALI and the underlying mechanism as well. We established LPS-induced ALI models to explore their interactive mechanisms in vitro and in vivo. Luciferase reporter assays were performed to determine that miR-138-5p could bind to lncRNA NLRP3 and NLRP3. We observed increased lncRNA NLRP3 expression, decreased miR-138-5p expression, NLRP3 inflammasome activation, and upregulated caspase-1, IL-1β, and IL-18 expression in the LPS-induced ALI model. Furthermore, lncRNA NLRP3 overexpression activated the NLRP3 inflammasome and promoted IL-1β and IL-18 secretion; the miR-138-5p mimic abolished these effects in vivo and in vitro. Consistently, miR-138-5p inhibition reversed the effects of lncRNA NLRP3 silencing on the expression of NLRP3-related molecules and inhibition of the NLRP3/caspase-1/IL-1β signalling pathway. Mechanistically, lncRNA NLRP3 sponging miR-138-5p facilitated NLRP3 activation through a competitive endogenous RNA (ceRNA) mechanism. In summary, our results suggested that lncRNA NLRP3 binding miR-138-5p promotes NLRP3-triggered inflammatory response via lncRNA NLRP3/miR-138-5p/NLRP3 ceRNA network (ceRNET) and provides insights into the treatment of early ALI.

Natural selection at the RASGEF1C (GGC) repeat in human and divergent genotypes in late-onset neurocognitive disorder

Expression dysregulation of the neuron-specific gene, RASGEF1C (RasGEF Domain Family Member 1C), occurs in late-onset neurocognitive disorders (NCDs), such as Alzheimer's disease. This gene contains a (GGC)13, spanning its core promoter and 5' untranslated region (RASGEF1C-201 ENST00000361132.9). Here we sequenced the (GGC)-repeat in a sample of human subjects (N = 269), consisting of late-onset NCDs (N = 115) and controls (N = 154). We also studied the status of this STR across various primate and non-primate species based on Ensembl 103. The 6-repeat allele was the predominant allele in the controls (frequency = 0.85) and NCD patients (frequency = 0.78). The NCD genotype compartment consisted of an excess of genotypes that lacked the 6-repeat (divergent genotypes) (Mid-P exact = 0.004). A number of those genotypes were not detected in the control group (Mid-P exact = 0.007). The RASGEF1C (GGC)-repeat expanded beyond 2-repeats specifically in primates, and was at maximum length in human. We conclude that there is natural selection for the 6-repeat allele of the RASGEF1C (GGC)-repeat in human, and significant divergence from that allele in late-onset NCDs. STR alleles that are predominantly abundant and genotypes that deviate from those alleles are underappreciated features, which may have deep evolutionary and pathological consequences.

Natural Selection at the NHLH2 Core Promoter Exceptionally Long CA-Repeat in Human and Disease-Only Genotypes in Late-Onset Neurocognitive Disorder

BACKGROUND

Approximately 2% of the human core promoter short tandem repeats (STRs) reach lengths of ≥6 repeats, which may in part be a result of adaptive evolutionary processes and natural selection. A single-exon transcript of the human nescient helix loop helix 2 (NHLH2) gene is flanked by the longest CA-repeat detected in a human protein-coding gene core promoter (Ensembl transcript ID: ENST00000369506.1). NHLH2 is involved in several biological and pathological pathways, such as motivated exercise, obesity, and diabetes.

METHODS

The allele and genotype distribution of the NHLH2 CA-repeat were investigated by sequencing in 655 Iranian subjects, consisting of late-onset neurocognitive disorder (NCD) as a clinical entity (n = 290) and matched controls (n = 365). The evolutionary trend of the CA-repeat was also studied across vertebrates.

RESULTS

The allele range was between 9 and 25 repeats in the NCD cases, and 12 and 24 repeats in the controls. At the frequency of 0.56, the 21-repeat allele was the predominant allele in the controls. While the 21-repeat was also the predominant allele in the NCD patients, we detected significant decline of the frequency (p < 0.0001) and homozygosity (p < 0.006) of this allele in this group. Furthermore, 12 genotypes were detected across 16 patients (5.5% of the entire NCD sample) and not in the controls (disease-only genotypes; p < 0.0003), consisting of at least one extreme allele. The extreme alleles were at 9, 12, 13, 18, and 19 repeats (extreme short end), and 23, 24, and 25 repeats (extreme long end), and their frequencies ranged between 0.001 and 0.04. The frequency of the 21-repeat allele significantly dropped to 0.09 in the disease-only genotype compartment (p < 0.0001). Evolutionarily, while the maximum length of the NHLH2 CA-repeat was 11 repeats in non-primates, this CA-repeat was ≥14 repeats in primates and reached maximum length in human.

CONCLUSION

We propose a novel locus for late-onset NCD at the NHLH2 core promoter exceptionally long CA-STR and natural selection at this locus. Furthermore, there was indication of genotypes at this locus that unambiguously linked to late-onset NCD. This is the first instance of natural selection in favor of a predominantly abundant STR allele in human and its differential distribution in late-onset NCD.

Link between short tandem repeats and translation initiation site selection

BACKGROUND

Despite their vast biological implication, the relevance of short tandem repeats (STRs)/microsatellites to the protein-coding gene translation initiation sites (TISs) remains largely unknown.

METHODS

We performed an Ensembl-based comparative genomics study of all annotated orthologous TIS-flanking sequences in human and 46 other species across vertebrates, on the genomic DNA and cDNA platforms (755,956 TISs), aimed at identifying human-specific STRs in this interval. The collected data were used to examine the hypothesis of a link between STRs and TISs. BLAST was used to compare the initial five amino acids (excluding the initial methionine), codons of which were flanked by STRs in human, with the initial five amino acids of all annotated proteins for the orthologous genes in other vertebrates (total of 5,314,979 pair-wise TIS comparisons on the genomic DNA and cDNA platforms) in order to compare the number of events in which human-specific and non-specific STRs occurred with homologous and non-homologous TISs (i.e., ≥ 50% and < 50% similarity of the five amino acids).

RESULTS

We detected differential distribution of the human-specific STRs in comparison to the overall distribution of STRs on the genomic DNA and cDNA platforms (Mann Whitney U test p = 1.4 × 10-11 and p < 7.9 × 10-11, respectively). We also found excess occurrence of non-homologous TISs with human-specific STRs and excess occurrence of homologous TISs with non-specific STRs on both platforms (p < 0.00001).

CONCLUSION

We propose a link between STRs and TIS selection, based on the differential co-occurrence rate of human-specific STRs with non-homologous TISs and non-specific STRs with homologous TISs.

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.

Support for "Disease-Only" Genotypes and Excess of Homozygosity at the CYTH4 Primate-Specific GTTT-Repeat in Schizophrenia

OBJECTIVE

The role of short tandem repeats (STRs) in the control of gene expression among species is being increasingly understood following the identification of several instances in which certain STRs occur identically, or expand differentially, in primates versus nonprimates. These STRs may regulate genes that participate in characteristics that are associated with the divergence of primates from sibling orders (e.g., brain higher order functions). The CYTH4 gene contains the longest tetranucleotide STR in its core promoter, at 7-repeats, and links to the evolution of human and nonhuman primates. Allele and genotype distribution of this STR were studied in patients affected by schizophrenia (SCZ) and controls.

METHODS

High-resolution data were obtained on the allele and genotype distribution of the CYTH4 STR and a novel C > T single-nucleotide polymorphism (SNP) at its immediate upstream sequence in 255 patients with SCZ and 249 controls. Each sample was sequenced twice using the fluorescent dye termination method.

RESULTS

Novel alleles were detected at the long extreme of the GTTT-repeat, at 10- and 11-repeats, in the SCZ cases and controls. Excess of homozygosity was observed for the entire range of alleles across the GTTT-repeat and the C > T SNP in the SCZ patients in comparison with the controls (Yates corrected p < 0.011). Three genotypes consisting of the 11-repeat allele (i.e., 11/11, 10/11, and 7/11) were detected only in the SCZ patients (i.e., disease-only genotypes), and contributed to 2.3% of the SCZ genotypes (Mid p exact <0.007). The frequency of the 11-repeat allele was estimated at 0.02 and 0.006 in the SCZ patients and controls, respectively (Mid p exact <0.006).

CONCLUSION

This indicates that STR genotypes that are absent in the control group may be risk factors for SCZ. Future studies are warranted to test the significance of our findings.

LincRNA-Gm4419 knockdown ameliorates NF-κB/NLRP3 inflammasome-mediated inflammation in diabetic nephropathy

Diabetic nephropathy (DN) as the primary cause of end-stage kidney disease is a common complication of diabetes. Recent researches have shown the activation of nuclear factor kappa light-chain enhancer of activated B cells (NF-κB) and NACHT, LRR and PYD domain-containing protein 3 (NLRP3) inflammasome are associated with inflammation in the progression of DN, but the exact mechanism is unclear. Long noncoding RNAs (lncRNAs) have roles in the development of many diseases including DN. However, the relationship between lncRNAs and inflammation in DN remains largely unknown. Our previous study has revealed that 14 lncRNAs are abnormally expressed in DN by RNA sequencing and real-time quantitative PCR (qRT-PCR) in the renal tissues of db/db DN mice. In this study, these lncRNAs were verified their expressions by qRT-PCR in mesangial cells (MCs) cultured under high- and low-glucose conditions. Twelve lncRNAs displayed the same expressional tendencies in both renal tissues and MCs. In particular, long intergenic noncoding RNA (lincRNA)-Gm4419 was the only one associating with NF-κB among these 12 lncRNAs by bioinformatics methods. Moreover, Gm4419 knockdown could obviously inhibit the expressions of pro-inflammatory cytokines and renal fibrosis biomarkers, and reduce cell proliferation in MCs under high-glucose condition, whereas overexpression of Gm4419 could increase the inflammation, fibrosis and cell proliferation in MCs under low-glucose condition. Interestingly, our results showed that Gm4419 could activate the NF-κB pathway by directly interacting with p50, the subunit of NF-κB. In addition, we found that p50 could interact with NLRP3 inflammasome in MCs. In conclusion, our findings suggest lincRNA-Gm4419 may participate in the inflammation, fibrosis and proliferation in MCs under high-glucose condition through NF-κB/NLRP3 inflammasome signaling pathway, and may provide new insights into the regulation of Gm4419 during the progression of DN.