Skip to main content

Relationship between polymorphism of receptor SCARB2 gene and clinical severity of enterovirus-71 associated hand-foot-mouth disease



To investigate the relationship between polymorphism of scavenger receptor class B member 2 (SCARB2) gene and clinical severity of enterovirus (EV)-71 associated hand-foot-mouth disease (HFMD).


Among the 100 recruited cases, 56 were in the severe HFMD group (case group) and 44 were in the general HFMD group (control group). By screening functional single nucleotide polymorphisms (SNPs) and hot SNPs, and performing SNP site optimization, some SNP sites of SCARB2 gene were selected for analysis. Genotyping was performed using a MassArray platform. PLINK software was used for statistical processing and analysis of the correlation differences between the mutant genotypes in the severe and general HFMD groups. The relationship between the SNPs and clinical severity of enterovirus (EV)-71 associated HFMD was assessed.


28 SNPs in SCARB2 were selected by site optimization. Then three loci were not in agreement with the minor allele frequency (MAF) in the 1000 Han Chinese in Beijing (CHB) dataset. Another three loci could not be detected. Nine loci were not suitable for further analysis (MAF < 0.01 and Hardy–Weinberg [HWE] P < 0.001). A total of 13 sites were subsequently analyzed. Through Fisher analysis, the frequency of the rs6812193 T allele was 0.134 and 0.034 in the severe and general HFMD groups, respectively (P 0.023 < 0.05, odds ratio [OR] 4.381 > 1). Logistic regression analysis of rs6812193 T alleles between the severe and general HFMD groups, respectively (P 0.023 < 0.05, OR 4.412 > 1, L95 1.210 > 1). Genotype logistic regression analysis of the rs6812193 alleles CT + TT versus CC gave an OR of 4.56 (95% confidence interval [95% CI] 1.22–17.04, P = 0.012).


The rs6812193 T allele was a susceptibility SNP for SHFMD, and the rs6812193 polymorphism might be significantly associated with the susceptibility to EV-71 infection.


Hand, foot, and mouth disease (HFMD) is an infectious disease caused by a variety of enteroviruses which belongs to the small RNA virus family [1]. HFMD is common in children under 5 years of age, and it is mainly manifested as herpes and maculopapules on the hands, feet, mouth, and other areas. A few patients progress rapidly and develop neurogenic pulmonary edema, circulatory disturbance, and even death at 1–5 days after disease onset [2]. The following indicators should alert the clinician of possible deterioration and impending critical type of severe case: persistent high fever, nervous system involvement, abnormal respiratory rate and rhythm, circulatory dysfunction, elevated peripheral white blood cell count, elevated blood glucose, elevated blood lactic acid [3]. HFMD is a global disease with a variety of causes. The main causes of HFMD are enterovirus (EV)-71 and Coxsackievirus A16 infection; however, EV-71 is responsible for most of the severe cases and fatal cases [4].

EV-71 is another important neuroenterophilic virus after poliovirus elimination. An analysis shows that EV-71 was circulating in the Netherlands as early as 1963 [5], but it was first reported in 1969 [6]. It was global distribution and occasionally concentrated outbreak. Since 1997, outbreaks of HFMD caused by EV-71 have occurred in the Asia–Pacific region, such as Malaysia [7], Taiwan [8], Singapore [9], etc. In 1998, the epidemic outbreaks began in Chinese mainland, especially in 2007, Shandong province [10] and 2008 Anhui province [11], which resulted in a large number of severe and dead children. Therefore, in 2008, the Ministry of Health of China listed HFMD in the list of Class C infectious diseases. In recent years, the incidence of HFMD has been 37.01–205.06 per 100,000, with a fatality rate of 6.46–51.00 per 100,000 [3]. HFMD has the highest number of cases and deaths of all Class C infectious diseases in China and represents a serious threat to the health of children.

As a major cause of severe and fatal cases, the pathogenesis of EV-71 has attracted more and more researchers' attention. In clinical work, it is not difficult to find that the severity of clinical symptoms and prognosis of different children with the same infection with EV-71 are significantly different. As the first important portal for the virus to enter the human body, the virus receptor determines the host range of a virus and tissue specificity. The influence of individual differences on the severity of clinical symptoms is worth further studying. Yamayoushi et al. [12, 13] confirmed that SCARB2 is the receptor of all EV-71 strains in cell experiments. The guidelines for the diagnosis and treatment of HFMD issued by the Ministry of Health of China (2018) clearly indicate that SCARB2 is the main receptor of EV-71 virus. Studies have found that people with different genotypes and alleles have different probability of disease and severity [14]. Single nucleotide polymorphisms (SNPs) are the most common form of variation in human genomic DNA. The SNPs of SCARB2 gene have naturally attracted great attention. At present, most studies focus on the relationship between the severity of EV-71 infection and the polymorphism of cytokines such as TNF-ɑ, IL-6, IL-10 [15, 16], chemokine IP-10, MCP-1(CCl2) [17], immune-related factors OAS1, OAS2, OAS3 and MXA [18, 19]. To the best of our knowledge, few studies on the relationship between SCARB2 SNP and EV-71 HFMD have been performed. In this study, 28 SNP sites in SCARB2 were selected as study loci, expecting to further clarify the pathogenesis of EV71 infection, and to provide a strong research basis for the early warning of critical disease and the reduction of case fatality rate.

Materials and methods

Clinical data and sample collection

We recruited 100 children with HFMD admitted to our hospital from April 2018 to October 2020 who were positive after EV-71 nucleic acid test. Diagnosis was based on the “Guidelines for the Diagnosis and Treatment of Hand, Foot, and Mouth Disease” (2018) [3]. According to the occurrence and development process of the disease, HFMD is divided into general HFMD and severe HFMD. General HFMD is usually in the eruption stage, but severe HFMD includes nervous system involvement stage, early cardiopulmonary failure stage, and cardiopulmonary failure stage, based on the degree of danger and heavy. There were 56 cases with severe HFMD (as the case group) and 44 cases with general HFMD (as the control group). This study was approved by the Ethics Committee of Tianjin Second People's Hospital, and informed consent was obtained from the patients’ parents or family members.

HFMD samples were collected before treatment on the day of admission; 3 mL peripheral venous blood was collected and stored at − 80 °C for later use.

SCARB2 SNP selection

Screening of SNPs

1. Screening of functional SNPs.

  1. a.

    In the National Center for Biotechnology Information (NCBI), SNP was searched using the name of the SCARB2 to identify functional SNPs including promoter proxy (upstream variant 2 KB), 5′-untranslated region (UTR), exons (missense, synonymous), 3′-UTR, etc. The relevant optimization parameter was a minor allele frequency (MAF) in Han Chinese in Beijing (CHB) > 0.05, according to the HapMap or 1000 Genomes databases.

  2. b.

    The screened SNPs were used for functional prediction.

  3. c.

    Linkage disequilibrium analysis was performed on the identified SNPs, and the linked sites with R2 = 1 were labeled.

2. Screening of hot SNPs.

  1. a.

    Through Google Scholar, the research literature of candidate gene-related polymorphisms was retrieved and susceptibility SNPs were screened out.

  2. b.

    The screened SNPs were verified using the CHB MAF values.

SNP site optimization method

According to the linkage disequilibrium analysis results, the completely linked loci with R2 = 1 were discarded, and the loci in the promoter region with R2 > 0.8 were retained (Haplotypes in this region are important for gene expression); however, the strong association sites with R2 > 0.8 identified in other regions and literature studies, are meaningless and omitted.

DNA extraction

DNA was extracted from the blood samples using a Tiangen kit. DNA samples were analyzed on a NanoDrop2000, and 1.25% agarose gel electrophoresis was performed. DNA was quantified and transferred to a 96-well plate for storage at − 20 °C for later use.

SNP typing

Primer design and synthesis

Assay Designer 3.1 software was used to design the primers, and the primers were synthesized by the company [BGI Tech Solutions (Beijing Liuhe) Co., Ltd]. The primer sequences are shown in Table 1.

Table 1 Primers used for SNP typing

Primer dilution and extension mix configuration

The single-tube PCR masters were diluted to 100 μM, and deionized water was added to achieve a final PCR master mix concentration of 0.5 μM. The single tube extension primers were diluted to a final concentration of 500 μM. Each primer was diluted to 8 μM, 10 μM, and 15 μM. According to the instructions of the DNA synthesis products, the molecular weight and number of moles, the amount of deionized water to be added were calculated according to the required concentration. According to the molecular weight of the mixed single-tube extension primers, 1 time (< 6300 Da), 1.2 times (6300–7200 Da), and 1.5 times (> 7200 Da) were taken for mixing.

MassArray reactions

PCR amplification was conducted in 5 µL reactions, containing 1.000 µL DNA (20 ng/μl), 1.000 µL PCR primers (500 nM each), 0.100 µL dNTP mix (25 mM each), 0.625 µL PCR buffer (15 mM MgCl2), 0.325 µL MgCl2 (25 mM), 1.850 µL water HPLC grade, and 0.100 µL Taq DNA polymerase (5 U/µL) (Agena Bioscience, San Diego, CA, USA). The PCR conditions were as follows: 94 °C for 5 min, 94 °C for 20 s, 45 cycles of 56 °C for 30 s, 72 °C for 1 min, and a final extension step at 72 °C for 3 min. Remaining unincorporated dNTPs were dephosphorylated and inactivated by treatment with 1 U shrimp alkaline phosphatase at 37 °C for 20 min and then 85 °C for 5 min. Finally, the single base extension reaction mix, including iPLEX Buffer Plus, iPLEX Termination mix, Extension Primers mix, and iPLEX enzyme (Agena Bioscience, San Diego, CA, USA), was added to the PCR amplification products. The single base extension reaction was carried out under the following conditions: 94 °C for 30 s, 40 cycles at 94 °C for 5 s (52 °C for 5 s and 80 °C for 5 s, repeated 5 times per cycle), and a final extension step at 72 °C for 3 min. The samples were spotted on a SpectroCHIP (Agena Bioscience, San Diego, CA, USA), and analyzed by mass spectrometry. The spectral profiles generated by matrix-assisted laser desorption/ionization-time of flight mass spectrometry were analyzed using Typer v.4.0 software (Agena Bioscience, San Diego, CA, USA).

Statistical analysis

PLINK software was used for statistical processing and analysis of the correlation differences between the mutant genotypes in the case and control groups. Case: severe HFMD group; Control: general HFMD group. A1: mutant; A2: wild-type (the default is the variant with the lowest allele frequency). A1 frequency is the MAF value. According to the Hardy–Weinberg equilibrium, the selected samples were from a random population. Fisher test was used to compare the genotype frequency between the case group and the control group. P-value represents the statistical difference between both groups, P < 0.05 indicates that there is a significant difference in A1 between the case and control groups; OR < 1 indicates that A1 is protective; OR = 1 indicates that A1 has no relationship with disease; OR > 1 indicates that A1 has a pathogenic effect. The differences of alleles and genotypes were compared by Logistic regression analysis. 95% confidence interval = L95 – U95. L95 and U95 represent the lower and upper limits of the confidence interval, respectively. OR > 1 and L95 > 1 indicate that the allele has a pathogenic effect, while OR < 1 and U95 < 1 indicate that the allele has a protective effect.



The characteristics of the subjects are shown in Table 2. There is no statistically significant difference between the case and control groups in terms of age and sex (P > 0.05).

Table 2 Characteristics of the subjects

Optimized SNP sites

The selected 28 SNP sites in SCARB2 were: rs1051326, rs11547135, rs121909118, rs121909119, rs1465922, rs1465923, rs1470194, rs200053119, rs2119733, rs2869851, rs3733255, rs3733256, rs57374265, rs6811781, rs6812193, rs6824953, rs6841815, rs727502772, rs727502781, rs72857048, rs75285019, rs755903502, rs7697073, rs78737354, rs8475, rs886041074, rs886041076, and rs886041078.

Among the 28 optimized sites, the MAFs for rs6811781, rs6841815, and rs72857048, were not in agreement with those in the 1000 CHB dataset, so it was discarded. rs11547135, rs1465922, and rs2869851 could not be detected so that they were excluded from further analysis. Of the 28 SNPs examined, consideration was given to the accuracy of the reaction system and primers. Therefore, SNPs with low detection rates were not used in further analysis. Therefore, a total of 22 SNPs were analyzed further.

The MAFs of the selected SNPs were greater than 0.01, and the P-values of the Hardy–Weinberg equilibrium test were greater than 0.001. Nine SNPs (rs121909118, rs121909119, rs200053119, rs727502772, rs727502781, rs755903502, rs886041074, rs886041076, and rs886041078) did not fulfill these criteria and were excluded from further analysis. At last, a total of 13 sites were subsequently analysed (Table 3).

Table 3 Genotyping results

Fisher analysis

As shown in Table 4, the frequencies of the rs6812193 T allele was 0.134 and 0.034 in the case and control group, respectively. P value 0.023 < 0.05, indicating a significant difference of A1 between the case and control groups; the OR of 4.381 > 1 indicates that A1 has a pathogenic effect. The remaining 12 SNPs may not be related to the pathogenicity of EV-71. Therefore, the rs6812193 T genotype is a susceptibility SNP.

Table 4 Fisher analysis

Allele logistic regression analysis

As shown in Table 5, the P value of the rs6812193 T allele was 0.0245 < 0.05, indicating a significant difference between the case and control groups; the OR of 4.412 > 1 and L95 value 1.210 > 1 indicate that the allele had a pathogenic effect. The remaining 12 SNPs may not be related to the pathogenicity of EV-71. Therefore, the rs6812193 T genotype is a susceptibility SNP.

Table 5 Allele logistic regression analysis

rs6812193 genotype logistic regression analysis

As shown in Table 6, in the dominant model, the rs6812193 T allele was associated with a risk of severe disease. CT + TT genotype carriers had an increased risk of severe disease compared with CC genotype carriers (OR = 4.56, 95% confidence interval = 1.22–17.04, P = 0.012).

Table 6 rs6812193 genotype logistic regression association with response status (n = 100)


The human SCARB2 gene is located on chromosome 4 and encodes a peptide chain containing 478 amino acids. SCARB2 is a transmembrane sialic acid glycoprotein with a relative molecular mass of 85 kDa, and belongs to the family of CD36 molecules [20]. SCARB2 is mainly located in lysosomes and endosome, and widely present on the membrane of most human cells including nerve cells [21]. This protein is also called lysosomal integral membrane protein 2. It is a type of specific glucose cerebral fat enzyme combined with ligands, involved in the lysosomal pathway. The related research fields are mostly Parkinson's disease with abnormal lysosomal metabolism [22, 23], Gaucher’s disease and myoclonic epilepsy [24, 25]. Yamayoshi et al. [12, 13] found that the tissue distribution of EV-71 virus antigen was well correlated with SCARB2, and further found that this receptor was involved in the endocytosis and membrane transport of pathogenic bacteria. This study speculated that the expression level of SCARB2 might be related to virus sensitivity and infection rate. Therefore, the 28 selected sites were all functional sites related to the expression level, including exons, promoters and introns.

Choi M et al. [26] found that exons contain the vast majority of protein coding synthesis, and about 85% of pathogenic mutations are located in the exon region. In 2009, Ng SB et al. [27] used exome sequencing for the first time to find point mutations located in MYH3 in 4 patients with Freeman Sheldon syndrome (autosomal dominant genetic disease), showing the powerful effect of exome sequencing in identifying pathogenic genes of Mendelian genetic disease. Many complex diseases have been identified by exome sequencing, such as genetic disease OHDO syndrome (KAT6B) [28], CTNNB1 mutation in craniopharyngioma patients [29], point mutation of dilated cardiomyopathy GATAD1 [30], etc. Jenny Do et al. [25] found that 3'-UTR mutations in SCARB2 may be associated with Gaucher disease and myoclonic epilepsy. Yock-Ping Chow et al. [31] found that SCARB2 exon mutation was associated with Pendred syndrome. Yamayoshi et al. [13] found that amino acids at position 142-204 of SCARB2 played an important role in promoting the binding of virus particles to cells and susceptibility to EV-71. However, the study of Ting-Yu Yen et al. did not find the correlation between amino acids at position 142-204 and clinical severity [32]. There are 12 exon sites selected in this study: rs1051326, rs3733255, rs3733256, rs6811781, rs6841815, and rs8475 belong to 3' UTR region, and its function was predicted as miRNA binding site. rs11547135 and rs1465922 belong to 5' UTR region, and its function was predicted to be a TFBS transcription factor binding region. rs121909118 rs200053119 rs755903502 and rs886041078 are NCBI pathogenic clinical significance sites. Among these exons, rs6811781 and rs6841815 (not in CHB); rs11547135 and rs1465922 (not detected); rs121909118, rs200053119, rs755903502, rs886041078 (not in line with MAF value). Finally, rs8475, rs1051326, rs3733255 and rs3733256 were included in the study, but no correlation was found between these four exon loci and the severity of clinical infection.

Promoter is an important cis-element in gene expression regulation and the core region of gene transcriptional regulation. In this study, two promoter loci were selected: rs1465923 and rs78737354. Finally, no correlation between these two promoters and clinical infection was found.

Previously, it was often believed that introns do not encode proteins and do not have biological functions in organisms. However, studies have found that the expression profiles of the same gene with and without introns are significantly different [33]. In many cases of transgenic expression, the addition of a universal intron to cDNA results in a significant increase in gene expression [34, 35]. The optimal expression of many endogenous genes has been demonstrated in mammalian tissue culture cells, transgenic mice, insects, and plant systems requiring the presence of one or more introns. Therefore, a variety of introns in organisms are an important part of eukaryotic genome and are closely related to the construction and dynamic changes of cytoskeleton of gene expression [36]. Ting-Yu Yen studied the relationship between SCARB2, PSGL-1, ANXA2 polymorphisms and clinical severity, and found that rs11097262 was associated with rs6824953 located in the intron region of SCARB2 gene, considering that it may regulate the function or expression of SCARB2 and thus affect the susceptibility to EV-71 [32]. In this study, 14 introns were selected: rs121909119, rs727502772, rs727502781, rs886041074, and rs886041076 were considered to be the sites with pathological clinical significance on NCBI website; rs6812193 is a hot spot site that can be found in the literatures [23, 37,38,39]; rs1470194, rs2119733, rs2869851, rs57374265, rs6824953, rs72857048, rs75285019, and rs7697073 are the TAGSNP sites. Among these introns of this study, rs72857048 is not in CHB; rs2869851 is not detected; rs121909119, rs727502772, rs727502781, rs886041074, and rs886041076 are not in line with MAF value. Finally, rs1470194, rs2119733, rs57374265, rs6812193, rs6824953, rs75285019 and rs7697073 were included in the final study, but only rs6812193 was correlated with the severity of clinical infection, and no correlation was found for the other 6 introns. By using Fisher analysis and allele logistic regression analysis, the rs6812193 T allele was shown to have a pathogenic effect. rs6812193 genotype logistic regression analysis in a dominant model showed that CT + TT genotype carriers had an increased risk of severe HFMD compared with CC genotype carriers. Therefore, the rs6812193 T genotype is considered to be a susceptibility SNP, and the rs6812193 polymorphism may be related to susceptibility to EV-71.

rs6812193 is actually a hot SNP close to SCARB2. A 2011 Web-based Genome-wide Association (GWA) study found that a nucleotide polymorphism rs6812193 close to SCARB2 was significantly associated with Parkinson’s disease (PD) in people of European ancestry [37]. In 2012, Shuai Chen et al. conducted a genotyping study on rs6812193 in 449 PD patients and 452 control patients in mainland China, and found that there is no statistically significant differences in allele and genotype distribution between the patients and the control group [38]. In 2013, Kallirhoe Kalinderi et al. studied 210 Greek patients with sporadic PD and 133 control subjects in Greece. It was found that there was no difference in genotype or allele frequency between PD patients and controls [39]. In 2021, T.S. Usenko et al. rs6812193 of the SCARB2 gene does not confer a significant risk for PD in Russian population [23]. At present, to the best of our knowledge, there is no research related to the role of rs6812193 in HFMD. How does rs6812193 affect the expression and function of SCARB2? This will be examined in the future by increasing the sample size to verify the association of rs6812193 with susceptibility to EV-71 HFMD, and to study the function of rs6812193.

The present study also has some limitations, such as the small number of cases, which may make it difficult to find significant differences between low-frequency SNPs. In addition, with sufficient funds and time, genome-wide tests can be performed to avoid screening for missing sites of interest. Clarifying the relationship between gene polymorphisms and disease will enable us to analyze disease pathogenesis further, to explore the nature of the diversity of disease phenotypes, and to develop more individualized treatment measures.


The rs6812193 T genotype was identified as a susceptibility SNP. CT + TT genotype carriers have an increased risk of severe HFMD compared with CC genotype carriers. Therefore, the rs6812193 polymorphism might be considerably related to clinical severity of enterovirus (EV)-71 associated HFMD, which can support doctors to make evidence-based health recommendations to patients.

Availability of data and materials

The data used and/or analyzed during this study are available from the corresponding author on request.



scavenger receptor class B member 2




hand-foot-mouth disease


severe hand-foot-mouth disease


single nucleotide polymorphism


minor allele frequency


Han Chinese in Beijing


ribonucleic acid


deoxyribonucleic acid




National Center for Biotechnology Information


untranslated region


polymerase chain reaction


polymerase chain reaction primers


high performance liquid chromatography


deoxy-ribonucleoside triphosphate






case group A1 allele frequency


control group A1 allele frequency


odds ratio


lower limits of the confidence interval


upper limits of the confidence interval


coefficient t-statistic


myosin heavy chain 3


lysine acetyltransferase 6B


GATA zinc finger domain containing 1


transcription factor binding site


P-selectin glycoprotein ligand 1


annexin II


Genome-wide Association


Parkinson’s disease


  1. 1.

    Solomon T, Lewthwaite P, Perera D, Cardosa MJ, McMinn P, Ooi MH. Virology, epidemiology, pathogenesis, and control of enterovirus71. Lancet Infect Dis. 2010;10:778–90.

    Article  Google Scholar 

  2. 2.

    Ooi MH, Wong SC, Lewthwaite P, Cardosa MJ, Solomon T. Clinical features, diagnosis, and management of enterovirus 71. Lancet Neurol. 2010;9:1097–105.

    Article  Google Scholar 

  3. 3.

    Li XW, Ni X, Qian SY, Wang Q, Jiang RM, Xu WB, Zhang YC, Yu GJ, Chen Q, Shang YX, Zhao CS, Yu H, Zhang T, Liu G, Deng HL, Gao J, Ran XG, Yang QZ, Xu BL, Huang XY, Wu XD, Bao YX, Chen YP, Chen ZH, Liu QQ, Lu GP, Liu CF, Wang RB, Zhang GL, Gu F, Xu HM, Li Y, Yang T. Chinese guidelines for the diagnosis and treatment of hand, foot and mouth disease, 2018. World J Pediatr. 2018;14:437–47.

    Article  Google Scholar 

  4. 4.

    Liu SL, Pan H, Liu P, Amer S, Chan TC, Zhan J, Huo XX, Liu YZ, Teng Z, Wang L, Zhuang H. Comparative epidemiology and virology of fatal and nonfatal cases of hand, foot and mouth disease in mainland China from 2008 to 2014. Rev Med Virol. 2015;25:115–28.

    Article  Google Scholar 

  5. 5.

    Sanden SVD, Koopmans M, Uslu GK, Avoort HVD. Epidemiology of enterovirus 71 in the Netherlands, 1963 to 2008. J Clin Microbiol. 2009;47:2826–33.

    Article  Google Scholar 

  6. 6.

    Schmidt NJ, Lennette EH, Ho HH. An apparently new enterovirus isolated from patients with disease of the central nervous system. J Infect Dis. 1974;129:304–9.

    CAS  Article  Google Scholar 

  7. 7.

    Chan LG, Parashar UD, Lye MS, Ong FGL, Zaki SR, Alexander JP, Ho KK, Han LL, Pallansch MA, Suleiman AB, Jegathesan M, Anderson LJ. Deaths of children during an outbreak of hand, foot, and mouth disease in Sarawak, Malaysia: clinical and pathological characteristics of the disease. Clin Infect Dis. 2000;31:678–83.

    CAS  Article  Google Scholar 

  8. 8.

    Ho M, Chen ER, Hsu KH, Twu SJ, Chen KT, Tsai SF, Wang JR, Shih SR. An epidemic of enterovirus 71 infection in Taiwan. New Engl J Med. 1999;341:929–35.

    CAS  Article  Google Scholar 

  9. 9.

    Chan KP, Goh KT, Chong CY, Teo ES, Lau G, Ling AE. Epidemic hand, foot and mouth disease caused by human enterovirus 71, Singapore. Emerg Infect Dis. 2003;9:78–85.

    Article  Google Scholar 

  10. 10.

    Zhang Y, Tan XJ, Wang HY, Yan DM, Zhu SL, Wang DY, Ji F, Wang XJ, Gao YJ, Chen L, An HQ, Li DX, Wang SW, Xu AQ, Wang ZJ, Xu WB. An outbreak of hand, foot, and mouth disease associated with subgenotype C4 of human enterovirus 71 in Shandong. China J Clin Virol. 2009;44:262–7.

    Article  Google Scholar 

  11. 11.

    Zhang Y, Zhu Z, Yang WZ, Ren J, Tan XJ, Wang Y, Mao NY, Xu ST, Zhu SL, Cui AL, Zhang Y, Yan DM, Li Q, Dong XP, Zhang J, Zhao YP, Wan JF, Feng ZJ, Sun JL, Wang SW, Li DX, Xu WB. An emerging recombinant human enterovirus 71 responsible for the 2008 outbreak of hand foot and mouth disease in Fuyang city of China. Virol J. 2010;7:94.

    Article  Google Scholar 

  12. 12.

    Yamayoshi S, Iizuka S, Yamashita T, Minagawa H, Mizuta K, Okamoto M, Nishimura H, Sanjoh K, Katsushima N, Itagaki T, Nagai Y, Fujii K, Koike S. Human SCARB2-dependent infection by coxsackievirus A7, A14, and A16 and enterovirus71. J Virol. 2012;86:5686–96.

    CAS  Article  Google Scholar 

  13. 13.

    Yamayoshi S, Ohka S, Fujii K, Koike S. Functional comparison of SCARB2 and PSGL1 as receptors for enterovirus 71. J Virol. 2013;87:3335–47.

    CAS  Article  Google Scholar 

  14. 14.

    Lv TG, Li J, Han ZL, Chen ZB. Association of interleukin-17F gene polymorphism with enterovirus 71 encephalitis in patients with hand, foot, and mouth disease. Inflammation. 2013;36:977–81.

    CAS  Article  Google Scholar 

  15. 15.

    Yuan AY, Li J, Liu PP, Chen ZB, Hou M, Wang JJ, Han ZL. Association of interleukin -6 -572C/G gene polymorphism and serum or cerebrospinal fluid interleukin-6 level with enterovirus 71 encephalitis in Chinese Han patients with hand, foot, and mouth disease. Inflammation. 2015;38:728–35.

    CAS  Article  Google Scholar 

  16. 16.

    Zhao N, Chen HL, Chen ZZ, Li J, Chen ZB. IL-10-592 polymorphism is associated with IL-10 expression and severity of enterovirus 71 infection in chinese children. J Clin Virol. 2017;95:42–6.

    CAS  Article  Google Scholar 

  17. 17.

    Li JA, Chen ZB, Lv TG, Han ZL. Genetic polymorphism of CCL2-2518, CXCL10-201, IL8+781 and susceptibility to severity of enterovirus-71 infection in a Chinese population. Inflamm Res. 2014;63:549–56.

    CAS  Article  Google Scholar 

  18. 18.

    Liu YD, Liu PP, Liu SH, Guo Y, He HF, Yang CQ, Song J, Zhang N, Cheng JG, Chen ZB. Oligoadenylate synthetase 3 S381R gene polymorphism is associated with severity of EV71 infection in Chinese children. J Clin Virol. 2018;101:29–33.

    CAS  Article  Google Scholar 

  19. 19.

    Zhang XA, Xu HM, Chen XD, Li XJ, Wang XJ, Ding SJ, Zhang RL, Liu LJ, He C, Zhuang L, Li H, Zhang PH, Yang H, Li TY, Liu W, Cao WC. Association of functional polymorphisms in the MxA gene with susceptibility to enterovirus 71 infection. Hum Genet. 2014;133:187–97.

    CAS  Article  Google Scholar 

  20. 20.

    Eskelinen EL, Tanaka Y, Saftig P. At the acidic edge: emerging functions for lysosomal membrane proteins. Trends Cell Biol. 2003;13:137–45.

    CAS  Article  Google Scholar 

  21. 21.

    Gonzalez A, Valeiras M, Sidransky E, Tayebi N. Lysosomal integral membrane protein-2: a new player in lysosome-related pathology. Mol Genet Metab. 2014;111:84–91.

    CAS  Article  Google Scholar 

  22. 22.

    Michelakakis H, Xiromerisiou G, Dardiotis E, Bozi M, Vassilatis D, Kountra PM, Patramani G, Moraitou M, Papadimitriou D, Stamboulis E, Stefanis L, Zintzaras E, Hadjigeorgiou DP. Evidence of an association between the scavenger receptor class B member 2 gene and Parkinson’s disease. Mov Disord. 2012;27:400–5.

    CAS  Article  Google Scholar 

  23. 23.

    Usenko TS, Bezrukova AI, Bogdanova DA, Kopytova AE, Senkevich KA, Gracheva EV, Timofeeva AA, Miliukhina IV, Zakharova EY, Emelyanov AK, Pchelina SN. Genetics variants and expression of the SCARB2 gene in the pathogenesis of Parkinson’s disease in Russia. Neurosci Lett. 2021;741:135509.

    CAS  Article  Google Scholar 

  24. 24.

    Mazzulli JR, Xu YH, Sun Y, Knight AL, McLean PJ, Caldwell GA, Sidransky E, Grabowski GA, Krainc D. Gaucher disease glucocerebrosidase and α-synuclein form a bidirectional pathogenic loop in synucleinopathies. Cell. 2011;146:37–52.

    CAS  Article  Google Scholar 

  25. 25.

    Do J, Gary S, Stubblefifield B, Ryan E, Lopez G, Sidransky E, Tayebi N. A 3’-UTR variant in SCARB2 modulates LIMP2 in patients with Gaucher disease and myoclonic epilepsy. Abstr Mol Genet and Metab. 2019;126:S49.

    Article  Google Scholar 

  26. 26.

    Choi M, Scholl UI, Ji WZ, Liu TW, Tikhonova IR, Zumbo P, Nayir A, Bakkaloğlu A, Özen S, Sanjad S, Nelson-Williams C, Farhi A, Mane S, Lifton RP. Genetic diagnosis by whole exome capture and massively parallel DNA sequencing. Proc Natl Acad Sci USA. 2009;106:19096–101.

    CAS  Article  Google Scholar 

  27. 27.

    Ng SB, Turner EH, Robertson PD, Flygare SD, Bigham AW, Lee C, Shaffer T, Wong M, Bhattacharjee A, Eichler EE, Bamshad M, Nickerson DA, Shendure J. Targeted capture and massively parallel sequencing of twelve human exomes. Nature. 2009;461:272–6.

    CAS  Article  Google Scholar 

  28. 28.

    Clayton-Smith J, O’Sullivan J, Daly S, Bhaskar S, Day R, Anderson B, Voss AK, Thomas T, Biesecker LG, Smith P, Fryer A, Chandler KE, Kerr B, Tassabehji M, Lynch SA, Krajewska-Walasek M, McKee S, Smith J, Sweeney E, Mansour S, Mohammed S, Donnai D, Black G. Whole-exome-sequencing identifies mutations in histone acetyltransferase gene KAT6B in individuals with the Say-Barber-Biesecker variant of Ohdo syndrome. Am J Hum Genet. 2011;89:675–81.

    CAS  Article  Google Scholar 

  29. 29.

    Brastianos PK, Taylor-Weiner A, Manley PE, Jones RT, Dias-Santagata D, Thorner AR, Rodriguez FJ, Bernardo LA, Schubert L, Sunkavalli A, Shillingford N, Calicchio ML, Lidov HGW, Taha H, Martinez-Lage M, Santi M, Storm PB, Lee JYK, Palmer JN, Adappa ND, Scott RM, Dunn IF, Laws ER, Stewart C, Ligon KL, Hoang MP, Hummelen PV, Hahn WC, Louis DN, Resnick AC, Kieran MW, Getz G, Santagata S. Exome sequencing identifies BRAF mutations in papillary craniopharyngiomas. Nat Genet. 2014;46:161–5.

    CAS  Article  Google Scholar 

  30. 30.

    Theis JL, Sharpe KM, Matsumoto ME, Chai HS, Nair AA, Theis JD, Andrade MD, Wieben ED, Michels VV, Olson TM. Homozygosity mapping and exome sequencing reveal GATAD1 mutation in autosomal recessive dilated cardiomyopathy. Circ Cardiovasc Genet. 2011;4:585–94.

    CAS  Article  Google Scholar 

  31. 31.

    Chow YP, Abdul Murad NA, Rani ZM, Khoo JS, Chong PS, Wu LL, Jamal R. Exome sequencing identifies SLC26A4, GJB2, SCARB2 and DUOX2 mutations in 2siblings with Pendred syndrome in a Malaysian family. Orphanet J Rare Dis. 2017;12:40.

    Article  Google Scholar 

  32. 32.

    Yen TY, Shih WL, Huang YC, Lee JT, Huang LM, Chang LY. Polymorphisms in enterovirus 71 receptors associated with susceptibility and clinical Severity. PLoS ONE. 2018;13:e0206769.

    Article  Google Scholar 

  33. 33.

    Nott A, Meislin SH, Moore MJ. A quantitative analysis of intron effects on mammalian gene expression. RNA. 2003;9(5):607–17.

    CAS  Article  Google Scholar 

  34. 34.

    Palmiter RD, Sandgren EP, Avarbock MR, Allen DD. Heterologous introns can enhance expression of transgenes in mice. Proc Natl Acad Sci USA. 1991;88:478–82.

    CAS  Article  Google Scholar 

  35. 35.

    Bourdon V, Harvey A, Lonsdale DM. Introns and their positions affect the translational activity of mRNA in plant cells. EMBO Rep. 2001;2:394–8.

    CAS  Article  Google Scholar 

  36. 36.

    Jeong YM, Mun JH, Lee I, Woo JC, Hong CB, Kim SG. Distinct roles of the first introns on the expression of arabidopsis profilin gene family members. Plant Physiol. 2006;140:196–209.

    CAS  Article  Google Scholar 

  37. 37.

    Do CB, Tung JY, Dorfman E, Kiefer AK, Drabant EM, Francke U, Mountain JL, Goldman SM, Tanner CM, Langston JW, Wojcicki A, Eriksson N. Web-based genome-wide association study identifies two novel loci and a substantial genetic component for Parkinson’s disease. PLoS Genet. 2011;7:e1002141.

    CAS  Article  Google Scholar 

  38. 38.

    Chen S, Zhang Y, Chen W, Wang Y, Liu J, Rong TY, Ma JF, Wang G, Zhang J, Pan J, Xiao Q, Chen SD. Association study of SCARB2 rs6812193 polymorphism with Parkinson’s disease in Han Chinese. Neurosci Lett. 2012;516:21–3.

    CAS  Article  Google Scholar 

  39. 39.

    Kalinderi K, Bostantjopoulou S, Katsarou Z, Fidani L. Association study of rs6812193 polymorphism with Parkinson’s disease in a Greek population. Neurosci Lett. 2013;541:190–2.

    CAS  Article  Google Scholar 

Download references


The authors gratefully acknowledge the editors and anonymous reviewers for their valuable comments on this manuscript.


This research was funded by the Key Project of Tianjin Second People’s Hospital (No. YS0019), and the 13th Five-year National Major Project for optimization of prevention and treatment plan for febrile rash syndrome based on HFMD (No. KY0112 2017ZX10103007-002).

Author information




Conceptualization and methodology, YL and RS; software and validation, XW and HL; formal analysis, XW; investigation, XW, HL, YL, RS, YL, KQ; resources, YL and RS; writing—original draft preparation, XW; writing—review and editing, XW, HL, YL, RS, YL, KQ; supervision, YL and RS; project administration, YL; funding acquisition, YL and XW. All authors read and approved the final manuscript.

Corresponding authors

Correspondence to Ying Li or Rui Su.

Ethics declarations

Ethics approval and consent to participate

The study was conducted according to the guidelines of the Declaration of Helsinki, and approved by the Institutional Review Board of Medical Ethics Committee of Tianjin Second People’s Hospital (Approval Certificate of Ethical Review No. 201825).

Consent for publication

No consent for publication applicable. This manuscript does not provide any patients’ data nor any animal studies or experiments. This manuscript does not contain any individual person’s data in any form. The authors declare no financial and non-financial conflict of interests. This manuscript has not been submitted to, nor is under review at, another journal or other publishing venue. The authors have no affiliation with any organization with a direct or indirect financial interest in the subject matter discussed in the manuscript.

Competing interests

The authors declare that they have no competing interests.

Additional information

Publisher's Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Rights and permissions

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit The Creative Commons Public Domain Dedication waiver ( applies to the data made available in this article, unless otherwise stated in a credit line to the data.

Reprints and Permissions

About this article

Verify currency and authenticity via CrossMark

Cite this article

Wang, X., Liu, H., Li, Y. et al. Relationship between polymorphism of receptor SCARB2 gene and clinical severity of enterovirus-71 associated hand-foot-mouth disease. Virol J 18, 132 (2021).

Download citation


  • Enterovirus 71
  • Human scavenger receptor B2
  • Single nucleotide polymorphism
  • Severe hand-foot-mouth disease