Next Article in Journal
Effects of Truffle Inoculation on Root Physiology and Mycorrhizosphere Microbial Communities of Carya illinoinensis Seedlings
Previous Article in Journal
Potential Distribution and Suitable Habitat for Chestnut (Castanea sativa)
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Fungi Inhabiting Stem Wounds of Quercus robur following Bark Stripping by Deer Animals

by
Adas Marčiulynas
1,*,
Vaida Sirgedaitė-Šėžienė
1 and
Audrius Menkis
2
1
Institute of Forestry, Lithuanian Research Centre for Agriculture and Forestry, Liepų Str. 1, Kaunas District, 53101 Girionys, Lithuania
2
Department of Forest Mycology and Plant Pathology, Uppsala BioCenter, Swedish University of Agricultural Sciences, P.O. Box 7026, SE-75007 Uppsala, Sweden
*
Author to whom correspondence should be addressed.
Forests 2023, 14(10), 2077; https://doi.org/10.3390/f14102077
Submission received: 27 August 2023 / Revised: 29 September 2023 / Accepted: 16 October 2023 / Published: 17 October 2023
(This article belongs to the Section Forest Biodiversity)

Abstract

:
We investigated fungal communities in oak wounds to determine how fungal species richness and community composition changes depending on the age of wounds. The sampling of wood cores was carried out from 10-, 20-, 30-, 40-, and 50-year-old wounds. The fungal community was analyzed using high-throughput sequencing of ITS2 rDNA. Sequence analysis showed the presence of 534 fungal OTUs, which were 83.4% Ascomycota, 16.3% Basidiomycota, and 0.3% Mucoromycota. The fungal OTU richness changed over time: it increased as compared between 10- and 20-year-old wounds, remained similar in 20- to 40-year-old wounds, and decreased in 50-year-old wounds. The fungal community composition also changed over time with the largest differences detected between 10-year-old and older wounds (p < 0.001). The most common representatives of Basidiomycota were Laetiporus sulphureus (34.7%), Mycena galericulata (17.0%), and Cylindrobasidium evolvens (6.5%), and the most common of Ascomycota were Aposphaeria corallinolutea (13.6%), Sclerostagonospora cycadis (7.6%), and Cadophora malorum (5.8%). In conclusion, oak wounds of different ages were colonized by a high diversity of fungi including oak-associated species. Fungal communities in oak wounds underwent qualitative and quantitative changes over time, which led to the gradual shift from fungal generalists in young wounds to oak specialists in older wounds.

1. Introduction

The common oak (Quercus robur) is a tree species with the highest diversity of associated species in Northern Europe [1,2]. Oaks are long-lived trees, and this longevity leads to many different micro-habitats being formed on and inside the trees, which promotes the establishment of vast biodiversity. According to some estimates, oaks in general can provide habitats to nearly 900 different species [3,4]. Among these, 252 species of phytophagous insects and mites have been identified as oak specialists. In the Nordic countries, about 770 oak-related species are included in the national red list, about 400 of which are in wood, including about 300 beetles and 50 fungi [5]. As oaks age, the activity of wood-decay fungi and wood-inhabiting insects results in decayed wood, which also provides a habitat for larger organisms such as small mammals (e.g., dormouse, bats) or birds [6]. However, the decline of oak trees over the last century in Europe has also resulted in the decline of many oak-associated species [7,8]. Although old oaks still exist in Europe, in many places, their populations are isolated, resulting in unfavorable conditions for the long-term survival of associated organisms [9,10,11,12]. In Lithuania, old oak trees make up a significant part of European oak habitats, where endangered and protected species of organisms are found. The starting point for the formation of valuable oak habitats is usually reached at about the age of 200 years [13]. In oaks older than 200 years, there is a rapid increase in both hollows [14] and colonization of trees by epiphytic lichens [15]. Hollows in trees younger than 200 years old are usually rare, which is likely due to factors such as a high tree vigor and/or occasional colonization of wood by wood-decay fungi [16]. In many areas, there is a lack of oak trees older than 100 years [17]. In Lithuania for example, the average age of oak stands is around 80 years. Moreover, trees growing in managed forests are often harvested at the age of 120 years, i.e., long before the time when valuable habitats appear on and in these trees. Therefore, the limited availability of old-growth oaks calls for alternative solutions to maintain associated biodiversity.
Tree veteranization is one means used for conservation and maintenance of tree-associated biodiversity [18,19,20]. The purpose of veteranization is to shorten the long waiting period until habitats suitable for protected and rare species are naturally formed. To achieve this goal, trees are damaged using tools to mimic natural processes in nature by creating damages such as a woodpecker hole or a broken branch, thereby creating habitats with exposed wood suitable for the establishment of new species [18]. Such measures are usually used where the number of old trees is not sufficient to ensure the preservation of biodiversity. For this, younger trees are usually used due to both their availability and the lack of exposed wood habitats needed for the natural establishment of particular species [21]. However, tree veteranization should be limited in areas with a relatively large number of large animals causing natural damages to younger oak trees.
Large herbivores, especially several species of deer (red deer or moose), cause significant damages to forestry due to browsing and bark stripping [22,23,24]. These animals generally prefer deciduous trees, and oaks are particularly frequently browsed [25,26]. Most often, these animals damage the stem of trees at a relatively early age, when oaks reach the age of 15–30 years, and when the bark is still relatively young and has not begun to roughen. In Lithuania, damages caused by deer animals encompass up to 1000 hectares of different forest stands each year. About 10%–15% of these damages are in young oak stands, while the remaining are in Picea abies and Pinus sylvestris stands [27]. Although severely damaged oak trees do not recover, the vast majority remain viable and are able to close wounds over time (Figure 1).
Tree vigor is a phenotypic trait that refers to the ability of a tree to grow and remain healthy in the prevailing environment [28]. Trees may also have well-functioning defense mechanisms that allow them to effectively resist attacks by native agents of damage [29]. Living and healthy trees are often colonized by fungal endophytes, which are fungi with different functional roles that colonize and persist in different tissues asymptomatically until environmental conditions become favorable for their development [30,31,32], i.e., when they can begin to respond to changes in substrate quality [33]. These endophytes can influence many biotic processes in trees, and under certain conditions, can shift to a pathogenic or saprophytic lifestyle [34]. They can act as first decomposers, thereby influencing the spatial and temporal active layer of biodiversity in their hosts [35]. However, wounds caused by animals and the presence of exposed wood facilitate the establishment of many different fungi due to new colonization opportunities [36,37]. In young wounds, there are often appropriate conditions for plant pathogens and/or wood-decay fungi to establish. Most often, these fungi are divided into three groups: obligate pathogens that interact only with living tree cells, facultative pathogens that live in both living and dead tree tissue, and characteristic saprophytes that only affect dead tree tissue [38]. Exposed wood can be infected by pathogenic or wood-decay fungi through fungal spores that are abundant in the forest [39]. In a common case, the transition of wood pathogens from the wound surface to deeper wood layers is limited due to the activity of wood preservatives and the limited availability of oxygen. Bacterial wet wood or a slow type of rot called “soft-rot” can develop in such wounds in the presence of low oxygen levels [40]. Although a relatively rapid wound closure may prevent further infections, closed wounds can eventually open due to wood decay caused by brown rot fungi, thereby in young oaks creating suitable habitats for many oak-associated organisms. Fungi such as Cadophora malorum, Panellus stipticus, and various species of Ophiostoma, Helotiales, or Cystofilabasidiales were often found in newly produced oak wounds, while Trametes versicolor, Peniophorella pubera, Phlebia radiata, Mollisia sp., or Stereum sp. were reported from older wounds [24,41,42].
The aim of this study was to investigate fungal communities in oak stem wounds caused by deer animals using high-throughput sequencing. We also aimed to determine how fungal species diversity changes depending on the wound age. We hypothesized that (a) wounds caused by deer animals in young oaks create conditions for the establishment of oak-associated fungal species, which can persist in colonized wood; and (b) the fungal species richness increases with the age of wounds.

2. Materials and Methods

2.1. Study Sites and Sampling

The study sites were at Kėdainiai (K), Prienai (P), and Telšiai (T), which were in three different regions of Lithuania (Figure 2). These regions have the highest density of Q. robur in forest stands. At each site, Q. robur trees, which were damaged by deer animals (Figure 1), were between 20 and 60 years old. The age of the trees and the age of the wounds were determined by drilling them with an increment borer (Haglöf, Sweden) and counting the annual tree rings. By comparing the difference between tree rings in the wound and in the intact wood, the age of the wound was determined. Wounds were often 15–20 cm long and 10 cm wide. Wounds that were 10 and 20 years old were not closed, but older wounds were either partially or completely closed. All 50-year-old wounds were closed. In each sampling site (K, P, or T), the distances between the individual oak stands, which were used for sampling, were up to 20 km. Different oak stands often differed in herbaceous vegetation, but other woody vegetation was often absent as pure stands of Q. robur prevailed.
The sampling of oak wood was carried out in 2021. Wood cores were taken from wounds that were either 10 (1 to 10 y), 20 (11–20 y), 30 (21–30 y), 40 (31–40 y), or 50 (41–50 y) years old. In each site, wood samples from wounds of the same age were collected from five trees growing within the same stand. Wood samples were taken ca. 1.2 m from the ground, i.e., at the height where the deer damage prevails. Prior to sampling, the upper layer (2–4 cm thick) of the wood, or bark covering the wound, was removed using an axe, and three replicate samples (ca. 7 cm in depth) were taken from each tree by drilling into the wound using an increment borer. Wood sampling tools were thoroughly cleaned between individual trees using 96% ethanol. A total of 75 oak trees were sampled (3 regions × 5 age classes × 5 trees). The collected wood cores were individually placed in sterile plastic tubes, labelled, transported the same day to the laboratory, and stored at −20 °C before further processing.

2.2. DNA Isolation, Amplification, and Sequencing

Before the isolation of DNA, three wood cores from the same tree were put together. The DNA work was based on the study by Marčiulynas et al. [43]. Firstly, wood samples were grinded in liquid nitrogen and ca. 0.5 g of wood was placed into a 2 mL screw-cap tube with two (2 mm in diameter) metal beads. Then, samples were homogenized using a Fast prep shaker (Montigny-le-Bretonneux, France), and the DNA was extracted using CTAB buffer and incubated at 65 °C for 1 h (vortex every 15 min). The supernatant was mixed with an equal volume of chloroform and cleaned with 2-propanol. The pellet was washed in 500 μL 70% ethanol, dried, and dissolved in 30 μL sterile milli-Q water. The DNA was further purified using a NucleoSpin®Soil kit (Macherey-Nagel GmbH & Co. Duren, Germany) according to recommendations by the producer. Following the extraction and purification of the DNA, its concentration was measured using a NanoDrop™ One spectrophotometer (Thermo Scientific, Rodchester, NY, USA) and adjusted to 10 ng/mL.
The amplification of the ITS2 rDNA region was achieved using a primer pair gITS7 [44] and ITS4 [45], both containing sample identification barcodes. The polymerase chain reaction (PCR) was performed in 50 μL reactions and consisted of the following final concentrations: 200 μM dNTPs, 750 μM MgCl2, 0.025 μM DreamTaq Green polymerase (5 U/μL) (Thermo Scientific, Waltham, MA, USA), 0.02 ng/μL template DNA, and 200 nM of each primer; sterile milli-Q water to a final volume of 50 μL. The amplifications were carried out using an Applied Biosystems 2720 thermal cycler (Applied Biosystems, Foster City, CA, USA). The following thermocycling pattern was used: 95 °C for 5 min; 95 °C for 30 s, 55 °C for 30 s, and 72 °C for 1 min (30 cycles); and 72 °C for 7 min. The PCR products were examined using gel electrophoresis on 1.5% agarose gels stained with GelRed (Biotium, Fremont, CA, USA). The PCR products were cleaned using a 96% ethanol and 3 M sodium acetate (pH 5.2) (Applichem GmbH, Darmstadt, Germany) mixture (20:1). After the quantification of PCR products using a Qubit fluorometer 4.0 (Life Technologies, Stockholm, Sweden), an equimolar mix was produced, which was sequenced using a PacBio platform and one Sequel II SMRT cell at a SciLifeLab facility in Uppsala, Sweden.

2.3. Bioinformatics

Obtained sequence reads were quality filtered in the Sequence Clustering and Analysis of Tagged Amplicons (SCATA) bioinformatics tool available at http://scata.mykopat.slu.se/ (accessed on 5 July 2023). Short sequences (<200 pb), and those with low read quality (Q < 20), and primer dimers were removed. Sequences without a tag or primer were also excluded. High-quality sequences were clustered into different OTUs using single linkage clustering based on 98.5% similarity. The GenBank (NCBI) database and Blastn algorithm were used to determine taxonomic identities of different OTUs. Taxon-delimiting ITS homology was 98%–100% at the species level, 94%–97% at the genus level, and 80%–94% at the higher level, corresponding to at least 90% of the sequence length [46]. Representative sequences of fungal OTUs were deposited in GenBank under accession numbers OR481117–OR481666. Fungal functional groups were assigned using the FUNGuild fungal database [47].

2.4. Statistical Analysis

