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.

Predicting antimicrobial mechanism-of-action from transcriptomes: A generalizable explainable artificial intelligence approach

To better combat the expansion of antibiotic resistance in pathogens, new compounds, particularly those with novel mechanisms-of-action [MOA], represent a major research priority in biomedical science. However, rediscovery of known antibiotics demonstrates a need for approaches that accurately identify potential novelty with higher throughput and reduced labor. Here we describe an explainable artificial intelligence classification methodology that emphasizes prediction performance and human interpretability by using a Hierarchical Ensemble of Classifiers model optimized with a novel feature selection algorithm called Clairvoyance; collectively referred to as a CoHEC model. We evaluated our methods using whole transcriptome responses from Escherichia coli challenged with 41 known antibiotics and 9 crude extracts while depositing 122 transcriptomes unique to this study. Our CoHEC model can properly predict the primary MOA of previously unobserved compounds in both purified forms and crude extracts at an accuracy above 99%, while also correctly identifying darobactin, a newly discovered antibiotic, as having a novel MOA. In addition, we deploy our methods on a recent E. coli transcriptomics dataset from a different strain and a Mycobacterium smegmatis metabolomics timeseries dataset showcasing exceptionally high performance; improving upon the performance metrics of the original publications. We not only provide insight into the biological interpretation of our model but also that the concept of MOA is a non-discrete heuristic with diverse effects for different compounds within the same MOA, suggesting substantial antibiotic diversity awaiting discovery within existing MOA.

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 ST7052^15B/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 bla_OXA-60, bla_OXA-22, and bla_BRO-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.

Genomic and phenotypic characterization of Burkholderia isolates from the potable water system of the International Space Station

The opportunistic pathogens Burkholderia cepacia and Burkholderia contaminans, both genomovars of the Burkholderia cepacia complex (BCC), are frequently cultured from the potable water dispenser (PWD) of the International Space Station (ISS). Here, we sequenced the genomes and conducted phenotypic assays to characterize these Burkholderia isolates. All recovered isolates of the two species fall within monophyletic clades based on phylogenomic trees of conserved single-copy core genes. Within species, the ISS-derived isolates all demonstrate greater than 99% average nucleotide identity (with 95-99% of genomes aligning) and share around 90% of the identified gene clusters from a pangenomic analysis-suggesting that the two groups are each composed of highly similar genomic lineages and their members may have all stemmed from the same two founding populations. The differences that can be observed between the recovered isolates at the pangenomic level are primarily located within putative plasmids. Phenotypically, macrophage intracellularization and lysis occurred at generally similar rates between all ISS-derived isolates, as well as with their respective type-terrestrial strain references. All ISS-derived isolates exhibited antibiotic sensitivity similar to that of the terrestrial reference strains, and minimal differences between isolates were observed. With a few exceptions, biofilm formation rates were generally consistent across each species. And lastly, though isolation date does not necessarily provide any insight into how long a given isolate had been aboard the ISS, none of the assayed physiology correlated with either date of isolation or distances based on nucleotide variation. Overall, we find that while the populations of Burkholderia present in the ISS PWS each maintain virulence, they are likely are not more virulent than those that might be encountered on planet and remain susceptible to clinically used antibiotics.

Whole genome comparison of Aspergillus flavus L-morphotype strain NRRL 3357 (type) and S-morphotype strain AF70

Aspergillus flavus is a saprophytic fungus that infects corn, peanuts, tree nuts and other agriculturally important crops. Once the crop is infected the fungus has the potential to secrete one or more mycotoxins, the most carcinogenic of which is aflatoxin. Aflatoxin contaminated crops are deemed unfit for human or animal consumption, which results in both food and economic losses. Within A. flavus, two morphotypes exist: the S strains (small sclerotia) and L strains (large sclerotia). Significant morphological and physiological differences exist between the two morphotypes. For example, the S-morphotypes produces sclerotia that are smaller (< 400 μm), greater in quantity, and contain higher concentrations of aflatoxin than the L-morphotypes (>400 μm). The morphotypes also differ in pigmentation, pH homeostasis in culture and the number of spores produced. Here we report the first full genome sequence of an A. flavus S morphotype, strain AF70. We provide a comprehensive comparison of the A. flavus S-morphotype genome sequence with a previously sequenced genome of an L-morphotype strain (NRRL 3357), including an in-depth analysis of secondary metabolic clusters and the identification SNPs within their aflatoxin gene clusters.

Secondary metabolite arsenal of an opportunistic pathogenic fungus

Aspergillus fumigatus is a versatile fungus able to successfully exploit diverse environments from mammalian lungs to agricultural waste products. Among its many fitness attributes are dozens of genetic loci containing biosynthetic gene clusters (BGCs) producing bioactive small molecules (often referred to as secondary metabolites or natural products) that provide growth advantages to the fungus dependent on environment. Here we summarize the current knowledge of these BGCs-18 of which can be named to product-their expression profiles in vivo, and which BGCs may enhance virulence of this opportunistic human pathogen. Furthermore, we find extensive evidence for the presence of many of these BGCs, or similar BGCs, in distantly related genera including the emerging pathogen Pseudogymnoascus destructans, the causative agent of white-nose syndrome in bats, and suggest such BGCs may be predictive of pathogenic potential in other fungi.This article is part of the themed issue 'Tackling emerging fungal threats to animal health, food security and ecosystem resilience'.

Mutual independence of alkaline- and calcium-mediated signalling in Aspergillus fumigatus refutes the existence of a conserved druggable signalling nexus

Functional coupling of calcium‐ and alkaline responsive signalling occurs in multiple fungi to afford efficient cation homeostasis. Host microenvironments exert alkaline stress and potentially toxic concentrations of Ca²⁺, such that highly conserved regulators of both calcium‐ (Crz) and pH‐ (PacC/Rim101) responsive signalling are crucial for fungal pathogenicity. Drugs targeting calcineurin are potent antifungal agents but also perturb human immunity thereby negating their use as anti‐infectives, abrogation of alkaline signalling has, therefore, been postulated as an adjunctive antifungal strategy. We examined the interdependency of pH‐ and calcium‐mediated signalling in Aspergillus fumigatus and found that calcium chelation severely impedes hyphal growth indicating a critical requirement for this ion independently of ambient pH. Transcriptomic responses to alkaline pH or calcium excess exhibited minimal similarity. Mutants lacking calcineurin, or its client CrzA, displayed normal alkaline tolerance and nuclear translocation of CrzA was unaffected by ambient pH. Expression of a highly conserved, alkaline‐regulated, sodium ATPase was tolerant of genetic or chemical perturbations of calcium‐mediated signalling, but abolished in null mutants of the pH‐responsive transcription factor PacC, and PacC proteolytic processing occurred normally during calcium excess. Taken together our data demonstrate that in A. fumigatus the regulatory hierarchy governing alkaline tolerance circumvents calcineurin signalling.

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.

Reducing the Bottleneck in Discovery of Novel Antibiotics

Most antibiotics were discovered by screening soil actinomycetes, but the efficiency of the discovery platform collapsed in the 1960s. By now, more than 3000 antibiotics have been described and most of the current discovery effort is focused on the rediscovery of known compounds, making the approach impractical. The last marketed broad-spectrum antibiotics discovered were daptomycin, linezolid, and fidaxomicin. The current state of the art in the development of new anti-infectives is a non-existent pipeline in the absence of a discovery platform. This is particularly troubling given the emergence of pan-resistant pathogens. The current practice in dealing with the problem of the background of known compounds is to use chemical dereplication of extracts to assess the relative novelty of a compound it contains. Dereplication typically requires scale-up, extraction, and often fractionation before an accurate mass and structure can be produced by MS analysis in combination with 2D NMR. Here, we describe a transcriptome analysis approach using RNA sequencing (RNASeq) to identify promising novel antimicrobial compounds from microbial extracts. Our pipeline permits identification of antimicrobial compounds that produce distinct transcription profiles using unfractionated cell extracts. This efficient pipeline will eliminate the requirement for purification and structure determination of compounds from extracts and will facilitate high-throughput screen of cell extracts for identification of novel compounds.

Ecology of aspergillosis: insights into the pathogenic potency of Aspergillus fumigatus and some other Aspergillus species

Fungi of the genus Aspergillus are widespread in the environment. Some Aspergillus species, most commonly Aspergillus fumigatus, may lead to a variety of allergic reactions and life-threatening systemic infections in humans. Invasive aspergillosis occurs primarily in patients with severe immunodeficiency, and has dramatically increased in recent years. There are several factors at play that contribute to aspergillosis, including both fungus and host-related factors such as strain virulence and host pulmonary structure/immune status, respectively. The environmental tenacity of Aspergilllus, its dominance in diverse microbial communities/habitats, and its ability to navigate the ecophysiological and biophysical challenges of host infection are attributable, in large part, to a robust stress-tolerance biology and exceptional capacity to generate cell-available energy. Aspects of its stress metabolism, ecology, interactions with diverse animal hosts, clinical presentations and treatment regimens have been well-studied over the past years. Here, we synthesize these findings in relation to the way in which some Aspergillus species have become successful opportunistic pathogens of human- and other animal hosts. We focus on the biophysical capabilities of Aspergillus pathogens, key aspects of their ecophysiology and the flexibility to undergo a sexual cycle or form cryptic species. Additionally, recent advances in diagnosis of the disease are discussed as well as implications in relation to questions that have yet to be resolved.

Candida albicans Transcriptional Profiling Within Biliary Fluid From a Patient With Cholangitis, Before and After Antifungal Treatment and Surgical Drainage

We used ribonucleic acid sequencing to profile Candida albicans transcription within biliary fluid from a patient with cholangitis; samples were collected before and after treatment with fluconazole and drainage. Candida albicans transcriptomes at the infection site distinguished treated from untreated cholangitis. After treatment, 1131 C. albicans genes were differentially expressed in biliary fluid. Up-regulated genes were enriched in hyphal growth, cell wall organization, adhesion, oxidation reduction, biofilm, and fatty acid and ergosterol biosynthesis. This is the first study to define Candida global gene expression during deep-seated human infection. Successful treatment of cholangitis induced C. albicans genes involved in fluconazole responses and pathogenesis.

Looking Beyond Respiratory Cultures: Microbiome-Cytokine Signatures of Bacterial Pneumonia and Tracheobronchitis in Lung Transplant Recipients

Bacterial pneumonia and tracheobronchitis are diagnosed frequently following lung transplantation. The diseases share clinical signs of inflammation and are often difficult to differentiate based on culture results. Microbiome and host immune-response signatures that distinguish between pneumonia and tracheobronchitis are undefined. Using a retrospective study design, we selected 49 bronchoalveolar lavage fluid samples from 16 lung transplant recipients associated with pneumonia (n = 8), tracheobronchitis (n = 12) or colonization without respiratory infection (n = 29). We ensured an even distribution of Pseudomonas aeruginosa or Staphylococcus aureus culture-positive samples across the groups. Bayesian regression analysis identified non-culture-based signatures comprising 16S ribosomal RNA microbiome profiles, cytokine levels and clinical variables that characterized the three diagnoses. Relative to samples associated with colonization, those from pneumonia had significantly lower microbial diversity, decreased levels of several bacterial genera and prominent multifunctional cytokine responses. In contrast, tracheobronchitis was characterized by high microbial diversity and multifunctional cytokine responses that differed from those of pneumonia-colonization comparisons. The dissimilar microbiomes and cytokine responses underlying bacterial pneumonia and tracheobronchitis following lung transplantation suggest that the diseases result from different pathogenic processes. Microbiomes and cytokine responses had complementary features, suggesting that they are closely interconnected in the pathogenesis of both diseases.

The In vitro Antibiotic Tolerant Persister Population in Burkholderia pseudomallei is Altered by Environmental Factors

Bacterial persistence due to antibiotic tolerance is a critical aspect of antibiotic treatment failure, disease latency, and chronic or reemergent infections. The levels of persisters is especially notable for the opportunistic Gram-negative pathogens from the Burkholderia and Pseudomonas genera. We examined the rate of drug tolerant persisters in Burkholderia pseudomallei, Burkholderia thailandensis, Burkholderia cepacia complex organisms, and Pseudomonas aeruginosa at mid-log growth in LB broth culture. We found that a fraction of the antibiotic-sensitive cells from every species were tolerant to a 24 h high-dose antibiotic challenge. All tested Burkholderia strains demonstrated a drug tolerant persister population at a rate that was at least 100-500 times higher than P. aeruginosa. When challenged with at least a 10X minimum inhibitory concentration (MIC) 24 h exposure to three different antibiotics with different modes of action we found that in B. pseudomallei Bp82 each of the tree antibiotics revealed different persister fractions at each of two different growth states. This observation suggests that our assay is detecting heterogeneous persister subpopulations. Persistence in B. pseudomallei Bp82 was highly dependent on growth stage, with a surprisingly high persister fraction of >64% of the late stationary phase cells being antibiotic tolerant to 100XMIC cefotaxime. Adaptation of B. pseudomallei to distilled water storage resulted in a population of drug tolerant cells up to 100% of the non-drug-challenged viable cell count in the same cefotaxime assay. Cultivation of B. pseudomallei with a sub-inhibitory concentration of several antibiotics resulted in altered persister fractions within the population relative to cultures lacking the antibiotic. Our study provides insight into the sensitivity of the persister fraction within the population of B. pseudomallei due to environmental variables and suggests diversity within the persister population revealed by different challenge antibiotics.

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.

Perturbations in small molecule synthesis uncovers an iron-responsive secondary metabolite network in Aspergillus fumigatus

