Gene copy-number variation and associated polymorphisms of complement component C4 in human systemic lupus erythematosus (SLE): low copy number is a risk factor for and high copy number is a protective factor against SLE susceptibility in European Americans

Interindividual gene copy-number variation (CNV) of complement component C4 and its associated polymorphisms in gene size (long and short) and protein isotypes (C4A and C4B) probably lead to different susceptibilities to autoimmune disease. We investigated the C4 gene CNV in 1,241 European Americans, including patients with systemic lupus erythematosus (SLE), their first-degree relatives, and unrelated healthy subjects, by definitive genotyping and phenotyping techniques. The gene copy number (GCN) varied from 2 to 6 for total C4, from 0 to 5 for C4A, and from 0 to 4 for C4B. Four copies of total C4, two copies of C4A, and two copies of C4B were the most common GCN counts, but each constituted only between one-half and three-quarters of the study populations. Long C4 genes were strongly correlated with C4A (R=0.695; P<.0001). Short C4 genes were correlated with C4B (R=0.437; P<.0001). In comparison with healthy subjects, patients with SLE clearly had the GCN of total C4 and C4A shifting to the lower side. The risk of SLE disease susceptibility significantly increased among subjects with only two copies of total C4 (patients 9.3%; unrelated controls 1.5%; odds ratio [OR] = 6.514; P=.00002) but decreased in those with > or =5 copies of C4 (patients 5.79%; controls 12%; OR=0.466; P=.016). Both zero copies (OR=5.267; P=.001) and one copy (OR=1.613; P=.022) of C4A were risk factors for SLE, whereas > or =3 copies of C4A appeared to be protective (OR=0.574; P=.012). Family-based association tests suggested that a specific haplotype with a single short C4B in tight linkage disequilibrium with the -308A allele of TNFA was more likely to be transmitted to patients with SLE. This work demonstrates how gene CNV and its related polymorphisms are associated with the susceptibility to a human complex disease.     

Discovery of new human beta-defensins using a genomics-based approach

Epithelial beta-defensins are broad-spectrum cationic antimicrobial peptides that also act as chemokines for adaptive immune cells. In the human genome, all known defensin genes cluster to a <1 Mb region of chromosome 8p22-p23. To identify new defensin genes, the DNA sequence from a contig of large-insert genomic clones from the region containing human beta-defensin-2 (HBD-2) was analyzed for the presence of defensin genes. This sequence survey identified a novel beta-defensin, termed HBD-3. The HBD-3 gene contains two exons, is located 13 kb upstream from the HBD-2 gene, and it is transcribed in the same direction. A partial HBD-3 cDNA clone was amplified from cDNA derived from IL-1beta induced fetal lung tissue. The cDNA sequence encodes for a 67 amino acid peptide that is approximately 43% identical to HBD-2 and shares the beta-defensin six cysteine motif. By PCR analysis of two commercial cDNA panels, HBD-3 expression was detected in adult heart, skeletal muscle, placenta and in fetal thymus. From RT-PCR experiments, HBD-3 expression was observed in skin, esophagus, gingival keratinocytes, placenta and trachea. Furthermore, in fetal lung explants and gingival keratinocytes, HBD-3 mRNA expression was induced by IL-1beta. Additional sequence analysis identified the HE2 (human epididymis secretory protein) gene 17 kb upstream from the HBD-3 gene. One splice variant of this gene (HE2beta1) encodes a beta-defensin consensus cysteine motif, suggesting it represents a defensin gene product. HE2beta1 mRNA expression was detected in gingival keratinocytes and bronchial epithelia using RT-PCR analysis. The discovery of these novel beta-defensin genes may allow further understanding of the role of defensins in host immunity at mucosal surfaces.     

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.     

A 450-kb contig of defensin genes on human chromosome 8p23.

