Meta-transcriptomics indicates biotic cross-tolerance in willow trees cultivated on petroleum hydrocarbon contaminated soil

Review: Current understanding on biological filtration for the removal of microcystins

Harmful algal blooms (HABs) have become a global problem not only in aquatic habitats but also in public health and safety due to the production of cyanotoxins as their secondary metabolites. Among the various identified cyanotoxin groups, microcystins (MCs) are one of the most prevalent cyanotoxin detected during HABs. Different strategies including advanced physical and chemical treatment processes have been developed to mitigate the threat of cyanotoxins in water utilities, but these have revealed certain limitations in terms of high operational costs, low removal efficacy, and harmful by-products formation. Recently, biological filtration systems (BFS) have gained attention for safe drinking water production as they can treat various natural organic matter (NOM) and emerging contaminants through a highly efficient and environmentally sustainable process. However, limited attention has been given to understand the current research progress, research challenges, and knowledge gaps for the successful implementation of BFS for MC removal. Therefore, in this review, currently identified MC biodegradation pathways and MC-degrading microorganisms with their degradation rates are summarized, which may be pivotal for studying bioaugmented BFS to enhance the MC removal during HABs. Moreover, both laboratory and field studies on BFS for MC removal are reviewed, followed by a discussion of current challenges and future research needs for the practical application of BFS.

Human milk microbiome is shaped by breastfeeding practices

There is evidence that breastfeeding practices may impact the milk microbiota diversity and differential abundance at the genera level; however, the possibility that distinct feeding practices, such as exclusive (EBF) and non-exclusive breastfeeding (non-EBF), might alter the milk microbiome at the species level has not been explored. This cross-sectional study analyzed the milk microbiome of 64 Mam-Mayan indigenous mothers from San Juan Ostuncalco in Guatemala. Two breastfeeding practices [exclusive (EBF) vs non-exclusive (non-EBF)] were analyzed at two stages of lactation [early (5–46 days post-partum) vs late (109–184 days post-partum)]. EBF was defined as offering only human milk and non-EBF was defined as feeding the infant herbal teas (agüitas) and/or complementary foods while continuing to breastfeed. Results identified four clusters with distinct microbial communities that segregated bacterial species by both breastfeeding practices and stage of lactation. Comparison among these clusters identified several notable patterns. First, during EBF, the microbiome differed by stage of lactation where there was a shift in differential abundance from Actinobacteria and Firmicutes in early to Bacteroidetes and Proteobacteria species in late lactation. Second, a similar comparison between non-EBF mothers by stage of lactation also identified a higher differential abundance of Actinobacteria and Firmicutes species in early lactation, but only Proteobacteria and not Bacteroidetes in late lactation, indicating a further shift in the milk microbial ecosystem with fewer oral bacteria present in late lactation. Third, comparisons between EBF and non-EBF mothers at both early and late lactation showed that mothers who exclusively breastfed had more differentially abundant species in early (11 vs 1) and late (13 vs 2) lactation. Fourth, EBF at early and late lactation had more commensal and lactic acid bacteria, including Lactobacillus gasseri, Granulicatella elegans, Streptococcus mitis, and Streptococcus parasanguinis, compared to those who did not exclusively breastfeed. Collectively, these results show that EBF has more differentially abundant bacteria, including commensal and lactic acid bacteria, and that the addition of agüitas (herbal teas) and/or complementary foods modify the milk microbiome composition by reducing the oral bacteria and introducing more environmentally sourced bacteria to the ecosystem.

Bacterial inoculant-assisted phytoremediation affects trace element uptake and metabolite content in Salix atrocinerea

Natural plant-associated microorganisms are of critical importance to plant growth and survival in field conditions under toxic concentrations of trace elements (TE) and these plant-microbial processes can be harnessed to enhance phytoremediation. The total bacterial diversity from grey willow (Salix atrocinerea) on a brownfield heavily-polluted with lead (Pb) and arsenic (As) was studied through pyrosequencing. Culturable bacteria were isolated and in vitro tested for plant growth-promotion (PGP) traits, arsenic (As) tolerance and impact on As speciation. Two of the most promising bacterial strains - the root endophyte Pantoea sp. AV62 and the rhizospheric strain Rhodococcus erythropolis AV96 - were inoculated in field to S. atrocinerea. This bioaugmentation resulted in higher As and Pb concentrations in both, roots and leaves of bacterial-inoculated plants as compared to non-inoculated plants. In consequence, bacterial bioaugmentation also affected parameters related to plant growth, oxidative stress, the levels of phytochelatins and phenylpropanoids, together with the differential expression of genes related to these tolerance mechanisms to TE in leaves. This study extends our understanding about plant-bacterial interactions and provides a solid basis for further bioaugmentation studies aiming to improve TE phytoremediation efficiency and predictability in the field.

