Susan M. Huse

The evolution of marine microbes over billions of years predicts that the composition of microbial communities should be much greater than the published estimates of a few thousand distinct kinds of microbes per liter of seawater. By adopting a massively parallel tag sequencing strategy, we show that bacterial communities of deep water masses of the North Atlantic and diffuse flow hydrothermal vents are one to two orders of magnitude more complex than previously reported for any microbial environment. A relatively small number of different populations dominate all samples, but thousands of low-abundance populations account for most of the observed phylogenetic diversity. This "rare biosphere" is very ancient and may represent a nearly inexhaustible source of genomic innovation. Members of the rare biosphere are highly divergent from each other and, at different times in earth's history, may have had a profound impact on shaping planetary processes.

The human intestinal microbiota is essential to the health of the host and plays a role in nutrition, development, metabolism, pathogen resistance, and regulation of immune responses. Antibiotics may disrupt these coevolved interactions, leading to acute or chronic disease in some individuals. Our understanding of antibiotic-associated disturbance of the microbiota has been limited by the poor sensitivity, inadequate resolution, and significant cost of current research methods. The use of pyrosequencing technology to generate large numbers of 16S rDNA sequence tags circumvents these limitations and has been shown to reveal previously unexplored aspects of the "rare biosphere." We investigated the distal gut bacterial communities of three healthy humans before and after treatment with ciprofloxacin, obtaining more than 7,000 full-length rRNA sequences and over 900,000 pyrosequencing reads from two hypervariable regions of the rRNA gene. A companion paper in PLoS Genetics (see Huse et al., doi: 10.1371/journal.pgen.1000255) shows that the taxonomic information obtained with these methods is concordant. Pyrosequencing of the V6 and V3 variable regions identified 3,300-5,700 taxa that collectively accounted for over 99% of the variable region sequence tags that could be obtained from these samples. Ciprofloxacin treatment influenced the abundance of about a third of the bacterial taxa in the gut, decreasing the taxonomic richness, diversity, and evenness of the community. However, the magnitude of this effect varied among individuals, and some taxa showed interindividual variation in the response to ciprofloxacin. While differences of community composition between individuals were the largest source of variability between samples, we found that two unrelated individuals shared a surprising degree of community similarity. In all three individuals, the taxonomic composition of the community closely resembled its pretreatment state by 4 weeks after the end of treatment, but several taxa failed to recover within 6 months. These pervasive effects of ciprofloxacin on community composition contrast with the reports by participants of normal intestinal function and with prior assumptions of only modest effects of ciprofloxacin on the intestinal microbiota. These observations support the hypothesis of functional redundancy in the human gut microbiota. The rapid return to the pretreatment community composition is indicative of factors promoting community resilience, the nature of which deserves future investigation.

Deep sequencing of PCR amplicon libraries facilitates the detection of low-abundance populations in environmental DNA surveys of complex microbial communities. At the same time, deep sequencing can lead to overestimates of microbial diversity through the generation of low-frequency, error-prone reads. Even with sequencing error rates below 0.005 per nucleotide position, the common method of generating operational taxonomic units (OTUs) by multiple sequence alignment and complete-linkage clustering significantly increases the number of predicted OTUs and inflates richness estimates. We show that a 2% single-linkage preclustering methodology followed by an average-linkage clustering based on pairwise alignments more accurately predicts expected OTUs in both single and pooled template preparations of known taxonomic composition. This new clustering method can reduce the OTU richness in environmental samples by as much as 30-60% but does not reduce the fraction of OTUs in long-tailed rank abundance curves that defines the rare biosphere.

Massively parallel pyrosequencing systems have increased the efficiency of DNA sequencing, although the published per-base accuracy of a Roche GS20 is only 96%. In genome projects, highly redundant consensus assemblies can compensate for sequencing errors. In contrast, studies of microbial diversity that catalogue differences between PCR amplicons of ribosomal RNA genes (rDNA) or other conserved gene families cannot take advantage of consensus assemblies to detect and minimize incorrect base calls.We performed an empirical study of the per-base error rate for the Roche GS20 system using sequences of the V6 hypervariable region from cloned microbial ribosomal DNA (tag sequencing). We calculated a 99.5% accuracy rate in unassembled sequences, and identified several factors that can be used to remove a small percentage of low-quality reads, improving the accuracy to 99.75% or better.By using objective criteria to eliminate low quality data, the quality of individual GS20 sequence reads in molecular ecological applications can surpass the accuracy of traditional capillary methods.

Most studies examining the commensal human oral microbiome are focused on disease or are limited in methodology. In order to diagnose and treat diseases at an early and reversible stage an in-depth definition of health is indispensible. The aim of this study therefore was to define the healthy oral microbiome using recent advances in sequencing technology (454 pyrosequencing).We sampled and sequenced microbiomes from several intraoral niches (dental surfaces, cheek, hard palate, tongue and saliva) in three healthy individuals. Within an individual oral cavity, we found over 3600 unique sequences, over 500 different OTUs or "species-level" phylotypes (sequences that clustered at 3% genetic difference) and 88 - 104 higher taxa (genus or more inclusive taxon). The predominant taxa belonged to Firmicutes (genus Streptococcus, family Veillonellaceae, genus Granulicatella), Proteobacteria (genus Neisseria, Haemophilus), Actinobacteria (genus Corynebacterium, Rothia, Actinomyces), Bacteroidetes (genus Prevotella, Capnocytophaga, Porphyromonas) and Fusobacteria (genus Fusobacterium).Each individual sample harboured on average 266 "species-level" phylotypes (SD 67; range 123 - 326) with cheek samples being the least diverse and the dental samples from approximal surfaces showing the highest diversity. Principal component analysis discriminated the profiles of the samples originating from shedding surfaces (mucosa of tongue, cheek and palate) from the samples that were obtained from solid surfaces (teeth).There was a large overlap in the higher taxa, "species-level" phylotypes and unique sequences among the three microbiomes: 84% of the higher taxa, 75% of the OTUs and 65% of the unique sequences were present in at least two of the three microbiomes. The three individuals shared 1660 of 6315 unique sequences. These 1660 sequences (the "core microbiome") contributed 66% of the reads. The overlapping OTUs contributed to 94% of the reads, while nearly all reads (99.8%) belonged to the shared higher taxa.We obtained the first insight into the diversity and uniqueness of individual oral microbiomes at a resolution of next-generation sequencing. We showed that a major proportion of bacterial sequences of unrelated healthy individuals is identical, supporting the concept of a core microbiome at health.

Here we describe, the longest microbial time-series analyzed to date using high-resolution 16S rRNA tag pyrosequencing of samples taken monthly over 6 years at a temperate marine coastal site off Plymouth, UK. Data treatment effected the estimation of community richness over a 6-year period, whereby 8794 operational taxonomic units (OTUs) were identified using single-linkage preclustering and 21 130 OTUs were identified by denoising the data. The Alphaproteobacteria were the most abundant Class, and the most frequently recorded OTUs were members of the Rickettsiales (SAR 11) and Rhodobacteriales. This near-surface ocean bacterial community showed strong repeatable seasonal patterns, which were defined by winter peaks in diversity across all years. Environmental variables explained far more variation in seasonally predictable bacteria than did data on protists or metazoan biomass. Change in day length alone explains >65% of the variance in community diversity. The results suggested that seasonal changes in environmental variables are more important than trophic interactions. Interestingly, microbial association network analysis showed that correlations in abundance were stronger within bacterial taxa rather than between bacteria and eukaryotes, or between bacteria and environmental variables.

Massively parallel pyrosequencing of hypervariable regions from small subunit ribosomal RNA (SSU rRNA) genes can sample a microbial community two or three orders of magnitude more deeply per dollar and per hour than capillary sequencing of full-length SSU rRNA. As with full-length rRNA surveys, each sequence read is a tag surrogate for a single microbe. However, rather than assigning taxonomy by creating gene trees de novo that include all experimental sequences and certain reference taxa, we compare the hypervariable region tags to an extensive database of rRNA sequences and assign taxonomy based on the best match in a Global Alignment for Sequence Taxonomy (GAST) process. The resulting taxonomic census provides information on both composition and diversity of the microbial community. To determine the effectiveness of using only hypervariable region tags for assessing microbial community membership, we compared the taxonomy assigned to the V3 and V6 hypervariable regions with the taxonomy assigned to full-length SSU rRNA sequences isolated from both the human gut and a deep-sea hydrothermal vent. The hypervariable region tags and full-length rRNA sequences provided equivalent taxonomy and measures of relative abundance of microbial communities, even for tags up to 15% divergent from their nearest reference match. The greater sampling depth per dollar afforded by massively parallel pyrosequencing reveals many more members of the "rare biosphere" than does capillary sequencing of the full-length gene. In addition, tag sequencing eliminates cloning bias and the sequences are short enough to be completely sequenced in a single read, maximizing the number of organisms sampled in a run while minimizing chimera formation. This technique allows the cost-effective exploration of changes in microbial community structure, including the rare biosphere, over space and time and can be applied immediately to initiatives, such as the Human Microbiome Project.

The analytical power of environmental DNA sequences for modeling microbial ecosystems depends on accurate assessments of population structure, including diversity (richness) and relative abundance (evenness). We investigated both aspects of population structure for microbial communities at two neighboring hydrothermal vents by examining the sequences of more than 900,000 microbial small-subunit ribosomal RNA amplicons. The two vent communities have different population structures that reflect local geochemical regimes. Descriptions of archaeal diversity were nearly exhaustive, but despite collecting an unparalleled number of sequences, statistical analyses indicated additional bacterial diversity at every taxonomic level. We predict that hundreds of thousands of sequences will be necessary to capture the vast diversity of microbial communities, and that different patterns of evenness for both high- and low-abundance taxa may be important in defining microbial ecosystem dynamics.

The genome of the eukaryotic protist Giardia lamblia, an important human intestinal parasite, is compact in structure and content, contains few introns or mitochondrial relics, and has simplified machinery for DNA replication, transcription, RNA processing, and most metabolic pathways. Protein kinases comprise the single largest protein class and reflect Giardia's requirement for a complex signal transduction network for coordinating differentiation. Lateral gene transfer from bacterial and archaeal donors has shaped Giardia's genome, and previously unknown gene families, for example, cysteine-rich structural proteins, have been discovered. Unexpectedly, the genome shows little evidence of heterozygosity, supporting recent speculations that this organism is sexual. This genome sequence will not only be valuable for investigating the evolution of eukaryotes, but will also be applied to the search for new therapeutics for this parasite.

Background Massively parallel pyrosequencing of amplicons from the V6 hypervariable regions of small-subunit (SSU) ribosomal RNA (rRNA) genes is commonly used to assess diversity and richness in bacterial and archaeal populations. Recent advances in pyrosequencing technology provide read lengths of up to 240 nucleotides. Amplicon pyrosequencing can now be applied to longer variable regions of the SSU rRNA gene including the V9 region in eukaryotes. Methodology/Principal Findings We present a protocol for the amplicon pyrosequencing of V9 regions for eukaryotic environmental samples for biodiversity inventories and species richness estimation. The International Census of Marine Microbes (ICoMM) and the Microbial Inventory Research Across Diverse Aquatic Long Term Ecological Research Sites (MIRADA-LTERs) projects are already employing this protocol for tag sequencing of eukaryotic samples in a wide diversity of both marine and freshwater environments. Conclusions/Significance Massively parallel pyrosequencing of eukaryotic V9 hypervariable regions of SSU rRNA genes provides a means of estimating species richness from deeply-sampled populations and for discovering novel species from the environment.

We present the Biological Observation Matrix (BIOM, pronounced “biome”) format: a JSON-based file format for representing arbitrary observation by sample contingency tables with associated sample and observation metadata. As the number of categories of comparative omics data types (collectively, the “ome-ome”) grows rapidly, a general format to represent and archive this data will facilitate the interoperability of existing bioinformatics tools and future meta-analyses. The BIOM file format is supported by an independent open-source software project (the biom-format project), which initially contains Python objects that support the use and manipulation of BIOM data in Python programs, and is intended to be an open development effort where developers can submit implementations of these objects in other programming languages. The BIOM file format and the biom-format project are steps toward reducing the “bioinformatics bottleneck” that is currently being experienced in diverse areas of biological sciences, and will help us move toward the next phase of comparative omics where basic science is translated into clinical and environmental applications. The BIOM file format is currently recognized as an Earth Microbiome Project Standard, and as a Candidate Standard by the Genomic Standards Consortium.

A good definition of commensal microflora and an understanding of its relation to health are essential in preventing and combating disease. We hypothesized that the species richness of human oral microflora is underestimated. Saliva and supragingival plaque were sampled from 71 and 98 healthy adults, respectively. Amplicons from the V6 hypervariable region of the small-subunit ribosomal RNA gene were generated by PCR, pooled into saliva and plaque pools, and sequenced by means of the Genome Sequencer 20 system at 454 Life Sciences. Data were evaluated by taxonomic and rarefaction analyses. The 197,600 sequences generated yielded about 29,000 unique sequences, representing 22 taxonomic phyla. Grouping the sequences in operational taxonomic units (6%) yielded 3621 and 6888 species-level phylotypes in saliva and plaque, respectively. This work gives a radically new insight into the diversity of human oral microflora, which, with an estimated number of 19,000 phylotypes, is considerably higher than previously reported.

Marine microbial communities have been essential contributors to global biomass, nutrient cycling, and biodiversity since the early history of Earth, but so far their community distribution patterns remain unknown in most marine ecosystems.The synthesis of 9.6 million bacterial V6-rRNA amplicons for 509 samples that span the global ocean's surface to the deep-sea floor shows that pelagic and benthic communities greatly differ, at all taxonomic levels, and share <10% bacterial types defined at 3% sequence similarity level. Surface and deep water, coastal and open ocean, and anoxic and oxic ecosystems host distinct communities that reflect productivity, land influences and other environmental constraints such as oxygen availability. The high variability of bacterial community composition specific to vent and coastal ecosystems reflects the heterogeneity and dynamic nature of these habitats. Both pelagic and benthic bacterial community distributions correlate with surface water productivity, reflecting the coupling between both realms by particle export. Also, differences in physical mixing may play a fundamental role in the distribution patterns of marine bacteria, as benthic communities showed a higher dissimilarity with increasing distance than pelagic communities.This first synthesis of global bacterial distribution across different ecosystems of the World's oceans shows remarkable horizontal and vertical large-scale patterns in bacterial communities. This opens interesting perspectives for the definition of biogeographical biomes for bacteria of ocean waters and the seabed.

We present a performance-optimized algorithm, subsampled open-reference OTU picking, for assigning marker gene (e.g., 16S rRNA) sequences generated on next-generation sequencing platforms to operational taxonomic units (OTUs) for microbial community analysis. This algorithm provides benefits over de novo OTU picking (clustering can be performed largely in parallel, reducing runtime) and closed-reference OTU picking (all reads are clustered, not only those that match a reference database sequence with high similarity). Because more of our algorithm can be run in parallel relative to "classic" open-reference OTU picking, it makes open-reference OTU picking tractable on massive amplicon sequence data sets (though on smaller data sets, "classic" open-reference OTU clustering is often faster). We illustrate that here by applying it to the first 15,000 samples sequenced for the Earth Microbiome Project (1.3 billion V4 16S rRNA amplicons). To the best of our knowledge, this is the largest OTU picking run ever performed, and we estimate that our new algorithm runs in less than 1/5 the time than would be required of "classic" open reference OTU picking. We show that subsampled open-reference OTU picking yields results that are highly correlated with those generated by "classic" open-reference OTU picking through comparisons on three well-studied datasets. An implementation of this algorithm is provided in the popular QIIME software package, which uses uclust for read clustering. All analyses were performed using QIIME's uclust wrappers, though we provide details (aided by the open-source code in our GitHub repository) that will allow implementation of subsampled open-reference OTU picking independently of QIIME (e.g., in a compiled programming language, where runtimes should be further reduced). Our analyses should generalize to other implementations of these OTU picking algorithms. Finally, we present a comparison of parameter settings in QIIME's OTU picking workflows and make recommendations on settings for these free parameters to optimize runtime without reducing the quality of the results. These optimized parameters can vastly decrease the runtime of uclust-based OTU picking in QIIME.

ABSTRACT Shifts in microbial communities are implicated in the pathogenesis of a number of gastrointestinal diseases, but we have limited understanding of the mechanisms that lead to altered community structures. One difficulty with studying these mechanisms in human subjects is the inherent baseline variability of the microbiota in different individuals. In an effort to overcome this baseline variability, we employed a mouse model to control the host genotype, diet, and other possible influences on the microbiota. This allowed us to determine whether the indigenous microbiota in such mice had a stable baseline community structure and whether this community exhibited a consistent response following antibiotic administration. We employed a tag-sequencing strategy targeting the V6 hypervariable region of the bacterial small-subunit (16S) rRNA combined with massively parallel sequencing to determine the community structure of the gut microbiota. Inbred mice in a controlled environment harbored a reproducible baseline community that was significantly impacted by antibiotic administration. The ability of the gut microbial community to recover to baseline following the cessation of antibiotic administration differed according to the antibiotic regimen administered. Severe antibiotic pressure resulted in reproducible, long-lasting alterations in the gut microbial community, including a decrease in overall diversity. The finding of stereotypic responses of the indigenous microbiota to ecologic stress suggests that a better understanding of the factors that govern community structure could lead to strategies for the intentional manipulation of this ecosystem so as to preserve or restore a healthy microbiota.

We explore the microbiota of 18 body sites in over 200 individuals using sequences amplified V1-V3 and the V3-V5 small subunit ribosomal RNA (16S) hypervariable regions as part of the NIH Common Fund Human Microbiome Project. The body sites with the greatest number of core OTUs, defined as OTUs shared amongst 95% or more of the individuals, were the oral sites (saliva, tongue, cheek, gums, and throat) followed by the nose, stool, and skin, while the vaginal sites had the fewest number of OTUs shared across subjects. We found that commonalities between samples based on taxonomy could sometimes belie variability at the sub-genus OTU level. This was particularly apparent in the mouth where a given genus can be present in many different oral sites, but the sub-genus OTUs show very distinct site selection, and in the vaginal sites, which are consistently dominated by the Lactobacillus genus but have distinctly different sub-genus V1-V3 OTU populations across subjects. Different body sites show approximately a ten-fold difference in estimated microbial richness, with stool samples having the highest estimated richness, followed by the mouth, throat and gums, then by the skin, nasal and vaginal sites. Richness as measured by the V1-V3 primers was consistently higher than richness measured by V3-V5. We also show that when such a large cohort is analyzed at the genus level, most subjects fit the stool "enterotype" profile, but other subjects are intermediate, blurring the distinction between the enterotypes. When analyzed at the finer-scale, OTU level, there was little or no segregation into stool enterotypes, but in the vagina distinct biotypes were apparent. Finally, we note that even OTUs present in nearly every subject, or that dominate in some samples, showed orders of magnitude variation in relative abundance emphasizing the highly variable nature across individuals.

Very few marine microbial communities are well characterized even with the weight of research effort presently devoted to it. Only a small proportion of this effort has been aimed at investigating temporal community structure. Here we present the first report of the application of high-throughput pyrosequencing to investigate intra-annual bacterial community structure. Microbial diversity was determined for 12 time points at the surface of the L4 sampling site in the Western English Channel. This was performed over 11 months during 2007. A total of 182 560 sequences from the V6 hyper-variable region of the small-subunit ribosomal RNA gene (16S rRNA) were obtained; there were between 11 327 and 17 339 reads per sample. Approximately 7000 genera were identified, with one in every 25 reads being attributed to a new genus; yet this level of sampling far from exhausted the total diversity present at any one time point. The total data set contained 17 673 unique sequences. Only 93 (0.5%) were found at all time points, yet these few lineages comprised 50% of the total reads sequenced. The most abundant phylum was Proteobacteria (50% of all sequenced reads), while the SAR11 clade comprised 21% of the ubiquitous reads and approximately 12% of the total sequenced reads. In contrast, 78% of all operational taxonomic units were only found at one time point and 67% were only found once, evidence of a large and transient rare assemblage. This time series shows evidence of seasonally structured community diversity. There is also evidence for seasonal succession, primarily reflecting changes among dominant taxa. These changes in structure were significantly correlated to a combination of temperature, phosphate and silicate concentrations.

The nasopharynx is the ecological niche for many commensal bacteria and for potential respiratory or invasive pathogens like Streptococcus pneumoniae, Haemophilus influenzae, and Neisseria meningitidis. Disturbance of a balanced nasopharyngeal (NP) microbiome might be involved in the onset of symptomatic infections with these pathogens, which occurs primarily in fall and winter. It is unknown whether seasonal infection patterns are associated with concomitant changes in NP microbiota. As young children are generally prone to respiratory and invasive infections, we characterized the NP microbiota of 96 healthy children by barcoded pyrosequencing of the V5-V6 hypervariable region of the 16S-rRNA gene, and compared microbiota composition between children sampled in winter/fall with children sampled in spring. The approximately 1,000,000 sequences generated represented 13 taxonomic phyla and approximately 250 species-level phyla types (OTUs). The 5 most predominant phyla were Proteobacteria (64%), Firmicutes (21%), Bacteroidetes (11%), Actinobacteria (3%) and Fusobacteria (1,4%) with Moraxella, Haemophilus, Streptococcus, Flavobacteria, Dolosigranulum, Corynebacterium and Neisseria as predominant genera. The inter-individual variability was that high that on OTU level a core microbiome could not be defined. Microbiota profiles varied strongly with season, with in fall/winter a predominance of Proteobacteria (relative abundance (% of all sequences): 75% versus 51% in spring) and Fusobacteria (absolute abundance (% of children): 14% versus 2% in spring), and in spring a predominance of Bacteroidetes (relative abundance: 19% versus 3% in fall/winter, absolute abundance: 91% versus 54% in fall/winter), and Firmicutes. The latter increase is mainly due to (Brevi)bacillus and Lactobacillus species (absolute abundance: 96% versus 10% in fall/winter) which are like Bacteroidetes species generally related to healthy ecosystems. The observed seasonal effects could not be attributed to recent antibiotics or viral co-infection.The NP microbiota of young children is highly diverse and appears different between seasons. These differences seem independent of antibiotic use or viral co-infection.

Summary The release of untreated sewage introduces non‐indigenous microbial populations of uncertain composition into surface waters. We used massively parallel 454 pyrosequencing of hypervariable regions in rRNA genes to profile microbial communities from eight untreated sewage influent samples of two wastewater treatment plants (WWTPs) in metropolitan Milwaukee. The sewage profiles included a discernible human faecal signature made up of several taxonomic groups including multiple Bifidobacteriaceae , Coriobacteriaceae , Bacteroidaceae , Lachnospiraceae and Ruminococcaceae genera. The faecal signature made up a small fraction of the taxa present in sewage but the relative abundance of these sequence tags mirrored the population structures of human faecal samples. These genera were much more prevalent in the sewage influent than standard indicators species. High‐abundance sequences from taxonomic groups within the Beta ‐ and Gammaproteobacteria dominated the sewage samples but occurred at very low levels in faecal and surface water samples, suggesting that these organisms proliferate within the sewer system. Samples from Jones Island (JI – servicing residential plus a combined sewer system) and South Shore (SS – servicing a residential area) WWTPs had very consistent community profiles, with greater similarity between WWTPs on a given collection day than the same plant collected on different days. Rainfall increased influent flows at SS and JI WWTPs, and this corresponded to greater diversity in the community at both plants. Overall, the sewer system appears to be a defined environment with both infiltration of rainwater and stormwater inputs modulating community composition. Microbial sewage communities represent a combination of inputs from human faecal microbes and enrichment of specific microbes from the environment to form a unique population structure.

In a series of studies of the gut microbiome, “enterotypes” have been used to classify gut microbiome samples that cluster together in ordination analyses. Initially, three distinct enterotypes were described, although later studies reduced this to two clusters, one dominated by Bacteroides or Clostridiales species found more commonly in Western (American and Western European) subjects and the other dominated by Prevotella more often associated with non-Western subjects. The two taxa, Bacteroides and Prevotella, have been presumed to represent consistent underlying microbial communities, but no one has demonstrated the presence of additional microbial taxa across studies that can define these communities. We analyzed the combined microbiome data from five previous studies with samples across five continents. We clearly demonstrate that there are no consistent bacterial taxa associated with either Bacteroides- or Prevotella-dominated communities across the studies. By increasing the number and diversity of samples, we found gradients of both Bacteroides and Prevotella and a lack of the distinct clusters in the principal coordinate plots originally proposed in the “enterotypes” hypothesis. The apparent segregation of the samples seen in many ordination plots is due to the differences in the samples’ Prevotella and Bacteroides abundances and does not represent consistent microbial communities within the “enterotypes” and is not associated with other taxa across studies. The projections we see are consistent with a continuum of values created from a simple mixture of Bacteroides and Prevotella; these two biomarkers are significantly correlated to the projection axes. We suggest that previous findings citing Bacteroides- and Prevotella-dominated clusters are the result of an artifact caused by the greater relative abundance of these two taxa over other taxa in the human gut and the sparsity of Prevotella abundant samples. We believe that the term “enterotypes” is misleading because it implies both an underlying consistency of community taxa and a clear separation of sets of human gut samples, neither of which is supported by the broader data. We propose the use of “biomarker” as a more accurate description of these and other taxa that correlate with diet, lifestyle, and disease state.

ABSTRACT The fecal microbiome of cattle plays a critical role not only in animal health and productivity but also in food safety, pathogen shedding, and the performance of fecal pollution detection methods. Unfortunately, most published molecular surveys fail to provide adequate detail about variability in the community structures of fecal bacteria within and across cattle populations. Using massively parallel pyrosequencing of a hypervariable region of the rRNA coding region, we profiled the fecal microbial communities of cattle from six different feeding operations where cattle were subjected to consistent management practices for a minimum of 90 days. We obtained a total of 633,877 high-quality sequences from the fecal samples of 30 adult beef cattle (5 individuals per operation). Sequence-based clustering and taxonomic analyses indicate less variability within a population than between populations. Overall, bacterial community composition correlated significantly with fecal starch concentrations, largely reflected in changes in the Bacteroidetes , Proteobacteria , and Firmicutes populations. In addition, network analysis demonstrated that annotated sequences clustered by management practice and fecal starch concentration, suggesting that the structures of bovine fecal bacterial communities can be dramatically different in different animal feeding operations, even at the phylum and family taxonomic levels, and that the feeding operation is a more important determinant of the cattle microbiome than is the geographic location of the feedlot.

To date, metagenomic studies have relied on the utilization and analysis of reads obtained using 454 pyrosequencing to replace conventional Sanger sequencing. After extensively scanning the 16S ribosomal RNA (rRNA) gene, we identified the V5 hypervariable region as a short region providing reliable identification of bacterial sequences available in public databases such as the Human Oral Microbiome Database. We amplified samples from the oral cavity of three healthy individuals using primers covering an ~ 82-base segment of the V5 loop, and sequenced using the Illumina technology in a single orientation. We identified 135 genera or higher taxonomic ranks from the resulting 1,373,824 sequences. While the abundances of the most common phyla (Firmicutes, Proteobacteria, Actinobacteria, Fusobacteria and TM7) are largely comparable to previous studies, Bacteroidetes were less present. Potential sources for this difference include classification bias in this region of the 16S rRNA gene, human sample variation, sample preparation and primer bias. Using an Illumina sequencing approach, we achieved a much greater depth of coverage than previous oral microbiota studies, allowing us to identify several taxa not yet discovered in these types of samples, and to assess that at least 30,000 additional reads would be required to identify only one additional phylotype. The evolution of high-throughput sequencing technologies, and their subsequent improvements in read length enable the utilization of different platforms for studying communities of complex flora. Access to large amounts of data is already leading to a better representation of sample diversity at a reasonable cost.

The Human Microbiome Project provided a census of bacterial populations in healthy individuals, but an understanding of the biomedical significance of this census has been hindered by limited taxonomic resolution. A high-resolution method termed oligotyping overcomes this limitation by evaluating individual nucleotide positions using Shannon entropy to identify the most information-rich nucleotide positions, which then define oligotypes. We have applied this method to comprehensively analyze the oral microbiome. Using Human Microbiome Project 16S rRNA gene sequence data for the nine sites in the oral cavity, we identified 493 oligotypes from the V1-V3 data and 360 oligotypes from the V3-V5 data. We associated these oligotypes with species-level taxon names by comparison with the Human Oral Microbiome Database. We discovered closely related oligotypes, differing sometimes by as little as a single nucleotide, that showed dramatically different distributions among oral sites and among individuals. We also detected potentially pathogenic taxa in high abundance in individual samples. Numerous oligotypes were preferentially located in plaque, others in keratinized gingiva or buccal mucosa, and some oligotypes were characteristic of habitat groupings such as throat, tonsils, tongue dorsum, hard palate, and saliva. The differing habitat distributions of closely related oligotypes suggest a level of ecological and functional biodiversity not previously recognized. We conclude that the Shannon entropy approach of oligotyping has the capacity to analyze entire microbiomes, discriminate between closely related but distinct taxa and, in combination with habitat analysis, provide deep insight into the microbial communities in health and disease.

Abstract Background An understanding of the relation of commensal microbiota to health is essential in preventing disease. Here we studied the oral microbial composition of children (N = 74, aged 3 - 18 years) in natural transition from their deciduous to a permanent dentition and related the microbial profiles to their oral health status. The microbial composition of saliva was assessed by barcoded pyrosequencing of the V5-V6 hypervariable regions of the 16 S rRNA, as well as by using phylogenetic microarrays. Results Pyrosequencing reads (126174 reads, 1045 unique sequences) represented 8 phyla and 113 higher taxa in saliva samples. Four phyla - Firmicutes, Bacteriodetes, Proteobacteria and Actinobacteria - predominated in all groups. The deciduous dentition harboured a higher proportion of Proteobacteria (Gammaproteobacteria, Moraxellaceae) than Bacteroidetes, while in all other groups Bacteroidetes were at least as abundant as Proteobacteria. Bacteroidetes (mainly genus Prevotella ), Veillonellaceae family, Spirochaetes and candidate division TM7 increased with increasing age, reflecting maturation of the microbiome driven by biological changes with age. Microarray analysis enabled further analysis of the individual salivary microbiota. Of 350 microarray probes, 156 gave a positive signal with, on average, 77 (range 48-93) probes per individual sample. A caries-free oral status significantly associated with the higher signal of the probes targeting Porphyromonas catoniae and Neisseria flavescens . Conclusions The potential role of P. catoniae and N. flavescens as oral health markers should be assessed in large-scale clinical studies. The combination of both, open-ended and targeted molecular approaches provides us with information that will increase our understanding of the interplay between the human host and its microbiome.

How microbial communities change over time in response to the environment is poorly understood. Previously a six-year time series of 16S rRNA V6 data from the Western English Channel demonstrated robust seasonal structure within the bacterial community, with diversity negatively correlated with day-length. Here we determine whether metagenomes and metatranscriptomes follow similar patterns. We generated 16S rRNA datasets, metagenomes (1.2 GB) and metatranscriptomes (157 MB) for eight additional time points sampled in 2008, representing three seasons (Winter, Spring, Summer) and including day and night samples. This is the first microbial 'multi-omic' study to combine 16S rRNA amplicon sequencing with metagenomic and metatranscriptomic profiling. Five main conclusions can be drawn from analysis of these data: 1) Archaea follow the same seasonal patterns as Bacteria, but show lower relative diversity; 2) Higher 16S rRNA diversity also reflects a higher diversity of transcripts; 3) Diversity is highest in winter and at night; 4) Community-level changes in 16S-based diversity and metagenomic profiles are better explained by seasonal patterns (with samples closest in time being most similar), while metatranscriptomic profiles are better explained by diel patterns and shifts in particular categories (i.e., functional groups) of genes; 5) Changes in key genes occur among seasons and between day and night (i.e., photosynthesis); but these samples contain large numbers of orphan genes without known homologues and it is these unknown gene sets that appear to contribute most towards defining the differences observed between times. Despite the huge diversity of these microbial communities, there are clear signs of predictable patterns and detectable stability over time. Renewed and intensified efforts are required to reveal fundamental deterministic patterns in the most complex microbial communities. Further, the presence of a substantial proportion of orphan sequences underscores the need to determine the gene products of sequences with currently unknown function.

Abstract Coastal sands filter and accumulate organic and inorganic materials from the terrestrial and marine environment, and thus provide a high diversity of microbial niches. Sands of temperate climate zones represent a temporally and spatially highly dynamic marine environment characterized by strong physical mixing and seasonal variation. Yet little is known about the temporal fluctuations of resident and rare members of bacterial communities in this environment. By combining community fingerprinting via pyrosequencing of ribosomal genes with the characterization of multiple environmental parameters, we disentangled the effects of seasonality, environmental heterogeneity, sediment depth and biogeochemical gradients on the fluctuations of bacterial communities of marine sands. Surprisingly, only 3–5% of all bacterial types of a given depth zone were present at all times, but 50–80% of them belonged to the most abundant types in the data set. About 60–70% of the bacterial types consisted of tag sequences occurring only once over a period of 1 year. Most members of the rare biosphere did not become abundant at any time or at any sediment depth, but varied significantly with environmental parameters associated with nutritional stress. Despite the large proportion and turnover of rare organisms, the overall community patterns were driven by deterministic relationships associated with seasonal fluctuations in key biogeochemical parameters related to primary productivity. The maintenance of major biogeochemical functions throughout the observation period suggests that the small proportion of resident bacterial types in sands perform the key biogeochemical processes, with minimal effects from the rare fraction of the communities.

The advent of next-generation DNA sequencing platforms has revolutionized molecular microbial ecology by making the detailed analysis of complex communities over time and space a tractable research pursuit for small research groups. However, the ability to generate 10⁵-10⁸ reads with relative ease brings with it many downstream complications. Beyond the computational resources and skills needed to process and analyze data, it is difficult to compare datasets in an intuitive and interactive manner that leads to hypothesis generation and testing.We developed the free web service VAMPS (Visualization and Analysis of Microbial Population Structures, http://vamps.mbl.edu) to address these challenges and to facilitate research by individuals or collaborating groups working on projects with large-scale sequencing data. Users can upload marker gene sequences and associated metadata; reads are quality filtered and assigned to both taxonomic structures and to taxonomy-independent clusters. A simple point-and-click interface allows users to select for analysis any combination of their own or their collaborators' private data and data from public projects, filter these by their choice of taxonomic and/or abundance criteria, and then explore these data using a wide range of analytic methods and visualizations. Each result is extensively hyperlinked to other analysis and visualization options, promoting data exploration and leading to a greater understanding of data relationships.VAMPS allows researchers using marker gene sequence data to analyze the diversity of microbial communities and the relationships between communities, to explore these analyses in an intuitive visual context, and to download data, results, and images for publication. VAMPS obviates the need for individual research groups to make the considerable investment in computational infrastructure and bioinformatic support otherwise necessary to process, analyze, and interpret massive amounts of next-generation sequence data. Any web-capable device can be used to upload, process, explore, and extract data and results from VAMPS. VAMPS encourages researchers to share sequence and metadata, and fosters collaboration between researchers of disparate biomes who recognize common patterns in shared data.

Verrucomicrobia is a bacterial phylum that is commonly detected in soil, but little is known about the distribution and diversity of this phylum in the marine environment. To address this, we analyzed the marine microbial community composition in 506 samples from the International Census of Marine Microbes as well as 11 coastal samples taken from the California Current. These samples from both the water column and sediments covered a wide range of environmental conditions. Verrucomicrobia were present in 98% of the analyzed samples, and thus appeared nearly ubiquitous in the ocean. Based on the occurrence of amplified 16S ribosomal RNA sequences, Verrucomicrobia constituted on average 2% of the water column and 1.4% of the sediment bacterial communities. The diversity of Verrucomicrobia displayed a biogeography at multiple taxonomic levels and thus, specific lineages appeared to have clear habitat preference. We found that subdivision 1 and 4 generally dominated marine bacterial communities, whereas subdivision 2 was more frequent in low salinity waters. Within the subdivisions, Verrucomicrobia community composition were significantly different in the water column compared with sediment as well as within the water column along gradients of salinity, temperature, nitrate, depth and overall water column depth. Although we still know little about the ecophysiology of Verrucomicrobia lineages, the ubiquity of this phylum suggests that it may be important for the biogeochemical cycle of carbon in the ocean.

We evaluated the population structure and temporal dynamics of the dominant community members within sewage influent from two wastewater treatment plants (WWTPs) in Milwaukee, WI. We generated > 1.1 M bacterial pyrotag sequences from the V6 hypervariable region of 16S rRNA genes from 38 influent samples and two samples taken upstream in the sanitary sewer system. Only a small fraction of pyrotags from influent samples (∼ 15%) matched sequences from human faecal samples. The faecal components of the sewage samples included enriched pyrotag populations from Lactococcus and Enterobacteriaceae relative to their fractional representation in human faecal samples. In contrast to the large number of distinct pyrotags that represent faecal bacteria such as Lachnospiraceae and Bacteroides, only one or two unique V6 sequences represented Acinetobacter, Aeromonas and Trichococcus, which collectively account for nearly 35% of the total sewage community. Two dominant Acinetobacter V6 pyrotags (designated Acineto tag 1 and Acineto tag 2) fluctuated inversely with a seasonal pattern over a 3-year period, suggesting two distinct Acinetobacter populations respond differently to ecological forcings in the system. A single nucleotide change in the V6 pyrotags accounted for the difference in these populations and corresponded to two phylogenetically distinct clades based on full-length sequences. Analysis of wavelet functions, derived from a mathematical model of temporal fluctuations, demonstrated that other abundant sewer associated populations including Trichococcus and Aeromonas had temporal patterns similar to either Acineto tag 1 or Acineto tag 2. Populations with related temporal fluctuations were found to significantly correlate with the same WWTP variables (5-day BOD, flow, ammonia, total phosphorous and suspended solids). These findings illustrate that small differences in V6 sequences can represent phylogenetically and ecologically distinct taxa. This work provides insight into microbial community composition and dynamics within the defined environment of urban sewer infrastructure.

Summary PCR‐based surveys of microbial communities commonly use regions of the small‐subunit ribosomal RNA (SSU rRNA) gene to determine taxonomic membership and estimate total diversity. Here we show that the length of the target amplicon has a significant effect on assessments of microbial richness and community membership. Using operational taxonomic unit (OTU)‐ and taxonomy‐based tools, we compared the V6 hypervariable region of the bacterial SSU rRNA gene of three amplicon libraries of c. 100, 400 and 1000 base pairs (bp) from each of two hydrothermal vent fluid samples. We found that the smallest amplicon libraries contained more unique sequences, higher diversity estimates and a different community structure than the other two libraries from each sample. We hypothesize that a combination of polymerase dissociation, cloning bias and mispriming due to secondary structure accounts for the differences. While this relationship is not linear, it is clear that the smallest amplicon libraries contained more different types of sequences, and accordingly, more diverse members of the community. Because divergent and lower abundant taxa can be more readily detected with smaller amplicons, they may provide better assessments of total community diversity and taxonomic membership than longer amplicons in molecular studies of microbial communities.

Recent advances in massively parallel sequencing technology have created new opportunities to probe the hidden world of microbes. Taxonomy-independent clustering of the 16S rRNA gene is usually the first step in analyzing microbial communities. Dozens of algorithms have been developed in the last decade, but a comprehensive benchmark study is lacking. Here, we survey algorithms currently used by microbiologists, and compare seven representative methods in a large-scale benchmark study that addresses several issues of concern. A new experimental protocol was developed that allows different algorithms to be compared using the same platform, and several criteria were introduced to facilitate a quantitative evaluation of the clustering performance of each algorithm. We found that existing methods vary widely in their outputs, and that inappropriate use of distance levels for taxonomic assignments likely resulted in substantial overestimates of biodiversity in many studies. The benchmark study identified our recently developed ESPRIT-Tree, a fast implementation of the average linkage-based hierarchical clustering algorithm, as one of the best algorithms available in terms of computational efficiency and clustering accuracy.

Background: Massively parallel pyrosequencing of amplicons from the V6 hypervariable regions of small-subunit (SSU) ribosomal RNA (rRNA) genes is commonly used to assess diversity and richness in bacterial and archaeal populations.Recent advances in pyrosequencing technology provide read lengths of up to 240 nucleotides.Amplicon pyrosequencing can now be applied to longer variable regions of the SSU rRNA gene including the V9 region in eukaryotes.Methodology/Principal Findings: We present a protocol for the amplicon pyrosequencing of V9 regions for eukaryotic environmental samples for biodiversity inventories and species richness estimation.The International Census of Marine Microbes (ICoMM) and the Microbial Inventory Research Across Diverse Aquatic Long Term Ecological Research Sites (MIRADA-LTERs) projects are already employing this protocol for tag sequencing of eukaryotic samples in a wide diversity of both marine and freshwater environments.Conclusions/Significance: Massively parallel pyrosequencing of eukaryotic V9 hypervariable regions of SSU rRNA genes provides a means of estimating species richness from deeply-sampled populations and for discovering novel species from the environment.

ABSTRACT Microbial sewage communities consist of a combination of human fecal microorganisms and nonfecal microorganisms, which may be residents of urban sewer infrastructure or flowthrough originating from gray water or rainwater inputs. Together, these different microorganism sources form an identifiable community structure that may serve as a signature for sewage discharges and as candidates for alternative indicators specific for human fecal pollution. However, the structure and variability of this community across geographic space remains uncharacterized. We used massively parallel 454 pyrosequencing of the V6 region in 16S rRNA genes to profile microbial communities from 13 untreated sewage influent samples collected from a wide range of geographic locations in the United States. We obtained a total of 380,175 high-quality sequences for sequence-based clustering, taxonomic analyses, and profile comparisons. The sewage profile included a discernible core human fecal signature made up of several abundant taxonomic groups within Firmicutes , Bacteroidetes , Actinobacteria , and Proteobacteria . DNA sequences were also classified into fecal, sewage infrastructure (i.e., nonfecal), and transient groups based on data comparisons with fecal samples. Across all sewage samples, an estimated 12.1% of sequences were fecal in origin, while 81.4% were consistently associated with the sewage infrastructure. The composition of feces-derived operational taxonomic units remained congruent across all sewage samples regardless of geographic locale; however, the sewage infrastructure community composition varied among cities, with city latitude best explaining this variation. Together, these results suggest that untreated sewage microbial communities harbor a core group of fecal bacteria across geographically dispersed wastewater sewage lines and that ambient water quality indicators targeting these select core microorganisms may perform well across the United States.

The operational taxonomic unit (OTU) is widely used in microbial ecology. Reproducibility in microbial ecology research depends on the reliability of OTU-based 16S ribosomal subunit RNA (rRNA) analyses.Here, we report that many hierarchical and greedy clustering methods produce unstable OTUs, with membership that depends on the number of sequences clustered. If OTUs are regenerated with additional sequences or samples, sequences originally assigned to a given OTU can be split into different OTUs. Alternatively, sequences assigned to different OTUs can be merged into a single OTU. This OTU instability affects alpha-diversity analyses such as rarefaction curves, beta-diversity analyses such as distance-based ordination (for example, Principal Coordinate Analysis (PCoA)), and the identification of differentially represented OTUs. Our results show that the proportion of unstable OTUs varies for different clustering methods. We found that the closed-reference method is the only one that produces completely stable OTUs, with the caveat that sequences that do not match a pre-existing reference sequence collection are discarded.As a compromise to the factors listed above, we propose using an open-reference method to enhance OTU stability. This type of method clusters sequences against a database and includes unmatched sequences by clustering them via a relatively stable de novo clustering method. OTU stability is an important consideration when analyzing microbial diversity and is a feature that should be taken into account during the development of novel OTU clustering methods.

ABSTRACT Dispersal limitation in phyllosphere communities was measured on the leaf surfaces of salt-excreting Tamarix trees, which offer unique, discrete habitats for microbial assemblages. We employed 16S rRNA gene pyrosequencing to measure bacterial community dissimilarity on leaves of spatially dispersed Tamarix specimens in sites with uniform climatic conditions across the Sonoran Desert in the Southwestern United States. Our analyses revealed diverse bacterial communities with four dominant phyla that exhibited differential effects of environmental and geographic variables. Geographical distance was the most important parameter that affected community composition, particularly that of betaproteobacteria, which displayed a statistically significant, distance-decay relationship.

Özok AR, Persoon IF, Huse SM, Keijser BJF, Wesselink PR, Crielaard W, Zaura E. Ecology of the microbiome of the infected root canal system: a comparison between apical and coronal root segments. International Endodontic Journal , 45 , 530–541, 2012. Abstract Aim To evaluate the microbial ecology of the coronal and apical segments of infected root canal systems using a complete sampling technique and next‐generation sequencing. Methodology The roots of 23 extracted teeth with apical periodontitis were sectioned in half, horizontally, and cryo‐pulverized. Bacterial communities were profiled using tagged 454 pyrosequencing of the 16S rDNA hypervariable V5‐V6 region. Results The sequences were classified into 606 taxa (species or higher taxon), representing 24 bacterial phyla or candidate divisions and one archaeal phylum. Proteobacteria were more abundant in the apical samples ( P &lt; 0.05), whilst Actinobacteria were in significantly higher proportions in the coronal samples. The apical samples harboured statistically significantly more taxa than the coronal samples ( P = 0.01) and showed a higher microbial diversity. Several taxa belonging to fastidious obligate anaerobes were significantly more abundant in the apical segments of the roots compared with their coronal counterparts. Conclusions Endodontic infections are more complex than reported previously. The apical part of the root canal system drives the selection of a more diverse and more anaerobic community than the coronal part. The presence of a distinct ecological niche in the apical region explains the difficulty of eradication of the infection and emphasizes the need for new treatment approaches to be developed.

Faecal pollution contains a rich and diverse community of bacteria derived from animals and humans, many of which might serve as alternatives to the traditional enterococci and Escherichia coli faecal indicators. We used massively parallel sequencing (MPS) of the 16S rRNA gene to characterize microbial communities from wastewater treatment plant (WWTP) influent sewage from 12 cities geographically distributed across the USA. We examined members of the Clostridiales, which included the families Clostridiaceae, Lachnospiraceae and Ruminococcaceae for their potential as sewage indicators. Lachnospiraceae was one of the most abundant groups of faecal bacteria in sewage, and several Lachnospiraceae high-abundance sewage pyrotags occurred in at least 46 of 48 human faecal samples. Clone libraries targeting Clostridium coccoides (C. coccoides) in sewage samples demonstrated that Lachnospiraceae-annotated V6 pyrotags encompassed the previously reported C. coccoides group. We used oligotyping to profile the genus Blautia within Lachnospiraceae and found oligotypes comprised of 24 entropy components that showed patterns of host specificity. These findings suggest that indicators based on Blautia might have the capacity to discriminate between different faecal pollution sources. Development of source-specific alternative indicators would enhance water quality assessments, which leads to improved ecosystem health and reduced human health risk due to waterborne disease.

Mucosal biopsy is the most common sampling technique used to assess microbial communities associated with the intestinal mucosa. Biopsies disrupt the epithelium and can be associated with complications such as bleeding. Biopsies sample a limited area of the mucosa, which can lead to potential sampling bias. In contrast to the mucosal biopsy, the mucosal brush technique is less invasive and provides greater mucosal coverage, and if it can provide equivalent microbial community data, it would be preferable to mucosal biopsies.We compared microbial samples collected from the intestinal mucosa using either a cytology brush or mucosal biopsy forceps. We collected paired samples from patients with ulcerative colitis (UC) who had previously undergone colectomy and ileal pouch anal anastomosis (IPAA), and profiled the microbial communities of the samples by sequencing V4-V6 or V4-V5 16S rRNA-encoding gene amplicons. Comparisons of 177 taxa in 16 brush-biopsy sample pairs had a mean R2 of 0.94. We found no taxa that varied significantly between the brush and biopsy samples after adjusting for multiple comparisons (false discovery rate ≤0.05). We also tested the reproducibility of DNA amplification and sequencing in 25 replicate pairs and found negligible variation (mean R2 = 0.99). A qPCR analysis of the two methods showed that the relative yields of bacterial DNA to human DNA were several-fold higher in the brush samples than in the biopsies.Mucosal brushing is preferred to mucosal biopsy for sampling the epithelial-associated microbiota. Although both techniques provide similar assessments of the microbial community composition, the brush sampling method has relatively more bacterial to host DNA, covers a larger surface area, and is less traumatic to the epithelium than the mucosal biopsy.

Massively parallel pyrosequencing of hypervariable regions from small subunit ribosomal RNA (SSU rRNA) genes can sample a microbial community two or three orders of magnitude more deeply per dollar and per hour than capillary sequencing of full-length SSU rRNA.As with full-length rRNA surveys, each sequence read is a tag surrogate for a single microbe.However, rather than assigning taxonomy by creating gene trees de novo that include all experimental sequences and certain reference taxa, we compare the hypervariable region tags to an extensive database of rRNA sequences and assign taxonomy based on the best match in a Global Alignment for Sequence Taxonomy (GAST) process.The resulting taxonomic census provides information on both composition and diversity of the microbial community.To determine the effectiveness of using only hypervariable region tags for assessing microbial community membership, we compared the taxonomy assigned to the V3 and V6 hypervariable regions with the taxonomy assigned to full-length SSU rRNA sequences isolated from both the human gut and a deep-sea hydrothermal vent.The hypervariable region tags and full-length rRNA sequences provided equivalent taxonomy and measures of relative abundance of microbial communities, even for tags up to 15% divergent from their nearest reference match.The greater sampling depth per dollar afforded by massively parallel pyrosequencing reveals many more members of the ''rare biosphere'' than does capillary sequencing of the full-length gene.In addition, tag sequencing eliminates cloning bias and the sequences are short enough to be completely sequenced in a single read, maximizing the number of organisms sampled in a run while minimizing chimera formation.This technique allows the cost-effective exploration of changes in microbial community structure, including the rare biosphere, over space and time and can be applied immediately to initiatives, such as the Human Microbiome Project.

Abstract Aim To test whether within‐species and among‐species patterns of abundance and latitudinal range in marine bacteria resemble those found for macro‐organisms, and whether these patterns differ along latitudinal clines. Location Global pelagic marine environments. Methods Taxon‐specific sequence abundance and location were retrieved from the open‐access V6‐ rRNA pyrotag sequence data base VAMPS ( http://vamps.mbl.edu/ ), which holds a massive collection of marine bacterial community data sets from the I nternational C ensus of M arine M icrobes sampling effort of global ocean water masses. Data were randomly subsampled to correct for spatial bias and for differences in sampling effort. Results We show that bacterial latitudinal ranges are narrower than expected by chance. When present in both N orthern and S outhern hemispheres, taxa occupy restricted ranges at similar latitudes on both sides of the equator. A significant and positive relationship exists between sequence abundance and latitudinal range, although this pattern contains a large amount of variance. Abundant taxa in the tropics and in the N orthern H emisphere generally have smaller ranges than those in the S outhern H emisphere. We show that the mean latitudinal range of bacterial taxa increases with latitude, supporting the existence of a R apoport effect in marine bacterioplankton. Finally, we show that bacterioplankton communities contain a higher proportion of abundant taxa as they approach the poles. Main conclusions Macroecological patterns such as the abundance–range relationship, in general, extend to marine bacteria. However, differences in the shape of these relationships between bacteria and macro‐organisms call into question whether the processes and their relative importance in shaping global marine bacteria and macro‐organism distributions are the same.

Low-temperature hydrothermal vent fluids represent access points to diverse microbial communities living in oceanic crust. This study examined the distribution, relative abundance, and diversity of Epsilonproteobacteria in 14 low-temperature vent fluids from five volcanically active seamounts of the Mariana Arc using a 454 tag sequencing approach. Most vent fluids were enriched in cell concentrations compared with background seawater, and quantitative PCR results indicated that all fluids were dominated by bacteria. Operational taxonomic unit-based statistical tools applied to 454 data show that all vents from the northern end of the Mariana Arc grouped together, to the exclusion of southern arc seamounts, which were as distinct from one another as they were from northern seamounts. Statistical analysis also showed a significant relationship between seamount and individual vent groupings, suggesting that community membership may be linked to geographical isolation and not geochemical parameters. However, while there may be large-scale geographic differences, distance is not the distinguishing factor in the microbial community composition. At the local scale, most vents host a distinct population of Epsilonproteobacteria, regardless of seamount location. This suggests that there may be barriers to exchange and dispersal for these vent endemic microorganisms at hydrothermal seamounts of the Mariana Arc.

Chapter 24 Microbial Diversity in the Deep Sea and the Underexplored “Rare Biosphere” David B. Mark Welch, David B. Mark Welch Marine Biological Laboratory at Woods Hole, Woods Hole, Massachusetts, USASearch for more papers by this authorSusan M. Huse, Susan M. Huse Marine Biological Laboratory at Woods Hole, Woods Hole, Massachusetts, USASearch for more papers by this author David B. Mark Welch, David B. Mark Welch Marine Biological Laboratory at Woods Hole, Woods Hole, Massachusetts, USASearch for more papers by this authorSusan M. Huse, Susan M. Huse Marine Biological Laboratory at Woods Hole, Woods Hole, Massachusetts, USASearch for more papers by this author Book Editor(s):Frans J. de Bruijn, Frans J. de Bruijn Laboratory of Plant Micro-organism Interaction, CNRS-INRA, Castanet Tolosan, FranceSearch for more papers by this author First published: 16 September 2011 https://doi.org/10.1002/9781118010549.ch24Citations: 17 AboutPDFPDF ToolsRequest permissionExport citationAdd to favoritesTrack citation ShareShareShare a linkShare onEmailFacebookTwitterLinkedInRedditWechat Summary This chapter contains sections titled: Introduction Methods Results and Discussion Internet Resources References INTERNET RESOURCES Primer sequences, the refSSU and refV6 databases, perl scripts used for SLP/PW-AN clustering, processed V6 tags, and other resources are available at http://vamps.mbl.edu. The RDP classifier is available at http://rdp.cme.msu.edu, and the Silva database of ribosomal sequences is available from www.arbsilva.de. Details of the ICoMM project can be found at http://icomm.mbl.edu. Google Scholar REFERENCES Google Scholar Acinas SG, Klepac-Ceraj V, Hunt DE, Pharino C, Ceraj I, Distel, DL. Polz, MF. 2004. Fine-scale phylogenetic architecture of a complex bacterial community. Nature 430: 551–554. 10.1038/nature02649 CASPubMedWeb of Science®Google Scholar Alain K, Querellou J, Lesongeur F, Pignet P, Crassous P, Raguenes G et al. 2002. Caminibacter hydrogeniphilus gen. nov., sp. nov., a novel thermophilic, hydrogen-oxidizing bacterium isolated from an East Pacific Rise hydrothermal vent. Int. J. Syst. Evol. Microbiol. 52: 1317–1323. 10.1099/ijs.0.02195-0 CASPubMedWeb of Science®Google Scholar Atlas RM, Bartha R. 1993. Microbial Ecology. Fundamentals and Applications. Redwood City, CA: Benjamin/Cummings. Google Scholar Béjà O, Suzuki MT, Koonin EV, Aravind L, Hadd A, Nguyen LP, et al. 2000. Construction and analysis of bacterial artificial chromo-some libraries from a marine microbial assemblage. Environ. Micro-biol. 2: 516–529. 10.1046/j.1462-2920.2000.00133.x CASPubMedWeb of Science®Google Scholar Campbell BJ, Jeanthon C, Kostka JE, Luther GW 3rd, Cary SC. 2001. Growth and phylogenetic properties of novel bacteria belonging to the epsilon subdivision of the Proteobacteria enriched from Alvinella pompejana and deep-sea hydrothermal vents. Appl. Environ. Microbiol. 67: 4566–4572. 10.1128/AEM.67.10.4566-4572.2001 CASPubMedWeb of Science®Google Scholar Chao A. 1984. Non-parametric estimation of the number of classes in a population. Scand. J. Stat. 11: 265–270. 10.1111/j.1462-2920.2006.01049.x PubMedWeb of Science®Google Scholar Chao A. 1992. Estimating the number of classes via sample coverage. J. Am. Stat. Assoc. 87: 210–217. 10.1080/01621459.1992.10475194 Web of Science®Google Scholar Chao A, Ma MC, Yang, K. 1993. Stopping rules and estimation for recapture debugging with unequal failure rates. Biometrika 80: 193–201. 10.1093/biomet/80.1.193 Web of Science®Google Scholar Cole JR, Chai B, Farris RJ, Wang Q, Kulam SA, McGarrell, DM, et al. 2005. The Ribosomal Database Project (RDP-II): sequences and tools for high-throughput rRNA analysis. Nucleic Acids Res. 33: D294–D296. 10.1093/nar/gki038 CASPubMedWeb of Science®Google Scholar Coleman ML, Sullivan MB, Martiny AC, Steglich C, Barry K, DeLong EF, Chisholm SW. 2006. Genomic islands and the ecology and evolution of Prochlorococcus. Science 311: 1768–1770. 10.1126/science.1122050 CASPubMedWeb of Science®Google Scholar Curtis TP, Sloan WT, Scannell, JW. 2002. Estimating prokaryotic diversity and its limits. Proc. Natl. Acad. Sci. USA 99: 10494–10499. 10.1073/pnas.142680199 CASPubMedWeb of Science®Google Scholar Curtis TP, Head IM, Lunn M, Woodcock S, Schloss PD, Sloan WT. 2006. What is the extent of microbial diversity? Philos. Trans. R. Soc. Lond. B 361: 2023–2037. 10.1098/rstb.2006.1921 PubMedWeb of Science®Google Scholar DeLong EF. 1992. Archaea in coastal marine environments. Proc. Natl. Acad. Sci. USA 89: 5685–5689. 10.1073/pnas.89.12.5685 CASPubMedWeb of Science®Google Scholar DeLong EF, Preston CM, Mincer T, Rich V, Hallam SJ, Frigaard NU, et al. 2006. Community genomics among stratified microbial assemblages in the ocean’s interior. Science 311: 496–503. 10.1126/science.1120250 CASPubMedWeb of Science®Google Scholar Frigaard NU, Martinez A, Mincer, TJ, DeLong, EF. 2006. Proteorhodopsin lateral gene transfer between marine planktonic Bacteria and Archaea. Nature 439: 847–850. 10.1038/nature04435 CASPubMedWeb of Science®Google Scholar Fuhrman JA, McCallum K, Davis, AA. 1993. Phylogenetic diversity of subsurface marine microbial communities from the Atlantic and Pacific Oceans. Appl. Environ. Microbiol. 59: 1294–1302. 10.1128/AEM.59.5.1294-1302.1993 CASPubMedWeb of Science®Google Scholar Giovannoni SJ, Britschgi TB, Moyer CL, Field KG. 1990. Genetic diversity in Sargasso Sea bacterioplankton. Nature 345: 60–63. 10.1038/345060a0 CASPubMedWeb of Science®Google Scholar Hallam SJ, Mincer TJ, Schleper C, Preston CM, Roberts K, Richardson PM, Delong, EF. 2006. Pathways of carbon assimilation and ammonia oxidation suggested by environmental genomic analyses of marine crenarchaeota. PLoS Biol. 4: e95. 10.1371/journal.pbio.0040095 CASPubMedWeb of Science®Google Scholar Hebert PD, Cywinska A, Ball, SL, deWaard, JR. 2003. Biological identifications through DNA barcodes. Proc. Biol. Sci. 270: 313–321. 10.1098/rspb.2002.2218 CASPubMedWeb of Science®Google Scholar Hewson I, Vargo GA, Fuhrman JA. 2003. Bacterial diversity in shallow oligotrophic marine benthos and overlying waters: Effects of virus infection, containment, and nutrient enrichment. Microb. Ecol. 46: 322–336. 10.1007/s00248-002-1067-3 CASPubMedWeb of Science®Google Scholar Hobbie JE, Daley RJ, Jasper S. 1977. Use of nuclepore filters for counting bacteria by fluorescence microscopy. Appl. Environ. Micro-biol. 33: 1225–1228. PubMedWeb of Science®Google Scholar Huber JA, Butterfield DA, Baross JA. 2003. Bacterial diversity in a subseafloor habitat following a deep-sea volcanic eruption. FEMS Microbiol. Ecol. 43: 393–409. 10.1111/j.1574-6941.2003.tb01080.x CASPubMedWeb of Science®Google Scholar Huber JA, Johnson HP, Butterfield DA, Baross, JA. 2006. Microbial life in ridge flank crustal fluids. Environ. Microbiol. 8: 88–99. 10.1111/j.1462-2920.2005.00872.x CASPubMedWeb of Science®Google Scholar Huber JA, Mark Welch DB, Morrison HG, Huse SM, Neal PR, Butterfield DA, Sogin, ML. 2007. Microbial population structures in the deep marine biosphere. Science 318: 97–100. 10.1126/science.1146689 CASPubMedWeb of Science®Google Scholar Hugler M, Wirsen CO, Fuchs G, Taylor CD, Sievert SM. 2005. Evidence for autotrophic CO2 fixation via the reductive tricarboxylic acid cycle by members of the epsilon subdivision of proteobacteria. J. Bacteriol. 187: 3020–3027. 10.1128/JB.187.9.3020-3027.2005 CASPubMedWeb of Science®Google Scholar Huse S, Huber J, Morrison H, Sogin M, Mark Welch D. 2007. Accuracy and quality of massively parallel DNA pyrosequencing. Genome Biol. 8: R143. 10.1186/gb-2007-8-7-r143 PubMedWeb of Science®Google Scholar Huse SM, Dethlefsen L, Huber JA, Mark Welch D, Relman DA, Sogin ML. 2008. Exploring microbial diversity and taxonomy using SSU rRNA hypervariable tag sequencing. PLoS Genetics 4: e1000255. 10.1371/journal.pgen.1000255 CASPubMedWeb of Science®Google Scholar Huse SM, Mark Welch DB, Morrison HG, Sogin ML. 2010. Ironing out the wrinkles in the rare biosphere through improved OTU clustering. Environ. Microbiol. in press. 10.1111/j.1462-2920.2010.02193.x Web of Science®Google Scholar Jeon S-O, Bunge J, Stoeck T, Barger KJ-A, Hong S-H, Epstein, SS. 2006. Synthetic statistical approach reveals a high degree of richness of microbial eukaryotes in an anoxic water column. Appl. Environ. Microbiol. 72: 6578–6583. 10.1128/AEM.00787-06 CASPubMedWeb of Science®Google Scholar Kemp PF, Aller, JY. 2003. Bacterial diversity in aquatic and other environments: What 16S rDNA libraries can tell us. GEMS Microbial. Ecol. 47: 161–177. 10.1016/S0168-6496(03)00257-5 CASWeb of Science®Google Scholar Knittel K, Losekann T, Boetius A, Kort R, Amann R. 2005. Diversity and distribution of methanotrophic archaea at cold seeps. Appl. Environ. Microbiol. 71: 467–479. 10.1128/AEM.71.1.467-479.2005 CASPubMedWeb of Science®Google Scholar Kunin V, Engelbrektson A, Ochman H, Hugenholtz P. 2010. Wrinkles in the rare biosphere: Pyrosequencing errors lead to artificial inflation of diversity estimates. Environ. Microbiol. 12: 118–123. 10.1111/j.1462-2920.2009.02051.x CASPubMedWeb of Science®Google Scholar Kysela DT, Palacios C, Sogin, M. 2005. Serial analysis of V6 ribosomal sequence tags (SARST-V6): A method for efficient, high-throughput analysis of microbial community composition. Environ. Microbiol. 7: 356–364. 10.1111/j.1462-2920.2004.00712.x CASPubMedWeb of Science®Google Scholar Lozupone C, Knight R. 2005. UniFrac: A new phylogenetic method for comparing microbial communities. Appl. Environ. Microbiol. 71: 8228–8235. 10.1128/AEM.71.12.8228-8235.2005 CASPubMedWeb of Science®Google Scholar Margulies M, Egholm M, Altman WE, Attiya S, Bader JS, Bemben LA, et al. 2005. Genome sequencing in microfabricated high-density picolitre reactors. Nature 437: 376–380. 10.1038/nature03959 CASPubMedWeb of Science®Google Scholar Pace NR. 1997. A molecular view of microbial diversity and the biosphere. Science 276: 734–740. 10.1126/science.276.5313.734 CASPubMedWeb of Science®Google Scholar Pedrós-Alió C. 2006. Marine microbial diversity: Can it be determined? Trends Microbiol. 14: 257–263. 10.1016/j.tim.2006.04.007 CASPubMedWeb of Science®Google Scholar Pedrós-Alió C. 2007. ECOLOGY: Dipping into the Rare Biosphere. Science 315: 192–193. 10.1126/science.1135933 CASPubMedWeb of Science®Google Scholar Porter KG, Feig YS. 1980. The use of DAPI for identifying and counting aquatic microflora Limnol. Oceanogr. 25: 943–948. 10.4319/lo.1980.25.5.0943 Web of Science®Google Scholar Pruesse E, Quast C, Knittel K, Fuchs BM, Ludwig W, Peplies J, Glockner, FO. 2007. SILVA: A comprehensive online resource for quality checked and aligned ribosomal RNA sequence data compatible with ARB. Nucleic Acids Res. 35: 7188–7196. 10.1093/nar/gkm864 CASPubMedWeb of Science®Google Scholar Quince C, Curtis TP, Sloan WT. 2008. The rational exploration of microbial diversity. ISME J. 2: 997–1006. 10.1038/ismej.2008.69 CASPubMedWeb of Science®Google Scholar Rappe MS, Giovannoni, SJ. 2003. The uncultured microbial majority. Annu. Rev. Microbiol. 57: 369–394. 10.1146/annurev.micro.57.030502.090759 CASPubMedWeb of Science®Google Scholar Ravenschlag K, Sahm K, Pernthaler J, Amann, R. 1999. High bacterial diversity in permanently cold marine sediments. Appl. Environ. Microbiol. 65: 3982–3989. 10.1128/AEM.65.9.3982-3989.1999 CASPubMedWeb of Science®Google Scholar Reysenbach AL, Longnecker K, Kirshtein J. 2000. Novel bacterial and archaeal lineages from an in situ growth chamber deployed at a Mid-Atlantic Ridge hydrothermal vent. Appl. Environ. Microbiol. 66: 3798–3806. 10.1128/AEM.66.9.3798-3806.2000 CASPubMedWeb of Science®Google Scholar Roux FL, Zouine M, Chakroun N, Binesse J, Saulnie, D, Bouchier C, et al. 2009. Genome sequence of Vibrio splendidus: An abundant planctonic marine species with a large genotypic diversity. Environ. Microbiol. 11: 1959–1970. 10.1111/j.1462-2920.2009.01918.x CASPubMedWeb of Science®Google Scholar Schloss PD, Handelsman, J. 2005. Introducing DOTUR, a computer program for defining operational taxonomic units and estimating species richness. Appl. Environ. Microbiol. 71: 1501–1506. 10.1128/AEM.71.3.1501-1506.2005 CASPubMedWeb of Science®Google Scholar Schloss PD, Westcott SL, Ryabin T, Hall JR, Hartmann M, Hollister EB, et al. 2009. Introducing mothur: Open source, platform-independent, community-supported software for describing and comparing microbial communities. Appl. Environ. Microbiol. 75: 7537–7541. 10.1128/AEM.01541-09 CASPubMedWeb of Science®Google Scholar Sloan WT, Woodcock S, Lunn M, Head IM, Curtis TP. 2007. Modeling taxa-abundance distributions in microbial communities using environmental sequence data. Microb. Ecol. 53: 443–455. 10.1007/s00248-006-9141-x PubMedWeb of Science®Google Scholar Sogin ML, Morrison HG, Huber JA, Mark Welch DB, Huse SM, Neal PR, et al. 2006. Microbial diversity in the deep sea and the underexplored “rare biosphere”. Proc. Natl. Acad. Sci. USA 103: 12115–12120. 10.1073/pnas.0605127103 CASPubMedWeb of Science®Google Scholar Takai K, Campbell BJ, Cary SC, Suzuki M, Oida H, Nunoura T, et al. 2005. Enzymatic and genetic characterization of carbon and energy metabolisms by deep-sea hydrothermal chemolithoautotrophic isolates of Epsilonproteobacteria. Appl. Environ. Microbiol. 71: 7310–7320. 10.1128/AEM.71.11.7310-7320.2005 CASPubMedWeb of Science®Google Scholar Tringe SG, von Mering C, Kobayashi A, Salamov AA, Chen K, Chang HW, et al. 2005. Comparative metagenomics of microbial communities. Science 308: 554–557. 10.1126/science.1107851 CASPubMedWeb of Science®Google Scholar Venter JC, Remington K, Heidelberg JF, Halpern AL, Rusch D, Eisen JA, et al. 2004. Environmental genome shotgun sequencing of the Sargasso Sea. Science 304: 66–74. 10.1126/science.1093857 CASPubMedWeb of Science®Google Scholar Whitman WB, Coleman DC, Wiebe WJ. 1998. Prokaryotes: The unseen majority. Proc. Natl. Acad. Sci. USA 95: 6578–6583. 10.1073/pnas.95.12.6578 CASPubMedWeb of Science®Google Scholar Citing Literature Handbook of Molecular Microbial Ecology II: Metagenomics in Different Habitats ReferencesRelatedInformation

Three-dimensional (3D) cultures are increasing in use because of their ability to represent in vivo human physiology when compared to monolayer two-dimensional (2D) cultures. When grown in 3D using scaffold-free agarose hydrogels, MCF-7 human breast cancer cells self-organize to form directionally-oriented microtissues that contain a luminal space, reminiscent of the in vivo structure of the mammary gland. When compared to MCF-7 cells cultured in 2D monolayer culture, MCF-7 microtissues exhibit increased mRNA expression of luminal epithelial markers keratin 8 and keratin 19 and decreased expression of basal marker keratin 14 and the mesenchymal marker vimentin. These 3D MCF-7 microtissues remain responsive to estrogens, as demonstrated by induction of known estrogen target mRNAs following exposure to 17β-estradiol. Culture of MCF-7 cells in scaffold-free conditions allows for the formation of more differentiated, estrogen-responsive structures that are a more relevant system for evaluation of estrogenic compounds than traditional 2D models.

Microorganisms associated with coastal sands serve as a natural biofilter, providing essential nutrient recycling in nearshore environments and acting to maintain coastal ecosystem health. Anthropogenic stressors often impact these ecosystems, but little is known about whether these disturbances can be identified through microbial community change. The blowout of the Macondo Prospect reservoir on April 20, 2010, which released oil hydrocarbons into the Gulf of Mexico, presented an opportunity to examine whether microbial community composition might provide a sensitive measure of ecosystem disturbance. Samples were collected on four occasions, beginning in mid-June, during initial beach oiling, until mid-November from surface sand and surf zone waters at seven beaches stretching from Bay St. Louis, MS to St. George Island, FL USA. Oil hydrocarbon measurements and NOAA shoreline assessments indicated little to no impact on the two most eastern beaches (controls). Sequence comparisons of bacterial ribosomal RNA gene hypervariable regions isolated from beach sands located to the east and west of Mobile Bay in Alabama demonstrated that regional drivers account for markedly different bacterial communities. Individual beaches had unique community signatures that persisted over time and exhibited spatial relationships, where community similarity decreased as horizontal distance between samples increased from one to hundreds of meters. In contrast, sequence analyses detected larger temporal and less spatial variation among the water samples. Superimposed upon these beach community distance and time relationships, was increased variability in bacterial community composition from oil hydrocarbon contaminated sands. The increased variability was observed among the core, resident, and transient community members, indicating the occurrence of community-wide impacts rather than solely an overprinting of oil hydrocarbon-degrading bacteria onto otherwise relatively stable sand population structures. Among sequences classified to genus, Alcanivorax, Alteromonas, Marinobacter, Winogradskyella, and Zeaxanthinibacter exhibited the largest relative abundance increases in oiled sands.

[This corrects the article DOI: 10.1186/s40168-015-0081-x.].

Spatial patterns of marine Synechococcus diversity across ocean domains have been reported on extensively.However, much less is known of seasonal and multiannual patterns of change in Synechococcus community composition.Here we report on the genotypic diversity of Synechococcus populations in the Gulf of Aqaba, Northern Red Sea, over seven annual cycles of deep mixing and stabile stratification, using ntcA as a phylogenetic marker.Synechococcus clone libraries were dominated by clade II and XII genotypes and a total of eight different clades were identified.Inclusion of ntcA sequences from the Global Ocean Sampling database in our analyses identified members of clade XII from beyond the Gulf of Aqaba, extending its known distribution.Most of the Synechococcus diversity was attributed to members of clade II during the spring bloom, while clade III contributed significantly to diversity during summer stratification.Clade XII diversity was most prevalent in fall and winter.Clade abundances were estimated from pyrosequencing of the V6 hypervariable region of 16S rRNA.Members of clade II dominated Synechococcus communities throughout the year, whereas the less frequent genotypes showed a pattern of seasonal succession.Based on the prevailing nutritional conditions we observed that clade I members thrive at higher nutrient concentrations during winter mixing.Clades V, VI and X became apparent during the transition periods between mixing and stratification.Clade III became prominent during sumeer stratification.We propose that members of clades V, VI, and X, and clade III are Synechococcus ecotypes that are adapted to intermediate and low nutrient levels respectively.This is the first time that molecular analyses have correlated population dynamics of Synechococcus genotypes with temporal fluctuations in nutrient regimes.Since these Synechococcus genotypes are routinely observed in the Gulf of Aqaba we suggest that seasonal fluctuations in nutrient levels create temporal niches that sustain their coexistence.

BackgroundProspective genetic evaluation of patients at this referral research hospital presents clinical research challenges.ObjectivesThis study sought not only a single-gene explanation for participants’ immune-related presentations, but viewed each participant holistically, with the potential to have multiple genetic contributions to their immune phenotype and other heritable comorbidities relevant to their presentation and health.MethodsThis study developed a program integrating exome sequencing, chromosomal microarray, phenotyping, results return with genetic counseling, and reanalysis in 1505 individuals from 1000 families with suspected or known inborn errors of immunity.ResultsProbands were 50.8% female, 71.5% were ≥18 years, and had diverse immune presentations. Overall, 327 of 1000 probands (32.7%) received 361 molecular diagnoses. These included 17 probands with diagnostic copy number variants, 32 probands with secondary findings, and 31 probands with multiple molecular diagnoses. Reanalysis added 22 molecular diagnoses, predominantly due to new disease-gene associations (9 of 22, 40.9%). One-quarter of the molecular diagnoses (92 of 361) did not involve immune-associated genes. Molecular diagnosis was correlated with younger age, male sex, and a higher number of organ systems involved. This program also facilitated the discovery of new gene-disease associations such as SASH3-related immunodeficiency. A review of treatment options and ClinGen actionability curations suggest that at least 251 of 361 of these molecular diagnoses (69.5%) could translate into ≥1 management option.ConclusionsThis program contributes to our understanding of the diagnostic and clinical utility whole exome analysis on a large scale. Prospective genetic evaluation of patients at this referral research hospital presents clinical research challenges. This study sought not only a single-gene explanation for participants’ immune-related presentations, but viewed each participant holistically, with the potential to have multiple genetic contributions to their immune phenotype and other heritable comorbidities relevant to their presentation and health. This study developed a program integrating exome sequencing, chromosomal microarray, phenotyping, results return with genetic counseling, and reanalysis in 1505 individuals from 1000 families with suspected or known inborn errors of immunity. Probands were 50.8% female, 71.5% were ≥18 years, and had diverse immune presentations. Overall, 327 of 1000 probands (32.7%) received 361 molecular diagnoses. These included 17 probands with diagnostic copy number variants, 32 probands with secondary findings, and 31 probands with multiple molecular diagnoses. Reanalysis added 22 molecular diagnoses, predominantly due to new disease-gene associations (9 of 22, 40.9%). One-quarter of the molecular diagnoses (92 of 361) did not involve immune-associated genes. Molecular diagnosis was correlated with younger age, male sex, and a higher number of organ systems involved. This program also facilitated the discovery of new gene-disease associations such as SASH3-related immunodeficiency. A review of treatment options and ClinGen actionability curations suggest that at least 251 of 361 of these molecular diagnoses (69.5%) could translate into ≥1 management option. This program contributes to our understanding of the diagnostic and clinical utility whole exome analysis on a large scale.

MIxS-BE: a MIxS extension defining a minimum information standard for sequence data from the built environment

Bacterial diversity in contaminated fuels has not been systematically investigated using cultivation-independent methods. The fuel industry relies on phenotypic cultivation-based contaminant identification, which may lack accuracy and neglect difficult-to-culture taxa. By the use of industry practice aerobic cultivation, 16S rRNA gene sequencing, and strain genotyping, a collection of 152 unique contaminant isolates from 54 fuel samples was assembled, and a dominance of Pseudomonas (21%), Burkholderia (7%), and Bacillus (7%) was demonstrated. Denaturing gradient gel electrophoresis (DGGE) of 15 samples revealed Proteobacteria and Firmicutes to be the most abundant phyla. When 16S rRNA V6 gene pyrosequencing of four selected fuel samples (indicated by "JW") was performed, Betaproteobacteria (42.8%) and Gammaproteobacteria (30.6%) formed the largest proportion of reads; the most abundant genera were Marinobacter (15.4%; JW57), Achromobacter (41.6%; JW63), Burkholderia (80.7%; JW76), and Halomonas (66.2%; JW78), all of which were also observed by DGGE. However, the Clostridia (38.5%) and Deltaproteobacteria (11.1%) identified by pyrosequencing in sample JW57 were not observed by DGGE or aerobic culture. Genotyping revealed three instances where identical strains were found: (i) a Pseudomonas sp. strain recovered from 2 different diesel fuel tanks at a single industrial site; (ii) a Mangroveibacter sp. strain isolated from 3 biodiesel tanks at a single refinery site; and (iii) a Burkholderia vietnamiensis strain present in two unrelated automotive diesel samples. Overall, aerobic cultivation of fuel contaminants recovered isolates broadly representative of the phyla and classes present but lacked accuracy by overrepresenting members of certain groups such as Pseudomonas.

In utero exposure to antiandrogenic xenobiotics such as di-n-butyl phthalate (DBP) has been linked to congenital defects of the male reproductive tract, including cryptorchidism and hypospadias, as well as later life effects such as testicular cancer and decreased sperm counts. Experimental evidence indicates that DBP has in utero antiandrogenic effects in the rat. However, it is unclear whether DBP has similar effects on androgen biosynthesis in human fetal testis. To address this issue, we developed a xenograft bioassay with multiple androgen-sensitive physiological endpoints, similar to the rodent Hershberger assay. Adult male athymic nude mice were castrated, and human fetal testis was xenografted into the renal subcapsular space. Hosts were treated with human chorionic gonadotropin for 4 weeks to stimulate testosterone production. During weeks 3 and 4, hosts were exposed to DBP or abiraterone acetate, a CYP17A1 inhibitor. Although abiraterone acetate (14 d, 75mg/kg/d po) dramatically reduced testosterone and the weights of androgen-sensitive host organs, DBP (14 d, 500mg/kg/d po) had no effect on androgenic endpoints. DBP did produce a near-significant trend toward increased multinucleated germ cells in the xenografts. Gene expression analysis showed that abiraterone decreased expression of genes related to transcription and cell differentiation while increasing expression of genes involved in epigenetic control of gene expression. DBP induced expression of oxidative stress response genes and altered expression of actin cytoskeleton genes.

Bisphenol A (BPA) is a ubiquitous industrial chemical that has been identified as an endocrine disrupting compound (EDC). There is growing concern that early life exposures to EDCs, such as BPA, can adversely affect the male reproductive tract and function. This study was conducted as part of the Consortium Linking Academic and Regulatory Insights on BPA Toxicity (CLARITY-BPA) to further delineate the toxicities associated with continuous exposure to BPA from early gestation, and to comprehensively examine the elicited effects on testes and sperm. NCTR Sprague Dawley rat dams were gavaged from gestational day (GD) 6 until parturition, and their pups were directly gavaged daily from postnatal day (PND) 1 to 90 with BPA (2.5, 25, 250, 2500, 25,000, 250,000 μg/kg/d) or vehicle control. At PND 90, the testes and sperm were collected for evaluation. The testes were histologically evaluated for altered germ cell apoptosis, sperm production, and altered spermiation. RNA and DNA isolated from sperm were assessed for elicited changes in global mRNA transcript abundance and altered DNA methylation. Effects of BPA were observed in changes in body, testis and epididymis weights only at the highest administered dose of BPA of 250,000 μg/kg/d. Genome-wide transcriptomic and epigenomic analyses failed to detect robust alterations in sperm mRNA and DNA methylation levels. These data indicate that prolonged exposure starting in utero to BPA over a wide range of levels has little, if any, impact on the testes and sperm molecular profiles of 90 day old rats as assessed by the histopathologic, morphometric, and molecular endpoints evaluated.

The indigenous gut microbiota are thought to play a crucial role in the development and maintenance of the abnormal inflammatory responses that are the hallmark of inflammatory bowel disease. Direct tests of the role of the gut microbiome in these disorders are typically limited by the fact that sampling of the microbiota generally occurs once disease has become manifest. This limitation could potentially be circumvented by studying patients who undergo total proctocolectomy with ileal pouch anal anastomosis (IPAA) for the definitive treatment of ulcerative colitis. A subset of patients who undergo IPAA develops an inflammatory condition known as pouchitis, which is thought to mirror the pathogenesis of ulcerative colitis. Following the development of the microbiome of the pouch would allow characterization of the microbial community that predates the development of overt disease. We monitored the development of the pouch microbiota in four patients who underwent IPAA. Mucosal and luminal samples were obtained prior to takedown of the diverting ileostomy and compared to samples obtained 2, 4 and 8 weeks after intestinal continuity had been restored. Through the combined analysis of 16S rRNA-encoding gene amplicons, targeted 16S amplification and microbial cultivation, we observed major changes in structure and function of the pouch microbiota following ileostomy. There is a relative increase in anaerobic microorganisms with the capacity for fermentation of complex carbohydrates, which corresponds to the physical stasis of intestinal contents in the ileal pouch. Compared to the microbiome structure encountered in the colonic mucosa of healthy individuals, the pouch microbial community in three of the four individuals was quite distinct. In the fourth patient, a community that was much like that seen in a healthy colon was established, and this patient also had the most benign clinical course of the four patients, without the development of pouchitis 2 years after IPAA. The microbiota that inhabit the ileal-anal pouch of patients who undergo IPAA for treatment of ulcerative colitis demonstrate significant structural and functional changes related to the restoration of fecal flow. Our preliminary results suggest once the pouch has assumed the physiologic role previously played by the intact colon, the precise structure and function of the pouch microbiome, relative to a normal colonic microbiota, will determine if there is establishment of a stable, healthy mucosal environment or the reinitiation of the pathogenic cascade that results in intestinal inflammation.

Amplicon sequencing of the hypervariable regions of the small subunit ribosomal RNA gene is a widely accepted method for identifying the members of complex bacterial communities. Several rRNA gene sequence reference databases can be used to assign taxonomic names to the sequencing reads using BLAST, USEARCH, GAST or the RDP classifier. Next-generation sequencing methods produce ample reads, but they are short, currently ∼100-450 nt (depending on the technology), as compared to the full rRNA gene of ∼1550 nt. It is important, therefore, to select the right rRNA gene region for sequencing. The primers should amplify the species of interest and the hypervariable regions should differentiate their taxonomy. Here, we introduce TaxMan: a web-based tool that trims reference sequences based on user-selected primer pairs and returns an assessment of the primer specificity by taxa. It allows interactive plotting of taxa, both amplified and missed in silico by the primers used. Additionally, using the trimmed sequences improves the speed of sequence matching algorithms. The smaller database greatly improves run times (up to 98%) and memory usage, not only of similarity searching (BLAST), but also of chimera checking (UCHIME) and of clustering the reads (UCLUST). TaxMan is available at http://www.ibi.vu.nl/programs/taxmanwww/.

DNA methylation is an important epigenetic control mechanism that has been shown to be associated with gene silencing through the course of development, maturation and aging. However, only limited data are available regarding the relationship between methylation and gene expression in human development. We analyzed the methylome and transcriptome of three human fetal liver samples (gestational age 20–22 weeks) and three adult human liver samples. Genes whose expression differed between fetal and adult numbered 7,673. Adult overexpression was associated with metabolic pathways and, in particular, cytochrome P450 enzymes while fetal overexpression reflected enrichment for DNA replication and repair. Analysis for DNA methylation using the Illumina Infinium 450 K HumanMethylation BeadChip showed that 42 % of the quality filtered 426,154 methylation sites differed significantly between adult and fetal tissue (q ≤ 0.05). Differences were small; 69 % of the significant sites differed in their mean methylation beta value by ≤0.2. There was a trend among all sites toward higher methylation in the adult samples with the most frequent difference in beta being 0.1. Characterization of the relationship between methylation and expression revealed a clear difference between fetus and adult. Methylation of genes overexpressed in fetal liver showed the same pattern as seen for genes that were similarly expressed in fetal and adult liver. In contrast, adult overexpressed genes showed fetal hypermethylation that differed from the similarly expressed genes. An examination of gene region-specific methylation showed that sites proximal to the transcription start site or within the first exon with a significant fetal-adult difference in beta (>0.2) showed an inverse relationship with gene expression. Nearly half of the CpGs in human liver show a significant difference in methylation comparing fetal and adult samples. Sites proximal to the transcription start site or within the first exon that show a transition from hypermethylation in the fetus to hypomethylation or intermediate methylation in the adult are associated with inverse changes in gene expression. In contrast, increases in methylation going from fetal to adult are not associated with fetal-to-adult decreased expression. These findings indicate fundamentally different roles for and/or regulation of DNA methylation in human fetal and adult liver.

Prostate cancer is the most frequent non-cutaneous malignancy in men. There is strong evidence in rodents that neonatal estrogen exposure plays a role in the development of this disease. However, there is little information regarding the effects of estrogen in human fetal prostate tissue. This study explored early life estrogen exposure, with and without a secondary estrogen and testosterone treatment in a human fetal prostate xenograft model. Histopathological lesions, proliferation, and serum hormone levels were evaluated at 7, 30, 90, and 200-day time-points after xenografting. The expression of 40 key genes involved in prostatic glandular and stromal growth, cell-cycle progression, apoptosis, hormone receptors and tumor suppressors was evaluated using a custom PCR array. Epigenome-wide analysis of DNA methylation was performed on whole tissue, and laser capture-microdissection (LCM) isolated epithelial and stromal compartments of 200-day prostate xenografts. Combined initial plus secondary estrogenic exposures had the most severe tissue changes as revealed by the presence of hyperplastic glands at day 200. Gene expression changes corresponded with the cellular events in the KEGG prostate cancer pathway, indicating that initial plus secondary exposure to estrogen altered the PI3K-Akt signaling pathway, ultimately resulting in apoptosis inhibition and an increase in cell cycle progression. DNA methylation revealed that differentially methylated CpG sites significantly predominate in the stromal compartment as a result of estrogen-treatment, thereby providing new targets for future investigation. By using human fetal prostate tissue and eliminating the need for species extrapolation, this study provides novel insights into the gene expression and epigenetic effects related to prostate carcinogenesis following early life estrogen exposure.

There is growing evidence that sperm DNA methylation is important in maintaining proper sperm health and function. Previous studies have associated sperm DNA methylation levels with sperm quality and function, however, little is known regarding the intra- and inter-individual variability in sperm methylation levels. This study characterizes this variation. Sperm epigenetic differences between successive semen samples from 12 patients were examined to identify the intra- and inter-individual differences globally across the genome, and in specifically defined genomic regions using the Illumina Infinium HumanMethylation450 BeadChips. Methylation analysis identified a bimodal distribution in the methylation levels that were non-uniformly distributed across the different genomic regions. The methylation levels were highly correlated in both the intra- and inter-individual comparisons. The intra-individual methylation levels were more highly correlated than the inter-individual comparison both globally and across the defined genomic regions, demonstrating that sperm DNA methylation levels are relatively stable between semen sample collections.

Exposure to the herbicide Agent Orange during the Vietnam War was widespread and is associated with numerous adverse health outcomes. A continuing concern of veterans is the possibility that exposure to the dioxin-containing herbicide might induce adverse reproductive outcomes. We sought to assess whether exposure to Agent Orange in Vietnam was associated with changes in DNA methylation in sperm in a subset of Vietnam veterans who participated in the Air Force Health Study (AFHS).We studied 37 members of the AFHS chosen to have no, low, medium or high exposure to Agent Orange, based upon serum dioxin levels obtained during a series of examinations. DNA from stored semen was extracted and DNA methylation assessed on the Illumina 450 K platform.Initial epigenome-wide analysis returned no loci that survived control for false discovery. However, the TEAD3 gene had four different CpG sites that showed loss of DNA methylation associated with dioxin exposure. Analysis assessing regional DNA methylation changes revealed 36 gene regions, including the region of the imprinted gene H19 to have altered DNA methylation associated with high exposure compared to the low exposure group. Additional comparison of our data with sperm DNA methylation data from Russian boys exposed to dioxin found an additional 5 loci that were altered in both studies and exhibited a consistent direction of association.Studying a small number of sperm samples from veterans enrolled in the AFHS, we did not find evidence of significant epigenome-wide alterations associated with exposure to Agent Orange. However, additional analysis showed that the H19 gene region is altered in the sperm of Agent Orange-exposed Ranch Hand veterans. Our study also replicated several findings of a prior study of dioxin-exposed Russian boys. These results provide additional candidate loci for further investigation and may have implications for the reproductive health of dioxin-exposed individuals.

Prostate cancer is the most commonly diagnosed nonskin cancer in men. The etiology of prostate cancer is unknown, although both animal and epidemiologic data suggest that early life exposures to various toxicants, may impact DNA methylation status during development, playing an important role.We have developed a xenograft model to characterize the growth and differentiation of human fetal prostate implants (gestational age 12-24 weeks) that can provide new data on the potential role of early life stressors on prostate cancer. The expression of key immunohistochemical markers responsible for prostate maturation was evaluated, including p63, cytokeratin 18, α-smooth muscle actin, vimentin, caldesmon, Ki-67, prostate-specific antigen, estrogen receptor-α, and androgen receptor. Xenografts were separated into epithelial and stromal compartments using laser capture microdissection (LCM), and the DNA methylation status was assessed in >480,000 CpG sites throughout the genome.Xenografts demonstrated growth and maturation throughout the 200 days of post-implantation evaluation. DNA methylation profiles of laser capture microdissected tissue demonstrated tissue-specific markers clustered by their location in either the epithelium or stroma of human prostate tissue. Differential methylated promoter region CpG-associated gene analysis revealed significantly more stromal than epithelial DNA methylation in the 30- and 90-day xenografts. Functional classification analysis identified CpG-related gene clusters in methylated epithelial and stromal human xenografts.This study of human fetal prostate tissue establishes a xenograft model that demonstrates dynamic growth and maturation, allowing for future mechanistic studies of the developmental origins of later life proliferative prostate disease.

In the article ‘Diversity and population structure of sewage-derived microorganisms in wastewater treatment plant influent’ (Mclellan et al., 2010), with reference to Table 2, we would like to clarify the following points: In column 4, the heading should be as follows: No. of sewage OTUs shared with humans. In column 5, the heading should be as follows: % of sewage tag data set shared with humans. In column 7, the heading should be as follows: % of sewage tag data set shared with surface water. We apologize for this mistake.

Abstract Marine bacterial diversity is vast, but seasonal variation in diversity is poorly understood. Here we present the longest bacterial diversity time series consisting of monthly (72) samples from the western English Channel over a 6 year period (2003-2008) using 747,494 16SrDNA-V6 amplicon-pyrosequences. Although there were characteristic cycles for each phylum, the overall community cycle was remarkably stable year after year. The majority of taxa were not abundant, although on occasion these rare bacteria could dominate the assemblage. Bacterial diversity peaked at the winter solstice and showed remarkable synchronicity with day-length, which had the best explanatory power compared to a combination of other variables (including temperature and nutrient concentrations). Day-length has not previously been recognised as a major force in structuring microbial communities.

In the development of human cell-based assays, 3-dimensional (3D) cell culture models are intriguing as they are able to bridge the gap between animal models and traditional two-dimensional (2D) cell culture. Previous work has demonstrated that MCF-7 human breast carcinoma cells cultured in a 3D scaffold-free culture system self-assemble and develop into differentiated microtissues that possess a luminal space. Exposure to estradiol for 7 days decreased lumen formation in MCF-7 microtissues, altered microtissue morphology and altered expression of genes involved in estrogen signaling, cell adhesion and cell cycle regulation. Exposure to receptor-specific agonists for estrogen receptor alpha, estrogen receptor beta and g-protein coupled estrogen receptor resulted in unique, receptor-specific phenotypes and gene expression signatures. The use of a differentiated scaffold-free 3D culture system offers a unique opportunity to study the phenotypic and molecular changes associated with exposure to estrogenic compounds.

The extremely high error rates reported by Keegan et al. in 'A platform-independent method for detecting errors in metagenomic sequencing data: DRISEE' (PLoS Comput Biol 2012; 8: :e1002541) for many next-generation sequencing datasets prompted us to re-examine their results. Our analysis reveals that the presence of conserved artificial sequences, e.g. Illumina adapters, and other naturally occurring sequence motifs accounts for most of the reported errors. We conclude that DRISEE reports inflated levels of sequencing error, particularly for Illumina data. Tools offered for evaluating large datasets need scrupulous review before they are implemented.

Over the past decade, laser capture microdissection (LCM) has grown as a tool for gene expression profiling of small numbers of cells from tumor samples and of specific cell populations in complex tissues. LCM can be used to study toxicant effects on selected cell populations within the testis at different stages of spermatogenesis. There are several LCM-related hurdles to overcome, including issues inherent to the method itself, as well as biases that result from amplifying the LCM-isolated RNA. Many technical issues associated with the LCM method are addressed here, including increasing RNA yield and obtaining more accurate quantification of RNA yields. We optimized the LCM method optimized to generate RNA quantities sufficient for quantitative reverse transcription polymerase chain reaction (qRT-PCR) array analysis without amplification and were able to validate the method through direct comparison of results from unamplified and amplified RNA from individual samples. The addition of an amplification step for gene expression studies using LCM RNA resulted in a bias, especially for low abundance transcripts. Although the amplification bias was consistent across samples, researchers should use caution when comparing results generated from amplified and unamplified LCM RNA. Here, we have validated the use of LCM-derived RNA with the qRT-PCR array, improving our ability to investigate cell-type and stage-specific responses to toxicant exposures.

Purified brush-border membrane vesicles (BBMV) of lobster antennal gland labyrinth and bladder were separately formed by a magnesium precipitation technique. L-[3H]proline uptake was stimulated by a transmembrane NaCl gradient [outside (o) greater than inside (i)] to a greater extent in BBMV from labyrinth than those from the bladder. Detailed study of the labyrinth proline-transport processes revealed a specific dependence on NaCl, with negligible stimulatory effects by NaSCN, Na-gluconate, or KCl. A transmembrane proton gradient (o greater than i) was without stimulatory effect on proline transport. A transmembrane potential difference alone, in the presence of equilibrated NaCl and L-[3H]proline, led to net influx of the labeled amino acid, suggesting that the uptake process was electrogenic and capable of bringing about the net transfer of positive charge to the vesicle interior. Although a transmembrane Na gradient alone, in the presence of equilibrated Cl and L-[3H]proline, was able to bring about the net influx of the amino acid, a transmembrane Cl gradient alone under Na- and L-[3H]proline-equilibrated conditions was not, suggesting that only the Na gradient could energize the carrier process through cotransport, while the anion served an essential activating role. Proline influx by these vesicles occurred by the combination of at least one saturable Michaelis-Menten carrier system (apparent Kt = 0.37 mM; apparent JM = 1.19 nmol.mg protein-1.10 s-1) and apparent diffusion (P = 0.33 nmol.mg protein-1.10 s-1.mM-1). Static head analysis of the transport process suggested a cotransport stoichiometry of 2 Na:1 proline with essential activation by Cl ion.

Abstract The nation's estuaries are at risk of further deterioration from land use change and intensification. These risks include direct impacts on wetland habitats and stream environments and indirect impacts from nonpoint source pollutant loading. This article reports on the methods, findings, and policy implications of a major study "The Effects of Land Use Change and Intensification on the San Francisco Estuary.”; By using a geographic information system (GIS), future growth scenarios were played out and the impacts on wetlands, streams, and water quality were estimated on a regionwide basis. The land use scenario developed from the General Plans of the Bay‐Delta Region's 12 counties shows that the total area planned as urban use outside existing incorporated cities is 331,530 acres, an increase of 37%. This land use change and intensification associated with increased growth will continue to stress an overtaxed estuarine system. Results are expressed according to 14 receiving water segments and the associated 34 watersheds. Direct impacts on wetlands and stream environment zones occur in every watershed containing these resources. We estimate that over 39,500 acres of wetlands may be potentially impacted. Of 377,000 acres of stream environment in the 12‐county study area, 28,000 acres are also subject to impacts of urbanization. Our analysis suggests that protection of farmed wetlands in the Delta and North Bay and the enhancement of biodiversity in the South Bay deserve special attention. The construction of land use scenarios for the estuary region has presented, for the first time, an opportunity to examine the cumulative contribution of nonpoint source urban runoff to the levels of pollutants in the bay and delta. To date, more modest studies in smaller urban watersheds have provided only a glimpse of the overall effect that urbanization has in a region the size of the estuary. We found that these impacts can be expected to decrease the overall water quality of the estuary. The existing system of land use planning delegates responsibility to local governments. However, of 111 jurisdictions within the estuary study region, only 18 have specific ordinances to protect streams and wetlands. We recommend that the existing system regulation and management be strengthened to protect, enhance, and restore the environmental well‐being of the estuary. The results of our study suggest that improvements are needed in the goals, management strategies, and institutional arrangements now in place for the San Francisco estuary. We advance several specific management options to help frame the debate over strengthening estuarine management. In particular, we urge that a specific focus on estuarine resource protection be incorporated in any new growth management legislation enacted in California. We identify several important national policy implications arising from our study. First, we believe the potential transferability of our methodology to other estuaries should be investigated. Second, we recommend that technical workshops be convened for estuarine managers who are addressing similar management issues. Third, we recommend that policy guidance to encourage the use of watersheds and receiving waters as the unit for analysis. We also recommend that GIS based analysis should be used to test the implications of alternative wetland definitions to inform the national policy debate.

Testicular effects of chemical mixtures may differ from those of the individual chemical constituents. This study assessed the co-exposure effects of the model germ cell- and Sertoli cell-specific toxicants, X-irradiation (x-ray), and 2,5-hexanedione (HD), respectively. In high-dose studies, HD has been shown to attenuate x-ray-induced germ cell apoptosis. Adult rats were exposed to different levels of x-ray (0.5 Gy, 1 Gy, and 2 Gy) or HD (0.33%), either alone or in combination. To assess cell type-specific attenuation of x-ray effects with HD co-exposure, we used laser capture microdissection (LCM) to enrich the targeted cell population and examine a panel of apoptosis-related transcripts using PCR arrays. The apoptosis PCR arrays identified significant dose-dependent treatment effects on several genes, with downregulation of death receptor 5 ( DR5), Naip2, Sphk2, Casp7, Aven, Birc3, and upregulation of Fas. The greatest difference in transcript response to exposure was seen with 0.5 Gy x-ray exposure, and the attenuation effect seen with the combined high-dose x-ray and HD did not persist into the low-dose range. Examination of protein levels in staged tubules revealed a significant upregulation in DR5, following high-dose co-exposure. These results provide insight into the testis cell-specific apoptotic response to low-dose co-exposures of model testicular toxicants.

We present a performance-optimized algorithm, subsampled open-reference OTU picking, for assigning marker gene (e.g., 16S rRNA) sequences generated on next-generation sequencing platforms to operational taxonomic units (OTUs) for microbial community analysis. This algorithm provides benefits over de novo OTU picking (clustering can be performed largely in parallel, reducing runtime) and closed-reference OTU picking (all reads are clustered, not only those that match a reference database sequence with high similarity). Because more of our algorithm can be run in parallel relative to “classic” open-reference OTU picking, it makes open-reference OTU picking tractable on massive amplicon sequence data sets (though on smaller data sets, “classic” open-reference OTU clustering is often faster). We illustrate that here by applying it to the first 15,000 samples sequenced for the Earth Microbiome Project (1.3 billion V4 16S rRNA amplicons). To the best of our knowledge, this is the largest OTU picking run ever performed, and we estimate that our new algorithm runs in less than 1/5 the time than would be required of “classic” open reference OTU picking. We show that subsampled open-reference OTU picking yields results that are highly correlated with those generated by “classic” open-reference OTU picking through comparisons on three well-studied datasets. An implementation of this algorithm is provided in the popular QIIME software package, which uses uclust for read clustering. All analyses were performed using QIIME’s uclust wrappers, though we provide details (aided by the open-source code in our GitHub repository) that will allow implementation of subsampled open-reference OTU picking independently of QIIME (e.g., in a compiled programming language, where runtimes should be further reduced). Our analyses should generalize to other implementations of these OTU picking algorithms. Finally, we present a comparison of parameter settings in QIIME’s OTU picking workflows and make recommendations on settings for these free parameters to optimize runtime without reducing the quality of the results. These optimized parameters can vastly decrease the runtime of uclust-based OTU picking in QIIME.

Long-term heavy alcohol abuse causes progressive degenerative diseases in liver and brain because of insulin resistance mediated by impairments in signaling at multiple levels within the insulin-signaling cascade. Alcohol-mediated insulin resistance dysregulates lipid metabolism, causing states of lipotoxicity, impairs growth and mitochondrial function, and increases oxidative and endoplasmic reticulum stress. Together, these effects worsen injury, inflammation, and insulin resistance, establishing a self-reinforcing vicious cycle. Alcohol-mediated degeneration of liver and brain are abetted by cofactors including nutritional deficiencies, particularly thiamine, tobacco smoke/tobacco nitrosamine exposures, and increased intestinal permeability coupled with pathogenic shifts in the gut microbiota. Furthermore, alcohol-related neurodegeneration is likely exacerbated by indirect mechanisms linked to steatohepatitis (i.e., liver-derived toxic lipids crossing the blood–brain barrier and damaging the brain).

We present a performance-optimized algorithm, subsampled open-reference OTU picking, for assigning marker gene (e.g., 16S rRNA) sequences generated on next-generation sequencing platforms to operational taxonomic units (OTUs) for microbial community analysis. This algorithm provides benefits over de novo OTU picking (clustering can be performed largely in parallel, reducing runtime) and closed-reference OTU picking (all reads are clustered, not only those that match a reference database sequence with high similarity). Because parts of our algorithm can be run in parallel, it makes open-reference OTU picking tractable on massive amplicon sequence data sets. We illustrate that here by applying it to the first 15,000 samples sequenced for the Earth Microbiome Project (1.3 billion V4 16S rRNA amplicons). To the best of our knowledge, this is the largest OTU picking run ever performed. We show that subsampled open-reference OTU picking yields results that are highly correlated with those generated by “legacy” open-reference OTU picking, where less of the process can be parallelized, through comparisons on three well-studied datasets. We therefore recommend that subsampled open-reference OTU picking always be applied in favor of “legacy” open-reference OTU picking. An implementation of this algorithm is provided in the popular QIIME software package. Finally, we present a comparison of parameter settings in QIIME’s OTU picking workflows and make recommendations on settings for these free parameters.

We present a performance-optimized algorithm, subsampled open-reference OTU picking, for assigning marker gene (e.g., 16S rRNA) sequences generated on next-generation sequencing platforms to operational taxonomic units (OTUs) for microbial community analysis. This algorithm provides benefits over de novo OTU picking (clustering can be performed largely in parallel, reducing runtime) and closed-reference OTU picking (all reads are clustered, not only those that match a reference database sequence with high similarity). Because more of our algorithm can be run in parallel relative to “classic” open-reference OTU picking, it makes open-reference OTU picking tractable on massive amplicon sequence data sets (though on smaller data sets, “classic” open-reference OTU clustering is often faster). We illustrate that here by applying it to the first 15,000 samples sequenced for the Earth Microbiome Project (1.3 billion V4 16S rRNA amplicons). To the best of our knowledge, this is the largest OTU picking run ever performed, and we estimate that our new algorithm runs in less than 1/5 the time than would be required of “classic” open reference OTU picking. We show that subsampled open-reference OTU picking yields results that are highly correlated with those generated by “classic” open-reference OTU picking through comparisons on three well-studied datasets. An implementation of this algorithm is provided in the popular QIIME software package, which uses uclust for read clustering. All analyses were performed using QIIME’s uclust wrappers, though we provide details (aided by the open-source code in our GitHub repository) that will allow implementation of subsampled open-reference OTU picking independently of QIIME (e.g., in a compiled programming language, where runtimes should be further reduced). Our analyses should generalize to other implementations of these OTU picking algorithms. Finally, we present a comparison of parameter settings in QIIME’s OTU picking workflows and make recommendations on settings for these free parameters to optimize runtime without reducing the quality of the results. These optimized parameters can vastly decrease the runtime of uclust-based OTU picking in QIIME.

RDATAand RMD.tar.gz is a zipped file of all the data and R markdown files to reproduce all the analyses done in the article and is the files can be found at : &lt; http://purl.stanford.edu/fs506ff9976 &gt;. (ZIP 11504Â kb)

Coyote Slough Lagoon was flooded as part of a local mitigation project in 1986. The design of the restored marsh included peninsulas to break up the local fetch and resulted in the development of 5 distinct tidal cells. Topographics maps were made just before and after flooding. Subsequent repeat bathymetric surveying and deployment of sedimentation plates allow us to examine the spatial and temporal variability of sedimentation within the lagoon. Below mean high water repeat surveying was used to determine sedimentation rates. In the 6 years following flooding sedimentation has been rapid in lagoon cell close to the main tidal inlet and a distinct channel has formed there. Sedimentation rates in the subtidal areas vary from greater than 60 cm/yr. near the main inlet to a minimum of 10 cm/yr. in the most remote area of the lagoon. On the mudflats above mean low water sediment plates were deployed over a two week period to monitor short term sedimentation rates. On these mudflats we calculate sedimentation rates varying from 5.67 ch/ yr in the central part of the lagoon to 1.64 cm/yr inn the most remote parts of the lagoon. Sedimentation rate varies with both distance from the tidal inlet and elevation. The deepest areas of the lagoon have the most rapid rates of sedimentation while the mudflats surrounding the lagoon are agrading more slowly.