Complete chloroplast genome sequences of Hordeum vulgare, Sorghum bicolor and Agrostis stolonifera, and comparative analyses with other grass genomes

Temperature-Caused Changes in Raman Pattern and Protein Profiles of Winter Triticale (x Triticosecale, Wittm.) Field-Grown Seedlings

Climate change, which causes periods with relatively high temperatures in winter in Poland, can lead to a shortening or interruption of the cold hardening of crops. Previous research indicates that cold acclimation is of key importance in the process of acquiring cereal tolerance to stress factors. The objective of this work was to verify the hypothesis that both natural temperature fluctuations and the plant genotype influence the content of metabolites as well as proteins, including antioxidant enzymes and photosystem proteins. The research material involved four winter triticale genotypes, differing in their tolerance to stress under controlled conditions. The values of chlorophyll a fluorescence parameters and antioxidant activity were measured in their seedlings. Subsequently, the contribution of selected proteins was verified using specific antibodies. In parallel, the profiling of the contents of chlorophylls, carotenoids, phenolic compounds, and proteins was carried out by Raman spectroscopy. The obtained results indicate that a better PSII performance along with a higher photosystem II proteins content and thioredoxin reductase abundance were accompanied by a higher antioxidant activity in the field-grown triticale seedlings. The Raman studies showed that the cold hardening led to a variation in photosynthetic dyes and an increase in the phenolic to carotenoids ratio in all DH lines.

Chloroplast genomes of Caragana tibetica and Caragana turkestanica: structures and comparative analysis

BACKGROUND

The genus Caragana encompasses multiple plant species that possess medicinal and ecological value. However, some species of Caragana are quite similar in morphology, so identifying species in this genus based on their morphological characteristics is considerably complex. In our research, illumina paired-end sequencing was employed to investigate the genetic organization and structure of Caragana tibetica and Caragana turkestanica, including the previously published chloroplast genome sequence of 7 Caragana plants.

RESULTS

The lengths of C. tibetica and C. turkestanica chloroplast genomes were 128,433 bp and 129,453 bp, respectively. The absence of inverted repeat sequences in these two species categorizes them under the inverted repeat loss clade (IRLC). They encode 110 and 111 genes (4 /4 rRNA genes, 30 /31tRNA genes, and 76 /76 protein-coding genes), respectively. Comparison of the chloroplast genomes of C. tibetica and C. turkestanica with 7 other Caragana species revealed a high overall sequence similarity. However, some divergence was observed between certain intergenic regions (matK-rbcL, psbD-psbM, atpA-psbI, and etc.). Nucleotide diversity (π) analysis revealed the detection of five highly likely variable regions, namely rps2-atpI, accD-psaI-ycf4, cemA-petA, psbN-psbH and rpoA-rps11. Phylogenetic analysis revealed that C. tibetica's sister species is Caragana jubata, whereas C. turkestanica's closest relative is Caragana arborescens.

CONCLUSIONS

The present study provides worthwhile information about the chloroplast genomes of C. tibetica and C. turkestanica, which aids in the identification and classification of Caragana species.

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.

The complete plastome sequences of invasive weed Parthenium hysterophorus: genome organization, evolutionary significance, structural features, and comparative analysis

Parthenium hysterophorus, a globally widespread weed, poses a significant threat to agricultural ecosystems due to its invasive nature. We investigated the chloroplast genome of P. hysterophorus in this study. Our analysis revealed that the chloroplast genome of P. hysterophorus spans a length of 151,881 base pairs (bp). It exhibits typical quadripartite structure commonly found in chloroplast genomes, including inverted repeat regions (IR) of 25,085 bp, a small single copy (SSC) region of 18,052 bp, and a large single copy (LSC) region of 83,588 bp. A total of 129 unique genes were identified in P. hysterophorus chloroplast genomes, including 85 protein-coding genes, 36 tRNAs, and eight rRNAs genes. Comparative analysis of the P. hysterophorus plastome with those of related species from the tribe Heliantheae revealed both conserved structures and intriguing variations. While many structural elements were shared among the species, we identified a rearrangement in the large single-copy region of P. hysterophorus. Moreover, our study highlighted notable gene divergence in several specific genes, namely matK, ndhF, clpP, rps16, ndhA, rps3, and ndhD. Phylogenetic analysis based on the 72 shared genes placed P. hysterophorus in a distinct clade alongside another species, P. argentatum. Additionally, the estimated divergence time between the Parthenium genus and Helianthus (sunflowers) was approximately 15.1 million years ago (Mya). These findings provide valuable insights into the evolutionary history and genetic relationships of P. hysterophorus, shedding light on its divergence and adaptation over time.

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.

A splicing site change between exon 5 and 6 of the nuclear-encoded chloroplast-localized HvYGL8 gene results in reduced chlorophyll content and plant height in barley

The chloroplast is an important cellular organelle and metabolic hub, which is not only responsible for plant photosynthesis but is also involved in the de novo biosynthesis of pigments, fatty acids, and hormone metabolisms. Several genes that are responsible for rice leaf color variations have been reported to be directly or indirectly involved in chlorophyll biosynthesis and chloroplast development, whereas a few genes have been functionally confirmed to be responsible for leaf color changes in barley at the molecular level. In this study, we obtained a yellow leaf and dwarf ygl8 mutant from the progeny of Morex (a variety of barley) seeds treated with EMS. We performed bulked-segregant analysis (BSA) and RNA-seq analysis and targeted a UMP kinase encoding gene, YGL8, which generated a splicing site change between exon 5 and 6 of YGL8 due to a G to A single-nucleotide transition in the 5th exon/intron junction in the ygl8 mutant. The splicing site change between exon 5 and 6 of YGL8 had no effects on chloroplast subcellular localization but resulted in an additional loop in the UMP kinase domain, which might disturb the access of the substrates. On one hand, the splicing site change between exon 5 and 6 of YGL8 downregulated the transcriptional expression of chloroplast-encoded genes and chlorophyll-biosynthesis-related genes in a temperature-dependent manner in the ygl8 mutant. On the other hand, the downregulation of bioactive GA-biosynthesis-related GA20ox genes and cell-wall-cellulose-biosynthesis-related CesA genes was also observed in the ygl8 mutant, which led to a reduction in plant height. Our study will facilitate the understanding of the regulation of leaf color and plant height in barley.

Development of InDels markers for the identification of cytoplasmic male sterility in Sorghum by complete chloroplast genome sequences analysis

Cytoplasmic male sterility (CMS) is predominantly used for F1 hybrid breeding and seed production in Sorghum. DNA markers to distinguish between normal fertile (CMS-N) and sterile (CMS-S) male cytoplasm can facilitate F1 hybrid cultivar development in Sorghum breeding programs. In this study, the complete chloroplast (cp) genome sequences of CMS-S and Korean Sorghum cultivars were obtained using next-generation sequencing. The de novo assembled genome size of ATx623, the CMS-S line of the chloroplast, was 140,644bp. When compared to the CMS-S and CMS-N cp genomes, 19 single nucleotide polymorphisms (SNPs) and 142 insertions and deletions (InDels) were identified, which can be used for marker development for breeding, population genetics, and evolution studies. Two InDel markers with sizes greater than 20 bp were developed to distinguish cytotypes based on the copy number variation of lengths as 28 and 22 bp tandem repeats, respectively. Using the newly developed InDel markers with five pairs of CMS-S and their near isogenic maintainer line, we were able to easily identify their respective cytotypes. The InDel markers were further examined and applied to 1,104 plants from six Korean Sorghum cultivars to identify variant cytotypes. Additionally, the phylogenetic analysis of seven Sorghum species with complete cp genome sequences, including wild species, indicated that CMS-S and CMS-N contained Milo and Kafir cytotypes that might be hybridized from S. propinquum and S. sudanese, respectively. This study can facilitate F1 hybrid cultivar development by providing breeders with reliable tools for marker-assisted selection to breed desirable Sorghum varieties.

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.

