skip to main content


Title: miRNA analysis with Prost! reveals evolutionary conservation of organ-enriched expression and post-transcriptional modifications in three-spined stickleback and zebrafish
Abstract

MicroRNAs (miRNAs) can have organ-specific expression and functions; they can originate from dedicated miRNA genes, from non-canonical miRNA genes, or from mirror-miRNA genes and can also experience post-transcriptional variation. It remains unclear, however, which mechanisms of miRNA production or modification are organ-specific and the extent of their evolutionary conservation. To address these issues, we developed the softwareProst!(PRocessing Of Short Transcripts), which, among other features, helps quantify mature miRNAs, accounts for post-transcriptional processing, such as nucleotide editing, and identifies mirror-miRNAs. Here, we appliedProst!to annotate and analyze miRNAs in three-spined stickleback (Gasterosteus aculeatus), a model fish for evolutionary biology reported to have a miRNome larger than most teleost fish. Zebrafish (Danio rerio), a distantly related teleost with a well-known miRNome, served as comparator. Our results provided evidence for the existence of 286 miRNA genes and 382 unique mature miRNAs (excludingmir430gene duplicates and the vaultRNA-derivedmir733), which doesn’t represent a miRNAome larger than other teleost miRNomes. In addition, small RNA sequencing data from brain, heart, testis, and ovary in both stickleback and zebrafish identified suites of mature miRNAs that display organ-specific enrichment, many of which are evolutionarily-conserved in the brain and heart in both species. These data also supported the hypothesis that evolutionarily-conserved, organ-specific mechanisms may regulate post-transcriptional variations in miRNA sequence. In both stickleback and zebrafish, miR2188-5p was edited frequently with similar nucleotide changes in the seed sequence with organ specific editing rates, highest in the brain. In summary,Prost!is a new tool to identify and understand small RNAs, to help clarify a species’ miRNA biology as shown here for an important model for the evolution of developmental mechanisms, and to provide insight into organ-enriched expression and the evolutionary conservation of miRNA post-transcriptional modifications.

 
more » « less
NSF-PAR ID:
10153489
Author(s) / Creator(s):
; ; ; ;
Publisher / Repository:
Nature Publishing Group
Date Published:
Journal Name:
Scientific Reports
Volume:
9
Issue:
1
ISSN:
2045-2322
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. Abstract

    microRNAs (miRNAs) are post‐transcriptional regulators of gene expression and can play an important role in modulating organismal development and physiology in response to environmental stress. However, the role of miRNAs in mediating adaptation to diverse environments in natural study systems remains largely unexplored. Here, we characterized miRNAs and their expression inPoecilia mexicana, a species of small fish that inhabits both normal streams and extreme environments in the form of springs rich in toxic hydrogen sulphide (H2S). We found thatP. mexicanahas a similar number of miRNA genes as other teleosts. In addition, we identified a large population of mature miRNAs that were differentially expressed between locally adapted populations in contrasting habitats, indicating that miRNAs may contribute toP. mexicanaadaptation to sulphidic environments. In silico identification of differentially expressed miRNA‐mRNA pairs revealed, in the sulphidic environment, the downregulation of miRNAs predicted to target mRNAs involved in sulphide detoxification and cellular homeostasis, which are pathways essential for life in H2S‐rich springs. In addition, we found that predicted targets of upregulated miRNAs act in the mitochondria (16.6% of predicted annotated targets), which is the main site of H2S toxicity and detoxification, possibly modulating mitochondrial function. Together, the differential regulation of miRNAs between these natural populations suggests that miRNAs may be involved in H2S adaptation by promoting functions needed for survival and reducing functions affected by H2S. This study lays the groundwork for further research to directly demonstrate the role of miRNAs in adaptation to H2S. Overall, this study provides a critical stepping‐stone towards a comprehensive understanding of the regulatory mechanisms underlying the adaptive variation in gene expression in a natural system.

     
    more » « less
  2. Social challenges like territorial intrusions evoke behavioral responses in widely diverging species. Recent work has showed that evolutionary “toolkits”—genes and modules with lineage‐specific variations but deep conservation of function—participate in the behavioral response to social challenge. Here, we develop a multispecies computational‐experimental approach to characterize such a toolkit at a systems level. Brain transcriptomic responses to social challenge was probed via RNA‐seq profiling in three diverged species—honey bees, mice and three‐spined stickleback fish—following a common methodology, allowing fair comparisons across species. Data were collected from multiple brain regions and multiple time points after social challenge exposure, achieving anatomical and temporal resolution substantially greater than previous work. We developed statistically rigorous analyses equipped to find homologous functional groups among these species at the levels of individual genes, functional and coexpressed gene modules, and transcription factor subnetworks. We identified six orthogroups involved in response to social challenge, including groups represented by mouse genesNpas4andNr4a1, as well as common modulation of systems such as transcriptional regulators, ion channels, G‐protein‐coupled receptors and synaptic proteins. We also identified conserved coexpression modules enriched for mitochondrial fatty acid metabolism and heat shock that constitute the shared neurogenomic response. Our analysis suggests a toolkit wherein nuclear receptors, interacting with chaperones, induce transcriptional changes in mitochondrial activity, neural cytoarchitecture and synaptic transmission after social challenge. It shows systems‐level mechanisms that have been repeatedly co‐opted during evolution of analogous behaviors, thus advancing the genetic toolkit concept beyond individual genes.

     
    more » « less
  3. Norwick, Katja (Ed.)
    Abstract In plants, miRNA production is orchestrated by a suite of proteins that control transcription of the pri-miRNA gene, post-transcriptional processing and nuclear export of the mature miRNA. Post-transcriptional processing of miRNAs is controlled by a pair of physically interacting proteins, hyponastic leaves 1 (HYL1) and Dicer-like 1 (DCL1). However, the evolutionary history and structural basis of the HYL1–DCL1 interaction is unknown. Here we use ancestral sequence reconstruction and functional characterization of ancestral HYL1 in vitro and in Arabidopsis thaliana to better understand the origin and evolution of the HYL1–DCL1 interaction and its impact on miRNA production and plant development. We found the ancestral plant HYL1 evolved high affinity for both double-stranded RNA (dsRNA) and its DCL1 partner before the divergence of mosses from seed plants (∼500 Ma), and these high-affinity interactions remained largely conserved throughout plant evolutionary history. Structural modeling and molecular binding experiments suggest that the second of two dsRNA-binding motifs (DSRMs) in HYL1 may interact tightly with the first of two C-terminal DCL1 DSRMs to mediate the HYL1–DCL1 physical interaction necessary for efficient miRNA production. Transgenic expression of the nearly 200 Ma-old ancestral flowering-plant HYL1 in A. thaliana was sufficient to rescue many key aspects of plant development disrupted by HYL1− knockout and restored near-native miRNA production, suggesting that the functional partnership of HYL1–DCL1 originated very early in and was strongly conserved throughout the evolutionary history of terrestrial plants. Overall, our results are consistent with a model in which miRNA-based gene regulation evolved as part of a conserved plant “developmental toolkit.” 
    more » « less
  4. Abstract Aims

    Dissecting complex interactions among transcription factors (TFs), microRNAs (miRNAs) and long noncoding RNAs (lncRNAs) are central for understanding heart development and function. Although computational approaches and platforms have been described to infer relationships among regulatory factors and genes, current approaches do not adequately account for how highly diverse, interacting regulators that include noncoding RNAs (ncRNAs) control cardiac gene expression dynamics over time.

    Methods

    To overcome this limitation, we devised an integrated framework, cardiac gene regulatory modeling (CGRM) that integrates LogicTRN and regulatory component analysis bioinformatics modeling platforms to infer complex regulatory mechanisms. We then used CGRM to identify and compare the TF-ncRNA gene regulatory networks that govern early- and late-stage cardiomyocytes (CMs) generated by in vitro differentiation of human pluripotent stem cells (hPSC) and ventricular and atrial CMs isolated during in vivo human cardiac development.

    Results

    Comparisons of in vitro versus in vivo derived CMs revealed conserved regulatory networks among TFs and ncRNAs in early cells that significantly diverged in late staged cells. We report that cardiac genes (“heart targets”) expressed in early-stage hPSC-CMs are primarily regulated by MESP1, miR-1, miR-23, lncRNAs NEAT1 and MALAT1, while GATA6, HAND2, miR-200c, NEAT1 and MALAT1 are critical for late hPSC-CMs. The inferred TF-miRNA-lncRNA networks regulating heart development and contraction were similar among early-stage CMs, among individual hPSC-CM datasets and between in vitro and in vivo samples. However, genes related to apoptosis, cell cycle and proliferation, and transmembrane transport showed a high degree of divergence between in vitro and in vivo derived late-stage CMs. Overall, late-, but not early-stage CMs diverged greatly in the expression of “heart target” transcripts and their regulatory mechanisms.

    Conclusions

    In conclusion, we find that hPSC-CMs are regulated in a cell autonomous manner during early development that diverges significantly as a function of time when compared to in vivo derived CMs. These findings demonstrate the feasibility of using CGRM to reveal dynamic and complex transcriptional and posttranscriptional regulatory interactions that underlie cell directed versus environment-dependent CM development. These results with in vitro versus in vivo derived CMs thus establish this approach for detailed analyses of heart disease and for the analysis of cell regulatory systems in other biomedical fields.

     
    more » « less
  5. Abstract Background

    Sexually dimorphic mating behaviors differ between sexes and involve gonadal hormones and possibly sexually dimorphic gene expression in the brain. However, the associations among the brain, gonad, and sexual behavior in teleosts are still unclear. Here, we utilized germ cells-freetdrd12knockout (KO) zebrafish, and steroid synthesis enzymecyp17a1-deficient zebrafish to investigate the differences and interplays in the brain–gonad–behavior axis, and the molecular control of brain dimorphism and male mating behaviors.

    Methods

    Tdrd12+/−;cyp17a1+/−double heterozygous parents were crossed to obtaintdrd12−/−;cyp17a1+/+(tdrd12 KO),tdrd12+/+;cyp17a1−/−(cyp17a1 KO), andtdrd12−/−;cyp17a1−/−(double KO) homozygous progenies. Comparative analysis of mating behaviors were evaluated using Viewpoint zebrafish tracking software and sexual traits were thoroughly characterized based on anatomical and histological experiments in these KOs and wild types. The steroid hormone levels (testosterone, 11-ketotestosterone and 17β-estradiol) in the brains, gonads, and serum were measured using ELISA kits. To achieve a higher resolution view of the differences in region-specific expression patterns of the brain, the brains of these KOs, and control male and female fish were dissected into three regions: the forebrain, midbrain, and hindbrain for transcriptomic analysis.

    Results

    Qualitative analysis of mating behaviors demonstrated thattdrd12−/−fish behaved in the same manner as wild-type males to trigger oviposition behavior, whilecyp17a1−/−and double knockout (KO) fish did not exhibit these behaviors. Based on the observation of sex characteristics, mating behaviors and hormone levels in these mutants, we found that the maintenance of secondary sex characteristics and male mating behavior did not depend on the presence of germ cells; rather, they depended mainly on the 11-ketotestosterone and testosterone levels secreted into the brain–gonad regulatory axis. RNA-seq analysis of different brain regions revealed that the brain transcript profile oftdrd12−/−fish was similar to that of wild-type males, especially in the forebrain and midbrain. However, the brain transcript profiles ofcyp17a1−/−and double KO fish were distinct from those of wild-type males and were partially biased towards the expression pattern of the female brain. Our results revealed important candidate genes and signaling pathways, such as synaptic signaling/neurotransmission, MAPK signaling, and steroid hormone pathways, that shape brain dimorphism and modulate male mating behavior in zebrafish.

    Conclusions

    Our results provide comprehensive analyses and new insights regarding the endogenous interactions in the brain–gonad–behavior axis. Moreover, this study revealed the crucial candidate genes and neural signaling pathways of different brain regions that are involved in modulating brain dimorphism and male mating behavior in zebrafish, which would significantly light up the understanding the neuroendocrine and molecular mechanisms modulating brain dimorphism and male mating behavior in zebrafish and other teleost fish.

    Graphical Abstract 
    more » « less