Search

Abstract 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

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.

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.

Aspergillus flavus and A. parasiticus infect peanut seeds and produce aflatoxins, which are associated with various diseases in domestic animals and humans throughout the world. The most cost-effective strategy to minimize aflatoxin contamination involves the development of peanut cultivars that are resistant to fungal infection and/or aflatoxin production. To identify peanut Aspergillus-interactive and peanut Aspergillus-resistance genes, we carried out a large scale peanut Expressed Sequence Tag (EST) project which we used to construct a peanut glass slide oligonucleotide microarray. The fabricated microarray represents over 40% of the protein coding genes in the peanut genome. For expression profiling, resistant and susceptible peanut cultivars were infected with a mixture of Aspergillus flavus and parasiticus spores. The subsequent microarray analysis identified 62 genes in resistant cultivars that were up-expressed in response to Aspergillus infection. In addition, we identified 22 putative Aspergillus-resistance genes that were constitutively up-expressed in the resistant cultivar in comparison to the susceptible cultivar. Some of these genes were homologous to peanut, corn, and soybean genes that were previously shown to confer resistance to fungal infection. This study is a first step towards a comprehensive genome-scale platform for developing Aspergillus-resistant peanut cultivars through targeted marker-assisted breeding and genetic engineering.

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.

ABSTRACT Dermatophytes are a uniquely pathogenic group of fungi that cause most common fungal infections globally. The major cause of athlete’s foot is Trichophyton rubrum, a pathogen of human skin. A recent paper in this journal reported the sequencing and analysis of five additional genome sequences, including that of Trichophyton rubrum. These five join the existing two additional genome sequences to bring the total to seven dermatophyte genome sequences, a notable milestone in the study of these fungi. These additional genomes set the stage for future genome-supported studies on the biology, pathogenicity, and host specificity of this important group of pathogens. To predict how this future might play out, we review the history of Aspergillus genomics since the initial publication of the first three Aspergillus genome sequences in 2005, an event that stimulated important studies of the pathogenic Aspergillus species. From these 7 years of Aspergillus history, we offer some speculation on the future of dermatophyte studies supported by the genome sequences given the similarities, differences, and relative levels of support for studies in these two groups of fungi and the diseases they cause.

ABSTRACT The mating efficiency of 50 Aspergillus fumigatus isolates from both clinical and environmental sources was analyzed. Forty isolates completed the sexual cycle in 4 weeks with variable levels of fertility designated high, medium, or low. Two opposite-mating-type strains exhibiting the highest fertility, AFB62 (MAT1-1), isolated from a case of invasive aspergillosis, and AFIR928 (MAT1-2), isolated from the environment, were chosen as the supermater pair. Single cleistothecia obtained from a cross of the two strains harbored a minimum of 1 × 104 ascospores. The viability of ascospores increased with the age of the fruiting body, 17% at 4 weeks and reaching 95% at 20 weeks. AFB62 and AFIR928 were equally virulent in two different murine models, despite differences in their sources. High recombination frequencies were observed when the closely linked genes alb1 (AFUA_2G17600) and abr2 (AFUA_2G17530) were used as genetic markers. Comparative genome hybridization analyses revealed that only 86 genes (ca. 0.86% of the genome) are significantly diverged between AFB62 and AFIR928. The high fertility in a relatively short period, combined with a high degree of virulence and a high recombination frequency, demonstrates that the mating pair AFB62 and AFIR928 provides an excellent tool for genetic studies of A. fumigatus. IMPORTANCE Aspergillus fumigatus is a heterothallic fungal pathogen that causes life-threatening infections in immunocompromised hosts. Although heterothallism facilitates genetic study via recombinational analysis, previous work showed that a 6-month incubation period is required for the completion of sexual reproduction in this species. Such a long incubation period impedes progress in genetic research. To discover a highly fertile (supermater) pair that can complete the sexual cycle in a considerably shorter period, we screened 50 strains collected from various geographic regions for mating efficiency. We identified a highly virulent pair of supermaters that can be an invaluable tool for genetic study.

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.

IntroductionNasopharyngeal (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.MethodsAs 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.ResultsAMR 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.ConclusionThis 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.

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.

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.

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.

Destruction of the pulmonary epithelium is a major feature of lung diseases caused by the mould pathogen Aspergillus fumigatus. Although it is widely postulated that tissue invasion is governed by fungal proteases, A. fumigatus mutants lacking individual or multiple enzymes remain fully invasive, suggesting a concomitant requirement for other pathogenic activities during host invasion. In this study we discovered, and exploited, a novel, tissue non-invasive, phenotype in A. fumigatus mutants lacking the pH-responsive transcription factor PacC. Our study revealed a novel mode of epithelial entry, occurring in a cell wall-dependent manner prior to protease production, and via the Dectin-1 β-glucan receptor. ΔpacC mutants are defective in both contact-mediated epithelial entry and protease expression, and significantly attenuated for pathogenicity in leukopenic mice. We combined murine infection modelling, in vivo transcriptomics, and in vitro infections of human alveolar epithelia, to delineate two major, and sequentially acting, PacC-dependent processes impacting epithelial integrity in vitro and tissue invasion in the whole animal. We demonstrate that A. fumigatus spores and germlings are internalised by epithelial cells in a contact-, actin-, cell wall- and Dectin-1 dependent manner and ΔpacC mutants, which aberrantly remodel the cell wall during germinative growth, are unable to gain entry into epithelial cells, both in vitro and in vivo. We further show that PacC acts as a global transcriptional regulator of secreted molecules during growth in the leukopenic mammalian lung, and profile the full cohort of secreted gene products expressed during invasive infection. Our study reveals a combinatorial mode of tissue entry dependent upon sequential, and mechanistically distinct, perturbations of the pulmonary epithelium and demonstrates, for the first time a protective role for Dectin-1 blockade in epithelial defences. Infecting ΔpacC mutants are hypersensitive to cell wall-active antifungal agents highlighting the value of PacC signalling as a target for antifungal therapy.

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.

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.

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.

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.

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.

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(i), 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.

Yersinia pestis is the causative agent of the plague. Y. pestis KIM 10+ strain was passaged and selected for loss of the 102 kb pgm locus, resulting in an attenuated strain, KIM D27. In this study, whole genome sequencing was performed on KIM D27 in order to identify any additional differences. Initial assemblies of 454 data were highly fragmented, and various bioinformatic tools detected between 15 and 465 SNPs and INDELs when comparing both strains, the vast majority associated with A or T homopolymer sequences. Consequently, Illumina sequencing was performed to improve the quality of the assembly. Hybrid sequence assemblies were performed and a total of 56 validated SNP/INDELs and 5 repeat differences were identified in the D27 strain relative to published KIM 10+ sequence. However, further analysis showed that 55 of these SNP/INDELs and 3 repeats were errors in the KIM 10+ reference sequence. We conclude that both 454 and Illumina sequencing were required to obtain the most accurate and rapid sequence results for Y. pestis KIMD27. SNP and INDELS calls were most accurate when both Newbler and CLC Genomics Workbench were employed. For purposes of obtaining high quality genome sequence differences between strains, any identified differences should be verified in both the new and reference genomes.

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.

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.

Certain environmental microorganisms can cause severe human infections, even in the absence of an obvious requirement for transition through an animal host for replication ("accidental virulence"). To understand this process, we compared eleven isolate genomes of Burkholderia pseudomallei (Bp), a tropical soil microbe and causative agent of the human and animal disease melioidosis. We found evidence for the existence of several new genes in the Bp reference genome, identifying 282 novel genes supported by at least two independent lines of supporting evidence (mRNA transcripts, database homologs, and presence of ribosomal binding sites) and 81 novel genes supported by all three lines. Within the Bp core genome, 211 genes exhibited significant levels of positive selection (4.5%), distributed across many cellular pathways including carbohydrate and secondary metabolism. Functional experiments revealed that certain positively selected genes might enhance mammalian virulence by interacting with host cellular pathways or utilizing host nutrients. Evolutionary modifications improving Bp environmental fitness may thus have indirectly facilitated the ability of Bp to colonize and survive in mammalian hosts. These findings improve our understanding of the pathogenesis of melioidosis, and establish Bp as a model system for studying the genetics of accidental virulence.

Aspergillus fumigatus is a common mould whose spores are a component of the normal airborne flora. Immune dysfunction permits developmental growth of inhaled spores in the human lung causing aspergillosis, a significant threat to human health in the form of allergic, and life-threatening invasive infections. The success of A. fumigatus as a pathogen is unique among close phylogenetic relatives and is poorly characterised at the molecular level. Recent genome sequencing of several Aspergillus species provides an exceptional opportunity to analyse fungal virulence attributes within a genomic and evolutionary context. To identify genes preferentially expressed during adaptation to the mammalian host niche, we generated multiple gene expression profiles from minute samplings of A. fumigatus germlings during initiation of murine infection. They reveal a highly co-ordinated A. fumigatus gene expression programme, governing metabolic and physiological adaptation, which allows the organism to prosper within the mammalian niche. As functions of phylogenetic conservation and genetic locus, 28% and 30%, respectively, of the A. fumigatus subtelomeric and lineage-specific gene repertoires are induced relative to laboratory culture, and physically clustered genes including loci directing pseurotin, gliotoxin and siderophore biosyntheses are a prominent feature. Locationally biased A. fumigatus gene expression is not prompted by in vitro iron limitation, acid, alkaline, anaerobic or oxidative stress. However, subtelomeric gene expression is favoured following ex vivo neutrophil exposure and in comparative analyses of richly and poorly nourished laboratory cultured germlings. We found remarkable concordance between the A. fumigatus host-adaptation transcriptome and those resulting from in vitro iron depletion, alkaline shift, nitrogen starvation and loss of the methyltransferase LaeA. This first transcriptional snapshot of a fungal genome during initiation of mammalian infection provides the global perspective required to direct much-needed diagnostic and therapeutic strategies and reveals genome organisation and subtelomeric diversity as potential driving forces in the evolution of pathogenicity in the genus Aspergillus.

We present the genome sequences of a new clinical isolate of the important human pathogen, Aspergillus fumigatus, A1163, and two closely related but rarely pathogenic species, Neosartorya fischeri NRRL181 and Aspergillus clavatus NRRL1. Comparative genomic analysis of A1163 with the recently sequenced A. fumigatus isolate Af293 has identified core, variable and up to 2% unique genes in each genome. While the core genes are 99.8% identical at the nucleotide level, identity for variable genes can be as low 40%. The most divergent loci appear to contain heterokaryon incompatibility (het) genes associated with fungal programmed cell death such as developmental regulator rosA. Cross-species comparison has revealed that 8.5%, 13.5% and 12.6%, respectively, of A. fumigatus, N. fischeri and A. clavatus genes are species-specific. These genes are significantly smaller in size than core genes, contain fewer exons and exhibit a subtelomeric bias. Most of them cluster together in 13 chromosomal islands, which are enriched for pseudogenes, transposons and other repetitive elements. At least 20% of A. fumigatus-specific genes appear to be functional and involved in carbohydrate and chitin catabolism, transport, detoxification, secondary metabolism and other functions that may facilitate the adaptation to heterogeneous environments such as soil or a mammalian host. Contrary to what was suggested previously, their origin cannot be attributed to horizontal gene transfer (HGT), but instead is likely to involve duplication, diversification and differential gene loss (DDL). The role of duplication in the origin of lineage-specific genes is further underlined by the discovery of genomic islands that seem to function as designated "gene dumps" and, perhaps, simultaneously, as "gene factories".

Burkholderia mallei is an understudied biothreat agent responsible for glanders which can be lethal in humans and animals. Research with this pathogen has been hampered in part by constraints of Select Agent regulations for safety reasons. Whole genomic sequencing (WGS) is an apt approach to characterize newly discovered or poorly understood microbial pathogens.We performed WGS on a strain of B. mallei, SAVP1, previously pathogenic, that was experimentally infected in 6 equids (4 ponies, 1 mule, 1 donkey), natural hosts, for purposes of producing antibodies. Multiple high inocula were used in some cases. Unexpectedly SAVP1 appeared to be avirulent in the ponies and mule, and attenuated in the donkey, but induced antibodies. We determined the genome sequence of SAVP1 and compared it to a strain that was virulent in horses and a human. In comparison, this phenotypic avirulent SAVP1 strain was missing multiple genes including all the animal type III secretory system (T3SS) complex of genes demonstrated to be essential for virulence in mice and hamster models. The loss of these genes in the SAVP1 strain appears to be the consequence of a multiple gene deletion across insertion sequence (IS) elements in the B. mallei genome. Therefore, the strain by itself is unlikely to revert naturally to its virulent phenotype. There were other genes present in one strain and not the other and vice-versa.The discovery that this strain of B. mallei was both avirulent in the natural host ponies, and did not possess T3SS associated genes may be fortuitous to advance biodefense research. The deleted virulence-essential T3SS is not likely to be re-acquired naturally. These findings may provide a basis for exclusion of SAVP1 from the Select Agent regulation or at least discussion of what else would be required for exclusion. This exclusion could accelerate research by investigators not possessing BSL-3 facilities and facilitate the production of reagents such as antibodies without the restraints of Select Agent regulation.

Secondary metabolites, including toxins and melanins, have been implicated as virulence attributes in invasive aspergillosis. Although not definitively proved, this supposition is supported by the decreased virulence of an Aspergillus fumigatus strain, DeltalaeA, that is crippled in the production of numerous secondary metabolites. However, loss of a single LaeA-regulated toxin, gliotoxin, did not recapitulate the hypovirulent DeltalaeA pathotype, thus implicating other toxins whose production is governed by LaeA. Toward this end, a whole-genome comparison of the transcriptional profile of wild-type, DeltalaeA, and complemented control strains showed that genes in 13 of 22 secondary metabolite gene clusters, including several A. fumigatus-specific mycotoxin clusters, were expressed at significantly lower levels in the DeltalaeA mutant. LaeA influences the expression of at least 9.5% of the genome (943 of 9,626 genes in A. fumigatus) but positively controls expression of 20% to 40% of major classes of secondary metabolite biosynthesis genes such as nonribosomal peptide synthetases (NRPSs), polyketide synthases, and P450 monooxygenases. Tight regulation of NRPS-encoding genes was highlighted by quantitative real-time reverse-transcription PCR analysis. In addition, expression of a putative siderophore biosynthesis NRPS (NRPS2/sidE) was greatly reduced in the DeltalaeA mutant in comparison to controls under inducing iron-deficient conditions. Comparative genomic analysis showed that A. fumigatus secondary metabolite gene clusters constitute evolutionarily diverse regions that may be important for niche adaptation and virulence attributes. Our findings suggest that LaeA is a novel target for comprehensive modification of chemical diversity and pathogenicity.

The complete sequence of the 1,267,782 bp genome of Wolbachia pipientis wMel, an obligate intracellular bacteria of Drosophila melanogaster, has been determined. Wolbachia, which are found in a variety of invertebrate species, are of great interest due to their diverse interactions with different hosts, which range from many forms of reproductive parasitism to mutualistic symbioses. Analysis of the wMel genome, in particular phylogenomic comparisons with other intracellular bacteria, has revealed many insights into the biology and evolution of wMel and Wolbachia in general. For example, the wMel genome is unique among sequenced obligate intracellular species in both being highly streamlined and containing very high levels of repetitive DNA and mobile DNA elements. This observation, coupled with multiple evolutionary reconstructions, suggests that natural selection is somewhat inefficient in wMel, most likely owing to the occurrence of repeated population bottlenecks. Genome analysis predicts many metabolic differences with the closely related Rickettsia species, including the presence of intact glycolysis and purine synthesis, which may compensate for an inability to obtain ATP directly from its host, as Rickettsia can. Other discoveries include the apparent inability of wMel to synthesize lipopolysaccharide and the presence of the most genes encoding proteins with ankyrin repeat domains of any prokaryotic genome yet sequenced. Despite the ability of wMel to infect the germline of its host, we find no evidence for either recent lateral gene transfer between wMel and D. melanogaster or older transfers between Wolbachia and any host. Evolutionary analysis further supports the hypothesis that mitochondria share a common ancestor with the alpha-Proteobacteria, but shows little support for the grouping of mitochondria with species in the order Rickettsiales. With the availability of the complete genomes of both species and excellent genetic tools for the host, the wMel-D. melanogaster symbiosis is now an ideal system for studying the biology and evolution of Wolbachia infections.