Alternative splicing associated with phenotypic plasticity in the bumble bee Bombus terrestris

Major changes in domain arrangements are associated with the evolution of termites

Domains as functional protein units and their rearrangements along the phylogeny can shed light on the functional changes of proteomes associated with the evolution of complex traits like eusociality. This complex trait is associated with sterile soldiers and workers, and long-lived, highly fecund reproductives. Unlike in Hymenoptera (ants, bees, and wasps), the evolution of eusociality within Blattodea, where termites evolved from within cockroaches, was accompanied by a reduction in proteome size, raising the question of whether functional novelty was achieved with existing rather than novel proteins. To address this, we investigated the role of domain rearrangements during the evolution of termite eusociality. Analysing domain rearrangements in the proteomes of three solitary cockroaches and five eusocial termites, we inferred more than 5,000 rearrangements over the phylogeny of Blattodea. The 90 novel domain arrangements that emerged at the origin of termites were enriched for several functions related to longevity, such as protein homeostasis, DNA repair, mitochondrial activity, and nutrient sensing. Many domain rearrangements were related to changes in developmental pathways, important for the emergence of novel castes. Along with the elaboration of social complexity, including permanently sterile workers and larger, foraging colonies, we found 110 further domain arrangements with functions related to protein glycosylation and ion transport. We found an enrichment of caste-biased expression and splicing within rearranged genes, highlighting their importance for the evolution of castes. Furthermore, we found increased levels of DNA methylation among rearranged compared to non-rearranged genes suggesting fundamental differences in their regulation. Our findings indicate the importance of domain rearrangements in the generation of functional novelty necessary for termite eusociality to evolve.

Identification of gene isoforms and their switching events between male and female embryos of the parthenogenetic crustacean Daphnia magna

The cladoceran crustacean Daphnia exhibits phenotypic plasticity, a phenomenon that leads to diverse phenotypes from one genome. Alternative usage of gene isoforms has been considered a key gene regulation mechanism for controlling different phenotypes. However, to understand the phenotypic plasticity of Daphnia, gene isoforms have not been comprehensively analyzed. Here we identified 25,654 transcripts derived from the 9710 genes expressed during environmental sex determination of Daphnia magna using the long-read RNA-Seq with PacBio Iso-Seq. We found that 14,924 transcripts were previously unidentified and 5713 genes produced two or more isoforms. By a combination of Illumina short-read RNA-Seq, we detected 824 genes that implemented switching of the highest expressed isoform between females and males. Among the 824 genes, we found isoform switching of an ortholog of CREB-regulated transcription coactivator, a major regulator of carbohydrate metabolism in animals, and a correlation of this switching event with the sexually dimorphic expression of carbohydrate metabolic genes. These results suggest that a comprehensive catalog of isoforms may lead to understanding the molecular basis for environmental sex determination of Daphnia. We also infer the applicability of the full-length isoform analyses to the elucidation of phenotypic plasticity in Daphnia.

Multidimensional plasticity jointly contributes to rapid acclimation to environmental challenges during biological invasions

Rapid plastic response to environmental changes, which involves extremely complex underlying mechanisms, is crucial for organismal survival during many ecological and evolutionary processes such as those in global change and biological invasions. Gene expression is among the most studied molecular plasticity, while co- or posttranscriptional mechanisms are still largely unexplored. Using a model invasive ascidian Ciona savignyi, we studied multidimensional short-term plasticity in response to hyper- and hyposalinity stresses, covering the physiological adjustment, gene expression, alternative splicing (AS), and alternative polyadenylation (APA) regulations. Our results demonstrated that rapid plastic response varied with environmental context, timescales, and molecular regulatory levels. Gene expression, AS, and APA regulations independently acted on different gene sets and corresponding biological functions, highlighting their nonredundant roles in rapid environmental adaptation. Stress-induced gene expression changes illustrated the use of a strategy of accumulating free amino acids under high salinity and losing/reducing them during low salinity to maintain the osmotic homoeostasis. Genes with more exons were inclined to use AS regulations, and isoform switches in functional genes such as SLC2a5 and Cyb5r3 resulted in enhanced transporting activities by up-regulating the isoforms with more transmembrane regions. The extensive 3'-untranslated region (3'UTR) shortening through APA was induced by both salinity stresses, and APA regulation predominated transcriptomic changes at some stages of stress response. The findings here provide evidence for complex plastic mechanisms to environmental changes, and thereby highlight the importance of systemically integrating different levels of regulatory mechanisms in studying initial plasticity in evolutionary trajectories.