Defensins are a large family of host defense peptides expressed in leukocytes and epithelia. Using P1 and BAC clones, we have determined the organization of the human alpha-defensin genes and the beta-defensin gene HDEFB1 on chromosome 8p23. From the telomere, the order of the genes (with encoded peptides in parentheses) is HDEFA5 (HD-5), HDEFA1/1A (HNP-1/3), HDEFA4 (HNP-4), HDEFA6 (HD-6), and HDEFB1 (HBD-1). These genes span a region of approximately 450kb. Genes encoding intestinal Paneth cell defensins (HDEFA5 and HDEFA6) flank the myeloid defensin gene cluster (HDEFA1, HDEFA1A, HDEFA4). Based on our previous studies, the remaining known defensin gene, HDEFB2 (HBD-2), is about 400kb centromeric to HDEFB1. This map supports the hypothesis, originally proposed because of sequence similarities, that myeloid alpha-defensin genes evolved by reduplication and divergence from Paneth cell defensin genes, and identifies regions and clones, which should be useful in the search for new defensin genes.     

IL-4 and IL-13 negatively regulate TNF-alpha- and IFN-gamma-induced beta-defensin expression through STAT-6, suppressor of cytokine signaling (SOCS)-1, and SOCS-3

Human beta-defensins (HBDs) are a major class of antimicrobial peptides that play an important role in the innate immune response, however, the induction and regulation of these antimicrobial peptides is not well understood. We demonstrate here that stimulation of keratinocytes with TNF-alpha/IFN-gamma induces HBD-2 and HBD-3 by activating STAT-1 and NF-kappaB signaling. We further demonstrate that IL-4 and IL-13 activate STAT-6 and induce the suppressors of cytokine signaling (SOCS)-1 and -3. This interferes with STAT-1 and NF-kappaB signaling, thereby inhibiting TNF-alpha/IFN-gamma-mediated induction of HBD-2 and HBD-3. These data suggest that targeting the STAT-1-signaling pathway or suppressor of cytokine signaling expression enhances beta-defensin expression and represents a new therapeutic strategy for reduction of infection in human diseases associated with beta-defensin deficiency.     

A novel murine beta -defensin expressed in tongue, esophagus, and trachea

beta-Defensins are broad spectrum antimicrobial peptides expressed at epithelial surfaces. Two human beta-defensins, HBD-1 and HBD-2, have been identified. In the lung, HBD-2 is an inducible product of airway epithelia and may play a role in innate mucosal defenses. We recently characterized rat homologs (RBD-1, RBD-2) of the human genes and used these sequences to identify novel mouse genes. Mouse beta-defensin-4 (MBD-4) was amplified from lung cDNA using polymerase chain reaction primers designed from conserved sequences of RBD-2 and HBD-2. A full-length cDNA was cloned which encodes a putative peptide with the sequence MRIHYLLFTFLLVLLSPLAAFTQIINNPITCMTNGAICWGPCPTAFRQIGNCGHFKVRCCKIR. The peptide shares approximately 40% identity with HBD-2. MBD-4 mRNA was expressed in the esophagus, tongue, and trachea but not in any of 20 other tissues surveyed. Cloning of the genomic sequence of MBD-4 revealed two nearly (>99%) identical sequences encoding MBD-4 and the presence of numerous additional highly similar genomic sequences. Radiation hybrid mapping localized this gene to a region of chromosome 8 near several other defensins, MBD-2, MBD-3, and alpha-defensins (cryptdins)-3 and -17, consistent with a gene cluster. Our genomic cloning and mapping data suggest that there is a large beta-defensin gene family in mice. Identification of murine beta-defensins provides an opportunity to understand further the role of these peptides in host defense through animal model studies and the generation of beta-defensin-deficient animals by gene targeting.     

Copy Number Variation of the Beta-Defensin Genes in Europeans: No Supporting Evidence for Association with Lung Function,Chronic Obstructive Pulmonary Disease or Asthma

