Skip to main content
Download PDF

Adaptive evolution and co-evolution of chloroplast genomes in Pteridaceae species occupying different habitats: overlapping residues are always highly mutated

BMC Plant Biology volume 23, Article number: 511 (2023) Cite this article

Abstract

Background

The evolution of protein residues depends on the mutation rates of their encoding nucleotides, but it may also be affected by co-evolution with other residues. Chloroplasts function as environmental sensors, transforming fluctuating environmental signals into different physiological responses. We reasoned that habitat diversity may affect their rate and mode of evolution, which might be evidenced in the chloroplast genome. The Pteridaceae family of ferns occupy an unusually broad range of ecological niches, which provides an ideal system for analysis.

Results

We conducted adaptive evolution and intra-molecular co-evolution analyses of Pteridaceae chloroplast DNAs (cpDNAs). The results indicate that the residues undergoing adaptive evolution and co-evolution were mostly independent, with only a few residues being simultaneously involved in both processes, and these overlapping residues tend to exhibit high mutations. Additionally, our data showed that Pteridaceae chloroplast genes are under purifying selection. Regardless of whether we grouped species by lineage (which corresponded with ecological niches), we determined that positively selected residues mainly target photosynthetic genes.

Conclusions

Our work provides evidence for the adaptive evolution of Pteridaceae cpDNAs, especially photosynthetic genes, to different habitats and sheds light on the adaptive evolution and co-evolution of proteins.

Peer Review reports

Background

The chloroplast genome of eukaryotes has been substantially reduced due to gene transfer to the host genome and simple gene loss [1]. Therefore, fewer than 5% of genes from the ancestral cyanobacteria remain in chloroplast DNA (cpDNA), mainly those encoding housekeeping and photosynthesis-associated proteins [2]. Nonetheless, chloroplasts are the sites of many important cellular processes, including photosynthesis and lipid, amino acid, chlorophyll, and carotenoid biosynthesis [3, 4]. Moreover, chloroplasts play important roles in plant adaptation to environmental stress, such as drought [5], salinity [6], extreme temperature [7], high light [8], and heavy metal stress [9]. Because chloroplasts act as the energy hubs of plant cells, their homeostasis is readily affected by environmental stress [3, 10]. Chloroplasts function as environmental sensors, transforming fluctuating environmental signals into physiological responses [11,12,13]. For instance, excessive reactive oxygen species (ROS) in plants can trigger retrograde signaling to subcellular compartments to readjust whole-cell metabolism in order to repair or turn over damaged macromolecules [7, 14]. In addition, the timely expression of chloroplast genes can promote adaptation to environmental fluctuations [15].

Genes encoding proteins with important functions are usually subject to strong evolutionary selection pressure. To measure selection pressure at the molecular level, the nucleotide substitution rate can be used to reflect the changes in protein-coding sequences: non-synonymous substitutions (dN) cause amino acid changes, while synonymous substitutions (dS) do not [16, 17]. dN/dS (ω) is a measure of natural selection that is widely used to detect genes related to environmental adaptation; a ω value greater than 1 is commonly believed to represent positive selection acting on a gene or gene lineage; a ω equal to 1 represents neutral selection, and a ω less than 1 represents purifying selection [18,19,20]. Many amino acids in functional proteins are conserved due to strong structural and functional constraints; thus, the average dN rate is rarely higher than the average dS rate [21, 22]. However, some genes may be positively selected under certain environmental conditions, resulting in higher ω values [23].

There may be an evolutionary association between amino acid residues of a given protein: co-evolution may occur between sequences within a protein that physically interact or are functionally relevant, so a change in residues at one site in the molecule may lead to a change in selection pressure at another site. However, amino acid residues located in crucial functional/structural regions will be subject to stronger selection restrictions because these sites may have a huge influence on the protein’s function [24]. Harmful mutations at such sites will be immediately eliminated from the population, but for a mutation that is actively selected for, compensation and substitution may occur to restore any fitness loss due to the mutation [25]. In other words, amino acid co-evolution or correlated mutation is a phenomenon whereby a deleterious substitution at one position is compensated for by another substitution elsewhere in the protein, so that the structure and function of the protein remain stable [26,27,28,29]. Notably, if the adaptive positively selected sites of a protein map to positions important for protein structure or function [30,31,32,33], and these sites show co-evolutionary relationships, this may indicate their functional/structural dependence [34]. Therefore, identifying the strategies employed by plants for adapting to different environments at the molecular level, along with their co-evolutionary dynamics, can shed light on the evolutionary pattern of species as well as the complex co-adaptation between residues in proteins.

Photosynthesis is particularly sensitive to environmental conditions and plays an important role in the unique niche occupied by land plants [35, 36]. Despite the wide variations in plant morphology and habitat throughout the plant kingdom, the structure and function of most cpDNA-encoded proteins are relatively well conserved. This conservation is largely believed to be driven by strong selection pressure related to the functional requirements of the proteins involved in photosynthesis [37]. Some functional genes with key roles in photosynthesis have undergone adaptive evolution with the “radiation”-type differentiation of species. For instance, genes encoding Rubisco components evolved under positive selection in most terrestrial plant lineages [38], and adaptive evolution of rbcL has also been observed in some aquatic plant lineages [39, 40], suggesting that RbcL subunits may have undergone continuous “fine-tuning” in different ecosystems. Similarly, atpB, psaB, and rbcL in Chlamydomonas sp. ICE-L show conserved adaptive evolution in extreme environments, which supports the notion that the adaptation of algae to extreme environments is related to the mode of selection of plastid proteins [41].

The Pteridaceae family of ferns contains five subfamilies, 53 genera, and an estimated 1,211 species, contributing to ~ 10% of extant leptosporangiate fern diversity [42, 43]. Pteridaceae have a diverse habitat distribution, ranging from wet tropical to arid regions, and they occupy an unusually broad range of ecological niches [44,45,46]. Most Pteridaceae species are terrestrial and inhabit open, often rocky environments, but representatives of some Adiantum and Pteris species are frequent in forests [47]. The Parkerioideae subfamily, represented by Ceratopteris, thrives in aquatic habitats, and Acrostichum is often associated with mangroves in the intertidal zone [48]. Species in the Cheilanthoideae subfamily are significant components of arid terrestrial habitats [49]. Vitarioid ferns, as sisters to the genus Adiantum L., are highly simplified and predominantly epiphytic species [50]. This provides an ideal system for studying the adaptive evolution and intra-molecular co-evolution of cpDNA in a wide range of niches. Previous studies on Pteridaceae have mainly focused on exploring their phylogenetic relationships based on their DNA sequences; however, few studies have focused on the molecular evolution of this family, especially between species inhabiting different environments. In this study, we divided 41 species of Pteridaceae into the terrestrial clade, the epiphytic clade, and the aquatic clade and analyzed the evolutionary dynamics of cpDNA of different habitat groups using common protein-coding sequence data sets to analyze (1) the molecular mechanism of adaptive evolution of Pteridaceae cpDNAs in different groups and (2) the intra-molecular co-evolution patterns of genes involved in this adaptive evolution.

Results

Construction of a phylogenetic tree for adaptive evolution and co-evolution analyses

We constructed a phylogenetic tree using protein sequences inferred from cpDNAs of 41 Pteridaceae species (Table S1). The topologies and support values from maximum likelihood (ML) and Bayesian inference (BI) analyses were similar, and the 41 species were resolved into five highly supported “subfamily clades”: Parkerioideae, Pteridoideae, Cryptogrammoideae, Cheilanthoideae, and Vittarioideae. Pteris and Adiantum were each resolved as monophyletic in this analysis, consistent with previous reports [43, 51,52,53]. We categorized this system into three clades based on their respective habitats: the aquatic clade belonging to the Parkerioideae subfamily, the epiphytic clade consisting of vittarioid ferns, while the terrestrial clade composing of the others (Fig. 1).

Fig. 1
figure 1

Phylogenetic frame of 41 Pteridaceae species constructed for adaptive evolution and co-evolution analysis. Out.: Outgroup; Park.: Parkerioideae; Pter.: Pteridoideae; Cryp.: Cryptogrammoideae; Chei.: Cheilanthoideae; Vitt.: Vittarioideae. + represents newly added cpDNA samples

Full size image

Analysis of selection pressure

We used the basic model to detect the selection pressure range in 41 Pteridaceae species. The ω values of Pteridaceae chloroplast protein-coding genes ranged from 0.0046 to 0.4432 (Table S2), and all 76 common protein-coding genes had ω less than 0.5, indicating that the protein-coding genes of Pteridaceae chloroplasts are dominated by purifying selection.

To test whether the evolutionary pattern of the chloroplast protein-coding genes in the Pteridaceae family is related to adaptability to different habitats, we set the terrestrial, aquatic, or epiphytic ancestral clade as the foreground clade successively when using the branch model and the other two non-target clades as background clade (Specific methods are detailed in Materials and Methods). The branch model will detect positive selection in a lineage only when the average dN of all sites is higher than the average dS [21]. We discarded five genes with extremely low dS values to avoid misestimating the ω value (999): accordingly, we eliminated infA from the terrestrial clade, psbJ from the epiphytic clade, and psaI, psaJ, and psbT from the aquatic clade. Using the remaining data sets, only when the terrestrial ancestral clade was used as the foreground clade did rbcL accept the M2 model, and the ω values became significant after p-value correction (ω = 0.0047, q-value = 0.0154). When the aquatic ancestral clade or the epiphytic ancestral clade was used as the foreground clade, we observed no significant positive selection after p-value correction. No matter which habitat was used as the foreground clade, the test results were mainly ω < 0.5 (Table S3), which also indicates that the Pteridaceae chloroplasts are dominated by purifying selection.

RELAX analysis

We identified a few genes with high but not statistically significant ω values in the above data. For instance, the ω values of psbL in the epiphytic ancestral clade and rps12 and rpl32 in the aquatic ancestral clade were greater than 1 and far higher than their respective background clades (Fig. 2). A higher ω value may be caused by the positive selection of specific genes or the relaxation of natural selection. Based on this result, we tested the selection intensity of Pteridaceae species from different habitats. We determined that the higher ω value of psbL in the epiphytic clade may be caused by the relaxation of selection pressure, while those of rps12 and rpl32 in the aquatic clade were caused by the intensification of selective pressure. In addition, rbcL with significant purifying selection in the terrestrial clade showed significant relaxation of selective pressure (Table S3).

Fig. 2
figure 2

The ω value of each retained protein-coding gene in different foreground clades under the branch model. Purple, blue, and green backgrounds represent genes related to the photosynthetic system, genetic system, and other functions respectively. * marked gene represents significance after p-value correction. Back.: Background clade; Ter.: Terrestrial clade; Aqu.: Aquatic clade; Epi.: Epiphytic clade

Full size image

Analysis of residues under positive selection

As most amino acids in a functional protein are under structural and functional constraints, and adaptive evolution likely affects only a few sites at a few time points [21], we used the site model, which ignores the variation in the ω value between branches, to determine whether any residues have undergone positive selection. Using the Bayes empirical Bayes (BEB) method in the site model with residue reference sites based on Pteris arisanensis, we detected positively selected residues within 20 chloroplast protein-coding genes that were statistically significant (Table 1).

Among these genes, atpE, matK, ndhF and rpoB rejected the M2a model in the M2a vs. M1a nested model, so we dismissed the results for these genes. Accordingly, 36 positively selected residues were retained in the M2a model (P > 95%), which are encoded within 16 genes. Of these genes, 10 are photosynthetic genes (atpA, atpB, atpH, ndhB, ndhK, ndhG, petB, psbF, psbL, and psbZ), four are genetic system genes (rpl16, rpoA, rpoC2, and rps7), and two are other functional genes (cemA and ycf3). Furthermore, atpA atpB and ndhK rejected the M8 model in the M8 vs. M8a nested model, so we dismissed the results for these genes. Accordingly, 59 positively selected residues were retained in the M8 model (P > 95%), which are encoded by 17 genes. Among these, nine are photosynthetic genes (atpE, atpH, ndhB, ndhF, ndhG, petB, psbF, psbL, and psbZ), five genetic system genes (rpl16, rpoA, rpoB, rpoC2, and rps7), and three are other functional genes (cemA, matK, and ycf3). Together, the results of the M2a and M8 models show that most of the positively selected residues are encoded in photosynthetic genes of Pteridaceae.

The branch-site model provided evidence of positive selection on different habitat groups of Pteridaceae. Using the BEB method, we did not detect positively selected residues in the terrestrial ancestral clade. Instead, we identified 3 positively selected residues in the epiphytic ancestral clade and 2 in the aquatic ancestral clade (P > 95%, with residue reference sites based on P. arisanensis) (Table 1). These residues are encoded by ndhI (S50), psaB (T511), and rps19 (R17) in the epiphytic ancestral clade. In the aquatic ancestral clade, they are encoded by ndhG (W121) and rpl33 (T16). These residues showed strong positive selection in the branch-site model, suggesting they may be critical for adaptive evolution in the epiphytic and aquatic ancestral clades of Pteridaceae.

Table 1 Parameter estimates for different model selection tests
Full size table

Intra-molecular co-evolution analysis

The evolution of residues depends on their intrinsic nucleotide mutation rates and the effects imposed by their complex co-evolutionary networks [24]. Co-evolution dynamics can highlight the intricate co-adaptive relationships between residues in a protein under an estimated timescale [34]. Using the protein sequences of P. arisanensis as a reference, we conducted co-evolution analysis of the genes encoding positively selected residues (i.e., adaptively evolving sites) identified using the site model (17 of the 20 protein sequences tested, except for PsbF, PsbL, and PsbZ). The 17 remaining protein sequences (AtpA, AtpB, AtpE, AtpH, CemA, MatK, NdhB, NdhF, NdhG, NdhK, PetB, Rpl16, RpoA, RpoB, RpoC2, Rps7, and Ycf3) contained co-evolved residue pairs (Table 2 displays proteins containing overlapping residues in the co-evolution and adaptive evolution analyses; other co-evolved residue pairs are shown in Table S4). Among these, RpoC2 had the most co-evolved residue pairs (363 pairs), followed by NdhB (90 pairs), AtpB (74), RpoB (50), NdhF (45), CemA (43), AtpA (38), RpoA (37), MatK (33), PetB (20), NdhG (6), NdhK (5), Rpl16 (5), AtpE (4), Ycf3 (4), AtpH (1), and Rps7 (1).

In complex co-evolution networks, some residues frequently appear in different co-evolved pairs. In particular, evolutionary dependence is identified between sites belonging to different domains, which plays an important role in the formation of residue co-evolution networks. For instance, based on the co-evolved pairs (S125 & S252 and S252 & P461) in the AtpB protein sequence, although S125 and P461 do not show direct co-evolution, there may be some indirect co-evolution pressure.

Table 2 Information about intra-molecular co-evolved residue pairs for protein sequences containing overlapping residues in co-evolution and adaptive evolution analyses
Full size table

Relationship between adaptive evolution and co-evolution

Having more positively selected residues makes it easier to observe the relationship between positively selected residues and co-evolved residues. Therefore, we used the union of the M2a and M8 model results in the following analyses. We analyzed the relationship between positively selected residues and co-evolved residues by predicting the tertiary structures of protein sequences. The co-evolved residues were scattered throughout the protein sequences (Fig. 3), and some sites were not in a direct co-evolutionary relationship even if they had a relatively close linear distance in the primary structure. For example, in the AtpA protein sequence, the co-evolved residues L82 and T83 were both predicted to interact with the distal sites in the tertiary structure (L82 & M220, distance: 35.4 Å; T83 & K255, distance: 43.8 Å) (Table 2). In addition, the positively selected residues and co-evolved residues within the protein sequence were mostly independent of each other, with only a few overlapping sites: AtpA (L308), AtpB (H123, S252, and S474), AtpH (L80), MatK (E118), NdhF (L629), PetB (S2 and R118), and RpoB (R575) (Fig. 3). After visualizing the co-evolved residues using Jalview [54], we determined that the co-evolved residues that overlapped with positively selected residues tended to be highly mutated, while the non-overlapping co-evolved residues had fewer mutations (Fig. 4). In this study, we defined a residue site as a high mutation site when the same residue was present in less than half (The results of multiple sequence alignment for specific co-evolved and positively selected residues can be found in Figure S1 and Figure S2, respectively). We employed a chi-square test to assess the conservativeness of residues in overlapping positions. The results indicated a statistical preference for overlapping residues to be located in positions with a low conservation ( χ2: 52.30, p-value: 0).

Fig. 3
figure 3

Distribution of positively selected residues and co-evolved residues in the AtpA, AtpB, AtpH, MatK, NdhF, PetB, and RpoB proteins of Pteridaceae. All protein tertiary structures were predicted by homology based on P. arisanensis. Purple represents co-evolution residues, magenta represents adaptive evolution residues, and firebrick red represents both and is labeled accordingly

Full size image
Fig. 4
figure 4

The proportion of differential residues in multiple sequence alignment of co-evolved sites, positively selected sites, and overlapping sites between the two

Full size image

Discussion

Adaptive evolutionary responses to changing environmental conditions result in accelerated evolution and the functional evolution of specific stress-response proteins that favor improved fitness in the new environment [55,56,57]. In this study, we observed that all protein-coding genes in the Pteridaceae chloroplast genome had ω < 0.5, indicating that Pteridaceae cpDNA has mainly undergone purifying selection. However, only rbcL showed a ω value (0.0047) in the terrestrial ancestral clade significantly different from that in the background clade, indicating that this gene has experienced significant purifying selection. Rubisco is crucial for photosynthesis, as it not only fixes CO2 [58, 59] but also improves plant growth performance under stress [60]. Truncations and mutations of the conserved N terminus of RbcL can substantially affect Rubisco activity [61]. RbcL has evolved adaptively in different environments and is thought to have undergone “continuous fine-tuning” in different ecosystems [38,39,40]. However, in the terrestrial clade of the Pteridaceae, rbcL is under significant purifying selection and is in a state of relaxed selection. On the other hand, the branch-site model did not detect any positively selected residues in the terrestrial ancestral clade. Therefore, we speculate that the ancestors of Pteridaceae likely grew on land and had completed adaptive evolution, so the terrestrial clade of Pteridaceae is now under purifying selection. Under purifying selection, most non-synonymous substitutions will be purified, so rbcL exhibits significant purifying selection.

When we applied the likelihood ratio test (LRT) and BEB test without distinguishing between lineages, 20 genes showed residues under positive selection (significant at the 95% level) (Table 1). These residues were mainly encoded within photosynthetic genes, including genes encoding subunits of ATP synthase (atpA, atpB, atpE, and atpH), the predominant site of photosynthetic flux control [62]; NADH dehydrogenase (ndhB, ndhF, ndhG, and ndhK), which is important for adaptation of the photosynthetic mechanism to abiotic stress [63]; the cytochrome b6f complex (petB); and photosystem II (psbF, psbL, and psbZ). These genes encode indispensable components of photosynthesis, and the finding that they are under positive selection indicates that photosynthetic genes play key roles in the adaptive evolution of Pteridaceae. The results of the site model also revealed positively selected genes related to plastid-encoded RNA polymerase (rpoA, rpoB, and rpoC2), which are involved in the transcription of photosynthesis-related genes [64, 65]; and small and large ribosome subunits (rpl16 and rps7). The other functional genes showing positive selection were cemA, matK, and ycf3.

To further elucidate the evolutionary pattern of Pteridaceae chloroplast-encoded proteins, we conducted a co-evolution analysis of the protein sequences encoded by 20 genes including positively selected residues, finding that 17 showed intra-molecular co-evolution. Among these, AtpB, NdhB, RpoB, and RpoC2 encoded abundant co-evolved residues (over 50 co-evolved pairs) and formed a complex co-evolution network (Figure S3); this complexity may be due to their large size, as larger proteins have more residue interactions and more possibilities for co-evolution [25]. The co-evolution approach emphasizes the role of spatial structure information in the study of protein evolution and provides strong evidence to reveal the structural and functional correlation of residue sites [66]. Co-evolved residue pairs often correspond to spatially proximal residues in protein structures [67]. In this study, in addition to the co-evolved pairs with relatively close linear and structural distances, we also observed many co-evolutionary relationships between distal residues at the tertiary structure level, which far outnumbered those of proximal residues (Table 2), suggesting that these co-evolved pairs may be more prone to functional co-evolution [68].

