Sex-biased expression of microRNAs in Drosophila melanogaster

The genome and sex-dependent responses to temperature in the common yellow butterfly, Eurema hecabe

BACKGROUND

Lepidoptera (butterflies and moths) is one of the most geographically widespread insect orders in the world, and its species play important and diverse ecological and applied roles. Climate change is one of the biggest challenges to biodiversity this century, and lepidopterans are vulnerable to climate change. Temperature-dependent gene expression differences are of relevance under the ongoing climate crisis. However, little is known about how climate affects gene expression in lepidopterans and the ecological consequences of this, particularly with respect to genes with biased expression in one of the sexes. The common yellow butterfly, Eurema hecabe (Family Pieridae), is one of the most geographically widespread lepidopterans that can be found in Asia, Africa, and Australia. Nevertheless, what temperature-dependent effects there may be and whether the effects differ between the sexes remain largely unexplored.

RESULTS

Here, we generated high-quality genomic resources for E. hecabe along with transcriptomes from eight developmental stages. Male and female butterflies were subjected to varying temperatures to assess sex-specific gene expression responses through mRNA and microRNA transcriptomics. We find that there are more temperature-dependent sex-biased genes in females than males, including genes that are involved in a range of biologically important functions, highlighting potential ecological impacts of increased temperatures. Further, by considering available butterfly data on sex-biased gene expression in a comparative genomic framework, we find that the pattern of sex-biased gene expression identified in E. hecabe is highly species-specific, rather than conserved across butterfly species, suggesting that sex-biased gene expression responses to climate change are complex in butterflies.

CONCLUSIONS

Our study lays the foundation for further understanding of differential responses to environmental stress in a widespread lepidopteran model and demonstrates the potential complexity of sex-specific responses of lepidopterans to climate change.

MicroRNAs reshape the immunity of insects in response to bacterial infection

The interaction between bacteria and insects can significantly impact a wide range of different areas because bacteria and insects are widely distributed around the globe. The bacterial-insect interactions have the potential to directly affect human health since insects are vectors for disease transmission, and their interactions can also have economic consequences. In addition, they have been linked to high mortality rates in economically important insects, resulting in substantial economic losses. MicroRNAs (miRNAs) are types of non-coding RNAs involved in regulating gene expression post-transcriptionally. The length of miRNAs ranges from 19 to 22 nucleotides. MiRNAs, in addition to their ability to exhibit dynamic expression patterns, have a diverse range of targets. This enables them to govern various physiological activities in insects, like innate immune responses. Increasing evidence suggests that miRNAs have a crucial biological role in bacterial infection by influencing immune responses and other mechanisms for resistance. This review focuses on some of the most recent and exciting discoveries made in recent years, including the correlation between the dysregulation of miRNA expression in the context of bacterial infection and the progression of the infection. Furthermore, it describes how they profoundly impact the immune responses of the host by targeting the Toll, IMD, and JNK signaling pathways. It also emphasizes the biological function of miRNAs in regulating immune responses in insects. Finally, it also discusses current knowledge gaps about the function of miRNAs in insect immunity, in addition to areas that require more research in the future.

The role of micro RNAs (miRNAs) in the regulation of Drosophila melanogaster's innate immunity

MicroRNAs (miRNAs) are a class of small non-coding RNAs ~19-22 nt long which post-transcriptionally regulate gene expression. Their ability to exhibit dynamic expression patterns coupled with their wide variety of targets allows miRNAs to regulate many processes, including the innate immune response of Drosophila melanogaster. Recent studies have identified miRNAs in Drosophila which are differentially expressed during infection with different pathogens as well as miRNAs that may affect immune signalling when differentially expressed. This review provides an overview of miRNAswhich have been identified to play a role in the immune response of Drosophila through targeting of the Toll and IMD signalling pathways and other immune processes. It will also explore the role of miRNAs in fine-tuning the immune response in Drosophila and highlight current gaps in knowledge regarding the role of miRNAs in immunity and areas for further research.

The small RNA landscape is stable with age and resistant to loss of dFOXO signaling in Drosophila

