Exploring Plastomic Resources in Sempervivum (Crassulaceae): Implications for Phylogenetics

The plastid organelle is vital for photosynthesis and energy production. Advances in sequencing technology have enabled the exploration of plastomic resources, offering insights into plant evolution, diversity, and conservation. As an important group of horticultural ornamentals in the Crassulaceae family, Sempervivum plants are known for their unique rosette-like structures and reproduction through offsets. Despite their popularity, the classification status of Sempervivum remains uncertain, with only a single plastome sequence currently available. Furthermore, codon usage bias (CUB) is a widespread phenomenon of the unbalanced usage of synonymous codons in the coding sequence (CDS). However, due to the limited available plastid data, there has been no research that focused on the CUB analysis among Sempervivum until now. To address these gaps, we sequenced and released the plastomes of seven species and one subspecies from Sempervivum, revealing several consistent patterns. These included a shared 110 bp extension of the rps19 gene, 14 hypervariable regions (HVRs) with distinct nucleotide diversity (π: 0.01173 to 0.02702), and evidence of selective pressures shaping codon usage. Notably, phylogenetic analysis robustly divided the monophyletic clade into two sections: Jovibarba and Sempervivum. In conclusion, this comprehensive plastomic resource provides valuable insights into Sempervivum evolution and offers potential molecular markers for DNA barcoding.

Comparative plastome analysis of the sister genera Ceratocephala and Myosurus (Ranunculaceae) reveals signals of adaptive evolution to arid and aquatic environments

BACKGROUND

Expansion and contraction of inverted repeats can cause considerable variation of plastid genomes (plastomes) in angiosperms. However, little is known about whether structural variations of plastomes are associated with adaptation to or occupancy of new environments. Moreover, adaptive evolution of angiosperm plastid genes remains poorly understood. Here, we sequenced the complete plastomes for four species of xerophytic Ceratocephala and hydrophytic Myosurus, as well as Ficaria verna. By an integration of phylogenomic, comparative genomic, and selection pressure analyses, we investigated evolutionary patterns of plastomes in Ranunculeae and their relationships with adaptation to dry and aquatic habitats.

RESULTS

Owing to the significant contraction of the boundary of IRA/LSC towards the IRA, plastome sizes and IR lengths of Myosurus and Ceratocephala are smaller within Ranunculeae. Compared to other Ranunculeae, the Myosurus plastome lost clpP and rps16, one copy of rpl2 and rpl23, and one intron of rpoC1 and rpl16, and the Ceratocephala plastome added an infA gene and lost one copy of rpl2 and two introns of clpP. A total of 11 plastid genes (14%) showed positive selection, two genes common to Myosurus and Ceratocephala, seven in Ceratocephala only, and two in Myosurus only. Four genes showed strong signals of episodic positive selection. The rps7 gene of Ceratocephala and the rpl32 and ycf4 genes of Myosurus showed an increase in the rate of variation close to 3.3 Ma.

CONCLUSIONS

The plastomic structure variations as well as the positive selection of two plastid genes might be related to the colonization of new environments by the common ancestor of Ceratocephala and Myosurus. The seven and two genes under positive selection might be related to the adaptation to dry and aquatic habitats in Ceratocephala and Myosurus, respectively. Moreover, intensified aridity and frequent sea-level fluctuations, as well as global cooling, might have favored an increased rate of change in some genes at about 3.3 Ma, associated with adaptation to dry and aquatic environments, respectively. These findings suggest that changing environments might have influenced structural variations of plastomes and fixed new mutations arising on some plastid genes owing to adaptation to specific habitats.

Plastid genome of Passiflora tripartita var. mollissima (poro-poro) from Huánuco, Peru

Passiflora tripartita var. mollissima, known locally as poro-poro, is an important native fruit used in traditional Peruvian medicine with relevant agro-industrial and pharmaceutical potential for its antioxidant capacity for human health. However, to date, only a few genetic data are available, which limits exploring its genetic diversity and developing new genetic studies for its improvement. We report the poro-poro plastid genome to expand the knowledge of its molecular markers, evolutionary studies, molecular pathways, and conservation genetics. The complete chloroplast (cp) genome is 163,451 bp in length with a typical quadripartite structure, containing a large single-copy region of 85,525 bp and a small single-copy region of 13,518 bp, separated by a pair of inverted repeat regions (IR) of 32,204 bp, and the overall GC content was 36.87%. This cp genome contains 128 genes (110 genes were unique and 18 genes were found duplicated in each IR region), including 84 protein-coding genes, 36 transfer RNA-coding genes, eight ribosomal RNA-coding genes, and 13 genes with introns (11 genes with one intron and two genes with two introns). The inverted repeat region boundaries among species were similar in organization, gene order, and content, with a few revisions. The phylogenetic tree reconstructed based on single-copy orthologous genes and maximum likelihood analysis demonstrates poro-poro is most closely related to Passiflora menispermifolia and Passiflora oerstedii. In summary, our study constitutes a valuable resource for studying molecular evolution, phylogenetics, and domestication. It also provides a powerful foundation for conservation genetics research and plant breeding programs. To our knowledge, this is the first report on the plastid genome of Passiflora tripartita var. mollissima from Peru.

Plastome phylogenomics and morphological traits analyses provide new insights into the phylogenetic position, species delimitation and speciation of Triplostegia (Caprifoliaceae)

BACKGROUND

The genus Triplostegia contains two recognized species, T. glandulifera and T. grandiflora, but its phylogenetic position and species delimitation remain controversial. In this study, we assembled plastid genomes and nuclear ribosomal DNA (nrDNA) cistrons sampled from 22 wild Triplostegia individuals, each from a separate population, and examined these with 11 recently published Triplostegia plastomes. Morphological traits were measured from herbarium specimens and wild material, and ecological niche models were constructed.

RESULTS

Triplostegia is a monophyletic genus within the subfamily Dipsacoideae comprising three monophyletic species, T. glandulifera, T. grandiflora, and an unrecognized species Triplostegia sp. A, which occupies much higher altitude than the other two. The new species had previously been misidentified as T. glandulifera, but differs in taproot, leaf, and other characters. Triplotegia is an old genus, with stem age 39.96 Ma, and within it T. glandulifera diverged 7.94 Ma. Triplostegia grandiflora and sp. A diverged 1.05 Ma, perhaps in response to Quaternary climate fluctuations. Niche overlap between Triplostegia species was positively correlated with their phylogenetic relatedness.

CONCLUSIONS

Our results provide new insights into the species delimitation of Triplostegia, and indicate that a taxonomic revision of Triplostegia is needed. We also identified that either rpoB-trnC or ycf1 could serve as a DNA barcode for Triplostegia.

Complete Plastid Genomes of Nine Species of Ranunculeae (Ranunculaceae) and Their Phylogenetic Inferences

The tribe Ranunculeae, Ranunculaceae, comprising 19 genera widely distributed all over the world. Although a large number of Sanger sequencing-based molecular phylogenetic studies have been published, very few studies have been performed on using genomic data to infer phylogenetic relationships within Ranunculeae. In this study, the complete plastid genomes of nine species (eleven samples) from Ceratocephala, Halerpestes, and Ranunculus were de novo assembled using a next-generation sequencing method. Previously published plastomes of Oxygraphis and other related genera of the family were downloaded from GenBank for comparative analysis. The complete plastome of each Ranunculeae species has 112 genes in total, including 78 protein-coding genes, 30 transfer RNA genes, and four ribosomal RNA genes. The plastome structure of Ranunculeae samples is conserved in gene order and arrangement. There are no inverted repeat (IR) region expansions and only one IR contraction was found in the tested samples. This study also compared plastome sequences across all the samples in gene collinearity, codon usage, RNA editing sites, nucleotide variability, simple sequence repeats, and positive selection sites. Phylogeny of the available Ranunculeae species was inferred by the plastome data using maximum-likelihood and Bayesian inference methods, and data partitioning strategies were tested. The phylogenetic relationships were better resolved compared to previous studies based on Sanger sequencing methods, showing the potential value of the plastome data in inferring the phylogeny of the tribe.

More than a spiny morphology: plastome variation in the prickly pear cacti (Opuntieae)

BACKGROUND

Plastid genomes (plastomes) have long been recognized as highly conserved in their overall structure, size, gene arrangement and content among land plants. However, recent studies have shown that some lineages present unusual variations in some of these features. Members of the cactus family are one of these lineages, with distinct plastome structures reported across disparate lineages, including gene losses, inversions, boundary movements or loss of the canonical inverted repeat (IR) region. However, only a small fraction of cactus diversity has been analysed so far.

METHODS

Here, we investigated plastome features of the tribe Opuntieae, the remarkable prickly pear cacti, which represent one of the most diverse and important lineages of Cactaceae. We assembled de novo the plastome of 43 species, representing a comprehensive sampling of the tribe, including all seven genera, and analysed their evolution in a phylogenetic comparative framework. Phylogenomic analyses with different datasets (full plastome sequences and genes only) were performed, followed by congruence analyses to assess signals underlying contentious nodes.

KEY RESULTS

Plastomes varied considerably in length, from 121 to 162 kbp, with striking differences in the content and size of the IR region (contraction and expansion events), including a lack of the canonical IR in some lineages and the pseudogenization or loss of some genes. Overall, nine different types of plastomes were reported, deviating in the presence of the IR region or the genes contained in the IR. Overall, plastome sequences resolved phylogenetic relationships within major clades of Opuntieae with high bootstrap values but presented some contentious nodes depending on the dataset analysed (e.g. whole plastome vs. genes only). Congruence analyses revealed that most plastidial regions lack phylogenetic resolution, while few markers are supporting the most likely topology. Likewise, alternative topologies are driven by a handful of plastome markers, suggesting recalcitrant nodes in the phylogeny.

CONCLUSIONS

Our study reveals a dynamic nature of plastome evolution across closely related lineages, shedding light on peculiar features of plastomes. Variation of plastome types across Opuntieae is remarkable in size, structure and content and can be important for the recognition of species in some major clades. Unravelling connections between the causes of plastome variation and the consequences for species biology, physiology, ecology, diversification and adaptation is a promising and ambitious endeavour in cactus research. Although plastome data resolved major phylogenetic relationships, the generation of nuclear genomic data is necessary to confront these hypotheses and assess the recalcitrant nodes further.

Comparative Plastid Genome and Phylogenomic Analyses of Potamogeton Species

Potamogetonaceae are aquatic plants divided into six genera. The largest genus in the family is Potamogeton, which is morphologically diverse with many hybrids and polyploids. Potamogetonaceae plastomes were conserved in genome size (155,863 bp-156,669 bp), gene contents (113 genes in total, comprising 79 protein-coding genes and 30 tRNA and 4 rRNA genes), and GC content (36.5%). However, we detected a duplication of the trnH gene in the IR region of the Potamogeton crispus and P. maakianus plastomes. A comparative analysis of Alismatales indicated that the plastomes of Potamogetonaceae, Cymodaceae, and Ruppiaceae have experienced a 6-kb inversion of the rbcL-trnV region and the ndh complex has been lost in the Najas flexilis plastome. Five divergent hotspots (rps16-trnQ, atpF intron, rpoB-trnC, trnC-psbM, and ndhF-rpl32) were identified among the Potamogeton plastomes, which will be useful for species identification. Phylogenetic analyses showed that the family Potamogetonaceae is a well-defined with 100% bootstrap support and divided into two different clades, Potamogeton and Stuckenia. Compared to the nucleotide substitution rates among Alismatales, we found neutral selection in all plastid genes of Potamogeton species. Our results reveal the complete plastome sequences of Potamogeton species, and will be helpful for taxonomic identification, the elucidation of phylogenetic relationships, and the plastome structural analysis of aquatic plants.

Comprehensive Comparative Analyses of Aspidistra Chloroplast Genomes: Insights into Interspecific Plastid Diversity and Phylogeny

Limestone karsts are renowned for extremely high species richness and endemism. Aspidistra (Asparagaceae) is among the highly diversified genera distributed in karst areas, making it an ideal group for studying the evolutionary mechanisms of karst plants. The taxonomy and identification of Aspidistra species are mainly based on their specialized and diverse floral structures. Aspidistra plants have inconspicuous flowers, and the similarity in vegetative morphology often leads to difficulties in species discrimination. Chloroplast genomes possess variable genetic information and offer the potential for interspecies identification. However, as yet there is little information about the interspecific diversity and evolution of the plastid genomes of Aspidistra. In this study, we reported chloroplast (cp) genomes of seven Aspidistra species (A. crassifila, A. dolichanthera, A. erecta, A. longgangensis, A. minutiflora, A. nankunshanensis, and A. retusa). These seven highly-conserved plastid genomes all have a typical quartile structure and include a total of 113 unique genes, comprising 79 protein-coding genes, 4 rRNA genes and 30 tRNA genes. Additionally, we conducted a comprehensive comparative analysis of Aspidistra cp genomes. We identified eight divergent hotspot regions (trnC-GCA-petN, trnE-UUC-psbD, accD-psaI, petA-psbJ, rpl20-rps12, rpl36-rps8, ccsA-ndhD and rps15-ycf1) that serve as potential molecular markers. Our newly generated Aspidistra plastomes enrich the resources of plastid genomes of karst plants, and an investigation into the plastome diversity offers novel perspectives on the taxonomy, phylogeny and evolution of Aspidistra species.

Characterization and development of a plastid genome base editor, ptpTALECD

The modification of photosynthesis-related genes in plastid genomes may improve crop yields. Recently, we reported that a plastid-targeting base editor named ptpTALECD, in which a cytidine deaminase DddA functions as the catalytic domain, can homoplasmically substitute a targeted C to T in plastid genomes of Arabidopsis thaliana. However, some target Cs were not substituted. In addition, although ptpTALECD could substitute Cs on the 3' side of T and A, it was unclear whether it could also substitute Cs on the 3' side of G and C. In this study, we identified the preferential positions of the substituted Cs in ptpTALECD-targeting sequences in the Arabidopsis plastid genome. We also found that ptpTALECD could substitute Cs on the 3' side of all four bases in plastid genomes of Arabidopsis. More recently, a base editor containing an improved version of DddA (DddA11) was reported to substitute Cs more efficiently, and to substitute Cs on the 3' side of more varieties of bases in human mitochondrial genomes than a base editor containing DddA. Here, we also show that ptpTALECD_v2, in which a modified version of DddA11 functions as the catalytic domain, more frequently substituted Cs than ptpTALECD in the Arabidopsis plastid genome. We also found that ptpTALECD_v2 tended to substitute Cs at more positions than ptpTALECD. Our results reveal that ptpTALECD can cause a greater variety of codon changes and amino acid substitutions than previously thought, and that ptpTALECD and ptpTALECD_v2 are useful tools for the targeted base editing of plastid genomes.

Phylogenomics and plastome evolution of a Brazilian mycoheterotrophic orchid, Pogoniopsis schenckii

PREMISE

Pogoniopsis likely represents an independent photosynthesis loss in orchids. We use phylogenomic data to better identify the phylogenetic placement of this fully mycoheterotrophic taxon, and investigate its molecular evolution.

METHODS

We performed likelihood analysis of plastid and mitochondrial phylogenomic data to localize the position of Pogoniopsis schenckii in orchid phylogeny, and investigated the evolution of its plastid genome.

RESULTS

All analyses place Pogoniopsis in subfamily Epidendroideae, with strongest support from mitochondrial data, which also place it near tribe Sobralieae with moderately strong support. Extreme rate elevation in Pogoniopsis plastid genes broadly depresses branch support; in contrast, mitochondrial genes are only mildly rate elevated and display very modest and localized reductions in bootstrap support. Despite considerable genome reduction, including loss of photosynthesis genes and multiple translation apparatus genes, gene order in Pogoniopsis plastomes is identical to related autotrophs, apart from moderately shifted inverted repeat (IR) boundaries. All cis-spliced introns have been lost in retained genes. Two plastid genes (accD, rpl2) show significant strengthening of purifying selection. A retained plastid tRNA gene (trnE-UUC) of Pogoniopsis lacks an anticodon; we predict that it no longer functions in translation but retains a secondary role in heme biosynthesis.

CONCLUSIONS

Slowly evolving mitochondrial genes clarify the placement of Pogoniopsis in orchid phylogeny, a strong contrast with analysis of rate-elevated plastome data. We documented the effects of the novel loss of photosynthesis: for example, despite massive gene loss, its plastome is fully colinear with other orchids, and it displays only moderate shifts in selective pressure in retained genes.

Resolution, conflict and rate shifts: insights from a densely sampled plastome phylogeny for Rhododendron (Ericaceae)

BACKGROUND AND AIMS

Rhododendron is a species-rich and taxonomically challenging genus due to recent adaptive radiation and frequent hybridization. A well-resolved phylogenetic tree would help to understand the diverse history of Rhododendron in the Himalaya-Hengduan Mountains where the genus is most diverse.

METHODS

We reconstructed the phylogeny based on plastid genomes with broad taxon sampling, covering 161 species representing all eight subgenera and all 12 sections, including ~45 % of the Rhododendron species native to the Himalaya-Hengduan Mountains. We compared this phylogeny with nuclear phylogenies to elucidate reticulate evolutionary events and clarify relationships at all levels within the genus. We also estimated the timing and diversification history of Rhododendron, especially the two species-rich subgenera Rhododendron and Hymenanthes that comprise >90 % of Rhododendron species in the Himalaya-Hengduan Mountains.

KEY RESULTS

The full plastid dataset produced a well-resolved and supported phylogeny of Rhododendron. We identified 13 clades that were almost always monophyletic across all published phylogenies. The conflicts between nuclear and plastid phylogenies suggested strongly that reticulation events may have occurred in the deep lineage history of the genus. Within Rhododendron, subgenus Therorhodion diverged first at 56 Mya, then a burst of diversification occurred from 23.8 to 17.6 Mya, generating ten lineages among the component 12 clades of core Rhododendron. Diversification in subgenus Rhododendron accelerated c. 16.6 Mya and then became fairly continuous. Conversely, Hymenanthes diversification was slow at first, then accelerated very rapidly around 5 Mya. In the Himalaya-Hengduan Mountains, subgenus Rhododendron contained one major clade adapted to high altitudes and another to low altitudes, whereas most clades in Hymenanthes contained both low- and high-altitude species, indicating greater ecological plasticity during its diversification.

CONCLUSIONS

The 13 clades proposed here may help to identify specific ancient hybridization events. This study will help to establish a stable and reliable taxonomic framework for Rhododendron, and provides insight into what drove its diversification and ecological adaption. Denser sampling of taxa, examining both organelle and nuclear genomes, is needed to better understand the divergence and diversification history of Rhododendron.

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

BACKGROUND

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

RESULTS

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

CONCLUSIONS

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

Plastid phylogenomic insights into relationships, divergence, and evolution of Apiales

