Succession and determinants of the early life nasopharyngeal microbiota in a South African birth cohort.

BACKGROUND: Bacteria colonizing the nasopharynx play a key role as gatekeepers of respiratory health. Yet, dynamics of early life nasopharyngeal (NP) bacterial profiles remain understudied in low- and middle-income countries (LMICs), where children have a high prevalence of risk factors for lower respiratory tract infection. We investigated longitudinal changes in NP bacterial profiles, and associated exposures, among healthy infants from low-income households in South Africa. METHODS: We used short fragment (V4 region) 16S rRNA gene amplicon sequencing to characterize NP bacterial profiles from 103 infants in a South African birth cohort, at monthly intervals from birth through the first 12 months of life and six monthly thereafter until 30 months. RESULTS: Corynebacterium and Staphylococcus were dominant colonizers at 1 month of life; however, these were rapidly replaced by Moraxella- or Haemophilus-dominated profiles by 4 months. This succession was almost universal and largely independent of a broad range of exposures. Warm weather (summer), lower gestational age, maternal smoking, no day-care attendance, antibiotic exposure, or low height-for-age z score at 12 months were associated with higher alpha and beta diversity. Summer was also associated with higher relative abundances of Staphylococcus, Streptococcus, Neisseria, or anaerobic gram-negative bacteria, whilst spring and winter were associated with higher relative abundances of Haemophilus or Corynebacterium, respectively. Maternal smoking was associated with higher relative abundances of Porphyromonas. Antibiotic therapy (or isoniazid prophylaxis for tuberculosis) was associated with higher relative abundance of anerobic taxa (Porphyromonas, Fusobacterium, and Prevotella) and with lower relative abundances of health associated-taxa Corynebacterium and Dolosigranulum. HIV-exposure was associated with higher relative abundances of Klebsiella or Veillonella and lower relative abundances of an unclassified genus within the family Lachnospiraceae. CONCLUSIONS: In this intensively sampled cohort, there was rapid and predictable replacement of early profiles dominated by health-associated Corynebacterium and Dolosigranulum with those dominated by Moraxella and Haemophilus, independent of exposures. Season and antibiotic exposure were key determinants of NP bacterial profiles. Understudied but highly prevalent exposures prevalent in LMICs, including maternal smoking and HIV-exposure, were associated with NP bacterial profiles. Video Abstract.

Sputum bacterial load and bacterial composition correlate with lung function and are altered by long-term azithromycin treatment in children with HIV-associated chronic lung disease.

BACKGROUND: Long-term azithromycin (AZM) treatment reduces the frequency of acute respiratory exacerbation in children and adolescents with HIV-associated chronic lung disease (HCLD). However, the impact of this treatment on the respiratory bacteriome is unknown. METHOD: African children with HCLD (defined as forced expiratory volume in 1 s z-score (FEV1z) less than - 1.0 with no reversibility) were enrolled in a placebo-controlled trial of once-weekly AZM given for 48-weeks (BREATHE trial). Sputum samples were collected at baseline, 48 weeks (end of treatment) and 72 weeks (6 months post-intervention in participants who reached this timepoint before trial conclusion). Sputum bacterial load and bacteriome profiles were determined using 16S rRNA gene qPCR and V4 region amplicon sequencing, respectively. The primary outcomes were within-participant and within-arm (AZM vs placebo) changes in the sputum bacteriome measured across baseline, 48 weeks and 72 weeks. Associations between clinical or socio-demographic factors and bacteriome profiles were also assessed using linear regression. RESULTS: In total, 347 participants (median age: 15.3 years, interquartile range [12.7-17.7]) were enrolled and randomised to AZM (173) or placebo (174). After 48 weeks, participants in the AZM arm had reduced sputum bacterial load vs placebo arm (16S rRNA copies/µl in log10, mean difference and 95% confidence interval [CI] of AZM vs placebo - 0.54 [- 0.71; - 0.36]). Shannon alpha diversity remained stable in the AZM arm but declined in the placebo arm between baseline and 48 weeks (3.03 vs. 2.80, p = 0.04, Wilcoxon paired test). Bacterial community structure changed in the AZM arm at 48 weeks compared with baseline (PERMANOVA test p = 0.003) but resolved at 72 weeks. The relative abundances of genera previously associated with HCLD decreased in the AZM arm at 48 weeks compared with baseline, including Haemophilus (17.9% vs. 25.8%, p < 0.05, ANCOM ω = 32) and Moraxella (1% vs. 1.9%, p < 0.05, ANCOM ω = 47). This reduction was sustained at 72 weeks relative to baseline. Lung function (FEV1z) was negatively associated with bacterial load (coefficient, [CI]: - 0.09 [- 0.16; - 0.02]) and positively associated with Shannon diversity (0.19 [0.12; 0.27]). The relative abundance of Neisseria (coefficient, [standard error]: (2.85, [0.7], q = 0.01), and Haemophilus (- 6.1, [1.2], q < 0.001) were positively and negatively associated with FEV1z, respectively. An increase in the relative abundance of Streptococcus from baseline to 48 weeks was associated with improvement in FEV1z (3.2 [1.11], q = 0.01) whilst an increase in Moraxella was associated with decline in FEV1z (-2.74 [0.74], q = 0.002). CONCLUSIONS: AZM treatment preserved sputum bacterial diversity and reduced the relative abundances of the HCLD-associated genera Haemophilus and Moraxella. These bacteriological effects were associated with improvement in lung function and may account for reduced respiratory exacerbations associated with AZM treatment of children with HCLD. Video Abstract.

The Influence of DNA Extraction and Lipid Removal on Human Milk Bacterial Profiles.

Culture-independent molecular techniques have advanced the characterization of environmental and human samples including the human milk (HM) bacteriome. However, extraction of high-quality genomic DNA that is representative of the bacterial population in samples is crucial. Lipids removal from HM prior to DNA extraction is common practice, but this may influence the bacterial population detected. The objective of this study was to compare four commercial DNA extraction kits and lipid removal in relation to HM bacterial profiles. Four commercial DNA extraction kits, QIAamp® DNA Microbiome Kit, ZR Fungal/Bacterial DNA MiniPrep™, QIAsymphony DSP DNA Kit and ZymoBIOMICS™ DNA Miniprep Kit, were assessed using milk collected from ten healthy lactating women. The kits were evaluated based on their ability to extract high quantities of pure DNA from HM and how well they extracted DNA from bacterial communities present in a commercial mock microbial community standard spiked into HM. Finally, the kits were evaluated by assessing their extraction repeatability. Bacterial profiles were assessed using Illumina MiSeq sequencing targeting the V4 region of the 16S rRNA gene. The ZR Fungal/Bacterial DNA MiniPrep™ and ZymoBIOMICS™ DNA Miniprep (Zymo Research Corp., Irvine, CA, USA) kits extracted the highest DNA yields with the best purity. DNA extracted using ZR Fungal/Bacterial DNA MiniPrep™ best represented the bacteria in the mock community spiked into HM. In un-spiked HM samples, DNA extracted using the QIAsymphony DSP DNA kit showed statistically significant differences in taxa prevalence from DNA extracted using ZR Fungal/Bacterial DNA MiniPrep™ and ZymoBIOMICS™ DNA Miniprep kits. The only difference between skim and whole milk is observed in bacterial profiles with differing relative abundances of Enhydrobacter and Acinetobacter. DNA extraction, but not lipids removal, substantially influences bacterial profiles detected in HM samples, emphasizing the need for careful selection of a DNA extraction kit to improve DNA recovery from a range of bacterial taxa.

Optimizing 16S rRNA gene profile analysis from low biomass nasopharyngeal and induced sputum specimens.

BACKGROUND: Careful consideration of experimental artefacts is required in order to successfully apply high-throughput 16S ribosomal ribonucleic acid (rRNA) gene sequencing technology. Here we introduce experimental design, quality control and "denoising" approaches for sequencing low biomass specimens. RESULTS: We found that bacterial biomass is a key driver of 16S rRNA gene sequencing profiles generated from bacterial mock communities and that the use of different deoxyribonucleic acid (DNA) extraction methods [DSP Virus/Pathogen Mini Kit® (Kit-QS) and ZymoBIOMICS DNA Miniprep Kit (Kit-ZB)] and storage buffers [PrimeStore® Molecular Transport medium (Primestore) and Skim-milk, Tryptone, Glucose and Glycerol (STGG)] further influence these profiles. Kit-QS better represented hard-to-lyse bacteria from bacterial mock communities compared to Kit-ZB. Primestore storage buffer yielded lower levels of background operational taxonomic units (OTUs) from low biomass bacterial mock community controls compared to STGG. In addition to bacterial mock community controls, we used technical repeats (nasopharyngeal and induced sputum processed in duplicate, triplicate or quadruplicate) to further evaluate the effect of specimen biomass and participant age at specimen collection on resultant sequencing profiles. We observed a positive correlation (r = 0.16) between specimen biomass and participant age at specimen collection: low biomass technical repeats (represented by < 500 16S rRNA gene copies/µl) were primarily collected at < 14 days of age. We found that low biomass technical repeats also produced higher alpha diversities (r = - 0.28); 16S rRNA gene profiles similar to no template controls (Primestore); and reduced sequencing reproducibility. Finally, we show that the use of statistical tools for in silico contaminant identification, as implemented through the decontam package in R, provides better representations of indigenous bacteria following decontamination. CONCLUSIONS: We provide insight into experimental design, quality control steps and "denoising" approaches for 16S rRNA gene high-throughput sequencing of low biomass specimens. We highlight the need for careful assessment of DNA extraction methods and storage buffers; sequence quality and reproducibility; and in silico identification of contaminant profiles in order to avoid spurious results.

Composition of Gut Microbiota of Children and Adolescents With Perinatal Human Immunodeficiency Virus Infection Taking Antiretroviral Therapy in Zimbabwe.

BACKGROUND: Human immunodeficiency virus (HIV) infection causes impairment of the gastrointestinal barrier, with substantial depletion of CD4+ T cells in the gut. Antiretroviral therapy (ART) restores CD4+ counts and may have beneficial effects on gut microbiota in adults. Little is known about effect of long-term ART on gut microbiome in HIV-infected children. We investigated composition of gut microbiota in HIV-infected and -uninfected children and assessed associations between gut microbiota and patient characteristics. METHODS: In a cross-sectional study, rectal swabs were collected from 177 HIV-infected and 103 HIV-uninfected controls. Gut microbial composition was explored using 16S ribosomal ribonucleic acid sequencing. RESULTS: Human immunodeficiency virus-infected children had significantly lower alpha-diversity and higher beta-diversity compared to HIV-uninfected. No association was observed between microbiome diversity and CD4+ T-cell count, HIV viral load, or HIV-associated chronic lung disease. We found enriched levels of Corynebacterium (P < .01), Finegoldia (P < .01), and Anaerococcus (P < .01) in HIV-infected participants and enrichment of Enterobacteriaceae (P = .02) in participants with low CD4+ counts (<400 cells/mm3). Prolonged ART-treatment (≥10 years) was significantly associated with a richer gut microbiota by alpha diversity. CONCLUSIONS: Human immunodeficiency virus-infected children have altered gut microbiota. Prolonged ART may restore the richness of the microbiota closer to that of HIV-uninfected children.

Influence of Socio-Economic and Psychosocial Profiles on the Human Breast Milk Bacteriome of South African Women.

The human breast milk (HBM) bacteriome is an important, continuous source of microbes to the neonate in early life, playing an important role in shaping the infant's intestinal bacteriome. Study of the composition of the HBM bacteriome is an emerging area of research, with little information available, particularly from low- and middle-income countries. The aim of this study was to characterize the diversity of bacterial communities in HBM samples collected between 6-10 weeks postpartum from lactating South African women and to study potential influencing factors of the bacteriome. Using 16S rRNA gene sequencing of samples from 554 women, we demonstrated that the HBM bacteriome was largely dominated by the phyla Firmicutes (mean relative abundance: 71.1%) and Actinobacteria (mean relative abundance: 16.4%). The most abundant genera identified from the HBM bacteriome were Streptococcus (mean relative abundance: 48.6%), Staphylococcus (mean relative abundance: 17.8%), Rothia (mean relative abundance: 5.8%), and Corynebacterium (mean relative abundance: 4.3%). "Core" bacterial genera including Corynebacterium, Streptococcus, Staphylococcus, Rothia, Veillonella, Gemella, Acinetobacter, Micrococcus and a genus belonging to the Enterobacteriaceae family were present in 80% of samples. HBM samples were classified, according to their bacteriome, into three major clusters, dominated by the genera Staphylococcus (cluster 1), a combination of Staphylococcus and Streptococcus (cluster 2), and Streptococcus (cluster 3). The cluster groups differed significantly for Shannon and chao1 richness indices. Bacterial interactions were studied using co-occurrence networks with positive associations observed between the abundances of Staphylococcus and Corynebacteria (members of the skin microflora) and between Streptococcus, Rothia, Veillonella, and Gemella (members of the oral microflora). HBM from older mothers had a higher Shannon diversity index. The study site was associated with differences in HBM bacteriome composition (permutational multivariate analysis of variance using distance matrices (PERMANOVA), p < 0.05). No other tested socio-demographic or psychosocial factors were associated with HBM bacterial composition.

Host and Microbiome Genome-Wide Association Studies: Current State and Challenges.

The involvement of the microbiome in health and disease is well established. Microbiome genome-wide association studies (mGWAS) are used to elucidate the interaction of host genetic variation with the microbiome. The emergence of this relatively new field has been facilitated by the advent of next generation sequencing technologies that enable the investigation of the complex interaction between host genetics and microbial communities. In this paper, we review recent studies investigating host-microbiome interactions using mGWAS. Additionally, we highlight the marked disparity in the sampling population of mGWAS carried out to date and draw attention to the critical need for inclusion of diverse populations.