Variation in mitogenome structural conformation in wild and cultivated lineages of sorghum corresponds with domestication history and plastome evolution

BACKGROUND

Mitochondria are organelles within eukaryotic cells that are central to the metabolic processes of cellular respiration and ATP production. However, the evolution of mitochondrial genomes (mitogenomes) in plants is virtually unknown compared to animal mitogenomes or plant plastids, due to complex structural variation and long stretches of repetitive DNA making accurate genome assembly more challenging. Comparing the structural and sequence differences of organellar genomes within and between sorghum species is an essential step in understanding evolutionary processes such as organellar sequence transfer to the nuclear genome as well as improving agronomic traits in sorghum related to cellular metabolism.

RESULTS

Here, we assembled seven sorghum mitochondrial and plastid genomes and resolved reticulated mitogenome structures with multilinked relationships that could be grouped into three structural conformations that differ in the content of repeats and genes by contig. The grouping of these mitogenome structural types reflects the two domestication events for sorghum in east and west Africa.

CONCLUSIONS

We report seven mitogenomes of sorghum from different cultivars and wild sources. The assembly method used here will be helpful in resolving complex genomic structures in other plant species. Our findings give new insights into the structure of sorghum mitogenomes that provides an important foundation for future research into the improvement of sorghum traits related to cellular respiration, cytonuclear incompatibly, and disease resistance.

SSR identification and phylogenetic analysis in four plant species based on complete chloroplast genome sequences

The effective utilization of traditional Chinese medicine (TCM) has been challenged by the difficulty to accurately distinguish between similar plant varieties. The stability and conservation of the chloroplast genome can aid in resolving genotypes. Previous studies using nuclear sequences and molecular markers have not effectively differentiated the species from related taxa, such as Machilus leptophylla, Hanceola exserta, Rubus bambusarum, and Rubus henryi. This study aimed to characterize the chloroplast genomes of these four plant species, and analyze their simple sequence repeats (SSRs) and phylogenetic positions. The results demonstrated the four chloroplast genomes consisted of 152.624 kb, 153.296 kb, 156.309 kb, and 158.953 kb in length, involving 124, 130, 129, and 131 genes, respectively. They also contained four specific regions with mononucleotide being the class with the most members. Moreover, these repeating types of SSR were various in individual class. Phylogenetic analysis showed that M. leptophylla was clustered with M. yunnanensis, and H. exserta was confirmed as belonging to the family Ocimeae. Additionally, R. bambusarum and R. henryi were grouped together but differed in their SSR features, indicating that they were not the same species. This research provides evidence for resolving species and contributes new genetic information for further studies.

Three Parts of the Plant Genome: On the Way to Success in the Production of Recombinant Proteins

Recombinant proteins are the most important product of current industrial biotechnology. They are indispensable in medicine (for diagnostics and treatment), food and chemical industries, and research. Plant cells combine advantages of the eukaryotic protein production system with simplicity and efficacy of the bacterial one. The use of plants for the production of recombinant proteins is an economically important and promising area that has emerged as an alternative to traditional approaches. This review discusses advantages of plant systems for the expression of recombinant proteins using nuclear, plastid, and mitochondrial genomes. Possibilities, problems, and prospects of modifications of the three parts of the genome in light of obtaining producer plants are examined. Examples of successful use of the nuclear expression platform for production of various biopharmaceuticals, veterinary drugs, and technologically important proteins are described, as are examples of a high yield of recombinant proteins upon modification of the chloroplast genome. Potential utility of plant mitochondria as an expression system for the production of recombinant proteins and its advantages over the nucleus and chloroplasts are substantiated. Although these opportunities have not yet been exploited, potential utility of plant mitochondria as an expression system for the production of recombinant proteins and its advantages over the nucleus and chloroplasts are substantiated.

What Do We Know about Barley miRNAs?

Plant miRNAs are powerful regulators of gene expression at the post-transcriptional level, which was repeatedly proved in several model plant species. miRNAs are considered to be key regulators of many developmental, homeostatic, and immune processes in plants. However, our understanding of plant miRNAs is still limited, despite the fact that an increasing number of studies have appeared. This systematic review aims to summarize our current knowledge about miRNAs in spring barley (Hordeum vulgare), which is an important agronomical crop worldwide and serves as a common monocot model for studying abiotic stress responses as well. This can help us to understand the connection between plant miRNAs and (not only) abiotic stresses in general. In the end, some future perspectives and open questions are summarized.

Complete Chloroplast Genome Sequence of the Endemic and Endangered Plant Dendropanax oligodontus: Genome Structure, Comparative and Phylogenetic Analysis

Dendropanax oligodontus, which belongs to the family Araliaceae, is an endemic and endangered species of Hainan Island, China. It has potential economic and medicinal value owing to the presence of phenylpropanoids, flavonoids, triterpenoids, etc. The analysis of the structure and characteristics of the D. oligodontus chloroplast genome (cpDNA) is crucial for understanding the genetic and phylogenetic evolution of this species. In this study, the cpDNA of D. oligodontus was sequenced for the first time using next-generation sequencing methods, assembled, and annotated. We observed a circular quadripartite structure comprising a large single-copy region (86,440 bp), a small single-copy region (18,075 bp), and a pair of inverted repeat regions (25,944 bp). The total length of the cpDNA was 156,403 bp, and the GC% was 37.99%. We found that the D. oligodontus chloroplast genome comprised 131 genes, with 86 protein-coding genes, 8 rRNA genes, and 37 tRNAs. Furthermore, we identified 26,514 codons, 13 repetitive sequences, and 43 simple sequence repeat sites in the D. oligodontus cpDNA. The most common amino acid encoded was leucine, with a strong A/T preference at the third position of the codon. The prediction of RNA editing sites in the protein-coding genes indicated that RNA editing was observed in 19 genes with a total of 54 editing sites, all of which involved C-to-T transitions. Finally, the cpDNA of 11 species of the family Araliaceae were selected for comparative analysis. The sequences of the untranslated regions and coding regions among 11 species were highly conserved, and minor differences were observed in the length of the inverted repeat regions; therefore, the cpDNAs were relatively stable and consistent among these 11 species. The variable hotspots in the genome included clpP, ycf1, rnK-rps16, rps16-trnQ, atpH-atpI, trnE-trnT, psbM-trnD, ycf3-trnS, and rpl32-trnL, providing valuable molecular markers for species authentication and regions for inferring phylogenetic relationships among them, as well as for evolutionary studies. Evolutionary selection pressure analysis indicated that the atpF gene was strongly subjected to positive environmental selection. Phylogenetic analysis indicated that D. oligodontus and Dendropanax dentiger were the most closely related species within the genus, and D. oligodontus was closely related to the genera Kalopanax and Metapanax in the Araliaceae family. Overall, the cp genomes reported in this study will provide resources for studying the genetic diversity and conservation of the endangered plant D. oligodontus, as well as resolving phylogenetic relationships within the family.

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.

The Plastome Sequences of Triticum sphaerococcum (ABD) and Triticum turgidum subsp. durum (AB) Exhibit Evolutionary Changes, Structural Characterization, Comparative Analysis, Phylogenomics and Time Divergence

