Oligonucleotide Array Comparative Genomic Hybridization (oaCGH) based characterization of genetic deficiencies as an aid to gene mapping in Caenorhabditis elegans
© Jones et al; licensee BioMed Central Ltd. 2007
Received: 11 August 2007
Accepted: 07 November 2007
Published: 07 November 2007
A collection of genetic deficiencies covering over 70% of the Caenorhabditis elegans genome exists, however the application of these valuable biological tools has been limited due to the incomplete correlation between their genetic and physical characterization.
We have applied oligonucleotide array Comparative Genomic Hybridization (oaCGH) to the high resolution, molecular characterization of several genetic deficiency and duplication strains in a 5 Mb region of Chromosome III. We incorporate this data into a physical deficiency map which is subsequently used to direct the positional cloning of essential genes within the region. From this analysis we are able to quickly determine the molecular identity of several previously unidentified mutations.
We have applied accurate, high resolution molecular analysis to the characterization of genetic mapping tools in Caenorhabditis elegans. Consequently we have generated a valuable physical mapping resource, which we have demonstrated can aid in the rapid molecular identification of mutations of interest.
A large resource of deletion strains (also known as genetic deficiencies) accounting for over 70% of the Caenorhabditis elegans genome has been generated by various research groups over the past three decades . These genetic deficiencies have proven advantageous for a variety of purposes including; characterization of mutant alleles , identification of specific loci affecting developmental processes , investigation of genome replication and stability [4, 5] and, most significantly, as tools for positional cloning of unmapped mutations to discrete regions of the genome [1, 6–8].
The full potential of these biological tools has however been limited due to the lack of high resolution characterization at a genome wide scale. The mapping of physical breakpoint positions within each deficiency strain is required to allow the deleted gene complement for that strain to be precisely defined. Additionally, genetic deficiencies may exhibit molecular complexity preventing their reliable use in mapping experiments .
Previously, characterization of genetic deficiencies has been performed by fairly low resolution or labor intensive techniques such as genetic linkage mapping, PCR analysis  and, more recently, by the application of snip-SNP [9, 10] and as a consequence many available deficiency strains remain poorly characterized.
Oligonucleotide array Comparative Genomic Analysis (oaCGH) is an emerging technology for high resolution mapping of chromosomal copy number changes at a genome wide scale through the comparison of the DNA ratio between two samples from the same organism [11, 12]. The recent development of a C. elegans specific oaCGH platform for identification of novel single gene deletions  represents a powerful technology that can be adapted to the rapid and precise characterization of deficiency mapping strains.
In this study we demonstrate the successful application of oaCGH to the physical characterization of deficiency strains in C. elegans. We use this data to annotate a physical deficiency map within a 5 Mb region of chromosome III and demonstrate the application of this map to aid in the molecular identification of previously generated mutations known to reside within this region.
Results and Discussion
oaCGH mapped deletions and duplications physically define 17 zones around the dpy-17 region of Chromosome III
7 deficiencies and 2 duplications lying in the region of dpy-17 on chromosome III were chosen for oaCGH analysis (Nimblegen) as they have been previously characterized by both genetic linkage and PCR analysis, and used to roughly position a large number of unidentified EMS generated lethal mutants . After oaCGH mapping had precisely defined the gene complement for each of these deficiencies a refined candidate gene approach was implemented to rapidly identify mutations in essential genes which map to this region.
Summary of mapping resolution for deficiencies using both genetic and molecular mapping techniques.
Genetic/PCR Mapping Data
oaCGH Mapping Data
a Sequenced Breakpoint Coordinates
Size of deficiency (kb)
Average resolution (kb)
Physical and genetic structure of zones 1A-14.
b # RNAi Lethal
# Lethal Mutants
7 c (+/- 8)
14 c (+/- 8)
The free translocation sDp3 undergoes significant modification within deficiency strains
This initial analysis revealed that the extent of sDp3 varies dramatically between each of the strains tested (Figure 2B–D). Since duplications maintained in balancer strains are known to breakdown both spontaneously and under the mutagenic conditions used to induce deficiencies [14–16] it is likely that the observed modifications to sDp3 are as a result of these factors.
In addition to the variability of sDp3 coverage, discrete duplications internal to sDp3 are observed within each of the strains tested (Figure 2). This lack of uniformity makes determination of deficiency structure problematic in some cases. For example, the presence of an internal duplication in BC4697 across the region containing the left breakpoint of sDf121 results in an increase of predicted copy number in this region (Figure 2F), while multiple variations within both sDp3 balanced strains containing sDf127 and sDf128 complicates the distinction between where DNA has been deleted or duplicated (Figure 2C–E and data not shown).
oaCGH of heterozygous deficiencies resolves copy number ambiguity and results in precise mapping of deletion breakpoints
Since heterozygous deletions can be easily resolved with oaCGH  we surmised that the creation of single copy deficiency strains would eliminate the ambiguity introduced through the use of the sDp3 balancer. By using this approach the expected ratio of reference DNA to test is 1:1 along the whole length of the chromosome except within the deleted region where a ratio of 2:1 will be expected (Figure 2B).
To implement this strategy all of the sDp3 balanced deficiency strains to be tested were out-crossed to JK2689, which contains the GFP marked heterozygous balancer hT2 (qIs48) known to cover the region under investigation in this study . For this analysis genomic JK2689 DNA was used as the reference sample to eliminate any genomic variability that may have been introduced into this strain through the heavy mutagenesis used for its original construction .
Combining data from multiple array experiments allows for bias removal resulting in reliable deficiency characterization
oaCGH data obtained from separate strains using the same chip design performs consistently enough to allow for the integration and direct comparison of the data generated from multiple experiments. Through this data integration, common inconsistencies that are not specific to the deletion of interest are identified and eliminated from the annotation, resulting in a reliable characterization of the whole deficiency genome (Figure 3B and 3C). This is an important step given that small mutagenic events which may have occurred outside the known boundaries of a given deficiency may remain unidentified by traditional mapping methods and could lead to conflicting mapping data.
oaCGH analysis reveals deficiency complexity on a genome-wide scale
We have analysed a number of deficiencies generated with both GRI and UV irradiation (Table 1 and data not shown), and have been unable to detect any significant mutagen specific deletion characteristics, the majority having a simple contiguous structure (Figure 3B and 3C, and data not shown). Exceptions to this observation are seen in the deficiencies sDf130, which contains two distinct deletions, the smaller of which defines zone 10 (Figure 1), and sDf128 which exhibits associated complexity.
Accurate identification of lethal mutations positioned within the deficiency map
Since 39 of the 40 molecularly identified lethal mutants in Wormbase  exhibit a variety of lethal or sterile RNAi phenotypes (Additional file 1), genes within each zone which display these RNAi phenotypes make good candidates for the lethal mutations mapping to the same region. Of the 21 genes within zone 1 B, only 6 exhibit lethal phenotypes by RNAi (Figure 5 and Additional file 1) and these 6 genes were consequently taken as initial candidates of the three lethal mutants in this region.
Summary of the mutations identified in essential genes as a result of this study.
Genomic Position (bp)
RRP5 Ribosomal RNA Processing Enzyme
3rd Intron Splice Site
Finally, through the incorporation of mapping data generated by our group from ongoing efforts to clone lethal mutations, we have been able to further identify the mutations let-716 and let-768 as alleles of the genes C16A3.3 and fum-1, respectively (our unpublished data and summarized in Table 3).
We have demonstrated that oaCGH can be successfully applied to the rapid and precise characterization of existing C. elegans genetic deficiencies. This study highlights how this type of analysis can transform low resolution genetic mapping tools into a precise physical mapping resource with which to accurately position molecularly unidentified mutations.
The implementation of oaCGH technology for this purpose is straightforward, requiring only the preparation of high quality genomic DNA from the deficiency strains of interest. Data is generated in a short time and can be verified experimentally by standard molecular techniques. Some caveats to the application of this methodology have been highlighted by the initial difficulties we describe resulting from the use of animals balanced by a free duplication, though we have demonstrated that by considered evaluation of the experimental approach these issues can be circumvented easily.
We envisage that this resource will be instrumental in resolving previously ambiguous genetic mapping data. Additionally, we propose that integration of a physical deficiency map into a high-throughput molecular mapping strategy such as snip-SNP will improve the efficiency of positional cloning in C. elegans. In such an approach snip-SNP would be employed to map the mutant of interest into a sub-chromosomal region, this region could then be further broken up by physically defined deficiencies, and a candidate sequencing pipeline, such as the one we describe here, employed to identify the molecular position of the mutation of interest.
Nematode culture, harvest, and DNA preparation
Nematodes were grown as previously described  on 150 mm NGM agar plates seeded with Escherichia coli strain OP50. Strains used: Bristol N2, JK2689 – pop-1(q645) dpy-5(e61) I/hT2 [bli-4(e937) let-?(q782) qIs48] (I;III), BC4697 – sDf121(s2098) unc-32(e189) III; sDp3 (III;f), BC7162 – dpy-17(e164) sDf125(s2424) unc-32(e189) III; hT2 [bli-4(e937) let-?(q782) qIs48], BC4638 – dpy-17(e164) sDf127(s2428) unc-32(e189) III; sDp3 (III;f), BC7163 – dpy-17(e164) sDf127(s2428) unc-32(e189) III; hT2 [bli-4(e937) let-?(q782) qIs48], BC4690 – dpy-17(e164) sDf128(s2786) unc-32(e189) III; sDp3 (III;f), BC7164 – dpy-17(e164) sDf128(s2786) unc-32(e189) III; hT2 [bli-4(e937) let-?(q782) qIs48], BC7244 – sDf130(s2427) unc-32(e189) III; hT2 [bli-4(e937) let-?(q782) qIs48], BC7181 – dpy-17(e164) sDf135(s2767) unc-32(e189) III; hT2 [bli-4(e937) let-?(q782) qIs48], BC3737 – sDp8(III;I); eT1(III;V), BC4158 sDp3(III;f); dpy-17(e164) let-722(s2448) unc-32(e189) III, BC4160 – sDp3(III;f); dpy-17(e164) let-722(s2450) unc-32(e189) III, BC4228 – sDp3(III;f); dpy-17(e164) let-722(s2588) unc-32(e189), BC4262 – sDp3(III;f); dpy-17(e164) let-754(s2622) unc-32(e189) III, BC4173 – dpy-17(e164) let-766(s2463) unc-32(e189) III; sDp3(III;f), BC4167 – sDp3(III;f); dpy-17(e164) let-716(s2457) unc-32(e189) III, BC4266 – sDp3(III;f); dpy-17(e164) let-716(s2626) unc-32(e189) III, BC4192 – sDp3(III;f); dpy-17(e164) let-768(s2482) unc-32(e189) III, BC4232 – sDp3(III;f); dpy-17(e164) let-768(s2592) unc-32(e189) III, BC4268 – sDp3(III;f); dpy-17(e164) let-768(s2628) unc-32(e189) III, BC2842 – dpy-18(e364)/eT1 III; unc-46(e177) let-444(s1418)/eT1 [let-500(s2165)]V, MT3022 – nDf20/sma-2(e502) unc-32(e189)III.
DNA preparation for oaCGH samples was done as previously described . For mutant sequencing ~50 homozygous lethal animals were picked from the progeny of sDp3 balanced strains into 5 ul of lysis buffer (50 mM KCl, 10 mM Tris pH 8.3, 2.5 mM MgCl2, 0.45% Tween-20, 0.01% Gelatin) and freeze cracked in Liquid N2 followed by digestion (1 hour at 60°C followed by 15 min at 95°C). 1 ul of this DNA preparation was subsequently used in a 25 ul PCR reaction.
oaCGH Data analysis
For the analysis of sDp3 balanced strains a C. elegans whole genome array was used, based on Wormbase release WS120, and available from NimbleGen Systems Inc. . Analysis of the hT2 (qIs48) balanced strains was performed with a whole genome C. elegans array designed with overlapping 50-mer probes targeting primarily annotated exons and micro-RNAs . oaCGH sample preparation, hybridization and analysis was done as previously described . Copy number aberrations were detected by visual inspection using the SignalMap™ browser software . The data discussed in this publication have been deposited in NCBIs Gene Expression Omnibus (GEO)  and are accessible through GEO Series accession numbers GSE9214 and GPL6047.
Molecular identification of deficiency breakpoints and mutations
PCR amplification across the region of the breakpoint was performed with nested primers (Primer sequences and amplification conditions available upon request) and sequenced using standard molecular methods.
We thank Dr Nigel O'Neil for his enthusiasm and advice concerning this project and his comments on the manuscript. This work was supported by grants from the Canadian Space Agency, the Canadian Institute of Health Research and from Genome Canada to DGM.
- Girard LR, Fiedler TJ, Harris TW, Carvalho F, Antoshechkin I, Han M, Sternberg PW, Stein LD, Chalfie M: WormBook: the online review of Caenorhabditis elegans biology. Nucleic Acids Res. 2007, 35 (Database issue): D472-5. 10.1093/nar/gkl894.PubMed CentralPubMedView ArticleGoogle Scholar
- Sigurdson DC, Spanier GJ, Herman RK: Caenorhabditis elegans deficiency mapping. Genetics. 1984, 108 (2): 331-345.PubMed CentralPubMedGoogle Scholar
- Chanal P, Labouesse M: A screen for genetic loci required for hypodermal cell and glial-like cell development during Caenorhabditis elegans embryogenesis. Genetics. 1997, 146 (1): 207-226.PubMed CentralPubMedGoogle Scholar
- Rosenbluth RE, Johnsen RC, Baillie DL: Pairing for recombination in LGV of Caenorhabditis elegans: a model based on recombination in deficiency heterozygotes. Genetics. 1990, 124 (3): 615-625.PubMed CentralPubMedGoogle Scholar
- McKim KS, Peters K, Rose AM: Two types of sites required for meiotic chromosome pairing in Caenorhabditis elegans. Genetics. 1993, 134 (3): 749-768.PubMed CentralPubMedGoogle Scholar
- Johnsen RC, Baillie DL: Genetic analysis of a major segment [LGV(left)] of the genome of Caenorhabditis elegans. Genetics. 1991, 129 (3): 735-752.PubMed CentralPubMedGoogle Scholar
- Stewart HI, O'Neil NJ, Janke DL, Franz NW, Chamberlin HM, Howell AM, Gilchrist EJ, Ha TT, Kuervers LM, Vatcher GP, Danielson JL, Baillie DL: Lethal mutations defining 112 complementation groups in a 4.5 Mb sequenced region of Caenorhabditis elegans chromosome III. Mol Gen Genet. 1998, 260 (2-3): 280-288. 10.1007/PL00013816.PubMedView ArticleGoogle Scholar
- Rogalski TM, Moerman DG, Baillie DL: Essential genes and deficiencies in the unc-22 IV region of Caenorhabditis elegans. Genetics. 1982, 102 (4): 725-736.PubMed CentralPubMedGoogle Scholar
- Kadandale P, Geldziler B, Hoffmann M, Singson A: Use of SNPs to determine the breakpoints of complex deficiencies, facilitating gene mapping in Caenorhabditis elegans. BMC Genet. 2005, 6 (1): 28-10.1186/1471-2156-6-28.PubMed CentralPubMedView ArticleGoogle Scholar
- Zhao Y, Lai K, Cheung I, Youds J, Tarailo M, Tarailo S, Rose A: A mutational analysis of Caenorhabditis elegans in space. Mutat Res. 2006, 601 (1-2): 19-29.PubMedView ArticleGoogle Scholar
- Dhami P, Coffey AJ, Abbs S, Vermeesch JR, Dumanski JP, Woodward KJ, Andrews RM, Langford C, Vetrie D: Exon array CGH: detection of copy-number changes at the resolution of individual exons in the human genome. Am J Hum Genet. 2005, 76 (5): 750-762. 10.1086/429588.PubMed CentralPubMedView ArticleGoogle Scholar
- Selzer RR, Richmond TA, Pofahl NJ, Green RD, Eis PS, Nair P, Brothman AR, Stallings RL: Analysis of chromosome breakpoints in neuroblastoma at sub-kilobase resolution using fine-tiling oligonucleotide array CGH. Genes Chromosomes Cancer. 2005, 44 (3): 305-319. 10.1002/gcc.20243.PubMedView ArticleGoogle Scholar
- Maydan JS, Flibotte S, Edgley ML, Lau J, Selzer RR, Richmond TA, Pofahl NJ, Thomas JH, Moerman DG: Efficient high-resolution deletion discovery in Caenorhabditis elegans by array comparative genomic hybridization. Genome Res. 2007, 17 (3): 337-347. 10.1101/gr.5690307.PubMed CentralPubMedView ArticleGoogle Scholar
- McKim KS, Rose AM: Chromosome I duplications in Caenorhabditis elegans. Genetics. 1990, 124 (1): 115-132.PubMed CentralPubMedGoogle Scholar
- McKim KS, Rose AM: Spontaneous duplication loss and breakage in Caenorhabditis elegans. Genome. 1994, 37 (4): 595-606.PubMedView ArticleGoogle Scholar
- Hedgecock EM, Herman RK: The ncl-1 gene and genetic mosaics of Caenorhabditis elegans. Genetics. 1995, 141 (3): 989-1006.PubMed CentralPubMedGoogle Scholar
- Siegfried KR, Kimble J: POP-1 controls axis formation during early gonadogenesis in C. elegans. Development. 2002, 129 (2): 443-453.PubMedGoogle Scholar
- Wormbase: release WS182 date 2007, [http://www.wormbase.org/]
- Brenner S: The genetics of Caenorhabditis elegans. Genetics. 1974, 77 (1): 71-94.PubMed CentralPubMedGoogle Scholar
- NimbleGen Systems Inc.: [http://www.nimblegen.com/]
- Gene Expression Omnibus: [http://www.ncbi.nlm.nih.gov/geo/]