The chloroplast genome of Hyoscyamus niger and a phylogenetic study of the tribe Hyoscyameae (Solanaceae)

Characterization of the complete chloroplast genome of the rare medicinal plant: Mandragora caulescens (Solanaceae)

In this study, we assembled high-quality chloroplast genomes of Mandragora caulescens through a reference-guided approach using high-throughput Illumina sequencing reads. The resulting chloroplast genome assembly displayed a typical quadripartite structural organization, comprising a large single-copy (LSC) region of 85,233 bp, two inverted repeat (IR) regions of 25,685 bp each, and a small single-copy (SSC) region of 18,207 bp. The chloroplast genome harbored 141 complete genes, and its overall GC content was 38.0%. In maximum-likelihood (ML) and Bayesian inference (BI) trees, the 19 Solanaceae species formed a monophyletic group, dividing into two main clades. M. caulescens and Nicandra physalodes formed a monophyletic group, suggesting a close relationship between the two species. The M. caulescens cp genome presented in this study lays a good foundation for further genetic and genomic studies of the Solanaceae.

Phylogenetic analysis and divergence time estimation of Lycium species in China based on the chloroplast genomes

BACKGROUND

Lycium is an economically and ecologically important genus of shrubs, consisting of approximately 70 species distributed worldwide, 15 of which are located in China. Despite the economic and ecological importance of Lycium, its phylogeny, interspecific relationships, and evolutionary history remain relatively unknown. In this study, we constructed a phylogeny and estimated divergence time based on the chloroplast genomes (CPGs) of 15 species, including subspecies, of the genus Lycium from China.

RESULTS

We sequenced and annotated 15 CPGs in this study. Comparative analysis of these genomes from these Lycium species revealed a typical quadripartite structure, with a total sequence length ranging from 154,890 to 155,677 base pairs (bp). The CPGs was highly conserved and moderately differentiated. Through annotation, we identified a total of 128-132 genes. Analysis of the boundaries of inverted repeat (IR) regions showed consistent positioning: the junctions of the IRb/LSC region were located in rps19 in all Lycium species, IRb/SSC between the ycf1 and ndhF genes, and SSC/IRa within the ycf1 gene. Sequence variation in the SSC region exceeded that in the IR region. We did not detect major expansions or contractions in the IR region or rearrangements or insertions in the CPGs of the 15 Lycium species. Comparative analyses revealed five hotspot regions in the CPG: trnR(UCU), atpF-atpH, ycf3-trnS(GGA), trnS(GGA), and trnL-UAG, which could potentially serve as molecular markers. In addition, phylogenetic tree construction based on the CPG indicated that the 15 Lycium species formed a monophyletic group and were divided into two typical subbranches and three minor branches. Molecular dating suggested that Lycium diverged from its sister genus approximately 17.7 million years ago (Mya) and species diversification within the Lycium species of China primarily occurred during the recent Pliocene epoch.

CONCLUSION

The divergence time estimation presented in this study will facilitate future research on Lycium, aid in species differentiation, and facilitate diverse investigations into this economically and ecologically important genus.

The complete chloroplast genome sequence of Castanopsis fordii Hance (Fagaceae)

Castanopsis fordii Hance 1884 is a typical evergreen broad-leaved forest plant in the south subtropical and middle subtropical regions of China. It has high utilization value in wood production and soil erosion protection. Here, we first reported and characterized the complete chloroplast (cp) genome sequence of C. fordii based on Illumina paired-end sequencing data. The complete cp genome sequence of C. fordii was 160,853 base pairs (bp) in length which contained two inverted repeats (IRs) of 25,699 bp separated by a large single-copy (LSC) and a small single copy (SSC) of 90,474 bp and 18,981 bp, respectively. The cpDNA contained 129 genes, comprising 85 protein-coding genes, 36 tRNA genes, 8 rRNA genes. The overall GC content of the plastome was 36.8%. Phylogenetic analysis base on 14 chloroplast genomes indicated that C. fordii was closely related to the species C. tibetana and C. concinna in Fagaceae.

