Emilio O. Casamayor

Using network analysis to explore co-occurrence patterns in soil microbial communities

Exploring large environmental datasets generated by high-throughput DNA sequencing technologies requires new analytical approaches to move beyond the basic inventory descriptions of the composition and diversity of natural microbial communities. In order to investigate potential interactions between microbial taxa, network analysis of significant taxon co-occurrence patterns may help to decipher the structure of complex microbial communities across spatial or temporal gradients. Here, we calculated associations between microbial taxa and applied network analysis approaches to a 16S rRNA gene barcoded pyrosequencing dataset containing >160 000 bacterial and archaeal sequences from 151 soil samples from a broad range of ecosystem types. We described the topology of the resulting network and defined operational taxonomic unit categories based on abundance and occupancy (that is, habitat generalists and habitat specialists). Co-occurrence patterns were readily revealed, including general non-random association, common life history strategies at broad taxonomic levels and unexpected relationships between community members. Overall, we demonstrated the potential of exploring inter-taxa correlations to gain a more integrated understanding of microbial community structure and the ecological rules guiding community assembly.

Identification of and Spatio-Temporal Differences between Microbial Assemblages from Two Neighboring Sulfurous Lakes: Comparison by Microscopy and Denaturing Gradient Gel Electrophoresis

ABSTRACT The microbial assemblages of Lake Cisó and Lake Vilar (Banyoles, northeast Spain) were analyzed in space and time by microscopy and by performing PCR-denaturing gradient gel electrophoresis (DGGE) and sequence analysis of 16S rRNA gene fragments. Samples obtained from different water depths and at two different times of the year (in the winter during holomixis and in the early spring during a phytoplankton bloom) were analyzed. Although the lakes have the same climatic conditions and the same water source, the limnological parameters were different, as were most of the morphologically distinguishable photosynthetic bacteria enumerated by microscopy. The phylogenetic affiliations of the predominant DGGE bands were inferred by performing a comparative 16S rRNA sequence analysis. Sequences obtained from Lake Cisó samples were related to gram-positive bacteria and to members of the divisionProteobacteria. Sequences obtained from Lake Vilar samples were related to members of theCytophaga-Flavobacterium-Bacteroides phylum and to cyanobacteria. Thus, we found that like the previously reported differences between morphologically distinct inhabitants of the two lakes, there were also differences among the community members whose morphologies did not differ conspicuously. The changes in the species composition from winter to spring were also marked. The two lakes both contained sequences belonging to phototrophic green sulfur bacteria, which is consistent with microscopic observations, but these sequences were different from the sequences of cultured strains previously isolated from the lakes. Euryarchaeal sequences (i.e., methanogen- and thermoplasma-related sequences) also were present in both lakes. These euryarchaeal group sequences dominated the archaeal sequences in Lake Cisó but not in Lake Vilar. In Lake Vilar, a new planktonic population related to the crenarchaeota produced the dominant archaeal band. The phylogenetic analysis indicated that new bacterial and archaeal lineages were present and that the microbial diversity of these assemblages was greater than previously known. We evaluated the correspondence between the abundances of several morphotypes and DGGE bands by comparing microscopy and sequencing results. Our data provide evidence that the sequences obtained from the DGGE fingerprints correspond to the microorganisms that are actually present at higher concentrations in the natural system.

Microbial Ecology of an Extreme Acidic Environment, the Tinto River

