Search

Abstract Background Adenosine-to-inosine (A-to-I) RNA editing is a co-/post-transcriptional modification introducing A-to-G variations in RNAs. There is extensive discussion on whether the flexibility of RNA editing exerts a proteomic diversification role, or it just acts like hardwired mutations to correct the genomic allele. Eusocial insects evolved the ability to generate phenotypically differentiated individuals with the same genome, indicating the involvement of epigenetic/transcriptomic regulation. Methods We obtained the genomes of 104 Hymenoptera insects and the transcriptomes of representative species. Comparative genomic analysis was performed to parse the evolutionary trajectory of a regulatory Ile > Met auto-recoding site in Adar gene. Results At genome level, the pre-editing Ile codon is conserved across a node containing all eusocial hymenopterans. At RNA level, the editing events are confirmed in representative species and shows considerable condition-specificity. Compared to random expectation, the editable Ile codon avoids genomic substitutions to Met or to uneditable Ile codons, but does not avoid mutations to other unrelated amino acids. Conclusions The flexibility of Adar auto-recoding site in Hymenoptera is selectively maintained, supporting the flexible RNA editing hypothesis. We proposed a new angle to view the adaptation of RNA editing, providing another layer to explain the great phenotypical plasticity of eusocial insects.

Adar-mediated adenosine-to-inosine (A-to-I) mRNA editing is a conserved mechanism that exerts diverse regulatory functions during the development, evolution, and adaptation of metazoans. The accurate detection of RNA editing sites helps us understand their biological significance. In this work, with an improved genome assembly of honeybee (Apis mellifera), we used a new orthology-based methodology to complement the traditional pipeline of (de novo) RNA editing detection. Compared to the outcome of traditional pipeline, we retrieved many novel editing sites in CDS that are deeply conserved between honeybee and other distantly related insects. The newly retrieved sites were missed by the traditional de novo identification due to the stringent criteria for controlling false-positive rate. Caste-specific editing sites are identified, including an Ile>Met auto-recoding site in Adar. This recoding was even conserved between honeybee and bumblebee, suggesting its putative regulatory role in shaping the phenotypic plasticity of eusocial Hymenoptera. In summary, we proposed a complementary approach to the traditional pipeline and retrieved several previously unnoticed CDS editing sites. From both technical and biological aspects, our works facilitate future researches on finding the functional editing sites and advance our understanding on the connection between RNA editing and the great phenotypic diversity of organisms.

Summary: Genetic basis underlying the biodiversity and phenotypic plasticity are fascinating questions in evolutionary biology. Such molecular diversity can be achieved at multi-omics levels. Here, we sequenced the first chromosome-level genome of assassin bug Rhynocoris fuscipes, a polyphagous generalist predator for biological control of agroecosystems. Compared to non-predatory true bugs Apolygus lucorum and Riptortus pedestris, the R. fuscipes-specific genes were enriched in diet-related genes (e.g., serine proteinase, cytochrome P450) which had higher expression level and more exons than non-diet genes. Extensive A-to-I RNA editing was identified in all three species and showed enrichment in genes associated with diet in R. fuscipes, diversifying the transcriptome. An extended analysis between five predaceous and 27 phytophagous hemipteran species revealed an expansion of diet-related genes in R. fuscipes. Our findings bridge the gap between genotype and phenotype, and also advance our understanding on genetic and epigenetic bases governing the diet shifts in ture bugs.

Abstract Aphidius gifuensis is a parasitoid wasp and primary endoparasitoid enemy of the peach potato aphid, Myzus persicae. Artificially reared, captive wasps of this species have been extensively and effectively used to control populations of aphids and limit crop loss. However, the consequences of large-scale releasing of captive A. gifuensis, such as genetic erosion and reduced fitness in wild populations of this species, remains unclear. Here, we sequence the genomes of 542 A. gifuensis individuals collected across China, including 265 wild and 277 human-intervened samples. Population genetic analyses on wild individuals recovered Yunnan populations as the ancestral group with the most complex genetic structure. We also find genetic signature of environmental adaptation during the dispersal of wild populations from Yunnan to other regions. While comparative genomic analyses of captive wasps revealed a decrease in genetic diversity during long-term rearing, population genomic analyses revealed signatures of natural selection by several biotic (host plants) or abiotic (climate) factors, which support maintenance of the gene pool of wild populations in spite of the introduction of captive wasps. Therefore, the impact of large-scale release is reduced. Our study suggests that A. gifuensis is a good system for exploring the genetic and evolutionary effects of mass rearing and release on species commonly used as biocontrol agents.

