Genetic risk for autoimmunity is associated with distinct changes in the human gut microbiome

Susceptibility to many human autoimmune diseases is under strong genetic control by class II human leukocyte antigen (HLA) allele combinations. These genes remain by far the greatest risk factors in the development of type 1 diabetes and celiac disease. Despite this, little is known about HLA influences on the composition of the human gut microbiome, a potential source of environmental influence on disease. Here, using a general population cohort from the All Babies in Southeast Sweden study, we report that genetic risk for developing type 1 diabetes autoimmunity is associated with distinct changes in the gut microbiome. Both the core microbiome and beta diversity differ with HLA risk group and genotype. In addition, protective HLA haplotypes are associated with bacterial genera Intestinibacter and Romboutsia. Thus, general population cohorts are valuable in identifying potential environmental triggers or protective factors for autoimmune diseases that may otherwise be masked by strong genetic control

Robust estimation of microbial diversity in theory and in practice

Quantifying diversity is of central importance for the study of structure, function and evolution of microbial communities. The estimation of microbial diversity has received renewed attention with the advent of large-scale metagenomic studies. Here, we consider what the diversity observed in a sample tells us about the diversity of the community being sampled. First, we argue that one cannot reliably estimate the absolute and relative number of microbial species present in a community without making unsupported assumptions about species abundance distributions. The reason for this is that sample data do not contain information about the number of rare species in the tail of species abundance distributions. We illustrate the difficulty in comparing species richness estimates by applying Chao's estimator of species richness to a set of in silico communities: they are ranked incorrectly in the presence of large numbers of rare species. Next, we extend our analysis to a general family of diversity metrics ("Hill diversities"), and construct lower and upper estimates of diversity values consistent with the sample data. The theory generalizes Chao's estimator, which we retrieve as the lower estimate of species richness. We show that Shannon and Simpson diversity can be robustly estimated for the in silico communities. We analyze nine metagenomic data sets from a wide range of environments, and show that our findings are relevant for empirically-sampled communities. Hence, we recommend the use of Shannon and Simpson diversity rather than species richness in efforts to quantify and compare microbial diversity.Comment: To be published in The ISME Journal. Main text: 16 pages, 5 figures. Supplement: 16 pages, 4 figure

Population dynamics and identification of efficient strains of Azospirillum in maize ecosystems of Bihar (India)

Information on inoculum load and diversity of native microbial community is an important prerequisite for crop management of microbial origin. Azospirillum has a proven role in benefiting the maize (Zea mays) crop in terms of nutrient (nitrogen) supply as well as plant growth enhancement. Bihar state has highest average national maize productivity although fertilizer consumption is minimum, indicating richness of Azospirillum both in terms of population and diversity in soils. An experiment was planned to generate basic information on Azospirillum population variation in maize soils under different agricultural practices and soil types of Bihar, to identify suitable agricultural practices supporting the target microorganism and efficient Azospirillum strain(s). No tillage, growing traditional maize cultivar, land use history (diara soil having history of maize cultivation), soil organic carbon (>1%) and intercrop with oat supported prevalence of Azospirillum in maize rhizosphere. Native Azospirillum population varied from 1 million to 1 billion/g soil under diverse agricultural practices and soil types. Such richness, however, does not necessarily mean that artificial inoculation of Azospirillum is not required in Bihar soils as 92% of Azospirillum isolates (50 isolates) were poor in nitrogen-fixing ability and 88% were poor on IAA production. Efficient strains of Azospirillum based on growth (three), acetylene reduction assay (three), IAA production (three), broad range of pH (two) and temperature tolerance were identified. The findings suggested that maize crop in Bihar should be inoculated in universal mode rather than site-specific mode

Conservative Fragments in Bacterial 16S rRNA Genes and Primer Design for 16S Ribosomal DNA Amplicons in Metagenomic Studies

Bacterial 16S ribosomal DNA (rDNA) amplicons have been widely used in the classification of uncultured bacteria inhabiting environmental niches. Primers targeting conservative regions of the rDNAs are used to generate amplicons of variant regions that are informative in taxonomic assignment. One problem is that the percentage coverage and application scope of the primers used in previous studies are largely unknown. In this study, conservative fragments of available rDNA sequences were first mined and then used to search for candidate primers within the fragments by measuring the coverage rate defined as the percentage of bacterial sequences containing the target. Thirty predicted primers with a high coverage rate (>90%) were identified, which were basically located in the same conservative regions as known primers in previous reports, whereas 30% of the known primers were associated with a coverage rate of <90%. The application scope of the primers was also examined by calculating the percentages of failed detections in bacterial phyla. Primers A519–539, E969–983, E1063–1081, U515 and E517, are highly recommended because of their high coverage in almost all phyla. As expected, the three predominant phyla, Firmicutes, Gemmatimonadetes and Proteobacteria, are best covered by the predicted primers. The primers recommended in this report shall facilitate a comprehensive and reliable survey of bacterial diversity in metagenomic studies

The Effects of Alignment Quality, Distance Calculation Method, Sequence Filtering, and Region on the Analysis of 16S rRNA Gene-Based Studies

Pyrosequencing of PCR-amplified fragments that target variable regions within the 16S rRNA gene has quickly become a powerful method for analyzing the membership and structure of microbial communities. This approach has revealed and introduced questions that were not fully appreciated by those carrying out traditional Sanger sequencing-based methods. These include the effects of alignment quality, the best method of calculating pairwise genetic distances for 16S rRNA genes, whether it is appropriate to filter variable regions, and how the choice of variable region relates to the genetic diversity observed in full-length sequences. I used a diverse collection of 13,501 high-quality full-length sequences to assess each of these questions. First, alignment quality had a significant impact on distance values and downstream analyses. Specifically, the greengenes alignment, which does a poor job of aligning variable regions, predicted higher genetic diversity, richness, and phylogenetic diversity than the SILVA and RDP-based alignments. Second, the effect of different gap treatments in determining pairwise genetic distances was strongly affected by the variation in sequence length for a region; however, the effect of different calculation methods was subtle when determining the sample's richness or phylogenetic diversity for a region. Third, applying a sequence mask to remove variable positions had a profound impact on genetic distances by muting the observed richness and phylogenetic diversity. Finally, the genetic distances calculated for each of the variable regions did a poor job of correlating with the full-length gene. Thus, while it is tempting to apply traditional cutoff levels derived for full-length sequences to these shorter sequences, it is not advisable. Analysis of β-diversity metrics showed that each of these factors can have a significant impact on the comparison of community membership and structure. Taken together, these results urge caution in the design and interpretation of analyses using pyrosequencing data

Effect of Cry1Ab Protein on Rhizobacterial Communities of Bt-Maize over a Four-Year Cultivation Period

Background: Bt-maize is a transgenic variety of maize expressing the Cry toxin from Bacillus turingiensis. The potential accumulation of the relative effect of the transgenic modification and the cry toxin on the rhizobacterial communities of Btmaize has been monitored over a period of four years. Methodology/Principal Findings: The accumulative effects of the cultivation of this transgenic plant have been monitored by means of high throughput DNA pyrosequencing of the bacterial DNA coding for the 16S rRNA hypervariable V6 region from rhizobacterial communities. The obtained sequences were subjected to taxonomic, phylogenetic and taxonomicindependent diversity studies. The results obtained were consistent, indicating that variations detected in the rhizobacterial community structure were possibly due to climatic factors rather than to the presence of the Bt-gene. No variations were observed in the diversity estimates between non-Bt and Bt-maize. Conclusions/Significance: The cultivation of Bt-maize during the four-year period did not change the maize rhizobacterial communities when compared to those of the non-Bt maize. This is the first study to be conducted with Bt-maize during such a long cultivation period and the first evaluation of rhizobacterial communities to be performed in this transgenic plant using Next Generation Sequencing

Amplification by PCR Artificially Reduces the Proportion of the Rare Biosphere in Microbial Communities

The microbial world has been shown to hold an unimaginable diversity. The use of rRNA genes and PCR amplification to assess microbial community structure and diversity present biases that need to be analyzed in order to understand the risks involved in those estimates. Herein, we show that PCR amplification of specific sequence targets within a community depends on the fractions that those sequences represent to the total DNA template. Using quantitative, real-time, multiplex PCR and specific Taqman probes, the amplification of 16S rRNA genes from four bacterial species within a laboratory community were monitored. Results indicate that the relative amplification efficiency for each bacterial species is a nonlinear function of the fraction that each of those taxa represent within a community or multispecies DNA template. Consequently, the low-proportion taxa in a community are under-represented during PCR-based surveys and a large number of sequences might need to be processed to detect some of the bacterial taxa within the ‘rare biosphere’. The structure of microbial communities from PCR-based surveys is clearly biased against low abundant taxa which are required to decipher the complete extent of microbial diversity in nature

New Detection Systems of Bacteria Using Highly Selective Media Designed by SMART: Selective Medium-Design Algorithm Restricted by Two Constraints

Culturing is an indispensable technique in microbiological research, and culturing with selective media has played a crucial role in the detection of pathogenic microorganisms and the isolation of commercially useful microorganisms from environmental samples. Although numerous selective media have been developed in empirical studies, unintended microorganisms often grow on such media probably due to the enormous numbers of microorganisms in the environment. Here, we present a novel strategy for designing highly selective media based on two selective agents, a carbon source and antimicrobials. We named our strategy SMART for highly Selective Medium-design Algorithm Restricted by Two constraints. To test whether the SMART method is applicable to a wide range of microorganisms, we developed selective media for Burkholderia glumae, Acidovorax avenae, Pectobacterium carotovorum, Ralstonia solanacearum, and Xanthomonas campestris. The series of media developed by SMART specifically allowed growth of the targeted bacteria. Because these selective media exhibited high specificity for growth of the target bacteria compared to established selective media, we applied three notable detection technologies: paper-based, flow cytometry-based, and color change-based detection systems for target bacteria species. SMART facilitates not only the development of novel techniques for detecting specific bacteria, but also our understanding of the ecology and epidemiology of the targeted bacteria

Fine-Scale Bacterial Beta Diversity within a Complex Ecosystem (Zodletone Spring, OK, USA): The Role of the Rare Biosphere

The adaptation of pyrosequencing technologies for use in culture-independent diversity surveys allowed for deeper sampling of ecosystems of interest. One extremely well suited area of interest for pyrosequencing-based diversity surveys that has received surprisingly little attention so far, is examining fine scale (e.g. micrometer to millimeter) beta diversity in complex microbial ecosystems.We examined the patterns of fine scale Beta diversity in four adjacent sediment samples (1mm apart) from the source of an anaerobic sulfide and sulfur rich spring (Zodletone spring) in southwestern Oklahoma, USA. Using pyrosequencing, a total of 292,130 16S rRNA gene sequences were obtained. The beta diversity patterns within the four datasets were examined using various qualitative and quantitative similarity indices. Low levels of Beta diversity (high similarity indices) were observed between the four samples at the phylum-level. However, at a putative species (OTU(0.03)) level, higher levels of beta diversity (lower similarity indices) were observed. Further examination of beta diversity patterns within dominant and rare members of the community indicated that at the putative species level, beta diversity is much higher within rare members of the community. Finally, sub-classification of rare members of Zodletone spring community based on patterns of novelty and uniqueness, and further examination of fine scale beta diversity of each of these subgroups indicated that members of the community that are unique, but non novel showed the highest beta diversity within these subgroups of the rare biosphere.The results demonstrate the occurrence of high inter-sample diversity within seemingly identical samples from a complex habitat. We reason that such unexpected diversity should be taken into consideration when exploring gamma diversity of various ecosystems, as well as planning for sequencing-intensive metagenomic surveys of highly complex ecosystems