Understanding Willow Transcriptional Response in the Context of Oil Sands Tailings Reclamation

One of the reclamation objectives for treated oil sands tailings (OST) is to establish boreal forest communities that can integrate with the surrounding area. Hence, selection of appropriate soil reclamation cover designs and plant species for revegetation are important aspects of tailings landform reclamation and closure. Research and monitoring of the long term and immediate impacts of capped OST on the growth and survival of native boreal plant species are currently underway. However, plant responses to OST-associated toxicity are not well known at the molecular level. Using RNA sequencing, we examined the effects of three types of OST on the willow transcriptome under different capping strategies. The transcriptomic data showed that some genes respond universally and others in a specific manner to different types of OST. Among the dominant and shared upregulated genes, we found some encoding protein detoxification (PD), Cytochrome P450 (CYPs), glutathione S-transferase regulatory process (GST), UDP-glycosyltransferase (UGT), and ABC transporter and regulatory process associated proteins. Moreover, genes encoding several stress-responsive transcription factors (bZIP, BHLH, ERF, MYB, NAC, WRKY) were upregulated with OST-exposure, while high numbers of transcripts related to photosynthetic activity and chloroplast structure and function were downregulated. Overall, the expression of 40 genes was found consistent across all tailings types and capping strategies. The qPCR analysis of a subset of these shared genes suggested that they could reliably distinguish plants exposed to different OST associated stress. Our results indicated that it is possible to develop OST stress exposure biosensors merely based on changes in the level of expression of a relatively small set of genes. The outcomes of this study will further guide optimization of OST capping and revegetation technology by using knowledge based plant stress adaptation strategies.

Stem metabolism under drought stress - a paradox of increasing respiratory substrates and decreasing respiratory rates

Metabolic changes underpinning drought-induced variations in stem respiration (R_s ) are unknown. We measured R_s rates and metabolite and gene expression profiles in Ulmus minor Mill. and Quercus ilex L. seedlings subjected to increasing levels of drought stress to better understand how carbon, nitrogen and energy metabolism interact during drought. In both species, only plants showing extreme stress symptoms - i.e. negligible rates of leaf stomatal conductance and photosynthesis, and high stem dehydration (30-50% of maximum water storage) and contraction (50-150 μm week^-1 ) - exhibited lower R_s rates than well-watered plants. Abundance of low-molecular weight sugars (e.g. glucose and fructose) and sugar alcohols (e.g. mannitol) increased with drought, at more moderate stress and to a higher extent in Q. ilex than U. minor. Abundance of amino acids increased at more severe stress, more abruptly, and to a higher extent in U. minor, coinciding with leaf senescence, which did not occur in Q. ilex. Organic acids changed less in response to drought: threonate and glycerate increased, and citrate decreased although slightly in both species. Transcripts of genes coding for enzymes of the Krebs cycle decreased in Q. ilex and increased in U. minor in conditions of extreme drought stress. The maintenance of R_s under severe growth and photosynthetic restrictions reveals the importance of stem mitochondrial activity in drought acclimation. The eventual decline in R_s diverts carbon substrates from entering the Krebs cycle that may help to cope with osmotic and oxidative stress during severe drought and to recover hydraulic functionality afterwards.

Arsenate-reducing bacteria affect As accumulation and tolerance in Salix atrocinerea

Arsenic (As)-reducing bacteria are able to influence As-speciation and, in this way, change As bio-availability. In consequence, this has an impact on As uptake by plants growing on polluted soil and on the effectiveness of the phytoremediation process. To be able to efficiently utilize these bacteria for As-phytoremediation in the field, a better understanding of the plant-bacterial interactions involved in As-tolerance or toxicity is needed. In this work, seedlings of a clone of Salix atrocinerea derived from a specimen naturally growing on an As-polluted brownfield were grown under gnotobiotic conditions exposed to As, and in the presence or absence of two of its field-associated and in vitro characterized plant growth-promoting (PGP) bacteria. The inoculation with Pantoea sp., induced a moderate reduction of AsV to AsIII in the exposure medium that, together with a coordinated plant response of As uptake, chelation and sequestration, increased As accumulation in roots; which is reflected into a higher phytostabilization. However, inoculation with Rhodococcus erythropolis due to a higher disproportionate reduction of AsV to AsIII in the medium caused less As accumulation in roots that non-bioaugmented plants and despite the lower As content, the concentrations of AsIII present in the medium and the damage suffered in roots and leaves, indicated that As tolerance mechanisms (such as prevention of AsIII uptake and efflux) did not occur in time to avoid physical disturbance and plants growth reduction. Interestingly, by two different metabolic pathways -coordinated by different key transporters mediating As uptake, tolerance, distribution and vacuolar accumulation at the roots- both bacteria limited As accumulation in Salix shoots. Our results provide for the first time a detailed insight in the plant-bacterial responses and physiological changes contributing to As tolerance in S. atrocinerea, that will facilitate the design of effective strategies for exploitation of plant-associated microorganisms for phytoremediation.