Adenosine-to-inosine (A-to-I) RNA editing recodes the genome and confers flexibility for the organisms to adapt to the environment. It is believed that RNA recoding sites are well suited for facilitating adaptive evolution by increasing the proteomic diversity in a temporal-spatial manner. The function and essentiality of a few conserved recoding sites are recognized. However, the experimentally discovered functional sites only make up a small corner of the total sites, and there is still the need to expand the repertoire of such functional sites with bioinformatic approaches. In this study, we define a new category of RNA editing sites termed ‘conserved editing with non-conserved recoding’ and systematically identify such sites in Drosophila editomes, figuring out their selection pressure and signals of adaptation at inter-species and intra-species levels. Surprisingly, conserved editing sites with non-conserved recoding are not suppressed and are even slightly overrepresented in Drosophila. DNA mutations leading to such cases are also favoured during evolution, suggesting that the function of those recoding events in different species might be diverged, specialized, and maintained. Finally, structural prediction suggests that such recoding in potassium channel Shab might increase ion permeability and compensate the effect of low temperature. In conclusion, conserved editing with non-conserved recoding might be functional as well. Our study provides novel aspects in considering the adaptive evolution of RNA editing sites and meanwhile expands the candidates of functional recoding sites for future validation.

ABSTRACTAlthough A-to-I RNA editing leads to similar effects to A-to-G DNA mutation, nonsynonymous RNA editing (recoding) is believed to confer its adaptiveness by ‘epigenetically’ regulating proteomic diversity in a temporospatial manner, avoiding the pleiotropic effect of genomic mutations. Recent discoveries on the evolutionary trajectory of Ser>Gly auto-editing site in insect Adar gene demonstrated a selective advantage to having an editable codon compared to uneditable ones. However, apart from pure observations, quantitative approaches for justifying the adaptiveness of individual RNA editing sites are still lacking. We performed a comparative genomic analysis on 113 Diptera species, focusing on the Adar Ser>Gly auto-recoding site in Drosophila. We only found one species having a derived Gly at the corresponding site, and this occurrence was significantly lower than genome-wide random expectation. This suggests that the Adar Ser>Gly site is unlikely to be genomically replaced with G during evolution, and thus indicating the advantage of editable status over hardwired genomic alleles. Similar trends were observed for the conserved Ile>Met recoding in gene Syt1. In the light of evolution, we established a comparative genomic approach for quantitatively justifying the adaptiveness of individual editing sites. Priority should be given to such adaptive editing sites in future functional studies.

ABSTRACTAdar-mediated adenosine-to-inosine (A-to-I) RNA editing mainly occurs in nucleus and diversifies the transcriptome in a flexible manner. It has been a challenging task to identify beneficial editing sites from the sea of total editing events. The functional Ser>Gly auto-recoding site in insect Adar gene has uneditable Ser codons in ancestral nodes, indicating the selective advantage to having an editable status. Here, we extended this case study to more metazoan species, and also looked for all Drosophila recoding events with potential uneditable synonymous codons. Interestingly, in D. melanogaster, the abundant nonsynonymous editing is enriched in the codons that have uneditable counterparts, but the Adar Ser>Gly case suggests that the editable orthologous codons in other species are not necessarily edited. The use of editable versus ancestral uneditable codon is a smart way to infer the selective advantage of RNA editing, and priority might be given to these editing sites for functional studies due to the feasibility to construct an uneditable allele. Our study proposes an idea to narrow down the candidates of beneficial recoding sites. Meanwhile, we stress that the matched transcriptomes are needed to verify the conservation of editing events during evolution.

