Exploring breast tissue microbial composition and the association with breast cancer risk factors.

BACKGROUND: Microbial dysbiosis has emerged as an important element in the development and progression of various cancers, including breast cancer. However, the microbial composition of the breast from healthy individuals, even relative to risk of developing breast cancer, remains unclear. Here, we performed a comprehensive analysis of the microbiota of the normal breast tissue, which was analyzed in relation to the microbial composition of the tumor and adjacent normal tissue. METHODS: The study cohorts included 403 cancer-free women (who donated normal breast tissue cores) and 76 breast cancer patients (who donated tumor and/or adjacent normal tissue samples). Microbiome profiling was obtained by sequencing the nine hypervariable regions of the 16S rRNA gene (V1V2, V2V3, V3V4, V4V5, V5V7, and V7V9). Transcriptome analysis was also performed on 190 normal breast tissue samples. Breast cancer risk score was assessed using the Tyrer-Cuzick risk model. RESULTS: The V1V2 amplicon sequencing resulted more suitable for the analysis of the normal breast microbiome and identified Lactobacillaceae (Firmicutes phylum), Acetobacterraceae, and Xanthomonadaceae (both Proteobacteria phylum) as the most abundant families in the normal breast. However, Ralstonia (Proteobacteria phylum) was more abundant in both breast tumors and histologically normal tissues adjacent to malignant tumors. We also conducted a correlation analysis between the microbiome and known breast cancer risk factors. Abundances of the bacterial taxa Acetotobacter aceti, Lactobacillus vini, Lactobacillus paracasei, and Xanthonomas sp. were associated with age (p < 0.0001), racial background (p < 0.0001), and parity (p < 0.0001). Finally, transcriptome analysis of normal breast tissues showed an enrichment in metabolism- and immune-related genes in the tissues with abundant Acetotobacter aceti, Lactobacillus vini, Lactobacillus paracasei, and Xanthonomas sp., whereas the presence of Ralstonia in the normal tissue was linked to dysregulation of genes involved in the carbohydrate metabolic pathway. CONCLUSIONS: This study defines the microbial features of normal breast tissue, thus providing a basis to understand cancer-related dysbiosis. Moreover, the findings reveal that lifestyle factors can significantly affect the normal breast microbial composition.

High-resolution phylogenetic and population genetic analysis of microbial communities with RoC-ITS.

Microbial communities are inter-connected systems of incredible complexity and dynamism that play crucial roles in health, energy, and the environment. To better understand microbial communities and how they respond to change, it is important to know which microbes are present and their relative abundances at the greatest taxonomic resolution possible. Here, we describe a novel protocol (RoC-ITS) that uses the single-molecule Nanopore sequencing platform to assay the composition of microbial communities at the subspecies designation. Using rolling-circle amplification, this methodology produces long-read sequences from a circular construct containing the complete 16S ribosomal gene and the neighboring internally transcribed spacer (ITS). These long reads can be used to generate a high-fidelity circular consensus sequence. Generally, the ribosomal 16S gene provides phylogenetic information down to the species-level, while the much less conserved ITS region contains strain-level information. When linked together, this combination of markers allows for the identification of individual ribosomal units within a specific organism and the assessment of their relative stoichiometry, as well as the ability to monitor subtle shifts in microbial community composition with a single generic assay. We applied RoC-ITS to an artificial microbial community that was also sequenced using the Illumina platform, to assess its accuracy in quantifying the relative abundance and identity of each species.

The Polar Flagellar Transcriptional Regulatory Network in Vibrio campbellii Deviates from Canonical Vibrio Species.

Swimming motility is a critical virulence factor in pathogenesis for numerous Vibrio species. Vibrio campbellii DS40M4 is a wild-type isolate that has been recently established as a highly tractable model strain for bacterial genetics studies. We sought to exploit the tractability and relevance of this strain for characterization of flagellar gene regulation in V. campbellii. Using comparative genomics, we identified homologs of V. campbellii flagellar and chemotaxis genes conserved in other members of the Vibrionaceae and determined the transcriptional profile of these loci using differential RNA-seq. We systematically deleted all 63 predicted flagellar and chemotaxis genes in V. campbellii and examined their effects on motility and flagellum production. We specifically focused on the core regulators of the flagellar hierarchy established in other vibrios: RpoN (σ54), FlrA, FlrC, and FliA. Our results show that V. campbellii transcription of flagellar and chemotaxis genes is governed by a multitiered regulatory hierarchy similar to other motile Vibrio species. However, there are several critical differences in V. campbellii: (i) the σ54-dependent regulator FlrA is dispensable for motility; (ii) the flgA, fliEFGHIJ, flrA, and flrBC operons do not require σ54 for expression; and (iii) FlrA and FlrC coregulate class II genes. Our model proposes that the V. campbellii flagellar transcriptional hierarchy has three classes of genes, in contrast to the four-class hierarchy in Vibrio cholerae. Our genetic and phenotypic dissection of the V. campbellii flagellar regulatory network highlights the differences that have evolved in flagellar regulation across the Vibrionaceae. IMPORTANCE Vibrio campbellii is a Gram-negative bacterium that is free-living and ubiquitous in marine environments and is an important global pathogen of fish and shellfish. Disruption of the flagellar motor significantly decreases host mortality of V. campbellii, suggesting that motility is a key factor in pathogenesis. Using this model organism, we identified >60 genes that encode proteins with predicted structural, mechanical, or regulatory roles in function of the single polar flagellum in V. campbellii. We systematically tested strains containing single deletions of each gene to determine the impact on motility and flagellum production. Our studies have uncovered differences in the regulatory network and function of several genes in V. campbellii compared to established systems in Vibrio cholerae and Vibrio parahaemolyticus.

