Evolutionary Rates of Bumblebee Genomes Are Faster at Lower Elevations

Genome-wide analysis reveals the contributors to fast molecular evolution of the Chinese hook snout carp (Opsariichthys bidens)

Variations in molecular evolutionary rate have been widely investigated among lineages and genes. However, it remains an open question whether fast rate of molecular evolution is driven by natural selection or random drift, and how the fast rate is linked to metabolic rate. Additionally, previous studies on fast molecular evolution have been largely restricted to concatenated matrix of genes or a few specifically selected genes, but less is known for individual genes at the genome-wide level. Here we addressed these questions using more than 5000 single-copy orthologous (SCO) genes through comparative genomic and phylogenetic analyses among fishes, with a special focus on a newly-sequenced clupeocephalan fish the Chinese hook snout carp Opsariichthys bidens. We showed O. bidens displays significantly higher mean substitution rate and more fast-evolving SCO genes (2172 genes) than most fishes studied here. The rapidly evolving genes are enriched in highly conserved and very basic functions such as translation and ribosome that are critical for biological fitness. We further revealed that ∼25 % of these fast-evolving genes exhibit a constant increase of substitution rate from the common ancestor down to the present, suggesting a neglected but important contribution from ancestral states. Model fitting showed that ∼85 % of fast-evolving genes exclusive to O. bidens and related species follow the adaptive evolutionary model rather than random-drift model, and 7.6 % of fast-evolving genes identified in O. bidens have experienced positive selection, both indicating the reflection of adaptive selection. Finally, metabolic rate was observed to be linked with substitution rate in a gene-specific manner. Overall, our findings reveal fast molecular evolution of SCO genes at genome-wide level in O. bidens, and uncover the evolutionary and ecological contributors to it.

Admixture of evolutionary rates across a butterfly hybrid zone

Hybridization is a major evolutionary force that can erode genetic differentiation between species, whereas reproductive isolation maintains such differentiation. In studying a hybrid zone between the swallowtail butterflies Papilio syfanius and Papilio maackii (Lepidoptera: Papilionidae), we made the unexpected discovery that genomic substitution rates are unequal between the parental species. This phenomenon creates a novel process in hybridization, where genomic regions most affected by gene flow evolve at similar rates between species, while genomic regions with strong reproductive isolation evolve at species-specific rates. Thus, hybridization mixes evolutionary rates in a way similar to its effect on genetic ancestry. Using coalescent theory, we show that the rate-mixing process provides distinct information about levels of gene flow across different parts of genomes, and the degree of rate-mixing can be predicted quantitatively from relative sequence divergence ([Formula: see text]) between the hybridizing species at equilibrium. Overall, we demonstrate that reproductive isolation maintains not only genomic differentiation, but also the rate at which differentiation accumulates. Thus, asymmetric rates of evolution provide an additional signature of loci involved in reproductive isolation.

Investigating the reliability of molecular estimates of evolutionary time when substitution rates and speciation rates vary

Background

An accurate timescale of evolutionary history is essential to testing hypotheses about the influence of historical events and processes, and the timescale for evolution is increasingly derived from analysis of DNA sequences. But variation in the rate of molecular evolution complicates the inference of time from DNA. Evidence is growing for numerous factors, such as life history and habitat, that are linked both to the molecular processes of mutation and fixation and to rates of macroevolutionary diversification. However, the most widely used methods rely on idealised models of rate variation, such as the uncorrelated and autocorrelated clocks, and molecular dating methods are rarely tested against complex models of rate change. One relationship that is not accounted for in molecular dating is the potential for interaction between molecular substitution rates and speciation, a relationship that has been supported by empirical studies in a growing number of taxa. If these relationships are as widespread as current evidence suggests, they may have a significant influence on molecular dates.

Results

