Mutational accessibility of essential genes on chromosome I(left) in Caenorhabditis elegans

Identification of Essential Genes in Caenorhabditis elegans with Lethal Mutations Maintained by Genetic Balancers

Genetic balancer systems, which allow effective capture and maintenance of lethal mutations stably, play an important role in identifying essential genes. Whole-genome sequencing (WGS) followed by bioinformatics analysis, combined with genetic mapping data analysis, allows for an efficient and economical means of identifying genomic mutations in essential genes. Using this approach, we successfully identified 104 essential genes on ChrI, ChrIII, and ChrV in C. elegans. In this report, we described a protocol that sequences the genome of prebalanced Caenorhabditis elegans (C. elegans) strains to carry lethal mutations and identifies candidate causal mutations and candidate essential genes using a robust bioinformatics procedure.

Evolutionarily conserved regulation of immunity by the splicing factor RNP-6/PUF60

Splicing is a vital cellular process that modulates important aspects of animal physiology, yet roles in regulating innate immunity are relatively unexplored. From genetic screens in C. elegans, we identified splicing factor RNP-6/PUF60 whose activity suppresses immunity, but promotes longevity, suggesting a tradeoff between these processes. Bacterial pathogen exposure affects gene expression and splicing in a rnp-6 dependent manner, and rnp-6 gain and loss-of-function activities reveal an active role in immune regulation. Another longevity promoting splicing factor, SFA-1, similarly exerts an immuno-suppressive effect, working downstream or parallel to RNP-6. RNP-6 acts through TIR-1/PMK-1/MAPK signaling to modulate immunity. The mammalian homolog, PUF60, also displays anti-inflammatory properties, and its levels swiftly decrease after bacterial infection in mammalian cells, implying a role in the host response. Altogether our findings demonstrate an evolutionarily conserved modulation of immunity by specific components of the splicing machinery.

The Enigmatic Canal-Associated Neurons Regulate Caenorhabditis elegans Larval Development Through a cAMP Signaling Pathway

Caenorhabditis elegans larval development requires the function of the two Canal-Associated Neurons (CANs): killing the CANs by laser microsurgery or disrupting their development by mutating the gene ceh-10 results in early larval arrest. How these cells promote larval development, however, remains a mystery. In screens for mutations that bypass CAN function, we identified the gene kin-29, which encodes a member of the Salt-Inducible Kinase (SIK) family and a component of a conserved pathway that regulates various C. elegans phenotypes. Like kin-29 loss, gain-of-function mutations in genes that may act upstream of kin-29 or growth in cyclic-AMP analogs bypassed ceh-10 larval arrest, suggesting that a conserved adenylyl cyclase/PKA pathway inhibits KIN-29 to promote larval development, and that loss of CAN function results in dysregulation of KIN-29 and larval arrest. The adenylyl cyclase ACY-2 mediates CAN-dependent larval development: acy-2 mutant larvae arrested development with a similar phenotype to ceh-10 mutants, and the arrest phenotype was suppressed by mutations in kin-29 ACY-2 is expressed predominantly in the CANs, and we provide evidence that the acy-2 functions in the CANs to promote larval development. By contrast, cell-specific expression experiments suggest that kin-29 acts in both the hypodermis and neurons, but not in the CANs. Based on our findings, we propose two models for how ACY-2 activity in the CANs regulates KIN-29 in target cells.

Genomic identification and functional analysis of essential genes in Caenorhabditis elegans

Background

Essential genes are required for an organism’s viability and their functions can vary greatly, spreading across many pathways. Due to the importance of essential genes, large scale efforts have been undertaken to identify the complete set of essential genes and to understand their function. Studies of genome architecture and organization have found that genes are not randomly disturbed in the genome.

Results

Using combined genetic mapping, Illumina sequencing, and bioinformatics analyses, we successfully identified 44 essential genes with 130 lethal mutations in genomic regions of C. elegans of around 7.3 Mb from Chromosome I (left). Of the 44 essential genes, six of which were genes not characterized previously by mutant alleles, let-633/let-638 (B0261.1), let-128 (C53H9.2), let-511 (W09C3.4), let-162 (Y47G6A.18), let-510 (Y47G6A.19), and let-131 (Y71G12B.6). Examine essential genes with Hi-C data shows that essential genes tend to cluster within TAD units rather near TAD boundaries. We have also shown that essential genes in the left half of chromosome I in C. elegans function in enzyme and nucleic acid binding activities during fundamental processes, such as DNA replication, transcription, and translation. From protein-protein interaction networks, essential genes exhibit more protein connectivity than non-essential genes in the genome. Also, many of the essential genes show strong expression in embryos or early larvae stages, indicating that they are important to early development.

