- Open Access
Landscape scale ecology of Tetracladium spp. fungal root endophytes
Environmental Microbiome volume 17, Article number: 40 (2022)
The genus Tetracladium De Wild. (Ascomycota) has been traditionally regarded as a group of Ingoldian fungi or aquatic hyphomycetes—a polyphyletic group of phylogenetically diverse fungi which grow on decaying leaves and plant litter in streams. Recent sequencing evidence has shown that Tetracladium spp. may also exist as root endophytes in terrestrial environments, and furthermore may have beneficial effects on the health and growth of their host. However, the diversity of Tetracladium spp. communities in terrestrial systems and the factors which shape their distribution are largely unknown.
Using a fungal community internal transcribed spacer amplicon dataset from 37 UK Brassica napus fields we found that soils contained diverse Tetracladium spp., most of which represent previously uncharacterised clades. The two most abundant operational taxonomic units (OTUs), related to previously described aquatic T. furcatum and T. maxilliforme, were enriched in roots relative to bulk and rhizosphere soil. For both taxa, relative abundance in roots, but not rhizosphere or bulk soil was correlated with B. napus yield. The relative abundance of T. furcatum and T. maxilliforme OTUs across compartments showed very similar responses with respect to agricultural management practices and soil characteristics. The factors shaping the relative abundance of OTUs homologous to T. furcatum and T. maxilliforme OTUs in roots were assessed using linear regression and structural equation modelling. Relative abundance of T. maxilliforme and T. furcatum in roots increased with pH, concentrations of phosphorus, and increased rotation frequency of oilseed rape. It decreased with increased soil water content, concentrations of extractable phosphorus, chromium, and iron.
The genus Tetracladium as a root colonising endophyte is a diverse and widely distributed part of the oilseed rape microbiome that positively correlates to crop yield. The main drivers of its community composition are crop management practices and soil nutrients.
Aquatic hyphomycetes or Ingoldian fungi are important decomposers in freshwater ecosystems . Spores of these fungi were first described from running freshwater streams in the 1940s , with species classified according to morphology—primarily sigmoid or tetraradiate . Sexual reproduction of these fungi has never been observed, and the members of this group don’t share common morphological or ecological characteristics . However, the common conidial shape suggests convergent evolution, and may contribute to spore dispersal via improved anchoring to the substrate or higher buoyancy for better aquatic dispersal . Recently, use of high-throughput sequencing has revealed the presence of aquatic hyphomycetes in fungal communities inhabiting soil and plants, although the ecological importance of these fungi in terrestrial habitats is unknown .
The genus Tetracladium is a group of common aquatic hyphomycete that was first described by de Wildeman in 1893 , and sits within the ascomycete class Leotiomycetes in the Han Clade 9/Stamnaria lineage/Vandijckellaceae clade as incertae sedis . The genus name was coined in response to its distinct ~ 60 × 100 µm tetra formatted conidiospores which have a central axis with three radiating branches . Since the initial description of Tetracladium, the genus has been found to be ubiquitous in aquatic environments [9,10,11,12,13,14]. The first terrestrial observations of Tetracladium were from forest litter [11, 15, 16], with fungal spores detected in the water film covering fallen leaves . However, most reports of Tetracladium in terrestrial environments came after the turn of the century as DNA sequencing techniques became more easily accessible. Most of this data comes from environmental metabarcoding studies, and there are only a few instances of Tetracladium ssp. being isolated in pure cultures. It has been hypothesised that there may have been under-reporting of Tetracladium in terrestrial habitats before the 2000s because of the strange nature of finding an aquatic organism in a terrestrial environment . It is unclear whether the species described based on spore morphology from aquatic habitats and the DNA sequences identified from terrestrial environmental samples belong the same organisms. However, nuclear ribosomal internal transcribed spacer (ITS) amplicon analysis has shown no sequence-based differences between aquatic and terrestrial strains of a number of species, indicating that some species may have diverse ecological functions .
One of the first observations of plant endophytic Tetracladium sp. came from riparian plant roots . These fungi do not appear to show host or habitat specificity as plant endophytes, and they have been found in roots of monocot species within Asparagales [18,19,20], Liliales , and Poales [22,23,24], as well as dicot species within Ericales , Brassicales [26, 27] and Vitales . Furthermore, they have been found associated with Equisetaceae [29, 30] and Bryophytes [31,32,33,34]. Tetracladium has most frequently been described in metabarcoding studies of soil from disturbed agricultural and grassland habitats [35,36,37,38]. There are numerous reports of the genus from the Arctics [39,40,41,42] and they have also been found in unvegetated habitats including glaciers and bedrock [43,44,45].
The dual ecology of Tetracladium sp., and particularly their importance as plant endophytes, are still debated. Anderson et al.  investigated an aquatic T. marchalianum population over time and space and found that the fungus maintains a high genotypic diversity throughout the year. They suggested that this could be attributed to their alternative lifestyles as terrestrial plant endophytes . It was proposed by Selosse et al.  that the terrestrial occurrence of aquatic hyphomycetes, and more specifically their endophytic nature, is attributed to the fungi precolonising plant tissues and building biomass, so in the event of abscission, they are already occupying the niche, and ready to decompose plant litter which reaches freshwater. It was also suggested that the tetraradiate spore morphology could aid them in becoming airborne . Based on this theory Tetracladium sp. should be the most common in aerial plant tissues, however there is currently no evidence to suggest that this is true.
There is conflicting evidence on whether Tetracladium infection provides benefits to the host. Glasshouse experiments have shown that inoculation with Tetracladium sp. can have beneficial effects on plant growth , while other studies have shown no effects . Importantly, Hilton et al.  found that Tetracladium sp. had a co-exclusion relationship with root pathogenic fungi, and relative abundance in roots was positively associated with crop yield. Furthermore, T. marchalianum showed an antagonistic effect against bacterial plant pathogens including Erwinia chrysanthemi and Xanthomonas phaseoli, although other Tetracladium sp. showed no such effects .
To date there have been no systematic studies that have investigated the diversity and distribution of Tetracladium in terrestrial habitats, and as a result the factors which shape Tetracladium spp. communities. Thus, the extent to which they interact with plants are unclear. In the current study we build on our earlier work  which characterised root fungal communities of Brassica napus across 37 UK fields to investigate (1) the diversity of Tetracladium spp. in soil and roots at the landscape scale (2) the extent to which different Tetracladium spp. are selectively recruited into roots and rhizosphere soil from bulk soil (3) the relationships between root, rhizosphere and soil populations of different Tetracladium spp. and crop yield and (4) to determine the importance of, and interactions between, soil nutrients, climate, soil physical properties, and crop management practices as drivers for the colonisation of roots by Tetracladium spp.
Sample collection and analyses
Root, rhizosphere soil and bulk soil samples were collected in March 2015 from 37 oilseed rape (B. napus) fields from 25 commercial farms in the UK. Five composite samples were taken from each field site. Loosely adhering soil was removed from the roots only leaving 2 mm rhizosphere soil. Roots with closely adhering soil were washed four times in sterile distilled water to release the rhizosphere soil which was then centrifuged (3250×g for 10 min) to leave a pellet of rhizosphere soil. Washed roots that were less than 2 mm diameter were cut into 5 mm pieces to obtain root samples. Bulk soil samples were sieved (7 mm followed by 2 mm sieve) then approximately 6 g was washed in sterile distilled water using the same washing method as the rhizosphere soil samples to ensure comparability. Following DNA extraction (PowerSoil-htp™ 96 Well Soil DNA Isolation Kit (MoBio Laboratories, Carlsbad, CA, USA)), the fungal community was amplified using internal transcribed spacer (ITS) primers fITS7-ITS4 . Sequencing was performed with Illumina MiSeq technology (Illumina MiSeq Reagent Kit v3), and taxonomy assigned using Quantitative Insights into Microbial Ecology (QIIME 1.8)  with the UNITE ITS database . Sequences were clustered to operational taxonomic units (OTU)  at 97% minimum identity threshold, and those OTUs assigned as Tetracladium spp. were selected for use in the current study. Metadata collected from each field included soil physico-chemical parameters (including C, N, P, micronutrient, pH, and soil type), climatic data, crop variety, rotation sequence and grain yield at the subsequent harvest. Full methodological details for sample preparation, DNA extraction, sequencing and bioinformatic analysis can be found in Hilton et al. . Information about farm locations, Tetracladium spp. OTUs and metadata can be found in Additional file 4 A and B.
To obtain more detailed information about the phylogenetic relatedness of recovered Tetracladium spp. sequences, the most closely related sequences to these OTUs were downloaded, including two representative ITS2 sequences from all described species. Sequences were accessed from the NCBI GenBank, and were aligned with our Tetracladium OTU sequences using MAFFT v.7 (e-ins-I algorithm) . The multiple sequence alignment can be found in Additional file 7. Maximum likelihood analyses were performed with RAxML on the CIPRES Science Gateway to build a phylogenetic tree using the default setting with 1000 bootstrap replicates [54, 55].
Observed species counts were used to generate richness plots to study OTU abundance differences in the three sampled compartments (bulk soil, rhizosphere, and root) using vegan (version 2.6–2) in R (version 4.12) . Rarefaction curves were constructed to assess the extent to which Tetracladium spp. richness was captured at the sequencing depth used. A heatmap was created to investigate the distribution of Tetracladium spp. OTUs across the sampling locations and between the soil, rhizosphere, and root compartments. Heatmap and rarefaction analyses were carried out using phyloseq (version 1.36.0) . Significance of differences in taxa richness between compartments and differences in OTU relative abundance between crop genotypes and previous cultivated crops were tested using the Kruskal–Wallis rank sum test. P values were corrected for multiple comparisons with a Dunn’s test using the false discovery rate with the Benjamini–Hochberg method. Linear regression was used to correlate relative abundance to yield and rotation. Zero values were introduced to accommodate for fields that never had oilseed rape sown before, therefore rotation length values are reciprocal. Ternary plots were created using ggtern (version 3.3.5)  to understand the compartment preference of the OTUs.
Out of all the OTUs found that resembled the genus Tetracladium we chose the two most abundant OTUs, which were also substantively enriched within the root compartment, for detailed ecological analyses. Data was normalised using modified Z-scores. To test for drivers of relative abundance in the roots, we created a piecewise structural equation model (PSEM). First, a correlogram was created to better understand relationships between metadata and OTU relative abundance (Additional file 1) using base R functions . Then, individual linear mixed-effect models were fitted with sampling location and soil compartment as random variables. Significant soil nutrient factors, relative to OTU abundance, taken from the correlogram output and the initial fitted model with all soil nutrients (modM) (Additional files 1 and 2), soil structure and climate were used as composite fixed variables. These composite variables were created so complicated constructs can be processed as simpler blocks that are easier to present and discuss . Variable reduction was done to unmask significant relationships that may be missed if too many factors were included based on the LMM fitted with all variables. Individual models were fitted using the R package lme4 . Fixed variables were reduced via assessing best model fit using the performance package . Finally, path models were fitted with the piecewiseSEM package , based on the findings of the individual models.
Across the 37 fields we found twelve OTUs that represented the genus Tetracladium (Fig. 1A). The closest match of the Tetracladium sequences to UNITE DOI is shown in Additional File 4 E. The accession numbers of the Tetracladium OTU sequences are shown in Additional File 4 F. Sequencing depth was sufficient to capture richness for most samples (Additional file 3). Higher abundance Tetracladium sp. OTUs (OTUs 19, 1088, 168, 6663, 359, 4156, and 5882) were found in all sampled fields, while lower abundance OTUs (813, 1055, 3952, 6656, and 10312) were found sporadically across sampling locations (Fig. 2, Additional file 4 C). The highest abundance OTUs, OTU 19 and 1088, grouped with T. maxilliforme and T. furcatum respectively, both of which have been described from water (Fig. 1B). OTU 813 clustered closely with, but was not identical to the species T. elipsoideum, which has been described from Arctic soil  and the closest uncultured environmental sequences to this OTU (MF181805.1 and MK627297.1) also originate from soil. OTUs 6663, 5882, 4156, and 168 clustered with uncultured Tetracladium sp. sequences in a clade close to T. ellipsoideum, T. psychrophilum and T. globosum. The closet uncultured environmental sequences for OTU 6663 (LR863329.1 and KX192428.1), OTU 5882 (MN660389.1), OTU 4156 (KF296960.1 and JX029127.1), and OTU 168 (GU055746.1) all originate from terrestrial samples. Environmental sequence MN660932.1—found in water—was a close match to OTU 168. OTUs 1055, 3952, 6656, 10312, and 359 formed a distinct clade with uncultured Tetracladium sp. sequences (MW050202.1, MH451254.1, KX193670.1, MH451294.1, MG756632.1, KM246272.1, MK246209.1, LR876862.1, and MK627349.1) found on land, likely representing uncharacterised species.
There was a significant difference between observed OTU richness of the three compartments (Fig. 3A), which was the lowest in the roots while the bulk soil and rhizosphere showed the same OTU richness (rhizosphere-bulk soil P = < 0.001, root-bulk soil P = < 0.001, root-rhizosphere P = < 0.001). Five of the OTUs (OTUs 19, 5882, 1055, 813, and 1088) were the most abundant in the roots, while OTU 6656 was found in greater abundance in the bulk soil, and OTU 10312 was only found in the rhizosphere. The rest of the OTUs did not show a specific preference for compartment (Fig. 3B). The most abundant Tetracladium OTUs (19 and 1088) both showed a strong preference for the roots. The mean relative abundance of OTU 19 was five times higher in the roots than in the bulk soil and it was over a thousand times higher in the case of OTU 1088 (Additional file 4 D). Based on these results we conclude that OTUs 19 and 1088 can be categorised as putative root colonising fungi and further analyses focussed on these two abundant and widely distributed OTUs. Other root associated OTUs (OTUs 813, 1055, and 5882) could only be found in high abundance at sporadic sampling locations and were present in low abundance or were absent at most sites (Fig. 2).
Correlating the metadata to OTU relative abundance
Relative abundance of OTUs 19 and 1088 in the root compartment had a significant positive linear correlation with oilseed rape (OSR, Brassica napus) yield (OTU 19 − R2 = 0.078, P = < 0.001, OTU 1088 − R2 = 0.067, P = < 0.001). The highest relative abundance of OTUs 19 and 1088 in the roots was associated with a yield increase of up to 25% relative to samples with the lowest relative abundance. No such relationship was seen in the bulk soil or rhizosphere (Fig. 3C).
We further investigated the relationships between previous crop cultivated at the site, OSR variety and rotation and the relative abundance of OTUs 19 and 1088. We found that relative abundance of OTU 19 and 1088 had a significant correlation with previous cultivated crop (R2 = 0.021 P = 0.013 and R2 = 0.021, P = 0.012 respectively). In roots, relative abundance of both OTU 19 and 1088 was significantly higher (P < 0.001) when rye rather than barley or wheat was the proceeding crop (Fig. 4A). In the rhizosphere and soil relative abundance of OTU 19 showed the same patterns in relation to previous crops as for the roots, but in addition relative abundance following fallow was significantly higher than barley and lower than wheat (Fig. 4A, B, Additional file 5).
OSR variety had a significant effect on relative abundance of OTU 19 (R2 = 0.615, P = < 0.001) and 1088 (R2 = 0.614, P = < 0.001) in roots, with the two OTUs showing very similar distribution patterns. Seven varieties (Quartz, Rocca, Incentive, Compass, Catena, Camelot and Cabernet) had very low mean relative abundances of these OTUs, while Nikita had substantively higher relative abundance than the other varieties (Fig. 5A, B, Additional file 6). Rotation also had a significant linear correlation with the relative abundance of OTUs 19 and 1088 in the roots (R2 = 0.032, P = 0.023 for OTU 19, R2 = 0.059, P = 0.002 for OTU 1088) (Fig. 5B), which increased as time since previous OSR crop increased.
We found a strong correlation between the relative abundances of the root specific Tetracladium OTUs (R2 = 0.887, P = < 0.001) so for further analyses we used the mean of the combined relative abundance of the two OTUs of each sample (Fig. 6). In the final PSEM, we correlated relative abundance of the combined OTUs to soil structure (soil moisture content, bulk density, and sand content), significant nutrients (Olsen P, iron, chromium, phosphorus, and manganese), pH, rotation, and climate (annual rainfall and minimum annual temperature) using sampled field as a random variable and used simple linear regression to correlate nutrients to soil structure and then pH to nutrients. Nutrient variables were chosen based on the results of the correlation matrix and the assessment of the model fit of the initial linear mixed-effect model modM (Additional file 1).
As a result, the PSEM showed the strongest significant direct positive effect of nutrients (R2 conditional = 0.48, standard estimate = 0.48, P = 0.001) followed by pH (R2 conditional = 0.07, standard estimate = 0.27, P = 0.002) and rotation (R2 conditional = 0.35, standard estimate = 0.17, P = 0.043) on the mean combined relative abundance of OTUs 19 and 1088 (Fig. 6). The effect of climate and soil structure was not significant in the path analyses. Soil structure and pH had significant positive correlations with nutrients (R2 = 0.12, standard estimate = 0.34, P = < 0.001 and R2 conditional = 0.07, standard estimate = 0.27, P = 0.005) thus indirectly affecting OTU relative abundance. To expand the composite variables of the PSEM, the mixed-effect linear models showed significant correlation between the combined mean relative abundance of OTUs 19 and 1088 (Fig. 6) with Olsen P (P = < 0.001), phosphorus (P = 0.003) and iron (P = 0.008). We also found a significant negative correlation between the combined mean relative abundance of the OTUs and soil water content (P = 0.045), while none of the climate variables were significant (Fig. 6).
Here we present the first systematic study of the landscape scale diversity and distribution of Tetracladium spp. within terrestrial systems and identify key factors controlling their occurrence as root endophytes. Tetracladium spp. were widely distributed, occurring in soil, rhizosphere, and roots in the 37 sampled sites. A total of 12 Tetracladium spp. OTUs were detected, and a subset of OTUs was specifically enriched in oilseed rape roots, including OTUs 19 and 1088 defined at 97% as T. maxilliforme and T. furcatum. Eight of the OTUs belonged to clades for which only environment sequences, largely from terrestrial habitats, have been recovered. There was a significant relationship between relative abundance of T. maxilliforme and T. furcatum within oilseed rape roots and crop genotype, previous cultivated crop, and oilseed rape rotation period. Linear mixed effects modelling and piecewise structural equation modelling showed that the most important environmental drivers of the relative abundance of Tetracladium spp. within plant roots were pH and select nutrients, including total phosphorus, extractable (Olsen) phosphorus, and iron.
Diversity and distribution
Tetracladium is generally considered to be a freshwater fungi. However, members of the genus have also been detected as root endophytes of terrestrial plants. In the sampled fields there was lower species richness in the roots than in the soil and rhizosphere which is a common feature of endophytes, and suggests selective recruitment into the mycobiome . Based on ITS sequences, the two most abundant OTUs (OTU 19 and 1088) clustered with T. maxilliforme and T. furcatum respectively and showed a strong root preference (Fig. 3B). T. maxilliforme and T. furcatum have been found several times in agriculture including the roots of crop plants [22, 25], though they are both traditionally regarded as aquatic organisms [14, 65,66,67,68,69]. OTU 813 clustered with the terrestrial T. ellipsoideum and it showed a strong root preference (Fig. 3B). Excluding OTUs 19, 1088, and 813, the closest GenBank sequence matches for all OTUs were from terrestrial habitats. Compartment preference did not have a correlation with taxonomic position as root-preferring OTUs were found across the main Tetracladium spp. clades. This suggests that the different species in the genus may not have one and the same defined lifestyle. Overall, our findings suggest that there is considerable diversity within terrestrial Tetracladium and that while some taxa may inhabit both aquatic and terrestrial habitats, others may inhabit terrestrial systems, with some occurring as plant endophytes.
Root endophytes have been shown to increase host resistance to environmental factors such as drought, heat and saline stress [70, 71], and can also increase plant health by inducing increased resistance via priming of the natural immune system to pathogens  or through competition with pathogens for nutrients and spatial niche exclusion . Plant growth and yield can also be increased by endophytes via direct nutrient transfer from fungus to plant . Furthermore, root colonising endophytes may shape plant community diversity and distribution .
Although Tetracladium species have previously been found in roots [19, 22, 25, 26], little is known about the process(es) by which Tetracladium spp. colonise roots, or the significance of infection for plant health. Sati and Arya  found that inoculation with T. nainitalense had no significant effect on the growth of Hibiscus esculentus or Solanum melongena following inoculation with T. nainitalense, although there was no evidence that the inoculant colonised plant tissues. In our previous work  we found a positive correlation between Tetracladium OTU relative abundance and oilseed rape yield on a landscape scale and in the current study we build on this to show a 25% yield increase from the lowest to the highest OTU relative abundance with both OTUs 19 and 1088. There is a clear need to understand the root infection process by Tetracladium spp. and to quantify benefits for plant health under controlled conditions, so that the significance of Tetracladium spp. root endophytes and their potential to act as beneficial symbionts can be established.
Drivers of relative abundance
The root associated OTUs homologous to T. maxilliforme and T. furcatum OTUs showed strong co-assembly patterns in oilseed rape roots at a landscape scale. These OTUs showed the same interactions with host genotype, crop management, and environmental factors. There is evidence of minimal competition between root colonising endophytes explaining their high diversity within a single host . This enables similar strains or species of fungi to colonise the same plant in high abundance. The data presented here could simply indicate similar adaptation of the two species without any ecological interaction.
According to our fitted models, the main drivers of Tetracladium sp. relative abundance in oilseed rape roots were soil nutrient content, crop rotation and pH. Relative abundance had a positive correlation with soil phosphorus and a negative correlation with iron content and Olsen P. Phosphorus and iron availability may limit plant and microbial growth in soil. For example, dark septate endophytes, which like Tetracladium spp. belong to the Helotiales, and occur as root endophytes, have been found to have iron phosphate solubilisation properties . Sati and Pant found phosphate solubilisation in T. setigerum isolated from riparian roots in agar and broth media . Moreover, it was shown that mineral fertilisers increased the relative abundance of Tetracladium sp. indicating that mineral fertiliser treatment might promote this plant-fungal symbiosis [23, 79]. In our study, total phosphorus content of the soil had a positive relationship with Tetracladium spp. OTU relative abundance in roots. In contrast, extractable phosphorus or Olsen P had a negative relationship with Tetracladium spp. OTU relative abundance in roots. This could suggest that Tetracladium spp. infection is promoted when bioavailability of P is low, in the same way that plants favour colonisation by arbuscular mycorrhizal symbioses under conditions of low P availability .
In contrast to bacteria, fungal communities are favoured by low soil pH . However, OTUs homologous to T. maxilliforme and T. furcatum showed increased relative abundance in roots as soil pH increased. Furthermore, pH was a key factor in determining fungal community structure in the soil in many cases where Tetracladium has been identified as a common genus [41, 82], however it was found to prefer neutral or slightly acidic soil in a long-term microplot experiment . In addition to pH, soil redox potential (Eh) is an important driver of microbial community growth, diversity and composition [83, 84]. Relative abundance of OTUs homologous to T. maxilliforme and T. furcatum increased as soil moisture decreased, which was surprising considering their dual ecology as aquatic taxa. Low soil moisture combined with high pH leads to lower Eh in the soil and results in slower rates of decomposition [85, 86]. In addition, high soil moisture and high Eh result in increased reducing conditions that limit extractable P availability in the soil through the solubilisation of Fe oxides that bind to available P . Putative root colonising Tetracladium OTUs showed higher relative abundance in low reducing conditions (low soil moisture and high pH) and had higher relative abundance under low P availability, however this data originates from bulk soil physicochemical measurements and therefore the results may differ when looking at the soil closely encapsulating the roots.
Finally, we found a strong correlation between OTU relative abundance and crop management practices. OTUs 19 and 1088 both had the highest relative abundances in all compartments when OSR was planted after rye and their relative abundance increased with OSR rotation length. Crop rotation is known to influence soil microbial community composition, and to influence the composition of plant associated microbiota including pathogens and symbionts [88, 89]. These changes may be attributed to a wide range of interactions. For endophytes, they may reflect plant species specificity, and the extent to which different crop species support proliferation of inoculum, either following recruitment in living root biomass or on organic material left in the field following harvest . Additionally, differences in management practices across crop types, such as fertiliser, tillage and pesticide use could also impact the inoculum . Enhanced colonisation following ryegrass relative to OSR could therefore indicate a preference of Tetracladium spp. for ryegrass as a host, or management practices associated with ryegrass.
The ecological interactions of Tetracladium spp. are a currently unknown; however, there is overwhelming evidence that some taxa within this group which were traditionally considered to be aquatic hyphomycetes can also occur as endophytes in terrestrial ecosystems, with several clades known only from environmental DNA, and which may represent terrestrial species. We found a correlation between crop yield and Tetracladium abundance, indicating that these fungi are signatures of a beneficial plant mycobiome. There is also indication that crop management practices, pH and nutrient enrichment are the main drivers of root colonisation of Tetracladium spp. in terrestrial environments. Further research is needed to determine their role in the plant’s life, particularly their effects on plant health and nutrition, to establish their potential value for utilisation in sustainable agricultural practices.
Availability of data and materials
All data generated or analysed during this study are included in this published article [and its supplementary information files] and in Hilton, S., Picot, E., Schreiter, S. et al. Identification of microbial signatures linked to oilseed rape yield decline at the landscape scale. Microbiome 9, 19 (2021). https://doi.org/10.1186/s40168-020-00972-0. The sequences used for this study are deposited in the NCBI GenBank database under accession numbers MN047188, MN047207, and ON782228-ON782237. The multiple sequence alignment is provided in Additional file 7.
Grossart H-P, Van den Wyngaert S, Kagami M, Wurzbacher C, Cunliffe M, Rojas-Jimenez K. Fungi in aquatic ecosystems. Nat Rev Microbiol. 2019;17:339–54.
Ingold CT. Aquatic hyphomycetes of decaying alder leaves. Trans Br Mycol Soc. 1942;25(4):339–417. https://doi.org/10.1016/S0007-1536(42)80001-7.
Ingold CT. Aquatic hyphomycetes from Switzerland. Trans Br Mycol Soc. 1949;32(3–4):341–5. https://doi.org/10.1016/s0007-1536(49)80026-x.
Johnston PR, Quijada L, Smith CA, Baral HO, Hosoya T, Baschien C, et al. A multigene phylogeny toward a new phylogenetic classification of leotiomycetes. Int Mycol Assoc Fungus. 2019;10:1. https://doi.org/10.1186/s43008-019-0002-x.
Webster J. Experiments with spores of aquatic hyphomycetes: I. Sedimentation, and impaction on smooth surfaces. Ann Bot. 1959;24(4):595–611.
Selosse M-A, Vohník M, Chauvet E. Out of the rivers: are some aquatic hyphomycetes plant endophytes? New Phytol. 2008;178(1):3–7. https://doi.org/10.1111/j.1469-8137.2008.02390.x.
de Wildeman É. Notes mycologiques IV. Annales de Société Belge de Microscopie. 1893;17(2):35–40.
Hibbett DS, Binder M, Bischoff JF, Blackwell M, Cannon PF, Eriksson OE, et al. A higher-level phylogenetic classification of the fungi. Mycol Res. 2007;111(Pt 5):509–47. https://doi.org/10.1016/j.mycres.2007.03.004.
Anderson JL, Shearer CA. Population genetics of the aquatic fungus Tetracladium marchalianum over space and time. PLoS ONE. 2011;6(1):e15908. https://doi.org/10.1371/journal.pone.0015908.
Conway KE. The aquatic hyphomycetes of Central New York. Mycologia. 1970;62(3):516–30. https://doi.org/10.1080/00275514.1970.12018992.
Makela K. Some aquatic hyphomycetes on grasses in Finland. Karstenia. 1973;13:16–22. https://doi.org/10.29203/ka.1973.78.
Chauvet E. Aquatic hyphomycete distribution in south-western France. J Biogeogr. 1991;18(6):699. https://doi.org/10.2307/2845551.
Bärlocher F. The ecology of aquatic hyphomycetes. Berlin: Springer; 1992.
Harrington TJ. Aquatic hyphomycetes of 21 rivers in southern Ireland. Biol Environ Proc R Ir Acad. 1997;97B(2):139–48.
Bandoni RJ. Terrestrial occurrence of some aquatic hyphomycetes. Can J Bot. 1972;50(11):2283–8. https://doi.org/10.1139/b72-297.
Bärlocher F, Oertli JJ. Colonization of conifer needles by aquatic hyphomycetes. Can J Bot. 1978;56(1):57–62. https://doi.org/10.1139/b78-005.
Sati SC, Belwal M. Aquatic hyphomycetes as endophytes of riparian plant roots. Mycologia. 2005;97(1):45–9. https://doi.org/10.1080/15572536.2006.11832837.
Abadie J-C, Puttsepp U, Gebauer G, Faccio A, Bonfante P, Selosse M-A. Cephalanthera longifolia (Neottieae, Orchidaceae) is mixotrophic: a comparative study between green and nonphotosynthetic individuals. Can J Bot. 2006;84(9):1462–77. https://doi.org/10.1139/b06-101.
Stark C, Babik W, Durka W. Fungi from the roots of the common terrestrial orchid Gymnadenia conopsea. Mycol Res. 2009;113(9):952–9. https://doi.org/10.1016/j.mycres.2009.05.002.
Vendramin E, Gastaldo A, Tondello A, Baldan B, Villan M, Squartini A. Identification of two fungal endophytes associated with the endangered orchid Orchis militaris L. J Microbiol Biotechnol. 2010;20(3):630–6.
Wang Y, Wang H, Cheng H, Chang F, Wan Y, She X. Niche differentiation in the rhizosphere and endosphere fungal microbiome of wild Paris polyphylla Sm. PeerJ. 2020;8:e8510. https://doi.org/10.7717/peerj.8510.
Grudzinska-Sterno M, Yuen J, Stenlid J, Djurle A. Fungal communities in organically grown winter wheat affected by plant organ and development stage. Eur J Plant Pathol. 2016;146(2):401–17. https://doi.org/10.1007/s10658-016-0927-5.
Chen Z, Wang Q, Ma J, Zou P, Yu Q, Jiang L. Fungal community composition change and heavy metal accumulation in response to the long-term application of anaerobically digested slurry in a paddy soil. Ecotoxicol Environ Saf. 2020;196:110453. https://doi.org/10.1016/j.ecoenv.2020.110453.
Zhao Y, Fu W, Hu C, Chen G, Xiao Z, Chen Y, et al. Variation of rhizosphere microbial community in continuous mono-maize seed production. Sci Rep. 2021. https://doi.org/10.1038/s41598-021-81228-1.
Bruzone MC, Fontenla SB, Vohnik M. Is the prominent ericoid mycorrhizal fungus Rhizoscyphus ericae absent in the Southern Hemisphere Ericaceae a case study on the diversity of root mycobionts in Gaultheria spp. from northwest Patagonia, Argentina. Mycorrhiza. 2014;25(1):25–40. https://doi.org/10.1007/s00572-014-0586-3.
Macia-Vicente JG, Nam B, Thines M. Root filtering, rather than host identity or age, determines the composition of root-associated fungi and oomycetes in three naturally co-occurring Brassicaceae. Soil Biol Biochem. 2020;146:107806. https://doi.org/10.1016/j.soilbio.2020.107806.
Chatterton S, Yang HE, Ortega Polo R, McAllister TA, Safarieskandari S, Lupwayi N. Bacterial and fungal communities, but not physicochemical properties, of soil differ according to root rot status of pea. Pedobiologia. 2021;84:150705. https://doi.org/10.1016/j.pedobi.2020.150705.
Ma W, Wu Y, Wei Y, Zou W, Yan Y, Xue J, et al. Microbial diversity analysis of vineyards in the Xinjiang region using high-throughput sequencing. J Inst Brew. 2018;124(3):276–83. https://doi.org/10.1002/jib.501.
Sati SC, Arya P, Belwal M. Tetracladium nainitalense sp. Nov., a root endophyte from Kumaun Himalaya, India. Mycologia. 2009;101(5):692–5. https://doi.org/10.3852/08-192.
Giesemann P, Eichenberg D, Stöckel M, Seifert LF, Gomes SIF, Merckx VSFT, et al. Dark septate endophytes and arbuscular mycorrhizal fungi (Paris morphotype) effect the stable isotope composition of ‘classically’ non-mycorrhizal plants. Funct Ecol. 2020;34(12):2453–66. https://doi.org/10.1111/1365-2435.13673.
Russell J, Bulman S. The liverwort Marchantia foliacea forms a specialized symbiosis with arbuscular mycorrhizal fungi in the genus glomus. New Phytol. 2004;165(2):567–79. https://doi.org/10.1111/j.1469-8137.2004.01251.x.
Rosa LH, Almeida-Vieira MDL, Santiago IF, Rosa CA. Endophytic fungi community associated with the dicotyledonous plant Colobanthus Quitensis (Kunth) Bartl. (Caryophyllaceae) in Antarctica. Fed Eur Microbiol Soc Microbiol Ecol. 2010;73(1):178–89. https://doi.org/10.1111/j.1574-6941.2010.00872.x.
Hirose D, Hobara S, Matsuoka S, Kato K, Tanabe Y, Uchida M, et al. Diversity and community assembly of moss-associated fungi in ice-free coastal outcrops of continental Antarctica. Fungal Ecol. 2016;24:94–101. https://doi.org/10.1016/j.funeco.2016.09.005.
Hirose D, Hobara S, Tanabe Y, Uchida M, Kudoh S, Osono T. Abundance, richness, and succession of microfungi in relation to chemical changes in Antarctic moss profiles. Polar Biol. 2017;40(12):2457–68. https://doi.org/10.1007/s00300-017-2157-5.
Klaubauf S, Inselsbacher E, Zechmeister-Boltenstern S, Wanek W, Gottsberger R, Strauss J, et al. Molecular diversity of fungal communities in agricultural soils from lower Austria. Fungal Diversity. 2010;44(1):65–75. https://doi.org/10.1007/s13225-010-0053-1.
Liu Y, Wu L, Wu X, Li H, Liao Q, Zhang X, et al. Analysis of microbial diversity in soil under ginger cultivation. Scientifica. 2017;2017:1–4. https://doi.org/10.1155/2017/8256865.
María T, Castano C, Rodríguez A, Ibanez M, Loboe A, Sebastia M-T. Fairy rings harbor distinct soil fungal communities and high fungal diversity in a montane grassland. Fungal Ecol. 2020;47:100962. https://doi.org/10.1016/j.funeco.2020.100962.
Hilton S, Picot E, Schreiter S, Bass D, Norman K, Oliver AE, et al. Identification of microbial signatures linked to oilseed rape yield decline at the landscape scale. Microbiome. 2021;9(1):19. https://doi.org/10.1186/s40168-020-00972-0.
Bridge PD, Newsham KK. Soil fungal community composition at Mars Oasis, a southern maritime Antarctic site, assessed by PCR amplification and cloning. Fungal Ecol. 2009;2(2):66–74. https://doi.org/10.1016/j.funeco.2008.10.008.
Tsuji M, Fujiu S, Xiao N, Hanada Y, Kudoh S, Kondo H, et al. Cold adaptation of fungi obtained from soil and lake sediment in the Skarvsnes ice-free area, Antarctica. Fed Eur Microbiol Soc Microbiol Lett. 2013;346(2):121–30. https://doi.org/10.1111/1574-6968.12217.
Zhang T, Wang N, Liu H, Zhang Y, Yu L. Soil pH is a key determinant of soil fungal community composition in the Ny-Ålesund Region, Svalbard (High Arctic). Front Microbiol. 2016. https://doi.org/10.3389/fmicb.2016.00227.
Durán P, Barra PJ, Jorquera MA, Viscardi S, Fernandez C, Paz C, et al. Occurrence of soil fungi in Antarctic pristine environments. Front Bioeng Biotechnol. 2019. https://doi.org/10.3389/fbioe.2019.00028.
Jiang N, Li Y, Zheng C, Chen L, Wei K, Feng J, et al. Characteristic microbial communities in the continuous permafrost beside the bitumen in Qinghai-Tibetan Plateau. Environ Earth Sci. 2015;74(2):1343–52. https://doi.org/10.1007/s12665-015-4124-1.
Wang M, Jiang X, Wu W, Hao Y, Su Y, Cai L, et al. Psychrophilic fungi from the world’s roof. Persoonia. 2015;34:100–12. https://doi.org/10.3767/003158515X685878.
Hu W, Schmidt SK, Sommers P, Darcy JL, Porazinska DL. Multiple-trophic patterns of primary succession following retreat of a high-elevation glacier. Ecosphere. 2021. https://doi.org/10.1002/ecs2.3400.
Sati SC, Pant P. Two root endophytic aquatic hyphomycetes Campylospora parvula and Tetracladium setigerum as plant growth promoters. Asian J Agric Res. 2020;4(1):28–33.
Sati SC, Arya P. Assessment of root endophytic aquatic hyphomycetous fungi on plant growth. Symbiosis. 2010;50(3):143–9. https://doi.org/10.1007/s13199-010-0054-6.
Arya P, Sati SC. Evaluation of endophytic aquatic hyphomycetes for their antagonistic activity against pathogenic bacteria. Int Res J Microbiol. 2011;2(9):343–7.
Ihrmark K, Bödeker ITM, Cruz-Martinez K, Friberg H, Kubartova A, Schenck J, et al. New primers to amplify the fungal ITS2 region—evaluation by 454-sequencing of artificial and natural communities. Fed Eur Microbiol Soc Microbiol Ecol. 2012;82(3):666–7.
Caporaso JG, Kuczynski J, Stombaugh J, Bittinger K, Bushman FD, Costello EK, et al. QIIME allows analysis of high-throughput community sequencing data. Nat Methods. 2010;7(5):335–6.
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(21):5271–7.
Blaxter M, Mann J, Chapman T, Thomas F, Whitton C, Floyd R, et al. Defining operational taxonomic units using DNA barcode data. Philos Trans R Soc B. 2005;360(1462):1935–43. https://doi.org/10.1098/rstb.2005.1725.
Katoh K, Rozewicki J, Yamada KD. MAFFT online service: multiple sequence alignment, interactive sequence choice and visualization. Brief Bioinform. 2019;20(4):1160–6. https://doi.org/10.1093/bib/bbx108.
Miller MA, Pfeiffer W, Schwartz T. Creating the CIPRES science gateway for inference of large phylogenetic trees. In: Proceedings of the gateway computing environments workshop (GCE). 2010. pp. 1–8. https://doi.org/10.1109/GCE.2010.5676129.
Stamatakis A. RAxML version 8: a tool for phylogenetic analysis and post-analysis of large phylogenies. Bioinformatics. 2014;30(9):1312–3.
Oksanen J, Blanchet FG, Friendly M, Kindt R, Legendre P, McGlinn D, et al. vegan: community ecology package. 2018. https://CRAN.R-project.org/package=vegan.
McMurdie PJ, Holmes S. phyloseq: An R package for reproducible interactive analysis and graphics of microbiome census data. Public Libr Sci ONE. 2013;8(4):e61217.
Hamilton NE, Ferry M. ggtern: ternary diagrams using ggplot2. J Stat Softw. 2018;87:1–17.
R Core Team. R: a language and environment for statistical computing. Vienna: R Foundation for Statistical Computing; 2021.
Grace JB, Keeley JE. A structural equation model snalysis of postfire plant diversity in California shrublands. Ecol Appl. 2006;16(2):503–14.
Bates D, Mächler M, Bolker B, Walker S. Fitting linear mixed-effects models using lme4. J Stat Softw. 2015;67:1–48.
Lüdecke D, Ben-Shachar MS, Patil I, Waggoner P, Makowski D. performance: an R package for assessment, comparison and testing of statistical models. J Open Source Softw. 2021;6(60):3139.
Lefcheck JS. piecewiseSEM: piecewise structural equation modeling in R for ecology, evolution, and systematics. Methods Ecol Evol. 2016;7(5):573–9.
Fernández-González AJ, Villadas PJ, Cabanás CG-L, Valverde-Corredor A, Belaj A, Mercado-Blanco J, et al. Defining the root endosphere and rhizosphere microbiomes from the world olive germplasm collection. Sci Rep. 2019;9:20423.
Engblom E, Lingdell P-E, Marvanova L, Miiller-Haeckel A. Foam spora in running waters of southern Greenland. Polar Res. 1986;4(1):47–51. https://doi.org/10.1111/j.1751-8369.1986.tb00517.x.
Gönczöl J, Révay Á. Fungal spores in rainwater: stemflow, throughfall and gutter conidial assemblages. Fungal Diversity. 2004;16:67–86.
Orłowska M, Lengiewicz I, Suszycka M. Hyphomycetes developing on water plants and bulrushes in fish ponds. Pol J Environ Stud. 2004;13(6):703–7.
Al-Riyami M, Victor R, Seena S, Elshafie AE, Bõrlocher F. Leaf decomposition in a mountain stream in the sultanate of Oman. Int Rev Hydrobiol. 2009;94(1):16–28. https://doi.org/10.1002/iroh.200811085.
Kravetz S, Gonzalez B, Marno A, Giorgi A. The genus Tetracladium in Pampean streams (Buenos Aires, Argentina). Phytotaxa. 2018;338(3):276. https://doi.org/10.11646/phytotaxa.338.3.5.
Hubbard M, Germida JJ, Vujanovic V. Fungal endophytes enhance wheat heat and drought tolerance in terms of grain yield and second-generation seed viability. J Appl Microbiol. 2013;116(1):109–22.
Molina-Montenegro MA, Acuña-Rodríguez IS, Torres-Díaz C, Gundel PE, Dreyer I. Antarctic root endophytes improve physiological performance and yield in crops under salt stress by enhanced energy production and Na+ sequestration. Sci Rep. 2020;10:5819.
Pańka D, Piesik D, Jeske M, Baturo-Cieśniewska A. Production of phenolics and the emission of volatile organic compounds by perennial ryegrass (Lolium perenne L.)/Neotyphodium lolii association as a response to infection by Fusarium poae. J Plant Physiol. 2013;170(11):1010–9.
Xia C, Li N, Zhang Y, Li C, Nan XZ. Role of Epichloë endophytes in defense responses of cool-season grasses to pathogens: a review. Plant Dis. 2018;102(11):2061–73.
Behie SW, Bidochka MJ. Nutrient transfer in plant-fungal symbioses. Trends Plant Sci. 2014;19(11):734–40. https://doi.org/10.1016/j.tplants.2014.06.007.
Abrego N, Roslin T, Huotari T, Tack AJM, Lindahl BD, Tikhonov G, et al. Accounting for environmental variation in co-occurrence modelling reveals the importance of positive interactions in root-associated fungal communities. Mol Ecol. 2020;29(14):2736–46.
Kia SH, Pallesch S, Piepenbring M, Macia-Vicente JG. Root endophytic fungi show low levels of interspecific competition in planta. Fungal Ecol. 2019;39:184–91.
Spagnoletti FN, Tobar NE, Pardo AFD, Chiocchio VM, Lavado RS. Dark septate endophytes present different potential to solubilize calcium, iron and aluminum phosphates. Appl Soil Ecol. 2017;111:25–32.
Sati SC, Pant P. Evaluation of phosphate solubilization by root endophytic aquatic hyphomycete Tetracladium setigerum. Symbiosis. 2018;77(2):141–5. https://doi.org/10.1007/s13199-018-0575-y.
Wang Y, Ji H, Hu Y, Wang R, Rui J, Guo S. Different selectivity in fungal communities between manure and mineral fertilizers: a study in an alkaline soil after 30 years fertilization. Front Microbiol. 2018. https://doi.org/10.3389/fmicb.2018.02613.
Gosling P, Hodge A, Goodlass G, Bending GD. Arbuscular mycorrhizal fungi and organic farming. Agric Ecosyst Environ. 2006;113(1–4):17–35.
Rousk J, Brookes PC, Bååth E. Contrasting soil pH effects on fungal and bacterial growth suggest functional redundancy in carbon mineralization. Appl Environ Microbiol. 2009;75(6):1589–696.
Grządziel J, Gałązka A. Fungal biodiversity of the most common types of polish soil in a long-term microplot experiment. Front Microbiol. 2019. https://doi.org/10.3389/fmicb.2019.00006.
Heintze SG. The use of the glass electrode in soil reaction and oxidation-reduction potential measurements. J Agric Sci. 1934;24(1):28–41.
Husson O. Redox potential (Eh) and pH as drivers of soil/plant/microorganism systems: a transdisciplinary overview pointing to integrative opportunities for agronomy. Plant Soil. 2012;362:389–417.
Włodarczyk T, Szarlip P, Brzezińska M, Kotowska U. Redox potential, nitrate content and pH in flooded eutric cambisol during nitrate reduction. Res Agric Eng. 2007;53(1):20–8.
Tano BF, Brou CY, Dossou-Yovo ER, Saito K, Futakuchi K, Wopereis MCS, et al. Spatial and temporal variability of soil redox potential, pH and electrical conductivity across a toposequence in the Savanna of West Africa. Agron J. 2020;10:1787.
Phillips IR. Phosphorus availability and sorption under alternating waterlogged and drying conditions. Commun Soil Sci Plant Anal. 1998;29(19–20):3045–59.
Taheri AE, Hamel C, Gan Y. Cropping practices impact fungal endophytes and pathogens in durum wheat roots. Appl Soil Ecol. 2016;100:104–11. https://doi.org/10.1016/j.apsoil.2015.12.007.
Benitez M-S, Osborne SL, Lehman RM. Previous crop and rotation history effects on maize seedling health and associated rhizosphere microbiome. Sci Rep. 2017;7:1–13.
Kirkegaard JA, Simpfendorfer S, Holland J, Bambach R, Moore KJ, Rebetzke GJ. Effect of previous crops on crown rot and yield of durum and bread wheat in northern NSW. Aust J Agric Res. 2004;55(3):321–34.
This project was founded by the UK Research and Innovation Natural Environment Research Council Grants 2433027 and NE/S010270/1 and Biotechnology and Biological Sciences Research Council grant BB/L025892/1.
Ethics approval and consent to participate
Consent for publication
The authors declare no competing financial interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Correlogram showing Pearson’s correlation of all metadata variables and OTU relative abundance. Boxes are coloured according to the R values blue indicating a negative, red indicating a positive relationship. Non-significant relations are shown with an x in the box.
Model indicators for LMMs. A is a visual representation of the measured model fit indices for the nutrient models. The final model included in Figure 6 is modMr. B is the actual values corresponding to the indicators of the nutrient models. C is the model fit indicators of the soil structure and climate models from Figure 6.
Sequencing efficacy of the samples for the Tetracladium sp. OTUs. Rarefaction curves showing Tetracladium sp. OTU richness across all samples in the bulk soil, the rhizosphere, and the roots.
Metadata tables. A - Metadata table from the 25 farms. B – Bulk soil properties from each five reps of the 37 field sites. C - Tetracladium OTU relative abundances across samples. D – Mean Tetracladium OTU relative abundances across the three sampled compartments. E – Closest sequence matches in the UNITE database with DOIs.
Dunn’s test results correcting for P values with Benjamini–Hochberg method. Comparing previous crop types across the different compartments in OTU 19 and OTU 1088.
Dunn’s test results correcting for P values with Benjamini–Hochberg method. Comparing OTU 19 and OTU 1088 relative abundance in different oilseed rape varieties in the roots.
Multiple sequence alignment of the twelve Tetracladium sp. OTUs and the reference sequences.
About this article
Cite this article
Lazar, A., Mushinski, R.M. & Bending, G.D. Landscape scale ecology of Tetracladium spp. fungal root endophytes. Environmental Microbiome 17, 40 (2022). https://doi.org/10.1186/s40793-022-00431-3
- Aquatic hyphomycete
- Root colonising