Complete sequence of the maize chloroplast genome: gene content, hotspots of divergence and fine tuning of genetic information by transcript editing

Removal of the large inverted repeat from the plastid genome reveals gene dosage effects and leads to increased genome copy number

The chloroplast genomes of most plants and algae contain a large inverted repeat (IR) region that separates two single-copy regions and harbours the ribosomal RNA operon. We have addressed the functional importance of the IR region by removing an entire copy of the 25.3-kb IR from the tobacco plastid genome. Using plastid transformation and subsequent selectable marker gene elimination, we precisely excised the IR, thus generating plants with a substantially reduced plastid genome size. We show that the lack of the IR results in a mildly reduced plastid ribosome number, suggesting a gene dosage benefit from the duplicated presence of the ribosomal RNA operon. Moreover, the IR deletion plants contain an increased number of plastid genomes, suggesting that genome copy number is regulated by measuring total plastid DNA content rather than by counting genomes. Together, our findings (1) demonstrate that the IR can enhance the translation capacity of the plastid, (2) reveal the relationship between genome size and genome copy number, and (3) provide a simplified plastid genome structure that will facilitate future synthetic biology applications.

Maize requires Embryo defective27 for embryogenesis and seedling development

The essential role of plastid translation in embryogenesis has been established in many plants, but a retrograde signal triggered by defective plastid translation machinery that may leads to embryogenesis arrest remains unknown. In this study, we characterized an embryo defective27 (emb27) mutant in maize (Zea mays), and cloning indicates that Emb27 encodes the plastid ribosomal protein S13. The null mutant emb27-1 conditions an emb phenotype with arrested embryogenesis; however, the leaky mutant emb27-2 exhibits normal embryogenesis but an albino seedling-lethal phenotype. The emb27-1/emb27-2 trans-heterozygotes display varying phenotypes from emb to normal seeds but albino seedlings. Analysis of the Emb27 transcription levels in these mutants revealed that the Emb27 expression level in the embryo corresponds with the phenotypic expression of the emb27 mutants. In the W22 genetic background, an Emb27 transcription level higher than 6% of the wild-type level renders normal embryogenesis, whereas lower than that arrests embryogenesis. Mutation of Emb27 reduces the level of plastid 16S rRNA and the accumulation of the plastid-encoded proteins. As a secondary effect, splicing of several plastid introns was impaired in emb27-1 and 2 other plastid translation-defective mutants, emb15 and emb16, suggesting that plastome-encoded factors are required for the splicing of these introns, such as Maturase K (MatK). Our results indicate that EMB27 is essential for plastid protein translation, embryogenesis, and seedling development in maize and reveal an expression threshold of Emb27 for maize embryogenesis.

Twelve newly assembled jasmine chloroplast genomes: unveiling genomic diversity, phylogenetic relationships and evolutionary patterns among Oleaceae and Jasminum species

BACKGROUND

Jasmine (Jasminum), renowned for its ornamental value and captivating fragrance, has given rise to numerous species and accessions. However, limited knowledge exists regarding the evolutionary relationships among various Jasminum species.

RESULTS

In the present study, we sequenced seven distinct Jasminum species, resulting in the assembly of twelve high-quality complete chloroplast (cp) genomes. Our findings revealed that the size of the 12 cp genomes ranged from 159 to 165 kb and encoded 134-135 genes, including 86-88 protein-coding genes, 38-40 tRNA genes, and 8 rRNA genes. J. nudiflorum exhibited a larger genome size compared to other species, mainly attributed to the elevated number of forward repeats (FRs). Despite the typically conservative nature of chloroplasts, variations in the presence or absence of accD have been observed within J. sambac. The calculation of nucleotide diversity (Pi) values for 19 cp genomes indicated that potential mutation hotspots were more likely to be located in LSC regions than in other regions, particularly in genes ycf2, rbcL, atpE, ndhK, and ndhC (Pi > 0.2). Ka/Ks values revealed strong selection pressure on the genes rps2, atpA, rpoA, rpoC1, and rpl33 when comparing J. sambac with the three most closely related species (J. auriculatum, J. multiflorum, and J. dichotomum). Additionally, SNP identification, along with the results of Structure, PCA, and phylogenetic tree analyses, divided the Jasminum cp genomes into six groups. Notably, J. polyanthum showed gene flow signals from both the G5 group (J. nudiflorum) and the G3 group (J. tortuosum and J. fluminense). Phylogenetic tree analysis reflected that most species from the same genus clustered together with robust support in Oleaceae, strongly supporting the monophyletic nature of cp genomes within the genus Jasminum.

CONCLUSION

Overall, this study provides comprehensive insights into the genomic composition, variation, and phylogenetic relationships among various Jasminum species. These findings enhance our understanding of the genetic diversity and evolutionary history of Jasminum.

The newly assembled chloroplast genome of Aeluropus littoralis: molecular feature characterization and phylogenetic analysis with related species

Aeluropus littoralis, a halophyte grass, is widely distributed from the Mediterranean to the Indian subcontinent through the Mongolian Gobi. This model halophyte has garnered increasing attention owing to its use as forage and its high tolerance to environmental stressors. The chloroplast genomes of many plants have been extensively examined for molecular, phylogenetic and transplastomic applications. However, no published research on the A. littoralis chloroplast (cp) genome was discovered. Here, the entire chloroplast genome of A. littoralis was assembled implementing accurate long-read sequences. The entire chloroplast genome, with an estimated length of 135,532 bp (GC content: 38.2%), has a quadripartite architecture and includes a pair of inverted repeat (IR) regions, IRa and IRb (21,012 bp each), separated by a large and a small single-copy regions (80,823 and 12,685 bp, respectively). The features of A. littoralis consist of 133 genes that synthesize 87 peptides, 38 transfer RNAs, and 8 ribosomal RNAs. Of these genes, 86 were unique, whereas 19 were duplicated in IR regions. Additionally, a total of forty-six simple sequence repeats, categorized into 32-mono, four-di, two-tri, and eight-tetranucleotides, were discovered. Furthermore, ten sets of repeats greater than 20 bp were located primarily in the LSC region. Evolutionary analysis based on chloroplast sequence data revealed that A. littoralis with A. lagopoides and A. sinensis belong to the Aeluropodinae subtribe, which is a sister to the Eleusininae in the tribe Cynodonteae and the subfamily Chloridoideae. This subfamily belongs to the PACMAD clade, which contains the majority of the C4 photosynthetic plants in the Poaceae. The newly constructed A. littoralis cp genome offers valuable knowledge for DNA barcoding, phylogenetic, transplastomic research, and other biological studies.

Molecular Characterization and Phylogenetic Analysis of Centipedegrass [Eremochloa ophiuroides (Munro) Hack.] Based on the Complete Chloroplast Genome Sequence

Centipedegrass (Eremochloa ophiuroides) is an important warm-season grass plant used as a turfgrass as well as pasture grass in tropical and subtropical regions, with wide application in land surface greening and soil conservation in South China and southern United States. In this study, the complete cp genome of E. ophiuroides was assembled using high-throughput Illumina sequencing technology. The circle pseudomolecule for E. ophiuroides cp genome is 139,107 bp in length, with a quadripartite structure consisting of a large single copyregion of 82,081 bp and a small single copy region of 12,566 bp separated by a pair of inverted repeat regions of 22,230 bp each. The overall A + T content of the whole genome is 61.60%, showing an asymmetric nucleotide composition. The genome encodes a total of 131 gene species, composed of 20 duplicated genes within the IR regions and 111 unique genes comprising 77 protein-coding genes, 30 transfer RNA genes, and 4 ribosome RNA genes. The complete cp genome sequence contains 51 long repeats and 197 simple sequence repeats, and a high degree of collinearity among E. ophiuroide and other Gramineae plants was disclosed. Phylogenetic analysis showed E. ophiuroides, together with the other two Eremochloa species, is closely related to Mnesithea helferi within the subtribe Rottboelliinae. These findings will be beneficial for the classification and identification of the Eremochloa taxa, phylogenetic resolution, novel gene discovery, and functional genomic studies for the genus Eremochloa.

Extreme Reconfiguration of Plastid Genomes in Papaveraceae: Rearrangements, Gene Loss, Pseudogenization, IR Expansion, and Repeats

The plastid genomes (plastomes) of angiosperms are typically highly conserved, with extreme reconfiguration being uncommon, although reports of such events have emerged in some lineages. In this study, we conducted a comprehensive comparison of the complete plastomes from twenty-two species, covering seventeen genera from three subfamilies (Fumarioideae, Hypecooideae, and Papaveroideae) of Papaveraceae. Our results revealed a high level of variability in the plastid genome size of Papaveraceae, ranging from 151,864 bp to 219,144 bp in length, which might be triggered by the expansion of the IR region and a large number of repeat sequences. Moreover, we detected numerous large-scale rearrangements, primarily occurring in the plastomes of Fumarioideae and Hypecooideae. Frequent gene loss or pseudogenization were also observed for ndhs, accD, clpP, infA, rpl2, rpl20, rpl32, rps16, and several tRNA genes, particularly in Fumarioideae and Hypecooideae, which might be associated with the structural variation in their plastomes. Furthermore, we found that the plastomes of Fumarioideae exhibited a higher GC content and more repeat sequences than those of Papaveroideae. Our results showed that Papaveroideae generally displayed a relatively conserved plastome, with the exception of Eomecon chionantha, while Fumarioideae and Hypecooideae typically harbored highly reconfigurable plastomes, showing high variability in the genome size, gene content, and gene order. This study provides insights into the plastome evolution of Papaveraceae and may contribute to the development of effective molecular markers.

The complete chloroplast genome of Hordeum marinum ssp. marinum

Hordeum marinum ssp. marinum (Huds.) R. J. Soreng (2003) is a halophyte wild relative of barley and wheat, which exhibits remarkable salt tolerance characteristics. In this study, we presented the first characterization of the complete chloroplast genome of H. marinum ssp. marinum. Our findings reveal that the chloroplast genome of H. marinum ssp. marinum consists of a small single-copy region (SSC: 12,715 bp), a large single-copy region (LSC: 81,130 bp), and a pair of inverted repeat regions (IRs: 21,517 bp), amounting to a total length of 136,989 bp. The chloroplast genome encodes 139 genes, including 91 protein-coding genes, 38 tRNA genes, and ten rRNA genes. By utilizing phylogenetic analysis, we determine the evolutionary position of H. marinum in Triticeae. Our study provides valuable insights into the chloroplast genome of H. marinum ssp. marinum, which may have important implications for the improvement of cereal crops.

Analysis of the RNA Editing Sites and Orthologous Gene Function of Transcriptome and Chloroplast Genomes in the Evolution of Five Deutzia Species

In this study, the chloroplast genomes and transcriptomes of five Deutzia genus species were sequenced, characterized, combined, and analyzed. A phylogenetic tree was constructed, including 32 other chloroplast genome sequences of Hydrangeoideae species. The results showed that the five Deutzia chloroplast genomes were typical circular genomes 156,860-157,025 bp in length, with 37.58-37.6% GC content. Repeat analysis showed that the Deutzia species had 41-45 scattered repeats and 199-201 simple sequence repeats. Comparative genomic and pi analyses indicated that the genomes are conservative and that the gene structures are stable. According to the phylogenetic tree, Deutzia species appear to be closely related to Kirengeshoma palmata and Philadelphus. By combining chloroplast genomic and transcriptomic analyses, 29-31 RNA editing events and 163-194 orthologous genes were identified. The ndh, rpo, rps, and atp genes had the most editing sites, and all RNA editing events were of the C-to-U type. Most of the orthologous genes were annotated to the chloroplast, mitochondria, and nucleus, with functions including energy production and conversion, translation, and protein transport. Genes related to the biosynthesis of monoterpenoids and flavonoids were also identified from the transcriptome of Deutzia spp. Our results will contribute to further studies of the genomic information and potential uses of the Deutzia spp.

The complete chloroplast genome sequence of Phyllostachys incarnata Wen, 1982

Phyllostachys incarnata Wen, 1982 is one of the important material and edible bamboo specie of high quality in China. We reported the complete chloroplast(cp) genome of P. incarnata in this study. The cp genome of P. incarnata (GenBank accession number: OL457160) was a typical tetrad structure with a full length of 139,689 bp, comprising a pair of inverted repeated (IR) regions (21,798 bp) separated by a large single-copy (LSC) region (83,221 bp) and a small single-copy (SSC) region (12,872 bp). And the cp genome contained 136 genes, including 90 protein-coding genes, 38 tRNA genes, and 8 rRNA genes. Phylogenetic analysis based on 19 cp genomes suggested that P. incarnata was relatively close to P. glauca among the species analyzed.

Backbone phylogeny and adaptive evolution of Pleurospermum s. l.: New insights from phylogenomic analyses of complete plastome data

Pleurospermum is a taxonomically challenging taxon of Apiaceae, as its circumscription and composition remain controversial for morphological similarities with several related genera, leading to a dispute between Pleurospermum in the broad sense and strict sense. While evidence from previous molecular studies recognized plural branching lineages within the Pleurospermum s. l., it did not support the latest delimitation of Pleurospermum s. str. by only two closely related northern species. So far, no proper delimitation for Pleurospermum has come up, and many of the plural taxa in Pleurospermum s. l. remain unresolved, which may be due to poor phylogenetic resolution yielded barely from ITS sequences. Herein, we newly assembled 40 complete plastomes from 36 species of Pleurospermum s. l. and related genera, 34 of which were first reported and generated a well-resolved backbone phylogeny in a framework of the subfamily Apioideae. From the phylogeny with greatly improved resolution, a total of six well-supported monophyletic lineages within Pleurospermum s. l. were recognized falling in different major clades of Apioideae. Combining morphological characteristics with phylogenetic inference, we suggested to re-delimit the Pleurospermum s. str. by introducing nine species mainly from the Himalayan regions and proposed its boundary features; the remaining species were suggested to be excluded from Pleurospermum to incorporate into their more related taxa being revealed. On this basis, the plastome comparison revealed not only the high conservatism but also the mild differences among lineages in plastome structure and gene evolution. Overall, our study provided a backbone phylogeny essential for further studies of the taxonomically difficult taxa within Pleurospermum s. l.

Complete chloroplast genome sequences of the medicinal plant Aconitum transsectum (Ranunculaceae): comparative analysis and phylogenetic relationships

BACKGROUND

Aconitum transsectum Diels. (Ranunculaceae) is an important medicinal plant that is widely used in traditional Chinese medicine, but its morphological traits make it difficult to recognize from other Aconitum species. No research has sequenced the chloroplast genome of A.transsectum, despite the fact that phylogenetic analysis based on chloroplast genome sequences provides essential evidence for plant classification.

RESULTS

In this study, the chloroplast (cp) genome of A. transsectum was sequenced, assembled, and annotated. A. transsectum cp genome is a 155,872 bp tetrameric structure including a large single copy (LSC, 87,671 bp) and a small single copy (SSC, 18,891 bp) section, as well as a pair of inverted repeat sequences (IRa and IRb, 25,894 bp each). 131 genes are encoded by the complete cp genome, comprising 86 protein-coding genes, 37 tRNAs, and 8 rRNAs. The most favored codon in the A. transsectum cp genome is AUG, and 46 repeats and 241 SSRs were also identified. The A. transsectum cp genome is similar in size, gene composition, and IR expansion and contraction to the cp genomes of seven Ranunculaceae species. Phylogenetic analysis of cp genomes of 28 plants from the Ranunculaceae family shows that A. transsectum is most closely related to A. vilmorinianum, A. episcopale, and A. forrestii of Subgen. Aconitum.

CONCLUSIONS

Overall, this study provides complete cp genome resources for A. transsectum that will be beneficial for identifying potential.

The population genomic analyses of chloroplast genomes shed new insights on the complicated ploidy and evolutionary history in Fragaria

The genus Fragaria consists of a rich diversity of ploidy levels with diploid (2x), tetraploid (4x), pentaploid (5x), hexaploidy (6x), octoploid (8x) and decaploid (10x) species. Only a few studies have explored the origin of diploid and octoploid strawberry, and little is known about the roles of tetraploidy and hexaploidy during the evolution of octoploid strawberry. The chloroplast genome is usually a stable circular genome and is widely used in investigating the evolution and matrilineal identification. Here, we assembled the chloroplast genomes of F. x ananassa cv. 'Benihoppe' (8x) using Illumina and HiFi data seperately. The genome alignment results showed that more InDels were located in the chloroplast genomes based on the PacBio HiFi data than Illumina data. We obtain highly accurate chloroplast genomes assembled through GetOrganelle using Illumina reads. We assembled 200 chloroplast genomes including 198 Fragaria (21 species) and 2 Potentilla samples. Analyses of sequence variation, phylogenetic and PCA analyses showed that Fragaria was divided into five groups. F. iinumae, F. nilgerrensis and all octoploid accessions formed Group A, C and E separately. Species native to western China were clustered into Group B. Group D consisted of F. virdis, F. orientalis, F. moschata, and F. vesca. STRUCTURE and haplotype network confirmed that the diploid F. vesca subsp. bracteata was the last maternal donator of octoploid strawberry. The dN/dS ratio estimated for the protein-coding genes revealed that genes involved in ATP synthase and photosystem function were under positive selection. These findings demonstrate the phylogeny of totally 21 Fragaria species and the origin of octoploid species. F. vesca was the last female donator of octoploid, which confirms the hypothesis that the hexaploid species F. moschata may be an evolutionary intermediate between the diploids and wild octoploid species.

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.

Chloroplast genomes of four Carex species: Long repetitive sequences trigger dramatic changes in chloroplast genome structure