The mechanism and course of Triticum plastome evolution is currently unknown; thus, it remains unclear how Triticum plastomes evolved during recent polyploidization. Here, we report the complete plastomes of two polyploid wheat species, Triticum sphaerococcum (AABBDD) and Triticum turgidum subsp. durum (AABB), and compare them with 19 available and complete Triticum plastomes to create the first map of genomic structural variation. Both T. sphaerococcum and T. turgidum subsp. durum plastomes were found to have a quadripartite structure, with plastome lengths of 134,531 bp and 134,015 bp, respectively. Furthermore, diploid (AA), tetraploid (AB, AG) and hexaploid (ABD, AGA^m) Triticum species plastomes displayed a conserved gene content and commonly harbored an identical set of annotated unique genes. Overall, there was a positive correlation between the number of repeats and plastome size. In all plastomes, the number of tandem repeats was higher than the number of palindromic and forward repeats. We constructed a Triticum phylogeny based on the complete plastomes and 42 shared genes from 71 plastomes. We estimated the divergence of Hordeum vulgare from wheat around 11.04-11.9 million years ago (mya) using a well-resolved plastome tree. Similarly, Sitopsis species diverged 2.8-2.9 mya before Triticum urartu (AA) and Triticum monococcum (AA). Aegilops speltoides was shown to be the maternal donor of polyploid wheat genomes and diverged ~0.2-0.9 mya. The phylogeny and divergence time estimates presented here can act as a reference framework for future studies of Triticum evolution.

Chloroplast Genes Are Involved in The Male-Sterility of K-Type CMS in Wheat

The utilization of crop heterosis can greatly improve crop yield. The sterile line is vital for the heterosis utilization of wheat (Triticum aestivum L.). The chloroplast genomes of two sterile lines and one maintainer were sequenced using second-generation high-throughput technology and assembled. The nonsynonymous mutated genes among the three varieties were identified, the expressed difference was further analyzed by qPCR, and finally, the function of the differentially expressed genes was analyzed by the barley stripe mosaic virus-induced gene silencing (BSMV-VIGS) method. A total of 16 genes containing 31 nonsynonymous mutations between K519A and 519B were identified. There were no base mutations in the protein-encoding genes between K519A and YS3038. The chloroplast genomes of 519B and K519A were closely related to the Triticum genus and Aegilops genus, respectively. The gene expression levels of the six selected genes with nonsynonymous mutation sites for K519A compared to 519B were mostly downregulated at the binucleate and trinucleate stages of pollen development. The seed setting rates of atpB-silenced or ndhH-silenced 519B plants by BSMV-VIGS method were significantly reduced. It can be concluded that atpB and the ndhH are likely to be involved in the reproductive transformation of 519B.

Evolutionary characteristics and phylogeny of cotton chloroplast tRNAs

The novel structural variations were identified in cotton chloroplast tRNAs and gene loss events were more obvious than duplications in chloroplast tRNAs. Transfer RNAs (tRNA) have long been believed an evolutionary-conserved molecular family, which play the key roles in the process of protein biosynthesis in plant life activities. In this study, we detected the evolutionary characteristics and phylogeny of chloroplast tRNAs in cotton plants, an economic and fibered important taxon in the world. We firstly annotated the chloroplast tRNAs of 27 Gossypium species to analyze their genetic composition, structural characteristics and evolution. Compared with the traditional view of evolutionary conservation of tRNA, some novel tRNA structural variations were identified in cotton plants. I.g., tRNA^Val-UAC and tRNA^Ile-GAU only contained one intron in the anti-condon loop region of tRNA secondary structure, respectively. In the variable region, some tRNAs contained a circle structure with a few nucleotides. Interestingly, the calculation result of free energy indicated that the variation of novel tRNAs contributed to the stability of tRNA structure. Phylogenetic analysis suggested that chloroplast tRNAs have evolved from multiple common ancestors, and the tRNA^Met seemed to be an ancestral tRNA, which can be duplicated and diversified to produce other tRNAs. The chloroplast tRNAs contained a group I intron in cotton plants, and the evolutionary analysis of introns indicated that group I intron of chloroplast tRNA originated from cyanobacteria. Analysis of gene duplication and loss events showed that gene loss events were more obvious than duplications in Gossypium chloroplast tRNAs. Additionally, we found that the rate of transition was higher than the ones of transversion in cotton chloroplast tRNAs. This study provided new insights into the structural characteristics and evolution of chloroplast tRNAs in cotton plants.

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.

Two NLR immune receptors acquired high-affinity binding to a fungal effector through convergent evolution of their integrated domain

A subset of plant NLR immune receptors carry unconventional integrated domains in addition to their canonical domain architecture. One example is rice Pik-1 that comprises an integrated heavy metal-associated (HMA) domain. Here, we reconstructed the evolutionary history of Pik-1 and its NLR partner, Pik-2, and tested hypotheses about adaptive evolution of the HMA domain. Phylogenetic analyses revealed that the HMA domain integrated into Pik-1 before Oryzinae speciation over 15 million years ago and has been under diversifying selection. Ancestral sequence reconstruction coupled with functional studies showed that two Pik-1 allelic variants independently evolved from a weakly binding ancestral state to high-affinity binding of the blast fungus effector AVR-PikD. We conclude that for most of its evolutionary history the Pik-1 HMA domain did not sense AVR-PikD, and that different Pik-1 receptors have recently evolved through distinct biochemical paths to produce similar phenotypic outcomes. These findings highlight the dynamic nature of the evolutionary mechanisms underpinning NLR adaptation to plant pathogens.

Complete chloroplast genome sequencing of Job's tears (Coix L.): genome structure, comparative analysis, and phylogenetic relationships

Job's tears, also known as adlay, is a valuable plant that has commonly been used in traditional Chinese medicine, as well as an edible food. Due to the lack of knowledge of its genetics and gaps in its evolutionary analysis, breeding of adlay has been hindered. Here, we report five complete chloroplast genomes of various species and varieties in the genus by Illumina sequencing, while their genome structure, comparative analysis, and phylogenetic relationships were conducted. Genome sizes ranged from 140,860 to 140,864 bp in length, GC contents were 38.43%, and genome architecture was of a typical quadripartite structure. We annotated 82～83 protein-coding genes and 46～47 non-coding RNA genes in each genome and they functionally associated with self-replication, photosynthesis, cytochrome synthesis and other unknown functions. Three codons that encoded tryptophan, arginine and leucine were used frequently at rates of 41.42, 37.98, and 32.28% respectively. The preferred codons consistently ended with A or T. A total of 146 simple sequence repeats (SSR), 9 insertions and deletions (InDels) and 143 single nucleotide polymorphisms (SNPs) were observed among genomes. The InDel and SNP variations were mostly distributed in intergenic regions. It confirmed that Coix, Sorghum, Saccharum, Zea, Tripsacum and Saccharum were closely genera and the genetic distance of Sorghum to Coix was closer than Zea to Coix. These results give us more insight into the evolution of Coix in a wide range of evolutionary studies.

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.

Identification of Vicia Species Native to South Korea Using Molecular and Morphological Characteristics

Recently, within the Fabaceae family, the Vicia genus has been recognized for its vital role in sustainable agriculture. Vicia species are economically important grain and forage crops. However, the presence of complex morphological characteristics makes identification and recognition of native species difficult. In this study, the possibility of using DNA barcoding regions (ITS2, matK, and rbcL) to distinguish among 19 Vicia taxa (59 accessions) found in South Korea was evaluated. The sequence alignment analysis revealed considerable nucleotide diversity (π) between the loci, in which ITS2 showed the highest mean interspecific distance, whereas there was no intraspecific variability among the barcode regions in 12 of the 19 taxa. Phylogenetic analysis of combined barcoding regions revealed well-resolved phylogeny with the highest species level discrimination. Combinations of barcode loci were also used in classification at the subgenera and section levels. The results revealed that the combined barcoding regions can be used effectively to differentiate the following species: Vicia angustifolia var. segetilis, Vicia bungei, Vicia villosa, Vicia cracca, Vicia dasycarpa, Vicia hirsuta, Vicia tetrasperma, Vicia amurensis, Vicia hirticalycina, and Vicia chosenensis. However, it is difficult to differentiate the species of Vicia unijuga, Vicia unijuga var. kaussanensis, Vicia linearifolia, Vicia unijuga f. angustifolia, Vicia nipponica, Vicia amoena, Vicia venosa var. cuspidata, Vicia pseudo-orobus, and Vicia japonica with the tested barcode regions. These species come under sect. Vicilla and are found to be closely related or species that have recently undergone speciation; thus, it has limitation to distinguish with recommended barcodes. Hence, to differentiate the unclassified species, 39 morphological characteristics were investigated, in which 16 useful characteristics were selected for efficient classification. Finally, the 16 selected morphological useful traits efficiently differentiated all the Vicia species. In conclusion, a combination of barcoding loci together with morphological characteristics of this study efficiently discriminated all the Korean Vicia species.