Human Milk Microbiota in an Indigenous Population Is Associated with Maternal Factors, Stage of Lactation, and Breastfeeding Practices

BACKGROUND

Human milk contains a diverse community of bacteria that are modified by maternal factors, but whether these or other factors are similar in developing countries has not been explored. Our objective was to determine whether the milk microbiota was modified by maternal age, BMI, parity, lactation stage, subclinical mastitis (SCM), and breastfeeding practices in the first 6 mo of lactation in an indigenous population from Guatemala.

METHODS

For this cross-sectional study, Mam-Mayan indigenous mothers nursing infants aged <6 mo were recruited. Unilateral human milk samples were collected (n = 86) and processed for 16S rRNA sequencing at the genus level. Microbial diversity and relative abundance were compared with maternal factors [age, BMI, parity, stage of lactation, SCM, and 3 breastfeeding practices (exclusive, predominant, mixed)] obtained through questionnaires.

RESULTS

Streptococcus was the most abundant genus (33.8%), followed by Pseudomonas (18.7%) and Sphingobium (10.7%) but relative abundance was associated with maternal factors. First, Lactobacillus and Streptococcus were more abundant in early lactation whereas the common oral (Leptotrichia) and environmental (Comamonas) bacteria were more abundant in established lactation. Second, Streptococcus,Lactobacillus,Lactococcus,Leuconostoc, and Micrococcus had a higher abundance in multiparous mothers compared with primiparous mothers. Third, a more diverse microbiota characterized by a higher abundance of lactic acid bacteria (Lactobacillus,Leuconostoc, and Lactococcus), Leucobacter, and Micrococcus was found in mothers with a healthy BMI. Finally, distinct microbial communities differed by stage of lactation and by exclusive, predominant, or mixed breastfeeding practices.

CONCLUSION

Milk bacterial communities in an indigenous community were associated with maternal factors. Higher microbial diversity was supported by having a healthy BMI, the absence of SCM, and by breastfeeding. Interestingly, breastfeeding practices when assessed by lactation stage were associated with distinct microbiota profiles.

Gene expression and biochemical response of giant reed under Ni and Cu stress

Giant reed (Arundo donax) has proved to be effective in detoxification, accumulation and tolerance of toxic metals. The present study explored the stress response of giant reed against Cu and Ni stress. The effect of metal stress was studied on dry weight, chlorophyll pigments antioxidant enzymes production and selected genes expression. The accumulation of heavy metals increased in a concentration-dependent manner and depicted toxicity symptoms in leaves beyond 75 mg/L of Cu or Ni. Oxidative stress was evident in giant reed under highest exposure of Ni and Cu which increased antioxidants activities (SOD, POD and CAT). It was observed that metal transport and detoxification were possible due to the expression of glutathione reductase, Natural Resistance-Associated Macrophage Protein (NRAMP) and Yellow Stripe-Like (YSL) genes. These insights into the genetic basis of a successful remediating plant species will be useful in understanding heavy metals tolerance in giant reed.

ANCHOR: a 16S rRNA gene amplicon pipeline for microbial analysis of multiple environmental samples

Analysis of 16S ribosomal RNA (rRNA) gene amplification data for microbial barcoding can be inaccurate across complex environmental samples. A method, ANCHOR, is presented and designed for improved species‐level microbial identification using paired‐end sequences directly, multiple high‐complexity samples and multiple reference databases. A standard operating procedure (SOP) is reported alongside benchmarking against artificial, single sample and replicated mock data sets. The method is then directly tested using a real‐world data set from surface swabs of the International Space Station (ISS). Simple mock community analysis identified 100% of the expected species and 99% of expected gene copy variants (100% identical). A replicated mock community revealed similar or better numbers of expected species than MetaAmp, DADA2, Mothur and QIIME1. Analysis of the ISS microbiome identified 714 putative unique species/strains and differential abundance analysis distinguished significant differences between the Destiny module (U.S. laboratory) and Harmony module (sleeping quarters). Harmony was remarkably dominated by human gastrointestinal tract bacteria, similar to enclosed environments on earth; however, Destiny module bacteria also derived from nonhuman microbiome carriers present on the ISS, the laboratory's research animals. ANCHOR can help substantially improve sequence resolution of 16S rRNA gene amplification data within biologically replicated environmental experiments and integrated multidatabase annotation enhances interpretation of complex, nonreference microbiomes.

