Biogeochemical Stoichiometry Reveals P and N Limitation Across the Post-glacial Landscape of Denali National Park, Alaska

Dual-Domain Primary Succession of Bacteria in Glacier Forefield Streams and Soils of a Maritime and Continental Glacier

Glaciers retreat rapidly and create newly exposed terrestrial and aquatic habitats in glacier forefields, where primary succession proceeds synchronously in glacier forefields. Here, we introduced the "Dual-Domain Primary Succession" concept to examine the parallel yet distinct primary succession processes in soil and stream ecosystems within glacier forefields, by focusing on Hailuogou Glacier and Urumqi Glacier No.1 in China. Findings showed that soil bacterial communities exhibited higher α-diversity with a decreasing pattern in Hailuogou Glacier, in contrast to Urumqi Glacier No.1, which displayed lower and unimodally distributed α-diversity along the glacier forefield chronosequence (GFC). A similar pattern emerged in streams, except for an increasing α-diversity trend in Urumqi Glacier No.1 stream along the GFC. Additionally, α-diversity in streams changed more rapidly than in soils for Hailuogou Glacier, but more slowly for Urumqi Glacier No.1. Along GFC, both soil and stream bacterial communities experienced spatial variations, primarily due to species turnover. The succession of community composition was evident at the OTU level, with each module in the co-occurrence network consisting of OTUs enriched at specific successional stages. A substantial number of OTUs shared between paired soil and stream samples showed a decreasing trend along the GFC, while β-diversity increased. The results suggested that bacterial communities have a similar succession pattern but in different pace between soil and stream while having distinct successional trajectories between the studied glaciers. This study highlighted the "Dual-Domain Primary Succession" in glacier forefields, but further studies with more glaciers are necessary to make broader generalizations.

Enzyme stoichiometry reveals microbial nitrogen limitation in stony soils

Resource limitation for soil microorganisms is the crucial factor in nutrient cycling and vegetation development, which are especially important in arid climate. Given that rock fragments strongly impact hydrologic and geochemical processes in arid areas, we hypothesized that microbial resource (C and N) limitation will increase along the rock fragment content (RFC) gradient. We conducted a field experiment in Minjiang river arid valleys with four RFC content (0 %, 25 %, 50 %, and 75 %, V V-1) and four vegetation types (Artemisia vestita, Bauhinia brachycarpa, Sophora davidii, and the soil without plants). Activities of C (β-1,4-glucosidase, BG), N (β-1,4-N-acetyl-glucosaminidase, NAG; L-leucine aminopeptidase, LAP), and P (acid phosphatase, ACP) acquiring enzymes were investigated to assess the limitations by C, N or P. In unplanted soil, the C acquiring enzyme activity decreased by 43 %, but N acquiring enzyme activity increased by 72 % in 75 % RFC than those in rock-free soils (0 % RFC). Increasing RFC reduced C:N and C:P enzymatic ratios, as well as vector length and vector angle (< 45°). Plants increased the activities of C and N acquiring enzymes in soils, as well as C:P and N:P enzyme activities, as well as vector length (by 5.6 %-25 %), but decreased vector angle (by 13 %-21 %). Enzyme stoichiometry was dependent on biotic and abiotic factors, such as soil water content, soil C:N, and total content of phospholipid fatty acids, reflecting microbial biomass content. Increased RFC shifted enzymatic stoichiometry toward lower C but stronger N limitation for microorganisms. Vegetation increased microbial C and N limitation, and impacted the enzymatic activities and stoichiometry depending on shrub functional groups. Consequently, the direct effects of vegetation, nutrient availability and microbial biomass content, as well as indirect effects of soil properties collectively increased microbial resource limitations along the RFC gradient.

Microbial nutrient limitation and carbon use efficiency changes under different degrees of litter decomposition

