- Copyright © 2002 American Society of Plant Physiologists
Abstract
Recombinant inbred lines (RILs) derived from B73 × M017 were screened for cold germination (CG) and desiccation tolerance (DT) phenotypes. Reciprocal F1 hybrids were made between divergent RILs, and hybrids that showed differential phenotypes (parent-of-origin effect) for CG or DT were selected for profiling mRNA and protein expression. mRNA and proteins were extracted from embryo axes of seed germinated for 11 d at 12.5°C in the dark and developing embryos at 40% seed moisture (R5 stage) for CG and DT, respectively. GeneCalling analysis, an open-ended mRNA profiling method, identified 336 of 32,496 and 656 of 32,940 cDNA fragments that showed ≥1.5-fold change in expression between the reciprocal F1 hybrids for CG and DT, respectively. Protein expression map (PEM) analysis, an open-ended two-dimensional polyacrylamide gel electrophoresis, identified 117 of 2,641 and 205 of 1,876 detected proteins to be differentially expressed with ≥1.5-fold change between the reciprocal F1 hybrids in CG and DT samples, respectively. A subset of these proteins was identified by tandem mass spectrometry followed by database query of the spectra. The differentially expressed genes/proteins were classified into various functional groups including carbohydrate and amino acid metabolism, ion transporters, stress and defense response, polyamine metabolism, chaperonins, cytoskeleton associated, etc. Phenotypic analysis of seed from self-pollinated ears of the reciprocal F1 hybrids displayed small differences compared with the reciprocal hybrids themselves, suggesting a negligible effect of cytoplasmic factors on CG and DT traits. The results provide leads to improving our understanding of the genes involved in stress response during seed maturation and germination.
There have been several reports where differential phenotypic expression was observed between reciprocal F1 hybrids in maize (Zea mays) for various kernel and germination traits such as whole-kernel growth rate (Groszmann and Sprague, 1948), embryo and endosperm dry weight (Bagnara and Daynard, 1983), embryo oil and protein (Miller and Brimhall, 1951), zein synthesis (Chaudhuri and Messing, 1994), seed germination at low temperatures (Pinnell, 1949), and tolerance to drying injury (Bdliya and Burris, 1988). These differential expressions can be attributed to epigenetic phenomena such as genomic imprinting (described below) and xenia (refers to the effect of pollen on the endosperm phenotype in the same generation), dosage effects (in case of triploid tissue such as endosperm), and cytoplasmic effects (so-called maternal effects from mitochondrial and chloroplast genomes or nuclear-encoded cytoplasmically stored gene products).
Genomic imprinting, also referred to as uniparental dominance or parent-of-origin effect, is a well-studied phenomenon in both mammals and plants (for reviews, see Alleman and Doctor, 2000; Reik and Walter, 2001). Imprinting refers to a reversible epigenetic modification of loci resulting in differential expression of genes depending on the parent of origin. Unlike Mendelian inheritance, a given allele of an imprinted gene is phenotypically expressed in the F1 depending on whether it is transmitted through the male or female parent. Imprinting in mammals occurs primarily via differential DNA methylation and alteration in chromatin structure resulting in silencing of specific allele (see Feil and Khosla, 1999;Tilghman, 1999). Possible roles of differential methylation in imprinting of zein and α-tubulin genes in maize (Lund et al., 1995a,1995b) and chromatin remodeling (structural alteration in chromatin) in regulating imprinting of MEDEA (MEA) locus in Arabidopsis (Jeddeloh et al., 1999; Vielle-Calzada et al., 1999) have been reported. However, the specific role of DNA methylation-mediated gene silencing in imprinting of genes during embryo and endosperm development in plants is not clear (Vielle-Calzada et al., 1999; Luo et al., 2000; Russinova and de Vries, 2000).
Comparison of reciprocal F1 hybrids provides an excellent system to study genes associated with parent-of-origin effects (some of which may be due to genomic imprinting) in diploid tissue that influences various traits of agronomic interest. In hybrid seed production systems such as maize, sorghum, and millets, identification of genes associated with imprinting provides an opportunity to select parents (male or female) through which to introduce the desired alleles.
Northern corn-growing regions have a short crop-growing season. This often requires planting of seed early in the season, when soil temperatures are low, followed by early season harvesting to avoid frost injury. Therefore, it is important for grain to be used as seed to tolerate artificial drying (desiccation) when harvested at high seed moisture (≥40%) and to retain the ability to germinate under cold conditions for better field stand establishment. Extensive efforts have been devoted to understanding the genetics and physiology of the cold germination (CG) and desiccation tolerance (DT) traits (Burris, 2000;Hoegemeyer and Gutormsen, 2000). Changes in various compositional (oligosaccharides, lipids, hormones, enzymes, free radicals, and metabolic intermediates) and organellar/suborganellar properties (mitochondria, oil and protein bodies, starch grains, membrane integrity, and structure of chromatin and proteins) in the seed and germinating seedlings of maize have been reported to be associated with CG (Yacoob and Filion, 1986; Janowiak and Markowski, 1987; Stewart et al., 1990a, 1990b; Schell et al., 1991; Prasad, 1996, 1997; Santis et al., 1999) and DT (Leprince et al., 1993, 1994, 1995; Bochicchio et al., 1994a, 1994b; Oishi and Bewley, 1990; Chen and Burris, 1991;Brenac et al., 1997; Obendorf, 1997; Perdomo and Burris, 1998; Wolkers et al., 1998).
Much of the literature relevant to cold stress response at the molecular level comes from Arabidopsis where several cold-regulated (COR) genes induced in response to abscisic acid (ABA) were characterized (for review, see Thomashow, 1999). An ABA-independent low-temperature responsive regulatory cascade involving a sequential expression of C-repeat (CRT)/dehydration responsive element (DRE) binding factors genes, followed by CRT/DRE-containing COR genes was recently characterized in Arabidopsis (Stockinger et al., 1997; Gilmour et al., 1998; Liu et al., 1998; Shinozaki and Yamaguchi-Shinozaki, 2000). Overexpression of CRT/DRE binding factor-1 or DREB1A under the control of cauliflower mosaic virus 35S or stress inducible rd29A promoters was found to induce COR genes enhancing freezing tolerance in Arabidopsis (Jaglo-Ottosen et al., 1998; Kasuga et al., 1999). Seki et al. (2001) reported a genome-wide survey of mRNA expression pattern of approximately 1,300 genes, using cDNA microarray, in Arabidopsis subjected to drought and cold stress at various stages from germination to maturity. They identified 44 and 19 cDNAs that were inducible by drought and cold stress, respectively, of which 12 were DREB1A target genes. Such detailed molecular studies of cold and desiccation stress responses have not been reported in maize.
The objective of this study was to identify differentially expressed genes related to CG and DT phenotypes. We have identified phenotypically divergent recombinant inbred lines (RILs) from the cross B73 × M017 (SX-19) for CG and DT and produced reciprocal F1 hybrids. The reciprocal F1 hybrids that showed differential phenotypes for CG and DT, an indication of potential imprinting effects (or nuclear cytoplasmic effects), were profiled for mRNA and protein expression using open-ended platforms. This approach allowed the profiling of genetically identical material thus enhancing the effectiveness of expression profiling.
RESULTS
Screening, Crossing, and Phenotypic Analysis
Five RILs that showed high germination (>90%) and five RILs that showed low germination (≤30%) under cold test conditions were identified. Similarly, five each of tolerant and sensitive RILs to high-temperature desiccation were identified. All genotypes germinated normally at ambient temperatures (data not shown). Reciprocal F1 hybrids produced in a five-by-five diallel-crossing scheme, showed varying degrees of responses within each trait (data not shown). The F1 hybrids displayed transgressive (heterosis) and additive responses compared with their self-pollinated parental phenotypes. In instances where the reciprocal F1 hybrids showed differential phenotypes, the phenotype of the F1s with the high parent (i.e. parent with high germination score) as the maternal parent was typically higher in their germination scores compared with those with low parents as the maternal parent. Maternal dominance was clearly evident among some, but not all, RIL combinations for both traits. Of all of the possible RIL combinations tested within each trait, at least six and five showed parent-of-origin effect for CG and GT, respectively (data not shown). Because of the complexity and resource requirements of the profiling methods used, only one representative RIL combination from each trait was selected for expression profiling (Figs.1–3). We commonly observed among the data, as demonstrated for CG in Figure1, that the phenotype of the F1 hybrid (i.e. M0023-syn0 × M0081-syn0), which involved the low-phenotype parent (i.e. M0023-syn0) as a female, was typically close to the mid-parent value. This pattern of maternal influence attributable to parent-of-origin effect was previously observed in other studies involving maize kernel growth and its components (Groszmann and Sprague, 1948; Bagnara and Daynard, 1983). Phenotypic differences between the reciprocal F1 hybrids could also be clearly seen when the seeds were germinated at 12.5°C for 11 d in the dark (Fig. 2).
CG response of self-pollinated RILs and their reciprocal F1 hybrids. The differential response of the reciprocal F1 hybrids indicates maternal dominance (parent-of-origin effect).
Seeds of reciprocal F1hybrids germinated for 11 d at 12.5°C in the dark. Top panel, F1 (M0081-syn0 × M0023-syn0); bottom panel, F1 (M0023-syn0 × M0081-syn0). Parents, M0081-syn0, and M0023-syn0 showed a similar response to that of F1-M0081-syn0 × M0023-syn0 and the F1-M0023-syn0 × M0081-syn0, respectively (data not shown).
When the seeds were harvested at high moisture (40%) and dried at high temperature (43°C), the phenotypic differences between the reciprocal F1s in the DT experiment was similar to those in the CG experiment, with one exception. The phenotype of the F1 M0337-syn4 × M0021-syn4 having the low-DT parent (M0337-syn4) as female did not approach the mid-parent value (Fig. 3). Although the reciprocal F1hybrids were significantly different in their phenotype after drying at high (43°C) temperature, the phenotypic differences were not significant after drying at the low (35°C) temperature (Fig. 3). When the reciprocal F1 hybrids of CG and DT traits were self-pollinated in the subsequent year, the magnitude of phenotypic differences between them was noticeably reduced (Fig. 4). This suggests that the nuclear genome played a major role in determining phenotypic differences observed in the F1 generation, rather than cytoplasmic (mitochondrial and chloroplast) genes.
CG response of self-pollinated RILs and their reciprocal F1s after high-moisture harvest and high- or low-temperature drying. The differential response of the reciprocal F1 hybrids at high-temperature (43°C) drying but not at low-temperature (35°C) drying indicates maternal dominance (parent-of-origin effect).
CG response of self-pollinated reciprocal F1 hybrids. The phenotypic differences between the self-pollinated reciprocal F1 hybrids were significantly reduced for both CG and DT traits. This suggests that the phenotypic differences between the reciprocal F1hybrids were mainly due to nuclear genes rather than cytoplasmic (mitochondrial or chloroplast genomes) factors.
Transcript Profiling
cDNAs synthesized from poly(A+) RNA of germinating embryo axes (CG) and embryos at 40% moisture content (DT) were digested with unique pairs of restriction enzymes in separate reactions. cDNA restriction fragments resolved by capillary gel electrophoresis varied in length from 40 to 450 nt. Use of 89 and 88 restriction enzyme pairs yielded a total of 32,496 and 32,940 detectable cDNA fragments in CG and DT samples, respectively. Quantitative expression analysis of the individual fragments, as shown by the peak heights in the electrophoretic scans, revealed that only 336 (1%) and 656 (2%) of the cDNA fragments were differentially expressed with ≥1.5-fold change between the reciprocal F1 hybrids in CG and DT traits, respectively.
GeneCalling analysis identified the genes that matched the cDNA restriction fragments against the predicted fragments in the public and proprietary databases. Most of the fragments matched more than one gene in the databases. Also, many genes were “called” by more than one cDNA restriction fragment. A subset of the cDNA fragments and their corresponding gene matches were selected for confirmation through cloning, sequencing, and competitive PCR analyses. The selection was based on whether the fragment showed a high degree of fold change between the reciprocal F1 hybrids and/or whether multiple cDNA fragments, with a consistent pattern of differential expression (measured by t test), were indicated for a given gene. In cases where the same gene matched multiple cDNA fragments with similar levels of expression, only one of the fragments was included for analysis. The confirmed genes with their function, fold change in expression, and cDNA restriction fragment size are listed in TableI for CG and in TableII for DT. The gene names and accession numbers were obtained using BLAST search of the public databases (GenBank and SwissProt/TrEMBL) with probability value of ≤1 × 10−07. Competitive PCR analysis using sequence-specific primers was used for sequence confirmation (by GeneCalling) or cloned sequence identification. Confirmation was revealed by the selective disappearance of the corresponding peak in the electrophoretic scans (Fig. 5).
Differentially expressed genes for the CG trait revealed by mRNA profiling
Differentially expressed genes for the DT trait revealed by mRNA profiling
Competitive PCR confirmation of the differentially expressed mRNA between the reciprocal F1 hybrids. The traces in red are controls and those in green are of competitive PCR using gene specific primers (see “Materials and Methods” for details). A and B show traces of cDNA fragments of the reciprocal F1 hybrids that were high (M0081-syn0 × M0023-syn0) and low (M0023-syn0 × M0081-syn0), respectively, for the CG phenotype. Competitive PCR with specific primers to abscisic acid- and stress-inducible protein (AF09573; see Table I) caused the red (control) peak (shown by the red arrow) to disappear, confirming the gene identity. Similarly, C and D show the traces of the desiccation tolerant (M0021-syn4 × M0337-syn4) and sensitive (M0337-syn4 × M0021-syn4) reciprocal F1 hybrids, respectively, showing the competitive PCR confirmation with alliinase (AL035356; see Table II) specific primers. The red traces in each of the panels represent an average of six electrophoretic traces of the digests from replicate samples. The differential expression (i.e. fold change) of the cDNA fragments can be gauged by comparing the peak heights of the red traces between A and B and also between C and D. The length of the cDNA fragment is shown along the x axis as number of nucleotides, and the height of the peaks (y axis) represent the abundance in arbitrary fluorescence units.
Protein Expression Profiles
The average protein yields were 14% and 19.5% (dry weight basis) for embryos at 40% moisture content and 11-d-old germinating embryo axes, respectively, suggesting that most of the proteins were extracted from the tissue. The principle steps involved in the protein expression profiling are illustrated in Figure6.
Schematic representation of the protein profiling process. Protein extracts from the samples were resolved on two-dimensional gels, first by their pI followed by their relative molecular mass. The gels were stained with a highly sensitive fluorescent dye and were scanned for image data processing. Each protein feature on the gel was cataloged by its position (pI andM r) and density (volume) in the gel, referred to as PEM. The PEMs of replicate samples of each F1 hybrid were compiled and the differentially expressed protein features between the reciprocal F1 hybrids were identified. These protein features were then robotically excised and subjected to trypsin digestion (site-specific cleavage at the C-terminal of Lys or Arg residue), and masses of the resulting oligopeptides were measured using matrix-assisted laser-desorption ionization time of flight mass spectrometer (MS). A portion of the digest was also subjected to MS/MS analysis using nano-electrospray ionization source (Z-spray) on a quadrupole time of flight MS. The uninterpreted fragmentation spectra in combination with the masses of the tryptic fragments were used to query the databases (GenPept, SwissProt, and Pioneer/DuPont and Oxford GlycoSciences' proprietary) using a computer program (SEQUEST program, Fennigan Corp., San Jose, CA) to obtain the protein identification (for details, see “Materials and Methods”).
A total of 7,419 and 6,585 unique features (proteins with different pI/M r combinations) were detected across samples within the CG and DT traits, respectively. Of these, 2,641 and 1,876 features were consistently detected above the background threshold level for CG and DT experiments, respectively. The F-test identified 117 features in CG and 205 features in DT samples that were significantly differentially expressed (P < 0.1) with ≥1.5-fold change between the reciprocal F1s. A subset of these features was further analyzed to obtain the protein identification.
The selected protein features were excised from the gels and subjected to in-gel trypsin digestion. The database search using tryptic-oligopeptide mass and tandem mass spectral (MS/MS) information with SEQUEST program (Fennigan Corp.) resulted in the identification (annotation) of many proteins. The representative PEMs with the annotated proteins that are differentially expressed between the reciprocal F1 hybrids are shown in Figure7. The reason why some of the protein features are difficult to see on the gels is that they are expressed at low levels and can only be detectable by the highly sensitive scanners.
PEMs of representative samples for the CG (A) and DT (B) experiments with annotations. The protein features indicated by red spots were differentially expressed between the reciprocal F1 hybrids. The identification of each of the proteins and their descriptions are provided in Tables III and IV for A and B, respectively.
Protein identifications obtained for 56 features in the CG experiment (Fig. 7A) and 54 features in the DT experiment (Fig. 7B) are listed in Tables III and IV, respectively. Identification of the remaining features in each trait could not be obtained because of the database limitations. The oligopetide masses and MS/MS spectral data of some of the features matched more than one gene. For example, globulin-2 precursor and nucleoside diphosphate kinase I with pI 6.5 and M r 12,010 (protein feature no. 4,413; TableIII) or globulin 1L and superoxide dismutase with pI 4.9 andM r 13,670 (protein feature no. 3,544; TableIV) had either the same or different oligopeptide sequence matches. On the other hand, two or more protein features had the same annotation. For example, several protein features in both CG (TableIII) and DT (Table IV) were annotated to be globulin-2 precursor proteins.
Differentially expressed genes for the CG trait revealed by protein expression profiling
Differentially expressed genes for the DT trait revealed by protein expression profiling
DISCUSSION
Screening and Phenotypic Analysis
Significant variation observed for CG and DT in this population of RILs enabled us to select the most divergent inbreds for evaluation of parent-of-origin effect in reciprocal F1 hybrids. Because the RIL population was derived from a single cross (B73 × M017), only two possible alleles can occur for any gene at a given locus. Furthermore, each RIL was self-pollinated for at least 11 generations, creating a population of inbreds each of which is homogeneous (one haplotype) and homozygous at most of the loci (Burr and Burr, 1991). Therefore, this RIL population was ideal for studying the genetics of many traits, especially complex ones such as CG and DT included in this study. Comparison of reciprocal F1 hybrids produced by crossing divergent RILs provided an ideal system (such as comparison of isogenic lines) to identify genes that influence the phenotypes considered.
Of the 25 RIL combinations (five high × five low) in each trait, the majority did not show significant differences between their reciprocal F1 hybrids (data not shown). However, reciprocal F1 hybrids between specific divergent RILs showed significant phenotypic differences, suggesting that parent-of-origin effects were influencing these traits in a genotype-specific manner. Differential phenotypic expression between reciprocal F1 hybrids, as was shown in other studies, is considered primary evidence of the imprinting phenomenon in plants (Kermicle, 1970; Lin, 1982; Ohad et al., 1996; Chadhury et al., 1997; Grossniklaus et al., 1998). Pinnell (1949) and Bdliya and Burris (1988) reported phenotypic differences between the reciprocal F1 hybrids of selected maize genotypes for CG and DT traits, respectively, but the underlying molecular bases were not explored.
Cytoplasmic and Dosage Effects
It is not uncommon to see phenotypic differences between reciprocal F1 hybrids for many morphological and agronomic traits in maize, especially those caused by cytoplasmic factors (mitochondrial and chloroplast genomes, commonly referred to as the maternal effect). To discern whether such factors were involved, reciprocal F1 hybrids were self-pollinated, and phenotypes of the resulting progeny were evaluated. Phenotypic differences between the self-pollinated progeny of the reciprocal F1 hybrids were significantly reduced compared with the reciprocal F1 hybrids themselves for both traits. This suggested that there was only a negligible effect of cytoplasmic factors on the phenotypes observed (Fig. 4). Unlike endosperm tissue, the tissues used in this study (embryos at R5 stage and germinating embryo axes) were diploid, and, therefore, dosage of any given gene should be equal from both female and male parent. The genotype of the endosperm can have a significant effect on early embryo development (Chang and Neuffer, 1994), but, to our knowledge, it has not been shown to have an effect on embryos at later stages of development.
Imprinting Effect
We hypothesize that the phenotypic differences between the reciprocal F1 hybrids observed in this study were most likely due to the imprinting phenomenon, and the differentially expressed genes identified were either directly or indirectly associated with such phenomenon. Imprinting was extensively studied in plants involving endosperm genes in maize (Kermicle and Alleman, 1990) and Arabidopsis (Luo et al., 2000; Russinova and de Vries, 2000). One of the best characterized imprinted loci of plants is the Rlocus that controls the anthocyanin pigmentation in maize aleurone. Using maize B-A translocation stocks with segmental duplications and deletions of chromosome 10L, Kermicle (1970, 1978) demonstrated that some of the R alleles condition mottled pigmentation in the aleurone after paternal transmission. Similarly, Lin (1982), also using maize B-A translocations, provided evidence for the imprinting effect of endosperm factors (Ef genes) in normal endosperm growth. The Ef genes, some of which are also located on chromosome 10L, were shown to support normal growth of endosperm only when inherited from the male parent.
In this study, we observed various α-tubulin genes to be differentially expressed in embryo and germinating-embryo axis tissues between the reciprocal crosses in both CG (Table I) and DT (Table IV) traits. Lund et al. (1995b) reported differential demethylation of α-tubulin genes between embryo and endosperm tissues in maize where DNA demethylation correlated with the increased accumulation of tubα2 and tubα4 RNA. Furthermore, they also reported that the differential demethylation occurred in the endosperm tissue between reciprocal crosses of maize inbred lines (W64A and A69Y) and that the demethylation state of tubα3 and tubα4 was maternally transmitted.
The fertilization independent seed (FIS) genes (MEA or FIS1, FIS2, and FIEor FIS3) that repress Arabidopsis endosperm development in the absence of pollination are perhaps the most studied imprinted genes in plants at the molecular level (Grossniklaus et al., 1998; Luo et al., 2000). Earlier observations of imprinting (parent-specific transmission) of these genes have come from studying the reciprocal F1 hybrids involving wild-type and mutant alleles of FIS genes (Ohad et al., 1996). Vielle-Calzada et al. (2000) recently reported that only maternal alleles were expressed during early embryogenesis (up to 64-cell stage embryo) and endosperm development in Arabidopsis. However, the authors also reported that the paternal alleles were expressed 3 to 4 d after fertilization. In this study, we demonstrate that embryos derived from seed of reciprocally cross-pollinated inbreds can display divergent phenotypes and differential patterns of gene expression, despite the fact that they are genetically identical.
mRNA Profiling
We used GeneCalling technology of CuraGen Corporation (New Haven, CT) to profile differences of mRNA expression between the reciprocal F1 hybrids that displayed significant phenotypic divergence in their ability to tolerate high-temperature desiccation or germinate under cold conditions. Shimkets et al. (1999) showed that the GeneCalling profiling procedure could yield at least one detectable cDNA fragment from more than 90% of expressed genes with an average of 2.3 fragments per gene when 72 restriction enzyme pairs were used. In our study, we used 89 and 88 enzyme pairs for CG and DT traits, respectively, yielding >32,000 cDNA fragments in each case. Only a subset (< 2%) of these showed ≥1.5-fold change in their expression between the reciprocal F1 hybrids. Comparison of reciprocal F1s, which are genetically identical (also, we sampled diploid tissue), helped eliminate the effect of DNA polymorphism that can significantly increase false positives (the number of cDNA fragments not relevant in influencing the trait) in the GeneCalling process.
A subset of the differentially expressed cDNA fragments was selected for further analyses via competitive PCR or cloning and sequencing. Twenty-four genes were positively confirmed from a total of 336 cDNA fragments that were differentially expressed (≥1.5-fold change) between the high- and low-reciprocal F1s for CG (Table I). Similarly, the identification of 64 genes was positively confirmed out of 656 differentially expressed cDNA fragments between desiccation-tolerant and -sensitive reciprocal F1hybrids (Table II). The maximum expression fold change observed for the CG trait was 14.8 (Table I) and that in DT was 8.0 (Table II). BLAST analyses revealed that all the genes were in public databases, either in GenBank composite/expressed sequence tag (EST) or in SwissProt/TrEMBL, and belong to wide functional categories (Tables Iand II).
Protein Profiling
Systematic survey of subcellular or whole-cell proteins by either N-terminal Edman sequencing or tandem mass spectrometry of individual spots resolved on a two-dimensional gel followed by database searching is fast becoming a common method of choice in developing maps of proteomes in various organisms, including plants (Peltier et al., 2000). To our knowledge, a survey of differentially expressed proteins that are attributable to a specific trait (phenotype) has never been reported in plants. In this study, we compared protein profiles of reciprocal F1 hybrids that showed differential phenotypes for CG and DT traits.
Protein expression map (PEM) analysis with Oxford GlycoSciences' (Oxfordshire, UK) technology involved high-resolution quantitative detection (femtomol levels) of polypeptides followed by identification of proteins by database query using MS/MS spectral data (Page et al., 1999). In addition, PEM analysis can also provide insight into posttranslational modifications, proteolytic derivation, or isoforms of the same protein. For example, enolase 1 was annotated for two different protein features, numbers 2,518 and 2,519 (Table III; Fig.7A), with a mass difference of 680 d and a pI difference of 0.1. These values suggest that the former protein is a proteolitic derivative of the latter or that they are isoforms of same enzyme, perhaps localized in different subcellular compartments. In maize, enolase 1 (eno1) is induced under anaerobic stress buteno2 is constitutive. Lal et al. (1998) speculated that the regulation of these genes involves post-translational modification. MS/MS of a few of the protein features in this study were matched to two different proteins each. Such annotations include protein feature numbers 4,413, 2,869, and 1,888 in Table III and 3,544 and 2,304 in Table IV. Some of these have common sequence motifs suggesting that they have functional similarity or that they belong to the same regulatory pathway.
Differentially Expressed Genes
The differentially expressed genes revealed by mRNA profiling and by protein profiling were different in both CG and DT traits. Nonetheless, several genes with similar functions and biochemical pathways were represented in both types of expression profiling. Examples include glycolytic/gluconeogenic pathway enzymes, members of cell wall and polyamine metabolism, ABA/stress-induced proteins, components of the cytoskeleton, chaperonins, etc. (Tables I-IV). Although we cannot know with certainty, perhaps the same genes were not observed by both profiling methods within each trait because only a subset of the total population of differentially expressed cDNA fragments and proteins were identified. A more comprehensive analysis of expressed proteins would require cellular fractionation and enrichment methods to identify low-abundance proteins.
Some of the genes observed in this study were previously shown to be associated with biotic and abiotic stress response in maize during seed development and germination. These include ABA/stress-inducible proteins, heat shock proteins, late embryogenesis abundant (LEA) proteins, and proteases (Yacoob and Filion, 1986; Prasad, 1996;Nieto-Sotelo et al., 1999; White et al., 2000). Upon comparing the accession numbers of drought and cold-induced genes in Arabidopsis reported by Seki et al. (2001) with the differentially expressed genes in this study (Tables I-IV), we found none of the genes to be in common. A closer examination by the gene description revealed that the expression of several genes, including LEA proteins, enolase, Cys protease, glyoxalase I, nodulin, heat shock protein dnaJ, and β-glucosidase were changed in both studies. Expression of these genes was modulated in the same direction in both the studies, except one of the LEA proteins, enolase, and glyoxaylase I (Seki et al., 2001; TablesI-IV). This could be explained by differential embryo morphology, temporal and developmental differences, mono- and dicotyledonous nature, differences in experimental conditions, and differential responses of various isoforms of these genes to stress.
Protein profiling of reciprocal F1s for the CG trait showed modulation in expression of genes involved in protein degradation such as proteases and proteins associated with trafficking of ubiquitinated proteins to proteosomes (Fu et al., 1999; Tables I andIII). These proteins may play a role in mobilization of proteins from the scutellum to the growing meristems of the primary root and shoot during germination. Proteosome (26S and 20S) component proteins and proteases were also differentially expressed for the DT trait (TableII) and, in an analogous way, may be important in preconditioning the embryo so as to render it more tolerant of high-temperature stress during the final phase of maturation.
Protein expression profiling of reciprocal F1hybrids divergent for DT showed 30 different globulin polypeptides. That is more than double the number of previously reported globulin fractions (Cross and Adams, 1983a, 1983b). Globulins are encoded byGlb1 and Glb2 genes, which are present in only one or two copies each in the maize genome (Belanger and Kriz, 1991;Wallace and Kriz, 1991). However, the details of post-translational processing and assembly of such a large number of polypeptides into protein bodies remain unclear. It was interesting to note that the majority of globulin 1 fractions were down-regulated, whereas globulin 2 fractions were up-regulated in the desiccation-tolerant F1 hybrid. Sequence alignment of Glb1and Glb2 revealed no significant similarity. Further examination is required to understand the effect of globulins on DT in maize.
Some of the differentially expressed genes observed in this study were implicated in more than one function, including structural, enzymatic, and regulatory with various subcellular localizations. True functions of many of the genes are currently being studied. It is important to understand that both of these traits are complex in that the allelic condition of many genes at different loci, rather than a single or a few genes, may determine the phenotype. It is not surprising that these traits were associated with a large number of genes encompassing a broad range of biochemical pathways. This is, in fact, borne out by the mere number of genes we measured as differentially expressed in this tightly controlled comparison. It is reasonable to assume that pleiotropic effects of a subset of the differentially expressed genes caused altered expression of some of the other genes identified in this study. It is expected that among the genes identified as “differential” there lie those that are intimately involved in controlling cellular processes that lead to the manifestation of differences in CG and DT phenotypes.
MATERIALS AND METHODS
Germplasm Screening and Reciprocal F1 Hybrid Production
Recombinant inbred populations were generated from a cross between B73 and MO17. The inbred lines were derived by self-pollinating the original F1 plants (syn0) or sib-mating the original F1 plants for two (syn2) or four (syn4) cycles, followed by self-pollinating each line of the F2 progeny for at least 11 generations. These lines were screened for their ability to tolerate high desiccation temperature when harvested at high seed-moisture content (for DT) and to germinate under cold conditions (for CG and DT). The lab-based CG test was performed according to Byrum and Copeland (1995) with minor modifications. Seeds were germinated in rolls of wet germination paper and incubated at 10°C for 7 d followed by 3 d at 27°C in the dark. Seedlings that had ≥2.5 cm shoot and ≥2.5 cm primary root lengths were considered germinated. In some cases, the primary root and shoot lengths were also recorded as a measure of seedling vigor (data not shown).
For the DT trait, 100 RILs that showed high germination under cold test conditions when harvested at low seed moisture (approximately 25%, R6 stage or physiological maturity) and dried at low (35°C/normal) temperature were screened. The screening was done by harvesting self-pollinated ears at low (≤25%) or high (≥40%, R5, or dent stage) seed-moisture content and drying at low (35°C, relative humidity [RH] 51%) or high (43°C, RH 28%) temperatures followed by performing the CG test. For the CG trait, self-pollinated ears harvested after physiological maturity from 780 RILs were screened using the CG test. Five contrasting RILs for each trait, DT (tolerant: M0002-syn0, M0021-syn4, M0212-syn4, M0241-syn0, M0363-syn4; and sensitive: M0079-syn0, M0213-syn0, M0247-syn4, M0337-syn4, M0367-syn2) and CG (high: M0029-syn2, M0062-syn0, M0081-syn0, M00324-syn0, M00367-syn2; and low: M0023-syn0, M0028-syn0, M00106-syn4, M0113-syn2, M0342-syn0), were selected for crossing. The divergent RILs were self-pollinated and cross-pollinated reciprocally in a diallel experimental design. Seeds from the selected reciprocally crossed ears (F1s) were grown and self-pollinated in the following season to test for cytoplasmic inheritance of the traits.
Phenotypic Analysis and Sample Collection
Self-pollinated and reciprocally cross-pollinated ears for the DT study were harvested at high (R5 stage) and low (R6 stage) seed moisture. After drying treatments at low (35°C, RH 51%) and high (43°C, RH 28%) temperatures, CG was evaluated on three individual ears for each treatment as described earlier. Similarly, five mature individual self- and reciprocally cross-pollinated ears were tested for the CG trait.
Specific combinations of RILs were identified based on the phenotypic analyses where the reciprocal F1 hybrids displayed significant differences. In case of the DT trait, a few rows of kernels were removed from each ear before drying treatment to estimate moisture content and to dissect embryos (only from seeds harvested at 40% moisture content) for RNA and protein extraction. In case of the CG trait, embryo axes tissue of seed germinated for 11 d at 12.5°C in the dark were used to extract RNA and protein.
GeneCalling Analysis of mRNA from Reciprocal F1Hybrids
The GeneCalling procedure involves a comprehensive quantitative analysis of differentially expressed transcripts. Isolated mRNA is converted to double-stranded cDNA, which in turn is subjected to digestion with specific pairs of restriction enzymes (restriction fragments). Gene identification is accomplished by matching fragment size(s) using a computer algorithm to query against virtual digests of genes in a database (see Rothberg et al., 1999; Bruce et al., 2000). In brief, total RNA was isolated from 1 to 2 g of ground tissue using the PUREscript kit (Gentra Systems, Inc., Minneapolis). Isolation of poly(A+) RNA, cDNA preparations, restriction digestion, adapter ligation and PCR amplification, restriction fragment separation, comparison and identification of differentially expressed fragments, and gene identification and confirmation were performed according to Shimkets et al. (1999). The cDNA preparations of CG and DT samples were digested with 89 and 88 restriction enzyme pairs, respectively. Fragments from each digest were ligated to and PCR amplified with restriction enzyme-specific oligonucleotide adaptor/primer pairs. One of the primers was labeled with fluorescamine and the other with biotin. Labeled fragments were purified using streptavidin beads and were resolved by high-resolution capillary gel electrophoresis to generate traces showing peaks whose position and height represented M r and abundance of cDNA fragment(s), respectively. Trace data were used for qualitative (M r) and quantitative (abundance) comparisons between the samples from reciprocal F1hybrids.
The cDNA fragments that showed ≥1.5-fold change between the reciprocal F1 hybrids were identified from the trace database and were queried against publicly available maize (Zea mays) sequences and ESTs from the proprietary Pioneer/DuPont database for gene identification. The cDNA fragments that did not match any genes in the database were cloned from the original PCR reaction and sequenced using standard protocols. Genes identified through the GeneCalling process or via cloning of fragments were confirmed by a competitive PCR method in which the original PCR reaction was re-amplified in the presence or absence of an excess quantity of unlabeled gene/sequence-specific primers (a few additional bases into the template). Gene identification was positively confirmed if competitive PCR caused the chromatographic peak of the particular cDNA species to completely disappear. On occasion, in cases where there was more than one cDNA species at a single point of the chromatogram (i.e. different cDNA fragments with the same M r), then a partial reduction of peak height would occur, also confirming the gene identification.
All the cloned and sequenced cDNA fragments were subjected to BLASTX and BLASTN analyses by querying against public and proprietary databases (Altschul et al., 1990). Initial searches were performed against GenBank and SwissProt/TrEMBL, and genes with ≥93% identity of the entire query fragment were recorded. Remaining sequences were searched in the EST databases, and the accession numbers were recorded.
Protein Expression Profiling
Proteomic analysis involved high-resolution separation on two-dimensional gels, gel image analysis, robotic excision, MS/MS analysis, and protein identification by querying virtual spectra of database proteins (Fig. 6). Page et al. (1999) described in detail various steps involved in this procedure. The same processes were followed in this experiment with minor modifications as described below.
Proteins from frozen ground tissue were extracted in a buffer containing 8 m urea, 2 m thiourea, 65 mm dithiothreitol, and 4% (w/v) CHAPS. After 4 h of extraction on a shaker at 4°C, the supernatants were subjected to two-dimensional gel electrophoresis consisting of isoelectric focusing in the first dimension followed by SDS-PAGE in the second dimension. Three gels were run for each sample. The gels were stained with fluorescent dye OGT MP17 that binds non-covalently to the SDS moiety attached to the protein, enabling features at sub-nanogram quantities to be visualized. High-resolution (200 μm) scans of the gels were electronically stored in a database for image analyses. Each protein feature found across replicate gels was indexed (with pI and relative M r) along with its average (across replicate gels) intensity, and recorded as percentage of the total feature intensity (percent volume). Such definition of all the features in a given sample across the replicate gels is referred to as PEM. Samples from three replicate ears from each reciprocal F1hybrid of a given trait were organized into a table by the protein feature index numbers using Melanie II software (Wilkins et al., 1996).
Differential protein expression between the reciprocal F1hybrids was analyzed using statistical methods. The raw percent volume data of each protein feature across replicate gels and samples of each F1 genotype were subjected to analysis of variance. Features that showed significant percent volume difference (P ≤ 0.1) between the reciprocal F1s and were clearly visible in the gel images of at least one genotype were used to calculate ratios indicating fold expression differences. Fold change in expression was derived by taking the ratio of mean percent volume of a given feature between the reciprocal F1samples. Proteins that showed ≥1.5-fold change in expression were excised from the gel by a software-driven robotic cutter and delivered in a 96-well plate for proteolysis and MS analyses.
Excised proteins were subjected to trypsin digestion, and the mass of the oligopeptides were estimated by matrix-assisted laser-desorption ionization time of flight mass spectrometer (ELITE, PerSeptive Biosystems, Framingham, MA). For MS/MS analyses, nano-electrospray ionization source (Z-spray) on a quadrupole time of flight instrument (Micromass, Manchester, UK) was used to generate fragmentation spectra. These uninterpreted MS/MS of the peptides were converted to centered spectra and used to query public (GenPept and SwissProt) and proprietary (Pioneer/DuPont and Oxford GlycoSciences) databases using a computer algorithm (SEQUEST program, Fennigan Corp.) based on the approach of Eng et al. (1994). The sequences of the oligopeptides were derived by matching the predicted ion series from a complete peptide sequence in the database consistent with the observed y-type fragmentation of peptide pools in the sample (Biemann, 1990).
ACKNOWLEDGMENTS
We thank Emily van de Mortel (Pioneer Seed Physiology group) for her extensive technical support in RI population screening, crossing, and phenotypic analyses; Dale Wilson (Pioneer Seed Physiology group) for his help with statistical analysis; the Pioneer winter nursery staff in Puerto Vallarta, Mexico, for field support; and Oswald Crasta (CuraGen Corporation) and Jeremy Hayter (Oxford GlycoSciences Plc) for technical support and helpful discussions.
Footnotes
-
↵* Corresponding author; e-mailKrishna.Kollipara{at}pioneer.com; fax 515–334–4778.
-
Article, publication date, and citation information can be found at www.plantphysiol.org/cgi/doi/10.1104/pp.000729.
- Received November 20, 2001.
- Revision received January 21, 2002.
- Accepted March 5, 2002.