Skip to main content

An accurate method for quantifying and analyzing copy number variation in porcine KIT by an oligonucleotide ligation assay



Aside from single nucleotide polymorphisms, copy number variations (CNVs) are the most important factors in susceptibility to genetic disorders because they affect expression levels of genes. In previous studies, pyrosequencing, mini-sequencing, real-time PCR, invader assays and other techniques have been used to detect CNVs. However, the higher the copy number in a genome, the more difficult it is to resolve the copies, so a more accurate method for measuring CNVs and assigning genotype is needed.


PCR followed by a quantitative oligonucleotide ligation assay (qOLA) was developed for quantifying CNVs. The accuracy and precision of the assay were evaluated for porcine KIT, which was selected as a model locus. Overall, the root mean squares of bias and standard deviation of qOLA were 2.09 and 0.45, respectively. These values are less than half of those in the published pyrosequencing assay for analyzing CNV in porcine KIT. Using a combined method of qOLA and another pyrosequencing for quantitative analysis of KIT copies with spliced forms, we confirmed the segregation of KIT alleles in 145 F1 animals with pedigree information and verified the correct assignment of genotypes. In a diagnostic test on 100 randomly sampled commercial pigs, there was perfect agreement between the genotypes obtained by grouping observations on a scatter plot and by clustering using the nearest centroid sorting method implemented in PROC FASTCLUS of the SAS package. In a test on 159 Large White pigs, there were only two discrepancies between genotypes assigned by the two clustering methods (98.7% agreement), confirming that the quantitative ligation assay established here makes genotyping possible through the accurate measurement of high KIT copy numbers (>4 per diploid genome). Moreover, the assay is sensitive enough for use on DNA from hair follicles, indicating that DNA from various sources could be used.


We have established a high resolution quantification method using an oligonucleotide ligation assay to measure CNVs, and verified the reliability of genotype assignment for random animal samples using the nearest centroid sorting method. This new method will make it more practical to determine KIT CNV and to genotype the complicated Dominant White/KIT locus in pigs. This procedure could have wide applications for studying gene or segment CNVs in other species.


Susceptibility to genetic disorders is known to be associated not only with single nucleotide polymorphisms (SNP), but also with structural and other genetic variations, including copy number variations (CNVs) [13]. Therefore, once identified, a CNV needs to be analyzed at the locus level, and ultimately, the genotype and haplotype must be determined to elucidate its relationship with a particular genetic alteration. Pyrosequencing, mini-sequencing, real-time PCR and invader assays are among the techniques that have been used to detect CNVs [46].

The porcine KIT was selected for this study because it is a well characterized and functionally important CNV. The Dominant White/KIT locus that determines white coat color is located in Sus scrofa chromosome 8 (SSC8) [7, 8]. Two KIT mutations cause the Dominant White phenotype in pigs: a gene duplication associated with a partially dominant phenotype, which is depicted as normal and duplicated in Figures 1a and 1b, and a splice mutation leading to the fully dominant allele [7, 9], which is marked in Figure 1a as an SNP(G/A) at the first nucleotide of intron 17. As shown in Figure 1c, there are four known major alleles at the KIT locus: the recessive i allele for the Color phenotype, the IPallele for the Patch phenotype, the dominant I allele for the White phenotype and IBefor the Belt phenotype [10, 11]. I allele diversity has been reported and classified in detail as I1, I2, I3 and IL[4]. All possible genotypes, which are derived from the alleles shown in Figure 1c, and theoretical ratios of spliced and duplicated copies corresponding to each genotype are presented in Table 1. The two ratios of each polymorphism were used as reference values when the genotypes of experimental pig samples were assigned in this study.

Figure 1
figure 1

(a) A schematic description of tandem duplication at the porcine KIT locus. A duplication unit is about 450 kb. A breakpoint at the junction between the two KIT copies is designated as duplicated, and another breakpoint at the distal end of the 2nd KIT copy is designated as normal. A splice donor mutation in intron 17, by which KIT becomes the fully dominant allele, is indicated by an arrow. The repeat elements around the breakpoint, L1MC1 and L1ME1, are described. The target points of pyrosequencing (Pyro_splice) for quantifying copies with the splice mutation and the quantitative oligonucleotide ligation assay (qOLA_CNV) for detecting total copy numbers are marked. (b) Nucleotide sequence around the breakpoint. While the breakpoint of the normal copy is on L1MC1, that of the duplicated copy is at the junction between L1MC1 and L1ME1. The common primer (BPT_Com) for the quantitative oligonucleotide ligation assay (qOLA_CNV) is marked in underlined italic letters and the two specific primers (BTP1 and BPT2) are indicated by underlined plain letters. The two nucleotides in a box, C and G, indicate the breakpoint and comparison point for qOLA_CNV in this study. Small letters g and c indicate the comparison point for measuring KIT CNV by pyrosequencing (Pyro_CNV) [4]. (c) Schematic descriptions of KIT alleles. The two target points of pyrosequencing (Pyro_splice) for quantifying copies with the splice mutation, and qOLA_CNV for detecting total copy numbers are marked with arrows. A question mark in the IBE allele means an unidentified polymorphism causing the Belt phenotype. Discrimination between i and IBE is not possible at present.

