Cellulosomics, a gene-centric approach to investigating the intraspecific diversity and adaptation of Ruminococcus flavefaciens within the rumen

Background

The bovine rumen maintains a diverse microbial community that serves to break down indigestible plant substrates. However, those bacteria specifically adapted to degrade cellulose, the major structural component of plant biomass, represent a fraction of the rumen microbiome. Previously, we proposed scaC as a candidate for phylotyping Ruminococcus flavefaciens, one of three major cellulolytic bacterial species isolated from the rumen. In the present report we examine the dynamics and diversity of scaC-types both within and between cattle temporally, following a dietary switch from corn-silage to grass-legume hay. These results were placed in the context of the overall bacterial population dynamics measured using the 16S rRNA.

Principal Findings

As many as 117 scaC-types were estimated, although just nineteen were detected in each of three rumens tested, and these collectively accounted for the majority of all types present. Variation in scaC populations was observed between cattle, between planktonic and fiber-associated fractions and temporally over the six-week survey, and appeared related to scaC phylogeny. However, by the sixth week no significant separation of scaC populations was seen between animals, suggesting enrichment of a constrained set of scaC-types. Comparing the amino-acid translation of each scaC-type revealed sequence variation within part of the predicted dockerin module but strong conservation in the N-terminus, where the cohesin module is located.

Conclusions

The R. flavefaciens species comprises a multiplicity of scaC-types in-vivo. Enrichment of particular scaC-types temporally, following a dietary switch, and between fractions along with the phylogenetic congruence suggests that functional differences exist between types. Observed differences in dockerin modules suggest at least part of the functional heterogeneity may be conferred by scaC. The polymorphic nature of scaC enables the relative distribution of R. flavefaciens strains to be examined and represents a gene-centric approach to investigating the intraspecific adaptation of an important specialist population.     

Diversity and Strain Specificity of Plant Cell Wall Degrading Enzymes Revealed by the Draft Genome of Ruminococcus flavefaciens FD-1

Background

Ruminococcus flavefaciens is a predominant cellulolytic rumen bacterium, which forms a multi-enzyme cellulosome complex that could play an integral role in the ability of this bacterium to degrade plant cell wall polysaccharides. Identifying the major enzyme types involved in plant cell wall degradation is essential for gaining a better understanding of the cellulolytic capabilities of this organism as well as highlighting potential enzymes for application in improvement of livestock nutrition and for conversion of cellulosic biomass to liquid fuels.

Methodology/Principal Findings

The R. flavefaciens FD-1 genome was sequenced to 29x-coverage, based on pulsed-field gel electrophoresis estimates (4.4 Mb), and assembled into 119 contigs providing 4,576,399 bp of unique sequence. As much as 87.1% of the genome encodes ORFs, tRNA, rRNAs, or repeats. The GC content was calculated at 45%. A total of 4,339 ORFs was detected with an average gene length of 918 bp. The cellulosome model for R. flavefaciens was further refined by sequence analysis, with at least 225 dockerin-containing ORFs, including previously characterized cohesin-containing scaffoldin molecules. These dockerin-containing ORFs encode a variety of catalytic modules including glycoside hydrolases (GHs), polysaccharide lyases, and carbohydrate esterases. Additionally, 56 ORFs encode proteins that contain carbohydrate-binding modules (CBMs). Functional microarray analysis of the genome revealed that 56 of the cellulosome-associated ORFs were up-regulated, 14 were down-regulated, 135 were unaffected, when R. flavefaciens FD-1 was grown on cellulose versus cellobiose. Three multi-modular xylanases ({"type":"entrez-protein","attrs":{"text":"ORF01222","term_id":"1178790230","term_text":"ORF01222"}}ORF01222, {"type":"entrez-protein","attrs":{"text":"ORF03896","term_id":"1178792974","term_text":"ORF03896"}}ORF03896, and {"type":"entrez-protein","attrs":{"text":"ORF01315","term_id":"1178790325","term_text":"ORF01315"}}ORF01315) exhibited the highest levels of up-regulation.

Conclusions/Significance

The genomic evidence indicates that R. flavefaciens FD-1 has the largest known number of fiber-degrading enzymes likely to be arranged in a cellulosome architecture. Functional analysis of the genome has revealed that the growth substrate drives expression of enzymes predicted to be involved in carbohydrate metabolism as well as expression and assembly of key cellulosomal enzyme components.     

