Skip to main content
Advertisement
Browse Subject Areas
?

Click through the PLOS taxonomy to find articles in your field.

For more information about PLOS Subject Areas, click here.

  • Loading metrics

Population dynamics of Phaius flavus in southeast China: Reproductive strategies and plants conservation

  • Jun Li,

    Roles Investigation, Writing – original draft

    Affiliations College of Landscape Architecture, Fujian Agriculture and Forestry University, Fuzhou, Fujian Province, China, National Forestry and Grassland Administration for Orchid Conservation and Utilization, Fujian Agriculture and Forestry University, Fuzhou, Fujian Province, China

  • Ya-ting Zhu,

    Roles Investigation

    Affiliations College of Landscape Architecture, Fujian Agriculture and Forestry University, Fuzhou, Fujian Province, China, National Forestry and Grassland Administration for Orchid Conservation and Utilization, Fujian Agriculture and Forestry University, Fuzhou, Fujian Province, China

  • Lun-yan Chen,

    Roles Data curation

    Affiliations College of Landscape Architecture, Fujian Agriculture and Forestry University, Fuzhou, Fujian Province, China, National Forestry and Grassland Administration for Orchid Conservation and Utilization, Fujian Agriculture and Forestry University, Fuzhou, Fujian Province, China

  • Ai-xian Lu,

    Roles Investigation

    Affiliations College of Landscape Architecture, Fujian Agriculture and Forestry University, Fuzhou, Fujian Province, China, National Forestry and Grassland Administration for Orchid Conservation and Utilization, Fujian Agriculture and Forestry University, Fuzhou, Fujian Province, China

  • Hong-yu Ji,

    Roles Software

    Affiliations College of Landscape Architecture, Fujian Agriculture and Forestry University, Fuzhou, Fujian Province, China, National Forestry and Grassland Administration for Orchid Conservation and Utilization, Fujian Agriculture and Forestry University, Fuzhou, Fujian Province, China

  • Hai-ping Liu,

    Roles Investigation

    Affiliations College of Landscape Architecture, Fujian Agriculture and Forestry University, Fuzhou, Fujian Province, China, National Forestry and Grassland Administration for Orchid Conservation and Utilization, Fujian Agriculture and Forestry University, Fuzhou, Fujian Province, China

  • Ze-xin Li,

    Roles Investigation

    Affiliations College of Landscape Architecture, Fujian Agriculture and Forestry University, Fuzhou, Fujian Province, China, National Forestry and Grassland Administration for Orchid Conservation and Utilization, Fujian Agriculture and Forestry University, Fuzhou, Fujian Province, China

  • Zuo-dong Lin,

    Roles Investigation

    Affiliations College of Landscape Architecture, Fujian Agriculture and Forestry University, Fuzhou, Fujian Province, China, National Forestry and Grassland Administration for Orchid Conservation and Utilization, Fujian Agriculture and Forestry University, Fuzhou, Fujian Province, China

  • Sha-sha Wu,

    Roles Conceptualization

    Affiliations College of Landscape Architecture, Fujian Agriculture and Forestry University, Fuzhou, Fujian Province, China, National Forestry and Grassland Administration for Orchid Conservation and Utilization, Fujian Agriculture and Forestry University, Fuzhou, Fujian Province, China

  • Jun-wen Zhai

    Roles Conceptualization, Funding acquisition, Investigation, Project administration

    zhai-jw@163.com

    Affiliations College of Landscape Architecture, Fujian Agriculture and Forestry University, Fuzhou, Fujian Province, China, National Forestry and Grassland Administration for Orchid Conservation and Utilization, Fujian Agriculture and Forestry University, Fuzhou, Fujian Province, China

Abstract

Because of species diversity and troubling conservation status in the wild, Orchidaceae has been one of the taxa with most concern in population ecological research for a long time. Although Orchidaceae is a group with high adaptability, they have become endangered for complex and various reasons such as the germination? difficulty and habitat loss, which makes it difficult to develop an accurate protection strategy. Phaius flavus is a terrestrial orchid which used to be widely distributed in central and southern Asia; however, large populations are difficult to find in the wild. Thus, the aim of this study was to provide a new perspective for conserving endangered P. flavus by investigating the mechanisms of its population decline; we established time-specific life and fertility tables, age pyramids, survival curves, and mortality curves for this plant and then conducted Leslie matrix model. We found that both of the populations from Wuhu Mount (WM) and Luohan Mount (LM) showed declining trends and exhibited pot-shaped age pyramids, low net reproductive rates, and negative intrinsic growth rates. The population from the Beikengding Mount (BM) showed a stable status with a bell-shaped age pyramid. However, it has a significant risk of decline because of the low net reproductive rate and intrinsic growth rate. This study use time-specific life and fertility tables, age pyramids, survival curves, and mortality curves, showed that the population decline of P. flavus could be attributed to 1) the shortage of seedlings caused by the low germination rate in the wild and 2) the loss of adult individuals caused by anthropogenic disturbances. To protect this species from extinction in these areas, we suggest that human activities in these habitats should be strictly forbidden and ex situ conservation of this plant in botanical gardens is also necessary.

