- Open Access
T. gondii rhoptry protein ROP18 induces apoptosis of neural cells via endoplasmic reticulum stress pathway
Parasites & Vectorsvolume 8, Article number: 554 (2015)
The neurotropic parasite T. gondii is widespread among mammalian hosts including humans. During the course of T. gondii infection, the central nervous system is the most commonly damaged of all invasive organs. The polymorphic rhoptry protein ROP18 has been identified as a key factor in the pathogenesis of T. gondii; however, the molecular mechanism by which this protein exerts neuropathogenesis remains elusive.
Immunofluorescence staining was performed to detect neuropathogenesis of the mouse brain tissues. The apoptosis of neural cells and the expressions of related proteins in the endoplasmic reticulum stress (ER Stress)-mediated apoptosis pathway were detected by flow cytometry and Western blotting.
Immunofluorescence staining reveals induction of the propidium iodide (PI) - positive neural cells in mouse cerebral cortex and hippocampus infected with ROP18 over-expressing transgenic tachyzoites. Western blotting analyses reveal that ROP18 increases the expressions of cleaved caspase-12, CHOP and cleaved caspase-3 when compared to the control groups. After the pretreatment of Z-ATAD-FMK (a specific caspase-12 inhibitor), the apoptotic level of neural cells had an apparent decline, and correspondingly, the expressions of those related proteins were notably decreased.
Our findings here highlight that the virulence factor ROP18 in T. gondii may contribute to neuronal apoptosis through the ER stress-mediated apoptosis pathway, which may be a potential molecular mechanism responsible for neurological disorders of toxoplasmosis.
T. gondii is a ubiquitous obligate intracellular parasite infecting a wide range of mammalian hosts including humans [1–3]. Infections with T. gondii are generally subclinical in healthy individuals, but toxoplasmosis remains a major problem in immunosuppressed adults or developing fetuses . In patients with severe immune dysfunction, a reactivation of the infection produces the neurological manifestations or even fatal Toxoplasmic encephalitis (TE) due to reactivation of Toxoplasma cysts in the brain . For those patients, dormant encysted bradyzoites can reactivate into fast replicating tachyzoites and cause severe damage to the brain . TE is the most vital outcome of toxoplasmosis in immunosuppressed individuals and the main symptoms including focal seizures, cranial nerve disturbances, altered mental state, ataxia, sensory abnormalities, hemiparesis, meningismus and neuropsychiatric disorders as well .
Although toxoplasmosis can be associated with multi-organ involvement, the central nervous system (CNS) is the most commonly affected of all invasive organs [7, 8]. As a neurotropic parasite, little is known about which effectors of the type I RH strain elicit neuropathogenesis during infection [7, 9]. Previous work revealed a polymorphic, parasite rhoptry protein, known as ROP18, is a key virulence determinant among different T. gondii clonal lineages [10–13]. The recent results have shown that ROP18 has several targets in the host cell, including IRGs (immunity-related GTPases) [14, 15] and NF-κB p65 . IRGs are strongly induced by interferon-γ (IFN-γ) and are important in the innate immune response against T. gondii . ROP18 contributes to avoidance of IRG recruitment to the parasitophorous vacuole membrane (PVM), thus protecting the parasite from clearance in interferon-activated macrophages. The nuclear factor NF-κB transcription factor has essential roles in immune and inflammatory responses. Its p65 subunit is regulated by several post-translational modifications, including phosphorylation, acetylation and ubiquitination. Our recent results show that ROP18 phosphorylates the host p65 and targets this protein to the ubiquitin-dependent degradation, thus inhibiting the NF-κB pathway in infected macrophages .
Our previous study also revealed that both the canonical type I RH strain and TgCtwh3, a representative strain prevalent in China, can induce the apoptosis of neural stem cells through ER stress signaling pathways [17, 18]. In this study, we investigated the effect of rhoptry protein ROP18 on the apoptosis of neural cells. Our results presented here showed that ROP18 can stimulate neural cell death through inducing ER stress-mediated apoptosis pathway.
The study protocol was approved from the Institutional Review Board of the Institute of Biomedicine at Anhui Medical University (Permit Number AMU 26–093628), which records and regulates all research activities in the school. All surgeries were performed under anesthesia, and all efforts were made to minimize animal suffering.
Parasite and cell lines
ROP18 over-expressing transgenic RH strain (ROP18-RH) was constructed as previously described. Briefly, the 5′-UTR-TUB promoter-ROP18 -Ty1-HXGPRT-3′ -UTR fragment was transfected into the Δku80Δhxgprt RH strain parasites (kindly provided by Professor John C. Boothroyd, Stanford University, USA) by electroporation. Stable integrants were selected in media with 50 μg/ml mycophenolic acid and 50 μg/ml xanthine and cloned by limiting dilution. Construction of ROP18-RH strain was confirmed by PCR, IF and Western-blotting . Then the parasites of ROP18-RH strain were harvested from the mouse peritoneal exudates by injection on the third day after infection, and were isolated by centrifugation to discard the contaminating host cells. Then the parasites were maintained by serial passage in the mouse N2A cells (ATCC, Neuro2a) for further experiments. The mouse neuroblastoma Neuro2a (N2a) cells were cultured in DMEM supplemented with 20 % FBS and 1 % penicillin/streptomycin in a humidified 5 % CO2 atmosphere. All parasite strains and cell lines were routinely assessed for mycoplasma contamination, and no contamination was detected.
Plasmid and reagents
Amplification of the open reading frame encoding T. gondii ROP18 (GenBank ID: AM075204.1) was achieved through RT-PCR of the RH tachyzoite RNA. Caspase-12, CHOP and caspase-3 antibodies were purchased from Cell Signaling Technology (USA); GAPDH antibody was purchased from Santa Cruz (USA).
Brain sections were hydrated and rinsed in PBS. Then they were placed in ACSF (artificial cerebrospinal fluid) containing 5 μg/ml PI for 30 min. After antigen retrieval, the brain sections were premeabilized, then incubated with primary antibody at 4 °C overnight. FITC-conjugated goat anti-mouse IgG, rhodamine-conjugated goat anti-rabbit IgG and DAPI dye were used for antigen and DNA visualization. The images were captured using a fluorescent microscopy (Olympus BX60, Tokyo, Japan).
Detection of apoptosis
The apoptosis of N2a cells were determined following the instruction of Annexin V-PE/7-AAD kit (BD, USA). Briefly, the N2a cells were harvested, then washed twice with cold PBS, and incubated in 1 × binding buffer (10 mM HEPES, 0.14 M NaCl, and 0.25 mM CaCl2) containing PE Annexin V and 7-AAD for 15 min. Then stained cells were analyzed using a Faces Calibur flow cytometer (BD Biosciences, USA) within 1 h, and the data were analyzed using FCS Express 4.0 software. Annexin V-PE+/7-AAD − cells represent the early apoptotic cells, and annexinV-PE+/7-AAD+ cells represent the late apoptotic cells.
Western blotting analysis
To further identify the apoptosis of N2a cells, the expressions of caspase-12, CHOP and caspase-3 were determined using Western blotting analysis. Western blotting was conducted as described previously . Briefly, after the N2a cells were infected with ROP18-RH or RH tachyzoites at an m.o.i of ~3 for 24 h, they were harvested, washed with PBS, and lysed in a lysis buffer (50 mM HEPES, pH7.4, 150 mM NaCl, 2 mM EGTA, 1 % Triton X-100, 1 mM phenylmethylsulfonyl fluoride, 10 g/ml leupeptin, and 10 g/ml pepstatin A) containing a protease inhibitor mixture (Sigma, USA). The cell lysates were separated through SDS-PAGE, transferred to nitrocellulose membranes, probed with the corresponding antibodies, and developed using an ECL kit.
All quantitative results were expressed as mean ± SD. Statistical analysis of variance used was a two-way ANOVA followed by the Scheffe’s test. P < 0.05 was considered statistically significant.
ROP18 induced the apoptosis of neural cells in vivo
To assay the effects of ROP18 on the survival of neurons during infection, we infected BALB/c mice intraperitoneally with 1000 tachyzoites from ROP18 over-expressing parasites (ROP18-RH) or RH tachyzoites . Pathology in the brain tissues of animals was examined as soon as obvious clinical manifestations were observed. Immunofluorescent staining NeuN (a specific neuron marker) and propidium iodide (PI) staining were used to reveal the neuronal damage. We found that in the un-infected group, the cerebral cortex was intact and there were many neurons there. However, the loss of neurons in the hippocampus of the mice infected with ROP18-RH or RH tachyzoites were 96.82 ± 2.62 % and 63.74 ± 7.38 %, significantly higher than those in the un-infected group (30.24 ± 5.2 %), accordingly (**P < 0.001, ***P < 0.0001). Consistently, the loss of neurons in the cerebral cortex of the mice infected with ROP18-RH or RH tachyzoites were 96.05 ± 3.92 and 77.28 ± 3.67 %, significantly higher than those in the un-infected group (1.88 ± 1.08 %), accordingly (***P < 0.0001). In addition, the numbers of dead neural cells in the cerebral cortex as well as in the hippocampus of the mice infected with ROP18-RH tachyzoites were significantly higher than RH tachyzoites infected group, accordingly (**P < 0.001) (Fig. 1). The results suggest ROP18 may induce the loss of neurons during infection in vivo.
ROP18 induced the apoptosis of neural cells in vitro
To investigate if the observed cell death of neurons from the mice infected with ROP18-RH or RH tachyzoites result from apoptosis, we performed FCM analysis with Annexin V-PE/7-AAD staining in a double-fluorescence assay. The mouse neuroblastoma line N2a cells were transfected with expression vectors for ROP18 or another parasite protein present in the parasitophorous vacuole, 14-3-3. We then ectopically expressed these parasite proteins and analyzed the effects on the apoptosis of N2a cells. Strikingly, the over-expression of ROP18, but not 14-3-3, in N2a cells led to high levels of apoptosis (Fig. 2a). To investigate this further, we infected N2a cells with ROP18-RH or RH tachyzoites. FCM analysis showed that ROP18-RH, compared with RH tachyzoites, obviously enhanced the apoptosis of neural cells (Fig. 2b). Collectively, these results provide that ROP18 induced the apoptosis of neural cells in vitro.
ROP18 induced the apoptosis of neural cells via ER stress pathway
Our previous results have shown that both the RH and TgCtwh3 strains can induce the apoptosis of neural stem cells through ER stress signaling pathways [17, 18]. To examine whether ROP18 was the effective molecule involved in the ER Stress-mediated apoptosis pathway, N2a cells transfected with ROP18-GFP or GFP vector were treated with 5 μmol/L Z-ATAD-FMK for 6 h and FCM analysis was performed. Consistently, the apoptotic rates of cells transfected with ROP18-GFP were significantly higher than the control groups (Fig. 3a). After treatment with Z-ATAD-FMK, the apoptotic rates of the cells transfected with ROP18-GFP were 8.87 ± 0.38 %, significantly lower than those in the Z-ATAD-FMK untreated ROP18-GFP group (15.04 ± 0.91 %), accordingly (*P < 0.05) (Fig. 3a). To further investigate it, the cells infected with ROP18-RH or RH tachyzoites were treated with 5 μmol/L Z-ATAD-FMK for 6 h. After treatment with Z-ATAD-FMK, the apoptotic rates of cells infected with ROP18-RH or RH tachyzoites were 12.98 ± 0.25 % and 8.05 ± 0.15 %, significantly lower than those in the Z-ATAD-FMK untreated ROP18-RH group (18.25 ± 0.62 %) and untreated RH group (11.53 ± 1.13 %), accordingly (*P < 0.05, **P < 0.01) (Fig. 3b). Consistently, the apoptotic rates of the cells infected with ROP18-RH tachyzoites (18.25 ± 0.62 %) were significantly higher than RH tachyzoites infected group (11.53 ± 1.13 %) (Fig. 3b). Immunoblotting analysis further displays that the expression levels of cleaved caspase-12, CHOP and cleaved caspase-3 in the N2a cells increased significantly in both the ROP18-RH and RH tachyzoites infected groups when compared to the control group. When the cells were pretreated with a caspase-12 inhibitor, Z-ATAD-FMK, the expressions of cleaved caspase-12 and caspase-3 in the inhibitor-pretreated groups significantly decreased (Fig. 4) whereas no significant differences were detected in the levels of CHOP when they were compared to the control groups.
Life-threatening toxoplasmic encephalitis is a common opportunistic infection of the CNS that usually results from reactivation of latent cysts into tachyzoites in the brain tissue of immunosuppressed individuals [19, 20]. After reaching the CNS, the parasites can invade all nucleated cells and initiate activation of resident microglia and astrocytes [21–23]. In the brain tissues, microglial cells, astrocytes and neurons are all susceptible to T. gondii infection [2, 24]. Microglia cells are implicated as essential mediators of the brain’s immune response to injury, inflammation or the presence of pathogens [25, 26]. During the course of T. gondii infection, microglia cells can secrete pro-inflammatory cytokines and lead to the neuronal cell damage of inflammatory surroundings . However, neurons are fundamental elements of the CNS, we hypothesized that there may be the possible direct mechanisms of T. gondii on neurons to disrupt their survival and function . Accumulated studies conducted in different animal experimental models have found that behavioural changes occur upon acute and chronic T. gondii infection . One of the most ‘convenient’ explanations for the neuropsychological and behavioural deficits could be that the parasites directly inject effective molecules into the neurons interfering with their activity and function .
As noted, ROP18 is a Ser/Thr kinase related to the ROP2 family, secreted from rhoptry organelles and relocalized to the PVM during invasion [14, 16, 29, 30]. During host cell invasion, ROP18 is secreted to the PVM where it is tethered to the cytosolic face of this host-pathogen interface. They are thought to modify the host cell environment, thus favoring intracellular parasitism . Although the functional role of rhoptry proteins in the biology of T. gondii has not been clearly elucidated, numerous studies have highlighted the role of ROP18 as the key virulence factor in the pathogenesis of T. gondii infection [10, 11, 31]. Indeed, we have previously demonstrated that ROP18 kinase activity is responsible for the phosphorylation and degradation of p65, thereby down-regulating Th1 responses and causing M2 phenotypes in macrophage . Here, we report that ROP18 is involved in the apoptosis of neural cells.
ER stress conditions have been observed in numerous diseases including ischemia/reperfusion injury, neurodegeneration, as well as infectious diseases [32–35]. Conditions that interfere with ER function lead to the accumulation and aggregation of unfolded proteins. ER trans-membrane receptors detect the onset of ER stress and initiate the unfolded protein response (UPR) to restore normal ER function. In eukaryotic cells, three ER trans-membrane proteins mediate the canonical UPR: the two kinases, IRE1 (inositol requiring enzyme 1) and PERK (PKR-like eukaryotic initiation factor 2a kinase), and the transcription factor precursor ATF6 (activating transcription factor 6). If the stress is prolonged, or the adaptive response fails, apoptotic cell death ensues through activation of caspase-12, CHOP and/or JNK (c-JUN NH2-terminal kinase) [36–38]. Many of these changes are induced by intracellular pathogens, so it is not surprising that ER stress is induced by cellular pathogens and viruses alike [32, 33]. Meanwhile, it has reported that T. gondii ROP18 is evolved mechanisms to mediate degradation of the host endoplasmic reticulum-localizing transcription factor, ATF6β, to down-regulate T cell-mediated type I adaptive immune responses .
In the present study, we detected the expressions of related proteins in the ER stress-mediated apoptosis pathway by Immunoblotting analysis, and the results showed that ROP18 enhanced the expressions of cleaved caspase-12, CHOP and cleaved caspase-3 in the neural cells (Fig. 4), and subsequently increased the apoptosis of the neural cells infected by ROP18 over-expressing parasite (Fig. 2). Moreover, with pretreatment of Z-ATAD-FMK, the protein level of cleaved caspase-12 and caspase-3 were significantly reduced and the apoptotic level of neural cells was down-regulated accordingly (Figs. 3 and 4). Taken together, our data indicate that the virulence factor ROP18 in T. gondii is involved in the apoptosis of neural cells through the ER stress pathway. Further studies focusing on ROP18-interacting host factors in ER stress pathway will help us to gain new insights into this pathological process.
Our findings here highlight that the virulence factor ROP18 in T. gondii may contribute to neuronal apoptosis through the ER stress-mediated apoptosis pathway, which may contribute to better understanding the possible mechanism of brain pathology during T. gondii infection.
Jones JL, Parise ME, Fiore AE. Neglected parasitic infections in the United States: toxoplasmosis. Am J Trop Med Hyg. 2014;90(5):794–9.
Parlog A, Schluter D, Dunay IR. Toxoplasma gondii-induced neuronal alterations. Parasite Immunol. 2015;37(3):159–70.
Li Z, Zhao ZJ, Zhu XQ, Ren QS, Nie FF, Gao JM, et al. Differences in iNOS and arginase expression and activity in the macrophages of rats are responsible for the resistance against T. gondii infection. PLoS One. 2012;7(4):e35834.
Gao XJ, Zhao ZJ, He ZH, Wang T, Yang TB, Chen XG, et al. Toxoplasma gondii infection in pregnant women in China. Parasitology. 2012;139(2):139–47.
Schmidt M, Sonneville R, Schnell D, Bige N, Hamidfar R, Mongardon N, et al. Clinical features and outcomes in patients with disseminated toxoplasmosis admitted to intensive care: a multicenter study. Clin Infect Dis. 2013;57(11):1535–41.
Reiter-Owona I, Seitz H, Gross U, Sahm M, Rockstroh JK, Seitz HM. Is stage conversion the initiating event for reactivation of Toxoplasma gondii in brain tissue of AIDS patients? J Parasitol. 2000;86(3):531–6.
Prandovszky E, Gaskell E, Martin H, Dubey JP, Webster JP, McConkey GA. The neurotropic parasite Toxoplasma gondii increases dopamine metabolism. PLoS One. 2011;6(9):e23866.
Holliman RE. Toxoplasmosis, behaviour and personality. J Infect. 1997;35(2):105–10.
Daniels BP, Sestito SR, Rouse ST. An expanded task battery in the Morris water maze reveals effects of Toxoplasma gondii infection on learning and memory in rats. Parasitol Int. 2015;64(1):5–12.
Taylor S, Barragan A, Su C, Fux B, Fentress SJ, Tang K, et al. A secreted serine-threonine kinase determines virulence in the eukaryotic pathogen Toxoplasma gondii. Science. 2006;314(5806):1776–80.
Saeij JP, Boyle JP, Coller S, Taylor S, Sibley LD, Brooke-Powell ET, et al. Polymorphic secreted kinases are key virulence factors in toxoplasmosis. Science. 2006;314(5806):1780–3.
Lei T, Wang H, Liu J, Nan H, Liu Q. ROP18 is a key factor responsible for virulence difference between Toxoplasma gondii and Neospora caninum. PLoS One. 2014;9(6):e99744.
Yuan ZG, Zhang XX, Lin RQ, Petersen E, He S, Yu M, et al. Protective effect against toxoplasmosis in mice induced by DNA immunization with gene encoding Toxoplasma gondii ROP18. Vaccine. 2011;29(38):6614–9.
Fentress SJ, Sibley LD. The secreted kinase ROP18 defends Toxoplasma’s border. Bioessays. 2011;33(9):693–700.
Steinfeldt T, Konen-Waisman S, Tong L, Pawlowski N, Lamkemeyer T, Sibley LD, et al. Phosphorylation of mouse immunity-related GTPase (IRG) resistance proteins is an evasion strategy for virulent Toxoplasma gondii. PLoS Biol. 2010;8(12):e1000576.
Du J, An R, Chen L, Shen Y, Chen Y, Cheng L, et al. Toxoplasma gondii virulence factor ROP18 inhibits the host NF-kappaB pathway by promoting p65 degradation. J Biol Chem. 2014;289(18):12578–92.
Zhou J, Gan X, Wang Y, Zhang X, Ding X, Chen L, et al. Toxoplasma gondii prevalent in China induce weaker apoptosis of neural stem cells C17.2 via endoplasmic reticulum stress (ERS) signaling pathways. Parasit Vectors. 2015;8(1):73.
Wang T, Zhou J, Gan X, Wang H, Ding X, Chen L, et al. Toxoplasma gondii induce apoptosis of neural stem cells via endoplasmic reticulum stress pathway. Parasitology. 2014;141(7):988–95.
Unno A, Suzuki K, Xuan X, Nishikawa Y, Kitoh K, Takashima Y. Dissemination of extracellular and intracellular Toxoplasma gondii tachyzoites in the blood flow. Parasitol Int. 2008;57(4):515–8.
Munoz M, Liesenfeld O, Heimesaat MM. Immunology of Toxoplasma gondii. Immunol Rev. 2011;240(1):269–85.
Hunter CA, Roberts CW, Alexander J. Kinetics of cytokine mRNA production in the brains of mice with progressive toxoplasmic encephalitis. Eur J Immunol. 1992;22(9):2317–22.
Sarciron ME, Gherardi A. Cytokines involved in Toxoplasmic encephalitis. Scand J Immunol. 2000;52(6):534–43.
Suzuki Y. Host resistance in the brain against Toxoplasma gondii. J Infect Dis. 2002;185 Suppl 1:S58–65.
Jung BK, Pyo KH, Shin KY, Hwang YS, Lim H, Lee SJ, et al. Toxoplasma gondii infection in the brain inhibits neuronal degeneration and learning and memory impairments in a murine model of Alzheimer’s disease. PLoS One. 2012;7(3):e33312.
Suzuki Y. Immunopathogenesis of cerebral toxoplasmosis. J Infect Dis. 2002;186 Suppl 2:S234–240.
Wilson EH, Hunter CA. The role of astrocytes in the immunopathogenesis of toxoplasmic encephalitis. Int J Parasitol. 2004;34(5):543–8.
Zhang YH, Chen H, Chen Y, Wang L, Cai YH, Li M, et al. Activated microglia contribute to neuronal apoptosis in Toxoplasmic Encephalitis. Parasit Vectors. 2014;7(1):372.
Gatkowska J, Wieczorek M, Dziadek B, Dzitko K, Dlugonska H. Behavioral changes in mice caused by Toxoplasma gondii invasion of brain. Parasitol Res. 2012;111(1):53–8.
Sanchez V, de-la-Torre A, Gomez-Marin JE. Characterization of ROP18 alleles in human toxoplasmosis. Parasitol Int. 2014;63(2):463–9.
Fentress SJ, Steinfeldt T, Howard JC, Sibley LD. The arginine-rich N-terminal domain of ROP18 is necessary for vacuole targeting and virulence of Toxoplasma gondii. Cell Microbiol. 2012;14(12):1921–33.
El Hajj H, Lebrun M, Arold ST, Vial H, Labesse G, Dubremetz JF. ROP18 is a rhoptry kinase controlling the intracellular proliferation of Toxoplasma gondii. PLoS Pathog. 2007;3(2):e14.
Kaushansky A, Kappe SH. Host ER stress during malaria parasite infection. EMBO Rep. 2015;16(883).
Inacio P, Zuzarte-Luis V, Ruivo MT, Falkard B, Nagaraj N, Rooijers K, et al. Parasite-induced ER stress response in hepatocytes facilitates Plasmodium liver stage infection. EMBO Rep 2015.;16(8):955–964.
Hrincius ER, Liedmann S, Finkelstein D, Vogel P, Gansebom S, Samarasinghe AE, et al. Acute Lung Injury Results from Innate Sensing of Viruses by an ER Stress Pathway. Cell Rep. 2015;11(10):1591–603.
Xu C, Bailly-Maitre B, Reed JC. Endoplasmic reticulum stress: cell life and death decisions. J Clin Invest. 2005;115(10):2656–64.
Ron D, Walter P. Signal integration in the endoplasmic reticulum unfolded protein response. Nat Rev Mol Cell Biol. 2007;8(7):519–29.
Zinszner H, Kuroda M, Wang X, Batchvarova N, Lightfoot RT, Remotti H, et al. CHOP is implicated in programmed cell death in response to impaired function of the endoplasmic reticulum. Genes Dev. 1998;12(7):982–95.
Nakagawa T, Zhu H, Morishima N, Li E, Xu J, Yankner BA, et al. Caspase-12 mediates endoplasmic-reticulum-specific apoptosis and cytotoxicity by amyloid-beta. Nature. 2000;403(6765):98–103.
Yamamoto M, Ma JS, Mueller C, Kamiyama N, Saiga H, Kubo E, et al. ATF6beta is a host cellular target of the Toxoplasma gondii virulence factor ROP18. J Exp Med. 2011;208(7):1533–46.
This work was supported by the National Natural Science Foundation of China (81271864, 30801329), Fok Ying Tung College Young Teachers Fund (131033) of the Chinese Ministry of Education, the Distinguished Young Scholar of Anhui Province (10040606Y19, 1308085JGD11) to Jian Du; the National Natural Science Foundation of China (81471983), the National Basic Research Program of China (973 program; Grant No. 2010CB530001) to Jilong Shen; the Natural Science Foundation of Anhui Educational Committee (No. KJ2015A330) to Wei Wang.
The authors declare that they have no competing interests.
JD and JS conceived and designed the study, and critically revised the manuscript. JD, LW, LG RA, WW and performed the experiments, analyzed the data and drafted the manuscript. JD, LW, RA, WW, LG, MZ, ZJ, YT, YZ, HC, LY and JS participated in the analysis and interpretation of data. All authors have read and approved the final manuscript.
Lijuan Wan and Lingli Gong contributed equally to this work.