ABSTRACT The Tinto River (Huelva, southwestern Spain) is an extreme environment with a rather constant acidic pH along the entire river and a high concentration of heavy metals. The extreme conditions of the Tinto ecosystem are generated by the metabolic activity of chemolithotrophic microorganisms thriving in the rich complex sulfides of the Iberian Pyrite Belt. Molecular ecology techniques were used to analyze the diversity of this microbial community. The communitys composition was studied by denaturing gradient gel electrophoresis (DGGE) using 16S rRNA and by 16S rRNA gene amplification. A good correlation between the two approaches was found. Comparative sequence analysis of DGGE bands showed the presence of organisms related to Leptospirillum spp., Acidithiobacillus ferrooxidans, Acidiphilium spp., “Ferrimicrobium acidiphilum,” Ferroplasma acidiphilum, and Thermoplasma acidophilum. The different phylogenetic groups were quantified by fluorescent in situ hybridization with a set of rRNA-targeted oligonucleotide probes. More than 80% of the cells were affiliated with the domain Bacteria, with only a minor fraction corresponding to Archaea. Members of Leptospirillum ferrooxidans, Acidithiobacillus ferrooxidans, and Acidiphilium spp., all related to the iron cycle, accounted for most of the prokaryotic microorganisms detected. Different isolates of these microorganisms were obtained from the Tinto ecosystem, and their physiological properties were determined. Given the physicochemical characteristics of the habitat and the physiological properties and relative concentrations of the different prokaryotes found in the river, a model for the Tinto ecosystem based on the iron cycle is suggested.

Bacterial Community Structure Associated with a Dimethylsulfoniopropionate-Producing North Atlantic Algal Bloom

ABSTRACT The bacteria associated with oceanic algal blooms are acknowledged to play important roles in carbon, nitrogen, and sulfur cycling, yet little information is available on their identities or phylogenetic affiliations. Three culture-independent methods were used to characterize bacteria from a dimethylsulfoniopropionate (DMSP)-producing algal bloom in the North Atlantic. Group-specific 16S rRNA-targeted oligonucleotides, 16S ribosomal DNA (rDNA) clone libraries, and terminal restriction fragment length polymorphism analysis all indicated that the marine Roseobacter lineage was numerically important in the heterotrophic bacterial community, averaging >20% of the 16S rDNA sampled. Two other groups of heterotrophic bacteria, the SAR86 and SAR11 clades, were also shown by the three 16S rRNA-based methods to be abundant in the bloom community. In surface waters, the Roseobacter, SAR86, and SAR11 lineages together accounted for over 50% of the bacterial rDNA and showed little spatial variability in abundance despite variations in the dominant algal species. Depth profiles indicated thatRoseobacter phylotype abundance decreased with depth and was positively correlated with chlorophyll a, DMSP, and total organic sulfur (dimethyl sulfide plus DMSP plus dimethyl sulfoxide) concentrations. Based on these data and previous physiological studies of cultured Roseobacter strains, we hypothesize that this lineage plays a role in cycling organic sulfur compounds produced within the bloom. Three other abundant bacterial phylotypes (representing a cyanobacterium and two members of the α Proteobacteria) were primarily associated with chlorophyll-rich surface waters of the bloom (0 to 50 m), while two others (representing Cytophagales and δProteobacteria) were primarily found in deeper waters (200 to 500 m).

Ecology of the rare microbial biosphere of the Arctic Ocean

Understanding the role of microbes in the oceans has focused on taxa that occur in high abundance; yet most of the marine microbial diversity is largely determined by a long tail of low-abundance taxa. This rare biosphere may have a cosmopolitan distribution because of high dispersal and low loss rates, and possibly represents a source of phylotypes that become abundant when environmental conditions change. However, the true ecological role of rare marine microorganisms is still not known. Here, we use pyrosequencing to describe the structure and composition of the rare biosphere and to test whether it represents cosmopolitan taxa or whether, similar to abundant phylotypes, the rare community has a biogeography. Our examination of 740,353 16S rRNA gene sequences from 32 bacterial and archaeal communities from various locations of the Arctic Ocean showed that rare phylotypes did not have a cosmopolitan distribution but, rather, followed patterns similar to those of the most abundant members of the community and of the entire community. The abundance distributions of rare and abundant phylotypes were different, following a log-series and log-normal model, respectively, and the taxonomic composition of the rare biosphere was similar to the composition of the abundant phylotypes. We conclude that the rare biosphere has a biogeography and that its tremendous diversity is most likely subjected to ecological processes such as selection, speciation, and extinction.

Global ecological patterns in uncultured Archaea.

We have applied a global analytical approach to uncultured Archaea that for the first time reveals well-defined community patterns along broad environmental gradients and habitat types. Phylogenetic patterns and the environmental factors governing the creation and maintenance of these patterns were analyzed for c. 2000 archaeal 16S rRNA gene sequences from 67 globally distributed studies. The sequences were dereplicated at 97% identity, grouped into seven habitat types, and analyzed with both Unifrac (to explore shared phylogenetic history) and multivariate regression tree (that considers the relative abundance of the lineages or taxa) approaches. Both phylogenetic and taxon-based approaches showed salinity and not temperature as one of the principal driving forces at the global scale. Hydrothermal vents and planktonic freshwater habitats emerged as the largest reservoirs of archaeal diversity and consequently are promising environments for the discovery of new archaeal lineages. Conversely, soils were more phylogenetically clustered and archaeal diversity was the result of a high number of closely related phylotypes rather than different lineages. Applying the ecological concept of ‘indicator species’, we detected up to 13 indicator archaeal lineages for the seven habitats prospected. Some of these lineages (that is, hypersaline MSBL1, marine sediment FCG1 and freshwater plSA1), for which ecological importance has remained unseen to date, deserve further attention as they represent potential key archaeal groups in terms of distribution and ecological processes. Hydrothermal vents held the highest number of indicator lineages, suggesting it would be the earliest habitat colonized by Archaea. Overall, our approach provided ecological support for the often arbitrary nomenclature within uncultured Archaea, as well as phylogeographical clues on key ecological and evolutionary aspects of archaeal biology.

DOES ECOSYSTEM SIZE DETERMINE AQUATIC BACTERIAL RICHNESS

With the advent of DNA-based molecular technologies, microbial ecologists now have the tools to test whether general ecological patterns apply to microorganisms. In this study, we selected 11 high-mountain lakes from Sierra Nevada (Spain) to test the predictions of island-biogeography theory in relation to ecosystem size and isolation, and to assess the influence of other factors (i.e., ecosystem productivity, resource richness, and biotic interactions) on bacterial community structure. Bacterial operational taxonomic units (OTUs), generated by denaturing-gradient gel electrophoresis of polymerase chain-reaction-amplified 16S rRNA genes, were used as a surrogate of predominant “biodiversity units.” OTU composition among lakes was heterogeneous, and the number of site-specific OTUs was near 50%. Lake remoteness did not affect the number of bacterial OTUs although the spatial distribution of the lakes significantly influenced bacterial composition. Lakes that were closer together had more similar bacterial ...

Distribution of prokaryotic genetic diversity in athalassohaline lakes of the Atacama Desert, Northern Chile

Athalassohaline lakes are inland saline aquatic environments with ionic proportions quite different from the dissolved salts in seawater. Prokaryotes inhabiting athalassohaline environments are poorly known and very few of such places have been surveyed for microbial diversity studies around the world. We analyzed the planktonic bacterial and archaeal assemblages inhabiting several of these evaporitic basins in a remote and vast area in northern Chile by PCR-denaturing gradient gel electrophoresis (DGGE) and sequencing of 16S rRNA gene fragments. Most systems were springs and athalassohaline ponds in different saltflats of the Atacama Desert region, including Salar de Llamará (in the Central Depression), Salar de Atacama (in the Pre-Andean Depression) and Salar de Ascotán (in the Altiplano). Overall, we analyzed more than 25 samples from 19 different environments with strong gradients of altitude, qualitative ionic compositions and UV influence. Between 4 and 25 well-defined DGGE bands were detected for Bacteria in each sample, whereas Archaea ranged between 1 and 5. Predominant DGGE bands (defined by intensity and frequency of appearance) were excised from the gel and sequenced. Bacterial assemblages were dominated by the Cytophaga-Flavobacterium-Bacteroides (CFB) phylum and a few Proteobacteria. There was a tendency for increasing contribution of CFB with higher salinities and altitude. Thus, CFB accounted for the major fraction of band intensity in the Ascotán samples and for lower percentages in Atacama and Llamará. When the distribution of particular CFB sequences was examined, there were several relatives of Psychroflexus torquis substituting each other as salinity changed in Ascotán. Another set of CFB sequences, very distantly related to Cytophaga marinovorus, was abundant in both Llamará and Atacama at salinities lower than 7%. Archaeal assemblages were dominated by uncultured haloarchaea distantly related to cultured strains mostly obtained from thalassohaline environments. Most of the archaeal sequences did not have a close match with environmental 16S rRNA genes deposited in the database either. Therefore, athalassohaline environments are excellent sources of new microorganisms different from their counterparts in thalassohaline sites and useful tools to relate microbial genetic diversity and environmental characteristics such as changes in salinity (both qualitative and quantitative) and altitude.

Microheterogeneity in 16S Ribosomal DNA-Defined Bacterial Populations from a Stratified Planktonic Environment Is Related to Temporal Changes and to Ecological Adaptations

ABSTRACT Temporal changes of the bacterioplankton from a meromictic lake (Lake Vilar, Banyoles, Spain) were analyzed with four culture-independent techniques: epifluorescence microscopy, PCR-denaturing gradient gel electrophoresis (DGGE) fingerprinting, fluorescence in situ whole-cell hybridization and flow cytometry sorting. Microscopically, blooms of one cyanobacterium (Synechococcus sp.-like), one green sulfur bacterium (Chlorobium phaeobacteroides-like), and one purple sulfur bacterium (Thiocystis minor-like) were observed at different depths and times. DGGE retrieved these populations and, additionally, populations related to the Cytophaga-Flavobacterium-Bacteroides phylum as predominant community members. The analyses of partial 16S ribosomal DNA sequences from the DGGE fingerprints (550 bp analyzed) revealed higher genetic diversity than expected from microscopic observation for most of these groups. Thus, the sequences of two Synechococcus spp. (both had a similarity of 97% to Synechococcus sp. strain PCC6307 in 16S rRNA), two Thiocystis spp. (similarities to Thiocystis minor of 93 and 94%, respectively), and three Cytophaga spp. (similarities to Cytophaga fermentans of 88 and 89% and to Cytophaga sp. of 93%, respectively) were obtained. The two populations of Synechococcus exhibited different pigment compositions and temporal distributions and their 16S rRNA sequences were 97.3% similar. The two Thiocystis populations differed neither in pigment composition nor in morphology, but their 16S rRNA sequences were only 92.3% similar and they also showed different distributions over time. Finally, two of the Cytophaga spp. showed 96.2% similarity between the 16S rRNA sequences, but one of them was found to be mostly attached to particles and only in winter. Thus, the identity of the main populations changed over time, but the function of the microbial guilds was maintained. Our data showed that temporal shifts in the identity of the predominant population is a new explanation for the environmental 16S rRNA microdiversity retrieved from microbial assemblages and support the hypothesis that clusters of closely related 16S rRNA environmental sequences may actually represent numerous closely related, yet ecologically distinct, populations.

Differential photoinhibition of bacterial and archaeal ammonia oxidation

Inhibition by light potentially influences the distribution of ammonia oxidizers in aquatic environments and is one explanation for nitrite maxima near the base of the euphotic zone of oceanic waters. Previous studies of photoinhibition have been restricted to bacterial ammonia oxidizers, rather than archaeal ammonia oxidizers, which dominate in marine environments. To compare the photoinhibition of bacterial and archaeal ammonia oxidizers, specific growth rates of two ammonia-oxidizing archaea (Nitrosopumilus maritimus and Nitrosotalea devanaterra) and bacteria (Nitrosomonas europaea and Nitrosospira multiformis) were determined at different light intensities under continuous illumination and light/dark cycles. All strains were inhibited by continuous illumination at the highest intensity (500 μE m(-2) s(-1)). At lower light intensities, archaeal growth was much more photosensitive than bacterial growth, with greater inhibition at 60 μE m(-2) s(-1) than at 15 μE m(-2) s(-1), where bacteria were unaffected. Archaeal ammonia oxidizers were also more sensitive to cycles of 8-h light/16-h darkness at two light intensities (60 and 15 μE m(-2) s(-1)) and, unlike bacterial strains, showed no evidence of recovery during dark phases. The findings provide evidence for niche differentiation in aquatic environments and reduce support for photoinhibition as an explanation of nitrite maxima in the ocean.