Skip to main content

Comparative genetic and epigenetic of the Sphagneticola trilobata (L.) Pruski from different regions in China



Sphagneticola trilobata (L.) Pruski is a prevalent and widely distributed invasive plant in South China. To investigate the molecular mechanisms underlying its rapid adaptation, we employed DNA methylation-sensitive amplified polymorphism (MSAP) and simple sequence repeat (SSR) analysis to study 60 S. trilobata individuals collected from Fuzhou (FZ), Haikou (HK), Jinghong (JH) and Guangzhou (GZ).


In this study, we computed the Shannon diversity index (I) of SSR and MSAP as 0.354 and 0.303, respectively. The UPGMA phylogenetic tree and PCoA analyses showed that MSAP had a better discriminatory power to distinguish populations from different regions. Notably, the GZ population was found to be the most distinct from the other three populations. Moreover, Mantel analysis revealed a significantly higher correlation between epigenetic distance and geographic distance as compared to genetic distance and geographic distance. Consequently, the correlation between epigenetic distance and geographic distance observed to be markedly stronger than that between genetic distance and geographical distance on Mantel analysis.


The S. trilobata populations in various regions displayed a high of complementary genetic and epigenetic diversity, which was a key feature contributing to their rapid invasion. Interestingly, the correlation between epigenetics and geographical distance was significantly stronger than that observed for genetics and geographical distance. These findings indicated that the epigenetic mechanism of S. trilobar exhibited high plasticity, leading to significant differences in methylation pattern across different populations.

Peer Review reports


With the globalization of the economy and trade, the rapid expansion in the scope and frequency of human activities has led to the ubiquity of invasive plants, resulting in significant economic losses and ecological damage on a global scale. Genetic and evolutionary processes are believed to be key determinants in the establishment and spread of invasive species. Invasive species must adapt to their new environment in order to survive and reproduce [1], and genetic variation plays an important role in plant adaptation [2]. Therefore, the ability of invasive species to rapidly evolve and adapt to new environmental conditions is a critical factor that determines their success in colonizing new areas and outcompeting native species. According to research, many invasive species exhibit a high level of genetic diversity, which is closely linked to their successful establishment in new environments [3]. For example, the South American weed Flaveria bidentis was introduced from Tianjin and quickly spread throughout in northern China. Its ability for sexual reproduction and multiple introductions have led to the maintained of high levels of genetic diversity, which in turn has facilitated local adaptation and invasion [4]. On the other hand, despite having very low genetic diversity, cloned invasive species are still capable of quickly adapting to new environments and becoming successful invaders. Therefore, the genome-level process that translates into phenotypic diversity needs to be re-evaluated. By studying successful plant invaders, we can gain a better understanding of this process.

Sphagneticola trilobata (L.) Pruski (synonyme: Wedelia trilobata, S. trilobata hereafter), belongs to the Asteraceae, a perennial herb native to tropical South America [5]. This plant has been listed as one of the 100 most invasive alien species in the world [6]. It was used widely in landscape greening because of its rapid reproduction, strong survival adaptability, and easy coverage. In the 1970s, S. trilobata was introduced to China as a ground cover plant, and S. trilobata quickly fled to the wild after settlement, and then developed and spread rapidly in southern China [7].Previous studies have found that S. trilobata has the ability to inhibit the growth and development of surrounding species, which poses a significant threat to local biodiversity and causes the imbalance of local ecological diversity. This is due to its rapid coverage of occupied open spaces through allelopathic therapy and asexual reproduction [8, 9]. Moreover, S. trilobata can invade other remote areas through water flow or human activities as its stem can directly develop into a complete plant [10, 11]. Currently, it is distributed in Guangdong, Guangxi, Fujian, Hainan, Yunnan and other provinces, mainly in parks, scenic spots, riversides, roadsides, fields, etc.

The extract from S. trilobata leaves has been shown to inhibit the germination of seeds and the growth of seedlings from both competitors and its own seeds, resulting in a high number of seed sterilizations and few seedlings in wild populations [12]. Consequently, S. trilobata primarily encroaches on new habitats through asexual clonal reproduction. In general, clonal reproduction typically leads to reduced gene flow between individuals within a species and a low level of genetic diversity, which is not conducive to adapting to changing environmental conditions. However, this was contradictory to the fact that S. trilobata was highly adaptable to diverse environments. Some researches have shown that epigenetics can provide clonal plants with more adaptive strategies. Invasive species tend to have higher phenotypic plasticity than non-invasive or native species, and even though invasive species populations in invasive areas were able to evolve greater plasticity than native populations [13]. For instance, the invasive plant Alternanthera philoxeroides, which predominated in clonal reproduction, responded rapidly to environmental stresses mainly through epigenetic variation [14]. Thus far, the molecular mechanisms of S. trilobata's environmental adaptations has been poorly reported. There was still unclearly about how S. trilobata itself adapted quickly in a short term. Therefore, it was hypothesized that S. trilobata might adapt quickly to the stresses of different environments in a similar manner.

Epigenetics provides clonal plants with more adaptive strategies. It is one of the best examples of how cloning can replace the slow mechanism of natural selection. Epigenetics depends to varying degrees on asexual reproduction and preserves or reverses gene ability to buffer plants against current challenges and future rapid environmental changes [15]. In addition, epigenetic modifications are genetically stable across generations [16], meaning that more favorable phenotypic variations can be transferred to future generations. Epigenetics refers to changes in gene expression levels based on covalent modifications on DNA or histones, rather than changes in the gene sequence, which are hereditary and inheritably reversible. Plants use epigenetic modifications of the genome to change the gene expression of individuals with the same genotype in response to changes in environmental conditions through phenotypic plasticity. The evolution of epigenetics may be driven by natural selection, indicating that in natural populations, epigenetics and genetic variation are two intertwined rather than independent evolutionary factors [17,18,19].

DNA methylation, RNA interference, and histone modification are widely recognized as essential components of epigenetic regulation. Among them, cytosine methylation stands out as the most important epigenetic regulation strategy and plays a central role in the epigenetic control of gene expression [20]. Furthermore, the hypothesis of epigenetic characteristics associated with invasion was largely based on that changes in DNA methylation were involved in gene expression regulation, particularly that methylation of gene promoters is associated with gene silencing [21]. It is worth noting that DNA methylation is one of the most common epigenetic modifications. On the contrary, in the contex of exploring genetic diversity research, due to the advantages of abundant multiple alleles, high polymorphism and dominance, simple sequence repeat (SSR) is a powerful tool for studying genetic diversity [22].

Since DNA methylation-sensitive amplified polymorphism (MSAP) does not require sequencing of a reference genome, it can be directly used to assess the methylation status of randomly distributed genomic cytosines. This method has the advantage of being easy to operate and inexpensive, and is currently the most widely used method for detecting cytosine methylation, which often occurs at the CG site in the promoter region of the DNA sequence [23]. MSAP is a method that can be used to directly assess the methylation status of randomly distributed genomic cytosines. Since it does not require sequencing of a reference genome, this method is widely used in the field of analytical epigenetics for detecting cytosine methylation changes in the ecological environment [24].

Previous studies on the invasion mechanism of S. trilobata have mainly focused on examining changes in phenotype or physiological indices in response to different environmental factors. However, there are few explorations based on genetics and epigenetics [25]. For better understanding the process of rapid adaptation of S. trilobata, this study explored the genetic and epigenetic structures and differences among populations of S. trilobata under different habitats by using molecular markers SSR and MSAP.

Materials and methods

Study species

Sphagneticola trilobata (L.) Pruski is a common invasive plant widely distributed in South China. It has been recorded in Guangdong, Guangxi, Fujian, Hainan and Yunnan province. It can adpapt to highly heterogeneous environments, such as parks, scenic spots, riversides, roadsides, fields, etc. S. trilobata, triangular-shaped leaves with toothed edges that are 5–10cm long and 2–5cm wide. The leaves are arranged alternately on the stem and are slightly hairy on both surfaces. The plant produces yellow daisy-like flowers with a diameter of approximately 2–3 cm, which bloom throughout the year.such as its triangular-shaped leaves and hairy surfaces. The plant's growth habit and morphology also vary depending on its environment, with plants growing in full sunlight typically being more compact and bushy than those growing in shade [10, 26].

Collection of meteorological data and materials

The current S. trilobata distribution records in China were obtained from literature reports and previous data collected by the laboratory. Maxent software and ArcGIS were utilized to predict contemporary suitable areas, and samples were gathered from highly suitable regions. As this species was introduced in the 1970s and has undergone a certain period of evolution, reference data from contemporary climate were chosen as the environmental factor data. The environmental data were collected from China Meteorological Administration.s website ( In this study, a total of 120 samples were collected solely from highly suitable areas in four different regions: Guangzhou, Guangdong Province (GZ), Haikou, Hainan Province (HK), Fuzhou, Fujian Province (FZ), and Jinghong, Yunnan Province (JH), respectively. Thirty samples were collected from each region (Table 1). The detailed sampling method we applied was as follows: 1) Three quadrats were randomly selected from each region, ensuring that each quadrat was more than 2 km apart (this was sample determined by the specific distribution location of S. trilobata in the area); 2) Ten samples were randomly selected from each quadrat and ensure each plant is more than 50 m apart; 3) Harvest the mature leaves after the third pair of leaves on the plant; 4) A total of 120 samples were collected from four populations, and the mature leaves were dried and numbered using silica gel from each region. DNA was extracted from each sample and stored at − 20 °C. Total genomic DNA extraction was extracted using the plant genome DNA extraction kit (Beijing Tiangen).

