Clusters of microRNAs emerge by new hairpins in existing transcripts

Plant bioactive compounds driven microRNAs (miRNAs): A potential source and novel strategy targeting gene and cancer therapeutics

Irrespective of medical technology improvements, cancer ranks among the leading causes of mortality worldwide. Although numerous cures and treatments exist, creating alternative cancer therapies with fewer adverse side effects is vital. Since ancient times, plant bioactive compounds have already been used as a remedy to heal cancer. These plant bioactive compounds and their anticancer activity can also deregulate the microRNAs (miRNAs) in the cancerous cells. Therefore, the deregulation of miRNAs in cancer cells by plant bioactive compounds and the usage of the related miRNA could be a promising approach for cancer cure, mainly to prevent cancer and overcome chemotherapeutic side effect problems. Hence, this review highlights the function of plant bioactive compounds as an anticancer agent through the underlying mechanism that alters the miRNA expression in cancer cells, ultimately leading to apoptosis. Moreover, this review provides insight into using plant bioactive compounds -driven miRNAs as an anticancer agent to develop miRNA-based cancer gene therapy. They can be the potential resource for gene therapy and novel strategies targeting cancer therapeutics.

Genome of tropical bed bug Cimex hemipterus (Cimicidae, Hemiptera) reveals tetraspanin expanded in bed bug ancestor

Cimex species are ectoparasites that exclusively feed on warm-blooded animals such as birds and mammals. Three cimicid species are known to be persistent pests for humans, including the tropical bed bug Cimex hemipterus, common bed bug Cimex lectularius, and Eastern bat bug Leptocimex boueti. To date, genomic information is restricted to the common bed bug C. lectularius, which limits understanding their biology and to provide controls of bed bug infestations. Here, a chromosomal-level genome assembly of C. hemipterus (495 Mb [megabase pairs]) contained on 16 pseudochromosomes (scaffold N50 = 34 Mb), together with 9 messenger RNA and small RNA transcriptomes were obtained. In comparison between hemipteran genomes, we found that the tetraspanin superfamily was expanded in the Cimex ancestor. This study provides the first genome assembly for the tropical bed bug C. hemipterus, and offers an unprecedented opportunity to address questions relating to bed bug infestations, as well as genomic evolution to hemipterans more widely.

Convergent and divergent evolution of microRNA-mediated regulation in metazoans

The evolution of microRNAs (miRNAs) has been studied extensively to understand their roles in gene regulation and evolutionary processes. This review focuses on how miRNA-mediated regulation has evolved in bilaterian animals, highlighting both convergent and divergent evolution. Since animals and plants display significant differences in miRNA biogenesis and target recognition, the 'independent origin' hypothesis proposes that miRNA pathways in these groups independently evolved from the RNA interference (RNAi) pathway, leading to modern miRNA repertoires through convergent evolution. However, recent evidence raises the alternative possibility that the miRNA pathway might have already existed in the last common ancestor of eukaryotes, and that the differences in miRNA pathway and miRNA repertoires among animal and plant lineages arise from lineage-specific innovations and losses of miRNA pathways, miRNA acquisition, and loss of miRNAs after eukaryotic divergence. The repertoire of miRNAs has considerably expanded during bilaterian evolution, primarily through de novo creation and duplication processes, generating new miRNAs. Although ancient functionally established miRNAs are rarely lost, many newly emerged miRNAs are transient and lineage specific, following a birth-death evolutionary pattern aligning with the 'out-of-the-testis' and 'transcriptional control' hypotheses. Our focus then shifts to the convergent molecular evolution of miRNAs. We summarize how miRNA clustering and seed mimicry contribute to this phenomenon, and we review how miRNAs from different sources converge to degrade maternal messenger RNAs (mRNAs) during animal development. Additionally, we describe how miRNAs evolve across species due to changes in sequence, seed shifting, arm switching, and spatiotemporal expression patterns, which can result in variations in target sites among orthologous miRNAs across distant strains or species. We also provide a summary of the current understanding regarding how the target sites of orthologous miRNAs can vary across strains or distantly related species. Although many paralogous miRNAs retain their seed or mature sequences after duplication, alterations can occur in the seed or mature sequences or expression patterns of paralogous miRNAs, leading to functional diversification. We discuss our current understanding of the functional divergence between duplicated miRNAs, and illustrate how the functional diversification of duplicated miRNAs impacts target site evolution. By investigating these topics, we aim to enhance our current understanding of the functions and evolutionary dynamics of miRNAs. Additionally, we shed light on the existing challenges in miRNA evolutionary studies, particularly the complexity of deciphering the role of miRNA-mediated regulatory network evolution in shaping gene expression divergence and phenotypic differences among species.

Parameters of clustered suboptimal miRNA biogenesis

The nuclear cleavage of a suboptimal primary miRNA hairpin by the Drosha/DGCR8 complex ("Microprocessor") can be enhanced by an optimal miRNA neighbor, a phenomenon termed cluster assistance. Several features and biological impacts of this new layer of miRNA regulation are not fully known. Here, we elucidate the parameters of cluster assistance of a suboptimal miRNA and also reveal competitive interactions amongst optimal miRNAs within a cluster. We exploit cluster assistance as a functional assay for suboptimal processing and use this to invalidate putative suboptimal substrates, as well as identify a "solo" suboptimal miRNA. Finally, we report complexity in how specific mutations might affect the biogenesis of clustered miRNAs in disease contexts. This includes how an operon context can buffer the effect of a deleterious processing variant, but reciprocally how a point mutation can have a nonautonomous effect to impair the biogenesis of a clustered, suboptimal, neighbor. These data expand our knowledge regarding regulated miRNA biogenesis in humans and represent a functional assay for empirical definition of suboptimal Microprocessor substrates.

ncOrtho: efficient and reliable identification of miRNA orthologs

MicroRNAs (miRNAs) are post-transcriptional regulators that finetune gene expression via translational repression or degradation of their target mRNAs. Despite their functional relevance, frameworks for the scalable and accurate detection of miRNA orthologs are missing. Consequently, there is still no comprehensive picture of how miRNAs and their associated regulatory networks have evolved. Here we present ncOrtho, a synteny informed pipeline for the targeted search of miRNA orthologs in unannotated genome sequences. ncOrtho matches miRNA annotations from multi-tissue transcriptomes in precision, while scaling to the analysis of hundreds of custom-selected species. The presence-absence pattern of orthologs to 266 human miRNA families across 402 vertebrate species reveals four bursts of miRNA acquisition, of which the most recent event occurred in the last common ancestor of higher primates. miRNA families are rarely modified or lost, but notable exceptions for both events exist. miRNA co-ortholog numbers faithfully indicate lineage-specific whole genome duplications, and miRNAs are powerful markers for phylogenomic analyses. Their exceptionally low genetic diversity makes them suitable to resolve clades where the phylogenetic signal is blurred by incomplete lineage sorting of ancestral alleles. In summary, ncOrtho allows to routinely consider miRNAs in evolutionary analyses that were thus far reserved to protein-coding genes.

microRNA Expression Dynamics in Culicoides sonorensis Biting Midges Following Blood-Feeding

Culicoides sonorensis midges vector multiple livestock arboviruses, resulting in significant economic losses worldwide. Due to the tight association between virus transmission, blood feeding, and egg development, understanding midge physiology is paramount to limiting pathogen transmission. Previous studies have demonstrated the importance of small non-coding RNAs (ncRNAs), specifically microRNAs (miRNAs), in multiple aspects of vector physiology. These small ncRNAs regulate gene expression at the post-transcriptional level and display differential expression during pathogen infection. Due to the lack of annotated miRNAs in the biting midge and associated expression profiles, we used small RNA-Seq and miRDeep2 analyses to determine the Culicoides miRNAs in whole females and midgut tissues in response to blood feeding. Our analyses revealed 76 miRNAs within C. sonorensis composed of 73 orthologous and three candidate novel miRNAs, as well as conserved miRNA clusters. miRNA conservation suggests an interesting evolutionary relationship between miRNA expression and hematophagy in the infraorder Culicomorpha. We also identified multiple blood meal-regulated and tissue-enriched miRNAs. Lastly, we further identified miRNAs with expression patterns potentially associated with virus infection by probing publicly available datasets. Together, our data provide a foundation for future ncRNA work to untangle the dynamics of gene regulation associated with midge physiology.