Social insect transcriptomics and the molecular basis of caste diversity

Studies of gene expression provide fundamentally important information on the molecular mechanisms underlying variation in phenotype. Recent technological advances have allowed for the robust study of gene expression through analysis of whole transcriptomes. Here, we review current advances in social insect transcriptomics and discuss their implications in understanding phenotypic diversity. Recent transcriptomic studies provide detailed inventories of the genes involved in producing distinct phenotypes in social species. These investigations have identified key genes and networks involved in producing distinct social insect castes. Nevertheless, questions concerning the evolution of gene expression patterns remain. We suggest a path forward for studying gene expression in future studies of biological systems.

Alternative splicing as a source of phenotypic diversity

A major goal of evolutionary genetics is to understand the genetic processes that give rise to phenotypic diversity in multicellular organisms. Alternative splicing generates multiple transcripts from a single gene, enriching the diversity of proteins and phenotypic traits. It is well established that alternative splicing contributes to key innovations over long evolutionary timescales, such as brain development in bilaterians. However, recent developments in long-read sequencing and the generation of high-quality genome assemblies for diverse organisms has facilitated comparisons of splicing profiles between closely related species, providing insights into how alternative splicing evolves over shorter timescales. Although most splicing variants are probably non-functional, alternative splicing is nonetheless emerging as a dynamic, evolutionarily labile process that can facilitate adaptation and contribute to species divergence.

Extent and complexity of RNA processing in honey bee queen and worker caste development

The distinct honeybee (Apis mellifera) worker and queen castes have become a model for the study of genomic mechanisms of phenotypic plasticity. Here we performed a nanopore-based direct RNA sequencing with exceptionally long reads to compare the mRNA transcripts between queen and workers at three points during their larval development. We found thousands of significantly differentially expressed transcript isoforms (DEIs) between queen and worker larvae. These DEIs were formatted by a flexible splicing system. We showed that poly(A) tails participated in this caste differentiation by negatively regulating the expression of DEIs. Hundreds of isoforms uniquely expressed in either queens or workers during their larval development, and isoforms were expressed at different points in queen and worker larval development demonstrating a dynamic relationship between isoform expression and developmental mechanisms. These findings show the full complexity of RNA processing and transcript expression in honey bee phenotypic plasticity.

•

Honeybee caste differentiation has a complexity of RNA processing

•

Isoforms differentially express between queens and workers during larval development

•

Isoforms are formatted by a flexible alternative splicing system

•

Poly(A) tails are negatively correlated with isoform expression

Alternative splicing in seasonal plasticity and the potential for adaptation to environmental change

Seasonal plasticity is accomplished via tightly regulated developmental cascades that translate environmental cues into trait changes. Little is known about how alternative splicing and other posttranscriptional molecular mechanisms contribute to plasticity or how these mechanisms impact how plasticity evolves. Here, we use transcriptomic and genomic data from the butterfly Bicyclus anynana, a model system for seasonal plasticity, to compare the extent of differential expression and splicing and test how these axes of transcriptional plasticity differ in their potential for evolutionary change. Between seasonal morphs, we find that differential splicing affects a smaller but functionally unique set of genes compared to differential expression. Further, we find strong support for the novel hypothesis that spliced genes are more susceptible than differentially expressed genes to erosion of genetic variation due to selection on seasonal plasticity. Our results suggest that splicing plasticity is especially likely to experience genetic constraints that could affect the potential of wild populations to respond to rapidly changing environments.

Little is known about how alternative splicing and other post-transcriptional molecular mechanisms impact plasticity. Steward et al. use transcriptomic and genomic data from the butterfly Bicyclus anynana, finding that splicing plasticity is likely to experience genetic constraints.

Independent variations in genome-wide expression, alternative splicing, and DNA methylation in brain tissues among castes of the buff-tailed bumblebee, Bombus terrestris