A high-quality genome assembly of the eggplant provides insights into the molecular basis of disease resistance and chlorogenic acid synthesis

The eggplant (Solanum melongena L.) is one of the most important Solanaceae crops, ranking third for total production and economic value in its genus. Herein, we report a high-quality, chromosome-scale eggplant reference genome sequence of 1155.8 Mb, with an N50 of 93.9 Mb, which was assembled by combining PacBio long reads and Hi-C sequencing data. Repetitive sequences occupied 70.1% of the assembly length, and 35,018 high-confidence protein-coding genes were annotated based on multiple sources. Comparative analysis revealed 646 species-specific families and 364 positive selection genes, conferring distinguishing traits on the eggplant. We performed genome-wide comparative identification of disease resistance genes and discovered an expanded gene family of bacterial spot resistance in eggplant and pepper, but not in tomato and potato. The genes involved in chlorogenic acid synthesis were comprehensively characterized. Highly similar chromosomal distribution patterns of polyphenol oxidase genes were observed in the eggplant, tomato, and potato genomes. The eggplant reference genome sequence will not only facilitate evolutionary studies of the Solanaceae but also facilitate their breeding and improvement.

Escalation in the host-pathogen arms race: A host resistance response corresponds to a heightened bacterial virulence response

The zig-zag model of host-pathogen interaction describes the relative strength of defense response across a spectrum of pathogen-induced plant phenotypes. A stronger defense response results in increased resistance. Here, we investigate the strength of pathogen virulence during disease and place these findings in the context of the zig-zag model. Xanthomonas vasicola pv. holcicola (Xvh) causes sorghum bacterial leaf streak. Despite being widespread, this disease has not been described in detail at the molecular level. We divided diverse sorghum genotypes into three groups based on disease symptoms: water-soaked lesions, red lesions, and resistance. Bacterial growth assays confirmed that these three phenotypes represent a range of resistance and susceptibility. To simultaneously reveal defense and virulence responses across the spectrum of disease phenotypes, we performed dual RNA-seq on Xvh-infected sorghum. Consistent with the zig-zag model, the expression of plant defense-related genes was strongest in the resistance interaction. Surprisingly, bacterial virulence genes related to the type III secretion system (T3SS) and type III effectors (T3Es) were also most highly expressed in the resistance interaction. This expression pattern was observed at multiple time points within the sorghum-Xvh pathosystem. Further, a similar expression pattern was observed in Arabidopsis infected with Pseudomonas syringae for effector-triggered immunity via AvrRps4 but not AvrRpt2. Specific metabolites were able to repress the Xvh virulence response in vitro and in planta suggesting a possible signaling mechanism. Taken together, these findings reveal multiple permutations of the continually evolving host-pathogen arms race from the perspective of host defense and pathogen virulence responses.

The arms race between plants and pathogens is a complex process. To dissect the plant defense and pathogen virulence responses simultaneously, we used sorghum and Xanthomonas vasicola pv. holcicola, as a model pathosystem. We performed dual RNA-seq on infected sorghum with a range of disease phenotypes. Our characterization of this pathosystem demonstrates that genes related to the plant defense and pathogen virulence responses are most highly induced during a resistance interaction. We observed a similar pattern of escalation in Arabidopsis infected with Pseudomonas syringae. These observations support a conceptual model of fluidity between the different stages of plant immunity and pathogen virulence.

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.

The Complete Plastid Genome of Artocarpus camansi: A High Degree of Conservation of the Plastome Structure in the Family Moraceae

Understanding the plastid genome is extremely important for the interpretation of the genetic mechanisms associated with essential physiological and metabolic functions, the identification of possible marker regions for phylogenetic or phylogeographic analyses, and the elucidation of the modes through which natural selection operates in different regions of this genome. In the present study, we assembled the plastid genome of Artocarpus camansi, compared its repetitive structures with Artocarpus heterophyllus, and searched for evidence of synteny within the family Moraceae. We also constructed a phylogeny based on 56 chloroplast genes to assess the relationships among three families of the order Rosales, that is, the Moraceae, Rhamnaceae, and Cannabaceae. The plastid genome of A. camansi has 160,096 bp, and presents the typical circular quadripartite structure of the Angiosperms, comprising a large single copy (LSC) of 88,745 bp and a small single copy (SSC) of 19,883 bp, separated by a pair of inverted repeat (IR) regions each with a length of 25,734 bp. The total GC content was 36.0%, which is very similar to Artocarpus heterophyllus (36.1%) and other moraceous species. A total of 23,068 codons and 80 SSRs were identified in the A. camansi plastid genome, with the majority of the SSRs being mononucleotide (70.0%). A total of 50 repeat structures were observed in the A. camansi plastid genome, in contrast with 61 repeats in A. heterophyllus. A purifying selection signal was found in 70 of the 79 protein-coding genes, indicating that they have all been highly conserved throughout the evolutionary history of the genus. The comparative analysis of the structural characteristics of the chloroplast among different moraceous species found a high degree of similarity in the sequences, which indicates a highly conserved evolutionary model in these plastid genomes. The phylogenetic analysis also recovered a high degree of similarity between the chloroplast genes of A. camansi and A. heterophyllus, and reconfirmed the hypothesis of the intense conservation of the plastome in the family Moraceae.

Barley's Second Spring as A Model Organism for Chloroplast Research

Barley (Hordeum vulgare) has been widely used as a model crop for studying molecular and physiological processes such as chloroplast development and photosynthesis. During the second half of the 20th century, mutants such as albostrians led to the discovery of the nuclear-encoded, plastid-localized RNA polymerase and the retrograde (chloroplast-to-nucleus) signalling communication pathway, while chlorina-f2 and xantha mutants helped to shed light on the chlorophyll biosynthetic pathway, on the light-harvesting proteins and on the organization of the photosynthetic apparatus. However, during the last 30 years, a large fraction of chloroplast research has switched to the more “user-friendly” model species Arabidopsis thaliana, the first plant species whose genome was sequenced and published at the end of 2000. Despite its many advantages, Arabidopsis has some important limitations compared to barley, including the lack of a real canopy and the absence of the proplastid-to-chloroplast developmental gradient across the leaf blade. These features, together with the availability of large collections of natural genetic diversity and mutant populations for barley, a complete genome assembly and protocols for genetic transformation and gene editing, have relaunched barley as an ideal model species for chloroplast research. In this review, we provide an update on the genomics tools now available for barley, and review the biotechnological strategies reported to increase photosynthesis efficiency in model species, which deserve to be validated in barley.

Assessment of Gene Flow to Wild Relatives and Nutritional Composition of Sugarcane in Brazil

