Pan-plastome approach empowers the assessment of genetic variation in cultivated Capsicum species

The pan-plastome of Hemerocallis citrina reveals new insights into the genetic diversity and cultivation history of an economically important food plant

BACKGROUND

Hemerocallis citrina Baroni (Huang hua cai in Chinese) is a perennial herbaceous plant grown for its flower buds that are eaten fresh or dried and is known as the vegetarian three treasures. The nuclear genome of H. citrina has been reported, but the intraspecific variation of the plastome (plastid genome) has not yet been studied. Therefore, the panplastome of this species collected from diverse locations is reported here for the first time.

RESULTS

In this study, 65 H. citrina samples were resequenced, de novo assembled, and aligned with the published plastome of H. citrina to resolve the H. citrina panplastome. The sizes of the 65 newly assembled complete plastomes of H. citrina ranged from 156,048 bp to 156,263 bp, and the total GC content ranged from 37.31 to 37.34%. The structure of the complete plastomes showed a typical tetrameric structure, including a large single copy (LSC), a small single copy (SSC), and a pair of inverted repeat regions (IRA and IRB). Many nucleotide variants were identified between plastomes, among which the variants in the intergenic spacer region were the most abundant, with the highest number of variants concentrated in the LSC region. Based on the phylogenetic tree constructed using the ML method, population structure analysis, and principal component analysis (PCA), the panplastome data were subdivided into five genetic clusters. The C5 genetic cluster was mostly represented by samples from Qidong, Hunan Province, while samples from Shanxi and Shaanxi Provinces were classified into the C4 genetic cluster. The greatest genetic diversity was found in the C1 genetic cluster, and the greatest genetic distance between any two clusters was found between the C4 and C5 clusters.

CONCLUSION

The resolution of the panplastome and the analysis of the population structure of H. citrina plastomes provide important data for future breeding projects and germplasm preservation.

Genomic divergence and demographic history of Quercus aliena populations

BACKGROUND

Quercus aliena is a major montane tree species of subtropical and temperate forests in China, with important ecological and economic value. In order to reveal the species' population dynamics, genetic diversity, genetic structure, and association with mountain habitats during the evolutionary process, we re-sequenced the genomes of 72 Q. aliena individuals.

RESULTS

The whole chloroplast and nuclear genomes were used for this study. Phylogenetic analysis using the chloroplast genome dataset supported four clades of Q. aliena, while the nuclear dataset supported three major clades. Sex-biased dispersal had a critical role in causing discordance between the chloroplast and nuclear genomes. Population structure analysis showed two groups in Q. aliena. The effective population size sharply declined 1 Mya, coinciding with the Poyang Glaciation in Eastern China. Using genotype-climate association analyses, we found a positive correlation between allele frequency variation in SNPs and temperature, suggesting the species has the capacity to adapt to changing temperatures.

CONCLUSION

Overall, this study illustrates the genetic divergence, genomic variation, and evolutionary processes behind the demographic history of Q. aliena.

The complete chloroplast genome of the halophyte flowering plant Suaeda monoica from Jeddah, Saudi Arabia

The complete chloroplast genome (plastome) of the annual flowering halophyte herb Suaeda monoica Forssk. ex J. F. Gmel. family (Amaranthaceae) that grows in Jeddah, Saudi Arabia, was identified for the first time in this study. Suaeda monoica is a medicinal plant species whose taxonomic classification remains controversial. Further, studying the species is useful for current conservation and management efforts. In the current study, the full chloroplast genome S. monoica was reassembled using whole-genome next-generation sequencing and compared with the previously published chloroplast genomes of Suaeda species. The chloroplast genome size of Suaeda monoica was 151,789 bp, with a single large copy of 83,404 bp, a small single copy of 18,007 bp and two inverted repeats regions of 25,189 bp. GC content in the whole genome was 36.4%. The cp genome included 87 genes that coded for proteins, 37 genes coding for tRNA, 8 genes coding for rRNA and one non-coding pseudogene. Five chloroplast genome features were compared between S. monoica and S. japonica, S. glauca, S. salsa, S. malacosperma and S. physophora. Among Suaeda genus and equal to most angiosperms chloroplast genomes, the RSCU values were conservative. Two pseudogenes (accD and ycf1), rpl16 intron and ndhF-rpl32 intergenic spacer, were highlighted as suitable DNA barcodes for different Suaeda species. Phylogenetic analyses show Suaeda cluster into three main groups; one in which S. monoica was closer to S. salsa. The obtained result provided valuable information on the characteristics of the S. monoica chloroplast genome and the phylogenetic relationships.

Evolutionary insights and expression dynamics of the CaNFYB transcription factor gene family in pepper (Capsicum annuum) under salinity stress

Introduction: The Capsicum annuum nuclear factor Y subunit B (CaNFYB) gene family plays a significant role in diverse biological processes, including plant responses to abiotic stressors such as salinity. Methods: In this study, we provide a comprehensive analysis of the CaNFYB gene family in pepper, encompassing their identification, structural details, evolutionary relationships, regulatory elements in promoter regions, and expression profiles under salinity stress. Results and discussion: A total of 19 CaNFYB genes were identified and subsequently characterized based on their secondary protein structures, revealing conserved domains essential for their functionality. Chromosomal distribution showed a non-random localization of these genes, suggesting potential clusters or hotspots for NFYB genes on specific chromosomes. The evolutionary analysis focused on pepper and comparison with other plant species indicated a complex tapestry of relationships with distinct evolutionary events, including gene duplication. Moreover, promoter cis-element analysis highlighted potential regulatory intricacies, with notable occurrences of light-responsive and stress-responsive binding sites. In response to salinity stress, several CaNFYB genes demonstrated significant temporal expression variations, particularly in the roots, elucidating their role in stress adaptation. Particularly CaNFYB01, CaNFYB18, and CaNFYB19, play a pivotal role in early salinity stress response, potentially through specific regulatory mechanisms elucidated by their cis-elements. Their evolutionary clustering with other Solanaceae family members suggests conserved ancestral functions vital for the family's survival under stress. This study provides foundational knowledge on the CaNFYB gene family in C. annuum, paving the way for further research to understand their functional implications in pepper plants and relative species and their potential utilization in breeding programs to enhance salinity tolerance.

Genome-wide identification of CaARR-Bs transcription factor gene family in pepper and their expression patterns under salinity stress

In plants, ARRs-B transcription factors play a crucial role in regulating cytokinin signal transduction, abiotic stress resistance, and plant development. A number of adverse environmental conditions have caused severe losses for the pepper (Capsicum annuum L.)-a significant and economically important vegetable. Among the transcription factors of the type B-ARRs family, multiple members have different functions. In pepper, only a few members of the ARRs-B family have been reported and characterized. The current study aimed to characterize ARRs-B transcription factors in C. annuum, including phylogenetic relationships, gene structures, protein motif arrangement, and RT-qPCR expression analyses and their role in salinity stress. In total, ten genes encode CaARRs-B transcription factors (CaARR1 to CaARR10) from the largest subfamily of type-B ARRs were identified in C. annum. The genome-wide analyses of the CaARRs-B family in C. annuum were performed based on the reported ARRs-B genes in Arabidopsis. An analysis of homologous alignments of candidate genes, including their phylogenetic relationships, gene structures, conserved domains, and qPCR expression profiles, was conducted. In comparison with other plant ARRs-B proteins, CaARRs-B proteins showed gene conservation and potentially specialized functions. In addition, tissue-specific expression profiles showed that CaARRs-B genes were differentially expressed, suggesting functionally divergent. CaARRs-B proteins had a typical conserved domain, including AAR-like (pfam: PF00072) and Myb DNA binding (pfam: PF00249) domains. Ten of the CaARRs-B genes were asymmetrically mapped on seven chromosomes in Pepper. Additionally, the phylogenetic tree of CaARRs-B genes from C. annuum and other plant species revealed that CaARRs-B genes were classified into four clusters, which may have evolved conservatively. Further, using quantitative real-time qRT-PCR, the study assessed the expression patterns of CaARRs-B genes in Capsicum annuum seedlings subjected to salt stress. The study used quantitative real-time qRT-PCR to examine CaARRs-B gene expression in Capsicum annuum seedlings under salt stress. Roots exhibited elevated expression of CaARR2 and CaARR9, while leaves showed decreased expression for CaARR3, CaARR4, CaARR7, and CaARR8. Notably, no amplification was observed for CaARR10. This research sheds light on the roles of CaARRs-B genes in pepper's response to salinity stress. These findings enrich our comprehension of the functional implications of CaARRs-B genes in pepper, especially in responding to salinity stress, laying a solid groundwork for subsequent in-depth studies and applications in the growth and development of Capsicum annuum.

Comparative Plastomes of Curcuma alismatifolia (Zingiberaceae) Reveal Diversified Patterns among 56 Different Cut-Flower Cultivars

Curcuma alismatifolia (Zingiberaceae) is an ornamental species with high economic value due to its recent rise in popularity among floriculturists. Cultivars within this species have mixed genetic backgrounds from multiple hybridization events and can be difficult to distinguish via morphological and histological methods alone. Given the need to improve identification resources, we carried out the first systematic study using plastomic data wherein genomic evolution and phylogenetic relationships from 56 accessions of C. alismatifolia were analyzed. The newly assembled plastomes were highly conserved and ranged from 162,139 bp to 164,111 bp, including 79 genes that code for proteins, 30 tRNA genes, and 4 rRNA genes. The A/T motif was the most common of SSRs in the assembled genomes. The Ka/Ks values of most genes were less than 1, and only two genes had Ka/Ks values above 1, which were rps15 (1.15), and ndhl (1.13) with petA equal to 1. The sequence divergence between different varieties of C. alismatifolia was large, and the percentage of variation in coding regions was lower than that in the non-coding regions. Such data will improve cultivar identification, marker assisted breeding, and preservation of germplasm resources.

The First Complete Chloroplast Genome of Cordia monoica: Structure and Comparative Analysis

Cordia monoica is a member of the Boraginaceae family. This plant is widely distributed in tropical regions and has a great deal of medical value as well as economic importance. In the current study, the complete chloroplast (cp) genome of C. monoica was sequenced, assembled, annotated, and reported. This circular chloroplast genome had a size of 148,711 bp, with a quadripartite structure alternating between a pair of repeated inverted regions (26,897-26,901 bp) and a single copy region (77,893 bp). Among the 134 genes encoded by the cp genome, there were 89 protein-coding genes, 37 transfer RNA (tRNA) genes, and 8 ribosomal RNA (rRNA) genes. A total of 1387 tandem repeats were detected, with the hexanucleotides class making up 28 percent of the repeats. Cordia monoica has 26,303 codons in its protein-coding regions, and leucine amino acid was the most frequently encoded amino acid in contrast to cysteine. In addition, 12 of the 89 protein-coding genes were found to be under positive selection. The phyloplastomic taxonomical clustering of the Boraginaceae species provides further evidence that chloroplast genome data are reliable not only at family level but also in deciphering the phylogeny at genus level (e.g., Cordia).

The pan-plastome of tartary buckwheat (fagopyrum tataricum): key insights into genetic diversity and the history of lineage divergence

BACKGROUND

Tartary buckwheat (Fagopyrum tataricum) is an important food and medicine crop plant, which has been cultivated for 4000 years. A nuclear genome has been generated for this species, while an intraspecific pan-plastome has yet to be produced. As such a detailed understanding of the maternal genealogy of Tartary buckwheat has not been thoroughly investigated.

RESULTS

In this study, we de novo assembled 513 complete plastomes of Fagopyrum and compared with 8 complete plastomes of Fagopyrum downloaded from the NCBI database to construct a pan-plastome for F. tartaricum and resolve genomic variation. The complete plastomes of the 513 newly assembled Fagopyrum plastome sizes ranged from 159,253 bp to 159,576 bp with total GC contents ranged from 37.76 to 37.97%. These plastomes all maintained the typical quadripartite structure, consisting of a pair of inverted repeat regions (IRA and IRB) separated by a large single copy region (LSC) and a small single copy region (SSC). Although the structure and gene content of the Fagopyrum plastomes are conserved, numerous nucleotide variations were detected from which population structure could be resolved. The nucleotide variants were most abundant in the non-coding regions of the genome and of those the intergenic regions had the most. Mutational hotspots were primarily found in the LSC regions. The complete 521 Fagopyrum plastomes were divided into five genetic clusters, among which 509 Tartary buckwheat plastomes were divided into three genetic clusters (Ft-I/Ft-II/Ft-III). The genetic diversity in the Tartary buckwheat genetic clusters was the greatest in Ft-III, and the genetic distance between Ft-I and Ft-II was the largest. Based on the results of population structure and genetic diversity analysis, Ft-III was further subdivided into three subgroups Ft-IIIa, Ft-IIIb, and Ft-IIIc. Divergence time estimation indicated that the genera Fagopyrum and Rheum (rhubarb) shared a common ancestor about 48 million years ago (mya) and that intraspecies divergence in Tartary buckwheat began around 0.42 mya.

