Hackflex library preparation enables low-cost metagenomic profiling.

Shotgun metagenomic sequencing provides valuable insights into microbial communities, but the high cost of library preparation with standard kits and protocols is a barrier for many. New methods such as Hackflex use diluted commercially available reagents to greatly reduce library preparation costs. However, these methods have not been systematically validated for metagenomic sequencing. Here, we evaluate Hackflex performance by sequencing metagenomic libraries from known mock communities as well as mouse fecal samples prepared by Hackflex, Illumina DNA Prep, and Illumina TruSeq methods. Hackflex successfully recovered all members of the Zymo mock community, performing best for samples with DNA concentrations <1 ng/uL. Furthermore, Hackflex was able to delineate microbiota of individual inbred mice from the same breeding stock at the same mouse facility, and statistical modeling indicated that mouse ID explained a greater fraction of the variance in metagenomic composition than did library preparation method. These results show that Hackflex is suitable for generating inventories of bacterial communities through metagenomic sequencing.

Female behavior drives the formation of distinct social structures in C57BL/6J versus wild-derived outbred mice in field enclosures.

BACKGROUND: Social behavior and social organization have major influences on individual health and fitness. Yet, biomedical research focuses on studying a few genotypes under impoverished social conditions. Understanding how lab conditions have modified social organizations of model organisms, such as lab mice, relative to natural populations is a missing link between socioecology and biomedical science. RESULTS: Using a common garden design, we describe the formation of social structure in the well-studied laboratory mouse strain, C57BL/6J, in replicated mixed-sex populations over 10-day trials compared to control trials with wild-derived outbred house mice in outdoor field enclosures. We focus on three key features of mouse social systems: (i) territory establishment in males, (ii) female social relationships, and (iii) the social networks formed by the populations. Male territorial behaviors were similar but muted in C57 compared to wild-derived mice. Female C57 sharply differed from wild-derived females, showing little social bias toward cage mates and exploring substantially more of the enclosures compared to all other groups. Female behavior consistently generated denser social networks in C57 than in wild-derived mice. CONCLUSIONS: C57 and wild-derived mice individually vary in their social and spatial behaviors which scale to shape overall social organization. The repeatable societies formed under field conditions highlights opportunities to experimentally study the interplay between society and individual biology using model organisms.

Assessing co-diversification in host-associated microbiomes.

When lineages of hosts and microbial symbionts engage in intimate interactions over evolutionary timescales, they can diversify in parallel (i.e., co-diversify), producing associations between the lineages' phylogenetic histories. Tests for co-diversification of individual microbial lineages and their hosts have been developed previously, and these have been applied to discover ancient symbioses in diverse branches of the tree of life. However, most host-microbe relationships are not binary but multipartite, in that a single host-associated microbiota can contain many microbial lineages, generating challenges for assessing co-diversification. Here, we review recent evidence for co-diversification in complex microbiota, highlight the limitations of prior studies, and outline a hypothesis testing approach designed to overcome some of these limitations. We advocate for the use of microbiota-wide scans for co-diversifying symbiont lineages and discuss tools developed for this purpose. Tests for co-diversification for simple host symbiont systems can be extended to entire phylogenies of microbial lineages (e.g., metagenome-assembled or isolate genomes, amplicon sequence variants) sampled from host clades, thereby providing a means for identifying co-diversifying symbionts present within complex microbiota. The relative ages of symbiont clades can corroborate co-diversification, and multi-level permutation tests can account for multiple comparisons and phylogenetic non-independence introduced by repeated sampling of host species. Discovering co-diversifying lineages will generate powerful opportunities for interrogating the molecular evolution and lineage turnover of ancestral, host-species specific symbionts within host-associated microbiota.

When microbes go missing: Understanding the impact of diversity loss within the gut microbiome.

Two recent papers published in Cell highlight the power of both top-down and bottom-up approaches to understanding the gut microbiome. The first uses ultra-deep sequencing to identify patterns across a gradient of human industrialization, while the second uses synthetic communities to determine how strain interactions impact microbiome structure and function.

Widespread extinctions of co-diversified primate gut bacterial symbionts from humans.

Humans and other primates harbour complex gut bacterial communities that influence health and disease, but the evolutionary histories of these symbioses remain unclear. This is partly due to limited information about the microbiota of ancestral primates. Here, using phylogenetic analyses of metagenome-assembled genomes (MAGs), we show that hundreds of gut bacterial clades diversified in parallel (that is, co-diversified) with primate species over millions of years, but that humans have experienced widespread losses of these ancestral symbionts. Analyses of 9,460 human and non-human primate MAGs, including newly generated MAGs from chimpanzees and bonobos, revealed significant co-diversification within ten gut bacterial phyla, including Firmicutes, Actinobacteriota and Bacteroidota. Strikingly, ~44% of the co-diversifying clades detected in African apes were absent from available metagenomic data from humans and ~54% were absent from industrialized human populations. In contrast, only ~3% of non-co-diversifying clades detected in African apes were absent from humans. Co-diversifying clades present in both humans and chimpanzees displayed consistent genomic signatures of natural selection between the two host species but differed in functional content from co-diversifying clades lost from humans, consistent with selection against certain functions. This study discovers host-species-specific bacterial symbionts that predate hominid diversification, many of which have undergone accelerated extinctions from human populations.

