Energy use efficiency is characterized by an ... - Semantic Scholar

1 downloads 0 Views 876KB Size Report
Nov 24, 2009 - aBayer BioScience N.V., 9052 Gent, Belgium; bDepartment of Plant Systems Biology, Flanders ..... vigour of oilseed rape lines and hybrids.
Energy use efficiency is characterized by an epigenetic component that can be directed through artificial selection to increase yield Miriam Haubena, Boris Haesendonckxa, Evi Standaerta, Katrien Van Der Kelenb,c, Abdelkrim Azmid, Herve´ Akpoa, Frank Van Breusegemb,c, Yves Guisezd, Marc Botsa, Bart Lamberta, Benjamin Lagaa, and Marc De Blocka,1 aBayer

BioScience N.V., 9052 Gent, Belgium; bDepartment of Plant Systems Biology, Flanders Institute for Biotechnology (VIB), 9052 Gent, Belgium; of Plant Biotechnology and Genetics, Ghent University, 9052 Gent, Belgium; and dLaboratory for Plant Biochemistry and Physiology, Department of Biology, University of Antwerp, 2020 Antwerpen, Belgium cDepartment

Edited by Marc C. E. Van Montagu, Ghent University, Ghent, Belgium, and approved September 29, 2009 (received for review August 4, 2009)

Brassica napus 兩 epigenome

T

he variability of a quantitative trait depends on genetic, epigenetic, and environmental factors as well as on their mutual interactions. Particularly, plants that cannot escape their environment have to adapt quickly. The importance of the epigenome in both short-term and long-term adaptations to stress was demonstrated by the elevated recombination frequencies in the offsprings of Arabidopsis thaliana populations that had been stressed with UV-B radiation or flagellin (1). These increased levels of recombination could be transmitted to successive generations that had not been exposed to the stress. Although there is no direct proof that recurrent selection can direct the epigenetic component of a quantitative trait, previous data suggest that it might be possible (2, 3). From a population obtained by a cross between two inbred flour beetles (Tribolium casteneum) with an average pupa weight of 2,400 ␮g, a population with an average pupa weight of 5,800 ␮g could be selected by continuous inbreeding for 120 generations (2, 3). Surprisingly, such a genetic variability was maintained over so many generations of inbreeding. Most probably, at a certain stage of the selection process, not genetic, but epigenetic variability, was dealt with, as further supported by identical experiments starting with one of the inbred lines homozygous for the genes affecting pupa weight (4). www.pnas.org兾cgi兾doi兾10.1073兾pnas.0908755106

For a long time, the epigenome has been believed to reset during meiosis. However, in numerous cases, epigenetic changes can be maintained during many generations. In plants, the flower-specific hypermethylated superman and agamous epialleles and the disease resistance hypomethylated late-flowering and ball epialleles are well-known (5–7). In fact, the transgenerational epigenetic inheritance holds true for many taxa, including animals, fungi, and bacteria (8). Recently, in Arabidopsis, epialleles across the genome were shown to be stably transmitted over consecutive generations (9). One of the most complex quantitative traits in plants is yield because it is determined by many interacting components. Although a major subject in plant research for many decennia, the understanding of the biological network controlling yield is only poorly understood. Here we show that in canola (Brassica napus), energy use efficiency (EUE) is a major yield component with an epigenetic feature that can be directed by recursive selection. The shaping of the epigenome allows the increase of seed yield in both open-pollinating lines and hybrids. Results and Discussion Selection of Sublines from an Isogenic Line Based on Respiration. A

