| HOME | HELP | FEEDBACK | SUBSCRIPTIONS | ARCHIVE | SEARCH | TABLE OF CONTENTS |

* Boyce Thompson Institute, Cornell University, Ithaca, New York 14853
Plant Biology Department, Cornell University, Ithaca, New York 14853
1 Corresponding author: Boyce Thompson Institute, Cornell University, Tower Rd., Ithaca, NY 14853.
E-mail-tpb8{at}cornell.edu
| ABSTRACT |
|---|
|
|
|---|
Over the past few years detailed genetic and biochemical studies of eudicots have greatly increased our understanding of how phytochrome signaling pathways regulate plant development (CASAL et al. 2003). In Arabidopsis, detailed mutant and molecular genetic analyses have defined roles for all five phytochrome apoproteins present in the genome. phyA is the primary photoreceptor controlling high-irradiance FR light responses such as hypocotyl elongation, cotyledon expansion, and seed germination (NAGATANI et al. 1993; PARKS and QUAIL 1993; WHITELAM et al. 1993). phyB is the primary red light photoreceptor mediating many of the shade-avoidance syndrome responses that result in petiole elongation, internode elongation, reduced leaf area, and early flowering (REED et al. 1993). phyC mediates seedling responses to red light and represses flowering under short-day conditions (FRANKLIN et al. 2003; MONTE et al. 2003). Although phyD and phyE are structurally similar to phyB, they are not functionally redundant (SHARROCK et al. 2003a). Responses mediated by phyD and phyE include petiole and internode elongation and the control of flowering time (AUKERMAN et al. 1997; DEVLIN et al. 1998, 1999).
Despite our detailed understanding of phytochrome signaling in eudicots, little is known of monocot light signal transduction. The monocots and, in particular, the grasses are by far the most agronomically important group of plants. Rice, maize, wheat, and sorghum are staple food crops throughout the world and have been under intense human selection for the last 500010,000 years. Furthermore, traits such as yield, flowering time, and stature are all influenced by red/far-red light (SMITH 1995). Early flowering varieties of barley (HANUMAPPA et al. 1999), sorghum (CHILDS et al. 1997), rice (IZAWA et al. 2000), and maize (SAWERS et al. 2002) have all recently been shown to contain lesions in the phytochrome signal transduction pathway. In maize, light response is a highly variable trait and a reduced light response is associated with early flowering North American and European inbred lines (MARKELZ et al. 2003). Together, these studies indicate that phytochrome response pathways contribute to the agronomic performance of these important crop species.
The phytochrome gene family structure in Arabidopsis contains five members (PHYAE); however, characterization of the gene family in other eudicots has revealed variation in gene family composition (HOWE et al. 1998; ALBA et al. 2000; LI and CHINNAPPA 2003). Unlike eudicots, the phytochrome apoprotein gene family in all monocots examined contains only three lineages: PhyA, PhyB/D, and PhyC (MATHEWS and SHARROCK 1996, 1997). In maize, a recent polyploidization event has resulted in large segmental duplications (GAUT and DOEBLEY 1997). Recent mapping experiments suggest that two PhyB sequences (PhyB1 and PhyB2) and two PhyA gene sequences (PhyA1 and PhyA2) arose from this genome duplication and are located in homeologous regions of the maize genome (CHRISTENSEN and QUAIL 1989; CHILDS et al. 1997; BASU et al. 2000). Preliminary studies also suggest that as many as three copies of PhyC-like genes may be present in the maize genome (BASU et al. 2000). However, only the PhyA2 gene of maize has been characterized in any detail (CHRISTENSEN and QUAIL 1989).
To examine evolutionary and functional relationships between phytochrome gene family members in maize, a detailed analysis of phytochrome gene family structure and expression has been performed. We report on the cloning and sequence analysis of full-length PhyA, PhyB, and PhyC genes in maize. Screens of bacterial artificial chromosome (BAC) clones generated from the standard maize inbred B73 indicate that the entire phytochrome gene family is duplicated in maize, resulting in closely related PhyA, PhyB, and PhyC gene pairs. These sequences have been mapped to homeologous regions of the maize genome, strongly suggesting that duplication events are a consequence of an ancestral tetraploidization event (GAUT and DOEBLEY 1997). To examine the expression of phytochrome gene family members, we developed a semiquantitative reverse transcriptase polymerase chain reaction (RT-PCR) assay to differentiate between each Phy homeolog. In conjunction with Northern gel blot analysis we show that all phytochrome genes are expressed in several maize tissues. The structure and function of the maize phytochrome gene family is discussed in relation to the evolutionary history of maize and current genome sequencing efforts.
| MATERIALS AND METHODS |
|---|
|
|
|---|
Maize (Z. mays) phytochrome chromophore attachment region and 3'-untranslated region probes:
The 456-bp PhyA chromophore attachment region (CAR) probe was generated from B73 genomic DNA using primers TBp23 (5'-GCTTGCTGCCAAGGCAATCTCCAA-3') and TBp24 (5'-GCTCGGGTTCATCATCCTCTTCA-3') designed to the published PhyA sequence (GenBank accession no. AY260865; CHRISTENSEN and QUAIL 1989). The 1055-bp PhyB CAR probe was generated from B73 genomic DNA with primers 19356F (5'-AACCCGATATGGATCCACTCCAGGGT-3') and 12417R (5'-GGAGTGACACCCAATGGATAGTATTTC-3') designed to a partial maize PhyB sequence (GenBank accession no. AF137332). A 293-bp PhyC CAR probe was generated from B73 genomic DNA using primers TBp29 (5'-TTTCTGTTTATGAAGAACAAAGTG-3') and TBp30 (5'-TTGTACTGTGCATGGCAACCATGGGAAG-3') designed to a partial maize PhyC sequence (GenBank accession no. U61220). PCR was conducted with 2.5 units rTaq (Panvera, Madison, WI), 0.04% DMSO, 500 mM Betaine, 1x PCR buffer. PCR reactions were performed on B73 genomic DNA under the following cycling conditions: 30 sec at 94°, 30 sec at 62°, 4 min at 72°, 35 cycles; 10-min final extension at 72° and 4° soak. PCR products were cloned into pCR4-TOPO vector and transformed into TOP10 Escherichia coli cells (Invitrogen, Carlsbad, CA) according to the manufacturer's protocol to generate the clones pPhyA-CAR (PhyA), pPhyB-CAR (PhyB), and pPhyC-CAR (PhyC). The identity of each clone was confirmed by sequence analysis using T3 or T7 sequencing primers.
Recombinant inbred mapping:
DNA from 94 individuals derived from the Intermated B73 x Mo17 (IBM) 94 (LEE et al. 2002) population were digested with EcoRI (PhyA and PhyB) or HindIII (PhyC) and hybridized to gene-specific probes. Initial hybridizations were performed with CAR probes, which resulted in two segregating polymorphisms. These polymorphisms were scored and submitted to CIMDE (http://www.maizemap.org/CIMDE/cimde.html) to determine the map positions of Phy genes. DNA probes derived from the 3'-untranslated regions (3'-UTRs) of PhyA, PhyB, and PhyC were then used to distinguish homeologs, on the basis of differential hybridization and restriction site polymorphisms. The 3'-UTRs were amplified from B73 genomic DNA for each gene and the products were cloned into pCR4-TOPO according to the manufacturer's specification. The 407-bp PhyA1 3'-UTR product was generated with ZmPA1_3'-UTR-F (5'-CCAGCTCCCTCAGAGTTCATTCGGAAAG-3') and ZmPA1_3'-UTR-R (5'-CCCATATTTAAACTCCACACTACTAAC-3') and then cloned into pCR4-TOPO to generate the clone pPhyA1-3'-UTR. The 356-bp PhyA2 3'-UTR product was generated with ZmPA2_3'-UTR-F (5'-CTCAGAGTTCATTCGGAAAGCCACGGTG-3') and ZmPA2_3'UTR-R (5'-GATGTGAAACTATGAACAAGGAGGCTG-3') and then cloned into pCR4-TOPO to generate the clone pPhyA2-3'-UTR. The 390-bp PhyB1 3'-UTR product was generated with ZmPB1_3'-UTR-F(5'-AGAGGGTGCCTAATTTACGAGAAG-3') and ZmPB1_3'-UTR-R (5'-CAACCGTGTAAGCTGAAACAACTC-3')and then cloned into pCR4-TOPO to generate the clone pPhyB1-3'-UTR. The 313-bp PhyB2 3'-UTR product was generated with ZmPB2_3'-UTR-F (5'-GGTCGATTAGGTAACTTAGGTGGTG3') and ZmPB2_3'-UTR-R (5'-GCCGCCTGTGTACATAACTATCTG-3') and then cloned into pCR4-TOPO to generate the clone pPhyB2-3'-UTR. The 410-bp PhyC1 3'-UTR product was generated with ZmPC1_3'-UTR-F (5'-GAAGTCTCGGTTGTAGCTAGCCAAG-3') and ZmPC1_3'-UTR-R (5'-GAGAGATCAAGAAAACAGACCTACCG-3') and then cloned into pCR4-TOPO to generate the clone pPhyC1-3'-UTR. The 353-bp PhyC2 3'-UTR product was generated with ZmPC2_3'-UTR-F (5'-CCGTGTCAAATTCTCGGGTGTAG-3') and ZmPC2_3'-UTR-R (5'-CAGCTAATCGATCGCGAAGAAAC-3') and then cloned into pCR4-TOPO to generate the clone pPhyC2-3'-UTR.
BAC library screens:
A maize HindIII BAC library, ZmMBBb, was screened to identify all potential Phy genes of maize (TOMKINS et al. 2000, 2002). Radiolabeled CAR probes were generated by Klenow polymerase (Promega, Madison, WI) incorporation of [
-32P]dCTP (New England Nuclear, Boston) using random hexamers (SAMBROOK et al. 1989). Hybridizations were performed as described (CHURCH and GILBERT 1984) at 65° in a bottle incubator (Thermo Hybaid, Needham Heights, MA) using extra-large bottles (300 x 70 mm, Continental Lab Products, San Diego). Probes were purified using Micro Bio-Spin P-30 Tris chromatography columns (Bio-Rad, Hercules, CA) and incorporation efficiency was estimated using a Beckman LS 5801 liquid scintillation counter (Beckman Coulter, Fullerton, CA). Filters were exposed to PhosphorImage screens (Amersham Biosciences, Sunnyvale, CA) overnight and scanned at 50-µm resolution using the Storm 840 scanner (Molecular Dynamics, Sunnyvale, CA). Clone scoring was performed manually aided by ImageQuant (IQMac v1.2) software (Molecular Dynamics). BAC clones were ordered from Clemson University Genomics Institute (Clemson, SC; http://www.genome.clemson.edu/orders/), and single colonies were isolated for overnight culture. Minipreps of putative positive BAC clones were grown in 4.0 ml Luria-Bertani medium (LB) with 12.5 µg ml1 chloramphenicol at 37° with 250-rpm continuous agitation overnight. Cells were pelleted and the BACs were isolated using a standard alkaline lysis protocol (SAMBROOK et al. 1989). DNA isolated from putative PhyA, PhyB, and PhyC clones was digested with restriction enzymes (New England Biolabs, Boston) EcoRI (PhyA) or HindIII (PhyB and PhyC) overnight at 37° and fractionated on 0.8% (w/v) agarose gels. DNA was transferred to Hybond-N+ nylon membranes (Amersham Biosciences). DNA gel blot analysis was performed as described above.
BAC sequencing:
Sequence analysis was performed directly on BAC clones to minimize cloning and Taq-generated errors. A single representative BAC clone served as a template for sequencing each of the six Phy gene family members. Phytochrome gene sequence was generated from the following BAC clones: PhyA1, ZMBBb0024P21; PhyA2, ZMBBb0020P08; PhyB1, ZMBBb0022H07; PhyB2, ZMBBb0001I24; PhyC1, ZMBBb0092I04; and PhyC2, ZMBBb0026K01. Midiprep cultures were generated for the desired BAC clones by direct inoculation from plates to 250 ml LB with 12.5 µg ml1 chloramphenicol. Cultures were grown overnight at 37° at 250 rpm. An alkaline lysis midiprep kit (QIAGEN, Valencia, CA) was used following the protocol for very low copy plasmids and following the modified protocol of KIRSCHNER and STRATAKIS (1999). BigDye terminator mix (Perkin-Elmer, Wellesley, MA) and subsequent purification were performed according to the manufacturer's protocol for BAC templates. Samples were fractionated on an ABI-3100 (Perkin-Elmer), and sequence analysis was performed using Sequencher 4.1.2 (Gene Codes, Ann Arbor, MI) and LaserGene 1.66 software (DNASTAR, Madison, WI). The majority of Phy gene sequence was obtained by primer walking initiated at the central CAR regions of the genes. Most of the Phy gene sequence was obtained with the standard BigDye protocols; however, several regions were refractory to sequence analysis and modifications were made to obtain sequence from CG- and AT-rich regions of the genes. Namely, 0.04% (v/v) DMSO and 500 mM Betaine were substituted for water in the linear amplification stage and eliminated the problem in many cases. Overnight digestion of the template DNA with frequent four-base cutting restriction enzymes (HinfI, MspI, and Sau3AI; New England Biolabs) prior to linear amplification was also utilized to obtain promoter sequence in half of the BAC clones. The accession numbers of the maize phytochrome genes are as follows: PhyA1, AY234825; PhyA2, AY260865 (CHRISTENSEN and QUAIL 1989); PhyB1, AY234827; PhyB2, AY234828; PhyC1, AY234829; and PhyC2, AY234830.
Predicted peptide alignments:
Protein sequence alignments were generated from full-length predicted sequences from Sorghum bicolor, Oryza sativa, and both Z. mays homeologs using MegAlign in LaserGene 1.66 (DNASTAR). MegAlign Clustal W parameters were as follows: Gap Penalty = 10, Gap Extension penalty = 0.2, protein weight matrix was the Gonnet Series. Accession numbers of sequences used are as follows: SbPHYA, AAB41397 SbPHYB, AAB41398 SbPHYC, U56731; OsPHYA, X14172; OsPHYB, X57563; and OsPHYC, AF141942. Domain architecture was defined using the Conserved Domain Architecture Retrieval Tool (CDART; http://www.ncbi.nlm.nih.gov/Structure/lexington/lexington.cgi?cmd=rps) with the default parameters.
RNA gel blot analysis:
Total RNA was extracted from
1 g of tissue as previously described (VAN TUNEN et al. 1988) and quantified using a Beckman DU 530 spectrophotometer. Approximately 20 µg total RNA per lane was fractionated on 6.8% (v/v) formaldehyde, 1.5% (w/v) agarose gel. Hybridizations and image analyses were performed as described for DNA blot analysis (see above).
Semiquantitative RT-PCR and DNA gel blot analysis:
A semiquantitative PCR assay was developed to monitor low abundance Phy gene transcripts and to identify transcripts encoded by homeologous genes. Fifty micrograms total RNA from light- and dark-grown sheath, blade, and mesocotyl tissues was first treated with DNaseI (Invitrogen) at 25° for 1 hr. RNA was purified using phenol:chloroform:isoamyl alcohol (25:24:1, v/v) and precipitated with 95% (v/v) ethanol and 1/10 volume 3 M sodium acetate. The RNA pellet was washed overnight at 4° in 70% (v/v) ethanol to remove excess salts and RNA samples were resuspended in pyrocarbonic acid diethyl ester (DEPC)-treated water at a concentration of 1.0 µg µl1. Five micrograms of total RNA was added to first strand cDNA synthesis reactions performed with SuperScript II kit (Invitrogen) according to the manufacturer's protocol along with "no-RT" controls. After a 2-hr incubation, all samples were treated with RNaseH (Invitrogen) according to the manufacturer's recommended protocol. The cDNA made from each sample was diluted 10-fold in DEPC-treated water and 10 µl was used in PCR amplifications of phytochrome genes (only 5 µl was used for Ubi control). To obtain a semiquantitative estimate of relative transcript levels, PCR reactions were performed for a limited number of cycles. HotStarTaq (QIAGEN) was used according to the suggested PCR protocol: 15-min initial denaturation at 95°; denaturation at 94° for 30 sec, annealing at 55° for 30 sec, 1-min extension at 72° for either 20 or 25 cycles; 10-min final extension at 72°; and 4° soak. A region of the maize Ubi gene was also amplified to control for equal addition of RNA templates. Amplification of Ubi was as described above except that 15 cycles of PCR were performed. The primers are as follows: PhyA-F (5'-GAGAGATCCATGAAGCAAAAGGTTTAC-3'), PhyA-R (5'-GAAGGTTGACATGCCCAGCTTCCCTGAG-3'), PhyB-F (5'-GTTTTGGCTGACTTCGCTAAGCATG-3'), PhyB-R (5'-GGACGATGAGGAAGAAACTCCGCTCTG-3'), PhyC-F (5'-GTCATGCAAGGCATTCCT TTGGGCAAG-3'), PhyC-R (5'-GCTGCTGAGCTGGGCGACCGGGAACTC-3'), UBI-F (5'-CTACAACATTCAGAAGGAGAGCAC-3'), and UBI-R (5'-TCTGCAAGGGTACGGCCATCC-3').
A restriction enzyme digest of PCR products was used to distinguish homeologs. The enzyme MaeIII was used to fractionate PhyA1 from PhyA2 products, AviI was used to distinguish PhyB1 and PhyB2, and AvaII was used to distinguish PhyC1 from PhyC2. Ten microliters of the PCR products was digested overnight according to manufacturer's recommendations (New England Biolabs). Products were fractionated on 2% (w/v) agarose gel and transferred by capillary action in 20x SSC to nylon membranes (Amersham). Hybridization and image analysis were as described above for DNA blot analysis. Bands were visualized with [
-32P]dCTP-labeled RT-PCR products derived from PhyA, PhyB, PhyC, or Ubi-specific amplification products.
| RESULTS |
|---|
|
|
|---|
|
To estimate phytochrome gene copy number, tallies of total Phy-containing clones were made. Six BAC filters were sequentially hybridized with PhyA-, PhyB-, and PhyC-CAR probes with the expectation that
16 clones (2.6 clones x 6 filters) would be recovered for each homeologous gene pair. DNA gel blot analysis of putative phy-containing BAC clones (see MATERIALS AND METHODS) confirmed the identities of 10 PhyA, 7 PhyB, and 10 PhyC clones. Restriction fragment length polymorphism (RFLP) and DNA blot analyses of BAC clones indicated that there were two classes of PhyA-, PhyB-, and PhyC-hybridizing clones. These results are consistent with the DNA gel blot analysis and indicate that duplicated copies of PhyA (PhyA1 and PhyA2), PhyB (PhyB1 and PhyB2), and PhyC (PhyC1 and PhyC2) are present in the B73 genome.
Sequence analysis of BAC clones containing Phy loci:
Full-length genomic sequences for PhyA1, PhyB1, PhyB2, PhyC1, and PhyC2 were obtained by direct sequencing of BAC clones as detailed in MATERIALS AND METHODS. Annotated Phy gene sequences were deposited in GenBank. The full-length PhyA2 sequence was previously reported (CHRISTENSEN and QUAIL 1989). Putative promoter regions and intron/exon junctions were identified through comparisons of previously published Phy gene family members of sorghum (CHILDS et al. 1997) and existing maize EST sequences. Comparison of the predicted genic structure of Phy homeologs reveals a striking conservation of gene organization. Figure 2 shows a schematic representation of the phytochrome gene family in maize. Despite 94% (PHYB and PHYC) to 96% (PHYA) shared amino acid identity between homeologs, several indels and nucleotide substitutions clearly differentiate each Phy gene pair. PhyA homeologs share 94% nucleotide identity across the coding region of the genes and extensive sequence similarity throughout putative 5'- and 3'-UTRs. However, the intron in the 5'-UTR of PhyA1 is approximately double the size of the intron in PhyA2. PhyB homeologs also share 94% sequence identity in coding regions but show greater heterogeneity in intron size. For instance, PhyB1 contains a 1.8-kb retroelement insertion (nucleotides 58767713) in intron II with similarity to a Ty3-gypsy-like retrotransposon. This sequence is not present in intron II of PhyB2 or in PhyB from sorghum (CHILDS et al. 1997), suggesting a relatively recent origin or that PhyB2 is the homolog of sorghum PhyB. Intron III of PhyB1 has also expanded relative to intron III of PhyB2 and PhyB in sorghum (CHILDS et al. 1997), although there is no obvious insertion element to account for this increase in size. PhyC homeologs have similar intron/exon structure and show little similarity in 5'-noncoding regions but not in the 3'-UTRs.
|
|
Despite the striking similarity among all Phy gene homeologs, a number of features differentiate the PhyB gene family homeologs (see Figure S1B at http://www.genetics.org/supplemental/). Several H repeats and a series of G repeats near the amino terminal region of PHYB define a hypervariable region among monocot PhyB sequences. This likely reflects the variable expansion of several microsatellite sequences near the 5' end of the PhyB genes. For instance, a variable number of CAC repeats account for the presence of three to seven H residues near the N terminus of ZmPHYB1, ZmPHYB2, SbPHYB, and OsPHYB. Similarly, slipped mispairing or polymerase stuttering accompanying DNA replication through long stretches of GC residues (FARABAUGH and MILLER 1978) may account for the variable number of G residues near the N terminus of monocot PHYB. The predicted PHYB1 protein is five amino acids shorter than PHYB2 due to additional G repeats at the N-terminal region of the protein and a +1-bp frameshift 9 bases before the predicted stop codon.
RNA gel blot analysis of the phytochrome gene family:
RNA gel blot analysis was used to examine the light regulation of Phy gene transcript accumulation (Figure 4). Due to the strong sequence similarity among homeologous genes, we were unable to develop gene-specific probes that could unambiguously discriminate each homeolog. Thus, the transcript pools examined by Northern analysis reflect the accumulation profiles of homeologous gene pairs and do not address the contributions by single gene family members. Nevertheless, hybridizations with gene-specific probes (see MATERIALS AND METHODS) detected transcripts of
4 kb for PhyA, PhyB, and PhyC, indicating that mature transcripts for PhyA, PhyB, and PhyC accumulate in seedling tissues. As observed previously (CHRISTENSEN and QUAIL 1989), transcripts encoded by PhyA were more abundant in D-grown tissues than in W-grown tissues. PhyB transcripts also accumulated to higher levels in D than in W. Interestingly, PhyC transcript levels were maintained at similar levels under light and dark growth conditions, indicating that the mechanisms regulating PhyC transcript accumulation may differ from those regulating PhyA and PhyB. In addition to light-induced changes in expression, tissue-specific differences in expression were also observed. Most notably, transcripts for PhyB and PhyC accumulated to similar levels in light- and dark-grown mesocotyl tissues and PhyA transcripts were much less abundant in dark-grown mesocotyl tissues relative to dark-grown blade and sheath tissues.
|
300450 bp) for each homeolog were also identical to prevent selective amplification of one or the other homeolog. A limited number of PCR cycles were then performed to enable semiquantitative analysis of gene expression profiles. cDNA products were digested with restriction enzymes that could discriminate between cDNA derived from homeologs due to the presence of single nucleotide polymorphisms (see schematic in Figure 5). Products were fractionated on agarose gels and transferred to nylon membranes and hybridized to gene-specific probes. The relative abundance was estimated from densitometric scans of radiolabeled filters.
|
Gene expression profiles for the PhyB genes were similar to those observed for the PhyA genes. Expression from only one of the PhyB homeologs predominated in all seedling tissues examined, and transcripts for PhyB1 and PhyB2 were more abundant in D-grown relative to W-grown blade and sheath tissues. However, both PhyB1 and PhyB2 were expressed at similar or slightly higher levels in W-grown mesocotyl tissue relative to D-grown mesocotyl tissue.
As observed for PhyA and PhyB homeologs, transcripts from the phyC homeolog located on chromosome 1 predominated in all seedling tissues examined. PhyC transcripts accumulated to slightly higher levels in D-relative to W-grown blade tissue but accumulated to similar levels in D- and W-grown sheath tissues. Together, these studies indicate that phytochrome gene family members are under both tissue- and light-dependent regulation during maize seedling development.
| DISCUSSION |
|---|
|
|
|---|
Although the function of any Phy gene has yet to be determined in maize, it is tempting to speculate on some of the consequences of maintaining such similar homeologs. Perhaps one of the most intriguing possibilities is that the homeologs could form active heterodimers that mediate responses different from those regulated by homodimers. Early studies of oat phytochrome proteins indicated that the biologically active forms of phytochrome function as homodimers (JONES and QUAIL 1986). Evidence strongly suggests that these dimerized forms interact with nuclear localized transcription factors to mediate phytochrome responses (SAKAMOTO and NAGATANI 1996; KIRCHER et al. 1999, 2002; YAMAGUCHI et al. 1999; HUQ and QUAIL 2002; HUQ et al. 2003). In the case of Arabidopsis phyB, dimerized and nuclear localized N-terminal sequences (amino acids 1651) are sufficient to mediate a number of red light responses in a phyB mutant background (MATSUSHITA et al. 2003). Although many regions of the mature protein assist or strengthen dimerization, the amino acids that are critical for dimerization of Arabidopsis phyA lie between 1116 and 1161 (QUAIL 1997). Due to the high degree of sequence identity between homeologous gene pairs, especially within this dimerization domain, it is likely that the maize PHY homeologs will also form heterodimers. Clearly, the consequences of heterodimerization are dependent on the, as yet, untested functional divergence of each protein. Nevertheless, the overlapping expression domains of maize Phy homeologs are consistent with the possibility of direct interaction of Phy homeolog proteins.
The high sequence conservation between homeologous phy gene pairs also highlights some of the challenges that face maize genome sequencing projects. Searches of The Institute for Genome Research (TIGR) Maize Gene Index (ZmGI) revealed tentative consensus (TC) sequences for PhyA, PhyB, and PhyC genes of maize. This database curates assemblies of maize ESTs (TIGR contigs or TCs) and identified two potential PhyA genes. For example, TC185264 incorporates several EST sequences derived from multiple cDNA libraries and includes the previously published PhyA2 sequence (CHRISTENSEN and QUAIL 1989). Comparisons of the individual EST clone sequences that constitute TC185264 with the PhyA1 and PhyA2 sequences reported here reveal that both PhyA1 and PhyA2 sequences are represented in the assembly. A tentative consensus sequence (TC171984) was also derived from ESTs that were all most similar to PhyA1 sequence. Thus, TC185264 is most similar to PhyA2 whereas TC171984 is most similar to PhyA1 sequence, but neither assembly is identical to the predicted PhyA1 gene reported here or to the previously published PhyA2 sequence (CHRISTENSEN and QUAIL 1989).
Attempts to distinguish highly similar gene homeologs such as the phytochrome genes are further complicated by the presence of repetitive elements within genes. As shown in Figure 2, a nearly 2-kb retroelement insertion is present in the second intron of PhyB1. Although gene islands are anticipated to be scattered throughout a sea of retrotransposons in maize (SANMIGUEL et al. 1996), the presence of large retroposon sequences within genes could prove problematic in genomic assemblies (RABINOWICZ et al. 2003). To enrich for gene-rich regions of the genome, current maize genome sequencing projects are utilizing methyl filtration technologies and high Cot sequences to reduce representation of repetitive LTR retroelements (RABINOWICZ et al. 2003). Since we did not examine the methylation status of this element, it is difficult to predict whether or not PhyB1 intron sequences would be represented in methylation-filtered genomic libraries (RABINOWICZ et al. 1999). Regardless, the repetitive nature of the element and its large size would likely lead to its exclusion from genomic assemblies. Recent methods developed to link gene-rich regions should prove particularly useful in assembling genes such as PhyB1 in maize (YUAN et al. 2002).
In addition to confounding bioinformatics-based approaches to genome annotation, the duplications of the maize Phy gene family also highlight challenges faced in defining syntenic relationships among grass genomes. The relatively recent polyploidization in the maize lineage has likely contributed to the duplication of >80% of the maize genome as colinear regions of homeologous maize chromosomes (GAUT 2001). These highly similar sequences have the potential to greatly complicate the analysis of orthologous relationships in the grasses. For instance, in a recent study, a PhyA-containing BAC clone of maize was characterized as a syntenic region with a PhyA gene sequence from sorghum (MORISHIGE et al. 2002). The authors characterized a single maize BAC clone and assigned this gene and all predicted gene sequences residing on this BAC to chromosome 1L, corresponding to the published map position of PhyA2 (CHRISTENSEN and QUAIL 1989). However, as we report here, the original assignment of this homeolog to chromosome 1 was incorrect and is instead located on chromosome 5S. Furthermore, because only one of the two PhyA orthologs was characterized, the previously reported synteny comparisons made between the sorghum and maize PhyA regions are incomplete (MORISHIGE et al. 2002). It is tempting to speculate that because all Phy genes located on chromosome 1 are more highly expressed in seedling tissues than are their homeologs located on chromosomes 5 and 9, genes tightly linked to these Phy genes may show similar homeolog-specific biases in gene expression.
An analysis of Phy transcript profiles was conducted to examine potential light- and tissue-specific regulation of Phy gene activity in maize. Transcripts for both PhyA1 and PhyA2 accumulated to much higher levels in dark relative to white light in blade, sheath, and mesocotyl tissues. This light-dependent decrease in phyA transcripts has been observed in several grass species (COLBERT et al. 1983; HERSHEY et al. 1984; CHRISTENSEN and QUAIL 1989; KAY et al. 1989; DEHESH et al. 1991) and likely reflects a common regulatory mechanism that acts to maintain high levels of phyA transcript in the dark (BRUCE et al. 1989). Furthermore, several elements have been defined in monocot PhyA promoters that likely mediate a high transcription rate of PhyA in the dark (BRUCE and QUAIL 1990; BRUCE et al. 1991). The presence of PE3 and RE1 elements in the promoter region of maize PhyA homeologs suggests a similar mechanism of transcriptional control.
Despite the common features shared by monocot phyA promoters, the kinetics of light-mediated reductions of PhyA transcript levels differs among grass species (COLBERT et al. 1985; LISSEMORE and QUAIL 1988; CHRISTENSEN and QUAIL 1989; KAY et al. 1989). In oat (LISSEMORE and QUAIL 1988) and rice (KAY et al. 1989) there is a rapid decrease in PhyA transcript levels that can largely be attributed to a decreased rate of PhyA transcription. In maize the reduction in PhyA transcript accumulation is less rapid (CHRISTENSEN and QUAIL 1989) and indicates that post-transcriptional regulation may also contribute to PhyA transcript accumulation. Here we have shown that this light-dependent reduction in PhyA transcript accumulation is intrinsic to PhyA1 and PhyA2 regulation and operates in several seedling tissues.
Although little is known of PhyB or PhyC gene regulation in monocots, detailed characterizations of organ- and tissue-specific regulation of phytochrome gene expression have been performed in several model eudicots. In Arabidopsis, PHYA, PHYB, PHYC, PHYD, and PHYE transcripts accumulate in several plant organs (CLACK et al. 1994), and promoter-Gus fusions indicate that the promoters determine tissue specificity, light induciblity, and magnitude of gene expression (SOMERS and QUAIL 1995a,b; GOOSEY et al. 1997). However, translational and/or post-translational mechanisms largely determine the accumulation of phytochrome proteins in Arabidopsis (CLACK et al. 1994; SHARROCK and CLACK 2002; SHARROCK et al. 2003b). In potato, a uniform pattern of PhyB expression was observed in several mature plant organs (HEYER and GATZ 1992). In tomato, however, temporal tissue- and organ-specific differences were observed for all phytochrome genes (HAUSER et al. 1997, 1998). Also, in general, PhyA transcripts are most abundant in several seedling and mature plant tissues, with the exception of ripening fruit, where PhyB2 and PhyF transcripts predominate (HAUSER et al. 1997). In tobacco, PHYB1-Gus promoter fusions also suggest tissue-specific regulation of PHYB1 transcription in developing seedlings and mature plants. In particular, PHYB1 expression appeared to be restricted to phloem cells of the stem, petiole, and midrib tissues (ADAM et al. 1996). The organ- and tissue-specific expression profiles of PHYB genes in tomato (HAUSER et al. 1997) and tobacco (ADAM et al. 1996) and maize (this study) contrast with the relatively uniform pattern of PHYB expression in Arabidopsis (SOMERS and QUAIL 1995b) and potato (HEYER and GATZ 1992) and suggest that there may be little conservation in the transcriptional regulation of PHYB and PHYC accumulation between monocots and dicots.
In rice (DEHESH et al. 1991; BASU et al. 2000) and sorghum (CHILDS et al. 1997) PhyB and PhyC transcripts are maintained at relatively constant levels under dark and white light conditions, but nothing is known of tissue- or organ-specific expression profiles. In maize, we have shown that tissue-specific accumulation of PhyB is light dependent. That is, PhyB transcripts are most abundant in blade and sheath tissues of dark-grown plants but accumulate to similar levels in mesocotyl tissues of light- and dark-grown plants (see Figure 4). Thus, light appears to play an important role in regulating the organ-specific patterns of PhyB gene expression in maize. PhyC transcripts were the least abundant of all Phy genes in dark-grown tissues and were maintained at low levels in light-grown tissues as well.
The identification and sequence analysis of all phytochrome gene family members reported here will now permit the utilization of association methods to examine the role of phytochrome gene function in this important crop plant. In addition, the elucidation of phytochrome gene family organization in maize will guide reverse genetic screens to disrupt all phytochrome gene family members in maize. The single, double, and higher order mutants will provide valuable insights into the functional divergence of the phytochrome gene family in maize and the role of phytochrome in the regulation of agronomically important traits.
| ACKNOWLEDGEMENTS |
|---|
|
|
|---|
| FOOTNOTES |
|---|
| LITERATURE CITED |
|---|
|
|
|---|
ADAM, E., L. KOZMA-BOGNAR, C. KOLAR, E. SCHAFER and F. NAGY, 1996 The tissue-specific expression of a tobacco phytochrome B gene. Plant Physiol. 110: 10811088.[Abstract]
ALBA, R., P. M. KELMENSON, M. M. CORDONNIER-PRATT and L. H. PRATT, 2000 The phytochrome gene family in tomato and the rapid differential evolution of this family in angiosperms. Mol. Biol. Evol. 17: 362373.
ARAVIND, L., and C. P. PONTING, 1997 The GAF domain: an evolutionary link between diverse phototransducing proteins. Trends Biochem. Sci. 22: 458459.[CrossRef][Medline]
AUKERMAN, M. J., M. HIRSCHFELD, L. WESTER, M. WEAVER, T. CLACK et al., 1997 A deletion in the PHYD gene of the Arabidopsis Wassilewskija ecotype defines a role for phytochrome D in red/far-red light sensing. Plant Cell 9: 13171326.[Abstract]
BALLARE, C. L., and J. J. CASAL, 2000 Light signals perceived by crop and weed plants. Field Crops Res. 67: 149160.[CrossRef]
BASU, D., K. DEHESH, H. J. SCHNEIDER-POETSCH, S. E. HARRINGTON, S. R. MCCOUCH et al., 2000 Rice PHYC gene: structure, expression, map position and evolution. Plant Mol. Biol. 44: 2742.[CrossRef][Medline]
BRUCE, W. B., and P. H. QUAIL, 1990 cis-acting elements involved in photoregulation of an oat phytochrome promoter in rice. Plant Cell 2: 10811089.
BRUCE, W. B., A. H. CHRISTENSEN, T. KLEIN, M. FROMM and P. H. QUAIL, 1989 Photoregulation of a phytochrome gene promoter from oat transferred into rice by particle bombardment. Proc. Natl. Acad. Sci. USA 86: 96929696.
BRUCE, W. B., X. W. DENG and P. H. QUAIL, 1991 A negatively acting DNA sequence element mediates phytochrome-directed repression of phyA gene transcription. EMBO J. 10: 30153024.[Medline]
CASAL, J. J., L. G. LUCCIONI, K. A. OLIVERIO and H. E. BOCCALANDRO, 2003 Light, phytochrome signalling and photomorphogenesis in Arabidopsis. Photochem. Photobiol. Sci. 2: 625636.[CrossRef][Medline]
CHILDS, K. L., F. R. MILLER, M. M. CORDONNIER-PRATT, L. H. PRATT, P. W. MORGAN et al., 1997 The sorghum photoperiod sensitivity gene, Ma3, encodes a phytochrome B. Plant Physiol. 113: 611619.[Abstract]
CHRISTENSEN, A. H., and P. H. QUAIL, 1989 Structure and expression of a maize phytochrome-encoding gene. Gene 85: 381390.[CrossRef][Medline]
CHURCH, G. M., and W. GILBERT, 1984 Genomic sequencing. Proc. Natl. Acad. Sci. USA 81: 19911995.
CLACK, T., S. MATHEWS and R. A. SHARROCK, 1994 The phytochrome apoprotein family in Arabidopsis is encoded by five genes: the sequences and expression of PHYD and PHYE. Plant Mol. Biol. 25: 413427.[CrossRef][Medline]
COLBERT, J. T., H. P. HERSHEY and P. H. QUAIL, 1983 Autoregulatory control of translatable phytochrome mRNA levels. Proc. Natl. Acad. Sci. USA 80: 22482252.
COLBERT, J. T., H. P. HERSHEY and P. H. QUAIL, 1985 Phytochrome regulation of phytochrome mRNA abundance. Plant Mol. Biol. 5: 91101.
DEHESH, K., J. TEPPERMAN, A. H. CHRISTENSEN and P. H. QUAIL, 1991 phyB is evolutionarily conserved and constitutively expressed in rice seedling shoots. Mol. Gen. Genet. 225: 305313.[CrossRef][Medline]
DEHESH, K., C. FRANCI, B. M. PARKS, K. A. SEELEY, T. W. SHORT et al., 1993 Arabidopsis HY8 locus encodes phytochrome A. Plant Cell 5: 10811088.
DEVLIN, P. F., S. R. PATEL and G. C. WHITELAM, 1998 Phytochrome E influences internode elongation and flowering time in Arabidopsis. Plant Cell 10: 14791487.
DEVLIN, P. F., P. R. ROBSON, S. R. PATEL, L. GOOSEY, R. A. SHARROCK et al., 1999 Phytochrome D acts in the shade-avoidance syndrome in Arabidopsis by controlling elongation growth and flowering time. Plant Physiol. 119: 909915.
FARABAUGH, P. J., and J. H. MILLER, 1978 Genetic studies of lac repressor VII. On the molecular nature of spontaneous hotspots in the lacI gene of Escherichia coli. J. Mol. Biol. 126: 847863.[CrossRef][Medline]
FRANKLIN, K. A., S. J. DAVIS, W. M. STODDART, R. D. VIERSTRA and G. C. WHITELAM, 2003 Mutant analyses define multiple roles for phytochrome C in Arabidopsis photomorphogenesis. Plant Cell 15: 19811989.
GAUT, B. S., 2001 Patterns of chromosomal duplication in maize and their implications for comparative maps of the grasses. Genome Res. 11: 5566.
GAUT, B. S., and J. F. DOEBLEY, 1997 DNA sequence evidence for the segmental allotetraploid origin of maize. Proc. Natl. Acad. Sci. USA 94: 68096814.
GOOSEY, L., L. PALECANDA and R. A. SHARROCK, 1997 Differential patterns of expression of the Arabidopsis PHYB, PHYD, and PHYE phytochrome genes. Plant Physiol. 115: 959969.[Abstract]
GYULA, P., E. SCHAFER and F. NAGY, 2003 Light perception and signalling in higher plants. Curr. Opin. Plant Biol. 6: 446452.[CrossRef][Medline]
HANUMAPPA, M., L. H. PRATT, M. M. CORDONNIER-PRATT and G. F. DEITZER, 1999 A photoperiod-insensitive barley line contains a light-labile phytochrome B. Plant Physiol. 119: 10331040.
HAUSER, B. A., L. H. PRATT and M. M. CORDONNIER-PRATT, 1997 Absolute quantification of five phytochrome transcripts in seedlings and mature plants of tomato (Solanum lycopersicum L.). Planta 201: 379387.[CrossRef][Medline]
HAUSER, B. A., M. M. CORDONNIER-PRATT and L. H. PRATT, 1998 Temporal and photoregulated expression of five tomato phytochrome genes. Plant J. 14: 431439.[CrossRef][Medline]
HERSHEY, H. P., J. T. COLBERT, J. L. LISSEMORE, R. F. BARKER and P. H. QUAIL, 1984 Molecular cloning of cDNA for Avena phytochrome. Proc. Natl. Acad. Sci. USA 81: 23322336.
HEYER, A., and C. GATZ, 1992 Isolation and characterization of a cDNA-clone coding for potato type B phytochrome. Plant Mol. Biol. 20: 589600.[CrossRef][Medline]
HO, Y. S., L. M. BURDEN and J. H. HURLEY, 2000 Structure of the GAF domain, a ubiquitous signaling motif and a new class of cyclic GMP receptor. EMBO J. 19: 52885299.[CrossRef][Medline]
HOWE, G. T., P. A. BUCCIAGLIA, W. P. HACKETT, G. R. FURNIER, M. M. CORDONNIER-PRATT et al., 1998 Evidence that the phytochrome gene family in black cottonwood has one PHYA locus and two PHYB loci but lacks members of the PHYC/F and PHYE subfamilies. Mol. Biol. Evol. 15: 160175.[Abstract]
HUQ, E., and P. H. QUAIL, 2002 PIF4, a phytochrome-interacting bHLH factor, functions as a negative regulator of phytochrome B signaling in Arabidopsis. EMBO J. 21: 24412450.[CrossRef][Medline]
HUQ, E., B. AL-SADY and P. H. QUAIL, 2003 Nuclear translocation of the photoreceptor phytochrome B is necessary for its biological function in seedling photomorphogenesis. Plant J. 35: 660664.[CrossRef][Medline]
IZAWA, T., T. OIKAWA, S. TOKUTOMI, K. OKUNO and K. SHIMAMOTO, 2000 Phytochromes confer the photoperiodic control of flowering in rice (a short-day plant). Plant J. 22: 391399.[CrossRef][Medline]
JONES, A. M., and P. H. QUAIL, 1986 Quaternary structure of 124 kiloDalton phytochrome from Avena sativa. Biochemistry 25: 29872995.[CrossRef]
KAY, S. A., B. KEITH, K. SHINOZAKI, M. L. CHYE and N. H. CHUA, 1989 The rice phytochrome gene: structure, autoregulated expression, and binding of GT-1 to a conserved site in the 5' upstream region. Plant Cell 1: 351360.
KIRCHER, S., L. KOZMA-BOGNAR, L. KIM, E. ADAM, K. HARTER et al., 1999 Light quality-dependent nuclear import of the plant photoreceptors phytochrome A and B. Plant Cell 11: 14451456.
KIRCHER, S., P. GIL, L. KOZMA-BOGNAR, E. FEJES, V. SPETH et al., 2002 Nucleocytoplasmic partitioning of the plant photoreceptors phytochrome A, B, C, D, and E is regulated differentially by light and exhibits a diurnal rhythm. Plant Cell 14: 15411555.
KIRSCHNER, L. S., and C. A. STRATAKIS, 1999 Large-scale preparation of sequence-ready bacterial artificial chromosome DNA using QIAGEN columns. Biotechniques 27: 7274.[Medline]
KRALL, L., and J. W. REED, 2000 The histidine kinase-related domain participates in phytochrome B function but is dispensable. Proc. Natl. Acad. Sci. USA 97: 81698174.
LEE, M., N. SHAROPOVA, W. D. BEAVIS, D. GRANT, M. KATT et al., 2002 Expanding the genetic map of maize with the intermated B73 x Mo17 (IBM) population. Plant Mol. Biol. 48: 453461.[CrossRef][Medline]
LI, W.-Z., and C.-C. CHINNAPPA, 2003 The phytochrome gene family in the Stellaria longipes complex. Int. J. Plant Sci. 164: 657673.[CrossRef]
LISSEMORE, J. L., and P. H. QUAIL, 1988 Rapid transcriptional regulation by phytochrome of the genes for phytochrome and chlorophyll a/b-binding protein in Avena sativa. Mol. Cell. Biol. 8: 48404850.
MARKELZ, N. H., D. E. COSTICH and T. P. BRUTNELL, 2003 Photomorphogenic responses in maize seedling development. Plant Physiol. 133: 15781591.
MATHEWS, S., and R. A. SHARROCK, 1996 The phytochrome gene family in grasses (Poaceae): a phylogeny and evidence that grasses have a subset of the loci found in dicot angiosperms. Mol. Biol. Evol. 13: 11411150.[Abstract]
MATHEWS, S., and R. A. SHARROCK, 1997 Phytochrome gene diversity. Plant Cell Environ. 20: 666671.[CrossRef]
MATHEWS, S., M. LAVIN and R. A. SHARROCK, 1995 Evolution of the phytochrome gene family in grasses (Poaceae): a phylogeny and evidence that grasses have a subset of the loci found in dicot angiosperms. Mol. Evol. Biol. 13: 11411150.
MATSUSHITA, T., N. MOCHIZUKI and A. NAGATANI, 2003 Dimers of the N-terminal domain of phytochrome B are functional in the nucleus. Nature 424: 571574.[CrossRef][Medline]
MONTE, E., J. M. ALONSO, J. R. ECKER, Y. ZHANG, X. LI et al., 2003 Isolation and characterization of phyC mutants in Arabidopsis reveals complex crosstalk between phytochrome signaling pathways. Plant Cell 15: 19621980.
MORISHIGE, D. T., K. L. CHILDS, L. D. MOORE and J. E. MULLET, 2002 Targeted analysis of orthologous phytochrome A regions of the sorghum, maize, and rice genomes using comparative gene-island sequencing. Plant Physiol. 130: 16141625.
NAGATANI, A., J. W. REED and J. CHORY, 1993 Isolation and initial characterization of Arabidopsis mutants that are deficient in phytochrome A. Plant Physiol. 102: 269277.[Abstract]
NEFF, M. M., C. FANKHAUSER and J. CHORY, 2000 Light: an indicator of time and place. Genes Dev. 14: 257271.
PARKS, B. M., and P. H. QUAIL, 1993 hy8, a new class of Arabidopsis long hypocotyl mutants deficient in functional phytochrome A. Plant Cell 5: 3948.
PELLEQUER, J. L., K. A. WAGER-SMITH, S. A. KAY and E. D. GETZOFF, 1998 Photoactive yellow protein: a structural prototype for the three-dimensional fold of the PAS domain superfamily. Proc. Natl. Acad. Sci. USA 95: 58845890.
PONTING, C. P., and L. ARAVIND, 1997 PAS: a multifunctional domain family comes to light. Curr. Biol. 7: R674R677.[Medline]
QUAIL, P. H., 1997 An emerging molecular map of the phytochromes. Plant Cell Environ. 20: 657665.[CrossRef]
QUAIL, P. H., 2002 Photosensory perception and signalling in plant