Conclusions

Our results confirmed that this work provided a more comprehensive picture of the essential gene and their functional characterization. These genetic resources will offer important tools for further heath and disease research.

Electronic supplementary material

The online version of this article (10.1186/s12864-018-5251-3) contains supplementary material, which is available to authorized users.

Rapid and Efficient Identification of Caenorhabditis elegans Legacy Mutations Using Hawaiian SNP-Based Mapping and Whole-Genome Sequencing

The production of viable embryos requires the coordination of many cellular processes, including protein synthesis, cytoskeletal reorganization, establishment of polarity, cell migration, cell division, and in Caenorhabditis elegans, eggshell formation. Defects in any of these processes can lead to embryonic lethality. We examined six temperature-sensitive mutants as well as one nonconditional mutant that were previously identified in genetic screens as either embryonic lethal (maternal-effect or zygotic lethal) or eggshell defective. The responsible molecular lesion for each had never been determined. After confirmation of temperature sensitivity and lethality, we performed whole-genome sequencing using a single-nucleotide polymorphism mapping strategy to pinpoint the molecular lesions. Gene candidates were confirmed by RNA interference phenocopy and/or complementation tests and one mutant was further validated by CRISPR (Clustered Regularly Interspaced Short Palidromic Repeats)/Cas9 gene editing. This approach identified new alleles of several genes that had only been previously studied by RNA interference depletion. Our identification of temperature-sensitive alleles for all of these essential genes provides an extremely useful tool for further investigation for the C. elegans community, such as the ability to address mutant phenotypes at various developmental stages and the ability to carry out suppressor/enhancer screens to identify other genes that function in a specific cellular process.

High-throughput capturing and characterization of mutations in essential genes of Caenorhabditis elegans

BACKGROUND

Essential genes are critical for the development of all organisms and are associated with many human diseases. These genes have been a difficult category to study prior to the availability of balanced lethal strains. Despite the power of targeted mutagenesis, there are limitations in identifying mutations in essential genes. In this paper, we describe the identification of coding regions for essential genes mutated using forward genetic screens in Caenorhabditis elegans. The lethal mutations described here were isolated and maintained by a wild-type allele on a rescuing duplication.

RESULTS

We applied whole genome sequencing to identify the causative molecular lesion resulting in lethality in existing C. elegans mutant strains. These strains are balanced and can be easily maintained for subsequent characterization. Our method can be effectively used to analyze mutations in a large number of essential genes. We describe here the identification of 64 essential genes in a region of chromosome I covered by the duplication sDp2. Of these, 42 are nonsense mutations, six are splice signal mutations, one deletion, and 15 are non-synonymous mutations. Many of the essential genes in this region function in cell cycle, transcriptional regulation, and RNA processing.

CONCLUSIONS

The essential genes identified here are represented by mutant strains, many of which have more than one mutant allele. The genetic resource can be utilized to further our understanding of essential gene function and will be applicable to the study of C. elegans development, conserved cellular function, and ultimately lead to improved human health.

Caenorhabditis elegans SWI/SNF subunits control sequential developmental stages in the somatic gonad

The Caenorhabditis elegans somatic gonadal precursors (SGPs) are multipotent progenitors that give rise to all somatic tissues of the adult reproductive system. The hunchback and Ikaros-like gene ehn-3 is expressed specifically in SGPs and is required for their development into differentiated tissues of the somatic gonad. To find novel genes involved in SGP development, we used a weak allele of ehn-3 as the basis for a reverse genetic screen. Feeding RNAi was used to screen ∼2400 clones consisting of transcription factors, signaling components, and chromatin factors. The screen identified five members of the C. elegans SWI/SNF chromatin remodeling complex as genetic enhancers of ehn-3. We characterized alleles of 10 SWI/SNF genes and found that SWI/SNF subunits are required for viability and gonadogenesis. Two conserved SWI/SNF complexes, PBAF and BAF, are defined by their unique array of accessory subunits around a common enzymatic core that includes a catalytic Swi2/Snf2-type ATPase. Tissue-specific RNAi experiments suggest that C. elegans PBAF and BAF complexes control different processes during somatic gonadal development: PBRM-1, a signature subunit of PBAF, is important for normal SGP development, whereas LET-526, the distinguishing subunit of BAF, is required for development of a differentiated cell type, the distal tip cell (DTC). We found that the SWSN-4 ATPase subunit is required for SGP and DTC development. Finally, we provide evidence that C. elegans PBAF subunits and hnd-1/dHand are important for the cell fate decision between SGPs and their differentiated sisters, the head mesodermal cells.

