Bacterial DnaK reduces the activity of anti-cancer drugs cisplatin and 5FU

BACKGROUND

Chemotherapy is a primary treatment for cancer, but its efficacy is often limited by cancer-associated bacteria (CAB) that impair tumor suppressor functions. Our previous research found that Mycoplasma fermentans DnaK, a chaperone protein, impairs p53 activities, which are essential for most anti-cancer chemotherapeutic responses.

METHODS

To investigate the role of DnaK in chemotherapy, we treated cancer cell lines with M. fermentans DnaK and then with commonly used p53-dependent anti-cancer drugs (cisplatin and 5FU). We evaluated the cells' survival in the presence or absence of a DnaK-binding peptide (ARV-1502). We also validated our findings using primary tumor cells from a novel DnaK knock-in mouse model. To provide a broader context for the clinical significance of these findings, we investigated human primary cancer sequencing datasets from The Cancer Genome Atlas (TCGA). We identified F. nucleatum as a CAB carrying DnaK with an amino acid composition highly similar to M. fermentans DnaK. Therefore, we investigated the effect of F. nucleatum DnaK on the anti-cancer activity of cisplatin and 5FU.

RESULTS

Our results show that both M. fermentans and F. nucleatum DnaKs reduce the effectiveness of cisplatin and 5FU. However, the use of ARV-1502 effectively restored the drugs' anti-cancer efficacy.

CONCLUSIONS

Our findings offer a practical framework for designing and implementing novel personalized anti-cancer strategies by targeting specific bacterial DnaKs in patients with poor response to chemotherapy, underscoring the potential for microbiome-based personalized cancer therapies.

Direct sequencing of insect symbionts via nanopore adaptive sampling

Insect symbionts can alter their host phenotype and their effects can range from beneficial to pathogenic. Moreover, many insects exhibit co-infections, making their study more challenging. Less than 1% of insect species have high-quality referenced genomes available and fewer still also have their symbionts sequenced. Two methods are commonly used to sequence symbionts: whole-genome sequencing to concomitantly capture the host and bacterial genomes, or isolation of the symbiont's genome before sequencing. These methods are limited when dealing with rare or poorly characterized symbionts. Long-read technology is an important tool to generate high-quality genomes as they can overcome high levels of heterozygosity, repeat content, and transposable elements that confound short-read methods. Oxford Nanopore (ONT) adaptive sampling allows a sequencing instrument to select or reject sequences in real time. We describe a method based on ONT adaptive sampling (subtractive) approach that readily permitted the sequencing of the complete genomes of mitochondria, Buchnera and its plasmids (pLeu, pTrp), and Wolbachia genomes in two aphid species, Aphis glycines and Pentalonia nigronervosa. Adaptive sampling is able to retrieve organelles such as mitochondria and symbionts that have high representation in their hosts such as Buchnera and Wolbachia, but is less successful at retrieving symbionts in low concentrations.

Comparative genomics of the niche-specific plant pathogen Streptomyces ipomoeae reveal novel genome content and organization

IMPORTANCE

While most plant-pathogenic Streptomyces species cause scab disease on a variety of plant hosts, Streptomyces ipomoeae is the sole causative agent of soil rot disease of sweet potato and closely related plant species. Here, genome sequencing of virulent and avirulent S. ipomoeae strains coupled with comparative genomic analyses has identified genome content and organization features unique to this streptomycete plant pathogen. The results here will enable future research into the mechanisms used by S. ipomoeae to cause disease and to persist in its niche environment.

Intestinal microbiota signatures of clinical response and immune-related adverse events in melanoma patients treated with anti-PD-1

Ample evidence indicates that the gut microbiome is a tumor-extrinsic factor associated with antitumor response to anti-programmed cell death protein-1 (PD-1) therapy, but inconsistencies exist between published microbial signatures associated with clinical outcomes. To resolve this, we evaluated a new melanoma cohort, along with four published datasets. Time-to-event analysis showed that baseline microbiota composition was optimally associated with clinical outcome at approximately 1 year after initiation of treatment. Meta-analysis and other bioinformatic analyses of the combined data show that bacteria associated with favorable response are confined within the Actinobacteria phylum and the Lachnospiraceae/Ruminococcaceae families of Firmicutes. Conversely, Gram-negative bacteria were associated with an inflammatory host intestinal gene signature, increased blood neutrophil-to-lymphocyte ratio, and unfavorable outcome. Two microbial signatures, enriched for Lachnospiraceae spp. and Streptococcaceae spp., were associated with favorable and unfavorable clinical response, respectively, and with distinct immune-related adverse effects. Despite between-cohort heterogeneity, optimized all-minus-one supervised learning algorithms trained on batch-corrected microbiome data consistently predicted outcomes to programmed cell death protein-1 therapy in all cohorts. Gut microbial communities (microbiotypes) with nonuniform geographical distribution were associated with favorable and unfavorable outcomes, contributing to discrepancies between cohorts. Our findings shed new light on the complex interaction between the gut microbiome and response to cancer immunotherapy, providing a roadmap for future studies.

Dietary fiber and probiotics influence the gut microbiome and melanoma immunotherapy response

Gut bacteria modulate the response to immune checkpoint blockade (ICB) treatment in cancer, but the effect of diet and supplements on this interaction is not well studied. We assessed fecal microbiota profiles, dietary habits, and commercially available probiotic supplement use in melanoma patients and performed parallel preclinical studies. Higher dietary fiber was associated with significantly improved progression-free survival in 128 patients on ICB, with the most pronounced benefit observed in patients with sufficient dietary fiber intake and no probiotic use. Findings were recapitulated in preclinical models, which demonstrated impaired treatment response to anti–programmed cell death 1 (anti–PD-1)–based therapy in mice receiving a low-fiber diet or probiotics, with a lower frequency of interferon-γ–positive cytotoxic T cells in the tumor microenvironment. Together, these data have clinical implications for patients receiving ICB for cancer.

Neonatal exposure to a wild-derived microbiome protects mice against diet-induced obesity

Obesity and its consequences are among the greatest challenges in healthcare. The gut microbiome is recognized as a key factor in the pathogenesis of obesity. Using a mouse model, we show here that a wild-derived microbiome protects against excessive weight gain, severe fatty liver disease and metabolic syndrome during a 10-week course of high-fat diet. This phenotype is transferable only during the first weeks of life. In adult mice, neither transfer nor severe disturbance of the wild-type microbiome modifies the metabolic response to a high-fat diet. The protective phenotype is associated with increased secretion of metabolic hormones and increased energy expenditure through activation of brown adipose tissue. Thus, we identify a microbiome that protects against weight gain and its negative consequences through metabolic programming in early life. Translation of these results to humans may identify early-life therapeutics that protect against obesity.

Abstract LB062: Efficacy of Responder-derived Fecal Microbiota Transplant (R-FMT) and Pembrolizumab in Anti-PD-1 Refractory Patients with Advanced Melanoma

Background: Monoclonal antibodies (mAb) targeting the programmed cell death protein 1 (PD-1) receptor provide durable long-term benefit in a subset of patients (pts) with advanced melanoma with response rates of 35-42% and 4-year progression-free survival (PFS) rate of 27%. Separately, the composition of the gut microbiota has been shown to correlate with anti-PD-1 efficacy in human cancer pts with melanoma, renal cell cancer and non-small cell lung cancer (NSCLC) although the precise organisms differ considerably across various studies. In preclinical models, responder-derived fecal microbiome and microbiome consortia produce anti-tumor responses. The effect of microbiome modulation in pts with anti-PD-1 refractory melanoma has not been evaluated. Methods: To evaluate whether primary resistance to anti-PD-1 immunotherapy could be overcome by intestinal microbiome modulation, we designed and conducted a phase II study (NCT03341143). We enrolled pts with primary refractory metastatic melanoma with best response of short-term stable disease (≤6 months) or progressive disease (PD) to prior anti-PD-1 based immunotherapy. Pts received single-administration of responder-derived fecal microbiota transplantation (R-FMT) together with pembrolizumab. Candidate donors were pts with advanced melanoma treated with anti-PD-1 immunotherapy with durable partial or complete response (PR, CR). Pembrolizumab was continued till intolerable toxicity or disease progression. Safety and clinical activity (based on RECIST v1.1) were main objectives; while progression-free survival (PFS) was a key secondary endpoint. Results: As of December 1, 2020, 16 pts with primary refractory melanoma were enrolled, of whom 15 were evaluable. LDH was elevated in 14/15 pts; and the median number of prior therapies was 2. Recipient pts were seromatched to receive a single R-FMT from one of eight candidate donors (5 CR; 3 PR; median PFS 58 months, range 43-70). R-FMT was administered via colonoscopy after bowel preparation with no use of antibiotics. Pembrolizumab was administered IV per label. R-FMT/pembrolizumab was well-tolerated, with no unusual toxicity signals. R-FMT induced rapid and durable microbiota perturbation in most pts; while 6 of 15 evaluable pts had evidence of clinical benefit. Response to R-FMT/pembrolizumab was associated with an increased abundance of taxa previously shown to be associated with response to anti-PD-1, increased CD8+ T cell activation, and decreased frequency of IL-8 expressing myeloid cells. Responders had distinct proteomic and metabolomic signatures, and transkingdom network analyses confirmed that the gut microbiome regulated these changes. Conclusions: In pts with anti-PD-1 primary refractory melanoma, R-FMT/pembrolizumab changed the gut microbiome and reprogrammed the tumor microenvironment to overcome resistance to anti-PD-1 immunotherapy. Response was associated with CD8 T cell induction and reduction of IL-8 expressing myeloid cells.