Lung function measures are heritable, predict mortality and are relevant in diagnosis of chronic obstructive pulmonary disease (COPD). COPD and asthma are diseases of the airways with major public health impacts and each have a heritable component. Genome-wide association studies of SNPs have revealed novel genetic associations with both diseases but only account for a small proportion of the heritability. Complex copy number variation may account for some of the missing heritability. A well-characterised genomic region of complex copy number variation contains beta-defensin genes (DEFB103, DEFB104 and DEFB4), which have a role in the innate immune response. Previous studies have implicated these and related genes as being associated with asthma or COPD. We hypothesised that copy number variation of these genes may play a role in lung function in the general population and in COPD and asthma risk. We undertook copy number typing of this locus in 1149 adult and 689 children using a paralogue ratio test and investigated association with COPD, asthma and lung function. Replication of findings was assessed in a larger independent sample of COPD cases and smoking controls. We found evidence for an association of beta-defensin copy number with COPD in the adult cohort (OR = 1.4, 95%CI:1.02–1.92, P = 0.039) but this finding, and findings from a previous study, were not replicated in a larger follow-up sample(OR = 0.89, 95%CI:0.72–1.07, P = 0.217). No robust evidence of association with asthma in children was observed. We found no evidence for association between beta-defensin copy number and lung function in the general populations. Our findings suggest that previous reports of association of beta-defensin copy number with COPD should be viewed with caution. Suboptimal measurement of copy number can lead to spurious associations. Further beta-defensin copy number measurement in larger sample sizes of COPD cases and children with asthma are needed.     

Involvement of Endoplasmic Reticulum Stress in Inflammatory Bowel Disease: A Different Implication for Colonic and Ileal Disease?

Background

Endoplasmic reticulum (ER) stress has been suggested to play a role in inflammatory bowel disease (IBD). The three branches (ATF6, IRE1 and PERK) of the unfolded protein response (UPR) have different roles and are not necessarily activated simultaneously.

Methodology/Principal Findings

Expression of UPR-related genes was investigated in colonic and ileal biopsies of 23 controls, 15 ulcerative colitis (UC) and 54 Crohn''s disease (CD) patients. This expression was confirmed at protein level in colonic and ileal samples of five controls, UC and CD patients. HSPA5, PDIA4 and XBP1s were significantly increased in colonic IBD at mRNA and/or protein levels, indicating activation of the ATF6 and IRE1 branch. Colonic IBD was associated with increased phosphorylation of EIF2A suggesting the activation of the PERK branch, but subsequent induction of GADD34 was not observed. In ileal CD, no differential expression of the UPR-related genes was observed, but our data suggested a higher basal activation of the UPR in the ileal mucosa of controls. This was confirmed by the increased expression of 16 UPR-related genes as 12 of them were significantly more expressed in ileal controls compared to colonic controls. Tunicamycin stimulation of colonic and ileal samples of healthy individuals revealed that although the ileal mucosa is exhibiting this higher basal UPR activation, it is still responsive to ER stress, even more than colonic mucosa.

Conclusions/Significance

Activation of the three UPR-related arms is seen in colonic IBD-associated inflammation. However, despite EIF2A activation, inflamed colonic tissue did not increase GADD34 expression, which is usually involved in re-establishment of ER homeostasis. This study also implies the presence of a constitutive UPR activation in healthy ileal mucosa, with no further activation during inflammation. Therefore, engagement of the UPR differs between colon and ileum and this could be a factor in the development of ileal or colonic disease.     

A flow cytometric assay to monitor antimicrobial activity of defensins and cationic tissue extracts

To determine the antibacterial activity of defensins and other antimicrobial peptides in biopsy extracts, we evaluated a flow cytometric method with the membrane potential sensitive dye bis-(1,3-dibutylbarbituric acid) trimethine oxonol [DiBAC4(3)]. This assay enables us to discriminate intact non-fluorescent and depolarized fluorescent bacteria after exposure to antimicrobial peptides by measurement at the direct target, the cytoplasmic membrane and the membrane potential. The feasibility of the flow cytometric assay was evaluated with recombinant human beta-defensin 3 (HBD-3) against 25 bacterial strains representing 12 species. HBD-3 showed a broad-spectrum dose dependent activity and the minimal dose to cause depolarization ranged from 1.25 to >15 microg/ml HBD-3, depending on the species tested. The antibacterial effect was diminished with sodium chloride or dithiothreitol and could be abrogated with a HBD-3 antibody. Additionally, isolated cationic extracts from human intestinal biopsies showed a strong bactericidal effect against Escherichia coli K12, E. coli ATCC 25922 and Staphylococcus aureus ATCC 25923, which was diminished towards E. coli at 150 mM NaCl, whereas the activity towards S. aureus ATCC 25923 remained unaffected at physiological salt concentrations. DTT blocked the bactericidal effect of biopsy extracts completely.     

