8p23 beta-defensin copy number determination by single-locus pseudogene-based paralog ratio tests risk bias due to low-frequency sequence variations

Development and Evaluation of a Droplet Digital PCR Assay for 8p23 β-Defensin Cluster Copy Number Determination

Introduction

It has been proposed that the copy number (CN) variation (CNV) in β-defensin genes (DEFB) on human chromosome 8p23 determines phenotypic differences in inflammatory diseases. However, no method for accurate and easy DEFB CN quantification is yet available.

Objective

Droplet digital polymerase chain reaction (ddPCR) is a novel method for CNV quantification and has been used for genes such as CCL4L, CCL3L1, AMY1, and HER2. However, to date, no ddPCR assay has been available for DEFB CN determination. In the present study, we aimed to develop and evaluate such a ddPCR assay.

Methods

The assay was designed using DEFB4 and RPP30 as the target and the reference gene, respectively. To evaluate the assay, 283 DNA samples with known CNs previously determined using the multiple ligation-dependent probe amplification (MLPA) method, the current gold standard, were used as standards. To discover the optimal DNA template amount, we tested 80 to 2.5 ng DNA by a serial of one to two dilutions of eight samples. To evaluate the reproducibility of the assay, 31 samples were repeated to calculate the intra- and inter-assay variations. To further validate the reliability of the assay, the CNs of all 283 samples were determined using ddPCR. To compare results with those using quantitative PCR (qPCR), DEFB CNs for 48 samples were determined using qPCR with the same primers and probes.

Results

In a one-dimensional plot, the positive and negative droplets were clearly separated in both DEFB4 and RPP30 detection channels. In a two-dimensional plot, four populations of droplets were observed. The 20 ng template DNA proved optimal, with either high (80 ng) or low (10, 5, 2.5 ng) template input leading to ambiguous or inaccurate results. For the 31 standard samples, DEFB CNs were accurately determined with small intra- and inter-assay variations (coefficient of variation < 0.04 for both). In the validation cohort, ddPCR provided the correct CN for all 283 samples with high confidence. qPCR measurements for the 48 samples produced noisy data with high uncertainty and low accuracy.

Conclusions

ddPCR is an accurate, reproducible, easy-to-use, cheap, high-throughput method for DEFB CN determination. ddPCR could be applied to DEFB CN quantification in large-scale case–control studies.

Supplementary Information

The online version contains supplementary material available at 10.1007/s40291-021-00546-2.

Analysis of Copy Number Variation Using the Paralogue Ratio Test (PRT)

Copy number variation (CNV), where a segment of DNA differs in copy number between different individuals, is an extensive and often underappreciated source of genetic variation within species. However, reliably determining copy number of a particular DNA sequence for a large number of samples can be challenging. Here, I describe and review the paralogue ratio test (PRT) in detail. PRT was developed to robustly type the CNV of the beta-defensin locus using small amounts of genomic DNA in a high-throughput manner, and has been applied successfully at many other loci. I discuss the strategies for designing successful PRT assays using both manual and bioinformatics methods, how to optimize experimental conditions, and approaches for analyzing the data. I discuss strengths and weaknesses of the approach, and how to troubleshoot results, as well as the range of problems to which PRT can be a potential solution.

Duplication of chicken defensin7 gene generated by gene conversion and homologous recombination

Defensins constitute an evolutionary conserved family of cationic antimicrobial peptides that play a key role in host innate immune responses to infection. Defensin genes generally reside in complex genomic regions that are prone to structural variation, and defensin genes exhibit extensive copy number variation in humans and in other species. Copy number variation of defensin genes was examined in inbred lines of Leghorn and Fayoumi chickens, and a duplication of defensin7 was discovered in the Fayoumi breed. Analysis of junction sequences confirmed the occurrence of a simple tandem duplication of defensin7 with sequence identity at the junction, suggesting nonallelic homologous recombination between defensin7 and defensin6 The duplication event generated two chimeric promoters that are best explained by gene conversion followed by homologous recombination. Expression of defensin7 was not elevated in animals with two genes despite both genes being transcribed in the tissues examined. Computational prediction of promoter regions revealed the presence of several putative transcription factor binding sites generated by the duplication event. These data provide insight into the evolution and possible function of large gene families and specifically, the defensins.

High Copy Numbers of β-Defensin Cluster on 8p23.1, Confer Genetic Susceptibility, and Modulate the Physical Course of Hidradenitis Suppurativa/Acne Inversa

Hidradenitis suppurativa/acne inversa (HS) has a multifactorial pathogenesis, with many patients reporting positive family history. Nine β-defensin genes (among them DEFB4 and DEFB103, encoding for proinflammatory mediators human β-defensin-2 and human β-defensin-3, respectively) exist as a cluster (DEFB) affected by copy number (CN). We hypothesized that CNs are greater in patients with HS and that they are linked to genetic susceptibility. CNs of DEFB were studied in two independent patient cohorts: 163 patients from Greece and 98 from Germany. CNs were greater in patients than control subjects in both studied cohorts. Carriage of more than six CNs was associated with a 7.53 odds ratio for HS in the Greek cohort and a 5.76 odds ratio for HS in the German cohort. The common odds ratio after meta-analysis was 6.72 (P < 0.0001). However, presence of fewer than six copies was linked with disease onset at an earlier age (P = 0.048), less frequent presentation of permanent purulence of the affected skin lesions (P = 0.036), and fewer skin localizations (P = 0.042). A robust genetic trait for susceptibility to HS is provided, and this is confirmed in two independent cohorts. Susceptibility arises from carriage of more than six DEFB copies, which interferes directly with the HS phenotype.

An MLPA-based approach for high-resolution genotyping of disease-related multi-allelic CNVs

Copy number variation has recently been recognized as an important type of genetic variation that modifies human phenotypes. Copy number variants (CNVs) are being increasingly associated with various human phenotypes and diseases. However, the lack of an appropriate method that allows fast, inexpensive and, most importantly, accurate CNVs genotyping significantly hampers CNV analysis. This limitation especially affects the analysis of multi-allelic CNVs that frequently modify various phenotypes. Recently, we developed a multiplex ligation-dependent probe amplification (MLPA)-based strategy for multiplex copy number genotyping and the validation of candidate CNV-miRNAs. Here we present the adaptation and optimization of this recently developed method for high-resolution genotyping of individual disease-related multi-allelic CNVs. We developed appropriate assays for three well-known and extensively studied CNVs: CNV-CCL3L1, CNV-DEFB, and CNV-UGT2B17, which have been associated with various human phenotypes including inflammation-related and infectious diseases. With the use of these assays we identified several general factors that allow to increase the resolution of the copy number genotyping. Performed experiments confirmed the high reproducibility and accuracy of the obtained genotyping results. The reliability of the results and relatively low per-genotype cost makes this strategy an attractive method for large-scale experiments such as genotype-phenotype association studies.