A worldwide analysis of beta-defensin copy number variation suggests recent selection of a high-expressing DEFB103 gene copy in East Asia

Beta-defensins are a family of multifunctional genes with roles in defense against pathogens, reproduction, and pigmentation. In humans, six beta-defensin genes are clustered in a repeated region which is copy-number variable (CNV) as a block, with a diploid copy number between 1 and 12. The role in host defense makes the evolutionary history of this CNV particularly interesting, because morbidity due to infectious disease is likely to have been an important selective force in human evolution, and to have varied between geographical locations. Here, we show CNV of the beta-defensin region in chimpanzees, and identify a beta-defensin block in the human lineage that contains rapidly evolving noncoding regulatory sequences. We also show that variation at one of these rapidly evolving sequences affects expression levels and cytokine responsiveness of DEFB103, a key inhibitor of influenza virus fusion at the cell surface. A worldwide analysis of beta-defensin CNV in 67 populations shows an unusually high frequency of high-DEFB103-expressing copies in East Asia, the geographical origin of historical and modern influenza epidemics, possibly as a result of selection for increased resistance to influenza in this region. Hum Mutat 32:743–750, 2011. © 2011 Wiley-Liss, Inc.

Determination of beta-defensin genomic copy number in different populations: a comparison of three methods

Background

There have been conflicting reports in the literature on association of gene copy number with disease, including CCL3L1 and HIV susceptibility, and β-defensins and Crohn's disease. Quantification of precise gene copy numbers is important in order to define any association of gene copy number with disease. At present, real-time quantitative PCR (QPCR) is the most commonly used method to determine gene copy number, however the Paralogue Ratio Test (PRT) is being used in more and more laboratories.

Findings

In this study we compare a Pyrosequencing-based Paralogue Ratio Test (PPRT) for determining beta-defensin gene copy number with two currently used methods for gene copy number determination, QPCR and triplex PRT by typing five different cohorts (UK, Danish, Portuguese, Ghanaian and Czech) of DNA from a total of 576 healthy individuals. We found a systematic measurement bias between DNA cohorts revealed by QPCR, but not by the PRT-based methods. Using PRT, copy number ranged from 2 to 9 copies, with a modal copy number of 4 in all populations.

Conclusions

QPCR is very sensitive to quality of the template DNA, generating systematic biases that could produce false-positive or negative disease associations. Both triplex PRT and PPRT do not show this systematic bias, and type copy number within the correct range, although triplex PRT appears to be a more precise and accurate method to type beta-defensin copy number.

Assessment of complement C4 gene copy number using the paralog ratio test

The complement C4 locus is in the class III region of the MHC, and exhibits copy number variation. Complement C4 null alleles have shown association with a number of diseases including systemic lupus erythematosus (SLE). However, most studies to date have used protein immunophenotyping and not direct interrogation of the genome to determine C4 null allele status. Moreover, a lack of accurate C4 gene copy number (GCN) estimation and tight linkage disequilibrium across the disease-associated MHC haplotypes has confounded attempts to establish whether or not these associations are causal. We have therefore developed a high throughput paralog ratio test (PRT) in association with two restriction enzyme digest variant ratio tests (REDVRs) to determine total C4 GCN, C4A GCN, and C4B GCN. In the densely genotyped CEU cohort we show that this method is accurate and reproducible when compared to gold standard Southern blot copy number estimation with a discrepancy rate of 9%. We find a broad range of C4 GCNs in the CEU and the 1958 British Birth Cohort populations under study. In addition, SNP-C4 CNV analyses show only moderate levels of correlation and therefore do not support the use of SNP genotypes as proxies for complement C4 GCN.

Beta-defensins and Crohn's disease: confusion from counting copies

Recent research has suggested that a family of antimicrobial and inflammatory molecules, called beta-defensins, are involved in the etiology of Crohn's disease. In this issue, Bentley and colleagues provide data that disagree with previous studies. They show that the beta-defensins that are copy number polymorphic have, on average, higher genomic copy number in patients when compared with healthy controls. This editorial places these new data in the context of previous research on the genetics of beta-defensin copy number polymorphism and association with Crohn's disease. It also suggests a path allowing this research area to move forward with confidence.

