Conservation Genomic Analysis of the Asian Honeybee in China Reveals Climate Factors Underlying Its Population Decline

Pan-genome analysis highlights the role of structural variation in the evolution and environmental adaptation of Asian honeybees

The Asian honeybee, Apis cerana, is an ecologically and economically important pollinator. Mapping its genetic variation is key to understanding population-level health, histories and potential capacities to respond to environmental changes. However, most efforts to date were focused on single nucleotide polymorphisms (SNPs) based on a single reference genome, thereby ignoring larger scale genomic variation. We employed long-read sequencing technologies to generate a chromosome-scale reference genome for the ancestral group of A. cerana. Integrating this with 525 resequencing data sets, we constructed the first pan-genome of A. cerana, encompassing almost the entire gene content. We found that 31.32% of genes in the pan-genome were variably present across populations, providing a broad gene pool for environmental adaptation. We identified and characterized structural variations (SVs) and found that they were not closely linked with SNP distributions; however, the formation of SVs was closely associated with transposable elements. Furthermore, phylogenetic analysis using SVs revealed a novel A. cerana ecological group not recoverable from the SNP data. Performing environmental association analysis identified a total of 44 SVs likely to be associated with environmental adaptation. Verification and analysis of one of these, a 330 bp deletion in the Atpalpha gene, indicated that this SV may promote the cold adaptation of A. cerana by altering gene expression. Taken together, our study demonstrates the feasibility and utility of applying pan-genome approaches to map and explore genetic feature variations of honeybee populations, and in particular to examine the role of SVs in the evolution and environmental adaptation of A. cerana.

Conservation genomics analysis reveals recent population decline and possible causes in bumblebee Bombus opulentus

Bumblebees are a genus of pollinators (Bombus) that play important roles in natural ecosystem and agricultural production. Several bumblebee species have been recorded as under population decline, and the proportion of species experiencing population decline within subgenus Thoracobombus is higher than average. Bombus opulentus is 1 species in Thoracobombus, but little is known about its recent population dynamics. Here, we employed conservation genomics methods to investigate the population dynamics of B. opulentus during the recent past and identify the likely environmental factors that may cause population decline. Firstly, we placed the scaffold-level of B. opulentus reference genome sequence onto chromosome-level using Hi-C technique. Then, based on this reference genome and whole-genome resequencing data for 51 B. opulentus samples, we reconstructed the population structure and effective population size (Ne ) trajectories of B. opulentus and identified genes that were under positive selection. Our results revealed that the collected B. opulentus samples could be divided into 2 populations, and 1 of them experienced a recent population decline; the declining population also exhibited lower genetic diversity and higher inbreeding levels. Genes related to high-temperature tolerance, immune response, and detoxication showed signals of positive selection in the declining population, suggesting that climate warming and pathogen/pesticide exposures may contribute to the decline of this B. opulentus population. Taken together, our study provided insights into the demography of B. opulentus populations and highlighted that populations of the same bumblebee species could have contrasting Ne trajectories and population decline could be caused by a combination of various stressors.

An empirical test of the estimation of historical effective population size using Drosophila melanogaster

The availability of a large number of high-density markers (SNPs) allows the estimation of historical effective population size (Ne ) from linkage disequilibrium between loci. A recent refinement of methods to estimate historical Ne from the recent past has been shown to be rather accurate with simulation data. The method has also been applied to real data for numerous species. However, the simulation data cannot encompass all the complexities of real genomes, and the performance of any estimation method with real data is always uncertain, as the true demography of the populations is not known. Here, we carried out an experimental design with Drosophila melanogaster to test the method with real data following a known demographic history. We used a population maintained in the laboratory with a constant census size of about 2800 individuals and subjected the population to a drastic decline to a size of 100 individuals. After a few generations, the population was expanded back to the previous size and after a few further generations again expanded to twice the initial size. Estimates of historical Ne were obtained with the software GONE both for autosomal and X chromosomes from samples of 17 individuals sequenced for the whole genome. Estimates of the historical effective size were able to infer the patterns of changes that occurred in the populations showing generally good performance of the method. We discuss the limitations of the method and the application of the software carried out so far.