A transcription map of the chloroplast genome from rice (Oryza sativa)

Characterization of the DNA accessibility of chloroplast genomes in grasses

Although the chloroplast genome (cpDNA) of higher plants is known to exist as a large protein-DNA complex called 'plastid nucleoid', researches on its DNA state and regulatory elements are limited. In this study, we performed the assay for transposase-accessible chromatin sequencing (ATAC-seq) on five common tissues across five grasses, and found that the accessibility of different regions in cpDNA varied widely, with the transcribed regions being highly accessible and accessibility patterns around gene start and end sites varying depending on the level of gene expression. Further analysis identified a total of 3970 putative protein binding footprints on cpDNAs of five grasses. These footprints were enriched in intergenic regions and co-localized with known functional elements. Footprints and their flanking accessibility varied dynamically among tissues. Cross-species analysis showed that footprints in coding regions tended to overlap non-degenerate sites and contain a high proportion of highly conserved sites, indicating that they are subject to evolutionary constraints. Taken together, our results suggest that the accessibility of cpDNA has biological implications and provide new insights into the transcriptional regulation of chloroplasts.

Impacts of phosphatidylglycerol on plastid gene expression and light induction of nuclear photosynthetic genes

Phosphatidylglycerol (PG) is the only major phospholipid in the thylakoid membrane of chloroplasts. PG is essential for photosynthesis, and loss of PG in Arabidopsis thaliana results in severe defects of growth and chloroplast development, with decreased chlorophyll accumulation, impaired thylakoid formation, and down-regulation of photosynthesis-associated genes encoded in nuclear and plastid genomes. However, how the absence of PG affects gene expression and plant growth remains unclear. To elucidate this mechanism, we investigated transcriptional profiles of a PG-deficient Arabidopsis mutant pgp1-2 under various light conditions. Microarray analysis demonstrated that reactive oxygen species (ROS)-responsive genes were up-regulated in pgp1-2. However, ROS production was not enhanced in the mutant even under strong light, indicating limited impacts of photooxidative stress on the defects of pgp1-2. Illumination to dark-adapted pgp1-2 triggered down-regulation of photosynthesis-associated nuclear-encoded genes (PhANGs), while plastid-encoded genes were constantly suppressed. Overexpression of GOLDEN2-LIKE1 (GLK1), a transcription factor gene regulating chloroplast development, in pgp1-2 up-regulated PhANGs but not plastid-encoded genes along with chlorophyll accumulation. Our data suggest a broad impact of PG biosynthesis on nuclear-encoded genes partially via GLK1 and a specific involvement of this lipid in plastid gene expression and plant development.

A rice mTERF protein V14 sustains photosynthesis establishment and temperature acclimation in early seedling leaves

BACKGROUND

Plant mitochondrial transcription termination factor (mTERF) family members play important roles in development and stress tolerance through regulation of organellar gene expression. However, their molecular functions have yet to be clearly defined.

RESULTS

Here an mTERF gene V14 was identified by fine mapping using a conditional albino mutant v14 that displayed albinism only in the first two true leaves, which was confirmed by transgenic complementation tests. Subcellular localization and real-time PCR analyses indicated that V14 encodes a chloroplastic protein ubiquitously expressed in leaves while spiking in the second true leaf. Chloroplastic gene expression profiling in the pale leaves of v14 through real-time PCR and Northern blotting analyses showed abnormal accumulation of the unprocessed transcripts covering the rpoB-rpoC1 and/or rpoC1-rpoC2 intercistronic regions accompanied by reduced abundance of the mature rpoC1 and rpoC2 transcripts, which encode two core subunits of the plastid-encoded plastid RNA polymerase (PEP). Subsequent immunoblotting analyses confirmed the reduced accumulation of RpoC1 and RpoC2. A light-inducible photosynthetic gene psbD was also found down-regulated at both the mRNA and protein levels. Interestingly, such stage-specific aberrant posttranscriptional regulation and psbD expression can be reversed by high temperatures (30 ~ 35 °C), although V14 expression lacks thermo-sensitivity. Meanwhile, three V14 homologous genes were found heat-inducible with similar temporal expression patterns, implicating their possible functional redundancy to V14.

CONCLUSIONS

These data revealed a critical role of V14 in chloroplast development, which impacts, in a stage-specific and thermo-sensitive way, the appropriate processing of rpoB-rpoC1-rpoC2 precursors and the expression of certain photosynthetic proteins. Our findings thus expand the knowledge of the molecular functions of rice mTERFs and suggest the contributions of plant mTERFs to photosynthesis establishment and temperature acclimation.

Complete chloroplast genome sequencing of Job's tears (Coix L.): genome structure, comparative analysis, and phylogenetic relationships

Job's tears, also known as adlay, is a valuable plant that has commonly been used in traditional Chinese medicine, as well as an edible food. Due to the lack of knowledge of its genetics and gaps in its evolutionary analysis, breeding of adlay has been hindered. Here, we report five complete chloroplast genomes of various species and varieties in the genus by Illumina sequencing, while their genome structure, comparative analysis, and phylogenetic relationships were conducted. Genome sizes ranged from 140,860 to 140,864 bp in length, GC contents were 38.43%, and genome architecture was of a typical quadripartite structure. We annotated 82～83 protein-coding genes and 46～47 non-coding RNA genes in each genome and they functionally associated with self-replication, photosynthesis, cytochrome synthesis and other unknown functions. Three codons that encoded tryptophan, arginine and leucine were used frequently at rates of 41.42, 37.98, and 32.28% respectively. The preferred codons consistently ended with A or T. A total of 146 simple sequence repeats (SSR), 9 insertions and deletions (InDels) and 143 single nucleotide polymorphisms (SNPs) were observed among genomes. The InDel and SNP variations were mostly distributed in intergenic regions. It confirmed that Coix, Sorghum, Saccharum, Zea, Tripsacum and Saccharum were closely genera and the genetic distance of Sorghum to Coix was closer than Zea to Coix. These results give us more insight into the evolution of Coix in a wide range of evolutionary studies.

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.

Complete Chloroplast Genome Sequence and Phylogenetic Analysis of Aster tataricus

We sequenced and analyzed the complete chloroplast genome of Aster tataricus (family Asteraceae), a Chinese herb used medicinally to relieve coughs and reduce sputum. The A. tataricus chloroplast genome was 152,992 bp in size, and harbored a pair of inverted repeat regions (IRa and IRb, each 24,850 bp) divided into a large single-copy (LSC, 84,698 bp) and a small single-copy (SSC, 18,250 bp) region. Our annotation revealed that the A. tataricus chloroplast genome contained 115 genes, including 81 protein-coding genes, 4 ribosomal RNA genes, and 30 transfer RNA genes. In addition, 70 simple sequence repeats (SSRs) were detected in the A. tataricus chloroplast genome, including mononucleotides (36), dinucleotides (1), trinucleotides (23), tetranucleotides (1), pentanucleotides (8), and hexanucleotides (1). Comparative chloroplast genome analysis of three Aster species indicated that a higher similarity was preserved in the IR regions than in the LSC and SSC regions, and that the differences in the degree of preservation were slighter between A. tataricus and A. altaicus than between A. tataricus and A. spathulifolius. Phylogenetic analysis revealed that A. tataricus was more closely related to A. altaicus than to A. spathulifolius. Our findings offer valuable information for future research on Aster species identification and selective breeding.

Complete chloroplast genome sequence and comparative analysis of loblolly pine (Pinus taeda L.) with related species

Pinaceae, the largest family of conifers, has a diversified organization of chloroplast (cp) genomes with two typical highly reduced inverted repeats (IRs). In the current study, we determined the complete sequence of the cp genome of an economically and ecologically important conifer tree, the loblolly pine (Pinus taeda L.), using Illumina paired-end sequencing and compared the sequence with those of other pine species. The results revealed a genome size of 121,531 base pairs (bp) containing a pair of 830-bp IR regions, distinguished by a small single copy (42,258 bp) and large single copy (77,614 bp) region. The chloroplast genome of P. taeda encodes 120 genes, comprising 81 protein-coding genes, four ribosomal RNA genes, and 35 tRNA genes, with 151 randomly distributed microsatellites. Approximately 6 palindromic, 34 forward, and 22 tandem repeats were found in the P. taeda cp genome. Whole cp genome comparison with those of other Pinus species exhibited an overall high degree of sequence similarity, with some divergence in intergenic spacers. Higher and lower numbers of indels and single-nucleotide polymorphism substitutions were observed relative to P. contorta and P. monophylla, respectively. Phylogenomic analyses based on the complete genome sequence revealed that 60 shared genes generated trees with the same topologies, and P. taeda was closely related to P. contorta in the subgenus Pinus. Thus, the complete P. taeda genome provided valuable resources for population and evolutionary studies of gymnosperms and can be used to identify related species.

Discovery of novel plastid phenylalanine (trnF) pseudogenes defines a distinctive clade in Solanaceae

BACKGROUND

The plastome of embryophytes is known for its high degree of conservation in size, structure, gene content and linear order of genes. The duplication of entire tRNA genes or their arrangement in a tandem array composed by multiple pseudogene copies is extremely rare in the plastome. Pseudogene repeats of the trnF gene have rarely been described from the chloroplast genome of angiosperms.

FINDINGS

We report the discovery of duplicated copies of the original phenylalanine (trnFGAA) gene in Solanaceae that are specific to a larger clade within the Solanoideae subfamily. The pseudogene copies are composed of several highly structured motifs that are partial residues or entire parts of the anticodon, T- and D-domains of the original trnF gene.

CONCLUSIONS

The Pseudosolanoid clade consists of 29 genera and includes many economically important plants such as potato, tomato, eggplant and pepper.

The Complete Plastid Genome Sequence of Madagascar Periwinkle Catharanthus roseus (L.) G. Don: Plastid Genome Evolution, Molecular Marker Identification, and Phylogenetic Implications in Asterids

The Madagascar periwinkle (Catharanthusroseus in the family Apocynaceae) is an important medicinal plant and is the source of several widely marketed chemotherapeutic drugs. It is also commonly grown for its ornamental values and, due to ease of infection and distinctiveness of symptoms, is often used as the host for studies on phytoplasmas, an important group of uncultivated plant pathogens. To gain insights into the characteristics of apocynaceous plastid genomes (plastomes), we used a reference-assisted approach to assemble the complete plastome of C. roseus, which could be applied to other C. roseus-related studies. The C. roseus plastome is the second completely sequenced plastome in the asterid order Gentianales. We performed comparative analyses with two other representative sequences in the same order, including the complete plastome of Coffeaarabica (from the basal Gentianales family Rubiaceae) and the nearly complete plastome of Asclepiassyriaca (Apocynaceae). The results demonstrated considerable variations in gene content and plastome organization within Apocynaceae, including the presence/absence of three essential genes (i.e., accD, clpP, and ycf1) and large size changes in non-coding regions (e.g., rps2-rpoC2 and IRb-ndhF). To find plastome markers of potential utility for Catharanthus breeding and phylogenetic analyses, we identified 41 C. roseus-specific simple sequence repeats. Furthermore, five intergenic regions with high divergence between C. roseus and three other euasterids I taxa were identified as candidate markers. To resolve the euasterids I interordinal relationships, 82 plastome genes were used for phylogenetic inference. With the addition of representatives from Apocynaceae and sampling of most other asterid orders, a sister relationship between Gentianales and Solanales is supported.

Complete plastid genome sequence of the basal asterid Ardisia polysticta Miq. and comparative analyses of asterid plastid genomes

Ardisia is a basal asterid genus well known for its medicinal values and has the potential for development of novel phytopharmaceuticals. In this genus of nearly 500 species, many ornamental species are commonly grown worldwide and some have become invasive species that caused ecological problems. As there is no completed plastid genome (plastome) sequence in related taxa, we sequenced and characterized the plastome of Ardisia polysticta to find plastid markers of potential utility for phylogenetic analyses at low taxonomic levels. The complete A. polysticta plastome is 156,506 bp in length and has gene content and organization typical of most asterids and other angiosperms. We identified seven intergenic regions as potentially informative markers with resolution for interspecific relationships. Additionally, we characterized the diversity of asterid plastomes with respect to GC content, plastome organization, gene content, and repetitive sequences through comparative analyses. The results demonstrated that the genome organizations near the boundaries between inverted repeats (IRs) and single-copy regions (SCs) are polymorphic. The boundary organization found in Ardisia appears to be the most common type among asterids, while six other types are also found in various asterid lineages. In general, the repetitive sequences in genic regions tend to be more conserved, whereas those in noncoding regions are usually lineage-specific. Finally, we inferred the whole-plastome phylogeny with the available asterid sequences. With the improvement in taxon sampling of asterid orders and families, our result highlights the uncertainty of the position of Gentianales within euasterids I.

Contradiction between plastid gene transcription and function due to complex posttranscriptional splicing: an exemplary study of ycf15 function and evolution in angiosperms

Plant chloroplast genes are usually co-transcribed while its posttranscriptional splicing is fairly complex and remains largely unsolved. On basis of sequencing the three complete Camellia (Theaceae) chloroplast genomes for the first time, we comprehensively analyzed the evolutionary patterns of ycf15, a plastid gene quite paradoxical in terms of its function and evolution, along the inferred angiosperm phylogeny. Although many species in separate lineages including the three species reported here contained an intact ycf15 gene in their chloroplast genomes, the phylogenetic mixture of both intact and obviously disabled ycf15 genes imply that they are all non-functional. Both intracellular gene transfer (IGT) and horizontal gene transfer (HGT) failed to explain such distributional anomalies. While, transcriptome analyses revealed that ycf15 was transcribed as precursor polycistronic transcript which contained ycf2, ycf15 and antisense trnL-CAA. The transcriptome assembly was surprisingly found to cover near the complete Camellia chloroplast genome. Many non-coding regions including pseudogenes were mapped by multiple transcripts, indicating the generality of pseudogene transcriptions. Our results suggest that plastid DNA posttranscriptional splicing may involve complex cleavage of non-functional genes.

The primary transcriptome of barley chloroplasts: numerous noncoding RNAs and the dominating role of the plastid-encoded RNA polymerase

Gene expression in plastids of higher plants is dependent on two different transcription machineries, a plastid-encoded bacterial-type RNA polymerase (PEP) and a nuclear-encoded phage-type RNA polymerase (NEP), which recognize distinct types of promoters. The division of labor between PEP and NEP during plastid development and in mature chloroplasts is unclear due to a lack of comprehensive information on promoter usage. Here, we present a thorough investigation into the distribution of PEP and NEP promoters within the plastid genome of barley (Hordeum vulgare). Using a novel differential RNA sequencing approach, which discriminates between primary and processed transcripts, we obtained a genome-wide map of transcription start sites in plastids of mature first leaves. PEP-lacking plastids of the albostrians mutant allowed for the unambiguous identification of NEP promoters. We observed that the chloroplast genome contains many more promoters than genes. According to our data, most genes (including genes coding for photosynthesis proteins) have both PEP and NEP promoters. We also detected numerous transcription start sites within operons, indicating transcriptional uncoupling of genes in polycistronic gene clusters. Moreover, we mapped many transcription start sites in intergenic regions and opposite to annotated genes, demonstrating the existence of numerous noncoding RNA candidates.

The mitochondrial genomes of the early land plants Treubia lacunosa and Anomodon rugelii: dynamic and conservative evolution

Early land plant mitochondrial genomes captured important changes of mitochondrial genome evolution when plants colonized land. The chondromes of seed plants show several derived characteristics, e.g., large genome size variation, rapid intra-genomic rearrangement, abundant introns, and highly variable levels of RNA editing. On the other hand, the chondromes of charophytic algae are still largely ancestral in these aspects, resembling those of early eukaryotes. When the transition happened has been a long-standing question in studies of mitochondrial genome evolution. Here we report complete mitochondrial genome sequences from an early-diverging liverwort, Treubia lacunosa, and a late-evolving moss, Anomodon rugelii. The two genomes, 151,983 and 104,239 base pairs in size respectively, contain standard sets of protein coding genes for respiration and protein synthesis, as well as nearly full sets of rRNA and tRNA genes found in the chondromes of the liverworts Marchantia polymorpha and Pleurozia purpurea and the moss Physcomitrella patens. The gene orders of these two chondromes are identical to those of the other liverworts and moss. Their intron contents, with all cis-spliced group I or group II introns, are also similar to those in the previously sequenced liverwort and moss chondromes. These five chondromes plus the two from the hornworts Phaeoceros laevis and Megaceros aenigmaticus for the first time allowed comprehensive comparative analyses of structure and organization of mitochondrial genomes both within and across the three major lineages of bryophytes. These analyses led to the conclusion that the mitochondrial genome experienced dynamic evolution in genome size, gene content, intron acquisition, gene order, and RNA editing during the origins of land plants and their major clades. However, evolution of this organellar genome has remained rather conservative since the origin and initial radiation of early land plants, except within vascular plants.

Identification and characterization of plastid trnF(GAA) pseudogenes in four species of Solanum (Solanaceae)

Non-functional trnF pseudogenes that rarely occur in embryophytes have been found in Solanaceae. We have sequenced the trnL-F intergenic spacer of four species of Solanum, and found duplicated regions of the original trnF gene. These repeats were 94-260 bp long causing large length variation in the trnL-F intergenic spacer resulting from differences in pseudogene copy number (2-4). The duplicated trnF regions are comprised of several highly structured motifs, which were partial residues, or entire parts of the Anticodon, T- and D-domains of the original gene, but all lacked the acceptor stems at the 5'- or 3'-end. Pseudogenes included several transitions and transversions in their sequences compared to the original trnF gene. Among pseudogene copies, T-domains were more frequent and fragmented than D-domain elements. Our results demonstrate that although chloroplast evolution is uniform such structural duplications in the sequences used for phylogenetic reconstructions should be treated with great caution.

Complete chloroplast DNA sequence from a Korean endemic genus, Megaleranthis saniculifolia, and its evolutionary implications

The chloroplast DNA sequences of Megaleranthis saniculifolia, an endemic and monotypic endangered plant species, were completed in this study (GenBank FJ597983). The genome is 159,924 bp in length. It harbors a pair of IR regions consisting of 26,608 bp each. The lengths of the LSC and SSC regions are 88,326 bp and 18,382 bp, respectively. The structural organizations, gene and intron contents, gene orders, AT contents, codon usages, and transcription units of the Megaleranthis chloroplast genome are similar to those of typical land plant cp DNAs. However, the detailed features of Megaleranthis chloroplast genomes are substantially different from that of Ranunculus, which belongs to the same family, the Ranunculaceae. First, the Megaleranthis cp DNA was 4,797 bp longer than that of Ranunculus due to an expanded IR region into the SSC region and duplicated sequence elements in several spacer regions of the Megaleranthis cp genome. Second, the chloroplast genomes of Megaleranthis and Ranunculus evidence 5.6% sequence divergence in the coding regions, 8.9% sequence divergence in the intron regions, and 18.7% sequence divergence in the intergenic spacer regions, respectively. In both the coding and noncoding regions, average nucleotide substitution rates differed markedly, depending on the genome position. Our data strongly implicate the positional effects of the evolutionary modes of chloroplast genes. The genes evidencing higher levels of base substitutions also have higher incidences of indel mutations and low Ka/Ks ratios. A total of 54 simple sequence repeat loci were identified from the Megaleranthis cp genome. The existence of rich cp SSR loci in the Megaleranthis cp genome provides a rare opportunity to study the population genetic structures of this endangered species. Our phylogenetic trees based on the two independent markers, the nuclear ITS and chloroplast matK sequences, strongly support the inclusion of the Megaleranthis to the Trollius. Therefore, our molecular trees support Ohwi's original treatment of Megaleranthis saniculiforia to Trollius chosenensis Ohwi.

Complete nucleotide sequence of the Cryptomeria japonica D. Don. chloroplast genome and comparative chloroplast genomics: diversified genomic structure of coniferous species

BACKGROUND

The recent determination of complete chloroplast (cp) genomic sequences of various plant species has enabled numerous comparative analyses as well as advances in plant and genome evolutionary studies. In angiosperms, the complete cp genome sequences of about 70 species have been determined, whereas those of only three gymnosperm species, Cycas taitungensis, Pinus thunbergii, and Pinus koraiensis have been established. The lack of information regarding the gene content and genomic structure of gymnosperm cp genomes may severely hamper further progress of plant and cp genome evolutionary studies. To address this need, we report here the complete nucleotide sequence of the cp genome of Cryptomeria japonica, the first in the Cupressaceae sensu lato of gymnosperms, and provide a comparative analysis of their gene content and genomic structure that illustrates the unique genomic features of gymnosperms.

RESULTS

The C. japonica cp genome is 131,810 bp in length, with 112 single copy genes and two duplicated (trnI-CAU, trnQ-UUG) genes that give a total of 116 genes. Compared to other land plant cp genomes, the C. japonica cp has lost one of the relevant large inverted repeats (IRs) found in angiosperms, fern, liverwort, and gymnosperms, such as Cycas and Gingko, and additionally has completely lost its trnR-CCG, partially lost its trnT-GGU, and shows diversification of accD. The genomic structure of the C. japonica cp genome also differs significantly from those of other plant species. For example, we estimate that a minimum of 15 inversions would be required to transform the gene organization of the Pinus thunbergii cp genome into that of C. japonica. In the C. japonica cp genome, direct repeat and inverted repeat sequences are observed at the inversion and translocation endpoints, and these sequences may be associated with the genomic rearrangements.

CONCLUSION

The observed differences in genomic structure between C. japonica and other land plants, including pines, strongly support the theory that the large IRs stabilize the cp genome. Furthermore, the deleted large IR and the numerous genomic rearrangements that have occurred in the C. japonica cp genome provide new insights into both the evolutionary lineage of coniferous species in gymnosperm and the evolution of the cp genome.

Dynamics and evolution of the inverted repeat-large single copy junctions in the chloroplast genomes of monocots

Background

Various expansions or contractions of inverted repeats (IRs) in chloroplast genomes led to fluxes in the IR-LSC (large single copy) junctions. Previous studies revealed that some monocot IRs contain a trnH-rps19 gene cluster, and it has been speculated that this may be an evidence of a duplication event prior to the divergence of monocot lineages. Therefore, we compared the organizations of genes flanking two IR-LSC junctions in 123 angiosperm representatives to uncover the evolutionary dynamics of IR-LSC junctions in basal angiosperms and monocots.

Results

The organizations of genes flanking IR-LSC junctions in angiosperms can be classified into three types. Generally each IR of monocots contains a trnH-rps19 gene cluster near the IR-LSC junctions, which differs from those in non-monocot angiosperms. Moreover, IRs expanded more progressively in monocots than in non-monocot angiosperms. IR-LSC junctions commonly occurred at polyA tract or A-rich regions in angiosperms. Our RT-PCR assays indicate that in monocot IR_A the trnH-rps19 gene cluster is regulated by two opposing promoters, S10_A and psbA.

Conclusion

Two hypotheses are proposed to account for the evolution of IR expansions in monocots. Based on our observations, the inclusion of a trnH-rps19 cluster in majority of monocot IRs could be reasonably explained by the hypothesis that a DSB event first occurred at IR_B and led to the expansion of IRs to trnH, followed by a successive DSB event within IR_A and lead to the expansion of IRs to rps19 or to rpl22 so far. This implies that the duplication of trnH-rps19 gene cluster was prior to the diversification of extant monocot lineages. The duplicated trnH genes in the IR_B of most monocots and non-monocot angiosperms have distinct fates, which are likely regulated by different expression levels of S10_A and S10_B promoters. Further study is needed to unravel the evolutionary significance of IR expansion in more recently diverged monocots.

The plastid sigma factor SIG1 maintains photosystem I activity via regulated expression of the psaA operon in rice chloroplasts

Sigma factors encoded by the nucleus of plants confer promoter specificity on the bacterial-type RNA polymerase in chloroplasts. We previously showed that transcripts of OsSIG1, which encodes one such sigma factor in rice, accumulate relatively late during leaf development. We have now isolated and characterized two allelic mutants of OsSIG1, in which OsSIG1 is disrupted by insertion of the retrotransposon Tos17, in order to characterize the functions of OsSIG1. The OsSIG1-/- plants were found to be fertile but they manifested an approximately one-third reduction in the chlorophyll content of mature leaves. Quantitative RT-PCR and northern blot analyses of chloroplast gene expression revealed that the abundance of transcripts derived from the psaA operon was markedly reduced in OsSIG1-/- plants compared with that in wild-type homozygotes. This effect was accompanied by a reduction in the abundance of the core protein complex (PsaA-PsaB) of photosystem I. Analysis of chlorophyll fluorescence also revealed a substantial reduction in the rate of electron transfer from photosystem II to photosystem I in the OsSIG1 mutants. Our results thus indicate that OsSIG1 plays an important role in the maintenance of photosynthetic activity in mature chloroplasts of rice by regulating expression of chloroplast genes for components of photosystem I.

Expression of matK: functional and evolutionary implications

Strong phylogenetic signal from matK has rendered it an invaluable gene in plant systematic and evolutionary studies at various evolutionary depths. Further, matK is proposed as the only chloroplast-encoded group II intron maturase, thus implicating MATK in chloroplast posttranscriptional processing. For a protein-coding gene, matK has an unusual evolutionary mode and tempo, including relatively high substitution rates at both the nucleotide and amino acids levels. These evolutionary features have raised questions about matK function. In the current study, we examined matK RNA and protein from representative land plant species to provide insight into functional aspects of this unusual gene. We report the first evidence of a transcript for matK separate from the trnK precursor and demonstrate that a full-length MATK protein exists in five angiosperm species. We also show that matK RNA and protein levels are regulated by light and developmental stage, suggesting functional roles for this putative maturase. Specifically, matK expression increased after etiolation and decreased at 4 weeks after germination. This work provides evidence for the expression of the only putative chloroplast-encoded group II intron maturase and insight into regulation mechanisms relating to plant development and, indirectly, to photosynthesis.

Supernetwork Identifies Multiple Events of Plastid trnF(GAA) Pseudogene Evolution in the Brassicaceae

The occurrence of nonfunctional trnF pseudogenes has been rarely described in flowering plants. However, we describe the first large-scale supernetwork for the Brassiccaeae built from gene trees for 5 loci (adh, chs, matK, trnL-F, and ITS) and report multiple independent origins for trnF pseudogenes in crucifers. The duplicated regions of the original trnF gene are comprised of its anticodon domain and several other highly structured motifs not related to the original gene. Length variation of the trnL-F intergenic spacer region in different taxa ranges from 219 to 900 bp as a result of differences in pseudocopy number (1-14). It is speculated that functional constraints favor 2-3 or 5-6 copies, as found in Arabidopsis and Boechera. The phylogenetic distribution of microstructural changes for the trnL-F region supports ancient patterns of divergence in crucifer evolution for some but not all gene loci.

Salt stress response in rice: genetics, molecular biology, and comparative genomics

Significant progress has been made in unraveling the molecular biology of rice in the past two decades. Today, rice stands as a forerunner amongst the cereals in terms of details known on its genetics. Evidence show that salt tolerance in plants is a quantitative trait. Several traditional cultivars, landraces, and wild types of rice like Pokkali, CSR types, and Porteresia coarctata appear as promising materials for donation of requisite salt tolerance genes. A large number of quantitative trait loci (QTL) have been identified for salt tolerance in rice through generation of recombinant inbred lines and are being mapped using different types of DNA markers. Salt-tolerant transgenic rice plants have been produced using a host of different genes and transcript profiling by micro- and macroarray-based methods has opened the gates for the discovery of novel salt stress mechanisms in rice, and comparative genomics is turning out to be a critical input in this respect. In this paper, we present a comprehensive review of the genetic, molecular biology, and comparative genomics effort towards the generation of salt-tolerant rice. From the data on comprehensive transcript expression profiling of clones representing salt-stress-associated genes of rice, it is shown that transcriptional and translational machineries are important determinants in controlling salt stress response, and gene expression response in tolerant and susceptible rice plants differs mainly in quantitative terms.

The chloroplast trnT-trnF region in the seed plant lineage Gnetales

The trnT-trnF region is located in the large single-copy region of the chloroplast genome. It consists of the trnL intron, a group I intron, and the trnT-trnL and trnL-trnF intergenic spacers. We analyzed the evolution of the region in the three genera of the gymnosperm lineage Gnetales (Gnetum, Welwitschia, and Ephedra), with especially dense sampling in Gnetum for which we sequenced 41 accessions, representing most of the 25-35 species. The trnL intron has a conserved secondary structure and contains elements that are homologous across land plants, while the spacers are so variable in length and composition that homology cannot be found even among the three genera. Palindromic sequences that form hairpin structures were detected in the trnL-trnF spacer, but neither spacer contained promoter elements for the tRNA genes. The absence of promoters, presence of hairpin structures in the trnL-trnF spacer, and high sequence variation in both spacers together suggest that trnT and trnF are independently transcribed. Our model for the expression and processing of the genes tRNA(Thr)(UGU), tRNA(Leu)(UAA), and tRNA(Phe) (GAA) therefore attributes the seemingly neutral evolution of the two spacers to their escape from functional constraints.

Evolution of the trnF(GAA) Gene in Arabidopsis Relatives and the Brassicaceae Family: Monophyletic Origin and Subsequent Diversification of a Plastidic Pseudogene

Recently, we used the 5'-trnL(UAA)-trnF(GAA) region of the chloroplast DNA for phylogeographic reconstructions and phylogenetic analysis among the genera Arabidopsis, Boechera, Rorippa, Nasturtium, and Cardamine. Despite the fact that extensive gene duplications are rare among the chloroplast genome of higher plants, within these taxa the anticodon domain of the trnF(GAA) gene exhibit extensive gene duplications with one to eight tandemly repeated copies in close 5' proximity of the functional gene. Interestingly, even in Arabidopsis thaliana we found six putative pseudogenic copies of the functional trnF gene within the 5'-intergenic trnL-trnF spacer. A reexamination of trnL(UAA)-trnF(GAA) regions from numerous published phylogenetic studies among halimolobine, cardaminoid, and other cruciferous taxa revealed not only extensive trnF gene duplications but also favor the hypothesis about a single origin of trnF pseudogene formation during evolution of the Brassicaceae family 16-21 MYA. Conserved sequence motifs from this tandemly repeated region are codistributed nonrandomly throughout the plastome, and we found some similarities with a DNA sequence duplication in the rps7 gene and its adjacent spacer. Our results demonstrate the potential evolutionary dynamics of a plastidic region generally regarded as highly conserved and probably cotranscribed and, as shown here for several genera among cruciferous plants, greatly characterized by parallel gains and losses of duplicated trnF copies.

Molecular evolution of the trnTUGU-trnFGAA region in Bryophytes

Structure, variability, and molecular evolution of the trnT-F region in the Bryophyta (mosses and liverworts) is analyzed based on about 200 sequences of the trnT-L spacer and trnL 5' exon, 1000 sequences of the trnL intron, and 800 sequences of the trnL 3' exon and trnL-F spacer, including comparisons of lengths, GC contents, sequence similarities, and functional elements. Mutations occurring in the trnL 5' and 3' exons, including compensatory base pair changes, and a transition in the trnL anticodon in Takakia lepidozioides, are discussed. All three non-coding regions display a mosaic structure of highly variable elements (V1 - V3 in the trnT-L spacer, V4/V5 corresponding to stem-loop regions P6/P8 in the trnL intron, and V6/V7 in the trnL-F spacer) and more conserved elements. In the trnL intron this structure is a consequence of the defined secondary structure necessary for correct splicing, whereas in both spacers conserved regions are restricted to promoter elements. At least the highly variable regions in the trnT-L spacer and stem-loop region P8 of the trnL intron seem to evolve independently in the major bryophyte lineages and are therefore not suitable for high taxonomic level phylogenetic reconstructions. In mosses, a trend of length reduction towards the more derived lineages is observed in all three non-coding regions. GC contents are mostly linked to sequence variability, with the conserved regions being more GC rich and the more variable AT rich. The lowest GC values (< 10 %) are found in the trnT-L spacer of mosses. In addition to two putative sigma (70)-type promoters in the trnT-L spacer, a third putative promoter is present in the trnL-F spacer, although trnL and trnF are assumed to be co-transcribed. Consensus sequences are provided for the -35 and -10 sequences of the major bryophyte lineages. The third promoter is part of a hairpin secondary structure, whose loop region is highly homoplastic in mosses due to an inversion occurring independently in non-related taxa, even at the intraspecific level.

A comparison of rice chloroplast genomes

Using high quality sequence reads extracted from our whole genome shotgun repository, we assembled two chloroplast genome sequences from two rice (Oryza sativa) varieties, one from 93-11 (a typical indica variety) and the other from PA64S (an indica-like variety with maternal origin of japonica), which are both parental varieties of the super-hybrid rice, LYP9. Based on the patterns of high sequence coverage, we partitioned chloroplast sequence variations into two classes, intravarietal and intersubspecific polymorphisms. Intravarietal polymorphisms refer to variations within 93-11 or PA64S. Intersubspecific polymorphisms were identified by comparing the major genotypes of the two subspecies represented by 93-11 and PA64S, respectively. Some of the minor genotypes occurring as intravarietal polymorphisms in one variety existed as major genotypes in the other subspecific variety, thus giving rise to intersubspecific polymorphisms. In our study, we found that the intersubspecific variations of 93-11 (indica) and PA64S (japonica) chloroplast genomes consisted of 72 single nucleotide polymorphisms and 27 insertions or deletions. The intersubspecific polymorphism rates between 93-11 and PA64S were 0.05% for single nucleotide polymorphisms and 0.02% for insertions or deletions, nearly 8 and 10 times lower than their respective nuclear genomes. Based on the total number of nucleotide substitutions between the two chloroplast genomes, we dated the divergence of indica and japonica chloroplast genomes as occurring approximately 86,000 to 200,000 years ago.

Noncoding plastid trnT-trnF sequences reveal a well resolved phylogeny of basal angiosperms

Recent contributions from DNA sequences have revolutionized our concept of systematic relationships in angiosperms. However, parts of the angiosperm tree remain unclear. Previous studies have been based on coding or rDNA regions of relatively conserved genes. A phylogeny for basal angiosperms based on noncoding, fast-evolving sequences of the chloroplast genome region trnT-trnF is presented. The recognition of simple direct repeats allowed a robust alignment. Mutational hot spots appear to be confined to certain sectors, as in two stem-loop regions of the trnL intron secondary structure. Our highly resolved and well-supported phylogeny depicts the New Caledonian Amborella as the sister to all other angiosperms, followed by Nymphaeaceae and an Austrobaileya-Illicium-Schisandra clade. Ceratophyllum is substantiated as a close relative of monocots, as is a monophyletic eumagnoliid clade consisting of Piperales plus Winterales sister to Laurales plus Magnoliales. Possible reasons for the striking congruence between the trnT-trnF based phylogeny and phylogenies generated from combined multi-gene, multi-genome data are discussed.

Post-transcriptional control of chloroplast gene expression. Accumulation of stable psaC mRNA is due to downstream RNA cleavages in the ndhD gene

Intergenic cleavages, intron splicing, and editing of primary transcripts of the plastid ndhH-D operon produce multiple overlapping RNAs, of which the most abundant by far is the monocistronic 400-nucleotide mRNA of psaC (encoding the PsaC protein of photosystem I), in contrast with the low level of transcripts of the six ndh genes. Like other plastid operons containing genes for functionally unrelated proteins, the contrasting accumulation of ndh and psaC transcripts provides a model to investigate the mechanisms of the post-transcriptional control of gene expression, a feature of chloroplast genetic machinery, with a minimum of interference by transcriptional control. In leek (Allium porrum L), the ndhD transcript (which follows the psaC gene and ends the ndhH-D operon) requires C --> U editing to restore its start codon and may be used as a marker for the processing of psaC and ndhD transcripts. By determining the editing state and 5' end sequences of specific transcripts, we demonstrated that stable monocistronic psaC mRNA results from downstream cleavages in the ndhD sequence, which renders non-functional ndhD transcripts as by-products. Alternative psaC-ndhD intergenic cleavages produce complete mRNAs for both genes, but only take place in precursors containing editing-restored ndhD start codons. Hence, post-transcriptional control acts by promoting the ndhD cleavage alternative, which allows the accumulation of psaC mRNA at the expense of ndhD mRNA levels.

Chloroplast ribosomal S14 protein transcript is edited to create a translation initiation codon in the moss Physcomitrella patens

The rps14 transcript is edited in the moss Physcomitrella patens chloroplast by a C-to-U transition, to create a translation initiation codon, AUG. The efficiency of RNA editing was low, with approximately 20% of rps14 transcripts edited. This suggests that the translation of rps14 mRNA is strictly regulated by RNA editing. This is the first report of RNA editing in P. patens and the creation of a translation initiation codon in rps14 mRNA in chloroplasts.