Probing biophysical sequence constraints within the transmembrane domains of rhodopsin by deep mutational scanning.

Membrane proteins must balance the sequence constraints associated with folding and function against the hydrophobicity required for solvation within the bilayer. We recently found the expression and maturation of rhodopsin are limited by the hydrophobicity of its seventh transmembrane domain (TM7), which contains polar residues that are essential for function. On the basis of these observations, we hypothesized that rhodopsin's expression should be less tolerant of mutations in TM7 relative to those within hydrophobic TM domains. To test this hypothesis, we used deep mutational scanning to compare the effects of 808 missense mutations on the plasma membrane expression of rhodopsin in HEK293T cells. Our results confirm that a higher proportion of mutations within TM7 (37%) decrease rhodopsin's plasma membrane expression relative to those within a hydrophobic TM domain (TM2, 25%). These results in conjunction with an evolutionary analysis suggest solvation energetics likely restricts the evolutionary sequence space of polar TM domains.

Unanticipated Geochemical and Microbial Community Structure under Seasonal Ice Cover in a Dilute, Dimictic Arctic Lake.

Despite most lakes in the Arctic being perennially or seasonally frozen for at least 40% of the year, little is known about microbial communities and nutrient cycling under ice cover. We assessed the vertical microbial community distribution and geochemical composition in early spring under ice in a seasonally ice-covered lake in southwest Greenland using amplicon-based sequencing that targeted 16S rRNA genes and using a combination of field and laboratory aqueous geochemical methods. Microbial communities changed consistently with changes in geochemistry. Composition of the abundant members responded strongly to redox conditions, shifting downward from a predominantly heterotrophic aerobic community in the suboxic waters to a heterotrophic anaerobic community in the anoxic waters. Operational taxonomic units (OTUs) of Sporichthyaceae, Comamonadaceae, and the SAR11 Clade had higher relative abundances above the oxycline and OTUs within the genus Methylobacter, the phylum Lentisphaerae, and purple sulfur bacteria (PSB) below the oxycline. Notably, a 13-fold increase in sulfide at the oxycline was reflected in an increase and change in community composition of potential sulfur oxidizers. Purple non-sulfur bacteria were present above the oxycline and green sulfur bacteria and PSB coexisted below the oxycline, however, PSB were most abundant. For the first time we show the importance of PSB as potential sulfur oxidizers in an Arctic dimictic lake.

Two temporal functions of Glass: Ommatidium patterning and photoreceptor differentiation.

Much progress has been made in elucidating the molecular networks required for specifying retinal cells, including photoreceptors, but the downstream mechanisms that maintain identity and regulate differentiation remain poorly understood. Here, we report that the transcription factor Glass has a dual role in establishing a functional Drosophila eye. Utilizing conditional rescue approaches, we confirm that persistent defects in ommatidium patterning combined with cell death correlate with the overall disruption of eye morphology in glass mutants. In addition, we reveal that Glass exhibits a separable role in regulating photoreceptor differentiation. In particular, we demonstrate the apparent loss of glass mutant photoreceptors is not only due to cell death but also a failure of the surviving photoreceptors to complete differentiation. Moreover, the late reintroduction of Glass in these developmentally stalled photoreceptors is capable of restoring differentiation in the absence of correct ommatidium patterning. Mechanistically, transcription profiling at the time of differentiation reveals that Glass is necessary for the expression of many genes implicated in differentiation, i.e. rhabdomere morphogenesis, phototransduction, and synaptogenesis. Specifically, we show Glass directly regulates the expression of Pph13, which encodes a transcription factor necessary for opsin expression and rhabdomere morphogenesis. Finally, we demonstrate the ability of Glass to choreograph photoreceptor differentiation is conserved between Drosophila and Tribolium, two holometabolous insects. Altogether, our work identifies a fundamental regulatory mechanism to generate the full complement of cells required for a functional rhabdomeric visual system and provides a critical framework to investigate the basis of differentiation and maintenance of photoreceptor identity.

Combined Flux Chamber and Genomics Approach Links Nitrous Acid Emissions to Ammonia Oxidizing Bacteria and Archaea in Urban and Agricultural Soil.

Nitrous acid (HONO) is a photochemical source of hydroxyl radical and nitric oxide in the atmosphere that stems from abiotic and biogenic processes, including the activity of ammonia-oxidizing soil microbes. HONO fluxes were measured from agricultural and urban soil in mesocosm studies aimed at characterizing biogenic sources and linking them to indigenous microbial consortia. Fluxes of HONO from agricultural and urban soil were suppressed by addition of a nitrification inhibitor and enhanced by amendment with ammonium (NH4(+)), with peaks at 19 and 8 ng m(-2) s(-1), respectively. In addition, both agricultural and urban soils were observed to convert (15)NH4(+) to HO(15)NO. Genomic surveys of soil samples revealed that 1.5-6% of total expressed 16S rRNA sequences detected belonged to known ammonia oxidizing bacteria and archaea. Peak fluxes of HONO were directly related to the abundance of ammonia-oxidizer sequences, which in turn depended on soil pH. Peak HONO fluxes under fertilized conditions are comparable in magnitude to fluxes reported during field campaigns. The results suggest that biogenic HONO emissions will be important in soil environments that exhibit high nitrification rates (e.g., agricultural soil) although the widespread occurrence of ammonia oxidizers implies that biogenic HONO emissions are also possible in the urban and remote environment.

Spatial soil heterogeneity has a greater effect on symbiotic arbuscular mycorrhizal fungal communities and plant growth than genetic modification with Bacillus thuringiensis toxin genes.

Maize, genetically modified with the insect toxin genes of Bacillus thuringiensis (Bt), is widely cultivated, yet its impacts on soil organisms are poorly understood. Arbuscular mycorrhizal fungi (AMF) form symbiotic associations with plant roots and may be uniquely sensitive to genetic changes within a plant host. In this field study, the effects of nine different lines of Bt maize and their corresponding non-Bt parental isolines were evaluated on AMF colonization and community diversity in plant roots. Plants were harvested 60 days after sowing, and data were collected on plant growth and per cent AMF colonization of roots. AMF community composition in roots was assessed using 454 pyrosequencing of the 28S rRNA genes, and spatial variation in mycorrhizal communities within replicated experimental field plots was examined. Growth responses, per cent AMF colonization of roots and AMF community diversity in roots did not differ between Bt and non-Bt maize, but root and shoot biomass and per cent colonization by arbuscules varied by maize cultivar. Plot identity had the most significant effect on plant growth, AMF colonization and AMF community composition in roots, indicating spatial heterogeneity in the field. Mycorrhizal fungal communities in maize roots were autocorrelated within approximately 1 m, but at greater distances, AMF community composition of roots differed between plants. Our findings indicate that spatial variation and heterogeneity in the field has a greater effect on the structure of AMF communities than host plant cultivar or modification by Bt toxin genes.

Concordance of bacterial communities of two tick species and blood of their shared rodent host.

High-throughput sequencing is revealing that most macro-organisms house diverse microbial communities. Of particular interest are disease vectors whose microbiome could potentially affect pathogen transmission and vector competence. We investigated bacterial community composition and diversity of the ticks Dermacentor variabilis (n = 68) and Ixodes scapularis (n = 15) and blood of their shared rodent host, Peromyscus leucopus (n = 45) to quantify bacterial diversity and concordance. The 16S rRNA gene was amplified from genomic DNA from field-collected tick and rodent blood samples, and 454 pyrosequencing was used to elucidate their bacterial communities. After quality control, over 300 000 sequences were obtained and classified into 118 operational taxonomic units (OTUs, clustered at 97% similarity). Analysis of rarefied communities revealed that the most abundant OTUs were tick species-specific endosymbionts, Francisella and Rickettsia, and the commonly flea-associated bacterium Bartonella in rodent blood. An Arsenophonus and additional Francisella endosymbiont were also present in D. variabilis samples. Rickettsia was found in both tick species but not in rodent blood, suggesting that it is not transmitted during feeding. Bartonella was present in larvae and nymphs of both tick species, even those scored as unengorged. Relatively, few OTUs (e.g. Bartonella, Lactobacillus) were found in all sample types. Overall, bacterial communities from each sample type were significantly different and highly structured, independent of their dominant OTUs. Our results point to complex microbial assemblages inhabiting ticks and host blood including infectious agents, tick-specific endosymbionts and environmental bacteria that could potentially affect arthropod-vectored disease dynamics.

Genome sequences of eight morphologically diverse Alphaproteobacteria.

