Plastid phylogenomics of the cool-season grass subfamily: clarification of relationships among early-diverging tribes

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.

The low level of plastome differentiation observed in some lineages of Poales hinders molecular species identification

Chloroplast genomes are a source of information successfully used in various fields of plant genetics, including molecular species identification. However, recent studies indicate an extremely low level of interspecific variability in the plastomes of some taxonomic groups of plants, including the genus Stipa L., which is a representative of the grass family. In this study we aimed to analyze the level of chloroplast genome diversity within particular genera as well as the effectiveness of identifying plant species in the Poaceae family and the other representatives of Poales order. Analysis of complete plastid genome alignments created for 96 genera comprising 793 species and 1707 specimens obtained from the GenBank database allowed defining and categorizing molecular diagnostic characters distinguishing the analyzed species from the other representatives of the genus. The results also demonstrate which species do not have any species-specific mutations, thereby they cannot be identified on the basis of differences between the complete chloroplast genomes. Our research showed a huge diversity of the analyzed species in terms of the number of molecular diagnostic characters and indicated which genera pose a particular challenge in terms of molecular species identification. The results show that a very low level of genetic diversity between plastomes is not uncommon in Poales. This is the first extensive research on super-barcoding that tests this method on a large data set and illustrates its effectiveness against the background of phylogenetic relationships.

Independent recruitment of FRUITFULL-like transcription factors in the convergent origins of vernalization-responsive grass flowering

Flowering in response to low temperatures (vernalization) has evolved multiple times independently across angiosperms as an adaptation to match reproductive development with the short growing season of temperate habitats. Despite the context of a generally conserved flowering time network, evidence suggests that the genes underlying vernalization responsiveness are distinct across major plant clades. Whether different or similar mechanisms underlie vernalization-induced flowering at narrower (e.g., family-level) phylogenetic scales is not well understood. To test the hypothesis that vernalization responsiveness has evolved convergently in temperate species of the grass family (Poaceae), we carried out flowering time experiments with and without vernalization in several representative species from different subfamilies. We then determined the likelihood that vernalization responsiveness evolved through parallel mechanisms by quantifying the response of Pooideae vernalization pathway FRUITFULL (FUL)-like genes to extended periods of cold. Our results demonstrate that vernalization-induced flowering has evolved multiple times independently in at least five grass subfamilies, and that different combinations of FUL-like genes have been recruited to this pathway on several occasions.

Insights into adaptive evolution of plastomes in Stipa L. (Poaceae)

BACKGROUND

The study presents results of research on the evolution of plastid genomes in Stipa L. which is a large genus of the Poaceae family, comprising species diverse in terms of geographic distribution, growing under highly variated habitat conditions. Complete plastome sequences of 43 taxa from Stipeae and Ampelodesmae tribes were analyzed for the variability of the coding regions against the background of phylogenetic relationships within the genus Stipa. The research hypothesis put forward in our research was that some of coding regions are affected by a selection pressure differentiated between individual phylogenetic lines of Stipa, potentially reducing the phylogenetic informativeness of these CDS. The study aimed to answer the question, which genes evolve in Stipa most rapidly and what kind of changes in the properties of encoded amino acids this entails. Another goal of this research was to find out whether individual genes are affected by positive selection and finally, whether selective pressure is uniform within the genus or does it vary between particular evolutionary lines within the genus.

RESULTS

Results of our study proved the presence of selective pressure in 11 genes: ccsA, matK, ndhC, ndhF, ndhK, rbcL, rpoA rpoC1, rpoC2, rps8 and rps11. For the first time the effect of positive selection on the rps8, rps11, and ndhK genes was documented in grasses. The varied pace of evolution, different intensity and effects of selective pressure have been demonstrated between particular phylogenetic lines of the genus tested.

CONCLUSIONS

Positive selection in plastid genome in Stipa mostly affects photosynthetic genes. The potential strongest adaptive pressure was observed in the rbcL gene, especially in the oldest evolutionary group comprising Central Asian high-mountain species: S. basiplumosa, S. klimesii, S. penicillata and S. purpurea, where adaptive pressure probably affected the amino acids directly related to the efficiency of CO₂ assimilation.

A phylogenetic analysis of Bromus (Poaceae: Pooideae: Bromeae) based on nuclear ribosomal and plastid data, with a focus on Bromus sect. Bromus

To investigate phylogenetic relationships among and within major lineages of Bromus, with focus on Bromus sect. Bromus, we analyzed DNA sequences from two nuclear ribosomal (ITS, ETS) and two plastid (rpl32-trnL_UAG , matK) regions. We sampled 103 ingroup accessions representing 26 taxa of B. section Bromus and 15 species of other Bromus sections. Our analyses confirm the monophyly of Bromus s.l. and identify incongruence between nuclear ribosomal and plastid data partitions for relationships within and among major Bromus lineages. Results support classification of B. pumilio and B. gracillimus within B. sect. Boissiera and B. sect. Nevskiella, respectively. These species are sister groups and are closely related to B. densus (B. sect. Mexibromus) in nrDNA trees and Bromus sect. Ceratochloa in plastid trees. Bromus sect. Bromopsis is paraphyletic. In nrDNA trees, species of Bromus sects. Bromopsis, Ceratochloa, Neobromus, and Genea plus B. rechingeri of B. sect. Bromus form a clade, in which B. tomentellus is sister to a B. sect. Genea-B. rechingeri clade. In the plastid trees, by contrast, B. sect. Bromopsis species except B. tomentosus form a clade, and B. tomentosus is sister to a clade comprising B. sect. Bromus and B. sect. Genea species. Affinities of B. gedrosianus, B. pulchellus, and B. rechingeri (members of the B. pectinatus complex), as well as B. oxyodon and B. sewerzowii, are discordant between nrDNA and plastid trees. We infer these species may have obtained their plastomes via chloroplast capture from species of B. sect. Bromus and B. sect. Genea. Within B. sect. Bromus, B. alopecuros subsp. caroli-henrici, a clade comprising B. hordeaceus and B. interruptus, and B. scoparius are successive sister groups to the rest of the section in the nrDNA phylogeny. Most relationships among the remaining species of B. sect. Bromus are unresolved in the nrDNA and plastid trees. Given these results, we infer that most B. sect. Bromus species likely diversified relatively recently. None of the subdivisional taxa proposed for Bromus sect. Bromus over the last century correspond to natural groups identified in our phylogenetic analyses except for a group including B. hordeaceus and B. interruptus.

A well-supported nuclear phylogeny of Poaceae and implications for the evolution of C4 photosynthesis

Poaceae (the grasses) includes rice, maize, wheat, and other crops, and is the most economically important angiosperm family. Poaceae is also one of the largest plant families, consisting of over 11 000 species with a global distribution that contributes to diverse ecosystems. Poaceae species are classified into 12 subfamilies, with generally strong phylogenetic support for their monophyly. However, many relationships within subfamilies, among tribes and/or subtribes, remain uncertain. To better resolve the Poaceae phylogeny, we generated 342 transcriptomic and seven genomic datasets; these were combined with other genomic and transcriptomic datasets to provide sequences for 357 Poaceae species in 231 genera, representing 45 tribes and all 12 subfamilies. Over 1200 low-copy nuclear genes were retrieved from these datasets, with several subsets obtained using additional criteria, and used for coalescent analyses to reconstruct a Poaceae phylogeny. Our results strongly support the monophyly of 11 subfamilies; however, the subfamily Puelioideae was separated into two non-sister clades, one for each of the two previously defined tribes, supporting a hypothesis that places each tribe in a separate subfamily. Molecular clock analyses estimated the crown age of Poaceae to be ∼101 million years old. Ancestral character reconstruction of C₃/C₄ photosynthesis supports the hypothesis of multiple independent origins of C₄ photosynthesis. These origins are further supported by phylogenetic analysis of the ppc gene family that encodes the phosphoenolpyruvate carboxylase, which suggests that members of three paralogous subclades (ppc-aL1a, ppc-aL1b, and ppc-B2) were recruited as functional C₄ppc genes. This study provides valuable resources and a robust phylogenetic framework for evolutionary analyses of the grass family.

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.

Unprecedented organelle genomic variations in morning glories reveal independent evolutionary scenarios of parasitic plants and the diversification of plant mitochondrial complexes

Background

The morning glories (Convolvulaceae) are distributed worldwide and produce economically important crops, medicinal herbs, and ornamentals. Members of this family are diverse in morphological characteristics and trophic modes, including the leafless parasitic Cuscuta (dodders). Organelle genomes were generally used for studying plant phylogeny and genomic variations. Notably, plastomes in parasitic plants always show non-canonical features, such as reduced size and accelerated rates. However, few organelle genomes of this group have been sequenced, hindering our understanding of their evolution, and dodder mitogenome in particular.

Results

We assembled 22 new mitogenomes and 12 new plastomes in Convolvulaceae. Alongside previously known ones, we totally analyzed organelle genomes of 23 species in the family. Our sampling includes 16 leafy autotrophic species and 7 leafless parasitic dodders, covering 8 of the 12 tribes. Both the plastid and mitochondrial genomes of these plants have encountered variations that were rarely observed in other angiosperms. All of the plastomes possessed atypical IR boundaries. Besides the gene and IR losses in dodders, some leafy species also showed gene and intron losses, duplications, structural variations, and insertions of foreign DNAs. The phylogeny reconstructed by plastid protein coding sequences confirmed the previous relationship of the tribes. However, the monophyly of ‘Merremieae’ and the sister group of Cuscuta remained uncertain. The mitogenome was significantly inflated in Cuscuta japonica, which has exceeded over 800 kb and integrated massive DNAs from other species. In other dodders, mitogenomes were maintained in small size, revealing divergent evolutionary strategies. Mutations unique to plants were detected in the mitochondrial gene ccmFc, which has broken into three fragments through gene fission and splicing shift. The unusual changes likely initially happened to the common ancestor of the family and were caused by a foreign insertion from rosids followed by double-strand breaks and imprecise DNA repairs. The coding regions of ccmFc expanded at both sides after the fission, which may have altered the protein structure.

Conclusions

Our family-scale analyses uncovered unusual scenarios for both organelle genomes in Convolvulaceae, especially in parasitic plants. The data provided valuable genetic resources for studying the evolution of Convolvulaceae and plant parasitism.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12915-022-01250-1.

Phylotranscriptomics Resolves the Phylogeny of Pooideae and Uncovers Factors for Their Adaptive Evolution

Adaptation to cool climates has occurred several times in different angiosperm groups. Among them, Pooideae, the largest grass subfamily with ∼3,900 species including wheat and barley, have successfully occupied many temperate regions and play a prominent role in temperate ecosystems. To investigate possible factors contributing to Pooideae adaptive evolution to cooling climates, we performed phylogenetic reconstruction using five gene sets (with 1,234 nuclear genes and their subsets) from 157 transcriptomes/genomes representing all 15 tribes and 24 of 26 subtribes. Our phylogeny supports the monophyly of all tribes (except Diarrheneae) and all subtribes with at least two species, with strongly supported resolution of their relationships. Molecular dating suggests that Pooideae originated in the late Cretaceous, with subsequent divergences under cooling conditions first among many tribes from the early middle to late Eocene and again among genera in the middle Miocene and later periods. We identified a cluster of gene duplications (CGD5) shared by the core Pooideae (with 80% Pooideae species) near the Eocene–Oligocene transition, coinciding with the transition from closed to open habitat and an upshift of diversification rate. Molecular evolutionary analyses homologs of CBF for cold resistance uncovered tandem duplications during the core Pooideae history, dramatically increasing their copy number and possibly promoting adaptation to cold habitats. Moreover, duplication of AP1/FUL-like genes before the Pooideae origin might have facilitated the regulation of the vernalization pathway under cold environments. These and other results provide new insights into factors that likely have contributed to the successful adaptation of Pooideae members to temperate regions.

The chloroplast genomes of four Bupleurum (Apiaceae) species endemic to Southwestern China, a diversity center of the genus, as well as their evolutionary implications and phylogenetic inferences

Background

As one of the largest genera in Apiaceae, Bupleurum L. is well known for its high medicinal value. The genus has frequently attracted the attention of evolutionary biologist and taxonomist for its distinctive characteristics in the Apiaceae family. Although some chloroplast genomes data have been now available, the changes in the structure of chloroplast genomes and selective pressure in the genus have not been fully understood. In addition, few of the species are endemic to Southwest China, a distribution and diversity center of Chinese Bupleurum. Endemic species are key components of biodiversity and ecosystems, and investigation of the chloroplast genomes features of endemic species in Bupleurum will be helpful to develop a better understanding of evolutionary process and phylogeny of the genus. In this study, we analyzed the sequences of whole chloroplast genomes of 4 Southwest China endemic Bupleurum species in comparison with the published data of 17 Bupleurum species to determine the evolutionary characteristics of the genus and the phylogenetic relationships of Asian Bupleurum.

Results

The complete chloroplast genome sequences of the 4 endemic Bupleurum species are 155,025 bp to 155,323 bp in length including a SSC and a LSC region separated by a pair of IRs. Comparative analysis revealed an identical chloroplast gene content across the 21 Bupleurum species, including a total of 114 unique genes (30 tRNA genes, 4 rRNA genes and 80 protein-coding genes). Chloroplast genomes of the 21 Bupleurum species showed no rearrangements and a high sequence identity (96.4–99.2%). They also shared a similar tendency of SDRs and SSRs, but differed in number (59–83). In spite of their high conservation, they contained some mutational hotspots, which can be potentially exploited as high-resolution DNA barcodes for species discrimination. Selective pressure analysis showed that four genes were under positive selection. Phylogenetic analysis revealed that the 21 Bupleurum formed two major clades, which are likely to correspond to their geographical distribution.

Conclusions

The chloroplast genome data of the four endemic Bupleurum species provide important insights into the characteristics and evolution of chloroplast genomes of this genu, and the phylogeny of Bupleurum.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-021-08008-z.

Variation in climatic tolerance, but not stomatal traits, partially explains Pooideae grass species distributions

BACKGROUND AND AIMS

Grasses in subfamily Pooideae live in some of the world's harshest terrestrial environments, from frigid boreal zones to the arid windswept steppe. It is hypothesized that the climate distribution of species within this group is driven by differences in climatic tolerance, and that tolerance can be partially explained by variation in stomatal traits.

METHODS

We determined the aridity index (AI) and minimum temperature of the coldest month (MTCM) for 22 diverse Pooideae accessions and one outgroup, and used comparative methods to assess predicted relationships for climate traits versus fitness traits, stomatal diffusive conductance to water (gw) and speed of stomatal closure following drought and/or cold.

KEY RESULTS

Results demonstrate that AI and MTCM predict variation in survival/regreening following drought/cold, and gw under drought/cold is positively correlated with δ 13C-measured water use efficiency (WUE). However, the relationship between climate traits and fitness under drought/cold was not explained by gw or speed of stomatal closure.

CONCLUSIONS

These findings suggest that Pooideae distributions are at least partly determined by tolerance to aridity and above-freezing cold, but that variation in tolerance is not uniformly explained by variation in stomatal traits.

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.

A 313 plastome phylogenomic analysis of Pooideae: Exploring relationships among the largest subfamily of grasses

In this study, we analyzed 313 plastid genomes (plastomes) of Poaceae with a focus on expanding our current knowledge of relationships among the subfamily Pooideae, which represented over half the dataset (164 representatives). In total, 47 plastomes were sequenced and assembled for this study. This is the largest study of its kind to include plastome-level data, to not only increase sampling at both the taxonomic and molecular levels with the aim of resolving complex and reticulate relationships, but also to analyze the effects of alignment gaps in large-scale analyses, as well as explore divergences in the subfamily with an expanded set of 14 accepted grass fossils for more accurate calibrations and dating. Incorporating broad systematic assessments of Pooideae taxa conducted by authors within the last five years, we produced a robust phylogenomic reconstruction for the subfamily, which included all but two supergeneric taxa (Calothecinae and Duthieeae). We further explored how including alignment gaps in plastome analyses oftentimes can produce incorrect or misinterpretations of complex or reticulate relationships among taxa of Pooideae. This presented itself as consistently changing relationships at specific nodes for different stripping thresholds (percentage-based removal of gaps per alignment column). Our summary recommendation for large-scale genomic plastome datasets is to strip alignment columns of all gaps to increase pairwise identity and reduce errant signal from poly A/T bias. To do this we used the "mask alignment" tool in Geneious software. Finally, we determined an overall divergence age for Pooideae of roughly 84.8 Mya, which is in line with, but slightly older than most recent estimates.

Molecular Phylogenetics and Micromorphology of Australasian Stipeae (Poaceae, Subfamily Pooideae), and the Interrelation of Whole-Genome Duplication and Evolutionary Radiations in This Grass Tribe

The mainly Australian grass genus Austrostipa (tribe Stipeae) comprising approximately 64 species represents a remarkable example of an evolutionary radiation. To investigate aspects of diversification, macro- and micromorphological variation in this genus, we conducted molecular phylogenetic and scanning electron microscopy (SEM) analyses including representatives from most of Austrostipa's currently accepted subgenera. Because of its taxonomic significance in Stipeae, we studied the lemma epidermal pattern (LEP) in 34 representatives of Austrostipa. Plastid DNA variation within Austrostipa was low and only few lineages were resolved. Nuclear ITS and Acc1 yielded comparable groupings of taxa and resolved subgenera Arbuscula, Petaurista, and Bambusina in a common clade and as monophyletic. In most of the Austrostipa species studied, the LEP was relatively uniform (typical maize-like), but six species had a modified cellular structure. The species representing subgenera Lobatae, Petaurista, Bambusina as well as A. muelleri from subg. Tuberculatae were well-separated from all the other species included in the analysis. We suggest recognizing nine subgenera in Austrostipa (with number of species): Arbuscula (4), Aulax (2), Austrostipa (36), Bambusina (2), Falcatae (10), Lobatae (5), Longiaristatae (2), Petaurista (2) and the new subgenus Paucispiculatae (1) encompassing A. muelleri. Two paralogous sequence copies of Acc1, forming two distinct clades, were found in polyploid Austrostipa and Anemanthele. We found analogous patterns for our samples of Stipa s.str. with their Acc1 clades strongly separated from those of Austrostipa and Anemanthele. This underlines a previous hypothesis of Tzvelev (1977) that most extant Stipeae are of hybrid origin. We also prepared an up-to-date survey and reviewed the chromosome number variation for our molecularly studied taxa and the whole tribe Stipeae. The chromosome base number patterns as well as dysploidy and whole-genome duplication events were interpreted in a phylogenetic framework. The rather coherent picture of chromosome number variation underlines the enormous phylogenetic and evolutionary significance of this frequently ignored character.

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.

Checklist of the grasses of India

A checklist of the grasses of India is presented, as compiled from survey of all available literature. Of the twelve subfamilies of grasses, ten are represented in India. Most subfamilies have been examined by taxonomic experts for up-to-date nomenclature. The list includes 1506 species plus infraspecific taxa and presents information on types, synonyms, distribution within India, and habit. Twelve new combinations are made, viz. Arctopoa tibetica (Munro ex Stapf) Prob. var. aristulata (Stapf) E.A. Kellogg, comb. nov.; Chimonocalamus nagalandianus (H.B. Naithani) L.G. Clark, comb. nov.; Chionachne digitata (L.f.) E.A. Kellogg, comb. nov.; Chionachne wallichiana (Nees) E.A. Kellogg, comb. nov.; Dinebra polystachyos (R. Br.) E.A. Kellogg, comb. nov.; Moorochloa eruciformis (Sm.) Veldkamp var. divaricata (Basappa & Muniv.) E.A. Kellogg, comb. nov.; Phyllostachys nigra (Lodd. ex Lindl.) Munro var. puberula (Miq.) Kailash, comb. & stat. nov.; Tzveleviochloa schmidii (Hook. f.) E.A. Kellogg, comb. nov.; Urochloa lata (Schumach.) C.E. Hubb. var. pubescens (C.E. Hubb.) E.A. Kellogg, comb. nov.; Urochloa ramosa (L.) T.Q. Nguyen var. pubescens (Basappa & Muniy.) E.A. Kellogg, comb. nov.; Urochloa semiundulata (Hochst. ex A. Rich.) Ashalatha & V.J. Nair var. intermedia (Basappa & Muniy.) E.A. Kellogg, comb. nov.

Plastome phylogenomics and phylogenetic diversity of endangered and threatened grassland species (Poaceae) in a North American tallgrass prairie

Native grasslands are one of the most endangered ecosystems in North America. In this study, we examined the ecological and evolutionary roles of endangered and threatened (e/t) grasses by establishing robust evolutionary relationships with other nonthreatened native and introduced grass species of the community. We hypothesized that the phylogenomic distribution of e/t species of grasses in Illinois would be phylogenetically clustered because closely related species would be vulnerable to the same threats and have similar requirements for survival. This study presents the first time a phylogeny based on complete plastome DNA of Poaceae was analyzed by phylogenetic diversity analysis. To avoid the disturbance of e/t populations, DNA was extracted from herbarium specimens. Next-generation sequencing (NGS) techniques were used to sequence DNA of plastid genomes (plastomes). The resulting phylogenomic tree was analyzed by phylogenetic diversity metrics. The extracted DNA successfully produced complete plastomes demonstrating that herbarium material is a practical source of DNA for genomic studies. The phylogenomic tree was strongly supported and defined Dichanthelium as a separate clade from Panicum. The phylogenetic metrics revealed phylogenetic clustering of e/t species, confirming our hypothesis.

Evolutionary patterns of plastome uncover diploid-polyploid maternal relationships in Triticeae

To investigate the diploid-polyploid relationships and the role of maternal progenitors in establishment of polyploid richness in Triticeae, 35 polyploids representing almost all genomic constitutions together with 48 diploid taxa representing 20 basic genomes in the tribe were analyzed. Phylogenomic reconstruction, genetic distance matrix, and nucleotide diversity patterns of plastome sequences indicated that (1) The maternal donor of the annual polyploid species with the U- and D-genome are related to extant Ae. umbellulata and Ae. tauschii, respectively. The maternal donor to the annual polyploid species with the S-, G-, and B-genome originated from the species of Sitopsis section of the genus Aegilops. The annual species with the Xe-containing polyploids were donated by Eremopyrum as the female parent; (2) Pseudoroegneria and Psathyrostachys were the maternal donor of perennial species with the St- and Ns-containing polyploids, respectively; (3) The Lophopyrum, Thinopyrum and Dasypyrum genomes contributed cytoplasm genome to Pseudoroegneria species as a result of incomplete lineage sorting and/or chloroplast captures, and these lineages were genetically transmitted to the St-containing polyploid species via polyploidization; (4) There is a reticulate relationship among the St-containing polyploid species. It can be suggested that genetic heterogeneity might associate with the richness of the polyploids in Triticeae.

Complete chloroplast genomes of all six Hosta species occurring in Korea: molecular structures, comparative, and phylogenetic analyses

BACKGROUND

The genus Hosta is a group of economically appreciated perennial herbs consisting of approximately 25 species that is endemic to eastern Asia. Due to considerable morphological variability, the genus has been well recognized as a group with taxonomic problems. Chloroplast is a cytoplasmic organelle with its own genome, which is the most commonly used for phylogenetic and genetic diversity analyses for land plants. To understand the genomic architecture of Hosta chloroplasts and examine the level of nucleotide and size variation, we newly sequenced four (H. clausa, H. jonesii, H. minor, and H. venusta) and analyzed six Hosta species (including the four, H. capitata and H. yingeri) distributed throughout South Korea.

RESULTS

The average size of complete chloroplast genomes for the Hosta taxa was 156,642 bp with a maximum size difference of ~ 300 bp. The overall gene content and organization across the six Hosta were nearly identical with a few exceptions. There was a single tRNA gene deletion in H. jonesii and four genes were pseudogenized in three taxa (H. capitata, H. minor, and H. jonesii). We did not find major structural variation, but there were a minor expansion and contractions in IR region for three species (H. capitata, H. minor, and H. venusta). Sequence variations were higher in non-coding regions than in coding regions. Four genic and intergenic regions including two coding genes (psbA and ndhD) exhibited the largest sequence divergence showing potential as phylogenetic markers. We found compositional codon usage bias toward A/T at the third position. The Hosta plastomes had a comparable number of dispersed and tandem repeats (simple sequence repeats) to the ones identified in other angiosperm taxa. The phylogeny of 20 Agavoideae (Asparagaceae) taxa including the six Hosta species inferred from complete plastome data showed well resolved monophyletic clades for closely related taxa with high node supports.

CONCLUSIONS

Our study provides detailed information on the chloroplast genome of the Hosta taxa. We identified nucleotide diversity hotspots and characterized types of repeats, which can be used for developing molecular markers applicable in various research area.

Structural variation of the complete chloroplast genome and plastid phylogenomics of the genus Asteropyrum (Ranunculaceae)

Two complete chloroplast genome sequences of Asteropyrum, as well as those of 25 other species from Ranunculaceae, were assembled using both Illumina and Sanger sequencing methods to address the structural variation of the cp genome and the controversial systematic position of the genus. Synteny and plastome structure were compared across the family. The cp genomes of the only two subspecies of Asteropyrum were found to be differentiated with marked sequence variation and different inverted repeat-single copy (IR-SC) borders. The plastomes of both subspecies contains 112 genes. However, the IR region of subspecies peltatum carries 27 genes, whereas that of subspecies cavaleriei has only 25 genes. Gene inversions, transpositions, and IR expansion-contraction were very commonly detected in Ranunculaceae. The plastome of Asteropyrum has the longest IR regions in the family, but has no gene inversions or transpositions. Non-coding regions of the cp genome were not ideal markers for inferring the generic relationships of the family, but they may be applied to interpret species relationship within the genus. Plastid phylogenomic analysis using complete cp genome with Bayesian method and partitioned modeling obtained a fully resolved phylogenetic framework for Ranunculaceae. Asteropyrum was detected to be sister to Caltha, and diverged early from subfamily Ranunculoideae.

The complete chloroplast genome sequence of Melica scabrosa (Poaceae)

Melica scabrosa Trin. is an important forage grass of Poaceae, wildly distributed in the Northeast Asia to Qinghai-Xizang Plateau. The complete chloroplast genome sequence of M. scabrosa was obtained by de novo assembly using whole genome sequence data. The chloroplast genome is 134,889 bp in length, containing 80,560 bp in a large single copy (LSC), 12,706 bp in a small single copy (SSC) and 20,810 bp in a pair of inverted repeats (IRs). A total of 129 genes including 83 protein-coding genes and 38 structural RNA genes were identified. Phylogenetic analysis represented close relationship among Melica species. This chloroplast genome sequencing offers a useful resource for future genetics and phylogenetic studies.

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.

Transcriptome Analysis and Differential Expression in Tall Fescue Harboring Different Endophyte Strains in Response to Water Deficit

Two tall fescue [Lolium arundinaceum (Schreb.) Darbysh. = Schedonorus arundinaceus (Schreb.) Dumort. = Festuca arundinacea var. arundinacea Schreb.] plant genotypes with an Epichloë coenophiala (Morgan-Jones & W. Gams) C.W. Bacon & Schardl common toxic endophyte (CTE), one with a nontoxic strain (NTE19) and one with another Epichloë species (FaTG-4) were evaluated and compared with their respective endophyte-free clones for responses to water-deficit stress in the greenhouse. One of the plant genotypes (P27) showed a positive effect of its CTE strain on tiller production after stress and resumed watering. In transcriptome analysis of the pseudostems (leaf sheath whorls), differentially expressed genes (DEGs) were defined as having at least twofold expression difference and false discovery rate (FDR) < 0.05 in comparisons of water treatment (stressed or watered), endophyte presence or absence, or both. Stress affected 38% of the plant transcripts including those for the expected stress-response pathways. The DEGs affected by endophyte in stressed plants were unique to individual plant genotypes. In unstressed plants, endophyte presence tended to reduce expression of genes putatively for defense against fungi, but in unstressed P27 endophyte presence there was enhanced expression of dehydrin and heat shock protein genes. Our results indicated subtle and variable effects of endophytes on tall fescue gene expression; where the endophyte confers protection, its effects on plant gene expression may help prime the plant for stress resistance.

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.

Evolution of six novel ORFs in the plastome of Mankyua chejuense and phylogeny of eusporangiate ferns

In this paper, three plastomes of Mankyua chejuense, Helminthostachys zeylanica, and Botrychium ternatum in Ophioglossaceae were completely sequenced in order to investigate the plastome evolution and phylogeny of eusporangiate ferns. They were similar to each other in terms of length and the gene orders; however, six unknown open reading frames (ORFs) were found between rps4 and trnL-UAA genes in M. chejuense. Similar sequence regions of six ORFs of M. chejuense were found at the plastomes of Ophioglossum californicum and H. zeylanica, as well as the mitochondrial genome (mitogenome) of H. zeylanica, but not in B. ternatum. Interestingly, the translated amino acid sequences of three ORFs were more similar to the proteins of distantly related taxa such as algae and bacteria than they were to proteins in land plants. It is likely that the six ORFs region arose from endosymbiotic gene transfer (EGT) or horizontal gene transfer (HGT), but further study is needed to verify this. Phylogenetic analyses suggested that Mankyua was resolved as the earliest diverging lineage and that Ophioglossum was subsequently diverged in Ophioglossaceae. This result supports why the plastome of M. chejuense have contained the most ancestral six ORFs in the family.

Monocot plastid phylogenomics, timeline, net rates of species diversification, the power of multi-gene analyses, and a functional model for the origin of monocots

PREMISE OF THE STUDY

We present the first plastome phylogeny encompassing all 77 monocot families, estimate branch support, and infer monocot-wide divergence times and rates of species diversification.

METHODS

We conducted maximum likelihood analyses of phylogeny and BAMM studies of diversification rates based on 77 plastid genes across 545 monocots and 22 outgroups. We quantified how branch support and ascertainment vary with gene number, branch length, and branch depth.

KEY RESULTS

Phylogenomic analyses shift the placement of 16 families in relation to earlier studies based on four plastid genes, add seven families, date the divergence between monocots and eudicots+Ceratophyllum at 136 Mya, successfully place all mycoheterotrophic taxa examined, and support recognizing Taccaceae and Thismiaceae as separate families and Arecales and Dasypogonales as separate orders. Only 45% of interfamilial divergences occurred after the Cretaceous. Net species diversification underwent four large-scale accelerations in PACMAD-BOP Poaceae, Asparagales sister to Doryanthaceae, Orchidoideae-Epidendroideae, and Araceae sister to Lemnoideae, each associated with specific ecological/morphological shifts. Branch ascertainment and support across monocots increase with gene number and branch length, and decrease with relative branch depth. Analysis of entire plastomes in Zingiberales quantifies the importance of non-coding regions in identifying and supporting short, deep branches.

CONCLUSIONS

We provide the first resolved, well-supported monocot phylogeny and timeline spanning all families, and quantify the significant contribution of plastome-scale data to resolving short, deep branches. We outline a new functional model for the evolution of monocots and their diagnostic morphological traits from submersed aquatic ancestors, supported by convergent evolution of many of these traits in aquatic Hydatellaceae (Nymphaeales).

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.

Investigation of mitochondrial-derived plastome sequences in the Paspalum lineage (Panicoideae; Poaceae)

BACKGROUND

The grass family (Poaceae), ca. 12,075 species, is a focal point of many recent studies that aim to use complete plastomes to reveal and strengthen relationships within the family. The use of Next Generation Sequencing technology has revealed intricate details in many Poaceae plastomes; specifically the trnI - trnL intergenic spacer region. This study investigates this region and the putative mitochondrial inserts within it in complete plastomes of Paspalum and other Poaceae.

RESULTS

Nine newly sequenced plastomes, seven of which contain an insert within the trnI - trnL intergenic spacer, were combined into plastome phylogenomic and divergence date analyses with 52 other species. A robust Paspalum topology was recovered, originating at 10.6 Ma, with the insert arising at 8.7 Ma. The alignment of the insert across Paspalum reveals 21 subregions with pairwise homology in 19. In an analysis of emergent self-organizing maps of tetranucleotide frequencies, the Paspalum insert grouped with mitochondrial DNA.

CONCLUSIONS

A hypothetical ancestral insert, 17,685 bp in size, was found in the trnI - trnL intergenic spacer for the Paspalum lineage. A different insert, 2808 bp, was found in the same region for Paraneurachne muelleri. Seven different intrastrand deletion events were found within the Paspalum lineage, suggesting selective pressures to remove large portions of noncoding DNA. Finally, a tetranucleotide frequency analysis was used to determine that the origin of the insert in the Paspalum lineage is mitochondrial DNA.

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.

Complete chloroplast genome of seven Fritillaria species, variable DNA markers identification and phylogenetic relationships within the genus

Fritillaria spp. constitute important traditional Chinese medicinal plants. Xinjiang is one of two diversity hotspots in China in which eight Fritillaria species occur, two of which are endemic to the region. Furthermore, the phylogenetic relationships of Xinjiang Fritillaria species (including F. yuminensis) within the genus are unclear. In the present study, we sequenced the chloroplast (cp) genomes of seven Fritillaria species in Xinjiang using the Illumina HiSeq platform, with the aim of assessing the global structural patterns of the seven cp genomes and identifying highly variable cp DNA sequences. These were compared to previously sequenced Fritillaria cp genomes. Phylogenetic analysis was then used to evaluate the relationships of the Xinjiang species and assess the evolution of an undivided stigma. The seven cp genomes ranged from 151,764 to 152,112 bp, presenting a traditional quadripartite structure. The gene order and gene content of the seven cp genomes were identical. A comparison of the 13 cp genomes indicated that the structure is highly conserved. Ten highly divergent regions were identified that could be valuable in phylogenetic and population genetic studies. The phylogenetic relationships of the 13 Fritillaria species inferred from the protein-coding genes, large single-copy, small single-copy, and inverted repeat regions were identical and highly resolved. The phylogenetic relationships of the species corresponded with their geographic distribution patterns, in that the north group (consisting of eight species from Xinjiang and Heilongjiang in North China) and the south group (including six species from South China) were basically divided at 40°N. Species with an undivided stigma were not monophyletic, suggesting that this trait might have evolved several times in the genus.

The complete chloroplast genome sequence of wild oat, Avena sterilis L. (Poaceae) and its phylogeny

Wild oat, Avena sterilis L. is a stout broad-leaved annual grass resembling cultivated oats in general appearance. In this study, we sequenced the complete chloroplast (cp) genome sequence of A. sterilis for the first time to investigate their phylogenetic relationship in the family Poaceae. The complete cp genome sequence is 135,887 bp in length with 38.5% overall GC content and exhibits a typical quadripartite structure comprising one pair of inverted repeats (21,603 bp) separated by a small single-copy region (12,575 bp) and a large single-copy region (80,106). The cp genome encodes 111 unique genes, 76 of which are protein-coding genes, four rRNA genes, 30 tRNA genes, and 18 duplicated genes in the inverted repeat region. The phylogenetic analysis indicated A. sterilis closely clustered with the cultivated oat, A. sativa L.

Brassicales phylogeny inferred from 72 plastid genes: A reanalysis of the phylogenetic localization of two paleopolyploid events and origin of novel chemical defenses

PREMISE OF THE STUDY

Previous phylogenetic studies employing molecular markers have yielded various insights into the evolutionary history across Brassicales, but many relationships between families remain poorly supported or unresolved. A recent phylotranscriptomic approach utilizing 1155 nuclear markers obtained robust estimates for relationships among 14 of 17 families. Here we report a complete family-level phylogeny estimated using the plastid genome.

METHODS

We conducted phylogenetic analyses on a concatenated data set comprising 44,926 bp from 72 plastid genes for species distributed across all 17 families. Our analysis includes three additional families, Tovariaceae, Salvadoraceae, and Setchellanthaceae, that were omitted in the previous phylotranscriptomic study.

KEY RESULTS

Our phylogenetic analyses obtained fully resolved and strongly supported estimates for all nodes across Brassicales. Importantly, these findings are congruent with the topology reported in the phylotranscriptomic study. This consistency suggests that future studies could utilize plastid genomes as markers for resolving relationships within some notoriously difficult clades across Brassicales. We used this new phylogenetic framework to verify the placement of the At-α event near the origin of Brassicaceae, with median date estimates of 31.8 to 42.8 million years ago and restrict the At-β event to one of two nodes with median date estimates between 85 to 92.2 million years ago. These events ultimately gave rise to novel chemical defenses and are associated with subsequent shifts in net diversification rates.

CONCLUSIONS

We anticipate that these findings will aid future comparative evolutionary studies across Brassicales, including selecting candidates for whole-genome sequencing projects.

Karyotype evolution in Phalaris (Poaceae): The role of reductional dysploidy, polyploidy and chromosome alteration in a wide-spread and diverse genus

Karyotype characteristics can provide valuable information on genome evolution and speciation, in particular in taxa with varying basic chromosome numbers and ploidy levels. Due to its worldwide distribution, remarkable variability in morphological traits and the fact that ploidy change plays a key role in its evolution, the canary grass genus Phalaris (Poaceae) is an excellent study system to investigate the role of chromosomal changes in species diversification and expansion. Phalaris comprises diploid species with two basic chromosome numbers of x = 6 and 7 as well as polyploids based on x = 7. To identify distinct karyotype structures and to trace chromosome evolution within the genus, we apply fluorescence in situ hybridisation (FISH) of 5S and 45S rDNA probes in four diploid and four tetraploid Phalaris species of both basic numbers. The data agree with a dysploid reduction from x = 7 to x = 6 as the result of reciprocal translocations between three chromosomes of an ancestor with a diploid chromosome complement of 2n = 14. We recognize three different genomes in the genus: (1) the exclusively Mediterranean genome A based on x = 6, (2) the cosmopolitan genome B based on x = 7 and (3) a genome C based on x = 7 and with a distribution in the Mediterranean and the Middle East. Both auto- and allopolyploidy of genomes B and C are suggested for the formation of tetraploids. The chromosomal divergence observed in Phalaris can be explained by the occurrence of dysploidy, the emergence of three different genomes, and the chromosome rearrangements accompanied by karyotype change and polyploidization. Mapping the recognized karyotypes on the existing phylogenetic tree suggests that genomes A and C are restricted to sections Phalaris and Bulbophalaris, respectively, while genome B occurs across all taxa with x = 7.

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.

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.

Plastome Sequencing of Ten Nonmodel Crop Species Uncovers a Large Insertion of Mitochondrial DNA in Cashew

In plant evolution, intracellular gene transfer (IGT) is a prevalent, ongoing process. While nuclear and mitochondrial genomes are known to integrate foreign DNA via IGT and horizontal gene transfer (HGT), plastid genomes (plastomes) have resisted foreign DNA incorporation and only recently has IGT been uncovered in the plastomes of a few land plants. In this study, we completed plastome sequences for l0 crop species and describe a number of structural features including variation in gene and intron content, inversions, and expansion and contraction of the inverted repeat (IR). We identified a putative in cinnamon ( J. Presl) and other sequenced Lauraceae and an apparent functional transfer of to the nucleus of quinoa ( Willd.). In the orchard tree cashew ( L.), we report the insertion of an ∼6.7-kb fragment of mitochondrial DNA into the plastome IR. BLASTn analyses returned high identity hits to mitogenome sequences including an intact open reading frame. Using three plastome markers for five species of , we generated a phylogeny to investigate the distribution and timing of the insertion. Four species share the insertion, suggesting that this event occurred <20 million yr ago in a single clade in the genus. Our study extends the observation of mitochondrial to plastome IGT to include long-lived tree species. While previous studies have suggested possible mechanisms facilitating IGT to the plastome, more examples of this phenomenon, along with more complete mitogenome sequences, will be required before a common, or variable, mechanism can be elucidated.

Peterson, Soreng RJ, Paszko B. Molecular phylogenetics of cool-season grasses in the subtribes Agrostidinae, Anthoxanthinae, Aveninae, Brizinae, Calothecinae, Koeleriinae and Phalaridinae (Poaceae, Pooideae, Poeae, Poeae chloroplast group 1)

Circumscriptions of and relationships among many genera and suprageneric taxa of the diverse grass tribe Poeae remain controversial. In an attempt to clarify these, we conducted phylogenetic analyses of >2400 new DNA sequences from two nuclear ribosomal regions (ITS, including internal transcribed spacers 1 and 2 and the 5.8S gene, and the 3'-end of the external transcribed spacer (ETS)) and five plastid regions (matK, trnL-trnF, atpF-atpH, psbK-psbI, psbA-rps19-trnH), and of more than 1000 new and previously published ITS sequences, focused particularly on Poeae chloroplast group 1 and including broad and increased species sampling compared to previous studies. Deep branches in the combined plastid and combined ITS+ETS trees are generally well resolved, the trees are congruent in most aspects, branch support across the trees is stronger than in trees based on only ITS and fewer plastid regions, and there is evidence of conflict between data partitions in some taxa. In plastid trees, a strongly supported clade corresponds to Poeae chloroplast group 1 and includes Agrostidinae p.p., Anthoxanthinae, Aveninae s.str., Brizinae, Koeleriinae (sometimes included in Aveninae s.l.), Phalaridinae and Torreyochloinae. In the ITS+ETS tree, a supported clade includes these same tribes as well as Sesleriinae and Scolochloinae. Aveninae s.str. and Sesleriinae are sister taxa and form a clade with Koeleriinae in the ITS+ETS tree whereas Aveninae s.str. and Koeleriinae form a clade and Sesleriinae is part of Poeae chloroplast group 2 in the plastid tree. All species of Trisetum are part of Koeleriinae, but the genus is polyphyletic. Koeleriinae is divided into two major subclades: one comprises Avellinia, Gaudinia, Koeleria, Rostraria, Trisetaria and Trisetum subg. Trisetum, and the other Calamagrostis/Deyeuxia p.p. (multiple species from Mexico to South America), Peyritschia, Leptophyllochloa, Sphenopholis, Trisetopsis and Trisetum subg. Deschampsioidea. Graphephorum, Trisetum cernuum, T. irazuense and T. macbridei fall in different clades of Koeleriinae in plastid vs. nuclear ribosomal trees, and are likely of hybrid origin. ITS and matK trees identify a third lineage of Koeleriinae corresponding to Trisetum subsect. Sibirica, and affinities of Lagurus ovatus with respect to Aveninae s.str. and Koeleriinae are incongruent in nuclear ribosomal and plastid trees, supporting recognition of Lagurus in its own subtribe. A large clade comprises taxa of Agrostidinae, Brizinae and Calothecinae, but neither Agrostidinae nor Calothecinae are monophyletic as currently circumscribed and affinities of Brizinae differ in plastid and nuclear ribosomal trees. Within this clade, one newly identified lineage comprises Calamagrostis coarctata, Dichelachne, Echinopogon (Agrostidinae p.p.) and Relchela (Calothecinae p.p.), and another comprises Chascolytrum (Calothecinae p.p.) and Deyeuxia effusa (Agrostidinae p.p.). Within Agrostidinae p.p., the type species of Deyeuxia and Calamagrostis s.str. are closely related, supporting classification of Deyeuxia as a synonym of Calamagrostis s.str. Furthermore, the two species of Ammophila are not sister taxa and are nested among different groups of Calamagrostis s.str., supporting their classification in Calamagrostis. Agrostis, Lachnagrostis and Polypogon form a clade and species of each are variously intermixed in plastid and nuclear ribosomal trees. Additionally, all but one species from South America classified in Deyeuxia sect. Stylagrostis resolve in Holcinae p.p. (Deschampsia). The current phylogenetic results support recognition of the latter species in Deschampsia, and we also demonstrate Scribneria is part of this clade. Moreover, Holcinae is not monophyletic in its current circumscription because Deschampsia does not form a clade with Holcus and Vahlodea, which are sister taxa. The results support recognition of Deschampsia in its own subtribe Aristaveninae. Substantial further changes to the classification of these grasses will be needed to produce generic circumscriptions consistent with phylogenetic evidence. The following 15 new combinations are made: Calamagrostis × calammophila, C. breviligulata, C. breviligulata subsp. champlainensis, C. × don-hensonii, Deschampsia aurea, D. bolanderi, D. chrysantha, D. chrysantha var. phalaroides, D. eminens, D. eminens var. fulva, D. eminens var. inclusa, D. hackelii, D. ovata, and D. ovata var. nivalis. D. podophora; the new name Deschampsia parodiana is proposed; the new subtribe Lagurinae is described; and a second-step lectotype is designated for the name Deyeuxia phalaroides.

DNA barcoding of invasive plants in China: A resource for identifying invasive plants

Invasive plants have aroused attention globally for causing ecological damage and having a negative impact on the economy and human health. However, it can be extremely challenging to rapidly and accurately identify invasive plants based on morphology because they are an assemblage of many different families and many plant materials lack sufficient diagnostic characteristics during border inspections. It is therefore urgent to evaluate candidate loci and build a reliable genetic library to prevent invasive plants from entering China. In this study, five common single markers (ITS, ITS2, matK, rbcL and trnH-psbA) were evaluated using 634 species (including 469 invasive plant species in China, 10 new records to China, 16 potentially invasive plant species around the world but not introduced into China yet and 139 plant species native to China) based on three different methods. Our results indicated that ITS2 displayed largest intra- and interspecific divergence (1.72% and 91.46%). Based on NJ tree method, ITS2, ITS, matK, rbcL and trnH-psbA provided 76.84%, 76.5%, 63.21%, 52.86% and 50.68% discrimination rates, respectively. The combination of ITS + matK performed best and provided 91.03% discriminatory power, followed by ITS2 + matK (85.78%). For identifying unknown individuals, ITS + matK had 100% correct identification rate based on our database, followed by ITS/ITS2 (both 93.33%) and ITS2 + matK (91.67%). Thus, we propose ITS/ITS2 + matK as the most suitable barcode for invasive plants in China. This study also demonstrated that DNA barcoding is an efficient tool for identifying invasive species.

Are diversification rates and chromosome evolution in the temperate grasses (Pooideae) associated with major environmental changes in the Oligocene-Miocene?

The Pooideae are a highly diverse C3 grass subfamily that includes some of the most economically important crops, nested within the highly speciose core-pooid clade. Here, we build and explore the phylogeny of the Pooideae within a temporal framework, assessing its patterns of diversification and its chromosomal evolutionary changes in the light of past environmental transformations. We sequenced five plastid DNA loci, two coding (ndhF, matk) and three non-coding (trnH-psbA, trnT-L and trnL-F), in 163 Poaceae taxa, including representatives for all subfamilies of the grasses and all but four ingroup Pooideae tribes. Parsimony and Bayesian phylogenetic analyses were conducted and divergence times were inferred in BEAST using a relaxed molecular clock. Diversification rates were assessed using the MEDUSA approach, and chromosome evolution was analyzed using the chromEvol software. Diversification of the Pooideae started in the Late-Eocene and was especially intense during the Oligocene-Miocene. The background diversification rate increased significantly at the time of the origin of the Poodae + Triticodae clade. This shift in diversification occurred in a context of falling temperatures that potentially increased ecological opportunities for grasses adapted to open areas around the world. The base haploid chromosome number n = 7 has remained stable throughout the phylogenetic history of the core pooids and we found no link between chromosome transitions and major diversification events in the Pooideae.

Dated tribe-wide whole chloroplast genome phylogeny indicates recurrent hybridizations within Triticeae

BACKGROUND

Triticeae, the tribe of wheat grasses, harbours the cereals barley, rye and wheat and their wild relatives. Although economically important, relationships within the tribe are still not understood. We analysed the phylogeny of chloroplast lineages among nearly all monogenomic Triticeae taxa and polyploid wheat species aiming at a deeper understanding of the tribe's evolution. We used on- and off-target reads of a target-enrichment experiment followed by Illumina sequencing.

RESULTS

The read data was used to assemble the plastid locus ndhF for 194 individuals and the whole chloroplast genome for 183 individuals, representing 53 Triticeae species and 15 genera. We conducted Bayesian and multispecies coalescent analyses to infer relationships and estimate divergence times of the taxa. We present the most comprehensive dated Triticeae chloroplast phylogeny and review previous hypotheses in the framework of our results. Monophyly of Triticeae chloroplasts could not be confirmed, as either Bromus or Psathyrostachys captured a chloroplast from a lineage closely related to a Bromus-Triticeae ancestor. The most recent common ancestor of Triticeae occurred approximately between ten and 19 million years ago.

CONCLUSIONS

The comparison of the chloroplast phylogeny with available nuclear data in several cases revealed incongruences indicating past hybridizations. Recent events of chloroplast capture were detected as individuals grouped apart from con-specific accessions in otherwise monopyhletic groups.

Evolution of the miR5200-FLOWERING LOCUS T flowering time regulon in the temperate grass subfamily Pooideae

Flowering time is a carefully regulated trait controlled primarily through the action of the central genetic regulator, FLOWERING LOCUS T (FT). Recently it was demonstrated that a microRNA, miR5200, targets the end of the second exon of FT under short-day photoperiods in the grass subfamily Pooideae, thus preventing FT transcripts from reaching threshold levels under non-inductive conditions. Pooideae are an interesting group in that they rapidly diversified from the tropics into the northern temperate region during a major global cooling event spanning the Eocene-Oligocene transition. We hypothesize that miR5200 photoperiod-sensitive regulation of Pooideae flowering time networks assisted their transition into northern seasonal environments. Here, we test predictions derived from this hypothesis that miR5200, originally found in bread wheat and later identified in Brachypodium distachyon, (1) was present in the genome of the Pooideae common ancestor, (2) is transcriptionally regulated by photoperiod, and (3) is negatively correlated with FT transcript abundance, indicative of miR5200 regulating FT. Our results demonstrate that miR5200 did evolve at or around the base of Pooideae, but only acquired photoperiod-regulated transcription within the Brachypodium lineage. Based on expression profiles and previous data, we posit that the progenitor of miR5200 was co-regulated with FT by an unknown mechanism.

The First Complete Plastid Genome from Joinvilleaceae (J. ascendens; Poales) Shows Unique and Unpredicted Rearrangements

Joinvilleaceae is a family of tropical grass-like monocots that comprises only the genus Joinvillea. Previous studies have placed Joinvilleaceae in close phylogenetic proximity to the well-studied grass family. A full plastome sequence was determined and characterized for J. ascendens. The plastome was sequenced with next generation methods, fully assembled de novo and annotated. The assembly revealed two novel inversions specific to the Joinvilleaceae lineage and at least one novel plastid inversion in the Joinvilleaceae-Poaceae lineage. Two previously documented inversions in the Joinvilleaceae-Poaceae lineage and one previously documented inversion in the Poaceae lineage were also verified. Inversion events were identified visually and verified computationally by simulation mutations. Additionally, the loss and subsequent degradation of the accD gene in order Poales was explored extensively in Poaceae and J. ascendens. The two novel inversions along with changes in gene composition between families better delimited lineages in the Poales. The presence of large inversions and subsequent reversals in this small family suggested a high potential for large-scale rearrangements to occur in plastid genomes.

Evolutionary relationships in Panicoid grasses based on plastome phylogenomics (Panicoideae; Poaceae)

BACKGROUND

Panicoideae are the second largest subfamily in Poaceae (grass family), with 212 genera and approximately 3316 species. Previous studies have begun to reveal relationships within the subfamily, but largely lack resolution and/or robust support for certain tribal and subtribal groups. This study aims to resolve these relationships, as well as characterize a putative mitochondrial insert in one linage.

RESULTS

35 newly sequenced Panicoideae plastomes were combined in a phylogenomic study with 37 other species: 15 Panicoideae and 22 from outgroups. A robust Panicoideae topology largely congruent with previous studies was obtained, but with some incongruences with previously reported subtribal relationships. A mitochondrial DNA (mtDNA) to plastid DNA (ptDNA) transfer was discovered in the Paspalum lineage.

CONCLUSIONS

The phylogenomic analysis returned a topology that largely supports previous studies. Five previously recognized subtribes appear on the topology to be non-monophyletic. Additionally, evidence for mtDNA to ptDNA transfer was identified in both Paspalum fimbriatum and P. dilatatum, and suggests a single rare event that took place in a common progenitor. Finally, the framework from this study can guide larger whole plastome sampling to discern the relationships in Cyperochloeae, Steyermarkochloeae, Gynerieae, and other incertae sedis taxa that are weakly supported or unresolved.

Phylogenomics and Plastome Evolution of Tropical Forest Grasses (Leptaspis, Streptochaeta: Poaceae)

Studies of complete plastomes have proven informative for our understanding of the molecular evolution and phylogenomics of grasses. In this study, a plastome phylogenomic analysis sampled species from lineages of deeply diverging grasses including Streptochaeta spicata (Anomochlooideae), Leptaspis banksii, and L. zeylanica (both Pharoideae). Plastomes from next generation sequences for three species were assembled by de novo methods. The unambiguously aligned coding and non-coding sequences of the entire plastomes were aligned with those from 43 other grasses and the outgroup Joinvillea ascendens. Outgroup sampling of grasses has previously posed a challenge for plastome phylogenomic studies because of major rearrangements of the plastome. Here, over 81,000 bases of homologous sequence were aligned for phylogenomic and divergence estimation analyses. Rare genomic changes, including persistently long ψycf1 and ψycf2 loci, the loss of the rpoC1 intron, and a 21 base tandem repeat insert in the coding sequence for rps19 defined branch points in the grass phylogeny. Marked differences were seen in the topologies inferred from the complete plastome and two gene matrices, and mean maximum likelihood support values for the former were 10% higher. In the full plastome phylogenomic analyses, the two species of Anomochlooideae were monophyletic. Leptaspis and Pharus were found to be reciprocally monophyletic, with the estimated divergence of two Leptaspis species preceding those of Pharus by over 14 Ma, consistent with historical biogeography. Our estimates for deep divergences among grasses were older than previous such estimates, likely influenced by more complete taxonomic and molecular sampling and the use of recently available or previously unused fossil calibration points.