1. Introduction

Organisms adjust their life history traits and reproductive strategies to adapt to the environmental conditions of their locations [1]. Plants face a variety of risks during their lifetimes, such as competition, predation, diseases and pests. To ensure their reproduction, species in different habitats have evolved a variety of population adaptation strategies [2]. The study of population dynamics can reflect the characteristics of quantitative variation, living status, and habitat fitness [3]. Gui et al. studied a population of Phoebe bournei in central Jiangxi Province and, based on survival curves and quantitative analysis, found its population was capable of self-development but was sensitive to external disturbances [4]. Zhu et al. analyzed the quantity characteristics of Ginkgo biloba by investigating population inventories and distributions, and argued that this plant had a long period what? under severe fragmentation of its distribution, and suggested improving its habitat quality and strengthening protections [5]. Xiao et al. studied a population of Michelia wilsonii, which is a kind of PSESP (Plants Species with Extremely Small Populations), in Ya-an and found a sharp decrease of the population in the early age period, stability in the middle age period, and a decline in the late age period [6]. The above studies show that quantitative analysis is a reliable reference to understand the living status of plant populations, and the findings would be helpful for future preservation strategy development.

Orchidaceae is one of the largest and most diverse flowering plant families and consists of approximately 800 genera and 27,800 species [7]. Orchids live in nearly every terrestrial ecosystem in the world except for polar regions and arid regions such as deserts, thus they are usually considered to be highly adaptive to various environments and particularly diversified in tropical areas [8]. However, the survival of wild orchids faces challenges due to various factors [9]. Thus, as a flagship group in plant conservation, orchids have been a hot topic in endangered plant conservation research for a long time. First, habitat loss is the most lethal variable for explaining extinction probabilities [10]. Human activities (including forest clear-cutting and the expansion of agricultural landscapes) have exacerbated the natural fragmentation of landscapes [11]. Habitat fragmentation decreases plant population sizes and increases spatial isolation, which hampers the exchange of seeds and pollen between fragmented populations [12]. Second, the physiological adaption of orchids makes them sensitive to changes in temperature and rainfall. Especially for perennial orchids, climate change may affect different parts of the life cycle, including germination [13], growth [14], flowering [15, 16], and survival [17]. Third, orchid survival has close associations with the biotic environment. Orchids have complicated pollination systems with specialized pollinators, and their survival strongly depends on pollination success [16, 18]. At the same time, the entire Orchidaceae family appears to consist of mycorrhizal “cheaters,” as orchids have short-circuited the reciprocal plant-fungal exchange [19]. The distribution and availability of suitable mycorrhiza usually determine the seeding recruitment of orchids, and successful seeding establishment is a critical life story stage [20]. However, we lack information necessary for the conservation of Orchidaceae till now, especially for species that are known to be endangered [16].

Phaius flavus (Bl.) Lindl. is a tall, tuberous, perennial orchid that is mainly distributed in Asia, including China, Japan, India and most Southeast Asian countries [21]. It grows on damp slopes under trees at altitudes from 300 to 2,500 m and has a long history of cultivation as an ornamental and medicinal plant in China. As scaled artificial cultivation systems have not been established and extensively applied, wild populations of P. flavus have been excessively harvested over a long period as commercial demand has increased. Meanwhile, habitat fragmentation caused by human activities has made it increasingly difficult for P. flavus populations to survive in the wild. As a result, it is now difficult to find large populations in the field. Wang et al. studied the tissue culture and aseptic germination of P. flavus [22]. They provided a technology to mass produce P. flavus seedlings, which has potential application prospects for the protection and development of P. flavus. Kawakami et al. isolated a flexuous virus from cultivated P. flavus in Japan and analyzed its complete nucleotide sequence, and the results showed that its open reading frames differed from those of previously reported potex viruses [23]. The flora of Chinese orchids are distinguished by rich diversity in geographical types, especially in the broad subtropical area between the Qinling Mountains and the Tropic of Cancer [24]. In China, the main distribution areas of P. flavus are in the subtropical monsoon climate zone, which is characterized by high temperatures and abundant rainfall in summer with warm and dry winters. It is a representative of wild P. flavus populations in Fujian because it is located in the southeast part of China and is in a typical zone of the subtropical monsoon climate.

Here, population dynamics and individual number predictions were studied in three populations of P. flavus in Fujian Province based on an investigation of its biological characteristics. The following questions were addressed: 1) What is the survival status for wild P. flavus? If the answer is not optimistic, what is the possible reason for this situation? 2) Is there a recruitment situation and development tendency of the three populations studied here? Additionally, if there is a gap of data between three populations, how did the gap come about? By solving these questions, we discuss the growth and reproduction of P. flavus and provide a theoretical basis for the conservation of orchid germplasm resources.

2 Materials and methods

2.1 Study sites