Aging can be defined as the progressive loss of physiological homeostasis that leads to a decline in cellular and organismal function. In recent years, it has become clear that small RNA pathways play a role in aging and aging related phenotypes. Small RNA pathways regulate many important processes including development, cellular physiology, and innate immunity. The pathways illicit a form of posttranscriptional gene regulation that relies on small RNAs bound by the protein components of the RNA-induced silencing complexes (RISCs), which inhibit the expression of complementary RNAs. In Drosophila melanogaster, Argonaute 1 (Ago1) is the core RISC component in microRNA (miRNA) silencing, while Argonaute 2 (Ago2) is the core RISC component in small interfering RNA (siRNA) silencing. The expression of Ago1 and Ago2 is regulated by stress response transcription factor Forkhead box O (dFOXO) increasing siRNA silencing efficiency. dFOXO plays a role in multiple stress responses and regulates pathways important for longevity. Here we use a next-generation sequencing approach to determine the effects of aging on small RNA abundance and RISC loading in male and female Drosophila. In addition, we examine the impact of the loss of dFOXO on these processes. We find that the relative abundance of the majority of small RNAs does not change with age. Additionally, under normal growth conditions, the loss of dFOXO has little effect on the small RNA landscape. However, we observed that age affects loading into RISC for a small number of miRNAs.

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.

The developmentally dynamic microRNA transcriptome of Glossina pallidipes tsetse flies, vectors of animal trypanosomiasis

Summary

MicroRNAs (miRNAs) are single stranded gene regulators of 18-25 bp in length. They play a crucial role in regulating several biological processes in insects. However, the functions of miRNA in Glossina pallidipes, one of the biological vectors of African animal trypanosomosis in sub-Saharan Africa, remain poorly characterized. We used a combination of both molecular biology and bioinformatics techniques to identify miRNA genes at different developmental stages (larvae, pupae, teneral and reproductive unmated adults, gravid females) and sexes of G. pallidipes. We identified 157 mature miRNA genes, including 12 novel miRNAs unique to G. pallidipes. Moreover, we identified 93 miRNA genes that were differentially expressed by sex and/or in specific developmental stages. By combining both miRanda and RNAhybrid algorithms, we identified 5550 of their target genes. Further analyses with the Gene Ontology term and KEGG pathways for these predicted target genes suggested that the miRNAs may be involved in key developmental biological processes. Our results provide the first repository of G. pallidipes miRNAs across developmental stages, some of which appear to play crucial roles in tsetse fly development. Hence, our findings provide a better understanding of tsetse biology and a baseline for exploring miRNA genes in tsetse flies.

Availability and implementation

Raw sequence data are available from NCBI Sequence Read Archives (SRA) under Bioproject accession number PRJNA590626.

Supplementary information

Supplementary data are available at Bioinformatics Advances online.

Evolutionary history of sexual selection affects microRNA profiles in Drosophila sperm

The presence of small RNAs in sperm is a relatively recent discovery and little is currently known about their importance and functions. Environmental changes including social conditions and dietary manipulations are known to affect the composition and expression of some small RNAs in sperm and may elicit a physiological stress response resulting in an associated change in gamete miRNA profiles. Here, we tested how microRNA profiles in sperm are affected by variation in both sexual selection and dietary regimes in Drosophila melanogaster selection lines. The selection lines were exposed to standard versus low yeast diet treatments and three different population sex ratios (male-biased, female-biased or equal sex) in a full-factorial design. After 38 generations of selection, all males were maintained on their selected diet and in a common garden male-only environment prior to sperm sampling. We performed transcriptome analyses on miRNAs in purified sperm samples. We found 11 differentially expressed miRNAs with the majority showing differences between male- and female-biased lines. Dietary treatment only had a significant effect on miRNA expression levels in interaction with sex ratio. Our findings suggest that long-term adaptation may affect miRNA profiles in sperm and that these may show varied interactions with short-term environmental changes. This article is protected by copyright. All rights reserved.

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.

Identification and Characterization of MicroRNAs in Gonads of Helicoverpa armigera (Lepidoptera: Noctuidae)

Simple Summary

For most insects, the development of the testis and ovary directly determines their reproductive ability. The cotton bollworm, Helicoverpa armigera (Hübner), is a polyphagous crop pest of the Lepidoptera Noctuidae. Owing to its broad range of host plants and strong fertility, H. armigera causes huge economic losses to agricultural production. Acting as a type of post-transcriptional regulatory factor, miRNAs participate in the gonadal development and reproductive regulation of H. arimgera. Our study uses H. armigera as a research object to identify and characterize the miRNAs and study their potential functions in the testis and ovary of this destructive crop pest. A total of 7,592,150 and 8,815,237 clean reads were obtained by constructing small RNA libraries of the testis and ovary, respectively. Length distribution analysis showed that the main types of small RNAs in the testis and ovary were different. Among the 74 known miRNAs, 60 miRNAs existed in the ovary, and 72 existed in the testis. Gene Ontology (GO) and KEGG pathway analyses indicated that the 8 gonad-biased differentially expressed miRNAs (miR-989a, miR-263-5p, miR-34, miR-2763, miR-998, miR-2c, miR-2765, and miR-252a-5p) had many target transcripts involved in the reproduction process.

Abstract

