Willemijn B. Breunis

A Genome-Wide Association Study Identifies Novel and Functionally Related Susceptibility Loci for Kawasaki Disease

Kawasaki disease (KD) is a pediatric vasculitis that damages the coronary arteries in 25% of untreated and approximately 5% of treated children. Epidemiologic data suggest that KD is triggered by unidentified infection(s) in genetically susceptible children. To investigate genetic determinants of KD susceptibility, we performed a genome-wide association study (GWAS) in 119 Caucasian KD cases and 135 matched controls with stringent correction for possible admixture, followed by replication in an independent cohort and subsequent fine-mapping, for a total of 893 KD cases plus population and family controls. Significant associations of 40 SNPs and six haplotypes, identifying 31 genes, were replicated in an independent cohort of 583 predominantly Caucasian KD families, with NAALADL2 (rs17531088, p combined = 1.13×10−6) and ZFHX3 (rs7199343, p combined = 2.37×10−6) most significantly associated. Sixteen associated variants with a minor allele frequency of >0.05 that lay within or close to known genes were fine-mapped with HapMap tagging SNPs in 781 KD cases, including 590 from the discovery and replication stages. Original or tagging SNPs in eight of these genes replicated the original findings, with seven genes having further significant markers in adjacent regions. In four genes (ZFHX3, NAALADL2, PPP1R14C, and TCP1), the neighboring markers were more significantly associated than the originally associated variants. Investigation of functional relationships between the eight fine-mapped genes using Ingenuity Pathway Analysis identified a single functional network (p = 10−13) containing five fine-mapped genes—LNX1, CAMK2D, ZFHX3, CSMD1, and TCP1—with functional relationships potentially related to inflammation, apoptosis, and cardiovascular pathology. Pair-wise blood transcript levels were measured during acute and convalescent KD for all fine-mapped genes, revealing a consistent trend of significantly reduced transcript levels prior to treatment. This is one of the first GWAS in an infectious disease. We have identified novel, plausible, and functionally related variants associated with KD susceptibility that may also be relevant to other cardiovascular diseases.

Genome-wide association study identifies variants in the CFH region associated with host susceptibility to meningococcal disease

Meningococcal disease is an infection caused by Neisseria meningitidis. Genetic factors contribute to host susceptibility and progression to disease, but the genes responsible for disease development are largely unknown. We report here a genome-wide association study for host susceptibility to meningococcal disease using 475 individuals with meningococcal disease (cases) and 4,703 population controls from the UK. We performed, in Western European and South European cohorts (consisting of 968 cases and 1,376 controls), two replication studies for the most significant SNPs. A cluster of complement factor SNPs replicated independently in both cohorts, including SNPs within complement factor H (CFH) (rs1065489 (p.936D<E), P = 2.2 × 10−11) and in CFH-related protein 3 (CFHR3)(rs426736, P = 4.6 × 10−13). N. meningitidis is known to evade complement-mediated killing by the binding of host CFH to the meningococcal factor H–binding protein (fHbp). Our study suggests that host genetic variation in these regulators of complement activation plays a role in determining the occurrence of invasive disease versus asymptomatic colonization by this pathogen.

Copy number variation at the FCGR locus includes FCGR3A, FCGR2C and FCGR3B but not FCGR2A and FCGR2B.

