- Open Access
Gene-expression profile comparisons distinguish seven organs of maize
© Cho et al., licensee BioMed Central Ltd 2002
- Received: 6 May 2002
- Accepted: 5 July 2002
- Published: 29 August 2002
A maize array was fabricated with 5,376 unique expressed sequence tag (EST) clones sequenced from 4-day-old roots, immature ears and adult organ cDNA libraries. To elucidate organ relationships, relative mRNA levels were quantified by hybridization with embryos, three maize vegetative organs (leaf blades, leaf sheaths and roots) from multiple developmental stages, husk leaves and two types of floral organs (immature ears and silks).
Clustering analyses of the hybridization data suggest that maize utilizes both the PEPCK and NADP-ME C4 photosynthetic routes as genes in these pathways are co-regulated. Husk RNA has a gene-expression profile more similar to floral organs than to vegetative leaves. Only 7% of the genes were highly organ specific, showing over a fourfold difference in at least one of 12 comparisons and 37% showed a two- to fourfold difference. The majority of genes were expressed in diverse organs with little difference in transcript levels. Cross-hybridization among closely related genes within multigene families could obscure tissue specificity. As a first step in elucidating individual gene-expression patterns, we show that 45-nucleotide oligo probes produce signal intensities and signal ratios comparable to PCR probes on the same matrix.
Gene-expression profile studies with cDNA microarrays provide a new molecular tool for defining plant organs and their relationships and for discovering new biological processes in silico. cDNA microarrays are insufficient for differentiating recently duplicated genes. Gene-specific oligo probes printed along with cDNA probes can query individual gene-expression profiles and gene families simultaneously.
- Leaf Blade
- Oligo Probe
- Floral Organ
- Leaf Sheath
- Signal Ratio
Flowering plants are composed of diverse cell types organized into tissues and organs. To achieve the morphological and functional specialization of cells and tissues, suites of genes are expressed in spatial and temporal patterns determined by regulatory hierarchies responding to environmental cues. Plants reiteratively produce photosynthetic organs such as individual leaves that often differ in morphology and physiology depending on environmental cues and their positions on the plant. The fine-tuning of gene expression to permit such diversity remains largely uncharacterized. Flower development provides one example where organs of distinctive morphologies (sepals, petals, stamens, carpels) are produced in rapid succession; specification of each floral organ requires temporally and spatially refined expression of specific genes .
Although the functions and expression patterns for dozens of genes have been characterized in particular plant tissues and organs, the number of well studied examples is meager compared to the total number of genes .
Using information and material from genomic and high-throughput expressed sequence tag (EST) sequencing projects, several approaches have been devised to investigate global gene-expression profiles. In particular, spotted microarrays of ESTs have been used to initiate functional analyses of thousands of genes simultaneously. The first microarray contained only 45 Arabidopsis thaliana genes , but the demonstrated success of the method was quickly followed by studies of human  and yeast  gene expression. Microarray analysis has been used in plants for such diverse purposes as discovering genes responsible for strawberry flavor , comparing mutant to wild-type plants , and monitoring organism-level responses to environmental stimuli [8,9,10]. In most studies, treated and untreated tissues of the same age were compared. To date, there are just a few studies comparing distinct developmental stages. For example, Ruan et al.  surveyed expression in the three fundamental organ types - leaves, roots and flowers - of Arabidopsis, using microarrays containing 1,400 EST cDNA clones. Fernandes et al.  compared the hybridization of maize 14-day endosperm and immature ear (1-2 cm) RNA populations on separate endosperm and ear microarrays containing approximately 2,800 and 2,500 distinct genes, respectively. They found that nearly all probes on the ear array hybridized to cDNA prepared from endosperm or ear, whereas the endosperm array contained many apparently tissue-specific elements.
Using the 152,746 maize ESTs available, 31,858 tentative unique genes (TUGs) have been assembled as reported by the Zea mays database (April 14, 2002 ). Most ESTs are from the Maize Gene Discovery project , and a UniGene set is being developed from this resource. From the UniGene1 EST assembly 5,376 cDNA gene probes were used to fabricate a spotted cDNA microarray for a suite of hybridizations. In addition, 384 synthetic oligonucleotides 30, 40 or 45 nucleotides in length were printed to explore whether hybridization and washing conditions compatible with signal retention for both oligos and cDNA clones could be devised. Thirteen RNA samples from 7 distinct organs were hybridized on this array to ask how many genes were expressed in all or most organs, which genes had discrete patterns of expression, and whether gene-expression profiles could be used to understand the relationships among organs. Among > 5,000 genes selected for analysis, 56% showed less than a twofold difference and 37% showed a two- to fourfold difference in mRNA amounts compared to the reference sample - maize seedling RNA. These results imply that the differentiated state of maize tissues and organs is characterized by combinations of small numbers of differentially expressed (> 4-fold) tissue- or organ-specific genes among the 5,376 genes in this study. A complication of this interpretation is that some members of gene families will cross-hybridize, making it difficult to resolve whether there are large numbers of organ-specific genes within such families. As a first step in resolving expression patterns among recently duplicated genes, oligo probes for well characterized genes were printed on the same microarray slides as the cDNAs. The oligo probes generated more accurate information about gene expression than did cDNA probes for the genes examined in this study. We discuss how rationally designed oligo probes can distinguish expression patterns of individual genes in a gene family with similar sequences.
Internal consistency of hybridization
To examine the consistency of experiments from the labeling reactions through to the scanning process, a pool of mRNA from 4-day-old roots (4DR) was used to synthesize two fluorescent cDNA targets using Cy3-dUTP or Cy5-dUTP. These labeled cDNAs were subsequently combined in equal proportion to perform control hybridization experiments. Signal intensities of the two fluorescence-measurement channels were linearly correlated for a majority of the 5,376 PCR probes, indicating that nearly all targets were labeled with each dye. The mean Cy5/Cy3 ratio of the hybridized group of 5,263 genes was 1.046 with a standard error of 0.0014 and standard deviation of ± 0.105. Nearly all probes (5,074 of 5,263, or 96%) were within the 0.75 to 1.25 range (-0.415 to 0.32 on the log2 scale), and only two probe ratios were slightly over 2 (1 on a log2 scale). Of 5,376 EST probes examined, 5,263 generated signal intensities exceeding 300 intensity units in each channel and > 1,500 in the summation of the two channels, indicating that the ESTs on the array correspond mainly to moderately expressed genes. The few ESTs with signal intensities below these values have been omitted from further analysis. Another data set selected with a lower standard (> 700 in the summation of two channels) generated very similar results (data not shown).
Expression-profile comparisons between 4-day-old roots and immature ears
Signal-ratio distribution according to the source libraries
-2.5 ~ -2.0
-2.0 ~ -1.5
-1.5 ~ -1.0
-1.0 ~ -0.5
-0.5 ~ 0.0
0.0 ~ 0.5
0.5 ~ 1.0
1.0 ~ 1.5
1.5 ~ 2.0
2.0 ~ 2.5
Hybridization results were further analyzed by considering the cDNA library of origin. There were 1,920 probes on the slides from cDNA library 614 (4DR), 768 probes from 606 (IME), and 2,683 probes from 707 and 945 (two samples of the same mixed adult tissue cDNA library). In the hybridization control experiment (Figure 1), the behavior of 4DR dye-labeled cDNA sample on the 4DR (library 614) section of the slide was examined. Some probes of library 614 showed strong red on the slide in Figure 1a and strong green on the slide in Figure 1b in the 614 section (Figure 1, arrowheads), which would be expected if these clones were expressed in a tissue-specific pattern. Other EST elements, however, presented the opposite color pattern (Figure 1, arrows). Of 1,920 probes originating from a 4DR cDNA library (614), 1,106 (57.6%) showed log2 signal ratios between -0.5 and 0.5. Only 28 (15 + 13, 1.5%) showed over fourfold higher (log2 ratio > 2) mRNA amount in 4-day-old roots compared to immature ears (Table 1). Furthermore, 713 (37%) of 1,920 ESTs originating from a 4DR library were expressed at lower levels in 4DR than in IME. Similarly, probes originating from the IME (606) and mixed adult tissue (707 + 945) libraries showed only a small fraction of organ-specific expression elements (Table 1). These data, in part, reflect EST and UniGene1 consolidation methods (see Methods and materials). When an EST is chosen for a UniGene set, this does not mean that it was found only in one cDNA library source; contigs assembled from maize ESTs typically have contributions from several cDNA libraries.
Expression-profile comparisons among thirteen samples from seven organs
Signal-ratio distribution for 12 samples
(-2) ~ (-1)
(-1) ~ 1
1 ~ 2
2-week leaf sheaths
8-day leaf sheaths
2-week leaf blades
8-day leaf blades
Inferred hierarchy of organ similarity
The relationship of the husk, a photosynthetic 'leaf-like' organ  surrounding the ears of maize, to other organs merits special attention. It appeared as an immediate sister to the silk (Figure 3, bottom panel) or sister to the inflorescence organ group including the silk (Figure 3, upper panel). After collapsing ambiguous nodes, the husk remained within the reproductive group, and therefore distinctive from both leaf blade and sheath. Adjacent positions with the long terminal branches implied a loose relationship between husks and silks. Of the flower-related genes, 32 are expressed at high levels in husks (Figure 2, node c).
Organ-specific gene expression
To gain greater confidence in deciphering organ-specific gene-expression patterns, gene-product classification focused on 326 probes with a more than fourfold ratio difference in signal intensities in at least one organ compared to the reference (Figure 2). According to the criterion of E-value < e-10 in BLASTX searches , only 136 of the 326 genes predict similar proteins in public databases, 71 of 326 genes matched an Arabidopsis genome sequence, and 119 genes did not have significantly similar sequence in the public databases.
Flower-related genes included hydroxyproline-rich glycoprotein (hrgp ), β-amylase (PID9294660), pollen allergen (PID4006978), a bZIP transcription factor (PID6288682), four MADS-box genes [1,14], and several unknown genes (Figure 2). A relatively high expression of diverse transcription factors in embryo, immature ear and silk is consistent with microscopic observations that many stages of organ differentiation are occurring within the immature inflorescence and embryos. Because the husks were morphologically fully expanded leaf sheaths surrounding the ear, it was surprising that they expressed the same genes, such as MADS-box genes that are associated with early stages of flower development. Heat-shock proteins (16.9 kDa, 82 kDa, and 101 kDa) were abundant in silks, husks and IME, but not in embryos. The genes for the 82 kDa and 101 kDa proteins are expressed at raised temperatures [19,20]; however, in this study they appear to be part of a developmental program. Embryos and flowers are distinguished by large quantitative differences in expression of these three heat-shock genes. Root-specific genes included a nodulin homolog (PID3482914), putative lipid-transfer protein gene (PID10140658), physical impedance induced protein gene (PID2226329), and four additional unknown genes. The organ-specific expression pattern of these genes may spark interest in defining their physiological functions.
Genes expressed preferentially in the leaf blade
Twenty-six transcripts were comparatively abundant only in leaf blades (Figure 2, node a). The expression ratio of these genes was more than twofold higher in 8-day, 2-week and adult leaf blades, and more than fourfold lower in roots, immature ears and embryos, in comparison to the reference 4-day shoots. Twenty-three of these leaf-enriched transcripts shared high sequence similarity to previously published (identified or putative) coding sequences. Gene products from 17 of these 23 genes were previously characterized as located in or predicted to locate to plastids. Two well characterized genes in this leaf-blade group are for Rubisco small subunit (rbcS) and phosphoribulokinase, key enzymes for converting CO2 into carbohydrate via the Calvin cycle. Other genes represented components of light-harvesting complexes (photosystems I and II), chloroplastic aldolase , and the phosphoenolpyruvate translocator gene . Most of these highly expressed leaf genes are encoded in the nuclear genome, and the proteins are imported into chloroplasts. Interestingly, three of the 26 'leaf-blade' genes are known to be encoded in the chloroplast genome in maize, as in other flowering plants . None of them contained a poly(A)+ tail track in the EST sequences. Transcripts for these genes are so abundant in leaf blades that poly (A)+ selection apparently failed to remove them during mRNA purification for cDNA library construction and hybridization target labeling. These chloroplast-encoded genes consistently showed a co-regulated expression pattern, clustering with photosynthesis genes encoded in the nuclear genome. Thus, although they are contaminants, their expression patterns confirm the co-regulation of plastid and nuclear-encoded genes required to construct photosynthetically competent organelles.
Verification of hybridization ratios and ratio interpretation
The blot hybridization also appears to report ratios over a wider range. The ratios from the microarrays for EMB and silks were about 7 for both organs, but they were 23 and 79 for EMB and silks in the blot hybridization, respectively. These observations may reflect two features of microarray analysis: the nonlinearity of fluorescence excitation and the saturation of signal intensities for abundant transcripts. The signal ratios from microarrays in this study probably underestimate the actual difference in amount for abundant transcripts. Therefore a small absolute ratio should be interpreted as a reliable indicator for the presence of the transcript type in both samples. In this experiment log2 ratios from -0.5 to 0.5 are interpreted as simply indicating transcript presence without ascribing a difference in absolute amount. Similarly, we conclude that a fourfold difference detected by microarray hybridization indicates more than a fourfold difference in RNA abundance and could indicate organ-specificity of expression.
Hybridization pattern comparisons between oligos and cDNA
There are five gene families on the 326-element cladogram shown in the middle of Figure 2: three αtubulin genes, three βtubulin genes, three carbonic anhydrase genes, four MADS box genes, and five putative cellulose synthase genes. Individual genes in each gene family displayed almost identical expression patterns in all 13 samples. A few individual α-tubulin , β-tubulin , MADS box [15,29], and cellulose synthase  genes have been reported to be differentially expressed in maize or other plants. For these gene families, individual gene expression profiles on microarrays are obscured by cross-hybridization among family members [12,31].
To test whether oligonucleotide probes can be utilized together with cDNA probes to resolve individual gene contributions, multiple oligos were printed on the same glass slide microarrays with the EST probes. We wished to determine whether multiple oligos designed to the same gene would exhibit a coherent hybridization pattern and whether the oligos from a particular gene would cluster with known examples of genes co-regulated in vivo, a powerful test of the microarray . For this analysis, 582 probes were selected, a combination of oligo and EST probes. Oligo probes from five genes (α-tub, hrgp, rbcS, eEF1-α, pepc) met the selection criteria (see legend to Figure 7) of demonstrating a high ratio in at least one hybridization. Multiple oligo probes for each of these genes were printed, as illustrated in the gene models in Figure 7. After cluster analysis, the multiple oligos for each gene established well-separated groups in only one restricted branch of the cladogram of the 582-element data set (Figure 7, left panel). Within the rbcS block as an example, both 45-nucleotide probes from each of two exons appear as close neighbors; no other probes separate them. Such tight groups are characteristic of all 45-nucleotide oligos present in this cladogram. Similar results were produced from data that were neither median-centered nor normalized.
The multiple cDNAs for rbcS, hrgp, and α-tub genes cluster with the corresponding oligo probes. It is notable, however, that two of three hrgp cDNA probes showed fairly reduced ratios in several organs (Figure 6, triangles and squares). The signal ratios from these two PCR probes differ significantly from the other probe and from the oligo probes (p < 0.01 in a paired t-test). Differences between a cDNA and multiple oligo probes are particularly evident for pepc. Six of nine oligos are shown on the cladogram with functionally related genes (Figures 2, 7). Three others were excluded during the selection process because of weak hybridization. On the other hand, the cDNA probe was not selected, because it exhibited an insufficient absolute ratio difference. We suspect that cross-hybridization among pepc family members (or other genes) obscured the authentic gene expression from the cDNA probe. The ratio patterns of oligo probes in 12 pairs of duplicated hybridization experiments suggest that 45-mer oligos are a good alternative to gene-specific RNA blot hybridization to measure expression patterns of specific genes. Oligos of 30 and 40 nucleotides were also used successfully, although signal strength was weaker (red dots on the gene list in Figure 7). We conclude that the oligo hybridization patterns reflect transcript levels relatively accurately for the five genes presented in Figure 7. In fact, the representation is likely to be more accurate than that based on PCR products based on the RNA blot hybridization comparisons (Figure 6).
Gene-expression profiles among thirteen samples from seven maize organs were analyzed using cDNA microarrays containing 5,376 unique genes. In addition, oligonucleotide elements included within the same microarrays yielded consistent hybridization patterns; oligo probes are a promising tool for resolving gene - or even allele-specific expression patterns. The majority of genes showed similar hybridization ratios among diverse maize organs, and only 326 (~ 7%) genes appeared highly organ-specific with > 4-fold ratio difference in comparison to the reference 4-day seedling sample. An organ hierarchy based on gene-expression profiles indicated a close relationship among silks, immature ears and embryos. These organs appeared distinct from vegetative organs such as leaf blades, leaf sheaths and roots. Surprisingly, husks were clustered in the floral organ group. In addition, analyses of coordinated expression patterns of photosynthetic genes strongly suggested the presence of two C4 pathways in maize leaf blades. As with other microarray experiments, the newly recognized patterns of gene expression are the springboard for additional genetic and molecular experiments.
Internal consistency of the microarray hybridization
Internal consistency of the array hybridization results was demonstrated by five pieces of evidence. First, a control hybridization with one type of mRNA for which aliquots were labeled with different dyes generated signal ratios within 0.75- to 1.25-fold for 96% of the genes. Second, a dye-swapping hybridization with two samples of mRNA hybridized on separate slides yielded very similar expression profiles (Figure 1). Third, multiple probes for each of several gene family members for five families deposited at random positions generated similar hybridization patterns (Figure 2), suggesting that local effects on hybridization were negligible. Fourth, functionally related genes clustered together, demonstrating a coherent pattern in 12 pairs of hybridization analysis (Figures 2, 7). Fifth, groups of oligos designed to match different positions within several genes generated similar signal ratios, and each oligo group clustered together (Figure 7). Collectively, these facts indicate that hybridization with these microarrays containing a mixture of cDNA and oligo probes was internally consistent.
Organ identity of husks
Each organ is expected to have a unique combination of expressed genes, allowing organ identification and assessment of similarity with other organs as shown by the cladograms in Figure 2. It is interesting that the highly expressed genes in husks parallel what is found in other floral organs. Anatomically, the husks around an ear are composed primarily of leaf sheath with a reduced ligule region subtending a highly reduced leaf blade in most maize inbred lines. Husks are usually classified as modified photosynthetic leaves, with the assumption that they are vegetative organs on a branch that terminates in an ear [16,33]. According the work of Langdale and colleagues , maize husks express mainly the C3 pathway of carbon fixation in contrast to leaf blades in which C4 fixation predominates. We found that husk gene-expression profiles are distinctive from both leaf blades and sheaths. For example, Rubisco subunit-binding protein (PID1345582) was expressed at a similar level in husks, but at > 2-fold lower levels in leaf blades compared to the reference 4-day-old shoot (Figure 2). On the other hand, all other photosynthetic genes expressed at high levels in leaf blades were at low levels in husks, relative to seedlings. Physiologically, photosynthetic rates in husks are consistently measured to be around 20-fold lower in leaf blades . Both the expression pattern of photosynthetic genes and the low rate of carbon fixation in husks suggest that these are distinctive organs. In contrast, those genes that are highly expressed in silks and immature ears were expressed at comparable levels in husks (Figure 2 node c). They included hrgp , β-amylase (PID9294660), pollen allergen (PID4006978), a bZIP transcription factor (PID6288682), four MADS-box [1,14], three heat-shock proteins, and a dozen uncharacterized genes. Consistent with the hybridization results, one MADS box gene, ZAP1, has been reported to be expressed in the sterile organs of maize florets and in husks . The Arabidopsis homolog AP1 is also expressed in non-reproductive organs such as sepal and petal primordia . The close relationship of husks to maize floral organs shown by gene-expression profiles suggests that husks could be considered as photosynthetic floral organs arising from an inflorescence meristem.
Two types of C4 photosynthesis pathways in maize
C4 plants have been classified into three subgroups on the basis of the distinctive enzymes that decarboxylate C4 acids in the bundle sheath cells. Maize is a classic NADP-ME type C4 plant . Interestingly, we found that the enzyme PEPCK is expressed in a pattern similar to NADP-ME and two additional universal C4 enzyme genes. PEPCK catalyzes the reversible decarboxylation of oxaloacetate (OAA) to PEP. This enzyme has several proposed functions, such as gluco-neogenesis in germinating seeds, carbon recovery during senescence, nitrogen assimilation during seed development and decarboxylation of OAA in PEPCK-type C4 photosynthesis [36,37]. The comparatively low expression level of the PEPCK gene in 4-day-old shoots weakens the hypothesis that its major function in maize is for gluconeogenesis in greening seedling parts. Similarly, high-level expression in seedling leaves and adult leaves cannot be for senescence-related carbon recovery. We concur with recent proposals that the major role of PEPCK in green tissues is decarboxylation of OAA during C4 photosynthesis in maize. Maize leaves have PEPCK activity equal to 45% of the activity levels of a 'pure' PEPCK-type C4 plant, Panicum maximum . Furthermore, the enzyme activity was localized in bundle-sheath cells where CO2 is released from OAA for refixation in the Calvin cycle . The cDNA was cloned from libraries enriched for maize bundle-sheath cells .
In addition to PEPCK, two additional genes for key enzymes of the PEPCK pathway were expressed coordinately: alanine aminotransferase and aspartate aminotransferase. Previously, the proteins were undetected by western analysis in purified maize bundle-sheath cells, using antibodies against Panicum maximum aspartate aminotransferase and Cucumis sativus alanine aminotransferase. Detection failure could reflect weak antibody cross-reactivity or enzyme degradation during bundle-sheath cell isolation . By EST sequencing and microarray hybridizations all three PEPCK C4 pathway-specific genes are expressed similarly to the NADP-ME pathway genes. We therefore propose that the PEPCK-type C4 pathway is active in addition to the NADP-ME type C4 pathway for CO2 fixation in maize leaf blades (Figure 5).
Extensive hybridization overlap by diverse organs
The UniGene microarray contained 5,376 elements. On this microarray, most genes hybridized to RNA from diverse organ samples, and very few genes hybridized to RNA from just one sample. According to the array hybridization results, over 60% of the genes produced similar signal ratios (|log2| < 0.5) between the reference and each of 13 samples examined. Thus most transcript types appear to be present at near equivalent levels in many organs of the plant. The interpretation of organ differences reported here reflects results based on a subset (17%) of the current tentative unique genes defined by maize EST collections; many of the EST elements queried are likely to be moderately expressed genes. However, microarray hybridization result is consistent with DNA-RNA reassociation kinetic studies using multiple organs of tobacco plants . About 40% of tobacco genes were expressed in all organs examined (leaf, petal, anther, ovary, root, stem) and 10-40% of the genes were tissue or organ-specific by the criterion of RNA complexity. The apparent low number of tissue- or organ-specific genes observed in maize is also consistent with other microarray studies indicating that only around 25% of Arabidopsis genes displayed significant (> 2-fold) difference in three organ comparisons: seed, root and leaf  or root, leaf and flower . Similarly, only 24% of tested genes were distinguishable at three stages of strawberry ripening . Studies of the same organ from different treatment regimes, such as dark-grown and light-grown seedlings of Arabidopsis  showed only a 16% difference. During a more complete study of the circadian cycle only 2% of genes examined showed differential expression with a circadian rhythm . From the data available, it appears that plant organs differentially express only a small subset of unique genes and that physiological perturbations result in induction or repression of an even smaller number.
Why do so many genes appear to be expressed in diverse plant organs including those of maize? Some housekeeping genes are constitutively expressed in similar amounts in all organs to insure the maintenance of basic cellular processes. Differential expressions might be cell-type dependent; such differences might not be detected in this study because most tissues were mixtures of multiple cell types. In a gene-expression study during Poplar wood development, > 40% genes were differentially expressed in different development zones within the vascular meristem . Some genes may be expressed at similar RNA levels but protein levels are controlled post-transcriptionally. On microarrays, the correspondence between individual gene expression and hybridization signal is not exact. Cross-hybridization among similar sequences is a major complication in microarrays fabricated with cDNAs. Substantial cross-hybridization has been reported among sequences showing 85% similarity over 30 nucleotides . Cross-hybridization between related genes will be a profound problem in the analysis of gene-expression patterns in plants. About 70% of genes are duplicated in A. thaliana through both chromosome and local duplications [2,43]. Maize has undergone an allo-tetraploid chromosome duplication event within the past 11.4 million years, preceded by other genome-wide duplications [44,45]. In the available studies, however, there are many examples of maize duplicated genes expressed in different organs . For example, two duplicated transcription factor genes (p1, p2) regulating pholaphene pigment synthesis are expressed fairly exclusively in two sets of organs. p1 and p2 arose following a local gene duplication and insertion of multiple retroelements between p1 and p2. Subsequently, p1 acquired a new regulatory sequence 5' of the gene, probably explaining its new expression pattern . There are many retroelements and DNA transposons flanking maize genes, and they may contribute to the rapid divergence of transcriptional regulation . Another example comes from duplicated chalcone synthase genes (C2, Whp). They share over 94% sequence similarity but are differentially expressed . It is also evident that some duplicated genes are expressed redundantly at the same time in the same organ. For example, five copies of cellulose synthase genes  and five copies of eEF1-α genes  are coexpressed in diverse organs.
Redundant expression among duplicated genes
Sequence comparisons among TUGs assembled within individual EST sequencing projects 614 (4DR) and 606 (IME) provide anecdotal information about the expression modes of duplicated genes. TUGs sharing > 90% sequence similarity over 100 nucleotides were identified by BLAST . Within library 606 (immature ear) 32% (963/3,032) of the TUGs are similar at this criterion, and 33% (1290/3,879) of the TUGs defined within library 614 (seedling root) appear to be duplicated genes. When TUGs assembled from libraries 606 and 614 are compared to each other, 20% appear to be duplicated. The sequence comparisons among TUGs probably underestimate the number of duplicated genes because sequence data are incomplete. Comparisons of full-length cDNA sequences of each gene would increase the fractions of gene families both within and between these two libraries. Because microarray hybridization conditions cannot resolve the precise expression patterns of gene > 90% similar, the true fraction of constitutively expressed genes cannot be calculated.
An important question is what fraction of closely related duplicated genes are expressed differentially during the maize life cycle. For the moderately expressed class of genes 'discovered' by EST sequencing of specific developmental stages, it is striking that so many gene families are expressed in all 13 samples examined here. Functional redundancy among individual genes within a gene family would produce no detectable phenotype until all functionally redundant genes are mutated (see examples in ). Yet, mutations in individual maize genes within a large gene family can produce a visible phenotype. This evidence indicates that functional specialization has occurred. By RNA blot hybridizations, it is often observed that the relative amount of transcripts varies among individual genes within a family, suggesting that promoter divergence produces quantitative differences [30,49]. In some cases, mutation that eliminates expression from one gene-family member may be compensated by higher expression of other members; even if there is no visible phenotype, a molecular phenotype is predicted. cDNA microarrays are not sensitive enough to detect minor changes of expression patterns or differential expression of recently duplicated genes with the current hybridization condition (65°C hybridization and 55°C washing). Oligo probes appear to be a good alternative for analysis of individual gene-expression profiles, either in conjunction with PCR products or by themselves. Oligo probes can be designed to represent individual genes by exploiting even small polymorphisms. Our results show that suitable hybridization and washing conditions can be used for the analysis of PCR and oligo probes on the same microarray slide.
The maize strain in this study has the genetic background K55 (75%), W23 (20%), Robertson's Mutator (5%). Seedlings were grown under 100 μE/m2/sec constant illumination conditions of cool-white fluorescent light in a 27°C growth room. For the 4-day seedling with coleoptile reference sample, the shoot was harvested; other seedling samples were taken at 8 days or 14 days after planting. Field-grown plants were the source of most organ samples. The same genotype was planted in mid-June at the Stanford University Plant Growth Facility. Three immature ears (3-5 cm) were harvested after the tips of the husks had emerged from a leaf sheath; silks were excluded from immature ears. Husks were collected from the same ear. The two outermost husk layers were excluded, and all other inside layers were collected. Mature but unpollinated silks were harvested from two ears; the ears had been shoot-bagged to prevent pollen contamination on the silks. Adult leaf blades were taken from fully expanded leaves.
cDNA probe preparation
cDNA clones were from three, non-normalized cDNA libraries: a mixture of adult tissue (projects 707 and 945, W23 inbred line with active Mutator transposons), 4-day-old roots (project 614, W23 inbred line), and immature ears (project 606, Oh43 inbred line). The 5,376 cDNA clones chosen represent approximately 17% of the tentative unique genes in the April 2002 EST assembly . They were designated as UniGene1 members after EST assembly of 73,000 available ESTs representing around 17,000 TUGs in September 2000. A clone for each TUG was selected on the basis of EST sequence length during UniGene1 consolidation. In many cases ESTs defining particular TUGs were recovered from multiple libraries. Clone identities were verified in UniGene1 by resequencing for approximately 50% of the clones to confirm well positions in the consolidation plates.
PCR amplifications of cDNA inserts were carried out at annealing temperatures of 50°C (614 and 606 libraries) or 60°C (707 library) in a 25 μl volume in a GeneAmp PCR system 9700 (Applied Biosystems, Foster City, CA) thermocycler for 35 cycles with a 2 min extension time. The reaction cocktail contained 1 ng EST plasmid DNA, 1.7 mM MgCl2, 1x reaction buffer (50 mM Tris pH 8.5, 20 mM KCl), and 0.1 units of Taq polymerase (GibcoBRL, Gaithersburg, MD). PCR-amplified products were purified with Gene Clean kits (BIO101, Carlsbad, CA), and eluted in 20 μl water. Samples of 3 μl were loaded onto a 1% agarose gel and electrophoresed to measure product size and yield. Of the 5,376 ESTs 197 produced multiple bands or smeared products and were excluded from the analyses. Samples of 10 μl were transferred to a 384-well plate and dried completely; the pellet was dissolved before printing in 5 μl of 150 mM phosphate pH 7.0 buffer to yield approximately 300 ng/μl DNA concentration. Probes were printed on 3D-link slides, followed by coupling and processing as recommended by the manufacturer (SurModics, Eden Prairie, MN).
Oligo probe preparation
A minimum of one oligo was synthesized within each exon, intron, and at exon/intron and exon1/exon2 junction regions from 17 selected maize genes. In most cases, two probes were synthesized in each exon and intron. Oligo design was based on double-stranded complete gene sequence, available from GenBank. Oligos were synthesized using phosphoramidite chemistry on an automated oligo-nucleotide synthesizer at the Stanford Genome Technology Center . Oligos were synthesized from the 3' to 5' direction, and the 5' end of each oligo was modified by addition of a C6-amide group. A total 184 oligos of 45 nucleotides were synthesized from the selected 17 genes. In addition, 96 oligos of 40 nucleotides and 96 oligos of 30 nucleotides were synthesized from 45 different genes for comparison of their hybridization behavior to PCR probes; short oligos were also prepared for the 17 genes for which the 45-nucleotides oligos were designed. The calculated Tm of exon probes ranged from 92-109°C, while the Tm of intron and intron/exon junction probes ranged from 89-95°C. Synthesized oligos were dissolved in 150 mM phosphate buffer, pH 8.5 at 40 μm. Multiple 45-nucleotide oligo probes for 17 genes were printed on the same arrays as the cDNA probes. Mean signal intensities were calculated from each of 10 genes represented by these 45-mers that showed consistent hybridization in all organ comparisons. They were used to examine the consistency of hybridization as positive control elements, and they were included in the cluster analyses.
RNA purification, labeling, and hybridization
Total RNA was extracted from 13 samples, using the Trizol method (GibcoBRL). mRNA was further purified from total RNA with Oligotex mini-columns (Qiagen, Valencia, CA). mRNA quantity and quality were examined by UV absorption at 260 and 280 nm. RNA quality was also examined by agarose gel electrophoresis to monitor loss of ribosomal RNA after mRNA purification. About 2 μg of poly (A)+ RNA was used to synthesize fluorescently labeled cDNA targets. The reaction cocktail contained ~2 μg poly(A)+ RNA, 1x reaction buffer (50 mM Tris-HCl pH 8, 75 mM KCl, 3 mM MgCL2, 50 μM dNTP, 10 mM DTT), 1 μg oligo dT, 3 μg random hexamer, and 400 units Superscript II (GibcoBRL).
Hybridization was performed as described at . Variations included hybridization at temperatures between 61-65°C and an initial washing at 55°C. mRNA from 4-day-old shoots with coleoptiles was labeled with Cy3-dUTP, which served as the common reference in all pairs of hybridization for the cluster analyses with Cy5-dUTP labeled samples from other stages. Microarray slides were scanned with an Axon400 scanner (Axon Instrument, Union City, CA). Signal was initially normalized during the image scanning process to adjust the average ratios between two channels. Grids were generated and adjusted automatically then refined manually to identify the microarray elements. Those probes whose signal intensity, subtracted by background, was lower than 300 in either channel or less than 1,500 in the sum of both channels were excluded from further analyses. Signal ratios for each probe element on each slide were calculated, using the mean intensity of pixels subtracted by median background for each channel. Array results are deposited at a public gene-expression database, Gene Expression Omnibus , and their accession numbers are GPL12 for the platform and GSM57-GSM80 for 24 samples.
Hierarchical clustering of the data was performed using the computer program Cluster . The output was visualized using the program TreeView (available at ). Cluster analyses were carried out before and after a secondary normalization process to make the sum of the squares 1.0 in each row and column. Although the results were very similar, we prefer the results from the unnormalized data for three reasons. First, the reference sample is identical in all hybridizations. Second, unnormalized data produced organ relationships consistent with organ identities and the relationships inferred from normalized data. Third, the normalization could compound variation by combining an uncertainty from a computation method on top of the variations from hybridization.
RNA blot hybridization
A 15 μg sample of total RNA was loaded onto a glyoxal gel as described in . Hybridization probes for RNA blots were prepared by the random primer labeling method to incorporate 32P. Blots were analyzed on a Phosphorlmager (Molecular Dynamics, Sunnyvale, CA).
We thank Brian Nakao, Gurpreet Randhawa and Khaled Sarsour for generating ESTs and UniGene1 verification sequencing, and ZmDB curators for database maintenance. We extend special thanks to Bret Schneider, Darren Morrow, Paula Casati, Mathew Fitzgerald and Dean Goodman for critical reading of the manuscript. The work is supported by National Science Foundation grant 98-72657 to V.W.
- Jack T: Plant development going MADS. Plant Mol Biol. 2001, 46: 515-520. 10.1023/A:1010689126632.PubMedView ArticleGoogle Scholar
- The Arabidopsis Genome Initiative: Analysis of the genome sequence of the flowering plant Arabidopsis thaliana. Nature. 2000, 408: 796-815. 10.1038/35048692.View ArticleGoogle Scholar
- Schena M, Shalon D, Davis RW, Brown PO: Quantitative monitoring of gene-expression patterns with a complementary-DNA microarray. Science. 1995, 270: 467-470.PubMedView ArticleGoogle Scholar
- Schena M, Shalon D, Heller R, Chai A, Brown PO, Davis RW: Parallel human genome analysis: microarray-based expression monitoring of 1000 genes. Proc Natl Acad Sci USA. 1996, 93: 10614-10619. 10.1073/pnas.93.20.10614.PubMedPubMed CentralView ArticleGoogle Scholar
- DeRisi JL, Iyer VR, Brown PO: Exploring the metabolic and genetic control of gene expression on a genomic scale. Science. 1997, 278: 680-686. 10.1126/science.278.5338.680.PubMedView ArticleGoogle Scholar
- Aharoni A, Keizer LC, Bouwmeester HJ, Sun Z, Alvarez-Huerta M, Verhoeven HA, Blaas J, van Houwelingen AM, De Vos RC, van der Voet H, et al: Identification of the SAAT gene involved in strawberry flavor biogenesis by use of DNA microarrays. Plant Cell. 2000, 12: 647-662. 10.1105/tpc.12.5.647.PubMedPubMed CentralView ArticleGoogle Scholar
- Helliwell CA, Chin-Atkins AN, Wilson IW, Chapple R, Dennis ES, Chaudhury A: The Arabidopsis amp1 gene encodes a putative glutamate carboxypeptidase. Plant Cell. 2001, 13: 2115-2125. 10.1105/tpc.13.9.2115.PubMedPubMed CentralView ArticleGoogle Scholar
- Reymond P, Weber H, Damond M, Farmer EE: Differential gene expression in response to mechanical wounding and insect feeding in Arabidopsis. Plant Cell. 2000, 12: 707-720. 10.1105/tpc.12.5.707.PubMedPubMed CentralView ArticleGoogle Scholar
- Seki M, Narusaka M, Abe H, Kasuga M, Yamaguchi-Shinozaki K, Carninci P, Hayashizaki Y, Shinozaki K: Monitoring the expression pattern of 1300 Arabidopsis genes under drought and cold stresses by using a full-length cDNA microarray. Plant Cell. 2001, 13: 61-72. 10.1105/tpc.13.1.61.PubMedPubMed CentralView ArticleGoogle Scholar
- Desprez T, Amselem J, Caboche M, Hofte H: Differential gene expression in Arabidopsis monitored using cDNA arrays. Plant J. 1998, 14: 643-652. 10.1046/j.1365-313X.1998.00160.x.PubMedView ArticleGoogle Scholar
- Ruan Y, Gilmore J, Conner T: Towards Arabidopsis genome analysis: monitoring expression profiles of 1400 genes using cDNA microarrays. Plant J. 1998, 15: 821-833. 10.1046/j.1365-313X.1998.00254.x.PubMedView ArticleGoogle Scholar
- Fernandes J, Brendel V, Gai X, Lal S, Chandler VL, Elumalai R, Galbraith DW, Pierson E, Walbot V: Comparison of RNA expression profiles based on maize EST frequency analysis and microarray hybridization. Plant Physiol. 2002, 128: 896-910. 10.1104/pp.010681.PubMedPubMed CentralView ArticleGoogle Scholar
- ZmDB: a maize genome database: class browser. [http://www.zmdb.iastate.edu]
- Heuer S, Hansen S, Bantin J, Brettschneider R, Kranz E, Lorz H, Dresselhaus T: The maize mads box gene zmmads3 affects node number and spikelet development and is co-expressed with zmmads1 during flower development, in egg cells, and early embryogenesis. Plant Physiol. 2001, 127: 33-45. 10.1104/pp.127.1.33.PubMedPubMed CentralView ArticleGoogle Scholar
- Mena M, Mandel MA, Lerner DR, Yanofsky MF, Schmidt RJ: A characterization of the MADS-box gene family in maize. Plant J. 1995, 8: 845-854.PubMedView ArticleGoogle Scholar
- Cheng P-C, Pareddy D: Morphology and development of the tassel and ear. In The Maize Handbook. Edited by: Freeling M, Walbot V. 1993, New York: Springer-Verlag, 37-47.Google Scholar
- Altschul SF, Madden TL, Schäffer AA, Zhang J, Zhang Z, Miller W, Lipman DJ: Gapped BLAST and PSI-BLAST: a new generation of protein database search programs. Nucleic Acids Res. 1997, 25: 3389-3402. 10.1093/nar/25.17.3389.PubMedPubMed CentralView ArticleGoogle Scholar
- Hood EE, Hood KR, Fritz SE: Hydroxyproline-rich glycoproteins in cell walls of pericarp from maize. Plant Sci. 1991, 79: 13-22. 10.1016/0168-9452(91)90063-E.View ArticleGoogle Scholar
- Nieto-Sotelo J, Kannan KB, Martinez LM, Segal C: Characterization of a maize heat-shock protein 101 gene, HSP101, encoding a ClpB/Hsp100 protein homologue. Gene. 1999, 230: 187-195. 10.1016/S0378-1119(99)00060-8.PubMedView ArticleGoogle Scholar
- Van Breusegem F, Dekeyser R, Garcia AB, Claes B, Gielen J, Van Montagu M, Caplan AB: Heat-inducible rice hsp82 and hsp70 are not always co-regulated. Planta. 1994, 193: 57-66.PubMedView ArticleGoogle Scholar
- Michelis R, Gepstein S: Identification and characterization of a heat-induced isoform of aldolase in oat chloroplast. Plant Mol Biol. 2000, 44: 487-498. 10.1023/A:1026528319769.PubMedView ArticleGoogle Scholar
- Fischer K, Arbinger B, Kammerer B, Busch C, Brink S, Wallmeier H, Sauer N, Eckerskorn C, Flugge UI: Cloning and in vivo expression of functional triose phosphate/phosphate translocators from C3- and C4-plants: evidence for the putative participation of specific amino acid residues in the recognition of phosphoenolpyruvate. Plant J. 1994, 5: 215-226. 10.1046/j.1365-313X.1994.05020215.x.PubMedView ArticleGoogle Scholar
- Maier RM, Neckermann K, Igloi GL, Kossel H: Complete sequence of the maize chloroplast genome: gene content, hotspots of divergence and fine tuning of genetic information by transcript editing. J Mol Biol. 1995, 251: 614-628. 10.1006/jmbi.1995.0460.PubMedView ArticleGoogle Scholar
- Ku MS, Kano-Murakami Y, Matsuoka M: Evolution and expression of C4 photosynthesis genes. Plant Physiol. 1996, 111: 949-957. 10.1104/pp.111.4.949.PubMedPubMed CentralView ArticleGoogle Scholar
- Galvez S, Hirsch AM, Wycoff KL, Hunt S, Layzell DB, Kondorosi A, Crespi M: Oxygen regulation of a nodule-located carbonic anhydrase in alfalfa. Plant Physiol. 2000, 124: 1059-1068. 10.1104/pp.124.3.1059.PubMedPubMed CentralView ArticleGoogle Scholar
- Furumoto T, Hata S, Izui K: cDNA cloning and characterization of maize phosphoenolpyruvate carboxykinase, a bundle sheath cell-specific enzyme. Plant Mol Biol. 1999, 41: 301-311. 10.1023/A:1006317120460.PubMedView ArticleGoogle Scholar
- Uribe X, Torres MA, Capellades M, Puigdomenech P, Rigau J: Maize alpha-tubulin genes are expressed according to specific patterns of cell differentiation. Plant Mol Biol. 1998, 37: 1069-1078. 10.1023/A:1006067710312.PubMedView ArticleGoogle Scholar
- Hussey PJ, Haas N, Hunsperger J, Larkin J, Snustad DP, Silflow CD: The beta-tubulin gene family in Zea mays: two differentially expressed beta-tubulin genes. Plant Mol Biol. 1990, 15: 957-972.PubMedView ArticleGoogle Scholar
- Cordts S, Bantin J, Wittich PE, Kranz E, Lorz H, Dresselhaus T: ZmES genes encode peptides with structural homology to defensins and are specifically expressed in the female gametophyte of maize. Plant J. 2001, 25: 103-114. 10.1046/j.0960-7412.2000.00944.x.PubMedView ArticleGoogle Scholar
- Holland N, Holland D, Helentjaris T, Dhugga KS, Xoconostle-Cazares B, Delmer DP: A comparative analysis of the plant cellulose synthase (CesA) gene family. Plant Physiol. 2000, 123: 1313-1324. 10.1104/pp.123.4.1313.PubMedPubMed CentralView ArticleGoogle Scholar
- Xu W, Bak S, Decker A, Paquette SM, Feyereisen R, Galbraith DW: Microarray-based analysis of gene expression in very large gene families: the cytochrome P450 gene superfamily of Arabidopsis thaliana. Gene. 2001, 272: 61-74. 10.1016/S0378-1119(01)00516-9.PubMedView ArticleGoogle Scholar
- Eisen MB, Spellman PT, Brown PO, Botstein D: Cluster analysis and display of genome-wide expression patterns. Proc Natl Acad Sci USA. 1998, 95: 14863-14868. 10.1073/pnas.95.25.14863.PubMedPubMed CentralView ArticleGoogle Scholar
- Langdale JA, Zelitch I, Miller E, Nelson T: Cell position and light influence C4 versus C3 patterns of photosynthetic gene expression in maize. EMBO J. 1998, 7: 3643-3651.Google Scholar
- Mandel MA, Gustafson-Brown C, Savidge B, Yanofsky MF: Molecular characterization of the Arabidopsis floral homeotic gene APETALA1. Nature. 1992, 360: 273-277. 10.1038/360273a0.PubMedView ArticleGoogle Scholar
- Martin R, Niyogi K: Photosynthesis. In Biochemistry and Molecular Biology of Plants. Edited by: Buchanan BB, Gruissem W, Jones RL. 2000, Rockville MD: American Society of Plant Physiologists, 619-628.Google Scholar
- Chen ZH, Walker RP, Acheson RM, Tecsi LI, Wingler A, Lea PJ, Leegood RC: Are isocitrate lyase and phosphoenolpyruvate carboxykinase involved in gluconeogenesis during senescence of barley leaves and cucumber cotyledons?. Plant Cell Physiol. 2000, 41: 960-967. 10.1093/pcp/pcd021.PubMedView ArticleGoogle Scholar
- Walker RP, Chen Z-H, Tecsi LI, Ramiani F, Lea PJ, Leegood RC: Phosphoenlpyruvate carboxykinase plays a role in interactions of carbon and nitrogen metabolism during grape seed development. Planta. 1999, 210: 9-18. 10.1007/s004250050648.PubMedView ArticleGoogle Scholar
- Wingler A, Walker RP, Chen ZH, Leegood RC: Phosphoenolpyruvate carboxykinase is involved in the decarboxylation of aspartate in the bundle sheath of maize. Plant Physiol. 1999, 120: 539-546. 10.1104/pp.120.2.539.PubMedPubMed CentralView ArticleGoogle Scholar
- Kamalay JC, Goldberg RB: Organ-specific nuclear RNAs in tobacco. Proc Natl Acad Sci USA. 1984, 81: 2801-2805.PubMedPubMed CentralView ArticleGoogle Scholar
- Girke T, Todd J, Ruuska S, White J, Benning C, Ohlrogge J: Microarray analysis of developing Arabidopsis seeds. Plant Physiol. 2000, 124: 1570-1581. 10.1104/pp.124.4.1570.PubMedPubMed CentralView ArticleGoogle Scholar
- Schaffer R, Landgraf J, Accerbi M, Simon V, Larson M, Wisman E: Microarray analysis of diurnal and circadian-regulated genes in Arabidopsis. Plant Cell. 2001, 13: 113-123. 10.1105/tpc.13.1.113.PubMedPubMed CentralView ArticleGoogle Scholar
- Hertzberg M, Aspeborg H, Schrader J, Andersson A, Erlandsson R, Blomqvist K, Bhalerao R, Uhlen M, Teeri TT, Lundeberg J, et al: A transcriptional roadmap to wood formation. Proc Natl Acad Sci USA. 2001, 98: 14732-14737. 10.1073/pnas.261293398.PubMedPubMed CentralView ArticleGoogle Scholar
- Vision TJ, Brown DG, Tanksley SD: The origins of genomic duplications in Arabidopsis. Science. 2000, 290: 2114-2117. 10.1126/science.290.5499.2114.PubMedView ArticleGoogle Scholar
- Gaut BS, Doebley JF: DNA sequence evidence for the segmental allotetraploid origin of maize. Proc Natl Acad Sci USA. 1997, 94: 6809-6814. 10.1073/pnas.94.13.6809.PubMedPubMed CentralView ArticleGoogle Scholar
- Wilson WA, Harrington SE, Woodman WL, Lee M, Sorrells ME, McCouch SR: Inferences on the genome structure of progenitor maize through comparative analysis of rice, maize and the domesticated panicoids. Genetics. 1999, 153: 453-473.PubMedPubMed CentralGoogle Scholar
- Franken P, Niesbach-Klosgen U, Weydemann U, Marechal-Drouard L, Saedler H, Wienand U: The duplicated chalcone synthase genes C2 and Whp (white pollen) of Zea mays are independently regulated: evidence for translational control of Whp expression by the anthocyanin intensifying gene in. EMBO J. 1991, 10: 2605-2612.PubMedPubMed CentralGoogle Scholar
- Zhang P, Chopra S, Peterson T: A segmental gene duplication generated differentially expressed myb-homologous genes in maize. Plant Cell. 2000, 12: 2311-2322. 10.1105/tpc.12.12.2311.PubMedPubMed CentralView ArticleGoogle Scholar
- Zhang Q, Arbuckle J, Wessler SR: Recent, extensive, and preferential insertion of members of the miniature inverted-repeat transposable element family Heartbreaker into genic regions of maize. Proc Natl Acad Sci USA. 2000, 97: 1160-1165. 10.1073/pnas.97.3.1160.PubMedPubMed CentralView ArticleGoogle Scholar
- Carneiro NP, Hughes PA, Larkins BA: The eEFIA gene family is differentially expressed in maize endosperm. Plant Mol Biol. 1999, 41: 801-813. 10.1023/A:1006391207980.PubMedView ArticleGoogle Scholar
- Lashkari DA, Hunicke-Smith SP, Norgren RM, Davis RW, Brennan T: An automated multiplex oligonucleotide synthesizer: development of high-throughput, low-cost DNA synthesis. Proc Natl Acad Sci USA. 1995, 92: 7912-7915.PubMedPubMed CentralView ArticleGoogle Scholar
- The Brown Lab: protocols. [http://cmgm.stanford.edu/pbrown/protocols/]
- Gene Expression Omnibus. [http://www.ncbi.nlm.nih.gov/geo/]
- Eisen lab. [http://rana.stanford.edu/software/]
- Sambrook J, Fritsch EF, Maniatis T: Molecular Cloning: A Laboratory Manual. 1989, Cold Spring Harbor, NY: Cold Spring Harbor Laboratory Press, 2Google Scholar