Oligoadenylate synthetase 3 S381R gene polymorphism is associated with severity of EV71 infection in Chinese children

EV71 3C protease cleaves host anti-viral factor OAS3 and enhances virus replication

The global spread of enteroviruses (EVs) has become more frequent, severe and life-threatening. Intereron (IFN) I has been proved to control EVs by regulating IFN-stimulated genes (ISG) expression. 2'-5'-oligoadenylate synthetases 3 (OAS3) is an important ISG in the OAS/RNase L antiviral system. The relationship between OAS3 and EVs is still unclear. Here, we reveal that OAS3, superior to OAS1 and OAS2, significantly inhibited EV71 replication in vitro. However, EV71 utilized autologous 3C protease (3Cpro) to cleave intracellular OAS3 and enhance viral replication. Rupintrivir, a human rhinovirus 3C protease inhibitor, completely abolished the cleavage of EV71 3Cpro on OAS3. And the proteolytically deficient mutants H40G, E71A, and C147G of EV71 3Cpro also lost the ability of OAS3 cleavage. Mechanistically, the Q982-G983 motif in C-terminal of OAS3 was identified as a crucial 3Cpro cutting site. Further investigation indicated that OAS3 inhibited not only EV71 but also Coxsackievirus B3 (CVB3), Coxsackievirus A16 (CA16), Enterovirus D68 (EVD68), and Coxsackievirus A6 (CA6) subtypes. Notably, unlike other four subtypes, CA16 3Cpro could not cleave OAS3. Two key amino acids variation Ile36 and Val86 in CA16 3Cpro might result in weak and delayed virus replication of CA16 because of failure of OAS and 3AB cleavage. Our works elucidate the broad anti-EVs function of OAS3, and illuminate a novel mechanism by which EV71 use 3Cpro to escape the antiviral effect of OAS3. These findings can be an important entry point for developing novel therapeutic strategies for multiple EVs infection.

Association of the IRAK4 rs4251545 genetic polymorphism with severity of enterovirus-71 infection in Chinese children

INTRODUCTION

This study aimed to explore the association between the IRAK4 polymorphism rs4251545 and the severity of enterovirus 71 (EV71) infection in Chinese children.

METHODS

We analyzed the IRAK4 polymorphism rs4251545 in 617 EV71-infected patients and 410 controls using the improved multiplex ligation detection reaction. IRAK4 mRNA expression was tested by qRT-PCR. Serum concentrations of IL-6 and NF-κB were detected using ELISA.

RESULTS

The frequencies of the GA + AA genotype and A allele in the mild EV71 infection group and in the severe EV71 infection group were significantly higher than those in the normal control group. The frequency of the GA + AA genotype and A allele in severely infected EV71 patients was markedly higher than that in mildly infected EV71 patients. IRAK4 mRNA expression in mildly infected EV71 patients and severely infected patients was significantly higher than that in the control group. IRAK4 mRNA expression in GA + AA genotypes in both mild and severe EV71 infection groups was significantly higher than that in patients with the GG genotype. IL-6 concentration and the ratio of IL-6/NF-κB in severe EV71 cases were significantly lower in patients with the GA + AA genotype than in those with the GG genotype. The ratio of IL-6/NF-κB was distinctly higher in severely infected EV71 patients than in mildly infected and control subjects.

CONCLUSIONS

The IRAK4 polymorphism rs4251545 was associated with the susceptibility and severity of EV71 infection. The A allele is a susceptible factor in the development of severe EV71 infection in Chinese children.

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

BACKGROUND

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).

METHODS

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.

RESULTS

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).

CONCLUSION

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.

Study on the differentially expressed genes and signaling pathways in dermatomyositis using integrated bioinformatics method

Dermatomyositis is a common connective tissue disease. The occurrence and development of dermatomyositis is a result of multiple factors, but its exact pathogenesis has not been fully elucidated. Here, we used biological information method to explore and predict the major disease related genes of dermatomyositis and to find the underlying pathogenic molecular mechanism.The gene expression data of GDS1956, GDS2153, GDS2855, and GDS3417 including 94 specimens, 66 cases of dermatomyositis specimens and 28 cases of normal specimens, were obtained from the Gene Expression Omnibus database. The 4 microarray gene data groups were combined to get differentially expressed genes (DEGs). The gene ontology and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway enrichments of DEGs were operated by the database for annotation, visualization and integrated discovery and KEGG orthology based annotation system databases, separately. The protein-protein interaction networks of the DEGs were built from the STRING website. A total of 4097 DEGs were extracted from the 4 Gene Expression Omnibus datasets, of which 2213 genes were upregulated, and 1884 genes were downregulated. Gene ontology analysis indicated that the biological functions of DEGs focused primarily on response to virus, type I interferon signaling pathway and negative regulation of viral genome replication. The main cellular components include extracellular space, cytoplasm, and blood microparticle. The molecular functions include protein binding, double-stranded RNA binding and MHC class I protein binding. KEGG pathway analysis showed that these DEGs were mainly involved in the toll-like receptor signaling pathway, cytosolic DNA-sensing pathway, RIG-I-like receptor signaling pathway, complement and coagulation cascades, arginine and proline metabolism, phagosome signaling pathway. The following 13 closely related genes, XAF1, NT5E, UGCG, GBP2, TLR3, DDX58, STAT1, GBP1, PLSCR1, OAS3, SP100, IGK, and RSAD2, were key nodes from the protein-protein interaction network.This research suggests that exploring for DEGs and pathways in dermatomyositis using integrated bioinformatics methods could help us realize the molecular mechanism underlying the development of dermatomyositis, be of actual implication for the early detection and prophylaxis of dermatomyositis and afford reliable goals for the curing of dermatomyositis.

Molecular epidemiology of enteroviruses associated with severe hand, foot and mouth disease in Shenzhen, China, 2014-2018

In this study, we investigated the epidemiology and molecular characteristics of enteroviruses associated with severe hand, foot and mouth disease (HFMD) in Shenzhen, China, during 2014-2018. A total of 137 fecal specimens from patients with severe HFMD were collected. Enterovirus (EV) types were determined using real-time reverse transcription polymerase chain reaction (RT-PCR), RT nested PCR, and sequencing. Sequences were analyzed using bioinformatics programs. Of 137 specimens tested, 97 (70.8%), 12 (8.8%), and 10 (7.3%) were positive for EV-A71, coxsackievirus A6 (CVA6), and CVA16, respectively. Other pathogens detected included CVA2 (2.9%, 4/137), CVA10 (2.9%, 4/137), CVA5 (0.7%, 1/137), echovirus 6 (E6) (0.7%, 1/137) and E18 (0.7%, 1/137). The most frequent complication in patients with proven EV infections was myoclonic jerk, followed by aseptic encephalitis, tachypnea, and vomiting. The frequencies of vomiting and abnormal eye movements were higher in EV-A71-infected patients than that in CVA6-infected or CVA16-infected patients. Molecular phylogeny based on the complete VP1 gene revealed no association between the subgenotype of the virus and disease severity. Nevertheless, 12 significant mutations that were likely to be associated with virulence or the clinical phenotype were observed in the 5’UTR, 2A^pro, 2C, 3A, 3D^pol and 3’UTR of CVA6. Eight significant mutations were observed in the 5’UTR, 2B, 3A, 3D^pol and 3’UTR of CVA16, and 10 significant mutations were observed in the 5’UTR, VP1, 3A and 3C^pro of CVA10. In conclusion, EV-A71 is still the main pathogen causing severe HFMD, although other EV types can also cause severe complications. Potential virulence or phenotype-associated sites were identified in the genomes of CVA6, CVA16, and CVA10.

Genome-wide association study identifies TPH2 variant as a novel locus for severe CV-A6-associated hand, foot, and mouth disease in Han Chinese

BACKGROUND

Little is known about the association between genetic susceptibility and the severity of hand, foot, and mouth disease (HFMD) infected with coxsackievirus A6 (CV-A6).

METHODS

Three hundred and sixty-four CV-A6 HFMD patients were enrolled, including 115 severe and 249 mild patients. A genome-wide association study (GWAS) was performed involving eight DNA pools of 115 severe and 115 mild CV-A6 HFMD patients pair-matched by age and gender. Differences in relative allele signal scores of SNPs in Illumina Human OmniZhongHua-8 BeadChips were compared between the two groups. The tag SNPS for potentially functional SNPs or their high linked SNPs were selected for individual genotyping in all 364 patients and assessed for their associations with severe CV-A6 HFMD using multivariable logistic regression analyses.

RESULTS

The top 30 significant SNPs obtained from pooled DNA GWAS analysis were checked for biological functions and their high linkage disequilibrium (LD) SNPs. Four tag SNPs (rs1558206, rs6927647, rs9375728 and rs10879355) were selected for further individual genotyping in 364 CV-A6 patients. Only SNP rs10879355 was associated with severe CV-A6 HFMD, with CC genotype having a greater risk of severe illness than TT+TC genotypes (OR=2.48, 95%CI: 1.34, 4.56). SNP rs4290270 is in complete LD with rs10879355 in Chinese Han children.

CONCLUSIONS

This is the first report that one potentially functional SNP rs4290270 in the TPH2 gene may be associated with the risk of severe CV-A6 HFMD.

Targeted exome sequencing identifies five novel loci at genome-wide significance for modulating antidepressant response in patients with major depressive disorder

In order to determine the role of single nucleotide variants (SNVs) in modulating antidepressant response, we conducted a study, consisting of 929 major depressive disorder (MDD) patients, who were treated with antidepressant drugs (drug-only) or in combination with a repetitive transcranial magnetic stimulation (plus-rTMS), followed by targeted exome sequencing analysis. We found that the "plus-rTMS" patients presented a more effective response to the treatment when compared to the 'drug-only' group. Our data firstly demonstrated that the SNV burden had a significant impact on the antidepressant response presented in the "drug-only" group, but was limited in the "plus-rTMS" group. Further, after controlling for overall SNV burden, seven single nucleotide polymorphisms (SNPs) at five loci, IL1A, GNA15, PPP2CB, PLA2G4C, and GBA, were identified as affecting the antidepressant response at genome-wide significance (P < 5 × 10^-08). Additional multiple variants achieved a level of correction for multiple testing, including GNA11, also shown as a strong signal for MDD risk. Our study showed some promising evidence on genetic variants that could be used as individualized therapeutic guides for MDD patients.