Functional diversity of microbial eukaryotes in a meromictic lake: Coupling between metatranscriptomic and a trait-based approach.

The advent of high-throughput sequencing has led to the discovery of a considerable diversity of microbial eukaryotes in aquatic ecosystems, nevertheless, their function and contribution to the trophic food web functioning remain poorly characterized especially in freshwater ecosystems. Based on metabarcoding data obtained from a meromictic lake ecosystem (Pavin, France), we performed a morpho-physio-phenological traits-based approach to infer functional groups of microbial eukaryotes. Metatranscriptomic data were also analysed to assess the metabolic potential of these groups across the diel cycle, size fraction, sampling depth, and periods. Our analysis highlights a huge microbial eukaryotic diversity in the monimolimnion characterized by numerous saprotrophs expressing transcripts related to sulfur and nitrate metabolism as well as dissolved and particulate organic matter degradation. We also describe strong seasonal variations of microbial eukaryotes in the mixolimnion, especially for parasites and mixoplankton. It appears that the water mixing (occurring during spring and autumn) which benefits photosynthetic host communities also promotes parasitic fungi dissemination and over-expression of genes involved in the zoospore phototaxis and stage transition in the parasitic cycle. Mixoplanktonic haptophytes over-expressing photosynthesis-, endocytosis- and phagosome-linked genes under nutrient limitation also suggest that phagotrophy may provide them an advantage over non-phagotrophic phytoplankton.

Reconstruction of plasmids by shotgun sequencing from environmental DNA: which bioinformatic workflow?

Plasmids play important roles in microbial evolution and also in the spread of antibiotic resistance. Plasmid sequences are extensively studied from clinical isolates but rarely from the environment with a metagenomic approach focused on the plasmid fraction referred to as the plasmidome. A clear challenge in this context is to define a workflow for discriminating plasmids from chromosomal contaminants existing in the plasmidome. For this purpose, we benchmarked existing tools from assembly to detection of the plasmids by reference-free methods (cBar and PlasFlow) and database-guided approaches. Our simulations took into account short-reads alone or combined with moderate long-reads like those actually generated in environmental genomics experiments. This benchmark allowed us to select the best tools for limiting false-positives associated to plasmid prediction tools and a combination of reference-guided methods based on plasmid and bacterial databases.

Temporal variations of Microsporidia diversity and discovery of new host-parasite interactions in a lake ecosystem.

Microsporidia are a large group of obligate intracellular eukaryotic parasites related to Fungi. Recent studies suggest that their diversity has been greatly underestimated and little is known about their hosts other than metazoans, and thus about their impact on the communities at the base of the food web. In this work, we therefore studied the diversity of Microsporidia over one year and identified potential new hosts in small-sized fractions (<150 µm) in a lake ecosystem using a metabarcoding approach coupled with co-occurrence networks and tyramide signal amplification-fluorescent in situ hybridization. Our analysis shows a great Microsporidia diversity (1 472 OTUs), with an important part of this diversity being unknown. Temporal variations of this diversity have been observed, which might follow temporal variations of their potential hosts such as protists and microzooplankton. New hosts among them were identified as well as associations with phytoplankton. Indeed, repeated infections were observed in Kellicottia (rotifers) with a prevalence of 38% (infected individuals). Microsporidia inside a Stentor (ciliate) were also observed. Finally, potential infections of the diatom Asterionella were identified (prevalence <0.1%). The microsporidian host spectrum could be therefore even more important than previously described, and their role in the functioning of lake ecosystems is undoubtedly largely unknown.

Plasmidome analysis of a hospital effluent biofilm: Status of antibiotic resistance.

Plasmids are widely involved in the dissemination of characteristics within bacterial communities. Their genomic content can be assessed by high-throughput sequencing of the whole plasmid fraction of an environment, the plasmidome. In this study, we analyzed the plasmidome of a biofilm formed in the effluents of the teaching hospital of Clermont-Ferrand (France). Our analysis discovered >350 new complete plasmids, with a length ranging from 1219 to 40,193 bp. Forty-two plasmid incompatibility (Inc) groups were found among all the plasmid contigs. Ten large plasmids, described here in detail, were reconstructed from plasmid contigs, seven of which carried antibiotic resistance genes. Four plasmids potentially confer resistance to numerous families of antibiotics, including carbapenems, aminoglycosides, colistin, and chloramphenicol. Most of these plasmids were affiliated to Proteobacteria, a phylum of Gram-negative bacteria. This study therefore illustrates the composition of an environmental mixed biofilm in terms of plasmids and antibiotic resistance genes.

