The complete nucleotide sequence of the tobacco chloroplast genome: its gene organization and expression

Rpl33, a nonessential plastid-encoded ribosomal protein in tobacco, is required under cold stress conditions

Plastid genomes contain a conserved set of genes encoding components of the translational apparatus. While knockout of plastid translation is lethal in tobacco (Nicotiana tabacum), it is not known whether each individual component of the plastid ribosome is essential. Here, we used reverse genetics to test whether several plastid genome-encoded ribosomal proteins are essential. We found that, while ribosomal proteins Rps2, Rps4, and Rpl20 are essential for cell survival, knockout of the gene encoding ribosomal protein Rpl33 did not affect plant viability and growth under standard conditions. However, when plants were exposed to low temperature stress, recovery of Rpl33 knockout plants was severely compromised, indicating that Rpl33 is required for sustaining sufficient plastid translation capacity in the cold. These findings uncover an important role for plastid translation in plant tolerance to chilling stress.

A plastid gene phylogeny of the non-photosynthetic parasitic Orobanche (Orobanchaceae) and related genera

The phylogenetic relationships of the non-photosynthetic Orobanche sensu lato (Orobanchaceae), which includes some of the economically most important parasitic weeds, remain insufficiently understood and controversial. This concerns both the phylogenetic relationships within the genus, in particular its monophyly or lack thereof, and the relationships to other holoparasitic genera such as Cistanche or Conopholis. Here we present the first comprehensive phylogenetic study of this group based on a region from the plastid genome (rps2 gene). Although substitution rates appear to be elevated compared to the photosynthetic members of Orobanchaceae, relationships among the major lineages Cistanche, Conopholis plus Epifagus, Boschniakia rossica (Cham. & Schltdl.) B. Fedtsch., B. himalaica Hook. f. & Thomson, B. hookeri Walp. plus B. strobilacea A. Gray, and Orobanche s. l. remain unresolved. Resolution within Orobanche, however, is much better. In agreement with morphological, cytological and other molecular phylogenetic evidence, five lineages, corresponding to the four traditionally recognised sections (Gymnocaulis, Myzorrhiza, Orobanche, Trionychon) and O. latisquama Reut. ex Boiss. (of sect. Orobanche), can be distinguished. A combined analysis of plastid rps2 and nuclear ITS sequences of the holoparasitic genera results in more resolved and better supported trees, although the relationships among Orobanche s. l., Cistanche, and the clade including the remaining genera is unresolved. Therefore, rps2 is a marker from the plastid genome that is well-suited to be used in combination with other already established nuclear markers for resolving generic relationships of Orobanche and related genera.

NDF6: a thylakoid protein specific to terrestrial plants is essential for activity of chloroplastic NAD(P)H dehydrogenase in Arabidopsis

NAD(P)H dehydrogenase (NDH) is a homolog of respiratory complex I and mediates one of the two pathways of cyclic electron flow around PSI (CEF I). Although 15 ndh subunits have been identified in the chloroplastic and nuclear genomes of higher plants, no electron accepter subunits have been identified to date. To identify the missing chloroplastic NDH subunits, we undertook an in silico approach based on co-expression analysis. In this report, we characterized the novel gene NDF6 (NDH-dependent flow 6; At1g18730) which encodes a protein that is essential for NDH activity. NDF6 has one transmembrane domain and is localized in the thylakoid membrane fraction. Homologous proteins of NDF6 were identified in the genomes of terrestrial plants; however, no homologs have been found in cyanobacteria, which are thought to be the origin of chloroplasts and have a minimal NDH complex unit. NDF6 is unstable in ndhB-impaired or disrupted mutants of higher plants in which the chloroplastic NDH complex is thought to be degraded. These results suggest that NDF6 is a novel subunit of chloroplastic NDH that was added to terrestrial plants during evolution.

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.

Comparative chloroplast genomics and phylogenetics of Fagopyrum esculentum ssp. ancestrale -a wild ancestor of cultivated buckwheat

BACKGROUND

Chloroplast genome sequences are extremely informative about species-interrelationships owing to its non-meiotic and often uniparental inheritance over generations. The subject of our study, Fagopyrum esculentum, is a member of the family Polygonaceae belonging to the order Caryophyllales. An uncertainty remains regarding the affinity of Caryophyllales and the asterids that could be due to undersampling of the taxa. With that background, having access to the complete chloroplast genome sequence for Fagopyrum becomes quite pertinent.

RESULTS

We report the complete chloroplast genome sequence of a wild ancestor of cultivated buckwheat, Fagopyrum esculentum ssp. ancestrale. The sequence was rapidly determined using a previously described approach that utilized a PCR-based method and employed universal primers, designed on the scaffold of multiple sequence alignment of chloroplast genomes. The gene content and order in buckwheat chloroplast genome is similar to Spinacia oleracea. However, some unique structural differences exist: the presence of an intron in the rpl2 gene, a frameshift mutation in the rpl23 gene and extension of the inverted repeat region to include the ycf1 gene. Phylogenetic analysis of 61 protein-coding gene sequences from 44 complete plastid genomes provided strong support for the sister relationships of Caryophyllales (including Polygonaceae) to asterids. Further, our analysis also provided support for Amborella as sister to all other angiosperms, but interestingly, in the bayesian phylogeny inference based on first two codon positions Amborella united with Nymphaeales.

CONCLUSION

Comparative genomics analyses revealed that the Fagopyrum chloroplast genome harbors the characteristic gene content and organization as has been described for several other chloroplast genomes. However, it has some unique structural features distinct from previously reported complete chloroplast genome sequences. Phylogenetic analysis of the dataset, including this new sequence from non-core Caryophyllales supports the sister relationship between Caryophyllales and asterids.

The complete plastid genome sequence of Welwitschia mirabilis: an unusually compact plastome with accelerated divergence rates

Background

Welwitschia mirabilis is the only extant member of the family Welwitschiaceae, one of three lineages of gnetophytes, an enigmatic group of gymnosperms variously allied with flowering plants or conifers. Limited sequence data and rapid divergence rates have precluded consensus on the evolutionary placement of gnetophytes based on molecular characters. Here we report on the first complete gnetophyte chloroplast genome sequence, from Welwitschia mirabilis, as well as analyses on divergence rates of protein-coding genes, comparisons of gene content and order, and phylogenetic implications.

Results

The chloroplast genome of Welwitschia mirabilis [GenBank: EU342371] is comprised of 119,726 base pairs and exhibits large and small single copy regions and two copies of the large inverted repeat (IR). Only 101 unique gene species are encoded. The Welwitschia plastome is the most compact photosynthetic land plant plastome sequenced to date; 66% of the sequence codes for product. The genome also exhibits a slightly expanded IR, a minimum of 9 inversions that modify gene order, and 19 genes that are lost or present as pseudogenes. Phylogenetic analyses, including one representative of each extant seed plant lineage and based on 57 concatenated protein-coding sequences, place Welwitschia at the base of all seed plants (distance, maximum parsimony) or as the sister to Pinus (the only conifer representative) in a monophyletic gymnosperm clade (maximum likelihood, bayesian). Relative rate tests on these gene sequences show the Welwitschia sequences to be evolving at faster rates than other seed plants. For these genes individually, a comparison of average pairwise distances indicates that relative divergence in Welwitschia ranges from amounts about equal to other seed plants to amounts almost three times greater than the average for non-gnetophyte seed plants.

Conclusion

Although the basic organization of the Welwitschia plastome is typical, its compactness, gene content and high nucleotide divergence rates are atypical. The current lack of additional conifer plastome sequences precludes any discrimination between the gnetifer and gnepine hypotheses of seed plant relationships. However, both phylogenetic analyses and shared genome features identified here are consistent with either of the hypotheses that link gnetophytes with conifers, but are inconsistent with the anthophyte hypothesis.

Molecular identification of twelve medicinal plants of Huperziaceae

Twelve species of Huperziaceae have been identified for resolving the resource depletion of Huperziaceae plants. The chloroplast genes rbcL, rpL16, and psbA-trnH intergenic spacers were used in this study. Identification of different species of Huperziaceae is possible by analyzing the information obtained from rbcL, rpL16, and psbA-trnH intergenic spacer sequences. The result also can provide information with regard to the use of plant tissue culture and cuttage propagation to solve the problem of the resource depletion of these plants.

The complete nucleotide sequence of the cassava (Manihot esculenta) chloroplast genome and the evolution of atpF in Malpighiales: RNA editing and multiple losses of a group II intron

The complete sequence of the chloroplast genome of cassava (Manihot esculenta, Euphorbiaceae) has been determined. The genome is 161,453 bp in length and includes a pair of inverted repeats (IR) of 26,954 bp. The genome includes 128 genes; 96 are single copy and 16 are duplicated in the IR. There are four rRNA genes and 30 distinct tRNAs, seven of which are duplicated in the IR. The infA gene is absent; expansion of IRb has duplicated 62 amino acids at the 3' end of rps19 and a number of coding regions have large insertions or deletions, including insertions within the 23S rRNA gene. There are 17 intron-containing genes in cassava, 15 of which have a single intron while two (clpP, ycf3) have two introns. The usually conserved atpF group II intron is absent and this is the first report of its loss from land plant chloroplast genomes. The phylogenetic distribution of the atpF intron loss was determined by a PCR survey of 251 taxa representing 34 families of Malpighiales and 16 taxa from closely related rosids. The atpF intron is not only missing in cassava but also from closely related Euphorbiaceae and other Malpighiales, suggesting that there have been at least seven independent losses. In cassava and all other sequenced Malphigiales, atpF gene sequences showed a strong association between C-to-T substitutions at nucleotide position 92 and the loss of the intron, suggesting that recombination between an edited mRNA and the atpF gene may be a possible mechanism for the intron loss.

Construction of a tobacco master line to improve Rubisco engineering in chloroplasts

The inability to assemble Rubisco from any photosynthetic eukaryote within Escherichia coli has hampered structure-function studies of higher plant Rubisco. Precise genetic manipulation of the tobacco chloroplast genome (plastome) by homologous recombination has facilitated the successful production of transplastomic lines that have either mutated the Rubisco large subunit (L) gene, rbcL, or replaced it with foreign variants. Here the capacity of a new tobacco transplastomic line, (cm)trL, to augment future Rubisco engineering studies is demonstrated. Initially the rbcL was replaced with the selectable marker gene, aadA, and an artificial codon-modified (cm)rbcM gene that codes for the structurally novel Rubisco dimer (L(2), approximately 100 kDa) from Rhodosprillum rubrum. To obtain (cm)trL, the aadA was excised by transiently introducing a T-DNA encoding CRE recombinase biolistically. Selection using aadA enabled transplantation of mutated and wild-type tobacco Rubisco genes into the (cm)trL plastome with an efficiency that was 3- to 10-fold higher than comparable transformations into wild-type tobacco. Transformants producing the re-introduced form I tobacco Rubisco variants (hexadecamers comprising eight L and eight small subunits, approximately 520 kDa) were identified by non-denaturing PAGE with fully segregated homoplasmic lines (where no L(2) Rubisco was produced) obtained within 6-9 weeks after transformation which enabled their Rubisco kinetics to be quickly examined. Here the usefulness of (cm)trL in more readily examining the production, folding, and assembly capabilities of both mutated tobacco and foreign form I Rubisco subunits in tobacco plastids is discussed, and the feasibility of quickly assessing the kinetic properties of those that functionally assemble is demonstrated.

Superwobbling facilitates translation with reduced tRNA sets

Some bacterial and most organelle genomes do not encode the full set of 32 tRNA species required to read all codons according to Crick's wobble rules. 'Superwobble', in which a tRNA species with an unmodified U in the wobble position reads all four nucleotides in the third codon position, represents one possible mechanism for how a reduced tRNA set could still suffice. We have tested the superwobble hypothesis by producing knockout mutants for the pair of plastid glycine tRNA genes. Here we show that, whereas the tRNA gene with U in the wobble position is essential, the gene with G in this position is nonessential, demonstrating that the U-containing anticodon can indeed read all four glycine triplets. We also show that the price for superwobbling is a reduced translational efficiency, which explains why most organisms prefer pairs of isoaccepting tRNAs over the superwobbling mechanism.

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.