positive correlation between yield and low cellular respiration rates has been demonstrated in canola, Lolium perenne, tomato (Solanum esculentum), and cucumber (Cucumis sativus) (10–16). Based on the published experimental concept (2–4), an artificial bidirectional selection procedure was established in canola for high and low cellular respiration rates. We started from an isogenic doubled haploid line to eliminate variation due to differences at the genetic level. As schematically presented (Fig. 1), out of 200 seedlings from the doubled haploid population ‘‘Simon,’’ seedlings with the lowest and highest cellular respiration were identified and retained. The two plants with the highest and lowest respiration were withheld for seed production by self-fertilization. These two populations were the starting material for four additional rounds of selection for sublines with high and low respiration rates. Seven lines were generated: four lines with lower (LR76, LR76-crashed, LR77, and LR79) and three with higher (HR80, HR82, and HR83) respiration than that of the control from which the selection had been initiated (Fig. 2). The ‘‘LR76-crashed’’ line was derived from a further Author contributions: M.D.B. designed research; M.H., B.H., E.S., K.V.D.K., A.A., and H.A. performed research; F.V.B., Y.G., M.B., B. Lambert, and B. Laga analyzed data; and M.D.B. wrote the paper. The authors declare no conflict of interest. This article is a PNAS Direct Submission. Freely available online through the PNAS open access option. 1To whom correspondence should be addressed. E-mail: [email protected].

This article contains supporting information online at www.pnas.org/cgi/content/full/ 0908755106/DCSupplemental.

PNAS 兩 November 24, 2009 兩 vol. 106 兩 no. 47 兩 20109 –20114

PLANT BIOLOGY

Quantitative traits, such as size and weight in animals and seed yield in plants, are distributed normally, even within a population of genetically identical individuals. For example, in plants, various factors, such as local soil quality, microclimate, and sowing depth, affect growth differences among individual plants of isogenic populations. Besides these physical factors, also epigenetic components contribute to differences in growth and yield. The network that regulates crop yield is still not well understood. Although this network is expected to have epigenetic elements, it is completely unclear whether it would be possible to shape the epigenome to increase crop yield. Here we show that energy use efficiency is an important factor in determining seed yield in canola (Brassica napus) and that it can be selected artificially through an epigenetic feature. From an isogenic canola population of which the individual plants and their self-fertilized progenies were recursively selected for respiration intensity, populations with distinct physiological and agronomical characteristics could be generated. These populations were found to be genetically identical, but epigenetically different. Furthermore, both the DNA methylation patterns as well as the agronomical and physiological characteristics of the selected lines were heritable. Hybrids derived from parent lines selected for high energy use efficiencies had a 5% yield increase on top of heterosis. Our results demonstrate that artificial selection allows the increase of the yield potential by selecting populations with particular epigenomic states.

The Selected Sublines Have a Distinct EUE. High cellular respiration

was expected to be linked to high energy production and vice versa. However, when the total amount of NAD(P)H, reflecting the energy content, was measured (Table 1), lines with a high and low respiration had a low and high NAD(P)H content, respectively, except for line ‘‘LR76-crashed.’’ When the ratio NAD(P)H versus respiration, defined as EUE, was calculated, respiration and EUE correlated inversely, except for line LR76crashed (Table 1). Probably, the respiration in line LR76-crashed dropped below a critical threshold to maintain sufficient energy production. Complex I Activity and Ascorbate Content Are High and Photorespiration Is Low in the Low-Respiration Sublines. During the selection

Fig. 1. Artificial selection for respiration and EUE in canola. The selection was initiated from an isogenic doubled haploid line. Approximately five selected plants with the lowest and highest respiration were self-fertilized, and the progenies were tested for respiration and EUE. The next rounds of selections were done in one direction for low or high respiration with the two lines with the lowest and highest respiration, respectively. Three to five rounds of selection were sufficient to generate lines with distinct respiration and EUE.