Complete Plastome of Physalis angulata var. villosa, Gene Organization, Comparative Genomics and Phylogenetic Relationships among Solanaceae

Physalis angulata var. villosa, rich in withanolides, has been used as a traditional Chinese medicine for many years. To date, few extensive molecular studies of this plant have been conducted. In the present study, the plastome of P. angulata var. villosa was sequenced, characterized and compared with that of other Physalis species, and a phylogenetic analysis was conducted in the family Solanaceae. The plastome of P. angulata var. villosa was 156,898 bp in length with a GC content of 37.52%, and exhibited a quadripartite structure typical of land plants, consisting of a large single-copy (LSC, 87,108 bp) region, a small single-copy (SSC, 18,462 bp) region and a pair of inverted repeats (IR: IRA and IRB, 25,664 bp each). The plastome contained 131 genes, of which 114 were unique and 17 were duplicated in IR regions. The genome consisted of 85 protein-coding genes, eight rRNA genes and 38 tRNA genes. A total of 38 long, repeat sequences of three types were identified in the plastome, of which forward repeats had the highest frequency. Simple sequence repeats (SSRs) analysis revealed a total of 57 SSRs, of which the T mononucleotide constituted the majority, with most of SSRs being located in the intergenic spacer regions. Comparative genomic analysis among nine Physalis species revealed that the single-copy regions were less conserved than the pair of inverted repeats, with most of the variation being found in the intergenic spacer regions rather than in the coding regions. Phylogenetic analysis indicated a close relationship between Physalis and Withania. In addition, Iochroma, Dunalia, Saracha and Eriolarynx were paraphyletic, and clustered together in the phylogenetic tree. Our study published the first sequence and assembly of the plastome of P. angulata var. villosa, reported its basic resources for evolutionary studies and provided an important tool for evaluating the phylogenetic relationship within the family Solanaceae.

Characterization of the Evolutionary Pressure on Anisodus tanguticus Maxim. with Complete Chloroplast Genome Sequence

Anisodus tanguticus Maxim. (Solanaceae), a traditional endangered Tibetan herb, is endemic to the Qinghai-Tibet Plateau. Here, we report the de novo assembled chloroplast (cp) genome sequences of A. tanguticus (155,765 bp). The cp contains a pair of inverted repeated (IRa and IRb) regions of 25,881 bp that are separated by a large single copy (LSC) region (86,516 bp) and a small single copy SSC (17,487 bp) region. A total of 132 functional genes were annotated in the cp genome, including 87 protein-coding genes, 37 tRNA genes, and 8 rRNA genes. Moreover, 199 simple sequence repeats (SSR) and 65 repeat structures were detected. Comparative plastome analyses revealed a conserved gene order and high similarity of protein-coding sequences. The A. tanguticus cp genome exhibits contraction and expansion, which differs from Przewalskia tangutica and other related Solanaceae species. We identified 30 highly polymorphic regions, mostly belonging to intergenic spacer regions (IGS), which may be suitable for the development of robust and cost-effective markers for inferring the phylogeny of the genus Anisodus and family Solanaceae. Analysis of the Ka/Ks ratios of the Hyoscyameae tribe revealed significant positive selection exerted on the cemA, rpoC2, and clpP genes, which suggests that protein metabolism may be an important strategy for A. tanguticus and other species in Hyoscyameae in adapting to the adverse environment on the Qinghai-Tibetan Plateau. Phylogenetic analysis revealed that A. tanguticus clustered closer with Hyoscyamus niger than P. tangutica. Our results provide reliable genetic information for future exploration of the taxonomy and phylogenetic evolution of the Hyoscyameae tribe and related species.

The first checklist of alien vascular plants of Kyrgyzstan, with new records and critical evaluation of earlier data. Contribution 2

Background

We continue the inventory of alien vascular plants of Kyrgyzstan, with emphasis on the time and pathways of introduction of the species and their current status in the territory. Each taxon is discussed in the context of plant invasions in Central Asia. This work is a further development of the preliminary checklist of alien plants of Kyrgyzstan, which was compiled for the Global Register of Introduced and Invasive Species in 2018.

