Contrasting Sex-and Caste-Dependent piRNA Profiles in the Transposon Depleted Haplodiploid Honeybee Apis mellifera

Evolution of piRNA-guided defense against transposable elements

Transposable elements (TEs) shape every aspect of genome biology, influencing genome stability, size, and organismal fitness. Following the 2007 discovery of the piRNA defense system, researchers have made numerous findings about organisms' defenses against these genomic invaders. TEs are suppressed by a 'genomic immune system', where TE insertions within specialized regions called PIWI-interacting RNA (piRNA) clusters produce small RNAs responsible for their suppression. The evolution of piRNA clusters and the piRNA system is only now being understood, largely because most research has been conducted in developmental biology labs using only one to two genotypes of Drosophila melanogaster. While piRNAs themselves were identified simultaneously in various organisms (flies, mice, rats, and zebrafish) in 2006-2007, detailed work on piRNA clusters has only recently expanded beyond D. melanogaster. By studying piRNA cluster evolution in various organisms from an evolutionary perspective, we are beginning to understand more about TE suppression mechanisms and organism-TE coevolution.

The eusocial non-code: Unveiling the impact of noncoding RNAs on Hymenoptera eusocial evolution

Eusociality, characterized by reproductive division of labor, cooperative brood care, and multi-generational cohabitation, represents a pinnacle of complex social evolution, most notably manifested within the Hymenoptera order including bees, ants, and wasps. The molecular underpinnings underlying these sophisticated social structures remain an enigma, with noncoding RNAs (ncRNAs) emerging as crucial regulatory players. This article delves into the roles of ncRNAs in exerting epigenetic control during the development and maintenance of Hymenopteran eusociality. We consolidate current findings on various classes of ncRNAs, underscoring their influence on gene expression regulation pertinent to caste differentiation, developmental plasticity, and behavioral modulation. Evidence is explored supporting the hypothesis that ncRNAs contribute to epigenetic landscapes fostering eusocial traits through genomic regulation. They are likely to play an important role in eusociality "point of no return". Critical analysis is provided on the functional insights garnered from ncRNA profiles correlated with caste-specific phenotypes, specifical for phylogenetic branches and transitional sociality models, drawing from comparative genomics and transcriptomics studies. Overall, ncRNA provides a missed understanding of both "genetic toolkit" and "unique genes" hypotheses of eusociality development. Moreover, it points to gaps in current knowledge, advocating for integrative approaches combining genomics, proteomics, and epigenetics to decipher the complexity of eusociality. Understanding the ncRNA contributions offers not only a window into the molecular intricacies of Hymenoptera sociality but also extends our comprehension of how complex biological systems evolve and function.

Analysis of the Expression Patterns of piRNAs in Response to Microsporidian Invasion in Midgut of Workers (Apis cerana cerana)

Piwi-interacting RNAs (piRNAs) play an essential part in transposon suppression, DNA methylation, and antiviral responses. The current understanding of the roles of piRNAs in honeybees is very limited. This study aims to analyze the expression pattern and regulatory role of piRNAs in the Asian honeybee (Apis cerana) responding to infection by Nosema ceranae, based on previously gained small RNA-seq data. Here, 450 and 422 piRNAs were respectively identified in the midgut tissues of Apis cerana cerana workers at 7 and 10 days post-inoculation (dpi) with N. ceranae, including 539 non-redundant ones. Additionally, one up-regulated (piR-ace-1216942) and one down-regulated (piR-ace-776728) piRNA were detected in the workers' midgut at 7 dpi, targeting 381 mRNAs involved in 31 GO terms, such as metabolic processes, catalytic activity, and organelles, as well as 178 KEGG pathways, including lysosome, MAPK signaling pathway, and purine metabolism. A total of 35 up-regulated and 11 down-regulated piRNAs were screened from the workers' midgut at 10 dpi, targeting 13,511 mRNAs engaged in 50 GO terms, such as biological regulation, transporter activity, and membrane, as well as 389 KEGG pathways, including the JAK-STAT signaling pathway, Hippo signaling pathway, and nitrogen metabolism. Further analysis indicated that 28 differentially expressed piRNAs (DEpiRNAs) in the midgut at 10 dpi could target 299 mRNAs annotated to three cellular immune pathways (lysosome, endocytosis, and phagosome), while 24 DEpiRNAs could target 205 mRNAs relevant to four humoral immune pathways (FoxO, JAK-STAT, NF-κB, and MAPK signaling pathway). Through Sanger sequencing and RT-qPCR, the expression of six randomly selected DEpiRNAs was verified. Moreover, the dual-luciferase reporter gene assay confirmed the binding relationships between piR-ace-446232 and CRT as well as between piR-ace-1008436 and EGFR. Our findings not only contribute to enrich our understanding of the role of piRNAs in honeybees but also provide a basis for exploring the host response to N. ceranae infection mediated by piRNAs.