Human Fcã receptors (FcãRs) are glycoproteins that bind the Fc region of IgG. The genes encoding the low‐affinity FcãRs are located on chromosome 1q23‐24. Beside single nucleotide polymorphisms (SNPs), gene copy number variation (CNV) is now being recognized as an important indicator for inter‐individual differences. Recent studies on identifying CNV in the human genome suggest large areas at chromosome 1q23‐24 to be involved, and CNV in this region has been associated with manifestations of systemic autoimmune disease. To study both SNPs and CNV of the low‐affinity FcãRs in one assay, we have developed a Multiplex Ligation‐dependent Probe Amplification (MLPA) assay. A novel CNV for FCGR3A was observed. Similar to FCGR3B and FCGR2C, a gene‐dosage effect of FCGR3A was found, that seemed to correlate nicely with the FcãRIIIa expression on NK cells. Next, we delineated the approximate boundaries of CNV at the FCGR locus. Variation in co‐segregation of neighboring FCGR genes was limited to four variants, with patterns of Mendelian inheritance. No CNV of the FCGR2A and FCGR2B genes was observed in over 600 individuals. In conclusion, we report a novel CNV of the FCGR3A gene that correlates with FcãRIIIa expression and function on NK cells. Only FCGR3A, FCGR2C and FCGR3B show CNV, in contrast to FCGR2A and FCGR2B.

Transforming Growth Factor-β Signaling Pathway in Patients With Kawasaki Disease

Background—Transforming growth factor (TGF)-&bgr; is a multifunctional peptide that is important in T-cell activation and cardiovascular remodeling, both of which are important features of Kawasaki disease (KD). We postulated that variation in TGF-&bgr; signaling might be important in KD susceptibility and disease outcome. Methods and Results—We investigated genetic variation in 15 genes belonging to the TGF-&bgr; pathway in a total of 771 KD subjects of mainly European descent from the United States, the United Kingdom, Australia, and the Netherlands. We analyzed transcript abundance patterns using microarray and reverse transcriptase–polymerase chain reaction for these same genes, and measured TGF-&bgr;2 protein levels in plasma. Genetic variants in TGFB2, TGFBR2, and SMAD3 and their haplotypes were consistently and reproducibly associated with KD susceptibility, coronary artery aneurysm formation, aortic root dilatation, and intravenous immunoglobulin treatment response in different cohorts. A SMAD3 haplotype associated with KD susceptibility replicated in 2 independent cohorts and an intronic single nucleotide polymorphism in a separate haplotype block was also strongly associated (A/G, rs4776338) (P=0.000022; odds ratio, 1.50; 95% confidence interval, 1.25 to 1.81). Pathway analysis using all 15 genes further confirmed the importance of the TGF-&bgr; pathway in KD pathogenesis. Whole-blood transcript abundance for these genes and TGF-&bgr;2 plasma protein levels changed dynamically over the course of the illness. Conclusions—These studies suggest that genetic variation in the TGF-&bgr; pathway influences KD susceptibility, disease outcome, and response to therapy, and that aortic root and coronary artery Z scores can be used for phenotype/genotype analyses. Analysis of transcript abundance and protein levels further support the importance of this pathway in KD pathogenesis.

Influence of Cytotoxic T Lymphocyte-associated Antigen 4 (CTLA4) Common Polymorphisms on Outcome in Treatment of Melanoma Patients With CTLA-4 Blockade

Blockade of the cytotoxic T lymphocyte-associated antigen 4 (CTLA-4), a down-regulator of T-cell activation, can cause cancer regression in patients with metastatic melanoma. However, not all patients respond well to the therapy and some develop severe autoimmune reactions. We hypothesized that common genetic variation in the CTLA4 gene could contribute to response to CTLA-4 blockade and the occurrence of autoimmune reactions. We investigated 7 common single nucleotide polymorphisms, SNPs, (rs733618, rs4553808, rs11571317, rs5742909, rs231775, rs3087243, and rs7565213) in 152 white melanoma patients who received CTLA-4 blockade. Three SNPs were associated with response to therapy: proximal promoter SNPs, rs4553808 [P=0.002; odds ratio (OR) 3.39; 95% confidence interval (CI), 1.62-7.10] and rs11571327 (P=0.02; OR 2.89; 95% CI, 1.23-6.83) and the nonsynonymous SNP rs231775 (Thr17Ala, P=0.009; OR 0.39; 95% CI, 0.18-0.82). A haplotype analysis including the 7 SNPs suggested that the common haplotype, TACCGGG could be associated with no response (P=0.02) whereas the haplotype TGCCAGG (P=0.06; OR 4.13; 95% CI, 1.17-14.5) could be associated with response to the treatment. No significant association was observed for occurrence of severe autoimmune reactions (grade III/IV) either by single SNP or haplotype analyses. Our results suggest that genetic variation in CTLA4 could influence response to CTLA-4 blockade therapy in metastatic melanoma patients, but further studies are necessary to confirm the observed associations.

Phenotypic Variation in IgG Receptors by Nonclassical FCGR2C Alleles

The balance between activating and inhibitory signals from the different FcγRs for IgG ensures homeostasis of many inflammatory responses. FCGR2C is the product of an unequal crossover of the FCGR2A and FCGR2B genes encoding the activating FcγRIIa (CD32a) and inhibitory FcγRIIb (CD32b), respectively. A single nucleotide polymorphism (SNP) in exon 3 of FCGR2C results in either expression of the activating FcγRIIc (CD32c) (FCGR2C-open reading frame [ORF]) or its absence because of a stop codon (FCGR2C-Stop). Two additional variations in FcγRIIb/c expression on leukocytes have now been identified. In case of “nonclassical” FCGR2C-ORF alleles, FcγRIIc expression was unexpectedly absent, because of novel splice site mutations near exon 7 leading to another stop codon. In some individuals with FCGR2C-Stop alleles FcγRIIb was detected on NK cells, which normally are devoid of this protein. Individuals with these nonclassical FCGR2C-Stop alleles carried a deletion of FCGR2C-FCGR3B that extends into the promoter region of the adjacent FCGR2B gene and probably deletes a negative regulatory element in the FCGR2B promoter in NK cells. FcγRIIb expression on NK cells effectively inhibited killing mediated by FcγRIIIa (CD16a) in Ab-dependent cytotoxicity tests. Our findings demonstrate a more extensive and previously unnoticed variation in FcγR expression with relevance to immunity and inflammation.

Polymorphisms in chemokine receptor genes and susceptibility to Kawasaki disease.

Kawasaki disease (KD) is an acute vasculitis occurring in young children. Its aetiology is unknown, but an infectious agent is assumed. Increased levels of proinflammatory cytokines and chemokines have been reported in KD. Genetic variation in these genes and the receptors for these genes could influence the regulation of cytokines and chemokines. In a case–control study of 170 Dutch Caucasian KD patients and 300 healthy Dutch Caucasian controls, common genetic variants in chemokine receptor genes CCR3, CCR2, CCR5, CX3CR1, CXCR1 and CXCR2 were analysed. Of the eight studied single nucleotide polymorphisms (SNPs) in the CCR3–CCR2–CCR5 gene cluster, four showed a significant association with susceptibility to KD. Moreover the CCR5‐Δ32 was observed with an allele frequency of 10·7% in the control population compared to 6·5% in the KD patients (P = 0·04). Two haplotypes of the CCR3–CCR2–CCR5 gene‐cluster appear to be at risk haplotypes for KD and one a protective haplotype. No association was observed with the studied SNPs in CX3CR1, CXCR1 and CXCR2. In conclusion, in a Dutch cohort of KD patients an association of KD occurrence with common genetic variants in the chemokine receptor gene‐cluster CCR3–CCR2–CCR5 was observed.

An analysis of genetic variation across the MBL2 locus in Dutch Caucasians indicates that 3¢ haplotypes could modify circulating levels of mannose-binding lectin

