Next Article in Journal
Chemical Composition and Nutritional Value of Flowers and Fruits of Cytisus striatus (Hill) Rothm
Previous Article in Journal
Identification of DREB Family Genes in Banana and Their Function under Drought and Cold Stress
Previous Article in Special Issue
Relationships within Bolbitis sinensis Species Complex Using RAD Sequencing
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Origin and Persistence of Lycopodium clavatum and Lycopodium annotinum (Lycopodiaceae) in Scots Pine Forests

by
Radvilė Rimgailė-Voicik
1,*,
Aleksandras Voicikas
1,
Julija Fediajevaitė
1,2,
Sigitas Juzėnas
1 and
Jolanta Patamsytė
1
1
Life Sciences Center, Institute of Biosciences, Vilnius University, Saulėtekio Ave. 7, 10257 Vilnius, Lithuania
2
School of Biological Sciences, University of Reading, Whiteknights, Reading RG6 6UR, UK
*
Author to whom correspondence should be addressed.
Plants 2024, 13(15), 2120; https://doi.org/10.3390/plants13152120
Submission received: 29 May 2024 / Revised: 5 July 2024 / Accepted: 23 July 2024 / Published: 31 July 2024
(This article belongs to the Special Issue Diversity and Evolution in Lycophytes and Ferns)

Abstract

:
Understanding the growth dynamics of spore-bearing clonal plant sporophytes and the influence of abiotic and biotic factors is crucial for predicting the persistence of club moss populations and implementing effective habitat management techniques. Despite this, the longevity and development of club-moss populations are rarely studied. This study adopted an integrated approach to assess the probability of repetitive young sporophyte recruitment via sexual propagation in Lycopodium annotinum L. and Lycopodium clavatum L. The size–age problem of clonal spore-bearing forest plants and their niche segregation were addressed. The canopy characteristics, insolation, small-scale disturbance, and genetic polymorphism were studied in temperate semi-natural Scots pine forests in Lithuania. Based on the size of the clones discovered, we hypothesize that initial sporophyte emergence occurred in 20-year-old pine stands, with subsequent sporophyte emergence continuing over time. The emergence was related to small-scale disturbances. High genetic polymorphism indicates that all sporophyte stands studied likely emerged via sexual reproduction. According to Ellenberg values, L. annotinum is related to shady habitats, but our findings show both species coexisting abundantly in the more open habitat, supposedly more suitable for L. clavatum.No significant differences in vegetation relevés and light availability was detected using hemispheric images.

1. Introduction

The balance between sexual recruitment and clonal offspring production is not well understood in land plants, and general explanations for low rates of seedling recruitment and clonality often point to malfunctions in sexual reproduction, such as seed or partner deficiency, trade-offs between clonal and sexual reproductive allocations [1], or environmental restrictions [2]. In woodlands, plants exhibit escape mechanisms to exploit spatially unpredictable patches, making them less dependent on seasonal regeneration [3]. Seedling establishment is less common compared to grasslands, and forest-floor species are usually adapted to long-distance dispersal [4]. Club mosses, in particular, exhibit low competitiveness [5,6] and tend to respond to unfavorable disturbances by growing away from such areas [7].
The life cycle of Lycopodium consists of two separate generations: gametophytes (n) develop in the soil, while sporophytes (2n) are present on the forest floor. Sexually, club-moss sporophytes originate from subterranean, achlorophyllous gametophytes which require a specific group of endophytic fungi to form [8,9]. Recent studies have demonstrated that successfully developing gametophytes are colonized by Densosporaceae (Endogonales and Mucoromycotina) fungi [10]. Continuous vegetative propagation likely enables club-moss sporophytes’ persistence despite low or no sexual reproduction, as the genet turnover rate is very slow. Clonal areas typically form from clumped and intermingled groups of genets and ramets. The club moss’s outgrowth formation and percentage of genets and ramets in a single outgrowth, and their spatial distribution patterns remain unknown. Physiological ramet integration allows for acropetal and basipetal resource translocations, as demonstrated in arctic [11] and temperate deciduous forest [12,13] clonal plants.
The genome size evolution and the degree of plasticity in long-lived clonal archaic spore-bearing vascular plants are of considerable biological interest, but they remain poorly understood [14,15]. Based on the frequency of sexual recruitment, all clonal plants can be separated into two major groups: (a) populations with repeated regular seedling recruitment and (b) populations with initial seedling recruitment and rare seedling addition [16]. It is unclear to which one lycophytes should be assigned due to the wide life-cycle timeframe and the rarity of gametophyte discovery [17]. Juvenile club-moss sporophytes are spatially isolated from mature Lycopodium annotinum L. and Lycopodium clavatum L. plants [8,18], which tend to form large outgrowths [5,7,19].
In clonal plants, mother plants produce ramets only in neighboring patches, and their performance is often regulated by local density rather than the whole population density. The persistence of genets in natural clonal populations largely depends on meristem demography [20]. It is believed that plant roots originated independently in lycophytes and euphyllophytes, but the processes of root evolution remain contentious, and little is known about the parallel evolution and molecular mechanisms governing lycophyte root functions [21]. It is generally accepted that Lycopod roots evolved around 400 million years ago and lack both the vascular cambium and organogenic properties outside the apical meristem [22]. Mature club moss sporophytes are opportunistic foragers that are well adapted to patchy nutrient resources [5,23]. The frequency of club-moss sexual reproduction is debatable and may depend on the level of habitat disturbance. Previously, club moss subterranean gametophytes and juvenile sporophytes were found on and near forest roads and tracks, near lines separating forest blocks and skiing tracks [8,17,24]. The presence of fire adaptations is also expected [19,25].
Lycopodium annotinum and L. clavatum (Figure 1) are crucial forest-floor species that are characteristic of the protected Boreal Europe Forest type Western Taiga (code 9010) [26]. Members of the Lycopodium genus are protected and considered to be of community interest. Their collection from the wild and exploitation may be subject to specific management measures at the national level within European Union member states [27]. Additionally, L. annotinum is listed on the International Union for Conservation of Nature’s (IUCN) Red List of Threatened Species, classified as being of ‘least concern’ in Europe [28]. Live specimens of L. clavatum, as well as dried and fresh plants, including leaves, roots/rootstocks, stems, seeds/spores, bark, and fruits, are subject to Annex D of the EU Wildlife Trade Regulation 318-2008 [29].
Stable understory vegetation species groups have been widely used in the major forest-habitat classifications [26,30] and the Ellenberg indicator value system (EIVs) [31,32]. Understory vegetation significantly influences the composition of future forests in response to disturbance [33] and plays a crucial role in nutrient cycling within forests [34]. Yet, there is still no universal agreement upon the concept of forest health with the respect to the whole forest ecosystem [35].
Club mosses are an important forest-floor component with significant indicator value for forest health assessment. It is unclear whether the origination of new club moss sporophytes is related to higher levels of forest floor disturbance during logging or forest clearing, or if it is linked to repetitive moderate to low disturbances. The impact of disturbance on successful club moss sporophyte emergence through sexual reproduction and maturation was never tested. In general, the life history and population dynamics of forest species are strongly shaped by timber extraction activities, which can disrupt the long-term resilience of forests and reduce the delivery of ecosystem services [36]. We hypothesize that the origination of club moss sporophytes from gametophytes is more closely associated with small-scale disturbances than with one-time catastrophic events such as forest fires or clear-cuts.
In this article, we compare L. annotinum and L. clavatum in Scots pine stands based on sporophyte outgrowth size, surrounding vegetation, light availability, and genetic polymorphism. The selected methods were not invasive and ensured low-to-no damage to club moss stands. This research highlights the importance of understanding the adaptive strategies of club moss species, which can inform forest management practices and conservation efforts. Further studies are needed to explore the long-term impacts of environmental disturbances on club moss populations and their role in forest ecosystem health.