Non-random distribution of extensive chromosome rearrangements in Brassica napus depends on genome organization

Chromosome rearrangements are common, but their dynamics over time, mechanisms of occurrence and the genomic features that shape their distribution and rate are still poorly understood. We used allohaploid Brassica napus (AC, n = 19) as a model to analyze the effect of genomic features on the formation and diversity of meiotically driven chromosome rearrangements. We showed that allohaploid B. napus meiosis leads to extensive new structural diversity. Almost every allohaploid offspring carried a unique combination of multiple rearrangements throughout the genome, and was thus structurally differentiated from both its haploid parent and its sister plants. This large amount of genome reshuffling was remarkably well-tolerated in the heterozygous state, as neither male nor female fertility were strongly reduced, and meiosis behavior was normal in most cases. We also used a quantitative statistical model, which accounted for 75% of the observed variation in rearrangement rates, to show that the distribution of meiotically driven chromosome rearrangements was not random but was shaped by three principal genomic features. In descending order of importance, the rate of marker loss increased strongly with genetic distance from the centromere, the degree of collinearity between chromosomes, and the genome of origin (A < C). Overall, our results demonstrate that B. napus accumulates a large number of genetic changes, but these rearrangements are not randomly distributed in the genome. The structural genetic diversity produced by the allohaploid pathway and its role in the evolution of polyploid species compared to diploid meiosis are discussed.

The near demise and subsequent revival of classical genetics for investigating Caenorhabditis elegans embryogenesis: RNAi meets next-generation DNA sequencing

Molecular genetic investigation of the early Caenorhabditis elegans embryo has contributed substantially to the discovery and general understanding of the genes, pathways, and mechanisms that regulate and execute developmental and cell biological processes. Initially, worm geneticists relied exclusively on a classical genetics approach, isolating mutants with interesting phenotypes after mutagenesis and then determining the identity of the affected genes. Subsequently, the discovery of RNA interference (RNAi) led to a much greater reliance on a reverse genetics approach: reducing the function of known genes with RNAi and then observing the phenotypic consequences. Now the advent of next-generation DNA sequencing technologies and the ensuing ease and affordability of whole-genome sequencing are reviving the use of classical genetics to investigate early C. elegans embryogenesis.

Allelic ratios and the mutational landscape reveal biologically significant heterozygous SNVs

The issue of heterozygosity continues to be a challenge in the analysis of genome sequences. In this article, we describe the use of allele ratios to distinguish biologically significant single-nucleotide variants from background noise. An application of this approach is the identification of lethal mutations in Caenorhabditis elegans essential genes, which must be maintained by the presence of a wild-type allele on a balancer. The h448 allele of let-504 is rescued by the duplication balancer sDp2. We readily identified the extent of the duplication when the percentage of read support for the lesion was between 70 and 80%. Examination of the EMS-induced changes throughout the genome revealed that these mutations exist in contiguous blocks. During early embryonic division in self-fertilizing C. elegans, alkylated guanines pair with thymines. As a result, EMS-induced changes become fixed as either G→A or C→T changes along the length of the chromosome. Thus, examination of the distribution of EMS-induced changes revealed the mutational and recombinational history of the chromosome, even generations later. We identified the mutational change responsible for the h448 mutation and sequenced PCR products for an additional four alleles, correlating let-504 with the DNA-coding region for an ortholog of a NFκB-activating protein, NKAP. Our results confirm that whole-genome sequencing is an efficient and inexpensive way of identifying nucleotide alterations responsible for lethal phenotypes and can be applied on a large scale to identify the molecular basis of essential genes.

Multiple functions of PBRM-1/Polybromo- and LET-526/Osa-containing chromatin remodeling complexes in C. elegans development