New information

This contribution includes all alien species of Kyrgyzstan belonging to Solanaceae and Asphodelaceae and one species of Asteraceae. Physalisphiladelphicus (syn. P.ixocarpa) is reported for the first time from Central Asia, as new to Kazakhstan, Kyrgyzstan and Uzbekistan, thus marking a recent invasion with a variety of imported grain and seed material. The old records of P.ixocarpa from Uzbekistan are based on misidentified specimens of P.angulata. Physalisangulata is an old cotton immigrant in Central Asia, whose invasion started in the 1920s; it is excluded from the alien flora of Kyrgyzstan as registered in error on the basis of cultivated plants. Alkekengiofficinarum is an archaeophyte of the Neolithic period in Central Asia, formerly used for food, now strongly declining and largely casual in Kyrgyzstan. The only historical record of Physalisviscosa from Uzbekistan was based on a technical error and belongs to A.officinarum. Daturastramonium and Hyoscyamusniger were introduced as medicinal plants during the period of the Arabic invasion of Central Asia, by the 11th century. Daturainnoxia is a newly recorded casual alien, recently escaped from ornamental cultivation. Nicandraphysalodes is a casual alien, which was cultivated by Russian colonists in the early 20^th century for culinary use and is currently used in ornamental cultivation. Hemerocallisfulva was a remnant of historical cultivation in the former Khanate of Buxoro, and its formerly established colonies are presumably extinct in the wild. Bidensfrondosa was seemingly introduced with contaminated forage and seed of American origin during the late Soviet period and started to spread in the period of independence; its invasion in the former USSR is analysed.

Plastid phylogenomics and biogeography of the medicinal plant lineage Hyoscyameae (Solanaceae)

The cosmopolitan family Solanaceae, which originated and first diversified in South America, is economically important. The tribe Hyoscyameae is one of the three clades in Solanaceae that occurs outside of the New World; Hyoscyameae genera are distributed mainly in Europe and Asia, and have centers of species diversity in the Qinghai-Tibet Plateau and adjacent regions. Although many phylogenetic studies have focused on Solanaceae, the phylogenetic relationships within the tribe Hyoscyameae and its biogeographic history remain obscure. In this study, we reconstructed the phylogeny of Hyoscyameae based on whole chloroplast genome data, and estimated lineage divergence times according to the newly reported fruit fossil from the Eocene Patagonia, Physalis infinemundi, the earliest known fossil of Solanaceae. We reconstructed a robust phylogeny of Hyoscyameae that reveals the berry fruit-type Atropa is sister to the six capsule-bearing genera (Hyoscyameae sensu stricto), Atropanthe is sister to the clade (Scopolia, Physochlaina, Przewalskia), and together they are sister to the robustly supported Anisodus-Hyoscyamus clade. The stem age of Hyoscyameae was inferred to be in the Eocene (47.11 Ma, 95% HPD: 36.75-57.86 Ma), and the crown ages of Hyoscyameae sensu stricto were estimated as the early Miocene (22.52 Ma, 95% HPD: 15.19-30.53 Ma), which shows a close correlation with the rapid uplift of the Qinghai-Tibet Plateau at the Paleogene/Neogene boundary. Our results provide insights into the phylogenetic relationships and the history of the biogeographic diversification of the tribe Hyoscyameae, as well as plant diversification on the Qinghai-Tibet Plateau.

Polyploidization is accompanied by synonymous codon usage bias in the chloroplast genomes of both cotton and wheat

Synonymous codon usage bias (SCUB) of both nuclear and organellar genes can mirror the evolutionary specialization of plants. The polyploidization process exposes the nucleus to genomic shock, a syndrome which promotes, among other genetic variants, SCUB. Its effect on organellar genes has not, however, been widely addressed. The present analysis targeted the chloroplast genomes of two leading polyploid crop species, namely cotton and bread wheat. The frequency of codons in the chloroplast genomes ending in either adenosine (NNA) or thymine (NNT) proved to be higher than those ending in either guanidine or cytosine (NNG or NNC), and this difference was conserved when comparisons were made between polyploid and diploid forms in both the cotton and wheat taxa. Preference for NNA/T codons was heterogeneous among genes with various numbers of introns and was also differential among the exons. SCUB patterns distinguished tetraploid cotton from its diploid progenitor species, as well as bread wheat from its diploid/tetraploid progenitor species, indicating that SCUB in the chloroplast genome partially mirrors the formation of polyploidies.