2. Results

Lithuanian forest cadastral data showed that trees in researched sites in Vilnius County were of similar age, with an average age of 91 ± 18 years, ranging from 52 to 127 years. The tallest trees were 24–28 m in height. The tree stands in all sites consisted of Pinus sylvestris L., with no subdominant or shrub layer being present. The most common understory species were Vaccinium vitis-idaea L., V. myrtillus L., and Melampyrum pratense L. The moss cover was almost 100%, consisting mostly of Pleurozium schreberi (Brid.) Mitt., Hylocomium splendens (Hedw.) Schimp., and Dicranum Hedw. sp. We compared vegetation relevés of mature L. annotinum and L. clavatum sites (Veržuva, Nemenčinė and Sakiškės) with juvenile Lycopodium growth sites (Varėna) relevés. Reyni diversity profiles were generated for forest sites, and ANOVA revealed significant effects of log species’ richness (α = 0), F (90) = 3.16, p < 0.03. A post hoc comparisons showed that the Sakiškės site differed from the Varėna site (Tukey HSD, p < 0.01). No significant differences among forest sites were found in Shannon (α = 1) and Simpson (α = 2). A t-test also showed that mature forests differ from juvenile forests in Varėna significantly (t (92) = 2.7, p < 0.01, d = 0.7). ANOVA revealed significant differences among forest sites in the proportion of the most abundant species ((α = Inf), F (3, 90) = 4.8, p < 0.03). The post hoc analysis showed that Varėna differed significantly both from Nemenčinė and Sakiškės (p < 0.01). Additionally, sites with mature sporophytes differed significantly from sites with juvenile sporophytes: t (92) = 3, p < 0.003, d = 0.8.
No significant correlation between forest stand age and clone size was found (Figure 2). According to the outgrowth size, the latest emergence of L. annotinum happened when the forest stand was 110 years of age, and for L. clavatum, when the stand was 100 years old. The age of sporophytes varied greatly, with a plausible 20-year gap between youngest and oldest sporophyte. Based on yearly growth rates, the age dispersion of L. annotinum sporophytes was 45 years, and for L. clavatum, it was 40 years. Lycopodium annotinum and L. clavatum populations emerged when trees were around 20 years of age, and population enrichment with new juvenile sporophytes was repetitive.
Weighted Ellenberg indicator values calculated for the grass–subshrub layer were used to test for differences among forest sites. ANOVA revealed a significant effect among forest sites based on species preference in temperature value (F (3, 90) = 3.14, p < 0.03). The post hoc analysis indicated a significant difference between the Nemenčinė and Varėna sites with juvenile sporophytes (Tukey HSD, p < 0.02). Additionally, the t-tests showed significant differences in temperature value (t (92) = 2.6, p < 0.01, d = 0.7) and in nutrients value (t (92) = 2.22, p < 0.03, d = 0.6) when comparing forests with mature sporophytes to those with juvenile sporophytes in Varėna.
In the Nemenčinė and Veržuva sites with L. annotinum stands, overall tree canopy shading ranged from 48 to 74%, and in sites with L. clavatum, shading ranged from 49 to 70%. The variation coefficient for L. annotinum sites was 12.5%, and for L. clavatum sites, it was 12.6%; thus, no statistically significant tree canopy shading differences were determined (Figure 3).
Also, within the outgrowths in the Nemenčinė and Veržuva sites, the determined vegetative shoot cover in 1 m2 was similar and varied from 5 to 50% both in L. annotinum and L. clavatum (Table 1). The average number of strobili in L. annotinum was 10.9 ± 10.3, and in L. clavatum, it was 7 ± 7.1. The maximum number of strobili in L. clavatum was 1077 and in L. annotinum was 2218.
The size of club moss outgrowths varied greatly, but no significant differences were determined between the species. In Nemenčinė and Veržuva sites, the smallest outgrowth of L. annotinum was 0.02 m2, and the largest was 300 m2, with an average of 35 m2 ± 50.8 (n = 45). For L. clavatum, the smallest outgrowth was 0.6 m2, and the largest was 400 m2, with an average of 38.9 m2 ± 73 (n = 31).
Within a ten-meter radius around the L. clavatum and L. annotinum outgrowths, no excessive forest damage was observed; specifically, no stumps, tree cutting, windthrow, or fire marks were found. Only forest paths with moss or grass cover were present near or went through the club-moss stands. The impact of forest paths on the horizontal L. annotinum and L. clavatum clone structure was noticed. The width of the path positively affected the development of young shoots and mature and dried-out shoots and negatively correlated with fertile shoots in L. annotinum. In L. clavatum, path width also positively correlated with young shoots and negatively with mature and dried-out shoots. Both club-moss species reacted to forest road-usage intensity by reducing their sporulation intensity.
The level of polymorphism differed between L. annotinum and L. clavatum species (Table A1). The polymorphism in L. clavatum was higher (74.21%) than in L. annotinum (68.02%). Nei genetic distances [37] differed significantly (PERMANOVA, F = 1.16, p < 0.05) among L. annotinum subpopulations in the Veržuva and Nemenčinė sites, while L. clavatum overlapped (Figure 4). The Mantel test showed a positive significant correlation (r = 0.582, p < 0.05) among Nei distances and geographical coordinates of L. annotinum clones, while the L. clavatum results were insignificant.
We conclude that, in the researched sites, mature L. annotinum and L. clavatum sporophytes occupy complementary niches, thriving in habitats characterized by similar species diversity, recurrent minor forest floor disturbances, and nearly identical light conditions. The ISSR analysis revealed that all outgrowths researched were not identical and, therefore, are a result of club moss sexual propagation. The results presented support the hypothesis that the origination of club moss sporophytes from gametophytes is more closely associated with understory vegetation diversity and small-scale disturbances.

3. Discussion

The balance between sexual and vegetative propagation is a key evolutionary feature of the life history strategies among clonal plants. Clonal spore-bearing plants are vital components of temperate forests, and yet few studies have addressed club moss sporophyte stands, their propagation strategy, and the abiotic and biotic factors influencing their successful sexual propagation. Bruchmann [8] was able to locate club moss gametophytes and juvenile sporophytes only in spruce plantations of 10–20 years of age and argued that only one age class, related to the time of forest planting, can be determined. Our findings align with those of other researchers [19,38], indicating that the emergence of new club moss sporophytes on the forest floor is linked to small-scale disturbances. The absence of these small-scale disturbances, along with interspecific competition, is the primary factor limiting the successful sexual propagation in club mosses.
Only a few studies have addressed the age–size problem in Lycopodium species (Table A2). Wittig [6] suggested that if growth rates during the past hundred years have remained unchanged, the average age of the population can be predicted quite accurately. Oinonen [39,40] constructed a timetable of vegetative spreading for woodland clonal species from parallel measurements and used it to determine plant emergence dates and their relationship with forest floor fires, other clonal growth species, and the age of the forest stand. However, as Oinonen emphasized, it has been difficult to determine the origin of the club-moss outgrowth confidently. After a fire, sporophytes can emerge from gametophytes or develop through vegetative spreading from relict plants that survived. Following partial destruction, new juvenile sporophytes can develop and form intermixed outgrowths [39]. Additionally, L. clavatum is known to form fairy rings, which are associated with nutrient withdrawal from the soil [41].
According to Oinonen [19], vegetative spreading of L. clavatum at normal speed begins after 21–22 years, while for L. annotinum, it starts after 17–18 years. Plagiotropic shoots survive longer than orthotropic shoots and begin to sporulate after three to seven years [5]. It has been proposed [6] that club mosses reach their maximum growth rates 18 years after the emergence of sporophyte. Considering that approximately eight years pass from spore germination, gametophyte maturation, and fertilization, it may take up to 25 years for L. clavatum and L. annotinum plants to begin spreading at their maximum potential. Wittig [6] suggested that, in Europe, L. annotinum can grow up to 20 cm/year. Yet, the growth rates determined vary (Appendix A Table A2). In Lithuania, the annual growth of young L. annotinum sporophytes and growth dynamics in Vaccinium myrtillosum Scots pine stands [38] showed that the populations were quite abundant: more than twenty young sporophytes were present and survived during a four-year study, and the average horizontal growth of young sporophytes from 1987 to 1991 was 10.6 cm (Appendix A Table A3).
The relationship among forest stand age, time since fire, and size of L. annotinum and L. clavatum outgrowths was investigated in Finland [19,39], showing that these factors are related. However, the outgrowth size was more closely related to the time since the last fire than to the age of the forest stand. We were unable to trace the history of forest fires in the study sites, and we found no significant correlation between forest age and sporophyte outgrowth diameter. According to the outgrowth size variation, we conclude that juvenile sporophytes originated from gametophytes repetitively, through different stages of forest development, starting with Scots pine stand of 20 years of age.
According to the Ellenberg values [31,32], L. annotinum has a light value of three and grows in sites where insolation is less than five % and vegetation cover exceeds 95%. Lycopodium clavatum has a light value of eight and grows in areas where insolation is at least 40% and vegetation cover is not more than 60%. In the dry pine forests of Eastern Lithuania, L. annotinum is often found in close proximity to L. clavatum, in habitats more suitable for L. clavatum. It was demonstrated that L. annotinum’s horizontal segment length and number of vertical apices significantly increased under the canopy of Vaccinium myrtillus compared with Vaccinium vitis-idaea; the canopies differed in regard to their red/far-red ratio, with V. myrtillus having the lowest ratio [42], but no similar research is present for L. clavatum. During research conducted in Canada, it was observed that three lycophyte species (L. annotinum, L. clavatum, and Lycopodium dendroideum Michx.) exhibited a notable affinity for old clear-cuts (aged 23–54 years) in mixed Acadian forests [43]. These species were evaluated as mid-seral invaders in the region, persisting longer in the community only if light and other understory conditions remained stable.
More than 95% of young sporophytes found in Scots pine forests in the Varėna District were L. annotinum [18]. It is plausible that L. annotinum exhibit sexual reproduction and, through the gametophytic life-cycle stage, expands the niche, becoming better adapted to environmental pressures related to global warming and anthropogenic contamination than L. clavatum. However, the determined polymorphism in L. clavatum was higher (74.21%) than in L. annotinum (68,.02%). Intensive vegetative propagation is usually associated with genetic monomorphism [44]. This has been shown to be true for L. annotinum stand [6] and population of Lycopodium lucidulum Michx. [45]. However, other genetic studies of clonal plants [46,47] have opposed this view, showing that populations rarely consist of one or few genets. If intragametophytic selfings were dominant, club moss populations would be expected to have low genetic diversity [48], but this was shown not to be the case [49]. The level of polymorphism discovered suggests that all club-moss sporophytes in the researched territory emerged from gametophytes rather than via vegetative propagation.
Scots pine exhibits high ecological plasticity but typically dominates nutrient-poor, dry sites, forming uneven-aged or cohort-structured stands driven by recurrent surface fires and/or gap-and-patch dynamics. In Europe, these stands are primarily shaped by site-specific disturbance regimes [50]. In Lithuania, clear-cutting accounted for 64% of the total area of felling in the state and private sector in 2021 [51]. The typical management cycle of Lithuanian boreal pine forests involves planting or sowing with soil preparation, two or three thinnings at 20-year intervals, and main (clear or selective) cutting at about 100–170 years of age, depending on the forest group [37]. Recently, management techniques oriented closer to natural forest processes, such as shelterwood cuttings, have been implemented [52]. Lycopodium clavatum was the most common species found in all forest types surveyed: clear-cuts, young stands, and mature pine forests; L. annotinum were only found in mature pine forests (unpublished research). Lycopodium clavatum coverage increases to more than 50% in Scots pine stands of 86–100 years of age (6 plots checked), and the first emergence was reported in the age group of 26–35 years (13 plots checked), with up to 20% coverage (Appendix A Table A4). The whole community-disturbance frequency indicator value for L. annotinum is −1.94, suggesting that the return time is nearly a century, while the herb layer disturbance-frequency indicator value is −0.92, indicating adaptations to a small-scale disturbance occurring less than ten years apart [53]. Because of the current forest management policy, forests older than 100 years in Lithuania are rare, and clear-cuts pose a threat to club-moss populations. The main recommendation proposed to keep club-moss sporophytes viable during clear-cut logging is to leave a fragment of tree stand understory [54,55]. Further research is needed to better understand how club mosses are adapting in changing environments and how the origin of club-moss sporophytes influences population vitality and longevity.

4. Materials and Methods

4.1. Study Area

We conducted assessments of dry pine forests in Nemenčinė and Sakiškės, Vilnius District; and Veržuva, Vilnius City. The overall distance between Nemenčinė and Sakiškės was five kilometers, and between Sakiškės and Veržuva, around eight kilometers. Additionally, sites with juvenile club mosses assessed in Varėna District were incorporated into the analysis.
The Scots pine (Pinus sylvestris L.) stands, predominantly of the Vacciniosa and Vaccinio myrtillosa forest types, grow on sandy, nutrient-poor dry podzols and are the most common in Lithuania, covering up to 37% of the total forested area [56]. These oligotrophic pine forests belong to the alliance Dicrano-Pinion sylvestris (Libbert 1933) W. Matuszkiewicz 1962 of the boreal coniferous forests (class Vaccinio-Piceetea).
Overall, the territory can be classified as recent forest (RF), with the age up to two hundred years [57], and has records of agricultural activity, including traditional even-aged (clear-cut) silviculture and salvage logging after natural disturbances. Most pine forests in Southeast Lithuania have a long history of moderate human disturbances and are often visited by locals and holidaymakers, so they are crisscrossed by numerous paths.
The annual gross solar irradiance was similar in all sites, approximately 3300–3400 MJ/m2, while the soil surface albedo was about 20% in Vilnius County sites, and in Varėna District it was lower—less than 17.5%. The mean annual temperature in Vilnius is +6–6.5 °C, with an absolute maximum of +34.9 °C and an absolute minimum of −32.2 °C, and the temperature remains above 0 °C for 248 days a year. The mean annual temperature in Varėna District is +6–7 °C, with an absolute maximum of +35.6 °C and an absolute minimum of −35.9 °C, and the temperature remains above 0 °C for 252 days a year. Snow cover is about 20–25 mm and lasts up to 95–100 days per year in Vilnius County and up to 90 days per year in Varėna District. Mean annual precipitation in all sites was 650 mm [58].

4.2. Sampling Design

From July to September in 2016 and 2017, a total of 76 adult club-moss sporophytes were analyzed by RRV and AV in Vilnius District (Supplementary Table S1): 25 in Nemenčinė site (15 L. annotinum and 10 L. clavatum), 28 in Veržuva site (15 L. annotinum and 13 L. clavatum), and 23 in Sakiškės site (15 L. annotinum and 8 L. clavatum). In 2016, sites with club moss outgrowths were identified. At each site, 10 × 10 m vegetation relevés were performed, the approximate size of each clone was evaluated, and samples for genetic analysis were collected. Only detached stands with a distance of 10 m or more between them were selected. Additionally, 18 vegetation relevés with juvenile club-moss sporophytes assessed in Varėna District in 2013–2015 by RRV were incorporated into the analysis (Figure 5).
We utilized the European Environment Agency (EEA) reference grid, with a resolution of 10 km × 10 km, to illustrate the distribution of club moss sporophyte-sampling sites in Lithuania (Figure 5). Additionally, Figure 5 incorporates detailed data on Scots pine-dominated forest cover in Lithuania, obtained from the Lithuanian Forest Cadastre [59].
In 2017, JF and SJ revisited the outgrowths in the Nemenčinė and Veržuva sites. They assessed the cover of healthy and dried vegetative shoots within the outgrowths, counted the number of strobili, and evaluated the effect of disturbances. For light availability analysis, images were taken at a height of 40 cm from the ground, using a Canon EOS 50D camera with a SIGMA 4.5 mm Circular FISHEYE lens. Hemispheric photos were processed with SideLook 1.1.01 [60] and analyzed with Gap Light Analyzer (GLA) 2.0 [61]. Observations of nearby disturbances, such as paths, pits, and dead trees, were also recorded.

4.3. Statistical Analysis

To evaluate and compare diversity, Renyi diversity profiles among the sites were generated in R using picante package [62]. Renyi diversity profiles were generated for the vegetation relevés of forest sites, and ANOVA was employed to test for significant effects, followed by Tukey HSD post hoc comparisons.
Sporophyte outgrowth size and relative age relationship with forest stand age were addressed using regression equations [19,39] and forest-age data. Forest stand-age data were obtained from the Lithuanian Forest Cadastre [59].
Regressive equations were formulated using the South Finland forest research data, [19], including a correction of eleven meters for secondary L. clavatum stands and eight meters for secondary L. annotinum stands:
dLC = −10.63 + 0.5a
dLA = −7.5 + 0.43a,
Later [39], the equations were adapted as follows:
dLC = −8.175 + 0.455a
dLA = −7.339 + 0.424a,
where d is the clone diameter, a is the age of the sporophyte, LA is L. annotinum, and LC is L. clavatum.
Circular or partially circular stands were included in the analysis as primary clones, while differently shaped outgrowths were attributed to the secondary clones, and their length was used as a diameter.
Only well-defined and reproducible DNA bands were included in the binary data matrix. The percentage of polymorphic loci, population genetic differentiation coefficient (GST), and expected heterozygosity were calculated using program POPGENE, version v. 1.31. Principal Coordinates Analysis (PCoA), Permutation AMOVA, and Mantel test based on Nei’s [63] genetic distances were conducted using GenAlEx v. 6.5 [64].
Data processing, statistical analysis, and illustrations were performed in Python, using pandas [65] and pingouin [66] libraries.

4.4. DNA Analysis

Outgrowths from Nemenčinė and Veržuva were analyzed. DNA was purified from fresh apical branches of L. annotinum and L. clavatum, using the modified CTAB method [67]. DNA concentrations were measured using a BioPhotometer, and then samples were frozen and stored at −20 °C. ISSR-PCR analyses were performed [68], and for PCR, we used 5 ng/μL DNA samples. Each PCR amplification was performed in a 48 10 μL reaction mixture (1 μL 10 × Taq MgCl2 buffer, 1.2 μL 25 mM MgCl2, 1 μL 2 mM dNTPs, 0.4 μL primer, 4.32 μL deionized H2O, 0.08 μL Taq DNA polymerase, and 2 μL of club moss DNA). ISSR reactions were carried out in an Eppendorf thermo cycler. The reaction program was set to 94 °C for 7 min, 32 cycles of 94 °C for 30 s, 55/39/46 °C for 5 min, 72 °C 2 min, and a final extension of 72 °C for 7 min. Fourteen primers were used, and seven were chosen for further analysis: ISSR-B, ISSR-C, ISSR I-28, ISSR I-39a, ISSR I-50a, WARD-2, and ARCADE-3. Sequences of these primers are provided in Appendix A, Table A1.
The amplification products were analyzed using electrophoresis in a 1.5% TBE agarose gel. The gel was run for about 3.5 h. Gel results were registered using a BioDocAnalyse documentation system (Biometra, Göttingen, Germany). GeneRulerTM DNA Ladder Mix (100–10,000 bp) was used as standard.

Supplementary Materials

The following supporting information can be downloaded at https://doi.org/10.5281/zenodo.12802736.

Author Contributions

Conceptualization, R.R.-V., A.V. and S.J.; methodology, R.R.-V., A.V., S.J. and J.P.; software, visualization, and data curation, A.V. and S.J.; field research, R.R.-V., A.V., J.F. and S.J.; molecular analysis, J.F. and J.P.; writing—original draft preparation, R.R.-V.; writing—review and editing, R.R.-V. and S.J. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the Research Council of Lithuania, National Science Program: ‘Agro-forest and water ecosystems sustainability’, Grant No. SIT-1/2015.

Data Availability Statement

Data is contained within the article or Supplementary Materials.

Conflicts of Interest

The authors declare no conflicts of interest.

Appendix A

Table A1. DNA polymorphism of club moss species determined by ISSR analysis.
Table A1. DNA polymorphism of club moss species determined by ISSR analysis.
PrimerSequence
5′→3′
Annealing Temperature, °CAnalyzed
Bands
Polymorphism,
%
DNA Band Length,
bp
LALCLALCLALC
ISSR B(AG)8CG516750.0057.14700–1150600–1250
ISSR C(AG)8TG4910750.0057.14620–1500600–1500
ISSR I-28(GT)6CG39232378.2669.57480–2000580–2100
ISSR I-39a(AGC)4AC42222122.7338.1480–2100580–2500
ISSR I-50aCCA(GCT)447262338.4621.74630–2700460–3000
WARD 2(AC)8G47221877.2766.67400–1800580–2000
ARCADE 3(CA)8TG49222245.4527.27400–2000400–1700
In total 131121
Average 51.74 ± 7.5948.23 ± 7.23
LA—L. annotinum, LC—L. clavatum.
Table A2. Summary of published growth data of two Lycopodium species.
Table A2. Summary of published growth data of two Lycopodium species.
MeasurementResult, cmReference
Lycopodium annotinum L.
Total growth/year44.2Oinonen [19]
Growth up front/year22.1[19]
Average and range of rhizome growth/year32 (22–45)[69]
Average number and range of aerial shoots/year7.9 (7–9)[69]
Root growth/year0.7[70]
Length of horizontal module with 1 root9.2[70]
Mean annual extension/year6.3[71]
Mean annual extension/year20
Lycopodium clavatum L.
Total growth/year50.2[19]
Growth up front/year25.1[19]
Average and range of rhizome growth/year70 (48–103)[69]
Average number and range of aerial shoots/year12.8 (12–14)[69]
Table A3. Growth dynamics of young sporophytes of Lycopodium annotinum in Vaccinium myrtillosum Scots pine stands near Beržupis rivulet, Varėna District *.
Table A3. Growth dynamics of young sporophytes of Lycopodium annotinum in Vaccinium myrtillosum Scots pine stands near Beržupis rivulet, Varėna District *.
Measurement1985198719891991
Number of sprouts23252116
New sprouts-93-
Died off sprouts-775
Average length of horizontal modules, cm8 ± 1.512.5 ± 1.523.3 ± 2.533.7 ± 3.4
Maximum length of horizontal modules, cm112760113
* Adapted from [38] Naujalis (1995).
Table A4. Lycopodium clavatum occurrence in Vaccinio-myrtilosa forests in Lithuania *.
Table A4. Lycopodium clavatum occurrence in Vaccinio-myrtilosa forests in Lithuania *.
Age GroupNo. of Plots EvaluatedOverall Coverage, %
26–3513Up to 20
36–4519Up to 20
46–5515Up to 20
66–758Up to 20
76–8520Up to 20
86–100641–60
101–12010Up to 20
121–140461–80
141>541–60
* Adapted from [56]. Karazija (2003). No L. clavatum sporophytes detected in pine stands younger than 26 years.

References

  1. Eriksson, O. Dynamics of genets in clonal plants. Trends Ecol. Evol. 1993, 8, 313–316. [Google Scholar] [CrossRef] [PubMed]
  2. Honnay, O.; Bossuyt, B. Prolonged clonal growth: Escape route or route to extinction? Oikos 2005, 108, 427–432. [Google Scholar] [CrossRef]
  3. Grime, J.P. Plant Strategies and Vegetation Processes, 2nd ed.; John Wiley & Son: Chichester, UK, 2001; pp. 138–177. [Google Scholar]
  4. Eriksson, O. Seedling dynamics and life histories in clonal plants. Oikos 1989, 55, 231–238. [Google Scholar] [CrossRef]
  5. Callaghan, T.V.; Headley, A.D.; Svensson, B.M.; Lixian, L.; Lee, J.A.; Lindley, D.K. Modular growth and function in the vascular cryptogam Lycopodium annotinum. Proc. R. Soc. Lond. Ser. B 1986, 228, 195–206. [Google Scholar] [CrossRef]
  6. Wittig, R.; Jungman, R.; Ballach, H.J. The extent of clonality in large stands of Lycopodium annotinum L. Flora 2007, 202, 98–105. [Google Scholar] [CrossRef]
  7. Svensson, B.M.; Rydin, H.; Carlsson, B.A. Clonal plants in the community. In Vegetation Ecology, 2nd ed.; van der Maarel, E., Franklin, J., Eds.; Wiley-Blackwell: Oxford, UK, 2013; pp. 141–163. [Google Scholar]
  8. Bruchmann, H. Über die Prothallien und die Keimpflanzen mehrerer europäischer Lycopodien, und zwar über die von Lycopodium clavatum, L. annotinum, L. complanatum und L. selago; Perthes: Gotha, Germany, 1898. [Google Scholar]
  9. Pressel, S.; Bidartondo, M.I.; Field, K.J.; Rimington, W.R.; Duckett, J.G. Pteridophyte fungal associations: Current knowledge and future perspectives. J. Syst. Evol. 2016, 54, 666–678. [Google Scholar] [CrossRef]
  10. Perez-Lamarque, B.; Laurent-Webb, L.; Bourceret, A.; Maillet, L.; Bik, F.; Cartier, D.; Labolle, F.; Holveck, P.; Epp, D.; Selosse, M.A. Fungal microbiomes associated with Lycopodiaceae during ecological succession. Environ. Microbiol. Rep. 2023, 15, 109–118. [Google Scholar] [CrossRef]
  11. Headley, A.D.; Callaghan, T.V.; Lee, J.A. Water uptake and movement in the clonal plants, Lycopodium annotinum L. and Diphasiastrum complanatum (L.) Holub. New Phytol. 1988, 110, 497–502. [Google Scholar] [CrossRef]
  12. Newell, S.Y. Estimating fungal biomass and productivity in decomposing litter. In The Fungal Community: Its Organization and Role in the Ecosystem, 2nd ed.; Carroll, G.C., Wicklow, D.T., Eds.; Marcel Dekker: New York, NY, USA, 1992; pp. 521–561. [Google Scholar]
  13. Landa, K.; Benner, B.; Watson, M.A.; Gartner, J. Physiological integration for carbon in mayapple (Podophyllum peltatum), a clonal perennial herb. Oikos 1992, 63, 348–356. [Google Scholar] [CrossRef]
  14. Bainard, J.D.; Henry, T.A.; Bainard, L.D.; Newmaster, S.G. DNA content variation in monilophytes and lycophytes: Large genomes that are not endopolyploidy. Chromosome Res. 2011, 19, 763–775. [Google Scholar] [CrossRef]
  15. Marchant, D.B.; Chen, G.; Cai, S.; Chen, F.; Schafran, P.; Jenkins, J.; Shu, S.; Plott, C.; Webber, J.; Lovell, J.T.; et al. Dynamic genome evolution in a model fern. Nat. Plants 2022, 8, 1038–1051. [Google Scholar] [CrossRef]
  16. Eriksson, O. Evolution of seed dispersal and recruitment in clonal plants. Oikos 1992, 63, 439–448. [Google Scholar] [CrossRef]
  17. Rimgailė-Voicik, R.; Naujalis, J.R. Techniques for locating and analyzing subterranean Lycopodium and Diphasiastrum gametophytes in the field. Appl. Plant Sci. 2022, 10, e11458. [Google Scholar] [CrossRef] [PubMed]
  18. Rimgailė-Voicik, R.; Naujalis, J.R. Presence of juvenile club moss (Lycopodiaceae) sporophytes and gametophytes in relation to vegetation cover in dry pine forests. Am. Fern J. 2016, 106, 242–257. [Google Scholar] [CrossRef]
  19. Oinonen, E. Lycopodium clavatum L. ja L. annotinum L.-kasvustojen laajuus rinnastettuna samanpaikkaisiin L. complanatum L.-ja Pteridium aquilinum (L.) Kuhn—Esiintymiiin sekä puuston ikään j paloaikoihin. Acta For. Fenn. 1968, 87, 5–53. [Google Scholar]
  20. Watson, M.A.; Casper, B.B. Morphogenetic constraints on patterns of carbon distribution in plants. Annu. Rev. Ecol. Evol. Syst. 1984, 15, 233–258. [Google Scholar] [CrossRef]
  21. Yang, X.; Poelmans, W.; Grones, C.; Lakehal, A.; Pevernagie, J.; Van Bel, M.; Njo, M.; Xu, L.; Nelissen, H.; De Rybel, B.; et al. Spatial transcriptomics of a lycophyte root sheds light on root evolution. Curr. Biol. 2023, 33, 4069–4084. [Google Scholar] [CrossRef]
  22. Hetherington, A.J.; Dolan, L. The evolution of lycopsid rooting structures: Conservatism and disparity. New Phytol. 2017, 215, 538–544. [Google Scholar] [CrossRef]
  23. Sutherland, W.J.; Stillman, R.A. The foraging tactics of plants. Oikos 1988, 52, 239–244. [Google Scholar] [CrossRef]
  24. Horn, K.; Franke, T.; Unterseher, M.; Schnittler, M.; Beenken, L. Morphological and molecular analyses of fungal endophytes of achlorophyllous gametophytes of Diphasiastrum alpinum (Lycopodiaceae). Am. J. Bot. 2013, 11, 2158–2174. [Google Scholar] [CrossRef]
  25. Eames, A.J. Illustrations of some Lycopodium gametophytes. Am. Fern J. 1942, 32, 1–12. [Google Scholar] [CrossRef]
  26. Chytrý, M.; Tichý, L.; Hennekens, S.M.; Knollová, I.; Janssen, J.A.M.; Rodwell, J.S.; Peterka, T.; Marceno, C.; Landucci, F.; Danihelka, J.; et al. EUNIS Habitat Classification: Expert system, characteristic species combinations and distribution maps of European habitats. Appl. Veg. Sci. 2020, 23, 648–675. [Google Scholar] [CrossRef]
  27. European Council. Council Directive 92/43/EEC of 21 May 1992 on the conservation of natural habitats and of wild fauna and flora. Off. J. Eur. Union 1992, 206, 50. [Google Scholar]
  28. IUCN. Lycopodium annotinum. 2017. Available online: https://www.iucnredlist.org/species/18821015/85433421#assessment-information (accessed on 6 June 2024).
  29. Commission Regulation (EC) No 318/2008 of 31 March 2008, Amending Council Regulation (EC) No 338/97 on the Protection of Species of Wild Fauna and Flora by Regulating Trade Therein. Available online: https://eur-lex.europa.eu/LexUriServ/LexUriServ.do?uri=OJ:L:2008:095:0003:0062:EN:PDF (accessed on 6 June 2024).
  30. Braun-Blanquet, J. Pflanzensoziologie: Grundzüge der Vegetationskunde; Springer: Berlin/Heidelber, Germany, 1928. [Google Scholar]
  31. Ellenberg, H. Zeigerwerte der Gefaesspflanzen Mitteleuropas. Scripta Geobotanica IX; Verlag Erich Goltz: Göttingen, Germany, 1974. [Google Scholar]
  32. Tichý, L.; Axmanová, T.; Dengler, J.; Guarino, R.; Jansen, F.; Midolo, G.; Nobis, M.P.; Van Meerbeek, K.; Acic, S.; Attore, F.; et al. Ellenberg-type indicator values for European vascular plant species. J. Veg. Sci. 2023, 34, e13168. [Google Scholar] [CrossRef]
  33. Pec, G.J.; Karst, J.; Sywenky, A.N.; Cigan, P.W.; Erbilgin, N.; Simard, S.W.; Cahill, J.F., Jr. Rapid Increases in Forest Understory Diversity and Productivity following a Mountain Pine Beetle (Dendroctonus ponderosae) Outbreak in Pine Forests. PLoS ONE 2015, 10, e0124691. [Google Scholar] [CrossRef] [PubMed]
  34. Gilliam, F.S. The ecological significance of the herbaceous layer in temperate forest ecosystems. BioScience 2007, 57, 845–858. [Google Scholar] [CrossRef]
  35. Cherubini, P.; Battipaglia, G.; Innes, J.L. Tree vitality and forest health: Can tree-ring stable isotopes be used as indicators? Curr. For. Rep. 2021, 7, 69–80. [Google Scholar] [CrossRef]
  36. Hooper, D.U.; Adair, E.C.; Cardinale, B.J.; Byrnes, J.E.K.; Hungate, B.A.; Matulich, K.L.; Gonzalez, A.; Duffy, J.E.; Gamfeldt, L.; O’Connoret, M.I. A global synthesis reveals biodiversity loss as a major driver of ecosystem change. Nature 2012, 486, 105–108. [Google Scholar] [CrossRef] [PubMed]
  37. Lithuanian Forest Cutting Rules. Available online: https://e-seimas.lrs.lt/portal/legalAct/lt/TAD/TAIS.364764/asr (accessed on 6 June 2024).
  38. Naujalis, J.R. Sporiniai Induočiai Kaip Augalų Bendrijų Komponentai; Baltic ECO: Vilnius, Lithuania, 1995; p. 295. [Google Scholar]
  39. Oinonen, E. The time table of vegetative spreading in oak fern (Carpogymnia dryopteris (L.) Löve & Löve) and may-lily (Maianthemum bifolium (L.) F. W. Schmidt in Southern Finland. Acta For. Fenn. 1971, 118, 4–37. [Google Scholar]
  40. Oinonen, E. Keltalieon (Lycopodium complanatum L.) itiöllinen uudistuminen Etelä-Suomessa kloonien laajuutta ja ikää koskevan tutkimuksen valossa. Acta For. Fenn. 1967, 83, 3–85. [Google Scholar] [CrossRef]
  41. Nemchinova, A.V.; Krinycyn, I.G.; Smirnova, L.A. Biometricheskye osobennosti biomorf plauna bulavovidnovo (Lycopodium clavatum L.) v ontogeneze. Vestnik KGU Im. N. A. Nekrasova 2014, 4, 45–48. [Google Scholar]
  42. Svensson, B.M.; Floderus, B.; Callaghan, T.V. Lycopodium annotinum and light quality: Growth responses under canopies of two Vaccinium species. Folia Geobot. 1994, 29, 159–166. [Google Scholar] [CrossRef]
  43. Moola, F.M.; Vasseur, L. Recovery of late-seral vascular plants in a chronosequence of post-clearcut forest stands in coastal Nova Scotia, Canada. Plant Ecol. 2004, 172, 183–197. [Google Scholar] [CrossRef]
  44. Harper, J.L. Population Biology of Plants; Academic Press: London, UK, 1977; p. 892. [Google Scholar]
  45. Levin, D.A.; Crepet, W.L. Genetic variation in Lycopodium lucidulum: A phylogenetic relic. Evolution 1973, 27, 622–632. [Google Scholar] [PubMed]
  46. Ellstrand, N.C.; Roose, M.L. Patterns of Genotypic Diversity in Clonal Plant Species. Am. J. Bot. 1987, 74, 123–131. [Google Scholar] [CrossRef]
  47. Reisch, C.; Poschlod, P. Clonal diversity and subpopulation structure in central European relict populations of Saxifraga paniculata Mill. (Saxifragaceae). Feddes Repert. 2004, 115, 239–247. [Google Scholar] [CrossRef]
  48. Ranker, T.A.; Geiger, J.M.O. Population genetics. In Biology and Evolution of Ferns and Lycophytes; Ranker, T.A., Haufler, C.H., Eds.; Cambridge University Press: Cambridge, UK, 2008; pp. 75–103. [Google Scholar]
  49. Schnittler, M.; Horn, K.; Kaufmann, R.; Rimgailė-Voicik, R.; Klahr, A.; Bog, M.; Fuchs, J.; Bennert, H.W. Genetic diversity and hybrid formation in Central European club-mosses (Diphasiastrum, Lycopodiaceae)—New insights from cp microsatellites, two nuclear markers and AFLP. Mol. Phylogenet. Evol. 2019, 131, 181–192. [Google Scholar] [CrossRef] [PubMed]
  50. Shorohova, E.; Kneeshaw, D.; Kuuluvainen, T.; Gauthier, S. Variability and dynamics of old-growth forests in the circumboreal zone: Implications for conservation, restoration and management. Silva Fenn. 2011, 45, 785–806. [Google Scholar] [CrossRef]
  51. State Forest Service. Forest Statistic. 2021. Available online: https://amvmt.lrv.lt/lt/atviri-duomenys-1/misku-statistikos-leidiniai/misku-ukio-statistika/2021-m-1 (accessed on 6 June 2024).
  52. Marozas, V.; Sasnauskienė, J. Changes of ground vegetation after shelter wood cuttings in pine forests, the hemiboreal zone, Lithuania. Balt. For. 2021, 27, 72–79. [Google Scholar] [CrossRef]
  53. Herben, T.; Chytrý, M.; Klimešová, J. A quest for species-level indicator values for disturbance. J. Veg. Sci. 2016, 27, 628–636. [Google Scholar] [CrossRef]
  54. Sokołowski, A.W.; Paluch, R. The effect of clear cutting and artificial regeneration with pine on species composition of fresh coniferous stands in the Białowieża Primeval Forest. Sylwan 2006, 150, 12–19. [Google Scholar]
  55. Obidziński, A. Disturbance as an element of forest dynamics. Sylwan 2001, 145, 51–59. [Google Scholar]
  56. Karazija, S. Age related dynamics of pine forest communities in Lithuania. Balt. For. 2003, 9, 50–62. [Google Scholar]
  57. Jogiste, K.; Frelich, L.E.; Laarmann, D.; Vodde, F.; Baders, E.; Donis, J.; Jansons, A.; Kangur, A.; Korjus, H.; Koster, K.; et al. Imprints of management history on hemiboreal forest ecosystems in the Baltic States. Ecosphere 2018, 9, e02503. [Google Scholar] [CrossRef]
  58. National Land Service under the Ministry of Agriculture of the Republic of Lithuania. Lithuanian National Atlas; Ministry of Agriculture of the Republic of Lithuania: Vilnius, Lithuania, 2016; Volume 1.
  59. State Forest Service under the Lithuanian Ministry of Environment. Forest Cadastre Data. 2023. Available online: https://www.geoportal.lt/map#portalAction=openService&serviceUrl=https%3A%2F%2Fwww.geoportal.lt%2Fmapproxy%2Fvmt_mkd%2FmapServer (accessed on 5 May 2024).
  60. Nobis, M.; Hunziker, U. Automatic thresholding for hemispherical canopy photographs based on edge detection. Agric. For. Meteorol. 2005, 128, 243–250. [Google Scholar] [CrossRef]
  61. Frazer, G.W.; Canham, C.D.; Lertzman, K.P. Gap Light Analyzer (GLA): Imaging Software to Extract Canopy Structure and Gap Light Transmission Indices from True-Colour Fish Eye Photographs; User’s Manual and Program Documentation; Simon Fraser University: Burnaby, BC, Canada; Institute of Ecosystem Studies: Millbrook, NY, USA, 1999; Available online: https://rem-main.rem.sfu.ca/downloads/Forestry/GLAV2UsersManual.pdf (accessed on 6 June 2024).
  62. Kembel, S.W.; Cowan, P.D.; Helmus, M.R.; Cornwell, W.K. Picante: R tools for integrating phylogenies and ecology. Bioinformatics 2010, 11, 1463–1464. [Google Scholar] [CrossRef] [PubMed]
  63. Nei, M. Estimation of average heterozigosity and genetic distance from a small number of individuals. Genetics 1978, 89, 583–590. [Google Scholar] [CrossRef] [PubMed]
  64. Peakall, R.; Smouse, P. GenAlEx v.6: Genetic analysis in Excel. Population genetic software for teaching and research—An update. Bioinformatics 2006, 78, 265–285. [Google Scholar] [CrossRef] [PubMed]
  65. McKinney, W. Data Structures for Statistical Computing in Python. In Proceedings of the 9th Python in Science Conference, Austin, TX, USA, 28 June–3 July 2010. [Google Scholar]
  66. Vallat, R. Pingouin: Statistics in Python. J. Open Source Softw. 2018, 3, 1026. [Google Scholar] [CrossRef]
  67. Doyle, J.J.; Doyle, J.L. Isolation of plant DNA from fresh tissue. Focus 1990, 12, 13–15. [Google Scholar]
  68. Patamsytė, J.; Čėsnienė, T.; Naugžemys, D.; Kleizaitė, V.; Vaitkūnienė, V.; Rančelis, V.; Žvingila, D. Genetic diversity of warty cabbage (Bunias orientalis L.) revealed by RAPD and ISSR markers. Zemdirbyste 1990, 98, 29–300. [Google Scholar]
  69. Primack, R.B. Growth patterns of five species of Lycopodium. Am. Fern J. 1973, 63, 3–7. [Google Scholar] [CrossRef]
  70. Callaghan, T.V.; Svensson, B.M.; Headley, A.D. Modular growth of Lycopodium annotinum. Fern Gaz. 1986, 13, 65–76. [Google Scholar]
  71. Svensson, B.M. Studies of the Metapopulation Dynamics of Lycopodium annotinum and Its Microenvironment. Ph.D. Thesis, University of Lund, Lund, Sweden, 1987. [Google Scholar]
Figure 1. Lycopodium annotinum and Lycopodium clavatum forming an intermixed outgrowth in Verkiai Regional Park, Vilnius city, Lithuania. Photo by Greta Valvonytė.
Figure 1. Lycopodium annotinum and Lycopodium clavatum forming an intermixed outgrowth in Verkiai Regional Park, Vilnius city, Lithuania. Photo by Greta Valvonytė.
Plants 13 02120 g001
Figure 2. Scatter plot of forest age and club moss outgrowth diameter. L. annotinumLycopodium annotinum; L. clavatumLycopodium clavatum. Primary clones, marked with dots, are circular or partially circular and most likely originated from gametophytes. Secondary clones, marked with diamonds, are differently shaped, and their origin is unclear.
Figure 2. Scatter plot of forest age and club moss outgrowth diameter. L. annotinumLycopodium annotinum; L. clavatumLycopodium clavatum. Primary clones, marked with dots, are circular or partially circular and most likely originated from gametophytes. Secondary clones, marked with diamonds, are differently shaped, and their origin is unclear.
Plants 13 02120 g002
Figure 3. Tree canopy shading in Nemenčinė and Veržuva sites for L. annotinum and L. clavatum sporophytes. Statistically significant p-values, indicated by asterisks (*), and non-significant p-values, indicated by ‘n.s.’ (not significant), were determined using a two-way ANOVA.
Figure 3. Tree canopy shading in Nemenčinė and Veržuva sites for L. annotinum and L. clavatum sporophytes. Statistically significant p-values, indicated by asterisks (*), and non-significant p-values, indicated by ‘n.s.’ (not significant), were determined using a two-way ANOVA.
Plants 13 02120 g003
Figure 4. Principal coordinates analysis (PCoA) reflecting differentiation among sporophytes of Lycopodium annotinum and L. clavatum, using ISSR method, in Veržuva and Nemenčinė. Results of PERMANOVA test using Nei genetic distances for each species provided. L. annotinum n = 14 in Veržuva and n = 15 in Nemenčinė. L. clavatum n = 13 in Veržuva and n = 10 in Nemenčinė.
Figure 4. Principal coordinates analysis (PCoA) reflecting differentiation among sporophytes of Lycopodium annotinum and L. clavatum, using ISSR method, in Veržuva and Nemenčinė. Results of PERMANOVA test using Nei genetic distances for each species provided. L. annotinum n = 14 in Veržuva and n = 15 in Nemenčinė. L. clavatum n = 13 in Veržuva and n = 10 in Nemenčinė.
Plants 13 02120 g004
Figure 5. Distribution of club moss sporophyte sampling sites and Scots pine-dominated forest cover in Lithuania.
Figure 5. Distribution of club moss sporophyte sampling sites and Scots pine-dominated forest cover in Lithuania.
Plants 13 02120 g005
Table 1. Growth data summary of two Lycopodium species.
Table 1. Growth data summary of two Lycopodium species.
Measurement, m2Result
Lycopodium annotinum L. (n = 30)
Average vegetative shoot cover23.9 ± 11.4%
Min vegetative shoot cover5%
Max vegetative shoot cover50%
Average dried-out shoot cover10.1 ± 8.3%
Average number of strobili10.9 ± 10.3
Total maximum number of strobili2218
Lycopodium clavatum L. (n = 23)
Average vegetative shoot cover17.9 ± 11.3%
Min vegetative shoot cover6%
Max vegetative shoot cover50%
Average dried-out shoot cover10%
Average number of strobili7 ± 7.1
Total maximum number of strobili1077
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

Rimgailė-Voicik, R.; Voicikas, A.; Fediajevaitė, J.; Juzėnas, S.; Patamsytė, J. Origin and Persistence of Lycopodium clavatum and Lycopodium annotinum (Lycopodiaceae) in Scots Pine Forests. Plants 2024, 13, 2120. https://doi.org/10.3390/plants13152120

AMA Style

Rimgailė-Voicik R, Voicikas A, Fediajevaitė J, Juzėnas S, Patamsytė J. Origin and Persistence of Lycopodium clavatum and Lycopodium annotinum (Lycopodiaceae) in Scots Pine Forests. Plants. 2024; 13(15):2120. https://doi.org/10.3390/plants13152120

Chicago/Turabian Style

Rimgailė-Voicik, Radvilė, Aleksandras Voicikas, Julija Fediajevaitė, Sigitas Juzėnas, and Jolanta Patamsytė. 2024. "Origin and Persistence of Lycopodium clavatum and Lycopodium annotinum (Lycopodiaceae) in Scots Pine Forests" Plants 13, no. 15: 2120. https://doi.org/10.3390/plants13152120

APA Style

Rimgailė-Voicik, R., Voicikas, A., Fediajevaitė, J., Juzėnas, S., & Patamsytė, J. (2024). Origin and Persistence of Lycopodium clavatum and Lycopodium annotinum (Lycopodiaceae) in Scots Pine Forests. Plants, 13(15), 2120. https://doi.org/10.3390/plants13152120

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