MBE Advance Access originally published online on July 8, 2008
Molecular Biology and Evolution 2008 25(9):2019-2030; doi:10.1093/molbev/msn149
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
Research Articles |
The Complete Nucleotide Sequences of the 5 Genetically Distinct Plastid Genomes of Oenothera, Subsection Oenothera: II. A Microevolutionary View Using Bioinformatics and Formal Genetic Data
,1

* Lehrstuhl für Botanik, Department Biology I, Botany, Ludwig-Maximilians-University, Munich, Germany
MIPS/IBI Institute for Bioinformatics and Systems Biology, Helmholtz Center Munich; German Research Center for Environmental Health GmbH, Neuherberg, Germany
E-mail: joerg.meurer{at}lrz.uni-muenchen.de.
| Abstract |
|---|
|
|
|---|
A unique combination of genetic features and a rich stock of information make the flowering plant genus Oenothera an appealing model to explore the molecular basis of speciation processes including nucleus–organelle coevolution. From representative species, we have recently reported complete nucleotide sequences of the 5 basic and genetically distinguishable plastid chromosomes of subsection Oenothera (I–V). In nature, Oenothera plastid genomes are associated with 6 distinct, either homozygous or heterozygous, diploid nuclear genotypes of the 3 basic genomes A, B, or C. Artificially produced plastome–genome combinations that do not occur naturally often display interspecific plastome–genome incompatibility (PGI). In this study, we compare formal genetic data available from all 30 plastome–genome combinations with sequence differences between the plastomes to uncover potential determinants for interspecific PGI. Consistent with an active role in speciation, a remarkable number of genes have high Ka/Ks ratios. Different from the Solanacean cybrid model Atropa/tobacco, RNA editing seems not to be relevant for PGIs in Oenothera. However, predominantly sequence polymorphisms in intergenic segments are proposed as possible sources for PGI. A single locus, the bidirectional promoter region between psbB and clpP, is suggested to contribute to compartmental PGI in the interspecific AB hybrid containing plastome I (AB-I), consistent with its perturbed photosystem II activity.
Key Words: Oenothera chloroplast genome cytonuclear coevolution interspecific plastome–genome incompatibility speciation Dobzhansky–Muller incompatibility microevolution
| Introduction |
|---|
|
|
|---|
The evolution of eukaryotic genomes originated in endosymbiotic cell conglomerates and was based on a conversion and an extension of the genetic potentials of initially free-living partner cells that coevolved into a single, integrated compartmentalized genetic system. To date, this machinery, with nucleus/cytosol and mitochondria in animals and fungi, and in addition with plastids in plants, is regulated spatiotemporally and quantitatively in its entirety. Mitochondria and chloroplasts possess only rudimentary genomes because they have lost a large fraction of their ancestral genes many of which by transfer to the nucleus (summarized in Herrmann [1997
Over the past century, the Onagracean genus Oenothera has developed as a model for studying plant evolution and represents an ideal taxon for application of modern molecular genetic strategies to elucidate mechanisms of speciation including aspects of the compartmentalized genetic machinery, notably of plastids. The genus represents a well-defined group of flowering plants, presumably the best-known plant taxon of its size (Raven 1988
; Levin et al. 2003
, 2004
), for which a rich stock of biosystematic, taxonomic, and genetic information is available. Oenothera genetics includes a unique combination of features, such as biparental transmission of plastids, a general interfertility of species, viable and fertile hybrid offspring, as well as permanent translocation heterozygotic genomes, generally operating in combination with a system of gametophytic or sporophytic lethal factors. Together they allow the exchange of plastids and nuclei as well as the substitution of entire haploid chromosome sets or of individual (or more) chromosome pairs between species. Such exchanges frequently result in developmentally impaired, though fertile, plastome–genome incompatible hybrids (Cleland 1972
; Stubbe and Raven 1979
; Stubbe 1989
; Harte 1994
; Dietrich et al. 1997
). Initially, Oenothera research was focused predominantly on an understanding of the nuclear genome, particularly of the phenomenon of translocation heterozygosity, its nature, origin, patterns, and impact for speciation (summarized in Cleland [1972
]; Harte [1994
]). The crucial role of the plastome in the evolution of the genus has become obvious since Renner's (1934)
fundamental genetic work. Together, these studies resulted in a greater insight into both the evolutionary structure and dynamics of the clade (Stubbe 1964
; Cleland 1972
; Stubbe and Raven 1979
; Harte 1994
; Dietrich et al. 1997
).
Subsection Oenothera (=Euoenothera) is 1 of the 5 subsections of the section Oenothera on which most of the experimental work on evolution in this taxon has been performed. It is comprised of somewhat more than a dozen species with numerous variants (Dietrich et al. 1997
). Their genotypes, 3 basic genomes (A, B, and C) in homozygous (AA, BB, CC) or heterozygous (AB, AC, BC) constitution associated with 5 genetically distinguishable plastomes (I–V) in distinct combinations, and habitats are well definable. Although all these species overlap geographically (fig. 1), are interbreeding, and their hybrids are fully fertile, the subsection has been convincingly demonstrated to represent a natural group of closely related species that can unambiguously be discerned morphologically, cytogenetically, and in particular on the basis of their relationships between their genomes and plastomes (Stubbe 1964
; Dietrich et al. 1997
; fig. 2).
|
|
To obtain a more precise understanding of speciation processes including the entire compartmentalized plant genome and of the taxonomic and evolutionary relationships of the species involved, we generated and applied an EST library, developed molecular markers of the nuclear and plastid genomes of subspecies for genotyping approaches (Mrá
ek et al. 2006| Materials and Methods |
|---|
|
|
|---|
Plant Material
The naturally occurring species Oenothera elata subsp. hookeri strain johansen (hjohansen·hjohansen Ijoh; AA-I) (Cleland 1935
Various Oenothera species and strains listed in Supplementary table 1 (Supplementary Material online) were investigated to monitor the specificity of the deletion in the clpP–psbB spacer region for plastome type I variants, that is, from 9 subplastomes I and 2 subplastomes each of plastomes II, III, and IV. The relevant sequence interval was PCR amplified with the primer pair VP9for 5'-catcttcttcgtcttcgtctcc-3' and VP10rev 5'-aatacacccaatgccagatagc-3' and sequenced.
Plastome Sequences
Sequences of the 5 basic plastomes of subsection Oenothera are deposited in GenBank (accession numbers: AJ271079.3, EU262887
[GenBank]
, EU262889
[GenBank]
, EU26890, and EU262891
[GenBank]
, respectively). Their design and differences in gene-coding and intergenic regions have been discussed previously (Greiner et al. 2008
).
Identification of Potential Editing Sites
Differences between deduced protein sequences of Oenothera and the corresponding components of Marchantia polymorpha (X04465
[GenBank]
) that lacks RNA editing (Freyer et al. 1997
) have been computationally analyzed whether a C-to-U substitution within the respective codon may restore the liverwort sequence. Such positions are designated potential editing sites in Oenothera genes.
Analysis of Variable Amino Acid Sites
To estimate the impact of a single amino acid exchange to PGI within Oenothera, 2 characteristics were considered. First, mutations exchanging amino acids with highly different biochemical properties have an increased likelihood to alter or destroy protein function. The skewness of biochemical properties at a variable site was therefore estimated within the 5 Oenothera plastomes. Second, highly conserved sites are less likely to undergo a drastic change. Hence, a background distribution of the biochemical properties of sites was derived from multiple alignments of sequences from 30 reference species covering dicots as well as monocots.
Biochemical property distributions were generated as follows: Let c be a specific column in the alignment of N species and C = {xi}, i = 1, 2, ... N, the set of letters x of the ith species at c. Given the Grantham (1974)
Matrix as a distance function d,
defines a distribution of pairwise differences of biochemical properties between all species at position c. To enrich for PGI candidate sites, we tested for significant differences (P
0.05) by a nonparametric Wilcoxon rank sum test (http://www.r-project.org/index.html) between the distributions derived from sequences of the 5 Oenothera and reference species, respectively. Note that the test excludes also sites that represent an Oenothera-specific adaptation, that is, sites that are similar within each data set but dissimilar between both sets. The Oenothera mutations were also checked whether they are located within known functional regions. Protein domains (e values
1 x10–10) were detected using the PFAM database (Bateman et al. 2004). Transmembrane domains were located by the InterPRO database (http://www.ebi.ac.uk/interpro) and analyzed with the online DAS server (Cserzo et al. 1997).
Estimation of Ka/Ks Substitution Rates
Nonsynonymous and synonymous substitution rates were determined using the yn00 program of the PAML package (Yang 1997
). F3x4 were selected as substitution matrix, and Ka and Ks were determined by the Nei–Gojobori method as implemented in yn00. Rates for protein-coding genes variable among at least 2 of the 5 Oenothera plastomes were estimated from pairwise codon-based alignments. There are 10 pairwise combinations for each gene, resulting in a total of 780 rates for all and 470 rates for variable genes. Note that the computation of
= Ka/Ks is not always applicable (e.g., for Ks = 0). Therefore,
could be determined for only 215 pairwise combinations. Mean
, Ka, and Ks values for a species were obtained from concatenated alignments of single gene alignments. Only orthologous genes present in Oenothera and all reference plastomes have been used for this approach.
Photosystem II Quantum Yield and Photosystem I Redox State
Efficiency and functional state of photosystem (PS) II reflected by chlorophyll a fluorescence parameters at room temperature (Schreiber et al. 1998
) were calculated from wild-type (AA-I and AB-III) and incompatible (AB-I) rosette leaves using a pulse-modulated fluorimeter (PAM 101, Walz, Effeltrich, Germany). The light intensity of the modulated measuring beam (1.6 kHz) was 0.5 µmol m–2 s–1. Leaves, dark adapted for 10 min, were used to detect the intrinsic (Fo) and maximal (Fm) fluorescence yields, the latter being determined by application of a saturating light pulse (0.8 s, 7,000 µmol m–2 s–1). The potential maximum quantum yield was calculated as (Fm – Fo)/Fm = Fv/Fm. Red actinic light (650 nm, 50 µmol m–2 s–1) was used for measurements of fluorescence quenching. Nonphotochemical quenching (NPQ) was determined by applying repetitive saturation pulses with 20-s intervals and calculated as (Fm' – F)/(Fm' – Fo) (Kooten and Snel 1990
).
The light-dependent redox state of PSI was measured on leaves as absorption changes at 830 nm in the absence or presence of actinic (650 nm, 50 µE m–2 s–1) and far-red light (12 Wm–2) using the PSI attachment of PAM101 (Klughammer and Schreiber 1998
). Saturating light pulses (0.8 s, 7,000 µmol m–2 s–1) were applied to follow PSII-dependent reduction of PSI in far-red background light.
| Results |
|---|
|
|
|---|
General Features
The geographical distribution and genomic constitution of the 11 genuine species of the North American subsection Oenothera, with their basic genomes A, B, and C associated with their respective plastomes I–V are illustrated in figure 1. In a preceding publication, we presented the complete nucleotide sequences of representatives of the 5 basic Oenothera plastomes, their comparison, evolutionary relationships, and temporal relation (Greiner et al. 2008
Estimation of Selection Pressure on Oenothera Plastomes
Genes causative for speciation have been suggested to be under positive selection for a limited period. Positive selection on plastid genes has not been reported, except recently for clpP and rbcL in Oenothera and other taxa (Kapralov and Filatov 2007
; Erixon and Oxelman 2008
). To estimate selection pressure on the Oenothera plastomes and to derive candidate genes for incompatibility factors, ratios of nonsynonymous (Ka) versus synonymous (Ks) substitutions were determined for genes, varying within the 5 plastomes, using alignments of their entire coding sequences. Out of 233 pairwise comparisons for which the method was applicable, 33 (14.1%) exhibited elevated Ka/Ks rates above 1.0. Notably, an excess of nonsynonymous substitutions was not equally distributed between the pairs under study but clustered predominantly to 5 genes, ycf1, ycf2, accD, clpP, and ndhA (Supplementary table 2, Supplementary Material online). For almost all plastome pairs, these genes displayed
values greater than 1, which indicates positive selection. The first 3 genes contain extended repetitive regions that are only weakly conserved in other species. It is therefore unclear whether divergence of these regions is functionally relevant or whether the increase of nonsynonymous substitution rates is simply the result of observed high variability. The highest rates (
= 4.1) calculated for clpP seem to differentiate plastomes I and II versus III, whereas maximal rates for accD (2.2 <
< 4) were observed between clade I/II and IV and V, respectively. A similar, though less pronounced, observation has been made for ndhA. Together, elevated rates in the 5 genes comprise 26 of 33 pairs with positive selection. For ccsA, petD, and matK each, only one pairwise
value exceeded the chosen criterion for positive selection. However, we note that several pairs of matK and ccsA indicate elevated evolutionary rates well above the median
= 0.405. Clustering of elevated rates has also been observed for ndhD, rps18, and rps3 (Supplementary table 2, Supplementary Material online). The mean Ka/Ks ratio
= 0.166 for pairwise comparisons of 30 angiosperm chloroplast genomes showed only a small variance (standard deviation: 0.03) in most branches of the mono- and eudicotyledonous plants investigated. Most ratios between Oenothera plastomes, however, were considerably higher (mean Ka/Ks ration
= 0,47 ± 0.11, Student's t-test: P < 10–6), with a minimum of 0.24 between plastomes IV and V and a maximum of 0.63 for plastomes I and II (table 1). The unusually high mean Ka/Ks ratios between most of the 5 plastomes and the relatively large number of genes indicating positive selection or fast evolutionary rates are consistent with a significant contribution of plastomes to speciation.
|
Search for Candidate Protein-Coding Loci Involved in PGI
Two attempts were made to estimate the impact of plastome-specific differences in genes. First, a putative functional impact of nonsynonymous sites in polypeptides was estimated from the degree of conservation and differences of biochemical properties between Oenothera and 30 reference species as described in Materials and Methods. Starting from 388 nonsynonymous replacements, excluding ycf1, ycf2, and the highly variable accD N-terminus (Greiner et al. 2008
Second, protein-coding genes with length polymorphisms are of intrinsic interest for PGI. Relevant are 2 genes, atpA and psbB, from which the electrophoretic mobility is known to differ in plastomes III, IV, and V, respectively (Herrmann and Possingham 1980
), and 8 loci that predictably should generate variant polypeptides between the evening primrose plastomes (Greiner et al. 2008
). The latter was arbitrarily divided into 2 groups, those with (ndhD, rpl22, and rps18) and those without reading frameshifts (ycf1, ycf2, accD, clpP, and ndhF). Pairwise as well as multiple protein alignments between these primrose size variants and their orthologs in 4–8 reference plastomes were manually edited and compared with the plastome–genome compatibility scheme of subsection Oenothera (fig. 2). Briefly, alternative stop codons caused by single-base pair indels were found for ndhD (plastome I) and rpl22, multiple base pair indels leading to different stop codons for accD, clpP, ndhD (plastome V), ndhF, and rps18; a point mutation causes a change of the stop codon in atpA of plastome III, and a single G (IV, V) to A (I, II, III) conversion at position 1195 in psbB leads to a conservative change from valine to isoleucine (Supplementary table 5, Supplementary Material online). The ycf1 and ycf2 genes are only moderately conserved in plastid genomes in general, and AccD is highly polymorphic in its N-terminal region as in the reference plastomes. Length polymorphisms in accD, atpA, clpP, ndhD, ndhF, rpl22, and rps18 could not be compared with distinct amino acid residues. However, none of these regions, except the altered ndhD 3' terminus, is located within a functional domain (Supplementary table 5, Supplementary Material online). Computational analysis of transmembrane domains was used also to check each candidate for aberrations in its transmembrane domain architecture. Marginal effects appeared only for components of the NAD(P)H dehydrogenase complex, for ndhA, ndhD, and ndhE (data not shown). Apart from the fact that it is unclear whether the predicted relatively small perturbations alter notably the stability of the transmembrane regions, disruption of ndh genes in tobacco did not result in a pronounced phenotype (Burrows et al. 1998
; Kofer et al. 1998
). Details of the changes in these loci including verification by PCR analysis as well as comparison with loci of reference species and their interpretation are therefore presented in Greiner et al. (2008
; supplementary figs. S1–S5, Supplementary Material online).
Supplementary tables 3 and 5 (Supplementary Material online) list the possible role of each of the significant mutations in coding regions for both single amino acid exchanges and variant polypeptides caused by indels for the distinct plastome–genome combinations in the compatibility scheme (fig. 2).
Role of RNA Editing in PGI
The findings that editotypes in plastids may differ between related species and even between ecotypes (Tsudzuki et al. 2001
; Tillich et al. 2005
) and that plant-specific editing sites often cannot be modified heterologously (Schmitz-Linneweber et al. 2005
; Shikanai 2006
) have been suggested that this kind of RNA maturation plays a crucial role not only in translation but also in speciation processes (Schmitz-Linneweber et al. 2001
, 2005
). Editing of the ndhD ACG start codon in plastomes I and IV as in spinach, tobacco, and Antirrhinum (Neckermann et al. 1994
) established the presence of an editing system in evening primroses. Comparison of the protein-coding sequences with those of the liverwort M. polymorpha, which lacks RNA editing uncovered 320 potential sites in coding regions in Oenothera (data not shown). These sites represent a comprehensive capture of all possible editing sites including true sites (in higher plants usually in the order of 30–35) and a large fraction of false positives. Despite this high number, only a single-nucleotide substitution differs among the potentially edited sites in the Onagracean plastomes. NdhA of plastomes I, II, and III contains a C-to-T conversion at amino acid position 309 compared with plastomes IV and V which would result in a T-to-I amino acid change. Because knockout lines of genes for nicotinamide adenine dinucleotide phosphate dehydrogenase (NDH) subunits in tobacco display no or only subtle phenotypes (Burrows et al. 1998
; Kofer et al. 1998
), editing—in contrast to the Solanacean model (Schmitz-Linneweber et al. 2001
, 2005
)—does not play a crucial role in speciation of the Oenothera clade nor in the generation of interspecific plastid–nuclear incompatibility.
Search for Candidate Loci for PGI in Intergenic Regions
The differences located in promoter sequences of the Oenothera plastomes were compared with the genetically determined compatibility relationships in interspecific hybrids (fig. 2) in order to classify them and to pinpoint potential determinants contributing to PGI.
Two types of RNA polymerases operate in plastids of higher plants, the ancient eubacterial type polymerase (PEP) and an acquired phage-type polymerase (NEP) (Shiina et al. 2005
; Liere and Börner 2006
). The polymerases recognize distinct sites in plastid promoters. In all, 38 putative PEP- and 25 NEP-binding sites that were computationally predicted were altered in at least one of the plastomes (Greiner et al. 2008
). Mutations were classified according to 3 criteria: 1) their similarity to an ideal consensus (Silhavy and Maliga 1998
; Kapoor and Sugiura 1999
; Homann and Link 2003
; Kanamaru and Tanaka 2004
; Shiina et al. 2005
; Liere and Börner 2006
), 2) number, and 3) position of predicted polymerase binding sites relative to a translational start site (supplementary table 4A and B, Supplementary Material online). With this selection scheme, 9 putative PEP promoters, notably of clpP, psbB, rpl16, rpl33, rps12, rps15, trnGGCC, trnLCAA, trnSUGA, and 7 predicted NEP promoters, namely of atpH, clpP, ndhG, psbB, psbK, rps4, and trnGGCC, were deduced as candidates causing PGI. Three promoter sequences indicated drastic changes for both polymerases: the promoter of trnGGCC harboring mutations specific for plastome V, and the bidirectional promoters for clpP (encoding a catalytic subunit of the protease Clp) and psbB (encoding the chlorophyll a–binding protein CP47 of the PSII core complex) that contain a large deletion specific for plastome I.
Delineation of an AB-I Incompatibility Determinant: The clpP–psbB Intergenic Region
Variant sequence regions in plastomes I–V were systematically compared with compatibility/incompatibility patterns (fig. 2) to filter loci relevant for PGI. This approach finally correlated the AB-I phenotype with a single major locus, the intergenic region psbB–clpP. AB-I displays a yellow–green (lutescent) phenotype (fig. 3A). The plastome–genome combinations AB-II and AB-III, which both occur in nature, and the hybrid containing plastome IV are green (fig. 2; Stubbe 1989
). In this distinct case, changes in plastome I shared with or similar to at least one in plastomes II–IV can be disregarded. Plastome V was excluded from our analysis because the combinations AB-V, AA-V, and BB-V are extremely disharmonic and differ substantially from the AB-I phenotype. They are fully bleached, largely pollen sterile, and display severe inhibition of cell division (Stubbe 1963
). The principal genetic determinants responsible for their phenotypes are presumably complex and different from those causing bleaching of AB-I individuals (see Discussion).
|
The outlined filtering strategy uncovered 22 regions, including minor changes such as single nucleic acid exchanges in noncoding regions. Six of them are polymorphic within all plastomes. A contribution of these regions to PGI is unlikely because plastomes II, III, and IV remain fully compatible in an AB background (fig. 2). This excludes the genes accD, ycf1, and ycf2, the intergenic regions rps16/rbcL and trnQUUG/accD, as well as the SSC/IRA junction. Sixteen regions specifically altered in plastome I remain. All regions involving the NADPH complex can be disregarded because knockouts of individual NDH subunits in tobacco lack a conspicuous phenotype (Burrows et al. 1998
|
First, sequence analysis of intraplastome variance (Herrmann et al. 1980
Second, a phylogenetic footprinting analysis of the intergenic region between both genes delimited distinct boxes in the promoter region described for clpP (Sriraman et al. 1998)
and psbB (Westhoff 1985
; Westhoff and Herrmann 1988
) and uncovered substantial changes to Oenothera and among its plastomes. Notably, a large deletion in plastome I directly upstream of the highly conserved psbB PEP promoter and the clpP 5' region eliminated a putative and a confirmed clpP promoter and 2 putative psbB promoters (data not shown) (fig. 4). Conservation as well as similarity to known sites in other species suggested that an element missing in plastome I caused the specific adaptation of plastome I to the AA genome and that the change of this region leads to an altered expression of psbB and/or clpP in the AB-I combination.
Third, if altered expression of psbB and/or clpP contributes to the bleaching of the AB-I phenotype, decreased PSII activity/levels and eventually pleiotropic effects due to clpP overexpression would be expected. Thus, the activity of PSII relative to that of PSI was monitored by spectroscopic analyses (fig. 5). Chlorophyll fluorescence induction was measured on dark-adapted leaves to adjust an oxidized plastoquinone pool of incompatible AB-I and compatible AB-III and AA-I plants. Maximum PSII quantum efficiency (Fv/Fm) was in fact reduced to 0.52 ± 0.04 in the incompatible AB-I leaf compared with 0.79 ± 0.03 in AB-III or AA-I consistent with a deficiency in PSII activity. The low Fv/Fm ratio observed was caused by an elevated Fo and decreased Fv fluorescence level in AB-I indicating inefficient exciton transfer to the reaction center and/or malfunction of PSII (fig. 5). The fluorescence level elicited by moderate actinic light (50 µE m–2 s–1) dropped far below the initial Fo level depending on the light intensity indicative defects in photosynthetic light utilization. Actinic light-induced quenching of chlorophyll fluorescence was faster in AB-I than in AB-III or AA-I, and a longer period was required to reach the steady state. NPQ increased dramatically from 0.39 ± 0.14 in AA-I and AB-III to 2.47 ± 0.42 in AB-I indicating an imbalance of photosynthetic electron transport. After switching off actinic light, the fluorescence decayed to the former Fo level in AA-I and AB-III but increased in AB-I within several minutes to finally reach the increased Fo already observed during the dark adaptation prior to the measurement, again indicating malfunction of PSII. Application of far-red light, which preferentially excites PSI had no notable effect on the half-life of the fluorescence rise upon light/dark switches nor on the elevated Fo levels suggesting that the AB-I incompatibility reflects a direct effect on PSII-driven electron transport, such as stable accumulation of QA species in the dark which induces an increased Fo level.
|
The extent of PSI oxidation in terms of balanced electron flow from PSII to PSI was monitored using absorption changes at 830 nm on actinic background light in order to further substantiate the primary lesion in the hybrid (fig. 4). Although PSI signal intensity in AB-I was generally slightly reduced, the results clearly showed a much higher actinic light-induced oxidation state of PSI in AB-I compared with AB-III. At 50 µE m–2 s–1, about 15% of P700 was oxidized in AB-I, whereas the reaction center of PSI remained almost completely reduced in AB-III, a response likely due to limited electron flow toward PSI of a deficient PSII rather than of downstream effects in the AB-I hybrid. Furthermore, short light pulses at far-red background light were sufficient to completely reduce PSI in AB-III and AA-I but not in AB-I, indicating a slower PSII-mediated electron transport. Taken together, chlorophyll fluorescence analysis is fully consistent with a decreased PSII activity relative to PSI and suggests that the lesions in the incompatible AB-I hybrid primarily affect PSII function. A detailed molecular analysis of the incompatible AB-1 phenotype will be presented in a forthcoming study.
All data obtained not shown in the manuscript are available upon request.
| Discussion |
|---|
|
|
|---|
Compartmental coevolution is accompanied by distinct changes in the 5 available Onagracean organelle chromosomes and the respective nuclear genomes. In subsection Oenothera, relationships between plastome and genome are crucial for the vitality of interspecific hybrids (Stubbe 1964
Disturbances caused by interspecific organelle exchanges normally affect a multitude of ontogenetic processes. Most conspicuously are lesions of the photosynthetic machinery (hybrid bleaching, hybrid variegation) and of the generative phase that can be impaired at various gametophytic and/or sporophytic stages, such as pollen and/or ovule development (summarized in Harte [1994
]). Both, photosynthetic performance and hybrid fertility are complex characters with a large number of potential states. Thus, PGI phenotypes are probably caused by processes that operate at many scales. This intricacy is not unexpected because about 30% of the nuclear coding capacity is required for the management of the organelles (Herrmann 1997
). Complexity also becomes apparent in various extremely disharmonic plastome–genome incompatible hybrids such as AA-V, AB-V, and BB-V, in which fertility and morphogenesis are severely affected (Stubbe 1963
; Stubbe et al. 1978
). In these cases, incompatibility may be caused by multiple Dobzhansky–Muller plastid–nuclear gene pairs. The profound knowledge of the photosynthetic process and of the biogenesis, maintenance, and adaptation of its underlying structures (Herrmann and Westhoff 2001
) can advantageously be used for evaluating phylogenetic questions and provide general and causal access to evolutionary relationships and speciation processes.
The relatively large number of plastome-encoded genes deduced to be under positive selection or displaying fast evolutionary rates, manifesting in high Ka/Ks values, clearly reflect an active contribution of plastomes to speciation. However, our knowledge to computationally predictable functional elements from primary sequences is limited, and Oenothera genome sequences—the complementary part on which prespeciation and coevolution processes are acting—are missing. Nevertheless, the alignment of sequence differences between plastomes of closely related, interbreeding species to predict genetic elements combined with filtering by their evolutionary and functional relevance as well as combinatorial logics is a promising strategy to pinpoint potential plastid-localized determinants involved in compartmental coevolution last not least because predictions can be tested. The case study of the hybrid AB-I has proven that the strategy of systematic filtering on genetically well-defined material is useful and attests to the power of the method. It correlated the bleached AB-I phenotype with a distinct major locus, a plastome I–specific deletion in the clpP–psbB intergenic region with reduced PSII activity in AB-I. Bioinformatic and biophysical data are consistent with a primary lesion in PSII and reminiscent to PSII downregulation in a bleached Arabidopsis mutant with severe deficiencies of psbB transcripts and CP47 proteins (Meurer et al. 1996
, 2002
). This mutant also displays comparable fluorescence characteristics. Loci such as the clpP–psbB intergenic region deduced by this approach are therefore potential candidates that deserve further study of underlying molecular mechanisms of PGI.
At a molecular scale, coevolution of polypeptides with their interaction partners, polypeptides with polypeptides or polypeptides with nucleic acid molecules, is a well-known phenomenon (Goh et al. 2000
). Basically, it could reflect a regulatory and/or a structural phenomenon. Diverging loci (supplementary table 5, Supplementary Material online) are therefore of intrinsic interest but the vast majority of such loci found in Oenothera either does not seem to be involved in interspecific compartmental incompatibility or does not contribute in a simple way. All changes occurred in parts of polypeptide chains, which are generally highly variable and do not affect conserved domains (supplementary figs. S1–S5 [Supplementary Material online] in Greiner et al. [2008
]). In terms of single amino acid exchanges of all 388 detected nonsynonymous substitutions, only 35 sites from which 25 are located in functional domains show biochemically significant different properties with respect to reference plastomes (supplementary table 3, Supplementary Material online). These sites are encoded by altogether 19 genes, namely accD, atpA, atpB, atpF, ccsA, clpP, matK, ndhA, ndhB, ndhC, ndhD, ndhE, ndhH, petB, rpoB, rpoC2, rps3, rps8, and rps15. Current knowledge decreases this small number of candidates even further. For instance, variance of RpoB between tobacco and Atropa, initially a candidate for PGI (Herrmann et al. 2003
), did not turn out to have a notable influence (Schmitz-Linneweber et al. 2005
), or gene deletions in tobacco suggest that disruption of genes for subunits of the NDH complex (Burrows et al. 1998
; Kofer et al. 1998
) are not of relevance either. Furthermore, a biochemically possibly significant amino acid change in petD fits with the compatibility chart in principle (fig. 2) but appears to be neutral because BB-II does not show a cytochrome phenotype (data not shown). Collectively, the available data strongly suggest that plastome-encoded genes for structural components of the photosynthetic machinery are usually highly conserved. Therefore, we suggest that compartmental coevolution in Oenothera influences predominantly gene expression and/or the transcript metabolism, as found also in the Solanacean model (Schmitz-Linneweber et al. 2001
, 2005
).
Species-specific editotype differences in plastid DNA and coevolving nuclear encoded trans factors have been shown to be important in compartmental coevolution between Atropa and tobacco. They play a crucial role in harmonious nucleo–plastid interaction between both species and also explain the pronounced phenotypic difference of their reciprocal cybrids (Schmitz-Linneweber et al. 2001
, 2002
, 2005
). RNA editing obviously does not influence compartmental coevolution in evening primroses. Thus, other aspects must cause or be involved in compartmental divergence, such as transcription and/or transcript stability of photosynthetic genes via interaction with trans factors of nuclear origin and corresponding cis elements in the psbB promoter or of stabilizing elements in 5' untranslated regions of its mRNA species in the AB-I hybrid. Collectively, these findings suggest that despite of photosynthetic defects in both, Solanaceae and Oenothera, the determinants and mechanisms causing plastid–nuclear incompatability are quite different. Obviously, the ways in which individual species or genera have evolved, their histories, and adaptation to their present-day habitats are diverse and include changes at quite different molecular levels.
Plastid chromosomes have well-defined structures, a limited coding potential and a conserved design. Their sequences have become a powerful research tool not only to probe into chloroplast biogenesis, function, and engineering but also into eukaryotic genome evolution. Classical data sets such as available for subsection Oenothera are of obvious value to study microevolutionary processes. The other subsections of the genus have not yet received a comparably intense genetic study, but quite substantial information is available regarding interfertility and plastid–nuclear compatibility relationships (Stubbe and Raven 1979
). As in subsection Oenothera, offspring of interspecific crosses within subsections Raimannia, Munzia, Nutantigemma, and Emersonia is generally fertile, and even the generation of intersubsectional hybrids is possible, but these are often extremely disharmonic and their fertility is usually reduced (Stubbe and Raven 1979
). Evening primroses represent an instructive, unique model because they offer ready, general, and causal access to study processes of microevolution not only for plastome/genome relationships but also of chromosomal evolution, pollination biology, and biogeographical radiation of this subsection in a continental dimension. With the substantial base of knowledge including taxonomic, biogeographical, genetic, and phylogenetic framework, nuclear and organelle genetics, and the current availability of molecular and cell biological strategies of analysis, Oenothera provides unique opportunities to overcome limitations of traditional phylogenetic analyses based on mere sequence data and to develop a "functional molecular phylogeny."
| Supplementary Material |
|---|
|
|
|---|
Supplementary tables 1–5 and figures S1–S5 are available at Molecular Biology and Evolution online (http://www.mbe.oxfordjournals.org/).
| Acknowledgements |
|---|
|
|
|---|
We thank Martina Silber for critical comments, Dario Leister for kindly reading the manuscript, and the Springer-Verlag for permission to adapt figure 2. The project was supported by the Deutsche Forschungsgemeinschaft (SFB-TR1) to R.G.H. and J.M. and the Hanns-Seidel-Stiftung of grants of the Bundesministerium für Bildung und Forschung to S.G. We dedicate this manuscript to the memory of Prof. Dr Rainer M. Maier.
| Footnotes |
|---|
1 These authors contributed equally to this work.
| References |
|---|
|
|
|---|
Bateman A, Coin L, Durbin R, et al, (13 co-authors). The Pfam protein families database. Nucleic Acids Res (2004) 32:D138–D141.
Burrows PA, Sazanov LA, Svab Z, Maliga P, Nixon PJ. Identification of a functional respiratory complex in chloroplasts through analysis of tobacco mutants containing disrupted plastid ndh genes. EMBO J (1998) 17:868–876.[CrossRef][Web of Science][Medline]
Chiu W-L, Stubbe W, Sears BB. Plastid inheritance in Oenothera: organelle genome modifies the extent of biparental plastid transmission. Curr Genet (1988) 13:181–189.[CrossRef][Web of Science]
Cleland RE. Cyto-taxonomic studies on certain Oenotheras from California. Proc Am Philos Soc (1935) 75:339–429.
Cleland RE. Oenothera cytogenetics and evolution. In: Experimental botany—Sutcliffe JF, Mahlberg P, eds. (1972) London: Academic Press. 43–299.
Cserzo M, Wallin E, Simon I, von Heijne G, Elofsson A. Prediction of transmembrane alpha-helices in prokaryotic membrane proteins: the dense alignment surface method. Protein Eng (1997) 10:673–676.
Dietrich W, Wagner WL, Raven PH. Systematics of Oenothera section Oenothera subsection Oenothera (Onagraceae). In: Systematic botany monographs—Anderson C, ed. (1997) Laramie (WY): The American Society of Plant Taxonomists. 1–123.
Dobzhansky T. Genetics and the origin of species (1937) New York: Columbia University Press.
Erixon P, Oxelman B. Whole-gene positive selection, elevated synonymous substitution rates, duplication, and indel evolution of the chloroplast clpP1 gene. PLoS ONE (2008) 3:e1386.[CrossRef]
Freyer R, Kiefer-Meyer MC, Kossel H. Occurrence of plastid RNA editing in all major lineages of land plants. Proc Natl Acad Sci USA (1997) 94:6285–6290.
Goh CS, Bogan AA, Joachimiak M, Walther D, Cohen FE. Co-evolution of proteins with their interaction partners. J Mol Biol (2000) 299:283–293.[CrossRef][Web of Science][Medline]
Grantham R. Amino acid difference formula to help explain protein evolution. Science (1974) 185:862–864.
Greiner S, Wang X, Rauwolf U, Silber MV, Mayer K, Meurer J, Haberer G, Herrmann RG, Forthcoming. The complete nucleotide sequences of the five genetically distinct plastid genomes of Oenothera, subsection Oenothera: i. Sequence evaluation and plastome evolution. Nucleic Acids Res (2008) 36:2366–2378.
Hajdukiewicz PT, Allison LA, Maliga P. The two RNA polymerases encoded by the nuclear and the plastid compartments transcribe distinct groups of genes in tobacco plastids. EMBO J (1997) 16:4041–4048.[CrossRef][Web of Science][Medline]
Harte C. Oenothera: contributions of a plant to biology. In: Monographs on theoretical and applied genetics—Frankel R, Grossman M, Linskens HF, Maliga P, Riley R, eds. (1994) Berlin (Germany): Springer. 1–261.
Heribert-Nilsson N. Die Variabilität der Oenothera Lamarckiana und das Problem der Mutation. Z indukt Abstamm Vererbungsl (1912) 9:89–231.
Herrmann RG. Eukaryotism, towards a new interpretation. In: Eukaryotism and symbiosis—Schenk HEA, Herrmann RG, Jeon KW, Müller NE, Schwemmler W, eds. (1997) Berlin (Germany): Springer. 73–118.
Herrmann RG, Maier RM, Schmitz-Linneweber C. Eukaryotic genome evolution: rearrangement and coevolution of compartmentalized genetic information. Philos Trans R Soc Lond B Biol Sci (2003) 358:87–97.
Herrmann RG, Possingham JV. Plastid DNA—the plastome. In: Results and problems in cell differentiation—Reinert J, ed. (1980) Berlin (Germany): Springer. 45–96.
Herrmann RG, Seyer P, Schedel R, et al, (11 co-authors). The plastid chromosomes of several dicotyledons. In: Biological chemistry of organelle formation—Bücher T, Sebald W, Weiß H, eds. (1980) Berlin (Germany): Springer. 97–112.
Herrmann RG, Westhoff P. Thylakoid biogenesis: the result of a complex phylogenetic puzzle. In: Regulation of photosynthesis—Aro E-M, Andersson B, eds. (2001) Dordrecht (The Netherlands): Kluwer Academic Publishing. 1–28.
Homann A, Link G. DNA-binding and transcription characteristics of three cloned sigma factors from mustard (Sinapis alba L.) suggest overlapping and distinct roles in plastid gene expression. Eur J Biochem (2003) 270:1288–1300.[Web of Science][Medline]
Kanamaru K, Tanaka K. Roles of chloroplast RNA polymerase sigma factors in chloroplast development and stress response in higher plants. Biosci Biotechnol Biochem (2004) 68:2215–2223.[CrossRef][Medline]
Kapoor S, Sugiura M. Identification of two essential sequence elements in the nonconsensus type II PatpB-290 plastid promoter by using plastid transcription extracts from cultured tobacco BY-2 cells. Plant Cell (1999) 11:1799–1810.
Kapralov MV, Filatov DA. Widespread positive selection in the photosynthetic Rubisco enzyme. BMC Evol Biol (2007) 7:73.[CrossRef][Medline]
Klughammer C, Schreiber U. Measuring P700 absorbance changes in the near infrared spectral region with a dual wavelength pulse modulation system. In: Photosynthesis: mechanisms and effects—Garab G, ed. (1998) Dordrecht (The Netherlands): Kluwer Academic Publishers. 4357–4360.
Kofer W, Koop HU, Wanner G, Steinmüller K. Mutagenesis of the genes encoding subunits A, C, H, I, J and K of the plastid NAD(P)H-plastoquinone-oxidoreductase in tobacco by polyethylene glycol-mediated plastome transformation. Mol Gen Genet (1998) 258:166–173.[CrossRef][Web of Science][Medline]
Kooten O, Snel JFH. The use of chlorophyll fluorescence nomenclature in plant stress physiology. Photosynth Res (1990) 25:147–150.[CrossRef][Web of Science]
Levin DA. The cytoplasmic factor in plant speciation. Syst Bot (2003) 28:5–11.
Levin RA, Wagner WL, Hoch PC, Hahn WJ, Rodriguez A, Baum DA, Katinas L, Zimmer EA, Sytsma KJ. Paraphyly in tribe Onagreae: insights into phylogenetic relationships of Onagraceae based on nuclear and chloroplast sequence data. Syst Bot (2004) 29:147–164.[CrossRef]
Levin RA, Wagner WL, Hoch PC, Nepokroff M, Pires JC, Zimmer EA, Sytsma KJ. Family-level relationships of Onagraceae based on chloroplast rbcL and ndhF data. Am J Bot (2003) 90:107–115.
Liere K, Börner T. Transcription of plastid genes. In: Regulation of transcription in plants—Grasser KD, ed. (2006) Oxford: Backwell Publishing. 184–233.
Martin W. Gene transfer from organelles to the nucleus: frequent and in big chunks. Proc Natl Acad Sci USA (2003) 100:8612–8614.
Martin W, Stoebe B, Goremykin V, Hapsmann S, Hasegawa M, Kowallik KV. Gene transfer to the nucleus and the evolution of chloroplasts. Nature (1998) 393:162–165.[CrossRef][Web of Science][Medline]
Meurer J, Lezhneva L, Amann K, Gödel M, Bezhani S, Sherameti I, Oelmüller R. A peptide chain release factor 2 affects the stability of UGA-containing transcripts in Arabidopsis chloroplasts. Plant Cell (2002) 14:3255–3269.
Meurer J, Meierhoff K, Westhoff P. Isolation of high-chlorophyll-fluorescence mutants of Arabidopsis thaliana and their characterisation by spectroscopy, immunoblotting and northern hybridisation. Planta (1996) 198:385–396.[CrossRef][Web of Science][Medline]
Mrá
ek J, Greiner S, Cho WK, et al, (16 co-authors). Construction, database integration, and application of an Oenothera EST library. Genomics (2006) 88:372–380.[CrossRef][Web of Science][Medline]
Neckermann K, Zeltz P, Igloi GL, Kossel H, Maier RM. The role of RNA editing in conservation of start codons in chloroplast genomes. Gene (1994) 146:177–182.[CrossRef][Web of Science][Medline]
Raven PH. Onagraceae as a model of plant evolution. In: Plant evolutionary biology: a symposium honouring G. Ledyard Stebbins—Gottlieb LD, Jain SK, eds. (1988) London: Chapman and Hall. 85–107.
Renner O. Die pflanzlichen Plastiden als selbständige Elemente der genetischen Konstitution. Ber Math Phys Kl Sächs Akad Wiss Leipzig (1934) 86:241–266.
Schmitz-Linneweber C, Kushnir S, Babiychuk E, Poltnigg P, Herrmann RG, Maier RM. Pigment deficiency in nightshade/tobacco cybrids is caused by the failure to edit the plastid ATPase alpha-subunit mRNA. Plant Cell (2005) 17:1815–1828.
Schmitz-Linneweber C, Regel R, Du TG, Hupfer H, Herrmann RG, Maier RM. The plastid chromosome of Atropa belladonna and its comparison with that of Nicotiana tabacum: the role of RNA editing in generating divergence in the process of plant speciation. Mol Biol Evol (2002) 19:1602–1612.
Schmitz-Linneweber C, Tillich M, Herrmann RG, Maier RM. Heterologous, splicing-dependent RNA editing in chloroplasts: allotetraploidy provides trans-factors. EMBO J (2001) 20:4874–4883.[CrossRef][Web of Science][Medline]
Schreiber U, Bilger W, Hormann H, Neubauer C. Chlorophyll fluorescence as a diagnostic tool: basics and some aspects of practical relevance. In: Photosynthesis: a comprehensive treatise—Rughavendra AS, ed. (1998) Cambridge: Cambridge University Press. 320–336.
Schwenkert S, Umate P, Dal Bosco C, Volz S, Mlçochová L, Zoryan M, Eichacker LA, Ohad I, Herrmann RG, Meurer J. PsbI affects the stability, function, and phosphorylation patterns of photosystem II assemblies in tobacco. J Biol Chem (2006) 281:34227–34238.
Shiina T, Tsunoyama Y, Nakahira Y, Khan MS. Plastid RNA polymerases, promoters, and transcription regulators in higher plants. Int Rev Cytol (2005) 244:1–68.[CrossRef][Web of Science][Medline]
Shikanai T. RNA editing in plant organelles: machinery, physiological function and evolution. Cell Mol Life Sci (2006) 63:698–708.[CrossRef][Web of Science][Medline]
Silhavy D, Maliga P. Mapping of promoters for the nucleus-encoded plastid RNA polymerase (NEP) in the iojap maize mutant. Curr Genet (1998) 33:340–344.[CrossRef][Web of Science][Medline]
Sriraman P, Silhavy D, Maliga P. The phage-type PclpP-53 plastid promoter comprises sequences downstream of the transcription initiation site. Nucleic Acids Res (1998) 26:4874–4879.
Steiner EE, Stubbe W. A contribution to the population biology of Oenothera grandiflora L'Her. Am J Bot (1984) 71:1293–1301.[CrossRef][Web of Science]
Stubbe W. Extrem disharmonische Genom-Plastom-Kombinationen und väterliche Plastidenvererbung bei Oenothera. Z Vererbungsl (1963) 94:392–411.[CrossRef][Web of Science]
Stubbe W. The role of the plastome in evolution of the genus Oenothera. Genetica (1964) 35:28–33.[Medline]
Stubbe W. Oenothera—an ideal system for studying the interaction of genome and plastome. Plant Mol Biol Rep (1989) 7:245–257.[CrossRef]
Stubbe W, Pietsch B, Kowallik KV. Cytologische Untersuchungen über väterliche Plastidenvererbung und plastomabhängige Degradation der Samenanlagen bei einem Oenothera-Bastard. Biol Zentralbl (1978) 97:39–52.
Stubbe W, Raven PH. A genetic contribution to the taxonomy of Oenothera sect. Oenothera (including subsection Euoenothera, Emersonia, Raimannia and Munzia). Plant Syst Evol (1979) 133:39–59.[CrossRef]
Swiatecka-Hagenbruch M, Liere K, Börner T. High diversity of plastidial promoters in Arabidopsis thaliana. Mol Genet Genomics (2007) 277:725–734.[CrossRef][Web of Science][Medline]
Tillich M, Funk HT, Schmitz-Linneweber C, Poltnigg P, Sabater B, Martin M, Maier RM. Editing of plastid RNA in Arabidopsis thaliana ecotypes. Plant J (2005) 43:708–715.[CrossRef][Web of Science][Medline]
Tsudzuki T, Wakasugi T, Sugiura M. Comparative analysis of RNA editing sites in higher plant chloroplasts. J Mol Evol (2001) 53:327–332.[CrossRef][Web of Science][Medline]
Westhoff P. Transcription of the gene encoding the 51 kd chlorophyll a-apoprotein of the photosystem II reaction centre from spinach. Mol Gen Genet (1985) 201:115–123.[CrossRef][Web of Science]
Westhoff P, Herrmann RG. Complex RNA maturation in chloroplasts. The psbB operon from spinach. Eur J Biochem (1988) 171:551–564.[Web of Science][Medline]
Yang Z. PAML: a program package for phylogenetic analysis by maximum likelihood. Comput Appl Biosci (1997) 13:555–556.
![]()
CiteULike
Connotea
Del.icio.us What's this?
This article has been cited by other articles:
![]() |
U. Rauwolf, H. Golczyk, J. Meurer, R. G. Herrmann, and S. Greiner Molecular Marker Systems for Oenothera Genetics Genetics, November 1, 2008; 180(3): 1289 - 1306. [Abstract] [Full Text] [PDF] |
||||
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||





