Skip to main content

Beta-defensin genomic copy number is not a modifier locus for cystic fibrosis

Abstract

Human beta-defensin 2 (DEFB4, also known as DEFB2 or hBD-2) is a salt-sensitive antimicrobial protein that is expressed in lung epithelia. Previous work has shown that it is encoded in a cluster of beta-defensin genes at 8p23.1, which varies in copy number between 2 and 12 in different individuals. We determined the copy number of this locus in 355 patients with cystic fibrosis (CF), and tested for correlation between beta-defensin cluster genomic copy number and lung disease associated with CF. No significant association was found.

Peer Review reports

Background

Cystic fibrosis (CF) is an autosomal recessive genetic disease in epithelia (Online Mendelian Inheritance in Man (OMIM) 219700, http://www.ncbi.nlm.nih.gov), with chronic lung infection being the major cause of morbidity and mortality. It has a a carrier frequency of about 1 in 25 in the United Kingdom, and the disease is due to mutations in the cystic fibrosis transmembrane conductance regulator gene (CFTR) gene, which encodes a membrane chloride channel (OMIM 602421). There is a wide spectrum of different mutations known to cause CF, but the most common is the ΔF508 mutation, which exists in the European population as a polymorphic allele with a frequency of between 2 and 3%.

The course of lung infection can vary greatly between patients with the same CFTR mutation, and despite the wide spectrum of disease causing mutations, only a few correlate with lung disease severity. The search for genetic loci modifying lung disease in CF patients has produced several candidates, with the gene for mannose binding lectin (MBL) being the most notable. MBL is a part of the innate immune system, and low MBL in serum has been associated with infections in both adults and children [1, 2]. Several studies have reported links with severity of CF lung disease, although the conclusion regarding the impact of the heterozygous deficient state differ amongst these [3–5].

Like MBL, beta-defensins are an important part of innate immunity. The peptide human beta-defensin 2 (encoded by the gene DEFB4, also known as DEFB2 or hBD-2) has been identified as a salt-sensitive antibacterial agent expressed in the lung airway [6]. Given the chronic bacterial colonisation of the lungs, this is a potential candidate for a CF modifier locus. It is possible that variation in DEFB4 levels could cause changes in the antibacterial effectiveness of airway surface fluid. Deterioration in lung function during clinical progress of CF is correlated with DEFB4 expression levels, suggesting again that DEFB4 has an important antibacterial role in airway surface fluid [7]. Additionally, DEFB4 expression is up-regulated by mucoid Pseudomonas aeruginosa, and has very effective bactericidal activity against it [8, 9], but another important pathogen in CF lungs, Burkholderia cepacia, is not killed by airway surface fluid expressing DEFB4 [10].

Three beta-defensins are expressed in significant levels in airway epithelia: human beta-defensin 1 (DEFB1) which is expressed in the lung [11], human beta-defensin 3 (DEFB103, also known as DEFB3 or hBD-3) which is expressed in the trachea [12], and DEFB4, also expressed in the lung [6]. All three genes are part of a beta-defensin cluster at chromosomal locus 8p23.1 [13]. This cluster of beta defensins, with the exception of DEFB1, has been shown to be on a repeat unit of at least 260 kb which varies in copy number in different people. Individuals can have between 2 and 12 copies, with most individuals having 3, 4 or 5 copies, and copy number is correlated with basal unstimulated expression level of DEFB4 in lymphoblastoid cell lines [14]. There have been several genome-wide studies confirming this copy number variation [15–17], and they highlight this type of large scale genomic polymorphism as an underappreciated source of phenotypic variation. Indeed, there is evidence that variable CCL3L1 gene copy number is associated with susceptibility to HIV infection [18].

The hypothesis that we wished to test in this study was that copy number variation of this beta-defensin cluster (and, by implication, DEFB4) is associated with lung disease in CF patients.

Results

We determined beta-defensin cluster copy number for 355 CF patients, of whom 343 had full clinical data. 159 patients were ΔF508 CFTR homozygotes, 141 were ΔF508 CFTR heterozygotes and 43 had non-ΔF508 mutations on both chromosomes. The mean age was 29.1 ± 8.9 years and 57% were male. Figure 1 shows the distribution of beta-defensin copy number phenotypes in the full CF patient cohort. There was no difference in the copy number phenotype distribution between this cohort and a random sample of 167 normal UK individuals (E.J.Hollox and J.A.L.Armour unpublished data).

Figure 1
figure 1

Beta-defensin diploid copy number distribution in the cohort of 355 cystic fibrosis patients analysed.

Each analysis was performed with the whole cohort, and with the ΔF508 CFTR homozygotes only. There were no significant correlations between beta-defensin copy number and each respiratory clinical parameter measured (mean and current FEV1, mean and current FVC). The patients were then arbitrarily divided into two groups based on copy number: one group comprising patients with 3 and fewer copies of the β-defensin region (n = 77), and another group with patients who have 4 copies and more (n = 266). Differences in lung function and inflammatory marker status were tested using two-tailed Mann-Whitney U-tests, differences in infection status were tested using χ2 tests. With this sample size, a difference of 9% in FEV1 median between the two groups could be detected with 80% power. No significant differences between any two groups were seen.

We also decided to determine the effectiveness of microsatellite dosage ratios in reinforcing, or possibly acting as a replacement for, MAPH in determining copy number. We analysed a subset of samples with copy numbers of three or four by determining the genotype of a microsatellite in the intron of DEFB4. The allelic distributions are shown in Figure 2. Out of 114 samples tested, 81 gave informative results, of which 71 (88%) agreed with the determination by MAPH. The samples that differed typically gave unrounded MAPH scores between the two values (range 3.42–3.54 or 4.5–4.54), and may be due either to error rate or, perhaps more likely given the small confidence intervals of each MAPH result (see Methods) and its reliability for other loci [19–23], genuine copy number heterogeneity between probes across the repeats. We should be able to distinguish between genuine copy number heterogeneity and MAPH error rate by analysing the copy number estimate given by the DEFB4 probe alone, which is 735 bp from the microsatellite in the intron of DEFB4, and should report the same copy number as the microsatellite. Analysis of this copy number estimate for the 10 discrepant results shows that it agrees with the microsatellite estimate for 5 samples. Therefore, the error rate in distinguishing 3 from 4 copies by MAPH is 5/81 = 6%. The result also suggests that copy number heterogeneity involving the DEFB4 probe (probe F) has a frequency of around 6%. In addition, copy number heterogeneity involving MAPH probe H has been found in several well-characterised individual DNAs (E.J.Hollox, J.A.L.Armour and J.C.K.Barber, unpublished data).

Figure 2
figure 2

Allelic distribution of DEFB4 microsatellite (EPEV-2) alleles in individuals with three and four copies of the beta-defensin gene cluster.

Discussion

This is the first study to test the hypothesis, raised in a previous paper [14], that variable β-defensin genomic copy number affects immune system function; in this case lung infection in a susceptible CF patient cohort. From the data we have presented, it is clear that either β-defensin genomic copy number has no effect on lung function in CF patients, or that any effect is too small for it to be detected using this cohort size. We cannot formally rule out an effect of single nucleotide polymorphisms (SNPs) within or surrounding the DEFB4 gene, given the apparent lack of association between copy number alleles and microsatellite alleles, and by inference any single nucleotide polymorphism (SNP). Any further study on these SNPs will have to carefully consider the copy number variation at this locus, given that a single nucleotide change could not only be allelic but differ between paralogues on the same chromosome [16].

There could be two reasons why β-defensin genomic copy number has no effect on lung function in CF patients. Firstly, unlike in cultured lymphocytes [14], there may be no correlation between DEFB4 copy number and gene expression in lung epithelia. Ideally a positive correlation between DEFB4 copy number and DEFB4 expression levels in non-inflamed lungs should be established, but realistically such studies are difficult to perform for both technical and ethical reasons. Even if such a direct link were proved, the lung in CF patients is characterised by a strong inflammatory response, and, given the dramatic increase in DEFB4 peptide levels in inflamed tissue [24], this could hide any differences in basal expression levels due to genomic copy number variation. Secondly, differences in DEFB4 expression are masked by other factors in the CF lung: for example, elastolytic cathepsins which breakdown DEFB4 and DEFB103 are present in the CF lung, but not in the healthy lung [25]. A recent report studying single nucleotide variation in the DEFB1 gene and lung function in CF patients has found no link, providing further evidence against beta-defensins modifying lung function in CF patients [26].

The key to further progress is to understand role of antimicrobial peptides in the inflamed lung, and exactly how they contribute to airway surface defence. A recent paper shows that DEFB4 is more than an antimicrobial peptide, and attracts neutrophils to sites of inflammation and infection [27], so it is possible that its chemoattractant role is as important in the lung as its antimicrobial role. How DEFB4 interacts with other components of the inflammatory response, and the innate and adaptive immune response, is only just beginning to be understood, and it is clear that much more work is needed before we understand these processes in lung disease.

Methods

Patient samples

DNA was extracted from peripheral blood from cystic fibrosis patients using the QIAamp DNA Blood Midi kit (Qiagen) according to the manufacturer's recommendations as part of a larger modifier gene study. Samples were coded and anonymised. All participants gave informed consent and the study was approved by the Royal Brompton Hospital and Harefield NHS Trust Ethics Committee.

Clinical parameters

Clinical data were obtained from the clinical databases, patients' hospital notes and computerised microbiology reports. Forced expiratory volume in 1 second (FEV1) and forced vital capacity (FVC) were obtained from annual lung function laboratory records, had been corrected for sex, age and height and are expressed as percent of predicted normal values. Annual values were recorded where available from 1990–2002 and annual rate of decline was calculated for all subjects with 2 or more values after linear regression of all available data. Transcutaneous oxygen saturation with the patient breathing room air was obtained at the annual visit, and a record of the requirement for supplementary oxygen was made. Additional data included the presence of infection with P. aeruginosa or B. cepacia and most recent total white blood cell count, C-reactive protein and liver function tests.

Multiplex amplifiable probe hybridisation (MAPH)

MAPH was carried out essentially as previously described [28, 29], with fluorescent detection of amplification products by an Applied Biosystems 3100 Genetic Analyzer. A probe set including six probes mapping to the beta-defensin cluster was constructed (Table 1). The values for these six probes were normalised first for lane intensity against control subtelomeric probes in the same probe set which do not show copy number variation [20]. These ratios were then calibrated against a weighted average of two control genomic samples in the same experiment, usually N8 (3 copies) and N25 (7 copies) to get a value for the number of copies of each probe per diploid cell [14]. In two instances, two other known controls were used on duplicate samples (European Collection of Animal Cell Cultures (ECACC) samples BO0183 3 copies – AF0105 4 copies). Of the 355 samples tested, 205 were tested in duplicate. Copy number was calculated by averaging the values of each probe mapping to this region and presenting the data both as a mean value +/- 95% confidence intervals (CIs), and as a rounded integer, representing the most likely copy number. Initially all six probes were used in the calculation, but it was found that DEFB104 probe was unusually sensitive to DNA quantity and removing the values for that probe in many cases improved the confidence intervals. All data in this paper are from five probes only. The 95% CIs for each copy number reading in the full cohort ranged from 0.05 to 1.2 with a median value of 0.24.

Table 1 MAPH probes that detect beta-defensin copy number variation

Microsatellite analysis

A microsatellite (EPEV-2) in the intron of DEFB4, was amplified using the following primers (EPEV-2F 5'-GCACCAGAGACCTCATGTTTTC-3' and EPEV-2R 5'-GTAACTTACAGTTGAAAACCAC-3'), with EPEV-2F 5' fluorescently-labelled with 6-FAM. The PCR conditions were as follows: 0.2 mM each of dATP, dCTP, dGTP, dTTP, 1 mM MgCl2, 75 mM Tris-HCl (pH 8.8 at 25°C), 20 mM (NH4)2SO4, 0.01% (v/v) Tween 20 in a final volume of 10 μl, with 5 pmol of each primer, 10 ng genomic DNA, 0.5 units Taq DNA polymerase (ABgene), followed by cycling for 25 cycles at 95°C 1 minute, 60°C 1 minute 72°C 1 minute, and a final extension incubation of 72°C 20 minutes. 1.5 μl of the PCR product was mixed with 10 μl Hi-Di formamide and run on an ABI 3100 using standard conditions, with G500 ROX marker. Allelic sizes ranged from 240 bp to 272 bp, and samples were omitted from analysis if stutter peaks made the interpretation of allele peak areas or allele lengths ambiguous.

References

  1. Summerfield JA, Ryder S, Sumiya M, Thursz M, Gorchein A, Monteil MA, Turner MW: Mannose binding protein gene mutations associated with unusual and severe infections in adults. Lancet. 1995, 345: 886-889. 10.1016/S0140-6736(95)90009-8.

    Article  CAS  PubMed  Google Scholar 

  2. Summerfield JA, Sumiya M, Levin M, Turner MW: Association of mutations in mannose binding protein gene with childhood infection in consecutive hospital series. Bmj. 1997, 314: 1229-1232.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  3. Gabolde M, Guilloud-Bataille M, Feingold J, Besmond C: Association of variant alleles of mannose binding lectin with severity of pulmonary disease in cystic fibrosis: cohort study. Bmj. 1999, 319: 1166-1167.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  4. Garred P, Pressler T, Madsen HO, Frederiksen B, Svejgaard A, Hoiby N, Schwartz M, Koch C: Association of mannose-binding lectin gene heterogeneity with severity of lung disease and survival in cystic fibrosis. J Clin Invest. 1999, 104: 431-437.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  5. Davies JC, Turner MW, Klein N: Impaired pulmonary status in cystic fibrosis adults with two mutated MBL-2 alleles. Eur Respir J. 2004, 24: 798-804. 10.1183/09031936.04.00055404.

    Article  CAS  PubMed  Google Scholar 

  6. Bals R, Wang X, Wu Z, Freeman T, Bafna V, Zasloff M, Wilson JM: Human beta-defensin 2 is a salt-sensitive peptide antibiotic expressed in human lung. J Clin Invest. 1998, 102: 874-880.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  7. Chen CI, Schaller-Bals S, Paul KP, Wahn U, Bals R: Beta-defensins and LL-37 in bronchoalveolar lavage fluid of patients with cystic fibrosis. J Cyst Fibros. 2004, 3: 45-50. 10.1016/j.jcf.2003.12.008.

    Article  CAS  PubMed  Google Scholar 

  8. Harder J, Meyer-Hoffert U, Teran LM, Schwichtenberg L, Bartels J, Maune S, Schroder JM: Mucoid Pseudomonas aeruginosa, TNF-alpha, and IL-1beta, but not IL-6, induce human beta-defensin-2 in respiratory epithelia. Am J Respir Cell Mol Biol. 2000, 22: 714-721.

    Article  CAS  PubMed  Google Scholar 

  9. Harder J, Bartels J, Christophers E, Schroder JM: A peptide antibiotic from human skin. Nature. 1997, 387: 861-10.1038/43088.

    Article  CAS  PubMed  Google Scholar 

  10. Baird RM, Brown H, Smith AW, Watson ML: Burkholderia cepacia is resistant to the antimicrobial activity of airway epithelial cells. Immunopharmacology. 1999, 44: 267-272. 10.1016/S0162-3109(99)00122-8.

    Article  CAS  PubMed  Google Scholar 

  11. Goldman MJ, Anderson GM, Stolzenberg ED, Kari UP, Zasloff M, Wilson JM: Human beta-defensin-1 is a salt-sensitive antibiotic in lung that is inactivated in cystic fibrosis. Cell. 1997, 88: 553-560. 10.1016/S0092-8674(00)81895-4.

    Article  CAS  PubMed  Google Scholar 

  12. Jia HP, Schutte BC, Schudy A, Linzmeier R, Guthmiller JM, Johnson GK, Tack BF, Mitros JP, Rosenthal A, Ganz T, McCray PBJ: Discovery of new human beta-defensins using a genomics-based approach. Gene. 2001, 263: 211-218. 10.1016/S0378-1119(00)00569-2.

    Article  CAS  PubMed  Google Scholar 

  13. Semple CA, Rolfe M, Dorin JR: Duplication and selection in the evolution of primate beta-defensin genes. Genome Biol. 2003, 4: R31-10.1186/gb-2003-4-5-r31.

    Article  PubMed  PubMed Central  Google Scholar 

  14. Hollox EJ, Armour JA, Barber JC: Extensive normal copy number variation of a beta-defensin antimicrobial-gene cluster. Am J Hum Genet. 2003, 73: 591-600. 10.1086/378157.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  15. Iafrate AJ, Feuk L, Rivera MN, Listewnik ML, Donahoe PK, Qi Y, Scherer SW, Lee C: Detection of large-scale variation in the human genome. Nat Genet. 2004, 36: 949-951. 10.1038/ng1416.

    Article  CAS  PubMed  Google Scholar 

  16. Fredman D, White SJ, Potter S, Eichler EE, Den Dunnen JT, Brookes AJ: Complex SNP-related sequence variation in segmental genome duplications. Nat Genet. 2004, 36: 861-866. 10.1038/ng1401.

    Article  CAS  PubMed  Google Scholar 

  17. Sebat J, Lakshmi B, Troge J, Alexander J, Young J, Lundin P, Maner S, Massa H, Walker M, Chi M, Navin N, Lucito R, Healy J, Hicks J, Ye K, Reiner A, Gilliam TC, Trask B, Patterson N, Zetterberg A, Wigler M: Large-scale copy number polymorphism in the human genome. Science. 2004, 305: 525-528. 10.1126/science.1098918.

    Article  CAS  PubMed  Google Scholar 

  18. 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 R J, 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: 1434-1440. 10.1126/science.1101160.

    Article  CAS  PubMed  Google Scholar 

  19. White S, Kalf M, Liu Q, Villerius M, Engelsma D, Kriek M, Vollebregt E, Bakker B, Van Ommen GJ, Breuning MH, Den Dunnen JT: Comprehensive Detection of Genomic Duplications and Deletions in the DMD Gene, by Use of Multiplex Amplifiable Probe Hybridization. Am J Hum Genet. 2002, 71: 365-374. 10.1086/341942.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  20. Hollox EJ, Atia T, Cross G, Parkin T, Armour JAL: High-throughput screening of human subtelomeric DNA for copy number changes using Multiplex Amplifiable Probe Hybridisation (MAPH). J Med Genet. 2002, 39: 790-795. 10.1136/jmg.39.11.790.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  21. Hollox EJ, Armour JA: No evidence for DNA copy number change associated with the DUP25 cytogenetic phenotype. Eur J Hum Genet. 2003, 11: 911-912. 10.1038/sj.ejhg.5201074.

    Article  CAS  PubMed  Google Scholar 

  22. Reid AG, Tarpey PS, Nacheva EP: High-resolution analysis of acquired genomic imbalances in bone marrow samples from chronic myeloid leukemia patients by use of multiple short DNA probes. Genes Chromosomes Cancer. 2003, 37: 282-290. 10.1002/gcc.10215.

    Article  CAS  PubMed  Google Scholar 

  23. 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 

  24. Singh PK, Jia HP, Wiles K, Hesselberth J, Liu L, Conway BA, Greenberg EP, Valore EV, Welsh MJ, Ganz T, Tack BF, McCray PBJ: Production of beta-defensins by human airway epithelia. Proc Natl Acad Sci U S A. 1998, 95: 14961-14966. 10.1073/pnas.95.25.14961.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  25. Taggart CC, Greene CM, Smith SG, Levine RL, McCray PBJ, O'Neill S, McElvaney NG: Inactivation of human beta-defensins 2 and 3 by elastolytic cathepsins. J Immunol. 2003, 171: 931-937.

    Article  CAS  PubMed  Google Scholar 

  26. Vankeerberghen A, Scudiero O, De Boeck K, Macek MJ, Pignatti PF, Van Hul N, Nuytten H, Salvatore F, Castaldo G, Zemkova D, Vavrova V, Cassiman JJ, Cuppens H: Distribution of human beta-defensin polymorphisms in various control and cystic fibrosis populations. Genomics. 2005, 85: 574-581. 10.1016/j.ygeno.2005.02.003.

    Article  CAS  PubMed  Google Scholar 

  27. Niyonsaba F, Ogawa H, Nagaoka I: Human beta-defensin-2 functions as a chemotactic agent for tumour necrosis factor-alpha-treated human neutrophils. Immunology. 2004, 111: 273-281. 10.1111/j.0019-2805.2004.01816.x.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  28. Armour JA, Sismani C, Patsalis PC, Cross G: Measurement of locus copy number by hybridisation with amplifiable probes. Nucleic Acids Res. 2000, 28: 605-609. 10.1093/nar/28.2.605.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  29. www.maph.info: .

Download references

Acknowledgements

We wish to thank Terence Low for DNA extraction, Ning Shen for help with clinical data collection, and Jess Tyson for helpful advice and comments. EJH is supported by a Wellcome Trust Bioarchaeology Postdoctoral Fellowship (grant no. 071024).

Author information

Authors and Affiliations

Authors

Corresponding author

Correspondence to Edward J Hollox.

Additional information

Authors' contributions

EJH carried out the molecular laboratory work, and participated in the design of the study. JD and EJH carried out the statistical analysis, and JD together with JB participated in clinical data collection and method design. UG, EA, and JALA participated in the design of the study. All authors helped to draft the manuscript, and have read and approved the final manuscript.

Authors’ original submitted files for images

Below are the links to the authors’ original submitted files for images.

Authors’ original file for figure 1

Authors’ original file for figure 2

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 ( https://creativecommons.org/licenses/by/2.0 ), 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

Hollox, E.J., Davies, J., Griesenbach, U. et al. Beta-defensin genomic copy number is not a modifier locus for cystic fibrosis. J Negat Results BioMed 4, 9 (2005). https://doi.org/10.1186/1477-5751-4-9

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/1477-5751-4-9

Keywords