Comparisons of the Fungal and Protistan Communities among Different Marine Sponge Holobionts by Pyrosequencing

To date, the knowledge of eukaryotic communities associated with sponges remains limited compared with prokaryotic communities. In a manner similar to prokaryotes, it could be hypothesized that sponge holobionts have phylogenetically diverse eukaryotic symbionts, and the eukaryotic community structures in different sponge holobionts were probably different. In order to test this hypothesis, the communities of eukaryota associated with 11 species of South China Sea sponges were compared with the V4 region of 18S ribosomal ribonucleic acid gene using 454 pyrosequencing. Consequently, 135 and 721 unique operational taxonomic units (OTUs) of fungi and protists were obtained at 97 % sequence similarity, respectively. These sequences were assigned to 2 phyla of fungi (Ascomycota and Basidiomycota) and 9 phyla of protists including 5 algal phyla (Chlorophyta, Haptophyta, Streptophyta, Rhodophyta, and Stramenopiles) and 4 protozoal phyla (Alveolata, Cercozoa, Haplosporidia, and Radiolaria) including 47 orders (12 fungi, 35 protists). Entorrhizales of fungi and 18 orders of protists were detected in marine sponges for the first time. Particularly, Tilletiales of fungi and Chlorocystidales of protists were detected for the first time in marine habitats. Though Ascomycota, Alveolata, and Radiolaria were detected in all the 11 sponge species, sponge holobionts have different fungi and protistan communities according to OTU comparison and principal component analysis at the order level. This study provided the first insights into the fungal and protistan communities associated with different marine sponge holobionts using pyrosequencing, thus further extending the knowledge on sponge-associated eukaryotic diversity.     

Groundtruthing next-gen sequencing for microbial ecology-biases and errors in community structure estimates from PCR amplicon pyrosequencing

Analysis of microbial communities by high-throughput pyrosequencing of SSU rRNA gene PCR amplicons has transformed microbial ecology research and led to the observation that many communities contain a diverse assortment of rare taxa-a phenomenon termed the Rare Biosphere. Multiple studies have investigated the effect of pyrosequencing read quality on operational taxonomic unit (OTU) richness for contrived communities, yet there is limited information on the fidelity of community structure estimates obtained through this approach. Given that PCR biases are widely recognized, and further unknown biases may arise from the sequencing process itself, a priori assumptions about the neutrality of the data generation process are at best unvalidated. Furthermore, post-sequencing quality control algorithms have not been explicitly evaluated for the accuracy of recovered representative sequences and its impact on downstream analyses, reducing useful discussion on pyrosequencing reads to their diversity and abundances. Here we report on community structures and sequences recovered for in vitro-simulated communities consisting of twenty 16S rRNA gene clones tiered at known proportions. PCR amplicon libraries of the V3-V4 and V6 hypervariable regions from the in vitro-simulated communities were sequenced using the Roche 454 GS FLX Titanium platform. Commonly used quality control protocols resulted in the formation of OTUs with >1% abundance composed entirely of erroneous sequences, while over-aggressive clustering approaches obfuscated real, expected OTUs. The pyrosequencing process itself did not appear to impose significant biases on overall community structure estimates, although the detection limit for rare taxa may be affected by PCR amplicon size and quality control approach employed. Meanwhile, PCR biases associated with the initial amplicon generation may impose greater distortions in the observed community structure.     

Changes in community structure of active protistan assemblages from the lower Pearl River to coastal Waters of the South China Sea

Protists make up an important component of aquatic ecosystems, playing crucial roles in biogeochemical processes on local and global scales. To reveal the changes of diversity and community structure of protists along the salinity gradients, community compositions of active protistan assemblages were characterized along a transect from the lower Pearl River estuary to the open waters of the South China Sea (SCS), using high-throughput sequencing of the hyper-variable V9 regions of 18S rRNA. This study showed that the alpha diversity of protists, both in the freshwater and in the coastal SCS stations was higher than that in the estuary. The protist community structure also changed along the salinity gradient. The relative sequence abundance of Stramenopiles was highest at stations with lower salinity and decreased with the increasing of salinity. By contrast, the contributions of Alveolata, Hacrobia and Rhizaria to the protistan communities generally increased with the increasing of salinity. The composition of the active protistan community was strongly correlated with salinity, indicating that salinity was the dominant factor among measured environmental parameters affecting protistan community composition and structure.     

