Human milk microbiome: associations with maternal diet and infant growth

Introduction

Ingestion of human milk (HM) is identified as a significant factor associated with early infant gut microbial colonization, which has been associated with infant health and development. Maternal diet has been associated with the HM microbiome (HMM). However, a few studies have explored the associations among maternal diet, HMM, and infant growth during the first 6 months of lactation.

Methods

For this cross-sectional study, Mam-Mayan mother-infant dyads (n = 64) were recruited from 8 rural communities in the Western Highlands of Guatemala at two stages of lactation: early (6-46 days postpartum, n = 29) or late (109-184 days postpartum, n = 35). Recruited mothers had vaginally delivered singleton births, had no subclinical mastitis or antibiotic treatments, and breastfed their infants. Data collected at both stages of lactation included two 24-h recalls, milk samples, and infant growth status indicators: head-circumference-for-age-z-score (HCAZ), length-for-age-z-score (LAZ), and weight-for-age-z-score (WAZ). Infants were divided into subgroups: normal weight (WAZ ≥ -1SD) and mildly underweight (WAZ < -1SD), non-stunted (LAZ ≥ -1.5SD) and mildly stunted (LAZ < -1.5SD), and normal head-circumference (HCAZ ≥ -1SD) and smaller head-circumference (HCAZ < -1SD). HMM was identified using 16S rRNA gene sequencing; amplicon analysis was performed with the high-resolution ANCHOR pipeline, and DESeq2 identified the differentially abundant (DA) HMM at the species-level between infant growth groups (FDR < 0.05) in both early and late lactation.

Results

Using both cluster and univariate analyses, we identified (a) positive correlations between infant growth clusters and maternal dietary clusters, (b) both positive and negative associations among maternal macronutrient and micronutrient intakes with the HMM at the species level and (c) distinct correlations between HMM DA taxa with maternal nutrient intakes and infant z-scores that differed between breast-fed infants experiencing growth faltering and normal growth in early and late lactation.

Conclusion

Collectively, these findings provide important evidence of the potential influence of maternal diet on the early-life growth of breastfed infants via modulation of the HMM.

Cytobacts: Abundant and Diverse Vertically Seed-Transmitted Cultivation-Recalcitrant Intracellular Bacteria Ubiquitous to Vascular Plants

We have recently described ‘Cytobacts’ as abundant intracellular endophytic bacteria inhabiting live plant cells based on the observations with callus and cell suspension cultures of grapevine and other plant species with the origin ascribable to field explants. In this study, we investigated the prevalence of such cytoplasmic bacterial associations in field plants across different taxa, their cultivability, and the extent of taxonomic diversity and explored the possibility of their embryo-mediated vertical transmission. Over 100 genera of field plants were surveyed for ‘Cytobacts’ through bright-field live-cell imaging as per our previous experience using fresh tissue sections from surface-sterilized shoot-tissues with parallel cultivation-based assessments. This revealed widespread cellular bacterial associations visualized as copious motile micro-particles in the cytoplasm with no or sparse colony forming units (CFU) from the tissue-homogenates indicating their general non-cultivability. Based on the ease of detection and the abundance of ‘Cytobacts’ in fresh tissue sections, the surveyed plants were empirically classified into three groups: (i) motile bacteria detected instantly in most cells; (ii) motility not so widely observed, but seen in some cells; and (iii) only occasional motile units observed, but abundant non-motile bacterial cells present. Microscopy versus 16S-rRNA V3–V4 amplicon profiling on shoot-tip tissues of four representative plants—tomato, watermelon, periwinkle, and maize—showed high bacterial abundance and taxonomic diversity (11–15 phyla) with the dominance of Proteobacteria followed by Firmicutes/Actinobacteria, and several other phyla in minor shares. The low CFU/absence of bacterial CFU from the tissue homogenates on standard bacteriological media endorsed their cultivation-recalcitrance. Intracellular bacterial colonization implied that the associated organisms are able to transmit vertically to the next generation through the seed-embryos. Microscopy and 16S-rRNA V3–V4 amplicon/metagenome profiling of mature embryos excised from fresh watermelon seeds revealed heavy embryo colonization by diverse bacteria with sparse or no CFU. Observations with grapevine fresh fruit-derived seeds and seed-embryos endorsed the vertical transmission by diverse cultivation-recalcitrant endophytic bacteria (CREB). By and large, Proteobacteria formed the major phylum in fresh seed-embryos with varying shares of diverse phyla. Thus, we document ‘Cytobacts’ comprising diverse and vertically transmissible CREBs as a ubiquitous phenomenon in vascular plants.

Fast and Accurate Classification of Meta-Genomics Long Reads With deSAMBA

There is still a lack of fast and accurate classification tools to identify the taxonomies of noisy long reads, which is a bottleneck to the use of the promising long-read metagenomic sequencing technologies. Herein, we propose de Bruijn graph-based Sparse Approximate Match Block Analyzer (deSAMBA), a tailored long-read classification approach that uses a novel pseudo alignment algorithm based on sparse approximate match block (SAMB). Benchmarks on real sequencing datasets demonstrate that deSAMBA enables to achieve high yields and fast speed simultaneously, which outperforms state-of-the-art tools and has many potentials to cutting-edge metagenomics studies.

Distinct Changes Occur in the Human Breast Milk Microbiome Between Early and Established Lactation in Breastfeeding Guatemalan Mothers

Human breast milk contains a diverse community of bacteria, but as breast milk microbiome studies have largely focused on mothers from high income countries where few women breastfeed to 6 months, the temporal changes in the breast milk microbiome that occur during later lactation stages have not been explored. For this cross-sectional study, microbiota from breast milk samples of Mam-Mayan mothers living in eight remote rural communities in the Western Highlands of Guatemala were analyzed. All mothers delivered vaginally and breastfed their infants for 6 months. Breast milk from 76 unrelated mothers was used to compare two lactation stages, either “early” (6–46 days post-partum, n = 33) or “late” (109–184 days post-partum, n = 43). Breast milk microbial communities were assessed using 16S ribosomal RNA gene sequencing and lactation stages were compared using DESeq2 differential abundance analysis. A total of 1,505 OTUs were identified, including 287 which could be annotated as putative species. Among several maternal factors, lactation stage explained microbiome variance and inertia in ordination with the most significance (p < 0.001). Differential abundance analysis identified 137 OTUs as significantly higher in either early or late lactation. These included a general shift from Staphylococcus and Streptococcus species in early lactation to Sphingobium and Pseudomonas species in late lactation. Species enriched in early lactation included putative commensal bacteria known to colonize the infant oral and intestinal tracts whereas species enriched in late lactation had a uniform functional trait associated with aromatic compound degradation. Differentially abundant species also included several species which have not previously been reported within breast milk, such as Janthinobacterium agaricidamnosum, Novosphingobium clariflavum, Ottowia beijingensis, and Flavobacterium cucumis. These discoveries describe temporal changes to the breast milk microbiome of healthy Guatemalan mothers from early to late lactation. Collectively, these findings illustrate how studying under-represented human populations might advance our understanding of factors that modulate the human milk microbiome in low and middle income countries (LMIC).

Can Targeting Non-Contiguous V-Regions With Paired-End Sequencing Improve 16S rRNA-Based Taxonomic Resolution of Microbiomes?: An In Silico Evaluation

Background: Next-generation sequencing (NGS) technologies have enabled probing of microbial diversity in different environmental niches with unprecedented sequencing depth. However, due to read-length limitations of popular NGS technologies, 16S amplicon sequencing-based microbiome studies rely on targeting short stretches of the 16S rRNA gene encompassing a selection of variable (V) regions. In most cases, such a short stretch constitutes a single V-region or a couple of V-regions placed adjacent to each other on the 16S rRNA gene. Given that different V-regions have different resolving ability with respect to various taxonomic groups, selecting the optimal V-region (or a combination thereof) remains a challenge.

Methods: The accuracy of taxonomic profiles generated from sequences encompassing 1) individual V-regions, 2) adjacent V-regions, and 3) pairs of non-contiguous V-regions were assessed and compared. Subsequently, the discriminating capability of different V-regions with respect to different taxonomic lineages was assessed. The possibility of using paired-end sequencing protocols to target combinations of non-adjacent V-regions was finally evaluated with respect to the utility of such an experimental design in providing improved taxonomic resolution.

Results: Extensive validation with simulated microbiome datasets mimicking different environmental and host-associated microbiome samples suggest that targeting certain combinations of non-contiguously placed V-regions might yield better taxonomic classification accuracy compared to conventional 16S amplicon sequencing targets. This work also puts forward a novel in silico combinatorial strategy that enables creation of consensus taxonomic profiles from experiments targeting multiple pair-wise combinations of V-regions to improve accuracy in taxonomic classification.

Conclusion: The study suggests that targeting non-contiguous V-regions with paired-end sequencing can improve 16S rRNA–based taxonomic resolution of microbiomes. Furthermore, employing the novel in silico combinatorial strategy can improve taxonomic classification without any significant additional experimental costs and/or efforts. The empirical observations obtained can potentially serve as a guideline for future 16S microbiome studies, and facilitate researchers in choosing the optimal combination of V-regions for a specific experiment/sampled environment.