An Arabidopsis sigma factor (SIG2)-dependent expression of plastid-encoded tRNAs in chloroplasts

A eubacteria-type RNA polymerase (PEP) plays crucial roles for chloroplast development in higher plants. The core subunits are encoded on plastid DNA (rpo genes) while the regulatory sigma factors are encoded on the nuclear DNA (SIG genes). However, the definite gene specificity of each sigma factor is unknown. We recently identified an Arabidopsis recessive pale-green mutant abc1 in which T-DNA is inserted in SIG2 (sigB). In this mutant, almost normal etioplasts were developed under dark conditions while the small chloroplasts with poor thylakoid membranes and stacked lamellar were developed under light conditions. The sig2-1 mutant was deficient in accumulating enough photosynthetic and photosynthesis-related proteins as well as chlorophyll. However, mRNAs of their structural genes were not significantly reduced. Further analyses revealed that several plastid-encoded tRNAs including trnE-UUC that has dual function for protein and ALA biosyntheses were drastically reduced in the sig2-1 mutant. In contrast, nucleus-encoded T7 phage-type RNA polymerase (NEP)-dependent gene transcripts were steadily accumulated in the mutant. These results indicate that progress of chloroplast development requires SIG2-dependent expression of plastid genes, particularly some of the tRNA genes.

Transcripts of the ndhH-D operon of barley plastids: possible role of unedited site III in splicing of the ndhA intron

The plastid ndhH-D operon produces several transcripts containing ndhA sequence with and without its group II intron. After sequencing an 8125 bp fragment of barley plastid DNA including the ndhH-D operon, we investigated the editing-splicing status of transcripts in the range 1.0-7.8 kb. Reverse transcription and sequencing of RNA bands separated by electrophoresis were used to determine C-->U editing sites. Sites I, II and IV of ndhA and site V of ndhD were edited in all transcripts analysed and, probably, were edited before any splicing had taken place. In contrast, site III of ndhA (13 bp from the 5'-end base of the second exon) was not edited in transcripts containing the intron (including the 1.7 kb intermediary transcript consisting of the intron and the second exon) but was edited in all transcripts lacking the ndhA intron. Comparison of the secondary structures of the ndhA intron and intron-second exon intermediate suggests that G pairing prevents editing of site III in transcripts containing the intron and maintains the secondary structure required for splicing. Splicing of the ndhA intron releases the site III C from pairing and, probably, brings it close to cis-acting elements for editing upstream in the first exon.

A chloroplast-derived sequence is utilized as a source of promoter sequences for the gene for subunit 9 of NADH dehydrogenase (nad9) in rice mitochondria

The chloroplast-derived sequence trnS-rps4/ 3'trnL-trnF-ndhJ-ndhK (4066 bases in length) is present in a region that starts 355 bases upstream of the gene for subunit 9 of NADH dehydrogenase (nad9) in the mitochondrial genome of rice. Northern blot hybridization revealed that three large transcripts of 3.05, 1.62 and 1.05 kb hybridized to strand-specific probes for both the nad9 gene and the chloroplast-derived sequence, indicating that the nad9 gene was transcribed together with the chloroplast-derived sequence. From the results of in vitro capping and ribonuclease protection experiments, as well as primer extension analysis, we identified at least seven sites for the initiation of transcription of nad9 in the chloroplast-derived sequence. All of the initiation sites for transcription of the nad9 gene were located in sequences homologous to chloroplast DNA. Two of seven initiation sites were flanked by a sequence homologous to the consensus promoter motif that includes the CRTA motif (where R is A or G) of the rice mitochondrion. However, the sequences surrounding the other five sites showed only limited similarity to the conserved sequence. It is suggested that all the promoters of the rice nad9 gene exist in a sequence that was transferred from the chloroplast during evolution. Thus, the chloroplast-derived sequence has a novel, significant function in the mitochondrial genome of this higher plant.

Impaired splicing of the rps12 transcript in ribosome-deficient plastids

Analysis of RNA maturation in ribosome-deficient plastids of four non-allelic barley mutants revealed an increased accumulation and altered processing of transcripts of the ribosomal protein gene CS12 (rps12) compared to normal chloroplasts. The three exons of rps12 are part of two different polycistronic transcription units. Generation of mature rps12-mRNA involves both cis- and trans-splicing. In ribosome-deficient plastids, the cis-intron separating exons 2 and 3 remains entirely unspliced whereas the splicing of the bipartite rps12 trans-intron between exon 1 and exon 2 occurs, but at a reduced level. A comparison of the 3' and 5' ends of the two RNAs that are generally assumed to interact during trans-splicing showed a difference in the processing pathways of 3' rps12 transcripts between mutated and normal chloroplasts. Nonetheless, the final products were identical.

Transcription, splicing and editing of plastid RNAs in the nonphotosynthetic plant Epifagus virginiana

Expression of the vestigial plastid genome of the nonphotosynthetic, parasitic flowering plant Epifagus virginiana was examined by northern analysis and by characterization of cDNAs. Probes for each of 12 plastid genes tested hybridized to all lanes of northern blots containing total RNA prepared from stems and fruits of Epifagus and from leaves of tobacco. Certain transcript patterns in Epifagus plastids are highly complex and similar to those of tobacco operons. In contrast, genes such as rps2, which have become orphaned in Epifagus as a result of evolutionary loss of formerly cotranscribed genes, show simpler transcript patterns in Epifagus than in tobacco. Sizing and sequencing of cDNAs generated by reverse transcriptase-PCR for three genes, rps12, rpl2, and clpP, show that their transcripts are properly cis- and/or trans-spliced at the same five group II intron insertion sites used in photosynthetic plants. A single, conventional C-->U edit in rps12 was found among the total of 1401 nucleotides of cDNA sequence that was determined for the three genes. An octanucleotide sequence identical to a putative guide RNA of plant organelles and perfectly complementary to the rps12 edit site itself was identified just 200 bp upstream of the edit site. These data, together with previous results from the complete sequencing of the Epifagus plastid genome, provide compelling evidence that this degenerate genome is nonetheless expressed and functional. Analysis of the putative maturase MatK, encoded by the group II intron of trnK in photosynthetic land plants but by a freestanding gene in Epifagus, leads us to hypothesize that it acts 'in trans' to assist the splicing of group II introns other than the one in which it is normally encoded.

ndhF sequence evolution and the major clades in the sunflower family

An extensive sequence comparison of the chloroplast ndhF gene from all major clades of the largest flowering plant family (Asteraceae) shows that this gene provides approximately 3 times more phylogenetic information than rbcL. This is because it is substantially longer and evolves twice as fast. The 5' region (1380 bp) of ndhF is very different from the 3' region (855 bp) and is similar to rbcL in both the rate and the pattern of sequence change. The 3' region is more A+T-rich, has higher levels of nonsynonymous base substitution, and shows greater transversion bias at all codon positions. These differences probably reflect different functional constraints on the 5' and 3' regions of ndhF. The two patterns of base substitutions of ndhF are particularly advantageous for phylogenetic reconstruction because the conserved and variable segments can be used for older and recent groups, respectively. Phylogenetic analyses of 94 ndhF sequences provided much better resolution of relationships than previous molecular and morphological phylogenies of the Asteraceae. The ndhF tree identified five major clades: (i) the Calyceraceae is the sister family of Asteraceae; (ii) the Barnadesioideae is monophyletic and is the sister group to the rest of the family; (iii) the Cichorioideae and its two basal tribes Mutisieae and Cardueae are paraphyletic; (iv) four tribes of Cichorioideae (Lactuceae, Arctoteae, Liabeae, and Vernonieae) form a monophyletic group, and these are the sister clade of the Asteroideae; and (v) the Asteroideae is monophyletic and includes three major clades.

An atpE-specific promoter within the coding region of the atpB gene in tobacco chloroplast DNA

The atpB and atpE genes encode beta and epsilon subunits, respectively, of chloroplast ATP synthase and are co-transcribed in the plant species so far studied. In tobacco, an atpB gene-specific probe hybridizes to 2.7- and 2.3-kb transcripts. In addition to these, a probe from the atpE coding region hybridizes also to a 1.0-kb transcript. The 5' end of the atpE-specific transcript has been mapped 430/431 nt upstream of the atpE translation initiation site, within the coding region of the atpB gene. In-vitro capping revealed that this transcript results from a primary transcriptional event and is also characterized by -10 and -35 canonical sequences in the 5' region. It has been found to share a common 3' end with the bi-cistronic transcripts that has been mapped within the coding region of the divergently transcribed trnM gene, approximately 236 nt downstream from the atpE termination codon. Interestingly, this transcript accumulates only in leaves and not in proplastid-containing cultured (BY-2) cells, indicating that, unless it is preferentially degraded in BY-2 cells, its expression might be transcriptionally controlled.

Developmental and light-dependent cues interact to establish steady-state levels of transcripts for photosynthesis-related genes (psbA, psbD, psaA and rbcL) in rice (Oryza sativa L.)

The steady-state transcript levels for psbA, psbD, psaA and rbcL are low in dark-grown rice seedlings as compared to those grown in light. Following seed germination, they accumulate in an age-dependent manner, in dark as well as light, reaching a maximal level on the 7th or 8th day, before a slow decline sets in. But transcripts for psbA and psbD continue to maintain relatively-high levels even after 10 days of growth in light. Exposure of 5-day-old dark-grown seedlings to light results in an approximately 25-60-fold increase in transcripts during a period of 72 h, followed by a decrease. An analysis of data from both lines of investigation reveals that the developmental programme increases the transcript levels for psbA, psbD, psaA and rbcL by about 10-, 2.3-, 7.0- and 8.0-fold, respectively, between 5-8 days after germination and it is independent of light. At the same time, exposure of the seedlings to light during this period further enhances transcript levels by 5-, 11.4-, 6.6- and 7.8-fold, respectively. Thus, both developmental and light-dependent cues contribute to establish steady-state levels of transcripts for the chloroplast genes investigated.

Separation of two classes of plastid DNA-dependent RNA polymerases that are differentially expressed in mustard (Sinapis alba L.) seedlings

Chloroplast and etioplast in vitro transcription systems from mustard have different functional properties, which is reflected in differences in phosphorylation status. Here we report another transcription control mechanism, which involves two plastid DNA-dependent RNA polymerases designated as peak A and peak B enzymes. Both are large multi-subunit complexes, but differ in their native molecular mass (> 700 kDa for peak A and ca. 420 kDa for peak B) and in their polypeptide composition. The A enzyme is composed of at least 13 polypeptides, while the B enzyme contains only four putative subunits. Peak B activity is inhibited by rifampicin, whereas that of peak A is resistant. RNA polymerase activity was compared for plastids from cotyledons of 4-day-old seedlings that were grown either under continuous light (chloroplasts) or in darkness (etioplasts), or were first dark-grown and then transferred to light for 16 h ('intermediate-type' plastids). While the total activity was approximately the same in all three cases, enzyme B was the predominant activity obtained from etioplasts and enzyme A that obtained from chloroplasts. Both had equal activity in preparations from the 'intermediate-type' plastid form. Both activation/inactivation and differential gene expression seem to play a role in the regulation of the plastid transcription machinery.

Variations in chloroplast DNA from rice (Oryza sativa): differences between deletions mediated by short direct-repeat sequences within a single species

In a previous study, we compared chloroplast DNAs (ctDNAs) from four species ofOryza and detected two independent deletions of DNA fragments in the ctDNAs (Kanno and Hirai 1992a). These deletions were genotype-specific variations. Since short direct-repeat sequences were detected at the borders of both deletions, the deletions were apparently the result of intramolecular recombination mediated by these direct-repeat sequences. In the present study, we examined whether or not this type of variation exists within a single species. Ishii et al. demonstrated three types of ctDNA inO. Sativa (1988), and the source of the variations that they identified seemed to be deletions. We determined the precise locations of the deletions and the sequences around them. As expected, our results showed that these variations were the results of deletions that were mediated by short direct-repeat sequences. While the deletions that had been found previously were located on spacer regions, those found in this study were located within open reading frames (ORFs). Northern hybridization analysis showed that one of the ORFs was-transcribed. In the case of this deletion, the amino acid sequence encoded by the C-terminal region of the ORF was altered and the short inverted-repeat sequences downstream of the ORF were deleted. In addition, there were other short inverted-repeat sequences downstream of the altered ORF.