key: cord-0897492-muf6mqoy
authors: Zhang, Xiaoai; Xu, Hongmei; Chen, Xiaodan; Li, Xiujun; Wang, Xianjun; Ding, Shujun; Zhang, Renli; Liu, Lijuan; He, Cui; Zhuang, Lu; Li, Hao; Zhang, Panhe; Yang, Hong; Li, Tingyu; Liu, Wei; Cao, Wuchun
title: Association of functional polymorphisms in the MxA gene with susceptibility to enterovirus 71 infection
date: 2013-10-02
journal: Hum Genet
DOI: 10.1007/s00439-013-1367-3
sha: 745215601ff67cccdeb7294ae116243da8a8b29c
doc_id: 897492
cord_uid: muf6mqoy

Myxovirus resistance A (MxA) is an antiviral protein induced by type I interferons α and β (IFN-α and IFN-β) that can inhibit virus replication. We examined whether the MxA polymorphisms were related to the risk and severity of enterovirus 71 (EV71) infection in Chinese populations. The MxA C-123A and G-88T polymorphisms were genotyped in two independent case–control populations in China by polymerase chain reaction-based restriction fragment length polymorphism (PCR-RFLP) analysis. Multivariate logistic regression analysis was used to calculate adjusted odds ratios (ORs) and 95 % confidence intervals (95 % CIs). MxA messenger RNA was quantified by real-time quantitative PCR in peripheral blood mononuclear cells (PBMCs) from 45 healthy children and 19 patients with EV71 infection. Significantly decreased susceptibility to EV71 infection was observed for the -123A allele and -88T allele carriers, with ORs (95 % CIs) estimated as 0.56 (0.39–0.81) and 0.64 (0.47–0.88), respectively, in the northern population. This association was confirmed in the southern population, with ORs (95 % CIs) estimated as 0.58 (0.38–0.89) and 0.67(0.47–0.95), respectively. The A(-) (123)T(-) (88) haplotype was also significantly associated with lower risk of EV71 infection in both the northern (OR = 0.62; 95 % CI = 0.44–0.85) and the southern population (OR = 0.63; 95 % CI = 0.43–0.92). Furthermore, we observed higher MxA messenger RNA levels in IFNβ1a-stimulated PBMCs from the -123A or -88T allele carriers compared with that from nocarriers. Our findings suggest that polymorphisms in the MxA promoter may play a role in mediating the susceptibility to EV71 infection in Chinese population.

Enterovirus 71 (EV71) is positive-stranded RNA virus belonging to the enterovirus genus of the Picornauiridae Abstract Myxovirus resistance A (MxA) is an antiviral protein induced by type I interferons α and β (IFN-α and IFN-β) that can inhibit virus replication. We examined whether the MxA polymorphisms were related to the risk and severity of enterovirus 71 (EV71) infection in Chinese populations. The MxA C-123A and G-88T polymorphisms were genotyped in two independent case-control populations in China by polymerase chain reaction-based restriction fragment length polymorphism (PCR-RFLP) analysis. Multivariate logistic regression analysis was used to calculate adjusted odds ratios (ORs) and 95 % confidence intervals (95 % CIs). MxA messenger RNA was quantified by real-time quantitative PCR in peripheral blood mononuclear cells (PBMCs) from 45 healthy children and 19 patients with EV71 infection. Significantly decreased susceptibility to EV71 infection was observed for the -123A allele and -88T allele carriers, with ORs (95 % CIs) estimated as 0.56 (0.39-0.81) and 0.64 (0.47-0.88), respectively, in the 1 3 family. EV71 has been regarded as the most important neurotropic enterovirus after the primary eradication of the poliovirus. There is no specific antiviral therapy or vaccine for EV71. Since its discovery in 1969, EV71 outbreaks have been reported periodically throughout the world. The size and frequency of EV71 outbreaks have greatly increased over the past decade in the Asia-Pacific region, especially in southeast Asia (Solomon et al. 2010) . In mainland China, EV71 infection has become more and more severe and led to 509 and 569 deaths in 2011 and 2012, respectively, according to the report of Ministry of Health, China.