selection for low respiration starting from the line LR76. The respiration of line LR76-crashed dropped below a critical threshold, resulting in reduced vigor (see below). At the individual plant level, the variation was high in the high-respiration lines [for instance, line HR82 with 125 ⫾ 74% (mean ⫾ standard deviation) compared to the control], and minimal in the lowrespiration lines (for instance, line LR76 with 90 ⫾ 18% versus the control). The high variation in respiration found between individual plants in the high-respiration lines is mainly due to the exponential increase of respiration once a certain threshold is passed (1.5 times the average respiration of the control line). This results in plants having up to six times higher respiration than the average respiration of the control line (supporting information (SI) Fig. S1). The high efficiency with which these lines were generated (from only 200 seedlings per selection cycle) suggests that no mutations are involved. Extensive amplified fragment-length polymorphism (AFLP) analyses did not identify differences between these lines and the control, except for the very high-respiration line HR83. The highest respiring plants of line HR83 have an abnormal development, are both male and female sterile, and have often chromosome deletions and inversions. Therefore, HR83 was not retained. Also the line LR79 was omitted in further analyses because its respiration was close to that of the control (97%). 20110 兩 www.pnas.org兾cgi兾doi兾10.1073兾pnas.0908755106

procedure, respiration was quantified by measuring the reduction of 2,3,5-triphenyltetrazolium chloride (TTC) of hypocotyl explants. TTC reduction occurs at the end of the mitochondrial respiratory chain at complex IV (17) and, therefore, reflects the total electron flow through the mitochondrial respiratory chain, including the alternative oxidative respiratory pathway. The inhibitor of the alternative oxidase SHAM inhibits the reduction of TTC. The electrons enter the mitochondrial electron transport chain through complex I, complex II, and the internal and external alternative NAD(P)H dehydrogenases. The activity of complex I, the main dehydrogenase, was measured in leaf 3 of control and selected lines (Table 1), and an inverse correlation with respiration and complex I activity was observed. Recently, l-galactono-1,4-lactone dehydrogenase, the last enzyme of the plant ascorbate biosynthesis pathway, has been found to associate with an 800-kDa subcomplex of complex I and to have an important function in the accumulation of complex I, because complex I does not accumulate in Arabidopsis null mutants (18, 19). This apparent interrelation between ascorbate biosynthesis and complex I prompted us to measure the ascorbate content in the leaves of the different lines (Table 1). The lines with the highest respiration and the lowest complex I activity had the lowest ascorbate content and inversely. In summary, these results suggest that lines with a high and low respiratory rate have a reduced and increased complex I content, respectively. The high-respiration lines resemble the cytoplasmic male-sterile II (cmsII) mutant of Nicotiana sylvestris that lacks the NAD7 subunit of complex I (20–21). The csmII mutant has no complex I activity, but its respiration is still high due to the compensatory enhanced activity of complex IV (cytochrome oxidase) and, probably, an increased activity of the alternative NAD(P)H dehydrogenases to meet the ATP needs of the cells. Lines with a high complex I content are more efficient in producing ATP, resulting in a low respiratory rate. In tobacco (Nicotiana tabacum) and cucumber (Cucumus sativus), the decreased complex I capacity resulted in an increase in photorespiration (16, 22). Similarly, we found that the highrespiration lines had a significantly increased photorespiration, while that of the low-respiration lines was strongly reduced (Table 1), implying that low respiration, but not too low (as in the LR76-crashed line), correlates with a favorable physiological state and high respiration with a less fitness. Yield Is Correlated with EUE. The specific characteristics of the

physiological states hint at a high or low stress tolerance (Fig. S2) and field performance. The selected and control lines were tested for seed yield in field trials (five locations with three to six replications of 10 m2/line) for 3 years (Table 1). Low-respiration lines with a high EUE had on average up to 8% higher seed yield than that of the control, while on average the seed yield was reduced by up to 10% in high-respiration lines with a low EUE. In fields with moderate drought stress, the line with the highest EUE had a 20% higher yield than that of the control, while the seed yield of the line with the highest respiration and lowest EUE Hauben et al.

Fig. 2. Respiration of six selected populations versus the original starting population. The values were percent-normalized versus the average of the starting line ‘‘Simon.’’ Each node represents an individual plant. Each population contains approximately 200 plants. The bar in the boxes represents the median of the population and the other bars, the standard deviation with indication of the mean.

dropped by 20%. With a Pearson correlation of 0.96, EUE and yield were very strongly positively correlated. The Selected Sublines Have Altered Global DNA Methylation and Histone Modification. The fact that the lines with high and low

respiration had been generated efficiently from an isogenic line and that, based on the AFLP results, the selected lines were indistinguishable, suggest that the distinct physiological characteristics of the lines have an epigenetic basis. Hence, differential cytosine methylation and overall histone acetylation and methylation were studied to explore the epigenetic changes at the molecular level. First, a methylation-sensitive AFLP (ms-AFLP) analysis with HpaII and MspI as methylation-sensitive restriction enzymes on genomic DNA of leaf 3 revealed distinct differences (⬇0.3%) between the selected lines and the control ‘‘Simon’’ (Fig. S3). Remarkably, all differentially methylated fragments that were sequenced (110 in total) corresponded to coding sequences. The ms-AFLP patterns were specific for each line and were maintained throughout different self-fertilizations for at least eight generations. Secondly, to evaluate differences in global methylation, cytosine extension assays were carried out on genomic DNA of leaf 3 in the different lines. In this assay, by means of the methylation-sensitive HpaII and MspI enzymes, it is possible to incorporate a single [3H]-labeled dCTP, after which the incorporated dCTP activity can be measured. An increase in

activity corresponds to a decrease in global methylation. In general, the genomic DNA of the selected lines was hypomethylated (Fig. 3B), but the hypomethylation was more prominent in the low-respiration lines, with line LR76-crashed, being the most hypomethylated. The distinct physiological state of the lines was also reflected at the histone level of the chromatin (Fig. 3A). In general, in the selected lines, except for HR82, the methylation at lysine 9 (Lys-9) and Lys-27 of histone 3 and acetylation at Lys-8, Lys-12, and Lys-16 of histone 4 was lower than that of the control line. The low-respiration line LR76 and its progeny LR76-crashed had a similarly low histone methylation and acetylation, whereas in the high-respiration line HR80 and line HR82, which is less vigorous than line HR80, the histone methylation and acetylation was the lowest and similar to that of the control, respectively. Hence, each of the selected lines is characterized by a specific degree of methylation and acetylation of histones 3 and 4. However, this overall picture does not reflect the physiological properties of the lines. More in-depth studies, such as chromatin immunoprecipitation (ChIP) analyses, might reveal a possible relation between histone modifications, gene expression, and vigor. Artificial selection for altered respiration and EUE clearly results in overall changes in DNA methylation and histone modifications through the plasticity of the epigenome (23). Extensive ms-AFLP analyses showed that none of the DNA

Table 1. Physiological and agronomical properties of selected canola lines (% normalized versus control†) Seed yield

LR76-crashed LR76 LR77 HR80 HR82

Respiration

NAD(P)H content

EUE‡

Complex I activity

Ascorbate content

Photorespiration

Average§

Field with stress¶

Optimal field储

74.5*** 90*** 92.5*** 125*** 139***

73** 98 103* 99 98

98 109 111 79 71

140** 127** 120** 100 95*

200*** 208*** 185*** 79* 74**

80** 81** 91* 105* 116**

97 104* 108** 92** 90**

106 114** 120** 81** 80**

97 97 103 93* 88** PLANT BIOLOGY

Line

†Original

line (variety ‘‘Simon’’) from which the selected lines were derived. use efficiency ⫽ NAD(P)H content/respiration. §Five locations, six repetitions/location, 3 years. ¶One location with sandy soil and moderate drought stress, 2 years. 储One location with loamy soil, 3 years. *, P ⬍ 0.05; **, P ⬍ 0.01; and ***, P ⬍ 0.001 versus control. ‡Energy

Hauben et al.

PNAS 兩 November 24, 2009 兩 vol. 106 兩 no. 47 兩 20111

