Next Article in Journal
Distribution of Exonic Variants in Glycogen Synthesis and Catabolism Genes in Late Onset Pompe Disease (LOPD)
Next Article in Special Issue
Genome-Wide Identification of Candidate Genes Associated with Heat Stress in Mulberry (Morus alba L.)
Previous Article in Journal
Apolipoprotein E (APOE) Haplotypes in Healthy Subjects from Worldwide Macroareas: A Population Genetics Perspective for Cardiovascular Disease, Neurodegeneration, and Dementia
Previous Article in Special Issue
Sucrose Transporter StSUT2 Affects Potato Plants Growth, Flowering Time, and Tuber Yield
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Characteristics, Comparative Analysis, and Phylogenetic Relationships of Chloroplast Genomes of Cultivars and Wild Relatives of Eggplant (Solanum melongena)

1
Institute of Vegetable Research, Guangxi Academy of Agricultural Sciences, Nanning 530003, China
2
Habin Academy of Agricultural Sciences, Harbin 150008, China
3
College of Horticulture, South China Agricultural University, Guangzhou 510642, China
4
Engineering Research Center for Horticultural Crop Germplasm Creation and New Variety Breeding, Ministry of Education, Changsha 410128, China
*
Author to whom correspondence should be addressed.
These authors contributed equally to this work.
Curr. Issues Mol. Biol. 2023, 45(4), 2832-2846; https://doi.org/10.3390/cimb45040185
Submission received: 27 February 2023 / Revised: 22 March 2023 / Accepted: 28 March 2023 / Published: 1 April 2023
(This article belongs to the Special Issue Functional Genomics and Comparative Genomics Analysis in Plants)

Abstract

:
The eggplant (Solanum melongena) is a popular vegetable around the world. However, the origin and evolution of eggplant has long been considered complex and unclear, which has become the barrier to improvements in eggplant breeding. Sequencing and comparative analyses of 13 complete chloroplast (cp) genomes of seven Solanum species were performed. Genome sizes were between 154,942 and 156,004 bp, the smallest genome was from S. torvum and the largest from S. macrocapon. Thirteen cp genomes showed highly conserved sequences and GC contents, particularly at the subgenus level. All genes in the 13 genomes were annotated. The cp genomes in this study comprised 130 genes (i.e., 80 protein-coding genes, 8 rRNA genes, and 42 tRNA genes), apart from S. sisymbriifolium, which had 129 (79 protein-coding genes, 8 rRNA genes, and 42 tRNA genes.). The rps16 was absent from the cp genome of S. sisymbriifolium, resulting in a nonsense mutation. Twelve hotspot regions of the cp genome were identified, which showed a series of sequence variations and differed significantly in the inverted repeat/single-copy boundary regions. Furthermore, phylogenetic analysis was conducted using 46 cp genomic sequences to determine interspecific genetic and phylogenetic relationships in Solanum species. All species formed two branches, one of which contained all cultivars of the subgenus Leptostemonum. The cp genome data and phylogenetic analysis provides molecular evidence revealing the origin and evolutionary relationships of S. melongena and its wild relatives. Our findings suggest precise intra- and interspecies relatedness within the subgenus Leptostemonum, which has positive implications for work on improvements in eggplant breeding, particularly in producing heterosis, expanding the source of species variation, and breeding new varieties.

1. Introduction

Solanaceae plants are medium-sized angiosperms, the largest group of vegetable crops, and the third-largest group of economic plants [1]. The taxa of Solanaceae plants are abundant and diverse, containing 90 genera and including 3000–4000 species. This genus includes many important crop species, raw industrial materials, and certain plant models used in research. Therefore, Solanaceae plants are often regarded as important research materials [2]. Due to the considerable economic and research significance of the Solanaceae family, taxonomic and phylogenetic studies on this family have received much attention [3,4,5]. The genus Solanum accounts for half of the Solanaceae family, with approximately 2000 species that are mostly distributed in tropical and subtropical regions [1,6]. Wild species of Solanum are found in Southeast Asia, Africa, South America, and southern China. Solanum is large and strongly monophyletic, which makes taxonomic analysis of this group difficult [1,6,7]. Isolation and differentiation of species are the basic processes involved in new species formation. However, there is no strict reproductive isolation in Solanum plants. Their pollen exhibits a certain affinity, and many species have the potential to cross-hybridize [8,9]. Further, during the process of artificial hybridization and variety selection, many intermediate hybrids are also produced. Therefore, there are many controversies regarding the identities of local breeds, intermediate hybrids, and related wild species, as well as their evolutionary relationships and classifications [10,11].
The chloroplast (cp) genome is the most prominent marker of green cells in leaves and occupies the second-largest area in mesophyll cells, besides the vacuole [12]. Plastids have diverse functions, including starch biosynthesis, nitrogen metabolism, sulfate reduction, fatty acid synthesis, and DNA and RNA synthesis [13]. The cp genomes of most terrestrial plants are circular chromosomes consisting of four parts: two inverted repeats (IR), a large single copy (LSC) region, and a small single copy (SSC) region [14]. In photosynthetic organisms, the genome length is 115–165 kb [14,15]. Most cp DNAs shows monoparental inheritance, with relatively few recombination and intraspecific mutation events. The phylogenetic tree can be constructed without relying on any other data, similar to investigations of the evolutionary history of plants [3,16]. The taxonomic and phylogenetic studies of the subgenus Leptostemonum have always been an area of great interest, and the species-level taxonomy of three cultivated aubergine species, the brinjal eggplant (Solanum melongena), the scarlet eggplant (S. aethiopicum) and the gboma eggplant (S. macrocarpon), which has long been considered complex [1,3,17,18,19].
With the rapid development of sequencing technology, bringing genomics into the large-scale, low-cost, high-throughput sequencing era has greatly promoted the research progress of the chloroplast genome in plants. The cp genomes of more than 1800 plants and 94 Solanum plants are available from the National Centre for Biotechnology Information (NCBI) Organelle Genome Resources. These cp genomes are widely used in comparative genomics and research on plant phylogeny and classification. The subgenus Leptostemonum contains a wide range of plants; however, its classification remains unclear.
In this study, we utilized Illumina HiSeq and analyzed 13 cp genomes of 7 Leptostemonum species. A phylogenetic tree was constructed to assess the interspecific differences and structural patterns of the cp genomes within the genus. Based on these results, the origin and evolutionary relationships of Leptostemonum species were analyzed to reveal the genetic resources available for use in eggplant breeding.

2. Materials and Methods

2.1. Materials

