Variation in copy number of a 24-base pair tandem repeat in the chloroplast DNA of Oenothera hookeri strain Johansen

A Dynamic Tandem Repeat in Monocotyledons Inferred from a Comparative Analysis of Chloroplast Genomes in Melanthiaceae

Chloroplast genomes (cpDNA) are highly valuable resources for evolutionary studies of angiosperms, since they are highly conserved, are small in size, and play critical roles in plants. Slipped-strand mispairing (SSM) was assumed to be a mechanism for generating repeat units in cpDNA. However, research on the employment of different small repeated sequences through SSM events, which may induce the accumulation of distinct types of repeats within the same region in cpDNA, has not been documented. Here, we sequenced two chloroplast genomes from the endemic species Heloniopsis tubiflora (Korea) and Xerophyllum tenax (USA) to cover the gap between molecular data and explore "hot spots" for genomic events in Melanthiaceae. Comparative analysis of 23 complete cpDNA sequences revealed that there were different stages of deletion in the rps16 region across the Melanthiaceae. Based on the partial or complete loss of rps16 gene in cpDNA, we have firstly reported potential molecular markers for recognizing two sections (Veratrum and Fuscoveratrum) of Veratrum. Melathiaceae exhibits a significant change in the junction between large single copy and inverted repeat regions, ranging from trnH_GUG to a part of rps3. Our results show an accumulation of tandem repeats in the rpl23-ycf2 regions of cpDNAs. Small conserved sequences exist and flank tandem repeats in further observation of this region across most of the examined taxa of Liliales. Therefore, we propose three scenarios in which different small repeated sequences were used during SSM events to generate newly distinct types of repeats. Occasionally, prior to the SSM process, point mutation event and double strand break repair occurred and induced the formation of initial repeat units which are indispensable in the SSM process. SSM may have likely occurred more frequently for short repeats than for long repeat sequences in tribe Parideae (Melanthiaceae, Liliales). Collectively, these findings add new evidence of dynamic results from SSM in chloroplast genomes which can be useful for further evolutionary studies in angiosperms. Additionally, genomics events in cpDNA are potential resources for mining molecular markers in Liliales.

Analysis of variation in chloroplast DNA sequences

This chapter introduces and reviews methods for analyzing variation in chloroplast DNA, mainly by polymerase chain reaction (PCR) and subsequent revelation of polymorphisms. Sources for chloroplast primers are discussed, as well as methods such as Sanger sequencing, PCR followed by restriction fragment length polymorphism (RFLP), gel electrophoresis, fragment analysis on automated DNA sequencers, denaturing high-performance liquid chromatography (dHPLC), and next-generation sequencing (NGS). A special section deals with peculiarities of chloroplast DNA variation, such as tandem repeats and mini- and microsatellites.

A frameshift mutation of the chloroplast matK coding region is associated with chlorophyll deficiency in the Cryptomeria japonica virescent mutant Wogon-Sugi

Wogon-Sugi has been reported as a cytoplasmically inherited virescent mutant selected from a horticultural variety of Cryptomeria japonica. Although previous studies of plastid structure and inheritance indicated that at least some mutations are encoded by the chloroplast genome, the causative gene responsible for the primary chlorophyll deficiency in Wogon-Sugi, has not been identified. In this study, we identified this gene by genomic sequencing of chloroplast DNA and genetic analysis. Chloroplast DNA sequencing of 16 wild-type and 16 Wogon-Sugi plants showed a 19-bp insertional sequence in the matK coding region in the Wogon-Sugi. This insertion disrupted the matK reading frame. Although an indel mutation in the ycf1 and ycf2 coding region was detected in Wogon-Sugi, sequence variations similar to that of Wogon-Sugi were also detected in several wild-type lines, and they maintained the reading frame. Genetic analysis of the 19 bp insertional mutation in the matK coding region showed that it was found only in the chlorophyll-deficient sector of 125 full-sibling seedlings. Therefore, the 19-bp insertion in the matK coding region is the most likely candidate at present for a mutation underlying the Wogon-Sugi phenotype.

Plastidic trnFUUC pseudogenes in North American genus Boechera (Brassicaceae): mechanistic aspects of evolution

The origin and maintenance of a plastidic tandem repeat next to the TRNF (UUC) gene were analyzed in the genus BOECHERA in a phylogenetic context and were compared to published analogous examples that emerged in parallel in the Asteraceae and Juncaceae, respectively. Although we identified some features common to these taxonomic groups with respect to structure and origin of the region, obvious differences were encountered, which argue against a specific mechanism or evolutionary principle underlying the parallel origin and maintenance of the TRNF-tandem repeats in those families. In contrast to the situation in the Asteraceae, no reciprocal recombinant repeat types have been observed in the Brassicaceae. Forty copy types, classified into three groups, were isolated from 103 chloroplast haplotypes of BOECHERA and it was demonstrated that they are composed of four subregions of various origins. We discuss various mutation mechanisms such as DNA replication slippage, and inter- and intrachromosomal recombination which were reported to mediate variation in copy numbers and other types of observed sequence length polymorphism. It is shown that the observed molecular structure of the tandem repeat region did not fully fit the particular patterns expected under a scenario of evolution including any of the known mechanisms. Nevertheless, it appeared that intermolecular unequal crossing-over is most likely the driving force in the evolution of this tandem repeat. However, it remains to be explained, why no reciprocal recombinant copy types have been observed. The reconstructed phylogenetic relationships among copies reflected different evolutionary scenarios as follows: (1) A single and ancient origin of copies pre-dates the radiation of BOECHERA. (2) Parallel expansion and shortening of the tandem repeat within different BOECHERA lineages. (3) Conservation of the first copy, as it was the only one present in all chloroplast haplotypes.

Hypervariable plastid locus variation and intron evolution in the Anacamptis palustris lineage

Data on the organization of a hypervariable chloroplast locus in the Anacamptis palustris (Orchidaceae) lineage are provided and used to infer the pattern of molecular evolution in this group. A large survey of sequence variation in A. palustris and allied taxa reveals several repeat types differing in number and organization that occur in the same plastid region. The resulting repeat type network suggests that at least seven different minisatellite loci evolved near each other in the tRNALEU intron and indicates the presence of at least three main phyletic lines in the A. palustris lineage.

Molecular evolution of a plastid tandem repeat locus in an orchid lineage

The molecular evolution of a chloroplast minisatellite locus in the Anacamptis palustris (Orchidaceae) lineage and haplotype variation in two Italian A. palustris populations were investigated. A phylogenetic analyses of the chloroplast tRNA(LEU) intron, where the minisatellite locus is located, revealed that a deletion in the ancestor of the A. palustris lineage led to the formation of two noncontiguous, complementary sequence motifs. We propose a model to explain the initial formation of the minisatellite repeat motif, starting with the two noncontiguous, complementary sequence motifs. A survey of minisatellite variation in four species of the A. palustris lineage revealed several haplotypes that differed not only in repeat number, but also in repeat organization. A haplotype network suggests that three different minisatellite loci evolved independently at the same position in the tRNA(LEU) intron. A secondary structure model revealed that the A. palustris minisatellite repeat forms a stem region of the tRNA(LEU) intron, which allows its notable expansion without negatively affecting splicing. Minisatellite variation was high in the two examined A. palustris populations where 20 haplotypes were detected, whereas no length variation was detected in a neighboring poly (A) microsatellite locus. We estimated a chloroplast minisatellite mutation rate of 3.2 x 10(-3) mutations per generation. Southern blot analyses did not find evidence for chloroplast heteroplasmy. Based on the analysis of the largest known, extant A. palustris population, a stepwise mutation model (SMM) was inferred.

Chloroplast mutations induced by 9-aminoacridine hydrochloride are independent of the plastome mutator in Oenothera

Oenothera plants homozygous for the recessive plastome mutator allele ( pm) show chloroplast DNA (cpDNA) mutation frequencies that are about 1,000-fold higher than spontaneous levels. The pm-encoded gene product has been hypothesized to have a function in cpDNA replication, repair and/or mutation avoidance. Previous chemical mutagenesis experiments with the alkylating agent nitroso-methyl urea (NMU) showed a synergistic effect of NMU on the induction of mutations in the pm line, suggesting an interaction between the pm-encoded gene product and one of the repair systems that corrects alkylation damage. The goal of the experiments described here was to examine whether the pm activity extends to the repair of damage caused by non-alkylating mutagens. To this end, the intercalating mutagen, 9-aminoacridine hydrochloride (9AA) was tested for synergism with the plastome mutator. A statistical analysis of the data reported here indicates that the pm-encoded gene product is not involved in the repair of the 9AA-induced mutations. However, the recovery of chlorotic sectors in plants derived from the mutagenized seeds shows that 9AA can act as a mutagen of the chloroplast genome.

A variable minisatellite sequence in the chloroplast genome of Sorbus L. (Rosaceae: Maloideae)