A genome-wide computational approach to define microRNA-Polycomb/trithorax gene regulatory circuits in Drosophila

Characterization of gene regulatory networks is fundamental to understanding homeostatic development. This process can be simplified by analyzing relatively simple genomes such as the genome of Drosophila melanogaster. In this work we have developed a computational framework in Drosophila to explore for the presence of gene regulatory circuits between two large groups of transcriptional regulators: the epigenetic group of the Polycomb/trithorax (PcG/trxG) proteins and the microRNAs (miRNAs). We have searched genome-wide for miRNA targets in PcG/trxG transcripts as well as for Polycomb Response Elements (PREs) in miRNA genes. Our results show that 10% of the analyzed miRNAs could be controlling PcG/trxG gene expression, while 40% of those miRNAs are putatively controlled by the selected set of PcG/trxG proteins. The integration of these analyses has resulted in the predicted existence of 3 classes of miRNA-PcG/trxG crosstalk interactions that define potential regulatory circuits. In the first class, miRNA-PcG circuits are defined by miRNAs that reciprocally crosstalk with PcG. In the second, miRNA-trxG circuits are defined by miRNAs that reciprocally crosstalk with trxG. In the third class, miRNA-PcG/trxG shared circuits are defined by miRNAs that crosstalk with both PcG and trxG regulators. These putative regulatory circuits may uncover a novel mechanism in Drosophila for the control of PcG/trxG and miRNAs levels of expression. The computational framework developed here for Drosophila melanogaster can serve as a model case for similar analyses in other species. Moreover, our work provides, for the first time, a new and useful resource for the Drosophila community to consult prior to experimental studies investigating the epigenetic regulatory networks of miRNA-PcG/trxG mediated gene expression.

Clustering pattern and evolution characteristic of microRNAs in grass carp (Ctenopharyngodon idella)

BACKGROUND

A considerable fraction of microRNAs (miRNAs) are highly conserved, and certain miRNAs correspond to genomic clusters. The clustering of miRNAs can be advantageous, possibly by allowing coordinated expression. However, little is known about the evolutionary forces responsible for the loss and acquisition of miRNA and miRNA clusters.

RESULTS

The results demonstrated that several novel miRNAs arose throughout grass carp evolution. Duplication and de novo production were critical strategies for miRNA cluster formation. Duplicates accounted for a smaller fraction of the expansion in the grass carp miRNA than de novo creation. Clustered miRNAs are more conserved and change slower, whereas unique miRNAs usually have high evolution rates and low expression levels. The expression level of miRNA expression in clusters is strongly correlated.

CONCLUSIONS

This study examines the genomic distribution, evolutionary background, and expression regulation of grass carp miRNAs. Our findings provide novel insights into the genesis and development of miRNA clusters in teleost.

The chromosomal distribution of sex-biased microRNAs in Drosophila is non-adaptive

Genes are often differentially expressed between males and females. In Drosophila melanogaster, the analysis of sex-biased microRNAs (short noncoding regulatory molecules) has revealed striking differences with protein-coding genes. Mainly, the X chromosome is enriched in male-biased microRNA genes, although it is depleted of male-biased protein-coding genes. The paucity of male-biased genes in the X chromosome is generally explained by an evolutionary process called demasculinization. I suggest that the excess of male-biased microRNAs in the X chromosome is due to high rates of de novo emergence of microRNAs (mostly in other neighboring microRNAs), a tendency of novel microRNAs in the X chromosome to be expressed in testis, and to a lack of a demasculinization process. To test this hypothesis, I analyzed the expression profile of microRNAs in males, females, and gonads in D. pseudoobscura, in which an autosome translocated into the X chromosome effectively becoming part of a sex chromosome (neo-X). I found that the pattern of sex-biased expression is generally conserved between D. melanogaster and D. pseudoobscura. Also, orthologous microRNAs in both species conserve their chromosomal location, indicating that there is no evidence of demasculinization or other interchromosomal movement of microRNAs. Drosophila pseudoobscura-specific microRNAs in the neo-X chromosome tend to be male-biased and particularly expressed in testis. In summary, the apparent paradox resulting from male-biased protein-coding genes depleted in the X chromosome and an enrichment in male-biased microRNAs is consistent with different evolutionary dynamics between coding genes and short RNAs.

A Rapid Evolving microRNA Cluster Rewires Its Target Regulatory Networks in Drosophila

New miRNAs are evolutionarily important but their functional evolution remains unclear. Here we report that the evolution of a microRNA cluster, mir-972C rewires its downstream regulatory networks in Drosophila. Genomic analysis reveals that mir-972C originated in the common ancestor of Drosophila where it comprises six old miRNAs. It has subsequently recruited six new members in the melanogaster subgroup after evolving for at least 50 million years. Both the young and the old mir-972C members evolved rapidly in seed and non-seed regions. Combining target prediction and cell transfection experiments, we found that the seed and non-seed changes in individual mir-972C members cause extensive target divergence among D. melanogaster, D. simulans, and D. virilis, consistent with the functional evolution of mir-972C reported recently. Intriguingly, the target pool of the cluster as a whole remains relatively conserved. Our results suggest that clustering of young and old miRNAs broadens the target repertoires by acquiring new targets without losing many old ones. This may facilitate the establishment of new miRNAs in existing regulatory networks.

Diagnostic and Prognostic Potential of MiR-379/656 MicroRNA Cluster in Molecular Subtypes of Breast Cancer

Introduction: Breast cancer is the most frequently diagnosed cancer globally and is one of the most important contributors to cancer-related deaths. Earlier diagnosis is known to reduce mortality, and better biomarkers are needed. MiRNA clusters often co-express and target mRNAs in a coordinated fashion, perturbing entire pathways; they thus merit further exploration for diagnostic or prognostic use. MiR-379/656, at chromosome 14q32, is the second largest miRNA cluster in the human genome and implicated in various malignancies including glioblastoma, melanoma, gastrointestinal tumors and ovarian cancer highlighting its potential importance. In this study, we focus on the diagnostic and prognostic potentials of MiR-379/656 in breast cancer and its molecular subtypes. Materials and Methods: We analyzed miRNA and mRNA next generation sequencing data from 903 primary tumors and 90 normal controls (source: The Cancer Genome Atlas). The differential expression profile between tumor and normal was analyzed using DeSEQ2. Penalized logistic regression modelling (lasso regression) was used to assess the predictive potential of MiR-379/656 expression for tumor and normal samples. The association between MiR-379/656 expression and overall patient survival was studied using Cox Proportional-Hazard Model. The target mRNAs (validated) of MiR-379/656 were annotated via pathway enrichment analysis to understand the biological significance of the cluster in breast cancer. Results: The differential expression analysis for 1390 miRNAs (miRnome) revealed 310 upregulated (22.3%) and 176 downregulated (12.66%) miRNAs in breast cancer patients compared with controls. For MiR-379/656, 32 miRNAs (32/42; 76%) were downregulated. The MiR-379/656 cluster was found to be the most differentially expressed cluster in the human genome (p < 10⁻³⁰). The Basal and Luminal B subtypes showed at least 83% (35/42) of the miRNAs to be downregulated. The binomial model prioritized 15 miRNAs, which distinguished breast cancer patients from controls with 99.15 ± 0.58% sensitivity and 77.78 ± 5.24% specificity. Overall, the Basal and Luminal B showed the most effective predictive power with respect to the 15 prioritized miRNAs at MiR-379/656 cluster. The decreased expression of MiR-379/656 was found to be associated with poorer clinical outcome in Basal and Luminal B subtypes, increasing tumor stage and tumor size/extent, and overall patient survival. Pathway enrichment for the validated targets of MiR-379/656 was significant for cancer-related pathways, especially DNA repair, transcriptional regulation by p53 and cell cycle checkpoints (adjusted p-value < 0.05). Conclusions: Genome informatics analysis of high throughput data for MiR-379/656 cluster has shown that a subset of 15 miRNAs from MiR-379/656 cluster can be used for the diagnostic and prognostic purpose of breast cancer and its subtypes—especially in Basal and Luminal B.

MicroRNAs and osteocytes