Arabidopsis (1), the epigenetic respiration component of canola can be transmitted in reciprocal backcrosses with a nonselected control line. For example, in the high-respiration line HR82 and similarly in the low-respiration line LR77, the self-fertilized progeny and the backcrosses to the original ‘‘Simon’’ line with HR82 and LR77 as females and males resulted in progenies with 127%, 124%, and 128% respiration and with 89%, 91%, and 93% respiration versus the control, respectively. The epigenetic inheritance was also reflected in the largely conserved DNA methylation patterns in the crosses compared to those of the self-fertilized selected lines. The methylation sites that were not conserved had a strong tendency toward hypomethylation, independently from the direction of the cross. The Epigenetic EUE Component Can Be Added on Top of Heterosis.

Although ‘‘Simon’’ is a rather good performing canola line, we wanted to investigate whether the epigenetic EUE component could increase yield in elite varieties and hybrids. To this end, we selected high EUE in the parental lines of a top-yielding canola hybrid. Both the female and male lines are doubled haploids, and for this, each parental line was isogenic. After three to four recurrent selections for high EUE, hybrids were made with the selected parental lines (Fig. S4). The yield in these hybrids was higher (102–106%) than that of the control hybrid of which the parental lines had not been selected. The physiological characteristics were similar to the high-EUE ‘‘Simon’’ sublines: decreased mitochondrial complex I activity, reduced photorespiration, and increased ascorbate, ATP, and NAD(P)H contents (Table S1). Also at the molecular level, the selected hybrids were distinct from the control. Besides the specific differences in global DNA methylation with hypermethylation and hypomethylation in the cotyledons and leaves of the selected hybrids, respectively (Fig. 3C), the transcription profiles as analyzed by cDNA-AFLP were very distinct from the control hybrid (Fig. S5).

Fig. 3. Characterization of the histone modifications and quantification of the global methylation of the genomic DNA by the cytosine extension assay. (A) Histones prepared from leaf 3 of the control line ‘‘Simon’’ and selected sublines. Below the immunoblots, the corresponding SYPRO Ruby-stained gel is shown to demonstrate the equal loading. (B) Cytosine extension assay using genomic DNA prepared from leaf 3 of ‘‘Simon’’ (control) and selected sublines. Leaves were harvested from the same population as used for studying the histone modifications. (C) Cytosine extension assay with genomic DNA prepared from cotyledons and the fourth leaf of a nonselected (control) and a selected hybrid. The cytosine extension assays were done on 10 individual plants. The error bars indicate the standard error.

methylation patterns found in the selected lines corresponded with methylation patterns of individual plants from the original population, indicating that recurrent selection for altered respiration at the individual plant level allows the selection for genes with specific and stable methylation states and distinct phenotypes, as a consequence. The Subline Specific Respiration Rate Is Transmitted in Self-Fertilized Seeds and Backcrosses. From the approximately eight generations of

seed upscalings in the greenhouse and 3 years of elaborate field trials, it became obvious that the epigenetic EUE component of the lines was stably inherited by self-fertilized siblings (Fig. S3). As already described for the epigenetic recombination component in 20112 兩 www.pnas.org兾cgi兾doi兾10.1073兾pnas.0908755106

Concluding Remarks. Our results show that EUE is a distinct feature of plant vigor and yield and that it possesses an epigenetic component that can be directed by artificial selection. In other words, it is feasible to select for a complex trait solely based on the epigenetic component. The main difference with the selection used in classical breeding schemes is that it is not only done at the population, but primary, at the plant level and in a recursive manner. However, the underlying mechanism that stabilizes the distinct epigenetic states is not understood. Starting from an isogenic population composed of plants with a variable epigenome, subpopulations with at least a partially fixed epigenome can be selected rather efficiently. It becomes even more puzzling when, during plant development, the epigenome of the selected sublines changes in opposite ways versus the control. For example, in hybrids generated by crossing parental lines that have been selected for higher EUE, the total genomic DNA is hypermethylated in the cotyledons, while in the fourth leaf, the genomic DNA is hypomethylated versus the control hybrid (Fig. 3C). It has to be mentioned that we found in our selections for higher EUE subpopulations that in the first rounds of selection returned to their original state. Only by recurrent selection the higher EUE could be fixed. It is not unlikely that the epigenetic stabilization can be tracked down to the expression of specific miRNAs. We are testing this possibility by making an inventory of the small RNAs in the selected and control lines. The ‘‘epigenetic EUE component’’ is stably inherited, allowing the creation of distinct isogenic sublines that can be used in breeding. Probably as importantly, these isogenic sublines with distinct yield characteristics have a low ‘‘genetic noise’’ and, therefore, are excellent material to study yield and EUE at the molecular level. Initial experiments with cDNA-AFLP show clear differences in transcription between the sublines (Fig. S5) Hauben et al.

