Gene Loss and Evolution of the Plastome

Insights into U-to-C RNA editing from the lycophyte Phylloglossum drummondii

The lycophyte Phylloglossum drummondii is the sole inhabitant of its genus in the Huperzioideae group and one of a small minority of plants which perform uridine to cytidine RNA editing. We assembled the P. drummondii chloroplast and mitochondrial genomes and used RNA sequence data to build a comprehensive profile of organellar RNA editing events. In addition to many C-to-U editing events in both organelles, we found just four U-to-C editing events in the mitochondrial transcripts cob, nad1, nad5 and rpl2. These events are conserved in related lycophytes in the genera Huperzia and Phlegmariurus. De novo transcriptomes for three of these lycophytes were assembled to search for putative U-to-C RNA editing enzymes. Four putative U-to-C editing factors could be matched to the four mitochondrial U-to-C editing sites. Due to the unusually few numbers of U-to-C RNA editing sites, P. drummondii and related lycophytes are useful models for studying this poorly understood mechanism.

Codon usage bias and phylogenetic analysis of chloroplast genome in 36 gracilariaceae species

Gracilariaceae is a group of marine large red algae and main source of agar with important economic and ecological value. The codon usage patterns of chloroplast genomes in 36 species from Graciliaceae show that GC range from 0.284 to 0.335, the average GC3 range from 0.135 to 0.243 and the value of ENC range from 35.098 to 42.327, which indicates these genomes are rich in AT and prefer to use codons ending with AT in these species. Nc plot, PR2 plot, neutrality plot analyses and correlation analysis indicate that these biases may be caused by multiple factors, such as natural selection and mutation pressure, but prolonged natural selection is the main driving force influencing codon usage preference. The cluster analysis and phylogenetic analysis show that the differentiation relationship of them is different and indicate that codons with weak or unbiased preferences may also play an irreplaceable role in these species' evolution. In addition, we identified 26 common high-frequency codons and 8-18 optimal codons all ending in A/U in these 36 species. Our results will not only contribute to carrying out transgenic work in Gracilariaceae species to maximize the protein yield in the future, but also lay a theoretical foundation for further exploring systematic classification of them.

Plastome diversity and evolution in mosses: Insights from structural characterization, comparative genomics, and phylogenetic analysis

Mosses play a significant role in ecology, evolution, and the economy. They belong to the nonvascular plant kingdom and are considered the closest living relatives of the first terrestrial plants. The circular chloroplast DNA molecules (plastomes) of mosses contain all the genetic information essential for chloroplast functions and represent the source of the evolutionary history of these organisms. This study comprehensively analyzed the plastomes of 47 moss species belonging to 14 orders, focusing on their size, GC content, gene loss, gene content, synteny, and evolution. The findings revealed great differences among plastome sizes, with Takakia lepidozioides (Takakiopsida) and Funaria hygrometrica (Funariales) having the largest and smallest plastomes, respectively. Moss plastomes included 69 to 89 protein-coding genes, 8 rRNA genes, and 34 to 42 tRNA genes, resulting in the total number of genes in a plastome ranging between 115 and 138. Various genes have been lost from the plastomes of different moss species, with Atrichum angustatum lacking the highest number of genes. This study also examined plastome synteny and moss evolution using comparative genomics and repeat sequence analysis. The results demonstrated that synteny and similarity levels varied across the 47 moss examined species, with some exhibiting structure similarity and others displaying structural inversions. Maximum likelihood and Bayesian approaches were used to construct a phylogenetic tree using 36 concatenated protein-coding genes, and the results revealed that the genera Sphagnum and Takakia are sister groups to the other mosses. Additionally, it was found that Tetraphidales, Polytrichales, Buxbaumiales, and Diphysciales are closely related. This research describes the evolutionary diversity of mosses and offers guidelines for future studies in this field. The findings also highlight the need for more investigations into the factors regulating plastome size variation in these plants.

Stochastic organelle genome segregation through Arabidopsis development and reproduction

Organelle DNA (oDNA) in mitochondria and plastids is vital for plant (and eukaryotic) life. Selection against damaged oDNA is mediated in part by segregation - sorting different oDNA types into different cells in the germline. Plants segregate oDNA very rapidly, with oDNA recombination protein MSH1 a key driver of this segregation, but we have limited knowledge of the dynamics of this segregation within plants and between generations. Here, we reveal how oDNA evolves through Arabidopsis thaliana development and reproduction. We combine stochastic modelling, Bayesian inference, and model selection with new and existing tissue-specific oDNA measurements from heteroplasmic Arabidopsis plant lines through development and between generations. Segregation proceeds gradually but continually during plant development, with a more rapid increase between inflorescence formation and the next generation. When MSH1 is compromised, the majority of observed segregation can be achieved through partitioning at cell divisions. When MSH1 is functional, mtDNA segregation is far more rapid; we show that increased oDNA gene conversion is a plausible mechanism quantitatively explaining this acceleration. These findings reveal the quantitative, time-dependent details of oDNA segregation in Arabidopsis. We also discuss the support for different models of the plant germline provided by these observations.

The complete chloroplast genome of Goodyera yunnanensis Schltr

The family Orchidaceae is renowned for its extensive diversity. Within this family, the genus Goodyera R. Br. is classified under the subtribe Goodyerinae, comprising approximately 99 species. In this study, a species Goodyera yunnanensis Schltr., its plastid genome was characterized. The plastid genome of G. yunnanensis is 146,197 bp in size and exhibits a typical quadripartite structure with a pair of inverted repeat regions (IRs) of 25,611 bp, a large single-copy region (LSC) of 81,300 bp and a small single-copy region (SSC) of 13,675 bp. A total of 126 genes were identified, containing 80 protein-coding genes, 38 tRNA genes and 8 rRNA genes. The overall GC content is 37.2%, with corresponding values of 43.3%, 34.7% and 29.1% in IR, LSC and SSC regions, respectively. Forty-seven simple sequence repeats (SSRs) are found in G. yunnanensis plastome, and the frequency of mononucleotide repeats is significantly higher than other repeat types. Phylogenetic analysis indicates that Goodyera is resolved into four clades. G. yunnanensis belongs to the monophyletic clade A, and its phylogenetic position can be reasonably supported by morphological and molecular data.

Complete Chloroplast Genomes of Pterodon emarginatus Vogel and Pterodon pubescens Benth: Comparative and Phylogenetic Analyses

Background

The species Pterodon emarginatus and P. pubescens, popularly known as white sucupira or faveira, are native to the Cerrado biome and have the potential for medicinal use and reforestation. They are sister species with evolutionary proximity.

Objective

Considering that the chloroplast genome exhibits a conserved structure and genes, the analysis of its sequences can contribute to the understanding of evolutionary, phylogenetic, and diversity issues.

Methods

The chloroplast genomes of P. emarginatus and P. pubescens were sequenced on the Illumina MiSeq platform. The genomes were assembled based on the de novo strategy. We performed the annotation of the genes and the repetitive regions of the genomes. The nucleotide diversity and phylogenetic relationships were analyzed using the gene sequences of these species and others of the Leguminosae family, whose genomes are available in databases.

Results

The complete chloroplast genome of P. emarginatus is 159,877 bp, and that of P. pubescens is 159,873 bp. The genomes of both species have circular and quadripartite structures. A total of 127 genes were predicted in both species, including 110 single-copy genes and 17 duplicated genes in the inverted regions. 141 microsatellite regions were identified in P. emarginatus and 140 in P. pubescens. The nucleotide diversity estimates of the gene regions in twenty-one species of the Leguminosae family were 0.062 in LSC, 0.086 in SSC, and 0.036 in IR. The phylogenetic analysis demonstrated the proximity between the genera Pterodon and Dipteryx, both from the clade Dipterygeae. Ten pairs of primers with potential for the development of molecular markers were designed.

Conclusion

The genetic information obtained on the chloroplast genomes of P. emarginatus and P. pubescens presented here reinforces the similarity and evolutionary proximity between these species, with a similarity percentage of 99.8%.

New Insights Into The Evolution of Chloroplast Genomes in Ochna Species (Ochnaceae, Malpighiales)

Ochnaceae DC. includes more than 600 species that exhibit potential values for environmental ecology, ornamental, pharmaceutical, and timber industries. Although studies on phylogeny and phytochemicals have been intensively conducted, chloroplast genome data of Ochnaceae species have not been fully explored. In this study, the next-generation sequencing method was used to sequence the chloroplast genomes of Ochna integerrima and Ochna serrulata which were 157 329 and 157 835 bp in length, respectively. These chloroplast genomes had a quadripartite structure and contained 78 protein-coding genes, 30 tRNAs, and 4 rRNAs. Comparative analysis revealed 8 hypervariable regions, including trnK_UUU-trnQ_UUG, rpoB-psbM, trnS_GGA-rps4, accD-psaI, rpl33-rps18, rpl14-rpl16, ndhF-trnL_UAG, and rps15-ycf1 among 6 Ochnaceae taxa. Additionally, there were shared and unique repeats among 6 examined chloroplast genomes. The notable changes were the loss of rpl32 in Ochna species and the deletion of rps16 exon 2 in O. integerrima compared to other taxa. This study is the first comprehensive comparative genomic analysis of complete chloroplast genomes of Ochna species and related taxa in Ochnaceae. Consequently, the current study provides initial results for further research on genomic evolution, population genetics, and developing molecular markers in Ochnaceae and related taxa.

Comparative analysis of the chloroplast genomes of Rosa species and RNA editing analysis

BACKGROUND

The genus Rosa (Rosaceae) contains approximately 200 species, most of which have high ecological and economic values. Chloroplast genome sequences are important for studying species differentiation, phylogeny, and RNA editing.

RESULTS

In this study, the chloroplast genomes of three Rosa species, Rosa hybrida, Rosa acicularis, and Rosa rubiginosa, were assembled and compared with other reported Rosa chloroplast genomes. To investigate the RNA editing sites in R. hybrida (commercial rose cultivar), we mapped RNA-sequencing data to the chloroplast genome and analyzed their post-transcriptional features. Rosa chloroplast genomes presented a quadripartite structure and had highly conserved gene order and gene content. We identified four mutation hotspots (ycf3-trnS, trnT-trnL, psbE-petL, and ycf1) as candidate molecular markers for differentiation in the Rosa species. Additionally, 22 chloroplast genomic fragments with a total length of 6,192 bp and > 90% sequence similarity with their counterparts were identified in the mitochondrial genome, representing 3.96% of the chloroplast genome. Phylogenetic analysis including all sections and all subgenera revealed that the earliest divergence in the chloroplast phylogeny roughly distinguished species of sections Pimpinellifoliae and Rosa and subgenera Hulthemia. Moreover, DNA- and RNA-sequencing data revealed 19 RNA editing sites, including three synonymous and 16 nonsynonymous, in the chloroplast genome of R. hybrida that were distributed among 13 genes.

CONCLUSIONS

The genome structure and gene content of Rosa chloroplast genomes are similar across various species. Phylogenetic analysis based on the Rosa chloroplast genomes has high resolution. Additionally, a total of 19 RNA editing sites were validated by RNA-Seq mapping in R. hybrida. The results provide valuable information for RNA editing and evolutionary studies of Rosa and a basis for further studies on genomic breeding of Rosa species.

Chloroplast genome characteristics and phylogeny of the sinodielsia clade (apiaceae: apioideae)

BACKGROUND

The Sinodielsia clade of the subfamily Apioideae (Apiacieae) was established in 2008, and it is composed of 37 species from 17 genera. Its circumscription is still poorly delimited and unstable, and interspecific relationships in the clade lack comprehensive analysis. Chloroplast (cp.) genomes provide valuable and informative data sources for evolutionary biology and have been widely used in studies on plant phylogeny. To infer the phylogenetic history of the Sinodielsia clade, we assembled complete cp. genomes of 39 species and then performed phylogenetic analysis based on these cp. genome sequence data combined with 66 published cp. genomes from 16 genera relative to the Sinodielsia clade.

RESULTS

These 39 newly assembled genomes had a typical quadripartite structure with two inverted repeat regions (IRs: 17,599-31,486 bp) separated by a large single-copy region (LSC: 82,048-94,046 bp) and a small single-copy region (SSC: 16,343-17,917 bp). The phylogenetic analysis showed that 19 species were clustered into the Sinodielsia clade, and they were divided into two subclades. Six mutation hotspot regions were detected from the whole cp. genomes among the Sinodielsia clade, namely, rbcL-accD, ycf4-cemA, petA-psbJ, ycf1-ndhF, ndhF-rpl32 and ycf1, and it was found that ndhF-rpl32 and ycf1 were highly variable in the 105 sampled cp. genomes.

CONCLUSION

The Sinodielsia clade was subdivided into two subclades relevant to geographical distributions, except for cultivated and introduced species. Six mutation hotspot regions, especially ndhF-rpl32 and ycf1, could be used as potential DNA markers in the identification and phylogenetic analyses of the Sinodielsia clade and Apioideae. Our study provided new insights into the phylogeny of the Sinodielsia clade and valuable information on cp. genome evolution in Apioideae.

Chloroplast genome characterization of Uncaria guianensis and Uncaria tomentosa and evolutive dynamics of the Cinchonoideae subfamily

Uncaria species are used in traditional medicine and are considered of high therapeutic value and economic importance. This work describes the assembly and annotation of the chloroplast genomes of U. guianensis and U. tomentosa, as well as a comparative analysis. The genomes were sequenced on MiSeq Illumina, assembled with NovoPlasty, and annotated using CHLOROBOX GeSeq. Addictionaly, comparative analysis were performed with six species from NCBI databases and primers were designed in Primer3 for hypervariable regions based on the consensus sequence of 16 species of the Rubiaceae family and validated on an in-silico PCR in OpenPrimeR. The genome size of U. guianensis and U. tomentosa was 155,505 bp and 156,390 bp, respectively. Both Species have 131 genes and GC content of 37.50%. The regions rpl32-ccsA, ycf1, and ndhF-ccsA showed the three highest values of nucleotide diversity within the species of the Rubiaceae family and within the Uncaria genus, these regions were trnH-psbA, psbM-trnY, and rps16-psbK. Our results indicates that the primer of the region ndhA had an amplification success for all species tested and can be promising for usage in the Rubiaceae family. The phylogenetic analysis recovered a congruent topology to APG IV. The gene content and the chloroplast genome structure of the analyzed species are conserved and most of the genes are under negative selection. We provide the cpDNA of Neotropical Uncaria species, an important genomic resource for evolutionary studies of the group.

Evolution of 101 Apocynaceae plastomes and phylogenetic implications

Apocynaceae are one of the ten species-richest angiosperm families. However, the backbone phylogeny of the family is yet less well supported, and the evolution of plastome structure has not been thoroughly studied for the whole family. Herein, a total of 101 complete plastomes including 35 newly sequenced, 24 reassembled from public raw data and the rest from the NCBI GenBank database, representing 26 of 27 tribes of Apocynaceae, were used for comparative plastome analysis. Phylogenetic analyses were conducted using a combined plastid data matrix of 77 protein-coding genes from 162 taxa, encompassing all tribes and 41 of 49 subtribes of Apocynaceae. Plastome lengths ranged from 150,897 bp in Apocynum venetum to 178,616 bp in Hoya exilis. Six types of boundaries between the inverted repeat (IR) regions and single copy (SC) regions were identified. Different sizes of IR expansion were found in three lineages, including Alyxieae, Ceropegieae and Marsdenieae, suggesting multiple expansion events of the IRs over the SC regions in Apocynaceae. The IR regions of Marsdenieae evolved in two ways: expansion towards the large single copy (LSC) region in Lygisma + Stephanotis + Ruehssia + Gymnema (Cosmopolitan clade), and expansion towards both LSC and small single copy (SSC) region in Dischidia-Hoya alliance and Marsdenia (Asia-Pacific clade). Six coding genes and five non-coding regions were identified as highly variable, including accD, ccsA-ndhD, clpP, matK, ndhF, ndhG-ndhI, trnG(GCC)-trnfM(CAU), trnH(GUG)-psbA, trnY(GUA)-trnE(UUC), ycf1, and ycf2. Maximum likelihood and Bayesian phylogenetic analyses resulted in nearly identical tree topologies and produced a well-resolved backbone comprising 15 consecutive dichotomies that subdivided Apocynaceae into 15 clades. The subfamily Periplocoideae were embedded in the Apocynoid grade and were sister to the Echiteae-Odontadenieae-Mesechiteae clade with high support values. Three tribes (Melodineae, Vinceae, and Willughbeieae), the subtribe Amphineuriinae, and four genera (Beaumontia, Ceropegia, Hoya, and Stephanotis) were not resolved as monophyletic. Our work sheds light on the backbone phylogenetic relationships in the family Apocynaceae and offers insights into the evolution of Apocynaceae plastomes using the most densely sampled plastome dataset to date.

A comparative study across the parasitic plants of Cuscuta subgenus Grammica (Convolvulaceae) reveals a possible loss of the plastid genome in its section Subulatae

Most species in Cuscuta subgenus Grammica retain many photosynthesis-related plastid genes, generally under purifying selection. A group of holoparasitic species in section Subulatae may have lost their plastid genomes entirely. The c. 153 species of plants belonging to Cuscuta subgenus Grammica are all obligate stem parasites. However, some have completely lost the ability to conduct photosynthesis while others retain photosynthetic machinery and genes. The plastid genome that primarily encodes key photosynthesis genes functions as a bellwether for how reliant plants are on primary production. This research assembles and analyses 17 plastomes across Cuscuta subgenus Grammica with the aim of characterizing the state of the plastome in each of its sections. By comparing the structure and content of plastid genomes across the subgenus, as well as by quantifying the selection acting upon each gene, we reconstructed the patterns of plastome change within the phylogenetic context for this group. We found that species in 13 of the 15 sections that comprise Grammica retain the bulk of plastid photosynthesis genes and are thus hemiparasitic. The complete loss of photosynthesis can be traced to two clades: the entire section Subulatae and a complex of three species within section Ceratophorae. We were unable to recover any significant plastome sequences from section Subulatae, suggesting that plastomes in these species are either drastically reduced or lost entirely.

Poaceae Chloroplast Genome Sequencing: Great Leap Forward in Recent Ten Years

The first complete chloroplast genome of rice (Oryza sativa) was published in 1989, ushering in a new era of studies of chloroplast genomics in Poaceae. Progresses in Next-Generation Sequencing (NGS) and Third-Generation Sequencing (TGS) technologiesand in the development of genome assembly software, have significantly advanced chloroplast genomics research. Poaceae is one of the most targeted families in chloroplast genome research because of its agricultural, ecological, and economic importance. Over the last 30 years, 2,050 complete chloroplast genome sequences from 40 tribes and 282 genera have been generated, most (97%) of them in the recent ten years. The wealth of data provides the groundwork for studies on species evolution, phylogeny, genetic transformation, and other aspects of Poaceae chloroplast genomes. As a result, we have gained a deeper understanding of the properties of Poaceae chloroplast genomes. Here, we summarize the achievements of the studies of the Poaceae chloroplast genomes and envision the challenges for moving the area ahead.