Alpine ecosystems are important terrestrial carbon (C) pools, and microbial decomposers play a key role in litter decomposition. Microbial metabolic limitations in these ecosystems, however, remain unclear. The objectives of this study aim to elucidate the characteristics of microbial nutrient limitation and their C use efficiency (CUE), and to evaluate their response to environmental factors. Five ecological indicators were utilized to assess and compare the degree of microbial elemental homeostasis and the nutrient limitations of the microbial communities among varying stages of litter decomposition (L, F, and H horizon) along an altitudinal gradient (2800, 3000, 3250, and 3500 m) under uniform vegetation (Abies fabri) on Gongga Mountain, eastern Tibetan Plateau. In this study, microorganisms in the litter reached a strictly homeostatic of C content exclusively during the middle stage of litter decomposition (F horizon). Based on the stoichiometry of soil enzymes, we observed that microbial N- and P-limitation increased during litter degradation, but that P-limitation was stronger than N-limitation at the late stages of degradation (H horizon). Furthermore, an increase in microbial CUE corresponded with a reduction in microbial C-limitation. Additionally, redundancy analysis (RDA) based on forward selection further showed that microbial biomass C (MBC) is closely associated with the enzyme activities and their ratios, and MBC was also an important factor in characterizing changes in microbial nutrient limitation and CUE. Our findings suggest that variations in MBC, rather than N- and P-related components, predominantly influence microbial metabolic processes during litter decomposition on Gongga Mountain, eastern Tibetan Plateau.

Principal role of fungi in soil carbon stabilization during early pedogenesis in the high Arctic

Climate warming is causing widespread deglaciation and pioneer soil formation over glacial deposits. Melting glaciers expose rocky terrain and glacial till sediment that is relatively low in biomass, oligotrophic, and depleted in nutrients. Following initial colonization by microorganisms, glacial till sediments accumulate organic carbon and nutrients over time. However, the mechanisms driving soil nutrient stabilization during early pedogenesis after glacial retreat remain unclear. Here, we traced amino acid uptake by microorganisms in recently deglaciated high-Arctic soils and show that fungi play a critical role in the initial stabilization of the assimilated carbon. Pioneer basidiomycete yeasts were among the predominant taxa responsible for carbon assimilation, which were associated with overall high amino acid use efficiency and reduced respiration. In intermediate- and late-stage soils, lichenized ascomycete fungi were prevalent, but bacteria increasingly dominated amino acid assimilation, with substantially decreased fungal:bacterial amino acid assimilation ratios and increased respiration. Together, these findings demonstrate that fungi are important drivers of pedogenesis in high-Arctic ecosystems that are currently subject to widespread deglaciation from global warming.

Mineral substrate quality determines the initial soil microbial development in front of the Nordenskiöldbreen, Svalbard

Substrate geochemistry is an important factor influencing early microbial development after glacial retreat on nutrient-poor geological substrates in the High Arctic. It is often difficult to separate substrate influence from climate because study locations are distant. Our study in the retreating Nordenskiöldbreen (Svalbard) is one of the few to investigate biogeochemical and microbial succession in two adjacent forefields, which share the same climatic conditions but differ in their underlying geology. The northern silicate forefield evolved in a classical chronosequence, where most geochemical and microbial parameters increased gradually with time. In contrast, the southern carbonate forefield exhibited high levels of nutrients and microbial biomass at the youngest sites, followed by a significant decline and then a gradual increase, which caused a rearrangement in the species and functional composition of the bacterial and fungal communities. This shuffling in the early stages of succession suggests that high nutrient availability in the bedrock could have accelerated early soil succession after deglaciation and thereby promoted more rapid stabilization of the soil and production of higher quality organic matter. Most chemical parameters and bacterial taxa converged with time, while fungi showed no clear pattern.

Diversity and co-occurrence networks of bacterial and fungal communities on two typical debris-covered glaciers, southeastern Tibetan Plateau

Debris-covered glaciers (DCGs) are globally distributed and thought to contain greater microbial diversity than clean surface continental glaciers, but the ecological characteristics of microbial communities on the surface of DCGs have remained underexplored. Here, we investigated bacterial and fungal diversity and co-occurrence networks on the supraglacial debris habitat of two DCGs (Hailuogou and Dagongba Glaciers) in the southeastern Tibetan Plateau. We found that the supraglacial debris harbored abundant microbes with Proteobacteria occupying more than half (51.5%) of the total bacteria operational taxonomic units. The composition, diversity, and co-occurrence networks of both bacterial and fungal communities in the debris were significantly different between Hailuogou Glacier and Dagongba Glacier even though the glaciers are geographically adjacent within the same mountain range. Bacteria were more diverse in the debris of the Dagongba Glacier, where a lower surface velocity and thicker debris layer allowed the supraglacial debris to continuously weather and accumulate nutrients. Fungi were more diverse in the debris of the Hailuogou Glacier, which experiences a wetter monsoonal climate, is richer in calcium, has greater debris instability, and greater ice velocity than the Dagongba Glacier. These factors may provide ideal conditions for the dispersal and propagation of fungi spores on the Hailuogou Glacier. In addition, we found an obvious gradient of bacterial diversity along the supraglacial debris transect on the Hailuogou Glacier. Bacterial diversity was lower where debris cover was thin and scattered and became more diverse near the glacial terminus in thick, slow-moving debris. No such increasing bacterial pattern was detected on the Dagongba Glacier, which implies a positive relationship of debris age, thickness, and weathering on bacterial diversity. Additionally, a highly connected bacterial co-occurrence network with low modularity was found in the debris of the Hailuogou Glacier. In contrast, debris from the Dagongba Glacier exhibited less connected but more modularized co-occurrence networks of both bacterial and fungal communities. These findings indicate that less disturbed supraglacial debris conditions are crucial for microbes to form stable communities on DCGs.