Comparing and phylogenetic analysis chloroplast genome of three Achyranthes species

In this study, the chloroplast genome sequencing of the Achyranthes longifolia, Achyranthes bidentata and Achyranthes aspera were performed by Next-generation sequencing technology. The results revealed that there were a length of 151,520 bp (A. longifolia), 151,284 bp (A. bidentata), 151,486 bp (A. aspera), respectively. These chloroplast genome have a highly conserved structure with a pair of inverted repeat (IR) regions (25,150 bp; 25,145 bp; 25,150 bp), a large single copy (LSC) regions (83,732 bp; 83,933 bp; 83,966 bp) and a small single copy (SSC) regions (17,252 bp; 17,263 bp; 17,254 bp) in A. bidentate, A. aspera and A. longifolia. There were 127 genes were annotated, which including 8 rRNA genes, 37 tRNA genes and 82 functional genes. The phylogenetic analysis strongly revealed that Achyranthes is monophyletic, and A. bidentata was the closest relationship with A. aspera and A. longifolia. A. bidentata and A. longifolia were clustered together, the three Achyranthes species had the same origin, then the gunes of Achyranthes is the closest relative to Alternanthera, and that forms a group with Alternanthera philoxeroides. The research laid a foundation and provided relevant basis for the identification of germplasm resources in the future.

RNA editing in the chloroplast of Asian Palmyra palm (Borassus flabellifer)

We have identified 46 RNA editing sites located in 20 chloroplast (cp) genes of Borassus flabellifer (Asian Palmyra palm), family Arecaceae, and tested these genes for supporting phylogenetic study among the commelinids. Among the 46 sites, 43 sites were found to cause amino acid alterations, which were predicted to increase the hydrophobicity and transmembrane regions of the proteins, and one site was to cause a premature stop codon. Analysis of these editing sites with data obtained from seed plants showed that a number of shared-editing sites depend on the evolutionary relationship between plants. We reconstructed a deep phylogenetic relationship among the commelinids using seven RNA edited genes that are orthologous among monocots. This tree could represent the relationship among subfamilies of Arecaceae family, but was insufficient to represent the relationship among the orders of the commelinid. After adding eight gene sequences with high parsimony-informative characters (PICs), the tree topology was improved and could support the topology for the commelinid orders ((Arecales,Dasypogenaceae) (Zingiberales+Commelinales,Poales)). The result provides support for inherent RNA editing along the evolution of seed plants, and we provide an alternative set of loci for the phylogenetic tree reconstruction of Arecaceae's subfamilies.

The complete chloroplast genome of Physochlaina physaloides (Solanaceae), an important medicinal plant

Physochlaina is an important perennial herbaceous genus with significant medicinal value, while the phylogeny of Physochlaina and tribe Hyoscyameae is not well resolved yet. In this study, we report the complete chloroplast genome sequences of Ph. physaloides, its complete chloroplast genome is 156,413 bp in length, which is a typical quadripartite structure that includes a large single-copy region of 86,659 bp, a small single-copy region of 18,012 bp, and its GC content was 37.7%. A total of 132 genes were identified, including 87 protein-coding genes, 37 tRNA genes, and 8 rRNA genes. Furthermore, a phylogenetic tree of the tribe Hyoscyameae was constructed based the complete chloroplast genome sequence, and a new topology of the tribe was obtained. This study provides valuable genetic information for the conservation and utilization of Ph. physaloides and also provide the potential for better understanding of the phylogeny of Hyoscyameae and Solanaceae.

Sequencing and Structural Analysis of the Complete Chloroplast Genome of the Medicinal Plant Lycium chinense Mill

Lycium chinense Mill, an important Chinese herbal medicine, is widely used as a dietary supplement and food. Here the chloroplast (CP) genome of L. chinense was sequenced and analyzed, revealing a size of 155,756 bp and with a 37.8% GC content. The L. chinense CP genome comprises a large single copy region (LSC) of 86,595 bp and a small single copy region (SSC) of 18,209 bp, and two inverted repeat regions (IRa and IRb) of 25,476 bp separated by the single copy regions. The genome encodes 114 genes, 16 of which are duplicated. Most of the 85 protein-coding genes (CDS) had standard ATG start codons, while 3 genes including rps12, psbL and ndhD had abnormal start codons (ACT and ACG). In addition, a strong A/T bias was found in the majority of simple sequence repeats (SSRs) detected in the CP genome. Analysis of the phylogenetic relationships among 16 species revealed that L. chinense is a sister taxon to Lycium barbarum. Overall, the complete sequence and annotation of the L. chinense CP genome provides valuable genetic information to facilitate precise understanding of the taxonomy, species and phylogenetic evolution of the Solanaceae family.

Sequencing and Comparative Analysis of the Chloroplast Genome of Angelica polymorpha and the Development of a Novel Indel Marker for Species Identification

The genus Angelica (Apiaceae) comprises valuable herbal medicines. In this study, we determined the complete chloroplast (CP) genome sequence of A. polymorpha and compared it with that of Ligusticum officinale (GenBank accession no. NC039760). The CP genomes of A. polymorpha and L. officinale were 148,430 and 147,127 bp in length, respectively, with 37.6% GC content. Both CP genomes harbored 113 unique functional genes, including 79 protein-coding, four rRNA, and 30 tRNA genes. Comparative analysis of the two CP genomes revealed conserved genome structure, gene content, and gene order. However, highly variable regions, sufficient to distinguish between A. polymorpha and L. officinale, were identified in hypothetical chloroplast open reading frame1 (ycf1) and ycf2 genic regions. Nucleotide diversity (Pi) analysis indicated that ycf4⁻chloroplast envelope membrane protein (cemA) intergenic region was highly variable between the two species. Phylogenetic analysis revealed that A. polymorpha and L. officinale were well clustered at family Apiaceae. The ycf4-cemA intergenic region in A. polymorpha carried a 418 bp deletion compared with L. officinale. This region was used for the development of a novel indel marker, LYCE, which successfully discriminated between A. polymorpha and L. officinale accessions. Our results provide important taxonomic and phylogenetic information on herbal medicines and facilitate their authentication using the indel marker.

The chloroplast genome sequence of bittersweet (Solanum dulcamara): Plastid genome structure evolution in Solanaceae

Bittersweet (Solanum dulcamara) is a native Old World member of the nightshade family. This European diploid species can be found from marshlands to high mountainous regions and it is a common weed that serves as an alternative host and source of resistance genes against plant pathogens such as late blight (Phytophthora infestans). We sequenced the complete chloroplast genome of bittersweet, which is 155,580 bp in length and it is characterized by a typical quadripartite structure composed of a large (85,901 bp) and small (18,449 bp) single-copy region interspersed by two identical inverted repeats (25,615 bp). It consists of 112 unique genes from which 81 are protein-coding, 27 tRNA and four rRNA genes. All bittersweet plastid genes including non-functional ones and even intergenic spacer regions are transcribed in primary plastid transcripts covering 95.22% of the genome. These are later substantially edited in a post-transcriptional phase to activate gene functions. By comparing the bittersweet plastid genome with all available Solanaceae sequences we found that gene content and synteny are highly conserved across the family. During genome comparison we have identified several annotation errors, which we have corrected in a manual curation process then we have identified the major plastid genome structural changes in Solanaceae. Interpreted in a phylogenetic context they seem to provide additional support for larger clades. The plastid genome sequence of bittersweet could help to benchmark Solanaceae plastid genome annotations and could be used as a reference for further studies. Such reliable annotations are important for gene diversity calculations, synteny map constructions and assigning partitions for phylogenetic analysis with de novo sequenced plastomes of Solanaceae.

