Multi-Omics Characterization of Host-Derived Bacillus spp. Probiotics for Improved Growth Performance in Poultry

Microbial feed ingredients or probiotics have been used widely in the poultry industry to improve production efficiency. Spore-forming Bacillus spp. offer advantages over traditional probiotic strains as Bacillus spores are resilient to high temperature, acidic pH, and desiccation. This results in increased strain viability during manufacturing and feed-pelleting processes, extended product shelf-life, and increased stability within the animal’s gastrointestinal tract. Despite numerous reports on the use of Bacillus spores as feed additives, detailed characterizations of Bacillus probiotic strains are typically not published. Insufficient characterizations can lead to misidentification of probiotic strains in product labels, and the potential application of strains carrying virulence factors, toxins, antibiotic resistance, or toxic metabolites. Hence, it is critical to characterize in detail the genomic and phenotypic properties of these strains to screen out undesirable properties and to tie individual traits to clinical outcomes and possible mechanisms. Here, we report a screening workflow and comprehensive multi-omics characterization of Bacillus spp. for use in broiler chickens. Host-derived Bacillus strains were isolated and screened for desirable probiotic properties. The phenotypic, genomic and metabolomic analyses of three probiotic candidates, two Bacillus amyloliquefaciens (Ba ATCC PTA126784 and ATCC PTA126785), and a Bacillus subtilis (Bs ATCC PTA126786), showed that all three strains had promising probiotic traits and safety profiles. Inclusion of Ba ATCC PTA12684 (Ba-PTA84) in the feed of broiler chickens resulted in improved growth performance, as shown by a significantly improved feed conversion ratio (3.3%), increased of European Broiler Index (6.2%), and increased average daily gain (ADG) (3.5%). Comparison of the cecal microbiomes from Ba PTA84-treated and control animals suggested minimal differences in microbiome structure, indicating that the observed growth promotion presumably was not mediated by modulation of cecal microbiome.

Genome analysis of a major urban malaria vector mosquito, Anopheles stephensi

Background

Anopheles stephensi is the key vector of malaria throughout the Indian subcontinent and Middle East and an emerging model for molecular and genetic studies of mosquito-parasite interactions. The type form of the species is responsible for the majority of urban malaria transmission across its range.

Results

Here, we report the genome sequence and annotation of the Indian strain of the type form of An. stephensi. The 221 Mb genome assembly represents more than 92% of the entire genome and was produced using a combination of 454, Illumina, and PacBio sequencing. Physical mapping assigned 62% of the genome onto chromosomes, enabling chromosome-based analysis. Comparisons between An. stephensi and An. gambiae reveal that the rate of gene order reshuffling on the X chromosome was three times higher than that on the autosomes. An. stephensi has more heterochromatin in pericentric regions but less repetitive DNA in chromosome arms than An. gambiae. We also identify a number of Y-chromosome contigs and BACs. Interspersed repeats constitute 7.1% of the assembled genome while LTR retrotransposons alone comprise more than 49% of the Y contigs. RNA-seq analyses provide new insights into mosquito innate immunity, development, and sexual dimorphism.

Conclusions

The genome analysis described in this manuscript provides a resource and platform for fundamental and translational research into a major urban malaria vector. Chromosome-based investigations provide unique perspectives on Anopheles chromosome evolution. RNA-seq analysis and studies of immunity genes offer new insights into mosquito biology and mosquito-parasite interactions.

Electronic supplementary material

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

Comparative phylogenomics and evolution of the Brucellae reveal a path to virulence

Brucella species include important zoonotic pathogens that have a substantial impact on both agriculture and human health throughout the world. Brucellae are thought of as "stealth pathogens" that escape recognition by the host innate immune response, modulate the acquired immune response, and evade intracellular destruction. We analyzed the genome sequences of members of the family Brucellaceae to assess its evolutionary history from likely free-living soil-based progenitors into highly successful intracellular pathogens. Phylogenetic analysis split the genus into two groups: recently identified and early-dividing "atypical" strains and a highly conserved "classical" core clade containing the major pathogenic species. Lateral gene transfer events brought unique genomic regions into Brucella that differentiated them from Ochrobactrum and allowed the stepwise acquisition of virulence factors that include a type IV secretion system, a perosamine-based O antigen, and systems for sequestering metal ions that are absent in progenitors. Subsequent radiation within the core Brucella resulted in lineages that appear to have evolved within their preferred mammalian hosts, restricting their virulence to become stealth pathogens capable of causing long-term chronic infections.

