Dual RNA-seq reveals large-scale non-conserved genotype × genotype-specific genetic reprograming and molecular crosstalk in the mycorrhizal symbiosis

Unraveling plant-microbe symbioses using single-cell and spatial transcriptomics

Plant-microbe symbioses require intense interaction and genetic coordination to successfully establish in specific cell types of the host and symbiont. Traditional RNA-seq methodologies lack the cellular resolution to fully capture these complexities, but single-cell and spatial transcriptomics (ST) are now allowing scientists to probe symbiotic interactions at an unprecedented level of detail. Here, we discuss the advantages that novel spatial and single-cell transcriptomic technologies provide in studying plant-microbe endosymbioses and highlight key recent studies. Finally, we consider the remaining limitations of applying these approaches to symbiosis research, which are mainly related to the simultaneous capture of both plant and microbial transcripts within the same cells.

The mycorrhizal symbiosis: research frontiers in genomics, ecology, and agricultural application

Mycorrhizal symbioses between plants and fungi are vital for the soil structure, nutrient cycling, plant diversity, and ecosystem sustainability. More than 250 000 plant species are associated with mycorrhizal fungi. Recent advances in genomics and related approaches have revolutionized our understanding of the biology and ecology of mycorrhizal associations. The genomes of 250+ mycorrhizal fungi have been released and hundreds of genes that play pivotal roles in regulating symbiosis development and metabolism have been characterized. rDNA metabarcoding and metatranscriptomics provide novel insights into the ecological cues driving mycorrhizal communities and functions expressed by these associations, linking genes to ecological traits such as nutrient acquisition and soil organic matter decomposition. Here, we review genomic studies that have revealed genes involved in nutrient uptake and symbiosis development, and discuss adaptations that are fundamental to the evolution of mycorrhizal lifestyles. We also evaluated the ecosystem services provided by mycorrhizal networks and discuss how mycorrhizal symbioses hold promise for sustainable agriculture and forestry by enhancing nutrient acquisition and stress tolerance. Overall, unraveling the intricate dynamics of mycorrhizal symbioses is paramount for promoting ecological sustainability and addressing current pressing environmental concerns. This review ends with major frontiers for further research.

Dual RNA-Seq Reveals Temperature-Mediated Gene Reprogramming and Molecular Crosstalk between Grapevine and Lasiodiplodia theobromae

High temperatures associated with a fluctuating climate profoundly accelerate the occurrence of a myriad of plant diseases around the world. A comprehensive insight into how plants respond to pathogenic microorganisms under high-temperature stress is required for plant disease management, whereas the underlying mechanisms behind temperature-mediated plant immunity and pathogen pathogenicity are still unclear. Here, we evaluated the effect of high temperature on the development of grapevine canker disease and quantified the contribution of temperature variation to the gene transcription reprogramming of grapevine and its pathogenic agent Lasiodiplodia theobromae using a dual RNA-seq approach. The results showed that both grapevine and the pathogen displayed altered transcriptomes under different temperatures, and even the transcription of a plethora of genes from the two organisms responded in different directions and magnitudes. The transcription variability that arose due to temperature oscillation allowed us to identify a total of 26 grapevine gene modules and 17 fungal gene modules that were correlated with more than one gene module of the partner organism, which revealed an extensive web of plant-pathogen gene reprogramming during infection. More importantly, we identified a set of temperature-responsive genes that were transcriptionally orchestrated within the given gene modules. These genes are predicted to be involved in multiple cellular processes including protein folding, stress response regulation, and carbohydrate and peptide metabolisms in grapevine and porphyrin- and pteridine-containing compound metabolisms in L. theobromae, implying that in response to temperature oscillation, a complex web of signaling pathways in two organism cells is activated during infection. This study describes a co-transcription network of grapevine and L. theobromae in the context of considering temperature variation, which provides novel insights into deciphering the molecular mechanisms underlying temperature-modulated disease development.

Differential Responses of Medicago truncatula NLA Homologs to Nutrient Deficiency and Arbuscular Mycorrhizal Symbiosis

NITROGEN LIMITATION ADAPTATION (NLA), a plasma-membrane-associated ubiquitin E3 ligase, plays a negative role in the control of the phosphate transporter family 1 (PHT1) members in Arabidopsis and rice. There are three NLA homologs in the Medicago truncatula genome, but it has been unclear whether the function of these homologs is conserved in legumes. Here we investigated the subcellular localization and the responses of MtNLAs to external phosphate and nitrate status. Similar to AtNLA1, MtNLA1/MtNLA2 was localized in the plasma membrane and nucleus. MtNLA3 has three alternative splicing variants, and intriguingly, MtNLA3.1, the dominant variant, was not able to target the plasma membrane, whereas MtNLA3.2 and MtNLA3.3 were capable of associating with the plasma membrane. In contrast with AtNLA1, we found that MtNLAs were not affected or even upregulated by low-phosphate treatment. We also found that MtNLA3 was upregulated by arbuscular mycorrhizal (AM) symbiosis, and overexpressing MtNLA3.1 in Medicago roots resulted in a decrease in the transcription levels of STR, an essential gene for arbuscule development. Taken together, our results highlight the difference between MtNLA homologs and AtNLA1. Further characterization will be required to reveal the regulation of these genes and their roles in the responses to external nutrient status and AM symbiosis.

Variation of growth and transcriptome responses to arbuscular mycorrhizal symbiosis in different foxtail millet lines

BACKGROUND

Arbuscular mycorrhizal fungi (AMF) have been applied to promote the growth of different crop species, but knowledge about the impacts of symbiosis on foxtail millet at the physiological and molecular levels have remained limited. In this study, we compared the mycorrhization phenotypes of one cultivar and three different landraces and performed a comprehensive transcriptomic analysis to assess the effects of genetic variation on the responses to symbiosis.

RESULTS

Our results showed that colonization by AMF did not enhance biomass accumulation but significantly increased grain production only in three lines. More than 2,000 genes were affected by AMF colonization in all lines. Most AM symbiosis-conserved genes were induced, but the induction levels varied between lines. Gene Ontology (GO) analysis showed that Biological Function terms related to nitrogen transport and assimilation were only enriched in TT8. Similarly, two of phosphate starvation-induced phosphate transporters were only simultaneously downregulated in TT8. In the other two lines, the enrichment of GO terms associated with cell wall reorganization and lignification was observed, though the effects were different.

CONCLUSION

This study reveals the impacts of genetic variation of millet lines on the responses to AM symbiosis and provides information regarding AMF application for millet production.

Cross-tissue omics analysis discovers ten adipose genes encoding secreted proteins in obesity-related non-alcoholic fatty liver disease

BACKGROUND

Non-alcoholic fatty liver disease (NAFLD) is a fast-growing, underdiagnosed, epidemic. We hypothesise that obesity-related inflammation compromises adipose tissue functions, preventing efficient fat storage, and thus driving ectopic fat accumulation into the liver.

METHODS

To identify adipose-based mechanisms and potential serum biomarker candidates (SBCs) for NAFLD, we utilise dual-tissue RNA-sequencing (RNA-seq) data in adipose tissue and liver, paired with histology-based NAFLD diagnosis, from the same individuals in a cohort of obese individuals. We first scan for genes that are differentially expressed (DE) for NAFLD in obese individuals' subcutaneous adipose tissue but not in their liver; encode proteins secreted to serum; and show preferential adipose expression. Then the identified genes are filtered to key adipose-origin NAFLD genes by best subset analysis, knockdown experiments during human preadipocyte differentiation, recombinant protein treatment experiments in human liver HepG2 cells, and genetic analysis.

FINDINGS

We discover a set of genes, including 10 SBCs, that may modulate NAFLD pathogenesis by impacting adipose tissue function. Based on best subset analysis, we further follow-up on two SBCs CCDC80 and SOD3 by knockdown in human preadipocytes and subsequent differentiation experiments, which show that they modulate crucial adipogenesis genes, LPL, SREBPF1, and LEP. We also show that treatment of the liver HepG2 cells with the CCDC80 and SOD3 recombinant proteins impacts genes related to steatosis and lipid processing, including PPARA, NFE2L2, and RNF128. Finally, utilizing the adipose NAFLD DE gene cis-regulatory variants associated with serum triglycerides (TGs) in extensive genome-wide association studies (GWASs), we demonstrate a unidirectional effect of serum TGs on NAFLD with Mendelian Randomization (MR) analysis. We also demonstrate that a single SNP regulating one of the SBC genes, rs2845885, produces a significant MR result by itself. This supports the conclusion that genetically regulated adipose expression of the NAFLD DE genes may contribute to NAFLD through changes in serum TG levels.

INTERPRETATION

Our results from the dual-tissue transcriptomics screening improve the understanding of obesity-related NAFLD by providing a targeted set of 10 adipose tissue-active genes as new serum biomarker candidates for the currently grossly underdiagnosed fatty liver disease.

FUNDING

The work was supported by NIH grants R01HG010505 and R01DK132775. The Genotype-Tissue Expression (GTEx) Project was supported by the Common Fund of the Office of the Director of the National Institutes of Health, and by NCI, NHGRI, NHLBI, NIDA, NIMH, and NINDS. The KOBS study (J. P.) was supported by the Finnish Diabetes Research Foundation, Kuopio University Hospital Project grant (EVO/VTR grants 2005-2019), and the Academy of Finland grant (Contract no. 138006). This study was funded by the European Research Council under the European Union's Horizon 2020 research and innovation program (Grant No. 802825 to M. U. K.). K. H. P. was funded by the Academy of Finland (grant numbers 272376, 266286, 314383, and 335443), the Finnish Medical Foundation, Gyllenberg Foundation, Novo Nordisk Foundation (grant numbers NNF10OC1013354, NNF17OC0027232, and NNF20OC0060547), Finnish Diabetes Research Foundation, Finnish Foundation for Cardiovascular Research, University of Helsinki, and Helsinki University Hospital and Government Research Funds. I. S. was funded by the Instrumentarium Science Foundation. Personal grants to U. T. A. were received from the Matti and Vappu Maukonen Foundation, Ella och Georg Ehrnrooths Stiftelse and the Finnish Foundation for Cardiovascular Research.

Ovarian transcriptional response to Wolbachia infection in D. melanogaster in the context of between-genotype variation in gene expression

Wolbachia is a maternally transmitted endosymbiotic bacteria that infects a wide variety of arthropod and nematode hosts. The effects of Wolbachia on host biology are far-reaching and include changes in host gene expression. However, previous work on the host transcriptional response has generally been investigated in the context of a single host genotype. Thus, the relative effect of Wolbachia infection versus vs. host genotype on gene expression is unknown. Here, we explicitly test the relative roles of Wolbachia infection and host genotype on host gene expression by comparing the ovarian transcriptomes of 4 strains of Drosophila melanogaster (D. melanogaster) infected and uninfected with Wolbachia. Our data suggest that infection explains a small amount of transcriptional variation, particularly in comparison to variation in gene expression among strains. However, infection specifically affects genes related to cell cycle, translation, and metabolism. We also find enrichment of cell division and recombination processes among genes with infection-associated differential expression. Broadly, the transcriptomic changes identified in this study provide novel understanding of the relative magnitude of the effect of Wolbachia infection on gene expression in the context of host genetic variation and also point to genes that are consistently differentially expressed in response to infection among multiple genotypes.

Analysis of the molecular and biochemical mechanisms involved in the symbiotic relationship between Arbuscular mycorrhiza fungi and Manihot esculenta Crantz

INTRODUCTION

Plants and arbuscular mycorrhizal fungi (AMF) mutualistic interactions are essential for sustainable agriculture production. Although it is shown that AMF inoculation improves cassava physiological performances and yield traits, the molecular mechanisms involved in AM symbiosis remain largely unknown. Herein, we integrated metabolomics and transcriptomics analyses of symbiotic (Ri) and asymbiotic (CK) cassava roots and explored AM-induced biochemical and transcriptional changes.

RESULTS

Three weeks (3w) after AMF inoculations, proliferating fungal hyphae were observable, and plant height and root length were significantly increased. In total, we identified 1,016 metabolites, of which 25 were differentially accumulated (DAMs) at 3w. The most highly induced metabolites were 5-aminolevulinic acid, L-glutamic acid, and lysoPC 18:2. Transcriptome analysis identified 693 and 6,481 differentially expressed genes (DEGs) in the comparison between CK (3w) against Ri at 3w and 6w, respectively. Functional enrichment analyses of DAMs and DEGs unveiled transport, amino acids and sugar metabolisms, biosynthesis of secondary metabolites, plant hormone signal transduction, phenylpropanoid biosynthesis, and plant-pathogen interactions as the most differentially regulated pathways. Potential candidate genes, including nitrogen and phosphate transporters, transcription factors, phytohormone, sugar metabolism-related, and SYM (symbiosis) signaling pathway-related, were identified for future functional studies.

DISCUSSION

Our results provide molecular insights into AM symbiosis and valuable resources for improving cassava production.

Synergetic effect of water deficit and arbuscular mycorrhizal symbiosis on the expression of aquaporins in wheat (Triticum aestivum L.) roots: insights from NGS RNA-sequencing

Wheat (Triticum aestivum) is one of the most important crops in the world. This investigation was attempted to evaluate the transcriptional responses of aquaporins (AQPs) to the mycorrhizal inoculation and/or water deficit conditions in wheat to clarify how the arbuscular mycorrhizal symbiosis can contribute to the modulation of water homeostasis. The wheat seedlings were subjected to the water deficiency, and mycorrhizal inoculation using arbuscular fungus Funneliformis mosseae and Illumina RNA-Seq analyses confirmed that aquaporins expressed differentially in response to both the irrigation levels and mycorrhizal colonization. Results of this study showed that only 13% of the studied AQPs were responsive to water deficit with a tiny fraction (3%) being up-regulated. Mycorrhizal inoculation had a greater impact on the expression of AQPs with ca. 26% being responsive, ca. 4% of which were up-regulated. The samples with arbuscular mycorrhizal inoculation yielded more root and stem biomass. Water deficit and mycorrhizal inoculation caused different AQPs to be up-regulated. The effect of mycorrhizal inoculation on the expression of AQPs was intensified by applying water deficiency with 32% of studied AQPs being responsive, 6% of which up-regulated. We also found that the overexpression of three genes TaNIP1-10, TaNIP3-3, and TaNIP3-4 was chiefly triggered by mycorrhizal inoculation. Our results show that water deficit has a lower impact on the expression of aquaporins compared to what the arbuscular mycorrhizal inoculation has; water deficit and arbuscular mycorrhizal inoculation mainly cause the down-regulation of the aquaporins, and water deficit and the arbuscular inoculation have synergetic effects. These findings could improve our knowledge of how arbuscular mycorrhizal symbiosis can contribute to the modulation of water homeostasis.

Supplementary Information

The online version contains supplementary material available at 10.1007/s12298-023-01285-w.

Stochastic nuclear organization and host-dependent allele contribution in Rhizophagus irregularis

BACKGROUND

Arbuscular mycorrhizal (AM) fungi are arguably the most important symbionts of plants, offering a range of benefits to their hosts. However, the provisioning of these benefits does not appear to be uniform among AM fungal individuals, with genetic variation between fungal symbionts having a substantial impact on plant performance. Interestingly, genetic variation has also been reported within fungal individuals, which contain millions of haploid nuclei sharing a common cytoplasm. In the model AM fungus, Rhizophagus irregularis, several isolates have been reported to be dikaryotes, containing two genetically distinct types of nuclei recognized based on their mating-type (MAT) locus identity. However, their extremely coenocytic nature and lack of a known single nucleus stage has raised questions on the origin, distribution and dynamics of this genetic variation.

RESULTS

Here we performed DNA and RNA sequencing at the mycelial individual, single spore and single nucleus levels to gain insight into the dynamic genetic make-up of the dikaryote-like R. irregularis C3 isolate and the effect of different host plants on its genetic variation. Our analyses reveal that parallel spore and root culture batches can have widely variable ratios of two main genotypes in C3. Additionally, numerous polymorphisms were found with frequencies that deviated significantly from the general genotype ratio, indicating a diverse population of slightly different nucleotypes. Changing host plants did not show consistent host effects on nucleotype ratio's after multiple rounds of subculturing. Instead, we found a major effect of host plant-identity on allele-specific expression in C3.

CONCLUSION

Our analyses indicate a highly dynamic/variable genetic organization in different isolates of R. irregularis. Seemingly random fluctuations in nucleotype ratio's upon spore formation, recombination events, high variability of non-tandemly repeated rDNA sequences and host-dependent allele expression all add levels of variation that may contribute to the evolutionary success of these widespread symbionts.

Dual RNA and 16S ribosomal DNA sequencing reveal arbuscular mycorrhizal fungi-mediated mitigation of selenate stress in Zea mays L. and reshaping of soil microbiota

Excessively high concentrations of selenium (Se) in soil are toxic to crop plants, and inoculation with arbuscular mycorrhizal fungi (AMF) can reverse Se stress in maize (Zea mays L.). To investigate the underlying mechanisms, maize seedlings were treated with sodium selenate (5 mg Se[VI] kg^-1) and/or AMF (Funneliformis mosseae and Claroideoglomus etunicatum). Dual RNA sequencing in mycorrhiza and 16 S ribosomal DNA sequencing in soil were performed. The results showed that Se(VI) application alone decreased plant dry weight, but increased plant Se concentration, total Se content (mainly selenocysteine), and root superoxide content. Inoculation with either F. mosseae or C. etunicatum increased plant dry weight, decreased Se accumulation and selenocysteine proportion, enhanced root peroxidase activity, and alleviated oxidative stress in Se(VI)-treated plants. Inoculation also downregulated the expression of genes encoding Se transporters, assimilation enzymes, and cysteine-rich receptor-like kinases in Se(VI)-stressed plants, similar to plant-pathogen interaction and glutathione metabolism related genes. Conversely, genes encoding selenium-binding proteins and those related to phenylpropanoid biosynthesis were upregulated in inoculated plants under Se(VI) stress. Compared with Se(VI)-free plants, Se tolerance index, symbiotic feedback percentage on plant dry weight, and root colonization rate were all increased in inoculated plants under Se(VI) stress, corresponding to upregulated expression of 'key genes' in symbiosis. AMF inoculation increased bacterial diversity, decreased the relative abundances of selenobacteria related to plant Se absorption (e.g., Proteobacteria and Firmicutes), and improved bacterial network complexity in Se(VI)-stressed soils. We suggest that stress-mediated enhancement of mycorrhizal symbiosis contributed to plant Se(VI) tolerance, whereas AMF-mediated reshaping of soil bacterial community structure prevented excessive Se accumulation in maize.

Integrating plant and fungal quantitative genetics to improve the ecological and agricultural applications of mycorrhizal symbioses

Finding and targeting genes that quantitatively contribute to agricultural and ecological processes progresses food production and conservation efforts. Typically, quantitative genetic approaches link variants in a single organism's genome with a trait of interest. Recently, genome-to-genome mapping has found genome variants interacting between species to produce the result of a multiorganism (including multikingdom) interaction. These were plant and bacterial pathogen genome interactions; plant-fungal coquantitative genetics have not yet been applied. Plant-mycorrhizae symbioses exist across most biomes, for a majority of land plants, including crop plants, and manipulate many traits from single organisms to ecosystems for which knowing the genetic basis would be useful. The availability of Rhizophagus irregularis mycorrhizal isolates, with genomic information, makes dual-genome methods with beneficial mutualists accessible and imminent.

A roadmap to understanding diversity and function of coral reef-associated fungi

Tropical coral reefs are hotspots of marine productivity, owing to the association of reef-building corals with endosymbiotic algae and metabolically diverse bacterial communities. However, the functional importance of fungi, well-known for their contribution to shaping terrestrial ecosystems and global nutrient cycles, remains underexplored on coral reefs. We here conceptualize how fungal functional traits may have facilitated the spread, diversification, and ecological adaptation of marine fungi on coral reefs. We propose that functions of reef-associated fungi may be diverse and go beyond their hitherto described roles of pathogens and bioeroders, including but not limited to reef-scale biogeochemical cycles and the structuring of coral-associated and environmental microbiomes via chemical mediation. Recent technological and conceptual advances will allow the elucidation of the physiological, ecological, and chemical contributions of understudied marine fungi to coral holobiont and reef ecosystem functioning and health and may help provide an outlook for reef management actions.

Nutrition vs association: plant defenses are altered by arbuscular mycorrhizal fungi association not by nutritional provisioning alone

BACKGROUND

While it is known that arbuscular mycorrhizal fungi (AMF) can improve nutrient acquisition and herbivore resistance in crops, the mechanisms by which AMF influence plant defense remain unknown. Plants respond to herbivory with a cascade of gene expression and phytochemical biosynthesis. Given that the production of defensive phytochemicals requires nutrients, a commonly invoked hypothesis is that the improvement to plant defense when grown with AMF is simply due to an increased availability of nutrients. An alternative hypothesis is that the AMF effect on herbivory is due to changes in plant defense gene expression that are not simply due to nutrient availability. In this study, we tested whether changes in plant defenses are regulated by nutritional provisioning alone or the response of plant to AMF associations. Maize plants grown with or without AMF and with one of three fertilizer treatments (standard, 2 × nitrogen, or 2 × phosphorous) were infested with fall armyworm (Spodoptera frugiperda; FAW) for 72 h. We measured general plant characteristics (e.g. height, number of leaves), relative gene expression (rtPCR) of three defensive genes (lox3, mpi, and pr5), total plant N and P nutrient content, and change in FAW mass per plant.

RESULTS

We found that AMF drove the defense response of maize by increasing the expression of mpi and pr5. Furthermore, while AMF increased the total phosphorous content of maize it had no impact on maize nitrogen. Fertilization alone did not alter upregulation of any of the 3 induced defense genes tested, suggesting the mechanism through which AMF upregulate defenses is not solely via increased N or P plant nutrition.

CONCLUSION

This work supports that maize defense may be optimized by AMF associations alone, reducing the need for artificial inputs when managing FAW.

RocTest: A standardized method to assess the performance of root organ cultures in the propagation of arbuscular mycorrhizal fungi

Over the past three decades, root organ cultures (ROCs) have been the gold standard method for studying arbuscular mycorrhizal fungi (AMF) under in vitro conditions, and ROCs derived from various plant species have been used as hosts for AM monoxenic cultures. While there is compelling evidence that host identity can significantly modify AMF fitness, there is currently no standardized methodology to assess the performance of ROCs in the propagation of their fungal symbionts. We describe RocTest, a robust methodological approach that models the propagation of AMF in symbiosis with ROCs. The development of extraradical fungal structures and the pattern of sporulation are modeled using cumulative link mixed models and linear mixed models. We demonstrate functionality of RocTest by evaluating the performance of three species of ROCs (Daucus carota, Medicago truncatula, Nicotiana benthamiana) in the propagation of three species of AMF (Rhizophagus clarus, Rhizophagus irregularis, Glomus sp.). RocTest produces a simple graphical output to assess the performance of ROCs and shows that fungal propagation depends on the three-way interaction between ROC, AMF, and time. RocTest makes it possible to identify the best combination of host/AMF for fungal development and spore production, making it an important asset for germplasm collections and AMF research.

Integrating genomics and transcriptomics to identify candidate genes for subcutaneous fat deposition in beef cattle

Fat deposition is a complex economic trait regulated by polygenic genetic basis and environmental factors. Therefore, integrating multi-omics data to uncover its internal regulatory mechanism has attracted extensive attention. Here, we performed genomics and transcriptomics analysis to detect candidates affecting subcutaneous fat (SCF) deposition in beef cattle. The association of 770K SNPs with the backfat thickness captured nine significant SNPs within or near 11 genes. Additionally, 13 overlapping genes regarding fat deposition were determined via the analysis of differentially expressed genes and weighted gene co-expression network analysis (WGCNA). We then calculated the correlations of these genes with BFT and constructed their interaction network. Finally, seven biomarkers including ACACA, SCD, FASN, ACOX1, ELOVL5, HACD2, and HSD17B12 were screened. Notably, ACACA, identified by the integration of genomics and transcriptomics, was more likely to exert profound effects on SCF deposition. These findings provided novel insights into the regulation mechanism underlying bovine fat accumulation.

The Costs and Benefits of Plant-Arbuscular Mycorrhizal Fungal Interactions

The symbiotic interaction between plants and arbuscular mycorrhizal (AM) fungi is often perceived as beneficial for both partners, though a large ecological literature highlights the context dependency of this interaction. Changes in abiotic variables, such as nutrient availability, can drive the interaction along the mutualism-parasitism continuum with variable outcomes for plant growth and fitness. However, AM fungi can benefit plants in more ways than improved phosphorus nutrition and plant growth. For example, AM fungi can promote abiotic and biotic stress tolerance even when considered parasitic from a nutrient provision perspective. Other than being obligate biotrophs, very little is known about the benefits AM fungi gain from plants. In this review, we utilize both molecular biology and ecological approaches to expand our understanding of the plant-AM fungal interaction across disciplines.

The hidden world within plants: metatranscriptomics unveils the complexity of wood microbiomes

The importance of plants as complex entities influenced by genomes of the associated microorganisms is now seen as a new source of variability for a more sustainable agriculture, also in the light of ongoing climate change. For this reason, we investigated through metatranscriptomics whether the taxa profile and behaviour of microbial communities associated with the wood of 20-year-old grapevine plants are influenced by the health status of the host. We report for the first time a metatranscriptome from a complex tissue in a real environment, highlighting that this approach is able to define the microbial community better than referenced transcriptomic approaches. In parallel, the use of total RNA enabled the identification of bacterial taxa in healthy samples that, once isolated from the original wood tissue, displayed potential biocontrol activities against a wood-degrading fungal taxon. Furthermore, we revealed an unprecedented high number of new viral entities (~120 new viral species among 180 identified) associated with a single and limited environment and with potential impact on the whole holobiont. Taken together, our results suggest a complex multitrophic interaction in which the viral community also plays a crucial role in raising new ecological questions for the exploitation of microbial-assisted sustainable agriculture.

Harnessing an arbuscular mycorrhizal fungus to improve the adaptability of a facultative metallophytic poplar (Populus yunnanensis) to cadmium stress: Physiological and molecular responses

Populus yunnanensis Dode, a facultative metallophytic poplar, exhibits afforestation potential in barren mine tailing areas. However, the interactions and functional roles of arbuscular mycorrhizal fungus (AMF) in P. yunnanensis adaptability to heavy metal stress remain unclear. Physiological and molecular responses of P. yunnanensis plantlets to AMF (Funneliformis mosseae) under cadmium (Cd) stress (50 mg kg^-1) were investigated. Results showed attenuation of Cd phytotoxicity effects on cell organelles upon AMF inoculation, which also reduced the Cd concentration in the poplar leaves, stems, and roots. Under Cd stress, AMF-blocking of metal transporter (e.g., Ca²⁺ channel) activity occurred, decreasing root cell Cd influx by reducing H⁺ efflux. Bioaugmentation of rhizosphere sediments by AMF to stabilize metals with a decreasing DTPA-extractable Cd also occurred. The AMF inoculation promoted Cd conversion into inactive, less phytotoxic forms, and helped to maintain ion homeostasis and relieve nutritional ion (e.g., Ca, Mg) disorders caused by excessive Cd. Leaf enzyme and non-enzyme antioxidant systems were triggered. Root and leaf physiological response patterns differed. The AMF regulated the poplar functional genes, and nine metal-responsive gene clusters were identified. We suggest that AMF is a functional component of P. yunnanensis phenotype extension, contributing to strong adaptability to unfavorable mine tailings conditions.

Symbiotic responses of Lotus japonicus to two isogenic lines of a mycorrhizal fungus differing in the presence/absence of an endobacterium

As other arbuscular mycorrhizal fungi, Gigaspora margarita contains unculturable endobacteria in its cytoplasm. A cured fungal line has been obtained and showed it was capable of establishing a successful mycorrhizal colonization. However, previous OMICs and physiological analyses have demonstrated that the cured fungus is impaired in some functions during the pre-symbiotic phase, leading to a lower respiration activity, lower ATP, and antioxidant production. Here, by combining deep dual-mRNA sequencing and proteomics applied to Lotus japonicus roots colonized by the fungal line with bacteria (B+) and by the cured line (B-), we tested the hypothesis that L. japonicus (i) activates its symbiotic pathways irrespective of the presence or absence of the endobacterium, but (ii) perceives the two fungal lines as different physiological entities. Morphological observations confirmed the absence of clear endobacteria-dependent changes in the mycorrhizal phenotype of L. japonicus, while transcript and proteomic datasets revealed activation of the most important symbiotic pathways. They included the iconic nutrient transport and some less-investigated pathways, such as phenylpropanoid biosynthesis. However, significant differences between the mycorrhizal B+/B- plants emerged in the respiratory pathways and lipid biosynthesis. In both cases, the roots colonized by the cured line revealed a reduced capacity to activate genes involved in antioxidant metabolism, as well as the early biosynthetic steps of the symbiotic lipids, which are directed towards the fungus. Similar to its pre-symbiotic phase, the intraradical fungus revealed transcripts related to mitochondrial activity, which were downregulated in the cured line, as well as perturbation in lipid biosynthesis.

A historical perspective on mycorrhizal mutualism emphasizing arbuscular mycorrhizas and their emerging challenges

Arbuscular mycorrhiza, one of the oldest interactions on earth (~ 450 million years old) and a first-class partner for plants to colonize emerged land, is considered one of the most pervasive ecological relationships on the globe. Despite how important and old this interaction is, its discovery was very recent compared to the long story of land plant evolution. The story of the arbuscular mycorrhiza cannot be addressed apart from the history, controversies, and speculations about mycorrhiza in its broad sense. The chronicle of mycorrhizal research is marked by multiple key milestones such as the initial description of a "persistent epiderm and pellicular wall structure" by Hartig; the introduction of the "Symbiotismus" and "Mycorrhiza" concepts by Frank; the description of diverse root-fungal morphologies; the first description of arbuscules by Gallaud; Mosse's pivotal statement of the beneficial nature of the arbuscular mycorrhizal symbiosis; the impact of molecular tools on the taxonomy of mycorrhizal fungi as well as the development of in vitro root organ cultures for producing axenic arbuscular mycorrhizal fungi (AMF). An appreciation of the story - full of twists and turns - of the arbuscular mycorrhiza, going from the roots of mycorrhiza history, along with the discovery of different mycorrhiza types such as ectomycorrhiza, can improve research to help face our days' challenge of developing sustainable agriculture that integrates the arbuscular mycorrhiza and its ecosystem services.

Differential responses of 23 maize cultivar seedlings to an arbuscular mycorrhizal fungus when grown in a metal-polluted soil

Modern breeding efforts have been accelerating crop improvement and yielding numerous cultivars with distinct genetic traits; however, interactions between different cultivars and their root-associated arbuscular mycorrhizal fungi (AMF) are not clear. Herein, we selected the 22 most common commercial maize (Zea mays) varieties in China and an inbred line (B73) to study the differential responses of these 23 cultivars to mycorrhizal inoculation when grown in an arable soil polluted by multiple metals (Pb, Zn, and Cd). We found that the different cultivars exhibited significant variations in plant metal accumulation, ranging from strong metal exclusion (ZYY9) to strong metal accumulation (B73). Mycorrhizal colonization substantially altered metal uptake and repartitioning, while bioaugmenting the inherent characteristics of metal accumulation; for example, the AMF enhanced leaf accumulation of the metal-accumulator B73, and markedly reduced the root uptake of the metal-excluder ZYY9. However, such AMF-induced alterations were also substantially dependent on plant organs (roots and shoots) and metal species. We found that the extent of the AMF-induced leaf alterations was substantially greater than that of the root alterations. Similarly, the number of instances where the AMF significantly altered the Zn and Cd accumulation was far higher than the number of instances where Pb accumulation was significantly altered by AMF. In addition, the presence of AMF appeared to trigger the maize antioxidant systems, which may have alleviated the toxicity of excessive Cd, increased the leaf chlorophyll content, augmented the net photosynthetic rate, and promoted the growth of 17.39% of the maize cultivars. Our results suggest that a future crop breeding challenge is to produce cultivars for safe production or phytoremediation, thereby optimizing the combinations of crop cultivars and their root-associated AMF in slightly to moderately metal-polluted arable soils.

Generation of unequal nuclear genotype proportions in Rhizophagus irregularis progeny causes allelic imbalance in gene transcription

Arbuscular mycorrhizal fungi (AMF) form mutualisms with most plant species. The model AMF Rhizophagus irregularis is common in many ecosystems and naturally forms homokaryons and dikaryons. Quantitative variation in allele frequencies in clonally dikaryon offspring suggests they disproportionately inherit two distinct nuclear genotypes from their parent. This is interesting, because such progeny strongly and differentially affect plant growth. Neither the frequency and magnitude of this occurrence nor its effect on gene transcription are known. Using reduced representation genome sequencing, transcriptomics, and quantitative analysis tools, we show that progeny of homokaryons and dikaryons are qualitatively genetically identical to the parent. However, dikaryon progeny differ quantitatively due to unequal inheritance of nuclear genotypes. Allele frequencies of actively transcribed biallelic genes resembled the frequencies of the two nuclear genotypes. More biallelic genes showed transcription of both alleles than monoallelic transcription, but biallelic transcription was less likely with greater allelic divergence. Monoallelic transcription levels of biallelic genes were reduced compared with biallelic gene transcription, a finding consistent with genomic conflict. Given that genetic variation in R. irregularis is associated with plant growth, our results establish quantitative genetic variation as a future consideration when selecting AMF lines to improve plant production.

Comparative Fungal Community Analyses Using Metatranscriptomics and Internal Transcribed Spacer Amplicon Sequencing from Norway Spruce

The health, growth, and fitness of boreal forest trees are impacted and improved by their associated microbiomes. Microbial gene expression and functional activity can be assayed with RNA sequencing (RNA-Seq) data from host samples. In contrast, phylogenetic marker gene amplicon sequencing data are used to assess taxonomic composition and community structure of the microbiome. Few studies have considered how much of this structural and taxonomic information is included in transcriptomic data from matched samples. Here, we described fungal communities using both host-derived RNA-Seq and fungal ITS1 DNA amplicon sequencing to compare the outcomes between the methods. We used a panel of root and needle samples from the coniferous tree species Picea abies (Norway spruce) growing in untreated (nutrient-deficient) and nutrient-enriched plots at the Flakaliden forest research site in boreal northern Sweden. We show that the relationship between samples and alpha and beta diversity indicated by the fungal transcriptome is in agreement with that generated by the ITS data, while also identifying a lack of taxonomic overlap due to limitations imposed by current database coverage. Furthermore, we demonstrate how metatranscriptomics data additionally provide biologically informative functional insights. At the community level, there were changes in starch and sucrose metabolism, biosynthesis of amino acids, and pentose and glucuronate interconversions, while processing of organic macromolecules, including aromatic and heterocyclic compounds, was enriched in transcripts assigned to the genus Cortinarius.

IMPORTANCE A deeper understanding of microbial communities associated with plants is revealing their importance for plant health and productivity. RNA extracted from plant field samples represents the host and other organisms present. Typically, gene expression studies focus on the plant component or, in a limited number of studies, expression in one or more associated organisms. However, metatranscriptomic data are rarely used for taxonomic profiling, which is currently performed using amplicon approaches. We created an assembly-based, reproducible, and hardware-agnostic workflow to taxonomically and functionally annotate fungal RNA-Seq data obtained from Norway spruce roots, which we compared to matching ITS amplicon sequencing data. While we identified some limitations and caveats, we show that functional, taxonomic, and compositional insights can all be obtained from RNA-Seq data. These findings highlight the potential of metatranscriptomics to advance our understanding of interaction, response, and effect between host plants and their associated microbial communities.

Lineage-Specific Genes and Cryptic Sex: Parallels and Differences between Arbuscular Mycorrhizal Fungi and Fungal Pathogens

Arbuscular mycorrhizal fungi (AMF) live as obligate root symbionts on almost all land plants. They have long been regarded as ancient asexuals that have propagated clonally for millions of years. However, genomic studies in Rhizophagus irregularis and other AMF revealed many features indicative of sex. Surprisingly, comparative genomics of conspecific isolates of R. irregularis revealed an unexpected interstrain diversity, suggesting that AMF carry a high number of lineage-specific (LS) genes. Intriguingly, cryptic sex and LS genomic regions have previously been reported in a number of fungal pathogens of plants and humans. Here, we discuss these genomic similarities and highlight their potential relevance for AMF adaptation to the environment and for symbiotic functioning.