A phylogenetic framework to investigate the microsporidian communities through metabarcoding and its application to lake ecosystems.

Microsporidia are obligate intracellular eukaryotic parasites known to parasitize many species of the animal kingdom as well as some protists. However, their diversity is underestimated, in part as a consequence of the failure of 'universal' primers to detect them in metabarcoding studies. Besides, due to the inconsistency between taxonomy and phylogenetic data, available databases may assign incorrectly sequences obtained with high-throughput sequencing. In this work, we developed a comprehensive reference database which positions microsporidian SSU rRNA gene sequences within a coherent ranked phylogenetic framework. We used this phylogenetic framework to study the microsporidian diversity in lacustrine ecosystems, focusing on < 150 µm planktonic size fractions. Our analysis shows a high diversity of Microsporidia, with the identification of 1531 OTUs distributed within seven clades, of which 76% were affiliated to clade IV2 and 20% to clade I (nomenclature presented hereby). About a quarter of the obtained sequences shared less than 85% identity to the closest known species, which might represent undescribed genera or families infecting small hosts. Variations in the abundance of Microsporidia were recorded between the two lakes sampled and across the sampling period, which might be explained by spatio-temporal variations of their potential hosts such as microeukaryotes and metazooplankton.

Honeybee gut microbiota dysbiosis in pesticide/parasite co-exposures is mainly induced by Nosema ceranae.

Honeybees ensure a key ecosystem service by pollinating many agricultural crops and wild plants. However, in the past few decades, managed bee colonies have been declining in Europe and North America. Researchers have emphasized both parasites and pesticides as the most important factors. Infection by the parasite Nosema ceranae and exposure to pesticides can contribute to gut dysbiosis, impacting the honeybee physiology. Here, we examined and quantified the effects of N. ceranae, the neonicotinoid thiamethoxam, the phenylpyrazole fipronil and the carboxamide boscalid, alone and in combination, on the honeybee gut microbiota. Chronic exposures to fipronil and thiamethoxam alone or combined with N. ceranae infection significantly decreased honeybee survival whereas the fungicide boscalid had no effect on uninfected bees. Interestingly, increased mortality was observed in N. ceranae-infected bees after exposure to boscalid, with synergistic negative effects. Regarding gut microbiota composition, co-exposure to the parasite and each pesticide led to decreased abundance of Alphaproteobacteria, and increased abundance of Gammaproteobacteria. The parasite also induced an increase of bacterial alpha-diversity (species richness). Our findings demonstrated that exposure of honeybees to N. ceranae and/or pesticides play a significant role in colony health and is associated with the establishment of a dysbiotic gut microbiota.

A Pediococcus strain to rescue honeybees by decreasing Nosema ceranae- and pesticide-induced adverse effects.

Honeybees ensure a key ecosystemic service by pollinating many agricultural crops and wild plants. However, since few decades, managed bee colonies have declined worldwide. This phenomenon is considered to be multifactorial, with a strong emphasis on both parasites and pesticides. Infection by the parasite Nosema ceranae and exposure to pesticides can contribute to adverse effects, resulting in a perturbation of the honeybee physiology. We thus hypothesized that probiotic treatment could be promising to treat or prevent these disturbances. The aim of this study was to evaluate the effects of probiotics on N. ceranae-infected and intoxicated honeybees (by the insecticide thiamethoxam and the fungicide boscalid). For this purpose, experiments were conducted with five probiotics. Among them, Pediococcus acidilactici (PA) showed the best protective effect against the parasite and pesticides. PA significantly improved the infected honeybee lifespan as prophylactic and curative treatments (respectively 2.3 fold and 1.7 fold). Furthermore, the exposure to pesticides induced an increase of honeybee mortality compared with the control group (p < .001) that was restored by the PA treatment. Despite its beneficial effect on honeybee lifespan, the PA administration did not induce changes in the gut bacterial communities (neither in abundance or diversity). N. ceranae and the pesticides were shown to deregulate genes involved in honeybee development (vitellogenin), immunity (serine protease 40, defensin) and detoxification system (glutathione peroxidase-like 2, catalase), and these effects were corrected by the PA treatment. This study highlights the promising use of PA to protect honeybees from both pathogens and pesticides.

Temporal covariation of epibacterial community and surface metabolome in the Mediterranean seaweed holobiont Taonia atomaria.

An integrative multi-omics approach allowed monthly variations for a year of the surface metabolome and the epibacterial community of the Mediterranean Phaeophyceae Taonia atomaria to be investigated. The LC-MS-based metabolomics and 16S rDNA metabarcoding data sets were integrated in a multivariate meta-omics analysis (multi-block PLS-DA from the MixOmic DIABLO analysis) showing a strong seasonal covariation (Mantel test: p < 0.01). A network based on positive and negative correlations between the two data sets revealed two clusters of variables, one relative to the 'spring period' and a second to the 'summer period'. The 'spring period' cluster was mainly characterized by dipeptides positively correlated with a single bacterial taxon of the Alteromonadaceae family (BD1-7 clade). Moreover, 'summer' dominant epibacterial taxa from the second cluster (including Erythrobacteraceae, Rhodospirillaceae, Oceanospirillaceae and Flammeovirgaceae) showed positive correlations with few metabolites known as macroalgal antifouling defences [e.g. dimethylsulphoniopropionate (DMSP) and proline] which exhibited a key role within the correlation network. Despite a core community that represents a significant part of the total epibacteria, changes in the microbiota structure associated with surface metabolome variations suggested that both environment and algal host shape the bacterial surface microbiota.

Monitoring human enteric viruses in wastewater and relevance to infections encountered in the clinical setting: a one-year experiment in central France, 2014 to 2015.

BackgroundHuman enteric viruses are resistant in the environment and transmitted via the faecal-oral route. Viral shedding in wastewater gives the opportunity to track emerging pathogens and study the epidemiology of enteric infectious diseases in the community. Aim: The aim of this study was to monitor the circulation of enteric viruses in the population of the Clermont-Ferrand area (France) by analysis of urban wastewaters. Methods: Raw and treated wastewaters were collected between October 2014 and October 2015 and concentrated by a two-step protocol using tangential flow ultrafiltration and polyethylene glycol precipitation. Processed samples were analysed for molecular detection of adenovirus, norovirus, rotavirus, parechovirus, enterovirus (EV), hepatitis A (HAV) and E (HEV) viruses. Results: All wastewater samples (n = 54) contained viruses. On average, six and four virus species were detected in, respectively, raw and treated wastewater samples. EV-positive samples were tested for EV-D68 to assess its circulation in the community. EV-D68 was detected in seven of 27 raw samples. We collected data from clinical cases of EV-D68 (n = 17), HAV (n = 4) and HEV infection (n = 16) and compared wastewater-derived sequences with clinical sequences. We showed the silent circulation of EV-D68 in September 2015, the wide circulation of HAV despite few notifications of acute disease and the presence in wastewater of the major HEV subtypes involved in clinical local cases. Conclusion: The environmental surveillance overcomes the sampling bias intrinsic to the study of infections associated with hospitalisation and allows the detection in real time of viral sequences genetically close to those reported in clinical specimens.

Tracking a century of changes in microbial eukaryotic diversity in lakes driven by nutrient enrichment and climate warming.

High-throughput sequencing of sedimentary DNA (sed-DNA) was utilized to reconstruct the temporal dynamics of microbial eukaryotic communities (MECs) at a centennial scale in two re-oligotrophicated lakes that were exposed to different levels of phosphorus enrichment. The temporal changes within the MECs were expressed in terms of richness, composition and community structure to investigate their relationships with two key forcing factors (i.e., nutrient enrichment and climate warming). Various groups, including Apicomplexa, Cercozoa, Chrysophyceae, Ciliophora, Chlorophyceae and Dinophyceae, responded to phosphorus enrichment levels with either positive or negative impacts on their richness and relative abundance. For both lakes, statistical modelling demonstrated that phosphorus concentration ([P]) was a dominant contributor to MECs modifications before the 1980s; after the mid-80s, the contribution of air temperature changes increased and potentially surpassed the contribution of [P]. Co-occurrence network analysis revealed that some clusters of taxa (i.e., modules) composed mainly of Dinophyceae and unclassified Alveolata were strongly correlated to air temperature in both lakes. Overall, our data showed that sed-DNA constitutes a precious archive of information on past biodiversity changes, allowing the study of the dynamics of numerous eukaryotic groups that were not traditionally considered in paleo-reconstructions.