We simulate phylogenies and molecular sequences under three different realistic rate variation models—one in which speciation rates and substitution rates both vary but are unlinked, one in which they covary continuously and one punctuated model in which molecular change is concentrated in speciation events, using empirical case studies to parameterise realistic simulations. We test three commonly used “relaxed clock” molecular dating methods against these realistic simulations to explore the degree of error in molecular dates under each model. We find average divergence time inference errors ranging from 12% of node age for the unlinked model when reconstructed under an uncorrelated rate prior using BEAST 2, to up to 91% when sequences evolved under the punctuated model are reconstructed under an autocorrelated prior using PAML.

Conclusions

We demonstrate the potential for substantial errors in molecular dates when both speciation rates and substitution rates vary between lineages. This study highlights the need for tests of molecular dating methods against realistic models of rate variation generated from empirical parameters and known relationships.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12862-022-02015-8.

Genomic Signatures of Recent Adaptation in a Wild Bumblebee

Environmental changes threaten insect pollinators, creating risks for agriculture and ecosystem stability. Despite their importance, we know little about how wild insects respond to environmental pressures. To understand the genomic bases of adaptation in an ecologically important pollinator, we analyzed genomes of Bombus terrestris bumblebees collected across Great Britain. We reveal extensive genetic diversity within this population, and strong signatures of recent adaptation throughout the genome affecting key processes including neurobiology and wing development. We also discover unusual features of the genome, including a region containing 53 genes that lacks genetic diversity in many bee species, and a horizontal gene transfer from a Wolbachia bacteria. Overall, the genetic diversity we observe and how it is distributed throughout the genome and the population should support the resilience of this important pollinator species to ongoing and future selective pressures. Applying our approach to more species should help understand how they can differ in their adaptive potential, and to develop conservation strategies for those most at risk.

Temperature predicts the rate of molecular evolution in Australian Eugongylinae skinks

Temperature differences over time and space have been hypothesized to cause variation in the rate of molecular evolution of species, but empirical evidence is mixed. To further test this hypothesis, we utilized a large exon-capture sequence data of Australian Eugongylinae skinks, exemplifying a radiation of temperature-sensitive ectotherms spanning a large latitudinal gradient. The association between temperature (and other species traits) and long-term substitution rate was assessed based on 1268 sequenced exons of 44 species pairs from the Eugongylinae subfamily using regression analyses. Temperature is the strongest, positively correlated predictor of variation in substitution rate across the Australian Eugongylinae. It explains 45% of variation in synonymous substitution rate, and 11% after controlling for all the other factors. Synonymous substitution rate is also negatively associated with body size, with a 6% variation explained by body size after controlling for the effects of temperature. Other factors are not associated with synonymous substitution rate after controlling for temperature. Overall, this study points to temperature as a strong predictor of the molecular evolution rate in the Eugongylinae subfamily, and demonstrates the power of large-scale exonic data to identify correlates of the rate of molecular evolution.

Molecular evolution of bumble bee vitellogenin and vitellogenin-like genes

Vitellogenin (Vg), a storage protein, has been significantly studied for its egg yolk precursor role in oviparous animals. Recent studies found that vitellogenin and its Vg-like homologs were fundamentally involved in many other biological processes in social insects such as female caste differences and oxidative stress resilience. In this study, we conducted the first large-scale molecular evolutionary analyses of vitellogenin coding genes (Vg) and Vg-like genes of bumble bees, a primitively eusocial insect belonging to the genus Bombus. We obtained sequences for each of the four genes (Vg, Vg-like-A, Vg-like-B, and Vg-like-C) from 27 bumble bee genomes (nine were newly sequenced in this study), and sequences from the two closest clades of Bombus, including five Apis species and five Tetragonula species. Our molecular evolutionary analyses show that in bumble bee, the conventional Vg experienced strong positive selection, while the Vg-like genes showed overall relaxation of purifying selection. In Apis and Tetragonula; however, all four genes were found under purifying selection. Furthermore, the conventional Vg showed signs of strong positive selection in most subgenera in Bombus, apart from the obligate parasitic subgenus Psithyrus which has no caste differentiation. Together, these results indicate that the conventional Vg, a key pleiotropic gene in social insects, is the most rapidly evolving copy, potentially due to its multiple known social functions for both worker and queen castes. This study shows that concerted evolution and purifying selection shaped the evolution of the Vg gene family following their ancient gene duplication and may be the leading forces behind the evolution of new potential protein function enabling functional social pleiotropy.

Genus-Wide Characterization of Bumblebee Genomes Provides Insights into Their Evolution and Variation in Ecological and Behavioral Traits

Bumblebees are a diverse group of globally important pollinators in natural ecosystems and for agricultural food production. With both eusocial and solitary life-cycle phases, and some social parasite species, they are especially interesting models to understand social evolution, behavior, and ecology. Reports of many species in decline point to pathogen transmission, habitat loss, pesticide usage, and global climate change, as interconnected causes. These threats to bumblebee diversity make our reliance on a handful of well-studied species for agricultural pollination particularly precarious. To broadly sample bumblebee genomic and phenotypic diversity, we de novo sequenced and assembled the genomes of 17 species, representing all 15 subgenera, producing the first genus-wide quantification of genetic and genomic variation potentially underlying key ecological and behavioral traits. The species phylogeny resolves subgenera relationships, whereas incomplete lineage sorting likely drives high levels of gene tree discordance. Five chromosome-level assemblies show a stable 18-chromosome karyotype, with major rearrangements creating 25 chromosomes in social parasites. Differential transposable element activity drives changes in genome sizes, with putative domestications of repetitive sequences influencing gene coding and regulatory potential. Dynamically evolving gene families and signatures of positive selection point to genus-wide variation in processes linked to foraging, diet and metabolism, immunity and detoxification, as well as adaptations for life at high altitudes. Our study reveals how bumblebee genes and genomes have evolved across the Bombus phylogeny and identifies variations potentially linked to key ecological and behavioral traits of these important pollinators.

Phylogenomics and evolutionary history of Oreobates (Anura: Craugastoridae) Neotropical frogs along elevational gradients

Mountain ranges offer opportunities for understanding how species evolved and diversified across different environmental conditions. Neotropical frogs of the genus Oreobates (Anura: Craugastoridae) are adapted to highland and lowland habitats along the Andes, but many aspects of their evolution remain unknown. We studied their evolutionary history using ~18,000 exons enriched by targeted sequence-capture. Since capture success was very variable across samples, we evaluated to what degree differing data filtering produced robust inferences. The inferred evolutionary framework evidenced phylogenetic discordances among lowland species that can be explained by taxonomic misidentification or admixture of ancestral lineages. Highland species showed smaller effective populations than lowland frogs, probably due to greater habitat fragmentation in montane environments. Stronger genetic drift likely decreased the power of purifying selection and led to an increased proportion of nonsynonymous mutations in highland populations that could play an important role in their adaptation. Overall, our work sheds light on the evolutionary history and diversification of this group of Neotropical frogs along elevational gradients in the Andes as well as on their patterns of intraspecific diversity.

De Novo Genome Assemblies for Three North American Bumble Bee Species: Bombus bifarius, Bombus vancouverensis, and Bombus vosnesenskii

Bumble bees are ecologically and economically important insect pollinators. Three abundant and widespread species in western North America, Bombus bifarius, Bombus vancouverensis, and Bombus vosnesenskii, have been the focus of substantial research relating to diverse aspects of bumble bee ecology and evolutionary biology. We present de novo genome assemblies for each of the three species using hybrid assembly of Illumina and Oxford Nanopore Technologies sequences. All three assemblies are of high quality with large N50s (> 2.2 Mb), BUSCO scores indicating > 98% complete genes, and annotations producing 13,325 - 13,687 genes, comparing favorably with other bee genomes. Analysis of synteny against the most complete bumble bee genome, Bombus terrestris, reveals a high degree of collinearity. These genomes should provide a valuable resource for addressing questions relating to functional genomics and evolutionary biology in these species.

Specialization of plant-pollinator interactions increases with temperature at Mt. Kilimanjaro

AIM

Species differ in their degree of specialization when interacting with other species, with significant consequences for the function and robustness of ecosystems. In order to better estimate such consequences, we need to improve our understanding of the spatial patterns and drivers of specialization in interaction networks.

METHODS

Here, we used the extensive environmental gradient of Mt. Kilimanjaro (Tanzania, East Africa) to study patterns and drivers of specialization, and robustness of plant-pollinator interactions against simulated species extinction with standardized sampling methods. We studied specialization, network robustness and other network indices of 67 quantitative plant-pollinator networks consisting of 268 observational hours and 4,380 plant-pollinator interactions along a 3.4 km elevational gradient. Using path analysis, we tested whether resource availability, pollinator richness, visitation rates, temperature, and/or area explain average specialization in pollinator communities. We further linked pollinator specialization to different pollinator taxa, and species traits, that is, proboscis length, body size, and species elevational ranges.

RESULTS

We found that specialization decreased with increasing elevation at different levels of biological organization. Among all variables, mean annual temperature was the best predictor of average specialization in pollinator communities. Specialization differed between pollinator taxa, but was not related to pollinator traits. Network robustness against simulated species extinctions of both plants and pollinators was lowest in the most specialized interaction networks, that is, in the lowlands.

CONCLUSIONS

Our study uncovers patterns in plant-pollinator specialization along elevational gradients. Mean annual temperature was closely linked to pollinator specialization. Energetic constraints, caused by short activity timeframes in cold highlands, may force ectothermic species to broaden their dietary spectrum. Alternatively or in addition, accelerated evolutionary rates might facilitate the establishment of specialization under warm climates. Despite the mechanisms behind the patterns have yet to be fully resolved, our data suggest that temperature shifts in the course of climate change may destabilize pollination networks by affecting network architecture.

Mating precedes selective immune priming which is maintained throughout bumblebee queen diapause

Background

Understanding the mechanisms by which organisms adapt to unfavourable conditions is a fundamental question in ecology and evolutionary biology. One such mechanism is diapause, a period of dormancy typically found in nematodes, fish, crustaceans and insects. This state is a key life-history event characterised by arrested development, suppressed metabolism and increased stress tolerance and allows an organism to avoid prolonged periods of harsh and inhospitable environmental conditions. For some species, diapause is preceded by mating which can have a profound effect on female behaviour, physiology and key biological processes, including immunity. However, our understanding of how mating impacts long-term immunity and whether these effects persist throughout diapause is currently limited. To address this, we explored molecular changes in the haemolymph of the ecologically important pollinator, the buff-tailed bumblebee Bombus terrestris. B. terrestris queens mate prior to entering diapause, a non-feeding period of arrested development that can last 6–9 months. Using mass-spectrometry-based proteomics, we quantified changes in the pre-diapause queen haemolymph after mating, as well as the subsequent protein expression of mated queens during and post-diapause.

Results

Our analysis identified distinct proteome profiles associated with diapause preparation, maintenance and termination. More specifically, mating pre-diapause was followed by an increase in the abundance of antimicrobial peptides, key effectors of the immune system. Furthermore, we identified the elevated abundance of these proteins to be maintained throughout diapause. This finding was in contrast to the general reduction observed in immune proteins during diapause suggestive of selective immune priming and expression during diapause. Diapause also affected the expression of proteins involved in cuticular maintenance, olfaction, as well as proteins of unknown function, which may have roles in diapause regulation.

Conclusions

Our results provide clear molecular evidence for the consequences and benefits of mating at the immune level as it precedes the selective increased abundance of antimicrobial peptides that are sustained throughout diapause. In addition, our results provide novel insights into the molecular mechanisms by which bumblebees prepare for, survive, and recover from diapause, insights that may have implications for our general understanding of these processes in other insect groups.