Non-gonadal expression of piRNAs is widespread across Arthropoda

PIWI-interacting RNAs (piRNAs) were discovered in the early 2000s and became known for their role in protecting the germline genome against mobile genetic elements. Successively, piRNAs were also detected in the somatic cells of gonads in multiple animal species. In recent years, piRNAs have been reported in non-gonadal tissues in various arthropods, contrary to the initial assumptions of piRNAs being exclusive to gonads. Here, we performed an extensive literature review, which revealed that reports on non-gonadal somatic piRNA expression are not limited to a few specific species. Instead, when multiple studies are considered collectively, it appears to be a widespread phenomenon across arthropods. Furthermore, we systematically analyzed 168 publicly available small RNA-seq datasets from diverse tissues in 17 species, which further supported the bibliographic reports that piRNAs are expressed across tissues and species in Arthropoda.

Embryonic piRNAs target horizontally transferred vertebrate transposons in assassin bugs

Introduction

Piwi proteins and the associated Piwi-interacting RNAs (piRNAs) coordinate a surveillance system that protects the animal genome from DNA damage induced by transposable element (TE) mobilization. While the pathway has been described in detail in the fruit fly Drosophila melanogaster, much less is known in more basal insects. Rhodnius prolixus is an hemipteran insect and one of the major vectors of Chagas disease. Rhodnius acquired specific classes of horizontally transferred transposons (HTTs) by feeding on bats, opossums and squirrel monkeys, thus providing the opportunity to investigate the piRNA-base response against HTTs in this species.

Methods

SmallRNA-Seq reads mapping to HTTs and resident transposable elements were quantified and checked for piRNA features like 1U a 10A biases, ping-pong and phasing signatures. Uniquely mapped piRNAs were used to identify piRNA clusters in Rhodnius' genome. RNA-Seq data was used to quantify transposon and Rp-PIWI genes expression levels and were validated by qRT-PCR.

Results

By analyzing the temporal dynamics of piRNA cluster expression and piRNA production during critical stages of Rhodnius development, we show that peak levels of ∼28 nt long piRNAs correlate with reduced HTT and resident TE expression primarily during embryogenesis. Strikingly, while resident TEs piRNAs seem to engage in a typical ping-pong amplification mechanism, sense and antisense HTT piRNAs instead overlap by ∼20 nt or do not display ping-pong signatures.

Discussion

Our data shed light on the biogenesis and functions of the piRNAs in Rhodnius prolixus and reveal that piRNAs, but not the siRNA pathway, responded to HTTs that were recently transferred from vertebrate tetrapods to a hematophagous insect of medical relevance.

Recombination Rate Variation in Social Insects: An Adaptive Perspective

Social insects have the highest rates of meiotic recombination among Metazoa, but there is considerable variation within the Hymenoptera. We synthesize the literature to investigate several hypotheses for these elevated recombination rates. We reexamine the long-standing Red Queen hypothesis, considering how social aspects of immunity could lead to increases in recombination. We examine the possibility of positive feedback between gene duplication and recombination rate in the context of caste specialization. We introduce a novel hypothesis that recombination rate may be driven up by direct selection on recombination activity in response to increases in lifespan. Finally, we find that the role of population size in recombination rate evolution remains opaque, despite the long-standing popularity of this hypothesis. Moreover, our review emphasizes how the varied life histories of social insect species provide an effective framework for advancing a broader understanding of adaptively driven variation in recombination rates.

Ago4-piRNA complex is a key component of genomic immune system against transposon expression in Penaeus monodon

Argonaute proteins are key constituents of small RNA-guided regulatory pathways. In crustaceans, members of the AGO subfamily of Argonaute proteins that play vital roles in immune defense are well studied, while proteins of the PIWI subfamily are less established. PmAgo4 of the black tiger shrimp, Penaeus monodon, though phylogenetically clustered with the AGO subfamily, has distinctive roles of the PIWI subfamily in safeguarding the genome from transposon invasion and controlling germ cell development. This study explored a molecular mechanism by which PmAgo4 regulates transposon expression in the shrimp germline. PmAgo4-associated small RNAs were co-immunoprecipitated from shrimp testis lysate using a PmAgo4-specific polyclonal antibody. RNA-seq revealed a majority of 26-27 nt long small RNAs in the PmAgo4-IP fraction suggesting that PmAgo4 is predominantly associated with piRNAs. Mapping of these piRNAs on nucleotide sequences of two gypsy and a mariner-like transposons of P. monodon suggested that most piRNAs were originated from the antisense strand of transposons. Suppression of PmAgo4 expression by a specific dsRNA elevated the expression levels of the three transposons while decreasing the levels of transposon-related piRNAs. Taken together, these results imply that PmAgo4 exerts its suppressive function on transposons by controlling the biogenesis of transposon-related piRNAs and thus, provides a defense mechanism against transposon invasion in shrimp germline cells.

Age- and caste-independent piRNAs in the germline and miRNA profiles linked to caste and fecundity in the ant Temnothorax rugatulus

Social insects are models for studies of phenotypic plasticity. Ant queens and workers vary in fecundity and lifespan, which are enhanced and extended in queens. Yet, the regulatory mechanisms underlying this variation are not well understood. Ant queens live and reproduce for years, so that they need to protect their germline from transposable element (TE) activity, which may be redundant in short-lived, often sterile workers. We analysed the expression of two protective classes of small RNAs, microRNAs (miRNAs) and Piwi-interacting RNAs (piRNAs), in various tissues, castes and age classes of the ant Temnothorax rugatulus. In queens, piRNAs were highly abundant in ovaries with TEs being their clear targets, with reduced but still detectable piRNA-specific ping-pong signatures in thorax and brains. piRNA pathway activity varied little with age in queens. Moreover, the reduced ovaries of workers also exhibited similar piRNA activity and this not only in young, fertile workers, but also in older foragers with regressed ovaries. Therefore, these ants protect their germline through piRNA activity, regardless of ovarian development, age or caste, even in sterile workers often considered the soma of the superorganism. Our tissue-specific miRNA analysis detected the expression of 304 miRNAs, of which 105 were expressed in all tissues, 10 enriched in the brain, three in the thorax, whereas 83 were ovarian-specific. We identified ovarian miRNAs whose expression was related to caste, fecundity and age, and which likely regulate group-specific gene expression. sRNA shifts in young- to middle-aged queens were minor, suggesting delayed senescence in this reproductive caste.

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

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

First Identification and Investigation of piRNAs in the Larval Gut of the Asian Honeybee, Apis cerana