IBD5 is a general risk factor for inflammatory bowel disease: replication of association with Crohn disease and identification of a novel association with ulcerative colitis

Inflammatory bowel disease (IBD) refers to complex chronic relapsing autoimmune disorders of the gastrointestinal tract that have been traditionally classified into Crohn disease (CD) and ulcerative colitis (UC). We have previously reported that genetic variation within a 250-kb haplotype (IBD5) in the 5q31 cytokine gene cluster confers susceptibility to CD in a Canadian population. In the current study, we first replicated this association by examining 368 German trios with CD and demonstrating, by transmission/disequilibrium testing (TDT), that the same haplotype is associated with CD (chi2=5.97; P=.007). Our original association study focused on the role of IBD5 in CD; we next explored the potential contribution of this locus to UC susceptibility in 187 German trios. Given the TDT results in the present cohort with UC, IBD5 may also act as a susceptibility locus for UC (chi2=8.10; P=.002). We then examined locus-locus interactions between IBD5 and CARD15, a locus reported elsewhere to confer risk exclusively to CD. Our current results indicate that the two loci act independently to confer risk to CD but that these two loci may behave in an epistatic fashion to promote the development of UC. Moreover, IBD5 was not associated with particular clinical manifestations upon phenotypic stratification in the current cohort with CD. Taken together, our results suggest that IBD5 may act as a general risk factor for IBD, with loci such as CARD15 modifying the clinical characteristics of disease.     

Dysregulation of Human β-Defensin-2 Protein in Inflammatory Bowel Disease

Background

Human β-defensin-2 (HBD2) is an antimicrobial peptide implicated in the pathogenesis of inflammatory bowel disease (IBD). Low copy number and concomitant low mRNA expression of the HBD2 gene have been implicated in susceptibility to colonic Crohn''s Disease (CD). We investigated the colonic distribution of HBD2 mRNA expression, and the contributions of genetic and environmental factors on HBD2 protein production.

Methodology/Principal Findings

We examined HBD2 mRNA expression at three colonic locations by microarray analysis of biopsies from 151 patients (53 CD, 67 ulcerative colitis [UC], 31 controls). We investigated environmental and genetic influences on HBD2 protein production using ex vivo cultured sigmoid colon biopsies from 69 patients (22 CD, 26 UC, 21 controls) stimulated with lipopolysaccharide (LPS) and/or nicotine for 24 hours. HBD2 and cytokines were measured in culture supernatants. Using DNA samples from these patients, regions in the HBD2 gene promoter were sequenced for NF-κB binding-sites and HBD2 gene copy number was determined. HBD2 mRNA expression was highest in inflamed (vs. uninflamed p = 0.0122) ascending colon in CD and in inflamed (vs. uninflamed p<0.0001) sigmoid colon in UC. HBD2 protein production was increased in inflamed UC biopsies (p = 0.0078). There was no difference in HBD2 protein production from unstimulated biopsies of CD, UC and controls. LPS-induced HBD2 production was significantly increased in CD (p = 0.0375) but not UC (p = 0.2017); this LPS-induced response was augmented by nicotine in UC (p = 0.0308) but not CD (p = 0.6872). Nicotine alone did not affect HBD2 production. HBD2 production correlated with IL8 production in UC (p<0.001) and with IL10 in CD (p<0.05). Variations in the HBD2 promoter and HBD2 gene copy number did not affect HBD2 production.

Significance/Conclusions

Colonic HBD2 was dysregulated at mRNA and protein level in IBD. Inflammatory status and stimulus but not germline variations influenced these changes.     

