The Carboxy Terminus of YCF1 Contains a Motif Conserved throughout >500 Myr of Streptophyte Evolution

Assembly and comparative analysis of the first complete mitochondrial genome of Setaria italica

MAIN CONCLUSION

In this study, we assembled the first complete mitochondrial genome of Setaria italica and confirmed the multi-branched architecture. The foxtail millet (Setaria italica) holds significant agricultural importance, particularly in arid and semi-arid regions. It plays a pivotal role in diversifying dietary patterns and shaping planting strategies. Although the chloroplast genome of S. italica has been elucidated in recent studies, the complete mitochondrial genome remains largely unexplored. In this study, we employed PacBio HiFi sequencing platforms to sequence and assemble the complete mitochondrial genome. The mitochondrial genome spans a total length of 446,614 base pairs and harbors a comprehensive set of genetic elements, including 33 unique protein-coding genes (PCGs), encompassing 24 unique mitochondrial core genes and 9 variable genes, along with 20 transfer RNA (tRNA) genes and 3 ribosomal RNA (rRNA) genes. Our analysis of mitochondrial PCGs revealed a pronounced codon usage preference. For instance, the termination codon exhibits a marked preference for UAA, while alanine (Ala) exhibits a preference for GCU, and glutamine (Gln) favors CAA. Notably, the maximum Relative Synonymous Codon Usage (RSCU) values for cysteine (Cys) and phenylalanine (Phe) are both below 1.2, indicating a lack of strong codon usage preference for these amino acids. Phylogenetic analyses consistently place S. italica in close evolutionary proximity to Chrysopogon zizanioides, relative to other Panicoideae plants. Collinearity analysis showed that a total of 39 fragments were identified to display homology with both the mitochondrial and chloroplast genomes. A total of 417 potential RNA-editing sites were discovered across the 33 mitochondrial PCGs. Notably, all these editing events involved the conversion of cytosine (C) to uracil (U). Through the employment of PCR validation coupled with Sanger sequencing for the anticipated editing sites of these codons, RNA-editing events were conclusively identified at two specific loci: nad4L-2 and atp6-1030. The results of this study provide a pivotal foundation for advanced genomic breeding research in foxtail millet. Furthermore, they impart essential insights that will be instrumental for forthcoming investigations into the evolutionary and molecular dynamics of Panicoideae species.

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 Genomic Analysis Uncovers the Chloroplast Genome Variation and Phylogenetic Relationships of Camellia Species

Camellia is the largest genus in the family Theaceae. Due to phenotypic diversity, frequent hybridization, and polyploidization, an understanding of the phylogenetic relationships between Camellia species remains challenging. Comparative chloroplast (cp) genomics provides an informative resource for phylogenetic analyses of Camellia. In this study, 12 chloroplast genome sequences from nine Camellia species were determined using Illumina sequencing technology via de novo assembly. The cp genome sizes ranged from 156,545 to 157,021 bp and were organized into quadripartite regions with the typical angiosperm cp genomes. Each genome harbored 87 protein-coding, 37 transfer RNA, and 8 ribosomal RNA genes in the same order and orientation. Differences in long and short sequence repeats, SNPs, and InDels were detected across the 12 cp genomes. Combining with the complete cp sequences of seven other species in the genus Camellia, a total of nine intergenic sequence divergent hotspots and 14 protein-coding genes with high sequence polymorphism were identified. These hotspots, especially the InDel (~400 bp) located in atpH-atpI region, had sufficient potential to be used as barcode markers for further phylogenetic analysis and species identification. Principal component and phylogenetic analysis suggested that regional constraints, rather than functional constraints, strongly affected the sequence evolution of the cp genomes in this study. These cp genomes could facilitate the development of new molecular markers, accurate species identification, and investigations of the phylogenomic relationships of the genus Camellia.

The complete plastomes of red fleshed pitaya (Selenicereus monacanthus) and three related Selenicereus species: insights into gene losses, inverted repeat expansions and phylogenomic implications

Selenicereus is a genus of perennial shrub from the family Cactaceae, and some of them play an important role in the food industry, pharmaceuticals, cosmetics and medicine. To date, there are few reports on Selenicereus plastomes, which limits our understanding of this genus. Here, we have reported the complete plastomes of four Selenicereus species (S. monacanthus, S. annthonyanus, S. grandifloras, and S. validus) and carried out a comprehensive comparative analysis. All four Selenicereus plastomes have a typical quartile structure. The plastome size ranged from 133,146 to 134,450 bp, and contained 104 unique genes, including 30 tRNA genes, 4 rRNA genes and 70 protein-coding genes. Comparative analysis showed that there were massive losses of ndh genes in Selenicereus. Besides, we observed the inverted repeat regions had undergone a dramatic expansion and formed a previously unreported small single copy/inverted repeat border in the intron region of the atpF gene. Furthermore, we identified 6 hypervariable regions (trnF-GAA-rbcL, ycf1, accD, clpP-trnS-GCU, clpP-trnT-CGU and rpl22-rps19) that could be used as potential DNA barcodes for the identification of Selenicereus species. Our study enriches the plastome in the family Cactaceae, and provides the basis for the reconstruction of phylogenetic relationships.

SUPPLEMENTARY INFORMATION

The online version contains supplementary material available at 10.1007/s12298-021-01121-z.

Comparative analysis of chloroplast genomes in Vasconcellea pubescens A.DC. and Carica papaya L

The chloroplast genome is an integral part of plant genomes in a species along with nuclear and mitochondrial genomes, contributing to adaptation, diversification, and evolution of plant lineages. In the family Caricaceae, only the Carica papaya chloroplast genome and its nuclear and mitochondrial genomes were sequenced, and no chloroplast genome-wide comparison across genera was conducted. Here, we sequenced and assembled the chloroplast genome of Vasconcellea pubescens A.DC. using Oxford Nanopore Technology. The size of the genome is 158,712 bp, smaller than 160,100 bp of the C. papaya chloroplast genome. And two structural haplotypes, LSC_IRa_SSCrc_IRb and LSC_IRa_SSC_IRb, were identified in both V. pubescens and C. papaya chloroplast genomes. The insertion-deletion mutations may play an important role in Ycf1 gene evolution in family Caricaceae. Ycf2 is the only one gene positively selected in the V. pubescens chloroplast genome. In the C. papaya chloroplast genome, there are 46 RNA editing loci with an average RNA editing efficiency of 63%. These findings will improve our understanding of the genomes of these two crops in the family Caricaceae and will contribute to crop improvement.

Sequencing and Analysis of the Complete Organellar Genomes of Prototheca wickerhamii

Of the Prototheca genus, Prototheca wickerhamii has the highest clinical significance in humans. However, neither nuclear nor organellar genomes of this species were sequenced until now. The hitherto determined and analyzed mitochondrial and plastid genomes of the alleged P. wickerhamii species belong in fact to another species, recently named Prototheca xanthoriae. This study provides a first insight into the organellar genomes of a true P. wickerhamii (type strain ATCC 16529). The P. wickerhamii mitochondrion had a 53.8-kb genome, which was considerably larger than that of Prototheca ciferrii (formerly Prototheca zopfii gen. 1) and Prototheca bovis (formerly Prototheca zopfii gen. 2), yet similarly functional, with the differences in size attributable to a higher number of introns and the presence of extra unique putative genes. The 48-kb plastid genome of P. wickerhamii, compared to autotrophic Trebouxiophyceae, was highly reduced due to the elimination of the photosynthesis-related genes. The gene content of the plastid genome of P. wickerhamii was, however, very similar to other colorless Prototheca species. Plastid genome-based phylogeny reinforced the polyphyly of the genus Prototheca, with Helicosporidium and Auxenochlorella branching within clades of Prototheca species. Phylogenetic reconstruction also confirmed the close relationship of P. wickerhamii and P. xanthoriae, which is reflected in the synteny of their organellar genomes. Interestingly, the entire set of atp genes was lost in P. wickerhamii plastid genome while being preserved in P. xanthoriae.

Comparative chloroplast genome analysis of Artemisia (Asteraceae) in East Asia: insights into evolutionary divergence and phylogenomic implications

BACKGROUND

Artemisia in East Asia includes a number of economically important taxa that are widely used for food, medicinal, and ornamental purposes. The identification of taxa, however, has been hampered by insufficient diagnostic morphological characteristics and frequent natural hybridization. Development of novel DNA markers or barcodes with sufficient resolution to resolve taxonomic issues of Artemisia in East Asia is significant challenge.

RESULTS

To establish a molecular basis for taxonomic identification and comparative phylogenomic analysis of Artemisia, we newly determined 19 chloroplast genome (plastome) sequences of 18 Artemisia taxa in East Asia, de novo-assembled and annotated the plastomes of two taxa using publicly available Illumina reads, and compared them with 11 Artemisia plastomes reported previously. The plastomes of Artemisia were 150,858-151,318 base pairs (bp) in length and harbored 87 protein-coding genes, 37 transfer RNAs, and 8 ribosomal RNA genes in conserved order and orientation. Evolutionary analyses of whole plastomes and 80 non-redundant protein-coding genes revealed that the noncoding trnH-psbA spacer was highly variable in size and nucleotide sequence both between and within taxa, whereas the coding sequences of accD and ycf1 were under weak positive selection and relaxed selective constraints, respectively. Phylogenetic analysis of the whole plastomes based on maximum likelihood and Bayesian inference analyses yielded five groups of Artemisia plastomes clustered in the monophyletic subgenus Dracunculus and paraphyletic subgenus Artemisia, suggesting that the whole plastomes can be used as molecular markers to infer the chloroplast haplotypes of Artemisia taxa. Additionally, analysis of accD and ycf1 hotspots enabled the development of novel markers potentially applicable across the family Asteraceae with high discriminatory power.

CONCLUSIONS

The complete sequences of the Artemisia plastomes are sufficiently polymorphic to be used as super-barcodes for this genus. It will facilitate the development of new molecular markers and study of the phylogenomic relationships of Artemisia species in the family Asteraceae.

Major Changes in Plastid Protein Import and the Origin of the Chloroplastida

Core components of plastid protein import and the principle of using N-terminal targeting sequences are conserved across the Archaeplastida, but lineage-specific differences exist. Here we compare, in light of plastid protein import, the response to high-light stress from representatives of the three archaeplastidal groups. Similar to land plants, Chlamydomonas reinhardtii displays a broad response to high-light stress, not observed to the same degree in the glaucophyte Cyanophora paradoxa or the rhodophyte Porphyridium purpureum. We find that only the Chloroplastida encode both Toc75 and Oep80 in parallel and suggest that elaborate high-light stress response is supported by changes in plastid protein import. We propose the origin of a phenylalanine-independent import pathway via Toc75 allowed higher import rates to rapidly service high-light stress, but with the cost of reduced specificity. Changes in plastid protein import define the origin of the green lineage, whose greatest evolutionary success was arguably the colonization of land.

•

Chloroplastida evolved a dual system, Toc75/Oep80, for high throughput protein import

•

Loss of F-based targeting led to dual organelle targeting using a single ambiguous NTS

•

Relaxation of functional constraints allowed a wider Toc/Tic modification

•

A broad response to high-light stress appears unique to Chloroplastida

Large Diversity of Nonstandard Genes and Dynamic Evolution of Chloroplast Genomes in Siphonous Green Algae (Bryopsidales, Chlorophyta)

Chloroplast genomes have undergone tremendous alterations through the evolutionary history of the green algae (Chloroplastida). This study focuses on the evolution of chloroplast genomes in the siphonous green algae (order Bryopsidales). We present five new chloroplast genomes, which along with existing sequences, yield a data set representing all but one families of the order. Using comparative phylogenetic methods, we investigated the evolutionary dynamics of genomic features in the order. Our results show extensive variation in chloroplast genome architecture and intron content. Variation in genome size is accounted for by the amount of intergenic space and freestanding open reading frames that do not show significant homology to standard plastid genes. We show the diversity of these nonstandard genes based on their conserved protein domains, which are often associated with mobile functions (reverse transcriptase/intron maturase, integrases, phage- or plasmid-DNA primases, transposases, integrases, ligases). Investigation of the introns showed proliferation of group II introns in the early evolution of the order and their subsequent loss in the core Halimedineae, possibly through RT-mediated intron loss.

Embryophyte stress signaling evolved in the algal progenitors of land plants

The evolution of land plants from algae is an age-old question in biology. The entire terrestrial flora stems from a grade of algae, the streptophyte algae. Recent phylogenomic studies have pinpointed the Zygnematophyceae as the modern-day streptophyte algal lineage that is most closely related to the algal land plant ancestor. Here, we provide insight into the biology of this ancestor that might have aided in its conquest of land. Specifically, we uncover the existence of stress-signaling pathways and the potential for intimate plastid-nucleus communication. Plastids act as environmental sensors in land plants; our data suggest that this feature was present in a common ancestor they shared with streptophyte algae.

Streptophytes are unique among photosynthetic eukaryotes in having conquered land. As the ancestors of land plants, streptophyte algae are hypothesized to have possessed exaptations to the environmental stressors encountered during the transition to terrestrial life. Many of these stressors, including high irradiance and drought, are linked to plastid biology. We have investigated global gene expression patterns across all six major streptophyte algal lineages, analyzing a total of around 46,000 genes assembled from a little more than 1.64 billion sequence reads from six organisms under three growth conditions. Our results show that streptophyte algae respond to cold and high light stress via expression of hallmark genes used by land plants (embryophytes) during stress–response signaling and downstream responses. Among the strongest differentially regulated genes were those associated with plastid biology. We observed that among streptophyte algae, those most closely related to land plants, especially Zygnema, invest the largest fraction of their transcriptional budget in plastid-targeted proteins and possess an array of land plant-type plastid-nucleus communication genes. Streptophyte algae more closely related to land plants also appear most similar to land plants in their capacity to respond to plastid stressors. Support for this notion comes from the detection of a canonical abscisic acid receptor of the PYRABACTIN RESISTANCE (PYR/PYL/RCAR) family in Zygnema, the first found outside the land plant lineage. We conclude that a fine-tuned response toward terrestrial plastid stressors was among the exaptations that allowed streptophytes to colonize the terrestrial habitat on a global scale.

How Embryophytic is the Biosynthesis of Phenylpropanoids and their Derivatives in Streptophyte Algae?

The origin of land plants from algae is a long-standing question in evolutionary biology. It is becoming increasingly clear that many characters that were once assumed to be 'embryophyte specific' can in fact be found in their closest algal relatives, the streptophyte algae. One such case is the phenylpropanoid pathway. While biochemical data indicate that streptophyte algae harbor lignin-like components, the phenylpropanoid core pathway, which serves as the backbone of lignin biosynthesis, has been proposed to have arisen at the base of the land plants. Here we revisit this hypothesis using a wealth of new sequence data from streptophyte algae. Tracing the biochemical pathway towards lignin biogenesis, we show that most of the genes required for phenylpropanoid synthesis and the precursors for lignin production were already present in streptophyte algae. Nevertheless, phylogenetic analyses and protein structure predictions of one of the key enzyme classes in lignin production, cinnamyl alcohol dehydrogenase (CAD), suggest that CADs of streptophyte algae are more similar to sinapyl alcohol dehydrogenases (SADs). This suggests that the end-products of the pathway leading to lignin biosynthesis in streptophyte algae may facilitate the production of lignin-like compounds and defense molecules. We hypothesize that streptophyte algae already possessed the genetic toolkit from which the capacity to produce lignin later evolved in vascular plants.