Piwi-interacting RNAs (piRNAs), a class of small non-coding RNAs (ncRNAs), play pivotal roles in maintaining the genomic stability and modulating biological processes such as growth and development via the regulation of gene expression. However, the piRNAs in the Asian honeybee (Apis cerana) are still largely unknown at present. In this current work, on the basis of previously gained high-quality small RNA-seq datasets, piRNAs in the larval gut of Apis cerana cerana, the nominated species of A. cerana, were identified for the first time, followed by an in-depth investigation of the regulatory roles of differentially expressed piRNAs (DEpiRNAs) in the developmental process of the A. c. cerana. Here, a total of 621 piRNAs were identified in A. c. cerana larval guts, among which 499 piRNAs were shared by 4-(Ac4 group), 5-(Ac5 group), and 6-day-old (Ac6 group) larval guts, while the numbers of unique ones equaled 79, 37, and 11, respectively. The piRNAs in each group ranged from 24 nucleotides (nt) to 33 nt in length, and the first base of the piRNAs had a cytosine (C) bias. Additionally, five up-regulated and five down-regulated piRNAs were identified in the Ac4 vs. Ac5 comparison group, nine of which could target 9011 mRNAs; these targets were involved in 41 GO terms and 137 pathways. Comparatively, 22 up-regulated piRNAs were detected in the Ac5 vs. Ac6 comparison group, 21 of which could target 28,969 mRNAs; these targets were engaged in 46 functional terms and 164 pathways. The results suggested an overall alteration of the expression pattern of piRNAs during the developmental process of A. c. cerana larvae. The regulatory network analysis showed that piR-ace-748815 and piR-ace-512574 in the Ac4 vs. Ac5 comparison group as well as piR-ace-716466 and piR-ace-828146 in the Ac5 vs. Ac6 comparison group were linked to the highest number of targets. Further investigation indicated that targets of DEpiRNAs in the abovementioned two comparison groups could be annotated to several growth and development-associated pathways, such as the Jak/STAT, TGF-β, and Wnt signaling pathways, indicating the involvement of DEpiRNAs in modulating larval gut development via these crucial pathways. Moreover, the expression trends of six randomly selected DEpiRNAs were verified using a combination of stem-loop RT-PCR and RT-qPCR. These results not only provide a novel insight into the development of the A. c. cerana larval gut, but also lay a foundation for uncovering the epigenetic mechanism underlying larval gut development.

Abundant small RNAs in the reproductive tissues and eggs of the honey bee, Apis mellifera

Background

Polyandrous social insects such as the honey bee are prime candidates for parental manipulation of gene expression in offspring. Although there is good evidence for parent-of-origin effects in honey bees the epigenetic mechanisms that underlie these effects remain a mystery. Small RNA molecules such as miRNAs, piRNAs and siRNAs play important roles in transgenerational epigenetic inheritance and in the regulation of gene expression during development.

Results

Here we present the first characterisation of small RNAs present in honey bee reproductive tissues: ovaries, spermatheca, semen, fertilised and unfertilised eggs, and testes. We show that semen contains fewer piRNAs relative to eggs and ovaries, and that piRNAs and miRNAs which map antisense to genes involved in DNA regulation and developmental processes are differentially expressed between tissues. tRNA fragments are highly abundant in semen and have a similar profile to those seen in the semen of other animals. Intriguingly we also find abundant piRNAs that target the sex determination locus, suggesting that piRNAs may play a role in honey bee sex determination.

Conclusions

We conclude that small RNAs may play a fundamental role in honey bee gametogenesis and reproduction and provide a plausible mechanism for parent-of-origin effects on gene expression and reproductive physiology.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-022-08478-9.

Transposable element repression using piRNAs, and its relevance to endogenous viral elements (EVEs) and immunity in insects

The piRNA system controls transposable element (TE) mobility by transcriptional gene silencing and post-transcriptional gene silencing. Dispersed in insect genomes, piRNA clusters contain TE copies, from which they produce piRNAs (specific small RNAs). These piRNAs can both target the nascent transcripts produced by active TE copies and directly repress them by heterochromatinization. They can also target mature transcripts and cleave them following amplification by the so-called 'ping-pong' loop mechanism. Moreover, piRNA clusters contain endogenous viral elements (EVEs), from which they produce piRNAs. The current idea is that these piRNAs could participate in the antiviral response against exogenous viral infection. In this review, we show that among insects, to date, this antiviral response by the piRNA system appears mainly restricted to mosquitoes, but this could be due to the focus of most studies on arboviruses.

The genomic basis of evolutionary differentiation among honey bees

