Strain level and comprehensive microbiome analysis in inflammatory bowel disease via multi-technology meta-analysis identifies key bacterial influencers of disease.

Objective: Inflammatory bowel disease (IBD) is a heterogenous disease in which the microbiome has been shown to play an important role. However, the precise homeostatic or pathological functions played by bacteria remain unclear. Most published studies report taxa-disease associations based on single-technology analysis of a single cohort, potentially biasing results to one clinical protocol, cohort, and molecular analysis technology. To begin to address this key question, precise identification of the bacteria implicated in IBD across cohorts is necessary. Methods: We sought to take advantage of the numerous and diverse studies characterizing the microbiome in IBD to develop a multi-technology meta-analysis (MTMA) as a platform for aggregation of independently generated datasets, irrespective of DNA-profiling technique, in order to uncover the consistent microbial modulators of disease. We report the largest strain-level survey of IBD, integrating microbiome profiles from 3,407 samples from 21 datasets spanning 15 cohorts, three of which are presented for the first time in the current study, characterized using three DNA-profiling technologies, mapping all nucleotide data against known, culturable strain reference data. Results: We identify several novel IBD associations with culturable strains that have so far remained elusive, including two genome-sequenced but uncharacterized Lachnospiraceae strains consistently decreased in both the gut luminal and mucosal contents of patients with IBD, and demonstrate that these strains are correlated with inflammation-related pathways that are known mechanisms targeted for treatment. Furthermore, comparative MTMA at the species versus strain level reveals that not all significant strain associations resulted in a corresponding species-level significance and conversely significant species associations are not always re-captured at the strain level. Conclusion: We propose MTMA for uncovering experimentally testable strain-disease associations that, as demonstrated here, are beneficial in discovering mechanisms underpinning microbiome impact on disease or novel targets for therapeutic interventions.

Some Gammaproteobacteria are enriched within CD14+ macrophages from intestinal lamina propria of Crohn's disease patients versus mucus.

Crohn's disease causes chronic inflammation in the gastrointestinal tract and its pathogenesis remains unclear. In the intestine of Crohn's disease patients, CD14+CD11+CD163low macrophages contribute to inflammation through the induction of Th17 cells and production of inflammatory cytokines; the CD14+CD11c+163high fraction is anti-inflammatory through the production of IL-10 in normal cases. In this report, the 16S rRNA gene amplicon sequencing method was used to identify bacteria that are specifically present in intestinal CD14+CD11c+ macrophages of Crohn's disease patients. Bacteria present in intestinal CD14+CD11c+ macrophages and mucus of Crohn's disease patients were separated into different clusters in principal coordinates analysis. There was a statistically significant increase in the relative composition of CD14+CD11c+ macrophages from mucus in two phyla (Proteobacteria [p = 0.01] and Actinobacteria [p = 0.02]) and two families (Moraxellaceae [p < 0.001] and Pseudomonadaceae [p = 0.01]). In addition, OTU-1: Acinetobacter and OTU-8: Pseudomonadaceae tended to concentrate in the CD14+CD11c+CD163low subset, whereas OTU-10: Proteus, OTU-15: Collinsella tended to concentrate more in the CD14+CD11c+CD163high subset than the other subset and mucus.

Impact of a probiotic product on bowel habits and microbial profile in participants with functional constipation: A randomized controlled trial.

OBJECTIVE: To investigate the clinical efficacy of a multi-strain probiotic product on bowel habits and microbial profile in participants with functional constipation. METHODS: This was a randomized, double-blind, placebo-controlled and parallel-arm study. Altogether 94 otherwise healthy adults aged 18 to 65 years with symptoms of functional constipation were randomized as part of the intention-to-treat population. The participants received a placebo or the probiotic product (1.5 × 1010 CFU/day), consisting of Lactobacillus acidophilus DDS-1, Bifidobacterium animalis subsp. lactis UABla-12, Bifidobacterium longum UABl-14 and Bifidobacterium bifidum UABb-10 over 4 weeks. Outcomes included the patient assessment of constipation-symptom (PAC-SYM) questionnaire, stool frequency and consistency, and microbial profile. RESULTS: There were no significant between-group differences in the PAC-SYM score, despite significant within-group differences (P < 0.001) over the study period. The probiotic group showed a faster normalization of stool frequency and consistency, with most participants achieving a normalized profile after 1 week. Fecal samples of the probiotic group exhibited higher relative abundance of Ruminococcaceae (P = 0.0047), including the Ruminococcus genus, and lower relative abundance of Erysipelotrichaceae (P = 0.0172) at end-point compared with baseline. Placebo group samples showed similar abundance profiles over the study, with the exception of Clostridiaceae, which was lower at the study end-point (P = 0.0033). Among treated participants, all four probiotic strains were significantly more abundant after the intervention. CONCLUSIONS: No significant differences were observed in symptomology, with both groups showing a more than 20% improvement. However, the probiotic helped modulate bowel function earlier than the placebo, with a corresponding shift to a more fibrolytic microbiota.