Anticodon table of the chloroplast genome and identification of putative quadruplet anticodons in chloroplast tRNAs

The chloroplast genome of 5959 species was analyzed to construct the anticodon table of the chloroplast genome. Analysis of the chloroplast transfer ribonucleic acid (tRNA) revealed the presence of a putative quadruplet anticodon containing tRNAs in the chloroplast genome. The tRNAs with putative quadruplet anticodons were UAUG, UGGG, AUAA, GCUA, and GUUA, where the GUUA anticodon putatively encoded tRNA^Asn. The study also revealed the complete absence of tRNA genes containing ACU, CUG, GCG, CUC, CCC, and CGG anticodons in the chloroplast genome from the species studied so far. The chloroplast genome was also found to encode tRNAs encoding N-formylmethionine (fMet), Ile2, selenocysteine, and pyrrolysine. The chloroplast genomes of mycoparasitic and heterotrophic plants have had heavy losses of tRNA genes. Furthermore, the chloroplast genome was also found to encode putative spacer tRNA, tRNA fragments (tRFs), tRNA-derived, stress-induced RNA (tiRNAs), and the group I introns. An evolutionary analysis revealed that chloroplast tRNAs had evolved via multiple common ancestors and the GC% had more influence toward encoding the tRNA number in the chloroplast genome than the genome size.

Decoding the Virtual 2D Map of the Chloroplast Proteomes

BACKGROUND

The chloroplast is a semi-autonomous organelle having its own genome and corresponding proteome. Although chloroplast genomes have been reported, no reports exist on their corresponding proteomes. Therefore, a proteome-wide analysis of the chloroplast proteomes of 2893 species was conducted, and a virtual 2D map was constructed.

RESULTS

The resulting virtual 2D map of the chloroplast proteome exhibited a bimodal distribution. The molecular mass of the chloroplast proteome ranged from 0.448 to 616.334 kDa, and the isoelectric point (pI) ranged from 2.854 to 12.954. Chloroplast proteomes were dominated by basic pI proteins with an average pI of 7.852. The molecular weight and isoelectric point of chloroplast proteome were found to show bimodal distribution. Leu was the most abundant and Cys the least abundant amino acid in the chloroplast proteome. Notably, Trp amino acid was absent in the chloroplast protein sequences of Pilostyles aethiopica. In addition, Selenocysteine (Sec) and Pyrrolysine (Pyl) amino acids were also found to be lacking in the chloroplast proteomes.

CONCLUSION

The virtual 2D map and amino acid composition of chloroplast proteome will enable the researchers to understand the biochemistry of chloroplast protein in detail. Further, the amino acid composition of the chloroplast proteome will also allow us to understand the codon usage bias. The codon usage bias and amino acid usage bias of chloroplast will be crucial to understanding their relationship.

Evidence from Phylogenomics and Morphology Provide Insights into the Phylogeny, Plastome Evolution, and Taxonomy of Kitagawia

Kitagawia Pimenov is one of the segregate genera of Peucedanum sensu lato within the Apiaceae. The phylogenetic position and morphological delimitation of Kitagawia have been controversial. In this study, we used plastid genome (plastome) and nuclear ribosomal DNA (nrDNA) sequences to reconstruct the phylogeny of Kitagawia, along with comparative plastome and morphological analyses between Kitagawia and related taxa. The phylogenetic results identified that all examined Kitagawia species were divided into Subclade I and Subclade II within the tribe Selineae, and they were all distant from the representative members of Peucedanum sensu stricto. The plastomes of Kitagawia and related taxa showed visible differences in the LSC/IRa junction (JLA) and several hypervariable regions, which separated Subclade I and Subclade II from other taxa. Fruit anatomical and micromorphological characteristics, as well as general morphological characteristics, distinguished the four Kitagawia species within Subclade I from Subclade II and other related genera. This study supported the separation of Kitagawia from Peucedanum sensu lato, confirmed that Kitagawia belongs to Selineae, and two species (K. praeruptora and K. formosana) within Subclade II should be placed in a new genus. We believe that the "core" Kitagawia should be limited to Subclade I, and this genus can be distinguished by the association of a series of morphological characteristics. Overall, our study provides new insights into the phylogeny, plastome evolution, and taxonomy of Kitagawia.

Evolutionary inference across eukaryotes identifies universal features shaping organelle gene retention

Mitochondria and plastids power complex life. Why some genes and not others are retained in their organelle DNA (oDNA) genomes remains a debated question. Here, we attempt to identify the properties of genes and associated underlying mechanisms that determine oDNA retention. We harness over 15k oDNA sequences and over 300 whole genome sequences across eukaryotes with tools from structural biology, bioinformatics, machine learning, and Bayesian model selection. Previously hypothesized features, including the hydrophobicity of a protein product, and less well-known features, including binding energy centrality within a protein complex, predict oDNA retention across eukaryotes, with additional influences of nucleic acid and amino acid biochemistry. Notably, the same features predict retention in both organelles, and retention models learned from one organelle type quantitatively predict retention in the other, supporting the universality of these features-which also distinguish gene profiles in more recent, independent endosymbiotic relationships. A record of this paper's transparent peer review process is included in the supplemental information.

Plastome characteristics and species identification of Chinese medicinal wintergreens (Gaultheria, Ericaceae)

Wintergreen oil is a folk medicine widely used in foods, pesticides, cosmetics and drugs. In China, nine out of 47 species within Gaultheria (Ericaceae) are traditionally used as Chinese medicinal wintergreens; however, phylogenetic approaches currently used to discriminating these species remain unsatisfactory. In this study, we sequenced and characterized plastomes from nine Chinese wintergreen species and identified candidate DNA barcoding regions for Gaultheria. Each Gaultheria plastome contained 110 unique genes (76 protein-coding, 30 tRNA, and four rRNA genes). Duplication of trnfM, rps14, and rpl23 genes were detected, while all plastomes lacked ycf1 and ycf2 genes. Gaultheria plastomes shared substantially contracted SSC regions that contained only the ndhF gene. Moreover, plastomes of Gaultheria leucocarpa var. yunnanensis contained an inversion in the LSC region and an IR expansion to cover the ndhF gene. Multiple rearrangement events apparently occurred between the Gaultheria plastomes and those from several previously reported families in Ericales. Our phylogenetic reconstruction using 42 plastomes revealed well-supported relationships within all nine Gaultheria species. Additionally, seven mutational hotspot regions were identified as potential DNA barcodes for Chinese medicinal wintergreens. Our study is the first to generate complete plastomes and describe the structural variations of the complicated genus Gaultheria. In addition, our findings provide important resources for identification of Chinese medicinal wintergreens.

Characterization of the plastome of Physaliscordata and comparative analysis of eight species of Physalis sensu stricto

In this study, we sequenced, assembled, and annotated the plastome of Physaliscordata Mill. and compared it with seven species of the genus Physalis sensu stricto. Sequencing, annotating, and comparing plastomes allow us to understand the evolutionary mechanisms associated with physiological functions, select possible molecular markers, and identify the types of selection that have acted in different regions of the genome. The plastome of P.cordata is 157,000 bp long and presents the typical quadripartite structure with a large single-copy (LSC) region of 87,267 bp and a small single-copy (SSC) region of 18,501 bp, which are separated by two inverted repeat (IRs) regions of 25,616 bp each. These values are similar to those found in the other species, except for P.angulata L. and P.pruinosa L., which presented an expansion of the LSC region and a contraction of the IR regions. The plastome in all Physalis species studied shows variation in the boundary of the regions with three distinct types, the percentage of the sequence identity between coding and non-coding regions, and the number of repetitive regions and microsatellites. Four genes and 10 intergenic regions show promise as molecular markers and eight genes were under positive selection. The maximum likelihood analysis showed that the plastome is a good source of information for phylogenetic inference in the genus, given the high support values and absence of polytomies. In the Physalis plastomes analyzed here, the differences found, the positive selection of genes, and the phylogenetic relationships do not show trends that correspond to the biological or ecological characteristics of the species studied.

Newly reported chloroplast genome of Sinosenecio albonervius Y. Liu & Q. E. Yang and comparative analyses with other Sinosenecio species

Background

Sinosenecio B. Nordenstam (Asteraceae) currently comprises 44 species. To investigate the interspecific relationship, several chloroplast markers, including ndhC-trnV, rpl32-trnL, matK, and rbcL, are used to analyze the phylogeny of Sinosenecio. However, the chloroplast genomes of this genus have not been thoroughly investigated. We sequenced and assembled the Sinosenecio albonervius chloroplast genome for the first time. A detailed comparative analysis was performed in this study using the previously reported chloroplast genomes of three Sinosenecio species.

Results

The results showed that the chloroplast genomes of four Sinosenecio species exhibit a typical quadripartite structure. There are equal numbers of total genes, protein-coding genes and RNA genes among the annotated genomes. Per genome, 49–56 simple sequence repeats and 99 repeat sequences were identified. Thirty codons were identified as RSCU values greater than 1 in the chloroplast genome of S. albonervius based on 54 protein-coding genes, indicating that they showed biased usage. Among 18 protein-coding genes, 46 potential RNA editing sites were discovered. By comparing these chloroplast genomes' structures, inverted repeat regions and coding regions were more conserved than single-copy and non-coding regions. The junctions among inverted repeat and single-copy regions showed slight difference. Several hot spots of genomic divergence were detected, which can be used as new DNA barcodes for species identification. Phylogenetic analysis of the whole chloroplast genome showed that the four Sinosenecio species have close interspecific relationships.

Conclusions

The complete chloroplast genome of Sinosenecio albonervius was revealed in this study, which included a comparison of Sinosenecio chloroplast genome structure, variation, and phylogenetic analysis for related species. These will help future research on Sinosenecio taxonomy, identification, origin, and evolution to some extent.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-022-08872-3.

Comparative Analysis of Bacillariophyceae Chloroplast Genomes Uncovers Extensive Genome Rearrangements Associated with Speciation

The Bacillariophyceae is a species-rich, ecologically significant class of Bacillariophyta. Despite their critical importance in marine ecosystems as primary producers and in the development of harmful algal blooms (HABs), taxonomic research on Bacillariophyceae species has been hindered because of their limited morphological features, plasticity of morphologies, and the low resolution of common molecular markers. Hence molecular markers with improved resolution are urgently needed. Organelle genomes, which can be constructed efficiently with the recent development of high throughput DNA sequencing technologies and the advancement of bioinformatics tools, have been proposed as super barcodes for their higher resolution for distinguishing different species and intra-species genomic variations. In this study, we tested the value of full-length chloroplast genomes (cpDNAs) as super barcodes for distinguishing diatom species, by constructing cpDNAs of 11 strains of the class Bacillariophyceae, including Nitzschia ovalis, Nitzschia traheaformis, Cylindrotheca spp., Psammodictyon constrictum, Bacillaria paxillifer, two strains of Haslea tsukamotoi, Haslea avium, Navicula arenaria, and Pleurosigma sp. Comparative analysis of cpDNAs revealed that cpDNAs were not only adequate for resolving different species, but also for enabling recognition of high levels of genome rearrangements between cpDNAs of different species, especially for species of the genera Nitzschia, Cylindrotheca, Navicula and Haslea. Additionally, comparative analysis suggested that the positioning of species in the genus Haslea should be transferred to the genus Navicula. Chloroplast genome-based evolutionary analysis suggested that the Bacillariophyceae species first appeared during the Cretaceous period and the diversity of species rose after the mass extinction about 65 Mya. This study highlighted the value of cpDNAs in research on the biodiversity and evolution of Bacillariophyceae species, and, with the construction of more cpDNAs representing additional genera, deeper insight into the biodiversity and evolutionary relationships of Bacillariophyceae species will be gained.

The complete chloroplast genome of Carex laevissima Nakai (Cyperaceae)

Carex laevissima Nakai 1914 (Cyperaceae) is vital for ecological conservation and land virescence, and has high ornamental value. Here the chloroplast genome of Carex laevissima was assembled and systematically analyzed for further genetic research of Carex plants. The chloroplast sequence of Carex laevissima was 188,029 bp in length, including two inverted repeat (IR) regions of 36,699 bp each, a large single-copy (LSC) region of 106,171 bp and a small single-copy region (SSC) of 8460 bp. The overall GC content is 34.0%. It contains 133 genes, including 89 protein-coding, 36 tRNA, and eight rRNA genes. Phylogenetic analysis showed that Carex laevissima is most closely related to Carex neurocarpa.

Strong plastid degradation is consistent within section Chondrophyllae, the most speciose lineage of Gentiana

Recovering phylogenetic relationships in lineages experiencing intense diversification has always been a persistent challenge in evolutionary studies, including in Gentiana section Chondrophyllae sensu lato (s.l.). Indeed, this subcosmopolitan taxon encompasses more than 180 mostly annual species distributed around the world. We sequenced and assembled 22 new plastomes representing 21 species in section Chondrophyllae s.l. In addition to previously released plastome data, our study includes all main lineages within the section. We reconstructed their phylogenetic relationships based on protein‐coding genes and recombinant DNA (rDNA) cistron sequences, and then investigated plastome structural evolution as well as divergence time. Despite an admittedly humble species cover overall, we recovered a well‐supported phylogenetic tree based on plastome data, and found significant discordance between phylogenetic relationships and taxonomic treatments. Our results show that G. capitata and G. leucomelaena diverged early within the section, which is then further divided into two clades. The divergence time estimation showed that section Chondrophyllae s.l. evolved in the second half of the Oligocene. We found that section Chondrophyllae s.l. had the smallest average plastome size (128 KB) in tribe Gentianeae (Gentianaceae), with frequent gene and sequence losses such as the ndh complex and its flanking regions. In addition, we detected both expansion and contraction of the inverted repeat (IR) regions. Our study suggests that plastome degradation parallels the diversification of this group, and illustrates the strong discordance between phylogenetic relationships and taxonomic treatments, which now need to be carefully revised.

Software Choice and Sequencing Coverage Can Impact Plastid Genome Assembly-A Case Study in the Narrow Endemic Calligonum bakuense

Most plastid genome sequences are assembled from short-read whole-genome sequencing data, yet the impact that sequencing coverage and the choice of assembly software can have on the accuracy of the resulting assemblies is poorly understood. In this study, we test the impact of both factors on plastid genome assembly in the threatened and rare endemic shrub Calligonum bakuense. We aim to characterize the differences across plastid genome assemblies generated by different assembly software tools and levels of sequencing coverage and to determine if these differences are large enough to affect the phylogenetic position inferred for C. bakuense compared to congeners. Four assembly software tools (FastPlast, GetOrganelle, IOGA, and NOVOPlasty) and seven levels of sequencing coverage across the plastid genome (original sequencing depth, 2,000x, 1,000x, 500x, 250x, 100x, and 50x) are compared in our analyses. The resulting assemblies are evaluated with regard to reproducibility, contig number, gene complement, inverted repeat length, and computation time; the impact of sequence differences on phylogenetic reconstruction is assessed. Our results show that software choice can have a considerable impact on the accuracy and reproducibility of plastid genome assembly and that GetOrganelle produces the most consistent assemblies for C. bakuense. Moreover, we demonstrate that a sequencing coverage between 500x and 100x can reduce both the sequence variability across assembly contigs and computation time. When comparing the most reliable plastid genome assemblies of C. bakuense, a sequence difference in only three nucleotide positions is detected, which is less than the difference potentially introduced through software choice.

Chloroplast Genome of Lithocarpus dealbatus (Hook.f. & Thomson ex Miq.) Rehder Establishes Monophyletic Origin of the Species and Reveals Mutational Hotspots with Taxon Delimitation Potential

There is phylogenetic ambiguity in the genus Lithocarpus and subfamily Quercoideae (Family: Fagaceae). Lithocarpus dealbatus, an ecologically important tree, is the dominant species among the Quercoideae in India. Although several studies have been conducted on the species' regeneration and ecological and economic significance, limited information is available on its phylo-genomics. To resolve the phylogeny in Quercoideae, we sequenced and assembled the 161,476 bp chloroplast genome of L. dealbatus, which has a large single-copy section of 90,732 bp and a small single-copy region of 18,987 bp, separated by a pair of inverted repeat regions of 25,879 bp. The chloroplast genome contained 133 genes, of which 86 were protein-coding genes, 39 were transfer RNAs, and eight were ribosomal RNAs. Analysis of repeat elements and RNA editing sites revealed interspecific similarities within the Lithocarpus genus. DNA diversity analysis identified five highly diverged coding and noncoding hotspot regions in the four genera, which can be used as polymorphic markers for species/taxon delimitation across the four genera of Quercoideae viz., Lithocarpus, Quercus, Castanea, and Castanopsis. The chloroplast-based phylogenetic analysis among the Quercoideae established a monophyletic origin of Lithocarpus, and a closer evolutionary lineage with a few Quercus species. Besides providing insights into the chloroplast genome architecture of L. dealbatus, the study identified five mutational hotspots having high taxon-delimitation potential across four genera of Quercoideae.

Complete Chloroplast Genome of Cnidium monnieri (Apiaceae) and Comparisons with Other Tribe Selineae Species

Cnidium monnieri is an economically important traditional Chinese medicinal plant. In this study, the complete chloroplast (cp) genome of C. monnieri was determined using the Illumina paired-end sequencing, the GetOrganelle de novo assembly strategy, as well as the GeSeq annotation method. Our results showed that the cp genome was 147,371 bp in length with 37.4% GC content and included a large single-copy region (94,361 bp) and a small single-copy region (17,552 bp) separated by a pair of inverted repeat regions (17,729 bp). A total of 129 genes were contained in the cp genome, including 85 protein-coding genes, 36 tRNA genes, and eight rRNA genes. We also investigated codon usage, RNA editing, repeat sequences, simple sequence repeats (SSRs), IR boundaries, and pairwise Ka/Ks ratios. Four hypervariable regions (trnD-trnY-trnE-trnT, ycf2, ndhF-rpl32-trnL, and ycf1) were identified as candidate molecular markers for species authentication. The phylogenetic analyses supported non-monophyly of Cnidium and C. monnieri located in tribe Selineae based on the cp genome sequences and internal transcribed spacer (ITS) sequences. The incongruence of the phylogenetic position of C. monnieri between ITS and cpDNA phylogenies suggested that C. monnieri might have experienced complex evolutions with hybrid and incomplete lineage sorting. All in all, the results presented herein will provide plentiful chloroplast genomic resources for studies of the taxonomy, phylogeny, and species authentication of C. monnieri. Our study is also conducive to elucidating the phylogenetic relationships and taxonomic position of Cnidium.

Plastid genomes reveal evolutionary shifts in elevational range and flowering time of Osmanthus (Oleaceae)

Species of Osmanthus are economically important ornamental trees, yet information regarding their plastid genomes (plastomes) have rarely been reported, thus hindering taxonomic and evolutionary studies of this small but enigmatic genus. Here, we performed comparative genomics and evolutionary analyses on plastomes of 16 of the 28 currently accepted species, with 11 plastomes newly sequenced. Phylogenetic studies identified four main lineages within the genus that are here designated the: “Caucasian Osmanthus” (corresponding to O. decorus), “Siphosmanthus” (corresponding to O. sect. Siphosmanthus), “O. serrulatus + O. yunnanensis,” and “Core Osmanthus: (corresponding to O. sect. Osmanthus + O. sect. Linocieroides). Molecular clock analysis suggested that Osmanthus split from its sister clade c. 15.83 Ma. The estimated crown ages of the lineages were the following: genus Osmanthus at 12.66 Ma; “Siphosmanthus” clade at 5.85 Ma; “O. serrulatus + O. yunnanensis” at 4.89 Ma; and “Core Osmanthus: clade at 6.2 Ma. Ancestral state reconstructions and trait mapping showed that ancestors of Osmanthus were spring flowering and originated at lower elevations. Phylogenetic principal component analysis clearly distinguished spring‐flowering species from autumn‐flowering species, suggesting that flowering time differentiation is related to the difference in ecological niches. Nucleotide substitution rates of 80 common genes showed slow evolutionary pace and low nucleotide variations, all genes being subjected to purifying selection.

Characterizing the complete chloroplast genome of the Impatiens davidii (Balsaminaceae)

Impatiens davidii Franch, 1886 is a rare ornamental flower used in gardens and has high economic value. In this study, we characterized the chloroplast genome of I. davidii and analyzed its phylogenetic relationship with other Impatiens species. The length of the complete chloroplast genome sequence of I. davidii is 152,214 bp, with a GC content of 36.9%. The chloroplast genome shows a typical quadripartite structure with a pair of inverted repeats (IRs) of 25,634 bp, separated by one large single copy (LSC) region of 83,128 bp and one small single copy (SSC) region of 17,818 bp. We annotated 125 genes, of which there were 85 protein-coding genes, 32 tRNA genes, and 8 rRNA genes. The Bayesian phylogenetic tree strongly supports that I. davidii has a close phylogenetic relationship with a group including I. piufanensis and I. alpicola.

Characterization of the complete chloroplast genome of Zephyranthes phycelloides (Amaryllidaceae, tribe Hippeastreae) from Atacama region of Chile

Sporadic rains in the Atacama Desert reveal a high biodiversity of plant species that only occur there. One of these rare species is the “Red añañuca” (Zephyranthes phycelloides), formerly known as Rhodophiala phycelloides. Many species of Zephyranthes in the Atacama Desert are dangerously threatened, due to massive extraction of bulbs and cutting of flowers. Therefore, studies of the biodiversity of these endemic species, which are essential for their conservation, should be conducted sooner rather than later. There are some chloroplast genomes available for Amaryllidaceae species, however there is no complete chloroplast genome available for any of the species of Zephyranthes subgenus Myostemma. The aim of the present work was to characterize and analyze the chloroplast of Z. phycelloides by NGS sequencing. The chloroplast genome of the Z. phycelloides consists of 158,107 bp, with typical quadripartite structures: a large single copy (LSC, 86,129 bp), a small single copy (SSC, 18,352 bp), and two inverted repeats (IR, 26,813 bp). One hundred thirty-seven genes were identified: 87 coding genes, 8 rRNA, 38 tRNA and 4 pseudogenes. The number of SSRs was 64 in Z. phycelloides and a total of 43 repeats were detected. The phylogenetic analysis of Z. phycelloides shows a distinct subclade with respect to Z. mesochloa. The average nucleotide variability (Pi) between Z. phycelloides and Z. mesochloa was of 0.02000, and seven loci with high variability were identified: psbA, trnS^GCU-trnG^UCC, trnD^GUC-trnY^GUA, trnL^UAA-trnF^GAA, rbcL, psbE-petL and ndhG-ndhI. The differences between the species are furthermore confirmed by the high amount of SNPs between these two species. Here, we report for the first time the complete cp genome of one species of the Zephyranthes subgenus Myostemma, which can be used for phylogenetic and population genomic studies.

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

The plastid organelle is essential for many vital cellular processes and the growth and development of plants. The availability of a large number of complete plastid genomes could be effectively utilized to understand the evolution of the plastid genomes and phylogenetic relationships among plants. We comprehensively analyzed the plastid genomes of Viridiplantae comprising 3,654 taxa from 298 families and 111 orders and compared the genomic organizations in their plastid genomic DNA among major clades, which include gene gain/loss, gene copy number, GC content, and gene blocks. We discovered that some important genes that exhibit similar functions likely formed gene blocks, such as the psb family presumably showing co-occurrence and forming gene blocks in Viridiplantae. The inverted repeats (IRs) in plastid genomes have doubled in size across land plants, and their GC content is substantially higher than non-IR genes. By employing three different data sets [all nucleotide positions (nt123), only the first and second codon positions (nt12), and amino acids (AA)], our phylogenomic analyses revealed Chlorokybales + Mesostigmatales as the earliest-branching lineage of streptophytes. Hornworts, mosses, and liverworts forming a monophylum were identified as the sister lineage of tracheophytes. Based on nt12 and AA data sets, monocots, Chloranthales and magnoliids are successive sister lineages to the eudicots + Ceratophyllales clade. The comprehensive taxon sampling and analysis of different data sets from plastid genomes recovered well-supported relationships of green plants, thereby contributing to resolving some long-standing uncertainties in the plant phylogeny.

Extensive reorganization of the chloroplast genome of Corydalis platycarpa: A comparative analysis of their organization and evolution with other Corydalis plastomes

Introduction

The chloroplast (cp) is an autonomous plant organelle with an individual genome that encodes essential cellular functions. The genome architecture and gene content of the cp is highly conserved in angiosperms. The plastome of Corydalis belongs to the Papaveraceae family, and the genome is comprised of unusual rearrangements and gene content. Thus far, no extensive comparative studies have been carried out to understand the evolution of Corydalis chloroplast genomes.

Methods

Therefore, the Corydalis platycarpa cp genome was sequenced, and wide-scale comparative studies were conducted using publicly available twenty Corydalis plastomes.

Results

Comparative analyses showed that an extensive genome rearrangement and IR expansion occurred, and these events evolved independently in the Corydalis species. By contrast, the plastomes of its closely related subfamily Papaveroideae and other Ranunculales taxa are highly conserved. On the other hand, the synapomorphy characteristics of both accD and the ndh gene loss events happened in the common ancestor of the Corydalis and sub-clade of the Corydalis lineage, respectively. The Corydalis-sub clade species (ndh lost) are distributed predominantly in the Qinghai-Tibetan plateau (QTP) region. The phylogenetic analysis and divergence time estimation were also employed for the Corydalis species.

Discussion

The divergence time of the ndh gene in the Corydalis sub-clade species (44.31 - 15.71 mya) coincides very well with the uplift of the Qinghai-Tibet Plateau in Oligocene and Miocene periods, and maybe during this period, it has probably triggered the radiation of the Corydalis species.

Conclusion

To the best of the authors' knowledge, this is the first large-scale comparative study of Corydalis plastomes and their evolution. The present study may provide insights into the plastome architecture and the molecular evolution of Corydalis species.

Complete chloroplast genome of a cultivated oil camellia species, Camellia gigantocarpa

Camellia gigantocarpa Hu et T. C. Huang, belonging to the Theaceae family, is an excellent landscape tree species with high ornamental value. It is also an important woody oil-bearing plant with high economic value. This study reports the first complete chloroplast genome sequence of C. gigantocarpa (GenBank accession number: MZ054232). Its whole chloroplast genome is 156,953 bp long with an overall GC content of 37.31%, which is composed of a large single copy region (86,631 bp), a small single copy region (18,402 bp), and a pair of inverted repeat regions (25,960 bp each). A total of 135 genes were predicted in this genome, namely eight ribosomal RNA genes, 37 transfer RNA genes, and 90 protein-coding genes. Based on maximum likelihood analysis results, we found that the Camellia species are clustered into a distinct branch, and the phylogenetic relationships among C. gigantocarpa, C. crapnelliana, and C. kissii were the closest.

Complete chloroplast genome sequences of five Bruguiera species (Rhizophoraceae): comparative analysis and phylogenetic relationships

Bruguiera is a genus of true mangroves that are mostly distributed in the Indo-West Pacific region. However, the number of published whole chloroplast genome sequences of Bruguiera species are limited. Here, the complete chloroplast sequences of five Bruguiera species were sequenced and assembled using Illumina data. The chloroplast genomes of B. gymnorhiza, B. hainesii, B. cylindrica, B. parviflora and B. sexangula were assembled into 161,195, 164,295, 164,297, 163,228 and 164,170 bp, respectively. All chloroplast genomes contain 37 tRNA and eight rRNA genes, with either 84 or 85 protein-coding genes. A comparative analysis of these genomes revealed high similarity in gene structure, gene order and boundary position of the LSC, SSC and two IR regions. Interestingly, B. gymnorhiza lost a rpl32 gene in the SSC region. In addition, a ndhF gene in B. parviflora straddles both the SSC and IRB boundary regions. These genes reveal differences in chloroplast evolution among Bruguiera species. Repeats and SSRs in the chloroplast genome sequences were found to be highly conserved between B. cylindrica and B. hainesii as well as B. gymnorhiza and B. sexangula indicating close genetic relationships based on maternal inheritance. Notably, B. hainesii, which is considered a hybrid between B. gymnorhiza and B. cylindrica, appears to have inherited the chloroplast from B. cylindrica. Investigating the effects of selection events on shared protein-coding genes showed a positive selection in rps7 and rpl36 genes in all species compared to land-plant species. A phylogenetic analysis, based on 59 conserved chloroplast protein-coding genes, showed strong support that all Bruguiera species are in the clade Rhizophoraceae. This study provides valuable genetic information for the study of evolutionary relationships and population genetics in Bruguiera and other mangrove species.

Lineage-specific plastid degradation in subtribe Gentianinae (Gentianaceae)

The structure and sequence of plastid genomes is highly conserved across most land plants, except for a minority of lineages that show gene loss and genome degradation. Understanding the early stages of plastome degradation may provide crucial insights into the repeatability and predictability of genomic evolutionary trends. We investigated these trends in subtribe Gentianinae of the Gentianaceae, which encompasses ca. 450 species distributed around the world, particularly in alpine and subalpine environments. We sequenced, assembled, and annotated the plastomes of 41 species, representing all six genera in subtribe Gentianinae and all main sections of the species-rich genus Gentiana L. We reconstructed the phylogeny, estimated divergence times, investigated the phylogenetic distribution of putative gene losses, and related these to substitution rate shifts and species' habitats. We obtained a strongly supported topology consistent with earlier studies, with all six genera in Gentianinae recovered as monophyletic and all main sections of Gentiana having full support. While closely related species have very similar plastomes in terms of size and structure, independent gene losses, particularly of the ndh complex, have occurred in multiple clades across the phylogeny. Gene loss was usually associated with a shift in the boundaries of the small single-copy and inverted repeat regions. Substitution rates were variable between clades, with evidence for both elevated and decelerated rate shifts. Independent lineage-specific loss of ndh genes occurred at a wide range of times, from Eocene to Pliocene. Our study illustrates that diverse degradation patterns shape the evolution of the plastid in this species-rich plant group.

The Genomic Impact of Mycoheterotrophy in Orchids

Mycoheterotrophic plants have lost the ability to photosynthesize and obtain essential mineral and organic nutrients from associated soil fungi. Despite involving radical changes in life history traits and ecological requirements, the transition from autotrophy to mycoheterotrophy has occurred independently in many major lineages of land plants, most frequently in Orchidaceae. Yet the molecular mechanisms underlying this shift are still poorly understood. A comparison of the transcriptomes of Epipogium aphyllum and Neottia nidus-avis, two completely mycoheterotrophic orchids, to other autotrophic and mycoheterotrophic orchids showed the unexpected retention of several genes associated with photosynthetic activities. In addition to these selected retentions, the analysis of their expression profiles showed that many orthologs had inverted underground/aboveground expression ratios compared to autotrophic species. Fatty acid and amino acid biosynthesis as well as primary cell wall metabolism were among the pathways most impacted by this expression reprogramming. Our study suggests that the shift in nutritional mode from autotrophy to mycoheterotrophy remodeled the architecture of the plant metabolism but was associated primarily with function losses rather than metabolic innovations.