Karst tiankeng create a unique habitat for the survival of soil microbes: Evidence from ecoenzymatic stoichiometry

Clarifying the soil microbial metabolism and resource limitations could help to understand the functions and processes of aboveground ecosystems, as well as to predict ecosystem stability under global climate change. Karst tiankeng is a kind of large-scale negative surface terrain on the surface which is similar to an oasis in degraded karst landscapes, but their soil microbial resource limitations still unclear. In this study, we evaluated and compared the soil microbial resource limitation in non-degraded tiankeng (NDT), moderately degraded tiankeng (MDT), heavily degraded tiankeng (HDT), and outside tiankeng (OT) by calculating soil ecoenzymatic stoichiometry. Overall, soil microbial communities were more limited by C and P in karst tiankeng ecosystem. The soil microbial C and P limitations significantly differed with the karst tiankeng degradation increased, and the lowest C and P limitations were observed in NDT. The higher microbial C and P limitations were observed in OT. Linear regression and redundancy analysis indicated that soil microbial C and P limitations were significantly influenced by soil nutrients. Karst tiankeng degradation influence the biogeochemical cycle and function of karst tiankeng systems. Our results highlight that karst tiankeng (especially the NDT) can provide a stable habitat for the survival of microorganisms in karst areas. Karst tiankeng is essential for regional ecological restoration and biodiversity conservation.

Abundant and rare soil fungi exhibit distinct succession patterns in the forefield of Dongkemadi glacier on the central Qinghai-Tibet Plateau

Glaciers are retreating rapidly, exposing extensive new soil habitats in glacier forefields and providing unique areas for studying primary succession. However, understanding the variation patterns and assembly mechanisms of abundant and rare fungi subcommunities along the glacier-retreating chronosequence remains a knowledge gap, especially true for the vast Qinghai-Tibet Plateau (QTP). Here, we investigated fungal communities in the glacier forefield in Dongkemadi Glaicer on the central QTP. The results showed that fungal alpha diversity exhibited a clear increasing pattern in response to increasing of distance to glacier. The percentage of abundant OTUs decreased while the percentage of rare OTUs increased, suggesting that soil development is more beneficial to the rare taxa. The distributions of both abundant and rare subcommunities exhibited a clear spatial pattern along the distance to glacier, and might be strongly controlled by multiple environmental variables, including pH, soil moisture, vegetation status, soil organic carbon, total nitrogen, and soluble reactive phosphorus. Abundant and rare fungal subcommunities were structured in different assembly regimes. Dispersal limitation processes were dominant for both abundant and rare subcommunities but with a stronger contribution to abundant subcommunity assembly. Heterogeneous selection processes contributed higher and non-dominant processes contributed lower to abundant subcommunities than to rare subcommunities. The modular structure of the fungal co-occurrence network was highly localized along the soil chronosequence. By revealing distinct diversity patterns and community assembly mechanisms of abundant and rare fungal subcommunities, our study improved our understanding of ecological succession along the glacier-retreating chronosequence.

Soil Microbial Resource Limitations and Community Assembly Along a Camellia oleifera Plantation Chronosequence

Understanding soil microbial element limitation and its relation with the microbial community can help in elucidating the soil fertility status and improving nutrient management of planted forest ecosystems. The stand age of a planted forest determines the aboveground forest biomass and structure and underground microbial function and diversity. In this study, we investigated 30 plantations of Camellia oleifera distributed across the subtropical region of China that we classified into four stand ages (planted <9 years, 9–20 years, 21–60 years, and >60 years age). Enzymatic stoichiometry analysis showed that microbial metabolism in the forests was mainly limited by C and P. P limitation significantly decreased and C limitation slightly increased along the stand age gradient. The alpha diversity of the soil microbiota remained steady along stand age, while microbial communities gradually converged from scattered to clustered, which was accompanied by a decrease in network complexity. The soil bacterial community assembly shifted from stochastic to deterministic processes, which probably contributed to a decrease in soil pH along stand age. Our findings emphasize that the stand age regulated the soil microbial metabolism limitation and community assembly, which provides new insight into the improvement of C and P management in subtropical planted forest.

Incorporating Biogeochemistry into Dryland Restoration

Dryland degradation is a persistent and accelerating global problem. Although the mechanisms initiating and maintaining dryland degradation are largely understood, returning productivity and function through ecological restoration remains difficult. Water limitation commonly drives slow recovery rates within drylands; however, the altered biogeochemical cycles that accompany degradation also play key roles in limiting restoration outcomes. Addressing biogeochemical changes and resource limitations may help improve restoration efforts within this difficult-to-restore biome. In the present article, we present a synthesis of restoration literature that identifies multiple ways biogeochemical understandings might augment dryland restoration outcomes, including timing restoration around resource cycling and uptake, connecting heterogeneous landscapes, manipulating resource pools, and using organismal functional traits to a restoration advantage. We conclude by suggesting ways to incorporate biogeochemistry into existing restoration frameworks and discuss research directions that may help improve restoration outcomes in the world's highly altered dryland landscapes.

Changes in leaf litter decomposition of primary Korean pine forests after degradation succession into secondary broad-leaved forests

Forest degradation succession often leads to changes in forest ecosystem functioning. Exactly how the decomposition of leaf litter is affected in a disturbed forest remains unknown. Therefore, in our study, we selected a primary Korean pine forest (PK) and a secondary broad-leaved forest (SF) affected by clear-cutting degradation, both in Northeast China. The aim was to explore the response to changes in the leaf litter decomposition converting PK to SF. The mixed litters of PK and SF were decomposed in situ (1 year). The proportion of remaining litter mass, main chemistry, and soil biotic and abiotic factors were assessed during decomposition, and then, we made an in-depth analysis of the changes in the leaf litter decomposition. According to our results, leaf litter decomposition rate was significantly higher in the PK than that in the SF. Overall, the remaining percent mass of leaf litter's main chemical quality in SF was higher than in PK, indicating that leaf litter chemical turnover in PK was relatively faster. PK had a significantly higher amount of total phospholipid fatty acids (PLFAs) than SF during decomposition. Based on multivariate regression trees, the forest type influenced the soil habitat factors related to leaf litter decomposition more than decomposition time. Structural equation modeling revealed that litter N was strongly and positively affecting litter decomposition, and the changes in actinomycetes PLFA biomass played a more important role among all the functional groups. Selected soil abiotic factors were indirectly driving litter decomposition through coupling with actinomycetes. This study provides evidence for the complex interactions between leaf litter substrate and soil physical-chemical properties in affecting litter decomposition via soil microorganisms.

Ecological stoichiometric characteristics of Carbon (C), Nitrogen (N) and Phosphorus (P) in leaf, root, stem, and soil in four wetland plants communities in Shengjin Lake, China

Ecological stoichiometric should be incorporated into management and nutrient impacted ecosystems dynamic to understand the status of ecosystems and ecological interaction. The present study focused on ecological stoichiometric characteristics of soil, and leaves, stems, and roots of different macrophytes after the banning of seine fishing in Shengjin Lake. For C, N, and P analysis from leaves, stems, roots, and soil to explore their stoichiometric ratio and deriving environmental forces, four dominant plant communities (Vallisneria natans, Zizania latifolia, Trapa natans and Carex schmidtii) were collected. The concentration of C, N, P and C: N: P ratio in leaves, stems, roots, and soil among the plant communities varied significantly. Along the depth gradient high C: N was measured in C.schmidtii soil (7.08±1.504) but not vary significantly (P >0.05). High C: P result was found in T.natans (81.14±43.88) and in V.natans soil (81.40±42.57) respectively with no significant difference (p>0.05). Besides, N: P ratio measured high in V. natans (13.7±4.05) and showed significant variation (P<0.05). High leaf C: N and N: P ratio was measured in C. schmidtii and V. natans respectively. Nevertheless, high leaf C: P ratio was measured in Z. latifolia. From the three studied organs, leaf C: N and N: P ratio showed high values compared to root and stems. The correlation analysis result showed that at 0-10cm depth soil organic carbon (SOC) correlated negatively with stem total phosphorus (STP), and root total nitrogen (RTN) (P<0.05) but positively strongly with leaf total phosphorus (LTP) and leaf total nitrogen (LTN) (P<0.01) respectively. Soil total nitrogen (STN) at 0-10cm strongly positively correlated with leaf total phosphorus (LTP) (P<0.01) and positively with RN: P and leaf total carbon (LTC) (P<0.05). Soil basic properties such as soil moisture content (SMC), bulky density (BD) and pH positively correlated with soil ecological stoichiometric characteristics. Redundancy analysis (RDA) result showed available nitrogen (AN), soil total nitrogen (STN), and available phosphorus (AP) were the potential determinants variables on plants stoichiometric characteristics.