Constitutional trisomy 8 and Behçet syndrome

The characteristic clinical features of constitutional trisomy 8 include varying degrees of developmental delay, joint contractures and deep palmar and plantar creases. There is an established literature, which describes features of Behçet syndrome occurring in phenotypically normal individuals with myelodysplastic syndromes and trisomy 8 in their bone marrow. In this article, we describe four patients with constitutional trisomy 8, all with varying clinical phenotypes, who developed features of Behçet, in particular but not exclusively mucocutaneous ulceration. In addition, we examined gene copy numbers of the variable-number neutrophil defensin genes DEFA1A3 in one of the cases (case 1) and her parents, together with 14 cases of Behçet syndrome in comparison with 121 normal controls. The gene copy number was highest in case 1 (copy number 14) and was also increased in her parents (both copy number 9). However the mean copy number for DEFA1A3 among the 14 Behçet syndrome patients was actually lower (5.1) than among the controls (mean of 6.8 copies). Thus, we conclude that patients with constitutional trisomy 8 and those with trisomy 8 confined to the bone marrow are both at increased risk of developing features of Behçet syndrome. The mechanism may relate to increased chromosome 8 gene dosage with further analysis of candidate genes on chromosome 8 required.

An integrated approach for measuring copy number variation at the FCGR3 (CD16) locus

Copy number variation (CNV) is an important source of genomic diversity in humans, and influences disease susceptibility. The immunoglobulin-receptor genes FCGR3A and FCGR3B on chromosome 1q23.3 show CNV, and CNV of the FCGR3B gene is associated with glomerulonephritis in systemic lupus erythematosus and organ-specific autoimmunity. Large-scale case-control association studies of CNV require technologies that are amenable to high-throughput analysis with low error rates. Here we propose an integrated suite of five assays, four of them duplexed to reduce DNA usage, that assays for CNV at FCGR3A and FCGR3B, and genotype the polymorphic neutrophil antigen HNA1. We show how a maximum-likelihood (ML) approach to combining the results from these five assays allows estimation of statistical confidence for each individual copy number, and therefore an appropriate significance threshold to be set, controlling the error rate. This approach results in a high-throughput copy number genotyping system, with demonstrable precision and accuracy, that can be applied to large case-control cohort studies. We demonstrate Mendelian inheritance of this CNV, variation in frequency between Europeans and East Asians, and a lack of strong association between the CNV and flanking SNP genotypes, with important consequences for genome-wide association studies.

Beta-defensin-2 protein is a serum biomarker for disease activity in psoriasis and reaches biologically relevant concentrations in lesional skin

BACKGROUND

Previous studies have extensively documented antimicrobial and chemotactic activities of beta-defensins. Human beta-defensin-2 (hBD-2) is strongly expressed in lesional psoriatic epidermis, and recently we have shown that high beta-defensin genomic copy number is associated with psoriasis susceptibility. It is not known, however, if biologically and pathophysiologically relevant concentrations of hBD-2 protein are present in vivo, which could support an antimicrobial and proinflammatory role of beta-defensins in lesional psoriatic epidermis.

METHODOLOGY/PRINCIPAL FINDINGS

We found that systemic levels of hBD-2 showed a weak but significant correlation with beta defensin copy number in healthy controls but not in psoriasis patients with active disease. In psoriasis patients but not in atopic dermatitis patients, we found high systemic hBD-2 levels that strongly correlated with disease activity as assessed by the PASI score. Our findings suggest that systemic levels in psoriasis are largely determined by secretion from involved skin and not by genomic copy number. Modelling of the in vivo epidermal hBD-2 concentration based on the secretion rate in a reconstructed skin model for psoriatic epidermis provides evidence that epidermal hBD-2 levels in vivo are probably well above the concentrations required for in vitro antimicrobial and chemokine-like effects.