Materials and Methods Artificial Selection for Respiration and EUE. Seedlings were grown in vitro for 2 weeks on agar in half-strength Murashige and Skoog medium supplemented with 2% sucrose. The shoot tips of the seedlings were placed on the medium for rooting, while five hypocotyl explants per seedling were cultured for 5 days on callus-inducing medium (Murashige and Skoog medium supplemented with 3% sucrose and 1 mg/L 2,4-dichlorophenoxyacetic acid, 0.25 mg/L ␣-naphthaleneacetic acid, and 1 mg/L 6-benzylaminopurine). Cellular respiration of the hypocotyl explants was measured. The rooted shoot tips of approximately five seedlings with the highest and lowest respiration were transferred to the greenhouse for seed production by self-fertilization. Both respiration and NAD(P)H content of approximately 35– 40 seedlings of the obtained progenies were measured. Lines with the lowest and highest respiration and highest and lowest EUE, respectively, were retained. The next rounds of selections were done in one direction for low or high respiration in lines with the lowest and highest respiration, respectively (Fig. 1). Three to five rounds of selection were sufficient to generate lines with distinct respiration and EUE. The same protocol was used for the production of more EUE hybrids. The male-sterile female, the maintainer, and the male restorer lines were each selected for low respiration and high EUE. After three cycles, the male-sterile female and restorer were crossed to produce hybrid seeds (Fig. S4).

1. Molinier J, Ries G, Zipfel C, Hohn B (2006) Transgeneration memory of stress in plants. Nature 442:1046 –1049. 2. Enfield FD, Comstock RE, Braskerud O (1966) Selection for pupa weight in Tribolium castaneum. I. Parameters in base populations. Genetics 54:523–533. 3. Kaufman PK, Enfield FD, Comstock RE (1977) Stabilizing selection for pupa weight in Tribolium castaneum. Genetics 87:324 –341. 4. Enfield FD, Braskerud O (1989) Mutational variance for pupa weight in Tribolium castaneum. Theor Appl Genet 77:416 – 420. 5. Jacobsen SE, Meyerowitz EM (1997) Hypermethylated SUPERMAN epigenetic alleles in Arabidopsis. Science 277:1100 –1103. 6. Jacobsen SE, Sakai H, Finnegan EJ, Cao X, Meyerowitz EM (2000) Ectopic hypermethylation of flower-specific genes in Arabidopsis. Curr Biol 10:179 –186. 7. Stokes TL, Kunkel BN, Richards EJ (2002) Epigenetic variation in Arabidopsis disease resistance. Genes Dev 16:171–182. 8. Jablonka E, Raz G (2009) Transgenerational epigenetic inheritance: Prevalence, mechanisms, and implications for the study of heredity and evolution. Q Rev Biol 84:131–176. 9. Johannes F, et al. (2009) Assessing the impact of transgenerational epigenetic variation on complex traits. PLoS Genet 5:e1000530. 10. Wilson D (1982) Response to selection of dark respiration rate of mature leaves in Lolium perenne and its effects on growth of young plants and simulated swards. Ann Bot 49:303–312. 11. Wilson D, Jones JG (1982) Effect of selection for dark respiration rate of mature leaves on crop yields of Lolium perenne cv. S23. Ann Bot 49:313–320. 12. Kraus E, Aydemir Y, Duin S, Kollo¨ffel C, Lambers H (1993) Yield advantage of a ‘slow-’over a ‘fast-’ respiring population of Lolium perenne cv. S23 depends on plant density. New Phytol 123:39 – 44. 13. Kraus E, Lambers H, Kollo¨ffel C (1993) The effect of handling on the yield of two populations of Lolium perenne selected for differences in mature leaf respiration rate. Physiol Plant 89:341–346.

Hauben et al.

Physiological and Biochemical Assays. Cellular respiration of hypocotyl explants was quantified by measuring the reduction of TTC as described (14). ATP and total NAD(P)H contents were quantified as described (24, 25). EUE was expressed as the ratio of the percent-normalized values versus the control of NAD(P)H content to the amount of reduced TTC. Ascorbate content and complex I activity were measured with the reflectometric ascorbic acid test (Merck) and the MitoProfile Dipstick Assay kit (MitoSciences), respectively. The photorespiration intensity was quantified in vitro by floating cotyledons or leaf pieces on ammonium-free medium containing 3–5 mg/L glufosinate, incubated for 24 h under continuous light (⬇70 ␮mol m⫺2 s⫺1). The ammonium production was evaluated as described (26). Photorespiration was expressed as percentage produced ammonium versus the control. Field Trials. For the field trials, five locations with different (from sandy to loamy) soil qualities were used. Each location was subdivided in six plots per line. Each plot was 10 m2 in size, with different arrangements of the plots for each location. Yield was expressed as kg seed/plot. Statistics. Data were statistically analyzed by one-way ANOVA with Dunnett’s posttest with the Prism version 5.00 (GraphPad Software). Molecular Analyses. Methylation-sensitive AFLP (27) was done with the restriction enzyme combinations EcoRI-MspI/HpaII. The methylation state of the total genomic DNA was quantified with the cytosine extension assay as described (28). In summary, genomic DNA was prepared from cotyledons, leaf 3, or leaf 4 and assayed with the Quant-iT dsDNA High-Sensitivity Assay kit (Invitrogen). Of 10 individual leaves per line, 150 ng were used for the cytosine extension assay and digested with a 10-fold excess of HpaII and MspI. After 24 h, again a 10-fold excess of the restriction enzymes was added for another 24-h incubation. After 48 h, the single-nucleotide extension reaction was done with the digested genomic DNA. Of the reaction mixture, 10 ␮L were spotted on DE81 Whatman filters and washed three times with phosphate buffer (0.5 M, pH 7.0). The radioactive incorporation was measured (cpm). Analyses of Histone Modifications. Histones were acid-extracted, separated in NuPAGE 4 –12% Bis-Tris gel (Invitrogen), and subjected to a standard procedure for immunoblotting with anti-H4-acetyl-Lys-8, anti-H4-acetyl-Lys-12 antibodies (Millipore), and with anti-H3-dimethyl-Lys-9 and anti-H3-trimethylLys-27 antibodies (kind gifts of Thomas Jenuwein, Max Planck-Institute of Immunobiology, Freiburg, Germany). ACKNOWLEDGMENTS. We thank T. Jenuwein (Max Planck-Institute of Immunobiology, Freiburg, Germany) for providing the antibodies against anti-H3dimethyl-Lys-9 and anti-H3-trimethyl-Lys-27 and Dr. Martine De Cock for help in preparing the manuscript. This work was supported by the Institute for the Promotion of Innovation by Science and Technology in Flanders, Belgium (Research and Development project ‘‘Phoenix’’) and an Industrial Agreement for Training through Research (CIFRE) fellowship (to H.A.).

14. De Block M, De Brouwer D (2002) A simple and robust in vitro assay to quantify the vigour of oilseed rape lines and hybrids. Plant Physiol Biochem 40:845– 852. 15. Nunes-Nesi A, et al. (2005) Enhanced photosynthetic performance and growth as a consequence of decreasing mitochondrial malate dehydrogenase activity in transgenic tomato plants. Plant Physiol 137:611– 622. 16. Juszczuk IM, et al. (2007) Effect of mitochondrial genome rearrangement on respiratory activity, photosynthesis, photorespiration and energy status of MSC16 cucumber (Cucumis sativus) mutant. Physiol Plant 131:527–541. 17. Rich PR, Mischis LA, Purton S, Wiskich JT (2001) The sites of interaction of triphenyltetrazolium chloride with mitochondrial respiratory chains. FEMS Microbiol Lett 202:181–187. 18. Millar AH, et al. (2003) Control of ascorbate synthesis by respiration and its implications for stress responses. Plant Physiol 133:443– 447. 19. Pineau B, Layoune O, Danon A, De Paepe R (2008) l-Galactono-1,4-lactone dehydrogenase is required for the accumulation of plant respiratory complex I. J Biol Chem 283:32500 –32505. 20. Gutierres S, et al. (1997) Lack of mitochondrial and nuclear-encoded subunits of complex I and alteration of the respiratory chain in Nicotiana sylvestris mitochondrial deletion mutants. Proc Natl Acad Sci USA 94:3436 –3441. 21. Sabar M, De Paepe R, de Kouchkovsky Y (2000) Complex I impairment, respiratory compensations, and photosynthetic decrease in nuclear and mitochondrial male sterile mutants of Nicotiana sylvestris. Plant Physiol 124:1239 –1249. 22. Priault P, et al. (2006) The lack of mitochondrial complex I in a CMSII mutant of Nicotiana sylvestris increases photorespiration through an increased internal resistance to CO2 diffusion. J Exp Bot 57:3195–3207. 23. Tsaftaris AS, Polidoros AN, Kapazoglou A, Tani E, Kovacevic NM (2008) in Plant Breeding Reviews, eds Janick J (Wiley, New York, NY), vol 30, pp 49 –178. 24. Rawyler A, Pavelic D, Gianinazzi C, Oberson J, Braendle R (1999) Membrane lipid integrity relies on a threshold of ATP production rate in potato cell cultures submitted to anoxia. Plant Physiol 120:293–300.

PNAS 兩 November 24, 2009 兩 vol. 106 兩 no. 47 兩 20113

PLANT BIOLOGY

and microarray experiments being carried out to obtain a broader picture. Making an inventory of the small RNAs of the sublines and ChIP experiments with differentially methylated and/or expressed genes will contribute to a better understanding of the relation between DNA methylation (as the discrepancy in DNA methylation between cotyledons and leaves in the hybrids), histone modifications, transcription, and the resulting phenotypes. The fact that the concept can be used to further improve yield in already top-performing hybrids open perspectives to go beyond heterosis. The central role of plant mitochondria in energy metabolism and their involvement in biochemical pathways, such as photorespiration, amino acid biosynthesis, and maintenance of redox status, is conserved in most plant species. Indeed, ongoing experiments indicate that analogous selections for EUE can be applied in rice (Oryza sativa) and tomato (Solanum esculentum). The implementation of the ‘‘epigenetic’’ selection for EUE has the potential to increase yield in many crops and this will further contribute to a better knowledge of the epigenetic mechanism, especially in crops, for which the genomics tools are well developed.

25. Nakamura J, et al. (2003) Quantitation of intracellular NAD(P)H can monitor an imbalance of DNA single strand break repair in base excision repair deficient cells in real time. Nucleic Acids Res 31:e104. 26. De Block M, De Sonville A, De Brouwer D (1995) The selection mechanism of phosphinothricin is influenced by the metabolic status of the tissue. Planta 197:619 – 626.

20114 兩 www.pnas.org兾cgi兾doi兾10.1073兾pnas.0908755106

27. Matthes M, Singh R, Cheah S-C, Karp A (2001) Variation in oil palm (Elaeis guineensis Jacq.) tissue culture-derived regenerants revealed by AFLPs with methylation-sensitive enzymes. Theor Appl Genet 102:971–979. 28. Boyko A, et al. (2007) Transgenerational changes in the genome stability and methylation in pathogen-infected plants (virus-induced plant genome instability). Nucleic Acids Res 35:1714 –1725.

Hauben et al.