Table 1 Theoretical genotype description of the KIT locus by the splice mutation and copy number variation.

To analyze the KIT locus, RFLP [9], minisequencing, real-time PCR [12], invader and pyrosequencing assays [4] have been used. Pyrosequencing has provided the best resolution for quantifying KIT CNV giving more accurate results than real-time PCR amplification and invader technologies. However, as the copy number increases, it gradually becomes more difficult to use the pyrosequencing method to accurately distinguish among genotype classes that differ by only one copy. This is because the relative increase in the signal from the duplicate breakpoint becomes smaller [4]. An underestimated CNV ratio may result in an ambiguous genotype assignment in samples for which family information, including parental genotypes, is not available.

We have therefore developed PCR followed by a quantitative oligonucleotide ligation assay (qOLA), which gives high resolution data for determining KIT CNV, especially if the copy number is high (>4). The development of qOLA is based on the strategy previously described in [4], but it improves on the pyrosequencing method [4] for analyzing CNV of the locus. We have also established a nearest centroid sorting procedure to verify the reliability of the genotype assignment for random animal samples. The qOLA used on a platform with an ABI sequencer is sensitive enough to analyze DNA from a few hair follicles, so DNA from various sources could be used for qOLA.


Verifying the specificity of the PCR primers used for analyzing KIT CNV

The PCR primers designed for the published pyrosequencing method [4] were used in this study. The primer sequences selected from the KIT duplication breakpoint are located on repetitive elements, L1MC1 and L1ME1 (Fig. 1a and 1b). The forward primer (KITBPF) shows 80% sequence identity with the L1MC1 consensus sequence and the two reverse primers, KIT1BPR for the normal copy and KIT2BPR for the duplicated copy, show 63.2% and 94.7% sequence identity with L1MC1 and L1ME1, respectively. This finding raised the question of whether the PCR products may contain nonspecific amplification products from other genomic regions. To evaluate the specificity of the PCR primers, somatic cell hybrid panel mapping was performed prior to the quantification assay. The two amplicons were located in SSC8p11, where the KIT locus exists (assignment probability/correlation: 0.8789/0.9250 for normal and 0.8791/0.9250 for duplicated), indicating that the amplifications of the primer sets were specific. As shown in Additional file 1, the primer sets were clearly amplified.

Evaluation of the established qOLA to measure the CNV of KIT (qOLA_CNV)

The amplicons of the duplicated and normal copies were cloned into the pCR®2.1-TOPO vector (Invitrogen, USA). The cloned amplicons were re-amplified using the M13 forward and reverse primers, and were then purified and serially diluted from 0% to 100% duplicated copy vs. normal copy. PCR followed by qOLA_CNV was performed on four replicates, and two standard curves were obtained for peak height and peak area (Additional file 2a and 2b). Correlation coefficients were 0.999 for both standard curves, indicating very good linearity. However, the correlation coefficient is only an index of the linearity of the standard curve. Bias, which indicates the amount of systematic error from the reference ratio, and the standard deviation (SD), which indicates the variation between replicates for a reference point, are more appropriate indices for the fit of the observed data to the expected results. As shown in Table 2, peak height values in qOLA_CNV fit the reference values better and show the least variation. In particular, for accurate genotyping of individuals with a total of more than 4 KIT copies, the assay needs better resolution in the zone between 60% and 90% in the standard curve. In this zone, the peak area values in qOLA_CNV showed root mean squares (RMS) of the bias and SD of 2.21 and 1.17, respectively. In contrast, the RMS of the bias and SD of the peak height values in qOLA_CNV were 0.86 and 0.51, respectively, in the same zone. qOLA_CNV was compared with the published pyrosequencing assay [4] for KIT CNV (Pyro_CNV). The same serial dilutions used for qOLA_CNV were used to obtain the standard curve for Pyro_CNV. As shown in Additional file 2c, the standard curve of Pyro_CNV showed good linearity (correlation coefficient 0.995). However, as shown in Table 2, the overall RMS of the bias (5.05) and SD (1.04) in Pyro_CNV were more than twice those for the peak height in qOLA_CNV (2.09 and 0.45). In conclusion, CNV estimation for porcine KIT using the peak height values in qOLA_CNV showed the lowest systematic errors and variations (Additional file 2d) of the studied methods, and therefore was used in further experiments to analyze KIT CNV and assign genotypes.

Table 2 Bias and standard deviation (SD) of each method for each duplicated copy ratio. Root mean square (RMS) of the bias and SD is calculated to compare accuracy and precision among the three standard curves.

Another pyrosequencing assay [12] (Pyro_Splice) for quantifying KIT copies with a splice donor mutation in intron 17 was combined with qOLA_CNV. qOLA_CNV gives information about the total KIT copy numbers in a sample, and Pyro_Splice gives additional information about the ratio of spliced KIT copies to the total copies estimated by qOLA_CNV. As shown in Table 1, several different genotypes have identical ratios in each polymorphism. Therefore, combined Information from the two polymorphisms should yield better discriminating power in assigning genotype.

CNV and genotyping tests using a combination of qOLA_CNV and Pyro_Splice (Fig. 2e) were performed to verify KIT allele segregation. One hundred and forty-five F1 animals produced by a cross between Korean native and Landrace pigs were used. By combining the qOLA_CNV and Pyro_Splice assays, we were able to resolve the genotypes in the founder and F1 populations. The 19 Korean native pigs, 8 boars and 11 sows were all homozygous i/i (Fig. 2a) for the Black phenotype (Additional file 3a), whereas the 17 Landrace pigs, 8 boars and 9 sows had three different genotypes (Fig. 2a) for the White phenotype (Additional file 3b), which consisted of combinations of the I1, I2 and i alleles. For genotype assignment we used two classification methods based on the clusters of measurements on a scatter plot (Fig. 2a and 2b) and the clusters of observations at 12 seed points using nearest centroid sorting implemented in PROC FASTCLUS of the SAS package [13, 14]. The genotyping results showed 100% agreement between the two methods and were in good agreement with both the theoretical genotype ratios (Table 1 and Fig. 2c) and phenotypes (Additional file 3c). The alleles of the F1 animals were consistent with the genotype of each Landrace founder, and were clustered at distinguishable ratios. This result indicates that genotype assignment based on qOLA_CNV in combination with Pyro_Splice is highly reliable and could be applied to random populations without pedigree information.

Figure 2
figure 2

Genotype assignment using a combination of qOLA_CNV to analyze KIT copy numbers and Pyro_Splice to quantify KIT copies with spliced forms. (a) Genotype assignment of parental animals. As the Patch and Belt phenotypes were not presented in the F1 population, all Korean native pigs were assigned to i/i, and Landrace pigs were assigned to I1/I1, I1/I2 and I1/i by the clustering measurements. (b) Genotype assignment of F1 animals. (c) Electrophoregram in qOLA_CNV and pyrogram in Pyro_Splice as representative examples of the three genotypes in the F1 population. In qOLA_CNV, the ratio was (duplicated copy/total copy) = [G/(G + C)]. In the Pyro_Splice assay, the ratio was (spliced copy/total copy) = [A/( A + G )].

Diagnostic tests for randomly collected samples

Twelve standard coordinates and centroids representing 15 genotypes are presented in Table 1, and are marked on the plot in Figure 3a. There are three pairs of genotypes that have the same ratios of spliced and duplicated copies.

Figure 3
figure 3

Diagnostic tests on random samples using the combined method of qOLA_CNV and Pyro_Splice. The x-axis is the ratio of duplicated copy to total copy measured by qOLA_CNV. The y-axis is the ratio of spliced copy to total copy measured by Pyro_Splice. Different genotypes are indicated by different symbols. Genotyping was performed in two ways: genotype assignment on the basis of clusters on the scatter plot, and statistical classification by clustering at 12 seed points using the FASTCLUS procedure. (a) Twelve standard coordinates and seed numbers corresponding to 15 genotypes derived from theoretical ratios in Table 1. Seed numbers used in the statistical analysis are indicated on the symbol. (b) Genotyping of 100 randomly sampled commercial pigs produced by a three-way cross using the Duroc × F1 sows (Landrace × Large White or vice versa). There was perfect agreement between the two genotyping procedures. (c) Genotyping of 159 unknown Large White pigs. The genotyping results by clustering on the scatter plot were the same with those by the statistical method except for two individuals (I3/I3 vs. I1/I3 and I1/I2 vs. I2/I2), which are indicated by arrows.