Induction of human beta-defensin-2 expression in human astrocytes by lipopolysaccharide and cytokines

Defensins are cationic peptides with broad-spectrum antimicrobial activity. They are members of a supergene family consisting of alpha and beta subtypes and each subtype is comprised of a number of different isoforms. For example, human alpha-defensin (HAD) has six isoforms, which are expressed by polymorphonuclear leukocytes and Paneth cells. In contrast, human beta-defensin (HBD) has two isoforms that are expressed by epithelial cells of the skin, gut, respiratory and urogenital tracts. Recently, HBD-1 was detected in human brain biopsy tissue. However, little is known about the expression of HBD-1 or HBD-2 in the CNS and whether neural cells can secrete these peptides. For the present study, human astrocyte, microglial, meningeal fibroblast and neuronal cultures were probed for the expression of HBD-1 and HBD-2 mRNA and protein. Each cell type was either maintained in tissue culture medium alone or in medium containing lipopolysaccharide (LPS) at concentrations ranging from 0.1 to 1 microgram/mL, interleukin-1 beta (IL-1beta) at 1-50 ng/mL, or tumor necrosis factor alpha (TNF-alpha) at the same concentrations. The expression of HBD-1 and HBD-2 mRNAs was monitored by RT-PCR. The cDNA products were sequenced to characterize the gene product. HBD-2 protein was detected by immunoblot, immunoprecipitation and immunocytochemistry. Results of these studies showed that HBD-1 mRNA was detected in all cell cultures except in those enriched for neurons. In contrast, HBD-2 mRNA was detected only in astrocyte cultures that were treated with LPS, IL-1beta or TNF-alpha. The detection of the respective proteins correlated positively with the mRNA results. As such, these data represent the first demonstration of HBD-2 expression by astrocytes and suggest that this peptide may play a role in host defense against bacterial CNS pathogenesis.     

Colonic expression and synthesis of interleukin 13 and interleukin 15 in inflammatory bowel disease

A dysregulated local immune reaction with unbalanced cytokine expression seems essential in inflammatory bowel disease (IBD), i.e. ulcerative colitis (UC) and Crohn's disease (CD). Since the roles of interleukin (IL-)13 and IL-15 remain unclear, this study aimed at studying intestinal expression of IL-13 and IL-15 in IBD. Methods: In colonic biopsies from 24 UC, 18 CD, and 12 controls IL-13 and IL-15 were measured using ELISA, and their gene expressions were assessed by RT-PCR. Leukocytes were visualised histochemically. Results: Concentrations of IL-13 were decreased in UC (median 56 pg/mg tissue; interquartile range 30-99 pg/mg) compared to CD (82 pg/mg tissue; 41-122;P=0.004) and controls (83 pg/mg tissue; 18-134;P>0.05), and lower in active UC (53 pg/mg tissue; 33-96) than in inactive UC (80 pg/mg tissue; 65-99;P=0.02). IL-15 concentrations were higher in CD patients (34 pg/mg tissue; 24-53) as compared to controls (20 pg/mg tissue; 15-21;P=0.001) whilst being 22 pg/mg tissue (15-32) in UC. IL-13 mRNA and IL-15 mRNA were detected in 20% and 15%, respectively. Infiltration of leukocytes correlated inversely with IL-13 levels (P=0.02). Conclusion: Active UC is associated with decreased colonic IL-13 suggesting that IL-13 levels are diminished as a part of UC exacerbations, or that exacerbations follow active downregulation of IL-13.     

Genetic evidence for interaction of the 5q31 cytokine locus and the CARD15 gene in Crohn disease

A common haplotype spanning 250 kb in the cytokine gene cluster on chromosome 5q31 has recently been reported to be strongly associated with Crohn disease (CD) in Canadian families. We have replicated this finding by both the transmission-disequilibrium test (TDT) (P=.016) and in a case-control association study (P=.008) in a large European cohort of patients with CD, although the increase in disease risk was small (odds ratio 1.49 for homozygotes, 95% CI 1.11-2.0). No association was detected in families or individuals with ulcerative colitis (UC). Stratification of offspring with CD in the TDT sample by mutation status in the CD susceptibility gene CARD15 showed that the association with the 5q31 risk haplotype was present only in offspring with at least one of the known CARD15 disease susceptibility alleles (P=.044). The 5q31 risk haplotype frequency was 53.1% in unrelated individuals with CD who had one or two CARD15 mutations versus 43.7% in control subjects (P=.0001) but was not significantly elevated in individuals with CD who had no CARD15 mutations (45.4%, P=.41). Kaplan-Meier survival analysis of age at disease onset showed a significantly earlier onset in homozygotes for the 5q31 risk haplotype (P=.0019). These findings suggest that genetic variants at the 5q31 (IBD5) locus may hasten the onset of Crohn disease and cooperate with CARD15 in disease causation.     

An ileal Crohn's disease gene signature based on whole human genome expression profiles of disease unaffected ileal mucosal biopsies

Previous genome-wide expression studies have highlighted distinct gene expression patterns in inflammatory bowel disease (IBD) compared to control samples, but the interpretation of these studies has been limited by sample heterogeneity with respect to disease phenotype, disease activity, and anatomic sites. To further improve molecular classification of inflammatory bowel disease phenotypes we focused on a single anatomic site, the disease unaffected proximal ileal margin of resected ileum, and three phenotypes that were unlikely to overlap: ileal Crohn's disease (ileal CD), ulcerative colitis (UC), and control patients without IBD. Whole human genome (Agilent) expression profiling was conducted on two independent sets of disease-unaffected ileal samples collected from the proximal margin of resected ileum. Set 1 (47 ileal CD, 27 UC, and 25 Control non-IBD patients) was used as the training set and Set 2 was subsequently collected as an independent test set (10 ileal CD, 10 UC, and 10 control non-IBD patients). We compared the 17 gene signatures selected by four different feature-selection methods to distinguish ileal CD phenotype with non-CD phenotype. The four methods yielded different but overlapping solutions that were highly discriminating. All four of these methods selected FOLH1 as a common feature. This gene is an established biomarker for prostate cancer, but has not previously been associated with Crohn's disease. Immunohistochemical staining confirmed increased expression of FOLH1 in the ileal epithelium. These results provide evidence for convergent molecular abnormalities in the macroscopically disease unaffected proximal margin of resected ileum from ileal CD subjects.     

Cytomegalovirus viremia as a risk factor for mortality prior to antiretroviral therapy among HIV-infected gold miners in South Africa

Background

Cytomegalovirus (CMV) viremia has been shown to be an independent risk factor for increased mortality among HIV-infected individuals in the developing world. While CMV infection is nearly ubiquitous in resource-poor settings, few data are available on the role of subclinical CMV reactivation on HIV.

Methods

Using a cohort of mineworkers with stored plasma samples, we investigated the association between CMV DNA concentration and mortality prior to antiretroviral therapy availability.

Results

Among 1341 individuals (median CD4 count 345 cells/µl, 70% WHO stage 1 or 2, median follow-up 0.9 years), 70 (5.2%) had CMV viremia at baseline; 71 deaths occurred. In univariable analysis CMV viremia at baseline was associated with a three-fold increase in mortality (hazard ratio [HR] 3.37; 95% confidence intervals [CI] 1.60, 7.10). After adjustment for CD4 count, WHO stage and HIV viral load (N = 429 with complete data), the association was attenuated (HR 2.27; 95%CI 0.88, 5.83). Mortality increased with higher CMV viremia (≥1,000 copies/ml vs. no viremia, adjusted HR 3.65, 95%CI: 1.29, 10.41). Results were similar using time-updated CMV viremia.

Conclusions

High copy number, subclinical CMV viremia was an independent risk factor for mortality among male HIV-infected adults in South Africa with relatively early HIV disease. Studies to determine whether anti-CMV therapy to mitigate high copy number viremia would increase lifespan are warranted.