Abstract
The present study aims to investigate the protective effect of quercetin against the joint toxic action induced by the mixture of four organophosphate pesticides (mixture-OPs) (dimethoate, acephate, dichlorvos, and phorate) at their corresponding no observed adverse effect level (NOAEL) using metabonomics. Rats were randomly divided into control, quercetin-treated, mixture-OPs-treated, and quercetin plus mixture-OPs-treated groups. Mixture-OPs and quercetin were given to the rats daily through drinking water and intragastric administration, respectively, for 90 days. The metabonomic profiles of rat urine were analyzed using ultra-performance liquid chromatography–mass spectrometry (UPLC/MS). The 14 metabolites significantly changed in the treatment groups compared with the control group, including the biomarkers of OPs exposure (dimethylphosphate, dimethyldithiophosphate, diethylphosphate) and the metabolites of quercetin (quercetin and isorhamnetina). The intensities of gentisic acid, creatinine, suberic acid, hippuric acid, uric acid, and citric acid significantly decreased, whereas the intensities of 7-methylguanine, estrone sulfate, and cholic acid significantly increased, in the mixture-OPs-treated group compared with the control group (
Introduction
Organophosphate pesticides (OPs) have been widely used in agriculture, industry, public health, and domestic application for several decades because of their highly effective insecticidal activities, low cost, and wide varieties. However, the wide use of OPs causes varying degrees of environmental pollution, including air, drinking water, and soil pollution. 1 People are inevitably exposed to OPs from various sources. The characteristic of exposure is long-term and low-dose oral intake. The primary route of OPs exposure for nonoccupational populations may be through residues in food and/or water and exposure to more than one kinds of OP simultaneously. Therefore, the influence of the mixture of OPs (mixture-OPs) on human health has gained the attention of researchers worldwide.
OPs is a well-known acetylcholinesterase inhibitor, which induced specific neurobehavioral effects, including cholinergic signs of toxicity and behavioral effects. The joint toxic action of OPs has been studied since the 1950s. 2 Until now, numerous studies have demonstrated that the joint toxic action of OPs includes interaction, independence, synergism, antagonism, and additivity. 3 Therefore, reducing the toxicity caused by mixture-OPs has become the focus of researchers.
Quercetin is the major flavonoid that belongs to the class called flavonols. Quercetin is found in many common foods including apples, tea, onions, nuts, berries, cauliflower, cabbage, and other foods. Numerous studies have shown that quercetin can elicit antioxidant, anti-inflammatory, anticancer, immunomodulatory, antiaggregatory, and vasodilating effects in vivo and in vitro. 4,5 Furthermore, epidemiological studies have found that increased intake of quercetin can reduce the incidence of cardiovascular disease, hypertensive disease, cancer, neurodegenerative diseases, and gastrointestinal diseases. 6,7 Some researchers have found that the biological activity of quercetin may be closely related to its antioxidant action. 8 In addition, animal experiments have proven that quercetin protects organisms against the toxicity of exogenous poisons, including OPs, polyaromatic hydrocarbons, and hazardous metals. 9 –12 Uzun et al. 9,10 have investigated the protective effect of quercetin against chlorpyrifos-induced lung toxicity, hepatotoxic, and hematologic changes in rats. The results indicated that quercetin has improved the cellular antioxidant status in the lung and liver. However, these studies have focused mainly on the protective property of quercetin on one or two kinds of pesticides, and mostly concerning the tissues at the organ level rather than at the organism level in rats. Therefore, a systemic approach is necessary to determine whether quercetin protects against the joint toxic action induced by mixture-OPs.
Metabonomics is defined as “the quantitative measurement of the dynamic multiparametric metabolic response of living systems to pathophysiological stimuli or genetic modification.” 13 Low-molecular-weight (typically < 1500 Da) metabolites (e.g. plasma, urine, synovial fluid) are measured to analyze metabolites and identify the biomarkers in biofluids or tissues. The metabolic pathways and networks altered by biological effects can also be identified. Metabonomics couples advanced spectroscopic techniques with multivariate or chemometric analysis to study the global alteration of metabolism in biofluids and tissues in response to external stimulus, offering a promising scientific platform for safety assessment. 13 In our previously studies, 14 –17 the toxicity of dichlorvos, acephate, dimethoate, and phorate were evaluated using metabonomics. The results showed that metabonomics can offer more toxicology information than the conventional toxicology method. Metabonomics has been widely used in many fields owing to its advantages. 18,19
Dichlorvos, phorate, dimethoate, and acephate are the most widely used OPs in agriculture and public health programs in China. Therefore, residues of these insecticides are commonly found in foods. 20 The toxicity of OPs is influenced by several factors, such as nature of OPs, exposure level, and exposure time. Humans are usually exposed to multiple chemicals simultaneously, which may interact with each other. Interaction between these chemicals may enhance or reduce their toxicity. Quercetin is contained in our diet, and may be metabolized by the same pathway as pesticides. 21 Therefore, we speculated that dietary quercetin may affect the toxicity of mixture-OPs at the body metabolism level. In our previously published articles, 22,23 the toxicity of the mixture of four aforementioned OPs at their corresponding no observed adverse effect level (NOAEL) on rats was evaluated using metabonomics. The results suggest that exposure to multipesticides with a common mechanism of toxicity causes joint toxic action. The current study aims to investigate the protection effect of quercetin against the joint toxic action induced by the mixture of the four OPs at the body metabolism level using metabonomics and elucidate the mechanism underlying such protective effect.
Materials and methods
Chemicals and reagents
Dimethoate (purity > 98%), acephate (purity > 95%), dichlorvos (purity > 95%), and the analytical standard of phorate were purchased from Hunan Haili Changde Pesticides Chemical Industry, Ltd. (Hunan, China), Nantong Weilike Chemical Industry, Ltd. (Nantong, China), Hebei Century Insecticides, Ltd. (Hebei, China), and Sigma–Aldrich (St. Louis, Missouri, USA), respectively. 7-Methylguanine (purity ≥ 98.0%), uric acid (purity ≥ 99%), suberic acid (purity ≥ 98.0%), estrone sulfate (purity ≥ 98.0%), citric acid (purity ≥ 99.5%), gentisic acid (purity ≥ 95%), isorhamnetin (purity ≥ 95%), quercetin (purity ≥ 95%), creatinine (purity > 99.7%), and leucine enkephalin (purity ≥ 97.0%) were produced from Sigma–Aldrich. Hippuric acid (purity ≥ 99.0%) and cholic acid (purity ≥ 98.0%) were produced from National Standard Substance Research Center (Beijing, China). The dimethylphosphate (DMP) and diethylphosphate (DEP) standards were obtained from AccuStandard, Inc. (New Haven, Connecticut, USA). Dimethyldithiophosphate (DMTP) (99%, purity) was purchased from Applichem (Germany). Ultra-performance liquid chromatography (UPLC)-grade acetonitrile and methanol were obtained from Dikma Science and Technology, Co. Ltd. (Canada). Distilled water was prepared using a Milli-Q system (Millipore, Billerica, Massachusetts, USA).
Lymphocyte separation medium (TBD) was purchased from Tianjin Haoyang Biological Manufacture Co., Ltd. (Tianjin, China). Dimethylsulfoxide (DMSO) was purchased from Tianjin Bodi Chemical Holding Co., Ltd. (Tianjin, China). Low melting point (LMP) and normal melting point (NMP) agarose were produced from Invitrogen (California, California, USA). Ethylendiaminetetracetice acid (EDTA-Na2), Triton X-100, and hydroxymethylaminomethane (Tris Base) were obtained from Amresco (Solon, OH, USA). All other chemicals, reagents, and buffers were analytical grade supplied by Amresco Llc. (Solon, OH, USA).
Animal treatment
Sixty male Wistar rats aged 6 weeks, 180–220 g, were procured from Vital River Laboratory Animal Technology Co. Ltd. (Beijing, China). The study was approved by the Institute of Zoology Animal and Medical Ethics Committee of Harbin Medical University. The animals were housed in stainless steel mesh cages individually and kept in rooms with 12-h light/12-h dark cycle, with a controlled humidity (50–60%) and temperature (22 ± 2°C). They had access freely to AIN-93 M diets and drinking water.
After acclimatization for 1 week, the rats were randomly divided into six groups (10 animals per group). Average initial body weight of each groups was no significantly different (
Sample collection and preparation
Urine samples were collected over 24 h in the metabolism cages over ice packs at each time point (24 h predose, 4-week, 8-week, and 12-week postdose) and centrifuged at 10,000 r/min for 10 min. The supernatants were transferred to polypropylene tubes and stored at −80°C for the metabonomics study. Prior to analysis, urine sample of each rat was thawed and centrifuged at 12,000 r/min for 10 min. The supernatants were diluted with distilled water in a 1:3 (v/v) ratio, and vortex for UPLC/MS analysis. A representative pooled quality control (QC) sample was prepared by mixing equal volumes of the urine samples from sixty rats to verify the reproducibility and reliability of the data; this was analyzed as every 15th sample throughout the analytical run.
After the rats were anesthetized with pentobarbital via intraperitoneal injection, the blood samples were obtained from the aorta abdominalis. The blood sample was added to a tube containing an anticoagulant (EDTA-K2) and shaken to mix the anticoagulant with the blood. Blood samples, 2 ml, diluted with phosphate-buffered saline (PBS) at 1:1. The diluted blood was layered on the top of TBD solution slowly at 1:1. The leukocytes was removed though sedimentation of erythrocytes after centrifuge at 2000 r/min for 10 min at room temperature. Peripheral blood lymphocyte were isolated at the interface of the TBD layers, and washed twice with PBS. The cell viability of lymphocyte was determined using Trypan Blue staining and found to be about 95%. The final concentration of the lymphocytes was adjusted to 1 −3 × 105 cells/ml for comet assay.
Single-cell gel electrophoresis (SCGE) or comet assay
The suspension of lymphocytes was used in comet assay. The frosted slide were prepared by coating 100 µl of 80% NMP agarose, and the 25 µl lymphocytes suspension was mixed with 75 µl LMP agarose were layered on it. Slides were covered with coverslips and cooled at 4°C for 30 min, then uncovered and immersed in the lysis buffer (2.5 M NaCl, 100 mM Na2EDTA, 10 mM Trizma, 10% DMSO, and 1% Triton X-100, pH 10) for 1.5 h at 4°C in the dark. The next step for comet assay followed by the method of Singh et al. 25
Chromatography
The separation system was performed on a 1.8 µm T3 column (ACQUITY (HSS); Waters Corp., Milford, Massachusetts, USA; 2.1 mm × 100 mm) equipped with a UPLC system (ACQUITY UPLC; Waters Corp., USA). A 2 µl aliquot of each urine supernatant sample was injected into the column at a flow rate at 0.35 ml/min for 18 min kept at 35°C. The UPLC mobile phase was as follows: water containing 0.1% formic acid (solution A) and acetonitrile (solution B). The duration of the elution gradient program was as follows: 0–5 min, 0–20% B; 5–9 min, 20–70% B; 9–12 min, 70–98% B; 12–13 min, 98–70% B; 13–16 min, 70–20% B; 16–18 min, 20–0%B. After chromatographic separation, the eluent was directed to the MS system in split mode.
Mass spectrometry
MS was performed on a Micromass Q-TOF MS (Waters Corp., Manchester, UK), which equipped with electrospray ionization (ESI). The ESI source operated in the negative ion mode (ESI−) and positive (ESI+), and in full scan mode from
Metabolite identification and metabolic pathway analysis
The raw data of UPLC/MS were imported into MarkerLynx Application Manager version 4.1 (Waters Corp., USA) for all the peaks detection, which allowed data deconvolution, alignment, and reduction to create a table of mass and retention time (RT). The data matrix including RT and
The pathways of biomarkers were interpreted using databases, including HMDB and Kyoto Encyclopedia of Genes and Genomes (KEGG, http://www.genome.jp/kegg/), to identify the top altered pathway analysis and visualization. 26 Then, other implicated pathways of biomarkers were further interpreted using references.
Data processing and analysis
Collection parameters of EZinfo software for multivariate statistical analysis were set as follows: RT window, 0.2 min; mass window, 0.05 Da; and mass tolerance, 10 mDa. The high and low mass ranges were set at 1000 Da and 50 Da, respectively, and the initial and final RTs to 0.5 min and 14 min, respectively, minimum intensity, 80; noise elimination level, 6.0; deisotope data, Yes. After being recognized and aligned, low-molecular weight metabolites were presented as chromatographic peaks in the base peak intensity (BPI) chromatograms. Pareto-scaled was used before multivariate statistical analysis to avoid chemical noise. To identify the difference in metabolites between each group, PCA and OPLS-DA models were constructed for all time points. Specificity and sensitivity of the potential biomarkers were assessed using the area under the curve (AUC) of the receiver operator characteristic (ROC) curve.
The data of SCGE were processed by CASP (version 1.2.1; Krzysztof Konca, Poland) software to measure the tail DNA% and tail length, which are the parameters of deoxyribonucleic acid (DNA) damage.
The analysis of one-way analysis of variance was used for statistical analysis by SPSS (version 13.0; Beijing Stats Data Mining Co. Ltd., China), and
Results
Body weight change
The BW of all rats was recorded at each time point and each group (see Additional file 1). The BW noticeably increased in the first 8 weeks, and then slowly increased in the last 4 weeks. The BW change in rats in the experimental groups was not significantly different compared with the time-matched control group at each time point (
Analysis of urine samples
Urine samples from the control and experimental groups were analyzed in both positive and negative ionization modes using UPLC/MS. Representative BPI chromatograms of urine in the positive mode in the control and mixture-OPs-treated groups at week 12 are shown in Figure 1. In this study, PCA was performed using UPLC coupled with electrospray ionization quadrupole time-of-flight MS (UPLC/ESI-Q-TOF/MS) to differentiate the distinct metabolite profiles of urine samples from all groups. An initial overview of the quality of the run, including all the QC injections, was obtained by core plot from PCA in the positive and negative mode (see Additional file 2). The result demonstrated that the data points of QC samples were clustered tightly, which provides some evidence that the UPLC/MS system was stable, and the reliable data were suitable for further statistical analysis. The score plot from PCA in the positive mode at each time point (0–12 weeks) is shown in Figure 2. The PCA scores collected from four different time points demonstrated time- and dose-dependent alterations. Figure 2(a) (0 weeks before treatment) shows that the data points did not separate from each other. Figure 2(b) (4 weeks after treatment) shows that the data points of the treatment groups tend to deviated from those of the control group; Figure 2(c) (8 weeks after treatment) shows that the data points of the treatment groups were completely separated from those of the control group. Figure 2(d) (12 weeks after treatment) shows that the data points of the treatment groups were clearly separated from those of the control group and the clusters of each group were close, however, overlaps were observed between groups PQ1 and P.

The representative BPI chromatograms of urine in positive mode in the control (a) and treatment group (b) at week 12. BPI: base peak intensity.

PCA scores plots derived from the UPLC/ESI-Q-TOF/MS of the urine in positive mode. C (box), control group; P (inverted triangle), mixture-OPs-treated group; PQ1 (dot), low-dose quercetin plus mixture-OPs-treated group; PQ2 (star), high-dose quercetin plus mixture-OPs-treated group; Q1 (triangle), low-dose quercetin-treated group; Q2 (diamond), high-dose quercetin-treated group; (a–d) 0, 4, 8, and 12 weeks after treatment. PCA: principal component analysis; UPLC/ESI-Q-TOF/MS: ultra-performance liquid chromatography coupled with electrospray ionization quadrupole time-of-flight mass spectrometry; OPs: organophosphate pesticides.
Fourteen metabolites exhibited statistically significant differences between the treated and the control groups (Table 1). These metabolites strongly influenced the patterns of each group in the score plot and were responsible for the separation between the control and the experiment groups. Eight of these metabolites were the loading plot of positive ions in the 12th week (Figure 3). The intensity value of these metabolites in the positive and negative modes at the 12th week is shown in Tables 2 and 3, respectively. Among these 14 metabolites, DMP, DMTP, uric acid, hippuric acid, and quercetin were identified both in positive and negative modes. DMP (ESI+ and ESI−), DMTP (ESI+ and ESI−), and DEP (ESI−) were only increased in groups P, PQ1, and PQ2 compared with the control group (
Urinary biomarkers identified after treatment in positive and negative mode.
RT: retention time; DMP: dimethylphosphate; DEP: diethylphosphate; DMTP: dimethyldithiophosphate.

Loading plot based on the urine profiling of the control and treatment group in positive mode at week 12.
Summary of intensity values of biomarkers detected in positive ESI mode in rats urine (mean ± SD,
ESI: electrospray ionization; DMP: dimethylphosphate; DMTP: dimethyldithiophosphate; C: control group; Q1: low-dose quercetin-treated group; Q2: high-dose quercetin-treated group; P: mixture-OPs-treated group; PQ1: low-dose quercetin plus mixture-OPs-treated group; PQ2: high-dose quercetin plus mixture-OPs-treated group; ANOVA: analysis of variance.
aSignificantly different from control rats at
bSignificantly different from control rats at
cSignificantly different from P rats at
dSignificantly different from P rats at
Summary of intensity values of biomarkers detected in negative ESI mode in rats urine (mean ± SD,
ESI: electrospray ionization; DMP: dimethylphosphate; DEP: diethylphosphate; DMTP: dimethyldithiophosphate; C: control group; Q1: low-dose quercetin-treated group; Q2: high-dose quercetin-treated group; P: mixture-OPs-treated group; PQ1: low-dose quercetin plus mixture-OPs-treated group; PQ2: high-dose quercetin plus mixture-OPs-treated group; ANOVA: analysis of variance.
aSignificantly different from control rats at
bSignificantly different from control rats at
cSignificantly different from P rats at
dSignificantly different from P rats at
In the positive mode (ESI+), the intensities of gentisic acid, creatinine, hippuric acid, and uric acid significantly decreased in group P compared with the control group (
In the negative mode (ESI−), the intensities of suberic acid, citric acid, hippuric acid, and uric acid significantly decreased in group P compared with the control group (
Comparing the intensities of the metabolites in the positive and negative modes, groups Q1 and Q2 had no significant difference compared with the control group (
ROC curve analysis
ROC curves were used to appraise the specificity and sensitivity. AUCs have been used as useful indexes for evaluating the diagnostic values of markers. If the AUC is greater than 0.7, the biomarkers are relatively exclusive. 27 ROC curves were plotted by entering data from control verses treated rats in urine. The ROC curve analysis of potential urine biomarker is shown in Additional file 3. The AUC of biomarkers were more than 0.7.
SCGE or comet assay
Comet assay is a versatile and sensitive method for measuring single strand breaks in DNA. The percentage of DNA in the tail (tail DNA%) and tail length, which are the most frequently used parameters in DNA strand breaks, provide a clear indication of the appearance of the comets and linearly related to DNA break frequency over a wide range of damage levels. The results of tail DNA% and tail length in each treatment group are shown in Table 4. These two parameters were not statistically significant change in groups Q1 and Q2 compared with the control group (
DNA damage in each group in rat lymphocyte.
C: control group; Q1: low-dose quercetin-treated group; Q2: high-dose quercetin-treated group; P: mixture-OPs-treated group; PQ1: low-dose quercetin plus mixture-OPs-treated group; PQ2: high-dose quercetin plus mixture-OPs-treated group; ANOVA: analysis of variance.
aSignificantly different from control rats at
bSignificantly different from control rats at
cSignificantly different from P rats at
Discussion
Fourteen metabolites in urine were detected by UPLC/Q-TOF/MS analysis, and the ROC curve analysis of potential urine biomarker indicated its specificity and sensitivity. The changes in the intensities of these metabolites in each group in the positive and negative ESI modes are shown in Tables 2 and 3, respectively. DMP, DMTP, and DEP are metabolites of OPs. 28 In the present study, DMP, DMTP, and DEP were only detected in groups exposed to mixture-OPs (P, PQ1, and PQ2), which showed high sensitivity and can be considered as biomarkers for rats exposed to mixture-OPs. Meanwhile, the intensities of quercetin and isorhamnetin significantly only increased in the quercetin-treated and quercetin plus mixture-OPs-treated groups. Their intensities positively correlated with the dosage of quercetin, in agreement with our previous research. 24
Tables 2 and 3 show that the intensities of estrone sulfate and cholic acid significantly increased, whereas gentisic acid, creatinine, suberic acid, and citric acid significantly decreased in group P compared with the control group (

The pathways in response to mixture-OPs and/or quercetin treatment. Upwards arrow or downwards arrow represents that the intensities of metabolites were significantly increased or decreased in the mixture-OPs-treated group compared with the control group; upwards dashed arrow or downwards dashed arrow represent the intensities of metabolites were significantly increased or decreased in the high-dose quercetin plus mixture-OPs-treated group compared with the mixture-OPs-treated group. OPs: organophosphate pesticides.
The first pathway is involved in fatty acids, energy, and sex hormone metabolism. Cholic acid is a major primary bile acid produced in the liver and usually conjugated with glycine or taurine. Bile acids are physiological detergents that facilitate excretion, absorption, and transport of fats and sterols in the intestine and liver. Bile acids are also steroidal amphipathic molecules derived from the catabolism of cholesterol. In the present study, cholic acid significantly increased in group P (
Citric acid is an organic acid that exists in various foods such as fruits, vegetables, and milk, and commonly found as a natural metabolite in living cells. Under aerobic conditions, citric acid is dissimilated by the tricarboxylic acid cycle (TCA cycle). Citric acid metabolism contributes to energy production by providing a major alternative pathway for NAD+ regeneration and allowing acetyl phosphate to yield acetate/adenosine triphosphate (ATP).
32
Creatine is an amino acid that exists in vertebrate tissues and urine. In muscle tissue, creatine occurs as phosphocreatine. Creatine is excreted as creatinine in urine, which functions as part of the cell energy shuttle. The high-energy phosphate group of ATP is transferred to creatine to form phosphocreatine through the following reaction: ATP + creatine ↔ ADP + phosphocreatine. This reaction is reversibly catalyzed by creatine kinase. In this study, the two metabolisms (citric acid and creatinine) significantly decreased in group P compared with the control group (
Estrone sulfate is a sulfated estrone derivative that can be converted to the more active estradiol as needed. Estradiol plays a critical role on reproductive and sexual functioning in women and affects other organs including the bones. A number of studies confirmed that concentrations of the biologically active estrogen, estradiol, are higher in malignant than in normal breast tissue.
34
In this study, the intensities of estrone sulfate significantly increased in group P compared with the control group (
The second pathway is involved in antioxidant system. Gentisic acid is an active metabolite of salicylic acid degradation, which has a broad spectrum of biological activities, including the function of efficiently scavenge hydroxyl radicals. Uric acid is the final oxidation product of purine metabolism and the end product of nitrogen metabolism. Uric acid has strong antioxidant capacity and has been suggested to function as an antioxidant in vivo. In this study, the intensities of uric acid and gentisic acid significantly decreased in group P compared with the control group (
The third pathway is involved in DNA damage. 7-Methylguanine, a metabolite of DNA methylation, is potentially the most useful marker of exposure to endogenous and exogenous methylation agent. The intensities of 7-methylguanine significantly increased in group P compared with the control group (
The fourth pathway concerns the liver and kidney function. Hippuric acid is the glycine-conjugated metabolite of benzoic acid in liver mitochondria used to evaluate liver function.
39
This study shows that hippuric acid significantly decreased in group P compared with the control group (
In this study, the abovementioned biomarker of liver and kidney dysfunction in group PQ2 was alleviated compared with group P (
Although the intensities of the identified metabolites (creatinine, citric acid, estrone sulfate, cholic acid, suberic acid, uric acid, gentisic acid, 7-methylguanine, and hippuric acid) were significantly ameliorated in group PQ2 compared with the group P (
Conclusion
A metabonomics-based systems approach was used to investigate the effect of quercetin on mixture-OPs-induced toxicity in rats. A high dose of quercetin (50 mg/(kg bw day)) has partial protective property on the toxicity induced by chronic exposure to low-level mixture-OPs, including metabolism of fatty acid, energy and sex hormone, antioxidant defense system, DNA damage, as well as liver and kidney function (Figure 4). Antioxidant activity is the most important among these protective effects of quercetin. This study has one limitation that the mechanism of interaction between quercetin and mixture-OPs was not studied. Further research is necessary, especially using molecular biology approaches, to address this point.
The proposed systems toxicity approach based on metabonomic profiling was able to reveal novel information of the protective effect of quercetin against the joint toxic action induced by the mixture of four OPs. Therefore, metabonomics has shown particular promise in the area of toxicology as a systems toxicity approach. However, bottlenecks of metabonomics, such as the lack of extensive and comprehensive metabolite databases, the nonautomated structure assignments, the limitations in data analysis interpretation restrict the throughput and hinder its broader use.
Footnotes
Acknowledgments
The Laboratory of Nutrition and Food Hygiene in Harbin Medical University is the key laboratory of Heilongjiang Province and Heilongjiang Higher Education Institutions.
Declaration of Conflicting Interests
The author(s) declared no potential conflicts of interest with respect to the research, authorship, and/or publication of this article.
Funding
The author(s) disclosed receipt of the following financial support for the research, authorship, and/or publication of this article: Financial support from the National Natural Science Foundation of China (81172672) is gratefully acknowledged.
References
Supplementary Material
Please find the following supplemental material available below.
For Open Access articles published under a Creative Commons License, all supplemental material carries the same license as the article it is associated with.
For non-Open Access articles published, all supplemental material carries a non-exclusive license, and permission requests for re-use of supplemental material or any part of supplemental material shall be sent directly to the copyright owner as specified in the copyright notice associated with the article.