The Alphaproteobacteria comprise morphologically diverse bacteria, including many species of stalked bacteria. Here we announce the genome sequences of eight alphaproteobacteria, including the first genome sequences of species belonging to the genera Asticcacaulis, Hirschia, Hyphomicrobium, and Rhodomicrobium.

An Ergatis-based prokaryotic genome annotation web server.

SUMMARY: Ergatis is a flexible workflow management system for designing and executing complex bioinformatics pipelines. However, its complexity restricts its usage to only highly skilled bioinformaticians. We have developed a web-based prokaryotic genome annotation server, Integrative Services for Genomics Analysis (ISGA), which builds upon the Ergatis workflow system, integrates other dynamic analysis tools and provides intuitive web interfaces for biologists to customize and execute their own annotation pipelines. ISGA is designed to be installed at genomics core facilities and be used directly by biologists. AVAILABILITY: ISGA is accessible at http://isga.cgb.indiana.edu/ and the system is also freely available for local installation.

ESTPiper--a web-based analysis pipeline for expressed sequence tags.

BACKGROUND: EST sequencing projects are increasing in scale and scope as the genome sequencing technologies migrate from core sequencing centers to individual research laboratories. Effectively, generating EST data is no longer a bottleneck for investigators. However, processing large amounts of EST data remains a non-trivial challenge for many. Web-based EST analysis tools are proving to be the most convenient option for biologists when performing their analysis, so these tools must continuously improve on their utility to keep in step with the growing needs of research communities. We have developed a web-based EST analysis pipeline called ESTPiper, which streamlines typical large-scale EST analysis components. RESULTS: The intuitive web interface guides users through each step of base calling, data cleaning, assembly, genome alignment, annotation, analysis of gene ontology (GO), and microarray oligonucleotide probe design. Each step is modularized. Therefore, a user can execute them separately or together in batch mode. In addition, the user has control over the parameters used by the underlying programs. Extensive documentation of ESTPiper's functionality is embedded throughout the web site to facilitate understanding of the required input and interpretation of the computational results. The user can also download intermediate results and port files to separate programs for further analysis. In addition, our server provides a time-stamped description of the run history for reproducibility. The pipeline can also be installed locally, allowing researchers to modify ESTPiper to suit their own needs. CONCLUSION: ESTPiper streamlines the typical process of EST analysis. The pipeline was initially designed in part to support the Daphnia pulex cDNA sequencing project. A web server hosting ESTPiper is provided at http://estpiper.cgb.indiana.edu/ to now support projects of all size. The software is also freely available from the authors for local installations.

BOV--a web-based BLAST output visualization tool.

BACKGROUND: The BLAST program is one of the most widely used sequence similarity search tools for genomic research, even by those biologists lacking extensive bioinformatics training. As the availability of sequence data increases, more researchers are downloading the BLAST program for local installation and performing larger and more complex tasks, including batch queries. In order to manage and interpret the results of batch queries, a host of software packages have been developed to assist with data management and post-processing. Among these programs, there is almost a complete lack of visualization tools to provide graphic representation of complex BLAST pair-wise alignments. We have developed a web-based program, BLAST Output Visualization Tool (BOV), that allows users to interactively visualize the matching regions of query and database hit sequences, thereby allowing the user to quickly and easily dissect complex matching patterns. RESULTS: Users can upload the standard BLAST output in pair-wise alignment format as input to the web server (including batch queries generated installing and running the stand-alone BLAST program on a local server). The program extracts the alignment coordinates of matching regions between the query and the corresponding database hit sequence. The coordinates are used to plot each matching region as colored lines or trapezoids. Using the straightforward control panels throughout the web site, each plotted matching region can be easily explored in detail by, for example, highlighting the region of interest or examining the raw pair-wise sequence alignment. Tutorials are provided at the website to guide users step-by-step through the functional features of BOV. CONCLUSION: BOV provides a user-friendly web interface to visualize the standard BLAST output for investigating wide-ranging genomic problems, including single query and batch query datasets. In particular, this software is valuable to users interested in identifying regions of co-linearity, duplication, translocation, and inversion among sequences. A web server hosting BOV is accessible via http://bioportal.cgb.indiana.edu/cgi-bin/BOV/index.cgi and the software is freely available for local installations.

A study of residue correlation within protein sequences and its application to sequence classification.

We investigate methods of estimating residue correlation within protein sequences. We begin by using mutual information (MI) of adjacent residues, and improve our methodology by defining the mutual information vector (MIV) to estimate long range correlations between nonadjacent residues. We also consider correlation based on residue hydropathy rather than protein-specific interactions. Finally, in experiments of family classification tests, the modeling power of MIV was shown to be significantly better than the classic MI method, reaching the level where proteins can be classified without alignment information.