Most chloroplast DNA of maize seedlings in linear molecules with defined ends and branched forms

Plastid Inheritance Revisited: Emerging Role of Organelle DNA Degradation in Angiosperms

Plastids are essential organelles in angiosperms and show non-Mendelian inheritance due to their evolution as endosymbionts. In approximately 80% of angiosperms, plastids are thought to be inherited from the maternal parent, whereas other species transmit plastids biparentally. Maternal inheritance can be generally explained by the stochastic segregation of maternal plastids after fertilization because the zygote is overwhelmed by the maternal cytoplasm. In contrast, biparental inheritance shows the transmission of organelles from both parents. In some species, maternal inheritance is not absolute and paternal leakage occurs at a very low frequency (∼10-5). A key process controlling the inheritance mode lies in the behavior of plastids during male gametophyte (pollen) development, with accumulating evidence indicating that the plastids themselves or their DNAs are eliminated during pollen maturation or at fertilization. Cytological observations in numerous angiosperm species have revealed several critical steps that mutually influence the degree of plastid transmission quantitatively among different species. This review revisits plastid inheritance from a mechanistic viewpoint. Particularly, we focus on a recent finding demonstrating that both low temperature and plastid DNA degradation mediated by the organelle exonuclease DEFECTIVE IN POLLEN ORGANELLE DNA DEGRADATION1 (DPD1) influence the degree of paternal leakage significantly in tobacco. Given these findings, we also highlight the emerging role of DPD1 in organelle DNA degradation.

Unprecedented variation pattern of plastid genomes and the potential role in adaptive evolution in Poales

BACKGROUND

The plastid is the photosynthetic organelle in plant cell, and the plastid genomes (plastomes) are generally conserved in evolution. As one of the most economically and ecologically important order of angiosperms, Poales was previously documented to exhibit great plastomic variation as an order of photoautotrophic plants.

RESULTS

We acquired 93 plastomes, representing all the 16 families and 5 major clades of Poales to reveal the extent of their variation and evolutionary pattern. Extensive variation including the largest one in monocots with 225,293 bp in size, heterogeneous GC content, and a wide variety of gene duplication and loss were revealed. Moreover, rare occurrences of three inverted repeat (IR) copies in angiosperms and one IR loss were observed, accompanied by short IR (sIR) and small direct repeat (DR). Widespread structural heteroplasmy, diversified inversions, and unusual genomic rearrangements all appeared in Poales, occasionally within a single species. Extensive repeats in the plastomes were found to be positively correlated with the observed inversions and rearrangements. The variation all showed a "small-large-moderate" trend along the evolution of Poales, as well as for the sequence substitution rate. Finally, we found some positively selected genes, mainly in C4 lineages, while the closely related lineages of those experiencing gene loss tended to have undergone more relaxed purifying selection.

CONCLUSIONS

The variation of plastomes in Poales may be related to its successful diversification into diverse habitats and multiple photosynthetic pathway transitions. Our order-scale analyses revealed unusual evolutionary scenarios for plastomes in the photoautotrophic order of Poales and provided new insights into the plastome evolution in angiosperms as a whole.

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.

Chloroplast Genomes and Phylogenetic Analysis of Three Carthamus (Asteraceae) Species

The genus Carthamus Linnaeus, which belongs to the tribe Cardueae in the Asteraceae family, originated in the Mediterranean region and consists of approximately 20 species worldwide. Understanding the phylogeny of the Carthamus is crucial for the cultivation of C. tinctorius. Although chloroplast genomes are widely used for species identification and evolutionary studies, there have been limited investigations on the chloroplast genomes of Carthamus species. In this study, we assembled the chloroplast genomes of C. persicus, C. tinctorius × C. persicus, and C. lanatus and combined them with the five chloroplast genomes of C. tinctorius for comparative genomic analysis. The sizes of the chloroplast genomes of C. lanatus, C. persicus, and C. tinctorius × C. persicus were 152,602 bp, 153,177 bp, and 153,177 bp, respectively. Comparative analysis showed that the chloroplast genome structures of the four Carthamus species were highly conserved. Additionally, the phylogenomic analysis demonstrated that the plastid genome and angiosperms353 dataset significantly improved the phylogenetic support of Carthamus species. This analysis supported Carthamus as a monophyletic taxon and its internal division into the sect. Carthamus and sect. Atractylis. The Carthamus was closely related to Carduncellus, Femeniasia, Phonus, and Centaurea. In conclusion, this study not only expands our understanding of the cp genomes of Carthamus species but also provides support for more comprehensive phylogenetic studies of Carthamus.

Chloroplast genome sequence of Nothapodytes nimmoniana (Icacinaceae), a vulnerable reservoir of camptothecin from Western Ghats

Nothapodytes nimmoniana belongs to family Icacinaceae and is a major source of compound Camptothecin. The global demand for Camptothecin has caused large-scale exploitation of N. nimmoniana from its wild habitat in Western Ghats of India, thereby making it vulnerable. The species is known to exhibit genetic diversity among the populations in Western Ghats. In this study, we report plastome sequence of N. nimmoniana, first for the genus. For the study, the species was collected from Western Ghats of Maharashtra. The plastome of N. nimmoniana was 150,726 bp in length and exhibited typical quadripartite structure with 83,771 bp LSC, 18,513 bp SSC and 24,221 IR region. The plastome was characterized by presence of 124 unique genes, 87 protein coding genes, 29 tRNA genes and four rRNA genes. Further, the plastome was compared with the available basal lamiid plastomes for gene order and composition. N. nimmoniana plastome exhibited SSC region in an inverted configuration. Phylogenomic study placed N. nimmoniana sister to Mappia mexicana. The SSR markers identified in this study, might help to distinguish genetically diverse populations, prioritizing the populations which need immediate conservation effects as well as for checking adulteration.

Loss of plastid ndh genes in an autotrophic desert plant

Plant plastid genomes are highly conserved with most flowering plants having the same complement of essential plastid genes. Here, we report the loss of five of the eleven NADH dehydrogenase subunit genes (ndh) in the plastid of a desert plant jojoba (Simmondsia chinensis). The plastid genome of jojoba was 156,496 bp with one large single copy region (LSC), a very small single copy region (SSC) and two expanded inverted repeats (IRA + IRB). The NADH dehydrogenase (NDH) complex is comprised of several protein subunits, encoded by the ndh genes of the plastome and the nucleus. The ndh genes are critical to the proper functioning of the photosynthetic electron transport chain and protection of plants from oxidative stress. Most plants are known to contain all eleven ndh genes. Plants with missing or defective ndh genes are often heterotrophs either due to their complete or holo- or myco- parasitic nature. Plants with a defective NDH complex, caused by the deletion/pseudogenisation of some or all the ndh genes, survive in milder climates suggesting the likely extinction of plant lineages lacking these genes under harsh climates. Interestingly, some autotrophic plants do exist without ndh gene/s and can cope with high or low light. This implies that these plants are protected from oxidative stress by mechanisms excluding ndh genes. Jojoba has evolved mechanisms to cope with a non-functioning NDH complex and survives in extreme desert conditions with abundant sunlight and limited water.

Multiple configurations of the plastid and mitochondrial genomes of Caragana spinosa

MAIN CONCLUSION

In this study, we assembled the complete plastome and mitogenome of Caragana spinosa and explored the multiple configurations of the organelle genomes. Caragana spinosa belongs to the Papilionoidea subfamily and has significant pharmaceutical value. To explore the possible interaction between the organelle genomes, we assembled and analyzed the plastome and mitogenome of C. spinosa using the Illumina and Nanopore DNA sequencing data. The plastome of C. spinosa was 129,995 bp belonging to the inverted repeat lacking clade (IRLC), which contained 77 protein-coding genes, 29 tRNA genes, and four rRNA genes. The mitogenome was 378,373 bp long and encoded 54 unique genes, including 33 protein-coding, three ribosomal RNA (rRNA), and 18 transfer RNA (tRNA) genes. In addition to the single circular conformation, alternative conformations mediated by one and four repetitive sequences in the plastome and mitogenome were identified and validated, respectively. The inverted repeat (PDR12, the 12th dispersed repeat sequence in C. spinosa plastome) of plastome mediating recombinant was conserved in the genus Caragana. Furthermore, we identified 14 homologous fragments by comparing the sequences of mitogenome and plastome, including eight complete tRNA genes. A phylogenetic analysis of protein-coding genes extracted from the plastid and mitochondrial genomes revealed congruent topologies. Analyses of sequence divergence found one intergenic region, trnN-GUU-ycf1, exhibiting a high degree of variation, which can be used to develop novel molecular markers to distinguish the nine Caragana species accurately. This plastome and mitogenome of C. spinosa could provide critical information for the molecular breeding of C. spinosa and be used as a reference genome for other species of Caragana. In this study, we assembled the complete plastome and mitogenome of Caragana spinosa and explored the multiple configurations of the organelle genomes.

Oxidative and Glycation Damage to Mitochondrial DNA and Plastid DNA during Plant Development

Oxidative damage to plant proteins, lipids, and DNA caused by reactive oxygen species (ROS) has long been studied. The damaging effects of reactive carbonyl groups (glycation damage) to plant proteins and lipids have also been extensively studied, but only recently has glycation damage to the DNA in plant mitochondria and plastids been reported. Here, we review data on organellar DNA maintenance after damage from ROS and glycation. Our focus is maize, where tissues representing the entire range of leaf development are readily obtained, from slow-growing cells in the basal meristem, containing immature organelles with pristine DNA, to fast-growing leaf cells, containing mature organelles with highly-fragmented DNA. The relative contributions to DNA damage from oxidation and glycation are not known. However, the changing patterns of damage and damage-defense during leaf development indicate tight coordination of responses to oxidation and glycation events. Future efforts should be directed at the mechanism by which this coordination is achieved.

Plastid phylogenomics and molecular evolution of Thismiaceae (Dioscoreales)

PREMISE

Species in Thismiaceae can no longer photosynthesize and instead obtain carbon from soil fungi. Here we infer Thismiaceae phylogeny using plastid genome data and characterize the molecular evolution of this genome.

METHODS

We assembled five Thismiaceae plastid genomes from genome skimming data, adding to previously published data for phylogenomic inference. We investigated plastid-genome structural changes, considering locally colinear blocks (LCBs). We also characterized possible shifts in selection pressure in retained genes by considering changes in the ratio of nonsynonymous to synonymous changes (ω).

RESULTS

Thismiaceae experienced two major pulses of gene loss around the early diversification of the family, with subsequent scattered gene losses across descendent lineages. In addition to massive size reduction, Thismiaceae plastid genomes experienced occasional inversions, and there were likely two independent losses of the plastid inverted repeat (IR) region. Retained plastid genes remain under generally strong purifying selection (ω << 1), with significant and sporadic weakening or strengthening in several instances. The bifunctional trnE-UUC gene of Thismia huangii may retain a secondary role in heme biosynthesis, despite a probable loss of functionality in protein translation. Several cis-spliced group IIA introns have been retained, despite the loss of the plastid intron maturase, matK.

CONCLUSIONS

We infer that most gene losses in Thismiaceae occurred early and rapidly, following the initial loss of photosynthesis in its stem lineage. As a species-rich, fully mycoheterotrophic lineage, Thismiaceae provide a model system for uncovering the unique and divergent ways in which plastid genomes evolve in heterotrophic plants.

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

PREMISE

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

METHODS

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

RESULTS

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

CONCLUSIONS

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

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

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

Analysis of the Plastid Genome Sequence During Maize Seedling Development

Shoot development in maize progresses from small, non-pigmented meristematic cells to expanded cells in the green leaf. During this transition, large plastid DNA (ptDNA) molecules in proplastids become fragmented in the photosynthetically-active chloroplasts. The genome sequences were determined for ptDNA obtained from Zea mays B73 plastids isolated from four tissues: base of the stalk (the meristem region); fully-developed first green leaf; first three leaves from light-grown seedlings; and first three leaves from dark-grown (etiolated) seedlings. These genome sequences were then compared to the Z. mays B73 plastid reference genome sequence that was previously obtained from green leaves. The assembled plastid genome was identical among these four tissues to the reference genome. Furthermore, there was no difference among these tissues in the sequence at and around the previously documented 27 RNA editing sites. There were, however, more sequence variants (insertions/deletions and single-nucleotide polymorphisms) for leaves grown in the dark than in the light. These variants were tightly clustered into two areas within the inverted repeat regions of the plastid genome. We propose a model for how these variant clusters could be generated by replication-transcription conflict.

Complete pan-plastome sequences enable high resolution phylogenetic classification of sugar beet and closely related crop wild relatives

Background

As the major source of sugar in moderate climates, sugar-producing beets (Beta vulgaris subsp. vulgaris) have a high economic value. However, the low genetic diversity within cultivated beets requires introduction of new traits, for example to increase their tolerance and resistance attributes – traits that often reside in the crop wild relatives. For this, genetic information of wild beet relatives and their phylogenetic placements to each other are crucial. To answer this need, we sequenced and assembled the complete plastome sequences from a broad species spectrum across the beet genera Beta and Patellifolia, both embedded in the Betoideae (order Caryophyllales). This pan-plastome dataset was then used to determine the wild beet phylogeny in high-resolution.

Results

We sequenced the plastomes of 18 closely related accessions representing 11 species of the Betoideae subfamily and provided high-quality plastome assemblies which represent an important resource for further studies of beet wild relatives and the diverse plant order Caryophyllales. Their assembly sizes range from 149,723 bp (Beta vulgaris subsp. vulgaris) to 152,816 bp (Beta nana), with most variability in the intergenic sequences. Combining plastome-derived phylogenies with read-based treatments based on mitochondrial information, we were able to suggest a unified and highly confident phylogenetic placement of the investigated Betoideae species.

Our results show that the genus Beta can be divided into the two clearly separated sections Beta and Corollinae. Our analysis confirms the affiliation of B. nana with the other Corollinae species, and we argue against a separate placement in the Nanae section. Within the Patellifolia genus, the two diploid species Patellifolia procumbens and Patellifolia webbiana are, regarding the plastome sequences, genetically more similar to each other than to the tetraploid Patellifolia patellaris. Nevertheless, all three Patellifolia species are clearly separated.

Conclusion

In conclusion, our wild beet plastome assemblies represent a new resource to understand the molecular base of the beet germplasm. Despite large differences on the phenotypic level, our pan-plastome dataset is highly conserved. For the first time in beets, our whole plastome sequences overcome the low sequence variation in individual genes and provide the molecular backbone for highly resolved beet phylogenomics. Hence, our plastome sequencing strategy can also guide genomic approaches to unravel other closely related taxa.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-022-08336-8.

Plastome structure, evolution, and phylogeny of Selaginella

As one of the earliest land plant lineages, Selaginella is important for studying land plant evolution. It is the largest genus of lycophytes containing 700-800 species. Some unique characters of Selaginella plastomes have been reported, but based only on 20 species. There have been no plastome phylogenies of Selaginella based on a relatively large sampling, and no efforts have been made to resolve the phylogeny of the enigmatic Sinensis group whose relationships have been unclear based on small datasets. Here we investigated the structures of 59 plastomes representing 51 species covering all six subgenera and 18 sections of Selaginella except two sections and including the intriguing Sinensis group for the first time. Our major results include: (1) the plastome size of Selaginella ranges tremendously from 78,492 bp to 187,632 bp; (2) there are numerous gene losses in Selaginella comparing with other lycophytes, Isoëtaceae and Lycopodiaceae; (3) the gene contents and plastome structures in Selaginella vary lineage-specifically and all infrageneric taxa are well supported in the plastome phylogeny; (4) the ndh gene family tends to lose or pseudogenize in those species with DR structure and without other short or medium repeats; (5) the short and medium repeat regions in SC mediate many conformations causing diverse and complex plastome structures, and six new conformations are discovered; (6) forty-eight species sampled have high GC content (>50%) but three species in the Sinensis group have ∼30% GC content in plastomes, similar to most vascular plants; (7) the Sinensis group is monophyletic, includes at least two subgroups, and has the smallest plastomes in land plants except some parasitic plants, and their plastomes do not contain any tRNAs; (8) the younger lineages in Selaginella tend to have higher GC content, whereas the older lineages tend to have lower GC content; and (9) because of incomplete genomic data and abnormal structures or some unknown reasons, even the concatenated plastomes could not well resolve the phylogenetic relationships in Selaginella with confidence, highlighting the difficulty in resolving the phylogeny and evolution of this particularly important land plant lineage.

Calophyllaceae plastomes, their structure and insights in relationships within the clusioids

A complete chloroplast genome is not yet available for numerous species of plants. Among the groups that lack plastome information is the clusioid clade (Malpighiales), which includes five families: Bonnetiaceae, Calophyllaceae, Clusiaceae, Hypericaceae, and Podostemaceae. With around 2200 species, it has few published plastomes and most of them are from Podostemaceae. Here we assembled and compared six plastomes from members of the clusioids: five from Calophyllaceae (newly sequenced) and one from Clusiaceae. Putative regions for evolutionary studies were identified and the newly assembled chloroplasts were analyzed with other available chloroplasts for the group, focusing on Calophyllaceae. Our results mostly agree with recent studies which found a general conserved structure, except for the two Podostemaceae species that have a large inversion (trnK-UUU–rbcL) and lack one intron from ycf3. Within Calophyllaceae we observed a longer LSC and reduced IRs in Mahurea exstipulata, resulting in some genic rearrangement, and a short inversion (psbJ–psbE) in Kielmeyera coriacea. Phylogenetic analyses recovered the clusioids and the five families as monophyletic and revealed that conflicts in relationships reported in the literature for the group agree with nodes concentrating uninformative or conflicting gene trees. Our study brings new insights about clusioid plastome architecture and its evolution.

Plastome Structural Evolution and Homoplastic Inversions in Neo-Astragalus (Fabaceae)

The plastid genomes of photosynthetic green plants have largely maintained conserved gene content and order as well as structure over hundreds of millions of years of evolution. Several plant lineages, however, have departed from this conservation and contain many plastome structural rearrangements, which have been associated with an abundance of repeated sequences both overall and near rearrangement endpoints. We sequenced the plastomes of 25 taxa of Astragalus L. (Fabaceae), a large genus in the inverted repeat-lacking clade of legumes, to gain a greater understanding of the connection between repeats and plastome inversions. We found plastome repeat structure has a strong phylogenetic signal among these closely related taxa mostly in the New World clade of Astragalus called Neo-Astragalus. Taxa without inversions also do not differ substantially in their overall repeat structure from four taxa each with one large-scale inversion. For two taxa with inversion endpoints between the same pairs of genes, differences in their exact endpoints indicate the inversions occurred independently. Our proposed mechanism for inversion formation suggests the short inverted repeats now found near the endpoints of the four inversions may be there as a result of these inversions rather than their cause. The longer inverted repeats now near endpoints may have allowed the inversions first mediated by shorter microhomologous sequences to propagate, something that should be considered in explaining how any plastome rearrangement becomes fixed regardless of the mechanism of initial formation.

The chicken or the egg? Plastome evolution and an independent loss of the inverted repeat in papilionoid legumes

The plastid genome (plastome), while surprisingly constant in gene order and content across most photosynthetic angiosperms, exhibits variability in several unrelated lineages. During the diversification history of the legume family Fabaceae, plastomes have undergone many rearrangements, including inversions, expansion, contraction and loss of the typical inverted repeat (IR), gene loss and repeat accumulation in both shared and independent events. While legume plastomes have been the subject of study for some time, most work has focused on agricultural species in the IR-lacking clade (IRLC) and the plant model Medicago truncatula. The subfamily Papilionoideae, which contains virtually all of the agricultural legume species, also comprises most of the plastome variation detected thus far in the family. In this study three non-papilioniods were included among 34 newly sequenced legume plastomes, along with 33 publicly available sequences, to assess plastome structural evolution in the subfamily. In an effort to examine plastome variation across the subfamily, approximately 20% of the sampling represents the IRLC with the remainder selected to represent the early-branching papilionoid clades. A number of IR-related and repeat-mediated changes were identified and examined in a phylogenetic context. Recombination between direct repeats associated with ycf2 resulted in intraindividual plastome heteroplasmy. Although loss of the IR has not been reported in legumes outside of the IRLC, one genistoid taxon was found to completely lack the typical plastome IR. The role of the IR and non-IR repeats in the progression of plastome change is discussed.

Unprecedented Intraindividual Structural Heteroplasmy in Eleocharis (Cyperaceae, Poales) Plastomes

Plastid genomes (plastomes) of land plants have a conserved quadripartite structure in a gene-dense unit genome consisting of a large inverted repeat that separates two single copy regions. Recently, alternative plastome structures were suggested in Geraniaceae and in some conifers and Medicago the coexistence of inversion isomers has been noted. In this study, plastome sequences of two Cyperaceae, Eleocharis dulcis (water chestnut) and Eleocharis cellulosa (gulf coast spikerush), were completed. Unlike the conserved plastomes in basal groups of Poales, these Eleocharis plastomes have remarkably divergent features, including large plastome sizes, high rates of sequence rearrangements, low GC content and gene density, gene duplications and losses, and increased repetitive DNA sequences. A novel finding among these features was the unprecedented level of heteroplasmy with the presence of multiple plastome structural types within a single individual. Illumina paired-end assemblies combined with PacBio single-molecule real-time sequencing, long-range polymerase chain reaction, and Sanger sequencing data identified at least four different plastome structural types in both Eleocharis species. PacBio long read data suggested that one of the four E. dulcis plastome types predominates.

Complete chloroplast genomes of Achnatherum inebrians and comparative analyses with related species from Poaceae

This article reports the complete chloroplast genome of Achnatherum inebrians, a poisonous herb that is widely distributed in the rangelands of Northern China. The genome is 137 714 bp in total and consists of a large single‐copy (81 758 bp) region and small single‐copy (12 682 bp) region separated by a pair of inverted repeats (21 637 bp). The genome contains 130 genes, including 84 protein‐coding genes, 38 tRNA genes and 8 ribosomal RNA genes, and the guanine + cytosine content is 36.17%. We subsequently performed comparative analysis of complete genomes from A. inebrians and other Poaceae‐related species from GenBank. Thirty‐eight simple sequence repeats were identified, further demonstrating rapid evolution in Poaceae. Finally, the phylogenetic trees of 37 species of Poaceae and 2 species of Amaranthaceae were constructed by using maximum likelihood and Bayesian inference methods, based on the genes of the complete chloroplast genome. We identified hotspots that can be used as molecular markers and barcodes for phylogenetic analysis, as well as for species identification. Phylogenetic analysis indicated that A. inebrians is a member of the genus Stipa rather than Achnatherum.

Advancing organelle genome transformation and editing for crop improvement

Plant cells contain three organelles that harbor DNA: the nucleus, plastids, and mitochondria. Plastid transformation has emerged as an attractive platform for the generation of transgenic plants, also referred to as transplastomic plants. Plastid genomes have been genetically engineered to improve crop yield, nutritional quality, and resistance to abiotic and biotic stresses, as well as for recombinant protein production. Despite many promising proof-of-concept applications, transplastomic plants have not been commercialized to date. Sequence-specific nuclease technologies are widely used to precisely modify nuclear genomes, but these tools have not been applied to edit organelle genomes because the efficient homologous recombination system in plastids facilitates plastid genome editing. Unlike plastid transformation, successful genetic transformation of higher plant mitochondrial genome transformation was tested in several research group, but not successful to date. However, stepwise progress has been made in modifying mitochondrial genes and their transcripts, thus enabling the study of their functions. Here, we provide an overview of advances in organelle transformation and genome editing for crop improvement, and we discuss the bottlenecks and future development of these technologies.

RHON1 Co-transcriptionally Resolves R-Loops for Arabidopsis Chloroplast Genome Maintenance

Preventing transcription-replication head-on conflict (HO-TRC)-triggered R-loop formation is essential for maintaining genome integrity in bacteria, plants, and mammals. The R-loop eraser RNase H can efficiently relax HO-TRCs. However, it is not clear how organisms resist HO-TRC-triggered R-loops when RNase H proteins are deficient. By screening factors that may relieve R-loop accumulation in the Arabidopsis atrnh1c mutant, we find that overexpression of the R-loop helicase RHON1 can rescue the defects of aberrantly accumulated HO-TRC-triggered R-loops co-transcriptionally. In addition, we find that RHON1 interacts with and orchestrates the transcriptional activity of plastid-encoded RNA polymerases to release the conflicts between transcription and replication. Our study illustrates that organisms employ multiple mechanisms to escape HO-TRC-triggered R-loop accumulation and thus maintain genome integrity.

Plastid Genomes of Flowering Plants: Essential Principles

The plastid genome (plastome ) has proved a valuable source of data for evaluating evolutionary relationships among angiosperms. Through basic and applied approaches, plastid transformation technology offers the potential to understand and improve plant productivity, providing food, fiber, energy, and medicines to meet the needs of a burgeoning global population. The growing genomic resources available to both phylogenetic and biotechnological investigations is allowing novel insights and expanding the scope of plastome research to encompass new species. In this chapter, we present an overview of some of the seminal and contemporary research that has contributed to our current understanding of plastome evolution and attempt to highlight the relationship between evolutionary mechanisms and the tools of plastid genetic engineering.

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

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

Plant Organellar DNA Polymerases Evolved Multifunctionality through the Acquisition of Novel Amino Acid Insertions

The majority of DNA polymerases (DNAPs) are specialized enzymes with specific roles in DNA replication, translesion DNA synthesis (TLS), or DNA repair. The enzymatic characteristics to perform accurate DNA replication are in apparent contradiction with TLS or DNA repair abilities. For instance, replicative DNAPs incorporate nucleotides with high fidelity and processivity, whereas TLS DNAPs are low-fidelity polymerases with distributive nucleotide incorporation. Plant organelles (mitochondria and chloroplast) are replicated by family-A DNA polymerases that are both replicative and TLS DNAPs. Furthermore, plant organellar DNA polymerases from the plant model Arabidopsis thaliana (AtPOLIs) execute repair of double-stranded breaks by microhomology-mediated end-joining and perform Base Excision Repair (BER) using lyase and strand-displacement activities. AtPOLIs harbor three unique insertions in their polymerization domain that are associated with TLS, microhomology-mediated end-joining (MMEJ), strand-displacement, and lyase activities. We postulate that AtPOLIs are able to execute those different functions through the acquisition of these novel amino acid insertions, making them multifunctional enzymes able to participate in DNA replication and DNA repair.

Reference Genome for the Highly Transformable Setaria viridis ME034V

Setaria viridis (green foxtail) is an important model system for improving cereal crops due to its diploid genome, ease of cultivation, and use of C4 photosynthesis. The S. viridis accession ME034V is exceptionally transformable, but the lack of a sequenced genome for this accession has limited its utility. We present a 397 Mb highly contiguous de novo assembly of ME034V using ultra-long nanopore sequencing technology (read N50 = 41kb). We estimate that this genome is largely complete based on our updated k-mer based genome size estimate of 401 Mb for S. viridis Genome annotation identified 37,908 protein-coding genes and >300k repetitive elements comprising 46% of the genome. We compared the ME034V assembly with two other previously sequenced Setaria genomes as well as to a diversity panel of 235 S. viridis accessions. We found the genome assemblies to be largely syntenic, but numerous unique polymorphic structural variants were discovered. Several ME034V deletions may be associated with recent retrotransposition of copia and gypsy LTR repeat families, as evidenced by their low genotype frequencies in the sampled population. Lastly, we performed a phylogenomic analysis to identify gene families that have expanded in Setaria, including those involved in specialized metabolism and plant defense response. The high continuity of the ME034V genome assembly validates the utility of ultra-long DNA sequencing to improve genetic resources for emerging model organisms. Structural variation present in Setaria illustrates the importance of obtaining the proper genome reference for genetic experiments. Thus, we anticipate that the ME034V genome will be of significant utility for the Setaria research community.

Methods of analysis of chloroplast genomes of C3, Kranz type C4 and Single Cell C4 photosynthetic members of Chenopodiaceae

BACKGROUND

Chloroplast genome information is critical to understanding forms of photosynthesis in the plant kingdom. During the evolutionary process, plants have developed different photosynthetic strategies that are accompanied by complementary biochemical and anatomical features. Members of family Chenopodiaceae have species with C3 photosynthesis, and variations of C4 photosynthesis in which photorespiration is reduced by concentrating CO2 around Rubisco through dual coordinated functioning of dimorphic chloroplasts. Among dicots, the family has the largest number of C4 species, and greatest structural and biochemical diversity in forms of C4 including the canonical dual-cell Kranz anatomy, and the recently identified single cell C4 with the presence of dimorphic chloroplasts separated by a vacuole. This is the first comparative analysis of chloroplast genomes in species representative of photosynthetic types in the family.

RESULTS

Methodology with high throughput sequencing complemented with Sanger sequencing of selected loci provided high quality and complete chloroplast genomes of seven species in the family and one species in the closely related Amaranthaceae family, representing C3, Kranz type C4 and single cell C4 (SSC4) photosynthesis six of the eight chloroplast genomes are new, while two are improved versions of previously published genomes. The depth of coverage obtained using high-throughput sequencing complemented with targeted resequencing of certain loci enabled superior resolution of the border junctions, directionality and repeat region sequences. Comparison of the chloroplast genomes with previously sequenced plastid genomes revealed similar genome organization, gene order and content with a few revisions. High-quality complete chloroplast genome sequences resulted in correcting the orientation the LSC region of the published Bienertia sinuspersici chloroplast genome, identification of stop codons in the rpl23 gene in B. sinuspersici and B. cycloptera, and identifying an instance of IR expansion in the Haloxylon ammodendron inverted repeat sequence. The rare observation of a mitochondria-to-chloroplast inter-organellar gene transfer event was identified in family Chenopodiaceae.

CONCLUSIONS

This study reports complete chloroplast genomes from seven Chenopodiaceae and one Amaranthaceae species. The depth of coverage obtained using high-throughput sequencing complemented with targeted resequencing of certain loci enabled superior resolution of the border junctions, directionality, and repeat region sequences. Therefore, the use of high throughput and Sanger sequencing, in a hybrid method, reaffirms to be rapid, efficient, and reliable for chloroplast genome sequencing.

Under the rug: Abandoning persistent misconceptions that obfuscate organelle evolution

The advent and advance of next generation sequencing over the past two decades made it possible to accumulate large quantities of sequence reads that could be used to assemble complete or nearly complete organelle genomes (plastome or mitogenome). The result has been an explosive increase in the availability of organelle genome sequences with over 4000 different species of green plants currently available on GenBank. During the same time period, plant molecular biologists greatly enhanced the understanding of the structure, repair, replication, recombination, transcription and translation, and inheritance of organelle DNA. Unfortunately many plant evolutionary biologists are unaware of or have overlooked this knowledge, resulting in misrepresentation of several phenomena that are critical for phylogenetic and evolutionary studies using organelle genomes. We believe that confronting these misconceptions about organelle genome organization, composition, and inheritance will improve our understanding of the evolutionary processes that underly organelle evolution. Here we discuss four misconceptions that can limit evolutionary biology studies and lead to inaccurate phylogenies and incorrect structure of the organellar DNA used to infer organelle evolution.

Long-Reads Reveal That the Chloroplast Genome Exists in Two Distinct Versions in Most Plants

The chloroplast genome usually has a quadripartite structure consisting of a large single copy region and a small single copy region separated by two long inverted repeats. It has been known for some time that a single cell may contain at least two structural haplotypes of this structure, which differ in the relative orientation of the single copy regions. However, the methods required to detect and measure the abundance of the structural haplotypes are labor-intensive, and this phenomenon remains understudied. Here, we develop a new method, Cp-hap, to detect all possible structural haplotypes of chloroplast genomes of quadripartite structure using long-read sequencing data. We use this method to conduct a systematic analysis and quantification of chloroplast structural haplotypes in 61 land plant species across 19 orders of Angiosperms, Gymnosperms, and Pteridophytes. Our results show that there are two chloroplast structural haplotypes which occur with equal frequency in most land plant individuals. Nevertheless, species whose chloroplast genomes lack inverted repeats or have short inverted repeats have just a single structural haplotype. We also show that the relative abundance of the two structural haplotypes remains constant across multiple samples from a single individual plant, suggesting that the process which maintains equal frequency of the two haplotypes operates rapidly, consistent with the hypothesis that flip-flop recombination mediates chloroplast structural heteroplasmy. Our results suggest that previous claims of differences in chloroplast genome structure between species may need to be revisited.

Plant organellar DNA polymerases repair double-stranded breaks by microhomology-mediated end-joining

Double-stranded breaks (DSBs) in plant organelles are repaired via genomic rearrangements characterized by microhomologous repeats. These microhomologous signatures predict the existence of an unidentified enzymatic machinery capable of repairing of DSBs via microhomology-mediated end-joining (MMEJ) in plant organelles. Here, we show that organellar DNA polymerases from Arabidopsis thaliana (AtPolIA and AtPolIB) perform MMEJ using microhomologous sequences as short as six nucleotides. AtPolIs execute MMEJ by virtue of two specialized amino acid insertions located in their thumb subdomains. Single-stranded binding proteins (SSBs) unique to plants, AtWhirly2 and organellar single-stranded binding proteins (AtOSBs), hinder MMEJ, whereas canonical mitochondrial SSBs (AtmtSSB1 and AtmtSSB2) do not interfere with MMEJ. Our data predict that organellar DNA rearrangements by MMEJ are a consequence of a competition for the 3'-OH of a DSBs. If AtWhirlies or AtOSBs gain access to the single-stranded DNA (ssDNA) region of a DSB, the reaction will shift towards high-fidelity routes like homologous recombination. Conversely MMEJ would be favored if AtPolIs or AtmtSSBs interact with the DSB. AtPolIs are not phylogenetically related to metazoan mitochondrial DNA polymerases, and the ability of AtPolIs to execute MMEJ may explain the abundance of DNA rearrangements in plant organelles in comparison to animal mitochondria.

Chloroplast genome sequences of Artemisia maritima and Artemisia absinthium: Comparative analyses, mutational hotspots in genus Artemisia and phylogeny in family Asteraceae

Artemisia L. is a complex genus of medicinal importance. Publicly available chloroplast genomes of few Artemisia species are insufficient to resolve taxonomic discrepancies at species level. We report chloroplast genome sequences of two further Artemisia species: A. maritima (151,061 bp) and A. absinthium (151,193 bp). Both genomes possess typical quadripartite structure comprising of a large single copy, a small single copy and a pair of long inverted repeats. The two genomes exhibited high similarities in genome sizes, gene synteny, GC content, synonymous and non-synonymous substitutions, codon usage, amino acids frequencies, RNA editing sites, microsatellites, and oligonucleotide repeats. Transition to transversion ratio was <1. Maximum likelihood tree showed Artemisia a monophyletic genus, sister to genus Chrysanthemum. We also identified 20 highly polymorphic regions including rpoC2-rps2, trnR-UCU-trnG-UCC, rps18-rpl20, and trnL-UAG-rpl32 that could be used to develop authentic and cost-effective markers to resolve taxonomic discrepancies and infer phylogenetic relationships among Artemisia species.

Integrated Genomic and Transcriptomic Analysis of the Peridinin Dinoflagellate Amphidinium carterae Plastid

The plastid genomes of peridinin-containing dinoflagellates are highly unusual, possessing very few genes, which are located on small chromosomal elements termed "minicircles". These minicircles may contain genes, or no recognisable coding information. Transcripts produced from minicircles may undergo unusual processing events, such as the addition of a 3' poly(U) tail. To date, little is known about the genetic or transcriptional diversity of non-coding sequences in peridinin dinoflagellate plastids. These sequences include empty minicircles, and regions of non-coding DNA in coding minicircles. Here, we present an integrated plastid genome and transcriptome for the model peridinin dinoflagellate Amphidinium carterae, identifying a previously undescribed minicircle. We also profile transcripts covering non-coding regions of the psbA and petB/atpA minicircles. We present evidence that antisense transcripts are produced within the A. carterae plastid, but show that these transcripts undergo different end cleavage events from sense transcripts, and do not receive 3' poly(U) tails. The difference in processing events between sense and antisense transcripts may enable the removal of non-coding transcripts from peridinin dinoflagellate plastid transcript pools.

Lost and Found: Return of the Inverted Repeat in the Legume Clade Defined by Its Absence

The plant genome comprises a coevolving, integrated genetic system housed in three subcellular compartments: the nucleus, mitochondrion, and the plastid. The typical land plant plastid genome (plastome) comprises the sum of repeating units of 130–160 kb in length. The plastome inverted repeat (IR) divides each plastome monomer into large and small single copy regions, an architecture highly conserved across land plants. There have been varying degrees of expansion or contraction of the IR, and in a few distinct lineages, including the IR-lacking clade of papilionoid legumes, one copy of the IR has been lost. Completion of plastome sequencing and assembly for 19 Medicago species and Trigonella foenum-graceum and comparative analysis with other IR-lacking clade taxa revealed modest divergence with regard to structural organization overall. However, one clade contained unique variation suggesting an ancestor had experienced repeat-mediated changes in plastome structure. In Medicago minima, a novel IR of ∼9 kb was confirmed and the role of repeat-mediated, recombination-dependent replication in IR reemergence is discussed.

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

BACKGROUND

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

RESULTS

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

CONCLUSIONS

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

Tension and Resolution: Dynamic, Evolving Populations of Organelle Genomes within Plant Cells

Mitochondria and plastids form dynamic, evolving populations physically embedded in the fluctuating environment of the plant cell. Their evolutionary heritage has shaped how the cell controls the genetic structure and the physical behavior of its organelle populations. While the specific genes involved in these processes are gradually being revealed, the governing principles underlying this controlled behavior remain poorly understood. As the genetic and physical dynamics of these organelles are central to bioenergetic performance and plant physiology, this challenges both fundamental biology and strategies to engineer better-performing plants. This article reviews current knowledge of the physical and genetic behavior of mitochondria and chloroplasts in plant cells. An overarching hypothesis is proposed whereby organelles face a tension between genetic robustness and individual control and responsiveness, and different species resolve this tension in different ways. As plants are immobile and thus subject to fluctuating environments, their organelles are proposed to favor individual responsiveness, sacrificing genetic robustness. Several notable features of plant organelles, including large genomes, mtDNA recombination, fragmented organelles, and plastid/mitochondrial differences may potentially be explained by this hypothesis. Finally, the ways that quantitative and systems biology can help shed light on the plethora of open questions in this field are highlighted.

Highly accelerated rates of genomic rearrangements and nucleotide substitutions in plastid genomes of Passiflora subgenus Decaloba

Plastid genomes (plastomes) of photosynthetic angiosperms are for the most part highly conserved in their organization, mode of inheritance and rates of nucleotide substitution. A small number of distantly related lineages share a syndrome of features that deviate from this general pattern, including extensive genomic rearrangements, accelerated rates of nucleotide substitution, biparental inheritance and plastome-genome incompatibility. Previous studies of plastomes in Passiflora with limited taxon sampling suggested that the genus exhibits this syndrome. To examine this phenomenon further, 15 new plastomes from Passiflora were sequenced and combined with previously published data to examine the phylogenetic relationships, genome organization and evolutionary rates across all five subgenera and the sister genus Adenia. Phylogenomic analyses using 68 protein-coding genes shared by Passiflora generated a fully resolved and strongly supported tree that is congruent with previous phylogenies based on a few plastid and nuclear loci. This phylogeny was used to examine the distribution of plastome rearrangements across Passiflora. Multiple gene and intron losses and inversions were identified in Passiflora with some occurring in parallel and others that extended across the Passifloraceae. Furthermore, extensive expansions and contractions of the inverted repeat (IR) were uncovered and in some cases this resulted in exclusion of all ribosomal RNA genes from the IR. The most highly rearranged lineage was subgenus Decaloba, which experienced extensive IR expansion that incorporated up to 25 protein-coding genes usually located in large single copy region. Nucleotide substitution rate analyses of 68 protein-coding genes across the genus showed lineage- and locus-specific acceleration. Significant increase in dS, dN and dN/dS was detected for clpP across the genus and for ycf4 in certain lineages. Significant increases in dN and dN/dS for ribosomal subunits and plastid-encoded RNA polymerase genes were detected in the branch leading to the expanded IR-clade in subgenus Decaloba. This subgenus displays the syndrome of unusual features, making it an ideal system to investigate the dynamic evolution of angiosperm plastomes.

Detecting useful genetic markers and reconstructing the phylogeny of an important medicinal resource plant, Artemisia selengensis, based on chloroplast genomics

Artemisia selengenesis is not only a health food, but also a well-known traditional Chinese medicine. Only a fraction of the chloroplast (cp) genome data of Artemisia has been reported and chloroplast genomic materials have been widely used in genomic evolution studies, molecular marker development, and phylogenetic analysis of the genus Artemisia, which makes evolutionary studies, genetic improvement, and phylogenetic identification very difficult. In this study, the complete chloroplast genome of A. selengensis was compared with that of other species within Artemisia and phylogenetic analyses was conducted with other genera in the Asteraceae family. The results showed that A. selengensis is an AT-rich species and has a typical quadripartite structure that is 151,215 bp in length. Comparative genome analyses demonstrated that the available chloroplast genomes of species of Artemisia were well conserved in terms of genomic length, GC contents, and gene organization and order. However, some differences, which may indicate evolutionary events, were found, such as a re-inversion event within the Artemisia genus, an unequal duplicate phenomenon of the ycf1 gene because of the expansion and contraction of the IR region, and the fast-evolving regions. Repeated sequences analysis showed that Artemisia chloroplast genomes presented a highly similar pattern of SSR or LDR distribution. A total of 257 SSRs and 42 LDRs were identified in the A. selengensis chloroplast genome. The phylogenetic analysis showed that A. selengensis was sister to A. gmelinii. The findings of this study will be valuable in further studies to understand the genetic diversity and evolutionary history of Asteraceae.

Assembly of chloroplast genomes with long- and short-read data: a comparison of approaches using Eucalyptus pauciflora as a test case

BACKGROUND

Chloroplasts are organelles that conduct photosynthesis in plant and algal cells. The information chloroplast genome contained is widely used in agriculture and studies of evolution and ecology. Correctly assembling chloroplast genomes can be challenging because the chloroplast genome contains a pair of long inverted repeats (10-30 kb). Typically, it is simply assumed that the gross structure of the chloroplast genome matches the most commonly observed structure of two single-copy regions separated by a pair of inverted repeats. The advent of long-read sequencing technologies should remove the need to make this assumption by providing sufficient information to completely span the inverted repeat regions. Yet, long-reads tend to have higher error rates than short-reads, and relatively little is known about the best way to combine long- and short-reads to obtain the most accurate chloroplast genome assemblies. Using Eucalyptus pauciflora, the snow gum, as a test case, we evaluated the effect of multiple parameters, such as different coverage of long-(Oxford nanopore) and short-(Illumina) reads, different long-read lengths, different assembly pipelines, with a view to determining the most accurate and efficient approach to chloroplast genome assembly.

RESULTS

Hybrid assemblies combining at least 20x coverage of both long-reads and short-reads generated a single contig spanning the entire chloroplast genome with few or no detectable errors. Short-read-only assemblies generated three contigs (the long single copy, short single copy and inverted repeat regions) of the chloroplast genome. These contigs contained few single-base errors but tended to exclude several bases at the beginning or end of each contig. Long-read-only assemblies tended to create multiple contigs with a much higher single-base error rate. The chloroplast genome of Eucalyptus pauciflora is 159,942 bp, contains 131 genes of known function.

CONCLUSIONS

Our results suggest that very accurate assemblies of chloroplast genomes can be achieved using a combination of at least 20x coverage of long- and short-reads respectively, provided that the long-reads contain at least ~5x coverage of reads longer than the inverted repeat region. We show that further increases in coverage give little or no improvement in accuracy, and that hybrid assemblies are more accurate than long-read-only or short-read-only assemblies.

Method for the Identification of Taxon-Specific k-mers from Chloroplast Genome: A Case Study on Tomato Plant (Solanum lycopersicum)

Polymerase chain reaction and different barcoding methods commonly used for plant identification from metagenomics samples are based on the amplification of a limited number of pre-selected barcoding regions. These methods are often inapplicable due to DNA degradation, low amplification success or low species discriminative power of selected genomic regions. Here we introduce a method for the rapid identification of plant taxon-specific k-mers, that is applicable for the fast detection of plant taxa directly from raw sequencing reads without aligning, mapping or assembling the reads. We identified more than 800 Solanum lycopersicum specific k-mers (32 nucleotides in length) from 42 different chloroplast genome regions using the developed method. We demonstrated that identified k-mers are also detectable in whole genome sequencing raw reads from S. lycopersicum. Also, we demonstrated the usability of taxon-specific k-mers in artificial mixtures of sequences from closely related species. Developed method offers a novel strategy for fast identification of taxon-specific genome regions and offers new perspectives for detection of plant taxa directly from sequencing raw reads.

Global RNA association with the transcriptionally active chromosome of chloroplasts

KEY MESSAGE

Processed chloroplast RNAs are co-enriched with preparations of the chloroplast transcriptionally active chromosome. Chloroplast genomes are organized as a polyploid DNA-protein structure called the nucleoid. Transcriptionally active chloroplast DNA together with tightly bound protein factors can be purified by gel filtration as a functional entity called the transcriptionally active chromosome (TAC). Previous proteomics analyses of nucleoids and of TACs demonstrated a considerable overlap in protein composition including RNA binding proteins. Therefore the RNA content of TAC preparations from Nicotiana tabacum was determined using whole genome tiling arrays. A large number of chloroplast RNAs was found to be associated with the TAC. The pattern of RNAs attached to the TAC consists of RNAs produced by different chloroplast RNA polymerases and differs from the pattern of RNA found in input controls. An analysis of RNA splicing and RNA editing of selected RNA species demonstrated that TAC-associated RNAs are processed to a similar extent as the RNA in input controls. Thus, TAC fractions contain a specific subset of the processed chloroplast transcriptome.

accD nuclear transfer of Platycodon grandiflorum and the plastid of early Campanulaceae

BACKGROUND

Campanulaceae species are known to have highly rearranged plastid genomes lacking the acetyl-CoA carboxylase (ACC) subunit D gene (accD), and instead have a nuclear (nr)-accD. Plastid genome information has been thought to depend on studies concerning Trachelium caeruleum and genome announcements for Adenophora remotiflora, Campanula takesimana, and Hanabusaya asiatica. RNA editing information for plastid genes is currently unavailable for Campanulaceae. To understand plastid genome evolution in Campanulaceae, we have sequenced and characterized the chloroplast (cp) genome and nr-accD of Platycodon grandiflorum, a basal member of Campanulaceae.

RESULTS

We sequenced the 171,818 bp cp genome containing a 79,061 bp large single-copy (LSC) region, a 42,433 bp inverted repeat (IR) and a 7840 bp small single-copy (SSC) region, which represents the cp genome with the largest IR among species of Campanulaceae. The genome contains 110 genes and 18 introns, comprising 77 protein-coding genes, four RNA genes, 29 tRNA genes, 17 group II introns, and one group I intron. RNA editing of genes was detected in 18 sites of 14 protein-coding genes. Platycodon has an IR containing a 3' rps12 operon, which occurs in the middle of the LSC region in four other species of Campanulaceae (T. caeruleum, A. remotiflora, C. takesimana, and H. asiatica), but lacks accD, clpP, infA, and rpl23, as has been found in these four species. Platycodon nr-accD contains about 3.2 kb intron between nr-accD.e1 and nr-accD.e2 at the same insertion point as in other Campanulaceae. The phylogenies of the plastid genomes and accD show that Platycodon is basal in the Campanulaceae clade, indicating that IR disruption in Campanulaceae occurred after the loss of accD, clpP, infA, and rpl23 in the cp genome, which occurred during plastid evolution in Campanulaceae.

CONCLUSIONS

The plastid genome of P. grandiflorum lacks the rearrangement of the IR found in T. caeruleum, A. remotiflora, C. takesimana, and H. asiatica. The absence of accD, clpP, infA, and rpl23 in the plastid genome is a synapomorphic characteristic of Campanulaceae. The chloroplast genome phylogeny supports the hypothesis that chloroplast genomic arrangement occurred after accD nuclear transfer and loss of the four genes in the plastid of early Campanulaceae as a lineage of asterids.

Recombination-dependent replication and gene conversion homogenize repeat sequences and diversify plastid genome structure

PREMISE OF THE STUDY

There is a misinterpretation in the literature regarding the variable orientation of the small single copy region of plastid genomes (plastomes). The common phenomenon of small and large single copy inversion, hypothesized to occur through intramolecular recombination between inverted repeats (IR) in a circular, single unit-genome, in fact, more likely occurs through recombination-dependent replication (RDR) of linear plastome templates. If RDR can be primed through both intra- and intermolecular recombination, then this mechanism could not only create inversion isomers of so-called single copy regions, but also an array of alternative sequence arrangements.

METHODS

We used Illumina paired-end and PacBio single-molecule real-time (SMRT) sequences to characterize repeat structure in the plastome of Monsonia emarginata (Geraniaceae). We used OrgConv and inspected nucleotide alignments to infer ancestral nucleotides and identify gene conversion among repeats and mapped long (>1 kb) SMRT reads against the unit-genome assembly to identify alternative sequence arrangements.

RESULTS

Although M. emarginata lacks the canonical IR, we found that large repeats (>1 kilobase; kb) represent ∼22% of the plastome nucleotide content. Among the largest repeats (>2 kb), we identified GC-biased gene conversion and mapping filtered, long SMRT reads to the M. emarginata unit-genome assembly revealed alternative, substoichiometric sequence arrangements.

CONCLUSION

We offer a model based on RDR and gene conversion between long repeated sequences in the M. emarginata plastome and provide support that both intra-and intermolecular recombination between large repeats, particularly in repeat-rich plastomes, varies unit-genome structure while homogenizing the nucleotide sequence of repeats.

Chloroplasts: state of research and practical applications of plastome sequencing

This review presents origins, structure and expression of chloroplast genomes. It also describes their sequencing, analysis and modification, focusing on potential practical uses and biggest challenges of chloroplast genome modification. During the evolution of eukaryotes, cyanobacteria are believed to have merged with host heterotrophic cell. Afterward, most of cyanobacterial genes from cyanobacteria were transferred to cell nucleus or lost in the process of endosymbiosis. As a result of these changes, a primary plastid was established. Nowadays, plastid genome (plastome) is almost always circular, has a size of 100-200 kbp (120-160 in land plants), and harbors 100-120 highly conserved unique genes. Plastids have their own gene expression system, which is similar to one of their cyanobacterial ancestors. Two different polymerases, plastid-derived PEP and nucleus-derived NEP, participate in transcription. Translation is similar to the one observed in cyanobacteria, but it also utilizes protein translation factors and positive regulatory mRNA elements absent from bacteria. Plastoms play an important role in genetic transformation. Transgenes are introduced into them either via gene gun (in undamaged tissues) or polyethylene glycol treatment (when protoplasts are targeted). Antibiotic resistance markers are the most common tool used for selection of transformed plants. In recent years, plastome transformation emerged as a promising alternative to nuclear transformation because of (1) high yield of target protein, (2) removing the risk of outcrossing with weeds, (3) lack of silencing mechanisms, and (4) ability to engineer the entire metabolic pathways rather than single gene traits. Currently, the main directions of such research regard: developing efficient enzyme, vaccine antigen, and biopharmaceutical protein production methods in plant cells and improving crops by increasing their resistance to a wide array of biotic and abiotic stresses. Because of that, the detailed knowledge of plastome structure and mechanism of functioning started to play a major role.

Chloroplast DNA Copy Number Changes during Plant Development in Organelle DNA Polymerase Mutants

Chloroplast genome copy number is very high in leaf tissue, with upwards of 10,000 or more copies of the chloroplast DNA (ctDNA) per leaf cell. This is often promoted as a major advantage for engineering the plastid genome, as it provides high gene copy number and thus is expected to result in high expression of foreign proteins from integrated genes. However, it is also known that ctDNA copy number and ctDNA integrity decrease as cells age. Quantitative PCR (qPCR) allows measurement of organelle DNA levels relative to a nuclear gene target. We have used this approach to determine changes in copy number of ctDNA relative to the nuclear genome at different ages of Arabidopsis plant growth and in organellar DNA polymerase mutants. The mutant plant lines have T-DNA insertions in genes encoding the two organelle localized DNA polymerases (PolIA and PolIB). Each of these mutant lines exhibits some delay in plant growth and development as compared to wild-type plants, with the PolIB plants having a more pronounced delay. Both mutant lines develop to maturity and produce viable seeds. Mutants for both proteins were observed to have a reduction in ctDNA and mtDNA copy number relative to wild type plants at all time points as measured by qPCR. Both DNA polymerase mutants had a fairly similar decrease in ctDNA copy number, while the PolIB mutant had a greater effect of reduction in mtDNA levels. However, despite similar decreases in genome copy number, RT-PCR analysis of PolIA mutants show that PolIB expression remains unchanged, suggesting that PolIA may not be essential to plant survival. Furthermore, genotypic analysis of plants from heterozygous parents display a strong pressure to maintain two functioning copies of PolIB. These results indicate that the two DNA polymerases are both important in ctDNA replication, and they are not fully redundant to each other, suggesting each has a specific function in plant organelles.

The linear plastid chromosomes of maize: terminal sequences, structures, and implications for DNA replication

The structure of a chromosomal DNA molecule may influence the way in which it is replicated and inherited. For decades plastid DNA (ptDNA) was believed to be circular, with breakage invoked to explain linear forms found upon extraction from the cell. Recent evidence indicates that ptDNA in vivo consists of linear molecules with discrete termini, although these ends were not characterized. We report the sequences of two terminal regions, End1 and End2, for maize (Zea mays L.) ptDNA. We describe structural features of these terminal regions and similarities found in other plant ptDNAs. The terminal sequences are within inverted repeat regions (leading to four genomic isomers) and adjacent to origins of replication. Conceptually, stem-loop structures may be formed following melting of the double-stranded DNA ends. Exonuclease digestion indicates that the ends in maize are unobstructed, but tobacco (Nicotiana tabacum L.) ends may have a 5'-protein. If the terminal structure of ptDNA molecules influences the retention of ptDNA, the unprotected molecular ends in mature leaves of maize may be more susceptible to degradation in vivo than the protected ends in tobacco. The terminal sequences and cumulative GC skew profiles are nearly identical for maize, wheat (Triticum aestivum L.) and rice (Oryza sativa L.), with less similarity among other plants. The linear structure is now confirmed for maize ptDNA and inferred for other plants and suggests a virus-like recombination-dependent replication mechanism for ptDNA. Plastid transformation vectors containing the terminal sequences may increase the chances of success in generating transplastomic cereals.

DNA maintenance in plastids and mitochondria of plants

The DNA molecules in plastids and mitochondria of plants have been studied for over 40 years. Here, we review the data on the circular or linear form, replication, repair, and persistence of the organellar DNA (orgDNA) in plants. The bacterial origin of orgDNA appears to have profoundly influenced ideas about the properties of chromosomal DNA molecules in these organelles to the point of dismissing data inconsistent with ideas from the 1970s. When found at all, circular genome-sized molecules comprise a few percent of orgDNA. In cells active in orgDNA replication, most orgDNA is found as linear and branched-linear forms larger than the size of the genome, likely a consequence of a virus-like DNA replication mechanism. In contrast to the stable chromosomal DNA molecules in bacteria and the plant nucleus, the molecular integrity of orgDNA declines during leaf development at a rate that varies among plant species. This decline is attributed to degradation of damaged-but-not-repaired molecules, with a proposed repair cost-saving benefit most evident in grasses. All orgDNA maintenance activities are proposed to occur on the nucleoid tethered to organellar membranes by developmentally-regulated proteins.

Short-range inversions: rethinking organelle genome stability: template switching events during DNA replication destabilize organelle genomes

In the organelles of plants and mammals, recent evidence suggests that genomic instability stems in large part from template switching events taking place during DNA replication. Although more than one mechanism may be responsible for this, some similarities exist between the different proposed models. These can be separated into two main categories, depending on whether they involve a single-strand-switching or a reciprocal-strand-switching event. Single-strand-switching events lead to intermediates containing Y junctions, whereas reciprocal-strand-switching creates Holliday junctions. Common features in all the described models include replication stress, fork stalling and the presence of inverted repeats, but no single element appears to be required in all cases. We review the field, and examine the ideas that several mechanisms may take place in any given genome, and that the presence of palindromes or inverted repeats in certain regions may favor specific rearrangements.

Changes in DNA damage, molecular integrity, and copy number for plastid DNA and mitochondrial DNA during maize development

The amount and structural integrity of organellar DNAs change during plant development, although the mechanisms of change are poorly understood. Using PCR-based methods, we quantified DNA damage, molecular integrity, and genome copy number for plastid and mitochondrial DNAs of maize seedlings. A DNA repair assay was also used to assess DNA impediments. During development, DNA damage increased and molecules with impediments that prevented amplification by Taq DNA polymerase increased, with light causing the greatest change. DNA copy number values depended on the assay method, with standard real-time quantitative PCR (qPCR) values exceeding those determined by long-PCR by 100- to 1000-fold. As the organelles develop, their DNAs may be damaged in oxidative environments created by photo-oxidative reactions and photosynthetic/respiratory electron transfer. Some molecules may be repaired, while molecules with unrepaired damage may be degraded to non-functional fragments measured by standard qPCR but not by long-PCR.

Nucleases in higher plants and their possible involvement in DNA degradation during leaf senescence

During leaf senescence, macromolecules such as proteins and lipids are known to be degraded for redistribution into upper tissues. Similarly, nucleic acids appear to undergo fragmentation or degradation during senescence, but the physiological role of nucleic acid degradation, particularly of genomic DNA degradation, remains unclear. To date, more than a dozen of plant deoxyribonucleases have been reported, whereas it remains to be verified whether any of them degrade DNA during leaf senescence. This review summarizes current knowledge related to the plant nucleases that are induced developmentally or in a tissue-specific manner and are known to degrade DNA biochemically. Of these, several endonucleases (BFN1, CAN1, and CAN2) and an exonuclease (DPD1) in Arabidopsis seem to act in leaf senescence because they were shown to be inducible at the transcript level. This review specifically examines DPD1, which is dual-targeted to chloroplasts and mitochondria. Results show that, among the exonuclease family to which DPD1 belongs, DPD1 expression is extraordinary when estimated using a microarray database. DPD1 is the only example among the nucleases in which DNA degradation has been confirmed in vivo in pollen by mutant analysis. These data imply a significant role of organelle DNA degradation during leaf senescence and implicate DPD1 as a potential target for deciphering nucleotide salvage in plants.

Replication intermediates of the linear mitochondrial DNA of Candida parapsilosis suggest a common recombination based mechanism for yeast mitochondria

Variation in the topology of mitochondrial DNA (mtDNA) in eukaryotes evokes the question if differently structured DNAs are replicated by a common mechanism. RNA-primed DNA synthesis has been established as a mechanism for replicating the circular animal/mammalian mtDNA. In yeasts, circular mtDNA molecules were assumed to be templates for rolling circle DNA-replication. We recently showed that in Candida albicans, which has circular mapping mtDNA, recombination driven replication is a major mechanism for replicating a complex branched mtDNA network. Careful analyses of C. albicans-mtDNA did not reveal detectable amounts of circular DNA molecules. In the present study we addressed the question of how the unit sized linear mtDNA of Candida parapsilosis terminating at both ends with arrays of tandem repeats (mitochondrial telomeres) is replicated. Originally, we expected to find replication intermediates diagnostic of canonical bi-directional replication initiation at the centrally located bi-directional promoter region. However, we found that the linear mtDNA of Candida parapsilosis also employs recombination for replication initiation. The most striking findings were that the mitochondrial telomeres appear to be hot spots for recombination driven replication, and that stable RNA:DNA hybrids, with a potential role in mtDNA replication, are also present in the mtDNA preparations.