Seven Leptostemonum species were collected from Africa, South America, Southeast Asia, and south China and transplanted in the Resource Garden of Vegetable Research Institute, Guangxi Academy of Agricultural Science, Nanning, China. Thirteen samples of the seven Leptostemonum species were obtained for cp genome sequencing.

2.2. DNA Extraction, Genome Sequencing, and Assembly

Approximately 5 g of fresh leaves was harvested for cp DNA isolation using an improved extraction method [20]. After DNA isolation, 1 μg of purified DNA was fragmented to construct short-insert libraries (insert size 430 bp) according to the manufacturer’s instructions (Illumina), followed by sequencing on the Illumina HiSeq 4000 [21]. Prior to assembly, low-quality raw reads with adaptors showing a quality score below 20 or containing N were filtered. The cp genome was reconstructed using a combination of de novo and reference-guided assemblies, and the following three steps were used to assemble the cp genomes [22]. First, filtered reads were assembled into contigs using SOAPdenovo2.04 [23]. Second, contigs were aligned to the reference genome of eggplant using BLAST and aligned contigs (≥80% similarity and query coverage) were ordered according to the reference genome. Third, clean reads were mapped to the assembled draft cp genome to correct wrong bases, and the majority of gaps were filled through local assembly by GapCloser. Insertions/deletions (InDels) of 13 cp genomes were assessed. The physical map of the circular cp genome of the 13 samples were drawn by Organellar GenomeDRAW (http://ogdraw.mpimp-golm.mpg.de/cgi-bin/ogdraw.pl (accessed on 26 March 2019)).

2.3. Genome Annotation

The cp protein-coding, transfer RNA (tRNA), and ribosomal RNA (rRNA) genes were annotated using the online DOGMA tool using default parameters [24]. A whole cp genome BLAST search (E-value 10−5, minimal alignment length percentage 40%) was performed against five databases: the Kyoto Encyclopedia of Genes and Genomes, Clusters of Orthologous Groups, Non-Redundant Protein Database databases, Swiss-Pro, and ChloroplastDB [25,26,27,28]. The circular chloroplast genome map of Solanum was drawn using Organellar Genome DRAW v1.2. The cp genomes of the 13 samples were compared using VISTA [29]. Genome, protein-coding gene, intron, and spacer sequence divergences were evaluated using DnaSP 5.10 after alignment [30]. The genomic sequences were aligned using MAFFT v5 and adjusted manually where necessary [31]. For protein-coding gene sequences, introns, and spacers, every gene or fragment was edited using the ClustalW multiple alignment option within BioEdit v7.0.9.0 [32].

2.4. IR Expansion and Contraction

IR expansion and contraction was determined using the site https://irscope.shinyapps.io/irapp/analysis (accessed on 26 March 2019).

2.5. Codon Usage

The EMBOSS-6.6.0 cusp module was used to analyze the coding sequence area [33].

2.6. Comparative Genome Analysis

MUMmer (3.23) and BLAST were used to conduct global and local alignments between the sample and reference genomes, determining potential single-nucleotide polymorphisms (SNPs). Subsequently, SNPs were filtered out in the repeat regions as detected using the software BLAST, Repeat Masker, and TRF. Finally, SNPs were annotated based on the position and interaction between genes. LASTZ software (Release 1.04.15) was used for global alignment between each sample sequence and the reference genome. The alignment result was corrected using axt_correction, axtSort, and axtBest to determine potential InDels with lengths less than 50 bp. Finally, BWA and SAM tools were used to map the reads to InDel sequences and filter out unreliable InDels [34,35].

2.7. Repeat and SSR Analysis

Repeats and simple sequence repeats (SSRs) were determined using the website https://bibiserv.cebitec.uni-bielefeld.de/reputer (accessed on 26 March 2019). Analysis, parameters of minimum repeat length 30 bp and hamming distance 3. SSRs are found widely in the genome. Generally, they are composed of 1–6 bp repeat sequences with a low degree, with 2–3 nucleotides as repeat units such as (GA)n, (AC)n, and (GAA)n. The MIcroSAtellite identification tool (MISA) was used to detect microsatellite loci. Definition (unit_size, min_repeats): 1–10, 2–6, 3–5, 4–5, 5–5, 6–5. The minimum distance between the two SSRs was set to 100 bp. Parameter meaning: 1 base repeat 10 time or more; 2 bases repeat 6 times or more; 3 base repeats 5 times or more; repeat 4 bases for 5 times or more; 5 base repeats 5 times or more; a sequence with 6 base repeats of 5 or more times is considered a microsatellite sequence. Concurrently, when the distance between the two microsatellites was less than 100 bp, they formed a composite microsatellite.

2.8. Identification of the Most Variable Regions

After aligning the sequence, the Pi value was calculated using a sliding window (window 300 bp, step size 200 bp). A window with a Pi value greater than or equal to 0.01 was selected and primers were designed at both ends of the window sequence.

2.9. Phylogenetic Analysis

ClustalW was used to align the cp DNA sequences with default parameters, and the alignment was checked manually [32]. Maximum-likelihood (ML) methods were used for gnome-wide phylogenetic analyses using PhyML 3.0 [36]. The nucleotide substitution model was selected using jModelTest 2.1.10 and Smart Model Selection in PhyML 3.0 [37]. TheGeneral Time Reversible model (GTR model) with default parameter was selected for ML analyses with 1000 bootstrap replicates to calculate the bootstrap values of the topology. The results were treated with iTOL 3.4.3 [38].

3. Results

3.1. Characteristics of Cp Genomes of Leptostemonum Species

The 13 cp DNA samples of seven Leptostemonum species were sequenced using the Illumina HiSeq 4000 (San Diego, CA, USA): 5120–10,165 Mb raw paired-end reads were generated, with an average organelle depth ranging from 636 to 3639 and a Q30 of 92.27–98.72% (Table 1). Data from the 13 complete cp genomes were submitted to NCBI (Table 2). All 13 complete cp genomes consisted of a single circular double-stranded DNA molecule with a classic four-part structure, as observed in most angiosperms, including the LSC, SSC, IRA, and IRB regions (Figure 1).
No prominent sequence inversions of genomic rearrangements were observed. Among the 13 cp genomes, the assembled length was 154,942–156,004 bp, LSC length was 85,646–86,542 bp, SSC length was 18,435–18,602 bp, IRA and IRB lengths were consistently 25,420–25,461 bp, and total GC content was 37.76–37.81% (Table 3). These results indicated that there were no significant differences in the genome size, GC content, LSC length, SSC length, IRA length, and IRB length of the seven Solanum species. These results also demonstrated the characteristics of Solanum cp genomes, which included short-length and conserved genes and genomic structures. Most Solanum species are shown together in a pie chart, as their numbers, orders, and gene names were the same (130 genes, i.e., 80 protein-coding genes, 8 rRNA genes, and 42 tRNA genes) (Figure 1A, Table 3 and Table 4). NN15 (Solanum sisymbriifolium) was annotated in a separate pie chart as it only contained 129 genes (79 protein-coding genes, 8 rRNA genes, and 42 tRNA genes), and the rps16 gene was divided into two parts in NN15 that resulted in a nonsense mutation (Figure 1B, Table 3 and Table 4).
The expansion and contraction of IR regions can cause differences in the genome size. Therefore, we compared the differences in the adjacent regions and their adjacent genes between the 13 cp genomes (Figure 2). The genes rps19, ndhF, ycf1, and trnH were located in the junctions of LSC/IRb, IRb/SSC, SSC/IRa, and IRa/LSC, respectively. Compared to the relatively conserved location of trnH and rps19, the SSC/IR boundary regions were more variable. trnH in 11 cp genomes was between the SSC/IRb region with 2215 and 2216 bp in the SSC region and 7 and 8 bp in the IRb region, respectively. Specifically, trnH of NN9 (S. wrightii) was in the SSC region and trnH of NN15 extended 7 bp into the SSC region. In all 13 cp genomes, ycf1 was located inthe SSC/IRa region, with 140–160 bp extending into the IRa region. No gene rearrangement or inversion events were observed.
Analysis of relative synonymous codon usage (RSCU) value showed that almost all amino acids had more than one synonymous codon, except for methionine and tryptophan (Figure 2A, Table S1). Nearly all protein-coding genes of the 13 cp genomes had a standard ATG start codon (RSCU = 1). Approximately half of the codons had an RSCU > 1, and most (29/31, 93.5%) ended with either A or T. RSCU = 1 of codon (TGG) indicated a balanced bias of the codon (Figure 2A). Leucine (2665, 10.82%) was the most common amino acid and cysteine (282, 1.14%) was the least common in the NN1 (S. aethiopicum) cp genome, as well as the other cp genomes (Figure 2B, Table S1).

3.2. Comparative Genome Analysis

Comparative analysis was performed on 13 cp genomes using KU682719 (S. melongena plastid) as a reference. The genomes exhibited a high degree of sequence synteny, suggesting a highly conserved evolutionary pattern. However, nucleotide substitutions, InDels, and length variations were observed among the species (Table 5). The number of SNPs in S. melongena (NN5, NN7, NN8) was smallest (9–10), with insignificant positions. The sample with the largest number of SNPs was S. sisymbriifolium (NN15; 259), which showed the mostdistinct genetic relationship with KU682719. Solanum aethiopicum, S. indicum, S. macrocapon, S. wrightii, and S. torvum showed moderate differences, with the first three closer to S. melongena and the latter two closer to S. sisymbriifolium; their ratios of intergenic SNPs were approximately 70%. Moreover, large differences were observed in the number, type, and positions of InDels among species, but little difference was observed among intraspecific species (Table 5). Solanum melongena contained the smallest number of InDels, ranging from 53 to 55. Eighty-five percent of the InDel loci were in the intergenic region, and more than 84% were insertions. The sample with the largest number of InDel loci was S. sisymbriifolium (NN15), which had the most distant relationship with KU682719 and largest number of differential InDel loci (up to 126). Ninety-six percent of the InDel loci were in the intergenic region, and 60.3% were insertions. Solanum aethiopicum, S. indicum, S. macrocapon, S. wrightii, and S. torvum were in the middle of the extremes. However, most InDels were in the intergenic region and the ratio of insertions in these five Solanum species was closer to that of S. sisymbriifolium than to KU682719. Most SNPs and InDels were found in the LSC and SSC regions, not in the IR regions (Figure 2C,D). This demonstrated the conserved characteristics of the IR region. At the whole-genome level, there were limited subgenus differences in the InDels of the Solanum species. The most common subgenus difference was single-base insertions or deletions. The number of single-base InDels varied in S. melongena compared to that in other species, with only 2–3 in S. melongena and 38 in S. sisymbriifolium. In other close relatives of Solanum, single-base InDel mutations accounted for more than 30% of the total number of InDels. Furthermore, the average length of the InDels in the 13 samples was 6–9 bp (Figure 2E).

3.3. Repeat and SSR Analysis

Four types of repeat sequences were determined: direct (forward), inverted (palindromic), complement, and reverse repeats. The numbers and distributions of the repeats in the 13 cp DNA were similar and conserved among these types. In total, 31–50 repeats were identified, including 16–30 forward repeats, 14–21 palindromic repeats, 0–9 reverse repeats, and 0–2 complement repeats (Figure 3A). Of the 13 cp genomes, forward and palindrome were the most abundant repeat types, with mostly consistent amounts. Reverse repeats were detected in S. aethiopium (NN1, NN2, NN3, NN4), S. macrocapon (NN6), S. wrightii (NN9), S. indicum (NN13), and S. sisymbriifolium (NN15) in amounts of 3, 3, 3, 3, 9, 5, 1, and 3, respectively (Figure 3B). Two complements were detected in NN1, NN2, NN3, and NN4 (Figure 3B). Most repeats were 30–39 bp, followed by 40–49 bp and over 50 bp (Figure 3B).
In the 13 cp genomes, 136–146 SSRs were identified: 107–121 mononucleotides (83.45%), 6–9 dinucleotides (5.52%), 4–7 trinucleotides (4.83%), 8–12 tetranucleotides (5.52%), and 1–3 pentanucleotides (0.69%) (Figure 3C,D). In NN1, 97.5% of the mononucleotides exhibited A or T types and 2.5% exhibited C or G types and all dinucleotides were composed of AT/TA (Figure 3E,F).

3.4. Identification of the Most Variable Regions

To determine the levels of sequence divergence, we calculated the nucleotide variability (Pi) values and investigated the levels of sequence divergence among genera in 46 Solanum genomes, including 13 Solanum cp genomes in this study and 33 other Solanum cp genomes downloaded from NCBI (Table S3). The Pi values within 300 bp across the 46 cp genomes ranged from 0 to 0.0695, with a mean value of 0.00351 (Figure 4A), indicating high similarity among species. However, hypervariable loci were also observed, including rps16-CDS1_trnQ-UUG, atpH_atpI, ccsA, petA-psbJ, psbA, rpl32_trnL-UAG, ndhD, atpI, psbJ, rpl32, rps15_ycf1-D2, and rbcL_psaI (Pi > 0.014). The rps15_ycf1-D2, ndhD, ccsA, rpl32_trnL-UAG, and rpl32 loci were present in the SSC region; psbJ, petA_psbJ, rbcL_psaI, atpI, atpH_atpI, rps16-CDS1_trnQ-UUG, and psbA were present in the LSC region (Figure 4A). These species shared the same order and orientation of syntenic blocks, indicating no rearrangement in gene organization and illustrating that cp genomes tend to be conserved and collineated, particularly in the same plant family.

3.5. Phylogenetic Analyses

To better determine the phylogenetic position of Leptostemonum and further clarify the evolutionary relationships within Leptostemonum, phylogenetic analysis was performed based on the 13 cp genomes in this study and an additional 33 Leptostemonum cp genomes downloaded from NCBI (Table S3). A phylogenetic tree was constructed using the maximum-likelihood (ML) method. All these species formed two branches. One branch mainly consisted of the most economically important cultivated aubergine worldwide. The three cultivated aubergine species, the brinjal eggplant (S. melongena), the scarlet eggplant (S. aethiopicum) and the gboma eggplant (S. macrocarpon), are closely related phylogenetically, and form a part of the subgenus Leptostemonum. Interestingly, high intraspecific genetic diversity was found among the three cultivated aubergine species, which formed different clades with their wild progenitors (Figure 4B). The African S. anguivi and the Asian S. violaceum are the sister species of S. aethiopicum. Moreover, three samples of S. melongena (NN5, NN7 and NN8) formed a phylogenetic tree with S. melongena (MH283708), but did not occur on the phylogenetic tree near S. melongena var. insanum (MH283711) (Figure 4B). The second branch consisted of S. torvum, S. sisymbriifolium, and S. wrightii (Figure 4B). Solanum torvum, S. sisymbriifolium, and S. wrightii are close relatives of aubergine, and are native to South America, Central America, and Africa. They are phylogenetically distant from the cultivated S. melongena.

4. Discussion

The Solanum subgenus Leptostemonum is the largest clade in genus Solanum and it contains three domesticated species (Solanum melongena, S. aethiopicum and S. macrocarpon). The taxonomy relationships and domestication of the cultivars and wild relatives is puzzling for breeders due to the large number of related species (Taher et al., 2017). In this study, we sequenced the cp genomes of seven subgenus Leptostemonum species, including three samples of S. melongena, four samples of S. aethiopicum, one sample of S. macrocarpon and five samples of their wild relatives, combining with 33 cp genome data of Leptostemonum species downloaded from NCBI to perform comparative genome and phylogenetic analysis. The results reveal precise intra- and interspecies relatedness within the subgenus Leptostemonum.
The cp genome of Solanum melongena, S. macrocarpon, S. anguivi, S. aethiopicum, S. violaceum, S. torvum and S. sisymbriifolium has been determined [18,39,40,41]. The cp genome of S. wrightii was the first report here. Chloroplast genomes have been successfully used in numerous phylogenetic studies of the subgenus Leptostemonum because of their high accuracy and resolution [18]. These analyses of cp genomes enhanced our understanding of the evolutionary trends and phylogenetic implications of Leptostemonum, and provide new molecular evidence for analyzing the genetic relationship of Leptostemonum species.
Interspecific hybridization is an effective way to create heterosis and species variation. Both S. aethiopicum and S. macrocarpon were partially interfertile to Solanum melongena [42,43]. Solanum aethiopicum exhibited higher cross-compatibility than S. macrocarpon when crossed with S. melongena [42,44]. This study strongly suggests that the genetic distances between S. aethiopicum and S. macrocarpon are smaller than that between S. melongena, respectively. Solanum aethiopicum consists of four cultivar groups (Gilo, Shum, Kumba, and Aculeatum), which exhibit high genetic diversity [45]. In our study, the four samples of S. aethiopicum were originally from different regions, and NN1 was the sister group of NN3 and NN4, while NN2 formed a split group. NN1, NN3 and NN4 were originally from Africa, whereas NN2 was from South America. This result was also present in the samples of S. melongena, which was comprised of four groups (E-H) [7,46]. These results suggested the genetic diversity of the subgenus Leptostemonum signatures associated with selection and domestication. Solanum anguivi is the ancestor of Solanum aethiopicum, which was consistent with the phylogenetic relationship of the two spices in this study [7,47,48]. Solanum violaceum occurs across Asia, but its close relatives were the African S. anguivi and S. aethiopicum, which suggested S. violaceum were originally from Africa [18]. The 2 samples of S. violaceum (NN13 and NN14) occur on a different group from the S. anguivi (MH283724) in our study, implying that was misidentification of species or hybridization. Solanum torvum, S. sisymbriifolium, and S. wrightii are phylogenetically distant from the cultivated S. melongena. S. torvum can be crossed with S. melongena to produce sterility F1, whereas S. sisymbriifolium cannot be crossed with S. melongena due to the distant genetic relationship [49,50]. There is no report about the interspecific hybridization between S. wrightii and S. melongena.

5. Conclusions

The cp genome sequences of 13 samples from seven Solanum species were analyzed in this study. Although the genomic structure and size were highly conserved, inverted repeat/single-copy boundary regions and variation among species were still detected, representing wide phylogenetic diversity in the genus. The availability of these cp genomes provides genetic information for identifying the species structure, taxonomy, phylogeny, and evolution in Solanum, along with insights into the utilization of Solanum plants.

Supplementary Materials

The following supporting information can be downloaded at https://www.mdpi.com/article/10.3390/cimb45040185/s1.

Author Contributions

Methodology, Q.Y. and P.W.; Software, Q.Y. and P.W.; Validation, Q.Y., L.C. and M.W.; Formal analysis, G.G.; Resources, C.X.; Data curation, Q.Y., W.L. (Weiliu Li) and Y.J.; Writing—original draft, Q.Y. and D.L.; Writing—review & editing, Y.L. and Y.W.; Supervision, W.L. (Wenjia Li) and R.C.; Project administration, Y.W.; Funding acquisition, Y.W. All authors have read and agreed to the published version of the manuscript.

Funding

This research was supported by the Science and Technology Planning Project of Guangxi (GuikeAD22035947, GuikeAA22068088-2), and Scientific Research Fund of Hunan Provincial Education Department (19B271).

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

The data used in this study have been deposited in the NCBI Genbank with the accession number in Tables S1 and S3.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Vorontsova, M.S.; Stern, S.; Bohs, L.; Knapp, S. African spiny Solanum (subgenus Leptostemonum, Solanaceae): A thorny phylogenetic tangle. Bot. J. Linn. Soc. 2013, 173, 176–193. [Google Scholar] [CrossRef] [Green Version]
  2. Rinaldi, R.; Van Deynze, A.; Portis, E.; Rotino, G.L.; Toppino, L.; Hill, T.; Ashrafi, H.; Barchi, L.; Lanteri, S. New Insights on Eggplant/Tomato/Pepper Synteny and Identification of Eggplant and Pepper Orthologous QTL. Front. Plant Sci. 2016, 7, 1031. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  3. Särkinen, T.; Bohs, L.; Olmstead, R.G.; Knapp, S. A phylogenetic framework for evolutionary study of the nightshades (Solanaceae): A dated 1000-tip tree. BMC Evol. Biol. 2013, 13, 214. [Google Scholar] [CrossRef] [Green Version]
  4. Olmstead, R.G.; Bohs, L.; Migid, H.A.; Santiago-Valentin, E.; Garcia, V.F.; Collier, S.M. A molecular phylogeny of the Solanaceae. TAXON 2008, 57, 1159–1181. [Google Scholar] [CrossRef]
  5. Wu, F.; Tanksley, S.D. Chromosomal evolution in the plant family Solanaceae. BMC Genom. 2010, 11, 182. [Google Scholar] [CrossRef] [Green Version]
  6. Knapp, S.; Vorontsova, M.S.; Prohens, J. Wild Relatives of the Eggplant (Solanum melongena L.: Solanaceae): New Understanding of Species Names in a Complex Group. PLoS ONE 2013, 8, e57039. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  7. Taher, D.; Solberg, S.; Prohens, J.; Chou, Y.-Y.; Rakha, M.; Wu, T.-H. World Vegetable Center Eggplant Collection: Origin, Composition, Seed Dissemination and Utilization in Breeding. Front. Plant Sci. 2017, 8, 1484. [Google Scholar] [CrossRef] [Green Version]
  8. Plazas, M.; Vilanova, S.; Gramazio, P.; Rodríguez-Burruezo, A.; Fita, A.; Herraiz, F.J.; Ranil, R.; Fonseka, R.; Niran, L.; Fonseka, H.; et al. Interspecific Hybridization between Eggplant and Wild Relatives from Different Genepools. J. Am. Soc. Hortic. Sci. 2016, 141, 34–44. [Google Scholar] [CrossRef] [Green Version]
  9. Afful, N.T.; Nyadanu, D.; Akromah, R.; Amoatey, H.M.; Annor, C.; Diawouh, R.G. Evaluation of crossability studies between selected eggplant accessions with wild relatives S. torvum, S. anguivi and S. aethopicum (Shum group). J. Plant Breed. Crop Sci. 2018, 10, 1–12. [Google Scholar] [CrossRef] [Green Version]
  10. Karihaloo, J.L.; Gottlieb, L.D. Allozyme variation in the eggplant, Solanum melongena L. (Solanaceae). Theor. Appl. Genet. 1995, 90, 578–583. [Google Scholar] [CrossRef]
  11. Ranil, R.H.G.; Prohens, J.; Aubriot, X.; Niran, H.M.L.; Plazas, M.; Fonseka, R.M.; Vilanova, S.; Fonseka, H.H.; Gramazio, P.; Knapp, S. Solanum insanum L. (subgenus Leptostemonum Bitter, Solanaceae), the neglected wild progenitor of eggplant (S. melongena L.): A review of taxonomy, characteristics and uses aimed at its enhancement for improved eggplant breeding. Genet. Resour. Crop. Evol. 2016, 64, 1707–1722. [Google Scholar] [CrossRef] [Green Version]
  12. Lopez-Juez, E.; Pyke, K. Plastids unleashed: Their development and their integration in plant development. Int. J. Dev. Biol. 2005, 49, 557–577. [Google Scholar] [CrossRef] [Green Version]
  13. Peeters, N.M.; Hanson, M.R.; Halter, C.P. RNA editing in ribosome-less plastids of iojap maize. Curr. Genet. 2004, 45, 331–337. [Google Scholar] [CrossRef] [PubMed]
  14. Jansen, R.K.; Raubeson, L.A.; Boore, J.L.; DePamphilis, C.W.; Chumley, T.W.; Haberle, R.C.; Wyman, S.K.; Alverson, A.J.; Peery, R.; Herman, S.J.; et al. Methods for Obtaining and Analyzing Whole Chloroplast Genome Sequences. In Methods in Enzymology; Academic Press: Cambridge, MA, USA, 2005; Volume 395, pp. 348–384. [Google Scholar]
  15. Wang, R.-J.; Cheng, C.-L.; Chang, C.-C.; Wu, C.-L.; Su, T.-M.; Chaw, S.-M. Dynamics and evolution of the inverted repeat-large single copy junctions in the chloroplast genomes of monocots. BMC Evol. Biol. 2008, 8, 36. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  16. Li, X.; Gao, H.; Wang, Y.; Song, J.; Henry, R.; Wu, H.; Hu, Z.; Yao, H.; Luo, H.; Luo, K.; et al. Complete chloroplast genome sequence of Magnolia grandiflora and comparative analysis with related species. Sci. China Life Sci. 2013, 56, 189–198. [Google Scholar] [CrossRef] [Green Version]
  17. Knapp, S.; Sagona, E.; Carbonell, A.K.; Chiarini, F. A revision of the Solanum elaeagnifolium clade (Elaeagnifolium clade; subgenus Leptostemonum, Solanaceae). Phytokeys 2017, 84, 1–104. [Google Scholar] [CrossRef] [Green Version]
  18. Aubriot, X.; Knapp, S.; Syfert, M.M.; Poczai, P.; Buerki, S. Shedding new light on the origin and spread of the brinjal eggplant (Solanum melongena L.) and its wild relatives. Am. J. Bot. 2018, 105, 1175–1187. [Google Scholar] [CrossRef] [Green Version]
  19. Vilanova, S.; Manzur, J.P.; Prohens, J. Development and characterization of genomic simple sequence repeat markers in eggplant and their application to the study of diversity and relationships in a collection of different cultivar types and origins. Mol. Breed. 2011, 30, 647–660. [Google Scholar] [CrossRef]
  20. McPherson, H.; Van Der Merwe, M.; Delaney, S.K.; Edwards, M.A.; Henry, R.J.; McIntosh, E.; Rymer, P.D.; Milner, M.L.; Siow, J.; Rossetto, M. Capturing chloroplast variation for molecular ecology studies: A simple next generation sequencing approach applied to a rainforest tree. BMC Ecol. 2013, 13, 8. [Google Scholar] [CrossRef] [Green Version]
  21. Borgström, E.; Lundin, S.; Lundeberg, J. Large Scale Library Generation for High Throughput Sequencing. PLoS ONE 2011, 6, e19119. [Google Scholar] [CrossRef] [Green Version]
  22. Cronn, R.; Liston, A.; Parks, M.; Gernandt, D.S.; Shen, R.; Mockler, T. Multiplex sequencing of plant chloroplast genomes using Solexa sequencing-by-synthesis technology. Nucleic Acids Res. 2008, 36, e122. [Google Scholar] [CrossRef] [Green Version]
  23. Luo, R.; Liu, B.; Xie, Y.; Li, Z.; Huang, W.; Yuan, J.; He, G.; Chen, Y.; Pan, Q.; Liu, Y.; et al. SOAPdenovo2: An empirically improved memory-efficient short-read de novo assembler. Gigascience 2012, 1, 18. [Google Scholar] [CrossRef] [PubMed]
  24. Wyman, S.K.; Jansen, R.K.; Boore, J.L. Automatic annotation of organellar genomes with DOGMA. Bioinformatics 2004, 20, 3252–3255. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  25. Kanehisa, M.; Goto, S.; Kawashima, S.; Okuno, Y.; Hattori, M. The KEGG resource for deciphering the genome. Nucleic Acids Res. 2004, 32, D277–D280. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  26. Tatusov, R.L.; Fedorova, N.D.; Jackson, J.D.; Jacobs, A.R.; Kiryutin, B.; Koonin, E.V.; Krylov, D.M.; Mazumder, R.; Mekhedov, S.L.; Nikolskaya, A.N.; et al. The COG database: An updated version includes eukaryotes. BMC Bioinform. 2003, 4, 41. [Google Scholar] [CrossRef] [Green Version]
  27. Magrane, M. Consortium UniProt Knowledgebase: A hub of integrated protein data. Database: J. Biol. Databases Curation 2011, 2011, bar009. [Google Scholar] [CrossRef] [Green Version]
  28. Cui, L.; Makałowska, I. ChloroplastDB: The Chloroplast Genome Database. Nucleic Acids Res. 2006, 34, D692–D696. [Google Scholar] [CrossRef] [Green Version]
  29. Lohse, M.; Drechsel, O.; Bock, R. OrganellarGenomeDRAW (OGDRAW): A tool for the easy generation of high-quality custom graphical maps of plastid and mitochondrial genomes. Curr. Genet. 2007, 52, 267–274. [Google Scholar] [CrossRef]
  30. Ogban, F.; Udensi, O.; Inyang, G.; Osang, F. Diversity in Content and Analytical Algorithms of Biologist-Centric Software(s) for DNA and Sequence Data: Mega, DNASp, GenAlex and ARLEQUIN. Int. J. Nat. Appl. Sci. 2019, 12, 47–60. [Google Scholar]
  31. Katoh, K.; Kuma, K.I.; Toh, H.; Miyata, T. MAFFT version 5: Improvement in accuracy of multiple sequence alignment. Nucleic Acids Res. 2005, 33, 511–518. [Google Scholar] [CrossRef]
  32. Larkin, M.A.; Blackshields, G.; Brown, N.P.; Chenna, R.; McGettigan, P.A.; McWilliam, H.; Valentin, F.; Wallace, I.M.; Wilm, A.; Lopez, R.; et al. Clustal W and Clustal X version 2.0. Bioinformatics 2007, 23, 2947–2948. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  33. Rice, P.; Longden, I.; Bleasby, A. EMBOSS: The European Molecular Biology Open Software Suite. Trends Genet. 2000, 16, 276–277. [Google Scholar] [CrossRef] [PubMed]
  34. Li, H.; Durbin, R. Fast and accurate short read alignment with Burrows—Wheeler transform. Bioinformatics 2009, 25, 1754–1760. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  35. Li, H.; Handsaker, B.; Wysoker, A.; Fennell, T.; Ruan, J.; Homer, N.; Marth, G.; Abecasis, G.; Durbin, R. The Sequence Alignment/Map format and SAMtools. Bioinformatics 2009, 25, 2078–2079. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  36. Guindon, S.; Dufayard, J.-F.; Lefort, V.; Anisimova, M.; Hordijk, W.; Gascuel, O. New Algorithms and Methods to Estimate Maximum-Likelihood Phylogenies: Assessing the Performance of PhyML 3.0. Syst. Biol. 2010, 59, 307–321. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  37. Darriba, D.; Taboada, G.L.; Doallo, R.; Posada, D. jModelTest 2: More models, new heuristics and parallel computing. Nat. Methods 2012, 9, 772. [Google Scholar] [CrossRef] [Green Version]
  38. Letunic, I.; Bork, P. Interactive tree of life (iTOL) v3: An online tool for the display and annotation of phylogenetic and other trees. Nucleic Acids Res. 2016, 44, W242–W245. [Google Scholar] [CrossRef]
  39. Ding, Q.-X.; Liu, J.; Gao, L.-Z. The complete chloroplast genome of eggplant (Solanum melongena L.). Mitochondrial DNA Part B 2016, 1, 843–844. [Google Scholar] [CrossRef] [Green Version]
  40. Li, D.; Gan, G.; Li, W.; Li, W.; Jiang, Y.; Liang, X.; Yu, N.; Chen, R.; Wang, Y. Inheritance of Solanum chloroplast genomic DNA in interspecific hybrids. Mitochondrial DNA Part B 2021, 6, 351–357. [Google Scholar] [CrossRef]
  41. Yin, M.; Yu, Y.; Gong, Y.; Gui, M.; Li, Z.; Bao, R.; Cheng, J.; Du, G.; Wu, L. The complete chloroplast genome of Solanum sisymbriifolium (Solanaceae), the wild eggplant. Mitochondrial DNA Part B 2022, 7, 886–888. [Google Scholar] [CrossRef]
  42. Devi, C.P.; Munshi, A.D.; Behera, T.K.; Choudhary, H.; Gurung, B.V.; Saha, P. Cross compatibility in interspecific hybridization of eggplant, Solanum melongena, with its wild relatives. Sci. Hortic. 2015, 193, 353–358. [Google Scholar] [CrossRef]
  43. Prohens, J.; Plazas, M.; Raigón, M.D.; Seguí-Simarro, J.M.; Stommel, J.R.; Vilanova, S. Characterization of interspecific hybrids and first backcross generations from crosses between two cultivated eggplants (Solanum melongena and S. aethiopicum Kumba group) and implications for eggplant breeding. Euphytica 2012, 186, 517–538. [Google Scholar] [CrossRef]
  44. Bletsos, F.; Roupakias, D.; Tsaktsira, M.; Scaltsoyjannes, A. Production and characterization of interspecific hybrids between three eggplant (Solanum melongena L.) cultivars and Solanum macrocarpon L. Sci. Hortic. 2004, 101, 11–21. [Google Scholar] [CrossRef]
  45. Stedje, B.; Bukenya-Ziraba, R. RAPD variation in Solanum anguivi Lam. and S. aethiopicum L. (Solanaceae) in Uganda. Euphytica 2003, 131, 293–297. [Google Scholar] [CrossRef]
  46. Hassan, S.M.; Lester, R.N. Origin and domestication of the Brinjal eggplant, Solanum melongena L., from Solanum incanum L., in Africa and Asia. In Proceedings of the III Congreso Mundial De Plantas Solanaceas International Solanaceae Congress, Bogota, Columbia, 25–30 July 1988; pp. 25–30. [Google Scholar]
  47. Anaso, H.U. Comparative cytological study of Solanum aethiopicum Gilo group, Solanum aethiopicum Shum group and Solanum anguivi. Euphytica 1991, 53, 81–85. [Google Scholar] [CrossRef]
  48. Rakha, M.; Prohens, J.; Taher, D.; Wu, T.-h.; Solberg, S.Ø. Advances in Plant Breeding Strategies: Vegetable Crops: Volume 9: Fruits and Young Shoots; Jameel, M., Al-Khayri, S., Mohan, J., Johnson, D.V., Eds.; Springer International Publishing: Berlin/Heidelberg, Germany, 2021; pp. 163–203. [Google Scholar]
  49. Bletsos, F.; Roupakias, D.; Tsaktsira, M.; Scaltsoyjannes, A.; Thanassoulopoulos, C.C. Interspecific hybrids between three eggplant (Solanum melongena L.) cultivars and two wild species (Solanum torvum Sw. and Solanum sisymbriifolium Lam.). Plant Breed. 1998, 117, 159–164. [Google Scholar] [CrossRef]
  50. Kumchai, J.; Wei, Y.-C.; Lee, C.-Y.; Chen, F.-C.; Chin, S.-W. Production of interspecific hybrids between commercial cultivars of the eggplant (Solanum melongena L.) and its wild relative S. torvum. Genet. Mol. Res.: GMR 2013, 12, 755–764. [Google Scholar] [CrossRef]
Figure 1. Structure of 13 cp genomes. (A) The circle characteristic of the reference genome cp genomes (S. melongena). (B) The physical map of the circular cp genome of S. sisymbriifolium. Genes shown inside the outer circle are transcribed counterclockwise, and those outside are transcribed clockwise. Genes belonging to different functional groups are color-coded. A gray area in the inner circle indicates the GC content. (C) Comparison of the borders of large single-copy (LSC), small single-copy (SSC), and inverted repeat (IR) regions among the cp genomes of 13 cp genomes.
Figure 1. Structure of 13 cp genomes. (A) The circle characteristic of the reference genome cp genomes (S. melongena). (B) The physical map of the circular cp genome of S. sisymbriifolium. Genes shown inside the outer circle are transcribed counterclockwise, and those outside are transcribed clockwise. Genes belonging to different functional groups are color-coded. A gray area in the inner circle indicates the GC content. (C) Comparison of the borders of large single-copy (LSC), small single-copy (SSC), and inverted repeat (IR) regions among the cp genomes of 13 cp genomes.
Cimb 45 00185 g001
Figure 2. Codon, amino acid and sequence variation analyses of the cp genome. (A) Codon content and RSCU value of 20 amino acids and the stop codon of NN1 (Solanum aethiopicum L.) cp genome. (B) Numbers of encoded amino acids of NN1 cp genome. (C,D) Percentage of SNPs and InDels in each region. (E) The distribution of InDel length.
Figure 2. Codon, amino acid and sequence variation analyses of the cp genome. (A) Codon content and RSCU value of 20 amino acids and the stop codon of NN1 (Solanum aethiopicum L.) cp genome. (B) Numbers of encoded amino acids of NN1 cp genome. (C,D) Percentage of SNPs and InDels in each region. (E) The distribution of InDel length.
Cimb 45 00185 g002
Figure 3. Repeat sequences and SSRs in 13 cp genomes. (A) F, P, R, and C indicate the repeat types F (forward), P (palindrome), R (reverse), and C (complement), respectively; (B) Repeats with different lengths. (C) Repeat unit and amounts of SSR in thirteen cp genomes. (D) Percentage of different SSR types in the NN1 cp genome. (E) SSR repeat sequences in the NN1 cp genome. (F) Percentage of repeated sequences in the NN1 cp genome.
Figure 3. Repeat sequences and SSRs in 13 cp genomes. (A) F, P, R, and C indicate the repeat types F (forward), P (palindrome), R (reverse), and C (complement), respectively; (B) Repeats with different lengths. (C) Repeat unit and amounts of SSR in thirteen cp genomes. (D) Percentage of different SSR types in the NN1 cp genome. (E) SSR repeat sequences in the NN1 cp genome. (F) Percentage of repeated sequences in the NN1 cp genome.
Cimb 45 00185 g003
Figure 4. Variable regions and phylogenetic tree of the Solanum species. (A) Comparison of variable regions among 46 Solanum cp genomes. (B) Phylogenetic analysis of 46 Solanum cp genomes.
Figure 4. Variable regions and phylogenetic tree of the Solanum species. (A) Comparison of variable regions among 46 Solanum cp genomes. (B) Phylogenetic analysis of 46 Solanum cp genomes.
Cimb 45 00185 g004
Table 1. List of test materials.
Table 1. List of test materials.
Sample NumberCodeSpeciesGene Bank Accession NumberOrigin
NN1166Solanum aethiopicum L.MN218076South Africa
NN253Solanum aethiopicum L.MN218077Brazil
NN3137Solanum aethiopicum L.MN218078Ethiopia
NN4Y11Solanum aethiopicum L.MN218079Ethiopia
NN5177Solanum melongena L.MN218080South China
NN6LwpqSolanum macrocapon L.MN218081Laos
NN7ShfSolanum melongena L.MN218082Viet Nam
NN8131Solanum melongena L.MN218083South China
NN9DhqSolanum wrightii L.MN218084South China
NN12SqSolanum torvumMN218087South China
NN13CtqSolanum anguivi L.MN218088Thailand
NN14ctq-BSolanum anguivi L.MN218089South China
NN15SjqSolanum sisymbriifoliumMN218090South China
Table 2. Summary of sequencing data for 13 cp genomes.
Table 2. Summary of sequencing data for 13 cp genomes.
Sample IDInsert Size (bp)Raw Data (Mb)Clean Data (Mb)Reads Length (bp)Clean Data GC(%)Clean Data Q30(%)Average Organelle Depth
NN143059815210(150:150)36.6593.811452
NN243091118248(150:150)36.3694.41521
NN343067546034(150:150)36.4794.211336
NN443063745638(150:150)36.9193.951824
NN543065155792(150:150)36.4894.131599
NN643086708302(150:150)36.9295.391477
NN743085848307(150:150)37.6495.741581
NN843098939557(150:150)37.1395.641798
NN943085048193(150:150)37.1295.86700
NN1243010,1659777(150:150)38.195.543639
NN1343084088139(150:150)36.4295.782337
NN1443057815193(150:150)36.2495.31769
NN1543052814503(150:150)39.2793.5636
Table 3. Information on 13 cp gene annotation.
Table 3. Information on 13 cp gene annotation.
Sample IDTotal Length (bp)LSC Length (bp)SSC Length (bp)IRA/IRB Length (bp)Total GC Content (%)Protein-
Coding Gene Number
tRNA NumberrRNA NumberTotal
Gene
Number
NN1155,60686,16418,56225,44037.780428130
NN2155,59886,16018,55825,44037.780428130
NN3155,60786,16518,56225,44037.780428130
NN4155,61486,16818,56425,44137.780428130
NN5155,58186,19418,50125,44337.7180428130
NN6156,00486,54218,53825,46237.6180428130
NN7155,58586,19718,50225,44337.7180428130
NN8155,58586,19718,50225,44337.7180428130
NN9155,50686,23118,42925,42337.6880428130
NN12154,94285,64618,43625,43037.8180428130
NN13155,75386,26118,60825,44237.6680428130
NN14155,75286,26218,60825,44137.6780428130
NN15155,77186,40418,52525,42137.7679428129
Table 4. List of genes encoded by 13 cp genomes.
Table 4. List of genes encoded by 13 cp genomes.
FunctionGene GroupGene Name
Photosynthesis pathwaysPhotosystem Ipsa (A, B, C *, I, J)
Photosystem I assemblyycf (3, 4)
Photosystem IIpsb (A-F, H-L, N, T, Z)
F-type ATP synthaseatp (A, B, E, F, H, I)
NDH complexndh (A *, B #, C, D *, E *, F *, G *, H *, I *, J, K)
Component of cytochrome b6/f complexpet (A, B, D, L)
Inner envelope membranecemA
Cytochrome biogenesis proteinccsA *
Large subunit of RubiscorbcL
Structural RNAsTransfer RNAstrnH-GUG; trnK-UUU; trnQ-UUG; trnS-GCU; trnG-UCC; trnR-UCU; trnF-GAA; trnD-GUC; trnY-GUA; trnE-UUC; trnT-GGU; trnM-CAU; trnS-UGA; trnC-GCA; trnG-GCC; trnfM-CAU; trnS-GGA; trnT-UGU; trnL-UAA; trnI-CAU #; trnV-UAC; trnW-CCA; trnP-UGG; trnP-GGG; trnI-GAU #; trnA-UGC #; trnN-GUU#; trnL-UAG; trnR-ACG #; trnV-GAC #; trnL-CAA #
Ribosomal RNAsrrn (4.5 #, 5 #, 16 #, 23 #)
Genetic apparatusLarge subunit of ribosomal proteinrpl (2 #, 14, 16, 20, 22, 23 #, 32 *, 33, 36)
Small subunit of ribosomal proteinrps (2, 3, 4, 7 #, 8, 11, 12 #&, 14, 15 *, 16 0, 18, 19)
Subunits of RNA polymeraserpo (A, B, C1, C2)
Post-transcriptional modificationMaturaseMatK
Protein-
modifying
ATP-
dependent Clp protease proteolytic subunit
clpP
Unknown Proteins ycf (1 #, 2 #, 15)
#, Gene repeated in IR region; * gene in SSC region; &, gene has 2 separate transcription units; 0 gene is not in material NN15.
Table 5. InDel and SNP types of 13 cp genomes.
Table 5. InDel and SNP types of 13 cp genomes.
Sample IDSpeciesCDS
InDel
Intergenic InDelInsertionDeletionTotal InDelCDS
SNPs
Intergenic_SNPsts/tvTotal_SNPs
NN1S. aethiopicum L.783573390982421.11340
NN2S. aethiopicum L.6865834921002391.08339
NN3S. aethiopicum L.783573390992431.12342
NN4S. aethiopicum L.7815632881002441.12344
NN5S. melongena L.846468549300.4439
NN7S. melongena L.847469559300.4439
NN8S. melongena L.8454585310310.4441
NN6S. melongena L.5906431951113440.95455
NN9S. wrightii L.410359481072676060.95873
NN12S. torvum511472471192435990.91842
NN13S. indicum L.6875736931032681.01371
NN14S. indicum L.6875736931032671.01370
NN15S. sisymbriifolium512176501262596370.94896
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Yang, Q.; Li, Y.; Cai, L.; Gan, G.; Wang, P.; Li, W.; Li, W.; Jiang, Y.; Li, D.; Wang, M.; et al. Characteristics, Comparative Analysis, and Phylogenetic Relationships of Chloroplast Genomes of Cultivars and Wild Relatives of Eggplant (Solanum melongena). Curr. Issues Mol. Biol. 2023, 45, 2832-2846. https://doi.org/10.3390/cimb45040185

AMA Style

Yang Q, Li Y, Cai L, Gan G, Wang P, Li W, Li W, Jiang Y, Li D, Wang M, et al. Characteristics, Comparative Analysis, and Phylogenetic Relationships of Chloroplast Genomes of Cultivars and Wild Relatives of Eggplant (Solanum melongena). Current Issues in Molecular Biology. 2023; 45(4):2832-2846. https://doi.org/10.3390/cimb45040185

Chicago/Turabian Style

Yang, Qihong, Ye Li, Liangyu Cai, Guiyun Gan, Peng Wang, Weiliu Li, Wenjia Li, Yaqin Jiang, Dandan Li, Mila Wang, and et al. 2023. "Characteristics, Comparative Analysis, and Phylogenetic Relationships of Chloroplast Genomes of Cultivars and Wild Relatives of Eggplant (Solanum melongena)" Current Issues in Molecular Biology 45, no. 4: 2832-2846. https://doi.org/10.3390/cimb45040185

APA Style

Yang, Q., Li, Y., Cai, L., Gan, G., Wang, P., Li, W., Li, W., Jiang, Y., Li, D., Wang, M., Xiong, C., Chen, R., & Wang, Y. (2023). Characteristics, Comparative Analysis, and Phylogenetic Relationships of Chloroplast Genomes of Cultivars and Wild Relatives of Eggplant (Solanum melongena). Current Issues in Molecular Biology, 45(4), 2832-2846. https://doi.org/10.3390/cimb45040185

Article Metrics

Back to TopTop