Adenosine-to-inosine (A-to-I) RNA editing, mediated by metazoan ADAR enzymes, is a prevalent post-transcriptional modification that diversifies the proteome and promotes adaptive evolution of organisms. The Drosophila Adar gene has an auto-recoding site (termed S>G site) that forms a negative-feedback loop and stabilizes the global editing activity. However, the evolutionary trajectory of Adar S>G site in many other insects remains largely unknown, preventing us from a deeper understanding on the significance of this auto-editing mechanism. In this study, we retrieved the well-annotated genomes of 375 arthropod species including the five major insect orders (Lepidoptera, Diptera, Coleoptera, Hymenoptera and Hemiptera) and several outgroup species. We performed comparative genomic analysis on the Adar auto-recoding S>G site. We found that the ancestral state of insect S>G site was an uneditable serine codon (unSer) and that this state was largely maintained in Hymenoptera. The editable serine codon (edSer) appeared in the common ancestor of Lepidoptera, Diptera and Coleoptera and was almost fixed in the three orders. Interestingly, Hemiptera species possessed comparable numbers of unSer and edSer codons, and a few ‘intermediate codons’, demonstrating a multi-step evolutionary trace from unSer-to-edSer with non-synchronized mutations at three codon positions. We argue that the evolution of Adar S>G site is the best genomic evidence supporting the ‘proteomic diversifying hypothesis’ of RNA editing. Our work deepens our understanding on the evolutionary significance of Adar auto-recoding site which stabilizes the global editing activity and controls transcriptomic diversity.

Abstract Background C-to-U RNA editing in plants is believed to confer its evolutionary adaptiveness by reversing unfavorable DNA mutations. This “restorative hypothesis” has not yet been tested genome-wide. In contrast, A-to-I RNA editing in insects like Drosophila and honeybee is already known to benefit the host by increasing proteomic diversity in a spatial-temporal manner (namely “diversifying hypothesis”). Methods We profiled the RNA editomes of multiple tissues of Arabidopsis thaliana, Drosophila melanogaster, and Apis melifera. We unprecedentedly defined the haplotype diversity (HD) of RNA molecules based on nonsynonymous editing events (recoding sites). Results Signals of adaptation is confirmed in Arabidopsis by observing higher frequencies and levels at nonsynonymous editing sites over synonymous sites. Compared to A-to-I recoding sites in Drosophila, the C-to-U recoding sites in Arabidopsis show significantly lower HD, presumably due to the stronger linkage between C-to-U events. Conclusions C-to-U RNA editing in Arabidopsis is adaptive but it is not designed for diversifying the proteome like A-to-I editing in Drosophila. Instead, C-to-U recoding sites resemble DNA mutations. Our observation supports the restorative hypothesis of plant C-to-U editing which claims that editing is used for fixing unfavorable genomic sequences.

Abstract Bean flower thrips Megalurothrips usitatus is a staple pest of cowpea and other legumes and causes dramatic economic losses. Its small size allows for easy concealment, and large reproductive capacity easily leads to infestations. Despite the importance of a genome in developing novel management strategies, genetic studies on M. usitatus remain limited. Thus, we generated a chromosome-level M. usitatus genome using a combination of PacBio long read and Hi-C technologies. The assembled genome was 238.14 Mb with a scaffold N50 of 13.85 Mb. The final genome was anchored into 16 pseudo-chromosomes containing 14,000 genes, of which 91.74% were functionally annotated. Comparative genomic analyses revealed that expanded gene families were enriched in fatty acid metabolism and detoxification metabolism (ABC transporters), and contracted gene families were strongly associated with chitin-based cuticle development and sensory perception of taste. In conclusion, this high-quality genome provides an invaluable resource for us to understand the thrips’ ecology and genetics, contributing to pest management.

Adenosine-to-inosine (A-to-I) RNA editing is the most prevalent RNA modification in the nervous systems of metazoans. To study the biological significance of RNA editing, we first have to accurately identify these editing events from the transcriptome. The genome-wide identification of RNA editing sites remains a challenging task. In this review, we will first introduce the occurrence, regulation, and importance of A-to-I RNA editing and then describe the established bioinformatic procedures and difficulties in the accurate identification of these sit esespecially in small sized non-model insects. In brief, (1) to obtain an accurate profile of RNA editing sites, a transcriptome coupled with the DNA resequencing of a matched sample is favorable; (2) the single-cell sequencing technique is ready to be applied to RNA editing studies, but there are a few limitations to overcome; (3) during mapping and variant calling steps, various issues, like mapping and base quality, soft-clipping, and the positions of mismatches on reads, should be carefully considered; (4) Sanger sequencing of both RNA and the matched DNA is the best verification of RNA editing sites, but other auxiliary evidence, like the nonsynonymous-to-synonymous ratio or the linkage information, is also helpful for judging the reliability of editing sites. We have systematically reviewed the understanding of the biological significance of RNA editing and summarized the methodology for identifying such editing events. We also raised several promising aspects and challenges in this field. With insightful perspectives on both scientific and technical issues, our review will benefit the researchers in the broader RNA editing community.