Citation Format: Diwakar Davar, Amiran Dzutsev, John A. McCulloch, Richard R. Rodrigues, Joe-Marc Chauvin, Robert M. Morrison, Richelle N. Deblasio, Carmine Menna, Quanquan Ding, Ornella Pagliano, Bochra Zidi, Shuowen Zhang, Jonathan H. Badger, Marie Vetizou, Alicia M. Cole, Miriam R. Fernandes, Stephanie Prescott, Raquel G. Costa, Ascharya K. Balaji, Andrey Morgun, Ivan Vujkovic-Cvijin, Hong Wang, Amir A. Borhani, Marc B. Schwartz, Howard M. Dubner, Scarlett J. Ernst, Amy Rose, Yana G. Najjar, Yasmine Belkaid, John M. Kirkwood, Giorgio Trinchieri, Hassane M. Zarour. Efficacy of Responder-derived Fecal Microbiota Transplant (R-FMT) and Pembrolizumab in Anti-PD-1 Refractory Patients with Advanced Melanoma [abstract]. In: Proceedings of the American Association for Cancer Research Annual Meeting 2021; 2021 Apr 10-15 and May 17-21. Philadelphia (PA): AACR; Cancer Res 2021;81(13_Suppl):Abstract nr LB062.

Fecal microbiota transplant overcomes resistance to anti-PD-1 therapy in melanoma patients

Anti-programmed cell death protein 1 (PD-1) therapy provides long-term clinical benefits to patients with advanced melanoma. The composition of the gut microbiota correlates with anti-PD-1 efficacy in preclinical models and cancer patients. To investigate whether resistance to anti-PD-1 can be overcome by changing the gut microbiota, this clinical trial evaluated the safety and efficacy of responder-derived fecal microbiota transplantation (FMT) together with anti-PD-1 in patients with PD-1-refractory melanoma. This combination was well tolerated, provided clinical benefit in 6 of 15 patients, and induced rapid and durable microbiota perturbation. Responders exhibited increased abundance of taxa that were previously shown to be associated with response to anti-PD-1, increased CD8⁺ T cell activation, and decreased frequency of interleukin-8-expressing myeloid cells. Responders had distinct proteomic and metabolomic signatures, and transkingdom network analyses confirmed that the gut microbiome regulated these changes. Collectively, our findings show that FMT and anti-PD-1 changed the gut microbiome and reprogrammed the tumor microenvironment to overcome resistance to anti-PD-1 in a subset of PD-1 advanced melanoma.

Soybean aphid biotype 1 genome: Insights into the invasive biology and adaptive evolution of a major agricultural pest

The soybean aphid, Aphis glycines Matsumura (Hemiptera: Aphididae) is a serious pest of the soybean plant, Glycine max, a major world-wide agricultural crop. We assembled a de novo genome sequence of Ap. glycines Biotype 1, from a culture established shortly after this species invaded North America. 20.4% of the Ap. glycines proteome is duplicated. These in-paralogs are enriched with Gene Ontology (GO) categories mostly related to apoptosis, a possible adaptation to plant chemistry and other environmental stressors. Approximately one-third of these genes show parallel duplication in other aphids. But Ap. gossypii, its closest related species, has the lowest number of these duplicated genes. An Illumina GoldenGate assay of 2380 SNPs was used to determine the world-wide population structure of Ap. Glycines. China and South Korean aphids are the closest to those in North America. China is the likely origin of other Asian aphid populations. The most distantly related aphids to those in North America are from Australia. The diversity of Ap. glycines in North America has decreased over time since its arrival. The genetic diversity of Ap. glycines North American population sampled shortly after its first detection in 2001 up to 2012 does not appear to correlate with geography. However, aphids collected on soybean Rag experimental varieties in Minnesota (MN), Iowa (IA), and Wisconsin (WI), closer to high density Rhamnus cathartica stands, appear to have higher capacity to colonize resistant soybean plants than aphids sampled in Ohio (OH), North Dakota (ND), and South Dakota (SD). Samples from the former states have SNP alleles with high F_ST values and frequencies, that overlap with genes involved in iron metabolism, a crucial metabolic pathway that may be affected by the Rag-associated soybean plant response. The Ap. glycines Biotype 1 genome will provide needed information for future analyses of mechanisms of aphid virulence and pesticide resistance as well as facilitate comparative analyses between aphids with differing natural history and host plant range.

Abstract B17: IKKα/STAT3 antagonistic signaling regulates fungi-bacteria endosymbiosis-associated carcinogenesis

Clinically, oral fungating lesions frequently occur in cancer patients following radiation and chemotherapy. However, their role in carcinogenesis has not been fully understood. Our previous finding showed that defective IKKα-mediated autoreactive T cells promote fungal infection and esophageal malignancy, but whether IKKα-mediated altered epithelial properties regulate fungal colonization and carcinogenesis remains to be tested. To investigate whether IKKα-mediated altered epithelial properties regulate fungal colonization and carcinogenesis, we orally infected mice lacking Ikkα (IkkαΔSOE) or both Ikkα and Stat3 (IkkαΔSOEStat3ΔSOE) in the stem cells of the skin and mouth with fungi. We found that oral fungal infection enhanced Ikkα ablation-mediated skin and oral carcinogenesis in association with increased oral and skin bacteria counts and altered oral bacteria phylum from gram-positive to gram-negative. Enhanced oral fungal/bacterial infection increased inflammasome activities and systemic IL-17A/IL-1β signals. Surprisingly, Stat3 ablation rescued these phenotypes in IkkαΔSOEStat3ΔSOE mice. In addition, amoxicillin treatment inhibited skin carcinogenesis and oral/skin bacteria counts and reversed oral bacteria phyla in fungal-inoculated IkkαΔSOE mice. These findings highlight that cancer cell IKKα/STAT3 antagonizing pathways determine susceptibility to fungal and bacterial colonization, which impacts carcinogenesis. Human head and neck squamous cell carcinomas also share oncogenic pathways and bacteria sequences/counts with mice, suggesting that oral fungating/ulcerating lesions enhance local and distal tumorigenesis in humans.

Citation Format: Na-Young Song, Xin Li, Jonathan H. Badger, Jami Willette Brown, Xhonghe Sun, Gongping Shi, Feng Zhu, Chengfei Jiang, Colum O'hUigin, Xiaolin Wu, Giorgio Trinchieri, Yinling Hu. IKKα/STAT3 antagonistic signaling regulates fungi-bacteria endosymbiosis-associated carcinogenesis [abstract]. In: Proceedings of the AACR Special Conference on the Microbiome, Viruses, and Cancer; 2020 Feb 21-24; Orlando, FL. Philadelphia (PA): AACR; Cancer Res 2020;80(8 Suppl):Abstract nr B17.

Microbiome as an Immunological Modifier

Humans are living ecosystems composed of human cells and microbes. The microbiome is the collection of microbes (microbiota) and their genes. Recent breakthroughs in the high-throughput sequencing technologies have made it possible for us to understand the composition of the human microbiome. Launched by the National Institutes of Health in USA, the human microbiome project indicated that our bodies harbor a wide array of microbes, specific to each body site with interpersonal and intrapersonal variabilities. Numerous studies have indicated that several factors influence the development of the microbiome including genetics, diet, use of antibiotics, and lifestyle, among others. The microbiome and its mediators are in a continuous cross talk with the host immune system; hence, any imbalance on one side is reflected on the other. Dysbiosis (microbiota imbalance) was shown in many diseases and pathological conditions such as inflammatory bowel disease, celiac disease, multiple sclerosis, rheumatoid arthritis, asthma, diabetes, and cancer. The microbial composition mirrors inflammation variations in certain disease conditions, within various stages of the same disease; hence, it has the potential to be used as a biomarker.

Description of soybean aphid (Aphis glycines Matsumura) mitochondrial genome and comparative mitogenomics of Aphididae (Hemiptera: Sternorrhyncha)

The complete mitochondrial genome of the soybean aphid (Aphis glycines Matsumura), a major agricultural pest in the world, is described for the first time, which consists of 13 protein-coding genes, 22 tRNA genes, 2 rRNA genes, as well as a large repeat region between tRNA-Glu and tRNA-Phe, and an AT-rich control region. The 17,954 bp mtgenome is the largest reported from the family Aphididae, and its gene order follows the ancestral insect mtgenome except for the repeat region, which contains a 195 bp unit repeated 11.9 times, representing the highest reported repeats among the known aphid mtgenomes to date. A new molecular phylogeny of Aphidae is reconstructed based on all available aphid mtgenomes, and it is shown that the mtgenome data can robustly resolve relationships at the subfamily level, but do not have sufficient phylogenetic information to resolve deep relationships. A phylogeny-based comparative analysis of mtgenomes has been performed to investigate the evolution of the repeat region between tRNA-Glu and tRNA-Phe. So far, among aphids, 13 species are known to have this repeat region of variable lengths, and a phylogenetic analysis of the repeat region shows that a large proportion of the sequences are conserved across the phylogeny, suggesting that the repeat region evolved in the most recent common ancestor of Aphidinae and Eriosomatinae, and that it has gone through numerous episodes of lineage-specific losses and expansions. Combined together, this study provides novel insights into how the repeat regions have evolved within aphids.

Laboratory mice born to wild mice have natural microbiota and model human immune responses

Laboratory mouse studies are paramount for understanding basic biological phenomena but also have limitations. These include conflicting results caused by divergent microbiota and limited translational research value. To address both shortcomings, we transferred C57BL/6 embryos into wild mice, creating "wildlings." These mice have a natural microbiota and pathogens at all body sites and the tractable genetics of C57BL/6 mice. The bacterial microbiome, mycobiome, and virome of wildlings affect the immune landscape of multiple organs. Their gut microbiota outcompete laboratory microbiota and demonstrate resilience to environmental challenges. Wildlings, but not conventional laboratory mice, phenocopied human immune responses in two preclinical studies. A combined natural microbiota- and pathogen-based model may enhance the reproducibility of biomedical studies and increase the bench-to-bedside safety and success of immunological studies.

An Interleukin-23-Interleukin-22 Axis Regulates Intestinal Microbial Homeostasis to Protect from Diet-Induced Atherosclerosis

Although commensal flora is involved in the regulation of immunity, the interplay between cytokine signaling and microbiota in atherosclerosis remains unknown. We found that interleukin (IL)-23 and its downstream target IL-22 restricted atherosclerosis by repressing pro-atherogenic microbiota. Inactivation of IL-23-IL-22 signaling led to deterioration of the intestinal barrier, dysbiosis, and expansion of pathogenic bacteria with distinct biosynthetic and metabolic properties, causing systemic increase in pro-atherogenic metabolites such as lipopolysaccharide (LPS) and trimethylamine N-oxide (TMAO). Augmented disease in the absence of the IL-23-IL-22 pathway was mediated in part by pro-atherogenic osteopontin, controlled by microbial metabolites. Microbiota transfer from IL-23-deficient mice accelerated atherosclerosis, whereas microbial depletion or IL-22 supplementation reduced inflammation and ameliorated disease. Our work uncovers the IL-23-IL-22 signaling as a regulator of atherosclerosis that restrains expansion of pro-atherogenic microbiota and argues for informed use of cytokine blockers to avoid cardiovascular side effects driven by microbiota and inflammation.

The Antimicrobial Peptide CRAMP Is Essential for Colon Homeostasis by Maintaining Microbiota Balance

Commensal bacteria are critical for physiological functions in the gut, and dysbiosis in the gut may cause diseases. In this article, we report that mice deficient in cathelin-related antimicrobial peptide (CRAMP) were defective in the development of colon mucosa and highly sensitive to dextran sulfate sodium (DSS)-elicited colitis, as well as azoxymethane-mediated carcinogenesis. Pretreatment of CRAMP^-/- mice with antibiotics markedly reduced the severity of DSS-induced colitis, suggesting CRAMP as a limiting factor on dysbiosis in the colon. This was supported by observations that wild-type (WT) mice cohoused with CRAMP^-/- mice became highly sensitive to DSS-induced colitis, and the composition of fecal microbiota was skewed by CRAMP deficiency. In particular, several bacterial species that are typically found in oral microbiota, such as Mogibacterium neglectum, Desulfovibrio piger, and Desulfomicrobium orale, were increased in feces of CRAMP^-/- mice and were transferred to WT mice during cohousing. When littermates of CRAMP^+/- parents were examined, the composition of the fecal microbiota of WT pups and heterozygous parents was similar. In contrast, although the difference in fecal microbiota between CRAMP^-/- and WT pups was small early on after weaning and single mouse housing, there was an increasing divergence with prolonged single housing. These results indicate that CRAMP is critical in maintaining colon microbiota balance and supports mucosal homeostasis, anti-inflammatory responses, and protection from carcinogenesis.

Wild Mouse Gut Microbiota Promotes Host Fitness and Improves Disease Resistance

Laboratory mice, while paramount for understanding basic biological phenomena, are limited in modeling complex diseases of humans and other free-living mammals. Because the microbiome is a major factor in mammalian physiology, we aimed to identify a naturally evolved reference microbiome to better recapitulate physiological phenomena relevant in the natural world outside the laboratory. Among 21 distinct mouse populations worldwide, we identified a closely related wild relative to standard laboratory mouse strains. Its bacterial gut microbiome differed significantly from its laboratory mouse counterpart and was transferred to and maintained in laboratory mice over several generations. Laboratory mice reconstituted with natural microbiota exhibited reduced inflammation and increased survival following influenza virus infection and improved resistance against mutagen/inflammation-induced colorectal tumorigenesis. By demonstrating the host fitness-promoting traits of natural microbiota, our findings should enable the discovery of protective mechanisms relevant in the natural world and improve the modeling of complex diseases of free-living mammals. VIDEO ABSTRACT.

Correction: Transcriptional Orchestration of the Global Cellular Response of a Model Pennate Diatom to Diel Light Cycling under Iron Limitation

[This corrects the article DOI: 10.1371/journal.pgen.1006490.].

Microbes and Cancer

Commensal microorganisms (the microbiota) live on all the surface barriers of our body and are particularly abundant and diverse in the distal gut. The microbiota and its larger host represent a metaorganism in which the cross talk between microbes and host cells is necessary for health, survival, and regulation of physiological functions locally, at the barrier level, and systemically. The ancestral molecular and cellular mechanisms stemming from the earliest interactions between prokaryotes and eukaryotes have evolved to mediate microbe-dependent host physiology and tissue homeostasis, including innate and adaptive resistance to infections and tissue repair. Mostly because of its effects on metabolism, cellular proliferation, inflammation, and immunity, the microbiota regulates cancer at the level of predisposing conditions, initiation, genetic instability, susceptibility to host immune response, progression, comorbidity, and response to therapy. Here, we review the mechanisms underlying the interaction of the microbiota with cancer and the evidence suggesting that the microbiota could be targeted to improve therapy while attenuating adverse reactions.

Transcriptional Orchestration of the Global Cellular Response of a Model Pennate Diatom to Diel Light Cycling under Iron Limitation

Environmental fluctuations affect distribution, growth and abundance of diatoms in nature, with iron (Fe) availability playing a central role. Studies on the response of diatoms to low Fe have either utilized continuous (24 hr) illumination or sampled a single time of day, missing any temporal dynamics. We profiled the physiology, metabolite composition, and global transcripts of the pennate diatom Phaeodactylum tricornutum during steady-state growth at low, intermediate, and high levels of dissolved Fe over light:dark cycles, to better understand fundamental aspects of genetic control of physiological acclimation to growth under Fe-limitation. We greatly expand the catalog of genes involved in the low Fe response, highlighting the importance of intracellular trafficking in Fe-limited diatoms. P. tricornutum exhibited transcriptomic hallmarks of slowed growth leading to prolonged periods of cell division/silica deposition, which could impact biogeochemical carbon sequestration in Fe-limited regions. Light harvesting and ribosome biogenesis transcripts were generally reduced under low Fe while transcript levels for genes putatively involved in the acquisition and recycling of Fe were increased. We also noted shifts in expression towards increased synthesis and catabolism of branched chain amino acids in P. tricornutum grown at low Fe whereas expression of genes involved in central core metabolism were relatively unaffected, indicating that essential cellular function is protected. Beyond the response of P. tricornutum to low Fe, we observed major coordinated shifts in transcript control of primary and intermediate metabolism over light:dark cycles which contribute to a new view of the significance of distinctive diatom pathways, such as mitochondrial glycolysis and the ornithine-urea cycle. This study provides new insight into transcriptional modulation of diatom physiology and metabolism across light:dark cycles in response to Fe availability, providing mechanistic understanding for the ability of diatoms to remain metabolically poised to respond quickly to Fe input and revealing strategies underlying their ecological success.

Oceanic diatoms live in constantly fluctuating environments to which they must adapt in order to survive. During sunlit hours, photosynthesis occurs allowing diatoms to store energy used at night to sustain energy demands. Cellular and molecular mechanisms for regulation of phytoplankton growth are important to understand because of their environmental roles at the base of food webs and in regulating carbon flux out of the atmosphere. In ocean ecosystems, the availability of iron (Fe) commonly limits phytoplankton growth and diatoms typically outcompete other phytoplankton when Fe is added, indicating they have adaptations allowing them to both survive at low Fe and rapidly respond to Fe additions. These adaptations may be unique depending on isolation from coastal or oceanic locations. To identify adaptive strategies, we characterized the response of a model diatom, Phaeodactylum tricornutum, to limiting Fe conditions over day:night cycles using a combination of gene expression analyses, metabolite, and physiology measurements. Major coordinated shifts in metabolism and growth were documented over diel cycles, with peak expression of low Fe expressed genes in the dark phase. Diatoms respond to limiting Fe by increasing Fe acquisition, while decreasing growth rate through slowed cell cycle progression, reduced energy acquisition, and subtle metabolic remodeling.

Structure of the germline genome of Tetrahymena thermophila and relationship to the massively rearranged somatic genome

The germline genome of the binucleated ciliate Tetrahymena thermophila undergoes programmed chromosome breakage and massive DNA elimination to generate the somatic genome. Here, we present a complete sequence assembly of the germline genome and analyze multiple features of its structure and its relationship to the somatic genome, shedding light on the mechanisms of genome rearrangement as well as the evolutionary history of this remarkable germline/soma differentiation. Our results strengthen the notion that a complex, dynamic, and ongoing interplay between mobile DNA elements and the host genome have shaped Tetrahymena chromosome structure, locally and globally. Non-standard outcomes of rearrangement events, including the generation of short-lived somatic chromosomes and excision of DNA interrupting protein-coding regions, may represent novel forms of developmental gene regulation. We also compare Tetrahymena's germline/soma differentiation to that of other characterized ciliates, illustrating the wide diversity of adaptations that have occurred within this phylum.

Draft genome sequences for the obligate bacterial predators Bacteriovorax spp. of four phylogenetic clusters

Bacteriovorax is the halophilic genus of the obligate bacterial predators, Bdellovibrio and like organisms. The predators are known for their unique biphasic life style in which they search for and attack their prey in the free living phase; penetrate, grow, multiply and lyse the prey in the intraperiplasmic phase. Bacteriovorax isolates representing four phylogenetic clusters were selected for genomic sequencing. Only one type strain genome has been published so far from the genus Bacteriovorax. We report the genomes from non-type strains isolated from aquatic environments. Here we describe and compare the genomic features of the four strains, together with the classification and annotation.

Functional tradeoffs underpin salinity-driven divergence in microbial community composition

Bacterial community composition and functional potential change subtly across gradients in the surface ocean. In contrast, while there are significant phylogenetic divergences between communities from freshwater and marine habitats, the underlying mechanisms to this phylogenetic structuring yet remain unknown. We hypothesized that the functional potential of natural bacterial communities is linked to this striking divide between microbiomes. To test this hypothesis, metagenomic sequencing of microbial communities along a 1,800 km transect in the Baltic Sea area, encompassing a continuous natural salinity gradient from limnic to fully marine conditions, was explored. Multivariate statistical analyses showed that salinity is the main determinant of dramatic changes in microbial community composition, but also of large scale changes in core metabolic functions of bacteria. Strikingly, genetically and metabolically different pathways for key metabolic processes, such as respiration, biosynthesis of quinones and isoprenoids, glycolysis and osmolyte transport, were differentially abundant at high and low salinities. These shifts in functional capacities were observed at multiple taxonomic levels and within dominant bacterial phyla, while bacteria, such as SAR11, were able to adapt to the entire salinity gradient. We propose that the large differences in central metabolism required at high and low salinities dictate the striking divide between freshwater and marine microbiomes, and that the ability to inhabit different salinity regimes evolved early during bacterial phylogenetic differentiation. These findings significantly advance our understanding of microbial distributions and stress the need to incorporate salinity in future climate change models that predict increased levels of precipitation and a reduction in salinity.

Candidate phylum TM6 genome recovered from a hospital sink biofilm provides genomic insights into this uncultivated phylum

The "dark matter of life" describes microbes and even entire divisions of bacterial phyla that have evaded cultivation and have yet to be sequenced. We present a genome from the globally distributed but elusive candidate phylum TM6 and uncover its metabolic potential. TM6 was detected in a biofilm from a sink drain within a hospital restroom by analyzing cells using a highly automated single-cell genomics platform. We developed an approach for increasing throughput and effectively improving the likelihood of sampling rare events based on forming small random pools of single-flow-sorted cells, amplifying their DNA by multiple displacement amplification and sequencing all cells in the pool, creating a "mini-metagenome." A recently developed single-cell assembler, SPAdes, in combination with contig binning methods, allowed the reconstruction of genomes from these mini-metagenomes. A total of 1.07 Mb was recovered in seven contigs for this member of TM6 (JCVI TM6SC1), estimated to represent 90% of its genome. High nucleotide identity between a total of three TM6 genome drafts generated from pools that were independently captured, amplified, and assembled provided strong confirmation of a correct genomic sequence. TM6 is likely a Gram-negative organism and possibly a symbiont of an unknown host (nonfree living) in part based on its small genome, low-GC content, and lack of biosynthesis pathways for most amino acids and vitamins. Phylogenomic analysis of conserved single-copy genes confirms that TM6SC1 is a deeply branching phylum.

Transcriptomic analysis of metabolic function in the giant kelp, Macrocystis pyrifera, across depth and season

To increase knowledge of transcript diversity for the giant kelp, Macrocystis pyrifera, and assess gene expression across naturally occurring depth gradients in light, temperature and nutrients, we sequenced four cDNA libraries created from blades collected at the sea surface and at 18 m depth during the winter and summer. Comparative genomics cluster analyses revealed novel gene families (clusters) in existing brown alga expressed sequence tag data compared with other related algal groups, a pattern also seen with the addition of M. pyrifera sequences. Assembly of 228 Mbp of sequence generated c. 9000 isotigs and c. 12,000 open reading frames. Annotations were assigned using families of hidden Markov models for c. 11% of open reading frames; M. pyrifera had highest similarity to other members of the Phaeophyceae, namely Ectocarpus siliculosus and Laminaria digitata. Quantitative polymerase chain reaction of transcript targets verified depth-related differences in gene expression; stress response and light-harvesting transcripts, especially members of the LI818 (also known as LHCSR) family, showed high expression in the surface compared with 18 m depth, while some nitrogen acquisition transcripts (e.g. nitrite reductase) were upregulated at depth compared with the surface, supporting a conceptual biological model of depth-dependent physiology.

The "most wanted" taxa from the human microbiome for whole genome sequencing

The goal of the Human Microbiome Project (HMP) is to generate a comprehensive catalog of human-associated microorganisms including reference genomes representing the most common species. Toward this goal, the HMP has characterized the microbial communities at 18 body habitats in a cohort of over 200 healthy volunteers using 16S rRNA gene (16S) sequencing and has generated nearly 1,000 reference genomes from human-associated microorganisms. To determine how well current reference genome collections capture the diversity observed among the healthy microbiome and to guide isolation and future sequencing of microbiome members, we compared the HMP's 16S data sets to several reference 16S collections to create a 'most wanted' list of taxa for sequencing. Our analysis revealed that the diversity of commonly occurring taxa within the HMP cohort microbiome is relatively modest, few novel taxa are represented by these OTUs and many common taxa among HMP volunteers recur across different populations of healthy humans. Taken together, these results suggest that it should be possible to perform whole-genome sequencing on a large fraction of the human microbiome, including the 'most wanted', and that these sequences should serve to support microbiome studies across multiple cohorts. Also, in stark contrast to other taxa, the 'most wanted' organisms are poorly represented among culture collections suggesting that novel culture- and single-cell-based methods will be required to isolate these organisms for sequencing.

Comparative genomics of the pathogenic ciliate Ichthyophthirius multifiliis, its free-living relatives and a host species provide insights into adoption of a parasitic lifestyle and prospects for disease control

Background

Ichthyophthirius multifiliis, commonly known as Ich, is a highly pathogenic ciliate responsible for 'white spot', a disease causing significant economic losses to the global aquaculture industry. Options for disease control are extremely limited, and Ich's obligate parasitic lifestyle makes experimental studies challenging. Unlike most well-studied protozoan parasites, Ich belongs to a phylum composed primarily of free-living members. Indeed, it is closely related to the model organism Tetrahymena thermophila. Genomic studies represent a promising strategy to reduce the impact of this disease and to understand the evolutionary transition to parasitism.

Results

We report the sequencing, assembly and annotation of the Ich macronuclear genome. Compared with its free-living relative T. thermophila, the Ich genome is reduced approximately two-fold in length and gene density and three-fold in gene content. We analyzed in detail several gene classes with diverse functions in behavior, cellular function and host immunogenicity, including protein kinases, membrane transporters, proteases, surface antigens and cytoskeletal components and regulators. We also mapped by orthology Ich's metabolic pathways in comparison with other ciliates and a potential host organism, the zebrafish Danio rerio.

Conclusions

Knowledge of the complete protein-coding and metabolic potential of Ich opens avenues for rational testing of therapeutic drugs that target functions essential to this parasite but not to its fish hosts. Also, a catalog of surface protein-encoding genes will facilitate development of more effective vaccines. The potential to use T. thermophila as a surrogate model offers promise toward controlling 'white spot' disease and understanding the adaptation to a parasitic lifestyle.

Efficient de novo assembly of single-cell bacterial genomes from short-read data sets

Whole genome amplification by the multiple displacement amplification (MDA) method allows sequencing of DNA from single cells of bacteria that cannot be cultured. Assembling a genome is challenging, however, because MDA generates highly nonuniform coverage of the genome. Here we describe an algorithm tailored for short-read data from single cells that improves assembly through the use of a progressively increasing coverage cutoff. Assembly of reads from single Escherichia coli and Staphylococcus aureus cells captures >91% of genes within contigs, approaching the 95% captured from an assembly based on many E. coli cells. We apply this method to assemble a genome from a single cell of an uncultivated SAR324 clade of Deltaproteobacteria, a cosmopolitan bacterial lineage in the global ocean. Metabolic reconstruction suggests that SAR324 is aerobic, motile and chemotaxic. Our approach enables acquisition of genome assemblies for individual uncultivated bacteria using only short reads, providing cell-specific genetic information absent from metagenomic studies.

Comparative genomic analysis of Vibrio parahaemolyticus: serotype conversion and virulence

BACKGROUND

Vibrio parahaemolyticus is a common cause of foodborne disease. Beginning in 1996, a more virulent strain having serotype O3:K6 caused major outbreaks in India and other parts of the world, resulting in the emergence of a pandemic. Other serovariants of this strain emerged during its dissemination and together with the original O3:K6 were termed strains of the pandemic clone. Two genomes, one of this virulent strain and one pre-pandemic strain have been sequenced. We sequenced four additional genomes of V. parahaemolyticus in this study that were isolated from different geographical regions and time points. Comparative genomic analyses of six strains of V. parahaemolyticus isolated from Asia and Peru were performed in order to advance knowledge concerning the evolution of V. parahaemolyticus; specifically, the genetic changes contributing to serotype conversion and virulence. Two pre-pandemic strains and three pandemic strains, isolated from different geographical regions, were serotype O3:K6 and either toxin profiles (tdh+, trh-) or (tdh-, trh+). The sixth pandemic strain sequenced in this study was serotype O4:K68.

RESULTS

Genomic analyses revealed that the trh+ and tdh+ strains had different types of pathogenicity islands and mobile elements as well as major structural differences between the tdh pathogenicity islands of the pre-pandemic and pandemic strains. In addition, the results of single nucleotide polymorphism (SNP) analysis showed that 94% of the SNPs between O3:K6 and O4:K68 pandemic isolates were within a 141 kb region surrounding the O- and K-antigen-encoding gene clusters. The "core" genes of V. parahaemolyticus were also compared to those of V. cholerae and V. vulnificus, in order to delineate differences between these three pathogenic species. Approximately one-half (49-59%) of each species' core genes were conserved in all three species, and 14-24% of the core genes were species-specific and in different functional categories.

CONCLUSIONS

Our data support the idea that the pandemic strains are closely related and that recent South American outbreaks of foodborne disease caused by V. parahaemolyticus are closely linked to outbreaks in India. Serotype conversion from O3:K6 to O4:K68 was likely due to a recombination event involving a region much larger than the O-antigen- and K-antigen-encoding gene clusters. Major differences between pathogenicity islands and mobile elements are also likely driving the evolution of V. parahaemolyticus. In addition, our analyses categorized genes that may be useful in differentiating pathogenic Vibrios at the species level.

The Ectocarpus genome and the independent evolution of multicellularity in brown algae

Brown algae (Phaeophyceae) are complex photosynthetic organisms with a very different evolutionary history to green plants, to which they are only distantly related. These seaweeds are the dominant species in rocky coastal ecosystems and they exhibit many interesting adaptations to these, often harsh, environments. Brown algae are also one of only a small number of eukaryotic lineages that have evolved complex multicellularity (Fig. 1). We report the 214 million base pair (Mbp) genome sequence of the filamentous seaweed Ectocarpus siliculosus (Dillwyn) Lyngbye, a model organism for brown algae, closely related to the kelps (Fig. 1). Genome features such as the presence of an extended set of light-harvesting and pigment biosynthesis genes and new metabolic processes such as halide metabolism help explain the ability of this organism to cope with the highly variable tidal environment. The evolution of multicellularity in this lineage is correlated with the presence of a rich array of signal transduction genes. Of particular interest is the presence of a family of receptor kinases, as the independent evolution of related molecules has been linked with the emergence of multicellularity in both the animal and green plant lineages. The Ectocarpus genome sequence represents an important step towards developing this organism as a model species, providing the possibility to combine genomic and genetic approaches to explore these and other aspects of brown algal biology further.

Metagenomes from high-temperature chemotrophic systems reveal geochemical controls on microbial community structure and function

The Yellowstone caldera contains the most numerous and diverse geothermal systems on Earth, yielding an extensive array of unique high-temperature environments that host a variety of deeply-rooted and understudied Archaea, Bacteria and Eukarya. The combination of extreme temperature and chemical conditions encountered in geothermal environments often results in considerably less microbial diversity than other terrestrial habitats and offers a tremendous opportunity for studying the structure and function of indigenous microbial communities and for establishing linkages between putative metabolisms and element cycling. Metagenome sequence (14–15,000 Sanger reads per site) was obtained for five high-temperature (>65°C) chemotrophic microbial communities sampled from geothermal springs (or pools) in Yellowstone National Park (YNP) that exhibit a wide range in geochemistry including pH, dissolved sulfide, dissolved oxygen and ferrous iron. Metagenome data revealed significant differences in the predominant phyla associated with each of these geochemical environments. Novel members of the Sulfolobales are dominant in low pH environments, while other Crenarchaeota including distantly-related Thermoproteales and Desulfurococcales populations dominate in suboxic sulfidic sediments. Several novel archaeal groups are well represented in an acidic (pH 3) Fe-oxyhydroxide mat, where a higher O₂ influx is accompanied with an increase in archaeal diversity. The presence or absence of genes and pathways important in S oxidation-reduction, H₂-oxidation, and aerobic respiration (terminal oxidation) provide insight regarding the metabolic strategies of indigenous organisms present in geothermal systems. Multiple-pathway and protein-specific functional analysis of metagenome sequence data corroborated results from phylogenetic analyses and clearly demonstrate major differences in metabolic potential across sites. The distribution of functional genes involved in electron transport is consistent with the hypothesis that geochemical parameters (e.g., pH, sulfide, Fe, O₂) control microbial community structure and function in YNP geothermal springs.

The complete genome sequence of Haloferax volcanii DS2, a model archaeon

BACKGROUND

Haloferax volcanii is an easily culturable moderate halophile that grows on simple defined media, is readily transformable, and has a relatively stable genome. This, in combination with its biochemical and genetic tractability, has made Hfx. volcanii a key model organism, not only for the study of halophilicity, but also for archaeal biology in general.

METHODOLOGY/PRINCIPAL FINDINGS

We report here the sequencing and analysis of the genome of Hfx. volcanii DS2, the type strain of this species. The genome contains a main 2.848 Mb chromosome, three smaller chromosomes pHV1, 3, 4 (85, 438, 636 kb, respectively) and the pHV2 plasmid (6.4 kb).

CONCLUSIONS/SIGNIFICANCE

The completed genome sequence, presented here, provides an invaluable tool for further in vivo and in vitro studies of Hfx. volcanii.

The complete genome of Teredinibacter turnerae T7901: an intracellular endosymbiont of marine wood-boring bivalves (shipworms)

Here we report the complete genome sequence of Teredinibacter turnerae T7901. T. turnerae is a marine gamma proteobacterium that occurs as an intracellular endosymbiont in the gills of wood-boring marine bivalves of the family Teredinidae (shipworms). This species is the sole cultivated member of an endosymbiotic consortium thought to provide the host with enzymes, including cellulases and nitrogenase, critical for digestion of wood and supplementation of the host's nitrogen-deficient diet. T. turnerae is closely related to the free-living marine polysaccharide degrading bacterium Saccharophagus degradans str. 2–40 and to as yet uncultivated endosymbionts with which it coexists in shipworm cells. Like S. degradans, the T. turnerae genome encodes a large number of enzymes predicted to be involved in complex polysaccharide degradation (>100). However, unlike S. degradans, which degrades a broad spectrum (>10 classes) of complex plant, fungal and algal polysaccharides, T. turnerae primarily encodes enzymes associated with deconstruction of terrestrial woody plant material. Also unlike S. degradans and many other eubacteria, T. turnerae dedicates a large proportion of its genome to genes predicted to function in secondary metabolism. Despite its intracellular niche, the T. turnerae genome lacks many features associated with obligate intracellular existence (e.g. reduced genome size, reduced %G+C, loss of genes of core metabolism) and displays evidence of adaptations common to free-living bacteria (e.g. defense against bacteriophage infection). These results suggest that T. turnerae is likely a facultative intracellular ensosymbiont whose niche presently includes, or recently included, free-living existence. As such, the T. turnerae genome provides insights into the range of genomic adaptations associated with intracellular endosymbiosis as well as enzymatic mechanisms relevant to the recycling of plant materials in marine environments and the production of cellulose-derived biofuels.

Green evolution and dynamic adaptations revealed by genomes of the marine picoeukaryotes Micromonas

Picoeukaryotes are a taxonomically diverse group of organisms less than 2 micrometers in diameter. Photosynthetic marine picoeukaryotes in the genus Micromonas thrive in ecosystems ranging from tropical to polar and could serve as sentinel organisms for biogeochemical fluxes of modern oceans during climate change. These broadly distributed primary producers belong to an anciently diverged sister clade to land plants. Although Micromonas isolates have high 18S ribosomal RNA gene identity, we found that genomes from two isolates shared only 90% of their predicted genes. Their independent evolutionary paths were emphasized by distinct riboswitch arrangements as well as the discovery of intronic repeat elements in one isolate, and in metagenomic data, but not in other genomes. Divergence appears to have been facilitated by selection and acquisition processes that actively shape the repertoire of genes that are mutually exclusive between the two isolates differently than the core genes. Analyses of the Micromonas genomes offer valuable insights into ecological differentiation and the dynamic nature of early plant evolution.

Complete genome sequence of the aerobic CO-oxidizing thermophile Thermomicrobium roseum

In order to enrich the phylogenetic diversity represented in the available sequenced bacterial genomes and as part of an “Assembling the Tree of Life” project, we determined the genome sequence of Thermomicrobium roseum DSM 5159. T. roseum DSM 5159 is a red-pigmented, rod-shaped, Gram-negative extreme thermophile isolated from a hot spring that possesses both an atypical cell wall composition and an unusual cell membrane that is composed entirely of long-chain 1,2-diols. Its genome is composed of two circular DNA elements, one of 2,006,217 bp (referred to as the chromosome) and one of 919,596 bp (referred to as the megaplasmid). Strikingly, though few standard housekeeping genes are found on the megaplasmid, it does encode a complete system for chemotaxis including both chemosensory components and an entire flagellar apparatus. This is the first known example of a complete flagellar system being encoded on a plasmid and suggests a straightforward means for lateral transfer of flagellum-based motility. Phylogenomic analyses support the recent rRNA-based analyses that led to T. roseum being removed from the phylum Thermomicrobia and assigned to the phylum Chloroflexi. Because T. roseum is a deep-branching member of this phylum, analysis of its genome provides insights into the evolution of the Chloroflexi. In addition, even though this species is not photosynthetic, analysis of the genome provides some insight into the origins of photosynthesis in the Chloroflexi. Metabolic pathway reconstructions and experimental studies revealed new aspects of the biology of this species. For example, we present evidence that T. roseum oxidizes CO aerobically, making it the first thermophile known to do so. In addition, we propose that glycosylation of its carotenoids plays a crucial role in the adaptation of the cell membrane to this bacterium's thermophilic lifestyle. Analyses of published metagenomic sequences from two hot springs similar to the one from which this strain was isolated, show that close relatives of T. roseum DSM 5159 are present but have some key differences from the strain sequenced.

IDEA: Interactive Display for Evolutionary Analyses

BACKGROUND

The availability of complete genomic sequences for hundreds of organisms promises to make obtaining genome-wide estimates of substitution rates, selective constraints and other molecular evolution variables of interest an increasingly important approach to addressing broad evolutionary questions. Two of the programs most widely used for this purpose are codeml and baseml, parts of the PAML (Phylogenetic Analysis by Maximum Likelihood) suite. A significant drawback of these programs is their lack of a graphical user interface, which can limit their user base and considerably reduce their efficiency.

RESULTS

We have developed IDEA (Interactive Display for Evolutionary Analyses), an intuitive graphical input and output interface which interacts with PHYLIP for phylogeny reconstruction and with codeml and baseml for molecular evolution analyses. IDEA's graphical input and visualization interfaces eliminate the need to edit and parse text input and output files, reducing the likelihood of errors and improving processing time. Further, its interactive output display gives the user immediate access to results. Finally, IDEA can process data in parallel on a local machine or computing grid, allowing genome-wide analyses to be completed quickly.

CONCLUSION

IDEA provides a graphical user interface that allows the user to follow a codeml or baseml analysis from parameter input through to the exploration of results. Novel options streamline the analysis process, and post-analysis visualization of phylogenies, evolutionary rates and selective constraint along protein sequences simplifies the interpretation of results. The integration of these functions into a single tool eliminates the need for lengthy data handling and parsing, significantly expediting access to global patterns in the data.

Genome sequencing and analysis of the filamentous fungus Penicillium chrysogenum

Industrial penicillin production with the filamentous fungus Penicillium chrysogenum is based on an unprecedented effort in microbial strain improvement. To gain more insight into penicillin synthesis, we sequenced the 32.19 Mb genome of P. chrysogenum Wisconsin54-1255 and identified numerous genes responsible for key steps in penicillin production. DNA microarrays were used to compare the transcriptomes of the sequenced strain and a penicillinG high-producing strain, grown in the presence and absence of the side-chain precursor phenylacetic acid. Transcription of genes involved in biosynthesis of valine, cysteine and alpha-aminoadipic acid-precursors for penicillin biosynthesis-as well as of genes encoding microbody proteins, was increased in the high-producing strain. Some gene products were shown to be directly controlling beta-lactam output. Many key cellular transport processes involving penicillins and intermediates remain to be characterized at the molecular level. Genes predicted to encode transporters were strongly overrepresented among the genes transcriptionally upregulated under conditions that stimulate penicillinG production, illustrating potential for future genomics-driven metabolic engineering.

Complete genome sequence of the N2-fixing broad host range endophyte Klebsiella pneumoniae 342 and virulence predictions verified in mice

We report here the sequencing and analysis of the genome of the nitrogen-fixing endophyte, Klebsiella pneumoniae 342. Although K. pneumoniae 342 is a member of the enteric bacteria, it serves as a model for studies of endophytic, plant-bacterial associations due to its efficient colonization of plant tissues (including maize and wheat, two of the most important crops in the world), while maintaining a mutualistic relationship that encompasses supplying organic nitrogen to the host plant. Genomic analysis examined K. pneumoniae 342 for the presence of previously identified genes from other bacteria involved in colonization of, or growth in, plants. From this set, approximately one-third were identified in K. pneumoniae 342, suggesting additional factors most likely contribute to its endophytic lifestyle. Comparative genome analyses were used to provide new insights into this question. Results included the identification of metabolic pathways and other features devoted to processing plant-derived cellulosic and aromatic compounds, and a robust complement of transport genes (15.4%), one of the highest percentages in bacterial genomes sequenced. Although virulence and antibiotic resistance genes were predicted, experiments conducted using mouse models showed pathogenicity to be attenuated in this strain. Comparative genomic analyses with the presumed human pathogen K. pneumoniae MGH78578 revealed that MGH78578 apparently cannot fix nitrogen, and the distribution of genes essential to surface attachment, secretion, transport, and regulation and signaling varied between each genome, which may indicate critical divergences between the strains that influence their preferred host ranges and lifestyles (endophytic plant associations for K. pneumoniae 342 and presumably human pathogenesis for MGH78578). Little genome information is available concerning endophytic bacteria. The K. pneumoniae 342 genome will drive new research into this less-understood, but important category of bacterial-plant host relationships, which could ultimately enhance growth and nutrition of important agricultural crops and development of plant-derived products and biofuels.

Genomic islands in the pathogenic filamentous fungus Aspergillus fumigatus

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”.

Aspergillus is an extremely diverse genus of filamentous ascomycetous fungi (molds) found ubiquitously in soil and decomposing vegetation. Being supreme opportunists, aspergilli have adapted to overcome various chemical, physical, and biological stresses found in heterogeneous environments. While most species in the genus are saprophytes, a surprising number are able to infect wounded plants and animals. Remarkably, the allergic human host also responds abnormally to the aspergilli with lung and sinus disease. The advent of immunosuppressive agents and other medical advances have created a large worldwide pool of human hosts susceptible to some Aspergillus species, including the world's most harmful mold and the causative agent of invasive aspergillosis, Aspergillus fumigatus. In this study, we have used the power of comparative genomics to gain insight into genetic mechanisms that may contribute to the metabolic versatility and pathogenicity of this important human pathogen. Comparison of the genomes of two A. fumigatus clinical isolates and two closely related, but rarely pathogenic species showed that their genomes contain several large isolate- and species-specific chromosomal islands. The metabolic capabilities encoded by these highly labile regions are likely to contribute to their rapid adaptation to heterogeneous environments such as soil or a living host.

Genome sequence and identification of candidate vaccine antigens from the animal pathogen Dichelobacter nodosus

Dichelobacter nodosus causes ovine footrot, a disease that leads to severe economic losses in the wool and meat industries. We sequenced its 1.4-Mb genome, the smallest known genome of an anaerobe. It differs markedly from small genomes of intracellular bacteria, retaining greater biosynthetic capabilities and lacking any evidence of extensive ongoing genome reduction. Comparative genomic microarray studies and bioinformatic analysis suggested that, despite its small size, almost 20% of the genome is derived from lateral gene transfer. Most of these regions seem to be associated with virulence. Metabolic reconstruction indicated unsuspected capabilities, including carbohydrate utilization, electron transfer and several aerobic pathways. Global transcriptional profiling and bioinformatic analysis enabled the prediction of virulence factors and cell surface proteins. Screening of these proteins against ovine antisera identified eight immunogenic proteins that are candidate antigens for a cross-protective vaccine.

Macronuclear genome sequence of the ciliate Tetrahymena thermophila, a model eukaryote

The ciliate Tetrahymena thermophila is a model organism for molecular and cellular biology. Like other ciliates, this species has separate germline and soma functions that are embodied by distinct nuclei within a single cell. The germline-like micronucleus (MIC) has its genome held in reserve for sexual reproduction. The soma-like macronucleus (MAC), which possesses a genome processed from that of the MIC, is the center of gene expression and does not directly contribute DNA to sexual progeny. We report here the shotgun sequencing, assembly, and analysis of the MAC genome of T. thermophila, which is approximately 104 Mb in length and composed of approximately 225 chromosomes. Overall, the gene set is robust, with more than 27,000 predicted protein-coding genes, 15,000 of which have strong matches to genes in other organisms. The functional diversity encoded by these genes is substantial and reflects the complexity of processes required for a free-living, predatory, single-celled organism. This is highlighted by the abundance of lineage-specific duplications of genes with predicted roles in sensing and responding to environmental conditions (e.g., kinases), using diverse resources (e.g., proteases and transporters), and generating structural complexity (e.g., kinesins and dyneins). In contrast to the other lineages of alveolates (apicomplexans and dinoflagellates), no compelling evidence could be found for plastid-derived genes in the genome. UGA, the only T. thermophila stop codon, is used in some genes to encode selenocysteine, thus making this organism the first known with the potential to translate all 64 codons in nuclear genes into amino acids. We present genomic evidence supporting the hypothesis that the excision of DNA from the MIC to generate the MAC specifically targets foreign DNA as a form of genome self-defense. The combination of the genome sequence, the functional diversity encoded therein, and the presence of some pathways missing from other model organisms makes T. thermophila an ideal model for functional genomic studies to address biological, biomedical, and biotechnological questions of fundamental importance.

Comparative genomic evidence for a close relationship between the dimorphic prosthecate bacteria Hyphomonas neptunium and Caulobacter crescentus

The dimorphic prosthecate bacteria (DPB) are alpha-proteobacteria that reproduce in an asymmetric manner rather than by binary fission and are of interest as simple models of development. Prior to this work, the only member of this group for which genome sequence was available was the model freshwater organism Caulobacter crescentus. Here we describe the genome sequence of Hyphomonas neptunium, a marine member of the DPB that differs from C. crescentus in that H. neptunium uses its stalk as a reproductive structure. Genome analysis indicates that this organism shares more genes with C. crescentus than it does with Silicibacter pomeroyi (a closer relative according to 16S rRNA phylogeny), that it relies upon a heterotrophic strategy utilizing a wide range of substrates, that its cell cycle is likely to be regulated in a similar manner to that of C. crescentus, and that the outer membrane complements of H. neptunium and C. crescentus are remarkably similar. H. neptunium swarmer cells are highly motile via a single polar flagellum. With the exception of cheY and cheR, genes required for chemotaxis were absent in the H. neptunium genome. Consistent with this observation, H. neptunium swarmer cells did not respond to any chemotactic stimuli that were tested, which suggests that H. neptunium motility is a random dispersal mechanism for swarmer cells rather than a stimulus-controlled navigation system for locating specific environments. In addition to providing insights into bacterial development, the H. neptunium genome will provide an important resource for the study of other interesting biological processes including chromosome segregation, polar growth, and cell aging.

Genome sequence of Synechococcus CC9311: Insights into adaptation to a coastal environment

Coastal aquatic environments are typically more highly productive and dynamic than open ocean ones. Despite these differences, cyanobacteria from the genus Synechococcus are important primary producers in both types of ecosystems. We have found that the genome of a coastal cyanobacterium, Synechococcus sp. strain CC9311, has significant differences from an open ocean strain, Synechococcus sp. strain WH8102, and these are consistent with the differences between their respective environments. CC9311 has a greater capacity to sense and respond to changes in its (coastal) environment. It has a much larger capacity to transport, store, use, or export metals, especially iron and copper. In contrast, phosphate acquisition seems less important, consistent with the higher concentration of phosphate in coastal environments. CC9311 is predicted to have differences in its outer membrane lipopolysaccharide, and this may be characteristic of the speciation of some cyanobacterial groups. In addition, the types of potentially horizontally transferred genes are markedly different between the coastal and open ocean genomes and suggest a more prominent role for phages in horizontal gene transfer in oligotrophic environments.

Comparative analysis of programmed cell death pathways in filamentous fungi

BACKGROUND

Fungi can undergo autophagic- or apoptotic-type programmed cell death (PCD) on exposure to antifungal agents, developmental signals, and stress factors. Filamentous fungi can also exhibit a form of cell death called heterokaryon incompatibility (HI) triggered by fusion between two genetically incompatible individuals. With the availability of recently sequenced genomes of Aspergillus fumigatus and several related species, we were able to define putative components of fungi-specific death pathways and the ancestral core apoptotic machinery shared by all fungi and metazoa.

RESULTS

Phylogenetic profiling of HI-associated proteins from four Aspergilli and seven other fungal species revealed lineage-specific protein families, orphan genes, and core genes conserved across all fungi and metazoa. The Aspergilli-specific domain architectures include NACHT family NTPases, which may function as key integrators of stress and nutrient availability signals. They are often found fused to putative effector domains such as Pfs, SesB/LipA, and a newly identified domain, HET-s/LopB. Many putative HI inducers and mediators are specific to filamentous fungi and not found in unicellular yeasts. In addition to their role in HI, several of them appear to be involved in regulation of cell cycle, development and sexual differentiation. Finally, the Aspergilli possess many putative downstream components of the mammalian apoptotic machinery including several proteins not found in the model yeast, Saccharomyces cerevisiae.

CONCLUSION

Our analysis identified more than 100 putative PCD associated genes in the Aspergilli, which may help expand the range of currently available treatments for aspergillosis and other invasive fungal diseases. The list includes species-specific protein families as well as conserved core components of the ancestral PCD machinery shared by fungi and metazoa.

Genomic analysis of Hyphomonas neptunium contradicts 16S rRNA gene-based phylogenetic analysis: implications for the taxonomy of the orders 'Rhodobacterales' and Caulobacterales

Hyphomonas neptunium is a marine prosthecate alpha-proteobacterium currently classified as a member of the order 'Rhodobacterales'. Although this classification is supported by 16S rRNA gene sequence phylogeny, 23S rRNA gene sequence analysis, concatenated ribosomal proteins, HSP70 and EF-Tu phylogenies all support classifying Hyphomonas neptunium as a member of the Caulobacterales instead. The possible reasons why the 16S rRNA gene sequence gives conflicting results in this case are also discussed.

Major structural differences and novel potential virulence mechanisms from the genomes of multiple campylobacter species

Sequencing and comparative genome analysis of four strains of Campylobacter including C. lari RM2100, C. upsaliensis RM3195, and C. coli RM2228 has revealed major structural differences that are associated with the insertion of phage- and plasmid-like genomic islands, as well as major variations in the lipooligosaccharide complex. Poly G tracts are longer, are greater in number, and show greater variability in C. upsaliensis than in the other species. Many genes involved in host colonization, including racR/S, cadF, cdt, ciaB, and flagellin genes, are conserved across the species, but variations that appear to be species specific are evident for a lipooligosaccharide locus, a capsular (extracellular) polysaccharide locus, and a novel Campylobacter putative licABCD virulence locus. The strains also vary in their metabolic profiles, as well as their resistance profiles to a range of antibiotics. It is evident that the newly identified hypothetical and conserved hypothetical proteins, as well as uncharacterized two-component regulatory systems and membrane proteins, may hold additional significant information on the major differences in virulence among the species, as well as the specificity of the strains for particular hosts.

An information-based sequence distance and its application to whole mitochondrial genome phylogeny

MOTIVATION

Traditional sequence distances require an alignment and therefore are not directly applicable to the problem of whole genome phylogeny where events such as rearrangements make full length alignments impossible. We present a sequence distance that works on unaligned sequences using the information theoretical concept of Kolmogorov complexity and a program to estimate this distance.

RESULTS

We establish the mathematical foundations of our distance and illustrate its use by constructing a phylogeny of the Eutherian orders using complete unaligned mitochondrial genomes. This phylogeny is consistent with the commonly accepted one for the Eutherians. A second, larger mammalian dataset is also analyzed, yielding a phylogeny generally consistent with the commonly accepted one for the mammals.

AVAILABILITY

The program to estimate our sequence distance, is available at http://www.cs.cityu.edu.hk/~cssamk/gencomp/GenCompress1.htm. The distance matrices used to generate our phylogenies are available at http://www.math.uwaterloo.ca/~mli/distance.html.

CRITICA: coding region identification tool invoking comparative analysis

Gene recognition is essential to understanding existing and future DNA sequence data. CRITICA (Coding Region Identification Tool Invoking Comparative Analysis) is a suite of programs for identifying likely protein-coding sequences in DNA by combining comparative analysis of DNA sequences with more common noncomparative methods. In the comparative component of the analysis, regions of DNA are aligned with related sequences from the DNA databases; if the translation of the aligned sequences has greater amino acid identity than expected for the observed percentage nucleotide identity, this is interpreted as evidence for coding. CRITICA also incorporates noncomparative information derived from the relative frequencies of hexanucleotides in coding frames versus other contexts (i.e., dicodon bias). The dicodon usage information is derived by iterative analysis of the data, such that CRITICA is not dependent on the existence or accuracy of coding sequence annotations in the databases. This independence makes the method particularly well suited for the analysis of novel genomes. CRITICA was tested by analyzing the available Salmonella typhimurium DNA sequences. Its predictions were compared with the DNA sequence annotations and with the predictions of GenMark. CRITICA proved to be more accurate than GenMark, and moreover, many of its predictions that would seem to be errors instead reflect problems in the sequence databases. The source code of CRITICA is freely available by anonymous FTP (rdp.life.uiuc.edu in/pub/critica) and on the World Wide Web (http:/(/)rdpwww.life.uiuc.edu).

Thermal adaptation analyzed by comparison of protein sequences from mesophilic and extremely thermophilic Methanococcus species

The genome sequence of the extremely thermophilic archaeon Methanococcus jannaschii provides a wealth of data on proteins from a thermophile. In this paper, sequences of 115 proteins from M. jannaschii are compared with their homologs from mesophilic Methanococcus species. Although the growth temperatures of the mesophiles are about 50 degrees C below that of M. jannaschii, their genomic G+C contents are nearly identical. The properties most correlated with the proteins of the thermophile include higher residue volume, higher residue hydrophobicity, more charged amino acids (especially Glu, Arg, and Lys), and fewer uncharged polar residues (Ser, Thr, Asn, and Gln). These are recurring themes, with all trends applying to 83-92% of the proteins for which complete sequences were available. Nearly all of the amino acid replacements most significantly correlated with the temperature change are the same relatively conservative changes observed in all proteins, but in the case of the mesophile/thermophile comparison there is a directional bias. We identify 26 specific pairs of amino acids with a statistically significant (P < 0.01) preferred direction of replacement.

The complete genome sequence of the hyperthermophilic, sulphate-reducing archaeon Archaeoglobus fulgidus

Archaeoglobus fulgidus is the first sulphur-metabolizing organism to have its genome sequence determined. Its genome of 2,178,400 base pairs contains 2,436 open reading frames (ORFs). The information processing systems and the biosynthetic pathways for essential components (nucleotides, amino acids and cofactors) have extensive correlation with their counterparts in the archaeon Methanococcus jannaschii. The genomes of these two Archaea indicate dramatic differences in the way these organisms sense their environment, perform regulatory and transport functions, and gain energy. In contrast to M. jannaschii, A. fulgidus has fewer restriction-modification systems, and none of its genes appears to contain inteins. A quarter (651 ORFs) of the A. fulgidus genome encodes functionally uncharacterized yet conserved proteins, two-thirds of which are shared with M. jannaschii (428 ORFs). Another quarter of the genome encodes new proteins indicating substantial archaeal gene diversity.