Analysis,Optimization and Verification of Illumina-Generated 16S rRNA Gene Amplicon Surveys

The exploration of microbial communities by sequencing 16S rRNA genes has expanded with low-cost, high-throughput sequencing instruments. Illumina-based 16S rRNA gene sequencing has recently gained popularity over 454 pyrosequencing due to its lower costs, higher accuracy and greater throughput. Although recent reports suggest that Illumina and 454 pyrosequencing provide similar beta diversity measures, it remains to be demonstrated that pre-existing 454 pyrosequencing workflows can transfer directly from 454 to Illumina MiSeq sequencing by simply changing the sequencing adapters of the primers. In this study, we modified 454 pyrosequencing primers targeting the V4-V5 hyper-variable regions of the 16S rRNA gene to be compatible with Illumina sequencers. Microbial communities from cows, humans, leeches, mice, sewage, and termites and a mock community were analyzed by 454 and MiSeq sequencing of the V4-V5 region and MiSeq sequencing of the V4 region. Our analysis revealed that reference-based OTU clustering alone introduced biases compared to de novo clustering, preventing certain taxa from being observed in some samples. Based on this we devised and recommend an analysis pipeline that includes read merging, contaminant filtering, and reference-based clustering followed by de novo OTU clustering, which produces diversity measures consistent with de novo OTU clustering analysis. Low levels of dataset contamination with Illumina sequencing were discovered that could affect analyses that require highly sensitive approaches. While moving to Illumina-based sequencing platforms promises to provide deeper insights into the breadth and function of microbial diversity, our results show that care must be taken to ensure that sequencing and processing artifacts do not obscure true microbial diversity.     

Investigating Microbial Eukaryotic Diversity from a Global Census: Insights from a Comparison of Pyrotag and Full-Length Sequences of 18S rRNA Genes

Next-generation DNA sequencing (NGS) approaches are rapidly surpassing Sanger sequencing for characterizing the diversity of natural microbial communities. Despite this rapid transition, few comparisons exist between Sanger sequences and the generally much shorter reads of NGS. Operational taxonomic units (OTUs) derived from full-length (Sanger sequencing) and pyrotag (454 sequencing of the V9 hypervariable region) sequences of 18S rRNA genes from 10 global samples were analyzed in order to compare the resulting protistan community structures and species richness. Pyrotag OTUs called at 98% sequence similarity yielded numbers of OTUs that were similar overall to those for full-length sequences when the latter were called at 97% similarity. Singleton OTUs strongly influenced estimates of species richness but not the higher-level taxonomic composition of the community. The pyrotag and full-length sequence data sets had slightly different taxonomic compositions of rhizarians, stramenopiles, cryptophytes, and haptophytes, but the two data sets had similarly high compositions of alveolates. Pyrotag-based OTUs were often derived from sequences that mapped to multiple full-length OTUs at 100% similarity. Thus, pyrotags sequenced from a single hypervariable region might not be appropriate for establishing protistan species-level OTUs. However, nonmetric multidimensional scaling plots constructed with the two data sets yielded similar clusters, indicating that beta diversity analysis results were similar for the Sanger and NGS sequences. Short pyrotag sequences can provide holistic assessments of protistan communities, although care must be taken in interpreting the results. The longer reads (>500 bp) that are now becoming available through NGS should provide powerful tools for assessing the diversity of microbial eukaryotic assemblages.     

Target region selection is a critical determinant of community fingerprints generated by 16S pyrosequencing

Pyrosequencing of 16S rRNA genes allows for in-depth characterization of complex microbial communities. Although it is known that primer selection can influence the profile of a community generated by sequencing, the extent and severity of this bias on deep-sequencing methodologies is not well elucidated. We tested the hypothesis that the hypervariable region targeted for sequencing and primer degeneracy play important roles in influencing the composition of 16S pyrotag communities. Subgingival plaque from deep sites of current smokers with chronic periodontitis was analyzed using Sanger sequencing and pyrosequencing using 4 primer pairs. Greater numbers of species were detected by pyrosequencing than by Sanger sequencing. Rare taxa constituted nearly 6% of each pyrotag community and less than 1% of the Sanger sequencing community. However, the different target regions selected for pyrosequencing did not demonstrate a significant difference in the number of rare and abundant taxa detected. The genera Prevotella, Fusobacterium, Streptococcus, Granulicatella, Bacteroides, Porphyromonas and Treponema were abundant when the V1-V3 region was targeted, while Streptococcus, Treponema, Prevotella, Eubacterium, Porphyromonas, Campylobacter and Enterococcus predominated in the community generated by V4-V6 primers, and the most numerous genera in the V7-V9 community were Veillonella, Streptococcus, Eubacterium, Enterococcus, Treponema, Catonella and Selenomonas. Targeting the V4-V6 region failed to detect the genus Fusobacterium, while the taxa Selenomonas, TM7 and Mycoplasma were not detected by the V7-V9 primer pairs. The communities generated by degenerate and non-degenerate primers did not demonstrate significant differences. Averaging the community fingerprints generated by V1-V3 and V7-V9 primers provided results similar to Sanger sequencing, while allowing a significantly greater depth of coverage than is possible with Sanger sequencing. It is therefore important to use primers targeted to these two regions of the 16S rRNA gene in all deep-sequencing efforts to obtain representational characterization of complex microbial communities.     

Short pyrosequencing reads suffice for accurate microbial community analysis

Pyrosequencing technology allows us to characterize microbial communities using 16S ribosomal RNA (rRNA) sequences orders of magnitude faster and more cheaply than has previously been possible. However, results from different studies using pyrosequencing and traditional sequencing are often difficult to compare, because amplicons covering different regions of the rRNA might yield different conclusions. We used sequences from over 200 globally dispersed environments to test whether studies that used similar primers clustered together mistakenly, without regard to environment. We then tested whether primer choice affects sequence-based community analyses using UniFrac, our recently-developed method for comparing microbial communities. We performed three tests of primer effects. We tested whether different simulated amplicons generated the same UniFrac clustering results as near-full-length sequences for three recent large-scale studies of microbial communities in the mouse and human gut, and the Guerrero Negro microbial mat. We then repeated this analysis for short sequences (100-, 150-, 200- and 250-base reads) resembling those produced by pyrosequencing. The results show that sequencing effort is best focused on gathering more short sequences rather than fewer longer ones, provided that the primers are chosen wisely, and that community comparison methods such as UniFrac are surprisingly robust to variation in the region sequenced.     

Analysis of the Cystic Fibrosis Lung Microbiota via Serial Illumina Sequencing of Bacterial 16S rRNA Hypervariable Regions

The characterization of bacterial communities using DNA sequencing has revolutionized our ability to study microbes in nature and discover the ways in which microbial communities affect ecosystem functioning and human health. Here we describe Serial Illumina Sequencing (SI-Seq): a method for deep sequencing of the bacterial 16S rRNA gene using next-generation sequencing technology. SI-Seq serially sequences portions of the V5, V6 and V7 hypervariable regions from barcoded 16S rRNA amplicons using an Illumina short-read genome analyzer. SI-Seq obtains taxonomic resolution similar to 454 pyrosequencing for a fraction of the cost, and can produce hundreds of thousands of reads per sample even with very high multiplexing. We validated SI-Seq using single species and mock community controls, and via a comparison to cystic fibrosis lung microbiota sequenced using 454 FLX Titanium. Our control runs show that SI-Seq has a dynamic range of at least five orders of magnitude, can classify >96% of sequences to the genus level, and performs just as well as 454 and paired-end Illumina methods in estimation of standard microbial ecology diversity measurements. We illustrate the utility of SI-Seq in a pilot sample of central airway secretion samples from cystic fibrosis patients.     

Illumina-based analysis of microbial community diversity

Microbes commonly exist in milieus of varying complexity and diversity. Although cultivation-based techniques have been unable to accurately capture the true diversity within microbial communities, these deficiencies have been overcome by applying molecular approaches that target the universally conserved 16S ribosomal RNA gene. The recent application of 454 pyrosequencing to simultaneously sequence thousands of 16S rDNA sequences (pyrotags) has revolutionized the characterization of complex microbial communities. To date, studies based on 454 pyrotags have dominated the field, but sequencing platforms that generate many more sequence reads at much lower costs have been developed. Here, we use the Illumina sequencing platform to design a strategy for 16S amplicon analysis (iTags), and assess its generality, practicality and potential complications. We fabricated and sequenced paired-end libraries of amplified hyper-variable 16S rDNA fragments from sets of samples that varied in their contents, ranging from a single bacterium to highly complex communities. We adopted an approach that allowed us to evaluate several potential sources of errors, including sequencing artifacts, amplification biases, non-corresponding paired-end reads and mistakes in taxonomic classification. By considering each source of error, we delineate ways to make biologically relevant and robust conclusions from the millions of sequencing reads that can be readily generated by this technology.     

Global biogeography of highly diverse protistan communities in soil

Protists are ubiquitous members of soil microbial communities, but the structure of these communities, and the factors that influence their diversity, are poorly understood. We used barcoded pyrosequencing to survey comprehensively the diversity of soil protists from 40 sites across a broad geographic range that represent a variety of biome types, from tropical forests to deserts. In addition to taxa known to be dominant in soil, including Cercozoa and Ciliophora, we found high relative abundances of groups such as Apicomplexa and Dinophyceae that have not previously been recognized as being important components of soil microbial communities. Soil protistan communities were highly diverse, approaching the extreme diversity of their bacterial counterparts across the same sites. Like bacterial taxa, protistan taxa were not globally distributed, and the composition of these communities diverged considerably across large geographic distances. However, soil protistan and bacterial communities exhibit very different global-scale biogeographical patterns, with protistan communities strongly structured by climatic conditions that regulate annual soil moisture availability.     

Next-generation phylogeography: a targeted approach for multilocus sequencing of non-model organisms

The field of phylogeography has long since realized the need and utility of incorporating nuclear DNA (nDNA) sequences into analyses. However, the use of nDNA sequence data, at the population level, has been hindered by technical laboratory difficulty, sequencing costs, and problematic analytical methods dealing with genotypic sequence data, especially in non-model organisms. Here, we present a method utilizing the 454 GS-FLX Titanium pyrosequencing platform with the capacity to simultaneously sequence two species of sea star (Meridiastra calcar and Parvulastra exigua) at five different nDNA loci across 16 different populations of 20 individuals each per species. We compare results from 3 populations with traditional Sanger sequencing based methods, and demonstrate that this next-generation sequencing platform is more time and cost effective and more sensitive to rare variants than Sanger based sequencing. A crucial advantage is that the high coverage of clonally amplified sequences simplifies haplotype determination, even in highly polymorphic species. This targeted next-generation approach can greatly increase the use of nDNA sequence loci in phylogeographic and population genetic studies by mitigating many of the time, cost, and analytical issues associated with highly polymorphic, diploid sequence markers.     

Channeling of bacterioplanktonic production toward phagotrophic flagellates and ciliates under different seasonal conditions in a river

The objective of this study was to analyze the flux of biomass through the communities of bacteria and phagotrophic protists in the cold and warm conditions occurring seasonally in Butrón River. Bacterial and heterotrophic protistan (flagellate and ciliate) abundance was determined by epifluorescence direct counts; protistan grazing on planktonic bacteria was measured from fluorescently labeled bacteria uptake rates; and the estimate of bacterial secondary production was obtained from [³H]thymidine incorporation rates. The abundance of bacterial, flagellate, and ciliate communities was similar during cold and warm situations. However, we observed that estimates of dynamic parameters, i.e., secondary bacterial production and protistan grazing, in both situations were noticeably different. In the warm situation, grazing rates of flagellates and ciliates (bacteria per protist per hour) were, respectively, 7 times and 18 times higher than those determined in the cold situation, and the grazing rates of the protistan communities (bacteria per protists present in 1 ml of water per hour) increased up to 5 times in the case of flagellates and 42 times in the case of ciliates. Estimates of bacterial secondary production were also higher during the warm situation, showing a ninefold increase. The percentage of bacterial production preyed upon by flagellates or ciliates was not significantly different between the two conditions. These results showed that in the different conditions of a system, the flux of biomass between the trophic levels may be quite different although this process may not be reflected in the abundance of each community of bacteria, flagellates, and ciliates. Offprint requests to: J. Iriberri.     

Cyanobacteria drive community composition and functionality in rock–soil interface communities

Most ecological research on hypoliths, significant primary producers in hyperarid deserts, has focused on the diversity of individual groups of microbes (i.e. bacteria). However, microbial communities are inherently complex, and the interactions between cyanobacteria, heterotrophic bacteria, protista and metazoa are likely to be very important for ecosystem functioning. Cyanobacterial and heterotrophic bacterial communities were analysed by pyrosequencing, while metazoan and protistan communities were assessed by T‐RFLP analysis. Microbial functionality was estimated using carbon substrate utilization. Cyanobacterial community composition was significant in shaping community structure and function in hypoliths. Ecological network analysis showed that most significant co‐occurrences were positive, representing potential synergistic interactions. There were several highly interconnected associations (modules), and specific cyanobacteria were important in driving the modular structure of hypolithic networks. Together, our results suggest that hypolithic cyanobacteria have strong effects on higher trophic levels and ecosystem functioning.     

A diverse bacterial community in an anoxic quinoline-degrading bioreactor determined by using pyrosequencing and clone library analysis

There is a concern of whether the structure and diversity of a microbial community can be effectively revealed by short-length pyrosequencing reads. In this study, we performed a microbial community analysis on a sample from a high-efficiency denitrifying quinoline-degrading bioreactor and compared the results generated by pyrosequencing with those generated by clone library technology. By both technologies, 16S rRNA gene analysis indicated that the bacteria in the sample were closely related to, for example, Proteobacteria, Actinobacteria, and Bacteroidetes. The sequences belonging to Rhodococcus were the most predominant, and Pseudomonas, Sphingomonas, Acidovorax, and Zoogloea were also abundant. Both methods revealed a similar overall bacterial community structure. However, the 622 pyrosequencing reads of the hypervariable V3 region of the 16S rRNA gene revealed much higher bacterial diversity than the 130 sequences from the full-length 16S rRNA gene clone library. The 92 operational taxonomic unit (OTUs) detected using pyrosequencing belonged to 45 families, whereas the 37 OTUs found in the clone library belonged to 25 families. Most sequences obtained from the clone library had equivalents in the pyrosequencing reads. However, 64 OTUs detected by pyrosequencing were not represented in the clone library. Our results demonstrate that pyrosequencing of the V3 region of the 16S rRNA gene is not only a powerful tool for discovering low-abundance bacterial populations but is also reliable for dissecting the bacterial community structure in a wastewater environment.     

Ultra-Deep Pyrosequencing (UDPS) Data Treatment to Study Amplicon HCV Minor Variants

We have investigated the reliability and reproducibility of HCV viral quasispecies quantification by ultra-deep pyrosequencing (UDPS) methods. Our study has been divided in two parts. First of all, by UDPS sequencing of clone mixes samples we have established the global noise level of UDPS and fine tuned a data treatment workflow previously optimized for HBV sequence analysis. Secondly, we have studied the reproducibility of the methodology by comparing 5 amplicons from two patient samples on three massive sequencing platforms (FLX+, FLX and Junior) after applying the error filters developed from the clonal/control study. After noise filtering the UDPS results, the three replicates showed the same 12 polymorphic sites above 0.7%, with a mean CV of 4.86%. Two polymorphic sites below 0.6% were identified by two replicates and one replicate respectively. A total of 25, 23 and 26 haplotypes were detected by GS-Junior, GS-FLX and GS-FLX+. The observed CVs for the normalized Shannon entropy (Sn), the mutation frequency (Mf), and the nucleotidic diversity (Pi) were 1.46%, 3.96% and 3.78%. The mean absolute difference in the two patients (5 amplicons each), in the GS-FLX and GS-FLX+, were 1.46%, 3.96% and 3.78% for Sn, Mf and Pi. No false polymorphic site was observed above 0.5%.Our results indicate that UDPS is an optimal alternative to molecular cloning for quantitative study of HCV viral quasispecies populations, both in complexity and composition. We propose an UDPS data treatment workflow for amplicons from the RNA viral quasispecies which, at a sequencing depth of at least 10,000 reads per strand, enables to obtain sequences and frequencies of consensus haplotypes above 0.5% abundance with no erroneous mutations, with high confidence, resistant mutants as minor variants at the level of 1%, with high confidence that variants are not missed, and highly confident measures of quasispecies complexity.     

Development of Single-Nucleotide Polymorphism- Based Phylum-Specific PCR Amplification Technique: Application to the Community Analysis Using Ciliates as a Reference Organism

Despite recent advance in mass sequencing technologies such as pyrosequencing, assessment of culture-independent microbial eukaryote community structures using universal primers remains very difficult due to the tremendous richness and complexity of organisms in these communities. Use of a specific PCR marker targeting a particular group would provide enhanced sensitivity and more in-depth evaluation of microbial eukaryote communities compared to what can be achieved with universal primers. We discovered that many phylum- or group-specific single-nucleotide polymorphisms (SNPs) exist in small subunit ribosomal RNA (SSU rRNA) genes from diverse eukaryote groups. By applying this discovery to a known simple allele-discriminating (SAP) PCR method, we developed a technique that enables the identification of organisms belonging to a specific higher taxonomic group (or phylum) among diverse types of eukaryotes. We performed an assay using two complementary methods, pyrosequencing and clone library screening. In doing this, specificities for the group (ciliates) targeted in this study in bulked environmental samples were 94.6% for the clone library and 99.2% for pyrosequencing, respectively. In particular, our novel technique showed high selectivity for rare species, a feature that may be more important than the ability to identify quantitatively predominant species in community structure analyses. Additionally, our data revealed that a target-specific library (or ciliate-specific one for the present study) can better explain the ecological features of a sampling locality than a universal library.     

Exploring Microbial Diversity and Taxonomy Using SSU rRNA Hypervariable Tag Sequencing

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.     

Evaluation of the Ion Torrent Personal Genome Machine for Gene-Targeted Studies Using Amplicons of the Nitrogenase Gene nifH

The sequencing chips and kits of the Ion Torrent Personal Genome Machine (PGM), which employs semiconductor technology to measure pH changes in polymerization events, have recently been upgraded. The quality of PGM sequences has not been reassessed, and results have not been compared in the context of a gene-targeted microbial ecology study. To address this, we compared sequence profiles across available PGM chips and chemistries and with 454 pyrosequencing data by determining error types and rates and diazotrophic community structures. The PGM was then used to assess differences in nifH-harboring bacterial community structure among four corn-based cropping systems. Using our suggested filters from mock community analyses, the overall error rates were 0.62, 0.36, and 0.39% per base for chips 318 and 314 with the 400-bp kit and chip 318 with the Hi-Q chemistry, respectively. Compared with the 400-bp kit, the Hi-Q kit reduced indel rates by 28 to 59% and produced one to seven times more reads acceptable for downstream analyses. The PGM produced higher frameshift rates than pyrosequencing that were corrected by the RDP FrameBot tool. Significant differences among platforms were identified, although the diversity indices and overall site-based conclusions remained similar. For the cropping system analyses, a total of 6,182 unique NifH operational taxonomic units at 5% amino acid dissimilarity were obtained. The current crop type, as well as the crop rotation history, significantly influenced the composition of the soil diazotrophic community detected.     

Bacterial Communities Associated with the Leaves and the Roots of Arabidopsis thaliana

Diverse communities of bacteria inhabit plant leaves and roots and those bacteria play a crucial role for plant health and growth. Arabidopsis thaliana is an important model to study plant pathogen interactions, but little is known about its associated bacterial community under natural conditions. We used 454 pyrosequencing to characterize the bacterial communities associated with the roots and the leaves of wild A. thaliana collected at 4 sites; we further compared communities on the outside of the plants with communities in the endophytic compartments. We found that the most heavily sequenced bacteria in A. thaliana associated community are related to culturable species. Proteobacteria, Actinobacteria, and Bacteroidetes are the most abundant phyla in both leaf and root samples. At the genus level, sequences of Massilia and Flavobacterium are prevalent in both samples. Organ (leaf vs root) and habitat (epiphytes vs endophytes) structure the community. In the roots, richness is higher in the epiphytic communities compared to the endophytic compartment (P = 0.024), while the reverse is true for the leaves (P = 0.032). Interestingly, leaf and root endophytic compartments do not differ in richness, diversity and evenness, while they differ in community composition (P = 0.001). The results show that although the communities associated with leaves and roots share many bacterial species, the associated communities differ in structure.