JGI Plant Gene Atlas: an updateable transcriptome resource to improve functional gene descriptions across the plant kingdom.

Gene functional descriptions offer a crucial line of evidence for candidate genes underlying trait variation. Conversely, plant responses to environmental cues represent important resources to decipher gene function and subsequently provide molecular targets for plant improvement through gene editing. However, biological roles of large proportions of genes across the plant phylogeny are poorly annotated. Here we describe the Joint Genome Institute (JGI) Plant Gene Atlas, an updateable data resource consisting of transcript abundance assays spanning 18 diverse species. To integrate across these diverse genotypes, we analyzed expression profiles, built gene clusters that exhibited tissue/condition specific expression, and tested for transcriptional response to environmental queues. We discovered extensive phylogenetically constrained and condition-specific expression profiles for genes without any previously documented functional annotation. Such conserved expression patterns and tightly co-expressed gene clusters let us assign expression derived additional biological information to 64 495 genes with otherwise unknown functions. The ever-expanding Gene Atlas resource is available at JGI Plant Gene Atlas (https://plantgeneatlas.jgi.doe.gov) and Phytozome (https://phytozome.jgi.doe.gov/), providing bulk access to data and user-specified queries of gene sets. Combined, these web interfaces let users access differentially expressed genes, track orthologs across the Gene Atlas plants, graphically represent co-expressed genes, and visualize gene ontology and pathway enrichments.

Predicted structural proteome of Sphagnum divinum and proteome-scale annotation.

MOTIVATION: Sphagnum-dominated peatlands store a substantial amount of terrestrial carbon. The genus is undersampled and under-studied. No experimental crystal structure from any Sphagnum species exists in the Protein Data Bank and fewer than 200 Sphagnum-related genes have structural models available in the AlphaFold Protein Structure Database. Tools and resources are needed to help bridge these gaps, and to enable the analysis of other structural proteomes now made possible by accurate structure prediction. RESULTS: We present the predicted structural proteome (25 134 primary transcripts) of Sphagnum divinum computed using AlphaFold, structural alignment results of all high-confidence models against an annotated nonredundant crystallographic database of over 90,000 structures, a structure-based classification of putative Enzyme Commission (EC) numbers across this proteome, and the computational method to perform this proteome-scale structure-based annotation. AVAILABILITY AND IMPLEMENTATION: All data and code are available in public repositories, detailed at https://github.com/BSDExabio/SAFA. The structural models of the S. divinum proteome have been deposited in the ModelArchive repository at https://modelarchive.org/doi/10.5452/ma-ornl-sphdiv.

Temperature and CO2 interactively drive shifts in the compositional and functional structure of peatland protist communities.

Microbes affect the global carbon cycle that influences climate change and are in turn influenced by environmental change. Here, we use data from a long-term whole-ecosystem warming experiment at a boreal peatland to answer how temperature and CO2 jointly influence communities of abundant, diverse, yet poorly understood, non-fungi microbial Eukaryotes (protists). These microbes influence ecosystem function directly through photosynthesis and respiration, and indirectly, through predation on decomposers (bacteria and fungi). Using a combination of high-throughput fluid imaging and 18S amplicon sequencing, we report large climate-induced, community-wide shifts in the community functional composition of these microbes (size, shape, and metabolism) that could alter overall function in peatlands. Importantly, we demonstrate a taxonomic convergence but a functional divergence in response to warming and elevated CO2 with most environmental responses being contingent on organismal size: warming effects on functional composition are reversed by elevated CO2 and amplified in larger microbes but not smaller ones. These findings show how the interactive effects of warming and rising CO2 levels could alter the structure and function of peatland microbial food webs-a fragile ecosystem that stores upwards of 25% of all terrestrial carbon and is increasingly threatened by human exploitation.

A research agenda for nonvascular photoautotrophs under climate change.

Nonvascular photoautotrophs (NVP), including bryophytes, lichens, terrestrial algae, and cyanobacteria, are increasingly recognized as being essential to ecosystem functioning in many regions of the world. Current research suggests that climate change may pose a substantial threat to NVP, but the extent to which this will affect the associated ecosystem functions and services is highly uncertain. Here, we propose a research agenda to address this urgent question, focusing on physiological and ecological processes that link NVP to ecosystem functions while also taking into account the substantial taxonomic diversity across multiple ecosystem types. Accordingly, we developed a new categorization scheme, based on microclimatic gradients, which simplifies the high physiological and morphological diversity of NVP and world-wide distribution with respect to several broad habitat types. We found that habitat-specific ecosystem functions of NVP will likely be substantially affected by climate change, and more quantitative process understanding is required on: (1) potential for acclimation; (2) response to elevated CO2 ; (3) role of the microbiome; and (4) feedback to (micro)climate. We suggest an integrative approach of innovative, multimethod laboratory and field experiments and ecophysiological modelling, for which sustained scientific collaboration on NVP research will be essential.

Climate drivers alter nitrogen availability in surface peat and decouple N2 fixation from CH4 oxidation in the Sphagnum moss microbiome.

Peat mosses (Sphagnum spp.) are keystone species in boreal peatlands, where they dominate net primary productivity and facilitate the accumulation of carbon in thick peat deposits. Sphagnum mosses harbor a diverse assemblage of microbial partners, including N2 -fixing (diazotrophic) and CH4 -oxidizing (methanotrophic) taxa that support ecosystem function by regulating transformations of carbon and nitrogen. Here, we investigate the response of the Sphagnum phytobiome (plant + constituent microbiome + environment) to a gradient of experimental warming (+0°C to +9°C) and elevated CO2 (+500 ppm) in an ombrotrophic peatland in northern Minnesota (USA). By tracking changes in carbon (CH4 , CO2 ) and nitrogen (NH4 -N) cycling from the belowground environment up to Sphagnum and its associated microbiome, we identified a series of cascading impacts to the Sphagnum phytobiome triggered by warming and elevated CO2 . Under ambient CO2 , warming increased plant-available NH4 -N in surface peat, excess N accumulated in Sphagnum tissue, and N2 fixation activity decreased. Elevated CO2 offset the effects of warming, disrupting the accumulation of N in peat and Sphagnum tissue. Methane concentrations in porewater increased with warming irrespective of CO2 treatment, resulting in a ~10× rise in methanotrophic activity within Sphagnum from the +9°C enclosures. Warming's divergent impacts on diazotrophy and methanotrophy caused these processes to become decoupled at warmer temperatures, as evidenced by declining rates of methane-induced N2 fixation and significant losses of keystone microbial taxa. In addition to changes in the Sphagnum microbiome, we observed ~94% mortality of Sphagnum between the +0°C and +9°C treatments, possibly due to the interactive effects of warming on N-availability and competition from vascular plant species. Collectively, these results highlight the vulnerability of the Sphagnum phytobiome to rising temperatures and atmospheric CO2 concentrations, with significant implications for carbon and nitrogen cycling in boreal peatlands.

Divergent selection and climate adaptation fuel genomic differentiation between sister species of Sphagnum (peat moss).

BACKGROUND AND AIMS: New plant species can evolve through the reinforcement of reproductive isolation via local adaptation along habitat gradients. Peat mosses (Sphagnaceae) are an emerging model system for the study of evolutionary genomics and have well-documented niche differentiation among species. Recent molecular studies have demonstrated that the globally distributed species Sphagnum magellanicum is a complex of morphologically cryptic lineages that are phylogenetically and ecologically distinct. Here, we describe the architecture of genomic differentiation between two sister species in this complex known from eastern North America: the northern S. diabolicum and the largely southern S. magniae. METHODS: We sampled plant populations from across a latitudinal gradient in eastern North America and performed whole genome and restriction-site associated DNA sequencing. These sequencing data were then analyzed computationally. KEY RESULTS: Using sliding-window population genetic analyses we find that differentiation is concentrated within 'islands' of the genome spanning up to 400 kb that are characterized by elevated genetic divergence, suppressed recombination, reduced nucleotide diversity and increased rates of non-synonymous substitution. Sequence variants that are significantly associated with genetic structure and bioclimatic variables occur within genes that have functional enrichment for biological processes including abiotic stress response, photoperiodism and hormone-mediated signalling. Demographic modelling demonstrates that these two species diverged no more than 225 000 generations ago with secondary contact occurring where their ranges overlap. CONCLUSIONS: We suggest that this heterogeneity of genomic differentiation is a result of linked selection and reflects the role of local adaptation to contrasting climatic zones in driving speciation. This research provides insight into the process of speciation in a group of ecologically important plants and strengthens our predictive understanding of how plant populations will respond as Earth's climate rapidly changes.

Extensive Genome-Wide Phylogenetic Discordance Is Due to Incomplete Lineage Sorting and Not Ongoing Introgression in a Rapidly Radiated Bryophyte Genus.

The relative importance of introgression for diversification has long been a highly disputed topic in speciation research and remains an open question despite the great attention it has received over the past decade. Gene flow leaves traces in the genome similar to those created by incomplete lineage sorting (ILS), and identification and quantification of gene flow in the presence of ILS is challenging and requires knowledge about the true phylogenetic relationship among the species. We use whole nuclear, plastid, and organellar genomes from 12 species in the rapidly radiated, ecologically diverse, actively hybridizing genus of peatmoss (Sphagnum) to reconstruct the species phylogeny and quantify introgression using a suite of phylogenomic methods. We found extensive phylogenetic discordance among nuclear and organellar phylogenies, as well as across the nuclear genome and the nodes in the species tree, best explained by extensive ILS following the rapid radiation of the genus rather than by postspeciation introgression. Our analyses support the idea of ancient introgression among the ancestral lineages followed by ILS, whereas recent gene flow among the species is highly restricted despite widespread interspecific hybridization known in the group. Our results contribute to phylogenomic understanding of how speciation proceeds in rapidly radiated, actively hybridizing species groups, and demonstrate that employing a combination of diverse phylogenomic methods can facilitate untangling complex phylogenetic patterns created by ILS and introgression.

amPEPpy 1.0: a portable and accurate antimicrobial peptide prediction tool.

SUMMARY: Antimicrobial peptides (AMPs) are promising alternative antimicrobial agents. Currently, however, portable, user-friendly and efficient methods for predicting AMP sequences from genome-scale data are not readily available. Here we present amPEPpy, an open-source, multi-threaded command-line application for predicting AMP sequences using a random forest classifier. AVAILABILITY AND IMPLEMENTATION: amPEPpy is implemented in Python 3 and is freely available through GitHub (https://github.com/tlawrence3/amPEPpy). SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

Plant carbohydrate storage: intra- and inter-specific trade-offs reveal a major life history trait.

Trade-offs among carbon sinks constrain how trees physiologically, ecologically, and evolutionarily respond to their environments. These trade-offs typically fall along a productive growth to conservative, bet-hedging continuum. How nonstructural carbohydrates (NSCs) stored in living tree cells (known as carbon stores) fit in this trade-off framework is not well understood. We examined relationships between growth and storage using both within species genetic variation from a common garden, and across species phenotypic variation from a global database. We demonstrate that storage is actively accumulated, as part of a conservative, bet-hedging life history strategy. Storage accumulates at the expense of growth both within and across species. Within the species Populus trichocarpa, genetic trade-offs show that for each additional unit of wood area growth (in cm2 yr-1 ) that genotypes invest in, they lose 1.2 to 1.7 units (mg g-1 NSC) of storage. Across species, for each additional unit of area growth (in cm2 yr-1 ), trees, on average, reduce their storage by 9.5% in stems and 10.4% in roots. Our findings impact our understanding of basic plant biology, fit storage into a widely used growth-survival trade-off spectrum describing life history strategy, and challenges the assumptions of passive storage made in ecosystem models today.

Phylogenomic structure and speciation in an emerging model: the Sphagnum magellanicum complex (Bryophyta).

Sphagnum magellanicum is one of two Sphagnum species for which a reference-quality genome exists to facilitate research in ecological genomics. Phylogenetic and comparative genomic analyses were conducted based on resequencing data from 48 samples and RADseq analyses based on 187 samples. We report herein that there are four clades/species within the S. magellanicum complex in eastern North America and that the reference genome belongs to Sphagnum divinum. The species exhibit tens of thousands (RADseq) to millions (resequencing) of fixed nucleotide differences. Two species, however, referred to informally as S. diabolicum and S. magni because they have not been formally described, are differentiated by only 100 (RADseq) to 1000 (resequencing) of differences. Introgression among species in the complex is demonstrated using D-statistics and f4 ratios. One ecologically important functional trait, tissue decomposability, which underlies peat (carbon) accumulation, does not differ between segregates in the S. magellanicum complex, although previous research showed that many closely related Sphagnum species have evolved differences in decomposability/carbon sequestration. Phylogenetic resolution and more accurate species delimitation in the S. magellanicum complex substantially increase the value of this group for studying the early evolutionary stages of climate adaptation and ecological evolution more broadly.

Habitat-adapted microbial communities mediate Sphagnum peatmoss resilience to warming.

Sphagnum peatmosses are fundamental members of peatland ecosystems, where they contribute to the uptake and long-term storage of atmospheric carbon. Warming threatens Sphagnum mosses and is known to alter the composition of their associated microbiome. Here, we use a microbiome transfer approach to test if microbiome thermal origin influences host plant thermotolerance. We leveraged an experimental whole-ecosystem warming study to collect field-grown Sphagnum, mechanically separate the associated microbiome and then transfer onto germ-free laboratory Sphagnum for temperature experiments. Host and microbiome dynamics were assessed with growth analysis, Chla fluorescence imaging, metagenomics, metatranscriptomics and 16S rDNA profiling. Microbiomes originating from warming field conditions imparted enhanced thermotolerance and growth recovery at elevated temperatures. Metagenome and metatranscriptome analyses revealed that warming altered microbial community structure in a manner that induced the plant heat shock response, especially the HSP70 family and jasmonic acid production. The heat shock response was induced even without warming treatment in the laboratory, suggesting that the warm-microbiome isolated from the field provided the host plant with thermal preconditioning. Our results demonstrate that microbes, which respond rapidly to temperature alterations, can play key roles in host plant growth response to rapidly changing environments.

The physiological acclimation and growth response of Populus trichocarpa to warming.

Plant metabolic acclimation to thermal stress remains underrepresented in current global climate models. Gaps exist in our understanding of how metabolic processes (i.e., photosynthesis, respiration) acclimate over time and how aboveground versus belowground acclimation differs. We measured the thermal acclimation of Populus trichocarpa, comparing aboveground versus belowground physiology over time. Ninety genetically identical ramets were propagated in mesocosms that separated root and microbial components. After establishment at 25°C for 6 weeks, 60 clones were warmed +4 or +8°C and monitored for 10 weeks, measuring photosynthesis (A), leaf respiration (R), soil respiration (Rs ), root plus soil respiration (Rs+r ), and root respiration (Rr ). We observed thermal acclimation in both A and R, with rates initially increasing, then declining as the thermal photosynthetic optimum (Topt ) and the temperature-sensitivity (Q10 ) of respiration adjusted to warmer conditions. Photosynthetic acclimation was constructive, based on an increase in both Topt and peak A. Belowground, Rs+r decreased linearly with warming, while Rs rates declined abruptly, then remained constant with additional warming. Plant biomass was greatest at +4°C, with 30% allocated belowground. Rates of mass-based Rr were similar among treatments; however, root nitrogen declined at +8°C leading to less mass nitrogen-based Rr in that treatment. The Q10 -temperature relationship of Rr was affected by warming, leading to differing values among treatments. Aboveground acclimation exceeded belowground acclimation, and plant nitrogen-use mediated the acclimatory response. Results suggest that moderate climate warming (+4°C) may lead to acclimation and increased plant biomass production but increases in production could be limited with severe warming (+8°C).

Phylogenomics reveals convergent evolution of red-violet coloration in land plants and the origins of the anthocyanin biosynthetic pathway.

The flavonoids, one of the largest classes of plant secondary metabolites, are found in lineages that span the land plant phylogeny and play important roles in stress responses and as pigments. Perhaps the most well-studied flavonoids are the anthocyanins that have human health benefits and help plants attract pollinators, regulate hormone production, and confer resistance to abiotic and biotic stresses. The canonical biochemical pathway responsible for the production of these pigments is well-characterized for flowering plants yet its conservation across deep divergences in land plants remains debated and poorly understood. Many early land plants such as mosses, liverworts, and ferns produce flavonoid pigments, but their biosynthetic origins and homologies to the anthocyanin pathway remain uncertain. We conducted phylogenetic analyses using full genome sequences representing nearly all major green plant lineages to reconstruct the evolutionary history of the anthocyanin biosynthetic pathway then test the hypothesis that genes in this pathway are present in early land plants. We found that the entire pathway was not intact until the most recent common ancestor of seed plants and that orthologs of many downstream enzymes are absent from seedless plants including mosses, liverworts, and ferns. Our results also highlight the utility of phylogenetic inference, as compared to pairwise sequence similarity, in orthology assessment within large gene families that have complex duplication-loss histories. We suggest that the production of red-violet flavonoid pigments widespread in seedless plants, including the 3-deoxyanthocyanins, requires the activity of novel, as-yet discovered enzymes, and represents convergent evolution of red-violet coloration across land plants.

Tradeoffs between phenotypic plasticity and local adaptation influence the ecophysiology of the moss, Sphagnum magellanicum.

Bryophytes are a diverse plant group and are functionally different from vascular plants. Yet, their peculiarities are rarely considered in the theoretical frameworks for plants. Currently, we lack information about the magnitude and the importance of intraspecific variability in the ecophysiology of bryophytes and how these might translate to local adaptation-a prerequisite for adaptive evolution. Capitalizing on two ecologically distinct (hummock and hollow) phenotypes of Sphagnum magellanicum, we explored the magnitude and pattern of intraspecific variability in this species and asked whether the environmental-mediated changes in shoot and physiological traits are due to phenotypic plasticity or local adaptation. Size, pigmentation, and habitat type that distinguished the phenotypes in the field did not influence the trait responses under a transplant and factorial experiment. In addition, the magnitude and pattern of trait variability (e.g., branch, stem and capitulum mass) changed with the treatments, which suggest that trait responses were due largely to phenotypic plasticity. The trait responses also suggest that the ecophysiological needs for mosses to grow in clumps, where they maintain a uniform growth may have an overriding effect over the potential for a fixed adaptive response to environmental heterogeneity, which would constrain local adaptation. We conclude that extending the trait-based framework to mosses or making comparisons between mosses and vascular plants under any theoretical framework would only be meaningful to the extent that growth form and dispersal strategies are considered.

Correction to: Diel rewiring and positive selection of ancient plant proteins enabled evolution of CAM photosynthesis in Agave.

Following publication of the original article.

Experimental warming alters the community composition, diversity, and N2 fixation activity of peat moss (Sphagnum fallax) microbiomes.

Sphagnum-dominated peatlands comprise a globally important pool of soil carbon (C) and are vulnerable to climate change. While peat mosses of the genus Sphagnum are known to harbor diverse microbial communities that mediate C and nitrogen (N) cycling in peatlands, the effects of climate change on Sphagnum microbiome composition and functioning are largely unknown. We investigated the impacts of experimental whole-ecosystem warming on the Sphagnum moss microbiome, focusing on N2 fixing microorganisms (diazotrophs). To characterize the microbiome response to warming, we performed next-generation sequencing of small subunit (SSU) rRNA and nitrogenase (nifH) gene amplicons and quantified rates of N2 fixation activity in Sphagnum fallax individuals sampled from experimental enclosures over 2 years in a northern Minnesota, USA bog. The taxonomic diversity of overall microbial communities and diazotroph communities, as well as N2 fixation rates, decreased with warming (p < 0.05). Following warming, diazotrophs shifted from a mixed community of Nostocales (Cyanobacteria) and Rhizobiales (Alphaproteobacteria) to predominance of Nostocales. Microbiome community composition differed between years, with some diazotroph populations persisting while others declined in relative abundance in warmed plots in the second year. Our results demonstrate that warming substantially alters the community composition, diversity, and N2 fixation activity of peat moss microbiomes, which may ultimately impact host fitness, ecosystem productivity, and C storage potential in peatlands.

Relatively rare root endophytic bacteria drive plant resource allocation patterns and tissue nutrient concentration in unpredictable ways.

PREMISE: Plant endophytic bacterial strains can influence plant traits such as leaf area and root length. Yet, the influence of more complex bacterial communities in regulating overall plant phenotype is less explored. Here, in two complementary experiments, we tested whether we can predict plant phenotype response to changes in microbial community composition. METHODS: In the first study, we inoculated a single genotype of Populus deltoides with individual root endophytic bacteria and measured plant phenotype. Next, data from this single inoculation were used to predict phenotypic traits after mixed three-strain community inoculations, which we tested in the second experiment. RESULTS: By itself, each bacterial endophyte significantly but weakly altered plant phenotype relative to noninoculated plants. In a mixture, bacterial strain Burkholderia BT03, constituted at least 98% of community relative abundance. Yet, plant resource allocation and tissue nutrient concentrations were disproportionately influenced by Pseudomonas sp. GM17, GM30, and GM41. We found a 10% increase in leaf mass fraction and an 11% decrease in root mass fraction when replacing Pseudomonas GM17 with GM41 in communities containing both Pseudomonas GM30 and Burkholderia BT03. CONCLUSIONS: Our results indicate that interactions among endophytic bacteria may drive plant phenotype over the contribution of each strain individually. Additionally, we have shown that low-abundance strains contribute to plant phenotype challenging the assumption that the dominant strains will drive plant function.

Microbial-type terpene synthase genes occur widely in nonseed land plants, but not in seed plants.

The vast abundance of terpene natural products in nature is due to enzymes known as terpene synthases (TPSs) that convert acyclic prenyl diphosphate precursors into a multitude of cyclic and acyclic carbon skeletons. Yet the evolution of TPSs is not well understood at higher levels of classification. Microbial TPSs from bacteria and fungi are only distantly related to typical plant TPSs, whereas genes similar to microbial TPS genes have been recently identified in the lycophyte Selaginella moellendorffii The goal of this study was to investigate the distribution, evolution, and biochemical functions of microbial terpene synthase-like (MTPSL) genes in other plants. By analyzing the transcriptomes of 1,103 plant species ranging from green algae to flowering plants, putative MTPSL genes were identified predominantly from nonseed plants, including liverworts, mosses, hornworts, lycophytes, and monilophytes. Directed searching for MTPSL genes in the sequenced genomes of a wide range of seed plants confirmed their general absence in this group. Among themselves, MTPSL proteins from nonseed plants form four major groups, with two of these more closely related to bacterial TPSs and the other two to fungal TPSs. Two of the four groups contain a canonical aspartate-rich "DDxxD" motif. The third group has a "DDxxxD" motif, and the fourth group has only the first two "DD" conserved in this motif. Upon heterologous expression, representative members from each of the four groups displayed diverse catalytic functions as monoterpene and sesquiterpene synthases, suggesting these are important for terpene formation in nonseed plants.

Diel rewiring and positive selection of ancient plant proteins enabled evolution of CAM photosynthesis in Agave.

BACKGROUND: Crassulacean acid metabolism (CAM) enhances plant water-use efficiency through an inverse day/night pattern of stomatal closure/opening that facilitates nocturnal CO2 uptake. CAM has evolved independently in over 35 plant lineages, accounting for ~ 6% of all higher plants. Agave species are highly heat- and drought-tolerant, and have been domesticated as model CAM crops for beverage, fiber, and biofuel production in semi-arid and arid regions. However, the genomic basis of evolutionary innovation of CAM in genus Agave is largely unknown. RESULTS: Using an approach that integrated genomics, gene co-expression networks, comparative genomics and protein structure analyses, we investigated the molecular evolution of CAM as exemplified in Agave. Comparative genomics analyses among C3, C4 and CAM species revealed that core metabolic components required for CAM have ancient genomic origins traceable to non-vascular plants while regulatory proteins required for diel re-programming of metabolism have a more recent origin shared among C3, C4 and CAM species. We showed that accelerated evolution of key functional domains in proteins responsible for primary metabolism and signaling, together with a diel re-programming of the transcription of genes involved in carbon fixation, carbohydrate processing, redox homeostasis, and circadian control is required for the evolution of CAM in Agave. Furthermore, we highlighted the potential candidates contributing to the adaptation of CAM functional modules. CONCLUSIONS: This work provides evidence of adaptive evolution of CAM related pathways. We showed that the core metabolic components required for CAM are shared by non-vascular plants, but regulatory proteins involved in re-reprogramming of carbon fixation and metabolite transportation appeared more recently. We propose that the accelerated evolution of key proteins together with a diel re-programming of gene expression were required for CAM evolution from C3 ancestors in Agave.

Multimodal MSI in Conjunction with Broad Coverage Spatially Resolved MS2 Increases Confidence in Both Molecular Identification and Localization.

One critical aspect of mass spectrometry imaging (MSI) is the need to confidently identify detected analytes. While orthogonal tandem MS (e.g., LC-MS2) experiments from sample extracts can assist in annotating ions, the spatial information about these molecules is lost. Accordingly, this could cause mislead conclusions, especially in cases where isobaric species exhibit different distributions within a sample. In this Technical Note, we employed a multimodal imaging approach, using matrix assisted laser desorption/ionization (MALDI)-MSI and liquid extraction surface analysis (LESA)-MS2I, to confidently annotate and localize a broad range of metabolites involved in a tripartite symbiosis system of moss, cyanobacteria, and fungus. We found that the combination of these two imaging modalities generated very congruent ion images, providing the link between highly accurate structural information onfered by LESA and high spatial resolution attainable by MALDI. These results demonstrate how this combined methodology could be very useful in differentiating metabolite routes in complex systems.