In contrast to the western honey bee, Apis mellifera, other honey bee species have been largely neglected despite their importance and diversity. The genetic basis of the evolutionary diversification of honey bees remains largely unknown. Here, we provide a genome-wide comparison of three honey bee species, each representing one of the three subgenera of honey bees, namely the dwarf (Apis florea), giant (A. dorsata), and cavity-nesting (A. mellifera) honey bees with bumblebees as an outgroup. Our analyses resolve the phylogeny of honey bees with the dwarf honey bees diverging first. We find that evolution of increased eusocial complexity in Apis proceeds via increases in the complexity of gene regulation, which is in agreement with previous studies. However, this process seems to be related to pathways other than transcriptional control. Positive selection patterns across Apis reveal a trade-off between maintaining genome stability and generating genetic diversity, with a rapidly evolving piRNA pathway leading to genomes depleted of transposable elements, and a rapidly evolving DNA repair pathway associated with high recombination rates in all Apis species. Diversification within Apis is accompanied by positive selection in several genes whose putative functions present candidate mechanisms for lineage-specific adaptations, such as migration, immunity, and nesting behavior.

Genome-wide profiling of piRNAs in the whitefly Bemisia tabaci reveals cluster distribution and association with begomovirus transmission

The whitefly Bemisia tabaci MEAM1 is a notorious vector capable of transmitting many plant viruses, resulting in serious crop loss and food shortage around the world. To investigate potential sRNA-mediated regulatory mechanisms in whiteflies that are affected by virus acquisition and transmission, we conducted small RNA (sRNA) deep sequencing and performed genome-wide profiling of piwi-interacting RNAs (piRNAs) in whiteflies that were fed on tomato yellow leaf curl virus (TYLCV)-infected or non-infected tomato plants for 24, 48, and 72 h. In the present study, piRNA reads ranging from 564,395 to 1,715,652 per library were identified and shown to distribute unevenly in clusters (57 to 96 per library) on the whitefly (B. tabaci MEAM1) genome. Among them, 53 piRNA clusters were common for all treatments. Comparative analysis between libraries generated from viruliferous and non-viruliferous whiteflies identified five TYLCV-induced and 24 TYLCV-suppressed piRNA clusters. Approximately 62% of piRNAs were derived from non-coding sequences including intergenic regions, introns, and untranslated regions (UTRs). The remaining 38% were derived from coding sequences (CDS) or repeat elements. Interestingly, six protein coding genes were targeted by the TYLCV-induced piRNAs. We identified a large number of piRNAs that were distributed in clusters across the whitefly genome, with 60% being derived from non-coding regions. Comparative analysis revealed that feeding on a virus-infected host caused induction and suppression of only a small number of piRNA clusters in whiteflies. Although piRNAs primarily regulate the activity of transposable elements, our results suggest that they may have additional functions in regulating protein coding genes and in insect-virus interactions.

Transcriptomic and functional analyses of the piRNA pathway in the Chagas disease vector Rhodnius prolixus

The piRNA pathway is a surveillance system that guarantees oogenesis and adult fertility in a range of animal species. The pathway is centered on PIWI clade Argonaute proteins and the associated small non-coding RNAs termed piRNAs. In this study, we set to investigate the evolutionary conservation of the piRNA pathway in the hemimetabolous insect Rhodnius prolixus. Our transcriptome profiling reveals that core components of the pathway are expressed during previtellogenic stages of oogenesis. Rhodnius’ genome harbors four putative piwi orthologs. We show that Rp-piwi2, Rp-piwi3 and Rp-ago3, but not Rp-piwi1 transcripts are produced in the germline tissues and maternally deposited in the mature eggs. Consistent with a role in Rhodnius oogenesis, parental RNAi against the Rp-piwi2, Rp-piwi3 and Rp-ago3 results in severe egg laying and female adult fertility defects. Furthermore, we show that the reduction of the Rp-piwi2 levels by parental RNAi disrupts oogenesis by causing a dramatic loss of trophocytes, egg chamber degeneration and oogenesis arrest. Intriguingly, the putative Rp-Piwi2 protein features a polyglutamine tract at its N-terminal region, which is conserved in PIWI proteins encoded in the genome of other Triatomine species. Together with R. prolixus, these hematophagous insects are primary vectors of the Chagas disease. Thus, our data shed more light on the evolution of the piRNA pathway and provide a framework for the development of new control strategies for Chagas disease insect vectors.