MicroRNAs, identified in the early 1990s, are believed to regulate approximately 30% of the human genome. The role of microRNA in bone cells was first reported in 2007 in a manuscript showing that microRNA-223 is essential for osteoclast differentiation in vitro, and a few studies reported a role of microRNAs in osteoblasts the same year. The first report of microRNA actions in osteocytes was published in 2010, in which it was demonstrated that the microRNA cluster 23a~27a~24-2 regulates osteocyte differentiation. Since then, few studies have described the role of these 18-25-nucleotide non-coding RNAs on osteocyte biology, reporting osteocytes both as producers and as targets of the actions of microRNAs. We review here the current knowledge on the effects of microRNAs on osteocyte biology.

Evolution and Phylogeny of MicroRNAs - Protocols, Pitfalls, and Problems

MicroRNAs are important regulators in many eukaryotic lineages. Typical miRNAs have a length of about 22nt and are processed from precursors that form a characteristic hairpin structure. Once they appear in a genome, miRNAs are among the best-conserved elements in both animal and plant genomes. Functionally, they play an important role in particular in development. In contrast to protein-coding genes, miRNAs frequently emerge de novo. The genomes of animals and plants harbor hundreds of mutually unrelated families of homologous miRNAs that tend to be persistent throughout evolution. The evolution of their genomic miRNA complement closely correlates with important morphological innovation. In addition, miRNAs have been used as valuable characters in phylogenetic studies. An accurate and comprehensive annotation of miRNAs is required as a basis to understand their impact on phenotypic evolution. Since experimental data on miRNA expression are limited to relatively few species and are subject to unavoidable ascertainment biases, it is inevitable to complement miRNA sequencing by homology based annotation methods. This chapter reviews the state of the art workflows for homology based miRNA annotation, with an emphasis on their limitations and open problems.

Clustering genomic organization of sea cucumber miRNAs impacts their evolution and expression

Echinoderms are marine deuterostomes with fascinating adaptation features such as aestivation and organ regeneration. However, post-transcriptional gene regulation by microRNAs (miRNAs) underlying these features are largely unexplored. Here, using homology-based and de novo approaches supported by expression data, we provided a comprehensive annotation of miRNA genes in the sea cucumber Apostichopus japonicus. By linkage and phylogenic analyses, we characterized miRNA genomic organization, evolutionary history and expression regulation. The results showed that sea cucumbers evolved a large number of new miRNAs, which tended to form polycistronic clusters via tandem duplication that had been especially active in the echinoderms. Most new miRNAs were weakly expressed, but miRNA clustering increased the expression level of clustered new miRNAs. The most abundantly expressed new miRNAs were organized in a single tandem cluster (cluster n2), which was activated during aestivation and intestine regeneration. Overall, our analyses suggest that clustering of miRNAs is important for their evolutionary origin, expression control, and functional cooperation.

Evolution after Whole-Genome Duplication: Teleost MicroRNAs

MicroRNAs (miRNAs) are important gene expression regulators implicated in many biological processes, but we lack a global understanding of how miRNA genes evolve and contribute to developmental canalization and phenotypic diversification. Whole-genome duplication events likely provide a substrate for species divergence and phenotypic change by increasing gene numbers and relaxing evolutionary pressures. To understand the consequences of genome duplication on miRNA evolution, we studied miRNA genes following the teleost genome duplication (TGD). Analysis of miRNA genes in four teleosts and in spotted gar, whose lineage diverged before the TGD, revealed that miRNA genes were retained in ohnologous pairs more frequently than protein-coding genes, and that gene losses occurred rapidly after the TGD. Genomic context influenced retention rates, with clustered miRNA genes retained more often than nonclustered miRNA genes and intergenic miRNA genes retained more frequently than intragenic miRNA genes, which often shared the evolutionary fate of their protein-coding host. Expression analyses revealed both conserved and divergent expression patterns across species in line with miRNA functions in phenotypic canalization and diversification, respectively. Finally, major strands of miRNA genes experienced stronger purifying selection, especially in their seeds and 3'-complementary regions, compared with minor strands, which nonetheless also displayed evolutionary features compatible with constrained function. This study provides the first genome-wide, multispecies analysis of the mechanisms influencing metazoan miRNA evolution after whole-genome duplication.

A genome-wide microRNA screen identifies the microRNA-183/96/182 cluster as a modulator of circadian rhythms

The regulatory mechanisms of circadian rhythms have been studied primarily at the level of the transcription-translation feedback loops of protein-coding genes. Regulatory modules involving noncoding RNAs are less thoroughly understood. In particular, emerging evidence has revealed the important role of microRNAs (miRNAs) in maintaining the robustness of the circadian system. To identify miRNAs that have the potential to modulate circadian rhythms, we conducted a genome-wide miRNA screen using U2OS luciferase reporter cells. Among 989 miRNAs in the library, 120 changed the period length in a dose-dependent manner. We further validated the circadian regulatory function of an miRNA cluster, miR-183/96/182, both in vitro and in vivo. We found that all three members of this miRNA cluster can modulate circadian rhythms. Particularly, miR-96 directly targeted a core circadian clock gene, PER2. The knockout of the miR-183/96/182 cluster in mice showed tissue-specific effects on circadian parameters and altered circadian rhythms at the behavioral level. This study identified a large number of miRNAs, including the miR-183/96/182 cluster, as circadian modulators. We provide a resource for further understanding the role of miRNAs in the circadian network and highlight the importance of miRNAs as a genome-wide layer of circadian clock regulation.

Differential expression of microRNAs and tRNA fragments mediate the adaptation of the liver fluke Fasciola gigantica to its intermediate snail and definitive mammalian hosts

The tropical liver fluke Fasciola gigantica affects livestock and humans in many Asian countries, large parts of Africa, and parts of Europe. Despite the public health and economic impacts of F. gigantica, understanding of F. gigantica biology and how the complex lifecycle of this liver fluke is transcriptionally regulated remain unknown. Here, we tested the hypothesis that the regulatory small non-coding RNAs (sncRNAs), microRNAs (miRNAs) and tRNA-derived fragments (tRFs) play roles in the adaptation of F. gigantica to its intermediate and definitive hosts. We sequenced sncRNAs of eight lifecycle stages of F. gigantica. In total, 56 miRNAs from 33 conserved families and four Fasciola-specific miRNAs were identified. Expression analysis of miRNAs suggested clear stage-related patterns. By leveraging the existing transcriptomic data, we predicted a miRNA-based regulation of metabolism, transport, growth and developmental processes. Also, by comparing miRNA complement of F. gigantica with that of Fasciola hepatica, we detected a high level of conservation and identified differences in some miRNAs, which can be used to distinguish the two species. Moreover, we found that tRFs at each lifecycle stage were predominantly derived by tRNA^-Lys and tRNA^-Gly at 5' half sites, but relatively high expression was related to the buffalo-infecting stages. Taken together, we provided a comprehensive overview of the dynamic transcriptional changes of small RNAs that occur during the developmental stages of F. gigantica. This global analysis of F. gigantica lifecycle stages revealed new roles of miRNAs and tRFs in parasite development and will facilitate future research into understanding of fasciolosis pathobiology.

Discovery and functional understanding of MiRNAs in molluscs: a genome-wide profiling approach

Small non-coding RNAs play a pivotal role in gene regulation, repression of transposable element and viral activity in various organisms. Among the various categories of these small non-coding RNAs, microRNAs (miRNAs) guide post-translational gene regulation in cellular development, proliferation, apoptosis, oncogenesis, and differentiation. Here, we performed a genome-wide computational prediction of miRNAs to improve the understanding of miRNA observation and function in molluscs. As an initial step, hundreds of conserved miRNAs were predicted in 35 species of molluscs through genome scanning. Afterwards, the miRNAs' population, isoforms, organization, and function were characterized in detail. Furthermore, the key miRNA biogenesis factors, including AGO2, DGCR8, DICER, DROSHA, TRABP2, RAN, and XPO5, were elucidated based on homologue sequence searching. We also summarized the miRNAs' function in biomineralization, immune and stress response, as well as growth and development in molluscs. Because miRNAs play a vital role in various lifeforms, this study will provide insight into miRNA biogenesis and function in molluscs, as well as other invertebrates.

Micro-RNA Clusters Integrate Evolutionary Constraints on Expression and Target Affinities: The miR-6/5/4/286/3/309 Cluster in Drosophila

A striking feature of micro-RNAs is that they are often clustered in the genomes of animals. The functional and evolutionary consequences of this clustering remain obscure. Here, we investigated a micro-RNA cluster miR-6/5/4/286/3/309 that is conserved across drosophilid lineages. Small RNA sequencing revealed expression of this micro-RNA cluster in Drosophila melanogaster leg discs, and conditional overexpression of the whole cluster resulted in leg appendage shortening. Transgenic overexpression lines expressing different combinations of micro-RNA cluster members were also constructed. Expression of individual micro-RNAs from the cluster resulted in a normal wild-type phenotype, but either the expression of several ancient micro-RNAs together (miR-5/4/286/3/309) or more recently evolved clustered micro-RNAs (miR-6-1/2/3) can recapitulate the phenotypes generated by the whole-cluster overexpression. Screening of transgenic fly lines revealed downregulation of leg-patterning gene cassettes in generation of the leg-shortening phenotype. Furthermore, cell transfection with different combinations of micro-RNA cluster members revealed a suite of downstream genes targeted by all cluster members, as well as complements of targets that are unique for distinct micro-RNAs. Considered together, the micro-RNA targets and the evolutionary ages of each micro-RNA in the cluster demonstrate the importance of micro-RNA clustering, where new members can reinforce and modify the selection forces on both the cluster regulation and the gene regulatory network of existing micro-RNAs. Key words: micro-RNA, cluster, evolution.

Analysis of conserved miRNAs in cynomolgus macaque genome using small RNA sequencing and homology searching

MicroRNAs (miRNAs) are important regulators that fine-tune diverse cellular activities. Cynomolgus macaques (Macaca fascicularis) are used extensively in biomedical and pharmaceutical research; however, substantially fewer miRNAs have been identified in this species than in humans. Consequently, we investigated conserved miRNA profiles in cynomolgus macaques by homology searching and small RNA sequencing. In total, 1,455 high-confidence miRNA gene loci were identified, 408 of which were also confirmed by RNA sequencing, including 73 new miRNA loci reported in cynomolgus macaques for the first time. Comparing miRNA expression with age, we found a positive correlation between sequence conservation and expression levels during miRNA evolution. Additionally, we found that the miRNA gene locations in cynomolgus macaque genome were very flexible. Most were embedded in intergenic spaces or introns and clustered together. Several miRNAs were found in certain gene locations, including 64 exon-resident miRNAs, six splice-site-overlapping miRNAs (SO-miRNAs), and two pairs of distinct mirror miRNAs. We also identified 78 miRNA clusters, 68 of which were conserved in the human genome, including 10 large miRNA clusters predicted to regulate diverse developmental and cellular processes in cynomolgus macaque. Thus, this study not only expands the number of identified miRNAs in cynomolgus macaques but also provides clues for future research on the differences in miRNA repertoire between macaques and humans.

Genomic Clustering Facilitates Nuclear Processing of Suboptimal Pri-miRNA Loci

Nuclear processing of most miRNAs is mediated by Microprocessor, comprised of RNase III enzyme Drosha and its cofactor DGCR8. Here, we uncover a hidden layer of Microprocessor regulation via studies of Dicer-independent mir-451, which is clustered with canonical mir-144. Although mir-451 is fully dependent on Drosha/DGCR8, its short stem and small terminal loop render it an intrinsically weak Microprocessor substrate. Thus, it must reside within a cluster for normal biogenesis, although the identity and orientation of its neighbor are flexible. We use DGCR8 tethering assays and operon structure-function assays to demonstrate that local recruitment and transfer of Microprocessor enhances suboptimal substrate processing. This principle applies more broadly since genomic analysis indicates suboptimal canonical miRNAs are enriched in operons, and we validate several of these experimentally. Proximity-based enhancement of suboptimal hairpin processing provides a rationale for genomic retention of certain miRNA operons and may explain preferential evolutionary emergence of miRNA operons.

Identification and characterization of miRNAs in an endoparasitoid wasp, Pteromalus puparum

MicroRNAs (miRNAs) are a form of endogenous small noncoding RNAs that regulate protein-coding gene expression at the posttranscriptional level. So far, knowledge of miRNAs in parasitoids remains rudimentary. We investigated miRNAs in Pteromalus puparum, a pupal endoparasitoid wasp with genome and transcriptome sequences completed. In this study, we constructed eight small RNA libraries from selected developmental stages and genders: male embryos, male larvae, male pupae, male adults, mixed-sex embryos, mixed-sex larvae, mixed-sex pupae, and female adults. We identified 254 mature miRNAs with 5p/3p arm features originated from 75 known and 119 novel miRNA genes in P. puparum, 88 of which reside in 26 clusters. The miRNAs in more than half of the clusters exhibit a consistent expression pattern, indicating they were co-transcribed from a long transcript. Comparing miRNA expression in the eight libraries, we found that 84 mature miRNAs were differentially expressed between embryos and larvae, 20 between larvae and pupae, and 26 between pupae and adults. We found some miRNAs were differentially expressed between sexes in embryos (10), larvae (29), pupae (8), and adults (14). Target predictions resulted in 211,571 miRNA-mRNA interactions for 254 different mature miRNAs. These miRNAs may be involved in sexual and developmental regulation of gene expression.

MicroRNAs in Daphnia magna identified and characterized by deep sequencing, genome mapping and manual curation

MicroRNAs (miRNAs) are small non-coding RNAs that function in RNA silencing and post-transcriptional regulation of gene expression in most organisms. The water flea, Daphnia magna is a key model to study phenotypic, physiological and genomic responses to environmental cues and miRNAs can potentially mediate these responses. By using deep sequencing, genome mapping and manual curations, we have characterised the miRNAome of D. magna. We identified 66 conserved miRNAs and 13 novel miRNAs; all of these were found in the three studied life stages of D. magna (juveniles, subadults, adults), but with variation in expression levels between stages. Forty-one of the miRNAs were clustered into 13 genome clusters also present in the D. pulex genome. Most miRNAs contained sequence variants (isomiRs). The highest expressed isomiRs were 3′ template variants with one nucleotide deletion or 3′ non-template variants with addition of A or U at the 3′ end. We also identified offset RNAs (moRs) and loop RNAs (loRs). Our work extends the base for further work on all species (miRNA, isomiRs, moRNAs, loRNAs) of the miRNAome of Daphnia as biomarkers in response to chemical substances and environment cues, and underline age dependency.

MicroRNA regulation of CTP synthase and cytoophidium in Drosophila melanogaster

CTPsyn is a crucial metabolic enzyme which synthesizes CTP nucleotides. It has the extraordinary ability to compartmentalize into filaments termed cytoophidia. Though the structure is evolutionarily conserved across kingdoms, the mechanisms behind their formation remain unknown. MicroRNAs (miRNAs) are short single-stranded RNA capable of directing mRNA silencing and degradation. D. melanogaster has a high total gene count to miRNA gene number ratio, alluding to the possibility that CTPsyn too may come under their regulation. A thorough miRNA overexpression involving 123 miRNAs was conducted, followed by CTPsyn-specific staining upon cytoophidia-rich egg chambers. This revealed a small group of candidates which confer either a lengthening or truncating effect on cytoophidia, suggesting they may play a role in regulating CTPsyn. MiR-975 and miR-1014 are both cytoophidia-elongating, whereas miR-190 and miR-932 are cytoophidia-shortening. Though target prediction shows that miR-975 and miR-932 do indeed have binding sites on CTPsyn mRNA, in vitro assays instead revealed a low probability of this being true, instead indicating that the effects asserted by overexpressed miRNAs indirectly reach CTPsyn and its cytoophidia through the actions of middling elements. In silico target prediction and qPCR quantification indicated that, at least for miR-932 and miR-1014, these undetermined elements may be players in fat metabolism. This is the first study to thoroughly investigate miRNAs in connection to CTPsyn expression and activity in any species. The findings presented could serve as a basis for further queries into not only the fundamental aspects of the enzyme's regulation, but may uncover new facets of closely related pathways as well.

Identification of microRNA clusters cooperatively acting on epithelial to mesenchymal transition in triple negative breast cancer

MicroRNAs play important roles in many biological processes. Their aberrant expression can have oncogenic or tumor suppressor function directly participating to carcinogenesis, malignant transformation, invasiveness and metastasis. Indeed, miRNA profiles can distinguish not only between normal and cancerous tissue but they can also successfully classify different subtypes of a particular cancer. Here, we focus on a particular class of transcripts encoding polycistronic miRNA genes that yields multiple miRNA components. We describe 'clustered MiRNA Master Regulator Analysis (ClustMMRA)', a fully redesigned release of the MMRA computational pipeline (MiRNA Master Regulator Analysis), developed to search for clustered miRNAs potentially driving cancer molecular subtyping. Genomically clustered miRNAs are frequently co-expressed to target different components of pro-tumorigenic signaling pathways. By applying ClustMMRA to breast cancer patient data, we identified key miRNA clusters driving the phenotype of different tumor subgroups. The pipeline was applied to two independent breast cancer datasets, providing statistically concordant results between the two analyses. We validated in cell lines the miR-199/miR-214 as a novel cluster of miRNAs promoting the triple negative breast cancer (TNBC) phenotype through its control of proliferation and EMT.

Integrated network analysis identifies hsa-miR-4756-3p as a regulator of FOXM1 in Triple Negative Breast Cancer

Both aberrantly expressed mRNAs and micro(mi)RNAs play important roles in cancer cell function, which makes integration analysis difficult. In this study, we first applied master regulator analysisalgorithm and confirmed hsa-miR-4756-3p as a candidate miRNA in triple negative breast cancer (TNBC) patients; hsa-miR-4756-3p could regulate TNBC cell line apoptosis, proliferation, migration, and cell cycle as well as suppress TGF-β1 signalling andtumour growth. In TNBC, forkhead box protein M1 (FOXM1)was found to be an hsa-miR-4756-3p target gene, and FOXM1 knockout completely inhibited hsa-miR-4756-3p-induced cell migration and metastasis, TGF-β1 signalling, and epithelial mesenchymal signal activation, which indicated that hsa-miR-4756-3p functions via the FOXM1-TGFβ1-EMT axis.

Genomic non-redundancy of the mir-183/96/182 cluster and its requirement for hair cell maintenance

microRNAs are important regulators of gene expression. In the retina, the mir-183/96/182 cluster is of particular interest due to its robust expression and studies in which loss of the cluster caused photoreceptor degeneration. However, it is unclear which of the three miRNAs in the cluster are ultimately required in photoreceptors, whether each may have independent, contributory roles, or whether a single miRNA from the cluster compensates for the loss of another. These are important questions that will not only help us to understand the role of these particular miRNAs in the retina, but will deepen our understanding of how clustered microRNAs evolve and operate. To that end, we have developed a complete panel of single, double, and triple mir-183/96/182 mutant zebrafish. While the retinas of all mutant animals were normal, the triple mutants exhibited acute hair cell degeneration which corresponded with impaired swimming and death at a young age. By measuring the penetrance of this phenotype in each mutant line, we determine which of the three miRNAs in the cluster are necessary and/or sufficient to ensure normal hair cell development and function.

The Origin and Evolution of Maternal Genes

Proteins and RNA molecules are deposited into the developing egg by the mother. These gene products will drive the first stages of development and are coded by maternal genes. Maternal genes are essential, yet, despite their importance, their evolutionary dynamics is largely unknown. Here I review the current knowledge of maternal gene evolution. The evolutionary origin of maternal genes tends to be more recent than that of zygotic genes. Some studies support the theoretical prediction that maternal genes evolve faster than zygotic genes. However, most studies were done on a limited set of species and genes. I also discuss the way forward to understand the evolution of maternal genes by combining high-throughput genomics and theoretical evolutionary approaches.

Small RNA-Seq Analysis Reveals miRNA Expression Dynamics Across Tissues in the Malaria Vector, Anopheles gambiae

Malaria continues to be a major global health problem, where disease transmission is deeply linked to the repeated blood feeding nature of the anautogenous mosquito. Given the tight link between blood feeding and disease transmission, understanding basic biology behind mosquito physiology is a requirement for developing effective vector-borne disease control strategies. In the mosquito, numerous loss of function studies with notable phenotypes demonstrate microRNAs (miRNAs) play significant roles in mosquito physiology. While the field appreciates the importance of a handful of miRNAs, we still need global mosquito tissue miRNA transcriptome studies. To address this need, our goal was to determine the miRNA transcriptome for multiple tissues of the pre-vitellogenic mosquito. To this end, by using small RNA-Seq analysis, we determined miRNA transcriptomes in tissues critical for mosquito reproduction and immunity including (i) fat body-abdominal wall enriched tissues, (ii) midguts, (iii) ovaries, and (iv) remaining tissues comprised of the head and thorax. We found numerous examples of miRNAs exhibiting pan-tissue high- or low- expression, tissue exclusion, and tissue enrichment. We also updated and consolidated the miRNA catalog and provided a detailed genome architecture map for the malaria vector, Anopheles gambiae. This study aims to build a foundation for future research on how miRNAs and potentially other small RNAs regulate mosquito physiology as it relates to vector-borne disease transmission.

Identification of the X-linked germ cell specific miRNAs (XmiRs) and their functions

MicroRNAs (miRNAs) play a critical role in multiple aspects of biology. Dicer, an RNase III endonuclease, is essential for the biogenesis of miRNAs, and the germ cell-specific Dicer1 knockout mouse shows severe defects in gametogenesis. How miRNAs regulate germ cell development is still not fully understood. In this study, we identified germ cell-specific miRNAs (miR-741-3p, miR-871-3p, miR-880-3p) by analyzing published RNA-seq data of mouse. These miRNA genes are contiguously located on the X chromosome near other miRNA genes. We named them X chromosome-linked miRNAs (XmiRs). To elucidate the functions of XmiRs, we generated knockout mice of these miRNA genes using the CRISPR/Cas9-mediated genome editing system. Although no histological abnormalities were observed in testes of F0 mice in which each miRNA gene was disrupted, a deletion covering miR-871 and miR-880 or covering all XmiRs (ΔXmiRs) resulted in arrested spermatogenesis in meiosis in a few seminiferous tubules, indicating their redundant functions in spermatogenesis. Among candidate targets of XmiRs, we found increased expression of a gene encoding a WNT receptor, FZD4, in ΔXmiRs testis compared with that in wildtype testis. miR-871-3p and miR-880-3p repressed the expression of Fzd4 via the 3′-untranslated region of its mRNA. In addition, downstream genes of the WNT/β-catenin pathway were upregulated in ΔXmiRs testis. We also found that miR-871, miR-880, and Fzd4 were expressed in spermatogonia, spermatocytes and spermatids, and overexpression of miR-871 and miR-880 in germ stem cells in culture repressed their increase in number and Fzd4 expression. Previous studies indicated that the WNT/β-catenin pathway enhances and represses proliferation and differentiation of spermatogonia, respectively, and our results consistently showed that stable β-catenin enhanced GSC number. In addition, stable β-catenin partially rescued reduced GSC number by overexpression of miR-871 and miR-880. The results together suggest that miR-871 and miR-880 cooperatively regulate the WNT/β-catenin pathway during testicular germ cell development.

Deciphering conserved identical sequences of mature miRNAs among six members of great apes

MicroRNAs (miRNAs) are a group of small RNA molecules which act as negative regulators of gene expression by controlling post-transcriptional regulation through binding to their corresponding mRNAs. Due to their small size, their nucleotide compositions are expected to be similar, but until now, the extent of similarity has not been reported in humans and their six phylogenetically closely related members of hominids. The present study allows direct comparison among six members of hominid species (Homosapiens, Gorillagorilla, Panpaniscus, Pongopygmaeus, Pantroglodytes and Symphalangussyndactylus) in terms of their miRNA repertoire, their evolutionary distance to human, as well as, the categorization of identical species-specific miRNAs. For this purpose, a total of 2694, 370, 157, 673, 590 and 10 mature miRNA sequences of Homosapiens, Gorillagorilla, Panpaniscus, Pongopygmaeus, Pantroglodytes and Symphalangussyndactylus respectively were retrieved from miRbase 22. A total of 12, 4, 4 and 3 conserved clusters with identical miRNA sequences that belong to the same gene families were found in Homosapiens, Gorillagorilla, Pongopygmaeus, Pantroglodytes respectively by neighbor-joining method using MEGA7 software. Interestingly, cross-species comparison has also shown a set of conserved identical miRNA sequences. Homologs of human mature miRNAs with 100% sequence identity are expected to have similar functions in the studied primates. Further in-vitro study is required to investigate common targets for identical miRNAs in the studied primates.

MiR-30e Attenuates Isoproterenol-induced Cardiac Fibrosis Through Suppressing Snai1/TGF-β Signaling

