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.

Transcriptome Analysis of Aspergillus flavus Reveals veA-Dependent Regulation of Secondary Metabolite Gene Clusters, Including the Novel Aflavarin Cluster.

The global regulatory veA gene governs development and secondary metabolism in numerous fungal species, including Aspergillus flavus. This is especially relevant since A. flavus infects crops of agricultural importance worldwide, contaminating them with potent mycotoxins. The most well-known are aflatoxins, which are cytotoxic and carcinogenic polyketide compounds. The production of aflatoxins and the expression of genes implicated in the production of these mycotoxins are veA dependent. The genes responsible for the synthesis of aflatoxins are clustered, a signature common for genes involved in fungal secondary metabolism. Studies of the A. flavus genome revealed many gene clusters possibly connected to the synthesis of secondary metabolites. Many of these metabolites are still unknown, or the association between a known metabolite and a particular gene cluster has not yet been established. In the present transcriptome study, we show that veA is necessary for the expression of a large number of genes. Twenty-eight out of the predicted 56 secondary metabolite gene clusters include at least one gene that is differentially expressed depending on presence or absence of veA. One of the clusters under the influence of veA is cluster 39. The absence of veA results in a downregulation of the five genes found within this cluster. Interestingly, our results indicate that the cluster is expressed mainly in sclerotia. Chemical analysis of sclerotial extracts revealed that cluster 39 is responsible for the production of aflavarin.

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.

Aflatoxin formation and gene expression in response to carbon source media shift in Aspergillus parasiticus.

Aflatoxins are toxic and carcinogenic polyketide metabolites produced by fungal species, including Aspergillus flavus and A. parasiticus. The biosynthesis of aflatoxins is modulated by many environmental factors, including the availability of a carbon source. The gene expression profile of A. parasiticus was evaluated during a shift from a medium with low concentration of simple sugars, yeast extract (YE), to a similar medium with sucrose, yeast extract sucrose (YES). Gene expression and aflatoxins (B1, B2, G1, and G2) were quantified from fungal mycelia harvested pre- and post-shifting. When compared with YE media, YES caused temporary reduction of the aflatoxin levels detected at 3-h post-shifting and they remained low well past 12 h post-shift. Aflatoxin levels did not exceed the levels in YE until 24 h post-shift, at which time point a tenfold increase was observed over YE. Microarray analysis comparing the RNA samples from the 48-h YE culture to the YES samples identified a total of 2120 genes that were expressed across all experiments, including most of the aflatoxin biosynthesis genes. One-way analysis of variance (ANOVA) identified 56 genes that were expressed with significant variation across all time points. Three genes responsible for converting norsolorinic acid to averantin were identified among these significantly expressed genes. The potential involvement of these genes in the regulation of aflatoxin biosynthesis is discussed.

Amino acid supplementation reveals differential regulation of aflatoxin biosynthesis in Aspergillus flavus NRRL 3357 and Aspergillus parasiticus SRRC 143.

Aflatoxins are toxic and carcinogenic secondary metabolites produced by the fungi Aspergillus flavus and Aspergillus parasiticus. To better understand the molecular mechanisms that regulate aflatoxin production, the biosynthesis of the toxin in A. flavus and A. parasticus grown in yeast extract sucrose media supplemented with 50 mM tryptophan (Trp) were examined. Aspergillus flavus grown in the presence of 50 mM tryptophan was found to have significantly reduced aflatoxin B(1) and B(2) biosynthesis, while A. parasiticus cultures had significantly increased B(1) and G(1) biosynthesis. Microarray analysis of RNA extracted from fungi grown under these conditions revealed 77 genes that are expressed significantly different between A. flavus and A. parasiticus, including the aflatoxin biosynthetic genes aflD (nor-1), aflE (norA), and aflO (omtB). It is clear that the regulatory mechanisms of aflatoxin biosynthesis in response to Trp in A. flavus and A. parasiticus are different. These candidate genes may serve as regulatory factors of aflatoxin biosynthesis.

What can comparative genomics tell us about species concepts in the genus Aspergillus?

Understanding the nature of species" boundaries is a fundamental question in evolutionary biology. The availability of genomes from several species of the genus Aspergillus allows us for the first time to examine the demarcation of fungal species at the whole-genome level. Here, we examine four case studies, two of which involve intraspecific comparisons, whereas the other two deal with interspecific genomic comparisons between closely related species. These four comparisons reveal significant variation in the nature of species boundaries across Aspergillus. For example, comparisons between A. fumigatus and Neosartorya fischeri (the teleomorph of A. fischerianus) and between A. oryzae and A. flavus suggest that measures of sequence similarity and species-specific genes are significantly higher for the A. fumigatus - N. fischeri pair. Importantly, the values obtained from the comparison between A. oryzae and A. flavus are remarkably similar to those obtained from an intra-specific comparison of A. fumigatus strains, giving support to the proposal that A. oryzae represents a distinct ecotype of A. flavus and not a distinct species. We argue that genomic data can aid Aspergillus taxonomy by serving as a source of novel and unprecedented amounts of comparative data, as a resource for the development of additional diagnostic tools, and finally as a knowledge database about the biological differences between strains and species.

