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

Expression Quantitative Trait Locus of Wood Formation-Related Genes in Salix suchowensis

Shrub willows are widely planted for landscaping, soil remediation, and biomass production, due to their rapid growth rates. Identification of regulatory genes in wood formation would provide clues for genetic engineering of willows for improved growth traits on marginal lands. Here, we conducted an expression quantitative trait locus (eQTL) analysis, using a full sibling F1 population of Salix suchowensis, to explore the genetic mechanisms underlying wood formation. Based on variants identified from simplified genome sequencing and gene expression data from RNA sequencing, 16,487 eQTL blocks controlling 5505 genes were identified, including 2148 cis-eQTLs and 16,480 trans-eQTLs. eQTL hotspots were identified, based on eQTL frequency in genomic windows, revealing one hotspot controlling genes involved in wood formation regulation. Regulatory networks were further constructed, resulting in the identification of key regulatory genes, including three transcription factors (JAZ1, HAT22, MYB36) and CLV1, BAM1, CYCB2;4, CDKB2;1, associated with the proliferation and differentiation activity of cambium cells. The enrichment of genes in plant hormone pathways indicates their critical roles in the regulation of wood formation. Our analyses provide a significant groundwork for a comprehensive understanding of the regulatory network of wood formation in S. suchowensis.

Soil fauna-microbial interactions shifts fungal and bacterial communities under a contamination disturbance

The aim of this study was to determine whether the soil faunal-microbial interaction complexity (SFMIC) is a significant factor influencing the soil microbial communities and the willow growth in the context of PAH contamination. The SFMIC treatment had eight levels: just the microbial community, or the microbial community with nematodes, springtails, earthworms and all the possible combinations. SFMIC affected the height and biomass of willows after eight weeks or growth. SFMIC affected the structure and the composition of the bacterial, archaeal and fungal communities, with significant effects of SFMIC on the relative abundance of fungal genera such as Sphaerosporella, a known willow symbiont during phytoremediation, and bacterial phyla such as Actinobacteriota, containing many polycyclic aromatic hydrocarbons (PAH) degraders. These SFMIC effects on microbial communities were not clearly reflected in the community structure and abundance of PAH degraders, even though some degraders related to Actinobacteriota and the diversity of Gram-negative degraders were affected by the SFMIC treatments. Over 95% of PAH was degraded in all pots at the end of the experiment. Overall, our results suggest that, under our experimental conditions, SFMIC changes willow phytoremediation outcomes.

The Evolutionary History, Diversity, and Ecology of Willows (Salix L.) in the European Alps

The genus Salix (willows), with 33 species, represents the most diverse genus of woody plants in the European Alps. Many species dominate subalpine and alpine types of vegetation. Despite a long history of research on willows, the evolutionary and ecological factors for this species richness are poorly known. Here we will review recent progress in research on phylogenetic relationships, evolution, ecology, and speciation in alpine willows. Phylogenomic reconstructions suggest multiple colonization of the Alps, probably from the late Miocene onward, and reject hypotheses of a single radiation. Relatives occur in the Arctic and in temperate Eurasia. Most species are widespread in the European mountain systems or in the European lowlands. Within the Alps, species differ ecologically according to different elevational zones and habitat preferences. Homoploid hybridization is a frequent process in willows and happens mostly after climatic fluctuations and secondary contact. Breakdown of the ecological crossing barriers of species is followed by introgressive hybridization. Polyploidy is an important speciation mechanism, as 40% of species are polyploid, including the four endemic species of the Alps. Phylogenomic data suggest an allopolyploid origin for all taxa analyzed so far. Further studies are needed to specifically analyze biogeographical history, character evolution, and genome evolution of polyploids.

Transcriptome Analysis of Melocactus glaucescens (Cactaceae) Reveals Metabolic Changes During in vitro Shoot Organogenesis Induction

Melocactus glaucescens is an endangered cactus highly valued for its ornamental properties. In vitro shoot production of this species provides a sustainable alternative to overharvesting from the wild; however, its propagation could be improved if the genetic regulation underlying its developmental processes were known. The present study generated de novo transcriptome data, describing in vitro shoot organogenesis induction in M. glaucescens. Total RNA was extracted from explants before (control) and after shoot organogenesis induction (treated). A total of 14,478 unigenes (average length, 520 bases) were obtained using Illumina HiSeq 3000 (Illumina Inc., San Diego, CA, USA) sequencing and transcriptome assembly. Filtering for differential expression yielded 2,058 unigenes. Pairwise comparison of treated vs. control genes revealed that 1,241 (60.3%) unigenes exhibited no significant change, 226 (11%) were downregulated, and 591 (28.7%) were upregulated. Based on database analysis, more transcription factor families and unigenes appeared to be upregulated in the treated samples than in controls. Expression of WOUND INDUCED DEDIFFERENTIATION 1 (WIND1) and CALMODULIN (CaM) genes, both of which were upregulated in treated samples, was further validated by real-time quantitative PCR (RT-qPCR). Differences in gene expression patterns between control and treated samples indicate substantial changes in the primary and secondary metabolism of M. glaucescens after the induction of shoot organogenesis. These results help to clarify the molecular genetics and functional genomic aspects underlying propagation in the Cactaceae family.

