
==== Front
Ecol Evol
Ecol Evol
10.1002/(ISSN)2045-7758
ECE3
Ecology and Evolution
2045-7758
John Wiley and Sons Inc. Hoboken

10.1002/ece3.70318
ECE370318
ECE-2024-01-00095.R2
Ecological Genetics
Evolutionary Ecology
Genomics
Research Article
Research Article
New insights on the phylogeny, evolutionary history, and ecological adaptation mechanism in cycle‐cup oaks based on chloroplast genomes
Li et al.
Li Yu 1 2
Zheng Si‐Si 1
Wang Tian‐Rui 1
Liu Mei‐Hua 2
Kozlowski Gregor https://orcid.org/0000-0003-4856-2005
1 3 4
Yi Li‐Ta 2 yilita@zafu.edu.cn

Song Yi‐Gang https://orcid.org/0000-0003-2584-2338
1 2 ygsong@cemps.ac.cn

1 Eastern China Conservation Centre for Wild Endangered Plant Resources Shanghai Chenshan Botanical Garden Shanghai China
2 College of Forestry and Biotechnology Zhejiang A&F University Hangzhou China
3 Department of Biology and Botanic Garden University of Fribourg Fribourg Switzerland
4 Natural History Museum Fribourg Fribourg Switzerland
* Correspondence
Yi‐Gang Song, Eastern China Conservation Centre for Wild Endangered Plant Resources, Shanghai Chenshan Botanical Garden, Shanghai 201602, China.
Email: ygsong@cemps.ac.cn
Li‐Ta Yi, College of Forestry and Biotechnology, Zhejiang A&F University, Lin'an, Hangzhou 311300, China.
Email: yilita@zafu.edu.cn

17 9 2024
9 2024
14 9 10.1002/ece3.v14.9 e7031828 8 2024
19 1 2024
03 9 2024
© 2024 The Author(s). Ecology and Evolution published by John Wiley & Sons Ltd.
https://creativecommons.org/licenses/by/4.0/ This is an open access article under the terms of the http://creativecommons.org/licenses/by/4.0/ License, which permits use, distribution and reproduction in any medium, provided the original work is properly cited.

Abstract

Cycle‐cup oaks (Quercus section Cyclobalanopsis) are one of the principal components of forests in the tropical and subtropical climates of East and Southeast Asia. They have experienced relatively recent increases in the diversification rate, driven by changing climates and the Himalayan orogeny. However, the evolutionary history and adaptive mechanisms at the chloroplast genome level in cycle‐cup oaks remain largely unknown. Therefore, we studied this problem by conducting chloroplast genomics on 50 of the ca. 90 species. Comparative genomics and other analyses showed that Quercus section Cyclobalanopsis had a highly conserved chloroplast genome structure. Highly divergent regions, such as the ndhF and ycf1 gene regions and the petN—psbM and rpoB—trnC‐GCA intergenic spacer regions, provided potential molecular markers for subsequent analysis. The chloroplast phylogenomic tree indicated that Quercus section Cyclobalanopsis was not monophyletic, which mixed with the other two sections of subgenus Cerris. The reconstruction of ancestral aera inferred that Palaeotropics was the most likely ancestral range of Quercus section Cyclobalanopsis, and then dispersed to Sino‐Japan and Sino‐Himalaya. Positive selection analysis showed that the photosystem genes had the lowest ω values among the seven functional gene groups. And nine protein‐coding genes containing sites for positive selection: ndhA, ndhD, ndhF, ndhH, rbcL, rpl32, accD, ycf1, and ycf2. This series of analyses together revealed the phylogeny, evolutionary history, and ecological adaptation mechanism of the chloroplast genome of Quercus section Cyclobalanopsis in the long river of earth history. These chloroplast genome data provide valuable information for deep insights into phylogenetic relationships and intraspecific diversity in Quercus.

Based on 50 plastomes of cycle‐cup oaks, we found a conserved structure while nine protein‐coding genes went through positive selection. The phylogenetic tree based on complete chloroplast genomes, coding DNA sequences, and highly divergent regions produced a poorly resolved genetic relationship of the plastid lineages of cycle‐cup oaks.

adaptive evolution
chloroplast genome
phylogenetic relationships
Quercus section Cyclobalanopsis
Special Fund for Scientific Research of Shanghai Landscaping & City Appearance Administrative BureauG192422 G222404 National Natural Science Foundation of China 10.13039/501100001809 31901217 source-schema-version-number2.0
cover-dateSeptember 2024
details-of-publishers-convertorConverter:WILEY_ML3GV2_TO_JATSPMC version:6.4.8 mode:remove_FC converted:17.09.2024
Li, Y. , Zheng, S.‐S. , Wang, T.‐R. , Liu, M.‐H. , Kozlowski, G. , Yi, L.‐T. , & Song, Y.‐G. (2024). New insights on the phylogeny, evolutionary history, and ecological adaptation mechanism in cycle‐cup oaks based on chloroplast genomes. Ecology and Evolution, 14 , e70318. 10.1002/ece3.70318
==== Body
pmc1 INTRODUCTION

Cycle‐cup oaks (Quercus section Cyclobalanopsis) are an important taxon of tropical and subtropical evergreen broad‐leaved forests in East and Southeast Asia (Denk et al., 2017; Denk & Grimm, 2010; Huang et al., 1999). Most cycle‐cup oaks are large trees, and only a few species living in montane cloud forests are shrubs (Li et al., 2022; Song et al., 2019; Wang et al., 2024). Morphologically, the leaf epidermal features reveal three main groups in cycle‐cup oaks that are the compound trichome base (CTB), branched uniseriate (BU), and single‐celled trichome base (STB) groups (Deng et al., 2014). Owing to the concentric circle formation of the shell bracts, cycle‐cup oaks (ca. 90 species) have been treated as a separate genus Cyclobalanopsis or Quercus subgenus Cyclobalanopsis for over a hundred years (Camus, 1936; Huang et al., 1999; Luo & Zhou, 2001; Menitsky, 2005; Nixon, 1993; Ørsted, 1871; Schwarz, 1936). In the recent unified systematic scheme for oaks based on morphological traits, molecular phylogenetic relationships, and evolutionary history, Cyclobalanopsis is downgraded to a section belonging to Quercus subgenus Cerris (Deng et al., 2018; Denk et al., 2017; Hipp et al., 2020).

With the development of molecular systematics, researchers have used genetics from different components (nDNA, cpDNA, and mtDNA) to better understand the phylogenetic relationships of Quercus (Li et al., 2021; Liu et al., 2022; Manos et al., 1999; McVay et al., 2017; Yang et al., 2016). However, most of these studies focus on the subgenus Quercus, while there are relatively few comprehensive phylogenetic studies on section Cyclobalanopsis. Based on restriction‐site‐associated DNA sequencing data of 35 species of Cyclobalanopsis, the phylogeny of cycle‐cup oaks was monophyletic and clearly divided into two main lineages (CTB and STB). Moreover, biogeographic reconstruction showed that cycle‐cup oaks originated in Palaeotropics, then spread and transferred to Sino‐Japan and Sino‐Himalaya during the Miocene (Deng et al., 2018; Hipp et al., 2020). This is the first comprehensive study on the phylogenetic and evolutionary history of section Cyclobalanopsis at the level of the nuclear genome inherited by both parents. However, the number of species selected in this study is still insufficient (only one‐third of the total number of species), and there is a lack of data on maternal inheritance to support the results.

Chloroplast DNA is a valid maternal genetic data for the phylogenetic analysis of Quercus (Petit et al., 2002; Simeone et al., 2016; Yang et al., 2017). Phylogenetic reconstruction of 147 individuals from 29 Cyclobalanopsis species was performed using four chloroplast DNA markers. The results showed that Quercus section Cyclobalanopsis was divided into one main clade and 11 small subclades, which were mixed with section Ilex (Yan et al., 2018). However, the results do not support monophyly in Cyclobalanopsis and the resolution is poor on some nodes. Chloroplast genomes have a typical circular quadripartite structure, generally ranging from 120 to 160 kb, and are highly conserved in terms of structure, size, and gene content (Marechal & Brisson, 2010; Palmer, 1985; Shaw et al., 2007). The whole chloroplast genome can improve the phylogenetic resolution of species to better resolve the phylogenetic relationships (Gitzendanner et al., 2018; He et al., 2021; Li, Yi et al., 2019; Tu et al., 2021; Zhai et al., 2019). Phylogenetic analysis of Quercus based on the whole chloroplast genomes also found that cycle‐cup oaks nested in Quercus section Ilex (Liu et al., 2021; Yang, Zhou et al., 2021; Yang, Qu et al, 2021). However, due to the minimal species coverage of cycle‐cup oaks, the phylogeny of the whole chloroplast genome in this section is not clear.

In this study, we sequenced and assembled the complete chloroplast genomes of 36 species of Quercus section Cyclobalanopsis. Additionally, 14 chloroplast genomes of Quercus section Cyclobalanopsis species were also obtained from the National Center for Biotechnology Information (NCBI) database. In total, our dataset comprised 50 of the approximate 90 species currently recognized, representing 55% of Quercus section Cyclobalanopsis. Using this dataset, we aimed to explore the structure and variations of the chloroplast genomes, reveal its phylogeny and evolution in the long river of earth history, and elucidate the adaptive mechanism of protein‐coding genes under positive selection. The findings will improve our overall understanding of the classification, phylogeny, and evolution of Quercus section Cyclobalanopsis.

2 MATERIALS AND METHODS

2.1 Taxon sampling, DNA extraction, and sequencing

Fresh and healthy plant leaves from 36 species were newly sampled from different distribution sites and desiccated in silica gel (Figure S1; Table S1). These specimens were formally identified by Dr. Yi‐Gang Song, and deposited in the Chenshan Herbarium (CSH) at Shanghai Chenshan Botanical Garden for subsequent DNA extraction. In addition, 14 chloroplast genomes of non‐duplicating species, previously published in NCBI, were also directly downloaded (Accessed on July 1 2022) (Table S1) (Cho et al., 2021; Ju et al., 2019; Li, Wang et al., 2019; Li, Wang, Liu et al.,2020; Li, Wang, Zhao et al.,2020; Li, Luo et al., 2021; Li et al., 2022; Su et al., 2019; Wang et al., 2021; Yang et al., 2018).

Total plant genomic DNA was extracted using a modified cetyl trimethyl ammonium bromide protocol (Doyle, 1987). The extracted DNA was used to construct DNA sequencing libraries with an average insert size of 350 bp, based on the whole‐genome shotgun strategy (Batzoglou et al., 1999). Paired‐end (PE) sequencing of the libraries was then performed using an Illumina NovaSeq 6000 (San Diego, CA, USA), according to the manufacturer's protocol, at Wuhan Benagen Technology Co., Ltd. (Wuhan, China). Raw data were converted using base calling analysis with reads of 150 bp in length (Ewing & Green, 1998).

2.2 Sequence assembly and annotation

Low‐quality reads from raw data were filtered using the SOAPnuke v1.3.0 software to obtain clean data (Chen et al., 2018). Based on clean data, the chloroplast genomes of 36 newly sequenced species in this study were de novo assembled using the GetOrganelle v1.7.6.1 software (Jin et al., 2020). All 36 newly sequenced and 14 downloaded chloroplast genomes were annotated using the online annotation tool GeSeq (https://chlorobox.mpimpgolm.mpg.de/geseq.html) (Tillich et al., 2017), and manually checked and adjusted. Chloroplast genome maps were drawn using the online program OGDRAW v1.3.1 (https://chlorobox.mpimp‐golm.mpg.de/OGDraw.html) (Greiner et al., 2019). The online program CPGView (http://47.96.249.172:16085/cpgview/home) was used to view the genetic structure of the spliced genes in the chloroplast genomes (Liu et al., 2023). Finally, the newly obtained genome sequences were uploaded and saved to the NCBI database under GenBank accession numbers OQ603420–OQ603455 (Table S1), and the raw sequence data were deposited in the Genome Sequence Archive in National Genomics Data Center, China National Center for Bioinformation/Beijing Institute of Genomics, Chinese Academy of Sciences (GSA: CRA017174) that were publicly accessible at https://ngdc.cncb.ac.cn/gsa (Chen et al., 2021; CNCB‐NGDC Members and Partners, 2022).

2.3 Genome structure analyses

The basic structural characteristics of the chloroplast genomes, such as genome length, GC content, and gene number, were calculated using the Geneious R9.0.2 software (Kearse et al., 2012). Simple sequence repeats (SSRs) with one to six base pairs as the basic repeat unit were identified using the MISA Perl script (https://webblast.ipk‐gatersleben.de/misa/) (Beier et al., 2017). The thresholds for repeat units of the SSRs were set to ten for mono‐; five for di‐; four for tri‐; and three for tetra‐, penta‐, and hexanucleotides. The online program Tandem Repeats Finder (TRF) v4.09 (https://tandem.bu.edu/trf/trf.html) was used to detect minisatellite sequence repeats (M) with repeat units greater than 10 bp (Benson, 1999). The alignment parameters for match, mismatch, and indels were set as two, seven, and seven, respectively, and the minimum alignment score and maximum period size were set to 80 and 500, respectively. REPuter website (https://bibiserv.cebitec.uni‐bielefeld.de/reputer) was used to predict four dispersed repeat sequences: forward (F), reverse (R), complementary (C), and palindromic (P) (Kurtz et al., 2001). Thirty and three were set as the minimum repeat size and hamming distance, respectively (Liang et al., 2019). For non‐repeating sequences and using the screening criteria of a length (≥300 bp) and the presence of both initial (ATG) and termination codons (TAA, TGA, or TAG), a total of 51 common coding sequences (CDS) were selected for codon bias analysis by the CodonW v1.4.2 software (Wright, 1990).

2.4 Comparative genome analyses

The boundary genes between the inverted repeat (IR) and single copy (SC) regions of each chloroplast genome were compared using the online tool IRscope (https://irscope.shinyapps.io/irapp/), which was specifically designed for the chloroplast genome (Amiryousefi et al., 2018). Fifty chloroplast genomes were aligned to assess possible structural changes, such as gene order rearrangements, inversions, and insertions, using the progressive mauve algorithm of the Mauve v2.4.0 software (Darling et al., 2004). Using the chloroplast genome of Quercus glauca as a reference genome, the mVISTA software (https://genome.lbl.gov/vista/mvista/submit.shtml) was used to compare interspecific variations in the 50 species at the chloroplast genome level in the Shuffle‐LAGAN mode (Brudno et al., 2003; Frazer et al., 2004). Sliding window analysis was performed using the DnaSP (DNA Sequence Polymorphism) v6.12.03 software to calculate the nucleotide variability (Pi) values and detect highly variable sites among the chloroplast genomes with step size of 200 bp and window length of 800 bp (Rozas et al., 2017).

2.5 Phylogenetic analyses

To infer the phylogenetic relationships of Quercus section Cyclobalanopsis, we performed phylogenetic analyses using three data sets: (1) the complete chloroplast genomes (CCGs); (2) the CDS of PCGs; and (3) the highly divergent regions detected in comparative genomic analyses (gene regions psbC, ndhF, and ycf1; and intergenic spacer regions petN—psbM, psbZ—trnG‐GCC, rpl32—trnL‐UAG, trnK‐UUU—rps16, trnH‐GUG—psbA, rpoB—trnC‐GCA, trnF‐GAA—ndhJ, and ndhF—rpl32). Quercus serrata (MK922350), Quercus acutissima (MH899015), and Quercus aquifolioides (KX911971) were chosen as outgroups to root the phylogenetic trees. Multiple sequences were aligned using the MAFFT v7.450 software with default parameters (Katoh & Standley, 2013). Maximum likelihood (ML) phylogenetic trees were constructed using the IQ‐tree v2.1.3 software with 1000 ultrafast bootstrap replicates (Minh et al., 2020). The best‐fitting models were selected using ModelFinder in IQ‐tree v2.1.3 (Kalyaanamoorthy et al., 2017). Finally, the constructed phylogenetic trees were visualized and beautified in the FigTree v.1.4.4 software (http://tree.bio.ed.ac.uk/software/figtree/).

2.6 Ancestral aera reconstruction

Based on Flora of China, the geographical distribution information of 50 species of Quercus section Cyclobalanopsis and outgroups were obtained in Table S2 (Huang et al., 1999). According to the floristic regions divided by Wu et al. (2011), it was divided into three biogeographical regions: (A) Sino‐Japan (including Eastern Asia subtropical areas east of the Tanaka Line); (B) Sino‐Himalaya (including the subtropical areas west of the Tanaka Line and the southeast edge of the Himalayas); and (C) Palaeotropics (including tropical regions of southern Yunnan, Guangxi, Hainan, and Indo‐China) (Deng et al., 2018; Wu et al., 2011). The statistical dispersal‐vicariance analysis (S‐DIVA) in the RASP v4.3 software was employed to reconstruct the ancestral geographical distribution of Quercus section Cyclobalanopsis (Yu et al., 2015, 2020). The analysis was performed on the phylogenetic tree constructed from CCGs with the highest support value.

2.7 Positive selection analyses

To explore the adaptive evolution of PCGs in the chloroplast genomes of Quercus section Cyclobalanopsis, positive selection analyses were performed using the Codeml program in the PAML v4.9j software (Yang, 2007). CDS were stripped of the termination codon and aligned using the ClustalW (Codons) plugins in the MEGA‐X software (Kumar et al., 2018; Thompson et al., 2002). An ML phylogenetic tree was generated based on the CDS of 50 chloroplast genomes using the IQ‐tree v2.1.3 software (Minh et al., 2020). Six site models (seqtype = 1; model = 0; and NSsites = 0, 1, 2, 3, 7, and 8) were used to identify the selection pressure and potential positive selection sites on the 79 common PCGs. The likelihood ratio test (LRT) was used to compare models between M1 (nearly neutral) and M2 (positive selection), M0 (one ratio) and M3 (discrete), and M7 (beta) and M8 (beta and ω). The chi‐square test (χ 2) was used to detect significance for the PCGs with positive selection sites, and genes with p < .05 were selected as positive selection genes. Finally, the posterior probability of the sites was calculated based on Bayes empirical Bayes (BEB) to evaluate the significance level of positive selection sites (p > 95%).

3 RESULTS

3.1 Basic structural characteristics of chloroplast genomes

The length of the assembled chloroplast genomes of Quercus section Cyclobalanopsis ranged from 160,098 bp (Quercus lamellosa) to 161,914 bp (Quercus brandisiana) (Figure 1). All chloroplast genomes were highly conserved and showed a typical circular quadripartite structure comprising a large single copy (LSC) region (89,725–90,577 bp) and a small single copy (SSC) region (18,175–18,996 bp) separated by a pair of IR regions (25,748–26,581 bp) (Figure 1c; Table S3). The average GC content (GC%) varied from 36.73% to 36.93%, with that of the IR region (42.32%–42.83%) being significantly higher than that of the LSC (34.57%–34.81%) and SSC regions (30.72%–31.24%) (Figure 1d; Table S3).

FIGURE 1 Schematic map and basic features of the chloroplast genomes of Quercus section Cyclobalanopsis. (a) Genome map of the chloroplast genomes. Genes outside the circle are transcribed in the counter‐clockwise direction, whereas those inside the circle are transcribed in the clockwise direction. The genes belonging to different functional groups are identified by different colors. The length and boundary of the LSC, SSC, and two IR regions are indicated in the inner circle. The dark gray area indicates the GC content while the lighter gray corresponds to the AT content of the genomes. IR, inverted repeat; LSC, large single copy; SSC, small single copy. (b) Gene structure map of the cis‐spliced genes and the trans‐spliced gene rps12 in the chloroplast genomes. The cis‐spliced genes are arranged from top to bottom in the order they are found in the chloroplast genomes. The gene names are shown on the left, and the gene structures are on the right. The exons are shown in black and the introns are shown in white. The arrows indicate the direction of the genes. (c, d) The length (c) and GC content (d) of each region of the chloroplast genomes.

The total number of annotated genes was uniform among the 50 chloroplast genomes, including 86 PCGs, 37 transfer RNA (tRNAs), and eight ribosomal RNA genes (rRNAs). Eighteen genes, comprising seven PCGs, seven tRNAs, and four rRNAs, were duplicated in the two IR regions (Figure 1a; Table 1). The rps12 and ycf1 genes spanned two regions between the IR and LSC/SSC, respectively (Figure 1a). The 113 unique genes were divided into four categories according to their functions, with transcription and translation, photosynthesis, biosynthesis, and unknown containing 60, 44, 5, and 4 genes, respectively. A total of 18 genes contained introns, among which 15 genes, trnK‐UUU, trnG‐GCC, trnL‐UAA, trnV‐UAC, trnI‐GAU, trnA‐UGC, rps16, rpl2, rpl16, rpoC1, atpF, ndhA, ndhB, petB, and petD, contained one intron, and the other three genes, ycf3, clpP, and rps12, contained two introns (Table 1). Of the 12 PCGs spliced by introns, 11 were cis‐spliced genes, whereas rps12 was a trans‐spliced gene (Figure 1b).

TABLE 1 Genetic classification of the chloroplast genomes of Quercus section Cyclobalanopsis.

Category	Group	Name	
Transcription and translation	Translational initiation factor	infA	
Ribosomal RNAs	rrn16S(×2), rrn4.5S(×2), rrn23S(×2), rrn5S(×2)	
Transfer RNAs	trnR‐UCU, trnfM‐CAU, trnD‐GUC, trnH‐GUG, trnM‐CAU, trnE‐UUC, trnS‐GGA, trnF‐GAA, trnP‐UGG, trnT‐UGU, trnG‐UCC, trnQ‐UUG, trnY‐GUA, trnW‐CCA, trnS‐UGA, trnC‐GCA, trnT‐GGU, trnL‐UAG, trnK‐UUU*, trnV‐UAC*, trnL‐UAA*, trnG‐GCC*, trnA‐UGC*(×2), trnI‐GAU*(×2), trnL‐CAA(×2), trnI‐CAU(×2), trnN‐GUU(×2), trnV‐GAC(×2), trnR‐ACG(×2)	
Small subunit of ribosome (SSU)	rps2, rps11, rps19, rps14, rps4, rps15, rps16*, rps8, rps18, rps3, rps12(×2)**, rps7(×2)	
Large subunit of ribosome (LSU)	rpl14, rpl20, rpl36, rpl33, rpl16*, rpl32, rpl22, rpl2*(×2), rpl23(×2)	
DNA‐dependent RNA polymerase	rpoC2, rpoB, rpoC1*, rpoA	
Photosynthesis	Photosystem I	psaB, psaJ, psaA, psaI, psaC	
Photosystem II	psbA, psbC, psbH, psbZ, psbI, psbJ, psbK, psbF, psbD, psbT, psbN, psbL, psbM, psbE, psbB	
Subunit of cytochrome	petB*, petN, petL, petG, petD*, petA	
ATP synthase	atpA, atpI, atpB, atpE, atpF*, atpH	
RubisCO large subunit	rbcL	
NADH dehydrogenase	ndhG, ndhD, ndhE, ndhK, ndhH, ndhI, ndhF, ndhA*, ndhJ, ndhC, ndhB*(×2)	
Biosynthesis	Maturase	matK	
ATP‐dependent protease	clpP**	
Acetyl‐CoA‐carboxylase	accD	
Envelop membrane protein	cemA	
C‐Type cytochrome synthesis	ccsA	
Unknown	Hypothetical chloroplast reading frames (ycf)	ycf4, ycf3**, ycf1(×2), ycf2(×2)	
Note: Genes marked with the asterisk are those with single (*) or double (**) introns. Duplicate genes located in the IR regions are marked as (×2).

3.2 Repeat sequences

The number and distribution of SSRs, minisatellites, and dispersed repeat sequences were analyzed in the chloroplast genomes of the 50 species. SSRs were most abundant, and all types of repeats were present. Quercus macrocalyx contained the most SSRs (134), whereas Quercus austrocochinchinensis had the lowest (110) (Figure 2a; Table S4). The amount of mononucleotide repeats was the highest (74–88), accounting for approximately 70% of the total, followed by dinucleotide repeats; with slightly more tetranucleotide than trinucleotide repeats. Hexanucleotide repeats were only present in Quercus gambleana, Quercus langbianensis, Q. macrocalyx, and Quercus augustinii (Figure 2b). Most SSRs comprised two complementary bases, adenine (A) and thymine (T), indicating a strong A/T bias (Figure 2c). SSRs were distributed mainly in the LSC and intergenic spacer regions (IGS), with only a few in the SSC, IRs, intron, and coding regions (Figure 2d). A total of 407 minisatellite repeat sequences (M) were detected, ranging from six in Quercus pachyloma to 11 in Quercus oxyodon and Quercus kiukiangensis (Figure 3a; Table S5). The copy number of the repeat sequences was mainly between two and four, and the period size of the repeat units ranged from 11 to 101 bp and were mainly between 20 and 29 bp (Figure 3b). Among the four types of dispersed repeat sequences, palindromic (P) and forward (F) sequences were the most abundant, whereas reverse (R) and complementary (C) sequences were rare and were not detected in some species (Figure 3a; Table S5). All dispersed repeat sequences were between 30 and 120 bp in size, and most were concentrated between 30 and 40 bp (Figure 3c).

FIGURE 2 Type, distribution, and number of SSRs in 50 chloroplast genomes of Quercus section Cyclobalanopsis. (a) Total SSR numbers in each chloroplast genome. (b) The percentage of SSRs with six repeat types (mono‐, di‐, tri‐, tetra‐, penta‐, and hexanucleotide). (c) The number of SSRs with different base types. (d) The number of SSRs in different regions of the chloroplast genomes. CDS, coding sequences; IGS, intergenic spacers.

FIGURE 3 The number and length distribution of various repeat types in 50 chloroplast genomes of Quercus section Cyclobalanopsis. (a) Number of different repeat types in each chloroplast genome. C, complementary repeat sequences; F, forward repeat sequences; M, minisatellite repeat sequences; P, palindromic repeat sequences; R, reverse repeat sequences. (b, c) The length distribution of minisatellite (b) and dispersed (c) repeat sequences.

3.3 Codon usage bias

The codon usage bias analysis of 51 selected coding sequences showed that each sequence contained 20,263 (Q. lamellosa) to 20,440 (Quercus acuta) codons, and the mean GC content of codons (GC_all) ranged from 37.72% (Quercus rex) to 37.89% (Quercus morii) (Table S6). The GC content of the first, second, and third codon sites, represented by GC1, GC2, and GC3, respectively, was less than 50% and showed a decreasing trend of GC1 > GC2 > GC3, further indicating that the chloroplast genomes were rich in A/T bases (Table S6). The effective number of codon values (ENc) were all significantly greater than 35, ranging from 49.79 (Q. macrocalyx) to 50.03 (Q. lamellosa), indicating weak codon usage bias in these chloroplast genomes (Table S6). Among the calculated relative synonymous codon usage values (RSCU) of the 59 synonymous codons, 30 codons had RSCU values greater than 1; among which only two ended in G/C (UCC and UUG), and the remaining 28 all ended in A/U with 16 ending in U and 11 ending in A. This showed that these codons have a bias in favor of the A/U endings (Table S7). The codon with the largest RSCU value was UUA, which encodes Leucine (Leu), followed by AGA, which encodes Arginine (Arg) (Figure 4).

FIGURE 4 The heat map of the relative synonymous codon usage (RSCU) in 50 chloroplast genomes of Quercus section Cyclobalanopsis. In the color scale, higher red values indicate higher RSCU values and lower blue values indicate lower RSCU values.

3.4 Comparative genome analyses of the chloroplast genomes

Fifteen different types of contraction and expansion of the IR regions were identified which resulted in variation in the chloroplast genomes of Quercus section Cyclobalanopsis (Figure 5; Table S8). The junction regions had the same relative position in all 50 chloroplast genomes. Both rps19 and trnH were located in the LSC region and were distributed near the LSC/IRb (JLB) and LSC/IRa (JLA) boundaries, respectively. The rps19 contracted by four bp and trnH expanded by 14 bp of the seven species of Type 2, whereas in the remaining species, rps19 contracted by 11 bp and trnH expanded by 16 bp. The total length of ndhF distributed at the SSC/IRb (JSB) boundary was 2243–2288 bp, with varying degrees of expansion or contraction at the boundary of the 11 species from Type 6 to Type 15. The ycf1 spanned the SSC and two IR regions, with the length of the repeats in the two IR regions being between 1045 and 1815 bp. The ycf1 at the JSB boundary expanded into the SSC region from eight to 72 bp, whereas the ycf1 at the SSC/IRa (JSA) boundary contracted into the SSC region from 3852 to 4622 bp.

FIGURE 5 Comparison of the junction regions (JLA, JLB, JSB, and JSA) of Quercus section Cyclobalanopsis. Genes are denoted by colored boxes. The numbers above the gene boxes indicate the distance between the end of the gene and the boundary sites.

Multiple alignment analysis of Quercus section Cyclobalanopsis by Mauve showed no genomic rearrangement events or inversion phenomena, indicating a good collinearity relationship (Figure S3). These results further indicate that the chloroplast genomes are highly conserved.

Using Q. glauca as the reference sequence, the mVISTA results showed that the chloroplast genomes of Quercus section Cyclobalanopsis had high sequence similarity (99.0%–99.9%). However, the noncoding and SC regions showed a higher level of divergence compared with the coding and IR regions (Figure 6 and Figure S4). High variation was observed in the exon regions of two PCGs, ndhF, and ycf1, and in the conserved noncoding regions of three intergenic spacer regions, petN—psbM, psbZ—trnG‐GCC, and rpl32—trnL‐UAG. The variation of ycf1 was particularly significant between the chloroplast genomes (Figure 6). These divergent regions contain abundant variation information which has the potential to be used to develop molecular markers as DNA barcodes for species authentication in Quercus section Cyclobalanopsis.

FIGURE 6 Visual local map (125–149 kb) of the alignment of Quercus section Cyclobalanopsis by the mVISTA. The gray arrows above show genes of the reference sequence, and the direction is forward or reverse. The position of the genome is shown on the horizontal axis at the bottom of each block. The alignment similarity percentages are shown on the right side of the graph (the vertical axis). Genome regions are color‐coded as exon, untranslated regions (UTR), mRNA, and conserved non‐coding sequences (CNS). The complete visualization comparison map is shown in Figure S4.

Using DnaSP, a total of 2445 polymorphic sites were detected in the genome, including 1403 singleton variable sites (SNP) and 1052 parsimony informative sites. The range of the nucleotide diversity value was 0–0.0717 and had an average of 0.00133. The nucleotide diversity values for the IR regions (0–0.0022, average = 0.00031) were lower than those for the LSC (0–0.0717, average = 0.00173) and SSC (0–0.0686, average = 0.00225) regions. Consistent with the results of the whole‐genome comparisons, the LSC and SSC regions were much more divergent than the IR regions, and the noncoding regions also varied more than the coding regions. Furthermore, the analysis detected eight highly divergent regions (Pi > 0.004), three of which were located in the gene regions, psbC, ndhF, and ycf1, and five in the intergenic spacer regions, trnK‐UUU—rps16, trnH‐GUG—psbA, rpoB—trnC‐GCA, trnF‐GAA—ndhJ, and ndhF—rpl32 (Figure 7). These highly divergent regions may have undergone rapid nucleotide substitutions during the evolution of the species, which is of great importance for phylogenetic analysis and species identification.

FIGURE 7 Sliding window analysis of 50 complete chloroplast genomes of Quercus section Cyclobalanopsis (window length: 600 bp; step size: 200 bp). The X‐axis represents nucleotide positions of the midpoint of the window and the Y‐axis represents the value of nucleotide diversity (Pi) per window.

3.5 Phylogenetic relationship of Quercus section Cyclobalanopsis

To infer the phylogenetic relationships within the Quercus section Cyclobalanopsis, we constructed phylogenetic trees using the ML method based on three datasets (Figure 8 and Figure S5). In all three trees, 10 (highly divergent regions), 13 (CDS), and 15 (CCGs) strongly supported clades were recognizable, and the phylogenetic relationships constructed from the CCGs data received relatively high support values. The topological structure of the three phylogenetic trees showed that Q. acutissima in section Cerris and Q. aquifolioides in section Ilex were mixed with section Cyclobalanopsis, and that the section Cyclobalanopsis was not monophyletic. The section Cyclobalanopsis differentiated step by step, and Q. macrocalyx, Q. lamellosa, and Q. brandisiana differentiated first and were located at the base of the phylogenetic trees. In the latest differentiation, section Cyclobalanopsis was clearly divided into two major evolutionary clades and the number of species included in this branch increased with the dataset (highly divergent regions > CDS > CCGs).

FIGURE 8 Phylogenetic tree constructed based on the CCGs and PCGs of Quercus section Cyclobalanopsis using ML methods. The bootstrap support values (BS) are labeled at the branch in the evolutionary tree, where BS less than 50% is represented by a “—.”

3.6 Ancestral aera reconstruction

Based on the S‐DIVA analysis, we reconstructed the ancestral distribution and revealed the origin and evolutionary history of the existing Quercus section Cyclobalanopsis (Figure 9). The results showed that Palaeotropics (C) was the most likely ancestral range of the species. It then dispersed to Sino‐Japan (A) and Sino‐Himalaya (B), with Sino‐Japan (A) being the most important diffusion region. The two major evolutionary clades in the latest differentiation based on the phylogenetic tree originated together in Sino‐Japan (A) and were mainly distributed in regions Sino‐Japan (A) and Palaeotropics (C) (except Quercus lobbii).

FIGURE 9 Ancestral distribution in Quercus section Cyclobalanopsis using S‐DIVA analysis based on RASP. (A) Sino‐Japan; (B) Sino‐Himalaya; and (C) Palaeotropics.

3.7 Positive selection evolution

Seventy‐nine common PCGs were detected in the selection pressure using the PAML site‐model. The results of the model M0 (model = 0, NSsites = 0) showed that the dS value of six genes was equal to 0, which meant that the ω value was meaningless. The PCGs were categorized into seven functional groups based on the classification: ATP synthase, NADPH dehydrogenase, cytochrome b/f complex, photosystem, ribosomal proteins, RNA polymerase, and other genes. Among them, the photosystem genes had the lowest ω values, while the ribosomal protein genes had the highest ω values (Figure 10; Table S9).

FIGURE 10 Boxplot of the values of the ω (dN/dS) in each functional gene group.

We also identified PCGs with positive selection sites in models M2 (28) and M8 (37); rejecting the null models M1 and M7, respectively. The LRT was performed for 37 PCGs with positive selection sites based on pair‐to‐pair comparisons of M0 versus M3, M1 versus M2, and M7 versus M8. Genes with p < .05 were selected as positive selection genes. Nine PCGs were affected by selection pressure, including four NADH dehydrogenase genes, ndhA, ndhD, ndhF, and ndhH; RubisCO large subunit gene rbcL; one ribosomal protein gene rpl32; acetyl‐CoA‐carboxylase gene accD; and two hypothetical chloroplast reading frame genes, ycf1, and ycf2 (Table 2). Based on BEB, there were 116 positive selection sites (M8) in these nine PCGs of which 33 had a p > 95% (Table 2 and Table S10). The ycf1 gene located in the IR region harbored the highest number of positive selection sites (51), including 14 significant positive selection sites.

TABLE 2 Likelihood ratio test (LRT) and positive selection sites under different models of chloroplast protein‐coding genes in Quercus section Cyclobalanopsis.

Gene	Model comparison	df	ΔlnL	2ΔlnL	LRT (p‐value)	Number of positively selected sites (P > 95%)	
accD	M0 vs. M3	4	14.456206	28.912412	8.14E‐06		
M1 vs. M2	2	7.046054	14.092108	8.71E‐04	10 (4)	
M7 vs. M8	2	7.290798	14.581596	6.82E‐04	11 (4)	
ndhA	M0 vs. M3	4	7.140073	14.280146	6.45E‐03		
M1 vs. M2	2	4.464701	8.929402	1.15E‐02	5 (1)	
M7 vs. M8	2	4.469385	8.938770	1.15E‐02	7 (1)	
ndhD	M0 vs. M3	4	15.944651	31.889302	2.02E‐06		
M1 vs. M2	2	8.332762	16.665524	2.41E‐04	7 (1)	
M7 vs. M8	2	8.590458	17.180916	1.86E‐04	10 (2)	
ndhF	M0 vs. M3	4	13.504137	27.008274	1.98E‐05		
M1 vs. M2	2	4.845863	9.691726	7.86E‐03	8 (2)	
M7 vs. M8	2	5.587439	11.174878	3.74E‐03	15 (2)	
ndhH	M0 vs. M3	4	6.127339	12.254678	1.56E‐02		
M1 vs. M2	2	3.169165	6.338330	4.20E‐02	3 (1)	
M7 vs. M8	2	3.198039	6.396078	4.01E‐02	5 (1)	
rbcL	M0 vs. M3	4	41.132773	82.265546	5.77E‐17		
M1 vs. M2	2	25.996025	51.992050	5.13E‐12	4 (4)	
M7 vs. M8	2	26.239243	52.478486	4.02E‐12	4 (4)	
rpl32	M0 vs. M3	4	9.891909	19.783818	5.50E‐04		
M1 vs. M2	2	3.854424	7.708848	2.12E‐02	1 (1)	
M7 vs. M8	2	4.135478	8.270956	1.56E‐02	1 (1)	
ycf1	M0 vs. M3	4	95.177045	190.354090	4.45E‐40		
M1 vs. M2	2	50.525122	101.050244	1.14E‐22	30 (10)	
M7 vs. M8	2	51.867464	103.734928	2.98E‐23	51 (14)	
ycf2	M0 vs. M3	4	55.836083	111.672166	3.20E‐23		
M1 vs. M2	2	42.528727	85.057454	3.39E‐19	12 (4)	
M7 vs. M8	2	38.003763	76.007526	3.13E‐17	12 (4)	

4 DISCUSSION

4.1 Chloroplast genome structure of Quercus section Cyclobalanopsis

In this study, we assembled and annotated the complete chloroplast genomes of 36 species of Quercus section Cyclobalanopsis, enriching the chloroplast genome database of section Cyclobalanopsis and genus Quercus. Similar to the chloroplast genomes of published species of section Cyclobalanopsis (Chen & Zhang, 2023; Li et al., 2022; Wang et al., 2021), the whole chloroplast genomes varied only slightly in size (160,098–161,914 bp) (Figure 1; Table S3) and were also highly similar in overall structure, gene order, and content (Figure 1). The GC content in the IR regions was significantly higher than that in the SC regions owing to the presence of unique rRNA genes (Yang et al., 2016; Zong et al., 2019). Compared with the chloroplast genomes of other sections of Quercus, the chloroplast genomes of the section Cyclobalanopsis were generally smaller in length, which may be one of the reasons why it has been distinguished from other Quercus species in the classification system (Liu et al., 2021; Pang et al., 2019; Yang et al., 2016; Yang, Zhou et al., 2021; Zhang et al., 2020).

A total of 131 genes, including 86 PCGs, 37 tRNAs, and eight rRNA genes, were annotated in the chloroplast genomes of 50 species (Table 2 and Table S3). Previous studies of the chloroplast genome sequences downloaded from the NCBI database have found that using different annotation software and reference genome sequences can lead to differences in the annotation results (Cho et al., 2021; Li, Wang, Liu et al., 2020; Yang et al., 2018). The ycf15 gene is defined as a pseudogene because it does not encode proteins; therefore, it is not annotated in most chloroplast genomes of Quercus. However, it has been annotated in a few species, such as Q. acutissima (MH607377) and Quercus fabri (MK693136) (Li et al., 2018; Li, Li et al., 2021). Therefore, it is necessary to compare and verify the results of genome annotation to reduce personal errors and ensure data reliability. Gene loss and transfer to the nucleus are the major feature of chloroplast genome evolution (Kleine et al., 2009; Stegemann et al., 2003). Consistent with most angiosperms, the results showed that the chloroplast genome of Quercus had almost no gene deletion or transfer. However, previous studies have found that some genes, such as ndhF, rpl22, and clpP, have been lost or transferred during the variation and evolution of chloroplast genomes in some species (Jansen et al., 2011; Yao et al., 2015).

4.2 Comparison and evolution of chloroplast genomes of Quercus section Cyclobalanopsis

Repeat sequences exist widely in plant genomes and may have played an important role in plant evolution (Timme et al., 2007; Weng et al., 2014). The repeat sequences detected in this study were composed of A and T bases and had a strong A/T bias (Figure 2), which is consistent with the results of other studies (Morton, 1997; Zhang et al., 2020). The length, quantity, and distribution of various repeat sequences were highly consistent and conserved, and these differences may be caused by the size and variation of each chloroplast genome (Wu et al., 2018). Different types of repeat sequences were detected in two gene regions, rpl12 and ycf2, which are likely involved in chloroplast genome rearrangement and sequence differentiation (Guisinger et al., 2010).

Codon usage bias is an important evolutionary characteristic that is common in many organisms. It is influenced by natural selection, base mutations, and other factors (Morton, 1998), and has been extensively studied in many plant species (Qin et al., 2013; Shi et al., 2022; Wang et al., 2018; Wen et al., 2021). In this study, the GC content decreased in the following order: GC1 > GC2 > GC3 (Table S6). GC content is a major factor in the formation of codon usage bias and may play an important role in the evolution of the genome structure (Yang et al., 2018). Among the 30 preferential synonymous codons, most of them terminated in A/T, except that UUG ended with G, and UCC ended with C (Figure 4). This demonstrated that the codon usages of the chloroplast genomes were biased towards A/T at the third position of codons, which was similar to reported results of other angiosperm, including the Bupleurum, Chamaesium, Polystachya, and Allium section Bromatorrhiza species (Chen et al., 2022; Guo et al., 2020; Jiang et al., 2022; Li, Xie et al., 2020).

The IR region plays an important role in stabilizing the chloroplast genome structure, and its expansion and contraction are the main causes of chloroplast genome length variation (Marechal & Brisson, 2010). In general, the distribution of boundary genes in Quercus section Cyclobalanopsis in this study did not show significant differences. However, the distribution of ndhF showed the greatest difference, being completely located in the SSC region in 39 species, whereas the remaining species had varying degrees of expansion or contraction at the JSB boundary (Figure 5). The same situation was observed in previous studies on Punica granatum, Allium, and other plants (Xie et al., 2020; Yan, Zhao et al.,2019). Chloroplast genomes have sequence homology; therefore, they have no gene rearrangement or inversion in most plant species (Li, Sylvester et al., 2019). In the present study, the chloroplast genome of Quercus section Cyclobalanopsis also showed good collinearity (Figure S3). Further comparison of the chloroplast genomes of Quercus section Cyclobalanopsis showed that the SC regions had a higher level of sequence diversity than the IR regions, which may be because the IR regions contained more conserved rRNA genes. This result was similar to that of other chloroplast genomes, such as Q. acutissima, Lonicera japonica, and Aceraceae (Han et al., 2016; He et al., 2017; Zhang et al., 2020). In addition, we identified highly divergent coding regions, psbC, ndhF, and ycf1, and noncoding regions, petN—psbM, trnK‐UUU—rps16, and rpoB—trnC‐GCA (Figures 6 and 7). These highly divergent regions are useful molecular resources for potential DNA barcoding and for subsequent population genetics studies of different species.

4.3 Phylogeny and evolutionary history of chloroplast genomes of Quercus section Cyclobalanopsis

Each angiosperm species contains three types of genomes, nuclear, chloroplast, and mitochondrial, all three of which have been used in phylogenetic analysis. Phylogenetic studies of Quercus are faced with great challenges, not only because of serious hybridization and introgression but also because of the inconsistent phylogenetic relationships constructed by different markers (Denk & Grimm, 2010; Manos et al., 1999; Simeone et al., 2013). In the present study, phylogenetic trees of Quercus section Cyclobalanopsis were constructed based on CCGs, PCGs, and highly divergent regions (Figure 8 and Figure S5). We also reconstructed the ancestral geographical distribution of Quercus section Cyclobalanopsis using the RASP (Figure 9).

First, some species had poor resolution between differentiation, which may be due to the low substitution rate and low genetic diversity of chloroplast genomes (Drouin et al., 2008; Yan et al., 2018; Yan, Liu, Li et al., 2019). Previous studies had shown that the species resolution of chloroplast markers in East Asian oak was poor (Ohyama et al., 2001; Yang et al., 2017), and the genetic diversity of section Cyclobalanopsis was the lowest in the subgenus Cerris (Yan et al., 2018; Yan,Liu, Li et al., 2019), suggesting that the chloroplast genome may not be able to infer a systematic classification of closely related species of Quercus section Cyclobalanopsis (Pham et al., 2017). Secondly, the species of Quercus section Cyclobalanopsis differentiated one by one during the evolutionary process and were obviously divided into two major evolutionary clades in the latest differentiation. This differed from the phylogenetic relationships previously constructed from the restriction‐site‐associated DNA sequencing data of 35 species of Quercus section Cyclobalanopsis which were divided into two main clades, named STB and CTB, respectively (Deng et al., 2018; Hipp et al., 2020). The species of each clade based on the chloroplast and nuclear genome were not consistent owing to hybridization and introgression. Maternally inherited chloroplast DNA was easily transferred during hybridization (Dumolin et al., 1995). In the reconstruction of ancestral distribution, it was revealed that the most probable ancestral distribution area of the existing species of Quercus section Cyclobalanopsis was Palaeotropics (C), and then dispersed to Sino‐Japan (A) and Sino‐Himalaya (B) (Figure 9). This result was consistent with the reconstruction of the ancestral region based on the nuclear genome (restriction‐site‐associated DNA) (Deng et al., 2018). However, the reconstruction of ancestral distribution based on the chloroplast genome of Quercus section Cyclobalanopsis showed a strong geographical structure. Therefore, the phylogeny of cpDNA can only reflect the maternal evolutionary history (Yan et al., 2018). Finally, Quercus section Cyclobalanopsis did not form a monophyletic group in our analysis. Species of the section Cyclobalanopsis were mixed with those from two other sections of the subgenus Cerris, which was consistent with the results of previous studies (Yan et al., 2018). In a basic morphological study, Quercus section Cyclobalanopsis was closely related to the section Cerris and also provided evidence for their nesting (Deng et al., 2014). The non‐monophyly of the chloroplast genomes and the close relationships among the sections may reflect ancient introgression or incomplete lineage sorting of the chloroplast genome in the ancestor lineages of subgenus Cerris (Yan et al., 2018; Yan, Liu, Li et al., 2019). In summary, the phylogenetic relationship of Quercus is complex and needs further analysis from multiple perspectives.

4.4 Adaptive evolution and mechanisms of chloroplast genomes of Quercus section Cyclobalanopsis

It is generally believed that PCGs in the organelle genomes of most plants have a lower substitution rate than those in the nuclear genomes, resulting in purifying selection (Xiang et al., 2023). The substitution rate of photosystem functional genes is the smallest in the chloroplast genome (Figure 10) (Zhang et al., 2024), suggesting that it is under stronger functional constraints (Wu et al., 2017; Yan, Liu, Yuan et al., 2019), whereas ribosomal protein genes and other functional genes have large substitution rates, possibly because they have duplicated genes in IR regions that may lead to evolution at a faster rate.

However, individual genes in the chloroplast genomes still undergo positive selection, particularly those involved in photosynthesis and other metabolic pathways (Silva et al., 2023). Species of Quercus section Cyclobalanopsis are widely distributed in the tropical and subtropical regions of East Asia (Huang et al., 1999). Over a long evolutionary period, they have adapted to diverse ecological environments and occupied different ecological niches (Jin et al., 2023). At the chloroplast genome level, nine PCGs with positive selection sites were identified in the Quercus section Cyclobalanopsis (Table 2). Most of these genes are involved in photosynthesis and other metabolic pathways, validating the previous statement. Four of these were NADH dehydrogenase genes, ndhA, ndhD, ndhF, and ndhH, which are involved in photosynthesis. The NADH complex in higher plants not only promotes ATP synthesis but also participates in chloroplast protection against photooxidative stress (Li et al., 2004; Martin et al., 1996). All species of Quercus section Cyclobalanopsis are evergreen trees with a large growth altitude span, therefore, adaptation to high ultraviolet radiation intensity may be one of the important genetic bases for the adaptive evolution of the chloroplast genome. The accD codes for a key enzyme involved in fatty acid biosynthesis (Slabas & Fawcett, 1992) and has undergone adaptive evolution, promoting the growth and development of Quercus section Cyclobalanopsis. Adaptive evolution of this gene has also been reported in other genera (Lee et al., 2007; Zhu et al., 2016). The rbcL has also been subject to positive selection in many species (Azarin et al., 2021; Kapralov & Filatov, 2007). The ycf1 and ycf2 are the two largest hypothetical chloroplast reading frame genes, but their functions are still unknown and their evolutionary significance remains to be studied.

5 CONCLUSION

In this study, we newly assembled and annotated the complete chloroplast genomes of 36 species of Quercus section Cyclobalanopsis. Combined with the existing complete chloroplast genomes of 14 other species, we carried out genome structure analysis, comparative genome analysis, phylogenetic analysis, ancestral aera reconstruction analysis, and selection pressure analysis for these 50 species of Quercus section Cyclobalanopsis. These results explored the diversity structure, phylogenetic relationship, evolutionary history, and adaptative evolutionary mechanism of the chloroplast genome of Quercus section Cyclobalanopsis. In conclusion, this study provides a wealth of valuable information that will help us to harvest deeper insights into the evolutionary history, phylogeny, and taxonomy of Quercus.

AUTHOR CONTRIBUTIONS

Yu Li: Data curation (equal); formal analysis (equal); methodology (equal); software (equal); validation (equal); visualization (equal); writing – original draft (equal). Si‐Si Zheng: Data curation (equal); resources (equal); visualization (equal); writing – review and editing (equal). Tian‐Rui Wang: Formal analysis (equal); software (equal); writing – review and editing (equal). Mei‐Hua Liu: Investigation (equal); supervision (equal); writing – review and editing (equal). Gregor Kozlowski: Project administration (equal); writing – review and editing (equal). Li‐Ta Yi: Conceptualization (equal); supervision (equal); writing – review and editing (equal). Yi‐Gang Song: Conceptualization (equal); funding acquisition (equal); project administration (equal); resources (equal); supervision (equal); writing – review and editing (equal).

FUNDING INFORMATION

This work was supported by the National Natural Science Foundation of China, grant number 31901217, and the Special Fund for Scientific Research of Shanghai Landscaping & City Appearance Administrative Bureau, grant number G192422 and G222404.

CONFLICT OF INTEREST STATEMENT

The authors declare that they have no conflict of interest. No approval/permission is required for the sample collection in accordance with relevant regulations and standards.

Supporting information

Figure S1.–S5.

Table S1.–S10.

ACKNOWLEDGEMENTS

We want to thank Wuhan Benagen Technology Co., Ltd. for their help with sequencing.

DATA AVAILABILITY STATEMENT

The data supporting the findings of this study are openly available in the NCBI GenBank repository at https://www.ncbi.nlm.nih.gov, reference number OQ603420–OQ60345.
==== Refs
REFERENCES

Amiryousefi, A. , Hyvonen, J. , & Poczai, P. (2018). IRscope: An online program to visualize the junction sites of chloroplast genomes. Bioinformatics, 34 (17 ), 3030–3031. 10.1093/bioinformatics/bty220 29659705
Azarin, K. , Usatov, A. , Makarenko, M. , Khachumov, V. , & Gavrilova, V. (2021). Comparative analysis of chloroplast genomes of seven perennial helianthus species. Gene, 774 , 145418. 10.1016/j.gene.2021.145418 33444687
Batzoglou, S. , Berger, B. , Mesirov, J. , & Lander, E. S. (1999). Sequencing a genome by walking with clone‐end sequences: A mathematical analysis. Genome Research, 9 (12 ), 1163–1174. 10.1101/gr.9.12.1163 10613838
Beier, S. , Thiel, T. , Münch, T. , Scholz, U. , & Mascher, M. (2017). MISA‐web: A web server for microsatellite prediction. Bioinformatics, 33 (16 ), 2583–2585. 10.1093/bioinformatics/btx198 28398459
Benson, G. (1999). Tandem repeats finder: A program to analyze DNA sequences. Nucleic Acids Research, 27 (2 ), 573–580. 10.1093/nar/27.2.573 9862982
Brudno, M. , Malde, S. , Poliakov, A. , Do, C. B. , Couronne, O. , Dubchak, I. , & Batzoglou, S. (2003). Glocal alignment: Finding rearrangements during alignment. Bioinformatics, 19 (Supplement 1 ), i54–i62. 10.1093/bioinformatics/btg1005 12855437
Camus, A. (1936–1938). Les Chênes. Monographie du genre Quercus. Tome I. Genre Quercus, sous‐genre Cyclobalanopsis, sous‐genre Euquercus (sections Cerris et Mesobalanus). Texte.
Chen, J. P. , Xie, D. F. , He, X. J. , Yang, Y. , & Li, X. F. (2022). Comparative analysis of the complete chloroplast genomes in Allium section Bromatorrhiza species (Amaryllidaceae): Phylogenetic relationship and adaptive evolution. Genes (Basel), 13 (7 ), 1279. 10.3390/genes13071279 35886061
Chen, T. T. , Chen, X. , Zhang, S. S. , Zhu, J. W. , Tang, B. X. , Wang, A. K. , Dong, Z. W. , Yu, C. X. , Sun, Y. L. , Chi, L. J. , Chen, H. X. , Zhai, S. , Sun, Y. B. , Lan, L. , Zhang, X. , Xiao, J. F. , Bao, Y. M. , Wang, Y. Q. , Zhang, Z. , & Zhao, W. M. (2021). The genome sequence archive family: Toward explosive data growth and diverse data types. Genomics, Proteomics & Bioinformatics, 19 (4 ), 578–583. 10.1016/j.gpb.2021.08.001
Chen, X. L. , & Zhang, X. M. (2023). Comparison of chloroplast genomes and phylogenetic analysis of four species in Quercus section Cyclobalanopsis . Scientific Reports, 13 , 18731. 10.1038/s41598-023-45421-8 37907468
Chen, Y. , Chen, Y. , Shi, C. , Huang, Z. , Zhang, Y. , Li, S. , Li, Y. , Ye, J. , Yu, C. , & Li, Z. (2018). SOAPnuke: A MapReduce acceleration‐supported software for integrated quality control and preprocessing of high‐throughput sequencing data. GigaScience, 7 (1 ), 1–6. 10.1093/gigascience/gix120
Cho, W. B. , Han, E. K. , Choi, I. S. , Kwak, W. , Kim, J. H. , Kim, B. Y. , & Lee, J. H. (2021). The complete plastid genome sequence of Quercus acuta (Fagaceae), an evergreen broad‐leaved oak endemic to East Asia. Mitochondrial DNA Part B Resources, 6 (2 ), 320–322. 10.1080/23802359.2020.1866449 33659663
CNCB‐NGDC Members and Partners . (2022). Database resources of the National Genomics Data Center, China National Center for bioinformation in 2022. Nucleic Acids Research, 50 (D1 ), D27–D38. 10.1093/nar/gkab951 34718731
Darling, A. C. E. , Mau, B. , Blattner, F. R. , & Perna, N. T. (2004). Mauve: Multiple alignment of conserved genomic sequence with rearrangements. Genome Research, 14 (7 ), 1394–1403. 10.1101/gr.2289704 15231754
Deng, M. , Hipp, A. , Song, Y. G. , Li, Q. S. , Coombes, A. , & Cotton, A. (2014). Leaf epidermal features of Quercus subgenus Cyclobalanopsis (Fagaceae) and their systematic significance. Botanical Journal of the Linnean Society, 176 (2 ), 224–259. 10.1111/boj.12207
Deng, M. , Jiang, X. L. , Hipp, A. L. , Manos, P. S. , & Hahn, M. (2018). Phylogeny and biogeography of east Asian evergreen oaks (Quercus section Cyclobalanopsis; Fagaceae): Insights into the Cenozoic history of evergreen broad‐leaved forests in subtropical Asia. Molecular Phylogenetics and Evolution, 119 , 170–181. 10.1016/j.ympev.2017.11.003 29175095
Denk, T. , & Grimm, G. W. (2010). The oaks of western Eurasia: Traditional classifications and evidence from two nuclear markers. Taxon, 59 (2 ), 351–366. 10.1002/tax.592002
Denk, T. , Grimm, G. W. , Manos, P. S. , Deng, M. , & Hipp, A. L. (2017). An updated infrageneric classification of the oaks: Review of previous taxonomic schemes and synthesis of evolutionary patterns. In E. Gil‐Pelegrín , J. J. Peguero‐Pina , & D. Sancho‐Knapik (Eds.), Oaks physiological ecology. Exploring the functional diversity of genus Quercus L. (pp. 13–38). Tree Physiology.
Doyle, J. J. (1987). A rapid DNA isolation procedure for small quantities of fresh leaf tissue. Phytochemical Bulletin, 19 , 11–15.
Drouin, G. , Daoud, H. , & Xia, J. (2008). Relative rates of synonymous substitutions in the mitochondrial, chloroplast and nuclear genomes of seed plants. Molecular Biology and Evolution, 49 (3 ), 827–831. 10.1016/j.ympev.2008.09.009
Dumolin, S. , Demesure, B. , & Petit, R. J. (1995). Inheritance of chloroplast and mitochondrial genomes in pedunculate oak investigated with an efficient PCR method. Theoretical and Applied Genetics, 91 (8 ), 1253–1256. 10.1007/BF00220937 24170054
Ewing, B. , & Green, P. (1998). Base‐calling of automated sequencer traces using phred. II. Error probabilities. Genome Research, 8 (3 ), 186–194. 10.1101/gr.8.3.186 9521922
Frazer, K. A. , Pachter, L. , Poliakov, A. , Rubin, E. M. , & Dubchak, I. (2004). VISTA: Computational tools for comparative genomics. Nucleic Acids Research, 32 (Web Server issue ), W273–W279. 10.1093/nar/gkh458 15215394
Gitzendanner, M. A. , Soltis, P. S. , Yi, T. S. , Li, D. Z. , & Soltis, D. E. (2018). Plastome phylogenetics: 30 years of inferences into plant evolution. Advances in Botanical Research, 85 , 293–313. 10.1016/bs.abr.2017.11.016
Greiner, S. , Lehwark, P. , & Bock, R. (2019). OrganellarGenomeDRAW (OGDRAW) version 1.3.1: Expanded toolkit for the graphical visualization of organellar genomes. Nucleic Acids Research, 47 (W1 ), W59–W64. 10.1093/nar/gkz238 30949694
Guisinger, M. M. , Kuehl, J. V. , Boore, J. L. , & Jansen, R. K. (2010). Extreme reconfiguration of plastid genomes in the angiosperm family Geraniaceae: Rearrangements, repeats, and codon usage. Molecular Biology and Evolution, 28 (1 ), 583–600. 10.1093/molbev/msq229 20805190
Guo, X. L. , Zheng, H. Y. , Price, M. , Zhou, S. D. , & He, X. J. (2020). Phylogeny and comparative analysis of Chinese Chamaesium species revealed by the complete plastid genome. Plants, 9 (8 ), 965. 10.3390/plants9080965 32751647
Han, Y. W. , Duan, D. , Ma, X. F. , Jia, Y. , Liu, Z. L. , Zhao, G. F. , & Li, Z. H. (2016). Efficient identification of the forest tree species in Aceraceae using DNA barcodes. Frontiers in Plant Science, 7 , 1707. 10.3389/fpls.2016.01707 27899929
He, L. , Qian, J. , Li, X. W. , Sun, Z. Y. , Xu, X. L. , & Chen, S. L. (2017). Complete chloroplast genome of medicinal plant Lonicera japonica: Genome rearrangement, intron gain and loss, and implications for phylogenetic studies. Molecules, 22 (2 ), 249. 10.3390/molecules22020249 28178222
He, W. C. , Chen, C. J. , Xiang, K. L. , Wang, J. , Zheng, P. , Tembrock, L. R. , Jin, D. M. , & Wu, Z. Q. (2021). The history and diversity of rice domestication as resolved from 1464 complete plastid genomes. Frontiers in Plant Science, 12 , 781793. 10.3389/fpls.2021.781793 34868182
Hipp, A. L. , Manos, P. S. , Hahn, M. , Avishai, M. , Bodenes, C. , Cavender‐Bares, J. , Crowl, A. A. , Deng, M. , Denk, T. , Fitz‐Gibbon, S. , Gailing, O. , Gonzalez‐Elizondo, M. S. , Gonzalez‐Rodriguez, A. , Grimm, G. W. , Jiang, X. L. , Kremer, A. , Lesur, I. , McVay, J. D. , Plomion, C. , … Valencia‐Avalos, S. (2020). Genomic landscape of the global oak phylogeny. New Phytologist, 226 (4 ), 1198–1212. 10.1111/nph.16162 31609470
Huang, C. J. , Zhang, Y. T. , Hu, Y. C. , & Jen, H. W. (1999). Fagaceae. In W. Y. Chun & C. J. Huang (Eds.), Flora of China (pp. 314–400). Science Press and Missouri Botanical Garden Press.
Jansen, R. K. , Saski, C. , Lee, S. B. , Hansen, A. K. , & Daniell, H. (2011). Complete plastid genome sequences of three Rosids (Castanea, Prunus, Theobroma): Evidence for at least two independent transfers of rpl22 to the nucleus. Molecular Biology and Evolution, 28 (1 ), 835–847. 10.1093/molbev/msq261 20935065
Jiang, H. , Tian, J. , Yang, J. X. , Dong, X. , Zhong, Z. X. , Mwahala, G. , Zhang, C. F. , Hu, G. W. , & Wang, Q. F. (2022). Comparative and phylogenetic analyses of six Kenya Polystachya (Orchidaceae) species based on the complete chloroplast genome sequences. BMC Plant Biology, 22 , 177. 10.1186/s12870-022-03529-5 35387599
Jin, D. M. , Yuan, Q. , Dai, X. L. , Kozlowski, G. , & Song, Y. G. (2023). Enhanced precipitation drove the evolution of subtropical evergreen broad‐leaved forests in eastern China since the early Miocene: Evidence from ring‐cupped oaks. Journal of Systematics and Evolution, 62 , 677–686. 10.1111/jse.13022
Jin, J. J. , Yu, W. B. , Yang, J. B. , Song, Y. , dePamphilis, C. W. , Yi, T. S. , & Li, D. Z. (2020). GetOrganelle: A fast and versatile toolkit for accurate de novo assembly of organelle genomes. Genome Biology, 21 (1 ), 241. 10.1186/s13059-020-02154-5 32912315
Ju, M. M. , Zhang, X. , Yang, Y. C. , Fan, W. B. , & Zhao, G. F. (2019). The complete chloroplast genome of a critically endangered tree species in China, Cyclobalanopsis obovatifolia (Fagaceae). Conservation Genetics Resources, 11 , 31–33. 10.1007/s12686-017-0958-z
Kalyaanamoorthy, S. , Minh, B. Q. , Wong, T. K. F. , von Haeseler, A. , & Jermiin, L. S. (2017). ModelFinder: Fast model selection for accurate phylogenetic estimates. Nature Methods, 14 (6 ), 587–589. 10.1038/nmeth.4285 28481363
Kapralov, M. V. , & Filatov, D. A. (2007). Widespread positive selection in the photosynthetic Rubisco enzyme. BMC Evolutionary Biology, 7 (1 ), 73. 10.1186/1471-2148-7-73 17498284
Katoh, K. , & Standley, D. M. (2013). MAFFT multiple sequence alignment software version 7: Improvements in performance and usability. Molecular Biology and Evolution, 30 (4 ), 772–780. 10.1093/molbev/mst010 23329690
Kearse, M. , Moir, R. , Wilson, A. , Stones‐Havas, S. , Cheung, M. , Sturrock, S. , Buxton, S. , Cooper, A. , Markowitz, S. , Duran, C. , Thierer, T. , Ashton, B. , Meintjes, P. , & Drummond, A. (2012). Geneious basic: An integrated and extendable desktop software platform for the organization and analysis of sequence data. Bioinformatics, 28 (12 ), 1647–1649. 10.1093/bioinformatics/bts199 22543367
Kleine, T. , Maier, U. G. , & Leister, D. (2009). DNA transfer from organelles to the nucleus: The idiosyncratic genetics of endosymbiosis. Annual Review of Plant Biology, 60 (1 ), 115–138. 10.1146/annurev.arplant.043008.092119
Kumar, S. , Stecher, G. , Li, M. , Knyaz, C. , & Tamura, K. (2018). MEGA X: Molecular evolutionary genetics analysis across computing platforms. Molecular Biology and Evolution, 35 (6 ), 1547–1549. 10.1093/molbev/msy096 29722887
Kurtz, S. , Choudhuri, J. V. , Ohlebusch, E. , Schleiermacher, C. , Stoye, J. , & Giegerich, R. (2001). REPuter: The manifold applications of repeat analysis on a genomic scale. Nucleic Acids Research, 29 (22 ), 4633–4642. 10.1093/nar/29.22.4633 11713313
Lee, H. L. , Jansen, R. K. , Chumley, T. W. , & Kim, K. J. (2007). Gene relocations within chloroplast genomes of Jasminum and Menodora (Oleaceae) are due to multiple, overlapping inversions. Molecular Biology and Evolution, 24 (5 ), 1161–1180. 10.1093/molbev/msm036 17329229
Li, H. T. , Yi, T. S. , Gao, L. M. , Ma, P. F. , Zhang, T. , Yang, J. B. , Gitzendanner, M. A. , Fritsch, P. W. , Cai, J. , Luo, Y. , Wang, H. , van der Bank, M. , Zhang, S. D. , Wang, Q. F. , Wang, J. , Zhang, Z. R. , Fu, C. N. , Yang, J. , Hollingsworth, P. M. , … Li, D. Z. (2019). Origin of angiosperms and the puzzle of the Jurassic gap. Nature Plants, 5 (5 ), 461–470. 10.1038/s41477-019-0421-0 31061536
Li, J. , Xie, D. F. , Guo, X. L. , Zheng, Z. Y. , He, X. J. , & Zhou, S. D. (2020). Comparative analysis of the complete plastid genome of five Bupleurum species and new insights into DNA barcoding and phylogenetic relationship. Plants, 9 (4 ), 543. 10.3390/plants9040543 32331381
Li, T. T. , Luo, C. S. , Mou, H. L. , Jiang, X. L. , & Xu, G. B. (2021). The complete chloroplast genome sequence of Quercus schottkyana, and comparative analysis with related species. Mitochondrial DNA Part B Resources, 6 (9 ), 2607–2609. 10.1080/23802359.2021.1961630 34395892
Li, X. , Li, Y. F. , Sylvester, S. P. , Zang, M. Y. , Ei‐Kassaby, Y. A. , & Fang, Y. M. (2021). Evolutionary patterns of nucleotide substitution rates in plastid genomes of Quercus . Ecology and Evolution, 11 (19 ), 13401–13414. 10.1002/ece3.8063 34646478
Li, X. , Li, Y. F. , Zang, M. Y. , Li, M. Z. , & Fang, Y. M. (2018). Complete chloroplast genome sequence and phylogenetic analysis of Quercus acutissima . International Journal of Molecular Sciences, 19 (8 ), 2443. 10.3390/ijms19082443 30126202
Li, X. G. , Wei, D. , Meng, Q. W. , Zou, Q. , & Zhao, S. J. (2004). The function of chloroplastic NAD(P)H dehydrogenase in tobacco during chilling stress under low irradiance. Plant and Cell Physiology, 45 (1 ), 103–108. 10.1093/pcp/pch011 14749491
Li, Y. , Wang, L. , & Fang, Y. (2019). The complete chloroplast genome sequence of Quercus myrsinifolia (Fagaceae). Mitochondrial DNA Part B Resources, 4 (2 ), 3139–3140. 10.1080/23802359.2019.1667276 33365888
Li, Y. , Wang, L. , Liu, Q. L. , & Fang, Y. M. (2020). The complete plastid genome sequence of Quercus ciliaris (Fagaceae). Mitochondrial DNA Part B Resources, 5 (2 ), 1954–1955. 10.1080/23802359.2020.1756955
Li, Y. , Wang, L. , Zhao, Y. , & Fang, Y. (2020). The complete chloroplast genome sequence of Quercus stewardiana (Fagaceae). Mitochondrial DNA Part B Resources, 5 (2 ), 1958–1959. 10.1080/23802359.2020.1756961
Li, Y. , Wang, T. R. , Kozlowski, G. , Liu, M. H. , Yi, L. T. , & Song, Y. G. (2022). Complete chloroplast genome of an endangered species Quercus litseoides, and its comparative, evolutionary, and phylogenetic study with other Quercus section Cyclobalanopsis species. Genes, 13 , 1184. 10.3390/genes13071184 35885967
Li, Y. F. , Sylvester, S. P. , Li, M. , Zhang, C. , Li, X. , Duan, Y. F. , & Wang, X. R. (2019). The complete plastid genome of Magnolia zenii and genetic comparison to Magnoliaceae species. Molecules, 24 (2 ), 261–278. 10.3390/molecules24020261 30641990
Liang, C. L. , Wang, L. , Lei, J. , Duan, B. Z. , Ma, W. S. , Xiao, S. M. , Qi, H. J. , Wang, Z. , Liu, Y. Q. , Shen, X. F. , Guo, S. , Hu, H. Y. , Xu, J. , & Chen, S. L. (2019). A comparative analysis of the chloroplast genomes of four Salvia medicinal plants. Engineering, 5 (5 ), 907–915. 10.1016/j.eng.2019.01.017
Liu, D. , Guo, H. , Zhu, J. , Qu, K. , Chen, Y. , Guo, Y. , Ding, P. , Yang, H. , Xu, T. , Jing, Q. , Han, S. , Li, W. , & Tong, B. (2022). Complex physical structure of complete mitochondrial genome of Quercus acutissima (Fagaceae): A significant energy plant. Genes (Basel), 13 (8 ), 1321. 10.3390/genes13081321 35893058
Liu, S. , Ni, Y. , Li, J. , Zhang, X. , Yang, H. , Chen, H. , & Liu, C. (2023). CPGView: A package for visualizing detailed chloroplast genome structures. Molecular Ecology Resources, 23 (3 ), 694–704. 10.1111/1755-0998.13729 36587992
Liu, X. , Chang, E. M. , Liu, J. F. , & Jiang, Z. P. (2021). Comparative analysis of the complete chloroplast genomes of six white oaks with high ecological amplitude in China. Journal of Forestry Research, 32 (5 ), 2203–2218. 10.1007/s11676-020-01288-3
Luo, Y. , & Zhou, Z. K. (2001). Phytogeography of Quercus subg. Cyclobalanopsis. (in Chinese.). Chinese Acta Botanica Yunnaica, 23 (1 ), 1–16.
Manos, P. S. , Doyle, J. J. , & Nixon, K. C. (1999). Phylogeny, biogeography, and processes of molecular differentiation in Quercus subgenus Quercus (Fagaceae). Molecular Phylogenetics and Evolution, 12 (3 ), 333–349. 10.1006/mpev.1999.0614 10413627
Marechal, A. , & Brisson, N. (2010). Recombination and the maintenance of plant organelle genome stability. New Phytologist, 186 (2 ), 299–317. 10.1111/j.1469-8137.2010.03195x 20180912
Martin, M. , Casano, L. M. , & Sabater, B. (1996). Identification of the product of ndhA gene as a thylakoid protein synthesized in response to photooxidative treatment. Plant and Cell Physiology, 37 (3 ), 293–298. 10.1093/oxfordjournals.pcp.a028945 8673340
McVay, J. D. , Hauser, D. , Hipp, A. L. , & Manos, P. S. (2017). Phylogenomics reveals a complex evolutionary history of lobed‐leaf white oaks in western North America. Genome, 60 (9 ), 733–742. 10.1139/gen-2016-0206 28727934
Menitsky, Y. L. (2005). Oaks of Asia (translated from the Russian original of 1984). Science Publishers.
Minh, B. Q. , Schmidt, H. A. , Chernomor, O. , Schrempf, D. , Woodhams, M. D. , von Haeseler, A. , & Lanfear, R. (2020). IQ‐TREE 2: New models and efficient methods for phylogenetic inference in the genomic era. Molecular Biology and Evolution, 37 (5 ), 1530–1534. 10.1101/849372 32011700
Morton, B. R. (1997). The influence of neighboring base composition on substitutions in plant chloroplast coding sequences. Molecular Biology and Evolution, 14 (2 ), 189–194. 10.1093/oxfordjournals.molbev.a025752
Morton, B. R. (1998). Selection on the codon bias of chloroplast and cyanelle genes in different plant and algal lineages. Journal of Molecular Evolution, 46 (4 ), 449–459. 10.1007/pl00006325 9541540
Nixon, K. C. (1993). Infrageneric classification of Quercus (Fagaceae) and typification of sectional names. Annales des Sciences Forestières, 36 (Supplement ), 25s–34s. 10.1051/forest:19930701
Ohyama, M. , Baba, K. , & Itoh, T. (2001). Wood identification of Japanese Cyclobalanopsis species (Fagaceae) based on DNA polymorphism of the intergenic spacer between trnT and trnL exon. Journal of Wood Science, 47 , 81–86. 10.1007/BF00780554
Ørsted, A. S. (1871). Bidrag til Kundskab om Egefamilien. Kongelige Danske Videnskabernes Selskabs Skrifter Naturvidenskabelig Og Mathematisk Afdeling, 6 , 331–538.
Palmer, J. D. (1985). Comparative organization of chloroplast genomes. Annual Review of Genetics, 19 (1 ), 325–354. 10.1146/annurev.ge.19.120185.001545
Pang, X. B. , Liu, H. S. , Wu, S. R. , Yuan, Y. C. , Li, H. J. , Dong, G. S. , Liu, Z. H. , An, C. Z. , Su, Z. H. , & Li, B. (2019). Species identification of oaks (Quercus L., Fagaceae) from gene to genome. International Journal of Molecular Sciences, 20 (23 ), 5940. 10.3390/ijms20235940 31779118
Petit, R. J. , Csaikl, U. M. , Bordács, S. , Burg, K. , Coart, E. , Cottrell, J. , van Dam, B. , Deans, J. D. , Dumolin‐Lapègue, S. , Fineschi, S. , Finkeldey, R. , Gillies, A. , Glaz, I. , Goicoechea, P. G. , Jensen, J. S. , König, A. O. , Lowe, A. J. , Madsen, S. F. , Mátyás, G. , … Kremer, A. (2002). Chloroplast DNA variation in European white oaks: Phylogeography and patterns of diversity based on data from over 2600 populations. Forest Ecology and Management, 156 , 5–26. 10.1016/S0378-1127(01)00645-4
Pham, K. K. , Hipp, A. L. , Manos, P. S. , & Cronn, R. C. (2017). A time and a place for everything: Phylogenetic history and geography as joint predictors of oak plastome phylogeny. Genome, 60 (9 ), 720–732. 10.1139/gen-2016-0191 28445658
Qin, Z. , Cai, Z. Q. , Xia, G. M. , & Wang, M. C. (2013). Synonymous codon usage bias is correlative to intron number and shows disequilibrium among exons in plants. BMC Genomics, 14 (1 ), 56. 10.1186/1471-2164-14-56 23350908
Rozas, J. , Ferrer‐Mata, A. , Sanchez‐DelBarrio, J. C. , Guirao‐Rico, S. , Librado, P. , Ramos‐Onsins, S. E. , & Sanchez‐Gracia, A. (2017). DnaSP 6: DNA sequence polymorphism analysis of large data sets. Molecular Biology and Evolution, 34 (12 ), 3299–3302. 10.1093/molbev/msx248 29029172
Schwarz, O. (1936). Entwurf zu einem natürlichen System der Cupuliferen und der Gattung Quercus L. Notizblatt des Botanischen Gartens und Museums zu Berlin‐Dahlem, 116 , 1–22.
Shaw, J. , Lickey, E. B. , Schilling, E. E. , & Small, R. L. (2007). Comparison of whole chloroplast genome sequences to choose noncoding regions for phylogenetic studies in angiosperms: The tortoise and the hare III. American Journal of Botany, 94 (3 ), 275–288. 10.3732/ajb.94.3.275 21636401
Shi, S. L. , Liu, Y. Q. , Xia, R. X. , & Qin, L. (2022). Comprehensive analysis of codon usage in Quercus chloroplast genome and focus on psbA gene. Genes, 13 , 2156. 10.3390/genes13112156 36421830
Silva, S. R. , Miranda, V. F. O. , Michael, T. P. , Płachno, B. J. , Matos, R. G. , Adamec, L. , Pond, S. L. K. , Lucaci, A. G. , Pinheiro, D. G. , & Varani, A. M. (2023). The phylogenomics and evolutionary dynamics of the organellar genomes in carnivorous Utricularia and Genlisea species (Lentibulariaceae). Molecular Phylogenetics and Evolution, 181 , 107711. 10.1016/j.ympev.2023.107711 36693533
Simeone, M. C. , Grimm, G. W. , Papini, A. , Vessella, F. , Cardoni, S. , Tordoni, E. , Piredda, R. , Franc, A. , & Denk, T. (2016). Plastome data reveal multiple geographic origins of Quercus group ilex . PeerJ, 4 , e1897. 10.7717/peerj.1897 27123376
Simeone, M. C. , Piredda, R. , Papini, A. , Vessella, F. , & Schirone, B. (2013). Application of plastid and nuclear markers to DNA barcoding of Euro‐Mediterranean oaks (Quercus, Fagaceae): Problems, prospects and phylogenetic implications. Botanical Journal of the Linnean Society, 172 (4 ), 478–499. 10.1111/boj.12059
Slabas, A. R. , & Fawcett, T. (1992). The biochemistry and molecular biology of plant lipid biosynthesis. Plant Molecular Biology, 19 (1 ), 169–191. 10.1007/BF00015613 1600168
Song, Y. G. , Petitpierre, B. , Deng, M. , Wu, J. P. , & Kozlowski, G. (2019). Predicting climate change impacts on the threatened Quercus arbutifolia in montane cloud forests in southern China and Vietnam: Conservation implications. Forest Ecology and Management, 444 , 269–279. 10.1016/j.foreco.2019.04.028
Stegemann, S. , Hartmann, S. , Ruf, S. , & Bock, R. (2003). High‐frequency gene transfer from the chloroplast genome to the nucleus. Proceedings of the National Academy of Sciences of the United States of America, 100 (15 ), 8828–8833. 10.1073/pnas.1430924100 12817081
Su, H. , Yang, Y. , Ju, M. , Li, H. , & Zhao, G. F. (2019). Characterization of the complete plastid genome of Quercus sichourensis . Conservation Genetics Resources, 11 , 129–131. 10.1007/s12686-017-0977-9
Thompson, J. D. , Gibson, T. J. , & Higgins, D. G. (2002). Multiple sequence alignment using ClustalW and ClustalX. Current Protocols in Bioinformatics, 2 (1 ), Unit 2.3. 10.1002/0471250953.bi0203s00
Tillich, M. , Lehwark, P. , Pellizzer, T. , Ulbricht‐Jones, E. S. , Fischer, A. , Bock, R. , & Greiner, S. (2017). GeSeq‐versatile and accurate annotation of organelle genomes. Nucleic Acids Research, 45 (W1 ), W6–W11. 10.1093/nar/gkx391 28486635
Timme, R. E. , Kuehl, J. V. , Boore, J. L. , & Jansen, R. K. (2007). A comparative analysis of the Lactuca and Helianthus (Asteraceae) plastid genomes: Identification of divergent regions and categorization of shared repeats. American Journal of Botany, 94 (3 ), 302–312. 10.3732/ajb.94.3.302 21636403
Tu, X. D. , Liu, D. K. , Xu, S. W. , Zhou, C. Y. , Gao, X. Y. , Zeng, M. Y. , Zhang, S. , Chen, J. L. , Ma, L. , Zhou, Z. , Huang, M. Z. , Chen, S. P. , Liu, Z. J. , Lan, S. R. , & Li, M. H. (2021). Plastid phylogenomics improves resolution of phylogenetic relationship in the Cheirostylis and Goodyera clades of Goodyerinae (Orchidoideae, Orchidaceae). Molecular Phylogenetics and Evolution, 164 , 107269. 10.1016/j.ympev.2021.107269 34324956
Wang, L. Y. , Xing, H. X. , Yuan, Y. C. , Wang, X. L. , Saeed, M. , Tao, J. C. , Feng, W. , Zhang, G. H. , Song, X. L. , & Sun, X. Z. (2018). Genome‐wide analysis of codon usage bias in four sequenced cotton species. PLoS One, 13 (3 ), e0194372. 10.1371/journal.pone.0194372 29584741
Wang, T. R. , Ning, X. , Zheng, S. S. , Li, Y. , Lu, Z. J. , Meng, H. H. , Ge, B. J. , Kozlowski, G. , Yan, M. X. , & Song, Y. G. (2024). Genomic insights into ecological adaptation of oaks revealed by phylogenomic analysis of multiple species. Plant Diversity. 10.1016/j.pld.2024.07.008
Wang, T. R. , Wang, Z. W. , Song, Y. G. , & Kozlowski, G. (2021). The complete chloroplast genome sequence of Quercus ningangensis and its phylogenetic implications. Plant and Fungal Systematics, 66 (2 ), 155–165. 10.35535/pfsyst-2021-0014
Wen, F. , Wu, X. Z. , Li, T. J. , Jia, M. L. , Liu, X. S. , & Liao, L. (2021). The complete chloroplast genome of Stauntonia chinensis and compared analysis revealed adaptive evolution of subfamily Lardizabaloideae species in China. BMC Genomics, 22 , 161. 10.1186/s12864-021-07484-7 33676415
Weng, M. L. , Blazier, J. C. , Govindu, M. , & Jansen, R. K. (2014). Reconstruction of the ancestral plastid genome in geraniaceae reveals a correlation between genome rearrangements, repeats, and nucleotide substitution rates. Molecular Biology and Evolution, 31 (3 ), 645–659. 10.1093/molbev/mst257 24336877
Wright, F. (1990). The ‘effective number of codons’ used in a gene. Gene, 87 (1 ), 23–29. 10.1016/0378-1119(90)90491-9 2110097
Wu, C. S. , Wang, T. J. , Wu, C. W. , Wang, Y. N. , & Chaw, S. M. (2017). Plastome evolution in the sole hemiparasitic genus Laurel dodder (Cassytha) and insights into the plastid phylogenomics of Lauraceae. Genome Biology and Evolution, 9 (10 ), 2604–2614. 10.1093/gbe/evx177 28985306
Wu, C. Y. , Sun, H. , Zhou, Z. K. , Li, D. Z. , & Peng, H. (2011). Floristics of seed plants from China. Science Press.
Wu, Y. , Liu, F. , Yang, D. G. , Li, W. , Zhou, X. J. , Pei, X. Y. , Liu, Y. G. , He, K. L. , Zhang, W. S. , Ren, Z. Y. , Zhou, K. H. , Ma, X. F. , & Li, Z. H. (2018). Comparative chloroplast genomics of Gossypium species: Insights into repeat sequence variations and phylogeny. Frontiers in Plant Science, 9 , 376. 10.3389/fpls.2018.00376 29619041
Xiang, K. L. , Wu, S. D. , Lian, L. , He, W. C. , Peng, D. , Peng, H. W. , Zhang, X. N. , Li, H. L. , Xue, J. Y. , Shan, H. Y. , Xu, G. X. , Liu, Y. , Wu, Z. Q. , & Wang, W. (2023). Genomic data and ecological niche modeling reveal an unusually slow rate of molecular evolution in the cretaceous Eupteleaceae. Science China Life Sciences, 67 , 803–816. 10.1007/s11427-023-2448-x 38087029
Xie, D. F. , Tan, J. B. , Yu, Y. , Gui, L. J. , Su, D. M. , Zhou, S. D. , & He, X. J. (2020). Insights into phylogeny, age and evolution of Allium (Amaryllidaceae) based on the whole plastome sequences. Annals of Botany, 125 (7 ), 1039–1055. 10.1093/aob/mcaa024 32239179
Yan, M. , Liu, R. , Li, Y. , Hipp, A. L. , Deng, M. , & Xiong, Y. (2019). Ancient events and climate adaptive capacity shaped distinct chloroplast genetic structure in the oak lineages. BMC Evolutionary Biology, 19 (1 ), 202. 10.1186/s12862-019-1523-z 31684859
Yan, M. , Zhao, X. Q. , Zhou, J. Q. , Huo, Y. , Ding, Y. , & Yuan, Z. H. (2019). The complete chloroplast genomes of Punica granatum and a comparison with other species in Lythraceae. International Journal of Molecular Sciences, 20 (12 ), 2886. 10.3390/ijms20122886 31200508
Yan, M. X. , Xiong, Y. S. , Liu, R. B. , Deng, M. , & Song, J. J. (2018). The application and limitation of universal chloroplast markers in discriminating east Asian evergreen oaks. Frontiers in Plant Science, 9 , 569. 10.3389/fpls.2018.00569 29868047
Yan, X. K. , Liu, T. J. , Yuan, X. , Xu, Y. , Yan, H. F. , & Hao, G. (2019). Chloroplast genomes and comparative analyses among thirteen taxa within Myrsinaceae s.str. Clade (Myrsinoideae, Primulaceae). International Journal of Molecular Sciences, 20 (18 ), 4534. 10.3390/ijms20184534 31540236
Yang, J. , Vázquez, L. , Chen, X. , Li, H. , Zhang, H. , Liu, Z. , & Zhao, G. (2017). Development of chloroplast and nuclear DNA markers for Chinese oaks (Quercus subgenus Quercus) and assessment of their utility as DNA barcodes. Frontiers in Plant Science, 8 , 816. 10.3389/fpls.2017.00816 28579999
Yang, Y. C. , Zhou, T. , Duan, D. , Yang, J. , Feng, L. , & Zhao, G. F. (2016). Comparative analysis of the complete chloroplast genomes of five Quercus species. Frontiers in Plant Science, 7 , 959. 10.3389/fpls.2016.00959 27446185
Yang, Y. C. , Zhou, T. , Qian, Z. Q. , & Zhao, G. F. (2021). Phylogenetic relationships in Chinese oaks (Fagaceae, Quercus): Evidence from plastid genome using low‐coverage whole genome sequencing. Genomics, 113 (3 ), 1438–1447. 10.1016/j.ygeno.2021.03.013 33744343
Yang, Y. C. , Zhu, J. , Feng, L. , Zhou, T. , Bai, G. Q. , Yang, J. , & Zhao, G. F. (2018). Plastid genome comparative and phylogenetic analyses of the key genera in Fagaceae: Highlighting the effect of codon composition bias in phylogenetic inference. Frontiers in Plant Science, 9 , 82. 10.3389/fpls.2018.00082 29449857
Yang, Y. Y. , Qu, X. J. , Zhang, R. , Stull, G. W. , & Yi, T. S. (2021). Plastid phylogenomic analyses of Fagales reveal signatures of conflict and ancient chloroplast capture. Molecular Phylogenetics and Evolution, 163 , 107232. 10.1016/j.ympev.2021.107232 34129935
Yang, Z. (2007). PAML 4: Phylogenetic analysis by maximum likelihood. Molecular Biology and Evolution, 24 (8 ), 1586–1591. 10.1093/molbev/msm088 17483113
Yao, X. H. , Tang, P. , Li, Z. Z. , Li, D. W. , Liu, Y. F. , & Huang, H. W. (2015). The first complete chloroplast genome sequences in Actinidiaceae: Genome structure and comparative analysis. PLoS One, 10 (6 ), e0129347. 10.1371/journal.pone.0129347 26046631
Yu, Y. , Blair, C. , & He, X. J. (2020). RASP 4: Ancestral state reconstruction tool for multiple genes and characters. Molecular Biology and Evolution, 37 (2 ), 604–606. 10.1093/molbev/msz257 31670774
Yu, Y. , Harris, A. J. , Blair, C. , & He, X. J. (2015). RASP (reconstruct ancestral state in phylogenies): A tool for historical biogeography. Molecular Phylogenetics and Evolution, 87 , 46–49.25819445
Zhai, W. , Duan, X. S. , Zhang, R. , Guo, C. C. , Li, L. , Xu, G. X. , Shan, H. Y. , Kong, H. Z. , & Ren, Y. (2019). Chloroplast genomic data provide new and robust insights into the phylogeny and evolution of the Ranunculaceae. Molecular Phylogenetics and Evolution, 135 , 12–21. 10.1016/j.ympev.2019.02.024 30826488
Zhang, R. S. , Yang, J. , Hu, H. L. , Xia, R. X. , Li, Y. P. , Su, J. F. , Li, Q. , Liu, Y. Q. , & Qin, L. (2020). A high level of chloroplast genome sequence variability in the Sawtooth oak Quercus acutissima . International Journal of Biological Macromolecules, 152 , 340–348. 10.1016/j.ijbiomac.2020.02.201 32109476
Zhang, S. Y. , Yan, H. F. , Wei, L. , Liu, T. J. , Chen, L. , Wu, X. , & Zhang, Q. L. (2024). Plastid genome and its phylogenetic implications of Asiatic Spiraea (Rosaceae). BMC Plant Biology, 24 , 23. 10.1186/s12870-023-04697-8 38166728
Zhu, A. D. , Guo, W. H. , Gupta, S. , Fan, W. S. , & Mower, J. P. (2016). Evolutionary dynamics of the plastid inverted repeat: The effects of expansion, contraction, and loss on substitution rates. New Phytologist, 209 (4 ), 1747–1756. 10.1111/nph.13743 26574731
Zong, D. , Zhou, A. P. , Zhang, Y. , Zou, X. L. , Li, D. , Duan, A. A. , & He, C. Z. (2019). Characterization of the complete chloroplast genomes of five Populus species from the western sichuan plateau, southwest China: Comparative and phylogenetic analyses. PeerJ, 7 (1 ), e6386. 10.7717/peerj.6386 30809432