One hundred commercial pigs produced by the cross between Duroc boars × F1 sows (Landrace sire × Large White dam or vice versa) were tested for CNV and then genotyped. The terminal sire, Duroc, is known to be a recessive homozygote, i/i [9]. F1 sows are phenotypically white, and are supposed to have the dominant I allele, either homozygous or heterozygous. Commercial pigs must have at least one i allele inherited from Duroc, and this feature allows us to predict the distribution of genotypes in the commercial pig population. Diagnostic testing of the population also facilitated evaluation of the established qOLA_CNV for low KIT copy numbers (<4). Estimated genotypes on the scatter plot were clustered into the four classes I1/i, I2/i, IP/i and i(IBe)/i (Fig. 3b), as expected, and this clustering were perfectly agreed with the genotypes obtained by statistical clustering.

One hundred and fifty-nine Large White pigs were tested. A wider range of CNV (3–6 copies) was expected in the Large White population than in the commercial pig population, so the resolution of the qOLA_CNV method could be evaluated. Furthermore, the robustness of qOLA_CNV could be addressed, because the DNA from Large White pigs was crudely prepared from hair follicles. First, we performed genotype classification on the basis of grouping observations on the scatter plot; the Large White pigs tested were clustered into 8 classes representing 11 genotypes (Fig. 3c). Secondly, they were genotyped by the statistical method. As shown in Additional file 4, both genotyping methods showed the same results, except for two individuals (IDs: T7183 and Y4410Q). The animal T7183 was scored 63.4% and 62.9% by the qOLA_CNV and Pyro_Splice methods, respectively, and the actual coordinate was clearly different from the I1/I3 standard coordinate (3:5 of duplicated copy to total copy, and 3:5 of splice copy to total copy). However, the Euclidean distance to the centroid No.12 (4:6 and 4:6) is 4.78, which is slightly further than 4.45 to the centroid No.11 (3:5 and 3:5). There was also a discrepancy in the assignment of animal Y4410Q between the two method: I2/I2(4:6 and 2:6) by clustering on the plot and the centroid No.8 (3:5 and 2:5) by the statistical method.

The two assays, qOLA_CNV and Pyro_CNV, are compared again in Table 3. The coefficient of variation (CV) of seventy I1/i commercial pigs was 2.9% by qOLA_CNV and 6.9% by Pyro_CNV. The CV of Pyro_CNV measurements was about twice that in qOLA_CNV, consistent with the comparison between the two standard curves. The Pyro_CNV data using hair DNA of the Large White pigs varied widely (CV range: 7.3–15.3%).

Table 3 Comparison of the accuracy and precision among qOLA_CNV, Pyro_CNV and Pyro_Splice in the genotyping of 100 randomly sampled commercial and 159 Large White pigs.


Ligation detection methods, such as the oligonucleotide ligation assay and ligase chain reaction, are widely used for detecting viral and microbial infections in clinical examination [1517] due to their high sensitivity. They are also widely used for typing the SNPs involved in genetic disorders [18] and in developing cancer cells [19]. Because neither approach is a quantitative assay, but is instead used for negative/positive screening, this study suggests another application for the ligation detection reaction.

In a previous study [4], Pyro_CNV was shown to be particularly useful for distinguishing individuals with two copies from those with three or more copies of the duplicated KIT region. However, Pyro_CNV has two drawbacks. First, it is over-influenced by the quality of DNA. The results within each class of copy number tended to form a line rather than a cluster in this study, which suggests that variations in DNA quality tend to influence the outcome of the test. This feature was most evident when crude Large White DNA samples prepared from hair follicles were used. The other drawback is that the two nucleotide (nt) positions used to estimate the duplication ratio are not the same distance from the pyrosequencing initiation point. As shown in Figure 1b, the 2nd nt, G, for the normal copy and 5th nt, C, for the duplicated copy are used as these positions. From the 1st to 3rd nt positions of the duplicated copy, there is a triple G. The pyrosequencing reaction for the triple G occurs simultaneously, which makes accurate comparison at the breakpoint difficult. There would be more variability in the efficiency of incorporation of dNTPs if the pyrosequencing reaction point were further from the initiation point. When Pyro_Splice was used for quantitative analysis of KIT copies with spliced forms, the comparison point SNP (G/A) was at the pyrosequencing initiation site, and the results showed much less deviation from the theoretical genotype description than Pyro_CNV (Table 2). This study verified that qOLA_CNV is suitable even for crudely prepared DNA samples. Moreover, the ligation point is the duplication breakpoint. Consequently, the assay established here overcomes the two difficulties in the Pyro_CNV assay. An additional advantage is that lower amounts of PCR amplicons are required for the assay than for pyrosequencing. We applied 25 cycles for blood DNA and 27 cycles for hair DNA, and used 1 μL of the PCR product for qOLA_CNV. Knowledge of PCR kinetics suggests that fewer PCR cycles should be better for post-PCR and end point analysis because the amplification of the two target amplicons is less biased.

Two methods were used to assign KIT genotype in this study. Genotyping through clustering observations on a scatter plot is simple and easy, but the subjective view of a researcher may enter into the genotyping of observations that are located on the boundary between neighboring clusters. Another additional assignment by statistical analysis using nearest centroid sorting is an objective method to verify clusters on the scatter plot. There were two animals showing different genotypes by the two assigning methods. In particular, the genotype disagreement of T7183 pig and given the CV figures of qOLA_CNV in Table 3 suggest that the distribution of I1/I3 and I3/I3would markedly overlap with means about 6% apart and SD about 2%, if a large number of animals with these genotypes were tested using the developed assay in this study. The samples that are not identically genotyped by the two assigning methods need to be reanalyzed. Incorrect ratios of qOLA_CNV and Pyro_Splice can be generated by mixed hair or blood samples, cross-contamination of purified DNAs, PCR and OLA errors, etc. As the possibility of having additional SNPs on the binding sites of the PCR and OLA primers is not excluded, sequence analysis of the sites would be required.

Genotypes assigned by quantifying porcine KIT CNV can be directly applied in the pig industry. White commercial pigs are commonly produced in some European and Asian countries, including Korea, where consumers prefer pork from white pigs. Most Korean pig producers perform a three-way cross using Duroc boars × F1 sows (Landrace sire × Large White dam or vice versa). In the cross, the state of fixation of the Dominant White/KIT allele in Large White and Landrace pigs is crucial for maintaining the White phenotype in the commercial pig population. Therefore, it was necessary to develop a high resolution assay to analyze a range of KIT CNV. The assay must be applicable for unknown samples that lack pedigree information and parental genotypes. The established qOLA_CNV successfully resolved KIT copy numbers from two to six and correctly genotyped unknown Large White samples.


We have established a reliable assay for measuring tandem CNV that could be applied for a variety of samples, such as those in a known pedigree, those with predictable segregation, those without pedigree information, and genomic DNA of poor quality. Combining this method with a verification procedure using statistical clustering, genotypes can be successfully assigned with high confidence. This development could be widely applicable to studies of the function and mechanism of CNV in other species, and may be particularly useful for tandemly repeated CNV.


Animals and DNA extraction

One hundred and forty-five F1 animals from a reciprocal cross between Korean native and Landrace pigs were used to verify the reliability of the qOLA_CNV through the analysis of KIT allele segregation. One hundred randomly selected commercial pigs produced by a three-way cross using the Duroc, Landrace and Large White breeds and 159 Large White pigs were tested for KIT CNV and genotyped for the KIT locus. Genomic DNA of the F1 and commercial pigs was extracted from whole blood by a red cell lysis-proteinase K method [20]. The DNA from the Large White pigs was prepared from hair follicles using 5% Chelex [21]. DNA recovered from blood was dissolved in TE buffer (pH 8.0) before use and the hair DNA was used directly as a PCR template.

Somatic cell hybrid panel mapping

In order to verify the specificity of the PCR primers, KITBPF, KIT1BPR and KIT2BPR [4], PCR analysis was performed using a porcine × rodent somatic cell hybrid panel [22]. Each PCR reaction across the clones in the panel was carried out in a total volume of 25 μL containing 25 ng of template DNA, 10 mM Tris-HCl (pH 9.0), 40 mM KCl, 2 mM MgCl2, 20 pmol of each forward and reverse primers (KITBPF and KIT1BPR; KITBPF and KIT2BPR), 100 μ M of each dNTP and 1 unit of Taq DNA polymerase (GenetBio, Korea). The PCR results were analyzed using the interpreting web pages at INRA [23].

qOLA for KIT duplication analysis (qOLA_CNV)

PCR was performed in a total volume of 25 uL with 20 ng of genomic DNA, 10 pmol of KITBPF [4] and tail primer [4], 0.1 pmol of each reverse primer (KIT1BPR and KIT2BPR) [4] containing a tail primer binding site, 200 μ M of each dNTP, 10 mM Tris-HCl (pH 9.0), 40 mM KCl, 2 mM MgCl2 and 1 unit of Taq DNA polymerase (GenetBio, Korea). The conditions were 27 cycles (for hair DNA) or 25 cycles (for blood DNA) of 20 s at 94°C, 20 s at 58°C and 30 s at 72°C. To remove the used Taq DNA polymerase activity, proteinase K was added to the PCR products (final conc. 10 μg/μL) and incubated at 55°C for 40 min. Then, the mixture was further incubated at 98°C for 40 min to get rid of proteinase K activity.

Oligonucleotides for qOLA_CNV were designed as shown in Figure 1b. A common oligonucleotide (BPT_Com: 5'-Fluorescein-GGC TAC ATA CTG TAT GAT TCC AA-3') spans from the -24 to the -1 nt position from the duplication breakpoint. One of each specific oligonucleotide for the duplicated or normal copy (BPT1: 5'-Phosohate-GGG TCA TGG CTT GAA AAA GAA AAA AAA AAA-3' and BPT2: 5'-Phosphate-CGA TAT GAC ATT CTG GAA ATA AAA AAA AAA AAA AA-3') was designed from the duplication breakpoint. The 5' end of BPT_Com was labeled with fluorescein, and a phosphate group was added at the 5' ends of BPT1 and BPT2. (dA)10 and (dA)15 were added at the 3' ends of BPT1 and BPT2, respectively, to distinguish the two OLA products and separate them further from the unused oligonucleotide peak at the size fractionation step.

qOLA was performed in 10 μL containing 1 μL of PCR product, 20 mM Tris-HCl (pH8.3), 25 mM KCl, 10 mM MgCl2, 0.5 mM NAD, 0.01% Triton X-100, 1.5 units of Ampligase (EPICENTRE Biotechnologies, USA), 1 pmol BPT_Com, and 0.5 pmol of each BPT1 and BPT2 primers. A cycled ligation reaction for qOLA can be applied because the Ampligase enzyme is a thermostable DNA ligase. The conditions for qOLA were 10 cycles at 94°C for 30 s and 50°C for 1.5 min. Ten μL of deionized formamide was added to the OLA and the products were resolved on an ABI Prism 3100 Genetic Analyzer (Applied Biosystems, USA). GeneScan software version 3.7 (Applied Biosystems, USA) was used to measure the height and area of each peak.

Pyrosequencing for analyzing KIT duplication (Pyro_CNV)

PCR and pyrosequencing to quantify KIT CNV was performed according to the published protocol [4] using PyroMark MD (Biotage, Sweden).

Quantitative analysis for KIT copies with spliced forms (Pyro_Splice)

The PCR primer set, KIT21 and KIT35, and the thermal conditions described in [9] were used. Pyrosequencing was performed as previously described [12].

Statistical analysis

Data were collected and basic statistics were calculated using the Microsoft Excel program (Microsoft, USA). Accuracy and precision were determined and compared among standard curves of qOLA_CNV, Pyro_CNV and Pyro_splice by the RMS of the bias and standard deviation, which were calculated using Minitab software (Minitab Inc., USA). The observations from the genotyping of 100 commercial and 159 Large White pigs were classified using nearest centroid sorting [13] implemented in PROC FASTCLUS of the SAS package [14]. The ratio of the spliced forms and the ratio of duplication in theoretical genotypes were selected as the 12 cluster seeds (the class centroids). Observations were assigned to the nearest class centroids on the basis of Euclidean distance. The clusters were then labeled as the class labels of the centroids.


  1. Redon R, Ishikawa S, Fitch KR, Feuk L, Perry GH, Andrews TD, Fiegler H, Shapero MH, Carson AR, Chen W, Cho EK, Dallaire S, Freeman JL, Gonzalez JR, Gratacos M, Huang J, Kalaitzopoulos D, Komura D, MacDonald JR, Marshall CR, Mei R, Montgomery L, Nishimura K, Okamura K, Shen F, Somerville MJ, Tchinda J, Valsesia A, Woodwark C, Yang F, Zhang J, Zerjal T, Zhang J, Armengol L, Conrad DF, Estivill X, Tyler-Smith C, Carter NP, Aburatani H, Lee C, Jones KW, Scherer SW, Hurles ME: Global variation in copy number in the human genome. Nature. 2006, 444: 444-454. 10.1038/nature05329.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  2. Stranger BE, Forrest MS, Dunning M, Ingle CE, Beazley C, Thorne N, Redon R, Bird commercial pigs de Grassi A, Lee C, Tyler-Smith C, Carter N, Scherer SW, Tavare S, Deloukas P, Hurles ME, Dermitzakis ET: Relative impact of nucleotide and copy number variation on gene expression phenotypes. Science. 2007, 315: 848-853. 10.1126/science.1136678.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  3. Kehrer-Sawatzki H: What a difference copy number variation makes. Bioessays. 2007, 29: 311-313. 10.1002/bies.20554.

    Article  CAS  PubMed  Google Scholar 

  4. Pielberg G, Day AE, Plastow GS, Andersson L: A sensitive method for detecting variation in copy numbers of duplicated genes. Genome Res. 2003, 13: 2171-2177. 10.1101/gr.1188203.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  5. Nevilie M, Selzer R, Aizenstein B, Maguire M, Hogan K, Walton R, Welsh K, Neri B, de Arruda M: Characterization of cytochrome P450 2D6 alleles using the Invader system. Biotechniques. 2002, : Suppl 34-Suppl 43.

    Google Scholar 

  6. Aldred PM, Hollox EJ, Armour JA: Copy number polymorphism and expression level variation of the human alpha-defensin genes DEFA1 and DEFA3. Hum Mol Genet. 2005, 14: 2045-2052. 10.1093/hmg/ddi209.

    Article  CAS  PubMed  Google Scholar 

  7. Johansson Moller M, Chaudhary R, Hellmen E, Hoyheim B, Chowdhary B, Andersson L: Pigs with the dominant white coat color phenotype carry a duplication of the KIT gene encoding the mast/stem cell growth factor receptor. Mamm Genome. 1996, 7: 822-830. 10.1007/s003359900244.

    Article  CAS  PubMed  Google Scholar 

  8. Hirooka H, de Koning DJ, van Arendonk JA, Harlizius B, de Groot PN, Bovenhuis H: Genome scan reveals new coat color loci in exotic pig cross. J Hered. 2002, 93: 1-8. 10.1093/jhered/93.1.1.

    Article  CAS  PubMed  Google Scholar 

  9. Marklund S, Kijas J, Rodriguez-Martinez H, Ronnstrand L, Funa K, Moller M, Lange D, Edfors-Lilja I, Andersson L: Molecular basis for the dominant white phenotype in the domestic pig. Genome Res. 1998, 8: 826-833.

    PubMed Central  CAS  PubMed  Google Scholar 

  10. Johansson M, Ellegren H, Marklund L, Gustavsson U, Ringmar-Cederberg E, Andersson K, Edfors-Lilja I, Andersson L: The gene for dominant white color in the pig is closely linked to ALB and PDGRFRA on chromosome 8. Genomics. 1992, 14: 965-969. 10.1016/S0888-7543(05)80118-1.

    Article  CAS  PubMed  Google Scholar 

  11. Giuffra E, Evans G, Tornsten A, Wales R, Day A, Looft H, Plastow G, Andersson L: The Belt mutation in pigs is an allele at the Dominant White (I/KIT) locus. Mamm Genome. 1999, 10: 1132-1136. 10.1007/s003359901178.

    Article  CAS  PubMed  Google Scholar 

  12. Pielberg G, Olsson C, Syvanen AC, Andersson L: Unexpectedly high allelic diversity at the KIT locus causing dominant white color in the domestic pig. Genetics. 2002, 160: 305-311.

    PubMed Central  CAS  PubMed  Google Scholar 

  13. Everitt B: Cluster Analysis. 1974, London: Heinemann Educational Books

    Google Scholar 

  14. SAS Institute Inc: SAS/STAT 9.1 user's guide. 2004, Cary, N.C.: SAS Institute, 6: 1377-1427.

    Google Scholar 

  15. Kar P, Polipalli SK, Chattopadhyay S, Hussain Z, Malik A, Husain SA, Medhi S, Begum N: Prevalence of hepatitis B virus genotype D in precore mutants among chronic liver disease patients from New Delhi, India. Dig Dis Sci. 2007, 52: 565-569. 10.1007/s10620-006-9625-z.

    Article  PubMed  Google Scholar 

  16. Kasehagen LJ, Mueller I, McNamara DT, Bockarie MJ, Kiniboro B, Rare L, Lorry K, Kastens W, Reeder JC, Kazura JW, Zimmerman PA: Changing patterns of Plasmodium blood-stage infections in the Wosera region of Papua New Guinea monitored by light microscopy and high throughput PCR diagnosis. Am J Trop Med Hyg. 2006, 75: 588-596.

    PubMed Central  PubMed  Google Scholar 

  17. Michel CE, Sonnex C, Carne CA, White JA, Magbanua JP, Nadala EC, Lee HH: Chlamydia trachomatis Load at Matched Anatomic Sites: Implications for Screening Strategies. J Clin Microbiol. 2007, 45: 1395-1402. 10.1128/JCM.00100-07.

    Article  PubMed Central  PubMed  Google Scholar 

  18. Milan D, Jeon JT, Looft C, Amarger V, Robic A, Thelander M, Rogel-Gaillard C, Paul S, Iannuccelli N, Rask L, Ronne H, Lundstrom K, Reinsch N, Gellin J, Kalm E, Roy PL, Chardon P, Andersson L: A mutation in PRKAG3 associated with excess glycogen content in pig skeletal muscle. Science. 2000, 288: 1248-1251. 10.1126/science.288.5469.1248.

    Article  CAS  PubMed  Google Scholar 

  19. Di Fiore F, Blanchard F, Charbonnier F, Le Pessot F, Lamy A, Galais MP, Bastit L, Killian A, Sesboue R, Tuech JJ, Queuniet AM, Paillot B, Sabourin JC, Michot F, Michel P, Frebourg T: Clinical relevance of KRAS mutation detection in metastatic colorectal cancer treated by Cetuximab plus chemotherapy. Br J Cancer. 2007, 96: 1166-1169. 10.1038/sj.bjc.6603685.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  20. Miller SA, Dykes DD, Polesky HF: A simple salting out procedure for extracting DNA from human nucleated cells. Nucleic Acids Res. 1988, 16: 1215-10.1093/nar/16.3.1215.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  21. Walsh PS, Metzger DA, Higuchi R: Chelex 100 as a medium for simple extraction of DNA for PCR-based typing from forensic material. Biotechniques. 1991, 10: 506-513.

    CAS  PubMed  Google Scholar 

  22. Yerle M, Echard G, Robic A, Mairal A, Dubut-Fontana C, Riquet J, Pinton P, Milan D, Lahbib-Mansais Y, Gellin J: A somatic cell hybrid panel for pig regional gene mapping characterized by molecular cytogenetics. Cytogenet Cell Genet. 1996, 73: 194-202.

    Article  CAS  PubMed  Google Scholar 

  23. Website title. []

Download references


We thank Dr. Martin Yerle (INRA, France) for supplying the somatic cell hybrid panel DNAs. This work was supported by a grant (no. 20050301034467) from the BioGreen 21 Program, Rural Development Administration, Korea. Bo-Young Seo, Hyun-Tae Im and Jae-Hwan Kim were supported by scholarships from the BK21 Program, Ministry of Education and Human Resources Development, Korea.

Author information

Authors and Affiliations


Corresponding author

Correspondence to Jin-Tae Jeon.

Additional information

Authors' contributions

JTJ conceived and designed the experiment, and drafted the manuscript. JHK did somatic cell hybrid mapping. BYS, EWP and HTI performed OLA and pyroseqencing. ICC, IKK and SHL collected blood and hair samples, and purified DNA. SJA and JHL performed the statistical analysis. All authors read and approved the final manuscript.

Electronic supplementary material


Additional file 1: Verification of the specificity of PCR primers used for amplifying breakpoints of KIT CNV using a porcine × rodent somatic cell hybrid panel. (a) The results from the primer set (KITBPF and KIT1BPR) for the normal copy. (b) The results from the primer set (KITBPF and KIT2BPR) for the duplicated copy. M, 100-bp size standard; numbers are positive clone numbers. (JPEG 66 KB)


Additional file 2: Standard curves for qOLA_CNV and Pyro_CNV. A serial dilution from 0% to 100% duplicated copy vs normal copy (PCR-amplified and cloned) was used for the estimation. (a) A curve estimated using peak height values from qOLA_CNV (correlation coefficient = 0.999). (b) A curve estimated using peak area values from qOLA_CNV (correlation coefficient = 0.999). (c) A curve estimated using Pyro_CNV (correlation coefficient = 0.995). (d) A comparison of the three curves by root mean square (RMS) of the bias to reference values and standard deviations. The qOLA_CNV using peak height measurements is the most accurate and precise of the three. (JPEG 126 KB)


Additional file 3: Typical coat colors of parental and F1 animals. (a) A Korean native boar (Black). (b) A Landrace sow (White). (c) An F1 littermate produced by a cross between a Korean native boar (i/i) and a Landrace sow (I1/i); four white pigs were genotyped as I1/i and three colored ones as i/i. (JPEG 191 KB)


Additional file 4: Comparison of genotyping results for 159 Large White pigs by the two genotyping methods. The genotypes for the clustering measurements on the plot are in the first column, and the numbers of class centroid for the statistical analysis are given in the first row. The two discrepancies between the assignment methods are indicated by italic and bold numbers. (DOC 38 KB)

Authors’ original submitted files for images

Rights and permissions

Open Access This article is published under license to BioMed Central Ltd. This is an Open Access article is distributed under the terms of the Creative Commons Attribution License ( ), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Reprints and permissions

About this article

Cite this article

Seo, BY., Park, EW., Ahn, SJ. et al. An accurate method for quantifying and analyzing copy number variation in porcine KIT by an oligonucleotide ligation assay. BMC Genet 8, 81 (2007).

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: