Plant mitochondrial genomes (mitogenomes) exhibit extensive structural variation yet extremely low nucleotide mutation rates, phenomena that remain only partially understood. The genus Gossypium, a globally important source of cotton, offers a wealth of long-read sequencing resources to explore mitogenome and plastome variation and dynamics accompanying the evolutionary divergence of its approximately 50 diploid and allopolyploid species. Here, we assembled 19 mitogenomes from Gossypium species, representing all genome groups (diploids A through G, K, and the allopolyploids AD) based on a uniformly applied strategy. A graph-based mitogenome assembly method revealed more alternative structural conformations than previously recognized, some of which confirmed the mitogenome structure reported in earlier studies on cotton. Using long-read data, we quantified alternative conformations mediated by recombination events between repeats, and phylogenetically informative structural variants were noted. Nucleotide substitution rate comparisons between coding and non-coding regions revealed low mutation rates across the entire mitogenome. Genome-wide mapping of nuclear organellar DNA transfers (NUOTs) in Gossypium revealed a nonrandom distribution of transfers in the nuclear genome. In cotton, the fate of NUOT events varied, with mitochondrion-to-nucleus transfer (NUMT) predominantly retained as short fragments in the nuclear genome, with more plastid sequences integrated into the nucleus. Phylogenetic relationships inferred using different data sets highlighted distinct evolutionary histories among these cellular compartments, providing ancillary evidence relevant to the evolutionary history of Gossypium. A comprehensive analysis of organellar genome variation demonstrates complex structural variation and low mutation rates across the entire mitogenome and reveals the history of organellar genome transfer among the three genomes throughout the cotton genus. The findings enhance our general understanding of mitogenome evolution, comparative organellar and nuclear evolutionary rates, and the history of inter-compartment genomic integration.

Deep relationships in the angiosperm tree of life remain highly controversial. To address this, we first assembled the complete mitochondrial genomes for Ceratophyllum demersum and Chloranthus sessilifolius, confirming a well-supported sister relationship that starkly conflicts with nuclear and plastid data. To dissect this classic cyto-nuclear conflict, we developed the ‘PhyloForensics’ framework, a novel diagnostic approach to systematically identify sources of phylogenetic instability. This framework revealed that signal heterogeneity (topological entropy variance) and information content (the proportion of informative sites) are the primary drivers of gene-tree conflict. Empirically validating this, we show that removing a small subset of “loudly conflicted” genes resolves deep-level incongruence, yielding a single, highly-supported topology previously obscured by noise. Finally, complementing this sequence-based resolution, we demonstrate that mitogenome architecture provides powerful phylogenetic signals, revealing predictable, mitogenome-wide evolutionary patterns, such as a significant negative correlation between branch length and both GC content and RNA editing sites. By integrating a validated conflict-resolution framework with architectural genomics, our study provides a comprehensive strategy for navigating the complexities of deep evolutionary histories.

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.

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.