Skip to main content

The endohyphal microbiome: current progress and challenges for scaling down integrative multi-omic microbiome research

Abstract

As microbiome research has progressed, it has become clear that most, if not all, eukaryotic organisms are hosts to microbiomes composed of prokaryotes, other eukaryotes, and viruses. Fungi have only recently been considered holobionts with their own microbiomes, as filamentous fungi have been found to harbor bacteria (including cyanobacteria), mycoviruses, other fungi, and whole algal cells within their hyphae. Constituents of this complex endohyphal microbiome have been interrogated using multi-omic approaches. However, a lack of tools, techniques, and standardization for integrative multi-omics for small-scale microbiomes (e.g., intracellular microbiomes) has limited progress towards investigating and understanding the total diversity of the endohyphal microbiome and its functional impacts on fungal hosts. Understanding microbiome impacts on fungal hosts will advance explorations of how “microbiomes within microbiomes” affect broader microbial community dynamics and ecological functions. Progress to date as well as ongoing challenges of performing integrative multi-omics on the endohyphal microbiome is discussed herein. Addressing the challenges associated with the sample extraction, sample preparation, multi-omic data generation, and multi-omic data analysis and integration will help advance current knowledge of the endohyphal microbiome and provide a road map for shrinking microbiome investigations to smaller scales.

Video Abstract

Background

Communities of microbes that exist in a particular environment, often referred to as microbiomes, have been universally observed across diverse ecosystems and biological niches including food products, soils, humans, and animals [1,2,3]. Research into the roles of microbiomes has revealed that they perform important ecological functions, through direct impacts on their biological hosts in the case of the human microbiome, and directly within their natural environment in the case of soil microbiomes [2, 4]. Microbiome research has rapidly expanded over the past decade, primarily due to advances in multi-omics and biotechnology, which have enabled investigations at various scales of physical size and complexity. Consequently, this has altered our understanding of microbiomes and the concept of a holobiont — the idea that an organism and the compendium of associated microbes should be considered as a singular entity [5, 6]. The term holobiont has been extensively used to describe humans, plants, and animals; however, progression within the field also led to the notion that microorganisms, often only considered constituents of larger microbiomes, can harbor their own microbiomes. The concept of “microbiomes within microbiomes” has steered the field towards considering the diversity and functional roles of smaller-scale microbiomes (e.g., intracellular scale or smaller physical scale, with lower symbiont biomass) and their impacts on host functioning as well as the dynamics of holobionts within larger and more complex microbial consortia [7].

Fungi are both common and integral members of environmental and host-associated microbiomes, and their presence and functional contributions are known to significantly affect microbiome dynamics and larger ecological processes [8, 9]. The biological complexity of fungi is often underestimated; however, fungi produce a number of complex structures as part of their life cycle, and they contain an array of organelles, lipid droplets, and other intracellular components. Furthermore, filamentous fungi are capable of harboring bacteria (including cyanobacteria), mycoviruses, and other fungi and can even internalize whole microalgal cells within their hyphae [10,11,12,13,14] (Fig. 1). As the diversity of observed endohyphal associates continues to expand, co-occurrences of these constituents have been more frequently observed, and it is becoming evident that many fungi harbor their own endohyphal microbiomes [15]. The endohyphal microbiome specifically refers to the microbiome found within living fungal hyphae and is the topic of this review. However, it is important to note that there are many known epihyphal associations where bacteria and other microorganisms live outside fungal hyphae in close physical proximity. More broadly, the mycosphere can be defined as the zone surrounding fungal mycelia that is characterized by elevated diversity and biological activity and is known to impact fungal physiology and ecosystem processes [16, 17]. As mentioned, this review will specifically discuss the endohyphal microbiome, but the ability to distinguish between epihyphal and endohyphal associates and their individual impacts on fungal host physiology and interactions in the mycosphere is an area of ongoing research in the field.

Fig. 1
figure 1

The known inhabitants of the endohyphal microbiome (microalgae, mycoviruses, bacteria [including cyanobacteria], and fungi)

Multi-omics techniques have rapidly gained popularity for investigating individual constituents and smaller subsets of fungal-associating communities but have yet to be used to interrogate entire endohyphal microbiomes. Due to their early discovery relative to other members of the endohyphal microbiome, bacterial endosymbionts represent the most frequently studied group of fungal associates. In line with canonical microbiome function, endohyphal bacterial associates have been shown to significantly alter fungal host function, development, and interactions with other organisms [18,19,20,21]. Certain endohyphal bacterial associates have been well studied, and thousands of bacterial-fungal interaction pairs (some of which were shown to be endohyphal relationships) have been described [10]. Mycoviruses represent the second most widely studied group within the endohyphal microbiome, and currently, members of 23 viral families have been observed associating with hosts spanning the fungal tree of life [22]. While the total diversity contained within the endohyphal microbiome and its impact on fungal host biology remains unclear, further investigation into these areas will provide a more complete understanding of fungal contributions, interactions, and roles within microbiomes.

Multi-omics approaches (genomics, transcriptomics, proteomics, and metabolomics) have been foundational in elucidating the diversity and functional roles of the endohyphal microbiome. However, these investigations have also highlighted several challenges associated with performing integrative multi-omics studies on small-scale microbiomes. The development of methods and standards relating to sample preparation, data generation, data analysis, and data integration at these small physical (micro to nano) scales will be critical towards enabling holistic investigations of the endohyphal and other small-scale microbiomes.

Many recent microbiome studies using multi-omics methods have centered around scaling up and capturing snapshots of entire environments or ecosystems. Interrogating large-scale microbiomes using multi-omics is informative for analyzing high-abundance sequences and biomolecules but can lack the resolution required for capturing contributions of potential underlying microbiomes of the constituents. In order to interrogate smaller, but potentially equally biologically relevant microbiomes, current techniques and tools for larger-scale investigations must be adjusted, and the development of novel protocols and methods is also warranted. Herein, progress on the utilization of integrative multi-omics for endohyphal microbiome investigations will be discussed along with current and future challenges relating to the application of these techniques to holistic studies of small-scale microbiomes.

Past and current usage of multi-omics for the endohyphal microbiome

Multi-omics investigations have provided key insights regarding the presence, taxonomic composition, and functional implications of members of the endohyphal microbiome on fungal hosts including impacts on host reproduction and pathogenesis [20, 23, 24] (Fig. 2). The publications portrayed in Fig. 2 can also be used as road maps for more specific methodological information for successfully conducting either individual or multi-omics-based investigations of members of the endohyphal microbiome.

Fig. 2
figure 2

Combinations of omics techniques utilized in endohyphal microbiome studies. The publications listed in each Venn diagram section are not exhaustive: one example was chosen for each omics combination that had multiple reference options. The Moebius et al. (2014) publication is highlighted in two sections, and no publications have yet used the three remaining combinations of omics types (genomics + proteomics + metabolomics, transcriptomics + proteomics + metabolomics, and genomics + transcriptomics + proteomics + metabolomics) to investigate the endohyphal microbiome

Endohyphal bacteria

Multi-omics investigations on endohyphal microbiomes have been most widely conducted on members of the Mucoromycota, which frequently harbor members from either one or both of two groups of highly studied bacterial endosymbionts: Mollicutes-related endobacteria (MRE) and Burkholderia-related endobacteria (BRE). Other fungal lineages, including Ascomycota and Basidiomycota, have been found to host endohyphal bacteria, but the distribution of these bacteria within fungal mycelial networks and their impacts on host physiology and functioning have been comparatively understudied [25,26,27].

Genomics- and metagenomics-based studies determined that MRE and BRE both have significantly reduced genome sizes and protein-coding content compared to closely related free-living bacteria, as well as reduced metabolic capacities [28,29,30,31]. These investigations have also revealed instances of horizontal gene transfer (HGT) between endobacteria and fungal hosts [29, 30]. Metagenomics has been utilized for simultaneous evolutionary and comparative examinations of both the bacterial endosymbiont and the fungal host, leading to insights into the origin of fungal-endobacteria interactions [28, 30, 32]. These genomic and metagenomic studies have provided fundamental knowledge around the endohyphal bacterial lifestyle and have demonstrated how these interactions can lead to genomic alterations of both the host and the endobacteria [33].

Transcriptomics has been used to uncover the genes and putative functions involved in interactions between endohyphal bacteria and their hosts. For example, comparative transcriptomic experiments between Rhizopus microsporus (Mucoromycota) isolates harboring BRE and those cured of their BRE partner revealed that the bacterial endosymbiont alters host gene expression, allowing it to control aspects of both asexual and sexual reproduction [34]. Comparative transcriptomic experiments have also provided insights into how a bacterial endosymbiont (Luteibacter, Gammaproteobacteria) interacts with its fungal host (Pestalotiopsis, Ascomycota) prior to endohyphal establishment [35]. Additionally, integrated analyses of transcriptome data with other omics data types have provided a more complete picture of the functional consequences of hosting endohyphal bacteria. For example, transcriptomic and proteomic data revealed that BRE in Gigaspora margarita (Mucoromycota) cause a shift in primary metabolism of the fungal host and an increase in fungal antioxidant production, which in turn alters interaction dynamics between plants and the fungal host [36]. Examination of the R. microsporus transcriptome and metabolome indicated differences in lipid metabolism and composition during interactions with its BRE partner [37]. To investigate the mechanisms by which bacteria can enter fungal host cells, genome mining, transcriptomics, and proteomics have been used to determine the critical role of chitinase during endosymbiosis establishment by Mycetohabitans rhizoxinica (Betaproteobacteria) in R. microsporus [38].

Metabolomics studies have identified broad shifts in the fungal metabolome related to harboring endobacteria. For example, the presence of BRE was found to induce substantial shifts in the metabolome of G. margarita [39]. Metabolomics has also been used in conjunction with genomics to determine that the presence of Paenibacillus (Bacilli) alters the lipid profiles of Thamnidium elegans (Mucoromycota) compared to other Mucoromycota fungi which do not host the bacterium [40]. Similarly, metabolomics was used in conjunction with genomics and transcriptomics in a study that determined that the presence/absence of a BRE (Mycoavidus cysteinexigens) substantially alters the Linnemannia elongata (Mucoromycota) transcriptome and metabolome, particularly with respect to fatty acid biosynthesis [28]. In addition to characterizing broad shifts in metabolomes, metabolomics has been used to identify important individual metabolites. For example, metabolomics was utilized to determine that rhizoxin, the causative toxin of rice seedling blight, is produced by a Rhizopus-associated BRE, not the fungus itself [41]. Proteomic and metabolomic investigations comparing spores of G. margarita with and without bacterial endosymbionts have shown that the endobacteria can significantly alter host protein expression and lipid profiles under several conditions and growth stages [42]. Investigations which integrated metabolomics and proteomics data identified key substrates exchanged between Linnemannia elongata (Mucoromycota) and its bacterial endosymbiont Mycoavidus sp. (Betaproteobacteria) [43]. Volatile emission measurements from L. elongata with and without Mycoavidus sp. showed that when the endobacterium is present, the production of fatty acids from the pyruvate pathway leads to high butyric acid and butyrate levels [44].

Other members of the endohyphal microbiome

Exploratory studies involving other members of the endohyphal microbiome have lagged behind in numbers, but these studies have led to key findings. Multi-omic studies have facilitated the identification and corroboration of novel members of the endohyphal microbiome and have started to provide insights into how members within the microbiome interact with one another in addition to their interactions with the fungal host. Genomics and transcriptomics were used to explore genes of interest, potential HGT, gene expression, and phylogeny of the fungal host Geosiphon pyriformis (Mucoromycota), the first fungus found to internalize Nostoc, a nitrogen-fixing cyanobacterium [45]. Stable isotope labeling was used in the first ever study that observed the internalization of whole algal cells into fungal hyphae, which demonstrated a reciprocal exchange of nitrogen and carbon between the symbionts [14].

In addition to beneficial symbiotic partnerships, other endohyphal microbiome constituents such as viruses or parasitic fungi can modulate fungal host function. Transcriptomic analyses were used to identify the molecular changes of the mycoparasite Ampelomyces quisqualis (Ascomycota) during the colonization of its host Podosphaera xanthii (Ascomycota) [46]. Mycoviruses are typically single-stranded or double-stranded RNA viruses; therefore, RNA sequencing is needed to identify and classify most mycoviruses [12]. Recent surveys of publicly available fungal transcriptome sequencing data have revealed that mycoviruses are taxonomically diverse and are common among members of several fungal phyla [47, 48]. Transcriptome sequencing of G. margarita has also been used to investigate potential modulation of mycovirus diversity as a result of bacterial presence in the endohyphal microbiome [49]. Transcriptomics, proteomics, and metabolomics have been used to elucidate phenotypic and functional outcomes of mycovirus infections on fungal hosts [50,51,52].

While multi-omics techniques have advanced current knowledge on the diversity and impacts of the broader endohyphal microbiome, many of these studies focus on individual microbiome members. Due to the demonstrated utility of multi-omics and the growing recognition of the complexity of the endohyphal microbiome, we suggest the routine use of multi-omics in examining the full diversity of endohyphal microbiome composition and function over time. To initiate these more complex studies, it is necessary to identify and work to overcome the current limitations of these methods when applied to investigations of small-scale microbiomes.

Challenges for integrative multi-omics of the endohyphal microbiome and other small-scale microbiomes

Sample harvesting and extraction

The physical isolation of holobionts from their environments enables small-scale investigations without concerns of contributions from the larger microbiomes. The ability to isolate a fungal host in culture reduces or eliminates external contamination and unassociated genetic and molecular signals, permitting high-resolution investigations of the endohyphal microbiome [53]. However, this approach biases comprehensive multi-omic investigations into endohyphal microbiomes to only culturable fungi. Recent advances in laboratory devices and techniques have increased the efficiency of isolating fungi from larger microbial communities. Selective media, media additives, and growth conditions (e.g., incubation temperature) can be used to culture specific fungal groups [8]. Technical advances in single-cell and low-input biomass techniques have also enabled fungal investigations at smaller physical scales, allowing interrogation of hyphal fragments and micro-scale spores such as conidia.

When endohyphal microbiome constituents are separated from their fungal host and the rest of the microbiome, multi-omic analyses can be more straightforward and can result in higher quality data [54, 55]. However, it may be challenging or impossible to isolate or culture certain members of the endohyphal microbiome, as some associates may be unculturable and others are known to be dependent on the fungal host, such as the auxotrophic endohyphal bacterium Mycoavidus cysteinexigens which relies on its host for production of cysteine [28, 56]. Additionally, the process of separating these components leads to perturbations of the native sampling environment which can alter what is captured by transcriptomic, metabolomic, and proteomic approaches and lead to results that are not indicative of endohyphal functions. Therefore, separation of the microbiome components prior to enrichment or extraction may not always be preferable. It can also be difficult to axenically isolate the fungal host, as the success of methods for curing fungal hosts of their endohyphal bacteria remains highly variable. Current methods can be host or associate specific, and no standard methods exist for curing a fungal host of its entire endohyphal microbiome [26, 27]. This in turn limits comparative multi-omic experiments between cured and uncured fungal hosts to directly measure and assess impacts of the endohyphal microbiome.

For sequence-based interrogations, the nucleic acid extraction method can impact extraction efficiency of the endohyphal microbiome. This is partly due to differences in cell lysis requirements for host fungi and their microbiome constituents. Extraction kits for various microbiomes can differentially affect results, particularly in regards to nucleic acid recovery and quality [57, 58]. To identify and characterize mycoviruses present within the endohyphal microbiome, RNA extractions are often required. Total nucleic acid kits for simultaneous extraction of DNA and RNA have become more commonplace, especially for integrative multi-omics experiments. Although these dual-extraction kits can be very effective, reports have also documented biases with some of these extraction methods [59]. Various protocols exist for performing high-molecular-weight extractions on fungi for long-read sequencing, but not all of these protocols work efficiently for all lineages and cell types, and nucleic acid yields and quality can be highly variable [60, 61].

Metabolite and bottom-up protein extraction methods and purification efficiency can also vary for small-scale microbiomes. The selection of cellular disruption (e.g., physical or chemical) and metabolome or proteome extraction methods must be consistent with the platform used for data generation, and different methods can bias results or impact biomolecule recovery [62, 63]. Protein extraction methods can have variable impacts on protein yields and accuracy of protein identification; thus, it is important to evaluate which method may be most appropriate based on the experimental design and scientific goals [64]. Protein extractions from fungi can also be especially challenging given the robustness of fungal cell walls and fungal secretion of proteases which can cause protein degradation [65]. This becomes even more challenging when considering lysis requirements for both the fungal host and its associated microbiota. Metabolite extractions also rarely capture all metabolites in the sample, which can lead to loss in overall diversity or the need for multiple extraction methods [63]. Several published techniques, for example, the MPLEx (metabolite, protein, and lipid extraction) method, have pioneered the parallel extraction of proteins and metabolites from the same starting sample, thus making multi-omic experiments more integrative, less time-consuming, and more informative on a spatiotemporal scale for each particular sample [66,67,68,69,70].

Sampling methods for extractions must also be considered when working with diverse fungal isolates. Filamentous fungi can vary substantially in their growth rates and hyphal division and production. All of these factors can impact the endohyphal microbiome, as microbiome constituents may replicate and disperse at different rates, creating differences in abundance across older and younger hyphae or hyphal segments. Certain members of the endohyphal microbiome may be present at very low levels, and not all members of the microbiome will have even distributions throughout the host’s hyphal network. Partial sampling of the fungal host, as opposed to collecting all available hyphae, could result in differential capture of microbiome components. Multiple extractions from total cultured fungal biomass or large portions of biomass at multiple locations across the hyphal network may be required to capture the diversity of heterogeneously distributed microbiome features. Spatiotemporal dynamics are important to consider as expression of transcripts, proteins, and metabolites will fluctuate on different time scales and differentially across individual hyphae. This can result in highly variable technical extraction replicates due to inconsistent capture of endohyphal constituent signatures. No extraction methods currently exist for the simultaneous extraction of DNA, RNA, proteins, and metabolites, making it important to consider splitting and preserving samples for each omic type. One of the biggest remaining challenges associated with performing multi-omics on the endohyphal microbiome is the consistently low yield of biomolecules produced by the microbiome constituents relative to the fungal host, which can be a significant hurdle for downstream multi-omics sample processing if the target signal is simply insufficient to provide robust or significant results. Figure 3 summarizes the aforementioned challenges as well as those associated with each step of the multi-omics experimental and analytical process.

Fig. 3
figure 3

Challenges, considerations, and limitations for conducting integrative multi-omics experiments on the endohyphal microbiome. The overall multi-omics workflow consists of sample harvesting, sample extraction, and sample preparation, as well as data generation, data analysis, and data integration steps. The blue arrow denotes progression through the experimental process (from sample isolation to integration of multi-omics data), and the orange triangle denotes the decrease in sample yield as the workflow progresses from sample harvesting and extraction to sample preparation. Created with BioRender.com; adapted from “Multi-Panel Horizontal Timeline (Layout 2 × 2),” by BioRender.com (2023). Retrieved from https://app.biorender.com/biorender-templates

Sample preparation for multi-omics

Some considerations for the endohyphal microbiome sample preparation stage are relevant for all omics types, such as the consistent challenge of low input of molecular targets from the microbiome, compared with the overwhelming molecular targets from the fungal host. For investigations of nucleic acids, this challenge can be addressed through sequence enrichment, host depletion, and/or target amplification. Amplification techniques are the most straightforward and cost-effective of these options, as many primer options exist that will generally target barcode sequences for identification and taxonomic classification of bacteria/cyanobacteria (16S rRNA gene), viruses (RdRp, RNA-dependent RNA polymerase), fungi (ITS, internal transcribed spacer), and algae (18S rRNA gene). The use of these somewhat “universal” primers also means there is no reliance on prior knowledge of the specific diversity of endohyphal microbiome constituents, and this type of analysis provides a general overview of the community profile of the endohyphal microbiome. However, this can require multiple sample preparations with multiple amplification steps which can introduce additional biases and noise, and the results will not provide any genomic information outside of these barcode sequences. Whole genome or whole transcriptome amplification (WGA or WTA) methods can also be used to amplify more than the standard amplicon barcode regions. While newer methods have minimized biases (e.g. GC content biases), chimera generation, and artifacts that can result from these protocols, these issues are still observed to occur and can significantly affect mixed community samples [71].

Alternatively, enrichment and/or depletion assays can be performed. Hybridization-based genomic enrichment strategies have been used in microbiome research to select for known sequences based on lineage-specific signatures from organisms of interest [72]. However, this strategy relies on prior knowledge of taxonomic and sequence diversity within the sample, as efficiency is driven by sequence similarity between the enrichment probes and sequences in the sample. The design of specific enrichment and host-depletion panels can be quite time consuming and cost-prohibitive, especially when trying to encompass the potentially extensive diversity of the endohyphal microbiome. Our current limited knowledge on the diversity within the endohyphal microbiome makes it challenging to design enrichment panels since non-targeted sequences will generally not be captured or sequenced.

For fungal transcriptome investigations, it is common to deplete rRNAs or to perform polyA selection, as total RNA extractions are dominated by fungal host rRNAs [73]. PolyA selection limits complete characterization of the microbiome as non-polyadenylated sequences from viral or prokaryotic members are discarded, and this method can be less efficient at reducing rRNA levels compared to direct rRNA depletion methods [74, 75]. Furthermore, rRNA, polyA, and other forms of depletion require large amounts of starting material, and they often result in significant losses of total RNA. Considering the low proportions of RNA from members of the endohyphal microbiome relative to the fungal host, the resulting RNA yields may be too low for downstream sample and library preparation and subsequent analysis [75]. If the sample also contains abundant endohyphal bacteria, bacterial rRNA will likely dominate the remaining transcript pool post-depletion. Bacterial rRNA can also be depleted but will once again result in an additional loss of total RNA yield [76]. Library preparation for genomic and transcriptomic sequencing has historically been severely hindered by low nucleic acid input amounts. Newer methods for single-cell sequencing and low-input library preparation have helped to alleviate this shortcoming; however, these methods are still prone to errors, technical noise, and biases and can have higher failure rates [77,78,79].

Similar to nucleic acids, proteins and metabolites produced by the fungal host will be much more prevalent than biomolecules from endohyphal microbiome constituents. Compared with nucleic acids, proteins and metabolites cannot be amplified in a similar manner, but differences in the physical and chemical properties, including charge, hydrophobicity, and molecular weight, of biomolecules of interest (e.g., lipids, pigments, organic acids) can be exploited for selective enrichment. Recent advances in enrichment for single-cell omics show promise for improving the detection of low abundance analytes from endohyphal microbiomes. Capillary electrophoresis (CE) is commonly used for low input metabolite and protein selection (for bottom-up proteomics) based on biomolecule size and charge [80, 81]. While these enrichment techniques can select for specific classes of biomolecules, they also reduce total biomolecule amounts. This is particularly concerning with proteomic sample preparation, as the steps involved often lead to significant losses of proteins and peptides, and this is further exacerbated with lower starting inputs [82]. Many low-input proteomics protocols involve minimizing the sample preparation steps to reduce the sample loss characterized in conventional proteomics methods, often during tube transfer [82,83,84].

Multi-omic data generation

The generation of sequencing data from endohyphal microbiome samples has unique considerations depending on the input sample type and research goals. Given the challenges listed above, metagenomic sequencing of a total fungal culture extract will often require deep sequencing to adequately assess the total diversity of the endohyphal microbiome. Deep sequencing can be very cost prohibitive, and most of the sequenced reads will belong to the fungal host. Rare members of the endohyphal microbiome may be captured at very low levels or not at all. Microbiome-enriched and/or fungal host-depleted metagenomic samples will reduce this need for ultra-deep sequencing. However, eukaryotic members of the endohyphal microbiome such as algae and other fungi can have much larger genomes compared to other members of the microbiome and require sufficient sequencing depth to acquire adequate genomic coverage.

The use of either short- or long-read sequencing technologies is also important to consider. Library preparation kits for short-read sequencing platforms generally require less sample input amounts than long-read platforms. However, single-cell/low-input library kits and workflows currently exist for both platforms for DNA, while only short-read platforms offer similar options for RNA samples. Short-read platforms are more commonly used for microbiome analyses, as their higher average sequencing depths and per-sample read counts are often superior for capturing rare microbiome members [85]. However, the use of long-read sequencing platforms in microbiome research is becoming more common. While long-read sequencing platforms may not be as efficient at capturing rare members of the microbiome, longer contigs collected from this platform aid in the generation of more complete genome references [86]. When possible, a combination of both long- and short-read sequencing platforms is ideal for these sample types [85].

Selection of an instrument or technique for performing metabolomic and proteomic analyses comes with its own set of challenges, as this requires considering methods best suited for investigating fungal hosts, microbiomes, and/or low-input (e.g., single cell) samples. Untargeted methods for both metabolomics and proteomics capture the most analyte diversity and capture unknown or undefined biomolecules produced by the fungal holobiont [87]. Techniques for performing untargeted low-input metabolomics and proteomics are still relatively nascent, but recent advancements have made this field more practical for small-scale microbiomes [88, 89]. The throughput of the technology must also be carefully considered (e.g., time per sample, amount of possible multiplexing), as rare analytes across samples may not be seen with lower throughput methods, even if the resolution per sample is ideal [90]. Mass spectrometry (MS) technologies are comparatively widely used for proteomic and metabolomic studies involving samples with low volumes and low biomolecule abundances, and MS instruments are continually improving in sensitivity and mass resolution [91].

Matrix-assisted laser desorption/ionization MS (MALDI-MS), secondary ion MS (SIMS), and electrospray ionization MS (ESI–MS) are currently some of the most popular low-input or single-cell metabolite and protein analysis platforms [92]. MALDI-MS has been widely used for single-cell metabolomics and proteomics, and it has shown promise for detecting analytes expressed from rare cells within a larger population [93]. However, MALDI-MS requires extensive sample preparation steps that may lead to additional sample loss and may compromise the natural metabolic states of the endohyphal microbiome constituents, and the technique does not perform well with low-molecular-weight organic acids [89]. SIMS provides one of the best possible cellular resolution ranges of currently available instrumentation, making it ideal for low input samples. However, it is still limited on mass resolution and can be a lower throughput technique. ESI–MS, especially when coupled with capillary electrophoresis (CE-MS), is also often used to perform low-input and small-scale metabolomics, even to the subcellular level [94, 95].

Omics techniques that also incorporate an imaging aspect can be critical contributions to investigations of microbiomes, including the endohyphal microbiome. Importantly, microscopy techniques are some of the best methods for distinguishing between members of the endohyphal and epihyphal communities. Imaging techniques also provide insights regarding the spatial heterogeneity within the fungal host and can be utilized to understand where biomolecules originate and end up within the fungal holobiont. Probes for fluorescence in situ hybridization (FISH) are traditionally designed from genomic and transcriptomic sequencing data. FISH imaging has been widely utilized for localization of endohyphal bacteria, and the technique could be expanded for the simultaneous detection and visualization of multiple members of the endohyphal microbiome [96, 97]. Single molecule FISH (smFISH) has also become much more reliable and applicable to rare microbiome constituents, and new methods have also increased the throughput of this approach [98]. Challenges and methods for visualizing endobacteria inside fungal hyphae using FISH have been previously discussed [97, 99]. Transmission electron microscopy has been used to identify MRE and BRE among other structures in the cytoplasm of Mortierellaceae and Glomeromycotina taxa to confirm their presence within fungal hyphae. Electron microscopy, in particular CryoEM techniques, continue to be invaluable for investigating the endohyphal microbiome [100, 101].

MALDI-MS can include an imaging aspect, which has even been used for 3D imaging within these small physical scales [102]. Recently developed metabolite imaging methods such as SpaceM directly correlate metabolomics data (e.g., MALDI data) with microscopy imaging to gain the spatial context of the cellular origin of the metabolites [103]. Various imaging techniques are difficult to integrate, as they each require unique sample preparation protocols (e.g., fixation); however, incorporating imaging methods into multi-omics studies will undoubtedly be indispensable for endohyphal microbiome studies. Imaging techniques also continue to be the most reliable way to validate which microbiome components are within the bounds of the fungal hyphae in order to define the endohyphal microbiome.

Multi-omic data processing and integration

Obtaining genome sequences for members of the endohyphal microbiome is often the first step in multi-omics data analysis. These sequences serve as a reference for transcriptomics and proteomics and provide information on taxonomic diversity and functional potential, such as the presence or absence of metabolic pathways, within the endohyphal microbiome. If members of the endohyphal microbiome can be separated from the fungal host through filtration or centrifugation and sequenced individually, the genomic sequencing datasets from each individual can be analyzed using bioinformatic tools and/or pipelines tailored to their specific taxonomy (e.g., prokaryotes or eukaryotes). In cases where separation is not possible, metagenomic sequencing of endohyphal microbiomes will result in mixed datasets containing sequences from all members at variable coverage, with lower abundance constituents having the least amount of sequence coverage. Sequences from metagenomic datasets can be assembled, taxonomically classified, and assessed for quality using a number of established and taxa-specific approaches [104, 105]. Tools have recently been developed to selectively identify eukaryotic and viral contigs from shotgun metagenomics data to aid in the separation and subsequent analyses of certain contigs [106,107,108]. More specific tools have also shown promise for these analyses: the Spore-associated Symbiotic Microbes (SeSaMe) bioinformatics tool was specifically developed for the sequence classification of microbial associates of arbuscular mycorrhizal fungi from metagenomic sequencing data [109]. For endohyphal microbiome constituents that are not adequately captured in shotgun metagenomic assemblies, more sensitive read-based taxonomic classification analyses can be used to detect their presence; however, the lack of contigs limits the analyses of other omics types to de novo (without a reference) methods. While a few genome assembly software packages have been developed for low input and single-cell applications, substantial optimization hurdles remain [110, 111].

Genome assemblies can be especially difficult for endohyphal bacteria due to their aforementioned rapid rates of evolution, divergence from available bacterial references, and the possible presence of multiple closely related populations within a single fungal host [29, 32, 112]. HGT can also occur as a result of the intimate interactions between endobacteria and their fungal hosts; furthermore, the mycosphere and fungal mycelia have been shown to be a niche experiencing increased bacterial HGT [16, 29, 30]. This creates additional challenges when assembling multiple low-coverage genomes simultaneously, as HGT sequences may be incorrectly assembled, and determining the origin of HGT sequences may not be trivial. Traditional techniques to evaluate genome completeness and quality may also not be appropriate for endosymbiont genomes, for the reasons stated above.

Adequate genome references are a traditional cornerstone for integrative multi-omics. Incomplete or low-quality reference genome assemblies can make it challenging to utilize transcriptomic analysis software and methods that rely on reference genomes. Tools for de novo transcriptome assembly such as the Trinity software can make transcriptomic analysis without reference genomes more feasible [113]. Trinity has been successfully utilized for de novo assemblies of transcriptomes and metatranscriptomes containing signatures from multiple kingdoms [114,115,116]. Other tools and methods for separating host and symbiont reads in holobiont systems have been designed that minimize the chimeras that can typically arise from de novo transcriptome analyses [117]. Assemblers specifically created for metatranscriptome samples may also be the most appropriate for diverse endohyphal microbiomes [118]. Importantly, the application of de novo transcriptome assembly methods to the endohyphal microbiome relies on sequencing datasets containing sufficient coverage of transcripts from each microbiome constituent. In the case of the endohyphal microbiome, it is also important to consider that samples obtained from different locations across the hyphal network or different timepoints could vary in their microbiome composition and the resulting endohyphal transcriptome. Performing metagenomic and metatranscriptomic analyses on the exact same sample is ideal; otherwise, artifacts or discrepancies may arise. Given that mycoviruses are typically RNA viruses, their assembly and classification from metatranscriptome data will require specific viral workflows and databases that contain sequences from other mycoviruses.

Proteomic analyses also heavily rely on adequate host genome and microbiome metagenome reference data. When targeting microbial consortia, and particularly poorly explored ones such as the endohyphal microbiome, it is common practice to construct a protein reference database from annotated assemblies generated from the community’s metagenomic and occasionally metatranscriptomic data. As discussed in previous sections, metagenomic analyses will often not capture complete genomes for all members of the endohyphal microbiome. To ensure a more complete taxonomic and/or functional assignment of proteomic data, reference genomes from similar taxa to those in the microbiome can be used as a substitute, or recently developed de novo peptide prediction methods can be employed. The use of reference genomes of species similar to those found in the microbiome may not adequately represent the expressed proteins which can lead to incorrect or uninterpretable results, and this method can be especially problematic when assessing rapidly evolving bacterial endosymbionts whose proteins may share very little homology with otherwise closely related taxa. De novo methods for proteomic sequence data analysis may be required, although these are relatively new and lack standardization in the field [119, 120].

Analysis of untargeted metabolomics data can be completed without genomic or metagenomic references; however, unique challenges remain. These untargeted analyses rely on databases for accurate metabolite identification against libraries of known metabolites. However, these databases are known to be largely incomplete, thus leading to large proportions of untargeted metabolomics data having “unknown” classifications [121]. Many databases are highly limited on representation of metabolites from fungi, and studies involving fungi often require analyses which query a number of separate databases [122]. The use of multiple libraries or databases will likely be required to account for all the potential taxonomic diversity found within the endohyphal microbiome, and this can make it challenging to perform holistic analyses.

Multi-omics-based studies are highly complementary, as examinations of the metabolome and proteome indicate changes in function, but genome and transcriptome investigations help reveal the genes and functional potential responsible for these changes. The combination of multiple omics types, and integration of the data, provides clearer interpretations and conclusions regarding the functions within the endohyphal microbiome and their impacts on the fungal host. Omics techniques and data generation have advanced tremendously in the past decade but still vary in their sensitivity and limits of detection. The utilization of multiple omics types helps to overcome the deficiencies of any single technique, and this approach is often required to determine the impacts of the endohyphal microbiome on the fungal host, as altered expression of biomolecules within the endohyphal microbiome can directly impact gene expression or functions of the fungal host. The process of multi-omic data integration still lacks optimization and standardization, although several tools and workflows exist which have been previously discussed and compared [123, 124].

Research into small-scale microbiomes and microorganism holobionts will require additional efforts towards new and optimized algorithms, bioinformatic tools and workflows, and standardization specifically centered around analyses of low-yield interkingdom samples. There has been a very active campaign towards increased standardization in multi-omics bioinformatics workflows; however, significant shortcomings and gaps remain, particularly for microbial community analyses, making it difficult for many researchers to easily perform end-to-end multi-omic analyses [89, 125]. Currently, no standardized methods exist for holistically investigating the endohyphal microbiome using multi-omics, and standardization will be key towards progressing this field and making these insights comparable across different studies (e.g., comparing molecular underpinning of different fungal hosts and their respective microbiomes) and applicable to other small-scale holobionts.

Conclusions

Many advances have been made in integrative multi-omic approaches for investigations into certain members of the endohyphal microbiome. However, significant challenges still remain at each experimental and analytical step that prevent fungal holobionts and other small-scale microbiomes from being routinely and holistically investigated. Integrative multi-omics is challenging even with abounding sample inputs into each omic analysis and with high-quality genome assemblies, and the challenges associated with integrative multi-omics are greatly compounded as additional members of the endohyphal microbiome are simultaneously interrogated and as additional omics types are integrated. The concept of a complex endohyphal microbiome is still relatively new, and the total taxonomic and functional diversity of the endohyphal microbiome has yet to be uncovered. Multi-omics has the potential to vastly expand these investigations, but it is imperative that the microbiome field considers the importance of these microbial holobionts and small-scale microbiomes when designing and advocating for new tools, methods, techniques, and standards. Understanding the functional contributions of each individual microbiome to their respective hosts as well as to larger microbial communities will unlock new opportunities and scientific questions in the field of microbiome research.

Availability of data and materials

Not applicable.

Abbreviations

BRE:

Burkholderia-related endobacteria

CE:

Capillary electrophoresis

ESI:

Electrospray ionization

FISH:

Fluorescence in situ hybridization

HGT:

Horizontal gene transfer

ITS:

Internal transcribed spacer

MALDI:

Matrix-assisted laser desorption/ionization

MRE:

Mollicutes-related endobacteria

MS:

Mass spectrometry

RdRp:

RNA-dependent RNA polymerase

SIMS:

Secondary ion mass spectrometry

smFISH:

Single-molecule fluorescence in situ hybridization

WGA:

Whole genome amplification

WTA:

Whole transcriptome amplification

References

  1. Walsh AM, Macori G, Kilcawley KN, Cotter PD. Meta-analysis of cheese microbiomes highlights contributions to multiple aspects of quality. Nat Food. 2020;1:500–10.

    Article  PubMed  Google Scholar 

  2. Jansson JK, Hofmockel KS. Soil microbiomes and climate change. Nat Rev Microbiol. 2020;18:35–46.

    Article  CAS  PubMed  Google Scholar 

  3. Turnbaugh PJ, Ley RE, Hamady M, Fraser-Liggett CM, Knight R, Gordon JI. The human microbiome project. Nature. 2007;449:804–10.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  4. The Integrative Human Microbiome Project. Nature. 2019;569:641–8.

    Article  Google Scholar 

  5. Berg G, Rybakova D, Fischer D, Cernava T, Vergès M-CC, Charles T, et al. Microbiome definition re-visited: old concepts and new challenges. Microbiome. 2020;8:103.

    Article  PubMed  PubMed Central  Google Scholar 

  6. Margulis L. Words as battle cries: symbiogenesis and the new field of endocytobiology. Bioscience. 1990;40:673–7.

    Article  CAS  PubMed  Google Scholar 

  7. Simon J-C, Marchesi JR, Mougel C, Selosse M-A. Host-microbiota interactions: from holobiont theory to analysis. Microbiome. 2019;7:5.

    Article  PubMed  PubMed Central  Google Scholar 

  8. Huffnagle GB, Noverr MC. The emerging world of the fungal microbiome. Trends Microbiol. 2013;21:334–41.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  9. Wagg C, Schlaeppi K, Banerjee S, Kuramae EE, van der Heijden MGA. Fungal-bacterial diversity and microbiome complexity predict ecosystem functioning. Nat Commun. 2019;10:4841.

    Article  PubMed  PubMed Central  Google Scholar 

  10. Robinson AJ, House GL, Morales DP, Kelliher JM, Gallegos-Graves LV, LeBrun ES, et al. Widespread bacterial diversity within the bacteriome of fungi. Commun Biol. 2021;4:1–13.

    Article  Google Scholar 

  11. Kluge M. A fungus eats a cyanobacterium: the story of the Geosiphon pyriformis endocyanosis. Biology and Environment: Proceedings of the Royal Irish Academy. 2002;102B:11–4.

  12. Kondo H, Botella L, Suzuki N. Mycovirus diversity and evolution revealed/inferred from recent studies. Annu Rev Phytopathol. 2022;60:307–36.

    Article  CAS  PubMed  Google Scholar 

  13. Sun JZ, Liu XZ, McKenzie EHC, Jeewon R, Liu JK, Zhang XL, et al. Fungicolous fungi: terminology, diversity, distribution, evolution, and species checklist. Fungal Diversity. 2019;95:337–430.

    Article  Google Scholar 

  14. Du Z-Y, Zienkiewicz K, Vande Pol N, Ostrom NE, Benning C, Bonito GM. Algal-fungal symbiosis leads to photosynthetic mycelium. Baldwin IT, Harrison MJ, Harrison MJ, Bonfante P, Hom E, editors. eLife. 2019;8:e47815.

  15. Desirò A, Salvioli A, Ngonkeu EL, Mondo SJ, Epis S, Faccio A, et al. Detection of a novel intracellular microbiome hosted in arbuscular mycorrhizal fungi. ISME J. 2014;8:257–70.

    Article  PubMed  Google Scholar 

  16. Zhang M, e Silva M de C, Chaib De Mares M, van Elsas JD. The mycosphere constitutes an arena for horizontal gene transfer with strong evolutionary implications for bacterial-fungal interactions. FEMS Microbiology Ecology. 2014;89:516–26.

    Article  CAS  PubMed  Google Scholar 

  17. Pratama AA, van Elsas JD. Gene mobility in microbiomes of the mycosphere and mycorrhizosphere –role of plasmids and bacteriophages. FEMS Microbiology Ecology. 2019;95:fiz053.

    Article  CAS  PubMed  Google Scholar 

  18. Pawlowska TE, Gaspar ML, Lastovetsky OA, Mondo SJ, Real-Ramirez I, Shakya E, et al. Biology of fungi and their bacterial endosymbionts. Annu Rev Phytopathol. 2018;56:289–309.

    Article  CAS  PubMed  Google Scholar 

  19. Büttner H, Niehs SP, Vandelannoote K, Cseresnyés Z, Dose B, Richter I, et al. Bacterial endosymbionts protect beneficial soil fungus from nematode attack. Proc Natl Acad Sci. 2021;118:e2110669118.

    Article  PubMed  PubMed Central  Google Scholar 

  20. Itabangi H, Sephton-Clark PCS, Tamayo DP, Zhou X, Starling GP, Mahamoud Z, et al. A bacterial endosymbiont of the fungus Rhizopus microsporus drives phagocyte evasion and opportunistic virulence. Curr Biol. 2022;32:1115–1130.e6.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  21. Richter I, Radosa S, Cseresnyés Z, Ferling I, Büttner H, Niehs SP, et al. Toxin-producing endosymbionts shield pathogenic fungus against micropredators. mBio. 2022;13:e01440–22.

    Article  PubMed  PubMed Central  Google Scholar 

  22. Hough B, Steenkamp E, Wingfield B, Read D. Fungal viruses unveiled: a comprehensive review of mycoviruses. Viruses. 2023;15:1202.

    Article  PubMed  PubMed Central  Google Scholar 

  23. Kotta-Loizou I. Mycoviruses and their role in fungal pathogenesis. Curr Opin Microbiol. 2021;63:10–8.

    Article  CAS  PubMed  Google Scholar 

  24. Partida-Martinez LP, Monajembashi S, Greulich K-O, Hertweck C. Endosymbiont-dependent host reproduction maintains bacterial-fungal mutualism. Curr Biol. 2007;17:773–7.

    Article  CAS  PubMed  Google Scholar 

  25. Arendt KR, Hockett KL, Araldi-Brondolo SJ, Baltrus DA, Arnold AE. Isolation of endohyphal bacteria from foliar Ascomycota and in vitro establishment of their symbiotic associations. Appl Environ Microbiol. 2016;82:2943–9.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  26. Sharma M, Schmid M, Rothballer M, Hause G, Zuccaro A, Imani J, et al. Detection and identification of bacteria intimately associated with fungi of the order Sebacinales. Cell Microbiol. 2008;10:2235–46.

    Article  CAS  PubMed  Google Scholar 

  27. Glaeser SP, Imani J, Alabid I, Guo H, Kumar N, Kämpfer P, et al. Non-pathogenic Rhizobium radiobacter F4 deploys plant beneficial activity independent of its host Piriformospora indica. ISME J. 2016;10:871–84.

    Article  PubMed  Google Scholar 

  28. Uehling J, Gryganskyi A, Hameed K, Tschaplinski T, Misztal PK, Wu S, et al. Comparative genomics of Mortierella elongata and its bacterial endosymbiont Mycoavidus cysteinexigens. Environ Microbiol. 2017;19:2964–83.

    Article  CAS  PubMed  Google Scholar 

  29. Torres-Cortés G, Ghignone S, Bonfante P, Schüßler A. Mosaic genome of endobacteria in arbuscular mycorrhizal fungi: transkingdom gene transfer in an ancient mycoplasma-fungus association. Proc Natl Acad Sci. 2015;112:7785–90.

    Article  PubMed  PubMed Central  Google Scholar 

  30. Naito M, Morton JB, Pawlowska TE. Minimal genomes of mycoplasma-related endobacteria are plastic and contain host-derived genes for sustained life within Glomeromycota. Proc Natl Acad Sci. 2015;112:7791–6.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  31. Lackner G, Moebius N, Hertweck C. Endofungal bacterium controls its host by an hrp type III secretion system. ISME J. 2011;5:252–61.

    Article  CAS  PubMed  Google Scholar 

  32. Sun X, Chen W, Ivanov S, MacLean AM, Wight H, Ramaraj T, et al. Genome and evolution of the arbuscular mycorrhizal fungus Diversispora epigaea (formerly Glomus versiforme) and its bacterial endosymbionts. New Phytol. 2019;221:1556–73.

    Article  CAS  PubMed  Google Scholar 

  33. Sharmin D, Guo Y, Nishizawa T, Ohshima S, Sato Y, Takashima Y, et al. Comparative genomic insights into endofungal lifestyles of two bacterial endosymbionts, Mycoavidus cysteinexigens and Burkholderia rhizoxinica. Microbes Environ. 2018;33:66–76.

    Article  PubMed  PubMed Central  Google Scholar 

  34. Mondo SJ, Lastovetsky OA, Gaspar ML, Schwardt NH, Barber CC, Riley R, et al. Bacterial endosymbionts influence host sexuality and reveal reproductive genes of early divergent fungi. Nat Commun. 2017;8:1843.

    Article  PubMed  PubMed Central  Google Scholar 

  35. Shaffer JP, Carter ME, Spraker JE, Clark M, Smith BA, Hockett KL, et al. Transcriptional profiles of a foliar fungal endophyte (Pestalotiopsis, Ascomycota) and its bacterial symbiont (Luteibacter, Gammaproteobacteria) reveal sulfur exchange and growth regulation during early phases of symbiotic interaction. mSystems. 2022;7:e00091–22.

  36. Vannini C, Carpentieri A, Salvioli A, Novero M, Marsoni M, Testa L, et al. An interdomain network: the endobacterium of a mycorrhizal fungus promotes antioxidative responses in both fungal and plant hosts. New Phytol. 2016;211:265–75.

    Article  CAS  PubMed  Google Scholar 

  37. Lastovetsky OA, Gaspar ML, Mondo SJ, LaButti KM, Sandor L, Grigoriev IV, et al. Lipid metabolic changes in an early divergent fungus govern the establishment of a mutualistic symbiosis with endobacteria. Proc Natl Acad Sci U S A. 2016;113:15102–7.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  38. Moebius N, Üzüm Z, Dijksterhuis J, Lackner G, Hertweck C. Active invasion of bacteria into living fungal cells. Nürnberger T, editor. eLife. 2014;3:e03007.

  39. Dearth SP, Castro HF, Venice F, Tague ED, Novero M, Bonfante P, et al. Metabolome changes are induced in the arbuscular mycorrhizal fungus Gigaspora margarita by germination and by its bacterial endosymbiont. Mycorrhiza. 2018;28:421–33.

    Article  CAS  PubMed  Google Scholar 

  40. Muszewska A, Okrasińska A, Steczkiewicz K, Drgas O, Orłowska M, Perlińska-Lenart U, et al. Metabolic potential, ecology and presence of associated bacteria is reflected in genomic diversity of Mucoromycotina. Frontiers in Microbiology. 2021;12. Available from: https://www.frontiersin.org/articles/https://doi.org/10.3389/fmicb.2021.636986. [cited 2023 Feb 27].

  41. Partida-Martinez LP, Hertweck C. Pathogenic fungus harbours endosymbiotic bacteria for toxin production. Nature. 2005;437:884–8.

    Article  CAS  PubMed  Google Scholar 

  42. Salvioli A, Chiapello M, Fontaine J, Hadj-Sahraoui AL, Grandmougin-Ferjani A, Lanfranco L, et al. Endobacteria affect the metabolic profile of their host Gigaspora margarita, an arbuscular mycorrhizal fungus. Environ Microbiol. 2010;12:2083–95.

    Article  CAS  PubMed  Google Scholar 

  43. Li Z, Yao Q, Dearth SP, Entler MR, Castro Gonzalez HF, Uehling JK, et al. Integrated proteomics and metabolomics suggests symbiotic metabolism and multimodal regulation in a fungal-endobacterial system. Environ Microbiol. 2017;19:1041–53.

    Article  CAS  PubMed  Google Scholar 

  44. Misztal PK, Lymperopoulou DS, Adams RI, Scott RA, Lindow SE, Bruns T, et al. Emission factors of microbial volatile organic compounds from environmental bacteria and fungi. Environ Sci Technol. 2018;52:8272–82.

    Article  CAS  PubMed  Google Scholar 

  45. Malar CM, Krüger M, Krüger C, Wang Y, Stajich JE, Keller J, et al. The genome of Geosiphon pyriformis reveals ancestral traits linked to the emergence of the arbuscular mycorrhizal symbiosis. Curr Biol. 2021;31:1570–1577.e4.

    Article  Google Scholar 

  46. Siozios S, Tosi L, Ferrarini A, Ferrari A, Tononi P, Bellin D, et al. Transcriptional reprogramming of the mycoparasitic fungus Ampelomyces quisqualis during the powdery mildew host-induced germination. Phytopathology®. 2015;105:199–209.

  47. Gilbert KB, Holcomb EE, Allscheid RL, Carrington JC. Hiding in plain sight: new virus genomes discovered via a systematic analysis of fungal public transcriptomes. PLoS ONE. 2019;14:e0219207.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  48. Myers JM, Bonds AE, Clemons RA, Thapa NA, Simmons DR, Carter-House D, et al. Survey of early-diverging lineages of fungi reveals abundant and diverse mycoviruses. mBio. 2020;11:e02027–20.

  49. Turina M, Ghignone S, Astolfi N, Silvestri A, Bonfante P, Lanfranco L. The virome of the arbuscular mycorrhizal fungus Gigaspora margarita reveals the first report of DNA fragments corresponding to replicating non-retroviral RNA viruses in fungi. Environ Microbiol. 2018;20:2012–25.

    Article  CAS  PubMed  Google Scholar 

  50. Lee K-M, Cho WK, Yu J, Son M, Choi H, Min K, et al. A comparison of transcriptional patterns and mycological phenotypes following infection of Fusarium graminearum by four mycoviruses. PLoS ONE. 2014;9:e100989.

    Article  PubMed  PubMed Central  Google Scholar 

  51. Zheng L, Shu C, Zhang M, Yang M, Zhou E. Molecular characterization of a novel endornavirus conferring hypovirulence in rice sheath blight fungus Rhizoctonia solani AG-1 IA strain GD-2. Viruses. 2019;11:178.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  52. Kwon S-J, Cho S-Y, Lee K-M, Yu J, Son M, Kim K-H. Proteomic analysis of fungal host factors differentially expressed by Fusarium graminearum infected with Fusarium graminearum virus-DK21. Virus Res. 2009;144:96–106.

    Article  CAS  PubMed  Google Scholar 

  53. Shi X-X, Qiu H-P, Wang J, Zhang Z, Wang Y-L, Sun G-C. A handy method to remove bacterial contamination from fungal cultures. PLoS ONE. 2019;14:e0224635.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  54. Baltrus DA, Dougherty K, Arendt KR, Huntemann M, Clum A, Pillay M, et al. Absence of genome reduction in diverse, facultative endohyphal bacteria. Microb Genom. 2017;3:e000101.

    PubMed  PubMed Central  Google Scholar 

  55. Uehling JK, Entler MR, Meredith HR, Millet LJ, Timm CM, Aufrecht JA, et al. Microfluidics and metabolomics reveal symbiotic bacterial–fungal interactions between Mortierella elongata and Burkholderia include metabolite exchange. Frontiers in Microbiology. 2019;10. Available from: https://www.frontiersin.org/articles/https://doi.org/10.3389/fmicb.2019.02163. [cited 2022 Jul 28].

  56. Ohshima S, Sato Y, Fujimura R, Takashima Y, Hamada M, Nishizawa T, et al. Mycoavidus cysteinexigens gen. nov., sp. nov., an endohyphal bacterium isolated from a soil isolate of the fungus Mortierella elongata. International Journal of Systematic and Evolutionary Microbiology. 2016;66:2052–7.

  57. Fiedorová K, Radvanský M, Němcová E, Grombiříková H, Bosák J, Černochová M, et al. The impact of DNA extraction methods on stool bacterial and fungal microbiota community recovery. Frontiers in Microbiology. 2019;10. Available from: https://www.frontiersin.org/article/https://doi.org/10.3389/fmicb.2019.00821. [cited 2022 Apr 27].

  58. Vesty A, Biswas K, Taylor MW, Gear K, Douglas RG. Evaluating the impact of DNA extraction method on the representation of human oral bacterial and fungal communities. PLoS ONE. 2017;12:e0169877.

    Article  PubMed  PubMed Central  Google Scholar 

  59. McCarthy A, Chiang E, Schmidt ML, Denef VJ. RNA preservation agents and nucleic acid extraction method bias perceived bacterial community composition. PLoS ONE. 2015;10:e0121659.

    Article  PubMed  PubMed Central  Google Scholar 

  60. Lee M-K, Park H-S, Han K-H, Hong S-B, Yu J-H. High molecular weight genomic DNA mini-prep for filamentous fungi. Fungal Genet Biol. 2017;104:1–5.

    Article  CAS  PubMed  Google Scholar 

  61. Petersen C, Sørensen T, Westphal KR, Fechete LI, Sondergaard TE, Sørensen JL, et al. High molecular weight DNA extraction methods lead to high quality filamentous ascomycete fungal genome assemblies using Oxford Nanopore sequencing. Microbial Genomics. 8:000816.

  62. Duportet X, Aggio RBM, Carneiro S, Villas-Bôas SG. The biological interpretation of metabolomic data can be misled by the extraction method used. Metabolomics. 2012;8:410–21.

    Article  CAS  Google Scholar 

  63. Chen F, Ma R, Chen X-L. Advances of metabolomics in fungal pathogen–plant interactions. Metabolites. 2019;9:169.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  64. Zhang X, Li L, Mayne J, Ning Z, Stintzi A, Figeys D. Assessing the impact of protein extraction methods for human gut metaproteomics. J Proteomics. 2018;180:120–7.

    Article  CAS  PubMed  Google Scholar 

  65. Bianco L, Perrotta G. Methodologies and perspectives of proteomics applied to filamentous fungi: from sample preparation to secretome analysis. Int J Mol Sci. 2015;16:5803–29.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  66. Nakayasu ES, Nicora CD, Sims AC, Burnum-Johnson KE, Kim Y-M, Kyle JE, et al. MPLEx: a robust and universal protocol for single-sample integrative proteomic, metabolomic, and lipidomic analyses. mSystems. 2016;1:e00043–16.

    Article  PubMed  PubMed Central  Google Scholar 

  67. Kang J, David L, Li Y, Cang J, Chen S. Three-in-one simultaneous extraction of proteins, metabolites and lipids for multi-omics. Front Genet. 2021;12: 635971.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  68. Salem MA, Jüppner J, Bajdzienko K, Giavalisco P. Protocol: a fast, comprehensive and reproducible one-step extraction method for the rapid preparation of polar and semi-polar metabolites, lipids, proteins, starch and cell wall polymers from a single sample. Plant Methods. 2016;12:45.

    Article  PubMed  PubMed Central  Google Scholar 

  69. Nicora CD, Burnum-Johnson KE, Nakayasu ES, Casey CP, Iii RAW, Chowdhury TR, et al. The MPLEx protocol for multi-omic analyses of soil samples. JoVE (Journal of Visualized Experiments). 2018;e57343.

  70. Coman C, Solari FA, Hentschel A, Sickmann A, Zahedi RP, Ahrends R. Simultaneous metabolite, protein, lipid extraction (SIMPLEX): a combinatorial multimolecular omics approach for systems biology *. Mol Cell Proteomics. 2016;15:1435–66.

    Article  Google Scholar 

  71. Paolillo C, Londin E, Fortina P. Single-cell genomics. Clin Chem. 2019;65:972–85.

    Article  CAS  PubMed  Google Scholar 

  72. Gasc C, Peyret P. Hybridization capture reveals microbial diversity missed using current profiling methods. Microbiome. 2018;6:61.

    Article  PubMed  PubMed Central  Google Scholar 

  73. O’Neil D, Glowatz H, Schlumpberger M. Ribosomal RNA depletion for efficient use of RNA-Seq capacity. Current Protocols in Molecular Biology. 2013;103:4.19.1–4.19.8.

  74. Chung M, Bruno VM, Rasko DA, Cuomo CA, Muñoz JF, Livny J, et al. Best practices on the differential expression analysis of multi-species RNA-seq. Genome Biol. 2021;22:121.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  75. Telzrow CL, Zwack PJ, Esher Righi S, Dietrich FS, Chan C, Owzar K, et al. Comparative analysis of RNA enrichment methods for preparation of Cryptococcus neoformans RNA sequencing libraries. G3 Genes|Genomes|Genetics. 2021;11:jkab301.

  76. Wahl A, Huptas C, Neuhaus K. Comparison of rRNA depletion methods for efficient bacterial mRNA sequencing. Sci Rep. 2022;12:1–11.

    Article  Google Scholar 

  77. Bhargava V, Head SR, Ordoukhanian P, Mercola M, Subramaniam S. Technical variations in low-input RNA-seq methodologies. Sci Rep. 2014;4:3678.

    Article  PubMed  PubMed Central  Google Scholar 

  78. Grün D, van Oudenaarden A. Design and analysis of single-cell sequencing experiments. Cell. 2015;163:799–810.

    Article  PubMed  Google Scholar 

  79. Kashima Y, Sakamoto Y, Kaneko K, Seki M, Suzuki Y, Suzuki A. Single-cell sequencing techniques from individual to multiomics analyses. Exp Mol Med. 2020;52:1419–27.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  80. Zhang W, Ramautar R. CE-MS for metabolomics: developments and applications in the period 2018–2020. Electrophoresis. 2021;42:381–401.

    Article  CAS  PubMed  Google Scholar 

  81. Zhang Z, Qu Y, Dovichi NJ. Capillary zone electrophoresis-mass spectrometry for bottom-up proteomics. TrAC, Trends Anal Chem. 2018;108:23–37.

    Article  Google Scholar 

  82. Wu R, Xing S, Badv M, Didar TF, Lu Y. Step-wise assessment and optimization of sample handling recovery yield for nanoproteomic analysis of 1000 mammalian cells. Anal Chem. 2019;91:10395–400.

    Article  CAS  PubMed  Google Scholar 

  83. Kassem S, van der Pan K, de Jager AL, Naber BAE, de Laat IF, Louis A, et al. Proteomics for low cell numbers: how to optimize the sample preparation workflow for mass spectrometry analysis. J Proteome Res. 2021;20:4217–30.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  84. Kraut A, Marcellin M, Adrait A, Kuhn L, Louwagie M, Kieffer-Jaquinod S, et al. Peptide storage: are you getting the best return on your investment? Defining optimal storage conditions for proteomics samples. J Proteome Res. 2009;8:3778–85.

    Article  CAS  PubMed  Google Scholar 

  85. Xu G, Zhang L, Liu X, Guan F, Xu Y, Yue H, et al. Combined assembly of long and short sequencing reads improve the efficiency of exploring the soil metagenome. BMC Genomics. 2022;23:37.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  86. Gehrig JL, Portik DM, Driscoll MD, Jackson E, Chakraborty S, Gratalo D, et al. Finding the right fit: evaluation of short-read and long-read sequencing approaches to maximize the utility of clinical microbiome data. Microb Genom. 2022;8:000794.

    PubMed  PubMed Central  Google Scholar 

  87. Zamboni N, Saghatelian A, Patti GJ. Defining the metabolome: size, flux, and regulation. Mol Cell. 2015;58:699–706.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  88. Duncan KD, Fyrestam J, Lanekoff I. Advances in mass spectrometry based single-cell metabolomics. Analyst. 2019;144:782–93.

    Article  CAS  PubMed  Google Scholar 

  89. Ali A, Abouleila Y, Shimizu Y, Hiyama E, Emara S, Mashaghi A, et al. Single-cell metabolomics by mass spectrometry: advances, challenges, and future applications. TrAC, Trends Anal Chem. 2019;120:115436.

    Article  Google Scholar 

  90. Slavov N. Scaling up single-cell proteomics. Molecular & cellular proteomics. 2022;21. Available from: https://www.mcponline.org/article/S1535-9476(21)00151-1/abstract. [cited 2023 Mar 1].

  91. Ibáñez AJ, Fagerer SR, Schmidt AM, Urban PL, Jefimovs K, Geiger P, et al. Mass spectrometry-based metabolomics of single yeast cells. Proc Natl Acad Sci. 2013;110:8790–4.

    Article  PubMed  PubMed Central  Google Scholar 

  92. Guo S, Zhang C, Le A. The limitless applications of single-cell metabolomics. Curr Opin Biotechnol. 2021;71:115–22.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  93. Ong T-H, Kissick DJ, Jansson ET, Comi TJ, Romanova EV, Rubakhin SS, et al. Classification of large cellular populations and discovery of rare cells using single cell matrix-assisted laser desorption/ionization time-of-flight mass spectrometry. Anal Chem. 2015;87:7036–42.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  94. Lapainis T, Rubakhin SS, Sweedler JV. Capillary electrophoresis with electrospray ionization mass spectrometric detection for single-cell metabolomics. Anal Chem. 2009;81:5858–64.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  95. Kawai T, Ota N, Okada K, Imasato A, Owa Y, Morita M, et al. Ultrasensitive single cell metabolomics by capillary electrophoresis–mass spectrometry with a thin-walled tapered emitter and large-volume dual sample preconcentration. Anal Chem. 2019;91:10564–72.

    Article  CAS  PubMed  Google Scholar 

  96. Okrasińska A, Bokus A, Duk K, Gęsiorska A, Sokołowska B, Miłobędzka A, et al. New endohyphal relationships between Mucoromycota and Burkholderiaceae representatives. Appl Environ Microbiol. 2021;87:e02707–e2720.

    Article  PubMed  PubMed Central  Google Scholar 

  97. Morales DP, Robinson AJ, Pawlowski AC, Ark C, Kelliher JM, Junier P, et al. Advances and challenges in fluorescence in situ hybridization for visualizing fungal endobacteria. Front Microbiol. 2022;13:892227.

    Article  PubMed  PubMed Central  Google Scholar 

  98. Safieddine A, Coleno E, Lionneton F, Traboulsi A-M, Salloum S, Lecellier C-H, et al. HT-smFISH: a cost-effective and flexible workflow for high-throughput single-molecule RNA imaging. Nat Protoc. 2023;18:157–87.

    Article  CAS  PubMed  Google Scholar 

  99. Desirò A, Naumann M, Epis S, Novero M, Bandi C, Genre A, et al. Mollicutes-related endobacteria thrive inside liverwort-associated arbuscular mycorrhizal fungi. Environ Microbiol. 2013;15:822–36.

    Article  PubMed  Google Scholar 

  100. Gundlach KA, Briegel A. Zooming in on host-symbiont interactions: advances in cryo-EM sample processing methods and future application to symbiotic tissues. Symbiosis. 2022;87:67–75.

    Article  Google Scholar 

  101. Albornoz FE, Hayes PE, Orchard S, Clode PL, Nazeri NK, Standish RJ, et al. First cryo-scanning electron microscopy images and X-ray microanalyses of Mucoromycotinian fine root endophytes in vascular plants. Frontiers in Microbiology. 2020;11. Available from: https://www.frontiersin.org/articles/https://doi.org/10.3389/fmicb.2020.02018. [cited 2023 Jul 18].

  102. Dueñas ME, Essner JJ, Lee YJ. 3D MALDI mass spectrometry imaging of a single cell: spatial mapping of lipids in the embryonic development of zebrafish. Sci Rep. 2017;7:14946.

    Article  PubMed  PubMed Central  Google Scholar 

  103. Rappez L, Stadler M, Triana S, Gathungu RM, Ovchinnikova K, Phapale P, et al. SpaceM reveals metabolic states of single cells. Nat Methods. 2021;18:799–805.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  104. Nayfach S, Camargo AP, Schulz F, Eloe-Fadrosh E, Roux S, Kyrpides NC. CheckV assesses the quality and completeness of metagenome-assembled viral genomes. Nat Biotechnol. 2021;39:578–85.

    Article  CAS  PubMed  Google Scholar 

  105. Parks DH, Imelfort M, Skennerton CT, Hugenholtz P, Tyson GW. CheckM: assessing the quality of microbial genomes recovered from isolates, single cells, and metagenomes. Genome Res. 2015;25:1043–55.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  106. Belliardo C, Koutsovoulos GD, Rancurel C, Clément M, Lipuma J, Bailly-Bechet M, et al. Improvement of eukaryotic protein predictions from soil metagenomes. Sci Data. 2022;9:311.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  107. Lind AL, Pollard KS. Accurate and sensitive detection of microbial eukaryotes from whole metagenome shotgun sequencing. Microbiome. 2021;9:58.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  108. Camargo A. geNomad. 2023. Available from: https://github.com/apcamargo/genomad. [cited 2023 Feb 20].

  109. Eun Kang J, Ciampi A, Hijri M. SeSaMe: metagenome sequence classification of arbuscular mycorrhizal fungi-associated microorganisms. Genomics Proteomics Bioinformatics. 2020;18:601–12.

    Article  PubMed  PubMed Central  Google Scholar 

  110. Bankevich A, Nurk S, Antipov D, Gurevich AA, Dvorkin M, Kulikov AS, et al. SPAdes: a new genome assembly algorithm and its applications to single-cell sequencing. J Comput Biol. 2012;19:455–77.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  111. Gatter T, von Löhneysen S, Fallmann J, Drozdova P, Hartmann T, Stadler PF. LazyB: fast and cheap genome assembly. Algorithms for Molecular Biology. 2021;16:8.

    Article  PubMed  PubMed Central  Google Scholar 

  112. Chang Y, Desirò A, Na H, Sandor L, Lipzen A, Clum A, et al. Phylogenomics of Endogonaceae and evolution of mycorrhizas within Mucoromycota. New Phytol. 2019;222:511–25.

    Article  PubMed  Google Scholar 

  113. Grabherr MG, Haas BJ, Yassour M, Levin JZ, Thompson DA, Amit I, et al. Full-length transcriptome assembly from RNA-Seq data without a reference genome. Nat Biotechnol. 2011;29:644–52.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  114. Gonzalez E, Pitre FE, Pagé AP, Marleau J, Guidi Nissim W, St-Arnaud M, et al. Trees, fungi and bacteria: tripartite metatranscriptomics of a root microbiome responding to soil contamination. Microbiome. 2018;6:53.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  115. Hölzer M, Marz M. De novo transcriptome assembly: a comprehensive cross-species comparison of short-read RNA-Seq assemblers. GigaScience. 2019;8:giz039.

    Article  PubMed  PubMed Central  Google Scholar 

  116. Shu L, Wang M, Xu H, Qiu Z, Li T. De novo transcriptome assembly and comprehensive assessment provide insight into fruiting body formation of Sparassis latifolia. Sci Rep. 2022;12:11075.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  117. Meng A, Marchet C, Corre E, Peterlongo P, Alberti A, Da Silva C, et al. A de novo approach to disentangle partner identity and function in holobiont systems. Microbiome. 2018;6:105.

    Article  PubMed  PubMed Central  Google Scholar 

  118. Shakya M, Lo C-C, Chain PSG. Advances and challenges in metatranscriptomic analysis. Frontiers in Genetics. 2019;10. Available from: https://www.frontiersin.org/articles/https://doi.org/10.3389/fgene.2019.00904. [cited 2022 Jul 19].

  119. Muth T, Hartkopf F, Vaudel M, Renard BY. A potential golden age to come—current tools, recent use cases, and future avenues for de novo sequencing in proteomics. Proteomics. 2018;18:1700150.

    Article  Google Scholar 

  120. Muth T, Renard BY. Evaluating de novo sequencing in proteomics: already an accurate alternative to database-driven peptide identification? Brief Bioinform. 2018;19:954–70.

    Article  CAS  PubMed  Google Scholar 

  121. Chaleckis R, Meister I, Zhang P, Wheelock CE. Challenges, progress and promises of metabolite annotation for LC–MS-based metabolomics. Curr Opin Biotechnol. 2019;55:44–50.

    Article  CAS  PubMed  Google Scholar 

  122. Li G, Jian T, Liu X, Lv Q, Zhang G, Ling J. Application of metabolomics in fungal research. Molecules. 2022;27:7365.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  123. Cambiaghi A, Ferrario M, Masseroli M. Analysis of metabolomic data: tools, current strategies and future challenges for omics data integration. Brief Bioinform. 2017;18:498–510.

    CAS  PubMed  Google Scholar 

  124. Subramanian I, Verma S, Kumar S, Jere A, Anamika K. Multi-omics data integration, interpretation, and its application: bioinformatics and biology insights. 2020; Available from: https://journals.sagepub.com/doi/full/https://doi.org/10.1177/1177932219899051. [cited 2022 Jul 27].

  125. Hu B, Canon S, Eloe-Fadrosh EA, Anubhav, Babinski M, Corilo Y, et al. Challenges in bioinformatics workflows for processing microbiome omics data at scale. Frontiers in Bioinformatics. 2022;1. Available from: https://www.osti.gov/pages/biblio/1840628-challenges-bioinformatics-workflows-processing-microbiome-omics-data-scale. [cited 2022 Apr 8].

Download references

Acknowledgements

Not applicable

Funding

This research was supported by a Science Focus Area Grant from the US Department of Energy (DOE), Biological and Environmental Research (BER), and Biological System Science Division (BSSD) under the grant number LANLF59T.

Author information

Authors and Affiliations

Authors

Contributions

JK conceptualized and wrote the majority of the manuscript. AR and PC were major contributors in writing the manuscript. RL, LJ, BH, DM, PJ, and GB each contributed substantially to the manuscript sections most relevant to their areas of expertise. GC significantly contributed to the figures and overall organization of the manuscript and its visualizations.

Corresponding authors

Correspondence to Julia M. Kelliher or Patrick S. G. Chain.

Ethics declarations

Ethics approval and consent to participate

Not applicable.

Consent for publication

Not applicable.

Competing interests

The authors declare no competing interests.

Additional information

Publisher’s Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Rights and permissions

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated in a credit line to the data.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Kelliher, J.M., Robinson, A.J., Longley, R. et al. The endohyphal microbiome: current progress and challenges for scaling down integrative multi-omic microbiome research. Microbiome 11, 192 (2023). https://doi.org/10.1186/s40168-023-01634-7

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s40168-023-01634-7

Keywords