Mannose-binding protein (MBL) is a critical component of innate immunity and provides first-line protection against pathogens. Both circulating MBL serum levels and functional activity have been correlated with common genetic variants in the MBL2 gene. Associations between MBL deficiency and severe infections have been reported in immuno-incompetent patients and for autoimmune disorders; however, measured MBL serum levels do not fully correlate with the ‘secretor haplotypes’. Previously, the MBL2 locus was resequenced and determined that a recombination hotspot divides MBL2 into two haplotype blocks. It was sought to investigate whether additional variants, in either block structure could associate with MBL serum levels. Therefore, 31 common variants were analysed across the locus in 212 DNA samples of healthy Caucasian individuals with known MBL serum concentrations. The additional 5′ variants were in strong linkage to the elements of the ‘secretor haplotypes’; functional alleles B, C and D also lie on restricted haplotypes. Four variants in the 3′ block (Ex4-1483T>C, Ex4-1067G>A, Ex4-901G>A and Ex4-710G>A) are components of a distinct haplotype block. The results of this study suggest that additional 5′ variants as well as markers of distinct 3′ haplotype blocks in MBL2 may contribute to circulating protein levels, but further studies are required to confirm these observations. Last, there could be a selective advantage for diversification of the 3′ region of the gene.

Genome-wide linkage and association mapping identify susceptibility alleles in ABCC4 for Kawasaki disease

Background Kawasaki disease (KD) is a self limited vasculitis in which host genetics plays a prominent role. To further the understanding of the role of host genetics in KD, a three-stage genetic study was conducted that began with a family linkage study and ultimately involved more than 3000 individuals to identify new genetic contributions to KD susceptibility. Methods and results A 26-family linkage study followed by fine mapping was performed in a cohort of 1284 KD subjects and their family members (total 3248 individuals). Suggestive evidence of disease linkage (logarithm of odds (LOD) ≥3.0, p<1.00×10−4) was found for five genomic locations (Chr 3q, 4q, 10p, 13q, 21q). Two of these loci (Chr 4q and Chr 13q) overlapped with validated findings from a recent KD genome-wide association study. Fine mapping analysis revealed three single nucleotide polymorphisms (SNPs) in ATP-binding cassette, subfamily C, member 4 (ABCC4) underlying the Chr 13q linkage peak showing evidence of association to KD (lowest p=8.82×10−5; combined OR 2.00, 95% CI 1.41 to 2.83). ABCC4 is a multifunctional cyclic nucleotide transporter that stimulates the migratory capacity of dendritic cells. It is also a mediator of prostaglandin efflux from human cells and is inhibited by non-steroidal anti-inflammatory medications such as aspirin. Conclusion These genetic data suggest that ABCC4 could play a fundamental role in KD pathogenesis with effects on immune activation and vascular response to injury.

Five years of Kawasaki disease in the Netherlands a national surveillance study

Background: The aim of this study was to evaluate the incidence, disease presentation, treatment and cardiac outcome of Kawasaki disease (KD) in The Netherlands. Methods: The national Dutch Pediatric Surveillance Unit was used to prospectively register new KD cases from 2008 through 2012. Questionnaires were sent to pediatricians to obtain clinical information. Results: Nationwide 341 cases were reported during the 5-year study period, of which 319 questionnaires (93.0%) were returned. The mean incidence of KD was estimated to be 5.8/100,000 children <5 years of age. The median age at disease onset was 2.4 years (range 0.1–14.6 years) and 79.2% of cases were <5 years of age. The male-to-female ratio was 1.5 to 1. Incomplete KD was diagnosed in 22.3% of cases and these cases were significantly younger than complete cases [median: 1.1 (0.1–13.7) vs. 2.8 (0.2–14.6) years, P < 0.001]. In total, 308 patients (96.6%) received intravenous immunoglobulins (IVIG). Retreatment with IVIG was given in 71 (23.1%) and additional steroid treatment in 17 patients (5.5%). During the acute phase, coronary artery aneurysms developed in 43 cases (13.5%). Multivariate logistic regression analysis showed that male gender, delay of treatment (>10 days) and IVIG retreatment were independent risk factors for coronary artery aneurysms development. Conclusions: This prospective study of KD in The Netherlands revealed a mean annual incidence of 5.8/100,000 children <5 years of age. Clinicians should consider the diagnosis of KD in young (male) children with persistent inexplicable fever to start IVIG treatment within 10 days to prevent development of coronary artery aneurysms.