1. Introduction
The origin of the Chloroplast (cp) can be traced back to more than one billion years ago as a result of Cyanobacterium endosymbiosis [1,2,3]. It is an organelle commonly found in the cytoplasmic matrix that is useful in the process of photosynthesis hence sustaining life on Earth [4,5]. The Chloroplast (cp) is a semi-autonomous organelle having its own genetic material, but some of its proteins are encoded in the nuclear genome [6]. The Chloroplast (cp) genome of angiosperms is mostly double-stranded circular structure containing four parts: a large single copy (LSC) region, a small single copy (SSC) region and a pair of inverted repeats (IRa/IRb) regions with the same sequence in opposite directions [7,8]. The first complete chloroplast genome to be sequenced was that of Nicotiana tabacum [9]. In higher plants, the plastome size is relatively smaller in size ranging between 120 and 180 kb in most terrestrial plants, and having highly conserved sequences that encodes approximately 110–130 genes [6,9,10]. These genes are involved in various functions including replication, translation and photosynthesis. In angiosperms, the plastome is maternally inherited. However, about 20% of the sequences may be inherited from patrilineage or from both parents [11,12]. Compared to nuclear genome, the plastome is relatively conserved and stable, with no recombination and low nucleotide substitutions. Therefore, Plastomes are very informative and valuable sources of genetic markers for molecular systematics and phylogenetic analysis [5,13,14]. Some genes have been used for DNA barcoding studies in plants, e.g rbcL, matK, and ycf1 [15,16]. In recent years, with the rapid development of sequencing technology and its affordability, more plastomes have been sequenced successfully [17,18,19,20,21]. Hence, chloroplast genomes have become a new and valuable tool for phylogenomic studies.
Fortunella venosa (Champ. ex Benth.) C.C.Huang is a perennial evergreen shrub in the flowering plant family Rutaceae. This species is endemic to China with its distribution area slightly overlapping to north with the tetraploid species F. hindsii (Champ. ex Hook.) Swingle. Its distribution range is relatively narrow occurring in Nanping (Fujian), Yongfeng (Jiangxi), and recently found in Ningyuan, Chaling, Guidong, and other counties in Hunan province [22,23]. Due to climate change and anthropogenic activities, the wild population of F. venosa is decreasing rapidly [24]. It was listed as an endangered species in the China Species Red List (vol.1 Red list) [25].
The genus Fortunella was described by Swingle in 1915, and currently there are six described species [26,27]. Flora Reipublicae Popularis Sinicae recorded five species and a few hybrids from China, F. venosa is one of them [22]. In the Flora of China (English edition), Fortunella was incorporated into the genus Citrus [28] and F. hindsii (Champ. ex Benth.) Swingle, F. japonica (Thunb.) Swingle, F. margarita (Lour.) Swingle, F. venosa (Champ. ex Hook.) C.C.Huang were treated as synonyms of C. japonica Thunb. Hence, the taxonomy and phylogeny of the genus is complex and controversial. In addition, due to the great number of varieties and easy hybridization of Citrus plants, the taxonomy of the genus and its phylogenetic relationships with Citrus L. has always been a concern for the taxonomists [26,29,30,31].
Members of Fortunella have high value in food, medicine and ornamental [24]. The introduction and cultivation of Fortunella species in different regions and the various types of highly hybrid germplasm resources has made it difficult to identify F. venosa [32]. A complete chloroplast genome sequencing will be helpful to solve the uncertainty of the taxonomic status among species. Currently, only the complete Chloroplast genome of F. japonica (GenBank accession no.: MN495932) has been sequenced and reported among the species in the genus Fortunella. Hence, in this study the initial complete Chloroplast genome of F. venosa was sequenced and reported, and a method for assembling, splicing and analysis of the complete Chloroplast genome of F. venosa was proposed. Furthermore, the Chloroplast genome of F. venosa was compared with other nine Rutaceae species from the NCBI database. Codon usage, repeat sequences, selection pressure, and phylogenetic relationships were analyzed. Sequencing of F. venosa plastome not only provides a theoretical basis for the phylogenetic relationship and its related taxonomic issues, but also provides an important foundation for the conservation and sustainable utilization of F. venosa resources.
2. Materials and Methods
2.1. Material Acquisition and Chloroplast Genome Sequencing
The plant materials were collected from Zhuzhou prefecture, Hunan province (Co-ordinates: 26°47′00.69″ N, 113°29′38.59″ E), China, on 9 April 2019. The fresh young leaves were collected and immediately placed in sealed bags and dried with silica gel. The voucher specimen (K.M. Liu, T. Wang 772949) was deposited in the Herbarium of Hunan Normal University (HNNU). Genomic DNA was extracted from 0.5 g of dried plant leaves with the conventional cetyltrimethylammonium bromide (CTAB) method [33], and sequenced using the second-generation sequencing platform illumina of the Novogene Company in Beijing, China.
2.2. Assembly of the Genome
The original quality of the sequences was evaluated using the FastQCv0.11.7 software [34]. Assembly was done using GetOrganellev1.6.22d [35] with default parameters. The GetOrganelle software is an advanced tool that provides a large number of scripts and libraries of Whole Genome Sequencing (WGS) read data, manipulating and disentangling assembly graphs, and generating reliable organelle genomes, accompanied by labeled assembly graphs for user-friendly manual completion and correction. Using the GetOrganelle software, we first filtered the plastid reads, then performed a de novo assembly, purified the assembly, and finally generated the complete Chloroplast genome [36,37,38]. Redundant sequences were then removed for subsequent genomic analysis. The final assembly map was visualized using Bandage [39] to identify automatically generated plastid genomes.
2.3. Annotation of the Genome
The assembled complete Chloroplast genome was annotated using the Plastid Genome Annotator (PGA) [40] and Strawberry Perl, using Amborella trichopoda (GenBank accession number: GCA_U000471905.1) as the initial reference. The published genomes of Citrus maxima (MN782007) and C. limon (MT880608) of the family Rutaceae were used as control for further annotation confirmation. Annotation tool in Geneious was used to manually correct and supplement problematic annotations. The whole Chloroplast genome circular map was drawn by using Organelle Genome Draw (OGDRAW) online software [41,42].
2.4. Repeat Sequence and Codon Usage
Dispersed repeats (forward, reverse, complementary, palindromic repeat sequences) in the complete Chloroplast genome sequence was analyzed using the REPuter online program (
2.5. Comparative Genome Analysis and Sequence Divergence
Using Fortunella venosa as a reference, the divergence within the ten Chloroplast genomes was analyzed using mVISTA tool [49,50]. The species sequences used included F. venosa and 9 other Rutaceae species; F. japonica (MN495932), Citrus aurantifolia (KJ865401), C. aurantium (MT702983), C. hongheensis (MT880607), C. cavaleriei (MT880606), C. limon (MT880608), C. maxima (MN782007), C. medica (MT106673), C. sinensis (DQ864733). To analyze the rearrangements and inversions within the boundary region of F. venosa, an insertion program Mauve in Geneious8 (Biomatcrs Ltd. Auckland, New Zealand) was used. The IRscope (IRscope.shinyapps.io/Chloroplot/) [51] software was used to analyze the expansion and contraction of IR boundary of the 10-representative species, and compared the differences within the IR boundaries. DnaSP v.5.0 [52] software was used to calculate nucleotide polymorphism (Pi), with the parameter set as follows: window length of 600 bp and the distance between each site (step size) was 200 bp. This was used to construct a polymorphic site line chart, and find fragments with high polymorphism among the Chloroplast genomes.
2.6. Adaptive Evolution and Substitution Rate
In order to analyze the rate of evolutionary changes in the Chloroplast genome of Fortunella venosa, the CDS sequence was extracted using geneious with Citrus aurantifolia as reference. The protein-coding sequences of the 10 Rutaceae species were extracted using PhyloSuite [53], MAFFT to automatically remove the stop codon. PhyloSuite was used to construct the maximum likelihood phylogenetic tree. GTR was selected as the best-fit model, and no outgroup was specified.1000 Bootstraps were performed to construct the phylogenetic unrooted tree. The PAML file and Newick file are imported into EasyCodeML for selective pressure analysis. Using the PAML v4.7 package of the EasyCodeML software [54,55], the positive selection pressure, non-synonymous (DN) and synonymous (DS) substitution rates, and their ratio (ω = DN/DS) of 10 Rutaceae species Plastomes were evaluated. The site-specific model in the software (M0 vs. M3, M1a vs. M2a, M7a vs. M8, and M8a vs. M8) were compared. In order to evaluate the adaptive evolution of Chloroplast genes, the computational likelihood ratio test (LRT) and ω were used to analyze the selection pressure of protein-coding genes in 10 plants.
2.7. Phylogeny
To determine the phylogenetic position and relationship of Fortunella venosa, a phylogenetic tree was reconstructed using 28 other species Chloroplast genome sequences downloaded from NCBI database with Melia azedarach Linn. as the outgroup. The outgroup was chosen according to the current APGIV system of classification (
3. Results
3.1. Analysis of Chloroplast Genome Structure
A genome paired-end sequencing obtained a total of 8,768,734 reads of 150 bp in length were obtained from Chloroplast genome sequencing, of which 3,244,455 bp reads were used for chloroplast genome assembly, accounting for 37% of the total. The base coverage reads used to assemble the Chloroplast genome was 793.65 times. The chloroplast genome of Fortunella venosa (GenBank accession No. MZ457935) has been submitted to the GenBank of the National Center for Biotechnology Information (NCBI). The complete chloroplast genome of Fortunella venosa had a size of 160,265 bp. The plastome of F. venosa is a typical circular four-part structure consisting of a large single copy region (LSC, 87,597 bp), a small single copy region (SSC, 18,732 bp) and two inverted repeat regions (IRa and IRb, 26,968 bp each) (Figure 1 and Table 1). A total of 134 functional genes, including 89 protein-coding genes (CDS), 8 ribosomal RNA (rRNA) genes, and 37 transfer RNA (tRNA) genes, were detected in F. venosa cpDNA (Table 1). The LSC region consists of 62 CDS, and 22 tRNAs, whereas, the SSC region is composed of 12 CDS and 1 tRNA. The IR regions is composed of 18 CDS, 14 tRNA and 8rRNA (Figure 1). The total GC-content of the F. venosa chloroplast genome was 38.4%. The IR region had the highest GC content (43.0%), while the LSC and SSC had 36.7% and 33.2%, respectively. The total length of the protein-coding region, tRNA and rRNA were 79,983 bp, 2792 bp, and 9044 bp, respectively, accounting for 49.9%, 1.7%, and 5.6% of the total length of the chloroplast genome. There were 21 genes duplicated in the IR region with two copies, including 10 protein coding genes (rpl22, rps19, rpl2, rpl23, ycf2, ycf15, nbhB, rps7, rps12, ycf68), seven tRNA genes (trnI-CAU, trnL-CAA, trnV-GAC, trnI-GAU, trnA-UGC, trnR-ACG, and trnN-GUU), and four rRNA genes (rrn16, rrn23, rrn4.5, and rrn5). There were 17 genes with introns, 15 genes had one intron (rps16, trnG-UCC, atpF, rpoC1, trnL-UAA, trnV-UAC, petB, petD, rpl16, trnI-GAU, trnA-UGC, ndhA, trnA-UGC, trnI-GAU, rps12, ndhB, rpl2) while two genes (ycf3, clpP) had two introns. Chloroplasts have maintained an autonomous genome that encodes important proteins required for their photosynthesis and different housekeeping functions. According to the function of genes, they can be divided into four categories, as shown in Table 2. There is a variation in the Chloroplast genomes of different species in terms of length, GC content and even the evolutionary rates. The comparison of Chloroplast genomes of ten species of Rutaceae is shown in Table 1.
3.2. Repeat Sequence Analysis
A total of 50 long repetitive sequences were detected in the Chloroplast genome of Fortunella venosa by REPuter, including 22 forward repeats (F), 7 reverse repeats (R), 19 palindromic repeats (P), and two complementary repeats (C). Forward repeats were the most abundant, followed by the palindromic repeats within all the species. The least abundant repeats were complementary repeats (Figure 2). Most of the repeat sites were located in the non-coding region of LSC, and some of them were located in rpoB, psaB, trnV-UAC, trnS-GCU, and trnL-UUA. Six repeat sites were located in the IR region and two in the SSC region. Analysis of the experimental data showed that most of the repeat sequences were 30–40 bp in length, with the longest being a palindrome repeat sequence with 54 bp. This repeat sequence was located between trnH-GUG and psbA section in the LSC region.
A total of 37 tandem repeats were detected by Tandem Repeats Finder, three repeats of which were longer than 30 bp in length and the others were between 1 bp and 26 bp. 20 repeat units reported mismatches and 10 had indels.
3.3. SSR Analysis
In this study, a total of 108 SSR loci were detected in the Chloroplast genome of Fortunella venosa. Among them, 74 were mononucleotides, five were dinucleotides, 15 were trinucleotide, 11 were tetranucleotide, two were pentanucleotides, and one was hexanucleotide (Figure 3). Most of these SSR loci were distributed in the Chloroplast genome, accounting for 74.1% of LSC region. The results are basically consistent with those of the other nine species (Figure 4 and Figure 5). Furthermore, 88.9% of the 108 SSRs were located in the non-coding region, and 11.1% of the rest in the coding region were located in the LSC region.
3.4. Codon Usage Analysis
A total of 89 CDS of the chloroplast genome of Fortunella venosa were used to estimate the relative frequency of synonymous codon usage. A total of 26,699 codons were detected, out of which 2844 (10.65%) encoded leucine and 315 (1.18%) encoded cysteine, which were the most and the least abundant amino acids in the Chloroplast genome of F. venosa, respectively. The most used codon was AUU (1071) encoding isoleucine and least used codon was AUG (1) that encoding methionine. From the analysis of the frequency of synonymous codon usage (RSCU) in the plastome, the codon usage was biased, among which 30 amino acids had a RSCU > 1. Three amino acids, methionine (AUG), serine (UCC), and tryptophan (UGG) do not have codon usage bias (RSCU = 1.00). Among the three stop codons, the use of the stop codon was biased towards UAA (RSCU > 1.00). The relative synonymous codon usage of F. venosa is shown in (Table 3).
3.5. Comparative Genome and Sequence Divergence Analysis
In general, the sequence sizes of these species are similar, ranging from 159,893 bp to 160,996 bp in length. As shown in Figure 6, the sequence similarity is very high, indicating that the Chloroplast genome is highly conserved having translocation and inversion of the genes (See File S2). In the 10 Plastomes, the IR regions were more conserved than LSC and SSC regions. Similarly, the coding regions were more conserved than non-coding regions. The regions that are relatively variable in non-coding section include; trnA(GUG)-psbA psbL-trnG(UGG), petN-psbM, psbE-trnM(CAU), trnL(UAA)-trnF(GAA), ndhC-trnV(UAC). These regions may have rapid nucleotide substitution at the species level, indicating that molecular markers have potential application value in phylogenetic analysis and plant identification (Figure 6).
In this study, the results showed that although Chloroplast genomes are generally conserved in length and genetic structure, they still show significant differences in the IR/SC boundary region (Figure 7). All genes at the border include rps3, rpl22, ndhF, ycf1, trnH. The expansion and contraction of the border region was analyzed for the 10 species. For example, the position of rpl22 gene in Citrus aurantium, C. cavaleriei, C. hongheensis, C. limon, and C. sinensis is located in the IRb region with a distance of 7 bp, 6 bp, 7 bp, 7 bp, and 7 bp from the boundary, respectively. The rpl22 in the other species spans the LSC and IRb regions, and the situation of rpl22 at the boundary of LSC and IRa is also different, the rpl22 gene is missing in C. maxima, C. medica, and Fortunella venosa. The gene ndhF located at the border between IRb and SSC is only 2 bp and 2200 bp in C. medica, and the rest are 31 bp and 2201 bp. The gene trnH located on the border of IRa and LSC is located on LSC but the length from the border varies from 2–65 bp. Ycf1 was lost at the boundary of IRb and SSC in C. medica and F. venosa, and ndhF crossed the boundary of LSC and IRb, but only 2 bp was located at IRb in C. medica, the rest was 31 bp. The length of ycf1 at the boundary between SSC and IRa is 5490 bp to 5505 bp. These results indicate that there is a contraction and expansion of IR region, which can be used for the study of species-specific gene loci.
The results of Dnaspv.5.0 showed that the regions with high nucleotide polymorphism were the LSC and SSC regions, which was basically consistent with the results of mVISTA (Figure 7). The highly polymorphic loci were trnG-GCC, trnfM-UAA, ndhJ, rpl2, rpl23, trnL-CAU, ccsA, ndhD, ycf1, trnN-GUU, and trnR-AGG. The highest value of Pi was 0.01563, recorded by the genes rpl2 and rpl23. The Pi value was more than 0.01, As shown in (Figure 8). Data on specific nucleotide polymorphisms are provided in File S3.
3.6. Adaptive Evolution Analysis
The dN/dS value can be used to measure the evolution rate of a specific gene [57]. This is the ratio of synonymous substitution rate (dS) to nonsynonymous substitution rate (dN) (ω = dN/dS). In the selection pressure analysis, when ω > l, it shows a positive selection, while, when ω = l, it is a neutral selection; if ω < 1, it is a purifying selection. In this study, we found that the model M7 vs. M8 is the most suitable model by EasyCodeML detection. A total of 344 positive selection sites were detected in 79 CDSs of the ten species (see File S4). Among them, the Naive Empirical Bayes (NEB) detected 54 genes loci, encoding 15 genes of selection pressure, accounting for 18. 99% of the total number of genes. The largest number of loci was rpoC2 with 27. Bayes Empirical Bayes (BEB) detected 290 positive selection sites, which encode the selection pressure of 53 genes, respectively accounting for 67.09% of the total number of genes, and rpoC2 has the most loci with 57. In NEB, photosynthesis-related genes ndhI (2) and self-replicating genes rpoC2 (8), rps2 (1), and rps18 (1) had p > 0.99%. In BEB, photosynthesis-related genes ndhB (1), ndhI (2), psbZ (1) and self-replicating genes rpoC2 (8), rps18 (1), rps19 (1), rps2 (1) had p > 0.99%. is shown in (Table 4). The results showed that the 10 species were under strong positive selection pressure, the nucleotide substitution rate was faster, and they showed strong adaptive variation to their environment.
3.7. Phylogenetic Analysis
The CDS phylogenetic tree results are shown in (Figure 9), Zanthoxylum was clustered into one branch. Glycosmis, Micromelum, Clausena, and Murraya showed a close relationship and hence formed single clade. Fortunella venosa and F. japonica were clustered together and showed a close relationship to genus Citrus. The two species were found to be closely related. They both show that genus Fortunella and genus Citrus are closely related. The results of the whole Chloroplast genome tree are shown in (Figure 10). In the phylogenetic tree reconstructed using the complete chloroplast genome, more than 95% of the branches have a support value of 100% which supports a close relationship among the species. However, one Citrus branch has a support value of 55.6% and its phylogenetic status needs to be further studied, which are basically consistent with the phylogenetic relationship constructed from CDS (Figure 9).
The Chloroplast genome sequence provides an important resource for phylogenetic research. In order to get a more detailed phylogenetic conclusion, more complete Chloroplast genomes of Fortunella are needed. As a highly primitive group of this genus, the complete Chloroplast genome characteristics of F. venosa is indispensable, which will be subsequently used for Citrus taxa phylogenetic study.
4. Discussion
In this study, we analyzed the complete chloroplast genome of Fortunella venosa and performed a comparative study with 10 Rutaceae species. The Chloroplast genome of Fortunella venosa is a circular structure with a size of 160,265 bp, which is similar to the size of other related species reported [58,59]. All the 10 complete chloroplast genomes of the Rutaceae species displayed attributes that is similar to other angiosperm Chloroplast genomes, with quadripartite structure including the LSC, SSC, and a pair of inverted repeats (IRa and IRb). Although there were no genomic rearrangements with gene order highly conserved, there were differences in the Chloroplast genomes ranging from 160,229 bp–160,265 bp in genus Fortunella, and 159,893–160, 996 bp in genus Citrus, hence suggesting some small genetic differences within the genomes. Previous studies have reported that loss of genes [60], variations in the inverted repeat regions [61], and the intergenic spacer region variation [62] are the three fundamental causes of the variations in the chloroplast genome sizes in plants. Additionally, the chloroplast gene length has also been associated with the genome size [63].
Repetitive sequences play an important role in genome rearrangement and are very helpful for phylogenetic studies [64]. In addition, analysis of various chloroplast genomes has shown that repetitive sequences are essential for the study of indels and substitutions [65]. All the ten Rutaceae species had reverse, compliment, forward and Palindromic repeats which varied in number among all the species. From the analysis, the number of repeats found within the chloroplast genomes indicate that Fortunella venosa and F. japonica are more similar than to the rest of the Citrus species. Studies have linked sequence variations and genome rearrangements to slipped-strand mispairing and improper recombination of the repeat sequences [66]. On the other hand, Simple Sequence Repeats (SSRs), also known as microsatellite sequences, are repeated DNA sequences, widely distributed in the whole Chloroplast genome, having lengths of about 1–6 bp [67]. The inheritance of cpDNA in higher plants is mostly maternal, and the structure of cpDNA is relatively conserved and simple, hence cpDNA SSR is widely used as molecular markers, variety identification and other molecular studies [68]. For example, SSR analysis is important for DNA markers used for population genetic and evolutionary studies [69,70,71]. In this study, we analyzed the SSRs in the Chloroplast genomes. Six categories of perfect SSRs (mono-, di-, tri-, tetra-, penta-, and hexa-nucleotide repeats) were detected in the Chloroplast genome of these ten species. In recent years, more evidence shows that the repetitive structure of genomic DNA is essential, not only important in plant molecular research [72], but also widely used in the study of population genetics of species [73,74,75]. SSR has the advantages of high mutation rate, site specificity and multiple alleles [76,77], which can be used for genetic diversity analysis [78,79].
The relative frequency of synonymous codon usage (RSCU) values in Chloroplast genomes have been shown to be as a result of mutation and selection [80,81], which are crucial in the study of the evolution of organisms. RSCU > 1 indicates a preference for the codon, RSCU < 1 indicates a low usage of the codon, and RSCU = 1 indicates no preference for the codon [82]. The codon usage was biased, among which 30 amino acids have RSCU >1. Three amino acids, methionine (AUG), serine (UCC), and tryptophan (UGG) do not have codon usage bias (RSCU = 1.00). Among the three stop codons, the use of the stop codon was biased towards UAA (RSCU > 1.00). This is basically consistent with the reports of other Chloroplast genomes in Rutaceae [58,59].
Comparative analysis in the 10 Plastomes showed that the IR regions were more conserved than LSC and SSC regions. Similarly, the coding regions were more conserved than non-coding regions. This is a common phenomenon in most angiosperms [83]. There is a variability in the size, position at the boundary regions among the species especially for genes such as rpl22, ndhF, and ycf1. This changes in the sizes and positions of the genes cause changes in the size of the genome, hence comparatively, there is a variation in length and number of genes as shown among the species. Expansion and contraction at the borders of the IR regions are considered important in the Chloroplast genome size and play a vital role in its evolution [84]. Nucleotide diversity among the 10 species of Rutaceae genomes was calculated, indicating that the average nucleotide diversity is 0.00252 (Supplementary File S4). This was comparatively higher as compared to the previous studies that compared the species level and the interspecific nucleotide diversity [85]. Most of the nucleotide diversity sites occurred in the LSC and SSC regions, with the highest peaks being rpl2/rpl23/trnL-CAU (0.016) and ycf1/trnN-GUU/trnR-AGG (0.015). This shows that there are low levels of nucleotide diversity throughout the chloroplast genome.
The genus Fortunella includes four species of the “kumquats” from eastern Asia (China, Hong Kong, and Malay Peninsula). It is traditionally separated from Citrus by quantitative characters, 3–7 (versus 8–18) locules in the ovary with two (vs. 4–12) ovules per locule, and by smaller fruits. In other vegetative, floral, and fruit characters, Fortunella is quite similar to Citrus, including the polyadelphous androecium (character 4) with numerous stamens cohering in bundles, a character more commonly found in Citrus subgenus Citrus. The results of this study (Figure 9 and Figure 10) indicate that Fortunella Swingle and Citrus L. are closely related, but do not support the incorporation of F. venosa into C. japonica. The complete chloroplast genomes have been shown to provide informative sites for resolving phylogenetic relationships of plants, and have been examined as well to be effective in the ability of differentiation in lower taxonomic levels [86]. The ML, BI, PhyML tree showed a very high level of support in our study. This study shows that F. venosa should be an independent species, which is significantly different from F. japonica in terms of morphology (habitat, leaf type, fruit size, etc.). F. venosa is a small shrub, usually no more than 1 m tall (the shortest mature plant is 0.28 m high); the leaves are single leaves (the joints of the petiole and the leaf are not joint); the leaves are usually small, 2–4 (−7) cm long, 1–2 cm wide, wedge-shaped base; petiole short, 1–3 (−5) mm long; flower solitary leaf axils, petals 3–5 mm long; ovary 2–3 compartments; fruit spherical or elliptical, diameter 6–8 mm, Orange-red when mature. On the other hand, F. japonica is a small tree or tree with a height of 2 to 8.5 m, and the main stem is usually slender; the leaf is a single leaflet with joints between the petiole and the leaf; the leaf is larger, 4–9 cm long, 1.5–3.5 cm wide, and a wide wedge-shaped base; The petiole is obviously longer, 6–10 (−15) mm long; the flower is single or 2–3 clusters with leaf axils, the petals are 6–8 mm long; the ovary is 4–6 compartments; the fruit is larger, spherical, 1.5–2.5 cm in diameter, Orange-yellow to orange-red when cooked (Figure 11). Due to the significant morphological difference between F. venosa and F. japonica, it is easy to distinguish the two in the wild. The molecular results obtained in this study provide strong support for the independent systematics status of F. venosa. In this paper, we still use F. japonica and F. venosa as the scientific names according to the Flora of China for better discussion. In addition, none of the research results done so far based on morphology, palynology and molecular biology supports the incorporation of F. venosa into C. japonica, showing that the two species are independent.
5. Conclusions
This paper reports the first complete Chloroplast genome sequence of Fortunella venosa. It provides a more detailed and complete information, laying a foundation for the identification of species in the genus Fortunella and the analysis of genetic differences at the individual level. In Rutaceae, Fortunella is phylogenetically related to Citrus, but the inter-species relationship is complicated. This study confirmed that the molecular phylogeny supports F. venosa as an independent species. Hence the Chloroplast genome proves an important basis for the study of systematic classification. In order to better solve the problem of systematic classification in Fortunella, we need to get more cpDNA sequence information of Fortunella. Furthermore, the variations among chloroplast genomes of both Fortunella and Citrus species provide a mechanism of distinguishing the species for future studies. The study of chloroplast genes is of great significance in revealing the mechanism and metabolic regulation of plant photosynthesis. At the same time, the in-depth study of the chloroplast genome helps understand the mutual regulation between the nuclear genome and the chloroplast genome, and the chloroplast as a semi-autonomous organelle is conducive to energy conversion research.
Supplementary Materials
Rearrangement and reversal,
Author Contributions
The collection of experimental materials was completed by T.W., R.-P.K., X.-L.L. and X.-H.W. Data analysis by T.W.; preparations for drafting the manuscript and diagrams were completed by T.W., X.-Z.C. and K.-M.L. The revision and manuscript editing were completed by T.W., K.-M.L., X.-Z.C., V.O.W. and G.-W.H. Proofreading of the English manuscript was completed by T.W., K.-M.L., X.-Z.C., V.O.W. and G.-W.H. Resources were provided by K.-M.L., G.-W.H. and X.-Z.C. The funds were provided by K.-M.L. and G.-W.H. All authors have read and agreed to the published version of the manuscript.
Funding
This work was supported by grants from the National Science and Technology Basic Resources Survey Project (2019FY101800), Investigation of Forest Tree Germplasm Resources in Hunan Province (Xiangcainongzhi (2015) 91), and International Partnership Program of Chinese Academy of Sciences (151853KYSB20190027).
Institutional Review Board Statement
Not applicable.
Informed Consent Statement
Not applicable.
Data Availability Statement
Not applicable.
Acknowledgments
We sincerely thank Hunan forestry bureau, Hainan forestry bureau, Chaling County Agriculture Bureau, Guidong County natural resources bureau, Rucheng County Forestry Bureau and other units, and relevant personnel for their strong support and help. Thanks to Ying Tan for revising the English version of this article and making valuable comments on the manuscript. Thanks to Yi-Yan Cong, Jing Tian, Liu Zhou, Qin-You Zhang, Cun-Zhong Huang and others for their assistance in field investigation and material collection. Thanks to Shuai Peng, Jia-Xin Yang, Xiang Dong, and Shi-Xiong Ding from Wuhan Botanical Garden, Chinese Academy of Sciences for their guidance in experiment and data analysis. Thanks to Hui-Juan Shu for her help in data processing. We also thank Hunan Province Key Laboratory of Crop Sterile Germplasm Resource Innovation and Application for their help.
Conflicts of Interest
The authors declare no conflict of interest.
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Figures and Tables
Enlarge this image.Figure 1. Gene map of the Fortunella venosa chloroplast genome. Genes drawn inside the circle are transcribed clockwise, and those outside are transcribed counterclockwise. Genes belonging to different functional groups are color coded. The darker gray color in the inner circle corresponds to the GC content, and the lighter gray color corresponds to the AT content.
Enlarge this image.Figure 2. Comparison of repeated sequences in ten Rutaceae chloroplast genomes (Type and abundance of long repetitive sequences). Note: C represents complementary repeats, F represents Forward repeats, R represents reverse repeats, P represents palindromic repeats.
Enlarge this image.Figure 3. Analysis of SSRs in Fortunella venosa chloroplast genomes. (Type and abundance of SSRs).
Enlarge this image.Figure 4. Analysis of SSRs in ten Rutaceae chloroplast genomes. (The number of different SSRs detected in 10 genomes). The Mono-, Di-, Tri-, Tetra-, Penta-, and Hexa- represents the nucleotide motifs of the SSRs present in the 10 species genome.
Enlarge this image.Figure 5. Analysis of SSRs in Fortunella venosa chloroplast genomes (frequency of identified SSRs in LSC, SSC, and IR regions). The Mono-, Di-, Tri-, Tetra-, Penta-, and Hexa- represents the nucleotide motifs of the SSRs present in the F. venosa genome.
Enlarge this image.Figure 6. DNA sequence comparison of the ten species of Rutaceae. VISTA-based identity plot showing sequence identity among ten Rutaceae species using Fortunella venosa as a reference.
Enlarge this image.Figure 7. Comparison of border distances between adjacent genes and junctions of LSC, SSC, and two IR regions among chloroplast genomes of the ten Rutaceae species. Boxes above or below the main line indicate the adjacent border genes. The figure is not to scale with regard to sequence length and only shows relative changes at or near IR/SC borders.
Enlarge this image.Figure 8. Nucleotide diversity of different regions of Rutaceae chloroplast genome (horizontal axis represents the midpoint of the window, and vertical axis represents the nucleotide diversity of the window Pi).
Enlarge this image.Figure 9. Phylogenetic tree generation using 76 CDS common in 27 Rutaceae species. Melia azedarach were used as outgroups. The numbers above the branch represent bootstrap support value for BI/ML/PhyML methods, where the asterisk signifies maximum support value of 100 in IQ and 1 BI. Blank branches signify 100% support value.
Enlarge this image.Figure 10. Phylogenetic trees based on BI of Rutaceae species based on whole chloroplast genome sequences, with one species from family Melia used as outgroup. The Bayesian inference (BI) tree with posterior probabilities values on the branches.
Enlarge this image.Figure 11. Main morphological forms of Fortunella japonica and F. venosa. (A). F. japonica; (B). F. venosa. 1. Plants and habitats; 2. Branches; 3. Leaves; 4. Flowers; 5. Fruit branches; 6. Fruits and fruits cross-cut.
Summary of complete chloroplast genomes for ten Rutaceae species.
| Species/Taxa | Fortunella venosa | Fortunella japonica | Citrus aurantifolia | Citrus aurantium | Citrus hongheensis | Citrus cavaleriei | Citrus limon | Citrus maxima | Citrus medica | Citrus sinensis | |
|---|---|---|---|---|---|---|---|---|---|---|---|
| Accession |
MZ457935 | MN495932 | KJ865401 | MT702983 | MT880607 | MT880606 | MT880608 | MN782007 | MT106673 | DQ864733 | |
| Total Number of Genes | 134 | 135 | 138 | 132 | 135 | 135 | 135 | 125 | 134 | 134 | |
| Genome | Total GC content (%) | 38.4 | 38.4 | 38.4 | 38.5 | 38.5 | 38.5 | 38.5 | 38.5 | 38.4 | 38.5 |
| Total Length(bp) | 160,265 | 160,229 | 159,893 | 160,140 | 160,275 | 160,996 | 160,141 | 160,186 | 160,031 | 160,129 | |
| CDS | number | 89 | 90 | 93 | 87 | 89 | 89 | 89 | 88 | 89 | 89 |
| Length(bp) | 79,983 | 80,568 | 81,363 | 80,097 | 79,509 | 80,097 | 79,509 | 79,971 | 80,370 | 79,971 | |
| GC (%) | 38.9 | 38.8 | 39 | 38.8 | 38.9 | 38.8 | 38.9 | 38.8 | 39 | 38.8 | |
| tRNA | number | 37 | 37 | 37 | 37 | 37 | 37 | 37 | 37 | 37 | 37 |
| Length(bp) | 2790 | 2792 | 2802 | 2793 | 2792 | 2792 | 2792 | 2800 | 2802 | 2802 | |
| GC (%) | 53.3 | 53.3 | 53.3 | 53.3 | 53.4 | 53.3 | 53.4 | 53.3 | 53.2 | 53.3 | |
| rRNA | number | 8 | 8 | 8 | 8 | 8 | 8 | 8 | 8 | 8 | 8 |
| Length(bp) | 9044 | 9048 | 9048 | 9044 | 9048 | 9050 | 9048 | 9046 | 9048 | 9048 | |
| GC (%) | 55.7 | 55.7 | 55.7 | 55.7 | 55.7 | 55.7 | 55.7 | 55.7 | 55.7 | 55.7 | |
Genes present and functional gene category in F. venosa chloroplast genome.
| Category | Group of Genes | Name of Genes |
|---|---|---|
| Self-replication | Ribosomal protein (LSU) | * rpl2, rpl14, * rpl16, rpl20, rpl22, rpl23, rpl33, rpl32, rpl36 |
| Ribosomal proteins (SSU) | rps2, rps3, rps4, rps7, rps8, rps11, * rps12, rps14, rps15, * rps16, rps18, rps19 | |
| DNA-dependent RNA polymerase | rpoA, rpoB, * rpoC1, rpoC2 | |
| rRNA genes | rrn4.5, rrn5, rrn16, rrn23 | |
| tRNA genes | * trnA-UGC, trnC-GCA, trnD-GUC, trnE-UUC, trnF-GAA, trnfM-CAU, | |
| trnG-GCC, * trnG-UCC, trnH-GUG, trnI-CAU, * trnI-GAU, trnK-UUU, | ||
| trnL-CAA, * trnL-UAA, trnL-UAG, trnM-CAU, trnN-GUU, trnP-UGG, | ||
| trnQ-UUG, trnR-ACG, trnR-UCU, trnS-GCU, trnS-GGA, trnS-UGA, trnT-GGU, | ||
| trnT-UGU, trnV-GAC, * trnV-UAC, trnW-CCA, trnY-GUA | ||
| Photosynthesis | Photosystem I | psaA, psaB, psaC, psaI, psaJ |
| Photosystem II | psbA, psbB, psbC, psbD, psbE, psbF, psbH, psbI, psbJ, psbK, psbL, psbM, | |
| psbN, psbT, psbZ | ||
| NADPH dehydrogenase | * ndhA, * ndhB, ndhC, ndhD, ndhE, ndhF, ndhG, ndhH, ndhI, ndhJ, ndhK | |
| ATP synthase | atpA, atpB, atpE, * atpF, atpH, atpI | |
| Cytochrome c-type synthesis | petA, * petB, * petD, petG, petL, petN | |
| Rubisco | rbcL | |
| Other genes | Maturase | matK |
| Cytochrome c-type synthesis | ccsA | |
| Carbon metabolism | cemA | |
| Fatty acid synthesis | accD | |
| Transfer initiation factor | infA | |
| Proteolysis | ** clpP | |
| unknown | Conserved open reading frames | ycf1, ycf2, ** ycf3, ycf4, ycf68, ycf15 |
* Genes have one intron. ** Genes have two introns.
Table 3Analysis of relative synonymous codon usage (RSCU)in Fortunella venosa chloroplast genome.
| Codon | Count | RSCU | Codon | Count | RSCU | Codon | Count | RSCU | Codon | Count | RSCU |
|---|---|---|---|---|---|---|---|---|---|---|---|
| UUU(F) | 975 | 1.27 | UCU(S) | 551 | 1.6 | UAU(Y) | 774 | 1.59 | UGU(C) | 233 | 1.48 |
| UUC(F) | 556 | 0.73 | UCC(S) | 343 | 1 | UAC(Y) | 197 | 0.41 | UGC(C) | 82 | 0.52 |
| UUA(L) | 832 | 1.76 | UCA(S) | 384 | 1.12 | UAA (*) | 51 | 1.72 | UGA (*) | 15 | 0.51 |
| UUG(L) | 586 | 1.24 | UCG(S) | 243 | 0.71 | UAG (*) | 23 | 0.78 | UGG(W) | 455 | 1 |
| CUU(L) | 590 | 1.24 | CCU(P) | 407 | 1.45 | CAU(H) | 461 | 1.43 | CGU(R) | 321 | 1.15 |
| CUC(L) | 222 | 0.47 | CCC(P) | 243 | 0.87 | CAC(H) | 183 | 0.57 | CGC(R) | 131 | 0.47 |
| CUA(L) | 395 | 0.83 | CCA(P) | 321 | 1.14 | CAA(Q) | 706 | 1.53 | CGA(R) | 386 | 1.38 |
| CUG(L) | 219 | 0.46 | CCG(P) | 152 | 0.54 | CAG(Q) | 215 | 0.47 | CGG(R) | 153 | 0.55 |
| AUU(I) | 1071 | 1.47 | ACU(T) | 528 | 1.56 | AAU(N) | 961 | 1.51 | AGU(S) | 398 | 1.16 |
| AUC(I) | 461 | 0.63 | ACC(T) | 264 | 0.78 | AAC(N) | 313 | 0.49 | AGC(S) | 147 | 0.43 |
| AUA(I) | 651 | 0.89 | ACA(T) | 391 | 1.16 | AAA(K) | 1041 | 1.47 | AGA(R) | 493 | 1.77 |
| AUG(M) | 632 | 1 | ACG(T) | 168 | 0.5 | AAG(K) | 371 | 0.53 | AGG(R) | 189 | 0.68 |
| GUU(V) | 527 | 1.45 | GCU(A) | 626 | 1.71 | GAU(D) | 848 | 1.58 | GGU(G) | 552 | 1.2 |
| GUC(V) | 181 | 0.5 | GCC(A) | 248 | 0.68 | GAC(D) | 226 | 0.42 | GGC(G) | 198 | 0.43 |
| GUA(V) | 536 | 1.48 | GCA(A) | 390 | 1.07 | GAA(E) | 1016 | 1.47 | GGA(G) | 698 | 1.52 |
| GUG(V) | 209 | 0.58 | GCG(A) | 200 | 0.55 | GAG(E) | 366 | 0.53 | GGG(G) | 394 | 0.86 |
* represents the stop codons.
Table 4Positively selected sites detected in the chloroplast genome of Rutaceae based of Bayes empirical Bayes (BEB) method.
| Gene Name | M8 | |
|---|---|---|
| Selected Sites | Pr (w > 1) | |
| ndhB | 4084R | 0.990 ** |
| ndhI | 6657M | 1.000 ** |
| 6658S | 1.000 ** | |
| psbZ | 11729L | 0.990 ** |
| rpoC2 | 16680Y | 0.999 ** |
| 16682C | 0.999 ** | |
| 16683I | 0.999 ** | |
| 16703T | 0.998 ** | |
| 16705R | 0.991 ** | |
| 16706A | 0.998 ** | |
| 16714G | 0.998 ** | |
| 16725Y | 0.997 ** | |
| rps18 | 17492N | 0.999 ** |
| rps19 | 17529A | 0.990 ** |
| rps2 | 17766Y | 1.000 ** |
** p < 0.01.
© 2021 by the authors.
Details
1 College of Life Sciences, Hunan Normal University, Changsha 410081, China;
2 Key Laboratory of Plant Germplasm Enhancement and Specialty Agriculture, Wuhan Botanical Garden, Chinese Academy of Sciences, Wuhan 430074, China;