Abundance and Diversity of Dockerin-Containing Proteins in the Fiber-Degrading Rumen Bacterium,Ruminococcus flavefaciens FD-1

Background

The cellulosome is a multi-enzyme machine, which plays a key role in the breakdown of plant cell walls in many anaerobic cellulose-degrading microorganisms. Ruminococcus flavefaciens FD-1, a major fiber-degrading bacterium present in the gut of herbivores, has the most intricate cellulosomal organization thus far described. Cellulosome complexes are assembled through high-affinity cohesin-dockerin interactions. More than two-hundred dockerin-containing proteins have been identified in the R. flavefaciens genome, yet the reason for the expansion of these crucial cellulosomal components is yet unknown.

Methodology/Principal Findings

We have explored the full spectrum of 222 dockerin-containing proteins potentially involved in the assembly of cellulosome-like complexes of R. flavefaciens. Bioinformatic analysis of the various dockerin modules showed distinctive conservation patterns within their two Ca²⁺-binding repeats and their flanking regions. Thus, we established the conceptual framework for six major groups of dockerin types, according to their unique sequence features. Within this framework, the modular architecture of the parent proteins, some of which are multi-functional proteins, was evaluated together with their gene expression levels. Specific dockerin types were found to be associated with selected groups of functional components, such as carbohydrate-binding modules, numerous peptidases, and/or carbohydrate-active enzymes. In addition, members of other dockerin groups were linked to structural proteins, e.g., cohesin-containing proteins, belonging to the scaffoldins.

Conclusions/Significance

This report profiles the abundance and sequence diversity of the R. flavefaciens FD-1 dockerins, and provides the molecular basis for future understanding of the potential for a wide array of cohesin-dockerin specificities. Conserved differences between dockerins may be reflected in their stability, function or expression within the context of the parent protein, in response to their role in the rumen environment.     

Magnesium requirement of some of the principal rumen cellulolytic bacteria

Information available on the role of Mg for growth and cellulose degradation by rumen bacteria is both limited and inconsistent. In this study, the Mg requirements for two strains each of the cellulolytic rumen species Fibrobacter succinogenes (A3c and S85), Ruminococcus albus (7 and 8) and Ruminococcus flavefaciens (B34b and C94) were investigated. Maximum growth, rate of growth and lag time were all measured using a complete factorial design, 2(3)×6; factors were: strains (2), within species (3) and Mg concentrations (6). R. flavefaciens was the only species that did not grow when Mg was singly deleted from the media, and both strains exhibited a linear growth response to increasing Mg concentrations (P<0.001). The requirement for R. flavefaciens B34b was estimated as 0.54 mM; whereas the requirement for R. flavefaciens C94 was >0.82 as there was no plateau in growth. Although not an absolute requirement for growth, strains of the two other species of cellulolytic bacteria all responded to increasing Mg concentrations. For F. succinogenes S85, R. albus 7 and R. albus 8, their requirement estimated from maximum growth was 0.56, 0.52 and 0.51, respectively. A requirement for F. succinogenes A3c could not be calculated because there was no solution for contrasts. Whether R. flavefaciens had a Mg requirement for cellulose degradation was determined in NH₃-free cellulose media, using a 2×4 factorial design, 2 strains and 4 treatments. Both strains of R. flavefaciens were found to have an absolute Mg requirement for cellulose degradation. Based on reported concentrations of Mg in the rumen, 1.0 to 10.1 mM, it seems unlikely that an in vivo deficiency of this element would occur.     

Metasecretome-selective phage display approach for mining the functional potential of a rumen microbial community

Background

In silico, secretome proteins can be predicted from completely sequenced genomes using various available algorithms that identify membrane-targeting sequences. For metasecretome (collection of surface, secreted and transmembrane proteins from environmental microbial communities) this approach is impractical, considering that the metasecretome open reading frames (ORFs) comprise only 10% to 30% of total metagenome, and are poorly represented in the dataset due to overall low coverage of metagenomic gene pool, even in large-scale projects.

Results

By combining secretome-selective phage display and next-generation sequencing, we focused the sequence analysis of complex rumen microbial community on the metasecretome component of the metagenome. This approach achieved high enrichment (29 fold) of secreted fibrolytic enzymes from the plant-adherent microbial community of the bovine rumen. In particular, we identified hundreds of heretofore rare modules belonging to cellulosomes, cell-surface complexes specialised for recognition and degradation of the plant fibre.

Conclusions