Comparative transcriptome analysis of Zea mays in response to petroleum hydrocarbon stress

The use of plants for the improvement of soils contaminated with hydrocarbons has been a primary research focus in phytoremediation studies. Obtaining insights regarding genes that are differentially induced by petroleum hydrocarbon stress and understanding plant response mechanisms against petroleum hydrocarbons at molecular level is essential for developing better phytoremediation strategies to remove these hazardous contaminants. The purpose of this study was to analyze the transcriptomal profile changes under hydrocarbon stress in maize plants and identify the genes associated with the phytoremediative capacity. Zea mays GeneChips were used to analyze the global transcriptome profiles of maize treated with different concentrations of petroleum hydrocarbons. In total, 883, 1281, and 2162 genes were differentially induced or suppressed in the comparisons of 0 (control) vs. 1% crude petroleum, 1 vs. 5% crude petroleum, and 0 vs. 5% crude petroleum, respectively. The differentially expressed genes were functionally associated with the osmotic stress response mechanism, likely preventing the uptake of water from the roots, and the phytoremediative capacity of plants, e.g., secretory pathway genes. The results presented here show the regulatory mechanisms in the response to petroleum hydrocarbon pollution in soil. Our study provides global gene expression data of Z. mays in response to petroleum hydrocarbon stress that could be useful for further studies investigating the biodegradation mechanism in maize and other plants.

Trees, fungi and bacteria: tripartite metatranscriptomics of a root microbiome responding to soil contamination

BACKGROUND

One method for rejuvenating land polluted with anthropogenic contaminants is through phytoremediation, the reclamation of land through the cultivation of specific crops. The capacity for phytoremediation crops, such as Salix spp., to tolerate and even flourish in contaminated soils relies on a highly complex and predominantly cryptic interacting community of microbial life.

METHODS

Here, Illumina HiSeq 2500 sequencing and de novo transcriptome assembly were used to observe gene expression in washed Salix purpurea cv. 'Fish Creek' roots from trees pot grown in petroleum hydrocarbon-contaminated or non-contaminated soil. All 189,849 assembled contigs were annotated without a priori assumption as to sequence origin and differential expression was assessed.

RESULTS

The 839 contigs differentially expressed (DE) and annotated from S. purpurea revealed substantial increases in transcripts encoding abiotic stress response equipment, such as glutathione S-transferases, in roots of contaminated trees as well as the hallmarks of fungal interaction, such as SWEET2 (Sugars Will Eventually Be Exported Transporter). A total of 8252 DE transcripts were fungal in origin, with contamination conditions resulting in a community shift from Ascomycota to Basidiomycota genera. In response to contamination, 1745 Basidiomycota transcripts increased in abundance (the majority uniquely expressed in contaminated soil) including major monosaccharide transporter MST1, primary cell wall and lamella CAZy enzymes, and an ectomycorrhiza-upregulated exo-β-1,3-glucanase (GH5). Additionally, 639 DE polycistronic transcripts from an uncharacterised Enterobacteriaceae species were uniformly in higher abundance in contamination conditions and comprised a wide spectrum of genes cryptic under laboratory conditions but considered putatively involved in eukaryotic interaction, biofilm formation and dioxygenase hydrocarbon degradation.

CONCLUSIONS

Fungal gene expression, representing the majority of contigs assembled, suggests out-competition of white rot Ascomycota genera (dominated by Pyronema), a sometimes ectomycorrhizal (ECM) Ascomycota (Tuber) and ECM Basidiomycota (Hebeloma) by a poorly characterised putative ECM Basidiomycota due to contamination. Root and fungal expression involved transcripts encoding carbohydrate/amino acid (C/N) dialogue whereas bacterial gene expression included the apparatus necessary for biofilm interaction and direct reduction of contamination stress, a potential bacterial currency for a role in tripartite mutualism. Unmistakable within the metatranscriptome is the degree to which the landscape of rhizospheric biology, particularly the important but predominantly uncharacterised fungal genetics, is yet to be discovered.

Soil contamination alters the willow root and rhizosphere metatranscriptome and the root-rhizosphere interactome