Iron plays a critical role in survival and virulence of the opportunistic pathogen Aspergillus fumigatus. Two transcription factors, the GATA-factor SreA and the bZip-factor HapX oppositely monitor iron homeostasis with HapX activating iron acquisition pathways (e.g., siderophores) and shutting down iron consumptive pathways (and SreA) during iron starvation conditions whereas SreA negatively regulates HapX and corresponding pathways during iron sufficiency. Recently the non-ribosomal peptide, hexadehydroastechrome (HAS; a tryptophan-derived iron (III)-complex), has been found important in A. fumigatus virulence. We found that HAS overproduction caused an iron starvation phenotype, from alteration of siderophore pools to regulation of iron homeostasis gene expression including sreA. Moreover, we uncovered an iron dependent secondary metabolism network where both SreA and HapX oppositely regulate multiple other secondary metabolites including HAS. This circuitry links iron-acquisition and consumption pathways with secondary metabolism-thus placing HAS as part of a metabolic feedback circuitry designed to balance iron pools in the fungus and presenting iron availability as one environmental trigger of secondary metabolism.

The tandem repeats enabling reversible switching between the two phases of β-lactamase substrate spectrum

Expansion or shrinkage of existing tandem repeats (TRs) associated with various biological processes has been actively studied in both prokaryotic and eukaryotic genomes, while their origin and biological implications remain mostly unknown. Here we describe various duplications (de novo TRs) that occurred in the coding region of a β-lactamase gene, where a conserved structure called the omega loop is encoded. These duplications that occurred under selection using ceftazidime conferred substrate spectrum extension to include the antibiotic. Under selective pressure with one of the original substrates (amoxicillin), a high level of reversion occurred in the mutant β-lactamase genes completing a cycle back to the original substrate spectrum. The de novo TRs coupled with reversion makes a genetic toggling mechanism enabling reversible switching between the two phases of the substrate spectrum of β-lactamases. This toggle exemplifies the effective adaptation of de novo TRs for enhanced bacterial survival. We found pairs of direct repeats that mediated the DNA duplication (TR formation). In addition, we found different duos of sequences that mediated the DNA duplication. These novel elements—that we named SCSs (same-strand complementary sequences)—were also found associated with β-lactamase TR mutations from clinical isolates. Both direct repeats and SCSs had a high correlation with TRs in diverse bacterial genomes throughout the major phylogenetic lineages, suggesting that they comprise a fundamental mechanism shaping the bacterial evolution.

β-lactamases can adapt to new antibiotics by mutations in their genes. The original and the extended substrate spectrums of β-lactamases define two phases of catalytic activity, and the conversion by point mutations is unidirectional from the initial to the new spectrum. We describe duplication mutations that enable reversible switching between the substrate spectrums, increasing the adaptability of the bacterium. We provide evidence supporting that two distinct groups of short sequences mediated the formation of DNA duplications in β-lactamases: direct repeats and novel elements that we named, SCSs (same-strand complementary sequences). Our study suggests that DNA duplication processes mediated by both direct repeats and SCSs are not just limited to the β-lactamase genes but comprise a fundamental mechanism in bacterial genome evolution.

Draft Genome Sequence of Penicillium expansum Strain R19, Which Causes Postharvest Decay of Apple Fruit

Among the species that cause blue mold, isolates of Penicillium expansum are the most prevalent and virulent species, causing more than 50 percent of postharvest decay. We report the draft genome sequence of P. expansum R19 in order to identify fungal virulence factors and to understand the mechanism of infection.

Standardized metadata for human pathogen/vector genomic sequences

High throughput sequencing has accelerated the determination of genome sequences for thousands of human infectious disease pathogens and dozens of their vectors. The scale and scope of these data are enabling genotype-phenotype association studies to identify genetic determinants of pathogen virulence and drug/insecticide resistance, and phylogenetic studies to track the origin and spread of disease outbreaks. To maximize the utility of genomic sequences for these purposes, it is essential that metadata about the pathogen/vector isolate characteristics be collected and made available in organized, clear, and consistent formats. Here we report the development of the GSCID/BRC Project and Sample Application Standard, developed by representatives of the Genome Sequencing Centers for Infectious Diseases (GSCIDs), the Bioinformatics Resource Centers (BRCs) for Infectious Diseases, and the U.S. National Institute of Allergy and Infectious Diseases (NIAID), part of the National Institutes of Health (NIH), informed by interactions with numerous collaborating scientists. It includes mapping to terms from other data standards initiatives, including the Genomic Standards Consortium's minimal information (MIxS) and NCBI's BioSample/BioProjects checklists and the Ontology for Biomedical Investigations (OBI). The standard includes data fields about characteristics of the organism or environmental source of the specimen, spatial-temporal information about the specimen isolation event, phenotypic characteristics of the pathogen/vector isolated, and project leadership and support. By modeling metadata fields into an ontology-based semantic framework and reusing existing ontologies and minimum information checklists, the application standard can be extended to support additional project-specific data fields and integrated with other data represented with comparable standards. The use of this metadata standard by all ongoing and future GSCID sequencing projects will provide a consistent representation of these data in the BRC resources and other repositories that leverage these data, allowing investigators to identify relevant genomic sequences and perform comparative genomics analyses that are both statistically meaningful and biologically relevant.

Polysome profiling reveals broad translatome remodeling during endoplasmic reticulum (ER) stress in the pathogenic fungus Aspergillus fumigatus

BACKGROUND

The unfolded protein response (UPR) is a network of intracellular signaling pathways that supports the ability of the secretory pathway to maintain a balance between the load of proteins entering the endoplasmic reticulum (ER) and the protein folding capacity of the ER lumen. Current evidence indicates that several pathogenic fungi rely heavily on this pathway for virulence, but there is limited understanding of the mechanisms involved. The best known functional output of the UPR is transcriptional upregulation of mRNAs involved in ER homeostasis. However, this does not take into account mechanisms of translational regulation that involve differential loading of ribosomes onto mRNAs. In this study, a global analysis of transcript-specific translational regulation was performed in the pathogenic mold Aspergillus fumigatus to determine the nature and scope of the translational response to ER stress.

RESULTS

ER stress was induced by treating the fungus with dithiothreitol, tunicamycin, or a thermal up-shift. The mRNAs were then fractionated on the basis of ribosome occupancy into an under-translated pool (U) and a well-translated pool (W). The mRNAs were used to interrogate microarrays and the ratio of the hybridization signal (W/U) was used as an indicator of the relative translational efficiency of a mRNA under each condition. The largest category of translationally upregulated mRNAs during ER stress encoded proteins involved in translation. Components of the ergosterol and GPI anchor biosynthetic pathways also showed increased polysome association, suggesting an important role for translational regulation in membrane and cell wall homeostasis. ER stress induced limited remodeling of the secretory pathway translatome. However, a select group of transcription factors was translationally upregulated, providing a link to subsequent modification of the transcriptome. Finally, we provide evidence that one component of the ER stress translatome is a novel mRNA isoform from the yvc1 gene that is induced by ER stress in a UPR-dependent manner.

CONCLUSIONS

Together, these findings define a core set of mRNAs subject to translational control during the adaptive response to acute ER stress in A. fumigatus and reveal a remarkable breadth of functions that are needed to resolve ER stress in this organism.

Mobile elements and mitochondrial genome expansion in the soil fungus and potato pathogen Rhizoctonia solani AG-3

