Skip to main content

Analysis of complete chloroplast genome sequences and insight into the phylogenetic relationships of Ferula L.



Ferula L. is one of the largest and most taxonomically complicated genera as well as being an important medicinal plant resource in the family Apiaceae. To investigate the plastome features and phylogenetic relationships of Ferula and its neighboring genera Soranthus Ledeb., Schumannia Kuntze., and Talassia Korovin, we sequenced 14 complete plastomes of 12 species. 


The size of the 14 complete chloroplast genomes ranged from 165,607 to 167,013 base pairs (bp) encoding 132 distinct genes (87 protein-coding, 37 tRNA, and 8 rRNA genes), and showed a typical quadripartite structure with a pair of inverted repeats (IR) regions. Based on comparative analysis, we found that the 14 plastomes were similar in codon usage, repeat sequence, simple sequence repeats (SSRs), and IR borders, and had significant collinearity. Based on our phylogenetic analyses, Soranthus, Schumannia, and Talassia should be considered synonymous with Ferula. Six highly divergent regions (rps16/trnQ-UUG, trnS-UGA/psbZ, psbH/petB, ycf1/ndhF, rpl32, and ycf1) were also detected, which may represent potential molecular markers, and combined with selective pressure analysis, the weak positive selection gene ccsA may be a discriminating DNA barcode for Ferula species.


Plastids contain abundant informative sites for resolving phylogenetic relationships. Combined with previous studies, we suggest that there is still much room for improvement in the classification of Ferula. Overall, our study provides new insights into the plastome evolution, phylogeny, and taxonomy of this genus.

Peer Review reports


Ferula L., a perennial single- or multi-bearing herb in the family Apiaceae, contains approximately 170 species mainly distributed in the Mediterranean region of southern Europe, northern Africa, Iran, Afghanistan, Central Asia, Siberia, Russia, India, and Pakistan [1]. Some Ferula species can secrete aromatic resins that have the aroma of onions and garlic, and these aromatic resins have insecticidal and fatigue-reducing properties, can be used to treat stomach diseases, dyspepsia, and abdominal pain, and is a plant resource with potentially important medicinal value [2,3,4,5].

Due to the similar morphologies and wide distribution of its constituent species, Ferula is recognized as one of the most taxonomically complicated genera within the Apiaceae [5,6,7,8]. Ferula was originally divided into three sections, Euferula Boiss., Peucedanoides Boiss., and Scorodosma Bunge [9], and later into four subgenera, Scorodosma (Bunge) Boiss., Narthex (Falc.) Drude, Soranthus Ledeb., and Euryangium (Kauffm.) Drude [10]. Fifty years later, Korovin [11] systematically divided the genus into six subgenera according to fruit, inflorescence, petals, and the number of vitta in fruits, namely Scorodosma (Bunge) Drude, Merwia (B. Fedtsch.) Korovin (including the sections Saprosmia Korovin, Phacocarpa Korovin, and Discicarpa Korovin), Narthex (Falc.) Drude (including the sections Paleonarthex Korovin and Neonarthex Korovin), Euferula (Boiss.) Korovin (including the sections Phyllites Korovin and Anatriches Korovin), Peucedanoides (Boiss.) Korovin (including the sections Xeronarthex Korovin and Macrorrhiza Korovin), and Dorematoides (Rgl. et Schmalh.) Korovin [11]. This arrangement was met with both approval [12] and opposition [13,14,15]. Safina and Pimenov [8] suggested that the genus Merwia was not naturally monophyletic and should be reduced as a section. Integrating the available research, a new classification system was subsequently proposed based on the nuclear ribosomal (nr) DNA internal transcribed spacer (ITS) and three plastid regions (rps16 intron, rpoC1 intron, and rpoB-trnC) with a total of four subgenera, namely Sinoferula Spalik, Puchałka & M.Panahi, Safinia Spalik, M.Panahi & Puchałka, Ferula (including the sections Ferula and Stenocarpa Puchałka & Spalik), and Narthex (Falc.) Drude (including the sections Glaucoselinum (Schischk.) Pimenov, Macrorrhiza Korovin, Soranthus (Ledeb.) Pimenov, Peucedanoides Boiss., Pachycarpa (Korovin) Banasiak, Euryangium (Kauffm.) Pimenov, Scorodosma (Bunge) Boiss., and Merwia (B. Fedtsch.) Koso-Pol.) [16]. However, this system still has many problems, such as the incongruence between nrDNA and plastid DNA as well as the chaotic interspecific relationship within the sections Merwia, Scorodosma, and Peucedanoides.

Moreover, the relationship between Ferula and some neighboring genera has been debated frequently, especially in the cases of Soranthus Ledeb., Schumannia Kuntz., and Talassia Korovin. Soranthus was established as a monotypic genus by Ledebour [17], with S. sibiricus (Willd.) Koso-Pol. considered a combination based on F. sibirica Willd. as published in 1798. However, this taxonomic treatment was not accepted, and Soranthus was subsequently merged into Ferula by Bunge [18], Drude [10], Safina and Pimenov [14], and Piminov [1]. Of specific note, Soranthus is recognized as a separate genus in the Flora of the Soviet Union [19], the Flora of China [20], and the Flora Xinjiangensis [21]. The same situation occurs in Schumannia, which was established as a monotypic genus with the type S. turcomanica Kuntz. [22]. S. turcomanica is a later homonym of Ferula karelinii Bunge, published by Bunge in 1851. In 1947, Korovin described the replacement Schumannia karelinii (Bunge) Korovin; however, Ferula karelinii was also listed within Ferula by Bunge [18], Drude [10], Safina and Pimenov [14], Piminov [1], and Tojibaev et al. [23], but not in the Flora of the Soviet Union [19], the Flora of China [20], or the Flora Xinjiangensis [21]. Talassia renardii (Regel & Schmalh.) Korovin and T. transiliensis (Herder) Korovin, which were isolated from Peucedanum transiliensis Regel & Herder from the genus Peucedanum L. [24], were recorded in the Flora of Kazakhstan and subsequently transferred to Ferula by Pimenov [25] and admitted by Govaerts et al. [26]. However, Talassia has also been listed as an independent genus in some Chinese floras [20, 21]. In addition, some studies have suggested that Schumannia should be merged with Soranthus based on their fruit, pollen morphology, and serological investigations [27, 28]. Recentlly, some molecular phylogeny based on the relatively limited number of nrDNA and cpDNA sequences indicated that Soranthus, Schumannia, and Talassia were embedded in Ferula, but show low support values [7, 16, 29].

Chloroplasts are independent organelles in plant cells that have their own complete set of genomes and typically covalently closed circular DNA, which exists in cells as multiple copies [30]. The chloroplast genomes of higher plants have a highly conserved tetrad structure involving inverted repeat sequences (IRs) and large single-copy (LSC) and small single-copy (SSC) regions [31]. Chloroplast genomes are relatively conserved in terms of gene number and sequence in terrestrial plants [32]. The sizes of chloroplast genomes are generally within the range of 115–165 kb, and genome size variation is mainly affected by reverse repeat length variation. Additionally, chloroplast genomes usually exhibit uniparental inheritance and low nucleotide substitution rates [33]. At present, chloroplast genome sequences and nuclear genome sequences can be obtained using shallow whole genome sequencing technology. This is considered an effective means of improving the rate of species identification and has been developed as a tool for plant phylogenetic studies at different taxonomic levels [34,35,36,37,38,39,40,41,42]. For example, the complete plastomes and nrDNA sequences obtained based on shallow genome sequencing have greatly improved the species identification rate of Rhododendron, which is also difficult to classify [43]. Thus, the complete plastomes might insight into the phylogenetic relationships of Ferula and its neighboring genera.

Here, we used plastomes to infer the phylogenetic relationships between Ferula and its confused neighboring genera. Fourteen newly sequenced plastomes of Ferula (including Soranthus, Schumannia, and Talassia) were analyzed to (1) conduct comprehensive research on the Ferula chloroplast genome; (2) identify hotspot regions, microsatellite types, and comparative genomic divergence; (3) analyze the relationships between Ferula, Soranthus, Schumannia, and Talassia based on their complete chloroplast genomes; and (4) serve as a reference for subsequent phylogenomic studies of the genus Ferula.


Chloroplast genome features

The 14 complete cp genomes ranged from 165,607 to 167,013 bp. Newly sequenced Ferula chloroplast genome maps are shown in Fig. 1. All cp genomes possessed the typical quadripartite structure of angiosperms, consisting of a pair of inverted repeat regions (IRs: 31,392–31,880 bp) and a circular molecular structure (Fig. 1; Table 1). All 14 cp genomes possessed 133 distinct genes arranged in the same order, including 87 protein-coding genes, 37 tRNA genes, and eight rRNA genes. Of these, 14 protein-coding genes and eight tRNAs contained at least one intron. The genes were classified into the following four groups based on their functions: (1) 74 self-replication genes; (2) 45 photosynthesis-related genes (in Rubisco, ATP synthase, Photosystem I, cytochrome b/f complex, photosystem II, and NADH dehydrogenase groups); and 13 other genes including (3) six genes with known functions (matK, cemA, accD, ccsA, infA, and clpP) and (4) seven genes with unknown functions (ycf1(2), ycf2(2), ycf3, ycf4, and ycf15) (Table 2). The total GC content for 12 sequenced species was 37.8–38.0% (Table 1).

Fig. 1
figure 1

Chloroplast genome maps for Ferula L. Genes on the inside of the circle are transcribed clockwise and those on the outside are transcribed counterclockwise. The darker gray inner circle corresponds to the GC content, whereas the lighter gray indicates the AT content. Different colors represent different functional genes

Table 1 Newly sequenced and complete chloroplast genomes of Ferula species
Table 2 List of genes in the chloroplast genomes of the examined Ferula species

Codon usage

The RSCU values of all codons are shown in Fig. 2 in the form of a heatmap; the red values indicate higher RSCU values, and the blue values indicate lower RSCU values. For Ferula species, the most commonly used transcription initiation codon was AUG, the most commonly used termination codon was UAA, and the initiation codon AUU only existed in F. olivacea. Except for the initiation codon and termination codon, the most used transcription codon was UTA, and AGC showed the lowest RSCU values; the most abundant amino acid (AA) was leucine, while cysteine was the lowest frequency AA. Except for tryptophan, all AAs had more than one synonymous codon, and three AAs (leucine, serine, and arginine) had the most (six) synonymous codons. The use of one codon, UGG, showed no bias (RSCU = 1) (Table S2).

Fig. 2
figure 2

Relative synonymous codon usage (RSCU) values of all merged protein-coding genes for 14 Ferula plastomes. Color key: red values indicate higher RSCU values, and blue values indicate lower RSCU values. M = initiation codon, * = termination codon, I = l-isoleucine codon, and V = Valine codon

Repeat structure analysis

Forward, palindromic, reverse, and complementary repeats were detected in 14 Ferula plastomes. Except for IR repeats, 837 repeats were identified in total; the numbers of forward repeats (398) and palindromic repeats (421) were much higher than the complement repeats (7) and reverse repeats (11). Reverse and complementary repeats were missing in four samples (F. sibirica 1, F. kelifi, F. ovina, and F. karelinii 3). F. kelifi contained the maximum number of repeats (94), whereas F. equisetacea and F. olivacea contained the least (46) (Table S3). A total of 1,061 SSRs were identified in the 14 species, six of which did not have pentanucleotides, and hexanucleotides were only found in F. olivacea. Additionally, mononucleotides were most frequent followed by dinucleotides, tetranucleotides, trinucleotides, pentanucleotides, and hexanucleotides. F. transiliensis-1 contained the highest number of SSRs (82), whereas F. oopoda contained the least (69). Poly (A/T) SSRs were typically most common, while poly (C/G) repeats were extremely rare (Table S4).

Comparisons of border and sequence identity

Single-copy and inverted repeat borders were examined; F. kelifi and F. equisetacea harbored the longest (31,880 bp) and shortest (31,392 bp) IR regions, respectively. Among all 14 Ferula species, rps19 is embedded in the LSC/IRb junction region and only 81 bp with the IRb overlap; ycf1 spans SSC/IRa and occupies a long section in both regions; and trnH occurs in the LSC region and is only 5 bp away from IRa, except for F. sibirica 3 (11 bp). The variety of IRb/SSC is relatively high, most (or all) of which occur in the SSC region, and the overlap with the IRb region varied from -18 to 16 bp (Fig. 3).

Fig. 3
figure 3

Comparison of the border regions of the 14 studied Ferula plastomes

According to the sequence identity plots, the 14 sequences were almost identical in their genetic structure and showed a very high degree of conservation (Fig. 4). To determine divergent hotspots, nucleotide diversity (Pi) values were calculated (Fig. 5, Table S5), yielding a maximum value of 0.01019 in ycf1. The SSC area showed the maximum nucleotide diversity followed by the LSC region, and the IR regions had the lowest Pi value. Additionally, six highly divergent regions (> 0.006) were detected in the LSC region (rps16/trnQ-UUG, trnS-UGA/psbZ, psbH/petB), SSC region (ycf1/ndhF, rpl32, ycf1), and IR region (0).

Fig. 4
figure 4

Sequence identity plots of the newly sequenced chloroplast genomes

Fig. 5
figure 5

Sliding window analysis of the newly sequenced chloroplast genomes of Ferula species

We calculated the Ka/Ks ratios of the 79 common protein-coding genes to reveal selection patterns among the protein-coding genes. The Ka/Ks ratios of most of the genes were less than 0.5 or could not be computed because either the Ka or Ks value was zero; three genes (ccsA, ndhC, and ycf2) had values greater than 1; and the total Ka/Ks ratio of all genes was 0.5331 (Table S6). In addition, we found several annotation errors (ndhH and ccsA) in the previously reported sequences of F. sinkiangensis (MW411057).

Phylogenetic analyses

To determine the phylogenetic relationship of Soranthus Ledeb., Schumannia Kuntz., Talassia Korovin, and Ferula L., 25 chloroplast genomes were used to construct maximum likelihood (ML) and Bayesian inference (BI) phylogenetic trees. These included 10 samples of 10 Ferula species (including F. sinkiangensis, GenBank accession no. MW411057), two samples of Soranthus, two samples of Schumannia, one sample of Talassia, and nine other Apiaceae genera, i.e., Caucalis L., Daucus L., Cuminum L., Anthriscus Pers., Aegopodium L., Cyclospermum Lag., Apium L., Cryptotaenia DC., and Oenanthe L. with an outgroup of Diplopanax stachyanthus Hand.-Mazz (Fig. 6).

