Large-Scale Annotation and Evolution Analysis of MiRNA in Insects

A microRNA is the effector gene of a classic evolutionary hotspot locus

In Lepidoptera (butterflies and moths), the genomic region around the gene cortex is a "hotspot" locus, repeatedly implicated in generating intraspecific melanic wing color polymorphisms across 100 million years of evolution. However, the identity of the effector gene regulating melanic wing color within this locus remains unknown. We show that none of the four candidate protein-coding genes within this locus, including cortex, serve as major effectors. Instead, a microRNA (miRNA), mir-193, serves as the major effector across three deeply diverged lineages of butterflies, and its role is conserved in Drosophila. In Lepidoptera, mir-193 is derived from a gigantic primary long noncoding RNA, ivory, and it functions by directly repressing multiple pigmentation genes. We show that a miRNA can drive repeated instances of adaptive evolution in animals.

The roles of different gene expression regulators in acoustic variation in the intermediate horseshoe bat revealed by long-read and short-read RNA sequencing data

Gene expression changes contribute greatly to phenotypic variations in nature. Studying patterns of regulators of gene expression is important to fully understand the molecular mechanism underlying phenotypic variations. In horseshoe bats, the cochleae are finely tuned to echoes of call frequency. Here, using 2 recently diverged subspecies of the intermediate horseshoe bat (Rhinolophus affinis hainanus and R. a. himalayanus) with great acoustic variations as the system, we aim to explore relative roles of different regulators of gene expression (differential gene expression, alternative splicing (AS) and long non-coding RNAs (lncRNAs)) in phenotypic variation with a combination of Illumina short-read and Nanopore long-read RNA-seq data from the cochlea. Compared to R. a. hainanus, R. a. himalayanus exhibited much more upregulated differentially expressed genes (DEGs) and multiple of them may play important roles in the maintenance and damage repair of auditory hair cells. We identified 411 differentially expressed lncRNAs and their target DEGs upregulated in R. a. himalayanus were also mainly involved in a protective mechanism for auditory hair cells. Using 3 different methods of AS analysis, we identified several candidate alternatively spliced genes (ASGs) that expressed different isoforms which may be associated with acoustic divergence of the 2 subspecies. We observed significantly less overlap than expected between DEGs and ASGs, supporting complementary roles of differential gene expression and AS in generating phenotypic variations. Overall, our study highlights the importance of a combination of short-read and long-read RNA-seq data in examining the regulation of gene expression changes responsible for phenotypic variations.

microRNA maintains nutrient homeostasis in the symbiont-host interaction

Endosymbionts provide essential nutrients for hosts, promoting growth, development, and reproduction. However, the molecular regulation of nutrient transport from endosymbiont to host is not well understood. Here, we used bioinformatic analysis, RNA-Sequencing, luciferase assays, RNA immunoprecipitation, and in situ hybridization to show that a bacteriocyte-distributed MRP4 gene (multidrug resistance-associated protein 4) is negatively regulated by a host (aphid)-specific microRNA (miR-3024). Targeted metabolomics, microbiome analysis, vitamin B6 (VB6) supplements, 3D modeling/molecular docking, in vitro binding assays (voltage clamp recording and microscale thermophoresis), and functional complementation of Escherichia coli were jointly used to show that the miR-3024/MRP4 axis controls endosymbiont (Serratia)-produced VB6 transport to the host. The supplementation of miR-3024 increased the mortality of aphids, but partial rescue was achieved by providing an external source of VB6. The use of miR-3024 as part of a sustainable aphid pest-control strategy was evaluated by safety assessments in nontarget organisms (pollinators, predators, and entomopathogenic fungi) using virus-induced gene silencing assays and the expression of miR-3024 in transgenic tobacco. The supplementation of miR-3024 suppresses MRP4 expression, restricting the number of membrane channels, inhibiting VB6 transport, and ultimately killing the host. Under aphids facing stress conditions, the endosymbiont titer is decreased, and the VB6 production is also down-regulated, while the aphid's autonomous inhibition of miR-3024 enhances the expression of MRP4 and then increases the VB6 transport which finally ensures the VB6 homeostasis. The results confirm that miR-3024 regulates nutrient transport in the endosymbiont-host system and is a suitable target for sustainable pest control.

Species-Specific miRNAs Contribute to the Divergence between Deciduous and Evergreen Species in Ilex