Adenosine-to-inosine (A-to-I) RNA editing leads to a similar effect to A-to-G mutations. RNA editing provides a temporo-spatial flexibility for organisms. Nonsynonymous (Nonsyn) RNA editing in insects is over-represented compared with synonymous (Syn) editing, suggesting adaptive signals of positive selection on Nonsyn editing during evolution. We utilized the brain RNA editome of Drosophila melanogaster to systematically study the LD (r2) between editing sites and infer its impact on the adaptive signals of RNA editing. Pairs of editing sites (PESs) were identified from the transcriptome. For CDS PESs of two consecutive editing sites, their occurrence was significantly biased to type-3 PES (Syn-Nonsyn). The haplotype frequency of type-3 PES exhibited a significantly higher abundance of AG than GA, indicating that the rear Nonsyn site is the driver that promotes the editing of the front Syn site (passenger). The exclusion of passenger Syn sites dramatically amplifies the adaptive signal of Nonsyn RNA editing. Our study for the first time quantitatively demonstrates that the linkage between RNA editing events comes from hitchhiking effects and leads to the underestimation of adaptive signals for Nonsyn editing. Our work provides novel insights for studying the evolutionary significance of RNA editing events.

Metazoan adenosine-to-inosine (A-to-I) RNA editing is a highly conserved mechanism that diversifies the transcriptome by post-transcriptionally converting adenosine to inosine. Millions of editing sites have been identified in different species and, based on abnormal editing observed in various disorders, it is intuitive to conclude that RNA editing is both functional and adaptive. In this review, we propose the following major points: (1) “Function/functional” only represents a molecular/phenotypic consequence and is not necessarily connected to “adaptation/adaptive”; (2) Adaptive editing should be judged in the light of evolution and emphasize advantages of temporal-spatial flexibility; (3) Adaptive editing could, in theory, be extended from nonsynonymous sites to all potentially functional sites. This review seeks to conceptually bridge the gap between molecular biology and evolutionary biology and provide a more objective understanding on the biological functions and evolutionary significance of RNA editing.

Coridius chinensis (Hemiptera: Dinidoridae) is a medicinal insect. Its mitochondrial gene arrangement is consistent with that of Drosophila melanogaster and Erthesina fullo, the two insects with well-studied mitochondrial transcription. To investigate whether the structural consistency of mitochondrial genes leads to similarities in transcription and post-transcriptional processing, we improved the gene annotation and constructed a quantitative transcription map for the C. chinensis mitochondrial genome (mitogenome) using full-length transcriptome sequencing. The size of this mitogenome was 16,214 bp and the proposed model of mitochondrial transcription was similar to that of Drosophila. Both strands were nearly entirely transcribed except for the antisense genes downstream of trnS2 on N strand. The expression of cytochrome c subunit genes is higher than that of NADH-dehydrogenase subunit genes. The post-transcriptional cleavage process followed the “tRNA punctuation” model, and both the “reverse cleavage” model in Drosophila and “forward cleavage” model in E. fullo were found in C. chinensis. In addition, we found that long non-coding RNAs from the control region contained tandem repeats. Polyadenylation was performed after CCA triplet at the 3′ end of tRNA. The isoform diversity of lrRNA was identified. Our study sheds light on the transcriptional regulation and RNA processing of insect mitogenomes with the putative ancestral gene arrangement.

Abstract Background The piwi-interacting RNAs (piRNAs) are small non-coding RNAs that specifically repress transposable elements (TEs) in the germline of Drosophila. Despite our expanding understanding of TE:piRNA interaction, whether there is an evolutionary arms race between TEs and piRNAs was unclear. Results Here, we studied the population genomics of TEs and piRNAs in the worldwide strains of D. melanogaster. By conducting a correlation analysis between TE contents and the abundance of piRNAs from ovaries of representative strains of D. melanogaster, we find positive correlations between TEs and piRNAs in six TE families. Our simulations further highlight that TE activities and the strength of purifying selection against TEs are important factors shaping the interactions between TEs and piRNAs. Our studies also suggest that the de novo generation of piRNAs is an important mechanism to repress the newly invaded TEs. Conclusions Our results revealed the existence of an evolutionary arms race between the copy numbers of TEs and the abundance of antisense piRNAs at the population level. Although the interactions between TEs and piRNAs are complex and many factors should be considered to impact their interaction dynamics, our results suggest the emergence, repression specificity and strength of piRNAs on TEs should be considered in studying the landscapes of TE insertions in Drosophila. These results deepen our understanding of the interactions between piRNAs and TEs, and also provide novel insights into the nature of genomic conflicts of other forms.

Summary: Social insects exhibit extensive phenotypic diversities among the genetically similar individuals, suggesting a role for the epigenetic regulations beyond the genome level. The ADAR-mediated adenosine-to-inosine (A-to-I) RNA editing, an evolutionarily conserved mechanism, facilitates adaptive evolution by expanding proteomic diversities. Here, we characterize the A-to-I RNA editome of honeybees (Apis mellifera), identifying 407 high-confidence A-to-I editing sites. Editing is most abundant in the heads and shows signatures for positive selection. Editing behavior differs between foragers and nurses, suggesting a role for editing in caste differentiation. Although only five sites are conserved between bees and flies, an unexpectedly large number of genes exhibit editing in both species, albeit at different locations, including the nonsynonymous auto-editing of Adar. This convergent evolution, where the same target genes independently acquire recoding events in distant diverged clades, together with the signals of adaptation observed in honeybees alone, further supports the notion of recoding being adaptive.

Adenosine-to-inosine (A-to-I) editing is hypothesized to facilitate adaptive evolution by expanding proteomic diversity through an epigenetic approach. However, it is challenging to provide evidences to support this hypothesis at the whole editome level. In this study, we systematically characterized 2,114 A-to-I RNA editing sites in female and male brains of D. melanogaster, and nearly half of these sites had events evolutionarily conserved across Drosophila species. We detected strong signatures of positive selection on the nonsynonymous editing sites in Drosophila brains, and the beneficial editing sites were significantly enriched in genes related to chemical and electrical neurotransmission. The signal of adaptation was even more pronounced for the editing sites located in X chromosome or for those commonly observed across Drosophila species. We identified a set of gene candidates (termed "PSEB" genes) that had nonsynonymous editing events favored by natural selection. We presented evidence that editing preferentially increased mutation sequence space of evolutionarily conserved genes, which supported the adaptive evolution hypothesis of editing. We found prevalent nonsynonymous editing sites that were favored by natural selection in female and male adults from five strains of D. melanogaster. We showed that temperature played a more important role than gender effect in shaping the editing levels, although the effect of temperature is relatively weaker compared to that of species effect. We also explored the relevant factors that shape the selective patterns of the global editomes. Altogether we demonstrated that abundant nonsynonymous editing sites in Drosophila brains were adaptive and maintained by natural selection during evolution. Our results shed new light on the evolutionary principles and functional consequences of RNA editing.

Background: Metazoan adenosine-to-inosine (A-to-I) RNA editing resembles A-to-G mutation and increases proteomic diversity in a temporal-spatial manner, allowing organisms adapting to changeable environment. The RNA editomes in many major animal clades remain unexplored, hampering the understanding on the evolution and adaptation of this essential post-transcriptional modification. Methods: We assembled the chromosome-level genome of Coridius chinensis belonging to Hemiptera, the fifth largest insect order where RNA editing has not been studied yet. We generated ten head RNA-Seq libraries with DNA-Seq from the matched individuals. Results: We identified thousands of high-confidence RNA editing sites in C. chinensis. Overrepresentation of nonsynonymous editing was observed, but conserved recoding across different orders was very rare. Under cold stress, the global editing efficiency was down-regulated and the general transcriptional processes were shut down. Nevertheless, we found an interesting site with "conserved editing but non-conserved recoding" in potassium channel Shab which was significantly up-regulated in cold, serving as a candidate functional site in response to temperature stress. Conclusions: RNA editing in C. chinensis largely recodes the proteome. The first RNA editome in Hemiptera indicates independent origin of beneficial recoding during insect evolution, which advances our understanding on the evolution, conservation, and adaptation of RNA editing. [ABSTRACT FROM AUTHOR]