The adaptation of bumblebees to extremely high elevation associated with their gut microbiota.

Bumblebees are among the most abundant and important pollinators for sub-alpine and alpine flowering plant species in the Northern Hemisphere, but little is known about their adaptations to high elevations. In this article, we focused on two bumblebee species, Bombus friseanus and Bombus prshewalskyi, and their respective gut microbiota. The two species, distributed through the Hengduan Mountains of southwestern China, show species replacement at different elevations. We performed genome sequencing based on 20 worker bee samples of each species. Applying evolutionary population genetics and metagenomic approaches, we detected genes under selection and analyzed functional pathways between bumblebees and their gut microbes. We found clear genetic differentiation between the two host species and significant differences in their microbiota. Species replacement occurred in both hosts and their bacteria (Snodgrassella) with an increase in elevation. These extremely high-elevation bumblebees show evidence of positive selection related to diverse biological processes. Positively selected genes involved in host immune systems probably contributed to gut microbiota changes, while the butyrate generated by gut microbiota may influence both host energy metabolism and immune systems. This suggests a close association between the genomes of the host species and their microbiomes based on some degree of natural selection.IMPORTANCETwo closely related and dominant bumblebee species, distributed at different elevations through the Hengduan Mountains of southwestern China, showed a clear genomic signature of adaptation to elevation at the molecular level and significant differences in their respective microbiota. Species replacement occurred in both hosts and their bacteria (Snodgrassella) with an increase in elevation. Bumblebees' adaptations to higher elevations are closely associated with their gut microbiota through two biological processes: energy metabolism and immune response. Information allowing us to understand the adaptive mechanisms of species to extreme conditions is implicit if we are to conserve them as their environments change.

Dioxin-like polychlorinated biphenyl 126 (PCB126) disrupts gut microbiota-host metabolic dysfunction in mice via aryl hydrocarbon receptor activation.

Exposure to environmental pollutants, including dioxin-like pollutants, can cause numerous health issues. A common exposure route to pollutants is through contaminated foods, and thus the gastrointestinal system and gut microbiota are often exposed to high amounts of pollutants. Multiple studies have focused on the imbalance in intestinal microbiota composition caused by dioxin-like pollutants. Here, we examined the effects of polychlorinated biphenyl 126 (PCB126) on the composition and functions of gut microbes through metagenomic sequencing, and explored the correlations between microflora dysbiosis and aryl hydrocarbon receptor (AHR) signaling. Adult male wild-type and Ahr-/- mice with a C57BL/6 background were weekly exposed to 50 µg/kg body weight of PCB126 for 8 weeks. Results showed that PCB126 had the opposite effect on gut microbiota composition and diversity in the wild-type and Ahr-/- mice. Functional prediction found that PCB126 exposure mainly altered carbon metabolism and signal regulatory pathways in wild-type mice but impacted DNA replication and lipopolysaccharide biosynthesis in Ahr-/- mice. In wild-type mice, PCB126 exposure induced liver injury, decreased serum lipid content, and delayed gastrointestinal motility, which were significantly correlated to several specific bacterial taxa, such as Helicobacter. Following AHR knockout, however, the holistic effects of PCB126 on the host were lessened or abolished. These results suggest that PCB126 may disrupt host metabolism and gut microbiota dynamics via AHR activation. Overall, our findings provide new insight into the complex interactions between host metabolism and gut microbiota, which may contribute to grouped assessment of environmental pollutants in the future.

Diet drives convergent evolution of gut microbiomes in bamboo-eating species.

Gut microbiota plays a critical role in host physiology and health. The coevolution between the host and its gut microbes facilitates animal adaptation to its specific ecological niche. Multiple factors such as host diet and phylogeny modulate the structure and function of gut microbiota. However, the relative contribution of each factor in shaping the structure of gut microbiota remains unclear. The giant (Ailuropoda melanoleuca) and red (Ailurus styani) pandas belong to different families of order Carnivora. They have evolved as obligate bamboo-feeders and can be used as a model system for studying the gut microbiome convergent evolution. Here, we compare the structure and function of gut microbiota of the two pandas with their carnivorous relatives using 16S rRNA and metagenome sequencing. We found that both panda species share more similarities in their gut microbiota structure with each other than each species shares with its carnivorous relatives. This indicates that the specialized herbivorous diet rather than host phylogeny is the dominant driver of gut microbiome convergence within Arctoidea. Metagenomic analysis revealed that the symbiotic gut microbiota of both pandas possesses a high level of starch and sucrose metabolism and vitamin B12 biosynthesis. These findings suggest a diet-driven convergence of gut microbiomes and provide new insight into host-microbiota coevolution of these endangered species.