CONCLUSIONS

The resolution of pan-plastome diversity in Tartary buckwheat provides an important resource for future projects such as marker-assisted breeding and germplasm preservation.

Poaceae Chloroplast Genome Sequencing: Great Leap Forward in Recent Ten Years

The first complete chloroplast genome of rice (Oryza sativa) was published in 1989, ushering in a new era of studies of chloroplast genomics in Poaceae. Progresses in Next-Generation Sequencing (NGS) and Third-Generation Sequencing (TGS) technologiesand in the development of genome assembly software, have significantly advanced chloroplast genomics research. Poaceae is one of the most targeted families in chloroplast genome research because of its agricultural, ecological, and economic importance. Over the last 30 years, 2,050 complete chloroplast genome sequences from 40 tribes and 282 genera have been generated, most (97%) of them in the recent ten years. The wealth of data provides the groundwork for studies on species evolution, phylogeny, genetic transformation, and other aspects of Poaceae chloroplast genomes. As a result, we have gained a deeper understanding of the properties of Poaceae chloroplast genomes. Here, we summarize the achievements of the studies of the Poaceae chloroplast genomes and envision the challenges for moving the area ahead.

Characterization and phylogenetic analysis of the complete mitochondrial genome of Mango tilapia (Sarotherodon galilaeus: Cichlidae)

BACKGROUND

Sarotherodon galilaeus (Linné, 1758) is a member of the family Cichlidae, which is considered the most important aquaculture freshwater species endemic to Africa and the Middle East. The genetics and molecular biology of this species are rare. This requires more comprehensive mitochondrial genomes-based phylogenetics to enhance understanding of the relationship and delineate this species.

METHODS AND RESULTS

Here, we assembled the complete mitogenome of S. galilaeus using Illumina high-throughput sequencing technology. The mango tilapia mitogenome was 16,631 bp in length with an AT composition of 53.4% and 46.4% GC content. It encodes 37 genes comprising two ribosomal RNA genes (rRNAs), 22 transfer RNA genes (tRNAs), and 13 protein-coding genes (PCGs) as well as the D-loop known as the control region. The phylogenetic tree was conducted to provide a relationship within the haplotilapiine lineage based on the maximum likelihood method, and the newly sequenced S. galilaeus was clustered with other Sarotherodon species.

CONCLUSION

Our results provide a new perception of the genetic basis of S. galilaeus species for further research on systematics, evolution, population genetics, and molecular ecology.

Pan-Plastome of Greater Yam (Dioscorea alata) in China: Intraspecific Genetic Variation, Comparative Genomics, and Phylogenetic Analyses

Dioscorea alata L. (Dioscoreaceae), commonly known as greater yam, water yam, or winged yam, is a popular tuber vegetable/food crop worldwide, with nutritional, health, and economical importance. China is an important domestication center of D. alata, and hundreds of cultivars (accessions) have been established. However, genetic variations among Chinese accessions remain ambiguous, and genomic resources currently available for the molecular breeding of this species in China are very scarce. In this study, we generated the first pan-plastome of D. alata, based on 44 Chinese accessions and 8 African accessions, and investigated the genetic variations, plastome evolution, and phylogenetic relationships within D. alata and among members of the section Enantiophyllum. The D. alata pan-plastome encoded 113 unique genes and ranged in size from 153,114 to 153,161 bp. A total of four whole-plastome haplotypes (Haps I-IV) were identified in the Chinese accessions, showing no geographical differentiation, while all eight African accessions shared the same whole-plastome haplotype (Hap I). Comparative genomic analyses revealed that all four whole plastome haplotypes harbored identical GC content, gene content, gene order, and IR/SC boundary structures, which were also highly congruent with other species of Enantiophyllum. In addition, four highly divergent regions, i.e., trnC-petN, trnL-rpl32, ndhD-ccsA, and exon 3 of clpP, were identified as potential DNA barcodes. Phylogenetic analyses clearly separated all the D. alata accessions into four distinct clades corresponding to the four haplotypes, and strongly supported that D. alata was more closely related to D. brevipetiolata and D. glabra than D. cirrhosa, D. japonica, and D. polystachya. Overall, these results not only revealed the genetic variations among Chinese D. alata accessions, but also provided the necessary groundwork for molecular-assisted breeding and industrial utilization of this species.

A comparative plastomic analysis of Ziziphus jujuba var. spinosa (Bunge) Hu ex H. F. Chow and implication of the origin of Chinese jujube

Comparative plastomics can be used to explicitly dissect various types of plastome variation. In the present study, the plastome variation pattern of Ziziphus jujuba var. spinosa (also called sour jujube) and its phylogenomic relationship with Chinese jujube were investigated. Plastomes of 21 sour jujube individuals were sequenced and assembled. The length of the sour jujube plastomes ranged between 159399 and 161279 bp. The plastomes exhibited collinearity of structure, gene order and content. The most divergent regions were located in the intergenic spacers, such as trnR-UCU-atpA and psbZ-trnG-UCC. Sliding window analysis demonstrated that the sequence variation among the sour jujube plastomes was relatively low. Sixty-two to 76 SSRs with 4 motif types were identified in the sour jujube plastomes with a predominant motif type of A/T. Three protein-coding genes exhibited higher nonsynonymous/synonymous substitution ratios, indicating that these genes may undergo positive selection. A total of 80 SNPs were detected and 1266 potential RNA editing sites of 23 protein-coding genes were predicted. In the phylogenomic tree constructed, sour jujube has a sister relationship to Chinese jujube, which indicates that Chinese jujube may have originated or been domesticated from sour jujube. The present study explicitly investigated the individual-level plastome variation of sour jujube and provides potential valuable molecular markers for future genetic-related study of this lineage.

Mitogenomic Features and Evolution of the Nile River Dominant Tilapiine Species (Perciformes: Cichlidae)

To better understand the diversity and evolution of cichlids, we sequenced, assembled, and annotated the complete mitochondrial genomes of three Nile tilapiine species (Coptodon zillii, Oreochromis niloticus, and Sarotherodon galilaeus) dominating the Nile River waters. Our results showed that the general mitogenomic features were conserved among the Nile tilapiine species. The genome length ranged from 16,436 to 16,631 bp and a total of 37 genes were identified (two ribosomal RNA genes (rRNAs), 22 transfer RNA genes (tRNAs), 13 protein-coding genes (PCGs), and 1 control region). The ND6 was the only CDS that presented a negative AT skew and a positive GC skew. The most extended repeat sequences were in the D-loop followed by the pseudogenes (trnS^GCU). The ND5 showed relatively high substitution rates whereas ATP8 had the lowest substitution rate. The codon usage bias displayed a greater quantity of NNA and NNC at the third position and anti-bias against NNG. The phylogenetic relationship based on the complete mitogenomes and CDS was able to differentiate the three species as previously reported. This study provides new insight into the evolutionary connections between various subfamilies within cichlids while providing new molecular data that can be applied to discriminate between Nile tilapiine species and their populations.

Chloroplast genome assemblies and comparative analyses of commercially important Vaccinium berry crops

Vaccinium is a large genus of shrubs that includes a handful of economically important berry crops. Given the numerous hybridizations and polyploidization events, the taxonomy of this genus has remained the subject of long debate. In addition, berries and berry-based products are liable to adulteration, either fraudulent or unintentional due to misidentification of species. The availability of more genomic information could help achieve higher phylogenetic resolution for the genus, provide molecular markers for berry crops identification, and a framework for efficient genetic engineering of chloroplasts. Therefore, in this study we assembled five Vaccinium chloroplast sequences representing the economically relevant berry types: northern highbush blueberry (V. corymbosum), southern highbush blueberry (V. corymbosum hybrids), rabbiteye blueberry (V. virgatum), lowbush blueberry (V. angustifolium), and bilberry (V. myrtillus). Comparative analyses showed that the Vaccinium chloroplast genomes exhibited an overall highly conserved synteny and sequence identity among them. Polymorphic regions included the expansion/contraction of inverted repeats, gene copy number variation, simple sequence repeats, indels, and single nucleotide polymorphisms. Based on their in silico discrimination power, we suggested variants that could be developed into molecular markers for berry crops identification. Phylogenetic analysis revealed multiple origins of highbush blueberry plastomes, likely due to the hybridization events that occurred during northern and southern highbush blueberry domestication.

Variations in genetic diversity in cultivated Pistacia chinensis

Identification of the evolution history and genetic diversity of a species is important in the utilization of novel genetic variation in this species, as well as for its conservation. Pistacia chinensis is an important biodiesel tree crop in China, due to the high oil content of its fruit. The aim of this study was to uncover the genetic structure of P. chinensis and to investigate the influence of intraspecific gene flow on the process of domestication and the diversification of varieties. We investigated the genetic structure of P. chinensis, as well as evolution and introgression in the subpopulations, through analysis of the plastid and nuclear genomes of 39 P. chinensis individuals from across China. High levels of variation were detected in the P. chinensis plastome, and 460 intraspecific polymorphic sites, 104 indels and three small inversions were identified. Phylogenetic analysis and population structure using the plastome dataset supported five clades of P. chinensis. Population structure analysis based on the nuclear SNPs showed two groups, clearly clustered together, and more than a third of the total individuals were classified as hybrids. Discordance between the plastid and nuclear genomes suggested that hybridization events may have occurred between highly divergent samples in the P. chinensis subclades. Most of the species in the P. chinensis subclade diverged between the late Miocene and the mid-Pliocene. The processes of domestication and cultivation have decreased the genetic diversity of P. chinensis. The extensive variability and structuring of the P. chinensis plastid together with the nuclear genomic variation detected in this study suggests that much unexploited genetic diversity is available for improvement in this recently domesticated species.

Diversification of chiles (Capsicum, Solanaceae) through time and space: New insights from genome-wide RAD-seq data

Capsicum L. encompasses 43 American species, including the five domesticated worldwide consumed sweet and hot chiles. This study presents new, updated and age-calibrated phylogenetic hypothesis on the genus evolution incorporating nearly all currently accepted Capsicum species. A new model of spatial diversification of Capsicum is proposed based on analyses of several thousands of genome-wide RAD-seq derived SNPs. Maximum likelihood approaches were used to reconstruct phylogenies and to estimate dates of all major splits. Ancestral ranges were inferred and diversification events were modeled in a time frame using a Bayesian approach. Nine clades corresponding to genetically and (mostly) geographically well-defined lineages, which diversified starting around mid-upper Miocene, were recovered with strong support. The Northern and Central Andes were inferred to represent the most likely ancestral range of the genus Capsicum. A few early vicariant and dispersal events were estimated to have driven the geographic divergence of the main Capsicum clades. Each lineage was inferred to have diversified within a distinct region of South America and expanded geographically to different extent. Extant species diversification was inferred to have begun at the beginning of Pliocene and continued through the Pleistocene. The Central Andes, represented mainly by the territory of present-day Bolivia, were proposed to play a central role in the diversification of lineages comprising domesticated capsicums and their wild allies. The genome-wide approach allowed for high resolution and support of deep phylogenetic nodes providing novel insights into the affinities of major lineages and clades as well as on the geographic expansion of Capsicum. This study provides the first dated evolutionary history of the genus encompassing most of the chile species diversity.

Characterization of the plastome of Physaliscordata and comparative analysis of eight species of Physalis sensu stricto

In this study, we sequenced, assembled, and annotated the plastome of Physaliscordata Mill. and compared it with seven species of the genus Physalis sensu stricto. Sequencing, annotating, and comparing plastomes allow us to understand the evolutionary mechanisms associated with physiological functions, select possible molecular markers, and identify the types of selection that have acted in different regions of the genome. The plastome of P.cordata is 157,000 bp long and presents the typical quadripartite structure with a large single-copy (LSC) region of 87,267 bp and a small single-copy (SSC) region of 18,501 bp, which are separated by two inverted repeat (IRs) regions of 25,616 bp each. These values are similar to those found in the other species, except for P.angulata L. and P.pruinosa L., which presented an expansion of the LSC region and a contraction of the IR regions. The plastome in all Physalis species studied shows variation in the boundary of the regions with three distinct types, the percentage of the sequence identity between coding and non-coding regions, and the number of repetitive regions and microsatellites. Four genes and 10 intergenic regions show promise as molecular markers and eight genes were under positive selection. The maximum likelihood analysis showed that the plastome is a good source of information for phylogenetic inference in the genus, given the high support values and absence of polytomies. In the Physalis plastomes analyzed here, the differences found, the positive selection of genes, and the phylogenetic relationships do not show trends that correspond to the biological or ecological characteristics of the species studied.

The Potential Role of Plastome Copy Number as a Quality Biomarker for Plant Products using Real-time Quantitative Polymerase Chain Reaction

Background: Plastids are plant-specific semi-autonomous self-replicating organelles, containing circular DNA molecules called plastomes. Plastids perform crucial functions, including photosynthesis, stress perception and response, synthesis of metabolites, and storage. The plastome and plastid numbers have been shown to be modulated by developmental stage and environmental stimuli and have been used as a biomarker (identification of plant species) and biosensor (an indicator of abiotic and biotic stresses). However, the determination of plastome sequence and plastid number is a laborious process requiring sophisticated equipment. Methods: This study proposes using plastome copy number (PCN), which can be determined rapidly by real-time quantitative polymerase chain reaction (RT-qPCR) as a plant product quality biomarker. This study shows that the PCN log₁₀ and range PCN log₁₀ values calculated from RT-qPCR data, which was obtained for two years from leaves and lint samples of cotton and seed samples of cotton, rice, soybean, maize, and sesame can be used for assessing the quality of the samples. Results: Observation of lower range PCN log₁₀ values for CS (0.31) and CR (0.58) indicated that the PCN showed little variance from the mean PCN log₁₀ values for CS (3.81) and CR (3.85), suggesting that these samples might have encountered ambient environmental conditions during growth and/ or post-harvest storage and processing. This conclusion was further supported by observation of higher range PCN log₁₀ values for RS (3.09) versus RP (0.05), where rice seeds in the RP group had protective hull covering compared to broken hull-less seeds in the RS group. To further support that PCN is affected by external factors, rice seeds treated with high temperatures and pathogens exhibited lower PCN values when compared to untreated seeds. Furthermore, the range PCN log₁₀ values were found to be high for cotton leaf (CL) and lint (Clt) sample groups, 4.11 and 3.63, respectively, where leaf and lint samples were of different sizes, indicating that leaf samples might be of different developmental stage and lint samples might have been processed differently, supporting that the PCN is affected by both internal and external factors, respectively. Moreover, PCN log₁₀ values were found to be plant specific, with oil containing seeds such as SeS (6.49) and MS (5.05) exhibiting high PCN log₁₀ values compared to non-oil seeds such as SS (1.96). Conclusion: In conclusion, it was observed that PCN log₁₀ values calculated from RT-qPCR assays were specific to plant species and the range of PCN log₁₀ values can be directly correlated to the internal and external factors and, therefore might be used as a potential biomarker for assessing the quality of plant products.

Phylogenomics and Genetic Diversity of Arnebiae Radix and Its Allies (Arnebia, Boraginaceae) in China

Arnebiae Radix is a traditional medicine with pleiotropic properties that has been used for several 100 years. There are five species of Arnebia in China, and the two species Arnebia euchroma and Arnebia guttata are the source plants of Arnebiae Radix according to the Chinese Pharmacopoeia. Molecular markers that permit species identification and facilitate studies of the genetic diversity and divergence of the wild populations of these two source plants have not yet been developed. Here, we sequenced the chloroplast genomes of 56 samples of five Arnebia species using genome skimming methods. The Arnebia chloroplast genomes exhibited quadripartite structures with lengths from 149,539 and 152,040 bp. Three variable markers (rps16-trnQ, ndhF-rpl32, and ycf1b) were identified, and these markers exhibited more variable sites than universal chloroplast markers. The phylogenetic relationships among the five Arnebia species were completely resolved using the whole chloroplast genome sequences. Arnebia arose during the Oligocene and diversified in the middle Miocene; this coincided with two geological events during the late Oligocene and early Miocene: warming and the progressive uplift of Tianshan and the Himalayas. Our analyses revealed that A. euchroma and A. guttata have high levels of genetic diversity and comprise two and three subclades, respectively. The two clades of A. euchroma exhibited significant genetic differences and diverged at 10.18 Ma in the middle Miocene. Three clades of A. guttata diverged in the Pleistocene. The results provided new insight into evolutionary history of Arnebia species and promoted the conservation and exploitation of A. euchroma and A. guttata.

The complete chloroplast genome of the Egyptian henbane (Hyoscyamus muticus L., Solanaceae)

Egyptian henbane (Hyoscyamus muticus L. Mant. Pl. 1767) is an important medicinal plant of the family Solanaceae. Its complete chloroplast (cp) genome was assembled using Illumina high-throughput sequencing technology to contribute to its conservation genetics studies. Here, we report the complete sequence of the chloroplast genome of H. muticus. The cp genome was 156,271 bp in length with an asymmetric base composition (30.9% A, 18.9% C, 18.6% G and 31.6% T). It encodes 132 genes comprising 87 protein-coding genes, 29 tRNA genes, and 8 rRNA genes. The overall GC content of the H. muticus chloroplast genome was 37.5%, and the corresponding values in the large single-copy (LSC), the small single-copy (SSC), and the inverted repeat (IR) regions were 35.5%, 31.5%, and 43.0%, respectively. The complete chloroplast genome sequence was deposited to the GenBank (NCBI, Accession number: MZ450974). The maximum likelihood phylogenetic analysis showed that the H. muticus and H. niger were clustered into one clade with strong support values, indicating their closer relationship.

Entailing the Next-Generation Sequencing and Metabolome for Sustainable Agriculture by Improving Plant Tolerance

Crop production is a serious challenge to provide food for the 10 billion individuals forecasted to live across the globe in 2050. The scientists' emphasize establishing an equilibrium among diversity and quality of crops by enhancing yield to fulfill the increasing demand for food supply sustainably. The exploitation of genetic resources using genomics and metabolomics strategies can help generate resilient plants against stressors in the future. The innovation of the next-generation sequencing (NGS) strategies laid the foundation to unveil various plants' genetic potential and help us to understand the domestication process to unmask the genetic potential among wild-type plants to utilize for crop improvement. Nowadays, NGS is generating massive genomic resources using wild-type and domesticated plants grown under normal and harsh environments to explore the stress regulatory factors and determine the key metabolites. Improved food nutritional value is also the key to eradicating malnutrition problems around the globe, which could be attained by employing the knowledge gained through NGS and metabolomics to achieve suitability in crop yield. Advanced technologies can further enhance our understanding in defining the strategy to obtain a specific phenotype of a crop. Integration among bioinformatic tools and molecular techniques, such as marker-assisted, QTLs mapping, creation of reference genome, de novo genome assembly, pan- and/or super-pan-genomes, etc., will boost breeding programs. The current article provides sequential progress in NGS technologies, a broad application of NGS, enhancement of genetic manipulation resources, and understanding the crop response to stress by producing plant metabolites. The NGS and metabolomics utilization in generating stress-tolerant plants/crops without deteriorating a natural ecosystem is considered a sustainable way to improve agriculture production. This highlighted knowledge also provides useful research that explores the suitable resources for agriculture sustainability.

Phylogeny of Crataegus (Rosaceae) based on 257 nuclear loci and chloroplast genomes: evaluating the impact of hybridization

Background

Hawthorn species (Crataegus L.; Rosaceae tribe Maleae) form a well-defined clade comprising five subgeneric groups readily distinguished using either molecular or morphological data. While multiple subsidiary groups (taxonomic sections, series) are recognized within some subgenera, the number of and relationships among species in these groups are subject to disagreement. Gametophytic apomixis and polyploidy are prevalent in the genus, and disagreement concerns whether and how apomictic genotypes should be recognized taxonomically. Recent studies suggest that many polyploids arise from hybridization between members of different infrageneric groups.

Methods

We used target capture and high throughput sequencing to obtain nucleotide sequences for 257 nuclear loci and nearly complete chloroplast genomes from a sample of hawthorns representing all five currently recognized subgenera. Our sample is structured to include two examples of intersubgeneric hybrids and their putative diploid and tetraploid parents. We queried the alignment of nuclear loci directly for evidence of hybridization, and compared individual gene trees with each other, and with both the maximum likelihood plastome tree and the nuclear concatenated and multilocus coalescent-based trees. Tree comparisons provided a promising, if challenging (because of the number of comparisons involved) method for visualizing variation in tree topology. We found it useful to deploy comparisons based not only on tree-tree distances but also on a metric of tree-tree concordance that uses extrinsic information about the relatedness of the terminals in comparing tree topologies.

Results

We obtained well-supported phylogenies from plastome sequences and from a minimum of 244 low copy-number nuclear loci. These are consistent with a previous morphology-based subgeneric classification of the genus. Despite the high heterogeneity of individual gene trees, we corroborate earlier evidence for the importance of hybridization in the evolution of Crataegus. Hybridization between subgenus Americanae and subgenus Sanguineae was documented for the origin of Sanguineae tetraploids, but not for a tetraploid Americanae species. This is also the first application of target capture probes designed with apple genome sequence. We successfully assembled 95% of 257 loci in Crataegus, indicating their potential utility across the genera of the apple tribe.

Comparative plastid genomics of Mazaceae: focusing on a new recognized genus, Puchiumazus

Six Mazaceae plastomes were assembled in this study, including the newly recognized genus, Puchiumazus. Comparative plastid genomic analysis provided new insights into Mazaceae. The phylogenetic categorization of Mazus lanceifolius (Mazaceae) has long been uncertain. In 2021, the scholars Bo Li, D. G. Zhang, and C. L. Xiang republished M. lanceifolius as a new species Puchiumazus lanceifolius, within a new genus Puchiumazus. However, there is little plastome information on Mazaceae. Following the publishing of the new genus Puchiumazus, it is now necessary to study the Mazaceae plastome features to comprehensively understand this young family. The Mazaceae plastomes all have a typical quartile structure. The plastomes have a size ranging from 152,388 to 154,252 bp, and each plastome contains 112 unique genes, including 78 protein-coding genes, 4 rRNA genes, and 30 tRNA genes. A comparative analysis showed that these plastome sequences are highly conserved. Furthermore, we identified four relatively hypervariable regions (trnQ-UUC-psbK, trnS-GCU- trnS-CGA, trnT-UGU-trnL-UAA and ycf1) that can be used as potential DNA barcodes for the identification of this clade. Phylogenetic relationships based on the whole plastome sequences of 25 samples of 14 genera of Lamiales placed M. lanceifolius in the basal clade of the family Mazaceae, with 100% bootstrap support. In summary, the M. lanceifolius results indicate that a new monotype genus (Puchiumazus) should be established at the whole-plastome level. This study provides plastid genomic resources for exploring the phylogeny of Mazaceae.

Plastome variation and phylogeny of Taxillus (Loranthaceae)

Several molecular phylogenetic studies of the mistletoe family Loranthaceae have been published such that now the general pattern of relationships among the genera and their biogeographic histories are understood. Less is known about species relationships in the larger (> 10 species) genera. This study examines the taxonomically difficult genus Taxillus composed of 35–40 Asian species. The goal was to explore the genetic diversity present in Taxillus plastomes, locate genetically variable hotspots, and test these for their utility as potential DNA barcodes. Using genome skimming, complete plastomes, as well as nuclear and mitochondrial rDNA sequences, were newly generated for eight species. The plastome sequences were used in conjunction with seven publicly available Taxillus sequences and three sequences of Scurrula, a close generic relative. The Taxillus plastomes ranged from 121 to 123 kbp and encoded 90–93 plastid genes. In addition to all of the NADH dehydrogenase complex genes, four ribosomal genes, infA and four intron-containing tRNA genes were lost or pseudogenized in all of the Taxillus and Scurrula plastomes. The topologies of the plastome, mitochondrial rDNA and nuclear rDNA trees were generally congruent, though with discordance at the position of T. chinensis. Several variable regions in the plastomes were identified that have sufficient numbers of parsimony informative sites as to recover the major clades seen in the complete plastome tree. Instead of generating complete plastome sequences, our study showed that accD alone or the concatenation of accD and rbcL can be used in future studies to facilitate identification of Taxillus samples and to generate a molecular phylogeny with robust sampling within the genus.

Genome Sequence of an Okra Leaf Curl Virus from Egypt

A complete Okra leaf curl virus DNA-A was sequenced from okra in Egypt. Here, we report the complete genome sequence of this monopartite virus, comprising 2,764 bp and encoding 6 open reading frames (ORFs) with a GC content of 44.6% and 88.3% similarity to a virus reported earlier from Cameroon.

Improved genome assembly and pan-genome provide key insights into eggplant domestication and breeding

Eggplant (Solanum melongena L.) is an important horticultural crop and one of the most widely grown vegetables from the Solanaceae family. It was domesticated from a wild, prickly progenitor carrying small, round, non-anthocyanic fruits. We obtained a novel, highly contiguous genome assembly of the eggplant '67/3' reference line, by Hi-C retrofitting of a previously released short read- and optical mapping-based assembly. The sizes of the 12 chromosomes and the fraction of anchored genes in the improved assembly were comparable to those of a chromosome-level assembly. We resequenced 23 accessions of S. melongena representative of the worldwide phenotypic, geographic, and genetic diversity of the species, and one each from the closely related species Solanum insanum and Solanum incanum. The eggplant pan-genome contained approximately 51.5 additional megabases and 816 additional genes compared with the reference genome, while the pan-plastome showed little genetic variation. We identified 53 selective sweeps related to fruit color, prickliness, and fruit shape in the nuclear genome, highlighting selection leading to the emergence of present-day S. melongena cultivars from its wild ancestors. Candidate genes underlying the selective sweeps included a MYBL1 repressor and CHALCONE ISOMERASE (for fruit color), homologs of Arabidopsis GLABRA1 and GLABROUS INFLORESCENCE STEMS2 (for prickliness), and orthologs of tomato FW2.2, OVATE, LOCULE NUMBER/WUSCHEL, SUPPRESSOR OF OVATE, and CELL SIZE REGULATOR (for fruit size/shape), further suggesting that selection for the latter trait relied on a common set of orthologous genes in tomato and eggplant.

Phylogenetics and comparative plastome genomics of two of the largest genera of angiosperms, Piper and Peperomia (Piperaceae)

Biological radiations provide unique opportunities to understand the evolution of biodiversity. One such radiation is the pepper plant family Piperaceae, an early-diverging and mega-diverse lineage that could serve as a model to study the diversification of angiosperms. However, traditional genetic markers lack sufficient variation for such studies, and testing hypotheses on poorly resolved phylogenetic frameworks becomes challenging. Limited genomic data is available for Piperaceae, which contains two of the largest genera of angiosperms, Piper (>2100 species) and Peperomia (>1300 species). To address this gap, we used genome skimming to assemble and annotate whole plastomes (152-161kbp) and >5kbp nuclear ribosomal DNA region from representatives of Piper and Peperomia. We conducted phylogenetic and comparative genomic analyses to study plastome evolution and investigate the role of hybridization in this group. Plastome phylogenetic trees were well resolved and highly supported, with a hard incongruence observed between plastome and nuclear phylogenetic trees suggesting hybridization in Piper. While all plastomes of Piper and Peperomia had the same gene content and order, there were informative structural differences between them. First, ycf1 was more variable and longer in Piper than Peperomia, extending well into the small single copy region by thousands of base pairs. We also discovered previously unknown structural variation in 14 out of 25 Piper taxa, tandem duplication of the trnH-GUG gene resulting in an expanded large single copy region. Other early-diverging angiosperms have a duplicated trnH-GUG, but the specific rearrangement we found is unique to Piper and serves to refine knowledge of relationships among early-diverging angiosperms. Our study demonstrates that genome skimming is an efficient approach to produce plastome assemblies for comparative genomics and robust phylogenies of species-rich plant genera.

Chloroplast genomes in Populus (Salicaceae): comparisons from an intensively sampled genus reveal dynamic patterns of evolution

The chloroplast is one of two organelles containing a separate genome that codes for essential and distinct cellular functions such as photosynthesis. Given the importance of chloroplasts in plant metabolism, the genomic architecture and gene content have been strongly conserved through long periods of time and as such are useful molecular tools for evolutionary inferences. At present, complete chloroplast genomes from over 4000 species have been deposited into publicly accessible databases. Despite the large number of complete chloroplast genomes, comprehensive analyses regarding genome architecture and gene content have not been conducted for many lineages with complete species sampling. In this study, we employed the genus Populus to assess how more comprehensively sampled chloroplast genome analyses can be used in understanding chloroplast evolution in a broadly studied lineage of angiosperms. We conducted comparative analyses across Populus in order to elucidate variation in key genome features such as genome size, gene number, gene content, repeat type and number, SSR (Simple Sequence Repeat) abundance, and boundary positioning between the four main units of the genome. We found that some genome annotations were variable across the genus owing in part from errors in assembly or data checking and from this provided corrected annotations. We also employed complete chloroplast genomes for phylogenetic analyses including the dating of divergence times throughout the genus. Lastly, we utilized re-sequencing data to describe the variations of pan-chloroplast genomes at the population level for P. euphratica. The analyses used in this paper provide a blueprint for the types of analyses that can be conducted with publicly available chloroplast genomes as well as methods for building upon existing datasets to improve evolutionary inference.

Comparative analysis of chloroplast genomes of cultivars and wild species of sweetpotato (Ipomoea batatas [L.] Lam)

Background

Sweetpotato (Ipomoea batatas [L.] Lam.) is an important food crop. However, the genetic information of the nuclear genome of this species is difficult to determine accurately because of its large genome and complex genetic background. This drawback has limited studies on the origin, evolution, genetic diversity and other relevant studies on sweetpotato.

Results

The chloroplast genomes of 107 sweetpotato cultivars were sequenced, assembled and annotated. The resulting chloroplast genomes were comparatively analysed with the published chloroplast genomes of wild species of sweetpotato. High similarity and certain specificity were found among the chloroplast genomes of Ipomoea spp. Phylogenetic analysis could clearly distinguish wild species from cultivars. Ipomoea trifida and Ipomoea tabascana showed the closest relationship with the cultivars, and different haplotypes of ycf1 could be used to distinguish the cultivars from their wild relatives. The genetic structure was analyzed using variations in the chloroplast genome. Compared with traditional nuclear markers, the chloroplast markers designed based on the InDels on the chloroplast genome showed significant advantages.

Conclusions

Comparative analysis of chloroplast genomes of 107 cultivars and several wild species of sweetpotato was performed to help analyze the evolution, genetic structure and the development of chloroplast DNA markers of sweetpotato.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-021-07544-y.

Integrative Analysis of the Metabolome and Transcriptome of a Cultivated Pepper and Its Wild Progenitor Chiltepin (Capsicum annuum L. var. glabriusculum) Revealed the Loss of Pungency During Capsicum Domestication

Pungency is a unique characteristic of chili peppers (Capsicum spp.) caused by capsaicinoids. The evolutionary emergence of pungency is thought to be a derived trait within the genus Capsicum. However, it is not well-known how pungency has varied during Capsicum domestication and specialization. In this study, we applied a comparative metabolomics along with transcriptomics analysis to assess various changes between two peppers (a mildly pungent cultivated pepper BB3 and its hot progenitor chiltepin) at four stages of fruit development, focusing on pungency variation. A total of 558 metabolites were detected in two peppers. In comparison with chiltepin, capsaicinoid accumulation in BB3 was almost negligible at the early stage. Next, 412 DEGs associated with the capsaicinoid accumulation pathway were identified through coexpression analysis, of which 18 genes (14 TFs, 3 CBGs, and 1 UGT) were deemed key regulators due to their high coefficients. Based on these data, we speculated that downregulation of these hub genes during the early fruit developmental stage leads to a loss in pungency during Capsicum domestication (from chiltepin to BB3). Of note, a putative UDP-glycosyltransferase, GT86A1, is thought to affect the stabilization of capsaicinoids. Our results lay the foundation for further research on the genetic diversity of pungency traits during Capsicum domestication and specialization.

Plastid genomes of the North American Rhus integrifolia-ovata complex and phylogenomic implications of inverted repeat structural evolution in Rhus L

Plastid genomes (plastomes) represent rich sources of information for phylogenomics, from higher-level studies to below the species level. The genus Rhus (sumac) has received a significant amount of study from phylogenetic and biogeographic perspectives, but genomic studies in this genus are lacking. Rhus integrifolia and R. ovata are two shrubby species of high ecological importance in the southwestern USA and Mexico, where they occupy coastal scrub and chaparral habitats. They hybridize frequently, representing a fascinating system in which to investigate the opposing effects of hybridization and divergent selection, yet are poorly characterized from a genomic perspective. In this study, complete plastid genomes were sequenced for one accession of R. integrifolia and one each of R. ovata from California and Arizona. Sequence variation among these three accessions was characterized, and PCR primers potentially useful in phylogeographic studies were designed. Phylogenomic analyses were conducted based on a robustly supported phylogenetic framework based on 52 complete plastomes across the order Sapindales. Repeat content, rather than the size of the inverted repeat, had a stronger relative association with total plastome length across Sapindales when analyzed with phylogenetic least squares regression. Variation at the inverted repeat boundary within Rhus was striking, resulting in major shifts and independent gene losses. Specifically, rps19 was lost independently in the R. integrifolia-ovata complex and in R. chinensis, with a further loss of rps22 and a major contraction of the inverted repeat in two accessions of the latter. Rhus represents a promising novel system to study plastome structural variation of photosynthetic angiosperms at and below the species level.

Development of nuclear SSR and chloroplast genome markers in diverse Liriodendron chinense germplasm based on low-coverage whole genome sequencing

Background

Liriodendron chinense ranges widely in subtropical China and northern Vietnam; however, it inhabits several small, isolated populations and is now an endangered species due to its limited seed production. The objective of this study was to develop a set of nuclear SSR (simple sequence repeats) and multiple chloroplast genome markers for genetic studies in L. chinense and their characterization in diverse germplasm.

Results

We performed low-coverage whole genome sequencing of the L. chinense from four genotypes, assembled the chloroplast genome and identified nuclear SSR loci by searching in contigs for SSR motifs. Comparative analysis of the four chloroplast genomes of L. chinense revealed 45 SNPs, 17 indels, 49 polymorphic SSR loci, and five small inversions. Most chloroplast intraspecific polymorphisms were located in the interspaces of single-copy regions. In total, 6147 SSR markers were isolated from low-coverage whole genome sequences. The most common SSR motifs were dinucleotide (70.09%), followed by trinucleotide motifs (23.10%). The motif AG/TC (33.51%) was the most abundant, followed by TC/AG (25.53%). A set of 13 SSR primer combinations were tested for amplification and their ability to detect polymorphisms in a set of 109 L. chinense individuals, representing distinct varieties or germplasm. The number of alleles per locus ranged from 8 to 28 with an average of 21 alleles. The expected heterozygosity (H_e) varied from 0.19 to 0.93 and the observed heterozygosity (H_o) ranged from 0.11 to 0.79.

Conclusions

The genetic resources characterized and tested in this study provide a valuable tool to detect polymorphisms in L. chinense for future genetic studies and breeding programs.

The Increase of Simple Sequence Repeats during Diversification of Marchantiidae, An Early Land Plant Lineage, Leads to the First Known Expansion of Inverted Repeats in the Evolutionarily-Stable Structure of Liverwort Plastomes

The chloroplast genomes of liverworts, an early land plant lineage, exhibit stable structure and gene content, however the known resources are very limited. The newly sequenced plastomes of Conocephalum, Riccia and Sphaerocarpos species revealed an increase of simple sequence repeats during the diversification of complex thalloid liverwort lineage. The presence of long TA motifs forced applying the long-read nanopore sequencing method for proper and dependable plastome assembly, since the length of dinucleotide repeats overcome the length of Illumina short reads. The accumulation of SSRs (simple sequence repeats) enabled the expansion of inverted repeats by the incorporation of rps12 and rps7 genes, which were part of large single copy (LSC) regions in the previously sequenced plastomes. The expansion of inverted repeat (IR) at the genus level is reported for the first time for non-flowering plants. Moreover, comparative analyses with remaining liverwort lineages revealed that the presence of SSR in plastomes is specific for simple thalloid species. Phylogenomic analysis resulted in trees confirming monophyly of Marchantiidae and partially congruent with previous studies, due to dataset-dependent results of Dumortiera-Reboulia relationships. Despite the lower evolutionary rate of Marchantiales plastomes, significant barcoding gap was detected, even for recently divergent holarctic Conocephalum species. The sliding window analyses revealed the presence of 18 optimal (500 bp long) barcodes that enable the molecular identification of all studied species.

The complete mitochondrial genome of the chiltepin pepper (Capsicum annuum var. glabriusculum), the wild progenitor of Capsicum annuum L

The complete mitochondrial genome of chiltepin pepper (Capsicum annuum var glabriusculum) was sequenced. The mitogenome of the American bird pepper was 505,190 bp, with 44.4% of GC content. A total of 218 genes were fully annotated, including 190 CDS (31 known genes and 158 open reading frames), three rRNA, and 25 tRNA genes. The gene synteny and number were equal to those of C. annuum var annuum, except for the partial annotation of ATP subunit 6 and the absence of ORF172 and ORF104b. The complete mt genome sequence was deposited to the GenBank (NCBI, Accession number: MN196478).