GenomeVx: simple web-based creation of editable circular chromosome maps

We describe GenomeVx, a web-based tool for making editable, publication-quality, maps of mitochondrial and chloroplast genomes and of large plasmids. These maps show the location of genes and chromosomal features as well as a position scale. The program takes as input either raw feature positions or GenBank records. In the latter case, features are automatically extracted and colored, an example of which is given. Output is in the Adobe Portable Document Format (PDF) and can be edited by programs such as Adobe Illustrator.

AVAILABILITY

GenomeVx is available at http://wolfe.gen.tcd.ie/GenomeVx

Evolution of early eukaryotic cells: genomes, proteomes, and compartments

Eukaryotes arose from an endosymbiotic association of an alpha-proteobacterium-like organism (the ancestor of mitochondria) with a host cell (lacking mitochondria or plastids). Plants arose by the addition of a cyanobacterium-like endosymbiont (the ancestor of plastids) to the two-member association. Each member of the association brought a unique internal environment and a unique genome. Analyses of recently acquired genomic sequences with newly developed algorithms have revealed (a) that the number of endosymbiont genes that remain in eukaryotic cells-principally in the nucleus-is surprisingly large, (b) that protein products of a large number of genes (or their descendents) that entered the association in the genome of the host are now directed to an organelle derived from an endosymbiont, and (c) that protein products of genes traceable to endosymbiont genomes are directed to the nucleo-cytoplasmic compartment. Consideration of these remarkable findings has led to the present suggestion that contemporary eukaryotic cells evolved through continual chance relocation and testing of genes as well as combinations of gene products and biochemical processes in each unique cell compartment derived from a member of the eukaryotic association. Most of these events occurred during about 300 million years, or so, before contemporary forms of eukaryotic cells appear in the fossil record; they continue today.

Molecular evolution of mitochondrial introns in the liverwort Marchantia polymorpha

We here describe in detail the characterization and molecular evolution of group II introns in the mitochondrial genome of the liverwort Marchantia polymorpha. We find that 18 introns of the 25 group II introns can be assigned by their similarities to six clusters, indicating an intra-genomic propagation of one ancestral intron each into the respective clusters in the liverwort mitochondrial genome. Interestingly, the intra-genomic propagation of some of these introns occurred only after the evolutionary separation of the bryophytes from the other clades of plants. Finally we report that the maturase-like sequences in the liverwort group II introns have further evolved by horizontal and independent transposition and substitution by analogous sequences from other fungal introns.

The catalytic properties of hybrid Rubisco comprising tobacco small and sunflower large subunits mirror the kinetically equivalent source Rubiscos and can support tobacco growth

Plastomic replacement of the tobacco (Nicotiana tabacum) Rubisco large subunit gene (rbcL) with that from sunflower (Helianthus annuus; rbcL(S)) produced tobacco(Rst) transformants that produced a hybrid Rubisco consisting of sunflower large and tobacco small subunits (L(s)S(t)). The tobacco(Rst) plants required CO(2) (0.5% v/v) supplementation to grow autotrophically from seed despite the substrate saturated carboxylation rate, K(m), for CO(2) and CO(2)/O(2) selectivity of the L(s)S(t) enzyme mirroring the kinetically equivalent tobacco and sunflower Rubiscos. Consequently, at the onset of exponential growth when the source strength and leaf L(s)S(t) content were sufficient, tobacco(Rst) plants grew to maturity without CO(2) supplementation. When grown under a high pCO(2), the tobacco(Rst) seedlings grew slower than tobacco and exhibited unique growth phenotypes: Juvenile plants formed clusters of 10 to 20 structurally simple oblanceolate leaves, developed multiple apical meristems, and the mature leaves displayed marginal curling and dimpling. Depending on developmental stage, the L(s)S(t) content in tobacco(Rst) leaves was 4- to 7-fold less than tobacco, and gas exchange coupled with chlorophyll fluorescence showed that at 2 mbar pCO(2) and growth illumination CO(2) assimilation in mature tobacco(Rst) leaves remained limited by Rubisco activity and its rate (approximately 11 micromol m(-2) s(-1)) was half that of tobacco controls. (35)S-methionine labeling showed the stability of assembled L(s)S(t) was similar to tobacco Rubisco and measurements of light transient CO(2) assimilation rates showed L(s)S(t) was adequately regulated by tobacco Rubisco activase. We conclude limitations to tobacco(Rst) growth primarily stem from reduced rbcL(S) mRNA levels and the translation and/or assembly of sunflower large with the tobacco small subunits that restricted L(s)S(t) synthesis.

Using plastid genome-scale data to resolve enigmatic relationships among basal angiosperms

Although great progress has been made in clarifying deep-level angiosperm relationships, several early nodes in the angiosperm branch of the Tree of Life have proved difficult to resolve. Perhaps the last great question remaining in basal angiosperm phylogeny involves the branching order among the five major clades of mesangiosperms (Ceratophyllum, Chloranthaceae, eudicots, magnoliids, and monocots). Previous analyses have found no consistent support for relationships among these clades. In an effort to resolve these relationships, we performed phylogenetic analyses of 61 plastid genes ( approximately 42,000 bp) for 45 taxa, including members of all major basal angiosperm lineages. We also report the complete plastid genome sequence of Ceratophyllum demersum. Parsimony analyses of combined and partitioned data sets varied in the placement of several taxa, particularly Ceratophyllum, whereas maximum-likelihood (ML) trees were more topologically stable. Total evidence ML analyses recovered a clade of Chloranthaceae + magnoliids as sister to a well supported clade of monocots + (Ceratophyllum + eudicots). ML bootstrap and Bayesian support values for these relationships were generally high, although approximately unbiased topology tests could not reject several alternative topologies. The extremely short branches separating these five lineages imply a rapid diversification estimated to have occurred between 143.8 +/- 4.8 and 140.3 +/- 4.8 Mya.

Phylogenetic and evolutionary implications of complete chloroplast genome sequences of four early-diverging angiosperms: Buxus (Buxaceae), Chloranthus (Chloranthaceae), Dioscorea (Dioscoreaceae), and Illicium (Schisandraceae)

We have determined the complete chloroplast genome sequences of four early-diverging lineages of angiosperms, Buxus (Buxaceae), Chloranthus (Chloranthaceae), Dioscorea (Dioscoreaceae), and Illicium (Schisandraceae), to examine the organization and evolution of plastid genomes and to estimate phylogenetic relationships among angiosperms. For the most part, the organization of these plastid genomes is quite similar to the ancestral angiosperm plastid genome with a few notable exceptions. Dioscorea has lost one protein-coding gene, rps16; this gene loss has also happened independently in four other land plant lineages, liverworts, conifers, Populus, and legumes. There has also been a small expansion of the inverted repeat (IR) in Dioscorea that has duplicated trnH-GUG. This event has also occurred multiple times in angiosperms, including in monocots, and in the two basal angiosperms Nuphar and Drimys. The Illicium chloroplast genome is unusual by having a 10 kb contraction of the IR. The four taxa sequenced represent key groups in resolving phylogenetic relationships among angiosperms. Illicium is one of the basal angiosperms in the Austrobaileyales, Chloranthus (Chloranthales) remains unplaced in angiosperm classifications, and Buxus and Dioscorea are early-diverging eudicots and monocots, respectively. We have used sequences for 61 shared protein-coding genes from these four genomes and combined them with sequences from 35 other genomes to estimate phylogenetic relationships using parsimony, likelihood, and Bayesian methods. There is strong congruence among the trees generated by the three methods, and most nodes have high levels of support. The results indicate that Amborella alone is sister to the remaining angiosperms; the Nymphaeales represent the next-diverging clade followed by Illicium; Chloranthus is sister to the magnoliids and together this group is sister to a large clade that includes eudicots and monocots; and Dioscorea represents an early-diverging lineage of monocots just internal to Acorus.

INTROGRESSIVE HYBRIDIZATION BETWEEN RHODODENDRON KIUSIANUM AND R. KAEMPFERI (ERICACEAE) IN KYUSHU, JAPAN BASED ON CHLOROPLAST DNA MARKERS

Wild evergreen azalea populations of Rhododendron kiusianum and R. kaempferi (Ericaceae) were analysed using a chloroplast DNA (cpDNA) PCR-RFLP marker that was used to detect introgressive hybridization in our previous study of the Kirishima Mts populations.

The populations of the intermediate region in the Unzen Mts, which show phenotypic variation, were demonstrated to result from interspecific hybridization between Rhododendron kiusianum and R. kaempferi, possessing cpDNA from either R. kiusianum (1030/420 bp) or R. kaempferi (950/420/80 bp).

Most individuals of Rhododendron kiusianum in the Kujyu Mts, the Aso Mts and the surrounding mountains exhibited the PCR-RFLP pattern of R. kaempferi. These results from the Kujyu Mts and the Aso Mts indicate that natural hybridization and cytoplasmic introgression from Rhododendron kaempferi to R. kiusianum have occurred in the relatively distant past. In the case of Mt Yufudake and Mt Haneyama, the Rhododendron kiusianum population retains the effects of natural hybridization with R. kaempferi in the cpDNA as well as in the variation in flower characteristics.

All individuals of Rhododendron kiusianum on Mt Onogaradake in the Takakuma Mts exhibit R. kiusianum cpDNA (1030/420 bp), in spite of variation in flower colour.

OrganellarGenomeDRAW (OGDRAW): a tool for the easy generation of high-quality custom graphical maps of plastid and mitochondrial genomes

Mitochondria and plastids are DNA-containing cell organelles whose genomes occur at high copy numbers per cell. Organellar genomes vary greatly in size ranging from approximately 15 kb for some animal mitochondrial genomes to more than 2 Mb for some plant mitochondrial genomes. The vast majority of organellar genomes map as circular molecules that are difficult to illustrate by available commercial or free software tools. Thus, published genome maps are extremely heterogeneous in design, often tediously drawn semi-manually and lack any consensus in display. Here, we present a new web-based tool, OrganellarGenomeDRAW (OGDRAW), which produces high-resolution custom graphical maps of DNA sequences as stored in standard GenBank format entries. GenBank data can be provided as either file uploads or accession numbers. The program is specially optimized for the display of chloroplast and mitochondrial genomes but can also be used to depict other circular DNA sequences. The design of the program core as a Perl module with an object-oriented interface allows easy integration into custom scripts.

Discoveries in Rubisco (Ribulose 1,5-bisphosphate carboxylase/oxygenase): a historical perspective

Historic discoveries and key observations related to Rubisco (Ribulose 1,5-bisphosphate carboxylase/oxygenase), from 1947 to 2006, are presented. Currently, around 200 papers describing Rubisco research are published each year and the literature contains more than 5000 manuscripts on the subject. While trying to ensure that all the major events over this period are recorded, this analysis will inevitably be incomplete and will reflect the areas of particular interest to the authors.

Photosynthetic H2 metabolism in Chlamydomonas reinhardtii (unicellular green algae)

Unicellular green algae have the ability to operate in two distinctly different environments (aerobic and anaerobic), and to photosynthetically generate molecular hydrogen (H2). A recently developed metabolic protocol in the green alga Chlamydomonas reinhardtii permitted separation of photosynthetic O2-evolution and carbon accumulation from anaerobic consumption of cellular metabolites and concomitant photosynthetic H2-evolution. The H2 evolution process was induced upon sulfate nutrient deprivation of the cells, which reversibly inhibits photosystem-II and O2-evolution in their chloroplast. In the absence of O2, and in order to generate ATP, green algae resorted to anaerobic photosynthetic metabolism, evolved H2 in the light and consumed endogenous substrate. This study summarizes recent advances on green algal hydrogen metabolism and discusses avenues of research for the further development of this method. Included is the mechanism of a substantial tenfold starch accumulation in the cells, observed promptly upon S-deprivation, and the regulated starch and protein catabolism during the subsequent H2-evolution. Also discussed is the function of a chloroplast envelope-localized sulfate permease, and the photosynthesis-respiration relationship in green algae as potential tools by which to stabilize and enhance H2 metabolism. In addition to potential practical applications of H2, approaches discussed in this work are beginning to address the biochemistry of anaerobic H2 photoproduction, its genes, proteins, regulation, and communication with other metabolic pathways in microalgae. Photosynthetic H2 production by green algae may hold the promise of generating a renewable fuel from nature's most plentiful resources, sunlight and water. The process potentially concerns global warming and the question of energy supply and demand.

Stable chloroplast transformation in cabbage (Brassica oleracea L. var. capitata L.) by particle bombardment

The objectives of this research were first to isolate plastid gene sequences from cabbage (Brassica oleracea L. var. capitata L.), and to establish the chloroplast transformation technology of Brassica. A universal transformation vector (pASCC201) for Brassica chloroplast was constructed with trnV-rrn16S (left) and trnI-trnA-rrn23S (right) of the IR(_A) region as a recombination site for the transformed gene. In transforming plasmid pASCC201, a chimeric aadA gene was cloned between the rrn16S and rrn23S plastid gene borders. Expression of aadA confers resistance to spectinomycin and streptomycin antibiotics. The uidA gene was also inserted into the pASCC201 and transferred into the leaf cells of cabbage via particle gun mediated transformation. Regenerated plantlets were selected by 200 mg/l spectinomycin and streptomycin. After antibiotic selection, the regeneration percentage of the two cabbage cultivars was about 2.7-3.3%. The results of PCR testing and Southern blot analysis confirmed that the uidA and aadA genes were present in the chloroplast genome via homologously recombined. Northern blot hybridizations, immunoblotting and GUS histochemical assays indicated that the uidA gene were stable integrated into the chloroplast genome. Foreign protein was accumulated at 3.2-5.2% of the total soluble protein in transgenic mature leaves. These results suggest that the expression of a variety of foreign genes in the chloroplast genome will be a powerful tool for use in future studies.

Translation of psbC mRNAs starts from the downstream GUG, not the upstream AUG, and requires the extended Shine-Dalgarno sequence in tobacco chloroplasts

The plastid gene psbC encodes the CP43 subunit of PSII. Most psbC mRNAs of many organisms possess two possible initiation codons, AUG and GUG, and their coding regions are generally annotated from the upstream AUG. Using a chloroplast in vitro translation system, we show here that translation of the tobacco plastid psbC mRNA initiates from the GUG. This mRNA possesses a long Shine-Dalgarno (SD)-like sequence, GAGGAGGU, nine nucleotides upstream of the GUG. Point mutations in this sequence abolished translation, suggesting that a strong interaction between this extended SD-like sequence and the 3' end of 16S rRNA facilitates translation initiation from the GUG.

Plastid tRNA genes trnC-GCA and trnN-GUU are essential for plant cell development

Higher plant chloroplast genomes code for a conserved set of 30 tRNAs. This set is believed to be sufficient to support translation, although import of cytosolic tRNA has been proposed to provide additional tRNA species to the chloroplast. Previous knock-outs of tRNA genes, or the pronounced reduction of the level of selected tRNAs, has not led to severe phenotypes. We deleted the two tRNA genes trnN-GUU and trnC-GCA independently from the plastid chromosome of tobacco. No homoplastomic tissue of either DeltatrnN or DeltatrnC plants could be isolated. Both mutants exhibit occasional loss of leaf sectors, and mutant plastid chromosomes are rapidly lost upon relief of selective pressure. This suggests that the knock-out of both trn genes is lethal, and that both tRNA species are required for cell survival. Surprisingly, the impact on chloroplast and cell development differs pronouncedly between the two mutants. Heteroplastomic DeltatrnC and DeltatrnN tissue exhibit different aberrations of the internal membrane systems and, more importantly, heteroplastomic DeltatrnN plants are variegated. Accumulation of tRNA-N and plastid-encoded proteins is reduced in white sectors of DeltatrnN plants, and differentiation of palisade cells is abolished. Our data demonstrate that plastid tRNAs are essential, i.e. not complemented by cytosolic tRNA, and have a differential impact on chloroplast and plant cell development.

Complete DNA sequences of the plastid genomes of two parasitic flowering plant species, Cuscuta reflexa and Cuscuta gronovii

BACKGROUND

The holoparasitic plant genus Cuscuta comprises species with photosynthetic capacity and functional chloroplasts as well as achlorophyllous and intermediate forms with restricted photosynthetic activity and degenerated chloroplasts. Previous data indicated significant differences with respect to the plastid genome coding capacity in different Cuscuta species that could correlate with their photosynthetic activity. In order to shed light on the molecular changes accompanying the parasitic lifestyle, we sequenced the plastid chromosomes of the two species Cuscuta reflexa and Cuscuta gronovii. Both species are capable of performing photosynthesis, albeit with varying efficiencies. Together with the plastid genome of Epifagus virginiana, an achlorophyllous parasitic plant whose plastid genome has been sequenced, these species represent a series of progression towards total dependency on the host plant, ranging from reduced levels of photosynthesis in C. reflexa to a restricted photosynthetic activity and degenerated chloroplasts in C. gronovii to an achlorophyllous state in E. virginiana.

RESULTS

The newly sequenced plastid genomes of C. reflexa and C. gronovii reveal that the chromosome structures are generally very similar to that of non-parasitic plants, although a number of species-specific insertions, deletions (indels) and sequence inversions were identified. However, we observed a gradual adaptation of the plastid genome to the different degrees of parasitism. The changes are particularly evident in C. gronovii and include (a) the parallel losses of genes for the subunits of the plastid-encoded RNA polymerase and the corresponding promoters from the plastid genome, (b) the first documented loss of the gene for a putative splicing factor, MatK, from the plastid genome and (c) a significant reduction of RNA editing.

CONCLUSION

Overall, the comparative genomic analysis of plastid DNA from parasitic plants indicates a bias towards a simplification of the plastid gene expression machinery as a consequence of an increasing dependency on the host plant. A tentative assignment of the successive events in the adaptation of the plastid genomes to parasitism can be inferred from the current data set. This includes (1) a loss of non-coding regions in photosynthetic Cuscuta species that has resulted in a condensation of the plastid genome, (2) the simplification of plastid gene expression in species with largely impaired photosynthetic capacity and (3) the deletion of a significant part of the genetic information, including the information for the photosynthetic apparatus, in non-photosynthetic parasitic plants.

Complete chloroplast genome sequences of Hordeum vulgare, Sorghum bicolor and Agrostis stolonifera, and comparative analyses with other grass genomes

Comparisons of complete chloroplast genome sequences of Hordeum vulgare, Sorghum bicolor and Agrostis stolonifera to six published grass chloroplast genomes reveal that gene content and order are similar but two microstructural changes have occurred. First, the expansion of the IR at the SSC/IRa boundary that duplicates a portion of the 5' end of ndhH is restricted to the three genera of the subfamily Pooideae (Agrostis, Hordeum and Triticum). Second, a 6 bp deletion in ndhK is shared by Agrostis, Hordeum, Oryza and Triticum, and this event supports the sister relationship between the subfamilies Erhartoideae and Pooideae. Repeat analysis identified 19-37 direct and inverted repeats 30 bp or longer with a sequence identity of at least 90%. Seventeen of the 26 shared repeats are found in all the grass chloroplast genomes examined and are located in the same genes or intergenic spacer (IGS) regions. Examination of simple sequence repeats (SSRs) identified 16-21 potential polymorphic SSRs. Five IGS regions have 100% sequence identity among Zea mays, Saccharum officinarum and Sorghum bicolor, whereas no spacer regions were identical among Oryza sativa, Triticum aestivum, H. vulgare and A. stolonifera despite their close phylogenetic relationship. Alignment of EST sequences and DNA coding sequences identified six C-U conversions in both Sorghum bicolor and H. vulgare but only one in A. stolonifera. Phylogenetic trees based on DNA sequences of 61 protein-coding genes of 38 taxa using both maximum parsimony and likelihood methods provide moderate support for a sister relationship between the subfamilies Erhartoideae and Pooideae.

Role of the low-molecular-weight subunits PetL, PetG, and PetN in assembly, stability, and dimerization of the cytochrome b6f complex in tobacco

The cytochrome b(6)f (Cyt b(6)f) complex in flowering plants contains nine conserved subunits, of which three, PetG, PetL, and PetN, are bitopic plastid-encoded low-molecular-weight proteins of largely unknown function. Homoplastomic knockout lines of the three genes have been generated in tobacco (Nicotiana tabacum 'Petit Havana') to analyze and compare their roles in assembly and stability of the complex. Deletion of petG or petN caused a bleached phenotype and loss of photosynthetic electron transport and photoautotrophy. Levels of all subunits that constitute the Cyt b(6)f complex were faintly detectable, indicating that both proteins are essential for the stability of the membrane complex. In contrast, DeltapetL plants accumulate about 50% of other Cyt b(6)f subunits, appear green, and grow photoautotrophically. However, DeltapetL plants show increased light sensitivity as compared to wild type. Assembly studies revealed that PetL is primarily required for proper conformation of the Rieske protein, leading to stability and formation of dimeric Cyt b(6)f complexes. Unlike wild type, phosphorylation levels of the outer antenna of photosystem II (PSII) are significantly decreased under state II conditions, although the plastoquinone pool is largely reduced in DeltapetL, as revealed by measurements of PSI and PSII redox states. This confirms the sensory role of the Cyt b(6)f complex in activation of the corresponding kinase. The reduced light-harvesting complex II phosphorylation did not affect state transition and association of light-harvesting complex II to PSI under state II conditions. Ferredoxin-dependent plastoquinone reduction, which functions in cyclic electron transport around PSI in vivo, was not impaired in DeltapetL.

Disruption of essential plastid gene expression caused by T7 RNA polymerase-mediated transcription of plastid transgenes during early seedling development

Transcription of plastid transgenes by plastid-targeted T7 RNA polymerase (ptT7RNAP) during early seedling development in tobacco was associated with a pale-green leaf phenotype, depletion of plastid rRNAs and arrest of shoot development. Extensive analysis of mutant seedlings at the transcript level using DNA microarrays and RNA gel blotting revealed severe disruption of plastid rRNA accumulation at 4-days post-germination and reduced transcript accumulation for the essential gene clpP. Several nuclear genes encoding plastid proteins were differentially regulated in mutant seedlings over time. Ef-Tu was upregulated at 4-days post-germination and then subsequently downregulated, while RbcS was already downregulated at this early time point. The downregulation of nuclear genes encoding plastid proteins suggests disruption of plastid-to-nucleus signalling. In contrast, transcripts of three plastid genes showed increased accumulation in mutant seedlings. Transcripts of ndhC and ndhK accumulated at high levels possibly due to T7RNAP-mediated enhancement of transcription, while ptT7RNAP-mediated transcription through the phage T7 Tphi terminator into the adjacent plastome increased the level of accD transcripts. The leakiness of the Tphi terminator has implications for the use of T7RNAP-based expression systems in plastid biotechnology.

Complete chloroplast genome sequences of Hordeum vulgare, Sorghum bicolor and Agrostis stolonifera, and comparative analyses with other grass genomes

Comparisons of complete chloroplast genome sequences of Hordeum vulgare, Sorghum bicolor and Agrostis stolonifera to six published grass chloroplast genomes reveal that gene content and order are similar but two microstructural changes have occurred. First, the expansion of the IR at the SSC/IRa boundary that duplicates a portion of the 5' end of ndhH is restricted to the three genera of the subfamily Pooideae (Agrostis, Hordeum and Triticum). Second, a 6 bp deletion in ndhK is shared by Agrostis, Hordeum, Oryza and Triticum, and this event supports the sister relationship between the subfamilies Erhartoideae and Pooideae. Repeat analysis identified 19-37 direct and inverted repeats 30 bp or longer with a sequence identity of at least 90%. Seventeen of the 26 shared repeats are found in all the grass chloroplast genomes examined and are located in the same genes or intergenic spacer (IGS) regions. Examination of simple sequence repeats (SSRs) identified 16-21 potential polymorphic SSRs. Five IGS regions have 100% sequence identity among Zea mays, Saccharum officinarum and Sorghum bicolor, whereas no spacer regions were identical among Oryza sativa, Triticum aestivum, H. vulgare and A. stolonifera despite their close phylogenetic relationship. Alignment of EST sequences and DNA coding sequences identified six C-U conversions in both Sorghum bicolor and H. vulgare but only one in A. stolonifera. Phylogenetic trees based on DNA sequences of 61 protein-coding genes of 38 taxa using both maximum parsimony and likelihood methods provide moderate support for a sister relationship between the subfamilies Erhartoideae and Pooideae.

Sequencing cucumber (Cucumis sativus L.) chloroplast genomes identifies differences between chilling-tolerant and -susceptible cucumber lines

Chilling injury in cucumber (Cucumis sativus L.) is conditioned by maternal factors, and the sequencing of its chloroplast genome could lead to the identification of economically important candidate genes. Complete sequencing of cucumber chloroplast (cp)DNA was facilitated by the development of 414 consensus chloroplast sequencing primers (CCSPs) from conserved cpDNA sequences of Arabidopsis (Arabidopsis thaliana L.), spinach (Spinacia oleracea L.), and tobacco (Nicotiana tabacum L.) cpDNAs, using degenerative primer technologies. Genomic sequence analysis led to the construction of 301 CCSPs and 72 cucumber chloroplast-specific sequencing primers (CSSPs), which were used for the complete sequencing of cpDNA of Gy14 (155 525 bp) and 'Chipper' (155 524 bp) cucumber lines, which are, respectively, susceptible and tolerant to chilling injury (4 degrees C for 5.5 h) in the first leaf stage. Comparative cpDNA sequence analyses revealed that 1 sequence span (located between genes trnK and rps16) and 2 nucleotides (located in genes atpB and ycf1) differed between chilling-susceptible and -tolerant lines. These sequence differences correspond to previously reported maternally inherited differences in chilling response between reciprocal F1 progeny derived from these lines. Sequence differences at these 3 cpDNA sites were also detected in a genetically diverse array of cucumber germplasm with different chilling responses. These and previously reported results suggest that 1 or several of these sequences could be responsible for the observed response to chilling injury in cucumber. The comprehensive sequencing of cpDNA of cucumber by CCSPs and CSSPs indicates that these primers have immediate applications in the analysis of cpDNAs from other dicotyledonous species and the investigation of evolutionary relationships.

The complete structure of the cucumber (Cucumis sativus L.) chloroplast genome: its composition and comparative analysis

The complete nucleotide sequence of the cucumber (C. sativus L. var. Borszczagowski) chloroplast genome has been determined. The genome is composed of 155,293 bp containing a pair of inverted repeats of 25,191 bp, which are separated by two single-copy regions, a small 18,222-bp one and a large 86,688-bp one. The chloroplast genome of cucumber contains 130 known genes, including 89 protein-coding genes, 8 ribosomal RNA genes (4 rRNA species), and 37 tRNA genes (30 tRNA species), with 18 of them located in the inverted repeat region. Of these genes, 16 contain one intron, and two genes and one ycf contain 2 introns. Twenty-one small inversions that form stem-loop structures, ranging from 18 to 49 bp, have been identified. Eight of them show similarity to those of other species, while eight seem to be cucumber specific. Detailed comparisons of ycf2 and ycf15, and the overall structure to other chloroplast genomes were performed.

Chloroplast DNA phylogeography of Pedicularis ser. Gloriosae (Orobanchaceae) in Japan

Phylogeographic analyses using chloroplast DNA (cpDNA) variation were performed for Pedicularis ser. Gloriosae (Orobanchaceae). Eighty-one plants of 18 populations of 6 species (P. gloriosa, P. iwatensis, P. nipponica, P. ochiaiana, P. sceptrum-carolinum and P. grandiflora) were analyzed. Fifteen distinct haplotypes were identified based on six cpDNA regions: the intergenic spacer between the trnT and trnL 3'exon, trnL 3'exon-trnF, atpB-rbcL, accD-psaI, the rpl16 intron and the trnK region (including the matK gene). Via phylogenetic analyses of the haplotypes, two continental species, P. sceptrum-carolinum and P. grandiflora, were placed at the most ancestral position in the trees. The former species is widely distributed in the Eurasian continent, and the latter is distributed in Far East Asia. Two robust major cpDNA clades (clades I and II) were revealed in the Japanese archipelago, although the statistical values of monophyly of these clades were weak. Clade I included the haplotypes (A-1, A-2, B-1, B-2 and J) of three species (P. gloriosa, P. iwatensis and P. ochiaiana), and Clade II included seven haplotypes (C-D, E-1, E-2 and F-H) of P. nipponica. These results suggest that this series originated on the Eurasian continent and that subsequently populations at the eastern edge of the continent differentiated into the two Japanese lineages.

Targeted inactivation of the tobacco plastome origins of replication A and B

According to the Kolodner and Tewari model [Kolodner, R.D. and Tewari, K.K. (1975) Nature, 256, 708.], plastid DNA replication involves displacement-loop and rolling-circle modes of replication, which are initiated on a pair of origins of replication (ori). In accordance with the model, such a pair of oris -oriA and oriB- was described in Nicotiana tabacum [Kunnimalaiyaan, M. and Nielsen B.L. (1997b) Nucl. Acids Res. 25, 3681.]. However, as reported previously, both copies of oriA can be deleted without abolishing replication. Deletion of both oriBs was not found [Mühlbauer, S.K. et al. (2002) Plant J. 32, 175.]. Here we describe new ori inactivation lines, in which one oriB is deleted and the other copy is strongly mutated. In addition, lines oriA and oriB were deleted from the same inverted repeat. In contrast to the expectations of the model, neither oriA nor oriB is essential. Some of the deletions led to reduced growth of plants and reduced plastid DNA copy number in later stages of leaf development. The gross structure of plastid DNA was unchanged; however, the location of the ends of branched plastid DNA complexes was different in the inactivation mutants. Taken together, the results indicate that there are additional mechanisms of plastid DNA replication and/or additional origins of replication. These mechanisms seem to be different from those found in eubacteria, which, according to the endosymbiont theory, are the progenitors of plastids.

Plastid marker gene excision by the phiC31 phage site-specific recombinase

Marker genes are essential for selective amplification of rare transformed plastid genome copies to obtain genetically stable transplastomic plants. However, the marker gene becomes dispensable when homoplastomic plants are obtained. Here we report excision of plastid marker genes by the phiC31 phage site-specific integrase (Int) that mediates recombination between bacterial (attB) and phage (attP) attachment sites. We tested marker gene excision in a two-step process. First we transformed the tobacco plastid genome with the pCK2 vector in which the spectinomycin resistance (aadA) marker gene is flanked with suitably oriented attB and attP sites. The transformed plastid genomes were stable in the absence of Int. We then transformed the nucleus with a gene encoding a plastid-targeted Int that led to efficient marker gene excision. The aadA marker free Nt-pCK2-Int plants were resistant to phosphinothricin herbicides since the pCK2 plastid vector also carried a bar herbicide resistance gene that, due to the choice of its promoter, causes a yellowish-golden (aurea) phenotype. Int-mediated marker excision reported here is an alternative to the currently used CRE/loxP plastid marker excision system and expands the repertoire of the tools available for the manipulation of the plastid genome.

Disruption of the psbA gene by the copy correction mechanism reveals that the expression of plastid-encoded genes is regulated by photosynthesis activity

The functional analysis of genes encoded by the chloroplast genome of tobacco by reverse genetics is routine. Nevertheless, for a small number of genes their deletion generates heteroplasmic genotypes, complicating their analysis. There is thus the need for additional strategies to develop deletion mutants for these genes. We have developed a homologous copy correction-based strategy for deleting/mutating genes encoded on the chloroplast genome. This system was used to produce psbA knockouts. The resulting plants are homoplasmic and lack photosystem II (PSII) activity. Further, the deletion mutants exhibit a distinct phenotype; young leaves are green, whereas older leaves are bleached, irrespective of light conditions. This suggests that senescence is promoted by the absence of psbA. Analysis of the transcript levels indicates that NEP (nuclear-encoded plastid RNA polymerase)-dependent plastid genes are up regulated in the psbA deletion mutants, whereas the bleached leaves retain plastid-encoded plastid RNA polymerase activity. Hence, the expression of NEP-dependent plastid genes may be regulated by photosynthesis, either directly or indirectly.

Exceptional transmission of plastids and mitochondria from the transplastomic pollen parent and its impact on transgene containment

Plastids in Nicotiana tabacum are normally transmitted to the progeny by the maternal parent only. However, low-frequency paternal plastid transmission has been reported in crosses involving parents with an alien cytoplasm. Our objective was to determine whether paternal plastids are transmitted in crosses between parents with the normal cytoplasm. The transplastomic father lines carried a spectinomycin resistance (aadA) transgene incorporated in the plastid genome. The mother lines in the crosses were either (i) alloplasmic, with the Nicotiana undulata cytoplasm that confers cytoplasmic male sterility (CMS92) or (ii) normal, with the fertile N. tabacum cytoplasm. Here we report that plastids from the transplastomic father were transmitted in both cases at low (10(-4)-10(-5)) frequencies; therefore, rare paternal pollen transmission is not simply due to breakdown of normal controls caused by the alien cytoplasm. Furthermore, we have found that the entire plastid genome was transmitted by pollen rather than small plastid genome (ptDNA) fragments. Interestingly, the plants, which inherited paternal plastids, also carried paternal mitochondrial DNA, indicating cotransmission of plastids and mitochondria in the same pollen. The detection of rare paternal plastid transmission described here was facilitated by direct selection for the transplastomic spectinomycin resistance marker in tissue culture; therefore, recovery of rare paternal plastids in the germline is less likely to occur under field conditions.

The chloroplast genome from a lycophyte (microphyllophyte), Selaginella uncinata, has a unique inversion, transpositions and many gene losses

We determined the complete nucleotide sequence of the chloroplast genome of Selaginella uncinata, a lycophyte belonging to the basal lineage of the vascular plants. The circular double-stranded DNA is 144,170 bp, with an inverted repeat of 25,578 bp separated by a large single copy region (LSC) of 77,706 bp and a small single copy region (SSC) of 40,886 bp. We assigned 81 protein-coding genes including four pseudogenes, four rRNA genes and only 12 tRNA genes. Four genes, rps15, rps16, rpl32 and ycf10, found in most chloroplast genomes in land plants were not present in S. uncinata. While gene order and arrangement of the chloroplast genome of another lycophyte, Hupertzia lucidula, are almost the same as those of bryophytes, those of S. uncinata differ considerably from the typical structure of bryophytes with respect to the presence of a unique 20 kb inversion within the LSC, transposition of two segments from the LSC to the SSC and many gene losses. Thus, the organization of the S. uncinata chloroplast genome provides a new insight into the evolution of lycophytes, which were separated from euphyllophytes approximately 400 million years ago.

The complete nucleotide sequence of the coffee (Coffea arabica L.) chloroplast genome: organization and implications for biotechnology and phylogenetic relationships amongst angiosperms

The chloroplast genome sequence of Coffea arabica L., the first sequenced member of the fourth largest family of angiosperms, Rubiaceae, is reported. The genome is 155 189 bp in length, including a pair of inverted repeats of 25,943 bp. Of the 130 genes present, 112 are distinct and 18 are duplicated in the inverted repeat. The coding region comprises 79 protein genes, 29 transfer RNA genes, four ribosomal RNA genes and 18 genes containing introns (three with three exons). Repeat analysis revealed five direct and three inverted repeats of 30 bp or longer with a sequence identity of 90% or more. Comparisons of the coffee chloroplast genome with sequenced genomes of the closely related family Solanaceae indicated that coffee has a portion of rps19 duplicated in the inverted repeat and an intact copy of infA. Furthermore, whole-genome comparisons identified large indels (> 500 bp) in several intergenic spacer regions and introns in the Solanaceae, including trnE (UUC)-trnT (GGU) spacer, ycf4-cemA spacer, trnI (GAU) intron and rrn5-trnR (ACG) spacer. Phylogenetic analyses based on the DNA sequences of 61 protein-coding genes for 35 taxa, performed using both maximum parsimony and maximum likelihood methods, strongly supported the monophyly of several major clades of angiosperms, including monocots, eudicots, rosids, asterids, eurosids II, and euasterids I and II. Coffea (Rubiaceae, Gentianales) is only the second order sampled from the euasterid I clade. The availability of the complete chloroplast genome of coffee provides regulatory and intergenic spacer sequences for utilization in chloroplast genetic engineering to improve this important crop.