Evolution of sensory complexity recorded in a myxobacterial genome.

Myxobacteria are single-celled, but social, eubacterial predators. Upon starvation they build multicellular fruiting bodies using a developmental program that progressively changes the pattern of cell movement and the repertoire of genes expressed. Development terminates with spore differentiation and is coordinated by both diffusible and cell-bound signals. The growth and development of Myxococcus xanthus is regulated by the integration of multiple signals from outside the cells with physiological signals from within. A collection of M. xanthus cells behaves, in many respects, like a multicellular organism. For these reasons M. xanthus offers unparalleled access to a regulatory network that controls development and that organizes cell movement on surfaces. The genome of M. xanthus is large (9.14 Mb), considerably larger than the other sequenced delta-proteobacteria. We suggest that gene duplication and divergence were major contributors to genomic expansion from its progenitor. More than 1,500 duplications specific to the myxobacterial lineage were identified, representing >15% of the total genes. Genes were not duplicated at random; rather, genes for cell-cell signaling, small molecule sensing, and integrative transcription control were amplified selectively. Families of genes encoding the production of secondary metabolites are overrepresented in the genome but may have been received by horizontal gene transfer and are likely to be important for predation.

Aspergillus flavus expressed sequence tags and microarray as tools in understanding aflatoxin biosynthesis.

Aflatoxins are the most toxic and carcinogenic naturally occurring mycotoxins. They are produced primarily byAspergillus flavus andA. parasiticus. In order to better understand the molecular mechanisms that control aflatoxin production, identification of genes usingA. flavus expressed sequence tags (ESTs) and microarrays is currently being performed. Sequencing and annotation ofA. flavus ESTs from a normalizedA. flavus cDNA library identified 7,218 unique EST sequences. Genes that are putatively involved in aflatoxin biosynthesis, regulation and signal transduction, fungal virulence or pathogenicity, stress response or antioxidation, and fungal development were identified from these ESTs. Microarrays containing over 5,000 uniqueA. flavus gene amplicons were constructed at The Institute for Genomic Research. Gene expression profiling under aflatoxin-producing and non-producing conditions using this microarray has identified hundreds of genes that are potentially involved in aflatoxin production. Further investigations on the functions of these genes by gene knockout experiments are underway. This research is expected to provide information for developing new strategies for controlling aflatoxin contamination of agricultural commodities.

Whole genome comparison of the A. fumigatus family.

The availability of the genome sequences of multiple Aspergillus spp. presents the research community with an unprecedented opportunity for discovery. The genomes of Neosartorya fischeri and Aspergillus clavatus have been sequenced in order to extend our knowledge of Aspergillus fumigatus, the primary cause of invasive aspergillosis. Through comparative genomic analysis, we hope to elucidate both obvious and subtle differences between genomes, developing new hypotheses that can be tested in the laboratory. A preliminary examination of the genomes and their predicted proteomes reveals extensive conservation between protein sequences and significant synteny, or conserved gene order. Comparative genomic analysis at the level of these closely related aspergilli should provide important insight into the evolutionary forces at play and their effect on gene content, regulation and expression.

Whole genome comparison of Aspergillus flavus and A. oryzae.

Aspergillus flavus is a plant and animal pathogen that also produces the potent carcinogen aflatoxin. Aspergillus oryzae is a closely related species that has been used for centuries in the food fermentation industry and is Generally Regarded As Safe (GRAS). Whole genome sequences for these two fungi are now complete, providing us with the opportunity to examine any genomic differences that may explain the different ecological niches of these two fungi, and perhaps to identify pathogenicity factors in A. flavus. These two fungi are very similar in genome size and number of predicted genes. The estimated genome size (36·8 Mb) and predicted number of genes (12 197) for A. flavus is similar to that of A. oryzae (36·7 Mb and 12 079, respectively). These two fungi have significantly larger genomes than Aspergillus nidulans (30·1) and Aspergillus fumigatus (29·4). The A. flavus and A. oryzae genomes are enriched in genes for secondary metabolism, but do not differ greatly from one another in the predicted number of polyketide synthases, nonribosomal peptide synthases or the number of genes coding for cytochrome P450 enzymes. A micro-scale analysis of the two fungi did show differences in DNA correspondence between the two species and in the number of transposable elements. Each species has approximately 350 unique genes. The high degree of sequence similarity between the two fungi suggests that they may be ecotypes of the same species and that A. oryzae has resulted from the domestication of A. flavus.

