New comparative genomic evidence supporting the proteomic diversification role of A-to-I RNA editing in insects.

Adenosine-to-inosine (A-to-I) RNA editing, resembling A-to-G mutation, confers adaptiveness by increasing proteomic diversity in a temporal-spatial manner. This evolutionary theory named "proteomic diversifying hypothesis" has only partially been tested in very few organisms like Drosophila melanogaster, mainly by observing the positive selection on nonsynonymous editing events. To find additional genome-wide evidences supporting this interesting assumption, we retrieved the genomes of four Drosophila species and collected 20 deep-sequenced transcriptomes of different developmental stages and neuron populations of D. melanogaster. We systematically profiled the RNA editomes in these samples and performed meticulous comparative genomic analyses. Further evidences were found to support the diversifying hypothesis. (1) None of the nonsynonymous editing sites in D. melanogaster had ancestral G-alleles, while the silent editing sites had an unignorable fraction of ancestral G-alleles; (2) Only very few nonsynonymous editing sites in D. melanogaster had corresponding G-alleles derived in the genomes of sibling species, and the fraction of such situation was significantly lower than that of silent editing sites; (3) The few nonsynonymous editing with corresponding G-alleles had significantly more variable editing levels (across samples) than other nonsynonymous editing sites in D. melanogaster. The proteomic diversifying nature of RNA editing in Drosophila excludes the restorative role which favors an ancestral G-allele. The few fixed G-alleles in sibling species might facilitate the adaptation to particular environment and the corresponding nonsynonymous editing in D. melanogaster would introduce stronger advantage of flexible proteomic diversification. With multi-Omics data, our study consolidates the nature of evolutionary significance of A-to-I RNA editing sites in model insects.

The first A-to-I RNA editome of hemipteran species Coridius chinensis reveals overrepresented recoding and prevalent intron editing in early-diverging insects.

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.

Comparative genomic analyses reveal evidence for adaptive A-to-I RNA editing in insect Adar gene.

Although 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.

Natural selection and genetic diversity maintenance in a parasitic wasp during continuous biological control application.

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.

Molecular insights into the adaptive evolution of SARS-CoV-2 spike protein.

The COVID-19 pandemic, caused by the severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2), has substantially damaged the global economy and human health. The spike (S) protein of coronaviruses plays a pivotal role in viral entry by binding to host cell receptors. Additionally, it acts as the primary target for neutralizing antibodies in those infected and is the central focus for currently utilized or researched vaccines. During the virus's adaptation to the human host, the S protein of SARS-CoV-2 has undergone significant evolution. As the COVID-19 pandemic has unfolded, new mutations have arisen and vanished, giving rise to distinctive amino acid profiles within variant of concern strains of SARS-CoV-2. Notably, many of these changes in the S protein have been positively selected, leading to substantial alterations in viral characteristics, such as heightened transmissibility and immune evasion capabilities. This review aims to provide an overview of our current understanding of the structural implications associated with key amino acid changes in the S protein of SARS-CoV-2. These research findings shed light on the intricate and dynamic nature of viral evolution, underscoring the importance of continuous monitoring and analysis of viral genomes. Through these molecular-level investigations, we can attain deeper insights into the virus's adaptive evolution, offering valuable guidance for designing vaccines and developing antiviral drugs to combat the ever-evolving viral threats.

Chromosome-level genome of the poultry shaft louse Menopon gallinae provides insight into the host-switching and adaptive evolution of parasitic lice.

BACKGROUND: Lice (Psocodea: Phthiraptera) are one important group of parasites that infects birds and mammals. It is believed that the ancestor of parasitic lice originated on the ancient avian host, and ancient mammals acquired these parasites via host-switching from birds. Here we present the first chromosome-level genome of Menopon gallinae in Amblycera (earliest diverging lineage of parasitic lice). We explore the transition of louse host-switching from birds to mammals at the genomic level by identifying numerous idiosyncratic genomic variations. RESULTS: The assembled genome is 155 Mb in length, with a contig N50 of 27.42 Mb. Hi-C scaffolding assigned 97% of the bases to 5 chromosomes. The genome of M. gallinae retains a basal insect repertoire of 11,950 protein-coding genes. By comparing the genomes of lice to those of multiple representative insects in other orders, we discovered that gene families of digestion, detoxification, and immunity-related are generally conserved between bird lice and mammal lice, while mammal lice have undergone a significant reduction in genes related to chemosensory systems and temperature. This suggests that mammal lice have lost some of these genes through the adaption to environment and temperatures after host-switching. Furthermore, 7 genes related to hematophagy were positively selected in mammal lice, suggesting their involvement in the hematophagous behavior. CONCLUSIONS: Our high-quality genome of M. gallinae provides a valuable resource for comparative genomic research in Phthiraptera and facilitates further studies on adaptive evolution of host-switching within parasitic lice.

Narrowing down the candidates of beneficial A-to-I RNA editing by comparing the recoding sites with uneditable counterparts.

Adar-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.

Genetic Modification of a Hox Locus Drives Mimetic Color Pattern Variation in a Highly Polymorphic Bumble Bee.

Müllerian mimicry provides natural replicates ideal for exploring mechanisms underlying adaptive phenotypic divergence and convergence, yet the genetic mechanisms underlying mimetic variation remain largely unknown. The current study investigates the genetic basis of mimetic color pattern variation in a highly polymorphic bumble bee, Bombus breviceps (Hymenoptera, Apidae). In South Asia, this species and multiple comimetic species converge onto local Müllerian mimicry patterns by shifting the abdominal setal color from orange to black. Genetic crossing between the orange and black phenotypes suggested the color dimorphism being controlled by a single Mendelian locus, with the orange allele being dominant over black. Genome-wide association suggests that a locus at the intergenic region between 2 abdominal fate-determining Hox genes, abd-A and Abd-B, is associated with the color change. This locus is therefore in the same intergenic region but not the same exact locus as found to drive red black midabdominal variation in a distantly related bumble bee species, Bombus melanopygus. Gene expression analysis and RNA interferences suggest that differential expression of an intergenic long noncoding RNA between abd-A and Abd-B at the onset setal color differentiation may drive the orange black color variation by causing a homeotic shift late in development. Analysis of this same color locus in comimetic species reveals no sequence association with the same color shift, suggesting that mimetic convergence is achieved through distinct genetic routes. Our study establishes Hox regions as genomic hotspots for color pattern evolution in bumble bees and demonstrates how pleiotropic developmental loci can drive adaptive radiations in nature.

Identification and Interpretation of A-to-I RNA Editing Events in Insect Transcriptomes.

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.

The Many Roles of A-to-I RNA Editing in Animals: Functional or Adaptive?

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.

Genome-Wide Analysis on Driver and Passenger RNA Editing Sites Suggests an Underestimation of Adaptive Signals in Insects.

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.

A full repertoire of Hemiptera genomes reveals a multi-step evolutionary trajectory of auto-RNA editing site in insect Adar gene.

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.

Differential adaptive RNA editing signals between insects and plants revealed by a new measurement termed haplotype diversity.

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.

Autorecoding A-to-I RNA editing sites in the Adar gene underwent compensatory gains and losses in major insect clades.

As one of the most prevalent RNA modifications in animals, adenosine-to-inosine (A-to-I) RNA editing facilitates the environmental adaptation of organisms by diversifying the proteome in a temporal-spatial manner. In flies and bees, the editing enzyme Adar has independently gained two different autorecoding sites that form an autofeedback loop, stabilizing the overall editing efficiency. This ensures cellular homeostasis by keeping the normal function of target genes. However, in a broader range of insects, the evolutionary dynamics and significance of this Adar autoregulatory mechanism are unclear. We retrieved the genomes of 377 arthropod species covering the five major insect orders (Hemiptera, Hymenoptera, Coleoptera, Diptera, and Lepidoptera) and aligned the Adar autorecoding sites across all genomes. We found that the two autorecoding sites underwent compensatory gains and losses during the evolution of two orders with the most sequenced species (Diptera and Hymenoptera), and that the two editing sites were mutually exclusive among them: One editable site is significantly linked to another uneditable site. This autorecoding mechanism of Adar could flexibly diversify the proteome and stabilize global editing activity. Many insects independently selected different autorecoding sites to achieve a feedback loop and regulate the global RNA editome, revealing an interesting phenomenon during evolution. Our study reveals the evolutionary force acting on accurate regulation of RNA editing activity in insects and thus deepens our understanding of the functional importance of RNA editing in environmental adaptation and evolution.

Chromosome-level genome assembly of bean flower thrips Megalurothrips usitatus (Thysanoptera: Thripidae).

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.

Chloroplast C-to-U RNA editing in vascular plants is adaptive due to its restorative effect: testing the restorative hypothesis.

The adaptiveness of nonsynonymous RNA editing (recoding) could be conferred by the flexibility of the temporal-spatially controllable proteomic diversity, or by its restorative effect which fixes unfavorable genomic mutations at the RNA level. These two complementary hypotheses, namely, the diversifying hypothesis and the restorative hypothesis, have distinct predictions on the landscape of RNA editing sites. We collected the chloroplast C-to-U RNA editomes of 21 vascular plants (11 angiosperms, four gymnosperms, and six ferns) from a previous study, aiming to testify whether the plant editomes typically conform to the restorative hypothesis. All predictions made by the restorative hypothesis are verified: (i) nonsynonymous editing sites are more frequent and have higher editing levels than synonymous sites; (ii) nonsynonymous editing levels are extremely high and show weak tissue-specificity in plants; (iii) on the inferred genomic sites with recent T-to-C mutations, nonsynonymous sites but not synonymous sites are compensated by C-to-U RNA editing. In conclusion, nonsynonymous C-to-U RNA editing in plants is adaptive due to its restorative effects. The recoding levels are high and are constantly required across the whole plant so that the recoding events could perfectly mimic DNA mutations. The evolutionary significance of plant RNA editing is systematically demonstrated at the genome-wide level.

Full-Length Transcriptome Profiling of Coridius chinensis Mitochondrial Genome Reveals the Transcription of Genes with Ancestral Arrangement in Insects.

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.

Evolutionary driving forces of A-to-I editing in metazoans.

Adenosine-to-inosine (A-to-I) RNA editing is an evolutionarily conserved mechanism that converts adenosines to inosines in metazoans' transcriptomes. However, the landscapes of editomes have considerably changed during evolution. Here, we review some of our current knowledge of A-to-I editing in the metazoan transcriptomes, focusing on the possible evolutionary driving forces underlying the editing events. First, we review the evolution of ADAR gene family in animals. Then, we summarize the recent advances in characterizing the editomes of various metazoan species. Next, we highlight several factors contributing to the interspecies differences in editomes, including variations in copy number and expression patterns of ADAR genes, the differences in genomic architectures and contents, and the differences in the efficacy of natural selection. After that, we review the possible diversifying and restorative effects of the editing (recoding) events that change the protein sequences. Finally, we discuss the possible convergent evolution of RNA editing in distantly related clades. This article is categorized under: RNA Evolution and Genomics > RNA and Ribonucleoprotein Evolution RNA Processing > RNA Editing and Modification.

A-to-I RNA editing in honeybees shows signals of adaptation and convergent evolution.

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.