{"paper_id":"caa8eda4-b320-4298-a442-ab002bb3fcb6","body_text":"DNA and RNA metabarcoding reveal distinct seed-borne mycobiota | Research Square window.SnipcartSettings = { analytics: { enabled: false } }; (function() { var accessVector = localStorage.getItem('access_vector') || ''; window.dataLayer = window.dataLayer || []; if (accessVector) { window.dataLayer.push({ user: { profile: { profileInfo: { snid: accessVector } } } }); } })(); (function(w,d,s,l,i){w[l]=w[l]||[];w[l].push({'gtm.start':new Date().getTime(),event:'gtm.js'});var f=d.getElementsByTagName(s)[0],j=d.createElement(s),dl=l!='dataLayer'?'&l='+l:'';j.async=true;j.src='https://www.googletagmanager.com/gtm.js?id='+i+dl;f.parentNode.insertBefore(j,f);})(window,document,'script','dataLayer','GTM-K279D39R'); Browse Preprints In Review Journals COVID-19 Preprints AJE Video Bytes Research Tools Research Promotion AJE Professional Editing AJE Rubriq About Preprint Platform In Review Editorial Policies Our Team Advisory Board Help Center Sign In Submit a Preprint Cite Share Download PDF Research Article DNA and RNA metabarcoding reveal distinct seed-borne mycobiota Iva Franić, Patrick Sherwood, Kinga Stolarek, René Eschen, Jana Orbach, and 2 more This is a preprint; it has not been peer reviewed by a journal. https://doi.org/ 10.21203/rs.3.rs-8095745/v1 This work is licensed under a CC BY 4.0 License Status: Under Review Version 1 posted 13 You are reading this latest preprint version Abstract Background Tree seeds harbor diverse fungal communities, including both pathogens and mutualists, that can influence plant health. These communities comprise living, metabolically active organisms as well as dormant or dead cells. Because only active fungi interact with their hosts, distinguishing active from inactive taxa is crucial, especially for environmental and phytosanitary monitoring. Traditional culturing methods capture living fungi but account for only a small fraction of the total fungal diversity. Currently, these methods are increasingly replaced by high-throughput DNA metabarcoding, which detects a broader range of taxa. However, DNA persists after cell death and occurs in dormant cells, preventing distinction between active and inactive fungi. In contrast, RNA metabarcoding detects metabolically active organisms and may better reflect living fungal communities than the other two methods, though its use in assessing plant-associated fungi remains underexplored. We used culturing, DNA-, and RNA-based metabarcoding to compare fungal communities associated with seeds of three key European tree species ( Fagus sylvatica , Abies alba , Pinus sylvestris ). Results DNA and RNA metabarcoding detected largely distinct, non-overlapping fungal communities, with differences primarily driven by rare active taxa in the RNA dataset. Several cultured genera—likely representing abundant and metabolically active taxa—were shared between both metabarcoding approaches. Conclusions These results highlight the complementary nature of the three methods for characterising seed-associated fungi. Combining culturing, DNA- and RNA-based metabarcoding may provide the most comprehensive assessment of fungal diversity, while RNA metabarcoding alone offers a promising opportunity to identify the active members of fungal communities for improved environmental and phytosanitary monitoring. microbiome high-throughput sequencing viable community metabolic activity relic DNA phytosanitary risk Figures Figure 1 Figure 2 Figure 3 Figure 4 Figure 5 Background Seeds are essential for plant reproduction and dispersal and act as selective microbial habitats. Packed with oils, starches, and sugars, their resources are protected by seed coats and chemical defences [ 1 ], allowing access only to specialized microorganisms [ 2 ]. Molecular studies show that tree seeds host diverse fungi [ 3 ], including plant pathogens [ 4 – 6 ], and fungi that protect plants against biotic and abiotic stresses [ 7 ]. However, few of these fungi may be metabolically active, while the remainder are dead or dormant. Differentiating these fractions is key to understanding the roles of seed-associated fungi in their hosts and their impact on tree health. Plant-associated fungi have been traditionally studied by growing them from plant tissues on nutrient media and identifying the isolates based on morphology and/or with molecular techniques [ 8 ]. This approach has been widely used to explore fungal biodiversity and ecology [ 9 – 11 ] and remains essential for obtaining living cultures, describing new species and conducting phenotyping assays (e.g., host range, virulence, optimal growth temperatures). For studying new host-pathogen associations, living cultures are also necessary to fulfil Koch’s postulates. However, culturing captures only a subset of fungal diversity—mainly abundant fungi—excluding taxa that are difficult or impossible to isolate [ 3 ], but which may affect tree health [ 12 ]. Moreover, culturing is time-consuming and costly [ 13 ], limiting its usefulness for large-scale applications such as environmental monitoring or phytosanitary screening. Recent advances in high-throughput sequencing (HTS) have enabled the rapid and sensitive analysis of many samples simultaneously. Amplicon-based HTS (i.e., metabarcoding) has become widely used for profiling fungal communities from various substrates, including plant tissues, to address questions related to their diversity and biosecurity [ 14 , 15 ]. This approach typically targets the nuclear ribosomal DNA (rDNA) internal transcribed spacer (ITS) region—a universal fungal barcode [ 16 ]. Sequencing platforms such as Illumina or PacBio are then used to characterise fungi from DNA extracted directly from bulk, host or environmental samples. However, DNA metabarcoding cannot distinguish living, metabolically active from dead or dormant organisms [ 14 ]. For example, Carini et al. [ 17 ] showed that over 40% of ITS sequences in soil belong to dead organisms, inflating diversity and skewing community profiles. This raises concerns when using DNA metabarcoding to study fungal responses to environmental change or in biosecurity [ 17 , 18 ], where only viable organisms pose a biosecurity threat that requires mitigation measures. Additionally, distinguishing dead from living organisms is essential for evaluating the effectiveness of phytosanitary treatments. Unlike DNA, RNA exists only in metabolically active cells and some viruses, and degrades rapidly after cell death [ 19 , 20 ]. In fungi, the rDNA operon is transcribed into a precursor ribosomal RNA (rRNA) that includes the ITS region. Targeting this precursor allows identification of fungi actively transcribing rRNA [ 21 ], resulting in a more accurate characterisation of the living community. While the potential of using RNA for improved biological monitoring has been recognised in environmental fields such as water research [ 22 ], the approach remains largely underexplored, including for detecting active fungi in plant samples. Studies using DNA and RNA metabarcoding revealed different total and active microbial communities in temperate soil, wood [ 23 , 24 ], ballast water [ 18 , 25 ] and ocean [ 26 ]. However, other studies found no differences. For example, 92% of fungal taxa in Picea abies needles overlapped between DNA and RNA datasets [ 27 ], and similar patterns were observed in Arctic soil [ 28 ]. The differences between total and active fungal communities may depend on environmental factors such as resource availability and temperature, which influence fungal metabolism [ 29 , 30 ] and DNA/RNA turnover [ 31 ]. Methodological biases may also contribute to the observed discrepancies, highlighting the need to assess how well these two approaches reflect fungal assemblages—especially when distinguishing active fungi from dead or dormant ones is critical for ecological inference and biosecurity. In this study, we characterised fungal communities associated with commercially traded seeds of three ecologically and economically important European tree species—common beech ( Fagus sylvatica L.), silver fir ( Abies alba Mill.), and Scots pine ( Pinus sylvestris L.)—using traditional culturing and DNA- and RNA-based metabarcoding. We assessed incidence-based and abundance-weighted alpha diversity (i.e., OTU richness and Inverse Simpson index, respectively) and beta diversity (i.e., Sørensen and Morisita-Horn dissimilarities, respectively) across tree species and molecular methods. We also compared taxonomic and lifestyle composition of fungal communities assessed with DNA and RNA metabarcoding within each tree species and compared fungal communities detected by each metabarcoding method with those from culturing. We hypothesised that DNA metabarcoding would reveal greater diversity than RNA metabarcoding, as RNA-based communities represent the active subset of the DNA-detected communities. Moreover, we hypothesised that fungi obtained by culturing would more closely resemble RNA-based communities than DNA-based ones, because culturing primarily recovers metabolically active fungi at the time of isolation. Methods Study material A total of 16 seed lots (i.e., batches of seeds from specific locations and tree species; F. sylvatica (n = 6), A. alba (n = 5), P. sylvestris (n = 5)) were obtained in winter 2019 from European commercial seed suppliers. Seeds were stored at 4°C until fungal assessment using culturing, DNA- and RNA-metabarcoding in spring 2020. Fungal assessment by culturing Seeds were surface sterilized by subsequent immersion in 70% ethanol (1 min), 1% sodium hypochlorite (5 min), 70% ethanol (30 sec), and sterile water (30 sec) [ 32 ], then dried in a laminar flow hood. For fungal assessment, 90 seeds per seed lot were crushed with a sterile pestle and placed on nutrient media containing streptomycin (0.05 mg/mL). To capture a broad range of fungi, three media with varying nutrient content were used—water agar (WA 15 g/L; VWR Chemicals, Solon, Ohio, USA), malt extract agar (MEA; ME 30g/L; Duchefa Biochemie, Haarlem, the Netherlands & WA 15 g/L; VWR Chemicals) and potato dextrose agar (PDA 39g/L; Merck KGaA, Darmstadt, Germany). Thirty seeds were plated per tree species x medium combination. Plates were checked every five days for 30 days. Emerging colonies were transferred to PDA (39 g/L; Merck KGaA) to obtain pure cultures. Isolates were grouped based on macromorphological traits (i.e., colour, texture, form, margin). One isolate per morphotype was selected for DNA extraction. Mycelia were scraped from agar, freeze-dried overnight, and DNA was extracted using the E.Z.N.A.® SP Plant DNA Kit (Omega Bio-Tek, Norcross, Georgia, USA) following the manufacturer’s instructions. DNA concentrations were measured using a DS-11 UV-Vis Spectrophotometer (DeNovix, Wilmington, Delaware, USA) and diluted to 5 ng/µl. DNA extracts with concentrations below 5 ng/µl were not diluted. The ITS region was amplified in 25 µl reactions, containing 2 µl DNA, 8.5 µl nuclease-free H 2 O, 12.5 µl DreamTaqPCR Master Mix (2X), and 1 µl each of ITS1 (F: TCCGTAGGTGAACCTGCGG) and ITS4 (R: TCCTCCGCTTATTGATATGC) primers [ 33 ]. PCRs were run on an Eppendorf Mastercycler (Eppendorf, Hamburg, Germany) under the following conditions: 95°C for 2 min; 35 cycles of 95°C for 1 min; 55°C for 45 sec; 72°C for 1.5 min; final extension at 72°C for 5 min. Amplification success was confirmed by gel electrophoresis. PCR products were purified and Sanger-sequenced by Macrogen-Europe (Amsterdam, the Netherlands) using the same primers as in PCRs. Forward and reverse sequences were trimmed and assembled into consensus sequences using CLC Main Workbench 26 (Qiagen, Aarhus, Denmark). Study sequences were compared against the National Centre for Biotechnology Information (NCBI) core nucleotide (nt) database using MEGABLAST [ 34 , 35 ], excluding sequences annotated as uncultured or environmental samples, to ensure reliable taxonomic assignment. Fungal assessment by metabarcoding Fungi were also assessed from 30 seeds per seed lot using DNA and RNA metabarcoding. Seeds were surface sterilised as for culturing, then ground in liquid nitrogen using RNase AWAY ™ Surface Decontaminant-treated mortars and pestles (Thermo Fisher Scientific, Waltham, USA). DNA was extracted using the E.Z.N.A.® SP Plant DNA Kit (Omega Bio-Tek) following the manufacturer’s instructions. DNA quantity was determined with a Qubit dsDNA BR Assay kit on a Qubit 3.0 Fluorometer (Thermo Fisher Scientific). Extracts were sent to Novogene (Cambridge, UK) for PCR, library preparation and sequencing. The ITS2 region was amplified using primers ITS3 (F: GCATCGATGAAGAACGCAGC) and ITS4 (R: TCCTCCGCTTATTGATATFC) [ 33 ], and libraries were sequenced on an Illumina MiSeq (2 x 250 bp). RNA was extracted using the E.Z.N.A.® Plant RNA kit (Omega Bio-Tek) following the manufacturer`s protocol for difficult samples, from the same ground seed pool used for DNA extractions. Multiple extractions per sample were pooled to obtain sufficient RNA for complementary DNA (i.e., cDNA) synthesis. After pooling, RNA quantity was measured using the Qubit™ RNA Broad Range (BR) Assay Kit with Qubit 3.0 Fluorometer (Thermo Fisher Scientific), and quality was assessed using the Agilent RNA 6000 Pico Kit and Agilent 2100 Bioanalyzer (Agilent, Santa Clara, California, USA). DNA contamination was removed using TURBO DNase (Invitrogen™, Thermo Fisher Scientific) with a 60 min incubation (4 U/reaction), confirmed by 35-cycle ITS PCR (ITS1/ITS4 primers; [ 33 ]) as described above and the absence of bands in 1.2% TAE agarose gel electrophoresis. RNA was cleaned and concentrated with the RNeasy MinElute Cleanup Kit (QIAGEN, Hilden, Germany), and cDNA synthesised using the SuperScript™ IV First-Strand System (Invitrogen™, Thermo Fisher Scientific) with random hexamer primers and RNase H treatment. The resulting cDNA was purified with AMPure XP beads (1:4 bead:sample ratio) and stored at − 80°C. Samples were sent to BMKGENE (Münster, Germany) for PCR, library preparation and Illumina NovaSeq 6000 sequencing (2 x 250 bp) of the ITS2 region using ITS3 and ITS4 primers [ 33 ]. Samples were demultiplexed by the sequencing facility, yielding separate forward and reverse FASTQ files for each sample and dataset (i.e., DNA and RNA). Bioinformatic processing was performed using a custom QIIME2-based pipeline [ 36 ]. Adapters and primers were removed with Cutadapt (i.e., error rate 0.1; [ 37 ]). Quality filtering, denoising, read merging and chimera removal were conducted using DADA2 [ 38 ]. Amplicon sequence variants of both DNA and RNA datasets were clustered into OTUs at 97% similarity using VSEARCH [ 39 ], and taxonomic assignment of the merged dataset was performed with a Naïve Bayes classifier [ 40 ] in QIIME2 against the UNITE v8.2 database [ 41 ]. Taxonomic assignments at species level were not considered, because species identification based on the ITS2 region is often unreliable, particularly for certain fungal groups [ 42 ]. Fungal lifestyle assignments Fungal lifestyles were assigned to OTUs based on genus-level matches using the FungalTraits database [ 43 ], referencing the “primary_lifestyle” column. The \"plant pathogen\" category was retained as is, while the following classification adjustments were made: 1) \"ectomycorrhizal\", \"arbuscular_mycorrhizal\" = \"mycorrhiza\", 2) \"foliar_endophyte\", \"root_endophyte\" = \"endophyte\", 3) \"soil_saprotroph\", \"dung_saprotroph\", \"litter_saprotroph\", \"nectar/tap_saprotroph\",\"wood_saprotroph”, “pollen_saprotroph\", \"unspecified_saprotroph\" = \"saprotroph\", 4) \"mycoparasite\",\"animal_parasite\",\"epiphyte\",\"lichen_parasite\", \"algal_parasite”, \"sooty_mold\", \"lichenized\" = \"other\". Statistical analyses All analyses were performed using R Statistical software [ 44 ]. To account for differences in sequencing depth (i.e., number of reads per sample ranged from 7,939 to 1,179,119 with a mean ± se of 48,663 ± 6,590) and retain all the samples, the dataset was rarefied to a minimum number of reads per sample (i.e., 7,939 reads in a DNA sample of P. sylvestris ; Rarefy function, GUniFrac package; [ 45 ]) prior to alpha and beta diversity analyses [ 46 ]. Alpha and beta diversity analyses were conducted using incidence-based and abundance weighted metrices. While incidence-based diversity measures capture overall diversity patterns by giving equal weight to all taxa regardless of abundance, abundance-weighted measures emphasize patterns driven by dominant taxa by assigning greater weight to abundant taxa. Alpha diversity was assessed using OTU richness (i.e., incidence-based metric) and Inverse Simpson Index (i.e., abundance weighted metric; [ 47 ], calculated using the hill_taxa function from the package hillR [ 48 , 49 ]). Alpha diversity metrics (i.e., response variables) were compared between explanatory variables—the methods (i.e., DNA and RNA metabarcoding) and tree species (i.e., F. sylvatica , A. alba and P. sylvestris ). To account for overdispersion, OTU richness was analysed using negative binomial models (i.e., glm.nb function, MASS package; [ 50 ]). Inverse Simpson Index was modelled with a gamma distribution generalized linear model (GLM) using log link (i.e., glm function, stats package; [ 50 ]) due to its non-integer nature. Models with and without the interaction term were compared based on Akaike information criterion (AIC). Models without the interaction term were retained (i.e., lower AIC) and further evaluated for factor significance (i.e., Anova function, car package; [ 51 ]). Estimated marginal means for significant factors and pairwise comparisons were computed with the function emmeans from the emmeans package [ 52 ]. Beta diversity was assessed using Sørensen dissimilarity (i.e., incidence-based metric), and Morisita-Horn dissimilarity (i.e., abundance-weighted metric). Pairwise distances between samples were calculated using the vegdist function from the vegan package [ 53 ]. Overall effects of method and tree species (i.e., explanatory variables) on fungal community composition (i.e., Sørensen and Morisita-Horn dissimilarity; response variables) were tested using Permutational multivariate ANOVA (i.e., PERMANOVA; adonis function, vegan package; [ 53 ]), without interaction term to match alpha diversity models. Differences in fungal community composition between the methods and tree species were visualized using non-metric multidimensional scaling (i.e., NMDS; metaMDS function, vegan package, [ 53 ]). As PERMANOVA revealed significant differences in fungal communities among tree species (see Results), taxonomic composition was assessed separately for each tree species and method, using a non-rarefied dataset to preserve full abundance information across samples [ 54 ]. We focused on reads assigned to fungal genera, enabling inference of ecological roles and host associations. First, the number of reads associated with each distinct fungal genus was calculated for every method and tree species. The genera were then ranked by read count, and the ten most abundant genera per method and tree species were selected for further analysis. All remaining reads were grouped as either “Other” (i.e., reads assigned to genera outside the top ten) or “Unassigned” (i.e., reads not assigned to any genus). Relative abundances of the top ten genera, and genera classified as “Other” or “Unassigned” were plotted across the methods and tree species. In addition, indicator species analysis (i.e., multipatt function, indicspecies package, [ 55 ]) was performed across all genera, irrespective of their read count, to identify those strongly associated with either the DNA or RNA dataset for each tree species. Variation in fungal lifestyles was evaluated using all non-rarefied reads assigned to a specific lifestyle category or classified as “Unassigned”. For each method and tree species, the relative abundance of reads associated with each lifestyle were calculated and plotted. Finally, we identified fungal genera shared between culturing and each metabarcoding dataset, as well as those unique to each dataset. This was first done for all tree species together to assess the overall trend, and for each tree species separately, to assess the differences across tree species. Results Of the 3,167,108 non-singleton reads assigned to 1,797 OTUs in the merged dataset comprising DNA and RNA reads, approximately 50% were fungal (i.e., 1,557,227 reads and 942 OTUs). Rarefaction curves indicated that the sequencing depth was sufficient to capture all fungal OTUs across all samples (Supplementary Fig. 1). More than 60% of the fungal reads originated from A. alba samples, while the remaining reads were similarly distributed between the other two species (Table 1 ). Moreover, the variation in the number of reads per sample was smaller for A. alba than for the other two species (Supplementary Fig. 1). Overall, a similar number of reads was obtained from both sequencing methods, although this varied among tree species. Read counts were comparable between methods in F. sylvatica and A. alba , whereas the P. sylvestris DNA dataset was almost twice the size of the RNA dataset (Table 1 ). Pinus sylvestris samples contained twice as many OTUs as A. alba , with F. sylvatica falling between the two (Table 1 ). The total number of OTUs in all samples of each of the three tree species was three to four times higher in the RNA dataset than in the DNA dataset (Table 1 ). The rarefied dataset that was used for the alpha and beta diversity analyses consisted of 254,048 non-singleton fungal reads equally distributed across all samples, as all samples were rarefied to the same number of reads (i.e., 7,939 reads). Overall, these reads were assigned to 872 OTUs (Table 1 ). Although 70 OTUs (i.e., around 7%) were lost during rarefaction, the number of OTUs in the RNA dataset remained three to four times higher than that in the DNA dataset (Table 1 ). Table 1 Fungal read and OTU counts across three tree species and DNA and RNA metabarcoding datasets . Number of fungal reads and OTUs in samples belonging to the three tree species (i.e., F. sylvatica , A. alba , P. sylvestris ; numbers of seed lots are indicated in brackets) analysed by DNA and RNA metabarcoding are indicated. Numbers are shown for the DNA and RNA datasets and for the combined dataset (i.e., all data) considering all fungal reads and rarefied fungal reads (in the brackets). Reads OTUs Tree species DNA RNA All data DNA RNA All data F. sylvatica (n = 6) 105,396 177,420 282,816 108 361 417 (47,634) (47,634) (95,268) (102) (342) (392) A. alba (n = 5) 444,834 522,881 967,715 93 232 277 (39,695) (39,695) (79,390) (72) (173) (202) P. sylvestris (n = 5) 212,820 93,876 306,696 152 519 628 (39,695) (39,695) (79,390) (139) (512) (610) Total 763,050 794,177 1,557,227 247 801 942 (127,024) (127,024) (254,048) (215) (512) (872) Alpha diversity RNA metabarcoding revealed approximately twice as many OTUs per sample compared to DNA metabarcoding (Fig. 1 A; χ² = 25.95, df = 1, p < 0.001). This pattern persisted for abundance-weighted alpha diversity (i.e., Inverse Simpson Index; χ² = 17.90, df = 1, p < 0.001), although the mean values of the Inverse Simpson Index were roughly ten times lower than those of OTU richness (Supplementary Fig. 2A). Overall, OTU richness per sample differed significantly between tree species (Fig. 1 B; χ² = 8.18, df = 1, p < 0.05), with the lowest richness observed in A. alba samples and the highest in P. sylvestris samples. A similar pattern was observed for the abundance-weighted diversity (χ² = 15.12, df = 2, p < 0.001) which was around ten times lower than OTU richness— F. sylvatica and P. sylvestris exhibited more than twice the diversity of A. alba (Supplementary Fig. 2B). Beta diversity When all OTUs were considered, regardless of their abundance, DNA and RNA metabarcoding revealed significantly different fungal communities associated with the same seeds (Fig. 2 A; F = 11.49, R 2 = 0.25, p < 0.001). Of the total 872 OTUs identified in the rarefied dataset, only 100 (~ 11%) were shared between the datasets, while around 13% and 75% were unique for the DNA and RNA dataset, respectively (Fig. 2 B). The proportion of OTUs shared between the two metabarcoding methods further differed across tree species: approximately 7% of OTUs were shared in P. sylvestris , compared to 13% in F. sylvatica and 21% in A. alba (Supplementary Fig. 3A–C). Comparable results were obtained using the non-rarefied dataset (Supplementary Fig. 4A-E). Significant differences in the composition of dominant fungal communities—assessed by Morisita-Horn dissimilarity, which gives greater weight to OTUs with higher read counts—were also observed between DNA and RNA datasets (Supplementary Fig. 5A; F = 1.97, R 2 = 0.05, p < 0.05). However, as indicated by R 2 values, the relative importance of the method in explaining the differences in fungal community composition decreased from 25% to 5% when an abundance-weighted beta diversity measure was used suggesting that the differences between methods were largely driven by rare OTUs. Fungal communities also differed among tree species, both when analysing the entire community (Fig. 2 C-D; F = 2.84, R 2 = 0.13, p < 0.001) and when focussing on dominant taxa (Supplementary Fig. 5B; F = 1.75, R 2 = 0.28, p < 0.001). The proportion of variation explained by tree species increased from 13% to nearly 28% when switching from incidence-based to abundance-weighted metrics, respectively, as indicated by R 2 values, suggesting that taxa with higher read counts were more host-specific than rare taxa. Taxonomic compositional differences between methods and tree species Out of a total of 1,557,227 fungal reads and 942 fungal OTUs in the non-rarefied dataset, the majority of reads and OTUs were assigned at the class, order, and family levels. At the genus level, around 72% of reads and 64% of OTUs were classified to 338 distinct genera (Supplementary Table 1). The proportion of fungal reads and OTUs assigned across taxonomic levels was broadly consistent between DNA and RNA datasets (Supplementary Table 1). The vast majority of reads and OTUs were assigned at high taxonomic levels (i.e., class, order, family) in both datasets. Approximately 74% of reads and 68% of OTUs in the DNA dataset, and 71% of reads and 65% of OTUs in the RNA dataset, were assigned to 110 and 305 unique genera, respectively. Seventy-seven genera (~ 23%) appeared in both datasets, while 33 (~ 10%) and 228 (~ 67%) were unique for DNA and RNA dataset, respectively. Approximately 85% of reads were assigned to genera in both DNA and RNA datasets of F. sylvatica , compared to 73% in A. alba and 61% in P. sylvestris (Fig. 3 ). The proportion of reads classified as “Other” (i.e., reads assigned to genera outside the top ten most abundant genera per method and tree species) was around only 3% in the A. alba DNA and RNA, and P. sylvestris DNA datasets. Higher proportion of reads classified as “Other” was found in the P. sylvestris RNA dataset (~ 14%) and, especially, in both F. sylvatica datasets (i.e., ~ 16% in the DNA and ~ 26% in the RNA dataset; Fig. 3 ). Taxonomic composition of dominant genera (i.e., the top ten most abundant genera per method and tree species) varied across tree species and sequencing methods (Fig. 3 ). The only genus consistently identified among the top ten genera across all methods and tree species was Aspergillus . Among the top ten genera identified for each method and tree species, five were shared between the DNA and RNA datasets in F. sylvatica (i.e., Alternaria , Apiognomonia , Aspergillus , Diaporthe , Gibberella ), and in A. alba (i.e., Aspergillus , Diplodia , Hormonema , Lambertella , Penicillium ), and two in P. sylvestris (i.e., Aspergillus , Hormonema ). Those shared genera accounted for a high proportion of reads in each dataset (i.e., F. sylvatica : 53% in DNA and 32% in RNA; A. alba : 60% in DNA and 55% in RNA; P. sylvestris : 44% in DNA and 40% in RNA). Indicator species analysis revealed significant, dataset-specific genus associations with either DNA or RNA datasets (Supplementary Table 2). Across tree species, the genera Mycocentrospora , Paraleptosphaeria , and Phoma were associated with the DNA datasets, whereas Fusarium , Kazachstania , and Vishniacozyma were associated with the RNA datasets, although with varying relative abundances. In F. sylvatica , four genera were associated with the DNA and 17 with the RNA dataset. Among these, several indicator genera ranked among the top ten genera per method in F. sylvatica — Paraleptosphaeria was an indicator of the DNA, whereas Vishniacozyma and Lophodermium were representative of the RNA dataset. In A. alba , nine genera were indicators of the DNA dataset and ten of the RNA dataset. Among the top ten genera per method in A. alba , Caloscypha was strongly linked to DNA, while Fusarium , Lophodermium and Neocatenulostroma emerged as indicators of the RNA dataset. For P. sylvestris , 12 genera were associated with the DNA dataset and nine with the RNA dataset. Within the top ten genera per method in P. sylvestris , Caloscypha , Fibulochlamys , Podospora , and Lambertella were DNA-associated, whereas Colletotrichum , Kazachstania , and Vishniacozyma were strongly associated with the RNA dataset. Additional genus-level associations were observed but fell outside the top ten most abundant genera per method and species (Supplementary Table 2). Lifestyle compositional differences between methods and tree species Of the total reads and OTUs in the non-rarified dataset (i.e., 1,557,227 reads and 942 OTUs), a primary lifestyle could be assigned to 1,128,634 reads (74%) representing 596 OTUs (64%). Among these, most fungal reads and OTUs belonged to saprotrophs (72% of reads, 59% of OTUs), followed by plant pathogens (25% of reads, 24% of OTUs), while mycorrhizal fungi, endophytes, and taxa categorized as “Other” each accounted for only a small fraction of the reads and OTUs (Supplementary Table 3). The proportion of fungal reads and OTUs assigned to a lifestyle was largely consistent across DNA and RNA datasets—around 74% and 71% out of 763,050 and 794,177 reads and 68% and 64% out of 247 and 801 OTUs in DNA and RNA dataset could be assigned to a lifestyle. Moreover, the proportions of reads and OTUs assigned to each lifestyle were largely similar between the metabarcoding datasets (Supplementary Table 3). Approximately 85%, 75%, and 65% of reads from F. sylvatica , A. alba , and P. sylvestris , respectively, were classified into fungal lifestyles, with varying proportions between the two datasets (Fig. 4 ). In F. sylvatica , plant pathogens accounted for 57% of DNA reads and 39% of RNA reads, while saprotrophs were represented with 29% of DNA reads and 38% of RNA reads. Saprotrophs were the most abundant guild in A. alba (60%) and P. sylvestris (52%) in both datasets. Plant pathogens constituted 15% and 12% of DNA and RNA reads in A. alba samples and 7% and 9% in P. sylvestris DNA and RNA dataset, respectively. In F. sylvatica and A. alba , mycorrhizal fungi (i.e., 740 and 392 reads) and endophytes (i.e., 217 and 233 reads) were exclusively detected in the RNA dataset. Unlike F. sylvatica and A. alba , P. sylvestris contained endophytes in both DNA and RNA datasets at similar read counts (i.e., 85 and 51, respectively). Mycorrhizal fungi in P. sylvestris were detected in both datasets but with only three reads in the DNA compared to 201 reads in the RNA dataset. Comparison of metabarcoding datasets with traditional culturing data A total of 1,553 fungal isolates were obtained across all seed samples and those were assigned to 203 morphotypes. Out of 203 representative isolates (i.e., one representative isolate per morphotype) which were selected for sequencing, 145 yielded good-quality sequences. Of those, 121 sequences were assigned to 31 distinct genera, compared to 110 and 305 genera identified through DNA and RNA metabarcoding, respectively. Of the 31 genera identified through culturing, ten (~ 32%) were found exclusively in the culturing datasets. In contrast, 17 (~ 55%) of the cultured genera were also detected in both DNA and RNA datasets and two additional cultured genera were detected in each of the metabarcoding datasets. Although the overlap between culturing and metabarcoding was broadly similar—with most cultured genera appearing in both metabarcoding datasets—the RNA dataset contained a higher number of unique genera, whereas the DNA dataset comprised of fewer unique taxa (Fig. 5 ). Most of the isolates originated from F. sylvatica seeds (i.e., 1,300 isolates representing 145 distinct morphotypes). Out of the 114 representative isolates which yielded good quality sequences 91 representative isolates could be assigned to one of 24 fungal genera (i.e., corresponding to 976 isolates and 91 morphotypes). Almost half of the cultured fungal genera (i.e., eleven) were detected by both DNA and RNA metabarcoding approaches, while an additional two (i.e., Oliveonia , Phacidium ) and three genera (i.e., Cladosporium , Chaetomium , Clonostachys ) were exclusive to the DNA or RNA datasets, respectively. Notably, eight genera were identified solely through culturing (i.e., Discosia , Epicoccum , Xylaria , Didymosphaeria , Plagiostoma , Biscogniauxia , Radulidium , Boeremia ; Table 2 ., Supplementary Fig. 6). From A. alba seeds, 231 isolates representing 44 morphotypes were obtained from seeds. Good quality sequences were obtained from 19 representative isolates, covering 97 isolates belonging to 19 morphotypes. All representative sequences could be assigned to one of the eight genera, five of which were detected in both DNA and RNA datasets (i.e., Penicillium , Trichoderma, Fusarium, Aspergillus, Talaromyces ), while three were exclusive to the culturing dataset (i.e., Sydowia , Mucor , Akanthomyces ; Table 2 ., Supplementary Fig. 6). From P. sylvestris seeds, 22 isolates representing 14 morphotypes were obtained. Good quality sequences were obtained from 12 representative isolates (morphotypes), covering 20 isolates. All but one representative isolate was assigned to one of seven genera. Three genera were shared between DNA and RNA datasets (i.e., Penicillium , Fusarium , Aspergillus ), two were shared with RNA dataset (i.e., Coniochaeta , Pseudopithomyces ) and two were detected only via culturing (i.e., Boeremia , Sydowia ; Table 2 ., Supplementary Fig. 6). Table 2 Fungal genera cultured from three tree species and their detection in DNA/RNA metabarcoding datasets Fungal genera detected in the culturing data of the three tree species (i.e., F. sylvatica , A.alba , P. sylvestris) are listed. Number of isolates belonging to each genus and number of morphotypes those isolates belong to is indicated, as well as the information about if the genera were detected in either of the two metabarcoding datasets (i.e., DNA, RNA). Genus Number of isolates Number of morphotypes Present in DNA dataset Present in RNA dataset Fagus sylvatica Alternaria 318 13 x x Apiognomonia 244 18 x x Fusarium 125 13 x x Discosia 57 5 Diaporthe 42 4 x x Epicoccum 33 5 Cladosporium 25 4 x Penicillium 24 5 x x Aureobasidium 14 3 x x Ceratobasidium 10 1 x x Oliveonia 10 1 x Clonostachys 9 4 x Trichothecium 9 1 x x Seimatosporium 8 2 x x Xylaria 8 2 Didymosphaeria 6 1 Muriphaeosphaeria 6 1 x x Chaetomium 5 1 x Biscogniauxia 5 1 Radulidium 5 1 Neosetophoma 4 2 x x Plagiostoma 3 1 Boeremia 3 1 Phacidium 3 1 x unidentified 324 54 Abies alba Penicillium 28 7 x x Trichoderma 27 2 x x Sydowia 25 2 Fusarium 7 2 x x Aspergillus 6 3 x x Mucor 2 1 Talaromyces 1 1 x x Akanthomyces 1 1 unidentified 134 25 Pinus sylvestris Penicillium 5 5 x x Fusarium 5 1 x x Boeremia 4 2 Coniochaeta 3 1 x Aspergillus 1 1 x x Sydowia 1 1 unidentified 3 3 Discussion Traditional culturing methods reveal limited fungal diversity, are resource-intensive, and have largely been replaced by the high-throughput metabarcoding approaches in studies where living cultures are not required. Metabarcoding studies are predominantly DNA-based, which does not allow differentiation between metabolically active and inactive taxa, including dormant or dead organisms. In contrast, RNA metabarcoding targets living, metabolically active organisms, potentially providing a more ecologically relevant snapshot of fungal communities compared to DNA metabarcoding. RNA metabarcoding has shown potential for studying active environmental microbiota [ 18 , 22 , 25 , 56 ]. However, the method has rarely been applied in studies of plant-associated microbiota. Given its high potential for ecological and phytosanitary monitoring, it is crucial to evaluate the RNA metabarcoding method against other commonly applied methods for studying plant-associated fungi before it can be widely implemented. In this study, we compared incidence-based and abundance-weighted seed-borne fungal diversity and community composition between DNA and RNA metabarcoding methods and tree species. Additionally, we assessed the overlap between fungal genera identified through culturing and those revealed by each metabarcoding method. Our results show that the two metabarcoding approaches capture distinct and largely non-overlapping fungal communities, with the observed differences primarily driven by rare OTUs in the RNA dataset. However, both approaches recovered similar cultured genera, which likely represent abundant and metabolically active taxa. Fungal diversity and community composition across metabarcoding and culturing datasets Contrary to our hypothesis, RNA metabarcoding revealed a higher number of fungal OTUs than DNA metabarcoding, including many unique OTUs, resulting in highly divergent community profiles. This unexpected outcome likely results from the high abundance of taxa represented by dead or dormant forms [ 57 ], as well as the large diversity of metabolically active rare or low-abundance taxa in dormant tree seeds. In complex or competitive environments—such as dormant seeds—metabolic activity often originates from low abundance taxa [ 58 ], as demonstrated by rare bacterial taxa showing high metabolic activity in atmospheric [ 56 , 59 ] and oceanic environments [ 60 ]. While DNA metabarcoding typically captures highly abundant taxa, RNA metabarcoding targets metabolically active fungi, leading to discrepancies in community profiles, especially when incidence-based measures are used. However, although abundance-weighted alpha diversity was also lower in the DNA than in the RNA dataset, the divergence between DNA- and RNA-based community profiles decreased when abundance-weighted metrics were applied. Among the most abundant genera identified per tree species and method, about half were shared between the DNA and RNA datasets in two out of three tree species. This indicates that many dominant fungi are generally metabolically active, and that discrepancies between methods largely reflect differences among rare active taxa. Methodological factors likely contributed to the observed differences in fungal communities between the DNA and RNA datasets. A first bias may have been introduced during nucleic acid extraction, as DNA and RNA were obtained from subsamples of one seed lot and using different extraction kits. Previous studies have demonstrated that the choice of extraction kit can significantly influence microbial community composition [ 61 ], a bias that could be mitigated by using co-extraction kits that allow simultaneous isolation of both DNA and RNA [ 62 ]. In addition, subsampling procedures differed between the extractions. While a single DNA extraction was performed per sample, RNA was extracted in duplicates and pooled to obtain sufficient input for reverse transcription. This variation in the amount of starting material may have led to higher alpha diversity in the RNA than DNA dataset and shifted community composition—an effect previously shown for seed mycobiota of P. sylvestris [ 63 ]. Nonetheless, our study accounted for differences in sampling effort through rarefaction and the use of abundance-weighted diversity metrics [ 64 , 65 ], helping to mitigate potential confounding effects. Furthermore, the reverse transcription step required to synthesize cDNA from RNA introduces biases such as uneven abundances of correctly transcribed molecules. The reverse transcription can also introduce artefacts due to incorrect primer binding or unpredictable reverse transcriptase behaviour [ 66 ]. These factors likely contributed to discrepancies between fungal community profiles derived from DNA and RNA. Additional bias may occur during PCR amplification of the ITS region originating from both DNA and cDNA, and library preparation, including primer-template mismatches and selective amplification of certain taxa [ 67 ]. Finally, DNA and RNA libraries were sequenced on separate runs, which may have introduced further variation, although, this effect is likely less pronounced than the other previously mentioned methodological sources of bias [ 61 ]. We found no evidence supporting our hypothesis that cultured fungi would more closely resemble the RNA-derived community than the DNA-derived one. Approximately 50% of cultured genera were detected by both metabarcoding approaches, while an additional 6% was uniquely recovered by either DNA- or RNA-based metabarcoding. These results suggest that both metabolic activity and abundance influence the success of isolation of the fungus. If culturing primarily captured metabolically active taxa, a stronger overlap with the RNA dataset—representing rare but active fungi—would have been expected. Instead, our culturing approach likely favoured abundant and competitive fungi [ 8 ]. Alternative strategies such as dilution-to-extinction culturing [ 68 ] or use of different media [ 69 ] could help recover a greater diversity of rare yet cultivable taxa and potentially increase overlap with RNA-based profiles. As expected, culturing recovered substantially lower diversity than metabarcoding, reflecting the well-known difficulty of growing many fungal taxa under standard laboratory conditions [ 8 ]. Moreover, taxa detected only by culturing may reflect stochastic recovery or primer bias between methods. As isolates are critical for species description, functional studies, and phenotypic assays, continued improvement of culture-dependent techniques remains essential. Developing more sensitive and high-throughput culturing approaches will be key to bridging molecular and culture-based views on fungal diversity. Taxonomic and functional fungal community composition across metabarcoding datasets and tree species Tree species identity was a major driver of the community composition of dominant fungi, explaining approximately 30% of the overall observed variation. The high host specificity of seed-borne fungi observed in our study aligns with findings from both woody and non-woody plants [ 70 , 71 ]. Although, several dominant fungal genera were associated with already known hosts, not all of them were previously known to occur in seeds. For example, the genus Apiognomonia was the most abundant in F. sylvatica seed, likely corresponding to A. errabunda , a common endophyte of F. sylvatica known to cause anthracnose under wet spring conditions or following insect damage [ 72 ], and was also identified in the culturing dataset. Abies alba seeds were dominated by OTUs assigned to the genus Lambertella which has been previously found in A. koreana needles [ 73 ]. Pinus sylvestris seeds were dominated by the genus Hormonema , likely corresponding to Sydowia polyspora (i.e., old name is H. dematioides ), a well-known endophyte and opportunistic pathogen in conifers and pre-emergent seed pathogen in P. ponderosa [ 6 , 74 ], which was also identified in the culturing dataset. Some fungal genera observed in seeds of specific tree species were found to be strongly associated with one of the metabarcoding datasets. For example, the genus Paraleptosphaeria was strongly associated with the DNA dataset of F. sylvatica . This genus is known to be saprobic, fungicolous, or pathogenic, and has been found in grasslands and on stems and leaves of various plants. However, Fagus has not been listed among confirmed hosts in recent studies [ 75 , 76 ]. The genus Caloscypha , likely represented by the species C. fulgens —a well-known conifer seed pathogen [ 77 ] was strongly associated with the A. alba and P. sylvestris DNA dataset. The genus Lophodermium showed strong association with the RNA data of F. sylvatica and A. alba and was represented by multiple OTUs identified to several species known as needle endophytes or pathogens of conifers [ 78 ]. Association of Lophodermium with F. sylvatica may stem from environmental or laboratory contamination. The genus Neocatenulostroma was strongly associated with the RNA dataset of A. alba and could correspond to N. abietis , species which was isolated from various substrates and was found as a saprobe or endophyte in pine needles, but also as a pathogen on a wide range of conifer hosts [ 79 ]. In P. sylvestris , Colletotrichum was strongly associated with the RNA dataset. Some Colletotrichum species are known pathogens of pines, although their occurrence in conifers is less frequent than in broadleaf hosts [ 80 ]. These results highlight the value of combining DNA and RNA-based metabarcoding to uncover both inactive and active fungal associations, some of which may have implication for tree health. Similar to findings by Mittelstrass et al. [ 3 ] most fungal reads in our study were assigned to saprotrophs and plant pathogens, which appeared in similar proportions across DNA and RNA metabarcoding datasets of all three host species. This suggests that these functional groups are not only abundant but also metabolically active in dormant tree seeds. Although dormant seeds may not seem like suitable habitats—given that saprotrophs rely on dead organic matter and plant pathogens require a susceptible host—research suggests that many potentially saprotrophic or pathogenic taxa can persist in asymptomatic plant tissues as endophytes (i.e., asymptomatic, commensal or weakly mutualistic inhabitants; [ 43 ]) and shift to saprotrophic or pathogenic lifestyle under favourable conditions [ 81 ]. Mycorrhizal fungi and endophytes were detected almost exclusively in RNA datasets, likely due to their low abundance, which may have prevented their detection by DNA metabarcoding. The sporadic occurrence across seed lots and low abundance of mycorrhizal fungi in seeds supports the notion that their seed-transmission is uncommon in trees as previously suggested [ 82 ]. However, the detection of taxa belonging to mycorrhizal fungi in seeds of various non-woody plant species [ 70 ] highlights the need to further investigate mycorrhizal transmission from seeds to seedlings. Low numbers of endophytic taxa are likely associated with their underrepresentation in databases such as the FungalTraits [ 43 ], underscoring the need for more targeted research into their diversity and function. Conclusions Our results show that DNA and RNA metabarcoding provide distinct but complementary insights into fungal communities and have valuable potential for high-throughput ecological and biosecurity monitoring. DNA metabarcoding has already been widely adopted for assessing biodiversity of different organisms in both environmental and host-associated samples. Its ability to detect living and non-living taxa provides a comprehensive perspective on total fungal community associated with the studied system. In contrast, RNA-based metabarcoding is still less commonly applied, despite its focus on metabolically active organisms which may enable real-time insights into ecosystem dynamics and functional responses to environmental change. Moreover, RNA metabarcoding might be particularly promising for biosecurity application where rapid detection and identification of living pests and pathogens is critical for decision making and timely intervention. In our study, RNA metabarcoding revealed a higher diversity, primarily of rare taxa, compared to DNA metabarcoding, making it particularly effective for detecting low-abundance pathogens, at least when dominant taxa are inactive or dormant (e.g., ensuring absence of pathogens in treated samples). However, its inability to detect metabolically inactive yet viable organisms represents a limitation as dormant pathogens may escape detection and later become active. This could be tackled with applying viability PCR (vPCR) using propidium or ethidium monoazide which enables selective suppression of DNA from membrane-compromised cells, allowing DNA-based assays to target only intact, viable fungi and thus reduce false positives from relic DNA [ 83 , 84 ]. Nevertheless, vPCR also has important limitations. Its effectiveness can vary depending on cell wall structure and environmental matrix composition, and incomplete dye penetration or binding may lead to false negatives. Moreover, its applicability to mixed environmental samples such as those that need to be ground prior to DNA/RNA extraction, remains uncertain. If technical obstacles are overcome, performing DNA metabarcoding with and without vPCR pretreatment could help differentiate living from dead organisms. Combining this with RNA-based methods can further help identify the metabolically active members of the community, leading to more accurate ecological interpretations and risk assessments. Looking ahead, advances in sequencing technologies and bioinformatics are expected to increase the sensitivity and resolution of metabarcoding methods, enabling more precise taxonomic identification and functional profiling. In the face of growing environmental pressures and accelerating species movement through trade, the combined use of DNA and RNA metabarcoding holds strong potential for ecosystem management and biosecurity, by revealing total and active communities. While culturing remains essential for studying fungal diversity, pathogenicity and functional traits, traditional approaches recover only a subset of the taxa detected by metabarcoding, typically the dominant and readily cultivable species. Future efforts should thus focus on developing more sensitive, high-throughput culturing techniques capable of capturing a broader spectrum of fungal diversity. Ultimately, integrating all three approaches—DNA- and RNA-based metabarcoding together with advanced culturing—offers the most comprehensive framework for studying fungal diversity and strengthening biosecurity, by capturing both active and inactive taxa as well as culturable and unculturable species. Declarations Ethics approval and consent to participate Not applicable. Consent for publication Not applicable. Data Availability Raw metabarcoding sequence reads have been deposited in the European Nucleotide Archive (ENA) at EMBL-EBI under accession number PRJEB96213. ITS barcode sequences generated from cultured specimens have been deposited in GenBank under accession numbers PX736467–PX736489, PX740711–PX740835, and PX745229–PX745242. All ITS sequences used in this study, together with their corresponding accession numbers, are provided in Supplementary Table 4. Competing Interests The authors declare that they have no competing interests. Funding This work was funded by the Swiss National Science Foundation SNSF (Grant: P500PB_211040) and by the Swiss Federal Office for the Environment (FOEN) (Finanzhilfevertrag betreffend wissenschaftlich-technischen Tätigkeiten im Bereich Waldschutz). René Eschen was supported by CABI. CABI is an international intergovernmental organisation, and we gratefully acknowledge the core financial support from our member countries (and lead agencies) including the UK (Foreign, Commonwealth & Development Office), China (Chinese Ministry of Agriculture and Rural Affairs), Australia (Australian Centre for International Agricultural Research), Canada (Agriculture and Agri-Food Canada), Netherlands (Directorate-General for International Cooperation) and Switzerland (Swiss Agency for Development and Cooperation). See https://www.cabi.org/about-cabi/who-we-work-with/key-donors/ for full details. Authors' contributions IF conceived and designed the study with input from RE, SP, and MC. Data collection was carried out by IF, PS, and KS. PS supervised all laboratory activities related to RNA metabarcoding. IF performed the data analysis and prepared the initial manuscript draft with input from MC. All co-authors reviewed the draft and contributed to the preparation of the final version of the manuscript Acknowledgments We thank Beatrice Tolio, Delnia Sepahvand, and Diana Marčiulynienė for their valuable assistance in the laboratory. We thank Beat Ruffner for his guidance with processing ITS sequences from cultured fungal isolates. References Hubert B, Leprince O, Buitink J. Sleeping but not defenceless: seed dormancy and protection. J Exp Bot. 2024;75:6110–24. https://doi.org/10.1093/jxb/erae213 . Nelson EB. The seed microbiome: Origins, interactions, and impacts. Plant Soil. 2018;422:7–34. https://doi.org/10.1007/s11104-017-3289-7 . Mittelstrass J, Heinzelmann R, Eschen R, Hartmann M, Kupper Q, Schneider S, et al. Metabarcoding with Illumina and Oxford Nanopore Technologies provides complementary insights into tree seed mycobiota. Environ Microbiome. 2025;20:53. https://doi.org/10.1186/s40793-025-00712-7 . Franić I, Prospero S, Hartmann M, Allan E, Auger-Rozenberg M-A, Grünwald NJ, et al. Are traded forest tree seeds a potential source of nonnative pests? Ecol Appl. 2019;29:e01971. https://doi.org/10.1002/eap.1971 . Franić I, Cleary M, Aday Kaya AG, Bragança H, Brodal G, Cech TL, et al. The Biosecurity Risks of International Forest Tree Seed Movements. Curr Forestry Rep. 2024;10:89–102. https://doi.org/10.1007/s40725-023-00211-3 . Cleary M, Oskay F, Doğmuş HT, Woodward S, Vettraino AM. Cryptic risks to forest biosecurity associated with the global movement of commercial seed. Forests. 2019;10:459. https://doi.org/10.3390/f10050459 . Rétif F, Kunz C, Calabro K, Duval C, Prado S, Bailly C, et al. Seed fungal endophytes as biostimulants and biocontrol agents to improve seed performance. Front Plant Sci. 2023;14:1260292. https://doi.org/10.3389/fpls.2023.1260292 . dos Reis JBA, Lorenzi AS, do Vale HMM. Methods used for the study of endophytic fungi: a review on methodologies and challenges, and associated tips. Arch Microbiol. 2022;204:675. https://doi.org/10.1007/s00203-022-03283-0 . Ganley RJ, Newcombe G. Fungal endophytes in seeds and needles of Pinus monticola . Mycol Res. 2006;110:318–27. https://doi.org/10.1016/j.mycres.2005.10.005 . Higgins KL, Arnold AE, Miadlikowska J, Sarvate SD, Lutzoni F. Phylogenetic relationships, host affinity, and geographic structure of boreal and arctic endophytes from three major plant lineages. Mol Phylogenet Evol. 2007;42:543–55. https://doi.org/10.1016/j.ympev.2006.07.012 . Zimmerman NB, Vitousek PM. Fungal endophyte communities reflect environmental structuring across a Hawaiian landscape. Proc Natl Acad Sci U S A. 2012;109:13022–7. https://doi.org/10.1073/pnas.1209872109 . Mapuranga J, Zhang N, Zhang L, Chang J, Yang W. Infection Strategies and Pathogenicity of Biotrophic Plant Fungal Pathogens. Front Microbiol. 2022;13:799396. https://doi.org/10.3389/fmicb.2022.799396 . Eschen R, De Groot M, Glavendekić M, Lacković N, Matošević D, Morales-Rodriguez C, et al. Spotting the pests of tomorrow—Sampling designs for detection of species associations with woody plants. J Biogeogr. 2019;46:2159–73. https://doi.org/10.1111/jbi.13670 . Nilsson RH, Anslan S, Bahram M, Wurzbacher C, Baldrian P, Tedersoo L. Mycobiome diversity: high-throughput sequencing and identification of fungi. Nat Rev Microbiol. 2019;17:95–109. https://doi.org/10.1038/s41579-018-0116-y . Tedersoo L, Drenkhan R, Anslan S, Morales-Rodriguez C, Cleary M. High-throughput identification and diagnostics of pathogens and pests: overview and practical recommendations. Mol Ecol Resour. 2018;19:47–76. https://doi.org/10.1111/1755-0998.12959 . Schoch CL, Seifert KA, Huhndorf S, Robert V, Spouge JL, Levesque CA, et al. Nuclear ribosomal internal transcribed spacer (ITS) region as a universal DNA barcode marker for Fungi. Proc Natl Acad Sci U S A. 2012;109:6241–6. https://doi.org/10.1073/pnas.1117018109 . Carini P, Marsden PJ, Leff JW, Morgan EE, Strickland MS, Fierer N. Relic DNA is abundant in soil and obscures estimates of soil microbial diversity. Nat Microbiol. 2017;2:16242. https://doi.org/10.1038/nmicrobiol.2016.242 . Pochon X, Zaiko A, Fletcher LM, Laroche O, Wood SA. Wanted dead or alive? Using metabarcoding of environmental DNA and RNA to distinguish living assemblages for biosecurity applications. PLoS ONE. 2017;12:e0187636. https://doi.org/10.1371/journal.pone.0187636 . Wong B, Leal I, Feau N, Dale A, Uzunovic A, Hamelin RC. Molecular assays to detect the presence and viability of Phytophthora ramorum and Grosmannia clavigera . PLoS ONE. 2020;15:e0221742. https://doi.org/10.1371/journal.pone.0221742 . Leal I, Feau N, Uzunovic A, Foord B, Hamelin RC. A Molecular Method to Assess Viability of Phytophthora in Infected Wood Following Phytosanitary Heat Treatment. PhytoFrontiersTM. 2024;4:690–9. https://doi.org/10.1094/phytofr-05-24-0056-r . Anderson IC, Parkin PI. Detection of active soil fungi by RT-PCR amplification of precursor rRNA molecules. J Microbiol Methods. 2007;68:248–53. https://doi.org/10.1016/j.mimet.2006.08.005 . Yates MC, Derry AM, Cristescu ME, Environmental RNA. A Revolution in Ecological Resolution? Trends in Ecology and Evolution. 2021;36:601–9. https://doi.org/10.1016/j.tree.2021.03.001 Baldrian P, Kolaiřík M, Štursová M, Kopecký J, Valášková V, Větrovský T, et al. Active and total microbial communities in forest soil are largely different and highly stratified during decomposition. ISME J. 2012;6:248–58. https://doi.org/10.1038/ismej.2011.95 . Rajala T, Peltoniemi M, Hantula J, Mäkipää R, Pennanen T. RNA reveals a succession of active fungi during the decay of Norway spruce logs. Fungal Ecol. 2011;4:437–48. https://doi.org/10.1016/j.funeco.2011.05.005 . Xue Z, Tian W, Han Y, Li S, Guo J, He H, et al. Environmental RNA metabarcoding for ballast water microbial diversity: Minimizing false positives. Sci Total Environ. 2024;955:176902. https://doi.org/10.1016/j.scitotenv.2024.176902 . Kong H, Yang EJ, Jiao N, Lee Y, Jung J, Cho KH, et al. RNA outperforms DNA-based metabarcoding in assessing the diversity and response of microeukaryotes to environmental variables in the Arctic Ocean. Sci Total Environ. 2023;876:162608. https://doi.org/10.1016/j.scitotenv.2023.162608 . Nguyen D, Boberg J, Ihrmark K, Stenström E, Stenlid J. Do foliar fungal communities of Norway spruce shift along a tree species diversity gradient in mature European forests? Fungal Ecol. 2016;23:97–108. https://doi.org/10.1016/j.funeco.2016.07.003 . Wutkowska M, Vader A, Mundra S, Cooper EJ, Eidesen PB. Dead or alive; Or does it really matter? Level of congruency between trophic modes in total and active fungal communities in high arctic soil. Front Microbiol. 2019;9:3243. https://doi.org/10.3389/fmicb.2018.03243 . Camenzind T, Weimershaus P, Lehmann A, Aguilar-Trigueros C, Rillig MC. Soil fungi invest into asexual sporulation under resource scarcity, but trait spaces of individual isolates are unique. Environ Microbiol. 2022;24:2962–78. https://doi.org/10.1111/1462-2920.16012 . Chen DV, Slowinski SP, Kido AK, Bruns EL. High temperatures reduce growth, infection, and transmission of a naturally occurring fungal plant pathogen. Ecology. 2024;105:e4373. https://doi.org/10.1002/ecy.4373 . Jo T, Tsuri K, Hirohara T, Yamanaka H. Warm temperature and alkaline conditions accelerate environmental RNA degradation. Environ DNA. 2023;5:836–48. https://doi.org/10.1002/edn3.334 . Gamboa MA, Laureano S, Bayman P. Measuring diversity of endophytic fungi in leaf fragments: Does size matter? Mycopathologia. 2003;156:41–5. https://doi.org/10.1023/A:1021362217723 White T, Bruns T, Lee S, Taylor J. Amplification and direct sequencing of fungal ribosomal RNA genes for phylogenetics. In: Innis MA, Gelfand DH, Sninsky JJ WT, editors. PCR Protocols: A Guide to Methods and Applications. New York: Academic Press Inc; 1990. pp. 315–22. Altschup SF, Gish W, Miller W, Myers EW, Lipman DJ. Basic Local Alignment Search Tool. J Mol Biol. 1990;215:403–10. Zhang Z, Schwartz S, Wagner L, Miller W. A Greedy Algorithm for Aligning DNA Sequences. J Comput Biol. 2000;7:203–14. Bolyen E, Rideout JR, Dillon MR, Bokulich NA, Abnet CC, Al-Ghalith GA, et al. Reproducible, interactive, scalable and extensible microbiome data science using QIIME 2. Nat Biotechnol. 2019;37:852–7. https://doi.org/10.1038/s41587-019-0209-9 . Martin M. Cutadapt removes adapter sequences from high-throughput sequencing reads. EMBnet J. 2011;17:10–2. https://doi.org/https://doi.org/10.14806/ej.17.1.200 . Callahan BJ, McMurdie PJ, Rosen MJ, Han AW, Johnson AJA, Holmes SP. DADA2: High-resolution sample inference from Illumina amplicon data. Nat Methods. 2016;13:581–3. https://doi.org/10.1038/nmeth.3869 . Rognes T, Flouri T, Nichols B, Quince C, Mahé F. VSEARCH: A versatile open source tool for metagenomics. PeerJ. 2016;10:e2584. https://doi.org/10.7717/peerj.2584 . Pedregosa F, Varoquaux G, Gramfort A, Michel V, Thirion B, Grisel O, et al. Scikit-learn: Machine Learning in Python. J Mach Learn Res. 2011;12:2825–30. Kõljalg U, Nilsson RH, Abarenkov K, Tedersoo L, Taylor AFS, Bahram M, et al. Towards a unified paradigm for sequence-based identification of fungi. Mol Ecol. 2013;22:5271–7. https://doi.org/10.1111/mec.12481 . Monard C, Gantner S, Stenlid J. Utilizing ITS1 and ITS2 to study environmental fungal diversity using pyrosequencing. FEMS Microbiol Ecol. 2013;84:165–75. https://doi.org/10.1111/1574-6941.12046 . Põlme S, Abarenkov K, Henrik Nilsson R, Lindahl BD, Clemmensen KE, Kauserud H, et al. FungalTraits: a user-friendly traits database of fungi and fungus-like stramenopiles. Fungal Divers. 2020;105:1–16. https://doi.org/10.1007/s13225-020-00466-2 . R Core Team. R: A Language and Environment for Statistical Computing. 2024. Chen J, Zhang X, Yang L, Zhang L, _GUniFrac. Generalized UniFrac Distances, Distance-Based Multivariate Methods and Feature-Based Univariate Methods for Microbiome Data Analysis_. 2025. Schloss PD. Rarefaction is currently the best approach to control for uneven sequencing effort in amplicon sequence analyses. mSphere. 2024;9:e00354–23. https://doi.org/10.1128/msphere.00354-23 . Cassol I, Ibañez M, Bustamante JP. Key features and guidelines for the application of microbial alpha diversity metrics. Sci Rep. 2025;15:622. https://doi.org/10.1038/s41598-024-77864-y . Jost L. Entropy and diversity. Oikos. 2006;113:363–75. https://doi.org/10.1111/j.2006.0030-1299.14714.x . Chao A, Chiu C-H, Jost L. Unifying Species Diversity, Phylogenetic Diversity, Functional Diversity, and Related Similarity and Differentiation Measures Through Hill Numbers. Annu Rev Ecol Evol Syst. 2014;45:297–324. https://doi.org/10.1146/annurev-ecolsys-120213-091540 . Venables WN, Ripley BD. Modern Applied Statistics with S. Fourth edition. New York: Springer; 2002. https://doi.org/10.1007/978-0-387-21706-2 . https://doi.org/. Fox J, Weisberg S, An R. Companion to Applied Regression. Third. Thousand Oaks: Sage; 2019. Lenth RV, Least-Squares, Means. The R Package lsmeans. J Stat Softw. 2016;69:1–33. https://doi.org/10.18637/jss.v069.i01 . Oksanen J, Blanchet FG, Friendly M, Kindt R, Legendre P, Dan, McGlinn et al. _vegan: Community Ecology Package_. R package version 2.5-3. 2025. McMurdie PJ, Holmes S. Waste Not, Want Not: Why Rarefying Microbiome Data Is Inadmissible. PLoS Comput Biol. 2014;10:e1003531. https://doi.org/10.1371/journal.pcbi.1003531 . Cáceres M, De, Legendre P. Associations between species and groups of sites: indices and statistical inference. Ecology. 2009;90:3566–74. https://doi.org/10.1890/08-1823.1 . Chen Y, Zhu X, Hou Z, Wang Y, Zhou Y, Wang L, et al. RNA-Based Analysis Reveals High Diversity of Plant-Associated Active Fungi in the Atmosphere. Front Microbiol. 2021;12:683266. https://doi.org/10.3389/fmicb.2021.683266 . Jones SE, Lennon JT. Dormancy contributes to the maintenance of microbial diversity. Proc Natl Acad Sci U S A. 2010;107:5881–6. https://doi.org/10.1073/pnas.0912765107 . Camenzind T, Lehmann A, Ahland J, Rumpel S, Rillig MC. Trait-based approaches reveal fungal adaptations to nutrient-limiting conditions. Environ Microbiol. 2020;22:3548–60. https://doi.org/10.1111/1462-2920.15132 . Klein AM, Bohannan BJM, Jaffe DA, Levin DA, Green JL. Molecular evidence for metabolically active bacteria in the atmosphere. Front Microbiol. 2016;7:772. https://doi.org/10.3389/fmicb.2016.00772 . Campbell BJ, Yu L, Heidelberg JF, Kirchman DL. Activity of abundant and rare bacteria in a coastal ocean. Proc Natl Acad Sci U S A. 2011;108:12776–81. https://doi.org/10.1073/pnas.1101405108 . Pu Y, Zhou X, Cai H, Lou T, Liu C, Kong M, et al. Impact of DNA Extraction Methods on Gut Microbiome Profiles: A Comparative Metagenomic Study. Phenomics. 2025;5:76–90. https://doi.org/10.1007/s43657-025-00232-x . 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. https://doi.org/10.1371/journal.pone.0121659 . Oskay F, Vettraino AM, Doğmuş HT, Lehtijärvi A, Woodward S, Cleary M. Seed quantity affects the fungal community composition detected using metabarcoding. Sci Rep. 2022;12:3060. https://doi.org/10.1038/s41598-022-06997-9 . Walker SC, Poos MS, Jackson DA. Functional rarefaction: Estimating functional diversity from field data. Oikos. 2008;117:286–96. https://doi.org/10.1111/j.2007.0030-1299.16171.x . McCoy CO, Matsen FA. Abundance-weighted phylogenetic diversity measures distinguish microbial community states and are robust to sampling depth. PeerJ. 2013;1:e157. https://doi.org/10.7717/peerj.157 . Verwilt J, Mestdagh P, Vandesompele JO. Artifacts and biases of the reverse transcription reaction in RNA sequencing. RNA. 2023;29:889–97. https://doi.org/10.1261/rna.079623.123 . Aird D, Ross MG, Chen WS, Danielsson M, Fennell T, Russ C, et al. Analyzing and minimizing PCR amplification bias in Illumina sequencing libraries. Genome Biol. 2011;12:R18. https://doi.org/10.1186/gb-2011-12-2-r18 . Unterseher M, Schnittler M. Dilution-to-extinction cultivation of leaf-inhabiting endophytic fungi in beech ( Fagus sylvatica L.) - Different cultivation techniques influence fungal biodiversity assessment. Mycol Res. 2009;113:645–54. https://doi.org/10.1016/j.mycres.2009.02.002 . Ordonez A, Hussain U, Cambon MC, Golyshin PN, Downie J, McDonald JE. Evaluating agar-plating and dilution-to-extinction isolation methods for generating oak-associated microbial culture collections. ISME Commun. 2025;5:ycaf019. https://doi.org/10.1093/ismeco/ycaf019 . Simonin M, Briand M, Chesneau G, Rochefort A, Marais C, Sarniguet A, et al. Seed microbiota revealed by a large-scale meta-analysis including 50 plant species. New Phytol. 2022;234:1448–63. https://doi.org/10.1111/nph.18037 . Franić I, Eschen R, Allan E, Hartmann M, Schneider S, Prospero S. Drivers of richness and community composition of fungal endophytes of tree seeds. FEMS Microbiol Ecol. 2020;96:fiaa166. https://doi.org/10.1093/femsec/fiaa166 . Boroń P, Grad B, Nawrot-Chorabik K, Kowalski T. The genetic relationship of Apiognomonia errabunda and related species. Mycologia. 2019;111:541–50. https://doi.org/10.2307/26756638 . Choi JW, Park E, Eo JK. An Unrecorded Genus Lambertella Höhn. (Rutstroemiaceae) and Its Unrecorded Species in Korea. Korean J Mycol. 2021;49:127–31. https://doi.org/10.4489/KJM.20210013 . Ridout M, Newcombe G. Sydowia polyspora is both a foliar endophyte and a preemergent seed pathogen in Pinus ponderosa . Plant Dis. 2018;102:640–4. https://doi.org/10.1094/PDIS-07-17-1074-RE . Gao Y, de Farias ARG, Jiang HB, Karunarathna SC, Xu JC, Tibpromma S, et al. Morphological and Phylogenetic Characterisations Reveal Four New Species in Leptosphaeriaceae (Pleosporales, Dothideomycetes). J Fungi. 2023;9. https://doi.org/10.3390/jof9060612 . Ariyawansa HA, Phukhamsakda C, Thambugala KM, Bulgakov TS, Wanasinghe DN, Perera RH, et al. Revision and phylogeny of Leptosphaeriaceae. Fungal Divers. 2015;74:19–51. https://doi.org/10.1007/s13225-015-0349-2 . Schröder T, Kehr R, Hüttermann A. First report of the seed-pathogen Geniculodendron pyriforme , the imperfect state of the ascomycete Caloscypha fulgens , on imported conifer seeds in Germany. Pathol. 2002;32:225–30. https://doi.org/10.1046/j.1439-0329.2002.00288.x . Salas-Lizana R, Oono R. Double-digest RADseq loci using standard Illumina indexes improve deep and shallow phylogenetic resolution of Lophodermium , a widespread fungal endophyte of pine needles. Ecol Evol. 2018;8:6638–51. https://doi.org/10.1002/ece3.4147 . Phukhamsakda C, Nilsson RH, Bhunjun CS, de Farias ARG, Sun YR, Wijesinghe SN, et al. The numbers of fungi: contributions from traditional taxonomic studies and challenges of metabarcoding. Fungal Divers. 2022;114:327–86. https://doi.org/10.1007/s13225-022-00502-3 . Talhinhas P, Baroncelli R. Colletotrichum species and complexes: geographic distribution, host range and conservation status. Fungal Divers. 2021;110:109–98. https://doi.org/10.1007/s13225-021-00491-9 . Liao C, Doilom M, Jeewon R, Hyde KD, Manawasinghe IS, Chethana KWT, et al. Challenges and update on fungal endophytes: classification, definition, diversity, ecology, evolution and functions. Fungal Divers. 2025;131:301–67. https://doi.org/10.1007/s13225-025-00550-5 . Cram MM, Kasten Dumroese R. Mycorrhizae in Forest Tree Nurseries. In: Cram MM, Frank MS, Mallams KM, editors. Forest Nursery Pests. Washington, DC, USA: USDA, Forest Service; 2012. pp. 20–3. Nocker A, Sossa KE, Camper AK. Molecular monitoring of disinfection efficacy using propidium monoazide in combination with quantitative PCR. J Microbiol Methods. 2007;70:252–60. https://doi.org/10.1016/j.mimet.2007.04.014 . Nocker A, Camper AK. Selective removal of DNA from dead cells of mixed bacterial communities by use of ethidium monoazide. Appl Environ Microbiol. 2006;72:1997–2004. https://doi.org/10.1128/AEM.72.3.1997-2004.2006 . Additional Declarations No competing interests reported. Supplementary Files Franicetal.2025SupplementaryMaterial.docx Cite Share Download PDF Status: Under Review Version 1 posted Editorial decision: Revision requested 27 Mar, 2026 Reviews received at journal 25 Mar, 2026 Reviews received at journal 17 Mar, 2026 Reviews received at journal 14 Mar, 2026 Reviews received at journal 14 Mar, 2026 Reviewers agreed at journal 16 Feb, 2026 Reviewers agreed at journal 12 Feb, 2026 Reviewers agreed at journal 12 Feb, 2026 Reviewers agreed at journal 11 Feb, 2026 Reviewers invited by journal 11 Feb, 2026 Editor assigned by journal 26 Jan, 2026 Submission checks completed at journal 23 Jan, 2026 First submitted to journal 19 Jan, 2026 You are reading this latest preprint version Research Square lets you share your work early, gain feedback from the community, and start making changes to your manuscript prior to peer review in a journal. As a division of Research Square Company, we’re committed to making research communication faster, fairer, and more useful. We do this by developing innovative software and high quality services for the global research community. Our growing team is made up of researchers and industry professionals working together to solve the most critical problems facing scientific publishing. Also discoverable on Platform About Our Team In Review Editorial Policies Advisory Board Help Center Resources Author Services Accessibility API Access RSS feed Manage Cookie Preferences © Research Square 2026 | ISSN 2693-5015 (online) Privacy Policy Terms of Service Do Not Sell My Personal Information {\"props\":{\"pageProps\":{\"initialData\":{\"identity\":\"rs-8095745\",\"acceptedTermsAndConditions\":true,\"allowDirectSubmit\":false,\"archivedVersions\":[],\"articleType\":\"Research Article\",\"associatedPublications\":[],\"authors\":[{\"id\":591443430,\"identity\":\"4bda9fa3-6af6-4794-9170-aba5e86a8cab\",\"order_by\":0,\"name\":\"Iva Franić\",\"email\":\"data:image/png;base64,iVBORw0KGgoAAAANSUhEUgAAAZAAAAAyAQMAAABI0h/eAAAABlBMVEX///8AAABVwtN+AAAACXBIWXMAAA7EAAAOxAGVKw4bAAAAyElEQVRIiWNgGAWjYBACNhBRAWYwH2BgKGBgMCBKyxkwgy0BrJ6gFgaYFgYGHgPitPBJJD9gOLjHLo+P/cw3yR8GDPLmBB0mkWbAcOBZcjEbT+42aaBFhjsbCGnhOWDA/OEAc2KbBO+220BXJRgcIKjl+AeGAwfqgVp4nt38QZQW9h6gww4cBmlhu8FDpJaCAwcOHE9s40kz/81jIGG4gZAW+Wb2jQ8OHKhOnN9++LHhjwobeYK2gACyGgki1I+CUTAKRsEoIAgAaL88KmR6W2wAAAAASUVORK5CYII=\",\"orcid\":\"\",\"institution\":\"Swiss Federal Institute for Forest, Snow and Landscape Research WSL\",\"correspondingAuthor\":true,\"prefix\":\"\",\"firstName\":\"Iva\",\"middleName\":\"\",\"lastName\":\"Franić\",\"suffix\":\"\"},{\"id\":591443431,\"identity\":\"69910f62-1e65-497a-9960-cc030ccfcb82\",\"order_by\":1,\"name\":\"Patrick Sherwood\",\"email\":\"\",\"orcid\":\"\",\"institution\":\"Swedish University of Agricultural Sciences, Southern Swedish Forest Research Centre\",\"correspondingAuthor\":false,\"prefix\":\"\",\"firstName\":\"Patrick\",\"middleName\":\"\",\"lastName\":\"Sherwood\",\"suffix\":\"\"},{\"id\":591443432,\"identity\":\"783f1a14-b042-461b-86b7-7f6e753eaff3\",\"order_by\":2,\"name\":\"Kinga Stolarek\",\"email\":\"\",\"orcid\":\"\",\"institution\":\"Swedish University of Agricultural Sciences, Southern Swedish Forest Research Centre\",\"correspondingAuthor\":false,\"prefix\":\"\",\"firstName\":\"Kinga\",\"middleName\":\"\",\"lastName\":\"Stolarek\",\"suffix\":\"\"},{\"id\":591443433,\"identity\":\"b8fec7f5-fd42-4950-995f-e5ca22d826d2\",\"order_by\":3,\"name\":\"René Eschen\",\"email\":\"\",\"orcid\":\"\",\"institution\":\"Federal Office for Agriculture\",\"correspondingAuthor\":false,\"prefix\":\"\",\"firstName\":\"René\",\"middleName\":\"\",\"lastName\":\"Eschen\",\"suffix\":\"\"},{\"id\":591443434,\"identity\":\"2eb125e4-ba51-49dd-a4ae-0604d39439de\",\"order_by\":4,\"name\":\"Jana Orbach\",\"email\":\"\",\"orcid\":\"\",\"institution\":\"Swiss Federal Institute for Forest, Snow and Landscape Research WSL\",\"correspondingAuthor\":false,\"prefix\":\"\",\"firstName\":\"Jana\",\"middleName\":\"\",\"lastName\":\"Orbach\",\"suffix\":\"\"},{\"id\":591443435,\"identity\":\"8069409a-da3d-4544-b796-2b19fdb6aa9d\",\"order_by\":5,\"name\":\"Simone Prospero\",\"email\":\"\",\"orcid\":\"\",\"institution\":\"Swiss Federal Institute for Forest, Snow and Landscape Research WSL\",\"correspondingAuthor\":false,\"prefix\":\"\",\"firstName\":\"Simone\",\"middleName\":\"\",\"lastName\":\"Prospero\",\"suffix\":\"\"},{\"id\":591443436,\"identity\":\"15ff801d-21f6-42f2-97ad-ba4c2581b390\",\"order_by\":6,\"name\":\"Michelle Cleary\",\"email\":\"\",\"orcid\":\"\",\"institution\":\"Swedish University of Agricultural Sciences, Southern Swedish Forest Research Centre\",\"correspondingAuthor\":false,\"prefix\":\"\",\"firstName\":\"Michelle\",\"middleName\":\"\",\"lastName\":\"Cleary\",\"suffix\":\"\"}],\"badges\":[],\"createdAt\":\"2025-11-12 11:23:53\",\"currentVersionCode\":1,\"declarations\":\"\",\"doi\":\"10.21203/rs.3.rs-8095745/v1\",\"doiUrl\":\"https://doi.org/10.21203/rs.3.rs-8095745/v1\",\"draftVersion\":[],\"editorialEvents\":[],\"editorialNote\":\"\",\"failedWorkflow\":false,\"files\":[{\"id\":102830952,\"identity\":\"6ef18288-e734-4fa0-94a2-e391167c5ff9\",\"added_by\":\"auto\",\"created_at\":\"2026-02-17 09:56:45\",\"extension\":\"png\",\"order_by\":1,\"title\":\"Figure 1\",\"display\":\"\",\"copyAsset\":false,\"role\":\"figure\",\"size\":340666,\"visible\":true,\"origin\":\"\",\"legend\":\"\\u003cp\\u003e\\u003cstrong\\u003eOTU richness per seed lot across metabarcoding methods and tree species.\\u003c/strong\\u003e OTU richness (i.e., number of observed OTUs) per seed lot comparing (A) different methods and (B) tree species is shown. Bars represent estimated means ± standard errors (SE). The model was run on the rarefied dataset. Different letters above the bars indicate statistically significant differences (p \\u0026lt; 0.05).\\u003c/p\\u003e\",\"description\":\"\",\"filename\":\"image1.png\",\"url\":\"https://assets-eu.researchsquare.com/files/rs-8095745/v1/7ef5f1be1ab40c5ae87bda14.png\"},{\"id\":102830932,\"identity\":\"a10aede6-cab7-4140-a533-67b9ebd6f989\",\"added_by\":\"auto\",\"created_at\":\"2026-02-17 09:56:30\",\"extension\":\"png\",\"order_by\":2,\"title\":\"Figure 2\",\"display\":\"\",\"copyAsset\":false,\"role\":\"figure\",\"size\":546820,\"visible\":true,\"origin\":\"\",\"legend\":\"\\u003cp\\u003e\\u003cstrong\\u003eDifferences in fungal community composition across metabarcoding methods and tree species.\\u003c/strong\\u003e\\u003cem\\u003e \\u003c/em\\u003eDifferences in fungal community composition between samples assessed with DNA and RNA metabarcoding (A, B) and samples belonging to different tree species (C, D) are shown. Non-metric multidimensional scaling ordination plots (A, C) are based on Sørensen distances and represent incidence-based community patterns. Stress value is 0.12. Points represent individual samples, coloured by method (A) or tree species (C). Ellipses represent 95% confidence intervals. Venn diagrams (B, D) are based on incidence data and illustrate the overlap in fungal operational taxonomic units between the two methods (B) or across three tree species (D), highlighting both shared and unique OTUs. All visualisations are based on the rarefied dataset.\\u003c/p\\u003e\",\"description\":\"\",\"filename\":\"image2.png\",\"url\":\"https://assets-eu.researchsquare.com/files/rs-8095745/v1/14bbd8d7c21759965a8195e1.png\"},{\"id\":102830956,\"identity\":\"eeee4843-4467-48b3-bcb6-689fd2999113\",\"added_by\":\"auto\",\"created_at\":\"2026-02-17 09:56:47\",\"extension\":\"png\",\"order_by\":3,\"title\":\"Figure 3\",\"display\":\"\",\"copyAsset\":false,\"role\":\"figure\",\"size\":810470,\"visible\":true,\"origin\":\"\",\"legend\":\"\\u003cp\\u003e\\u003cstrong\\u003eDominant fungal genera in seeds of three tree species from DNA and RNA metabarcoding.\\u003c/strong\\u003e Taxonomic profiles of dominant fungal communities in seeds of \\u003cem\\u003eFagus sylvatica\\u003c/em\\u003e, \\u003cem\\u003eAbies alba\\u003c/em\\u003e, and \\u003cem\\u003ePinus sylvestris\\u003c/em\\u003e, based on DNA and RNA metabarcoding are shown. Bars show the relative abundance of the ten most abundant fungal genera per tree species and metabarcoding method. Reads assigned to genera outside the top ten are grouped under “Other,” while those not classified at genus level are labelled “Unassigned.” Calculations are based on non-rarefied dataset.\\u003c/p\\u003e\",\"description\":\"\",\"filename\":\"image3.png\",\"url\":\"https://assets-eu.researchsquare.com/files/rs-8095745/v1/3c7e85fb4d4e03832fd49623.png\"},{\"id\":102830954,\"identity\":\"c11fa057-3fe8-4b33-9ebc-31cacecdfdb4\",\"added_by\":\"auto\",\"created_at\":\"2026-02-17 09:56:45\",\"extension\":\"png\",\"order_by\":4,\"title\":\"Figure 4\",\"display\":\"\",\"copyAsset\":false,\"role\":\"figure\",\"size\":815685,\"visible\":true,\"origin\":\"\",\"legend\":\"\\u003cp\\u003e\\u003cstrong\\u003eFungal lifestyle composition in seeds of three tree species from DNA and RNA metabarcoding\\u003c/strong\\u003e\\u003cem\\u003e\\u003cstrong\\u003e.\\u003c/strong\\u003e\\u003c/em\\u003e\\u003cem\\u003e \\u003c/em\\u003eLifestyle profiles of fungal communities in seeds of \\u003cem\\u003eFagus sylvatica\\u003c/em\\u003e, \\u003cem\\u003eAbies alba\\u003c/em\\u003e, and \\u003cem\\u003ePinus sylvestris\\u003c/em\\u003e, based on DNA and RNA metabarcoding are shown. Fungal lifestyles were assigned using FungalTraits database. Bars show the relative abundance of different fungal lifestyles per tree species and metabarcoding method. Percentages are not shown for lifestyles representing less than 1% of reads per method and tree species. Calculations are based on non-rarefied dataset.\\u003c/p\\u003e\",\"description\":\"\",\"filename\":\"image4.png\",\"url\":\"https://assets-eu.researchsquare.com/files/rs-8095745/v1/c27435acfb1de52d3a5c4b0e.png\"},{\"id\":102830944,\"identity\":\"7e8379af-e3f8-46b1-b8b8-58219e44fe08\",\"added_by\":\"auto\",\"created_at\":\"2026-02-17 09:56:34\",\"extension\":\"png\",\"order_by\":5,\"title\":\"Figure 5\",\"display\":\"\",\"copyAsset\":false,\"role\":\"figure\",\"size\":598193,\"visible\":true,\"origin\":\"\",\"legend\":\"\\u003cp\\u003e\\u003cstrong\\u003eOverlap of fungal genera detected by culturing and DNA/RNA metabarcoding methods\\u003c/strong\\u003e\\u003cem\\u003e\\u003cstrong\\u003e.\\u003c/strong\\u003e\\u003c/em\\u003e\\u003cem\\u003e \\u003c/em\\u003eVenn diagram illustrates the overlap in fungal genera detected by culturing and two metabarcoding methods (i.e., DNA and RNA), highlighting shared and unique OTUs. The visualization is based on the non-rarefied dataset.\\u003c/p\\u003e\",\"description\":\"\",\"filename\":\"image5.png\",\"url\":\"https://assets-eu.researchsquare.com/files/rs-8095745/v1/7cfbccb3004d5ddb5496a2ee.png\"},{\"id\":102830970,\"identity\":\"d2f214c7-65d8-44aa-9b8a-4db8069186e5\",\"added_by\":\"auto\",\"created_at\":\"2026-02-17 09:57:01\",\"extension\":\"pdf\",\"order_by\":0,\"title\":\"\",\"display\":\"\",\"copyAsset\":false,\"role\":\"manuscript-pdf\",\"size\":3950239,\"visible\":true,\"origin\":\"\",\"legend\":\"\",\"description\":\"\",\"filename\":\"manuscript.pdf\",\"url\":\"https://assets-eu.researchsquare.com/files/rs-8095745/v1/12631af5-932a-45fb-8382-58ebd391a7a7.pdf\"},{\"id\":102830962,\"identity\":\"bcb3420b-f3ac-44d0-a83f-5e94002d9149\",\"added_by\":\"auto\",\"created_at\":\"2026-02-17 09:56:53\",\"extension\":\"docx\",\"order_by\":0,\"title\":\"\",\"display\":\"\",\"copyAsset\":false,\"role\":\"supplement\",\"size\":1093851,\"visible\":true,\"origin\":\"\",\"legend\":\"\",\"description\":\"\",\"filename\":\"Franicetal.2025SupplementaryMaterial.docx\",\"url\":\"https://assets-eu.researchsquare.com/files/rs-8095745/v1/ddebaad6d7bde9b81fe4e330.docx\"}],\"financialInterests\":\"No competing interests reported.\",\"formattedTitle\":\"DNA and RNA metabarcoding reveal distinct seed-borne mycobiota\",\"fulltext\":[{\"header\":\"Background\",\"content\":\"\\u003cp\\u003eSeeds are essential for plant reproduction and dispersal and act as selective microbial habitats. Packed with oils, starches, and sugars, their resources are protected by seed coats and chemical defences [\\u003cspan citationid=\\\"CR1\\\" class=\\\"CitationRef\\\"\\u003e1\\u003c/span\\u003e], allowing access only to specialized microorganisms [\\u003cspan citationid=\\\"CR2\\\" class=\\\"CitationRef\\\"\\u003e2\\u003c/span\\u003e]. Molecular studies show that tree seeds host diverse fungi [\\u003cspan citationid=\\\"CR3\\\" class=\\\"CitationRef\\\"\\u003e3\\u003c/span\\u003e], including plant pathogens [\\u003cspan additionalcitationids=\\\"CR5\\\" citationid=\\\"CR4\\\" class=\\\"CitationRef\\\"\\u003e4\\u003c/span\\u003e\\u0026ndash;\\u003cspan citationid=\\\"CR6\\\" class=\\\"CitationRef\\\"\\u003e6\\u003c/span\\u003e], and fungi that protect plants against biotic and abiotic stresses [\\u003cspan citationid=\\\"CR7\\\" class=\\\"CitationRef\\\"\\u003e7\\u003c/span\\u003e]. However, few of these fungi may be metabolically active, while the remainder are dead or dormant. Differentiating these fractions is key to understanding the roles of seed-associated fungi in their hosts and their impact on tree health.\\u003c/p\\u003e \\u003cp\\u003ePlant-associated fungi have been traditionally studied by growing them from plant tissues on nutrient media and identifying the isolates based on morphology and/or with molecular techniques [\\u003cspan citationid=\\\"CR8\\\" class=\\\"CitationRef\\\"\\u003e8\\u003c/span\\u003e]. This approach has been widely used to explore fungal biodiversity and ecology [\\u003cspan additionalcitationids=\\\"CR10\\\" citationid=\\\"CR9\\\" class=\\\"CitationRef\\\"\\u003e9\\u003c/span\\u003e\\u0026ndash;\\u003cspan citationid=\\\"CR11\\\" class=\\\"CitationRef\\\"\\u003e11\\u003c/span\\u003e] and remains essential for obtaining living cultures, describing new species and conducting phenotyping assays (e.g., host range, virulence, optimal growth temperatures). For studying new host-pathogen associations, living cultures are also necessary to fulfil Koch\\u0026rsquo;s postulates. However, culturing captures only a subset of fungal diversity\\u0026mdash;mainly abundant fungi\\u0026mdash;excluding taxa that are difficult or impossible to isolate [\\u003cspan citationid=\\\"CR3\\\" class=\\\"CitationRef\\\"\\u003e3\\u003c/span\\u003e], but which may affect tree health [\\u003cspan citationid=\\\"CR12\\\" class=\\\"CitationRef\\\"\\u003e12\\u003c/span\\u003e]. Moreover, culturing is time-consuming and costly [\\u003cspan citationid=\\\"CR13\\\" class=\\\"CitationRef\\\"\\u003e13\\u003c/span\\u003e], limiting its usefulness for large-scale applications such as environmental monitoring or phytosanitary screening.\\u003c/p\\u003e \\u003cp\\u003eRecent advances in high-throughput sequencing (HTS) have enabled the rapid and sensitive analysis of many samples simultaneously. Amplicon-based HTS (i.e., metabarcoding) has become widely used for profiling fungal communities from various substrates, including plant tissues, to address questions related to their diversity and biosecurity [\\u003cspan citationid=\\\"CR14\\\" class=\\\"CitationRef\\\"\\u003e14\\u003c/span\\u003e, \\u003cspan citationid=\\\"CR15\\\" class=\\\"CitationRef\\\"\\u003e15\\u003c/span\\u003e]. This approach typically targets the nuclear ribosomal DNA (rDNA) internal transcribed spacer (ITS) region\\u0026mdash;a universal fungal barcode [\\u003cspan citationid=\\\"CR16\\\" class=\\\"CitationRef\\\"\\u003e16\\u003c/span\\u003e]. Sequencing platforms such as Illumina or PacBio are then used to characterise fungi from DNA extracted directly from bulk, host or environmental samples. However, DNA metabarcoding cannot distinguish living, metabolically active from dead or dormant organisms [\\u003cspan citationid=\\\"CR14\\\" class=\\\"CitationRef\\\"\\u003e14\\u003c/span\\u003e]. For example, Carini et al. [\\u003cspan citationid=\\\"CR17\\\" class=\\\"CitationRef\\\"\\u003e17\\u003c/span\\u003e] showed that over 40% of ITS sequences in soil belong to dead organisms, inflating diversity and skewing community profiles. This raises concerns when using DNA metabarcoding to study fungal responses to environmental change or in biosecurity [\\u003cspan citationid=\\\"CR17\\\" class=\\\"CitationRef\\\"\\u003e17\\u003c/span\\u003e, \\u003cspan citationid=\\\"CR18\\\" class=\\\"CitationRef\\\"\\u003e18\\u003c/span\\u003e], where only viable organisms pose a biosecurity threat that requires mitigation measures. Additionally, distinguishing dead from living organisms is essential for evaluating the effectiveness of phytosanitary treatments. Unlike DNA, RNA exists only in metabolically active cells and some viruses, and degrades rapidly after cell death [\\u003cspan citationid=\\\"CR19\\\" class=\\\"CitationRef\\\"\\u003e19\\u003c/span\\u003e, \\u003cspan citationid=\\\"CR20\\\" class=\\\"CitationRef\\\"\\u003e20\\u003c/span\\u003e]. In fungi, the rDNA operon is transcribed into a precursor ribosomal RNA (rRNA) that includes the ITS region. Targeting this precursor allows identification of fungi actively transcribing rRNA [\\u003cspan citationid=\\\"CR21\\\" class=\\\"CitationRef\\\"\\u003e21\\u003c/span\\u003e], resulting in a more accurate characterisation of the living community. While the potential of using RNA for improved biological monitoring has been recognised in environmental fields such as water research [\\u003cspan citationid=\\\"CR22\\\" class=\\\"CitationRef\\\"\\u003e22\\u003c/span\\u003e], the approach remains largely underexplored, including for detecting active fungi in plant samples.\\u003c/p\\u003e \\u003cp\\u003eStudies using DNA and RNA metabarcoding revealed different total and active microbial communities in temperate soil, wood [\\u003cspan citationid=\\\"CR23\\\" class=\\\"CitationRef\\\"\\u003e23\\u003c/span\\u003e, \\u003cspan citationid=\\\"CR24\\\" class=\\\"CitationRef\\\"\\u003e24\\u003c/span\\u003e], ballast water [\\u003cspan citationid=\\\"CR18\\\" class=\\\"CitationRef\\\"\\u003e18\\u003c/span\\u003e, \\u003cspan citationid=\\\"CR25\\\" class=\\\"CitationRef\\\"\\u003e25\\u003c/span\\u003e] and ocean [\\u003cspan citationid=\\\"CR26\\\" class=\\\"CitationRef\\\"\\u003e26\\u003c/span\\u003e]. However, other studies found no differences. For example, 92% of fungal taxa in \\u003cem\\u003ePicea abies\\u003c/em\\u003e needles overlapped between DNA and RNA datasets [\\u003cspan citationid=\\\"CR27\\\" class=\\\"CitationRef\\\"\\u003e27\\u003c/span\\u003e], and similar patterns were observed in Arctic soil [\\u003cspan citationid=\\\"CR28\\\" class=\\\"CitationRef\\\"\\u003e28\\u003c/span\\u003e]. The differences between total and active fungal communities may depend on environmental factors such as resource availability and temperature, which influence fungal metabolism [\\u003cspan citationid=\\\"CR29\\\" class=\\\"CitationRef\\\"\\u003e29\\u003c/span\\u003e, \\u003cspan citationid=\\\"CR30\\\" class=\\\"CitationRef\\\"\\u003e30\\u003c/span\\u003e] and DNA/RNA turnover [\\u003cspan citationid=\\\"CR31\\\" class=\\\"CitationRef\\\"\\u003e31\\u003c/span\\u003e]. Methodological biases may also contribute to the observed discrepancies, highlighting the need to assess how well these two approaches reflect fungal assemblages\\u0026mdash;especially when distinguishing active fungi from dead or dormant ones is critical for ecological inference and biosecurity.\\u003c/p\\u003e \\u003cp\\u003eIn this study, we characterised fungal communities associated with commercially traded seeds of three ecologically and economically important European tree species\\u0026mdash;common beech (\\u003cem\\u003eFagus sylvatica\\u003c/em\\u003e L.), silver fir (\\u003cem\\u003eAbies alba\\u003c/em\\u003e Mill.), and Scots pine (\\u003cem\\u003ePinus sylvestris\\u003c/em\\u003e L.)\\u0026mdash;using traditional culturing and DNA- and RNA-based metabarcoding. We assessed incidence-based and abundance-weighted alpha diversity (i.e., OTU richness and Inverse Simpson index, respectively) and beta diversity (i.e., S\\u0026oslash;rensen and Morisita-Horn dissimilarities, respectively) across tree species and molecular methods. We also compared taxonomic and lifestyle composition of fungal communities assessed with DNA and RNA metabarcoding within each tree species and compared fungal communities detected by each metabarcoding method with those from culturing. We hypothesised that DNA metabarcoding would reveal greater diversity than RNA metabarcoding, as RNA-based communities represent the active subset of the DNA-detected communities. Moreover, we hypothesised that fungi obtained by culturing would more closely resemble RNA-based communities than DNA-based ones, because culturing primarily recovers metabolically active fungi at the time of isolation.\\u003c/p\\u003e \"},{\"header\":\"Methods\",\"content\":\"\\u003cp\\u003eStudy material\\u003c/p\\u003e \\u003cp\\u003eA total of 16 seed lots (i.e., batches of seeds from specific locations and tree species; \\u003cem\\u003eF. sylvatica\\u003c/em\\u003e (n\\u0026thinsp;=\\u0026thinsp;6), \\u003cem\\u003eA. alba\\u003c/em\\u003e (n\\u0026thinsp;=\\u0026thinsp;5), \\u003cem\\u003eP. sylvestris\\u003c/em\\u003e (n\\u0026thinsp;=\\u0026thinsp;5)) were obtained in winter 2019 from European commercial seed suppliers. Seeds were stored at 4\\u0026deg;C until fungal assessment using culturing, DNA- and RNA-metabarcoding in spring 2020.\\u003c/p\\u003e \\u003cp\\u003eFungal assessment by culturing\\u003c/p\\u003e \\u003cp\\u003eSeeds were surface sterilized by subsequent immersion in 70% ethanol (1 min), 1% sodium hypochlorite (5 min), 70% ethanol (30 sec), and sterile water (30 sec) [\\u003cspan citationid=\\\"CR32\\\" class=\\\"CitationRef\\\"\\u003e32\\u003c/span\\u003e], then dried in a laminar flow hood. For fungal assessment, 90 seeds per seed lot were crushed with a sterile pestle and placed on nutrient media containing streptomycin (0.05 mg/mL). To capture a broad range of fungi, three media with varying nutrient content were used\\u0026mdash;water agar (WA 15 g/L; VWR Chemicals, Solon, Ohio, USA), malt extract agar (MEA; ME 30g/L; Duchefa Biochemie, Haarlem, the Netherlands \\u0026amp; WA 15 g/L; VWR Chemicals) and potato dextrose agar (PDA 39g/L; Merck KGaA, Darmstadt, Germany). Thirty seeds were plated per tree species x medium combination. Plates were checked every five days for 30 days. Emerging colonies were transferred to PDA (39 g/L; Merck KGaA) to obtain pure cultures. Isolates were grouped based on macromorphological traits (i.e., colour, texture, form, margin).\\u003c/p\\u003e \\u003cp\\u003eOne isolate per morphotype was selected for DNA extraction. Mycelia were scraped from agar, freeze-dried overnight, and DNA was extracted using the E.Z.N.A.\\u0026reg; SP Plant DNA Kit (Omega Bio-Tek, Norcross, Georgia, USA) following the manufacturer\\u0026rsquo;s instructions. DNA concentrations were measured using a DS-11 UV-Vis Spectrophotometer (DeNovix, Wilmington, Delaware, USA) and diluted to 5 ng/\\u0026micro;l. DNA extracts with concentrations below 5 ng/\\u0026micro;l were not diluted. The ITS region was amplified in 25 \\u0026micro;l reactions, containing 2 \\u0026micro;l DNA, 8.5 \\u0026micro;l nuclease-free H\\u003csub\\u003e2\\u003c/sub\\u003eO, 12.5 \\u0026micro;l DreamTaqPCR Master Mix (2X), and 1 \\u0026micro;l each of ITS1 (F: TCCGTAGGTGAACCTGCGG) and ITS4 (R: TCCTCCGCTTATTGATATGC) primers [\\u003cspan citationid=\\\"CR33\\\" class=\\\"CitationRef\\\"\\u003e33\\u003c/span\\u003e]. PCRs were run on an Eppendorf Mastercycler (Eppendorf, Hamburg, Germany) under the following conditions: 95\\u0026deg;C for 2 min; 35 cycles of 95\\u0026deg;C for 1 min; 55\\u0026deg;C for 45 sec; 72\\u0026deg;C for 1.5 min; final extension at 72\\u0026deg;C for 5 min. Amplification success was confirmed by gel electrophoresis. PCR products were purified and Sanger-sequenced by Macrogen-Europe (Amsterdam, the Netherlands) using the same primers as in PCRs.\\u003c/p\\u003e \\u003cp\\u003eForward and reverse sequences were trimmed and assembled into consensus sequences using CLC Main Workbench 26 (Qiagen, Aarhus, Denmark). Study sequences were compared against the National Centre for Biotechnology Information (NCBI) core nucleotide (nt) database using MEGABLAST [\\u003cspan citationid=\\\"CR34\\\" class=\\\"CitationRef\\\"\\u003e34\\u003c/span\\u003e, \\u003cspan citationid=\\\"CR35\\\" class=\\\"CitationRef\\\"\\u003e35\\u003c/span\\u003e], excluding sequences annotated as uncultured or environmental samples, to ensure reliable taxonomic assignment.\\u003c/p\\u003e \\u003cp\\u003eFungal assessment by metabarcoding\\u003c/p\\u003e \\u003cp\\u003eFungi were also assessed from 30 seeds per seed lot using DNA and RNA metabarcoding. Seeds were surface sterilised as for culturing, then ground in liquid nitrogen using RNase \\u003cem\\u003eAWAY\\u003c/em\\u003e\\u0026trade; Surface Decontaminant-treated mortars and pestles (Thermo Fisher Scientific, Waltham, USA). DNA was extracted using the E.Z.N.A.\\u0026reg; SP Plant DNA Kit (Omega Bio-Tek) following the manufacturer\\u0026rsquo;s instructions. DNA quantity was determined with a Qubit dsDNA BR Assay kit on a Qubit 3.0 Fluorometer (Thermo Fisher Scientific). Extracts were sent to Novogene (Cambridge, UK) for PCR, library preparation and sequencing. The ITS2 region was amplified using primers ITS3 (F: GCATCGATGAAGAACGCAGC) and ITS4 (R: TCCTCCGCTTATTGATATFC) [\\u003cspan citationid=\\\"CR33\\\" class=\\\"CitationRef\\\"\\u003e33\\u003c/span\\u003e], and libraries were sequenced on an Illumina MiSeq (2 x 250 bp).\\u003c/p\\u003e \\u003cp\\u003eRNA was extracted using the E.Z.N.A.\\u0026reg; Plant RNA kit (Omega Bio-Tek) following the manufacturer`s protocol for difficult samples, from the same ground seed pool used for DNA extractions. Multiple extractions per sample were pooled to obtain sufficient RNA for complementary DNA (i.e., cDNA) synthesis. After pooling, RNA quantity was measured using the Qubit\\u0026trade; RNA Broad Range (BR) Assay Kit with Qubit 3.0 Fluorometer (Thermo Fisher Scientific), and quality was assessed using the Agilent RNA 6000 Pico Kit and Agilent 2100 Bioanalyzer (Agilent, Santa Clara, California, USA). DNA contamination was removed using TURBO DNase (Invitrogen\\u0026trade;, Thermo Fisher Scientific) with a 60 min incubation (4 U/reaction), confirmed by 35-cycle ITS PCR (ITS1/ITS4 primers; [\\u003cspan citationid=\\\"CR33\\\" class=\\\"CitationRef\\\"\\u003e33\\u003c/span\\u003e]) as described above and the absence of bands in 1.2% TAE agarose gel electrophoresis. RNA was cleaned and concentrated with the RNeasy MinElute Cleanup Kit (QIAGEN, Hilden, Germany), and cDNA synthesised using the SuperScript\\u0026trade; IV First-Strand System (Invitrogen\\u0026trade;, Thermo Fisher Scientific) with random hexamer primers and RNase H treatment. The resulting cDNA was purified with AMPure XP beads (1:4 bead:sample ratio) and stored at \\u0026minus;\\u0026thinsp;80\\u0026deg;C. Samples were sent to BMKGENE (M\\u0026uuml;nster, Germany) for PCR, library preparation and Illumina NovaSeq 6000 sequencing (2 x 250 bp) of the ITS2 region using ITS3 and ITS4 primers [\\u003cspan citationid=\\\"CR33\\\" class=\\\"CitationRef\\\"\\u003e33\\u003c/span\\u003e].\\u003c/p\\u003e \\u003cp\\u003eSamples were demultiplexed by the sequencing facility, yielding separate forward and reverse FASTQ files for each sample and dataset (i.e., DNA and RNA). Bioinformatic processing was performed using a custom QIIME2-based pipeline [\\u003cspan citationid=\\\"CR36\\\" class=\\\"CitationRef\\\"\\u003e36\\u003c/span\\u003e]. Adapters and primers were removed with Cutadapt (i.e., error rate 0.1; [\\u003cspan citationid=\\\"CR37\\\" class=\\\"CitationRef\\\"\\u003e37\\u003c/span\\u003e]). Quality filtering, denoising, read merging and chimera removal were conducted using DADA2 [\\u003cspan citationid=\\\"CR38\\\" class=\\\"CitationRef\\\"\\u003e38\\u003c/span\\u003e]. Amplicon sequence variants of both DNA and RNA datasets were clustered into OTUs at 97% similarity using VSEARCH [\\u003cspan citationid=\\\"CR39\\\" class=\\\"CitationRef\\\"\\u003e39\\u003c/span\\u003e], and taxonomic assignment of the merged dataset was performed with a Na\\u0026iuml;ve Bayes classifier [\\u003cspan citationid=\\\"CR40\\\" class=\\\"CitationRef\\\"\\u003e40\\u003c/span\\u003e] in QIIME2 against the UNITE v8.2 database [\\u003cspan citationid=\\\"CR41\\\" class=\\\"CitationRef\\\"\\u003e41\\u003c/span\\u003e]. Taxonomic assignments at species level were not considered, because species identification based on the ITS2 region is often unreliable, particularly for certain fungal groups [\\u003cspan citationid=\\\"CR42\\\" class=\\\"CitationRef\\\"\\u003e42\\u003c/span\\u003e].\\u003c/p\\u003e \\u003cp\\u003eFungal lifestyle assignments\\u003c/p\\u003e \\u003cp\\u003eFungal lifestyles were assigned to OTUs based on genus-level matches using the \\u003cem\\u003eFungalTraits\\u003c/em\\u003e database [\\u003cspan citationid=\\\"CR43\\\" class=\\\"CitationRef\\\"\\u003e43\\u003c/span\\u003e], referencing the \\u0026ldquo;primary_lifestyle\\u0026rdquo; column. The \\\"plant pathogen\\\" category was retained as is, while the following classification adjustments were made:\\u003c/p\\u003e \\u003cp\\u003e1) \\\"ectomycorrhizal\\\", \\\"arbuscular_mycorrhizal\\\" = \\\"mycorrhiza\\\",\\u003c/p\\u003e \\u003cp\\u003e2) \\\"foliar_endophyte\\\", \\\"root_endophyte\\\" = \\\"endophyte\\\",\\u003c/p\\u003e \\u003cp\\u003e3) \\\"soil_saprotroph\\\", \\\"dung_saprotroph\\\", \\\"litter_saprotroph\\\", \\\"nectar/tap_saprotroph\\\",\\\"wood_saprotroph\\u0026rdquo;, \\u0026ldquo;pollen_saprotroph\\\", \\\"unspecified_saprotroph\\\" = \\\"saprotroph\\\",\\u003c/p\\u003e \\u003cp\\u003e4) \\\"mycoparasite\\\",\\\"animal_parasite\\\",\\\"epiphyte\\\",\\\"lichen_parasite\\\", \\\"algal_parasite\\u0026rdquo;, \\\"sooty_mold\\\", \\\"lichenized\\\" = \\\"other\\\".\\u003c/p\\u003e \\u003cp\\u003eStatistical analyses\\u003c/p\\u003e \\u003cp\\u003eAll analyses were performed using R Statistical software [\\u003cspan citationid=\\\"CR44\\\" class=\\\"CitationRef\\\"\\u003e44\\u003c/span\\u003e].\\u003c/p\\u003e \\u003cp\\u003eTo account for differences in sequencing depth (i.e., number of reads per sample ranged from 7,939 to 1,179,119 with a mean\\u0026thinsp;\\u0026plusmn;\\u0026thinsp;se of 48,663\\u0026thinsp;\\u0026plusmn;\\u0026thinsp;6,590) and retain all the samples, the dataset was rarefied to a minimum number of reads per sample (i.e., 7,939 reads in a DNA sample of \\u003cem\\u003eP. sylvestris\\u003c/em\\u003e; \\u003cem\\u003eRarefy\\u003c/em\\u003e function, GUniFrac package; [\\u003cspan citationid=\\\"CR45\\\" class=\\\"CitationRef\\\"\\u003e45\\u003c/span\\u003e]) prior to alpha and beta diversity analyses [\\u003cspan citationid=\\\"CR46\\\" class=\\\"CitationRef\\\"\\u003e46\\u003c/span\\u003e]. Alpha and beta diversity analyses were conducted using incidence-based and abundance weighted metrices. While incidence-based diversity measures capture overall diversity patterns by giving equal weight to all taxa regardless of abundance, abundance-weighted measures emphasize patterns driven by dominant taxa by assigning greater weight to abundant taxa.\\u003c/p\\u003e \\u003cp\\u003eAlpha diversity was assessed using OTU richness (i.e., incidence-based metric) and Inverse Simpson Index (i.e., abundance weighted metric; [\\u003cspan citationid=\\\"CR47\\\" class=\\\"CitationRef\\\"\\u003e47\\u003c/span\\u003e], calculated using the \\u003cem\\u003ehill_taxa\\u003c/em\\u003e function from the package hillR [\\u003cspan citationid=\\\"CR48\\\" class=\\\"CitationRef\\\"\\u003e48\\u003c/span\\u003e, \\u003cspan citationid=\\\"CR49\\\" class=\\\"CitationRef\\\"\\u003e49\\u003c/span\\u003e]). Alpha diversity metrics (i.e., response variables) were compared between explanatory variables\\u0026mdash;the methods (i.e., DNA and RNA metabarcoding) and tree species (i.e., \\u003cem\\u003eF. sylvatica\\u003c/em\\u003e, \\u003cem\\u003eA. alba\\u003c/em\\u003e and \\u003cem\\u003eP. sylvestris\\u003c/em\\u003e). To account for overdispersion, OTU richness was analysed using negative binomial models (i.e., \\u003cem\\u003eglm.nb\\u003c/em\\u003e function, MASS package; [\\u003cspan citationid=\\\"CR50\\\" class=\\\"CitationRef\\\"\\u003e50\\u003c/span\\u003e]). Inverse Simpson Index was modelled with a gamma distribution generalized linear model (GLM) using log link (i.e., \\u003cem\\u003eglm\\u003c/em\\u003e function, stats package; [\\u003cspan citationid=\\\"CR50\\\" class=\\\"CitationRef\\\"\\u003e50\\u003c/span\\u003e]) due to its non-integer nature. Models with and without the interaction term were compared based on Akaike information criterion (AIC). Models without the interaction term were retained (i.e., lower AIC) and further evaluated for factor significance (i.e., \\u003cem\\u003eAnova\\u003c/em\\u003e function, car package; [\\u003cspan citationid=\\\"CR51\\\" class=\\\"CitationRef\\\"\\u003e51\\u003c/span\\u003e]). Estimated marginal means for significant factors and pairwise comparisons were computed with the function \\u003cem\\u003eemmeans\\u003c/em\\u003e from the emmeans package [\\u003cspan citationid=\\\"CR52\\\" class=\\\"CitationRef\\\"\\u003e52\\u003c/span\\u003e].\\u003c/p\\u003e \\u003cp\\u003eBeta diversity was assessed using S\\u0026oslash;rensen dissimilarity (i.e., incidence-based metric), and Morisita-Horn dissimilarity (i.e., abundance-weighted metric). Pairwise distances between samples were calculated using the \\u003cem\\u003evegdist\\u003c/em\\u003e function from the vegan package [\\u003cspan citationid=\\\"CR53\\\" class=\\\"CitationRef\\\"\\u003e53\\u003c/span\\u003e]. Overall effects of method and tree species (i.e., explanatory variables) on fungal community composition (i.e., S\\u0026oslash;rensen and Morisita-Horn dissimilarity; response variables) were tested using Permutational multivariate ANOVA (i.e., PERMANOVA; adonis function, vegan package; [\\u003cspan citationid=\\\"CR53\\\" class=\\\"CitationRef\\\"\\u003e53\\u003c/span\\u003e]), without interaction term to match alpha diversity models. Differences in fungal community composition between the methods and tree species were visualized using non-metric multidimensional scaling (i.e., NMDS; \\u003cem\\u003emetaMDS\\u003c/em\\u003e function, vegan package, [\\u003cspan citationid=\\\"CR53\\\" class=\\\"CitationRef\\\"\\u003e53\\u003c/span\\u003e]).\\u003c/p\\u003e \\u003cp\\u003eAs PERMANOVA revealed significant differences in fungal communities among tree species (see Results), taxonomic composition was assessed separately for each tree species and method, using a non-rarefied dataset to preserve full abundance information across samples [\\u003cspan citationid=\\\"CR54\\\" class=\\\"CitationRef\\\"\\u003e54\\u003c/span\\u003e]. We focused on reads assigned to fungal genera, enabling inference of ecological roles and host associations. First, the number of reads associated with each distinct fungal genus was calculated for every method and tree species. The genera were then ranked by read count, and the ten most abundant genera per method and tree species were selected for further analysis. All remaining reads were grouped as either \\u0026ldquo;Other\\u0026rdquo; (i.e., reads assigned to genera outside the top ten) or \\u0026ldquo;Unassigned\\u0026rdquo; (i.e., reads not assigned to any genus). Relative abundances of the top ten genera, and genera classified as \\u0026ldquo;Other\\u0026rdquo; or \\u0026ldquo;Unassigned\\u0026rdquo; were plotted across the methods and tree species. In addition, indicator species analysis (i.e., \\u003cem\\u003emultipatt\\u003c/em\\u003e function, indicspecies package, [\\u003cspan citationid=\\\"CR55\\\" class=\\\"CitationRef\\\"\\u003e55\\u003c/span\\u003e]) was performed across all genera, irrespective of their read count, to identify those strongly associated with either the DNA or RNA dataset for each tree species.\\u003c/p\\u003e \\u003cp\\u003eVariation in fungal lifestyles was evaluated using all non-rarefied reads assigned to a specific lifestyle category or classified as \\u0026ldquo;Unassigned\\u0026rdquo;. For each method and tree species, the relative abundance of reads associated with each lifestyle were calculated and plotted.\\u003c/p\\u003e \\u003cp\\u003eFinally, we identified fungal genera shared between culturing and each metabarcoding dataset, as well as those unique to each dataset. This was first done for all tree species together to assess the overall trend, and for each tree species separately, to assess the differences across tree species.\\u003c/p\\u003e\"},{\"header\":\"Results\",\"content\":\"\\u003cp\\u003eOf the 3,167,108 non-singleton reads assigned to 1,797 OTUs in the merged dataset comprising DNA and RNA reads, approximately 50% were fungal (i.e., 1,557,227 reads and 942 OTUs). Rarefaction curves indicated that the sequencing depth was sufficient to capture all fungal OTUs across all samples (Supplementary Fig.\\u0026nbsp;1).\\u003c/p\\u003e \\u003cp\\u003eMore than 60% of the fungal reads originated from \\u003cem\\u003eA. alba\\u003c/em\\u003e samples, while the remaining reads were similarly distributed between the other two species (Table\\u0026nbsp;\\u003cspan refid=\\\"Tab1\\\" class=\\\"InternalRef\\\"\\u003e1\\u003c/span\\u003e). Moreover, the variation in the number of reads per sample was smaller for \\u003cem\\u003eA. alba\\u003c/em\\u003e than for the other two species (Supplementary Fig.\\u0026nbsp;1). Overall, a similar number of reads was obtained from both sequencing methods, although this varied among tree species. Read counts were comparable between methods in \\u003cem\\u003eF. sylvatica\\u003c/em\\u003e and \\u003cem\\u003eA. alba\\u003c/em\\u003e, whereas the \\u003cem\\u003eP. sylvestris\\u003c/em\\u003e DNA dataset was almost twice the size of the RNA dataset (Table\\u0026nbsp;\\u003cspan refid=\\\"Tab1\\\" class=\\\"InternalRef\\\"\\u003e1\\u003c/span\\u003e). \\u003cem\\u003ePinus sylvestris\\u003c/em\\u003e samples contained twice as many OTUs as \\u003cem\\u003eA. alba\\u003c/em\\u003e, with \\u003cem\\u003eF. sylvatica\\u003c/em\\u003e falling between the two (Table\\u0026nbsp;\\u003cspan refid=\\\"Tab1\\\" class=\\\"InternalRef\\\"\\u003e1\\u003c/span\\u003e). The total number of OTUs in all samples of each of the three tree species was three to four times higher in the RNA dataset than in the DNA dataset (Table\\u0026nbsp;\\u003cspan refid=\\\"Tab1\\\" class=\\\"InternalRef\\\"\\u003e1\\u003c/span\\u003e).\\u003c/p\\u003e \\u003cp\\u003eThe rarefied dataset that was used for the alpha and beta diversity analyses consisted of 254,048 non-singleton fungal reads equally distributed across all samples, as all samples were rarefied to the same number of reads (i.e., 7,939 reads). Overall, these reads were assigned to 872 OTUs (Table\\u0026nbsp;\\u003cspan refid=\\\"Tab1\\\" class=\\\"InternalRef\\\"\\u003e1\\u003c/span\\u003e). Although 70 OTUs (i.e., around 7%) were lost during rarefaction, the number of OTUs in the RNA dataset remained three to four times higher than that in the DNA dataset (Table\\u0026nbsp;\\u003cspan refid=\\\"Tab1\\\" class=\\\"InternalRef\\\"\\u003e1\\u003c/span\\u003e).\\u003c/p\\u003e \\u003cp\\u003e \\u003cdiv class=\\\"gridtable\\\"\\u003e\\u003ctable float=\\\"Yes\\\" id=\\\"Tab1\\\" border=\\\"1\\\"\\u003e \\u003ccaption language=\\\"En\\\"\\u003e \\u003cdiv class=\\\"CaptionNumber\\\"\\u003eTable 1\\u003c/div\\u003e \\u003cdiv class=\\\"CaptionContent\\\"\\u003e \\u003cp\\u003e\\u003cb\\u003eFungal read and OTU counts across three tree species and DNA and RNA metabarcoding datasets\\u003c/b\\u003e. Number of fungal reads and OTUs in samples belonging to the three tree species (i.e., \\u003cem\\u003eF. sylvatica\\u003c/em\\u003e, \\u003cem\\u003eA. alba\\u003c/em\\u003e, \\u003cem\\u003eP. sylvestris\\u003c/em\\u003e; numbers of seed lots are indicated in brackets) analysed by DNA and RNA metabarcoding are indicated. Numbers are shown for the DNA and RNA datasets and for the combined dataset (i.e., all data) considering all fungal reads and rarefied fungal reads (in the brackets).\\u003c/p\\u003e \\u003c/div\\u003e \\u003c/caption\\u003e \\u003ccolgroup cols=\\\"7\\\"\\u003e \\u003cdiv align=\\\"left\\\" class=\\\"colspec\\\" colname=\\\"c1\\\" colnum=\\\"1\\\"\\u003e\\u003c/div\\u003e \\u003cdiv align=\\\"left\\\" class=\\\"colspec\\\" colname=\\\"c2\\\" colnum=\\\"2\\\"\\u003e\\u003c/div\\u003e \\u003cdiv align=\\\"left\\\" class=\\\"colspec\\\" colname=\\\"c3\\\" colnum=\\\"3\\\"\\u003e\\u003c/div\\u003e \\u003cdiv align=\\\"left\\\" class=\\\"colspec\\\" colname=\\\"c4\\\" colnum=\\\"4\\\"\\u003e\\u003c/div\\u003e \\u003cdiv align=\\\"left\\\" class=\\\"colspec\\\" colname=\\\"c5\\\" colnum=\\\"5\\\"\\u003e\\u003c/div\\u003e \\u003cdiv align=\\\"left\\\" class=\\\"colspec\\\" colname=\\\"c6\\\" colnum=\\\"6\\\"\\u003e\\u003c/div\\u003e \\u003cdiv align=\\\"left\\\" class=\\\"colspec\\\" colname=\\\"c7\\\" colnum=\\\"7\\\"\\u003e\\u003c/div\\u003e \\u003cthead\\u003e \\u003ctr\\u003e \\u003cth align=\\\"left\\\" colname=\\\"c1\\\"\\u003e\\u0026nbsp;\\u003c/th\\u003e \\u003cth align=\\\"left\\\" colspan=\\\"3\\\" nameend=\\\"c4\\\" namest=\\\"c2\\\"\\u003e \\u003cp\\u003eReads\\u003c/p\\u003e \\u003c/th\\u003e \\u003cth align=\\\"left\\\" colspan=\\\"3\\\" nameend=\\\"c7\\\" namest=\\\"c5\\\"\\u003e \\u003cp\\u003eOTUs\\u003c/p\\u003e \\u003c/th\\u003e \\u003c/tr\\u003e \\u003c/thead\\u003e \\u003ctbody\\u003e \\u003ctr\\u003e \\u003ctd align=\\\"left\\\" colname=\\\"c1\\\"\\u003e \\u003cp\\u003eTree species\\u003c/p\\u003e \\u003c/td\\u003e \\u003ctd align=\\\"left\\\" colname=\\\"c2\\\"\\u003e \\u003cp\\u003eDNA\\u003c/p\\u003e \\u003c/td\\u003e \\u003ctd align=\\\"left\\\" colname=\\\"c3\\\"\\u003e \\u003cp\\u003eRNA\\u003c/p\\u003e \\u003c/td\\u003e \\u003ctd align=\\\"left\\\" colname=\\\"c4\\\"\\u003e \\u003cp\\u003eAll data\\u003c/p\\u003e \\u003c/td\\u003e \\u003ctd align=\\\"left\\\" colname=\\\"c5\\\"\\u003e \\u003cp\\u003eDNA\\u003c/p\\u003e \\u003c/td\\u003e \\u003ctd align=\\\"left\\\" colname=\\\"c6\\\"\\u003e \\u003cp\\u003eRNA\\u003c/p\\u003e \\u003c/td\\u003e \\u003ctd align=\\\"left\\\" colname=\\\"c7\\\"\\u003e \\u003cp\\u003eAll data\\u003c/p\\u003e \\u003c/td\\u003e \\u003c/tr\\u003e \\u003ctr\\u003e \\u003ctd align=\\\"left\\\" colname=\\\"c1\\\" morerows=\\\"1\\\" rowspan=\\\"2\\\"\\u003e \\u003cp\\u003e\\u003cem\\u003eF. sylvatica\\u003c/em\\u003e (n\\u0026thinsp;=\\u0026thinsp;6)\\u003c/p\\u003e \\u003c/td\\u003e \\u003ctd align=\\\"left\\\" colname=\\\"c2\\\"\\u003e \\u003cp\\u003e105,396\\u003c/p\\u003e \\u003c/td\\u003e \\u003ctd align=\\\"left\\\" colname=\\\"c3\\\"\\u003e \\u003cp\\u003e177,420\\u003c/p\\u003e \\u003c/td\\u003e \\u003ctd align=\\\"left\\\" colname=\\\"c4\\\"\\u003e \\u003cp\\u003e282,816\\u003c/p\\u003e \\u003c/td\\u003e \\u003ctd align=\\\"left\\\" colname=\\\"c5\\\"\\u003e \\u003cp\\u003e108\\u003c/p\\u003e \\u003c/td\\u003e \\u003ctd align=\\\"left\\\" colname=\\\"c6\\\"\\u003e \\u003cp\\u003e361\\u003c/p\\u003e \\u003c/td\\u003e \\u003ctd align=\\\"left\\\" colname=\\\"c7\\\"\\u003e \\u003cp\\u003e417\\u003c/p\\u003e \\u003c/td\\u003e \\u003c/tr\\u003e \\u003ctr\\u003e \\u003ctd align=\\\"left\\\" colname=\\\"c2\\\"\\u003e \\u003cp\\u003e(47,634)\\u003c/p\\u003e \\u003c/td\\u003e \\u003ctd align=\\\"left\\\" colname=\\\"c3\\\"\\u003e \\u003cp\\u003e(47,634)\\u003c/p\\u003e \\u003c/td\\u003e \\u003ctd align=\\\"left\\\" colname=\\\"c4\\\"\\u003e \\u003cp\\u003e(95,268)\\u003c/p\\u003e \\u003c/td\\u003e \\u003ctd align=\\\"left\\\" colname=\\\"c5\\\"\\u003e \\u003cp\\u003e(102)\\u003c/p\\u003e \\u003c/td\\u003e \\u003ctd align=\\\"left\\\" colname=\\\"c6\\\"\\u003e \\u003cp\\u003e(342)\\u003c/p\\u003e \\u003c/td\\u003e \\u003ctd align=\\\"left\\\" colname=\\\"c7\\\"\\u003e \\u003cp\\u003e(392)\\u003c/p\\u003e \\u003c/td\\u003e \\u003c/tr\\u003e \\u003ctr\\u003e \\u003ctd align=\\\"left\\\" colname=\\\"c1\\\" morerows=\\\"1\\\" rowspan=\\\"2\\\"\\u003e \\u003cp\\u003e\\u003cem\\u003eA. alba\\u003c/em\\u003e (n\\u0026thinsp;=\\u0026thinsp;5)\\u003c/p\\u003e \\u003c/td\\u003e \\u003ctd align=\\\"left\\\" colname=\\\"c2\\\"\\u003e \\u003cp\\u003e444,834\\u003c/p\\u003e \\u003c/td\\u003e \\u003ctd align=\\\"left\\\" colname=\\\"c3\\\"\\u003e \\u003cp\\u003e522,881\\u003c/p\\u003e \\u003c/td\\u003e \\u003ctd align=\\\"left\\\" colname=\\\"c4\\\"\\u003e \\u003cp\\u003e967,715\\u003c/p\\u003e \\u003c/td\\u003e \\u003ctd align=\\\"left\\\" colname=\\\"c5\\\"\\u003e \\u003cp\\u003e93\\u003c/p\\u003e \\u003c/td\\u003e \\u003ctd align=\\\"left\\\" colname=\\\"c6\\\"\\u003e \\u003cp\\u003e232\\u003c/p\\u003e \\u003c/td\\u003e \\u003ctd align=\\\"left\\\" colname=\\\"c7\\\"\\u003e \\u003cp\\u003e277\\u003c/p\\u003e \\u003c/td\\u003e \\u003c/tr\\u003e \\u003ctr\\u003e \\u003ctd align=\\\"left\\\" colname=\\\"c2\\\"\\u003e \\u003cp\\u003e(39,695)\\u003c/p\\u003e \\u003c/td\\u003e \\u003ctd align=\\\"left\\\" colname=\\\"c3\\\"\\u003e \\u003cp\\u003e(39,695)\\u003c/p\\u003e \\u003c/td\\u003e \\u003ctd align=\\\"left\\\" colname=\\\"c4\\\"\\u003e \\u003cp\\u003e(79,390)\\u003c/p\\u003e \\u003c/td\\u003e \\u003ctd align=\\\"left\\\" colname=\\\"c5\\\"\\u003e \\u003cp\\u003e(72)\\u003c/p\\u003e \\u003c/td\\u003e \\u003ctd align=\\\"left\\\" colname=\\\"c6\\\"\\u003e \\u003cp\\u003e(173)\\u003c/p\\u003e \\u003c/td\\u003e \\u003ctd align=\\\"left\\\" colname=\\\"c7\\\"\\u003e \\u003cp\\u003e(202)\\u003c/p\\u003e \\u003c/td\\u003e \\u003c/tr\\u003e \\u003ctr\\u003e \\u003ctd align=\\\"left\\\" colname=\\\"c1\\\" morerows=\\\"1\\\" rowspan=\\\"2\\\"\\u003e \\u003cp\\u003e\\u003cem\\u003eP. sylvestris\\u003c/em\\u003e (n\\u0026thinsp;=\\u0026thinsp;5)\\u003c/p\\u003e \\u003c/td\\u003e \\u003ctd align=\\\"left\\\" colname=\\\"c2\\\"\\u003e \\u003cp\\u003e212,820\\u003c/p\\u003e \\u003c/td\\u003e \\u003ctd align=\\\"left\\\" colname=\\\"c3\\\"\\u003e \\u003cp\\u003e93,876\\u003c/p\\u003e \\u003c/td\\u003e \\u003ctd align=\\\"left\\\" colname=\\\"c4\\\"\\u003e \\u003cp\\u003e306,696\\u003c/p\\u003e \\u003c/td\\u003e \\u003ctd align=\\\"left\\\" colname=\\\"c5\\\"\\u003e \\u003cp\\u003e152\\u003c/p\\u003e \\u003c/td\\u003e \\u003ctd align=\\\"left\\\" colname=\\\"c6\\\"\\u003e \\u003cp\\u003e519\\u003c/p\\u003e \\u003c/td\\u003e \\u003ctd align=\\\"left\\\" colname=\\\"c7\\\"\\u003e \\u003cp\\u003e628\\u003c/p\\u003e \\u003c/td\\u003e \\u003c/tr\\u003e \\u003ctr\\u003e \\u003ctd align=\\\"left\\\" colname=\\\"c2\\\"\\u003e \\u003cp\\u003e(39,695)\\u003c/p\\u003e \\u003c/td\\u003e \\u003ctd align=\\\"left\\\" colname=\\\"c3\\\"\\u003e \\u003cp\\u003e(39,695)\\u003c/p\\u003e \\u003c/td\\u003e \\u003ctd align=\\\"left\\\" colname=\\\"c4\\\"\\u003e \\u003cp\\u003e(79,390)\\u003c/p\\u003e \\u003c/td\\u003e \\u003ctd align=\\\"left\\\" colname=\\\"c5\\\"\\u003e \\u003cp\\u003e(139)\\u003c/p\\u003e \\u003c/td\\u003e \\u003ctd align=\\\"left\\\" colname=\\\"c6\\\"\\u003e \\u003cp\\u003e(512)\\u003c/p\\u003e \\u003c/td\\u003e \\u003ctd align=\\\"left\\\" colname=\\\"c7\\\"\\u003e \\u003cp\\u003e(610)\\u003c/p\\u003e \\u003c/td\\u003e \\u003c/tr\\u003e \\u003ctr\\u003e \\u003ctd align=\\\"left\\\" colname=\\\"c1\\\" morerows=\\\"1\\\" rowspan=\\\"2\\\"\\u003e \\u003cp\\u003eTotal\\u003c/p\\u003e \\u003c/td\\u003e \\u003ctd align=\\\"left\\\" colname=\\\"c2\\\"\\u003e \\u003cp\\u003e763,050\\u003c/p\\u003e \\u003c/td\\u003e \\u003ctd align=\\\"left\\\" colname=\\\"c3\\\"\\u003e \\u003cp\\u003e794,177\\u003c/p\\u003e \\u003c/td\\u003e \\u003ctd align=\\\"left\\\" colname=\\\"c4\\\"\\u003e \\u003cp\\u003e1,557,227\\u003c/p\\u003e \\u003c/td\\u003e \\u003ctd align=\\\"left\\\" colname=\\\"c5\\\"\\u003e \\u003cp\\u003e247\\u003c/p\\u003e \\u003c/td\\u003e \\u003ctd align=\\\"left\\\" colname=\\\"c6\\\"\\u003e \\u003cp\\u003e801\\u003c/p\\u003e \\u003c/td\\u003e \\u003ctd align=\\\"left\\\" colname=\\\"c7\\\"\\u003e \\u003cp\\u003e942\\u003c/p\\u003e \\u003c/td\\u003e \\u003c/tr\\u003e \\u003ctr\\u003e \\u003ctd align=\\\"left\\\" colname=\\\"c2\\\"\\u003e \\u003cp\\u003e(127,024)\\u003c/p\\u003e \\u003c/td\\u003e \\u003ctd align=\\\"left\\\" colname=\\\"c3\\\"\\u003e \\u003cp\\u003e(127,024)\\u003c/p\\u003e \\u003c/td\\u003e \\u003ctd align=\\\"left\\\" colname=\\\"c4\\\"\\u003e \\u003cp\\u003e(254,048)\\u003c/p\\u003e \\u003c/td\\u003e \\u003ctd align=\\\"left\\\" colname=\\\"c5\\\"\\u003e \\u003cp\\u003e(215)\\u003c/p\\u003e \\u003c/td\\u003e \\u003ctd align=\\\"left\\\" colname=\\\"c6\\\"\\u003e \\u003cp\\u003e(512)\\u003c/p\\u003e \\u003c/td\\u003e \\u003ctd align=\\\"left\\\" colname=\\\"c7\\\"\\u003e \\u003cp\\u003e(872)\\u003c/p\\u003e \\u003c/td\\u003e \\u003c/tr\\u003e \\u003c/tbody\\u003e \\u003c/colgroup\\u003e \\u003c/table\\u003e\\u003c/div\\u003e \\u003c/p\\u003e \\u003cp\\u003eAlpha diversity\\u003c/p\\u003e \\u003cp\\u003eRNA metabarcoding revealed approximately twice as many OTUs per sample compared to DNA metabarcoding (Fig.\\u0026nbsp;\\u003cspan refid=\\\"Fig1\\\" class=\\\"InternalRef\\\"\\u003e1\\u003c/span\\u003eA; χ\\u0026sup2; = 25.95, df\\u0026thinsp;=\\u0026thinsp;1, \\u003cem\\u003ep\\u003c/em\\u003e\\u0026thinsp;\\u0026lt;\\u0026thinsp;0.001). This pattern persisted for abundance-weighted alpha diversity (i.e., Inverse Simpson Index; χ\\u0026sup2; = 17.90, df\\u0026thinsp;=\\u0026thinsp;1, \\u003cem\\u003ep\\u003c/em\\u003e\\u0026thinsp;\\u0026lt;\\u0026thinsp;0.001), although the mean values of the Inverse Simpson Index were roughly ten times lower than those of OTU richness (Supplementary Fig.\\u0026nbsp;2A).\\u003c/p\\u003e \\u003cp\\u003eOverall, OTU richness per sample differed significantly between tree species (Fig.\\u0026nbsp;\\u003cspan refid=\\\"Fig1\\\" class=\\\"InternalRef\\\"\\u003e1\\u003c/span\\u003eB; χ\\u0026sup2; = 8.18, df\\u0026thinsp;=\\u0026thinsp;1, \\u003cem\\u003ep\\u003c/em\\u003e\\u0026thinsp;\\u0026lt;\\u0026thinsp;0.05), with the lowest richness observed in \\u003cem\\u003eA. alba\\u003c/em\\u003e samples and the highest in \\u003cem\\u003eP. sylvestris\\u003c/em\\u003e samples. A similar pattern was observed for the abundance-weighted diversity (χ\\u0026sup2; = 15.12, df\\u0026thinsp;=\\u0026thinsp;2, \\u003cem\\u003ep\\u003c/em\\u003e\\u0026thinsp;\\u0026lt;\\u0026thinsp;0.001) which was around ten times lower than OTU richness\\u0026mdash;\\u003cem\\u003eF. sylvatica\\u003c/em\\u003e and \\u003cem\\u003eP. sylvestris\\u003c/em\\u003e exhibited more than twice the diversity of \\u003cem\\u003eA. alba\\u003c/em\\u003e (Supplementary Fig.\\u0026nbsp;2B).\\u003c/p\\u003e \\u003cp\\u003e \\u003c/p\\u003e \\u003cp\\u003eBeta diversity\\u003c/p\\u003e \\u003cp\\u003eWhen all OTUs were considered, regardless of their abundance, DNA and RNA metabarcoding revealed significantly different fungal communities associated with the same seeds (Fig.\\u0026nbsp;\\u003cspan refid=\\\"Fig2\\\" class=\\\"InternalRef\\\"\\u003e2\\u003c/span\\u003eA; F\\u0026thinsp;=\\u0026thinsp;11.49, R\\u003csup\\u003e2\\u003c/sup\\u003e\\u0026thinsp;=\\u0026thinsp;0.25, p\\u0026thinsp;\\u0026lt;\\u0026thinsp;0.001). Of the total 872 OTUs identified in the rarefied dataset, only 100 (~\\u0026thinsp;11%) were shared between the datasets, while around 13% and 75% were unique for the DNA and RNA dataset, respectively (Fig.\\u0026nbsp;\\u003cspan refid=\\\"Fig2\\\" class=\\\"InternalRef\\\"\\u003e2\\u003c/span\\u003eB). The proportion of OTUs shared between the two metabarcoding methods further differed across tree species: approximately 7% of OTUs were shared in \\u003cem\\u003eP. sylvestris\\u003c/em\\u003e, compared to 13% in \\u003cem\\u003eF. sylvatica\\u003c/em\\u003e and 21% \\u003cem\\u003ein A. alba\\u003c/em\\u003e (Supplementary Fig.\\u0026nbsp;3A\\u0026ndash;C). Comparable results were obtained using the non-rarefied dataset (Supplementary Fig.\\u0026nbsp;4A-E). Significant differences in the composition of dominant fungal communities\\u0026mdash;assessed by Morisita-Horn dissimilarity, which gives greater weight to OTUs with higher read counts\\u0026mdash;were also observed between DNA and RNA datasets (Supplementary Fig.\\u0026nbsp;5A; F\\u0026thinsp;=\\u0026thinsp;1.97, R\\u003csup\\u003e2\\u003c/sup\\u003e\\u0026thinsp;=\\u0026thinsp;0.05, p\\u0026thinsp;\\u0026lt;\\u0026thinsp;0.05). However, as indicated by R\\u003csup\\u003e2\\u003c/sup\\u003e values, the relative importance of the method in explaining the differences in fungal community composition decreased from 25% to 5% when an abundance-weighted beta diversity measure was used suggesting that the differences between methods were largely driven by rare OTUs.\\u003c/p\\u003e \\u003cp\\u003eFungal communities also differed among tree species, both when analysing the entire community (Fig.\\u0026nbsp;\\u003cspan refid=\\\"Fig2\\\" class=\\\"InternalRef\\\"\\u003e2\\u003c/span\\u003eC-D; F\\u0026thinsp;=\\u0026thinsp;2.84, R\\u003csup\\u003e2\\u003c/sup\\u003e\\u0026thinsp;=\\u0026thinsp;0.13, p\\u0026thinsp;\\u0026lt;\\u0026thinsp;0.001) and when focussing on dominant taxa (Supplementary Fig.\\u0026nbsp;5B; F\\u0026thinsp;=\\u0026thinsp;1.75, R\\u003csup\\u003e2\\u003c/sup\\u003e\\u0026thinsp;=\\u0026thinsp;0.28, p\\u0026thinsp;\\u0026lt;\\u0026thinsp;0.001). The proportion of variation explained by tree species increased from 13% to nearly 28% when switching from incidence-based to abundance-weighted metrics, respectively, as indicated by R\\u003csup\\u003e2\\u003c/sup\\u003e values, suggesting that taxa with higher read counts were more host-specific than rare taxa.\\u003c/p\\u003e \\u003cp\\u003e \\u003c/p\\u003e \\u003cp\\u003eTaxonomic compositional differences between methods and tree species\\u003c/p\\u003e \\u003cp\\u003eOut of a total of 1,557,227 fungal reads and 942 fungal OTUs in the non-rarefied dataset, the majority of reads and OTUs were assigned at the class, order, and family levels. At the genus level, around 72% of reads and 64% of OTUs were classified to 338 distinct genera (Supplementary Table\\u0026nbsp;1).\\u003c/p\\u003e \\u003cp\\u003eThe proportion of fungal reads and OTUs assigned across taxonomic levels was broadly consistent between DNA and RNA datasets (Supplementary Table\\u0026nbsp;1). The vast majority of reads and OTUs were assigned at high taxonomic levels (i.e., class, order, family) in both datasets. Approximately 74% of reads and 68% of OTUs in the DNA dataset, and 71% of reads and 65% of OTUs in the RNA dataset, were assigned to 110 and 305 unique genera, respectively. Seventy-seven genera (~\\u0026thinsp;23%) appeared in both datasets, while 33 (~\\u0026thinsp;10%) and 228 (~\\u0026thinsp;67%) were unique for DNA and RNA dataset, respectively.\\u003c/p\\u003e \\u003cp\\u003eApproximately 85% of reads were assigned to genera in both DNA and RNA datasets of \\u003cem\\u003eF. sylvatica\\u003c/em\\u003e, compared to 73% in \\u003cem\\u003eA. alba\\u003c/em\\u003e and 61% in \\u003cem\\u003eP. sylvestris\\u003c/em\\u003e (Fig.\\u0026nbsp;\\u003cspan refid=\\\"Fig3\\\" class=\\\"InternalRef\\\"\\u003e3\\u003c/span\\u003e). The proportion of reads classified as \\u0026ldquo;Other\\u0026rdquo; (i.e., reads assigned to genera outside the top ten most abundant genera per method and tree species) was around only 3% in the \\u003cem\\u003eA. alba\\u003c/em\\u003e DNA and RNA, and \\u003cem\\u003eP. sylvestris\\u003c/em\\u003e DNA datasets. Higher proportion of reads classified as \\u0026ldquo;Other\\u0026rdquo; was found in the \\u003cem\\u003eP. sylvestris\\u003c/em\\u003e RNA dataset (~\\u0026thinsp;14%) and, especially, in both \\u003cem\\u003eF. sylvatica\\u003c/em\\u003e datasets (i.e., ~\\u0026thinsp;16% in the DNA and ~\\u0026thinsp;26% in the RNA dataset; Fig.\\u0026nbsp;\\u003cspan refid=\\\"Fig3\\\" class=\\\"InternalRef\\\"\\u003e3\\u003c/span\\u003e).\\u003c/p\\u003e \\u003cp\\u003eTaxonomic composition of dominant genera (i.e., the top ten most abundant genera per method and tree species) varied across tree species and sequencing methods (Fig.\\u0026nbsp;\\u003cspan refid=\\\"Fig3\\\" class=\\\"InternalRef\\\"\\u003e3\\u003c/span\\u003e). The only genus consistently identified among the top ten genera across all methods and tree species was \\u003cem\\u003eAspergillus\\u003c/em\\u003e. Among the top ten genera identified for each method and tree species, five were shared between the DNA and RNA datasets in \\u003cem\\u003eF. sylvatica\\u003c/em\\u003e (i.e., \\u003cem\\u003eAlternaria\\u003c/em\\u003e, \\u003cem\\u003eApiognomonia\\u003c/em\\u003e, \\u003cem\\u003eAspergillus\\u003c/em\\u003e, \\u003cem\\u003eDiaporthe\\u003c/em\\u003e, \\u003cem\\u003eGibberella\\u003c/em\\u003e), and in \\u003cem\\u003eA. alba\\u003c/em\\u003e (i.e., \\u003cem\\u003eAspergillus\\u003c/em\\u003e, \\u003cem\\u003eDiplodia\\u003c/em\\u003e, \\u003cem\\u003eHormonema\\u003c/em\\u003e, \\u003cem\\u003eLambertella\\u003c/em\\u003e, \\u003cem\\u003ePenicillium\\u003c/em\\u003e), and two in \\u003cem\\u003eP. sylvestris\\u003c/em\\u003e (i.e., \\u003cem\\u003eAspergillus\\u003c/em\\u003e, \\u003cem\\u003eHormonema\\u003c/em\\u003e). Those shared genera accounted for a high proportion of reads in each dataset (i.e., \\u003cem\\u003eF. sylvatica\\u003c/em\\u003e: 53% in DNA and 32% in RNA; \\u003cem\\u003eA. alba\\u003c/em\\u003e: 60% in DNA and 55% in RNA; \\u003cem\\u003eP. sylvestris\\u003c/em\\u003e: 44% in DNA and 40% in RNA).\\u003c/p\\u003e \\u003cp\\u003e \\u003c/p\\u003e \\u003cp\\u003eIndicator species analysis revealed significant, dataset-specific genus associations with either DNA or RNA datasets (Supplementary Table\\u0026nbsp;2). Across tree species, the genera \\u003cem\\u003eMycocentrospora\\u003c/em\\u003e, \\u003cem\\u003eParaleptosphaeria\\u003c/em\\u003e, and \\u003cem\\u003ePhoma\\u003c/em\\u003e were associated with the DNA datasets, whereas \\u003cem\\u003eFusarium\\u003c/em\\u003e, \\u003cem\\u003eKazachstania\\u003c/em\\u003e, and \\u003cem\\u003eVishniacozyma\\u003c/em\\u003e were associated with the RNA datasets, although with varying relative abundances. In \\u003cem\\u003eF. sylvatica\\u003c/em\\u003e, four genera were associated with the DNA and 17 with the RNA dataset. Among these, several indicator genera ranked among the top ten genera per method in \\u003cem\\u003eF. sylvatica\\u003c/em\\u003e\\u0026mdash;\\u003cem\\u003eParaleptosphaeria\\u003c/em\\u003e was an indicator of the DNA, whereas \\u003cem\\u003eVishniacozyma\\u003c/em\\u003e and \\u003cem\\u003eLophodermium\\u003c/em\\u003e were representative of the RNA dataset. In \\u003cem\\u003eA. alba\\u003c/em\\u003e, nine genera were indicators of the DNA dataset and ten of the RNA dataset. Among the top ten genera per method in \\u003cem\\u003eA. alba\\u003c/em\\u003e, \\u003cem\\u003eCaloscypha\\u003c/em\\u003e was strongly linked to DNA, while \\u003cem\\u003eFusarium\\u003c/em\\u003e, \\u003cem\\u003eLophodermium\\u003c/em\\u003e and \\u003cem\\u003eNeocatenulostroma\\u003c/em\\u003e emerged as indicators of the RNA dataset. For \\u003cem\\u003eP. sylvestris\\u003c/em\\u003e, 12 genera were associated with the DNA dataset and nine with the RNA dataset. Within the top ten genera per method in \\u003cem\\u003eP. sylvestris\\u003c/em\\u003e, \\u003cem\\u003eCaloscypha\\u003c/em\\u003e, \\u003cem\\u003eFibulochlamys\\u003c/em\\u003e, \\u003cem\\u003ePodospora\\u003c/em\\u003e, and \\u003cem\\u003eLambertella\\u003c/em\\u003e were DNA-associated, whereas \\u003cem\\u003eColletotrichum\\u003c/em\\u003e, \\u003cem\\u003eKazachstania\\u003c/em\\u003e, and \\u003cem\\u003eVishniacozyma\\u003c/em\\u003e were strongly associated with the RNA dataset. Additional genus-level associations were observed but fell outside the top ten most abundant genera per method and species (Supplementary Table\\u0026nbsp;2).\\u003c/p\\u003e \\u003cp\\u003eLifestyle compositional differences between methods and tree species\\u003c/p\\u003e \\u003cp\\u003eOf the total reads and OTUs in the non-rarified dataset (i.e., 1,557,227 reads and 942 OTUs), a primary lifestyle could be assigned to 1,128,634 reads (74%) representing 596 OTUs (64%). Among these, most fungal reads and OTUs belonged to saprotrophs (72% of reads, 59% of OTUs), followed by plant pathogens (25% of reads, 24% of OTUs), while mycorrhizal fungi, endophytes, and taxa categorized as \\u0026ldquo;Other\\u0026rdquo; each accounted for only a small fraction of the reads and OTUs (Supplementary Table\\u0026nbsp;3).\\u003c/p\\u003e \\u003cp\\u003eThe proportion of fungal reads and OTUs assigned to a lifestyle was largely consistent across DNA and RNA datasets\\u0026mdash;around 74% and 71% out of 763,050 and 794,177 reads and 68% and 64% out of 247 and 801 OTUs in DNA and RNA dataset could be assigned to a lifestyle. Moreover, the proportions of reads and OTUs assigned to each lifestyle were largely similar between the metabarcoding datasets (Supplementary Table\\u0026nbsp;3).\\u003c/p\\u003e \\u003cp\\u003eApproximately 85%, 75%, and 65% of reads from \\u003cem\\u003eF. sylvatica\\u003c/em\\u003e, \\u003cem\\u003eA. alba\\u003c/em\\u003e, and \\u003cem\\u003eP. sylvestris\\u003c/em\\u003e, respectively, were classified into fungal lifestyles, with varying proportions between the two datasets (Fig.\\u0026nbsp;\\u003cspan refid=\\\"Fig4\\\" class=\\\"InternalRef\\\"\\u003e4\\u003c/span\\u003e). In \\u003cem\\u003eF. sylvatica\\u003c/em\\u003e, plant pathogens accounted for 57% of DNA reads and 39% of RNA reads, while saprotrophs were represented with 29% of DNA reads and 38% of RNA reads. Saprotrophs were the most abundant guild in \\u003cem\\u003eA. alba\\u003c/em\\u003e (60%) and \\u003cem\\u003eP. sylvestris\\u003c/em\\u003e (52%) in both datasets. Plant pathogens constituted 15% and 12% of DNA and RNA reads in \\u003cem\\u003eA. alba\\u003c/em\\u003e samples and 7% and 9% in \\u003cem\\u003eP. sylvestris\\u003c/em\\u003e DNA and RNA dataset, respectively. In \\u003cem\\u003eF. sylvatica\\u003c/em\\u003e and \\u003cem\\u003eA. alba\\u003c/em\\u003e, mycorrhizal fungi (i.e., 740 and 392 reads) and endophytes (i.e., 217 and 233 reads) were exclusively detected in the RNA dataset. Unlike \\u003cem\\u003eF. sylvatica\\u003c/em\\u003e and \\u003cem\\u003eA. alba\\u003c/em\\u003e, \\u003cem\\u003eP. sylvestris\\u003c/em\\u003e contained endophytes in both DNA and RNA datasets at similar read counts (i.e., 85 and 51, respectively). Mycorrhizal fungi in \\u003cem\\u003eP. sylvestris\\u003c/em\\u003e were detected in both datasets but with only three reads in the DNA compared to 201 reads in the RNA dataset.\\u003c/p\\u003e \\u003cp\\u003e \\u003c/p\\u003e \\u003cp\\u003eComparison of metabarcoding datasets with traditional culturing data\\u003c/p\\u003e \\u003cp\\u003eA total of 1,553 fungal isolates were obtained across all seed samples and those were assigned to 203 morphotypes. Out of 203 representative isolates (i.e., one representative isolate per morphotype) which were selected for sequencing, 145 yielded good-quality sequences. Of those, 121 sequences were assigned to 31 distinct genera, compared to 110 and 305 genera identified through DNA and RNA metabarcoding, respectively. Of the 31 genera identified through culturing, ten (~\\u0026thinsp;32%) were found exclusively in the culturing datasets. In contrast, 17 (~\\u0026thinsp;55%) of the cultured genera were also detected in both DNA and RNA datasets and two additional cultured genera were detected in each of the metabarcoding datasets. Although the overlap between culturing and metabarcoding was broadly similar\\u0026mdash;with most cultured genera appearing in both metabarcoding datasets\\u0026mdash;the RNA dataset contained a higher number of unique genera, whereas the DNA dataset comprised of fewer unique taxa (Fig.\\u0026nbsp;\\u003cspan refid=\\\"Fig5\\\" class=\\\"InternalRef\\\"\\u003e5\\u003c/span\\u003e).\\u003c/p\\u003e \\u003cp\\u003e \\u003c/p\\u003e \\u003cp\\u003eMost of the isolates originated from \\u003cem\\u003eF. sylvatica\\u003c/em\\u003e seeds (i.e., 1,300 isolates representing 145 distinct morphotypes). Out of the 114 representative isolates which yielded good quality sequences 91 representative isolates could be assigned to one of 24 fungal genera (i.e., corresponding to 976 isolates and 91 morphotypes). Almost half of the cultured fungal genera (i.e., eleven) were detected by both DNA and RNA metabarcoding approaches, while an additional two (i.e., \\u003cem\\u003eOliveonia\\u003c/em\\u003e, \\u003cem\\u003ePhacidium\\u003c/em\\u003e) and three genera (i.e., \\u003cem\\u003eCladosporium\\u003c/em\\u003e, \\u003cem\\u003eChaetomium\\u003c/em\\u003e, \\u003cem\\u003eClonostachys\\u003c/em\\u003e) were exclusive to the DNA or RNA datasets, respectively. Notably, eight genera were identified solely through culturing (i.e., \\u003cem\\u003eDiscosia\\u003c/em\\u003e, \\u003cem\\u003eEpicoccum\\u003c/em\\u003e, \\u003cem\\u003eXylaria\\u003c/em\\u003e, \\u003cem\\u003eDidymosphaeria\\u003c/em\\u003e, \\u003cem\\u003ePlagiostoma\\u003c/em\\u003e, \\u003cem\\u003eBiscogniauxia\\u003c/em\\u003e, \\u003cem\\u003eRadulidium\\u003c/em\\u003e, \\u003cem\\u003eBoeremia\\u003c/em\\u003e; Table\\u0026nbsp;\\u003cspan refid=\\\"Tab2\\\" class=\\\"InternalRef\\\"\\u003e2\\u003c/span\\u003e., Supplementary Fig.\\u0026nbsp;6).\\u003c/p\\u003e \\u003cp\\u003eFrom \\u003cem\\u003eA. alba\\u003c/em\\u003e seeds, 231 isolates representing 44 morphotypes were obtained from seeds. Good quality sequences were obtained from 19 representative isolates, covering 97 isolates belonging to 19 morphotypes. All representative sequences could be assigned to one of the eight genera, five of which were detected in both DNA and RNA datasets (i.e., \\u003cem\\u003ePenicillium\\u003c/em\\u003e, \\u003cem\\u003eTrichoderma, Fusarium, Aspergillus, Talaromyces\\u003c/em\\u003e), while three were exclusive to the culturing dataset (i.e., \\u003cem\\u003eSydowia\\u003c/em\\u003e, \\u003cem\\u003eMucor\\u003c/em\\u003e, \\u003cem\\u003eAkanthomyces\\u003c/em\\u003e; Table\\u0026nbsp;\\u003cspan refid=\\\"Tab2\\\" class=\\\"InternalRef\\\"\\u003e2\\u003c/span\\u003e., Supplementary Fig.\\u0026nbsp;6).\\u003c/p\\u003e \\u003cp\\u003eFrom \\u003cem\\u003eP. sylvestris\\u003c/em\\u003e seeds, 22 isolates representing 14 morphotypes were obtained. Good quality sequences were obtained from 12 representative isolates (morphotypes), covering 20 isolates. All but one representative isolate was assigned to one of seven genera. Three genera were shared between DNA and RNA datasets (i.e., \\u003cem\\u003ePenicillium\\u003c/em\\u003e, \\u003cem\\u003eFusarium\\u003c/em\\u003e, \\u003cem\\u003eAspergillus\\u003c/em\\u003e), two were shared with RNA dataset (i.e., \\u003cem\\u003eConiochaeta\\u003c/em\\u003e, \\u003cem\\u003ePseudopithomyces\\u003c/em\\u003e) and two were detected only via culturing (i.e., \\u003cem\\u003eBoeremia\\u003c/em\\u003e, \\u003cem\\u003eSydowia\\u003c/em\\u003e; Table\\u0026nbsp;\\u003cspan refid=\\\"Tab2\\\" class=\\\"InternalRef\\\"\\u003e2\\u003c/span\\u003e., Supplementary Fig.\\u0026nbsp;6).\\u003c/p\\u003e \\u003cp\\u003e \\u003cdiv class=\\\"gridtable\\\"\\u003e\\u003ctable float=\\\"Yes\\\" id=\\\"Tab2\\\" border=\\\"1\\\"\\u003e \\u003ccaption language=\\\"En\\\"\\u003e \\u003cdiv class=\\\"CaptionNumber\\\"\\u003eTable 2\\u003c/div\\u003e \\u003cdiv class=\\\"CaptionContent\\\"\\u003e \\u003cp\\u003e\\u003cb\\u003eFungal genera cultured from three tree species and their detection in DNA/RNA metabarcoding datasets\\u003c/b\\u003e Fungal genera detected in the culturing data of the three tree species (i.e., \\u003cem\\u003eF. sylvatica\\u003c/em\\u003e, \\u003cem\\u003eA.alba\\u003c/em\\u003e, \\u003cem\\u003eP. sylvestris)\\u003c/em\\u003e are listed. Number of isolates belonging to each genus and number of morphotypes those isolates belong to is indicated, as well as the information about if the genera were detected in either of the two metabarcoding datasets (i.e., DNA, RNA).\\u003c/p\\u003e \\u003c/div\\u003e \\u003c/caption\\u003e \\u003ccolgroup cols=\\\"5\\\"\\u003e \\u003cdiv align=\\\"left\\\" class=\\\"colspec\\\" colname=\\\"c1\\\" colnum=\\\"1\\\"\\u003e\\u003c/div\\u003e \\u003cdiv align=\\\"left\\\" class=\\\"colspec\\\" colname=\\\"c2\\\" colnum=\\\"2\\\"\\u003e\\u003c/div\\u003e \\u003cdiv align=\\\"left\\\" class=\\\"colspec\\\" colname=\\\"c3\\\" colnum=\\\"3\\\"\\u003e\\u003c/div\\u003e \\u003cdiv align=\\\"left\\\" class=\\\"colspec\\\" colname=\\\"c4\\\" colnum=\\\"4\\\"\\u003e\\u003c/div\\u003e \\u003cdiv align=\\\"left\\\" class=\\\"colspec\\\" colname=\\\"c5\\\" colnum=\\\"5\\\"\\u003e\\u003c/div\\u003e \\u003cthead\\u003e \\u003ctr\\u003e \\u003cth align=\\\"left\\\" colname=\\\"c1\\\"\\u003e \\u003cp\\u003eGenus\\u003c/p\\u003e \\u003c/th\\u003e \\u003cth align=\\\"left\\\" colname=\\\"c2\\\"\\u003e \\u003cp\\u003eNumber of isolates\\u003c/p\\u003e \\u003c/th\\u003e \\u003cth align=\\\"left\\\" colname=\\\"c3\\\"\\u003e \\u003cp\\u003eNumber of morphotypes\\u003c/p\\u003e \\u003c/th\\u003e \\u003cth align=\\\"left\\\" colname=\\\"c4\\\"\\u003e \\u003cp\\u003ePresent in DNA dataset\\u003c/p\\u003e \\u003c/th\\u003e \\u003cth align=\\\"left\\\" colname=\\\"c5\\\"\\u003e \\u003cp\\u003ePresent in RNA dataset\\u003c/p\\u003e \\u003c/th\\u003e \\u003c/tr\\u003e \\u003c/thead\\u003e \\u003ctbody\\u003e \\u003ctr\\u003e \\u003ctd align=\\\"left\\\" colspan=\\\"5\\\" nameend=\\\"c5\\\" namest=\\\"c1\\\"\\u003e \\u003cp\\u003e\\u003cb\\u003eFagus sylvatica\\u003c/b\\u003e\\u003c/p\\u003e \\u003c/td\\u003e \\u003c/tr\\u003e \\u003ctr\\u003e \\u003ctd align=\\\"left\\\" colname=\\\"c1\\\"\\u003e \\u003cp\\u003e\\u003cem\\u003eAlternaria\\u003c/em\\u003e\\u003c/p\\u003e \\u003c/td\\u003e \\u003ctd align=\\\"left\\\" colname=\\\"c2\\\"\\u003e \\u003cp\\u003e318\\u003c/p\\u003e \\u003c/td\\u003e \\u003ctd align=\\\"left\\\" colname=\\\"c3\\\"\\u003e \\u003cp\\u003e13\\u003c/p\\u003e \\u003c/td\\u003e \\u003ctd align=\\\"left\\\" colname=\\\"c4\\\"\\u003e \\u003cp\\u003ex\\u003c/p\\u003e \\u003c/td\\u003e \\u003ctd align=\\\"left\\\" colname=\\\"c5\\\"\\u003e \\u003cp\\u003ex\\u003c/p\\u003e \\u003c/td\\u003e \\u003c/tr\\u003e \\u003ctr\\u003e \\u003ctd align=\\\"left\\\" colname=\\\"c1\\\"\\u003e \\u003cp\\u003e\\u003cem\\u003eApiognomonia\\u003c/em\\u003e\\u003c/p\\u003e \\u003c/td\\u003e \\u003ctd align=\\\"left\\\" colname=\\\"c2\\\"\\u003e \\u003cp\\u003e244\\u003c/p\\u003e \\u003c/td\\u003e \\u003ctd align=\\\"left\\\" colname=\\\"c3\\\"\\u003e \\u003cp\\u003e18\\u003c/p\\u003e \\u003c/td\\u003e \\u003ctd align=\\\"left\\\" colname=\\\"c4\\\"\\u003e \\u003cp\\u003ex\\u003c/p\\u003e \\u003c/td\\u003e \\u003ctd align=\\\"left\\\" colname=\\\"c5\\\"\\u003e \\u003cp\\u003ex\\u003c/p\\u003e \\u003c/td\\u003e \\u003c/tr\\u003e \\u003ctr\\u003e \\u003ctd align=\\\"left\\\" colname=\\\"c1\\\"\\u003e \\u003cp\\u003e\\u003cem\\u003eFusarium\\u003c/em\\u003e\\u003c/p\\u003e \\u003c/td\\u003e \\u003ctd align=\\\"left\\\" colname=\\\"c2\\\"\\u003e \\u003cp\\u003e125\\u003c/p\\u003e \\u003c/td\\u003e \\u003ctd align=\\\"left\\\" colname=\\\"c3\\\"\\u003e \\u003cp\\u003e13\\u003c/p\\u003e \\u003c/td\\u003e \\u003ctd align=\\\"left\\\" colname=\\\"c4\\\"\\u003e \\u003cp\\u003ex\\u003c/p\\u003e \\u003c/td\\u003e \\u003ctd align=\\\"left\\\" colname=\\\"c5\\\"\\u003e \\u003cp\\u003ex\\u003c/p\\u003e \\u003c/td\\u003e \\u003c/tr\\u003e \\u003ctr\\u003e \\u003ctd align=\\\"left\\\" colname=\\\"c1\\\"\\u003e \\u003cp\\u003e\\u003cem\\u003eDiscosia\\u003c/em\\u003e\\u003c/p\\u003e \\u003c/td\\u003e \\u003ctd align=\\\"left\\\" colname=\\\"c2\\\"\\u003e \\u003cp\\u003e57\\u003c/p\\u003e \\u003c/td\\u003e \\u003ctd align=\\\"left\\\" colname=\\\"c3\\\"\\u003e \\u003cp\\u003e5\\u003c/p\\u003e \\u003c/td\\u003e \\u003ctd align=\\\"left\\\" colname=\\\"c4\\\"\\u003e\\u0026nbsp;\\u003c/td\\u003e \\u003ctd align=\\\"left\\\" colname=\\\"c5\\\"\\u003e\\u0026nbsp;\\u003c/td\\u003e \\u003c/tr\\u003e \\u003ctr\\u003e \\u003ctd align=\\\"left\\\" colname=\\\"c1\\\"\\u003e \\u003cp\\u003e\\u003cem\\u003eDiaporthe\\u003c/em\\u003e\\u003c/p\\u003e \\u003c/td\\u003e \\u003ctd align=\\\"left\\\" colname=\\\"c2\\\"\\u003e \\u003cp\\u003e42\\u003c/p\\u003e \\u003c/td\\u003e \\u003ctd align=\\\"left\\\" colname=\\\"c3\\\"\\u003e \\u003cp\\u003e4\\u003c/p\\u003e \\u003c/td\\u003e \\u003ctd align=\\\"left\\\" colname=\\\"c4\\\"\\u003e \\u003cp\\u003ex\\u003c/p\\u003e \\u003c/td\\u003e \\u003ctd align=\\\"left\\\" colname=\\\"c5\\\"\\u003e \\u003cp\\u003ex\\u003c/p\\u003e \\u003c/td\\u003e \\u003c/tr\\u003e \\u003ctr\\u003e \\u003ctd align=\\\"left\\\" colname=\\\"c1\\\"\\u003e \\u003cp\\u003e\\u003cem\\u003eEpicoccum\\u003c/em\\u003e\\u003c/p\\u003e \\u003c/td\\u003e \\u003ctd align=\\\"left\\\" colname=\\\"c2\\\"\\u003e \\u003cp\\u003e33\\u003c/p\\u003e \\u003c/td\\u003e \\u003ctd align=\\\"left\\\" colname=\\\"c3\\\"\\u003e \\u003cp\\u003e5\\u003c/p\\u003e \\u003c/td\\u003e \\u003ctd align=\\\"left\\\" colname=\\\"c4\\\"\\u003e\\u0026nbsp;\\u003c/td\\u003e \\u003ctd align=\\\"left\\\" colname=\\\"c5\\\"\\u003e\\u0026nbsp;\\u003c/td\\u003e \\u003c/tr\\u003e \\u003ctr\\u003e \\u003ctd align=\\\"left\\\" colname=\\\"c1\\\"\\u003e \\u003cp\\u003e\\u003cem\\u003eCladosporium\\u003c/em\\u003e\\u003c/p\\u003e \\u003c/td\\u003e \\u003ctd align=\\\"left\\\" colname=\\\"c2\\\"\\u003e \\u003cp\\u003e25\\u003c/p\\u003e \\u003c/td\\u003e \\u003ctd align=\\\"left\\\" colname=\\\"c3\\\"\\u003e \\u003cp\\u003e4\\u003c/p\\u003e \\u003c/td\\u003e \\u003ctd align=\\\"left\\\" colname=\\\"c4\\\"\\u003e\\u0026nbsp;\\u003c/td\\u003e \\u003ctd align=\\\"left\\\" colname=\\\"c5\\\"\\u003e \\u003cp\\u003ex\\u003c/p\\u003e \\u003c/td\\u003e \\u003c/tr\\u003e \\u003ctr\\u003e \\u003ctd align=\\\"left\\\" colname=\\\"c1\\\"\\u003e \\u003cp\\u003e\\u003cem\\u003ePenicillium\\u003c/em\\u003e\\u003c/p\\u003e \\u003c/td\\u003e \\u003ctd align=\\\"left\\\" colname=\\\"c2\\\"\\u003e \\u003cp\\u003e24\\u003c/p\\u003e \\u003c/td\\u003e \\u003ctd align=\\\"left\\\" colname=\\\"c3\\\"\\u003e \\u003cp\\u003e5\\u003c/p\\u003e \\u003c/td\\u003e \\u003ctd align=\\\"left\\\" colname=\\\"c4\\\"\\u003e \\u003cp\\u003ex\\u003c/p\\u003e \\u003c/td\\u003e \\u003ctd align=\\\"left\\\" colname=\\\"c5\\\"\\u003e \\u003cp\\u003ex\\u003c/p\\u003e \\u003c/td\\u003e \\u003c/tr\\u003e \\u003ctr\\u003e \\u003ctd align=\\\"left\\\" colname=\\\"c1\\\"\\u003e \\u003cp\\u003e\\u003cem\\u003eAureobasidium\\u003c/em\\u003e\\u003c/p\\u003e \\u003c/td\\u003e \\u003ctd align=\\\"left\\\" colname=\\\"c2\\\"\\u003e \\u003cp\\u003e14\\u003c/p\\u003e \\u003c/td\\u003e \\u003ctd align=\\\"left\\\" colname=\\\"c3\\\"\\u003e \\u003cp\\u003e3\\u003c/p\\u003e \\u003c/td\\u003e \\u003ctd align=\\\"left\\\" colname=\\\"c4\\\"\\u003e \\u003cp\\u003ex\\u003c/p\\u003e \\u003c/td\\u003e \\u003ctd align=\\\"left\\\" colname=\\\"c5\\\"\\u003e \\u003cp\\u003ex\\u003c/p\\u003e \\u003c/td\\u003e \\u003c/tr\\u003e \\u003ctr\\u003e \\u003ctd align=\\\"left\\\" colname=\\\"c1\\\"\\u003e \\u003cp\\u003e\\u003cem\\u003eCeratobasidium\\u003c/em\\u003e\\u003c/p\\u003e \\u003c/td\\u003e \\u003ctd align=\\\"left\\\" colname=\\\"c2\\\"\\u003e \\u003cp\\u003e10\\u003c/p\\u003e \\u003c/td\\u003e \\u003ctd align=\\\"left\\\" colname=\\\"c3\\\"\\u003e \\u003cp\\u003e1\\u003c/p\\u003e \\u003c/td\\u003e \\u003ctd align=\\\"left\\\" colname=\\\"c4\\\"\\u003e \\u003cp\\u003ex\\u003c/p\\u003e \\u003c/td\\u003e \\u003ctd align=\\\"left\\\" colname=\\\"c5\\\"\\u003e \\u003cp\\u003ex\\u003c/p\\u003e \\u003c/td\\u003e \\u003c/tr\\u003e \\u003ctr\\u003e \\u003ctd align=\\\"left\\\" colname=\\\"c1\\\"\\u003e \\u003cp\\u003e\\u003cem\\u003eOliveonia\\u003c/em\\u003e\\u003c/p\\u003e \\u003c/td\\u003e \\u003ctd align=\\\"left\\\" colname=\\\"c2\\\"\\u003e \\u003cp\\u003e10\\u003c/p\\u003e \\u003c/td\\u003e \\u003ctd align=\\\"left\\\" colname=\\\"c3\\\"\\u003e \\u003cp\\u003e1\\u003c/p\\u003e \\u003c/td\\u003e \\u003ctd align=\\\"left\\\" colname=\\\"c4\\\"\\u003e \\u003cp\\u003ex\\u003c/p\\u003e \\u003c/td\\u003e \\u003ctd align=\\\"left\\\" colname=\\\"c5\\\"\\u003e\\u0026nbsp;\\u003c/td\\u003e \\u003c/tr\\u003e \\u003ctr\\u003e \\u003ctd align=\\\"left\\\" colname=\\\"c1\\\"\\u003e \\u003cp\\u003e\\u003cem\\u003eClonostachys\\u003c/em\\u003e\\u003c/p\\u003e \\u003c/td\\u003e \\u003ctd align=\\\"left\\\" colname=\\\"c2\\\"\\u003e \\u003cp\\u003e9\\u003c/p\\u003e \\u003c/td\\u003e \\u003ctd align=\\\"left\\\" colname=\\\"c3\\\"\\u003e \\u003cp\\u003e4\\u003c/p\\u003e \\u003c/td\\u003e \\u003ctd align=\\\"left\\\" colname=\\\"c4\\\"\\u003e\\u0026nbsp;\\u003c/td\\u003e \\u003ctd align=\\\"left\\\" colname=\\\"c5\\\"\\u003e \\u003cp\\u003ex\\u003c/p\\u003e \\u003c/td\\u003e \\u003c/tr\\u003e \\u003ctr\\u003e \\u003ctd align=\\\"left\\\" colname=\\\"c1\\\"\\u003e \\u003cp\\u003e\\u003cem\\u003eTrichothecium\\u003c/em\\u003e\\u003c/p\\u003e \\u003c/td\\u003e \\u003ctd align=\\\"left\\\" colname=\\\"c2\\\"\\u003e \\u003cp\\u003e9\\u003c/p\\u003e \\u003c/td\\u003e \\u003ctd align=\\\"left\\\" colname=\\\"c3\\\"\\u003e \\u003cp\\u003e1\\u003c/p\\u003e \\u003c/td\\u003e \\u003ctd align=\\\"left\\\" colname=\\\"c4\\\"\\u003e \\u003cp\\u003ex\\u003c/p\\u003e \\u003c/td\\u003e \\u003ctd align=\\\"left\\\" colname=\\\"c5\\\"\\u003e \\u003cp\\u003ex\\u003c/p\\u003e \\u003c/td\\u003e \\u003c/tr\\u003e \\u003ctr\\u003e \\u003ctd align=\\\"left\\\" colname=\\\"c1\\\"\\u003e \\u003cp\\u003e\\u003cem\\u003eSeimatosporium\\u003c/em\\u003e\\u003c/p\\u003e \\u003c/td\\u003e \\u003ctd align=\\\"left\\\" colname=\\\"c2\\\"\\u003e \\u003cp\\u003e8\\u003c/p\\u003e \\u003c/td\\u003e \\u003ctd align=\\\"left\\\" colname=\\\"c3\\\"\\u003e \\u003cp\\u003e2\\u003c/p\\u003e \\u003c/td\\u003e \\u003ctd align=\\\"left\\\" colname=\\\"c4\\\"\\u003e \\u003cp\\u003ex\\u003c/p\\u003e \\u003c/td\\u003e \\u003ctd align=\\\"left\\\" colname=\\\"c5\\\"\\u003e \\u003cp\\u003ex\\u003c/p\\u003e \\u003c/td\\u003e \\u003c/tr\\u003e \\u003ctr\\u003e \\u003ctd align=\\\"left\\\" colname=\\\"c1\\\"\\u003e \\u003cp\\u003e\\u003cem\\u003eXylaria\\u003c/em\\u003e\\u003c/p\\u003e \\u003c/td\\u003e \\u003ctd align=\\\"left\\\" colname=\\\"c2\\\"\\u003e \\u003cp\\u003e8\\u003c/p\\u003e \\u003c/td\\u003e \\u003ctd align=\\\"left\\\" colname=\\\"c3\\\"\\u003e \\u003cp\\u003e2\\u003c/p\\u003e \\u003c/td\\u003e \\u003ctd align=\\\"left\\\" colname=\\\"c4\\\"\\u003e\\u0026nbsp;\\u003c/td\\u003e \\u003ctd align=\\\"left\\\" colname=\\\"c5\\\"\\u003e\\u0026nbsp;\\u003c/td\\u003e \\u003c/tr\\u003e \\u003ctr\\u003e \\u003ctd align=\\\"left\\\" colname=\\\"c1\\\"\\u003e \\u003cp\\u003e\\u003cem\\u003eDidymosphaeria\\u003c/em\\u003e\\u003c/p\\u003e \\u003c/td\\u003e \\u003ctd align=\\\"left\\\" colname=\\\"c2\\\"\\u003e \\u003cp\\u003e6\\u003c/p\\u003e \\u003c/td\\u003e \\u003ctd align=\\\"left\\\" colname=\\\"c3\\\"\\u003e \\u003cp\\u003e1\\u003c/p\\u003e \\u003c/td\\u003e \\u003ctd align=\\\"left\\\" colname=\\\"c4\\\"\\u003e\\u0026nbsp;\\u003c/td\\u003e \\u003ctd align=\\\"left\\\" colname=\\\"c5\\\"\\u003e\\u0026nbsp;\\u003c/td\\u003e \\u003c/tr\\u003e \\u003ctr\\u003e \\u003ctd align=\\\"left\\\" colname=\\\"c1\\\"\\u003e \\u003cp\\u003e\\u003cem\\u003eMuriphaeosphaeria\\u003c/em\\u003e\\u003c/p\\u003e \\u003c/td\\u003e \\u003ctd align=\\\"left\\\" colname=\\\"c2\\\"\\u003e \\u003cp\\u003e6\\u003c/p\\u003e \\u003c/td\\u003e \\u003ctd align=\\\"left\\\" colname=\\\"c3\\\"\\u003e \\u003cp\\u003e1\\u003c/p\\u003e \\u003c/td\\u003e \\u003ctd align=\\\"left\\\" colname=\\\"c4\\\"\\u003e \\u003cp\\u003ex\\u003c/p\\u003e \\u003c/td\\u003e \\u003ctd align=\\\"left\\\" colname=\\\"c5\\\"\\u003e \\u003cp\\u003ex\\u003c/p\\u003e \\u003c/td\\u003e \\u003c/tr\\u003e \\u003ctr\\u003e \\u003ctd align=\\\"left\\\" colname=\\\"c1\\\"\\u003e \\u003cp\\u003e\\u003cem\\u003eChaetomium\\u003c/em\\u003e\\u003c/p\\u003e \\u003c/td\\u003e \\u003ctd align=\\\"left\\\" colname=\\\"c2\\\"\\u003e \\u003cp\\u003e5\\u003c/p\\u003e \\u003c/td\\u003e \\u003ctd align=\\\"left\\\" colname=\\\"c3\\\"\\u003e \\u003cp\\u003e1\\u003c/p\\u003e \\u003c/td\\u003e \\u003ctd align=\\\"left\\\" colname=\\\"c4\\\"\\u003e\\u0026nbsp;\\u003c/td\\u003e \\u003ctd align=\\\"left\\\" colname=\\\"c5\\\"\\u003e \\u003cp\\u003ex\\u003c/p\\u003e \\u003c/td\\u003e \\u003c/tr\\u003e \\u003ctr\\u003e \\u003ctd align=\\\"left\\\" colname=\\\"c1\\\"\\u003e \\u003cp\\u003e\\u003cem\\u003eBiscogniauxia\\u003c/em\\u003e\\u003c/p\\u003e \\u003c/td\\u003e \\u003ctd align=\\\"left\\\" colname=\\\"c2\\\"\\u003e \\u003cp\\u003e5\\u003c/p\\u003e \\u003c/td\\u003e \\u003ctd align=\\\"left\\\" colname=\\\"c3\\\"\\u003e \\u003cp\\u003e1\\u003c/p\\u003e \\u003c/td\\u003e \\u003ctd align=\\\"left\\\" colname=\\\"c4\\\"\\u003e\\u0026nbsp;\\u003c/td\\u003e \\u003ctd align=\\\"left\\\" colname=\\\"c5\\\"\\u003e\\u0026nbsp;\\u003c/td\\u003e \\u003c/tr\\u003e \\u003ctr\\u003e \\u003ctd align=\\\"left\\\" colname=\\\"c1\\\"\\u003e \\u003cp\\u003e\\u003cem\\u003eRadulidium\\u003c/em\\u003e\\u003c/p\\u003e \\u003c/td\\u003e \\u003ctd align=\\\"left\\\" colname=\\\"c2\\\"\\u003e \\u003cp\\u003e5\\u003c/p\\u003e \\u003c/td\\u003e \\u003ctd align=\\\"left\\\" colname=\\\"c3\\\"\\u003e \\u003cp\\u003e1\\u003c/p\\u003e \\u003c/td\\u003e \\u003ctd align=\\\"left\\\" colname=\\\"c4\\\"\\u003e\\u0026nbsp;\\u003c/td\\u003e \\u003ctd align=\\\"left\\\" colname=\\\"c5\\\"\\u003e\\u0026nbsp;\\u003c/td\\u003e \\u003c/tr\\u003e \\u003ctr\\u003e \\u003ctd align=\\\"left\\\" colname=\\\"c1\\\"\\u003e \\u003cp\\u003e\\u003cem\\u003eNeosetophoma\\u003c/em\\u003e\\u003c/p\\u003e \\u003c/td\\u003e \\u003ctd align=\\\"left\\\" colname=\\\"c2\\\"\\u003e \\u003cp\\u003e4\\u003c/p\\u003e \\u003c/td\\u003e \\u003ctd align=\\\"left\\\" colname=\\\"c3\\\"\\u003e \\u003cp\\u003e2\\u003c/p\\u003e \\u003c/td\\u003e \\u003ctd align=\\\"left\\\" colname=\\\"c4\\\"\\u003e \\u003cp\\u003ex\\u003c/p\\u003e \\u003c/td\\u003e \\u003ctd align=\\\"left\\\" colname=\\\"c5\\\"\\u003e \\u003cp\\u003ex\\u003c/p\\u003e \\u003c/td\\u003e \\u003c/tr\\u003e \\u003ctr\\u003e \\u003ctd align=\\\"left\\\" colname=\\\"c1\\\"\\u003e \\u003cp\\u003e\\u003cem\\u003ePlagiostoma\\u003c/em\\u003e\\u003c/p\\u003e \\u003c/td\\u003e \\u003ctd align=\\\"left\\\" colname=\\\"c2\\\"\\u003e \\u003cp\\u003e3\\u003c/p\\u003e \\u003c/td\\u003e \\u003ctd align=\\\"left\\\" colname=\\\"c3\\\"\\u003e \\u003cp\\u003e1\\u003c/p\\u003e \\u003c/td\\u003e \\u003ctd align=\\\"left\\\" colname=\\\"c4\\\"\\u003e\\u0026nbsp;\\u003c/td\\u003e \\u003ctd align=\\\"left\\\" colname=\\\"c5\\\"\\u003e\\u0026nbsp;\\u003c/td\\u003e \\u003c/tr\\u003e \\u003ctr\\u003e \\u003ctd align=\\\"left\\\" colname=\\\"c1\\\"\\u003e \\u003cp\\u003e\\u003cem\\u003eBoeremia\\u003c/em\\u003e\\u003c/p\\u003e \\u003c/td\\u003e \\u003ctd align=\\\"left\\\" colname=\\\"c2\\\"\\u003e \\u003cp\\u003e3\\u003c/p\\u003e \\u003c/td\\u003e \\u003ctd align=\\\"left\\\" colname=\\\"c3\\\"\\u003e \\u003cp\\u003e1\\u003c/p\\u003e \\u003c/td\\u003e \\u003ctd align=\\\"left\\\" colname=\\\"c4\\\"\\u003e\\u0026nbsp;\\u003c/td\\u003e \\u003ctd align=\\\"left\\\" colname=\\\"c5\\\"\\u003e\\u0026nbsp;\\u003c/td\\u003e \\u003c/tr\\u003e \\u003ctr\\u003e \\u003ctd align=\\\"left\\\" colname=\\\"c1\\\"\\u003e \\u003cp\\u003e\\u003cem\\u003ePhacidium\\u003c/em\\u003e\\u003c/p\\u003e \\u003c/td\\u003e \\u003ctd align=\\\"left\\\" colname=\\\"c2\\\"\\u003e \\u003cp\\u003e3\\u003c/p\\u003e \\u003c/td\\u003e \\u003ctd align=\\\"left\\\" colname=\\\"c3\\\"\\u003e \\u003cp\\u003e1\\u003c/p\\u003e \\u003c/td\\u003e \\u003ctd align=\\\"left\\\" colname=\\\"c4\\\"\\u003e \\u003cp\\u003ex\\u003c/p\\u003e \\u003c/td\\u003e \\u003ctd align=\\\"left\\\" colname=\\\"c5\\\"\\u003e\\u0026nbsp;\\u003c/td\\u003e \\u003c/tr\\u003e \\u003ctr\\u003e \\u003ctd align=\\\"left\\\" colname=\\\"c1\\\"\\u003e \\u003cp\\u003eunidentified\\u003c/p\\u003e \\u003c/td\\u003e \\u003ctd align=\\\"left\\\" colname=\\\"c2\\\"\\u003e \\u003cp\\u003e324\\u003c/p\\u003e \\u003c/td\\u003e \\u003ctd align=\\\"left\\\" colname=\\\"c3\\\"\\u003e \\u003cp\\u003e54\\u003c/p\\u003e \\u003c/td\\u003e \\u003ctd align=\\\"left\\\" colname=\\\"c4\\\"\\u003e\\u0026nbsp;\\u003c/td\\u003e \\u003ctd align=\\\"left\\\" colname=\\\"c5\\\"\\u003e\\u0026nbsp;\\u003c/td\\u003e \\u003c/tr\\u003e \\u003ctr\\u003e \\u003ctd align=\\\"left\\\" colspan=\\\"5\\\" nameend=\\\"c5\\\" namest=\\\"c1\\\"\\u003e \\u003cp\\u003e\\u003cb\\u003eAbies alba\\u003c/b\\u003e\\u003c/p\\u003e \\u003c/td\\u003e \\u003c/tr\\u003e \\u003ctr\\u003e \\u003ctd align=\\\"left\\\" colname=\\\"c1\\\"\\u003e \\u003cp\\u003e\\u003cem\\u003ePenicillium\\u003c/em\\u003e\\u003c/p\\u003e \\u003c/td\\u003e \\u003ctd align=\\\"left\\\" colname=\\\"c2\\\"\\u003e \\u003cp\\u003e28\\u003c/p\\u003e \\u003c/td\\u003e \\u003ctd align=\\\"left\\\" colname=\\\"c3\\\"\\u003e \\u003cp\\u003e7\\u003c/p\\u003e \\u003c/td\\u003e \\u003ctd align=\\\"left\\\" colname=\\\"c4\\\"\\u003e \\u003cp\\u003ex\\u003c/p\\u003e \\u003c/td\\u003e \\u003ctd align=\\\"left\\\" colname=\\\"c5\\\"\\u003e \\u003cp\\u003ex\\u003c/p\\u003e \\u003c/td\\u003e \\u003c/tr\\u003e \\u003ctr\\u003e \\u003ctd align=\\\"left\\\" colname=\\\"c1\\\"\\u003e \\u003cp\\u003e\\u003cem\\u003eTrichoderma\\u003c/em\\u003e\\u003c/p\\u003e \\u003c/td\\u003e \\u003ctd align=\\\"left\\\" colname=\\\"c2\\\"\\u003e \\u003cp\\u003e27\\u003c/p\\u003e \\u003c/td\\u003e \\u003ctd align=\\\"left\\\" colname=\\\"c3\\\"\\u003e \\u003cp\\u003e2\\u003c/p\\u003e \\u003c/td\\u003e \\u003ctd align=\\\"left\\\" colname=\\\"c4\\\"\\u003e \\u003cp\\u003ex\\u003c/p\\u003e \\u003c/td\\u003e \\u003ctd align=\\\"left\\\" colname=\\\"c5\\\"\\u003e \\u003cp\\u003ex\\u003c/p\\u003e \\u003c/td\\u003e \\u003c/tr\\u003e \\u003ctr\\u003e \\u003ctd align=\\\"left\\\" colname=\\\"c1\\\"\\u003e \\u003cp\\u003e\\u003cem\\u003eSydowia\\u003c/em\\u003e\\u003c/p\\u003e \\u003c/td\\u003e \\u003ctd align=\\\"left\\\" colname=\\\"c2\\\"\\u003e \\u003cp\\u003e25\\u003c/p\\u003e \\u003c/td\\u003e \\u003ctd align=\\\"left\\\" colname=\\\"c3\\\"\\u003e \\u003cp\\u003e2\\u003c/p\\u003e \\u003c/td\\u003e \\u003ctd align=\\\"left\\\" colname=\\\"c4\\\"\\u003e\\u0026nbsp;\\u003c/td\\u003e \\u003ctd align=\\\"left\\\" colname=\\\"c5\\\"\\u003e\\u0026nbsp;\\u003c/td\\u003e \\u003c/tr\\u003e \\u003ctr\\u003e \\u003ctd align=\\\"left\\\" colname=\\\"c1\\\"\\u003e \\u003cp\\u003e\\u003cem\\u003eFusarium\\u003c/em\\u003e\\u003c/p\\u003e \\u003c/td\\u003e \\u003ctd align=\\\"left\\\" colname=\\\"c2\\\"\\u003e \\u003cp\\u003e7\\u003c/p\\u003e \\u003c/td\\u003e \\u003ctd align=\\\"left\\\" colname=\\\"c3\\\"\\u003e \\u003cp\\u003e2\\u003c/p\\u003e \\u003c/td\\u003e \\u003ctd align=\\\"left\\\" colname=\\\"c4\\\"\\u003e \\u003cp\\u003ex\\u003c/p\\u003e \\u003c/td\\u003e \\u003ctd align=\\\"left\\\" colname=\\\"c5\\\"\\u003e \\u003cp\\u003ex\\u003c/p\\u003e \\u003c/td\\u003e \\u003c/tr\\u003e \\u003ctr\\u003e \\u003ctd align=\\\"left\\\" colname=\\\"c1\\\"\\u003e \\u003cp\\u003e\\u003cem\\u003eAspergillus\\u003c/em\\u003e\\u003c/p\\u003e \\u003c/td\\u003e \\u003ctd align=\\\"left\\\" colname=\\\"c2\\\"\\u003e \\u003cp\\u003e6\\u003c/p\\u003e \\u003c/td\\u003e \\u003ctd align=\\\"left\\\" colname=\\\"c3\\\"\\u003e \\u003cp\\u003e3\\u003c/p\\u003e \\u003c/td\\u003e \\u003ctd align=\\\"left\\\" colname=\\\"c4\\\"\\u003e \\u003cp\\u003ex\\u003c/p\\u003e \\u003c/td\\u003e \\u003ctd align=\\\"left\\\" colname=\\\"c5\\\"\\u003e \\u003cp\\u003ex\\u003c/p\\u003e \\u003c/td\\u003e \\u003c/tr\\u003e \\u003ctr\\u003e \\u003ctd align=\\\"left\\\" colname=\\\"c1\\\"\\u003e \\u003cp\\u003e\\u003cem\\u003eMucor\\u003c/em\\u003e\\u003c/p\\u003e \\u003c/td\\u003e \\u003ctd align=\\\"left\\\" colname=\\\"c2\\\"\\u003e \\u003cp\\u003e2\\u003c/p\\u003e \\u003c/td\\u003e \\u003ctd align=\\\"left\\\" colname=\\\"c3\\\"\\u003e \\u003cp\\u003e1\\u003c/p\\u003e \\u003c/td\\u003e \\u003ctd align=\\\"left\\\" colname=\\\"c4\\\"\\u003e\\u0026nbsp;\\u003c/td\\u003e \\u003ctd align=\\\"left\\\" colname=\\\"c5\\\"\\u003e\\u0026nbsp;\\u003c/td\\u003e \\u003c/tr\\u003e \\u003ctr\\u003e \\u003ctd align=\\\"left\\\" colname=\\\"c1\\\"\\u003e \\u003cp\\u003e\\u003cem\\u003eTalaromyces\\u003c/em\\u003e\\u003c/p\\u003e \\u003c/td\\u003e \\u003ctd align=\\\"left\\\" colname=\\\"c2\\\"\\u003e \\u003cp\\u003e1\\u003c/p\\u003e \\u003c/td\\u003e \\u003ctd align=\\\"left\\\" colname=\\\"c3\\\"\\u003e \\u003cp\\u003e1\\u003c/p\\u003e \\u003c/td\\u003e \\u003ctd align=\\\"left\\\" colname=\\\"c4\\\"\\u003e \\u003cp\\u003ex\\u003c/p\\u003e \\u003c/td\\u003e \\u003ctd align=\\\"left\\\" colname=\\\"c5\\\"\\u003e \\u003cp\\u003ex\\u003c/p\\u003e \\u003c/td\\u003e \\u003c/tr\\u003e \\u003ctr\\u003e \\u003ctd align=\\\"left\\\" colname=\\\"c1\\\"\\u003e \\u003cp\\u003e\\u003cem\\u003eAkanthomyces\\u003c/em\\u003e\\u003c/p\\u003e \\u003c/td\\u003e \\u003ctd align=\\\"left\\\" colname=\\\"c2\\\"\\u003e \\u003cp\\u003e1\\u003c/p\\u003e \\u003c/td\\u003e \\u003ctd align=\\\"left\\\" colname=\\\"c3\\\"\\u003e \\u003cp\\u003e1\\u003c/p\\u003e \\u003c/td\\u003e \\u003ctd align=\\\"left\\\" colname=\\\"c4\\\"\\u003e\\u0026nbsp;\\u003c/td\\u003e \\u003ctd align=\\\"left\\\" colname=\\\"c5\\\"\\u003e\\u0026nbsp;\\u003c/td\\u003e \\u003c/tr\\u003e \\u003ctr\\u003e \\u003ctd align=\\\"left\\\" colname=\\\"c1\\\"\\u003e \\u003cp\\u003eunidentified\\u003c/p\\u003e \\u003c/td\\u003e \\u003ctd align=\\\"left\\\" colname=\\\"c2\\\"\\u003e \\u003cp\\u003e134\\u003c/p\\u003e \\u003c/td\\u003e \\u003ctd align=\\\"left\\\" colname=\\\"c3\\\"\\u003e \\u003cp\\u003e25\\u003c/p\\u003e \\u003c/td\\u003e \\u003ctd align=\\\"left\\\" colname=\\\"c4\\\"\\u003e\\u0026nbsp;\\u003c/td\\u003e \\u003ctd align=\\\"left\\\" colname=\\\"c5\\\"\\u003e\\u0026nbsp;\\u003c/td\\u003e \\u003c/tr\\u003e \\u003ctr\\u003e \\u003ctd align=\\\"left\\\" colspan=\\\"5\\\" nameend=\\\"c5\\\" namest=\\\"c1\\\"\\u003e \\u003cp\\u003e\\u003cb\\u003ePinus sylvestris\\u003c/b\\u003e\\u003c/p\\u003e \\u003c/td\\u003e \\u003c/tr\\u003e \\u003ctr\\u003e \\u003ctd align=\\\"left\\\" colname=\\\"c1\\\"\\u003e \\u003cp\\u003e\\u003cem\\u003ePenicillium\\u003c/em\\u003e\\u003c/p\\u003e \\u003c/td\\u003e \\u003ctd align=\\\"left\\\" colname=\\\"c2\\\"\\u003e \\u003cp\\u003e5\\u003c/p\\u003e \\u003c/td\\u003e \\u003ctd align=\\\"left\\\" colname=\\\"c3\\\"\\u003e \\u003cp\\u003e5\\u003c/p\\u003e \\u003c/td\\u003e \\u003ctd align=\\\"left\\\" colname=\\\"c4\\\"\\u003e \\u003cp\\u003ex\\u003c/p\\u003e \\u003c/td\\u003e \\u003ctd align=\\\"left\\\" colname=\\\"c5\\\"\\u003e \\u003cp\\u003ex\\u003c/p\\u003e \\u003c/td\\u003e \\u003c/tr\\u003e \\u003ctr\\u003e \\u003ctd align=\\\"left\\\" colname=\\\"c1\\\"\\u003e \\u003cp\\u003e\\u003cem\\u003eFusarium\\u003c/em\\u003e\\u003c/p\\u003e \\u003c/td\\u003e \\u003ctd align=\\\"left\\\" colname=\\\"c2\\\"\\u003e \\u003cp\\u003e5\\u003c/p\\u003e \\u003c/td\\u003e \\u003ctd align=\\\"left\\\" colname=\\\"c3\\\"\\u003e \\u003cp\\u003e1\\u003c/p\\u003e \\u003c/td\\u003e \\u003ctd align=\\\"left\\\" colname=\\\"c4\\\"\\u003e \\u003cp\\u003ex\\u003c/p\\u003e \\u003c/td\\u003e \\u003ctd align=\\\"left\\\" colname=\\\"c5\\\"\\u003e \\u003cp\\u003ex\\u003c/p\\u003e \\u003c/td\\u003e \\u003c/tr\\u003e \\u003ctr\\u003e \\u003ctd align=\\\"left\\\" colname=\\\"c1\\\"\\u003e \\u003cp\\u003e\\u003cem\\u003eBoeremia\\u003c/em\\u003e\\u003c/p\\u003e \\u003c/td\\u003e \\u003ctd align=\\\"left\\\" colname=\\\"c2\\\"\\u003e \\u003cp\\u003e4\\u003c/p\\u003e \\u003c/td\\u003e \\u003ctd align=\\\"left\\\" colname=\\\"c3\\\"\\u003e \\u003cp\\u003e2\\u003c/p\\u003e \\u003c/td\\u003e \\u003ctd align=\\\"left\\\" colname=\\\"c4\\\"\\u003e\\u0026nbsp;\\u003c/td\\u003e \\u003ctd align=\\\"left\\\" colname=\\\"c5\\\"\\u003e\\u0026nbsp;\\u003c/td\\u003e \\u003c/tr\\u003e \\u003ctr\\u003e \\u003ctd align=\\\"left\\\" colname=\\\"c1\\\"\\u003e \\u003cp\\u003e\\u003cem\\u003eConiochaeta\\u003c/em\\u003e\\u003c/p\\u003e \\u003c/td\\u003e \\u003ctd align=\\\"left\\\" colname=\\\"c2\\\"\\u003e \\u003cp\\u003e3\\u003c/p\\u003e \\u003c/td\\u003e \\u003ctd align=\\\"left\\\" colname=\\\"c3\\\"\\u003e \\u003cp\\u003e1\\u003c/p\\u003e \\u003c/td\\u003e \\u003ctd align=\\\"left\\\" colname=\\\"c4\\\"\\u003e\\u0026nbsp;\\u003c/td\\u003e \\u003ctd align=\\\"left\\\" colname=\\\"c5\\\"\\u003e \\u003cp\\u003ex\\u003c/p\\u003e \\u003c/td\\u003e \\u003c/tr\\u003e \\u003ctr\\u003e \\u003ctd align=\\\"left\\\" colname=\\\"c1\\\"\\u003e \\u003cp\\u003e\\u003cem\\u003eAspergillus\\u003c/em\\u003e\\u003c/p\\u003e \\u003c/td\\u003e \\u003ctd align=\\\"left\\\" colname=\\\"c2\\\"\\u003e \\u003cp\\u003e1\\u003c/p\\u003e \\u003c/td\\u003e \\u003ctd align=\\\"left\\\" colname=\\\"c3\\\"\\u003e \\u003cp\\u003e1\\u003c/p\\u003e \\u003c/td\\u003e \\u003ctd align=\\\"left\\\" colname=\\\"c4\\\"\\u003e \\u003cp\\u003ex\\u003c/p\\u003e \\u003c/td\\u003e \\u003ctd align=\\\"left\\\" colname=\\\"c5\\\"\\u003e \\u003cp\\u003ex\\u003c/p\\u003e \\u003c/td\\u003e \\u003c/tr\\u003e \\u003ctr\\u003e \\u003ctd align=\\\"left\\\" colname=\\\"c1\\\"\\u003e \\u003cp\\u003e\\u003cem\\u003eSydowia\\u003c/em\\u003e\\u003c/p\\u003e \\u003c/td\\u003e \\u003ctd align=\\\"left\\\" colname=\\\"c2\\\"\\u003e \\u003cp\\u003e1\\u003c/p\\u003e \\u003c/td\\u003e \\u003ctd align=\\\"left\\\" colname=\\\"c3\\\"\\u003e \\u003cp\\u003e1\\u003c/p\\u003e \\u003c/td\\u003e \\u003ctd align=\\\"left\\\" colname=\\\"c4\\\"\\u003e\\u0026nbsp;\\u003c/td\\u003e \\u003ctd align=\\\"left\\\" colname=\\\"c5\\\"\\u003e\\u0026nbsp;\\u003c/td\\u003e \\u003c/tr\\u003e \\u003ctr\\u003e \\u003ctd align=\\\"left\\\" colname=\\\"c1\\\"\\u003e \\u003cp\\u003eunidentified\\u003c/p\\u003e \\u003c/td\\u003e \\u003ctd align=\\\"left\\\" colname=\\\"c2\\\"\\u003e \\u003cp\\u003e3\\u003c/p\\u003e \\u003c/td\\u003e \\u003ctd align=\\\"left\\\" colname=\\\"c3\\\"\\u003e \\u003cp\\u003e3\\u003c/p\\u003e \\u003c/td\\u003e \\u003ctd align=\\\"left\\\" colname=\\\"c4\\\"\\u003e\\u0026nbsp;\\u003c/td\\u003e \\u003ctd align=\\\"left\\\" colname=\\\"c5\\\"\\u003e\\u0026nbsp;\\u003c/td\\u003e \\u003c/tr\\u003e \\u003c/tbody\\u003e \\u003c/colgroup\\u003e \\u003c/table\\u003e\\u003c/div\\u003e \\u003c/p\\u003e\"},{\"header\":\"Discussion\",\"content\":\"\\u003cp\\u003eTraditional culturing methods reveal limited fungal diversity, are resource-intensive, and have largely been replaced by the high-throughput metabarcoding approaches in studies where living cultures are not required. Metabarcoding studies are predominantly DNA-based, which does not allow differentiation between metabolically active and inactive taxa, including dormant or dead organisms. In contrast, RNA metabarcoding targets living, metabolically active organisms, potentially providing a more ecologically relevant snapshot of fungal communities compared to DNA metabarcoding. RNA metabarcoding has shown potential for studying active environmental microbiota [\\u003cspan citationid=\\\"CR18\\\" class=\\\"CitationRef\\\"\\u003e18\\u003c/span\\u003e, \\u003cspan citationid=\\\"CR22\\\" class=\\\"CitationRef\\\"\\u003e22\\u003c/span\\u003e, \\u003cspan citationid=\\\"CR25\\\" class=\\\"CitationRef\\\"\\u003e25\\u003c/span\\u003e, \\u003cspan citationid=\\\"CR56\\\" class=\\\"CitationRef\\\"\\u003e56\\u003c/span\\u003e]. However, the method has rarely been applied in studies of plant-associated microbiota. Given its high potential for ecological and phytosanitary monitoring, it is crucial to evaluate the RNA metabarcoding method against other commonly applied methods for studying plant-associated fungi before it can be widely implemented.\\u003c/p\\u003e \\u003cp\\u003eIn this study, we compared incidence-based and abundance-weighted seed-borne fungal diversity and community composition between DNA and RNA metabarcoding methods and tree species. Additionally, we assessed the overlap between fungal genera identified through culturing and those revealed by each metabarcoding method. Our results show that the two metabarcoding approaches capture distinct and largely non-overlapping fungal communities, with the observed differences primarily driven by rare OTUs in the RNA dataset. However, both approaches recovered similar cultured genera, which likely represent abundant and metabolically active taxa.\\u003c/p\\u003e \\u003cp\\u003eFungal diversity and community composition across metabarcoding and culturing datasets\\u003c/p\\u003e \\u003cp\\u003eContrary to our hypothesis, RNA metabarcoding revealed a higher number of fungal OTUs than DNA metabarcoding, including many unique OTUs, resulting in highly divergent community profiles. This unexpected outcome likely results from the high abundance of taxa represented by dead or dormant forms [\\u003cspan citationid=\\\"CR57\\\" class=\\\"CitationRef\\\"\\u003e57\\u003c/span\\u003e], as well as the large diversity of metabolically active rare or low-abundance taxa in dormant tree seeds. In complex or competitive environments\\u0026mdash;such as dormant seeds\\u0026mdash;metabolic activity often originates from low abundance taxa [\\u003cspan citationid=\\\"CR58\\\" class=\\\"CitationRef\\\"\\u003e58\\u003c/span\\u003e], as demonstrated by rare bacterial taxa showing high metabolic activity in atmospheric [\\u003cspan citationid=\\\"CR56\\\" class=\\\"CitationRef\\\"\\u003e56\\u003c/span\\u003e, \\u003cspan citationid=\\\"CR59\\\" class=\\\"CitationRef\\\"\\u003e59\\u003c/span\\u003e] and oceanic environments [\\u003cspan citationid=\\\"CR60\\\" class=\\\"CitationRef\\\"\\u003e60\\u003c/span\\u003e]. While DNA metabarcoding typically captures highly abundant taxa, RNA metabarcoding targets metabolically active fungi, leading to discrepancies in community profiles, especially when incidence-based measures are used. However, although abundance-weighted alpha diversity was also lower in the DNA than in the RNA dataset, the divergence between DNA- and RNA-based community profiles decreased when abundance-weighted metrics were applied. Among the most abundant genera identified per tree species and method, about half were shared between the DNA and RNA datasets in two out of three tree species. This indicates that many dominant fungi are generally metabolically active, and that discrepancies between methods largely reflect differences among rare active taxa.\\u003c/p\\u003e \\u003cp\\u003eMethodological factors likely contributed to the observed differences in fungal communities between the DNA and RNA datasets. A first bias may have been introduced during nucleic acid extraction, as DNA and RNA were obtained from subsamples of one seed lot and using different extraction kits. Previous studies have demonstrated that the choice of extraction kit can significantly influence microbial community composition [\\u003cspan citationid=\\\"CR61\\\" class=\\\"CitationRef\\\"\\u003e61\\u003c/span\\u003e], a bias that could be mitigated by using co-extraction kits that allow simultaneous isolation of both DNA and RNA [\\u003cspan citationid=\\\"CR62\\\" class=\\\"CitationRef\\\"\\u003e62\\u003c/span\\u003e]. In addition, subsampling procedures differed between the extractions. While a single DNA extraction was performed per sample, RNA was extracted in duplicates and pooled to obtain sufficient input for reverse transcription. This variation in the amount of starting material may have led to higher alpha diversity in the RNA than DNA dataset and shifted community composition\\u0026mdash;an effect previously shown for seed mycobiota of \\u003cem\\u003eP. sylvestris\\u003c/em\\u003e [\\u003cspan citationid=\\\"CR63\\\" class=\\\"CitationRef\\\"\\u003e63\\u003c/span\\u003e]. Nonetheless, our study accounted for differences in sampling effort through rarefaction and the use of abundance-weighted diversity metrics [\\u003cspan citationid=\\\"CR64\\\" class=\\\"CitationRef\\\"\\u003e64\\u003c/span\\u003e, \\u003cspan citationid=\\\"CR65\\\" class=\\\"CitationRef\\\"\\u003e65\\u003c/span\\u003e], helping to mitigate potential confounding effects.\\u003c/p\\u003e \\u003cp\\u003eFurthermore, the reverse transcription step required to synthesize cDNA from RNA introduces biases such as uneven abundances of correctly transcribed molecules. The reverse transcription can also introduce artefacts due to incorrect primer binding or unpredictable reverse transcriptase behaviour [\\u003cspan citationid=\\\"CR66\\\" class=\\\"CitationRef\\\"\\u003e66\\u003c/span\\u003e]. These factors likely contributed to discrepancies between fungal community profiles derived from DNA and RNA. Additional bias may occur during PCR amplification of the ITS region originating from both DNA and cDNA, and library preparation, including primer-template mismatches and selective amplification of certain taxa [\\u003cspan citationid=\\\"CR67\\\" class=\\\"CitationRef\\\"\\u003e67\\u003c/span\\u003e]. Finally, DNA and RNA libraries were sequenced on separate runs, which may have introduced further variation, although, this effect is likely less pronounced than the other previously mentioned methodological sources of bias [\\u003cspan citationid=\\\"CR61\\\" class=\\\"CitationRef\\\"\\u003e61\\u003c/span\\u003e].\\u003c/p\\u003e \\u003cp\\u003eWe found no evidence supporting our hypothesis that cultured fungi would more closely resemble the RNA-derived community than the DNA-derived one. Approximately 50% of cultured genera were detected by both metabarcoding approaches, while an additional 6% was uniquely recovered by either DNA- or RNA-based metabarcoding. These results suggest that both metabolic activity and abundance influence the success of isolation of the fungus. If culturing primarily captured metabolically active taxa, a stronger overlap with the RNA dataset\\u0026mdash;representing rare but active fungi\\u0026mdash;would have been expected. Instead, our culturing approach likely favoured abundant and competitive fungi [\\u003cspan citationid=\\\"CR8\\\" class=\\\"CitationRef\\\"\\u003e8\\u003c/span\\u003e]. Alternative strategies such as dilution-to-extinction culturing [\\u003cspan citationid=\\\"CR68\\\" class=\\\"CitationRef\\\"\\u003e68\\u003c/span\\u003e] or use of different media [\\u003cspan citationid=\\\"CR69\\\" class=\\\"CitationRef\\\"\\u003e69\\u003c/span\\u003e] could help recover a greater diversity of rare yet cultivable taxa and potentially increase overlap with RNA-based profiles. As expected, culturing recovered substantially lower diversity than metabarcoding, reflecting the well-known difficulty of growing many fungal taxa under standard laboratory conditions [\\u003cspan citationid=\\\"CR8\\\" class=\\\"CitationRef\\\"\\u003e8\\u003c/span\\u003e]. Moreover, taxa detected only by culturing may reflect stochastic recovery or primer bias between methods. As isolates are critical for species description, functional studies, and phenotypic assays, continued improvement of culture-dependent techniques remains essential. Developing more sensitive and high-throughput culturing approaches will be key to bridging molecular and culture-based views on fungal diversity.\\u003c/p\\u003e \\u003cp\\u003eTaxonomic and functional fungal community composition across metabarcoding datasets and tree species\\u003c/p\\u003e \\u003cp\\u003eTree species identity was a major driver of the community composition of dominant fungi, explaining approximately 30% of the overall observed variation. The high host specificity of seed-borne fungi observed in our study aligns with findings from both woody and non-woody plants [\\u003cspan citationid=\\\"CR70\\\" class=\\\"CitationRef\\\"\\u003e70\\u003c/span\\u003e, \\u003cspan citationid=\\\"CR71\\\" class=\\\"CitationRef\\\"\\u003e71\\u003c/span\\u003e]. Although, several dominant fungal genera were associated with already known hosts, not all of them were previously known to occur in seeds. For example, the genus \\u003cem\\u003eApiognomonia\\u003c/em\\u003e was the most abundant in \\u003cem\\u003eF. sylvatica\\u003c/em\\u003e seed, likely corresponding to \\u003cem\\u003eA. errabunda\\u003c/em\\u003e, a common endophyte of \\u003cem\\u003eF. sylvatica\\u003c/em\\u003e known to cause anthracnose under wet spring conditions or following insect damage [\\u003cspan citationid=\\\"CR72\\\" class=\\\"CitationRef\\\"\\u003e72\\u003c/span\\u003e], and was also identified in the culturing dataset. \\u003cem\\u003eAbies alba\\u003c/em\\u003e seeds were dominated by OTUs assigned to the genus \\u003cem\\u003eLambertella\\u003c/em\\u003e which has been previously found in \\u003cem\\u003eA. koreana\\u003c/em\\u003e needles [\\u003cspan citationid=\\\"CR73\\\" class=\\\"CitationRef\\\"\\u003e73\\u003c/span\\u003e]. \\u003cem\\u003ePinus sylvestris\\u003c/em\\u003e seeds were dominated by the genus \\u003cem\\u003eHormonema\\u003c/em\\u003e, likely corresponding to \\u003cem\\u003eSydowia polyspora\\u003c/em\\u003e (i.e., old name is \\u003cem\\u003eH. dematioides\\u003c/em\\u003e), a well-known endophyte and opportunistic pathogen in conifers and pre-emergent seed pathogen in \\u003cem\\u003eP. ponderosa\\u003c/em\\u003e [\\u003cspan citationid=\\\"CR6\\\" class=\\\"CitationRef\\\"\\u003e6\\u003c/span\\u003e, \\u003cspan citationid=\\\"CR74\\\" class=\\\"CitationRef\\\"\\u003e74\\u003c/span\\u003e], which was also identified in the culturing dataset.\\u003c/p\\u003e \\u003cp\\u003eSome fungal genera observed in seeds of specific tree species were found to be strongly associated with one of the metabarcoding datasets. For example, the genus \\u003cem\\u003eParaleptosphaeria\\u003c/em\\u003e was strongly associated with the DNA dataset of \\u003cem\\u003eF. sylvatica\\u003c/em\\u003e. This genus is known to be saprobic, fungicolous, or pathogenic, and has been found in grasslands and on stems and leaves of various plants. However, \\u003cem\\u003eFagus\\u003c/em\\u003e has not been listed among confirmed hosts in recent studies [\\u003cspan citationid=\\\"CR75\\\" class=\\\"CitationRef\\\"\\u003e75\\u003c/span\\u003e, \\u003cspan citationid=\\\"CR76\\\" class=\\\"CitationRef\\\"\\u003e76\\u003c/span\\u003e]. The genus \\u003cem\\u003eCaloscypha\\u003c/em\\u003e, likely represented by the species \\u003cem\\u003eC. fulgens\\u003c/em\\u003e\\u0026mdash;a well-known conifer seed pathogen [\\u003cspan citationid=\\\"CR77\\\" class=\\\"CitationRef\\\"\\u003e77\\u003c/span\\u003e] was strongly associated with the \\u003cem\\u003eA. alba\\u003c/em\\u003e and \\u003cem\\u003eP. sylvestris\\u003c/em\\u003e DNA dataset. The genus \\u003cem\\u003eLophodermium\\u003c/em\\u003e showed strong association with the RNA data of \\u003cem\\u003eF. sylvatica\\u003c/em\\u003e and \\u003cem\\u003eA. alba\\u003c/em\\u003e and was represented by multiple OTUs identified to several species known as needle endophytes or pathogens of conifers [\\u003cspan citationid=\\\"CR78\\\" class=\\\"CitationRef\\\"\\u003e78\\u003c/span\\u003e]. Association of \\u003cem\\u003eLophodermium\\u003c/em\\u003e with \\u003cem\\u003eF. sylvatica\\u003c/em\\u003e may stem from environmental or laboratory contamination. The genus \\u003cem\\u003eNeocatenulostroma\\u003c/em\\u003e was strongly associated with the RNA dataset of \\u003cem\\u003eA. alba\\u003c/em\\u003e and could correspond to \\u003cem\\u003eN. abietis\\u003c/em\\u003e, species which was isolated from various substrates and was found as a saprobe or endophyte in pine needles, but also as a pathogen on a wide range of conifer hosts [\\u003cspan citationid=\\\"CR79\\\" class=\\\"CitationRef\\\"\\u003e79\\u003c/span\\u003e]. In \\u003cem\\u003eP. sylvestris\\u003c/em\\u003e, \\u003cem\\u003eColletotrichum\\u003c/em\\u003e was strongly associated with the RNA dataset. Some \\u003cem\\u003eColletotrichum\\u003c/em\\u003e species are known pathogens of pines, although their occurrence in conifers is less frequent than in broadleaf hosts [\\u003cspan citationid=\\\"CR80\\\" class=\\\"CitationRef\\\"\\u003e80\\u003c/span\\u003e]. These results highlight the value of combining DNA and RNA-based metabarcoding to uncover both inactive and active fungal associations, some of which may have implication for tree health.\\u003c/p\\u003e \\u003cp\\u003eSimilar to findings by Mittelstrass et al. [\\u003cspan citationid=\\\"CR3\\\" class=\\\"CitationRef\\\"\\u003e3\\u003c/span\\u003e] most fungal reads in our study were assigned to saprotrophs and plant pathogens, which appeared in similar proportions across DNA and RNA metabarcoding datasets of all three host species. This suggests that these functional groups are not only abundant but also metabolically active in dormant tree seeds. Although dormant seeds may not seem like suitable habitats\\u0026mdash;given that saprotrophs rely on dead organic matter and plant pathogens require a susceptible host\\u0026mdash;research suggests that many potentially saprotrophic or pathogenic taxa can persist in asymptomatic plant tissues as endophytes (i.e., asymptomatic, commensal or weakly mutualistic inhabitants; [\\u003cspan citationid=\\\"CR43\\\" class=\\\"CitationRef\\\"\\u003e43\\u003c/span\\u003e]) and shift to saprotrophic or pathogenic lifestyle under favourable conditions [\\u003cspan citationid=\\\"CR81\\\" class=\\\"CitationRef\\\"\\u003e81\\u003c/span\\u003e]. Mycorrhizal fungi and endophytes were detected almost exclusively in RNA datasets, likely due to their low abundance, which may have prevented their detection by DNA metabarcoding. The sporadic occurrence across seed lots and low abundance of mycorrhizal fungi in seeds supports the notion that their seed-transmission is uncommon in trees as previously suggested [\\u003cspan citationid=\\\"CR82\\\" class=\\\"CitationRef\\\"\\u003e82\\u003c/span\\u003e]. However, the detection of taxa belonging to mycorrhizal fungi in seeds of various non-woody plant species [\\u003cspan citationid=\\\"CR70\\\" class=\\\"CitationRef\\\"\\u003e70\\u003c/span\\u003e] highlights the need to further investigate mycorrhizal transmission from seeds to seedlings. Low numbers of endophytic taxa are likely associated with their underrepresentation in databases such as the FungalTraits [\\u003cspan citationid=\\\"CR43\\\" class=\\\"CitationRef\\\"\\u003e43\\u003c/span\\u003e], underscoring the need for more targeted research into their diversity and function.\\u003c/p\\u003e\"},{\"header\":\"Conclusions\",\"content\":\"\\u003cp\\u003eOur results show that DNA and RNA metabarcoding provide distinct but complementary insights into fungal communities and have valuable potential for high-throughput ecological and biosecurity monitoring. DNA metabarcoding has already been widely adopted for assessing biodiversity of different organisms in both environmental and host-associated samples. Its ability to detect living and non-living taxa provides a comprehensive perspective on total fungal community associated with the studied system. In contrast, RNA-based metabarcoding is still less commonly applied, despite its focus on metabolically active organisms which may enable real-time insights into ecosystem dynamics and functional responses to environmental change. Moreover, RNA metabarcoding might be particularly promising for biosecurity application where rapid detection and identification of living pests and pathogens is critical for decision making and timely intervention.\\u003c/p\\u003e \\u003cp\\u003eIn our study, RNA metabarcoding revealed a higher diversity, primarily of rare taxa, compared to DNA metabarcoding, making it particularly effective for detecting low-abundance pathogens, at least when dominant taxa are inactive or dormant (e.g., ensuring absence of pathogens in treated samples). However, its inability to detect metabolically inactive yet viable organisms represents a limitation as dormant pathogens may escape detection and later become active. This could be tackled with applying viability PCR (vPCR) using propidium or ethidium monoazide which enables selective suppression of DNA from membrane-compromised cells, allowing DNA-based assays to target only intact, viable fungi and thus reduce false positives from relic DNA [\\u003cspan citationid=\\\"CR83\\\" class=\\\"CitationRef\\\"\\u003e83\\u003c/span\\u003e, \\u003cspan citationid=\\\"CR84\\\" class=\\\"CitationRef\\\"\\u003e84\\u003c/span\\u003e]. Nevertheless, vPCR also has important limitations. Its effectiveness can vary depending on cell wall structure and environmental matrix composition, and incomplete dye penetration or binding may lead to false negatives. Moreover, its applicability to mixed environmental samples such as those that need to be ground prior to DNA/RNA extraction, remains uncertain. If technical obstacles are overcome, performing DNA metabarcoding with and without vPCR pretreatment could help differentiate living from dead organisms. Combining this with RNA-based methods can further help identify the metabolically active members of the community, leading to more accurate ecological interpretations and risk assessments.\\u003c/p\\u003e \\u003cp\\u003eLooking ahead, advances in sequencing technologies and bioinformatics are expected to increase the sensitivity and resolution of metabarcoding methods, enabling more precise taxonomic identification and functional profiling. In the face of growing environmental pressures and accelerating species movement through trade, the combined use of DNA and RNA metabarcoding holds strong potential for ecosystem management and biosecurity, by revealing total and active communities. While culturing remains essential for studying fungal diversity, pathogenicity and functional traits, traditional approaches recover only a subset of the taxa detected by metabarcoding, typically the dominant and readily cultivable species. Future efforts should thus focus on developing more sensitive, high-throughput culturing techniques capable of capturing a broader spectrum of fungal diversity. Ultimately, integrating all three approaches\\u0026mdash;DNA- and RNA-based metabarcoding together with advanced culturing\\u0026mdash;offers the most comprehensive framework for studying fungal diversity and strengthening biosecurity, by capturing both active and inactive taxa as well as culturable and unculturable species.\\u003c/p\\u003e\"},{\"header\":\"Declarations\",\"content\":\"\\u003ch2\\u003eEthics approval and consent to participate\\u003c/h2\\u003e\\n\\u003cp\\u003eNot applicable.\\u003c/p\\u003e\\n\\u003ch2\\u003eConsent for publication\\u003c/h2\\u003e\\n\\u003cp\\u003eNot applicable.\\u003c/p\\u003e\\n\\u003ch2\\u003eData Availability\\u003c/h2\\u003e\\n\\u003cp\\u003eRaw metabarcoding sequence reads have been deposited in the European Nucleotide Archive (ENA) at EMBL-EBI under accession number PRJEB96213.\\u003c/p\\u003e\\n\\u003cp\\u003eITS barcode sequences generated from cultured specimens have been deposited in GenBank under accession numbers PX736467\\u0026ndash;PX736489, PX740711\\u0026ndash;PX740835, and PX745229\\u0026ndash;PX745242. All ITS sequences used in this study, together with their corresponding accession numbers, are provided in Supplementary Table\\u0026nbsp;4.\\u003c/p\\u003e\\n\\u003cp\\u003eCompeting Interests\\u003c/p\\u003e\\n\\u003cp\\u003eThe authors declare that they have no competing interests.\\u003c/p\\u003e\\n\\u003ch2\\u003eFunding\\u003c/h2\\u003e\\n\\u003cp\\u003eThis work was funded by the Swiss National Science Foundation SNSF (Grant: P500PB_211040) and by the Swiss Federal Office for the Environment (FOEN) (Finanzhilfevertrag betreffend wissenschaftlich-technischen T\\u0026auml;tigkeiten im Bereich Waldschutz).\\u003c/p\\u003e\\n\\u003cp\\u003eRen\\u0026eacute; Eschen was supported by CABI. CABI is an international intergovernmental organisation, and we gratefully acknowledge the core financial support from our member countries (and lead agencies) including the UK (Foreign, Commonwealth \\u0026amp; Development Office), China (Chinese Ministry of Agriculture and Rural Affairs), Australia (Australian Centre for International Agricultural Research), Canada (Agriculture and Agri-Food Canada), Netherlands (Directorate-General for International Cooperation) and Switzerland (Swiss Agency for Development and Cooperation). See https://www.cabi.org/about-cabi/who-we-work-with/key-donors/ for full details.\\u003c/p\\u003e\\n\\u003ch2\\u003eAuthors\\u0026apos; contributions\\u003c/h2\\u003e\\n\\u003cp\\u003eIF conceived and designed the study with input from RE, SP, and MC. Data collection was carried out by IF, PS, and KS. PS supervised all laboratory activities related to RNA metabarcoding. IF performed the data analysis and prepared the initial manuscript draft with input from MC. All co-authors reviewed the draft and contributed to the preparation of the final version of the manuscript\\u003c/p\\u003e\\n\\u003ch2\\u003eAcknowledgments\\u003c/h2\\u003e\\n\\u003cp\\u003eWe thank Beatrice Tolio, Delnia Sepahvand, and Diana Marčiulynienė for their valuable assistance in the laboratory. We thank Beat Ruffner for his guidance with processing ITS sequences from cultured fungal isolates.\\u003c/p\\u003e\"},{\"header\":\"References\",\"content\":\"\\u003col\\u003e\\u003cli\\u003e\\u003cspan\\u003eHubert B, Leprince O, Buitink J. Sleeping but not defenceless: seed dormancy and protection. J Exp Bot. 2024;75:6110\\u0026ndash;24. \\u003cspan class=\\\"ExternalRef\\\"\\u003e\\u003cspan class=\\\"RefSource\\\"\\u003ehttps://doi.org/10.1093/jxb/erae213\\u003c/span\\u003e\\u003cspan address=\\\"10.1093/jxb/erae213\\\" targettype=\\\"DOI\\\" class=\\\"RefTarget\\\"\\u003e\\u003c/span\\u003e\\u003c/span\\u003e.\\u003c/span\\u003e\\u003c/li\\u003e \\u003cli\\u003e\\u003cspan\\u003eNelson EB. The seed microbiome: Origins, interactions, and impacts. Plant Soil. 2018;422:7\\u0026ndash;34. \\u003cspan class=\\\"ExternalRef\\\"\\u003e\\u003cspan class=\\\"RefSource\\\"\\u003ehttps://doi.org/10.1007/s11104-017-3289-7\\u003c/span\\u003e\\u003cspan address=\\\"10.1007/s11104-017-3289-7\\\" targettype=\\\"DOI\\\" class=\\\"RefTarget\\\"\\u003e\\u003c/span\\u003e\\u003c/span\\u003e.\\u003c/span\\u003e\\u003c/li\\u003e \\u003cli\\u003e\\u003cspan\\u003eMittelstrass J, Heinzelmann R, Eschen R, Hartmann M, Kupper Q, Schneider S, et al. Metabarcoding with Illumina and Oxford Nanopore Technologies provides complementary insights into tree seed mycobiota. Environ Microbiome. 2025;20:53. \\u003cspan class=\\\"ExternalRef\\\"\\u003e\\u003cspan class=\\\"RefSource\\\"\\u003ehttps://doi.org/10.1186/s40793-025-00712-7\\u003c/span\\u003e\\u003cspan address=\\\"10.1186/s40793-025-00712-7\\\" targettype=\\\"DOI\\\" class=\\\"RefTarget\\\"\\u003e\\u003c/span\\u003e\\u003c/span\\u003e.\\u003c/span\\u003e\\u003c/li\\u003e \\u003cli\\u003e\\u003cspan\\u003eFranić I, Prospero S, Hartmann M, Allan E, Auger-Rozenberg M-A, Gr\\u0026uuml;nwald NJ, et al. Are traded forest tree seeds a potential source of nonnative pests? Ecol Appl. 2019;29:e01971. \\u003cspan class=\\\"ExternalRef\\\"\\u003e\\u003cspan class=\\\"RefSource\\\"\\u003ehttps://doi.org/10.1002/eap.1971\\u003c/span\\u003e\\u003cspan address=\\\"10.1002/eap.1971\\\" targettype=\\\"DOI\\\" class=\\\"RefTarget\\\"\\u003e\\u003c/span\\u003e\\u003c/span\\u003e.\\u003c/span\\u003e\\u003c/li\\u003e \\u003cli\\u003e\\u003cspan\\u003eFranić I, Cleary M, Aday Kaya AG, Bragan\\u0026ccedil;a H, Brodal G, Cech TL, et al. The Biosecurity Risks of International Forest Tree Seed Movements. Curr Forestry Rep. 2024;10:89\\u0026ndash;102. \\u003cspan class=\\\"ExternalRef\\\"\\u003e\\u003cspan class=\\\"RefSource\\\"\\u003ehttps://doi.org/10.1007/s40725-023-00211-3\\u003c/span\\u003e\\u003cspan address=\\\"10.1007/s40725-023-00211-3\\\" targettype=\\\"DOI\\\" class=\\\"RefTarget\\\"\\u003e\\u003c/span\\u003e\\u003c/span\\u003e.\\u003c/span\\u003e\\u003c/li\\u003e \\u003cli\\u003e\\u003cspan\\u003eCleary M, Oskay F, Doğmuş HT, Woodward S, Vettraino AM. Cryptic risks to forest biosecurity associated with the global movement of commercial seed. Forests. 2019;10:459. \\u003cspan class=\\\"ExternalRef\\\"\\u003e\\u003cspan class=\\\"RefSource\\\"\\u003ehttps://doi.org/10.3390/f10050459\\u003c/span\\u003e\\u003cspan address=\\\"10.3390/f10050459\\\" targettype=\\\"DOI\\\" class=\\\"RefTarget\\\"\\u003e\\u003c/span\\u003e\\u003c/span\\u003e.\\u003c/span\\u003e\\u003c/li\\u003e \\u003cli\\u003e\\u003cspan\\u003eR\\u0026eacute;tif F, Kunz C, Calabro K, Duval C, Prado S, Bailly C, et al. Seed fungal endophytes as biostimulants and biocontrol agents to improve seed performance. Front Plant Sci. 2023;14:1260292. \\u003cspan class=\\\"ExternalRef\\\"\\u003e\\u003cspan class=\\\"RefSource\\\"\\u003ehttps://doi.org/10.3389/fpls.2023.1260292\\u003c/span\\u003e\\u003cspan address=\\\"10.3389/fpls.2023.1260292\\\" targettype=\\\"DOI\\\" class=\\\"RefTarget\\\"\\u003e\\u003c/span\\u003e\\u003c/span\\u003e.\\u003c/span\\u003e\\u003c/li\\u003e \\u003cli\\u003e\\u003cspan\\u003edos Reis JBA, Lorenzi AS, do Vale HMM. Methods used for the study of endophytic fungi: a review on methodologies and challenges, and associated tips. Arch Microbiol. 2022;204:675. \\u003cspan class=\\\"ExternalRef\\\"\\u003e\\u003cspan class=\\\"RefSource\\\"\\u003ehttps://doi.org/10.1007/s00203-022-03283-0\\u003c/span\\u003e\\u003cspan address=\\\"10.1007/s00203-022-03283-0\\\" targettype=\\\"DOI\\\" class=\\\"RefTarget\\\"\\u003e\\u003c/span\\u003e\\u003c/span\\u003e.\\u003c/span\\u003e\\u003c/li\\u003e \\u003cli\\u003e\\u003cspan\\u003eGanley RJ, Newcombe G. Fungal endophytes in seeds and needles of \\u003cem\\u003ePinus monticola\\u003c/em\\u003e. Mycol Res. 2006;110:318\\u0026ndash;27. \\u003cspan class=\\\"ExternalRef\\\"\\u003e\\u003cspan class=\\\"RefSource\\\"\\u003ehttps://doi.org/10.1016/j.mycres.2005.10.005\\u003c/span\\u003e\\u003cspan address=\\\"10.1016/j.mycres.2005.10.005\\\" targettype=\\\"DOI\\\" class=\\\"RefTarget\\\"\\u003e\\u003c/span\\u003e\\u003c/span\\u003e.\\u003c/span\\u003e\\u003c/li\\u003e \\u003cli\\u003e\\u003cspan\\u003eHiggins KL, Arnold AE, Miadlikowska J, Sarvate SD, Lutzoni F. Phylogenetic relationships, host affinity, and geographic structure of boreal and arctic endophytes from three major plant lineages. Mol Phylogenet Evol. 2007;42:543\\u0026ndash;55. \\u003cspan class=\\\"ExternalRef\\\"\\u003e\\u003cspan class=\\\"RefSource\\\"\\u003ehttps://doi.org/10.1016/j.ympev.2006.07.012\\u003c/span\\u003e\\u003cspan address=\\\"10.1016/j.ympev.2006.07.012\\\" targettype=\\\"DOI\\\" class=\\\"RefTarget\\\"\\u003e\\u003c/span\\u003e\\u003c/span\\u003e.\\u003c/span\\u003e\\u003c/li\\u003e \\u003cli\\u003e\\u003cspan\\u003eZimmerman NB, Vitousek PM. Fungal endophyte communities reflect environmental structuring across a Hawaiian landscape. Proc Natl Acad Sci U S A. 2012;109:13022\\u0026ndash;7. \\u003cspan class=\\\"ExternalRef\\\"\\u003e\\u003cspan class=\\\"RefSource\\\"\\u003ehttps://doi.org/10.1073/pnas.1209872109\\u003c/span\\u003e\\u003cspan address=\\\"10.1073/pnas.1209872109\\\" targettype=\\\"DOI\\\" class=\\\"RefTarget\\\"\\u003e\\u003c/span\\u003e\\u003c/span\\u003e.\\u003c/span\\u003e\\u003c/li\\u003e \\u003cli\\u003e\\u003cspan\\u003eMapuranga J, Zhang N, Zhang L, Chang J, Yang W. Infection Strategies and Pathogenicity of Biotrophic Plant Fungal Pathogens. Front Microbiol. 2022;13:799396. \\u003cspan class=\\\"ExternalRef\\\"\\u003e\\u003cspan class=\\\"RefSource\\\"\\u003ehttps://doi.org/10.3389/fmicb.2022.799396\\u003c/span\\u003e\\u003cspan address=\\\"10.3389/fmicb.2022.799396\\\" targettype=\\\"DOI\\\" class=\\\"RefTarget\\\"\\u003e\\u003c/span\\u003e\\u003c/span\\u003e.\\u003c/span\\u003e\\u003c/li\\u003e \\u003cli\\u003e\\u003cspan\\u003eEschen R, De Groot M, Glavendekić M, Lacković N, Matošević D, Morales-Rodriguez C, et al. Spotting the pests of tomorrow\\u0026mdash;Sampling designs for detection of species associations with woody plants. J Biogeogr. 2019;46:2159\\u0026ndash;73. \\u003cspan class=\\\"ExternalRef\\\"\\u003e\\u003cspan class=\\\"RefSource\\\"\\u003ehttps://doi.org/10.1111/jbi.13670\\u003c/span\\u003e\\u003cspan address=\\\"10.1111/jbi.13670\\\" targettype=\\\"DOI\\\" class=\\\"RefTarget\\\"\\u003e\\u003c/span\\u003e\\u003c/span\\u003e.\\u003c/span\\u003e\\u003c/li\\u003e \\u003cli\\u003e\\u003cspan\\u003eNilsson RH, Anslan S, Bahram M, Wurzbacher C, Baldrian P, Tedersoo L. Mycobiome diversity: high-throughput sequencing and identification of fungi. Nat Rev Microbiol. 2019;17:95\\u0026ndash;109. \\u003cspan class=\\\"ExternalRef\\\"\\u003e\\u003cspan class=\\\"RefSource\\\"\\u003ehttps://doi.org/10.1038/s41579-018-0116-y\\u003c/span\\u003e\\u003cspan address=\\\"10.1038/s41579-018-0116-y\\\" targettype=\\\"DOI\\\" class=\\\"RefTarget\\\"\\u003e\\u003c/span\\u003e\\u003c/span\\u003e.\\u003c/span\\u003e\\u003c/li\\u003e \\u003cli\\u003e\\u003cspan\\u003eTedersoo L, Drenkhan R, Anslan S, Morales-Rodriguez C, Cleary M. High-throughput identification and diagnostics of pathogens and pests: overview and practical recommendations. Mol Ecol Resour. 2018;19:47\\u0026ndash;76. \\u003cspan class=\\\"ExternalRef\\\"\\u003e\\u003cspan class=\\\"RefSource\\\"\\u003ehttps://doi.org/10.1111/1755-0998.12959\\u003c/span\\u003e\\u003cspan address=\\\"10.1111/1755-0998.12959\\\" targettype=\\\"DOI\\\" class=\\\"RefTarget\\\"\\u003e\\u003c/span\\u003e\\u003c/span\\u003e.\\u003c/span\\u003e\\u003c/li\\u003e \\u003cli\\u003e\\u003cspan\\u003eSchoch CL, Seifert KA, Huhndorf S, Robert V, Spouge JL, Levesque CA, et al. Nuclear ribosomal internal transcribed spacer (ITS) region as a universal DNA barcode marker for Fungi. Proc Natl Acad Sci U S A. 2012;109:6241\\u0026ndash;6. \\u003cspan class=\\\"ExternalRef\\\"\\u003e\\u003cspan class=\\\"RefSource\\\"\\u003ehttps://doi.org/10.1073/pnas.1117018109\\u003c/span\\u003e\\u003cspan address=\\\"10.1073/pnas.1117018109\\\" targettype=\\\"DOI\\\" class=\\\"RefTarget\\\"\\u003e\\u003c/span\\u003e\\u003c/span\\u003e.\\u003c/span\\u003e\\u003c/li\\u003e \\u003cli\\u003e\\u003cspan\\u003eCarini P, Marsden PJ, Leff JW, Morgan EE, Strickland MS, Fierer N. Relic DNA is abundant in soil and obscures estimates of soil microbial diversity. Nat Microbiol. 2017;2:16242. \\u003cspan class=\\\"ExternalRef\\\"\\u003e\\u003cspan class=\\\"RefSource\\\"\\u003ehttps://doi.org/10.1038/nmicrobiol.2016.242\\u003c/span\\u003e\\u003cspan address=\\\"10.1038/nmicrobiol.2016.242\\\" targettype=\\\"DOI\\\" class=\\\"RefTarget\\\"\\u003e\\u003c/span\\u003e\\u003c/span\\u003e.\\u003c/span\\u003e\\u003c/li\\u003e \\u003cli\\u003e\\u003cspan\\u003ePochon X, Zaiko A, Fletcher LM, Laroche O, Wood SA. Wanted dead or alive? Using metabarcoding of environmental DNA and RNA to distinguish living assemblages for biosecurity applications. PLoS ONE. 2017;12:e0187636. \\u003cspan class=\\\"ExternalRef\\\"\\u003e\\u003cspan class=\\\"RefSource\\\"\\u003ehttps://doi.org/10.1371/journal.pone.0187636\\u003c/span\\u003e\\u003cspan address=\\\"10.1371/journal.pone.0187636\\\" targettype=\\\"DOI\\\" class=\\\"RefTarget\\\"\\u003e\\u003c/span\\u003e\\u003c/span\\u003e.\\u003c/span\\u003e\\u003c/li\\u003e \\u003cli\\u003e\\u003cspan\\u003eWong B, Leal I, Feau N, Dale A, Uzunovic A, Hamelin RC. Molecular assays to detect the presence and viability of \\u003cem\\u003ePhytophthora ramorum\\u003c/em\\u003e and \\u003cem\\u003eGrosmannia clavigera\\u003c/em\\u003e. PLoS ONE. 2020;15:e0221742. \\u003cspan class=\\\"ExternalRef\\\"\\u003e\\u003cspan class=\\\"RefSource\\\"\\u003ehttps://doi.org/10.1371/journal.pone.0221742\\u003c/span\\u003e\\u003cspan address=\\\"10.1371/journal.pone.0221742\\\" targettype=\\\"DOI\\\" class=\\\"RefTarget\\\"\\u003e\\u003c/span\\u003e\\u003c/span\\u003e.\\u003c/span\\u003e\\u003c/li\\u003e \\u003cli\\u003e\\u003cspan\\u003eLeal I, Feau N, Uzunovic A, Foord B, Hamelin RC. A Molecular Method to Assess Viability of \\u003cem\\u003ePhytophthora\\u003c/em\\u003e in Infected Wood Following Phytosanitary Heat Treatment. PhytoFrontiersTM. 2024;4:690\\u0026ndash;9. \\u003cspan class=\\\"ExternalRef\\\"\\u003e\\u003cspan class=\\\"RefSource\\\"\\u003ehttps://doi.org/10.1094/phytofr-05-24-0056-r\\u003c/span\\u003e\\u003cspan address=\\\"10.1094/phytofr-05-24-0056-r\\\" targettype=\\\"DOI\\\" class=\\\"RefTarget\\\"\\u003e\\u003c/span\\u003e\\u003c/span\\u003e.\\u003c/span\\u003e\\u003c/li\\u003e \\u003cli\\u003e\\u003cspan\\u003eAnderson IC, Parkin PI. Detection of active soil fungi by RT-PCR amplification of precursor rRNA molecules. J Microbiol Methods. 2007;68:248\\u0026ndash;53. \\u003cspan class=\\\"ExternalRef\\\"\\u003e\\u003cspan class=\\\"RefSource\\\"\\u003ehttps://doi.org/10.1016/j.mimet.2006.08.005\\u003c/span\\u003e\\u003cspan address=\\\"10.1016/j.mimet.2006.08.005\\\" targettype=\\\"DOI\\\" class=\\\"RefTarget\\\"\\u003e\\u003c/span\\u003e\\u003c/span\\u003e.\\u003c/span\\u003e\\u003c/li\\u003e \\u003cli\\u003e\\u003cspan\\u003eYates MC, Derry AM, Cristescu ME, Environmental RNA. A Revolution in Ecological Resolution? Trends in Ecology and Evolution. 2021;36:601\\u0026ndash;9. \\u003cspan class=\\\"ExternalRef\\\"\\u003e\\u003cspan class=\\\"RefSource\\\"\\u003ehttps://doi.org/10.1016/j.tree.2021.03.001\\u003c/span\\u003e\\u003cspan address=\\\"10.1016/j.tree.2021.03.001\\\" targettype=\\\"DOI\\\" class=\\\"RefTarget\\\"\\u003e\\u003c/span\\u003e\\u003c/span\\u003e\\u003c/span\\u003e\\u003c/li\\u003e \\u003cli\\u003e\\u003cspan\\u003eBaldrian P, Kolaiř\\u0026iacute;k M, Štursov\\u0026aacute; M, Kopeck\\u0026yacute; J, Val\\u0026aacute;škov\\u0026aacute; V, Větrovsk\\u0026yacute; T, et al. Active and total microbial communities in forest soil are largely different and highly stratified during decomposition. ISME J. 2012;6:248\\u0026ndash;58. \\u003cspan class=\\\"ExternalRef\\\"\\u003e\\u003cspan class=\\\"RefSource\\\"\\u003ehttps://doi.org/10.1038/ismej.2011.95\\u003c/span\\u003e\\u003cspan address=\\\"10.1038/ismej.2011.95\\\" targettype=\\\"DOI\\\" class=\\\"RefTarget\\\"\\u003e\\u003c/span\\u003e\\u003c/span\\u003e.\\u003c/span\\u003e\\u003c/li\\u003e \\u003cli\\u003e\\u003cspan\\u003eRajala T, Peltoniemi M, Hantula J, M\\u0026auml;kip\\u0026auml;\\u0026auml; R, Pennanen T. RNA reveals a succession of active fungi during the decay of Norway spruce logs. Fungal Ecol. 2011;4:437\\u0026ndash;48. \\u003cspan class=\\\"ExternalRef\\\"\\u003e\\u003cspan class=\\\"RefSource\\\"\\u003ehttps://doi.org/10.1016/j.funeco.2011.05.005\\u003c/span\\u003e\\u003cspan address=\\\"10.1016/j.funeco.2011.05.005\\\" targettype=\\\"DOI\\\" class=\\\"RefTarget\\\"\\u003e\\u003c/span\\u003e\\u003c/span\\u003e.\\u003c/span\\u003e\\u003c/li\\u003e \\u003cli\\u003e\\u003cspan\\u003eXue Z, Tian W, Han Y, Li S, Guo J, He H, et al. Environmental RNA metabarcoding for ballast water microbial diversity: Minimizing false positives. Sci Total Environ. 2024;955:176902. \\u003cspan class=\\\"ExternalRef\\\"\\u003e\\u003cspan class=\\\"RefSource\\\"\\u003ehttps://doi.org/10.1016/j.scitotenv.2024.176902\\u003c/span\\u003e\\u003cspan address=\\\"10.1016/j.scitotenv.2024.176902\\\" targettype=\\\"DOI\\\" class=\\\"RefTarget\\\"\\u003e\\u003c/span\\u003e\\u003c/span\\u003e.\\u003c/span\\u003e\\u003c/li\\u003e \\u003cli\\u003e\\u003cspan\\u003eKong H, Yang EJ, Jiao N, Lee Y, Jung J, Cho KH, et al. RNA outperforms DNA-based metabarcoding in assessing the diversity and response of microeukaryotes to environmental variables in the Arctic Ocean. Sci Total Environ. 2023;876:162608. \\u003cspan class=\\\"ExternalRef\\\"\\u003e\\u003cspan class=\\\"RefSource\\\"\\u003ehttps://doi.org/10.1016/j.scitotenv.2023.162608\\u003c/span\\u003e\\u003cspan address=\\\"10.1016/j.scitotenv.2023.162608\\\" targettype=\\\"DOI\\\" class=\\\"RefTarget\\\"\\u003e\\u003c/span\\u003e\\u003c/span\\u003e.\\u003c/span\\u003e\\u003c/li\\u003e \\u003cli\\u003e\\u003cspan\\u003eNguyen D, Boberg J, Ihrmark K, Stenstr\\u0026ouml;m E, Stenlid J. Do foliar fungal communities of Norway spruce shift along a tree species diversity gradient in mature European forests? Fungal Ecol. 2016;23:97\\u0026ndash;108. \\u003cspan class=\\\"ExternalRef\\\"\\u003e\\u003cspan class=\\\"RefSource\\\"\\u003ehttps://doi.org/10.1016/j.funeco.2016.07.003\\u003c/span\\u003e\\u003cspan address=\\\"10.1016/j.funeco.2016.07.003\\\" targettype=\\\"DOI\\\" class=\\\"RefTarget\\\"\\u003e\\u003c/span\\u003e\\u003c/span\\u003e.\\u003c/span\\u003e\\u003c/li\\u003e \\u003cli\\u003e\\u003cspan\\u003eWutkowska M, Vader A, Mundra S, Cooper EJ, Eidesen PB. Dead or alive; Or does it really matter? Level of congruency between trophic modes in total and active fungal communities in high arctic soil. Front Microbiol. 2019;9:3243. \\u003cspan class=\\\"ExternalRef\\\"\\u003e\\u003cspan class=\\\"RefSource\\\"\\u003ehttps://doi.org/10.3389/fmicb.2018.03243\\u003c/span\\u003e\\u003cspan address=\\\"10.3389/fmicb.2018.03243\\\" targettype=\\\"DOI\\\" class=\\\"RefTarget\\\"\\u003e\\u003c/span\\u003e\\u003c/span\\u003e.\\u003c/span\\u003e\\u003c/li\\u003e \\u003cli\\u003e\\u003cspan\\u003eCamenzind T, Weimershaus P, Lehmann A, Aguilar-Trigueros C, Rillig MC. Soil fungi invest into asexual sporulation under resource scarcity, but trait spaces of individual isolates are unique. Environ Microbiol. 2022;24:2962\\u0026ndash;78. \\u003cspan class=\\\"ExternalRef\\\"\\u003e\\u003cspan class=\\\"RefSource\\\"\\u003ehttps://doi.org/10.1111/1462-2920.16012\\u003c/span\\u003e\\u003cspan address=\\\"10.1111/1462-2920.16012\\\" targettype=\\\"DOI\\\" class=\\\"RefTarget\\\"\\u003e\\u003c/span\\u003e\\u003c/span\\u003e.\\u003c/span\\u003e\\u003c/li\\u003e \\u003cli\\u003e\\u003cspan\\u003eChen DV, Slowinski SP, Kido AK, Bruns EL. High temperatures reduce growth, infection, and transmission of a naturally occurring fungal plant pathogen. Ecology. 2024;105:e4373. \\u003cspan class=\\\"ExternalRef\\\"\\u003e\\u003cspan class=\\\"RefSource\\\"\\u003ehttps://doi.org/10.1002/ecy.4373\\u003c/span\\u003e\\u003cspan address=\\\"10.1002/ecy.4373\\\" targettype=\\\"DOI\\\" class=\\\"RefTarget\\\"\\u003e\\u003c/span\\u003e\\u003c/span\\u003e.\\u003c/span\\u003e\\u003c/li\\u003e \\u003cli\\u003e\\u003cspan\\u003eJo T, Tsuri K, Hirohara T, Yamanaka H. Warm temperature and alkaline conditions accelerate environmental RNA degradation. Environ DNA. 2023;5:836\\u0026ndash;48. \\u003cspan class=\\\"ExternalRef\\\"\\u003e\\u003cspan class=\\\"RefSource\\\"\\u003ehttps://doi.org/10.1002/edn3.334\\u003c/span\\u003e\\u003cspan address=\\\"10.1002/edn3.334\\\" targettype=\\\"DOI\\\" class=\\\"RefTarget\\\"\\u003e\\u003c/span\\u003e\\u003c/span\\u003e.\\u003c/span\\u003e\\u003c/li\\u003e \\u003cli\\u003e\\u003cspan\\u003eGamboa MA, Laureano S, Bayman P. Measuring diversity of endophytic fungi in leaf fragments: Does size matter? Mycopathologia. 2003;156:41\\u0026ndash;5. \\u003cspan class=\\\"ExternalRef\\\"\\u003e\\u003cspan class=\\\"RefSource\\\"\\u003ehttps://doi.org/10.1023/A:1021362217723\\u003c/span\\u003e\\u003cspan address=\\\"10.1023/A:1021362217723\\\" targettype=\\\"DOI\\\" class=\\\"RefTarget\\\"\\u003e\\u003c/span\\u003e\\u003c/span\\u003e\\u003c/span\\u003e\\u003c/li\\u003e \\u003cli\\u003e\\u003cspan\\u003eWhite T, Bruns T, Lee S, Taylor J. Amplification and direct sequencing of fungal ribosomal RNA genes for phylogenetics. In: Innis MA, Gelfand DH, Sninsky JJ WT, editors. PCR Protocols: A Guide to Methods and Applications. New York: Academic Press Inc; 1990. pp. 315\\u0026ndash;22.\\u003c/span\\u003e\\u003c/li\\u003e \\u003cli\\u003e\\u003cspan\\u003eAltschup SF, Gish W, Miller W, Myers EW, Lipman DJ. Basic Local Alignment Search Tool. J Mol Biol. 1990;215:403\\u0026ndash;10.\\u003c/span\\u003e\\u003c/li\\u003e \\u003cli\\u003e\\u003cspan\\u003eZhang Z, Schwartz S, Wagner L, Miller W. A Greedy Algorithm for Aligning DNA Sequences. J Comput Biol. 2000;7:203\\u0026ndash;14.\\u003c/span\\u003e\\u003c/li\\u003e \\u003cli\\u003e\\u003cspan\\u003eBolyen E, Rideout JR, Dillon MR, Bokulich NA, Abnet CC, Al-Ghalith GA, et al. Reproducible, interactive, scalable and extensible microbiome data science using QIIME 2. Nat Biotechnol. 2019;37:852\\u0026ndash;7. \\u003cspan class=\\\"ExternalRef\\\"\\u003e\\u003cspan class=\\\"RefSource\\\"\\u003ehttps://doi.org/10.1038/s41587-019-0209-9\\u003c/span\\u003e\\u003cspan address=\\\"10.1038/s41587-019-0209-9\\\" targettype=\\\"DOI\\\" class=\\\"RefTarget\\\"\\u003e\\u003c/span\\u003e\\u003c/span\\u003e.\\u003c/span\\u003e\\u003c/li\\u003e \\u003cli\\u003e\\u003cspan\\u003eMartin M. Cutadapt removes adapter sequences from high-throughput sequencing reads. EMBnet J. 2011;17:10\\u0026ndash;2. \\u003cspan class=\\\"ExternalRef\\\"\\u003e\\u003cspan class=\\\"RefSource\\\"\\u003ehttps://doi.org/https://doi.org/10.14806/ej.17.1.200\\u003c/span\\u003e\\u003cspan address=\\\"10.14806/ej.17.1.200\\\" targettype=\\\"DOI\\\" class=\\\"RefTarget\\\"\\u003e\\u003c/span\\u003e\\u003c/span\\u003e.\\u003c/span\\u003e\\u003c/li\\u003e \\u003cli\\u003e\\u003cspan\\u003eCallahan BJ, McMurdie PJ, Rosen MJ, Han AW, Johnson AJA, Holmes SP. DADA2: High-resolution sample inference from Illumina amplicon data. Nat Methods. 2016;13:581\\u0026ndash;3. \\u003cspan class=\\\"ExternalRef\\\"\\u003e\\u003cspan class=\\\"RefSource\\\"\\u003ehttps://doi.org/10.1038/nmeth.3869\\u003c/span\\u003e\\u003cspan address=\\\"10.1038/nmeth.3869\\\" targettype=\\\"DOI\\\" class=\\\"RefTarget\\\"\\u003e\\u003c/span\\u003e\\u003c/span\\u003e.\\u003c/span\\u003e\\u003c/li\\u003e \\u003cli\\u003e\\u003cspan\\u003eRognes T, Flouri T, Nichols B, Quince C, Mah\\u0026eacute; F. VSEARCH: A versatile open source tool for metagenomics. PeerJ. 2016;10:e2584. \\u003cspan class=\\\"ExternalRef\\\"\\u003e\\u003cspan class=\\\"RefSource\\\"\\u003ehttps://doi.org/10.7717/peerj.2584\\u003c/span\\u003e\\u003cspan address=\\\"10.7717/peerj.2584\\\" targettype=\\\"DOI\\\" class=\\\"RefTarget\\\"\\u003e\\u003c/span\\u003e\\u003c/span\\u003e.\\u003c/span\\u003e\\u003c/li\\u003e \\u003cli\\u003e\\u003cspan\\u003ePedregosa F, Varoquaux G, Gramfort A, Michel V, Thirion B, Grisel O, et al. Scikit-learn: Machine Learning in Python. J Mach Learn Res. 2011;12:2825\\u0026ndash;30.\\u003c/span\\u003e\\u003c/li\\u003e \\u003cli\\u003e\\u003cspan\\u003eK\\u0026otilde;ljalg U, Nilsson RH, Abarenkov K, Tedersoo L, Taylor AFS, Bahram M, et al. Towards a unified paradigm for sequence-based identification of fungi. Mol Ecol. 2013;22:5271\\u0026ndash;7. \\u003cspan class=\\\"ExternalRef\\\"\\u003e\\u003cspan class=\\\"RefSource\\\"\\u003ehttps://doi.org/10.1111/mec.12481\\u003c/span\\u003e\\u003cspan address=\\\"10.1111/mec.12481\\\" targettype=\\\"DOI\\\" class=\\\"RefTarget\\\"\\u003e\\u003c/span\\u003e\\u003c/span\\u003e.\\u003c/span\\u003e\\u003c/li\\u003e \\u003cli\\u003e\\u003cspan\\u003eMonard C, Gantner S, Stenlid J. Utilizing ITS1 and ITS2 to study environmental fungal diversity using pyrosequencing. FEMS Microbiol Ecol. 2013;84:165\\u0026ndash;75. \\u003cspan class=\\\"ExternalRef\\\"\\u003e\\u003cspan class=\\\"RefSource\\\"\\u003ehttps://doi.org/10.1111/1574-6941.12046\\u003c/span\\u003e\\u003cspan address=\\\"10.1111/1574-6941.12046\\\" targettype=\\\"DOI\\\" class=\\\"RefTarget\\\"\\u003e\\u003c/span\\u003e\\u003c/span\\u003e.\\u003c/span\\u003e\\u003c/li\\u003e \\u003cli\\u003e\\u003cspan\\u003eP\\u0026otilde;lme S, Abarenkov K, Henrik Nilsson R, Lindahl BD, Clemmensen KE, Kauserud H, et al. FungalTraits: a user-friendly traits database of fungi and fungus-like stramenopiles. Fungal Divers. 2020;105:1\\u0026ndash;16. \\u003cspan class=\\\"ExternalRef\\\"\\u003e\\u003cspan class=\\\"RefSource\\\"\\u003ehttps://doi.org/10.1007/s13225-020-00466-2\\u003c/span\\u003e\\u003cspan address=\\\"10.1007/s13225-020-00466-2\\\" targettype=\\\"DOI\\\" class=\\\"RefTarget\\\"\\u003e\\u003c/span\\u003e\\u003c/span\\u003e.\\u003c/span\\u003e\\u003c/li\\u003e \\u003cli\\u003e\\u003cspan\\u003eR Core Team. R: A Language and Environment for Statistical Computing. 2024.\\u003c/span\\u003e\\u003c/li\\u003e \\u003cli\\u003e\\u003cspan\\u003eChen J, Zhang X, Yang L, Zhang L, _GUniFrac. Generalized UniFrac Distances, Distance-Based Multivariate Methods and Feature-Based Univariate Methods for Microbiome Data Analysis_. 2025.\\u003c/span\\u003e\\u003c/li\\u003e \\u003cli\\u003e\\u003cspan\\u003eSchloss PD. Rarefaction is currently the best approach to control for uneven sequencing effort in amplicon sequence analyses. mSphere. 2024;9:e00354\\u0026ndash;23. \\u003cspan class=\\\"ExternalRef\\\"\\u003e\\u003cspan class=\\\"RefSource\\\"\\u003ehttps://doi.org/10.1128/msphere.00354-23\\u003c/span\\u003e\\u003cspan address=\\\"10.1128/msphere.00354-23\\\" targettype=\\\"DOI\\\" class=\\\"RefTarget\\\"\\u003e\\u003c/span\\u003e\\u003c/span\\u003e.\\u003c/span\\u003e\\u003c/li\\u003e \\u003cli\\u003e\\u003cspan\\u003eCassol I, Iba\\u0026ntilde;ez M, Bustamante JP. Key features and guidelines for the application of microbial alpha diversity metrics. Sci Rep. 2025;15:622. \\u003cspan class=\\\"ExternalRef\\\"\\u003e\\u003cspan class=\\\"RefSource\\\"\\u003ehttps://doi.org/10.1038/s41598-024-77864-y\\u003c/span\\u003e\\u003cspan address=\\\"10.1038/s41598-024-77864-y\\\" targettype=\\\"DOI\\\" class=\\\"RefTarget\\\"\\u003e\\u003c/span\\u003e\\u003c/span\\u003e.\\u003c/span\\u003e\\u003c/li\\u003e \\u003cli\\u003e\\u003cspan\\u003eJost L. Entropy and diversity. Oikos. 2006;113:363\\u0026ndash;75. \\u003cspan class=\\\"ExternalRef\\\"\\u003e\\u003cspan class=\\\"RefSource\\\"\\u003ehttps://doi.org/10.1111/j.2006.0030-1299.14714.x\\u003c/span\\u003e\\u003cspan address=\\\"10.1111/j.2006.0030-1299.14714.x\\\" targettype=\\\"DOI\\\" class=\\\"RefTarget\\\"\\u003e\\u003c/span\\u003e\\u003c/span\\u003e.\\u003c/span\\u003e\\u003c/li\\u003e \\u003cli\\u003e\\u003cspan\\u003eChao A, Chiu C-H, Jost L. Unifying Species Diversity, Phylogenetic Diversity, Functional Diversity, and Related Similarity and Differentiation Measures Through Hill Numbers. Annu Rev Ecol Evol Syst. 2014;45:297\\u0026ndash;324. \\u003cspan class=\\\"ExternalRef\\\"\\u003e\\u003cspan class=\\\"RefSource\\\"\\u003ehttps://doi.org/10.1146/annurev-ecolsys-120213-091540\\u003c/span\\u003e\\u003cspan address=\\\"10.1146/annurev-ecolsys-120213-091540\\\" targettype=\\\"DOI\\\" class=\\\"RefTarget\\\"\\u003e\\u003c/span\\u003e\\u003c/span\\u003e.\\u003c/span\\u003e\\u003c/li\\u003e \\u003cli\\u003e\\u003cspan\\u003eVenables WN, Ripley BD. Modern Applied Statistics with S. Fourth edition. New York: Springer; 2002. \\u003cspan class=\\\"ExternalRef\\\"\\u003e\\u003cspan class=\\\"RefSource\\\"\\u003ehttps://doi.org/10.1007/978-0-387-21706-2\\u003c/span\\u003e\\u003cspan address=\\\"10.1007/978-0-387-21706-2\\\" targettype=\\\"DOI\\\" class=\\\"RefTarget\\\"\\u003e\\u003c/span\\u003e\\u003c/span\\u003e. https://doi.org/.\\u003c/span\\u003e\\u003c/li\\u003e \\u003cli\\u003e\\u003cspan\\u003eFox J, Weisberg S, An R. Companion to Applied Regression. Third. Thousand Oaks: Sage; 2019.\\u003c/span\\u003e\\u003c/li\\u003e \\u003cli\\u003e\\u003cspan\\u003eLenth RV, Least-Squares, Means. The R Package lsmeans. J Stat Softw. 2016;69:1\\u0026ndash;33. \\u003cspan class=\\\"ExternalRef\\\"\\u003e\\u003cspan class=\\\"RefSource\\\"\\u003ehttps://doi.org/10.18637/jss.v069.i01\\u003c/span\\u003e\\u003cspan address=\\\"10.18637/jss.v069.i01\\\" targettype=\\\"DOI\\\" class=\\\"RefTarget\\\"\\u003e\\u003c/span\\u003e\\u003c/span\\u003e.\\u003c/span\\u003e\\u003c/li\\u003e \\u003cli\\u003e\\u003cspan\\u003eOksanen J, Blanchet FG, Friendly M, Kindt R, Legendre P, Dan, McGlinn et al. _vegan: Community Ecology Package_. R package version 2.5-3. 2025.\\u003c/span\\u003e\\u003c/li\\u003e \\u003cli\\u003e\\u003cspan\\u003eMcMurdie PJ, Holmes S. Waste Not, Want Not: Why Rarefying Microbiome Data Is Inadmissible. PLoS Comput Biol. 2014;10:e1003531. \\u003cspan class=\\\"ExternalRef\\\"\\u003e\\u003cspan class=\\\"RefSource\\\"\\u003ehttps://doi.org/10.1371/journal.pcbi.1003531\\u003c/span\\u003e\\u003cspan address=\\\"10.1371/journal.pcbi.1003531\\\" targettype=\\\"DOI\\\" class=\\\"RefTarget\\\"\\u003e\\u003c/span\\u003e\\u003c/span\\u003e.\\u003c/span\\u003e\\u003c/li\\u003e \\u003cli\\u003e\\u003cspan\\u003eC\\u0026aacute;ceres M, De, Legendre P. Associations between species and groups of sites: indices and statistical inference. Ecology. 2009;90:3566\\u0026ndash;74. \\u003cspan class=\\\"ExternalRef\\\"\\u003e\\u003cspan class=\\\"RefSource\\\"\\u003ehttps://doi.org/10.1890/08-1823.1\\u003c/span\\u003e\\u003cspan address=\\\"10.1890/08-1823.1\\\" targettype=\\\"DOI\\\" class=\\\"RefTarget\\\"\\u003e\\u003c/span\\u003e\\u003c/span\\u003e.\\u003c/span\\u003e\\u003c/li\\u003e \\u003cli\\u003e\\u003cspan\\u003eChen Y, Zhu X, Hou Z, Wang Y, Zhou Y, Wang L, et al. RNA-Based Analysis Reveals High Diversity of Plant-Associated Active Fungi in the Atmosphere. Front Microbiol. 2021;12:683266. \\u003cspan class=\\\"ExternalRef\\\"\\u003e\\u003cspan class=\\\"RefSource\\\"\\u003ehttps://doi.org/10.3389/fmicb.2021.683266\\u003c/span\\u003e\\u003cspan address=\\\"10.3389/fmicb.2021.683266\\\" targettype=\\\"DOI\\\" class=\\\"RefTarget\\\"\\u003e\\u003c/span\\u003e\\u003c/span\\u003e.\\u003c/span\\u003e\\u003c/li\\u003e \\u003cli\\u003e\\u003cspan\\u003eJones SE, Lennon JT. Dormancy contributes to the maintenance of microbial diversity. Proc Natl Acad Sci U S A. 2010;107:5881\\u0026ndash;6. \\u003cspan class=\\\"ExternalRef\\\"\\u003e\\u003cspan class=\\\"RefSource\\\"\\u003ehttps://doi.org/10.1073/pnas.0912765107\\u003c/span\\u003e\\u003cspan address=\\\"10.1073/pnas.0912765107\\\" targettype=\\\"DOI\\\" class=\\\"RefTarget\\\"\\u003e\\u003c/span\\u003e\\u003c/span\\u003e.\\u003c/span\\u003e\\u003c/li\\u003e \\u003cli\\u003e\\u003cspan\\u003eCamenzind T, Lehmann A, Ahland J, Rumpel S, Rillig MC. Trait-based approaches reveal fungal adaptations to nutrient-limiting conditions. Environ Microbiol. 2020;22:3548\\u0026ndash;60. \\u003cspan class=\\\"ExternalRef\\\"\\u003e\\u003cspan class=\\\"RefSource\\\"\\u003ehttps://doi.org/10.1111/1462-2920.15132\\u003c/span\\u003e\\u003cspan address=\\\"10.1111/1462-2920.15132\\\" targettype=\\\"DOI\\\" class=\\\"RefTarget\\\"\\u003e\\u003c/span\\u003e\\u003c/span\\u003e.\\u003c/span\\u003e\\u003c/li\\u003e \\u003cli\\u003e\\u003cspan\\u003eKlein AM, Bohannan BJM, Jaffe DA, Levin DA, Green JL. Molecular evidence for metabolically active bacteria in the atmosphere. Front Microbiol. 2016;7:772. \\u003cspan class=\\\"ExternalRef\\\"\\u003e\\u003cspan class=\\\"RefSource\\\"\\u003ehttps://doi.org/10.3389/fmicb.2016.00772\\u003c/span\\u003e\\u003cspan address=\\\"10.3389/fmicb.2016.00772\\\" targettype=\\\"DOI\\\" class=\\\"RefTarget\\\"\\u003e\\u003c/span\\u003e\\u003c/span\\u003e.\\u003c/span\\u003e\\u003c/li\\u003e \\u003cli\\u003e\\u003cspan\\u003eCampbell BJ, Yu L, Heidelberg JF, Kirchman DL. Activity of abundant and rare bacteria in a coastal ocean. Proc Natl Acad Sci U S A. 2011;108:12776\\u0026ndash;81. \\u003cspan class=\\\"ExternalRef\\\"\\u003e\\u003cspan class=\\\"RefSource\\\"\\u003ehttps://doi.org/10.1073/pnas.1101405108\\u003c/span\\u003e\\u003cspan address=\\\"10.1073/pnas.1101405108\\\" targettype=\\\"DOI\\\" class=\\\"RefTarget\\\"\\u003e\\u003c/span\\u003e\\u003c/span\\u003e.\\u003c/span\\u003e\\u003c/li\\u003e \\u003cli\\u003e\\u003cspan\\u003ePu Y, Zhou X, Cai H, Lou T, Liu C, Kong M, et al. Impact of DNA Extraction Methods on Gut Microbiome Profiles: A Comparative Metagenomic Study. Phenomics. 2025;5:76\\u0026ndash;90. \\u003cspan class=\\\"ExternalRef\\\"\\u003e\\u003cspan class=\\\"RefSource\\\"\\u003ehttps://doi.org/10.1007/s43657-025-00232-x\\u003c/span\\u003e\\u003cspan address=\\\"10.1007/s43657-025-00232-x\\\" targettype=\\\"DOI\\\" class=\\\"RefTarget\\\"\\u003e\\u003c/span\\u003e\\u003c/span\\u003e.\\u003c/span\\u003e\\u003c/li\\u003e \\u003cli\\u003e\\u003cspan\\u003eMcCarthy 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. \\u003cspan class=\\\"ExternalRef\\\"\\u003e\\u003cspan class=\\\"RefSource\\\"\\u003ehttps://doi.org/10.1371/journal.pone.0121659\\u003c/span\\u003e\\u003cspan address=\\\"10.1371/journal.pone.0121659\\\" targettype=\\\"DOI\\\" class=\\\"RefTarget\\\"\\u003e\\u003c/span\\u003e\\u003c/span\\u003e.\\u003c/span\\u003e\\u003c/li\\u003e \\u003cli\\u003e\\u003cspan\\u003eOskay F, Vettraino AM, Doğmuş HT, Lehtij\\u0026auml;rvi A, Woodward S, Cleary M. Seed quantity affects the fungal community composition detected using metabarcoding. Sci Rep. 2022;12:3060. \\u003cspan class=\\\"ExternalRef\\\"\\u003e\\u003cspan class=\\\"RefSource\\\"\\u003ehttps://doi.org/10.1038/s41598-022-06997-9\\u003c/span\\u003e\\u003cspan address=\\\"10.1038/s41598-022-06997-9\\\" targettype=\\\"DOI\\\" class=\\\"RefTarget\\\"\\u003e\\u003c/span\\u003e\\u003c/span\\u003e.\\u003c/span\\u003e\\u003c/li\\u003e \\u003cli\\u003e\\u003cspan\\u003eWalker SC, Poos MS, Jackson DA. Functional rarefaction: Estimating functional diversity from field data. Oikos. 2008;117:286\\u0026ndash;96. \\u003cspan class=\\\"ExternalRef\\\"\\u003e\\u003cspan class=\\\"RefSource\\\"\\u003ehttps://doi.org/10.1111/j.2007.0030-1299.16171.x\\u003c/span\\u003e\\u003cspan address=\\\"10.1111/j.2007.0030-1299.16171.x\\\" targettype=\\\"DOI\\\" class=\\\"RefTarget\\\"\\u003e\\u003c/span\\u003e\\u003c/span\\u003e.\\u003c/span\\u003e\\u003c/li\\u003e \\u003cli\\u003e\\u003cspan\\u003eMcCoy CO, Matsen FA. Abundance-weighted phylogenetic diversity measures distinguish microbial community states and are robust to sampling depth. PeerJ. 2013;1:e157. \\u003cspan class=\\\"ExternalRef\\\"\\u003e\\u003cspan class=\\\"RefSource\\\"\\u003ehttps://doi.org/10.7717/peerj.157\\u003c/span\\u003e\\u003cspan address=\\\"10.7717/peerj.157\\\" targettype=\\\"DOI\\\" class=\\\"RefTarget\\\"\\u003e\\u003c/span\\u003e\\u003c/span\\u003e.\\u003c/span\\u003e\\u003c/li\\u003e \\u003cli\\u003e\\u003cspan\\u003eVerwilt J, Mestdagh P, Vandesompele JO. Artifacts and biases of the reverse transcription reaction in RNA sequencing. RNA. 2023;29:889\\u0026ndash;97. \\u003cspan class=\\\"ExternalRef\\\"\\u003e\\u003cspan class=\\\"RefSource\\\"\\u003ehttps://doi.org/10.1261/rna.079623.123\\u003c/span\\u003e\\u003cspan address=\\\"10.1261/rna.079623.123\\\" targettype=\\\"DOI\\\" class=\\\"RefTarget\\\"\\u003e\\u003c/span\\u003e\\u003c/span\\u003e.\\u003c/span\\u003e\\u003c/li\\u003e \\u003cli\\u003e\\u003cspan\\u003eAird D, Ross MG, Chen WS, Danielsson M, Fennell T, Russ C, et al. Analyzing and minimizing PCR amplification bias in Illumina sequencing libraries. Genome Biol. 2011;12:R18. \\u003cspan class=\\\"ExternalRef\\\"\\u003e\\u003cspan class=\\\"RefSource\\\"\\u003ehttps://doi.org/10.1186/gb-2011-12-2-r18\\u003c/span\\u003e\\u003cspan address=\\\"10.1186/gb-2011-12-2-r18\\\" targettype=\\\"DOI\\\" class=\\\"RefTarget\\\"\\u003e\\u003c/span\\u003e\\u003c/span\\u003e.\\u003c/span\\u003e\\u003c/li\\u003e \\u003cli\\u003e\\u003cspan\\u003eUnterseher M, Schnittler M. Dilution-to-extinction cultivation of leaf-inhabiting endophytic fungi in beech (\\u003cem\\u003eFagus sylvatica\\u003c/em\\u003e L.) - Different cultivation techniques influence fungal biodiversity assessment. Mycol Res. 2009;113:645\\u0026ndash;54. \\u003cspan class=\\\"ExternalRef\\\"\\u003e\\u003cspan class=\\\"RefSource\\\"\\u003ehttps://doi.org/10.1016/j.mycres.2009.02.002\\u003c/span\\u003e\\u003cspan address=\\\"10.1016/j.mycres.2009.02.002\\\" targettype=\\\"DOI\\\" class=\\\"RefTarget\\\"\\u003e\\u003c/span\\u003e\\u003c/span\\u003e.\\u003c/span\\u003e\\u003c/li\\u003e \\u003cli\\u003e\\u003cspan\\u003eOrdonez A, Hussain U, Cambon MC, Golyshin PN, Downie J, McDonald JE. Evaluating agar-plating and dilution-to-extinction isolation methods for generating oak-associated microbial culture collections. ISME Commun. 2025;5:ycaf019. \\u003cspan class=\\\"ExternalRef\\\"\\u003e\\u003cspan class=\\\"RefSource\\\"\\u003ehttps://doi.org/10.1093/ismeco/ycaf019\\u003c/span\\u003e\\u003cspan address=\\\"10.1093/ismeco/ycaf019\\\" targettype=\\\"DOI\\\" class=\\\"RefTarget\\\"\\u003e\\u003c/span\\u003e\\u003c/span\\u003e.\\u003c/span\\u003e\\u003c/li\\u003e \\u003cli\\u003e\\u003cspan\\u003eSimonin M, Briand M, Chesneau G, Rochefort A, Marais C, Sarniguet A, et al. Seed microbiota revealed by a large-scale meta-analysis including 50 plant species. New Phytol. 2022;234:1448\\u0026ndash;63. \\u003cspan class=\\\"ExternalRef\\\"\\u003e\\u003cspan class=\\\"RefSource\\\"\\u003ehttps://doi.org/10.1111/nph.18037\\u003c/span\\u003e\\u003cspan address=\\\"10.1111/nph.18037\\\" targettype=\\\"DOI\\\" class=\\\"RefTarget\\\"\\u003e\\u003c/span\\u003e\\u003c/span\\u003e.\\u003c/span\\u003e\\u003c/li\\u003e \\u003cli\\u003e\\u003cspan\\u003eFranić I, Eschen R, Allan E, Hartmann M, Schneider S, Prospero S. Drivers of richness and community composition of fungal endophytes of tree seeds. FEMS Microbiol Ecol. 2020;96:fiaa166. \\u003cspan class=\\\"ExternalRef\\\"\\u003e\\u003cspan class=\\\"RefSource\\\"\\u003ehttps://doi.org/10.1093/femsec/fiaa166\\u003c/span\\u003e\\u003cspan address=\\\"10.1093/femsec/fiaa166\\\" targettype=\\\"DOI\\\" class=\\\"RefTarget\\\"\\u003e\\u003c/span\\u003e\\u003c/span\\u003e.\\u003c/span\\u003e\\u003c/li\\u003e \\u003cli\\u003e\\u003cspan\\u003eBoroń P, Grad B, Nawrot-Chorabik K, Kowalski T. The genetic relationship of \\u003cem\\u003eApiognomonia errabunda\\u003c/em\\u003e and related species. Mycologia. 2019;111:541\\u0026ndash;50. \\u003cspan class=\\\"ExternalRef\\\"\\u003e\\u003cspan class=\\\"RefSource\\\"\\u003ehttps://doi.org/10.2307/26756638\\u003c/span\\u003e\\u003cspan address=\\\"10.2307/26756638\\\" targettype=\\\"DOI\\\" class=\\\"RefTarget\\\"\\u003e\\u003c/span\\u003e\\u003c/span\\u003e.\\u003c/span\\u003e\\u003c/li\\u003e \\u003cli\\u003e\\u003cspan\\u003eChoi JW, Park E, Eo JK. An Unrecorded Genus \\u003cem\\u003eLambertella\\u003c/em\\u003e H\\u0026ouml;hn. (Rutstroemiaceae) and Its Unrecorded Species in Korea. Korean J Mycol. 2021;49:127\\u0026ndash;31. \\u003cspan class=\\\"ExternalRef\\\"\\u003e\\u003cspan class=\\\"RefSource\\\"\\u003ehttps://doi.org/10.4489/KJM.20210013\\u003c/span\\u003e\\u003cspan address=\\\"10.4489/KJM.20210013\\\" targettype=\\\"DOI\\\" class=\\\"RefTarget\\\"\\u003e\\u003c/span\\u003e\\u003c/span\\u003e.\\u003c/span\\u003e\\u003c/li\\u003e \\u003cli\\u003e\\u003cspan\\u003eRidout M, Newcombe G. \\u003cem\\u003eSydowia polyspora\\u003c/em\\u003e is both a foliar endophyte and a preemergent seed pathogen in \\u003cem\\u003ePinus ponderosa\\u003c/em\\u003e. Plant Dis. 2018;102:640\\u0026ndash;4. \\u003cspan class=\\\"ExternalRef\\\"\\u003e\\u003cspan class=\\\"RefSource\\\"\\u003ehttps://doi.org/10.1094/PDIS-07-17-1074-RE\\u003c/span\\u003e\\u003cspan address=\\\"10.1094/PDIS-07-17-1074-RE\\\" targettype=\\\"DOI\\\" class=\\\"RefTarget\\\"\\u003e\\u003c/span\\u003e\\u003c/span\\u003e.\\u003c/span\\u003e\\u003c/li\\u003e \\u003cli\\u003e\\u003cspan\\u003eGao Y, de Farias ARG, Jiang HB, Karunarathna SC, Xu JC, Tibpromma S, et al. Morphological and Phylogenetic Characterisations Reveal Four New Species in Leptosphaeriaceae (Pleosporales, Dothideomycetes). J Fungi. 2023;9. \\u003cspan class=\\\"ExternalRef\\\"\\u003e\\u003cspan class=\\\"RefSource\\\"\\u003ehttps://doi.org/10.3390/jof9060612\\u003c/span\\u003e\\u003cspan address=\\\"10.3390/jof9060612\\\" targettype=\\\"DOI\\\" class=\\\"RefTarget\\\"\\u003e\\u003c/span\\u003e\\u003c/span\\u003e.\\u003c/span\\u003e\\u003c/li\\u003e \\u003cli\\u003e\\u003cspan\\u003eAriyawansa HA, Phukhamsakda C, Thambugala KM, Bulgakov TS, Wanasinghe DN, Perera RH, et al. Revision and phylogeny of Leptosphaeriaceae. Fungal Divers. 2015;74:19\\u0026ndash;51. \\u003cspan class=\\\"ExternalRef\\\"\\u003e\\u003cspan class=\\\"RefSource\\\"\\u003ehttps://doi.org/10.1007/s13225-015-0349-2\\u003c/span\\u003e\\u003cspan address=\\\"10.1007/s13225-015-0349-2\\\" targettype=\\\"DOI\\\" class=\\\"RefTarget\\\"\\u003e\\u003c/span\\u003e\\u003c/span\\u003e.\\u003c/span\\u003e\\u003c/li\\u003e \\u003cli\\u003e\\u003cspan\\u003eSchr\\u0026ouml;der T, Kehr R, H\\u0026uuml;ttermann A. First report of the seed-pathogen \\u003cem\\u003eGeniculodendron pyriforme\\u003c/em\\u003e, the imperfect state of the ascomycete \\u003cem\\u003eCaloscypha fulgens\\u003c/em\\u003e, on imported conifer seeds in Germany. Pathol. 2002;32:225\\u0026ndash;30. \\u003cspan class=\\\"ExternalRef\\\"\\u003e\\u003cspan class=\\\"RefSource\\\"\\u003ehttps://doi.org/10.1046/j.1439-0329.2002.00288.x\\u003c/span\\u003e\\u003cspan address=\\\"10.1046/j.1439-0329.2002.00288.x\\\" targettype=\\\"DOI\\\" class=\\\"RefTarget\\\"\\u003e\\u003c/span\\u003e\\u003c/span\\u003e.\\u003c/span\\u003e\\u003c/li\\u003e \\u003cli\\u003e\\u003cspan\\u003eSalas-Lizana R, Oono R. Double-digest RADseq loci using standard Illumina indexes improve deep and shallow phylogenetic resolution of \\u003cem\\u003eLophodermium\\u003c/em\\u003e, a widespread fungal endophyte of pine needles. Ecol Evol. 2018;8:6638\\u0026ndash;51. \\u003cspan class=\\\"ExternalRef\\\"\\u003e\\u003cspan class=\\\"RefSource\\\"\\u003ehttps://doi.org/10.1002/ece3.4147\\u003c/span\\u003e\\u003cspan address=\\\"10.1002/ece3.4147\\\" targettype=\\\"DOI\\\" class=\\\"RefTarget\\\"\\u003e\\u003c/span\\u003e\\u003c/span\\u003e.\\u003c/span\\u003e\\u003c/li\\u003e \\u003cli\\u003e\\u003cspan\\u003ePhukhamsakda C, Nilsson RH, Bhunjun CS, de Farias ARG, Sun YR, Wijesinghe SN, et al. The numbers of fungi: contributions from traditional taxonomic studies and challenges of metabarcoding. Fungal Divers. 2022;114:327\\u0026ndash;86. \\u003cspan class=\\\"ExternalRef\\\"\\u003e\\u003cspan class=\\\"RefSource\\\"\\u003ehttps://doi.org/10.1007/s13225-022-00502-3\\u003c/span\\u003e\\u003cspan address=\\\"10.1007/s13225-022-00502-3\\\" targettype=\\\"DOI\\\" class=\\\"RefTarget\\\"\\u003e\\u003c/span\\u003e\\u003c/span\\u003e.\\u003c/span\\u003e\\u003c/li\\u003e \\u003cli\\u003e\\u003cspan\\u003eTalhinhas P, Baroncelli R. \\u003cem\\u003eColletotrichum\\u003c/em\\u003e species and complexes: geographic distribution, host range and conservation status. Fungal Divers. 2021;110:109\\u0026ndash;98. \\u003cspan class=\\\"ExternalRef\\\"\\u003e\\u003cspan class=\\\"RefSource\\\"\\u003ehttps://doi.org/10.1007/s13225-021-00491-9\\u003c/span\\u003e\\u003cspan address=\\\"10.1007/s13225-021-00491-9\\\" targettype=\\\"DOI\\\" class=\\\"RefTarget\\\"\\u003e\\u003c/span\\u003e\\u003c/span\\u003e.\\u003c/span\\u003e\\u003c/li\\u003e \\u003cli\\u003e\\u003cspan\\u003eLiao C, Doilom M, Jeewon R, Hyde KD, Manawasinghe IS, Chethana KWT, et al. Challenges and update on fungal endophytes: classification, definition, diversity, ecology, evolution and functions. Fungal Divers. 2025;131:301\\u0026ndash;67. \\u003cspan class=\\\"ExternalRef\\\"\\u003e\\u003cspan class=\\\"RefSource\\\"\\u003ehttps://doi.org/10.1007/s13225-025-00550-5\\u003c/span\\u003e\\u003cspan address=\\\"10.1007/s13225-025-00550-5\\\" targettype=\\\"DOI\\\" class=\\\"RefTarget\\\"\\u003e\\u003c/span\\u003e\\u003c/span\\u003e.\\u003c/span\\u003e\\u003c/li\\u003e \\u003cli\\u003e\\u003cspan\\u003eCram MM, Kasten Dumroese R. Mycorrhizae in Forest Tree Nurseries. In: Cram MM, Frank MS, Mallams KM, editors. Forest Nursery Pests. Washington, DC, USA: USDA, Forest Service; 2012. pp. 20\\u0026ndash;3.\\u003c/span\\u003e\\u003c/li\\u003e \\u003cli\\u003e\\u003cspan\\u003eNocker A, Sossa KE, Camper AK. Molecular monitoring of disinfection efficacy using propidium monoazide in combination with quantitative PCR. J Microbiol Methods. 2007;70:252\\u0026ndash;60. \\u003cspan class=\\\"ExternalRef\\\"\\u003e\\u003cspan class=\\\"RefSource\\\"\\u003ehttps://doi.org/10.1016/j.mimet.2007.04.014\\u003c/span\\u003e\\u003cspan address=\\\"10.1016/j.mimet.2007.04.014\\\" targettype=\\\"DOI\\\" class=\\\"RefTarget\\\"\\u003e\\u003c/span\\u003e\\u003c/span\\u003e.\\u003c/span\\u003e\\u003c/li\\u003e \\u003cli\\u003e\\u003cspan\\u003eNocker A, Camper AK. Selective removal of DNA from dead cells of mixed bacterial communities by use of ethidium monoazide. Appl Environ Microbiol. 2006;72:1997\\u0026ndash;2004. \\u003cspan class=\\\"ExternalRef\\\"\\u003e\\u003cspan class=\\\"RefSource\\\"\\u003ehttps://doi.org/10.1128/AEM.72.3.1997-2004.2006\\u003c/span\\u003e\\u003cspan address=\\\"10.1128/AEM.72.3.1997-2004.2006\\\" targettype=\\\"DOI\\\" class=\\\"RefTarget\\\"\\u003e\\u003c/span\\u003e\\u003c/span\\u003e.\\u003c/span\\u003e\\u003c/li\\u003e\\u003c/ol\\u003e\"}],\"fulltextSource\":\"\",\"fullText\":\"\",\"funders\":[],\"hasAdminPriorityOnWorkflow\":false,\"hasManuscriptDocX\":true,\"hasOptedInToPreprint\":true,\"hasPassedJournalQc\":\"\",\"hasAnyPriority\":false,\"hideJournal\":false,\"highlight\":\"\",\"institution\":\"\",\"isAcceptedByJournal\":false,\"isAuthorSuppliedPdf\":false,\"isDeskRejected\":\"\",\"isHiddenFromSearch\":false,\"isInQc\":false,\"isInWorkflow\":false,\"isPdf\":false,\"isPdfUpToDate\":true,\"isWithdrawnOrRetracted\":false,\"journal\":{\"display\":true,\"email\":\"info@researchsquare.com\",\"identity\":\"environmental-microbiome\",\"isNatureJournal\":false,\"hasQc\":true,\"allowDirectSubmit\":false,\"externalIdentity\":\"sigs\",\"sideBox\":\"Learn more about [Environmental Microbiome](https://environmentalmicrobiome.biomedcentral.com)\",\"snPcode\":\"40793\",\"submissionUrl\":\"https://submission.nature.com/new-submission/40793/3\",\"title\":\"Environmental Microbiome\",\"twitterHandle\":\"@bmc\",\"acdcEnabled\":true,\"dfaEnabled\":true,\"editorialSystem\":\"em\",\"reportingPortfolio\":\"BMC/SO AJ\",\"inReviewEnabled\":true,\"inReviewRevisionsEnabled\":true},\"keywords\":\"microbiome, high-throughput sequencing, viable community, metabolic activity, relic DNA, phytosanitary risk\",\"lastPublishedDoi\":\"10.21203/rs.3.rs-8095745/v1\",\"lastPublishedDoiUrl\":\"https://doi.org/10.21203/rs.3.rs-8095745/v1\",\"license\":{\"name\":\"CC BY 4.0\",\"url\":\"https://creativecommons.org/licenses/by/4.0/\"},\"manuscriptAbstract\":\"\\u003ch2\\u003eBackground\\u003c/h2\\u003e \\u003cp\\u003eTree seeds harbor diverse fungal communities, including both pathogens and mutualists, that can influence plant health. These communities comprise living, metabolically active organisms as well as dormant or dead cells. Because only active fungi interact with their hosts, distinguishing active from inactive taxa is crucial, especially for environmental and phytosanitary monitoring. Traditional culturing methods capture living fungi but account for only a small fraction of the total fungal diversity. Currently, these methods are increasingly replaced by high-throughput DNA metabarcoding, which detects a broader range of taxa. However, DNA persists after cell death and occurs in dormant cells, preventing distinction between active and inactive fungi. In contrast, RNA metabarcoding detects metabolically active organisms and may better reflect living fungal communities than the other two methods, though its use in assessing plant-associated fungi remains underexplored. We used culturing, DNA-, and RNA-based metabarcoding to compare fungal communities associated with seeds of three key European tree species (\\u003cem\\u003eFagus sylvatica\\u003c/em\\u003e, \\u003cem\\u003eAbies alba\\u003c/em\\u003e, \\u003cem\\u003ePinus sylvestris\\u003c/em\\u003e).\\u003c/p\\u003e\\u003ch2\\u003eResults\\u003c/h2\\u003e \\u003cp\\u003eDNA and RNA metabarcoding detected largely distinct, non-overlapping fungal communities, with differences primarily driven by rare active taxa in the RNA dataset. Several cultured genera\\u0026mdash;likely representing abundant and metabolically active taxa\\u0026mdash;were shared between both metabarcoding approaches.\\u003c/p\\u003e\\u003ch2\\u003eConclusions\\u003c/h2\\u003e \\u003cp\\u003eThese results highlight the complementary nature of the three methods for characterising seed-associated fungi. Combining culturing, DNA- and RNA-based metabarcoding may provide the most comprehensive assessment of fungal diversity, while RNA metabarcoding alone offers a promising opportunity to identify the active members of fungal communities for improved environmental and phytosanitary monitoring.\\u003c/p\\u003e\",\"manuscriptTitle\":\"DNA and RNA metabarcoding reveal distinct seed-borne mycobiota\",\"msid\":\"\",\"msnumber\":\"\",\"nonDraftVersions\":[{\"code\":1,\"date\":\"2026-02-17 09:54:46\",\"doi\":\"10.21203/rs.3.rs-8095745/v1\",\"editorialEvents\":[{\"type\":\"communityComments\",\"content\":0},{\"type\":\"decision\",\"content\":\"Revision requested\",\"date\":\"2026-03-27T04:17:00+00:00\",\"index\":\"\",\"fulltext\":\"\"},{\"type\":\"editorInvitedReview\",\"content\":\"\",\"date\":\"2026-03-25T10:40:11+00:00\",\"index\":\"hide\",\"fulltext\":\"\"},{\"type\":\"editorInvitedReview\",\"content\":\"\",\"date\":\"2026-03-17T22:23:18+00:00\",\"index\":\"hide\",\"fulltext\":\"\"},{\"type\":\"editorInvitedReview\",\"content\":\"\",\"date\":\"2026-03-14T15:50:33+00:00\",\"index\":\"hide\",\"fulltext\":\"\"},{\"type\":\"editorInvitedReview\",\"content\":\"\",\"date\":\"2026-03-14T15:41:20+00:00\",\"index\":\"hide\",\"fulltext\":\"\"},{\"type\":\"reviewerAgreed\",\"content\":\"139325054428239318688772311727467733990\",\"date\":\"2026-02-17T03:40:58+00:00\",\"index\":\"hide\",\"fulltext\":\"\"},{\"type\":\"reviewerAgreed\",\"content\":\"317470364710376684768060229454784206253\",\"date\":\"2026-02-12T15:43:34+00:00\",\"index\":\"hide\",\"fulltext\":\"\"},{\"type\":\"reviewerAgreed\",\"content\":\"62352256862363872486189041166491381722\",\"date\":\"2026-02-12T15:12:39+00:00\",\"index\":\"hide\",\"fulltext\":\"\"},{\"type\":\"reviewerAgreed\",\"content\":\"236384891441576557720689319848852214520\",\"date\":\"2026-02-12T04:52:50+00:00\",\"index\":\"hide\",\"fulltext\":\"\"},{\"type\":\"reviewersInvited\",\"content\":\"\",\"date\":\"2026-02-12T04:49:11+00:00\",\"index\":\"\",\"fulltext\":\"\"},{\"type\":\"editorAssigned\",\"content\":\"\",\"date\":\"2026-01-26T12:35:30+00:00\",\"index\":\"\",\"fulltext\":\"\"},{\"type\":\"checksComplete\",\"content\":\"\",\"date\":\"2026-01-23T09:17:24+00:00\",\"index\":\"\",\"fulltext\":\"\"},{\"type\":\"submitted\",\"content\":\"Environmental Microbiome\",\"date\":\"2026-01-19T13:27:52+00:00\",\"index\":\"\",\"fulltext\":\"\"}],\"status\":\"published\",\"journal\":{\"display\":true,\"email\":\"info@researchsquare.com\",\"identity\":\"environmental-microbiome\",\"isNatureJournal\":false,\"hasQc\":true,\"allowDirectSubmit\":false,\"externalIdentity\":\"sigs\",\"sideBox\":\"Learn more about [Environmental Microbiome](https://environmentalmicrobiome.biomedcentral.com)\",\"snPcode\":\"40793\",\"submissionUrl\":\"https://submission.nature.com/new-submission/40793/3\",\"title\":\"Environmental Microbiome\",\"twitterHandle\":\"@bmc\",\"acdcEnabled\":true,\"dfaEnabled\":true,\"editorialSystem\":\"em\",\"reportingPortfolio\":\"BMC/SO AJ\",\"inReviewEnabled\":true,\"inReviewRevisionsEnabled\":true}}],\"origin\":\"\",\"ownerIdentity\":\"8323dbf6-d8a9-420c-a369-d8a76e558e0f\",\"owner\":[],\"postedDate\":\"February 17th, 2026\",\"published\":true,\"recentEditorialEvents\":[],\"rejectedJournal\":[],\"revision\":\"\",\"amendment\":\"\",\"status\":\"under-review\",\"subjectAreas\":[],\"tags\":[],\"updatedAt\":\"2026-04-27T04:24:00+00:00\",\"versionOfRecord\":[],\"versionCreatedAt\":\"2026-02-17 09:54:46\",\"video\":\"\",\"vorDoi\":\"\",\"vorDoiUrl\":\"\",\"workflowStages\":[]},\"version\":\"v1\",\"identity\":\"rs-8095745\",\"journalConfig\":\"researchsquare\"},\"__N_SSP\":true},\"page\":\"/article/[identity]/[[...version]]\",\"query\":{\"redirect\":\"/article/rs-8095745\",\"identity\":\"rs-8095745\",\"version\":[\"v1\"]},\"buildId\":\"XKTyCvWXoU3ODBz1xrDgd\",\"isFallback\":false,\"isExperimentalCompile\":false,\"dynamicIds\":[84888],\"gssp\":true,\"scriptLoader\":[]}","source_license":"CC-BY-4.0","license_restricted":false}