<?xml version="1.0" encoding="UTF-8" standalone="no"?>
<!DOCTYPE article PUBLIC "-//NLM//DTD Journal Publishing DTD v2.3 20070202//EN" "journalpublishing.dtd">
<article xml:lang="EN" xmlns:mml="http://www.w3.org/1998/Math/MathML" xmlns:xlink="http://www.w3.org/1999/xlink" article-type="research-article">
<front>
<journal-meta>
<journal-id journal-id-type="publisher-id">Front. Microbiol.</journal-id>
<journal-title>Frontiers in Microbiology</journal-title>
<abbrev-journal-title abbrev-type="pubmed">Front. Microbiol.</abbrev-journal-title>
<issn pub-type="epub">1664-302X</issn>
<publisher>
<publisher-name>Frontiers Media S.A.</publisher-name>
</publisher>
</journal-meta>
<article-meta>
<article-id pub-id-type="doi">10.3389/fmicb.2023.1218560</article-id>
<article-categories>
<subj-group subj-group-type="heading">
<subject>Microbiology</subject>
<subj-group>
<subject>Hypothesis and Theory</subject>
</subj-group>
</subj-group>
</article-categories>
<title-group>
<article-title>Potential role of the gut microbiota of bumblebee <italic>Bombus pyrosoma</italic> in adaptation to high-altitude habitats</article-title>
</title-group>
<contrib-group>
<contrib contrib-type="author">
<name><surname>Zhang</surname> <given-names>Zhengyi</given-names></name>
<xref ref-type="aff" rid="aff1"><sup>1</sup></xref>
<uri xlink:href="http://loop.frontiersin.org/people/2289093/overview"/>
</contrib>
<contrib contrib-type="author">
<name><surname>Guo</surname> <given-names>Yulong</given-names></name>
<xref ref-type="aff" rid="aff1"><sup>1</sup></xref>
<uri xlink:href="http://loop.frontiersin.org/people/1107944/overview"/>
</contrib>
<contrib contrib-type="author">
<name><surname>Zhuang</surname> <given-names>Mingsheng</given-names></name>
<xref ref-type="aff" rid="aff1"><sup>1</sup></xref>
<xref ref-type="aff" rid="aff2"><sup>2</sup></xref>
</contrib>
<contrib contrib-type="author">
<name><surname>Liu</surname> <given-names>Fugang</given-names></name>
<xref ref-type="aff" rid="aff1"><sup>1</sup></xref>
</contrib>
<contrib contrib-type="author">
<name><surname>Xia</surname> <given-names>Zhongyan</given-names></name>
<xref ref-type="aff" rid="aff1"><sup>1</sup></xref>
</contrib>
<contrib contrib-type="author">
<name><surname>Zhang</surname> <given-names>Zhihao</given-names></name>
<xref ref-type="aff" rid="aff1"><sup>1</sup></xref>
<uri xlink:href="http://loop.frontiersin.org/people/1592173/overview"/>
</contrib>
<contrib contrib-type="author">
<name><surname>Yang</surname> <given-names>Fan</given-names></name>
<xref ref-type="aff" rid="aff1"><sup>1</sup></xref>
<uri xlink:href="http://loop.frontiersin.org/people/883635/overview"/>
</contrib>
<contrib contrib-type="author">
<name><surname>Zeng</surname> <given-names>Huayan</given-names></name>
<xref ref-type="aff" rid="aff3"><sup>3</sup></xref>
</contrib>
<contrib contrib-type="author">
<name><surname>Wu</surname> <given-names>Yueguo</given-names></name>
<xref ref-type="aff" rid="aff3"><sup>3</sup></xref>
</contrib>
<contrib contrib-type="author" corresp="yes">
<name><surname>Huang</surname> <given-names>Jiaxing</given-names></name>
<xref ref-type="aff" rid="aff1"><sup>1</sup></xref>
<xref ref-type="corresp" rid="c001"><sup>&#x0002A;</sup></xref>
<uri xlink:href="http://loop.frontiersin.org/people/1233741/overview"/>
</contrib>
<contrib contrib-type="author" corresp="yes">
<name><surname>Li</surname> <given-names>Jilian</given-names></name>
<xref ref-type="aff" rid="aff1"><sup>1</sup></xref>
<xref ref-type="corresp" rid="c002"><sup>&#x0002A;</sup></xref>
<uri xlink:href="http://loop.frontiersin.org/people/1110771/overview"/>
</contrib>
</contrib-group>
<aff id="aff1"><sup>1</sup><institution>State Key Laboratory of Resource Insects, Institute of Apicultural Research, Chinese Academy of Agricultural Science</institution>, <addr-line>Beijing</addr-line>, <country>China</country></aff>
<aff id="aff2"><sup>2</sup><institution>Shanghai Suosheng Biotechnology Co., Ltd.</institution>, <addr-line>Shanghai</addr-line>, <country>China</country></aff>
<aff id="aff3"><sup>3</sup><institution>Luoping Yunling Bee Industry and Trade Co., Ltd., Qujing</institution>, <addr-line>Yunnan</addr-line>, <country>China</country></aff>
<author-notes>
<fn fn-type="edited-by"><p>Edited by: Kashif Ur Rehman, Islamia University of Bahawalpur, Pakistan</p></fn>
<fn fn-type="edited-by"><p>Reviewed by: Nazish Roy, Forman Christian College, Pakistan; Raees Khan, National University of Medical Sciences (NUMS), Pakistan</p></fn>
<corresp id="c001">&#x0002A;Correspondence: Jiaxing Huang <email>huangjiaxing&#x00040;caas.cn</email></corresp>
<corresp id="c002">Jilian Li <email>bumblebeeljl&#x00040;hotmail.com</email></corresp>
</author-notes>
<pub-date pub-type="epub">
<day>02</day>
<month>08</month>
<year>2023</year>
</pub-date>
<pub-date pub-type="collection">
<year>2023</year>
</pub-date>
<volume>14</volume>
<elocation-id>1218560</elocation-id>
<history>
<date date-type="received">
<day>07</day>
<month>05</month>
<year>2023</year>
</date>
<date date-type="accepted">
<day>11</day>
<month>07</month>
<year>2023</year>
</date>
</history>
<permissions>
<copyright-statement>Copyright &#x000A9; 2023 Zhang, Guo, Zhuang, Liu, Xia, Zhang, Yang, Zeng, Wu, Huang and Li.</copyright-statement>
<copyright-year>2023</copyright-year>
<copyright-holder>Zhang, Guo, Zhuang, Liu, Xia, Zhang, Yang, Zeng, Wu, Huang and Li</copyright-holder>
<license xlink:href="http://creativecommons.org/licenses/by/4.0/"><p>This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.</p></license> </permissions>
<abstract>
<p>The gut microbiota affects the health and overall fitness of bumblebees. It can enhance the host&#x00027;s ecological range by leveraging their metabolic capacities. However, the diversity of the gut microbiota and adaptive functional evolution in high-altitude regions remain unclear. To explore how the gut microbiota helps the host adapt to high-altitude environments, we analyzed the differences in diversity and function of the gut microbiota between high- and low-altitude regions through full-length 16S rRNA sequencing. Our results show that high-altitude regions have a lower abundance of <italic>Fructobacillus</italic> and <italic>Saccharibacter</italic> compared to low-altitude regions. Additionally, some individuals in low-altitude regions were invaded by opportunistic pathogens. The gut microbiota in high-altitude regions has a greater number of pathways involved in &#x0201C;Protein digestion and absorption&#x0201D; and &#x0201C;Biosynthesis of amino acids,&#x0201D; while fewer carbohydrate pathways are involved in &#x0201C;digestion and absorption&#x0201D; and &#x0201C;Salmonella infection.&#x0201D; Our finding suggests that plateau hosts typically reduce energy metabolism and enhance immunity in response to adverse environments. Correspondingly, the gut microbiota also makes changes, such as reducing carbohydrate degradation and increasing protein utilization in response to the host. Additionally, the gut microbiota regulates their abundance and function to help the host adapt to adverse high-altitude environments.</p></abstract>
<kwd-group>
<kwd>gut microbiota</kwd>
<kwd>functional evolution</kwd>
<kwd>high-altitude</kwd>
<kwd>energy metabolism</kwd>
<kwd>immunity</kwd>
</kwd-group>
<counts>
<fig-count count="5"/>
<table-count count="0"/>
<equation-count count="0"/>
<ref-count count="51"/>
<page-count count="8"/>
<word-count count="4915"/>
</counts>
<custom-meta-wrap>
<custom-meta>
<meta-name>section-at-acceptance</meta-name>
<meta-value>Microbial Symbioses</meta-value>
</custom-meta>
</custom-meta-wrap>
</article-meta>
</front>
<body>
<sec id="s1">
<title>1. Introduction</title>
<p>Animals harbor gut microorganisms that have evolved along with host lineages, including humans, chimpanzees, gorillas, and orangutans. This provides evidence for the co-diversification of hosts and some lineages of the gut bacteria, implying a long-term vertical association (Moeller et al., <xref ref-type="bibr" rid="B30">2016</xref>; Moran et al., <xref ref-type="bibr" rid="B32">2019</xref>). Similarly, in social corbiculate bees (honey bees, bumble bees, and stingless bees), five core lineages of the gut microbiota show phylogenies mostly matching their hosts, supporting co-diversification over approximately 80 million years (Kwong et al., <xref ref-type="bibr" rid="B22">2017</xref>). The five core gut bacteria, <italic>Gilliamella, Snodgrassella, Bifidobacterium, Lactobacillus Firm-4</italic>, and <italic>Lactobacillus Firm-5</italic>, have undergone functional evolution to adapt to the host. By utilizing their own metabolic capabilities, the gut microbiota help hosts expand their ecological range, which benefits them (Kwong and Moran, <xref ref-type="bibr" rid="B24">2016</xref>). A previous study has indicated that symbiotic bacteria in insects are essential for host health by contributing to food digestion, detoxifying toxic molecules, providing essential nutrients, and protecting against pathogens and parasites (Engel and Moran, <xref ref-type="bibr" rid="B15">2013</xref>; Zheng et al., <xref ref-type="bibr" rid="B51">2016</xref>; Mockler et al., <xref ref-type="bibr" rid="B29">2018</xref>; Nishida and Ochman, <xref ref-type="bibr" rid="B34">2018</xref>).</p>
<p>Bumblebees are important pollinators of wild plants and crops. China has &#x0007E;50% of the world&#x00027;s species of bumblebees and is the country with the richest bumblebee species (Huang and An, <xref ref-type="bibr" rid="B19">2018</xref>). China is known for its abundant and distinct species of bumblebees, such as <italic>Bombus pyrosoma</italic>, which is found across a wide range, stretching from western Liaoning to eastern Qinghai provinces, and has also been reported in provinces such as Shanxi, Hebei, and Gansu (An et al., <xref ref-type="bibr" rid="B3">2008</xref>; Wu et al., <xref ref-type="bibr" rid="B44">2009</xref>). Their habitats span different ecosystems, from low elevations in the North China Plain to high elevations in the Qinghai&#x02013;Tibet Plateau (An et al., <xref ref-type="bibr" rid="B2">2014</xref>; Williams et al., <xref ref-type="bibr" rid="B43">2016</xref>). Bumblebees prefer cooler climates, leading to greater diversity in the temperate habitats of the Northern Hemisphere, particularly across Eurasia. In contrast, bumblebee species inhabiting tropical regions tend to be restricted to higher elevations due to the more favorable conditions found there (Engel and Rasmussen, <xref ref-type="bibr" rid="B14">2020</xref>).</p>
<p>High-altitude areas such as the Qinghai-Tibetan Plateau are often accompanied by extreme weather conditions of low temperature, hypoxia, and strong ultraviolet radiation, which are usually unsuitable for the growth and development of biological organisms (Yu et al., <xref ref-type="bibr" rid="B47">2023</xref>). However, organisms such as mammals and ruminants have evolved some strategies to adapt to high-altitude environments (Friedrich and Wiener, <xref ref-type="bibr" rid="B16">2020</xref>). One strategy for dealing with low oxygen levels is to increase oxygen availability in mammals (Storz et al., <xref ref-type="bibr" rid="B41">2010</xref>). For example, populations at high altitudes have lower ventilation than populations at low altitudes (Brutsaert, <xref ref-type="bibr" rid="B6">2007</xref>). In addition, oxygen affinity is enhanced by increasing red blood cell count and hemoglobin level (Weber, <xref ref-type="bibr" rid="B42">2007</xref>). Another strategy is to reduce the body&#x00027;s need for oxygen by lowering its overall metabolic rate (Hochachka et al., <xref ref-type="bibr" rid="B18">1996</xref>). Studies have shown that the energy expenditure of native Andean llamas is significantly lower compared to other ruminants (Riek et al., <xref ref-type="bibr" rid="B39">2019</xref>). High-altitude areas not only affect mammals but also their symbiotic microbiota. Compared to low-altitude individuals, high-altitude species such as humans, Tibetan sheep, and pigs have different compositions and diversity of their symbiotic microbiota (Yang et al., <xref ref-type="bibr" rid="B46">2017</xref>; Ma et al., <xref ref-type="bibr" rid="B27">2019</xref>). However, the metabolic role of the gut microbiota in animals living in high-altitude environments and the potential co-evolution between gut microbiota and host in high-altitude adaptation have not yet been fully understood.</p>
<p>The interaction between gut microbiota not only regulates the organism&#x00027;s health but also forms a crucial bridge between the environment and the host, helping the host to better adapt to the environment (Gao et al., <xref ref-type="bibr" rid="B17">2020</xref>). The unique habitat conditions and simple gut patterns of <italic>Bombus pyrosoma</italic> provide an excellent opportunity to study how gut microbiota help the host adapt to high-altitude evolutionary pressures (Zhang and Zheng, <xref ref-type="bibr" rid="B50">2022</xref>). We hypothesize that the composition and abundance of gut microbiota from <italic>Bombus pyrosoma</italic> in high-altitude areas may differ and help the host improve its ability to adapt to high-altitude environments.</p></sec>
<sec id="s2">
<title>2. Materials and methods</title>
<sec>
<title>2.1. <italic>B. pyrosoma</italic> sample collection and DNA extraction</title>
<p>A total of 41 workers of <italic>B. pyrosoma</italic> were collected from six different sites in September 2019 (<xref ref-type="fig" rid="F1">Figure 1</xref>). Three sites in Gansu belong to the high-altitude (HA, <italic>n</italic> = 21) Tibet Plateau (103&#x000B0;9&#x02032;29.79&#x02033;, 34&#x000B0;14&#x02032;23.92&#x02033;, 3,269 m, 102&#x000B0;39&#x02032;24.12&#x02033;, 36&#x000B0;58&#x02032;14.77&#x02033;, 2,942 m, 103&#x000B0;41&#x02032;21.48&#x02033;, 34&#x000B0;54&#x02032;34.06&#x02033;, 2,918 m), with each site containing four to ten worker bees as replicate samples. Another three sites in Hebei belong to the low-altitude (LA, <italic>n</italic> = 20) North China Plain (114&#x000B0;57&#x02032;5.40&#x02033;, 39&#x000B0;56&#x02032;42.04&#x02033;, 1,192 m, 117&#x000B0;35&#x02032;25.22&#x02033;, 42&#x000B0;7&#x02032;26.50&#x02033;, 1,027 m, 117&#x000B0;36&#x02032;53.70&#x02033;, 42&#x000B0;6&#x02032;54.95&#x02033;, 1,005 m) with five to eight replicate samples. The gut DNA of 41 workers of <italic>B. pyrosoma</italic> was extracted according to the protocol of the Wizard<sup>&#x000AE;</sup> Genomic DNA Purification Kit (Item number: A1125, Promega, Madison, Wisconsin, USA) (Paulos et al., <xref ref-type="bibr" rid="B36">2016</xref>). DNA extraction was dissolved in 30 ml Tris-EDTA (TE) buffer, quantified using a NanoDrop 2000 UV-visible (UV-vis) spectrophotometer (NanoDrop, DE, USA), and qualitatively evaluated by gel electrophoresis.</p>
<fig id="F1" position="float">
<label>Figure 1</label>
<caption><p>A topographic map of mainland China shows the sampling sites (red dots) of <italic>B. pyrosoma</italic> in this study. Three sites in Gansu represent <italic>B. pyrosoma</italic> individuals at high altitudes (HA, <italic>n</italic> = 21), with each site containing four to ten worker bees as replicate samples. Another three sites in Hebei represent <italic>B. pyrosoma</italic> individuals at low altitudes (LA, <italic>n</italic> = 20), with each site containing five to eight worker bees as replicate samples. The map was created with R (v 4.1.1).</p></caption>
<graphic mimetype="image" mime-subtype="tiff" xlink:href="fmicb-14-1218560-g0001.tif"/>
</fig></sec>
<sec>
<title>2.2. Full-length 16S rRna sequencing</title>
<p>PCR amplification of the nearly full-length bacterial 16S rRNA genes was conducted using the forward primer 27F (5&#x02032;-AGRGTTYGATYMTGGCTCAG-3&#x02032;) and the reverse primer 1492R (5&#x02032;-RGYTACCTTGTTACGACTT-3&#x02032;). The PCR reaction conditions were 95&#x000B0;C for 5 min, 95&#x000B0;C for 30 s, 50&#x000B0;C for 30 s, and 72&#x000B0;C for 1 min, with 25 cycles in a reaction volume of 10 &#x003BC;l. A total of PCR amplicons was purified with Agencourt AMPure Beads (Beckman Coulter, Indianapolis, IN) and quantified using the PicoGreen dsDNA Assay Kit (Invitrogen, Carlsbad, CA, USA). After the individual quantification step, amplicons were pooled in equal amounts, and Single Molecule Real Time (SMRT) sequencing technology was performed using the PacBio Sequel platform at Beijing Personal Biotechnology Co., Ltd. (Beijing, China). PacBio circular consensus sequencing (CCS) reads were derived from the multiple alignments of sub-reads to decrease the sequencing error rate.</p></sec>
<sec>
<title>2.3. Bioinformatics and statistical analysis</title>
<p>We used the complete DADA2 workflow (<ext-link ext-link-type="uri" xlink:href="https://github.com/benjjneb/LRASManuscript">https://github.com/benjjneb/LRASManuscript</ext-link>) to obtain amplicon sequence variants (ASVs) and feature abundance from PacBio CCS reads (Callahan et al., <xref ref-type="bibr" rid="B7">2019</xref>). The obtained ASVs were taxonomically annotated in the SILVA 138 Database (Quast et al., <xref ref-type="bibr" rid="B38">2013</xref>). The calculation of alpha diversity and beta diversity was based on the randomly extracted 4,530 sequencing abundances per sample. Analysis of similarities (ANOSIM) and permutational multivariate analysis of variance (PERMANOVA) were performed with R (v 4.1.1) scripts (Liu et al., <xref ref-type="bibr" rid="B26">2021</xref>). PICRUSt2 (Douglas et al., <xref ref-type="bibr" rid="B12">2020</xref>) contained the Kyoto Encyclopedia of Genes and Genomes (KEGG) orthologs (KO) database and was performed to predict categories and the abundance of Kyoto Encyclopedia of Genes and Genomes (KEGG) pathways (Kanehisa et al., <xref ref-type="bibr" rid="B20">2012</xref>). Specifically, we first used PICRUSt2 to predict KO with the sequence and abundance of ASVs as input and then converted KO to KEGG pathway categories by executing scripts. In addition, we used LDA effect size (LEfSe) (Chang et al., <xref ref-type="bibr" rid="B8">2022</xref>) to conduct different taxonomic hierarchical analyses and STAMP (Parks et al., <xref ref-type="bibr" rid="B35">2014</xref>) for functional differential analysis. The Kruskal&#x02013;Wallis test was conducted to test the difference in variables (e.g., alpha diversity and KEGG pathways) of gut microbiota found in high-altitude and of those found in low-altitude areas.</p></sec></sec>
<sec id="s3">
<title>3. Results</title>
<sec>
<title>3.1. Diversity differences in the gut microbiota between bumblebees in high-altitude and those in low-altitude areas</title>
<p>The rarefaction curves of bacterial ASVs were detected in the gut microbiota of high- and low-altitude bumblebees and reached a saturation phase by rarefying the samples to a minimum abundance of 4,530, which satisfied the subsequent analyses of gut microbiota alpha and beta diversity (<xref ref-type="fig" rid="F2">Figure 2A</xref>). The two groups have 42 shared ASVs, with 154 unique ASVs at high altitude regions and 236 unique ASVs at low altitude regions (<xref ref-type="fig" rid="F2">Figure 2B</xref>). The richness (Chao1) and phylogenetic diversity (faith_pd) of high-altitude bumblebee gut bacterial communities significantly decreased (<italic>P</italic> &#x0003C; 0.05) compared to those in low-altitude areas (<xref ref-type="fig" rid="F2">Figures 2C</xref>, <xref ref-type="fig" rid="F2">D</xref>). Analysis of similarities (ANOSIM) and permutational multivariate analysis of variance (PERMANOVA) were used to statistically test gut bacterial community similarities. The intestinal microbial community composition of bumblebees at high altitudes was significantly different from that in low-altitude areas (<italic>P</italic> &#x0003C; 0.05) (<xref ref-type="fig" rid="F3">Figure 3</xref>).</p>
<fig id="F2" position="float">
<label>Figure 2</label>
<caption><p>Alpha diversities of the gut microbiota differ between the HA (high-altitude, <italic>n</italic> = 21) and the LA (low-altitude, <italic>n</italic> = 20). <bold>(A)</bold> Rarefaction curves of detected bacterial ASVs of HA and LA; <bold>(B)</bold> A Venn diagram based on ASVS of HA and LA; <bold>(C)</bold> Chao1 index of the gut microbiota between HA and LA (<italic>P</italic> &#x0003C; 0.05); <bold>(D)</bold> Faith_pd index of the gut microbiota between HA and LA (<italic>P</italic> &#x0003C; 0.05). a and b are methods of representation in statistics, and they represent significant differences.</p></caption>
<graphic mimetype="image" mime-subtype="tiff" xlink:href="fmicb-14-1218560-g0002.tif"/>
</fig>
<fig id="F3" position="float">
<label>Figure 3</label>
<caption><p>Beta diversities of the gut microbiota differ between the HA (high-altitude, <italic>n</italic> = 21) and the LA (low-altitude, <italic>n</italic> = 20) based on unweighted UniFrac distance. The horizontal axis represents the number of pairwise distances between intestinal individuals at HA and LA, respectively.</p></caption>
<graphic mimetype="image" mime-subtype="tiff" xlink:href="fmicb-14-1218560-g0003.tif"/>
</fig></sec>
<sec>
<title>3.2. Taxonomy difference of the gut microbiota of bumblebees between high- and low-altitude areas</title>
<p>The top six average relative abundances of the gut microbiota of <italic>B. pyrosoma</italic> are the core symbiotic bacteria <italic>Snodgrassella, Lactobacillus, Gilliamella, Bombiscardovi, Apibacter</italic>, and <italic>Serratia</italic> (<xref ref-type="fig" rid="F4">Figure 4A</xref>). LEfSe (LDA score &#x0003E; 2) analysis of the different classification hierarchies showed that the gut community of high-altitude bumblebees has fewer Leuconostocaceae (Family), <italic>Fructobacillus</italic> (Genus), Alphaproteobacteria (Class), Acetobacterales (Order), Acetobacteraceae (Family), and <italic>Saccharibacter</italic> (Genus) compared with that of low-altitude bumblebees (<xref ref-type="fig" rid="F4">Figure 4B</xref>). Furthermore, the boxplot showed that the bumblebee gut at high-altitude regions had less <italic>Fructobacillus</italic> (<italic>P</italic> &#x0003C; 0.05) and <italic>Saccharibacter</italic> (<italic>P</italic> &#x0003C; 0.05) compared to that at low-altitude regions (<xref ref-type="fig" rid="F4">Figures 4C</xref>, <xref ref-type="fig" rid="F4">D</xref>).</p>
<fig id="F4" position="float">
<label>Figure 4</label>
<caption><p>Variation of microbial taxonomy between the HA (high altitude, <italic>n</italic> = 21) and the LA (low altitude, <italic>n</italic> = 20). <bold>(A)</bold> Genus-level relative abundance of gut microbiota of the HA and LA; <bold>(B)</bold> LEfSe (LDA score &#x0003E; 2) analysis of the different classification hierarchy; <bold>(C)</bold> <italic>Fructobacillus</italic> relative abundance of gut microbiota at the HA and LA (<italic>P</italic> &#x0003C; 0.05); <bold>(D)</bold> <italic>Saccharibacter</italic> relative abundance of gut microbiota at the HA and LA (<italic>P</italic> &#x0003C; 0.05). a and b are methods of representation in statistics, and they represent significant differences.</p></caption>
<graphic mimetype="image" mime-subtype="tiff" xlink:href="fmicb-14-1218560-g0004.tif"/>
</fig></sec>
<sec>
<title>3.3. KEGG pathway difference in the gut microbiota of bumblebees between high- and low-altitude areas</title>
<p>The intestinal microbes of bumblebees at high altitudes had three significant KEGG secondary classification pathways compared with those at low altitudes (<italic>P</italic> &#x0003C; 0.05). The gut microbiota of bumblebees at high altitudes has a greater number of pathways involving &#x0201C;Drug resistance: antineoplastic&#x0201D; and &#x0201C;Cell growth and death,&#x0201D; while the gut microbiota of bumblebees in low-altitude areas has a greater number of &#x0201C;Infectious disease: bacterial&#x0201D; (<xref ref-type="fig" rid="F5">Figure 5A</xref>). For third-level KEGG pathways, the gut microbiota of bumblebee individuals in high-altitude regions exhibits a higher abundance of pathways related to &#x0201C;One carbon pool by folate,&#x0201D; &#x0201C;Antifolate resistance,&#x0201D; &#x0201C;Protein digestion and absorption,&#x0201D; and &#x0201C;Biosynthesis of amino acids.&#x0201D; Conversely, these high-altitude regions show a lower abundance of pathways associated with &#x0201C;Carbohydrate digestion and absorption&#x0201D; and &#x0201C;Salmonella infection&#x0201D; compared to the gut microbiota of bumblebees in low-altitude regions (<xref ref-type="fig" rid="F5">Figure 5B</xref>).</p>
<fig id="F5" position="float">
<label>Figure 5</label>
<caption><p>Functionally predicted KEGG pathways differing in gut microbiota between HA (high-altitude, <italic>n</italic> = 21) and LA (low-altitude, <italic>n</italic> = 20). The difference in proportions between the two groups is shown with 95% confidence intervals. Only a <italic>P</italic>-value of &#x0003C;0.05 (FDR adjusted) is shown and composition. <bold>(A)</bold> Secondary-level KEGG pathways; <bold>(B)</bold> Third-level KEGG pathways.</p></caption>
<graphic mimetype="image" mime-subtype="tiff" xlink:href="fmicb-14-1218560-g0005.tif"/>
</fig></sec></sec>
<sec id="s4">
<title>4. Discussion</title>
<p>Compared with bumblebees in low-altitude regions, the gut bacteria of <italic>B. pyrosoma</italic> in high-altitude areas have lower alpha diversity and a relatively uniform community composition, which may result from collecting only a single type of plant pollen and nectar. Apart from host attributes (e.g., species, age, and caste) determining the gut microbiota of bumblebees, selective forces (e.g., antibiotic-resistant genes and habitat ecological environment) may also play an important role in shaping the formation of the gut microbial community (Colman et al., <xref ref-type="bibr" rid="B10">2012</xref>; Kwong et al., <xref ref-type="bibr" rid="B21">2014</xref>; Kwong and Moran, <xref ref-type="bibr" rid="B23">2015</xref>, <xref ref-type="bibr" rid="B24">2016</xref>; Zhang et al., <xref ref-type="bibr" rid="B49">2021</xref>). For example, differences in external environments lead to different microbial community structures in wild and artificially reared bee colonies, with wild colonies having more non-core microbial groups (Zhang and Zheng, <xref ref-type="bibr" rid="B50">2022</xref>). Owing to the obvious difference in eco-climate between high- and low-altitude areas, low temperature, hypoxia, ultraviolet radiation, and other climatic characteristics in low-altitude areas lead to the decline of plant diversity (Bruns and Kennedy, <xref ref-type="bibr" rid="B5">2009</xref>; Murray, <xref ref-type="bibr" rid="B33">2016</xref>). However, plant diversity is higher in low-altitude areas, where the gut microbiota displays higher diversity and has a more dispersed community composition.</p>
<p>The high-altitude regions have fewer <italic>Fructobacillus</italic> and <italic>Saccharibacter</italic>. It may be because there are fewer carbon sources at higher altitudes. The gut microbiota, such as <italic>Fructobacillus</italic> and <italic>Saccharibacter</italic>, generally benefit the health of their hosts. Members of the genus <italic>Frugibacillus</italic> are classified as frugiophilic lactic acid bacteria (Endo et al., <xref ref-type="bibr" rid="B13">2018</xref>). The fructose content found in the nectar of quince cultivars from undamaged varieties was higher than that from cultivars found under frozen conditions (Ak&#x00161;i&#x00107; et al., <xref ref-type="bibr" rid="B1">2014</xref>). To adapt to fructose-rich environments such as fruit and honey-rich flowers, <italic>Fructobacillus</italic> underwent a specific reductive evolution similar to that of <italic>Lactobacillus kunkeei</italic> (Lee et al., <xref ref-type="bibr" rid="B25">2015</xref>; Maeno et al., <xref ref-type="bibr" rid="B28">2016</xref>). <italic>Frutobacillus</italic> and <italic>Saccharibacter</italic>, present in the bumblebee gut, help the host degrade fructose through carbohydrate catabolism (Praet et al., <xref ref-type="bibr" rid="B37">2018</xref>). The genus <italic>Saccharibacter</italic> is an acetic acid bacterium that has been shown to be associated with a variety of insects that rely on a high-sugar diet (Crotti et al., <xref ref-type="bibr" rid="B11">2010</xref>; Chouaia et al., <xref ref-type="bibr" rid="B9">2014</xref>). Since the genus <italic>Saccharibacter</italic> is isolated from flowers, it may also come into contact with insect pollinators and play a role in their use of nectar as a food source (Smith et al., <xref ref-type="bibr" rid="B40">2020</xref>). However, some individuals in low-altitude areas even have a predominance of opportunistic bacteria such as <italic>Serratia</italic>, which may indicate that bumblebees in low-altitude areas face more severe environments such as ecological pollution.</p>
<p>Plateau-adapted hosts have developed strategies to cope with adverse environments, such as lowering their overall metabolic rate to reduce the body&#x00027;s demand for oxygen (Moore, <xref ref-type="bibr" rid="B31">2017</xref>). Under low temperatures and hypoxic conditions, insects upregulate the expression of immune-related genes, such as those responsible for encoding antimicrobial peptides and the activation of the IKK/NF-kB signaling pathway, which in turn stimulates the immune response (Zhang et al., <xref ref-type="bibr" rid="B48">2011</xref>; Xu and James, <xref ref-type="bibr" rid="B45">2012</xref>; Bandarra et al., <xref ref-type="bibr" rid="B4">2014</xref>). Compared to gut bacteria of bumblebee individuals in low-altitude areas, those in high-altitude regions have a higher abundance of pathways associated with &#x0201C;One carbon pool by folate,&#x0201D; &#x0201C;Antifolate resistance,&#x0201D; &#x0201C;Protein digestion and absorption,&#x0201D; and &#x0201C;Biosynthesis of amino acids.&#x0201D; Conversely, there is a decrease in pathways related to &#x0201C;Carbohydrate digestion and absorption&#x0201D; and &#x0201C;Salmonella infection.&#x0201D; This variation may be due to adverse environmental conditions such as low temperatures, hypoxia, and pollen and nectar food shortages at high altitudes. The Tibetan <italic>B. pyrosoma</italic> reduces energy metabolism and enhances immunity, which necessitates fewer carbohydrates and more protein to cope with these challenges. Correspondingly, the gut microbiota regulates members&#x00027; species and abundance (e.g., <italic>Fructobacillus</italic> and <italic>Saccharibacter</italic>), exerting less demand for carbohydrates and improving protein utilization to help the host cope with an adverse high-altitude environment.</p></sec>
<sec id="s5">
<title>5. Conclusion</title>
<p>The gut microbiota co-evolved with their hosts for over 80 million years. Intestinal bacteria expand the host&#x00027;s ecological range through their metabolic capacities. Similarly, the host&#x00027;s ecological environment can also affect the composition and abundance of intestinal bacteria. Hosts at high altitudes, such as <italic>Bombus pyrosoma</italic>, adapt to the plateau environment by reducing energy metabolism and improving immunity. In turn, the gut microbiota helps the host adapt to the adverse environment by regulating their members&#x00027; abundance and functional metabolism, such as reducing carbohydrate degradation and improving the utilization of protein. Our findings provide insights into the interactions between the bumblebee gut microbiota and the host.</p></sec>
<sec sec-type="data-availability" id="s6">
<title>Data availability statement</title>
<p>The PacBio circular consensus sequencing (CCS) reads have been deposited in the Genome Sequence Archive in the BIG Data Center, Beijing Institute of Genomics (BIG), Chinese Academy of Sciences, under accession numbers CRA010937 that are publicly accessible at <ext-link ext-link-type="uri" xlink:href="http://bigd.big.ac.cn/gsa">http://bigd.big.ac.cn/gsa</ext-link>, accessed on 7 May 2023.</p></sec>
<sec sec-type="author-contributions" id="s7">
<title>Author contributions</title>
<p>ZheZ and FY conceived the research project. ZheZ, MZ, FL, ZX, ZhiZ, HZ, YW, and JH collected the samples, analyzed the data, and sequenced the samples. ZheZ, YG, and JL wrote the manuscript. All authors read and approved the final version of the manuscript.</p></sec>
</body>
<back>
<sec sec-type="funding-information" id="s8">
<title>Funding</title>
<p>This study was supported by the Agricultural Science and Technology Innovation Program (CAAS&#x02013;ASTIP-2016&#x02013;IAR) and the Agriculture Research System Bee (CARS-44-KXJ5).</p>
</sec>
<ack><p>We thank Yuanjian Li and Yingping Qu for collecting samples and providing help and advice in the laboratory.</p>
</ack>
<sec sec-type="COI-statement" id="conf1">
<title>Conflict of interest</title>
<p>MZ was employed by Shanghai Suosheng Biotechnology Co., Ltd. HZ and YW were employed by Luoping Yunling Bee Industry and Trade Co., Ltd. The remaining authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.</p>
</sec>
<sec sec-type="disclaimer" id="s9">
<title>Publisher&#x00027;s note</title>
<p>All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.</p>
</sec>
<ref-list>
<title>References</title>
<ref id="B1">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Ak&#x00161;i&#x00107;</surname> <given-names>M. F.</given-names></name> <name><surname>Tosti</surname> <given-names>T.</given-names></name> <name><surname>Nedi&#x00107;</surname> <given-names>N.</given-names></name> <name><surname>Markovi&#x00107;</surname> <given-names>M.</given-names></name> <name><surname>Li&#x0010D;ina</surname> <given-names>V.</given-names></name> <name><surname>Milojkovi&#x00107;-Opsenica</surname> <given-names>D.</given-names></name> <etal/></person-group>. (<year>2014</year>). <article-title>Influence of frost damage on the sugars and sugar alcohol composition in quince (Cydonia oblonga Mill.) floral nectar</article-title>. <source>Acta. Physiologiae Plantarum</source> 37. <pub-id pub-id-type="doi">10.1007/s11738-014-1701-y</pub-id></citation>
</ref>
<ref id="B2">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>An</surname> <given-names>J.</given-names></name> <name><surname>Huang</surname> <given-names>J.</given-names></name> <name><surname>Shao</surname> <given-names>Y.</given-names></name> <name><surname>Zhang</surname> <given-names>S.</given-names></name> <name><surname>Wang</surname> <given-names>B.</given-names></name> <name><surname>Liu</surname> <given-names>X.</given-names></name> <etal/></person-group>. (<year>2014</year>). <article-title>The bumblebees of North China (Apidae, Bombus Latreille)</article-title>. <source>Zootaxa</source> <volume>3830</volume>, <fpage>1</fpage>&#x02013;<lpage>89</lpage>. <pub-id pub-id-type="doi">10.11646/zootaxa.3830.1.1</pub-id><pub-id pub-id-type="pmid">25081273</pub-id></citation></ref>
<ref id="B3">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>An</surname> <given-names>J.</given-names></name> <name><surname>Yao</surname> <given-names>J.</given-names></name> <name><surname>Huang</surname> <given-names>J.</given-names></name> <name><surname>Shao</surname> <given-names>Y.</given-names></name> <name><surname>Wu</surname> <given-names>J.</given-names></name> <name><surname>Li</surname> <given-names>J.</given-names></name> <etal/></person-group>. (<year>2008</year>). <article-title><italic>Bombus fauna</italic> (Hymenoptera, Apidae) in Shanxi, China</article-title>. <source>Acta. Zootaxonomica Sinica</source>, <volume>33</volume>, <fpage>80</fpage>&#x02013;<lpage>88</lpage>.</citation>
</ref>
<ref id="B4">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Bandarra</surname> <given-names>D.</given-names></name> <name><surname>Biddlestone</surname> <given-names>J.</given-names></name> <name><surname>Mudie</surname> <given-names>S.</given-names></name> <name><surname>Muller</surname> <given-names>H. A.</given-names></name> <name><surname>Rocha</surname> <given-names>S.</given-names></name></person-group> (<year>2014</year>). <article-title>Hypoxia activates IKK-NF-kappaB and the immune response in Drosophila melanogaster</article-title>. <source>Biosci. Rep.</source> <volume>34</volume>, <fpage>e00127</fpage>. <pub-id pub-id-type="doi">10.1042/BSR20140095</pub-id><pub-id pub-id-type="pmid">24993778</pub-id></citation></ref>
<ref id="B5">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Bruns</surname> <given-names>T. D.</given-names></name> <name><surname>Kennedy</surname> <given-names>P. G.</given-names></name></person-group> (<year>2009</year>). <article-title>Individuals, populations, communities and function: the growing field of ectomycorrhizal ecology</article-title>. <source>New Phytol.</source> <volume>182</volume>, <fpage>12</fpage>&#x02013;<lpage>14</lpage>. <pub-id pub-id-type="doi">10.1111/j.1469-8137.2009.02788.x</pub-id><pub-id pub-id-type="pmid">19291070</pub-id></citation></ref>
<ref id="B6">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Brutsaert</surname> <given-names>T. D.</given-names></name></person-group> (<year>2007</year>). <article-title>Population genetic aspects and phenotypic plasticity of ventilatory responses in high altitude natives</article-title>. <source>Respir. Physiol. Neurobiol.</source> <volume>158</volume>, <fpage>151</fpage>&#x02013;<lpage>160</lpage>. <pub-id pub-id-type="doi">10.1016/j.resp.2007.03.004</pub-id><pub-id pub-id-type="pmid">17400521</pub-id></citation></ref>
<ref id="B7">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Callahan</surname> <given-names>B. J.</given-names></name> <name><surname>Wong</surname> <given-names>J.</given-names></name> <name><surname>Heiner</surname> <given-names>C.</given-names></name> <name><surname>Oh</surname> <given-names>S.</given-names></name> <name><surname>Theriot</surname> <given-names>C. M.</given-names></name> <name><surname>Gulati</surname> <given-names>A. S.</given-names></name> <etal/></person-group>. (<year>2019</year>). <article-title>High-throughput amplicon sequencing of the full-length 16S rRNA gene with single-nucleotide resolution</article-title>. <source>Nucleic Acids Res.</source> <volume>47</volume>, <fpage>e103</fpage>. <pub-id pub-id-type="doi">10.1093/nar/gkz569</pub-id><pub-id pub-id-type="pmid">31269198</pub-id></citation></ref>
<ref id="B8">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Chang</surname> <given-names>F.</given-names></name> <name><surname>He</surname> <given-names>S.</given-names></name> <name><surname>Dang</surname> <given-names>C.</given-names></name></person-group> (<year>2022</year>). <article-title>Assisted selection of biomarkers by linear discriminant analysis effect size (LEfSe) in microbiome data</article-title>. <source>J. Vis. Exp</source>. e61715. <pub-id pub-id-type="doi">10.3791/61715</pub-id><pub-id pub-id-type="pmid">35635468</pub-id></citation></ref>
<ref id="B9">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Chouaia</surname> <given-names>B.</given-names></name> <name><surname>Gaiarsa</surname> <given-names>S.</given-names></name> <name><surname>Crotti</surname> <given-names>E.</given-names></name> <name><surname>Comandatore</surname> <given-names>F.</given-names></name> <name><surname>Degli Esposti</surname> <given-names>M.</given-names></name> <name><surname>Ricci</surname> <given-names>I.</given-names></name> <etal/></person-group>. (<year>2014</year>). <article-title>Acetic acid bacteria genomes reveal functional traits for adaptation to life in insect guts</article-title>. <source>Genome Biol. Evol.</source> <volume>6</volume>, <fpage>912</fpage>&#x02013;<lpage>920</lpage>. <pub-id pub-id-type="doi">10.1093/gbe/evu062</pub-id><pub-id pub-id-type="pmid">24682158</pub-id></citation></ref>
<ref id="B10">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Colman</surname> <given-names>D. R.</given-names></name> <name><surname>Toolson</surname> <given-names>E. C.</given-names></name> <name><surname>Takacs-Vesbach</surname> <given-names>C. D.</given-names></name></person-group> (<year>2012</year>). <article-title>Do diet and taxonomy influence insect gut bacterial communities?</article-title> <source>Mol. Ecol.</source> <volume>21</volume>, <fpage>5124</fpage>&#x02013;<lpage>5137</lpage>. <pub-id pub-id-type="doi">10.1111/j.1365-294X.2012.05752.x</pub-id><pub-id pub-id-type="pmid">22978555</pub-id></citation></ref>
<ref id="B11">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Crotti</surname> <given-names>E.</given-names></name> <name><surname>Rizzi</surname> <given-names>A.</given-names></name> <name><surname>Chouaia</surname> <given-names>B.</given-names></name> <name><surname>Ricci</surname> <given-names>I.</given-names></name> <name><surname>Favia</surname> <given-names>G.</given-names></name> <name><surname>Alma</surname> <given-names>A.</given-names></name> <etal/></person-group>. (<year>2010</year>). <article-title>Acetic acid bacteria, newly emerging symbionts of insects</article-title>. <source>Appl. Environ. Microbiol.</source> <volume>76</volume>, <fpage>6963</fpage>&#x02013;<lpage>6970</lpage>. <pub-id pub-id-type="doi">10.1128/AEM.01336-10</pub-id><pub-id pub-id-type="pmid">20851977</pub-id></citation></ref>
<ref id="B12">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Douglas</surname> <given-names>G. M.</given-names></name> <name><surname>Maffei</surname> <given-names>V. J.</given-names></name> <name><surname>Zaneveld</surname> <given-names>J. R.</given-names></name> <name><surname>Yurgel</surname> <given-names>S. N.</given-names></name> <name><surname>Brown</surname> <given-names>J. R.</given-names></name> <name><surname>Taylor</surname> <given-names>C. M.</given-names></name> <etal/></person-group>. (<year>2020</year>). <article-title>PICRUSt2 for prediction of metagenome functions</article-title>. <source>Nat. Biotechnol.</source> <volume>38</volume>, <fpage>685</fpage>&#x02013;<lpage>688</lpage>. <pub-id pub-id-type="doi">10.1038/s41587-020-0548-6</pub-id><pub-id pub-id-type="pmid">32483366</pub-id></citation></ref>
<ref id="B13">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Endo</surname> <given-names>A.</given-names></name> <name><surname>Maeno</surname> <given-names>S.</given-names></name> <name><surname>Tanizawa</surname> <given-names>Y.</given-names></name> <name><surname>Kneifel</surname> <given-names>W.</given-names></name> <name><surname>Arita</surname> <given-names>M.</given-names></name> <name><surname>Dicks</surname> <given-names>L.</given-names></name> <etal/></person-group>. (<year>2018</year>). <article-title>Fructophilic lactic acid bacteria, a unique group of fructose-fermenting microbes</article-title>. <source>Appl. Environ. Microbiol.</source> <volume>84</volume>:<fpage>e01290</fpage>-<lpage>18</lpage>. <pub-id pub-id-type="doi">10.1128/AEM.01290-18</pub-id><pub-id pub-id-type="pmid">30054367</pub-id></citation></ref>
<ref id="B14">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Engel</surname> <given-names>M. S.</given-names></name> <name><surname>Rasmussen</surname> <given-names>C.</given-names></name></person-group> (<year>2020</year>). <article-title>Corbiculate bees</article-title>. <source>Encyclopedia of Social Insects</source> <pub-id pub-id-type="doi">10.1007/978-3-319-90306-4_30-1</pub-id></citation>
</ref>
<ref id="B15">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Engel</surname> <given-names>P.</given-names></name> <name><surname>Moran</surname> <given-names>N. A.</given-names></name></person-group> (<year>2013</year>). <article-title>Functional and evolutionary insights into the simple yet specific gut microbiota of the honey bee from metagenomic analysis</article-title>. <source>Gut Microbes</source> <volume>4</volume>, <fpage>60</fpage>&#x02013;<lpage>65</lpage>. <pub-id pub-id-type="doi">10.4161/gmic.22517</pub-id><pub-id pub-id-type="pmid">23060052</pub-id></citation></ref>
<ref id="B16">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Friedrich</surname> <given-names>J.</given-names></name> <name><surname>Wiener</surname> <given-names>P.</given-names></name></person-group> (<year>2020</year>). <article-title>Selection signatures for high-altitude adaptation in ruminants</article-title>. <source>Anim. Genet.</source> <volume>51</volume>, <fpage>157</fpage>&#x02013;<lpage>165</lpage>. <pub-id pub-id-type="doi">10.1111/age.12900</pub-id><pub-id pub-id-type="pmid">31943284</pub-id></citation></ref>
<ref id="B17">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Gao</surname> <given-names>H.</given-names></name> <name><surname>Chi</surname> <given-names>X.</given-names></name> <name><surname>Li</surname> <given-names>G.</given-names></name> <name><surname>Qin</surname> <given-names>W.</given-names></name> <name><surname>Song</surname> <given-names>P.</given-names></name> <name><surname>Jiang</surname> <given-names>F.</given-names></name> <etal/></person-group>. (<year>2020</year>). <article-title>Gut microbial diversity and stabilizing functions enhance the plateau adaptability of Tibetan wild ass (Equus kiang)</article-title>. <source>Microbiologyopen</source> <volume>9</volume>, <fpage>1150</fpage>&#x02013;<lpage>1161</lpage>. <pub-id pub-id-type="doi">10.1002/mbo3.1025</pub-id><pub-id pub-id-type="pmid">32157819</pub-id></citation></ref>
<ref id="B18">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Hochachka</surname> <given-names>P. W.</given-names></name> <name><surname>Buck</surname> <given-names>L. T.</given-names></name> <name><surname>Doll</surname> <given-names>C. J.</given-names></name> <name><surname>Land</surname> <given-names>S. C.</given-names></name></person-group> (<year>1996</year>). <article-title>Unifying theory of hypoxia tolerance: molecular/metabolic defense and rescue mechanisms for surviving oxygen lack</article-title>. <source>Proc. Natl. Acad. Sci. USA.</source> <volume>93</volume>, <fpage>9493</fpage>&#x02013;<lpage>9498</lpage>. <pub-id pub-id-type="doi">10.1073/pnas.93.18.9493</pub-id><pub-id pub-id-type="pmid">8790358</pub-id></citation></ref>
<ref id="B19">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Huang</surname> <given-names>J. X.</given-names></name> <name><surname>An</surname> <given-names>J.</given-names></name></person-group> (<year>2018</year>). <article-title>Species diversity, pollination application and strategy for conservation of the bumblebees of China</article-title>. <source>Biodivers. Sci.</source> <volume>26</volume>, <fpage>486</fpage>&#x02013;<lpage>497</lpage>. <pub-id pub-id-type="doi">10.17520/biods.2018068</pub-id></citation>
</ref>
<ref id="B20">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Kanehisa</surname> <given-names>M.</given-names></name> <name><surname>Goto</surname> <given-names>S.</given-names></name> <name><surname>Sato</surname> <given-names>Y.</given-names></name> <name><surname>Furumichi</surname> <given-names>M.</given-names></name> <name><surname>Tanabe</surname> <given-names>M.</given-names></name></person-group> (<year>2012</year>). <article-title>KEGG for integration and interpretation of large-scale molecular data sets</article-title>. <source>Nucleic Acids Res.</source> <volume>40</volume>, <fpage>D109</fpage>&#x02013;<lpage>D114</lpage>. <pub-id pub-id-type="doi">10.1093/nar/gkr988</pub-id><pub-id pub-id-type="pmid">22080510</pub-id></citation></ref>
<ref id="B21">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Kwong</surname> <given-names>W. K.</given-names></name> <name><surname>Engel</surname> <given-names>P.</given-names></name> <name><surname>Koch</surname> <given-names>H.</given-names></name> <name><surname>Moran</surname> <given-names>N. A.</given-names></name></person-group> (<year>2014</year>). <article-title>Genomics and host specialization of honey bee and bumble bee gut symbionts</article-title>. <source>Proc. Natl. Acad. Sci. USA.</source> <volume>111</volume>, <fpage>11509</fpage>&#x02013;<lpage>11514</lpage>. <pub-id pub-id-type="doi">10.1073/pnas.1405838111</pub-id><pub-id pub-id-type="pmid">25053814</pub-id></citation></ref>
<ref id="B22">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Kwong</surname> <given-names>W. K.</given-names></name> <name><surname>Medina</surname> <given-names>L. A.</given-names></name> <name><surname>Koch</surname> <given-names>H.</given-names></name> <name><surname>Sing</surname> <given-names>K. W.</given-names></name> <name><surname>Soh</surname> <given-names>E. J. Y.</given-names></name> <name><surname>Ascher</surname> <given-names>J. S.</given-names></name> <etal/></person-group>. (<year>2017</year>). <article-title>Dynamic microbiome evolution in social bees</article-title>. <source>Sci. Adv.</source> 3,e1600513. <pub-id pub-id-type="doi">10.1126/sciadv.1600513</pub-id><pub-id pub-id-type="pmid">28435856</pub-id></citation></ref>
<ref id="B23">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Kwong</surname> <given-names>W. K.</given-names></name> <name><surname>Moran</surname> <given-names>N. A.</given-names></name></person-group> (<year>2015</year>). <article-title>Evolution of host specialization in gut microbes: the bee gut as a model</article-title>. <source>Gut Microbes</source> <volume>6</volume>, <fpage>214</fpage>&#x02013;<lpage>220</lpage>. <pub-id pub-id-type="doi">10.1080/19490976.2015.1047129</pub-id><pub-id pub-id-type="pmid">26011669</pub-id></citation></ref>
<ref id="B24">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Kwong</surname> <given-names>W. K.</given-names></name> <name><surname>Moran</surname> <given-names>N. A.</given-names></name></person-group> (<year>2016</year>). <article-title>Gut microbial communities of social bees</article-title>. <source>Nat. Rev. Microbiol.</source> <volume>14</volume>, <fpage>374</fpage>&#x02013;<lpage>384</lpage>. <pub-id pub-id-type="doi">10.1038/nrmicro.2016.43</pub-id><pub-id pub-id-type="pmid">27140688</pub-id></citation></ref>
<ref id="B25">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Lee</surname> <given-names>F. J.</given-names></name> <name><surname>Rusch</surname> <given-names>D. B.</given-names></name> <name><surname>Stewart</surname> <given-names>F. J.</given-names></name> <name><surname>Mattila</surname> <given-names>H. R.</given-names></name> <name><surname>Newton</surname> <given-names>I. L.</given-names></name></person-group> (<year>2015</year>). <article-title>Saccharide breakdown and fermentation by the honey bee gut microbiome</article-title>. <source>Environ. Microbiol.</source> <volume>17</volume>, <fpage>796</fpage>&#x02013;<lpage>815</lpage>. <pub-id pub-id-type="doi">10.1111/1462-2920.12526</pub-id><pub-id pub-id-type="pmid">24905222</pub-id></citation></ref>
<ref id="B26">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Liu</surname> <given-names>Y. X.</given-names></name> <name><surname>Qin</surname> <given-names>Y.</given-names></name> <name><surname>Chen</surname> <given-names>T.</given-names></name> <name><surname>Lu</surname> <given-names>M.</given-names></name> <name><surname>Qian</surname> <given-names>X.</given-names></name> <name><surname>Guo</surname> <given-names>X.</given-names></name> <etal/></person-group>. (<year>2021</year>). <article-title>A practical guide to amplicon and metagenomic analysis of microbiome data</article-title>. <source>Protein Cell</source> <volume>12</volume>, <fpage>315</fpage>&#x02013;<lpage>330</lpage>. <pub-id pub-id-type="doi">10.1007/s13238-020-00724-8</pub-id><pub-id pub-id-type="pmid">32394199</pub-id></citation></ref>
<ref id="B27">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Ma</surname> <given-names>Y.</given-names></name> <name><surname>Ma</surname> <given-names>S.</given-names></name> <name><surname>Chang</surname> <given-names>L.</given-names></name> <name><surname>Wang</surname> <given-names>H.</given-names></name> <name><surname>Ga</surname> <given-names>Q.</given-names></name> <name><surname>Ma</surname> <given-names>L.</given-names></name> <etal/></person-group>. (<year>2019</year>). <article-title>Gut microbiota adaptation to high altitude in indigenous animals</article-title>. <source>Biochem. Biophys. Res. Commun.</source> <volume>516</volume>, <fpage>120</fpage>&#x02013;<lpage>126</lpage>. <pub-id pub-id-type="doi">10.1016/j.bbrc.2019.05.085</pub-id><pub-id pub-id-type="pmid">31196622</pub-id></citation></ref>
<ref id="B28">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Maeno</surname> <given-names>S.</given-names></name> <name><surname>Tanizawa</surname> <given-names>Y.</given-names></name> <name><surname>Kanesaki</surname> <given-names>Y.</given-names></name> <name><surname>Kubota</surname> <given-names>E.</given-names></name> <name><surname>Kumar</surname> <given-names>H.</given-names></name> <name><surname>Dicks</surname> <given-names>L.</given-names></name> <etal/></person-group>. (<year>2016</year>). <article-title>Genomic characterization of a fructophilic bee symbiont Lactobacillus kunkeei reveals its niche-specific adaptation</article-title>. <source>Syst. Appl. Microbiol.</source> <volume>39</volume>, <fpage>516</fpage>&#x02013;<lpage>526</lpage>. <pub-id pub-id-type="doi">10.1016/j.syapm.2016.09.006</pub-id><pub-id pub-id-type="pmid">27776911</pub-id></citation></ref>
<ref id="B29">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Mockler</surname> <given-names>B. K.</given-names></name> <name><surname>Kwong</surname> <given-names>W. K.</given-names></name> <name><surname>Moran</surname> <given-names>N. A.</given-names></name> <name><surname>Koch</surname> <given-names>H.</given-names></name></person-group> (<year>2018</year>). <article-title>Microbiome structure influences infection by the parasite crithidia bombi in bumble bees</article-title>. <source>Appl. Environ. Microbiol.</source> <volume>84</volume>, <fpage>e02335</fpage>-<lpage>17</lpage>. <pub-id pub-id-type="doi">10.1128/AEM.02335-17</pub-id><pub-id pub-id-type="pmid">29374030</pub-id></citation></ref>
<ref id="B30">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Moeller</surname> <given-names>A. H.</given-names></name> <name><surname>Foerster</surname> <given-names>S.</given-names></name> <name><surname>Wilson</surname> <given-names>M. L.</given-names></name> <name><surname>Pusey</surname> <given-names>A. E.</given-names></name> <name><surname>Hahn</surname> <given-names>B. H.</given-names></name> <name><surname>Ochman</surname> <given-names>H.</given-names></name></person-group> (<year>2016</year>). <article-title>Social behavior shapes the chimpanzee pan-microbiome</article-title>. <source>Sci. Adv.</source> <volume>2</volume>, <fpage>e1500997</fpage>. <pub-id pub-id-type="doi">10.1126/sciadv.1500997</pub-id><pub-id pub-id-type="pmid">26824072</pub-id></citation></ref>
<ref id="B31">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Moore</surname> <given-names>L. G.</given-names></name></person-group> (<year>2017</year>). <article-title>Measuring high-altitude adaptation</article-title>. <source>J. Appl. Physiol. (1985)</source> <volume>123</volume>, <fpage>1371</fpage>&#x02013;<lpage>1385</lpage>. <pub-id pub-id-type="doi">10.1152/japplphysiol.00321.2017</pub-id><pub-id pub-id-type="pmid">28860167</pub-id></citation></ref>
<ref id="B32">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Moran</surname> <given-names>N. A.</given-names></name> <name><surname>Ochman</surname> <given-names>H.</given-names></name> <name><surname>Hammer</surname> <given-names>T. J.</given-names></name></person-group> (<year>2019</year>). <article-title>Evolutionary and ecological consequences of gut microbial communities</article-title>. <source>Annu. Rev. Ecol. Evol. Syst.</source> <volume>50</volume>, <fpage>451</fpage>&#x02013;<lpage>475</lpage>. <pub-id pub-id-type="doi">10.1146/annurev-ecolsys-110617-062453</pub-id><pub-id pub-id-type="pmid">32733173</pub-id></citation></ref>
<ref id="B33">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Murray</surname> <given-names>A. J.</given-names></name></person-group> (<year>2016</year>). <article-title>Energy metabolism and the high-altitude environment</article-title>. <source>Exp. Physiol.</source> <volume>101</volume>, <fpage>23</fpage>&#x02013;<lpage>27</lpage>. <pub-id pub-id-type="doi">10.1113/EP085317</pub-id><pub-id pub-id-type="pmid">26315373</pub-id></citation></ref>
<ref id="B34">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Nishida</surname> <given-names>A. H.</given-names></name> <name><surname>Ochman</surname> <given-names>H.</given-names></name></person-group> (<year>2018</year>). <article-title>Rates of gut microbiome divergence in mammals</article-title>. <source>Mol. Ecol.</source> <volume>27</volume>, <fpage>1884</fpage>&#x02013;<lpage>1897</lpage>. <pub-id pub-id-type="doi">10.1111/mec.14473</pub-id><pub-id pub-id-type="pmid">29290090</pub-id></citation></ref>
<ref id="B35">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Parks</surname> <given-names>D. H.</given-names></name> <name><surname>Tyson</surname> <given-names>G. W.</given-names></name> <name><surname>Hugenholtz</surname> <given-names>P.</given-names></name> <name><surname>Beiko</surname> <given-names>R. G.</given-names></name></person-group> (<year>2014</year>). <article-title>STAMP: statistical analysis of taxonomic and functional profiles</article-title>. <source>Bioinformatics</source> <volume>30</volume>, <fpage>3123</fpage>&#x02013;<lpage>3124</lpage>. <pub-id pub-id-type="doi">10.1093/bioinformatics/btu494</pub-id><pub-id pub-id-type="pmid">25061070</pub-id></citation></ref>
<ref id="B36">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Paulos</surname> <given-names>S.</given-names></name> <name><surname>Mateo</surname> <given-names>M.</given-names></name> <name><surname>De Lucio</surname> <given-names>A.</given-names></name> <name><surname>Hernandez-De Mingo</surname> <given-names>M.</given-names></name> <name><surname>Bailo</surname> <given-names>B.</given-names></name> <name><surname>Saugar</surname> <given-names>J. M.</given-names></name> <etal/></person-group>. (<year>2016</year>). <article-title>Evaluation of five commercial methods for the extraction and purification of DNA from human faecal samples for downstream molecular detection of the enteric protozoan parasites Cryptosporidium spp., Giardia duodenalis, and Entamoeba spp</article-title>. <source>J. Microbiol. Methods</source> <volume>127</volume>, <fpage>68</fpage>&#x02013;<lpage>73</lpage>. <pub-id pub-id-type="doi">10.1016/j.mimet.2016.05.020</pub-id><pub-id pub-id-type="pmid">27241828</pub-id></citation></ref>
<ref id="B37">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Praet</surname> <given-names>J.</given-names></name> <name><surname>Parmentier</surname> <given-names>A.</given-names></name> <name><surname>Schmid-Hempel</surname> <given-names>R.</given-names></name> <name><surname>Meeus</surname> <given-names>I.</given-names></name> <name><surname>Smagghe</surname> <given-names>G.</given-names></name> <name><surname>Vandamme</surname> <given-names>P.</given-names></name></person-group> (<year>2018</year>). <article-title>Large-scale cultivation of the bumblebee gut microbiota reveals an underestimated bacterial species diversity capable of pathogen inhibition</article-title>. <source>Environ. Microbiol.</source> <volume>20</volume>, <fpage>214</fpage>&#x02013;<lpage>227</lpage>. <pub-id pub-id-type="doi">10.1111/1462-2920.13973</pub-id><pub-id pub-id-type="pmid">29076622</pub-id></citation></ref>
<ref id="B38">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Quast</surname> <given-names>C.</given-names></name> <name><surname>Pruesse</surname> <given-names>E.</given-names></name> <name><surname>Yilmaz</surname> <given-names>P.</given-names></name> <name><surname>Gerken</surname> <given-names>J.</given-names></name> <name><surname>Schweer</surname> <given-names>T.</given-names></name> <name><surname>Yarza</surname> <given-names>P.</given-names></name> <etal/></person-group>. (<year>2013</year>). <article-title>The SILVA ribosomal RNA gene database project: improved data processing and web-based tools</article-title>. <source>Nucleic Acids Res.</source> <volume>41</volume>, <fpage>D590</fpage>&#x02013;<lpage>D596</lpage>. <pub-id pub-id-type="doi">10.1093/nar/gks1219</pub-id><pub-id pub-id-type="pmid">23193283</pub-id></citation></ref>
<ref id="B39">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Riek</surname> <given-names>A.</given-names></name> <name><surname>Stolzl</surname> <given-names>A.</given-names></name> <name><surname>Marquina Bernedo</surname> <given-names>R.</given-names></name> <name><surname>Ruf</surname> <given-names>T.</given-names></name> <name><surname>Arnold</surname> <given-names>W.</given-names></name> <name><surname>Hambly</surname> <given-names>C.</given-names></name> <etal/></person-group>. (<year>2019</year>). <article-title>Energy expenditure and body temperature variations in llamas living in the high andes of Peru</article-title>. <source>Sci. Rep.</source> <volume>9</volume>, <fpage>4037</fpage>. <pub-id pub-id-type="doi">10.1038/s41598-019-40576-9</pub-id><pub-id pub-id-type="pmid">30858417</pub-id></citation></ref>
<ref id="B40">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Smith</surname> <given-names>E. A.</given-names></name> <name><surname>Vuong</surname> <given-names>H. Q.</given-names></name> <name><surname>Miller</surname> <given-names>D. L.</given-names></name> <name><surname>Parish</surname> <given-names>A. J.</given-names></name> <name><surname>Mcfrederick</surname> <given-names>Q. S.</given-names></name> <name><surname>Newton</surname> <given-names>I. L. G.</given-names></name></person-group> (<year>2020</year>). <article-title>Draft genome sequences of four saccharibacter sp. strains isolated from native bees</article-title>. <source>Microbiol. Resour. Announc</source>. <volume>9</volume>, <fpage>e00022</fpage>-<lpage>20</lpage>. <pub-id pub-id-type="doi">10.1128/MRA.00022-20</pub-id><pub-id pub-id-type="pmid">32139579</pub-id></citation></ref>
<ref id="B41">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Storz</surname> <given-names>J. F.</given-names></name> <name><surname>Scott</surname> <given-names>G. R.</given-names></name> <name><surname>Cheviron</surname> <given-names>Z. A.</given-names></name></person-group> (<year>2010</year>). <article-title>Phenotypic plasticity and genetic adaptation to high-altitude hypoxia in vertebrates</article-title>. <source>J. Exp. Biol.</source> <volume>213</volume>, <fpage>4125</fpage>&#x02013;<lpage>4136</lpage>. <pub-id pub-id-type="doi">10.1242/jeb.048181</pub-id><pub-id pub-id-type="pmid">33301891</pub-id></citation></ref>
<ref id="B42">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Weber</surname> <given-names>R. E.</given-names></name></person-group> (<year>2007</year>). <article-title>High-altitude adaptations in vertebrate hemoglobins</article-title>. <source>Respir. Physiol. Neurobiol.</source> <volume>158</volume>, <fpage>132</fpage>&#x02013;<lpage>142</lpage>. <pub-id pub-id-type="doi">10.1016/j.resp.2007.05.001</pub-id><pub-id pub-id-type="pmid">17561448</pub-id></citation></ref>
<ref id="B43">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Williams</surname> <given-names>P. H.</given-names></name> <name><surname>Huang</surname> <given-names>J.</given-names></name> <name><surname>Rasmont</surname> <given-names>P.</given-names></name> <name><surname>An</surname> <given-names>J.</given-names></name></person-group> (<year>2016</year>). <article-title>Early-diverging bumblebees from across the roof of the world: the high-mountain subgenus Mendacibombus revised from species&#x00027; gene coalescents and morphology (Hymenoptera, Apidae)</article-title>. <source>Zootaxa</source> <volume>4204</volume>, <fpage>1</fpage>&#x02013;<lpage>72</lpage>. <pub-id pub-id-type="doi">10.11646/zootaxa.4204.1.1</pub-id><pub-id pub-id-type="pmid">27988613</pub-id></citation></ref>
<ref id="B44">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Wu</surname> <given-names>J.</given-names></name> <name><surname>An</surname> <given-names>J. D.</given-names></name> <name><surname>Yao</surname> <given-names>J.</given-names></name> <name><surname>Huang</surname> <given-names>J. X.</given-names></name> <name><surname>Feng</surname> <given-names>X. Q.</given-names></name></person-group> (<year>2009</year>). <article-title><italic>Bombus fauna</italic> (Hymenoptera, Apidae) in Hebei, China</article-title>. <source>Acta. Zootaxonomica Sinica</source>, <volume>34</volume>, <fpage>87</fpage>&#x02013;<lpage>97</lpage>.</citation>
</ref>
<ref id="B45">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Xu</surname> <given-names>J.</given-names></name> <name><surname>James</surname> <given-names>R. R.</given-names></name></person-group> (<year>2012</year>). <article-title>Temperature stress affects the expression of immune response genes in the alfalfa leafcutting bee, Megachile rotundata</article-title>. <source>Insect Mol. Biol.</source> <volume>21</volume>, <fpage>269</fpage>&#x02013;<lpage>280</lpage>. <pub-id pub-id-type="doi">10.1111/j.1365-2583.2012.01133.x</pub-id><pub-id pub-id-type="pmid">22356318</pub-id></citation></ref>
<ref id="B46">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Yang</surname> <given-names>J.</given-names></name> <name><surname>Jin</surname> <given-names>Z. B.</given-names></name> <name><surname>Chen</surname> <given-names>J.</given-names></name> <name><surname>Huang</surname> <given-names>X. F.</given-names></name> <name><surname>Li</surname> <given-names>X. M.</given-names></name> <name><surname>Liang</surname> <given-names>Y. B.</given-names></name> <etal/></person-group>. (<year>2017</year>). <article-title>Genetic signatures of high-altitude adaptation in Tibetans</article-title>. <source>Proc. Natl. Acad. Sci. USA.</source> <volume>114</volume>, <fpage>4189</fpage>&#x02013;<lpage>4194</lpage>. <pub-id pub-id-type="doi">10.1073/pnas.1617042114</pub-id><pub-id pub-id-type="pmid">28373541</pub-id></citation></ref>
<ref id="B47">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Yu</surname> <given-names>X.</given-names></name> <name><surname>Wei</surname> <given-names>P.</given-names></name> <name><surname>Zhao</surname> <given-names>S.</given-names></name> <name><surname>Chen</surname> <given-names>Z.</given-names></name> <name><surname>Li</surname> <given-names>X.</given-names></name> <name><surname>Zhang</surname> <given-names>W.</given-names></name> <etal/></person-group>. (<year>2023</year>). <article-title>Population transcriptomics uncover the relative roles of positive selection and differential expression in Batrachium bungei adaptation to the Qinghai-Tibetan plateau</article-title>. <source>Plant Cell Rep</source>. <volume>42</volume>, <fpage>879</fpage>&#x02013;<lpage>893</lpage>. <pub-id pub-id-type="doi">10.1007/s00299-023-03005-w</pub-id><pub-id pub-id-type="pmid">36973418</pub-id></citation></ref>
<ref id="B48">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Zhang</surname> <given-names>J.</given-names></name> <name><surname>Marshall</surname> <given-names>K. E.</given-names></name> <name><surname>Westwood</surname> <given-names>J. T.</given-names></name> <name><surname>Clark</surname> <given-names>M. S.</given-names></name> <name><surname>Sinclair</surname> <given-names>B. J.</given-names></name></person-group> (<year>2011</year>). <article-title>Divergent transcriptomic responses to repeated and single cold exposures in Drosophila melanogaster</article-title>. <source>J. Exp. Biol.</source> <volume>214</volume>, <fpage>4021</fpage>&#x02013;<lpage>4029</lpage>. <pub-id pub-id-type="doi">10.1242/jeb.059535</pub-id><pub-id pub-id-type="pmid">22071194</pub-id></citation></ref>
<ref id="B49">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Zhang</surname> <given-names>Z. J.</given-names></name> <name><surname>Huang</surname> <given-names>M. F.</given-names></name> <name><surname>Qiu</surname> <given-names>L. F.</given-names></name> <name><surname>Song</surname> <given-names>R. H.</given-names></name> <name><surname>Zhang</surname> <given-names>Z. X.</given-names></name> <name><surname>Ding</surname> <given-names>Y. W.</given-names></name> <etal/></person-group>. (<year>2021</year>). <article-title>Diversity and functional analysis of Chinese bumblebee gut microbiota reveal the metabolic niche and antibiotic resistance variation of Gilliamella</article-title>. <source>Insect Sci.</source> <volume>28</volume>, <fpage>302</fpage>&#x02013;<lpage>314</lpage>. <pub-id pub-id-type="doi">10.1111/1744-7917.12770</pub-id><pub-id pub-id-type="pmid">32101381</pub-id></citation></ref>
<ref id="B50">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Zhang</surname> <given-names>Z. J.</given-names></name> <name><surname>Zheng</surname> <given-names>H.</given-names></name></person-group> (<year>2022</year>). <article-title>Bumblebees with the socially transmitted microbiome: a novel model organism for gut microbiota research</article-title>. <source>Insect Sci.</source> <volume>29</volume>, <fpage>958</fpage>&#x02013;<lpage>976</lpage>. <pub-id pub-id-type="doi">10.1111/1744-7917.13040</pub-id><pub-id pub-id-type="pmid">35567381</pub-id></citation></ref>
<ref id="B51">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Zheng</surname> <given-names>H.</given-names></name> <name><surname>Nishid,A</surname> <given-names>A.</given-names></name> <name><surname>Kwong</surname> <given-names>W. K.</given-names></name> <name><surname>Koch</surname> <given-names>H.</given-names></name> <name><surname>Engel</surname> <given-names>P.</given-names></name> <name><surname>Steele</surname> <given-names>M. I.</given-names></name> <etal/></person-group>. (<year>2016</year>). <article-title>Metabolism of toxic sugars by strains of the bee gut symbiont gilliamella apicola</article-title>. <source>mBio</source>, <volume>7</volume>, <fpage>e01326</fpage>-<lpage>16</lpage>. <pub-id pub-id-type="doi">10.1128/mBio.01326-16</pub-id><pub-id pub-id-type="pmid">27803186</pub-id></citation></ref>
</ref-list> 
</back>
</article> 