Table1 The information and environmental factors in sampling sites of Sphagneticola trilobata from different regions in China

SSR method

Based on the transcriptome data, six pairs of SSR primers were selected for capillary electrophoresis. The volume of Polymerase Chain Reaction (PCR) system was 10 μL, consisting of 1 μL of DNA template, 5 μL of 2 × Taq PCR Master Mix, 0.5 μL of upstream primers (the 5’ end of the upstream primer was labeled with phosphorescent dye during fluorescence PCR), 0.5 μL of downstream primers and 3 μL of ddH2O. The PCR conditions were set as follows: 95 ℃, 3 min; 34 cycles of (95 ℃, 30 s; The annealing temperature is the annealing temperature of each primer, 30 s; 72 ℃, 60 s); and a final extension at 72 ℃ for 5 min, and finally stored at -20 °C. In the end, PCR products were separated by 6% non-denaturing polyacrylamide gel electrophoresis [27].

MSAP method

Using EcoRI and two restriction sites 5′-CCGG-3′, the restriction endonucleases MspI and HpaII with different sensitivity to methylation were cut at 37 °C for 12 h. The resulting DNA digested products were connected and amplified by PCR. After 2% agarose gel detection, 8 pairs of primers were selected for PCR selective amplification. The selective PCR products were detected with 6% polypropylene gel, and then the color was observed by silver staining (detailed steps can be found in the supplementary materials). There were four cytosine methylation patterns in DNA (Table 2). Type I: sites with the presence of bands in both EcoRI/HpaII and EcoRI/MspI, namely (1.1), it represents CCGG site was unmethylated. Type II: bands present in EcoRI/HpaII, but absent in EcoRI/MspI, namely (1.0), indicating that the external cytosine site of only one DNA chain was methylated, which was also known as the Hemi-methylation of CCGG site. Type III: bands present in EcoRI/HpaII, but absent in EcoRI/MspI, namely (0.1), stands for the internal cytosine of the double DNA chain was methylated, which was also known as the complete methylation of CCGG site. Type IV: bands has not present in EcoRI/HpaII and EcoRI/MspI, namely (0.0). The markers for type IV mode (0.0) were treated as hypermethylated sites in this study followed [28,29,30].

Table 2 DNA methylation band types

Data analysis

The SSR bands of S. trilobata were quantified by numbers, where “1” denotes the presence of a band and “0” denotes the presence of no band. Similarly, the MSAP bands were counted by numbers, where “1” indicating DNA methylation, and “0” indicating non-methylation. The genetic diversity was calculated using Popgene software, and the phylogenetic trees of SSR with MSAP were constructed using MEGA7. NTSYS software was used to conduct evolutionary tree analysis; PCoA analysis was performed on the the resulting data. Furthermore, GenAlex software was used for mantel test analysis and AMOVA analysis. The R package (4.1.2) and Python (3.10.5) were utilized for data analysis and figure drawing.

Result and analysis

Genetic diversity and epigenetic genetic diversity analysis

The numbers of observed alleles (Na) in SSR ranged from 1.656 to 1.750, while the number of effective alleles (Ne) range from 1.319 to 1.530. The average values of Na and Ne were 1.703 and 1.406, respectively. In MSAP, Na ranges from 1.603 to 2.000, while the Ne ranges from 1.148 to 1.344, with the average values of 1.798 and 1.275, respectively. The average values of gene diversity (H) and the Shannon diversity index (I) of both SSR and MSAP were 0.236, 0.184 and 0.354, 0.303, respectively, with little variance (Table 3). In summary, compared with other populations, the genetic diversity of the GZ population was the highest, with slightly lower epigenetic diversity than that of the FZ population, while the genetic and epigenetic diversities of the HK population were significantly lower than those of the other three populations. In terms of genetic distance (Table 4), HK and FZ were the closest, while GZ and JH were the farthest. The epigenetics of JH and FZ were the closest, while the epigenetics of JH and HK were the farthest.

Table 3 Genetic and epigenetic diversity analysis on the four Sphagneticola trilobata populations from different invasive regions in china
Table 4 The distance of genetic and the distance of epigenetic genetic between different Sphagneticola trilobata populations in china

Cluster analysis of genetics and epigenetics of different populations

In Fig. 1, it can be observed that most samples from HK and GZ regions were clearly clustered into independent branches, while samples from JH and FZ regions were partially intertwined. The cluster of MSAP molecular markers displays a greater degree of independence compared to the SSR cluster in each region. The Principal Coordinate Analysis (PCoA) results reveal that the JH and GZ populations were closely related based on SSR, but they were separated from the FZ population (Fig. 2a). The PCoA findings was consistent with the results from UPGMA clustering. Based on MSAP, FZ, HK and JH populations were close, the second axis distinguishing these three populations from the GZ population (Fig. 2b).

Fig. 1
figure 1

UPGMA clustering analysis reveals differences between four populations of Sphagneticola trilobata in China based on SSR (a) and MSAP (b) analysis

Fig. 2
figure 2

Principal Coordinate Analysis (PCoA) showing the divergence in four Sphagneticola trilobata populations in China based on SSR (a) and MSAP (b) data

Correlation analysis of genetic and epigenetic distance and geographic distance for all populations

Mantel analysis of the relationship between geographical distance and genetic or epigenetic distance indicated that the R values were 0.0093 and 0.475, respectively, both with the 0.01 P values (Fig. 3a and b). In addition, the Mantel analysis examining the relationship between genetic and epigenetic distances yielded an R value of 0.2 (Fig. 3c), suggesting a stronger correlation in contrast to genetic distance.

Fig. 3
figure 3

The correlation between the genetic distance (a) and epigenetic distance (b) of all populations of Sphagneticola trilobata in China and geographical changes, as well as the correlation between the calculated genetic distance and epigenetic distance (c)

AMOVA analysis of SSR and MSAP

AMOVA analysis (Table 5) revealed that the epigenetic fixation index (Fst) of S. trilobata in different regions was 0.475, which was higher than the genetic fixation index (0.277). This suggested that the differentiation among populations was mainly attributed to epigenetics. Genetic and epigenetic gene flow (Nm) were 0.652 and 0.276, respectively, which were less than 1, indicating that gene flow among S. trilobata populations in different regions was limited by geographical distance.

Table 5 AMOVA analysis of genetic and epigenetic genetic between different Sphagneticola trilobata populations in China

Mode structure of epigenetic regions

The total methylated rates in the four populations (Fig. 4) were 23.82%, 18.25%, 32.92% and 27.55%, respectively; the fully methylated rates were 16.25%, 11.51%, 25.75% and 24.43%, respectively; the hemi-methylation rates were 7.57%, 6.75%, 7.16% and 3.12%, respectively; moreover, the non-methylated rates were 76.18%, 81.75%, 67.08% and 72.45%, respectively.

Fig. 4
figure 4

Comparation of Sphagneticola trilobata methylation status from different invasive areas in China Note: Each gray dot represents a sample collected from each region for MSAP detection, with a total of 15 samples collected from each region; Total methylated bands =  I + II + III + IV; Full-methylated bands = III + IV; Total methylated rate/% = [(II + III + IV)/(I + II + III + IV)] × 100%; Fully methylated rate(%) = [(III + IV)/(I + II + III + IV)] × 100%; Hemi-methylated rate (%) = [(II)/(I + II + III + IV)] × 100%; Non-methylated rate (%) = [(I)/(I + II + III + IV)] × 100%


Differences in genetic and epigenetic composition in S. trilobata rapid invasion

In this study, the Shannon Diversity Indexes (I) of genetic diversity and epigenetic diversity in S. trilobata were 0.472 and 0.440, respectively. In Alternanthera philoxeroides, the Shannon Diversity Indexes (I) of genetic diversity and epigenetic diversity were 0.002 and 0.046, respectively [14].These results indicated S. trilobata has relatively higher genetic and epigenetic diversity. The Shannon Diversity Index of genetic and epigenetic diversity in HK was low. The reason could be HK population was located in Hainan Province, which was an independent island, resulting in less frequent gene exchange between HK population with the other populations. As the earliest region in China to introduce S. trilobata, GZ has the longest time for the formation of its escape zone and natural community, which may be the reason have led to its highest genetic diversity. During the process of clonal reproduction, plants evolve complex adaptation with defense mechanisms through potential epigenetic traits due to the stress they feel from the environment, and are genetically mediated to produce phenotypic changes called phenotypic plasticity. Studies had shown that invasive species tend to had higher phenotypic plasticity than non-invasive or native species, even that populations of invasive species in invasive regions were able to evolve greater plasticity than native populations [13]. The FZ population had the highest epigenetic diversity, it may be due to the unique distribution of samples in FZ area. Farther distribution between different populations and denser the distribution within the population, would bring on greater epigenetic differentiation between different populations. However, the distribution of samples in other areas was more uniform. At the same time, GZ also has high epigenetic diversity, with only a very slight gap is lower than Fuzhou.

