Comparative Analyses of 3,654 Plastid Genomes Unravel Insights Into Evolutionary Dynamics and Phylogenetic Discordance of Green Plants

Intraspecific differentiation of Lindera obtusiloba as revealed by comparative plastomic and evolutionary analyses

Lindera obtusiloba Blume is the northernmost tree species in the family Lauraceae, and it is a key species in understanding the evolutionary history of this family. The species of L. obtusiloba in East Asia has diverged into the Northern and Southern populations, which are geographically separated by an arid belt. Though the morphological differences between populations have been observed and well documented, intraspecific variations at the plastomic level have not been systematically investigated to date. Here, ten chloroplast genomes of L. obtusiloba individuals were sequenced and analyzed along with three publicly available plastomes. Comparative plastomic analysis suggests that both the Northern and the Southern populations share similar overall structure, gene order, and GC content in their plastomes although the size of the plasome and the level of intraspecific variability do vary between the two populations. The Northern have relatively larger plastomes while the Southern population possesses higher intraspecific variability, which could be attributed to the complexity of the geological environments in the South. Phylogenomic analyses also support the split of the Northern and Southern clades among L. obtusiloba individuals. However, there is no obvious species boundary between var. obtusiloba and var. heterophylla in the Southern population, indicating that gene flow could still occur between these two varieties, and this could be used as a good example of reticulate evolution. It is also found that a few photosynthesis-related genes are under positive selection, which is mainly related to the geological and environmental differences between the Northern and the Southern regions. Our results provide a reference for phylogenetic analysis within species and suggest that phylogenomic analyses with a sufficient number of nuclear and chloroplast genomic target loci from widely distributed individuals could provide a deeper understanding of the population evolution of the widespread species.

Recurrent evolutionary switches of mitochondrial cytochrome c maturation systems in Archaeplastida

Mitochondrial cytochrome c maturation (CCM) requires heme attachment via distinct pathways termed systems I and III. The mosaic distribution of these systems in Archaeplastida raises questions about the genetic mechanisms and evolutionary forces promoting repeated evolution. Here, we show a recurrent shift from ancestral system I to the eukaryotic-specific holocytochrome c synthase (HCCS) of system III in 11 archaeplastid lineages. Archaeplastid HCCS is sufficient to rescue mutants of yeast system III and Arabidopsis system I. Algal HCCS mutants exhibit impaired growth and respiration, and altered biochemical and metabolic profiles, likely resulting from deficient CCM and reduced cytochrome c-dependent respiratory activity. Our findings demonstrate that archaeplastid HCCS homologs function as system III components in the absence of system I. These results elucidate the evolutionary trajectory and functional divergence of CCM pathways in Archaeplastida, providing insight into the causes, mechanisms, and consequences of repeated cooption of an entire biological pathway.

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

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

Characterization and phylogenetic analysis of the complete plastid genome of Theobroma bicolor (Malvaceae) from Peru

Theobroma bicolor Bonpl. 1806 is distributed in the Neotropics from southern Mexico to the Peruvian and Brazilian Amazon. High-throughput sequencing of T. bicolor from Peru (KUELAP2926) resulted in the assembly of its complete plastid genome (GenBank accession number OQ557154). The chloroplast genome of T. bicolor is A + T-rich (62.97%), having 160,317 bp in size and containing 130 genes; including a pair of inverted repeat regions (IRs) of 25,462 bp separated by a large single copy region (LSC) of 89,221 bp and a small single copy region (SSC) of 20,172 bp. This plastid genome is similar in length, content, and organization to other members of the genus Theobroma. Phylogenetic analyses of T. bicolor support its sistership to the clade comprising T. cacao and T. grandiflorum. This study may contribute valuable information to the phylogenetic relationships within the genus Theobroma.

Comprehensive analysis of the chloroplast genome and phylogenetic relationships of Sasa quelpaertensis Nakai

Jeju-Joritdae (Sasa quelpaertensis Nakai) is a broad-leaved bamboo grass endemic to Mount Halla, Jeju Island, South Korea. In this study, we report the complete chloroplast genome sequence of S. quelpaertensis. Its chloroplast genome is 139,730 bp in size and consists of a large single-copy (LSC, 83,351 bp) region, one small single-copy (SSC, 12,788 bp) region, and two inverted repeats (IRs, 21,796 bp each). The chloroplast genome of S. quelpaertensis encodes 131 genes, including 86 protein-coding, 37 tRNA, and 8 rRNA genes. The overall GC content of the S. quelpaertensis chloroplast genome is 38.86%. Phylogenetic analysis using the chloroplast genome sequence showed that S. quelpaertensis is closely related to Sasa veitchii and Sasella kogasensis. These findings provide valuable genomic resources for future studies of the Sasa genus in South Korea and other countries encompassing its distribution area.

Characterizing conflict and congruence of molecular evolution across organellar genome sequences for phylogenetics in land plants

Chloroplasts and mitochondria each contain their own genomes, which have historically been and continue to be important sources of information for inferring the phylogenetic relationships among land plants. The organelles are predominantly inherited from the same parent, and therefore should exhibit phylogenetic concordance. In this study, we examine the mitochondrion and chloroplast genomes of 226 land plants to infer the degree of similarity between the organelles' evolutionary histories. Our results show largely concordant topologies are inferred between the organelles, aside from four well-supported conflicting relationships that warrant further investigation. Despite broad patterns of topological concordance, our findings suggest that the chloroplast and mitochondrial genomes evolved with significant differences in molecular evolution. The differences result in the genes from the chloroplast and the mitochondrion preferentially clustering with other genes from their respective organelles by a program that automates selection of evolutionary model partitions for sequence alignments. Further investigation showed that changes in compositional heterogeneity are not always uniform across divergences in the land plant tree of life. These results indicate that although the chloroplast and mitochondrial genomes have coexisted for over 1 billion years, phylogenetically, they are still evolving sufficiently independently to warrant separate models of evolution. As genome sequencing becomes more accessible, research into these organelles' evolution will continue revealing insight into the ancient cellular events that shaped not only their history, but the history of plants as a whole.

Genome-Wide Analysis of the MADS-Box Gene Family in Holoparasitic Plants (Balanophora subcupularis and Balanophora fungosa var. globosa)

MADS-box is an important transcription factor family that is involved in the regulation of various stages of plant growth and development, especially flowering regulation and flower development. Being a holoparasitic plant, the body structure of Balanophoraceae has changed dramatically over time, and its vegetative and reproductive organs have been extensively modified, with rudimentary flower organs. Meanwhile, extraordinary gene losses have been identified in holoparasitic plants compared with autotrophs. Our study reveals that the MADS-box gene family contracted sharply in Balanophora subcupularis and Balanophora fungosa var. globosa, and some subfamilies were lost, exhibiting reduced redundancy in both. The genes that functioned in the transition from the vegetative to floral production stages suffered a significant loss, but the ABCE model genes remained intact. We further investigated genes related to flowering regulation in B. subcupularis and B. fungosa var. globosa, vernalization and autonomous ways of regulating flowering time remained comparatively integrated, while genes in photoperiod and circadian clock pathways were almost lost. Convergent gene loss in flowering regulation occurred in Balanophora and another holoparasitic plant Sapria himalayana (Rafflesiaceae). The genome-wide analysis of the MADS-box gene family in Balanophora species provides valuable information for understanding the classification, gene loss pattern, and flowering regulation mechanism of MADS-box gene family in parasitic plants.