Phylogeographic evidence of cognate recognition site patterns and transformation efficiency differences in H. pylori: theory of strain dominance

BACKGROUND

Helicobacter pylori has diverged in parallel to its human host, leading to distinct phylogeographic populations. Recent evidence suggests that in the current human mixing in Latin America, European H. pylori (hpEurope) are increasingly dominant at the expense of Amerindian haplotypes (hspAmerind). This phenomenon might occur via DNA recombination, modulated by restriction-modification systems (RMS), in which differences in cognate recognition sites (CRS) and in active methylases will determine direction and frequency of gene flow. We hypothesized that genomes from hspAmerind strains that evolved from a small founder population have lost CRS for RMS and active methylases, promoting hpEurope's DNA invasion. We determined the observed and expected frequencies of CRS for RMS in DNA from 7 H. pylori whole genomes and 110 multilocus sequences. We also measured the number of active methylases by resistance to in vitro digestion by 16 restriction enzymes of genomic DNA from 9 hpEurope and 9 hspAmerind strains, and determined the direction of DNA uptake in co-culture experiments of hspAmerind and hpEurope strains.

RESULTS

Most of the CRS were underrepresented with consistency between whole genomes and multilocus sequences. Although neither the frequency of CRS nor the number of active methylases differ among the bacterial populations (average 8.6 ± 2.6), hspAmerind strains had a restriction profile distinct from that in hpEurope strains, with 15 recognition sites accounting for the differences. Amerindians strains also exhibited higher transformation rates than European strains, and were more susceptible to be subverted by larger DNA hpEurope-fragments than vice versa.

CONCLUSIONS

The geographical variation in the pattern of CRS provides evidence for ancestral differences in RMS representation and function, and the transformation findings support the hypothesis of Europeanization of the Amerindian strains in Latin America via DNA recombination.

Contrasting patterns of evolution following whole genome versus tandem duplication events in Populus

Comparative analysis of multiple angiosperm genomes has implicated gene duplication in the expansion and diversification of many gene families. However, empirical data and theory suggest that whole-genome and small-scale duplication events differ with respect to the types of genes preserved as duplicate pairs. We compared gene duplicates resulting from a recent whole genome duplication to a set of tandemly duplicated genes in the model forest tree Populus trichocarpa. We used a combination of microarray expression analyses of a diverse set of tissues and functional annotation to assess factors related to the preservation of duplicate genes of both types. Whole genome duplicates are 700 bp longer and are expressed in 20% more tissues than tandem duplicates. Furthermore, certain functional categories are over-represented in each class of duplicates. In particular, disease resistance genes and receptor-like kinases commonly occur in tandem but are significantly under-retained following whole genome duplication, while whole genome duplicate pairs are enriched for members of signal transduction cascades and transcription factors. The shape of the distribution of expression divergence for duplicated pairs suggests that nearly half of the whole genome duplicates have diverged in expression by a random degeneration process. The remaining pairs have more conserved gene expression than expected by chance, consistent with a role for selection under the constraints of gene balance. We hypothesize that duplicate gene preservation in Populus is driven by a combination of subfunctionalization of duplicate pairs and purifying selection favoring retention of genes encoding proteins with large numbers of interactions.

Microbial diversity in saliva of oral squamous cell carcinoma

In the oral cavity, chronic inflammation has been observed at various stages of oral squamous cell carcinomas (OSCC). Such inflammation could result from persistent mucosal or epithelial cell colonization by microorganisms. There is increasing evidence of the involvement of oral bacteria in inflammation, warranting further studies on the association of bacteria with the progression of OSCC. The objective of this study was to evaluate the diversity and relative abundance of bacteria in the saliva of subjects with OSCC. Using 454 parallel DNA sequencing, ∼58,000 PCR amplicons that span the V4-V5 hypervariable region of rRNAs from five subjects were sequenced. Members of eight phyla (divisions) of bacteria were detected. The majority of classified sequences belonged to the phyla Firmicutes (45%) and Bacteroidetes (25%). Further, 52 different genera containing approximately 860 (16.51%) known species were identified and 1077 (67%) sequences belonging to various uncultured bacteria or unclassified groups. The species diversity estimates obtained with abundance-based coverage estimators and Chao1 were greater than published analyses of other microbial profiles from the oral cavity. Fifteen unique phylotypes were present in all three OSCC subjects.