CONCLUSIONS/SIGNIFICANCE

Serum hBD-2 appears to be a useful surrogate marker for disease activity in psoriasis. The discrepancy between hBD-2 levels in psoriasis and atopic dermatitis could explain the well known differences in infection rate between these two diseases.

Defensins and the dynamic genome: what we can learn from structural variation at human chromosome band 8p23.1

Over the past four years, genome-wide studies have uncovered numerous examples of structural variation in the human genome. This includes structural variation that changes copy number, such as deletion and duplication, and structural variation that does not change copy number, such as orientation and positional polymorphism. One region that contains all these types of variation spans the chromosome band 8p23.1. This region has been studied in some depth, and the focus of this review is to examine our current understanding of the variation of this region. We also consider whether this region is a good model for other structurally variable regions in the genome and what the implications of this variation are for clinical studies. Finally, we discuss the bioinformatics challenges raised, discuss the evolution of the region, and suggest some future priorities for structural variation research.

Accurate, high-throughput typing of copy number variation using paralogue ratios from dispersed repeats

Recent work has demonstrated an unexpected prevalence of copy number variation in the human genome, and has highlighted the part this variation may play in predisposition to common phenotypes. Some important genes vary in number over a high range (e.g. DEFB4, which commonly varies between two and seven copies), and have posed formidable technical challenges for accurate copy number typing, so that there are no simple, cheap, high-throughput approaches suitable for large-scale screening. We have developed a simple comparative PCR method based on dispersed repeat sequences, using a single pair of precisely designed primers to amplify products simultaneously from both test and reference loci, which are subsequently distinguished and quantified via internal sequence differences. We have validated the method for the measurement of copy number at DEFB4 by comparison of results from >800 DNA samples with copy number measurements by MAPH/REDVR, MLPA and array-CGH. The new Paralogue Ratio Test (PRT) method can require as little as 10 ng genomic DNA, appears to be comparable in accuracy to the other methods, and for the first time provides a rapid, simple and inexpensive method for copy number analysis, suitable for application to typing thousands of samples in large case-control association studies.

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

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.

Duplications and copy number variants of 8p23.1 are cytogenetically indistinguishable but distinct at the molecular level

It has been proposed that duplications of 8p23.1 are either euchromatic variants of the 8p23.1 defensin domain with no phenotypic consequences or true duplications associated with developmental delay and heart defects. Here, we provide evidence for both alternatives in two new families. A duplication of most of band 8p23.1 (circa 5 Mb) was found in a girl of 8 years with pulmonary stenosis and mild language delay. BAC fluorescence in situ hybridisation (FISH) and multiplex amplifiable probe hybridisation (MAPH) showed that the two copies of the duplicated segment were sited, in an alternating fashion, between three copies of a circa 300-450 kb segment from 8p23.1 distal to REPD. Copy number of the variable 8p23.1 defensin domain was consistent with duplication but within the normal range. Duplication of the GATA-binding protein 4 gene (GATA4) in this patient and others with and without heart defects, suggests it is a dosage-sensitive gene with variable penetrance. A cytogenetically similar duplication of 8p23.1 was found at prenatal diagnosis in a fetus, father and grandmother. There was no duplication using BAC FISH but MAPH showed 11 copies of the 360 kb variable defensin domain which is within the expanded range found in previous euchromatic variant carriers. Semiquantitative FISH (SQ-FISH) was consistent with a simultaneous expansion of the adjacent olfactory receptor repeats. These results distinguish duplications of 8p23.1 with clinically significant consequences from benign copy number variants, which have not yet been associated with qualitative or quantitative traits.

Copy number polymorphism and expression level variation of the human alpha-defensin genes DEFA1 and DEFA3