Rhodnius prolixus together with other blood-feeding bugs of the Triatominae family are primary vectors of the protozoan Trypanosoma cruzi, the causative agent of the Chagas disease. It has been estimated that 7–8 million people are affected by this life-threatening illness worldwide, which makes the Chagas disease one of the most neglected tropical diseases. In this study, we describe the transcriptome of previtellogenic stages of Rhodnius oogenesis. Furthermore, by using a combination of molecular biology techniques and functional analyses we show that central components of the piRNA pathway are conserved in this species. The piRNA pathway guarantees genomic stability in the germ cells of organisms as distant as flies and mice. In accordance, we find that the knock-down of the piwi genes, which form the backbone of the pathway, results in partial or complete female adult sterility in Rhodnius. Our data will help improve the annotation of the Rhodnius genome and provide a framework for the development of novel techniques aiming at the eradication of Rhodnius prolixus and other Triatomine species from the infested areas. The achievement of this goal will ultimately prevent the transmission of trypanosomes to humans and reduce or eliminate the diffusion of the Chagas disease.

Detecting rare asymmetrically methylated cytosines and decoding methylation patterns in the honeybee genome

Context-dependent gene expression in eukaryotes is controlled by several mechanisms including cytosine methylation that primarily occurs in the CG dinucleotides (CpGs). However, less frequent non-CpG asymmetric methylation has been found in various cell types, such as mammalian neurons, and recent results suggest that these sites can repress transcription independently of CpG contexts. In addition, an emerging view is that CpG hemimethylation may arise not only from deregulation of cellular processes but also be a standard feature of the methylome. Here, we have applied a novel approach to examine whether asymmetric CpG methylation is present in a sparsely methylated genome of the honeybee, a social insect with a high level of epigenetically driven phenotypic plasticity. By combining strand-specific ultra-deep amplicon sequencing of illustrator genes with whole-genome methylomics and bioinformatics, we show that rare asymmetrically methylated CpGs can be unambiguously detected in the honeybee genome. Additionally, we confirm differential methylation between two phenotypically and reproductively distinct castes, queens and workers, and offer new insight into the heterogeneity of brain methylation patterns. In particular, we challenge the assumption that symmetrical methylation levels reflect symmetry in the underlying methylation patterns and conclude that hemimethylation may occur more frequently than indicated by methylation levels. Finally, we question the validity of a prior study in which most of cytosine methylation in this species was reported to be asymmetric.

Detecting rare asymmetrically methylated cytosines and decoding methylation patterns in the honeybee genome

Context-dependent gene expression in eukaryotes is controlled by several mechanisms including cytosine methylation that primarily occurs in the CG dinucleotides (CpGs). However, less frequent non-CpG asymmetric methylation has been found in various cell types, such as mammalian neurons, and recent results suggest that these sites can repress transcription independently of CpG contexts. In addition, an emerging view is that CpG hemimethylation may arise not only from deregulation of cellular processes but also be a standard feature of the methylome. Here, we have applied a novel approach to examine whether asymmetric CpG methylation is present in a sparsely methylated genome of the honeybee, a social insect with a high level of epigenetically driven phenotypic plasticity. By combining strand-specific ultra-deep amplicon sequencing of illustrator genes with whole-genome methylomics and bioinformatics, we show that rare asymmetrically methylated CpGs can be unambiguously detected in the honeybee genome. Additionally, we confirm differential methylation between two phenotypically and reproductively distinct castes, queens and workers, and offer new insight into the heterogeneity of brain methylation patterns. In particular, we challenge the assumption that symmetrical methylation levels reflect symmetry in the underlying methylation patterns and conclude that hemimethylation may occur more frequently than indicated by methylation levels. Finally, we question the validity of a prior study in which most of cytosine methylation in this species was reported to be asymmetric.