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.

Protein and Microbial Biomarkers in Sputum Discern Acute and Latent Tuberculosis in Investigation of Pastoral Ethiopian Cohort.

Differential diagnosis of tuberculosis (TB) and latent TB infection (LTBI) remains a public health priority in high TB burden countries. Pulmonary TB is diagnosed by sputum smear microscopy, chest X-rays, and PCR tests for distinct Mycobacterium tuberculosis (Mtb) genes. Clinical tests to diagnose LTBI rely on immune cell stimulation in blood plasma with TB-specific antigens followed by measurements of interferon-γ concentrations. The latter is an important cytokine for cellular immune responses against Mtb in infected lung tissues. Sputum smear microscopy and chest X-rays are not sufficiently sensitive while both PCR and interferon-γ release assays are expensive. Alternative biomarkers for the development of diagnostic tests to discern TB disease states are desirable. This study's objective was to discover sputum diagnostic biomarker candidates from the analysis of samples from 161 human subjects including TB patients, individuals with LTBI, negative community controls (NCC) from the province South Omo, a pastoral region in Ethiopia. We analyzed 16S rRNA gene-based bacterial taxonomies and proteomic profiles. The sputum microbiota did not reveal statistically significant differences in α-diversity comparing the cohorts. The genus Mycobacterium, representing Mtb, was only identified for the TB group which also featured reduced abundance of the genus Rothia in comparison with the LTBI and NCC groups. Rothia is a respiratory tract commensal and may be sensitive to the inflammatory milieu generated by infection with Mtb. Proteomic data supported innate immune responses against the pathogen in subjects with pulmonary TB. Ferritin, an iron storage protein released by damaged host cells, was markedly increased in abundance in TB sputum compared to the LTBI and NCC groups, along with the α-1-acid glycoproteins ORM1 and ORM2. These proteins are acute phase reactants and inhibit excessive neutrophil activation. Proteomic data highlight the effector roles of neutrophils in the anti-Mtb response which was not observed for LTBI cases. Less abundant in the sputum of the LTBI group, compared to the NCC group, were two immunomodulatory proteins, mitochondrial TSPO and the extracellular ribonuclease T2. If validated, these proteins are of interest as new biomarkers for diagnosis of LTBI.

Genetic diversity of clinical and environmental Mucorales isolates obtained from an investigation of mucormycosis cases among solid organ transplant recipients.

Mucormycoses are invasive infections by Rhizopus species and other Mucorales. Over 10 months, four solid organ transplant (SOT) recipients at our centre developed mucormycosis due to Rhizopus microsporus (n=2), R. arrhizus (n=1) or Lichtheimia corymbifera (n=1), at a median 31.5 days (range: 13-34) post-admission. We performed whole genome sequencing (WGS) on 72 Mucorales isolates (45 R. arrhizus, 19 R. delemar, six R. microsporus, two Lichtheimia species) from these patients, from five patients with community-acquired mucormycosis, and from hospital and regional environments. Isolates were compared by core protein phylogeny and global genomic features, including genome size, guanine-cytosine percentages, shared protein families and paralogue expansions. Patient isolates fell into six core phylogenetic lineages (clades). Phylogenetic and genomic similarities of R. microsporus isolates recovered 7 months apart from two SOT recipients in adjoining hospitals suggested a potential common source exposure. However, isolates from other patients and environmental sites had unique genomes. Many isolates that were indistinguishable by core phylogeny were distinct by one or more global genomic comparisons. Certain clades were recovered throughout the study period, whereas others were found at particular time points. In conclusion, mucormycosis cases could not be genetically linked to a definitive environmental source. Comprehensive genomic analyses eliminated false associations between Mucorales isolates that would have been assigned using core phylogenetic or less extensive genomic comparisons. The genomic diversity of Mucorales mandates that multiple isolates from individual patients and environmental sites undergo WGS during epidemiological investigations. However, exhaustive surveillance of fungal populations in a hospital and surrounding community is probably infeasible.

Characterization of Pneumococcal Colonization Dynamics and Antimicrobial Resistance Using Shotgun Metagenomic Sequencing in Intensively Sampled South African Infants.