MicroRNAs (miRNAs) are pivotal regulators of gene expression, playing crucial roles in plant developmental processes and environmental responses. However, the function of miRNAs in influencing deciduous traits has been little explored. Here, we utilized sRNA-seq on two deciduous species, Ilex polyneura (Hand.-Mazz.) S. Y. Hu and Ilex asprella Champ. ex Benth., along with an evergreen species, Ilex latifolia Thunb., to identify and annotate miRNAs within these species. Our analysis revealed 162 species-specific miRNAs (termed SS-miRNAs) from 120 families, underscoring the fundamental roles and potential influence of SS-miRNAs on plant phenotypic diversity and adaptation. Notably, three SS-miRNAs in I. latifolia were found to target crucial genes within the abscission signaling pathway. Analysis of cis-regulatory elements suggested a novel regulatory relationship that may contribute to the evergreen phenotype of I. latifolia by modulating the abscission process in a light-independent manner. These findings propose a potential mechanism by which SS-miRNAs can influence the conserved abscission pathway, contributing to the phenotypic divergence between deciduous and evergreen species within the genus Ilex.

A micro-RNA is the effector gene of a classic evolutionary hotspot locus

In Lepidoptera (butterflies and moths), the genomic region around the gene cortex is a 'hotspot' locus, repeatedly used to generate intraspecific melanic wing color polymorphisms across 100-million-years of evolution. However, the identity of the effector gene regulating melanic wing color within this locus remains unknown. Here, we show that none of the four candidate protein-coding genes within this locus, including cortex, serve as major effectors. Instead, a micro-RNA (miRNA), mir-193, serves as the major effector across three deeply diverged lineages of butterflies, and its function is conserved in Drosophila. In Lepidoptera, mir-193 is derived from a gigantic long non-coding RNA, ivory, and it functions by directly repressing multiple pigmentation genes. We show that a miRNA can drive repeated instances of adaptive evolution in animals.

Mitochondrially mediated RNA interference, a retrograde signaling system affecting nuclear gene expression

Several functional classes of short noncoding RNAs are involved in manifold regulatory processes in eukaryotes, including, among the best characterized, miRNAs. One of the most intriguing regulatory networks in the eukaryotic cell is the mito-nuclear crosstalk: recently, miRNA-like elements of mitochondrial origin, called smithRNAs, were detected in a bivalve species, Ruditapes philippinarum. These RNA molecules originate in the organelle but were shown in vivo to regulate nuclear genes. Since miRNA genes evolve easily de novo with respect to protein-coding genes, in the present work we estimate the probability with which a newly arisen smithRNA finds a suitable target in the nuclear transcriptome. Simulations with transcriptomes of 12 bivalve species suggest that this probability is high and not species specific: one in a hundred million (1 × 10-8) if five mismatches between the smithRNA and the 3' mRNA are allowed, yet many more are allowed in animals. We propose that novel smithRNAs may easily evolve as exaptation of the pre-existing mitochondrial RNAs. In turn, the ability of evolving novel smithRNAs may have played a pivotal role in mito-nuclear interactions during animal evolution, including the intriguing possibility of acting as speciation trigger.

Rapid adaptation of cellular metabolic rate to the MicroRNA complements of mammals and its relevance to the evolution of endothermy

The metabolic efficiency of mammalian cells depends on the attenuation of intrinsic translation noise by microRNAs. We devised a metric of cellular metabolic rate (cMR), rMR/Mexp optimally fit to the number of microRNA families (mirFam), that is robust to variation in mass and sensitive to body temperature (Tb), consistent with the heat dissipation limit theory of Speakman and Król (2010). Using mirFam as predictor, an Ornstein-Uhlenbeck process of stabilizing selection, with an adaptive shift at the divergence of Boreoeutheria, accounted for 95% of the variation in cMR across mammals. Branchwise rates of evolution of cMR, mirFam and Tb concurrently increased 6- to 7-fold at the divergence of Boreoeutheria, independent of mass. Cellular MR variation across placental mammals was also predicted by the sum of model conserved microRNA-target interactions, revealing an unexpected degree of integration of the microRNA-target apparatus into the energy economy of the mammalian cell.

microRNAs in Syrista parreyssi (Hymenoptera) and Lepisma saccharina (Zygentoma) possibly involved in the mitochondrial function