Home-site advantage for host species-specific gut microbiota.

Mammalian species harbor compositionally distinct gut microbial communities, but the mechanisms that maintain specificity of symbionts to host species remain unclear. Here, we show that natural selection within house mice (Mus musculus domesticus) drives deterministic assembly of the house-mouse gut microbiota from mixtures of native and non-native microbiotas. Competing microbiotas from wild-derived lines of house mice and other mouse species (Mus and Peromyscus spp.) within germ-free wild-type (WT) and Rag1-knockout (Rag1-/-) house mice revealed widespread fitness advantages for native gut bacteria. Native bacterial lineages significantly outcompeted non-native lineages in both WT and Rag1-/- mice, indicating home-site advantage for native microbiota independent of host adaptive immunity. However, a minority of native Bacteriodetes and Firmicutes favored by selection in WT hosts were not favored or disfavored in Rag1-/- hosts, indicating that Rag1 mediates fitness advantages of these strains. This study demonstrates home-site advantage for native gut bacteria, consistent with local adaptation of gut microbiota to their mammalian species.

Influential factors of saliva microbiota composition.

The oral microbiota is emerging as an influential factor of host physiology and disease state. Factors influencing oral microbiota composition have not been well characterised. In particular, there is a lack of population-based studies. We undertook a large hypothesis-free study of the saliva microbiota, considering potential influential factors of host health (frailty; diet; periodontal disease), demographics (age; sex; BMI) and sample processing (storage time), in a sample (n = 679) of the TwinsUK cohort of adult twins. Alpha and beta diversity of the saliva microbiota was associated most strongly with frailty (alpha diversity: ß = -0.16, Q = 0.003, Observed; ß = -0.16, Q = 0.002, Shannon; ß = -0.16, Q = 0.003, Simpson; Beta diversity: Q = 0.002, Bray Curtis dissimilarity) and age (alpha diversity: ß = 0.15, Q = 0.006, Shannon; ß = 0.12, Q = 0.003, Simpson; beta diversity: Q = 0.002, Bray Curtis dissimilarity; Q = 0.032, Weighted UniFrac) in multivariate models including age, frailty, sex, BMI, frailty and diet, and adjustment for multiple testing. Those with a more advanced age were more likely to be dissimilar in the saliva microbiota composition than younger participants (P = 5.125e-06, ANOVA). In subsample analyses, including consideration of periodontal disease (total n = 138, periodontal disease n = 66), the association with frailty remained for alpha diversity (Q = 0.002, Observed ASVs; Q = 0.04 Shannon Index), but not beta diversity, whilst age was not demonstrated to associate with alpha or beta diversity in this subsample, potentially due to insufficient statistical power. Length of time that samples were stored prior to sequencing was associated with beta diversity (Q = 0.002, Bray Curtis dissimilarity). Six bacterial taxa were associated with age after adjustment for frailty and diet. Of the factors studied, frailty and age emerged as the most influential with regards to saliva microbiota composition. Whilst age and frailty are correlates, the associations were independent of each other, giving precedence to both biological and chronological ageing as processes of potential importance when considering saliva microbiota composition.

Microbiomes: Infant Chimps Crawling with Bacteria.

Human guts are colonized at birth by a limited set of microbes that gradually increases in diversity throughout infancy. A new study reports the opposite pattern in infant chimpanzees, raising questions about how host-microbiota relationships have changed during ape evolution.

Microbiota assembly, structure, and dynamics among Tsimane horticulturalists of the Bolivian Amazon.

Selective and neutral forces shape human microbiota assembly in early life. The Tsimane are an indigenous Bolivian population with infant care-associated behaviors predicted to increase mother-infant microbial dispersal. Here, we characterize microbial community assembly in 47 infant-mother pairs from six Tsimane villages, using 16S rRNA gene amplicon sequencing of longitudinal stool and tongue swab samples. We find that infant consumption of dairy products, vegetables, and chicha (a fermented drink inoculated with oral microbes) is associated with stool microbiota composition. In stool and tongue samples, microbes shared between mothers and infants are more abundant than non-shared microbes. Using a neutral model of community assembly, we find that neutral processes alone explain the prevalence of 79% of infant-colonizing microbes, but explain microbial prevalence less well in adults from river villages with more regular access to markets. Our results underscore the importance of neutral forces during microbiota assembly. Changing lifestyle factors may alter traditional modes of microbiota assembly by decreasing the role of neutral processes.

Comparison of pectin-degrading fungal communities in temperate forests using glycosyl hydrolase family 28 pectinase primers targeting Ascomycete fungi.

Fungi have developed a wide assortment of enzymes to break down pectin, a prevalent polymer in plant cell walls that is important in plant defense and structure. One enzyme family used to degrade pectin is the glycosyl hydrolase family 28 (GH28). In this study we developed primers for the amplification of GH28 coding genes from a database of 293 GH28 sequences from 40 fungal genomes. The primers were used to successfully amplify GH28 pectinases from all Ascomycota cultures tested, but only three out of seven Basidiomycota cultures. In addition, we further tested the primers in PCRs on metagenomic DNA extracted from senesced tree leaves from different forest ecosystems, followed by cloning and sequencing. Taxonomic specificity for Ascomycota GH28 genes was tested by comparing GH28 composition in leaves to internal transcribed spacer (ITS) amplicon composition using pyrosequencing. All sequences obtained from GH28 primers were classified as Ascomycota; in contrast, ITS sequences indicated that fungal communities were up to 39% Basidiomycetes. Analysis of leaf samples indicated that both forest stand and ecosystem type were important in structuring fungal communities. However, site played the prominent role in explaining GH28 composition, whereas ecosystem type was more important for ITS composition, indicating possible genetic drift between populations of fungi. Overall, these primers will have utility in understanding relationships between fungal community composition and ecosystem processes, as well as detection of potentially pathogenic Ascomycetes.

Use of 16S rRNA sequencing and quantitative PCR to correlate venous leg ulcer bacterial bioburden dynamics with wound expansion, antibiotic therapy, and healing.

Clinical diagnosis of infection in chronic wounds is currently limited to subjective clinical signs and culture-based methods that underestimate the complexity of wound microbial bioburden as revealed by DNA-based microbial identification methods. Here, we use 16S rRNA next generation sequencing and quantitative polymerase chain reaction to characterize weekly changes in bacterial load, community structure, and diversity associated with a chronic venous leg ulcer over the 15-week course of treatment and healing. Our DNA-based methods and detailed sampling scheme reveal that the bacterial bioburden of the wound is unexpectedly dynamic, including changes in the bacterial load and community structure that correlate with wound expansion, antibiotic therapy, and healing. We demonstrate that these multidimensional changes in bacterial bioburden can be summarized using swabs taken prior to debridement, and therefore, can be more easily collected serially than debridement or biopsy samples. Overall, this case illustrates the importance of detailed clinical indicators and longitudinal sampling to determine the pathogenic significance of chronic wound microbial dynamics and guide best use of antimicrobials for improvement of healing outcomes.

Evolutionary analysis of glycosyl hydrolase family 28 (GH28) suggests lineage-specific expansions in necrotrophic fungal pathogens.

Glycosyl hydrolase family 28 (GH28) is a set of structurally related enzymes that hydrolyze glycosidic bonds in pectin, and are important extracellular enzymes for both pathogenic and saprotrophic fungi. Yet, very little is understood about the evolutionary forces driving the diversification of GH28s in fungal genomes. We reconstructed the evolutionary history of family GH28 in fungi by examining the distribution of GH28 copy number across the phylogeny of fungi, and by reconstructing the phylogeny of GH28 genes. We also examined the relationship between lineage-specific GH28 expansions and fungal ecological strategy, testing the hypothesis that GH28 evolution in fungi is driven by ecological strategy (pathogenic vs. non-pathogenic) and pathogenic niche (necrotrophic vs. biotrophic). Our results showed that GH28 phylogeny of Ascomycota and Basidiomycota sequences was structured by specific biochemical function, with endo-polygalacturonases and endo-rhamnogalacturonases forming distinct, apparently ancient clades, while exo-polygalacturonases are more widely distributed. In contrast, Mucoromycotina and Stramenopile sequences formed taxonomically-distinct clades. Large, lineage-specific variation in GH28 copy number indicates that the evolution of this gene family is consistent with the birth-and-death model of gene family evolution, where diversity of GH28 loci within genomes was generated through multiple rounds of gene duplication followed by functional diversification and loss of some gene family members. Although GH28 copy number was correlated with genome size, our findings suggest that ecological strategy also plays an important role in determining the GH28 repertoire of fungi. Both necrotrophic and biotrophic fungi have larger genomes than non-pathogens, yet only necrotrophs possess more GH28 enzymes than non-pathogens. Hence, lineage-specific GH28 expansion is the result of both variation in genome size across fungal species and diversifying selection within the necrotrophic plant pathogen ecological niche. GH28 evolution among necrotrophs has likely been driven by a co-evolutionary arms race with plants, whereas the need to avoid plant immune responses has resulted in purifying selection within biotrophic fungi.