Background:

MicroRNAs are a class of small RNA molecules that inhibit protein expression through either degradation of messenger RNA or interference with protein translation. Our previous work suggested an involvement of miR-30e in myocardial fibrosis; however, the exact role of miR-30e in the pathogenesis of cardiac fibrosis and the underlying mechanisms are not known.

Methods:

Male Sprague Dawley rats were treated with isoproterenol (ISO) to induce cardiac remodeling and fibrosis and treated with either miR-30e agomir (AG) or antagomir and respective controls. The expression of miR-30e was evaluated by reverse transcription and quantitative polymerase chain reaction. Myocardial fibrosis was assessed by Masson's trichrome staining, and the level of oxidative stress and the expression of Snai1 and transforming growth factor-beta (TGF-β) were detected using Western blots.

Results:

A significant downregulation of miR-30e was found in the hearts of ISO-treated rats with cardiac fibrosis compared with nontreated controls. In vivo administration of miR-30e AG increased the survival of ISO-treated rats compared with AG-negative control administration, which was associated with reduced oxidative stress. We further identified Snai1 as a novel miR-30e target. Snai1 expression was significantly increased in hearts from ISO-treated rats, which coincided with decreased miR-30e expression and increased TGF-β expression. An miR-30e putative target sequence was identified in the 3′-untranslated region (UTR) Snai1. In a reporter assay, miR-30e greatly suppressed the activity of wild-type 3′-UTR–fused luciferase reporter, but showed no significant effect with the mutated 3′-UTR–fused reporter.

Conclusion:

MiR-30e attenuated ISO-induced cardiac dysfunction and cardiac fibrosis in a rat cardiac remodeling model. Mechanistically, miR-30e suppressed Snai1/TGF-β pathway which was involved in ISO-induced cardiac remodeling.

MicroRNA duplication accelerates the recruitment of new targets during vertebrate evolution

The repertoire of miRNAs has considerably expanded during metazoan evolution, and duplication is an important mechanism for generating new functional miRNAs. However, relatively little is known about the functional divergence between paralogous miRNAs and the possible coevolution between duplicated miRNAs and the genomic contexts. By systematically examining small RNA expression profiles across various human tissues and interrogating the publicly available miRNA:mRNA pairing chimeras, we found that changes in expression patterns and targeting preferences are widespread for duplicated miRNAs in vertebrates. Both the empirical interactions and target predictions suggest that evolutionarily conserved homo-seed duplicated miRNAs pair with significantly higher numbers of target sites compared to the single-copy miRNAs. Our birth-and-death evolutionary analysis revealed that the new target sites of miRNAs experienced frequent gains and losses during function development. Our results suggest that a newly emerged target site has a higher probability to be functional and maintained by natural selection if it is paired to a seed shared by multiple paralogous miRNAs rather than being paired to a single-copy miRNA. We experimentally verified the divergence in target repression between two paralogous miRNAs by transfecting let-7a and let-7b mimics into kidney-derived cell lines of four mammalian species and measuring the resulting transcriptome alterations by extensive high-throughput sequencing. Our results also suggest that the gains and losses of let-7 target sites might be associated with the evolution of repressiveness of let-7 across mammalian species.

microRNA regulation in an ancient obligate endosymbiosis

Although many insects are associated with obligate bacterial endosymbionts, the mechanisms by which these host/endosymbiont associations are regulated remain mysterious. While microRNAs (miRNAs) have been recently identified as regulators of host/microbe interactions, including host/pathogen and host/facultative endosymbiont interactions, the role miRNAs may play in mediating host/obligate endosymbiont interactions is virtually unknown. Here, we identified conserved miRNAs that potentially mediate symbiotic interactions between aphids and their obligate endosymbiont, Buchnera aphidicola. Using small RNA sequence data from Myzus persicae and Acyrthosiphon pisum, we annotated 93 M. persicae and 89 A. pisum miRNAs, among which 69 were shared. We found 14 miRNAs that were either highly expressed in aphid bacteriome, the Buchnera-housing tissue, or differentially expressed in bacteriome vs. gut, a non-Buchnera-housing tissue. Strikingly, 10 of these 14 miRNAs have been implicated previously in other host/microbe interaction studies. Investigating the interaction networks of these miRNAs using a custom computational pipeline, we identified 103 miRNA::mRNA interactions shared between M. persicae and A. pisum. Functional annotation of the shared mRNA targets revealed only two over-represented cluster of orthologous group categories: amino acid transport and metabolism, and signal transduction mechanisms. Our work supports a role for miRNAs in mediating host/symbiont interactions between aphids and their obligate endosymbiont Buchnera. In addition, our results highlight the probable importance of signal transduction mechanisms to host/endosymbiont coevolution.

Nonprotein-Coding RNAs as Regulators of Development in Tunicates

Tunicates, or urochordates, are a group of small marine organisms that are found widely throughout the seas of the world. As most plausible sister group of the vertebrates, they are of utmost importance for a comprehensive understanding of chordate evolution; hence, they have served as model organisms for many aspects of the developmental biology. Current genomic analysis of tunicates indicates that their genomes evolved with a fast rate not only at the level of nucleotide substitutions but also in terms of genomic organization. The latter involves genome reduction, rearrangements, as well as the loss of some important coding and noncoding RNA (ncRNAs) elements and even entire genomic regions that are otherwise well conserved. These observations are largely based on evidence from comparative genomics resulting from the analysis of well-studied gene families such as the Hox genes and their noncoding elements. In this chapter, the focus lies on the ncRNA complement of tunicates, with a particular emphasis on microRNAs, which have already been studied extensively for other animal clades. MicroRNAs are known as important regulators of key genes in animal development, and they are intimately related to the increase morphological complexity in higher metazoans. Here we review the discovery, evolution, and genome organization of the miRNA repertoire, which has been drastically reduced and restructured in tunicates compared to the chordate ancestor. Known functions of microRNAs as regulators of development in tunicates are a central topic. For instance, we consider the role of miRNAs as regulators of the muscle development and their importance in the regulation of the differential expression during the oral siphon regeneration. Beyond microRNAs, we touch upon the functions of some other ncRNAs such as yellow crescent RNA, moRNAs, RMST lncRNAs, or spliced-leader (SL) RNAs, which have diverse functions associated with the embryonic development, neurogenesis, and mediation of mRNA stability in general.

A miRNA catalogue and ncRNA annotation of the short-living fish Nothobranchius furzeri

BACKGROUND

The short-lived fish Nothobranchius furzeri is the shortest-lived vertebrate that can be cultured in captivity and was recently established as a model organism for aging research. Small non-coding RNAs, especially miRNAs, are implicated in age dependent control of gene expression.

RESULTS

Here, we present a comprehensive catalogue of miRNAs and several other non-coding RNA classes (ncRNAs) for Nothobranchius furzeri. Analyzing multiple small RNA-Seq libraries, we show most of these identified miRNAs are expressed in at least one of seven Nothobranchius species. Additionally, duplication and clustering of N. furzeri miRNAs was analyzed and compared to the four fish species Danio rerio, Oryzias latipes, Gasterosteus aculeatus and Takifugu rubripes. A peculiar characteristic of N. furzeri, as compared to other teleosts, was a duplication of the miR-29 cluster.

CONCLUSION

The completeness of the catalogue we provide is comparable to that of the zebrafish. This catalogue represents a basis to investigate the role of miRNAs in aging and development in this species.

MicroRNAs That Contribute to Coordinating the Immune Response in Drosophila melanogaster

Small noncoding RNAs called microRNAs (miRNAs) have emerged as post-transcriptional regulators of gene expression related to host defenses. Here, we have used Drosophila melanogaster to explore the contribution of individual or clusters of miRNAs in countering systemic Candida albicans infection. From a total of 72 tested, we identify 6 miRNA allelic mutant backgrounds that modulate the survival response to infection and the ability to control pathogen number. These mutants also exhibit dysregulation of the Toll pathway target transcripts Drosomycin (Drs) and Immune-Induced Molecule 1 (IM1). These are characteristics of defects in Toll signaling, and consistent with this, we demonstrate dependency for one of the miRNA mutants on the NF-κΒ homolog Dif. We also quantify changes in the miRNA expression profile over time in response to three pathogen types, and identify 13 mature miRNA forms affected by pathogens that stimulate Toll signaling. To complement this, we provide a genome-wide map of potential NF-κB sites in proximity to miRNA genes. Finally, we demonstrate that systemic C. albicans infection contributes to a reduction in the total amount of branch-chained amino acids, which is miRNA-regulated. Overall, our data reveal a new layer of miRNA complexity regulating the fly response to systemic fungal infection.

Chronic low-dose exposure to a mixture of environmental endocrine disruptors induces microRNAs/isomiRs deregulation in mouse concomitant with intratesticular estradiol reduction

Humans are environmentally exposed not only to single endocrine-disrupting chemicals (EDCs) but to mixtures that affect their reproductive health. In reproductive tissues, microRNAs (miRNAs) are emerging as key targets of EDCs. Here, we analysed changes in the testis "miRNome" (and their biogenesis mechanism) in chronically exposed adult mice to a cocktail of five EDCs containing 0.3 mg/kg-body weight (BW)/day of each phthalate (DEHP, DBP, BBP) and 0.05 mg/kg-BW/day of each alkylphenol (NP, OP), from conception to adulthood. The testis "miRNome" was characterised using next-generation sequencing (NGS). Expression levels of genes involved in miRNA biogenesis were measured by RT-qPCR, as well as several physiological and cytological parameters. We found two up-regulated, and eight down-regulated miRNAs and thirty-six differentially expressed isomiRs along with an over-expression of Drosha, Adar and Zcchc11. A significant decrease of intratesticular estradiol but not testosterone was detected. Functional analysis showed altered spermatogenesis, germ cell apoptosis and negative correlation of miR-18a-5p with Nr1h2 involved in the deregulation of the steroidogenesis pathway. Here, we present the first association between miRNA/isomiRs deregulation, their mechanisms of biogenesis and histopathological and hormonal alterations in testes of adult mice exposed to a mixture of low-dose EDCs, which can play a role in male infertility.

microRNAs in the Same Clusters Evolve to Coordinately Regulate Functionally Related Genes

MicroRNAs (miRNAs) are endogenously expressed small noncoding RNAs. The genomic locations of animal miRNAs are significantly clustered in discrete loci. We found duplication and de novo formation were important mechanisms to create miRNA clusters and the clustered miRNAs tend to be evolutionarily conserved. We proposed a "functional co-adaptation" model to explain how clustering helps newly emerged miRNAs survive and develop functions. We presented evidence that abundance of miRNAs in the same clusters were highly correlated and those miRNAs exerted cooperative repressive effects on target genes in human tissues. By transfecting miRNAs into human and fly cells and extensively profiling the transcriptome alteration with deep-sequencing, we further demonstrated the functional co-adaptation between new and old miRNAs in the miR-17-92 cluster. Our population genomic analysis suggest that positive Darwinian selection might be the driving force underlying the formation and evolution of miRNA clustering. Our model provided novel insights into mechanisms and evolutionary significance of miRNA clustering.

Pervasive microRNA Duplication in Chelicerates: Insights from the Embryonic microRNA Repertoire of the Spider Parasteatoda tepidariorum

MicroRNAs are small (∼22 nt) noncoding RNAs that repress translation and therefore regulate the production of proteins from specific target mRNAs. microRNAs have been found to function in diverse aspects of gene regulation within animal development and many other processes. Among invertebrates, both conserved and novel, lineage specific, microRNAs have been extensively studied predominantly in holometabolous insects such as Drosophila melanogaster However little is known about microRNA repertoires in other arthropod lineages such as the chelicerates. To understand the evolution of microRNAs in this poorly sampled subphylum, we characterized the microRNA repertoire expressed during embryogenesis of the common house spider Parasteatoda tepidariorum We identified a total of 148 microRNAs in P. tepidariorum representing 66 families. Approximately half of these microRNA families are conserved in other metazoans, while the remainder are specific to this spider. Of the 35 conserved microRNAs families 15 had at least two copies in the P. tepidariorum genome. A BLAST-based approach revealed a similar pattern of duplication in other spiders and a scorpion, but not among other chelicerates and arthropods, with the exception of a horseshoe crab. Among the duplicated microRNAs we found examples of lineage-specific tandem duplications, and the duplication of entire microRNA clusters in three spiders, a scorpion, and in a horseshoe crab. Furthermore, we found that paralogs of many P. tepidariorum microRNA families exhibit arm switching, which suggests that duplication was often followed by sub- or neofunctionalization. Our work shows that understanding the evolution of microRNAs in the chelicerates has great potential to provide insights into the process of microRNA duplication and divergence and the evolution of animal development.

Differences in molecular evolutionary rates among microRNAs in the human and chimpanzee genomes

Background

The rise of the primate lineage is accompanied by an outstanding emergence of microRNAs, small non-coding RNAs with a prominent role in gene regulation. In spite of their biological importance little is known about the way in which natural selection has influenced microRNAs in the human lineage. To study the recent evolutionary history of human microRNAs and to analyze the signatures of natural selection in genomic regions harbouring microRNAs we have investigated the nucleotide substitution rates of 1,872 human microRNAs in the human and chimpanzee lineages.

Results

We produced a depurated set of microRNA alignments of human, chimpanzee and orang-utan orthologs combining BLAT and liftOver and selected 1,214 microRNA precursors presenting optimal secondary structures. We classified microRNAs in categories depending on their genomic organization, duplication status and conservation along evolution. We compared substitution rates of the aligned microRNAs between human and chimpanzee using Tajima’s Relative Rate Test taking orang-utan as out-group and found several microRNAs with particularly high substitution rates in either the human or chimpanzee branches. We fitted different models of natural selection on these orthologous microRNA alignments and compared them using a likelihood ratio test that uses ancestral repeats and microRNA flanking regions as neutral sequences. We found that although a large fraction of human microRNAs is highly conserved among the three species studied, significant differences in rates of molecular evolution exist among microRNA categories. Particularly, primate-specific microRNAs, which are enriched in isolated and single copy microRNAs, more than doubled substitution rates of those belonging to older, non primate-specific microRNA families.

Conclusions

Our results corroborate the remarkable conservation of microRNAs, a proxy of their functional relevance, and indicate that a subset of human microRNAs undergo nucleotide substitutions at higher rates, which may be suggestive of the action of positive selection.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-016-2863-3) contains supplementary material, which is available to authorized users.

Microsynteny and phylogenetic analysis of tandemly organised miRNA families across five members of Brassicaceae reveals complex retention and loss history

Plant genomes are characterized by the presence of large miRNA gene families which are few in number. The expansion of miRNA families is thought to be driven by gene and genome duplication. Some members of these miRNA gene families are tandemly arranged and their analysis is of interest because such organisation may indicate origin through tandem duplication and also to investigate whether some such tandem clusters have similar expression patterns, and whether these are regulated through a common set of cis-regulatory elements (eg. promoters and enhancers). As a first step, we undertake a comprehensive study using micro-synteny analyses of tandemly organised miRNA families across the Brassicaceae spanning an evolutionary time scale of ca. 45 million years, among Arabidopsis, Capsella, Brassica and Thellungiella species, to address the following questions: Are most miRNA gene families present as tandem clusters? To what extent are these tandem patterns retained? To what extent can family sizes be ascribed to genome duplication? Our analysis of thirteen tandemly organised miRNA families revealed that synteny is largely conserved among Arabidopsis thaliana, A. lyrata and Capsella rubella, which form a clade spanning approximately between 6.2-9.8 my (Acarkan et al., 2000) [1]. On the other hand, comparison of sequences from these species with Brassica rapa, B. oleracea and Thellungiella halophila, which form a separate clade spanning 31 my (Franzke et al., 2011)[2] reveals many differences. The latter clade reveals several paralogous duplications that probably resulted from whole genome duplication, as well as disrupted synteny. Phylogenetic analyses of precursor sequences generally support the history inferred from synteny analysis. Synteny and phylogenetic analysis of six members of the tandemly organised miR169 family suggest that the Brassicaceae ancestral state consisted of a "dimer as a unit" which may have undergone direct local duplication to retain the transcriptional orientation followed by lineage specific changes. MiR169, to the best of our knowledge, is one of the largest tandemly organised miRNA gene family across plant kingdom and further analysis should reveal the generality of this pattern of evolution. The conserved organisation of miR395A-B-C and miR395 D-E-F as two clusters on same chromosome/scaffold across A. thaliana, B. rapa and salsuginea demonstrates retention of the large chromosomal segment across the two lineages. MiRNA family miR845 was detected only in Arabidopsis species and Thellungiella indicating a complex loss and retention history. MiR447A-B family was only found in A. thaliana indicating that it is a species-specific gene family of recent origin.

Identification and characterization of microRNAs in the white-backed planthopper, Sogatella furcifera

MicroRNAs (miRNAs) are a novel class of small, non-coding endogenous RNAs that play critical regulatory roles in many metabolic activities in eukaryotes. Reports of the identification of miRNAs in Sogatella furcifera (white-backed planthopper), the insect that acts as the only confirmed vector of the southern rice black-streaked dwarf virus (SRBSDV), are limited. In this study, a total of 382 miRNAs were identified in S. furcifera, including 106 conserved and 276 novel miRNAs, using high-throughput sequencing based on two small RNA libraries from viruliferous and non-viruliferous S. furcifera, and these miRNAs belonged to 52 conserved miRNA families and 58 S. furcifera-specific families, respectively. Comparison with miRNAs from 26 insect species and five other species in miRBase showed that more than half of the conserved miRNA families are highly conserved in Hexapoda, while other miRNAs are only conserved in non-dipterans. Furthermore, 4 117 target genes predicted for the 382 identified miRNAs could be categorized into 45 functional groups annotated by Gene Ontology. Compared with non-viruliferous cells, eight up-regulated miRNAs and four down-regulated miRNAs were identified in cells inoculated with SRBSDV, among which miR-14 and miR-n98a may be involved in the immune response to SRBSDV infection. Analyses of the identified miRNAs will provide insights into the roles of these miRNAs in the regulation and expression of genes involved in the metabolism, development and viral infection of S. furcifera.

Specific down-regulation of spermatogenesis genes targeted by 22G RNAs in hybrid sterile males associated with an X-Chromosome introgression

Hybrid incompatibility (HI) prevents gene flow between species, thus lying at the heart of speciation genetics. One of the most common HIs is male sterility. Two superficially contradictory observations exist for hybrid male sterility. First, an introgression on the X Chromosome is more likely to produce male sterility than on autosome (so-called large-X theory); second, spermatogenesis genes are enriched on the autosomes but depleted on the X Chromosome (demasculinization of X Chromosome). Analysis of gene expression in Drosophila hybrids suggests a genetic interaction between the X Chromosome and autosomes that is essential for male fertility. However, the prevalence of such an interaction and its underlying mechanism remain largely unknown. Here we examine the interaction in nematode species by contrasting the expression of both coding genes and transposable elements (TEs) between hybrid sterile males and its parental nematode males. We use two lines of hybrid sterile males, each carrying an independent introgression fragment from Caenorhabditis briggsae X Chromosome in an otherwise Caenorhabditis nigoni background, which demonstrate similar defects in spermatogenesis. We observe a similar pattern of down-regulated genes that are specific for spermatogenesis between the two hybrids. Importantly, the down-regulated genes caused by the X Chromosome introgressions show a significant enrichment on the autosomes, supporting an epistatic interaction between the X Chromosome and autosomes. We investigate the underlying mechanism of the interaction by measuring small RNAs and find that a subset of 22G RNAs specifically targeting the down-regulated spermatogenesis genes is significantly up-regulated in hybrids, suggesting that perturbation of small RNA-mediated regulation may contribute to the X-autosome interaction.

Promoter-based identification of novel non-coding RNAs reveals the presence of dicistronic snoRNA-miRNA genes in Arabidopsis thaliana

Background

In the past few decades, non-coding RNAs (ncRNAs) have emerged as important regulators of gene expression in eukaryotes. Most studies of ncRNAs in plants have focused on the identification of silencing microRNAs (miRNAs) and small interfering RNAs (siRNAs). Another important family of ncRNAs that has been well characterized in plants is the small nucleolar RNAs (snoRNAs) and the related small Cajal body-specific RNAs (scaRNAs). Both target chemical modifications of ribosomal RNAs (rRNAs) and small nuclear RNAs (snRNAs). In plants, the snoRNA genes are organized in clusters, transcribed by RNA Pol II from a common promoter and subsequently processed into mature molecules. The promoter regions of snoRNA polycistronic genes in plants are highly enriched in two conserved cis-regulatory elements (CREs), Telo-box and Site II, which coordinate the expression of snoRNAs and ribosomal protein coding genes throughout the cell cycle.

Results

In order to identify novel ncRNA genes, we have used the snoRNA Telo-box/Site II motifs combination as a functional promoter indicator to screen the Arabidopsis genome. The predictions generated by this process were tested by detailed exploration of available RNA-Seq and expression data sets and experimental validation. As a result, we have identified several snoRNAs, scaRNAs and 'orphan' snoRNAs. We also show evidence for 16 novel ncRNAs that lack similarity to any reported RNA family. Finally, we have identified two dicistronic genes encoding precursors that are processed to mature snoRNA and miRNA molecules. We discuss the evolutionary consequences of this result in the context of a tight link between snoRNAs and miRNAs in eukaryotes.

Conclusions

We present an alternative computational approach for non-coding RNA detection. Instead of depending on sequence or structure similarity in the whole genome screenings, we have explored the properties of promoter regions of well-characterized ncRNAs. Interestingly, besides expected ncRNAs predictions we were also able to recover single precursor arrangement for snoRNA-miRNA. Accompanied by analyses performed on rice sequences, we conclude that such arrangement might have interesting functional and evolutionary consequences and discuss this result in the context of a tight link between snoRNAs and miRNAs in eukaryotes.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-015-2221-x) contains supplementary material, which is available to authorized users.

Identification of Novel and Conserved microRNAs in Homalodisca vitripennis, the Glassy-Winged Sharpshooter by Expression Profiling

The glassy-winged sharpshooter (GWSS) Homalodisca vitripennis (Hemiptera: Cicadellidae), is a xylem-feeding leafhopper and an important vector of the bacterium Xylella fastidiosa; the causal agent of Pierce’s disease of grapevines. MicroRNAs are a class of small RNAs that play an important role in the functional development of various organisms including insects. In H. vitripennis, we identified microRNAs using high-throughput deep sequencing of adults followed by computational and manual annotation. A total of 14 novel microRNAs that are not found in the miRBase were identified from adult H. vitripennis. Conserved microRNAs were also found in our datasets. By comparison to our previously determined transcriptome sequence of H. vitripennis, we identified the potential targets of the microRNAs in the transcriptome. This microRNA profile information not only provides a more nuanced understanding of the biological and physiological mechanisms that govern gene expression in H. vitripennis, but may also lead to the identification of novel mechanisms for biorationally designed management strategies through the use of microRNAs.

Protein and small non-coding RNA-enriched extracellular vesicles are released by the pathogenic blood fluke Schistosoma mansoni

Background

Penetration of skin, migration through tissues and establishment of long-lived intravascular partners require Schistosoma parasites to successfully manipulate definitive host defences. While previous studies of larval schistosomula have postulated a function for excreted/secreted (E/S) products in initiating these host-modulatory events, the role of extracellular vesicles (EVs) has yet to be considered. Here, using preparatory ultracentrifugation as well as methodologies to globally analyse both proteins and small non-coding RNAs (sncRNAs), we conducted the first characterization of Schistosoma mansoni schistosomula EVs and their potential host-regulatory cargos.

Results

Transmission electron microscopy analysis of EVs isolated from schistosomula in vitro cultures revealed the presence of numerous, 30–100 nm sized exosome-like vesicles. Proteomic analysis of these vesicles revealed a core set of 109 proteins, including homologs to those previously found enriched in other eukaryotic EVs, as well as hypothetical proteins of high abundance and currently unknown function. Characterization of E/S sncRNAs found within and outside of schistosomula EVs additionally identified the presence of potential gene-regulatory miRNAs (35 known and 170 potentially novel miRNAs) and tRNA-derived small RNAs (tsRNAs; nineteen 5′ tsRNAs and fourteen 3′ tsRNAs).

Conclusions

The identification of S. mansoni EVs and the combinatorial protein/sncRNA characterization of their cargo signifies that an important new participant in the complex biology underpinning schistosome/host interactions has now been discovered. Further work defining the role of these schistosomula EVs and the function/stability of intra- and extra-vesicular sncRNA components presents tremendous opportunities for developing novel schistosomiasis diagnostics or interventions.