Mitochondria are essential organelles for maintaining vital cellular functions, and microRNAs (miRNAs) regulate gene expression posttranscriptionally. miRNAs exhibit tissue and time-specific patterns in mitochondria and specifically mitochondrial miRNAs (mitomiRs) can regulate the mRNA expression both originating from mitochondrial and nuclear transcription which affect mitochondrial metabolic activity and cell homeostasis. In this study, miRNAs of two insect species, Syrista parreyssi (Hymenoptera) and Lepisma saccharina (Zygentoma), were investigated for the first time. The known and possible novel miRNAs were predicted and characterized and their potential effects on mitochondrial transcription were investigated in these insect species using deep sequencing. The previously reported mitomiRs were also investigated and housekeeping miRNAs were characterized. miRNAs that are involved in mitochondrial processes such as apoptosis and signaling and that affect genes encoding the subunits of OXPHOS complexes have been identified in each species. Here, 81 and 161 novel mature miRNA candidates were bioinformatically predicted and 9 and 24 of those were aligned with reference mitogenomes of S. parreyssi and L. saccharina, respectively. As a result of RNAHybrid analysis, 51 and 69 potential targets of miRNAs were found in the mitogenome of S. parreyssi and L. saccharina, respectively. cox1 gene was the most targeted gene and cytB, rrnS, and rrnL genes were highly targeted in both of the species by novel miRNAs, hypothetically. We speculate that these novel miRNAs, originating from or targeting mitochondria, influence on rRNA genes or positively selected mitochondrial protein-coding genes. These findings may provide a new perspective in evaluating miRNAs for maintaining mitochondrial function and transcription.

Identification of Species-Specific MicroRNAs Provides Insights into Dynamic Evolution of MicroRNAs in Plants

MicroRNAs (miRNAs) are an important class of regulatory small RNAs that program gene expression, mainly at the post-transcriptional level. Although sporadic examples of species-specific miRNAs (termed SS-miRNAs) have been reported, a genome-scale study across a variety of distant species has not been assessed. Here, by comprehensively analyzing miRNAs in 81 plant species phylogenetically ranging from chlorophytes to angiosperms, we identified 8048 species-specific miRNAs from 5499 families, representing over 61.2% of the miRNA families in the examined species. An analysis of the conservation from different taxonomic levels supported the high turnover rate of SS-miRNAs, even over short evolutionary distances. A comparison of the intrinsic features between SS-miRNAs and NSS-miRNAs (non-species-specific miRNAs) indicated that the AU content of mature miRNAs was the most striking difference. Our data further illustrated a significant bias of the genomic coordinates towards SS-miRNAs lying close to or within genes. By analyzing the 125,267 putative target genes for the 7966 miRNAs, we found the preferentially regulated functions of SS-miRNAs related to diverse metabolic processes. Collectively, these findings underscore the dynamic evolution of miRNAs in the species-specific lineages.

The microRNAs in the antennae of Apolygus lucorum (Hemiptera: Miridae): Expression properties and targets prediction

MicroRNAs (miRNAs) regulate gene expression and contribute to numerous physiological processes. However, little is known about the functions of miRNAs in insect chemosensation. In this study, nine small RNA libraries were constructed and sequenced from the antennae of nymphs, adult males, and adult females of Apolygus lucorum. In total, 399 (275 known and 124 novel) miRNAs were identified. miR-7-5p_1 was the most abundant miRNA. Altogether, 69,708 target genes related to biogenesis, membrane, and binding activities were predicted. In particular, 15 miRNAs targeted 16 olfactory genes. Comparing the antennae of nymphs and adult males and females, 94 miRNAs were differentially expressed. Alternatively, a subset of differentially expressed miRNAs was verified by qPCR, supporting the reliability of the sequencing results. This study provides a global miRNA transcriptome for the antennae of A. lucorum and valuable information for further investigations of the functions of miRNAs in the regulation of chemosensation.

A high-quality genome of the dobsonfly Neoneuromus ignobilis reveals molecular convergences in aquatic insects

Neoneuromus ignobilis is an archaic holometabolous aquatic predatory insect. However, a lack of genomic resources hinders the use of whole genome sequencing to explore their genetic basis and molecular mechanisms for adaptive evolution. Here, we provided a high-contiguity, chromosome-level genome assembly of N. ignobilis using high coverage Nanopore and PacBio reads with the Hi-C technique. The final assembly is 480.67 MB in size, containing 12 telomere-ended pseudochromosomes with only 17 gaps. We compared 42 hexapod species genomes including six independent lineages comprising 11 aquatic insects, and found convergent expansions of long wavelength-sensitive and blue-sensitive opsins, thermal stress response TRP channels, and sulfotransferases in aquatic insects, which may be related to their aquatic adaptation. We also detected strong nonrandom signals of convergent amino acid substitutions in aquatic insects. Collectively, our comparative genomic analysis revealed the evidence of molecular convergences in aquatic insects during both gene family evolution and convergent amino acid substitutions.

A transcriptomic atlas underlying developmental plasticity of seasonal forms of Bicyclus anynana butterflies

Organisms residing in regions with alternating seasons often develop different phenotypes, or forms, in each season. These forms are often adaptations to each season and result from an altered developmental response to specific environmental cues such as temperature. While multiple studies have examined form-specific gene expression profiles in a diversity of species, little is known about how environments and developmental transitions, cued by hormone pulses, alter post-transcriptional patterns. In this study, we examine how gene expression, alternative splicing, and miRNA-mediated gene silencing in Bicyclus anynana butterfly hindwing tissue, varies across two rearing temperatures at four developmental timepoints. These timepoints flank two temperature-sensitive periods that coincide with two pulses of the insect hormone 20E. Our results suggest that developmental transitions, coincident with 20E pulses, elicit a greater impact on all these transcriptomic patterns than rearing temperatures per se. More similar transcriptomic patterns are observed pre-20E pulses than those observed post-20E pulses. We also found functionally distinct sets of differentially expressed and differentially spliced genes in the seasonal forms. Furthermore, around 10% of differentially expressed genes are predicted to be direct targets of, and regulated by, differentially expressed miRNAs between the seasonal forms. Many differentially expressed genes, miRNAs, or differentially spliced genes potentially regulate eyespot size plasticity, and we validated the differential splicing pattern of one such gene, daughterless. We present a comprehensive and interactive transcriptomic atlas of the hindwing tissue of both seasonal forms of B. anynana throughout development, a model organism of seasonal plasticity.

Advances in Editing Silkworms (Bombyx mori) Genome by Using the CRISPR-Cas System

CRISPR (clustered regularly interspaced short palindromic repeats)-Cas (CRISPR-associated) represents a powerful genome editing technology that revolutionized in a short period of time numerous natural sciences branches. Therefore, extraordinary progress was made in various fields, such as entomology or biotechnology. Bombyx mori is one of the most important insects, not only for the sericulture industry, but for numerous scientific areas. The silkworms play a key role as a model organism, but also as a bioreactor for the recombinant protein production. Nowadays, the CRISPR-Cas genome editing system is frequently used in order to perform gene analyses, to increase the resistance against certain pathogens or as an imaging tool in B. mori. Here, we provide an overview of various studies that made use of CRISPR-Cas for B. mori genome editing, with a focus on emphasizing the high applicability of this system in entomology and biological sciences.

Identification and functional characterization of olfactory indolergic receptors in Drosophila melanogaster

Indole-sensitive odorant receptors or indolORs belong to a mosquito-specific expansion as ancient as the Culicidae lineage. Brachyceran flies appeared to lack representative members of this group despite the importance of indolics in this important group of dipterans. To explore whether indolORs occur in other brachyceran species, we searched for candidate indolORs in Drosophila melanogaster. Using phylogenetic tools, we show that D. melanogaster OR30a, OR43a, and OR49b form a distinct monophyletic lineage with mosquito indolORs. To explore a potential functional orthology with indolORs, we expressed these three Drosophila ORs in Xenopus laevis oocytes and measured their responses to a panel of indolic compounds. We provide evidence that OR30a, OR43a, and OR49b exhibit high sensitivity to indoles. Along with the recent discovery of indolORs in the housefly Musca domestica, our findings suggest that indolORs are a widespread feature of the peripheral olfactory systems of Diptera.

MirGeneDB 2.1: toward a complete sampling of all major animal phyla

We describe an update of MirGeneDB, the manually curated microRNA gene database. Adhering to uniform and consistent criteria for microRNA annotation and nomenclature, we substantially expanded MirGeneDB with 30 additional species representing previously missing metazoan phyla such as sponges, jellyfish, rotifers and flatworms. MirGeneDB 2.1 now consists of 75 species spanning over ∼800 million years of animal evolution, and contains a total number of 16 670 microRNAs from 1549 families. Over 6000 microRNAs were added in this update using ∼550 datasets with ∼7.5 billion sequencing reads. By adding new phylogenetically important species, especially those relevant for the study of whole genome duplication events, and through updating evolutionary nodes of origin for many families and genes, we were able to substantially refine our nomenclature system. All changes are traceable in the specifically developed MirGeneDB version tracker. The performance of read-pages is improved and microRNA expression matrices for all tissues and species are now also downloadable. Altogether, this update represents a significant step toward a complete sampling of all major metazoan phyla, and a widely needed foundation for comparative microRNA genomics and transcriptomics studies. MirGeneDB 2.1 is part of RNAcentral and Elixir Norway, publicly and freely available at http://www.mirgenedb.org/.