The commercial release of genetically modified organisms (GMO) requires a prior environmental and human/animal health risk assessment. In Brazil, the National Biotechnology Technical Commission (CTNBio) requires a survey of the area of natural occurrence of wild relatives of the GMO in the Brazilian ecosystems to evaluate the possibility of introgressive hybridization between sexually compatible species. Modern sugarcane cultivars, the focus of this study, derive from a series of hybridization and backcrossing events among Saccharum species. The so-called "Saccharum broad sense" group includes around 40 species from a few genera, including Erianthus, found in various tropical regions, particularly South-Eastern Asia. In Brazil, three native species, originally considered to belong to Erianthus, were reclassified as S. angustifolium (Nees) Trin., S. asperum (Nees) Steud., and S. villosum Steud., based on inflorescence morphology. Thus, we have investigated the potential occurrence of gene flow among the Brazilian Saccharum native species and commercial hybrids as a requisite for GMO commercial release. A comprehensive survey was carried out to map the occurrence of the three native Saccharum species in Brazil, concluding that they are sympatric with sugarcane cultivation only from around 14°S southwards, which precludes most Northeastern sugarcane-producing states from undergoing introgression. Based on phenology, we concluded that the Brazilian Saccharum species are unable to outcross naturally with commercial sugarcane since the overlap between the flowering periods of sugarcane and the native species is limited. A phylogenomic reconstruction based on the full plastid genome sequence showed that the three native Saccharum species are the taxa closest to sugarcane in Brazil, being closer than introduced Erianthus or Miscanthus. A 2-year study on eight nutritional composition traits of the 20 main sugarcane cultivars cultivated in Brazil was carried out in six environments. The minimum and maximum values obtained were, in percent: moisture (62.6-82.5); sucrose (9.65-21.76); crude fiber (8.06-21.03); FDN (7.20-20.68); FDA (4.55-16.90); lipids (0.06-1.59); ash (0.08-2.67); and crude protein (0.18-1.18). Besides a considerable amount of genetic variation and plastic responses, many instances of genotype-by-environment interaction were detected.

Gene transfer and nucleotide sequence evolution by Gossypium cytoplasmic genomes indicates novel evolutionary characteristics

The DNA fragments transferred among cotton cytoplasmic genomes are highly differentiated. The wild D group cotton species have undergone much greater evolution compared with cultivated AD group. Cotton (Gossypium spp.) is one of the most economically important fiber crops worldwide. Gene transfer, nucleotide evolution, and the codon usage preferences in cytoplasmic genomes are important evolutionary characteristics of high plants. In this study, we analyzed the nucleotide sequence evolution, codon usage, and transfer of cytoplasmic DNA fragments in Gossypium chloroplast (cp) and mitochondrial (mt) genomes, including the A genome group, wild D group, and cultivated AD group of cotton species. Our analyses indicated that the differences in the length of transferred cytoplasmic DNA fragments were not significant in mitochondrial and chloroplast sequences. Analysis of the transfer of tRNAs found that trnQ and nine other tRNA genes were commonly transferred between two different cytoplasmic genomes. The Codon Adaptation Index values showed that Gossypium cp genomes prefer A/T-ending codons. Codon preference selection was higher in the D group than the other two groups. Nucleotide sequence evolution analysis showed that intergenic spacer sequences were more variable than coding regions and nonsynonymous mutations were clearly more common in cp genomes than mt genomes. Evolutionary analysis showed that the substitution rate was much higher in cp genomes than mt genomes. Interestingly, the D group cotton species have undergone much faster evolution compared with cultivated AD groups, possibly due to the selection and domestication of diverse cotton species. Our results demonstrate that gene transfer and differential nucleotide sequence evolution have occurred frequently in cotton cytoplasmic genomes.

The study of organelle DNA variability in alloplasmic barley lines in the NGS era

Alloplasmic lines are a suitable model for studying molecular coevolution and interrelations between genetic systems of plant cells. Whole chloroplast (cp) and mitochondrial (mt) genome sequences were obtained by the MiSeq System (Illumina). Organelle DNA samples were prepared from a set of 12 alloplasmic barley lines with different cytoplasms of Hordeum vulgare ssp. spontaneum and H. vulgare ssp. vulgare, as well as from their paternal varieties. A bioinformatic approach for analysis of NGS data obtained on an organellar DNA mix has been developed and verified. A comparative study of Hordeum organelle genomes’ variability and disposition of polymorphic loci was conducted. Eight types of chloroplast DNA and 5 types of mitochondrial DNA were distinguished for the barley sample set examined. These results were compared with the previous data of a restriction fragment length polymorphism (RFLP) study of organelle DNAs for the same material. Formerly established data about a field evaluation of alloplasmic barley lines were revised in the light of information about organelle genomes gained after NGS. Totally 17 polymorphic loci were found at exons of chloroplast genomes. Seven of the SNPs were located in the genes of the Ndh complex. The nonsynonymous changes of nucleotides were detected in the matK, rpoC1, ndhK, ndhG and infA genes. Some of the SNPs detected are very similar in codon position and in the type of amino acid substitution to the places where RNA editing can occur. Thus, these results outline new perspectives for the future study of nuclear-cytoplasmic interactions in alloplasmic lines.

The draft genome of a wild barley genotype reveals its enrichment in genes related to biotic and abiotic stresses compared to cultivated barley

Wild barley (Hordeum spontaneum) is the progenitor of cultivated barley (Hordeum vulgare) and provides a rich source of genetic variations for barley improvement. Currently, the genome sequences of wild barley and its differences with cultivated barley remain unclear. In this study, we report a high-quality draft assembly of wild barley accession (AWCS276; henceforth named as WB1), which consists of 4.28 Gb genome and 36 395 high-confidence protein-coding genes. BUSCO analysis revealed that the assembly included full lengths of 95.3% of the 956 single-copy plant genes, illustrating that the gene-containing regions have been well assembled. By comparing with the genome of the cultivated genotype Morex, it is inferred that the WB1 genome contains more genes involved in resistance and tolerance to biotic and abiotic stresses. The presence of the numerous WB1-specific genes indicates that, in addition to enhance allele diversity for genes already existing in the cultigen, exploiting the wild barley taxon in breeding should also allow the incorporation of novel genes. Furthermore, high levels of genetic variation in the pericentromeric regions were detected in chromosomes 3H and 5H between the wild and cultivated genotypes, which may be the results of domestication. This H. spontaneum draft genome assembly will help to accelerate wild barley research and be an invaluable resource for barley improvement and comparative genomics research.

The complete chloroplast genome of Dendrobium nobile, an endangered medicinal orchid from north-east India and its comparison with related Dendrobium species

The medicinal orchid genus Dendrobium belonging to the Orchidaceae family is a huge genus comprising about 800-1,500 species. To better illustrate the species status in the genus Dendrobium, a comparative analysis of 33 available chloroplast genomes retrieved from NCBI RefSeq database was compared with that of the first complete chloroplast genome of D. nobile from north-east India based on next-generation sequencing methods (Illumina HiSeq 2500-PE150). Our results provide comparative chloroplast genomic information for taxonomical identification, alignment-free phylogenomic inference and other statistical features of Dendrobium plastomes, which can also provide valuable information on their mutational events and sequence divergence.

Comparative Analysis of two Sugarcane Ancestors Saccharum officinarum and S. spontaneum based on Complete Chloroplast Genome Sequences and Photosynthetic Ability in Cold Stress