Spatio-Temporal Variations of Marine Biofilm Communities Colonizing Artificial Substrata Including Antifouling Coatings in Contrasted French Coastal Environments.

Surface colonization in seawater first corresponds to the selection of specific microbial biofilm communities. By coupling flow cytometry, microscopy and high throughput sequencing (HTS, 454 pyrosequencing) with artificial surfaces and environmental analyses, we intend to identify the contribution of biofilm community drivers at two contrasted French sites, one temperate and eutrophic (Lorient, Atlantic coast) and the other at a mesotrophic but highly contaminated bay (Toulon, North-Western Mediterranean Sea). Microbial communities were shaped by high temperatures, salinity and lead at Toulon by but nutrients and DOC at Lorient. Coatings including pyrithione exhibited a significant decrease of their microbial densities except for nanoeukaryotes. Clustering of communities was mainly based on the surface type and secondly the site, whereas seasons appeared of less importance. The in-depth HTS revealed that γ- and α-proteobacteria, but also Bacteroidetes, dominated highly diversified bacterial communities with a relative low ß-diversity. Sensitivity to biocides released by the tested antifouling coatings could be noticed at different taxonomic levels: the percentage of Bacteroidetes overall decreased with the presence of pyrithione, whereas the α/γ-proteobacteria ratio decreased at Toulon when increased at Lorient. Small diatom cells (Amphora and Navicula spp.) dominated on all surfaces, whereas site-specific sub-dominant taxa appeared clearly more sensitive to biocides. This overall approach exhibited the critical significance of surface characteristics in biofilm community shaping.

Analysis of metagenomic data reveals common features of halophilic viral communities across continents.

Microbial communities from hypersaline ponds, dominated by halophilic archaea, are considered specific of such extreme conditions. The associated viral communities have accordingly been shown to display specific features, such as similar morphologies among different sites. However, little is known about the genetic diversity of these halophilic viral communities across the Earth. Here, we studied viral communities in hypersaline ponds sampled on the coast of Senegal (8-36% of salinity) using metagenomics approach, and compared them with hypersaline viromes from Australia and Spain. The specificity of hyperhalophilic viruses could first be demonstrated at a community scale, salinity being a strong discriminating factor between communities. For the major viral group detected in all samples (Caudovirales), only a limited number of halophilic Caudovirales clades were highlighted. These clades gather viruses from different continents and display consistent genetic composition, indicating that they represent related lineages with a worldwide distribution. Non-tailed hyperhalophilic viruses display a greater rate of gene transfer and recombination, with uncharacterized genes conserved across different kind of viruses and plasmids. Thus, hypersaline viral communities around the world appear to form a genetically consistent community that are likely to harbour new genes coding for enzymes specifically adapted to these environments.

Long-term dynamics in microbial eukaryotes communities: a palaeolimnological view based on sedimentary DNA.

Assessing the extent to which changes in lacustrine biodiversity are affected by anthropogenic or climatic forces requires extensive palaeolimnological data. We used high-throughput sequencing to generate time-series data encompassing over 2200 years of microbial eukaryotes (protists and Fungi) diversity changes from the sedimentary DNA record of two lakes (Lake Bourget in French Alps and Lake Igaliku in Greenland). From 176 samples, we sequenced a large diversity of microbial eukaryotes, with a total 16 386 operational taxonomic units distributed within 50 phylogenetic groups. Thus, microbial groups, such as Chlorophyta, Dinophyceae, Haptophyceae and Ciliophora, that were not previously considered in lacustrine sediment record analyses appeared to be potential biological markers of trophic status changes. Our data suggest that shifts in relative abundance of extant species, including shifts between rare and abundant taxa, drive ecosystem responses to local and global environmental changes. Community structure shift events were concomitant with major climate variations (more particularly in Lake Igaliku). However, this study shows that the impacts of climatic fluctuations may be overpassed by the high-magnitude eutrophication impacts, as observed in the eutrophicated Lake Bourget. Overall, our data show that DNA preserved in sediment constitutes a precious archive of information on past biodiversity changes.

Structure of the rare archaeal biosphere and seasonal dynamics of active ecotypes in surface coastal waters.

Marine Archaea are important players among microbial plankton and significantly contribute to biogeochemical cycles, but details regarding their community structure and long-term seasonal activity and dynamics remain largely unexplored. In this study, we monitored the interannual archaeal community composition of abundant and rare biospheres in northwestern Mediterranean Sea surface waters by pyrosequencing 16S rDNA and rRNA. A detailed analysis of the rare biosphere structure showed that the rare archaeal community was composed of three distinct fractions. One contained the rare Archaea that became abundant at different times within the same ecosystem; these cells were typically not dormant, and we hypothesize that they represent a local seed bank that is specific and essential for ecosystem functioning through cycling seasonal environmental conditions. The second fraction contained cells that were uncommon in public databases and not active, consisting of aliens to the studied ecosystem and representing a nonlocal seed bank of potential colonizers. The third fraction contained Archaea that were always rare but actively growing; their affiliation and seasonal dynamics were similar to the abundant microbes and could not be considered a seed bank. We also showed that the major archaeal groups, Thaumarchaeota marine group I and Euryarchaeota group II.B in winter and Euryarchaeota group II.A in summer, contained different ecotypes with varying activities. Our findings suggest that archaeal diversity could be associated with distinct metabolisms or life strategies, and that the rare archaeal biosphere is composed of a complex assortment of organisms with distinct histories that affect their potential for growth.

Evidence for an active rare biosphere within freshwater protists community.

Studies on the active rare biosphere at the RNA level are mainly focused on Bacteria and Archaea and fail to include the protists, which are involved in the main biogeochemical cycles of the earth. In this study, the richness, composition and activity of the rare protistan biosphere were determined from a temporal survey of two lakes by pyrosequencing. In these ecosystems, the always rare OTUs represented 77.2% of the total OTUs and 76.6% of the phylogenetic diversity. From the various phylogenetic indices computed, the phylogenetic units (PUs) constituted exclusively by always rare OTUs were discriminated from the other PUs. Therefore, the rare biosphere included mainly taxa that are distant from the reference databases compared to the dominant ones. In addition, the rarest OTUs represented 59.8% of the active biosphere depicted by rRNA and the activity (rRNA:rDNA ratio) increased with the rarity. The high rRNA:rDNA ratio determined in the rare fraction highlights that some protists were active at low abundances and contribute to ecosystem functioning. Interestingly, the always rare and active OTUs were characterized by seasonal changes in relation with the main environmental parameters measured. In conclusion, the rare eukaryotes represent an active, dynamic and overlooked fraction in the lacustrine ecosystems.

Technical challenges in metatranscriptomic studies applied to the bacterial communities of freshwater ecosystems.

Metatranscriptome analysis relates to the transcriptome of microbial communities directly sampled in the environment. Accessing the mRNA pool in natural bacterial communities presents some technical challenges such as the RNA extraction, rRNA depletion, and the choice of the high-throughput sequencing technique. The lack of technical details in scientific articles is a major problem to correctly obtained mRNA from a microbial community and thus the corresponding sequencing data. In our study, we present the methodological procedure that was developed in order to access to the metatranscriptome of the microbial communities during two cyanobacterial blooms successively occurring in a freshwater eutrophic lake. Each procedure step was detailed and discussed with regard to the choices and difficulties encountered and to the recent literature. Finally, the two major limits for metatranscriptomic approaches targeting bacterial communities from natural environments were (i) the removal of rRNA in order to increase the putative mRNA reads number after sequencing, and (ii) for most of the bacterial communities living in natural environments, the lack of reference genomes in databases that leads to the non-assignation of numerous reads. Once these challenges overcome, we managed to access putative mRNA of dominant species, i.e. cyanobacteria (from 6 to 72 % of mRNA assigned), and of the surrounding bacteria (from 1 to 5 % of mRNA assigned).

Is Planktonic Diversity Well Recorded in Sedimentary DNA? Toward the Reconstruction of Past Protistan Diversity.

Studies based on the coupling of a paleolimnological approach and molecular tools (e.g., sequencing of sedimentary DNA) present a promising opportunity to obtain long-term data on past lacustrine biodiversity. However, certain validations are still required, such as the evaluation of DNA preservation in sediments for various planktonic taxa that do not leave any morphological diagnostic features. In this study, we focused on the diversity of planktonic unicellular eukaryotes and verified the presence of their DNA in sediment archives. We compared the molecular inventories (high-throughput sequencing of 18S ribosomal DNA) obtained from monitoring the water column with those obtained for DNA archived in the first 30 cm of sediment. Seventy-one percent of taxonomic units found in the water samples were detected in sediment samples, including pigmented taxa, such as Chlorophyta, Dinophyceae, and Chrysophyceae, phagotrophic taxa, such as Ciliophora, parasitic taxa, such as Apicomplexa and Chytridiomycota, and saprotrophs, such as Cryptomycota. Parallel analysis of 18S ribosomal RNA (rRNA) transcripts revealed the presence of living eukaryotic taxa only in the top 2 cm of sediment; although some limits exist in using RNA/DNA ratio as indicator of microbial activity, these results suggested that the sedimentary DNA mostly represented DNA from past and inactive communities. Only the diversity of a few groups, such as Cryptophyta and Haptophyta, seemed to be poorly preserved in sediments. Our overall results showed that the application of sequencing techniques to sedimentary DNA could be used to reconstruct past diversity for numerous planktonic eukaryotic groups.

Temporal Dynamics of Active Prokaryotic Nitrifiers and Archaeal Communities from River to Sea.

To test if different niches for potential nitrifiers exist in estuarine systems, we assessed by pyrosequencing the diversity of archaeal gene transcript markers for taxonomy (16S ribosomal RNA (rRNA)) during an entire year along a salinity gradient in surface waters of the Charente estuary (Atlantic coast, France). We further investigated the potential for estuarine prokaryotes to oxidize ammonia and hydrolyze urea by quantifying thaumarchaeal amoA and ureC and bacterial amoA transcripts. Our results showed a succession of different nitrifiers from river to sea with bacterial amoA transcripts dominating in the freshwater station while archaeal transcripts were predominant in the marine station. The 16S rRNA sequence analysis revealed that Thaumarchaeota marine group I (MGI) were the most abundant overall but other archaeal groups like Methanosaeta were also potentially active in winter (December-March) and Euryarchaeota marine group II (MGII) were dominant in seawater in summer (April-August). Each station also contained different Thaumarchaeota MGI phylogenetic clusters, and the clusters' microdiversity was associated to specific environmental conditions suggesting the presence of ecotypes adapted to distinct ecological niches. The amoA and ureC transcript dynamics further indicated that some of the Thaumarchaeota MGI subclusters were involved in ammonia oxidation through the hydrolysis of urea. Our findings show that ammonia-oxidizing Archaea and Bacteria were adapted to contrasted conditions and that the Thaumarchaeota MGI diversity probably corresponds to distinct metabolisms or life strategies.

Metavir 2: new tools for viral metagenome comparison and assembled virome analysis.

BACKGROUND: Metagenomics, based on culture-independent sequencing, is a well-fitted approach to provide insights into the composition, structure and dynamics of environmental viral communities. Following recent advances in sequencing technologies, new challenges arise for existing bioinformatic tools dedicated to viral metagenome (i.e. virome) analysis as (i) the number of viromes is rapidly growing and (ii) large genomic fragments can now be obtained by assembling the huge amount of sequence data generated for each metagenome. RESULTS: To face these challenges, a new version of Metavir was developed. First, all Metavir tools have been adapted to support comparative analysis of viromes in order to improve the analysis of multiple datasets. In addition to the sequence comparison previously provided, viromes can now be compared through their k-mer frequencies, their taxonomic compositions, recruitment plots and phylogenetic trees containing sequences from different datasets. Second, a new section has been specifically designed to handle assembled viromes made of thousands of large genomic fragments (i.e. contigs). This section includes an annotation pipeline for uploaded viral contigs (gene prediction, similarity search against reference viral genomes and protein domains) and an extensive comparison between contigs and reference genomes. Contigs and their annotations can be explored on the website through specifically developed dynamic genomic maps and interactive networks. CONCLUSIONS: The new features of Metavir 2 allow users to explore and analyze viromes composed of raw reads or assembled fragments through a set of adapted tools and a user-friendly interface.