Transcriptomic analysis to uncover genes affecting cold resistance in the Chinese honey bee (Apis cerana cerana)

Drivers of genomic differentiation landscapes in populations of disparate ecological and geographical settings within mainland Apis cerana

Elucidating the evolutionary processes that drive population divergence can enhance our understanding of the early stages of speciation and inform conservation management decisions. The honeybee Apis cerana displays extensive population divergence, providing an informative natural system for exploring these processes. The mainland lineage A. cerana includes several peripheral subspecies with disparate ecological and geographical settings radiated from a central ancestor. Under this evolutionary framework, we can explore the patterns of genome differentiation and the evolutionary models that explain them. We can also elucidate the contribution of non-genomic spatiotemporal mechanisms (extrinsic features) and genomic mechanisms (intrinsic features) that influence these genomic differentiation landscapes. Based on 293 whole genomes, a small part of the genome is highly differentiated between central-peripheral subspecies pairs, while low and partial parallelism partly reflects idiosyncratic responses to environmental differences. Combined elements of recurrent selection and speciation-with-gene-flow models generate the heterogeneous genome landscapes. These elements weight differently between central-island and other central-peripheral subspecies pairs, influenced by glacial cycles superimposed on different geomorphologies. Although local recombination rates exert a significant influence on patterns of genomic differentiation, it is unlikely that low-recombination rates regions were generated by structural variation. In conclusion, complex factors including geographical isolation, divergent ecological selection and non-uniform genome features have acted concertedly in the evolution of reproductive barriers that could reduce gene flow in part of the genome and facilitate the persistence of distinct populations within mainland lineage of A. cerana.

A rapid return to normal: temporal gene expression patterns following cold exposure in the bumble bee Bombus impatiens

Bumble bees are common in cooler climates and many species likely experience periodic exposure to very cold temperatures, but little is known about the temporal dynamics of cold response mechanisms following chill exposure, especially how persistent effects of cold exposure may facilitate tolerance of future events. To investigate molecular processes involved in the temporal response by bumble bees to acute cold exposure, we compared mRNA transcript abundance in Bombus impatiens workers exposed to 0°C for 75 min (inducing chill coma) and control bees maintained at a constant ambient temperature (28°C). We sequenced the 3' end of mRNA transcripts (TagSeq) to quantify gene expression in thoracic tissue of bees at several time points (0, 10, 30, 120 and 720 min) following cold exposure. Significant differences from control bees were only detectable within 30 min after the treatment, with most occurring at the 10 min recovery time point. Genes associated with gluconeogenesis and glycolysis were most notably upregulated, while genes related to lipid and purine metabolism were downregulated. The observed patterns of expression indicate a rapid recovery after chill coma, suggesting an acute differential transcriptional response during recovery from chill coma and return to baseline expression levels within an hour, with no long-term gene expression markers of this cold exposure. Our work highlights the functions and pathways important for acute cold recovery, provides an estimated time frame for recovery from cold exposure in bumble bees, and suggests that cold hardening may be less important for these heterothermic insects.

Genomic selection pressure discovery using site-frequency spectrum and reduced local variability statistics in Pakistani Dera-Din-Panah goat

BACKGROUND

Population geneticists have long sought to comprehend various selection traces accumulated in the goat genome due to natural or human driven artificial selection through breeding practices, which led the wild animals to domestication, so understanding evolutionary process may helpful to utilize the full genetic potential of goat genome.

METHODS AND RESULTS

As a step forward to pinpoint the selection signals in Pakistani Dera-Din-Panah (DDP) goat, whole-genome pooled sequencing (n = 12) was performed, and 618,236,192 clean paired-end reads were mapped against ARS1 reference goat assembly. Five different selection signature statistics were applied using four site-frequency spectrum (SFS) methods (Tajima's D ([Formula: see text]), Fay and Wu's H ([Formula: see text]), Zeng's E ([Formula: see text]), [Formula: see text]) and one reduced local variability approach named pooled heterozygosity ([Formula: see text]). The under-selection regions were annotated with significant threshold values of [Formula: see text]≥4.7, [Formula: see text]≥6, [Formula: see text]≥2.5, Pool-HMM ≥ 12, and [Formula: see text]≥5 that resulted in accumulative 364 candidate gene hits. The highest genomic selection signals were observed on Chr. 4, 6, 10, 12, 15, 16, 18, 20, and 27 and harbor ADAMTS6, CWC27, RELN, MYCBP2, FGF14, STIM1, CFAP74, GNB1, CALML6, TMEM52, FAM149A, NADK, MMP23B, OPN3, FH, MFHAS1, KLKB1, RRM1, KMO, SPEF2, F11, KIT, KMO, ERI1, ATP8B4, and RHOG genes. Next, the validation of our captured genomic hits was also performed by more than one applied statistics which harbor meat production, immunity, and reproduction associated genes to strengthen our hypothesis of under-selection traits in this Pakistani goat breed. Furthermore, common candidate genes captured by more than one statistical method were subjected to gene ontology and KEGG pathway analysis to get insights of particular biological processes associated with this goat breed.

CONCLUSION

Current perception of genomic architecture of DDP goat provides a better understanding to improve its genetic potential and other economically important traits of medium to large body size, milk, and fiber production by updating the genomic insight driven breeding strategies to boost the livestock and agriculture-based economy of the country.

The impact of geography and climate on the population structure and local adaptation in a wild bee

Deciphering processes that contribute to genetic differentiation and divergent selection of natural populations is useful for evaluating the adaptive potential and resilience of organisms faced with various anthropogenic stressors. Insect pollinator species, including wild bees, provide critical ecosystem services but are highly susceptible to biodiversity declines. Here, we use population genomics to infer the genetic structure and test for evidence of local adaptation in an economically important native pollinator, the small carpenter bee (Ceratina calcarata). Using genome-wide SNP data (n = 8302), collected from specimens across the species' entire distribution, we evaluated population differentiation and genetic diversity and identified putative signatures of selection in the context of geographic and environmental variation. Results of the analyses of principal component and Bayesian clustering were concordant with the presence of two to three genetic clusters, associated with landscape features and inferred phylogeography of the species. All populations examined in our study demonstrated a heterozygote deficit, along with significant levels of inbreeding. We identified 250 robust outlier SNPs, corresponding to 85 annotated genes with known functional relevance to thermoregulation, photoperiod, and responses to various abiotic and biotic stressors. Taken together, these data provide evidence for local adaptation in a wild bee and highlight genetic responses of native pollinators to landscape and climate features.

A key gene for the climatic adaptation of Apis cerana populations in China according to selective sweep analysis

BACKGROUND

Apis cerana is widely distributed in China and, prior to the introduction of western honeybees, was the only bee species kept in China. During the long-term natural evolutionary process, many unique phenotypic variations have occurred among A. cerana populations in different geographical regions under varied climates. Understanding the molecular genetic basis and the effects of climate change on the adaptive evolution of A. cerana can promote A. cerana conservation in face of climate change and allow for the effective utilization of its genetic resources.

RESULT

To investigate the genetic basis of phenotypic variations and the impact of climate change on adaptive evolution, A. cerana workers from 100 colonies located at similar geographical latitudes or longitudes were analyzed. Our results revealed an important relationship between climate types and the genetic variation of A. cerana in China, and a greater influence of latitude compared with longitude was observed. Upon selection and morphometry analyses combination for populations under different climate types, we identified a key gene RAPTOR, which was deeply involved in developmental processes and influenced the body size.

CONCLUSION

The selection of RAPTOR at the genomic level during adaptive evolution could allow A. cerana to actively regulate its metabolism, thereby fine-tuning body sizes in response to harsh conditions caused by climate change, such as food shortages and extreme temperatures, which may partially elucidate the size differences of A. cerana populations. This study provides crucial support for the molecular genetic basis of the expansion and evolution of naturally distributed honeybee populations.

Molecular Mechanisms of Winter Survival

Winter provides many challenges for insects, including direct injury to tissues and energy drain due to low food availability. As a result, the geographic distribution of many species is tightly coupled to their ability to survive winter. In this review, we summarize molecular processes associated with winter survival, with a particular focus on coping with cold injury and energetic challenges. Anticipatory processes such as cold acclimation and diapause cause wholesale transcriptional reorganization that increases cold resistance and promotes cryoprotectant production and energy storage. Molecular responses to low temperature are also dynamic and include signaling events during and after a cold stressor to prevent and repair cold injury. In addition, we highlight mechanisms that are subject to selection as insects evolve to variable winter conditions. Based on current knowledge, despite common threads, molecular mechanisms of winter survival vary considerably across species, and taxonomic biases must be addressed to fully appreciate the mechanistic basis of winter survival across the insect phylogeny.

Full length transcriptomes analysis of cold-resistance of Apis cerana in Changbai Mountain during overwintering period

Apis cerana in Changbai Mountain is an ecological type of Apis cerana, which is an excellent breeding material with cold-resistant developed by long-term natural selection under the ecological conditions. However, the physiological and molecular mechanisms of Changbai Mountain population under cold stress are still unclear. In this study, the Nanopore sequencing was carried out for the transcriptome of Apis cerana in Changbai Mountain in the coldest period of overwintering, which will provide a reference to the cold-resistant mechanism. We determined 5,941 complete ORF sequences, 1,193 lncRNAs, 619 TFs, 10,866 SSRs and functional annotations of 11,599 new transcripts. Our results showed that the myosin family and the C2H2 zinc finger protein transcription factor family possibly have significant impacts on the response mechanism of cold stress during overwintering. In addition, the cold environment alters genes expression profiles in honeybees via different AS and APA mechanisms. These altered genes in Hippo, Foxo, and MARK pathways help them counter the stress of cold in overwinter period. Our results might provide clues about the response of eastern honeybees to extreme cold, and reflect the possible genetic basis of physiological changes.

Identification and Functional Characterization of Antifreeze Protein and Its Mutants in Dendroctonus armandi (Coleoptera: Curculionidae: Scolytinae) Larvae Under Cold Stress

Dendroctonus armandi (Tsai and Li) (Coleoptera: Curculionidae: Scolytinae) is considered to be the most destructive forest pest in the Qinling and Bashan Mountains of China. Low winter temperatures limit insect's populations, distribution, activity, and development. Insects have developed different strategies such as freeze-tolerance and freeze-avoidance to survive in low temperature conditions. In the present study, we used gene cloning, real-time polymerase chain reaction (PCR), RNA interference (RNAi), and heterologous expression to study the function of the D. armandi antifreeze protein gene (DaAFP). We cloned the 800 bp full-length cDNA encoding 228 amino acids of DaAFP and analyzed its structure using bioinformatics analysis. The DaAFP amino acid sequence exhibited 24-86% similarity with other insect species. The expression of DaAFP was high in January and in the larvae, head, and midgut of D. armandi. In addition, the expression of DaAFP increased with decreasing temperature and increasing exposure time. RNAi analysis also demonstrated that AFP plays an important role in the cold tolerance of overwintering larvae. The thermal hysteresis and antifreeze activity assay of DaAFP and its mutants indicated that the more regular the DaAFP threonine-cystine-threonine (TXT) motif, the stronger the antifreeze activity. These results suggest that DaAFP plays an essential role as a biological cryoprotectant in overwintering D. armandi larvae and provides a theoretical basis for new pest control methods.

Geometric morphology and population genomics provide insights into the adaptive evolution of Apis cerana in Changbai Mountain

Background

Exploration of adaptive evolution of organisms in response to environmental change can help to understand the evolutionary history of species and the underlying mechanisms of adaptation to local environments, thus guiding future conservation programmes. Before the introduction of Apis mellifera in China, eastern honey bees (Apis cerana) were the only species used for beekeeping in this region. In the mountains of Changbai, populations of A. cerana are considered a distinct ecotype of the species which formed through the distinct selective pressures in this area over time.

Result

We performed a measure of 300 wing specimens of eastern honey bees and obtained the geometric morphological variation in the wing of A. cerana in Changbai Mountain. A total of 3,859,573 high-quality SNP loci were yielded via the whole-genome resequencing of 130 individuals in 5 geographic regions.

Conclusion

Corresponding geometric morphology and population genomics confirmed the particularity of the A. cerana in Changbai Mountain. Genetic differentiation at the subspecies level exists between populations in Changbai Mountain and remaining geographic regions, and a significant reduction in the effective population size and an excessive degree of inbreeding may be responsible for a substantial loss of population genetic diversity. Candidate genes potentially associated with cold environmental adaptations in populations under natural selection were identified, which may represent local adaptations in populations. Our study provided insights into the evolutionary history and adaptation of A. cerana in Changbai Mountain, as well as its conservation.

Identification of quantitative trait loci associated with calmness and gentleness in honey bees using whole-genome sequences

The identification of quantitative trait loci (QTL) through genome-wide association studies (GWAS) is a powerful method for unravelling the genetic background of selected traits and improving early-stage predictions. In honey bees (Apis mellifera), past genetic analyses have particularly focused on individual queens and workers. In this study, we used pooled whole-genome sequences to ascertain the genetic variation of the entire colony. In total, we sampled 216 Apis mellifera mellifera and 28 Apis mellifera carnica colonies. Different experts subjectively assessed the gentleness and calmness of the colonies using a standardised protocol. Conducting a GWAS for calmness on 211 purebred A. m. mellifera colonies, we identified three QTL, on chromosomes 8, 6, and 12. The two first QTL correspond to LOC409692 gene, coding for a disintegrin and metalloproteinase domain-containing protein 10, and to Abscam gene, coding for a Dscam family member Abscam protein, respectively. The last gene has been reported to be involved in the domestication of A. mellifera. The third QTL is located 13 kb upstream of LOC102655631, coding for a trehalose transporter. For gentleness, two QTL were identified on chromosomes 4 and 3. They are located within gene LOC413669, coding for a lap4 protein, and gene LOC413416, coding for a bicaudal C homolog 1-B protein, respectively. The identified positional candidate genes of both traits mainly affect the olfaction and nervous system of honey bees. Further research is needed to confirm the results and to better understand the genetic and phenotypic basis of calmness and gentleness.

Expression analysis of genes related to cold tolerance in Dendroctonus valens

Pine beetles are well known in North America for their widespread devastation of pine forests. However, Dendroctonus valens LeConte is an important invasive forest pest in China also. Adults and larvae of this bark beetle mainly winter at the trunks and roots of Pinus tabuliformis and Pinus sylvestris; larvae, in particular, result in pine weakness or even death. Since the species was introduced from the United States to Shanxi in 1998, its distribution has spread northward. In 2017, it invaded a large area at the junction of Liaoning, Inner Mongolia and Hebei provinces, showing strong cold tolerance. To identify genes relevant to cold tolerance and the process of overwintering, we sequenced the transcriptomes of wintering and non-wintering adult and larval D. valens using the Illumina HiSeq platform. Differential expression analysis methods for other non-model organisms were used to compare transcript abundances in adults and larvae at two time periods, followed by the identification of functions and metabolic pathways related to genes associated with cold tolerance. We detected 4,387 and 6,091 differentially expressed genes (DEGs) between sampling dates in larvae and adults, respectively, and 1,140 common DEGs, including genes encoding protein phosphatase, very long-chain fatty acids protein, cytochrome P450, and putative leucine-rich repeat-containing proteins. In a Gene Ontology (GO) enrichment analysis, 1,140 genes were assigned to 44 terms, with significant enrichment for cellulase activity, hydrolase activity, and carbohydrate metabolism. Kyoto Encyclopedia of Genes and Genomes (KEGG) classification and enrichment analyses showed that the lysosomal and purine metabolism pathways involved the most DEGs, the highly enriched terms included autophagy-animal, pentose and glucuronate interconversions and lysosomal processes. We identified 140 candidate genes associated with cold tolerance, including genes with established roles in this trait (e.g., genes encoding trehalose transporter, fructose-1,6-bisphosphatase, and trehalase). Our comparative transcriptome analysis of adult and larval D. valens in different conditions provides basic data for the discovery of key genes and molecular mechanisms underlying cold tolerance.

The role of melatonin and Tryptophan-5-hydroxylase-1 in different abiotic stressors in Apis cerana cerana

Tryptophan-5-hydroxylase-1 (T5H-1) is the rate-limiting enzyme in the biosynthesis of serotonin, which is involved in the biosynthesis of melatonin (Mel). Mel, a biological hormone, plays crucial roles in stressors tolerance, such as cold, hot, Ultraviolet (UV) and pesticide tolerance. However, the direct correlation between T5H-1 and Mel and the underlying mechanism in organisms remains elusive. Mel-mediated cold tolerance was studied extensively in plants and somewhat in insects, including bees. The present study isolated the Mel synthesis gene T5H-1 from Apis cerana cerana for the first time. qRT-PCR analysis indicated that AccT5H-1 played vital roles during some adverse conditions, including 4 °C, 8 °C, 10 °C, 45 °C, UV, cyhalothrin, abamectin, paraquat and bifenthrin exposure. Knockdown of AccT5H-1 using RNA interference (RNAi) technology upregulated most antioxidant genes. Additionally, an enzyme activity assay revealed higher contents of Malondialdehyde (MDA) and Hydrogen peroxide (H₂O₂), lower content of Vitamin C (VC), and higher activities of Glutathione S-transferase (GST), Superoxide dismutase (SOD), Catalase (CAT) and Peroxidase (POD) in the AccT5H-1 silenced group than the control group. These results suggest that AccT5H-1 is involved in the response to different oxidative stressors in A. cerana cerana. The survival rate of A. cerana cerana exposed to low temperature treatment revealed that the optimal concentration of Mel in the diet was 10 µg/mL. We also found that the antioxidant enzyme (GST, SOD, POD and CAT) concentrations at 10 µg/mL Mel increased to different degrees, and the content of oxidizing substances (MDA and H₂O₂) decreased, the content of VC increased, and the content of substances that promote cold resistance (glycerol and glycogen) increased. Mel increased the resistance of A. cerana cerana exposed to low temperatures. The expression of AccT5H-1 decreased after the feeding of exogenous Mel to bees. These results provide a reference for other insect studies on Mel and T5H-1.

Genomic analysis of Asian honeybee populations in China reveals evolutionary relationships and adaptation to abiotic stress

The geographic and biological diversity of China has resulted in the differential adaptation of the eastern honeybee, Apis cerana, to these varied habitats. A. cerana were collected from 14 locations in China. Their genomes were sequenced, and nucleotide polymorphisms were identified at more than 9 million sites. Both STRUCTURE and principal component analysis placed the bees into seven groups. Phylogenomic analysis groups the honeybees into many of the same clusters with high bootstrap values (91%-100%). Populations from Tibet and South Yunnan are sister taxa and together represent the earliest diverging lineage included in this study. We propose that the evolutionary origin of A. cerana in China was in the southern region of Yunnan Province and expanded from there into the southeastern regions and into the northeastern mountain regions. The Cold-Temperate West Sichuan Plateau and Tropical Diannan populations were compared to identify genes under adaptive selection in these two habitats. Pathway enrichment analysis showing genes under selection, including the Hippo signaling pathway, GABAergic pathway, and trehalose-phosphate synthase, indicates that most genes under selection pressure are involved in the process of signal transduction and energy metabolism. qRT-PCR analysis reveals that one gene under selection, the AcVIAAT gene, involved in the GABAergic pathway, is responding to cold temperature stress. Through homologous recombination, we show that the AcVIAAT gene is able to replace the CNAG_01904 gene in the fungus Cryptococcus neoformans and that it makes the fungus less sensitive to conditions of oxidative stress and variations in temperature. Our results contribute to our understanding of the evolutionary origin of A. cerana in China and the molecular basis of environmental adaptation.

Adaptations to thermal stress in social insects: recent advances and future directions

Thermal stress is a major driver of population declines and extinctions. Shifts in thermal regimes create new environmental conditions, leading to trait adaptation, population migration, and/or species extinction. Extensive research has examined thermal adaptations in terrestrial arthropods. However, little is known about social insects, despite their major role in ecosystems. It is only within the last few years that the adaptations of social insects to thermal stress have received attention. Herein, we discuss what is currently known about thermal tolerance and thermal adaptation in social insects - namely ants, termites, social bees, and social wasps. We describe the behavioural, morphological, physiological, and molecular adaptations that social insects have evolved to cope with thermal stress. We examine individual and collective responses to both temporary and persistent changes in thermal conditions and explore the extent to which individuals can exploit genetic variability to acclimatise. Finally, we consider the costs and benefits of sociality in the face of thermal stress, and we propose some future research directions that should advance our knowledge of individual and collective thermal adaptations in social insects.

Transcriptomic analysis to elucidate the response of honeybees (Hymenoptera: Apidae) to amitraz treatment

Amitraz is an acaricide that is widely used in apiculture. Several studies have reported that in honeybees (Apis mellifera Linnaeus; Hymenoptera: Apidae), amitraz affects learning, memory, behavior, immunity, and various other physiological processes. Despite this, few studies have explored the molecular mechanisms underlying the action of amitraz on honeybees. Here, we investigated the transcriptome of honeybees after exposure to 9.4 mg/L amitraz for 10 d, a subchronic dose. Overall, 279 differentially expressed genes (DEGs) were identified (237 upregulated, 42 downregulated). Several, including Pla2, LOC725381, LOC413324, LOC724386, LOC100577456, LOC551785, and P4504c3, were validated by quantitative PCR. According to gene ontology, DEGs were mainly involved in metabolism, biosynthesis, and translation. Kyoto Encyclopedia of Genes and Genomes pathway analyses revealed that amitraz treatment affected the relaxin signaling pathway, platelet activation, and protein digestion and absorption.

Within species expressed genetic variability and gene expression response to different temperatures in the rotifer Brachionus calyciflorus sensu stricto

Genetic divergence is impacted by many factors, including phylogenetic history, gene flow, genetic drift, and divergent selection. Rotifers are an important component of aquatic ecosystems, and genetic variation is essential to their ongoing adaptive diversification and local adaptation. In addition to coding sequence divergence, variation in gene expression may relate to variable heat tolerance, and can impose ecological barriers within species. Temperature plays a significant role in aquatic ecosystems by affecting species abundance, spatio-temporal distribution, and habitat colonization. Recently described (formerly cryptic) species of the Brachionus calyciflorus complex exhibit different temperature tolerance both in natural and in laboratory studies, and show that B. calyciflorus sensu stricto (s.s.) is a thermotolerant species. Even within B. calyciflorus s.s., there is a tendency for further temperature specializations. Comparison of expressed genes allows us to assess the impact of stressors on both expression and sequence divergence among disparate populations within a single species. Here, we have used RNA-seq to explore expressed genetic diversity in B. calyciflorus s.s. in two mitochondrial DNA lineages with different phylogenetic histories and differences in thermotolerance. We identify a suite of candidate genes that may underlie local adaptation, with a particular focus on the response to sustained high or low temperatures. We do not find adaptive divergence in established candidate genes for thermal adaptation. Rather, we detect divergent selection among our two lineages in genes related to metabolism (lipid metabolism, metabolism of xenobiotics).

Long-Distance Transportation Causes Temperature Stress in the Honey Bee, Apis mellifera (Hymenoptera: Apidae)

Pollination services provided by the honey bee, Apis mellifera (Hymenoptera: Apidae, Linnaeus, 1758) have broad economic impacts and are necessary for production of a diversity of important crops. Hives may be transported multiple times per year to provide pollination. To test how temperature may contribute to transportation stress, temperature sensors were placed in hives in different locations and orientations on the trailer during shipping. Colony size prior to shipping significantly contributed to loss of population immediately after shipping which contributed to colony failure with smaller colonies more likely to fail and fail faster. Colony size also affects thermoregulation and temperature stress. Internal hive temperature varies significantly based on location and orientation. While colonies near the front and rear of the trailer and those oriented toward the center aisle had significantly different average internal temperatures, colony size best predicts loss of thermoregulation. Additionally, we profiled gene expression at departure, on arrival, and after a recovery period to identify transcriptional responses to transportation. Functional and enrichment analysis identified increased methylation and decreased ribosomal and protein-folding activity. Pheromone and odorant-binding transcripts were up-regulated after transportation. After recovery, transcripts associated with defense response, immune activity, and heat shock decreased, while production of antibiotic peptides increased. We conclude that hives experience considerable temperature stress possibly caused by turbulent airflow in exposed locations. Transportation stress should be considered an important component of annual colony losses which can be mitigated with improved management strategies.

Proteomic Responses Under Cold Stress Reveal Unique Cold Tolerance Mechanisms in the Pacific White Shrimp (Litopenaeus vannamei)

The Pacific white shrimp (Litopenaeus vannamei), one of the most widely cultured shrimp species in the world, often suffers from cold stress. To understand the molecular mechanism of cold tolerance in Pacific white shrimp, we conducted a proteomic analysis on two contrasting shrimp cultivars, namely, cold-tolerant Guihai2 (GH2) and cold-sensitive Guihai1 (GH1), under normal temperature (28°C), under cold stress (16°C), and during recovery to 28°C. In total, 3,349 proteins were identified, among which 2,736 proteins were quantified. Based on gene ontology annotations, differentially expressed proteins largely belonged to biological processes, cellular components, and molecular functions. KEGG pathway annotations indicated that the main changes were observed in the lysosome, ribosomes, and oxidative phosphorylation. Subcellular localization analysis showed a significant increase in proteins present in cytosol, extracellular regions, and mitochondria. Combining enrichment-based clustering analysis and qRT-PCR analysis, we found that glutathione S-transferase, zinc proteinase, m7GpppX diphosphatase, AP2 transcription complex, and zinc-finger transcription factors played a major role in the cold stress response in Pacific white shrimp. Moreover, structure proteins, including different types of lectin and DAPPUDRAFT, were indispensable for cold stress tolerance of the Pacific white shrimp. Results indicate the molecular mechanisms of the Pacific white shrimp in response to cold stress and provide new insight into breeding new cultivars with increased cold tolerance.

Sequencing and Expression Characterization of Antifreeze Protein Maxi-Like in Apis cerana cerana

Antifreeze proteins (AFPs) are biological cryoprotectants with unique properties that play a crucial role in regulating the molecular mechanisms governing cold resistance in insects. To identify and characterize Apis cerana cerana AFP (AcerAFP), we cloned the full-length cDNA of AcerAFP and examined its expression patterns in A. cerana cerana. A nucleotide alignment analysis showed that the entire coding region of AcerAFP is 1095 bp and encodes a polypeptide of 365 amino acids. The amino acid sequence of this protein exhibits 63–96% homology with AFP homologs from other hymenopterans. α-helices form the main secondary and tertiary structures of AcerAFP, which is similar to the molecular structure of fish AFP type-I. The expression profiles of AcerAFP revealed that expression was tissue, sex, and developmentally specific. In response to cold stress, the mRNA and protein expression of AcerAFP were both induced by low temperatures, and were also related to the concentrations of several cryoprotectants, including glucose, glycerin, glutamic acid, cysteine, histidine, alanine, and methionine. In addition, we found that the knockdown of AcerAFP by RNA interference remarkably increased the total freezing temperature of hemolymph in A. cerana cerana, where levels of AcerAFP mRNA were correlated with the expression of most antifreeze-related proteins. Taken together, these results suggest that AcerAFP plays an essential role as a biological cryoprotectant in honeybees, and is in turn regulated by small cryoprotectants and antifreeze-related proteins.