Species of the genus Nicotiana (Solanaceae), commonly referred to as tobacco plants, are often cultivated as non-food crops and garden ornamentals. In addition to the worldwide production of tobacco leaves, they are also used as evolutionary model systems due to their complex development history tangled by polyploidy and hybridization. Here, we assembled the plastid genomes of five tobacco species: N. knightiana, N. rustica, N. paniculata, N. obtusifolia and N. glauca. De novo assembled tobacco plastid genomes had the typical quadripartite structure, consisting of a pair of inverted repeat (IR) regions (25,323–25,369 bp each) separated by a large single-copy (LSC) region (86,510–86,716 bp) and a small single-copy (SSC) region (18,441–18,555 bp). Comparative analyses of Nicotiana plastid genomes with currently available Solanaceae genome sequences showed similar GC and gene content, codon usage, simple sequence and oligonucleotide repeats, RNA editing sites, and substitutions. We identified 20 highly polymorphic regions, mostly belonging to intergenic spacer regions (IGS), which could be suitable for the development of robust and cost-effective markers for inferring the phylogeny of the genus Nicotiana and family Solanaceae. Our comparative plastid genome analysis revealed that the maternal parent of the tetraploid N. rustica was the common ancestor of N. paniculata and N. knightiana, and the later species is more closely related to N. rustica. Relaxed molecular clock analyses estimated the speciation event between N. rustica and N. knightiana appeared 0.56 Ma (HPD 0.65–0.46). Biogeographical analysis supported a south-to-north range expansion and diversification for N. rustica and related species, where N. undulata and N. paniculata evolved in North/Central Peru, while N. rustica developed in Southern Peru and separated from N. knightiana, which adapted to the Southern coastal climatic regimes. We further inspected selective pressure on protein-coding genes among tobacco species to determine if this adaptation process affected the evolution of plastid genes. These analyses indicate that four genes involved in different plastid functions, including DNA replication (rpoA) and photosynthesis (atpB, ndhD and ndhF), came under positive selective pressure as a result of specific environmental conditions. Genetic mutations in these genes might have contributed to better survival and superior adaptations during the evolutionary history of tobacco species.

The genus Nicotiana of the family Solanaceae, commonly referred to as tobacco plants, are a group cultivated as garden ornamentals. Besides their use in the worldwide production of tobacco leaves, they are also used as evolutionary model systems due to their complex development history, which is tangled by polyploidy and hybridization. Here, we assembled the plastid genomes of five tobacco species, namely N. knightiana, N. rustica, N. paniculata, N. obtusifolia and N. glauca. De novo assembled tobacco plastid genomes showed typical quadripartite structure, consisting of a pair of inverted repeats (IR) regions (25,323–25,369 bp each) separated by a large single copy (LSC) region (86,510 –86,716 bp) and a small single copy (SSC) region (18,441–18,555 bp). Comparative analyses of Nicotiana plastid genomes showed similar GC content, gene content, codon usage, simple sequence repeats, oligonucleotide repeats, RNA editing sites and substitutions with currently available Solanaceae genomes sequences. We identified twenty highly polymorphic regions mostly belonging to intergenic spacer regions (IGS), which could be appropriate for the development of robust and cost-effective markers to infer the phylogeny of genus Nicotiana and family Solanaceae. Our comparative plastid genome analysis revealed that the maternal parent of the tetraploid N. rustica was the common ancestor of N. paniculata and N. knightiana, and the later species is more closely related to N. rustica. The relaxed molecular clock analyses estimated that the speciation event between N. rustica and knightiana appeared 0.56 Ma (HPD 0.65–0.46). The biogeographical analysis showed a south-to-north range expansion and diversification for N. rustica and related species, where N. undulata and N. paniculata evolved in North/Central Peru, while N. rustica developed in Southern Peru and separated from N. knightiana, which adapted to the Southern coastal climatic regimes. We further inspected selective pressure on protein-coding genes among tobacco species to determine if this adaptation process affected the evolution of plastid genes. These analyses indicated that four genes involved in different plastid functions, such as DNA replication (rpoA) and photosynthesis (atpB, ndhD and ndhF), came under positive selective pressure as a result of specific environmental conditions. Genetic mutations of the following genes might have contributed to the survival and better adaptation during the evolutionary history of tobacco species.

Within the family Solanaceae, Withania is a small genus belonging to the Solanoideae subfamily. Here, we report the de novo assembled chloroplast genome sequences of W. coagulans, W. adpressa, and W. riebeckii. The length of these genomes ranged from 154,162 to 154,364 base pairs (bp). These genomes contained a pair of inverted repeats (IRa and IRb) ranging from 25,029 to 25,071 bp that were separated by a large single-copy (LSC) region of 85,635–85,765 bp and a small single-copy (SSC) region of 18,457–18,469 bp. We analyzed the structural organization, gene content and order, guanine-cytosine content, codon usage, RNA-editing sites, microsatellites, oligonucleotide and tandem repeats, and substitutions of Withania plastomes, which revealed high similarities among the species. Comparative analysis among the Withania species also highlighted 10 divergent hotspots that could potentially be used for molecular marker development, phylogenetic analysis, and species identification. Furthermore, our analyses showed that even three mutational hotspots (rps4-trnT, trnM-atpE, and rps15) were sufficient to discriminate the Withania species included in current study.

Within the family Solanaceae, Withania is a small genus belonging to the Solanoideae subfamily. Here, we report the de novo assembled, complete, plastomed genome sequences of W. coagulans, W. adpressa, and W. riebeckii. The length of these genomes ranged from 154,198 base pairs (bp) to 154,361 bp and contained a pair of inverted repeats (IRa and IRb) of 25,027--25,071 bp that were separated by a large single-copy (LSC) region of 85,675--85,760 bp and a small single-copy (SSC) region of 18,457--18,469 bp. We analyzed the structural organization, gene content and order, guanine-cytosine content, codon usage, RNA-editing sites, microsatellites, oligonucleotide and tandem repeats, and substitutions of Withania plastid genomes, which revealed close resemblance among the species. Both the substitution and insertion and deletion analyses confirmed that the IR region was significantly conserved compared with the LSC and SSC regions. Further comparative analysis among the Withania species highlighted 30 divergent hotspots that could potentially be used for molecular marker development, phylogenetic analysis, and species identification.