Rarefaction analysis was implemented using Analytical Rarefaction v.1.3, (http://www.uga.edu/strata/software/index.html) (accessed on 9 July 2023). Differences in fungal OTU richness (the number of different OTUs) in Q. robur wounds of different ages were compared using the non-parametric chi-square test [43]. The Shannon diversity index (α-diversity) and Sørensen qualitative similarity (β-diversity) index were calculated using SAS v. 9.4 (Cary, NC, USA) [48,49]. The non-parametric Mann–Whitney test in SAS was used to assess whether the Shannon diversity index differed among wounds of different ages. The composition of fungal communities in wounds of different ages was analyzed using nonmetric multidimensional scaling (NMDS) based on the Bray–Curtis similarity index (β-diversity). One-way ANOSIM analysis was carried out to assess for significant differences between different samples. These analyses were carried out using Vegan 2.5.7 and Stats 3.6.2 in R 4.1.1 (https://www.r-project.org, accessed on 19 July 2023) [50,51].

3. Results

PacBio sequencing and quality filtering revealed the presence of 141,671 high-quality sequences, representing 629 OTUs. Taxonomic identification showed that 534 OTUs were fungal (Table S1) and 95 were non-fungal, which were removed from further analyses. The detected fungi were 83.4% Ascomycota, 16.3% Basidiomycota, and 0.3% Mucoromycota. Although 10-year-old wounds had the highest absolute richness of fungal OTUs (all samples taken together, Table 1), when the number of fungal sequences had been taken into consideration, the OTU richness was lower as compared to older wounds (e.g., 10 y wounds had 213 OTUs among 62,406 sequences vs. 20 y wounds had 159 OTUs among 10,427 sequences) (p < 0.05) (Table 1). The OTU richness in 50-year-old wounds was significantly lower than in 20-, 30-, or 40-year-old wounds (p < 0.05), which did not differ significantly from each other (p < 0.05). Differences in OTU richness were also demonstrated by species accumulation curves (Figure 3). The species accumulation curve from 10-year-old wounds was approaching the asymptote, while those from 20–50-year-old wounds did not reach the asymptote, showing that in these wounds, a higher OTU richness could be detected with deeper sequencing (Figure 3). The Shannon diversity index (α-diversity) was moderate in wounds of different ages and ranged between 2.51 and 2.93 (Table 1). Consequently, no significant differences in this respect were found between wounds of different ages (p > 0.05).
There were certain variations in OTU richness among different sampling sites and wounds of different ages (Figure 4). In general, the highest richness of fungal OTUs was in wounds from Kėdainiai, while the lowest was in wounds from Prienai. Regarding the age of wounds, the highest variation in OTU richness was in 30- and 40-year-old wounds from Kėdainiai and in 10-, 20-, and 50-year-old wounds from Telšiai (Figure 4).
The most common fungi in 10-year-old wounds were Aposphaeria corallinolutea (23.7%), Sclerostagonospora cycadis (13.3%), and Diatrype stigma (9.6%); in 20-year-old wounds, these were Cadophora malorum (22.3%), Diaporthe eres (15.6%), and Unidentified sp. 5631_30 (15.4%); in 30-year-old wounds, these were Laetiporus sulphureus (22.0%), Menispora sp. 5631_29 (12.2%), and Phialemonium inflatum (10.6%); in 40-year-old wounds, these were Fimetariella rabenhorstii (30.7%), Leptosillia macrospora (18.3%), and Mycena galericulata (16.5%); and in 50-year-old, these were Phaeostalagmus cyclosporus (20.1%), L. sulphureus (16.6%), and C. malorum (10.3%) (Table 2).
The most common Basidiomycota fungi were L. sulphureus (34.7%), M. galericulata (17.0%), and Cylindrobasidium evolvens (6.5%), and the most common Ascomycota fungi were A. corallinolutea (13.6%), S. cycadis (7.6%), and C. malorum (5.8%) (Figure 5).
The number of unique OTUs (found only in a particular age class) varied substantially in stem wounds of different ages (Table S1). The most frequent unique OTUs in 10-year-old wounds were Cylindrobasidium evolvens (2.2%), Peniophora incarnata (0.6%), and Hormodochis aggregata (0.4%); in 20-year-old wounds, these were Candida norvegica (1.2%), Unidentified sp. 5631_140 (0.4%), and Hyphodiscus luxurians (0.3%); in 30-year-old wounds, these were Phlebia acerina (7.5%), Unidentified sp. 5631_39 (5.8%), and Stereum hirsutum (5.3%); in 40-year-old wounds, these were Toniniopsis subincompta (2.6%), Zygosaccharomyces parabailii (2.1%), and Unidentified sp. 5631_56 (2.0%); and in 50-year-old wounds, these were Lophiostoma sp. 5631_40 (4.0%), Phleogena sp. 5631_45 (2.2%), and Mortierella hypsicladia (1.2%) (Table S1). Among all fungi, unique OTUs in wounds of different ages constituted between 24.0% and 52.3% of all OTUs. The highest number of shared OTUs was between 30- and 40-year-old (86 OTUs) wounds, while the least number was between 10- and 50-year-old (37 OTUs) wounds. Most of the unique OTUs were in 10-year-old (112) and 50-year-old (71) wounds (Figure 6).
The Sørensen similarity index (β-diversity) of fungal communities between wounds of different ages was moderate to low. The lowest Sørensen similarity index of fungal communities was between 10- and 50-year-old (0.19) wounds, while the highest was between 30- and 40-year-old (0.46) wounds (Table 3).
Fungal functional groups were determined for 91.4% of fungal sequences, and their relative abundance in wounds of different ages is shown in Figure 7. In all samples, the most abundant fungal functional groups were undefined saprotrophs (37.9%), wood saprotrophs (24.2%), plant pathogens (14.2%), and endophytes (12.7%) (Figure 7). Undefined saprotrophs showed the highest relative abundance in 50-year-old (55.8%), 10-year-old (37.8%), and 30-year-old (37.7%) wounds, and their relative abundance differed significantly from those in 20-year-old wounds (p < 0.05). The highest relative abundance of wood saprotrophs was in 30-year-old (32.7%), 10-year-old (27.8%), and 50-year-old (23.2%) wounds, and their relative abundance differed significantly from those in 20-year-old wounds (p < 0.05) (Figure 7). The relative abundance of plant pathogens was highest in 10-year-old (24.6%) and 20-year-old (22.7%) wounds, which was significantly higher than in wounds of other ages (p < 0.05) (Figure 7). The relative abundance of endophytes varied substantially in wounds of different ages, with the highest difference being between 10-year-old (4.5%) and 20-year-old (38.3%) wounds (p < 0.05).
The detected wood-decay fungi were L. sulphureus (5.7%) (detected in ≥30-year-old wounds), Armillaria gallica (0.9%) (only in 50-year-old wounds), Phlebia acerina (0.9%) (only in 30-year-old wounds), Stereum hirstum (0.6%) (in ≥30-year-old wounds), Stereum rugosum (0.5%) (in ≥30-year-old wounds), and Tremetes versicolor (0.03%) (in ≥40-year-old wounds). Fungi strongly associated with oaks such as Phlebia acerina (1.2%) (only in 30-year-old wounds) and Pleurotus dryinus (0.7%) (mostly in 50-year-old wounds) were also detected (Table S1).
NMDS analysis revealed a partial overlap of fungal communities in oak wounds of different ages (Figure 8). The permutation test showed that fungal communities in 10-year-old wounds differed significantly from those in older wounds (20 y, 30 y, 40 y, and 50 y) (p < 0.001). There was also a significant difference between fungal communities in 20-year-old and 50-year-old wounds (p < 0.01), and in 40-year-old and 50-year-old wounds (p < 0.05), while in wounds of other ages, fungal communities were similar (p > 0.05) (Figure 8).

4. Discussion

In most cases, living tree tissues beneath the intact bark represent an unfavorable environment for colonization by many fungi, but these tissues become readily available for colonization when the tree is injured, or in the case of bark dieback. Post-damage fungal colonization in angiosperms was extensively studied in North America [43,52], showing that fungal colonization is often determined by both the nature of the wound and its location on the stem. If the wound is only superficial, i.e., only the bark is removed in a small area, then conditions are often limited for fungal colonization. However, if the wound is relatively large, and especially if it penetrates several or more annual rings, the decay can spread more rapidly into the inner wood. The primary colonization of wounds starts with latent fungi that were already in the wood, but it may also include opportunistic pathogens, or fungi whose spores were on the damaged surface [13]. Following the expansion of these fungi in wounds of living trees, complex fungal communities are formed. In addition to latent fungi, wounds are also colonized by those with ruderal characteristics, even though their abundance is usually not high [42].
Fungal communities colonizing living tree stems were shown to persist in deadwood for many years after the tree injury, but they are also constantly changing both temporally and spatially [42,53]. In agreement, the results of the present study showed that the richness and composition of fungal communities changed over time. Consequently, certain trends in richness of fungal OTUs in oak wounds of different ages (10 y–50 y) were revealed, showing that there was an increase in OTU richness between 10- and 20-year-old wounds, the OTU richness remained similar between 20- and 40-year-old wounds, and the OTU richness decreased in 50-year-old wounds (Table 1, Figure 3), thereby only partly supporting the hypothesis that the species richness increases with the age of wounds. Although it was shown that the diversity of fungal species changes strongly during different stages of wood decomposition [54], it may be different in wounds. For example, the increase in richness of fungal OTUs in relatively young wounds can be associated with the availability of new substrate for colonization, while competition among fungi already established in oak wood and the ongoing process of wound closure could contribute to minor changes in the richness of fungal OTUs in mid-age (20- to 40-year-old) wounds. Furthermore, the complete wound closure in old wounds can drive the decline of fungal OTU richness due to both the lack of access to exposed wood and reduced supply of oxygen. However, variation in richness of fungal OTUs was high at different study sites (Figure 4), showing a strong impact of local environmental conditions.
Differences in fungal community composition in wounds of different ages were also notable (Figure 8), and these differences were further supported by low values of the Sørensen similarity index (Table 3) and a low overall number of shared OTUs (21) (Figure 6). These changes, which were likely driven by changes in habitat availability and quality, led to the establishment and persistence of oak-associated fungi, thereby supporting the hypothesis that wounds in young oaks create conditions for the establishment and long-term survival of such fungi. Among these, there were fungi such as L. sulphureus, which causes a heart rot in living oak trees [55,56]; a widespread fungus P. dryinus, which causes a white rot [57]; and the wood-decay fungi S. hirstum, S. rogosum, and Peniophora incarnata [58,59,60]. A tree pathogen A. gallica was also detected in old (especially in 50-year-old) wounds, showing that it is a secondary plant parasite that usually causes infections only after the defense of host trees is weakened, e.g., after the insect defoliation, drought, or infection by other fungi. Ophiostoma quercus, which was suggested to be associated with the decline of oaks in Central Europe [61,62], was also detected. All these oak-associated wood-decay and/or pathogenic fungi were found to be more abundant in older wounds (≥30 years old) (Table 2 and Table S1). By contrast, endophytic and/or generalist pathogenic fungi were more common in 10-year-old wounds (Table 2). For example, in 10-year-old wounds, there were fungal endophytes such as S. cycadis and Paraphaeosphaeria sporulosa [63,64,65], and plant pathogens Neoleptosphaeria rubefaciens [66,67,68] and Sporocadus trimorphus [69] (Table 2), but their functional role in colonized oak wood is obscure.
The abundance and composition of fungal functional groups has also changed over time (Figure 7), as tree wounding enabled latent colonizers of intact wood to switch to saprophytic or pathogenic lifestyles while at the same time having a competitive advantage over other later arriving fungi due to readily available wood resources [70,71]. For example, Lophiostoma corticola, which is a widely distributed saprotroph, showed a high relative abundance in 10-year-old wounds [72,73]. Although it is not common in oaks, it is well characterized from leaves, petioles, or branches of Fraxinus excelsior [73,74,75]. Diatrype stigma was also found among the most common fungi associated with 10-year-old wounds. It is a saprotrophic species that is commonly associated with the genus Quercus but can also grow on various plants within Rosaceae and Betulaceae [76]. Cytospora predappioensis was also discovered in relatively young wounds. Some Cytospora species are known to cause canker diseases in woody plants [77,78]. Cylindrobasidium evolvens is a widespread saprotroph that is common in young wounds of hardwood tree species [79]. Such a composition of fungi in 10-year-old wounds indicates that newly created wounds provide favorable conditions for the establishment and expansion of primary colonizers, which are widely distributed in the environment [32]. As a possible result of expansion of such fungi, the relative abundance of endophytes has peaked in 20-year-old wounds (Figure 7) with Diaporthe eres and Unidentified sp. 5631_30 being the most common and unique fungi (Table 2). Although D. eres can cause stem canker, stem necrosis, dead branch, shoot blight, fruit rot, leaf spot, leaf necrosis, and crown browning [80,81,82], in Mediterranean oaks it was characterized as an endophyte [83]. In 30-year-old and older wounds, there was a notable change from endophytic to saprotrophic fungi (Figure 7), which can probably be attributed to the activity of secondary wood decomposers [84,85]. For example, M. galericulata and Ascocoryne sarcoides were commonly detected in older wounds (Table 2). Mycena galericulata is a saprotrophic wood-decay fungus, which occurs at later stages of wood decomposition, as it requires a substrate prepared by other fungi [86,87]. Ascocoryne sarcoides is commonly found in high humidity temperate forests, which are rich in beech and oak [88,89]. Although the fungus usually occurs on fallen trees or large pieces of wood, suggesting that it can be a saprotroph, the precise trophic mode is obscure [90].

5. Conclusions

Oak wounds of different ages are inhabited by a high diversity of fungi including oak-associated species such as L. sulphureus, P. dryinus, and P. incarnata, which often persist in old wounds after their closure, demonstrating that wounds in young oaks create conditions for the establishment and long-term survival of oak-associated fungi. The fungal OTU richness changed over time: it increased between 10- and 20-year-old wounds, remained similar between 20- and 40-year-old wounds, and decreased in 50-year-old wounds. Consequently, fungal communities in oak wounds undergo qualitative and quantitative changes over time, leading to the gradual shift from fungal generalists in young wounds to oak specialists in older wounds.

Supplementary Materials

The following Supporting Information can be downloaded at: https://www.mdpi.com/article/10.3390/f14102077/s1, Table S1: Distribution and relative abundance of fungal OTUs (shown as a proportion of all high-quality fungal sequences) in Quercus robur wounds of different ages (10–50 years old). The data from different study sites are combined.

Author Contributions

Conceptualization, A.M. (Adas Marčiulynas) and A.M. (Audrius Menkis); methodology, A.M. (Adas Marčiulynas), V.S.-Š. and A.M. (Audrius Menkis); software, A.M. (Adas Marčiulynas), V.S.-Š. and A.M. (Audrius Menkis); validation, A.M. (Adas Marčiulynas) and A.M. (Audrius Menkis); formal analysis, A.M. (Adas Marčiulynas), V.S.-Š. and A.M. (Audrius Menkis); investigation, A.M. (Adas Marčiulynas) and V.S.-Š.; resources, A.M. (Adas Marčiulynas) and A.M. (Audrius Menkis); data curation, A.M. (Adas Marčiulynas) and A.M. (Audrius Menkis); writing—original draft preparation, A.M. (Adas Marčiulynas) and V.S.-Š.; writing—review and editing, A.M. (Audrius Menkis); visualization, A.M. (Adas Marčiulynas) and A.M. (Audrius Menkis); supervision, A.M. (Audrius Menkis). All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the European Social Fund under the 09.3.3-LMT-K-712 “Development of Competences of Scientists other Researchers and Students through Practical Research Activities” measure.

Data Availability Statement

Not applicable.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Leuschner, C.; Ellenberg, H. Ecology of Central European Non-Forest Vegetation: Coastal to Alpine, Natural to Man-Made Habitats: Vegetation Ecology of Central Europe; Springer: Berlin/Heidelberg, Germany, 2017; Volume 2, 1094p. [Google Scholar] [CrossRef]
  2. Mölder, A.; Meyer, P.; Nagel, R.V. Integrative management to sustain biodiversity and ecological continuity in Central European temperate oak (Quercus robur, Q. petraea) forests: An overview. For. Ecol. Manag. 2019, 437, 324–339. [Google Scholar] [CrossRef]
  3. Jansson, N.; Bergman, K.O.; Jonsell, M.; Milberg, P. An indicator system for identification of sites of high conservation value for saproxylic oak (Quercus spp.) beetles in southern Sweden. J. Insect Conserv. 2009, 13, 399–412. [Google Scholar] [CrossRef]
  4. Sundberg, S.; Carlberg, T.; Sandström, J.; Thor, G. Värdväxters Betydelse för Andra Organismer–Med Fokus på Vedartade Värdväxter; SLU Artdatabanken: Uppsala, Sweden, 2019; 52p. [Google Scholar]
  5. Tingstad, L.; Grytnes, J.A.; Felde, V.A.; Juslén, A.; Hyvärinen, E.; Dahlberg, A. The potential to use documentation in national Red Lists to characterize red-listed forest species in Fennoscandia and to guide conservation. GECCO 2018, 15, e00410. [Google Scholar] [CrossRef]
  6. Kõrkjas, M.; Remm, L.; Lõhmus, A. Development rates and persistence of the microhabitats initiated by disease and injuries in live trees: A review. For. Ecol. Manag. 2021, 482, 118833. [Google Scholar] [CrossRef]
  7. Sunhede, S.; Vasiliauskas, R. Wood and bark inhabiting fungi on oak in Lithuania. Balt. For. 1996, 2, 23–27. [Google Scholar]
  8. Sunhede, S.; Vasiliauskas, R. Ecology and decay pattern of Inocutis dryophila on Quercus robur. Karstenia 2003, 43, 45–53. [Google Scholar] [CrossRef]
  9. Bergman, K.-O. Bedömning av långsiktig överlevnad för hotade arter knutna till ekar på Händelö. Rapport till Norrköpings kommun. Nat. Norrköping 2003, 3, 1–32. [Google Scholar]
  10. Hedin, J. Metapopulation Ecology of Osmoderma eremita—Dispersal, Habitat Quality and Habitat History. Ph.D. Thesis, Lunds Universitet, Lund, Sweden, 2003. [Google Scholar]
  11. Ranius, T. Extinction risk in metapopulations of a beetle inhabiting hollow trees predicted from time series. Ecography 2007, 30, 716–726. [Google Scholar] [CrossRef]
  12. Bengtsson, V.; Fay, L. 2009 Veteran Tree Survey, Ashtead and Epsom Commons; City of London Report; City of London Open Spaces Department: London, UK, 2009. [Google Scholar]
  13. Menkis, A.; Redr, D.; Bengtsson, V.; Hedin, J.; Niklasson, M.; Nordén, B.; Dahlberg, A. Endophytes dominate fungal communities in six-year-old veteranisation wounds in living oak trunks. Fun. Ecol. 2022, 59, 101020. [Google Scholar] [CrossRef]
  14. Ranius, T.; Niklasson, M.; Berg, N. Development of tree hollows in pedunculate oak (Quercus robur). For. Ecol. Manag. 2009, 257, 303–310. [Google Scholar] [CrossRef]
  15. Ranius, T.; Eliasson, P.; Johansson, P. Large-scale occurrence patterns of red-listed lichens and fungi on old oaks are influenced both by current and historical habitat density. Biodivers. Conserv. 2008, 17, 2371–2381. [Google Scholar] [CrossRef]
  16. Gauthier, N.W.; Fountain, W.E.; Missun, T. Tree Wounds—Invitations to Wood Decay Fungi; Plant Pathology Fact Sheet, Collage of Agriculture, Food and Environment, University of Kentucky: Lexington, KY, USA, 2015; Available online: https://plantpathology.ca.uky.edu/files/ppfs-or-w-01.pdf (accessed on 22 July 2023).
  17. Niklasson, M.; Nilsson, S.G. Skogsdynamik och arters bevarande. Studentlitteratur 2005. Available online: https://www.studentlitteratur.se/kurslitteratur/naturvetenskap-och-miljo/miljo/skogsdynamik-och-arters-bevarande (accessed on 25 July 2023).
  18. Read, H. Veteran Trees—A Guide to Good Management; English Nature: Peterborough, UK, 2000. [Google Scholar]
  19. Lonsdale, D. Ancient and other Veteran Trees: Further Guidance on Management; The Tree Council: London, UK, 2013. [Google Scholar]
  20. Hedin, J.; Niklasson, M.; Bengtsson, V. Veteranisering–Verktyg i Stället för Tid. Fauna Flora 2018, 113, 12–25. Available online: https://lifebridgingthegap.se/wp-content/uploads/VeteraniseringFochF2018Hedin-et-al.pdf (accessed on 20 July 2023).
  21. Bengtsson, V.; Hedin, J.; Niklasson, M. Veteranisation of oak—Managing trees to speed up habitat production. In Trees Beyond the Wood: An Exploration of Concepts of Woods, Forests and Trees; Rotherham, I.D., Handley, C., Agnoletti, M., Samojlik, T., Eds.; Wildtrack Publishing: Sheffield, UK, 2013; pp. 61–68. [Google Scholar]
  22. Gill, R.M.A. A review of damage by mammals in north temperate forests. J. For. Res. 1992, 65, 145–169. [Google Scholar] [CrossRef]
  23. Reimoser, F.; Putman, R.J. Impact of large ungulates on agriculture, forestry and conservation habitats in Europe. In Ungulate Management in Europe: Problems and Practices; Putman, R.J., Apollonio, M., Andersen, R., Eds.; Cambridge University Press: Cambridge, UK, 2011; pp. 144–191. [Google Scholar]
  24. Borkowski, J.; Dobrowolska, D.; Dąbrowski, W. Young conifer stands form a deer browsing refuge for an oak admixture: Silvicultural implications for forest regeneration under herbivore pressure. Eur. J. Forest. Res. 2017, 136, 787–800. [Google Scholar] [CrossRef]
  25. Götmark, F.; Berglund, Å.; Wiklander, K. Browsing damage on broadleaved trees in semi-natural temperate forest in Sweden, with a focus on oak regeneration. Scand. J. For. Res. 2005, 20, 223–234. [Google Scholar] [CrossRef]
  26. Bergquist, J.; Löf, M.; Örlander, G. Effects of roe deer browsing and site preparation on performance of planted broadleaved and conifer seedlings when using temporary fences. Scand. J. For. Res. 2009, 24, 308–317. [Google Scholar] [CrossRef]
  27. Lithuanian Statistical Yearbook of Forestry. Ministry of Environment, State Forest Service. 2019. Available online: https://amvmt.lrv.lt/lt/atviri-duomenys-1/misku-statistikos-leidiniai/misku-ukio-statistika/2021-m-1 (accessed on 14 August 2023).
  28. Shigo, A.L. Modern arboriculture: A systems approach to the care of trees and their associates. In Shigo and Trees; Associates Durham: Durham, UK, 1991; 424p. [Google Scholar]
  29. Wargo, P.M. Consequences of environmental stress on oak: Predisposition to pathogens. Ann. For. Sci. 1996, 53, 359–368. [Google Scholar] [CrossRef]
  30. Griffith, G.S.; Boddy, L. Fungal decomposition of attached angiosperm twigs: III. Effect of water potential and temperature on fungal growth, survival and decay of wood. New Phyt. 1991, 117, 259–269. Available online: https://www.jstor.org/stable/2556881 (accessed on 27 July 2023). [CrossRef]
  31. Sieber, T.N. Endophytic fungi in forest trees: Are they mutualists? Fungal Biol. Rev. 2007, 21, 75–89. [Google Scholar] [CrossRef]
  32. Parfitt, D.; Hunt, J.; Dockrell, D.; Rogers, H.J.; Boddy, L. Do all trees carry the seeds of their own destruction? PCR reveals numerous wood decay fungi latently present in sapwood of a wide range of angiosperm trees. Fungal Ecol. 2010, 3, 338–346. [Google Scholar] [CrossRef]
  33. Cline, L.C.; Schilling, J.S.; Menke, J.; Groenhof, E.; Kennedy, P.G. Ecological and functional effects of fungal endophytes on wood decomposition. Funct. Ecol. 2018, 32, 181–191. [Google Scholar] [CrossRef]
  34. Kogel, K.H.; Franken, P.; Hückelhoven, R. Endophyte or parasite–what decides? Curr. Opin. 2006, 9, 358–363. [Google Scholar] [CrossRef] [PubMed]
  35. Agostinelli, M.; Cleary, M.; Martin, J.A.; Albrectsen, B.R.; Witzell, J. Pedunculate oaks (Quercus robur L.) differing in vitality as reservoirs for fungal biodiversity. Front. Microb. 2018, 9, 1758. [Google Scholar] [CrossRef]
  36. Boddy, L. Fungal community ecology and wood decomposition processes in angiosperms: From standing tree to complete decay of coarse woody debris. Ecol. Bull. 2001, 49, 43–56. Available online: http://www.jstor.org/stable/20113263 (accessed on 27 July 2023).
  37. Bässler, C.; Müller, J.; Dziock, F.; Brandl, R. Effects of resource availability and climate on the diversity of wood-decaying fungi. J. Ecol. 2010, 98, 822–832. [Google Scholar] [CrossRef]
  38. Shortle, W.C.; Dudzik, K.R. Wood Decay in Living and Dead Trees: A Pictorial Overview; US Department of Agriculture, Forest Service, Northern Research Station: Newtown Square, PA, USA, 2012; pp. 10–21.
  39. Marčiulynas, A.; Lynikienė, J.; Marčiulynienė, D.; Gedminas, A.; Menkis, A. Seasonal and Site-Specific Patterns of Airborne Fungal Diversity Revealed Using Passive Spore Traps and High-Throughput DNA Sequencing. Diversity 2023, 15, 539. [Google Scholar] [CrossRef]
  40. Xu, Z.; Leininger, T.D.; Lee, A.W.; Tainter, F.H. Chemical properties associated with bacterial wetwood in red oaks. Wood Fiber Sci. 2001, 33, 76–83. [Google Scholar]
  41. van der Wal, A.; Ottosson, E.; De Boer, W. Neglected role of fungal community composition in explaining variation in wood decay rates. Ecology 2015, 96, 124–133. [Google Scholar] [CrossRef]
  42. Rayner, A.D.M.; Boddy, L. Fungal communities in the decay of wood. Adv. Microb. Ecol. 1988, 10, 115–166. [Google Scholar] [CrossRef]
  43. Marčiulynas, A.; Marčiulynienė, D.; Lynikienė, J.; Bakys, R.; Menkis, A. Fungal Communities in Leaves and Roots of Healthy-Looking and Diseased Ulmus glabra. Microorganisms 2022, 10, 2228. [Google Scholar] [CrossRef] [PubMed]
  44. Ihrmark, K.; Bodeker, I.T.M.; Cruz-Martinez, K.; Friberg, H.; Kubartova, A.; Schenck, J.; Strid, Y.; Stenlid, J.; Brandstrom-Durling, M.; Clemmensen, K.E.; et al. New primers to amplify the fungal ITS2—Evaluation by 454-sequencing of artificial and natural communities. FEMS Microbiol. Ecol. 2012, 82, 666–677. [Google Scholar] [CrossRef] [PubMed]
  45. White, T.J.; Bruns, T.; Lee, S.; Taylor, J. Amplification and direct sequencing of fungal ribosomal RNA genes for phylogenetics. In PCR Protocols: A Guide to Methods and Applications; Innis, M.A., Gelfand, D.H., Sninsky, J.J., White, T.J., Eds.; Academic Press: San Diego, CA, USA, 1990; pp. 315–322. [Google Scholar]
  46. Kubartová, A.; Ottosson, E.; Dahlberg, A.; Stenlid, J. Patterns of fungal communities among and within decaying logs, revealed by 454 sequencing. Molec. Ecol. 2012, 21, 4514–4532. [Google Scholar] [CrossRef] [PubMed]
  47. Nguyen, N.H.; Song, Z.; Bates, S.T.; Branco, S.; Tedersoo, L.; Menke, J.; Schilling, J.S.; Kennedy, P.G. FUNGuild: An open annotation tool for parsing fungal community datasets by ecological guild. Fung. Ecol. 2016, 20, 241–248. [Google Scholar] [CrossRef]
  48. Magurran, A.E. Ecological Diversity and Its Measurement; Princeton University Press: Princeton, NJ, USA, 1988; 192p. [Google Scholar]
  49. Shannon, C.E. A mathematical theory of communication. Bell Syst. Tech. J. 1948, 27, 379–423. [Google Scholar] [CrossRef]
  50. Oksanen, J.; Blanchet, F.G.; Kindt, R.; Legendre, P.; Minchin, P.R.; O’Hara, R.B.; Simpson, G.L.; Solymos, P.; Stevens, M.H.H.; Wagner, H. Community Ecology Package. R Package Version 2. 2013. Available online: https://cran.r-project.org/web/packages/vegan/index.html (accessed on 27 July 2023).
  51. R Core Team. R: A Language and Environment for Statistical Computing; R Foundation for Statistical Computing: Vienna, Austria, 2021. [Google Scholar]
  52. Blanchette, R.A. Anatomical responses of xylem to injury and invasion by fungi. In Defense Mechanisms of Woody Plants against Fungi; Springer: Berlin/Heidelberg, Germany, 1992; pp. 77–95. [Google Scholar]
  53. Boddy, L. Microenvironmental aspects of xylem defenses to wood decay fungi. In Defense Mechanisms of Woody Plants against Fungi; Springer: Berlin/Heidelberg, Germany, 1992; pp. 96–132. [Google Scholar]
  54. Boddy, L. Effect of temperature and water potential on growth rate of wood-rotting basidiomycetes. Trans. Br. Mycol. Soc. 1983, 80, 141–149. [Google Scholar] [CrossRef]
  55. Ryvarden, L.; Gilbertson, R.L. European Polypores: Part 1: Abortiporus-Lindtneria; Department of Botany, University of Oslo: Blindern/Oslo, Norway, 1993; 387p. [Google Scholar]
  56. Vasaitis, R.; Menkis, A.; Lim, Y.W.; Seok, S.; Tomsovsky, M.; Jankovsky, L.; Lygis, V.; Slippers, B.; Stenlid, J. Genetic variation and relationships in Laetiporus sulphureus s. lat., as determined by ITS rDNA sequences and in vitro growth rate. Mycol. Res. 2009, 113, 326–336. [Google Scholar] [CrossRef]
  57. Elisashvili, V.; Penninckx, M.; Kachlishvili, E.; Asatiani, M.; Kvesitadze, G. Use of Pleurotus dryinus for lignocellulolytic enzymes production in submerged fermentation of mandarin peels and tree leaves. Enzyme Microb. Technol. 2006, 38, 998–1004. [Google Scholar] [CrossRef]
  58. Swift, M.J. Growth of Stereum hirsutum during the long-term decomposition of oak branch-wood. Soil Biol. Biochem. 1978, 10, 335–337. [Google Scholar] [CrossRef]
  59. Vasiliauskas, R.; Stenlid, J. Fungi inhabiting stems of Picea abies in a managed stand in Lithuania. For. Ecol. Manag. 1998, 109, 119–126. [Google Scholar] [CrossRef]
  60. Mirić, M.; Stefanović, M. The spread of four wood-decaying fungi through artificially infected healthy trees of pedunculate oak (Quercus robur L.) in vivo. Šumarstvo 2018, 1, 79–90. [Google Scholar]
  61. Brasier, C.M.; Kirk, S.A.; Mehrotra, M.D. Ophiostoma himalulmi sp. nov., a new species of Dutch elm disease fungus endemic to the Himalayas. Mycol. Res. 1995, 99, 205–215. [Google Scholar] [CrossRef]
  62. Harrington, T.C.; McNew, D.; Steimel, J.; Hofstra, D.; Farrell, R. Phylogeny and taxonomy of the Ophiostoma piceae complex and the Dutch elm disease fungi. Mycologia 2001, 93, 111–136. [Google Scholar] [CrossRef]
  63. de Amorim, M.R.; Hilario, F.; Dos Santos, F.M., Jr.; Batista, J.M., Jr.; Bauab, T.M.; Araujo, A.R.; Carlos, I.Z.; Vilegas, W.; Dos Santos, L.C. New benzaldehyde and benzopyran compounds from the endophytic fungus Paraphaeosphaeria sp. F03 and their antimicrobial and cytotoxic activities. Planta Med. 2019, 85, 957–964. [Google Scholar] [CrossRef] [PubMed]
  64. Liu, T.; Greenslade, A.; Yang, S. Levels of rhizome endophytic fungi fluctuatein Paris polyphylla var. yunnanensis as plants age. Plant Divers. 2017, 39, 60–64. [Google Scholar] [CrossRef]
  65. Verkley, G.J.M.; Dukik, K.; Renfurm, R.; Göker, M.; Stielow, J.B. Novel genera and species of coniothyrium-like fungi in Montagnulaceae (Ascomycota). Persoonia 2014, 32, 25–51. [Google Scholar] [CrossRef]
  66. Rooney, S.N.; Eskalen, A.; Gubler, W.D. Recovery of Phaeomoniella chlamydospora and Phaeoacremonium in-flatipes from soil and grapevine tissues. Phytopathol. Mediterr. 2001, 40, 351–356. Available online: https://www.jstor.org/stable/44981641 (accessed on 7 August 2023).
  67. Ariyawansa, H.A.; Phukhamsakda, C.; Thambugala, K.M.; Bulgakov, T.S.; Wanasinghe, D.N.; Perera, R.H.; Mapook, A.; Camporesi, E.; Kang, J.-C.; Jones, E.B.G.; et al. Revision and phylogeny of Leptosphaeriaceae. Fun. Divers. 2015, 74, 19–51. [Google Scholar] [CrossRef]
  68. El-Demerdash, A.; El-Demerdash, A. Chemical diversity and biological activities of Phaeosphaeria fungi genus: A systematic review. J. Fungi 2018, 4, 130. [Google Scholar] [CrossRef]
  69. Khaledi, E.; Moghadam, J.N.; Abdollahzadeh, J.; Amini, J. Fungi associated with grapevine trunk diseases (GTDs) with emphasize on pestalotioid species in Kurdistan Province, Iran. Res. Sq. 2021. [Google Scholar] [CrossRef]
  70. Niemelä, T.; Renvall, P.; Pentillä, R. Interactions of fungi at late stages of wood decomposition. Ann. Bot. Fenn. 1995, 32, 141–152. Available online: https://www.jstor.org/stable/23726315 (accessed on 8 August 2023).
  71. Heilmann-Clausen, J.; Boddy, L. Inhibition and stimulation effects in communities of wood decay fungi: Exudates from colonized wood influence growth by other species. Microb. Ecol. 2005, 49, 399–406. [Google Scholar] [CrossRef] [PubMed]
  72. Liew, E.C.Y.; Aptroot, A.; Hyde, K.D. An evaluation of the monophyly of Massarina based on ribosomal DNA sequences. Mycologia 2002, 94, 803–813. [Google Scholar] [CrossRef] [PubMed]
  73. Chen, J. Fungal Community Survey of Fraxinus excelsior in New Zealand. Dissertation’s Thesis, Swedish University of Agricultural Sciences, Uppsala, Switzerland, 2012. [Google Scholar]
  74. Davydenko, K.; Vasaitis, R.; Stenlid, J.; Menkis, A. Fungi in foliage and shoots of Fraxinus excelsior in eastern Ukraine: A first report on Hymenoscyphus pseudoalbidus. For. Pathol. 2013, 43, 462–467. [Google Scholar] [CrossRef]
  75. Gherghel, F.; Fussi, B.; Donges, K.; Haustein, M.; Jakob, K.M.; Müller, K.; Piškur, B.; Hauptman, T.; Lenz, H.D.; Konnert, M.; et al. The development of a species-specific test to detect Hymenoscyphus pseudoalbidus in ash tissues. For. Pathol. 2014, 44, 137–144. [Google Scholar] [CrossRef]
  76. Rappaz, F. Taxonomie et nomenclature des Diatrypacees a asques octospores. Mycol. Helv. 1987, 2, 285–648. Available online: https://folia.unifr.ch/unifr/documents/302452 (accessed on 7 August 2023).
  77. Kepley, J.B.; Reeves, F.B.; Jacobi, W.R.; Adams, G.C. Species associated with Cytospora canker on Populus tremuloides. Mycotaxon 2015, 130, 783–805. [Google Scholar] [CrossRef]
  78. Wang, Y.L.; Lu, Q.; Decock, C.; Li, Y.X.; Zhang, X.Y. Cytospora species from Populus and Salix in China with C. davidiana sp. nov. Fungal Biol. 2015, 119, 420–432. [Google Scholar] [CrossRef]
  79. Jacobs, K.; Plaschkies, K.; Scheiding, W.; Weiß, B.; Melcher, E.; Conti, E.; Le Bayon, I. Natural durability of important European wood species against wood decay fungi. Part 2: Field tests and fungal community. Int. Biodeterior. Biodegrad. 2019, 137, 118–126. [Google Scholar] [CrossRef]
  80. Anagnostakis, S.L. Diaporthe eres (Phomopsis oblonga) as a pathogen of butternut (Juglans cinerea) in Connecticut. Plant Dis. 2007, 91, 1198. [Google Scholar] [CrossRef]
  81. Gao, Y.H.; Su, Y.Y.; Sun, W.; Cai, L. Diaporthe species occurring on Lithocarpus glabra in China, with descriptions of five new species. Fungal Biol. 2015, 119, 295–309. [Google Scholar] [CrossRef] [PubMed]
  82. Dissanayake, A.J.; Chen, Y.Y.; Liu, J.K. Unravelling Diaporthe species associated with woody hosts from Karst Formations (Guizhou) in China. J. Fungi 2020, 6, 251. [Google Scholar] [CrossRef]
  83. Moricca, S.; Ginetti, B.; Ragazzi, A. Species- and organ-specificity in endophytes colonizing healthy and declining Mediterranean oaks. Phytopathol. Mediterr. 2012, 51, 587–598. Available online: http://www.jstor.org/stable/43872345 (accessed on 10 August 2023).
  84. Coates, D.; Rayner, A.D.M. Fungal population and community development in cut beech logs: I. Establishment via the aerial cut surface. New Phytol. 1985, 101, 153–171. [Google Scholar] [CrossRef] [PubMed]
  85. Fukasawa, Y.; Osono, T.; Takeda, H. Wood decomposing abilities of diverse lignicolous fungi on nondecayed and decayed beech wood. Mycologia 2011, 103, 474–482. [Google Scholar] [CrossRef] [PubMed]
  86. Iršėnaitė, R.; Kutorga, E. Diversity of fungi on decaying common oak coarse woody debris. Ekologija 2006, 4, 22–30. Available online: http://elibrary.lt/resursai/LMA/Ekologija/Eko64/Eko64_07.pdf (accessed on 10 August 2023).
  87. Županić, M.; Matošević, D.; Pernek, M.; Diminić, D. Lignicolous fungi on Pedunculate oak in lowland forests of Central Croatia. Period. Biol. 2009, 111, 397–403. Available online: https://hrcak.srce.hr/47874 (accessed on 11 August 2023).
  88. Mallette, N.D.; Pankrantz, E.M.; Busse, S.; Strobel, G.A.; Carlson, R.P.; Peyton, B.M. Evaluation of cellulose as a substrate for hydrocarbon fuel production by Ascocoryne sarcoides (NRRL 50072). J. Sustain. Bioenergy Syst. 2014, 4, 33–49. [Google Scholar] [CrossRef]
  89. Lindhe, A.; Åsenblad, N.; Toresson, H.G. Cut logs and high stumps of spruce, birch, aspen and oak–nine years of saproxylic fungi succession. Biol. Conserv. 2004, 119, 443–454. [Google Scholar] [CrossRef]
  90. Bunyard, B.A.; Wang, Z.; Malloch, D.; Clayden, S.; Voitk, A. New North American Records for Ascocoryne turficola (Ascomycota: Helotiales). Fungi 2008, 1, 23–31. [Google Scholar]
Figure 1. Wounds of different ages made by deer animals on the stems of Quercus robur: (A) a 20-year-old oak with an eight-year-old wound with exposed wood; (B) a 50-year-old oak with a more than 30-year-old wound, which was closed.
Figure 1. Wounds of different ages made by deer animals on the stems of Quercus robur: (A) a 20-year-old oak with an eight-year-old wound with exposed wood; (B) a 50-year-old oak with a more than 30-year-old wound, which was closed.
Forests 14 02077 g001
Figure 2. Map of Lithuania showing the distribution of common oak (Quercus robur) in forest stands (in green). The intensity of green color shows the percentage of Q. robur in the composition of forest stands. Sampling sites are denoted by K—Kėdainiai, P—Prienai, and T—Telšiai. The map was reproduced with permission from FORESTGEN, www.forestgen.mi.lt (accessed on 13 July 2023).
Figure 2. Map of Lithuania showing the distribution of common oak (Quercus robur) in forest stands (in green). The intensity of green color shows the percentage of Q. robur in the composition of forest stands. Sampling sites are denoted by K—Kėdainiai, P—Prienai, and T—Telšiai. The map was reproduced with permission from FORESTGEN, www.forestgen.mi.lt (accessed on 13 July 2023).
Forests 14 02077 g002
Figure 3. Rarefaction curves showing the relationship between the accumulated number of fungal OTUs and the number of ITS2 rDNA sequences from oak stem wounds of different ages (10 y, 20 y, 30 y, 40 y, and 50 y). Data from different sites are combined.
Figure 3. Rarefaction curves showing the relationship between the accumulated number of fungal OTUs and the number of ITS2 rDNA sequences from oak stem wounds of different ages (10 y, 20 y, 30 y, 40 y, and 50 y). Data from different sites are combined.
Forests 14 02077 g003
Figure 4. The richness of fungal OTUs in oak wounds of different ages (10 y, 20 y, 30 y, 40 y, and 50 y) at three different sites (Prienai, Telšiai, and Kėdainiai) in Lithuania. Within each site, columns followed by the same letter do not differ significantly at p > 0.05.
Figure 4. The richness of fungal OTUs in oak wounds of different ages (10 y, 20 y, 30 y, 40 y, and 50 y) at three different sites (Prienai, Telšiai, and Kėdainiai) in Lithuania. Within each site, columns followed by the same letter do not differ significantly at p > 0.05.
Forests 14 02077 g004
Figure 5. Relative abundance (%) of most common fungal OTUs in Quercus robur wounds of different ages (10 y–50 y). Others include fungal OTUs with a relative abundance of <0.7% for Basidiomycota and of <1.5% for Ascomycota.
Figure 5. Relative abundance (%) of most common fungal OTUs in Quercus robur wounds of different ages (10 y–50 y). Others include fungal OTUs with a relative abundance of <0.7% for Basidiomycota and of <1.5% for Ascomycota.
Forests 14 02077 g005
Figure 6. Venn diagrams showing the diversity and overlap of fungal OTUs in oak wounds of different ages (10 y–50 y). Data from different sites are combined.
Figure 6. Venn diagrams showing the diversity and overlap of fungal OTUs in oak wounds of different ages (10 y–50 y). Data from different sites are combined.
Forests 14 02077 g006
Figure 7. Relative abundance (%) of fungal functional groups in Q. robur stem wounds of different ages, estimated based on fungal sequences. Others represent fungi, which are not associated with plants (e.g., animal pathogens).
Figure 7. Relative abundance (%) of fungal functional groups in Q. robur stem wounds of different ages, estimated based on fungal sequences. Others represent fungi, which are not associated with plants (e.g., animal pathogens).
Forests 14 02077 g007
Figure 8. Non-metric multidimensional scaling (NMDS) of fungal communities associated with Quercus robur stem wounds of different ages (10, 20, 30, 40, or 50 years old).
Figure 8. Non-metric multidimensional scaling (NMDS) of fungal communities associated with Quercus robur stem wounds of different ages (10, 20, 30, 40, or 50 years old).
Forests 14 02077 g008
Table 1. Generated high-quality fungal sequences and detected diversity of fungal OTUs in oak wounds of different ages in three sampling sites in Lithuania.
Table 1. Generated high-quality fungal sequences and detected diversity of fungal OTUs in oak wounds of different ages in three sampling sites in Lithuania.
SiteSequences/fungal OTUs in oak wounds of different ages
10 y20 y30 y40 y50 y
Prienai118/31738/3168/244/31261/21
564/154/41714/17179/201092/24
1759/258/8208/2050/18128/20
1542/272633/1616/1137/275154/18
4731/2487/2617/1226/13654/18
Prienai total8714/713450/592023/57296/548289/73
Telšiai88/131/1400/1064/61283/31
2/2621/25202/1436/64008/60
243/142179/1022/4477/655/9
20/91011/4024/8352/8139/6
6832/463/3536/173162/9-/-
Telšiai total7185/723815/651184/444091/265485/92
Kėdainiai8565/622338/393514/8124/483300/8
8527/6855/2460/14372/12161/8
10,894/7816/111571/597848/261067/8
8965/58662/263789/515770/745105/36
9556/5991/233804/33488/33380/16
Kėdainiai total46,507/1523162/8312,738/12814,602/14110,013/56
All total62,406/21310,427/15915,945/17918,989/19323,787/179
Shannon diversity index (H)2.742.892.902.512.93
Table 2. Occurrence and relative abundance of the 30 most common fungal OTUs (shown as a proportion of all high-quality fungal sequences) in Quercus robur stem wounds of different ages (10–50 years old). Data from different sites are combined.
Table 2. Occurrence and relative abundance of the 30 most common fungal OTUs (shown as a proportion of all high-quality fungal sequences) in Quercus robur stem wounds of different ages (10–50 years old). Data from different sites are combined.
Phylum *OTUGenbank ReferenceSimilarity, (%)Age of Wound Total, %
10 y, %20 y, %30 y, %40 y, %50 y, %
AAposphaeria corallinoluteaMT177916245/245 (100)23.741.38-0.04-11.38
ASclerostagonospora cycadisKR611890248/248 (100)13.290.330.040.01-6.34
BLaetiporus sulphureusMH321898296/296 (100)--22.040.0116.625.68
ACadophora malorumMT561395241/241 (100)0.0122.288.131.6010.254.84
ADiatrype stigmaKT004563251/252 (99)9.570.130.01--4.55
AFimetariella rabenhorstiiMN547388248/250 (99)-0.090.4730.69-4.49
ANeosetophoma italicaLC206635249/249 (100)8.980.07-0.01-4.27
APhaeostalagmus cyclosporusON989633233/238 (98)--0.780.3920.983.94
ANeocucurbitaria quercinaOP896095249/249 (100)4.374.630.035.400.163.25
BMycena galericulataKJ705178314/314 (100)--2.0516.540.782.78
ALeptosillia macrosporaNR164064265/265 (100)0.140.65-18.320.012.76
ACytospora predappioensisMT595364260/260 (100)5.340.090.030.01-2.54
ACladosporium cladosporioidesOP963820243/243 (100)5.170.190.090.170.032.51
AParaphaeosphaeria sporulosaKX664338244/245 (99)5.170.47---2.49
ALophiostoma corticolaKU712227244/244 (100)3.760.290.030.04-1.81
ANeoleptosphaeria rubefaciensKT804116250/250 (100)3.350.010.31--1.63
ASporocadus trimorphusMT223846248/248 (100)3.350.01---1.59
AAscocoryne sarcoidesMH857562240/240 (100)0.0030.700.240.068.041.55
AMenispora sp. 5631_29MH859920222/233 (95)--12.200.06-1.49
APhialocephala compactaMH862480234/239 (98)0.0037.364.951.380.351.45
ADiaporthe eresMT111115255/255 (100)0.3115.63---1.39
APhialemonium inflatumMT221573242/246 (98)-0.0210.570.040.0041.29
AUnidentified sp. 5631_30MN096587254/265 (96)-15.380.03--1.22
APhlyctema vagabundaMK174720239/239 (100)2.500.030.19--1.21
BCylindrobasidium evolvensMN947592387/387 (100)2.24----1.06
AHyaloscypha fuckeliiMT231692238/238 (100)---0.015.761.04
BArmillaria gallicaKY474051476/477 (99)--0.01-5.250.95
BPhlebia acerinaMN945144289/289 (100)--7.46--0.90
BPleurotus dryinusMK169240284/284 (100)-0.40-0.243.660.73
ALophiostoma sp. 5631_40MH178565240/252 (95)----3.960.72
Total 91.2870.1469.6475.0275.8481.84
* A—Ascomycota, B—Basidiomycota.
Table 3. The Sørensen similarity index (β-diversity) of fungal communities among wounds of different ages (10–50 years old). Data from different sites are combined.
Table 3. The Sørensen similarity index (β-diversity) of fungal communities among wounds of different ages (10–50 years old). Data from different sites are combined.
Wound Age10 y20 y30 y40 y50 y
10 y-0.320.260.350.19
20 y0.32-0.430.420.36
30 y0.260.43-0.460.40
40 y0.350.420.46-0.35
50 y0.190.360.400.35-
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Marčiulynas, A.; Sirgedaitė-Šėžienė, V.; Menkis, A. Fungi Inhabiting Stem Wounds of Quercus robur following Bark Stripping by Deer Animals. Forests 2023, 14, 2077. https://doi.org/10.3390/f14102077

AMA Style

Marčiulynas A, Sirgedaitė-Šėžienė V, Menkis A. Fungi Inhabiting Stem Wounds of Quercus robur following Bark Stripping by Deer Animals. Forests. 2023; 14(10):2077. https://doi.org/10.3390/f14102077

Chicago/Turabian Style

Marčiulynas, Adas, Vaida Sirgedaitė-Šėžienė, and Audrius Menkis. 2023. "Fungi Inhabiting Stem Wounds of Quercus robur following Bark Stripping by Deer Animals" Forests 14, no. 10: 2077. https://doi.org/10.3390/f14102077

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop