Mitochondrial DNA of Vitis vinifera and the issue of rampant horizontal gene transfer

Substoichiometrically different mitotypes coexist in mitochondrial genomes of Brassica napus L

Cytoplasmic male sterility (CMS) has been identified in numerous plant species. Brassica napus CMS plants, such as Polima (pol), MI, and Shaan 2A, have been identified independently by different researchers with different materials in conventional breeding processes. How this kind of CMS emerges is unclear. Here, we report the mitochondrial genome sequence of the prevalent mitotype in the most widely used pol-CMS line, which has a length of 223,412 bp and encodes 34 proteins, 3 ribosomal RNAs, and 18 tRNAs, including two near identical copies of trnH. Of these 55 genes, 48 were found to be identical to their equivalents in the “nap” cytoplasm. The nap mitotype carries only one copy of trnH, and the sequences of five of the six remaining genes are highly similar to their equivalents in the pol mitotype. Forty-four open reading frames (ORFs) with unknown function were detected, including two unique to the pol mitotype (orf122 and orf132). At least five rearrangement events are required to account for the structural differences between the pol and nap sequences. The CMS-related orf224 neighboring region (∼5 kb) rearranged twice. PCR profiling based on mitotype-specific primer pairs showed that both mitotypes are present in B. napus cultivars. Quantitative PCR showed that the pol cytoplasm consists mainly of the pol mitotype, and the nap mitotype is the main genome of nap cytoplasm. Large variation in the copy number ratio of mitotypes was found, even among cultivars sharing the same cytoplasm. The coexistence of mitochondrial mitotypes and substoichiometric shifting can explain the emergence of CMS in B. napus.

The mitochondrial genome of the legume Vigna radiata and the analysis of recombination across short mitochondrial repeats

The mitochondrial genomes of seed plants are exceptionally fluid in size, structure, and sequence content, with the accumulation and activity of repetitive sequences underlying much of this variation. We report the first fully sequenced mitochondrial genome of a legume, Vigna radiata (mung bean), and show that despite its unexceptional size (401,262 nt), the genome is unusually depauperate in repetitive DNA and "promiscuous" sequences from the chloroplast and nuclear genomes. Although Vigna lacks the large, recombinationally active repeats typical of most other seed plants, a PCR survey of its modest repertoire of short (38–297 nt) repeats nevertheless revealed evidence for recombination across all of them. A set of novel control assays showed, however, that these results could instead reflect, in part or entirely, artifacts of PCR-mediated recombination. Consequently, we recommend that other methods, especially high-depth genome sequencing, be used instead of PCR to infer patterns of plant mitochondrial recombination. The average-sized but repeat- and feature-poor mitochondrial genome of Vigna makes it ever more difficult to generalize about the factors shaping the size and sequence content of plant mitochondrial genomes.

Role of intercompartmental DNA transfer in producing genetic diversity

In eukaryotic cells, genes are found in three compartments-the nucleus, mitochondria, and plastids-and extensive gene transfer has occurred between them. Most organellar genes in the nucleus migrated there long ago, but transfer is ongoing and ubiquitous. It now generates mostly noncoding nuclear DNA, can also disrupt gene functions, and reshape genes by adding novel exons. Plastid or nuclear sequences have also contributed to the formation of mitochondrial tRNA genes. It is now clear that organelle-to-nucleus DNA transfer involves the escape of DNA molecules from the organelles at times of stress or at certain developmental stages, and their subsequent incorporation at sites of double-stranded breaks in nuclear DNA by nonhomologous recombination. Intercompartmental DNA transfer thus appears to be an inescapable phenomenon that has had a broad impact on eukaryotic evolution, affecting DNA repair, gene and genome evolution, and redirecting proteins to different target compartments.

Horizontal acquisition of multiple mitochondrial genes from a parasitic plant followed by gene conversion with host mitochondrial genes

BACKGROUND

Horizontal gene transfer (HGT) is relatively common in plant mitochondrial genomes but the mechanisms, extent and consequences of transfer remain largely unknown. Previous results indicate that parasitic plants are often involved as either transfer donors or recipients, suggesting that direct contact between parasite and host facilitates genetic transfer among plants.

RESULTS

In order to uncover the mechanistic details of plant-to-plant HGT, the extent and evolutionary fate of transfer was investigated between two groups: the parasitic genus Cuscuta and a small clade of Plantago species. A broad polymerase chain reaction (PCR) survey of mitochondrial genes revealed that at least three genes (atp1, atp6 and matR) were recently transferred from Cuscuta to Plantago. Quantitative PCR assays show that these three genes have a mitochondrial location in the one species line of Plantago examined. Patterns of sequence evolution suggest that these foreign genes degraded into pseudogenes shortly after transfer and reverse transcription (RT)-PCR analyses demonstrate that none are detectably transcribed. Three cases of gene conversion were detected between native and foreign copies of the atp1 gene. The identical phylogenetic distribution of the three foreign genes within Plantago and the retention of cytidines at ancestral positions of RNA editing indicate that these genes were probably acquired via a single, DNA-mediated transfer event. However, samplings of multiple individuals from two of the three species in the recipient Plantago clade revealed complex and perplexing phylogenetic discrepancies and patterns of sequence divergence for all three of the foreign genes.

CONCLUSIONS

This study reports the best evidence to date that multiple mitochondrial genes can be transferred via a single HGT event and that transfer occurred via a strictly DNA-level intermediate. The discovery of gene conversion between co-resident foreign and native mitochondrial copies suggests that transferred genes may be evolutionarily important in generating mitochondrial genetic diversity. Finally, the complex relationships within each lineage of transferred genes imply a surprisingly complicated history of these genes in Plantago subsequent to their acquisition via HGT and this history probably involves some combination of additional transfers (including intracellular transfer), gene duplication, differential loss and mutation-rate variation. Unravelling this history will probably require sequencing multiple mitochondrial and nuclear genomes from Plantago. See Commentary: http://www.biomedcentral.com/1741-7007/8/147.

Bioinformatic analysis of ESTs collected by Sanger and pyrosequencing methods for a keystone forest tree species: oak

BACKGROUND

The Fagaceae family comprises about 1,000 woody species worldwide. About half belong to the Quercus family. These oaks are often a source of raw material for biomass wood and fiber. Pedunculate and sessile oaks, are among the most important deciduous forest tree species in Europe. Despite their ecological and economical importance, very few genomic resources have yet been generated for these species. Here, we describe the development of an EST catalogue that will support ecosystem genomics studies, where geneticists, ecophysiologists, molecular biologists and ecologists join their efforts for understanding, monitoring and predicting functional genetic diversity.

RESULTS

We generated 145,827 sequence reads from 20 cDNA libraries using the Sanger method. Unexploitable chromatograms and quality checking lead us to eliminate 19,941 sequences. Finally a total of 125,925 ESTs were retained from 111,361 cDNA clones. Pyrosequencing was also conducted for 14 libraries, generating 1,948,579 reads, from which 370,566 sequences (19.0%) were eliminated, resulting in 1,578,192 sequences. Following clustering and assembly using TGICL pipeline, 1,704,117 EST sequences collapsed into 69,154 tentative contigs and 153,517 singletons, providing 222,671 non-redundant sequences (including alternative transcripts). We also assembled the sequences using MIRA and PartiGene software and compared the three unigene sets. Gene ontology annotation was then assigned to 29,303 unigene elements. Blast search against the SWISS-PROT database revealed putative homologs for 32,810 (14.7%) unigene elements, but more extensive search with Pfam, Refseq_protein, Refseq_RNA and eight gene indices revealed homology for 67.4% of them. The EST catalogue was examined for putative homologs of candidate genes involved in bud phenology, cuticle formation, phenylpropanoids biosynthesis and cell wall formation. Our results suggest a good coverage of genes involved in these traits. Comparative orthologous sequences (COS) with other plant gene models were identified and allow to unravel the oak paleo-history. Simple sequence repeats (SSRs) and single nucleotide polymorphisms (SNPs) were searched, resulting in 52,834 SSRs and 36,411 SNPs. All of these are available through the Oak Contig Browser http://genotoul-contigbrowser.toulouse.inra.fr:9092/Quercus_robur/index.html.