The urinary microbiome associated with bladder cancer.

Recent findings suggest that human microbiome can influence the development of cancer, but the role of microorganisms in bladder cancer pathogenesis has not been explored yet. The aim of this study was to characterize and compare the urinary microbiome of bladder cancer patients with those of healthy controls. Bacterial communities present in urine specimens collected from 12 male patients diagnosed with bladder cancer, and from 11 healthy, age-matched individuals were analysed using 16S sequencing. Our results show that the most abundant phylum in both groups was Firmicutes, followed by Actinobacteria, Bacteroidetes and Proteobacteria. While microbial diversity and overall microbiome composition were not significantly different between groups, we could identify operational taxonomic units (OTUs) that were more abundant in either group. Among those that were significantly enriched in the bladder cancer group, we identified an OTU belonging to genus Fusobacterium, a possible protumorigenic pathogen. In an independent sample of 42 bladder cancer tissues, 11 had Fusobacterium nucleatum sequences detected by PCR. Three OTUs from genera Veillonella, Streptococcus and Corynebacterium were more abundant in healthy urines. However, due to the limited number of participants additional studies are needed to determine if urinary microbiome is associated with bladder cancer.

The kinetoplastid-infecting Bodo saltans virus (BsV), a window into the most abundant giant viruses in the sea.

Giant viruses are ecologically important players in aquatic ecosystems that have challenged concepts of what constitutes a virus. Herein, we present the giant Bodo saltans virus (BsV), the first characterized representative of the most abundant group of giant viruses in ocean metagenomes, and the first isolate of a klosneuvirus, a subgroup of the Mimiviridae proposed from metagenomic data. BsV infects an ecologically important microzooplankton, the kinetoplastid Bodo saltans. Its 1.39 Mb genome encodes 1227 predicted ORFs, including a complex replication machinery. Yet, much of its translational apparatus has been lost, including all tRNAs. Essential genes are invaded by homing endonuclease-encoding self-splicing introns that may defend against competing viruses. Putative anti-host factors show extensive gene duplication via a genomic accordion indicating an ongoing evolutionary arms race and highlighting the rapid evolution and genomic plasticity that has led to genome gigantism and the enigma that is giant viruses.

Dietary Gluten-Induced Gut Dysbiosis Is Accompanied by Selective Upregulation of microRNAs with Intestinal Tight Junction and Bacteria-Binding Motifs in Rhesus Macaque Model of Celiac Disease.

The composition of the gut microbiome reflects the overall health status of the host. In this study, stool samples representing the gut microbiomes from 6 gluten-sensitive (GS) captive juvenile rhesus macaques were compared with those from 6 healthy, age- and diet-matched peers. A total of 48 samples representing both groups were studied using V4 16S rRNA gene DNA analysis. Samples from GS macaques were further characterized based on type of diet administered: conventional monkey chow, i.e., wheat gluten-containing diet (GD), gluten-free diet (GFD), barley gluten-derived diet (BOMI) and reduced gluten barley-derived diet (RGB). It was hypothesized that the GD diet would lower the gut microbial diversity in GS macaques. This is the first report illustrating the reduction of gut microbial alpha-diversity (p < 0.05) following the consumption of dietary gluten in GS macaques. Selected bacterial families (e.g., Streptococcaceae and Lactobacillaceae) were enriched in GS macaques while Coriobacteriaceae was enriched in healthy animals. Within several weeks after the replacement of the GD by the GFD diet, the composition (beta-diversity) of gut microbiome in GS macaques started to change (p = 0.011) towards that of a normal macaque. Significance for alpha-diversity however, was not reached by the day 70 when the feeding experiment ended. Several inflammation-associated microRNAs (miR-203, -204, -23a, -23b and -29b) were upregulated (p < 0.05) in jejunum of 4 biopsied GS macaques fed GD with predicted binding sites on 16S ribosomal RNA of Lactobacillus reuteri (accession number: NR_025911), Prevotella stercorea (NR_041364) and Streptococcus luteciae (AJ297218) that were overrepresented in feces. Additionally, claudin-1, a validated tight junction protein target of miR-29b was significantly downregulated in jejunal epithelium of GS macaques. Taken together, we predict that with the introduction of effective treatments in future studies the diversity of gut microbiomes in GS macaques will approach those of healthy individuals. Further studies are needed to elucidate the regulatory pathways of inflammatory miRNAs in intestinal mucosa of GS macaques and to correlate their expression with gut dysbiosis.

Combining genomic sequencing methods to explore viral diversity and reveal potential virus-host interactions.

Viral diversity and virus-host interactions in oxygen-starved regions of the ocean, also known as oxygen minimum zones (OMZs), remain relatively unexplored. Microbial community metabolism in OMZs alters nutrient and energy flow through marine food webs, resulting in biological nitrogen loss and greenhouse gas production. Thus, viruses infecting OMZ microbes have the potential to modulate community metabolism with resulting feedback on ecosystem function. Here, we describe viral communities inhabiting oxic surface (10 m) and oxygen-starved basin (200 m) waters of Saanich Inlet, a seasonally anoxic fjord on the coast of Vancouver Island, British Columbia using viral metagenomics and complete viral fosmid sequencing on samples collected between April 2007 and April 2010. Of 6459 open reading frames (ORFs) predicted across all 34 viral fosmids, 77.6% (n = 5010) had no homology to reference viral genomes. These fosmids recruited a higher proportion of viral metagenomic sequences from Saanich Inlet than from nearby northeastern subarctic Pacific Ocean (Line P) waters, indicating differences in the viral communities between coastal and open ocean locations. While functional annotations of fosmid ORFs were limited, recruitment to NCBI's non-redundant "nr" database and publicly available single-cell genomes identified putative viruses infecting marine thaumarchaeal and SUP05 proteobacteria to provide potential host linkages with relevance to coupled biogeochemical cycling processes in OMZ waters. Taken together, these results highlight the power of coupled analyses of multiple sequence data types, such as viral metagenomic and fosmid sequence data with prokaryotic single cell genomes, to chart viral diversity, elucidate genomic and ecological contexts for previously unclassifiable viral sequences, and identify novel host interactions in natural and engineered ecosystems.

Seasonal and interannual variability of the marine bacterioplankton community throughout the water column over ten years.

Microbial activities that affect global oceanographic and atmospheric processes happen throughout the water column, yet the long-term ecological dynamics of microbes have been studied largely in the euphotic zone and adjacent seasonally mixed depths. We investigated temporal patterns in the community structure of free-living bacteria, by sampling approximately monthly from 5 m, the deep chlorophyll maximum (∼15-40 m), 150, 500 and 890 m, in San Pedro Channel (maximum depth 900 m, hypoxic below ∼500 m), off the coast of Southern California. Community structure and biodiversity (inverse Simpson index) showed seasonal patterns near the surface and bottom of the water column, but not at intermediate depths. Inverse Simpson's index was highest in the winter in surface waters and in the spring at 890 m, and varied interannually at all depths. Biodiversity appeared to be driven partially by exchange of microbes between depths and was highest when communities were changing slowly over time. Meanwhile, communities from the surface through 500 m varied interannually. After accounting for seasonality, several environmental parameters co-varied with community structure at the surface and 890 m, but not at the intermediate depths. Abundant and seasonally variable groups included, at 890 m, Nitrospina, Flavobacteria and Marine Group A. Seasonality at 890 m is likely driven by variability in sinking particles, which originate in surface waters, pass transiently through the middle water column and accumulate on the seafloor where they alter the chemical environment. Seasonal subeuphotic groups are likely those whose ecology is strongly influenced by these particles. This surface-to-bottom, decade-long, study identifies seasonality and interannual variability not only of overall community structure, but also of numerous taxonomic groups and near-species level operational taxonomic units.

Biogeography of Viruses in the Sea.

Viral ecology is a rapidly progressing area of research, as molecular methods have improved significantly for targeted research on specific populations and whole communities. To interpret and synthesize global viral diversity and distribution, it is feasible to assess whether macroecology concepts can apply to marine viruses. We review how viral and host life history and physical properties can influence viral distribution in light of biogeography and metacommunity ecology paradigms. We highlight analytical approaches that can be applied to emerging global data sets and meta-analyses to identify individual taxa with global influence and drivers of emergent properties that influence microbial community structure by drawing on examples across the spectrum of viral taxa, from RNA to ssDNA and dsDNA viruses.

Top-down controls on bacterial community structure: microbial network analysis of bacteria, T4-like viruses and protists.

Characterizing ecological relationships between viruses, bacteria and protists in the ocean are critical to understanding ecosystem function, yet these relationships are infrequently investigated together. We evaluated these relationships through microbial association network analysis of samples collected approximately monthly from March 2008 to January 2011 in the surface ocean (0-5 m) at the San Pedro Ocean Time series station. Bacterial, T4-like myoviral and protistan communities were described by Automated Ribosomal Intergenic Spacer Analysis and terminal restriction fragment length polymorphism of the gene encoding the major capsid protein (g23) and 18S ribosomal DNA, respectively. Concurrent shifts in community structure suggested similar timing of responses to environmental and biological parameters. We linked T4-like myoviral, bacterial and protistan operational taxonomic units by local similarity correlations, which were then visualized as association networks. Network links (correlations) potentially represent synergistic and antagonistic relationships such as viral lysis, grazing, competition or other interactions. We found that virus-bacteria relationships were more cross-linked than protist-bacteria relationships, suggestive of increased taxonomic specificity in virus-bacteria relationships. We also found that 80% of bacterial-protist and 74% of bacterial-viral correlations were positive, with the latter suggesting that at monthly and seasonal timescales, viruses may be following their hosts more often than controlling host abundance.

Temporal variability and coherence of euphotic zone bacterial communities over a decade in the Southern California Bight.

Time-series are critical to understanding long-term natural variability in the oceans. Bacterial communities in the euphotic zone were investigated for over a decade at the San Pedro Ocean Time-series station (SPOT) off southern California. Community composition was assessed by Automated Ribosomal Intergenic Spacer Analysis (ARISA) and coupled with measurements of oceanographic parameters for the surface ocean (0-5 m) and deep chlorophyll maximum (DCM, average depth ≈ 30 m). SAR11 and cyanobacterial ecotypes comprised typically more than one-third of the measured community; diversity within both was temporally variable, although a few operational taxonomic units (OTUs) were consistently more abundant. Persistent OTUs, mostly Alphaproteobacteria (SAR11 clade), Actinobacteria and Flavobacteria, tended to be abundant, in contrast to many rarer yet intermittent and ephemeral OTUs. Association networks revealed potential niches for key OTUs from SAR11, cyanobacteria, SAR86 and other common clades on the basis of robust correlations. Resilience was evident by the average communities drifting only slightly as years passed. Average Bray-Curtis similarity between any pair of dates was ≈ 40%, with a slight decrease over the decade and obvious near-surface seasonality; communities 8-10 years apart were slightly more different than those 1-4 years apart with the highest rate of change at 0-5 m between communities <4 years apart. The surface exhibited more pronounced seasonality than the DCM. Inter-depth Bray-Curtis similarities repeatedly decreased as the water column stratified each summer. Environmental factors were better predictors of shifts in community composition than months or elapsed time alone; yet, the best predictor was community composition at the other depth (that is, 0-5 m versus DCM).

Short-term observations of marine bacterial and viral communities: patterns, connections and resilience.

Observation of short-term temporal variation in bacterial and viral communities is important for understanding patterns of aquatic microbial diversity. We collected surface seawater once daily for 38 consecutive days with seven more samples interspersed over 40 more days at one location ∼2 km from Santa Catalina Island, California. Bacterial communities were analyzed by automated ribosomal intergenic spacer analysis (ARISA) and viral communities were analyzed by terminal restriction fragment length polymorphism (TRFLP) of the conserved T4-like myoviral gene encoding the major capsid protein (g23). Common bacterial and viral taxa were consistently dominant, and relatively few displayed dramatic increases/decreases or 'boom/bust' patterns that might be expected from dynamic predator-prey interactions. Association network analysis showed most significant covariations (associations) occurred among bacterial taxa or among viral taxa and there were several modular (highly-interconnected) associations (P≤0.005). Associations observed between bacteria and viruses (P≤0.005) occurred with a median time lag of 2 days. Regression of all pairwise Bray-Curtis similarities between samples indicated a rate of bacterial community change that slows from 2.1%-0.18% per day over a week to 2 months; the rate stays around 0.4% per day for viruses. Our interpretation is that, over the scale of days, individual bacterial and viral OTUs can be dynamic and patterned; resulting in statistical associations regarded as potential ecological interactions. However, over the scale of weeks, average bacterial community variation is slower, suggesting that there is strong community-level ecological resilience, that is, a tendency to converge towards a 'mean' microbial community set by longer-term controlling factors.

Seasonality and monthly dynamics of marine myovirus communities.

Marine myoviruses (i.e. bacteriophages with a contractile tail sheath) are numerically abundant and genetically diverse. We developed a terminal restriction fragment length polymorphism assay (TRFLP) for g23, the conserved gene encoding the major capsid protein, to investigate T4-like myovirus communities at USC's Microbial Observatory at the San Pedro Ocean Time-series (SPOT), where we previously reported bacterial seasonality. Between 71 and 154 operational taxonomic units (OTUs) were observed monthly over 3 years. Roughly 25% of OTUs were detected in 31 or more months. T4-like myoviral community structure varied seasonally with some OTUs peaking repeatedly in spring-summer and others in fall-winter, while moderately abundant OTUs persisted year-round. Recurring community structure was demonstrated using discriminant function analysis (DFA, selecting taxa that best predict months) and average Bray-Curtis similarity. DFA showed communities from adjacent months or 12 months apart were positively auto-correlated, while communities 3-7 months apart were negatively auto-correlated. Bray-Curtis similarity was highest between adjacent months - with a local maximum at 12-month and local minima at 6- and 18- to 20-month lags. The T4-like virus community at SPOT exhibited seasonality, yet the somewhat unexpected persistence of moderately abundant OTUs and predictability of the community add new twists to existing conceptual models of marine viruses.

Marine bacterial, archaeal and protistan association networks reveal ecological linkages.

Microbes have central roles in ocean food webs and global biogeochemical processes, yet specific ecological relationships among these taxa are largely unknown. This is in part due to the dilute, microscopic nature of the planktonic microbial community, which prevents direct observation of their interactions. Here, we use a holistic (that is, microbial system-wide) approach to investigate time-dependent variations among taxa from all three domains of life in a marine microbial community. We investigated the community composition of bacteria, archaea and protists through cultivation-independent methods, along with total bacterial and viral abundance, and physico-chemical observations. Samples and observations were collected monthly over 3 years at a well-described ocean time-series site of southern California. To find associations among these organisms, we calculated time-dependent rank correlations (that is, local similarity correlations) among relative abundances of bacteria, archaea, protists, total abundance of bacteria and viruses and physico-chemical parameters. We used a network generated from these statistical correlations to visualize and identify time-dependent associations among ecologically important taxa, for example, the SAR11 cluster, stramenopiles, alveolates, cyanobacteria and ammonia-oxidizing archaea. Negative correlations, perhaps suggesting competition or predation, were also common. The analysis revealed a progression of microbial communities through time, and also a group of unknown eukaryotes that were highly correlated with dinoflagellates, indicating possible symbioses or parasitism. Possible 'keystone' species were evident. The network has statistical features similar to previously described ecological networks, and in network parlance has non-random, small world properties (that is, highly interconnected nodes). This approach provides new insights into the natural history of microbes.

Separating the post-Glacial coancestry of European and Asian Y chromosomes within haplogroup R1a.

Human Y-chromosome haplogroup structure is largely circumscribed by continental boundaries. One notable exception to this general pattern is the young haplogroup R1a that exhibits post-Glacial coalescent times and relates the paternal ancestry of more than 10% of men in a wide geographic area extending from South Asia to Central East Europe and South Siberia. Its origin and dispersal patterns are poorly understood as no marker has yet been described that would distinguish European R1a chromosomes from Asian. Here we present frequency and haplotype diversity estimates for more than 2000 R1a chromosomes assessed for several newly discovered SNP markers that introduce the onset of informative R1a subdivisions by geography. Marker M434 has a low frequency and a late origin in West Asia bearing witness to recent gene flow over the Arabian Sea. Conversely, marker M458 has a significant frequency in Europe, exceeding 30% in its core area in Eastern Europe and comprising up to 70% of all M17 chromosomes present there. The diversity and frequency profiles of M458 suggest its origin during the early Holocene and a subsequent expansion likely related to a number of prehistoric cultural developments in the region. Its primary frequency and diversity distribution correlates well with some of the major Central and East European river basins where settled farming was established before its spread further eastward. Importantly, the virtual absence of M458 chromosomes outside Europe speaks against substantial patrilineal gene flow from East Europe to Asia, including to India, at least since the mid-Holocene.

Paleolithic Y-haplogroup heritage predominates in a Cretan highland plateau.

The island of Crete, credited by some historical scholars as a central crucible of western civilization, has been under continuous archeological investigation since the second half of the nineteenth century. In the present work, the geographic stratification of the contemporary Cretan Y-chromosome gene pool was assessed by high-resolution haplotyping to investigate the potential imprints of past colonization episodes and the population substructure. In addition to analyzing the possible geographic origins of Y-chromosome lineages in relatively accessible areas of the island, this study includes samples from the isolated interior of the Lasithi Plateau--a mountain plain located in eastern Crete. The potential significance of the results from the latter region is underscored by the possibility that this region was used as a Minoan refugium. Comparisons of Y-haplogroup frequencies among three Cretan populations as well as with published data from additional Mediterranean locations revealed significant differences in the frequency distributions of Y-chromosome haplogroups within the island. The most outstanding differences were observed in haplogroups J2 and R1, with the predominance of haplogroup R lineages in the Lasithi Plateau and of haplogroup J lineages in the more accessible regions of the island. Y-STR-based analyses demonstrated the close affinity that R1a1 chromosomes from the Lasithi Plateau shared with those from the Balkans, but not with those from lowland eastern Crete. In contrast, Cretan R1b microsatellite-defined haplotypes displayed more resemblance to those from Northeast Italy than to those from Turkey and the Balkans.

Polarity and temporality of high-resolution y-chromosome distributions in India identify both indigenous and exogenous expansions and reveal minor genetic influence of Central Asian pastoralists.

Although considerable cultural impact on social hierarchy and language in South Asia is attributable to the arrival of nomadic Central Asian pastoralists, genetic data (mitochondrial and Y chromosomal) have yielded dramatically conflicting inferences on the genetic origins of tribes and castes of South Asia. We sought to resolve this conflict, using high-resolution data on 69 informative Y-chromosome binary markers and 10 microsatellite markers from a large set of geographically, socially, and linguistically representative ethnic groups of South Asia. We found that the influence of Central Asia on the pre-existing gene pool was minor. The ages of accumulated microsatellite variation in the majority of Indian haplogroups exceed 10,000-15,000 years, which attests to the antiquity of regional differentiation. Therefore, our data do not support models that invoke a pronounced recent genetic input from Central Asia to explain the observed genetic variation in South Asia. R1a1 and R2 haplogroups indicate demographic complexity that is inconsistent with a recent single history. Associated microsatellite analyses of the high-frequency R1a1 haplogroup chromosomes indicate independent recent histories of the Indus Valley and the peninsular Indian region. Our data are also more consistent with a peninsular origin of Dravidian speakers than a source with proximity to the Indus and with significant genetic input resulting from demic diffusion associated with agriculture. Our results underscore the importance of marker ascertainment for distinguishing phylogenetic terminal branches from basal nodes when attributing ancestral composition and temporality to either indigenous or exogenous sources. Our reappraisal indicates that pre-Holocene and Holocene-era--not Indo-European--expansions have shaped the distinctive South Asian Y-chromosome landscape.