Gut microbial diversity across a contact zone for California voles: Implications for lineage divergence of hosts and mitonuclear mismatch in the assembly of the mammalian gut microbiome

Gut microbial diversity is thought to reflect the co‐evolution of microbes and their hosts as well as current host‐specific attributes such as genetic background and environmental setting. To explore interactions among these parameters, we characterized variation in gut microbiome composition of California voles (Microtus californicus) across a contact zone between two recently diverged lineages of this species. Because this contact zone contains individuals with mismatched mitochondrial‐nuclear genomes (cybrids), it provides an important opportunity to explore how different components of the genotype contribute to gut microbial diversity. Analyses of bacterial 16S rRNA sequences and joint species distribution modelling revealed that host genotypes and genetic differentiation among host populations together explained more than 50% of microbial community variation across our sampling transect. The ranked importance (most to least) of factors contributing to gut microbial diversity in our study populations were: genome‐wide population differentiation, local environmental conditions, and host genotypes. However, differences in microbial communities among vole populations (β‐diversity) did not follow patterns of lineage divergence (i.e., phylosymbiosis). Instead, among‐population variation was best explained by the spatial distribution of hosts, as expected if the environment is a primary source of gut microbial diversity (i.e., dispersal limitation hypothesis). Across the contact zone, several bacterial taxa differed in relative abundance between the two parental lineages as well as among individuals with mismatched mitochondrial and nuclear genomes. Thus, genetic divergence among host lineages and mitonuclear genomic mismatches may also contribute to microbial diversity by altering interactions between host genomes and gut microbiota (i.e., hologenome speciation hypothesis).     

The origin of malarial parasites in orangutans

Background

Recent findings of Plasmodium in African apes have changed our perspectives on the evolution of malarial parasites in hominids. However, phylogenetic analyses of primate malarias are still missing information from Southeast Asian apes. In this study, we report molecular data for a malaria parasite lineage found in orangutans.

Methodology/Principal Findings

We screened twenty-four blood samples from Pongo pygmaeus (Kalimantan, Indonesia) for Plasmodium parasites by PCR. For all the malaria positive orangutan samples, parasite mitochondrial genomes (mtDNA) and two antigens: merozoite surface protein 1 42 kDa (MSP-1₄₂) and circumsporozoite protein gene (CSP) were amplified, cloned, and sequenced. Fifteen orangutans tested positive and yielded 5 distinct mitochondrial haplotypes not previously found. The haplotypes detected exhibited low genetic divergence among them, indicating that they belong to one species. We report phylogenetic analyses using mitochondrial genomes, MSP-1₄₂ and CSP. We found that the orangutan malaria parasite lineage was part of a monophyletic group that includes all the known non-human primate malaria parasites found in Southeast Asia; specifically, it shares a recent common ancestor with P. inui (a macaque parasite) and P. hylobati (a gibbon parasite) suggesting that this lineage originated as a result of a host switch. The genetic diversity of MSP-1₄₂ in orangutans seems to be under negative selection. This result is similar to previous findings in non-human primate malarias closely related to P. vivax. As has been previously observed in the other Plasmodium species found in non-human primates, the CSP shows high polymorphism in the number of repeats. However, it has clearly distinctive motifs from those previously found in other malarial parasites.

Conclusion

The evidence available from Asian apes indicates that these parasites originated independently from those found in Africa, likely as the result of host switches from other non-human primates.     

Genetic markers in Andean Puya species (Bromeliaceae) with implications on plastome evolution and phylogeny

The Andean plant endemic Puya is a striking example of recent and rapid diversification from central Chile to the northern Andes, tracking mountain uplift. This study generated 12 complete plastomes representing nine Puya species and compared them to five published plastomes for their features, genomic evolution, and phylogeny. The total size of the Puya plastomes ranged from 159,542 to 159,839 bp with 37.3%–37.4% GC content. The Puya plastomes were highly conserved in organization and structure with a typical quadripartite genome structure. Each of the 17 consensus plastomes harbored 133 genes, including 87 protein‐coding genes, 38 tRNA (transfer RNA) genes, and eight rRNA (ribosomal RNA) genes; we found 69–78 tandem repeats, 45–60 SSRs (simple sequence repeats), and 8–22 repeat structures among 13 species. Four protein‐coding genes were identified under positive site‐specific selection in Puya. The complete plastomes and hypervariable regions collectively provided pronounced species discrimination in Puya and a practical tool for future phylogenetic studies. The reconstructed phylogeny and estimated divergence time for the lineage suggest that the diversification of Puya is related to Andean orogeny and Pleistocene climatic oscillations. This study provides plastome resources for species delimitation and novel phylogenetic and biogeographic studies.     

Phylogenetic analysis and evolutionary studies of plant carotenoid cleavage dioxygenase gene

The oxidative breakdown of carotenoid evidences the formation of apocarotenoids through carotenoid cleavage dioxygenases (CCDs). Numerous CCDs and apocarotenoids have been identified and characterized in plants. Using available sequence data, a study was performed to investigate the phylogenetic relationship among CCD genes and to statistically estimate the sequence conservation and functional divergence. In total, 77 genes were identified from 39 species belonging to 21 families. Our result of phylogenetic analysis indicated the existence of well-conserved subfamilies. Moreover, comparative genomic analysis showed that the gene structures of the CCDs were highly conserved across some different lineage species. Through functional divergence analysis, a substantial divergence was found between CCD subfamilies. In addition, examination of the site-specific profile revealed the critical amino acid residues accounting for functional divergence. This study mainly focused on the evolution of CCD genes and their functional divergence which may deliver an initial step for further experimental verifications.     

Cryptic species and co‐diversification in sand scorpions from the Karakum and Kyzylkum deserts of Central Asia

Little is known about species diversification within the deserts of Central Asia. For example, the degree of lineage divergence and timing of population differentiation, as well as potential biogeographic barriers driving diversification, are nearly unknown. Here, we analysed a multi‐locus data set for a widespread sand scorpion (Mesobuthus gorelovi) to evaluate cryptic species diversity and phylogeographic patterns across the Karakum and Kyzylkum deserts. We also combined these data with previously published sequence data to test for a signal of co‐diversification. A consensus species delimitation approach indicated that the widespread M. gorelovi is likely composed of up to five distinct species that began to diversify at the Miocene–Pliocene boundary. We observed shared patterns of lineage divergence across the Amu Darya River region in three scorpion taxa and found support for a shared history of assemblage diversification across this biogeographic barrier. Thus, major river systems appear to facilitate diversification among desert scorpions.     

Phylogenomics of the Maverick Virus-Like Mobile Genetic Elements of Vertebrates

Mavericks are virus-like mobile genetic elements found in the genomes of eukaryotes. Although Mavericks encode capsid morphogenesis homologs, their viral particles have not been observed. Here, we provide new evidence supporting the viral nature of Mavericks and the potential existence of virions. To this end, we conducted a phylogenomic analysis of Mavericks in hundreds of vertebrate genomes, discovering 134 elements with an intact coding capacity in 17 host species. We reveal an extensive genomic fossil record in 143 species and date three groups of elements to the Late Cretaceous. Bayesian phylogenetic analysis using genomic fossil orthologs suggests that Mavericks have infected osteichthyans for ∼419 My. They have undergone frequent cross-species transmissions in cyprinid fish and all core genes are subject to strong purifying selection. We conclude that vertebrate Mavericks form an ancient lineage of aquatic dsDNA viruses which are probably still functional in some vertebrate lineages.     

Hidden Diversity in Honey Bee Gut Symbionts Detected by Single-Cell Genomics

Microbial communities in animal guts are composed of diverse, specialized bacterial species, but little is known about how gut bacteria diversify to produce genetically and ecologically distinct entities. The gut microbiota of the honey bee, Apis mellifera, presents a useful model, because it consists of a small number of characteristic bacterial species, each showing signs of diversification. Here, we used single-cell genomics to study the variation within two species of the bee gut microbiota: Gilliamella apicola and Snodgrassella alvi. For both species, our analyses revealed extensive variation in intraspecific divergence of protein-coding genes but uniformly high levels of 16S rRNA similarity. In both species, the divergence of 16S rRNA loci appears to have been curtailed by frequent recombination within populations, while other genomic regions have continuously diverged. Furthermore, gene repertoires differ markedly among strains in both species, implying distinct metabolic capabilities. Our results show that, despite minimal divergence at 16S rRNA genes, in situ diversification occurs within gut communities and generates bacterial lineages with distinct ecological niches. Therefore, important dimensions of microbial diversity are not evident from analyses of 16S rRNA, and single cell genomics has potential to elucidate processes of bacterial diversification.     

Multiple Single-Cell Genomes Provide Insight into Functions of Uncultured Deltaproteobacteria in the Human Oral Cavity

Despite a long history of investigation, many bacteria associated with the human oral cavity have yet to be cultured. Studies that correlate the presence or abundance of uncultured species with oral health or disease highlight the importance of these community members. Thus, we sequenced several single-cell genomic amplicons from Desulfobulbus and Desulfovibrio (class Deltaproteobacteria) to better understand their function within the human oral community and their association with periodontitis, as well as other systemic diseases. Genomic data from oral Desulfobulbus and Desulfovibrio species were compared to other available deltaproteobacterial genomes, including from a subset of host-associated species. While both groups share a large number of genes with other environmental Deltaproteobacteria genomes, they encode a wide array of unique genes that appear to function in survival in a host environment. Many of these genes are similar to virulence and host adaptation factors of known human pathogens, suggesting that the oral Deltaproteobacteria have the potential to play a role in the etiology of periodontal disease.     

Genomic and genetic analyses of diversity and plant interactions of Pseudomonas fluorescens

Background

Pseudomonas fluorescens are common soil bacteria that can improve plant health through nutrient cycling, pathogen antagonism and induction of plant defenses. The genome sequences of strains SBW25 and Pf0-1 were determined and compared to each other and with P. fluorescens Pf-5. A functional genomic in vivo expression technology (IVET) screen provided insight into genes used by P. fluorescens in its natural environment and an improved understanding of the ecological significance of diversity within this species.

Results

Comparisons of three P. fluorescens genomes (SBW25, Pf0-1, Pf-5) revealed considerable divergence: 61% of genes are shared, the majority located near the replication origin. Phylogenetic and average amino acid identity analyses showed a low overall relationship. A functional screen of SBW25 defined 125 plant-induced genes including a range of functions specific to the plant environment. Orthologues of 83 of these exist in Pf0-1 and Pf-5, with 73 shared by both strains. The P. fluorescens genomes carry numerous complex repetitive DNA sequences, some resembling Miniature Inverted-repeat Transposable Elements (MITEs). In SBW25, repeat density and distribution revealed ''repeat deserts'' lacking repeats, covering approximately 40% of the genome.

Conclusions

P. fluorescens genomes are highly diverse. Strain-specific regions around the replication terminus suggest genome compartmentalization. The genomic heterogeneity among the three strains is reminiscent of a species complex rather than a single species. That 42% of plant-inducible genes were not shared by all strains reinforces this conclusion and shows that ecological success requires specialized and core functions. The diversity also indicates the significant size of genetic information within the Pseudomonas pan genome.     

Uncovering Cryptic Diversity in Two Amoebozoan Species Using Complete Mitochondrial Genome Sequences

The Amoebozoa are a major eukaryotic lineage that encompasses a wide range of amoeboid organisms. The group is understudied from a systematic perspective: molecular tools have only been applied in the last 15 yr. Hence, there is an undersampling of both genes and taxa in the group especially compared to plants, animals, and fungi. Here, we present the complete mitochondrial genomes of two ubiquitous and abundant morpho‐species (Acanthamoeba castellanii and Vermamoeba vermiformis). Both have mitochondrial genomes of close relatives previously available, enabling insights into recent divergences at a genomic scale, while simultaneously offering comparisons with divergence estimates obtained from traditionally used single genes, SSU rDNA and cox1. The newly sequenced mt genomes are significantly divergent from their previously sequenced conspecifics (A. castellannii 16.4% divergence at nucleotide level and 10.4% amino acid; V. vermiformis 21.6% and 13.1%, respectively), while divergence at the small subunit ribosomal DNA is below 1% within both species. Morphological analyses determined that these lineages are indistinguishable from their previously sequenced counterparts. Phylogenetic reconstructions using 26 mt genes also indicate a level of divergence that is comparable to divergence among species, while reconstructions using the small subunit ribosomal DNA (SSU rDNA) do not. In addition, we demonstrate that between closely related taxa, there are high levels of synteny, which can be explored for primer design to obtain larger fragments than the traditional barcoding genes. We conclude that, although most systematic work has relied on SSU, this gene alone can severely underestimate diversity. Thus, we suggest that the mt genome emerges as an alternative for unraveling the lower level phylogenetic relationships of Amoebozoa.     

Conserved features of cohesin binding along fission yeast chromosomes

Background

Pseudomonas fluorescens are common soil bacteria that can improve plant health through nutrient cycling, pathogen antagonism and induction of plant defenses. The genome sequences of strains SBW25 and Pf0-1 were determined and compared to each other and with P. fluorescens Pf-5. A functional genomic in vivo expression technology (IVET) screen provided insight into genes used by P. fluorescens in its natural environment and an improved understanding of the ecological significance of diversity within this species.

Results

Comparisons of three P. fluorescens genomes (SBW25, Pf0-1, Pf-5) revealed considerable divergence: 61% of genes are shared, the majority located near the replication origin. Phylogenetic and average amino acid identity analyses showed a low overall relationship. A functional screen of SBW25 defined 125 plant-induced genes including a range of functions specific to the plant environment. Orthologues of 83 of these exist in Pf0-1 and Pf-5, with 73 shared by both strains. The P. fluorescens genomes carry numerous complex repetitive DNA sequences, some resembling Miniature Inverted-repeat Transposable Elements (MITEs). In SBW25, repeat density and distribution revealed 'repeat deserts' lacking repeats, covering approximately 40% of the genome.

Conclusions

P. fluorescens genomes are highly diverse. Strain-specific regions around the replication terminus suggest genome compartmentalization. The genomic heterogeneity among the three strains is reminiscent of a species complex rather than a single species. That 42% of plant-inducible genes were not shared by all strains reinforces this conclusion and shows that ecological success requires specialized and core functions. The diversity also indicates the significant size of genetic information within the Pseudomonas pan genome.     

First reported chloroplast genome sequence of <Emphasis Type="Italic">Punica granatum</Emphasis> (cultivar Helow) from Jabal Al-Akhdar,Oman: phylogenetic comparative assortment with <Emphasis Type="Italic">Lagerstroemia</Emphasis>

Pomegranate (Punica granatum L.) is one of the oldest known edible fruits. It has grown in popularity and is a profitable fruit crop due to its attractive features including a bright red appearance and its biological activities. Scientific exploration of the genetics and evolution of these beneficial traits has been hampered by limited genomic information. In this study, we sequenced the complete chloroplast (cp) genome of the native P. granatum (cultivar Helow) cultivated in the mountains of Jabal Al-Akhdar, Oman. The results revealed a P. granatum cp genome length of 158,630 bp, characterized by a relatively conserved structure containing 2 inverted repeat regions of 25,466 bp, an 18,686 bp small single copy regions, and an 89,015 bp large single copy region. The 86 protein-coding genes included 37 transfer RNA genes and 8 ribosomal RNA genes. Comparison of the P. granatum whole cp genome with seven Lagerstroemia species revealed an overall high degree of sequence similarity with divergence among intergenic spacers. The location, distribution, and divergence of repeat sequences and shared genes of the Punica and Lagerstroemia species were highly similar. Analyses of nucleotide substitution, insertion/deletions, and highly variable regions in these cp genomes identified potential plastid markers for taxonomic and phylogenetic studies in Myrtales. A phylogenetic study of the cp genomes and 76 shared coding regions generated similar cladograms. The complete cp genome of P. granatum will aid in taxonomical studies of the family Lythraceae.     

Prevotella copri ameliorates cholestasis and liver fibrosis in primary sclerosing cholangitis by enhancing the FXR signalling pathway

Primary sclerosing cholangitis (PSC) is a chronic cholestatic liver disease characterized by bile duct inflammation, fibrosis, bile acid (BA) metabolism disorders and gut microbiota dysbiosis. At present, the aetiology and pathogenesis of PSC are not clear, and there is no specific or effective treatment available. Therefore, new research perspectives are needed to explore effective methods to treat PSC and improve symptoms. The intestinal microbiota of patients with PSC is known to be significantly different from that of healthy people. By comparing differentially abundant bacterial genera in PSC patients, it was found that the abundance of Prevotella copri (P. copri) was significantly decreased, suggesting that this species may have a protective effect against PSC disease. Therefore, comprehensively exploring the role and possible function of P. copri in the disease process is worthwhile. In this study, a PSC mouse model was established by feeding mice a customized diet supplemented with 0.1% (w/w) 3,5-diethoxycarbonyl-1,4-dihydrocollidine (DDC) for one week, and the abundance of P. copri was confirmed to be decreased in this model. Previous studies in patients and animal models have demonstrated that gut microbiota intervention is an acceptable treatment for some diseases. We found that intervention with P. copri could significantly improve cholestasis and liver fibrosis by enhancing the FXR-related signalling pathway in PSC mice. Together, through the overall effect of P. copri on intestinal microbiota structure and its association with BAs, we speculate that P. copri intervention might be as potential biological treatment of PSC.     

Phylogeography of Phrynocephalus erythrurus from the Qiangtang Plateau of the Tibetan Plateau

Phrynocephalus erythrurus of the Qiangtang Plateau occupies the highest regions of any reptile on earth. Here, we report mitochondrial DNA haplotypes sampled throughout the distribution of P. erythrurus and analyze patterns of genetic divergence among populations. The species diverged into two major lineages/subspecies at 3.67 mya corresponding to the Northern and Southern Qiangtang Plateau. The Northern Qiangtang lineage diverged into two subpopulations at 2.76 mya separated by the Beilu River Region and Wulanwula Mountains. Haplotypes from the southern Qiangtang lineage diverged 0.98 mya as a star-shaped pattern. Analyses of molecular variance indicated that most of the observed genetic variation occurred among populations/regions implying long-term interruptions to gene flow. There was no evidence of sudden recent range expansions within any of the clades/lineages. NCPA infers allopatric fragmentation and restricted gene flow as the most likely mechanisms of population differentiation. Our results also indicate the presence of at least three refugia since the Hongya glaciation. Mountain movement and glaciations since mid-Pliocene are considered to have shaped phylogenetic patterns of P. erythrurus. P. erythrurus parva is suggested as a valid subspecies of P. erythrurus. Using four calibration points, we estimate an evolutionary rate of 0.762% divergence per lineage per million years for a mitochondrial genomic segment comprising the genes encoding ND2, tRNA^Trp and tRNA^Ala.     

Identification of a New P335 Subgroup through Molecular Analysis of Lactococcal Phages Q33 and BM13

Lactococcal dairy starter strains are under constant threat from phages in dairy fermentation facilities, especially by members of the so-called 936, P335, and c2 species. Among these three phage groups, members of the P335 species are the most genetically diverse. Here, we present the complete genome sequences of two P335-type phages, Q33 and BM13, isolated in North America and representing a novel lineage within this phage group. The Q33 and BM13 genomes exhibit homology, not only to P335-type, but also to elements of the 936-type phage sequences. The two phage genomes also have close relatedness to phages infecting Enterococcus and Clostridium, a heretofore unknown feature among lactococcal P335 phages. The Q33 and BM13 genomes are organized in functionally related clusters with genes encoding functions such as DNA replication and packaging, morphogenesis, and host cell lysis. Electron micrographic analysis of the two phages highlights the presence of a baseplate more reminiscent of the baseplate of 936 phages than that of the majority of members of the P335 group, with the exception of r1t and LC3.     

History cleans up messes: The impact of time in driving divergence and introgression in a tropical suture zone

Contact zones provide an excellent arena in which to address questions about how genomic divergence evolves during lineage divergence. They allow us to both infer patterns of genomic divergence in allopatric populations isolated from introgression and to characterize patterns of introgression after lineages meet. Thusly motivated, we analyze genome‐wide introgression data from four contact zones in three genera of lizards endemic to the Australian Wet Tropics. These contact zones all formed between morphologically cryptic lineage‐pairs within morphologically defined species, and the lineage‐pairs meeting in the contact zones diverged anywhere from 3.1 to 5.8 million years ago. By characterizing patterns of molecular divergence across an average of 11K genes and fitting geographic clines to an average of 7.5K variants, we characterize how patterns of genomic differentiation and introgression change through time. Across this range of divergences, we find that genome‐wide differentiation increases but becomes no less heterogeneous. In contrast, we find that introgression heterogeneity decreases dramatically, suggesting that time helps isolated genomes “congeal.” Thus, this work emphasizes the pivotal role that history plays in driving lineage divergence.     

Divergent evolution and molecular adaptation in the Drosophila odorant-binding protein family: inferences from sequence variation at the OS-E and OS-F genes

Background

Lineage-specific genes, the genes that are restricted to a limited subset of related organisms, may be important in adaptation. In parasitic organisms, lineage-specific gene products are possible targets for vaccine development or therapeutics when these genes are absent from the host genome.

Results

In this study, we utilized comparative approaches based on a phylogenetic framework to characterize lineage-specific genes in the parasitic protozoan phylum Apicomplexa. Genes from species in two major apicomplexan genera, Plasmodium and Theileria, were categorized into six levels of lineage specificity based on a nine-species phylogeny. In both genera, lineage-specific genes tend to have a higher level of sequence divergence among sister species. In addition, species-specific genes possess a strong codon usage bias compared to other genes in the genome. We found that a large number of genus- or species-specific genes are putative surface antigens that may be involved in host-parasite interactions. Interestingly, the two parasite lineages exhibit several notable differences. In Plasmodium, the (G + C) content at the third codon position increases with lineage specificity while Theileria shows the opposite trend. Surface antigens in Plasmodium are species-specific and mainly located in sub-telomeric regions. In contrast, surface antigens in Theileria are conserved at the genus level and distributed across the entire lengths of chromosomes.

Conclusion

Our results provide further support for the model that gene duplication followed by rapid divergence is a major mechanism for generating lineage-specific genes. The result that many lineage-specific genes are putative surface antigens supports the hypothesis that lineage-specific genes could be important in parasite adaptation. The contrasting properties between the lineage-specific genes in two major apicomplexan genera indicate that the mechanisms of generating lineage-specific genes and the subsequent evolutionary fates can differ between related parasite lineages. Future studies that focus on improving functional annotation of parasite genomes and collection of genetic variation data at within- and between-species levels will be important in facilitating our understanding of parasite adaptation and natural selection.     

Identification of Genes to Differentiate Closely Related Salmonella Lineages

Background

Salmonella are important human and animal pathogens. Though highly related, the Salmonella lineages may be strictly adapted to different hosts or cause different diseases, from mild local illness like gastroenteritis to fatal systemic infections like typhoid. Therefore, rapid and accurate identification of Salmonella is essential for timely and correct diagnosis of Salmonella infections. The current identification methods such as 16S rRNA sequencing and multilocus sequence typing are expensive and time consuming. Additionally, these methods often do not have sufficient distinguishing resolution among the Salmonella lineages.

Methodologies/Principal Findings

We compared 27 completely sequenced Salmonella genomes to identify possible genomic features that could be used for differentiation of individual lineages. We concatenated 2372 core genes in each of the 27 genomes and constructed a neighbor-joining tree. On the tree, strains of each serotype were clustered tightly together and different serotypes were unambiguously separated with clear genetic distances, demonstrating systematic genomic divergence among the Salmonella lineages. We made detailed comparisons among the 27 genomes and identified distinct sets of genomic differences, including nucleotide variations and genomic islands (GIs), among the Salmonella lineages. Two core genes STM4261 and entF together could unambiguously distinguish all Salmonella lineages compared in this study. Additionally, strains of a lineage have a common set of GIs and closely related lineages have similar sets of GIs.

Conclusions

Salmonella lineages have accumulated distinct sets of mutations and laterally acquired DNA (e.g., GIs) in evolution. Two genes entF and STM4261 have diverged sufficiently among the Salmonella lineages to be used for their differentiation. Further investigation of the distinct sets of mutations and GIs will lead to novel insights into genomic evolution of Salmonella and greatly facilitate the elucidation of pathogeneses of Salmonella infections.