Caste differentiation in social hymenopterans is an intriguing example of phenotypic plasticity. However, the co-ordination among gene regulatory factors to mediate caste differentiation remains inconclusive. In this study, we determined the role of gene regulation and related epigenetic processes in pre-imaginal caste differentiation in the primitively eusocial bumblebee Bombus terrestris. By combining RNA-Seq data from Illumina and PacBio and accurately quantifying methylation at whole-genomic base pair resolution, we found that queens, workers, and drones mainly differentiate in gene expression but not in alternative splicing and DNA methylation. Gynes are the most distinct with the lowest global level of whole-genomic methylation and with the largest number of caste-specific transcripts and alternative splicing events. By contrast, workers exhibit few uniquely expressed or alternatively spliced genes. Moreover, several genes involved in hormone and neurotransmitter metabolism are related to caste differentiation, whereas several neuropeptides are linked with sex differentiation. Despite little genome-wide association among differential gene expression, splicing, and differential DNA methylation, the overlapped gene ontology (GO) terms point to nutrition-related activity. Therefore, variations in gene regulation correlate with the behavioral differences among castes and highlight the specialization of toolkit genes in bumblebee gynes at the beginning of the adult stage.

Comparative Transcriptome Analysis Reveals Sex-Based Differences during the Development of the Adult Parasitic Wasp Cotesia vestalis (Hymenoptera: Braconidae)

The endoparasitic wasp Cotesia vestalis is an important biological agent for controlling the population of Plutella xylostella, a major pest of cruciferous crops worldwide. Though the genome of C. vestalis has recently been reported, molecular mechanisms associated with sexual development have not been comprehensively studied. Here, we combined PacBio Iso-Seq and Illumina RNA-Seq to perform genome-wide profiling of pharate adult and adult development of male and female C. vestalis. Taking advantage of Iso-Seq full-length reads, we identified 14,466 novel transcripts as well as 8770 lncRNAs, with many lncRNAs showing a sex- and stage-specific expression pattern. The differentially expressed gene (DEG) analyses showed 2125 stage-specific and 326 sex-specific expressed genes. We also found that 4819 genes showed 11,856 alternative splicing events through combining the Iso-Seq and RNA-Seq data. The results of comparative analyses showed that most genes were alternatively spliced across developmental stages, and alternative splicing (AS) events were more prevalent in females than in males. Furthermore, we identified six sex-determining genes in this parasitic wasp and verified their sex-specific alternative splicing profiles. Specifically, the characterization of feminizer and doublesex splicing between male and female implies a conserved regulation mechanism of sexual development in parasitic wasps.

Sex-specific expression and DNA methylation in a species with extreme sexual dimorphism and paternal genome elimination

Phenotypic differences between sexes are often mediated by differential expression and alternative splicing of genes. However, the mechanisms that regulate these expression and splicing patterns remain poorly understood. The mealybug, Planococcus citri, displays extreme sexual dimorphism and exhibits an unusual instance of sex-specific genomic imprinting, paternal genome elimination (PGE), in which the paternal chromosomes in males are highly condensed and eliminated from the sperm. Planococcus citri has no sex chromosomes and both sexual dimorphism and PGE are predicted to be under epigenetic control. We recently showed that P. citri females display a highly unusual DNA methylation profile for an insect species, with the presence of promoter methylation associated with lower levels of gene expression. Here, we therefore decided to explore genome-wide differences in DNA methylation between male and female P. citri using whole-genome bisulphite sequencing. We identified extreme differences in genome-wide levels and patterns between the sexes. Males display overall higher levels of DNA methylation which manifest as more uniform low levels across the genome. Whereas females display more targeted high levels of methylation. We suggest these unique sex-specific differences are due to chromosomal differences caused by PGE and may be linked to possible ploidy compensation. Using RNA-Seq, we identify extensive sex-specific gene expression and alternative splicing, but we find no correlation with cis-acting DNA methylation.

Gene expression during larval caste determination and differentiation in intermediately eusocial bumblebees, and a comparative analysis with advanced eusocial honeybees