Individual susceptibility to EV71 infection seems to be variable (Chang et al. 2008; Yang et al. 2001 Yang et al. , 2012 . Patients with EV71 infection displayed wide clinical spectrum, ranging from asymptomatic or mild hand-foot-mouth disease (HFMD) to severe neurological disease. Similar to other infectious diseases, EV71 infection may also be the result of complex interactions between pathogen, environmental factors and host genes. The identification of susceptibility genes contributing to EV71 infection would assist in predicting individual and population risks of EV71 infection, as well as clarifying pathogenesis of this infectious disease.

Type I interferons α and β (IFN-α and IFN-β) system is a key component of the innate immune response and the first line of defense against viral infection including EV71 infection (Liu et al. 2005; Yi et al. 2011 ). The Myxovirus resistance (Mx) proteins are dynamin-like GTPases that are induced by IFN-α and IFN-β and have antiviral activity against RNA viruses (Sadler and Williams 2008) . In humans, two distinct Mx (MxA and MxB) are encoded on chromosome 21, while only MxA has exhibited antiviral activity (Sadler and Williams 2008) . Transfected cells (Chieux et al. 2001; Kochs and Haller 1999) and transgenic mice (Pavlovic et al. 1995) overexpressing MxA acquire a high degree of antiviral resistance against several RNA viruses. In humans, synthesis of MxA protein is induced during viral infections and may protect humans from infection (Roers et al. 1994; von Wussow et al. 1990 ). As for EV71, it has been shown that type I IFN subtypes can inhibit EV71 replication and the antiviral effect was caused by a strong induction of the downstream antiviral effectors especially MxA (Yi et al. 2011 ). On the basis of the functional role of MxA in the pathogenesis of viral infections in vivo and in vitro, we hypothesize that MxA might be a potential candidate susceptibility gene for EV71 infection.

Among the validated polymorphisms in MxA gene, two functional single nucleotide polymorphisms (SNPs) (C-123A, rs17000900; G-88T, rs2071430) in promoter region have been extensively studied. The electrophoretic mobility shift assays shown that these two polymorphisms could influence the binding of nuclear protein(s) (Ching et al. 2010) . Furthermore, by luciferase reporter assay, the basal promoter activity was demonstrated to be higher in promoter constructs containing the minor alleles of the two polymorphisms (Furuyama et al. 2006; Hijikata et al. 2001; Torisu et al. 2004 ). We and other investigators have reported the association between the aforementioned MxA functional SNPs and risk of infections of viruses, including severe acute respiratory syndrome-coronavirus (SARS-CoV) (Ching et al. 2010; Hamano et al. 2005; He et al. 2006) , hepatitis B virus (HBV) (Cao et al. 2009 ), hepatitis C virus (HCV) (Hijikata et al. 2000) and measles virus (Torisu et al. 2004) . The role of the MxA polymorphisms in EV71 infection, however, has never been investigated. In the present study, we examined whether polymorphisms in the MxA gene were related to the risk and severity of EV71 infection among Chinese populations.

Patients were recruited between April and November 2010 from 7 hospitals located in either northern or southern China. All the hospitals were surveillance sentinels which were set by China CDC. HFMD was clinically defined as a patient having oral ulcers and vesicular rash on the hands, feet, knees, or buttocks. Complicated patients were defined as cases with central nervous system (CNS) involvement (including encephalitis, poliomyelitis-like syndrome, and encephalomyelitis), myocarditis or pulmonary edema/hemorrhage. In each hospital, all patients who met the clinical definition and willing to participate into the study were recruited.

Stool samples were collected for EV71 RNA detection by real-time PCR using previously described primers (Verstrepen et al. 2002) . Viral isolation was performed as described previously (Chen et al. 2010) . Serum samples were collected to test the IgM antibody against EV71 by a commercial enzyme-linked immunosorbent assay (ELISA) kit (Beijing Wantai Biological Pharmacy Enterprise Co, Ltd, China). Confirmation of EV71 infections was recognized by one of the following criteria: isolation of EV71 from stool specimens; positive EV71 viral RNA quantitation or positive EV71 IgM results during acute phase of infection.

The controls were randomly selected from healthy children who underwent routine physical examination in the same region during the same period that the cases were recruited. After interviewing their guardians, healthy children with no self-reported or medical-recorded history of HFMD or herpangina were recruited. Their sera samples were collected for EV71-specific IgG antibody test by the Diagnostic Kit for IgG Antibody to Human EV71 (Beijing Beier Bioengineering Co, Ltd, China). Only those with negative EV71 IgG results were recruited as eligible controls. By checking the medical records or by interviewing the participants' guardians, we determined that all cases and controls were genetically unrelated Han Chinese.

For MxA expression assay, blood samples were obtained from a group of 20 healthy children (13 males and 7 females) aged 26-50 months (mean ± SD, 39.6 ± 6.6) from the northern population, 25 healthy children (16 males and 9 females) aged 24-52 months (mean ± SD, 38.2 ± 8.4) from southern population, and 19 patients with EV71 infection (12 males and 7 females) aged 24-50 months (mean ± SD, 37.8 ± 7.8). At recruitment, informed consent was obtained from all participants' guardians, and personal information on demographic factors and medical history were collected via structured questionnaire. This study was performed with the approval of the Ethical Committee of Beijing Institute of Microbiology and Epidemiology and conducted according to the principles expressed in the Declaration of Helsinki.

Genomic DNA was extracted from peripheral blood by Gene Relax DNA extraction kit (TianGen, Beijing, China) according to the manufacturer's instructions. The polymorphisms were genotyped using polymerase chain reaction-based restriction fragment length polymorphism (PCR-RFLP) analysis. The primers 5′-GGGAAGGACAT GCCTAGGTTC-3′ and 5′-CCGCGAAGAAATGAAACT CAC-3′ were used for amplifying the target region. The a PCR-based RFLP analysis for the C-123A polymorphism. Lanes 2-7 were CC genotypes; lane 1 was CA genotype; lane 8 was AA genotype; M was DL2000 marker. b PCR-based RFLP analysis for the G-88T polymorphism. Lanes 4 and 8 were GG genotypes; lanes 1, 2, and 5-7 were GT genotypes; lane 3 was TT genotype; M was DL2000 marker. c Sequencing analysis for genotypes of the C-123A polymorphism. The 3 charts represent the CC, CA and AA genotypes, respectively. d Sequencing analysis for genotypes of the G-88T polymorphism. The 3 charts represent the GG, GT and TT genotypes, respectively. The arrows localize the base changes at the nucleotide positions 1 3 PCR conditions were performed as described previously except for an annealing temperature of 57.5 °C (Zhang et al. 2012) . The reaction yielded a 263-base pair amplicon. PCR products of 3 μL were then digested with 4U of PstI and 4U HhaI (Takara BioTech, Dalian, China) for the C-123A and G-88T polymorphisms, respectively. Genotyping was done by staff blinded to the subjects' case/control status. The accuracy of genotyping data obtained from PCR-RFLP analyses was validated by direct sequencing of a 15 % masked, randomly selected sample of cases and controls, and the results were 100 % concordant ( Fig. 1 ).

Real-time analysis of MxA mRNA Total RNA of peripheral blood mononuclear cells (PBMCs) from patients with EV71 infection and IFNβ1a (1,000 U/mL, R&D Systems) stimulated PBMCs from healthy children were extracted using Trizol reagent (Invitrogen, USA). Reverse transcription was done using reverse transcriptase system (Promega, USA) according to the manufacturer's instructions. Quantitative realtime PCR for MxA, with GAPDH as an internal reference gene, was performed using the Roche Light Cycler 2.0 (Roche Applied Science, USA) based on the SYBR-Green method. The primers used for MxA (Fernandez-Arcas et al. 2004 ) were 5′-ACAGCAAATCAAGGCACTGGA-3′ and 5′-CGGATCAGCTTCTCACCTTCTC-3′and for GAPDH (Xie et al. 2011) were 5′-GGGAAGGTGAAGGTCGG AGT-3′ and 5′-TTGAGGTCAATGAAGGGGTCA-3′. The PCR reaction mixture consisted of 0.1 μmol/L of each primer, 1 X SYBR Green Master mix (Takara), and 50 ng complementary DNA to a final volume of 20 μL. Cycling conditions were 95 °C for 5 min, followed by 40 cycles at 95 °C for 15 s and 57 °C for 15 s. Data were analyzed using the 2 −ΔΔCt method (Livak and Schmittgen 2001) . All analyses were performed in a blinded fashion with the laboratory persons unaware of genotyping data.

Genotype and allele frequencies for each polymorphism were determined by gene counting. The fitness to Hardy-Weinberg equilibrium was tested using the random-permutation procedure implemented in the Arlequin package (http://lgb.unige.ch/arlequin). Associations between polymorphisms and risk of EV71 infection were estimated by use of logistic regression analyses. Odds ratios (ORs) and 95 % confidence intervals (CIs) were used to measure the strength of association. The Mantel-Haenszel (MH) test was performed for the combined analysis of the data from both populations. In view of the multiple comparisons, the correction factor n (m − 1) (n loci with m alleles each) was applied to correct the significance level. These analyses were performed using SPSS software (version 15.0). The pairwise linkage disequilibrium (LD) index (Lewontin's D' and r 2 ) calculation was performed using the Arlequin package. The haplotype analysis was performed using PLINK software (version 1.07) (http://pngu.mgh.harvard.edu/ purcell/plink/) (Purcell et al. 2007) .

A total of 441 laboratory-confirmed patients with EV71 infection and 958 controls were included in the study, including 253 cases and 550 controls from northern China, and 188 cases and 408 controls from southern China. In both populations, controls were comparable with cases with regard to mean age and sex distribution. A total of 150 (59.3 %) patients of northern origin and 151 (80.3 %) cases of southern origin were determined as complicated patients. The characteristics of selected subjects are shown in Table 1 .

The genotyping results were presented in Table 2 . The observed genotype frequencies of the C-123A and G-88T polymorphisms were in Hardy-Weinberg equilibrium (all P > 0.05, data not shown). The -123A-positive genotypes were significantly associated with decreased susceptibility to EV71 infection either in the northern population (OR, Table 2) . Similarly, the -88T-positive genotypes were also significantly associated with decreased susceptibility to EV71 infection with an OR of 0.64 (95 % CI, 0.47-0.87; P = 0.004) in the northern population and 0.67 (95 % CI, 0.47-0.95; P = 0.02) in the southern population (Table 2) . Then multivariate logistic regression analyses were performed by adjusting for age and sex. For the C-123A polymorphism, subjects bearing at least one -123A allele had a decreased susceptibility to EV71 infection compared with those without -123A allele in both independent populations (in the northern population: adjusted OR, 0.56; 95 % CI, 0.39-0.81; P = 0.002; in the southern population: adjusted OR, 0.58; 95 % CI, 0.38-0.89; P = 0.01) ( Table 2) . For the G-88T polymorphism, subjects bearing at least one -88T allele had a decreased susceptibility to EV71 infection compared with those without -88T allele in both independent populations (in the northern population: adjusted OR, 0.64; 95 % CI, 0.47-0.88; P = 0.005; in the southern population: adjusted OR, 0.67; 95 % CI, 0.47-0.95; P = 0.02) ( Table 2 ). The associations remained significant after correction for multiple comparisons.

When northern and southern population were combined, both the -123A-positive and -88T-positive genotypes were significantly associated with decreased susceptibility to EV71 infection, for the -123A-positive Table 2 ). There was no significant heterogeneity for both the -123A-positive and -88T-positive genotypes associated ORs among the study populations (Breslow-Day P = 0.87 for C-123A; Breslow-Day P = 0.85 for G-88T).

The effects of C-123A and G-88T genotypes on clinical severity of EV71 infection were also assessed. In both populations, there was no association between G-88T polymorphism and severity of EV71 infection. For C-123A polymorphism, no association was observed with severity of EV71 infection in the southern population. While in the northern population, subjects bearing at least one -123A allele had a decreased susceptibility to EV71 infection (Table 3) .

By LD analysis, the D' values for the MxA C-123A and G-88T polymorphisms in the northern and southern populations were both 1.0, whereas r 2 was 0.48 and 0.44, respectively. Three haplotypes (C -123 -G -88 , C -123 -T -88 and A -123 -T -88 ) were observed based on these two polymorphisms. Consistent with the findings of published studies (Ching et al. 2010; Weinstock-Guttman et al. 2007) , the A -123 -G -88 haplotype was absent in our two populations. The overall haplotype analysis gave a P value of 0.01 in the northern and 0.02 in the southern population (Table 4) . Compared with all other haplotypes, the A -123 -T -88 was associated with a significantly lower risk to EV71 infection in the two independent populations (in the northern population: OR, 0.62; 95 % CI, 0.44-0.85; P = 0.003; in the southern population: OR, 0.63; 95 % CI, 0.43-0.92; P = 0.01) as well as in the combined population (OR, 0.62; 95 % CI, 0.48-0.80; P < 0.001) ( Table 4 ). There was a significantly association between the MxA A -123 -T -88 haplotype and severity of EV71 infection in the northern population (OR, 0.44; 95 % CI, 0.25-0.79; P = 0.005), which however failed to be detected in the southern population (Table 4) .

Quantitative real-time RT-PCR for MxA expression As shown in Fig. 2 , the relative MxA mRNA level from individuals with -123A allele (AA + AC genotypes) was significantly higher than that from individuals with -123CC genotype in both northern and southern populations (t test; P < 0.05). With regard to the G-88T polymorphism, subjects with the -88T allele (TT + GT genotypes) appeared to express higher MxA mRNA than those with the -88GG genotype in both northern and southern populations (t test; P < 0.05). In addition, those subjects that are homozygous for the major C-123A and G-88T haplotype (CC/GG) had lower levels of MxA mRNA than subjects that are heterozygous (CA/GT) in both populations (t test; P < 0.05). The relative MxA mRNA level was also compared between uncomplicated and complicated cases (Fig. 3) . There were no significant differences in the relative MxA mRNA level between uncomplicated and complicated cases in both populations (t test; P > 0.05). 

In this study, we found that two functional polymorphisms in the promoter of MxA gene, C-123A and G-88T, were significantly associated with decreased susceptibility to EV71 infection in two independent Chinese subpopulations.Haplotype A -123 -T -88 , derived from these two polymorphisms was also shown to be associated with the decreased risk of EV71 infection. Furthermore, consistent with the results of population-based association study, higher MxA messenger RNA levels in IFNβ1a-stimulated PBMCs from the -123A or -88T allele carriers were observed. To the best of our knowledge, this study is the first to report a genetic association between MxA and susceptibility to EV71 infection, confirming the initial hypothesis that MxA may play a role in the pathogenesis of this disorder. The genetic association between MxA polymorphism and EV71 infection is biologically plausible. MxA is an antiviral protein induced by type I IFNs that can inhibit viral replication. Following stimulation with type I IFNs, MxA protein accumulates to high levels in the cytoplasm on intracellular membranes. At the point of infection, MxA are released and bind viral nucleocapsids or other viral components, to trap and then degrade them (Sadler and Williams 2008) . With regard to EV71, it has been shown that the superior anti-EV71 effect of type I IFN was caused by a strong induction of the downstream antiviral effectors, especially MxA (Sadler and Williams 2008) . In addition, Fig. 2 Quantitative real-time RT-PCR for MxA in healthy children. Expression of MxA for three different genotypes of C-123A and G-88T was measured in RNA from 45 healthy children (20 from the northern population and 25 from the southern population) stimulated with IFNβ1a. Normalization for mRNA quantity was performed with human GAPDH control primers for each sample and final abundance figures adjusted to yield an arbitrary value of 1 for -123CC or -88GG genotype. Columns mean from triplicate measurements; the vertical bars indicate the standard deviation (SD). a Correlation of MxA mRNA expression with C-123A genotypes in the northern population. Compared to the CC carriers, the A allele carriers had a markedly elevated MxA transcription (P < 0.05; t test). b Correlation of MxA mRNA expression with G-88T genotypes in the northern population. Compared to the GG carriers, the T allele carriers had a markedly elevated MxA transcription (P < 0.05; t test). c Correlation of MxA mRNA expression with MxA haplotypes in the northern population. The subjects that are homozygous for the major C-123A and G-88T haplotype (CC/GG) had lower levels of MxA mRNA than subjects that are heterozygous (CA/GT) (t test; P < 0.05). d Correlation of MxA mRNA expression with C-123A genotypes in the southern population. Compared to the CC carriers, the A allele carriers had a markedly elevated MxA transcription (P < 0.05; t test). e Correlation of MxA mRNA expression with G-88T genotypes in the southern population. Compared to the GG carriers, the T allele carriers had a markedly elevated MxA transcription (P < 0.05; t test). f Correlation of MxA mRNA expression with MxA haplotypes in the southern population. The subjects that are homozygous for the major C-123A and G-88T haplotype (CC/GG) had lower levels of MxA mRNA than subjects that are heterozygous (CA/GT) (t test; P < 0.05) it has recently been reported that the antiviral activity of MxA extended to coxsackieviruses B virus, a member of the enterovirus family (Chieux et al. 2001) . Thus, it is highly possible that MxA can bind to the components of EV71 and exert its role as an antivirus, although additional studies are needed to clarify this possibility.

The MxA C-123A and G-88T polymorphisms are both functional polymorphisms. Previous in vitro studies showed that -123A and -88T alleles provided stronger binding affinity to nuclear proteins than the wild-type alleles (Ching et al. 2010) . Furthermore, the C-123A and G-88T polymorphisms are in strong LD and the A -123 -T -88 haplotype displays higher promoter activity compared with the C -123 -G -88 haplotype in luciferase reporter assay (Furuyama et al. 2006; Hijikata et al. 2001; Torisu et al. 2004 ). In the present study, we also detected higher MxA transcription in IFNβ1a-stimulated PBMCs from the -123A or -88T allele carriers. Given the well-characterized antiviral role of MxA, one might expect that individuals who carry the -123A and/or -88T alleles, and thus have increased expression of MxA, may be at a lower susceptibility to EV71 infection.

In this study, the association with risk of EV71 infection was stronger for the C-123A polymorphism than the G-88T polymorphism. This finding is biologically rational. Although administration of type I IFN can inhibit EV71 replication in vitro and in mouse model, a recent study found that EV71 inhibits the expression of IFN-β in infected cells (Lei et al. 2010) , which suggests that the virus evades the antiviral activity of IFN-β. As the C-123A polymorphism affects basal MxA expression without IFN-α and IFN-β induction, whereas the G-88T polymorphism influenced MxA expression after IFN-β stimulation, the C-123A SNP probably exerts greater effect on basal MxA expression in innate immune response against EV71, in which endogenous IFN induction is suppressed. In addition, we found markedly elevated MxA transcription in response to IFNβ1a stimulation for genotypes of the G-88T polymorphism, while the response is not so marked for the C-123A polymorphism. Our results are also consistent with previous findings that the C-123A polymorphism plays a more important role in modifying basal MxA expression, and contributes more significantly to the innate immune response against SARS (Ching et al. 2010) .

In this study, we observe associations between MxA C-123A polymorphism and A -123 -T -88 haplotype with disease severity only in northern population. This result may be due to the limited sample size of patients with uncomplicated cases (n = 37) in the southern population in our study. However, no association between C-123A polymorphism and severity of EV71 infection was observed when northern and southern populations were combined for analysis. Moreover, there were no significant differences in the relative MxA mRNA level between uncomplicated and complicated cases. Therefore, our results indicate that MxA C-123A polymorphism might not play a role in disease severity of EV71 infection. Still, this result warrants further confirmation in future studies.

To our knowledge, only three genetic association studies have addressed susceptibility genes that related to EV71 infection (Chang et al. 2008; Yang et al. 2001 Yang et al. , 2012 . Some of the results, however, could not be replicated in patients with EV71 infection (12 from the northern population and 7 from the southern population). Normalization for mRNA quantity was performed with human GAPDH control primers for each sample and final abundance figures adjusted to yield an arbitrary value of 1 for complicated cases. Columns mean from triplicate measurements; the vertical bars indicate the standard deviation (SD). a MxA mRNA expression in patients with EV71 infection in the northern population. Compared to the complicated cases, the uncomplicated cases had a comparable MxA transcription level (P > 0.05; t test). b MxA mRNA expression in patients with EV71 infection in the southern population. Compared to the complicated cases, the uncomplicated cases had a comparable MxA transcription level (P > 0.05; t test) subsequent studies in other populations. The inconsistent results could be attributable to the small sample size and/or inadequate adjustment for confounding factors such as age and sex. Age and sex are well-known confounding factors for EV71 infection, because individuals aged 1-4 years are associated with a higher risk (Chan et al. 2003; Mao et al. 2010) , and the proportion of male was higher in patients with EV71 infection (Chen et al. 2007 ). This tendency was also identified in our case patients. In the present study, the association was derived from two different subpopulations of large sample size, with potential confounding factors (sex and age) adjusted. Our study, therefore, could probably pose as one of the largest case-control studies demonstrating valid and replicable association of EV71 infection.

Recently, several association studies have shown that the MxA gene was related to the susceptibility to various specific infections (Cao et al. 2009; Ching et al. 2010; Hamano et al. 2005; He et al. 2006; Hijikata et al. 2000) . Some of the results, however, could not be replicated in subsequent studies because of multiple factors, such as the different ethnicities of study populations and/or different pathogenesis of the studied diseases. Although the highly significant association between MxA and susceptibility to EV71 derived from a biologically based a priori hypothesis, our findings should be independently verified in populations of different ancestry, such as white subjects and African subjects.

In conclusion, our study shows an association between the MxA polymorphisms and decreased susceptibility to EV71 infection in two Chinese subpopulations. These findings may provide support for the importance of MxA in the pathogenesis of EV71 infection. If confirmed by other studies, knowledge of genetic factors contributing to the pathogenesis of EV71 infection as presented here would be important for the assessment of host susceptibility to EV71 infection and other infectious diseases, especially those sharing a mode of action similar to that of EV71 infection.

The haplotype of the MxA gene promoter is associated with hepatitis B virus infection in a Chinese population

Epidemic hand, foot and mouth disease caused by human enterovirus 71

HLA-A33 is associated with susceptibility to enterovirus 71 infection

An eightyear study of epidemiologic features of enterovirus 71 infection in Taiwan

Comparison of clinical features between coxsackievirus A2 and enterovirus 71 during the enterovirus outbreak in Taiwan, 2008: a children's hospital experience

Inhibition of coxsackievirus B4 replication in stably transfected cells expressing human MxA protein

Significance of the myxovirus resistance A (MxA) gene −123C >a single-nucleotide polymorphism in suppressed interferon beta induction of severe acute respiratory syndrome coronavirus infection

Differential transcriptional expresion of the polymorphic myxovirus resistance protein A in response to interferon-alpha treatment

Single nucleotide polymorphisms and functional analysis of MxA promoter region in multiple sclerosis

Polymorphisms of interferon-inducible genes OAS-1 and MxA associated with SARS in the Vietnamese population

Association of SARS susceptibility with single nucleic acid polymorphisms of OAS1 and MxA genes: a case-control study

Identification of a single nucleotide polymorphism in the MxA gene promoter (G/T at nt −88) correlated with the response of hepatitis C patients to interferon

Genetic polymorphism of the MxA gene promoter and interferon responsiveness of hepatitis C patients: revisited by analyzing two SNP sites (−123 and −88) in vivo and in vitro

Interferon-induced human MxA GTPase blocks nuclear import of Thogoto virus nucleocapsids

The 3C protein of enterovirus 71 inhibits retinoid acid-inducible gene I-mediated interferon regulatory factor 3 activation and type I interferon responses

Type I interferons protect mice against enterovirus 71 infection

Analysis of relative gene expression data using real-time quantitative PCR and the 2(-Delta Delta C(T)) Method

Epidemiology of hand, foot, and mouth disease and genotype characterization of Enterovirus 71 in Jiangsu

Enhanced virus resistance of transgenic mice expressing the human MxA protein

PLINK: a tool set for whole-genome association and population-based linkage analyses

MxA gene expression after live virus vaccination: a sensitive marker for endogenous type I interferon

Interferon-inducible antiviral effectors

Virology, epidemiology, pathogenesis, and control of enterovirus 71

Functional MxA promoter polymorphism associated with subacute sclerosing panencephalitis

Evaluation of a rapid real-time RT-PCR assay for detection of enterovirus RNA in cerebrospinal fluid specimens

The interferoninduced Mx-homologous protein in people with symptomatic HIV-1 infection

Pharmacogenetics of MXA SNPs in interferon-beta treated multiple sclerosis patients

Smurf1 ubiquitin ligase targets Kruppel-like factor KLF2 for ubiquitination and degradation in human lung cancer H1299 cells

Association of interleukin 10 and interferon gamma gene polymorphisms with enterovirus 71 encephalitis in patients with hand, foot and mouth disease

Altered cellular but not humoral reactions in children with complicated enterovirus 71 infections in Taiwan

Potent inhibition of human enterovirus 71 replication by type I interferon subtypes

The p73 G4C14-to-A4T14 polymorphism is associated with risk of lung cancer in the Han nationality of North China

Acknowledgments We thank all the tested individuals, their families, and collaborating clinicians for their participation. 

The authors declare that they have no conflict of interest.