CONCLUSIONS

This genomic resource provides a unique tool to discover genes of interest, study the oak transcriptome, and develop new markers to investigate functional diversity in natural populations.

Identification of new polymorphic regions and differentiation of cultivated olives (Olea europaea L.) through plastome sequence comparison

BACKGROUND

The cultivated olive (Olea europaea L.) is the most agriculturally important species of the Oleaceae family. Although many studies have been performed on plastid polymorphisms to evaluate taxonomy, phylogeny and phylogeography of Olea subspecies, only few polymorphic regions discriminating among the agronomically and economically important olive cultivars have been identified. The objective of this study was to sequence the entire plastome of olive and analyze many potential polymorphic regions to develop new inter-cultivar genetic markers.

RESULTS

The complete plastid genome of the olive cultivar Frantoio was determined by direct sequence analysis using universal and novel PCR primers designed to amplify all overlapping regions. The chloroplast genome of the olive has an organisation and gene order that is conserved among numerous Angiosperm species and do not contain any of the inversions, gene duplications, insertions, inverted repeat expansions and gene/intron losses that have been found in the chloroplast genomes of the genera Jasminum and Menodora, from the same family as Olea.The annotated sequence was used to evaluate the content of coding genes, the extent, and distribution of repeated and long dispersed sequences and the nucleotide composition pattern. These analyses provided essential information for structural, functional and comparative genomic studies in olive plastids. Furthermore, the alignment of the olive plastome sequence to those of other varieties and species identified 30 new organellar polymorphisms within the cultivated olive.

CONCLUSIONS

In addition to identifying mutations that may play a functional role in modifying the metabolism and adaptation of olive cultivars, the new chloroplast markers represent a valuable tool to assess the level of olive intercultivar plastome variation for use in population genetic analysis, phylogenesis, cultivar characterisation and DNA food tracking.

Extensive loss of translational genes in the structurally dynamic mitochondrial genome of the angiosperm Silene latifolia

Background

Mitochondrial gene loss and functional transfer to the nucleus is an ongoing process in many lineages of plants, resulting in substantial variation across species in mitochondrial gene content. The Caryophyllaceae represents one lineage that has experienced a particularly high rate of mitochondrial gene loss relative to other angiosperms.

Results

In this study, we report the first complete mitochondrial genome sequence from a member of this family, Silene latifolia. The genome can be mapped as a 253,413 bp circle, but its structure is complicated by a large repeated region that is present in 6 copies. Active recombination among these copies produces a suite of alternative genome configurations that appear to be at or near "recombinational equilibrium". The genome contains the fewest genes of any angiosperm mitochondrial genome sequenced to date, with intact copies of only 25 of the 41 protein genes inferred to be present in the common ancestor of angiosperms. As observed more broadly in angiosperms, ribosomal proteins have been especially prone to gene loss in the S. latifolia lineage. The genome has also experienced a major reduction in tRNA gene content, including loss of functional tRNAs of both native and chloroplast origin. Even assuming expanded wobble-pairing rules, the mitochondrial genome can support translation of only 17 of the 61 sense codons, which code for only 9 of the 20 amino acids. In addition, genes encoding 18S and, especially, 5S rRNA exhibit exceptional sequence divergence relative to other plants. Divergence in one region of 18S rRNA appears to be the result of a gene conversion event, in which recombination with a homologous gene of chloroplast origin led to the complete replacement of a helix in this ribosomal RNA.

Conclusions

These findings suggest a markedly expanded role for nuclear gene products in the translation of mitochondrial genes in S. latifolia and raise the possibility of altered selective constraints operating on the mitochondrial translational apparatus in this lineage.

Lineage-specific group II intron gains and losses of the mitochondrial rps3 gene in gymnosperms

According to PCR assays and sequencing, we now report the shared presence of two rps3 introns, namely the rps3i74 and the rps3i249, in the mitochondria of all the classes representing the surviving lineages of gymnosperms, and unveil several lineages experiencing intron loss. Interestingly, the rps3 intron gains and losses within the four groups of gymnosperms let us sort out the Pinaceae and the non-Pinaceae into intron (+)- and intron (-)-lineages, respectively. Worthy of mention is also the finding that only Gnetum within the Gnetales harbours both the rps3 introns. This intron distribution pattern is consistent with the hypothesis that the two rps3 introns were likely present in the common ancestor of the seed plants and, then, independently lost in the non-Pinaceae during gymnosperm evolution. The derived secondary structural model of the novel group IIA intron improves our understanding of the significance and origin of the extraordinary length polymorphisms observed among rps3i249 orthologs. Despite the remarkable structural plasticity to adopt and reject introns, the rps3 mRNAs undergo accurate processing by splicing and extensive editing in gymnosperm mitochondria. This study provides additional insights into the evolutionarily high dynamics of mitochondrial introns which may come and go in closely related plant species. The turnover of the mitochondrial rps3 group II introns seen among lineages of seed plants further suggests that these introns might be an additional signature to discriminate between particularly cryptical taxonomic groups for which there is a need of a further evaluation of their evolutionary affiliation.

Cost of Having the Largest Mitochondrial Genome: Evolutionary Mechanism of Plant Mitochondrial Genome

The angiosperm mitochondrial genome is the largest and least gene-dense among the eukaryotes, because its intergenic regions are expanded. There seems to be no functional constraint on the size of the intergenic regions; angiosperms maintain the large mitochondrial genome size by a currently unknown mechanism. After a brief description of the angiosperm mitochondrial genome, this review focuses on our current knowledge of the mechanisms that control the maintenance and alteration of the genome. In both processes, the control of homologous recombination is crucial in terms of site and frequency. The copy numbers of various types of mitochondrial DNA molecules may also be controlled, especially during transmission of the mitochondrial genome from one generation to the next. An important characteristic of angiosperm mitochondria is that they contain polypeptides that are translated from open reading frames created as byproducts of genome alteration and that are generally nonfunctional. Such polypeptides have potential to evolve into functional ones responsible for mitochondrially encoded traits such as cytoplasmic male sterility or may be remnants of the former functional polypeptides.

Large-scale detection and analysis of RNA editing in grape mtDNA by RNA deep-sequencing

RNA editing is a widespread post-transcriptional molecular phenomenon that can increase proteomic diversity, by modifying the sequence of completely or partially non-functional primary transcripts, through a variety of mechanistically and evolutionarily unrelated pathways. Editing by base substitution has been investigated in both animals and plants. However, conventional strategies based on directed Sanger sequencing are time-consuming and effectively preclude genome wide identification of RNA editing and assessment of partial and tissue-specific editing sites. In contrast, the high-throughput RNA-Seq approach allows the generation of a comprehensive landscape of RNA editing at the genome level. Short reads from Solexa/Illumina GA and ABI SOLiD platforms have been used to investigate the editing pattern in mitochondria of Vitis vinifera providing significant support for 401 C-to-U conversions in coding regions and an additional 44 modifications in non-coding RNAs. Moreover, 76% of all C-to-U conversions in coding genes represent partial RNA editing events and 28% of them were shown to be significantly tissue specific. Solexa/Illumina and SOLiD platforms showed different characteristics with respect to the specific issue of large-scale editing analysis, and the combined approach presented here reduces the false positive rate of discovery of editing events.

OrgConv: detection of gene conversion using consensus sequences and its application in plant mitochondrial and chloroplast homologs

Background

The ancestry of mitochondria and chloroplasts traces back to separate endosymbioses of once free-living bacteria. The highly reduced genomes of these two organelles therefore contain very distant homologs that only recently have been shown to recombine inside the mitochondrial genome. Detection of gene conversion between mitochondrial and chloroplast homologs was previously impossible due to the lack of suitable computer programs. Recently, I developed a novel method and have, for the first time, discovered recurrent gene conversion between chloroplast mitochondrial genes. The method will further our understanding of plant organellar genome evolution and help identify and remove gene regions with incongruent phylogenetic signals for several genes widely used in plant systematics. Here, I implement such a method that is available in a user friendly web interface.

Results

OrgConv (Organellar Conversion) is a computer package developed for detection of gene conversion between mitochondrial and chloroplast homologous genes. OrgConv is available in two forms; source code can be installed and run on a Linux platform and a web interface is available on multiple operating systems. The input files of the feature program are two multiple sequence alignments from different organellar compartments in FASTA format. The program compares every examined sequence against the consensus sequence of each sequence alignment rather than exhaustively examining every possible combination. Making use of consensus sequences significantly reduces the number of comparisons and therefore reduces overall computational time, which allows for analysis of very large datasets. Most importantly, with the significantly reduced number of comparisons, the statistical power remains high in the face of correction for multiple tests.

Conclusions

Both the source code and the web interface of OrgConv are available for free from the OrgConv website http://www.indiana.edu/~orgconv. Although OrgConv has been developed with main focus on detection of gene conversion between mitochondrial and chloroplast genes, it may also be used for detection of gene conversion between any two distinct groups of homologous sequences.

Molecular poltergeists: mitochondrial DNA copies (numts) in sequenced nuclear genomes

The natural transfer of DNA from mitochondria to the nucleus generates nuclear copies of mitochondrial DNA (numts) and is an ongoing evolutionary process, as genome sequences attest. In humans, five different numts cause genetic disease and a dozen human loci are polymorphic for the presence of numts, underscoring the rapid rate at which mitochondrial sequences reach the nucleus over evolutionary time. In the laboratory and in nature, numts enter the nuclear DNA via non-homolgous end joining (NHEJ) at double-strand breaks (DSBs). The frequency of numt insertions among 85 sequenced eukaryotic genomes reveal that numt content is strongly correlated with genome size, suggesting that the numt insertion rate might be limited by DSB frequency. Polymorphic numts in humans link maternally inherited mitochondrial genotypes to nuclear DNA haplotypes during the past, offering new opportunities to associate nuclear markers with mitochondrial markers back in time.

Insights into the evolution of mitochondrial genome size from complete sequences of Citrullus lanatus and Cucurbita pepo (Cucurbitaceae)

The mitochondrial genomes of seed plants are unusually large and vary in size by at least an order of magnitude. Much of this variation occurs within a single family, the Cucurbitaceae, whose genomes range from an estimated 390 to 2,900 kb in size. We sequenced the mitochondrial genomes of Citrullus lanatus (watermelon: 379,236 nt) and Cucurbita pepo (zucchini: 982,833 nt)--the two smallest characterized cucurbit mitochondrial genomes--and determined their RNA editing content. The relatively compact Citrullus mitochondrial genome actually contains more and longer genes and introns, longer segmental duplications, and more discernibly nuclear-derived DNA. The large size of the Cucurbita mitochondrial genome reflects the accumulation of unprecedented amounts of both chloroplast sequences (>113 kb) and short repeated sequences (>370 kb). A low mutation rate has been hypothesized to underlie increases in both genome size and RNA editing frequency in plant mitochondria. However, despite its much larger genome, Cucurbita has a significantly higher synonymous substitution rate (and presumably mutation rate) than Citrullus but comparable levels of RNA editing. The evolution of mutation rate, genome size, and RNA editing are apparently decoupled in Cucurbitaceae, reflecting either simple stochastic variation or governance by different factors.

Discovery of the rpl10 gene in diverse plant mitochondrial genomes and its probable replacement by the nuclear gene for chloroplast RPL10 in two lineages of angiosperms

Mitochondrial genomes of plants are much larger than those of mammals and often contain conserved open reading frames (ORFs) of unknown function. Here, we show that one of these conserved ORFs is actually the gene for ribosomal protein L10 (rpl10) in plant. No rpl10 gene has heretofore been reported in any mitochondrial genome other than the exceptionally gene-rich genome of the protist Reclinomonas americana. Conserved ORFs corresponding to rpl10 are present in a wide diversity of land plant and green algal mitochondrial genomes. The mitochondrial rpl10 genes are transcribed in all nine land plants examined, with five seed plant genes subject to RNA editing. In addition, mitochondrial-rpl10-like cDNAs were identified in EST libraries from numerous land plants. In three lineages of angiosperms, rpl10 is either lost from the mitochondrial genome or a pseudogene. In two of them (Brassicaceae and monocots), no nuclear copy of mitochondrial rpl10 is identifiably present, and instead a second copy of nuclear-encoded chloroplast rpl10 is present. Transient assays using green fluorescent protein indicate that this duplicate gene is dual targeted to mitochondria and chloroplasts. We infer that mitochondrial rpl10 has been functionally replaced by duplicated chloroplast counterparts in Brassicaceae and monocots.

Ribosomal protein L10 is encoded in the mitochondrial genome of many land plants and green algae

Background

The mitochondrial genomes of plants generally encode 30-40 identified protein-coding genes and a large number of lineage-specific ORFs. The lack of wide conservation for most ORFs suggests they are unlikely to be functional. However, an ORF, termed orf-bryo1, was recently found to be conserved among bryophytes suggesting that it might indeed encode a functional mitochondrial protein.

Results

From a broad survey of land plants, we have found that the orf-bryo1 gene is also conserved in the mitochondria of vascular plants and charophycean green algae. This gene is actively transcribed and RNA edited in many flowering plants. Comparative sequence analysis and distribution of editing suggests that it encodes ribosomal protein L10 of the large subunit of the ribosome. In several lineages, such as crucifers and grasses, where the rpl10 gene has been lost from the mitochondrion, we suggest that a copy of the nucleus-encoded chloroplast-derived rpl10 gene may serve as a functional replacement.

Conclusion

Despite the fact that there are now over 20 mitochondrial genome sequences for land plants and green algae, this gene has remained unidentified and largely undetected until now because of the unlikely coincidence that most of the earlier sequences were from the few lineages that lack the intact gene. These results illustrate the power of comparative sequencing to identify novel genomic features.

Fine-scale mergers of chloroplast and mitochondrial genes create functional, transcompartmentally chimeric mitochondrial genes

The mitochondrial genomes of flowering plants possess a promiscuous proclivity for taking up sequences from the chloroplast genome. All characterized chloroplast integrants exist apart from native mitochondrial genes, and only a few, involving chloroplast tRNA genes that have functionally supplanted their mitochondrial counterparts, appear to be of functional consequence. We developed a novel computational approach to search for homologous recombination (gene conversion) in a large number of sequences and applied it to 22 mitochondrial and chloroplast gene pairs, which last shared common ancestry some 2 billion years ago. We found evidence of recurrent conversion of short patches of mitochondrial genes by chloroplast homologs during angiosperm evolution, but no evidence of gene conversion in the opposite direction. All 9 putative conversion events involve the atp1/atpA gene encoding the alpha subunit of ATP synthase, which is unusually well conserved between the 2 organelles and the only shared gene that is widely sequenced across plant mitochondria. Moreover, all conversions were limited to the 2 regions of greatest nucleotide and amino acid conservation of atp1/atpA. These observations probably reflect constraints operating on both the occurrence and fixation of recombination between ancient homologs. These findings indicate that recombination between anciently related sequences is more frequent than previously appreciated and creates functional mitochondrial genes of chimeric origin. These results also have implications for the widespread use of mitochondrial atp1 in phylogeny reconstruction.

Phylogenomic analysis demonstrates a pattern of rare and ancient horizontal gene transfer between plants and fungi

Horizontal gene transfer (HGT) describes the transmission of genetic material across species boundaries and is an important evolutionary phenomenon in the ancestry of many microbes. The role of HGT in plant evolutionary history is, however, largely unexplored. Here, we compare the genomes of six plant species with those of 159 prokaryotic and eukaryotic species and identify 1689 genes that show the highest similarity to corresponding genes from fungi. We constructed a phylogeny for all 1689 genes identified and all homolog groups available from the rice (Oryza sativa) genome (3177 gene families) and used these to define 14 candidate plant-fungi HGT events. Comprehensive phylogenetic analyses of these 14 data sets, using methods that account for site rate heterogeneity, demonstrated support for nine HGT events, demonstrating an infrequent pattern of HGT between plants and fungi. Five HGTs were fungi-to-plant transfers and four were plant-to-fungi HGTs. None of the fungal-to-plant HGTs involved angiosperm recipients. These results alter the current view of organismal barriers to HGT, suggesting that phagotrophy, the consumption of a whole cell by another, is not necessarily a prerequisite for HGT between eukaryotes. Putative functional annotation of the HGT candidate genes suggests that two fungi-to-plant transfers have added phenotypes important for life in a soil environment. Our study suggests that genetic exchange between plants and fungi is exceedingly rare, particularly among the angiosperms, but has occurred during their evolutionary history and added important metabolic traits to plant lineages.

Retention of the virus-derived sequences in the nuclear genome of grapevine as a potential pathway to virus resistance

BACKGROUND

Previous studies have revealed a wide-spread occurence of the partial and complete genomes of the reverse-transcribing pararetroviruses in the nuclear genomes of herbaceous plants. Although the absence of the virus-encoded integrases attests to the random and incidental incorporation of the viral sequences, their presence could have functional implications for the virus-host interactions.

HYPOTHESIS

Analyses of two nuclear genomes of grapevine revealed multiple events of horizontal gene transfer from pararetroviruses. The approximately 200-800 bp inserts that corresponded to partial ORFs encoding reverse transcriptase apparently derived from unknown or extinct caulimoviruses and tungroviruses, were found in 11 grapevine chromosomes. In contrast to the previous reports, no reliable cases of the inserts derived from the positive-strand RNA viruses were found. Because grapevine is known to be infected by the diverse positive-strand RNA viruses, but not pararetroviruses, we hypothesize that pararetroviral inserts have conferred host resistance to these viruses. Furthermore, we propose that such resistance involves RNA interference-related mechanisms acting via small RNA-mediated methylation of pararetroviral DNAs and/or via degradation of the viral mRNAs.

CONCLUSION

The pararetroviral sequences in plant genomes may be maintained due to the benefits of virus resistance to this class of viruses conferred by their presence. Such resistance could be particularly significant for the woody plants that must withstand years- to centuries-long virus assault. Experimental research into the RNA interference pathways involving the integrated pararetroviral inserts is required to test this hypothesis.

REVIEWERS

This article was reviewed by Arcady R. Mushegian, I. King Jordan, and Eugene V. Koonin.

Evolution of reduced and compact chloroplast genomes (cpDNAs) in gnetophytes: selection toward a lower-cost strategy

The cpDNA of Welwitschia mirabilis (the only species of Welwitschiales) was recently reported to be the most reduced and compact among photosynthetic land plants. However, cpDNAs of the other two gnetophyte lineages (viz. Ephedrales and Gnetales) have not yet been studied. It remains unclear what underlining mechanisms have downsized the cpDNA. To pin down major factors for cpDNA reduction and compaction in gnetophytes, we have determined 4 complete cpDNAs, including one from each of the 3 gnetophyte orders, Ephedra equisetina, Gnetum parvifolium, and W. mirabilis, and one from the non-Pinus Pinaceae, Keteleeria davidiana. We report that the cpDNAs of E. equisetina (109,518bp) and G.parvifolium (114,914bp) are not only smaller but more compact than that of W. mirabilis (118,919bp). The gnetophyte cpDNAs have commonly lost at least 18 genes that are retained in other seed plants. Furthermore, they have significantly biased usages of AT-rich codons and shorter introns and intergenic spaces, which are largely due to more deletions at inter-operon than intra-operon spaces and removal of segment sequences rather than single-nucleotides. We show that the reduced gnetophyte cpDNAs clearly resulted from selection for economy by deletions of genes and non-coding sequences, which then led to the compactness and the accelerated substitution rates. The smallest C-values in gnetophyte nuclear DNAs and the competitive or resource-poor situations encountered by gnetophytes further suggest a critical need for an economic strategy.