Transcriptome sequencing of the Microarray Quality Control (MAQC) RNA reference samples using next generation sequencing

Background

Transcriptome sequencing using next-generation sequencing platforms will soon be competing with DNA microarray technologies for global gene expression analysis. As a preliminary evaluation of these promising technologies, we performed deep sequencing of cDNA synthesized from the Microarray Quality Control (MAQC) reference RNA samples using Roche's 454 Genome Sequencer FLX.

Results

We generated more that 3.6 million sequence reads of average length 250 bp for the MAQC A and B samples and introduced a data analysis pipeline for translating cDNA read counts into gene expression levels. Using BLAST, 90% of the reads mapped to the human genome and 64% of the reads mapped to the RefSeq database of well annotated genes with e-values ≤ 10^-20. We measured gene expression levels in the A and B samples by counting the numbers of reads that mapped to individual RefSeq genes in multiple sequencing runs to evaluate the MAQC quality metrics for reproducibility, sensitivity, specificity, and accuracy and compared the results with DNA microarrays and Quantitative RT-PCR (QRTPCR) from the MAQC studies. In addition, 88% of the reads were successfully aligned directly to the human genome using the AceView alignment programs with an average 90% sequence similarity to identify 137,899 unique exon junctions, including 22,193 new exon junctions not yet contained in the RefSeq database.

Conclusion

Using the MAQC metrics for evaluating the performance of gene expression platforms, the ExpressSeq results for gene expression levels showed excellent reproducibility, sensitivity, and specificity that improved systematically with increasing shotgun sequencing depth, and quantitative accuracy that was comparable to DNA microarrays and QRTPCR. In addition, a careful mapping of the reads to the genome using the AceView alignment programs shed new light on the complexity of the human transcriptome including the discovery of thousands of new splice variants.

Molecular and physiological adaptation to prolonged drought stress in the leaves of two Andean potato genotypes

Responses to prolonged drought and recovery from drought of two South American potato (Solanum tuberosum L. ssp. andigena (Juz & Buk) Hawkes) landraces, Sullu and Ccompis were compared under field conditions. Physiological and biomass measurements, yield analysis, the results of hybridisation to a potato microarray platform (44 000 probes) and metabolite profiling were used to characterise responses to water deficit. Drought affected shoot and root biomass negatively in Ccompis but not in Sullu, whereas both genotypes maintained tuber yield under water stress. Ccompis showed stronger reduction in maximum quantum yield under stress than Sullu, and less decrease in stomatal resistance. Genes associated with PSII functions were activated during recovery in Sullu only. Evidence for sucrose accumulation in Sullu only during maximum stress and recovery was observed, in addition to increases in cell wall biosynthesis. A depression in the abundance of plastid superoxide dismutase transcripts was observed under maximum stress in Ccompis. Both sucrose and the regulatory molecule trehalose accumulated in the leaves of Sullu only. In contrast, in Ccompis, the raffinose oligosaccharide family pathway was activated, whereas low levels of sucrose and minor stress-mediated changes in trehalose were observed. Proline, and expression of the associated genes, rose in both genotypes under drought, with a 3-fold higher increase in Sullu than in Ccompis. The results demonstrate the presence of distinct molecular and biochemical drought responses in the two potato landraces leading to yield maintenance but differential biomass accumulation in vegetative tissues.

Physiological and molecular adaptations to drought in Andean potato genotypes

The drought stress tolerance of two Solanum tuberosum subsp. andigena landraces, one hybrid (adgxtbr) and Atlantic (S. tuberosum subsp. tuberosum) has been evaluated. Photosynthesis in the Andigena landraces during prolonged drought was maintained significantly longer than in the Tuberosum (Atlantic) line. Among the Andigena landraces, 'Sullu' (SUL) was more drought resistant than 'Negra Ojosa' (NOJ). Microarray analysis and metabolite data from leaf samples taken at the point of maximum stress suggested higher mitochondrial metabolic activity in SUL than in NOJ. A greater induction of chloroplast-localized antioxidant and chaperone genes in SUL compared with NOJ was evident. ABA-responsive TFs were more induced in NOJ compared with SUL, including WRKY1, mediating a response in SA signalling that may give rise to increased ROS. NOJ may be experiencing higher ROS levels than SUL. Metabolite profiles of NOJ were characterized by compounds indicative of stress, for example, proline, trehalose, and GABA, which accumulated to a higher degree than in SUL. The differences between the Andigena lines were not explained by protective roles of compatible solutes; hexoses and complex sugars were similar in both landraces. Instead, lower levels of ROS accumulation, greater mitochondrial activity and active chloroplast defences contributed to a lower stress load in SUL than in NOJ during drought.

Early PLDalpha-mediated events in response to progressive drought stress in Arabidopsis: a transcriptome analysis

Phospholipase D (PLD) has been implicated in a variety of stresses including osmotic stress and wounding. PLDalpha1-derived phosphatidic acid interacts with ABI1 phosphatase 2C and promotes abscisic acid signalling. It has also been shown to regulate proline biosynthesis negatively. Plants with abrogated PLDalpha show insensitivity to abscisic acid (ABA) and impaired stomatal conductance. The goal in the present study was to identify early PLDalpha-mediated events in response to progressive drought stress in Arabidopsis. Water was withheld from 7-week-old Arabidopsis thaliana (Col-0) and antisense-PLDalpha1 (anti-PLDalpha) in a controlled environment chamber. Diurnal leaf water potential (LWP) and photosynthesis measurements were recorded five and three times a day, respectively. Quantitative reverse transcription-polymerase chain reaction (qRT-PCR) and microarray analyses were conducted using RNA from shoots collected at the fourth LWP time point on the ninth day after stress imposition. Anti-PLDalpha experienced severe water stress (-1.28 MPa) at the same time period that Col-0 experienced less water stress (-0.31 MPa). Diurnal LWP measurements showed that anti-PLDalpha had a lower LWP than Col-0 in both control and drought-stress conditions. Photosynthesis was also more affected in anti-PLDalpha than in Col-0. Anti-PLDalpha plants recovered fully following rehydration after 10 d of stress. qRT-PCR revealed up to 18-fold lower values for PLDalpha transcripts in stressed anti-PLDalpha plants when compared with stressed Col-0. Microarray expression profiles revealed distinct gene expression patterns in Col-0 and anti-PLDalpha. No differences in gene expression were detected between the two genotypes in the absence of drought stress. ROP8, PLDdelta, and lipid transfer proteins were among the differentially expressed genes between the two genotypes.

PATRIC: the VBI PathoSystems Resource Integration Center

The PathoSystems Resource Integration Center (PATRIC) is one of eight Bioinformatics Resource Centers (BRCs) funded by the National Institute of Allergy and Infection Diseases (NIAID) to create a data and analysis resource for selected NIAID priority pathogens, specifically proteobacteria of the genera Brucella, Rickettsia and Coxiella, and corona-, calici- and lyssaviruses and viruses associated with hepatitis A and E. The goal of the project is to provide a comprehensive bioinformatics resource for these pathogens, including consistently annotated genome, proteome and metabolic pathway data to facilitate research into counter-measures, including drugs, vaccines and diagnostics. The project's curation strategy has three prongs: ‘breadth first’ beginning with whole-genome and proteome curation using standardized protocols, a ‘targeted’ approach addressing the specific needs of researchers and an integrative strategy to leverage high-throughput experimental data (e.g. microarrays, proteomics) and literature. The PATRIC infrastructure consists of a relational database, analytical pipelines and a website which supports browsing, querying, data visualization and the ability to download raw and curated data in standard formats. At present, the site warehouses complete sequences for 17 bacterial and 332 viral genomes. The PATRIC website () will continually grow with the addition of data, analysis and functionality over the course of the project.

Response diversity of Arabidopsis thaliana ecotypes in elevated [CO2] in the field

Free Air [CO(2)] Enrichment (FACE) allows for plant growth under fully open-air conditions of elevated [CO(2)] at concentrations expected to be reached by mid-century. We used Arabidopsis thaliana ecotypes Col-0, Cvi-0, and WS to analyze changes in gene expression and metabolite profiles of plants grown in "SoyFACE" (http://www.soyface.uiuc.edu/), a system of open-air rings within which [CO(2)] is elevated to approximately 550 ppm. Data from multiple rings, comparing plants in ambient air and elevated [CO(2)], were analyzed by mixed model ANOVA, linear discriminant analysis (LDA) and data-mining tools. In elevated [CO(2)], decreases in the expression of genes related to chloroplast functions characterized all lines but individual members of distinct multi-gene families were regulated differently between lines. Also, different strategies distinguished the lines with respect to the regulation of genes related to carbohydrate biosynthesis and partitioning, N-allocation and amino acid metabolism, cell wall biosynthesis, and hormone responses, irrespective of the plants' developmental status. Metabolite results paralleled reactions seen at the level of transcript expression. Evolutionary adaptation of species to their habitat and intrinsic genetic plasticity seem to determine the nature of responses to elevated [CO(2)]. Irrespective of their underlying genetic diversity, and evolutionary adaptation to different habitats, a small number of common, predominantly stress-responsive, signature transcripts appear to characterize responses of the Arabidopsis ecotypes in FACE.

Effects of chronic ozone exposure on gene expression in Arabidopsis thaliana ecotypes and in Thellungiella halophila

Arabidopsis thaliana (At) ecotypes Columbia-0 (Col-0), Wassilewskija (WS), Cape Verde Islands (Cvi-0) and a relative, Thellungiella halophila (Th), were exposed to 20-25% over ambient ozone [O3] in a free air concentration enrichment (FACE) experiment (http://www.soyFACE. uiuc.edu), mirroring increases expected in the near future. Col-0 and WS accelerated development and developed lesions within 10 d under increased ozone, while Cvi-0 and Th grew slowly. RNAs were used in microarray hybridizations (Col-0-based 26 000 elements, 70-mer oligonucleotides). A two-step analysis of variance (ANOVA) model, including comparison with values obtained under [O3], was used for analyses. WS showed the greatest number of changes in gene expression in response to ozone. Th showed the least changes, suggesting that its expression state at [O3] was sufficient for resistance at increased ozone. Patterns observed in ambient air controls for Cvi-0 and Col-0 were most similar, while Th showed the greatest number of differences compared with the other controls. Compared with Col-0, however, Cvi-0 showed higher levels of expression of chaperones, receptor kinase-like and photosynthesis-related genes in ambient air. Cvi-0 exhibited ozone-mediated changes in a pathway involving AtSR, a homologue of the mammalian NF kappa B family of redox-sensitive transcription factors, changes in chaperones, WRKY and C2H2 proteins and antioxidants. WS displayed ozone-mediated decreases in the expression of two AtSR/NF kappa B family members, C2-domain proteins and genes associated with cell wall growth and changes in the expression of marker genes for programmed cell death (PCD), among them RCD1, a key regulator in this pathway. Microarray data were verified by reverse transcriptase (RT)-PCR. We relate O3-response diversity across the four lines to different responses among signaling and transcriptional response networks and differences in gene expression at [O3] levels.

The statistics of identifying differentially expressed genes in Expresso and TM4: a comparison

Background

Analysis of DNA microarray data takes as input spot intensity measurements from scanner software and returns differential expression of genes between two conditions, together with a statistical significance assessment. This process typically consists of two steps: data normalization and identification of differentially expressed genes through statistical analysis. The Expresso microarray experiment management system implements these steps with a two-stage, log-linear ANOVA mixed model technique, tailored to individual experimental designs. The complement of tools in TM4, on the other hand, is based on a number of preset design choices that limit its flexibility. In the TM4 microarray analysis suite, normalization, filter, and analysis methods form an analysis pipeline. TM4 computes integrated intensity values (IIV) from the average intensities and spot pixel counts returned by the scanner software as input to its normalization steps. By contrast, Expresso can use either IIV data or median intensity values (MIV). Here, we compare Expresso and TM4 analysis of two experiments and assess the results against qRT-PCR data.

Results

The Expresso analysis using MIV data consistently identifies more genes as differentially expressed, when compared to Expresso analysis with IIV data. The typical TM4 normalization and filtering pipeline corrects systematic intensity-specific bias on a per microarray basis. Subsequent statistical analysis with Expresso or a TM4 t-test can effectively identify differentially expressed genes. The best agreement with qRT-PCR data is obtained through the use of Expresso analysis and MIV data.

Conclusion

The results of this research are of practical value to biologists who analyze microarray data sets. The TM4 normalization and filtering pipeline corrects microarray-specific systematic bias and complements the normalization stage in Expresso analysis. The results of Expresso using MIV data have the best agreement with qRT-PCR results. In one experiment, MIV is a better choice than IIV as input to data normalization and statistical analysis methods, as it yields as greater number of statistically significant differentially expressed genes; TM4 does not support the choice of MIV input data. Overall, the more flexible and extensive statistical models of Expresso achieve more accurate analytical results, when judged by the yardstick of qRT-PCR data, in the context of an experimental design of modest complexity.