| HOME | HELP | FEEDBACK | SUBSCRIPTIONS | ARCHIVE | SEARCH | TABLE OF CONTENTS |
Previous Article | Next Article ![]()
Eukaryotic Cell, January 2005, p. 121-133, Vol. 4, No. 1
1535-9778/05/$08.00+0 doi:10.1128/EC.4.1.121-133.2005
Copyright © 2005, American Society for Microbiology. All Rights Reserved.
Department for Organismic Interactions, Max Planck Institute for Terrestrial Microbiology, Marburg, Germany
Received 26 August 2004/ Accepted 26 October 2004
| ABSTRACT |
|---|
|
|
|---|
strains were affected in filamentous growth and pathogenicity. Rrm4 is an RRM protein with a so far unique domain organization consisting of three N-terminal RRMs as well as a domain found in the C terminus of poly(A)-binding proteins. These results indicate a role for RNA-binding proteins in regulation of morphology as well as in pathogenic development in U. maydis. | INTRODUCTION |
|---|
|
|
|---|
Mating compatibility is determined by a tetrapolar mating system consisting of two mating-type loci designated a and b. The biallelic a locus (a1 and a2) encodes the precursor and receptor of lipopeptide pheromones regulating intercellular recognition (9). The multiallelic b locus encodes a pair of homeodomain proteins, bE and bW, that are only active as heterodimeric transcription factors with subunits derived from different allelic origins (29). The formation of an active bE/bW heterodimer is necessary and sufficient for the developmental switch to fast polar growth. The concerted action of mating-type gene products is orchestrated by cross talk between protein kinase A (PKA) and mitogen-activated protein kinase (MAPK) signaling (17, 36, 54).
To investigate regulation of cell polarity and morphology during pathogenic development, we tested the involvement of RNA-binding proteins, since this class of proteins determines cell polarity in higher eukaryotes during developmental programs such as oogenesis and early embryogenesis in Drosophila melanogaster and sex determination in Caenorhabditis elegans, as well as mammalian stem cell proliferation (34). RNA-binding proteins regulate mRNA activity at the level of maturation, nuclear export, cytoplasmic localization, stability, and translation (15, 40, 42). A well-studied example is Pumilio from D. melanogaster, which recognizes a specific response element in the 3'-untranslated region of hunchback mRNA. Thereby translation is repressed, resulting in a protein gradient throughout the embryo that determines the posterior polarity of the embryo (44, 69). Also regarding lower eukaryotes, the number of reports that describe connections between RNA-binding proteins and development is growing. PufA from Dictyostelium discoideum is involved in determining sporulation in response to starvation by regulating translation of mRNA encoding the catalytic subunit of PKA (60). Nrd1 and Msa1 are negative regulators of sexual differentiation in Schizosaccharomyces pombe (26), and Scp160p regulates mating in Saccharomyces cerevisiae by interacting with the heterotrimeric G protein Gpa1p during pheromone-responsive MAPK signaling (22).
To study the function of RNA-binding proteins in U. maydis, we chose a combinatorial approach based on bioinformatics and reverse genetics, using the recently released genome sequence (http://www.broad.mit.edu/annotation/fungi/ustilago_maydis/) and an efficient gene replacement strategy (10, 28). Genes encoding potential RNA-binding proteins can be easily identified due to the presence of at least one RNA-binding domain. In eukaryotes, the majority of RNA-binding domains belong to the following four classes that are well defined on primary, secondary, and tertiary structure: the Pumilio homology domain (PUM), the K homology domain (KHD), the double-stranded RNA-binding domain (DSRM), and the RNA recognition motif (RRM) (12, 41, 49).
The PUM was named after its founding member, Pumilio (69, 73). The domain recognizes RNA sequence specifically and is composed of eight tandem repeats, each
40 amino acids (aa) in length. Characteristic of these repeats is a core consensus containing aromatic and basic amino acids (70). The KHD was first identified in hnRNP K (heterogeneous nuclear ribonucleoprotein K), a major pre-mRNA-binding protein (58). The domain comprises
60 aa and contains a characteristic Gly-X-X-Gly loop that grips RNA (38, 45). The double-stranded RNA-binding motif is an
65-aa domain found in various proteins that bind double-stranded RNAs such as Staufen from D. melanogaster or RNase III from Escherichia coli (31, 61). The RRM is by far the best-characterized and most widespread domain. At present there are 4,160 proteins with 7,079 domains known (http://smart.embl.de/). The RRM is 90 to 100 aa in length, containing two well-conserved sequences, the RNP1 octamer and the RNP2 hexamer (6, 12). Charged and aromatic side chains of RNP1 and RNP2 are solvent exposed and contact RNA (12).
Here, we describe the identification of 27 open reading frames (ORFs) in U. maydis encoding putative RNA-binding proteins of the PUM, KHD, DSRM, and RRM types. A subset of 18 gene replacement mutants were generated and subjected to a detailed phenotypic analysis revealing that KHD and RRM candidates are important for pathogenic development.
| MATERIALS AND METHODS |
|---|
|
|
|---|
(Bethesda Research Laboratories) and Top10 (Invitrogen) were used for cloning purposes. U. maydis strains FB1, FB2, UM521, SG200, FBD12-17, and FBD11-7 have been published previously (5, 8). Bioinformatic analysis. In order to identify PUM-type RNA-binding proteins, the following protein sequences (http://smart.embl-heidelberg.de/smart/do_annotation.pl?DOMAIN=Pumilio&BLAST=DUMMY&EVOLUTION=Show#Evolution) have been compared to the genomic sequence from U. maydis (http://www.broad.mit.edu/annotation/fungi/ustilago_maydis/) by using BLAST (2): S. cerevisiae (7), Q12221 [GenBank] , P47077 [GenBank] , P39016 [GenBank] , Q07807 [GenBank] , P25339 [GenBank] , Q04373 [GenBank] , and P47135 [GenBank] ; S. pombe (9), Q10238 [GenBank] , Q09829 [GenBank] , Q92347 [GenBank] , O94462, O60059, Q92359 [GenBank] , Q9P789, Q9UU76, and O74438 [GenBank] ; Neurospora crassa (2), Q8X035, Q9P6D4; D. discoideum (1), Q9Y1J5; Arabidopsis thaliana (6), Q9LRZ3, Q9SVZ7, BAB08273 [GenBank] Q95547, Q9SFL0, and Q9LVC3; H. sapiens (3), AAG31807 [GenBank] O00234 [SwissProt] , and AAG31806 [GenBank] C. elegans (6), O01322, QZ3191, Q9U2G4, Q09487, Q9N3Q2, and Q9N5M6; and D. melanogaster (1), Q9VHH6.
To identify KHD-containing proteins, we used the following protein sequences from the SMART (Simple Modular Architecture Research Tool) database (http://smart.embl-heidelberg.de/smart/do_annotation.pl?DOMAIN=KH&BLAST=DUMMY&EVOLUTION=Show#Evolution): S. cerevisiae (6), Q07834 [GenBank] , Q12186 [GenBank] , YBD2, PBP2, RS3, and Q99216 [GenBank] ; S. pombe (7), O59810, O74359 [GenBank] , O74919, O74555 [GenBank] , O74777 [GenBank] , YAJE_Schpo, and O14044 [GenBank] ; N. crassa (1), Q9P6B5; H. sapiens (5), Q13601 [GenBank] , AAG02184 [GenBank] AAG45476 [GenBank] Q9Y307, and Vigilin; C. elegans (1), Q9U3B2; and D. melanogaster (1), Q9V8H6.
In order to identify DSRM-type proteins, we used the following protein sequences (http://smart.embl-heidelberg.de/smart/do_annotation.pl?DOMAIN=DSRM&BLAST=DUMMY&EVOLUTION=Show#Evolution): S. cerevisiae (2), RM03 and RNT1; S. pombe (2), O43042 and PAC1; N. crassa (2), Q7S6N0 and Q9P6D6; and D. melanogaster (23), AAQ23615 [GenBank] Q8MRC7, Q8MRY4, Q95YG3, Q960Y4, Q961S7, Q9GNJ2, Q9NHW9, Q9NII2, Q9U6N4, Q9V847, Q9V9V7, Q9VCU9, Q9VHB0, Q9VJY9, Q9VLW8, Q9VRL8, Q9VY45, Q9VZ88, Q9W5S7, Q9XYN5, MLE, and STAU.
To identify RRM-type proteins, we used the following sequences (http://smart.embl-heidelberg.de/smart/do_annotation.pl?DOMAIN=RRM&BLAST=DUMMY&EVOLUTION=Show#Evolution): S. cerevisiae (6), P04147 [GenBank] , P32588 [GenBank] , P38922 [GenBank] , P32831 [GenBank] , P27476 [GenBank] , and P34217 [GenBank] ; S. pombe (2), Q09702 [GenBank] and P08965 [GenBank] ; H. sapiens (3), Q9UNP9, Q15717 [GenBank] , and Q9Y3B4; and D. melanogaster (5), P19018 [GenBank] , P19339 [GenBank] , P16914 [GenBank] , O02374, and AAF54070 [GenBank]
For the phylogenetic analysis, sequence alignments were performed with CLUSTAL_X (63). Dendrograms were constructed with MEGA 2.1 (http://www.megasoftware.net) (35) with the minimum evolution algorithm and an initial tree obtained by the neighbor-joining method. Two KHD proteins each have been predicted from C. neoformans (Cn; http://www.broad.mit.edu/annotation/fungi/cryptococcus_neoformans/) and Coprinus cinereus (Cc; http://www.broad.mit.edu/annotation/fungi/coprinus_cinereus/). KHD protein Cn1, 1,113 aa in length, is located on contig 1.77 between nucleotide positions 175691 and 172233. Two introns were predicted (intron 1, 173461 to 1733960; and intron 2, 172828 to 172779). KHD protein Cn2, 891 aa in length, is located on contig 1.12 between nucleotide positions 65877 and 61787. Eight introns were predicted (intron 1, 65225 to 65028; intron 2, 64929 to 64881; intron 3, 64745 to 64107; intron 4, 63950 to 63919; intron 5, 63720 to 63383; intron 6, 62923 to 62877; intron 7, 62675 to 62618; intron 8, 62371 to 62318). KHD protein Cc1, 1,228 aa in length, is located on contig 1.44 between nucleotide positions 260066 and 264281. Seven introns were predicted (intron 1, 260264 to 260745; intron 2, 261695 to 264754; intron 3, 262019 to 262074; intron 4, 262244 to 262303; intron 5, 262418 to 262528; intron 6, 262595 to 262654; intron 7, 263275 to 263334). KHD protein Cc2, 797 aa in length, is located on contig 1.9 between nucleotide positions 41512 and 44651. Eight introns were predicted (intron 1, 41834 to 41904; intron 2, 42041 to 42094; intron 3, 42237 to 42284; intron 4, 42550 to 42610; intron 5, 42821 to 43159; intron 6, 43619 to 43671; intron 7, 43873 to 43920; intron 8, 44306 to 44376).
Generation of gene replacement strains in U. maydis.
Gene replacement mutants were generated by a PCR-based approach (10, 28). For each gene replacement construct, eight PCR primers were designed: u1, u2, and u3 (upstream flanking sequence); d1, d2, and d3 (downstream flanking sequence); and p1 and p2 (10) (see Fig. 2A and Table 2; see Table S1 in the supplemental material). Primer u3 and d1 introduced SfiI(u) and SfiI(d) recognition sites at the 3' end of the upstream flank and the 5' end of the downstream flank, respectively. Initially
1-kb upstream and downstream flanking sequences of the ORF region of interest were amplified with genomic DNA of wild-type strain UM521 (a1b1) as a template and primer combinations u2/u3 and d1/d2, respectively (Fig. 2B, lanes 1 and 2; Table 2). PCR products were cleaved with SfiI and ligated in the presence of a compatible 1.8-kb SfiI(u)/SfiI(d) fragment containing the hygromycin resistance cassette (28). Gel-purified ligation products were inserted in pCR2.1-TOPO (Invitrogen; Table 2, pUMa). The resulting plasmids were used as a template for PCR with primers u2 and d2 (Fig. 2B, lane 3), and the respective products were transformed in U. maydis wild-type strains FB1 and FB2 (a1b1 and a2b2, respectively). Hygromycin-resistant colonies were subjected to whole-cell PCR analysis using primers p1 and p2 (Fig. 2B, lanes 4 to 6; Table 2). Amplification of an
200-bp-long fragment was diagnostic for the presence of the wild-type allele indicative for a nonhomologous integration event (Fig. 2B, lane 4). Only those transformants in which the diagnostic PCR product was absent were analyzed further (Fig. 2B, lanes 5 and 6). Successful gene replacements were verified by detection of an
1-kb PCR product using primers u1 and hu (Fig. 2B, lanes 7 to 9) and by amplification of
1-kb PCR products with primer combinations d3 and hd (Fig. 2B, lanes 10 to 12). To validate this approach, we verified successful gene replacements in the cases of pum1-3
and khd4
, as well as rrm4
, by Southern analysis (10).
|
|
Measuring doubling time during exponential growth. Starter cultures were grown overnight on a rotary shaker at 220 rpm in liquid CM at 28°C (24). A 1-µl cell suspension of stationary-phase starter cultures was transferred in 1 ml of CM and incubated on a thermomixer compact (Eppendorf) at 400 rpm in a 24-well plate at 28°C. At appropriate time intervals, optical density (OD) was measured as ODA600 with a TECAN Saphire fluorescence reader. Doubling time was calculated from the slope of a linear regression line derived from values during the exponential phase of growth (see Fig. S1 in the supplemental material).
Quantifying response to external cAMP.
Starter cultures (1 ml) were grown overnight on a thermomixer compact (Eppendorf) at 1,100 rpm in liquid CM (24) at 28°C using 2-ml tubes sealed with LidBac (Eppendorf). Two-microliter cell suspensions of stationary-phase starter cultures were transferred in 1 ml of CM in the presence or absence of 15 mM cyclic AMP (cAMP) and incubated on a thermomixer compact (Eppendorf) at 1,100 rpm for 18 h at 28°C. An aliquot was transferred to a Neubauer counting chamber, and the percentage of cells in aggregates of more than two cells was counted (for each quantification,
100 cells; see Fig. S2 in the supplemental material).
Mating, pheromone production, and filamentation assays. Starter cultures were grown overnight on a rotary shaker at 220 rpm in PD medium (2.4% potato dextrose broth; Difco). Stationary-phase cultures were diluted in PD medium to an OD at 600 nm (OD600) of 0.2. After growth at 28°C on a rotary shaker at 220 rpm to an OD600 of 0.5, a 0.5-ml aliquot of cell suspension was centrifuged at 3,500 rpm (5 min at room temperature). Cells were resuspended in 0.5 ml of water. For filamentation assays, cells were applied as spots and, for mating as well as pheromone production assays, respective strains were applied as cospots on charcoal-containing PD plates that were sealed with Parafilm and incubated at 22°C for 24 to 48 h.
Cell morphology and conjugation hypha formation. Starter cultures were grown overnight on a rotary shaker at 220 rpm in CM (24) at 28°C. Stationary-phase cultures were diluted in CM to an OD600 of 0.2. After incubation as mentioned above to an OD600 of 0.5, a 1-ml aliquot of cells was centrifuged at 3,500 rpm (5 min at room temperature) and cells were resuspended in CM with or without synthetic a1 or a2 pheromone dissolved in dimethyl sulfoxide (2.5 µg/ml) (62). Cells were incubated for 5 h at 28°C in a 15-ml plastic tube rotating with 20 rpm. For microscopic observation, we used a Zeiss Axiophot with differential interference contrast optics. Pictures were taken with a charge-coupled device camera (AxioCam HRm; Zeiss).
Pathogenicity assay. Plant infections of corn variety Early Golden Bantam (Olds Seeds, Madison, Wis.) were performed as described previously (11). Tumor formation was scored after 14 days. In the case of rrm4 deletion mutants, a more detailed disease rating was adopted from earlier studies (20) using the following categories: no tumor, minor disease symptoms at leaves (disease area with tumors is not larger than 1.5 cm and the number of tumors does not exceed a cutoff of 20), major disease symptoms at leaves (disease area with tumors is wider than 1.5 cm and the number of tumors exceeds 20), minor disease symptoms at stems (disease area with tumors is not larger than 4 cm and heights of tumors do not exceed 0.6 cm), major disease symptoms at stems (disease with tumors is larger than 4 cm and heights of tumors exceed 0.6 cm), as well as wilted or dead plants.
| RESULTS |
|---|
|
|
|---|
|
|
To identify KHD proteins, we chose 21 candidates from SMART: 14 from fungi, 5 from H. sapiens, and 1 each from D. melanogaster and C. elegans. Nine candidate ORFs were found and designated khd1 to -9 (Fig. 1; Table 1). Khd1 and Khd2 share the same domain organization consisting of three KHD regions and are similar in sequence to the Pab1p-binding protein Pbp2p that is required for proper chromosome segregation during meiosis (39, 51). Khd4 and Khd5 are KHD proteins that contain 5 and 11 KHD regions, respectively. Both are more than 1,000 aa in length. Khd5 is most likely the homologue of Scp160p from S. cerevisiae (18, 67), since it shares substantial sequence similarity over the entire length of the protein (BLAST expect value of 9 x e42; Table 1). Khd6 is related to the branch point bridging protein Msl5p involved in splicing (1). According to their similarity to proteins from S. cerevisiae, Khd7, Khd8, and Khd9 are likely to be involved in ribosome maturation and function. Khd7 exhibits significant sequence similarity to Krr1p, a protein involved in ribosomal small subunit assembly and maintenance (55). The Khd8-related protein Pno1p is a component of the 90S preribosome involved in ribosome maturation (57). Khd9 displays high sequence similarity to the ribosomal protein S3 (Rps3p), an integral part of the small subunit of ribosomes (Fig. 1; Table 1) (66).
In order to identify DSRM proteins in U. maydis, we chose 29 proteins containing the DSRM: 6 from fungi and 23 from D. melanogaster. Only one candidate ORF, termed dsr1, encodes a DSRM domain (Fig. 1; Table 1). Dsr1 exhibits high sequence similarity to the mitochondrial precursor of the 60S ribosomal protein L3 from S. cerevisiae (Mrpl3p, Table 1) (21). Both proteins share a characteristic RNase III RIBOc domain in N-terminal proximity to the DSRM (http://smart.embl.de/). The presence of only few DSRM proteins was not surprising since S. cerevisiae, S. pombe, and N. crassa also contain only two each (http://smart.embl.de/). Apparently, U. maydis lacks a clear homologue of Rnt1p, the second DSRM protein from S. cerevisiae involved in small nucleolar RNA metabolism (13).
To detect RRM proteins in U. maydis, we chose 16 RRM-containing proteins from SMART: 8 from fungi, 3 from H. sapiens, and 5 from D. melanogaster. Since the number of RRM proteins is substantially higher than the number of proteins containing the other domains (e.g., S. cerevisiae contains 53), we chose to list only 12 candidates (Fig. 1; Table 1). Rrm1 and Rrm5 contain two central RRM domains, while Rrm11 harbors two such motifs near the N terminus. Rrm2 exhibits a distinct domain organization consisting of an RRM domain adjacent to a SWAP (suppressor-of-white-apricot splicing regulator) domain and an RPR (regulation of pre-mRNA) domain. The latter two domains are found in proteins involved in regulation of processing nuclear pre-mRNA (http://smart.embl.de/). Proteins sharing this particular domain organization are also found in C. elegans, H. sapiens, D. melanogaster, and A. thaliana but are absent in S. pombe or S. cerevisiae (http://smart.embl.de/). Rrm12 is most likely the functional homologue of the poly(A)-binding Pab1p. It contains the characteristic domain organization of poly(A)-binding proteins from other species consisting of four N-terminal RRM domains in combination with a C-terminal PABC domain [poly(A)-binding protein C-terminal domain; Fig. 1; Table 1; BLAST expect value of e132]. Rrm4 shares a similar domain composition. However, this potential RNA-binding protein exhibits a novel domain architecture, since it contains only three N-terminal RRM domains that are combined with a C-terminal PABC domain (http://smart.embl.de/). Rrm8 contains four central RRM domains and is similar to Ylp184p from S. cerevisiae and to the negative regulator of development Nrd1 from S. pombe (64). Rrm9 carries in addition to its RRM region an R3H domain, which exhibits a characteristic spacing of conserved arginine (R) and histidine (H) residues and is implicated in binding single-stranded nucleic acids. This domain organization has not been found in proteins of the aforementioned model organisms (http://smart.embl.de). In summary, using known RNA-binding proteins as a template, we were able to identify 27 ORFs in U. maydis encoding 5 proteins of the PUM type, 9 of the KHD type, a single protein of the DSRM type, and 12 of the RRM type.
Phenotyping of mutants carrying deletions in genes encoding proteins of the PUM, KHD, and RRM types. In order to investigate the function of the identified ORFs, we chose a subset of 18 ORFs for a reverse genetic approach: pum1 to -3, khd1 to -4, and rrm1 to -11 (Tables 1 and 2). We focused on those genes that were not highly similar in sequence to essential genes from S. cerevisiae. Therefore, proteins such as Khd6 to Khd9, potentially involved in ribosomal processes, or the putative poly(A)-binding protein Rrm12 were omitted (Table 1). Respective gene replacement mutants were obtained by a PCR-based approach (Materials and Methods; Table 2 and Fig. 2) (10, 28). For each gene, two independent deletion strains were generated in strain FB1 as well as FB2 (a1b1 and a2b2, respectively; Table 2).
All deletion strains were viable and were subjected to a detailed phenotypic analysis assaying morphology and pigmentation of colonies, growth on solid minimal medium at various temperatures, doubling time in liquid medium, cell morphology, cAMP-dependent cell aggregation, mating on plates containing activated charcoal, as well as disease symptoms on corn seedlings (Materials and Methods) (Tables 3 and 4). None of the deletions caused auxotrophy, heat sensitivity, abnormal morphology, or pigmentation of colonies (Table 3). To score the cAMP response, deletion strains were grown in the presence of 15 mM cAMP (Materials and Methods). Wild-type strains respond to the presence of external cAMP with a characteristic cell separation defect resulting in cell aggregates (Materials and Methods) (20, 33). Only khd4
strains differed from the wild type by containing a large proportion of cell aggregates in the absence of external cAMP. This phenotype was not enhanced by cAMP addition (see Fig. S2 in the supplemental material).
|
|
, khd4
, and rrm4
strains (Tables 3 and 4). In comparison to the wild type, all khd1
strains were reduced in growth at 15°C on plates as well as in liquid culture, indicating a cold-sensitive growth phenotype (Table 3) (data not shown). khd4
and rrm4
strains were disturbed in mating and pathogenicity. In addition, khd4
strains differed from the wild type in growth rate and cell morphology (Table 3).
khd4
strains are affected in growth, pheromone response, filamentation, and cell morphology, as well as pathogenicity.
Analysis of khd4
strains revealed a number of aberrant phenotypes. First, mutants showed a twofold longer doubling time of 330 min with respect to wild-type strains (Materials and Methods) (see Fig. S1 in the supplemental material). Second, mating assays on activated charcoal-containing plates revealed an almost complete absence of dikaryotic hyphae in crosses of compatible khd4 deletion strains and a significant reduction when khd4
strains were cospotted with wild-type strains, indicative of defects in pheromone response and filamentation (Fig. 3A). Third, to test whether pheromone secretion was disturbed in khd4
strains, serial dilutions were cospotted with the diploid pheromone tester strain FBD12-17 (a2a2b1b2) on plates containing activated charcoal. In response to a1 pheromone, FBD12-17 switches to filament formation without prior cell fusion (5). khd4 deletion mutants were significantly reduced in elicitation of filament formation in FBD12-17, indicating reduced pheromone production (Fig. 3B). Fourth, to assay conjugation tube formation, wild-type as well as khd4
strains were incubated for 6 h in liquid CM with synthetic a1 or a2 pheromone (Materials and Methods) (62). Untreated khd4
cells were defective in cytokinesis, often shaped like lemons or droplets, and showed an increased diameter with respect to the wild type (Fig. 3C). Upon pheromone stimulation, only a few khd4
cells responded by formation of conjugation tubes, and these were substantially shorter and thicker than the wild type (Fig. 3C).
|
Sixth, to test for virulence of khd4
strains, 7-day-old corn seedlings were infected with a cross of compatible deletion strains (Materials and Methods). Formation of tumors was scored after 2 weeks and compared to that of the wild type. Mixtures of FB1khd4
#1 with FB2khd4
#1 and FB1khd4
#2 with FB2khd4
#2 showed 16 and 27% tumor formation, respectively. This was approximately fourfold lower than the tumor formation of a mixture of wild-type strains (80% tumor formation in FB1 x FB2 crosses; Table 4). In summary, loss of Khd4 results in reduced pheromone secretion, conjugation tube formation, and b-dependent filamentation, as well as pathogenicity.
rrm4
strains are reduced in filamentous growth and virulence.
According to our initial phenotypic analysis, rrm4
strains were slightly affected in mating (Table 3). Since the mating defects were only subtle, we performed more elaborate assays. Cospotting of strain FB2rrm4
#1 with serial dilutions of either FB1 (a1b1) or FB1rrm4
#1 on activated charcoal-containing plates confirmed that rrm4 deletion strains were impaired in mating (Fig. 4A, top), which could be due to defects in a-dependent pheromone response or in b-dependent filamentation. To test pheromone response, respective strains were cospotted in serial dilutions with diploid pheromone tester strain FBD12-17 (a2a2b1b2) (5). We observed no difference between elicitation of pheromone-induced filamentation of the diploid tester strain by either the wild type or FB1rrm4
#1 (Fig. 4B). We also tested conjugation tube formation by treatment of FB1rrm4
#1/#2 or FB2rrm4
#1/#2 with synthetic a1 or a2 pheromone (62) and could not observe any difference from the wild type (data not shown). To assay for postfusion defects in filament formation, we addressed whether wild-type strains were able to complement the mating defect of rrm4
strains. Therefore, we cospotted FB2 (a2b2) with serial dilutions of either FB1 (a1b1) or FB1rrm4
#1. The observation that there was no difference between the wild type and FB1rrm4
#1 (Fig. 4C) indicated that wild-type strains are able to rescue the mating defect of rrm4
strains, suggesting a role for Rrm4 in postfusion events. For verification, we deleted rrm4 in strain SG200 (a1mfa2bW2bE1) (8). SG200rrm4
#1 and -#2 strains produced only short filaments compared to the SG200 progenitor (Fig. 4D). These data illustrate that Rrm4 is not drastically affected in pheromone response but is crucial for b-dependent filamentation.
|
strains were mixed (FB1rrm4
#1 with FB2rrm4
#1 and FB1rrm4
#2 with FB2rrm4
#2) and more than 360 corn seedlings each were infected in five independent infection experiments (Materials and Methods; Fig. 5B). We observed that, in mixtures of rrm4
strains, the number of wilted or dead plants decreased approximately 1 order of magnitude in comparison to the wild type, whereas the number of plants without symptoms or with minor leaf tumors increased approximately threefold (Fig. 5B). Thus, rrm4
strains are still able to form tumors, but virulence is strongly affected.
|
| DISCUSSION |
|---|
|
|
|---|
RNA-binding proteins of the PUM type in U. maydis. U. maydis contains five PUM-type proteins, and it is very likely that these are all members of this particular family, since the primary sequence and iterative nature of Pumilio-like repeats are very conserved (70). This notion is supported by a phylogenetic analysis of fungal PUM-type proteins. Representatives from U. maydis are present in all four clusters and in an additional region of nonclustering sequences. All five proteins fall into separate branches of the dendrogram, indicating that they did not arise by recent gene duplication. This is likely the case for other proteins such as Puf1p and Puf2p from S. cerevisiae (Fig. 6A).
|
RNA-binding proteins of the KHD type in U. maydis.
According to our analysis, nine KHD-type proteins are present in U. maydis. This number is comparable to the numbers of KHD proteins from S. cerevisiae and S. pombe (10 and 8 proteins, respectively; http://smart.embl.de/), and most of the KHD proteins from U. maydis have related counterparts in S. cerevisiae (Table 1). The phenotypic analysis of khd1 to -4 deletion strains revealed that khd1
strains exhibit a cold-sensitive growth phenotype, while growth of khd2
strains was not affected by temperature. This indicates that Khd1 and Khd2 have distinct biological roles, which is in accordance with the fact that the region of high sequence similarity is restricted to their three KHD RNA-binding domains (20% overall sequence identity, but the first, second, and third KHD domains are 35, 43, and 52% identical in sequence, respectively). Loss of Khd4 affected growth rate, morphology, pheromone response, filamentous growth, and pathogenicity. At present it is not clear whether these pleiotropic defects are due to altered regulation of a number of target mRNAs or whether these phenotypes are connected through the same key component regulated by Khd4.
Comparison to sequences from other model fungi such as A. nidulans, N. crassa, Fusarium graminearum, Magnaporthe grisea, C. neoformans, and C. cinereus revealed that Khd4 is part of a novel group of KHD proteins that are related to multi-KHD proteins of the Vigilin type. Interestingly, members of this group are absent in S. pombe and S. cerevisiae (Fig. 6B). The closest relatives from S. cerevisiae and mammals are the multi-KHD proteins Scp160p and Vigilin, respectively, containing 14 KHDs each. This type of protein is also found in U. maydis (Khd5), as well as in the other fungi mentioned above (Fig. 6B). Interestingly, loss of Scp160p from S. cerevisiae results in abnormal morphology (72), and recently this RNA-binding protein has been described as being involved in mating, since it interacts with the heterotrimeric G protein Gpa1p during the pheromone response (22). It will be interesting to determine whether in U. maydis multi-KHD proteins such as Khd4 and Khd5 share functions in regulating morphogenesis and pheromone signaling.
RNA-binding proteins of the RRM type in U. maydis.
An automated feature search predicted the presence of 43 RRM proteins in U. maydis (http://www.broad.mit.edu/cgi-bin/annotation/fungi/ustilago_maydis/findfeatures.cgi), a number that is close to the 53 RRM proteins from S. cerevisiae (http://smart.embl.de/). Phenotyping strains carrying a deletion in a subset of 11 RRM protein-encoding genes revealed that only the loss of Rrm4 resulted in discernible phenotypes. This was surprising since we would have expected more RRM proteins to be important for developmental processes in U. maydis. For example, Rrm8 shares the same domain organization as well as similar flanking sequences to the negative regulator of development Nrd1 from S. pombe. Loss of Nrd1 function causes cells to initiate sexual development in the absence of nutrient starvation (64). In contrast, rrm8
strains from U. maydis were not drastically disturbed in mating or pathogenicity.
In the case of Rrm4, we showed that neither pheromone secretion nor pheromone-responsive conjugation tube formation was significantly disturbed. However, loss of Rrm4 negatively affects dikaryotic filaments as well as filamentation of a solopathogenic strain. This defect is likely to be the cause for the reduced disease symptoms, since it is conceivable that shorter filaments have difficulties in finding an appropriate entry site for plant penetration.
Rrm4 displays a novel domain organization that is so far unique in eukaryotes (Fig. 6C). At its C terminus, it contains a PABC domain that was initially found in poly(A)-binding proteins. However, this domain is also associated with a HECT domain (homologous to the E6-AP carboxyl terminus) found in ubiquitin-protein ligases and implicated as a protein-protein interaction interface (32, 39). The PABC domain in human poly(A)-binding protein contains four amino acids (F22, I25, A33, and K35) shown to form a deep hydrophobic pocket that could function as a peptide binding domain (32). Strikingly, all four residues are conserved in Rrm4. In its N terminus, Rrm4 contains three RRM domains with typical spacing known from ELAV-like proteins (3), such as Elav or Bruno from D. melanogaster, HuR from H. sapiens, and Pub1p from S. cerevisiae. All of these proteins contain two RRM domains juxtaposed with each other and separated from a third by a hinge region of variable length (Fig. 6C). Members of this family of proteins have been described to function as sequence-specific RNA-binding proteins regulating mRNA translation or stability (3). For example, Bruno represses translation of Oskar, a determinant for posterior body patterning in embryos from D. melanogaster, by binding to Bruno response elements in the 3'-untranslated region of the osk mRNA (47, 68, 71). HuR from humans has been implicated in regulating AUUUA-mediated mRNA decay (16, 46). Pub1p from S. cerevisiae recognizes stabilizer elements and thereby protects mRNAs that contain upstream ORFs from degradation by the nonsense-mediated decay pathway (52).
According to sequence similarity, we hypothesize that Rrm4 is a sequence-specific RNA-binding protein that might be involved in regulation of mRNA splicing, transport, translation, or stability. For two reasons, we do not believe that Rrm4 functions as poly(A)-binding protein. Rrm4 does not contain the typical four RRM domains in its N terminus, and in vitro binding experiments revealed that Rrm4 does not bind to poly(A) or poly(U) but recognizes specific sequences enriched by SELEX (C. Julius and M. Feldbrügge, unpublished observation). A much better candidate for a poly(A)-binding protein in U. maydis is Rrm12, since it contains the characteristic domain organization (Fig. 6C) and exhibits high sequence similarity to other PABPs over its entire length (comparison to Pab1p from S. cerevisiae results in a BLAST expect value of e132).
Conclusions. A great surprise from our study of RNA-binding proteins in U. maydis was the fact that, despite our extensive phenotypic analysis, the majority of gene replacement mutants did not exhibit any discernible phenotypes. A possible explanation would be functional redundancy of related RNA-binding proteins. Alternatively, it could be possible that our assays were not sensitive enough or that we did not test the affected phenotype. Nevertheless, we were able to identify 2 out of 18 RNA-binding proteins implicated in pathogenic development of U. maydis. One of them is a founding member of a new group of KHD proteins, and the other is a so far unique ELAV-like RRM protein. Investigation of the subcellular localization of these candidates in combination with identification of specifically recognized target RNAs will allow us to unravel the connection between the mode of action of these RNA-binding proteins and pathogenic development. For example, nuclear localization would indicate that they function in mRNA maturation or export. We expect to find within the group of target RNAs important determinants of pathogenic development whose expression is manifested at the posttranscriptional level. The involvement of RNA-binding proteins in regulation of pathogenic development supports our initial assumption predicting fundamental similarities to developmental programs of higher eukaryotes.
| ACKNOWLEDGMENTS |
|---|
This work was supported by Bayer CropScience.
| FOOTNOTES |
|---|
Supplemental material for this article may be found at http://ec.asm.org/. ![]()
| REFERENCES |
|---|
|
|
|---|