The Complete Chloroplast Genomes of Six Ipomoea Species and Indel Marker Development for the Discrimination of Authentic Pharbitidis Semen (Seeds of I. nil or I. purpurea)

Ipomoea L. is the largest genus within the Convolvulaceae and contains 600-700 species. Ipomoea species (morning glories) are economically valuable as horticultural species and scientifically valuable as ecological model plants to investigate mating systems, molecular evolution, and both plant-herbivore and plant-parasite interactions. Furthermore, the dried seeds of I. nil or I. purpurea are used in Korean traditional herbal medicines. In this study, chloroplast (cp) genomes were sequenced from six Ipomoea species, namely, I. nil and I. purpurea and, for the first time, I. triloba, I. lacunosa, I. hederacea, and I. hederacea var. integriuscula. The cp genomes were 161,354-161,750 bp in length and exhibited conserved quadripartite structures. In total, 112 genes were identified, including 78 protein-coding regions, 30 transfer RNA genes, and 4 ribosomal RNA genes. The gene order, content, and orientation of the six Ipomoea cp genomes were highly conserved and were consistent with the general structure of angiosperm cp genomes. Comparison of the six Ipomoea cp genomes revealed locally divergent regions, mainly within intergenic spacer regions (petN-psbM, trnI-CAU-ycf2, ndhH-ndhF, psbC-trnS, and ccsA-ndhD). In addition, the protein-coding genes accD, cemA, and ycf2 exhibited high sequence variability and were under positive selection (Ka/Ks > 1), indicating adaptive evolution to the environment within the Ipomoea genus. Phylogenetic analysis of the six Ipomoea species revealed that these species clustered according to the APG IV system. In particular, I. nil and I. hederacea had monophyletic positions, with I. purpurea as a sister. I. triloba and I. lacunosa in the section Batatas and I. hederacea and I. hederacea var. integriuscula in the section Quamoclit were supported in this study with strong bootstrap values and posterior probabilities. We uncovered high-resolution phylogenetic relationships between Ipomoeeae. Finally, indel markers (IPOTY and IPOYCF) were developed for the discrimination of the important herbal medicine species I. nil and I. purpurea. The cp genomes and analyses in this study provide useful information for taxonomic, phylogenetic, and evolutionary analysis of the Ipomoea genome, and the indel markers will be useful for authentication of herbal medicines.

The Complete Chloroplast Genome Sequences of Fritillaria ussuriensis Maxim. and Fritillaria cirrhosa D. Don, and Comparative Analysis with Other Fritillaria Species

The genus Fritillaria belongs to the widely distributed Liliaceae. The bulbs of Fritillaria, F. ussuriensis and F. cirrhosa are valuable herbaceous medicinal ingredients. However, they are still used indiscriminately in herbal medicine. Identification and molecular phylogenic analysis of Fritillaria species are therefore required. Here, we report the complete chloroplast (CP) genome sequences of F. ussuriensis and F. cirrhosa. The two Fritillaria CP genomes were 151,524 and 151,083 bp in length, respectively, and each included a pair of inverted repeated regions (52,678 and 52,156 bp) that was separated by a large single copy region (81,732 and 81,390 bp), and a small single copy region (17,114 and 17,537 bp). A total of 111 genes in F. ussuriensis and 112 in F. cirrhosa comprised 77 protein-coding regions in F. ussuriensis and 78 in F. cirrhosa, 30 transfer RNA (tRNA) genes, and four ribosomal RNA (rRNA) genes. The gene order, content, and orientation of the two Fritillaria CP genomes exhibited the general structure of flowering plants, and were similar to those of other Fritillaria species. Comparison of the six Fritillaria species’ CP genomes indicated seven highly divergent regions in intergenic spacers and in the matK, rpoC1, rpoC2, ycf1, ycf2, ndhD, and ndhF coding regions. We established the position of the six species through phylogenic analysis. The complete chloroplast genome sequences of the two Fritillaria species and a comparison study are useful genomic information for identifying and for studying the phylogenetic relationship among Fritillaria species within the Liliaceae.

Chloroplast genomes: diversity, evolution, and applications in genetic engineering

Chloroplasts play a crucial role in sustaining life on earth. The availability of over 800 sequenced chloroplast genomes from a variety of land plants has enhanced our understanding of chloroplast biology, intracellular gene transfer, conservation, diversity, and the genetic basis by which chloroplast transgenes can be engineered to enhance plant agronomic traits or to produce high-value agricultural or biomedical products. In this review, we discuss the impact of chloroplast genome sequences on understanding the origins of economically important cultivated species and changes that have taken place during domestication. We also discuss the potential biotechnological applications of chloroplast genomes.

Synonymous codon usage bias in plant mitochondrial genes is associated with intron number and mirrors species evolution

Synonymous codon usage bias (SCUB) is a common event that a non-uniform usage of codons often occurs in nearly all organisms. We previously found that SCUB is correlated with both intron number and exon position in the plant nuclear genome but not in the plastid genome; SCUB in both nuclear and plastid genome can mirror the evolutionary specialization. However, how about the rules in the mitochondrial genome has not been addressed. Here, we present an analysis of SCUB in the mitochondrial genome, based on 24 plant species ranging from algae to land plants. The frequencies of NNA and NNT (A- and T-ending codons) are higher than those of NNG and NNC, with the strongest preference in bryophytes and the weakest in land plants, suggesting an association between SCUB and plant evolution. The preference for NNA and NNT is more evident in genes harboring a greater number of introns in land plants, but the bias to NNA and NNT exhibits even among exons. The pattern of SCUB in the mitochondrial genome differs in some respects to that present in both the nuclear and plastid genomes.

Synonymous Codon Usage Bias in the Plastid Genome is Unrelated to Gene Structure and Shows Evolutionary Heterogeneity

Synonymous codon usage bias (SCUB) is the nonuniform usage of codons, occurring often in nearly all organisms. Our previous study found that SCUB is correlated with intron number, is unequal among exons in the plant nuclear genome, and mirrors evolutionary specialization. However, whether this rule exists in the plastid genome has not been addressed. Here, we present an analysis of SCUB in the plastid genomes of 25 species from lower to higher plants (algae, bryophytes, pteridophytes, gymnosperms, and spermatophytes). We found NNA and NNT (A- and T-ending codons) are preferential in the plastid genomes of all plants. Interestingly, this preference is heterogeneous among taxonomies of plants, with the strongest preference in bryophytes and the weakest in pteridophytes, suggesting an association between SCUB and plant evolution. In addition, SCUB frequencies are consistent among genes with varied introns and among exons, indicating that the bias of NNA and NNT is unrelated to either intron number or exon position. Further, SCUB is associated with DNA methylation–induced conversion of cytosine to thymine in the vascular plants but not in algae or bryophytes. These data demonstrate that these SCUB profiles in the plastid genome are distinctly different compared with the nuclear genome.

Homologous recombination and retention of a single form of most genes shape the highly chimeric mitochondrial genome of a cybrid plant

The structure and evolution of angiosperm mitochondrial genomes are driven by extremely high rates of recombination and rearrangement. An excellent experimental system for studying these events is offered by cybrid plants, in which parental mitochondria usually fuse and their genomes recombine. Little is known about the extent, nature and consequences of mitochondrial recombination in these plants. We conducted the first study in which the organellar genomes of a cybrid - between Nicotiana tabacum and Hyoscyamus niger - were sequenced and compared to those of its parents. This cybrid mitochondrial genome is highly recombinant, reflecting at least 30 crossovers and five gene conversions between its parental genomes. It is also surprisingly large (41% and 64% larger than the parental genomes), yet contains single alleles for 90% of mitochondrial genes. Recombination produced a remarkably chimeric cybrid mitochondrial genome and occurred entirely via homologous mechanisms involving the double-strand break repair and/or break-induced replication pathways. Retention of a single form of most genes could be advantageous to minimize intracellular incompatibilities and/or reflect neutral forces that preferentially eliminate duplicated regions. We discuss the relevance of these findings to the surprisingly frequent occurrence of horizontal gene - and genome - transfer in angiosperm mitochondrial DNAs.