Phytoremediation using willows is thought to be a sustainable alternative to traditional remediation techniques involving excavation, transport, and landfilling. However, the complexity of the interaction between the willow and its associated highly diverse microbial communities makes the optimization of phytoremediation very difficult. Here, we have sequenced the rhizosphere metatranscriptome of four willow species and the plant root metatranscriptome for two willow species growing in petroleum hydrocarbon-contaminated and non-contaminated soils on a former petroleum refinery site. Significant differences in the abundance of transcripts related to different bacterial and fungal taxa were observed between willow species, mostly in contaminated soils. When comparing transcript abundance in contaminated vs. non-contaminated soil for each willow species individually, transcripts for many microbial taxa and functions were significantly more abundant in contaminated rhizosphere soil for Salix eriocephala, S. miyabeana and S. purpurea, in contrast to what was observed in the rhizosphere of S. caprea. This agrees with the previously reported sensitivity of S. caprea to contamination, and the superior tolerance of S. miyabeana and S. purpurea to soil contamination at that site. The root metatranscriptomes of two species were compared and revealed that plants transcripts are mainly influenced by willow species, while microbial transcripts mainly responded to contamination. A comparison of the rhizosphere and root metatranscriptomes in the S. purpurea species revealed a complete reorganization of the linkages between root and rhizosphere pathways when comparing willows growing in contaminated and non-contaminated soils, mainly because of large shifts in the rhizosphere metatranscriptome.

Transcriptomic Response of Purple Willow (Salix purpurea) to Arsenic Stress

Arsenic (As) is a toxic element for plants and one of the most common anthropogenic pollutants found at contaminated sites. Despite its severe effects on plant metabolism, several species can accumulate substantial amounts of arsenic and endure the associated stress. However, the genetic mechanisms involved in arsenic tolerance remains obscure in many model plant species used for land decontamination (phytoremediation), including willows. The present study assesses the potential of Salix purpurea cv. 'Fish Creek' for arsenic phytoextraction and reveals the genetic responses behind arsenic tolerance, phytoextraction and metabolism. Four weeks of hydroponic exposure to 0, 5, 30 and 100 mg/L revealed that plants were able to tolerate up to 5 mg/L arsenic. Concentrations of 0 and 5 mg/L of arsenic treatment were then used to compare alterations in gene expression of roots, stems and leaves using RNA sequencing. Differential gene expression revealed transcripts encoding proteins putatively involved in entry of arsenic into the roots, storage in vacuoles and potential transport through the plant as well as primary and secondary (indirect) toxicity tolerance mechanisms. A major role for tannin as a compound used to relieve cellular toxicity is implicated as well as unexpected expression of the cadmium transporter CAX2, providing a potential means for internal arsenic mobility. These insights into the underpinning genetics of a successful phytoremediating species present novel opportunities for selection of dedicated arsenic tolerant crops as well as the potential to integrate such tolerances into a wider Salix ideotype alongside traits including biomass yield, biomass quality, low agricultural inputs and phytochemical production.

Comparative Transcriptomic Approaches Exploring Contamination Stress Tolerance in Salix sp. Reveal the Importance for a Metaorganismal de Novo Assembly Approach for Nonmodel Plants

Metatranscriptomic study of nonmodel organisms requires strategies that retain the highly resolved genetic information generated from model organisms while allowing for identification of the unexpected. A real-world biological application of phytoremediation, the field growth of 10 Salix cultivars on polluted soils, was used as an exemplar nonmodel and multifaceted crop response well-disposed to the study of gene expression. Sequence reads were assembled de novo to create 10 independent transcriptomes, a global transcriptome, and were mapped against the Salix purpurea 94006 reference genome. Annotation of assembled contigs was performed without a priori assumption of the originating organism. Global transcriptome construction from 3.03 billion paired-end reads revealed 606,880 unique contigs annotated from 1588 species, often common in all 10 cultivars. Comparisons between transcriptomic and metatranscriptomic methodologies provide clear evidence that nonnative RNA can mistakenly map to reference genomes, especially to conserved regions of common housekeeping genes, such as actin, α/β-tubulin, and elongation factor 1-α. In Salix, Rubisco activase transcripts were down-regulated in contaminated trees across all 10 cultivars, whereas thiamine thizole synthase and CP12, a Calvin Cycle master regulator, were uniformly up-regulated. De novo assembly approaches, with unconstrained annotation, can improve data quality; care should be taken when exploring such plant genetics to reduce de facto data exclusion by mapping to a single reference genome alone. Salix gene expression patterns strongly suggest cultivar-wide alteration of specific photosynthetic apparatus and protection of the antenna complexes from oxidation damage in contaminated trees, providing an insight into common stress tolerance strategies in a real-world phytoremediation system.