We have defined unexpectedly extensive copy number variation at the human anti-microbial alpha-defensin genes DEFA1 and DEFA3, encoding human neutrophil peptides HNP-1, HNP-2 and HNP-3. There was variation in both number and position of DEFA1/DEFA3 genes in arrays of 19 kb tandem repeats on 8p23.1, so that the DEFA1 and DEFA3 genes appear to be interchangeable variant cassettes within tandem gene arrays. For this reason, the official symbol for this locus has been revised to DEFA1A3. The total number of gene copies per diploid genome varied between four and 11 in a sample of 111 control individuals from the UK, with approximately 10% (11/111) of people lacking DEFA3 completely. DEFA1 appeared to be at high copy number in all great apes studied; at one variable site in the repeat unit, both variants have persisted in humans, chimpanzees and gorillas since their divergence. Analysis of expression levels in human white blood cells showed a clear correlation between the relative proportions of DEFA1:DEFA3 mRNA and corresponding gene numbers. However, there was no relationship between total (DEFA1+DEFA3) mRNA levels and total gene copy number, suggesting the superimposed influence of trans-acting factors. The persistence of DEFA1 at high copy number in other apes suggests an alternative model for the early stages of the evolution of novel genes by duplication and divergence. Duplicated genes present in variant tandem arrays may have greater potential than simple duplications for the combinatorial creation of new functions by recombination and gene conversion, while still preserving pre-existing functions on the same haplotype.

A 4q35.2 subtelomeric deletion identified in a screen of patients with co-morbid psychiatric illness and mental retardation

Background

Cryptic structural abnormalities within the subtelomeric regions of chromosomes have been the focus of much recent research because of their discovery in a percentage of people with mental retardation (UK terminology: learning disability). These studies focused on subjects (largely children) with various severities of intellectual impairment with or without additional physical clinical features such as dysmorphisms. However it is well established that prevalence of schizophrenia is around three times greater in those with mild mental retardation. The rates of bipolar disorder and major depressive disorder have also been reported as increased in people with mental retardation. We describe here a screen for telomeric abnormalities in a cohort of 69 patients in which mental retardation co-exists with severe psychiatric illness.

Methods

We have applied two techniques, subtelomeric fluorescence in situ hybridisation (FISH) and multiplex amplifiable probe hybridisation (MAPH) to detect abnormalities in the patient group.

Results

A subtelomeric deletion was discovered involving loss of 4q in a patient with co-morbid schizoaffective disorder and mental retardation.

Conclusion

The precise region of loss has been defined allowing us to identify genes that may contribute to the clinical phenotype through hemizygosity. Interestingly, the region of 4q loss exactly matches that linked to bipolar affective disorder in a large multiply affected Australian kindred.

Extensive normal copy number variation of a beta-defensin antimicrobial-gene cluster

Using a combination of multiplex amplifiable probe hybridization and semiquantitative fluorescence in situ hybridization (SQ-FISH), we analyzed DNA copy number variation across chromosome band 8p23.1, a region that is frequently involved in chromosomal rearrangements. We show that a cluster of at least three antimicrobial beta-defensin genes (DEFB4, DEFB103, and DEFB104) at 8p23.1 are polymorphic in copy number, with a repeat unit >/=240 kb long. Individuals have 2-12 copies of this repeat per diploid genome. By segregation, microsatellite dosage, and SQ-FISH chromosomal signal intensity ratio analyses, we deduce that individual chromosomes can have one to eight copies of this repeat unit. Chromosomes with seven or eight copies of this repeat unit are identifiable by cytogenetic analysis as a previously described 8p23.1 euchromatic variant. Analysis of RNA from different individuals by semiquantitative reverse-transcriptase polymerase chain reaction shows a significant correlation between genomic copy number of DEFB4 and levels of its messenger RNA (mRNA) transcript. The peptides encoded by these genes are potent antimicrobial agents, especially effective against clinically important pathogens, such as Pseudomonas aeruginosa and Staphylococcus aureus, and DEFB4 has been shown to act as a cytokine linking the innate and adaptive immune responses. Therefore, a copy number polymorphism involving these genes, which is reflected in mRNA expression levels, is likely to have important consequences for immune system function.

High throughput screening of human subtelomeric DNA for copy number changes using multiplex amplifiable probe hybridisation (MAPH)

BACKGROUND

Subtelomeric regions of the human genome are gene rich, with a high level of sequence polymorphism. A number of clinical conditions, including learning disability, have been attributed to subtelomeric deletions or duplications, but screening for deletion in these regions using conventional cytogenetic methods and fluorescence in situ hybridisation (FISH) is laborious. Here we report that a new method, multiplex amplifiable probe hybridisation (MAPH), can be used to screen for copy number at subtelomeric regions.

METHODS

We have constructed a set of MAPH probes with each subtelomeric region represented at least once, so that one gel lane can assay copy number at all chromosome ends in one person. Each probe has been sequenced and, where possible, its position relative to the telomere determined by comparison with mapped clones.

RESULTS

The sensitivity of the probes has been characterised on a series of cytogenetically verified positive controls and 83 normal controls were used to assess the frequency of polymorphic copy number with no apparent phenotypic effect. We have also used MAPH to test a cohort of 37 people selected from males referred for fragile X syndrome testing and found six changes that were confirmed by dosage PCR.

CONCLUSIONS

MAPH can be used to screen subtelomeric regions of chromosomes for deletions and duplications before confirmation by FISH or dosage PCR. The high throughput nature of this technique allows it to be used for large scale screening of subtelomeric copy number, before confirmation by FISH. In practice, the availability of a rapid and efficient screen may allow subtelomeric analysis to be applied to a wider selection of patients than is currently possible using FISH alone.

DNA copy number analysis by MAPH: molecular diagnostic applications

DNA copy number variation is an important cause of genetic disease. There are several techniques available to detect copy number changes of various sizes, each with their limitations in resolution and cost. Here we outline the development of multiplex amplifiable probe hybridization (MAPH) into a high-throughput diagnostic technique for detecting copy number variation of almost any size. Its application in testing for genetic mutations causing diseases, such as familial breast cancer, Charcot-Marie-Tooth disease Type 1A, Duchenne/Becker muscular dystrophy and familial colorectal cancer is described, as well as its use in identifying chromosomal changes in some individuals with mental retardation. The analysis of the data produced by MAPH is also considered, along with its potential for automation and development of microarray-based MAPH.

Intolerance to lactose and other dietary sugars

Intolerance of dietary carbohydrate and sugars can result from a variety of genetically determined enzyme and transporter deficiencies. This article reviews this topic and discusses in more detail the current state of our own research on lactase.

Lactase haplotype diversity in the Old World

Lactase persistence, the genetic trait in which intestinal lactase activity persists at childhood levels into adulthood, varies in frequency in different human populations, being most frequent in northern Europeans and certain African and Arabian nomadic tribes, who have a history of drinking fresh milk. Selection is likely to have played an important role in establishing these different frequencies since the development of agricultural pastoralism approximately 9,000 years ago. We have previously shown that the element responsible for the lactase persistence/nonpersistence polymorphism in humans is cis-acting to the lactase gene and that lactase persistence is associated, in Europeans, with the most common 70-kb lactase haplotype, A. We report here a study of the 11-site haplotype in 1,338 chromosomes from 11 populations that differ in lactase persistence frequency. Our data show that haplotype diversity was generated both by point mutations and recombinations. The four globally common haplotypes (A, B, C, and U) are not closely related and have different distributions; the A haplotype is at high frequencies only in northern Europeans, where lactase persistence is common; and the U haplotype is virtually absent from Indo-European populations. Much more diversity is seen in sub-Saharan Africans than in non-Africans, consistent with an "Out of Africa" model for peopling of the Old World. Analysis of recent recombinant haplotypes by allele-specific PCR, along with deduction of the root haplotype from chimpanzee sequence, allowed construction of a haplotype network that assisted in evaluation of the relative roles of drift and selection in establishing the haplotype frequencies in the different populations. We suggest that genetic drift was important in shaping the general pattern of non-African haplotype diversity, with recent directional selection in northern Europeans for the haplotype associated with lactase persistence.

Common polymorphism in a highly variable region upstream of the human lactase gene affects DNA-protein interactions

In most mammals lactase activity declines after weaning when lactose is no longer part of the diet, but in many humans lactase activity persists into adult life. The difference responsible for this phenotypic polymorphism has been shown to be cis-acting to the lactase gene. The causal sequence difference has not been found so far, but a number of polymorphic sites have been found within and near to the lactase gene. We have shown previously that in Europeans there are two polymorphic sites in a small region between 974 bp and 852 bp upstream from the start of transcription, which are detectable by denaturing gradient gel electrophoresis (DGGE). In this study, analysis of individuals from five other population groups by the same DGGE method reveals four new alleles resulting from three additional nucleotide changes within this very small region. Analysis of sequence in four primate species and comparison with the published pig sequence shows that the overall sequence of this highly variable human region is conserved in pigs as well as primates, and that it lies within a 1kb region which has been shown to control lactase downregulation in pigs. Electrophoretic mobility shift assay (EMSA) studies were carried out to determine whether common variation affected protein-DNA binding and several binding activities were found using this technique. A novel two base-pair deletion that is common in most populations tested, but is not present in Europeans, caused no change in binding activity. However, a previously published C to T transition at -958bp dramatically reduced binding activity, although the functional significance of this is not clear.

The genetically programmed down-regulation of lactase in children

BACKGROUND & AIMS

Intestinal lactase activity is high in all healthy human babies, but in adults a genetic polymorphism, which acts in cis to the lactase gene, determines high or low messenger RNA (mRNA) expression and activity (lactase persistence and nonpersistence, respectively). Our aim was to investigate the onset of expression of this polymorphism in children.

METHODS

Activities were analyzed in relation to age in normal biopsy specimens from a 20-year collection of diagnostic specimens. In a smaller set of 32 samples, aged 2-132 months, RNA was extracted for semiquantitative reverse-transcription polymerase chain reaction. Marker polymorphisms were used to determine the allelic origin of lactase mRNA transcripts.

RESULTS

Analysis of 866 children showed evidence that the lactase persistence/nonpersistence polymorphism began before 5 years of age. The 32 children tested had high lactase mRNA and activity. Six children aged 2-16 months showed equal expression of two alleles, 2 children aged 7 and 14 months showed slightly asymmetric expression, and 7 children aged 22-132 months showed very asymmetric expression, the second allele being undetectable in the 11-year-old, as previously seen in lactase-persistent heterozygote adults.

CONCLUSIONS

Genetically programmed down-regulation of the lactase gene is detectable in children from the second year of life, although the onset and extent are somewhat variable.

Lactase haplotype frequencies in Caucasians: association with the lactase persistence/non-persistence polymorphism

A genetic polymorphism is responsible for determining that some humans express lactase at high levels throughout their lives and are thus lactose tolerant, while others lose lactase expression during childhood and are lactose intolerant. We have previously shown that this polymorphism is controlled by an element or elements which act in cis to the lactase gene. We have also reported that 7 polymorphisms in the lactase gene are highly associated and lead to only 3 common haplotypes (A, B and C) in individuals of European extraction. Here we report the frequencies of these polymorphisms in Caucasians from north and south Europe and also from the Indian sub-continent, and show that the alleles differ in frequency, the B and C haplotypes being much more common in southern Europe and India. Allelic association studies with lactase persistence and non-persistence phenotypes show suggestive evidence of association of lactase persistence with certain alleles. This association was rather more clear in the analysis of small families, where haplotypes could be determined. Furthermore haplotype and RNA transcript analysis of 11 unrelated lactase persistent individuals shows that the persistence (highly expressed) allele is almost always on the A haplotype background. Non-persistence is found on a variety of haplotypes including A. Thus it appears that lactase persistence arose more recently than the DNA marker polymorphisms used here to define the main Caucasian haplotypes, possibly as a single mutation on the A haplotype background. The high frequency of the A haplotype in northern Europeans is consistent with the high frequency of lactase persistence.