Department of Biology, Massachusetts Institute of Technology, Cambridge, Massachusetts c

Transcription

1 The Plant Cell, Vol. 17, , May 2005, ª 2005 American Society of Plant Biologists MicroRNA-Directed Regulation of Arabidopsis AUXIN RESPONSE FACTOR17 Is Essential for Proper Development and Modulates Expression of Early Auxin Response Genes W Allison C. Mallory, a David P. Bartel, a,b,1 and Bonnie Bartel c,1 a Whitehead Institute for Biomedical Research, Cambridge, Massachusetts b Department of Biology, Massachusetts Institute of Technology, Cambridge, Massachusetts c Department of Biochemistry and Cell Biology, Rice University, Houston, Texas The phytohormone auxin plays critical roles during plant growth, many of which are mediated by the auxin response transcription factor (ARF) family. MicroRNAs (mirnas), endogenous 21-nucleotide riboregulators, target several mrnas implicated in auxin responses. mir160 targets ARF10, ARF16, and ARF17, three of the 23 Arabidopsis thaliana ARF genes. Here, we describe roles of mir160-directed ARF17 posttranscriptional regulation. Plants expressing a mirna-resistant version of ARF17 have increased ARF17 mrna levels and altered accumulation of auxin-inducible GH3-like mrnas, YDK1/GH3.2, GH3.3, GH3.5, and DFL1/GH3.6, which encode auxin-conjugating proteins. These expression changes correlate with dramatic developmental defects, including embryo and emerging leaf symmetry anomalies, leaf shape defects, premature inflorescence development, altered phyllotaxy along the stem, reduced petal size, abnormal stamens, sterility, and root growth defects. These defects demonstrate the importance of mir160-directed ARF17 regulation and implicate ARF17 as a regulator of GH3-like early auxin response genes. Many of these defects resemble phenotypes previously observed in plants expressing viral suppressors of RNA silencing and plants with mutations in genes important for mirna biogenesis or function, providing a molecular rationale for phenotypes previously associated with more general disruptions of mirna function. INTRODUCTION MicroRNAs (mirnas) are endogenous 21-nucleotide riboregulators that modulate gene expression in plants and animals (Bartel and Bartel, 2003; Carrington and Ambros, 2003; Bartel, 2004; Mallory and Vaucheret, 2004). In plants, mirnas are processed from imperfectly complementary stem-loop precursor RNAs (Reinhart et al., 2002), and proper accumulation of some mirnas depends on the nuclear activity of DICER-LIKE1 (DCL1) as well as HYL1, HEN1, and AGO1 (Park et al., 2002; Reinhart et al., 2002; Boutet et al., 2003; Kasschau et al., 2003; Han et al., 2004; Vaucheret et al., 2004; Vazquez et al., 2004a). In both plants and animals, mirnas can direct cleavage of target mrnas and/or inhibition of productive translation, although plant mirnas appear to act primarily via mrna cleavage (Lee et al., 1993; Wightman et al., 1993; Olsen and Ambros, 1999; Hutvágner and Zamore, 2002; Llave et al., 2002; Aukerman and Sakai, 2003; Kasschau et al., 2003; Tang et al., 2003; Chen, 2004; Yekta et al., 2004). Studies examining pairing requirements between mir165/166:phabulosa and mir171:scl6 mirna: 1 To whom correspondence should be addressed. bartel@rice. edu or dbartel@wi.mit.edu; fax or The authors responsible for distribution of materials integral to the findings presented in this article in accordance with the policy described in the Instructions for Authors ( are: David P. Bartel (dbartel@wi.mit.edu) and Bonnie Bartel (bartel@rice.edu). W Online version contains Web-only data. Article, publication date, and citation information can be found at target pairs have revealed that mismatches in the mirna 59 region disrupt cleavage more effectively than mismatches in the mirna 39 region (Mallory et al., 2004b; Parizotto et al., 2004). Many plant and animal mirnas do not accumulate ubiquitously, but instead are restricted to specific tissue types and developmental stages, suggesting spatial and temporal regulation of mirna accumulation and, thus, of target expression (Carrington and Ambros, 2003; Ambros, 2004; Bartel, 2004). Many mirnas isolated from the dicot Arabidopsis thaliana are conserved in the monocot rice (Oryza sativa) and in other plants, implying conserved evolutionary roles for plant mirnas (Reinhart et al., 2002; Floyd and Bowman, 2004; Jones-Rhoades and Bartel, 2004; Sunkar and Zhu, 2004; Wang et al., 2004a; Axtell and Bartel, 2005). Plant mirnas often have extensive, evolutionarily conserved complementarity to plant mrnas (Rhoades et al., 2002; Axtell and Bartel, 2005). This observation has enabled numerous regulatory targets to be confidently predicted (Park et al., 2002; Rhoades et al., 2002; Xie et al., 2003; Jones- Rhoades and Bartel, 2004; Sunkar and Zhu, 2004; Wang et al., 2004b), more than 50 of which have been experimentally validated in plants (Dugas and Bartel, 2004). Plant mirna complementary sites are usually present as single copies in the open reading frame of the target mrna, although complementary sites in both 59 and 39 untranslated regions (UTRs) have been reported. This trend contrasts with the initially characterized animal mirna complementary sites, which typically are present as multiple sites in the 39 UTR of target mrnas (Bartel, 2004). However, recent studies indicate that targeting of single sites and targeting within open reading frames are

2 mir160 Regulates Auxin Response Factors 1361 widespread for mammalian mirnas (Lewis et al., 2005; Lim et al., 2005). Many plant mirna targets encode transcription factors involved in cell fate determination (Rhoades et al., 2002), supporting the idea that mirnas regulate plant development. In further support of this idea, plants impaired in mirna accumulation, such as dcl1, ago1, hyl1, and hen1 mutants (Park et al., 2002; Reinhart et al., 2002; Boutet et al., 2003; Kasschau et al., 2003; Han et al., 2004; Vaucheret et al., 2004; Vazquez et al., 2004a), and plants expressing viral suppressors of gene silencing that alter mirna accumulation (Mallory et al., 2002; Kasschau et al., 2003; Chapman et al., 2004; Chen et al., 2004; Dunoyer et al., 2004) display dramatic anomalies during vegetative and floral development. The overexpression of mirnas and the expression of mirnaresistant targets in vivo have allowed assignment of developmental roles to mir159/mir319/jaw, mir164, mir165/166, and mir172 mirna families. These mirnas regulate rosette leaf expansion and curvature (Palatnik et al., 2003), embryonic, vegetative, and floral organ boundary formation (Laufs et al., 2004; Mallory et al., 2004a), radial patterning (Emery et al., 2003; Juarez et al., 2004; Mallory et al., 2004b; McHale and Koning, 2004; Zhong and Ye, 2004), and floral organ identity and flowering time (Aukerman and Sakai, 2003; Chen, 2004). mir162 and mir168 likely influence development through feedback regulation of mirna pathway components, DCL1 and AGO1, respectively (Xie et al., 2003; Vaucheret et al., 2004). Particular mirnas accumulate in response to sulfur starvation (Jones- Rhoades and Bartel, 2004), abiotic stresses (Sunkar and Zhu, 2004), or phytohormones (Achard et al., 2004; Sunkar and Zhu, 2004), suggesting that mirnas are important not only for development but also for responses to environmental stimuli. The phytohormone auxin influences many aspects of plant development (Rogg and Bartel, 2001; Liscum and Reed, 2002; Swarup et al., 2002; Friml, 2003; Jürgens, 2003), and the identity of several mirna targets suggests roles for mirnas in auxin signaling. mir160 is complementary to AUXIN RESPONSE FACTOR10 (ARF10), ARF16, and ARF17 (Rhoades et al., 2002), and mir167 is complementary to ARF6 and ARF8 (Rhoades et al., 2002; Bartel and Bartel, 2003). Thus, at least five of the 23 Arabidopsis ARF transcription factors are potentially mirna regulated. ARFs are a plant-specific family of DNA binding proteins that control auxin-regulated transcription (Guilfoyle et al., 1998). They bind to auxin-responsive promoter elements (AuxREs), which are found in early auxin response genes, including Auxin/Indole-3-Acetic Acid (Aux/IAA), SAUR, and GH3, and can either enhance or repress transcription (Abel and Theologis, 1996; Ulmasov et al., 1999a, 1999b; Hagen and Guilfoyle, 2002). Most ARFs have a conserved N-terminal DNA binding domain (DBD), a nonconserved middle region conferring transcriptional repression or activation, and a conserved C-terminal dimerization domain (CTD) that can mediate ARF homodimerization and heterodimerization with Aux/IAA repressors (Guilfoyle et al., 1998; Ulmasov et al., 1999a, 1999b; Guilfoyle and Hagen, 2001; Liscum and Reed, 2002; Tiwari et al., 2003). Aux/IAA proteins, short-lived nuclear proteins, can heterodimerize with activating ARF proteins, preventing early auxinresponse gene expression (Kim et al., 1997; Rouse et al., 1998; Reed, 2001; Dharmasiri and Estelle, 2002; Kepinski and Leyser, 2002; Liscum and Reed, 2002; Tiwari et al., 2004). C-terminal domains of Aux/IAA proteins mediate heterodimerization and are conserved with the CTD of most ARF proteins. Increased auxin levels accelerate proteolysis of Aux/IAA proteins, which would allow ARF proteins to homodimerize and impose their regulatory functions on early auxin-response gene expression. The transcriptional repressors ARF1-ARF4 and ARF9 have Pro-Ser-Thr rich middle regions, whereas the transcriptional activators ARF5-ARF8 have Gln-Leu-Ser rich middle regions (Ulmasov et al., 1999b; Tiwari et al., 2003). These observations suggest that ARFs can be classified as activators or repressors based on the amino acid composition of their middle regions. In addition, most tested ARFs require CTD-mediated dimerization to bind stably to AuxREs in vitro (Ulmasov et al., 1997, 1999a); however, ARF1 does not require the CTD to bind AuxREs in vitro (Ulmasov et al., 1997, 1999a), and ARF3 lacks a CTD. ARF10 and ARF16 display typical ARF sequence characteristics: each has a conserved DBD and CTD and a nonconserved middle region. By contrast, ARF17 is unusual because it lacks a conserved CTD, which is present in 21 of 23 Arabidopsis ARF proteins (Ulmasov et al., 1997, 1999a; Hagen and Guilfoyle, 2002). mir160 and target ARFs are conserved in dicots and monocots (Rhoades et al., 2002; Bartel and Bartel, 2003), but the importance of mir160-directed regulation of ARF17, ARF16, and ARF10 has not been explored. These three ARFs comprise a subgroup of Arabidopsis ARF proteins (Remington et al., 2004). arf10 and arf16 loss-of-function mutants do not display obvious developmental anomalies (Okushima et al., 2005), and ARF17 has not been characterized in vivo. Here, we demonstrate that disrupting mir160 regulation of ARF17 increases ARF17 mrna levels, leads to severe developmental abnormalities, including defects in embryonic, root, vegetative, and floral development, and alters GH3-like gene expression. These results indicate that mir160-directed regulation is critical for the developmental functions of ARF17 and expose ARF17 as a possible transcriptional regulator of GH3-like early auxin-response genes. RESULTS mir160 Regulates ARF10, ARF16,andARF17 Expression We detected fragments of ARF10, ARF16, and ARF17 mrnas from Arabidopsis tissues and mapped the end of each fragment to precisely the nucleotide predicted for mir160-directed cleavage (Figure 1A) using a modified form of rapid amplification of cdna ends (59 RACE) (Llave et al., 2002). This analysis confirmed what had been reported for ARF10 and ARF17 (Kasschau et al., 2003) and indicated that ARF16 also is posttranscriptionally regulated by mir160. To determine the tissues in which mir160 regulates ARF10, ARF16, and ARF17 mrnas, we assayed RNA accumulation in a variety of wild-type Arabidopsis tissues using gel blot analysis. We found mir160 and ARF10, ARF16, and ARF17 full-length mrnas in all tissues analyzed (Figure 1B), indicating that these genes are widely expressed. In addition to full-length transcripts, mrnas the size of predicted ARF10, ARF16, and ARF17 39 cleavage products were detected in all

3 1362 The Plant Cell tissues analyzed (Figure 1B), suggesting that mir160 is posttranscriptionally regulating these messages in many cells of the plant. The accumulation of mir160 and ARF10, ARF16, and ARF17 39 cleavage products was highest in roots (Figure 1B), but there was no simple correlation among the level of mir160 and full-length ARF mrna accumulation in the tissues analyzed (Figure 1B). The relative abundance of ARF17 39 cleavage product compared with full-length ARF17 mrna in many tissues (Figure 1B) suggested extensive mirna-directed cleavage of this transcript and prompted us to explore ARF17 regulation in more detail. 35S:ARF17 Plants Overaccumulate ARF17 39 Cleavage Product To assess the importance of ARF17 transcriptional regulation during development, we expressed an ARF17 genomic clone in wild-type Arabidopsis under the control of the 35S promoter of Cauliflower mosaic virus (CaMV), which drives strong constitutive expression (Cary et al., 2002). Among 23 35S:ARF17 primary transformants, 21 displayed no obvious developmental defects, whereas two displayed a slight reduction in rosette leaf size and overall stature (data not shown). RNA gel blot analysis revealed a threefold to fourfold increase in full-length ARF17 mrna accumulation in rosette leaves of four out of five 35S:ARF17 primary transformants, including the two smaller plants Figure 1. mir160 Regulates ARF10, ARF16, and ARF17 in Vivo. (A) mir160 cleavage sites in ARF10, ARF16, andarf17 mrnas determined by RNA ligase-mediated 59 RACE. The agarose gel shows the nested PCR products representing the 39 cleavage fragments that were cloned and sequenced for each gene. The frequency of 59 RACE clones corresponding to each cleavage site (arrows) is shown as a fraction, with the number of clones matching the target message in the numerator. As expected for mirna-directed cleavage (Llave et al., 2002), the clones ended at the nucleotide that pairs with the 10th nucleotide of mir160. Mismatched nucleotides are in lowercase letters. (B) mir160 and ARF10, ARF16, andarf17 mrnas accumulate differentially in Arabidopsis (Col-0) tissues. Shown is an RNA gel blot analysis of 30 mg total RNA prepared from root, rosette leaf, stem, cauline leaf, buds and inflorescence meristems (bud), flower, seedpods and embryos (silique), and 12-d-old seedling tissues with a DNA probe complementary to mir160. The blot was stripped and reprobed with an oligonucleotide complementary to U6 as a loading control. Also shown is RNA gel blot analysis of 7 mg of the same RNA with RNA probes complementary to 39 ARF10, 39 ARF16, or39 ARF17. The positions of full-length mrnas and 39 cleavage products are noted at the right. Ethidium bromide staining of 25S rrna is shown as a loading control. Ratios of mir160 to U6 RNA or ARF full-length mrna to 39 cleavage product are indicated below each lane. (C) RNA gel blot analysis of 10 mg total RNA from rosette leaves of wildtype Col-0, 35S:ARF17, and 35S:5mARF17 T1 plants with a 39 ARF17 probe. The absence (ÿ) or presence of a moderate (þ) or severe (þþ) rosette leaf phenotype and individual transformant numbers are indicated above the lanes. The positions of full-length ARF17 mrna and ARF17 39 cleavage product are noted at the right. The 25S rrna is shown as a loading control. The blot was rehybridized with an ACTIN2 probe, and normalized values of ARF17 full-length mrna and ARF17 39 cleavage product to ACTIN2 (At3g18780) RNA (with Col-0 levels set at 1.0) and ratios of ARF17 full-length mrna to 39 cleavage product are indicated.

4 mir160 Regulates Auxin Response Factors 1363 (35S:ARF17-1 and 35S:ARF17-4) (Figure 1C). The level of mir160-directed ARF17 39 cleavage product also was increased approximately threefold in these four 35S:ARF17 plants compared with wild-type Columbia-0 (Col-0) (Figure 1C). The modest increase in ARF17 mrna expression in 35S:ARF17 Arabidopsis plants prompted us to test the efficacy of our 35S:ARF17 construct with an Agrobacterium tumefaciens mediated transient expression assay in Nicotiana benthamiana leaves. Although we did not detect endogenous ARF17 mrna in N. benthamiana leaves, high levels of ARF17 full-length mrna accumulated after Agrobacterium inoculation with the 35S:ARF17 construct (data not shown), indicating that the 35S:ARF17 construct was capable of driving high level expression. The relatively low level of ARF17 mrna and the increased level of ARF17 39 cleavage product in 35S:ARF17 Arabidopsis plants (Figure 1C), along with the wide mir160 expression profile in wild-type plants (Figure 1B), suggested that endogenous mir160 may be sufficiently abundant to direct cleavage of much of the excess ARF17 mrna that results from 35S driven expression, at least in the relevant tissues. Plants Expressing mir160-resistant ARF17 Have Dramatic Developmental Defects The above results suggested that posttranscriptional regulation limits ARF17 mrna accumulation. To examine the importance of mir160 in this regulation, we constructed a mir160-resistant version of ARF17. This construct, designated 5mARF17, had five silent mutations within the mir160-complementary domain of an ARF17 genomic clone, thereby increasing the number of mismatches between mir160 and the ARF17 mrna from one in wild type to six without altering the amino acid sequence of the encoded ARF17 protein (Figure 2A). These substitutions also created an ApaLI restriction site (Figure 2A) that allowed us to distinguish RT-PCR products of 5mARF17 transcripts from those of ARF17. To determine whether mir160-directed cleavage of 5mARF17 RNA was reduced, we used an in vitro assay that relies on mirnaprogrammed silencing complexes endogenously present in wheat germ extract (Tang et al., 2003), which can direct cleavage of wild-type ARF17 RNA (G. Tang, M. Jones-Rhoades, D.P. Bartel, and P. Zamore, unpublished data), consistent with the presence of a mir160 locus in wheat (Triticum aestivum) (Jones- Rhoades and Bartel, 2004) and the cloning of mir160 from wheat germ extract (R. Rajagopalan and D.P. Bartel, unpublished data). For wild-type ARF17 RNA, but not for 5mARF17 RNA, a product of the size expected for mir160-directed cleavage was detected, demonstrating that 5mARF17 mismatches interfered with mir160-directed ARF17 cleavage (Figure 2B). To assess the in vivo consequences of disrupting mir160 regulation, we transformed wild-type Arabidopsis plants with 5mARF17 and control ARF17 genomic constructs under the control of the native ARF17 59 and 39 regulatory sequences. To preserve endogenous ARF17 transcriptional regulation, we included 1.9 kb of 59 flanking sequence, which extends ;120 bp into the annotated 39 UTR of the upstream gene At1g77840 (;60 bp downstream of the At1g77840 stop codon), and 1.7 kb of 39 flanking sequence, which stops 170 bp upstream of the stop Figure 2. 5mARF17 RNA Is Resistant to mir160-directed Cleavage. (A) Illustration of the ARF17 genomic clone used to transform Arabidopsis. Locations of the intergenic sequences (solid lines), introns (dashed lines), 59 and 39 UTRs (black boxes), DBD, remainder of the ARF17 open reading frame (open boxes), start (ATG) and stop (TAA) codons, mir160 complementary site (asterisk), and the EcoRI and SpeI sites used for cloning are shown. mrna segments from ARF17 and 5mARF17 are shown paired with mir160, with the ApaLI restriction site unique to 5mARF17 and the protein sequence encoded by the mir160 complementary sites of both mrnas indicated. (B) ARF17 RNA cleavage directed by mir160 in wheat germ extracts. The 59-radiolabeled transcripts prepared from wild-type ARF17 and 5mARF17 constructs were introduced into wheat germ extracts, and the time course of cleavage was examined on a denaturing polyacrylamide gel. The positions of ARF17 RNA and 59 ARF17 cleavage product are noted at the left. codon of the downstream, reverse-oriented gene At1g77855 (Figure 2A). None of the 101 control ARF17 primary transformants displayed obvious or consistent developmental anomalies (Figures 3 and 4). By contrast, 51 of 90 5mARF17 transformants displayed prominent vegetative and floral defects, including rosette and cauline leaf margin serration (Figures 3E to 3I), upward curling of the leaf margins (Figures 3E to 3I), reduced plant size (Figures 3E to 3I), accelerated flowering time (Figures 3G to 3I), altered phyllotaxy along the primary and lateral stems (Figures 4I and 4J), reduced petal size (Figures 4A to 4E), abnormal stamen structure (Figures 4F to 4H), and reduced fertility (Figures 4I and 4J). Eight of the 51 5mARF17 transformants with developmental abnormalities died before the transition to flowering, and five plants flowered successfully but were sterile.

5 1364 The Plant Cell 5mARF17 progeny plants displayed presumed embryonic defects, including symmetry anomalies in which the normally bilateral symmetry of germinating seedlings was instead trilateral or quadrilateral, resulting in one or two extra cotyledons (Figures 3A to 3C). Two individual hemizygous 5mARF17 plants segregated 10% (9/89) and 14% (11/81) progeny with three cotyledons and 3% (3/89) and 4% (3/81) progeny with four cotyledons. Plants with three cotyledons had three simultaneously emerging first true leaves (Figure 3B) at a frequency of 6/9 and 7/11, whereas plants with four cotyledons always had four simultaneously emerging first true leaves (Figure 3C). In addition to cotyledon-number defects, 5mARF17 cotyledons were often moderately or deeply lobed (heart shaped) and downwardly curled (Figures 3A to 3C). Abaxial cotyledon epidermal cells were altered in shape, size, and distribution (Figure 3D), perhaps contributing to the downward curling of the cotyledons. By contrast, adaxial epidermal cells of 5mARF17 cotyledons appeared relatively normal (data not shown). 5mARF17 plants displaying mild embryonic phenotypes (moderately lobed cotyledons) had mild vegetative and floral organ defects, whereas plants displaying more severe embryonic phenotypes (heartshaped or extra cotyledons) had more severe vegetative and floral organ defects. 5mARF17 seedlings also displayed reduced root and hypocotyl lengths and decreased root branching (Figure 4M). Like the primary transformants, 5mARF17 progeny plants flowered slightly earlier than control plants. Five independent Figure 3. Embryonic and Vegetative Developmental Defects of 5mARF17 Plants. (A) Eight-day-old control ARF17 and 5mARF17 T2 seedlings displaying a range of embryonic phenotypes, including lobed cotyledons (bottom left of 5mARF17 panel) and extra cotyledons indicative of trilateral and quadrilateral embryonic symmetry (right 5mARF17 panels). (B) Scanning electron micrographs of a 6-d-old control ARF17 (left) plant with the normal two cotyledons (asterisks) and two emerging true leaves (arrowheads) and 10-d-old 5mARF17 seedlings (middle and right). The middle plant displays a lobed cotyledon and a downwardly curled cotyledon, and the right plant has three cotyledons (asterisks) and three emerging true leaves (trilateral symmetry). (C) Scanning electron micrographs of a 10-d-old 5mARF17 seedling (left), 6-d-old 5mARF17 seedling (middle), and magnified view of the 6-dold 5mARF17 shoot apex (right). Indicative of quadrilaterally symmetric embryos, the 5mARF17 seedlings had four cotyledons (asterisks) and four emerging rosette leaves (arrowheads) rather than two of each. (D) Scanning electron micrographs of abaxial surfaces of 6-d-old control ARF17 and 5mARF17 cotyledons showing altered epidermal cell shape and organization in the 5mARF17 cotyledon. (E) Twenty-two-day-old control ARF17 (left) and 5mARF17 (top right) T2 plants and a 30-d-old 5mARF17 T2 plant with severe leaf serration (bottom right). (F) Rosette leaves from 30-d-old control ARF17 (left) and 5mARF17 (right) T2 plants. (G) Thirty-two-day-old control ARF17 (left) and 5mARF17 (right) T1 plants. The 5mARF17 plant has bolted. (H) Magnified view of 5mARF17 plant in (G). (I) Seventy-two-day-old control ARF17 and three 5mARF17 T2 plants grown in short days (8 h light, 16 h dark). The 5mARF17 developmental defects are increasing in severity from left to right; only the plant with the most severe phenotype has bolted. (J) Col-0 and 35S:5mARF17 rosette leaves with defects increasing in severity from left to right.

6 mir160 Regulates Auxin Response Factors mARF17 lines with similar rosette leaf defects flowered on day (n ¼ 34), (n ¼ 35), (n ¼ 30), (n ¼ 29), and (n ¼ 30), whereas nontransformed plants flowered on day (n ¼ 44), and four independent ARF17 control lines flowered on day (n ¼ 24), (n ¼ 24), (n ¼ 20), and (n ¼ 20). Despite extensive efforts, we were unable to recover a homozygous 5mARF17 plant with rosette leaf defects that was fertile. To observe the effects of disrupting both transcriptional and mir160 regulation of ARF17, we transformed wild-type plants with the 5mARF17 construct under the control of the CaMV 35S promoter (35S:5mARF17). Like the 5mARF17 primary transformants, nine out of 19 35S:5mARF17 primary transformants displayed embryonic (data not shown), rosette leaf (Figure 3J), and floral defects (Figures 4K and 4L). Although many of the phenotypes were qualitatively similar to those of 5mARF17 plants, the severity of leaf curling (Figure 3J) and floral organ defects (Figures 4K and 4L) and the frequency of premature death was greater in 35S:5mARF17 plants than 5mARF17 plants; two 35S:5mARF17 plants died as seedlings, four plants never transitioned to flowering but died after producing 8 to 10 upwardly curled and serrated rosette leaves, and three plants produced flowers with stamen and petal defects accompanied by reduced fertility or sterility. 5mARF17 and 35S:5mARF17 Plants Overaccumulate ARF17 mrna RNA gel blot analysis revealed that 12- and 16-d-old 5mARF17 seedlings displaying cotyledon phenotypes (Figures 5A and 5C), Figure 4. Floral Organ and Seedling Defects of 5mARF17 Plants. (A) Side view of the inflorescence bud cluster of ARF17 and 5mARF17 T1 plants. (B) Mature ARF17 and 5mARF17 T2 flowers with 5mARF17 petal and stamen defects increasing in severity from left to right. (C) to (E) Scanning electron micrographs of an ARF17 T2 mature flower (C) and 5mARF17 T2 mature flowers with moderate floral organ defects and reduced fertility ([D] and [E]). Petal (P) length and width was reduced in 5mARF17 flowers. (F) to (H) Scanning electron micrographs of ARF17 T2 (F) and 5mARF17 T2 (G) reproductive structures. Arrowheads point to the stamen. (H) is a magnified view of 5mARF17 stamen. Note the reduced filament length, abnormal stamen orientation, and reduced anther size in 5mARF17 flowers. The 5mARF17 flowers in (G) and (H) were from infertile plants. (I) Wild-type Col-0 plant floral phyllotaxy along the primary stem. (J) 5mARF17 floral phyllotaxy along the primary and lateral stems. 5mARF17 plants have reduced fertility, illustrated by the failure to produce filled seedpods. The inset shows a magnified view of altered 5mARF17 floral phyllotaxy. (K) Mature wild-type Col-0 flower (left) and Col-0 flower with sepals and petals removed to expose reproductive organs (right). (L) Mature 35S:5mARF17 flower (left) and 35S:5mARF17 flower with sepals and petals removed to expose reproductive organs (right). (M) Average hypocotyl length, primary root length, and lateral root number (expressed as a fraction of the primary root length) of 7-d-old untransformed control seedlings (C; black bars) and 7-d-old progeny of two ARF17 T3 homozygous lines (gray bars) and two 5mARF17 T3 heterozygous lines (white bars). The 5mARF17 progeny with wild-type phenotypes (ÿ) lacked the transgene; those with mild (þ) or severe (þþ) phenotypes scored positive for the 5mARF17 transgene. The number of plants analyzed is indicated above each bar, and error bars represent standard deviations.

7 1366 The Plant Cell as well as roots, rosette, and cauline leaves, stems, inflorescence meristems and buds, flowers, and siliques of older 5mARF17 plants with aberrant phenotypes (Figure 5C), all accumulated more ARF17 mrna than control plants. Moreover, 35S:5mARF17 transformants with aberrant phenotypes accumulated more full-length ARF17 mrna than 35S:ARF17 or untransformed wild-type Col-0 plants (Figure 1C). The increased ARF17 mrna levels in 5mARF17 and 35S:5mARF17 plants correlated with the severity of the rosette leaf phenotype (Figures 1C and 5B). The increase was not accompanied by a similar increase in 39 cleavage product accumulation (Figures 1 and 5), suggesting that the excess ARF17 mrna did not derive from the endogenous ARF17 gene but instead derived from the 5mARF17 or 35S:5mARF17 transgene, which both produce mrna refractory to mir160-directed cleavage (Figure 2B). Indeed, the fraction of ARF17 RT-PCR product that corresponded to 5mARF17, as monitored by ApaLI digestion, correlated with the presence of the rosette leaf phenotype (Figure 5D). mir160 levels were unchanged in control ARF17 and 5mARF17 plants (data not shown), indicating that an additional copy of ARF17 or 5mARF17 does not affect mir160 accumulation. Because developmental defects consistent with those observed in 5mARF17 and 35S:5mARF17 plants were not observed in 101 control ARF17 primary transformants or numerous progeny of these plants or 23 35S:ARF17 primary transformants, we conclude that the developmental phenotypes and increased Figure 5. Phenotypic 5mARF17 Plants Accumulate Increased ARF17 mrna. (A) RNA gel blot analysis of 7 mg total RNA prepared from 16-d-old T2 seedlings with a 39 ARF17 probe. The positions of full-length ARF17 mrna and ARF17 39 cleavage product are noted on the right. The 25S rrna is shown as a loading control. The blot was rehybridized with an ACTIN2 probe, and normalized values of ARF17 full-length mrna to ACTIN2 mrna (with Col-0 levels set at 1.0) and ratios of ARF17 fulllength mrna to 39 cleavage product are indicated. (B) RNA gelblot analysis of10 mg total RNA prepared from rosette leaves of 32-d-old T2 plants with a 39 ARF17 probe. The absence (ÿ) or presence of a mild (þ) or severe (þþ) rosette leaf phenotype is indicated. The positions of full-length ARF17 mrna and ARF17 39 cleavage products are noted at the right. Normalization was as in (A). (C) RNA gel blot analysis of 5 mg total RNA prepared from Col-0 and 5mARF17 root (r), rosette leaf (rl), stem (st), cauline leaf (cl), buds and inflorescence meristems (b), flower (f), silique and embryonic tissues (si), and 12-d-old seedling (se) tissues with a 39 ARF17 probe. The positions of full-length ARF17 mrna and ARF17 39 cleavage product are noted at the right. Two uncharacterized ARF17-hybridizing RNAs, one longer and one shorter than the full-length ARF17 transcript, accumulated in siliques of Col-0 and 5mARF17 plants. Ratios of ARF17 full-length mrna to 39 cleavage product are indicated, along with the relative levels of ARF17 mrna in 5mARF17 plants compared with wild-type Col-0 plants in each tissue (5mARF17/Col-0 ARF17). (D) Expression of 5mARF17 mrna in untransformed (untxf) and T2 transgenic plants that were individual progeny of the indicated primary transformant. Rosette leaf RNA was reverse transcribed, and the resulting cdna was PCR-amplified to completion, then digested with ApaLI, which cuts the 5mARF17 amplicon but not the ARF17 amplicon or heteroduplex molecules resulting from annealing ARF17 and 5mARF17 strands. Agarose gel separation and ethidium bromide staining revealed the full-length PCR product (330 bp) and an ApaLI digestion fragment (250 bp). Each lane is an analysis of an individual plant, and the presence (þ) or absence (ÿ) of a rosette leaf phenotype in that plant is indicated. DNA gel blot analysis (data not shown) was used to quantitate the 330- and 250-bp DNAs, and the amount of 5mARF17 relative to ARF17 after correcting for heteroduplex formation is shown below each lane.

8 mir160 Regulates Auxin Response Factors 1367 ARF17 mrna accumulation in 5mARF17 and 35S:5mARF17 plants resulted from disrupting mir160-directed ARF17 regulation, rather than from expressing an extra copy of ARF17. Together, these results show that mir160 is crucial for the posttranscriptional regulation of ARF17 expression and that this regulation is necessary for the proper growth and development of many Arabidopsis organs. 5mARF17 Plants Respond to Auxin Treatment 5mARF17 plants have fewer lateral roots and shorter hypocotyls than control plants (Figure 4M), traits characteristic of auxinresistant mutants (Estelle and Somerville, 1987; Hobbie and Estelle, 1995; Monroe-Augustus et al., 2003). 5mARF17 plants also have shorter primary roots than control plants (Figure 4M), whereas auxin-resistant mutants typically have increased root length (Estelle and Somerville, 1987; Hobbie and Estelle, 1995; Monroe-Augustus et al., 2003). To determine if the response to exogenous auxin was altered in 5mARF17 plants, we measured root and hypocotyl length after exposing plants to the auxin IAA. Even though the primary root and hypocotyl lengths of 5mARF17 plants were shorter than those of control plants before IAA treatment, IAA treatment decreased primary root and hypocotyl growth in 5mARF17 plants (Figures 6A and 6B), indicating that 5mARF17 plants were not dramatically impaired in these responses to exogenous IAA. IAA Treatment Does Not Appreciably Alter mir160, mir164, and mir167 Accumulation in Seedlings In addition to mir160 and mir167, which regulate ARF genes (Figure 1; Rhoades et al., 2002; Kasschau et al., 2003), mir164 regulates NAC-domain genes (Laufs et al., 2004; Mallory et al., 2004a), some of which are also implicated in auxin signaling (Xie et al., 2000; Aida et al., 2002; Furutani et al., 2004). To determine whether levels of mirnas implicated in auxin signaling are affected by exogenous IAA, we analyzed mir160, mir164, and mir167 accumulation in wild-type seedlings after IAA treatment. Although IAA treatment dramatically increased IAA1 and IAA19 early auxin-response transcript levels, mir160, mir164, and mir167 accumulation was not substantially altered at the times analyzed (Figure 6C), indicating that seedling levels of these mirnas are not greatly influenced by IAA application. Similarly, ARF17 mrna and ARF17 39 cleavage product levels appeared unaffected by IAA treatment (Figure 6C). These results do not provide evidence for a role of auxin in mediating ARF17, mir160, mir164, or mir167 expression in seedlings, suggesting that regulation by these three mirnas may instead be needed to set components of the auxin-response machinery to proper levels so that tissues can respond appropriately to auxin. all have Pro-Ser-Thr rich middle regions (see Supplemental Figure 1 online), suggesting that they may repress transcription of early auxin-response genes. To identify molecular changes in 5mARF17 plants, we monitored mrna accumulation profiles of five GH3-like transcripts and three Aux/IAA transcripts in seedlings and rosette leaves. In 16-d-old seedlings, we found levels of GH3.3 increased approximately fourfold in 5mARF17 seedlings compared with wild-type or control ARF17 seedlings (Figure 7A), whereas levels of IAA1 Aux/IAA mrna and other GH3-like mrnas, YDK1/GH3.2, GH3.5, DFL1/GH3.6, and DFL2/GH3.10, were similar in wild-type, control ARF17, and 5mARF17 seedlings (Figure 7). In addition to an approximately fourfold increase in GH3.3 mrna levels, YDK1/GH3.2 mrna levels were increased in 30-d-old 5mARF17 rosette leaves (Figure 8A). By contrast, levels of GH3.5 and DFL1/GH3.6 mrnas were decreased to undetectable levels in rosette leaves of 30-d-old 5mARF17 plants (Figure 8B), whereas levels of DFL2/GH3.10, IAA1, IAA17, and IAA19 mrnas were similar in rosette leaves of 30-d-old wildtype, control ARF17, and 5mARF17 plants (Figure 8A; data not shown). These results indicate that mir160-directed ARF17 regulation is necessary for proper expression of a subset of GH3-like mrnas and establish a molecular connection between ARF17 and early auxin responses. YDK1/GH3.2, GH3.3, GH3.5, and DFL1/GH3.6 are group II GH3-like proteins that conjugate IAA to amino acids in vitro (Staswick et al., 2005) and therefore are predicted to decrease active IAA levels in the cell. To determine if responsiveness to endogenous auxin is changed in 5mARF17 plants, we expressed control ARF17 and 5mARF17 in the widely used DR5- GUS reporter line (Guilfoyle, 1999). The DR5-GUS line expresses b-glucuronidase (GUS) under the control of a minimal CaMV 35S promoter fused to multiple copies of the TGTCTC AuxRE, first defined in the auxin responsive promoter of the soybean (Glycine max) GH3 gene. Rosette leaves of control DR5-GUS plants transformed with control ARF17 accumulated similar levels of GUS mrna as those of the parental line (Figure 8C), and control ARF17 DR5-GUS and DR5-GUS cotyledons and rosette leaves had a similar GUS expression pattern (data not shown). By contrast, phenotypic 5mARF17 DR5-GUS plants with increased ARF17 mrna levels accumulated approximately twofold more GUS mrna in rosette leaves (Figure 8C). These increases correlated with more expanded GUS expression in cotyledons and rosette leaves (Figures 8D and 8E), suggesting that auxin responsiveness or perhaps auxin is more widely distributed in cotyledons and rosette leaves of 5mARF17 plants. No change in overall GUS expression was observed in roots (data not shown), suggesting that root auxin responsiveness is not greatly influenced by 5mARF17 expression. Together, these results suggest that mir160-directed regulation of ARF17 is necessary for proper GH3-like gene expression and perhaps auxin distribution. mir160-directed ARF17 Regulation Is Necessary for Proper Expression of Certain GH3-Like Early Auxin Response Genes Although no molecular connections between ARF10, ARF16, and ARF17 and auxin responses have been reported, these ARFs DISCUSSION Here, we describe in vivo roles of mir160-directed ARF17 posttranscriptional regulation. Plants expressing a mirnaresistant version of ARF17 (5mARF17) display increased ARF17

9 1368 The Plant Cell mrna accumulation and exhibit dramatic developmental defects. These phenotypes correlate with reduced accumulation of GH3.5 and DFL1/GH3.6, two closely related mrnas (Staswick et al., 2002), and increased accumulation of YDK1/GH3.2 and GH3.3, two other closely related mrnas (Staswick et al., 2002), as well as DR5-GUS. These results indicate that mir160- directed regulation of ARF17 is critical for proper development, establish a molecular link between ARF17 and the auxin response pathway, and add another posttranscriptional regulatory dimension to ARF-mediated regulation. The Developmental Abnormalities of 5mARF17 Plants Overlap with Those of Plants Impaired in mirna Functioning In Arabidopsis, mutations in DCL1, AGO1, HYL1, and HEN1 impair the mirna pathway and lead to developmental defects that overlap with those exhibited by 5mARF17 plants. In particular, hypomorphic ago1 rosette leaves are serrated and ago1, hyl1, and hen1 null mutants exhibit upwardly curled rosette leaves and a dwarfed stature. Indeed, mir160 accumulation is reduced and ARF17 mrna accumulation is increased in dcl1, ago1, hyl1, and hen1 mutants (Kasschau et al., 2003; Vaucheret et al., 2004; Vazquez et al., 2004a), consistent with the possibility that reduced mir160-directed ARF17 regulation contributes to the developmental abnormalities of these mutants. Viral proteins can interfere with the mirna pathway and affect development when expressed in plants (Mallory et al., 2002; Kasschau et al., 2003; Chapman et al., 2004; Dunoyer et al., 2004). For example, Arabidopsis plants expressing viral proteins P1/HC-Pro, p19, p15, and p21 display reduced mirna-directed Figure 6. IAA Treatment of Wild-Type and 5mARF17 Plants. (A) Primary root length of 8-d-old Col-0 and T2 seedlings mock treated (white bars) or treated with 50 nm IAA (black bars). 5mARF17 progeny with wild-type phenotypes (ÿ) lacked the transgene; those with mild (þ) or severe (þþ) phenotypes scored positive for the 5mARF17 transgene. The number of plants analyzed is indicated above each bar, and standard deviations are shown (except in cases where fewer than three plants were analyzed). (B) Hypocotyl length of 7-d-old Col-0 and T2 seedlings mock treated (white bars) or treated with 0.1 mm IAA (gray bars) and 10 mm IAA (black bars). 5mARF17 progeny with wild-type phenotypes (ÿ) lacked the transgene; those with mild (þ) or severe (þþ) phenotypes scored positive for the 5mARF17 transgene. The number of plants analyzed is indicated above each bar, and standard deviations are shown (except in cases where fewer than three plants were analyzed). (C) RNA gel blot analysis of 10 mg total RNA prepared from IAA treated (þ) or untreated (ÿ) Col-0 seedlings with DNA probes complementary to mir160, mir164, mir167, or U6 (loading control). Blots were stripped between each hybridization. The positions of 32 P-labeled RNA oligonucleotides are noted at the right. Also shown are RNA gel blot analyses of 10 mg of the same RNA with RNA probes complementary to IAA1, IAA19, or 39 ARF17. The positions of each full-length transcript and ARF17 39 cleavage product are noted at the left. The 25S rrna is shown as a loading control. Normalized values of mirna to U6 RNA and IAA1, IAA19, and ARF17 full-length mrnas to ACTIN2 mrna (with untreated Col-0 levels set at 1.0) and ratios of ARF17 full-length mrna to 39 cleavage product are indicated.

10 mir160 Regulates Auxin Response Factors 1369 regulation can largely explain many of the developmental defects of plants expressing viral-encoded silencing suppressors. ARF10, ARF16, and ARF17 Are Similar to Repressing ARFs Figure 7. GH3-Like Expression in 16-d-Old 5mARF17 Seedlings. (A) RNA gel blot analysis of 7 mg of total RNA from 16-d-old seedlings with RNA probes complementary to YDK1/GH3.2 or GH3.3. The 25S rrna is shown as a loading control. Blots were rehybridized with an ACTIN2 probe, and normalized values of GH3-like mrnas to ACTIN2 mrna are indicated, with mrna levels in Col-0 control plants set at 1.0. (B) Similar analysis with a probe complementary to DFL2/GH3.10. Normalized values of DFL2 mrna to 25S rrna are indicated, with the mrna level in Col-0 control plants set at 1.0. (C) Similar analysis with a probe complementary to GH3.5, DFL1/GH3.6, IAA1, and 39 ARF17. Normalizations were as in (A). mrna cleavage and developmental abnormalities (Kasschau et al., 2003; Chapman et al., 2004; Dunoyer et al., 2004), whereas plants expressing viral proteins 2b, p38, and p25 are not impaired in mirna-directed target mrna cleavage and lack dramatic developmental defects (Chapman et al., 2004; Dunoyer et al., 2004). The small, serrated rosette leaf phenotypes of P1/ HC-Pro, p19, p15, and p21 expressing plants (Kasschau et al., 2003; Chapman et al., 2004; Dunoyer et al., 2004) are strikingly similar to those of 5mARF17 plants (Figures 3E to 3I). In addition, P1/HC-Pro plants exhibit reduced stamen size (Kasschau et al., 2003; Chapman et al., 2004; Dunoyer et al., 2004), and p19 and p15 plants have smaller petals (Chapman et al., 2004; Dunoyer et al., 2004), also reminiscent of 5mARF17 plants (Figures 4A to 4H), suggesting that disrupted mir160-directed ARF17 ARF proteins can either activate or repress transcription, depending on the nature of the middle domain (Ulmasov et al., 1999b). ARF5, ARF6, ARF7, ARF8, and ARF19 are activating ARFs with Gln-Leu-Ser rich middle regions, whereas ARF1, ARF2, ARF3, ARF4, and ARF9 are repressing ARFs with Pro- Ser-Thr rich middle regions (Ulmasov et al., 1999b; Tiwari et al., 2003). ARF10, ARF16, and ARF17 and their four rice homologs have Pro-Ser-Thr rich middle regions (see Supplemental Figure 1 online), suggesting that they might be repressors. However, the five known repressing ARFs are more closely related to the activating ARFs than either group is to ARF10, ARF16, and ARF17 (Remington et al., 2004), so ARF10, ARF16, and ARF17 may define a specialized ARF class. arf10 and arf16 lossof-function mutants do not display obvious developmental anomalies (Okushima et al., 2005). In an attempt to identify a loss-of-function arf17 mutant, we searched the Salk Institute Genomic Analysis Laboratory collection (Alonso et al., 2003) for plants with disruptions in ARF17. No mutants were found with insertions in the open reading frame of ARF17, but one mutant, SALK_062511, had a T-DNA inserted ;210 bp upstream of the ARF17 start codon. Plants homozygous for this insertion did not display obvious developmental defects (data not shown); however, RT-PCR revealed that these plants still accumulated ARF17 mrna (data not shown), indicating that this mutant was not a null allele. As has been demonstrated for other members of the ARF family (Okushima et al., 2005), functional redundancy among ARF10, ARF16, and ARF17 may preclude informative analyses of single arf mutants in this class. ARF17 Regulates GH3-Like Expression There are 20 Arabidopsis GH3 homologs, which fall into three clades. Group II GH3 proteins, including YDK1/GH3.2, GH3.3, GH3.5, and DFL1/GH3.6, conjugate IAA to amino acids in vitro (Staswick et al., 2005). Regulating IAA conjugation is important for maintaining endogenous IAA levels (Ljung et al., 2002); these GH3 proteins likely play an important role in auxin responsiveness by reducing active auxin levels and thus negatively regulating auxin signaling. We observed decreased GH3.5 and DFL1/GH3.6 mrna levels in 5mARF17 rosette leaves (Figures 8B). Because the Pro-Ser-Thr rich middle region of ARF17 is consistent with transcriptional repression, an appealing model is that ARF17 acts directly to repress GH3.5 and DFL1/GH3.6 transcription, and the consequent GH3.5 and DFL1/GH3.6 reduction leads to increased IAA, which is known to stimulate the expression of YDK1/GH3.2, GH3.3, and DR5-GUS. However, because it is currently impossible to discern whether the changes in GH3- like transcript levels were direct or indirect consequences of disrupting mir160-directed ARF17 regulation, there are many other possibilities. For example, ARF17 might indirectly regulate GH3-like and DR5-GUS expression by repressing an

11 1370 The Plant Cell activator of GH3.5 and DFL1/GH3.6 or by repressing a repressor of YDK1/GH3.2, GH3.3, and DR5-GUS. Although the exact mechanism by which ARF17 regulates GH3-like gene expression is unclear, the observation that mirna repression of ARF17 is important for proper GH3-like mrna accumulation provides an entry point for the study of early auxin-responsive gene expression. Gain-of-Function ydk1-d Mutants Mimic 5mARF17 Hypocotyl and Root Phenotypes YDK1/GH3.2 mrna is increased in 5mARF17 rosette leaves (Figure 8A) but is not obviously changed in 5mARF17 16-d-old seedlings (Figure 7A). Dominant ydk1-d mutants, which overexpress YDK1, have hypocotyl and root phenotypes (Takase et al., 2004) similar to 5mARF17 plants (Figure 4M). Both 5mARF17 and ydk1-d plants display reduced primary root length, lateral root number, hypocotyl length, and stature. ydk1-d plants also exhibit reduced apical dominance, which we did not observe in 5mARF17 plants. Because the hypocotyl and root phenotypes of 5mARF17 and ydk1-d plants are similar, it is possible that increased YDK1 expression contributes to 5mARF17 phenotypes. Gain-of-function dfl1-d plants are resistant to exogenous auxin and exhibit reduced lateral root number and hypocotyl length and a dwarf stature. Antisense DFL1 plants have increased lateral roots (Nakazawa et al., 2001), whereas 5mARF17 plants, which show reduced DFL1/GH3.6 mrna accumulation in rosette leaves (Figure 8B) but not in seedlings (Figure 7C), have fewer lateral roots (Figure 4M). The expression of at least two other GH3-like genes, GH3.3 and GH3.5, also is altered in 5mARF17 plants (Figures 7 and 8), but mutants in these genes have not been described. Because the expression of at least four GH3-like genes is altered in 5mARF17 plants, the Figure 8. GH3-Like Gene Expression Changes in 5mARF17 Rosette Leaves. (A) RNA gel blot analysis of 7 mg of total RNA from 30-d-old rosette leaves with RNA probes complementary to GH3.3, YDK1/GH3.2, DFL2/ GH3.10, IAA1, or39 ARF17. The absence (ÿ) or presence of a severe (þþ) rosette leaf phenotype is indicated. The positions of each full-length transcript and ARF17 39 cleavage product are noted at the right. The 25S rrna is shown as a loading control. Normalized values of mrnas to ACTIN2 mrna are indicated, with mrna levels in Col-0 control plants set at 1.0. Because YDK1 transcripts were not detected (ND) in Col-0 or ARF17 control rosette leaves, the approximately fourfold increase in YDK1 signal observed represents a minimum. (B) Similar analysis with probes complementary to GH3.5 or DFL1/ GH3.6. These transcripts were not detected in 5mARF17 rosette leaves, so the fold changes observed represent a minimum. (C) RNA gel blot analysis of 7 mg of total RNA from 30-d-old rosette leaves with RNA probes complementary to GUS or 39 ARF17. The absence (ÿ)or presence of a mild (þ) or severe (þþ) rosette leaf phenotype is indicated. The positions of each full-length transcript and ARF17 39 cleavage product are noted at the right. Normalization is as in (A). (D) Fifteen-day-old ARF17 DR5-GUS and 5mARF17 DR5-GUS cotyledons. (E) Twenty-eight-day-old ARF17 DR5-GUS and 5mARF17 DR5-GUS rosette leaves with moderate (top) or more severe (bottom) leaf defects.

12 mir160 Regulates Auxin Response Factors 1371 Table 1. Primer Sequences Primer Primer Sequence (59 39) 39 ARF10 for CAACTCTTCGGATCACCATCTCCGTC 39 ARF10 rev GGCCGTAATACGACTCACTATAGGCCGAGAGATCGAGTGTGCGTCC 39 ARF16 for GGCCGTAATACGACTCACTATAGGGGCTCCTGATGCATCCCTATAGAG 39 ARF16 rev GACAATGGGAACAACACCATGC 39 ARF17 for GAGATGATGAACTTTGGCAGTCCG 39 ARF17 rev GGCCGTAATACGACTCACTATAGGCCTGCGTTGCTGTTGGAATGTT GH3.5 for CGTGAAGTTTGCTCCAATTATCGAGC GH3.5 rev GGCCGTAATACGACTCACTATAGGCATGATAAAAGTCCAAACAATTCCGC GH3.17 for CAAGCAGAGCATGTACTGTCAGCTTC GH3.17 rev GGCCGTAATACGACTCACTATAGGGCCATCGAACCAGTCACAATCACC GH3.6 for GGCCGTAATACGACTCACTATAGGCCTTGGTGTCTTGTACTGATTGATCG GH3.6 rev GCGCCTCAGTTCAGCTTCATATGC GH3.10 for GGCCGTAATACGACTCACTATAGGGCATGTGAAACAATAGGCACCAAGG GH3.10 rev TTTGTCTGAAGACAAAGGGAGTGGTC GH3.2 for CCATCCTCTGCTCCGACTCGTCC GH3.2 rev GGCCGTAATACGACTCACTATAGGGCTCCAGTAACAATCACGTCGAGG GH3.3 for GGCCGTAATACGACTCACTATAGGGACAAGATCATTGACCGGTCACC GH3.3 rev GCTCTGCGATCTCCGATGATGC IAA1 for GGAAGTCACCAATGGGCTTAACC IAA1 rev GGCCGTAATACGACTCACTATAGGGCAGGAGGAGGAGCAGATTCTTCTG IAA17 for GGCAGTGTCGAGCTGAATCTGAGG IAA17 rev GGCCGTAATACGACTCACTATAGGGGCCGGAGGTTTGGCTGGATC IAA19 for GAGAAGGAAGGACTCGGGCTTGAG IAA19 rev GGCCGTAATACGACTCACTATAGGCGATGCCACGGAAACCGAAGAG GUS for GCGTATCGTGCTGCGTTTCGATG GUS rev GGCCGTAATACGACTCACTATAGGGCCTTGTCCAGTTGCAACCACC ARF17 mut for CGTATTCTACATTTCCTGCTGGAATGCAAGGTGCACGGCAATATGATTTTGGGTCTTTCAATCC ARF17 mut rev GGATTGAAAGACCCAAAATCATATTGCCGTGCACCTTGCATTCCAGCAGGAAATGTAGAATACG ARF17 wg for GGCCGTAATACGACTCACTATAGGCCTCTACTTATCAGGAGACCG ARF17 wg rev CCAGAGGACAGATTAGTGGTG ARF10.1 RACE CCGAGAGATCGAGTGTGCGTCC ARF10.2 RACE CTCCTCCGCTTCCGCCTCTTCTTCC ARF16.1 RACE GGCTCCTGATGCATCCCTATAGAG ARF16.2 RACE CCAAACCTGATGCATCATGAACTTTCTTAG ARF17.1 RACE CTTGGGAGCTAGAACCTGCGTTGC ARF17.2 RACE CCTGCGTTGCTGTTGGAATGTTTAGG ARF17-RT for GGATGAACCTGAGATTCTGC ARF17-RT rev TGGTGAATAGCTGGGGAGGATTTC ARF17 screen CTCTCTGTGGATCACTCTGTCTTGC pgii far SacI GCTATGACCATGATTACGCC ACTIN2 AACCCTCGTAGATTGGCACA for, forward; rev, reverse; mut, mutagenesis; wg, wheat germ. contribution of individual GH3-like genes to the 5mARF17 developmental defects is difficult to determine. The expression of GH3.5, DFL1, YDK1, and GH3.3 is altered in 5mARF17 rosette leaves (Figures 8A and 8B), but among the five GH3-like transcripts we monitored, GH3.3 was the only transcript noticeably changed in both 5mARF17 seedlings and rosette leaves (Figures 7A and 8A). Because the increase in ARF17 expression is consistently greater in 5mARF17 rosette leaves than in 5mARF17 seedlings (Figures 5, 7, and 8), it is possible that ARF17 levels in seedlings are not sufficient to detectably change expression of the other GH3-like genes when RNA from entire seedlings is pooled for RNA gel blot analysis. Indeed, we observe expanded DR5-GUS expression domains not only in rosette leaves, but also in 5mARF17 cotyledons (Figures 8D and 8E), and the changes visualized by histochemical staining are obvious, whereas the changes in GUS mrna levels are modest (approximately twofold). mir160 and mir167 May Coordinately Modulate GH3-Like Expression ARF8 appears to negatively regulate free IAA levels by controlling GH3-like gene expression (Tian et al., 2004). Levels of three GH3 mrnas, GH3.5, DFL1/GH3.6, and GH3.17, are reduced in arf8 loss-of-function mutants and increased in ARF8 overexpressing plants (Tian et al., 2004), suggesting that ARF8 activates GH3-like