Chloroplast genome variation in upland and lowland switchgrass

Characterization of the organellar genomes of Mazzaella laminarioides and Mazzaella membranacea (Gigartinaceae, Rhodophyta)

Mazzaella, a genus with no genomic resources available, has extensive distribution in the cold waters of the Pacific, where they represent ecologically and economically important species. In this study, we aimed to sequence, assemble, and annotate the complete mitochondrial and chloroplast genomes from two Mazzaella spp. and characterize the intraspecific variation among them. We report for the first time seven whole organellar genomes (mitochondria: OR915856, OR947465, OR947466, OR947467, OR947468, OR947469, OR947470; chloroplast: OR881974, OR909680, OR909681, OR909682, OR909683, OR909684, OR909685) obtained through high-throughput sequencing for six M. laminarioides sampled from three Chilean regions and one M. membranacea. Sequenced Mazzaella mitogenomes have identical gene number, gene order, and genome structure. The same results were observed for assembled plastomes. A total of 52 genes were identified in mitogenomes, and a total of 235 genes were identified in plastomes. Although the M. membranacea plastome included a full-length pbsA gene, in all M. laminarioides samples, the pbsA gene was split in three open reading frames (ORFs). Within M. laminarioides, we observed important plastome lineage-specific variations, such as the pseudogenization of the two hypothetical protein-coding genes, ycf23 and ycf45. Nonsense mutations in the ycf23 and ycf45 genes were only detected in the northern lineage. These results are consistent with phylogenetic reconstructions and divergence time estimation using concatenated coding sequences that not only support the monophyly of M. laminarioides but also underscore that the three M. laminarioides lineages are in an advanced stage of divergence. These new results open the question of the existence of still undisclosed species in M. laminarioides.

The complete Chloroplast genome of Stachys geobombycis and comparative analysis with related Stachys species

Herb genomics, at the forefront of traditional Chinese medicine research, combines genomics with traditional practices, facilitating the scientific validation of ancient remedies. This integration enhances public understanding of traditional Chinese medicine's efficacy and broadens its scope in modern healthcare. Stachys species encompass annual or perennial herbs or small shrubs, exhibiting simple petiolate or sessile leaves. Despite their wide-ranging applications across various fields, molecular data have been lacking, hindering the precise identification and taxonomic elucidation of Stachys species. To address this gap, we assembled the complete chloroplast (CP) genome of Stachys geobombycis and conducted reannotation and comparative analysis of seven additional species within the Stachys genus. The findings demonstrate that the CP genomes of these species exhibit quadripartite structures, with lengths ranging from 14,523 to 150,599 bp. Overall, the genome structure remains relatively conserved, hosting 131 annotated genes, including 87 protein coding genes, 36 tRNA genes, and 8 rRNA genes. Additionally, 78 to 98 SSRs and long repeat sequences were detected , and notably, 6 highly variable regions were identified as potential molecular markers in the CP genome through sequence alignment. Phylogenetic analysis based on Bayesian inference and maximum likelihood methods strongly supported the phylogenetic position of the genus Stachys as a member of Stachydeae tribe. Overall, this comprehensive bioinformatics study of Stachys CP genomes lays the groundwork for phylogenetic classification, plant identification, genetic engineering, evolutionary studies, and breeding research concerning medicinal plants within the Stachys genus.

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.

A Comprehensive Evolutionary Study of Chloroplast RNA Editing in Gymnosperms: A Novel Type of G-to-A RNA Editing Is Common in Gymnosperms

Although more than 9100 plant plastomes have been sequenced, RNA editing sites of the whole plastome have been experimentally verified in only approximately 21 species, which seriously hampers the comprehensive evolutionary study of chloroplast RNA editing. We investigated the evolutionary pattern of chloroplast RNA editing sites in 19 species from all 13 families of gymnosperms based on a combination of genomic and transcriptomic data. We found that the chloroplast C-to-U RNA editing sites of gymnosperms shared many common characteristics with those of other land plants, but also exhibited many unique characteristics. In contrast to that noted in angiosperms, the density of RNA editing sites in ndh genes was not the highest in the sampled gymnosperms, and both loss and gain events at editing sites occurred frequently during the evolution of gymnosperms. In addition, GC content and plastomic size were positively correlated with the number of chloroplast RNA editing sites in gymnosperms, suggesting that the increase in GC content could provide more materials for RNA editing and facilitate the evolution of RNA editing in land plants or vice versa. Interestingly, novel G-to-A RNA editing events were commonly found in all sampled gymnosperm species, and G-to-A RNA editing exhibits many different characteristics from C-to-U RNA editing in gymnosperms. This study revealed a comprehensive evolutionary scenario for chloroplast RNA editing sites in gymnosperms, and reported that a novel type of G-to-A RNA editing is prevalent in gymnosperms.

Haplotype Analysis of Chloroplast Genomes for Jujube Breeding

Jujube (family Rhamnaceae) is an important economic fruit tree in China. In this study, we reported 26 chloroplast (cp) sequences of jujube using Illumina paired-end sequencing. The sequence length of cp genome was 161, 367-161, 849 bp, which was composed of a large single-copy region (89053-89437 bp) and a small single-copy region (19356-19362 bp) separated by a pair of reverse repeat regions (26478-26533 bp). Each cp genome encodes the same 130 genes, including 112 unique genes, being quite conserved in genome structure and gene sequence. A total of 118 single base substitutions (SNPs) and 130 InDels were detected in 65 jujube accessions. Phylogenetic and haplotype network construction methods were used to analyze the origin and evolution of jujube and its sour-tasting relatives. We detected 32 effective haplotypes, consisting of 20 unique jujube haplotypes and 9 unique sour-jujube haplotypes. Compared with sour-jujube, jujube showed greater haplotype diversity at the chloroplast DNA level. To cultivate crisp and sweet fruit varieties featuring strong resistance, by combining the characteristics of sour-jujube and cultivated jujube, three hybrid combinations were suggested for reciprocal crosses: "Dongzao" × "Jingzao39," "Dongzao" × "Jingzao60," "Dongzao" × "Jingzao28." This study provides the basis for jujube species' identification and breeding, and lays the foundation for future research.

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

BACKGROUND AND AIMS

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

METHODS

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

KEY RESULTS

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

CONCLUSIONS

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

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.

Linkage mapping evidence for a syntenic QTL associated with flowering time in perennial C4 rhizomatous grasses Miscanthus and switchgrass

Flowering in perennial species is directed via complex signalling pathways that adjust to developmental regulations and environmental cues. Synchronized flowering in certain environments is a prerequisite to commercial seed production, and so the elucidation of the genetic architecture of flowering time in Miscanthus and switchgrass could aid breeding in these underdeveloped species. In this context, we assessed a mapping population in Miscanthus and two ecologically diverse switchgrass mapping populations over 3 years from planting. Multiple flowering time quantitative trait loci (QTL) were identified in both species. Remarkably, the most significant Miscanthus and switchgrass QTL proved to be syntenic, located on linkage groups 4 and 2, with logarithm of odds scores of 17.05 and 21.8 respectively. These QTL regions contained three flowering time transcription factors: Squamosa Promoter-binding protein-Like, MADS-box SEPELLATA2 and gibberellin-responsive bHLH137. The former is emerging as a key component of the age-related flowering time pathway.

The complete chloroplast genome of Papaver setigerum and comparative analyses in Papaveraceae

Papaver setigerum is an annual herb that is closely related to the opium poppy, P. somniferum. Genetic resources for P. setigerum are scarce. In the present study, we assembled the complete chloroplast (cp) genome of P. setigerum based on genome skimming data, and we conducted comparative cp genome analyses to study the evolutionary pattern in Papaveraceae. The cp genome of P. setigerum is 152,862 bp in length with a typical quadripartite structure. Comparative analyses revealed no gene rearrangement in the Papaveraceae family, although differences were evident in genome size, gene losses, as well as inverted repeats (IR) region expansion and contraction. The rps15 gene has been lost from the genomes of Meconopsis racemosa, Coreanomecon hylomeconoides, P. orientale, P. somniferum, and P. setigerum, and the ycf15 gene is found only in C. hylomeconoides. Moreover, 13 cpDNA markers, including psbA-trnH, rps16-trnQ, trnS-trnG, trnC-petN, trnE-trnT, trnL-trnF, trnF-ndhJ, petA-psbJ, ndhF-rpl32, rpl32-trnL, ccsA-ndhD, ndhE-ndhG, and rps15-ycf1, were identified with relatively high levels of variation within Papaver, which will be useful for species identification in this genus. Among those markers, psbA-trnH is the best one to distinguish P. somniferum and P. setigerum.

The dynamic evolution of mobile open reading frames in plastomes of Hymenophyllum Sm. and new insight on Hymenophyllum coreanum Nakai

In this study, four plastomes of Hymenophyllum, distributed in the Korean peninsula, were newly sequenced and phylogenomic analysis was conducted to reveal (1) the evolutionary history of plastomes of early-diverging fern species at the species level, (2) the importance of mobile open reading frames in the genus, and (3) plastome sequence divergence providing support for H. coreanum to be recognized as an independent species distinct from H. polyanthos. In addition, 1C-values of H. polyanthos and H. coreanum were measured to compare the genome size of both species and to confirm the diversification between them. The rrn16-trnV intergenic regions in the genus varied in length caused by Mobile Open Reading Frames in Fern Organelles (MORFFO). We investigated enlarged noncoding regions containing MORFFO throughout the fern plastomes and found that they were strongly associated with tRNA genes or palindromic elements. Sequence identity between plastomes of H. polyanthos and H. coreanum is quite low at 93.35% in the whole sequence and 98.13% even if the variation in trnV-rrn16 intergenic spacer was ignored. In addition, different genome sizes were found for these species based on the 1C-value. Consequently, there is no reason to consider them as a conspecies.

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.

The Complete Chloroplast Genome Sequencing and Comparative Analysis of Reed Canary Grass (Phalaris arundinacea) and Hardinggrass (P. aquatica)

There are 22 species in the Phalaris genera that distribute almost all over the temperate regions of the world. Among them, reed canary grass (Phalaris arundinacea, tetraploid and hexaploid) and hardinggrass (P. aquatica, tetraploid) have been long cultivated as forage grass and have received attention as bio-energy materials in recent years. We aimed to facilitate inter-species/ploidies comparisons, and to illuminate the degree of sequence variation within existing gene pools, chloroplast (cp) genomes of three Phalaris cytotypes (P. aquatica/4x, P. arundinacea/4x and P. arundinacea/6x) were sequenced and assembled. The result indicated that certain sequence variations existed between the cp genomes of P. arundinacea and P. aquatica. Several hotspot regions (atpI~atpH, trnT-UGU~ndhJ, rbcL~psaI, and ndhF~rpl32) were found, and variable genes (infA, psaI, psbK, etc.) were detected. SNPs (single nucleotide polymorphisms) and/or indels (insertions and deletions) were confirmed by the high Ka/Ks and Pi value. Furthermore, distribution and presence of cp simple sequence repeats (cpSSRs) were identified in the three Phalaris cp genomes, although little difference was found between hexaploid and tetraploid P. arundinacea, and no rearrangement was detected among the three Phalaris cp genomes. The evolutionary relationship and divergent time among these species were discussed. The RNA-seq revealed several differentially expressed genes (DEGs), among which psaA, psaB, and psbB related to leaf color were further verified by leaf color differences.

Transcriptome divergence during leaf development in two contrasting switchgrass (Panicum virgatum L.) cultivars

The genetics and responses to biotic stressors of tetraploid switchgrass (Panicum virgatum L.) lowland cultivar 'Kanlow' and upland cultivar Summer are distinct and can be exploited for trait improvement. In general, there is a paucity of data on the basal differences in transcription across tissue developmental times for switchgrass cultivars. Here, the changes in basal and temporal expression of genes related to leaf functions were evaluated for greenhouse grown 'Kanlow', and 'Summer' plants. Three biological replicates of the 4th leaf pooled from 15 plants per replicate were harvested at regular intervals beginning from leaf emergence through senescence. Increases and decreases in leaf chlorophyll and N content were similar for both cultivars. Likewise, multidimensional scaling (MDS) analysis indicated both cultivar-independent and cultivar-specific gene expression. Cultivar-independent genes and gene-networks included those associated with leaf function, such as growth/senescence, carbon/nitrogen assimilation, photosynthesis, chlorophyll biosynthesis, and chlorophyll degradation. However, many genes encoding nucleotide-binding leucine rich repeat (NB-LRRs) proteins and wall-bound kinases associated with detecting and responding to environmental signals were differentially expressed. Several of these belonged to unique cultivar-specific gene co-expression networks. Analysis of genomic resequencing data provided several examples of NB-LRRs genes that were not expressed and/or apparently absent in the genomes of Summer plants. It is plausible that cultivar (ecotype)-specific genes and gene-networks could be one of the drivers for the documented differences in responses to leaf-borne pathogens between these two cultivars. Incorporating broad resistance to plant pathogens in elite switchgrass germplasm could improve sustainability of biomass production under low-input conditions.

Pervasive admixture between eucalypt species has consequences for conservation and assisted migration

Conservation management often uses information on genetic population structure to assess the importance of local provenancing for ecological restoration and reintroduction programs. For species that do not exhibit complete reproductive isolation, the estimation of population genetic parameters may be influenced by the extent of admixture. Therefore, to avoid perverse outcomes for conservation, genetically informed management strategies must determine whether hybridization between species is relevant, and the extent to which observed population genetic patterns are shaped by interspecific versus intraspecific gene flow. We used genotyping by sequencing to identify over 2,400 informative single nucleotide polymorphisms across 18 populations of Eucalyptus regnans F. Muell., a foundation tree species of montane forests in south-eastern Australia. We used these data to determine the extent of hybridization with another species, Eucalyptus obliqua L'Hér., and investigate how admixture influences genetic diversity parameters, by estimating metrics of genetic diversity and examining population genetic structure in datasets with and without admixed individuals. We found hybrid individuals at all sites and two highly introgressed populations. Hybrid individuals were not distributed evenly across environmental gradients, with logistic regression identifying hybrids as being associated with temperature. Removal of hybrids resulted in increases in genetic differentiation (F ST), expected heterozygosity, observed heterozygosity and the inbreeding coefficient, and different patterns of isolation by distance. After removal of hybrids and introgressed populations, mountain ash showed very little population genetic structure, with a small effect of isolation by distance, and very low global F ST(0.03). Our study shows that, in plants, decisions around provenancing of individuals for restoration depend on knowledge of whether hybridization is influencing population genetic structure. For species in which most genetic variation is held within populations, there may be little benefit in planning conservation strategies around environmental adaptation of seed sources. The possibility for adaptive introgression may also be relevant when species regularly hybridize.

Characterization of chloroplast genomes of Alnus rubra and Betula cordifolia, and their use in phylogenetic analyses in Betulaceae

OBJECTIVE

Betulaceae is a relatively small birch family that comprises about 160 deciduous trees and shrubs. Chloroplast (cp) genome sequencing of Alnus rubra and Betula cordifolia was carried out to elucidate their molecular features and phylogenetic relationship among species in Betulaceae family.

METHODS

Chloroplast genome sequencing was carried out using next generation sequencing method. Molecular and genomic features of the two cp genomes were characterized with other cp genomes in Betulaceae. Also, molecular phylogenetic analysis was performed using the whole cp genome sequences.

RESULTS

The average cp genome length was 160,136 bp among the Betulaceae species. Base compositions of the cp genomes were skewed toward a high AT ratio, with an average of 63.4%. We identified 117 different genes 83 with protein coding, 4 with ribosomal RNA, and 30 with tRNA. Eighteen genes contained introns which were conserved among the cp genomes of all Betulaceae. We mined 82 SSRs from the cp genomes of A. rubra, A. cordifolia, and A. nana. The SSRs were variable in motif repeat numbers and presence/absence among the cp genomes.

CONCLUSION

Chloroplast genome-wide sequence comparison from 11 Betulaceae species and one cp genome of evergreen oak revealed that the patterns of sequence variations were congruent with two subfamily classification Betuloideae (Alnus and Betula) and Corylaceae (Corylus, Ostrya, and Carpinus). Subsequent phylogenetic analysis also supports the sub-classifications of these species.

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.

Natural variation in genes potentially involved in plant architecture and adaptation in switchgrass (Panicum virgatum L.)

Background

Advances in genomic technologies have expanded our ability to accurately and exhaustively detect natural genomic variants that can be applied in crop improvement and to increase our knowledge of plant evolution and adaptation. Switchgrass (Panicum virgatum L.), an allotetraploid (2n = 4× = 36) perennial C4 grass (Poaceae family) native to North America and a feedstock crop for cellulosic biofuel production, has a large potential for genetic improvement due to its high genotypic and phenotypic variation. In this study, we analyzed single nucleotide polymorphism (SNP) variation in 372 switchgrass genotypes belonging to 36 accessions for 12 genes putatively involved in biomass production to investigate signatures of selection that could have led to ecotype differentiation and to population adaptation to geographic zones.

Results

A total of 11,682 SNPs were mined from ~ 15 Gb of sequence data, out of which 251 SNPs were retained after filtering. Population structure analysis largely grouped upland accessions into one subpopulation and lowland accessions into two additional subpopulations. The most frequent SNPs were in homozygous state within accessions. Sixty percent of the exonic SNPs were non-synonymous and, of these, 45% led to non-conservative amino acid changes. The non-conservative SNPs were largely in linkage disequilibrium with one haplotype being predominantly present in upland accessions while the other haplotype was commonly present in lowland accessions. Tajima’s test of neutrality indicated that PHYB, a gene involved in photoperiod response, was under positive selection in the switchgrass population. PHYB carried a SNP leading to a non-conservative amino acid change in the PAS domain, a region that acts as a sensor for light and oxygen in signal transduction.

Conclusions

Several non-conservative SNPs in genes potentially involved in plant architecture and adaptation have been identified and led to population structure and genetic differentiation of ecotypes in switchgrass. We suggest here that PHYB is a key gene involved in switchgrass natural selection. Further analyses are needed to determine whether any of the non-conservative SNPs identified play a role in the differential adaptation of upland and lowland switchgrass.

Electronic supplementary material

The online version of this article (10.1186/s12862-018-1193-2) contains supplementary material, which is available to authorized users.

Comparative chloroplast genomics between the invasive weed Mikania micrantha and its indigenous congener Mikania cordata: Structure variation, identification of highly divergent regions, divergence time estimation, and phylogenetic analysis

Mikania micrantha and Mikania cordata are the only two species in genus Mikania (Asteraceae) in China. They share very similar morphological and life-history characteristics but occupy quite different habitats. Most importantly, they generate totally different ecological consequences. While M. micrantha has become an exotic invasive weed, M. cordata exists as an indigenous species with no harmful effects on native plants or habitats. As a continuous study of our previously reported M. micrantha chloroplast (cp) genome, in this study we have further sequenced the M. cordata cp genome to (1) conduct a comparative genome analysis to gain insights into the mechanism of invasiveness; (2) develop cp markers to examine the population genetic adaptation of M. micrantha; and (3) screen variable genome regions of phylogenetic utility. The M. cordata chloroplast genome is 151,984 bp in length and displays a typical quadripartite structure. The number and distribution of protein coding genes, tRNA genes, and rRNA genes of M. cordata are identical to those of M. micrantha. The main difference lays in that the pseudogenization of ndhF and a 118-bp palindromic repeat only arises in M. cordata. Fourteen highly divergent regions, 235 base substitutions, and 58 indels were identified between the two cp genomes. Phylogenetic inferences revealed a sister relationship between M. micrantha and M. cordata whose divergence was estimated to occur around 1.78 million years ago (MYA). Twelve cpSSR loci were detected to be polymorphic and adopted to survey the genetic adaptation of M. micrantha populations. No cpSSR loci were found to undergo selection. Our results build a foundation to examine the invasive mechanism of Mikania weed.

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.

Molecular dating of phylogenetic divergence between Urochloa species based on complete chloroplast genomes

BACKGROUND

Forage species of Urochloa are planted in millions of hectares of tropical and subtropical pastures in South America. Most of the planted area is covered with four species (U. ruziziensis, U. brizantha, U. decumbens and U. humidicola). Breeding programs rely on interspecific hybridizations to increase genetic diversity and introgress traits of agronomic importance. Knowledge of phylogenetic relationships is important to optimize compatible hybridizations in Urochloa, where phylogeny has been subject of some controversy. We used next-generation sequencing to assemble the chloroplast genomes of four Urochloa species to investigate their phylogenetic relationships, compute their times of divergence and identify chloroplast DNA markers (microsatellites, SNPs and InDels).

RESULTS

Whole plastid genome sizes were 138,765 bp in U. ruziziensis, 138,945 bp in U. decumbens, 138,946 bp in U. brizantha and 138,976 bp in U. humidicola. Each Urochloa chloroplast genome contained 130 predicted coding regions and structural features that are typical of Panicoid grasses. U. brizantha and U. decumbens chloroplast sequences are highly similar and show reduced SNP, InDel and SSR polymorphism as compared to U. ruziziensis and U. humidicola. Most of the structural and sequence polymorphisms were located in intergenic regions, and reflected phylogenetic distances between species. Divergence of U. humidicola from a common ancestor with the three other Urochloa species was estimated at 9.46 mya. U. ruziziensis, U. decumbens, and U. brizantha formed a clade where the U. ruziziensis lineage would have diverged by 5.67 mya, followed by a recent divergence event between U. decumbens and U. brizantha around 1.6 mya.

CONCLUSION

Low-coverage Illumina sequencing allowed the successful sequence analysis of plastid genomes in four species of Urochloa used as forages in the tropics. Pairwise sequence comparisons detected multiple microsatellite, SNP and InDel sites prone to be used as molecular markers in genetic analysis of Urochloa. Our results placed the origin of U. humidicola and U. ruziziensis divergence in the Miocene-Pliocene boundary, and the split between U. brizantha and U. decumbens in the Pleistocene.

Complete Chloroplast Genomes of Erianthus arundinaceus and Miscanthus sinensis: Comparative Genomics and Evolution of the Saccharum Complex

The genera Erianthus and Miscanthus, both members of the Saccharum complex, are of interest as potential resources for sugarcane improvement and as bioenergy crops. Recent studies have mainly focused on the conservation and use of wild accessions of these genera as breeding materials. However, the sequence data are limited, which hampers the studies of phylogenetic relationships, population structure, and evolution of these grasses. Here, we determined the complete chloroplast genome sequences of Erianthus arundinaceus and Miscanthus sinensis by using 454 GS FLX pyrosequencing and Sanger sequencing. Alignment of the E. arundinaceus and M. sinensis chloroplast genome sequences with the known sequence of Saccharum officinarum demonstrated a high degree of conservation in gene content and order. Using the data sets of 76 chloroplast protein-coding genes, we performed phylogenetic analysis in 40 taxa including E. arundinaceus and M. sinensis. Our results show that S. officinarum is more closely related to M. sinensis than to E. arundinaceus. We estimated that E. arundinaceus diverged from the subtribe Sorghinae before the divergence of Sorghum bicolor and the common ancestor of S. officinarum and M. sinensis. This is the first report of the phylogenetic and evolutionary relationships inferred from maternally inherited variation in the Saccharum complex. Our study provides an important framework for understanding the phylogenetic relatedness of the economically important genera Erianthus, Miscanthus, and Saccharum.

The Complete Chloroplast Genome of Chinese Bayberry (Morella rubra, Myricaceae): Implications for Understanding the Evolution of Fagales

Morella rubra (Myricaceae), also known as Chinese bayberry, is an economically important, subtropical, evergreen fruit tree. The phylogenetic placement of Myricaceae within Fagales and the origin of Chinese bayberry's domestication are still unresolved. In this study, we report the chloroplast (cp) genome of M. rubra and take advantage of several previously reported chloroplast genomes from related taxa to examine patterns of evolution in Fagales. The cp genomes of three M. rubra individuals were 159,478, 159,568, and 159.586 bp in length, respectively, comprising a pair of inverted repeat (IR) regions (26,014-26,069 bp) separated by a large single-copy (LSC) region (88,683-88,809 bp) and a small single-copy (SSC) region (18,676-18,767 bp). Each cp genome encodes the same 111 unique genes, consisting of 77 different protein-coding genes, 30 transfer RNA genes and four ribosomal RNA genes, with 18 duplicated in the IRs. Comparative analysis of chloroplast genomes from four representative Fagales families revealed the loss of infA and the pseudogenization of ycf15 in all analyzed species, and rpl22 has been pseudogenized in M. rubra and Castanea mollissima, but not in Juglans regia or Ostrya rehderiana. The genome size variations are detected mainly due to the length of intergenic spacers rather than gene loss, gene pseudogenization, IR expansion or contraction. The phylogenetic relationships yielded by the complete genome sequences strongly support the placement of Myricaceae as sister to Juglandaceae. Furthermore, seven cpDNA markers (trnH-psbA, psbA-trnK, rps2-rpoC2, ycf4-cemA, petD-rpoA, ndhE-ndhG, and ndhA intron) with relatively high levels of variation and variable cpSSR loci were identified within M. rubra, which will be useful in future research characterizing the population genetics of M. rubra and investigating the origin of domesticated Chinese bayberry.

Transcriptome analysis in switchgrass discloses ecotype difference in photosynthetic efficiency

BACKGROUND

Switchgrass, a warm-season perennial grass studied as a potential dedicated biofuel feedstock, is classified into two main taxa - lowland and upland ecotypes - that differ in morphology and habitat of adaptation. But there is limited information on their inherent molecular variations.

RESULTS

Transcriptome analysis by RNA-sequencing (RNA-Seq) was conducted for lowland and upland ecotypes to document their gene expression variations. Mapping of transcriptome to the reference genome (Panicum virgatum v1.1) revealed that the lowland and upland ecotypes differ substantially in sets of genes transcribed as well as levels of expression. Differential gene expression analysis exhibited that transcripts related to photosynthesis efficiency and development and photosystem reaction center subunits were upregulated in lowlands compared to upland genotype. On the other hand, catalase isozymes, helix-loop-helix, late embryogenesis abundant group I, photosulfokinases, and S-adenosyl methionine synthase gene transcripts were upregulated in the upland compared to the lowlands. At ≥100x coverage and ≥5% minor allele frequency, a total of 25,894 and 16,979 single nucleotide polymorphism (SNP) markers were discovered for VS16 (upland ecotype) and K5 (lowland ecotype) against the reference genome. The allele combination of the SNPs revealed that the transition mutations are more prevalent than the transversion mutations.

CONCLUSIONS

The gene ontology (GO) analysis of the transcriptome indicated lowland ecotype had significantly higher representation for cellular components associated with photosynthesis machinery controlling carbon fixation. In addition, using the transcriptome data, SNP markers were detected, which were distributed throughout the genome. The differentially expressed genes and SNP markers detected in this study would be useful resources for traits mapping and gene transfer across ecotypes in switchgrass breeding for increased biomass yield for biofuel conversion.

Chloroplast genome structures in Gentiana (Gentianaceae), based on three medicinal alpine plants used in Tibetan herbal medicine

The genus Gentiana is the largest in the Gentianaceae family with ca. 400 species. However, with most species growing on the Qinghai-Tibet plateau, the processes of adaptive evolution and speciation within the genus is not clear. Also, the genomic analyses could provide important information. So far, the complete chloroplast (cp) genome data of the genus are still deficient. As the second and third sequenced members within Gentianaceae, we report the construction of complete cp sequences of Gentiana robusta King ex Hook. f. and Gentiana crassicaulis Duthie ex Burk., and describe a comparative study of three Gentiana cp genomes, including the cp genome of Gentiana straminea Maxim. published previously. These cp genomes are highly conserved in gene size, gene content, and gene order and the rps16 pseudogene with exon2 missing was found common. Three repeat types and five SSR types were investigated, and the number and distribution are similar among the three genomes. Sixteen genome divergent hotspot regions were identified across these cp genomes that could provide potential molecular markers for further phylogenetic studies in Gentiana. The IR/SC boundary organizations in Gentianales cp genomes were compared and three different types of boundaries were observed. Six data partitions of cp genomes in Gentianales were used for phylogenetic analyses and different data partitions were largely congruent with each other. The ML phylogenetic tree was constructed based on the fragments in cp genomes commonly available in 33 species from Lamiids, including 12 species in Gentianales, 1 in Boraginaceae, 10 in Solanales, and 10 in Lamiales. The result strongly supports the position of Boraginaceae (Ehretia acuminata) as the sister of Solanales, with the bootstrap values of 97 %. This study provides a platform for further research into the molecular phylogenetics of species in the order Gentianales (family Gentianaceae) notably in respect of speciation and species identification.

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

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

Genome-Wide Characterization of Major Intrinsic Proteins in Four Grass Plants and Their Non-Aqua Transport Selectivity Profiles with Comparative Perspective

Major intrinsic proteins (MIPs), commonly known as aquaporins, transport not only water in plants but also other substrates of physiological significance and heavy metals. In most of the higher plants, MIPs are divided into five subfamilies (PIPs, TIPs, NIPs, SIPs and XIPs). Herein, we identified 68, 42, 38 and 28 full-length MIPs, respectively in the genomes of four monocot grass plants, specifically Panicum virgatum, Setaria italica, Sorghum bicolor and Brachypodium distachyon. Phylogenetic analysis showed that the grass plants had only four MIP subfamilies including PIPs, TIPs, NIPs and SIPs without XIPs. Based on structural analysis of the homology models and comparing the primary selectivity-related motifs [two NPA regions, aromatic/arginine (ar/R) selectivity filter and Froger's positions (FPs)] of all plant MIPs that have been experimentally proven to transport non-aqua substrates, we predicted the transport profiles of all MIPs in the four grass plants and also in eight other plants. Groups of MIP subfamilies based on ar/R selectivity filter and FPs were linked to the non-aqua transport profiles. We further deciphered the substrate selectivity profiles of the MIPs in the four grass plants and compared them with their counterparts in rice, maize, soybean, poplar, cotton, Arabidopsis thaliana, Physcomitrella patens and Selaginella moellendorffii. In addition to two NPA regions, ar/R filter and FPs, certain residues, especially in loops B and C, contribute to the functional distinctiveness of MIP groups. Expression analysis of transcripts in different organs indicated that non-aqua transport was related to expression of MIPs since most of the unexpressed MIPs were not predicted to facilitate the transport of non-aqua molecules. Among all MIPs in every plant, TIP (BdTIP1;1, SiTIP1;2, SbTIP2;1 and PvTIP1;2) had the overall highest mean expression. Our study generates significant information for understanding the diversity, evolution, non-aqua transport profiles and insight into comparative transport selectivity of plant MIPs, and provides tools for the development of transgenic plants.

The complete chloroplast genome of Gentiana straminea (Gentianaceae), an endemic species to the Sino-Himalayan subregion

Endemic to the Sino-Himalayan subregion, the medicinal alpine plant Gentiana straminea is a threatened species. The genetic and molecular data about it is deficient. Here we report the complete chloroplast (cp) genome sequence of G. straminea, as the first sequenced member of the family Gentianaceae. The cp genome is 148,991bp in length, including a large single copy (LSC) region of 81,240bp, a small single copy (SSC) region of 17,085bp and a pair of inverted repeats (IRs) of 25,333bp. It contains 112 unique genes, including 78 protein-coding genes, 30 tRNAs and 4 rRNAs. The rps16 gene lacks exon2 between trnK-UUU and trnQ-UUG, which is the first rps16 pseudogene found in the nonparasitic plants of Asterids clade. Sequence analysis revealed the presence of 13 forward repeats, 13 palindrome repeats and 39 simple sequence repeats (SSRs). An entire cp genome comparison study of G. straminea and four other species in Gentianales was carried out. Phylogenetic analyses using maximum likelihood (ML) and maximum parsimony (MP) were performed based on 69 protein-coding genes from 36 species of Asterids. The results strongly supported the position of Gentianaceae as one member of the order Gentianales. The complete chloroplast genome sequence will provide intragenic information for its conservation and contribute to research on the genetic and phylogenetic analyses of Gentianales and Asterids.

Plastid Genotyping Reveals the Uniformity of Cytoplasmic Male Sterile-T Maize Cytoplasms

Cytoplasmic male-sterile (CMS) lines in maize (Zea mays) have been classified by their response to specific restorer genes into three categories: cms-C, cms-S, and cms-T. A mitochondrial genome representing each of the CMS cytotypes has been sequenced, and male sterility in the cms-S and cms-T cytotypes is linked to chimeric mitochondrial genes. To identify markers for plastid genotyping, we sequenced the plastid genomes of three fertile maize lines (B37, B73, and A188) and the B37 cms-C, cms-S, and cms-T cytoplasmic substitution lines. We found that the plastid genomes of B37 and B73 lines are identical. Furthermore, the fertile and CMS plastid genomes are conserved, differing only by zero to three single-nucleotide polymorphisms (SNPs) in coding regions and by eight to 22 SNPs and 10 to 21 short insertions/deletions in noncoding regions. To gain insight into the origin and transmission of the cms-T trait, we identified three SNPs unique to the cms-T plastids and tested the three diagnostic SNPs in 27 cms-T lines, representing the HA, I, Q, RS, and T male-sterile cytoplasms. We report that each of the tested 27 cms-T group accessions have the same three diagnostic plastid SNPs, indicating a single origin and maternal cotransmission of the cms-T mitochondria and plastids to the seed progeny. Our data exclude exceptional pollen transmission of organelles or multiple horizontal gene transfer events as the source of the mitochondrial urf13-T (unidentified reading frame encoding 13-kD cms-T protein) gene in the cms-T cytoplasms. Plastid genotyping enables a reassessment of the evolutionary relationships of cytoplasms in cultivated maize.

Physiological Evaluation of Alkali-Salt Tolerance of Thirty Switchgrass (Panicum virgatum) Lines

Soil salt-alkalization is a major limiting factor for crop production in many regions. Switchgrass (Panicum virgatum L.) is a warm-season C₄ perennial rhizomatous bunchgrass and a target lignocellulosic biofuel species. The objective of this study was to evaluate relative alkali-salt tolerance among 30 switchgrass lines. Tillers of each switchgrass line were transplanted into pots filled with fine sand. Two months after transplanting, plants at E5 developmental stage were grown in either half strength Hoagland’s nutrient solution with 0 mM Na⁺ (control) or half strength Hoagland’s nutrient solution with 150 mM Na⁺ and pH of 9.5 (alkali-salt stress treatment) for 20 d. Alkali-salt stress damaged cell membranes [higher electrolyte leakage (EL) ], reduced leaf relative water content (RWC), net photosynthetic rate (Pn), stomatal conductance (g_s), and transpiration rate (Tr). An alkali-salt stress tolerance trait index (ASTTI) for each parameter was calculated based on the ratio of the value under alkali-salt stress and the value under non-stress conditions for each parameter of each line. Relative alkali-salt tolerance was determined based on principal components analysis and cluster analysis of the physiological parameters and their ASTTI values. Significant differences in alkali-salt stress tolerance were found among the 30 lines. Lowland lines TEM-SEC, Alamo, TEM-SLC and Kanlow were classified as alkali-salt tolerant. In contrast, three lowland lines (AM-314/MS-155, BN-13645-64) and two upland lines (Caddo and Blackwell-1) were classified as alkali-salt sensitive. The results suggest wide variations exist in alkali-salt stress tolerance among the 30 switchgrass lines. The approach of using a combination of principal components and cluster analysis of the physiological parameters and related ASTTI is feasible for evaluating alkali-salt tolerance in switchgrass.

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

Whole plastid genomes are being sequenced rapidly from across the green plant tree of life, and phylogenetic analyses of these are increasing resolution and support for relationships that have varied among or been unresolved in earlier single- and multi-gene studies. Pooideae, the cool-season grass lineage, is the largest of the 12 grass subfamilies and includes important temperate cereals, turf grasses and forage species. Although numerous studies of the phylogeny of the subfamily have been undertaken, relationships among some 'early-diverging' tribes conflict among studies, and some relationships among subtribes of Poeae have not yet been resolved. To address these issues, we newly sequenced 25 whole plastomes, which showed rearrangements typical of Poaceae. These plastomes represent 9 tribes and 11 subtribes of Pooideae, and were analysed with 20 existing plastomes for the subfamily. Maximum likelihood (ML), maximum parsimony (MP) and Bayesian inference (BI) robustly resolve most deep relationships in the subfamily. Complete plastome data provide increased nodal support compared with protein-coding data alone at nodes that are not maximally supported. Following the divergence of Brachyelytrum, Phaenospermateae, Brylkinieae-Meliceae and Ampelodesmeae-Stipeae are the successive sister groups of the rest of the subfamily. Ampelodesmeae are nested within Stipeae in the plastome trees, consistent with its hybrid origin between a phaenospermatoid and a stipoid grass (the maternal parent). The core Pooideae are strongly supported and include Brachypodieae, a Bromeae-Triticeae clade and Poeae. Within Poeae, a novel sister group relationship between Phalaridinae and Torreyochloinae is found, and the relative branching order of this clade and Aveninae, with respect to an Agrostidinae-Brizinae clade, are discordant between MP and ML/BI trees. Maximum likelihood and Bayesian analyses strongly support Airinae and Holcinae as the successive sister groups of a Dactylidinae-Loliinae clade.

Echinochloa chloroplast genomes: insights into the evolution and taxonomic identification of two weedy species

The genus Echinochloa (Poaceae) includes numerous problematic weeds that cause the reduction of crop yield worldwide. To date, DNA sequence information is still limited in the genus Echinochloa. In this study, we completed the entire chloroplast genomes of two Echinochloa species (Echinochloa oryzicola and Echinochloa crus-galli) based on high-throughput sequencing data from their fresh green leaves. The two Echinochloa chloroplast genomes are 139,891 and 139,800 base pairs in length, respectively, and contain 131 protein-coding genes, 79 indels and 466 substitutions helpful for discrimination of the two species. The divergence between the genus Echinochloa and Panicum occurred about 21.6 million years ago, whereas the divergence between E. oryzicola and E. crus-galli chloroplast genes occurred about 3.3 million years ago. The two reported Echinochloa chloroplast genome sequences contribute to better understanding of the diversification of this genus.

Chloroplast genome differences between Asian and American Equisetum arvense (Equisetaceae) and the origin of the hypervariable trnY-trnE intergenic spacer

Comparative analyses of complete chloroplast (cp) DNA sequences within a species may provide clues to understand the population dynamics and colonization histories of plant species. Equisetum arvense (Equisetaceae) is a widely distributed fern species in northeastern Asia, Europe, and North America. The complete cp DNA sequences from Asian and American E. arvense individuals were compared in this study. The Asian E. arvense cp genome was 583 bp shorter than that of the American E. arvense. In total, 159 indels were observed between two individuals, most of which were concentrated on the hypervariable trnY-trnE intergenic spacer (IGS) in the large single-copy (LSC) region of the cp genome. This IGS region held a series of 19 bp repeating units. The numbers of the 19 bp repeat unit were responsible for 78% of the total length difference between the two cp genomes. Furthermore, only other closely related species of Equisetum also show the hypervariable nature of the trnY-trnE IGS. By contrast, only a single indel was observed in the gene coding regions: the ycf1 gene showed 24 bp differences between the two continental individuals due to a single tandem-repeat indel. A total of 165 single-nucleotide polymorphisms (SNPs) were recorded between the two cp genomes. Of these, 52 SNPs (31.5%) were distributed in coding regions, 13 SNPs (7.9%) were in introns, and 100 SNPs (60.6%) were in intergenic spacers (IGS). The overall difference between the Asian and American E. arvense cp genomes was 0.12%. Despite the relatively high genetic diversity between Asian and American E. arvense, the two populations are recognized as a single species based on their high morphological similarity. This indicated that the two regional populations have been in morphological stasis.

Complete chloroplast genome sequences of Praxelis (Eupatorium catarium Veldkamp), an important invasive species

Praxelis (Eupatorium catarium Veldkamp) is a new hazardous invasive plant species that has caused serious economic losses and environmental damage in the Northern hemisphere tropical and subtropical regions. Although previous studies focused on detecting the biological characteristics of this plant to prevent its expansion, little effort has been made to understand the impact of Praxelis on the ecosystem in an evolutionary process. The genetic information of Praxelis is required for further phylogenetic identification and evolutionary studies. Here, we report the complete Praxelis chloroplast (cp) genome sequence. The Praxelis chloroplast genome is 151,410 bp in length including a small single-copy region (18,547 bp) and a large single-copy region (85,311 bp) separated by a pair of inverted repeats (IRs; 23,776 bp). The genome contains 85 unique and 18 duplicated genes in the IR region. The gene content and organization are similar to other Asteraceae tribe cp genomes. We also analyzed the whole cp genome sequence, repeat structure, codon usage, contraction of the IR and gene structure/organization features between native and invasive Asteraceae plants, in order to understand the evolution of organelle genomes between native and invasive Asteraceae. Comparative analysis identified the 14 markers containing greater than 2% parsimony-informative characters, indicating that they are potential informative markers for barcoding and phylogenetic analysis. Moreover, a sister relationship between Praxelis and seven other species in Asteraceae was found based on phylogenetic analysis of 28 protein-coding sequences. Complete cp genome information is useful for plant phylogenetic and evolutionary studies within this invasive species and also within the Asteraceae family.

Synergistic and antagonistic effects of salinity and pH on germination in switchgrass (Panicum virgatum L.)

The effects of salt-alkaline mixed stress on switchgrass were investigated by evaluating seed germination and the proline, malondialdehyde (MDA) and soluble sugar contents in three switchgrass (Panicum virgatum L.) cultivars in order to identify which can be successfully produced on marginal lands affected by salt-alkaline mixed stress. The experimental conditions consisted of four levels of salinity (10, 60, 110 and 160 mM) and four pH levels (7.1, 8.3, 9.5 and 10.7). The effects of salt-alkaline mixed stress with equivalent coupling of the salinity and pH level on the switchgrass were explored via model analyses. Switchgrass was capable of germinating and surviving well in all treatments under low-alkaline pH (pH≤8.3), regardless of the salinity. However, seed germination and seedling growth were sharply reduced at higher pH values in conjunction with salinity. The salinity and pH had synergetic effects on the germination percentage, germination index, plumular length and the soluble sugar and proline contents in switchgrass. However, these two factors exhibited antagonistic effects on the radicular length of switchgrass. The combined effects of salinity and pH and the interactions between them should be considered when evaluating the strength of salt-alkaline mixed stress.

Insights from the complete chloroplast genome into the evolution of Sesamum indicum L

Sesame (Sesamum indicum L.) is one of the oldest oilseed crops. In order to investigate the evolutionary characters according to the Sesame Genome Project, apart from sequencing its nuclear genome, we sequenced the complete chloroplast genome of S. indicum cv. Yuzhi 11 (white seeded) using Illumina and 454 sequencing. Comparisons of chloroplast genomes between S. indicum and the 18 other higher plants were then analyzed. The chloroplast genome of cv. Yuzhi 11 contains 153,338 bp and a total of 114 unique genes (KC569603). The number of chloroplast genes in sesame is the same as that in Nicotiana tabacum, Vitis vinifera and Platanus occidentalis. The variation in the length of the large single-copy (LSC) regions and inverted repeats (IR) in sesame compared to 18 other higher plant species was the main contributor to size variation in the cp genome in these species. The 77 functional chloroplast genes, except for ycf1 and ycf2, were highly conserved. The deletion of the cp ycf1 gene sequence in cp genomes may be due either to its transfer to the nuclear genome, as has occurred in sesame, or direct deletion, as has occurred in Panax ginseng and Cucumis sativus. The sesame ycf2 gene is only 5,721 bp in length and has lost about 1,179 bp. Nucleotides 1–585 of ycf2 when queried in BLAST had hits in the sesame draft genome. Five repeats (R10, R12, R13, R14 and R17) were unique to the sesame chloroplast genome. We also found that IR contraction/expansion in the cp genome alters its rate of evolution. Chloroplast genes and repeats display the signature of convergent evolution in sesame and other species. These findings provide a foundation for further investigation of cp genome evolution in Sesamum and other higher plants.

Complete chloroplast genome of Sedum sarmentosum and chloroplast genome evolution in Saxifragales

Comparative chloroplast genome analyses are mostly carried out at lower taxonomic levels, such as the family and genus levels. At higher taxonomic levels, chloroplast genomes are generally used to reconstruct phylogenies. However, little attention has been paid to chloroplast genome evolution within orders. Here, we present the chloroplast genome of Sedum sarmentosum and take advantage of several available (or elucidated) chloroplast genomes to examine the evolution of chloroplast genomes in Saxifragales. The chloroplast genome of S. sarmentosum is 150,448 bp long and includes 82,212 bp of a large single-copy (LSC) region, 16.670 bp of a small single-copy (SSC) region, and a pair of 25,783 bp sequences of inverted repeats (IRs).The genome contains 131 unique genes, 18 of which are duplicated within the IRs. Based on a comparative analysis of chloroplast genomes from four representative Saxifragales families, we observed two gene losses and two pseudogenes in Paeonia obovata, and the loss of an intron was detected in the rps16 gene of Penthorum chinense. Comparisons among the 72 common protein-coding genes confirmed that the chloroplast genomes of S. sarmentosum and Paeonia obovata exhibit accelerated sequence evolution. Furthermore, a strong correlation was observed between the rates of genome evolution and genome size. The detected genome size variations are predominantly caused by the length of intergenic spacers, rather than losses of genes and introns, gene pseudogenization or IR expansion or contraction. The genome sizes of these species are negatively correlated with nucleotide substitution rates. Species with shorter duration of the life cycle tend to exhibit shorter chloroplast genomes than those with longer life cycles.

Transcriptome analysis of heat stress response in switchgrass (Panicum virgatum L.)

BACKGROUND

Global warming predictions indicate that temperatures will increase by another 2-6°C by the end of this century. High temperature is a major abiotic stress limiting plant growth and productivity in many areas of the world. Switchgrass (Panicum virgatum L.) is a model herbaceous bioenergy crop, due to its rapid growth rate, reliable biomass yield, minimal requirements of water and nutrients, adaptability to grow on marginal lands and widespread distribution throughout North America. The effect of high temperature on switchgrass physiology, cell wall composition and biomass yields has been reported. However, there is void in the knowledge of the molecular responses to heat stress in switchgrass.

RESULTS

We conducted long-term heat stress treatment (38°/30°C, day/night, for 50 days) in the switchgrass cultivar Alamo. A significant decrease in the plant height and total biomass was evident in the heat stressed plants compared to controls. Total RNA from control and heat stress samples were used for transcriptome analysis with switchgrass Affymetrix genechips. Following normalization and pre-processing, 5365 probesets were identified as differentially expressed using a 2-fold cutoff. Of these, 2233 probesets (2000 switchgrass unigenes) were up-regulated, and 3132 probesets (2809 unigenes) were down-regulated. Differential expression of 42 randomly selected genes from this list was validated using RT-PCR. Rice orthologs were retrieved for 78.7% of the heat stress responsive switchgrass probesets. Gene ontology (GOs) enrichment analysis using AgriGO program showed that genes related to ATPase regulator, chaperone binding, and protein folding was significantly up-regulated. GOs associated with protein modification, transcription, phosphorus and nitrogen metabolic processes, were significantly down-regulated by heat stress.

CONCLUSIONS

Plausible connections were identified between the identified GOs, physiological responses and heat response phenotype observed in switchgrass plants. Comparative transcriptome analysis in response to heat stress among four monocots - switchgrass, rice, wheat and maize identified 16 common genes, most of which were associated with protein refolding processes. These core genes will be valuable biomarkers for identifying heat sensitive plant germplasm since they are responsive to both short duration as well as chronic heat stress treatments, and are also expressed in different plant growth stages and tissue types.

Advances in biotechnology and genomics of switchgrass

Switchgrass (Panicum virgatum L.) is a C4 perennial warm season grass indigenous to the North American tallgrass prairie. A number of its natural and agronomic traits, including adaptation to a wide geographical distribution, low nutrient requirements and production costs, high water use efficiency, high biomass potential, ease of harvesting, and potential for carbon storage, make it an attractive dedicated biomass crop for biofuel production. We believe that genetic improvements using biotechnology will be important to realize the potential of the biomass and biofuel-related uses of switchgrass. Tissue culture techniques aimed at rapid propagation of switchgrass and genetic transformation protocols have been developed. Rapid progress in genome sequencing and bioinformatics has provided efficient strategies to identify, tag, clone and manipulate many economically-important genes, including those related to higher biomass, saccharification efficiency, and lignin biosynthesis. Application of the best genetic tools should render improved switchgrass that will be more economically and environmentally sustainable as a lignocellulosic bioenergy feedstock.

Comparative analyses of plastid sequences between native and introduced populations of aquatic weeds Elodea canadensis and E. nuttallii

Non-indigenous species (NIS) are species living outside their historic or native range. Invasive NIS often cause severe environmental impacts, and may have large economical and social consequences. Elodea (Hydrocharitaceae) is a New World genus with at least five submerged aquatic angiosperm species living in fresh water environments. Our aim was to survey the geographical distribution of cpDNA haplotypes within the native and introduced ranges of invasive aquatic weeds Elodea canadensis and E. nuttallii and to reconstruct the spreading histories of these invasive species. In order to reveal informative chloroplast (cp) genome regions for phylogeographic analyses, we compared the plastid sequences of native and introduced individuals of E. canadensis. In total, we found 235 variable sites (186 SNPs, 47 indels and two inversions) between the two plastid sequences consisting of 112,193 bp and developed primers flanking the most variable genomic areas. These 29 primer pairs were used to compare the level and pattern of intraspecific variation within E. canadensis to interspecific variation between E. canadensis and E. nuttallii. Nine potentially informative primer pairs were used to analyze the phylogeographic structure of both Elodea species, based on 70 E. canadensis and 25 E. nuttallii individuals covering native and introduced distributions. On the whole, the level of variation between the two Elodea species was 53% higher than that within E. canadensis. In our phylogeographic analysis, only a single haplotype was found in the introduced range in both species. These haplotypes H1 (E. canadensis) and A (E. nuttallii) were also widespread in the native range, covering the majority of native populations analyzed. Therefore, we were not able to identify either the geographic origin of the introduced populations or test the hypothesis of single versus multiple introductions. The divergence between E. canadensis haplotypes was surprisingly high, and future research may clarify mechanisms that structure native E. canadensis populations.

Complete chloroplast genome sequences of Mongolia medicine Artemisia frigida and phylogenetic relationships with other plants

Background

Artemisia frigida Willd. is an important Mongolian traditional medicinal plant with pharmacological functions of stanch and detumescence. However, there is little sequence and genomic information available for Artemisia frigida, which makes phylogenetic identification, evolutionary studies, and genetic improvement of its value very difficult. We report the complete chloroplast genome sequence of Artemisia frigida based on 454 pyrosequencing.

Methodology/Principal Findings

The complete chloroplast genome of Artemisia frigida is 151,076 bp including a large single copy (LSC) region of 82,740 bp, a small single copy (SSC) region of 18,394 bp and a pair of inverted repeats (IRs) of 24,971 bp. The genome contains 114 unique genes and 18 duplicated genes. The chloroplast genome of Artemisia frigida contains a small 3.4 kb inversion within a large 23 kb inversion in the LSC region, a unique feature in Asteraceae. The gene order in the SSC region of Artemisia frigida is inverted compared with the other 6 Asteraceae species with the chloroplast genomes sequenced. This inversion is likely caused by an intramolecular recombination event only occurred in Artemisia frigida. The existence of rich SSR loci in the Artemisia frigida chloroplast genome provides a rare opportunity to study population genetics of this Mongolian medicinal plant. Phylogenetic analysis demonstrates a sister relationship between Artemisia frigida and four other species in Asteraceae, including Ageratina adenophora, Helianthus annuus, Guizotia abyssinica and Lactuca sativa, based on 61 protein-coding sequences. Furthermore, Artemisia frigida was placed in the tribe Anthemideae in the subfamily Asteroideae (Asteraceae) based on ndhF and trnL-F sequence comparisons.

Conclusion

The chloroplast genome sequence of Artemisia frigida was assembled and analyzed in this study, representing the first plastid genome sequenced in the Anthemideae tribe. This complete chloroplast genome sequence will be useful for molecular ecology and molecular phylogeny studies within Artemisia species and also within the Asteraceae family.

Switchgrass genomic diversity, ploidy, and evolution: novel insights from a network-based SNP discovery protocol

Switchgrass (Panicum virgatum L.) is a perennial grass that has been designated as an herbaceous model biofuel crop for the United States of America. To facilitate accelerated breeding programs of switchgrass, we developed both an association panel and linkage populations for genome-wide association study (GWAS) and genomic selection (GS). All of the 840 individuals were then genotyped using genotyping by sequencing (GBS), generating 350 GB of sequence in total. As a highly heterozygous polyploid (tetraploid and octoploid) species lacking a reference genome, switchgrass is highly intractable with earlier methodologies of single nucleotide polymorphism (SNP) discovery. To access the genetic diversity of species like switchgrass, we developed a SNP discovery pipeline based on a network approach called the Universal Network-Enabled Analysis Kit (UNEAK). Complexities that hinder single nucleotide polymorphism discovery, such as repeats, paralogs, and sequencing errors, are easily resolved with UNEAK. Here, 1.2 million putative SNPs were discovered in a diverse collection of primarily upland, northern-adapted switchgrass populations. Further analysis of this data set revealed the fundamentally diploid nature of tetraploid switchgrass. Taking advantage of the high conservation of genome structure between switchgrass and foxtail millet (Setaria italica (L.) P. Beauv.), two parent-specific, synteny-based, ultra high-density linkage maps containing a total of 88,217 SNPs were constructed. Also, our results showed clear patterns of isolation-by-distance and isolation-by-ploidy in natural populations of switchgrass. Phylogenetic analysis supported a general south-to-north migration path of switchgrass. In addition, this analysis suggested that upland tetraploid arose from upland octoploid. All together, this study provides unparalleled insights into the diversity, genomic complexity, population structure, phylogeny, phylogeography, ploidy, and evolutionary dynamics of switchgrass.

Recent advances in sequencing technologies have enabled large-scale surveys of genetic diversity in model species with a wholly or partly sequenced reference genome. However, thousands of key species, which are essential for food, health, energy, and ecology, do not have reference genomes. To accelerate their breeding cycle via marker assisted selection, high-throughput genotyping is required for these valuable species, in spite of the absence of reference genomes. Based on genotyping by sequencing (GBS) technology, we developed a new single nucleotide polymorphism (SNP) discovery protocol, the Universal Network-Enabled Analysis Kit (UNEAK), which can be widely used in any species, regardless of genome complexity or the availability of a reference genome. Here we test this protocol on switchgrass, currently the prime energy crop species in the United States of America. In addition to the discovery of over a million SNPs and construction of high-density linkage maps, we provide novel insights into the genome complexity, ploidy, phylogeny, and evolution of switchgrass. This is only the beginning: we believe UNEAK offers the key to the exploration and exploitation of the genetic diversity of thousands of non-model species.

Switchgrass genomic diversity, ploidy, and evolution: novel insights from a network-based SNP discovery protocol

Switchgrass (Panicum virgatum L.) is a perennial grass that has been designated as an herbaceous model biofuel crop for the United States of America. To facilitate accelerated breeding programs of switchgrass, we developed both an association panel and linkage populations for genome-wide association study (GWAS) and genomic selection (GS). All of the 840 individuals were then genotyped using genotyping by sequencing (GBS), generating 350 GB of sequence in total. As a highly heterozygous polyploid (tetraploid and octoploid) species lacking a reference genome, switchgrass is highly intractable with earlier methodologies of single nucleotide polymorphism (SNP) discovery. To access the genetic diversity of species like switchgrass, we developed a SNP discovery pipeline based on a network approach called the Universal Network-Enabled Analysis Kit (UNEAK). Complexities that hinder single nucleotide polymorphism discovery, such as repeats, paralogs, and sequencing errors, are easily resolved with UNEAK. Here, 1.2 million putative SNPs were discovered in a diverse collection of primarily upland, northern-adapted switchgrass populations. Further analysis of this data set revealed the fundamentally diploid nature of tetraploid switchgrass. Taking advantage of the high conservation of genome structure between switchgrass and foxtail millet (Setaria italica (L.) P. Beauv.), two parent-specific, synteny-based, ultra high-density linkage maps containing a total of 88,217 SNPs were constructed. Also, our results showed clear patterns of isolation-by-distance and isolation-by-ploidy in natural populations of switchgrass. Phylogenetic analysis supported a general south-to-north migration path of switchgrass. In addition, this analysis suggested that upland tetraploid arose from upland octoploid. All together, this study provides unparalleled insights into the diversity, genomic complexity, population structure, phylogeny, phylogeography, ploidy, and evolutionary dynamics of switchgrass.

The organelle genomes of Hassawi rice (Oryza sativa L.) and its hybrid in saudi arabia: genome variation, rearrangement, and origins

Hassawi rice (Oryza sativa L.) is a landrace adapted to the climate of Saudi Arabia, characterized by its strong resistance to soil salinity and drought. Using high quality sequencing reads extracted from raw data of a whole genome sequencing project, we assembled both chloroplast (cp) and mitochondrial (mt) genomes of the wild-type Hassawi rice (Hassawi-1) and its dwarf hybrid (Hassawi-2). We discovered 16 InDels (insertions and deletions) but no SNP (single nucleotide polymorphism) is present between the two Hassawi cp genomes. We identified 48 InDels and 26 SNPs in the two Hassawi mt genomes and a new type of sequence variation, termed reverse complementary variation (RCV) in the rice cp genomes. There are two and four RCVs identified in Hassawi-1 when compared to 93–11 (indica) and Nipponbare (japonica), respectively. Microsatellite sequence analysis showed there are more SSRs in the genic regions of both cp and mt genomes in the Hassawi rice than in the other rice varieties. There are also large repeats in the Hassawi mt genomes, with the longest length of 96,168 bp and 96,165 bp in Hassawi-1 and Hassawi-2, respectively. We believe that frequent DNA rearrangement in the Hassawi mt and cp genomes indicate ongoing dynamic processes to reach genetic stability under strong environmental pressures. Based on sequence variation analysis and the breeding history, we suggest that both Hassawi-1 and Hassawi-2 originated from the Indonesian variety Peta since genetic diversity between the two Hassawi cultivars is very low albeit an unknown historic origin of the wild-type Hassawi rice.

Karyotype variation is indicative of subgenomic and ecotypic differentiation in switchgrass

BACKGROUND

Karyotypes can provide information about taxonomic relationships, genetic aberrations, and the evolutionary origins of species. However, differentiation of the tiny chromosomes of switchgrass (Panicum virgatum L.) and creation of a standard karyotype for this bioenergy crop has not been accomplished due to lack of distinguishing features and polyploidy.

RESULTS

A cytogenetic study was conducted on a dihaploid individual (2n = 2X = 18) of switchgrass to establish a chromosome karyotype. Size differences, condensation patterns, and arm-length ratios were used as identifying features and fluorescence in-situ hybridization (FISH) assigned 5S and 45S rDNA loci to chromosomes 7 and 2 respectively. Both a maize CentC and a native switchgrass centromeric repeat (PviCentC) that shared 73% sequence identity demonstrated a strong signal on chromosome 3. However, only the PviCentC probe labeled the centromeres of all chromosomes. Unexpected PviCentC and 5S rDNA hybidization patterns were consistent with severe reduction or total deletion of these repeats in one subgenome. These patterns were maintained in tetraploid and octoploid individuals. The 45S rDNA repeat produced the expected number of loci in dihaploid, tetraploid and octoploid individuals. Differences observed at the 5S rDNA loci between the upland and lowland ecotypes of switchgrass provided a basis for distinguishing these subpopulations.

CONCLUSION

Collectively, these results provide a quantitative karyotype of switchgrass chromosomes. FISH analyses indicate genetic divergence between subgenomes and allow for the classification of switchgrass plants belonging to divergent genetic pools. Furthermore, the karyotype structure and cytogenetic analysis of switchgrass provides a framework for future genetic and genomic studies.