- Open Access
Study on expression of lncRNA RGMB-AS1 and repulsive guidance molecule b in non-small cell lung cancer
Diagnostic Pathologyvolume 10, Article number: 63 (2015)
The relationships between lncRNAs and tumors have currently become one of the focuses on cancer studies. However, there are a few studies about lncRNAs in non-small cell lung cancer (NSCLC) at present.
Microarray analysis was designed to study the expression patterns of lncRNAs in three pairs of NSCLC tissues. The expression of lncRNA RGMB-AS1 and repulsive guidance molecule b (RGMB) were detected in 72 paired NSCLC tissues and adjacent normal tissues by qRT-PCR assay. The relations of lncRNA RGMB-AS1 and RGMB expression with clinicopathological factors of NSCLC patients were explored. A549 and SPC-A-1 cells were transfected with siRNA of lncRNA RGMB-AS1 and negative control. RGMB expression level was detected by qRT-PCR assay and western blot analysis.
The results of microarray found that 571 lncRNAs were differentially expressed in NSCLC tissues (Fold change cut-off: 5.0, P < 0.05), including 304 upregulated and 267 downregulated lncRNAs. The results of qRT-PCR showed that lncRNA RGMB-AS1 expression was significantly higher in NSCLC tissues than in adjacent normal tissues (P < 0.05), while RGMB mRNA showed an opposite trend (P < 0.05). Correlation analysis indicated that the expression of lncRNA RGMB-AS1and RGMB mRNA were inversely correlated (R2 = 0.590, P < 0.05). While lncRNA RGMB-AS1 and RGMB expression levels in NSCLC tissues were associated with the occurrence of differentiation status, lymph node metastases and TNM stage (P < 0.05). Transfection with siRNA of lncRNA RGMB-AS1, subsequent results showed that RGMB mRNA and protein expression were upregulated (P < 0.05) in A549 and SPC-A-1 cells compared to the control groups.
We identified lncRNA RGMB-AS1 was upregulated and RGMB was downregulated in NSCLC patients. Both were related to differentiation status, lymph node metastases and TNM stage. Studies also indicated that lncRNA RGMB-AS1and RGMB were inversely correlated.
The virtual slide(s) for this article can be found here: http://www.diagnosticpathology.diagnomx.eu/vs/7911587521528276
Lung cancer, one of the most frequent malignancies in the world, is rapidly becoming the main cause of cancer related death nowadays . Approximately 85 % of patients are diagnosed with non-small cell lung cancer (NSCLC). Modifications of chemotherapy combinations, the addition of monoclonal antibodies, including bevacizumab [2, 3] and cetuximab [4, 5], and the incorporation of histologic subtype  into treatment decisions have added to the survival of patients with advanced NSCLC. Unfortunately, therapeutic outcomes appear to have reached a plateau, with response rates of 20 to 35 % and median survival of 8 to 12 months . Currently, encouraging new targeted agents have led to a better understanding of the molecular pathogenesis that causes NSCLC, especially the non-protein-coding portion (~90 %) as non-coding RNA (ncRNA) of the genome . The ncRNAs characterize as three types, long ncRNAs, mid-size ncRNAs and short ncRNAs . Although most studies on ncRNAs are focused on short ncRNAs, such as microRNAs (miRNAs) , long non-coding RNAs (lncRNAs) are rapidly gaining prominence recently.
LncRNAs, tentatively defined as non-coding RNAs more than 200 nucleotides (nt) , are usually divided into exonic lncRNAs, intronic lncRNAs, intergenic lncRNAs and overlapping lncRNAs in accordance with their location relative to the protein-coding transcripts . In recent years, a large number of lncRNAs have been identified and the abnormal expression of lncRNAs has been implicated in imprinting , enhancing various biological functions , X chromosome inactivation, charomatin structure . Thus, lncRNAs are critical for normal development and, in many cases, are deregulated in diseases, such as cancer [16–19].
Currently, studies on lncRNAs in NSCLC are limited. In this study, we analyzed the expression patterns of lncRNAs in three pairs of NSCLC tissues. With performing qRT-PCR assay, we chose some kinds of lncRNAs to test and verify the results of microarray analysis. Combined with the UCSC data-base and bioinformatics analysis, not all lncRNAs from microarray showed potential significance. With the screening tests, we found lncRNA RGMB-AS1 may have potential function in the development of NSCLC. So we focused on the expression of lncRNA RGMB-AS1, and analyzed expression of its potential protein on the base of their location. Then the relationships between the expression of lncRNA RGMB-AS1 and repulsive guidance molecule b (RGMB) and clinicopathological factors of patients with NSCLC were explored. Our results may provide a new perspective for the mechanisms of NSCLC.
Patients and tissue samples
This study was approved by the Ethics Committee of Zhengzhou University. Seventy-two paired NSCLC tissues and adjacent normal tissues (≥3 cm away from tumor) were obtained from patients who received surgical resection of NSCLC between 2012 and 2014 in the First Affiliated Hospital of Zhengzhou University. Of the seventy-two pairs of samples, three were used in lncRNA microarray analysis and all pairs were used for additional evaluations. Patients were diagnosed with NSCLC was confirmed by histopathology. The tumor samples and matched adjacent normal tissues were snap-frozen in liquid nitrogen immediately after resection until total RNA extraction.
Cell culture and transfection
Human NSCLC cell lines (A549 and SPC-A-1) were obtained from the Type Culture Collection of the Chinese Academy of Sciences (Shanghai, China). Both were cultured in RPMI-1640 medium (Gibco BRL, Grand Island, NY) supplemented with 10 % fetal bovine serum, 100 U/mL penicillin, and 100 μg/mL streptomycin at 37 °C in a humidified 5 % CO2 atmosphere. In all experiments, exponentially growing cells were used.
For transfection, cells were seeded into six-well plates at a density of 5 × 104 cells/well. When cells viability reached approximately 80 %, transient transfection was performed using Lipofectamine™ 2000 (Invitrogen, Carlsbad, CA, USA) following the manufacturer’s instructions. siRNA, working on lncRNA RGMB-AS1, and control oligonucleotide (negative control) were synthesized by Shanghai GenePharma Co. Ltd.
For NSCLC tissues, if the proportion of NSCLC cells in a tissue section was >80 % then the frozen block was subjected to RNA extraction. According to the manufacturer's protocol, total RNA was extracted from 72 pairs of snap-frozen NSCLC tissues and adjacent normal tissues using TRIzol reagent (Invitrogen, CA, USA). For NSCLC cell lines, total RNA was extracted with an RNA Extraction Kit (Qiagen). The integrity of the RNA was evaluated by a Nanodrop ND-1000 (Thermo Scientific, Worcester, USA). The value of OD260/280 is around 1.8 as a criterion of acceptable purity.
Of the seventy-two pairs of samples, three pairs of matched samples were used for microarray analysis. RNA was purified from total RNA after removal of rRNA (mRNA-ONLY™ Eukaryotic mRNA Isolation Kit, Epicentre). Then, each sample was transcribed into fluorescent cRNA along the entire length of the transcripts without 3’ bias utilizing a random priming method. The labeled cRNAs were hybridized onto the Human lncRNA Array v2.0 (8 × 60 K, Arraystar).
For qRT-PCR assay, RNA was reverse transcribed to cDNA from 1.0 μg of total RNA using a Reverse Transcription Kit (Thermo Scientific),. Real-time PCR analyses were conducted using the Power SYBR Green (Thermo Scientific). All protocols were carried out according to the manufacturer's instructions. LncRNA RGMB-AS1 (NR_033932) and RGMB mRNA (NM_001012761) expression were determined by qRT-PCR. GAPDH mRNA (NM_002046) served as an endogenous control for normalization. The correlating primers are the following sequences: lncRNA RGMB-AS1, forward: 5’AGTGGGCAAACTTCAACGTTC 3’, and reverse: 5’ GAGCTGCCATTGAATTAATCCG 3’. RGMB, forward: 5’TTCAGGTTCAAGTGA CAAACG-3’ and reverse: 5’-ACTGAACCTGACCGTACATCATCTGTCACAGCT TGGTA -3’. GAPDH, forward: 5’-AGAGGCAGGGATGATGTTCTG-3’, and reverse: 5’-GACTCATGACCACAGTCCATGC-3’. The PCR reaction was executed on an ABI 7500 Fast Real-Time PCR System (Applied Biosystems, CA, USA) at 95 °C for 180 s, followed by 40 cycles of 95 °C for 15 s and 65 °C for 30 s. All experiments were performed in triplicate, and melting curve analysis was performed to ensure the specificity of the products after amplification. The median in each triplicate was used to calculate relative LncRNA RGMB-AS1 and RGMB mRNA concentrations (ΔCt = Ct median lncRNA or mRNA - Ct median GAPDH), which was then converted to x-fold changes(2-ΔCt) .
Western blot analysis
Total proteins were extracted from the transfected cells at 48 h after transfection and then were subjected to SDS-PAGE (SDS/polyacrylamide (10 %) gel electrophoresis) and transferred electrophoretically onto polyvinylidene difluoride (PVDF) filter membranes (Whatman). The membranes were blocked in 5 % skim milk for 1 h, washed four times with Tris-buffered saline containing Tween 20 (TBST) at room temperature, then incubated overnight at 4 °C with diluted primary antibody (rabbit anti-human RGMB antibody, 1:500, Santa Cruz Biotechnology). Following extensive washing with TBST, the membranes were incubated with secondary antibody (goat anti-rabbit IgG, 1:2000, Santa Cruz Biotechnology) for 1 h. After four washes (15 min each) with TBST at room temperature, the immunoreactivity was visualized by enhanced chemiluminescence. The antibody against GAPDH (Santa Cruz Biotechnology) served as an endogenous reference.
SPSS17.0 soft-ware (SPSS, Inc, Chicago, IL) was performed for statistical analysis. Data were expressed as the mean ± the standard deviation (SD). The student’s t-test and a one-way analysis of variance (ANOVA) were used in the comparison of means from different samples. P values of less than 0.05 were considered statistically significant.
Overview of different lncRNA expression profilings in NSCLC tissues relative to adjacent normal tissues
We profiled lncRNAs expression in tumors from NSCLC patients. Comparing to adjacent normal tissues on the basis of lncRNA array results, we found that 571 lncRNAs were differentially expressed in NSCLC tissues (Fold change cut-off: 5.0, P < 0.05), including 304 upregulated and 267 downregulated lncRNAs. Table 1 showed 14 lncRNAs randomly selected among the differential lncRNAs (Log2 fold change, P < 0.05).
Expression of lncRNA RGMB-AS1 and RGMB mRNA in NSCLC tissues
In 72 NSCLC tissues, we performed qRT-PCR to check lncRNA RGMB-AS1 and RGMB mRNA expression, and found that lncRNA RGMB-AS1 expression was significantly higher in NSCLC tissues than in adjacent normal tissues (P < 0.05, Fig. 1a), while RGMB mRNA showed an opposite trend (P < 0.05, Fig. 1b). These data indicated that the expression of lncRNA RGMB-AS1 and RGMB mRNA were inversely correlated (R2 = 0.590, P < 0.05, Fig. 1c).
Association of lncRNA RGMB-AS1 and RGMB expression with clinicopathologic features of NSCLC patients
Using the qRT–PCR analysis and with adjacent normal tissues as reference, the clinicopathologic features of the 72 NSCLC samples included in this study are presented in Table 2. We found that lncRNA RGMB-AS1 and RGMB expression levels in NSCLC tissues were associated with the occurrence of differentiation status, lymph node metastases and TNM stage (P < 0.05; Table 2, Fig. 2). No significant differences were observed between lncRNA RGMB-AS1, RGMB expression and either gender or age (P >0.05; Table 2).
siRNA of lncRNA RGMB-AS1 can promote the expression of RGMB in A549 and SPC-A-1 cells
To further study that lncRNA RGMB-AS1 expression has an opposite correlation with RGMB expression, we studied whether knockdown of lncRNA RGMB-AS1 would have an effect on the expression of RGMB. Transfection with siRNA of lncRNA RGMB-AS1, subsequent qRT-PCR and western blot analysis indeed showed that RGMB mRNA (P < 0.05, Fig. 3a) and protein expression (P < 0.05, Fig. 3b) were upregulated in A549 and SPC-A-1 cells compared to the control groups.
The relationship between lncRNAs and tumors has currently become one of the focuses of cancer studies. In digestive system tumors, lncRNA HNF1A-AS1 has been shown to regulate proliferation and migration in oesophageal adenocarcinoma cells . LncRNA HOTAIR was found overexpression in human hepatocellular carcinomas (HCC) . The highly upregulated in liver cancer (HULC) gene expression is confined to colorectal carcinomas that metastasize to the liver . In hematological system tumors, lncRNA NEAT1 expression has been revealed to impair myeloid differentiation in acute promyelocytic leukemia cells . In urinary system tumors, there have been reports in recent years about lncRNA H19 associating with bladder cancer [25, 26]. In the male reproductive system, prostate cancer gene expression marker 1 (PCGEM1) has been demonstrated as a prostate tissue-specific, and prostate cancer-associated noncoding RNA (ncRNA) gene . In respiratory system tumors, MALAT1 was originally identified as a prognostic marker for metastasis and patient survival in NSCLC, specifically in early stages of lung adenocarcinoma . In this study, we have identified lncRNA RGMB-AS1 is aberrantly expressed in human NSCLC tissues compared to paired adjacent normal tissues. We also found that altered lncRNA RGMB-AS1 expression levels are associated with the occurrence of lymph node metastases, and the differentiation status and TNM stage in NSCLC patients, which may be a hint for NSCLC occurrence and be a potential biomarker for the diagnosis of early NSCLC.
Simultaneously, we also identified RGMB is downregulated in human NSCLC tissues via qRT-PCR analysis. Statistical analysis showed RGMB expression was also related to the occurrence of lymph node metastases, and the differentiation status and TNM stage in NSCLC patients. RGMB, also known as DRAGON, is a member of the repulsive guidance molecules (RGMs) family which consists of RGMA, RGMB, and RGMC . RGMs are a group of cysteine rich 33 kDa proteins, including an N-terminal signal peptide, proteolytic cleavage site, partial von Willebrand factor type D domain, and glycophosphatidylinositol (GPI) anchor [30, 31]. RGM proteins can function as bone morphogenetic protein (BMP) coreceptors [32–34], which are members of the transforming growth factor beta (TGF-β) family of ligands and play a role in many biological activities. Specifically, RGMB directly interacts with BMP receptors for BMP-2 and BMP-4 enhancing binding to their ligands . RGMs work
as BMP co-receptors and by participating in BMP signalling pathway may also be involved in cancer development and progression. In prostate cancer, the knockdown of RGMB significantly enhanced the prostate cancer cell capacity, namely increased growth, adhesive, motility and mobility . Knockdown of RGMB also was studied in breast cancer. The results showed the promotion of growth, survival, adhesion, and migration of breast cancer cells [37, 38].
The UCSC data-base results showed that lncRNA RGMB-AS1 was localized in human chromosome 5 between 98105322 and 98108829 base sites. To check the location of RGMB, UCSC data-base also showed that it was localized in human chromosome 5 between 98109606 and 98129556 base sites. The locations of both lead us to infer that lncRNA RGMB-AS1 and RGMB may have a certain expression control. Applying statistical analysis and performing RNA interference technique in A549 and SPC-A-1 cells, we revealed that lncRNA RGMB-AS1 expression was inversely associated with RGMB expression. The detailed molecular analysis of lncRNA RGMB-AS1 and RGMB still need to be further studied.
In summary, we identified upregulation of lncRNA RGMB-AS1 and downregulation of RGMB in 72 NSCLC patients, and both were related to differentiation status, lymph node metastases and TNM stage. We also provided evidence demonstrating the probable association between lncRNA RGMB-AS1 and RGMB. Based on these findings, we propose that lncRNA RGMB-AS1 and RGMB may serve as a new direction in NSCLC research. But more elaborate studies will be necessary for further exploration of the potential role of lncRNA RGMB-AS1 and RGMB in development of NSCLC in the future.
Non-small-cell lung cancers
Long non-coding RNAs
Repulsive guidance molecule b
Squamous cell carcinoma
Human hepatocellular carcinomas
The highly upregulated in liver cancer
Prostate cancer gene expression marker 1
Bone morphogenetic protein
Transforming growth factor beta
Siegel R, Naishadham D, Jemal A. Cancer statistics, 2013. CA Cancer J Clin. 2013;63(1):11–30.
Le Moulec S, Hadoux J, Gontier E, Chargari C, Helissey C, Lamand V, et al. Combination of paclitaxel and bevacizumab in heavily pre-treated non-small-cell lung cancer (NSCLC) patients: a case series study on 15 patients. Bull Cancer. 2013;100(12):30–7.
Zarza V, Couraud S, Bosc C, Toffart AC, Moro-Sibilot D, Souquet PJ. Paclitaxel-bevacizumab is a possible alternative as salvage chemotherapy in advanced non-small cell bronchial carcinoma. Rev Mal Respi. 2014;31(7):601–7.
Lin CH, Lin MT, Kuo YW, Ho CC. Afatinib combined with cetuximab for lung adenocarcinoma with leptomeningeal carcinomatosis. Lung Cancer. 2014;85(3):479–80.
Sgambato A, Casaluce F, Maione P, Rossi A, Ciardiello F, Gridelli C. Cetuximab in advanced non-small cell lung cancer (NSCLC): the showdown? J Thorac Dis. 2014;6(6):578–80.
Scagliotti GV, Parikh P, von Pawel J, Biesma B, Vansteenkiste J, Manegold C, et al. Phase III study comparing cisplatin plus gemcitabine with cisplatin plus pemetrexed in chemotherapy-naive patients with advanced-stage non-small-cell lung cancer. J Clin Oncol. 2008;26(21):3543–51.
Schiller JH, Harrington D, Belani CP, Langer C, Sandler A, Krook J, et al. Comparison of four chemotherapy regimens for advanced non-small-cell lung cancer. N Eng J Med. 2002;346(2):92–8.
Cufer T, Ovcaricek T, O’Brien ME. Systemic therapy of advanced non-small cell lung cancer: major-developments of the last 5-years. Eur J Cancer. 2013;49(6):1216–25.
Haserlat SM, Supko JG, Haluska FG, Louis DN, Christiani DC, Settleman J, et al. Activating mutations in the epidermal growth factor receptor underlying responsiveness of non-small-cell lung cancer to gefitinib. N Engl J Med. 2004;350(21):2129–39.
Boggon TJ, Naoki K, Sasaki H, Fujii Y, Eck MJ, Sellers WR, et al. EGFR mutations in lung cancer: correlation with clinical response to gefitinib therapy. Science. 2004;304(5676):1497–500.
Soda M, Choi YL, Enomoto M, Takada S, Yamashita Y, Ishikawa S, et al. Identification of the transforming EML4-ALK fusion gene in non-small-cell lung cancer. Nature. 2007;448(7153):561–6.
Derrien T, Johnson R, Bussotti G, Tanzer A, Djebali S, Tilgner H, et al. The GENCODE v7 catalog of human long noncoding RNAs: analysis of their gene structure, evolution, and expression. Genome Res. 2012;22(9):1775–89.
Nagano T, Fraser P. Emerging similarities in epigenetic gene silencing by long noncoding RNAs. Mamm Genome. 2009;20(9–10):557–62.
Kim A, Zhao H, Ifrim I, Dean A. Beta-globin intergenic transcription and histone acetylation dependent on an enhancer. Mol Cell Biol. 2007;27(8):2980–6.
Rinn JL, Kertesz M, Wang JK, Squazzo SL, Xu X, Brugmann SA, et al. Functional demarcation of active and silent chromatin domains in human HOX loci by noncoding RNAs. Cell. 2007;129(7):1311–23.
Silva JM, Perez DS, Pritchett JR, Halling ML, Tang H, Smith DI. Identification of long stress-induced non-coding transcripts that have altered expression in cancer. Genomics. 2010;95(6):355–62.
You QY, Tao H, Ling B. Long Noncoding RNA HOX Transcript Antisense Intergenic RNA (HOTAIR) as a Foe and Novel Potential Therapeutic Target for Endometrial Carcinoma. Int J Gynecol Cancer. 2014;24(9):1536.
Zhao W, Luo J, Jiao S. Comprehensive characterization of cancer subtype associated long non-coding RNAs and their clinical implications. Sci Rep. 2014;4:6591.
Ding X, Zhu L, Ji T, Zhang X, Wang F, Gan S, et al. Long intergenic non-coding RNAs (LincRNAs) identified by RNA-seq in breast cancer. PLoS One. 2014;9(8), e103270.
Schmittgen TD, Livak KJ. Analyzing real-time PCR data by the comparative C (T) method. Nat Protoc. 2008;3(6):1101–8.
Yang X, Song JH, Cheng Y, Wu W, Bhagat T, Yu Y, et al. Long non-coding RNA HNF1A-AS1 regulates proliferation and migration in oesophageal adenocarcinoma cells. Gut. 2014;63(6):881–90.
Ishibashi M, Kogo R, Shibata K, Sawada G, Takahashi Y, Kurashige J, et al. Clinical significance of the expression of long non-coding RNA HOTAIR in primary hepatocellular carcinoma. Oncol Rep. 2013;29(3):946–50.
Matouk IJ, Abbasi I, Hochberg A, Galun E, Dweik H, Akkawi M. Highly upregulated in liver cancer noncoding RNA is overexpressed in hepatic colorectal metastasis. Eur J Gastroenterol Hepatol. 2009;21(6):688–92.
Zeng C, Xu Y, Xu L, Yu X, Cheng J, Yang L, et al. Inhibition of long non-coding RNA NEAT1 impairs myeloid differentiation in acute promyelocytic leukemia cells. BMC Cancer. 2014;14:693.
Luo M, Li Z, Wang W, Zeng Y, Liu Z, Qiu J. Upregulated H19 contributes to bladder cancer cell proliferation by regulating ID2 expression. FEBS J. 2013;280(7):1709–16.
Luo M, Li Z, Wang W, Zeng Y, Liu Z, Qiu J. Long non-coding RNA H19 increases bladder cancer metastasis by associating with EZH2 and inhibiting E-cadherin expression. Cancer Lett. 2013;333(2):213–21.
Fu X, Ravindranath L, Tran N, Petrovics G, Srivastava S. Regulation of apoptosis by a prostate-specific and prostate cancer-associated noncoding gene, PCGEM1. DNA Cell Biol. 2006;25(3):135–41.
Ji P, Diederichs S, Wang W, Böing S, Metzger R, Schneider PM, et al. MALAT-1, a novel noncoding RNA, and thymosin beta4 predict metastasis and survival in early-stage non-small cell lung cancer. Oncogene. 2003;22(39):8031–41.
Severyn CJ, Shinde U, Rotwein P. Molecular biology, genetics and biochemistry of the repulsive guidance molecule family. Biochem J. 2009;422(3):393–403.
Monnier PP, Sierra A, Macchi P, Deitinghoff L, Andersen JS, Mann M, et al. RGM is a repulsive guidance molecule for retinal axons. Nature. 2002;419(6905):392–5.
Niederkofler V, Salie R, Sigrist M, Arber S. Repulsive guidance molecule (RGM) gene function is required for neural tube closure but not retinal topography in the mouse visual system. J Neurosci. 2004;24(4):808–18.
Babitt JL, Zhang Y, Samad TA, Xia Y, Tang J, Campagna JA, et al. Repulsive guidance molecule (RGMa), a DRAGON homologue, is a bone morphogenetic protein co-receptor. J Biol Chem. 2005;280(33):29820–7.
Babitt JL, Huang FW, Wrighting DM, Xia Y, Sidis Y, Samad TA, et al. Bone morphogenetic protein signaling by hemojuvelin regulates hepcidin expression. Nat Genet. 2006;38(5):531–9.
Samad TA, Srinivasan A, Karchewski LA, Jeong SJ, Campagna JA, Ji RR, et al. DRAGON: a member of the repulsive guidance molecule-related family of neuronal- and muscle-expressed membrane proteins is regulated by DRG11 and has neuronal adhesive properties. J Neurosci. 2004;24(8):2027–36.
Samad TA, Rebbapragada A, Bell E, Zhang Y, Sidis Y, Jeong SJ, et al. DRAGON, a bone morphogenetic protein co-receptor. J Biol Chem. 2005;280(14):14122–9.
Li J, Ye L, Kynaston HG, Jiang WG. Repulsive guidance molecules, novel bone morphogenetic protein co-receptors, are key regulators of the growth and aggressiveness of prostate cancer cells. Int J Oncol. 2012;40(2):544–50.
Li J, Ye L, Sanders AJ, Jiang WG. Repulsive guidance molecule B (RGMB) plays negative roles in breast cancer by coordinating BMP signaling. J Cell Biochem. 2012;113(7):2523–31.
Li J, Ye L, Mansel RE, Jiang WG. Potential prognostic value of repulsive guidance molecules in breast cancer. Anticancer Res. 2011;31(5):1703–11.
The authors are grateful to all staff at the study centre who contributed to this study. This study was supported by Ministry of Major Science & Technology of Henan (201302005).
The authors declare that they have no conflicts of interest.
GJZ and GQZ, and PL: conceived of the study, and participated in its design and coordination and helped to draft the manuscript. PL, JL, RY, FRZ, HQW, HYC, YL, and SZD: collected the samples. PL, JL, YYW, WQZ, YWD and XNC: carried out part of experiments and wrote the manuscript. YYW, GJZ and GQZ performed the statistical analysis. All authors read and approved the final manuscript.