A New Alignment-Free Whole Metagenome Comparison Tool and Its Application on Gut Microbiomes of Wild Giant Pandas.

The comparison of metagenomes is crucial for studying the relationship between microbial communities and environmental factors. One recently published alignment-free whole metagenome comparison method based on k-mer frequencies, Libra, showed higher resolutions than the present fastest method, Mash, on whole metagenomic sequencing reads, but it did not perform as well on the assembled contigs. Here, we developed a new alignment-free tool, KmerFreqCalc, for the comparison of the whole metagenomic data, which first calculated the frequencies of both forward and reverse complementary sequences of k-mers like Mash and then computed the cosine distance between the samples based on k-mer frequency vectors like Libra. We applied KmerFreqCalc on the assembled contigs of the gut microbiomes of wild giant pandas and compared the results to Libra and Mash. The results indicated that KmerFreqCalc was able to detect the subtle difference between giant panda samples caused by seasonal diet change, showing better clustering than Libra and Mash. Therefore, KmerFreqCalc has high resolution and accuracy for assembled contigs, being very suitable for comparison of samples with low dissimilarity.

Detection and typing of human-infecting influenza viruses in China by using a multiplex DNA biochip assay.

Rapid identification of the infections of specific subtypes of influenza viruses is critical for patient treatment and pandemic control. Here we report the application of multiplex reverse transcription polymerase chain reaction (RT-PCR) coupled with membrane-based DNA biochip to the detection and discrimination of the type (A and B) and subtype (human H1N1, human H3N2, avian H5N1 and avian H7N9) of influenza viruses in circulation in China. A multiplex one-step RT-PCR assay was designed to simultaneously amplify the HA and NA genes of the four subtypes of influenza A viruses and NS genes to discriminate type A and B viruses. PCR products were analyzed by a membrane-based biochip. The analytical sensitivity of the assay was determined at a range of 2-100 copies/reactions for each of the gene transcripts. Eighty one clinical samples, containing 66 positive samples with evident seasonal influenza virus infections, were tested, which gives the clinical sensitivity and specificity of 95.5% and 100% respectively. For the avian influenza samples, 3 out of 4 H5N1 samples and 2 out of 2 H7N9 avian samples were correctly identified. We argue this method could allow a rapid, reliable and inexpensive detection and differentiation of human-infecting influenza viruses.

Integrating computational modeling and functional assays to decipher the structure-function relationship of influenza virus PB1 protein.

The influenza virus PB1 protein is the core subunit of the heterotrimeric polymerase complex (PA, PB1 and PB2) in which PB1 is responsible for catalyzing RNA polymerization and binding to the viral RNA promoter. Among the three subunits, PB1 is the least known subunit so far in terms of its structural information. In this work, by integrating template-based structural modeling approach with all known sequence and functional information about the PB1 protein, we constructed a modeled structure of PB1. Based on this model, we performed mutagenesis analysis for the key residues that constitute the RNA template binding and catalytic (TBC) channel in an RNP reconstitution system. The results correlated well with the model and further identified new residues of PB1 that are critical for RNA synthesis. Moreover, we derived 5 peptides from the sequence of PB1 that form the TBC channel and 4 of them can inhibit the viral RNA polymerase activity. Interestingly, we found that one of them named PB1(491-515) can inhibit influenza virus replication by disrupting viral RNA promoter binding activity of polymerase. Therefore, this study has not only deepened our understanding of structure-function relationship of PB1, but also promoted the development of novel therapeutics against influenza virus.

Coiled-coil networking shapes cell molecular machinery.

The highly abundant α-helical coiled-coil motif not only mediates crucial protein-protein interactions in the cell but is also an attractive scaffold in synthetic biology and material science and a potential target for disease intervention. Therefore a systematic understanding of the coiled-coil interactions (CCIs) at the organismal level would help unravel the full spectrum of the biological function of this interaction motif and facilitate its application in therapeutics. We report the first identified genome-wide CCI network in Saccharomyces cerevisiae, which consists of 3495 pair-wise interactions among 598 predicted coiled-coil regions. Computational analysis revealed that the CCI network is specifically and functionally organized and extensively involved in the organization of cell machinery. We further show that CCIs play a critical role in the assembly of the kinetochore, and disruption of the CCI network leads to defects in kinetochore assembly and cell division. The CCI network identified in this study is a valuable resource for systematic characterization of coiled coils in the shaping and regulation of a host of cellular machineries and provides a basis for the utilization of coiled coils as domain-based probes for network perturbation and pharmacological applications.

