Full-length transcriptome analysis of Ophioglossum vulgatum: effects of experimentally identified chloroplast gene clusters on expression and evolutionary patterns.

Genes with similar or related functions in chloroplasts are often arranged in close proximity, forming clusters on chromosomes. These clusters are transcribed coordinated to facilitate the expression of genes with specific function. Our previous study revealed a significant negative correlation between the chloroplast gene expression level of the rare medicinal fern Ophioglossum vulgatum and its evolutionary rates as well as selection pressure. Therefore, in this study, we employed a combination of SMRT and Illumina sequencing technology to analyze the full-length transcriptome sequencing of O. vulgatum for the first time. In particular, we experimentally identified gene clusters based on transcriptome data and investigated the effects of chloroplast gene clustering on expression and evolutionary patterns. The results revealed that the total sequenced data volume of the full-length transcriptome of O. vulgatum amounted to 71,950,652,163 bp, and 110 chloroplast genes received transcript coverage. Nine different types of gene clusters were experimentally identified in their transcripts. The chloroplast cluster genes may cause a decrease in non-synonymous substitution rate and selection pressure, as well as a reduction in transversion rate, transition rate, and their ratio. While expression levels of chloroplast cluster genes in leaf, sporangium, and stem would be relatively elevated. The Mann-Whitney U test indicated statistically significant in the selection pressure, sporangia and leaves groups (P < 0.05). We have contributed novel full-length transcriptome data resources for ferns, presenting new evidence on the effects of chloroplast gene clustering on expression land evolutionary patterns, and offering new theoretical support for transgenic research through gene clustering.

Chloroplast gene control: unlocking RNA thermometer mechanisms in photosynthetic systems.

RNA thermometers offer straightforward, protein-independent methods to regulate gene expression at the post-transcriptional level. In this context, Chung and colleagues have discovered a revolutionary RNA thermometer in the chloroplast genome of Chlamydomonas reinhardtii. This will facilitate temperature-driven control of inducible transgene expression for biotechnology applications in plant and algal systems.

Context-Dependent Mutation Dynamics, Not Selection, Explains the Codon Usage Bias of Most Angiosperm Chloroplast Genes.

Two competing proposals about the degree to which selection affects codon usage of angiosperm chloroplast genes are examined. The first, based on observations that codon usage does not match expectations under the naïve assumption that base composition will be identical at all neutral sites, is that selection plays a significant role. The second is that codon usage is determined almost solely by mutation bias and drift, with selection influencing only one or two highly expressed genes, in particular psbA. First it is shown that, as a result of an influence of neighboring base composition on mutation dynamics, compositional biases are expected to be widely divergent at different sites in the absence of selection. The observed mutation properties are then used to predict expected neutral codon usage biases and to show that observed deviations from the naïve expectations are in fact expected given the context-dependent mutational dynamics. It is also shown that there is a match between the observed and expected codon usage when context effects are taken into consideration, with psbA being a notable exception. Overall, the data support the model that selection is not a widespread factor affecting the codon usage of angiosperm chloroplast genes and highlight the need to have an accurate model of mutational dynamics.

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.

KEY MESSAGE: 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 trnPGGG. 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.

On the Edge of Dispensability, the Chloroplast ndh Genes.

The polypeptides encoded by the chloroplast ndh genes and some nuclear genes form the thylakoid NADH dehydrogenase (Ndh) complex, homologous to the mitochondrial complex I. Except for Charophyceae (algae related to higher plants) and a few Prasinophyceae, all eukaryotic algae lack ndh genes. Among vascular plants, the ndh genes are absent in epiphytic and in some species scattered among different genera, families, and orders. The recent identification of many plants lacking plastid ndh genes allows comparison on phylogenetic trees and functional investigations of the ndh genes. The ndh genes protect Angiosperms under various terrestrial stresses, maintaining efficient photosynthesis. On the edge of dispensability, ndh genes provide a test for the natural selection of photosynthesis-related genes in evolution. Variable evolutionary environments place Angiosperms without ndh genes at risk of extinction and, probably, most extant ones may have lost ndh genes recently. Therefore, they are evolutionary endpoints in phylogenetic trees. The low number of sequenced plastid DNA and the long lifespan of some Gymnosperms lacking ndh genes challenge models about the role of ndh genes protecting against stress and promoting leaf senescence. Additional DNA sequencing in Gymnosperms and investigations into the molecular mechanisms of their response to stress will provide a unified model of the evolutionary and functional consequences of the lack of ndh genes.

A PPR Protein ACM1 Is Involved in Chloroplast Gene Expression and Early Plastid Development in Arabidopsis.

Chloroplasts cannot develop normally without the coordinated action of various proteins and signaling connections between the nucleus and the chloroplast genome. Many questions regarding these processes remain unanswered. Here, we report a novel P-type pentatricopeptide repeat (PPR) factor, named Albino Cotyledon Mutant1 (ACM1), which is encoded by a nuclear gene and involved in chloroplast development. Knock-down of ACM1 transgenic plants displayed albino cotyledons but normal true leaves, while knock-out of the ACM1 gene in seedlings was lethal. Fluorescent protein analysis showed that ACM1 was specifically localized within chloroplasts. PEP-dependent plastid transcript levels and splicing efficiency of several group II introns were seriously affected in cotyledons in the RNAi line. Furthermore, denaturing gel electrophoresis and Western blot experiments showed that the accumulation of chloroplast ribosomes was probably damaged. Collectively, our results indicate ACM1 is indispensable in early chloroplast development in Arabidopsis cotyledons.

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

KEY MESSAGE: 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.

Analysis of codon usage bias of chloroplast genes in Oryza species : Codon usage of chloroplast genes in Oryza species.

MAIN CONCLUSION: The codon usage bias in chloroplast genes of Oryza species was low and AT rich. The pattern of codon usage was different among Oryza species and mainly influenced by mutation pressure and natural selection. Codon usage bias (CUB) is the unequal usage of synonymous codons in which some codons are more preferred to others in the coding sequences of genes. It shows a species-specific property. We studied the patterns of codon usage and the factors that influenced the CUB of protein-coding chloroplast (cp) genes in 18 Oryza species as no work was yet reported. The nucleotide composition analysis revealed that the overall GC content of cp genes in different species of Oryza was lower than 50%, i.e., Oryza cp genes were AT rich. Synonymous codon usage order (SCUO) suggested that CUB was weak in the cp genes of different Oryza species. A highly significant correlation was observed between overall nucleotides and its constituents at the third codon position suggesting that both, mutation pressure and natural selection, might influence the CUB. Correspondence analysis (COA) revealed that codon usage pattern differed across Oryza species. In the neutrality plot, a narrow range of GC3 distribution was recorded and some points were diagonally distributed in all the plots, suggesting that natural selection and mutation pressure might have influenced the CUB. The slope of the regression line was < 0.5, augmenting our inference that natural selection might have played a major role, while mutation pressure had a minor role in shaping the CUB of cp genes. The magnitudes of mutation pressure and natural selection on cp genes varied across Oryza species.

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.

EMB-7L is required for embryogenesis and plant development in maize involved in RNA splicing of multiple chloroplast genes.

Embryo and endosperm originate from the double fertilization, but they have different developmental fates and biological functions. We identified a previously undescribed maize seed mutant, wherein the embryo appears to be more severely affected than the endosperm (embryo-specific, emb). In the W22 background, the emb embryo arrests at the transition stage whereas its endosperm appears nearly normal in size. At maturity, the embryo in W22-emb is apparently small or even invisible. In contrast, the emb endosperm develops into a relative normal size. We cloned the mutant gene on the Chromosome 7L and designated it emb-7L. This gene is generally expressed, but it has a relatively higher expression level in leaves. Emb-7L encodes a chloroplast-localized P-type pentatricopeptide repeat (PPR) protein, consistent with the severe chloroplast deficiency in emb-7L albino seedling leaves. Full transcriptome analysis of the leaves of WT and emb-7L seedlings reveals that transcription of chloroplast protein-encoding genes are dramatically variable with pre-mRNA intron splicing apparently affected in a tissue-dependent pattern and the chloroplast structure and activity were dramatically affected including chloroplast membrane and photosynthesis machinery component and synthesis of metabolic products (e.g., fatty acids, amino acids, starch).

Mutation pressure and natural selection on codon usage in chloroplast genes of two species in Pisum L. (Fabaceae: Faboideae).

This study was attempted to focus on the pattern of codon usage bias (CUB) of chloroplast genes in two species of Pisum viz. P. fulvum and P. sativum and to identify the factors which influence CUB. Bioinformatic tools were used to understand codon usage pattern in the protein-coding sequences of Pisum chloroplast genomes. It was found that GC content was lower than AT content in the genes. Low synonymous codon usage order (SCUO) values of genes indicated low CUB in chloroplast genes of Pisum species. Heatmaps showed positive correlations of GC3 with all the GC and AT ending codons. Neutrality plot analysis revealed that natural selection might have played a prominent role over mutation pressure in sculpturing the CUB of chloroplast genes in these two taxa. Positive correlation between SCUO and mRNA free energy (mFE) suggested that higher energy release by entire mRNA was related to high degree of CUB. Further, highly significant (p < .01) negative correlation was found between parameters in pair i.e. mFE-GC, mFE-GC1, mFE-GC2 and mFE for entire mRNA-GC3. This pointed out that higher GC content might have influenced lesser energy release by mRNA molecules of chloroplast genes.

Evidence of mitochondrial DNA in the chloroplast genome of Convallaria keiskei and its subsequent evolution in the Asparagales.

DNA transfer between internal organelles such as the nucleus, mitochondrion, and plastid is a well-known phenomenon in plant evolution, and DNA transfer from the plastid and mitochondrion to the nucleus, from the plastid to the mitochondrion, and from the nucleus to the mitochondrion has been well-documented in angiosperms. However, evidence of the transfer of mitochondrial DNA (mtDNA) to the plastid has only been found in three dicotyledons and one monocotyledon. In the present study, we characterised and analysed two chloroplast (cp) genome sequences of Convallaria keiskei and Liriope spicata, and found that C. keiskei has the largest cp genome (162,109 bp) in the Asparagaceae. Interestingly, C. keiskei had a ~3.3-kb segment of mtDNA in its cp genome and showed similarity with the mt gene rpl10 as a pseudogene. Further analyses revealed that mtDNA transfer only occurred in C. keiskei in the Nolinoideae, which diverged very recently (7.68 million years ago (mya); 95% highest posterior density (HPD): 14.55-2.97 mya). These findings indicate that the C. keiskei cp genome is unique amongst monocotyledon land plants, but further work is necessary to understand the direction and mechanism involved in the uptake of mtDNA by the plastid genome of C. keiskei.

A nuclear-encoded protein, mTERF6, mediates transcription termination of rpoA polycistron for plastid-encoded RNA polymerase-dependent chloroplast gene expression and chloroplast development.

The expression of plastid genes is regulated by two types of DNA-dependent RNA polymerases, plastid-encoded RNA polymerase (PEP) and nuclear-encoded RNA polymerase (NEP). The plastid rpoA polycistron encodes a series of essential chloroplast ribosome subunits and a core subunit of PEP. Despite the functional importance, little is known about the regulation of rpoA polycistron. In this work, we show that mTERF6 directly associates with a 3'-end sequence of rpoA polycistron in vitro and in vivo, and that absence of mTERF6 promotes read-through transcription at this site, indicating that mTERF6 acts as a factor required for termination of plastid genes' transcription in vivo. In addition, the transcriptions of some essential ribosome subunits encoded by rpoA polycistron and PEP-dependent plastid genes are reduced in the mterf6 knockout mutant. RpoA, a PEP core subunit, accumulates to about 50% that of the wild type in the mutant, where early chloroplast development is impaired. Overall, our functional analyses of mTERF6 provide evidence that it is more likely a factor required for transcription termination of rpoA polycistron, which is essential for chloroplast gene expression and chloroplast development.

Differential impact of heat stress on the expression of chloroplast-encoded genes.

Heat shock is one of the major abiotic factors that causes severe retardation in plant growth and development. To dissect the principal effects of hyperthermia on chloroplast gene expression, we studied the temporal dynamics of transcript accumulation for chloroplast-encoded genes in Arabidopsis thaliana and genes for the chloroplast transcription machinery against a background of changes in physiological parameters. A marked reduction in the transcript amounts of the majority of the genes at the early phases of heat shock (HS) was followed by a return to the baseline levels of rbcL and the housekeeping genes clpP, accD, rps14 and rrn16. The decline in the mRNA levels of trnE (for tRNAglu) and the PSI genes psaA and psaB was opposed by the transient increase in the transcript accumulation of ndhF and the PSII genes psbA, psbD, and psbN and their subsequent reduction with the development of stress. However, the up-regulation of PSII genes in response to elevated temperature was absent in the heat stress-sensitive mutants abi1 and abi2 with the impaired degradation of D2 protein. The expression of rpoA and rpoB, which encode subunits of PEP, was strongly down-regulated throughout the duration of the heat treatment. In addition, heat stress-induced PEP deficiency caused the compensatory up-regulation of the genes for the nuclear-encoded RNA polymerases RPOTp and RPOTmp, the PEP-associated proteins PAP6 and PAP8, the Ser/Thr protein kinase cPCK2, and the stress-inducible sigma factor gene SIG5. Thus, heat stress differentially modulates the transcript accumulation of plastid-encoded genes in A. thaliana at least in part via the expression of HS-responsive nuclear genes for the plastid transcription machinery.

Genome and evolution of the shade-requiring medicinal herb Panax ginseng.

Panax ginseng C. A. Meyer, reputed as the king of medicinal herbs, has slow growth, long generation time, low seed production and complicated genome structure that hamper its study. Here, we unveil the genomic architecture of tetraploid P. ginseng by de novo genome assembly, representing 2.98 Gbp with 59 352 annotated genes. Resequencing data indicated that diploid Panax species diverged in association with global warming in Southern Asia, and two North American species evolved via two intercontinental migrations. Two whole genome duplications (WGD) occurred in the family Araliaceae (including Panax) after divergence with the Apiaceae, the more recent one contributing to the ability of P. ginseng to overwinter, enabling it to spread broadly through the Northern Hemisphere. Functional and evolutionary analyses suggest that production of pharmacologically important dammarane-type ginsenosides originated in Panax and are produced largely in shoot tissues and transported to roots; that newly evolved P. ginseng fatty acid desaturases increase freezing tolerance; and that unprecedented retention of chlorophyll a/b binding protein genes enables efficient photosynthesis under low light. A genome-scale metabolic network provides a holistic view of Panax ginsenoside biosynthesis. This study provides valuable resources for improving medicinal values of ginseng either through genomics-assisted breeding or metabolic engineering.

A novel chloroplast gene reported for flagellate plants.

PREMISE OF THE STUDY: Gene space in plant plastid genomes is well characterized and annotated, yet we discovered an unrecognized open reading frame (ORF) in the fern lineage that is conserved across flagellate plants. METHODS: We initially detected a putative uncharacterized ORF by the existence of a highly conserved region between rps16 and matK in a series of matK alignments of leptosporangiate ferns. We mined available plastid genomes for this ORF, which we now refer to as ycf94, to infer evolutionary selection pressures and assist in functional prediction. To further examine the transcription of ycf94, we assembled the plastid genome and sequenced the transcriptome of the leptosporangiate fern Adiantum shastense Huiet & A.R. Sm. KEY RESULTS: The ycf94 predicted protein has a distinct transmembrane domain but with no sequence homology to other proteins with known function. The nonsynonymous/synonymous substitution rate ratio of ycf94 is on par with other fern plastid protein-encoding genes, and additional homologs can be found in a few lycophyte, moss, hornwort, and liverwort plastid genomes. Homologs of ycf94 were not found in seed plants. In addition, we report a high level of RNA editing for ycf94 transcripts-a hallmark of protein-coding genes in fern plastomes. CONCLUSIONS: The degree of sequence conservation, together with the presence of a distinct transmembrane domain and RNA-editing sites, suggests that ycf94 is a protein-coding gene of functional significance in ferns and, potentially, bryophytes and lycophytes. However, the origin and exact function of this gene require further investigation.

Comparative chloroplast genomes of Paris Sect. Marmorata: insights into repeat regions and evolutionary implications.

BACKGROUND: Species of Paris Sect. Marmorata are valuable medicinal plants to synthesize steroidal saponins with effective pharmacological therapy. However, the wild resources of the species are threatened by plundering exploitation before the molecular genetics studies uncover the genomes and evolutionary significance. Thus, the availability of complete chloroplast genome sequences of Sect. Marmorata is necessary and crucial to the understanding the plastome evolution of this section and facilitating future population genetics studies. Here, we determined chloroplast genomes of Sect. Marmorata, and conducted the whole chloroplast genome comparison. RESULTS: This study presented detailed sequences and structural variations of chloroplast genomes of Sect. Marmorata. Over 40 large repeats and approximately 130 simple sequence repeats as well as a group of genomic hotspots were detected. Inverted repeat contraction of this section was inferred via comparing the chloroplast genomes with the one of P. verticillata. Additionally, almost all the plastid protein coding genes were found to prefer ending with A/U. Mutation bias and selection pressure predominately shaped the codon bias of most genes. And most of the genes underwent purifying selection, whereas photosynthetic genes experienced a relatively relaxed purifying selection. CONCLUSIONS: Repeat sequences and hotspot regions can be scanned to detect the intraspecific and interspecific variability, and selected to infer the phylogenetic relationships of Sect. Marmorata and other species in subgenus Daiswa. Mutation and natural selection were the main forces to drive the codon bias pattern of most plastid protein coding genes. Therefore, this study enhances the understanding about evolution of Sect. Marmorata from the chloroplast genome, and provide genomic insights into genetic analyses of Sect. Marmorata.

Comparative assessment of chloroplast transcriptional responses highlights conserved and unique patterns across Triticeae members under salt stress.

Chloroplast functional genomics, in particular understanding the chloroplast transcriptional response is of immense importance mainly due to its role in oxygenic photosynthesis. As a photosynthetic unit, its efficiency and transcriptional activity is directly regulated by reactive oxygen species during abiotic and biotic stress and subsequently affects carbon assimilation, and plant biomass. In crops, understanding photosynthesis is crucial for crop domestication by identifying the traits that could be exploited for crop improvement. Transcriptionally and translationally active chloroplast plays a key role by regulating the PSI and PSII photo-reaction centres, which ubiquitously affects the light harvesting. Using a comparative transcriptomics mapping approach, we identified differential regulation of key chloroplast genes during salt stress across Triticeae members with potential genes involved in photosynthesis and electron transport system such as CytB6f. Apart from differentially regulated genes involved in PSI and PSII, we found widespread evidence of intron splicing events, specifically uniquely spliced petB and petD in Triticum aestivum and high proportion of RNA editing in ndh genes across the Triticeae members during salt stress. We also highlight the role and differential regulation of ATP synthase as member of CF0CF1 and also revealed the effect of salt stress on the water-splitting complex under salt stress. It is worthwhile to mention that the observed conserved down-regulation of psbJ across the Triticeae is limiting the assembly of water-splitting complexes and thus making the BEP clade Triticeae members more vulnerable to high light during the salt stress. Comparative understanding of the chloroplast transcriptional dynamics and photosynthetic regulation will improve the approaches for improved crop domestication.

The Linum usitatissimum L. plastome reveals atypical structural evolution, new editing sites, and the phylogenetic position of Linaceae within Malpighiales.

KEY MESSAGE: The plastome of Linum usitatissimum was completely sequenced allowing analyses of evolution of genome structure, RNA editing sites, molecular markers, and indicating the position of Linaceae within Malpighiales. Flax (Linum usitatissimum L.) is an economically important crop used as food, feed, and industrial feedstock. It belongs to the Linaceae family, which is noted by high morphological and ecological diversity. Here, we reported the complete sequence of flax plastome, the first species within Linaceae family to have the plastome sequenced, assembled and characterized in detail. The plastome of flax is a circular DNA molecule of 156,721 bp with a typical quadripartite structure including two IRs of 31,990 bp separating the LSC of 81,767 bp and the SSC of 10,974 bp. It shows two expansion events from IRB to LSC and from IRB to SSC, and a contraction event in the IRA-LSC junction, which changed significantly the size and the gene content of LSC, SSC and IRs. We identified 109 unique genes and 2 pseudogenes (rpl23 and ndhF). The plastome lost the conserved introns of clpP gene and the complete sequence of rps16 gene. The clpP, ycf1, and ycf2 genes show high nucleotide and aminoacid divergence, but they still possibly retain the functionality. Moreover, we also identified 176 SSRs, 20 tandem repeats, and 39 dispersed repeats. We predicted in 18 genes a total of 53 RNA editing sites of which 32 were not found before in other species. The phylogenetic inference based on 63 plastid protein-coding genes of 38 taxa supports three major clades within Malpighiales order. One of these clades has flax (Linaceae) sister to Chrysobalanaceae family, differing from earlier studies that included Linaceae into the euphorbioid clade.

Complete chloroplast genome sequences of Lilium: insights into evolutionary dynamics and phylogenetic analyses.

Lilium is a large genus that includes approximately 110 species distributed throughout cold and temperate regions of the Northern Hemisphere. The species-level phylogeny of Lilium remains unclear; previous studies have found universal markers but insufficient phylogenetic signals. In this study, we present the use of complete chloroplast genomes to explore the phylogeny of this genus. We sequenced nine Lilium chloroplast genomes and retrieved seven published chloroplast genomes for comparative and phylogenetic analyses. The genomes ranged from 151,655 bp to 153,235 bp in length and had a typical quadripartite structure with a conserved genome arrangement and moderate divergence. A comparison of sixteen Lilium chloroplast genomes revealed ten mutation hotspots. Single nucleotide polymorphisms (SNPs) for any two Lilium chloroplast genomes ranged from 8 to 1,178 and provided robust data for phylogeny. Except for some of the shortest internodes, phylogenetic relationships of the Lilium species inferred from the chloroplast genome obtained high support, indicating that chloroplast genome data will be useful to help resolve the deeper branches of phylogeny.