High biomass yield increases in a primary effluent wastewater phytofiltration are associated to altered leaf morphology and stomatal size in Salix miyabeana

Municipal wastewater treatment using willow 'phyto'-filtration has the potential for reduced environmental impact compared to conventional treatment practices. However, the physiological adaptations underpinning tolerance to high wastewater irrigation in willow are unknown. A one-hectare phytofiltration plantation established using the Salix miyabeana cultivar 'SX67' in Saint-Roch-de-l'Achigan, Quebec, Canada, tested the impact of unirrigated, potable water or two loads of primary effluent wastewater 19 and 30 ML ha^-1 yr^-1. A nitrogen load of 817 kg N ha^-1 from wastewater did not increase soil pore water nitrogen concentrations beyond Quebec drinking water standards. The willow phytofiltration phenotype had increased leaf area (+106-142%) and leaf nitrogen (+94%) which were accompanied by significant increases in chlorophyll a + b content. Wastewater irrigated trees had higher stomatal sizes and a higher stomatal pore index, despite lower stomatal density, resulting in increased stomatal conductance (+42-78%). These developmental responses led to substantial increases in biomass yields of 56-207% and potable water controls revealed the nitrogen load to be necessary for the high productivity of 28-40 t ha^-1 yr^-1 in wastewater irrigated trees. Collectively, this study suggests phytofiltration plantations could treat primary effluent municipal wastewater at volumes of at least 19 million litres per hectare and benefit from increased yields of sustainable biomass over a two-year coppice cycle. Added-value cultivation practices, such as phytofiltration, have the potential to mitigate negative local and global environmental impact of wastewater treatment while providing valuable services and sustainable bioproducts.

Co-cropping with three phytoremediation crops influences rhizosphere microbiome community in contaminated soil

Human industrial activities have left millions of hectares of land polluted with trace element metals and persistent organic pollutants (POPs) around the world. Although contaminated sites are environmentally damaging, high economic costs often discourage soil remediation efforts. Phytoremediation is a potential green technology solution but can be challenging due to the diversity of anthropogenic contaminants. Co-cropping could provide improved tolerance to diverse soil challenges by taking advantage of distinct crop capabilities. Co-cropping of three species with potentially complementary functions, Festuca arundinacea, Salix miyabeana and Medicago sativa, perform well on diversely contaminated soils. Here, rhizosphere microbiomes of each crop in monoculture and in all co-cropping combinations were compared using 16S rRNA gene amplification, sequencing and differential abundance analysis. The hyperaccumulating F. arundinacea rhizosphere microbiome included putative plant growth promoting bacteria (PGPB) and metal tolerance species, such as Rhizorhapis suberifaciens, Cellvibrio fibrivorans and Pseudomonas lini. The rhizosphere microbiome of the fast-growing tree S. miyabeana included diverse taxa involved in POP degradation, including the species Phenylobacterium panacis. The well-characterised nitrogen-fixing M. sativa microbiome species, Sinorhizobium meliloti, was identified alongside others involved in nutrient acquisition and putative yet-to-be-cultured Candidatus saccharibacteria (TM7-1 group). The majority of differentially abundant rhizosphere-associated bacterial species were maintained in co-cropping pairs, with pairs having higher numbers of differentially abundant taxa than monocultures in all cases. This was not the case when all three crops were co-cropped, where most host-specific bacterial species were not detected as differentially abundant, indicating the potential for reduced rhizosphere functionality. The crops cultivated in pairs here retained rhizosphere microbiome bacteria involved in these monoculture ecosystem services of plant growth promotion, POP tolerance and degradation, and improved nutrient acquisition. These findings provide a promising outlook of the potential for complementary co-cropping strategies for phytoremediation of the multifaceted anthropogenic pollution which can disastrously affect soils around the world.

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 genomics and transcriptomics analysis reveals evolution patterns of selection in the Salix phylogeny

BACKGROUND

Willows are widely distributed in the northern hemisphere and have good adaptability to different living environment. The increasing of genome and transcriptome data provides a chance for comparative analysis to study the evolution patterns with the different origin and geographical distributions in the Salix phylogeny.

RESULTS

Transcript sequences of 10 Salicaceae species were downloaded from public databases. All pairwise of orthologues were identified by comparative analysis in these species, from which we constructed a phylogenetic tree and estimated the rate of diverse. Divergence times were estimated in the 10 Salicaceae using comparative transcriptomic analysis. All of the fast-evolving positive selection sequences were identified, and some cold-, drought-, light-, universal-, and heat- resistance genes were discovered.

CONCLUSIONS

The divergence time of subgenus Vetrix and Salix was about 17.6-16.0 Mya during the period of Middle Miocene Climate Transition (21-14 Mya). Subgenus Vetrix diverged to migratory and resident groups when the climate changed to the cool and dry trend by 14 Mya. Cold- and light- stress genes were involved in positive selection among the resident Vetrix, and which would help them to adapt the cooling stage. Universal- stress genes exhibited positive selection among the migratory group and subgenus Salix. These data are useful for comprehending the adaptive evolution and speciation in the Salix lineage.

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.