The SWI/SNF-like chromatin remodeling complexes consist of two evolutionarily conserved subclasses, which are characterized by specific accessory components, the OSA/BAF250 and Polybromo proteins. These complexes regulate the expressions of distinct sets of target genes, with some overlap, and the regulatory components are thought to determine the target specificity for each complex. Here we isolated C. elegans mutants of the genes for the OSA/BAF250 homolog, LET-526, and the Polybromo homolog, PBRM-1, in a screen for the abnormal asymmetric cell division phenotype. In the asymmetric division of the T cell, both LET-526 and PBRM-1 regulated the asymmetric expression of psa-3/Meis between the T cell daughters, suggesting that the two subclasses share the same target. In the gonad, PBRM-1 regulated gonad primordium formation during embryogenesis, whereas LET-526 was required post-embryonically for distal tip cell (DTC) production from the gonad primordium, suggesting that these proteins have distinct targets for DTC development. Thus, the same cellular process is regulated by LET-526 and PBRM-1 in the asymmetric division of the T cell, but they regulate distinct cellular processes in the gonad morphogenesis. Although disruption of the core component PSA-1 or PSA-4 caused similar defects in the gonad and T cell, it also caused early embryonic arrest, which was not observed in the let-526, pbrm-1, or let-526 pbrm-1 double mutants, suggesting that some targets of SWI/SNF-like complexes do not require LET-526 or PBRM-1 for their transcription. Our results show that the target selection by SWI/SNF-like complexes during C. elegans development is intricately regulated by accessory components.

Rapid mapping and identification of mutations in Caenorhabditis elegans by restriction site-associated DNA mapping and genomic interval pull-down sequencing

Forward genetic screens provide a powerful approach for inferring gene function on the basis of the phenotypes associated with mutated genes. However, determining the causal mutation by traditional mapping and candidate gene sequencing is often the rate-limiting step, especially when analyzing many mutants. We report two genomic approaches for more rapidly determining the identity of the affected genes in Caenorhabditis elegans mutants. First, we report our use of restriction site-associated DNA (RAD) polymorphism markers for rapidly mapping mutations after chemical mutagenesis and mutant isolation. Second, we describe our use of genomic interval pull-down sequencing (GIPS) to selectively capture and sequence megabase-sized portions of a mutant genome. Together, these two methods provide a rapid and cost-effective approach for positional cloning of C. elegans mutant loci, and are also applicable to other genetic model systems.

Genomic sequence of a mutant strain of Caenorhabditis elegans with an altered recombination pattern

Background

The original sequencing and annotation of the Caenorhabditis elegans genome along with recent advances in sequencing technology provide an exceptional opportunity for the genomic analysis of wild-type and mutant strains. Using the Illumina Genome Analyzer, we sequenced the entire genome of Rec-1, a strain that alters the distribution of meiotic crossovers without changing the overall frequency. Rec-1 was derived from ethylmethane sulfonate (EMS)-treated strains, one of which had a high level of transposable element mobility. Sequencing of this strain provides an opportunity to examine the consequences on the genome of altering the distribution of meiotic recombination events.

Results

Using Illumina sequencing and MAQ software, 83% of the base pair sequence reads were aligned to the reference genome available at Wormbase, providing a 21-fold coverage of the genome. Using the software programs MAQ and Slider, we observed 1124 base pair differences between Rec-1 and the reference genome in Wormbase (WS190), and 441 between the mutagenized Rec-1 (BC313) and the wild-type N2 strain (VC2010). The most frequent base-substitution was G:C to A:T, 141 for the entire genome most of which were on chromosomes I or X, 55 and 31 respectively. With this data removed, no obvious pattern in the distribution of the base differences along the chromosomes was apparent. No major chromosomal rearrangements were observed, but additional insertions of transposable elements were detected. There are 11 extra copies of Tc1, and 8 of Tc2 in the Rec-1 genome, most likely the remains of past high-hopper activity in a progenitor strain.

Conclusion

Our analysis of high-throughput sequencing was able to detect regions of direct repeat sequences, deletions, insertions of transposable elements, and base pair differences. A subset of sequence alterations affecting coding regions were confirmed by an independent approach using oligo array comparative genome hybridization. The major phenotype of the Rec-1 strain is an alteration in the preferred position of the meiotic recombination event with no other significant phenotypic consequences. In this study, we observed no evidence of a mutator effect at the nucleotide level attributable to the Rec-1 mutation.

Evolution of the Caenorhabditis elegans genome

A fundamental problem in genome biology is to elucidate the evolutionary forces responsible for generating nonrandom patterns of genome organization. As the first metazoan to benefit from full-genome sequencing, Caenorhabditis elegans has been at the forefront of research in this area. Studies of genomic patterns, and their evolutionary underpinnings, continue to be augmented by the recent push to obtain additional full-genome sequences of related Caenorhabditis taxa. In the near future, we expect to see major advances with the onset of whole-genome resequencing of multiple wild individuals of the same species. In this review, we synthesize many of the important insights to date in our understanding of genome organization and function that derive from the evolutionary principles made explicit by theoretical population genetics and molecular evolution and highlight fertile areas for future research on unanswered questions in C. elegans genome evolution. We call attention to the need for C. elegans researchers to generate and critically assess nonadaptive hypotheses for genomic and developmental patterns, in addition to adaptive scenarios. We also emphasize the potential importance of evolution in the gonochoristic (female and male) ancestors of the androdioecious (hermaphrodite and male) C. elegans as the source for many of its genomic and developmental patterns.

Retroviral insertions in the VISION database identify molecular pathways in mouse lymphoid leukemia and lymphoma

AKXD recombinant inbred (RI) strains develop a variety of leukemias and lymphomas due to somatically acquired insertions of retroviral DNA into the genome of hematopoetic cells that can mutate cellular proto-oncogenes and tumor suppressor genes. We generated a new set of tumors from nine AKXD RI strains selected for their propensity to develop B-cell tumors, the most common type of human hematopoietic cancers. We employed a PCR technique called viral insertion site amplification (VISA) to rapidly isolate genomic sequence at the site of provirus insertion. Here we describe 550 VISA sequence tags (VSTs) that identify 74 common insertion sites (CISs), of which 21 have not been identified previously. Several suspected proto-oncogenes and tumor suppressor genes lie near CISs, providing supportive evidence for their roles in cancer. Furthermore, numerous previously uncharacterized genes lie near CISs, providing a pool of candidate disease genes for future research. Pathway analysis of candidate genes identified several signaling pathways as common and powerful routes to blood cancer, including Notch, E-protein, NFκB, and Ras signaling. Misregulation of several Notch signaling genes was confirmed by quantitative RT-PCR. Our data suggest that analyses of insertional mutagenesis on a single genetic background are biased toward the identification of cooperating mutations. This tumor collection represents the most comprehensive study of the genetics of B-cell leukemia and lymphoma development in mice. We have deposited the VST sequences, CISs in a genome viewer, histopathology, and molecular tumor typing data in a public web database called VISION (Viral Insertion Sites Identifying Oncogenes), which is located at http://www.mouse-genome.bcm.tmc.edu/vision.

Efficient target-selected mutagenesis in Caenorhabditis elegans: toward a knockout for every gene

Reverse genetic or gene-driven knockout approaches have contributed significantly to the success of model organisms for fundamental and biomedical research. Although various technologies are available for C. elegans, none of them scale very well for genome-wide application. To address this, we implemented a target-selected knockout approach that is based on random chemical mutagenesis and detection of single nucleotide mutations in genes of interest using high-throughput resequencing. A clonal library of 6144 EMS-mutagenized worms was established and screened, resulting in the identification of 1044 induced mutations in 109 Mbp, which translates into an average spacing between exonic mutations in the library of only 17 bp. We covered 25% of the open reading frames of 32 genes and identified one or more inactivating mutations (nonsense or splice site) in 84% of them. Extrapolation of our results indicates that nonsense mutations for >90% of all C. elegans genes are present in the library. To identify all of these mutations, one only needs to inspect those positions that--given the known specificity of the mutagen--can result in the introduction of a stop codon. We define these positions as nonsense introducing mutations (NIMs). The genome-wide collection of possible NIMs can be calculated for any organism with a sequenced genome and reduces the screening complexity by 200- to 2000-fold, depending on the organism and mutagen. For EMS-mutagenized C. elegans, there are only approximately 500,000 NIMs. We show that a NIM genotyping approach employing high-density microarrays can, in principle, be used for the genome-wide identification of C. elegans knockouts.

Regional variation in the density of essential genes in mice

In most species, and particularly in vertebrates, the percentage of genes absolutely required for survival, the essential genes, has not been estimated. To obtain this estimation, we used the mouse as an experimental model to carry out high-efficiency N-ethyl-N-nitrosourea (ENU) mutagenesis screens in two balancer chromosome regions, and compared our results to a third previously published screen. The number of essential genes in each region was predicted based on allele frequencies. We determined that the density of essential genes differs by up to an order of magnitude among genomic regions. This indicates that extrapolating from regional estimates to genome-wide estimates of essential genes has a huge variance. A particularly high density of essential genes on mouse Chromosome 11 coincides with a high degree of regional linkage conservation, providing a possible causal explanation for the density variation. This is the first demonstration of regional variation in essential gene density in the mouse genome.

Deletion polymorphism in wheat chromosome regions with contrasting recombination rates

Polymorphism for deletions was investigated in 1027 lines of tetraploid and hexaploid wheat and 420 lines of wheat diploid ancestors. A total of 26 deletions originating during the evolution of polyploid wheat were discovered among 155 investigated loci. Wheat chromosomes were divided into a proximal, low-recombination interval containing 69 loci and a distal, high-recombination interval containing 86 loci. A total of 23 deletions involved loci in the distal, high-recombination interval and only 3 involved loci in the proximal, low-recombination interval. The rates of DNA loss differed by several orders of magnitude in the two intervals. The rate of diploidization of polyploid wheat by deletions was estimated and was shown to have proceeded faster in the distal, high-recombination interval than in the proximal, low-recombination interval.

Bioinformatical assay of human gene morbidity

Only a fraction of eukaryotic genes affect the phenotype drastically. We compared 18 parameters in 1273 human morbid genes, known to cause diseases, and in the remaining 16 580 unambiguous human genes. Morbid genes evolve more slowly, have wider phylogenetic distributions, are more similar to essential genes of Drosophila melanogaster, code for longer proteins containing more alanine and glycine and less histidine, lysine and methionine, possess larger numbers of longer introns with more accurate splicing signals and have higher and broader expressions. These differences make it possible to classify as non-morbid 34% of human genes with unknown morbidity, when only 5% of known morbid genes are incorrectly classified as non-morbid. This classification can help to identify disease-causing genes among multiple candidates.

The organization and rate of evolution of wheat genomes are correlated with recombination rates along chromosome arms

Genes detected by wheat expressed sequence tags (ESTs) were mapped into chromosome bins delineated by breakpoints of 159 overlapping deletions. These data were used to assess the organizational and evolutionary aspects of wheat genomes. Relative gene density and recombination rate increased with the relative distance of a bin from the centromere. Single-gene loci present once in the wheat genomes were found predominantly in the proximal, low-recombination regions, while multigene loci tended to be more frequent in distal, high-recombination regions. One-quarter of all gene motifs within wheat genomes were represented by two or more duplicated loci (paralogous sets). For 40 such sets, ancestral loci and loci derived from them by duplication were identified. Loci derived by duplication were most frequently located in distal, high-recombination chromosome regions whereas ancestral loci were most frequently located proximal to them. It is suggested that recombination has played a central role in the evolution of wheat genome structure and that gradients of recombination rates along chromosome arms promote more rapid rates of genome evolution in distal, high-recombination regions than in proximal, low-recombination regions.

Large-scale analysis of gene function in Caenorhabditis elegans by high-throughput RNAi

Genome-wide analysis of gene function is essential for the post-genome era, and development of efficient and economical technology suitable for it has been in demand. Here we report a large-scale inactivation of the expressed genes in the nematode Caenorhabditis elegans. For this purpose, we have established a high-throughput "RNAi-by-soaking" methodology by modifying the conventional RNAi method [1, 2]. A set of tag-sequenced, nonredundant cDNAs corresponding to approximately 10,000 genes [3] (representing half of the predicted genes [4]) was used for the systematic RNAi analysis. We have processed approximately 2500 genes to date. In development, 27% of them showed detectable phenotypes, such as embryonic lethality, post-embryonic lethality, sterility, and morphological abnormality. Of these, we analyzed the phenotypes of F1 sterility in detail, and we have identified 24 genes that might play important roles in germline development. Combined with the ongoing analysis of expression patterns of these cDNAs [3, 5], the functional information obtained in this work will provide a starting point for the further analysis of each gene. Another finding from this screening is that the incidence of essential genes is significantly lower in the X chromosome than in the autosomes.