The chloroplast genome is now known to be more variable than was once thought. Reports of RFLP (restriction fragment length polymorphism) and sequence variation, as well as variation in chloroplast microsatellites, are common. Here, data are presented on the variability of a minisatellite sequence in the chloroplast genome of Sorbus species. RFLP analysis of a PCR product comprising the region between the trnM and rbcL genes of nine Sorbus species identified seven size variants. Sequencing revealed the observed size polymorphism to be due to differences in the number of copies of an imperfect 9-bp motif. A more intensive survey of the variability of the minisatellite was undertaken in populations of Sorbus aucuparia. The potential uses of such regions in chloroplast DNA are discussed and a possible mechanism for the evolution of the minisatellite is presented.

Plastome mutator-induced alterations arise in Oenothera chloroplast DNA through template slippage

The plastome mutator of Oenothera hookeri strain Johansen causes deletions and duplications at target sites defined by direct repeats in the plastid genome. Previous studies characterized the mutations long after they had occurred and could not discriminate between the possibilities that the plastome mutator acted through unequal homologous recombination or template slippage. From the known hotspots, the rRNA spacer in the large inverted repeat was chosen for this study because it contains both direct and indirect repeats. Identical deletions were recovered from independently derived plants; the altered regions were always flanked by direct repeats. The regions in which the deletions occurred have the potential to form secondary structures that would stabilize the intervening sequence. Of the two affected regions, the one with the stronger potential secondary structure was altered more frequently. Because no duplication products or inversions were recovered, it is proposed that the plastome mutator acts through template slippage rather than through a recombination mechanism.

Coordination of nuclear and chloroplast gene expression in plant cells

Plastid proteins are encoded in two genomes, one in the nucleus and the other in the organelle. The expression of genes in these two compartments in coordinated during development and in response to environmental parameters such as light. Two converging approaches reveal features of this coordination: the biochemical analysis of proteins involved in gene expression, and the genetic analysis of mutants affected in plastid function or development. Because the majority of proteins implicated in plastid gene expression are encoded in the nucleus, regulatory processes in the nucleus and in the cytoplasm control plastid gene expression, in particular during development. Many nucleus-encoded factors involved in transcriptional and posttranscriptional steps of plastid gene expression have been characterized. We are also beginning to understand whether and how certain developmental or environmental signals perceived in one compartment may be transduced to the other.

A mutation hotspot in the chloroplast genome of a conifer (Douglas-fir: Pseudotsuga) is caused by variability in the number of direct repeats derived from a partially duplicated tRNA gene

We determined the DNA sequence of a 2.7-kb cpDNA XbaI fragment from douglas-fir [Pseudotsuga menziesii (Mirb.) Franco]. RFLPs revealed by the 2.7-kb XbaI clone were observed to vary up to 1 kb among species within the genus Pseudotsuga and up to 200 bp among trees of P. menziesii. The polymerase chain reaction (PCR) allowed the locus of polymorphism to be identified, and the variable region was then sequenced in a second Douglas-fir tree, a single tree of a related species, Japanese Douglas-fir (P. japonica), and in a species lacking a mutation hotspot in the region, Pinus radiata (Monterey pine). The locus of polymorphism is characterized by hundreds of base pairs of imperfect, tandem direct repeats flanked by a partially duplicated and an intact trn Y-GUA gene. The duplication is direct in orientation and consists of 43 bp of the 3' end of trnY and 25 bp of its 3' flanking sequence. Tandem repeats show high sequence similarity to a 27-bp region of the trnY gene that overlaps one end of the duplication. The two trees of Douglas-fir sequenced differed by a single tandem repeat unit, whereas these trees differed from the Japanese Douglas-fir sequenced by approximately 34 repeat units. Repetitive DNA in the Pseudotsuga cpDNA hotspot was most likely generated at the time of the partial trnY gene duplication and these sequences expanded by slipped-strand mispairing and unequal crossing-over.

Structure and evolution of the largest chloroplast gene (ORF2280): internal plasticity and multiple gene loss during angiosperm evolution

We have determined the nucleotide sequence of the Pelargonium x hortorum ORF2280 homolog, the largest gene in the plastid genome of most land plants, and compared it to published homologs from Nicotiana tabacum, Epifagus virginiana, Spinacia oleracea, and Marchantia polymorpha. Multiple alignment of protein sequences requires an extraordinary number of gaps, indicating a very high frequency of insertion/deletion events during the evolution of the protein; however, the overall predicted size of the protein varies relatively little among the five species. At 2,109 codons, the Pelargonium gene is smaller than other land plant ORF2280 homologs and exhibits a rate of nucleotide substitution several times higher relative to Nicotiana, Epifagus, and Spinacia. Southern-blot and restriction-mapping studies were carried out to uncover length variation in ORF2280 homologs from 279 species (representing 111 families) of angiosperms. In many independent angiosperm lineages, this gene has sustained deletions ranging in size from 200 bp to almost 6 kb. Based on the severity of deletions, we postulate that the chloroplast homolog of ORF2280 has become nonfunctional in at least four independent lineages of angiosperms.

Similarity between putative ATP-binding sites in land plant plastid ORF2280 proteins and the FtsH/CDC48 family of ATPases

Plastid ORF2280 proteins from five species of land plant are shown to have limited amino-acid sequence similarity to a family of proteins that includes the yeast CDC48, SEC18, PAS1 and SUG1 proteins, three subunits of the mammalian 26S protease, and the Escherichia coli FtsH protein. These proteins all contain one or two ATPase domains and many are involved in cell division, transport of proteins across membranes, or proteolysis. Similarity with the ORF2280 proteins is restricted to a single region of about 130 amino acids that contains: (1) sequences resembling a nucleotide binding site but lacking two normally conserved residues, and (2) a downstream conserved motif with the consensus sequence VIX2TX2PX3DPALX2P. Most of the rest of ORF2280 is very poorly conserved among land plants, even though other family members such as CDC48 have slow rates of protein sequence evolution. In contrast, a protein encoded by plastid DNA of the rhodophyte alga Porphyra purpurea is very similar to E. coli FtsH. Phylogenetic analysis suggests that the red and green plastid genes are not true homologues (orthologues) but distinct members of an ancient gene family.

In-frame length mutations associated with short tandem repeats are located in unassigned open reading frames of Oenothera chloroplast DNA

Chloroplast DNAs were compared between two closely related species in the subsection Munzia of the genus Oenothera. A restriction fragment length dimorphism (273 bp) within the large inverted repeats was localized to an unassigned open reading frame that is homologous to ORF 2280 of tobacco chloroplast DNA. This dimorphism is due to different copy numbers of various short tandem repeated sequences, with each repeat unit specifying an in-frame addition or deletion. Other small length mutations were detected within an unassigned reading frame that appears to be homologous to the tobacco ORF 1244, and in the non-coding sequence upstream of that frame. These insertions and/or deletions are all associated with short direct repeats that lie in tandem.

Large unidentified open reading frame in plastid DNA (ORF2280) is expressed in chloroplasts

The chloroplast DNA encodes genes for components of photosynthesis and the transcription-translation machinery; a number of unidentified open reading frames (ORFs) are also present. To determine whether a large ORF in the inverted repeat of chloroplast DNA of tobacco (ORF2280) encodes a chloroplast protein, a conserved region of the ORF was expressed in Escherichia coli. An antibody against the ORF protein was prepared using the purified fusion protein as an antigen. When incubated with proteins from the soluble fraction of tobacco, spinach and Oenothera chloroplasts, the antiserum detects relatively labile polypeptides, which have apparent molecular weights of 170 to 180 kDa. The ORF in tobacco and spinach is large enough to encode a protein of 240-250 kDa, thus it is possible that post-transcriptional or post-translational processing reduces the size of the expression product. Analysis of Oenothera chloroplasts representing four different plastome types revealed endonuclease restriction fragment length polymorphisms in chloroplast DNA indicative of insertion/deletion events in a region of the chloroplast DNA that shared significant sequence similarity with ORF2280. The ORF2280 antiserum was used to demonstrate that there are qualitative differences in the ORF proteins from different Oenothera plastome types.

Structure and expression of cytochrome f in an Oenothera plastome mutant

The chloroplast mutant pm7 is one of a number of mutants derived from the plastome mutator (pm) line of Oenothera hookeri, strain Johansen. Immunoblotting showed that this mutant accumulates a protein that is cross-antigenic with cytochrome f, but five kilodaltons larger than the mature wild-type protein. Since cytochrome f is known to be translated on plastid ribosomes as a precursor with an amino-terminal extension, it is proposed that the unprocessed cytochrome f precursor accumulates in pm7. In addition to this precursor-sized cytochrome f protein, some mature-sized cytochrome f was also found in the mutant plastids. The pm7 mutation is inherited in a non-Mendelian fashion; but no alterations in chloroplast DNA restriction patterns, or differences in DNA sequence in the region encoding cytochrome f, were found in a comparison of the wild-type and pm7 chloroplast DNAs. Although the mutant was capable of synthesizing heme, no covalently-bound heme, normally found associated with mature, functional, cytochrome f was detected in the mutant at sizes expected for the presumed precursor, or for mature cytochrome f. These results indicate that the aberrant accumulation of a precursor-sized cytochrome f in pm7 is not due to a lesion directly in the plastid gene encoding cytochrome f, petA, or to a deficiency in the ability of the mutant plastids to synthesize or accumulate heme.