Clonorchis sinensis omega-class glutathione transferases play major roles in the protection of the reproductive system during maturation and the response to oxidative stress
© The Author(s). 2016
Received: 6 April 2016
Accepted: 2 June 2016
Published: 13 June 2016
Clonorchis sinensis causes a major food-borne helminthic infection. This species locates in mammalian hepatobiliary ducts, where oxidative stressors and hydrophobic substances are profuse. To adapt to the hostile micromilieu and to ensure its long-term survival, the parasite continuously produces a diverse repertoire of antioxidant enzymes including several species of glutathione transferases (GSTs). Helminth GSTs play pertinent roles during sequestration of harmful xenobiotics since most helminths lack the cytochrome P-450 detoxifying enzyme.
We isolated and analyzed the biochemical properties of two omega-class GSTs of C. sinensis (CsGSTo1 and CsGSTo2). We observed spatiotemporal expression patterns in accordance with the maturation of the worm’s reproductive system. Possible biological protective roles of CsGSTos in these organs under oxidative stress were investigated.
The full-length cDNAs of CsGSTo1 and 2 constituted 965 bp and 1,061 bp with open reading frames of 737 bp (246 amino acids) and 669 bp (223 amino acids). They harbored characteristic N-terminal thioredoxin-like and C-terminal α-helical domains. A cysteine residue, which constituted omega-class specific active site, and the glutathione-binding amino acids, were recognized in appropriate positions. They shared 44 % sequence identity with each other and 14.8–44.8 % with orthologues/homologues from other organisms. Bacterially expressed recombinant proteins (rCsGSTo1 and 2) exhibited dehydroascorbate reductase (DHAR) and thioltransferase activities. DHAR activity was higher than thioltransferase activity. They showed weak canonical GST activity toward 1-chloro-2,4-dinitrobenzene. S-hexylglutathione potently and competitively inhibited the active-site at nanomolar concentrations (0.63 and 0.58 nM for rCsGSTo1 and 2). Interestingly, rCsGSTos exhibited high enzyme activity toward mu- and theta-class GST specific substrate, 4-nitrobenzyl chloride. Expression of CsGSTo transcripts and proteins increased beginning in 2-week-old juveniles and reached their highest levels in 4-week-old adults. The proteins were mainly expressed in the elements of the reproductive system, such as vitelline follicles, testes, seminal receptacle, sperm and eggs. Oxidative stressors induced upregulated expression of CsGSTos in these organs. Regardless of oxidative stresses, CsGSTos continued to be highly expressed in eggs. CsGSTo1 or 2 overexpressing bacteria demonstrated high resistance under oxidative killing.
CsGSTos might be critically involved in protection of the reproductive system during maturation of C. sinensis worms and in response to oxidative conditions, thereby contributing to maintenance of parasite fecundity.
Clonorchis sinensis causes one of the major fish-borne-zoonotic trematodiases. It is prevalent in several countries of Asia, especially where aquaculture systems associated with paddy field are widely used . Approximately 35 million people are infected and another 600 million people are at risk worldwide . Humans are infected by eating raw/undercooked freshwater fish infected with metacercariae. Light infections usually are asymptomatic. However, chronic cumulative infections invoke several hepatobiliary symptoms including cholecystitis, obstructive jaundice, cholangitis and ascites . Pathological alterations like adenomatous hyperplasia and/or dysplasia of the biliary epithelium, mucin-secreting metaplasia and ductal dilatation with fibrosis are frequently observed in those patients . Epidemiological evidence indicates that approximately 10 % of cholangiocarcinoma is associated with chronic C. sinensis infections [5, 6]. Long-standing inflammations accompanied by clonorchiasis might result in oxidative damage of the biliary ductal epithelium and malignant transformation. C. sinensis is classified as a Group 1 biocarcinogen .
Clonorchis sinensis survives more than ten years within the biliary lumen, where free oxygen radicals generated by lipid peroxidation and several hydrophobic substances derived from liver metabolism prevail . In order to adapt to the hostile micromilieu, C. sinensis continuously produces diverse antioxidant enzymes, among which several species of glutathione transferases (GSTs: E.C. 126.96.36.199) are the major components [9, 10]. At least eight proteoforms of mu- and sigma-class GST isozymes have been described. Some are intimately involved in protection of the worm during oxidative stress as well as in neutralization of cytopathic host bile . Nucleotide sequences coding for kappa- (GAA51086) and zeta-type (GAA44819) GSTs have also been identified, but their protein identity and biological properties remain elusive.
GSTs are ubiquitously expressed in almost all known organisms . Typical catalytic activity of GSTs is refined by the conjugation of glutathione (GSH; γ-Glu-Cys-Gly) to a wide variety of non-polar electrophilic, endogenous and exogenous toxic compounds . GSTs play crucial roles against various toxicants, especially in helminth parasites that lack the cytochrome P-450 (CYP450) phase II detoxification enzyme. Most helminth GSTs can be classified into mu- and sigma-types [10, 13], although some GSTs demonstrate mosaic patterns of substrate/inhibitor specificity .
Omega-class GST (GSTo) is a relatively ancient cytosolic enzyme, but is the most recently characterized [11, 15]. A RNA polymerase-related protein designated stringent starvation protein A (SspA) represents a bacterial GST-like molecule due to its highly comparable structural property with GSTo, but lacks GST activity . GSTo has interesting features compared with the other types of GSTs. GSTo has distinct enzymatic properties, e.g. GSH-dependent thioltransferase and dehydroascorbate reductase activity (DHAR), which might be attributable to its structural similarity to glutaredoxin . GSTo shows high affinity toward S-hexylglutathione and 4-nitrophenyl acetate, but low affinity to 1-chloro-2,4-dinitrobenzne. GSTo utilizes cysteine residue to form a mixed disulfide bond with GSH, while most GSTs possess typical tyrosine or serine residues in the active site . GSTo participates in modulation of calcium channels, interaction with cytokine inhibitory drugs, multistep biotransformation, signaling pathway during c-Jun N-terminal kinase (JNK)-mediated apoptosis, sequestration of hydrophobic substances/byproducts generated through diverse hepatic metabolisms and cellular protection from oxidative damages [17–19].
Helminth GSTos have been characterized from Onchocerca volvulus and Schistosoma mansoni, Fasciola spp. and free-living Caenorhabditis elegans [13, 20–22]. Onchocerca volvulus GSTo (Ov-GST3) identified by differential display RT-PCR demonstrates stress resistant effects . Introduction of double-stranded RNA of the Ov-GST3 into mutant C. elegans induces resistance to oxidative stress . Transgenic C. elegans that overexpress GSTo (GSTO-1) exhibits increased resistance during oxidative injuries .
In our previous study involving proteome analysis of C. sinensis GSTs, we observed that CsGSTos were inducible during stimulation of the worm with bile juice . This result prompted us to further characterize biochemical features and biological functions relevant to the CsGSTos in response to oxidative stress. In this study, we characterized biochemical properties of two species of C. sinensis GSTos. We demonstrated that expression profiles of the CsGSTos were spatiotemporally regulated in accordance with the maturation of the worm’s reproductive system. We subsequently investigated possible biological protective roles of CsGSTos in these organs under oxidative stressful conditions.
Clonorchis sinensis metacercariae were collected from naturally infected freshwater fish (Pseudorasbora parva) in an endemic area in Korea. Each of 100 metacercariae was orally infect to Sprague–Dawley rats. Worms were collected from the bile ducts at 1-, 2-, 3- and 4-weeks post-infection. Worms were washed more than 10 times with phosphate buffered saline (PBS, 100 mM, pH 7.4) at 4 °C and were stored at -80 °C. Fresh intact worms were immediately used in ex vivo stimulation experiments (see later section). Animals were housed in accordance with guidelines from the Association for the Assessment and Accreditation of Laboratory Animal Care.
Cloning of C. sinensis omega-class GSTs
Expressed sequence tags (ESTs) were constructed through a screening of randomly selected clones of an adult C. sinensis cDNA library [25, 26]. The similarity patterns of the EST sequences were analyzed against the non-redundant database using BLASTX at the NCBI (http://www.ncbi.nlm.nih.gov). Clones showing high-level sequence identities with Schistosoma mansoni (AAO49385) and Fasciola hepatica (JX156880) GSTo were selected. The adult C. sinensis cDNA library was screened by PCR using vectors (T3 and T7 promoter primers) and gene-specific primers, which contained BamHI (forward) and XhoI (reverse) restriction sites: CsGSTo1-forward (5′-CCG GAT CCA TGC CAA CCT GTT CCA AGC ATT TGC-3′); CsGSTo1-revese (5′-GGC TCG AGT TAC ATG TCC CAG TCA GGA TGA CCA-3′); CsGSTo2-forward (5′-ATG GAT CCA TGT GCT ATC TGG GAG ACG CAG GGA-3′); and CsGSTo2-reverse (5′-GGC TCG AGC TAG GCA ATT TCA AGA TTT GGC TTT CCA GC-3′). T7 promoter and forward primers were used to amplify the 3′-region; the T3 promoter and reverse primers were used for amplification of the 5′-region. The PCR thermal cycler profile included 35 cycles at 94 °C (50 s), 58 °C (50 s) and 72 °C (90 s), followed by 10 min extension at 72 °C. Amplicons were purified using a QIAquick PCR purification kit (Qiagen, Valencia, CA, USA), digested with respective enzymes and cloned into the pET-28a(+) vector (Novagen, Madison, WI, USA). The plasmids were transformed into Escherichia coli DH5α. Nucleotide sequences were determined from both strands. Two full-length cDNAs were obtained by overlapping the 5′- and 3′-region sequences.
The coding profiles and the homology patterns were analyzed using the ORF Finder and BLAST programs. The functionally/structurally conserved domains were searched using ProfileScan (http://myhits.isb-sib.ch/cgi-bin/motif_scan). The secondary structure elements were predicted by the Jpred (www.compbio.dundee.ac.uk/jpred/). The Expasy-Sib Bioinformatics Resource Portal (http://web.expasy.org/compute_pi/) was used to predict the theoretical molecular mass (M r) and isoelectric point (pI). Tertiary structures were simulated using Swiss-PdbViewer (ver4.1) based on the human omega-class GST1 and 2 (pdb 1EEM and 3QAG). The amino acid sequences were employed as queries during sequence analyses using Hidden Markov Models (InterProScan, http://www.ebi.ac.uk/InterProScan/). The retrieved amino acid sequences were aligned with ClustalX 2.1 and optimized with GeneDoc (ver2.7) . The phylogenetic tree was constructed using the neighbor-joining method and the molecular evolution genetics analysis (MEGA) ver5.1 software . Statistical significance of each branching node was examined by a bootstrap analysis of 1000 replicates using SEQBOOT in the PHYLIP package .
Expression of recombinant proteins
The full-length CsGSTo1 and 2 cloned into the pET-28a(+) vector were introduced into E. coli BL21 (DE3). Expression of the recombinant proteins was induced with 0.1 mM isopropyl-β-D-thiogalactopyranoside (IPTG) for 4 h at 37 °C. Bacterial cells were sonicated and rCsGSTos were purified by Ni-nitrilotriacetic acid (NTA) affinity column (Qiagen) using Tris-HCl (50 mM, pH 8.0) supplemented with NaCl (200 mM) and imidazole (250 mM). His-tag was removed by a Thrombin CleanCleave kit (Sigma-Aldrich, St. Louis, MO, USA). The purified proteins were dialyzed against PBS (100 mM, pH 7.4) for 4 h at 4 °C, concentrated by lyophilization and analyzed by 12 % SDS-PAGE under reducing conditions.
Specific antibodies and immunoblotting
Specific-pathogen free 6-week-old female BALB/c mice were immunized with rCsGSTos (100 μg each) emulsified with 2 % ammonium hydroxide gel adjuvant (InvivoGen, San Diego, CA, USA). Two weeks later, proteins (100 μg) mixed with emulsifier were boosted three times at one-week intervals. One-week later, the sera were collected and IgG fractions were purified using Protein G affinity column.
Proteins were separated by 12 % reducing SDS-PAGE and/or isoelectrically focused using IPG strips (pH 3–10; GE Healthcare, Piscataway, NJ, USA) for 30 kVh followed by 12 % SDS-PAGE (2-dimensional electrophoresis; 2-DE). The proteins transferred onto nitrocellulose membranes (Santa Cruz Biotechnology, Santa Cruz, CA, USA) were blocked in Tris-buffered saline (100 mM, pH 8.8) containing 0.01 % Tween 20 and 3 % skim milk for 1 h. The membrane was incubated with anti-rCsGSTo1 or 2 antibody (1:1,000 dilution) overnight at 4 °C and subsequently with horseradish peroxidase (HRP)-conjugated goat anti-mouse IgG (1:4,000 dilution; Cappel, West Chester, PA, USA) for 2 h. Signals were detected using West-Q Pico enhanced chemiluminescence (ECL) kit (GenDEPOT, Dallas, TX, USA). All images were obtained after 2 min exposure for quantitative analysis.
Binding characteristics of native CsGSTos against S-hexylglutathione (SHG) and GSH
We determined the binding specificity of CsGSTos toward SHG and GSH. Adult C. sinensis (4-week-old) were homogenized with a Teflon-pestle homogenizer in PBS (100 mM, pH 7.4) containing a protease inhibitor cocktail (Roche, Basel, Switzerland). The supernatant was obtained by centrifugation at 20,000 g for 30 min at 4 °C. Proteins (200 μg protein per column) were loaded onto a SHG-agarose column (Sigma-Aldrich) or a glutathione-Sepharose 4B column (GE Healthcare). The columns were washed with 20 bed volumes of Tris-HCl buffer (50 mM, pH 7.8) containing 200 mM NaCl. Bound proteins were eluted using Tris-HCl buffer (50 mM, pH 7.8) with 0, 2 and 4 mM step-wise gradient fashions of SHG or GSH. Purified proteins resolved by 12 % SDS-PAGE/2-DE were stained with Coomassie brilliant G-250 (CBB) or further processed by immunoblotting using anti-rCsGSTo1 and 2 antibodies.
GST activity was spectrophotometrically determined employing a panel of substrates (Sigma-Aldrich); 1-chloro-2,4-dinitrobenzene (CDNB; pH 6.5, 340 nm), 1,2-dichloro-4-nitrobenzene (DCNB; pH 7.5, 345 nm), 4-nitrobenzyl chloride (4-NBC; pH 6.5, 310 nm), 4-nitrophenyl acetate (4-NPA; pH 7.0, 400 nm), 4-hydroxy nonenal (pH 7.5, 340 nm), cumene hydroperoxide (CHP; pH 6.5, 340 nm) and ethacrynic acid (pH 7.5, 340 nm). The reactions were recorded for 5 min at 25 °C in 100 mM potassium phosphate buffer (pH 7.2, 200 μl) containing each 4 mM substrate and 4 mM GSH. The formation of ascorbate by the glutathione-dependent DHAR was detected in potassium phosphate buffer (50 mM, pH 7.2) supplemented with 1 mM GSH and 0.25 mM dehydroascorbate (DHA) at 265 nm. Thioltransferase activity was assayed using hydroxylethyl disulfide (HEDS, 2 mM) in potassium phosphate buffer (50 mM, pH 7.2) containing 0.2 mM NADPH, 0.5 mM GSH and 0.5 units of glutathione reductase for 2 min at 340 nm. One unit of enzyme activity was defined as the amount of enzyme that catalyzed the formation of one micromole of product per min in the presence of respective substrates. Vmax and appKm were determined by one site saturation assays in ranges of 0.01–5 mM DHA with 5 mM GSH (saturating concentration). We also used variable concentrations of GSH ranging from 0.01–5.0 mM with 5 mM DHA (saturating concentration). Enzyme activity was monitored by changes of absorbance and was converted to specific activity using a molar extinction coefficient (Δε = 5.3). Non-enzymatic reaction was concomitantly monitored and subtracted from the entire reaction rate. All enzyme assays were independently performed in triplicate at 25 °C. Data were analyzed by best fit algorithm in SigmaPlot10.0.1 (Systat, San Diego, CA, USA).
We determined the inhibition mode of CsGSTos employing SHG and the anthelminthic drug praziquantel (PZQ; Shinpoong, Seoul, Korea). rCsGSTos (each 100 ng) were preincubated with Dulbecco’s PBS supplemented with 10–1,000 μM PZQ or 10–500 nM SHG for 2 min, after which the reaction was initiated by adding 1 mM GSH and 1 mM DHA. The increase in absorbance of the resulting GSH conjugate was recorded spectrophotometrically at 265 nm. A set of reactions under identical conditions was done for each inhibitor concentration and for controls. To examine the inhibition mode of the specific active sites, rCsGSTos were preincubated with saturating concentrations of GSH (5 mM) for 2 min before conjugating reaction with varying concentrations of DHA (0.01–5 mM). To determine the inhibition mode of inhibitors against the ligand-binding sites, the initial velocity of the enzyme reactions was observed in the presence of the respective inhibitors. rCsGSTos were incubated with saturating concentration of DHA (5 mM) for 2 min prior to the conjugating reaction with variable concentrations of GSH (0.01–5 mM). All measurements were independently done in triplicate. Data were analyzed by Lineweaver-Burk plots.
In order to observe tissue distribution pattern of CsGSTos, immunohistochemical staining was done on adult worm sections. Clonorchis sinensis adult worms were fixed in 4 % neutral paraformaldehyde, embedded and cut into 4 μm-thick pieces. Sections were treated with 3 % H2O2 for 10 min, subsequently with Tris buffered saline (100 mM, pH 8.0) containing 3 % BSA and 0.1 % Tween 20 (TBS/T-BSA) for 1 h. The slides were incubated overnight at 4 °C with anti-rCsGSTo1 or 2 antibody (1:400 dilution in TBS/T-BSA) and further incubated with HRP-conjugated goat anti-mouse IgG (1:1,000 dilution; Cappel). Color reactions were developed using HighDef blue chromogen (Enzo Life Sciences, Farmingdale, NY, USA) with PBS (100 mM, pH 7.4) supplemented with 0.05 % 3,3'-diaminobenzidine blue and 0.015 % H2O2 for 5 min. The images were photographed under a TissueFAXS plus (TissueGnostics, Vienna, Austria).
In vitro induction of CsGSTos under oxidative stresses
To assess biological reactivity of CsGSTos under oxidative stressful conditions, we observed induction profile of CsGSTos upon treatment with oxidizing chemicals. Fresh intact worms were stabilized for 1 h at 37 °C in 5 % CO2 atmosphere in serum- and phenol red-free RPMI medium. Worms (10 worms per group per 1 ml of medium) were transferred into fresh medium containing different doses of 5-hydroxy-1,4-naphthoquinone (Juglone; Sigma-Aldrich) (25–100 μM) or CHP (1–4 mM) and incubated for 1 h at 37 °C. The worms were harvested and fractionated into individual compartments, such as seminal receptacle, vitelline follicle-enriched parenchyma and eggs under a dissecting microscope. The conditioned medium containing excretory-secretory products (ESP) was also harvested and sperm was separately collected under a dissecting microscope. Proteins of the respective compartments were extracted in PBS (100 mM, pH 7.4) containing protease inhibitor cocktail (one tablet/25 ml PBS; Complete; Roche) and centrifugation at 12,000 g for 30 min at 4 °C. The proteins prepared from individual compartments of the worms incubated without oxidative treatment were used as controls. Cs tubulin (CsTub; DF143021), whose constant expression was verified by RT-PCR in association with the respective stimuli , was used as an internal control. The induction profiles of CsGSTo were examined by immunoblotting probed each with anti-rCsGSTo antibody and by quantitative real-time RT-PCR (qRT-PCR). At least three independent experiments were done with freshly prepared worms.
Reverse-transcription PCR (RT-PCR) and qRT-PCR
Expressional changes of CsGSTos during worm’s maturation and in response to oxidative stress were determined. Total RNA was extracted from the experimental worms or different developmental stages of worms using a RNeasy Mini kit (Qiagen). RT-PCR was done using RT-PCR PreMix kit (iNtRON, Seongnam, Korea). CsGSTo transcripts in the RNA (1 μg) were amplified by PCR with the following primers: CsGSTo1-forward, 5′-GTT TCC ATT TGT GGA C-3′ and -reverse, 5′-TGG TAG CTG CAA TAC G-3′; CsGSTo2-forward, 5′-TCG TTT GAG CGA ATC G-3′ and -reverse, 5′-CAG CGA GAC TGA GTT G-3′. The thermal cycler profile included pre-heating for 30 min (50 °C) and 10 min (94 °C), 30 cycles of 50 s (94 °C), 50 s (58 °C) and 90 s (72 °C) with 10 min final extension (72 °C). CsTub gene was amplified by PCR using gene specific forward primer (5′-ATT CAG CTG TCC TGG GAA AC-3′) and reverse primer (5′-ACT GCA TTG ATA ACG AAG CG-3′). Thermal cycler profile included 25, 30 and 35 cycles at 94 °C (50 s), 58 °C (50 s) and 72 °C (90 s), followed by a 10 min final extension (72 °C). The PCR products were analyzed on 1 % agarose gel with ethidium bromide staining.
qRT-PCR was conducted using the Rotor-Gene SYBR Green PCR kit and the Rotor-Gene Q Real-time PCR (Qiagen). Total RNA (200 ng) treated with DNase was reverse-transcribed into cDNA using the SuperScript First-Strand Synthesis system (Thermo Fisher Scientific, Waltham, MA, USA). The cDNA was used to examine transcriptional activities of CsGSTo genes with specific primers. The qRT-PCR program included pre-denaturation for 5 min (94 °C), 40 cycles of amplification (94 °C for 15 s, 60 °C for 30 s and 72 °C for 30 s) and a melt cycle from 65 °C to 95 °C. Control reactions were done with RNAs that had not been reverse-transcribed. mRNA abundancy was evaluated in three independent samples for each group with three technical repeats. Data were normalized against those obtained with the CsTrop (ΔCT). Fold inductions (ΔΔCT) were calculated by comparison of the non-stimulated controls. Data were analyzed with the Rotor-Gene Q ScreenClust HRM software using the 2-ΔΔC T method .
Disc diffusion, cell growth and survival assays under oxidative stress
We observed effects of CsGSTo during oxidative stress employing CsGSTo overexpressing E. coli. Escherichia coli BL21 cells transformed with CsGSTo1 or 2 expression plasmids and control cells transformed with mock vector were induced for expression of recombinant proteins by adding 0.1 mM IPTG for 4 h at 37 °C. The cells (5 × 108) were cultured on LB-kanamycin agar plates for 1 h at 37 °C. Discs (6-mm diameter) soaked with 10, 50, 100 and 200 mM CHP or Juglone were placed on the surface of the top agar. The cells were grown for 24 h at 37 °C and the inhibition zones were measured. For the cell growth assay, stationary-phase cultures of CsGSTo overexpressing bacterial cells and control cells were diluted to an optical density at 600 nm (OD600) of 0.004. The cells were grown in LB broth at 37 °C until exponential phase (OD600 = 0.13–0.14). Aliquots were treated with 0.5, 1, 2 and 4 mM CHP. Growth curves were obtained in Erlenmeyer flasks at 37 °C and 225 rpm. The cultures were diluted to an OD600 of 0.01 in LB broth and the OD600 was measured every 1 h for 25 h. CsGSTo overexpressing E. coli and control cells were subjected to survival assay against oxidative injury. Escherichia coli cells were grown (OD600 = 0.5), after which CHP or Juglone was added to exponential bacterial suspension (0.49 ml) to final concentrations of 1, 2 and 4 mM, and incubated from 20–60 min at 37 °C with shaking (225 rpm). The cells diluted in PBS (0.1 ml) were plated on LB agar plates and grown for 48 h at 37 °C. Cell viability was determined by counting colony-forming units per ml (CFU/ml) as percentage of surviving cells compared to untreated cells. The limit of detection was 100 CFU/ml. All assays were done independently in triplicate.
Data are expressed as mean ± standard deviation (SD) of 3–5 independent experiments. Statistical significance was evaluated by a one-way analysis of variance (ANOVA) using the Statistical Package for the Social Sciences (SPSS; ver20.0) software (SPSS, Chicago, IL, USA), or Student’s t-test followed by a Bonferroni correction, as appropriate. Differences in mean values were considered statistically significant at P < 0.05.
Molecular properties of C. sinensis omega-class GSTs
When we simulated the tertiary structure of these proteins, a cysteine residue that constituted the omega-class specific active site (C30/C36) and glutathione binding amino acids were recognized in appropriate positions (K57/K63, V70/V76, E83/E88 and S84/S89 for CsGSTo1 and 2, respectively). Despite their low sequence identity (14.8–18.8 %), the general topology of CsGSTos was comparable with human GSTos. However, additional amino acids between α4 and α5 helices (CsGSTo1) and the N-terminal extension (CsGSTo2) were observed. These appeared to form an unstructured loop-like domain, which could not be readily determined (Additional file 1: Figure S1a, b).
The genomic structure of CsGSTo1 spanned 11.2 kbp with five exons and four intervening introns. The 19.1 kbp-coding DNA sequence of CsGSTo2 was split into four exons by three introns. In contrast to C. sinensis, only a single gene orthologous to CsGSTo was retrieved from another platyhelminth examined. Comparison of CsGSTo structures with other trematode orthologues revealed that a single exon corresponding to the fourth exon of CsGSTo1 was deleted in CsGSTo2. The first three introns of the trematode genes were also detected in the orthologous positions of human genes. However, the fourth intron (and fifth intron of human genes) appeared to be introduced after the divergence of Chordata (Fig. 1b). A phylogenetic tree suggested that the GSTo gene has undergone donor organism-specific duplication event(s), at least in C. sinensis and some kinds of ecdysozoan invertebrates (Fig. 1c).
Identification of the native CsGSTos by affinity binding to SHG and GSH
Enzymatic properties of CsGSTos
Substrate specificity of recombinant CsGSTo1 and 2
Specific activity (μmol/min/mg)
0.13 ± 0.01
0.10 ± 0.02
μ and θ
1.42 ± 0.09
0.91 ± 0.24
0.84 ± 0.06
0.73 ± 0.04
α and θ
1.16 ± 0.04
1.08 ± 0.02
1.08 ± 0.02
0.76 ± 0.02
Kinetic parameters for recombinant CsGSTo1 and 2
1.16 ± 0.02 μmol/min/mg
1.08 ± 0.02 μmol/min/mg
0.21 ± 0.02 mM
0.17 ± 0.02 mM
2.29 × 103/s/M
2.53 × 103/s/M
0.42 ± 0.04 μmol/min/mg
0.56 ± 0.08 μmol/min/mg
0.19 ± 0.02 mM
0.16 ± 0.02 mM
1.95 × 103/s/M
2.06 × 103/s/M
Inhibition characteristics of SHG and PZQ against hydrophobic ligand- and glutathione-binding sites
Spatiotemporal expression patterns of CsGSTos
Induction profile of CsGSTo transcripts and proteins under oxidative stress
Worms treated with different dosages and for different times were harvested and fractionated into individual compartments including eggs, vitelline follicle-enriched parenchyma (see also Fig. 5a), seminal receptacle and sperm. When Western blots were probed with anti-rCsGSTo1 and 2 antibodies, expression of CsGSTo1 and 2 proteins appeared to gradually increase in a dose-dependent manner (Fig. 6c). Expression of CsGSTo proteins was also upregulated with time-lapse, although the increasing levels could not be precisely determined (Fig. 6d). Interestingly, expression of GSTos in eggs was maintained at high levels regardless of oxidative conditions and was not changed by the oxidizing chemicals. Expressional changes of CsμGST2 (negative control) were not observed as previously reported . Cs tubulin, employed as an internal control, also did not show expressional changes (Fig. 6c, d).
CsGSTo overexpressing E. coli are resistant to oxidant-mediated killing
In this study, we characterized two GSTos of the liver fluke C. sinensis, which showed differential and common biochemical and biological properties compared to GSTos characterized from other organisms including helminths. Analysis of expression patterns of CsGSTos in different developmental stages and in the different tissues, as well as induction profile under oxidative stress might provide information regarding functional relevance of CsGSTos. Uniquely, C. sinensis possessed two paralogous GSTos, which exhibited high activity toward 4-NBC, a mu- and theta-specific substrate. Expression of CsGSTos was spatiotemporally regulated along with maturation of the worm’s reproductive system. They showed augmented expression in response to oxidative stress and suggested their biological impacts in those organs. CsGSTos shared common properties with other homologous enzymes, such as harboring a cysteine residue in the active sites to form mixed disulfide bond , DHAR/thioltransferase activities and a weak conjugation activity of CDNB. They also have high binding affinity to SHG [15, 20]. Finally CsGSTos are involved in the cellular protection from oxidative injuries [22, 24].
Interestingly, rCsGSTos demonstrated a significant enzyme activity toward the mu- and theta-class specific substrate, 4-NBC (Table 1). GSTos characterized from humans (HsGSTO1), Anopheles cracens (AcGSTO1-1) and S. mansoni (SmGSTO) do not show conjugating activity toward 4-NBC [20, 33, 34]. A recent genome-wide survey of C. sinensis revealed that the parasite harbors 12 GST species, which are segregated to omega-, zeta-, mu-, sigma-, kappa- and membrane-associated protein in eicosanoid and glutathione metabolism-like protein (MAPEG), while theta-class GST(s) is not recognized [10, 20, 35]. Previous studies involving proteome analyses of F. hepatica and F. gigantica GSTs also did not identify theta-class GST [13, 36]. During a survey of S. mansoni, F. hepatica and F. gigantica genomes in the wormbase (http://parasite.wormbase.org/), we could not retrieve theta-class GST(s). Since GSTos are thought to be intermediate forms, which have occurred between the ancient glutaredoxin and the later mu- and/or theta-class GSTs , it is possible that CsGSTos might retain additional roles substitute for theta-class GST. Alternatively, CsGSTos might acquire a unique catabolic activity toward 4-NBC to ensure its physiological demand during maturation of the reproductive system. To address this intriguing phenomenon, identification of theta-class GST(s) and 4-NBC catalyzing activity of mu-class GSTs of diverse platyhelminths deserve further studies.
The differential expression of organ-specific transcripts of GSTos in higher mammals and insects has been well described [17, 37–40]. The expression of helminth GSTos also appeared to be determined by expression of tissue- and/or stage-specific transcripts of different helminths including a free-living nematode. Onchocerca volvulus GSTo was expressed only in the egg shell at the morula stage of the developing embryo . Caenorhabditis elegans GSTo was expressed in the intestinal cells of the late embryo and adult hermaphrodite, where intestine specific GATA-type transcription factor (Elt-2) was predominantly expressed . The expression of CsGSTos was detected to begin in 2-week-old juveniles and rapidly increased up to the 4-week-old adult stage and the highest expression was observed in eggs. CsGSTos were abundantly distributed in the reproductive system, such as vitelline follicles, testes, seminal receptacle and eggs (Fig. 5). This result demonstrated that expression regulation of CsGSTos is associated with the development of the reproductive system. The genuine sexual maturation of C. sinensis reportedly begins in juveniles between 1- and 2-weeks . It further implies that functional roles of CsGSTos might also be engaged in a relatively narrow target specific adaptive protection of specialized protein(s)/organ(s) compared to other types of GSTs, e.g. for protection of the worm’s reproductive system.
Clonorchis sinensis is equipped with several antioxidant enzymes to respond to and modulate exogenously- and endogenously-derived harmful byproducts. We previously reported that different species of CsGSTs might have evolved for each of the multiple specialized functions, i.e. Cs28σGST1 and its paralogues might be specialized for detoxification of endogenous toxicants, while Cs28σGST3 and Cs26μGST2 conjugate xenobiotics/hydrophobic substances in extracellular environments [9, 10]. The mu- and sigma-class CsGSTs are distributed in the tegument and parenchymal tissues [9, 41] and the phospholipid hydroperoxide glutathione peroxidases, which also revealed protective functions, are localized on the vitellocytes and eggs . We assessed induction profile of CsGSTs by in vitro exposing the worms with oxidizing chemicals. CsGSTo expression increased in a dose- and a time-dependent fashion in the sperm, seminal receptacle and vitelline follicle-enriched parenchyma, but different induction patterns were evident between the two genes and between the chemicals used. CsGSTo1 showed a response relatively sensitive to Juglone compared to CsGSTo2 (Figs. 6, 7). This result might reflect greater activity of CsGSTo1 in scavenging intracellular reactive oxygen species compared to CsGSTo2, because Juglone is an internal inducer of the generation of superoxide anion from molecular oxygen during aerobic metabolism .
We could not demonstrate whether or not induction of CsGSTos was definitively related with their protective roles. We transformed CsGSTo expressing plasmids into E. coli and observed their effects during oxidative killing activity. CsGSTo overexpressing bacteria exhibited significant resistance, while control cells were vulnerable under oxidative harmful states. This result indicates a pivotal role of CsGSTos in protection of E. coli during oxidative killing. Molecular modeling of the active site of CsGSTos was highly comparable to those of humans (Additional file 1: Figure S1) and suggested strongly that Cys-30 (CsGSTo1) and -36 (CsGSTo2) located nearby the helix α1 might form a mixed disulfide bond with the thiol group of GSH as other GSTos do [15, 43]. Induction of CsGSTos might result in increased glutathionylation and this redox detoxification activity might be responsible for the neutralization of toxic components, thus influencing their protective activity .
CsGSTos expressed in eggs were not induced by oxidative stresses, but were maintained at high levels regardless of environmental conditions (Fig. 6). Eggs play critical roles for preservation and expansion of the species, therefore maintenance of cellular viability is important. However, eggs will contact an extremely hostile external environment featuring aridity/humidity, temperature fluctuation and other physicochemical toxicants when expelled from the parasite. It seems reasonable to consider that sufficient amounts of GSTos might be expressed and accumulate in the egg during intrauterine maturation.
It is noteworthy that C. sinensis expresses two paralogous GSTo genes, while other platyhelminths examined possess a single GSTo orthologous gene (Fig. 1). The gene is also multiplied in mammalians and insects. The duplication event appears to be a lineage- or species-specific exclusion in mammalians, in which the gene might be multiplied into paralogues during an early stage of Chordata (or mammalian) evolution. Currently we have no information regarding the selective pressure that drives the lineage-specific duplication(s) of GSTo genes. Whether the presence of two GSTo isotypes and biological functions specialized to these molecules are unique characteristics of C. sinensis or common in the Opisthorchiidae, is a question that requires further studies.
CsGSTo1 and 2 might constitute target toward development of novel chemotherapeutics by inhibiting worm’s resistance within the biliary lumen, where oxidative stress is harsh. Impediment of egg viability by control CsGSTo activity might lead to interruption of C. sinensis life-cycle. Therefore, CsGSTos might also be exploited as transmission blocking vaccines. It would be especially useful for management of reservoir hosts in the fields, which may ultimately contribute to control of human clonorchiasis.
Our data demonstrate that functional roles of CsGSTos are specialized for protection of the reproductive system during maturation and in response to oxidative stress, thereby contributing to maintenance of parasite fecundity. CsGSTos might be involved in the cellular defense against hostile environments and/or in the target specific adaptive response to maintain cellular redox conditions. The detailed understanding on the parasitic bioactive proteins including their unique biochemical features and action mechanisms might be helpful to establish control strategies, where the parasitic diseases are prevalent.
During the revision process, we recognized that genes showing high sequence identity with CsGSTo1 and 2 were registered in the GenBank database under accession numbers GAA34234 and GAA51230 as glutathione S-transferase omega-1 during analysis of C. sinensis draft genome . These genes share 92 and 98 % sequence identity with CsGSTo1 (KX163088) and 2 (KX163089) characterized in this study. Verification of actual relationships of these genes needs further elucidation.
2-DE, two-dimensional electrophoresis; 4-NBC, 4-nitrobenzyl chloride; 4-NPA, 4-nitrophenyl acetate; CDNB, 1-chloro-2,4-dinitrobenzene; CHP, cumene hydroperoxide; CsGSTo1, Clonorchis sinensis omega-class GST1; CsGSTo2, Clonorchis sinensis omega-class GST2; CsTub, Cs tubulin; DCNB, 1,2-dichloro-4-nitrobenzene; DHA, dehydroascorbate; DHAR, dehydroascorbate reductase; ECL, enhanced chemiluminescence; GSH, reduced glutathione; GST, glutathione transferase; HEDS, hydroxylethyl disulfide; Juglone, 5-hydroxy-1,4-naphthoquinone; PBS, phosphate buffered saline; qRT-PCR, quantitative real-time reverse transcription-PCR; SHG, S-hexylglutathione
This work was supported by the Korea National Research Foundation (2015R1D1A1A02061866).
Experimental conception and design: YK, JGK, YAB; sample collections: JGK, CSA, NYK, HJY, WMS; cloning and expression: NYK, YAB, JGK, CSA; biochemistry: JGK, CSA, NYK, SHK, IK; immunohistochemistry: SHK, JGK, YAB, NYK; data analysis: YK, JGK, IK, YAB; paper writing: YK, JGK. All authors read and approved the final version of the manuscript.
The authors declare that they have no competing interests.
Ethics approval and consent to participate
All protocols including parasite infection, recovery of the worms under intraperitoneal injection of 500 μl of a 7:3 mixture of alfaxlone (10 mg/ml) and xylazine (23.32 mg/ml) and euthanasia of the moribund rats, and production of specific antibodies BALB/c mice, were approved by the Institutional Review Board and conducted in the Laboratory Animal Research Center of Sungkyunkwan University (protocol 2015–3).
Open AccessThis article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.
- Madsen H, Dung BT, The DT, Viet NK, Dalsgaard A, Van PT. The role of rice fields, fish ponds and water canals for transmission of fish-borne zoonotic trematodes in aquaculture ponds in Nam Dinh Province, Vietnam. Parasit Vectors. 2015;8:625.View ArticlePubMedPubMed CentralGoogle Scholar
- Keiser J, Utzinger J. Food-borne trematodiases. Clin Microbiol Rev. 2009;22:466–83.View ArticlePubMedPubMed CentralGoogle Scholar
- Marcos LA, Terashima A, Gotuzzo E. Update on hepatobiliary flukes: fascioliasis, opisthorchiasis and clonorchiasis. Curr Opin Infect Dis. 2008;21:523–30.View ArticlePubMedGoogle Scholar
- Jang KT, Hong SM, Lee KT, Lee JG, Choi SH, Heo JS, et al. Intraductal papillary neoplasm of the bile duct associated with Clonorchis sinensis infection. Virchows Arch. 2008;453:589–98.View ArticlePubMedGoogle Scholar
- Choi D, Lim JH, Lee KT, Lee JK, Choi SH, Heo JS, et al. Cholangiocarcinoma and Clonorchis sinensis infection: a case–control study in Korea. J Hepatol. 2006;44:1066–73.View ArticlePubMedGoogle Scholar
- Shin HR, Oh JK, Masuyer E, Curado MP, Bouvard V, Fang YY, et al. Epidemiology of cholangiocarcinoma: an update focusing on risk factors. Cancer Sci. 2010;101:579–85.View ArticlePubMedGoogle Scholar
- Bouvard V, Baan R, Straif K, Grosse Y, Secretan B, El Ghissassi F, et al. WHO International Agency for Research on Cancer Monograph Working Group. A review of human carcinogens - Part B: biological agents. Lancet Oncol. 2009;10:321–2.View ArticlePubMedGoogle Scholar
- Lechner S, Müller-Ladner U, Schlottmann K, Jung B, McClelland M, Rüschoff J, et al. Bile acids mimic oxidative stress induced upregulation of thioredoxin reductase in colon cancer cell lines. Carcinogenesis. 2002;23:1281–8.View ArticlePubMedGoogle Scholar
- Bae YA, Ahn DW, Lee EG, Kim SH, Cai GB, Kang I, et al. Differential activation of diverse glutathione transferases of Clonorchis sinensis in response to the host bile and oxidative stressors. PLoS Negl Trop Dis. 2013;7:e2211.View ArticlePubMedPubMed CentralGoogle Scholar
- Bae YA, Kim JG, Kong Y. Phylogenetic characterization of Clonorchis sinensis proteins homologous to the sigma-class glutathione transferase and their differential expression profiles. Mol Biochem Parasitol. 2016; 206:46-55. doi:10.1016/j.molbiopara.2016.01.002.
- Sheehan D, Meade G, Foley VM, Dowd CA. Structure, function and evolution of glutathione transferases: implications for classification of non-mammalian members of an ancient enzyme superfamily. Biochem J. 2001;360:1–16.View ArticlePubMedPubMed CentralGoogle Scholar
- Oakley AJ. Glutathione transferase: new functions. Curr Opin Struc Biol. 2005;15:2978–84.View ArticleGoogle Scholar
- Morphew RM, Eccleston N, Wilkinson TJ, McGarry J, Perally S, Prescott M, et al. Proteomics and in silico approaches to extend understanding of the glutathione transferase superfamily of the tropical liver fluke Fasciola gigantica. J Proteome Res. 2012;11:5876–89.PubMedGoogle Scholar
- Torres-Rivera A, Landa A. Glutathione transeferases from parasites: a biochemical view. Acta Trop. 2008;105:99–112.View ArticlePubMedGoogle Scholar
- Board PG, Coggan M, Chelvanayagam G, Easteal S, Jermiin LS, Schulte GK, et al. Identification, characterization, and crystal structure of the omega class glutathione transferase. J Biol Chem. 2000;275:24798–806.View ArticlePubMedGoogle Scholar
- Hansen AM, Gu Y, Li M, Andrykovitch M, Waugh DS, Jin DJ, Ji X. Structural basis for the function of stringent starvation protein a as a transcription factor. J Biol Chem. 2005;80:17380–91.View ArticleGoogle Scholar
- Whitbread AK, Masoumi A, Tetlow N, Schmuck E, Coggan M, Board PG. Characterization of the omega class of glutathione transferases. Methods Enzymol. 2005;401:78–99.View ArticlePubMedGoogle Scholar
- Schmuck EM, Board PG, Whitbread AK, Tetlow N, Cavanaugh JA, Blackburn AC, Masoumi A. Characterization of the monomethylarsonate reductase and dehydroascorbate reductase activities of omega class glutathione transferase variants: implications for arsenic metabolism and the age-at-onset of Alzheimer’s and Parkinson’s diseases. Pharmacogenet Genomics. 2005;15:493–501.View ArticlePubMedGoogle Scholar
- Paul S, Jakhar R, Bhardwaj M, Kang SC. Glutathione-S-transferase omega 1 (GSTO1-1) acts as mediator of signaling pathways involved in aflatoxin B1-induced apoptosis-autophagy crosstalk in macrophages. Free Radic Biol Med. 2015;89:1218–30.View ArticlePubMedGoogle Scholar
- Girardini J, Amirante A, Zemzoumi K, Serra E. Characterization of an omega-class glutathione S-transferase from Schistosoma mansoni with glutaredoxin-like dehydroascorbate reductase and thiol transferase activities. Eur J Biochem. 2002;269:5512–21.View ArticlePubMedGoogle Scholar
- Liebau E, Eschbach ML, Tawe W, Sommer A, Fischer P, Walter RD, et al. Identification of a stress-responsive Onchocerca volvulus glutathione S-transferase (Ov-GST-3) by RT-PCR differential display. Mol Biochem Parasitol. 2000;109:101–10.View ArticlePubMedGoogle Scholar
- Liebau E, Höppner J, Mühlmeister M, Burmeister C, Lüersen K, Perbandt M, et al. The secretory omega-class glutathione transferase OvGST3 from the human pathogenic parasite Onchocerca volvulus. FEBS J. 2008;275:3438–53.View ArticlePubMedGoogle Scholar
- Kampkötter A, Volkmann TE, de Castro SH, Leiers B, Klotz LO, Johnson TE, et al. Functional analysis of the glutathione S-transferase 3 from Onchocerca volvulus (Ov-GST-3): a parasite GST confers increased resistance to oxidative stress in Caenorhabditis elegans. J Mol Biol. 2003;325:25–37.View ArticlePubMedGoogle Scholar
- Burmeister C, Lüersen K, Heinick A, Hussein A, Domagalski M, Walter RD, et al. Oxidative stress in Caenorhabditis elegans: protective effects of the omega class glutathione transferase (GSTO-1). FASEB J. 2008;22:343–54.View ArticlePubMedGoogle Scholar
- Cai GB, Bae YA, Kim SH, Sohn WM, Lee YS, Jiang MS, et al. Vitellocyte-specific expression of phospholipid hydroperoxide glutathione peroxidases in Clonorchis sinensis. Int J Parasitol. 2008;38:1613–23.View ArticlePubMedGoogle Scholar
- Bae YA, Cai GB, Kim SH, Sohn WM, Kong Y. Expression pattern and substrate specificity of Clonorchis sinensis tyrosinases. Int J Parasitol. 2013;43:891–900.View ArticlePubMedGoogle Scholar
- Jeanmougin F, Thompson JD, Gouy M, Higgins DG, Gibson TJ. Multiple sequence alignment with Clustal X. Trends Biochem Sci. 1998;23:403–5.View ArticlePubMedGoogle Scholar
- Nicholas KB, Nicholas HBJ, Deerfield DW. GeneDoc: analysis and visualization of genetic variation. Embnew News. 1997;4:14.Google Scholar
- Felsenstein J. PHYLIP - Phylogeny Inference Package (Version 3.2). Cladistics. 1989;5:164–6.Google Scholar
- Livak KJ, Schmittgen TD. Analysis of relative gene expression data using real-time quantitative PCR and the 2-ΔΔC T method. Methods. 2001;25:402–8.View ArticlePubMedGoogle Scholar
- Swope MD, Sun HW, Klockow B, Blake P, Lolis E. Macrophage migration inhibitory factor interactions with glutathione and S-hexylglutathione. J Biol Chem. 1998;273:14877–84.View ArticlePubMedGoogle Scholar
- Brock J, Board PG, Oakley AJ. Structural insights into omega-class glutathione transferases: a snapshot of enzyme reduction and identification of a non-catalytic ligandin site. PLoS One. 2013;8:e60324.View ArticlePubMedPubMed CentralGoogle Scholar
- Wongtrakul W, Pongjaroenkit S, Leelapat P, Nachaiwieng W, Prapanthadara L, Ketterman AJ. Expression and characterization of three new glutathione transferases, an epsilon (AcGSTE2-2), omega (AcGSTO1-1), and theta (AcGSTT1-1) from Anopheles cracens (Diptera: Culicidae), a major Thai malaria vector. J Med Entomol. 2010;47:162–71.View ArticlePubMedGoogle Scholar
- Board PG. The omega-class glutathione transferases: structure, function, and genetics. Drug Metab Rev. 2011;43:226–35.View ArticlePubMedGoogle Scholar
- Wang X, Chen W, Huang Y, Sun J, Men J, Liu H, et al. The draft genome of the carcinogenic human liver fluke Clonorchis sinensis. Genome Biol. 2011;12:R107.View ArticlePubMedPubMed CentralGoogle Scholar
- Chemale G, Morphew R, Moxon JV, Morassuti AL, LaCourse EJ, Barrett J, et al. Proteomic analysis of glutathione transferases from the liver fluke parasite, Fasciola hepatica. Proteomics. 2006;6:6263–73.View ArticlePubMedGoogle Scholar
- Rouimi P, Anglade P, Benzekri A, Costet P, Debrauwer L, Pineau T, et al. Purification and characterization of a glutathione S-transferase omega in pig: evidence for two distinct organ-specific transcripts. Biochem J. 2001;358:257–62.View ArticlePubMedPubMed CentralGoogle Scholar
- Whitbread AK, Tetlow N, Eyre HJ, Sutherland GR, Board PG. Characterization of the human omega class glutathione transferase genes and associated polymorphisms. Pharmacogenetics. 2003;13:131–44.View ArticlePubMedGoogle Scholar
- Walters KB, Grant P, Johnson DL. Evolution of the GST omega gene family in 12 Drosophila species. J Hered. 2009;100:742–53.View ArticlePubMedGoogle Scholar
- Zhang Y, Yan H, Lu W, Li Y, Guo X, Xu B. A novel omega-class glutathione S-transferase gene in Apis cerana cerana: molecular characterisation of GSTO2 and its protective effects in oxidative stress. Cell Stress Chaperon. 2013;18:503–16.View ArticleGoogle Scholar
- Kang SY, Ahn IY, Park CY, Chung YB, Hong ST, Kong Y, et al. Clonorchis sinensis: molecular cloning and characterization of 28-kDa glutathione S-transferase. Exp Parasitol. 2001;97:186–95.View ArticlePubMedGoogle Scholar
- Inbaraj JJ, Chignell CF. Cytotoxic action of juglone and plumbagin: a mechanistic study using HaCaT keratinocytes. Chem Res Toxicol. 2004;17:55–62.View ArticlePubMedGoogle Scholar
- Winayanuwattikun P, Ketterman AJ. An electron-sharing network involved in the catalytic mechanism is functionally conserved in different glutathione transferase classes. J Biol Chem. 2005;280:31776–82.View ArticlePubMedGoogle Scholar
- Menon D, Board PG. A role for glutathione transferase Omega 1 (GSTO1-1) in the glutathionylation cycle. J Biol Chem. 2013;288:25769–79.View ArticlePubMedPubMed CentralGoogle Scholar