The queen‐worker caste system of eusocial insects represents a prime example of developmental polyphenism (environmentally‐induced phenotypic polymorphism) and is intrinsic to the evolution of advanced eusociality. However, the comparative molecular basis of larval caste determination and subsequent differentiation in the eusocial Hymenoptera remains poorly known. To address this issue within bees, we profiled caste‐associated gene expression in female larvae of the intermediately eusocial bumblebee Bombus terrestris. In B. terrestris, female larvae experience a queen‐dependent period during which their caste fate as adults is determined followed by a nutrition‐sensitive period also potentially affecting caste fate but for which the evidence is weaker. We used mRNA‐seq and qRT‐PCR validation to isolate genes differentially expressed between each caste pathway in larvae at developmental stages before and after each of these periods. We show that differences in gene expression between caste pathways are small in totipotent larvae, then peak after the queen‐dependent period. Relatively few novel (i.e., taxonomically‐restricted) genes were differentially expressed between castes, though novel genes were significantly enriched in late‐instar larvae in the worker pathway. We compared sets of caste‐associated genes in B. terrestris with those reported from the advanced eusocial honeybee, Apis mellifera, and found significant but relatively low levels of overlap of gene lists between the two species. These results suggest both the existence of low numbers of shared toolkit genes and substantial divergence in caste‐associated genes between Bombus and the advanced eusocial Apis since their last common eusocial ancestor.

High-quality chromosome-level genome assembly and full-length transcriptome analysis of the pharaoh ant Monomorium pharaonis

Background

Ants with complex societies have fascinated scientists for centuries. Comparative genomic and transcriptomic analyses across ant species and castes have revealed important insights into the molecular mechanisms underlying ant caste differentiation. However, most current ant genomes and transcriptomes are highly fragmented and incomplete, which hinders our understanding of the molecular basis for complex ant societies.

Findings

By hybridizing Illumina, Pacific Biosciences, and Hi-C sequencing technologies, we de novo assembled a chromosome-level genome for Monomorium pharaonis, with a scaffold N50 of 27.2 Mb. Our new assembly provides better resolution for the discovery of genome rearrangement events at the chromosome level. Analysis of full-length isoform sequencing (ISO-seq) suggested that ∼15 Gb of ISO-seq data were sufficient to cover most expressed genes, but the number of transcript isoforms steadily increased with sequencing data coverage. Our high-depth ISO-seq data greatly improved the quality of gene annotation and enabled the accurate detection of alternative splicing isoforms in different castes of M. pharaonis. Comparative transcriptome analysis across castes based on the ISO-seq data revealed an unprecedented number of transcript isoforms, including many caste-specific isoforms. We also identified a number of conserved long non-coding RNAs that evolved specifically in ant lineages and several that were conserved across insect lineages.

Conclusions

We produced a high-quality chromosome-level genome for M. pharaonis, which significantly improved previous short-read assemblies. Together with full-length transcriptomes for all castes, we generated a highly accurate annotation for this ant species. These long-read sequencing results provide a useful resource for future functional studies on the genetic mechanisms underlying the evolution of social behaviors and organization in ants.

Evidence of capacitation in the parasitoid wasp, Nasonia vitripennis, and its potential role in sex allocation

The allocation of resources to the production of one sex or another has been observed in a large variety of animals. Its theoretical basis allows accurate predictions of offspring sex ratios in many species, but the mechanisms by which sex allocation is controlled are poorly understood. Using previously published data, we investigated whether alternative splicing, combined with differential gene expression, was involved with sex allocation in the parasitoid wasp, Nasonia vitripennis. We found that sex allocation is not controlled by alternative splicing but changes in gene and transcript-specific expression, which were identified to be involved with oviposition, were shown to be similar to those involved in sperm motility and capacitation. Genes involved in cholesterol efflux, a key component of capacitation, along with calcium transport, neurotransmission, trypsin, and MAPKinase activity were regulated in ovipositing wasps. The results show evidence for regulation of sperm motility and of capacitation in an insect which, in the context of the physiology of the N. vitripennis spermatheca, could be important for sex allocation.

Multiple plasticity regulators reveal targets specifying an induced predatory form in nematodes

The ability to translate a single genome into multiple phenotypes, or developmental plasticity, defines how phenotype derives from more than just genes. However, to study the evolutionary targets of plasticity and their evolutionary fates, we need to understand how genetic regulators of plasticity control downstream gene expression. Here, we have identified a transcriptional response specific to polyphenism (i.e., discrete plasticity) in the nematode Pristionchus pacificus. This species produces alternative resource-use morphs - microbivorous and predatory forms, differing in the form of their teeth, a morphological novelty - as influenced by resource availability. Transcriptional profiles common to multiple polyphenism-controlling genes in P. pacificus reveal a suite of environmentally sensitive loci, or ultimate target genes, that make up an induced developmental response. Additionally, in vitro assays show that one polyphenism regulator, the nuclear receptor (NR) NHR-40, physically binds to promoters with putative HNF4⍺ (the NR class including NHR-40) binding sites, suggesting this receptor may directly regulate genes that describe alternative morphs. Among differentially expressed genes were morph-limited genes, highlighting factors with putative "on-off" function in plasticity regulation. Further, predatory morph-biased genes included candidates - namely, all four P. pacificus homologs of Hsp70, which have HNF4⍺ motifs - whose natural variation in expression matches phenotypic differences among P. pacificus wild isolates. In summary, our study links polyphenism regulatory loci to the transcription producing alternative forms of a morphological novelty. Consequently, our findings establish a platform for determining how specific regulators of morph-biased genes may influence selection on plastic phenotypes.

Morphological characterization and staging of bumble bee pupae

Bumble bees (Hymenoptera: Apidae, Bombus) are important pollinators and models for studying mechanisms underlying developmental plasticity, such as factors influencing size, immunity, and social behaviors. Research on such processes, as well as expanding use of gene-manipulation and gene expression technologies, requires a detailed understanding of how these bees develop. Developmental research often uses time-staging of pupae, however dramatic size differences in these bees can generate variation in developmental timing. To study developmental mechanisms in bumble bees, appropriate staging of developing bees using morphology is necessary. In this study, we describe morphological changes across development in several bumble bee species and use this to establish morphology-based staging criteria, establishing 20 distinct illustrated stages. These criteria, defined largely by eye and cuticle pigmentation patterns, are generalizable across members of the subgenus Pyrobombus, and can be used as a framework for study of other bumble bee subgenera. We examine the effects of temperature, caste, size, and species on pupal development, revealing that pupal duration shifts with each of these factors, confirming the importance of staging pupae based on morphology rather than age and the need for standardizing sampling.

Extensive Differential Splicing Underlies Phenotypically Plastic Aphid Morphs

Phenotypic plasticity results in a diversity of phenotypes from a single genotype in response to environmental cues. To understand the molecular basis of phenotypic plasticity, studies have focused on differential gene expression levels between environmentally determined phenotypes. The extent of alternative splicing differences among environmentally determined phenotypes has largely been understudied. Here, we study alternative splicing differences among plastically produced morphs of the pea aphid using RNA-sequence data. Pea aphids express two separate polyphenisms (plasticity with discrete phenotypes): a wing polyphenism consisting of winged and wingless females and a reproduction polyphenism consisting of asexual and sexual females. We find that pea aphids alternatively splice 34% of their genes, a high percentage for invertebrates. We also find that there is extensive use of differential spliced events between genetically identical, polyphenic females. These differentially spliced events are enriched for exon skipping and mutually exclusive exon events that maintain the open reading frame, suggesting that polyphenic morphs use alternative splicing to produce phenotype-biased proteins. Many genes that are differentially spliced between polyphenic morphs have putative functions associated with their respective phenotypes. We find that the majority of differentially spliced genes is not differentially expressed genes. Our results provide a rich candidate gene list for future functional studies that would not have been previously considered based solely on gene expression studies, such as ensconsin in the reproductive polyphenism, and CAKI in the wing polyphenism. Overall, this study suggests an important role for alternative splicing in the expression of environmentally determined phenotypes.

Vitellogenin-like A-associated shifts in social cue responsiveness regulate behavioral task specialization in an ant

Division of labor and task specialization explain the success of human and insect societies. Social insect colonies are characterized by division of labor, with workers specializing in brood care early and foraging later in life. Theory posits that this task switching requires shifts in responsiveness to task-related cues, yet experimental evidence is weak. Here, we show that a Vitellogenin (Vg) ortholog identified in an RNAseq study on the ant T. longispinosus is involved in this process: using phylogenetic analyses of Vg and Vg-like genes, we firstly show that this candidate gene does not cluster with the intensively studied honey bee Vg but falls into a separate Vg-like A cluster. Secondly, an experimental knockdown of Vg-like A in the fat body caused a reduction in brood care and an increase in nestmate care in young ant workers. Nestmate care is normally exhibited by older workers. We demonstrate experimentally that this task switch is at least partly based on Vg-like A-associated shifts in responsiveness from brood to worker cues. We thus reveal a novel mechanism leading to early behavioral maturation via changes in social cue responsiveness mediated by Vg-like A and associated pathways, which proximately play a role in regulating division of labor.