Background: There remains a significant proportion of deaths due to pneumococcal pneumonia in infants from low- and middle-income countries despite the marginal global declines recorded in the past decade. Monitoring changes in pneumococcal carriage is key to understanding vaccination-induced shifts in the ecology of carriage, patterns of antimicrobial resistance, and impact on health. We longitudinally investigated pneumococcal carriage dynamics in PCV-13 vaccinated infants by collecting nasopharyngeal (NP) samples at 2-weekly intervals from birth through the first year of life from 137 infants. As a proof of concept, 196 NP samples were retrieved from a subset of 23 infants to explore strain-level pneumococcal colonization patterns and associated antimicrobial-resistance determinants. These were selected on the basis of changes in serotype and antibiogram over time. NP samples underwent short-term enrichment for streptococci prior to total nucleic acid extraction and whole metagenome shotgun sequencing (WMGS). Reads were assembled and aligned to pneumococcal reference genomes for the extraction of pneumococcal and non-pneumococcal bacterial reads. Pneumococcal contigs were aligned to the Antibiotic Resistance Gene-ANNOTation database of acquired AMR genes. In silico pneumococcal capsular and multilocus sequence typing were performed. Results: Of the 196 samples sequenced, 174 had corresponding positive cultures for pneumococci, of which, 152 were assigned an in silico serotype. Metagenomic sequencing detected a single pneumococcal serotype in 85% (129/152), and co-colonization in 15% (23/152) of the samples. Twenty-two different pneumococcal serotypes were identified, with 15B/15C and 16F being the most common non-PCV13 serotypes, while 23F and 19A were the most common PCV13 serotypes. Twenty-six different sequence types (STs), including four novel STs were identified in silico. Mutations in the folA and folP genes, associated with cotrimoxazole resistance, were detected in 89% (87/98) of cotrimoxazole-non-susceptible pneumococci, as well as in the pbp1a and pbp2x genes, in penicillin non-susceptible ST705215B/15C isolates. Conclusions: Metagenomic sequencing of NP samples is a valuable culture-independent technique for a detailed evaluation of the pneumococcal component and resistome of the NP microbiome. This method allowed for the detection of novel STs, as well as co-colonization, with a predominance of non-PCV13 serotypes in this cohort. Forty-eight resistance genes, as well as mutations associated with resistance were detected, but the correlation with phenotypic non-susceptibility was lower than expected.

Longitudinal changes in the nasopharyngeal resistome of South African infants using shotgun metagenomic sequencing.

INTRODUCTION: Nasopharyngeal (NP) colonization with antimicrobial-resistant bacteria is a global public health concern. Antimicrobial-resistance (AMR) genes carried by the resident NP microbiota may serve as a reservoir for transfer of resistance elements to opportunistic pathogens. Little is known about the NP antibiotic resistome. This study longitudinally investigated the composition of the NP antibiotic resistome in Streptococcus-enriched samples in a South African birth cohort. METHODS: As a proof of concept study, 196 longitudinal NP samples were retrieved from a subset of 23 infants enrolled as part of broader birth cohort study. These were selected on the basis of changes in serotype and antibiogram over time. NP samples underwent short-term enrichment for streptococci prior to total nucleic acid extraction and whole metagenome shotgun sequencing (WMGS). Reads were assembled and aligned to pneumococcal reference genomes for the extraction of streptococcal and non-streptococcal bacterial reads. Contigs were aligned to the Antibiotic Resistance Gene-ANNOTation database of acquired AMR genes. RESULTS: AMR genes were detected in 64% (125/196) of the samples. A total of 329 AMR genes were detected, including 36 non-redundant genes, ranging from 1 to 14 genes per sample. The predominant AMR genes detected encoded resistance mechanisms to beta-lactam (52%, 172/329), macrolide-lincosamide-streptogramin (17%, 56/329), and tetracycline antibiotics (12%, 38/329). MsrD, ermB, and mefA genes were only detected from streptococcal reads. The predominant genes detected from non- streptococcal reads included blaOXA-60, blaOXA-22, and blaBRO-1. Different patterns of carriage of AMR genes were observed, with only one infant having a stable carriage of mefA, msrD and tetM over a long period. CONCLUSION: This study demonstrates that WMGS can provide a broad snapshot of the NP resistome and has the potential to provide a comprehensive assessment of resistance elements present in this niche.

Annotated Genome Sequence of Aspergillus tanneri NIH1004.

The annotated genome of Aspergillus tanneri, a recently discovered drug-resistant pathogen, was determined by employing the Oxford Nanopore MinION platform and the Funannotate pipeline. The genome size and the number of protein-coding genes are notably larger than those of the most common etiological agent of aspergillosis, Aspergillus fumigatus.

HIV-exposure, early life feeding practices and delivery mode impacts on faecal bacterial profiles in a South African birth cohort.

There are limited data on meconium and faecal bacterial profiles from African infants and their mothers. We characterized faecal bacterial communities of infants and mothers participating in a South African birth cohort. Stool and meconium specimens were collected from 90 mothers and 107 infants at birth, and from a subset of 72 and 36 infants at 4-12 and 20-28 weeks of age, respectively. HIV-unexposed infants were primarily exclusively breastfed at 4-12 (49%, 26/53) and 20-28 weeks (62%, 16/26). In contrast, HIV-exposed infants were primarily exclusively formula fed at 4-12 (53%; 10/19) and 20-28 weeks (70%, 7/10). Analysis (of the bacterial 16S rRNA gene sequences of the V4 hypervariable region) of the 90 mother-infant pairs showed that meconium bacterial profiles [dominated by Proteobacteria (89%)] were distinct from those of maternal faeces [dominated by Firmicutes (66%) and Actinobacteria (15%)]. Actinobacteria predominated at 4-12 (65%) and 20-28 (50%) weeks. HIV-exposed infants had significantly higher faecal bacterial diversities at both 4-12 (p = 0.026) and 20-28 weeks (p = 0.002). HIV-exposed infants had lower proportions of Bifidobacterium (p = 0.010) at 4-12 weeks. Maternal faecal bacterial profiles were influenced by HIV status, feeding practices and mode of delivery. Further longitudinal studies are required to better understand how these variables influence infant and maternal faecal bacterial composition.

Strain Level Streptococcus Colonization Patterns during the First Year of Life.

Pneumococcal pneumonia has decreased significantly since the implementation of the pneumococcal conjugate vaccine (PCV), nevertheless, in many developing countries pneumonia mortality in infants remains high. We have undertaken a study of the nasopharyngeal (NP) microbiome during the first year of life in infants from The Philippines and South Africa. The study entailed the determination of the Streptococcus sp. carriage using a lytA qPCR assay, whole metagenomic sequencing, and in silico serotyping of Streptococcus pneumoniae, as well as 16S rRNA amplicon based community profiling. The lytA carriage in both populations increased with infant age and lytA+ samples ranged from 24 to 85% of the samples at each sampling time point. We next developed informatic tools for determining Streptococcus community composition and pneumococcal serotype from metagenomic sequences derived from a subset of longitudinal lytA-positive Streptococcus enrichment cultures from The Philippines (n = 26 infants, 50% vaccinated) and South African (n = 7 infants, 100% vaccinated). NP samples from infants were passaged in enrichment media, and metagenomic DNA was purified and sequenced. In silico capsular serotyping of these 51 metagenomic assemblies assigned known serotypes in 28 samples, and the co-occurrence of serotypes in 5 samples. Eighteen samples were not typeable using known serotypes but did encode for capsule biosynthetic cluster genes similar to non-encapsulated reference sequences. In addition, we performed metagenomic assembly and 16S rRNA amplicon profiling to understand co-colonization dynamics of Streptococcus sp. and other NP genera, revealing the presence of multiple Streptococcus species as well as potential respiratory pathogens in healthy infants. A range of virulence and drug resistant elements were identified as circulating in the NP microbiomes of these infants. This study revealed the frequent co-occurrence of multiple S. pneumoniae strains along with Streptococcus sp. and other potential pathogens such as S. aureus in the NP microbiome of these infants. In addition, the in silico serotype analysis proved powerful in determining the serotypes in S. pneumoniae carriage, and may lead to developing better targeted vaccines to prevent invasive pneumococcal disease (IPD) in these countries. These findings suggest that NP colonization by S. pneumoniae during the first years of life is a dynamic process involving multiple serotypes and species.

Genetic Analysis Using an Isogenic Mating Pair of Aspergillus fumigatus Identifies Azole Resistance Genes and Lack of MAT Locus's Role in Virulence.

Invasive aspergillosis (IA) due to Aspergillus fumigatus is a major cause of mortality in immunocompromised patients. The discovery of highly fertile strains of A. fumigatus opened the possibility to merge classical and contemporary genetics to address key questions about this pathogen. The merger involves sexual recombination, selection of desired traits, and genomics to identify any associated loci. We constructed a highly fertile isogenic pair of A. fumigatus strains with opposite mating types and used them to investigate whether mating type is associated with virulence and to find the genetic loci involved in azole resistance. The pair was made isogenic by 9 successive backcross cycles of the foundational strain AFB62 (MAT1-1) with a highly fertile (MAT1-2) progeny. Genome sequencing showed that the F9 MAT1-2 progeny was essentially identical to the AFB62. The survival curves of animals infected with either strain in three different animal models showed no significant difference, suggesting that virulence in A. fumigatus was not associated with mating type. We then employed a relatively inexpensive, yet highly powerful strategy to identify genomic loci associated with azole resistance. We used traditional in vitro drug selection accompanied by classical sexual crosses of azole-sensitive with resistant isogenic strains. The offspring were plated under varying drug concentrations and pools of resulting colonies were analyzed by whole genome sequencing. We found that variants in 5 genes contributed to azole resistance, including mutations in erg11A (cyp51A), as well as multi-drug transporters, erg25, and in HMG-CoA reductase. The results demonstrated that with minimal investment into the sequencing of three pools from a cross of interest, the variation(s) that contribute any phenotype can be identified with nucleotide resolution. This approach can be applied to multiple areas of interest in A. fumigatus or other heterothallic pathogens, especially for virulence associated traits.

