Chloroplast genome sequence of the moss Tortula ruralis: gene content, polymorphism, and structural arrangement relative to other green plant chloroplast genomes
© Oliver et al; licensee BioMed Central Ltd. 2010
Received: 27 August 2009
Accepted: 27 February 2010
Published: 27 February 2010
Tortula ruralis, a widely distributed species in the moss family Pottiaceae, is increasingly used as a model organism for the study of desiccation tolerance and mechanisms of cellular repair. In this paper, we present the chloroplast genome sequence of T. ruralis, only the second published chloroplast genome for a moss, and the first for a vegetatively desiccation-tolerant plant.
The Tortula chloroplast genome is ~123,500 bp, and differs in a number of ways from that of Physcomitrella patens, the first published moss chloroplast genome. For example, Tortula lacks the ~71 kb inversion found in the large single copy region of the Physcomitrella genome and other members of the Funariales. Also, the Tortula chloroplast genome lacks petN, a gene found in all known land plant plastid genomes. In addition, an unusual case of nucleotide polymorphism was discovered.
Although the chloroplast genome of Tortula ruralis differs from that of the only other sequenced moss, Physcomitrella patens, we have yet to determine the biological significance of the differences. The polymorphisms we have uncovered in the sequencing of the genome offer a rare possibility (for mosses) of the generation of DNA markers for fine-level phylogenetic studies, or to investigate individual variation within populations.
Tortula ruralis (Hedw.) Gaertn., also known as Syntrichia ruralis (Hedw.) F. Weber & D. Mohr (Pottiaceae) is a moss with a cosmopolitan distribution in relatively dry habitats. In North America the species is widespread in northern latitudes but is more common in the Western U.S., south into Mexico . Tortula ruralis has received considerable attention over the last forty years as a model for the study of vegetative desiccation tolerance, i.e., the ability to equilibrate to the water potential of dry air and survive, regaining growth and development upon rehydration. Tortula ruralis offers much as an experimental model for the study of environmental impacts on plants: it grows easily in culture, has a limited number of cell types, and, because of its morphology, experimental treatments act directly at the cellular level [2, 3]. It is the latter property that also makes it an ideal choice for an indicator species in air pollution studies [4, 5].
Tortula ruralis is among the most desiccation-tolerant of land plants and it can recover from desiccation even after at least three years in the dried state [6, 7]. Physcomitrella patens is relatively tolerant of dehydration but cannot tolerate the levels of drying that T. ruralis can survive . It is well established that the chloroplast plays a central role in the recovery of vegetative plants cells from desiccation  and it is possible that differences between the chloroplast genomes of T. ruralis and P. patens may relate to this fundamental difference between the two mosses.
The rapid recovery of photosynthesis is critical in order to recover and re-establish growth when water is available, thus maximizing the time available to the moss for carbon fixation and productivity . Following slow drying to -100 Mpa, photosystem II (PSII) activity in T. ruralis recovers within minutes after rewetting , with normal rates of carbon fixation returning within an hour . Photosynthesis is essential for the production of the energy required for repair and protein synthesis following the desiccation event. Obviously, the integrity and metabolic capacity of the chloroplast is central to the speed of recovery of photosynthesis. It is clear from electron microscopic observation of freeze-fracture preparations that chloroplast membranes, both the envelope and thylakoid membranes, in T. ruralis are unaltered by desiccation , which supports the idea that desiccation does not damage the photosynthetic apparatus. Such protection of chloroplast structure has also been demonstrated for gametophytes of Polytrichum formosum, which also appear to be unaltered by the imposition of desiccation and the rigors of rehydration . Thus it is clear that the chloroplast holds a central role in the response of T. ruralis to desiccation and rehydration and it is important to study the nature of its genome in this plant, the first vegetatively desiccation-tolerant plant to have its chloroplast genome sequenced. The genome sequence of T. ruralis, and its comparison to other chloroplast genomes, is critical if we want to understand how the interplay between the nuclear and chloroplast genomes plays a role in desiccation tolerance.
In addition to the relevance of the T. ruralis chloroplast genome to the important trait of desiccation tolerance, the genome sequence has considerable relevance to our current understanding of evolutionary history of the land plants. Current evidence suggests that mosses are the sister group of hornworts plus tracheophytes, diverging at least 450 million years ago [13, 14]. As an early diverging lineage, mosses hold a place in the phylogeny of land plants that is important for comparative purposes to seed plants , although comparisons are currently hampered because only one published chloroplast genome is available for mosses, while hundreds are available for its sister group. Several chloroplast genome sequences will be required to estimate the ancestral genome sequence for mosses, which will in turn allow comparisons with tracheophyte genomes.
The interest in Tortula ruralis as a model desiccation-tolerant organism has increased as our need to understand how plants survive dehydration stress grows and the global impact of climate change becomes more critical. This impetus and the need for increasing sampling within the mosses for phylogenetic comparative purposes led to the choice of T. ruralis for chloroplast genome sequencing. The assembly and annotation of the T. ruralis chloroplast genome sequence is presented here, only the second chloroplast genome sequenced for a moss and the first for a desiccation tolerant plant. The first chloroplast genome for a moss, Physcomitrella patens, was completed in 2003 . Because P. patens was found to have a major rearrangement in the chloroplast compared to what PCR-based methods show for most other moss lineages , the T. ruralis chloroplast genome will serve as an important point of comparison and will assist in ongoing efforts to utilize whole-genomic sequences and structural characters in a comparative phylogenetic framework [13, 18].
Whole chloroplast genome description
List of all sequence polymorphisms located on the Tortula ruralis chloroplast genome sequence.
GTA or GTG - synonymous
TTA or TTC - nonsynonymous
GTA or TTA - nonsynonymous
AGT or ATT - nonsynonymous
TAA or TAG - synonymous
TTG or TTT - nonsynonymous
CGA or CGG - synonymous
TGT or TTT - nonsynonymous
CGC or CGT - synonymous
CCA or TCA nonsynonymous
rps12 exon2/rps7 IGR
GTA or GTC - synonymous
Whole chloroplast genome description
Compared to other published chloroplast genomes, the chloroplast genome of Tortula ruralis is most similar to that of Physcomitrella patens. However, T. ruralis lacks the large ~71 kb inversion in the LSC region of the genome that is found in P. patens. This large inversion is found only in the moss order Funariales, whereas other moss lineages have a more plesiomorphic gene order similar to the liverwort Marchantia polymorpha [16, 17].
It is interesting that the Tortula ruralis chloroplast genome lacks the petN gene because the protein encoded by this gene plays an important role in photosynthesis electron transport. The PetN protein is a subunit of the cytochrome b6f (Cyt b6f) complex of cyanobacteria and plants. The Cytb6f complex is a plastoquinol-plastocyanin oxidoreductase within thylakoid membranes that functions as the photosynthetic redox control of energy distribution between the two photosystems, PSII and PSI, and gene expression . With such a critical role in chloroplast function it is clear that T. ruralis cannot be devoid of the PetN polypeptide and it is very likely, given the conservation of gene order in plant chloroplasts, that the petN gene does not reside in the region of the genome for which we do not yet have a sequence (e.g., within ycf2). What remains is the probability that in T. ruralis the petN gene has moved into the nuclear genome or there is a nuclear-encoded gene product that serves the same function as the PetN protein. This hypothesis remains to be tested but it offers an exciting speculation that the translocation of the petN gene from the chloroplast to nuclear genome may be related to the stability of the chloroplast in desiccation-tolerant mosses. Obviously, a much more detailed survey of moss chloroplast genomes needs to be accomplished and further experimental evidence needs to be collected before such speculation can be tested.
Few studies have detected nucleotide polymorphism for plastid DNA within plant populations . Our results suggest a considerable level of population polymorphism in the Tortula ruralis chloroplast DNA. Mosses have monoplastidic cell division [21, 22], thus it is unlikely that the polymorphism occurs within individual plants. Because the sequenced material was grown from wild-collected spores from multiple parents, the polymorphism is more likely to be due to genotypic races or otherwise cryptic lineages within the source population. One alternative source of polymorphism that is difficult to discount completely is that the variation is caused by inclusion of nuclear DNA which can occasionally contain fragments of chloroplast DNA . We consider this unlikely because our source of plastid DNA was from isolated chloroplasts, so nuclear DNA levels would be very low and the chances of capturing rare inserts would be negligible. Also, the 29 polymorphisms were distributed across the genome rather than clustered in one region which is what would be expected if one or a few nuclear regions had been captured. Many recently sequenced plastid genomes have used shotgun or PCR-based approaches, and completely eliminating nuclear DNA is impossible.
The polymorphisms that give rise to non-synonymous changes in amino acid at a defined site in a protein coding region are possible targets for post transcriptional editing should the non-synonymous change result in a loss or detrimental change in protein function. For example, the polymorphic position 744 in the matK gene (A or C) could result in either leucine (UUA) or phenylalanine (UUC) in the resultant peptide. Given the importance of this group II intron maturase gene to chloroplast function  this change either does not result in a significant change in the normal activity of the polypeptide or it would be corrected by the RNA editing activity within the chloroplast prior to translation. Testing this hypothesis, and estimating the extent of RNA editing will require sequencing of cDNA from Tortula chloroplast genes.
The shotgun sequencing method used for assembling this genome makes a determination of the number of haplotypes in the sequenced sample impossible. However, there are some indications that only two haplotypes may be present: (1) all variable sites were bimodal (e.g., either an A or G); and (2) in 5 instances, multiple polymorphic sites were located within a single DNA fragment and in all 5 cases only two versions of the sequence were present. Interestingly, one third of the variable sites occurred in coding regions of the chloroplast, including two sites in rbcL and three in ndhF. Polymorphic sites that can be easily sequenced from chloroplast DNA have potential for future use in population-level studies in T. ruralis. Future studies that include cloning and sequencing chloroplasts genomic regions from single individuals, gathered from across the range of this species, would be very productive. It would also be of interest to determine if the polymorphisms within the coding regions of the chloroplast genes generate functional allelic variation in the gene products involved.
Among the known chloroplast genomes of land plants, the trnPGGG gene shows an unusual evolutionary pattern. In Physcomitrella patens and Tortula ruralis the sequence data indicate that this gene has become a pseudogene in these lineages (Figure 2); further substantiated by the low Cove scores for the tRNA structure prediction from tRNAscan-SE . In the liverwort Marchantia polymorpha, the anticodon is intact but there are stem mismatches and a deletion of 5 bases, indicating that this gene may also be pseudogenized in M. polymorpha. However, in the hornwort Anthoceros formosae and in vascular plants trnPGGG appears to be functional. It is of course possible that the change in the anticodon sequence seen in T. ruralis and P. patens are modified post-transcriptionally to allow non-wobble interactions and thus rendering the tRNA functional. However, if these alterations in the trnPGGG seen in T. ruralis and P. patens do indicate a loss of function then given the current understanding of land plant phylogeny, this would imply at least two independent deactivation events for the trnPGGG gene in the chloroplast. Functional copies of trnPGGG occur in the chloroplast genomes of the green algae Chaetosphaeridium and Zygnema, otherwise this gene appears to have been lost in Chara and is absent in other examined green-algal chloroplast genomes . The GGG codon is used in coding regions of the chloroplast with roughly equal frequency in all of these lineages.
Relationships to desiccation-tolerance
Protection of chloroplast structure appears to be a major aspect of the mechanism for desiccation tolerance in bryophytes (for review see Oliver ). It is unclear at the moment as to how much the chloroplast genome, and the genes it encodes, influences the stability of the chloroplast during dehydration, or how much it contributes to the rapid resumption of photosynthetic electron transport. However, the importance of the chloroplast in the phenotype of desiccation tolerance cannot be underestimated. This is evident in the percentage of transcripts encoding chloroplast-directed proteins that are present in the transcriptome of recovering rapid-dried gametophytes of Tortula ruralis . Of the transcripts that can be annotated in the rehydration transcriptome, 12.5% are classified as chloroplastic in the Gene Ontology (GO) functional classification scheme. One of the more prominent being the transcript that encode the Early Light Inducible Protein (ELIP) which is described as being an important protein for chloroplast structural protection during desiccation and upon rehydration of T. ruralis gametophytes . Unfortunately, because of the nature of plastid transcripts (lack of a polyA tail) they are generally not present in the cDNA preparations used for transcriptome sequencing projects and thus we know little of how their synthesis responds to desiccation and rehydration or what is the role of the chloroplast transcriptome in the mechanism of desiccation tolerance we see in T. ruralis. Sequencing the chloroplast genome is an important first step in answering these questions about the interplay between the two genomes, nuclear and plastid, in this important trait. The role of the chloroplast genome and the expression of its genes remains a fertile area for study.
The Tortula ruralis genome differs from that of the only other published chloroplast genome sequence, that of Physcomitrella patens. The differences in chloroplast genome structure and gene content offer some tantalizing hypotheses in relation to one of the fundamental differences between the two mosses in their ability to tolerate the stresses associated with extreme water loss. The most intriguing observation is the loss of the petN gene, presumably to the nuclear genome, given its important role in photosynthetic electron transport and the significance of this with regards to desiccation tolerance.
The polymorphisms we have uncovered in the sequencing of the genome offer the possibility of the generating chloroplast DNA markers for future fine-level phylogenetic studies, or for future population genetic studies examining individual variation within populations.
DNA Isolation and Sequencing
Material of Tortula ruralis for DNA extraction was grown in sterile culture from spores collected from a wild population located under a pine canopy along the southern bank of the Bow River west of Calgary, Alberta, Canada, approximately 51° 06' 04° N, 114° 17' 10° W (voucher specimen deposited in the University and Jepson Herbaria, University of California, Berkeley, CA). Isolated sporophytes with intact capsules were carefully removed from the parent gametophyte, washed twice in 5% bleach solution each followed by a sterile water wash, and placed on sterile minimal media and the spores removed by breaking open the capsules. The moss cultures used in the preparation of chloroplasts for genome sequencing were derived from five different capsules. Isolation of intact chloroplasts was achieved using fluorescence-activated cell sorting (FACS) wherein fluorescently stained organelles can be visualized, separated, and collected using flow cytometry. Total chloroplast DNA was amplified using rolling circle amplification (RCA) with random hexamer primers . Amplified products were sheared in approximately equal-sized fragments of ~3000 bp by repeated passage through a narrow aperture using a Hydroshear (Genemachines, San Carlos, CA). The resulting fragments were subcloned and sequenced by the Department of Energy Joint Genome Institute (JGI). Detailed methods for the complete process can be found in Wolf et al., .
Genome Assembly, Finishing & Annotation
Because the sequences covered random sections of the chloroplast, sufficient numbers of DNA fragments were sequenced to provide an average of eight times coverage for the majority of the genome. Sequence fragments were assembled using the Phrap software package, and Consed was used to visualize assembly of contigs (see http://www.phrap.org). Remaining gaps, low-coverage regions, and areas of questionable assembly were manually resequenced from specific clones or from genomic DNA. The finishing sequences were obtained by cloning, and sequencing PCR-generated fragments obtained using flanking primers from known sequence or by use of a GenomeWalker™ kit (Clontech, Mountain View, CA) using either purified chloroplast genomic DNA or total genomic DNA extractions that contain significant levels of chloroplast DNA.
The resulting genome was annotated using the on-line Dual Organellar GenoMe Annotator (DOGMA)  and tRNAscan-SE , and gene content was compared to published annotated chloroplast genomes and in particular to the chloroplast genome of Physcomitrella patens (NC_005087). Location of introns, pseudogenes, and beginning and endpoints of many genes was aided by the libraries of plastid DNA and annotated genomes available through NCBI GenBank http://www.ncbi.nlm.nih.gov.
The inverted repeat (IR) characteristic of the chloroplast genome in embryophytes causes problems for automatic alignment programs such as Phrap, because the two IR copies cannot be differentiated by the program. However, this apparent problem can be turned into an advantage and can be used to determine the ends of the IR. Consed (an assembly viewing program) visually indicates when paired sequences from the same clone are placed too far apart in a contig, making the general area of the IR visually apparent if sufficient reads are presented. A consensus sequence for this region with an additional ~2 kb on either side was excised, and all reads were then aligned to this sequence. Using the "color means match" option in Consed, the ends of the IR were visually determined by finding the points at which there was a transition from all reads agreeing to the presence of two distinct sequence motifs.
This research was supported in part by a collaborative grant from the US National Science Foundation: ATOL: Collaborative Research: Deep Green Plant Phylogenetics: Novel Analytical Methods for Scaling from Genomics to Morphology http://ucjeps.berkeley.edu/TreeofLife/, NSF grant numbers 0228729 to BDM and JLB, 0228660 to DFM, and 0228432 to PGW. Additional funding was provided by the US National Institutes of Health Interdisciplinary Training in Genomic Sciences Grant T32-HG00035 to KGK and USDA-CSREES NRI grant 2007-02007 to MJO. The authors would like to thank Dr Aru K. Arumuganathan, Director of the Flow Cytometry Core, Benaroya Research Institute at Virginia Mason for his excellent technical assistance in the isolation of intact chloroplasts using fluorescence-activated cell sorting (FACS). We would also like to acknowledge the excellent technical assistance of Dean Kelch, Jeremy Hudgeons, and Jim Elder. Work presented here was also supported in part by USDA-ARS CRIS project 3622-21000-027-00 (MJO). This work was also partly performed under the auspices of the US Department of Energy's Office of Science, Biological and Environmental Research Program, and by the University of California, Lawrence Berkeley National Laboratory under Contract No. DE-AC02-05CH11231. Mention of a trademark or proprietary product does not constitute a guarantee or warranty of the product by the United States Department of Agriculture, and does not imply its approval to the exclusion of other products that may also be suitable.
- Mishler BD:Syntrichia. Flora of North America North of Mexico. 2007, the Flora of North America Editorial Committee: Oxford University Press, USA, 27: 618-627.Google Scholar
- Bewley JD: Physiological aspects of desiccation-tolerance. Ann Rev Plant Physiol. 1979, 30: 195-238.View ArticleGoogle Scholar
- Oliver MJ, Bewley JD: Desiccation-tolerance of plant tissues: A mechanistic overview. Hort Reviews. 1997, 18: 171-214.Google Scholar
- Winner WE, Bewley JD: Terrestrial mosses as bioindicators of SO2 pollution stress. Oecologia. 1978, 35: 221-230.View ArticleGoogle Scholar
- Naszradi T, Badacsonyi A, Németh N, Tuba Z, Bati F: Zinc, Lead, and Cadmium content in meadow plants and mosses along the M3 motorway (Hungary). J Atmos Chem. 2004, 49: 593-603.View ArticleGoogle Scholar
- Alpert P, Oliver MJ: Drying without Dying. Desiccation and survival in plants: Drying without Dying. Edited by: Black M, Pritchard HW. 2002, CABI Publishing, Wallingford, Oxon, 3-43.View ArticleGoogle Scholar
- Oliver MJ, Mishler BD, Quisenberry JE: Comparative measures of desiccation-tolerance in the Tortula ruralis complex. I. Variation in damage control and repair. Am J Bot. 1993, 80: 127-136.View ArticleGoogle Scholar
- Wood AJ: The nature and distribution of vegetative desiccation-tolerance in hornworts, liverworts and mosses. The Bryologist. 2007, 110: 163-177.View ArticleGoogle Scholar
- Proctor MCF, Oliver MJ, Wood AJ, Alpert P, Stark LR, Cleavitt N, Mishler BD: Desiccation Tolerance in Bryophytes. The Bryologist. 2007, 110: 595-621.View ArticleGoogle Scholar
- Oliver MJ, Velten J, Mishler BD: Desiccation tolerance in bryophytes: a reflection of the primitive strategy for plant survival in dehydrating habitats?. Integ Comp Biol. 2005, 45: 788-799.View ArticleGoogle Scholar
- Tuba Z, Csintalan Z, Proctor MCF: Photosynthetic responses of a moss, Tortula ruralis, ssp. ruralis, and the lichens Cladonia convoluta and C. furcata to water deficit and short periods of desiccation, and their ecophysiological significance: a baseline study at present-day CO2 concentration. New Phytol. 1996, 133: 353-361.View ArticleGoogle Scholar
- Platt KA, Oliver MJ, Thomson WW: Membranes and organelles of dehydrated Selaginella and Tortula retain their normal configuration and structural integrity: freeze fracture evidence. Protoplasma. 1994, 178: 57-65.View ArticleGoogle Scholar
- Kelch DG, Driskell A, Mishler BD: Inferring phylogeny using genomic characters: a case study using land plant plastomes. Molecular Systematics of Bryophytes [Monographs in Systematic Botany 98]. Edited by: Goffinet B, Hollowell V, Magill R. 2004, Missouri Botanical Garden Press, St. Louis, 3-12.Google Scholar
- Qiu Y-L, Li L, Wang B, Chen Z, Knoop V, Groth-Malonek M, Dombrovska O, Lee J, Kent L, Rest J, Estabrook GF, Hendry TA, Taylor DW, Testa CM, AMbros M, Crandall-Stotler B, Duff RJ, Stech M, Frey W, Quandt D, Davis CC: The deepest divergences in land plants inferred from phylogenomic evidence. PNAS. 2006, 103: 15511-1551.PubMed CentralPubMedView ArticleGoogle Scholar
- Mishler BD, Oliver MJ: Putting Physcomitrella patens on the Tree of Life: The evolution and ecology of mosses. The Moss Physcomitrella patens: Annual Plant Reviews. Edited by: Knight C, Perroud P-F, Cove D. 2009, Wiley-Blackwell, New Jersey, 36: 1-15.View ArticleGoogle Scholar
- Sugiura C, Kobayashi Y, Aoki S, Sugita C, Sugita M: Complete chloroplast DNA sequence of the moss Physcomitrella patens : evidence for the loss and relocation of rpoA from the chloroplast to the nucleus. Nuc Acids Res. 2003, 31: 5324-5331.View ArticleGoogle Scholar
- Goffinet B, Wickett N, Werner O, Ros RM, Shaw AJ, Cox CJ: Distribution and phylogenetic significance of a 71 kb inversion in the chloroplast genome of the Funariidae (Bryophyta). Annals of Botany. 2007, 99: 747-753.PubMed CentralPubMedView ArticleGoogle Scholar
- Mishler BD, Kelch DG: Phylogenomics and early land plant evolution. Bryophyte Biology. Edited by: Shaw AJ, Goffinet B. 2009, Cambridge University Press, New York, 173-197. SecondGoogle Scholar
- Allen JF: Cytochrome b6f: structure for signalling and vectorial metabolism. Trends Plant Sci. 2004, 3: 130-137.View ArticleGoogle Scholar
- Mccauley DE, Stevens JE, Peroni PA, Raveill JA: The spatial distribution of chloroplast DNA and allozyme polymorphisms within a population of Silene alba (Caryophyllaceae). Am J Bot. 1996, 83: 727-731.View ArticleGoogle Scholar
- Brown RC, Lemmon BL: Monoplastidic cell division in lower land plants. Am J Bot. 1990, 77: 559-571.View ArticleGoogle Scholar
- Renzaglia KS, Brown RC, Lemmon BE, Duckett JG, Ligrone R: The occurrence and phylogenetic significance of monoplastidic meiosis in liverworts. Can J Bot. 1994, 72: 65-7.View ArticleGoogle Scholar
- Guo XY, Ruan SL, Hu WM, Ca DG, Fan LJ: Chloroplast DNA insertions into the nuclear genome of rice: the genes, sites and ages of insertion involved. Functional & Integrative Genomics. 2008, 8: 101-108.View ArticleGoogle Scholar
- Barthet MM, Hilu KW: Expression of matK : functional and evolutionary implications. Am J Bot. 2007, 94: 1402-1412.PubMedView ArticleGoogle Scholar
- Lowe TM, Eddy SR: tRNAscan-SE: a program for improved detection of transfer RNA genes in genomic sequence. Nucleic Acids Res. 1997, 25: 955-964.PubMed CentralPubMedView ArticleGoogle Scholar
- Turmel M, Otis C, Lemieux C: The chloroplast genome sequence of Chara vulgaris sheds new light into the closest green algal relatives of land plants. Mol Biol Evol. 2006, 23: 1324-1338.PubMedView ArticleGoogle Scholar
- Oliver MJ: Lessons on dehydration tolerance from desiccation tolerant plants. Plant Desiccation Tolerance. Edited by: Jenks MA, Wood AJ. 2007, Wiley-Blackwell, New Jersey, 11-50.View ArticleGoogle Scholar
- Oliver MJ, Dowd SE, Zaragoza J, Mauget SA, Payton PR: The rehydration transcriptome of the desiccation-tolerant bryophyte Tortula ruralis : Transcript classification and analysis. BMC Genomics. 2004, 5 (89): 1-19.Google Scholar
- Zeng Q, Chen X, Wood AJ: Two early light-inducible protein (ELIP) cDNAs from the resurrection plant Tortula ruralis are differentially expressed in response to desiccation, rehydration, salinity, and high light. J Exp Bot. 2002, 53: 1197-1205.PubMedView ArticleGoogle Scholar
- Dean FB, Nelson JR, Giesler TL, Lasken RS: Rapid amplification of plasmid and phage DNA using Phi 29 DNA polymerase and multiply-primed rolling circle. Genome Res. 2001, 11: 1095-1099.PubMed CentralPubMedView ArticleGoogle Scholar
- Wolf PG, Karol KG, Mandoli DF, Kuehl J, Arumuganathan K, Ellis MW, Mishler BD, Kelch DG, Olmstead RG, Boore JL: The first complete chloroplast genome sequence of a lycophyte, Huperzia lucidula (Lycopodiaceae). Gene. 2005, 350: 117-128.PubMedView ArticleGoogle Scholar
- Wyman SK, Boore JL, Jansen RK: Automatic annotation of organellar genomes with DOGMA. Bioinformatics. 2004, 20: 3252-3255.PubMedView ArticleGoogle Scholar
- National Center for Biotechnology Information's Genbank Database. [http://www.ncbi.nlm.nih.gov]