What the Aspergillus genomes have told us.

The sequencing and annotation of the genomes of the first strains of Aspergillus nidulans, Aspergillus oryzae, and Aspergillus fumigatus will be seen in retrospect as a transformational event in Aspergillus biology. With this event the entire genetic composition of A. nidulans, the sexual experimental model organism of the genus Aspergillus, A. oryzae, the food biotechnology organism which is the product of centuries of cultivation, and A. fumigatus, the most common causative agent of invasive aspergillosis is now revealed to the extent that we are at present able to understand. Each genome exhibits a large set of genes common to the three as well as a much smaller set of genes unique to each. Moreover, these sequences serve as resources providing the major tool to expanding our understanding of the biology of each. Transcription profiling of A. fumigatus at high temperatures and comparative genomic hybridization between A. fumigatus and a closely related Aspergillus species provides microarray based examples of the beginning of functional analysis of the genomes of these organisms going forward from the genome sequence.

Mouse BAC ends quality assessment and sequence analyses.

A large-scale BAC end-sequencing project at The Institute for Genomic Research (TIGR) has generated one of the most extensive sets of sequence markers for the mouse genome to date. With a sequencing success rate of >80%, an average read length of 485 bp, and ABI3700 capillary sequencers, we have generated 449,234 nonredundant mouse BAC end sequences (mBESs) with 218 Mb total from 257,318 clones from libraries RPCI-23 and RPCI-24, representing 15x clone coverage, 7% sequence coverage, and a marker every 7 kb across the genome. A total of 191,916 BACs have sequences from both ends providing 12x genome coverage. The average Q20 length is 406 bp and 84% of the bases have phred quality scores > or = 20. RPCI-24 mBESs have more Q20 bases and longer reads on average than RPCI-23 sequences. ABI3700 sequencers and the sample tracking system ensure that > 95% of mBESs are associated with the right clone identifiers. We have found that a significant fraction of mBESs contains L1 repeats and approximately 48% of the clones have both ends with > or = 100 bp contiguous unique Q20 bases. About 3% mBESs match ESTs and > 70% of matches were conserved between the mouse and the human or the rat. Approximately 0.1% mBESs contain STSs. About 0.2% mBESs match human finished sequences and > 70% of these sequences have EST hits. The analyses indicate that our high-quality mouse BAC end sequences will be a valuable resource to the community.

Complete genome sequence of Caulobacter crescentus.

The complete genome sequence of Caulobacter crescentus was determined to be 4,016,942 base pairs in a single circular chromosome encoding 3,767 genes. This organism, which grows in a dilute aquatic environment, coordinates the cell division cycle and multiple cell differentiation events. With the annotated genome sequence, a full description of the genetic network that controls bacterial differentiation, cell growth, and cell cycle progression is within reach. Two-component signal transduction proteins are known to play a significant role in cell cycle progression. Genome analysis revealed that the C. crescentus genome encodes a significantly higher number of these signaling proteins (105) than any bacterial genome sequenced thus far. Another regulatory mechanism involved in cell cycle progression is DNA methylation. The occurrence of the recognition sequence for an essential DNA methylating enzyme that is required for cell cycle regulation is severely limited and shows a bias to intergenic regions. The genome contains multiple clusters of genes encoding proteins essential for survival in a nutrient poor habitat. Included are those involved in chemotaxis, outer membrane channel function, degradation of aromatic ring compounds, and the breakdown of plant-derived carbon sources, in addition to many extracytoplasmic function sigma factors, providing the organism with the ability to respond to a wide range of environmental fluctuations. C. crescentus is, to our knowledge, the first free-living alpha-class proteobacterium to be sequenced and will serve as a foundation for exploring the biology of this group of bacteria, which includes the obligate endosymbiont and human pathogen Rickettsia prowazekii, the plant pathogen Agrobacterium tumefaciens, and the bovine and human pathogen Brucella abortus.

DNA sequence of both chromosomes of the cholera pathogen Vibrio cholerae.

Here we determine the complete genomic sequence of the gram negative, gamma-Proteobacterium Vibrio cholerae El Tor N16961 to be 4,033,460 base pairs (bp). The genome consists of two circular chromosomes of 2,961,146 bp and 1,072,314 bp that together encode 3,885 open reading frames. The vast majority of recognizable genes for essential cell functions (such as DNA replication, transcription, translation and cell-wall biosynthesis) and pathogenicity (for example, toxins, surface antigens and adhesins) are located on the large chromosome. In contrast, the small chromosome contains a larger fraction (59%) of hypothetical genes compared with the large chromosome (42%), and also contains many more genes that appear to have origins other than the gamma-Proteobacteria. The small chromosome also carries a gene capture system (the integron island) and host 'addiction' genes that are typically found on plasmids; thus, the small chromosome may have originally been a megaplasmid that was captured by an ancestral Vibrio species. The V. cholerae genomic sequence provides a starting point for understanding how a free-living, environmental organism emerged to become a significant human bacterial pathogen.

Genome data: what do we learn?

Genome sequence information has continued to accumulate at a spectacular pace during the past year. Details of the sequence and gene content of human chromosome 22 were published. The sequencing and annotation of the first two Arabidopsis thaliana chromosomes was completed. The sequence of chromosome 3 from Plasmodium falciparum, the second sequenced malaria chromosome, was reported, as was that of chromosome 1 from Leishmania major. The complete genomic sequences of five microbes were reported. Approaches to using data from completely sequenced microbial genomes in phylogenetic studies are being explored, as is the application of microarrays to whole genome expression analysis.

Sequence and analysis of chromosome 2 of the plant Arabidopsis thaliana.

Arabidopsis thaliana (Arabidopsis) is unique among plant model organisms in having a small genome (130-140 Mb), excellent physical and genetic maps, and little repetitive DNA. Here we report the sequence of chromosome 2 from the Columbia ecotype in two gap-free assemblies (contigs) of 3.6 and 16 megabases (Mb). The latter represents the longest published stretch of uninterrupted DNA sequence assembled from any organism to date. Chromosome 2 represents 15% of the genome and encodes 4,037 genes, 49% of which have no predicted function. Roughly 250 tandem gene duplications were found in addition to large-scale duplications of about 0.5 and 4.5 Mb between chromosomes 2 and 1 and between chromosomes 2 and 4, respectively. Sequencing of nearly 2 Mb within the genetically defined centromere revealed a low density of recognizable genes, and a high density and diverse range of vestigial and presumably inactive mobile elements. More unexpected is what appears to be a recent insertion of a continuous stretch of 75% of the mitochondrial genome into chromosome 2.

The human proteinase-activated receptor-3 (PAR-3) gene. Identification within a Par gene cluster and characterization in vascular endothelial cells and platelets.

Proteolytically activated receptors (PARs) represent an emerging subset of seven transmembrane G protein-coupled receptors that mediate cell activation events by receptor cleavage at distinct scissile bonds located within receptor amino termini. Differential genomic blotting using a yeast artificial chromosome known to contain the PAR-1 and PAR-2 genes identified the PAR-3 gene within a PAR gene cluster spanning approximately 100 kilobases at 5q13. The PAR-3 gene is relatively small (approximately 12 kilobases); and, like the PAR-1 and PAR-2 genes, it displays a two-exon structure, with the majority of the coding sequence and the proteolytic cleavage site contained within the larger second exon. Sequence analysis of the 5'-flanking region demonstrates that the promoter is TATA-less, similar to that seen with PAR-1, with the identification of nucleic acid motifs potentially involved in transcriptional gene regulation, including AP-1, GATA, and octameric sequences. PAR-3 transcripts were apparent in human vascular endothelial cells, although at considerably lower levels than those of PAR-1 and not significantly modulated by the endothelial cell stimulus tumor necrosis factor-alpha. Likewise, although PAR-3 mRNA was evident in human platelets, receptor cell surface expression was modest (approximately 10%) compared with that of PAR-1. Thus, although PAR-3 is postulated to represent a second thrombin receptor, its modest endothelial cell and platelet expression suggest that PAR-3 activation by alpha-thrombin is less relevant for physiological responses in these mature cells. Rather, given its disparately greater expression in megakaryocytes (and megakaryocyte-like human erythroleukemia cells), a regulatory role in cellular development (by protease activation) could be postulated.

The human thrombin receptor and proteinase activated receptor-2 genes are tightly linked on chromosome 5q13.

The thrombin receptor (TR) and proteinase activated receptor-2 (PAR-2) may represent the prototypes of an emerging family of cell-surface receptors that effect cell activation events mediated by serine proteases generated during inflammatory, fibrinolytic or haemostatic-regulated pathways. To further characterize the molecular genetics of these receptors, we have refined the genetic and physical mapping of both PAR-2 and TR. Utilization of two distinct radiation hybrid mapping panels with different levels of resolution demonstrated that both genes are tightly linked to the microsatellite markers D5S424, D5S1977, D5S2529 and D5S2596 (in order of decreasing LOD scores, from 13.7 for D5S424 to 7.7 for D5S2596). Physical mapping using yeast artificial chromosomes (YACs) and inversion field gel electrophoresis demonstrated that they are maximally separate by 90 kb. If the association of TR and PAR-2 genes resulted from a relatively recent gene duplication event from a common ancestral gene, these observations provide a general framework for the identification of gene transcripts representing alternative proteolytically activated receptors which may be clustered within this region of the human genome. These observations are especially relevant given recent evidence that murine and human platelets express alternative signalling mechanisms or receptors for thrombin.