The soil fungus Rhizoctonia solani is an economically important pathogen of agricultural and forestry crops. Here, we present the complete sequence and analysis of the mitochondrial genome of R. solani, field isolate Rhs1AP. The genome (235 849 bp) is the largest mitochondrial genome of a filamentous fungus sequenced to date and exhibits a rich accumulation of introns, novel repeat sequences, homing endonuclease genes, and hypothetical genes. Stable secondary structures exhibited by repeat sequences suggest that they comprise functional, possibly catalytic RNA elements. RNA-Seq expression profiling confirmed that the majority of homing endonuclease genes and hypothetical genes are transcriptionally active. Comparative analysis suggests that the mitochondrial genome of R. solani is an example of a dynamic history of expansion in filamentous fungi.

Aspergillus galactosaminogalactan mediates adherence to host constituents and conceals hyphal β-glucan from the immune system

Aspergillus fumigatus is the most common cause of invasive mold disease in humans. The mechanisms underlying the adherence of this mold to host cells and macromolecules have remained elusive. Using mutants with different adhesive properties and comparative transcriptomics, we discovered that the gene uge3, encoding a fungal epimerase, is required for adherence through mediating the synthesis of galactosaminogalactan. Galactosaminogalactan functions as the dominant adhesin of A. fumigatus and mediates adherence to plastic, fibronectin, and epithelial cells. In addition, galactosaminogalactan suppresses host inflammatory responses in vitro and in vivo, in part through masking cell wall β-glucans from recognition by dectin-1. Finally, galactosaminogalactan is essential for full virulence in two murine models of invasive aspergillosis. Collectively these data establish a role for galactosaminogalactan as a pivotal bifunctional virulence factor in the pathogenesis of invasive aspergillosis.

Invasive aspergillosis is the most common mold infection in humans, predominately affecting immunocompromised patients. The mechanisms by which the mold Aspergillus fumigatus adheres to host tissues and causes disease are poorly understood. In this report, we compared mutants of Aspergillus with different adhesive properties to identify fungal factors involved in adherence to host cells. This approach identified a cell wall associated polysaccharide, galactosaminogalactan, which is required for adherence to a wide variety of substrates. Galactosaminogalactan was also observed to suppress inflammation by concealing β-glucans, key pattern associated microbial pattern molecules in Aspergillus hyphae, from recognition by the innate immune system. Mutants that were deficient in galactosaminogalactan were less virulent in mouse models of invasive aspergillosis. These data identify a bifunctional role for galactosaminogalactan in the pathogenesis of invasive aspergillosis, and suggest that it may serve as a useful target for antifungal therapy.

Sequencing of mitochondrial genomes of nine Aspergillus and Penicillium species identifies mobile introns and accessory genes as main sources of genome size variability

BACKGROUND

The genera Aspergillus and Penicillium include some of the most beneficial as well as the most harmful fungal species such as the penicillin-producer Penicillium chrysogenum and the human pathogen Aspergillus fumigatus, respectively. Their mitochondrial genomic sequences may hold vital clues into the mechanisms of their evolution, population genetics, and biology, yet only a handful of these genomes have been fully sequenced and annotated.

RESULTS

Here we report the complete sequence and annotation of the mitochondrial genomes of six Aspergillus and three Penicillium species: A. fumigatus, A. clavatus, A. oryzae, A. flavus, Neosartorya fischeri (A. fischerianus), A. terreus, P. chrysogenum, P. marneffei, and Talaromyces stipitatus (P. stipitatum). The accompanying comparative analysis of these and related publicly available mitochondrial genomes reveals wide variation in size (25-36 Kb) among these closely related fungi. The sources of genome expansion include group I introns and accessory genes encoding putative homing endonucleases, DNA and RNA polymerases (presumed to be of plasmid origin) and hypothetical proteins. The two smallest sequenced genomes (A. terreus and P. chrysogenum) do not contain introns in protein-coding genes, whereas the largest genome (T. stipitatus), contains a total of eleven introns. All of the sequenced genomes have a group I intron in the large ribosomal subunit RNA gene, suggesting that this intron is fixed in these species. Subsequent analysis of several A. fumigatus strains showed low intraspecies variation. This study also includes a phylogenetic analysis based on 14 concatenated core mitochondrial proteins. The phylogenetic tree has a different topology from published multilocus trees, highlighting the challenges still facing the Aspergillus systematics.

CONCLUSIONS

The study expands the genomic resources available to fungal biologists by providing mitochondrial genomes with consistent annotations for future genetic, evolutionary and population studies. Despite the conservation of the core genes, the mitochondrial genomes of Aspergillus and Penicillium species examined here exhibit significant amount of interspecies variation. Most of this variation can be attributed to accessory genes and mobile introns, presumably acquired by horizontal gene transfer of mitochondrial plasmids and intron homing.

Clonality despite sex: the evolution of host-associated sexual neighborhoods in the pathogenic fungus Penicillium marneffei

Molecular genetic approaches typically detect recombination in microbes regardless of assumed asexuality. However, genetic data have shown the AIDS-associated pathogen Penicillium marneffei to have extensive spatial genetic structure at local and regional scales, and although there has been some genetic evidence that a sexual cycle is possible, this haploid fungus is thought to be genetically, as well as morphologically, asexual in nature because of its highly clonal population structure. Here we use comparative genomics, experimental mixed-genotype infections, and population genetic data to elucidate the role of recombination in natural populations of P. marneffei. Genome wide comparisons reveal that all the genes required for meiosis are present in P. marneffei, mating type genes are arranged in a similar manner to that found in other heterothallic fungi, and there is evidence of a putatively meiosis-specific mutational process. Experiments suggest that recombination between isolates of compatible mating types may occur during mammal infection. Population genetic data from 34 isolates from bamboo rats in India, Thailand and Vietnam, and 273 isolates from humans in China, India, Thailand, and Vietnam show that recombination is most likely to occur across spatially and genetically limited distances in natural populations resulting in highly clonal population structure yet sexually reproducing populations. Predicted distributions of three different spatial genetic clusters within P. marneffei overlap with three different bamboo rat host distributions suggesting that recombination within hosts may act to maintain population barriers within P. marneffei.

Twelve positions in a β-lactamase that can expand its substrate spectrum with a single amino acid substitution

The continuous evolution of β-lactamases resulting in bacterial resistance to β-lactam antibiotics is a major concern in public health, and yet the underlying molecular basis or the pattern of such evolution is largely unknown. We investigated the mechanics of the substrate fspectrum expansion of the class A β-lactamase using PenA of Burkholderia thailandensis as a model. By analyzing 516 mutated enzymes that acquired the ceftazidime-hydrolyzing activity, we found twelve positions with single amino acid substitutions (altogether twenty-nine different substitutions), co-localized at the active-site pocket area. The ceftazidime MIC (minimum inhibitory concentration) levels and the relative frequency in the occurrence of substitutions did not correlate well with each other, and the latter appeared be largely influenced by the intrinsic mutational biases present in bacteria. Simulation studies suggested that all substitutions caused a congruent effect, expanding the space in a conserved structure called the omega loop, which in turn increased flexibility at the active site. A second phase of selection, in which the mutants were placed under increased antibiotic pressure, did not result in a second mutation in the coding region, but a mutation that increased gene expression arose in the promoter. This result suggests that the twelve amino acid positions and their specific substitutions in PenA may represent a comprehensive repertoire of the enzyme's adaptability to a new substrate. These mapped substitutions represent a comprehensive set of general mechanical paths to substrate spectrum expansion in class A β-lactamases that all share a functional evolutionary mechanism using common conserved residues.

Development stage-specific proteomic profiling uncovers small, lineage specific proteins most abundant in the Aspergillus Fumigatus conidial proteome

Background

The pathogenic mold Aspergillus fumigatus is the most frequent infectious cause of death in severely immunocompromised individuals such as leukemia and bone marrow transplant patients. Germination of inhaled conidia (asexual spores) in the host is critical for the initiation of infection, but little is known about the underlying mechanisms of this process.

Results

To gain insights into early germination events and facilitate the identification of potential stage-specific biomarkers and vaccine candidates, we have used quantitative shotgun proteomics to elucidate patterns of protein abundance changes during early fungal development. Four different stages were examined: dormant conidia, isotropically expanding conidia, hyphae in which germ tube emergence has just begun, and pre-septation hyphae. To enrich for glycan-linked cell wall proteins we used an alkaline cell extraction method. Shotgun proteomic resulted in the identification of 375 unique gene products with high confidence, with no evidence for enrichment of cell wall-immobilized and secreted proteins. The most interesting discovery was the identification of 52 proteins enriched in dormant conidia including 28 proteins that have never been detected in the A. fumigatus conidial proteome such as signaling protein Pil1, chaperones BipA and calnexin, and transcription factor HapB. Additionally we found many small, Aspergillus specific proteins of unknown function including 17 hypothetical proteins. Thus, the most abundant protein, Grg1 (AFUA_5G14210), was also one of the smallest proteins detected in this study (M.W. 7,367). Among previously characterized proteins were melanin pigment and pseurotin A biosynthesis enzymes, histones H3 and H4.1, and other proteins involved in conidiation and response to oxidative or hypoxic stress. In contrast, expanding conidia, hyphae with early germ tubes, and pre-septation hyphae samples were enriched for proteins responsible for housekeeping functions, particularly translation, respiratory metabolism, amino acid and carbohydrate biosynthesis, and the tricarboxylic acid cycle.

Conclusions

The observed temporal expression patterns suggest that the A. fumigatus conidia are dominated by small, lineage-specific proteins. Some of them may play key roles in host-pathogen interactions, signal transduction during conidial germination, or survival in hostile environments.

Carbohydrate-active enzymes revealed in Coptotermes formosanus (Isoptera: Rhinotermitidae) transcriptome

Coptotermes formosanus is one of the most destructive wood-feeding termites. To understand the molecular mechanisms that regulate the development of the termite, a normalized C. formosanus cDNA library was constructed using mixed RNA isolated from workers, soldiers, nymphs and alates of both sexes. The sequencing of this library generated 131 636 expressed sequence tags (ESTs) and 25 939 assembled unigenes. The carbohydrate-active enzymes (CAZymes) revealed in this library were analysed in the present report. A total of 509 putative CAZymes were identified. Diverse cellulolytic enzymes were uncovered from both the host termite and from symbionts harboured by the termite, which were possibly the result of the high efficiency of cellulose utilization. CAZymes associated with trehalose biosynthetic and metabolic pathways were also identified, which are potential regulators of the physiological activities of trehalose, an important insect blood sugar. Representative CAZyme coding genes in glycoside hydrolase family 1 (GH1) were quantitatively analysed. The results showed that the five GH1 β-glucosidase genes were expressed differentially among different castes and one of them was female alate-specific. Overall, the normalized EST library provides a comprehensive genetic resource of C. formosanus and will serve a diverse range of research areas. The CAZymes represent one of the repositories of enzymes useful for physiological studies and applications in sugar-based biofuel production.

Cryptococcus gattii comparative genomics and transcriptomics: a NIH/NIAID White Paper

Cryptococcus gattii is an emerging global pathogen. Recent reports suggest that C. gattii cryptococcosis is more common in immunocompetent as well as HIV-infected AIDS patients than earlier estimated. An ongoing outbreak of C. gattii in Vancouver, Canada, and the US Pacific Northwest has heightened public health awareness in North America. We have few clues as to what causes emergence or re-emergence of highly pathogenic strains, why C. gattii split up from its sibling pathogen C. neoformans, why it thrives in trees instead, and why immunocompetent individuals are vulnerable to this pathogen? C. gattii comprises of four distinct lineages, but the information on the genome of C. gattii is inadequate and unrepresentative as it is limited to two strains, R265 and WM276, which are MATα, serotype B, genotype VGII/VGI from Canada and Australia, respectively. There is a wide gap in knowledge about the genomes of VGIII and VGIV strains, serotype C strains, and MATa strains. The geographical representation is inadequate in the absence of strains from California, South America, Asia, and Africa. Additional obstacles to work with this pathogen are the following: (a) complex molecular typing schemes and (b) lack of functional genomics analyses. We propose to complete genome sequencing of 12 reference strains by next-generation sequencing technology and to map their transcriptomes by RNA-Seq technology. This effort would lead to new resources for the scientific community including (1) insight from additional C. gattii genomes to anchor future research studies, (2) validation of single-nucleotide polymorphisms (SNPs) for molecular typing to improve epidemiology studies, and (3) transcript analyses from strains under relevant pathogenic and non-pathogenic conditions to accelerate the discovery of proteins for diagnostics, drug targets, and vaccines.

HacA-independent functions of the ER stress sensor IreA synergize with the canonical UPR to influence virulence traits in Aspergillus fumigatus

Endoplasmic reticulum (ER) stress is a condition in which the protein folding capacity of the ER becomes overwhelmed by an increased demand for secretion or by exposure to compounds that disrupt ER homeostasis. In yeast and other fungi, the accumulation of unfolded proteins is detected by the ER-transmembrane sensor IreA/Ire1, which responds by cleaving an intron from the downstream cytoplasmic mRNA HacA/Hac1, allowing for the translation of a transcription factor that coordinates a series of adaptive responses that are collectively known as the unfolded protein response (UPR). Here, we examined the contribution of IreA to growth and virulence in the human fungal pathogen Aspergillus fumigatus. Gene expression profiling revealed that A. fumigatus IreA signals predominantly through the canonical IreA-HacA pathway under conditions of severe ER stress. However, in the absence of ER stress IreA controls dual signaling circuits that are both HacA-dependent and HacA-independent. We found that a ΔireA mutant was avirulent in a mouse model of invasive aspergillosis, which contrasts the partial virulence of a ΔhacA mutant, suggesting that IreA contributes to pathogenesis independently of HacA. In support of this conclusion, we found that the ΔireA mutant had more severe defects in the expression of multiple virulence-related traits relative to ΔhacA, including reduced thermotolerance, decreased nutritional versatility, impaired growth under hypoxia, altered cell wall and membrane composition, and increased susceptibility to azole antifungals. In addition, full or partial virulence could be restored to the ΔireA mutant by complementation with either the induced form of the hacA mRNA, hacAⁱ, or an ireA deletion mutant that was incapable of processing the hacA mRNA, ireA^Δ10. Together, these findings demonstrate that IreA has both HacA-dependent and HacA-independent functions that contribute to the expression of traits that are essential for virulence in A. fumigatus.