Coexistence of genetically different Rhizophagus irregularis isolates induces genes involved in a putative fungal mating response

Arbuscular mycorrhizal fungi (AMF) are of great ecological importance because of their effects on plant growth. Closely related genotypes of the same AMF species coexist in plant roots. However, almost nothing is known about the molecular interactions occurring during such coexistence. We compared in planta AMF gene transcription in single and coinoculation treatments with two genetically different isolates of Rhizophagus irregularis in symbiosis independently on three genetically different cassava genotypes. Remarkably few genes were specifically upregulated when the two fungi coexisted. Strikingly, almost all of the genes with an identifiable putative function were known to be involved in mating in other fungal species. Several genes were consistent across host plant genotypes but more upregulated genes involved in putative mating were observed in host genotype (COL2215) compared with the two other host genotypes. The AMF genes that we observed to be specifically upregulated during coexistence were either involved in the mating pheromone response, in meiosis, sexual sporulation or were homologs of MAT-locus genes known in other fungal species. We did not observe the upregulation of the expected homeodomain genes contained in a putative AMF MAT-locus, but observed upregulation of HMG-box genes similar to those known to be involved in mating in Mucoromycotina species. Finally, we demonstrated that coexistence between the two fungal genotypes in the coinoculation treatments explained the number of putative mating response genes activated in the different plant host genotypes. This study demonstrates experimentally the activation of genes involved in a putative mating response and represents an important step towards the understanding of coexistence and sexual reproduction in these important plant symbionts.

Spatial and morphological reorganization of endosymbiosis during metamorphosis accommodates adult metabolic requirements in a weevil

Bacterial intracellular symbiosis (endosymbiosis) is widespread in nature and impacts many biological processes. In holometabolous symbiotic insects, metamorphosis entails a complete and abrupt internal reorganization that creates a constraint for endosymbiont transmission from larvae to adults. To assess how endosymbiosis copes-and potentially evolves-throughout this major host-tissue reorganization, we used the association between the cereal weevil Sitophilus oryzae and the bacterium Sodalis pierantonius as a model system. S. pierantonius are contained inside specialized host cells, the bacteriocytes, that group into an organ, the bacteriome. Cereal weevils require metabolic inputs from their endosymbiont, particularly during adult cuticle synthesis, when endosymbiont load increases dramatically. By combining dual RNA-sequencing analyses and cell imaging, we show that the larval bacteriome dissociates at the onset of metamorphosis and releases bacteriocytes that undergo endosymbiosis-dependent transcriptomic changes affecting cell motility, cell adhesion, and cytoskeleton organization. Remarkably, bacteriocytes turn into spindle cells and migrate along the midgut epithelium, thereby conveying endosymbionts to midgut sites where future mesenteric caeca will develop. Concomitantly, endosymbiont genes encoding a type III secretion system and a flagellum apparatus are transiently up-regulated while endosymbionts infect putative stem cells and enter their nuclei. Infected cells then turn into new differentiated bacteriocytes and form multiple new bacteriomes in adults. These findings show that endosymbiosis reorganization in a holometabolous insect relies on a synchronized host-symbiont molecular and cellular "choreography" and illustrates an adaptive feature that promotes bacteriome multiplication to match increased metabolic requirements in emerging adults.

The combined analysis as the best strategy for Dual RNA-Seq mapping

In Dual RNA-Seq experiments the simultaneous extraction of RNA and analysis of gene expression data from both interacting organisms could be a challenge. One alternative is separating the reads during in silico data analysis. There are two main mapping methods used: sequential and combined. Here we present a combined approach in which the libraries were aligned to a concatenated genome to sort the reads before mapping them to the respective annotated genomes. A comparison of this method with the sequential analysis was performed. Two RNA-Seq libraries available in public databases consisting of a eukaryotic (Zea mays) and a prokaryotic (Herbaspirillum seropediceae) organisms were mixed to simulate a Dual RNA-Seq experiment. Libraries from real Dual RNA-Seq experiments were also used. The sequential analysis consistently attributed more reads to the first reference genome used in the analysis (due to cross-mapping) than the combined approach. More importantly, the combined analysis resulted in lower numbers of cross-mapped reads. Our results highlight the necessity of combining the reference genomes to sort reads previously to the counting step to avoid losing information in Dual RNA-Seq experiments. Since most studies first map the RNA-Seq libraries to the eukaryotic genome, much prokaryotic information has probably been lost.

Bioinformatic Methods for the Analysis of High-Throughput RNA Sequencing in Arbuscular Mycorrhizal Fungi

RNA-seq is a powerful method for transcriptome profiling that allows the detection of total RNA present in a single cell, tissues, or organs. mRNA-seq is focused on protein-coding RNAs, and results in large datasets of reads, or portion of sequenced mRNA that can be assembled back to the original transcripts to reconstruct a virtual gene catalog. Studies on the biology of arbuscular mycorrhizal fungi (AMF) often took great advantage of mRNA-seq, and several attempts to decipher their coding potential relied on de novo transcriptome assembly. As the transcriptional profile of an organism is modulated depending on cell types, and in response to specific biological conditions, mRNA-seq is an attractive approach to study the physiology of AMF, which are axenically unculturable and genetically intractable. mRNA-seq analyses require bioinformatic workflows to manipulate the huge amount of raw data generated by the sequencing run, with several crucial steps (e.g., library trimming, reads mapping, normalization, and differential expression calculation) which can strongly affect the final results. Here, we propose a standard workflow for de novo transcriptome assembly and differential expression calculation for AMF, which considers the most common technical issues of working in the absence of reference sequences and with mixed biological samples.

Gene expression in Rhizoglomus irregulare at two different time points of mycorrhiza establishment in Helianthus annuus roots, as revealed by RNA-seq analysis

Arbuscular mycorrhizal fungi (AMF) play a fundamental role in plant growth and nutrition in natural and agricultural ecosystems. Despite the importance of such symbionts, the different developmental changes occurring during the AMF life cycle have not been fully elucidated at the molecular level. Here, the RNA-seq approach was used to investigate Rhizoglomus irregulare specific and common transcripts at two different time points of mycorrhizal establishment in Helianthus annuus in vivo. Four days after inoculation, transcripts related to cellular remodeling (actin and tubulin), cellular signaling (calmodulin, serine/threonine protein kinase, 14-3-3 protein, and calcium transporting ATPase), lipid metabolism (fatty acid desaturation, steroid hormone, and glycerophospholipid biosynthesis), and biosynthetic processes were detected. In addition to such transcripts, 16 days after inoculation, expressed genes linked to binding and catalytic activities; ion (K⁺, Ca²⁺, Fe²⁺, Zn²⁺, Mn²⁺, Pi, ammonia), sugar, and lipid transport; and those involved in vacuolar polyphosphate accumulation were found. Knowledge of transcriptomic changes required for symbiosis establishment and performance is of great importance to understand the functional role of AMF symbionts in food crop nutrition and health, and in plant diversity in natural ecosystems.

Next-generation technologies for studying host-pathogen interactions: a focus on dual transcriptomics, CRISPR/Cas9 screening and organs-on-chips

Pathogens constantly interact with their hosts and the environment, and therefore have evolved unique virulence mechanisms to target and breach host defense barriers and manipulate host immune response to establish an infection. Advances in technologies that allow genome mining, gene editing such as CRISPR/Cas9, genomic, epigenomic and transcriptomic studies such as dual RNA-seq, coupled with bioinformatics, have accelerated the field of host-pathogen interactions within a broad range of infection models. Underpinning of the molecular changes that accompany invasion of eukaryotic cells with pathogenic microorganisms at the intersection of host, pathogen and their local environment has provided a better understanding of infectious disease mechanisms and antimicrobial strategies. The recent evolution of physiologically relevant three-dimensional (3-D) tissue/organ models and microfluidic organ-on-chip devices also provided a window to a more predictive framework of infectious disease processes. These approaches combined hold the potential to highly impact discovery of novel drug targets and vaccine candidates of the future. Here, we review three of the available and emerging technologies-dual RNA-seq, CRISPR/Cas9 screening and organs-on-chips, applicable to the high throughput study and deciphering of interaction networks between pathogens and their hosts that are critical for the development of novel therapeutics.

Tomato RNA-seq Data Mining Reveals the Taxonomic and Functional Diversity of Root-Associated Microbiota

Next-generation approaches have enabled researchers to deeply study the plant microbiota and to reveal how microbiota associated with plant roots has key effects on plant nutrition, disease resistance, and plant development. Although early "omics" experiments focused mainly on the species composition of microbial communities, new "meta-omics" approaches such as meta-transcriptomics provide hints about the functions of the microbes when interacting with their plant host. Here, we used an RNA-seq dataset previously generated for tomato (Solanum lycopersicum) plants growing on different native soils to test the hypothesis that host-targeted transcriptomics can detect the taxonomic and functional diversity of root microbiota. Even though the sequencing throughput for the microbial populations was limited, we were able to reconstruct the microbial communities and obtain an overview of their functional diversity. Comparisons of the host transcriptome and the meta-transcriptome suggested that the composition and the metabolic activities of the microbiota shape plant responses at the molecular level. Despite the limitations, mining available next-generation sequencing datasets can provide unexpected results and potential benefits for microbiota research.

Transcriptome Analysis of Wheat Roots Reveals a Differential Regulation of Stress Responses Related to Arbuscular Mycorrhizal Fungi and Soil Disturbance

Symbioses with soil microorganisms are central in shaping the diversity and productivity of land plants and provide protection against a diversity of stresses, including metal toxicity. Arbuscular mycorrhizal fungi (AMF) can form extensive extraradical mycelial networks (ERM), which are very efficient in colonizing a new host. We quantified the responses of transcriptomes of wheat and one AMF partner, Rhizoglomus irregulare, to soil disturbance (Undisturbed vs. Disturbed) and to two different preceding mycotrophic species (Ornithopus compressus and Lolium rigidum). Soil disturbance and preceding plant species engender different AMF communities in wheat roots, resulting in a differential tolerance to soil manganese (Mn) toxicity. Soil disturbance negatively impacted wheat growth under manganese toxicity, probably due to the disruption of the ERM, and activated a large number of stress and starvation-related genes. The O. compressus treatment, which induces a greater Mn protection in wheat than L. rigidum, activated processes related to cellular division and growth, and very few related to stress. The L. rigidum treatment mostly induced genes that were related to oxidative stress, disease protection, and metal ion binding. R. irregulare cell division and molecular exchange between nucleus and cytoplasm were increased by O. compressus. These findings are highly relevant for sustainable agricultural systems, when considering a fit-for-purpose symbiosis.

In vitro Propagation of Arbuscular Mycorrhizal Fungi May Drive Fungal Evolution

Transformed root cultures (TRC) are used to mass produce arbuscular mycorrhizal (AM) fungal propagules in vitro. These propagules are then used in research, agriculture, and ecological restoration. There are many examples from other microbial systems that long-term in vitro propagation leads to domesticated strains that differ genetically and functionally. Here, we discuss potential consequences of in TRC propagation on AM fungal traits, and how this may affect their functionality. We examine weather domestication of AM fungi has already happened and finally, we explore whether it is possible to overcome TRC-induced domestication.

Biocontrol of Cereal Crop Diseases Using Streptomycetes

A growing world population and an increasing demand for greater food production requires that crop losses caused by pests and diseases are dramatically reduced. Concurrently, sustainability targets mean that alternatives to chemical pesticides are becoming increasingly desirable. Bacteria in the plant root microbiome can protect their plant host against pests and pathogenic infection. In particular, Streptomyces species are well-known to produce a range of secondary metabolites that can inhibit the growth of phytopathogens. Streptomyces are abundant in soils and are also enriched in the root microbiomes of many different plant species, including those grown as economically and nutritionally valuable cereal crops. In this review we discuss the potential of Streptomyces to protect against some of the most damaging cereal crop diseases, particularly those caused by fungal pathogens. We also explore factors that may improve the efficacy of these strains as biocontrol agents in situ, as well as the possibility of exploiting plant mechanisms, such as root exudation, that enable the recruitment of microbial species from the soil to the root microbiome. We argue that a greater understanding of these mechanisms may enable the development of protective plant root microbiomes with a greater abundance of beneficial bacteria, such as Streptomyces species.