The chloroplast genomes of angiosperms usually have a stable circular quadripartite structure that exhibits high consistency in genome size and gene order. As one of the most diverse genera of angiosperms, Carex is of great value for the study of evolutionary relationships and speciation within its genus, but the study of the structure of its chloroplast genome is limited due to its highly expanded and restructured genome with a large number of repeats. In this study, we provided a more detailed account of the chloroplast genomes of Carex using a hybrid assembly of second- and third-generation sequencing and examined structural variation within this genus. The study revealed that chloroplast genomes of four Carex species are significantly longer than that of most angiosperms and are characterized by high sequence rearrangement rates, low GC content and gene density, and increased repetitive sequences. The location of chloroplast genome structural variation in the species of Carex studied is closely related to the positions of long repeat sequences; this genus provides a typical example of chloroplast structural variation and expansion caused by long repeats. Phylogenetic relationships constructed based on the chloroplast protein-coding genes support the latest taxonomic system of Carex, while revealing that structural variation in the chloroplast genome of Carex may have some phylogenetic significance. Moreover, this study demonstrated a hybrid assembly approach based on long and short reads to analyze complex chloroplast genome assembly and also provided an important reference for the analysis of structural rearrangements of chloroplast genomes in other taxa.

The chloroplast genome of Lappula myosotis V. Wolf, a medicinal species

Lappula myosotis V. Wolf 1776 is an annual or biennial plant with important medicinal value. In the present study, we report the complete chloroplast genome data of L. myosotis, which has a length of 146,668 bp, including a small single-copy (SSC) region of 17,059 bp, a large single-copy (LSC) region of 79,691 bp, and a pair of inverted repeats (IRs) of 24,959 bp. A total of 127 genes encoding tRNA and rRNA were annotated. The total CG content of the chloroplast genome was 37.7%. The maximum-likelihood (ML) phylogenetic tree strongly supported that L. myosotis is closely related to Trigonotis peduncularis. The complete chloroplast genome of L. myosotis provides useful information on the evolution and phylogenetic relationship among Boraginaceae plants.

Plastid phylogenomics and morphological character evolution of Chloridoideae (Poaceae)

Chloridoideae is one of the largest subfamilies of Poaceae, containing many species of great economic and ecological value; however, phylogenetic relationships among the subtribes and genera of Cynodonteae are controversial. In the present study, we combined 111 plastomes representing all five tribes, including 25 newly sequenced plastomes that are mostly from Cynodonteae. Phylogenetic analyses supported the five monophyletic tribes of Chloridoideae, including Centropodieae, Triraphideae, Eragrostideae, Zoysieae and Cynodonteae. Simultaneously, nine monophyletic lineages were revealed in Cynodonteae: supersubtribe Boutelouodinae, subtribes Tripogoninae, Aeluropodinae, Eleusininae, Dactylocteniinae, supersubtribe Gouiniodinae, Cleistogenes and Orinus, and subtribe Triodiinae. Within the tribe of Cynodonteae, the basal lineage is supersubtribe Boutelouodinae and Tripogoninae is sister to the remaining lineages. The clade formed of Aeluropodinae and Eleusininae is sister to the clade composed of Dactylocteniinae, supersubtribe Gouiniodinae, Cleistogenes and Orinus, and subtribe Triodiinae. The clade comprising Dactylocteniinae and supersubtribe Gouiniodinae is sister to the clade comprising Cleistogenes, Orinus, and Triodiinae. Acrachne is a genus within Eleusininae but not within Dactylocteniinae. Molecular evidence determined that Diplachne is not clustered with Leptochloa, which indicated that Diplachne should not be combined into Leptochloa. Cleistogenes is sister to a clade composed of Orinus and Triodia, whereas the recently proposed subtribe Orininae was not supported. Cynodonteae was estimated to have experienced rapid divergence within a short period, which could be a major obstacle in resolving its phylogenetic relationships. Ancestral state reconstructions of morphological characters showed that the most recent common ancestor (MRCA) of Chloridoideae has a panicle, multiple florets in each spikelet, the peaked type of stomatal subsidiary cells, and a saddle-shaped phytoliths, while the ancestral morphological characters of Cynodonteae are the panicle, peaked type of stomatal subsidiary cells, sharp-cap cell typed and equal-base-cell microhair, and square-shaped phytoliths. Overall, plastome phylogenomics provides new insights into the phylogenetic relationships and morphological character evolution of Chloridoideae.

Phylogenomic Analysis of the Plastid Genome of the Peruvian Purple Maize Zea mays subsp. mays cv. 'INIA 601'

Peru is an important center of diversity for maize; its different cultivars have been adapted to distinct altitudes and water availability and possess an array of kernel colors (red, blue, and purple), which are highly appreciated by local populations. Specifically, Peruvian purple maize is a collection of native landraces selected and maintained by indigenous cultures due to its intense purple color in the seed, bract, and cob. This color is produced by anthocyanin pigments, which have gained interest due to their potential use in the food, agriculture, and pharmaceutical industry. It is generally accepted that the Peruvian purple maize originated from a single ancestral landrace 'Kculli', but it is not well understood. To study the origin of the Peruvian purple maize, we assembled the plastid genomes of the new cultivar 'INIA 601' with a high concentration of anthocyanins, comparing them with 27 cultivars/landraces of South America, 9 Z. mays subsp. parviglumis, and 5 partial genomes of Z. mays subsp. mexicana. Using these genomes, plus four other maize genomes and two outgroups from the NCBI database, we reconstructed the phylogenetic relationship of Z. mays. Our results suggest a polyphyletic origin of purple maize in South America and agree with a complex scenario of domestication with recurrent gene flow from wild relatives. Additionally, we identify 18 plastid positions that can be used as high-confidence genetic markers for further studies. Altogether, these plastid genomes constitute a valuable resource to study the evolution and domestication of Z. mays in South America.

C-to-U RNA Editing: A Site Directed RNA Editing Tool for Restoration of Genetic Code

The restoration of genetic code by editing mutated genes is a potential method for the treatment of genetic diseases/disorders. Genetic disorders are caused by the point mutations of thymine (T) to cytidine (C) or guanosine (G) to adenine (A), for which gene editing (editing of mutated genes) is a promising therapeutic technique. In C-to-Uridine (U) RNA editing, it converts the base C-to-U in RNA molecules and leads to nonsynonymous changes when occurring in coding regions; however, for G-to-A mutations, A-to-I editing occurs. Editing of C-to-U is not as physiologically common as that of A-to-I editing. Although hundreds to thousands of coding sites have been found to be C-to-U edited or editable in humans, the biological significance of this phenomenon remains elusive. In this review, we have tried to provide detailed information on physiological and artificial approaches for C-to-U RNA editing.

Complete Chloroplast Genome of Gladiolus gandavensis (Gladiolus) and Genetic Evolutionary Analysis

Gladiolus is an important ornamental plant that is one of the world’s four most-grown cut flowers. Gladiolus gandavensis has only been found in the Cangnan County (Zhejiang Province) of China, which is recorded in the “Botanical”. To explore the origin of G. gandavensis, chloroplast genome sequencing was conducted. The results indicated that a total of 151,654 bp of circular DNA was obtained. The chloroplast genome of G. gandavensis has a quadripartite structure (contains a large single-copy (LSC) region (81,547 bp), a small single-copy region (SSC) (17,895 bp), and two inverted repeats (IRs) (IRa and IRb, 52,212 bp)), similar to that of other species. In addition, a total of 84 protein-coding genes, 8 rRNA-encoding genes, and 38 tRNA-encoding genes were present in the chloroplast genome. To further study the structural characteristics of the chloroplast genome in G. gandavensis, a comparative analysis of eight species of the Iridaceae family was conducted, and the results revealed higher similarity in the IR regions than in the LSC and SSC regions. In addition, 265 simple sequence repeats (SSRs) were detected in this study. The results of the phylogenetic analysis indicated that the chloroplast genome of G. gandavensis has high homology with the Crocus cartwrightianus and Crocus sativus chloroplast genomes. Genetic analysis based on the rbcl sequence among 49 Gladiolus species showed that samples 42, 49, 50, and 54 had high homology with the three samples from China (64, 65, and 66), which might be caused by chance similarity in genotypes. These results suggest that G. gandavensis may have originated from South Africa.

Comparative chloroplast genome analyses of 23 species in Swertia L. (Gentianaceae) with implications for its phylogeny

Swertia L. is a large genus in the family Gentianaceae. Different chloroplast gene segments have been used to study systematic evolutionary relationships between species of Swertia L. However, as gene fragment–based phylogenies lack sufficient resolution, the systematic evolutionary relationships between Swertia L. species have remained unclear. We sequenced and annotated the complete chloroplast genomes of four Swertia species, namely, S. bifolia, S. tetraptera, S. franchetian, and S. przewalskii, using next generation sequencing and the plastid genome annotator tool. The chloroplast genome sequences of 19 additional species of Swertia L. were downloaded from the NCBI database and also assessed. We found that all 23 Swertia L. species had a similar genetic structure, that is, a ring tetrad structure, but with some clear differences. The chloroplast genomes of the 23 Swertia L. species were 149036–153691 bp long, averaging 152385 bp; the genomes contained 134 functional genes: 38 tRNA, eight rRNA, and 88 protein-encoding genes. A comparative analysis showed that chloroplasts genome of Swertia was conserved in terms of genome structure, codon preference, and repeat sequences, but it differed in terms of genome sizes, gene contents, and SC/IR boundary. Using Swertia wolfangiana as a reference, we found clear divergences in most of the non-coding and intergenic regions of the complete chloroplast genomes of these species; we also found that rpoC1, ccsA, ndhI, ndhA, and rps15 protein-coding genes had large variations. These highly variable hotspots will be useful for future phylogenetic and population genetic studies. Phylogenetic analysis with high bootstrap support showed that Swertia L. was not monophyletic. The classification of subgen. Swertia and subgen. Ophelia was supported by molecular data, which also partly supported the division of sect. Ophelia, sect. Platynema, sect. Poephila, sect. Swertia, and sect. Macranthos. However, the systematic positions of other groups and species require further exploration. The Swertia L formed at 29.60 Ma. Speciation of 10 species occurred in succession after 12 Ma and 13 species occurred in succession after 2.5 Ma. Our analysis provides insight into the unresolved evolutionary relationships of Swertia L. species.

Comparative and Phylogenetic Analysis of Complete Chloroplast Genomes in Leymus (Triticodae, Poaceae)

Leymus is a perennial genus that belongs to the tribe Triticeae (Poaceae) which has an adaptive capacity to ecological conditions and strong resistance to cold, drought, and salinity. Most Leymus species are fine herbs that can be used for agriculture, conservation, and landscaping. Due to confusion taxonomy within genera, the complete chloroplast (cp) genome of 13 Leymus species was sequenced, assembled, and compared with those of three other previously published Leymus species (Leymus condensatus, Leymus angustus, and Leymus mollis) to clarify the issue. Overall, the whole cp genome size ranged between 135,057 (L. condensatus) and 136,906 bp (Leymus coreanus) and showed a typical quadripartite structure. All studied species had 129 genes, including 83 protein-coding genes, 38 transfer RNAs, and 8 ribosomal RNAs. In total, 800 tandem repeats and 707 SSR loci were detected, most of which were distributed in the large single-copy region, followed by the inverted repeat (IR) and small single-copy regions. The sequence identity of all sequences was highly similar, especially concerning the protein-coding and IR regions; in particular, the protein-coding regions were significantly similar to those in the IR regions, regardless of small sequence differences in the whole cp genome. Moreover, the coding regions were more conserved than the non-coding regions. Comparisons of the IR boundaries showed that IR contraction and expansion events were reflected in different locations of rpl22, rps19, ndhH, and psbA genes. The close phylogenetic relationship of Leymus and Psathyrostachys indicated that Psathyrostachys possibly is the donor of the Ns genome sequence identified in Leymus. Altogether, the complete cp genome sequence of Leymus will lay a solid foundation for future population genetics and phylogeography studies, as well as for the analysis of the evolution of economically valuable plants.

Plastome comparison and phylogenomics of Fagopyrum (Polygonaceae): insights into sequence differences between Fagopyrum and its related taxa

BACKGROUND

Fagopyrum (Polygonaceae) is a small plant lineage comprised of more than fifteen economically and medicinally important species. However, the phylogenetic relationships of the genus are not well explored, and the characteristics of Fagopyrum chloroplast genomes (plastomes) remain poorly understood so far. It restricts the comprehension of species diversity in Fagopyrum. Therefore, a comparative plastome analysis and comprehensive phylogenomic analyses are required to reveal the taxonomic relationship among species of Fagopyrum.

RESULTS

In the current study, 12 plastomes were sequenced and assembled from eight species and two varieties of Fagopyrum. In the comparative analysis and phylogenetic analysis, eight previously published plastomes of Fagopyrum were also included. A total of 49 plastomes of other genera in Polygonaceae were retrieved from GenBank and used for comparative analysis with Fagopyrum. The variation of the Fagopyrum plastomes is mainly reflected in the size and boundaries of inverted repeat/single copy (IR/SC) regions. Fagopyrum is a relatively basal taxon in the phylogenomic framework of Polygonaceae comprising a relatively smaller plastome size (158,768-159,985 bp) than another genus of Polygonaceae (158,851-170,232 bp). A few genera of Polygonaceae have nested distribution of the IR/SC boundary variations. Although most species of Fagopyrum show the same IRb/SC boundary with species of Polygonaceae, only a few species show different IRa/SC boundaries. The phylogenomic analyses of Fagopyrum supported the cymosum and urophyllum groups and resolved the systematic position of subclades within the urophyllum group. Moreover, the repeat sequence types and numbers were found different between groups of Fagopyrum. The plastome sequence identity showed significant differences between intra-group and inter-group.

CONCLUSIONS

The deletions of intergenic regions cause a short length of Fagopyrum plastomes, which may be the main reason for plastome size diversity in Polygonaceae species. The phylogenomic reconstruction combined with the characteristics comparison of plastomes supports grouping within Fagopyrum. The outcome of these genome resources may facilitate the taxonomy, germplasm resources identification as well as plant breeding of Fagopyrum.

The complete chloroplast genome of critically endangered Chimonobambusa hirtinoda (Poaceae: Chimonobambusa) and phylogenetic analysis

Chimonobambusa hirtinoda, a threatened species, is only naturally distributed in Doupeng Mountain, Duyun, Guizhou, China. Next-generation sequencing (NGS) is used to obtain the complete chloroplast (cp) genome sequence of C. hirtinoda. The sequence was assembled and analyzed for phylogenetic and evolutionary studies. Additionally, we compared the cp genome of C. hirtinoda with previously published Chimonobambusa species. The cp genome of C. hirtinoda has a total length of 139, 561 bp and 38.90% GC content. This genome included a large single -copy (LSC) region of 83, 166 bp, a small single-copy (SSC) region of 20, 811 bp and a pair of inverted repeats of 21,792 bp each. We discovered 130 genes in the cp genome, including 85 protein-coding genes, 37 tRNA, and 8 rRNA genes. A total of 48 simple sequence repeats (SSRs) were detected. The A/U preference of the third nucleotide in the cp genome of C. hirtinoda was obtained by measuring the codon usage frequency of amino acids. Furthermore, phylogenetic analysis using complete cp sequences and matK gene revealed a genetic relationship within the Chimonobambusa genus. This study reported the chloroplast genome of the C. hirtinoda.

Comparative analysis of two Korean irises (Iris ruthenica and I. uniflora, Iridaceae) based on plastome sequencing and micromorphology

Iris ruthenica Ker Gawl. and I. uniflora Pall. ex Link, which are rare and endangered species in Korea, possess considerable horticultural and medicinal value among Korean irises. However, discrimination of the species is hindered by extensive morphological similarity. Thus, the aim of the present study was to identify discriminating features by comparing the species’ complete plastid genome (i.e., plastome) sequences and micromorphological features, including leaf margins, stomatal complex distribution (hypostomatic vs. amphistomatic leaves), anther stomata density, and tepal epidermal cell patterns. Plastome comparison revealed slightly divergent regions within intergenic spacer regions, and the most variable sequences, which were distributed in non-coding regions, could be used as molecular markers for the discrimination of I. ruthenica and I. uniflora. Phylogenetic analysis of the Iris species revealed that I. ruthenica and I. uniflora formed a well-supported clade. The comparison of plastomes and micromorphological features performed in this study provides useful information for elucidating taxonomic, phylogenetic, and evolutionary relationships in Iridaceae. Further studies, including those based on molecular cytogenetic approaches using species specific markers, will offer insights into species delimitation of the two closely related Iris species.

Analysis of the Plastid Genome Sequence During Maize Seedling Development

Shoot development in maize progresses from small, non-pigmented meristematic cells to expanded cells in the green leaf. During this transition, large plastid DNA (ptDNA) molecules in proplastids become fragmented in the photosynthetically-active chloroplasts. The genome sequences were determined for ptDNA obtained from Zea mays B73 plastids isolated from four tissues: base of the stalk (the meristem region); fully-developed first green leaf; first three leaves from light-grown seedlings; and first three leaves from dark-grown (etiolated) seedlings. These genome sequences were then compared to the Z. mays B73 plastid reference genome sequence that was previously obtained from green leaves. The assembled plastid genome was identical among these four tissues to the reference genome. Furthermore, there was no difference among these tissues in the sequence at and around the previously documented 27 RNA editing sites. There were, however, more sequence variants (insertions/deletions and single-nucleotide polymorphisms) for leaves grown in the dark than in the light. These variants were tightly clustered into two areas within the inverted repeat regions of the plastid genome. We propose a model for how these variant clusters could be generated by replication-transcription conflict.

Comprehensive analysis of plastid gene expression during fruit development and ripening of kiwifruit

Global survey of plastid gene expression during fruit ripening in kiwifruit provides cis-elements for the future engineering of the plastid genome of kiwifruit. A limitation in the application of plastid biotechnology for molecular farming is the low-level expression of transgenes in non-green plastids compared with photosynthetically active chloroplasts. Unlike other fruits, not all chloroplasts are transformed into chromoplasts during ripening of red-fleshed kiwifruit (Actinidia chinensis cv. Hongyang) fruits, which may make kiwifruit an ideal horticultural plant for recombinant protein production by plastid engineering. To identify cis-elements potentially triggering high-level transgene expression in edible tissues of the 'Hongyang' kiwifruit, here we report a comprehensive analysis of kiwifruit plastid gene transcription in green leaves and fruits at three different developmental stages. While transcripts of a few photosynthesis-related genes and most genetic system genes were substantially upregulated in green fruits compared with leaves, nearly all plastid genes were significantly downregulated at the RNA level during fruit development. Expression of a few genes remained unchanged, including psbA, the gene encoding the D1 polypeptide of photosystem II. However, PsbA protein accumulation decreased continuously during chloroplast-to-chromoplast differentiation. Analysis of post-transcriptional steps in mRNA maturation, including intron splicing and RNA editing, revealed that splicing and editing may contribute to regulation of plastid gene expression. Altogether, 40 RNA editing sites were verified, and 5 of them were newly discovered. Taken together, this study has generated a valuable resource for the analysis of plastid gene expression and provides cis-elements for future efforts to engineer the plastid genome of kiwifruit.

The complete chloroplast genome sequencing analysis revealed an unusual IRs reduction in three species of subfamily Zygophylloideae

Tetraena mongolica, Zygophyllum xanthoxylon, and Z. fabago are three typical dryland plants with important ecological values in subfamily Zygophylloideae of Zygophyllaceae. Studies on the chloroplast genomes of them are favorable for understanding the diversity and phylogeny of Zygophyllaceae. Here, we sequenced and assembled the whole chloroplast genomes of T. mongolica, Z. xanthoxylon, and Z. fabago, and performed comparative genomic and phylogenetic analysis. The total size, structure, gene content and orders of these three chloroplast genomes were similar, and the three chloroplast genomes exhibited a typical quadripartite structure with a large single-copy region (LSC; 79,696–80,291 bp), a small single-copy region (SSC; 16,462–17,162 bp), and two inverted repeats (IRs; 4,288–4,413 bp). A total of 107 unique genes were identified from the three chloroplast genomes, including 70 protein-coding genes, 33 tRNAs, and 4 rRNAs. Compared with other angiosperms, the three chloroplast genomes were significantly reduced in overall length due to an unusual 16–24 kb shrinkage of IR regions and loss of the 11 genes which encoded subunits of NADH dehydrogenase. Genome-wide comparisons revealed similarities and variations between the three species and others. Phylogenetic analysis based on the three chloroplast genomes supported the opinion that Zygophyllaceae belonged to Zygophyllales in Fabids, and Z. xanthoxylon and Z. fabago belonged to Zygophyllum. The genome-wide comparisons revealed the similarity and variations between the chloroplast genomes of the three Zygophylloideae species and other plant species. This study provides a valuable molecular biology evidence for further studies of phylogenetic status of Zygophyllaceae.

Beyond RuBisCO: convergent molecular evolution of multiple chloroplast genes in C4 plants

BACKGROUND

The recurrent evolution of the C4 photosynthetic pathway in angiosperms represents one of the most extraordinary examples of convergent evolution of a complex trait. Comparative genomic analyses have unveiled some of the molecular changes associated with the C4 pathway. For instance, several key enzymes involved in the transition from C3 to C4 photosynthesis have been found to share convergent amino acid replacements along C4 lineages. However, the extent of convergent replacements potentially associated with the emergence of C4 plants remains to be fully assessed. Here, we conducted an organelle-wide analysis to determine if convergent evolution occurred in multiple chloroplast proteins beside the well-known case of the large RuBisCO subunit encoded by the chloroplast gene rbcL.

METHODS

Our study was based on the comparative analysis of 43 C4 and 21 C3 grass species belonging to the PACMAD clade, a focal taxonomic group in many investigations of C4 evolution. We first used protein sequences of 67 orthologous chloroplast genes to build an accurate phylogeny of these species. Then, we inferred amino acid replacements along 13 C4 lineages and 9 C3 lineages using reconstructed protein sequences of their reference branches, corresponding to the branches containing the most recent common ancestors of C4-only clades and C3-only clades. Pairwise comparisons between reference branches allowed us to identify both convergent and non-convergent amino acid replacements between C4:C4, C3:C3 and C3:C4 lineages.

RESULTS

The reconstructed phylogenetic tree of 64 PACMAD grasses was characterized by strong supports in all nodes used for analyses of convergence. We identified 217 convergent replacements and 201 non-convergent replacements in 45/67 chloroplast proteins in both C4 and C3 reference branches. C4:C4 branches showed higher levels of convergent replacements than C3:C3 and C3:C4 branches. Furthermore, we found that more proteins shared unique convergent replacements in C4 lineages, with both RbcL and RpoC1 (the RNA polymerase beta' subunit 1) showing a significantly higher convergent/non-convergent replacements ratio in C4 branches. Notably, more C4:C4 reference branches showed higher numbers of convergent vs. non-convergent replacements than C3:C3 and C3:C4 branches. Our results suggest that, in the PACMAD clade, C4 grasses experienced higher levels of molecular convergence than C3 species across multiple chloroplast genes. These findings have important implications for our understanding of the evolution of the C4 photosynthesis pathway.

Extensive reorganization of the chloroplast genome of Corydalis platycarpa: A comparative analysis of their organization and evolution with other Corydalis plastomes

Introduction

The chloroplast (cp) is an autonomous plant organelle with an individual genome that encodes essential cellular functions. The genome architecture and gene content of the cp is highly conserved in angiosperms. The plastome of Corydalis belongs to the Papaveraceae family, and the genome is comprised of unusual rearrangements and gene content. Thus far, no extensive comparative studies have been carried out to understand the evolution of Corydalis chloroplast genomes.

Methods

Therefore, the Corydalis platycarpa cp genome was sequenced, and wide-scale comparative studies were conducted using publicly available twenty Corydalis plastomes.

Results

Comparative analyses showed that an extensive genome rearrangement and IR expansion occurred, and these events evolved independently in the Corydalis species. By contrast, the plastomes of its closely related subfamily Papaveroideae and other Ranunculales taxa are highly conserved. On the other hand, the synapomorphy characteristics of both accD and the ndh gene loss events happened in the common ancestor of the Corydalis and sub-clade of the Corydalis lineage, respectively. The Corydalis-sub clade species (ndh lost) are distributed predominantly in the Qinghai-Tibetan plateau (QTP) region. The phylogenetic analysis and divergence time estimation were also employed for the Corydalis species.

Discussion

The divergence time of the ndh gene in the Corydalis sub-clade species (44.31 - 15.71 mya) coincides very well with the uplift of the Qinghai-Tibet Plateau in Oligocene and Miocene periods, and maybe during this period, it has probably triggered the radiation of the Corydalis species.

Conclusion

To the best of the authors' knowledge, this is the first large-scale comparative study of Corydalis plastomes and their evolution. The present study may provide insights into the plastome architecture and the molecular evolution of Corydalis species.

On the Edge of Dispensability, the Chloroplast ndh Genes

The polypeptides encoded by the chloroplast ndh genes and some nuclear genes form the thylakoid NADH dehydrogenase (Ndh) complex, homologous to the mitochondrial complex I. Except for Charophyceae (algae related to higher plants) and a few Prasinophyceae, all eukaryotic algae lack ndh genes. Among vascular plants, the ndh genes are absent in epiphytic and in some species scattered among different genera, families, and orders. The recent identification of many plants lacking plastid ndh genes allows comparison on phylogenetic trees and functional investigations of the ndh genes. The ndh genes protect Angiosperms under various terrestrial stresses, maintaining efficient photosynthesis. On the edge of dispensability, ndh genes provide a test for the natural selection of photosynthesis-related genes in evolution. Variable evolutionary environments place Angiosperms without ndh genes at risk of extinction and, probably, most extant ones may have lost ndh genes recently. Therefore, they are evolutionary endpoints in phylogenetic trees. The low number of sequenced plastid DNA and the long lifespan of some Gymnosperms lacking ndh genes challenge models about the role of ndh genes protecting against stress and promoting leaf senescence. Additional DNA sequencing in Gymnosperms and investigations into the molecular mechanisms of their response to stress will provide a unified model of the evolutionary and functional consequences of the lack of ndh genes.

Sample Sequence Analysis Uncovers Recurrent Horizontal Transfers of Transposable Elements among Grasses

Limited genome resources are a bottleneck in the study of horizontal transfer (HT) of DNA in plants. To solve this issue, we tested the usefulness of low-depth sequencing data generated from 19 previously uncharacterized panicoid grasses for HT investigation. We initially searched for horizontally transferred LTR-retrotransposons by comparing the 19 sample sequences to 115 angiosperm genome sequences. Frequent HTs of LTR-retrotransposons were identified solely between panicoids and rice (Oryza sativa). We consequently focused on additional Oryza species and conducted a nontargeted investigation of HT involving the panicoid genus Echinochloa, which showed the most HTs in the first set of analyses. The comparison of nine Echinochloa samples and ten Oryza species identified recurrent HTs of diverse transposable element (TE) types at different points in Oryza history, but no confirmed cases of HT for sequences other than TEs. One case of HT was observed from one Echinochloa species into one Oryza species with overlapping geographic distributions. Variation among species and data sets highlights difficulties in identifying all HT, but our investigations showed that sample sequence analyses can reveal the importance of HT for the diversification of the TE repertoire of plants.

Complete chloroplast genome of the medicinal plant Evolvulus alsinoides: comparative analysis, identification of mutational hotspots and evolutionary dynamics with species of Solanales

Evolvulus alsinoides, belonging to the family Convolvulaceae, is an important medicinal plant widely used as a nootropic in the Indian traditional medicine system. In the genus Evolvulus, no research on the chloroplast genome has been published. Hence, the present study focuses on annotation, characterization, identification of mutational hotspots, and phylogenetic analysis in the complete chloroplast genome (cp) of E. alsinoides. Genome comparison and evolutionary dynamics were performed with the species of Solanales. The cp genome has 114 genes (80 protein-coding genes, 30 transfer RNA, and 4 ribosomal RNA genes) that were unique with total genome size of 157,015 bp. The cp genome possesses 69 RNA editing sites and 44 simple sequence repeats (SSRs). Predicted SSRs were randomly selected and validated experimentally. Six divergent hotspots such as trnQ-UUG, trnF-GAA, psaI, clpP, ndhF, and ycf1 were discovered from the cp genome. These microsatellites and divergent hot spot sequences of the Taxa 'Evolvulus' could be employed as molecular markers for species identification and genetic divergence investigations. The LSC area was found to be more conserved than the SSC and IR region in genome comparison. The IR contraction and expansion studies show that nine genes rpl2, rpl23, ycf1, ycf2, ycf1, ndhF, ndhA, matK, and psbK were present in the IR-LSC and IR-SSC boundaries of the cp genome. Fifty-four protein-coding genes in the cp genome were under negative selection pressure, indicating that they were well conserved and were undergoing purifying selection. The phylogenetic analysis reveals that E. alsinoides is closely related to the genus Cressa with some divergence from the genus Ipomoea. This is the first time the chloroplast genome of the genus Evolvulus has been published. The findings of the present study and chloroplast genome data could be a valuable resource for future studies in population genetics, genetic diversity, and evolutionary relationship of the family Convolvulaceae.

SUPPLEMENTARY INFORMATION

The online version contains supplementary material available at 10.1007/s12298-021-01051-w.

The chicken or the egg? Plastome evolution and an independent loss of the inverted repeat in papilionoid legumes

The plastid genome (plastome), while surprisingly constant in gene order and content across most photosynthetic angiosperms, exhibits variability in several unrelated lineages. During the diversification history of the legume family Fabaceae, plastomes have undergone many rearrangements, including inversions, expansion, contraction and loss of the typical inverted repeat (IR), gene loss and repeat accumulation in both shared and independent events. While legume plastomes have been the subject of study for some time, most work has focused on agricultural species in the IR-lacking clade (IRLC) and the plant model Medicago truncatula. The subfamily Papilionoideae, which contains virtually all of the agricultural legume species, also comprises most of the plastome variation detected thus far in the family. In this study three non-papilioniods were included among 34 newly sequenced legume plastomes, along with 33 publicly available sequences, to assess plastome structural evolution in the subfamily. In an effort to examine plastome variation across the subfamily, approximately 20% of the sampling represents the IRLC with the remainder selected to represent the early-branching papilionoid clades. A number of IR-related and repeat-mediated changes were identified and examined in a phylogenetic context. Recombination between direct repeats associated with ycf2 resulted in intraindividual plastome heteroplasmy. Although loss of the IR has not been reported in legumes outside of the IRLC, one genistoid taxon was found to completely lack the typical plastome IR. The role of the IR and non-IR repeats in the progression of plastome change is discussed.

Plastome genomics in South American maize landraces: chloroplast lineages parallel the geographical structuring of nuclear gene pools

BACKGROUND AND AIMS

The number of plastome sequences has increased exponentially during the last decade. However, there is still little knowledge of the levels and distribution of intraspecific variation. The aims of this study were to estimate plastome diversity within Zea mays and analyse the distribution of haplotypes in connection with the landrace groups previously delimited for South American maize based on nuclear markers.

METHODS

We obtained the complete plastomes of 30 South American maize landraces and three teosintes by means of next-generation sequencing (NGS) and used them in combination with data from public repositories. After quality filtering, the curated data were employed to search for single-nucleotide polymorphisms, indels and chloroplast simple sequence repeats. Exact permutational contingency tests were performed to assess associations between plastome and nuclear variation. Network and Bayesian phylogenetic analyses were used to infer evolutionary relationships among haplotypes.

KEY RESULTS

Our analyses identified a total of 124 polymorphic plastome loci, with the intergenic regions psbE-rps18, petN-rpoB, trnL_UAG-ndhF and rpoC2-atpI exhibiting the highest marker densities. Although restricted in number, these markers allowed the discrimination of 27 haplotypes in a total of 51 Zea mays individuals. Andean and lowland South American landraces differed significantly in haplotype distribution. However, overall differentiation patterns were not informative with respect to subspecies diversification, as evidenced by the scattered distribution of maize and teosinte plastomes in both the network and Bayesian phylogenetic reconstructions.

CONCLUSIONS

Knowledge of intraspecific plastome variation provides the framework for a more comprehensive understanding of evolutionary processes at low taxonomic levels and may become increasingly important for future plant barcoding efforts. Whole-plastome sequencing provided useful variability to contribute to maize phylogeographic studies. The structuring of haplotype diversity in the maize landraces examined here clearly reflects the distinction between the Andean and South American lowland gene pools previously inferred based on nuclear markers.

Plastid Phylogenomic Data Offers Novel Insights Into the Taxonomic Status of the Trichosanthes kirilowii Complex (Cucurbitaceae) in South Korea

Trichosanthes is a genus in Cucurbitaceae comprising 90-100 species. Trichosanthes species are valuable as herbaceous medicinal ingredients. The fruits, seeds, and roots of species such as T. kirilowii and T. rosthornii are used in Korean traditional herbal medicines. T. rosthornii is only found in China, whereas in South Korea two varieties, T. kirilowii var. kirilowii and T. kirilowii var. japonica, are distributed. T. kirilowii var. kirilowii and T. kirilowii var. japonica have different fruit and leaf shapes but are recognized as belonging to the same species. Furthermore, although its members have herbal medicine applications, genomic information of the genus is still limited. The broad goals of this study were (i) to evaluate the taxonomy of Trichosanthes using plastid phylogenomic data and (ii) provide molecular markers specific for T. kirilowii var. kirilowii and T. kirilowii var. japonica, as these have differences in their pharmacological effectiveness and thus should not be confused and adulterated. Comparison of five Trichosanthes plastid genomes revealed locally divergent regions, mainly within intergenic spacer regions (trnT-UGU-trnL-UAA: marker name Tri, rrn4.5-rrn5: TRr, trnE-UUC-trnT-GGU: TRtt). Using these three markers as DNA-barcodes for important herbal medicine species in Trichosanthes, the identity of Trichosanthes material in commercial medicinal products in South Korea could be successfully determined. Phylogenetic analysis of the five Trichosanthes species revealed that the species are clustered within tribe Sicyoeae. T. kirilowii var. kirilowii and T. rosthornii formed a clade with T. kirilowii var. japonica as their sister group. As T. kirilowii in its current circumscription is paraphyletic and as the two varieties can be readily distinguished morphologically (e.g., in leaf shape), T. kirilowii var. japonica should be treated (again) as an independent species, T. japonica.

The Chloroplast Phylogenomics and Systematics of Zoysia (Poaceae)

The genus Zoysia Willd. (Chloridoideae) is widely distributed from the temperate regions of Northeast Asia—including China, Japan, and Korea—to the tropical regions of Southeast Asia. Among these, four species—Zoysia japonica Steud., Zoysia sinica Hance, Zoysia tenuifolia Thiele, and Zoysia macrostachya Franch. & Sav.—are naturally distributed in the Korean Peninsula. In this study, we report the complete plastome sequences of these Korean Zoysia species (NCBI acc. nos. MF953592, MF967579~MF967581). The length of Zoysia plastomes ranges from 135,854 to 135,904 bp, and the plastomes have a typical quadripartite structure, which consists of a pair of inverted repeat regions (20,962~20,966 bp) separated by a large (81,348~81,392 bp) and a small (12,582~12,586 bp) single-copy region. In terms of gene order and structure, Zoysia plastomes are similar to the typical plastomes of Poaceae. The plastomes encode 110 genes, of which 76 are protein-coding genes, 30 are tRNA genes, and four are rRNA genes. Fourteen genes contain single introns and one gene has two introns. Three evolutionary hotspot spacer regions—atpB~rbcL, rps16~rps3, and rpl32~trnL-UAG—were recognized among six analyzed Zoysia species. The high divergences in the atpB~rbcL spacer and rpl16~rpl3 region are primarily due to the differences in base substitutions and indels. In contrast, the high divergence between rpl32~trnL-UAG spacers is due to a small inversion with a pair of 22 bp stem and an 11 bp loop. Simple sequence repeats (SSRs) were identified in 59 different locations in Z. japonica, 63 in Z. sinica, 62 in Z. macrostachya, and 63 in Z. tenuifolia plastomes. Phylogenetic analysis showed that the Zoysia (Zoysiinae) forms a monophyletic group, which is sister to Sporobolus (Sporobolinae), with 100% bootstrap support. Within the Zoysia clade, the relationship of (Z. sinica, Z japonica), (Z. tenuifolia, Z. matrella), (Z. macrostachya, Z. macrantha) was suggested.

Unprecedented Intraindividual Structural Heteroplasmy in Eleocharis (Cyperaceae, Poales) Plastomes

Plastid genomes (plastomes) of land plants have a conserved quadripartite structure in a gene-dense unit genome consisting of a large inverted repeat that separates two single copy regions. Recently, alternative plastome structures were suggested in Geraniaceae and in some conifers and Medicago the coexistence of inversion isomers has been noted. In this study, plastome sequences of two Cyperaceae, Eleocharis dulcis (water chestnut) and Eleocharis cellulosa (gulf coast spikerush), were completed. Unlike the conserved plastomes in basal groups of Poales, these Eleocharis plastomes have remarkably divergent features, including large plastome sizes, high rates of sequence rearrangements, low GC content and gene density, gene duplications and losses, and increased repetitive DNA sequences. A novel finding among these features was the unprecedented level of heteroplasmy with the presence of multiple plastome structural types within a single individual. Illumina paired-end assemblies combined with PacBio single-molecule real-time sequencing, long-range polymerase chain reaction, and Sanger sequencing data identified at least four different plastome structural types in both Eleocharis species. PacBio long read data suggested that one of the four E. dulcis plastome types predominates.

Phylogenetic position of Bupleurum sikangense inferred from the complete chloroplast genome sequence

Bupleurum sikangense is an endemic species to China distributed in Xizang (Tibet), which has high saikosaponin content and potential medicinal value. Morphologically, it extremely resembles B. commelynoideum. In order to get a better understanding of the relationship between B. sikangense and B. commelynoideum, and on the phylogenetic status of the two species in the genus, the complete chloroplast (cp) genomes of them were sequenced. The genome organization, repeat sequences, codon usage, RNA-editing sites, and variation of their cp genomes revealed high similarity between the species. Some highly variable regions like trnK-UUU_rps16, rps16_trnQ-UUG, ndhC_trnV-UAC, petA_psbJ, accD_psaI, and petL_psbE were identified, providing potential molecular markers for differentiating the two species. Phylogenetic analysis indicated that B. commelynoideum has a closer relationship to B. chinese than that to B. sikangense. Overall, this study will not only improve our knowledge about cp genomes of these two species, and but also provide data for further research on species identification, safe medical application, conservation genetics, etc., of Bupleurum plants.

ZmPPR26, a DYW-type pentatricopeptide repeat protein, is required for C-to-U RNA editing at atpA-1148 in maize chloroplasts

Pentatricopeptide repeat (PPR) proteins are involved in the C-to-U RNA editing of organellar transcripts. The maize genome contains over 600 PPR proteins and few have been found to function in the C-to-U RNA editing in chloroplasts. Here, we report the function of ZmPPR26 in the C-to-U RNA editing and chloroplast biogenesis in maize. ZmPPR26 encodes a DYW-type PPR protein targeted to chloroplasts. The zmppr26 mutant exhibits albino seedling-lethal phenotype. Loss of function of ZmPPR26 abolishes the editing at atpA-1148 site, and decreases the editing at ndhF-62, rpl20-308, rpl2-2, rpoC2-2774, petB-668, rps8-182, and ndhA-50 sites. Overexpression of ZmPPR26 in zmppr26 restores the editing efficiency and rescues the albino seedling-lethal phenotype. Abolished editing at atpA-1148 causes a Leu to Ser change at AtpA-383 that leads to a reduction in the abundance of chloroplast ATP synthase in zmppr26. The accumulation of photosynthetic complexes are also markedly reduced in zmppr26, providing an explanation for the albino seedling-lethal phenotype. These results indicate that ZmPPR26 is required for the editing at atpA-1148 and is important for editing at the other seven sites in maize chloroplasts. The editing at atpA-1148 is critical for AtpA function, assembly of ATP synthase complex, and chloroplast biogenesis in maize.

Validation of a Triplex Quantitative Polymerase Chain Reaction Assay for Detection and Quantification of Traditional Protein Sources, Pisum sativum L. and Glycine max (L.) Merr., in Protein Powder Mixtures

Several botanicals have been traditionally used as protein sources, including the leguminous Pisum sativum L. and Glycine max (L.) Merr. While a rich history exists of cultivating these plants for their whole, protein-rich grain, modern use as powdered supplements present a new challenge in material authentication. The absence of clear morphological identifiers of an intact plant and the existence of long, complex supply chains behoove industry to create quick, reliable analytical tools to identify the botanical source of a protein product (many of which contain multiple sources). The utility of molecular tools for plant-based protein powder authentication is gaining traction, but few validated tools exist. Multiplex quantitative polymerase chain reaction (qPCR) can provide an economical means by which sources can be identified and relative proportions quantified. We followed established guidelines for the design, optimization, and validation of qPCR assay, and developed a triplex qPCR assay that can amplify and quantify pea and soy DNA targets, normalized by a calibrator. The assay was evaluated for analytical specificity, analytical sensitivity, efficiency, precision, dynamic range, repeatability, and reproducibility. We tested the quantitative ability of the assay using pea and soy DNA mixtures, finding exceptional quantitative linearity for both targets - 0.9983 (p < 0.0001) for soy and 0.9915 (p < 0.0001) for pea. Ratios based on mass of protein powder were also tested, resulting in non-linear patterns in data that suggested the requirement of further sample preparation optimization or algorithmic correction. Variation in fragment size within different lots of commercial protein powder samples was also analyzed, revealing low SD among lots. Ultimately, this study demonstrated the utility of qPCR in the context of protein powder mixtures and highlighted key considerations to take into account for commercial implementation.

Complete chloroplast genomes of Achnatherum inebrians and comparative analyses with related species from Poaceae

This article reports the complete chloroplast genome of Achnatherum inebrians, a poisonous herb that is widely distributed in the rangelands of Northern China. The genome is 137 714 bp in total and consists of a large single‐copy (81 758 bp) region and small single‐copy (12 682 bp) region separated by a pair of inverted repeats (21 637 bp). The genome contains 130 genes, including 84 protein‐coding genes, 38 tRNA genes and 8 ribosomal RNA genes, and the guanine + cytosine content is 36.17%. We subsequently performed comparative analysis of complete genomes from A. inebrians and other Poaceae‐related species from GenBank. Thirty‐eight simple sequence repeats were identified, further demonstrating rapid evolution in Poaceae. Finally, the phylogenetic trees of 37 species of Poaceae and 2 species of Amaranthaceae were constructed by using maximum likelihood and Bayesian inference methods, based on the genes of the complete chloroplast genome. We identified hotspots that can be used as molecular markers and barcodes for phylogenetic analysis, as well as for species identification. Phylogenetic analysis indicated that A. inebrians is a member of the genus Stipa rather than Achnatherum.

Green giant-a tiny chloroplast genome with mighty power to produce high-value proteins: history and phylogeny

Free-living cyanobacteria were entrapped by eukaryotic cells ~2 billion years ago, ultimately giving rise to chloroplasts. After a century of debate, the presence of chloroplast DNA was demonstrated in the 1960s. The first chloroplast genomes were sequenced in the 1980s, followed by ~100 vegetable, fruit, cereal, beverage, oil and starch/sugar crop chloroplast genomes in the past three decades. Foreign genes were expressed in isolated chloroplasts or intact plant cells in the late 1980s and stably integrated into chloroplast genomes, with typically maternal inheritance shown in the 1990s. Since then, chloroplast genomes conferred the highest reported levels of tolerance or resistance to biotic or abiotic stress. Although launching products with agronomic traits in important crops using this concept has been elusive, commercial products developed include enzymes used in everyday life from processing fruit juice, to enhancing water absorption of cotton fibre or removal of stains as laundry detergents and in dye removal in the textile industry. Plastid genome sequences have revealed the framework of green plant phylogeny as well as the intricate history of plastid genome transfer events to other eukaryotes. Discordant historical signals among plastid genes suggest possible variable constraints across the plastome and further understanding and mitigation of these constraints may yield new opportunities for bioengineering. In this review, we trace the evolutionary history of chloroplasts, status of autonomy and recent advances in products developed for everyday use or those advanced to the clinic, including treatment of COVID-19 patients and SARS-CoV-2 vaccine.

Differentiation of Hedyotis diffusa and Common Adulterants Based on Chloroplast Genome Sequencing and DNA Barcoding Markers

Chinese herbal tea, also known as Liang Cha or cooling beverage, is popular in South China. It is regarded as a quick-fix remedy to relieve minor health problems. Hedyotis diffusa Willd. (colloquially Baihuasheshecao) is a common ingredient of cooling beverages. H. diffusa is also used to treat cancer and bacterial infections. Owing to the high demand for H. diffusa, two common adulterants, Hedyotis brachypoda (DC.) Sivar and Biju (colloquially Nidingjingcao) and Hedyotis corymbosa (L.) Lam. (colloquially Shuixiancao), are commonly encountered in the market. Owing to the close similarity of their morphological characteristics, it is difficult to differentiate them. Here, we sequenced the complete chloroplast genomes of the three species of Hedyotis using next-generation sequencing (NGS). By comparing the complete chloroplast genomes, we found that they are closely related in the subfamily Rubioideae. We also discovered that there are significant differences in the number and repeating motifs of microsatellites and complex repeats and revealed three divergent hotspots, rps16-trnQ intergenic spacer, ndhD and ycf1. By using these species-specific sequences, we propose new DNA barcoding markers for the authentication of H. diffusa and its two common adulterants.

Plastid Genomes of Flowering Plants: Essential Principles

The plastid genome (plastome ) has proved a valuable source of data for evaluating evolutionary relationships among angiosperms. Through basic and applied approaches, plastid transformation technology offers the potential to understand and improve plant productivity, providing food, fiber, energy, and medicines to meet the needs of a burgeoning global population. The growing genomic resources available to both phylogenetic and biotechnological investigations is allowing novel insights and expanding the scope of plastome research to encompass new species. In this chapter, we present an overview of some of the seminal and contemporary research that has contributed to our current understanding of plastome evolution and attempt to highlight the relationship between evolutionary mechanisms and the tools of plastid genetic engineering.

Variation in Plastome Sizes Accompanied by Evolutionary History in Monogenomic Triticeae (Poaceae: Triticeae)

To investigate the pattern of chloroplast genome variation in Triticeae, we comprehensively analyzed the indels in protein-coding genes and intergenic sequence, gene loss/pseudonization, intron variation, expansion/contraction in inverted repeat regions, and the relationship between sequence characteristics and chloroplast genome size in 34 monogenomic Triticeae plants. Ancestral genome reconstruction suggests that major length variations occurred in four-stem branches of monogenomic Triticeae followed by independent changes in each genus. It was shown that the chloroplast genome sizes of monogenomic Triticeae were highly variable. The chloroplast genome of Pseudoroegneria, Dasypyrum, Lophopyrum, Thinopyrum, Eremopyrum, Agropyron, Australopyrum, and Henradia in Triticeae had evolved toward size reduction largely because of pseudogenes elimination events and length deletion fragments in intergenic. The Aegilops/Triticum complex, Taeniatherum, Secale, Crithopsis, Herteranthelium, and Hordeum in Triticeae had a larger chloroplast genome size. The large size variation in major lineages and their subclades are most likely consequences of adaptive processes since these variations were significantly correlated with divergence time and historical climatic changes. We also found that several intergenic regions, such as petN-trnC and psbE-petL containing unique genetic information, which can be used as important tools to identify the maternal relationship among Triticeae species. Our results contribute to the novel knowledge of plastid genome evolution in Triticeae.

Dissection for Floral Micromorphology and Plastid Genome of Valuable Medicinal Borages Arnebia and Lithospermum (Boraginaceae)

The genera Arnebia and Lithospermum (Lithospermeae-Boraginaceae) comprise 25-30 and 50-60 species, respectively. Some of them are economically valuable, as their roots frequently contain a purple-red dye used in the cosmetic industry. Furthermore, dried roots of Arnebia euchroma, A. guttata, and Lithospermum erythrorhizon, which have been designated Lithospermi Radix, are used as traditional Korean herbal medicine. This study is the first report on the floral micromorphology and complete chloroplast (cp) genome sequences of A. guttata (including A. tibetana), A. euchroma, and L. erythrorhizon. We reveal great diversity in floral epidermal cell patterns, gynoecium, and structure of trichomes. The cp genomes were 149,361-150,465 bp in length, with conserved quadripartite structures. In total, 112 genes were identified, including 78 protein-coding regions, 30 tRNA genes, and four rRNA genes. Gene order, content, and orientation were highly conserved and were consistent with the general structure of angiosperm cp genomes. Comparison of the four cp genomes revealed locally divergent regions, mainly within intergenic spacer regions (atpH-atpI, petN-psbM, rbcL-psaI, ycf4-cemA, ndhF-rpl32, and ndhC-trnV-UAC). To facilitate species identification, we developed molecular markers psaA- ycf3 (PSY), trnI-CAU- ycf2 (TCY), and ndhC-trnV-UAC (NCTV) based on divergence hotspots. High-resolution phylogenetic analysis revealed clear clustering and a close relationship of Arnebia to its Lithospermum sister group, which was supported by strong bootstrap values and posterior probabilities. Overall, gynoecium characteristics and genetic distance of cp genomes suggest that A. tibetana, might be recognized as an independent species rather than a synonym of A. guttata. The present morphological and cp genomic results provide useful information for future studies, such as taxonomic, phylogenetic, and evolutionary analysis of Boraginaceae.

Setaria viridis chlorotic and seedling-lethal mutants define critical functions for chloroplast gene expression

Deep insights into chloroplast biogenesis have been obtained by mutant analysis; however, in C₄ plants a relevant mutant collection has only been developed and exploited for maize. Here, we report the initial characterization of an ethyl methyl sulfonate-induced mutant population for the C₄ model Setaria viridis. Approximately 1000 M₂ families were screened for the segregation of pale-green seedlings in the M₃ generation, and a subset of these was identified to be deficient in post-transcriptional steps of chloroplast gene expression. Causative mutations were identified for three lines using deep sequencing-based bulked segregant analysis, and in one case confirmed by transgenic complementation. Using chloroplast RNA-sequencing and other molecular assays, we describe phenotypes of mutants deficient in PSRP7, a plastid-specific ribosomal protein, OTP86, an RNA editing factor, and cpPNP, the chloroplast isozyme of polynucleotide phosphorylase. The psrp mutant is globally defective in chloroplast translation, and has varying deficiencies in the accumulation of chloroplast-encoded proteins. The otp86 mutant, like its Arabidopsis counterpart, is specifically defective in editing of the rps14 mRNA; however, the conditional pale-green mutant phenotype contrasts with the normal growth of the Arabidopsis mutant. The pnp mutant exhibited multiple defects in 3' end maturation as well as other qualitative changes in the chloroplast RNA population. Overall, our collection opens the door to global analysis of photosynthesis and early seedling development in an emerging C₄ model.

Chloroplast Genome Analysis of Two Medicinal Coelogyne spp. (Orchidaceae) Shed Light on the Genetic Information, Comparative Genomics, and Species Identification

Although the medicinal properties of Coelogyne spp. have been previously studied, there is little genomic information providing a valuable tool for the plant taxonomy, conservation, and utilization of this genus. This study used the next-generation MiSeq sequencing platform to characterize the chloroplast (cp) genomes of Coelogyne fimbriata and Coelogyne ovalis. The Maximum Likelihood (ML) and Bayesian (BI) methods were employed to confirm the phylogenetic position of two Coelogyne species based on the whole chloroplast genome sequences. Additionally, we developed eight new primers based on the two cp genomes' medium variable regions and evaluated the transferability to another 16 Coelogyne species. We constructed phylogenetic trees including 18 Coelogyne species and four outgroup species using the chloroplast fragments with the ML method. Our results showed that the cp genomes of C. fimbriata and C. ovalis contained a small single-copy region (18,839 and 18,851 bp, respectively) and a large single-copy region (87,606 and 87,759 bp, respectively), separated by two same-length inverted-repeat regions (26,675 bp in C. fimbriata and 26,715 bp C. ovalis, respectively). They all contained 86 protein-coding genes, 38 tRNA genes, and eight rRNA genes, revealing strong structure and gene content similarities. The phylogenetic analysis indicated a close relationship between the genera Coelogyne and Pleione. The newly developed primers revealed good transferability among the Coelogyne taxa and provided enough variable sites to distinguish C. fimbriata and C. ovalis. The two complete cp genomes and the eight new primers of Coelogyne provide new genomic data for further studies on phylogenomics, population genetics, and evolutionary history of Coelogyne taxa.