Variation in denitrifying bacterial communities along a primary succession in the Hailuogou Glacier retreat area, China

Background

The Hailuogou Glacier is located at the Gongga Mountain on the southeastern edge of the Tibetan Plateau, and has retreated continuously as a result of global warming. The retreat of the Hailuogou Glacier has left behind a primary succession along soil chronosequences. Hailuogou Glacier’s retreated area provides an excellent living environment for the colonization of microbes and plants, making it an ideal model to explore plant successions, microbial communities, and the interaction of plants and microbes during the colonization process. However, to date, the density of the nitrogen cycling microbial communities remain unknown, especially for denitrifiers in the primary succession of the Hailuogou Glacier. Therefore, we investigated the structural succession and its driving factors for denitrifying bacterial communities during the four successional stages (0, 20, 40, and 60 years).

Methods

The diversity, community composition, and abundance of nosZ-denitrifiers were determined using molecular tools, including terminal restriction fragment length polymorphism and quantitative polymerase chain reactions (qPCR).

Results

nosZ-denitrifiers were more abundant and diverse in soils from successional years 20–60 compared to 0–5 years, and was highest in Site3 (40 years). The denitrifying bacterial community composition was more complex in older soils (40–60 years) than in younger soils (≤20 years). The terminal restriction fragments (T-RFs) of Azospirillum (90 bp) and Rubrivivax (95 bp) were dominant in soisl during early successional stages (0–20 years) and in the mature phase (40–60 years), respectively. Specific T-RFs of Bradyrhizobium (100 bp) and Pseudomonas (275 bp) were detected only in Site3 and Site4, respectively. Moreover, the unidentified 175 bp T-RFs was detected only in Site3. Of the abiotic factors that were measured in this study, soil available phosphorus, available potassium and denitrifying enzyme activity (DEA) correlated significantly with the community composition of nosZ-denitrifiers (P < 0.05 by Monte Carlo permutation test within RDA analysis).

High nitrogen contribution by Gunnera magellanica and nitrogen transfer by mycorrhizas drive an extraordinarily fast primary succession in sub-Antarctic Chile

Chronosequences at the forefront of retreating glaciers provide information about colonization rates of bare surfaces. In the northern hemisphere, forest development can take centuries, with rates often limited by low nutrient availability. By contrast, in front of the retreating Pia Glacier (Tierra del Fuego, Chile), a Nothofagus forest is in place after only 34 yr of development, while total soil nitrogen (N) increased from near zero to 1.5%, suggesting a strong input of this nutrient. We measured N-fixation rates, carbon fluxes, leaf N and phosphorus contents and leaf δ¹⁵ N in the dominant plants, including the herb Gunnera magellanica, which is endosymbiotically associated with a cyanobacterium, in order to investigate the role of N-fixing and mycorrhizal symbionts in N-budgets during successional transition. G. magellanica presented some of the highest nitrogenase activities yet reported (potential maximal contribution of 300 kg N ha^-1  yr^-1 ). Foliar δ¹⁵ N results support the framework of a highly efficient N-uptake and transfer system based on mycorrhizas, with c. 80% of N taken up by the mycorrhizas potentially transferred to the host plant. Our results suggest the symbiosis of G. magellanica with cyanobacteria, and trees and shrubs with mycorrhizas, to be the key processes driving this rapid succession.

Resource limitation of soil microbes in karst ecosystems