Aspergillus fumigatus is the predominant mold pathogen of humans, responsible for life-threatening infections in patients with depressed immunity. The fungus is highly adapted for secretion, a feature that it uses to extract nutrients from the host environment. High rates of protein secretion can overwhelm the protein folding capacity of the endoplasmic reticulum (ER). The resulting ER stress is alleviated by the unfolded protein response (UPR), a signaling pathway that is triggered by the ER-membrane sensor IreA and executed by the downstream transcription factor HacA. This paper uncovers a novel role for IreA in the expression of multiple adaptive traits that allow the fungus to cope with stress conditions that are encountered during infection. Gene expression profiling of ΔireA and ΔhacA mutants revealed that IreA signals predominantly through the canonical IreA-HacA UPR pathway under extreme conditions of ER stress, but has unexpected HacA-dependent and HacA-independent functions even in the absence of ER stress. These findings establish IreA as an important regulator of A. fumigatus pathogenicity and suggest that therapeutic targeting of the dual functions of this protein could be an effective antifungal strategy.

Tight control of mycotoxin biosynthesis gene expression in Aspergillus flavus by temperature as revealed by RNA-Seq

To better understand the effect of temperature on mycotoxin biosynthesis, RNA-Seq technology was used to profile the Aspergillus flavus transcriptome under different temperature conditions. This approach allowed us to quantify transcript abundance for over 80% of fungal genes including 1153 genes that were differentially expressed at 30 and 37 °C. Eleven of the 55 secondary metabolite clusters were upregulated at the lower temperature, including aflatoxin biosynthesis genes, which were among the most highly upexpressed genes. On average, transcript abundance for the 30 aflatoxin biosynthesis genes was 3300 times greater at 30 °C as compared with 37 °C. The results are consistent with the view that high temperature negatively affects aflatoxin production by turning down transcription of the two key transcriptional regulators, aflR and aflS. Subtle changes in the expression levels of aflS to aflR appear to control transcription activation of the aflatoxin cluster.

Transcript amplification from single bacterium for transcriptome analysis

Total transcript amplification (TTA) from single eukaryotic cells for transcriptome analysis is established, but TTA from a single prokaryotic cell presents additional challenges with much less starting material, the lack of poly(A)-tails, and the fact that the messages can be polycistronic. Here, we describe a novel method for single-bacterium TTA using a model organism, Burkholderia thailandensis, exposed to a subinhibitory concentration of the antibacterial agent, glyphosate. Utilizing a B. thailandensis microarray to assess the TTA method showed low fold-change bias (less than twofold difference and Pearson correlation coefficient R ≈ 0.87-0.89) and drop-outs (4%-6% of 2842 detectable genes), compared with data obtained from the larger-scale nonamplified RNA samples. Further analysis of the microarray data suggests that B. thailandensis, when exposed to the aromatic amino acid biosynthesis inhibitor glyphosate, induces (or represses) genes to possibly recuperate and balance the intracellular amino acid pool. We validated our single-cell microarray data at the multi-cell and single-cell levels with lacZ and gfp reporter-gene fusions, respectively. Sanger sequencing of 192 clones generated from the TTA product of a single cell, with and without enrichment by elimination of rRNA and tRNA, detected only B. thailandensis sequences with no contamination. These data indicate that RNA-seq of TTA from a single cell is possible using this novel method.

Characterization of a 5-azacytidine-induced developmental Aspergillus fumigatus variant

The hypomethylating agent 5-azacytidine (5AC) is widely used in patients at risk of invasive mycoses. We sought to determine whether 5AC affects the developmental competence and virulence of Aspergillus fumigatus. Incubation of A. fumigatus strain 293 with 5AC induced high-frequency conversion to a fluffy-variant (Af293 (FL) ). The conidiation defect was bypassed by exposing Af293 (FL) to light during the initial 18 hours of growth on solid media. Transcriptional profiling revealed differential expression of multiple genes involved in G-protein signaling, including a putative G-protein coupled photoreceptor (opsin), suggesting that impaired signaling through a light-responsive pathway upstream of brlA is responsible for this phenotype. Af293 (FL) was fully virulent in fruit fly and murine models of invasive aspergillosis. Moreover, Af293 (FL) overexpressed aspergillopepsin F, had increased elastase activity and was more angioinvasive than the parental wild-type strain. The 5AC-induced A. fumigatus fluffy variant illustrates the potential effects of chemotherapeutic agents on the developmental and pathobiologic characteristics of opportunistic fungi.

The temporal dynamics of differential gene expression in Aspergillus fumigatus interacting with human immature dendritic cells in vitro

Dendritic cells (DC) are the most important antigen presenting cells and play a pivotal role in host immunity to infectious agents by acting as a bridge between the innate and adaptive immune systems. Monocyte-derived immature DCs (iDC) were infected with viable resting conidia of Aspergillus fumigatus (Af293) for 12 hours at an MOI of 5; cells were sampled every three hours. RNA was extracted from both organisms at each time point and hybridised to microarrays. iDC cell death increased at 6 h in the presence of A. fumigatus which coincided with fungal germ tube emergence; >80% of conidia were associated with iDC. Over the time course A. fumigatus differentially regulated 210 genes, FunCat analysis indicated significant up-regulation of genes involved in fermentation, drug transport, pathogenesis and response to oxidative stress. Genes related to cytotoxicity were differentially regulated but the gliotoxin biosynthesis genes were down regulated over the time course, while Aspf1 was up-regulated at 9 h and 12 h. There was an up-regulation of genes in the subtelomeric regions of the genome as the interaction progressed. The genes up-regulated by iDC in the presence of A. fumigatus indicated that they were producing a pro-inflammatory response which was consistent with previous transcriptome studies of iDC interacting with A. fumigatus germ tubes. This study shows that A. fumigatus adapts to phagocytosis by iDCs by utilising genes that allow it to survive the interaction rather than just up-regulation of specific virulence genes.

Comparative genomic analysis of fruiting body formation in Myxococcales

Genetic programs underlying multicellular morphogenesis and cellular differentiation are most often associated with eukaryotic organisms, but examples also exist in bacteria such as the formation of multicellular, spore-filled fruiting bodies in the order Myxococcales. Most members of the Myxococcales undergo a multicellular developmental program culminating in the formation of spore-filled fruiting bodies in response to starvation. To gain insight into the evolutionary history of fruiting body formation in Myxococcales, we performed a comparative analysis of the genomes and transcriptomes of five Myxococcales species, four of these undergo fruiting body formation (Myxococcus xanthus, Stigmatella aurantiaca, Sorangium cellulosum, and Haliangium ochraceum) and one does not (Anaeromyxobacter dehalogenans). Our analyses show that a set of 95 known M. xanthus development-specific genes--although suffering from a sampling bias--are overrepresented and occur more frequently than an average M. xanthus gene in S. aurantiaca, whereas they occur at the same frequency as an average M. xanthus gene in S. cellulosum and in H. ochraceum and are underrepresented in A. dehalogenans. Moreover, genes for entire signal transduction pathways important for fruiting body formation in M. xanthus are conserved in S. aurantiaca, whereas only a minority of these genes are conserved in A. dehalogenans, S. cellulosum, and H. ochraceum. Likewise, global gene expression profiling of developmentally regulated genes showed that genes that upregulated during development in M. xanthus are overrepresented in S. aurantiaca and slightly underrepresented in A. dehalogenans, S. cellulosum, and H. ochraceum. These comparative analyses strongly indicate that the genetic programs for fruiting body formation in M. xanthus and S. aurantiaca are highly similar and significantly different from the genetic program directing fruiting body formation in S. cellulosum and H. ochraceum. Thus, our analyses reveal an unexpected level of plasticity in the genetic programs for fruiting body formation in the Myxococcales and strongly suggest that the genetic program underlying fruiting body formation in different Myxococcales is not conserved. The evolutionary implications of this finding are discussed.

Aspergillus fumigatus AcuM regulates both iron acquisition and gluconeogenesis

Relatively few transcription factors that govern the virulence of Aspergillus fumigatus are known. We constructed 11 A. fumigatus transcription factor mutants and screened them for altered virulence in Galleria mellonella larvae. We discovered that the zinc cluster transcription factor, AcuM, is essential for maximal virulence in this model, as well as in murine models of haematogenously disseminated and invasive pulmonary aspergillosis. Transcriptional profiling experiments suggested that AcuM suppresses sreA and induces hapX to stimulate expression of genes involved in both reductive iron assimilation and siderophore-mediated iron uptake. Consistent with these results, a ΔacuM mutant had reduced iron incorporation, decreased extracellular siderophore production and impaired capacity to grow under iron-limited conditions. Interestingly, an Aspergillus nidulansΔacuM mutant had normal extracellular siderophore production and growth under iron-limited conditions, indicating that AcuM does not govern iron acquisition in this organism. A. fumigatus AcuM also regulated genes involved in gluconeogenesis, and the ΔacuM mutant had impaired growth on gluconeogenic carbon sources. Deletion of sreA in the ΔacuM mutant restored iron uptake, extracellular siderophore production and virulence, but not the defect in gluconeogenesis. Thus, AcuM represses SreA and thereby induces iron acquisition, a process that is essential for the maximal virulence of A. fumigatus.

SMURF: Genomic mapping of fungal secondary metabolite clusters

Fungi produce an impressive array of secondary metabolites (SMs) including mycotoxins, antibiotics and pharmaceuticals. The genes responsible for their biosynthesis, export, and transcriptional regulation are often found in contiguous gene clusters. To facilitate annotation of these clusters in sequenced fungal genomes, we developed the web-based software SMURF (www.jcvi.org/smurf/) to systematically predict clustered SM genes based on their genomic context and domain content. We applied SMURF to catalog putative clusters in 27 publicly available fungal genomes. Comparison with genetically characterized clusters from six fungal species showed that SMURF accurately recovered all clusters and detected additional potential clusters. Subsequent comparative analysis revealed the striking biosynthetic capacity and variability of the fungal SM pathways and the correlation between unicellularity and the absence of SMs. Further genetics studies are needed to experimentally confirm these clusters.

Genetic and phenotypic diversity in Burkholderia: contributions by prophage and phage-like elements

BACKGROUND

Burkholderia species exhibit enormous phenotypic diversity, ranging from the nonpathogenic, soil- and water-inhabiting Burkholderia thailandensis to the virulent, host-adapted mammalian pathogen B. mallei. Genomic diversity is evident within Burkholderia species as well. Individual isolates of Burkholderia pseudomallei and B. thailandensis, for example, carry a variety of strain-specific genomic islands (GIs), including putative pathogenicity and metabolic islands, prophage-like islands, and prophages. These GIs may provide some strains with a competitive advantage in the environment and/or in the host relative to other strains.

RESULTS

Here we present the results of analysis of 37 prophages, putative prophages, and prophage-like elements from six different Burkholderia species. Five of these were spontaneously induced to form bacteriophage particles from B. pseudomallei and B. thailandensis strains and were isolated and fully sequenced; 24 were computationally predicted in sequenced Burkholderia genomes; and eight are previously characterized prophages or prophage-like elements. The results reveal numerous differences in both genome structure and gene content among elements derived from different species as well as from strains within species, due in part to the incorporation of additional DNA, or 'morons' into the prophage genomes. Implications for pathogenicity are also discussed. Lastly, RNAseq analysis of gene expression showed that many of the genes in varphi1026b that appear to contribute to phage and lysogen fitness were expressed independently of the phage structural and replication genes.

CONCLUSIONS

This study provides the first estimate of the relative contribution of prophages to the vast phenotypic diversity found among the Burkholderiae.

The sequence of a 1.8-mb bacterial linear plasmid reveals a rich evolutionary reservoir of secondary metabolic pathways

Plasmids are mobile genetic elements that play a key role in the evolution of bacteria by mediating genome plasticity and lateral transfer of useful genetic information. Although originally considered to be exclusively circular, linear plasmids have also been identified in certain bacterial phyla, notably the actinomycetes. In some cases, linear plasmids engage with chromosomes in an intricate evolutionary interplay, facilitating the emergence of new genome configurations by transfer and recombination or plasmid integration. Genome sequencing of Streptomyces clavuligerus ATCC 27064, a Gram-positive soil bacterium known for its production of a diverse array of biotechnologically important secondary metabolites, revealed a giant linear plasmid of 1.8 Mb in length. This megaplasmid (pSCL4) is one of the largest plasmids ever identified and the largest linear plasmid to be sequenced. It contains more than 20% of the putative protein-coding genes of the species, but none of these is predicted to be essential for primary metabolism. Instead, the plasmid is densely packed with an exceptionally large number of gene clusters for the potential production of secondary metabolites, including a large number of putative antibiotics, such as staurosporine, moenomycin, β-lactams, and enediynes. Interestingly, cross-regulation occurs between chromosomal and plasmid-encoded genes. Several factors suggest that the megaplasmid came into existence through recombination of a smaller plasmid with the arms of the main chromosome. Phylogenetic analysis indicates that heavy traffic of genetic information between Streptomyces plasmids and chromosomes may facilitate the rapid evolution of secondary metabolite repertoires in these bacteria.

The early stage of bacterial genome-reductive evolution in the host

The equine-associated obligate pathogen Burkholderia mallei was developed by reductive evolution involving a substantial portion of the genome from Burkholderia pseudomallei, a free-living opportunistic pathogen. With its short history of divergence (approximately 3.5 myr), B. mallei provides an excellent resource to study the early steps in bacterial genome reductive evolution in the host. By examining 20 genomes of B. mallei and B. pseudomallei, we found that stepwise massive expansion of IS (insertion sequence) elements ISBma1, ISBma2, and IS407A occurred during the evolution of B. mallei. Each element proliferated through the sites where its target selection preference was met. Then, ISBma1 and ISBma2 contributed to the further spread of IS407A by providing secondary insertion sites. This spread increased genomic deletions and rearrangements, which were predominantly mediated by IS407A. There were also nucleotide-level disruptions in a large number of genes. However, no significant signs of erosion were yet noted in these genes. Intriguingly, all these genomic modifications did not seriously alter the gene expression patterns inherited from B. pseudomallei. This efficient and elaborate genomic transition was enabled largely through the formation of the highly flexible IS-blended genome and the guidance by selective forces in the host. The detailed IS intervention, unveiled for the first time in this study, may represent the key component of a general mechanism for early bacterial evolution in the host.

Continuing evolution of Burkholderia mallei through genome reduction and large-scale rearrangements

Burkholderia mallei (Bm), the causative agent of the predominately equine disease glanders, is a genetically uniform species that is very closely related to the much more diverse species Burkholderia pseudomallei (Bp), an opportunistic human pathogen and the primary cause of melioidosis. To gain insight into the relative lack of genetic diversity within Bm, we performed whole-genome comparative analysis of seven Bm strains and contrasted these with eight Bp strains. The Bm core genome (shared by all seven strains) is smaller in size than that of Bp, but the inverse is true for the variable gene sets that are distributed across strains. Interestingly, the biological roles of the Bm variable gene sets are much more homogeneous than those of Bp. The Bm variable genes are found mostly in contiguous regions flanked by insertion sequence (IS) elements, which appear to mediate excision and subsequent elimination of groups of genes that are under reduced selection in the mammalian host. The analysis suggests that the Bm genome continues to evolve through random IS-mediated recombination events, and differences in gene content may contribute to differences in virulence observed among Bm strains. The results are consistent with the view that Bm recently evolved from a single strain of Bp upon introduction into an animal host followed by expansion of IS elements, prophage elimination, and genome rearrangements and reduction mediated by homologous recombination across IS elements.

Aspergillus flavus genomics as a tool for studying the mechanism of aflatoxin formation

Aspergillus flavus is a weak pathogen that infects plants, animals and humans. When it infects agricultural crops, however, it produces one of the most potent carcinogens known (aflatoxins). To devise strategies to control aflatoxin contamination of pre-harvest agricultural crops and post-harvest grains during storage, we launched the A. flavus genomics program. The major objective of this program is the identification of genes involved in aflatoxin biosynthesis and regulation, as well as in pathogenicity, to gain a better understanding of the mechanism of aflatoxin formation. The sequencing of A. flavus whole genome has been completed. Initial annotation of the sequence revealed that there are about 13,071 genes in the A. flavus genome. Genes which potentially encode for enzymes involved in secondary metabolite production in the A. flavus genome have been identified. Preliminary comparative genome analysis of A. flavus with A. oryzae is summarized here.

Genomics of industrial Aspergilli and comparison with toxigenic relatives

Aspergillus oryzae has been used in Japanese fermentation industries for more than a thousand years. The species produces large amounts of various hydrolytic enzymes and has been successfully applied to modern biotechnology. The size of the A. oryzae genome (37.5 Mb) is very close to that of A. flavus and A. niger, and 20-30% larger than that of either A. nidulans or A. fumigatus. A. oryzae and A. flavus have exactly the same number of aspartic proteinase genes, of which each orthologous pair shares highly conserved amino acid sequences. Synteny analysis with A. fumigatus and A. nidulans showed that the A. oryzae genome has a mosaic structure consisting of syntenic and non-syntenic blocks. In the microorganisms to be compared, the density of the genes having homologs was obviously higher on the syntenic than on the non-syntenic blocks. Expression analysis by the DNA microarray supported the significantly lower expression of genes on the non-syntenic than on the syntenic blocks.

Effect of competition on the production and activity of secondary metabolites in Aspergillus species

Secondary metabolites are of intense interest to humans due to their pharmaceutical and/or toxic properties. Also, these metabolites are clinically relevant because of their importance in fungal pathogenesis. Aspergillus species secrete secondary metabolites when grown individually and in the presence of other fungal species. However, it is not known whether secreted secondary metabolites provide a competitive advantage over other fungal species, or whether competition has any effect on the production of those metabolites. Here, we have performed co-cultivation competition assays among different species of Aspergillus to determine relative species fitness in culture, and to analyze the presence of possible antifungal activity of secondary metabolites in extracts. The results show that, for the most part, at 30 degrees C only one species is able to survive direct competition with a second species. In contrast, survival of both competitors was often observed at 37 degrees C. Consistent with these observations, antifungal activity of extracts from cultures grown at 30 degrees C was greater than that of extract from cultures at 37 degrees C. Interestingly, culture extracts from all species studied had some degree of antifungal activity, but in general, the extracts had greater antifungal activity when species were grown in the presence of a competitor. Using gas chromatography it was determined that the composition of extracts changed due to competition and a shift in temperature. These findings indicate that co-cultivation could be a very promising method for inducing and characterizing novel antifungal compounds produced by species of Aspergillus.

Simple sequence repeat (SSR)-based gene diversity in Burkholderia pseudomallei and Burkholderia mallei

Pathogens Burkholderia pseudomallei (Bp) and Burkholderia mallei (Bm) contain a large number (> 12,000) of Simple Sequence Repeats (SSRs). To study the extent to which these features have contributed to the diversification of genes, we have conducted comparative studies with nineteen genomes of these bacteria. We found 210 genes with characteristic types of SSR variations. SSRs with nonamer repeat units were the most abundant, followed by hexamers and trimers. Amino acids with smaller and nonpolar R-groups are preferred to be encoded by the variant SSRs, perhaps due to their minimal impacts to protein functionality. A majority of these genes appears to code for surface or secreted proteins that may directly interact with the host factors during pathogenesis or other environmental factors. There also are others that encode diverse functions in the cytoplasm, and this protein variability may reflect an extensive involvement of phase variation in survival and adaptation of these pathogens.

Exogenous Yersinia pestis quorum sensing molecules N-octanoyl-homoserine lactone and N-(3-oxooctanoyl)-homoserine lactone regulate the LcrV virulence factor

LcrV is a key Yersinia pestis antigen, immune regulator, and component of the type III secretion system (T3SS). Researchers have shown that N-acyl-homoserine lactones (AHLs) can down-regulate the expression of the LcrV homolog, PcrV, in Pseudomonas aeruginosa. Using ELISA, western blot, DNA microarray analysis, and real time PCR we demonstrate that the addition of AHL molecules N-octanoyl-homoserine lactone (C8) or N-(3-oxooctanoyl)-homoserine lactone (oxo-C8) to Y. pestis cultures down-regulates LcrV protein expression. DNA microarray analysis shows 10 additional T3SS genes are consistently down-regulated by C8 or oxo-C8. This is the first report demonstrating that AHLs regulate Y. pestis virulence factor expression.