The nematode C. briggsae serves as a useful model organism for comparative analysis of developmental and behavioral processes. The amenability of C. briggsae to genetic manipulations and the availability of its genome sequence have prompted researchers to study evolutionary changes in gene function and signaling pathways. These studies rely on the availability of forward genetic tools such as mutants and mapping markers.
We have computationally identified more than 30,000 polymorphisms (SNPs and indels) in C. briggsae strains AF16 and HK104. These include 1,363 SNPs that change restriction enzyme recognition sites (snip-SNPs) and 638 indels that range between 7 bp and 2 kb. We established bulk segregant and single animal-based PCR assay conditions and used these to test 107 polymorphisms. A total of 75 polymorphisms, consisting of 14 snip-SNPs and 61 indels, were experimentally confirmed with an overall success rate of 83%. The utility of polymorphisms in genetic studies was demonstrated by successful mapping of 12 mutations, including 5 that were localized to sub-chromosomal regions. Our mapping experiments have also revealed one case of a misassembled contig on chromosome 3.
We report a comprehensive set of polymorphisms in C. briggsae wild-type strains and demonstrate their use in mapping mutations. We also show that molecular markers can be useful tools to improve the C. briggsae genome sequence assembly. Our polymorphism resource promises to accelerate genetic and functional studies of C. briggsae genes.
Comparative analysis of developmental and behavioral processes in closely related species is a powerful approach to understand the mechanisms of evolution. It facilitates identification of molecular components that are conserved over millions of years due to their role in specifying common features as well as those that are variable because they confer species-specific features. The model organism Caenorhabditis elegans (a nematode) and its congener, C. briggsae, are particularly suitable for such investigations. Their many experimental advantages include rapid growth, small size, transparency, ease of culture and genetic manipulations, and the availability of fully sequenced genomes [1-3].
C. briggsae is phenotypically almost indistinguishable from C. elegans and has a similar (hermaphroditic) reproductive mode. The last common ancestor of these two species lived about 30 million years ago , and despite the rapid molecular evolution typical of the family Rhabditidae, more than half (~52%) of the C. elegans genome aligns with the C. briggsae genome assembly . This includes two-thirds of all C. briggsae genes (13,107 or 67.8%) with reciprocal orthologs in C. elegans . Thus C. elegans-C. briggsae comparative genomic and genetic studies promise powerful new tools for the identification of genes and pathways and the study of both conservation and divergence.
Like C. elegans, C. briggsae has six chromosomes that display extensive conservation of synteny, but not exact colinearity relative to C. elegans . While C. briggsae shares many of the experimental advantages of C. elegans, it has the further advantage of increased natural variability for single nucleotide polymorphisms (SNPs) and insertion-deletions (indels) [7,8]. This elevated natural variation potentially enhances its use for genotype-phenotype association studies, and is also very useful for the mapping aspects of forward genetics projects.
Initial work on gene function in C. briggsae employed cross-species transgene rescue of C. elegans mutants (e.g. [9-13]) and RNA interference (RNAi; e.g. [13-17]). However, a number of laboratories are now generating true mutations in C. briggsae, using both forward mutagenesis screens [18-20] (R.E. Ellis, personal communication; B.P. Gupta, unpublished results; H. Chamberlin, unpublished results) and PCR-based deletion mutation screens . Positional cloning of C. briggsae mutations without relying upon obvious candidate genes requires a set of mapping tools. Development of such tools is facilitated by a high-quality whole-genome shotgun assembly  and the organization of most of the resulting contigs into chromosomes via a SNP-based recombination map .
Among the tools needed to facilitate forward genetics in C. briggsae, a set of easily scored DNA polymorphisms is especially important. Experimentally validated polymorphisms can serve as useful markers for mapping mutations that cause visible phenotypes. Additionally, these markers can be integrated with the phenotype-based genetic linkage map (e.g., dpy and unc mutants ) to further enhance their utility. Integration of polymorphisms and phenotype-based maps increases map density and anchors the relative locations of molecular and phenotypic markers. With this goal in mind we have discovered a large set of genome-wide polymorphisms (SNPs and indels) in wild-type strains, using AF16 as a reference strain and four other natural isolate strains: HK104, HK105, VT847, and PB800.
The indels were placed into three classes: small (7-49 bp), medium (50-2,000 bp), and large (>2 kb). We have focused on medium and small indels (212 and 7,530, respectively), which offer the greatest utility as genetic markers. In the case of SNPs (23,829), we found that 4,700 modify restriction enzyme sites (termed snip-SNPs) and therefore can be easily detected as restriction fragment length polymorphisms (RFLPs). We established assay conditions for bulk segregant analysis (BSA) and used these to experimentally validate 14 snip-SNPs, 28 medium and 32 small indels. The validated polymorphisms were used to genetically map known mutations causing visible phenotypes thus demonstrating the effectiveness of the polymorphisms in linkage mapping studies. We also developed single animal-based PCR assay to determine map distance. Five mutations were successfully localized to sub-chromosomal regions by 3 or more indels, greatly facilitating the search for each candidate gene. These results demonstrate the utility of our mapping toolkit in genetic linkage and gene identification studies.
We performed SNP discovery in four C. briggsae strains by aligning paired shotgun sequence reads to the AF16-based reference sequence (cb25 assembly, ). These sequences were obtained by capillary gel electrophoresis at Washington University Genome Center (see Methods). To build on previous SNP discovery efforts , we applied the ssahaSNP algorithm, which detects SNPs and small indels based on SSAHA alignments to a reference sequence (see Methods). Compared to AF16, ssahaSNP detected 23,829 unique SNP loci in HK104 DNA, or one substitution per 163 bp on average (Table 1 and additional file 1). Consistent with C. briggsae clade structure , SNP density was slightly lower in strains HK105 (1/168 bp) and PB800 (1/197 bp) and much lower in strain VT847 (1/475 bp). In HK104, the most common substitution by far was A(T) to G(C), which accounted for 57.1% of all substitutions (Figure 1A).
Format: XLS Size: 4.9MB Download file
This file can be viewed with: Microsoft Excel Viewer
Table 1. SNPs in various C. briggsae strains identified by ssahaSNP (in comparison to AF16).
Figure 1. SNP discovery results for C. briggsae strain HK104. (A) Incidence of SNPs and small indels identified by ssahaSNP, by polymorphism type. (B) Distribution of small (<50 bp) indels, by size. (C) Distribution of moderate (50-2,000 bp) indels, by size.
RFLP genotyping assays and validation of snip-SNPs
We screened the SNPs predicted for HK104 for variants that altered the recognition site of a restriction enzyme, and thus might be amenable to restriction fragment length polymorphism (RFLP) genotyping. To make this a practical resource, we limited the analysis to 30 restriction enzymes from REBASE  that are reliable and inexpensive. Of 23,829 HK104 SNPs, some 4,700 (19.72%) were predicted to alter the recognition site of at least one of the 30 restriction enzymes. To develop restriction fragment length polymorphism (RFLP) assays from these snip-SNPs, we designed PCR primers with a standard protocol and performed in silico digests of the resulting amplicons to infer the banding patterns for each strain. RFLPs not easily distinguishable on a gel, or SNPs on ultracontigs not yet included in the genetic map, were removed. Finally, we used assembly AGP information and BLAST alignment to obtain coordinates for each snip-SNP on the cb3 sequence assembly. Our set contained 1,987 predicted RFLP assays from 1,362 snip-SNPs (some SNPs alter multiple RE sites) positioned on both the genetic and physical maps (see additional file 2). Another snip-SNP, bdP3, was identified in a separate study (see Methods and additional files 1 and 2).
Additional file 2. List of snip-SNPs and RFLP assays. An Excel file containing 1988 SNPs that alter restriction enzyme site (snip-SNPs) along with PCR primers and predicted digestion patterns in AF16 and HK104.
Format: XLS Size: 523KB Download file
This file can be viewed with: Microsoft Excel Viewer
We selected a total of 20 RFLP assays (between 3 and 4 for each chromosome) based on HindIII, DraI and SalI snip-SNPs for validation in AF16 and HK104 parental DNA (Table 2). Roughly a third of the assays (6) failed PCR in one or both strains in repeated attempts. Although we did not investigate the issue of PCR failure, it is possible that redesigning primers (by moving them out or in) and testing different PCR conditions may produce desired products in some cases. All of the 14 assays successfully gave rise to strain-specific RFLP banding patterns, validating the predicted snip-SNP (Table 2, Figure 2A for two examples). Interestingly, two of these assays (cb55670 and cb20723) exhibited HK104 fragments that varied from in silico predictions, another possible consequence of unknown variants in this highly divergent strain. Consistent with C. briggsae clade structure , VT847 was not polymorphic (from AF16) for the snip-SNPs we examined.
Table 2. List of validated snip-SNPs by RFLP assays.
Figure 2. Validation of polymorphism assays in C. briggsae. (A) RFLP assays for snip-SNPs in three parental strains, namely AF16 (A), HK104 (H), and VT847 (V) by HindIII restriction digestion. (B) Medium indels in three parental strains showing a 100 bp deletion on fpc4171, a 145 bp insertion on fpc4140, and a 536 bp deletion on fpc0011. (C) Six additional medium indels (one for each chromosome) in AF16 and HK104 genotypes. The indels are (Chr. 1-5 and X, from left to right): cb-m142 (250 bp), cb-m21 (200 bp), cb-m205 (300 bp), cb-m172 (250 bp), cb-m103 (200 bp), and cb-m204 (250 bp). (D) Validation of small indels in AF16, HK104 and F1 heterozygotes (AH). Chr 1: bhP19 (32 bp), Chr 2: bhP21 (20 bp), Chr 3: bhP12 (39 bp), Chr 4: bhP9 (16 bp), Chr 5: bhP48 (22 bp) and Chr X: bhP24 (16 bp).
The ssahaSNP program  was also able to detect insertion/deletion (indel) polymorphisms of 1-50 bp (Figure 1B). We used the parse_indel utility to extract 7,530 candidate indels (4,686 deletions and 2,844 insertions compared to AF16) for the HK104 strain. Most insertions and deletions detected by ssahaSNP were single base pair events; the largest was 49 bp. To identify larger insertion/deletion events, we developed a customized algorithm called BreakPointRead that detects indels based on BLAST alignments between read sequences and the reference genome. First, the algorithm identifies "breakpoint reads" with alignment gaps of 10 bp or larger compared to the reference sequence. Next, it analyzes the gap size and alignment orientations to infer the nature (insertion, deletion, inversion, etc.) and approximate size of sequence variation. When we applied BreakPointRead to the HK104 sequence traces, it identified 689 breakpoint reads suggestive of 635 underlying variants. We filtered the results to obtain insertion/deletion events between 50 and 2,000 bp. The resulting set contained 212 predicted indels (144 deletions and 68 insertions), the largest of which was a 1,707 bp deletion on chromosome IV (Figure 1C and additional file 3).
Validation of indels
We used the BreakPointRead algorithm to develop PCR fragment length polymorphism (PLP) assays for medium indels to facilitate high-throughput gene mapping. Of the 212 putative indels, we selected 40 for validation in AF16 and HK104 parental DNA (Figure 2B, C). Two other indels (bdP1 and bdP4) were also chosen for a similar analysis (see Methods and additional file 3). Ten of the 42 PLP assays failed PCR in one or both strains (Table 3). As mentioned earlier, some of these errors may be eliminated by redoing PCR using new primers. One assay (cb-m16) resulted in PCR products that were quite large (>1,500 bp) in both strains and therefore could not be accurately resolved on the gel. Of the remaining 31 assays that were successfully amplified, 29 confirmed the presence of polymorphism, 22 of which were similar in type and size as predicted (Table 3).
Table 3. List of medium indels tested by PCR.
We developed PLP assays for small (<50 bp) indels in AF16 and HK104 that were identified by ssahaSNP. To allow for gel resolution, we excluded indels smaller than 7 bp. This resulted in 436 assays that had 7-49 bp band size differences between AF16 and HK104 (see additional file 4). We tested 45 indels (between 4 and 9 for each chromosome) by PCR and found that except one (bhP44), for which HK104 amplification failed, all others could be successfully amplified (Table 4). A total of 32 indels showed bands of predicted sizes (Table 4, Figure 2D). Of the remaining 12, 1 showed no indel (i.e., identical PCR products in AF16 and HK104), 4 had multiple products (either due to PCR error, incorrect in silico predictions, or misassembly), and 7 showed PCR products that were inconsistent and unreliable (Table 4).
Format: XLS Size: 133KB Download file
This file can be viewed with: Microsoft Excel Viewer
Table 4. List of small indels tested by PCR.
Altogether we experimentally confirmed 75 polymorphisms (14 snip-SNPs, 29 medium indels, and 32 small indels) on all six chromosomes (Figure 3, Table 5). The utility of these 'working' markers in genetic studies is demonstrated by successful mapping of several mutations that cause visible phenotypes. In two cases these mapping experiments also helped improve the genome sequence assembly. Specifically, the bhP42 contig fpc4184 was placed near the center of chromosome 1 and bhP18 contig fpc4010 was reassigned to the left arm of chromosome 3 (see below).
Table 5. Summary of polymorphisms experimentally tested in this study.
Figure 3. A polymorphism-based genetic linkage map of C. briggsae. The map consists of 13 snip-SNPs (italics), 22 medium indels, and 32 small (bold) indels.
Linkage mapping studies using snip-SNPs and indels
The snip-SNPs and indels were used to map a set of 12 mutations with easily recognizable phenotypes (e.g., Uncoordinated or Unc and Dumpy or Dpy) that were previously isolated in different laboratories (Table 6 and Methods). Except lin(bh25) and unc(sy5415), all other loci were uniquely assigned to linkage groups by standard 2-factor mapping using known mutations that serve as reference (Table 6, also see ). The dpy(s1272), unc(s1270), and unc(sa997) are reference markers for LGIII, LGIV, and LGV, respectively. The remaining autosomal loci are linked to C. elegans orthologs Cbr-lin-11 (LGI) and Cbr-unc-4 (LGII). Not only did the polymorphism-based mapping agree with phenotypic marker-based mapping (see dpy(sy5001) and dpy(sy5148) in Figure 4A, B), it also helped to identify linkage groups of lin(bh25) (LGI, see Figure 4D and Table 6) and unc(sy5415) (LGV, see Table 6). In each of these cases a single cross with HK104 provided enough genomic DNA and usually one PCR per chromosome was sufficient to establish the linkage (using bulk-segregant approach, BSA). We also quantified DNA band intensities to determine linkages as
Table 6. List of mutations used in polymorphism mapping experiments.
Figure 4. Mutation mapping by polymorphisms. The indels were used to map dpy(sy5001) (A) and dpy(sy5148) (B, C) and snip-SNP (cb650) was used to map lin(bh25) (D). (A) Mapping of X-linked mutation dpy(sy5001) using six medium indels (one per chromosome). W, non-mutant (phenotypically wild type) pool; M, mutant pool. (B) dpy(sy5148) localization on chromosome 2 by small indels (Chr 1: bhP19, Chr 2: bhP21, Chr 3: bhP12, Chr 4: bhP11 and Chr X: bhP26). (C) ULVs for dpy(sy5148) show linkage to bhP21. X-axis shows chromosomes whereas Y-axis linkage values. The dotted line shows the baseline for unlinked chromosomes. (D) SacI digested PCR amplified genomic DNA of wild-type controls (A: AF16, H: HK104) and lin(bh25) mutant (M) and non-mutant (W) categories. There is a clear bias towards AF16 DNA (uncut) in the mutant pool compared to the non-mutant pool, demonstrating that lin(bh25) is linked to cb650 (chr. 1).
In addition to linking mutations to chromosomes we also investigated whether polymorphisms could be used in more precise mapping i.e., placing mutations in specific chromosomal regions (left, right arms, or middle). We reasoned that by narrowing down genetic intervals of mutations it should be possible to identify potential candidates, including C. elegans orthologs, thereby facilitating gene cloning by RNAi and transgene rescue approaches. To this end we used three medium indels on chromosome X to map dpy(sy5001). The ULVs for dpy(sy5001) suggest weak linkages to indels cb-m204 (left arm) and cb-m136 (right arm) and tight linkage to the middle indel cb-m197 (Figure 5A). Similar result was also obtained with the small indel bhP26 that is located close to cb-m197 and is strongly linked to dpy(sy5001) (data not shown).
Figure 5. Sub-chromosomal localization of mutations by medium and small indel-based mapping. The gel images show non-mutant (wild type, W) and mutant (M) pool of PCR amplified DNA from F2 worms. The histograms show ULVs for various indels. (A) dpy(sy5001) is most strongly linked to the medium indel cb-m197 (located roughly in the middle of chromosome X, ChrX-M) compared to flanking indels cb-m204 (left arm, ChrX-L) and cb-m136 (right arm, ChrX-R). (B) dpy(s1272), unc(sa972), and lin(bh20) are located on chromosome 3. While dpy(s1272) and unc(sa972) are strongly linked to bhP14 and bhP18 and appear to be on the left arm, lin(bh20) maps closer to bhP40 (center right region). (C) unc(sy5422) is tightly linked to indels bhP9 and bhP16 on the right arm of chromosome 4.
Besides dpy(sy5001), we also mapped 4 autosomal mutations to sub-chromosomal regions using small indels. The results showed that dpy(s1272) is most strongly linked to bhP18, unc(sa972) to bhP14 and bhP18, and lin(bh20) to bhP40 (all on chromosome 3) (Figure 5B). The unc(sy5422) appears to be located on the right arm of chromosome 4 (closer to the middle) since it shows strongest linkage to bhP9 (Figure 5C).
The C. briggsae sequence assembly 'cb3' had placed the bhP18 contig fpc4010 on the right arm of chromosome 3 (~36.5 mu). We found that this location was inconsistent with ULVs for dpy(s1272) and unc(sa972) since both mutations are also linked to the left arm indel bhP14 (Figure 5B). This suggested that there could be a possible error in the sequence assembly. In a separate study Zhao et al.  used the SNP-based oligonucleotide array to map dpy(s1272) and found that the region corresponding to bhP18 is misassembled. Therefore, we have assigned bhP18, as well as dpy(s1272) and unc(sa972), to the very end of the left arm of chromosome 3. Additional mapping using polymorphisms and visible markers will resolve the extent of misassembly.
In addition to the above BSA approach we also analyzed single F2 mutants segregating from a cross (termed single recombinant analysis or SRA) to determine recombination distances between polymorphisms and mutations. For this we used a multivulva mutation lin(sy5353) and an Unc mutation unc(sy5506). The lin(sy5353) mutation is located on chromosome 1 since it is strongly linked to three small indels bhP1, bhP7 and bhP42 (B.P. Gupta, unpublished results). From a cross between lin(sy5353) and HK104 we picked 23 F2 mutant animals (46 chromosomes) and analyzed their DNA for the presence of bhP1. A single recombinant chromosome was recovered suggesting a recombination frequency of ~2% between lin(sy5353) and bhP1. In the case of unc(sy5506) mutation, located on chromosome X, we analyzed 20 F2 mutant animals for the presence of indel bhP26. A total of 4 recombinant chromosomes were recovered (see additional file 5) suggesting that the two loci are 10% apart. These results demonstrate that SRA mapping protocol can be used in C. briggsae to localize mutations to sub-chromosomal regions and narrow down their genetic interval.
Additional file 5. Mapping of unc(sy5506) mutation by single recombinant analysis (SRA). Twenty single F2 mutant animals were individually examined by PCR for the presence of indel bhP26. Four of these (#7, #10, #16, and #17) were found to be recombinants, as judged by the presence of two bands on the agarose gel (corresponding to AF16 and HK104 DNA).
Format: PNG Size: 89KB Download file
Discussion and Conclusions
We took a bioinformatics approach to identify polymorphisms in the C. briggsae genome and experimentally validated a small set of these to facilitate mapping of mutations. Comparison of AF16 (reference strain) to four other isolates (HK104, HK105, VT847 and PB800) revealed that HK104 is most polymorphic since it has the highest density of SNPs among all the strains. Altogether we identified ~31,300 polymorphisms (23,800 SNPs and 7,500 medium and small indels) between AF16 and HK104 that promise to be a valuable resource for mutation mapping and genome evolution studies. Roughly 20% of the SNPs are predicted to alter restriction enzyme sites (snip-SNPs) that could be detected by PCR followed by restriction digestion and agarose gel electrophoresis.
A total of 107 polymorphisms (20 snip-SNPs, 42 medium indels, and 45 small indels) that were experimentally tested, 66 (12 snip-SNPs, 22 medium indels, and 32 small indels) showed DNA fragments identical (or close) to in silico predictions (Table 5). Another 9 cases (2 snip-SNPs and 7 medium indels) were significantly different but nonetheless showed the presence of underlying variants. In 15 cases no polymorphism could be detected. Thus, excluding PCR failures (total 17), the success rate of correctly predicted polymorphisms was 73% (69-86% range) (Table 5). This suggests that both ssahaSNP and BreakPointRead algorithms work equally efficiently regardless of the type of polymorphism in question. A similar study in C. elegans  showed that greater than 95% of the polymorphisms predicted by the Polybayes program  are true. It remains to be seen whether the lower success rate in C. briggsae is due to intrinsic differences between the programs alone or if the quality of sequence data and assembly are additional contributing factors.
We used snip-SNPs and indels to map 12 mutations with visible phenotypes, and found that polymorphism-based mapping agreed with phenotypic marker-based results. Furthermore, it helped map two mutations, lin(bh25) and unc(sy5415), for which no prior genetic linkage data was available. Five mutations were also localized to sub-chromosomal regions. Thus our mapping resource can be used to rapidly map new mutations in C. briggsae. It is also relatively easy to validate additional polymorphisms if one needs a greater resolution. It should be pointed out that Hillier et al.  have validated another set of 9 snip-SNPs by sequencing during the process of C. briggsae genome sequence assembly (see additional file 6). Given the high density of such markers (>2,000), it should be possible to map a mutation within a small genetic interval to facilitate molecular cloning (e.g., see ).
In addition to mapping mutations, SNPs and indels could also be used to improve the genetic linkage map of C. briggsae. The current C. briggsae sequence assembly, cb3, incorporates 90.2% (91.2 Mb) of the genome united into six chromosomes . The remaining 9.8% of sequences are tentatively associated with chromosomes. These unmapped regions could be integrated into chromosomes by polymorphism-based recombination mapping. We have successfully used this approach to place the contig fpc4184 in the vicinity of fpc3441 (chromosome 1) based on the recombination distance of 5% between bhP42 (fpc4184) and bhP1 (fpc3441) (Figure 3) (A. Seetharaman, P. Cumbo, B. Nagagireesh and B. P. Gupta, manuscript submitted). In the other case, we have reassigned the bhP18 contig fpc4010 to the left arm of chromosome 3 based on its tight linkage to dpy(s1272) and unc(sa972). Additional snip-SNPs and indels could further refine the locations of these contigs.
Strains and culture conditions
All strains were maintained at 22°C. The general methods of culturing nematodes are previously described . AF16 is the wild type genetic background for all strains. The four other wild-type C. briggsae isolates that were used for polymorphism discovery are HK104, HK105, VT847, and PB800. The HK104 strain was used in all mapping experiments. Various mutations used in this study are: csp-1(sa972), dpy(nm4), dpy(s1272), dpy(sy5001), dpy(sy5148), lev(sy5440), lin(bh20), lin(bh25), lin(sy5353), unc(sa997), unc(sy5415), unc(sy5422), and unc(sy5506).
In silico predictions of polymorphisms
SNP discovery was performed on 13,632 shotgun sequence traces from strains HK104, VT847, HK105, and PB800. The ssahaSNP program (version SSAHA2)  was used to call SNPs due to its robust and efficient performance; only polymorphisms with quality scores above the minimum threshold were accepted. We also tested Polyphred (v5.04)  and PolyBayes (v3.0)  programs but found that only ssahaSNP could efficiently handle the entire read set and reference genome sequences as input. For a reference sequence the cb25 genome sequence assembly, which is based on strain AF16 and organized into ultra (fingerprint) contigs, was obtained from Wormbase. Flanking sequences for predicted SNPs were repeat-masked to lower case using the RepeatMasker program (v3.1.5)  with a customized C. briggsae repeat library.
The HK104 SNPs were positioned on the cb25 sequence assembly during SNP discovery. To position them on the newer cb3 sequence assembly, which is by chromosome, we obtained the assembly AGP files from Wormbase. SNP positions were inferred based on the coordinates and orientation of their cb25 ultracontig. SNPs on cb3-unmapped ultracontigs were mapped by WU-BLAST v2.0 (Gish, W., personal communication) alignment of their flanking sequences. Some 699 SNPs could not be positioned on the cb3 assembly by either method.
Candidate AF16-HK104 indels were extracted from HK104 sequence traces using the parse_indel utility of ssahaSNP. In the HK104 set, the largest indel event identified by ssahaSNP was 49 bp. To identify larger insertion/deletion variants we implemented BreakPointRead, a custom algorithm that detects structural variations (insertions, deletions, inversions, and copy number variants) spanned by individual sequence traces. Traces were aligned to the reference genome (cb25) using WU-BLAST v2.0 (Gish, W., personal communication), and screened for alignments with "gaps" of > = 10 bp. The alignment patterns of such "breakpoint reads" were used to infer the type and size of polymorphism. Predicted insertions and deletions were set aside for assay development.
The bdP polymorphisms described in this study (snip-SNP bdP3 and medium indels bdP1 and bdP4) were identified in the laboratory of SEB. The snip-SNP bdP3 was earlier used in a study involving ray pattern variation in C. briggsae .
Development of RFLP and PCR Length Polymorphism (PLP) assays
SNPs were screened for substitutions that altered the recognition sequence of restriction enzymes using the Bio::Restriction::Analysis library of BioPerl . The analysis was limited to 30 restriction enzymes from REBASE  known to be reliable and inexpensive. PCR assays were designed (amplicon sizes of 500 to 1000 bp, primer Tm's of 54-56 °C) using a local installation of the primer3 program . In silico fragment analysis of the PCR products was performed to predict band sizes for AF16 and HK104; assays with more than 4 bands in either strain were removed.
In the case of small indels (7-49 bp), primers were selected to generate AF16 amplicon sizes within the range of 200 and 400 bp. For medium indels (50-2,000 bp), primers flanking each indel and specifying an AF16 amplicon size of 300-800 bp were selected.
In all experiments the genomic DNA from F2 worms (derived from a cross between AF16 and HK104 animals) was used as a PCR template. In some control experiments genomic DNA from F1 heterozygous animals was also used. PCR results that gave rise to unexpected or no products were repeated at least twice. In some cases we also tested different annealing temperatures. Those that consistently failed were termed as "PCR failure".
For RFLP and medium indel assays, 10 ng dry PCR primers (IDT, Coralville, IA) were resuspended into 40 μM in a 96-well format. Our PCR mixture consisted of 300 pg genomic DNA, 0.1 μM PCR of each up and down primer, 0.02 U Platinum Taq (Invitrogen), 83.3 μM per base dNTPs, 2.92 μM MgCl2, 10× Buffer (16.7 mM Tris-HCl pH 8.4, 41.67 mM KCl - Invitrogen PCR Kit), and 4.2% DMSO. Amplification was carried out in a Perkin-Elmer 384 PCR plate containing 12 μl of 1× PCR mixture. The 384-well plate was sealed with Microseal A (MJ Research) before carrying out the PCR. After initial heat denaturation step (95°C - 2 min.) we used a fixed 35 cycles PCR (94°C - 10 sec., 58°C - 20 sec., 68°C - 30 sec.) and a final extension (68°C for 10 min.). PCR products were analyzed on 10% polyacrylamide gels that consisted of 33.3% of 29:1 acrylamide:bis (Biorad), 10% of 10×TBE, 15% of glycerol, 40.2% H20, 3.4 × 10-2% TEMED (Int'l Biotech), and 1.42% of 10%APS. Gels were stained with SYBR green (Invitrogen) and inspected over a UV light box at 254 nm.
To detect small indels (<50 bp) we used a standard 35 cycles PCR (94°C - 10 sec., 48°C - 30 sec., 72°C - 60 sec.) for all amplifications. For BSA-based mapping the genomic DNA was prepared from 25 adults in a 10 μl lysis buffer (consisting of 50 mM KCl, 10 mM Tris pH 8.2, 2.5 mM MgCl, 0.45% Tween 20, 0.45% NP40, 0.001% Gelatin, and 30 μg Proteinase K). The mixture was frozen at -80°C (at least 30 min) and then placed in a thermal cycler for 1 hr incubation at 60°C followed by 15 min heat inactivation at 95°C. The resulting genomic DNA was diluted to 25 μl using sterile distilled water and stored at -20°C. The PCR mixture (25 μl) contained 1 μl genomic DNA, 1 μl each of up and down primers, 1 μl dNTPs, 2.5 μl NEB ThermoPol 10× PCR buffer, 0.2 μl NEB Taq Enzyme, and 18.3 μl sterile distilled PCR grade water. For SRA-based mapping, single worms were placed in 5 μl lysis buffer and processed as above. The lysed samples were used as DNA templates in PCR experiments. The amplified products were first analyzed on 1% agarose gel (Invitrogen UltraPure, Catalog #15510-027). Successful amplifications were subsequently examined on a 4% high-resolution agarose gel (Invitrogen UltaPure Agarose-1000, Catalog #10975-035) to determine the presence of indels.
We picked 12 strains for linkage mapping studies (Table 6). The mutations were obtained from EMS (ethyl methane sulfonate) mutagenesis screens in an AF16 genetic background in various laboratories. The strains were outcrossed several times (3 or more). For mapping, mutant hermaphrodites were crossed with HK104 males and the genomic DNA from 20 F2 animals (wild type and mutant separately) was prepared as described in the previous section. The linkage was determined by PCR using protocols established for control experiments.
Linkage and ULV analysis
To determine the linkage of a mutation to a chromosome, we initially relied on the visual inspection of DNA band intensities on Ethidium bromide-stained agarose gels. Subsequently, in indel-based mapping experiments, we calculated linkages as unitless values (ULVs) for an unbiased analysis. The mean intensities of DNA bands were measured by NIH ImageJ software (version 1.41o; ) using Measure tool under Analyze menu. For each mutation a ratio of band intensities in the "mutant" lane was calculated by dividing the mean intensities of the AF16 bands by the mean intensities of the HK104 bands. This ratio was termed as the ULV. As expected, ULVs were one for unlinked mutations and higher for linked mutations.
Genetic positions of polymorphisms
The genetic positions of snip-SNPs and indels in this study correspond to nearest SNPs that were experimentally validated (D.C.K. and R.D.M., unpublished). These 'verified' SNPs (400 in total) were genotyped in RILs derived from two independent crosses between AF16, HK104 and VT847 (AF16 × HK104 and AF16 × VT847). The details are available on the C. briggsae SNP Research Facility website .
Polymorphisms and sequence data availability
All identified snip-SNPs, indels and PCR primers, confirmed or otherwise, are accessible via the C. briggsae resource website (see "Polymorphism" link) . The website also contains sequence reads of all four species (HK104, VT847, HK105, and PB800). The sequence directories are organized into Polyphred structure and contain additional files (such as read quality). Additional information on polymorphism discovery using sequence data can be obtained from the Washington University C. briggsae SNP Research Facility website . The confirmed polymorphisms have also been submitted to Wormbase .
List of abbreviations
BSA: Bulk segregant analysis; Indel: Insertion-deletion; PLP: PCR fragment length polymorphism; RIL: Recombinant inbred line; RFLP: Restriction fragment length polymorphism; SNP: Single nucleotide polymorphism; SRA: Single recombinant analysis; ULV: Unitless linkage value.
The C. briggsae Advisory Group (BPG, ESH, HMC, RDM, and SEB) planned the project and collected the genetic and financial resources to perform the work. DCK and RDM carried out bioinformatics searches. DCK designed snip-SNP and PLP assays. JES, BT and KH performed PCR experiments to validate polymorphisms and map mutations. BT and BPG calculated ULVs for mutations linked to medium and small indels. bdP polymorphisms were identified and validated in the laboratory of SEB. BPG wrote the paper with input from ESH, HMC, SEB, DCK, and JES. All authors read and approved the final manuscript.
The authors dedicate this work to the memory of our dear colleague and friend, Dr. Raymond D. Miller. We thank Prateek Goyal and Nagagireesh Bojanala for assistance in mapping some of the mutants. We also thank Paul Sternberg (California Institute of Technology), David Baillie (Simon Fraser), and James Thomas (University of Washington Seattle) for sharing strains and unpublished data on mutant isolation and genetic mapping. Many of the oligos were obtained from the MOBIX lab at McMaster University. This work was supported by a grant from National Institutes of Health (R24 GM075101) to RDM and BPG.
Hillier LW, Miller RD, Baird SE, Chinwalla A, Fulton LA, Koboldt DC, Waterston RH: Comparison of C. elegans and C. briggsae Genome Sequences Reveals Extensive Conservation of Chromosome Organization and Synteny.
Methods Mol Biol 2000, 132:365-386. PubMed Abstract
C. briggsae SNP Research Facility [http://snp.wustl.edu/snp-research/c-briggsae/genetic-map-of-c-briggsae.html] webcite