Polyploid Saccharum with complex genomes hindered the progress of sugarcane improvement, while their chloroplast genomes are much smaller and simpler. Chloroplast (cp), the vital organelle, is the site of plant photosynthesis, which also evolves other functions, such as tolerance to environmental stresses. In this study, the cp genome of two sugarcane ancestors Saccharum officinarum and S. spontaneum were sequenced, and genome comparative analysis between these two species was carried out, together with the photosynthetic ability. The length is 141,187 bp for S. officinarum and that is 7 bp longer than S. spontaneum, with the same GC content (38.44%) and annotated gene number (134), 13 with introns among them. There is a typical tetrad structure, including LSC, SSC, IRb and IRa. Of them, LSC and IRa/IRb are 18 bp longer and 6 bp shorter than those in S. spontaneum (83,047 bp and 22,795 bp), respectively, while the size of SSC is same (12,544 bp). Five genes exhibit contraction and expansion at the IR junctions, but only one gene ndhF with 29 bp expansion at the border of IRb/SSC. Nucleotide diversity (Pi) based on sliding window analysis showed that the single copy and noncoding regions were more divergent than IR- and coding regions, and the variant hotspots trnG-trnM, psbM-petN, trnR-rps14, ndhC-trnV and petA-psbJ in the LSC and trnL-ccsA in the SSC regions were detected, and petA-psbJ with the highest divergent value of 0.01500. Genetic distances of 65 protein genes vary from 0.00000 to 0.00288 between two species, and the selective pressure on them indicated that only petB was subjected to positive selection, while more genes including rpoC2, rps3, ccsA, ndhA, ndhA, psbI, atpH and psaC were subjected to purifying or very strong purifying selection. There are larger number of codons in S. spontaneum than that in S. officinarum, while both species have obvious codon preference and the codons with highest-(AUG) and lowest frequency (AUA) are same. Whilst, the most abundant amino acid is leucine in both S. officinarum and S. spontaneum, with number of 2175 (10.88% of total) and 2228 (10.90% of total) codons, respectively, and the lowest number is cysteine, with only 221 (1.105%) and 224 (1.096%), respectively. Protein collinearity analysis showed the high collinearity though several divergences were present in cp genomes, and identification of simple sequence repeats (SSRs) were included in this study. In addition, in order to compare cold tolerance and explore the expanding function of this environmental stress, the chlorophyll relative content (SPAD) and chlorophyll fluorescence Fv/Fm were measured. The significantly higher SPAD were observed in S. spontaneum than those in S. officinarum, no matter what the control conditions, exposure to low temperature or during recovery, and so was for Fv/Fm under exposure to low temperature, together with higher level of SPAD in S. spontaneum in each measurement. Aforementioned results suggest much stronger photosynthetic ability and cold tolerance in S. spontaneum. Our findings build a foundation to investigate the biological mechanism of two sugarcane ancestor chloroplasts and retrieve reliable molecular resources for phylogenetic and evolutionary studies, and will be conducive to genetic improvement of photosynthetic ability and cold resistance in modern sugarcane.

The rpl23 gene and pseudogene are hotspots of illegitimate recombination in barley chloroplast mutator seedlings

Previously, through a TILLING (Targeting Induced Local Lesions in Genomes) approach applied on barley chloroplast mutator (cpm) seedlings a high frequency of polymorphisms in the rpl23 gene was detected. All the polymorphisms corresponded to five differences already known to exist in nature between the rpl23 gene located in the inverted repeats (IRs) and the rpl23 pseudogene located in the large single copy region (LSC). In this investigation, polymorphisms in the rpl23 gene were verified and besides, a similar situation was found for the pseudogene in cpm seedlings. On the other hand, no polymorphisms were found in any of those loci in 40 wild type barley seedlings. Those facts and the independent occurrence of polymorphisms in the gene and pseudogene in individual seedlings suggest that the detected polymorphisms initially arose from gene conversion between gene and pseudogene. Moreover, an additional recombination process involving small recombinant segments seems to occur between the two gene copies as a consequence of their location in the IRs. These and previous results support the hypothesis that the CPM protein is a component of the plastome mismatch repair (MMR) system, whose failure of the anti-recombination activity results in increased illegitimate recombination between the rpl23 gene and pseudogene.

Leaf Variegation and Impaired Chloroplast Development Caused by a Truncated CCT Domain Gene in albostrians Barley

Chloroplasts fuel plant development and growth by converting solar energy into chemical energy. They mature from proplastids through the concerted action of genes in both the organellar and the nuclear genome. Defects in such genes impair chloroplast development and may lead to pigment-deficient seedlings or seedlings with variegated leaves. Such mutants are instrumental as tools for dissecting genetic factors underlying the mechanisms involved in chloroplast biogenesis. Characterization of the green-white variegated albostrians mutant of barley (Hordeum vulgare) has greatly broadened the field of chloroplast biology, including the discovery of retrograde signaling. Here, we report identification of the ALBOSTRIANS gene HvAST (also known as Hordeum vulgare CCT Motif Family gene 7, HvCMF7) by positional cloning as well as its functional validation based on independently induced mutants by Targeting Induced Local Lesions in Genomes (TILLING) and RNA-guided clustered regularly interspaced short palindromic repeats-associated protein 9 endonuclease-mediated gene editing. The phenotypes of the independent HvAST mutants imply residual activity of HvCMF7 in the original albostrians allele conferring an imperfect penetrance of the variegated phenotype even at homozygous state of the mutation. HvCMF7 is a homolog of the Arabidopsis (Arabidopsis thaliana) CONSTANS, CO-like, and TOC1 (CCT) Motif transcription factor gene CHLOROPLAST IMPORT APPARATUS2, which was reported to be involved in the expression of nuclear genes essential for chloroplast biogenesis. Notably, in barley we localized HvCMF7 to the chloroplast, without any clear evidence for nuclear localization.

Plastome phylogenomics and characterization of rare genomic changes as taxonomic markers in plastome groups 1 and 2 Poeae (Pooideae; Poaceae)

A phylogenomic analysis of 42 complete plastid genomes (plastomes), including 16 that were newly sequenced, was conducted. Plastomes were sampled from 19 subtribes of Pooideae, to investigate relationships within and between Chloroplast Group 1 (Aveneae) and Group 2 (Poeae) species. Two data partitions: complete plastomes, and a combined plastome and rare genomic change (RGC) data matrix, were analyzed. Overall, 156 non-ambiguous RGC were identified, of which homology was inferred for 38 RGC. Among the 38 RGC identified, six were synapomorphic among the Group 1 subtribes: Aveninae, Agrostidinae, and Anthoxanthinae, (Phalaridinae + Torreyochloinae), and 27 were synapomorphic among the Group 2 subtribes: Loliinae, (Ammochloinae + Parapholiinae + Dactylidinae), Parapholiinae, Dactylidinae, Poinae, and Coleanthinae. Four RGC were determined to be homoplasious in Groups 1 and 2. Two other RGC originated through intrastrand deletion events. The remaining RGC events likely originated through recombination given their size and lack of sequence evidence for other types of mutations. This study also determined that relationships between taxa, even those only weakly supported in previous studies, could be inferred with strong support when utilizing complete plastomes.

Comparative Chloroplast Genomes of Sorghum Species: Sequence Divergence and Phylogenetic Relationships

Sorghum comprises 31 species that exhibit considerable morphological and ecological diversity. The phylogenetic relationships among Sorghum species still remain unresolved due to lower information on the traditional DNA markers, which provides a limited resolution for identifying Sorghum species. In this study, we sequenced the complete chloroplast genomes of Sorghum sudanense and S. propinquum and analyzed the published chloroplast genomes of S. bicolor and S. timorense to retrieve valuable chloroplast molecular resources for Sorghum. The chloroplast genomes ranged in length from 140,629 to 140,755 bp, and their gene contents, gene orders, and GC contents were similar to those for other Poaceae species but were slightly different in the number of SSRs. Comparative analyses among the four chloroplast genomes revealed 651 variable sites, 137 indels, and nine small inversions. Four highly divergent DNA regions (rps16-trnQ, trnG-trnM, rbcL-psaI, and rps15-ndhF), which were suitable for phylogenetic and species identification, were detected in the Sorghum chloroplast genomes. A phylogenetic analysis strongly supported that Sorghum is a monophyletic group in the tribe Andropogoneae. Overall, the genomic resources in this study could provide potential molecular markers for phylogeny and species identification in Sorghum.

Genomic and evolutionary aspects of chloroplast tRNA in monocot plants

BACKGROUND

Chloroplasts are one of the most indispensable organelles that make life forms on the earth possible by their capacity to photosynthesize. These organelles possess a circular genome with a number of coding genes responsible for self-regulation. tRNAs are an important evolutionary-conserved gene family that are responsible for protein translation. However, within the chloroplast genome, tRNA machinery are poorly understood.