After analyzing the UPGMA phylogenetic tree and PCoA results, it was found that the MSAP had a pronounced effect than the SSR clustering. The epigenetic molecular markers seemed to be able to better distinguish regions, which indicated that epigenetics was more malleable during adaptation to different habitats. It was even more obvious in the PCoA result plot of MSAP that the GZ population was clearly separated from the rest of the region, which was a very interesting phenomenon. As the GZ region was the first to be introduced and cultivated in China, and its population has undergone a long-term stage of adaptive evolution, whereas other regions were introduced and cultivated at a later time. This observation also supports the notion that genetic and epigenetic diversity indices of species increase proportionally with the stage of adaptive evolution.This result was consistent with the findings of Shi et al., indicating that during clonal reproduction, plants evolved complex adaptations with defense mechanisms through underlying epigenetic traits due to environmental pressures, and improved adaptive ability through gene mediation [14]. Moreover, our research reveals a positive trend in the Shannon diversity index of genetic and epigenetic diversity in different populations. Previous studies had shown that incidental genetic variation can greatly shape epigenetic variation in order to truly enrich epigenetic mechanisms interact with genetics to promote heredity and ultimately lead to the transformation of genetic and other non-genetic mechanisms into body performance [31,32,33]. This also showed that in the invasion mechanism of S. trilobata, heredity and epigenetics were mutually supportive.

Comparison of the different Sphagneticola trilobata populations in genetic distances and epigenetic distance with their geographical

Mantel’s analysis revealed that the correlation between epigenetic distance and geographical distance was significantly stronger than that between genetic distance and geographical distance, indicating that epigenetic distance was more affected by geographical distance. Additionally, the Fst results suggested that epigenetics exhibited greater differentiation than genetics. Herrera et al. discovered that the differences in functional traits of Helleborus foetidus were related to environment and epigenetic distance [34]. Similarly, Pagel et al. observed a robust correlation between genetic distance and epigenetic distance in Linum catharticum from both dry and wet grasslands [35]. These findings provide evidence that environmental factors play a significant role in shaping genetic and epigenetic levels. In this study, we found that the genetic and epigenetic distances in S. trilobata just showed a significantly weakly correlated, where the habitats of our four regions were not in extreme opposition. Besides, Ne, H, I and altitude were significantly positively correlated in genetics, as well as PIC, Na, Ne, H, and I in epigenetics were all positively correlated with the highest temperature (Fig. S2).As a consequence, Adaptation pressures from various habitats are driving the genetic and epigenetic diversity in the species. Environmental factors had a certain degree of influence on the genetic and epigenetic levels, with a more significant effect on epigenetics.

Differences epigenetic composition in S. trilobata rapid invasion

The methylation rate in plants typically falls within the range of 20% to 40% [36]. In this study, the average methylation rate of S. trilobata in four regions ranged from 18.25%—32.92%. The results manifested that S. trilobata was significant differences in methylation among the four populations. The methylated DNA can cause chromatin extremely to become highly compact [37], resulting in reduce genetic exchange and a decrease in polymorphic bands. Furthermore, since type IV (0 0) was an unintuitive amplification site, the number of methylation sites detected by MSAP may be underestimated, implying that the actual level of methylation in S. trilobata may be higher than reported. The total methylated rate (%) and fully methylated rate (%) were both the highest in HK population, whereas those in GZ population were the lowest. The hemi-methylated rate (%) was the highest in FZ population. The non-methylated rate was the highest in GZ population and the lowest in HK population. Overall, the methylation level of GZ population was significantly lower than that of other populations, indicating that there were some differences in methylation patterns between different populations. This could be attributed to variations in local temperature and habitats at the sampling sites. In cloned plants, epigenetic modifications can created new modules over time and accumulate them for adaptations [15, 38]. There were abundant evidences suggesting that sexual and clonal reproduction can effectively superimposed, with epigenetic regulation participating in and regulating this clonal reproduction [39, 40].

For example, Shi et al. (2019) discovered a substantial amount of epigenetic differences within the clone population of Alternanthera philoxeroides in the same environment, which can be maintained until 10 generations before gradually disappearing, indicating that epigenetic variation was likely caused by environmental induction and spontaneous mutations [14]. Similarly, Guarino et al. (2019) denoted on the invading clonal population Arundo donax, which showed that epigenetic diversity was higher than genomic diversity, suggesting its invasion under different environmental conditions was attributed to higher epigenetic variability [41].

The differences in methylation patterns in different regions of S. trilobata indicated that the methylation patterns of cloned plants produced significant differences in invasion regions. These results suggest that epigenetic mechanisms partially compensate for the limited gene exchange between populations due to clonal reproduction. This indicates that epigenetics plays a supportive role in the process of adaptation and evolution, and exhibits greater plasticity in aiding clonal plants to adapt to invasion.


As can be concluded from this study, S. trilobata populations in different regions had high genetic and epigenetic diversity, which have played a mutually reinforcing role in facilitating its rapid invasion. The underlying drivers for these changes are attributed to diverse environmental pressures, which have been one of the primary factors responsible for S. trilobata's successful invasion. Notably, the correlation between epigenetics and geographic distance was significantly higher than that observed for genetic and geographic distance. By virtue of the significant difference in methylation patterns among different populations, the epigenetic mechanism of S. trilobata have endowed it with greater plasticity, thereby enabling it to adapt rapidly to changing environments. Overall, this study lays a foundation for a deeper understanding of how cloned plants undergo rapid adaptation.

Availability of data and materials

All relevant data can be found within the manuscript and its supporting materials. All data analyzed during this study are included in the supplementary information fles.


  1. Sakai AK, Allendorf FW, Holt JH, Lodge DM, Molofsky J, Baughman S, Cabin RJ, Cohen JE, Ellstrand NC, Mccauley DE. The population biology of invasive species. Annu Rev Ecol Syst. 2001;32:305–32.

    Article  Google Scholar 

  2. Henderson IR, Salt DE. Natural genetic variation and hybridization in plants. J Exp Bot. 2017;68(20):5415–7.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  3. Müller-Schärer H, Schaffner U, Steinger T. Evolution in invasive plants: implications for biological control. Trends Ecol Evol. 2004;19(8):417–22.

    Article  PubMed  Google Scholar 

  4. Ma J, Geng S, Wang S, Zhang G, Fu W, Shu B. Genetic diversity of the newly invasive weed Flaveria bidentis (Asteraceae) reveals consequences of its rapid range expansion in northern China. Weed Res. 2011;51(4):363–72.

    Article  Google Scholar 

  5. Zhang Q, Chen G, Huang J, Peng C. Comparison of the Ability to Control Water Loss in the Detached Leaves of Wedelia trilobata, Wedelia chinensis, and Their Hybrid. Plants (Basel). 2020;9(9):1227.

    Article  CAS  PubMed  Google Scholar 

  6. Invasive Species Specialist Group. 100 of the world’s worst invasive alien species. Auckland: Invasive Species Specialist Group; 2001.

  7. Song LY, Li CH, Peng SL. Elevated CO2 increases energy use efficiency of invasive Wedelia trilobata over its indigenous Congener. Biol Invasions. 2010;12(5):1221–30.

    Article  Google Scholar 

  8. Qi SS, Dai ZC, Miao SL, Zhai DL, Si CC, Huang P, Wang RP, Du DL. Light limitation and litter of an invasive clonal plant, Wedelia trilobata, inhibit its seedling recruitment. Ann Bot. 2014;114(2):425–33.

    Article  PubMed  PubMed Central  Google Scholar 

  9. Dai ZC, Wang XY, Qi SS, Cai HH, Sun JF, Huang P, Du DL. Effects of leaf litter on inter-specific competitive ability of the invasive plant Wedelia trilobata. Ecol Res. 2016;31(3):367–74.

    Article  Google Scholar 

  10. Sun Z, Chen Y, Schaefer V, Liang H, Li W, Huang S, Peng C. Responses of the Hybrid between Sphagneticola trilobata and Sphagneticola calendulacea to Low Temperature and Weak Light Characteristic in South China. Sci Rep. 2015;5:16906.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  11. Javed Q, Sun J, Azeem A, Jabran K, Du D. Competitive ability and plasticity of Wedelia trilobata (L.) under wetland hydrological variations. Sci Rep. 2020;10(1):9431.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  12. Azizan KA, Ibrahim S, Ghani NHA, Nawawi MF. Metabolomics approach to investigate phytotoxic effects of Wedelia trilobata leaves, litter and soil. Plant Biosystems. 2019;153(5):1–9.

    Article  Google Scholar 

  13. Davidson AM, Jennions M, Nicotra AB. Do invasive species show higher phenotypic plasticity than native species and if so, is it adaptive? Ecol Lett. 2011;14(4):419–31.

    Article  PubMed  Google Scholar 

  14. Shi W, Chen X, Gao L, Xu CY, Ou X, Bossdorf O, Yang J, Geng Y. Transient Stability of Epigenetic Population Differentiation in a Clonal Invader. Front Plant Sci. 2019;9:1851.

    Article  PubMed  PubMed Central  Google Scholar 

  15. Douhovnikoff V, Dodd RS. Epigenetics:a potential mechanism for clonal plant success. Plant Ecol. 2015;216(2):227–33.

    Article  Google Scholar 

  16. Gáspár B, Bossdorf O, Durka W. Structure, stability and ecological significance of natural epigenetic variation: a large-scale survey in Plantago lanceolata. New Phytol. 2019;221(3):1585–96.

    Article  PubMed  Google Scholar 

  17. Herrera CM, Medrano M, Bazaga P. Variation in DNA methylation transmissibility, genetic heterogeneity and fecundity-related traits in natural populations of the perennial herb Helleborus foetidus. Mol Ecol. 2014;23(5):1085–95.

    Article  CAS  PubMed  Google Scholar 

  18. Bossdorf O, Richards CL, Pigliucci M. Epigenetics for ecologists. Ecol Lett. 2008;11(2):106–15.

    PubMed  Google Scholar 

  19. Fujimoto R, Sasaki T, Ishikawa R, Osabe K, Kawanabe T, Dennis ES. Molecular mechanisms of epigenetic variation in plants. Int J Mol Sci. 2012;13(8):9900–22.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  20. Law JA, Jacobsen SE. Establishing, maintaining and modifying DNA methylation patterns in plants and animals. Nat Rev Genet. 2010;11(3):204–20.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  21. Paun O, Verhoeven K, Richards CL. Opportunities and limitations of reduced representation bisulfite sequencing in plant ecological epigenomics. New Phytol. 2019;221:738–42.

    Article  CAS  PubMed  Google Scholar 

  22. Powell W, Morgante M, Andre C, Hanafey M, Vogelet J, Tingey S, Rafalski A. The comparison of RFLP, RAPD, AFLP and SSR (microsatellite) markers for germplasm analysis. Mol Breed. 1996;2(3):225–38.

  23. Richards CL, Alonso C, Becker C, Bossdorf O, Bucher E, Colomé-Tatché M, Durka W, Engelhardt J, Gaspar B, Gogol-Döring A, Grosse I, van Gurp TP, Heer K, Kronholm I, Lampei C, Latzel V, Mirouze M, Opgenoorth L, Paun O, Prohaska SJ, Rensing SA, Stadler PF, Trucchi E, Ullrich K, Verhoeven KJF. Ecological plant epigenetics: Evidence from model and non-model species, and the way forward. Ecol Lett. 2017;20(12):1576–90.

    Article  PubMed  Google Scholar 

  24. Schrey AW, Mariano A, Foust CM, Kilvitis HJ, Lee JD, Liebl AL, Martin LB, Richards CL, Marta R. Ecological epigenetics:beyond MS-AFLP. Integr Comp Biol. 2013;53(2):340–50.

    Article  PubMed  Google Scholar 

  25. Huang P, Shen F, Abbas A, Wang H, Du Y, Du D, Hussain S, Javed T, Alamri S. Effects of Different Nitrogen Forms and Competitive Treatments on the Growth and Antioxidant System of Wedelia trilobata and Wedelia chinensis Under High Nitrogen Concentrations. Front Plant Sci. 2022;13:851099.

    Article  PubMed  PubMed Central  Google Scholar 

  26. Li Q. A new natural hybrid of Sphagneticola (Asteraceae, Heliantheae) from Guangdong, China. Phytotaxa. 2015;221(1):71.

    Article  Google Scholar 

  27. Liu T, Li Y, Duan W, Huang F, Hou X. Cold acclimation alters DNA methylation patterns and confers tolerance to heat and increases growth rate in Brassica rapa. J Exp Bot. 2017;68(5):1213–24.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  28. Gupta V, Bijo AJ, Kumar M, et al. Detection of epigenetic variations in the protoplast-derived germlings of ulvareticulata using methylation sensitive amplification polymorphism (MSAP). Marine Biotechnol. 2012;14(6):692–700.

    Article  CAS  Google Scholar 

  29. Schulz B, Eckstein RL, Durka W. Scoring and analysis of methylation-sensitive amplification polymorphismsfor epigenetic population studies. Mol Ecol Resour. 2013;13(4):642–53.

    Article  CAS  PubMed  Google Scholar 

  30. Ashikawa I. Surveying CpG methylation at 5’-CCGG in the genomes of rice cultivars. Plant Mol Biol. 2001;45(1):31–9.

    Article  CAS  PubMed  Google Scholar 

  31. Becker C, Hagmann J, Müller J, Koenig D, Stegle O, Borgwardt K, Weigel D. Spontaneous epigenetic variation in the Arabidopsis thaliana methylome. Nature. 2011;480:245–9.

    Article  CAS  PubMed  Google Scholar 

  32. Dubin MJ, Zhang P, Meng D, Remigereau MS, Osborne EJ, Paolo Casale F, Drewe P, Kahles A, Jean G, Vilhjálmsson B, Jagoda J, Irez S, Voronin V, Song Q, Long Q, Rätsch G, Stegle O, Clark RM, Nordborg M. DNA methylation in Arabidopsis has a genetic basis and shows evidence of local adaptation. Elife. 2015;4:e05255.

    Article  PubMed  PubMed Central  Google Scholar 

  33. Sasaki E, Kawakatsu T, Ecker JR, Nordborg M. Common alleles of CMT2 and NRPE1 are major determinants of CHH methylation variation in Arabidopsis thaliana. PLoS Genet. 2019;15(12):e1008492.

    Article  PubMed  PubMed Central  Google Scholar 

  34. Herrera CM, Medrano M, Bazaga P. Comparative epigenetic and genetic spatial structure of the perennial herb Helleborus foetidus: Isolation by environment, isolation by distance, and functional trait divergence. Am J Bot. 2017;104(8):1195–204.

    Article  PubMed  Google Scholar 

  35. Pagel E, Poschlod P, Reisch C. Habitat matters - Strong genetic and epigenetic differentiation in Linum catharticum from dry and wet grasslands. Ecol Evol. 2020;10(18):10271–80.

    Article  PubMed  PubMed Central  Google Scholar 

  36. Portis E, Acquadro A, Comino C, Lanteri S. Analysis of DNA methylation during germination of pepper (Capsicum annuum L.) seeds using methylation-sensitive amplification polymorphism (MSAP). Plant Sci. 2004;166(1):169–78.

  37. Fletcher H, Hickey I. BIOS Instant Notes in Genetics. 4th Ed. Boca Raton: CRC Press; 2013.

  38. Guarino F, Cicatelli A, Brundu G, Improta G, Triassi M, Castiglione S. The use of MSAP reveals epigenetic diversity of the invasive clonal populations of Arundo donax L. PLoS ONE. 2019;14(4):e0215096.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  39. Hand ML, Koltunow AMG. The genetic control of apomixis: asexual seed formation. Genetics. 2014;197:441–50.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  40. Mounger J, Ainouche ML, Bossdorf O, Cavé-Radet A, Li B, Parepa M, Salmon A, Yang J, Richards CL. Epigenetics and the success of invasive plants. Philos Trans R Soc Lond B Biol Sci. 1826;2021(376):20200117.

    Google Scholar 

  41. Verhoeven KJF, Preite V. Epigenetic variation in asexually reproducing organisms. Evolution. 2014;68(3):644–55.

    Article  PubMed  Google Scholar 

Download references


We are grateful to Xishuang Banna Tropical Botanical Garden of the CAS, South China Normal University and Prof. Guowen Xie for supporting the the collection of plant materials for this study.


This work was supported by the Natural Science Foundation of Guangdong Province [2018A030313478] and Science and Technology Program of GuangZhou [201707010257].

Author information

Authors and Affiliations



L.S., Y.X. and X.C. designed the study. Y.X., X.C.,Y.Y.,J.Z.,H.Y. collected the samples and performed the experiment. Y.X., X.C. performed data analysis. Y.X. drafted the manuscript. L.S. revised the manuscript and edited the language of the manuscript. All authors have read and agreed to the published version of the manuscript. The author(s) read and approved the final manuscript.

Corresponding author

Correspondence to Liying Song.

Ethics declarations

Ethics approval and consent to participate

We promise that Sphagneticola trilobata (L.) Pruski material from all regions has been licensed for collection. The authors confirm that all methods comply with relevant institutional, national, and international guidelines and legislations, including the collection of plant materials for this study. The plant materials were collected from the regions of Fuzhou city, Fujian; Haikou city, Hainan; Jinghong city, Yunan; Guangzhou city, Guangdong by the authors (Y.X.and X.C.), and identified as Sphagneticola trilobata (L.) Pruski by Prof. Guowen Xie of GZHU. The voucher specimen (No# GZHU-SLY-2020–07-FZ; No# GZHU-SLY-2020–07-HK; No# GZHU-SLY-2020–07-JH; No# GZHU-SLY-2020–07-GZ) were deposited at the Herbarium of School of Life Sciences, GZHU. L.S. undertook the formal identification of the plant material used in our study.

Consent for publication

Not applicable.

Competing interests

The authors declare that they have no competing interests.

Additional information

Publisher’s Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Supplementary Information

Additional file 1: Table S1.

Components of enzyme digestion system for MSAP (20 μl). Table S2. Formula of adapters for MSAP. Table S3. System of ligation for MSAP (20 μl). Table S4. Reaction system of pre-amplification (20 μl). Table S5. Reaction system of selection PCR amplification. Table S6. Sequences of adapters and primers for MSAP. Table S7. Reaction system of SSR amplification (10μl). Figure S1. Sampling sites distribution of Sphagneticola trilobata. Figure S2. Correlation between epigenetic (a) and genetic diversity (b) of Sphagneticola trilobata and environmental factors.

Rights and permissions

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit The Creative Commons Public Domain Dedication waiver ( applies to the data made available in this article, unless otherwise stated in a credit line to the data.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Xiao, Y., Chen, X., Yin, Y. et al. Comparative genetic and epigenetic of the Sphagneticola trilobata (L.) Pruski from different regions in China. BMC Plant Biol 23, 289 (2023).

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: