Technical considerations for genotyping multi-allelic copy number variation (CNV), in regions of segmental duplication
© Cantsilieris et al.; licensee BioMed Central Ltd. 2014
Received: 14 October 2013
Accepted: 22 April 2014
Published: 1 May 2014
Intrachromosomal segmental duplications provide the substrate for non-allelic homologous recombination, facilitating extensive copy number variation in the human genome. Many multi-copy gene families are embedded within genomic regions with high levels of sequence identity (>95%) and therefore pose considerable analytical challenges. In some cases, the complexity involved in analyzing such regions is largely underestimated. Rapid, cost effective analysis of multi-copy gene regions have typically implemented quantitative approaches, however quantitative data are not an absolute means of certainty. Therefore any technique prone to degrees of measurement error can produce ambiguous results that may lead to spurious associations with complex disease.
In this study we have focused on testing the accuracy and reproducibility of quantitative analysis techniques. With reference to the C-C Chemokine Ligand-3-like-1 (CCL3L1) gene, we performed analysis using real-time Quantitative PCR (QPCR), Multiplex Ligation-dependent Probe Amplification (MLPA) and Paralogue Ratio Test (PRT). After controlling for potential outside variables on assay performance, including DNA concentration, quality, preparation and storage conditions, we find that real-time QPCR produces data that does not cluster tightly around copy number integer values, with variation substantially greater than that of the MLPA or PRT systems. We find that the method of rounding real-time QPCR measurements can potentially lead to mis-scoring of copy number genotypes and suggest caution should be exercised in interpreting QPCR data.
We conclude that real-time QPCR is inherently prone to measurement error, even under conditions that would seem favorable for association studies. Our results indicate that potential variability in the physicochemical properties of the DNA samples cannot solely explain the poor performance exhibited by the real-time QPCR systems. We recommend that more robust approaches such as PRT or MLPA should be used to genotype multi-allelic copy number variation in disease association studies and suggest several approaches which can be implemented to ensure the quality of the copy number typing using quantitative methods.
KeywordsCopy number variation Multi-allelic CCL3L1 Real-time quantitative PCR Multiplex ligation-dependent probe amplification Paralogue ratio test
Structurally complex regions of the human genome contain a wealth of information that is yet to be fully understood within the context of evolution and disease. Here we use the term “complex” to describe genomic regions containing highly duplicated sequence of >95% sequence identity (eg. segmental duplications) demonstrating a range of diploid copies (>3 copy number classes). The C-C Chemokine Ligand-3-like-1 (CCL3L1) gene is one such example .
The quantitative methods commonly used to genotype multi-allelic copy number polymorphisms (CNPs) are not absolute, therefore accuracy and reproducibly are required to generate robust disease associations . Single or multi locus genotyping methods including Multiplex Ligation-dependent Probe Amplification (MLPA), Paralogue Ratio Test (PRT) and real-time Quantitative PCR (QPCR) are flexible, and cost effective ways to analyze multi-allelic CNPs . Although studies have shown that in some cases these techniques can provide an accurate assessment across a range of diploid copy numbers, their reproducibility can be influenced by a number of factors including assay design, sample quality and operator error.
A challenging aspect of multi-allelic copy number analysis is that the relative difference between copy number classes becomes smaller and more difficult to discriminate as copy number increases . Despite this challenge, CNP analysis should produce data that clusters tightly around integer values , unless there is evidence for tissue mosaicism in which different cell types contain different copy numbers . As such, tissue mosaicism can occur over a continuous spectrum (1-99%), and can result in extensive overlap between integer classes resulting in data that seemingly indicates a continuous copy number distribution. However, similar data can also be produced when measurements of non-mosaic tissues are imprecise and unable to discriminate CNVs with a high level of accuracy, potentially leading to erroneous conclusions.
Indeed, differences in these observations made in a number of studies support the conclusion that such measurement error can complicate data interpretation  and these effects are clearly evident in studies of the CCL3L1 gene and disease susceptibility [1, 8, 9]. In some cases, there could be misconceptions surrounding assay validation and/or data interpretation, and if this is coupled with a lack of transparency regarding the distribution of raw unrounded copy number measurements, it may lead to spurious results that are difficult to explain. The CCL3L1 gene represents a particularly intriguing locus for study, most notably for its implication in HIV AIDS susceptibility , but also because this region demonstrates several layers of complexity. In addition to being copy number variable, this locus is highly population stratified and the presence of extreme levels of sequence duplication can complicate assay design. As such due to the presence of the CC3L3 pseudogene, the analysis of certain CCL3L1 sequences differs with respect to copy number content and should be considered when attempting to compare methods and studies.
In this study we aim to demonstrate some simple approaches for assessing the quality of genotyping data when using methods such as real-time QPCR, PRT and MLPA for the analysis of multi-allelic copy number. These methods are evaluated by examining the copy number distribution of the CCL3L1 gene in DNA samples of European ethnicity.
We selected 150 DNA samples from patients of European origin. Genomic DNA was isolated from venous blood leukocytes using a phenol/chloroform extraction procedure as previously described . The concentration and purity of the DNA samples was determined using the NanoDrop Specrophotometer (Thermo Scientific). We considered the genomic DNA samples to be free of contaminants or solvents if the 260/280 ratio resulted in a range between 1.7-2.1 and the 260/230 ratio was >1.5. The samples were collected as part of ongoing eye studies at the Royal Victorian Eye and Ear Hospital (RVEEH), Melbourne. Subjects were given a standard risk factor questionnaire and following clinical examination are considered healthy control individuals. The study was conducted in accordance with the Declaration of Helsinki and according to the National Health and Medical Research Council of Australia’s statement on ethical conduct in research involving humans, revised in 2000. Written informed consent was obtained from all individuals, and ethics approval for the project was provided by the Human Research and Ethics Committee of the Royal Victorian Eye and Ear Hospital (RVEEH), Melbourne.
Multiplex Ligation-dependent Probe Amplification (MLPA)
MLPA probe sequences used for targeting the CCL3L1 gene and two reference loci EP300 and CREBBP
Real-time quantitative PCR (QPCR)
Real-time Quantitative PCR was performed according to the TaqMan® Copy Number Assays protocol (Applied Biosystems). We used 20 ng/μL of genomic DNA in a 10 μL reaction, consisting of 5 μL of 2× TaqMan® Genotyping Master Mix, 0.5 μL of RNaseP TaqMan assay 20× working stock, 0.625 μL of CCL3L1 Forward 10 μM, 0.625 μL of CCL3L1 Reverse 10 μM, 0.3 μL of CCL3L1 MGB probe 10 μM and 1 μL DH2O. The CCL3L1 MGB probe sequence is as follows MGB-6FA: TTCGAGGCCCAGCGACCTCA. The PCR primer sequences are as follows CCL3L1 Forward: GGGTCCAGAAATACGTCAGT and CCL3L1 Reverse: CATGTTCCCAAGGCTCAG as previously described . Reactions were performed in duplicate, the cycling times were 1 cycle of 95°C for 10 minutes, and 40 cycles of 95°C for 15 seconds and 60°C for 1 minute. To assess the PCR efficiencies of the test and reference locus, we performed a series of 10-fold serial dilutions for both the CCL3L1 and RNaseP PCR reactions. We found that the PCR efficiencies were approximately equal, differing by less than 1%. We exported the CT values for both and analyzed them using CopyCaller™ software from Applied Biosystems. CopyCaller™ performs a comparative CT (ΔΔCT) relative quantitation analysis whereby the difference between the threshold cycles of the target and reference gene is calculated, the ΔCT of the test sample is compared to a calibrator of known copy number. We chose the calibrator sample which represented two copies of the CCL3L1 gene as tested by MLPA.
Paralogue Ratio Test (PRT)
We performed Single Plex Paralogue Ratio Test (PRT) using the method described previously . We used 20 ng/μL of genomic DNA in a 25 μL reaction, comprising of 5 μL of 5× Bioline My Taq Buffer, 2 μL of 10 pmol/μL CCL3C Forward: GGC TAA GAC CCC TTC TAG AG and CCL3C Reverse: AAT CAT GCA GGT CTC CAC T Primer, 1 μL Dimethyl sulfoxide (DMSO), 0.25 μL Bioline My Taq Polymerase and 14.75 μL DH2O. The cycling times were 1 cycle of 95°C for 1 minute, 24 cycles of 95°C for 30 seconds, 55°C for 30 seconds and 70°C for 1 minute, followed by 1 cycle of 70°C for 40 minutes. The PCR primers were FAM labeled and for each PCR reaction, 1 μl was mixed with 8.9 μl of HiDi Formamide and 0.1 μl of ROX500 size standard (Applied Biosystems). PCR product separation was performed on the ABI 3130 Electrophoresis 16 capillary sequencer (Applied Biosystems). Data Analysis was performed using the GeneMapper Software (Applied Biosystems) and the peak heights were exported to Microsoft Excel. Copy number was calculated as the ratio of the CCL3L1 peak height against the reference peak height of the diploid CCL3 gene.
Validation of CNV genotyping methods
In order to validate each of the three genotyping methodologies, we compared 12 samples previously identified as having 1, 2, 3 or 4 CCL3L1 copies (three samples of each copy number class) using the MLPA assay, and compared these with results generated by both the PRT and real-time QPCR systems. We found that the genotype calls between all three methodologies were highly concordant, demonstrating that each assay represented data specific to each exon.
Assessment of assay performance and separation of integer copy number
CNV genotyping calls from three independent techniques, QPCR, PRT and MLPA, when analyzing the CCL3L1 gene
Measurement error has the power to drastically influence disease association studies, and this is particularly evident when performing analysis of multi-allelic CNPs. Studies have shown that even a small degree of differential bias, applied systematically to large case/control data sets, is likely to result in spurious association between multi-allelic CNPs and disease [2, 5] Our results demonstrate quite conclusively that real-time QPCR is prone to measurement error and performs poorly when applied to the CCL3L1 multi-allelic CNP region, when compared to alternative techniques such as MLPA or PRT.
The performance of real-time QPCR has previously been assessed for genotyping multi-allelic CNP regions such as the CCL3L1 gene  and the Beta-defensin gene cluster [2, 15, 16]. In each case, two main conclusions concerning the results have been reported 1. There is a lack of clear integer clustering 2. There are discrepancies between some copy number calls when compared to alternative approaches. Both of these findings are supported by the data presented here. Previous studies have attempted to compare the performance of CCL3L1 genotyping using the real-time QPCR and PRT systems . However as the assays targeted different exons, the presence of a pseudogene excluded a more direct comparison of the techniques. Additionally, the cases and controls used in the study by Field et al. were extracted using different techniques potentially introducing variation in template quality as a contributor to the poor performance of real-time QPCR.
Several studies have suggested that the biological properties of the DNA including protein content, shearing, methods of extraction and storage conditions are potential contributors to the poor performance of real-time QPCR for analyzing multi-allelic CNPs [2, 14, 15, 17]. Additionally, it has been suggested that real-time QPCR will perform comparably to PRT if the concentration of DNA samples is normalized to a uniform concentration, and the efficiency between the test and reference loci differ by less than 5% . With knowledge of these potential variables on assay performance, we felt that it would be unfair to perform a method comparison on DNA samples that were subject to such limitations. Therefore, we reasoned that if DNA quality is crucial to the success of real-time QPCR, then using high quality input DNA samples of similar concentration, generated using a single technique, and stored under uniform conditions, should produce high quality data. However, the results we present here do not support the conclusion that sample quality or PCR efficiency are the only reason for the variability observed with the real-time QPCR system, and thus we present an alternative hypothesis.
Reported number of technical replicates used in real-time QPCR studies for disease association
Technical replicates performed
Gonzalez et al. 
HIV related Outcomes
Urban et al. 
HIV related Outcomes
Bhattacharya et al. 
HIV related Outcomes
Field et al. 
Mamtani et al. 
Systemic Lupus Erythematosus
Fellerman et al. 
Bentley et al. 
Mehlotra et al. 
HIV related Outcomes
Previous studies have highlighted the critical importance of implementing copy number reference samples, plotting the raw unrounded copy number measurement data and observing the clustering around integer values in disease association studies [2, 5, 25]. We also demonstrate that assay validation should not consist solely in genotyping positive controls as real-time QPCR performed comparably to PRT and MLPA when a small number of positive controls were genotyped. Our data shows that because the variance for real-time QPCR measurements surrounding integer values is substantially greater than that of the PRT or MLPA systems, only by plotting and ordering the raw unrounded measurements from a larger sample size on a histogram can the short-comings of QPCR be identified. It is when a smaller sample size is utilized that the results can potentially be misleading. We would like to point out that this can be done by genotyping as few as 40 test samples, well before an entire cohort is genotyped. Therefore, we believe the use of positive controls should be strictly limited to assessing the ability of an assay to detect copy number variation, and then calibrating inter-run variation.
In some cases, studies using real-time QPCR have systematically removed uncertain data that fails to cluster around integer values in an attempt to reduce the impact of measurement error [26, 27]. However, calling samples with only a high degree of certainty may result in bias, as the failure to cluster is unlikely to be independent of genotype . A more robust approach using real-time QPCR would be to firstly; perform a sufficient number of replicates as dictated by the copy number content of the given region; secondly repeat samples in which the inference of the most likely copy number has extended error bars, utilize the z-scores to highlight data with low confidence ; and thirdly, implement a statistical model such as CNVtools which specifically controls for differential bias in the case/control analysis while also allowing one to input data as probability weighted CNV scores as well as integers .
One of the more troubling aspects regarding the use of real-time QPCR in our study is the failure to form integer clusters on a sample set of high quality DNA of similar concentration. This suggests that real-time QPCR is inherently prone to measurement error. Moreover this error may be exacerbated by varying sample quality in large case/control cohorts potentially leading to false positive associations. However, it is also noteworthy that real-time QPCR is not substantially incorrect in its assessment of copy number genotype, as demonstrated in Table 2. Our data shows that even though under optimal conditions the variance of real-time QPCR is considerable across large sample numbers, however it still may have a place as purely a confirmatory assay for sequencing projects or MLPA and PRT data sets.
The consistent level of clustering observed by both the PRT and MLPA systems performed without replicates or repeat testing is unequivocal evidence of their accuracy and reliability on a large scale. Additionally, we have not fully utilized the multiplexing capacity of MLPA  or Triplex PRT  which would further improve the copy number typing of these systems. One important aspect to consider is that combining measurements from two independent techniques, such as PRT and MLPA, is very powerful with respect to accuracy and reproducibility. We found that MLPA and PRT calls differed by a single copy in five samples, notably in the 4–5 copy range. However, by combining the two techniques using a concordance plot we were able to unambiguously assign copy number to these samples (Figure 4A).
We argue that many multi-allelic CNP regions have been largely underestimated with regard to the complexity of analysis, as evidenced by the large number of contradictory evidence concerning disease association studies (reviewed in ). There also appears to be a fundamental need to understand how these techniques work, along with which variables can impact the reproducibility of the assay. We also highlight that real-time QPCR appears to be satisfactory for detecting differences between 1 and 2 copies, thus may be appropriate for the detection of di-allelic CNPs where the differences between the integers are at least 2-fold. Recently, water in oil digital droplet PCR (ddPCR) has been described for genotyping multi-allelic loci with an impressive level of reproducibility and accuracy [30, 31]. While this technology has not yet been widely implemented for association studies, it is likely that ddPCR will resolve a number of issues regarding measurement error in large cohorts. The development of a low cost genotyping assay that can provide accurate copy number content, and also possesses the ability to discriminate sequence content, including smaller insertion/deletions and paralogue specific variation, is required to fully understand the complex nature of multi-copy genomic regions. Recent work combining molecular inversion probes (MIPs) together with massively parallel DNA sequencing  appears perfectly suited for this type of analysis and will bridge the gap between locus specific quantitative measures and whole genome sequencing technologies. Although MIPs have immense promise for accurately genotyping multi-allelic CNPs such as CCL3L1, it is noteworthy that MIP technology is heavily reliant on high quality reference sequences and availability of singly unique nucleotides (SUNs) to accurately distinguish gene family paralogs. Therefore at present MIP- based genotyping cannot be universally applied to all multi-allelic CNP regions.
In conclusion the data presented here clearly demonstrates that real-time QPCR is prone to measurement error when applied to genotyping the CCL3L1 multi-allelic CNP region. We recommend that real-time QPCR data be treated with a high degree of caution and suggest that if it is used for the analysis of multi-allelic CNPs, that it is accompanied by independent genotyping measures such as MLPA or PRT. With reference to quantitative measures, ddPCR technology, MIPs or MLPA and PRT, represent far more reliable and robust approaches for accurately resolving multi-allelic CNPs and should be used in preference if possible.
Availability of supporting data
The data supporting the results of this article are included within the article.
Paul N Baird and Stefan J White are recognized as joint senior authors.
This work is supported by funding from Monash University awarded to PW and SW, and the National Health and Medical Research Council Centre for Clinical Research Excellence Grant Number #529923. Translational Clinical Research in Major Eye Diseases. NHMRC CJ Martin Biomedical Fellowship awarded to SC; Grant Number #1073726. NHMRC Senior Research Fellowship awarded to PNB Grant Number #1028444. MIMR-PHI Institute and The Centre for Eye Research Australia (CERA) receive operational infrastructure support from the Victorian Government.
- Gonzalez E, Kulkarni H, Bolivar H, Mangano A, Sanchez R, Catano G, Nibbs RJ, Freedman BI, Quinones MP, Bamshad MJ, Murthy KK, Rovin BH, Bradley W, Clark RA, Anderson SA, O'Connell RJ, Agan BK, Ahuja SS, Bologna R, Sen L, Dolan MJ, Ahuja SK: The influence of CCL3L1 gene-containing segmental duplications on HIV-1/AIDS susceptibility. Science. 2005, 307 (5714): 1434-1440. 10.1126/science.1101160.PubMedView ArticleGoogle Scholar
- Aldhous M, Abu Bakar S, Prescott N, Palla R, Soo K, Mansfield J, Mathew C, Satsangi J, Armour J: Measurement methods and accuracy in copy number variation: failure to replicate associations of beta-defensin copy number with Crohn’s disease. Hum Mol Genet. 2010, 19 (24): 4930-4938. 10.1093/hmg/ddq411.PubMed CentralPubMedView ArticleGoogle Scholar
- Cantsilieris S, Baird PN, White SJ: Molecular methods for genotyping complex copy number polymorphisms. Genomics. 2012, 101 (2): 86-93.PubMedView ArticleGoogle Scholar
- Hollox EJ: The challenges of studying complex and dynamic regions of the human genome. Methods Mol Biol. 2012, 838: 187-207. 10.1007/978-1-61779-507-7_9.PubMedView ArticleGoogle Scholar
- Carpenter D, Walker S, Prescott N, Schalkwijk J, Armour J: Accuracy and differential bias in copy number measurement of CCL3L1 in association studies with three auto-immune disorders. BMC Genomics. 2011, 12 (1): 418-10.1186/1471-2164-12-418.PubMed CentralPubMedView ArticleGoogle Scholar
- Notini AJ, Craig JM, White SJ: Copy number variation and mosaicism. Cytogenet Genome Res. 2009, 123 (1–4): 270-277.Google Scholar
- Cantsilieris S, White SJ: Correlating multi-allelic copy number polymorphisms with disease susceptibility. Hum Mutat. 2013, 34 (1): 1-13.PubMedView ArticleGoogle Scholar
- Bhattacharya T, Stanton J, Kim EY, Kunstman KJ, Phair JP, Jacobson LP, Wolinsky SM: CCL3L1 and HIV/AIDS susceptibility. Nat Med. 2009, 15 (10): 1112-1115. 10.1038/nm1009-1112.PubMedView ArticleGoogle Scholar
- Urban TJ, Weintrob AC, Fellay J, Colombo S, Shianna KV, Gumbs C, Rotger M, Pelak K, Dang KK, Detels R, Martinson JJ: CCL3L1 and HIV/AIDS susceptibility. Nat Med. 2009, 15 (10): 1110-1112. 10.1038/nm1009-1110.PubMed CentralPubMedView ArticleGoogle Scholar
- Sambrook J, Fritsch EF, Maniatis T: Molecular cloning: a Laboratory manual. Cold Spring Harbor Laboratory Press. 1989, 1: E3-E7.Google Scholar
- White SJ, Vink GR, Kriek M, Wuyts W, Schouten J, Bakker B, Breuning MH, den Dunnen JT: Two-color multiplex ligation-dependent probe amplification: detecting genomic rearrangements in hereditary multiple exostoses. Hum Mutat. 2004, 24 (1): 86-92. 10.1002/humu.20054.PubMedView ArticleGoogle Scholar
- White SJ, Breuning MH, den Dunnen JT: Detecting copy number changes in genomic DNA: MAPH and MLPA. Methods Mol Biol. 2004, 75: 751-768.Google Scholar
- Sudmant P, Kitzman JO, Antonacci F, Alkan C, Malig M, Tsalenko A, Sampas N, Bruhn L, Shendure J, Eichler EE: 1000 Genomes Project. Diversity of human copy number variation and multicopy genes. Science. 2010, 330 (6004): 641-646. 10.1126/science.1197005.PubMed CentralPubMedView ArticleGoogle Scholar
- Field SF, Howson JMM, Maier LM, Walker S, Walker NM, Smyth DJ, Armour JAL, Clayton DG, Todd JA: Experimental aspects of copy number variant assays at CCL3L1. Nat Med. 2009, 15 (10): 1115-1117. 10.1038/nm1009-1115.PubMed CentralPubMedView ArticleGoogle Scholar
- Fode P, Jespersgaard C, Hardwick R, Bogle H, Theisen M, Dodoo D, Lenicek M, Vitek L, Vieira A, Freitas J, Andersen PS, Hollox EJ: Determination of beta-defensin genomic copy number in different populations: a comparison of three methods. PLoS One. 2011, 6 (2): e16768-10.1371/journal.pone.0016768.PubMed CentralPubMedView ArticleGoogle Scholar
- Perne A, Zhang X, Lehmann LE, Groth M, Stuber F, Book M: Comparison of Multiplex Ligation-dependant Probe Amplification and real time PCR accuracy for gene copy number quantification using the Beta Defensin Locus. Biotechniques. 2009, 47 (6): 1023-1026. 10.2144/000113300.PubMedView ArticleGoogle Scholar
- Cukier HN, Pericak-Vance MA, Gilbert JR, Hedges DJ: Sample degradation leads to false-positive copy number variation calls in multiplex real-time polymerase chain reaction assays. Anal Biochem. 2009, 386 (2): 288-290. 10.1016/j.ab.2008.11.040.PubMedView ArticleGoogle Scholar
- Fernandez-Jimenez N, Castellanos-Rubio A, Plaza-Izurieta L, Gutierrez G, Irastorza I, Castan L, Vitoria J, Bilbao J: Accuracy in copy number calling by qPCR and PRT: a matter of DNA. PLoS One. 2011, 6 (12): e28910-10.1371/journal.pone.0028910.PubMed CentralPubMedView ArticleGoogle Scholar
- Weaver S, Dube S, Mir A, Qin J, Gang S, Ramakrishnan R, Jones RC, Livak K: Taking qPCR to a higher level: analysis of CNV reveals the power of high throughpout qPCR to enhance quantitative resolution. Methods. 2010, 50: 271-276. 10.1016/j.ymeth.2010.01.003.PubMedView ArticleGoogle Scholar
- Hollox EJ, Armour JA, Barber JCK: Extensive normal copy number variation of a β-defensin antimicrobial-gene cluster. A J Hum Genet. 2003, 73 (3): 591-600. 10.1086/378157.View ArticleGoogle Scholar
- Mamtani M, Rovin B, Brey R, Camargo JF, Kulkarni H, Herrera M, Correa P, Holliday S, Anaya JM, Ahuja SK: CCL3L1 gene-containing segmental duplications and polymorphisms in CCR5 affect risk of systemic lupus erythaematosus. Ann Rheum Dis. 2008, 67 (8): 1076-1083. 10.1136/ard.2007.078048.PubMed CentralPubMedView ArticleGoogle Scholar
- Fellermann K, Stange DE, Schaeffeler E, Schmalzl H, Wehkamp J, Bevins CL, Reinisch W, Teml A, Schwab M, Lichter P, Radlwimmer B, Stange EF: A chromosome 8 gene-cluster polymorphism with Low human beta-defensin 2 gene copy number predisposes to Crohn disease of the colon. Am J Hum Genet. 2006, 79 (3): 439-448. 10.1086/505915.PubMed CentralPubMedView ArticleGoogle Scholar
- Bentley RW, Pearson J, Gearry RB, Barclay ML, McKinney C, Merriman TR, Roberts RL: Association of higher DEFB4 genomic copy number with Crohn’s disease. Am J Gastroenterol. 2010, 105 (2): 354-359. 10.1038/ajg.2009.582.PubMedView ArticleGoogle Scholar
- Mehlotra R, Dazard J, John B, Zimmerman P, Weinberg A, Jurevic R: Copy number variation within human β-defensin gene cluster influences progression to AIDS in the multicenter AIDS cohort study. J Aids Clin Res. 2012, 3 (10): 1000184-PubMed CentralPubMedView ArticleGoogle Scholar
- Hollox EJ: Beta-Defensins and Crohn’s disease: confusion from counting copies. Am J Gastroenterol. 2010, 105 (2): 360-362. 10.1038/ajg.2009.573.PubMedView ArticleGoogle Scholar
- McKinney C, Broen JCA, Vonk MC, Beretta L, Hesselstrand R, Hunzelmann N, Riemekasten G, Scorza R, Simeon CP, Fonollosa V, Carreira PE, Ortego-Centeno N, Gonzalez-Gay MA, Airo P, Coenen M, Martin J, Radstake TRDJ, Merriman TR: Evidence that deletion at FCGR3B is a risk factor for systemic sclerosis. Genes Immun. 2012, 13 (6): 458-460. 10.1038/gene.2012.15.PubMedView ArticleGoogle Scholar
- McKinney C, Merriman ME, Chapman PT, Gow PJ, Harrison AA, Highton J, Jones PBB, McLean L, O’Donnell JL, Pokorny V, Spellerberg M, Stamp LK, Willis J, Steer S, Merriman TR: Evidence for an influence of chemokine ligand 3-like 1 (CCL3L1) gene copy number on susceptibility to rheumatoid arthritis. Ann Rheum Dis. 2008, 67 (3): 409-413.PubMedView ArticleGoogle Scholar
- Barnes C, Plagnol V, Fitzgerald T, Redon R, Marchini J, Clayton D, Hurles ME: A robust statistical method for case–control association testing with copy number variation. Nat Genet. 2008, 40 (10): 1245-1252. 10.1038/ng.206.PubMed CentralPubMedView ArticleGoogle Scholar
- D’haene B, Vandesompele J, Hellemans J: Accurate and objective copy number profiling using real-time quantitative PCR. Methods. 2010, 50: 262-270. 10.1016/j.ymeth.2009.12.007.PubMedView ArticleGoogle Scholar
- Boettger L, Handsaker R, Zody M, McCarroll S: Structural haplotypes and recent evolution of the human 17q21.31 region. Nat Genet. 2012, 44 (8): 881-886. 10.1038/ng.2334.PubMed CentralPubMedView ArticleGoogle Scholar
- Hindson BJ, Ness KD, Masquelier DA, Belgrader P, Heredia NJ, Makarewicz AJ, Bright IJ, Lucero MY, Hiddessen AL, Legler TC, Kitano TK, Hodel MR, Petersen JF, Wyatt PW, Steenblock ER, Shah PH, Bousse LJ, Troup CB, Mellen JC, Wittmann DK, Erndt NG, Cauley TH, Koehler RT, So AP, Dube S, Rose KA, Montesclaros L, Wang S, Stumbo DP, Hodges SP, et al: High-throughput droplet digital PCR system for absolute quantitation of DNA copy number. Anal Chem. 2011, 83 (22): 8604-8610. 10.1021/ac202028g.PubMed CentralPubMedView ArticleGoogle Scholar
- Nuttle X, Huddleston J, O’Roak BJ, Antonacci F, Fichera M, Romano C, Shendure J, Eichler EE: Rapid and accurate large-scale genotyping of duplicated genes and discovery of interlocus gene conversions. Nat Meth. 2013, 10 (9): 903-909. 10.1038/nmeth.2572.View ArticleGoogle Scholar
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly credited. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.