INFORMATION ABOUT PROJECT,
SUPPORTED BY RUSSIAN SCIENCE FOUNDATION

The information is prepared on the basis of data from the information-analytical system RSF, informative part is represented in the author's edition. All rights belong to the authors, the use or reprinting of materials is permitted only with the prior consent of the authors.

 

COMMON PART


Project Number22-24-01063

Project titleMetagenomic analysis of xylobiont fungi and bacteria in deadwoods of deciduous and coniferous trees at different stages of decomposition after a catastrophic windthrow in a mesic broad-leaved forest

Project LeadKhanina Larisa

AffiliationKeldysh Institute of Applied Mathematics of Russian Academy of Sciences,

Implementation period 2022 - 2023 

Research area 04 - BIOLOGY AND LIFE SCIENCES, 04-107 - Ecology of biological systems

Keywordsbiodiversity, DNA sequencing, mycobiota, fungi, macromycetes, bacteria, DNA, metagenomic analysis, coarse woody debris, dead wood, physicochemical properties of wood, ecology of microorganisms, broadleaved deciduous forests, windthrow


 

PROJECT CONTENT


Annotation
The species diversity of xylobiont fungi and bacteria in a catastrophic windthrow of a multi-species broad-leaved forest in the Kaluzhskie Zaseki State Nature Reserve will be assessed by means of metagenomic analysis. Forests in the Reserve are unique in their multi-species composition, high tree age and relatively large area; they are "reference areas" for both forest biodiversity protection and ecosystem-based forest management in the temperate forest region. Coarse woody debris (deadwood) plays a major role in nutrient cycling as a temporary storage stock of carbon and macro- and micronutrients, which become available to biota during the deadwood decomposition. Understanding the factors affecting wood decomposition is necessary to assess the carbon sequestration capacity of forests, which is a key issue addressing climate change. Fungi and bacteria are the main destructors of deadwood; xylotrophic fungal communities being highly diverse and including rare and endangered taxa. Despite the high importance of research on wood decomposition processes, their detailed analysis has only recently begun. It becomes possible due to the emergence of the newest methods for analyzing natural ecosystems, methods for assessing their diversity and functioning. One of these newest methods is metagenomic research, which allows analysis of fungal diversity not only by fruiting bodies, but also by their dark diversity and by the bacterial community that develops in dead wood. The novelty and originality of the research is determined both by the integrated approach (ecology of biosystem, mycology and structural genomics) and by the uniquely wide range of research objects. For the first time, metagenomics methods will be applied to identify the diversity of xylobiont fungi and bacteria on deadwood samples of eight tree species (Quercus robur, Fraxinus excelsior, Tilia cordata, Ulmus glabra, Acer platanoides, Populus tremula, Betula pendula, and Picea abies) at five decay stages; all the deadwood was formed due to the natural disturbance by catastrophic windthrow happened in 2006. A significant objective of the project is DNA barcoding of macromycetes based on sequencing of fruit bodies collected in habitats where sampling for metagenomic analysis will be carried out and for which there is no information in reference databases on nucleotide sequences. Sanger sequencing of freshly collected fruit body samples of fungi that develop directly in the habitats under study will allow the generation of an array of reference nucleotide sequences and reduce the number of 'nameless' OTUs when interpreting metagenomic sequencing results. The project will conclude with an assessment of fungal and bacterial diversity and an analysis of the relationships between the diversity and composition of xylobionts and the taxonomic identity, decomposition stage of the wood samples, and their physical and chemical characteristics as well.

Expected results
For the first time, the diversity of fungi and bacteria developing on fallen logs of eight tree species (Quercus robur, Fraxinus excelsior, Tilia cordata, Ulmus glabra, Acer platanoides, Populus tremula, Betula pendula, and Picea abies) at five decay stages and located in natural conditions of a multi-species broadleaved forest will be identified by performing metagenomic analysis. The results of the fungi occurrences will be available through international repositories. Reference DNA sequences will be obtained from fruit body samples of fungi for which no data are available in nucleotide databases. Information on the new nucleotide sequences will be deposited in the international GenBank database. The diversity of fungi and bacteria will be assessed; the relationships between the diversity and composition of fungal and bacterial communities with the taxonomic identity of the woody samples and their physico-chemical characteristics, including the stage of deadwood decomposition, will be analyzed. The expected results are highly significant for international and Russian science; they are important for solving the urgent problem of climate change mitigation on the Earth.


 

REPORTS


Annotation of the results obtained in 2023
In the second year of the project, we deposited in the international GenBank database information on 20 new nucleotide sequences of fungi whose fruiting bodies were collected directly in the sites where sampling for metagenomic analysis was performed. The complete nucleotide sequences of ITS1-5.8S-ITS2 nrDNA obtained and deposited in NCBI GenBank are represented mainly by taxa about which there is no information in reference databases. In particular, information on the primary structure of the ITS region of nrDNA was obtained for the first time for corticioid fungi from the genera Peniophorella (Hymenochaetales), Hyphoderma, Phlebia (Polyporales), and Amylocorticium (Amylocorticiales). In addition, nucleotide sequences were obtained from authentic samples for little-known and rarely found in nature representatives of the genera of xylotrophic basidiomycetes Athelia (Atheliales), Leucogyrophana (Boletales), Kneiffiella (Hymenochaetales), Yuchengia (Polyporales). Information on 146 fungal findings from a mass windthrow area in the mixed broadleaved forest was published in GBIF. A total of 4717 different ITS1 sequence variants of fungi and 92262 sequence variants of the V4 regions of the 16S rRNA gene of bacteria were obtained from deadwood metagenomes. After downsampling, about 2500 and about 22500 unique sequences remained in the fungal and bacterial community, respectively (with only 643 and 665 of these representing one per cent or more of all community ASVs in at least one sample). All tree species examined showed an increase in the alpha diversity of the fungal and bacterial communities between the initial and final stages of decomposition, with the bacterial community being more diverse than the fungal community at all stages of decomposition. In the first stages, the fungal community consisted on average of about 25 species with strongly dominant representatives. The exception was Picea abies, on which more than 60 species were found. The bacterial community of each tree species consisted of over 200 species, except for Acer platanoides and Betula pendula, in which 45 and 86 species were observed, respectively. By the fourth and fifth stages of decomposition, an average of over 100 fungal species and between approximately 400 and over 600 bacterial species were observed in all tree species. Some tree species showed a decline in the diversity of the investigated destructors at the second or third stage of decomposition. With the overall high noisiness of the data, the concentration of amplicons from the bacterial community was more often higher at the last stage of decomposition than at the first, whereas the fungal community often showed a peak at earlier stages (1-3) and a decline at the last. According to our data, the highest concentration of amplicons from fungi and bacteria was characteristic for Tilia cordata and the lowest for Quercus rubra. In oak wood, bacteria appeared in noticeable numbers only by the third stage of decomposition, which probably correlates with its resistance to decomposition and, consequently, its active use in the economy. Interestingly, in some cases (e.g. at the first stage of decomposition in Fraxinus excelsior and the second in Picea) fungi significantly prevailed over bacteria, which may be due to the release of bactericidal substances by them, and the number of bacteria increased with the transition to later stages of decomposition. The fungal community was characterised by very high beta diversity (mean Bray-Curtis distance 0.96 ± 0.08), with the distance between samples of different tree species of the same decomposition stage being slightly smaller than the distance between samples of the same species but different stages, resulting in a more pronounced clustering of samples by decomposition stage than by tree species. During wood decomposition, the fungal community for some tree species changed dramatically, showing a change of dominant representatives at the phylum level. Of the 47 classes found in the fungal community, 13 were observed in at least one sample of each species and each stage of decomposition. At the same time, most fungal genera were observed only on samples of one species of one stage, suggesting that there is no single common pathway of wood decomposition by the fungal community. In the bacterial community, 13 out of 104 classes and 146 out of 781 genera were observed in samples of all types, and, as in the fungal community, frequent classes included dominant representatives. When beta diversity was assessed using the Bray-Curtis distance at the ASV level, bacterial communities also varied greatly between samples. In terms of similarity of bacterial communities (as with fungi, but more so for bacteria), samples clustered more by stage than by tree species. Based on the records of fruiting bodies of aphyllophoroid fungi, 127 species of basidial fungi were identified on dead tree trunks of 8 tree species. The fungal species belonged to 78 genera, 50 families, and 14 orders. The maximum number of fungal species (46) occurred on Picea and minimum (11) on Quercus while 38, 24, 22, 19, 18, and 13 were surveyed on Acer, Fraxinus, Populus tremula, Ulmus glabra, and Tilia, respectively. White rot fungi dominated (74%), while brown rot fungi accounted for 19%; morphologically, corticoid and popypore comprised 55 and 32%, respectively, and clavarioid and hydnoid together accounted for 12% of sporocarps. The mean number of fungal species was 1.8 ± 0.06 in the 1-m plot and 10.4 ± 0.98 in the trunk. High gamma diversity was supported by high beta diversity of the fungal community: mean Bray-Curtis distance was 0.95 ± 0.001 and 0.92 ± 0.006 for plots and trunks, respectively. Trunks of different tree species did not differ in the number of fungal species found on them, while groups of tree species differed significantly: the number of fungal species found on ring-porous angiosperms (10.6 ± 1.2) was significantly lower than those found on the gymnosperm Picea (16.8 ± 2.6); Kruskal–Wallis test χ2 = 9.48, P = 0.009. The composition of basidiomycetes on lying trunks of various tree species differed rather well, as shown by the results of non-metric multidimensional scaling (NMDS) and permutational multivariate analysis (PERMANOVA). Tree species identity determined 45% of the variation in basidiomycete species composition (P = 0.0001); together with mean stem decomposition stage (R2 = 5%, P = 0.05), dbh and number of non-zero sites, the model determined 58% of the variation in fungal species composition. Additive partitioning of biodiversity showed that the contributions of all diversity components were significant, but the maximum contribution (81%) to the total gamma diversity of basidiomycetes was made by compositional diversity among tree species, which was significantly higher than that expected under the null model. Thus, we obtained overall consistent estimates of the diversity of fungal communities identified from fruiting bodies and metagenomic analyses. The two analyses showed an increase in community diversity with increasing decomposition stage and the importance of the contribution of woody substrate diversity to the total gamma diversity of the fungal population of the analyzed site. Metagenomic analysis allowed a detailed assessment of the quantitative and qualitative change in the species diversity of destructor communities (bacteria and fungi) during wood decomposition in different tree species; specific features of the decomposition process for each of the eight tree species analyzed were suggested. The results of the research were presented at 4 scientific conferences. Nine scientific articles were published, among them one in Q1 journal (IF 4.384); one in a journal registered in Scopus database; also 4 articles in conference proceedings and 2 datasets were published. Web pages of the project were updated https://www.impb.ru/index.php?id=div/lce/22-24-01063/index and https://www.impb.ru/index.php?id=div/lce/22-24-01063/index&lang=eng

 

Publications

1. Khanina L., Bobrovsky M., Smirnov V., Romanov M. Wood decomposition, carbon, nitrogen, and pH values in logs of 8 tree species 14 and 15 years after a catastrophic windthrow in a mesic broad-leaved forest in the East European plain Forest Ecology and Management, Vol. 545. ID 121275 (year - 2023) https://doi.org/10.1016/j.foreco.2023.121275

2. Khanina L.G., Bobrovsky M.V., Smirnov V.E. Динамика запасов углерода и азота в валеже и почве после массового ветровала в широколиственном лесу на песках Эволюция биосферы, биогеохимические циклы и биогеохимические технологии: связь фундаментальных и прикладных исследований. Материалы конференции. Пущино: КМК., С. 167–170. (year - 2023)

3. Khanina L.G., Bobrovsky M.V., Smirnov V.E. Динамика запасов биофильных элементов в валеже и почве после массового ветровала в широколиственном лесу на флювиогляциальных песках Вестник Томского государственного университета. Биология., № 62. С. 29–52. (year - 2023) https://doi.org/10.17223/19988591/62/2

4. Khanina L.G., Smirnov V.E., Bobrovskii M.V. The elemental composition of dead wood in different tree species and stages of decay in the broad-leaved forests of the Kaluzhskie Zaseki State Nature Reserve Russian Journal of Ecology, Vol. 54. No. 7. P. 658–670. (year - 2023) https://doi.org/10.1134/S106741362307007X

5. Khanina L.G., Volobuev S.V., Shelyakin P.V., Tutukina M.N., Bessonova T.A., Bobrovsky M.V. Разнообразие ксилобионтных грибов на валеже в смешанном широколиственном лесу зависит от стадии разложения и видовой принадлежности упавших деревьев Математическое моделирование в экологии / Материалы Восьмой Национальной научной конференции с международным участием, 9–11 ноября 2023 г. Пущино, ФИЦ ПНЦБИ РАН, С. 114-116. (year - 2023)

6. Volobuev S., Shelyakin P., Tutkina M., Bessonova T., Bobrovsky M., Khanina L. Метагеномный анализ ксилобионтных грибов и бактерий в валежных стволах лиственных и хвойных деревьев разной степени разложения после массового ветровала в полидоминантном широколиственном лесу Proceedings of 11th Moscow Conference on Computational Molecular Biology MCCMB'23. М.: ИППИ РАН., 158 (year - 2023)

