key: cord-0684698-9i8zpdfy
authors: Hoh, B.P.; Umi-Shakina, H.; Zuraihan, Z.; Zaiharina, M.Z.; Rafidah-Hanim, S.; Mahiran, M.; Nik Khairudin, N.Y.; Benedict, L.H. Sim; Masliza, Z.; Christopher, K.C. Lee; Sazaly, A.B.
title: Common variants of chemokine receptor gene CXCR3 and its ligands CXCL10 and CXCL11 associated with vascular permeability of dengue infection in peninsular Malaysia
date: 2015-04-07
journal: Hum Immunol
DOI: 10.1016/j.humimm.2015.03.019
sha: f61afe136befbbe424ddd37c5b92411741d705d8
doc_id: 684698
cord_uid: 9i8zpdfy

Dengue causes significantly more human disease than any other arboviruses. It causes a spectrum of illness, ranging from mild self-limited fever, to severe and fatal dengue hemorrhagic fever, as evidenced by vascular leakage and multifactorial hemostatic abnormalities. There is no specific treatment available till date. Evidence shows that chemokines CXCL10, CXCL11 and their receptor CXCR3 are involved in severity of dengue, but their genetic association with the susceptibility of vascular leakage during dengue infection has not been reported. We genotyped 14 common variants of these candidate genes in 176 patients infected with dengue. rs4859584 and rs8878 (CXCL10) were significantly associated with vascular permeability of dengue infection (P < 0.05); while variants of CXCL11 showed moderate significance of association (P = 0.0527). Haplotype blocks were constructed for genes CXCL10 and CXCL11 (5 and 7 common variants respectively). Haplotype association tests performed revealed that, “CCCCA” of gene CXCL10 and “AGTTTAC” of CXCL11 were found to be significantly associated with vascular leakage (P = 0.0154 and 0.0366 respectively). In summary, our association study further strengthens the evidence of the involvement of CXCL10 and CXCL11 in the pathogenesis of dengue infection.

Dengue virus (DENV) is a mosquito-borne virus in the genus of Flavivirus, family Flaviviridae, and consists of four related serotypes: DENV1 -DENV4. This disease causes a spectrum of illness ranging from asymptomatic infection or mild febrile illness to severe and fatal hemorrhagic disease. While majority experience uncomplicated Dengue Fever (DF), patients could potentially progress to severe clinical manifestations, which include occurrence of vascular permeability defect resulting in plasma leakage, and multifactorial hemostatic abnormalities [1] . There is no specific treatment available till date.

Dengue causes significantly more human disease than any other aboviruses, and it is endemic in more than 100 countries, including most of the Southeast Asia, South and Central America, the Caribbean and South Pacific regions; while dengue hemorrhagic fever/dengue shock syndrome (DHF/DSS) has been reported in more than 60 countries [2] . Annually, approximately 50-100 million cases of severe dengue patients require hospitalization in which 500,000 resulted in DHF/DSS, with more than 20,000 deaths worldwide [2] . In Malaysia, 46,171 cases with 134 deaths were reported in 2010 [3] .

The severe form of dengue infection has been associated with various factors including a robust host inflammatory immune response. ''Original antigenic sin'' and antibody-dependent enhancement (ADE) of viral replication is the most widely accepted explanation for the association between DHF and pre-existing http antibody [4] [5] [6] [7] . However, it remains uncertain as to how the virushost interaction triggers the inflammatory response resulting in vascular leakage, the hallmark of DHF/DSS, as evidenced by haemoconcentration and/or effusions in the pleural or peritoneal spaces. The fact that the Cuban, Caribbean, and African black populations which share a common ancestry are protected against DHF explains at least partly the role of host genetic variability in this disease susceptibility [8] .

A study reported by Fink et al. [9] suggested the involvement of some host candidate genes associated with dengue infection namely, CXCR3 chemokine receptor gene and its ligands chemokine CXCL10 (IP10) and CXCL11 (I-TAC). These genes are commonly involved in the NFkB pathway. Production of these chemokines lead to recruitment of CXCR3 expressing T-cells and NK cells. On top of the in vivo study, earlier studies on gene expression also suggest the involvement of these chemokines in SARS, West Nile virus encephalitis and chronic Hepatitis C patients [10] [11] [12] . We hypothesized that the differential clinical expression of the selected chemokines in different individuals are related to the genetic variation of the chemokine receptor CXCR3 gene, and its ligands CXCL10 and CXCL11. We therefore, attempted to investigate if the risk of acquiring vascular leakage -the hallmark symptom of progression to severe dengue -is attributed to genetic variations in these candidate genes.

