Genetic differentiation and local adaptation of the Japanese honeybee, Apis cerana japonica.

We examine the population genetic structure and divergence among the regional populations of the Japanese honeybee, Apis cerana japonica, by re-sequencing the genomes of 105 individuals from the three main Japanese islands with diverse climates. The genetic structure results indicated that these individuals are distinct from the mainland Chinese A. cerana samples. Furthermore, population structure analyses have identified three genetically distinct geographic regions in Japan: Northern (Tohoku-Kanto-Chubu districts), Central (Chugoku district), and Southern (Kyushu district). In some districts, "possible non-native" individuals, likely introduced from other regions in recent years, were discovered. Then, genome-wide scans were conducted to detect candidate genes for adaptation by two different approaches. We performed a population branch statistics (PBS) analysis to identify candidate genes for population-specific divergence. A latent factor mixed model (LFMM) was used to identify genes associated with climatic variables along a geographic gradient. The PBSmax analysis identified 25 candidate genes for population-specific divergence whereas the LFMM analysis identified 73 candidate genes for adaptation to climatic variables along a geographic gradient. However, no common genes were identified by both methods.

Does selection occur at the intermediate zone of two insufficiently isolated populations? A whole-genome analysis along an altitudinal gradient.

Adaptive divergence occurs even between insufficiently isolated populations when there is a great difference in environments between their habitats. Individuals present in an intermediate zone of the two divergent populations are expected to have an admixed genetic structure due to gene flow. A selective pressure that acts on the genetically admixed individuals may limit the gene flow and maintain the adaptive divergence. Here, we addressed a question whether selection occurs in the genetically admixed individuals between two divergent populations. Arabidopsis halleri is a perennial montane plant, which has clear phenotypic dimorphisms between highland and lowland habitats in Mt. Ibuki, central Japan. We obtained the whole-genome sequences of Arabidopsis halleri plants along an altitudinal gradient of 359-1,317 m with a high spatial resolution (mean altitudinal interval of 20 m). We found a zone where the highland and lowland genes were mixing (intermediate subpopulation). In the intermediate subpopulation, we identified 5 and 13 genome regions, which included 3 and 8 genes, that had a high frequency of alleles that are accumulated in highland and lowland subpopulations, respectively. In addition, we also found that the frequency of highland alleles of these selected genome regions was smaller in the lowland subpopulation compared with that of the non-selected regions. These results suggest that the selection in the intermediate subpopulation might limit the gene flow and contribute to the adaptive divergence between altitudes. We also identified 7 genome regions that had low heterozygote frequencies in the intermediate subpopulation. We conclude that different types of selection in addition to gene flow occur at the intermediate altitude and shape the genetic structure across altitudes.

Meta-Analysis of the Public RNA-Seq Data of the Western Honeybee Apis mellifera to Construct Reference Transcriptome Data.

The Western honeybee (Apis mellifera) is valuable in biological research and agriculture. Its genome sequence was published before those for other insect species. RNA-Seq data for A. mellifera have been applied in several recently published studies. Nevertheless, these data have not been prepared for use in subsequent meta-analyses. To promote A. mellifera transcriptome analysis, we constructed reference transcriptome data using the reference genome sequence and RNA-Seq data curated from about 1,000 runs of public databases. The new reference transcriptome data construct comprised 149,685 transcripts, and 194,174 protein sequences were predicted. Approximately 50-60% of the predicted protein sequences were functionally annotated using the protein sequence data for several model and insect species. Novel candidate immune-related transcripts were searched by meta-analysis using immune-response-related RNA-Seq and reference transcriptome data. Three to twenty candidate transcripts including autophagy-related protein 3 were upregulated or downregulated in response to both viral and bacterial infections. The constructed reference transcriptome data may facilitate future transcriptome analyses of A. mellifera.

The complete mitochondrial genome of the cavity-nesting honeybee, Apis koschevnikovi (Insecta: Hymenoptera: Apidae).

In this study, we analyzed the complete mitochondrial genome of the cavity-nesting honeybee, A. koschevnikovi. The mitochondrial genome of A. koschevnikovi was observed to be a circular molecule of 15,278 bp and was similar to that of the other cavity-nesting honeybee species. The average AT content in the A. koschevnikovi mitochondrial genome was 84%. It was predicted to contain 13 protein-coding, 24 tRNA and two rRNA genes, along with one A + T-rich control region, besides three tRNA-Met repeats.

The complete mitochondrial genome of the Japanese honeybee, Apis cerana japonica (Insecta: Hymenoptera: Apidae).

In this study, we analyzed the complete mitochondrial genome of the Japanese honeybee Apis cerana japonica. The mitochondrial genome of A. c. japonica is a circular molecule of 15 917 bp and is similar to that of A. c. cerana. It contains 13 protein-coding genes, 22 tRNA genes, two rRNA genes and one A + T-rich control region. All protein-coding genes are initiated by ATT and ATG codons and are terminated by the typical stop codon TAA or TAG, except for the start codon of ATP8 which ends with C. All tRNA genes typically form a cloverleaf secondary structure, except for tRNA-Ser (AGN).

Complete mitochondrial genome of the Japanese bumblebee, Bombus hypocrita sapporensis (Insecta: Hymenoptera: Apidae).

In the present study, we describe the complete mitochondrial genome of the bumblebee, Bombus hypocrita sapporensis from the Rebun Island, in Hokkaido, Japan. The mitochondrial genome of B. hypocrita sapporensis includes a circular molecule of 15 700 bp. It contains 13 protein-coding genes, 22 tRNA genes, two rDNA genes and an A + T-rich control region. All protein-coding genes are initiated by ATA, ATG, and ATT codons and are terminated by the typical stop codon TAA or T, except for ND4L, which ends with TA. All tRNA genes typically form a cloverleaf secondary structure.