RESULTS

In the present study, the chloroplast genome of six monocot plants, Oryza nivara (NC_005973), Oryza sativa (NC_001320), Sachharum officinarum (NC_006084), Sorghum bicolor (NC_008602), Triticum aestivum (NC_002762), and Zea mays (NC_001666) were downloaded and analyzed to identify tRNA sequences. Further analysis of the tRNA sequences in the chloroplast genomes of the monocot plants resulted in the identification of several novel features. The length of tRNAs in the chloroplast genome of the monocot plants ranged from 59 to 155 nucleotides. Pair-wise sequence alignment revealed the presence of a conserved A-C-x-U-A-x-U-A-x-U-x5-U-A-A nucleotide consensus sequence. In addition, the tRNAs in chloroplast genomes of the monocot plants also contain 21-28 anti-codons against 61 sense codons in the genome. They also contain a group I intron and a C-A-U anti-codon for tRNAIle, which is a common anti-codon of tRNAMet. Evolutionary analysis indicates that tRNAs in the chloroplast genome have evolved from multiple common ancestors, and tRNAMet appears to be the ancestral tRNA that underwent duplication and diversification to give rise to other tRNAs.

CONCLUSION

The results obtained from the study of chloroplast tRNA will greatly help to increase our understanding of tRNA biology at a new level. Functional studies of the reported novel aspects of the chloroplast tRNA of the monocot plants will greatly help to decipher their roles in diverse cellular processes.

The Plastid Genome of Deschampsia cespitosa (Poaceae)

Plastid genome analysis of non-model organisms provides valuable information for basic research e.g., molecular evolutionary genomics, phylogeny and phylogeography. Deschampsia cespitosa is the most widespread species of the genus and it is a common grass that is found across Eurasia and North America. Scattered populations in regions of appropriate ecological conditions are also found in Australia, New Zealand and southern South America, where it is sympatric with D. antarctica. We analyzed the plastid genome of a sample of Deschampsia cespitosa of the Austrian Alps using high-throughput sequencing. The plastid (cp) genome shows the typical quadripartite structure with a length of 135,340 bp, comprising a large single-copy (LSC) region of 79,992 bp, a small single-copy (SSC) region of 12,572 bp and two inverted repeats (IR) regions of 21,388 bp each. It contains 115 genes, including 85 protein-coding genes, four ribosomal RNA genes and 30 transfer RNA genes. The GC content (%), number of repeats and microsatellites, RNA editing sites and codon usage were highly similar to those of D. antarctica. The results of this present study highlight the extremely conserved nature of the cp genome in this group, since the comparison involved individuals separated by about 13,000 km, from the Alps to Antarctica.

Green Biotechnology: A Brief Update on Plastid Genome Engineering

Plant genetic engineering has become an inevitable tool in the molecular breeding of crops. Significant progress has been made in the generation of novel plastid transformation vectors and optimized transformation protocols. There are several advantages of plastid genome engineering over conventional nuclear transformation. Some of the advantages include multigene engineering by expression of biosynthetic pathway genes as operons, extremely high-level expression of protein accumulation, lack of transgene silencing, etc. Transgene containment owing to maternal inheritance is another important advantage of plastid genome engineering. Chloroplast genome modification usually results in alteration of several thousand plastid genome copies in a cell. Several therapeutic proteins, edible vaccines, antimicrobial peptides, and industrially important enzymes have been successfully expressed in chloroplasts so far. Here, we critically recapitulate the latest developments in plastid genome engineering. Latest advancements in plastid genome sequencing are briefed. In addition, advancement of extending the toolbox for plastid engineering for selected applications in the area of molecular farming and production of industrially important enzyme is briefed.

Complete chloroplast genome sequence of Dryopteris fragrans (L.) Schott and the repeat structures against the thermal environment

Dryopteris fragrans (L.) Schott is a fern growing on the surface of hot rocks and lava. It is exposed to sunlight directly and bears local hot environment. We sequenced the complete nucleotide sequence of its chloroplast (cp) genome. The cp genome was 151,978 bp in length, consisting of a large single-copy region (85,332 bp), a small single-copy region (31,947 bp) and a pair of inverted repeats (17,314 bp). The cp genome contained 112 genes and 345 RNA editing sites in protein-coding genes. Simple sequence repeats (SSRs) and long repeat structure pairs (30–55 bp) were identified. The number and percent of repeat structures are extremely high in ferns. Thermal denaturation experiments showed its cp genome to have numerous, dispersed and high GC percent repeat structures, which conferred the strongest thermal stability. This repeat-heavy genome may provide the molecular basis of how D. fragrans cp survives its hot environment.

The complete chloroplast genome of the endangered plant Paphiopedilum wardii (Orchidaceae)

Chloroplast (cp) genome sequences provide a valuable source for phylogenetic analysis, which becomes a popular tool for population and phylogeny in a recent report. Here, the complete chloroplast genome of the Paphiopedilum wardii has been reconstructed from the whole-genome Illumina sequencing data. The circular genome is 169,820 bp in size and comprises a pair of inverted repeat (IR) regions of 33,641 bp each, a large single-copy (LSC) region of 98,872 bp and a small single-copy (SSC) region of 3666 bp. The total GC content is 33.8%, while the corresponding values of the IR regions, LSC, and SSC and are 38.7%, 31.7% and 26.7%, respectively. The chloroplast genome contains 151 genes, including 91 protein-coding genes, eight ribosomal RNA genes and 52 transfer RNA genes. Two genes (ndhA and ndhB) are pseudogenized or lost in its cpDNA. The maximum-likelihood phylogenetic analysis showed a close relationship with P. armeniacum in Orchidaceae. Our findings provide useful information for phylogenetic and evolutionary research of Paphiopedilum species.

First reported chloroplast genome sequence of Punica granatum (cultivar Helow) from Jabal Al-Akhdar, Oman: phylogenetic comparative assortment with Lagerstroemia

Pomegranate (Punica granatum L.) is one of the oldest known edible fruits. It has grown in popularity and is a profitable fruit crop due to its attractive features including a bright red appearance and its biological activities. Scientific exploration of the genetics and evolution of these beneficial traits has been hampered by limited genomic information. In this study, we sequenced the complete chloroplast (cp) genome of the native P. granatum (cultivar Helow) cultivated in the mountains of Jabal Al-Akhdar, Oman. The results revealed a P. granatum cp genome length of 158,630 bp, characterized by a relatively conserved structure containing 2 inverted repeat regions of 25,466 bp, an 18,686 bp small single copy regions, and an 89,015 bp large single copy region. The 86 protein-coding genes included 37 transfer RNA genes and 8 ribosomal RNA genes. Comparison of the P. granatum whole cp genome with seven Lagerstroemia species revealed an overall high degree of sequence similarity with divergence among intergenic spacers. The location, distribution, and divergence of repeat sequences and shared genes of the Punica and Lagerstroemia species were highly similar. Analyses of nucleotide substitution, insertion/deletions, and highly variable regions in these cp genomes identified potential plastid markers for taxonomic and phylogenetic studies in Myrtales. A phylogenetic study of the cp genomes and 76 shared coding regions generated similar cladograms. The complete cp genome of P. granatum will aid in taxonomical studies of the family Lythraceae.

Complete Chloroplast Genome Sequence of Broomcorn Millet (Panicum miliaceum L.) and Comparative Analysis with Other Panicoideae Species

Broomcorn millet (Panicum miliaceum L.) is one of the earliest domesticated cereals worldwide, holding significant agricultural, historical, and evolutionary importance. However, our genomic knowledge of it is rather limited at present, hampering further genetic and evolutionary studies. Here, we sequenced and assembled the chloroplast genome (cp) of broomcorn millet and compared it with five other Panicoideae species. Results showed that the cp genome of broomcorn millet was 139,826 bp in size, with a typical quadripartite structure. In total, 108 genes were annotated and 18 genes were duplicated in the IR (inverted region) region, which was similar to other Panicoideae species. Comparative analysis showed a rather conserved genome structure between them, with three common regions. Furthermore, RNA editing, codon usage, and expansion of the IR, as well as simple sequence repeat (SSR) elements, were systematically investigated and 13 potential DNA markers were developed for Panicoideae species identification. Finally, phylogenetic analysis implied that broomcorn millet was a sister species to Panicum virgatum within the tribe Paniceae, and supported a monophyly of the Panicoideae. This study has reported for the first time the genome organization, gene content, and structural features of the chloroplast genome of broomcorn millet, which provides valuable information for genetic and evolutionary studies in the genus Panicum and beyond.

First chloroplast genomics study of Phoenix dactylifera (var. Naghal and Khanezi): A comparative analysis

Date palm (Phoenix dactylifera L.) is one of the oldest fruit crops in the arid regions of the Middle East. However, little information is available regarding its plastid genomes. In this study, we sequenced the chloroplast (cp) genomes of two economically important but genomically unexplored date palm cultivars of Phoenix dactylifera (var. Naghal and Khanezi). The data assembly and genome annotation revealed a typical quadripartite structure similar to Arecaceae, and the genome sizes of Naghal and Khanezi were 158,210 bp and 158,211 bp, respectively. Structurally, both cp genomes were comprised of four regions: a pair of inverted repeats (27,273 bp for Khanezi and for Naghal 27,272 bp), a large single-copy region (86,090 bp and 86,092 bp) and a small single-copy region (17,575 bp and 17,574 bp). Both genomes had 138 representative genes, whereas 227 and 229 randomly distributed microsatellites were also observed in Khanezi and Naghal, respectively. Phylogenetic analysis based on the whole cp genomes and 68 shared genes showed identical phylogenetic trees of Khanezi and Naghal forming clades with Khalas and Aseel cultivars, respectively. The current study showed detailed comparative cp genome analysis, which could be essential for broader population genetics and molecular studies of these four date palm cultivars.

Comparative plastome genomics and phylogenomics of Brachypodium: flowering time signatures, introgression and recombination in recently diverged ecotypes

Few pan-genomic studies have been conducted in plants, and none of them have focused on the intraspecific diversity and evolution of their plastid genomes. We address this issue in Brachypodium distachyon and its close relatives B. stacei and B. hybridum, for which a large genomic data set has been compiled. We analyze inter- and intraspecific plastid comparative genomics and phylogenomic relationships within a family-wide framework. Major indel differences were detected between Brachypodium plastomes. Within B. distachyon, we detected two main lineages, a mostly Extremely Delayed Flowering (EDF+) clade and a mostly Spanish (S+) - Turkish (T+) clade, plus nine chloroplast capture and two plastid DNA (ptDNA) introgression and micro-recombination events. Early Oligocene (30.9 million yr ago (Ma)) and Late Miocene (10.1 Ma) divergence times were inferred for the respective stem and crown nodes of Brachypodium and a very recent Mid-Pleistocene (0.9 Ma) time for the B. distachyon split. Flowering time variation is a main factor driving rapid intraspecific divergence in B. distachyon, although it is counterbalanced by repeated introgression between previously isolated lineages. Swapping of plastomes between the three different genomic groups, EDF+, T+, S+, probably resulted from random backcrossing followed by stabilization through selection pressure.

A 250 plastome phylogeny of the grass family (Poaceae): topological support under different data partitions

The systematics of grasses has advanced through applications of plastome phylogenomics, although studies have been largely limited to subfamilies or other subgroups of Poaceae. Here we present a plastome phylogenomic analysis of 250 complete plastomes (179 genera) sampled from 44 of the 52 tribes of Poaceae. Plastome sequences were determined from high throughput sequencing libraries and the assemblies represent over 28.7 Mbases of sequence data. Phylogenetic signal was characterized in 14 partitions, including (1) complete plastomes; (2) protein coding regions; (3) noncoding regions; and (4) three loci commonly used in single and multi-gene studies of grasses. Each of the four main partitions was further refined, alternatively including or excluding positively selected codons and also the gaps introduced by the alignment. All 76 protein coding plastome loci were found to be predominantly under purifying selection, but specific codons were found to be under positive selection in 65 loci. The loci that have been widely used in multi-gene phylogenetic studies had among the highest proportions of positively selected codons, suggesting caution in the interpretation of these earlier results. Plastome phylogenomic analyses confirmed the backbone topology for Poaceae with maximum bootstrap support (BP). Among the 14 analyses, 82 clades out of 309 resolved were maximally supported in all trees. Analyses of newly sequenced plastomes were in agreement with current classifications. Five of seven partitions in which alignment gaps were removed retrieved Panicoideae as sister to the remaining PACMAD subfamilies. Alternative topologies were recovered in trees from partitions that included alignment gaps. This suggests that ambiguities in aligning these uncertain regions might introduce a false signal. Resolution of these and other critical branch points in the phylogeny of Poaceae will help to better understand the selective forces that drove the radiation of the BOP and PACMAD clades comprising more than 99.9% of grass diversity.

Complete structure and variation of the chloroplast genome of Agropyron cristatum (L.) Gaertn

Agropyron cristatum (L.) Gaertner, a perennial grass in the tribe Triticeae (Poaceae), is a wild relative of cereal crops that is suitable for genetic improvement. In this study, we first sequenced the complete chloroplast (cp) genome of Ag. cristatum using Hiseq4000 PE150. The Ag. cristatum chloroplast genome is 135,554bp in length, has a typical quadripartite structure and contains 76 protein-coding genes, 29 tRNA genes and four rRNA genes. The cp genome of Ag. cristatum was used for comparison with other seven Triticeae species. One large variable region (800bp), which primarily contained the rpl23 (non-reciprocally translocated from IRs) and accD genes, was detected between rbcL gene and psaI gene within LSC region. The deletion of the accD and translocated rpl23 genes in Ag. cristatum indicated an independent gene-loss events or an additional divergence in Triticeae. Analyses of the dn/ds ratio and K2-P's genetic distance for 76 protein-coding genes showed that genes with evolutionary divergence might suffer from the effect of sequence regional constraints or gene functional constraints in Triticeae species. Our research will generally contribute to the knowledge of plastid genome evolution in Triticeae.

The complete chloroplast genome sequence of the CAM epiphyte Spanish moss (Tillandsia usneoides, Bromeliaceae) and its comparative analysis

Spanish moss (Tillandsia usneoides) is an epiphytic bromeliad widely distributed throughout tropical and warm temperate America. This plant is highly adapted to extreme environmental conditions. Striking features of this species include specialized trichomes (scales) covering the surface of its shoots aiding the absorption of water and nutrients directly from the atmosphere and a specific photosynthesis using crassulacean acid metabolism (CAM). Here we report the plastid genome of Spanish moss and present the comparison of genome organization and sequence evolution within Poales. The plastome of Spanish moss has a quadripartite structure consisting of a large single copy (LSC, 87,439 bp), two inverted regions (IRa and IRb, 26,803 bp) and short single copy (SSC, 18,612 bp) region. The plastid genome had 37.2% GC content and 134 genes with 88 being unique protein-coding genes and 20 of these are duplicated in the IR, similar to other reported bromeliads. Our study shows that early diverging lineages of Poales do not have high substitution rates as compared to grasses, and plastid genomes of bromeliads show structural features considered to be ancestral in graminids. These include the loss of the introns in the clpP and rpoC1 genes and the complete loss or partial degradation of accD and ycf genes in the Graminid clade. Further structural rearrangements appeared in the graminids lacking in Spanish moss, which include a 28-kb inversion between the trnG-UCC-rps14 region and 6-kb in the trnG-UCC-psbD, followed by a third <1kb inversion in the trnT sequence.