The chloroplast is one of two organelles containing a separate genome that codes for essential and distinct cellular functions such as photosynthesis. Given the importance of chloroplasts in plant metabolism, the genomic architecture and gene content have been strongly conserved through long periods of time and as such are useful molecular tools for evolutionary inferences. At present, complete chloroplast genomes from over 4000 species have been deposited into publicly accessible databases. Despite the large number of complete chloroplast genomes, comprehensive analyses regarding genome architecture and gene content have not been conducted for many lineages with complete species sampling. In this study, we employed the genus Populus to assess how more comprehensively sampled chloroplast genome analyses can be used in understanding chloroplast evolution in a broadly studied lineage of angiosperms. We conducted comparative analyses across Populus in order to elucidate variation in key genome features such as genome size, gene number, gene content, repeat type and number, SSR (Simple Sequence Repeat) abundance, and boundary positioning between the four main units of the genome. We found that some genome annotations were variable across the genus owing in part from errors in assembly or data checking and from this provided corrected annotations. We also employed complete chloroplast genomes for phylogenetic analyses including the dating of divergence times throughout the genus. Lastly, we utilized re-sequencing data to describe the variations of pan-chloroplast genomes at the population level for P. euphratica. The analyses used in this paper provide a blueprint for the types of analyses that can be conducted with publicly available chloroplast genomes as well as methods for building upon existing datasets to improve evolutionary inference.

The genus Dalbergia contains more than 200 species, several of which are trees that produce traditional medicines and extremely high-value timber commonly referred to as rosewood. Due to the rarity of these species in the wild, the high value of the timber, and a growing international illicit trade, CITES (Convention on International Trade in Endangered Species of Wild Fauna and Flora) has listed the entire genus in appendix II and the species Dalbergia nigra in appendix I because species in this genus are considered at risk of extinction. Given this, and the fact that species or even genus level determination is nearly impossible from cut timber morphology, alternative molecular methods are needed to identify and track intercepted rosewood. To better identify rosewood using molecular methods, we sequenced and assembled eight chloroplast genomes including D. nigra as well as conducted comparative analyses with all other available chloroplast genomes in Dalbergia and closely related lineages. From these analyses, numerous repeats including simple sequence repeats (SSR) and conserved nucleotide polymorphisms unique to subclades within the genus were detected. From phylogenetic analysis based on the CDS from 77 chloroplast genes, the groups Siam rosewood and scented rosewood resolved as monophyletic, supporting the morphological traits used to delimit these species. In addition, several instances of paraphyly and polyphyly resulting from mismatches between taxonomic determinations and phylogenetic tree topology were identified. Ultimately, the highly variable regions in the chloroplast genomes will provide useful plastid markers for further studies regarding the identification, phylogeny, and population genetics of Dalbergia species, including those frequently intercepted in illegal trade.

Bougainvillea is a perennial ornamental shrub that is highly regarded in ornamental horticulture around the world. However, the absence of genome data limits our understanding of the pathways involved in bract coloration and breeding. Here, we report a chromosome-level assembly of the giga-genome of Bougainvillea × buttiana ‘Mrs Butt’, a cultivar thought to be the origin of many other Bougainvillea cultivars. The assembled genome is ~5 Gb with a scaffold N50 of 151 756 278 bp and contains 86 572 genes which have undergone recent whole-genome duplication. We confirmed that multiple rounds of whole-genome multiplication have occurred in the evolutionary history of the Caryophyllales, reconstructed the relationship in the Caryophyllales at whole genome level, and found discordance between species and gene trees as the result of complex introgression events. We investigated betalain and anthocyanin biosynthetic pathways and found instances of independent evolutionary innovations in the nine different Caryophyllales species. To explore the potential formation mechanism of diverse bract colors in Bougainvillea, we analyzed the genes involved in betalain and anthocyanin biosynthesis and found extremely low expression of ANS and DFR genes in all cultivars, which may limit anthocyanin biosynthesis. Our findings indicate that the expression pattern of the betalain biosynthetic pathway did not directly correlate with bract color, and a higher expression level in the betalain biosynthetic pathway is required for colored bracts. This improved understanding of the correlation between gene expression and bract color allows plant breeding outcomes to be predicted with greater certainty.

Compared with most flowers where the showy part comprises specialized leaves (petals) directly subtending the reproductive structures, most Zingiberaceae species produce showy “flowers” through modifications of leaves (bracts) subtending the true flowers throughout an inflorescence. Curcuma alismatifolia, belonging to the Zingiberaceae family, a plant species originating from Southeast Asia, has become increasingly popular in the flower market worldwide because of its varied and esthetically pleasing bracts produced in different cultivars. Here, we present the chromosome-scale genome assembly of C. alismatifolia “Chiang Mai Pink” and explore the underlying mechanisms of bract pigmentation. Comparative genomic analysis revealed C. alismatifolia contains a residual signal of whole-genome duplication. Duplicated genes, including pigment-related genes, exhibit functional and structural differentiation resulting in diverse bract colors among C. alismatifolia cultivars. In addition, we identified the key genes that produce different colored bracts in C. alismatifolia, such as F3′5'H, DFR, ANS and several transcription factors for anthocyanin synthesis, as well as chlH and CAO in the chlorophyll synthesis pathway by conducting transcriptomic analysis, bulked segregant analysis using both DNA and RNA data, and population genomic analysis. This work provides data for understanding the mechanism of bract pigmentation and will accelerate breeding in developing novel cultivars with richly colored bracts in C. alismatifolia and related species. It is also important to understand the variation in the evolution of the Zingiberaceae family.

Curcuma alismatifolia (Zingiberaceae) is an ornamental species with high economic value due to its recent rise in popularity among floriculturists. Cultivars within this species have mixed genetic backgrounds from multiple hybridization events and can be difficult to distinguish via morphological and histological methods alone. Given the need to improve identification resources, we carried out the first systematic study using plastomic data wherein genomic evolution and phylogenetic relationships from 56 accessions of C. alismatifolia were analyzed. The newly assembled plastomes were highly conserved and ranged from 162,139 bp to 164,111 bp, including 79 genes that code for proteins, 30 tRNA genes, and 4 rRNA genes. The A/T motif was the most common of SSRs in the assembled genomes. The Ka/Ks values of most genes were less than 1, and only two genes had Ka/Ks values above 1, which were rps15 (1.15), and ndhl (1.13) with petA equal to 1. The sequence divergence between different varieties of C. alismatifolia was large, and the percentage of variation in coding regions was lower than that in the non-coding regions. Such data will improve cultivar identification, marker assisted breeding, and preservation of germplasm resources.

There is a paradox in the plant mitochondrial genome, that is, the genic region evolves slowly while the intergenic region evolves rapidly. Thus, the intergenic regions of the plant mitochondrial genome are difficult to align across different species, even in closely related species. Here, to character the mechanism of this paradox, we identified interspecific variations in the Ginkgo biloba, Oryza sativa, and Arabidopsis thaliana mitochondrial and plastid genome at a genome-wide level. The substitution rate of synonymous sites in genic regions was similar to the substitution rate of intergenic regions, while the substitution rate of nonsynonymous sites in genic regions was lower than that in intergenic regions, suggesting the mutation inputs were the same among different categories within the organelle genome, but the selection pressure varied. The substitution rate of single-copy regions was higher than that of IR (inverted repeats) in the plastid genome at an intraspecific level. The substitution rate of single-copy regions was higher than that of repeats in the G. biloba and A. thaliana mitochondrial genomes, but lower in that of O. sativa. This difference may be related to the length and distribution of repeats. Copy number variations that existed in the G. biloba and O. sativa mitochondrial genomes were confirmed. This study reveals the intraspecific variation pattern of organelle genomes at a genome-wide level, and that copy number variations were common in plant mitochondrial genomes.