The high fecundity of the most destructive pest Helicoverpa armigera and its great resistance risk to insecticides and Bt crops make the reproductive-destruction-based control of this pest extremely appealing. To find suitable targets for disruption of its reproduction, we observed the testis and ovary development of H. armigera and conducted deep sequencing of the ovary and testis small RNAs of H. armigera and quantitative RT-PCR (RT-qPCR) validation to identify reproduction-related micro RNAs (miRNAs). A total of 7,592,150 and 8,815,237 clean reads were obtained from the testis and ovary tissue, respectively. After further analysis, we obtained 173 novel and 74 known miRNAs from the two libraries. Among the 74 known miRNAs, 60 miRNAs existed in the ovary and 72 existed in the testis. Further RT-qPCR validation of 5 miRNAs from the ovary and 6 miRNAs from the testis confirmed 8 of them were indeed ovary- (miR-989a, miR-263-5p, miR-34) or testis-biased (miR-2763, miR-998, miR-2c, miR-2765, miR-252a-5p). The 8 ovary- or testis-biased miRNAs had a total of 30,172 putative non-redundant target transcripts, as predicted by miRanda and RNAhybrid. Many of these target transcripts are assigned to reproduction-related GO terms (e.g., oocyte maturation, vitellogenesis, spermatogenesis) and are members of multiple reproduction-related KEGG pathways, such as the JAK-STAT signaling pathway, oocyte meiosis, the insulin signaling pathway, and insect hormone biosynthesis. These results suggest that the 8 gonad-biased miRNAs play important roles in reproduction and may be used as the targets for the development of reproductive-destruction-based control of H. armigera and, possibly, other lepidopteran pests.

Genome-wide identification of miRNAs and target regulatory network in the invasive ectoparasitic mite Varroa destructor

Varroa destructor is an ectoparasite mite that attacks bees leading to colony disorders worldwide. microRNAs (miRNAs) are key molecules used by eukaryotes to post-transcriptional control of gene expression. Nevertheless, still lack information aboutV. destructor miRNAs and its regulatory networks. Here, we used an integrative strategy to characterize the miRNAs in the V. destructor mite. We identified 310 precursors that give rise to 500 mature miRNAs, which 257 are likely mite-specific elements. miRNAs showed canonical length ranging between 18 and 25 nucleotides and 5' uracil preference. Top 10 elements concentrated over 80% of total miRNA expression, with bantam alone representing ~50%. We also detected non-templated bases in precursor-derived small RNAs, indicative of miRNA post-transcriptional regulatory mechanisms. Finally, we note that conserved miRNAs control similar processes in different organisms, suggesting a conservative role. Altogether, our findings contribute to the better understanding of the mite biology that can assist future studies on varroosis control.

miRNAs of Aedes aegypti (Linnaeus 1762) conserved in six orders of the class Insecta

Aedes aegypti L. is the most important vector of arboviruses such as dengue, Zika, chikungunya, Mayaro, and yellow fever, which impact millions of people's health per year. MicroRNA profile has been described in some mosquito species as being important for biological processes such as digestion of blood, oviposition, sexual differentiation, insecticide resistance, and pathogens dissemination. We identified the miRNAs of Ae. aegypti females, males and eggs of a reference insecticide susceptible strain New Orleans and compared them with those other insects to determine miRNA fingerprint by new-generation sequencing. The sequences were analyzed using data mining tools and categorization, followed by differential expression analysis and conservation with other insects. A total of 55 conserved miRNAs were identified, of which 34 were of holometabolous insects and 21 shared with hemimetabolous insects. Of these miRNAs, 32 had differential expression within the stages analyzed. Three predominant functions of miRNA were related to embryonic development regulation, metamorphosis, and basal functions. The findings of this research describe new information on Ae. aegypti physiology which could be useful for the development of new control strategies, particularly in mosquito development and metamorphosis processes.

A comparative genomic approach using mouse and fruit fly data to discover genes involved in testis function in hymenopterans with a focus on Nasonia vitripennis

Background

Spermatogenesis appears to be a relatively well-conserved process even among distantly related animal taxa such as invertebrates and vertebrates. Although Hymenopterans share many characteristics with other organisms, their complex haplodiploid reproduction system is still relatively unknown. However, they serve as a complementary insect model to Drosophila for studying functional male fertility. In this study, we used a comparative method combining taxonomic, phenotypic data and gene expression to identify candidate genes that could play a significant role in spermatogenesis in hymenopterans.

Results

Of the 546 mouse genes predominantly or exclusively expressed in the mouse testes, 36% had at least one ortholog in the fruit fly. Of these genes, 68% had at least one ortholog in one of the six hymenopteran species we examined. Based on their gene expression profiles in fruit fly testes, 71 of these genes were hypothesized to play a marked role in testis function. Forty-three of these 71 genes had an ortholog in at least one of the six hymenopteran species examined, and their enriched GO terms were related to the G2/M transition or to cilium organization, assembly, or movement. Second, of the 379 genes putatively involved in male fertility in Drosophila, 224 had at least one ortholog in each of the six Hymenoptera species. Finally, we showed that 199 of these genes were expressed in early pupal testis in Nasonia vitripennis; 86 exhibited a high level of expression, and 54 displayed modulated expression during meiosis.

Conclusions

In this study combining phylogenetic and experimental approaches, we highlighted genes that may have a major role in gametogenesis in hymenopterans; an essential prerequisite for further research on functional importance of these genes.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12862-021-01825-6.

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.

Divergence in Transcriptional and Regulatory Responses to Mating in Male and Female Fruitflies

Mating induces extensive physiological, biochemical and behavioural changes in female animals of many taxa. In contrast, the overall phenotypic and transcriptomic consequences of mating for males, hence how they might differ from those of females, are poorly described. Post mating responses in each sex are rapidly initiated, predicting the existence of regulatory mechanisms in addition to transcriptional responses involving de novo gene expression. That post mating responses appear different for each sex also predicts that the genome-wide signatures of mating should show evidence of sex-specific specialisation. In this study, we used high resolution RNA sequencing to provide the first direct comparisons of the transcriptomic responses of male and female Drosophila to mating, and the first comparison of mating-responsive miRNAs in both sexes in any species. As predicted, the results revealed the existence of sex- and body part-specific mRNA and miRNA expression profiles. More genes were differentially expressed in the female head-thorax than the abdomen following mating, whereas the opposite was true in males. Indeed, the transcriptional profile of male head-thorax tissue was largely unaffected by mating, and no differentially expressed genes were detected at the most stringent significance threshold. A subset of ribosomal genes in females were differentially expressed in both body parts, but in opposite directions, consistent with the existence of body part-specific resource allocation switching. Novel, mating-responsive miRNAs in each sex were also identified, and a miRNA-mRNA interactions analysis revealed putative targets among mating-responsive genes. We show that the structure of genome-wide responses by each sex to mating is strongly divergent, and provide new insights into how shared genomes can achieve characteristic distinctiveness.

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.

Gonadal microRNA Expression Profiles and Their Potential Role in Sex Differentiation and Gonadal Maturation of Mud Crab Scylla paramamosain

Although the sexual dimorphism in terms of gonadal development and gametogenesis of mud crab has been described, the internal regulating mechanism and sex differentiation process remain unclear. A comparative gonadal miRNA transcriptomic study was conducted to identify miRNAs that are differentially expressed between testes and ovaries, and potentially uncover miRNAs that might be involved in sex differentiation and gonadal maturation mechanisms of mud crabs (Scylla paramamosain). A total of 10 known miRNAs and 130 novel miRNAs were identified, among which 54 were differentially expressed. Target gene prediction revealed a significant enrichment in 30 KEGG pathways, including some reproduction-related pathways, e.g. phosphatidylinositol signalling system and inositol phosphate metabolism pathways. Further analysis on six differentially expressed known miRNAs, six differentially expressed novel miRNAs and their reproduction-related putative target genes shows that both miRNAs and putative target genes showed stage-specific expression during gonadal maturation, suggesting their potential regulatory roles in sex differentiation and reproductive development. This study reveals the sex-biased miRNA profile and establishes a solid foundation for understanding the sex differentiation and gonadal maturation mechanisms of S. paramamosain.

Characterization of the salmon louse Lepeophtheirus salmonis miRNome: Sex-biased differences related to the coding and non-coding RNA interplay

The salmon louse Lepeophtheirus salmonis is a marine ectoparasite that has a detrimental impact on salmon farms. Genomic knowledge of adult stages is critical to understand the reproductive success and lifecycle completion of this species. Here, we report a comprehensive characterization of the L. salmonis miRNome with emphasis on the sex-differences of the parasite. Small-RNA sequencing was conducted on males and females, and mRNA-sequencing was also conducted to identify miRNA-targets at these stages. Based on bioinformatics analyses, 3101 putative miRNAs were found in L. salmonis, including precursors and variants. The most abundant and over-expressed miRNAs belonged to the bantam, mir-100, mir-1, mir-263a and mir-276 families, while the most differentially expressed mRNAs corresponded to genes related to reproduction and other biological processes involved in cell-differentiation. Target analyses revealed that the most up-regulated miRNAs in males can act by inhibiting the expression of genes related to female differentiation such as vitellogenin genes. Target prediction and expression patterns suggested a pivotal role of miRNAs in the reproductive development of L. salmonis.

Characterization and comparative analysis of microRNAs in the rice pest Sogatella furcifera

MicroRNAs (miRNAs) are a class of endogenous regulatory RNA molecules 21-24 nucleotides in length that act as functional regulators of post-transcriptional repression of messenger RNA. We report the identification and characterization of a conserved miRNA and 171 novel miRNAs in the migratory rice pest Sogatella furcifera by deep sequencing, which were observed to be biased towards female adults of the insect, modulating the functionality and targets of the miRNAs in sex differentiation. A switch in arm usage was also observed in 9 miRNA when compared to the insect ancestor during insect evolution. The miRNA loci showed high 5’ fidelity in both miRNA and star species and about 93.4% of WBPH miRNAs conserved within non-planthopper species were homologous with planthopper species. The novel miRNAs identified in this study provide a better understanding of the sRNA and the regulatory role of miRNA in sexual dimorphism and alteration in the expression or function of miRNAs in the rice pest.

Differentially expressed microRNAs associated with changes of transcript levels in detoxification pathways and DDT-resistance in the Drosophila melanogaster strain 91-R

Dichloro-diphenyl-trichloroethane (DDT) resistance among arthropod species is a model for understanding the molecular adaptations in response to insecticide exposures. Previous studies reported that DDT resistance may involve one or multiple detoxification genes, such as cytochrome P450 monooxygenases (P450s), glutathione S-transferases (GSTs), esterases, and ATP binding cassette (ABC) transporters, or changes in the voltage-sensitive sodium channel. However, the possible involvement of microRNAs (miRNAs) in the post-transcriptional regulation of genes associated with DDT resistance in the Drosophila melanogaster strain 91-R remains poorly understood. In this study, the majority of the resulting miRNAs discovered in small RNA libraries from 91-R and the susceptible control strain, 91-C, ranged from 16-25 nt, and contained 163 precursors and 256 mature forms of previously-known miRNAs along with 17 putative novel miRNAs. Quantitative analyses predicted the differential expression of ten miRNAs between 91-R and 91-C, and, based on Gene Ontology and pathway analysis, these ten miRNAs putatively target transcripts encoding proteins involved in detoxification mechanisms. RT-qPCR validated an inverse correlation between levels of differentially-expressed miRNAs and their putatively targeted transcripts, which implies a role of these miRNAs in the differential regulation of detoxification pathways in 91-R compared to 91-C. This study provides evidence associating the differential expression of miRNAs in response to multigenerational DDT selection in Drosophila melanogaster and provides important clues for understanding the possible roles of miRNAs in mediating insecticide resistance traits.

miRNAome expression profiles in the gonads of adult Melopsittacus undulatus

The budgerigar (Melopsittacus undulatus) is one of the most widely studied parrot species, serving as an excellent animal model for behavior and neuroscience research. Until recently, it was unknown how sexual differences in the behavior, physiology, and development of organisms are regulated by differential gene expression. MicroRNAs (miRNAs) are endogenous short non-coding RNA molecules that can post-transcriptionally regulate gene expression and play a critical role in gonadal differentiation as well as early development of animals. However, very little is known about the role gonadal miRNAs play in the early development of birds. Research on the sex-biased expression of miRNAs in avian gonads are limited, and little is known about M. undulatus. In the current study, we sequenced two small non-coding RNA libraries made from the gonads of adult male and female budgerigars using Illumina paired-end sequencing technology. We obtained 254 known and 141 novel miRNAs, and randomly validated five miRNAs. Of these, three miRNAs were differentially expressed miRNAs and 18 miRNAs involved in sexual differentiation as determined by functional analysis with GO annotation and KEGG pathway analysis. In conclusion, this work is the first report of sex-biased miRNAs expression in the budgerigar, and provides additional sequences to the avian miRNAome database which will foster further functional genomic research.

Dynamic miRNA-mRNA regulations are essential for maintaining Drosophila immune homeostasis during Micrococcus luteus infection

Pathogen bacteria infections can lead to dynamic changes of microRNA (miRNA) and mRNA expression profiles, which may control synergistically the outcome of immune responses. To reveal the role of dynamic miRNA-mRNA regulation in Drosophila innate immune responses, we have detailedly analyzed the paired miRNA and mRNA expression profiles at three time points during Drosophila adult males with Micrococcus luteus (M. luteus) infection using RNA- and small RNA-seq data. Our results demonstrate that differentially expressed miRNAs and mRNAs represent extensively dynamic changes over three time points during Drosophila with M. luteus infection. The pathway enrichment analysis indicates that differentially expressed genes are involved in diverse signaling pathways, including Toll and Imd as well as orther signaling pathways at three time points during Drosophila with M. luteus infection. Remarkably, the dynamic change of miRNA expression is delayed by compared to mRNA expression change over three time points, implying that the "time" parameter should be considered when the function of miRNA/mRNA is further studied. In particular, the dynamic miRNA-mRNA regulatory networks have shown that miRNAs may synergistically regulate gene expressions of different signaling pathways to promote or inhibit innate immune responses and maintain homeostasis in Drosophila, and some new regulators involved in Drosophila innate immune response have been identified. Our findings strongly suggest that miRNA regulation is a key mechanism involved in fine-tuning cooperatively gene expressions of diverse signaling pathways to maintain innate immune response and homeostasis in Drosophila. Taken together, the present study reveals a novel role of dynamic miRNA-mRNA regulation in immune response to bacteria infection, and provides a new insight into the underlying molecular regulatory mechanism of Drosophila innate immune responses.

MiR-184 expression is regulated by AMPK in pancreatic islets

AMPK is a critical energy sensor and target for widely used antidiabetic drugs. In β cells, elevated glucose concentrations lower AMPK activity, and the ablation of both catalytic subunits [β-cell–specific AMPK double-knockout (βAMPKdKO) mice] impairs insulin secretion in vivo and β-cell identity. MicroRNAs (miRNAs) are small RNAs that silence gene expression that are essential for pancreatic β-cell function and identity and altered in diabetes. Here, we have explored the miRNAs acting downstream of AMPK in mouse and human β cells. We identified 14 down-regulated and 9 up-regulated miRNAs in βAMPKdKO vs. control islets. Gene ontology analysis of targeted transcripts revealed enrichment in pathways important for β-cell function and identity. The most down-regulated miRNA was miR-184 (miR-184-3p), an important regulator of β-cell function and compensatory expansion that is controlled by glucose and reduced in diabetes. We demonstrate that AMPK is a potent regulator and an important mediator of the negative effects of glucose on miR-184 expression. Additionally, we reveal sexual dimorphism in miR-184 expression in mouse and human islets. Collectively, these data demonstrate that glucose-mediated changes in AMPK activity are central for the regulation of miR-184 and other miRNAs in islets and provide a link between energy status and gene expression in β cells.—Martinez-Sanchez, A., Nguyen-Tu, M.-S., Cebola, I., Yavari, A., Marchetti, P., Piemonti, L., de Koning, E., Shapiro, A. M. J., Johnson, P., Sakamoto, K., Smith, D. M., Leclerc, I., Ashrafian, H., Ferrer, J., Rutter, G. A. MiR-184 expression is regulated by AMPK in pancreatic islets.

Sex-biased microRNA expression in mammals and birds reveals underlying regulatory mechanisms and a role in dosage compensation

Sexual dimorphism depends on sex-biased gene expression, but the contributions of microRNAs (miRNAs) have not been globally assessed. We therefore produced an extensive small RNA sequencing data set to analyze male and female miRNA expression profiles in mouse, opossum, and chicken. Our analyses uncovered numerous cases of somatic sex-biased miRNA expression, with the largest proportion found in the mouse heart and liver. Sex-biased expression is explained by miRNA-specific regulation, including sex-biased chromatin accessibility at promoters, rather than piggybacking of intronic miRNAs on sex-biased protein-coding genes. In mouse, but not opossum and chicken, sex bias is coordinated across tissues such that autosomal testis-biased miRNAs tend to be somatically male-biased, whereas autosomal ovary-biased miRNAs are female-biased, possibly due to broad hormonal control. In chicken, which has a Z/W sex chromosome system, expression output of genes on the Z Chromosome is expected to be male-biased, since there is no global dosage compensation mechanism that restores expression in ZW females after almost all genes on the W Chromosome decayed. Nevertheless, we found that the dominant liver miRNA, miR-122-5p, is Z-linked but expressed in an unbiased manner, due to the unusual retention of a W-linked copy. Another Z-linked miRNA, the male-biased miR-2954-3p, shows conserved preference for dosage-sensitive genes on the Z Chromosome, based on computational and experimental data from chicken and zebra finch, and acts to equalize male-to-female expression ratios of its targets. Unexpectedly, our findings thus establish miRNA regulation as a novel gene-specific dosage compensation mechanism.

Redundant and incoherent regulations of multiple phenotypes suggest microRNAs' role in stability control

Each microRNA (miRNA) represses a web of target genes and, through them, controls multiple phenotypes. The difficulties inherent in such controls cast doubt on how effective miRNAs are in driving phenotypic changes. A "simple regulation" model posits "one target-one phenotype" control under which most targeting is nonfunctional. In an alternative "coordinate regulation" model, multiple targets are assumed to control the same phenotypes coherently, and most targeting is functional. Both models have some empirical support but pose different conceptual challenges. Here, we concurrently analyze multiple targets and phenotypes associated with the miRNA-310 family (miR310s) of Drosophila Phenotypic rescue in the mir310s knockout background is achieved by promoter-directed RNA interference that restores wild-type expression. For one phenotype (eggshell morphology), we observed redundant regulation, hence rejecting "simple regulation" in favor of the "coordinate regulation" model. For other phenotypes (egg-hatching and male fertility), however, one gene shows full rescue, but three other rescues aggravate the phenotype. Overall, phenotypic controls by miR310s do not support either model. Like a thermostat that controls both heating and cooling elements to regulate temperature, redundancy and incoherence in regulation generally suggest some capacity in stability control. Our results therefore support the published view that miRNAs play a role in the canalization of transcriptome and, hence, phenotypes.

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.

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.

Identification of male-biased microRNA-107 as a direct regulator for nuclear receptor subfamily 5 group A member 1 based on sexually dimorphic microRNA expression profiling from chicken embryonic gonads

Several studies indicate that sexual dimorphic microRNAs (miRNAs) in chicken gonads are likely to have important roles in sexual development, but a more global understanding of the roles of miRNAs in sexual differentiation is still needed. To this end, we performed miRNA expression profiling in chicken gonads at embryonic day 5.5 (E5.5). Among the sex-biased miRNAs validated by qRT-PCR, twelve male-biased and six female-biased miRNAs were consistent with the sequencing results. Bioinformatics analysis revealed that some sex-biased miRNAs were potentially involved in gonadal development. Further functional analysis found that over-expression of miR-107 directly inhibited nuclear receptor subfamily 5 group A member 1 (NR5a1), and its downstream cytochrome P450 family 19 subfamily A, polypeptide 1 (CYP19A1). However, anti-Mullerian hormone (AMH) was not directly or indirectly regulated by miR-107. Overall results indicate that miR-107 may specifically mediate avian ovary-development by post-transcriptional regulation of NR5a1 and CYP19A1 in estrogen signaling pathways.

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.

Characterization of the genome and transcriptome of the blue tit Cyanistes caeruleus: polymorphisms, sex-biased expression and selection signals

Decoding genomic sequences and determining their variation within populations has potential to reveal adaptive processes and unravel the genetic basis of ecologically relevant trait variation within a species. The blue tit Cyanistes caeruleus--a long-time ecological model species--has been used to investigate fitness consequences of variation in mating and reproductive behaviour. However, very little is known about the underlying genetic changes due to natural and sexual selection in the genome of this songbird. As a step to bridge this gap, we assembled the first draft genome of a single blue tit, mapped the transcriptome of five females and five males to this reference, identified genomewide variants and performed sex-differential expression analysis in the gonads, brain and other tissues. In the gonads, we found a high number of sex-biased genes, and of those, a similar proportion were sex-limited (genes only expressed in one sex) in males and females. However, in the brain, the proportion of female-limited genes within the female-biased gene category (82%) was substantially higher than the proportion of male-limited genes within the male-biased category (6%). This suggests a predominant on-off switching mechanism for the female-limited genes. In addition, most male-biased genes were located on the Z-chromosome, indicating incomplete dosage compensation for the male-biased genes. We called more than 500,000 SNPs from the RNA-seq data. Heterozygote detection in the single reference individual was highly congruent between DNA-seq and RNA-seq calling. Using information from these polymorphisms, we identified potential selection signals in the genome. We list candidate genes which can be used for further sequencing and detailed selection studies, including genes potentially related to meiotic drive evolution. A public genome browser of the blue tit with the described information is available at http://public-genomes-ngs.molgen.mpg.de.

microRNAs and the evolution of complex multicellularity: identification of a large, diverse complement of microRNAs in the brown alga Ectocarpus

There is currently convincing evidence that microRNAs have evolved independently in at least six different eukaryotic lineages: animals, land plants, chlorophyte green algae, demosponges, slime molds and brown algae. MicroRNAs from different lineages are not homologous but some structural features are strongly conserved across the eukaryotic tree allowing the application of stringent criteria to identify novel microRNA loci. A large set of 63 microRNA families was identified in the brown alga Ectocarpus based on mapping of RNA-seq data and nine microRNAs were confirmed by northern blotting. The Ectocarpus microRNAs are highly diverse at the sequence level with few multi-gene families, and do not tend to occur in clusters but exhibit some highly conserved structural features such as the presence of a uracil at the first residue. No homologues of Ectocarpus microRNAs were found in other stramenopile genomes indicating that they emerged late in stramenopile evolution and are perhaps specific to the brown algae. The large number of microRNA loci in Ectocarpus is consistent with the developmental complexity of many brown algal species and supports a proposed link between the emergence and expansion of microRNA regulatory systems and the evolution of complex multicellularity.

Dynamic evolution and biogenesis of small RNAs during sex reversal

Understanding origin, evolution and functions of small RNA (sRNA) genes has been a great challenge in the past decade. Molecular mechanisms underlying sexual reversal in vertebrates, particularly sRNAs involved in this process, are largely unknown. By deep-sequencing of small RNA transcriptomes in combination with genomic analysis, we identified a large amount of piRNAs and miRNAs including over 1,000 novel miRNAs, which were differentially expressed during gonad reversal from ovary to testis via ovotesis. Biogenesis and expressions of miRNAs were dynamically changed during the reversal. Notably, phylogenetic analysis revealed dynamic expansions of miRNAs in vertebrates and an evolutionary trajectory of conserved miR-17-92 cluster in the Eukarya. We showed that the miR-17-92 cluster in vertebrates was generated through multiple duplications from ancestor miR-92 in invertebrates Tetranychus urticae and Daphnia pulex from the Chelicerata around 580 Mya. Moreover, we identified the sexual regulator Dmrt1 as a direct target of the members miR-19a and -19b in the cluster. These data suggested dynamic biogenesis and expressions of small RNAs during sex reversal and revealed multiple expansions and evolutionary trajectory of miRNAs from invertebrates to vertebrates, which implicate small RNAs in sexual reversal and provide new insight into evolutionary and molecular mechanisms underlying sexual reversal.

MicroRNAs of two medically important mosquito species: Aedes aegypti and Anopheles stephensi

MicroRNAs (miRNAs) are endogenous, single-stranded small RNAs that have important regulatory functions at the post-transcriptional level. In the present study, we characterize miRNAs in two divergent mosquito species, Aedes aegypti and Anopheles stephensi, through deep sequencing of small RNAs spanning all developmental stages. We discovered eight novel miRNAs in Ae. aegypti and 20 novel miRNAs in An. stephensi, which enabled the first systematic analysis of miRNA evolution in mosquitos. We traced the phylogenetic history of all miRNAs in both species and report a rate of 0.055-0.13 miRNA net gain per million years. Most novel miRNAs originate de novo. Duplications that produced miRNA clusters and families are more common in Ae. aegypti than in An. stephensi. We also identified arm-switch as a source of new miRNAs. Expression profile analysis identified mosquito-specific miRNAs that showed strong stage-specific expression in one or both lineages. For example, the aae-miR-2941/2946 family represents the most abundant maternally deposited and zygotically transcribed miRNAs in Ae. aegypti. miR-2943 is a highly expressed zygotic miRNA in both Ae. aegypti and An. stephensi. Such information provides the basis from which to study the function of these miRNAs in biology common to all mosquitos or unique to one particular lineage.

The germline of the malaria mosquito produces abundant miRNAs, endo-siRNAs, piRNAs and 29-nt small RNAs

BACKGROUND

Small RNAs include different classes essential for endogenous gene regulation and cellular defence against genomic parasites. However, a comprehensive analysis of the small RNA pathways in the germline of the mosquito Anopheles gambiae has never been performed despite their potential relevance to reproductive capacity in this malaria vector.

RESULTS

We performed small RNA deep sequencing during larval and adult gonadogenesis and find that they predominantly express four classes of regulatory small RNAs. We identified 45 novel miRNA precursors some of which were sex-biased and gonad-enriched , nearly doubling the number of previously known miRNA loci. We also determine multiple genomic clusters of 24-30 nt Piwi-interacting RNAs (piRNAs) that map to transposable elements (TEs) and 3'UTR of protein coding genes. Unusually, many TEs and the 3'UTR of some endogenous genes produce an abundant peak of 29-nt small RNAs with piRNA-like characteristics. Moreover, both sense and antisense piRNAs from TEs in both Anopheles gambiae and Drosophila melanogaster reveal novel features of piRNA sequence bias. We also discovered endogenous small interfering RNAs (endo-siRNAs) that map to overlapping transcripts and TEs.

CONCLUSIONS

This is the first description of the germline miRNome in a mosquito species and should prove a valuable resource for understanding gene regulation that underlies gametogenesis and reproductive capacity. We also provide the first evidence of a piRNA pathway that is active against transposons in the germline and our findings suggest novel piRNA sequence bias. The contribution of small RNA pathways to germline TE regulation and genome defence in general is an important finding for approaches aimed at manipulating mosquito populations through the use of selfish genetic elements.