Knowledge about resource limitation to soil microbes is crucial for understanding ecosystem functions and processes, and for predicting ecosystem responses to global changes as well. Karst ecosystems are widespread in the world, and play a key role in regulating the global climate, however, the patterns of and mechanisms underlying microbial resource limitation in karst ecosystems remain poorly known. Here we investigated the microbial resource limitation in a karst region, by selecting four main land-use types, i.e. cropland, grassland, shrubland and secondary forest, in areas underlain by two lithology types, i.e. dolomite and limestone, in southwest China. Ecoenzymatic stoichiometry was used as an indicator of microbial resource limitation. Overall, soil microbes in karst ecosystems were more limited by carbon and phosphorus, rather than by nitrogen. Further analyses revealed that the patterns of carbon and phosphorus limitation were different among land-use or lithology types. Microbial carbon limitation was greatest in cropland and forest but lowest in grassland, and was greater under dolomite than under limestone. Microbial phosphorus limitation decreased from secondary forest to cropland under dolomite areas, but showed no difference among ecosystem types under limestone areas, indicating that lithology controls the pattern of microbial phosphorus limitation along the post-agriculture succession. Our study describes a general pattern of microbial resource limitation in karst ecosystems, and we suggest that lithology may provide a new mechanism for explaining the variations of microbial resource limitation along the post-agriculture succession in different regions.

Climatic and topographic variables control soil nitrogen, phosphorus, and nitrogen: Phosphorus ratios in a Picea schrenkiana forest of the Tianshan Mountains

Knowledge about soil nitrogen (N) and phosphorus (P) concentrations, stocks, and stoichiometric ratios is crucial for understanding the biogeochemical cycles and ecosystem function in arid mountainous forests. However, the corresponding information is scarce, particularly in arid mountainous forests. To fill this gap, we investigated the depth and elevational patterns of the soil N and P concentrations and the N: P ratios in a Picea schrenkiana forest using data from soil profiles collected during 2012-2017. Our results showed that the soil N and P concentrations and the N: P ratios varied from 0.15 g kg-1 to 0.56 g kg-1 (average of 0.31 g kg-1), from 0.09 g kg-1 to 0.16 g kg-1 (average of 0.12 g kg-1), and from 2.42 g kg-1 to 4.36 g kg-1 (average of 3.42 g kg-1), respectively; additionally, values significantly and linearly decreased with soil depth. We did not observe a significant variation in the soil N and P concentrations and the N: P ratios with the elevational gradient. In contrast, our results revealed that the mean annual temperature and mean annual precipitation exhibited a more significant influence on the soil N and P concentrations and the N: P ratios than did elevation. This finding indicated that climatic variables might have a more direct impact on soil nutrient status than elevation. The observed relationship among the soil N and P concentrations and the N: P ratios demonstrated that the soil N was closely coupled with the soil P in the P. schrenkiana forest.

Phosphorus, not nitrogen, limits plants and microbial primary producers following glacial retreat

Current models of ecosystem development hold that low nitrogen availability limits the earliest stages of primary succession, but these models were developed from studies conducted in areas with temperate or wet climates. Global warming is now causing rapid glacial retreat even in inland areas with cold, dry climates, areas where ecological succession has not been adequately studied. We combine field and microcosm studies of both plant and microbial primary producers and show that phosphorus, not nitrogen, is the nutrient most limiting to the earliest stages of primary succession along glacial chronosequences in the Central Andes and central Alaska. We also show that phosphorus addition greatly accelerates the rate of succession for plants and for microbial phototrophs, even at the most extreme deglaciating site at over 5000 meters above sea level in the Andes of arid southern Peru. These results challenge the idea that nitrogen availability and a severe climate limit the rate of plant and microbial succession in cold-arid regions and will inform conservation efforts to mitigate the effects of global change on these fragile and threatened ecosystems.

Spatial autocorrelation of microbial communities atop a debris-covered glacier is evidence of a supraglacial chronosequence

Although microbial communities from many glacial environments have been analyzed, microbes living in the debris atop debris-covered glaciers represent an understudied frontier in the cryosphere. The few previous molecular studies of microbes in supraglacial debris have either had limited phylogenetic resolution, limited spatial resolution (e.g. only one sample site on the glacier) or both. Here, we present the microbiome of a debris-covered glacier across all three domains of life, using a spatially-explicit sampling scheme to characterize the Middle Fork Toklat Glacier's microbiome from its terminus to sites high on the glacier. Our results show that microbial communities differ across the supraglacial transect, but surprisingly these communities are strongly spatially autocorrelated, suggesting the presence of a supraglacial chronosequence. This pattern is dominated by phototrophic microbes (both bacteria and eukaryotes) which are less abundant near the terminus and more abundant higher on the glacier. We use these data to refute the hypothesis that the inhabitants of the glacier are randomly deposited atmospheric microbes, and to provide evidence that succession from a predominantly photosynthetic to a more heterotrophic community is occurring on the glacier.