7. Volobuev S.V. Предпочтения и исключения: особенности субстратной приуроченности ксилотрофных афиллофороидных грибов в полидоминантном широколиственном лесу Экология грибов и грибоподобных организмов: факты, гипотезы, тенденции: сборник материалов конференции / Ярославль. Ярославский государственный педагогический университет им. К. Д. Ушинского, С. 15-16 (year - 2023)

8. Khanina L.G., Bobrovsky M.V., Ivanova N.V. Element content, bulk density, and pH in living and dead wood of 8 tree species in a mesic broad-leaved forest in the East European Plain Institute of Mathematical Problems of Biology RAS – the Branch of Keldysh Institute of Applied Mathematics of Russian Academy of Sciences. Sampling event dataset https://doi.org/10.15468/epgqxk accessed via GBIF.org, Sampling event dataset https://doi.org/10.15468/epgqxk accessed via GBIF.org on 2023-11-25 (year - 2023) https://doi.org/10.15468/epgqxk

9. Volobuev S.V. Aphyllophoroid fungi (Basidiomycota) at deadwood of eight tree species in a mesic broad-leaved forest after a catastrophic windthrow in the Kaluzhskie Zaseki Reserve. Institute of Mathematical Problems of Biology RAS – the Branch of Keldysh Institute of Applied Mathematics of Russian Academy of Sciences. Occurrence dataset. https://doi.org/10.15468/bq3m99 accessed via GBIF.org, Occurrence dataset https://doi.org/10.15468/bq3m99 accessed via GBIF.org (year - 2023) https://doi.org/10.15468/bq3m99


Annotation of the results obtained in 2022
The aim of the project is to study the diversity of xylobiont fungi and bacteria in natural conditions of multi-species mesic Quercus broad-leaved forest at the sites of mass windthrow, which occurred in 2006 in the Kaluzhskie Zaseki State Nature Reserve. In 2022 the following work on the project was done. For metagenomic analysis, 150 wood samples of eight tree species at five decomposition stages were collected. Samples were taken from the deadwood (coarse woody debris) of seven angiosperms: maple (Acer platanoides), birch (Betula pendula), common ash (Fraxinus excelsior), European aspen (Populus tremula), English oak (Quercus robur), linden (Tilia cordata) and elm (Ulmus glabra), and one gymnosperm European spruce (Picea abies). A protocol for the isolation of total DNA and amplification of ITS variant sites of fungi was optimized. Total DNA of bacteria and fungi was isolated from the wood samples. More than half of the samples were prepared for subsequent sequencing. An analysis of the nucleotide sequences in the available international database NCBI Nucleotide was performed. It was found that 12 of the 175 species (6.9%) of aphyllophoroid fungi previously identified in the study area lacked any information on the nucleotide sequences of the ITS1-5.8S-ITS2 region of rDNA. Additionally, a list of 11 taxa represented in the NCBI Nucleotide database for one and two target nucleotide sequences was compiled. DNA barcoding of fruit body samples of xylobiont aphyllophoroid fungi collected from deadwood of eight tree species at the sampling sites for metagenomic analysis was performed. Fifty nucleotide sequences of the ITS1-5.8S-ITS2 region of rDNA were obtained. Preservation of isolated genomic DNA samples in freezer conditions was ensured, determining the possibility of amplification of additional genomic regions. There was created a library of 50 reference macromycete DNA sequences prepared for deposit into the international NCBI GenBank database. For each nucleotide sequence, information is given on the DNA extraction method, the primers used, and details about the voucher herbarium sample. Based on the sequencing of pilot samples, we identified bacterial and fungal communities containing representatives with known wood decomposition ability. In the process of wood decomposition, a pronounced increase in the diversity of the bacterial community was observed, with an increase in both the number of bacterial groups observed and the Shannon index (no individual dominant groups emerged). The relative number of representatives of the Acidobacteriaceae family and Bacili class decreased and the number of representatives of the Actinobacteria and Thermoleophilia classes increased. In the fungal community, there was a change in the dominance of representatives of Ascomycota at the first stage of decomposition by representatives of Basidiomycota at the third stage. The fungal community looked less diverse than the bacterial community: during wood decomposition, the diversity of the fungal community significantly decreased, both in terms of number of observed groups and the Shannon index (dominant groups in the community appeared). Physical and chemical properties of deadwood of 8 tree species at 5 stages of decomposition were analyzed. Analysis of variance showed that tree species differed most clearly in the content and mass concentration of Mn, Zn, Mg, Ca and K. The leaders in the content of these elements were: Mn – Acer, Betula, Picea, and Tilia; Zn – Betula and Populus; Mg – Acer and Ulmus; Ca - Ulmus; K – Tilia and Ulmus. Wood decomposition stages proved to be a significant factor of variation for the content of Mn, P, Cu, Zn, and Ca. In the course of deadwood degradation, there was a significant increase in the content of these elements. It is concluded that the maintenance of cycles of biophilic elements is more successfully implemented in the presence of dead wood of different species at different stages of decomposition. The results of the study were reported in oral presentations at 6 scientific conferences. There were four articles in scientific issues. One of the articles was accepted for publication in the journal Lesovedeine indexed by the SCOPUS database (to be published in 2023). Web pages on the project have been developed; they are located at: https://www.impb.ru/index.php?id=div/lce/22-24-01063/index https://www.impb.ru/index.php?id=div/lce/22-24-01063/index&lang=en

 

Publications

1. Khanina L.G., Bobrovsky M.V., Smirnov V.E., Ivashchenko K.V., Zhuravleva A.I., Romanov M.S. МАЛОНАРУШЕННЫЕ ШИРОКОЛИСТВЕННЫЕ ЛЕСА: УГЛЕРОД И АЗОТ В ВАЛЕЖЕ И ПОЧВАХ НА МАССОВЫХ ВЕТРОВАЛАХ Научные основы устойчивого управления лесами. Материалы Всероссийской научной конференции с международным участием, посвященной 30-летию ЦЭПЛ РАН. Москва., С. 130-133 (year - 2022)

2. Khanina, L., Volobuev, S., Smirnov, V., Tutukina, M., Shelyakin, P., Bobrovsky, M. Разнообразие сообществ грибов и бактерий и динамика физико-химических свойств древесины широколиственных видов деревьев и ели европейской в ходе деструкции после массового ветровала Доклады Международной конференции “Математическая биология и биоинформатика”. Под ред. В.Д. Лахно. Том 9. Пущино: ИМПБ РАН, 2022., Том 9. Статья № e47. doi: 10.17537/icmbb22.48 (year - 2022) https://doi.org/10.17537/icmbb22.48

3. Khanina, L.G., Smirnov, V.E., Bobrovsky, M.V. Влияние древесной породы и стадии разложения на элементный состав валежа в широколиственном лесу на Восточно-Европейской равнине. Лесоведение, № 4. С. 353-368. (year - 2023) https://doi.org/10.31857/S0024114823040034

4. Volobuev S.V. POLYPORE FUNGI SPECIES (AGARICOMYCETES, BASIDIOMYCOTA), NEW TO THE KALUGA REGION, RUSSIA Разнообразие растительного мира, № 4 (15). С. 43–49. (year - 2022) https://doi.org/10.22281/2686-9713-2022-4-43-49

5. Volobuev, S.V., Khanina, L.G. РАЗНООБРАЗИЕ КСИЛОТРОФНЫХ АФИЛЛОФОРОИДНЫХ ГРИБОВ НА МАССОВОМ ВЕТРОВАЛЕ В ПОЛИДОМИНАНТНОМ ШИРОКОЛИСТВЕННОМ ЛЕСУ (НА ПРИМЕРЕ ЗАПОВЕДНИКА «КАЛУЖСКИЕ ЗАСЕКИ») МАТЕРИАЛЫ V (ХIII) МЕЖДУНАРОДНОЙ БОТАНИЧЕСКОЙ КОНФЕРЕНЦИИ МОЛОДЫХ УЧЁНЫХ В САНКТ-ПЕТЕРБУРГЕ Сборник тезисов конференции. Санкт-Петербург, 2022 Издательство: Ботанический институт им. В.Л. Комарова РАН (Санкт-Петербург), с. 130 (year - 2022)

6. Bobrovsky M.V., Khanina L.G., Ivanova N.V. Sampling-event data of vascular plant species on microsites caused by catastrophic windthrow in a mesic broad-leaved forest (Kaluzhskie Zaseki State Nature Reserve). Version 1.6. Institute of Mathematical Problems of Biology RAS – the Branch of Keldysh Institute of Applied Mathematics of Russian Academy of Sciences. Sampling event dataset https://doi.org/10.15468/zkbjxv accessed via GBIF.org on 2022-12-07., Sampling event dataset https://doi.org/10.15468/zkbjxv accessed via GBIF.org on 2022-12-07. (year - 2022) https://doi.org/10.15468/zkbjxv