MicroRNA evolution, expression, and function during short germband development in Tribolium castaneum

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.

Identification and Characterization of Development-Related microRNAs in the Red Flour Beetle, Tribolium castaneum

MicroRNAs (miRNAs) play important roles in insect growth and development, but they were poorly studied in insects. In this study, a total of 883 miRNAs were detected from the early embryo (EE), late larva (LL), early pupa (EP), late pupa (LP), and early adult (EA) of Tribolium castaneum by microarray assay. Further analysis identified 179 differentially expressed unique miRNAs (DEmiRNAs) during these developmental stages. Of the DEmiRNAs, 102 DEmiRNAs exhibited stage-specific expression patterns during development, including 53 specifically highly expressed miRNAs and 20 lowly expressed miRNAs in EE, 19 highly expressed miRNAs in LL, 5 weakly expressed miRNAs in EP, and 5 abundantly expressed miRNAs in EA. These miRNAs were predicted to target 747, 265, 472, 234, and 121 genes, respectively. GO enrichment analysis indicates that the targets were enriched by protein phosphorylation, calcium ion binding, sequence-specific DNA binding transcription factor activity, and cytoplasm. An RNA interference-mediated knockdown of the DEmiRNAs tca-miR-6-3p, tca-miR-9a-3p, tca-miR-9d-3p, tca-miR-11-3p, and tca-miR-13a-3p led to defects in metamorphosis and wing development of T. castaneum. This study has completed the identification and characterization of development-related miRNAs in T. castaneum, and will enable us to investigate their roles in the growth and development of insect.

What Are the Functional Roles of Piwi Proteins and piRNAs in Insects?

Research on Piwi proteins and piRNAs in insects has focused on three experimental models: oogenesis and spermatogenesis in Drosophila melanogaster, the antiviral response in Aedes mosquitoes and the molecular analysis of primary and secondary piRNA biogenesis in Bombyx mori-derived BmN4 cells. Significant unique and complementary information has been acquired and has led to a greater appreciation of the complexity of piRNA biogenesis and Piwi protein function. Studies performed in other insect species are emerging and promise to add to the current state of the art on the roles of piRNAs and Piwi proteins. Although the primary role of the piRNA pathway is genome defense against transposons, particularly in the germline, recent findings also indicate an expansion of its functions. In this review, an extensive overview is presented of the knowledge of the piRNA pathway that so far has accumulated in insects. Following a presentation of the three major models, data from other insects were also discussed. Finally, the mechanisms for the expansion of the function of the piRNA pathway from transposon control to gene regulation were considered.

Maternal body condition and season influence RNA deposition in the oocytes of alfalfa leafcutting bees (Megachile rotundata)

Maternal effects are an important source of phenotypic variance, whereby females influence offspring developmental trajectory beyond direct genetic contributions, often in response to changing environmental conditions. However, relatively little is known about the mechanisms by which maternal experience is translated into molecular signals that shape offspring development. One such signal may be maternal RNA transcripts (mRNAs and miRNAs) deposited into maturing oocytes. These regulate the earliest stages of development of all animals, but are understudied in most insects. Here we investigated the effects of female internal (body condition) and external (time of season) environmental conditions on maternal RNA in the maturing oocytes and 24-h-old eggs (24-h eggs) of alfalfa leafcutting bees. Using gene expression and WGCNA analysis, we found that females adjust the quantity of mRNAs related to protein phosphorylation, transcriptional regulation, and nuclease activity deposited into maturing oocytes in response to both poor body condition and shorter day lengths that accompany the late season. However, the magnitude of these changes was higher for time of season. Females also adjusted miRNA deposition in response to seasonal changes, but not body condition. We did not observe significant changes in maternal RNAs in response to either body condition or time of season in 24-h eggs, which were past the maternal-to-zygotic transition. Our results suggest that females adjust the RNA transcripts they provide for offspring to regulate development in response to both internal and external environmental cues. Variation in maternal RNAs may, therefore, be important for regulating offspring phenotype in response to environmental change.

The embryonic transcriptome of Parhyale hawaiensis reveals different dynamics of microRNAs and mRNAs during the maternal-zygotic transition

Parhyale hawaiensis has emerged as the crustacean model of choice due to its tractability, ease of imaging, sequenced genome, and development of CRISPR/Cas9 genome editing tools. However, transcriptomic datasets spanning embryonic development are lacking, and there is almost no annotation of non-protein-coding RNAs, including microRNAs. We have sequenced microRNAs, together with mRNAs and long non-coding RNAs, in Parhyale using paired size-selected RNA-seq libraries at seven time-points covering important transitions in embryonic development. Focussing on microRNAs, we annotate 175 loci in Parhyale, 88 of which have no known homologs. We use these data to annotate the microRNAome of 37 crustacean genomes, and suggest a core crustacean microRNA set of around 61 sequence families. We examine the dynamic expression of microRNAs and mRNAs during the maternal-zygotic transition. Our data suggest that zygotic genome activation occurs in two waves in Parhyale with microRNAs transcribed almost exclusively in the second wave. Contrary to findings in other arthropods, we do not predict a general role for microRNAs in clearing maternal transcripts. These data significantly expand the available transcriptomics resources for Parhyale, and facilitate its use as a model organism for the study of small RNAs in processes ranging from embryonic development to regeneration.

Yak miR-2285o-3p attenuates hypoxia-induced apoptosis by targeting caspase-3

miRNAs are a class of hairpin‐derived RNAs, 21–24 nucleotides in length, which are involved in a range of biological processes. The bta‐miR‐2285 family has over 40 members spanning the entire bovine genome. We previously found that bta‐miR‐2285o‐3p was highly expressed in yak heart and lung when compared with cattle, which prompted us to investigate its potential function in high‐altitude adaptation of yaks. In this study, we detected wide‐spread high expression of bta‐miR‐2285o‐3p in yak tissues. Further experiments revealed that the protein tyrosine phosphatase receptor type M (PTPRM) gene was the host gene of bta‐miR‐2285o‐3p and that two linked SNPs in bta‐mir‐2285o precursor affected the biogenesis of mature miRNA (bta‐miR‐2285o‐3p). Functional analysis in vitro indicated that bta‐miR‐2285o‐3p attenuated hypoxia‐induced apoptosis by targeting very low‐density lipoprotein receptor (VLDLR), phosphatase and tensin homolog (PTEN) and caspase‐3. Expression level analysis in vivo revealed the high negative Pearson’s correlation between bta‐miR‐2285o‐3p and caspase3 in yak, highlighting the potential important roles of bta‐miR‐2285o‐3p in yak high‐altitude adaptation. Our study provides a typical model for deciphering the function of miRNAs in environmental adaptation.

Large-Scale Annotation and Evolution Analysis of MiRNA in Insects

Insects are among the most diverse and successful groups of animals and exhibit great morphological diversity and complexity. The innovation of wings and metamorphosis are some examples of the fascinating biological evolution of insects. Most microRNAs (miRNAs) contribute to canalization by conferring robustness to gene networks and thus increase the heritability of important phenotypes. Though previous studies have demonstrated how miRNAs regulate important phenotypes, little is still known about miRNA evolution in insects. Here, we used both small RNA-seq data and homology searching methods to annotate the miRNA repertoires of 152 arthropod species, including 135 insects and 17 noninsect arthropods. We identified 16,212 miRNA genes, and classified them into highly conserved (62), insect-conserved (90), and lineage-specific (354) miRNA families. The phylogenetic relationship of miRNA binary presence/absence dynamics implies that homoplastic loss of conserved miRNA families tends to occur in far-related morphologically simplified taxa, including scale insects (Coccoidea) and twisted-wing insects (Strepsiptera), leading to inconsistent phylogenetic tree reconstruction. The common ancestor of Insecta shares 62 conserved miRNA families, of which five were rapidly gained in the early winged-insects (Pterygota). We also detected extensive miRNA losses in Paraneoptera that are correlated with morphological reduction, and miRNA gains in early Endopterygota around the time holometabolous metamorphosis appeared. This was followed by abundant miRNA gains in Hymenoptera and Lepidoptera. In summary, we provide a comprehensive data set and a detailed evolutionary analysis of miRNAs in insects. These data will be important for future studies on miRNA functions associated with insect morphological innovation and trait biodiversity.

Effects of Dinotefuran on Brain miRNA Expression Profiles in Young Adult Honey Bees (Hymenopptera: Apidae)

Honey bees are important pollinators of wild plants and crops. MicroRNAs (miRNAs) are endogenous regulators of gene expression. In this study, we initially determined that the lethal concentration 50 (LC50) of dinotefuran was 0.773 mg/l. Then, the expression profiles and differentially expressed miRNAs (DE miRNAs) in honey bee brains after 1, 5, and 10 d of treatment with the lethal concentration 10 (LC10) of dinotefuran were explored via deep small-RNA sequencing and bioinformatics. In total, 2, 23, and 27 DE miRNAs were identified after persistent exposure to the LC10 of dinotefuran for 1, 5, and 10 d, respectively. Some abundant miRNAs, such as ame-miR-375-3p, ame-miR-281-5p, ame-miR-3786-3p, ame-miR-10-5p, and ame-miR-6037-3p, were extremely significantly differentially expressed. Enrichment analysis suggested that the candidate target genes of the DE miRNAs are involved in the regulation of biological processes, cellular processes, and behaviors. These results expand our understanding of the regulatory roles of miRNAs in honey bee Apis mellifera (Hymenopptera: Apidae) responses to neonicotinoid insecticides and facilitate further studies on the functions of miRNAs in honey bees.

Distinct Regulation of the Expression of Satellite DNAs in the Beetle Tribolium castaneum

In the flour beetle, Tribolium castaneum (peri)centromeric heterochromatin is mainly composed of a major satellite DNA TCAST1 interspersed with minor satellites. With the exception of heterochromatin, clustered satellite repeats are found dispersed within euchromatin. In order to uncover a possible satellite DNA function within the beetle genome, we analysed the expression of the major TCAST1 and a minor TCAST2 satellite during the development and upon heat stress. The results reveal that TCAST1 transcription was strongly induced at specific embryonic stages and upon heat stress, while TCAST2 transcription is stable during both processes. TCAST1 transcripts are processed preferentially into piRNAs during embryogenesis and into siRNAs during later development, contrary to TCAST2 transcripts, which are processed exclusively into piRNAs. In addition, increased TCAST1 expression upon heat stress is accompanied by the enrichment of the silent histone mark H3K9me3 on the major satellite, while the H3K9me3 level at TCAST2 remains unchanged. The transcription of the two satellites is proposed to be affected by the chromatin state: heterochromatin and euchromatin, which are assumed to be the prevalent sources of TCAST1 and TCAST2 transcripts, respectively. In addition, distinct regulation of the expression might be related to diverse roles that major and minor satellite RNAs play during the development and stress response.

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.

A Day in the Life: Identification of Developmentally Regulated MicroRNAs in the Colorado Potato Beetle (Leptinotarsa decemlineata; Coleoptera: Chrysomelidae)

The Colorado potato beetle (Leptinotarsa decemlineata (Say)) is an important pest of the cultivated potato (Solanum tuberosum (L.) [Solanales: Solanaceae]). With its broad resistance toward commonly used insecticides, it is clear that more sophisticated control strategies are needed. Due to their importance in insect development, microRNAs (miRNAs) represent a potential tool to employ in insect control strategies. However, most studies conducted in this area have focused on model species with well-annotated genomes. In this study, next-generation sequencing was used to catalogue the miRNAs produced by L. decemlineata across all eight stages of its development, from eggs to adults. For most stages, the length of miRNAs peaked between 21 and 22 nt, though it was considerably longer for the egg stage (26 nt). Global profiling of miRNAs revealed three distinct developmental clusters: 1) egg stage; 2) early stage (first, second, and third instar); and 3) late stage (fourth instar, prepupae, pupae, and adult). We identified 86 conserved miRNAs and 33 bonafide novel miRNAs, including stage-specific miRNAs and those not previously identified in L. decemlineata. Most of the conserved miRNAs were found in multiple developmental stages, whereas the novel miRNAs were often stage specific with the bulk identified in the egg stage. The identified miRNAs have a myriad of putative functions, including growth, reproduction, and insecticide resistance. We discuss the putative roles of some of the most notable miRNAs in the regulation of L. decemlineata development, as well as the potential applications of this research in Colorado potato beetle management.

MicroRNA-574-5p in gastric cancer cells promotes angiogenesis by targeting protein tyrosine phosphatase non-receptor type 3 (PTPN3)

We elucidate in this study that up-regulation of miR-574-5p in gastric cancer cells under hypoxic conditions contributed to angiogenesis. We found that miR-574-5p and HIF-1α were up-regulated in gastric cancer cells cultured under 2% O₂ or in medium containing CoCl_2, and in muscle tissues of mice injected with NaNO₂, indicating up-regulation of miR-574-5p in vitro or in vivo in response to hypoxic conditions. We hypothesized that up-regulation of miR-574-5p could promote angiogenesis. Transfection of gastric cancer cells with miR-574-5p mimics or inhibitor resulted in increase or decrease in the expression of VEGFA. Viability, migration, invasion and tube formation of HUVECs cultured with conditioned medium from SGC/574 cells transfected with miR-574-5p inhibitor were reduced. Tube formation of HUVECs cultured with conditioned medium from SGC-7901 cells transfected with miR-574-5p mimics was increased. An in vivo study demonstrated that inhibition of miR-574-5p in the tumor xenografts of mice reduced the expression of CD31 one of the endothelial cell markers. We identified PTPN3 a tyrosine phosphatase as a target of miR-574-5p that bound to the 3'UTR of PTPN3 mRNA to inhibit the expression of PTPN3. Furthermore, the data in this study demonstrated that inhibition of PTPN3 in gastric cancer cells enhanced phosphorylation of p44/42 MAPKs and promoted angiogenesis. We conclude that miR-574-5p in gastric cancer cells promoted angiogenesis via enhancing phosphorylation of p44/42 MAPKs by miR-574-5p inhibition of PTPN3 expression.

Enhanced genome assembly and a new official gene set for Tribolium castaneum

Background

The red flour beetle Tribolium castaneum has emerged as an important model organism for the study of gene function in development and physiology, for ecological and evolutionary genomics, for pest control and a plethora of other topics. RNA interference (RNAi), transgenesis and genome editing are well established and the resources for genome-wide RNAi screening have become available in this model. All these techniques depend on a high quality genome assembly and precise gene models. However, the first version of the genome assembly was generated by Sanger sequencing, and with a small set of RNA sequence data limiting annotation quality.

Results

Here, we present an improved genome assembly (Tcas5.2) and an enhanced genome annotation resulting in a new official gene set (OGS3) for Tribolium castaneum, which significantly increase the quality of the genomic resources. By adding large-distance jumping library DNA sequencing to join scaffolds and fill small gaps, the gaps in the genome assembly were reduced and the N50 increased to 4753kbp. The precision of the gene models was enhanced by the use of a large body of RNA-Seq reads of different life history stages and tissue types, leading to the discovery of 1452 novel gene sequences. We also added new features such as alternative splicing, well defined UTRs and microRNA target predictions. For quality control, 399 gene models were evaluated by manual inspection. The current gene set was submitted to Genbank and accepted as a RefSeq genome by NCBI.

Conclusions

The new genome assembly (Tcas5.2) and the official gene set (OGS3) provide enhanced genomic resources for genetic work in Tribolium castaneum. The much improved information on transcription start sites supports transgenic and gene editing approaches. Further, novel types of information such as splice variants and microRNA target genes open additional possibilities for analysis.

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.

Differentially and Co-expressed Genes in Embryo, Germ-Line and Somatic Tissues of Tribolium castaneum

Transcriptomic studies of Tribolium castaneum have led to significant advances in our understanding of co-regulation and differential expression of genes in development. However, previously used microarray approaches have covered only a subset of known genes. The aim of this study was to investigate gene expression patterns of beetle embryo, germ-line and somatic tissues. We identified 12,302 expressed genes and determined differentially expressed up and down-regulated genes among all samples. For example, 1624 and 3639 genes were differentially increased in expression greater than or equal to twofold change (FDR < 0.01) in testis vs. ovary (virgin female) and ovary vs. embryo (0-5 hr), respectively. Of these, many developmental, somatic and germ-line differentially expressed genes were identified. Furthermore, many maternally deposited transcripts were identified, whose expression either decreased rapidly or persisted during embryogenesis. Genes with the largest change in expression were predominantly decreased during early embryogenesis as compared to ovary or were increased in testis compared to embryo. We also identify zygotic genes induced after fertilization. The genome wide variation in transcript regulation in maternal and zygotic genes could provide additional information on how the anterior posterior axis formation is established in Tribolium embryos as compared to Drosophila Together, our data will facilitate studies of comparative developmental biology as well as help identify candidate genes for identifying cis-elements to drive transgenic constructs.

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.

Zelda and the maternal-to-zygotic transition in cockroaches

In the endopterygote Drosophila melanogaster, Zelda is an activator of the zygotic genome during the maternal-to-zygotic transition (MZT). Zelda binds cis-regulatory elements (TAGteam heptamers), making chromatin accessible for gene transcription. Zelda has been studied in other endopterygotes: Apis mellifera and Tribolium castaneum, and the paraneopteran Rhodnius prolixus. We studied Zelda in the cockroach Blattella germanica, a hemimetabolan, short germ-band, and polyneopteran species. B. germanica Zelda has the complete set of functional domains, which is typical of species displaying ancestral features concerning embryogenesis. Interestingly, we found D. melanogaster TAGteam heptamers in the B. germanica genome. The canonical one, CAGGTAG, is present at a similar proportion in the genome of these two species and in the genome of other insects, suggesting that the genome admits as many CAGGTAG motifs as its length allows. Zelda-depleted embryos of B. germanica show defects involving blastoderm formation and abdomen development, and genes contributing to these processes are down-regulated. We conclude that in B. germanica, Zelda strictly activates the zygotic genome, within the MZT, a role conserved in more derived endopterygote insects. In B. germanica, zelda is expressed during MZT, whereas in D. melanogaster and T. castaneum it is expressed beyond this transition. In these species and A. mellifera, Zelda has functions even in postembryonic development. The expansion of zelda expression beyond the MZT in endopterygotes might be related with the evolutionary innovation of holometabolan metamorphosis. DATABASES: The RNA-seq datasets of B. germanica, D. melanogaster, and T. castaneum are accessible at the GEO databases GSE99785, GSE18068, GSE63770, and GSE84253. In addition, the RNA-seq library from T. castaneum adult females is available at SRA: SRX021963. The B. germanica reference genome is available as BioProject PRJNA203136.

miRNAs and their roles in KSHV pathogenesis

Kaposi's sarcoma-associated herpesvirus (KSHV) is the etiological agent of Kaposi's sarcoma (KS), primary effusion lymphoma (PEL), and multicentric Castleman Disease (MCD). Recent mechanistic advances have discerned the importance of microRNAs in the virus-host relationship. KSHV has two modes of replication: lytic and latent phase. KSHV entry into permissive cells, establishment of infection, and maintenance of latency are contingent upon successful modulation of the host miRNA transcriptome. Apart from host cell miRNAs, KSHV also encodes viral miRNAs. Among various cellular and molecular targets, miRNAs are appearing to be key players in regulating viral pathogenesis. Therefore, the use of miRNAs as novel therapeutics has gained considerable attention as of late. This innovative approach relies on either mimicking miRNA species by identical oligonucleotides, or selective silencing of miRNA with specific oligonucleotide inhibitors. Here, we provide an overview of KSHV pathogenesis at the molecular level with special emphasis on the various roles miRNAs play during virus infection.

Multiple functions of miR-8-3p in the development and metamorphosis of the red flour beetle, Tribolium castaneum

The microRNA miR-8-3p is conserved among insects and closely involved in development and immunity, but its functions in vivo are unexplored in the red flour beetle, Tribolium castaneum. Here, we show that miR-8-3p was highly expressed in late larva and early adult stages, as determined by quantitative real-time PCR. It was enriched in the fat body and cuticle in late larval tissues and abundant in the head and cuticle in early adult tissues, indicating this microRNA plays important roles during T. castaneum development. Specific inhibition of miR-8-3p in late larvae led to metamorphosis defects in the development of wings, eyes, legs and embryo. Moreover, a series of genes related to organism development were identified as miR-8-3p targets by computational prediction and microRNA-messenger RNA interaction validation, including Wingless, Eyg, Fpps and Sema-1a. These genes were critical for the regulation of the larva-to-adult transition. Eyg, as a functional target of miR-8-3p, participates in eye development, which was further confirmed by luciferase assay and loss-of-function analyses. In brief, miR-8-3p is broadly involved in the development of wings, eyes and legs through its target genes and has extensive regulatory roles during T. castaneum development.

microRNAs in Macrobrachium olfersii embryos: Identification, their biogenesis components and potential targets

In embryonic development, microRNAs (miRNAs) regulate the complex gene expression associated with the complexity of embryogenesis. Today, few studies have been conducted on the identification of miRNAs and components of miRNA biogenesis on embryonic development in crustaceans, especially in prawns. In this context, the aim of this study was to identify in silico components of miRNA biogenesis, and miRNAs and potential target genes during embryonic development in the prawn Macrobrachium olfersii through small RNAs and transcriptome analyses. Using the miRDeep2 program, we identified 17 miRNA precursors in M. olfersii, which seven (miR-9, miR-10, miR-92, miR-125, miR-305, miR-1175, and miR-2788) were reported in the miRBase database, indicating high evolutionary conservation of these sequences among animals. The other 10 miRNAs of M. olfersii were novel miRNAs and only similar to Macrobrachium niponnense miRNAs, indicating genus-specific miRNAs. In addition, eight key components of miRNA biogenesis (DROSHA, PASHA/DGCR8, XPO5, RAN, DICER, TRBP2, AGO, and PIWI) were identified in M. olfersii embryos unigenes. In the annotation of miRNA targets, 516 genes were similar to known sequences in the GenBank database. Regarding the conserved miRNAs, we verified that they were differentially expressed during embryonic development in M. olfersii. In conclusion, this is the first study that identifies conserved and novel miRNAs in the prawn M. olfersii with some miRNA target genes involved in embryonic development. Our results will allow further studies on the function of these miRNAs and miRNA biogenesis components during embryonic development in M. olfersii and other prawns of commercial interest.

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.

A six-microRNA signature predicts survival of patients with uterine corpus endometrial carcinoma

Uterine corpus endometrial carcinoma (UCEC) is one of the most common female gynecological malignant tumors that threaten women health seriously. MicroRNAs (miRNAs) has been proved to play critical roles in tumor pathogenesis and malignant progression. In this study, we aimed to explore a novel signature of microRNA expression for predicting the overall survival (OS) of patients with UCEC. The genome-wide miRNA expression profiles and relevant clinical characteristics of 348 patients with UCEC were downloaded from the Cancer Genome Atlas (TCGA) data portal and analyzed comprehensively. A total of 144 miRNAs were confirmed to be expressed differentially in tumor tissues. Among them, 6 miRNAs (hsa-mir-15a.MIMAT0000068, hsa-mir-142.MIMAT0000433, hsa-mir-142.MIMAT0000434, hsa-mir-3170.MIMAT0015045, hsa-mir-1976.MIMAT0009451, and hsa-mir-146a.MIMAT0000449) were validated to be significantly correlated with the OS of patients with UCEC. The risk indictor established by the 6-microRNA signature was proved be an independent prognostic factor (Hazard ratio = 0.391; 95% CI: 0.195-0.783; P = 0.008). In conclusion, we identified miRNAs that were correlated with the occurrence and progression of UCEC and established a 6-microRNA expression signature as a predictor for the OS of patients with UCEC.

Daphnia magna microRNAs respond to nutritional stress and ageing but are not transgenerational

Maternal effects, where the performance of offspring is determined by the condition of their mother, are widespread and may in some cases be adaptive. The crustacean Daphnia magna shows strong maternal effects: offspring size at birth and other proxies for fitness are altered when their mothers are older or when mothers have experienced dietary restriction. The mechanisms for this transgenerational transmission of maternal experience are unknown, but could include changes in epigenetic patterning. MicroRNAs (miRNAs) are regulators of gene expression that have been shown to play roles in intergenerational information transfer, and here, we test whether miRNAs are involved in D. magna maternal effects. We found that miRNAs were differentially expressed in mothers of different ages or nutritional state. We then examined miRNA expression in their eggs, their adult daughters and great granddaughters, which did not experience any treatments. The maternal (treatment) generation exhibited differential expression of miRNAs, as did their eggs, but this was reduced in adult daughters and lost by great granddaughters. Thus, miRNAs are a component of maternal provisioning, but do not appear to be the cause of transgenerational responses under these experimental conditions. MicroRNAs may act in tandem with egg provisioning (e.g., with carbohydrates or fats), and possibly other small RNAs or epigenetic modifications.

LncRNA SNHG12 ameliorates brain microvascular endothelial cell injury by targeting miR-199a

Long non-coding RNAs regulate brain microvascular endothelial cell death, the inflammatory response and angiogenesis during and after ischemia/reperfusion and oxygen-glucose deprivation/reoxygenation (OGD/R) insults. The long non-coding RNA, SNHG12, is upregulated after ischemia/reperfusion and OGD/R in microvascular endothelial cells of the mouse brain. However, its role in ischemic stroke has not been studied. We hypothesized that SNHG12 positively regulates ischemic stroke, and therefore we investigated its mechanism of action. We established an OGD/R mouse cell model to mimic ischemic stroke by exposing brain microvascular endothelial cells to OGD for 0, 2, 4, 8, 16 or 24 hours and reoxygenation for 4 hours. Quantitative real-time polymerase chain reaction showed that SNHG12 levels in brain microvascular endothelial cells increased with respect to OGD exposure time. Brain microvascular endothelial cells were transfected with pcDNA-control, pcDNA-SNHG12, si-control, or si-SNHG12. After exposure to OGD for 16 hours, these cells were then analyzed by 3-(4,5-dimethyl-2-thiazolyl)-2,5-diphenyl-2-H-tetrazolium bromide, trypan blue exclusion, western blot, and capillary-like tube formation assays. Overexpression of SNHG12 inhibited brain microvascular endothelial cell death and the inflammatory response but promoted angiogenesis after OGD/R, while SNHG12 knockdown had the opposite effects. miR-199a was identified as a target of SNHG12, and SNHG12 overexpression reversed the effect of miR-199a on brain microvascular endothelial cell death, the inflammatory response, and angiogenesis. These findings suggest that SNHG12 suppresses endothelial cell injury induced by OGD/R by targeting miR-199a.

Comparative analysis of miRNA expression during the development of insects of different metamorphosis modes and germ-band types

Background

Do miRNAs contribute to specify the germ-band type and the body structure in the insect embryo? Our goal was to address that issue by studying the changes in miRNA expression along the ontogeny of the German cockroach Blattella germanica, which is a short germ-band and hemimetabolan species.

Results

We sequenced small RNA libraries representing 11 developmental stages of B. germanica ontogeny (with especial emphasis on embryogenesis) and the changes in miRNA expression were examined. Data were compared with equivalent data for two long germ-band holometabolan species Drosophila melanogaster and Drosophila virilis, and the short germ-band holometabolan species Tribolium castaneum. The identification of B. germanica embryo small RNA sequences unveiled miRNAs not detected in previous studies, such as those of the MIR-309 family and 54 novel miRNAs. Four main waves of miRNA expression were recognized (with most miRNA changes occurring during the embryonic stages): the first from day 0 to day 1 of embryogenesis, the second during mid-embryogenesis (days 0–6), the third (with an acute expression peak) on day 2 of embryonic development, and the fourth during post-embryonic development. The second wave defined the boundaries of maternal-to-zygotic transition, with maternal mRNAs being cleared, presumably by Mir-309 and associated scavenger miRNAs.

Conclusion

miRNAs follow well-defined patterns of expression over hemimetabolan ontogeny, patterns that are more diverse during embryonic development than during the nymphal stages. The results suggest that miRNAs play important roles in the developmental transitions between the embryonic stages of development (starting with maternal loading), during which they might influence the germ-band type and metamorphosis mode.

Electronic supplementary material

The online version of this article (10.1186/s12864-017-4177-5) contains supplementary material, which is available to authorized users.

Abundant expression of somatic transposon-derived piRNAs throughout Tribolium castaneum embryogenesis

BACKGROUND

Piwi-interacting RNAs (piRNAs) are a class of short (~26-31-nucleotide) non-protein-coding RNAs expressed in the metazoan germline. The piRNA pathway in arthropods is best understood in the ovary of Drosophila melanogaster, where it acts to silence active transposable elements (TEs). Maternal loading of piRNAs in oocytes is further required for the inheritance of piRNA-mediated transposon defence. However, our understanding of the diversity, evolution and function of the piRNA complement beyond drosophilids is limited. The red flour beetle, Tribolium castaneum, is an emerging model organism separated from Drosophila by ~ 350 million years of evolution that displays a number of features ancestral to arthropods, including short germ embryogenesis. Here, we characterize the maternally deposited and zygotically expressed small RNA and mRNA complements throughout T. castaneum embryogenesis.

RESULTS

We find that beetle oocytes and embryos of all stages are abundant in heterogeneous ~ 28-nucleotide RNAs. These small RNAs originate from discrete genomic loci enriched in TE sequences and display the molecular signatures of transposon-derived piRNAs. In addition to the maternally loaded primary piRNAs, Tribolium embryos produce secondary piRNAs by the cleavage of zygotically activated TE transcripts via the ping-pong mechanism. The two Tribolium piRNA pathway effector proteins, Tc-Piwi/Aub and Tc-Ago3, are also expressed throughout the soma of early embryos.

CONCLUSIONS

Our results show that the piRNA pathway in Tribolium is not restricted to the germline, but also operates in the embryo and may act to antagonize zygotically activated transposons. Taken together, these data highlight a functional divergence of the piRNA pathway between insects.

Evaluation of high-throughput isomiR identification tools: illuminating the early isomiRome of Tribolium castaneum

Background

MicroRNAs carry out post-transcriptional gene regulation in animals by binding to the 3' untranslated regions of mRNAs, causing their degradation or translational repression. MicroRNAs influence many biological functions, and dysregulation can therefore disrupt development or even cause death. High-throughput sequencing and the mining of animal small RNA data has shown that microRNA genes can yield differentially expressed isoforms, known as isomiRs. Such isoforms are particularly relevant during early development, and the extension or truncation of the 5' end can change the profile of mRNA targets compared to the original mature sequence. We used the publicly available small RNA dataset of the model beetle Tribolium castaneum to create the first comparative isomiRome of early developmental stages in this species. Standard microRNA analysis software does not specifically account for isomiRs. We therefore carried out the first comparative evaluation of the specialized tools isomiRID, isomiR-SEA and miraligner, which can be downloaded for local use and can handle next generation sequencing data.

Results

We compared the performance of isomiRID, isomiR-SEA and miraligner using simulated Illumina HiSeq2000 and MiSeq data to test the impact of technical errors. We also created artificial microRNA isoforms to determine the effect of biological variants on the performance of each algorithm. We found that isomiRID achieved the best true positive rate among the three algorithms, but only accounted for one mutation at a time. In contrast, miraligner reported all variations simultaneously but with 78% sensitivity, yielding isomiRs with 3' or 5' deletions. Finally, isomiR-SEA achieved a sensitivity of 25–33% when the seed region was mutated or partly deleted, but was the only tool that could accommodate more than one mismatch. Using the best tool, we performed a complete isomiRome analysis of the early developmental stages of T. castaneum.

Conclusions

Our findings will help researchers to select the most suitable isomiR analysis tools for their experiments. We confirmed the dynamic expression of 3′ non-template isomiRs and expanded the isomiRome by all known isomiR modifications during the early development of T. castaneum.

Electronic supplementary material

The online version of this article (doi:10.1186/s12859-017-1772-z) contains supplementary material, which is available to authorized users.

What Drives Positive Selection in the Drosophila piRNA Machinery? The Genomic Autoimmunity Hypothesis

In animals, PIWI-interacting RNAs (piRNAs) play a crucial role in genome defense. Moreover, because piRNAs can be maternally transmitted, they contribute to the epigenetic profile of inheritance. Multiple studies, especially in Drosophila, have demonstrated that the machinery of piRNA biogenesis is often the target of positive selection. Because transposable elements (TEs) are a form of genetic parasite, positive selection in the piRNA machinery is often explained by analogy to the signatures of positive selection commonly observed in genes that play a role in host-parasite dynamics. However, the precise mechanisms that drive positive selection in the piRNA machinery are not known. In this review, we outline several mechanistic models that might explain pervasive positive selection in the piRNA machinery of Drosophila species. We propose that recurrent positive selection in the piRNA machinery can be partly explained by an ongoing tension between selection for sensitivity required by genome defense and selection for specificity to avoid the off-target effects of maladaptive genic silencing by piRNA.

The microRNA toolkit of insects

Is there a correlation between miRNA diversity and levels of organismic complexity? Exhibiting extraordinary levels of morphological and developmental complexity, insects are the most diverse animal class on earth. Their evolutionary success was in particular shaped by the innovation of holometabolan metamorphosis in endopterygotes. Previously, miRNA evolution had been linked to morphological complexity, but astonishing variation in the currently available miRNA complements of insects made this link unclear. To address this issue, we sequenced the miRNA complement of the hemimetabolan Blattella germanica and reannotated that of two other hemimetabolan species, Locusta migratoria and Acyrthosiphon pisum, and of four holometabolan species, Apis mellifera, Tribolium castaneum, Bombyx mori and Drosophila melanogaster. Our analyses show that the variation of insect miRNAs is an artefact mainly resulting from poor sampling and inaccurate miRNA annotation, and that insects share a conserved microRNA toolkit of 65 families exhibiting very low variation. For example, the evolutionary shift toward a complete metamorphosis was accompanied only by the acquisition of three and the loss of one miRNA families.

Identification and Characterization of Sex-Biased MicroRNAs in Bactrocera dorsalis (Hendel)

MicroRNAs (miRNAs) are a class of endogenous small non-coding RNAs that regulate various biological processes including sexual dimorphism. The oriental fruit fly Bactrocera dorsalis is one of the most destructive agricultural insect pests in many Asian countries. However, no miRNAs have been identified from the separate sex and gonads to elucidate sex gonad differentiation in B. dorsalis. In this study, we constructed four small RNA libraries from whole body of females, males (except ovaries and testes) and ovaries, testes of B. dorsalis for deep sequencing. The data analysis revealed 183 known and 120 novel miRNAs from these libraries. 18 female-biased and 16 male-biased miRNAs that may be involved in sexual differentiation were found by comparing the miRNA expression profiles in the four libraries. Using a bioinformatic approach, we predicted doublesex (dsx) as a target gene of the female-biased miR-989-3p which is considered as the key switch gene in the sex determination of tephritid insects. This study reveals the first miRNA profile related to the sex differentiation and gives a first insight into sex differences in miRNA expression of B. dorsalis which could facilitate studies of the reproductive organ specific roles of miRNAs.

Selection Against Maternal microRNA Target Sites in Maternal Transcripts

In animals, before the zygotic genome is expressed, the egg already contains gene products deposited by the mother. These maternal products are crucial during the initial steps of development. In Drosophila melanogaster, a large number of maternal products are found in the oocyte, some of which are indispensable. Many of these products are RNA molecules, such as gene transcripts and ribosomal RNAs. Recently, microRNAs (small RNA gene regulators) have been detected early during development and are important in these initial steps. The presence of some microRNAs in unfertilized eggs has been reported, but whether they have a functional impact in the egg or early embryo has not being explored. I have extracted and sequenced small RNAs from Drosophila unfertilized eggs. The unfertilized egg is rich in small RNAs and contains multiple microRNA products. Maternal microRNAs often are encoded within the intron of maternal genes, suggesting that many maternal microRNAs are the product of transcriptional hitchhiking. Comparative genomics analyses suggest that maternal transcripts tend to avoid target sites for maternal microRNAs. I also developed a microRNA target mutation model to study the functional impact of polymorphisms at microRNA target sites. The analysis of Drosophila populations suggests that there is selection against maternal microRNA target sites in maternal transcripts. A potential role of the maternal microRNA mir-9c in maternal-to-zygotic transition is also discussed. In conclusion, maternal microRNAs in Drosophila have a functional impact in maternal protein-coding transcripts.