This study was conducted on three natural populations of P. flavus located in Fujian, China, namely, the Wuhu Mount population (WM population), the Luohan Mount population (LM population) and the Beikengding Mount population (BM population). The WM population is present in the Sandiejing National Forest Park in the western part of Fuzhou, southern subtropical region, while the LM and BM populations are present in Sanming City, central subtropical region [25]. Other information for these sites is shown in Table 1. The dominant species in these study areas are Sloanea sinensis, Castanopsis sclerophylla, Schima superba, Ficus erecta, Aspidistra elatior, Oplismenus undulatifolius, Colysis elliptica and Gardneria multiflora.

2.2 Investigation of biological characteristics

Three transects were established along three valleys with the distribution of P. flavus. Based on a complete survey along these transects, we defined one site each in Wuhu Mount and Luohan Mount and two in Beikengding Mount. Because the distance between the two study sites in Beikengding Mount was less than 1 km and there were no habitat breaks or barriers between individuals, all P. flavus plants at that location were collectively placed into the BM population. Each transect consisted of adjacent 50 cm × 50 cm plots. A total of 146 P. flavus plants along three transects were mapped and labeled, and a grid method was used to obtain their relative coordinates. Pseudobulbs linked by flashy rhizomes were counted as a single individual plant and genotype. These pseudobulbs consisted of leafed pseudobulbs and leafless pseudobulbs, and their numbers were recorded, separately. Because every individual produces a new pseudobulb each year, we determine the age of each individual based on the number of pseudobulbs, this method is recorded in cultivation references and some previous studies. This regularity is shared by Trias verrucose [26], Cymbidium sinense [8], Dendrobium sinense [27], Calanthe tsoongiana [28], Phaius flavus and many other orchids. The above studies provide more ideas for wild orchid conservation, and also provide a reference for our study. To confirm this rule, we counted the number of buds germinated by every individual between Feb. 2019 and Apr. 2019 and further analyzed these results. The numbers of infructescences and persistent fruits were recorded.

2.3 Data analysis

2.3.1 Establishment of time-specific life tables.

Pseudobulbs of P. flavus plants form a single-chain structure during the growth process. We used an approach of spatial sequences instead temporal sequences [8, 29, 30], and the specific indicators and their data processing are as follows. We defined X as the age of a sample in years, ax as the number of plants surviving at the start of the age interval x, and lx as the proportion of plants surviving at the start of age interval x.

(1)

We then defined dx as the proportion of plants that died during the age interval from x to x + 1, qx as the death rate during the age interval from x to x + 1, and Kx as the rate of disappearance.

(2)(3)(4)

Lx is the average number of plants that survived during the age interval x and x + 1, Tx is the total number of plants that survived from age interval x, and ex is the life expectancy at the start of age interval x.

(5)(6)(7)

2.3.2 Age pyramids, survival curves and mortality curves.

Age pyramids were plotted with lx as the horizontal axis and the corresponding age level (X) as the vertical axis. The survival curve was plotted with ln lx as the horizontal axis and the corresponding age level (X) as the vertical axis. The mortality curve was plotted with the qx as the horizontal axis and the corresponding age level (X) as the vertical axis.

2.3.3 Fertility table establishment.

Based on data from the time-specific life tables, specific indicators for the fertility tables and their data processing are as follows. We defined X as the age of a sample in years, lx as the proportion of plants surviving at the start of age interval x, mx as the average number of progeny produced by individuals in X-age, R0 as the net reproductive rate, rm as the intrinsic rate of increase, λ as the increment rate per unit, and Tx as the average period from the birth of parents to the birth of progeny.

(8)(9)(10)(11)

2.3.4 Leslie matrix construction.

The Leslie matrix is one of the best-known methods for describing the growth of populations; it is used in ecology to model changes in populations of organisms over a period of time [31]. Px (the sum of survival rates from age interval x) was calculated on the basis of survival rate (lx) from the specific-time life table and fx (the average number of progeny that were produced by X-age adults and survived to the next age) was calculated by using data from the fertility table.

(12)(13)

Then, the population massive and age structure after one unit time interval were determined from the data of the previous time interval. Nt is the population quantity in the t time interval and N0 is the population quantity at the time of this survey.

(14)

M is the population projection matrix.

(15)

3 Results and analysis

3.1 Biological characteristics of P. flavus plants

P. flavus plants grow in gullies and/or near creeks with deep moist humus soil, good drainage and sunlight-scattered forest slopes or stone surfaces. Their oval pseudobulbs are linked by short, fleshy stems for which pseudobulbs are 2 to 8 cm tall and 1 to 6 cm in diameter. Pseudobulbs produce 1–6 leaves at the apex, and they remain green and functional for 1–14 years; at the same time, pseudobulbs that have lost all their leaves still retain the possibility of sprouting new buds to resume their vegetative growth and reproductive process. The leaf lengths of different individuals in the population varied widely, and the adults were stable at approximately 48 cm. New pseudobulbs usually germinated yearly from the bottom of the last generation of pseudobulbs in the middle of March. The persistence of leafless pseudobulbs and the annual proliferation of new pseudobulbs worked together to form the chain-shaped P. flavus plants.

According to our survey, the average numbers of buds that each individual generated in Feb. 2019 were 0.87 (WM population), 0.91 (LM population) and 0.83 (BM population). Based on these data, we estimate that each individual produced nearly one additional pseudobulb per year, and the number of pseudobulbs is further assumed to be this individual’s age grade. According to the statistics of persistent flower stems, P. flavus recruitments take at least 4 years from germination to flowering. Once these adult individuals enter the sexual reproduction stage, they remain fertile until their maximum lifespan has been reached. The sums of the number of persistent fruits in the 3 populations were 3 (WM population), 11 (LM population) and 50 (BM population), and the average numbers of fruits per individual in the 3 populations were 0.05 (WM population), 0.26 (LM population) and 2.63 (BM population) (Fig 1).

thumbnail
Fig 1. The habitat and growth characteristics of Phaius flavus.

A. Habitat; B. pseudobulbs and fleshy stems; C. persistent dehiscent fruits; D. a budding bud.

https://doi.org/10.1371/journal.pone.0272929.g001

3.2 Time-specific life tables

The time-specific life tables for the 3 populations are shown in Tables 24. With 62 individuals and 18 age classes, the WM population is the largest but youngest population among these three, while the BM population is the smallest and oldest and has 41 individuals and 39 age classes. Abruption phenomena frequently appear for age classes over 11 and cause many formulas in the tables to have errors in the counting process. This is why there are many blanks in the tables. Especially for the BM population, there are only 9 age classes that have surviving individuals between age classes 12–39, which account for 32.1% of the total.

thumbnail
Table 2. Time-specific life table of P. flavus WM population.

https://doi.org/10.1371/journal.pone.0272929.t002

thumbnail
Table 3. Time-specific life table of P. flavus LM population.

https://doi.org/10.1371/journal.pone.0272929.t003

thumbnail
Table 4. Time-specific life table of P. flavus BM population.

https://doi.org/10.1371/journal.pone.0272929.t004

As shown in Tables 13, negative mortality (qx) is common and similar in the three populations, and the frequency of appearance was 50% (WM population), 41.7% (LM population) and 52.6% (BM population). Viewed from age class divisions, negative mortality was widely distributed across every age stage but is more obviously distributed in the seedling and aging stages and less so in the middle stage, especially when the plants are just reaching puberty.

3.3 Age pyramid, survival curve and mortality curve

3.3.1 Age pyramid.

The WM and LM populations have narrower age spans and similar age structures, and as Figs 2 and 3 show, their age pyramids are spindly. The BM population has a wider age span and more significant abruption, and as shown in Fig 4, its age pyramid is roughly bell-shaped and is widest at the bottom and diminishes upward. Therefore, the WM and LM populations are decreasing populations with scarce seedlings. Although most individuals in this kind of population are fertile, it is dangerous to continue lacking seedling recruitment. The BM population has equal numbers of adolescents and adults, which is characteristic of a stable population. In addition, the individuals in these populations were unevenly distributed across age classes, which means the occurrence of unusual deaths. This conclusion is similar to what we found in the time-specific life tables.

3.3.2 Survival curve and mortality curve.

By analyzing the survival curves (Fig 5) and mortality curves (Fig 6) of three P. flavus populations, the following conclusions can be drawn: 1) As there is a higher survival rate of seedlings (e.g., individuals in the 1–4 age class) in the BM population than in the other two populations, recruitments in the BM population have a higher possibility of entering the breeding period; 2) The peak survival rate appears in the 5–8 age class, and plants in this stage have just begun to flower, while their highest survival rate causes them to become mainstays of population quantities; 3) In the aging stage, the survival curve of the BM population declines steadily, while other survival curves decline steeply, which causes the quantities of the WM and LM populations to be more dependent on their middle-aged individuals.

According to Fig 6, the mortality of P. flavus plants as a whole increases with increased age, and the peak appears later in life. From the perspective of single age classes, there are significant differences in mortality between one age class and another, which cause the mortality curve to fluctuate widely. Compared with the BM population, the WM and LM populations have higher death rates in the seedling and aging stages, which means that the living status for these two populations are worse than that for the BM population, which is similar to what we found in the survival curves.

Taken together, wild P. flavus populations currently face a lack of recruitment and difficulty in seedling survival, and their individuals in the middle age class have high living quality and large quantities and bear the burden of sustaining the populations’ existence. Thus, for the WM and LM populations, sharp decreases in seedling numbers cause decreases in their long-term stability; although the mortality curve and survival curve of the BM population showed jagged changes, the BM population still has higher stability, as its overall trend matches the Deevey Ⅰ curve.

3.4 Fertility table parameter analysis

r0 (net reproductive rate) and rm (intrinsic rate of increase) are calculated from the data in Tables 57. The r0 values of the three populations are 0.075 (WM population), 0.081 (LM population) and 1.24 (BM population). The rm values are -0.336 (WM population), -0.392 (LM population) and 0.07 (BM population).

Parameter r0 is the increase multiple per generation, and r0 values below 1 mean that the increase quantity is smaller than the decrease quantity and this population cannot achieve growth in numbers and is in danger of dying out; r0 values greater than 1 mean that this population is showing an upward trend. The r0 values for the WM and LM populations are far less than 1, which means that, in these two populations, the instantaneous birth rate is far less than the instantaneous death rate and the tendency for increase is very weak compared with tendency for decrease. Although the r0 of the BM population is greater than 1, it is close to 1. Therefore, the BM population can slowly maintain its growth under present conditions.

Parameter rm is similar to r0, but its critical value is 0. As we found in the analysis of parameter r0, the WM and LM populations show downward trends, while the BM population is stable.

3.5 Leslie matrix and dynamic quantitative prediction model construction

Based on the investigation of biological characteristics, the average number of seedlings produced by fertile individuals in each population was calculated as parameter mx (average fertility), and the population density was then simulated separately with a matrix model that combined the time-specific life table and fertility table.

Overall, the three populations in this study are declining populations. Population numbers of the WM and LM populations will sharply decrease according to Tables 9 and 11. The continuous reduction of seedling stage plants and, at the same time, the loss of adults aggravates this situation by reducing seed sources and causing a vicious circle until the populations finally vanish. The BM population is different from the two populations mentioned above; its numbers will increase in the mid-early term and decrease later, as shown in Table 13. Before N20, high survival rates of mid-aged individuals become the cornerstone of population numbers; due to the stable supply of recruits offered by fertile plants, the numbers of the BM population will even increase slightly in a few years. The turning point will appear after N20 when the high quality of living of middle stage-plants cannot offset the downward tendency caused by low fertility, and the BM population may begin to decrease gradually until disappearing completely.

3.5.1 WM population.

mx for the WM population is 0.02. The Leslie matrix model and dynamic quantitative prediction model are shown in Tables 8 and 9. The formula for the dynamic quantitative prediction model is Nt = Nt-1- Nt-1e-1.56, where Nt is the number of individuals at time t. As predicted, the WM population will vanish after 18 years.

thumbnail
Table 9. Dynamic quantitative prediction model of P. flavus WM population.

https://doi.org/10.1371/journal.pone.0272929.t009

3.5.2 LM population.

mx of the LM population is 0.02. The Leslie matrix model and dynamic quantitative prediction model are shown in Tables 10 and 11. The formula for the dynamic quantitative prediction model is Nt = Nt-1- Nt-1e0.75, where Nt is the number of individuals at t time. As predicted, the LM population will vanish after 22 years.

thumbnail
Table 10. Leslie matrix model of P. flavus LM population.

https://doi.org/10.1371/journal.pone.0272929.t010

thumbnail
Table 11. Dynamic quantitative prediction model of P. flavus LM population.

https://doi.org/10.1371/journal.pone.0272929.t011

3.5.3 BM population.

mx of the BM population is 0.24. The Leslie matrix model and dynamic quantitative prediction model are shown in Tables 12 and 13. The formula of the dynamic quantitative prediction model is Nt = Nt-1- Nt-1e0.71, where Nt is the number of individuals at t time. As predicted, the BM population will vanish after 40 years.

thumbnail
Table 12. Leslie matrix model of P. flavus BM population.

https://doi.org/10.1371/journal.pone.0272929.t012

thumbnail
Table 13. Dynamic quantitative prediction model of P. flavus BM population.

https://doi.org/10.1371/journal.pone.0272929.t013

4 Discussion

4.1 The formation of a negative mortality phenomenon

In the analysis of time-specific life tables, we found that negative mortality (qx) was common in every period of each population. This is unusual. As a kind of ecological study method that works backward from spatial sequence to temporal sequence, it is usually used to study long-lived species (their complete life history is hard to monitor) whose ages can be accurately evaluated. This is why specific-time life tables are generally used for tree research [32]. P. flavus plants have age indicators, and the age of the oldest plants in three populations was 39 years, so this method applies to P. flavus. In recent years, wild populations of Cypripedium lentiginosum [30], Cymbidium sinense [8], Pleione formosana [29], Dendrobium sinense [27] and other orchids have been studied by life table methods. In the research of D.sinense, the negative mortality issue was raised and was explained by “unusual deaths”. CHEN used data-smoothing techniques to correct the negative death rate problem in his research on Cinnamomum micranthum [33]. This technique is widely recognized and has been applied in some studies in recent years [34, 35]. Throughout the studies that used data-smoothing techniques, their objects of study not only had far greater quantities than this study but also had more balanced distributions of age classes; these are probably the characteristics that made them suitable for this technique. This method is more applicable for repairing errors that do not affect overall trends, so this study was excluded from its application. We posed two hypotheses regarding the high frequency of negative mortality: 1) seedling numbers have fluctuated and 2) P. flavus plants have been overcollected.

4.1.1 Speculation on the cause of recruitment fluctuations.

The WM and LM populations only produced 1 recruitment last year, while the BM population had none. This is far below the maximum number of survivors in each population. P. flavus plants use both sexual and asexual reproduction as a kind of clonal plant. Clonal plants produce ramets by breaking the connecting structures between the ramet and genet [36]. Most of the time, one P. flavus individual will produce one or more pseudobulb per year but these potential ramets have difficulty separating from their genets because the short fleshy stems that link them together are hard to break or rot under natural conditions. From another perspective, the long-term existence of connecting structures may decrease ramet mortality by nutrition storage and absorption of resources [37]. Therefore, it is shown in this study that asexual reproduction of P. flavus plants has more value for the extension of mother plants, stabilization of population numbers, and accumulation of dry matter, rather than the production of new individuals. The increase in seedling numbers mainly depends on sexual reproduction.

According to our survey of persistent fruits, there are significant differences between quantities, and each population had ripe and dehiscent fruits. This evidence suggests that P. flavus plants can complete their sexual reproduction process without any apparent difficulties. Additionally, Wang et al. studied the aseptic germination of P. flavus plants, and they proved that the seed germination rate could be above 70% and that the rooting rate could be above 90% without special treatment under experimental conditions [22]. It can be inferred from this report that P. flavus seeds have no internal structural barriers to germination when compared with other orchids that do have difficulty germinating, such as Paphiopedilum. Our study supports the notion that germination of P. flavus seeds is highly dependent on environmental suitability and on the availability of mycorrhiza, which means that the seedling complements of P. flavus populations in the wild are easily affected by environmental change, thus causing fluctuations in recruitment quantities.

4.1.2 The effect of wild harvesting.

Wild harvesting of orchids, such as P. flavus, is assumed to be unsustainable because their life story characteristics cause many orchids to be vulnerable to disturbances [38]. Particularly in the flowering period, their large, bright, yellow flowers may attract more attention from passers-by and this increases the risk of individual blossoms been harvested. This not only directly causes the loss of adult individuals but also indirectly causes the loss of potential seedlings. Lin et al. believe that maintaining human interference at a moderate level can provide more opportunities for other plants to settle in the community [39]. A similar conclusion was drawn by Niu Li-qin in her essay, and she believes that moderate human disturbance is beneficial for increasing plant species richness [40]. Both of their studies are based on the idea that species that have been disturbed can still maintain their growth in population numbers, while the overall trend of P. flavus wild populations has been falling. Therefore, protecting these populations should come first, as wild harvesting still does more harm than good.

4.2 Development gaps between three P. flavus populations

Compared with the BM population, the population structure of the LM population is more similar to that of the WM population in the southern subtropical region. The LM and BM populations are both in the middle subtropical region. Therefore, the development of differences between these three populations is not entirely dependent on natural environmental factors but also influenced by internal factors and human interference.

On the one hand, nutrient accumulations by plants are important for reproductive success. The phenomenon that fruit and seedling production in orchids differs between different populations within species has also been found in other research [41]. Compared with other populations, the BM population is significantly higher than the others in the number of age classes and average age, which means that plants in the BM population made fuller matter preparation. In the statistics of persistent fruits, we found that the output due to sexual reproduction in the BM population was highest in the three populations, as the fruit average of the BM population (2.63) was far greater than that for the others (e.g., 0.05 for the WM population and 0.26 for the LM population). Therefore, for the case of slight differences in other natural environmental factors, it is speculated that the sexual reproduction activity of P. flavus is likely to be positively correlated with the degree of material accumulation. On the other hand, the WM population is near a footpath in Sandiejing Forest Park, so its habitat highly overlaps with areas of human activity; the distance between the LM population and a local tea garden is less than 1 km; the BM population is located in an inaccessible valley. Based on accessibility analysis, the WM population is the most likely to be disturbed by humans, the LM population is second most likely, and the BM population is least likely to be disturbed.

4.3 Reproduction strategies of P. flavus

In a study of Trias verrucosa, LIU et al. found that this orchid has a reproduction strategy of strengthening asexual reproduction and weakening sexual reproduction as a response to the low rewards of sexual reproduction [26]. In their matrix analyses, Miguel Franco and Jonathan Silvertown also suggested that seed production and seedling recruitment are secondary to the long-term viability of populations of long-lived plant species, including orchids [42]. We assume that a similar phenomenon occurs in P. flavus populations. Although there is no barrier to its sexual reproduction, the development focus of young populations may shift to asexual reproduction to ensure the survival of those populations. The adult stages of orchids are profoundly different from the juvenile stages. Jacquemyn et al. suggested that adult orchids generally exhibit much lower annual mortality rates and, in many species, have the potential to achieve long lifespans [43]. Therefore, investment in asexual reproduction appears to be more worthwhile for young populations to use less energy and gain more profits [43]. After long-term accumulation of nutrients, mature populations are more competitive in their surroundings, which gives them the ability to continue surviving while pursuing sexual reproduction to seek more potential opportunities for future development and expansion [26]. Zhang suggested that most resources should be allocated toward sexual reproduction in habitats with fluctuating environmental conditions and strong competition [44]. Obviously, the development of this reproduction strategy is related to the relatively stable growing environment of P. flavus populations.

4.4 Protective and restoring measures

As we found in this study, the survival prognosis for wild P. flavus populations is poor; the WM and LM populations are at risk of dying out, and the BM population can barely maintain stable numbers. Based on the above analysis, the fundamental cause of the decline of P. flavus wild populations lies in recruitment shortages caused by low fertility, while the loss of adults triggered by human harvesting aggravates this difficult position and speeds up the decrease in the number of individuals. The seedling shortage issue is probably due to a combination of P. flavus own reproduction strategies and the germination difficulties that are commonly seen in orchids. Therefore, ensuring survival of adults is a critical part of promoting the recovery and renewal of P. flavus populations. According to the analysis of endangering mechanisms, it is suggested that human activities around P. flavus populations be strictly limited to reduce over-collection and habitat breakage. These populations should be revisited every year to confirm their status so that we can conduct ex situ conservation in time.

References

  1. 1. Willems JH, Ellers J. Plant performance and population characteristics of Orchis simia (Orchidaceae) in two extremes of its distribution area. FLORA. 1996;191(1):41–8.
  2. 2. Zhong ZC. The Reproduction Strategies of Plants. Chinese Journal of Ecology. 1995(01):37–42.
  3. 3. Kang HJ, Chen ZL, Liu P, Hao ZY, Wei FM. Population structure and distribution pattern of Emmenopterys henryi in Dapanshan National Nature Reserve. Acta Ecolohica Sinica. 2007(01):389–96.
  4. 4. Gui YK, Pan P, Ouyang XZ, Zang H, Guo R, Li Y. Quantitative Characteristics and Distribution Pattern of Phoebe Bournei Natural Population in Central Jiangxi Province. Acta Ecologica Sinica. 2019;38(10):2918–24.
  5. 5. Zhu W, Zhang GF, Gui ZY. Population Structure and Dynamics of Ancient Ginkgo Biloba in Nantong, Jiangsu Province. Acta Ecologica Sinica. 2019;12(38):3609–16.
  6. 6. Xiao SL, Fu MY, Yang K, Chen XH. Population Structure and Quantitative Dynamics of an Extremely Small Population, Michelia wilsonii. Acta Botanica Boreali-Occidentalia Sinica. 2019;39(07):1279–88.
  7. 7. Christenhusz M, Byng JW. The number of known plants species in the world and its annual increase. PHYTOTAXA. 2016; 261:201–17.
  8. 8. Liu ZJ, Chen LJ, Liu KW, Li LQ, Zhang YT, Huang LQ. Climate warming brings about extinction tendency in wild population of Cymbidium sinense. Acta Ecologica Sinica. 2009;29(7):3443–55.
  9. 9. Kull T, Hutchings MJ. A comparative analysis of decline in the distribution ranges of orchid species in Estonia and the United Kingdom. BIOL CONSERV. 2005;129(1):39.
  10. 10. Crain BJ, Tremblay RL, Rouget M. Do richness and rarity hotspots really matter for orchid conservation in light of anticipated habitat loss? Divers Distrib. 2014;20(6):652–62.
  11. 11. Bender DJ, Fahrig CL. Habitat Loss and Population Decline: A Meta-Analysis of the Patch Size Effect. ECOLOGY. 1998;79(2):517–33.
  12. 12. Meekers T, Honnay O. Effects of habitat fragmentation on the reproductive success of the nectar-producing orchid Gymnadenia conopsea and the nectarless Orchis mascula. Plant Ecol. 2011;212(11):1791–801.
  13. 13. Mondoni A, Rossi G, Orsenigo S, Probert RJ. Climate warming could shift the timing of seed germination in alpine plants. Ann Bot. 2012;110(1):155–64. pmid:22596094
  14. 14. Reina-Rodríguez GA M.JER, Llanos FAC, Otero JT. Spatial distribution of dry forest orchids in the Cauca River Valley and Dagua Canyon: Towards a conservation strategy to climate change. J Nat Conserv. 2016; 30:32–43.
  15. 15. Attila MV, Lyi JTK, Végvári Z, Sramkó G, Sulyok J, Barta Z. Pollination mode predicts phenological response to climate change in terrestrial orchids: a case study from central Europe. J Ecol. 2012;100(5):1141–52.
  16. 16. Zuzana S, Spyros T, Pavel K. Pollination Mechanisms are Driving Orchid Distribution in Space. Sci Rep-Uk. 2020;10(1).
  17. 17. Liu H, Feng CL, Luo YB, Chen BS, Wang ZS, Gu HY. Potential Challenges of Climate Change to Orchid Conservation in a Wild Orchid Hotspot in Southwestern China. Bot Rev. 2010;76(2):174–92.
  18. 18. Damo A, Hernández-Ramírez F, Riggi L, Verspoor R, Bertolini V, Walker ML, et al. Pollination of euglossinophylic epiphytic orchids in agroecosystems and forest fragments in southeast Mexico. Eur J Agron. 2012;1(2):5–14.
  19. 19. Taylor DL, Bruns TD. Population, habitat and genetic correlates of mycorrhizal specialization in the ‘cheating’ orchids Corallorhiza maculata and C. mertensiana. Mol Ecol. 1999;8(10):1719–32. pmid:10583834
  20. 20. Jacquemyn H, Brys R, Vandepitte K, Honnay O, Ruiz IR, Wiegand T. A spatially explicit analysis of seedling recruitment in the terrestrial orchid Orchis purpurea. New Phytol. 2007;176(2):448–59. pmid:17888122
  21. 21. Committee FOCE, Kai-yong L, Xin-Qi C, Yi-bo L. Flora of China: Orchidaceae.: Science Press; 1999. p. 259.
  22. 22. Wang LH, Zhu YQ, Jiang YL, Feng YC. Tissue Culture and Rapid Propagation of Phaius flavus (Bl.) Lindl. Plant Physiology Communications, 2007(05):899–900.
  23. 23. Kawakami K, Fuji S, Miyoshi K. Complete nucleotide sequence of a new potexvirus, “Phaius virus X”, isolated from Phaius flavus Lindl. Arch Virol. 2008; Vol.153(NO.3):527–31.
  24. 24. Luo YB, Jia JS, Wang CL. A general review of the conservation status of Chinese orchids. Biodiversity Science. 2003;11(1):70–7.
  25. 25. Zheng JY, Bian JJ, Ge QS, Hao ZX, Yin YH, Liao YM. The climate regionalization in China for 1981–2010 (in Chinese). Chinese Science Bulletin. 2013;58(30).
  26. 26. Liu ZJ, Chen LJ, Lei SP, Rao WH, Li LQ. The reproduction stratagy of Trias verrucosa from China. Acta Ecologica Sinica. 2007;27(11):4460–8.
  27. 27. Wang TX, Qi SJ, Song XQ, Meng QW, Yu XD. Correlation analysis between population dynamics and reproductive strategy of Dendrobium sinense. Journal of Tropical Biology. 2018;9(2):189–97.
  28. 28. Lian JJ. Reproductive Biology of Calanthe tsoongiana. Chinese Academy of Forestry. 2013.
  29. 29. Wu SS, Zhou YZ, Li SX, Ma XY, Peng DH. Population dynamics and reproductive modes of Pleione feomosana Hayata. Journal of Fujian College of Forestry. 2014;34(04):297–303.
  30. 30. Liu ZJ, Chen LJ, Rao WH, Li LQ, Zhang YT. Correlation between numeric dynamics and reproductive behaviour in Cypripedium lentiginosum. Acta Ecologica Sinica. 2008;28(1):111–21.
  31. 31. Miller FP, Vandome AF, Mcbrewster J. Leslie matrix. Germany: Alphascript Publishing; 2010.
  32. 32. Wu CZ, Hong W, Xie JS, Wu JL. Life table analysis of Tsuga longibracteata population. Chinese Journal of Applied Ecology. 2000;11(3):333–6. pmid:11767626
  33. 33. Chen YZ, Ma QX, Feng LZ, Huang YL, Zheng RQ. The population life table and periodic fluctuation of Cinnamomun micranthum, an endangered plant. Acta Ecologica Sinica. 2006;26(12):4267–72.
  34. 34. Wang DP, Peng SL, Chen FP, Ji SY. Population dynamics and considerations for the conservation of the rare Cycas fairylakea in China. For Ecosyst. 2012;14(2):118–23.
  35. 35. Liu ZQ, Dong HG, Liang QL, Chen WM. Life table and survival analysis of relict plant Malus Sieversii population. Journal of Arid Land Resources and Environment. 2016;30(10):98–103.
  36. 36. Dong M. Clonal growing in plants in relation to resource heterogeneity: foraging behavior. Acta Botanica Sinica. 1996;10(38):828–35.
  37. 37. Dong M. Clonal growing in plants in relation to resource heterogeneity: risk spreading. Acta Botanica Sinica. 1996;20(6):543–8.
  38. 38. Koopowitz H. Orchids and their conservation. London, UK: Batsford Press; 2001.
  39. 39. Lin JY, Tang FC, He YM, Li J, Liang RL. Effects of Human Distribance on Community Structure and Species Diversity of Phoebe bournei Community. Journal of West China Forestry Science. 2019;48(4):72–8.
  40. 40. Niu LQ. Effects of anthropogenic disturbance on forest community structure in Wutai Mountain, China. Chinese Journal of Applied and Environmental Biology. 2019;25(02):300–12.
  41. 41. Chapurlat E, Agren J, Sletvold N. Spatial variation in pollinator-mediated selection on phenology, floral display and spur length in the orchid Gymnadenia conopsea. New Phytol. 2015;4(208):1264–75. pmid:26183369
  42. 42. Franco M, Silvertown J. A comparative demography of plants based upon elasticities of vital rates. Ecology. 2004(85):531–8.
  43. 43. Jacquemyn H, Brys R, Jongejans E. Size-dependent flowering and costs of reproduction affect population dynamics in a tuberous perennial woodland orchid. Journal of Ecology. 2010;5(98):1204–15.
  44. 44. Zhang YF, Zhang DY. Asexual and sexual reproductive strategies in clonal plants. J Plant Ecol. 2006;30(1):174–83.