One hundred and seventy-six hospitalized suspected adult dengue-patients, aged 13 years old and above were recruited from three different hospitals namely, Hospital Universiti Sains Malaysia (HUSM), Hospital Kota Bharu (HKB) and Hospital Sungai Buloh (HSB). Patients were from HUSM and HKB were recruited between 2008 and 2010. Both hospitals are located in the state of Kelantan; while patients from HSB located in the state of Selangor, were recruited between 2010 and 2012. Blood samples were collected on the first day of admission while the clinical data were retrieved on standardized case report forms, after the patients being discharged. Laboratory confirmation, namely IgM serological test was carried out at the respected hospitals. These tests were independently repeated at least twice to confirm the results in our laboratory in addition to IgG serological test. Dengue specific captured IgG and IgM ELISA kit (PanBio Diagnostics, Brisbane, Australia) were used in our laboratory. Diagnosis of dengue was determined by the expert clinicians at each study centre respectively based on the WHO 2009 criteria [2] .

Informed consents were collected from all recruited patients upon hospitalization. Clinically diagnosed patients with positive serological test either for IgG, IgM or both were included into this study. Patients who had been co-infected with other pathogens, or were negative for both IgG and IgM, were excluded.

Complete clinical history, laboratory and other parameters pertaining to dengue infection diagnosis were taken. All personal information and clinical/laboratory results were handled confidentially as mentioned in consent form.

Ethics approval was obtained earlier from the Ethics Committee 

Since dengue patients may show a vast spectrum of clinical manifestations, we decided to select vascular permeability as the trait of interest in this study.

Dengue patients who manifested a significant increase in hematocrit (Hct) and/or any sign of vascular permeability (eg clinical presence of ascites and/or pleural effusion) were classified as cases, while dengue patients who did not show any of these characteristics were classified as controls. According to the Ministry of Health Malaysia [13] , the increase of hematocrit (Hct) was defined as more than 46% and 40% for male and female respectively. This is based on the local normal range of Hct in adults and due to the unavailability of the baseline Hct level in the respective study centres.

Selection of SNPs was carried out by tagSNP selection calculated by Tagger in the HAPLOVIEW [14] , based on the 90 HapMap CHB/ JPT genotype data. Additional SNPs with putative functional variants resulting in changes in protein sequence were also included.

Genomic DNA was extracted either from whole, or clotted blood using commercially available kit, with slight modification, according to Zuraihan et al. [15] .

Genotyping of SNPs was performed using the MassARRAY System (Sequenom Inc, San Diego, USA). Following PCR amplification, primer extension products were analyzed by chip-based MALDI-TOF MS. Extension primers were designed to extend beyond the SNP site by one or two bases. Primer extension and PCR were performed according to the manufacturer's instructions. The MassEXTEND reaction product was performed using groups of identical termination mixtures provided by the manufacturer. Desalting method was carried out and products were loaded into a SpectroCHIP (Sequenom Inc) preloaded with patches of crystalline matrix. The SpectroCHIPs were analyzed with MALDI-TOF MassArray System.

The genotype and allele frequencies of each variant were determined. Hardy-Weinberg Equilibrium (HWE) was performed to compute the deviation of the variants tested. Fisher's exact test was used to determine the significance of the difference in genotype distributions.

Chi-square test with Yates' correction and/or Fisher's exact test (when applicable) was performed on the polymorphic variants for association with vascular permeability of dengue using GraphPad QuickCals (http://graphpad.com/quickcalcs/contingency1/). Twosided P values were calculated. Odds ratio (OR) and a Cornfield's 95% confidence interval (95% CI) were calculated.

Analyses of haplotypes were carried out using HAPLOVIEW [14] . Haplotype blocks were constructed according to definition by Gabriel et al. [16] . Haplotype blocks and the haplotype frequencies were estimated. Haplotype association tests were performed using chi square test in HAPLOVIEW. Permutation test (1000Â) was performed to examine the association significance.

Out of the 176 recruited subjects, 162 dengue patients were included into this study. This includes 103 cases and 59 controls. The demographic and clinical characteristics of the subjects are listed in Table 1 . Since the recruited samples were hospitalized patients, they were expected to have more severe clinical manifestations, therefore more cases were recruited than the controls. In most cases, subjects clinically diagnosed as dengue with warning sign (DW) based on WHO2009 presented with significantly increased of Hct (according to the definition by Ministry of Health Malaysia). Subjects clinically diagnosed as DW but did not show marked increased of Hct were excluded from this analysis to minimize potentially confounding factors to the genetic finding. On the other hand, all patients diagnosed as uncomplicated DF based on WHO2009 did not show increased of Hct therefore, all controls in this study were clinically diagnosed as uncomplicated DF. Sample classifications were blinded prior to genotyping experiment to avoid any potential bias. Fourteen samples were excluded from further analysis due to various factors including co-infections (such as leptospirosis and H1N1), incomplete clinical data and poor genotyping quality. Mean Hct was found to be significantly different between cases and controls (P < 0.0001), but no significant difference was observed in platelet count, ethnicity, gender, age and infection status.

A total of 25 SNPs were selected in this study, of which 14 were polymorphic. For all tested coding SNPs the MAF was = 0 in our dengue cohort. All variants did not deviate from the Hardy-Weinberg Equilibrium (P > 0.05). Allele frequencies obtained were compared with the datasets from 1000Genomes [17] , as shown in Table S1 .

CXCR3 SNPs were analyzed separately for gender since this gene is hemizygous, ie located in chromosome X. Fisher exact test revealed both CXCR3 SNPs were not significantly associated with vascular permeability (Table 2) .

Two out of the four polymorphic SNPs from CXCL10 were significantly associated with case group namely, rs4859584 (P = 0.0230; OR = 2.0; 95% CI, 1.1343-3.5262) and rs8878 (P = 0.0460; OR = 1.8; 95% CI, 1.0503-3.3176) ( Table 3) . These variants showed a slightly lower r 2 value with the tagSNP rs3921 (0.938 and 0.979 respectively) (Fig. 1) . CXCL11 trended towards the association with vascular permeability (P = 0.0527; OR = 1.86; 95% CI, 1.0338-3.3162) ( Table 4) (Fig. 2) .

However, genotypes of all candidate genes did not show any significant association with vascular permeability. Since CXCL10 and CXCL11 are physically located close to each other, a combined LD analysis was performed (Table 5 ). Analysis revealed that variants of these genes are in strong LD. Slight lower D 0 estimation was observed between rs9994667 and rs8878 (D 0 = 0.9) (Fig. 3) . This is in accordance with the HapMap CHB LD block.

Two haplotypes were formed in CXCL10; while three were formed in CXCL11, each with the frequencies as indicated in Figs. S1 and S2. Fisher's exact were performed on haplotype associations of genes CXCL10 and CXCL11. Haplotype ''CCCCA'' (CXCL10) was found to be significantly associated with the case (P = 0.0154); while ''AGTTTAC'' (CXCL11) was found to be significantly associated with the control (P = 0.0366) ( Table 5 ). Haplotype ''CCCCA'' (CXCL10) remained significant after permutation test (1000 times) was performed (permutation P value = 0.0340; OR = 2; 95% CI 1.1343-3.5262).

Haplotypes generated from the combination of the two chemokine CXCL10 and CXCL11 genes revealed three haplotypes with frequencies of 0.173, 0.302 and 0.513 respectively (Fig. S4) . Association analysis revealed that haplotype ''TGGAGAGTTTAC'' was significantly associated with vascular permeability of dengue (P = 0.0437; OR of 0.522, 95% CI [0.3086-0.9882]) ( Table 5 ).

CXCL10 (IP-10) and CXCL11 (I-TAC) are both IFN-gamma induced chemokines. The production of these inflammatory chemokines leads to the recruitment of CXCR3 expressing T cells and NK cells to the site of infection or inflammation [18, 19] . Increased CXCL10, CXCL11 and CXCR3 expression has been implicated in number viral infections including viral meningitis and dengue [9, 20, 21] , inflammatory and autoimmune diseases [22] [23] [24] . However, the impact of genetic variations of these genes in infectious diseases in particular dengue has yet to be widely investigated.

CXCL10 (IP-10) is believed to compete with dengue virus for binding to cell surface heparan sulfate, consequently reduces the DENV uptake and infection of cells [25] [26] [27] . CXCL10 variant namely rs8878, located at the 3 0 -UTR, is believed to affect RNA-binding protein important for efficient translation of this gene via RNA stabilization [28] . It has been associated with autoimmune diseases like type-I diabetes [29] , and severe acute respiratory syndrome amongst Chinese [30] . As for another SNP with significant association, rs4859584, a search in dbSNP though revealed no known disease association was reported [31] . Function prediction analysis performed using FastSNP [32] suggested that both rs8878 and rs4859584 might play a role as an intronic enhancer.

CXCL11 is known to be involved in the inflammatory process of HCV infection, and its expression is upregulated in chronic hepatitis C infection (CHC) [33] . However, in the current study, all variants only reached a moderate significance level, most likely due to small number of samples. Analysis of haplotypes from a LD block has been a powerful tool to localize the candidate region(s) underlying complex diseases as it increases the statistical power to assess an association over an individual marker [34] . In this study, haplotype analysis further assured the influence of CXCL10 and CXCL11 in susceptibility of vascular permeability in dengue.

Further functional and epidemiological studies involving larger number of samples and markers are crucial to elucidate the underlying mechanisms of these genes, in particular associated variants and to define the magnitude of the genetic factor on the vascular permeability of dengue.

In summary, to the best of our knowledge, this study marks the first report on the association of common variants of chemokine genes CXCL10 and CXCL11 with vascular permeability of dengue infection. The finding of this study would shed light to the understanding of the pathogenesis of vascular leakage in dengue infection, which would be invaluable in improving diagnosis, treatment and prevention of the disease.

World Health Organization. Dengue, Dengue Hemorrhagic Fever and Dengue Shock Syndrome in the Context of the Integrated Management of Childhood Illness. Geneva: WHO

Dengue: Guidelines for Diagnosis, Treatment, Prevention and Control

Original antigenic sin in dengue

Original antigenic sin and apoptosis in the pathogenesis of dengue hemorrhagic fever

Dengue hemorrhagic fever with special emphasis on immunopathogenesis

Dengue: defining protective versus pathologic immunity

Race: a risk factor for dengue hemorrhagic fever

Host gene expression profiling of dengue virus infection in cell lines and patients

Both CXCR3 and CXCL10/IFN-Inducible protein 10 are required for resistance to primary infection by dengue virus

Peripheral CXCR3-associated chemokines as biomarkers of fibrosis in chronic hepatitis C virus infection

CXCR3 mediates regionspecific antiviral T cell trafficking within the central nervous system during West Nile virus encephalitis

Clinical Practice Guideline: Management of Dengue Infection in Adult. Malaysia: second ed

Haploview: analysis and visualization of LD and haplotype maps

An alternate method for DNA and RNA extraction from clotted blood

The structure of haplotype blocks in the human genome

Chemokine receptor specific for IP10 and mig: structure, function, and expression in activated T-lymphocytes

Interferon-inducible T cell alpha chemoattractant (I-TAC): a novel non-ELR CXC chemokine with potent activity on activated T cells through selective high affinity binding to CXCR3

Chemotactic activity on mononuclear cells in the cerebrospinal fluid of patients with viral meningitis is mediated by interferon-gamma inducible protein-10 and monocyte chemotactic protein-1

Expression and agonist responsiveness of CXCR3 variants in human T lymphocytes

The chemokine receptors CXCR3 and CCR5 mark subsets of T cells associated with certain inflammatory reactions

Differential expression of three T lymphocyte-activating CXC chemokines by human atheroma-associated cells

Expression of specific chemokines and chemokine receptors in the central nervous system of multiple sclerosis patients

Dengue virus induces expression of CXC chemokine ligand 10/IFN-c-inducible protein 10, which competitively inhibits viral binding to cell surface heparan sulfate

Resistance to dengue virus infection in mice is potentiated by CXCL10 and is independent of CXCL10-mediated leukocyte recruitment

Susceptibility and response of human blood monocyte subsets to primary dengue virus infection

Interferon-gamma-inducible protein (IP)-10 mRNA stabilized by RNA-binding proteins in monocytes treated with S100b

Effect of the IP10 (CXCL10) and HLA genotype on the risk of type 1 diabetes in children

The association of RANTES polymorphism with severe acute respiratory syndrome in Hong Kong and Beijing Chinese

Database for Short Genetic Variation

Expression of the CXCR3 ligand I-TAC by hepatocytes in chronic hepatitis C and its correlation with hepatic inflammation

Haplotype analysis in population genetics and association studies

HBP and SAB conceptualized the study; HBP prepared the manuscript; USH, ZZ, ZMZ performed the laboratory experiments and data analysis; RHS, MM, NKNY collected and diagnosed samples from Kota Bharu, BS, MZ and CL collected and diagnosed samples from Hospital Sungai Buloh.

Supplementary data associated with this article can be found, in the online version, at http://dx.doi.org/10.1016/j.humimm.2015. 03.019.