As a method, metasecretome phage display combined with next-generation sequencing has a power to sample the diversity of low-abundance surface and secreted proteins that would otherwise require exceptionally large metagenomic sequencing projects. As a resource, metasecretome display library backed by the dataset obtained by next-generation sequencing is ready for i) affinity selection by standard phage display methodology and ii) easy purification of displayed proteins as part of the virion for individual functional analysis.

Electronic supplementary material

The online version of this article (doi:10.1186/1471-2164-15-356) contains supplementary material, which is available to authorized users.     

Cellulosilyticum ruminicola,a Newly Described Rumen Bacterium That Possesses Redundant Fibrolytic-Protein-Encoding Genes and Degrades Lignocellulose with Multiple Carbohydrate- Borne Fibrolytic Enzymes

Cellulosilyticum ruminicola H1 is a newly described bacterium isolated from yak (Bos grunniens) rumen and is characterized by its ability to grow on a variety of hemicelluloses and degrade cellulosic materials. In this study, we performed the whole-genome sequencing of C. ruminicola H1 and observed a comprehensive set of genes encoding the enzymes essential for hydrolyzing plant cell wall. The corresponding enzymatic activities were also determined in strain H1; these included endoglucanases, cellobiohydrolases, xylanases, mannanase, pectinases, and feruloyl esterases and acetyl esterases to break the interbridge cross-link, as well as the enzymes that degrade the glycosidic bonds. This bacterium appears to produce polymer hydrolases that act on both soluble and crystal celluloses. Approximately half of the cellulytic activities, including cellobiohydrolase (50%), feruloyl esterase (45%), and one third of xylanase (31%) and endoglucanase (36%) activities were bound to cellulosic fibers. However, only a minority of mannase (6.78%) and pectinase (1.76%) activities were fiber associated. Strain H1 seems to degrade the plant-derived polysaccharides by producing individual fibrolytic enzymes, whereas the majority of polysaccharide hydrolases contain carbohydrate-binding module. Cellulosome or cellulosomelike protein complex was never isolated from this bacterium. Thus, the fibrolytic enzyme production of strain H1 may represent a different strategy in cellulase organization used by most of other ruminal microbes, but it applies the fungal mode of cellulose production.The ruminant rumens are long believed to be the anaerobic environments efficiently degrading the plant-derived polysaccharides, which is attributed to the inhabited abundant rumen microorganisms. They implement the fibrolytic degradation by the combination of the enzymes comprising of cellulases, hemicellulases, and to a lesser extent pectinases and ligninases (12). The rumen bacteria are outnumbered of the other rumen microbes; however, only a few of cellulolytic bacteria have been isolated from rumens. Ruminococcus flavefaciens, Ruminococcus albus, and Fibrobacter succinogenes are considered to be the most important cellulose-degrading bacteria in the rumen (18), and they produce a set of cellulolytic enzymes, including endoglucanases, exoglucanases (generally cellobiohydrolase), and β-glucosidases, as well as hemicellulases. In addition, the predominant ruminal hemicellulose-digesting bacteria such as Butyrivibrio fibrisolvens and Prevotella ruminicola lack the ability to digest cellulose but degrade xylan and pectin and utilize the degraded soluble sugars as substrates (10, 14). Although the robust cellulolytic species F. succinogenes degrades xylan, it cannot use the pentose product as a carbon source (24). Culture-independent approaches indicate that the three cellulolytic bacterial species represent only ∼2% of the ruminal bacterial 16S rRNA (43). Therefore, many varieties of rumen microbes remain uncultured (2). In recent years, rumen metagenomics studies have revealed the vast diversity of fibrolytic enzymes, multiple domain proteins, and the complexity of microbial composition in the ecosystem (9, 17). Hence, it is likely that the entire microbial community is necessary for the implementation of an efficient fibrolytic process in the rumen, including the uncultured species.In the rumen and other fibrolytic ecosystems, cellulolytic bacteria have to cope with the structural complexity of lignocelluloses and the interspecies competition; thus, not only a variety of plant polymer-degrading enzymes but also a noncatalytic assistant strategy, such as including adhesion of cells to substrates by a variety of anchoring domains, is required (8, 33, 38, 39). The (hemi)cellulolytic enzyme systems have been intensively studied for nonrumen anaerobic bacteria, including Clostridium thermocellum (19, 40), Clostridium cellulolyticum (6), Clostridium cellulovorans (13), and Clostridium stercorarium (47), as well as the rumen species, Rumicoccocus albus (35), Ruminococcus flavefaciens (32), and Fibrobacter succinogenes (4). The results indicate that most of them, except for Fibrobacter succinogenes, produce multiple cellulolytic enzymes integrated in a complex, cellulosome, and free individual proteins.The yak (Bos grunniens) is a large ruminant (∼1,000 kg) in the bovine family that lives mainly on the Qinghai-Tibetan Plateau in China at an altitude of 3,000 m above sea level. It is a local species that lives mainly on the world''s highest plateau. Yaks live in a full-grazing style with grasses, straws, and lichens as their exclusive feed, so the yak rumen can harbor a microbial flora distinct from those of other ruminants due to their fiber-component diet, since diet can be a powerful factor in regulating mammalian gut microbiome (27). A very different prokaryote community structure was revealed for yak rumen in our previous work based on the 16S rRNA diversity, which showed fewer phyla than for cattle but that a higher ratio of sequences was related to uncultured bacteria (2).We previously isolated a novel anaerobic fibrolytic bacterium, Cellulosilyticum ruminicola H1, from the rumen of a domesticated yak (11). Strain H1 grew robustly on natural plant fibers such as corn cob, alfalfa, and ryegrass as the sole carbon and energy sources, as well as on a variety of polysaccharides, including cellulose, xylan, mannan, and pectin, but not monosaccharides such as glucose, which is preferred by most ruminal bacteria. In the present study, using a draft of its genome and enzymatic characterization, we analyzed the enzymatic activities and the structures of the polymer hydrolases of strain H1 that were involved in the hydrolysis of complex polysaccharides.     

Comparative genome sequencing reveals chemotype-specific gene clusters in the toxigenic black mold Stachybotrys

Background

The fungal genus Stachybotrys produces several diverse toxins that affect human health. Its strains comprise two mutually-exclusive toxin chemotypes, one producing satratoxins, which are a subclass of trichothecenes, and the other producing the less-toxic atranones. To determine the genetic basis for chemotype-specific differences in toxin production, the genomes of four Stachybotrys strains were sequenced and assembled de novo. Two of these strains produce atranones and two produce satratoxins.

Results

Comparative analysis of these four 35-Mbp genomes revealed several chemotype-specific gene clusters that are predicted to make secondary metabolites. The largest, which was named the core atranone cluster, encodes 14 proteins that may suffice to produce all observed atranone compounds via reactions that include an unusual Baeyer-Villiger oxidation. Satratoxins are suggested to be made by products of multiple gene clusters that encode 21 proteins in all, including polyketide synthases, acetyltransferases, and other enzymes expected to modify the trichothecene skeleton. One such satratoxin chemotype-specific cluster is adjacent to the core trichothecene cluster, which has diverged from those of other trichothecene producers to contain a unique polyketide synthase.

Conclusions

The results suggest that chemotype-specific gene clusters are likely the genetic basis for the mutually-exclusive toxin chemotypes of Stachybotrys. A unified biochemical model for Stachybotrys toxin production is presented. Overall, the four genomes described here will be useful for ongoing studies of this mold’s diverse toxicity mechanisms.

Electronic supplementary material

The online version of this article (doi:10.1186/1471-2164-15-590) contains supplementary material, which is available to authorized users.     

Identification of Ruminococcus flavefaciens as the Predominant Cellulolytic Bacterial Species of the Equine Cecum

Detection and quantification of cellulolytic bacteria with oligonucleotide probes showed that Ruminococcus flavefaciens was the predominant species in the pony and donkey cecum. Fibrobacter succinogenes and Ruminococcus albus were present at low levels. Four isolates, morphologically resembling R. flavefaciens, differed from ruminal strains by their carbohydrate utilization and their end products of cellobiose fermentation.     

Analysis of the Pantoea ananatis pan-genome reveals factors underlying its ability to colonize and interact with plant,insect and vertebrate hosts

Background

Pantoea ananatis is found in a wide range of natural environments, including water, soil, as part of the epi- and endophytic flora of various plant hosts, and in the insect gut. Some strains have proven effective as biological control agents and plant-growth promoters, while other strains have been implicated in diseases of a broad range of plant hosts and humans. By analysing the pan-genome of eight sequenced P. ananatis strains isolated from different sources we identified factors potentially underlying its ability to colonize and interact with hosts in both the plant and animal Kingdoms.

Results

The pan-genome of the eight compared P. ananatis strains consisted of a core genome comprised of 3,876 protein coding sequences (CDSs) and a sizeable accessory genome consisting of 1,690 CDSs. We estimate that ~106 unique CDSs would be added to the pan-genome with each additional P. ananatis genome sequenced in the future. The accessory fraction is derived mainly from integrated prophages and codes mostly for proteins of unknown function. Comparison of the translated CDSs on the P. ananatis pan-genome with the proteins encoded on all sequenced bacterial genomes currently available revealed that P. ananatis carries a number of CDSs with orthologs restricted to bacteria associated with distinct hosts, namely plant-, animal- and insect-associated bacteria. These CDSs encode proteins with putative roles in transport and metabolism of carbohydrate and amino acid substrates, adherence to host tissues, protection against plant and animal defense mechanisms and the biosynthesis of potential pathogenicity determinants including insecticidal peptides, phytotoxins and type VI secretion system effectors.

Conclusions

P. ananatis has an ‘open’ pan-genome typical of bacterial species that colonize several different environments. The pan-genome incorporates a large number of genes encoding proteins that may enable P. ananatis to colonize, persist in and potentially cause disease symptoms in a wide range of plant and animal hosts.

Electronic supplementary material

The online version of this article (doi: 10.1186/1471-2164-15-404) contains supplementary material, which is available to authorized users.     

Influence of Forage Phenolics on Ruminal Fibrolytic Bacteria and In Vitro Fiber Degradation

In vitro cultures of ruminal microorganisms were used to determine the effect of cinnamic acid and vanillin on the digestibility of cellulose and xylan. Cinnamic acid and vanillin depressed in vitro dry matter disappearance of cellulose 14 and 49%, respectively, when rumen fluid was the inoculum. The number of viable Bacteroides succinogenes cells, the predominant cellulolytic organism, was threefold higher for fermentations which contained vanillin than for control fermentations. When xylan replaced cellulose as the substrate, a 14% decrease in the digestibility of xylan was observed with vanillin added; however, the number of viable xylanolytic bacteria cultured from the batch fermentation was 10-fold greater than that of control fermentations. The doubling time of B. succinogenes was increased from 2.32 to 2.58 h when vanillin was added to cellobiose medium, and absorbance was one-half that of controls after 18 h. The growth rate of Ruminococcus albus and Ruminococcus flavefaciens was inhibited more by p-coumaric acid than by vanillin, although no reduction of final absorbance was observed in their growth cycles. Vanillin, and to a lesser extent cinnamic acid, appeared to prevent the attachment of B. succinogenes cells to cellulose particles, but did not affect dissociation of cells from the particles. B. succinogenes, R. albus, R. flavefaciens, and Butyrivibrio fibrisolvens all modified the parent monomers cinnamic acid, p-coumaric acid, ferulic acid, and vanillin, with B. fibrisolvens causing the most extensive modification. These results suggest that phenolic monomers can inhibit digestibility of cellulose and xylan, possibly by influencing attachment of the fibrolytic microorganisms to fiber particles. The reduced bacterial attachment to structural carbohydrates in the presence of vanillin may generate more free-floating fibrolytic organisms, thus giving a deceptively higher viable count.     

Effects of phosphorus supply on growth, phosphate concentration and cluster-root formation in three Lupinus species

Background and Aims

In some lupin species, phosphate deficiency induces cluster-root formation, which enhances P uptake by increasing root surface area and, more importantly, the release of root exudates which enhances P availability.

Methods

Three species of Lupinus, L. albus, L. atlanticus and L. micranthus, with inherently different relative growth rates were cultivated under hydroponics in a greenhouse at four phosphate concentrations (1, 10, 50 and 150 µm) to compare the role of internal P in regulating cluster-root formation.

Key Results

The highest growth rate was observed in L. atlanticus, followed by L. albus and L. micranthus. At 1 µm P, cluster-root formation was markedly induced in all three species. The highest P uptake and accumulation was observed in L. micranthus, followed by L. atlanticus and then L. albus. Inhibition of cluster-root formation was severe at 10 µm P in L. atlanticus, but occurred stepwise with increasing P concentration in the root medium in L. albus.

Conclusions

In L. atlanticus and L. albus cluster-root formation was suppressed by P treatments above 10 µm, indicating a P-inducible regulating system for cluster-root formation, as expected. By contrast, production of cluster roots in L. micranthus, in spite of a high internal P concentration, indicated a lower sensitivity to P status, which allowed P-toxicity symptoms to develop.     

Comparative genomics of the Bifidobacterium breve taxon

Background

Bifidobacteria are commonly found as part of the microbiota of the gastrointestinal tract (GIT) of a broad range of hosts, where their presence is positively correlated with the host’s health status. In this study, we assessed the genomes of thirteen representatives of Bifidobacterium breve, which is not only a frequently encountered component of the (adult and infant) human gut microbiota, but can also be isolated from human milk and vagina.

Results

In silico analysis of genome sequences from thirteen B. breve strains isolated from different environments (infant and adult faeces, human milk, human vagina) shows that the genetic variability of this species principally consists of hypothetical genes and mobile elements, but, interestingly, also genes correlated with the adaptation to host environment and gut colonization. These latter genes specify the biosynthetic machinery for sortase-dependent pili and exopolysaccharide production, as well as genes that provide protection against invasion of foreign DNA (i.e. CRISPR loci and restriction/modification systems), and genes that encode enzymes responsible for carbohydrate fermentation. Gene-trait matching analysis showed clear correlations between known metabolic capabilities and characterized genes, and it also allowed the identification of a gene cluster involved in the utilization of the alcohol-sugar sorbitol.

Conclusions

Genome analysis of thirteen representatives of the B. breve species revealed that the deduced pan-genome exhibits an essentially close trend. For this reason our analyses suggest that this number of B. breve representatives is sufficient to fully describe the pan-genome of this species. Comparative genomics also facilitated the genetic explanation for differential carbon source utilization phenotypes previously observed in different strains of B. breve.

Electronic supplementary material

The online version of this article (doi:10.1186/1471-2164-15-170) contains supplementary material, which is available to authorized users.     

LOCAS--a low coverage assembly tool for resequencing projects

Motivation

Next Generation Sequencing (NGS) is a frequently applied approach to detect sequence variations between highly related genomes. Recent large-scale re-sequencing studies as the Human 1000 Genomes Project utilize NGS data of low coverage to afford sequencing of hundreds of individuals. Here, SNPs and micro-indels can be detected by applying an alignment-consensus approach. However, computational methods capable of discovering other variations such as novel insertions or highly diverged sequence from low coverage NGS data are still lacking.

Results

We present LOCAS, a new NGS assembler particularly designed for low coverage assembly of eukaryotic genomes using a mismatch sensitive overlap-layout-consensus approach. LOCAS assembles homologous regions in a homology-guided manner while it performs de novo assemblies of insertions and highly polymorphic target regions subsequently to an alignment-consensus approach. LOCAS has been evaluated in homology-guided assembly scenarios with low sequence coverage of Arabidopsis thaliana strains sequenced as part of the Arabidopsis 1001 Genomes Project. While assembling the same amount of long insertions as state-of-the-art NGS assemblers, LOCAS showed best results regarding contig size, error rate and runtime.

Conclusion

LOCAS produces excellent results for homology-guided assembly of eukaryotic genomes with short reads and low sequencing depth, and therefore appears to be the assembly tool of choice for the detection of novel sequence variations in this scenario.     

Whole genome sequences of three Treponema pallidum ssp. pertenue strains: yaws and syphilis treponemes differ in less than 0.2% of the genome sequence

Background

The yaws treponemes, Treponema pallidum ssp. pertenue (TPE) strains, are closely related to syphilis causing strains of Treponema pallidum ssp. pallidum (TPA). Both yaws and syphilis are distinguished on the basis of epidemiological characteristics, clinical symptoms, and several genetic signatures of the corresponding causative agents.

Methodology/Principal Findings

To precisely define genetic differences between TPA and TPE, high-quality whole genome sequences of three TPE strains (Samoa D, CDC-2, Gauthier) were determined using next-generation sequencing techniques. TPE genome sequences were compared to four genomes of TPA strains (Nichols, DAL-1, SS14, Chicago). The genome structure was identical in all three TPE strains with similar length ranging between 1,139,330 bp and 1,139,744 bp. No major genome rearrangements were found when compared to the four TPA genomes. The whole genome nucleotide divergence (d_A) between TPA and TPE subspecies was 4.7 and 4.8 times higher than the observed nucleotide diversity (π) among TPA and TPE strains, respectively, corresponding to 99.8% identity between TPA and TPE genomes. A set of 97 (9.9%) TPE genes encoded proteins containing two or more amino acid replacements or other major sequence changes. The TPE divergent genes were mostly from the group encoding potential virulence factors and genes encoding proteins with unknown function.

Conclusions/Significance

Hypothetical genes, with genetic differences, consistently found between TPE and TPA strains are candidates for syphilitic treponemes virulence factors. Seventeen TPE genes were predicted under positive selection, and eleven of them coded either for predicted exported proteins or membrane proteins suggesting their possible association with the cell surface. Sequence changes between TPE and TPA strains and changes specific to individual strains represent suitable targets for subspecies- and strain-specific molecular diagnostics.     

Influence of Plant Phenolic Acids on Growth and Cellulolytic Activity of Rumen Bacteria

Isolated rumen bacteria were examined for growth and, where appropriate, for their ability to degrade cellulose in the presence of the hydroxycinnamic acids trans-p-coumaric acid and trans-ferulic acid and the hydroxybenzoic acids vanillic acid and 4-hydroxybenzoic acid. Ferulic and p-coumaric acids proved to be the most toxic of the acids examined and suppressed the growth of the cellulolytic strains Ruminococcus albus, Ruminococcus flavefaciens, and Bacteroides succinogenes when included in a simple sugars medium at concentrations of >5 mM. The extent of cellulose digestion by R. flavefaciens and B. succinogenes but not R. albus was also substantially reduced. Examination of rumen fluid from sheep maintained on dried grass containing 0.51% phenolic acids showed the presence of phloretic acid (0.1 mM) and 3-methoxyphloretic acid (trace) produced by hydrogenation of the 2-propenoic side chain of p-coumaric and ferulic acids, respectively. The parent acids were found in trace amounts only, although they represented the major phenolic acids ingested. Phloretic and 3-methoxyphloretic acids proved to be considerably less toxic than their parent acids. All of the cellulolytic strains (and Streptococcus bovis) showed at least a limited ability to hydrogenate hydroxycinnamic acids, with Ruminococcus spp. proving the most effective. No further modification of hydroxycinnamic acids was produced by the single strains of bacteria examined. However, a considerable shortfall in the recovery of added phenolic acids was noted in media inoculated with rumen fluid. It is suggested that hydrogenation may serve to protect cellulolytic strains from hydroxycinnamic acids.     

Comparative analyses of Legionella species identifies genetic features of strains causing Legionnaires’ disease

Background

The genus Legionella comprises over 60 species. However, L. pneumophila and L. longbeachae alone cause over 95% of Legionnaires’ disease. To identify the genetic bases underlying the different capacities to cause disease we sequenced and compared the genomes of L. micdadei, L. hackeliae and L. fallonii (LLAP10), which are all rarely isolated from humans.

Results

We show that these Legionella species possess different virulence capacities in amoeba and macrophages, correlating with their occurrence in humans. Our comparative analysis of 11 Legionella genomes belonging to five species reveals highly heterogeneous genome content with over 60% representing species-specific genes; these comprise a complete prophage in L. micdadei, the first ever identified in a Legionella genome. Mobile elements are abundant in Legionella genomes; many encode type IV secretion systems for conjugative transfer, pointing to their importance for adaptation of the genus. The Dot/Icm secretion system is conserved, although the core set of substrates is small, as only 24 out of over 300 described Dot/Icm effector genes are present in all Legionella species. We also identified new eukaryotic motifs including thaumatin, synaptobrevin or clathrin/coatomer adaptine like domains.

Conclusions

Legionella genomes are highly dynamic due to a large mobilome mainly comprising type IV secretion systems, while a minority of core substrates is shared among the diverse species. Eukaryotic like proteins and motifs remain a hallmark of the genus Legionella. Key factors such as proteins involved in oxygen binding, iron storage, host membrane transport and certain Dot/Icm substrates are specific features of disease-related strains.

Electronic supplementary material

The online version of this article (doi:10.1186/s13059-014-0505-0) contains supplementary material, which is available to authorized users.     

Characterization of the core and accessory genomes of Pseudomonas aeruginosa using bioinformatic tools Spine and AGEnt

Background

Pseudomonas aeruginosa is an important opportunistic pathogen responsible for many infections in hospitalized and immunocompromised patients. Previous reports estimated that approximately 10% of its 6.6 Mbp genome varies from strain to strain and is therefore referred to as “accessory genome”. Elements within the accessory genome of P. aeruginosa have been associated with differences in virulence and antibiotic resistance. As whole genome sequencing of bacterial strains becomes more widespread and cost-effective, methods to quickly and reliably identify accessory genomic elements in newly sequenced P. aeruginosa genomes will be needed.

Results

