Evolutionary origins and interactomes of human, young microproteins and small peptides translated from short open reading frames.

All species continuously evolve short open reading frames (sORFs) that can be templated for protein synthesis and may provide raw materials for evolutionary adaptation. We analyzed the evolutionary origins of 7,264 recently cataloged human sORFs and found that most were evolutionarily young and had emerged de novo. We additionally identified 221 previously missed sORFs potentially translated into peptides of up to 15 amino acids-all of which are smaller than the smallest human microprotein annotated to date. To investigate the bioactivity of sORF-encoded small peptides and young microproteins, we subjected 266 candidates to a mass-spectrometry-based interactome screen with motif resolution. Based on these interactomes and additional cellular assays, we can associate several candidates with mRNA splicing, translational regulation, and endocytosis. Our work provides insights into the evolutionary origins and interaction potential of young and small proteins, thereby helping to elucidate this underexplored territory of the human proteome.

Splicing and editing of ionotropic glutamate receptors: a comprehensive analysis based on human RNA-Seq data.

Ionotropic glutamate receptors (iGluRs) play key roles for signaling in the central nervous system. Alternative splicing and RNA editing are well-known mechanisms to increase iGluR diversity and to provide context-dependent regulation. Earlier work on isoform identification has focused on the analysis of cloned transcripts, mostly from rodents. We here set out to obtain a systematic overview of iGluR splicing and editing in human brain based on RNA-Seq data. Using data from two large-scale transcriptome studies, we established a workflow for the de novo identification and quantification of alternative splice and editing events. We detected all canonical iGluR splice junctions, assessed the abundance of alternative events described in the literature, and identified new splice events in AMPA, kainate, delta, and NMDA receptor subunits. Notable events include an abundant transcript encoding the GluA4 amino-terminal domain, GluA4-ATD, a novel C-terminal GluD1 (delta receptor 1) isoform, GluD1-b, and potentially new GluK4 and GluN2C isoforms. C-terminal GluN1 splicing may be controlled by inclusion of a cassette exon, which shows preference for one of the two acceptor sites in the last exon. Moreover, we identified alternative untranslated regions (UTRs) and species-specific differences in splicing. In contrast, editing in exonic iGluR regions appears to be mostly limited to ten previously described sites, two of which result in silent amino acid changes. Coupling of proximal editing/editing and editing/splice events occurs to variable degree. Overall, this analysis provides the first inventory of alternative splicing and editing in human brain iGluRs and provides the impetus for further transcriptome-based and functional investigations.

Evidence for miR-548c-5p regulation of FOXC2 transcription through a distal genomic target site in human podocytes.

Podocytes are highly differentiated epithelial cells outlining the glomerular vessels. FOXC2 is a transcription factor essential for inducing podocyte differentiation, development and maturation, and is considered to be the earliest podocyte marker. miRNA prediction analysis revealed a full-length target site for the primate-specific miR-548c-5p at a genomic region > 8 kb upstream of FOXC2. We hypothesised that the transcription rates of FOXC2 during podocyte differentiation might be tuned by miR-548c-5p through this target site. Experiments were performed with cultured human podocytes, transfected with luciferase reporter constructs bearing this target site region within an enhancer element of the native plasmid. The results confirmed a seed region-driven targeting potential by the miRNA, with mimics downregulating and inhibitors enhancing luciferase activity. Introducing mutations into the miRNA target seed region abolished the expected response. In cultured podocytes, FOXC2 mRNA and protein levels responded to miR-548c-5p abundance in a coordinated manner before and after induction of differentiation, with high statistical significance. Ago-ChIP experiments revealed occupancy of the miRNA target site by miRNA/RISC in undifferentiated cells and its release when differentiation is initiated, allowing its interaction with the gene's promoter region to amplify FOXC2 expression, as shown by chromosome conformation capture and qRT-PCR. Moreover, the expression pattern of FOXC2 during podocyte differentiation seems to be affected by miR-548c-5p, as removal of either endogenous or mimic miR-548c-5p results in increased FOXC2 protein levels and cells resembling those undergoing differentiation. Collectively, results indicate a well-orchestrated regulatory model of FOXC2 expression by a remote upstream target site for miR-548c-5p.

Two distinct nuclear stress bodies containing different sets of RNA-binding proteins are formed with HSATIII architectural noncoding RNAs upon thermal stress exposure.

Nuclear stress bodies (nSBs) are thermal stress-inducible membrane-less nuclear bodies that are formed on highly repetitive satellite III architectural noncoding RNAs (HSATIII arcRNAs). Upon thermal stress exposure, HSATIII expression is induced to sequestrate specific sets of RNA-binding proteins and form nSBs. The major population of nSBs contain SAFB as a marker, whereas the minor population are SAFB-negative. Here, we found that HNRNPM, which was previously reported to localize in nuclear foci adjacent to SAFB-positive foci upon thermal stress, localizes in a minor population of HSATIII-dependent nSBs. Hence, we used the terms nSB-S and nSB-M to distinguish the SAFB foci and HNRNPM foci, respectively. Analysis of the components of the nSBs revealed that each set contains distinct RNA-binding proteins, including SLTM and NCO5A in nSB-Ss and HNRNPA1 and HNRNPH1 in nSB-Ms. Overall, our findings indicate that two sets of nSBs containing HSATIII arcRNAs and distinct sets of RNA-binding proteins are formed upon thermal stress exposure.

Massive contribution of transposable elements to mammalian regulatory sequences.

Barbara McClintock discovered the existence of transposable elements (TEs) in the late 1940s and initially proposed that they contributed to the gene regulatory program of higher organisms. This controversial idea gained acceptance only much later in the 1990s, when the first examples of TE-derived promoter sequences were uncovered. It is now known that half of the human genome is recognizably derived from TEs. It is thus important to understand the scope and nature of their contribution to gene regulation. Here, we provide a timeline of major discoveries in this area and discuss how transposons have revolutionized our understanding of mammalian genomes, with a special emphasis on the massive contribution of TEs to primate evolution. Our analysis of primate-specific functional elements supports a simple model for the rate at which new functional elements arise in unique and TE-derived DNA. Finally, we discuss some of the challenges and unresolved questions in the field, which need to be addressed in order to fully characterize the impact of TEs on gene regulation, evolution and disease processes.

PRINS, a primate-specific long non-coding RNA, plays a role in the keratinocyte stress response and psoriasis pathogenesis.

In the last few years with the recent emergence of high-throughput technologies, thousands of long non-coding RNAs (lncRNAs) have been identified in the human genome. However, assigning functional annotation and determining cellular contexts for these RNAs are still in its infancy. As information gained about lncRNA structure, interacting partners, and roles in human diseases may be helpful in the characterization of novel lncRNAs, we review our knowledge on a selected group of lncRNAs that were identified serendipitously years ago by large-scale gene expression methods used to study human diseases. In particular, we focus on the Psoriasis-susceptibility-Related RNA Gene Induced by Stress (PRINS) lncRNA, first identified by our research group as a transcript highest expressed in psoriatic non-lesional epidermis. Results gathered for PRINS in the last 10 years indicate that it is conserved in primates and plays a role in keratinocyte stress response. Elevated levels of PRINS expression in psoriatic non-lesional keratinocytes alter the stress response of non-lesional epidermis and contribute to disease pathogenesis. Finally, we propose a categorization for the PRINS lncRNA based on a recently elaborated system for lncRNA classification.

MicroRNAs in placental health and disease.

MicroRNAs (miRNAs) constitute a large family of small noncoding RNAs that are encoded by the genomes of most organisms. They regulate gene expression through posttranscriptional mechanisms to attenuate protein output in various genetic networks. The discovery of miRNAs has transformed our understanding of gene regulation and sparked intense efforts intended to harness their potential as diagnostic markers and therapeutic tools. Over the last decade, a flurry of studies has shed light on placental miRNAs but has also raised many questions regarding the scope of their biologic action. Moreover, the recognition that miRNAs of placental origin are released continually in the maternal circulation throughout pregnancy suggested that circulating miRNAs might serve as biomarkers for placental function during pregnancy. Although this generated much enthusiasm, recently recognized challenges have delayed the application of miRNA-based biomarkers and therapeutics in clinical practice. In this review, we summarize key findings in the field and discuss current knowledge related to miRNAs in the context of placental biology.

Leveraging primate-specific genomic information for genetic studies of complex diseases.

Genomic changes specific to higher primates are regarded as primate-specific genomic information (PSI). Using PSI to inform genetic studies is highly desirable but hampered by three factors: heterogeneity among PSI studies, lack of integrated profiles of the identified PSI elements and dearth of relevant functional information. We report a database of 19,767 PSI elements collated from nine types of brain-related studies, which form 19,473 non-overlapping PSI regions that distribute unevenly but jointly cover only 0.81% of the genome. About 2.5% of the PSI regions colocalized with variants identified in genome-wide association studies, with disease loci more likely colocalized than quantitative trait loci (p = 1.6 × 10-5), particularly in regions without obvious regulatory roles. We further showed an LRP4 exemplar region with PSI elements orchestrated with common and rare disease variants and other functional elements. Our results render PSI elements as a valuable source to inform genetic studies of complex diseases.

Pan-cancer surveys indicate cell cycle-related roles of primate-specific genes in tumors and embryonic cerebrum.

BACKGROUND: Despite having been extensively studied, it remains largely unclear why humans bear a particularly high risk of cancer. The antagonistic pleiotropy hypothesis predicts that primate-specific genes (PSGs) tend to promote tumorigenesis, while the molecular atavism hypothesis predicts that PSGs involved in tumors may represent recently derived duplicates of unicellular genes. However, these predictions have not been tested. RESULTS: By taking advantage of pan-cancer genomic data, we find the upregulation of PSGs across 13 cancer types, which is facilitated by copy-number gain and promoter hypomethylation. Meta-analyses indicate that upregulated PSGs (uPSGs) tend to promote tumorigenesis and to play cell cycle-related roles. The cell cycle-related uPSGs predominantly represent derived duplicates of unicellular genes. We prioritize 15 uPSGs and perform an in-depth analysis of one unicellular gene-derived duplicate involved in the cell cycle, DDX11. Genome-wide screening data and knockdown experiments demonstrate that DDX11 is broadly essential across cancer cell lines. Importantly, non-neutral amino acid substitution patterns and increased expression indicate that DDX11 has been under positive selection. Finally, we find that cell cycle-related uPSGs are also preferentially upregulated in the highly proliferative embryonic cerebrum. CONCLUSIONS: Consistent with the predictions of the atavism and antagonistic pleiotropy hypotheses, primate-specific genes, especially those PSGs derived from cell cycle-related genes that emerged in unicellular ancestors, contribute to the early proliferation of the human cerebrum at the cost of hitchhiking by similarly highly proliferative cancer cells.

Intragenic Transcriptional cis-Antagonism Across SLC6A3.

A promoter can be regulated by various cis-acting elements so that delineation of the regulatory modes among them may help understand developmental, environmental and genetic mechanisms in gene activity. Here we report that the human dopamine transporter gene SLC6A3 carries a 5' distal 5-kb super enhancer (5KSE) which upregulated the promoter by 5-fold. Interestingly, 5KSE is able to prevent 3' downstream variable number tandem repeats (3'VNTRs) from silencing the promoter. This new enhancer consists of a 5'VNTR and three repetitive sub-elements that are conserved in primates. Two of 5KSE's sub-elements, E-9.7 and E-8.7, upregulate the promoter, but only the later could continue doing so in the presence of 3'VNTRs. Finally, E-8.7 is activated by novel dopaminergic transcription factors including SRP54 and Nfe2l1. Together, these results reveal a multimodal regulatory mechanism in SLC6A3.

Primate-specific miR-944 activates p53-dependent tumor suppression in human colorectal cancers.

As cancers with a high incidence rate, colorectal cancers are a main cause of cancer-related death. MicroRNAs are often deregulated in cancers. The primate-specific miR-944, located in a p63 intron, is known to be highly expressed in patients exhibiting low colorectal cancer recurrence rates. However, the biological functions of miR-944 in colorectal cancers remain unclear. In this study, we found that miR-944 was downregulated in colorectal cancer tissues, and inhibited cancer cell growth in a xenograft mouse model. The overexpression of miR-944 caused G1 phase arrest and increased p53 expression in cancer cells. p53 stability was enhanced by miR-944s targeting E3 ligases COP1 and MDM2. Overexpression of COP1 and MDM2 restored cell growth inhibition caused by miR-944. Taken together, our results suggest that miR-944 acts as a potential tumor suppressor in colorectal cancers through the ubiquitin-proteasome system.

Evolutionary Gain of Dbx1 Expression Drives Subplate Identity in the Cerebral Cortex.

Changes in transcriptional regulation through cis-regulatory elements are thought to drive brain evolution. However, how this impacts the identity of primate cortical neurons is still unresolved. Here, we show that primate-specific cis-regulatory sequences upstream of the Dbx1 gene promote human-like expression in the mouse embryonic cerebral cortex, and this imparts cell identity. Indeed, while Dbx1 is expressed in highly restricted cortical progenitors in the mouse ventral pallium, it is maintained in neurons in primates. Phenocopy of the primate-like Dbx1 expression in mouse cortical progenitors induces ectopic Cajal-Retzius and subplate (SP) neurons, which are transient populations playing crucial roles in cortical development. A conditional expression solely in neurons uncouples mitotic and postmitotic activities of Dbx1 and exclusively promotes a SP-like fate. Our results highlight how transcriptional changes of a single fate determinant in postmitotic cells may contribute to the expansion of neuronal diversity during cortical evolution.

Evolution and cell-type specificity of human-specific genes preferentially expressed in progenitors of fetal neocortex.

Understanding the molecular basis that underlies the expansion of the neocortex during primate, and notably human, evolution requires the identification of genes that are particularly active in the neural stem and progenitor cells of the developing neocortex. Here, we have used existing transcriptome datasets to carry out a comprehensive screen for protein-coding genes preferentially expressed in progenitors of fetal human neocortex. We show that 15 human-specific genes exhibit such expression, and many of them evolved distinct neural progenitor cell-type expression profiles and levels compared to their ancestral paralogs. Functional studies on one such gene, NOTCH2NL, demonstrate its ability to promote basal progenitor proliferation in mice. An additional 35 human genes with progenitor-enriched expression are shown to have orthologs only in primates. Our study provides a resource of genes that are promising candidates to exert specific, and novel, roles in neocortical development during primate, and notably human, evolution.

Primate-specific Long Non-coding RNAs and MicroRNAs.

Non-coding RNAs (ncRNAs) are critical regulators of gene expression in essentially all life forms. Long ncRNAs (lncRNAs) and microRNAs (miRNAs) are two important RNA classes possessing regulatory functions. Up to date, many primate-specific ncRNAs have been identified and investigated. Their expression specificity to primate lineage suggests primate-specific roles. It is thus critical to elucidate the biological significance of primate or even human-specific ncRNAs, and to develop potential ncRNA-based therapeutics. Here, we have summarized the studies regarding regulatory roles of some key primate-specific lncRNAs and miRNAs.

De novo methylation in male germ cells of the common marmoset monkey occurs during postnatal development and is maintained in vitro.

The timing of de novo DNA methylation in male germ cells during human testicular development is yet unsolved. Apart from that, the stability of established imprinting patterns in vitro is controversially discussed. This study aimed at determining the timing of DNA de novo methylation and at assessing the stability of the methylation status in vitro. We employed the marmoset monkey (Callithrix jacchus) as it is considered the best non-human primate model for human testicular development. We selected neonatal, pre-pubertal, pubertal, and adult animals (n = 3, each) and assessed germ cell global DNA methylation levels by 5-methyl cytosine staining, and Alu elements and gene-specific methylation (H19, LIT1, SNRPN, MEST, OCT4, MAGE-A4, and DDX-4) by pyrosequencing. De novo methylation is progressively established during postnatal primate development and continues until adulthood, a process that is different in most other species. Importantly, once established, methylation patterns remained stable, as demonstrated using in vitro cultures. Thus, the marmoset monkey is a unique model for the study of postnatal DNA methylation mechanisms in germ cells and for the identification of epimutations and their causes.

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.

Dysregulation of an X-linked primate-specific epididymal microRNA cluster in unexplained asthenozoospermia.

Asthenoszoopermia, characterized by reduced sperm motility, is one of the primary forms of male infertility. Whereas most cases were diagnosed into unexplained asthenozoospermia (UA) because the etiology cannot be identified. In animal models, epigenetic dysregulation in epididymis can impair sperm maturation and result in asthenozoospermia. However, researches of epididymal epigenetic regulation on humans are impeded by the difficulty in obtaining epididymal tissues. We previously identified cell-free seminal microRNAs predominately derived from epididymis in human ejaculate. In the present study, these microRNAs were used to screen and validate the microRNA dysregulation in men with UA, which were divided into screening set and validation set. The expression of five miRNAs (miR-891b, miR-892b, miR-892a, miR-888 and miR-890) was found and confirmed to be dysregulated in men with UA. Interestingly, these five miRNAs belong to a primate-specific miRNA cluster located on the X chromosome with epididymis specific expression. Moreover, obvious coherent dysregulation of these miRNAs were observed in 13% men with UA. Regression analysis demonstrated that levels of these miRNAs were significantly correlated with progressive sperm motility. Functions and pathways of predicted target genes of this cluster suggested its role in sperm maturation. Dysregulation of this miRNA cluster might be an epigenetic basis for some patients with UA. We also showed a noninvasive and feasible approach to get epigenetic information of human epididymis.

miRNA Enriched in Human Neuroblast Nuclei Bind the MAZ Transcription Factor and Their Precursors Contain the MAZ Consensus Motif.

While the cytoplasmic function of microRNA (miRNA) as post-transcriptional regulators of mRNA has been the subject of significant research effort, their activity in the nucleus is less well characterized. Here we use a human neuronal cell model to show that some mature miRNA are preferentially enriched in the nucleus. These molecules were predominantly primate-specific and contained a sequence motif with homology to the consensus MAZ transcription factor binding element. Precursor miRNA containing this motif were shown to have affinity for MAZ protein in nuclear extract. We then used Ago1/2 RIP-Seq to explore nuclear miRNA-associated mRNA targets. Interestingly, the genes for Ago2-associated transcripts were also significantly enriched with MAZ binding sites and neural function, whereas Ago1-transcripts were associated with general metabolic processes and localized with SC35 spliceosomes. These findings suggest the MAZ transcription factor is associated with miRNA in the nucleus and may influence the regulation of neuronal development through Ago2-associated miRNA induced silencing complexes. The MAZ transcription factor may therefore be important for organizing higher order integration of transcriptional and post-transcriptional processes in primate neurons.

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

Recent reports of a link between human- and primate-specific genetic factors and human/primate-specific characteristics and diseases necessitate genome-wide identification of those factors. We have previously reported core promoter short tandem repeats (STRs) of extreme length (≥6-repeats) that have expanded exceptionally in primates vs. non-primates, and may have a function in adaptive evolution. In the study reported here, we extended our study to the human STRs of ≥3-repeats in the category of penta and hexaucleotide STRs, across the entire human protein coding gene core promoters, and analyzed their status in several superorders and orders of vertebrates, using the Ensembl database. The ConSite software was used to identify the transcription factor (TF) sets binding to those STRs. STR specificity was observed at different levels of human and non-human primate (NHP) evolution. 73% of the pentanucleotide STRs and 68% of the hexanucleotide STRs were found to be specific to human and NHPs. AP-2alpha, Sp1, and MZF were the predominantly selected TFs (90%) binding to the human-specific STRs. Furthermore, the number of TF sets binding to a given STR was found to be a selection factor for that STR. Our findings indicate that selected STRs, the cognate binding TFs, and the number of TF set binding to those STRs function as switch codes at different levels of human and NHP evolution and speciation.

ZSCAN5B and primate-specific paralogs bind RNA polymerase III genes and extra-TFIIIC (ETC) sites to modulate mitotic progression.

Mammalian genomes contain hundreds of genes transcribed by RNA Polymerase III (Pol III), encoding noncoding RNAs and especially the tRNAs specialized to carry specific amino acids to the ribosome for protein synthesis. In addition to this well-known function, tRNAs and their genes (tDNAs) serve a variety of other critical cellular functions. For example, tRNAs and other Pol III transcripts can be cleaved to yield small RNAs with potent regulatory activities. Furthermore, from yeast to mammals, active tDNAs and related "extra-TFIIIC" (ETC) loci provide the DNA scaffolds for the most ancient known mechanism of three-dimensional chromatin architecture. Here we identify the ZSCAN5 TF family - including mammalian ZSCAN5B and its primate-specific paralogs - as proteins that occupy mammalian Pol III promoters and ETC sites. We show that ZSCAN5B binds with high specificity to a conserved subset of Pol III genes in human and mouse. Furthermore, primate-specific ZSCAN5A and ZSCAN5D also bind Pol III genes, although ZSCAN5D preferentially localizes to MIR SINE- and LINE2-associated ETC sites. ZSCAN5 genes are expressed in proliferating cell populations and are cell-cycle regulated, and siRNA knockdown experiments suggested a cooperative role in regulation of mitotic progression. Consistent with this prediction, ZSCAN5A knockdown led to increasing numbers of cells in mitosis and the appearance of cells. Together, these data implicate the role of ZSCAN5 genes in regulation of Pol III genes and nearby Pol II loci, ultimately influencing cell cycle progression and differentiation in a variety of tissues.