- THIS ARTICLE
-
Abstract
- Full Text (PDF)
- Alert me when this article is cited
- Alert me if a correction is posted
- SERVICES
- Similar articles in this journal
- Similar articles in PubMed
- Alert me to new issues of the journal
- Download to citation manager
- Reprints & Permissions
- CITING ARTICLES
- Citing Articles via HighWire
- Citing Articles via Google Scholar
- GOOGLE SCHOLAR
- Articles by Colaiácovo, M. P.
- Articles by Villeneuve, A. M.
- Search for Related Content
- PUBMED
- PubMed Citation
- Articles by Colaiácovo, M. P.
- Articles by Villeneuve, A. M.
A Targeted RNAi Screen for Genes Involved in Chromosome Morphogenesis and Nuclear Organization in the Caenorhabditis elegans Germline
M. P. Colaiácovoa, G. M. Stanfielda, K. C. Reddya, V. Reinke1,a, S. K. Kima, and A. M. Villeneuveaa Departments of Developmental Biology and Genetics, Stanford University School of Medicine, Stanford, California 94305-5329
Corresponding author: A. M. Villeneuve, Stanford University School of Medicine, Beckman Center, Rm. B309, 279 Campus Dr., Stanford, CA 94305-5329., villen{at}cmgm.stanford.edu (E-mail)
Communicating editor: B. MEYER
| ABSTRACT |
|---|
We have implemented a functional genomics strategy to identify genes involved in chromosome morphogenesis and nuclear organization during meiotic prophase in the Caenorhabditis elegans germline. This approach took advantage of a gene-expression survey that used DNA microarray technology to identify genes preferentially expressed in the germline. We defined a subset of 192 germline-enriched genes whose expression profiles were similar to those of previously identified meiosis genes and designed a screen to identify genes for which inhibition by RNA interference (RNAi) elicited defects in function or development of the germline. We obtained strong germline phenotypes for 27% of the genes tested, indicating that this targeted approach greatly enriched for genes that function in the germline. In addition to genes involved in key meiotic prophase events, we identified genes involved in meiotic progression, germline proliferation, and chromosome organization and/or segregation during mitotic growth.
IN metazoans, genetic information is transmitted from one generation to the next via a specialized cell lineage known as the germline. The germ cell lineage is set aside early during development, and its chromosomes are insulated from events that can occur in somatic lineages that disturb or interfere with genomic integrity such as programmed gene rearrangements, chromatin diminution, and telomere shortening. Further, most animals reproduce sexually, so after germ cells exit a proliferative state they undergo a specialized nuclear division program called meiosis, which enables diploid germ cells to generate haploid gametes.
The nematode Caenorhabditis elegans is an especially favorable system for investigating the mechanisms that govern the faithful transmission of the genome through the germline (![]()
![]()
![]()
![]()
![]()
![]()
Further advances in the identification of germline genes have been made possible by the completion of the C. elegans genome sequence (C. ELEGANS SEQUENCING CONSORTIUM 1998) and the development of DNA microarrays that can be used to survey gene expression on a genome-wide scale (![]()
![]()
![]()
![]()
Here we report the use of a targeted functional genomics strategy using RNA-mediated interference (RNAi) to investigate the function of a selected pool of the germline-intrinsic genes identified by ![]()
![]()
![]()
![]()
| MATERIALS AND METHODS |
|---|
Worm strains:
Only wild-type C. elegans from the Bristol N2 strain were utilized. Worms were cultured at 20° under standard conditions as described by ![]()
Primer pair design and dsRNA synthesis:
For 185 of the genes tested, PCR products amplified from genomic DNA were used as templates for dsRNA synthesis. The primers used were adapted from those used by ![]()
1 kb of coding sequence from the corresponding gene. For most primer pairs, T7 promoter sequences (TAATACGACTCACTATAG) were added to the 5' ends of each primer, allowing subsequent transcription and annealing steps to be performed as single-tube reactions. For three genes, C08F8.3, C49C3.7, and T07C4.3, one primer of the pair included the T7 promoter sequence while the other included a T3 promoter sequence (ATTAACCCTCACTAAAG). Comparison of these reactions with the yield obtained from the same reactions where both primers carried a T7 promoter and were synthesized as single-tube reactions indicated no difference in yield.
Primary PCR reactions (50 µl) were done in 96-well plates using
0.1 µg template genomic DNA and 0.5 µM primers. To optimize the yield of products with full-length T7 promoters, we used 2 µl of the primary PCR reaction as template for a secondary amplification using T7 primers (0.5 µM). Both primary and secondary PCR products were assessed on 1% agarose gels to confirm size and yield. Products were purified with a 96-well QIAGEN (Valencia, CA) PCR purification kit and eluted in 30 µl. All the eluted DNA was used in a 100-µl transcription reaction with T7 RNA polymerase (Promega, Madison, WI), followed by DNAse I treatment. RNA samples were then purified (QIAGEN) and eluted in 30 µl; an aliquot (ssRNA) was removed prior to annealing samples 10 min at 68° followed by 30 min at 37°. dsRNA concentrations ranged from 1 to 5 µg/µl. To score for shifts in mobility for the dsRNA, ssRNA was run parallel to dsRNA on 1.4% agarose gels.
dsRNA synthesis for seven additional genes (F57B10.4, R12B2.4, T06E4.1, ZK1055.1, F26D2.2, F39H2.4, and F56A3.4) tested in our pilot screen and included in this article was done following the protocol described by FIRE et al. (1998) using cDNA clones yk176b1, yk428e9, yk204f9, yk362d7, yk414g2, yk252b7, and yk295b7, respectively, provided by Dr. Yuji Kohara, NIG (Japan).
Injections and phenotypic analysis:
Each dsRNA was injected into one or both gonad arms of 15 young adult hermaphrodites (P0's), which were plated individually and transferred serially to fresh plates to collect broods of F1's (Fig 1). In most cases, four F1's per injected P0 (arising from embryos laid 1248 hr post-injection) were picked to individual plates, left to lay eggs for 2436 hr, and transferred to new plates. P0's were kept for cytological analysis. In some cases, RNAi led to death of a substantial fraction of F1's, and <60 viable F1's were available. Because our primary goal was to identify genes involved in meiosis, plates containing progeny produced by P0's (during 1248 hr post-injection) and the initial plates containing progeny of F1's were scored for (1) a high incidence of males (or Him) phenotype (![]()
|
P0's or F1's that displayed any of the above phenotypes were subjected to cytological analysis. Undissected whole worms were fixed with Carnoy's fixative and stained with 4',6-diamidino-2-phenylindole (DAPI) as described in ![]()
List of genes that elicited no defect upon RNAi:
3R5.1, B0001.2, B0001.3, C01G5.8, C08C3.2, C08F8.3, C13F10.6, C16A11.3, C16A11.4, C16C8.4, C18H2.2, C29H12.5, C34G6.5, C50C3.8, C56A3.5, D1081.7, D2030.8, F01G4.4, F07H5.10, F10B5.5, F13G3.6, F14D2.8, F23B12.8 (klp-14), F26A1.1, F26H9.4, F28F8.6, F30F8.3, F32E10.2, F33H2.1, F35G12.12, F38A5.13, F38B7.7 a.k.a. H12C20.2a (pms-2), F39H2.1, F45E4.10 (gfi-4), F49E8.7, F52C9.7, F53F4.14, F54D5.9, F57A10.4, F57B10.4, F57B10.6, K07A1.1, K07H8.10, K08F4.2, K08F4.3, K10D2.1, K12D12.5, M03E7.5, R05D11.5, R05D3.11, R06C7.2, R06C7.9, R07B7.2, R09B3.1 (exo-3), R10D12.14, R11A8.2, R74.8, T01G9.4 (kup-2), T05H10.2 (apn-1), T06E4.1 (hcp-2), T07C12.3, T07C4.3a and b, T08B2.11, T09A5.8, T13F2.6, T19B10.8, T23B12.4, T24A6.1, T24D1.3, T24G10.2, T25G3.3, T26A5.2, T26A5.5, W02D3.10, W05F2.2a and b, Y102E9.2, Y17G7A.1 (hmg-12), Y32B12B.2, ZC155.3, ZC410.3, ZK1055.1 (hcp-1), ZK1307.9, ZK856.12.
Statistical analysis:
Statistical analyses (Fisher's exact test, chi-square test for independence, and Mann-Whitney test) were performed using the InStat software package (GraphPad.com).
| RESULTS AND DISCUSSION |
|---|
Rationale of the RNAi screen:
A major goal of our research program is to investigate the mechanisms responsible for faithful segregation of homologous chromosomes during meiosis. Because previously identified components of the meiotic machinery are expressed preferentially in the germline, we reasoned that it should be possible to discover additional components by screening among genes exhibiting germline-enriched gene expression to identify those for which inhibition by RNAi elicits defects in meiotic chromosome morphology or behavior. Meiotic defects can be detected indirectly by examining the products of meiosis through progeny testing of affected individuals; in C. elegans, defects in meiotic chromosome segregation are manifested by a high frequency of XO male progeny and inviable aneuploid embryos among the self-progeny of XX hermaphrodites (the Him phenotype). Alternatively, meiotic defects can be detected more directly, by cytological examination of DAPI-stained meiotic prophase chromosomes in the germlines of affected individuals (Fig 2). Any defects leading to an absence or reduced frequency of crossover recombination are readily detected at diakinesis, the last stage of meiotic prophase: whereas wild-type nuclei at this stage contain six discrete DAPI-stained bodies, each corresponding to a pair of homologs attached by a chiasma, a deficit in crossing over results in the presence of up to 12 univalent chromosomes that are unattached to their homologous partners. Further, failure of chromosomes to reorganize or align lengthwise with their homologs at earlier stages of meiotic prophase can be detected by an altered appearance of the DAPI-stained chromatin in more distal parts of the germline.
|
Because RNAi effects are often most pronounced in the F1 progeny of treated animals, the former means of assessment of meiotic defects requires examination of progeny two generations after dsRNA administration, whereas the latter requires microscopic imaging of the germlines of F1 animals after they have reached adulthood. Thus since an effective screening procedure for identifying meiotic defects would involve a nontrivial investment of effort for each gene tested, we also chose to document the effects of RNAi on several additional aspects of germline development, organization, and function.
Selection of 192 candidate genes:
A subset of genes from the germline-intrinsic list of ![]()
![]()
![]()
![]()
![]()
To a large extent, we excluded from our list genes for which the biological role was already known or which we knew to be under investigation by other laboratories. We also used sequence information to make strategic decisions to exclude genes we thought unlikely to play a direct or specific role in meiotic chromosome behavior (e.g., RNA polymerase, ribosomal proteins). Including both the initial small-scale screen and the main screen, we screened a total of 192 genes.
Design and validation of the screening strategy:
To develop a screening strategy that would be successful in identifying genes involved in diverse meiotic prophase events, we conducted pilot RNAi experiments using six genes previously implicated in various aspects of the meiotic prophase program. We used the following genes: him-3, which encodes a meiosis-specific component of the chromosome axis similar to Hop1, an axial element/lateral element protein from budding yeast (![]()
![]()
![]()
![]()
In our initial experiments with the above-described positive control genes, we injected dsRNA at concentrations of
50250 ng/µl. Although this concentration elicited a robust meiotic-defective phenotype at high frequency for some genes, for other genes this concentration either elicited only weak phenotypes or produced a strong phenotype only at a low frequency. By increasing the concentration of the dsRNA to 15 µg/µl (as in ![]()
In addition to defining an appropriate dsRNA concentration and the numbers of worms we would be handling during the screen, it was important to define a time window during which the F1 progeny produced by the injected P0 animals would be most likely to exhibit germline defects. By transferring the injected P0's to new plates at 12-hr intervals, we were able to monitor F1's generated at successive intervals both for their plate phenotypes (frequencies of males and dead embryos produced) and for their cytological phenotypes (altered chromosome morphology or organization in germline nuclei). We observed the strongest effects of RNAi in F1's derived from embryos laid between 2448 hr post-injection, while F1's laid before or after that time period exhibited a weaker RNAi effect for some of the targeted genes.
For the genes used in our design and validation phase, reduction or loss of gene function leads to chromosome missegregation but does not prevent completion of the meiotic program, so embryos are produced (albeit many are aneuploid and inviable). We also wished to identify genes that are crucial for initiating or completing the meiotic program, so in addition to the Him and cytological phenotypes described above, in the screening phase F1's were also scored for sterility (failure to produce any embryos).
Classification of genes:
For 57% of the 192 genes tested, phenotypic defects were observed following RNAi. These fall into four broad classes, as summarized in Table 1. We have used various descriptor terms to indicate the phenotypes observed. For many of the genes, RNAi elicited a combination of phenotypes, so more than one descriptor term is applicable. Where we have used multiple descriptors, the first descriptor generally denotes what we considered to be either the primary defect or that which most directly reflects the germline function for the gene in question. This may indicate a phenotype that preceded all others observed and/or the phenotype most frequently observed among the adult F1's. In some cases, additional descriptors refer to phenotypes present among subpopulations of the affected F1's. In other cases more general descriptors (e.g., sterile) were also applied after more specific primary descriptors (e.g., tumorous germline) both to indicate a phenotype that could be assessed at lower resolution and to facilitate searching the database containing a compilation of our results (the full database is available by ftp at http://villeneuveRNAi.stanford.edu; user ID and password: villeneuveRNAi; see CONCLUSION).
|
Strong germline class
Fifty-one genes (27%) were classified as "strong germline" genes because RNAi elicited clear defects in function and/or development of the germline. Genes in this class were listed under the following primary descriptors.
Meiotic: This descriptor was applied when RNAi elicited defects in meiotic chromosome segregation generating anywhere from medium to high levels of the Him phenotype (see MATERIALS AND METHODS) and/or resulted in the presence of achiasmate chromosomes at the diakinesis stage of meiotic prophase. In some cases the appearance of nuclei at earlier stages of meiotic prophase was also affected.
Included under this primary descriptor are three genes (syp-1, syp-2, and syp-3) encoding proteins predicted to have extended coiled-coil domains. RNAi for each of these genes led to a high frequency of univalent chromosomes at the diakinesis stage of meiotic prophase (Fig 2K) and defects in chromosomal organization earlier in prophase, including a persistence of the polarized nuclear organization normally seen only at the onset of homologous chromosome pairing in early meiotic prophase (extended transition zone phenotype; Fig 2F). All three genes have now been matched with meiotic mutants identified in our genetic screens (A. J. MACQUEEN, M. P. COLAIÁCOVO, J. ENGEBRECHT, K. C. REDDY and A. M. VILLENEUVE, unpublished data), and an in-depth analysis of their functional roles will be reported elsewhere.
The rad-50 gene encodes a conserved protein whose fungal orthologs play important roles in both meiotic recombination and DNA repair, acting in a complex together with the Mre11 protein (![]()
![]()
Chromosome segregation defects were also accompanied by variable numbers of achiasmate chromosomes in diakinesis-stage oocytes for several other genes tested. RNAi for one gene, T09E8.2, appeared to preferentially affect the segregation of the X chromosomes; we have since found that T09E8.2 corresponds to the him-17 gene (defined by multiple mutant alleles; K. REDDY, J. HODGKIN and A. M. VILLENEUVE, unpublished data) and are currently investigating its function in more detail. For F59A1.7, affected worms exhibiting a high frequency of achiasmate chromosomes at diakinesis also had abnormally large nuclei at the pachytene stage, earlier in meiotic prophase. For C05D2.5, a medium Him phenotype in affected F1's was accompanied by defects in meiotic progression (see below) and achiasmate chromosomes at diakinesis in a subset of animals; a Him phenotype had been seen previously for C05D2.5 when gene function was inhibited by transgene-mediated cosuppression or by RNAi (L. KUERVERS and D. BAILLIE, personal communication).
Meiotic defects were also observed following RNAi for M04F3.1, which encodes the C. elegans ortholog of RPA2/Rfa2, a subunit of the eukaryotic single-stranded DNA binding protein shown in other systems to function in DNA replication, repair, and recombination (![]()
![]()
For all but one of the genes to which the "meiotic" descriptor was applied, defects were evident at one or more stage(s) of meiotic prophase in at least a subset of affected animals. The exception was klp-16, which encodes the C. elegans ortholog of kinesin-related motor proteins Drosophila melanogaster Ncd and Saccharomyces cerevisiae Kar3 (![]()
![]()
![]()
![]()
![]()
We also detected a meiotic phenotype for K12H4.8, recently named dcr-1 (![]()
![]()
![]()
Meiotic progression:
For several genes the spatial/temporal gradient of meiotic prophase was altered by RNAi, as evidenced either by a pachytene arrest or by the presence of reduced numbers of nuclei representing particular substages in meiotic prophase as if meiotic progression were accelerated (see Fig 2G). While affected P0's laid low-to-medium levels of dead eggs, F1's were sterile. The pachytene arrest and accompanying sterility phenotype are similar to phenotypes reported for mutations in genes encoding components of the mitogen-activated protein kinase (MAPK) signaling pathway required for exit from pachytene (![]()
![]()
Abnormal gonad structure: For several genes, RNAi led to a variety of abnormalities in the shape or size of the gonads in the affected individuals (Fig 2H). In addition, we observed intermixing of nuclei at different stages of meiotic prophase rather than a clear spatial/temporal gradient of meiotic nuclei. In such cases it was not possible to discern from our analysis whether the abnormality in gonad structure was the cause of the aberrant organization of meiotic stages or whether these two phenotypes represent separate effects of the RNAi. Further, in some cases a subset of nuclei appeared to have degenerated or fragmented. Such phenotypes were sometimes associated with embryonic lethality or sterility in the F1.
Mitotic:
This descriptor was applied to a broad spectrum of phenotypes in which the appearance of nuclei was abnormal in the premeiotic region of the germline. (This region includes both mitotically proliferating nuclei as well as nuclei in the G1 and S phases immediately preceding meiotic prophase.) For many genes in this class, we observed variation in the size, shape, and brightness of the DAPI signals in germline nuclei of affected adults, likely indicating defects in chromosome segregation during mitotic growth (Fig 2I). Examples of genes in this subclass include R12B2.4, recently reported by ![]()
![]()
![]()
![]()
![]()
![]()
![]()
![]()
(involved in both replication and repair), and C54G10.2/rfc-1, which encodes the large subunit of replication factor C (the clamp-loading protein; ![]()
For most genes in the mitotic class, the majority of F1's produced by treated P0's died as embryos, arrested as larvae, or became thin, sterile uncoordinated adults (a syndrome that reflects defects in postembryonic cell divisions; ![]()
Rudimentary gonads:
For some genes, the treated P0 animals produced sterile F1 progeny that had very few or no germline nuclei and gonad arms that were very small or missing. In most cases, this likely reflected a defect in early proliferation of the germline, as seen in mutants such as glp-4 (![]()
![]()
![]()
Tumorous germline:
RNAi for some genes generated a tumorous germline phenotype. For these genes, affected F1's were usually sterile, and occasional escapers that presumably avoided the full effect of RNAi laid very small broods. Cytologically, the gonad arms of most affected F1's had normally shaped, mitotically dividing nuclei throughout their entire lengths, with no evidence of entry into meiotic prophase. In a few cases some apparent meiotic prophase nuclei were also observed, but no diakinesis nuclei were evident. These tumorous germline phenotypes are reminiscent of those observed in loss-of-function gld-1 and gld-1; gld-2 double mutants and in glp-1 gain-of-function mutants, which are defective in regulating exit from the mitotic cell cycle and entrance into meiosis (![]()
![]()
![]()
Low-frequency phenotypes
For a significant fraction of genes tested in our screen (19%), we observed F1 animals exhibiting robust germline-defective phenotypes only at very low frequencies (17 of 60 F1's scored were affected). In the rare affected animals, however, the phenotypes observed were quite strong. For example, for most of the genes in the "low-frequency meiotic" subclass, a high frequency of achiasmate chromosomes was seen in all oocytes of the affected animals.
For nine genes for which low-frequency germline phenotypes were detected, essentially all of the F1's produced during the later broods of the injected P0's died as embryos, indicating that these genes provide essential functions. In these cases, it seems likely that most of the F1's that survived to adulthood were in fact escapers from the effects of RNAi. In such cases, detection of a small subset of the viable F1's that have strong germline-defective phenotypes may reflect a small "window of opportunity" in which reduction of gene activity by RNAi was not sufficient (either in time or in amount) to preclude viability but was sufficient to interfere with germ cell function. This category includes several genes encoding proteins with orthologs or paralogs that have been implicated in chromosome segregation or cell-cycle progression: F08H9.1/coh-3, which encodes a protein related to the Scc1 subunit of cohesins (![]()
![]()
![]()
For 29 genes, a "low-frequency" phenotype was the only phenotype observed following RNAi. For 28 of these genes, the defects observed were in the germlines of affected F1's, while for one (C49C3.7), the phenotype was embryonic lethality in the progeny of the rare affected F1's. For the largest single subclass (17 genes), affected F1's exhibited a high frequency of achiasmate chromosomes in diakinesis-stage oocytes, a phenotype that is normally a reliable diagnostic of defects in meiotic prophase. It was unexpected that we would see a very strong germline phenotype in only a few F1's for such a large fraction of the genes tested. It is possible that these low-frequency phenotypes reflect a nonspecific effect of the RNAi procedure or a baseline spontaneous occurrence of these phenotypes detected because of the large numbers of animals analyzed. If either of these were the case, then it would be appropriate to group the "low-frequency-only" genes into a single class together with the genes for which no defect was observed following RNAi. While the basis for these low- frequency phenotypes remains unclear, several considerations argue for a specific effect of RNAi. First, we found that the low-frequency-only and "no defect" classes differed significantly with respect to the fraction of genes in the class having probable orthologs in other species (see below); this finding is inconsistent with the idea that these two groups constitute a single class. Second, a further prediction of a "single-class" model (in which the incidence of the phenotype reflects random occurrence of a low-frequency event) is that the proportions of genes for which zero, one, two, or more F1's exhibit a particular low-frequency phenotype should follow the Poisson distribution. This expectation is not borne out by the existing data set. We examined the distribution of worms exhibiting the "achiasmate diakinesis" phenotype among genes in the consolidated low-frequency-only or no defect class and found that it differed significantly from the Poisson distribution (P < 0.04).
Thus we infer that for at least a subset of genes falling into the low-frequency-only category, the phenotype observed was not merely random or nonspecific, but rather provided a clue about the biological role for the gene tested. Some of these genes may be refractory to RNAi and, for some genes, a more reliable phenotype might be obtained by delivering dsRNA by another means or at a different time or by targeting a different part of the gene. Further, although RNAi may not prove to be an effective approach to investigate their functions, the low-frequency class of genes may serve as useful candidates to facilitate molecular identification of genes identified by mutational analysis.
While we have just argued that at least a subset of genes in the low-frequency class probably have legitimate functions in meiosis, we do not exclude the possibility that some of these low-frequency events could have been a nonspecific consequence of the RNAi treatment. It has been previously shown that loss of function of ego-1, which encodes an RNA-dependent RNA polymerase involved in the mechanism of RNAi in the germline, also causes pleiotropic effects including defective oogenesis, altered meiotic progression, and some achiasmate chromosomes (![]()
Nongermline class
Of the genes tested, 10% showed a "nongermline" phenotype such as larval arrest, embryonic lethality, or unco-ordinated movement (Unc) without any evident germline-associated phenotype. For such genes it is possible there is no role for the gene product in the germline per se and that the germline expression of the gene reflects transcription solely for the purpose of deposition of the corresponding mRNA or protein in the embryo. Alternatively, it is possible that some of these gene products also function in the germline, but that the early arrest or lethality elicited by RNAi precluded our ability to uncover such roles in F1 animals. In some cases the embryonic lethal descriptor was accompanied by the secondary descriptor "mitotic-emb"; this designation indicated that DAPI staining revealed evidence of abnormal cell or nuclear divisions in the embryos, such as nuclei with abnormal DNA content and/or chromatin bridges visible between nuclei.
Sequence conservation and detection of RNAi phenotypes: Genetic redundancy is regularly cited as a possible explanation for failure to detect a phenotype following targeting of a gene. The simplest form of redundancy is the presence in the genome of a second gene whose sequence is closely related to that of the gene in question. We assessed whether this type of redundancy might have contributed significantly to cases where no defect was detected or where phenotypes were detected only at low frequency. Specifically, we asked whether genes with close paralogs were overrepresented among genes in the no defect or low-frequency-only classes compared with genes for which RNAi elicited strong phenotypes. For this analysis we considered a gene to have a close worm paralog when one or more additional C. elegans genes encoded a protein with an overall level of extended similarity greater than or equal to that of the closest homologs from other species. Using these criteria, we found that genes with close paralogs accounted for 27% of no defect genes, 21% of low-frequency-only genes, 22% of "strong germline" genes, and 23% of "all strong phenotype" genes; there were no significant differences between any of these groups regarding the proportion of genes with and without close paralogs (pairwise comparisons were conducted using Fisher's exact test).
We then took into account whether paralog genes were likely to be subject to co-RNAi in our experiments. For this analysis a paralog gene was considered "likely" to be subject to co-RNAi if it contained multiple stretches considerably greater than 23 nucleotides (nt) in length with 100% nucleotide identity with the input dsRNA; co-RNAi was considered "plausible" (but uncertain) for paralog genes that contained only one to three stretches of identity between 23 and 35 nt in length or for one short gene that had multiple close paralogs with significant stretches of identity. (These criteria were based on experimental data from M. MONTGOMERY, personal communication). We discovered that genes for which co-RNAi was likely or plausible were significantly overrepresented among the paralogs of strong germline genes and all strong phenotype genes compared with those of no defect genes (P = 0.013 for the strong germline vs. no defect comparison and P = 0.007 for the all strong phenotype vs. no defect comparison using chi-square test for independence). When we then examined the whole gene list with regard to whether genes had a C. elegans paralog that would not be subject to co-RNAi, we found that such genes were significantly overrepresented among the no defect genes (23%) compared with the strong germline genes (8%) and all strong phenotype genes (9%) (P = 0.033 for the strong germline vs. no defect comparison and P = 0.018 for the all strong phenotype vs. no defect comparison using Fisher's exact test). These analyses suggest that genetic redundancy probably does account for a subset of the cases in which no phenotype was elicited by RNAi.
We next considered whether there might be differences among the classes with respect to the fraction of genes with probable orthologs in other species. We found that RNAi was significantly more likely to elicit a detectable phenotype for genes with apparent functional conservation across species than for genes lacking probable orthologs. The proportion of genes with probable orthologs in one or more of the species referenced in the Proteome database (WORMPD Bioknowledge Library from Incyte Genomics; ![]()
![]()
![]()
![]()
Gene-expression profiles and RNAi phenotypes:
The genes tested in our screen were chosen from a list of genes exhibiting germline-enriched expression in the microarray analysis of ![]()
To ask whether there might be any correlation between the expression level of a gene and the phenotypic class to which it belonged, we compared the distributions of mean gene expression (MGE) values between the different classes. The MGE value for a given gene is a normalized representation of its expression level. For each gene in a given experiment, we calculated the ratio of the raw expression level for the individual gene to the average per gene-expression level for all genes on the microarray; MGE values for each gene were derived from the published data set (![]()
When a nongermline phenotype (particularly embryonic lethality) was the major or only phenotype detected for a gene, it might be the case that transcripts and/or protein products expressed in the germline were there primarily for deposition in the oocyte to support development of the embryo. It seemed plausible that the ratio of expression in worms producing only oocytes to worms producing only sperm might tend to be higher for such genes compared with genes for which strong phenotypes were observed in the germline itself, but we failed to find any statistical support for this possibility. Comparison of the distributions of oocyte/sperm ratios between the strong germline and nongermline classes did not reveal any significant difference between these classes (Mann-Whitney test).
Finally, we asked whether the level of induction of a gene's expression in the germline might be a useful predictor of a gene's phenotypic class. We conducted pairwise comparisons of the distributions of "fold-induction" values between the strong germline class and the low-frequency, nongermline, and no defect classes, but once again found no significant differences. Thus for the set of genes tested in our screen, we failed to find any expression criteria that could be used to further subdivide the genes in a way that would be predictive of biological function.
Comparison with previous RNAi screens:
The results of several RNAi screens have been reported to date. In two of these screens, the set of genes tested was defined by chromosomal position; screens of 96% of the predicted genes on chromosome III (![]()
![]()
![]()
![]()
![]()
In our screen of 192 genes representing a defined subset of germline-enriched genes, we detected some type of phenotype for 57% of genes tested. Even if we eliminate from consideration those genes for which a phenotype was detected only at low frequency, we still saw phenotypes for 42% of genes tested, a "hit rate" significantly higher than that for any of the previously reported screens. Moreover, we detected strong germline phenotypes for 27% of tested genes, indicating that our focused strategy was highly effective in identifying genes required for normal development and function of the germline. We attribute our high success rate in identifying germline phenotypes to at least two key factors. First, we focused on genes expressed in the relevant tissue at the relevant time. Second, our multigenerational screening strategy allowed us both to identify phenotypes revealed preferentially in later generations and/or later time points than those examined in some screens and to uncover germline phenotypes that may have been masked by embryonic lethality among F1 animals.
The efficacy of our screening strategy is emphasized when we examine the 105 genes that were tested both in our screen and in one or more of the other screens mentioned above (Table 1). For 69 of these genes, similar phenotypes were detected in our screen and the other screens. However, we detected RNAi phenotypes for a total of 27 genes that were classified as "no phenotype" in the other screens; most of these were in the low-frequency-only class, but 12 showed strong phenotypes. Further, we also observed additional phenotypes associated with embryonic lethality for 9 genes identified only as embryonic lethal in the other screens; for 7 of these, RNAi elicited strong germline phenotypes in our hands. In contrast, there were only 3 genes among the 38 for which a phenotype that was not detected in our screen was detected in other screens.
| CONCLUSION |
|---|
Our work has shown that the combination of a preselected pool of target genes and a multigenerational screening strategy has a high success rate for identifying genes that function in the C. elegans germline. We achieved our initial goal of identifying several new genes involved in key meiotic prophase events. Moreover, these efforts, in synergy with other work in the laboratory, have greatly expedited the molecular identification of genes we had previously defined by mutational analysis (A. J. MACQUEEN, M. P. COLAIÁCOVO, J. ENGEBRECHT, K. C. REDDY and A. M. VILLENEUVE, unpublished data). We also identified genes involved in many other aspects of germline development and function, including meiotic progression, germline proliferation, and chromosome organization and/or segregation during mitotic growth. Thus the information obtained in our screen will be useful not only for understanding germline function per se, but also for understanding chromosome metabolism and cell division in general. To provide broad access to detailed information about each of the genes analyzed in the screen, a FilemakerPro file containing our full database has been made available by ftp at ftp://villeneuveRNAi.stanford.edu (user ID and password: villeneuveRNAi). The database includes extensive descriptions of our phenotypic observations, sample images of DAPI-stained worms, and links to Wormbase, WormPD, and the germline microarray data from ![]()
| FOOTNOTES |
|---|
1 Present address: Department of Genetics, Yale University School of Medicine, New Haven, CT 06520. ![]()
| ACKNOWLEDGMENTS |
|---|
We are grateful to J. Wang for assistance with primer synthesis and J. Lund for assistance with installing database links. We thank Y. Kohara of the National Institute of Genetics (Japan) for cDNA clones. This work was supported by National Institutes of Health grant GM-53804 to A.M.V., a Bank of America-Giannini Foundation Postdoctoral Fellowship to M.P.C., a Katharine McCormick Postdoctoral Fellowship and a Cancer Research Fund of the Damon Runyon-Walter Winchell Foundation Postdoctoral Fellowship (DRG-1636) to G.M.S., and a National Science Foundation Predoctoral Fellowship to K.C.R.
Manuscript received February 22, 2002; Accepted for publication May 24, 2002.
| LITERATURE CITED |
|---|
ALBERTSON, D. G., J. E. SULSTON, and J. G. WHITE, 1978 Cell cycling and DNA replication in a mutant blocked in cell division in the nematode Caenorhabditis elegans.. Dev. Biol. 63:165-178.[Medline]
ALBERTSON, D. G., A. M. ROSE and A. M. VILLENEUVE, 1997 Chromosome organization, mitosis, and meiosis, pp. 4778 in C. elegans II, edited by D. L. RIDDLE, T. BLUMENTHAL, B. J. MEYER and J. R. PRIESS. Cold Spring Harbor Laboratory Press, Cold Spring Harbor, NY.
BAI, C., P. SEN, K. HOFMANN, L. MA, and M. GOEBL et al., 1996 SKP1 connects cell cycle regulators to the ubiquitin proteolysis machinery through a novel motif, the F-box. Cell 86:263-274.[Medline]
BEANAN, M. J. and S. STROME, 1992 Characterization of a germline proliferation mutation in C. elegans.. Development 116:755-766.[Abstract]
BERGER, B., D. B. WILSON, E. WOLF, T. TONCHEV, and M. MILLA et al., 1995 Predicting coiled coils by use of pairwise residue correlation. Proc. Natl. Acad. Sci. USA 92:8259-8263.
BERRY, L. W., B. WESTLUND, and T. SCHEDL, 1997 Germ-line tumor formation caused by activation of glp-1, a Caenorhabditis elegans member of the Notch family of receptors. Development 124:925-936.[Abstract]
BLELLOCH, R. and J. KIMBLE, 1999 Control of organ shape by a secreted metalloprotease in the nematode Caenorhabditis elegans.. Nature 399:586-590.[Medline]
BRENNER, S., 1974 The genetics of Caenorhabditis elegans.. Genetics 77:71-94.
Genome sequence of the nematode C. elegans: a platform for investigating biology. (1998) Science 282:2012-2018.
CHIN, G. M. and A. M. VILLENEUVE, 2001 C. elegans mre-11 is required for meiotic recombination and DNA repair but is dispensable for the meiotic G(2) DNA damage checkpoint. Genes Dev. 15:522-534.
CHURCH, D. L., K. L. GUAN, and E. J. LAMBIE, 1995 Three genes of the MAP kinase cascade, mek-2, mpk-1/sur-1 and let-60 ras, are required for meiotic cell cycle progression in Caenorhabditis elegans.. Development 121:2525-2535.[Abstract]
COHEN-FIX, O., 2001 The making and breaking of sister chromatid cohesion. Cell 106:137-140.[Medline]
COSTANZO, M. C., M. E. CRAWFORD, J. E. HIRSCHMAN, J. E. KRANZ, and P. OLSEN et al., 2001 YPDTM, PombePDTM, and WormPDTM: model organism volumes of the BioKnowledge® library, an integrated resource for protein information. Nucleic Acids Res. 29:75-79.
DERNBURG, A. F., K. MCDONALD, G. MOULDER, R. BARSTEAD, and M. DRESSER et al., 1998 Meiotic recombination in C. elegans initiates by a conserved mechanism and is dispensable for homologous chromosome synapsis. Cell 94:387-398.[Medline]
DERNBURG, A. F., J. ZALEVSKY, M. P. COLAIÁCOVO, and A. M. VILLENEUVE, 2000 Transgene-mediated cosuppression in the C. elegans germline. Genes Dev. 14:1578-1583.
ENDOW, S. A., S. HENIKOFF, and L. SOLER-NIEDZIELA, 1990 Mediation of meiotic and early mitotic chromosome segregation in Drosophila by a protein related to kinesin. Nature 345:81-83.[Medline]
FLORES-ROZAS, H. and R. D. KOLODNER, 2000 Links between replication, recombination and genome instability in eukaryotes. Trends Biochem. Sci. 25:196-200.[Medline]
FRANCIS, R., M. K. BARTON, J. KIMBLE, and T. SCHEDL, 1995 gld-1, a tumor suppressor gene required for oocyte development in Caenorhabditis elegans.. Genetics 139:579-606.[Abstract]
FRASER, A. G., R. S. KAMATH, P. ZIPPERLEN, M. MARTINEZ-CAMPOS, and M. SOHRMANN et al., 2000 Functional genomic analysis of C. elegans chromosome I by systematic RNA interference. Nature 408:325-330.[Medline]
FRIEDMAN, L., S. SANTA ANNA-ARRIOLA, J. HODGKIN, and J. KIMBLE, 2000 gon-4, a cell lineage regulator required for gonadogenesis in Caenorhabditis elegans.. Dev. Biol. 228:350-362.[Medline]
GARTNER, A., S. MILSTEIN, S. AHMED, J. HODGKIN, and M. O. HENGARTNER, 2000 A conserved checkpoint pathway mediates DNA damage-induced apoptosis and cell cycle arrest in C. elegans.. Mol. Cell 5:435-443.[Medline]
GÖNCZY, P., G. ECHEVERRI, K. OEGEMA, A. COULSON, and S. J. JONES et al., 2000 Functional genomic analysis of cell division in C. elegans using RNAi of genes on chromosome III. Nature 408:331-336.[Medline]
GRISHOK, A., H. TABARA, and C. C. MELLO, 2000 Genetic requirements for inheritance of RNAi in C. elegans.. Science 287:2494-2497.
GRISHOK, A., A. E. PASQUINELLI, D. CONTE, N. LI, and S. PARRISH et al., 2001 Genes and mechanisms related to RNA interference regulate expression of the small temporal RNAs that control C. elegans developmental timing. Cell 106:23-34.[Medline]
HABER, J. E., 1998 The many interfaces of Mre11. Cell 95:583-586.[Medline]
HANAZAWA, M., M. MOCHII, N. UENO, Y. KOHARA, and Y. IINO, 2001 Use of cDNA subtraction and RNA interference screens in combination reveals genes required for germline development in Caenorhabditis elegans.. Proc. Natl. Acad. Sci. USA 98:8686-8691.
HODGKIN, J., H. R. HORVITZ, and S. BRENNER, 1979 Nondisjunction mutants of the nematode Caenorhabditis elegans.. Genetics 91:67-94.
HOWE, M., K. L. MCDONALD, D. G. ALBERTSON, and B. J. MEYER, 2001 HIM-10 is required for kinetochore structure and function on Caenorhabditis elegans holocentric chromosomes. J. Cell Biol. 153:1227-1238.
HUBBARD, E. J. and D. GREENSTEIN, 2000 The Caenorhabditis elegans gonad: a test tube for cell and developmental biology. Dev. Dyn. 218:2-22.[Medline]
JANTSCH-PLUNGER, V., P. GÖNCZY, A. ROMANO, H. SCHNABEL, and D. HAMILL et al., 2000 CYK-4: a Rho family gtpase activating protein (GAP) required for central spindle formation and cytokinesis. J. Cell Biol. 149:1391-1404.
KADYK, L. C. and J. KIMBLE, 1998 Genetic regulation of entry into meiosis in Caenorhabditis elegans.. Development 125:1803-1813.[Abstract]
KELLY, K. O., A. F. DERNBURG, G. M. STANFIELD, and A. M. VILLENEUVE, 2000 Caenorhabditis elegans msh-5 is required for both normal and radiation-induced meiotic crossing over but not for completion of meiosis. Genetics 156:617-630.
KETTING, R. F., S. E. FISCHER, E. BERNSTEIN, T. SIJEN, and G. J. HANNON et al., 2001 Dicer functions in RNA interference and in synthesis of small RNA involved in developmental timing in C. elegans.. Genes Dev. 15:2654-2659.
KNIGHT, S. W. and B. L. BASS, 2001 A role for the RNase III enzyme DCR-1 in RNA interference and germline development in Caenorhabditis elegans.. Science 293:2269-2271.
KOWALCZYKOWSKI, S. C., 2000 Initiation of genetic recombination and recombination-dependent replication. Trends Biochem. Sci. 25:156-165.[Medline]
LUPAS, A., M. VAN DYKE, and J. STOCK, 1991 Predicting coiled coils from protein sequences. Science 252:1162-1164.[Medline]
MACQUEEN, A. J. and A. M. VILLENEUVE, 2001 Nuclear reorganization and homologous chromosome pairing during meiotic propha