Homologs of eukaryotic Ras superfamily proteins in prokaryotes and their novel phylogenetic correlation with their eukaryotic analogs.

Ras superfamily proteins are key regulators in a wide variety of cellular processes. Previously, they were considered to be specific to eukaryotes, and MglA, a group of obviously different prokaryotic proteins, were recognized as their only prokaryotic analogs or even ancestors. Here, taking advantage of quite a current accumulation of prokaryotic genomic databases, we have investigated the existence and taxonomic distribution of Ras superfamily protein homologs in a much wider prokaryotic range, and analyzed their phylogenetic correlation with their eukaryotic analogs. Thirteen unambiguous prokaryotic homologs, which possess the GDP/GTP-binding domain with all the five characteristic motifs of their eukaryotic analogs, were identified in 12 eubacteria and one archaebacterium, respectively. In some other archaebacteria, including four methanogenic archaebacteria and three Thermoplasmales, homologs were also found, but with the GDP/GTP-binding domains not containing all the five characteristic motifs. Many more MglA orthologs were identified than in previous studies mainly in delta-proteobacteria, and all were shown to have common unique features distinct from the Ras superfamily proteins. Our phylogenetic analysis indicated eukaryotic Rab, Ran, Ras, and Rho families have the closest phylogenetic correlation with the 13 unambiguous prokaryotic homologs, whereas the other three eukaryotic protein families (SRbeta, Sar1, and Arf) branch separately from them, but have a relatively close relationship with the methanogenic archaebacterial homologs and MglA. Although homologs were identified in a relative minority of prokaryotes with genomic databases, their presence in a relatively wide variety of lineages, their unique sequence characters distinct from those of eukaryotic analogs, and the topology of our phylogenetic tree altogether do not support their origin from eukaryotes as a result of lateral gene transfer. Therefore, we argue that Ras superfamily proteins might have already emerged at least in some prokaryotic lineages, and that the seven eukaryotic protein families of the Ras superfamily may have two independent prokaryotic origins, probably reflecting the 'fusion' evolutionary history of the eukaryotic cell.

Phylogenetic positions of several amitochondriate protozoa--evidence from phylogenetic analysis of DNA topoisomerase II.

Several groups of parasitic protozoa, as represented by Giardia, Trichomonas, Entamoeba and Microsporida, were once widely considered to be the most primitive extant eukaryotic group--Archezoa. The main evidence for this is their 'lacking mitochondria' and possessing some other primitive features between prokaryotes and eukaryotes, and being basal to all eukaryotes with mitochondria in phylogenies inferred from many molecules. Some authors even proposed that these organisms diverged before the endosymbiotic origin of mitochondria within eukaryotes. This view was once considered to be very significant to the study of origin and evolution of eukaryotic cells (eukaryotes). However, in recent years this has been challenged by accumulating evidence from new studies. Here the sequences of DNA topoisomerase II in G lamblia, T. vaginalis and E. histolytica were identified first by PCR and sequencing, then combining with the sequence data of the microsporidia Encephalitozoon cunicul and other eukaryotic groups of different evolutionary positions from GenBank, phylogenetic trees were constructed by various methods to investigate the evolutionary positions of these amitochondriate protozoa. Our results showed that since the characteristics of DNA topoisomerase II make it avoid the defect of 'long-branch attraction' appearing in the previous phylogenetic analyses, our trees can not only reflect effectively the relationship of different major eukaryotic groups, which is widely accepted, but also reveal phylogenetic positions for these amitochondriate protozoa, which is different from the previous phylogenetic trees. They are not the earliest-branching eukaryotes, but diverged after some mitochondriate organisms such as kinetoplastids and mycetozoan; they are not a united group but occupy different phylogenetic positions. Combining with the recent cytological findings of mitochondria-like organelles in them, we think that though some of them (e.g. diplomonads, as represented by Giardia) may occupy a very low evolutionary position, generally these organisms are not as extremely primitive as was thought before; they should be polyphyletic groups diverging after the endosymbiotic origin of mitochondrion to adapt themselves to anaerobic parasitic life.