Using Bayesian modelling to investigate factors governing antibiotic-induced Candida albicans colonization of the GI tract.

Receipt of broad-spectrum antibiotics enhances Candida albicans colonization of the GI tract, a risk factor for haematogenously-disseminated candidiasis. To understand how antibiotics influence C. albicans colonization, we treated mice orally with vancomycin or a combination of penicillin, streptomycin, and gentamicin (PSG) and then inoculated them with C. albicans by gavage. Only PSG treatment resulted in sustained, high-level GI colonization with C. albicans. Furthermore, PSG reduced bacterial diversity in the colon much more than vancomycin. Both antibiotic regimens significantly reduced IL-17A, IL-21, IL-22 and IFN-γ mRNA levels in the terminal ileum but had limited effect on the GI fungal microbiome. Through a series of models that employed Bayesian model averaging, we investigated the associations between antibiotic treatment, GI microbiota, and host immune response and their collective impact on C. albicans colonization. Our analysis revealed that bacterial genera were typically associated with either C. albicans colonization or altered cytokine expression but not with both. The only exception was Veillonella, which was associated with both increased C. albicans colonization and reduced IL-21 expression. Overall, antibiotic-induced changes in the bacterial microbiome were much more consistent determinants of C. albicans colonization than either the GI fungal microbiota or the GI immune response.

A systematic evaluation of high-dimensional, ensemble-based regression for exploring large model spaces in microbiome analyses.

BACKGROUND: Microbiome studies incorporate next-generation sequencing to obtain profiles of microbial communities. Data generated from these experiments are high-dimensional with a rich correlation structure but modest sample sizes. A statistical model that utilizes these microbiome profiles to explain a clinical or biological endpoint needs to tackle high-dimensionality resulting from the very large space of variable configurations. Ensemble models are a class of approaches that can address high-dimensionality by aggregating information across large model spaces. Although such models are popular in fields as diverse as economics and genetics, their performance on microbiome data has been largely unexplored. RESULTS: We developed a simulation framework that accurately captures the constraints of experimental microbiome data. Using this setup, we systematically evaluated a selection of both frequentist and Bayesian regression modeling ensembles. These are represented by variants of stability selection in conjunction with elastic net and spike-and-slab Bayesian model averaging (BMA), respectively. BMA ensembles that explore a larger space of models relative to stability selection variants performed better and had lower variability across simulations. However, stability selection ensembles were able to match the performance of BMA in scenarios of low sparsity where several variables had large regression coefficients. CONCLUSIONS: Given a microbiome dataset of interest, we present a methodology to generate simulated data that closely mimics its characteristics in a manner that enables meaningful evaluation of analytical strategies. Our evaluation demonstrates that the largest ensembles yield the strongest performance on microbiome data with modest sample sizes and high-dimensional measurements. We also demonstrate the ability of these ensembles to identify microbiome signatures that are associated with opportunistic Candida albicans colonization during antibiotic exposure. As the focus of microbiome research evolves from pilot to translational studies, we anticipate that our strategy will aid investigators in making evaluation-based decisions for selecting appropriate analytical methods.

Large-scale transcriptional response to hypoxia in Aspergillus fumigatus observed using RNAseq identifies a novel hypoxia regulated ncRNA.

We utilized RNAseq analysis of the Aspergillus fumigatus response to early hypoxic condition exposure. The results show that more than 89% of the A. fumigatus genome is expressed under normoxic and hypoxic conditions. Replicate samples were highly reproducible; however, comparisons between normoxia and hypoxia revealed that >23 and 35% of genes were differentially expressed after 30 and 120 min of hypoxia exposure, respectively. Consistent with our previous report detailing transcriptomic and proteomic responses at later time points, the results here show major repression of ribosomal function and induction of ergosterol biosynthesis, as well as activation of alternate respiratory mechanisms at the later time point. RNAseq data were used to define 32 hypoxia-specific genes, which were not expressed under normoxic conditions. Transcripts of a C6 transcription factor and a histidine kinase-response regulator were found only in hypoxia. In addition, several genes involved in the phosphoenylpyruvate and D-glyceraldehyde-3-phosphate metabolism were only expressed in hypoxia. Interestingly, a 216-bp ncRNA Afu-182 in the 3' region of insA (AFUB_064770) was significantly repressed under hypoxia with a 40-fold reduction in expression. A detailed analysis of Afu-182 showed similarity with several genes in the genome, many of which were also repressed in hypoxia. The results from this study show that hypoxia induces very early and widely drastic genome-wide responses in A. fumigatus that include expression of protein-coding and ncRNA genes. The role of these ncRNA genes in regulating the fungal hypoxia response is an exciting future research direction.