We developed a bioinformatic method for identifying the accessory genome of P. aeruginosa. First, the core genome was determined based on sequence conserved among the completed genomes of twelve reference strains using Spine, a software program developed for this purpose. The core genome was 5.84 Mbp in size and contained 5,316 coding sequences. We then developed an in silico genome subtraction program named AGEnt to filter out core genomic sequences from P. aeruginosa whole genomes to identify accessory genomic sequences of these reference strains. This analysis determined that the accessory genome of P. aeruginosa ranged from 6.9-18.0% of the total genome, was enriched for genes associated with mobile elements, and was comprised of a majority of genes with unknown or unclear function. Using these genomes, we showed that AGEnt performed well compared to other publically available programs designed to detect accessory genomic elements. We then demonstrated the utility of the AGEnt program by applying it to the draft genomes of two previously unsequenced P. aeruginosa strains, PA99 and PA103.

Conclusions

The P. aeruginosa genome is rich in accessory genetic material. The AGEnt program accurately identified the accessory genomes of newly sequenced P. aeruginosa strains, even when draft genomes were used. As P. aeruginosa genomes become available at an increasingly rapid pace, this program will be useful in cataloging the expanding accessory genome of this bacterium and in discerning correlations between phenotype and accessory genome makeup. The combination of Spine and AGEnt should be useful in defining the accessory genomes of other bacterial species as well.

Electronic supplementary material

The online version of this article (doi:10.1186/1471-2164-15-737) contains supplementary material, which is available to authorized users.     

Comparative genomic analysis of Ralstonia solanacearum reveals candidate genes for host specificity

Background

Ralstonia solanacearum is a vascular soil-borne plant pathogen with an unusually broad host range. This economically destructive and globally distributed bacterium has thousands of distinct lineages within a heterogeneous and taxonomically disputed species complex. Some lineages include highly host-adapted strains (ecotypes), such as the banana Moko disease-causing strains, the cold-tolerant potato brown rot strains (also known as R3bv2) and the recently emerged Not Pathogenic to Banana (NPB) strains.

Results

These distinct ecotypes offer a robust model to study host adaptation and the emergence of ecotypes because the polyphyletic Moko strains include lineages that are phylogenetically close to the monophyletic brown rot and NPB strains. Draft genomes of eight new strains belonging to these three model ecotypes were produced to complement the eleven publicly available R. solanacearum genomes. Using a suite of bioinformatics methods, we searched for genetic and evolutionary features that distinguish ecotypes and propose specific hypotheses concerning mechanisms of host adaptation in the R. solanacearum species complex. Genome-wide, few differences were identified, but gene loss events, non-synonymous polymorphisms, and horizontal gene transfer were identified among type III effectors and were associated with host range differences.

Conclusions

This extensive comparative genomics analysis uncovered relatively few divergent features among closely related strains with contrasting biological characteristics; however, several virulence factors were associated with the emergence of Moko, NPB and brown rot and could explain host adaptation.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-015-1474-8) contains supplementary material, which is available to authorized users.     

Comparative analysis of plant carbohydrate active enZymes and their role in xylogenesis

Background

Carbohydrate metabolism is a key feature of vascular plant architecture, and is of particular importance in large woody species, where lignocellulosic biomass is responsible for bearing the bulk of the stem and crown. Since Carbohydrate Active enZymes (CAZymes) in plants are responsible for the synthesis, modification and degradation of carbohydrate biopolymers, the differences in gene copy number and regulation between woody and herbaceous species have been highlighted previously. There are still many unanswered questions about the role of CAZymes in land plant evolution and the formation of wood, a strong carbohydrate sink.

Results

Here, twenty-two publically available plant genomes were used to characterize the frequency, diversity and complexity of CAZymes in plants. We find that a conserved suite of CAZymes is a feature of land plant evolution, with similar diversity and complexity regardless of growth habit and form. In addition, we compared the diversity and levels of CAZyme gene expression during wood formation in trees using mRNA-seq data from two distantly related angiosperm tree species Eucalyptus grandis and Populus trichocarpa, highlighting the major CAZyme classes involved in xylogenesis and lignocellulosic biomass production.

Conclusions

CAZyme domain ratio across embryophytes is maintained, and the diversity of CAZyme domains is similar in all land plants, regardless of woody habit. The stoichiometric conservation of gene expression in woody and non-woody tissues of Eucalyptus and Populus are indicative of gene balance preservation.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-015-1571-8) contains supplementary material, which is available to authorized users.