Primate-specific Alus harbor different regulatory features, including miRNA targets. In this study, we provide evidence for miRNA-mediated modulation of transcript isoform levels during heat-shock response through exaptation of Alu-miRNA sites in mature mRNA. We performed genome-wide expression profiling coupled with functional validation of miRNA target sites within exonized Alus, and analyzed conservation of these targets across primates. We observed that two miRNAs (miR-15a-3p and miR-302d-3p) elevated in stress response, target RAD1, GTSE1, NR2C1, FKBP9 and UBE2I exclusively within Alu. These genes map onto the p53 regulatory network. Ectopic overexpression of miR-15a-3p downregulates GTSE1 and RAD1 at the protein level and enhances cell survival. This Alu-mediated fine-tuning seems to be unique to humans as evident from the absence of orthologous sites in other primate lineages. We further analyzed signatures of selection on Alu-miRNA targets in the genome, using 1000 Genomes Phase-I data. We found that 198 out of 3177 Alu-exonized genes exhibit signatures of selection within Alu-miRNA sites, with 60 of them containing SNPs supported by multiple evidences (global-FST > 0.3, pair-wise-FST > 0.5, Fay-Wu’s H < −20, iHS > 2.0, high ΔDAF) and implicated in p53 network. We propose that by affecting multiple genes, Alu-miRNA interactions have the potential to facilitate population-level adaptations in response to environmental challenges.

Transforming growth factor beta3 (TGFB3) gene mutations in patients of arrhythmogenic right ventricular dysplasia/cardiomyopathy (ARVD1) and Loeys-Dietz syndrome-5 (LDS5)/Rienhoff syndrome are associated with cardiomyopathy, cardiac arrhythmia, cardiac fibrosis, cleft palate, aortic aneurysms, and valvular heart disease. Although the developing heart of embryos express Tgfb3, its overarching role remains unclear in cardiovascular development and disease. We used histological, immunohistochemical, and molecular analyses of Tgfb3−/− fetuses and compared them to wildtype littermate controls. The cardiovascular phenotypes were diverse with approximately two thirds of the Tgfb3−/− fetuses having one or more cardiovascular malformations, including abnormal ventricular myocardium (particularly of the right ventricle), outflow tract septal and alignment defects, abnormal aortic and pulmonary trunk walls, and thickening of semilunar and/or atrioventricular valves. Ventricular septal defects (VSD) including the perimembranous VSDs were observed in Tgfb3−/− fetuses with myocardial defects often accompanied by the muscular type VSD. In vitro studies using TGFβ3-deficient fibroblasts in 3-D collagen lattice formation assays indicated that TGFβ3 was required for collagen matrix reorganization. Biochemical studies indicated the ‘paradoxically’ increased activation of canonical (SMAD-dependent) and noncanonical (MAP kinase-dependent) pathways. TGFβ3 is required for cardiovascular development to maintain a balance of canonical and noncanonical TGFβ signaling pathways.

Alu repeats contribute to phylogenetic novelties in conserved regulatory networks in primates. Our study highlights how exonized Alus could nucleate large-scale mRNA-miRNA interactions. Using a functional genomics approach, we characterize a transcript isoform of an orphan gene, CYP20A1 (CYP20A1_Alu-LT) that has exonization of 23 Alus in its 3’UTR. CYP20A1_Alu-LT, confirmed by 3’RACE, is an outlier in length (9 kb 3’UTR) and widely expressed. Using publically available datasets, we demonstrate its expression in higher primates and presence in single nucleus RNA-seq of 15928 human cortical neurons. miRanda predicts ∼4700 miRNA recognition elements (MREs) for ∼1000 miRNAs, primarily originated within these 3’UTR-Alus. CYP20A1_Alu-LT could be a potential multi-miRNA sponge as it harbors ≥10 MREs for 140 miRNAs and has cytosolic localization. We further tested whether expression of CYP20A1_Alu-LT correlates with mRNAs harboring similar MRE targets. RNA-seq with conjoint miRNA-seq analysis was done in primary human neurons where we observed CYP20A1_Alu-LT to be downregulated during heat shock response and upregulated in HIV1-Tat treatment. 380 genes were positively correlated with its expression (significantly downregulated in heat shock and upregulated in Tat) and they harbored MREs for nine expressed miRNAs which were also enriched in CYP20A1_Alu-LT. MREs were significantly enriched in these 380 genes compared to random sets of differentially expressed genes (p = 8.134e-12). Gene ontology suggested involvement of these genes in neuronal development and hemostasis pathways thus proposing a novel component of Alu-miRNA mediated transcriptional modulation that could govern specific physiological outcomes in higher primates.