Fig. 6
figure 6

a Branch length diagram of the phylogenetic tree. b Phylogenetic tree of the 25 species inferred from maximum likelihood (ML) and Bayesian inference (BI) analyses based on the complete plastomes. The Shimodaira-Hasegawa-like support values approximate the likelihood ratio test (only F. oopoda and F. gigantea had SH-aLRT values below 80 in the terminal branch), and ultrafast bootstrap values (UFBS ≥ 95%, on the right) are shown on the branches. Green indicates two sequences of S. meyeri (F. sibiraca), blue indicates one sequence of S. karelinii (F. karelinii), and red indicates two sequences of T. transiliensis (F. transiliensis)

The ML and BI topologies were highly supported. Ten selected genera formed 10 monophyletic groups, all of which had support values of 100 or 1 in the ML and BI trees, respectively. Ferula was divided into three main lineages (A, B and C) with maximal support (PP = 1, BS ≥ 97%), and three genera (Soranthus meyeri, Schumannia karelinii, Talassia transiliensis) were clustered into Ferula. Lineage A contained 11 Ferula species, S. meyeri, T. transiliensis, and S. karelinii. Within this lineage, S. sibirica and S. karelinii are sister species, and F. sinkiangensis and F. litwinowiana are sister species. Lineages B and C contained only F. equisetacea and F. olivacea, respectively. Moreover, Ferula and four genera of Apiaceae formed a monophyletic group.


Comparison of Ferula plastid genome

Plastomes are considered an effective means used in taxonomic and evolutionary studies to assess evolutionary relationships and compare genome structure at different taxonomic levels [34,35,36,37,38,39,40,41,42]. Generally, the plastomes are highly conserved in genome structure, gene order, and gene content [32]. In this study, all 14 plastomes are divided into four regions consisting of an LSC (84,839–85,598 bp), an SSC (17,559–17,687 bp), and two IRs (31,392–31,880 bp). The comparative analysis of 14 complete plastomes showed great similarities in terms of genome length (165,607–167,013 bp), structure, IR/SC borders and GC content (37.8–38.0), the equal number of CDs, rRNA, and tRNA genes, and no rearrangement or a good collinearity relationship among them (Fig. 1; Table S1), indicated that the Ferula are relatively conserved.

Although the IR region is thought to be the most conserved region in the chloroplast genome, contraction and expansion of the IR region is common, and is the main reason for the variation in chloroplast genome size [44,45,46]. The junction of IRb/LSC located at ycf2 gene is defined as the type without any expansion or contraction [47]. In this study, we observed that 14 sequenced complete plastomes exhibited significant IR expansion (Fig. 3). All the species expanded into rps19 at the IRb/LSC junction region, contributing to rps19 fragment in the IRa/LSC region, and they also expanded into ycf1 at the IRb/SSC junction region, leading to an overlap between the ycf1 pseudo-gene and ndhF. This was consistent with previous studies, in which the pseudogenes ycf1 and rps19 were produced by contraction and expansion of the IR region in angiosperms [48,49,50].

RSCU value is the ratio of specific codon usage frequency to desired frequency, which can eradicate the influence of amino acid composition on codon usage and promotes the detection of synonymous codons [51, 52]. Generally, the content of A/T was higher than that of G/C in plastomes codons and A/T is preferred in the third codon position [53], the bias also showed in the Ferula plastomes (Fig. 2). Leucine was encoded by 6 codons, the order of codon preference was UTA > CUT > UTG > CUA > CUC > CUG, which following previous studies [54, 55]. The analysis of RSCU can provide a basis for studying the specific mechanism of synonymous codon bias preference in different species, which plays a crucial role in molecular biology basis research [56, 57].

As a primary source of molecular markers, SSRs have been widely used in Ferula genetic diversity studies because of their high polymorphism rate and abundant variation at the species level [58, 59]. In our study, we identified 837 repeats (Table S3) and 1,061 SSRs (Table S4) in the 14 Ferula samples. In which, the single nucleotide and dinucleotide repeats were common, which is consistent with the results of previous studies [55, 60]. In general, during the evolutionary process of species, most repeated sequences in the genome are distributed in the non-coding region and retain as little genetic information as possible to improve its genetic efficiency. Therefore, repeat sequences play an important role in species evolution [61,62,63]. The repeats found in the 12 analyzed species indicate genetic variation among the Ferula species. In addition, we also observed that the poly (A/T) SSRs were typically most common, while poly (C/G) repeats were extremely rare. These results are consistent with those of a previous study and verify the hypothesis that cpSSRs generally consist of short polyadenine (polyA) or polythymine (polyT) repeats and rarely contain tandem guanine (G) or cytosine (C) repeats [64,65,66].

Divergent hotspots play a significant role in species identification and phylogenetic information. Moreover, IR regions often show lower sequence divergence than SSC and LSC regions [67], this probably due to higher mutation rates lead to rapid genome evolution compared to other regions [68]. In our study, this phenomenon was evident that the SSC area showed the maximum nucleotide diversity followed by the LSC region, and the IR regions had the lowest Pi value (Fig. 5, Table S5). And rps16/trnQ-UUG, trnS-UGA/psbZ, psbH/petB, ycf1/ndhF, rpl32, ycf1 were detected as the most divergent regions (Pi > 0.006) across all tested plastomes, suggesting that these variable loci can be used as important references and potential molecular markers for future studies on the evolution and diversity in Ferula. Generally, the Ka/Ks ratio is used to divide genes into positive selection, neutral evolution, and purification, with a limit of one [69]. Previously studies indicated that Ka/Ks ratios mostly are lower due to synonymous nucleotide substitutions rates that occur more often compared to nonsynonymous substitutions rates [70]. The genes with the highest Ka/Ks variability can be used as candidate barcodes to diferentiate species and in the future applied to perform phylogenetic and phylogeographic analyses [71]. Our study suggests that 76 common protein-coding genes were under purifying selection, which indicates the typical evolutionary conservation of plant plastid genes [55, 72, 73], and three genes (ccsA, ndhC, ycf2) were under weak positive selection (Table S6), ycf2 have been proved to be pseudogenized in many studies [74] and ccsA was located in one of the most divergent regions, possibly as a discriminating DNA barcode for Ferula species.

The relationships between Soranthus, Schumannia, Talassia and Ferula

Based on the anatomical morphological characteristics of sclerosing cell layers in the mesocarp, the genera Soranthus, Schumannia and Talassia have been proposed to be located under the genus Ferula [25, 75], all of which are recognized in the Flora of China [76]. It is easily distinguished Ferula from Soranthus and Schumannia by gross morphology and inflorescence structure, combined with the presence of luteolin 7-glycosides in the leaves, that seems reasonable to combine the two genera into Soranthus [77]. Also, Talassia tends to be incorporated into Ferula because insignificant morphological differences, although a large extent similarity between T. transiliensis and F. conocaula in the spectrum of leaf flavonoids [77]. Through a comparative study of plant external morphology, fruit anatomy, and pollen morphology, Qin and Shen [27] suggest that Talassia should be an independent genus and agreed to combine the other two monotypic genera. However, the above four genera have been suggested to merge into one genus according to the presence or absence of coumarins [78]. Recentlly, the molecular phylogeny of Ferula constructed Kurzyna-Młynik et al. [7] and Panahi et al. [16, 29] based on nrDNA ITS and cpDNA sequences (the rps16 intron, the rpoC1 intron and the rpoB-trnC) indicated that Soranthus, Schumannia, and Talassia were embedded in Ferula with low support values. In our study, 15 sequences (including S. meyeri, S. karelinii and T. transiliensis) covered all of the branches except the subgenera Ferula (including section Ferula and section Stenocarpa) according to the latest Ferula phylogenetic tree [16]. Our results show that all those three species representing the genera Soranthus, Schumannia and Talassia were embedded in Ferula based on phylogenetic trees with high bootstrap values (Fig. 6). The species S. meyeri and S. karelinii were clustered into section Soranthus (PP = 1, BS = 100%), and T. transiliensis and F. renardii clustered into section Glaucoselinum (PP = 1, BS = 100%), which was coincident with Panahi [16] while with higher support values. Therefore, we support the standpoint of sinking Soranthus Ledeb., Schumannia Kuntz., Talassia Korovin into synonymy of Ferula L.

Plastomes might provide new insight on phylogenetic relationships in Ferula

As one of a complex taxonomic genus within Apiaceae, the system of Ferula is paid attention at the morphological and molecular levels [7, 9,10,11, 13,14,15,16]. All those efforts on taxonomic systems have contributed greatly to understanding of the genus Ferula. Kurzyna-Młynik et al. [7] published the first molecular phylogeny for Ferula to solve the relationship among Dorema, Ferula and Leutea, in which nrDNA ITS sequences were used to construct a phylogenetic tree revising Dorema and Leutea to Ferula and transferring Ferula to Scandiceae from Peucedaneae. Later, nrDNA ITS sequences and three fragments of cpDNA (the rps16 intron, the rpoC1 intron and the rpoB-trnC) were used to explore the relationship among the three genera, and it was found that Dorema was incorporated into Ferula and Leutea independently [16, 29]. Although these results provide an important foundation for the identification and classification of Ferula species, all previous studies have been based on relatively short sequences with low support values owing to the relatively limited number of nuclear/chloroplast genes. In addition, nrDNA and plastid DNA are highly incongruent, and intense reticulate evolution in Ferula means that proposing an unambiguous hierarchical classification system is almost impossible [16]. Furthermore, many species of Ferula have not been specifically addressed, and many only broadly grouped into branches.

Notably, studies based on plastomes can provide new insights into the phylogenetic relationships between species. For example, Clerodendranthus spicatus is closely related to two Lamiacea species, Tectona grandis L.f. and Glechoma longituba (Nakai) Kuprian. [79]; Juglandaceae is monophyletic, and Carya cathayensis Sarg. is a sister to C. kweichowensis Kuang & A.M.Lu and C. illinoinensis (Wangenh.) K.Koch [66]; and Fagus longipetiolata Seemen and F. engleriana Seemen ex Diels form a close relationship [41]. Here, we performed phylogenetic analyses for Ferula and other genera of Apiaceae using complete plastomes, and we recognized Ferula as a monophyletic group with the integration of Soranthus, Schumannia, Talassia (PP = 1, BS = 100%). Within Ferula, we recovered three main lineages in agreement with Panahi et al. [16], who proposed a new classification based on morphological characteristics and sequence data (nrDNA ITS sequences and three cpDNA fragments). This classification divides Ferula into four subgenera and 10 sections. In addition, Caucalis, Daucus, Cuminum, and Anthriscus were all typical of Scandiceae and formed a monophyletic system with Ferula. This provides strong evidence and support for the transfer of Ferula from the Peucedaneae to the Scandiceae [7]. However, we also observed some differences. When added into Panahi et al.’s phylogenetic tree, F. sinkiangensis was clustered into the Scorodosma branch with the sister species F. kelifi. Based on our results, F. sinkiangensis is separated from F. kelifi, being clustered with F. litwinowiana in the Merwia branch. Further research is needed to confirm this phenomenon. Overall, our work demonstrates that plastome studies can provide highly useful information for future phylogenetic, taxonomic, and evolutionary studies of Ferula.


We obtained 14 complete cp genome sequences from 12 Ferula species (including Soranthus, Schumannia, and Talassia) and compared them based on genome structure, gene content, and gene sequences. Some hotspots in the LSC and SSC regions were identified, which may provide useful markers for phylogenetic analysis. Notably, the Gene ccsA can be used as a DNA barcode for Ferula species. Our phylogenetic analysis showed a tight connection between Soranthus Ledb., Schumannia Kuntz., Talassia Korov., and Ferula L., indicating that treatment as separate genera is unreasonable. Instead, their phylogenetic relationship, which is now well resolved, strongly supports that they can be considered synonymous with Ferula. This new genomic information not only contributes to the better development and utilization of Ferula but also provides a basis for further understanding the evolutionary, genetic, and phylogenetic relationships of this important genera.

Materials and methods

Plant materials and DNA extraction

Fourteen samples were collected from the field and herbaria (Table S1). Of these, five specimens were taken from the specimen museum of the Xinjiang Institute of Ecology and Geography, Chinese Academy of Sciences (XJBI), one was obtained from the Komarov Botanical Institute of RAS (LE), five specimens were taken from the National Herbarium of Uzbekistan (TASH), and three were collected from the field in Tajikistan. Leaf samples were dried in silica gel and stored at -20 °C for DNA extraction. DNA extraction was performed using a plant genome extraction kit (DP320) from Tiangen Biochemical Technology (Beijing) according to the manufacturer’s instructions.

DNA sequencing and genome assembly and annotation

The extracted DNA was sent to a sequencing company for automatic sequencing using the NEBNext Ultra II DNA Library Prep Kit for Illumina (New England BIolabs) [80]. DNA extracts were quantified and sheared into approximately 500 base pair (bp) fragments for library construction using standard protocols (NEBNext Ultra IITMDNA Library Prep Kit for Illumina). Paired-end sequencing from both ends of 150 bp fragments was performed on the Illumina HiSeq X Ten platform at the Molecular Biology Experiment Center, Germplasm Bank of Wild Species in Southwest China, to generate no less than 2 GB data for each individual.

The paired-end reads were filtered using the GetOrganelle pipeline ( to obtain plastid-like reads [81] and then assembled using SPAdes version 3.10 [82]. A complete circular assembly graph was checked and further extracted using Bandage version 0.8.1 [83]. The genomes were automatically annotated using CpGAVAS [84], PGA (, and then manually adjusted using Geneious version 9.1.7 [85]. The chloroplast sequences generated in this study have been submitted to GenBank (Table S1). Circular genome maps of all 14 plastomes were also obtained using the Organellar Genome DRAW (OGDRAW) tool [86].

Codons, repeat sequences, and simple sequences repeat analysis

The protein-coding genes were extracted for codon analysis. The final dataset included 86 protein-coding genes from each species. Codon usage and relative synonymous codon usage (RSCU) values were calculated using JSHYCloud ( A heatmap of all the RSCU values of the 14 plastomes was produced using ClustVis [87]. Using the parameters of a Hamming distance of 3, a minimum repeat size of 30 bp, and a maximum repeat size of 5,000 bp, REPuter was used to identify the size and location of four types of repeat sequences (i.e., forward, palindromic, reverse, and complement) [88]. Simple sequence repeats (SSRs) were detected using the online MISA software ( with minimum repeat number settings of 10, 5, 4, 3, 3, and 3 for mononucleotides, dinucleotides, trinucleotides, tetranucleotides, pentanucleotides, and hexanucleotides, respectively.

Genome comparison with other Ferula species and selective pressure analysis

Sequence divergence among the 14 chloroplast (cp) genomes was compared using Mafft (version 7.0) [89], IRscope ( and Mauve [90]. DnaSP [91] was used to calculate nucleotide divergence values using the sliding window method, with a window length of 800 bp and a step size of 200 bp. Selective pressures were analyzed for 79 common protein-coding genes among 15 Ferula species (including one published plastome). The ratio of nonsynonymous to synonymous nucleotide substitution rates (Ka/Ks) was calculated using DnaSP.

Phylogenetic analysis

We used 25 complete plastome sequences to infer the phylogenetic relationships of Ferula. After comparison with Mafft, Trimal [92], and Phylosuite [93] were used to trim areas with poor quality. The phylogenetic tree was then constructed using RaxML-HPC v.8 [94] and the maximum likelihood method with 1,000 replicates and the GTRGAMMA model. After screening for the best model using jModelTest2 [95], MrBayes 3.2.7a [96] was used to construct a Bayes tree, and the selected models for the complete plastome sequences in BI analyses were TPM1uf + I + G, and iTOL [97] and FigTree 1.4.2 [98] were used to construct the phylogenetic tree.

Availability of data and materials

Fourteen newly sequenced and annotated plastomes have been submitted into NCBI ( Accession numbers: ON324036-ON324048, OM993535.


  1. Pimenov MG, Leonov MV. The genera of the Umbelliferae: a nomenclator. Kew: Royal Botanic Gardens; 1993.

    Google Scholar 

  2. Sina AA. Kanon vrachebnoy nauki (Canon). In: Zaxidov TZ, editor. Book. 1. Tashkent: Academy of Sciences of the UzSSR; 1954. p. 549.

  3. Shen GM. Chinese herbal medicine series: Ferula. Urumqi: Xinjiang people’s publishing house; 1986.

    Google Scholar 

  4. Nazari ZE, Iranshahi M. Biologically active sesquiterpene coumarins from Ferula species. Phytother Res. 2011;25(3):315–23.

    Article  CAS  PubMed  Google Scholar 

  5. Mahendra P, Bisht S. Ferula asafoetida: Traditional uses and pharmacological activity. Pharmacogn Rev. 2012;6(12):141–6.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  6. Ajani Y, Ajani A, Cordes JM, Watson MF, Downie SR. Phylogenetic analysis of nrDNA ITS sequences reveals relationships within five groups of Iranian Apiaceae subfamily Apioideae. Taxon. 2008;57(2):383–401.

    Google Scholar 

  7. Kurzyna-Młynik R, Oskolski AA, Downie SR, Kopacz R, Wojewódzka A, Spalik K. Phylogenetic position of the genus Ferula (Apiaceae) and its placement in tribe Scandiceae as inferred from nrDNA ITS sequence variation. Plant Syst Evol. 2008;274(1–2):47–66.

    Article  CAS  Google Scholar 

  8. Safina LK, Ostroumova TA, Pimenov MG. Carpology of the species of Ferula subgen. Merwia(Umbelliferae-Apioideae) and some taxonomic implications. Nord J Bot. 2015;33(2):140–50.

    Article  Google Scholar 

  9. Boissier PE. Flora orientalis sive, enumeratio plantarum in Oriente a Graecia et Aegypto ad Indiae fines hucusque observatae 2. Genève, Basel & Lyon: H.Georg; 1872. p. 1159.

  10. Drude CGO. Umbelliferae. In: Engler A, Prantl K, editors. Die natürlichen Pflanzenfamilien, vol. 3. Leipzig: Verlag von Wilhelm Engelman (Druck von Breitkopf & Härtel in Leipzig); 1898. p. 63–250.

  11. Korovin EP. Generis Ferula (Tourn.) L. monographia illustrata. Tashkent: Academiae Scientiarum UzRSS; 1947. p. 91.

  12. Chamberlain DF, Rechinger KH. Ferula L. In: Hedge IC, Lamond JM, Rechinger KH, editors. Umbelliferae, Flora Iranica, vol. 162. Graz: Akademische Druck- und Verlagsanstalt; 1987. p. 387–426.

  13. Safina LK, Pimenov MG. The carpoanatomical features of the species of the genus Ferula of the subgenus Peucedanoides (Apiaceae) in connection with the systematics of the genus. Bot Zhurn (Leningrad). 1983;68:730–9.

    Google Scholar 

  14. Safina LK, Pimenov MG. Feruly Kazakhstana. Alma-ata: Nauka Kazakhskoĭ SSR; 1984. p. 110.

  15. Safina LK, Pimenov MG. Carpology of the species of type subgenus of the genus Ferula and some problems of their systematics. Feddes Repertorium. 2008;101(3–4):135–51.

    Google Scholar 

  16. Panahi M, Banasiak L, Piwczyński M, Puchałka R, Kanani MR, Oskolski AA, Modnicki D, Miłobędzka A, Spalik K. Taxonomy of the traditional medicinal plant genus Ferula (Apiaceae) is confounded by incongruence between nuclear rDNA and plastid DNA. Bot J Linn Soc. 2018;188(2):173–89.

    Article  Google Scholar 

  17. Ledebour CF, Bunge A, Meyer CA. Flora Altaica. Berolini: G. Reimeri; 1829. p. 197–206.

  18. von Bunge A. Beitrag zur kenntniss der flor Russlands und der steppen Central-Asiens. St. Petersburg: Kaiserliche Akademie der Wissenschaften; 1851. p. 359.

  19. Korovin EP. Ferula L. In: Schischkin BK, editor. Flora of the USSR. Moscow and Leningrad: Akad. nauk SSSR; 1951. p. 62–214.

    Google Scholar 

  20. She ML, Pu FD, Pan ZH, Watson MF, Cannon JFM, Holmes-Smith I, et al. Apiaceae (Umbelliferae). In: Wu ZY, Raven RH, editors. Flora of China. Beijing and St. Louis: Science Press and Missouri Botanical garden Press; 2005. p. 1–205.

  21. Shen GM. Apiaceae (Umbelliferae). In: Shen GM, editor. Flora Xinjiangensis. Urumqi: Xinjiang Science & Technology Publishing House; 2011. p. 464–621.

  22. Kuntze O. Plantae orientali-rossicae. Trudy Imp: S.-Peterburgsk. Bot. Sada. 1887;10:35–262.

  23. Tojibaev KSh, Sennikov AN, Lazkov GA, Jang GG, Choi HJ, Chang KS, et al. Checklist of vascular plants of the Tian-Shan Mountain System. Pocheon: Korea National Arboretum; 2021. p. 607.

  24. Korovin EP. The new genera and species of Umbelliferae from Kazakhstan flora. Trudy Instituta Botaniki: Akademiya Nauk Kazakhskoi SSR. 1962;13:242–62.

  25. Pimenov MG, Kirillina NA. The carpology of Soranthus, Ladyginia, Eriosynaphe and Schumannia in connection with the problem of the taxonomic limits of the genus Ferula (Apiaceae). Botanicheskii Zhurnal. 1980;65:1756–66.

    Google Scholar 

  26. Govaerts R, Nic Lughadha E, Black N, Turner R, Paton A. The World Checklist of Vascular Plants, a continuously updated resource for exploring global plant diversity. Sci Data. 2021;8(1):215.

    Article  PubMed  PubMed Central  Google Scholar 

  27. Qin XM, Shen KM. Taxonomic studies on the Genus Ferula and its close genera in Xinjiang. Arid Zone Res. 1990;7(4):23–33.

    Google Scholar 

  28. Hui H, Liu QX, Liu MH. Study on serum classification and genetic relationship of Ferula of Peucedaneae subtribe Ferulinae of Apiaceae in China. J Syst Evol. 2003;41(4):369–80.

    Google Scholar 

  29. Panahi M, Banasiak Ł, Piwczyński M, Puchałka R, Oskolski AA, Spalik K. Phylogenetic relationships among Dorema, Ferula and Leutea (Apiaceae: Scandiceae: Ferulinae) inferred from nrDNA ITS and cpDNA noncoding sequences. Taxon. 2015;64(4):770–83.

    Article  Google Scholar 

  30. Corriveau JL, Coleman AW. Rapid Screening Method to Detect Potential Biparental Inheritance of Plastid DNA and Results for over 200 Angiosperm Species. Am J Bot. 1988;75(10):1443–58.

    Article  Google Scholar 

  31. Jansen RK, Raubeson LA, Boore JL, dePamphilis CW, Chumley TW, Haberle RC, Wyman SK, Alverson AJ, Peery R, Herman SJ, et al. Methods for obtaining and analyzing whole chloroplast genome sequences. Methods Enzymol. 2005;395:348–84.

    Article  CAS  PubMed  Google Scholar 

  32. Ravi V, Khurana JP, Tyagi AK, Khurana P. An update on chloroplast genomes. Plant Syst Evol. 2007;271(1–2):101–22.

    Google Scholar 

  33. Wicke S, Schneeweiss GM, dePamphilis CW, Muller KF, Quandt D. The evolution of the plastid chromosome in land plants: gene content, gene order, gene function. Plant Mol Biol. 2011;76(3–5):273–97.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  34. Yang JB, Tang M, Li HT, Zhang ZR, Li DZ. Complete chloroplast genome of the genus Cymbidium: lights into the species identification, phylogenetic implications and population genetic analyses. BMC Evol Biol. 2013;13:84.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  35. Dong W, Liu H, Xu C, Zuo Y, Chen Z, Zhou S. A chloroplast genomic strategy for designing taxon specific DNA mini-barcodes: a case study on ginsengs. BMC Genet. 2014;15:138.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  36. Ma PF, Zhang YX, Zeng CX, Guo ZH, Li DZ. Chloroplast phylogenomic analyses resolve deep-level relationships of an intractable bamboo tribe Arundinarieae (poaceae). Syst Biol. 2014;63(6):933–50.

    Article  PubMed  Google Scholar 

  37. Coissac E, Hollingsworth PM, Lavergne S, Taberlet P. From barcodes to genomes: extending the concept of DNA barcoding. Mol Ecol. 2016;25(7):1423–8.

    Article  CAS  PubMed  Google Scholar 

  38. Hollingsworth PM, Li DZ, van der Bank M, Twyford AD. Telling plant species apart with DNA: from barcodes to genomes. Philos Trans R Soc Lond B Biol Sci. 2016;371(1702):20150338.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  39. Huang Y, Li X, Yang Z, Yang C, Yang J, Ji Y. Analysis of Complete Chloroplast Genome Sequences Improves Phylogenetic Resolution in Paris (Melanthiaceae). Front Plant Sci. 2016;7:1797.

    Article  PubMed  PubMed Central  Google Scholar 

  40. Xie DF, Yu Y, Deng YQ, Li J, Liu HY, Zhou SD, He XJ. Comparative Analysis of the Chloroplast Genomes of the Chinese Endemic Genus Urophysa and Their Contribution to Chloroplast Phylogeny and Adaptive Evolution. Int J Mol Sci. 2018;19(7):1847.

    Article  PubMed Central  CAS  Google Scholar 

  41. Liang D, Wang H, Zhang J, Zhao Y, Wu F. Complete Chloroplast Genome Sequence of Fagus longipetiolata Seemen (Fagaceae): Genome Structure, Adaptive Evolution, and Phylogenetic Relationships. Life (Basel). 2022;12(1):92.

    CAS  PubMed Central  Google Scholar 

  42. Wang N, Chen S, Xie L, Wang L, Feng Y, Lv T, Fang Y, Ding H. The complete chloroplast genomes of three Hamamelidaceae species: Comparative and phylogenetic analyses. Ecol Evol. 2022;12(2):e8637.

    Article  PubMed  PubMed Central  Google Scholar 

  43. Fu CN, Mo ZQ, Yang JB, Cai J, Ye LJ, Zou JY, Qin HT, Zheng W, Hollingsworth PM, Li DZ, et al. Testing genome skimming for species discrimination in the large and taxonomically difficult genus Rhododendro. Mol Ecol Resour. 2022;22(1):404–14.

    Article  CAS  PubMed  Google Scholar 

  44. Goulding SE, Olmstead RG, Morden CW, Wolfe KH. Ebb and flow of the chloroplast inverted repeat. Mol Gen Genet. 1996;252(1–2):195–206.

    Article  CAS  PubMed  Google Scholar 

  45. Huang J-L, Sun G-L, Zhang D-M. Molecular evolution and phylogeny of the angiosperm ycf2 gene. J Syst Evol. 2010;48(4):240–8.

    Article  Google Scholar 

  46. Zhu A, Guo W, Gupta S, Fan W, Mower JP. Evolutionary dynamics of the plastid inverted repeat: the effects of expansion, contraction, and loss on substitution rates. New Phytol. 2016;209(4):1747–56.

    Article  CAS  PubMed  Google Scholar 

  47. Wen J, Xie DF, Price M, Ren T, Deng YQ, Gui LJ, Guo XL, He XJ. Backbone phylogeny and evolution of Apioideae (Apiaceae): New insights from phylogenomic analyses of plastome data. Mol Phylogenet Evol. 2021;161:107183.

    Article  PubMed  Google Scholar 

  48. Gu C, Ma L, Wu Z, Chen K, Wang Y. Comparative analyses of chloroplast genomes from 22 Lythraceae species: inferences for phylogenetic relationships and genome evolution within Myrtales. BMC Plant Biol. 2019;19(1):281.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  49. Liu X, Chang E-M, Liu J-F, Huang Y-N, Wang Y, Yao N, Jiang Z-P. Complete Chloroplast Genome Sequence and Phylogenetic Analysis of Quercus bawanglingensis Huang, Li et Xing, a Vulnerable Oak Tree in China. Forests. 2019;10(7):587.

    Article  Google Scholar 

  50. Kim S-C, Lee J-W, Choi B-K. Seven Complete Chloroplast Genomes from Symplocos: Genome Organization and Comparative Analysis. Forests. 2021;12(5):608.

    Article  Google Scholar 

  51. Rono PC, Dong X, Yang JX, Mutie FM, Oulo MA, Malombe I, Kirika PM, Hu GW, Wang QF. Initial Complete Chloroplast Genomes of Alchemilla (Rosaceae): Comparative Analysis and Phylogenetic Relationships. Front Genet. 2020;11:560368.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  52. Wanga VO, Dong X, Oulo MA, Mkala EM, Yang JX, Onjalalaina GE, Gichua MK, Kirika PM, Gituru RW, Hu GW, et al. Complete Chloroplast Genomes of Acanthochlamys bracteata (China) and Xerophyta (Africa) (Velloziaceae): Comparative Genomics and Phylogenomic Placement. Front Plant Sci. 2021;12:691833.

    Article  PubMed  PubMed Central  Google Scholar 

  53. Morton BR. Selection on the codon bias of chloroplast and cyanelle genes in different plant and algal lineages. J Mol Evol. 1998;46(4):449–59.

    Article  CAS  PubMed  Google Scholar 

  54. Guisinger MM, Kuehl JV, Boore JL, Jansen RK. Extreme reconfiguration of plastid genomes in the angiosperm family Geraniaceae: rearrangements, repeats, and codon usage. Mol Biol Evol. 2011;28(1):583–600.

    Article  CAS  PubMed  Google Scholar 

  55. Ren T, Li ZX, Xie DF, Gui LJ, Peng C, Wen J, He XJ. Plastomes of eight Ligusticum species: characterization, genome evolution, and phylogenetic relationships. BMC Plant Biol. 2020;20(1):519.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  56. Lu H, Zhao WM, Zheng Y, Wang H, Qi M, Yu XP. Analysis of synonymous codon usage bias in Chlamydia. Acta Biochim Biophys Sin (Shanghai). 2005;37(1):1–10.

    Article  CAS  Google Scholar 

  57. Hassan S, Mahalingam V, Kumar V. Synonymous codon usage analysis of thirty two mycobacteriophage genomes. Adv Bioinformatics. 2009;2009:316936.

    Article  CAS  Google Scholar 

  58. Li WJ, Su ZH, Yang L, Cao QM, Fengi Y. Genetic diversity of the critically endangered Ferula sinkiangensis KM Shen (Apiaceae) and the implications for conservation. Turk J Bot. 2020;44(2):145–52.

    Google Scholar 

  59. Yang L, Hisoriev H, Kurbonova P, Boboev M, Bobokalonov K, Feng Y, Li W. High genetic diversity and low differentiation of endangered Ferula tadshikorum Pimenov in Tajikistan. GECCO. 2021;28:e01627.

    Google Scholar 

  60. Ren T, Yang Y, Zhou T, Liu ZL. Comparative Plastid Genomes of Primula Species: Sequence Divergence and Phylogenetic Relationships. Int J Mol Sci. 2018;19(4):1050.

    Article  PubMed Central  CAS  Google Scholar 

  61. Li B, Zheng Y. Dynamic evolution and phylogenomic analysis of the chloroplast genome in Schisandraceae. Sci Rep. 2018;8(1):9285.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  62. Chen Y, Hu N, Wu H. Analyzing and Characterizing the Chloroplast Genome of Salix wilsonii. Biomed Res Int. 2019;2019:5190425.

    PubMed  PubMed Central  Google Scholar 

  63. Khan A, Asaf S, Khan AL, Al-Harrasi A, Al-Sudairy O, AbdulKareem NM, Khan A, Shehzad T, Alsaady N, Al-Lawati A, et al. First complete chloroplast genomics and comparative phylogenetic analysis of Commiphora gileadensis and C foliacea: Myrrh producing trees. PLoS One. 2019;14(1):e0208511.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  64. Wang L, Wuyun T-n, Du H, Wang D, Cao D. Complete chloroplast genome sequences of Eucommia ulmoides: genome structure and evolution. Tree Genet Genomes. 2016;12(1):15.

    Article  Google Scholar 

  65. Shen X, Wu M, Liao B, Liu Z, Bai R, Xiao S, Li X, Zhang B, Xu J, Chen S. Complete Chloroplast Genome Sequence and Phylogenetic Analysis of the Medicinal Plant Artemisia annua. Molecules. 2017;22(8):1330.

    Article  PubMed Central  CAS  Google Scholar 

  66. Shen J, Li X, Chen X, Huang X, Jin S. The Complete Chloroplast Genome of Carya cathayensis and Phylogenetic Analysis. Genes (Basel). 2022;13(2):369.

    Article  CAS  Google Scholar 

  67. Li W, Zhang C, Guo X, Liu Q, Wang K. Complete chloroplast genome of Camellia japonica genome structures, comparative and phylogenetic analysis. PLoS ONE. 2019;14(5):e0216645.

    Article  PubMed  PubMed Central  Google Scholar 

  68. Tyagi S, Jung JA, Kim JS, Won SY. A comparative analysis of the complete chloroplast genomes of three Chrysanthemum boreale strains. PeerJ. 2020;8:e9448.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  69. Hurst LD. The Ka/Ks ratio: diagnosing the form of sequence evolution. Trends Genet. 2002;18(9):486–7.

    Article  PubMed  Google Scholar 

  70. Yang J, Kang GH, Pak JH, Kim SC. Characterization and Comparison of Two Complete Plastomes of Rosaceae Species (Potentilla dickinsii var. glabrata and Spiraea insularis) Endemic to Ulleung Island, Korea. Int J Mol Sci. 2020;21(14):4933.

    Article  CAS  PubMed Central  Google Scholar 

  71. Dong X, Mkala EM, Mutinda ES, Yang JX, Wanga VO, Oulo MA, Onjolo VO, Hu GW, Wang QF. Taxonomy, comparative genomics of Mullein (Verbascum, Scrophulariaceae), with implications for the evolution of Verbascum and Lamiales. BMC Genomics. 2022;23(1):566.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  72. Lee-Yaw JA, Grassa CJ, Joly S, Andrew RL, Rieseberg LH. An evaluation of alternative explanations for widespread cytonuclear discordance in annual sunflowers (Helianthus). New Phytol. 2019;221(1):515–26.

    Article  CAS  PubMed  Google Scholar 

  73. Zhang X, Deng T, Moore MJ, Ji Y, Lin N, Zhang H, Meng A, Wang H, Sun Y, Sun H. Plastome phylogenomics of Saussurea (Asteraceae: Cardueae). BMC Plant Biol. 2019;19(1):290.

    Article  PubMed  PubMed Central  Google Scholar 

  74. Timme RE, Kuehl JV, Boore JL, Jansen RK. A comparison of the first two sequenced chloroplast genomes in Asteraceae: lettuce and sunflower. United States: Lawrence Berkeley National Laboratory; 2006. p. 1–33.

  75. Pimenov MG. Glaucoselinum section (Schischk.) M. Pimen of genus Ferula L. (Umbelliferae). Moscow, Biologicheskie nauki: Nauchnye doklady vysshei shkoly. 1983;12:74–9.

  76. Shan RH, She ML. Flora Reipublcae Popularis Sinicae, vol. 55. Beijing: Science Press; 1979.

    Google Scholar 

  77. Chen XY, Liu QX. Luteolin glycosides as taxonomic markers in Ferula and related genera. Biochem Syst Ecol. 1989;17(4):309–10.

    Article  CAS  Google Scholar 

  78. Liu QX, Wu MY, Rao GX, Ye JS, Hui H. H-NMR detection of coumarin and its application in the chemical classification of Ferula. J Plant Resour Environ. 1999;8(1):46–51.

    CAS  Google Scholar 

  79. Du Q, Jiang M, Sun S, Wang L, Liu S, Jiang C, Gao H, Chen H, Li Y, Wang B, et al. The complete chloroplast genome sequence of Clerodendranthus spicatus, a medicinal plant for preventing and treating kidney diseases from Lamiaceae family. Mol Biol Rep. 2022;49(4):3073–83.

    Article  CAS  PubMed  Google Scholar 

  80. Meyer M, Kircher M. Illumina sequencing library preparation for highly multiplexed target capture and sequencing. Cold Spring Harb Protoc. 2010;2010(6):pdb prot5448.

    Article  PubMed  Google Scholar 

  81. 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.

    Article  PubMed  PubMed Central  Google Scholar 

  82. Bankevich A, Nurk S, Antipov D, Gurevich AA, Dvorkin M, Kulikov AS, Lesin VM, Nikolenko SI, Pham S, Prjibelski AD, et al. SPAdes: a new genome assembly algorithm and its applications to single-cell sequencing. J Comput Biol. 2012;19(5):455–77.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  83. Wick RR, Schultz MB, Zobel J, Holt KE. Bandage: interactive visualization of de novo genome assemblies. Bioinformatics. 2015;31(20):3350–2.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  84. Liu C, Shi L, Zhu Y, Chen H, Zhang J, Lin X, Guan X. CpGAVAS, an integrated web server for the annotation, visualization, analysis, and GenBank submission of completely sequenced chloroplast genome sequences. BMC Genomics. 2012;13:715.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  85. Kearse M, Moir R, Wilson A, Stones-Havas S, Cheung M, Sturrock S, Buxton S, Cooper A, Markowitz S, Duran C, et al. Geneious Basic: an integrated and extendable desktop software platform for the organization and analysis of sequence data. Bioinformatics. 2012;28(12):1647–9.

    Article  PubMed  PubMed Central  Google Scholar 

  86. Lohse M, Drechsel O, Kahlau S, Bock R. OrganellarGenomeDRAW–a suite of tools for generating physical maps of plastid and mitochondrial genomes and visualizing expression data sets. Nucleic Acids Res. 2013;41(Web server issue):W575-581.

    Article  PubMed  PubMed Central  Google Scholar 

  87. Metsalu T, Vilo J. ClustVis: a web tool for visualizing clustering of multivariate data using Principal Component Analysis and heatmap. Nucleic Acids Res. 2015;43(W1):W566-570.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  88. Kurtz S, Choudhuri JV, Ohlebusch E, Schleiermacher C, Stoye J, Giegerich R. REPuter: the manifold applications of repeat analysis on a genomic scale. Nucleic Acids Res. 2001;29(22):4633–42.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  89. Katoh K, Rozewicki J, Yamada KD. MAFFT online service: multiple sequence alignment, interactive sequence choice and visualization. Brief Bioinform. 2019;20(4):1160–6.

    Article  CAS  PubMed  Google Scholar 

  90. Darling AC, Mau B, Blattner FR, Perna NT. Mauve: multiple alignment of conserved genomic sequence with rearrangements. Genome Res. 2004;14(7):1394–403.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  91. Rozas J, Ferrer-Mata A, Sanchez-DelBarrio JC, Guirao-Rico S, Librado P, Ramos-Onsins SE, Sanchez-Gracia A. DnaSP 6: DNA Sequence Polymorphism Analysis of Large Data Sets. Mol Biol Evol. 2017;34(12):3299–302.

    Article  CAS  PubMed  Google Scholar 

  92. Capella-Gutierrez S, Silla-Martinez JM, Gabaldon T. trimAl: a tool for automated alignment trimming in large-scale phylogenetic analyses. Bioinformatics. 2009;25(15):1972–3.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  93. 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.

    Article  PubMed  Google Scholar 

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

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  95. Darriba D, Taboada GL, Doallo R, Posada D. jModelTest 2: more models, new heuristics and parallel computing. Nat Methods. 2012;9(8):772.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

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

    Article  PubMed  PubMed Central  Google Scholar 

  97. Letunic I, Bork P. Interactive Tree Of Life (iTOL) v5: an online tool for phylogenetic tree display and annotation. Nucleic Acids Res. 2021;49(W1):W293–6.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  98. Rambaut A. FigTree 1.4.2 software, a graphical viewer of phylogenetic trees. Edinburgh: Institute of Evolutionary Biology University of Edinburgh; 2014.

Download references


We thank Yu Wenbin, Liu Rong, Xin Yaxuan and Li Jinyue in Xishuangbanna Tropical Botanical Garden of Chinese Academy of Sciences for their assistance to assemble and annotation sequences. Our sincere thanks also go to the anonymous reviewers for their comments and suggestions.


The materials used in this article are not related to plant protection and are within the limits of national laws. And we have permission from The Komarov Botanical Institute of RAS (LE) and The National Herbarium of Uzbekistan (TASH) to enter and collect plant specimens.


The work was supported by the Third Xinjiang Scientific Expedition Program (Grant No.2021xjkk0600), the National Science and Technology Basic Program of China (2019FY100204), the Shanghai cooperation organization partnership and international technology cooperation plan of science and technology projects (2021E01020), the Youth Innovation Promotion Association Foundation of the Chinese Academy of Sciences, China (No. 2019429), and Taxonomic revision of polymorphic plant families of the flora of Uzbekistan (FZ-20200929321). 

Author information

Authors and Affiliations



WJL, KHT and YMZ conceived and designed the work. LY and WJL analyzed the sequence data. WJL, OA and KFS provided the materials. LY wrote the manuscript. WJL and KHT revised the manuscript. All authors gave final approval of the paper.

Corresponding author

Correspondence to Wen-Jun Li.

Ethics declarations

Ethics approval and consent to participate

Not applicable. No specific permits were required for the collection of specimens for this study. This research was carried out in compliance with the relevant laws and all methods were performed in accordance with the relevant guidelines and regulations.

Consent for publication

Not applicable.

Competing interests

The authors declare that they have no competing interests.

Additional information

Publisher’s Note

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

Supplementary Information

Additional file 1:

Table S1. All the materials used in this article.

Additional file 2:

Table S2. Codon usage and relative synonymous codon usage (RSCU) values of protein-coding genes of the 14 Ferula plastomes.

Additional file 3:

Table S3. Distribution of repeat sequences in the 14 studied Ferula plastomes.

Additional file 4:

Table S4. Distribution of simple sequence repeats (SSRs) in the 14 studied Ferula plastomes.

Additional file 5:

Table S5. Nucleotide variability (Pi) of Ferula species.

Additional file 6:

Table S6. Non-synonymous to synonymous nucleotide substitution rates (Ka/Ks) of Ferula species.

Rights and permissions

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

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Yang, L., Abduraimov, O., Tojibaev, K. et al. Analysis of complete chloroplast genome sequences and insight into the phylogenetic relationships of Ferula L.. BMC Genomics 23, 643 (2022).

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: