Toward molecular trait-based ecology through integration of biogeochemical, geographical and metagenomic data

Using metagenomic ‘parts lists' to study microbial ecology remains a significant challenge. This work proposes a molecular trait-based approach to biogeography by integrating metagenomic data with external metadata and using functional community composition as readout

Local Expansion of a Panmictic Lineage of Water Bloom-Forming Cyanobacterium Microcystis aeruginosa

In previous studies, we have demonstrated that the population structure of the bloom-forming cyanobacterium Microcystis aeruginosa is clonal. Expanded multilocus sequence typing analysis of M. aeruginosa using 412 isolates identified five intraspecific lineages suggested to be panmictic while maintaining overall clonal structure probably due to a reduced recombination rate between lineages. Interestingly, since 2005 most strains belonging to one of these panmictic clusters (group G) have been found in a particular locality (Lake Kasumigaura Basin) in Japan. In this locality, multiple, similar but distinct genotypes of this lineage predominated in the bloom, a pattern that is unprecedented for M. aeruginosa. The population structure underlying blooms associated with this lineage is comparable to epidemics of pathogens. Our results may reveal an expansion of the possible adaptive lineage in a localized aquatic environment, providing us with a unique opportunity to investigate its ecological and biogeographical consequences

Bacterial diversity and community composition from seasurface to subseafloor

© The International Society for Microbial Ecology, 2015. This article is distributed under the terms of the Creative Commons Attribution License. The definitive version was published in ISME Journal 10 (2016): 979–989, doi:10.1038/ismej.2015.175.We investigated compositional relationships between bacterial communities in the water column and those in deep-sea sediment at three environmentally distinct Pacific sites (two in the Equatorial Pacific and one in the North Pacific Gyre). Through pyrosequencing of the v4–v6 hypervariable regions of the 16S ribosomal RNA gene, we characterized 450 104 pyrotags representing 29 814 operational taxonomic units (OTUs, 97% similarity). Hierarchical clustering and non-metric multidimensional scaling partition the samples into four broad groups, regardless of geographic location: a photic-zone community, a subphotic community, a shallow sedimentary community and a subseafloor sedimentary community (greater than or equal to1.5 meters below seafloor). Abundance-weighted community compositions of water-column samples exhibit a similar trend with depth at all sites, with successive epipelagic, mesopelagic, bathypelagic and abyssopelagic communities. Taxonomic richness is generally highest in the water-column O2 minimum zone and lowest in the subseafloor sediment. OTUs represented by abundant tags in the subseafloor sediment are often present but represented by few tags in the water column, and represented by moderately abundant tags in the shallow sediment. In contrast, OTUs represented by abundant tags in the water are generally absent from the subseafloor sediment. These results are consistent with (i) dispersal of marine sedimentary bacteria via the ocean, and (ii) selection of the subseafloor sedimentary community from within the community present in shallow sediment.This study was funded by the Biological Oceanography Program of the US National Science Foundation (grant OCE-0752336) and by the NSF-funded Center for Dark Energy Biosphere Investigations (grant NSF-OCE-0939564)

Ocean currents shape the microbiome of Arctic marine sediments

Prokaryote communities were investigated on the seasonally stratified Alaska Beaufort Shelf (ABS). Water and sediment directly underlying water with origin in the Arctic, Pacific or Atlantic oceans were analyzed by pyrosequencing and length heterogeneity-PCR in conjunction with physicochemical and geographic distance data to determine what features structure ABS microbiomes. Distinct bacterial communities were evident in all water masses. Alphaproteobacteria explained similarity in Arctic surface water and Pacific derived water. Deltaproteobacteria were abundant in Atlantic origin water and drove similarity among samples. Most archaeal sequences in water were related to unclassified marine Euryarchaeota. Sediment communities influenced by Pacific and Atlantic water were distinct from each other and pelagic communities. Firmicutes and Chloroflexi were abundant in sediment, although their distribution varied in Atlantic and Pacific influenced sites. Thermoprotei dominated archaea in Pacific influenced sediments and Methanomicrobia dominated in methane-containing Atlantic influenced sediments. Length heterogeneity-PCR data from this study were analyzed with data from methane-containing sediments in other regions. Pacific influenced ABS sediments clustered with Pacific sites from New Zealand and Chilean coastal margins. Atlantic influenced ABS sediments formed another distinct cluster. Density and salinity were significant structuring features on pelagic communities. Porosity co-varied with benthic community structure across sites and methane did not. This study indicates that the origin of water overlying sediments shapes benthic communities locally and globally and that hydrography exerts greater influence on microbial community structure than the availability of methane

Preservation of organic matter in sediments promoted by iron

The biogeochemical cycles of iron and organic carbon are strongly interlinked. In oceanic waters, organic ligands have been shown to control the concentration of dissolved iron. In soils, solid iron phases shelter and preserve organic carbon, but the role of iron in the preservation of organic matter in sediments has not been clearly established. Here we use an iron reduction method previously applied to soils to determine the amount of organic carbon associated with reactive iron phases in sediments of various mineralogies collected from a wide range of depositional environments. Our findings suggest that 21.5 ± 8.6 per cent of the organic carbon in sediments is directly bound to reactive iron phases. We further estimate that a global mass of (19–45) × 1015 grams of organic carbon is preserved in surface marine sediments as a result of its association with iron. We propose that these associations between organic carbon and iron, which are formed primarily through co-precipitation and/or direct chelation, promote the preservation of organic carbon in sediments. Because reactive iron phases are metastable over geological timescales, we suggest that they serve as an efficient ‘rusty sink’ for organic carbon, acting as a key factor in the long-term storage of organic carbon and thus contributing to the global cycles of carbon, oxygen and sulphur

Taxa-area relationship of aquatic fungi on deciduous leaves

One of the fundamental patterns in macroecology is the increase in the number of observed taxa with size of sampled area. For microbes, the shape of this relationship remains less clear. The current study assessed the diversity of aquatic fungi, by the traditional approach based on conidial morphology (captures reproducing aquatic hyphomycetes) and next generation sequencing (NGS; captures other fungi as well), on graded sizes of alder leaves (0.6 to 13.6 cm2). Leaves were submerged in two streams in geographically distant locations: the Oliveira Stream in Portugal and the Boss Brook in Canada. Decay rates of alder leaves and fungal sporulation rates did not differ between streams. Fungal biomass was higher in Boss Brook than in Oliveira Stream, and in both streams almost 100% of the reads belonged to active fungal taxa. In general, larger leaf areas tended to harbour more fungi, but these findings were not consistent between techniques. Morphospecies-based diversity increased with leaf area in Boss Brook, but not in Oliveira Stream; metabarcoding data showed an opposite trend. The higher resolution of metabarcoding resulted in steeper taxa-accumulation curves than morphospecies-based assessments (fungal conidia morphology). Fungal communities assessed by metabarcoding were spatially structured by leaf area in both streams. Metabarcoding promises greater resolution to assess biodiversity patterns in aquatic fungi and may be more accurate for assessing taxa-area relationships and local to global diversity ratios.This work was supported by the strategic programme UID/BIA/04050/2013 (POCI-01-0145-FEDER-007569), funded by national funds through the Portuguese Foundation for Science and Technology (FCT) I.P. (http://www.fct.pt/) and by the ERDF through the COMPETE2020 - Programa Operacional Competitividade e Internacionalizacao (POCI) and by the project PTDC/AAC-AMB/117068/2010, funded by national funds through FCT I.P. and the European Regional Development Funds through the Operational Competitiveness Program (FEDER-COMPETE). Support from FCT to SD (SFRH/BPD/47574/2008 and SFRH/BPD/109842/2015) and from NSERC Discovery grant program (http://www.nserc-crsng.gc.ca/index_eng.asp) to FB is also acknowledged. The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript.info:eu-repo/semantics/publishedVersio

Phylogeography of the Microcoleus vaginatus (Cyanobacteria) from Three Continents – A Spatial and Temporal Characterization

It has long been assumed that cyanobacteria have, as with other free-living microorganisms, a ubiquitous occurrence. Neither the geographical dispersal barriers nor allopatric speciation has been taken into account. We endeavoured to examine the spatial and temporal patterns of global distribution within populations of the cyanobacterium Microcoleus vaginatus, originated from three continents, and to evaluate the role of dispersal barriers in the evolution of free-living cyanobacteria. Complex phylogeographical approach was applied to assess the dispersal and evolutionary patterns in the cyanobacterium Microcoleus vaginatus (Oscillatoriales). We compared the 16S rRNA and 16S-23S ITS sequences of strains which had originated from three continents (North America, Europe, and Asia). The spatial distribution was investigated using a phylogenetic tree, network, as well as principal coordinate analysis (PCoA). A temporal characterization was inferred using molecular clocks, calibrated from fossil DNA. Data analysis revealed broad genetic diversity within M. vaginatus. Based on the phylogenetic tree, network, and PCoA analysis, the strains isolated in Europe were spatially separated from those which originated from Asia and North America. A chronogram showed a temporal limitation of dispersal barriers on the continental scale. Dispersal barriers and allopatric speciation had an important role in the evolution of M. vaginatus. However, these dispersal barriers did not have a permanent character; therefore, the genetic flow among populations on a continental scale was only temporarily present. Furthermore, M. vaginatus is a recently evolved species, which has been going through substantial evolutionary changes

Korarchaeota Diversity, Biogeography, and Abundance in Yellowstone and Great Basin Hot Springs and Ecological Niche Modeling Based on Machine Learning

Over 100 hot spring sediment samples were collected from 28 sites in 12 areas/regions, while recording as many coincident geochemical properties as feasible (>60 analytes). PCR was used to screen samples for Korarchaeota 16S rRNA genes. Over 500 Korarchaeota 16S rRNA genes were screened by RFLP analysis and 90 were sequenced, resulting in identification of novel Korarchaeota phylotypes and exclusive geographical variants. Korarchaeota diversity was low, as in other terrestrial geothermal systems, suggesting a marine origin for Korarchaeota with subsequent niche-invasion into terrestrial systems. Korarchaeota endemism is consistent with endemism of other terrestrial thermophiles and supports the existence of dispersal barriers. Korarchaeota were found predominantly in >55°C springs at pH 4.7–8.5 at concentrations up to 6.6×106 16S rRNA gene copies g−1 wet sediment. In Yellowstone National Park (YNP), Korarchaeota were most abundant in springs with a pH range of 5.7 to 7.0. High sulfate concentrations suggest these fluids are influenced by contributions from hydrothermal vapors that may be neutralized to some extent by mixing with water from deep geothermal sources or meteoric water. In the Great Basin (GB), Korarchaeota were most abundant at spring sources of pH<7.2 with high particulate C content and high alkalinity, which are likely to be buffered by the carbonic acid system. It is therefore likely that at least two different geological mechanisms in YNP and GB springs create the neutral to mildly acidic pH that is optimal for Korarchaeota. A classification support vector machine (C-SVM) trained on single analytes, two analyte combinations, or vectors from non-metric multidimensional scaling models was able to predict springs as Korarchaeota-optimal or sub-optimal habitats with accuracies up to 95%. To our knowledge, this is the most extensive analysis of the geochemical habitat of any high-level microbial taxon and the first application of a C-SVM to microbial ecology

A population genomics approach shows widespread geographical distribution of cryptic genomic forms of the symbiotic fungus Rhizophagus irregularis.

Arbuscular mycorrhizal fungi (AMF; phylum Gomeromycota) associate with plants forming one of the most successful microbe-plant associations. The fungi promote plant diversity and have a potentially important role in global agriculture. Plant growth depends on both inter- and intra-specific variation in AMF. It was recently reported that an unusually large number of AMF taxa have an intercontinental distribution, suggesting long-distance gene flow for many AMF species, facilitated by either long-distance natural dispersal mechanisms or human-assisted dispersal. However, the intercontinental distribution of AMF species has been questioned because the use of very low-resolution markers may be unsuitable to detect genetic differences among geographically separated AMF, as seen with some other fungi. This has been untestable because of the lack of population genomic data, with high resolution, for any AMF taxa. Here we use phylogenetics and population genomics to test for intra-specific variation in Rhizophagus irregularis, an AMF species for which genome sequence information already exists. We used ddRAD sequencing to obtain thousands of markers distributed across the genomes of 81 R. irregularis isolates and related species. Based on 6 888 variable positions, we observed significant genetic divergence into four main genetic groups within R. irregularis, highlighting that previous studies have not captured underlying genetic variation. Despite considerable genetic divergence, surprisingly, the variation could not be explained by geographical origin, thus also supporting the hypothesis for at least one AMF species of widely dispersed AMF genotypes at an intercontinental scale. Such information is crucial for understanding AMF ecology, and how these fungi can be used in an environmentally safe way in distant locations