When a species is faced with an extreme environment, some residues within proteins undergo adaptive changes [30,31,32,33]. To maintain protein stability, this effect may be mitigated by compensatory changes in other amino acid residues [25]. Therefore, the adaptive evolution of residues in a protein and the intra-molecular co-evolution will tend to be linked. However, our results were contrary to this expectation: residues undergoing adaptive evolution and co-evolution are scattered in the three-dimensional structure of proteins, and their overlap is limited (overlapping residues/adaptively evolving residues: 10/64; overlapping residues/co-evolved residues: 10/741), indicating that the positions of these two categories of sites in the protein sequence are nearly independent. After visualizing multiple amino acid sequences, whereas adaptively evolving residues generally exhibit lower conservation (Figure S2), co-evolved residues tend to display relatively higher conservation (Figure S1). Interestingly, the co-evolved residues that overlap with adaptively evolving residues generally show low conservation: that is, the sequences of the overlapping sites were always highly mutated. On the one hand, non-synonymous substitutions are more strongly affected by natural selection than are synonymous substitutions [69]. If adaptive evolution leads to an increase in the non-synonymous substitution rates, it may indicate that the species is actively adapting to the new environment. Therefore, adaptive evolution is expected to increase the rate of non-synonymous substitutions. Notably, functionally/structurally important residues are often subject to strong selection constraints and thus tend to be conserved, even these are subject to mutation under extreme environmental stress. On the other hand, because the structure and function of a macromolecule depend on complex interactions between residues, such changes can alter crucial interactions with other residues. Thus, the fitness effect of a mutation at a given location may depend on the state of interacting residues, resulting in non-independent evolution (co-evolution) [70]. This co-evolution should generally be quite slow, with these residues showing greater conservation [25]. Moreover, the geometry of interactions in the tertiary structures of proteins may vary in time and sequence, with the residues involved potentially changing over time during the evolutionary history of a molecule [25]. Thus, co-evolved residue pairs and adaptively evolved residue sites may vary in time and space. Therefore, one possible explanation for the presence of these overlapping residues is that adaptive evolution and intra-molecular co-evolution in species may be driven by distinct mechanisms, and they happen to overlap. Notably, these overlapping residues tend to exhibit lower conservation.

We also detected genes encoding positively selected residues in the epiphytic and the aquatic clades (95% level, BEB method): ndhG (W121) and rpl33 (T16) in the aquatic ancestral clade and ndhI (S50), psaB (T511), and rps19 (R17) in the epiphytic ancestral clade. The habitats of epiphytic species in a tree canopy are diverse, with some limited to smaller branches and others to larger branches or trunks [71]. The aquatic A. speciosum usually grows in the shady understory of a mangrove forest and is frequently flooded by tides [72], while Ceratopteris species are often restricted to aquatic habitats, such as those near ponds and streams [73]. Such shade-living species are especially vulnerable to damage to the photosynthetic machinery due to excessive light. The genes that were positively selected in these clades encode proteins crucial to this machinery: NdhG, NdhI, and PsaB are involved in photosystem I (PSI) cyclic electron transport [74,75,76], and the NDH complex appears to be particularly important for enabling the photosynthetic machinery to adapt to stress conditions [63] and/or alleviating stress [75]. Rpl33 and Rps19 are ribosomal proteins, with Rpl33 being necessary to maintain sufficient cpDNA translation under cold conditions [77]. However, we did not detect positively selected residues in the terrestrial ancestral clade, possibly because the Pteridaceae are descended from a terrestrial ancestor that was previously adapted to such environments. In summary, some residues encoded within Pteridaceae chloroplast genes showed positive selection, indicating that they play an important role in the adaptive evolution to different habitats.

Conclusion

Pteridaceae cpDNAs have mainly undergone purifying selection. Only rbcL in the terrestrial clade showed significant purification selection and is in a state of relaxed selection. Regardless of whether the lineage was distinguished, the positively selected residues were mostly encoded by photosynthetic genes, indicating that photosynthetic genes play an important role in the adaptive evolution of Pteridaceae species. We carried out a co-evolution analysis of 20 genes encoding adaptively evolved residues to explore the complex evolutionary pattern of proteins. These adaptively evolved sites and co-evolved sites were mostly independent, with only a few overlapping sites, and the amino acid sequences of these overlapping sites were always highly mutated. These overlapping sites may be due to the different mechanisms between adaptive evolution and co-evolution, and they may overlap in different spaces and times. Obtaining more structural/functional information for these overlapping sites will be crucial for a deeper understanding of the relationship between adaptive evolution and co-evolution. Here, we present evidence at the molecular level about the adaptive evolution of Pteridaceae cpDNAs to different habitats and provide insight into the adaptive evolution and co-evolution of proteins.

Materials and methods

Sample collection

Fresh leaves of Pteris ensiformis, P. arisanensis, and Taenitis blechnoides were sampled from the campus of Shenzhen Fairy Lake Botanical Garden (location: E114°09′, N22°34′; altitude: 944 m), quickly frozen in liquid nitrogen, and stored at − 80 °C until use. The plant materials used in the study were identified by Ting Wang, and specimens were stored in the Herbarium of the College of Life Sciences, South China Agricultural University (Specimen numbers of P. ensiformis, P. arisanensis, and T. blechnoides are GXL20210903, GXL20210904 and GXL20210905, respectively).

Library preparation, sequencing and genome assembly

DNA was extracted from the samples using a Tiangen Plant Genome DNA Kit (Tiangen Biotech Co., Ltd., Beijing, China) according to the manufacturer’s instructions. The Illumina NovaSeq6000 platform was used for sequencing. The complete chloroplast genomes of P. arisanensis and T. blechnoides were assembled using GetOrganelle [78]. However, the complete chloroplast genome of P. ensiformis failed to assemble into a circle using GetOrganelle and was instead assembled using Novoplasty [79]. The sequences were submitted to the National Center for Biotechnology Information (NCBI) under GenBank accession numbers OP441371 (P. arisanensis), OP743918 (P. ensiformis), and OP743919 (T. blechnoides).

DNA sequence alignment and phylogenetic analysis

The complete chloroplast genomes of 38 Pteridaceae species and Alsophila denticulata (outgroup) were downloaded from GenBank (Table S1). Combining our three newly sequenced species, a total of 41 Pteridaceae species were examined, covering all subfamilies. Seventy-six common but different protein-coding sequences of these species were retained (Table S1), MAFFT [80] was used to perform sequence alignment, and the gap area was deleted to exclude poorly aligned positions. PhyloSuite [81] was used to concatenate these sequences into a dataset for phylogenetic analysis. The maximum likelihood (ML) tree was inferred using RAxML [82], GTRGAMMAI was selected as the nucleotide substitution model, and bootstrap values for each branch were obtained by performing 1,000 bootstrap replicates. The Bayesian inference (BI) tree was established by MrBayes [83] and was estimated by running 2,000,000 generations (Nst = 6, rates = invgamma).

Analysis of selection pressure and adaptive evolution

The 76 common non-repeating protein-coding sequences were used to build independent data sets (Table S1). The 41 Pteridaceae species were classified into three clades based on their habitat: the terrestrial, aquatic, and epiphyte clades. Before performing sequence alignment, the stop codon caused by RNA editing inside the sequences was modified, and the tail stop codon was removed. Alignment gaps and uncertainties were deleted to avoid false positives.

Selection pressure was analyzed using the CODEML [16] program. The ω values of each common protein-coding gene were calculated under the basic model (model = 0, Nsites = 0, which assumes no site-wise or branch-wise dN/dS variation). The branch model was applied by comparing the single-ratio model (M0) and two-ratio model (M2) (using the likelihood ratio test (LRT) with a χ2 distribution); the M0 hypothesis was rejected if P < 0.05. The M0 model assumes that all Pteridaceae clades have the same ratio of non-synonymous to synonymous substitution rates; the M2 model assumes different ω between the foreground clade and background clade. The false discovery rate (FDR) correction was applied to the P values calculated above [84].

Positive selection models (M2a and M8) and null hypothesis models (M1a, M7, and M8a) provided by PAML [16] were used to perform site adaptive evolutionary analyses on a dataset of shared genes from 41 Pteridaceae species. Three sets of nested models (M2a vs. M1a, M8 vs. M7, and M8a vs. M8) were used to infer the genes experiencing positive selection. If the positive selection model significantly outperformed the null hypothesis model (p < 0.05), the genes were assumed to be under positive selection, where M8a vs. M8 resulted in lower false positive results.

The branch-site model allows different site-encoding genes to have different values of ω in different branches of the phylogenetic tree in order to test whether positive selection acts on certain sites in the foreground clade. Model A (model = 2, Nsites = 2, fixed omega = 0, omega = 2) assumes that only the foreground clade undergoes positive selection; Model A null (model = 2, Nsites = 2, fixed omega = 1, omega = 1) fixes the ω of the foreground clade in model A to 1. If Model A is significantly better than Model A null (p < 0.05), then the gene has undergone positive selection in the foreground clade.

In the above models, the LRT was compared to a χ2 null distribution with the corresponding degrees of freedom. All sites under positive selection were retrieved using the Bayes empirical Bayes (BEB) method.

Natural selection pressure analyses

Parameter K was calculated for the data by running RELAX (implemented in HYPHY v2.5.42) [85] to test the relaxation of natural selection, the selection pressure in this context includes both purifying selection and positive selection. The K parameter is related to the value of ω as (ω background)K = (ω foreground) [85]. K < 1 is indicative of relaxed natural selection and K > 1 suggests intensification in the test compared to the background. RELAX was used for LRT analysis by comparing the model with K = 1 to the model with K < 1 (or K > 1). If K < 1, P < 0.05 indicates that the test clade was under significantly relaxed selection, and if K > 1, P < 0.05 indicates that the test clade was under significantly intensified selection.

Co-evolution analysis

To further explore adaptive evolution at the molecular level, co-evolution analysis was conducted on the protein sequences encoded by genes that have undergone adaptive evolution in a specific lineage, and the structures of proteins with co-evolved residue pairs were predicted in order to understand their co-evolutionary mechanism. Co-evolution analysis was performed using the program CAPS [24], which reveals structural and functional correlations between sites by detecting whether amino acid site variants are associated. The alpha-value in the program was set to 0.01; the random sampling value was set to 1000. Protein 3D structure was predicted using Phyre2 [86] online in INTENSIVE mode and visualized using PyMOL (http://www.pymol.org). Multiple sequences of co-evolved residues were checked using Jalview [54].

Data availability

These sequence data have been submitted to the GenBank databases under accession numbers OP441371, OP743918 and OP743919. The data underlying this article are available in the article and in its online supplementary material.

References

  1. Timmis JN, Ayliffe MA, Huang CY, Martin W. Endosymbiotic gene transfer: organelle genomes forge eukaryotic chromosomes. NAT REV GENET. 2004;5(2):123–35.

    CAS  PubMed  Google Scholar 

  2. Li M, Kim C. Chloroplast ROS and stress signaling. Plant Commun. 2022;3(1):100264.

    CAS  PubMed  Google Scholar 

  3. Chan KX, Crisp PA, Estavillo GM, Pogson BJ. Chloroplast-to-nucleus communication: current knowledge, experimental strategies and relationship to drought stress signaling. Plant Signal Behav. 2010;5(12):1575–82.

    CAS  PubMed Central  PubMed  Google Scholar 

  4. Gu X, Zhu M, Su Y, Wang T. A large intergenic spacer leads to the increase in genome size and sequential gene movement around IR/SC boundaries in the chloroplast genome of Adiantum malesianum (Pteridaceae). INT J MOL SCI. 2022;23(24):15616.

    CAS  PubMed Central  PubMed  Google Scholar 

  5. Zhao C, Haigh AM, Holford P, Chen ZH. Roles of chloroplast retrograde signals and ion transport in plant drought tolerance. INT J MOL SCI. 2018;19(4):963.

    PubMed Central  PubMed  Google Scholar 

  6. Krasensky J, Broyart C, Rabanal FA, Jonak C. The redox-sensitive chloroplast trehalose-6-phosphate phosphatase AtTPPD regulates salt stress tolerance. Antioxid Redox Signal. 2014;21(9):1289–304.

    CAS  PubMed Central  PubMed  Google Scholar 

  7. Luo S, Kim C. Current understanding of temperature stress-responsive chloroplast FtsH metalloproteases. INT J MOL SCI. 2021;22(22):12106.

    CAS  PubMed Central  PubMed  Google Scholar 

  8. Phee BK, Cho JH, Park S, Jung JH, Lee YH, Jeon JS, Bhoo SH, Hahn TR. Proteomic analysis of the response of Arabidopsis chloroplast proteins to high light stress. Proteomics. 2004;4(11):3560–8.

    CAS  PubMed  Google Scholar 

  9. Khan M, Nawaz N, Ali I, Azam M, Rizwan M, Ahmad P, Ali S. Regulation of photosynthesis under metal stress. Photosynthesis Productivity and Environmental Stress 2019:95–105.

  10. Zhu JK. Abiotic stress signaling and responses in plants. Cell. 2016;167(2):313–24.

    CAS  PubMed Central  PubMed  Google Scholar 

  11. Chan KX, Phua SY, Crisp P, McQuinn R, Pogson BJ. Learning the languages of the chloroplast: retrograde signaling and beyond. ANNU REV PLANT BIOL. 2016;67:25–53.

    CAS  PubMed  Google Scholar 

  12. Woodson JD. Chloroplast quality control - balancing energy production and stress. NEW PHYTOL. 2016;212(1):36–41.

    CAS  PubMed  Google Scholar 

  13. Zhao C, Wang Y, Chan KX, Marchant DB, Franks PJ, Randall D, Tee EE, Chen G, Ramesh S, Phua SY, et al. Evolution of chloroplast retrograde signaling facilitates green plant adaptation to land. Proc Natl Acad Sci. 2019;116(11):5015–20.

    CAS  PubMed Central  PubMed  Google Scholar 

  14. Mamaeva A, Taliansky M, Filippova A, Love AJ, Golub N, Fesenko I. The role of chloroplast protein remodeling in stress responses and shaping of the plant peptidome. NEW PHYTOL. 2020;227(5):1326–34.

    CAS  PubMed  Google Scholar 

  15. Zhang Y, Zhang A, Li X, Lu C. The role of chloroplast gene expression in plant responses to environmental stress. INT J MOL SCI. 2020;21(17):6082.

    CAS  PubMed Central  PubMed  Google Scholar 

  16. Yang Z. PAML 4: phylogenetic analysis by maximum likelihood. MOL BIOL EVOL. 2007;24(8):1586–91.

    CAS  PubMed  Google Scholar 

  17. Garud NR, Pollard KS. Population genetics in the human microbiome. TRENDS GENET. 2020;36(1):53–67.

    CAS  PubMed  Google Scholar 

  18. Scobeyeva VA, Artyushin IV, Krinitsina AA, Nikitin PA, Antipin MI, Kuptsov SV, Belenikin MS, Omelchenko DO, Logacheva MD, Konorov EA, et al. Gene loss, pseudogenization in plastomes of genus Allium (Amaryllidaceae), and putative selection for adaptation to environmental conditions. FRONT GENET. 2021;12:674783.

    CAS  PubMed Central  PubMed  Google Scholar 

  19. Zhang L, Yan HF, Wu W, Yu H, Ge XJ. Comparative transcriptome analysis and marker development of two closely related primrose species (Primula Poissonii and Primula wilsonii). BMC Genomics. 2013;14:329.

    CAS  PubMed Central  PubMed  Google Scholar 

  20. Wang B, Jiang B, Zhou Y, Su Y, Wang T. Higher substitution rates and lower dN/dS for the plastid genes in Gnetales than other gymnosperms. BIOCHEM SYST ECOL. 2015;59:278–87.

    CAS  Google Scholar 

  21. Yang Z, Nielsen R. Codon-substitution models for detecting molecular adaptation at individual sites along specific lineages. MOL BIOL EVOL. 2002;19(6):908–17.

    CAS  PubMed  Google Scholar 

  22. Smith DR. Mutation rates in plastid genomes: they are lower than you might think. GENOME BIOL EVOL. 2015;7(5):1227–34.

    CAS  PubMed Central  PubMed  Google Scholar 

  23. Zhang ZR, Yang X, Li WY, Peng YQ, Gao J. Comparative chloroplast genome analysis of Ficus (Moraceae): insight into adaptive evolution and mutational hotspot regions. FRONT PLANT SCI. 2022;13:965335.

    PubMed Central  PubMed  Google Scholar 

  24. Fares MA, McNally D. CAPS: coevolution analysis using protein sequences. Bioinformatics. 2006;22(22):2821–2.

    CAS  PubMed  Google Scholar 

  25. Chaurasia S, Dutheil JY. The structural determinants of intra-protein compensatory substitutions. MOL BIOL EVOL 2022, 39(4).

  26. Jacob E, Unger R, Horovitz A. Codon-level information improves predictions of inter-residue contacts in proteins by correlated mutation analysis. ELIFE. 2015;4:e8932.

    Google Scholar 

  27. Basu S, Bahadur RP. Conservation and coevolution determine evolvability of different classes of disordered residues in human intrinsically disordered proteins. Proteins. 2022;90(3):632–44.

    CAS  PubMed  Google Scholar 

  28. Dos SR, Jiang X, Martinez L, Morcos F. Coevolutionary signals and structure-based models for the prediction of protein native conformations. Methods Mol Biol. 2019;1851:83–103.

    Google Scholar 

  29. Codoner FM, Fares MA, Elena SF. Adaptive covariation between the coat and movement proteins of prunus necrotic ringspot virus. J VIROL. 2006;80(12):5833–40.

    CAS  PubMed Central  PubMed  Google Scholar 

  30. Song W, Qin Y, Zhu Y, Yin G, Wu N, Li Y, Hu Y. Delineation of plant caleosin residues critical for functional divergence, positive selection and coevolution. BMC EVOL BIOL. 2014;14:124.

    PubMed Central  PubMed  Google Scholar 

  31. Lynch VJ. Inventing an arsenal: adaptive evolution and neofunctionalization of snake venom phospholipase A2 genes. BMC EVOL BIOL. 2007;7:2.

    PubMed Central  PubMed  Google Scholar 

  32. Wang M, Zhang X, Zhao H, Wang Q, Pan Y. FoxO gene family evolution in vertebrates. BMC EVOL BIOL. 2009;9:222.

    PubMed Central  PubMed  Google Scholar 

  33. Zhao F, Qin S. Evolutionary analysis of phycobiliproteins: implications for their structural and functional relationships. J MOL EVOL. 2006;63(3):330–40.

    CAS  PubMed  Google Scholar 

  34. Sen L, Fares MA, Liang B, Gao L, Wang B, Wang T, Su YJ. Molecular evolution of rbcL in three gymnosperm families: identifying adaptive and coevolutionary patterns. BIOL DIRECT. 2011;6:29.

    PubMed Central  PubMed  Google Scholar 

  35. Kapoor B, Kumar P, Verma V, Irfan M, Sharma R, Bhargava B. How plants conquered land: evolution of terrestrial adaptation. J Evol Biol 2022.

  36. Wellman CH, Strother PK. The terrestrial biota prior to the origin of land plants (embryophytes): a review of the evidence. Palaeontology. 2015;58(4):601–27.

    Google Scholar 

  37. Bungard RA. Photosynthetic evolution in parasitic plants: insight from the chloroplast genome. BioEssays. 2004;26(3):235–47.

    CAS  PubMed  Google Scholar 

  38. Kapralov MV, Filatov DA. Widespread positive selection in the photosynthetic rubisco enzyme. BMC EVOL BIOL. 2007;7:73.

    PubMed Central  PubMed  Google Scholar 

  39. Iida S, Miyagi A, Aoki S, Ito M, Kadono Y, Kosuge K. Molecular adaptation of rbcL in the heterophyllous aquatic plant Potamogeton. PLoS ONE. 2009;4(2):e4633.

    PubMed Central  PubMed  Google Scholar 

  40. Miwa H, Odrzykoski IJ, Matsui A, Hasegawa M, Akiyama H, Jia Y, Sabirov R, Takahashi H, Boufford DE, Murakami N. Adaptive evolution of rbcL in Conocephalum (Hepaticae, bryophytes). Gene. 2009;441(1–2):169–75.

    CAS  PubMed  Google Scholar 

  41. Zhang Z, An M, Miao J, Gu Z, Liu C, Zhong B. The Antarctic Sea ice alga Chlamydomonas sp. ICE-L provides insights into adaptive patterns of chloroplast evolution. BMC PLANT BIOL. 2018;18(1):53.

    PubMed Central  PubMed  Google Scholar 

  42. Schuettpelz E, Schneider H, Smith AR, Hovenkamp P, Prado J, Rouhan G, Salino A, Sundue M, Almeida TIE, Parris B, et al. A community-derived classification for extant lycophytes and ferns. J SYST EVOL. 2016;54(6):563–603.

    Google Scholar 

  43. Zhang L, Zhou XM, Lu NT, Zhang LB. Phylogeny of the fern subfamily Pteridoideae (Pteridaceae; Pteridophyta), with the description of a new genus: Gastoniella. MOL PHYLOGENET EVOL. 2017;109:59–72.

    PubMed  Google Scholar 

  44. Schuettpelz E, Schneider H, Huiet L, Windham MD, Pryer KM. A molecular phylogeny of the fern family Pteridaceae: assessing overall relationships and the affinities of previously unsampled genera. MOL PHYLOGENET EVOL. 2007;44(3):1172–85.

    CAS  PubMed  Google Scholar 

  45. Rothfels CJ, Windham MD, Grusz AL, Gastony GJ, Pryer KM. Toward a monophyletic Notholaena (Pteridaceae): resolving patterns of evolutionary convergence in xeric-adapted ferns. Taxon. 2008;57(3):712–24.

    Google Scholar 

  46. Abotsi KE, Radji AR, Rouhan G, Dubuisson JY, Kokou K. The Pteridaceae family diversity in Togo. Biodivers Data J. 2015;3:e5078.

    Google Scholar 

  47. Prado J, Del Nero Rodrigues C, Salatino A, Salatino MLF. Phylogenetic relationships among Pteridaceae, including Brazilian species, inferred from rbcL sequences. Taxon. 2007;56(2):355–68.

    Google Scholar 

  48. Wei X, He M, Jin X. The first complete chloroplast genome of a mangrove fern, Acrostichum speciosum. Mitochondrial DNA. Part B. Resources 2020, 5(2):1413–1414.

  49. Kirkpatrick RE. Investigating the monophyly of Pellaea (Pteridaceae) in the context of a phylogenetic analysis of cheilanthoid ferns. SYST BOT. 2007;32(3):504–18.

    Google Scholar 

  50. Schuettpelz E, Chen C, Kessler M, Pinson JB, Johnson G, Davila A, Cochran AT, Huiet L, Pryer KM. A revised generic classification of vittarioid ferns (Pteridaceae) based on molecular, micromorphological, and geographic data. Taxon. 2016;65(4):708–22.

    Google Scholar 

  51. Zhang L, Zhang LB. Phylogeny and systematics of the brake fern genus Pteris (Pteridaceae) based on molecular (plastid and nuclear) and morphological evidence. MOL PHYLOGENET EVOL. 2018;118:265–85.

    PubMed  Google Scholar 

  52. Schuettpelz E, Davila A, Prado J, Hirai RY, Yatskievych G. Molecular phylogenetic and morphological affinities of Adiantum senae (Pteridaceae). TAXON 2014, 63(2):258–264.

  53. Lu JM, Wen J, Lutz S, Wang YP, Li DZ. Phylogenetic relationships of Chinese Adiantum based on five plastid markers. J PLANT RES. 2012;125(2):237–49.

    PubMed  Google Scholar 

  54. Waterhouse AM, Procter JB, Martin DM, Clamp M, Barton GJ. Jalview Version 2 - a multiple sequence alignment editor and analysis workbench. Bioinformatics. 2009;25(9):1189–91.

    CAS  PubMed Central  PubMed  Google Scholar 

  55. Nevo E. Evolution of genome-phenome diversity under environmental stress. Proc Natl Acad Sci U S A. 2001;98(11):6233–40.

    CAS  PubMed Central  PubMed  Google Scholar 

  56. Wright BE. Does selective gene activation direct evolution? FEBS LETT. 1997;402(1):4–8.

    CAS  PubMed  Google Scholar 

  57. Yang J, Wang ZL, Zhao XQ, Wang DP, Qi DL, Xu BH, Ren YH, Tian HF. Natural selection and adaptive evolution of leptin in the Ochotona family driven by the cold environmental stress. PLoS ONE. 2008;3(1):e1472.

    PubMed Central  PubMed  Google Scholar 

  58. Iniguez C, Aguilo-Nicolau P, Galmes J. Improving photosynthesis through the enhancement of Rubisco carboxylation capacity. Biochem Soc Trans. 2021;49(5):2007–19.

    CAS  PubMed  Google Scholar 

  59. Galmes J, Capo-Bauca S, Niinemets U, Iniguez C. Potential improvement of photosynthetic CO2 assimilation in crops by exploiting the natural variation in the temperature response of Rubisco catalytic traits. CURR OPIN PLANT BIOL. 2019;49:60–7.

    CAS  PubMed  Google Scholar 

  60. Qu Y, Sakoda K, Fukayama H, Kondo E, Suzuki Y, Makino A, Terashima I, Yamori W. Overexpression of both Rubisco and Rubisco activase rescues rice photosynthesis and biomass under heat stress. PLANT CELL ENVIRON. 2021;44(7):2308–20.

    CAS  PubMed  Google Scholar 

  61. Ng J, Guo Z, Mueller-Cajar O. Rubisco activase requires residues in the large subunit N terminus to remodel inhibited plant Rubisco. J BIOL CHEM. 2020;295(48):16427–35.

    CAS  PubMed Central  PubMed  Google Scholar 

  62. Neupane P, Bhuju S, Thapa N, Bhattarai HK. ATP synthase: structure, function and inhibition. Biomol Concepts. 2019;10(1):1–10.

    CAS  PubMed  Google Scholar 

  63. Rumeau D, Peltier G, Cournac L. Chlororespiration and cyclic electron flow around PSI during photosynthesis and plant stress response. PLANT CELL ENVIRON. 2007;30(9):1041–51.

    CAS  PubMed  Google Scholar 

  64. Krause K, Maier RM, Kofer W, Krupinska K, Herrmann RG. Disruption of plastid-encoded RNA polymerase genes in Tobacco: expression of only a distinct set of genes is not based on selective transcription of the plastid chromosome. Mol Gen Genet. 2000;263(6):1022–30.

    CAS  PubMed  Google Scholar 

  65. Hajdukiewicz PT, Allison LA, Maliga P. The two RNA polymerases encoded by the nuclear and the plastid compartments transcribe distinct groups of genes in Tobacco plastids. EMBO J. 1997;16(13):4041–8.

    CAS  PubMed Central  PubMed  Google Scholar 

  66. Oteri F, Sarti E, Nadalin F, Carbone A. iBIS2Analyzer: a web server for a phylogeny-driven coevolution analysis of protein families. NUCLEIC ACIDS RES 2022.

  67. de Juan D, Pazos F, Valencia A. Emerging methods in protein co-evolution. NAT REV GENET. 2013;14(4):249–61.

    PubMed  Google Scholar 

  68. Fares MA, Travers SA. A novel method for detecting intramolecular coevolution: adding a further dimension to selective constraints analyses. Genetics. 2006;173(1):9–23.

    CAS  PubMed Central  PubMed  Google Scholar 

  69. Shirai K, Inomata N, Mizoiri S, Aibara M, Terai Y, Okada N, Tachida H. High prevalence of non-synonymous substitutions in mtDNA of cichlid fishes from Lake Victoria. GENE 2014, 552(2):239–245.

  70. Starr TN, Thornton JW. Epistasis in protein evolution. PROTEIN SCI. 2016;25(7):1204–18.

    CAS  PubMed Central  PubMed  Google Scholar 

  71. Rasmussen HN, Rasmussen FN. The epiphytic habitat on a living host: reflections on the orchid–tree relationship. BOT J LINN SOC. 2018;186(4):456–72.

    Google Scholar 

  72. Zhang Z, He Z, Xu S, Li X, Guo W, Yang Y, Zhong C, Zhou R, Shi S. Transcriptome analyses provide insights into the phylogeny and adaptive evolution of the mangrove fern genus Acrostichum. Sci Rep. 2016;6:35634.

    CAS  PubMed Central  PubMed  Google Scholar 

  73. Lloyd RM. Systematics of the genus Ceratopteris (Parkeriaceae), I. sexual and vegetative reproduction in hawaiian Ceratopteris thalictroides. AM FERN J. 1973;69(1):12–8.

    Google Scholar 

  74. Peng L, Yamamoto H, Shikanai T. Structure and biogenesis of the chloroplast NAD(P)H dehydrogenase complex. Biochim Biophys Acta. 2011;1807(8):945–53.

    CAS  PubMed  Google Scholar 

  75. Wang P, Duan W, Takabayashi A, Endo T, Shikanai T, Ye JY, Mi H. Chloroplastic NAD(P)H dehydrogenase in Tobacco leaves functions in alleviation of oxidative damage caused by temperature stress. PLANT PHYSIOL. 2006;141(2):465–74.

    CAS  PubMed Central  PubMed  Google Scholar 

  76. Sommer F, Drepper F, Haehnel W, Hippler M. The hydrophobic recognition site formed by residues PsaA-Trp651 and PsaB-Trp627 of photosystem I in Chlamydomonas reinhardtii confers distinct selectivity for binding of plastocyanin and cytochrome c6. J BIOL CHEM. 2004;279(19):20009–17.

    CAS  PubMed  Google Scholar 

  77. Rogalski M, Schöttler MA, Thiele W, Schulze WX, Bock R. Rpl33, a nonessential plastid-encoded ribosomal protein in Tobacco, is required under cold stress conditions. Plant Cell. 2008;20(8):2221–37.

    CAS  PubMed Central  PubMed  Google Scholar 

  78. Jin JJ, Yu WB, Yang JB, Song Y, DePamphilis CW, Yi TS, Li DZ. GetOrganelle: a fast and versatile toolkit for accurate de novo assembly of organelle genomes. GENOME BIOL. 2020;21(1):241.

    PubMed Central  PubMed  Google Scholar 

  79. Dierckxsens N, Mardulyn P, Smits G. NOVOPlasty: de novo assembly of organelle genomes from whole genome data. NUCLEIC ACIDS RES. 2017;45(4):e18.

    PubMed  Google Scholar 

  80. Katoh K, Standley DM. MAFFT multiple sequence alignment software version 7: improvements in performance and usability. MOL BIOL EVOL. 2013;30(4):772–80.

    CAS  PubMed Central  PubMed  Google Scholar 

  81. Zhang D, Gao F, Jakovlic I, Zou H, Zhang J, Li WX, Wang GT. PhyloSuite: an integrated and scalable desktop platform for streamlined molecular sequence data management and evolutionary phylogenetics studies. MOL ECOL RESOUR. 2020;20(1):348–55.

    PubMed  Google Scholar 

  82. Stamatakis A. RAxML version 8: a tool for phylogenetic analysis and post-analysis of large phylogenies. Bioinformatics. 2014;30(9):1312–3.

    CAS  PubMed Central  PubMed  Google Scholar 

  83. Ronquist F, Teslenko M, van der Mark P, Ayres DL, Darling A, Hohna S, Larget B, Liu L, Suchard MA, Huelsenbeck JP. MrBayes 3.2: efficient bayesian phylogenetic inference and model choice across a large model space. SYST BIOL. 2012;61(3):539–42.

    PubMed Central  PubMed  Google Scholar 

  84. Storey JD, Tibshirani R. Statistical significance for genomewide studies. Proc Natl Acad Sci U S A. 2003;100(16):9440–5.

    CAS  PubMed Central  PubMed  Google Scholar 

  85. Wertheim JO, Murrell B, Smith MD, Kosakovsky Pond SL, Scheffler K. RELAX: detecting relaxed selection in a phylogenetic framework. MOL BIOL EVOL. 2015;32(3):820–32.

    CAS  PubMed  Google Scholar 

  86. Kelley LA, Mezulis S, Yates CM, Wass MN, Sternberg MJ. The Phyre2 web portal for protein modeling, prediction and analysis. NAT PROTOC. 2015;10(6):845–58.

    CAS  PubMed Central  PubMed  Google Scholar 

Download references

Acknowledgements

We would like to thank Shenzhen Fairy Lake Botanical Garden for their support in providing the plant materials.

Funding

This work was supported by the National Natural Science Foundation of China (31872670 and 32071781), Guangdong Basic and Applied Basic Research Foundation (2021A1515010911), Science and Technology Projects in Guangzhou (202206010107), and Project of Department of Science and Technology of Shenzhen City, Guangdong, China (JCYJ20190813172001780 and JCYJ20210324141000001).

Author information

Authors and Affiliations

  1. College of Life Sciences, South China Agricultural University, Guangzhou, 510642, China

    Xiaolin Gu, Lingling Li, Wanxin Shi, Xiaona Zhong & Ting Wang

  2. College of Natural Resources and Environment, South China Agricultural University, Guangzhou, 510642, China

    Sicong Li

  3. School of Life Sciences, Sun Yat-sen University, Guangzhou, 510275, China

    Yingjuan Su

  4. Research Institute of Sun Yat-sen University in Shenzhen, Shenzhen, 518057, China

    Yingjuan Su

Authors
  1. Xiaolin Gu
    View author publications

    You can also search for this author in PubMed Google Scholar

  2. Lingling Li
    View author publications

    You can also search for this author in PubMed Google Scholar

  3. Sicong Li
    View author publications

    You can also search for this author in PubMed Google Scholar

  4. Wanxin Shi
    View author publications

    You can also search for this author in PubMed Google Scholar

  5. Xiaona Zhong
    View author publications

    You can also search for this author in PubMed Google Scholar

  6. Yingjuan Su
    View author publications

    You can also search for this author in PubMed Google Scholar

  7. Ting Wang
    View author publications

    You can also search for this author in PubMed Google Scholar

Contributions

Data analyses, visualization and curation were conducted by X.G. Sample collection were conducted by X.G., L.L and S.L. Funding and supervision were contributed by Y.S. and T.W. X.G., Y.S., T.W., L.L., S.L., W.S. and X.Z. all contributed towards investigation and manuscript complete. X.G. wrote the manuscript together with Y.S. and T.W. The study was conceptualized by X.G. and T.W. All authors read and approved the final manuscript.

Corresponding authors

Correspondence to Yingjuan Su or Ting Wang.

Ethics declarations

Ethics approval and consent to participate

Not applicable.

Consent for publication

Not applicable.

Competing interests

The authors declare no competing interests.

Additional information

Publisher’s Note

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

Electronic supplementary material

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 http://creativecommons.org/licenses/by/4.0/. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) 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

Gu, X., Li, L., Li, S. et al. Adaptive evolution and co-evolution of chloroplast genomes in Pteridaceae species occupying different habitats: overlapping residues are always highly mutated. BMC Plant Biol 23, 511 (2023). https://doi.org/10.1186/s12870-023-04523-1

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s12870-023-04523-1

Keywords

BMC Plant Biology

ISSN: 1471-2229

Contact us