Members of Apiales are monophyletic and radiated in the Late Cretaceous. Fruit morphologies are critical for Apiales evolution and negative selection and mutation pressure play important roles in environmental adaptation. Apiales include many foods, spices, medicinal, and ornamental plants, but the phylogenetic relationships, origin and divergence, and adaptive evolution remain poorly understood. Here, we reconstructed Apiales phylogeny based on 72 plastid genes from 280 species plastid genomes representing six of seven families of this order. Highly supported phylogenetic relationships were detected, which revealed that each family of Apiales is monophyletic and confirmed that Pennanticeae is a member of Apiales. Genera Centella and Dickinsia are members of Apiaceae, and the genus Hydrocotyle previously classified into Apiaceae is confirmed to belong to Araliaceae. Besides, coalescent phylogenetic analysis and gene trees cluster revealed ten genes that can be used for distinguishing species among families of Apiales. Molecular dating suggested that the Apiales originated during the mid-Cretaceous (109.51 Ma), with the families' radiation occurring in the Late Cretaceous. Apiaceae species exhibit higher differentiation compared to other families. Ancestral trait reconstruction suggested that fruit morphological evolution may be related to shifts in plant types (herbaceous or woody), which in turn is related to the distribution areas and species numbers. Codon bias and positive selection analyses suggest that negative selection and mutation pressure may play important roles in environmental adaptation of Apiales members. Our results improve the phylogenetic framework of Apiales and provide insights into the origin, divergence, and adaptive evolution of this order and its members.

Phylogenomics, plastome degradation and mycoheterotrophy evolution of Neottieae (Orchidaceae), with emphasis on the systematic position and Loess Plateau-Changbai Mountains disjunction of Diplandrorchis

BACKGROUND

Mycoheterotrophy is a unique survival strategy adapted to dense forests and has attracted biologists' attention for centuries. However, its evolutionary origin and related plastome degradation are poorly understood. The tribe Neottieae contains various nutrition types, i.e., autotrophy, mixotrophy, and mycoheterotrophy. Here, we present a comprehensive phylogenetic analysis of the tribe based on plastome and nuclear ITS data. We inferred the evolutionary shift of nutrition types, constructed the patterns of plastome degradation, and estimated divergence times and ancestral ranges. We also used an integration of molecular dating and ecological niche modeling methods to investigate the disjunction between the Loess Plateau and Changbai Mountains in Diplandrorchis, a mycoheterotrophic genus endemic to China that was included in a molecular phylogenetic study for the first time.

RESULTS

Diplandrorchis was imbedded within Neottia and formed a clade with four mycoheterotrophic species. Autotrophy is the ancestral state in Neottieae, mixotrophy independently originated at least five times, and three shifts from mixotrophy to mycoheterotrophy independently occurred. The five mixotrophic lineages possess all plastid genes or lost partial/all ndh genes, whereas each of the three mycoheterotroph lineages has a highly reduced plastome: one lost part of its ndh genes and a few photosynthesis-related genes, and the other two lost almost all ndh, photosynthesis-related, rpo, and atp genes. These three mycoheterotrophic lineages originated at about 26.40 Ma, 25.84 Ma, and 9.22 Ma, respectively. Diplandrorchis had presumably a wide range in the Pliocene and migrated southward in the Pleistocene.

CONCLUSIONS

The Pleistocene climatic fluctuations and the resultant migration resulted in the Loess Plateau-Changbai Mountains disjunction of Diplandrorchis. In the evolution of mycoheterotrophic lineages, the loss of plastid-encoded genes and plastome degradation are staged and irreversible, constraining mycoheterotrophs to inhabit understories with low light levels. Accordingly, the rise of local forests might have promoted the origin of conditions in which mycoheterotrophy is advantageous.

Comparing complete organelle genomes of holoparasitic Christisonia kwangtungensis (Orabanchaceae) with its close relatives: how different are they?

BACKGROUND

Orobanchaceae is the only flowering plant family with species from free-living nonparasite, hemi-parasite to holoparasite, making it an ideal system for studying the evolution of parasitism. However, both plastid and mitochondrial genome have been sequenced in only few parasitic species in Orobanchaceae. Therefore, further comparative study is wanted to investigate the impact of holoparasitism on organelle genomes evolution between close relatives. Here, we sequenced organelle genomes and transcriptome of holoparasitic Christisonia kwangtungensis and compared it with its closely related groups to analyze similarities and differences in adaption strategies to the holoparasitic lifestyle.

RESULTS

The plastid genome of C. kwangtungensis has undergone extensive pseudogenization and gene loss, but its reduction pattern is different from that of Aeginetia indica, the close relative of C. kwangtungensis. Similarly, the gene expression detected in the photosynthetic pathway of these two genera is different. In Orobanchaceae, holoparasites in Buchnereae have more plastid gene loss than Rhinantheae, which reflects their longer history of holoparasitism. Distinct from severe degradation of the plastome, protein-coding genes in the mitochondrial genome of C. kwangtungensis are relatively conserved. Interestingly, besides intracellularly transferred genes which are still retained in its plastid genome, we also found several horizontally transferred genes of plastid origin from diverse donors other than their current hosts in the mitochondrial genome, which probably indicate historical hosts.

CONCLUSION

Even though C. kwangtungensis and A. indica are closely related and share severe degradation of plastome, they adapt organelle genomes to the parasitic lifestyle in different ways. The difference between their gene loss and gene expression shows they ultimately lost photosynthetic genes but through different pathways. Our study exemplifies how parasites part company after achieving holoparasitism.

Genome diploidization associates with cladogenesis, trait disparity, and plastid gene evolution

Angiosperm genome evolution was marked by many clade-specific whole-genome duplication events. The Microlepidieae is one of the monophyletic clades in the mustard family (Brassicaceae) formed after an ancient allotetraploidization. Postpolyploid cladogenesis has resulted in the extant c. 17 genera and 60 species endemic to Australia and New Zealand (10 species). As postpolyploid genome diploidization is a trial-and-error process under natural selection, it may proceed with different intensity and be associated with speciation events. In Microlepidieae, different extents of homoeologous recombination between the two parental subgenomes generated clades marked by slow ("cold") versus fast ("hot") genome diploidization. To gain a deeper understanding of postpolyploid genome evolution in Microlepidieae, we analyzed phylogenetic relationships in this tribe using complete chloroplast sequences, entire 35S rDNA units, and abundant repetitive sequences. The four recovered intra-tribal clades mirror the varied diploidization of Microlepidieae genomes, suggesting that the intrinsic genomic features underlying the extent of diploidization are shared among genera and species within one clade. Nevertheless, even congeneric species may exert considerable morphological disparity (e.g. in fruit shape), whereas some species within different clades experience extensive morphological convergence despite the different pace of their genome diploidization. We showed that faster genome diploidization is positively associated with mean morphological disparity and evolution of chloroplast genes (plastid-nuclear genome coevolution). Higher speciation rates in perennials than in annual species were observed. Altogether, our results confirm the potential of Microlepidieae as a promising subject for the analysis of postpolyploid genome diploidization in Brassicaceae.

Structural mutations of small single copy (SSC) region in the plastid genomes of five Cistanche species and inter-species identification

BACKGROUND

Cistanche is an important genus of Orobanchaceae, with critical medicinal, economic, and desertification control values. However, the phylogenetic relationships of Cistanche genus remained obscure. To date, no effective molecular markers have been reported to discriminate effectively the Cistanche closely related species reported here. In this study, we obtained and characterized the plastomes of four Cistanche species from China, to clarify the phylogenetic relationship within the genus, and to develop molecular markers for species discrimination.  RESULTS: Four Cistanche species (Cistanche deserticola, Cistanche salsa, Cistanche tubulosa and Cistanche sinensis), were deep-sequenced with Illumina. Their plastomes were assembled using SPAdes and annotated using CPGAVAS2. The plastic genomes were analyzed in detail, finding that all showed the conserved quadripartite structure (LSC-IR-SSC-IR) and with full sizes ranging from 75 to 111 Kbp. We observed a significant contraction of small single copy region (SSC, ranging from 0.4-29 Kbp) and expansion of inverted repeat region (IR, ranging from 6-30 Kbp), with C. deserticola and C. salsa showing the smallest SSCs with only one gene (rpl32). Compared with other Orobanchaceae species, Cistanche species showed extremely high rates of gene loss and pseudogenization, as reported for other parasitic Orobanchaceae species. Furthermore, analysis of sequence divergence on protein-coding genes showed the three genes (rpl22, clpP and ycf2) had undergone positive selection in the Cistanche species under study. In addition, by comparison of all available Cistanche plastomes we found 25 highly divergent intergenic spacer (IGS) regions that were used to predict two DNA barcode markers (Cis-mk01 and Cis-mk02 based on IGS region trnR-ACG-trnN-GUU) and eleven specific DNA barcode markers using Ecoprimer software. Experimental validation showed 100% species discrimination success rate with both type of markers.

CONCLUSION

Our findings have shown that Cistanche species are an ideal model to investigate the structure variation, gene loss and pseudogenization during the process of plastome evolution in parasitic species, providing new insights into the evolutionary relationships among the Cistanche species. In addition, the developed DNA barcodes markers allow the proper species identification, ensuring the effective and safe use of Cistanche species as medicinal products.

The evolution of extremely diverged plastomes in Selaginellaceae (lycophyte) is driven by repeat patterns and the underlying DNA maintenance machinery

Two factors are proposed to account for the unusual features of organellar genomes: the disruptions of organelle-targeted DNA replication, repair, and recombination (DNA-RRR) systems in the nuclear genome and repetitive elements in organellar genomes. Little is known about how these factors affect organellar genome evolution. The deep-branching vascular plant family Selaginellaceae is known to have a deficient DNA-RRR system and convergently evolved organellar genomes. However, we found that the plastid genome (plastome) of Selaginella sinensis has extremely accelerated substitution rates, a low GC content, pervasive repeat elements, a dynamic network structure, and it lacks direct or inverted repeats. Unexpectedly, its organelle DNA-RRR system is short of a plastid-targeted Recombinase A1 (RecA1) and a mitochondrion-targeted RecA3, in line with other explored Selaginella species. The plastome contains a large collection of short- and medium-sized repeats. Given the absence of RecA1 surveillance, we propose that these repeats trigger illegitimate recombination, accelerated mutation rates, and structural instability. The correlations between repeat quantity and architectural complexity in the Selaginella plastomes support these conclusions. We, therefore, hypothesize that the interplay of the deficient DNA-RRR system and the high repeat content has led to the extraordinary divergence of the S. sinensis plastome. Our study not only sheds new light on the mechanism of plastome divergence by emphasizing the power of cytonuclear integration, but it also reconciles the longstanding contradiction on the effects of DNA-RRR system disruption on genome structure evolution.

Comparative analysis of the complete plastid genomes in Prunus subgenus Cerasus (Rosaceae): Molecular structures and phylogenetic relationships

Prunus subgenus Cerasus (cherry) is an economically important group that distributed in temperate regions of the northern hemisphere. However, shared interspecific morphological traits and variability across taxa of Cerasus are among the impediments to taxonomic efforts to correctly delimit taxa. This is further complicated by a lack of genetic information on these taxa, with no focused genomic or phylogenetic studies being done on Cerasus. In this study, we conducted comparative analysis on the complete plastid genomes (plastomes) of 20 Cerasus species to gain a greater understanding of the attributes of the plastome of these taxa while helping resolve their phylogenetic placement in Prunus sensu lato and interspecific relationships within the subgenus. Our results displayed that (1) the plastomes of the 20 Cerasus species studied exhibited a typical quadripartite structure with conversed genome arrangement, structure, and moderate divergence. (2) The average size of complete plastomes for the Cerasus taxa studied was 157,861 bp, ranging from 157,458 to 158,024 bp. A total of 134 genes were annotated, including 86 protein-coding genes, 40 tRNAs, and 8 rRNAs across all species. In simple sequence repeat analysis, we found Cerasus had a comparable number of dispersed and tandem repeats to those identified in other angiosperm taxa, with only P. pseudocerasus found to contain trinucleotide repeats. Nucleotide diversity analysis revealed that the trnG-GCC gene and rpl32-trnL region had the highest Pi value showing potential as phylogenetic markers. (3) Two phylogenetic trees of the plastomes verified the monophyletic relationship of Cerasus and provided a more resolved species-level phylogeny. Our study provides detailed plastome information for exploring the phylogeny of subg. Cerasus taxa. We identified various types of repeats and nucleotide diversity hotspots, which can be a reference for species identification and reconstruction of phylogenetic relationships.

Towards the plastome evolution and phylogeny of Cycas L. (Cycadaceae): molecular-morphology discordance and gene tree space analysis

BACKGROUND

Plastid genomes (plastomes) present great potential in resolving multiscale phylogenetic relationship but few studies have focused on the influence of genetic characteristics of plastid genes, such as genetic variation and phylogenetic discordance, in resolving the phylogeny within a lineage. Here we examine plastome characteristics of Cycas L., the most diverse genus among extant cycads, and investigate the deep phylogenetic relationships within Cycas by sampling 47 plastomes representing all major clades from six sections.

RESULTS

All Cycas plastomes shared consistent gene content and structure with only one gene loss detected in Philippine species C. wadei. Three novel plastome regions (psbA-matK, trnN-ndhF, chlL-trnN) were identified as containing the highest nucleotide variability. Molecular evolutionary analysis showed most of the plastid protein-coding genes have been under purifying selection except ndhB. Phylogenomic analyses that alternatively included concatenated and coalescent methods, both identified four clades but with conflicting topologies at shallow nodes. Specifically, we found three species-rich Cycas sections, namely Stangerioides, Indosinenses and Cycas, were not or only weakly supported as monophyly based on plastomic phylogeny. Tree space analyses based on different tree-inference methods both revealed three gene clusters, of which the cluster with moderate genetic properties showed the best congruence with the favored phylogeny.

CONCLUSIONS

Our exploration in plastomic data for Cycas supports the idea that plastid protein-coding genes may exhibit discordance in phylogenetic signals. The incongruence between molecular phylogeny and morphological classification reported here may largely be attributed to the uniparental attribute of plastid, which cannot offer sufficient information to resolve the phylogeny. Contrasting to a previous consensus that genes with longer sequences and a higher proportion of variances are superior for phylogeny reconstruction, our result implies that the most effective phylogenetic signals could come from loci that own moderate variation, GC content, sequence length, and underwent modest selection.

The complete plastid genome sequence of the enigmatic moss, Takakia lepidozioides (Takakiopsida, Bryophyta): evolutionary perspectives on the largest collection of genes in mosses and the intensive RNA editing

Complete chloroplast genome sequence of a moss, Takakia lepidozioides (Takakiopsida) is reported. The largest collection of genes in mosses and the intensive RNA editing were discussed from evolutionary perspectives. We assembled the entire plastid genome sequence of Takakia lepidozioides (Takakiopsida), emerging from the first phylogenetic split among extant mosses. The genome sequences were assembled into a circular molecule 149,016 bp in length, with a quadripartite structure comprising a large and a small single-copy region separated by inverted repeats. It contained 88 genes coding for proteins, 32 for tRNA, four for rRNA, two open reading frames, and at least one pseudogene (tufA). This is the largest number of genes of all sequenced plastid genomes in mosses and Takakia is the only moss that retains the seven coding genes ccsA, cysA, cysT, petN rpoA, rps16 and trnP^GGG. Parsimonious interpretation of gene loss suggests that the last common ancestor of bryophytes had all seven genes and that mosses lost at least three of them during their diversification. Analyses of the plastid transcriptome identified the extraordinary frequency of RNA editing with more than 1100 sites. We indicated a close correlation between the monoplastidy of vegetative tissue and the intensive RNA editing sites in the plastid genome in land plant lineages. Here, we proposed a hypothesis that the small population size of plastids in each vegetative cell of some early diverging land plants, including Takakia, might cause the frequent fixation of mutations in plastid genome through the intracellular genetic drift and that deleterious mutations might be continuously compensated by RNA editing during or following transcription.

Nuclear-cytoplasmic balance: whole genome duplications induce elevated organellar genome copy number

The plant genome is partitioned across three distinct subcellular compartments: the nucleus, mitochondria, and plastids. Successful coordination of gene expression among these organellar genomes and the nuclear genome is critical for plant function and fitness. Whole genome duplication (WGD) events in the nucleus have played a major role in the diversification of land plants and are expected to perturb the relative copy number (stoichiometry) of nuclear, mitochondrial, and plastid genomes. Thus, elucidating the mechanisms whereby plant cells respond to the cytonuclear stoichiometric imbalance that follows WGDs represents an important yet underexplored question in understanding the evolutionary consequences of genome doubling. We used droplet digital PCR to investigate the relationship between nuclear and organellar genome copy numbers in allopolyploids and their diploid progenitors in both wheat and Arabidopsis. Polyploids exhibit elevated organellar genome copy numbers per cell, largely preserving the cytonuclear stoichiometry observed in diploids despite the change in nuclear genome copy number. To investigate the timescale over which cytonuclear stoichiometry may respond to WGD, we also estimated the organellar genome copy number in Arabidopsis synthetic autopolyploids and in a haploid-induced diploid line. We observed corresponding changes in organellar genome copy number in these laboratory-generated lines, indicating that at least some of the cellular response to cytonuclear stoichiometric imbalance is immediate following WGD. We conclude that increases in organellar genome copy numbers represent a common response to polyploidization, suggesting that maintenance of cytonuclear stoichiometry is an important component in establishing polyploid lineages.

Extensive genomic rearrangements mediated by repetitive sequences in plastomes of Medicago and its relatives

BACKGROUND

Although plastomes are highly conserved with respect to gene content and order in most photosynthetic angiosperms, extensive genomic rearrangements have been reported in Fabaceae, particularly within the inverted repeat lacking clade (IRLC) of Papilionoideae. Two hypotheses, i.e., the absence of the IR and the increased repeat content, have been proposed to affect the stability of plastomes. However, this is still unclear for the IRLC species. Here, we aimed to investigate the relationships between repeat content and the degree of genomic rearrangements in plastomes of Medicago and its relatives Trigonella and Melilotus, which are nested firmly within the IRLC.

RESULTS

We detected abundant repetitive elements and extensive genomic rearrangements in the 75 newly assembled plastomes of 20 species, including gene loss, intron loss and gain, pseudogenization, tRNA duplication, inversion, and a second independent IR gain (IR ~ 15 kb in Melilotus dentata) in addition to the previous first reported cases in Medicago minima. We also conducted comparative genomic analysis to evaluate plastome evolution. Our results indicated that the overall repeat content is positively correlated with the degree of genomic rearrangements. Some of the genomic rearrangements were found to be directly linked with repetitive sequences. Tandem repeated sequences have been detected in the three genes with accelerated substitution rates (i.e., accD, clpP, and ycf1) and their length variation could be explained by the insertions of tandem repeats. The repeat contents of the three localized hypermutation regions around these three genes with accelerated substitution rates are also significantly higher than that of the remaining plastome sequences.

CONCLUSIONS

Our results suggest that IR reemergence in the IRLC species does not ensure their plastome stability. Instead, repeat-mediated illegitimate recombination is the major mechanism leading to genome instability, a pattern in agreement with recent findings in other angiosperm lineages. The plastome data generated herein provide valuable genomic resources for further investigating the plastome evolution in legumes.

Plastid phylogenomics of Pleurothallidinae (Orchidaceae): Conservative plastomes, new variable markers, and comparative analyses of plastid, nuclear, and mitochondrial data

We present the first comparative plastome study of Pleurothallidinae with analyses of structural and molecular characteristics and identification of the ten most-variable regions to be incorporated in future phylogenetic studies. We sequenced complete plastomes of eight species in the subtribe and compared phylogenetic results of these to parallel analyses of their nuclear ribosomal DNA operon (26S, 18S, and 5.8S plus associated spacers) and partial mitochondrial genome sequences (29–38 genes and partial introns). These plastomes have the typical quadripartite structure for which gene content is similar to those of other orchids, with variation only in the composition of the ndh genes. The independent loss of ndh genes had an impact on which genes border the inverted repeats and thus the size of the small single-copy region, leading to variation in overall plastome length. Analyses of 68 coding sequences indicated the same pattern of codon usage as in other orchids, and 13 protein-coding genes under positive selection were detected. Also, we identified 62 polymorphic microsatellite loci and ten highly variable regions, for which we designed primers. Phylogenomic analyses showed that the top ten mutational hotspots represent well the phylogenetic relationships found with whole plastome sequences. However, strongly supported incongruence was observed among plastid, nuclear ribosomal DNA operon, and mitochondrial DNA trees, indicating possible occurrence of incomplete lineage sorting and/or introgressive hybridization. Despite the incongruence, the mtDNA tree retrieved some clades found in other analyses. These results, together with performance in recent studies, support a future role for mitochondrial markers in Pleurothallidinae phylogenetics.

airpg: automatically accessing the inverted repeats of archived plastid genomes

BACKGROUND

In most flowering plants, the plastid genome exhibits a quadripartite genome structure, comprising a large and a small single copy as well as two inverted repeat regions. Thousands of plastid genomes have been sequenced and submitted to public sequence repositories in recent years. The quality of sequence annotations in many of these submissions is known to be problematic, especially regarding annotations that specify the length and location of the inverted repeats: such annotations are either missing or portray the length or location of the repeats incorrectly. However, many biological investigations employ publicly available plastid genomes at face value and implicitly assume the correctness of their sequence annotations.

RESULTS

We introduce airpg, a Python package that automatically assesses the frequency of incomplete or incorrect annotations of the inverted repeats among publicly available plastid genomes. Specifically, the tool automatically retrieves plastid genomes from NCBI Nucleotide under variable search parameters, surveys them for length and location specifications of inverted repeats, and confirms any inverted repeat annotations through self-comparisons of the genome sequences. The package also includes functionality for automatic identification and removal of duplicate genome records and accounts for taxa that genuinely lack inverted repeats. A survey of the presence of inverted repeat annotations among all plastid genomes of flowering plants submitted to NCBI Nucleotide until the end of 2020 using airpg, followed by a statistical analysis of potential associations with record metadata, highlights that release year and publication status of the genome records have a significant effect on the frequency of complete and equal-length inverted repeat annotations.

CONCLUSION

The number of plastid genomes on NCBI Nucleotide has increased dramatically in recent years, and many more genomes will likely be submitted over the next decade. airpg enables researchers to automatically access and evaluate the inverted repeats of these plastid genomes as well as their sequence annotations and, thus, contributes to increasing the reliability of publicly available plastid genomes. The software is freely available via the Python package index at http://pypi.python.org/pypi/airpg .

Evidence of horizontal gene transfer between land plant plastids has surprising conservation implications

BACKGROUND AND AIMS

Horizontal gene transfer (HGT) is an important evolutionary mechanism because it transfers genetic material that may code for traits or functions between species or genomes. It is frequent in mitochondrial and nuclear genomes but has not been demonstrated between plastid genomes of different green land plant species.

METHODS

We Sanger-sequenced the nuclear internal transcribed spacers (ITS1 and 2) and the plastid rpl16 G2 intron (rpl16). In five individuals with foreign rpl16 we also sequenced atpB-rbcL and trnLUAA-trnFGAA.

KEY RESULTS

We discovered 14 individuals of a moss species with typical nuclear ITSs but foreign plastid rpl16 from a species of a distant lineage. None of the individuals with three plastid markers sequenced contained all foreign markers, demonstrating the transfer of plastid fragments rather than the entire plastid genome, i.e. entire plastids were not transferred. The two lineages diverged 165-185 Myr BP. The extended time interval since lineage divergence suggests that the foreign rpl16 is more likely explained by HGT than by hybridization or incomplete lineage sorting.

CONCLUSIONS

We provide the first conclusive evidence of interspecific plastid-to-plastid HGT among land plants. Two aspects are critical: it occurred at several localities during the massive colonization of recently disturbed open habitats that were created by large-scale liming as a freshwater biodiversity conservation measure; and it involved mosses whose unique life cycle includes spores that first develop a filamentous protonema phase. We hypothesize that gene transfer is facilitated when protonema filaments of different species intermix intimately when colonizing disturbed early succession habitats.

Plastid genome evolution in Amazonian açaí palm (Euterpe oleracea Mart.) and Atlantic forest açaí palm (Euterpe edulis Mart.)

The plastomes of E. edulis and E. oleracea revealed several molecular markers useful for genetic studies in natural populations and indicate specific evolutionary features determined by vicariant speciation. Arecaceae is a large and diverse family occurring in tropical and subtropical ecosystems worldwide. E. oleracea is a hyperdominant species of the Amazon forest, while E. edulis is a keystone species of the Atlantic forest. It has reported that E. edulis arose from vicariant speciation after the emergence of the belt barrier of dry environment (Cerrado and Caatinga biomes) between Amazon and Atlantic forests, isolating the E. edulis in the Atlantic forest. We sequenced the complete plastomes of E. edulis and E. oleracea and compared them concerning plastome structure, SSRs, tandem repeats, SNPs, indels, hotspots of nucleotide polymorphism, codon Ka/Ks ratios and RNA editing sites aiming to investigate evolutionary traits possibly affected by distinct environments. Our analyses revealed 303 SNPs, 91 indels, and 82 polymorphic SSRs among both species. Curiously, the narrow correlation among localization of repetitive sequences and indels strongly suggests that replication slippage is involved in plastid DNA mutations in Euterpe. Moreover, most non-synonymous substitutions represent amino acid variants in E. edulis that evolved specifically or in a convergent manner across the palm phylogeny. Amino acid variants observed in several plastid proteins in E. edulis were also identified as positive signatures across palm phylogeny. The higher incidence of specific amino acid changes in plastid genes of E. edulis in comparison with E. oleracea probably configures adaptive genetic variations determined by vicariant speciation. Our data indicate that the environment generates a selective pressure on the plastome making it more adapted to specific conditions.

Distinctive evolutionary pattern of organelle genomes linked to the nuclear genome in Selaginellaceae

Plastids and mitochondria are endosymbiotic organelles that store genetic information. The genomes of these organelles generally exhibit contrasting patterns regarding genome architecture and genetic content. However, they have similar genetic features in Selaginellaceae, and little is known about what causes parallel evolution. Here, we document the multipartite plastid genomes (plastomes) and the highly divergent mitochondrial genomes (mitogenomes) from spikemoss obtained by combining short- and long-reads. The 188-kb multipartite plastome has three ribosomal operon copies in the master genomic conformation, creating the alternative subgenomic conformation composed of 110- and 78-kb subgenomes. The long-read data indicated that the two different genomic conformations were present in almost equal proportions in the plastomes of Selaginella nipponica. The mitogenome of S. nipponica was assembled into 27 contigs with a total size of 110 kb. All contigs contained directly arranged repeats at both ends, which introduced multiple conformations. Our results showed that plastomes and mitogenomes share high tRNA losses, GC-biased nucleotides, elevated substitution rates and complicated organization. The exploration of nuclear-encoded organelle DNA replication, recombination and repair proteins indicated that, several single-targeted proteins, particularly plastid-targeted recombinase A1, have been lost in Selaginellaceae; conversely, the dual-targeted proteins remain intact. According to the reported function of recombinase A1, we propose that the plastomes of spikemoss often fail to pair homologous sequences during recombination, and the dual-targeted proteins play a key role in the convergent genetic features of plastomes and mitogenomes. Our results provide a distinctive evolutionary pattern of the organelle genomes in Selaginellaceae and evidence of their convergent evolution.

Plastome phylogenomic study of Gentianeae (Gentianaceae): widespread gene tree discordance and its association with evolutionary rate heterogeneity of plastid genes

BACKGROUND

Plastome-scale data have been prevalent in reconstructing the plant Tree of Life. However, phylogenomic studies currently based on plastomes rely primarily on maximum likelihood inference of concatenated alignments of plastid genes, and thus phylogenetic discordance produced by individual plastid genes has generally been ignored. Moreover, structural and functional characteristics of plastomes indicate that plastid genes may not evolve as a single locus and are experiencing different evolutionary forces, yet the genetic characteristics of plastid genes within a lineage remain poorly studied.

RESULTS

We sequenced and annotated 10 plastome sequences of Gentianeae. Phylogenomic analyses yielded robust relationships among genera within Gentianeae. We detected great variation of gene tree topologies and revealed that more than half of the genes, including one (atpB) of the three widely used plastid markers (rbcL, atpB and matK) in phylogenetic inference of Gentianeae, are likely contributing to phylogenetic ambiguity of Gentianeae. Estimation of nucleotide substitution rates showed extensive rate heterogeneity among different plastid genes and among different functional groups of genes. Comparative analysis suggested that the ribosomal protein (RPL and RPS) genes and the RNA polymerase (RPO) genes have higher substitution rates and genetic variations among plastid genes in Gentianeae. Our study revealed that just one (matK) of the three (matK, ndhB and rbcL) widely used markers show high phylogenetic informativeness (PI) value. Due to the high PI and lowest gene-tree discordance, rpoC2 is advocated as a promising plastid DNA barcode for taxonomic studies of Gentianeae. Furthermore, our analyses revealed a positive correlation of evolutionary rates with genetic variation of plastid genes, but a negative correlation with gene-tree discordance under purifying selection.

CONCLUSIONS

Overall, our results demonstrate the heterogeneity of nucleotide substitution rates and genetic characteristics among plastid genes providing new insights into plastome evolution, while highlighting the necessity of considering gene-tree discordance into phylogenomic studies based on plastome-scale data.

Phylogenetic and evolutionary features of the plastome of Tropaeolum pentaphyllum Lam. (Tropaeolaceae)

Complete plastome sequence of Tropaeolum pentaphyllum revealed molecular markers, hotspots of nucleotide polymorphism, RNA editing sites and phylogenetic aspects Tropaeolaceae Juss. ex DC. comprises approximately 95 species across North and South Americas. Tropaeolum pentaphyllum Lam. is an unconventional and endangered species with occurrence in some countries of South America. Although this species presents nutritional, medicinal and ornamental uses, genetic studies involving natural populations or promising genotypes are practically non-existent. Here, we report the nucleotide sequence of T. pentaphyllum plastome. It represents the first complete plastome sequence of the family Tropaeolaceae to be fully sequenced and analyzed in detail. The sequencing data revealed that the T. pentaphyllum plastome is highly similar to the plastomes of other Brassicales. Notwithstanding, our analyses detected some specific features concerning events of IR expansion and structural changes in some genes such as matK, rpoA, and rpoC2. We also detected 251 SSR loci, nine hotspots of nucleotide polymorphism, and two specific RNA editing sites in the plastome of T. pentaphyllum. Moreover, plastid phylogenomic inference indicated a closed relationship between the families Tropaeolaceae and Akaniaceae, which formed a sister group to Moringaceae-Caricaceae. Finally, our data bring new molecular markers and evolutionary features to be applied in the natural population, germplasm collection, and genotype selection aiming conservation, genetic diversity evaluation, and exploitation of this endangered species.

Degradation of key photosynthetic genes in the critically endangered semi-aquatic flowering plant Saniculiphyllum guangxiense (Saxifragaceae)

BACKGROUND

Plastid gene loss and pseudogenization has been widely documented in parasitic and mycoheterotrophic plants, which have relaxed selective constraints on photosynthetic function. More enigmatic are sporadic reports of pseudogenization and loss of important photosynthesis genes in lineages thought to be fully photosynthetic. Here we report the complete plastid genome of Saniculiphyllum guangxiense, a critically endangered and phylogenetically isolated plant lineage, along with genomic evidence of reduced chloroplast function. We also report 22 additional plastid genomes representing the diversity of its containing clade Saxifragales, characterizing gene content and placing variation in a broader phylogenetic context.

RESULTS

We find that the plastid genome of Saniculiphyllum has experienced pseudogenization of five genes of the ndh complex (ndhA, ndhB, ndhD, ndhF, and ndhK), previously reported in flowering plants with an aquatic habit, as well as the surprising pseudogenization of two genes more central to photosynthesis (ccsA and cemA), contrasting with strong phylogenetic conservatism of plastid gene content in all other sampled Saxifragales. These genes participate in photooxidative protection, cytochrome synthesis, and carbon uptake. Nuclear paralogs exist for all seven plastid pseudogenes, yet these are also unlikely to be functional.

CONCLUSIONS

Saniculiphyllum appears to represent the greatest degree of plastid gene loss observed to date in any fully photosynthetic lineage, perhaps related to its extreme habitat specialization, yet plastid genome length, structure, and substitution rate are within the variation previously reported for photosynthetic plants. These results highlight the increasingly appreciated dynamism of plastid genomes, otherwise highly conserved across a billion years of green plant evolution, in plants with highly specialized life history traits.

Exploration of Plastid Phylogenomic Conflict Yields New Insights into the Deep Relationships of Leguminosae

Phylogenomic analyses have helped resolve many recalcitrant relationships in the angiosperm tree of life, yet phylogenetic resolution of the backbone of the Leguminosae, one of the largest and most economically and ecologically important families, remains poor due to generally limited molecular data and incomplete taxon sampling of previous studies. Here, we resolve many of the Leguminosae's thorniest nodes through comprehensive analysis of plastome-scale data using multiple modified coding and noncoding data sets of 187 species representing almost all major clades of the family. Additionally, we thoroughly characterize conflicting phylogenomic signal across the plastome in light of the family's complex history of plastome evolution. Most analyses produced largely congruent topologies with strong statistical support and provided strong support for resolution of some long-controversial deep relationships among the early diverging lineages of the subfamilies Caesalpinioideae and Papilionoideae. The robust phylogenetic backbone reconstructed in this study establishes a framework for future studies on legume classification, evolution, and diversification. However, conflicting phylogenetic signal was detected and quantified at several key nodes that prevent the confident resolution of these nodes using plastome data alone. [Leguminosae; maximum likelihood; phylogenetic conflict; plastome; recalcitrant relationships; stochasticity; systematic error.].

Plastid-targeted forms of restriction endonucleases enhance the plastid genome rearrangement rate and trigger the reorganization of its genomic architecture

Plant cells have acquired chloroplasts (plastids) with a unique genome (ptDNA), which developed during the evolution of endosymbiosis. The gene content and genome structure of ptDNAs in land plants are considerably stable, although those of algal ptDNAs are highly varied. Plant cells seem, therefore, to be intolerant of any structural or organizational changes in the ptDNA. Genome rearrangement functions as a driver of genomic evolutionary divergence. Here, we aimed to create various types of rearrangements in the ptDNA of Arabidopsis genomes using plastid-targeted forms of restriction endonucleases (pREs). Arabidopsis plants expressing each of the three specific pREs, i.e., pTaqI, pHinP1I, and pMseI, were generated; they showed the leaf variegation phenotypes associated with impaired chloroplast development. We confirmed that these pREs caused double-stranded breaks (DSB) at their recognition sites in ptDNAs. Genome-wide analysis of ptDNAs revealed that the transgenic lines exhibited a large number of rearrangements such as inversions and deletions/duplications, which were dominantly repaired by microhomology-mediated recombination and microhomology-mediated end-joining, and less by non-homologous end-joining. Notably, pHinP1I, which recognized a small number of sites in ptDNA, induced drastic structural changes, including regional copy number variations throughout ptDNAs. In contrast, the transient expression of either pTaqI or pMseI, whose recognition site numbers were relatively larger, resulted in small-scale changes at the whole genome level. These results indicated that DSB frequencies and their distribution are major determinants in shaping ptDNAs.

Leveraging plastomes for comparative analysis and phylogenomic inference within Scutellarioideae (Lamiaceae)

Scutellaria, or skullcaps, are medicinally important herbs in China, India, Japan, and elsewhere. Though Scutellaria is the second largest and one of the more taxonomically challenging genera within Lamiaceae, few molecular systematic studies have been undertaken within the genus; in part due to a paucity of available informative markers. The lack of informative molecular markers for Scutellaria hinders our ability to accurately and robustly reconstruct phylogenetic relationships, which hampers our understanding of the diversity, phylogeny, and evolutionary history of this cosmopolitan genus. Comparative analyses of 15 plastomes, representing 14 species of subfamily Scutellarioideae, indicate that plastomes within Scutellarioideae contain about 151,000 nucleotides, and possess a typical quadripartite structure. In total, 590 simple sequence repeats, 489 longer repeats, and 16 hyper-variable regions were identified from the 15 plastomes. Phylogenetic relationships among the 14 species representing four of the five genera of Scutellarioideae were resolved with high support values, but the current infrageneric classification of Scutellaria was not supported in all analyses. Complete plastome sequences provide better resolution at an interspecific level than using few to several plastid markers in phylogenetic reconstruction. The data presented here will serve as a foundation to facilitate DNA barcoding, species identification, and systematic research within Scutellaria, which is an important medicinal plant resource worldwide.

Chloroplast nucleoids are highly dynamic in ploidy, number, and structure during angiosperm leaf development

Chloroplast nucleoids are large, compact nucleoprotein structures containing multiple copies of the plastid genome. Studies on structural and quantitative changes of plastid DNA (ptDNA) during leaf development are scarce and have produced controversial data. We have systematically investigated nucleoid dynamics and ptDNA quantities in the mesophyll of Arabidopsis, tobacco, sugar beet, and maize from the early post-meristematic stage until necrosis. DNA of individual nucleoids was quantified by DAPI-based supersensitive epifluorescence microscopy. Nucleoids occurred in scattered, stacked, or ring-shaped arrangements and in recurring patterns during leaf development that was remarkably similar between the species studied. Nucleoids per organelle varied from a few in meristematic plastids to >30 in mature chloroplasts (corresponding to about 20-750 nucleoids per cell). Nucleoid ploidies ranged from haploid to >20-fold even within individual organelles, with average values between 2.6-fold and 6.7-fold and little changes during leaf development. DNA quantities per organelle increased gradually from about a dozen plastome copies in tiny plastids of apex cells to 70-130 copies in chloroplasts of about 7 μm diameter in mature mesophyll tissue, and from about 80 plastome copies in meristematic cells to 2600-3300 copies in mature diploid mesophyll cells without conspicuous decline during leaf development. Pulsed-field electrophoresis, restriction of high-molecular-weight DNA from chloroplasts and gerontoplasts, and CsCl equilibrium centrifugation of single-stranded and double-stranded ptDNA revealed no noticeable fragmentation of the organelle DNA during leaf development, implying that plastid genomes in mesophyll tissues are remarkably stable until senescence.

Comparison of the abilities of universal, super, and specific DNA barcodes to discriminate among the original species of Fritillariae cirrhosae bulbus and its adulterants

Fritillariae cirrhosae bulbus is a famous type of traditional Chinese medicine used for cough relief and eliminating phlegm. The medicine originates from dried bulbs of five species and one variety of Fritillaria. Recently, immature bulbs from other congeneric species, such as F. ussuriensis, have been sold as adulterants of Fritillariae cirrhosae bulbus in medicine markets owing to the high price and limited availability of the genuine medicine. However, it is difficult to accurately identify the bulbs from different original species of Fritillariae cirrhosae bulbus and its adulterants based on traditional methods, although such medicines have different prices and treatment efficacies. The present study adopted DNA barcoding to identify these different species and compared the discriminatory power of super, universal, and specific barcodes in Fritillaria. The results revealed that the super-barcode had strong discriminatory power (87.5%). Among universal barcodes, matK provided the best species resolution (87.5%), followed by ITS (62.5%), rbcL (62.5%), and trnH-psbA (25%). The combination of these four universal barcodes provided the highest discriminatory power (87.5%), which was equivalent to that of the super-barcode. Two plastid genes, ycf1 and psbM-psbD, had much better discriminatory power (both 87.5%) than did other plastid barcodes, and were suggested as potential specific barcodes for identifying Fritillaria species. Phylogenetic analyses indicated that F. cirrhosa was not a "good" species that was composed of multiple lineages, which might have affected the evaluation of the discriminatory ability. This study revealed that the complete plastid genome, as super barcode, was an efficient and reliable tool for identifying the original species of Fritillariae cirrhosae bulbus and its adulterants.

Plastid genomes of Elaeagnus mollis: comparative and phylogenetic analyses

Plastomes, which are maternally inherited and show a moderate rate of evolution, play a critical role in phylogenetic reconstruction and assignment of plant species. However, little is known about the sequence divergence and molecular evolutionary patterns of plastid genomes in Elaeagnus mollis, a plant of great economic, medicinal, edible and ecological values. The plastid genome of E. mollis is 152,224-bp long and has 47 repeat sequences, including tandem (17), dispersed (12), and palindromic (18) types of repeat variations. Here, we reported six divergence hotspots (atpH-atpI, petN-psbM, trnT-psbD, trnP-psaJ, rpl32-trnL and ycf1) that could potentially be used as molecular genetic markers for population genetics and phylogenetic studies of E. mollis. A comparison of plastid genomes in the order Rosales showed that the trnH gene was duplicated only in Elaeagnaceae; therefore, it is an important marker in Elaeagnaceae. Phylogenetic analyses based on whole plastid genome sequences in 33 species revealed that Rosales is divided into two strongly supported clades and that the families Elaeagnaceae and Barbeyaceae are closely related.

Variation in plastid genomes in the gynodioecious species Silene vulgaris

BACKGROUND

Gynodioecious species exist in two sexes - male-sterile females and hermaphrodites. Male sterility in higher plants often results from mitonuclear interaction between the CMS (cytoplasmic male sterility) gene(s) encoded by mitochondrial genome and by nuclear-encoded restorer genes. Mitochondrial and nuclear-encoded transcriptomes in females and hermaphrodites are intensively studied, but little is known about sex-specific gene expression in plastids. We have compared plastid transcriptomes between females and hermaphrodites in two haplotypes of a gynodioecious species Silene vulgaris with known CMS candidate genes.

RESULTS

We generated complete plastid genome sequences from five haplotypes S. vulgaris including the haplotypes KRA and KOV, for which complete mitochondrial genome sequences were already published. We constructed a phylogenetic tree based on plastid sequences of S. vulgaris. Whereas lowland S. vulgaris haplotypes including KRA and KOV clustered together, the accessions from high European mountains diverged early in the phylogram. S. vulgaris belongs among Silene species with slowly evolving plastid genomes, but we still detected 212 substitutions and 112 indels between two accessions of this species. We estimated elevated Ka/Ks in the ndhF gene, which may reflect the adaptation of S. vulgaris to high altitudes, or relaxed selection. We compared depth of coverage and editing rates between female and hermaphrodite plastid transcriptomes and found no significant differences between the two sexes. We identified 51 unique C to U editing sites in the plastid genomes of S. vulgaris, 38 of them in protein coding regions, 2 in introns, and 11 in intergenic regions. The editing site in the psbZ gene was edited only in one of two plastid genomes under study.

CONCLUSIONS

We revealed no significant differences between the sexes in plastid transcriptomes of two haplotypes of S. vulgaris. It suggests that gene expression of plastid genes is not affected by CMS in flower buds of S. vulgaris, although both sexes may still differ in plastid gene expression in specific tissues. We revealed the difference between the plastid transcriptomes of two S. vulgaris haplotypes in editing rate and in the coverage of several antisense transcripts. Our results document the variation in plastid genomes and transcriptomes in S. vulgaris.

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

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

OrganellarGenomeDRAW (OGDRAW) version 1.3.1: expanded toolkit for the graphical visualization of organellar genomes

Organellar (plastid and mitochondrial) genomes play an important role in resolving phylogenetic relationships, and next-generation sequencing technologies have led to a burst in their availability. The ongoing massive sequencing efforts require software tools for routine assembly and annotation of organellar genomes as well as their display as physical maps. OrganellarGenomeDRAW (OGDRAW) has become the standard tool to draw graphical maps of plastid and mitochondrial genomes. Here, we present a new version of OGDRAW equipped with a new front end. Besides several new features, OGDRAW now has access to a local copy of the organelle genome database of the NCBI RefSeq project. Together with batch processing of (multi-)GenBank files, this enables the user to easily visualize large sets of organellar genomes spanning entire taxonomic clades. The new OGDRAW server can be accessed at https://chlorobox.mpimp-golm.mpg.de/OGDraw.html.

Directed Repeats Co-occur with Few Short-Dispersed Repeats in Plastid Genome of a Spikemoss, Selaginella vardei (Selaginellaceae, Lycopodiopsida)

BACKGROUND

It is hypothesized that the highly conserved inverted repeats (IR) structure of land plant plastid genomes (plastomes) is beneficial for stabilizing plastome organization, whereas the mechanism of the occurrence and stability maintenance of the recently reported direct repeats (DR) structure is yet awaiting further exploration. Here we describe the DR structure of the Selaginella vardei (Selaginellaceae) plastome, to elucidate the mechanism of DR occurrence and stability maintenance.

RESULTS

The plastome of S. vardei is 121,254 bp in length and encodes 76 genes, of which 62 encode proteins, 10 encode tRNAs, and four encode rRNAs. Unexpectedly, the two identical rRNA gene regions (13,893 bp) are arranged in a direct orientation (DR), rather than inverted. Comparing to the IR organization in Isoetes flaccida (Isoetaceae, Lycopodiopsida) plastome, a ca. 50-kb trnN-trnF inversion that spans one DR copy was found in the plastome of S. vardei, which might cause the orientation change. In addition, we find extremely rare short dispersed repeats (SDRs) in the plastomes of S. vardei and its closely related species S. indica.

CONCLUSIONS

We suggest that the ca. 50-kb inversion resulted in the DR structure, and the reduction in SDRs plays a key role in maintaining the stability of plastomes with DR structure by avoiding potential secondary recombination. We further confirmed the presence of homologous recombination between DR regions, which are able to generate subgenomes and form diverse multimers. Our study deepens the understanding of Selaginella plastomes and provides new insights into the diverse plastome structures in land plants.

Targeting mutations to the plastidial psbA gene of Chlamydomonas reinhardtii without direct positive selection

Targeting mutations to specific genomic loci is invaluable for assessing in vivo the effect of these changes on the biological role of the gene in study. Here, we attempted to introduce a mutation that was previously implicated in an increased heat stability of the mesophilic cyanobacterium Synechocystis sp. PCC6803 via homologous recombination to the psbA gene of Chlamydomonas reinhardtii. For that, we established a strategy for targeted mutagenesis that was derived from the efficient genome-wide homologous-recombination-based methodology that was used to target individual genes of Saccharomyces cerevisiae. While the isolated mutants did not show any benefit under elevated temperature conditions, the new strategy proved to be efficient for C. reinhardtii even in the absence of direct positive selection.

Cytonuclear interactions remain stable during allopolyploid evolution despite repeated whole-genome duplications in Brassica

Several plastid macromolecular protein complexes are encoded by both nuclear and plastid genes. Therefore, cytonuclear interactions are held in place to prevent genomic conflicts that may lead to incompatibilities. Allopolyploidy resulting from hybridization and genome doubling of two divergent species can disrupt these fine-tuned interactions, as newly formed allopolyploid species confront biparental nuclear chromosomes with a uniparentally inherited plastid genome. To avoid any deleterious effects of unequal genome inheritance, preferential transcription of the plastid donor over the other donor has been hypothesized to occur in allopolyploids. We used Brassica as a model to study the effects of paleopolyploidy in diploid parental species, as well as the effects of recent and ancient allopolyploidy in Brassica napus, on genes implicated in plastid protein complexes. We first identified redundant nuclear copies involved in those complexes. Compared with cytosolic protein complexes and with genome-wide retention rates, genes involved in plastid protein complexes show a higher retention of genes in duplicated and triplicated copies. Those redundant copies are functional and are undergoing strong purifying selection. We then compared transcription patterns and sequences of those redundant gene copies between resynthesized allopolyploids and their diploid parents. The neopolyploids showed no biased subgenome expression or maternal homogenization via gene conversion, despite the presence of some non-synonymous substitutions between plastid genomes of parental progenitors. Instead, subgenome dominance was observed regardless of the maternal progenitor. Our results provide new insights on the evolution of plastid protein complexes that could be tested and generalized in other allopolyploid species.

Recombination and intraspecific polymorphism for the presence and absence of entire chromosomes in mitochondrial genomes

Although mitochondrial genomes are typically thought of as single circular molecules, these genomes are fragmented into multiple chromosomes in many eukaryotes, raising intriguing questions about inheritance and (in)stability of mtDNA in such systems. A previous comparison of mitochondrial genomes from two different individuals of the angiosperm species Silene noctiflora found variation in the presence of entire mitochondrial chromosomes. Here, we expand on this work with a geographically diverse sampling of 25 S. noctiflora populations and the closely related species S. turkestanica and S. undulata. Using a combination of deep sequencing and PCR-based screening for the presence of 22 different mitochondrial chromosomes, we found extensive variation in the complement of chromosomes across individuals. Much of this variation could be attributed to recent chromosome loss events, suggesting that the massively expanded and fragmented mitochondrial genomes of S. noctiflora may have entered a phase of genome reduction in which they are losing entire chromosomes at a rapid rate. Sequence analysis of mitochondrial and plastid genomes revealed genealogical differences both between these organelles and within the mitochondrial genome, indicating a history of recombination. Evidence that recombination has generated novel combinations of alleles was more frequent between loci on different mitochondrial chromosomes than it was within chromosomes. Therefore, the fragmentation of mitochondrial genomes and the assortment of chromosomes during mitochondrial inheritance appears to have contributed to a history of sexual-like recombination in the mtDNA of this species.

On the postglacial spread of human commensal Arabidopsis thaliana: journey to the East

With more sequenced genomes, our understanding of the demographic history of Arabidopsis thaliana is rapidly expanding. However, no-one has yet compiled previous data to investigate patterns of genetic variation across Eurasia. While sub-Saharan accessions have been reported to be the most divergent group, in the nuclear genome we found accessions from Yunnan, China to be genetically closest to the sub-Saharan group. In chloroplast, several deeply diverged haplogroups exist only in Eurasia, and African populations have lower variation in many haplogroups that they share with the Eurasian populations. These patterns cannot be easily explained by a single out-of-Africa event suggested previously. For more recent demographic history, we dated the nonrelict expansion to 10 ka. In the Chinese Yangtze nonrelicts, we found clear traces of gene flow with local relicts, and genes under strong selection were enriched for traces of relict introgression, especially those related to biotic and immune responses. The results suggest the ability of nonrelicts to obtain locally adaptive alleles through admixture with relicts is important for the expansion across environmental gradients of Eurasia. Our re-analyses provide another model for the early history as well as elucidating factors contributing to the recent demographic turnover event of this species.

Chloroplast competition is controlled by lipid biosynthesis in evening primroses

In most eukaryotes, organellar genomes are transmitted preferentially by the mother, but molecular mechanisms and evolutionary forces underlying this fundamental biological principle are far from understood. It is believed that biparental inheritance promotes competition between the cytoplasmic organelles and allows the spread of so-called selfish cytoplasmic elements. Those can be, for example, fast-replicating or aggressive chloroplasts (plastids) that are incompatible with the hybrid nuclear genome and therefore maladaptive. Here we show that the ability of plastids to compete against each other is a metabolic phenotype determined by extremely rapidly evolving genes in the plastid genome of the evening primrose Oenothera Repeats in the regulatory region of accD (the plastid-encoded subunit of the acetyl-CoA carboxylase, which catalyzes the first and rate-limiting step of lipid biosynthesis), as well as in ycf2 (a giant reading frame of still unknown function), are responsible for the differences in competitive behavior of plastid genotypes. Polymorphisms in these genes influence lipid synthesis and most likely profiles of the plastid envelope membrane. These in turn determine plastid division and/or turnover rates and hence competitiveness. This work uncovers cytoplasmic drive loci controlling the outcome of biparental chloroplast transmission. Here, they define the mode of chloroplast inheritance, as plastid competitiveness can result in uniparental inheritance (through elimination of the "weak" plastid) or biparental inheritance (when two similarly "strong" plastids are transmitted).

Phylogeny, character evolution and spatiotemporal diversification of the species-rich and world-wide distributed tribe Rubieae (Rubiaceae)

The Rubiaceae tribe Rubieae has a world-wide distribution with up to 1,000 species. These collectively exhibit an enormous ecological and morphological diversity, making Rubieae an excellent group for macro- and microevolutionary studies. Previous molecular phylogenetic analyses used only a limited sampling within the tribe or missed lineages crucial for understanding character evolution in this group. Here, we analyze sequences from two plastid spacer regions as well as morphological and biogeographic data from an extensive and evenly distributed sampling to establish a sound phylogenetic framework. This framework serves as a basis for our investigation of the evolution of important morphological characters and the biogeographic history of the Rubieae. The tribe includes three major clades, the Kelloggiinae Clade (Kelloggia), the Rubiinae Clade (Didymaea, Rubia) and the most species-rich Galiinae Clade (Asperula, Callipeltis, Crucianella, Cruciata, Galium, Mericarpaea, Phuopsis, Sherardia, Valantia). Within the Galiinae Clade, the largest genera Galium and Asperula are para- and polyphyletic, respectively. Smaller clades, however, usually correspond to currently recognized taxa (small genera or sections within genera), which may be used as starting points for a refined classification in this clade. Life-form (perennial versus annual), flower shape (long versus short corolla tube) and fruit characters (dry versus fleshy, with or without uncinate hairs) are highly homoplasious and have changed multiple times independently. Inference on the evolution of leaf whorls, a characteristic feature of the tribe, is sensitive to model choice. Multi-parted leaf whorls appear to have originated from opposite leaves with two small interpetiolar stipules that are subsequently enlarged and increased in number. Early diversification of Rubieae probably started during the Miocene in western Eurasia. Disjunctions between the Old and the New World possibly are due to connections via a North Atlantic land bridge. Diversification of the Galiineae Clade started later in the Miocene, probably in the Mediterranean, from where lineages reached, often multiple times, Africa, eastern Asia and further on the Americas and Australia.

Efficient Replication of the Plastid Genome Requires an Organellar Thymidine Kinase

Thymidine kinase (TK) is a key enzyme of the salvage pathway that recycles thymidine nucleosides to produce deoxythymidine triphosphate. Here, we identified the single TK of maize (Zea mays), denoted CPTK1, as necessary in the replication of the plastidial genome (cpDNA), demonstrating the essential function of the salvage pathway during chloroplast biogenesis. CPTK1 localized to both plastids and mitochondria, and its absence resulted in an albino phenotype, reduced cpDNA copy number and a severe deficiency in plastidial ribosomes. Mitochondria were not affected, indicating they are less reliant on the salvage pathway. Arabidopsis (Arabidopsis thaliana) TKs, TK1A and TK1B, apparently resulted from a gene duplication after the divergence of monocots and dicots. Similar but less-severe effects were observed for Arabidopsis tk1a tk1b double mutants in comparison to those in maize cptk1 TK1B was important for cpDNA replication and repair in conditions of replicative stress but had little impact on the mitochondrial phenotype. In the maize cptk1 mutant, the DNA from the small single-copy region of the plastidial genome was reduced to a greater extent than other regions, suggesting preferential abortion of replication in this region. This was accompanied by the accumulation of truncated genomes that resulted, at least in part, from unfaithful microhomology-mediated repair. These and other results suggest that the loss of normal cpDNA replication elicits the mobilization of new replication origins around the rpoB (beta subunit of plastid-encoded RNA polymerase) transcription unit and imply that increased transcription at rpoB is associated with the initiation of cpDNA replication.

Dynamic evolution of inverted repeats in Euglenophyta plastid genomes

Photosynthetic euglenids (Euglenophyta) are a monophyletic group of unicellular eukaryotes characterized by the presence of plastids, which arose as the result of the secondary endosymbiosis. Many Euglenophyta plastid (pt) genomes have been characterized recently, but they represented mainly one family - Euglenaceae. Here, we report a comparative analysis of plastid genomes from eight representatives of the family Phacaceae. Newly sequenced plastid genomes share a number of features including synteny and gene content, except for genes mat2 and mat5 encoding maturases. The observed diversity of intron number and presence/absence of maturases corroborated previously suggested correlation between the number of maturases in the pt genome and intron proliferation. Surprisingly, pt genomes of taxa belonging to Discoplastis and Lepocinclis encode two inverted repeat (IR) regions containing the rDNA operon, which are absent from the Euglenaceae. By mapping the presence/absence of IR region on the obtained phylogenomic tree, we reconstructed the most probable events in the evolution of IRs in the Euglenophyta. Our study highlights the dynamic nature of the Euglenophyta plastid genome, in particular with regards to the IR regions that underwent losses repeatedly.