- Primary research
- Open Access
Genetic variations in E6, E7 and the long control region of human papillomavirus type 16 among patients with cervical lesions in Xinjiang, China
Cancer Cell Internationalvolume 19, Article number: 65 (2019)
Xinjiang is one of the areas with the highest incidence of cervical cancer in China. Genetic variation in Human papillomavirus type 16 (HPV16) may increase the ability of the virus to mediate carcinogenesis and immune escape, which are risk factors for the progression of cervical cancer. We investigated polymorphism in HPV16 and the distribution of its sub-lineages in the region by analyzing the E6, E7 and long control region (LCR) gene sequences from women with HPV16-positive cervical samples in Xinjiang.
A total of 138 cases of cervical lesions and squamous cell carcinoma with infection of HPV16 virus were collected. The E6 and E7 genes and LCR of HPV16 virus were sequenced and compared with the HPV16 European prototype reference and other HPV16 mutants for single nucleotide polymorphisms. Neighbor-joining phylogenetic trees were constructed using E6, E7 and LCR sequences.
Fourteen missense mutations were found in the E6 gene; the loci with the highest mutation frequency were T350G (36/75, 48%) and T178G (19/75, 25.3%). In the E7 gene, the locus with the highest mutation frequency was A647G (18/75, 24%). A total of 33 polymorphic sites were found in the LCR, of which T7447C (39/95, 40.1%) was the most frequent.
HPV16 in Xinjiang is mainly of the European variant, followed by the Asian variant type; no Africa 1, 2 or Asia–America variant types were found.
Cervical cancer is the second most common malignant disease after breast cancer among women globally and has been ranked fourth among cancer-related deaths. The newly diagnosed cases in China (130,000) account for about a quarter of those in the world (approximately half a million). Cervical cancer is the most common malignant tumor of the female reproductive tract and has a major impact on women’s health and life expectancy . Epidemiological data have confirmed that persistent infection with high-risk human papillomavirus (HR-HPV) (such as HPV16 and HPV18) can lead to atypical hyperplasia and cancer in the cervix, thus confirming HPV as a central cause of this malignant tumor . Xinjiang is one of the areas with the highest incidence of cervical cancer in China, most of which is closely related to HPV16 infection . The production of cervical cancer by HPV is mainly mediated by the E6 and E7 genes . In addition, the E1 protein promotes viral genome replication and the E2 protein is negatively correlated with oncogene expression. The integration and mutation of E1 and E2 may promote the expression of virus genes E6 and E7, thus leading to the occurrence and development of cervical cancer . E7 is an oncoprotein with high carcinogenic risk containing 98 amino acids, which makes epithelial cells immortal with the cooperation of HPV E6 protein .
Many researchers have sequenced the HPV16 genes isolated from cervical lesions. It has been found that the variation in HPV16 E6 and E7 is correlated with the progression of cervical lesions and that HPV16 variant and HPV16 E6 and E7 mutations vary by racial and geographic area [7,8,9].
A meta-analysis was reported in 2013 by Comet et al. , using the scope of cancer caused by HPV16 all over the world, which showed that the European (Eur) variant is the most common sub-lineage and may be associated with persistent infection by HPV and increased risk of the progression of cervical lesions. By contrast, the results of Sun et al.  in the Liaoning province of China show that the Eur variants are negatively correlated with the severity of cervical disease at a level of cervical intraepithelial neoplasia 2/3 (CIN 2/3) or higher. In that study, Asian (As) variants were found in precancerous lesions and cases of cervical cancer, and were associated with disease progression. In addition, epidemiological studies have suggested that Asia–America (AA) variants can promote continuous infection with HPV and disease progression . Smith  and Berumen , in Costa Rica and Mexico, have reported that both AA and North America (NA) variants can increase the risk of CIN 3 cervical cancer in women. At the same time, they opined that variation can increase the ability of keratinocytes to achieve in vitro transformation. Overall, different HPV variants are associated with severity of cervical lesion. The variations of the HPV16 gene sequence may increase the virulence of the virus and mediate immune escape, which are risk factors for cervical cancer progression. However, different studies have inconsistent results, which may be related to differences in the genetic background of local populations.
Our study was designed to analyze the E6, E7, and long control region (LCR) gene sequences of HPV16-positive cases of cervical disease in the female population of Xinjiang, to investigate polymorphism of the HPV16 sub-lineages and their distribution in the local area. We explored the relationships between HPV16 E6, E7, and LCR gene mutations and the incidence of cervical cancer in Xinjiang, to accumulate molecular epidemiological data for the study of HPV16 gene polymorphism in cervical lesions.
Materials and methods
Primers for HPV16 E6 were used to screen samples from women with cervical lesions or those who had undergone surgery for cervical squamous cell carcinoma, and biopsy samples. The participants attended the Friendship Hospital (in 2016), or the People’s Hospital of Kashi (southern Xinjiang, China) and the People’s Hospital of Autonomous region (northern Xinjiang, China) during 2011–2014. Cases were confirmed by pathological examination and identified as cervical squamous cell carcinoma. A total of 138 cases of cervical lesions and squamous cell carcinoma with HPV16 infection were obtained. None of the patients had a history of long-term residence in other places. The samples were stored at − 80 °C for further processing. Informed consent was obtained from all patients and the study protocol was reviewed and approved by the ethics committees of the hospitals.
DNA extraction and HPV16 identification
DNA extraction was carried out using an SK1252 genomic DNA extraction kit (Shanghai Sangon Biological Engineering Technology and Services Company, Shangai, China) according to the instructions of the kit. The HPV16 E6 primers were designed and polymerase chain reaction (PCR) was used to detect HPV16 virus. The primer sequences for identification of HPV16 E6 were: HPV16 E6-F: 5′-gacccagaaagttaccacag-3′, HPV16 E6-R: 5′-cacaacggtttgttgtattg-3′ (F, forward primer; R, reverse primer).
PCR amplification and sequencing
The HPV16 E6, E7 and LCR gene fragments were amplified by PCR. Each 50-µL PCR reaction mixture contained 20 pmol of each primer, 50 mM KCl, 2.5 mM MgCl2, 100 mM Tris–HCl, pH 8.3, 0.1% Triton X-100, 50 µM of each dNTP, 1.8 U of HotStar Taq polymerase (QIAGEN, Hilden, Germany) and 5 µL template DNA. The PCR reaction conditions were as follows: 94 °C for 5 min; 30 cycles of 55 °C for 45 s, 72 °C for 60 s, 94 °C for 15 s; 55 °C for 45 s, 72 °C for 5 min. The sequencing primers are shown in Table 1.
Phylogenetic analysis of HPV16 variants
The PCR products were purified using SAP (Promega) and Exo I (Epicentre) and sequenced directly using an ABI Big-Dye Terminator v3.1 Cycle Sequencing Kit on a DNA analyzer (ABI3130XL) at Genesky Biotechnologies Inc (Shanghai, China). Single nucleotide polymorphisms (SNP) were analyzed using Polyphred and were aligned with the European prototype (GenBank: NC_001526.2)  and other typical HPV16 variants including As (GenBank: AF534061; AB889492), Af-1 (GenBank: AF472508; HQ644238), AA (GenBank: AF402678) and AA1 (GenBank: HQ644247) for comparison. A phylogenetic tree was built using MEGA 6  (Fig. 1).
Mutation analysis of HPV16 E6, E7 and LCR genes
In total, 138 HPV16 positive samples were sequenced to give full-length HPV16 LCR sequences, and 90 cases were chosen to carry out full-length HPV16 E6 and E7 gene sequencing. Finally, 75 E6 and E7 sequences and 95 LCR sequences were obtained. The HPV16 prototype (European prototype, GenBank accession: NC_001526.2) was used as the standard for comparison, and the gene polymorphism sites are shown in Tables 2 and 3.
Nineteen single nucleotide changes were identified among the sequencing results for E6 and E7, including 14 missense mutations and five synonymous mutations: 11/14 missense mutations were distributed in the E6 gene (C168G, G176A, T178G, T178A, T183G, A275G, A276G, T295G, T350G, C360A and A442C) and 3/14 missense mutations (A645C, A646C, A647G) were located in the E7 gene. These caused the amino acid to change from threonine to serine (T22S), aspartic acid to asparagine (D25N), aspartic acid to glutamic acid (D25E), aspartic acid to glutamine (D25Q), isoleucine to arginine (I27R), asparagine to aspartic acid (N58D), asparagine to serine (N58S), aspartic acid to glutamic acid (D64E), leucine to valine (L83V), threonine to lysine (T86K), glutamic acid to aspartic acid (E113D), leucine to phenylalanine (L28F), asparagine to histidine (N29H) and asparagine to serine (N29S), respectively. The two synonymous mutations of the E6 gene were nt109 (T to C) and nt 131 (A to C). The most common mutation sites in HPV16 E6 were T350G (36/75, 48%) and T178G (19/75, 25.3%). Compared with E6, the E7 gene is more conserved; the most common mutation site in E7 was A647G (18/75, 24%), and only one sample of A645C and A646C showed mutation. Notably, 18 cases of the A647G mutation were combined with the T178G mutation, and we considered that the A647G of E7 and T178G of E6 were co-variations.
Among the results of LCR sequencing, 33 polymorphic sites were found (the Table lists only 29 sites), of which four polymorphic sites (G7191T, nt7432–7433:GC to CGG, G7518A and A7861 deletions) showed 100% mutation in 95 samples. The remaining 29 mutations were mainly T7447C (39/95, 40.1%), followed by T7199C (16/95, 16.8%), C7268T (16/95, 16.8%), A7285C (16/95, 16.8%), T7711G (16/95, 16.8%), A7727C (16/95, 16.8%), G7839A (16/95, 16.8%) and C24T (16/95, 16.8%). We also found that six loci, T7199C, C7268T, A7285C, A7727C, G7839A, and C24T, in the LCR sequence were possibly co-variations.
Phylogenetic tree analysis
Phylogenetic analysis used the neighbor-joining method and the Kimura 2-Parameter model with bootstrap resampling (1000 replicates) of Mega 6 to construct the phylogenetic tree. Bootstrap values > 95% indicate very high confidence; 75–94%: high confidence; 50–74%: low confidence; and < 50%: no confidence. The bootstrap values (≥ 50%) are indicated in internal nodes.
The results of the phylogenetic analysis of E6 and E7 (Fig. 1) showed that, among the 75 samples, 18 samples (bootstrap values = 77%) belonged to As and the remaining 57 samples (bootstrap values = 99%) belonged to Eur strains.
The LCR phylogenetic analysis result (Fig. 2) showed that, among the 95 samples, only 16 of the samples (bootstrap values = 86%) belonged to As; Af-1, Af-2 and AA strains were not found.
Distribution of HPV16 sub-lineages in different histopathological grades
According to the E6 and E7 sequences, the HPV16 in our study could be divided into the Ep, E-350G and As sub-lineages. Use of the Pearson Chi squared test or Fisher’s exact test on the distribution of the three sub-lineages in different histological grades (considering P < 0.05 as statistically significant), showed that in different histological grades, the distributions of the three sub-lineages were significantly different when comparing normal with CIN 1 and CIN 2/3 (P = 0.001), or CIN 2/3 with cervical cancer (P = 0.023) (Table 4). The analysis showed no significant difference in the distribution of the sub-lineages of the HPV16 LCR sequence (P > 0.05) (Table 5).
A large number of experiments have been conducted to study the mutation of a variety of high-risk human papillomaviruses, concentrating mainly on HPV16, 18, 33, 35, and 45. However, HPV16 is the primary pathogenic factor in the occurrence of cervical cancer, and therefore most reports are related to HPV16 variation. They indicate that mutation of the virus may lead to the substitution of amino acids in the corresponding encoded proteins, which can alter the biological characteristics and immunogenicity of the virus, and affect the carcinogenic ability of HPV16 . HPV16 gene variation shows a degree of regionality and it is divided it into six main branches according to geographic location: Eur, As, AA, Af-1, Af-2 and NA [10, 18]. In addition, Eur viruses have been divided into a few small branches, such as E-T350 (Ep), E-G350 (T350G), E-G131 (A131G) and E-C109 (T109C). Each study of genetic variation in HPV16 has shown that a particular gene locus mutation will allow the virus to escape more easily from the monitoring of the host immune system or increase the opportunity of secondary virus infection, or cause cells to be malignant.
In this study, a total of 13 polymorphic sites were identified in 75 HPV16 E6 sequences, and T350G (36/75, 48%) and T178G (19/75, 25.3%) comprised the majority of 11 missense mutations, this result is consistent with that of Cai et al. . Studies have shown that HPV16 E6 T350G (L83V) variants are prevalent in high cervical lesions in Moroccan women and are closely related to the progression of cervical cancer . T178G (D25E) variations are mainly distributed in the Asian population (such as China, Japan and South Korea) [21, 22]; the mutation can interact with Human Lymphocyte Antigen (HLA) gene polymorphisms and promote the development of cervical cancer. Chansaenroj et al.  also suggest this mutation can increase the likelihood of persistent viral infection and cervical cancer progression.
The results of this study show that the common polymorphism site of the E7 gene is A647G, which is similar to the results of Yang et al. . Interactions between the E7 gene and tumor suppressor protein Rb are widely considered to be one of the main causes of cervical cancer. Amino acids 21–34 of the E7 protein form a region that combines with tumor suppressor protein Rb , and the A647G mutation may block the physiological function of Rb, thereby maintaining long-term infection with HPV. In addition, most of the mutations in A647G occur in the As sub-lineage. Further studies are needed to determine whether the mutation leads to the higher carcinogenicity of the As sub-lineage. It was found that A645C (L28F) has an incidence of 19% in cervical cancer tissues from Korean women  and that this mutation also occurs in Italy and Japan [27, 28]. However, only one mutation was found in our study, further emphasizing that HPV mutation has regional characteristics.
Most studies aim only to investigate whether there is a combination of mutations within individual genes, and there are relatively few reports of multiple genes with joint mutations. However, it has been reported that joint mutation of E6 T350G and E7 A647G may be Chinese specific . In contrast to the above report, 36 cases of T350G polymorphism were found in this study, and there was no polymorphism of the joint site mentioned above. Among 19 samples with polymorphism in E6 T178G or E7 A647G, 18 cases (94.7%) had two loci changed at the same time, which was consistent with the studies of Ding . In addition, three cases showed the E7 T846C mutation in association with the combined mutation of the above two loci. The effect of these joint mutations on the carcinogenicity of HPV remains to be further studied.
The LCR is the most variable region of the HPV16 genome, and may exert a vital function in persistent virus infection and the progression of cervical cancer. It contains the sequences associated with transcriptional regulation, and it is also the replication origin of HPV [8, 30]. Accounting for 100% of all infections, G7191T and G7518A mutations, were predicted by Xi et al.  to be, respectively, the binding sites for FOXA1 (Forkhead box protein A1), which is involved in the regulation of breast cancer, liver cancer, prostate cancer, lung cancer and endometrial cancer, and for SOX9 (sex-determining region Y-box 9), which is the potential cervical cancer tumor suppressor that, through the activation of p21WAF1/CIP1, inhibits cervical cancer cell growth. The frequency of A7727C mutation is relatively high. It may be the binding site of transcription factor PHOX2A (paired-like homeobox 2a), which is involved in cell proliferation and lung cancer metastasis.
The phylogenetic tree of the E6 and E7 sequences constructed using MEGA 6 shows that the prevalent strains of HPV16 in the Xinjiang area were mainly of the Eur (57/75, 76%), followed by the As (18/75, 24%); no Af-1 or Af-2 and AA variant types were found. The same conclusion was obtained in the analysis of the LCR sequences. Some of these branches had lower bootstrap values and may be associated with mutations in certain sites. As discussed above, the Eur variant of HPV16 is the most common sub-lineage in Xinjiang. In this study, we found that the association of HPV16 sub-lineages with different histopathological grades was statistically significant (P = 0.007). Future studies will be conducted to explore the role of HPV16 variants in the development of cervical cancer, based on a larger sample size, and will involve in vitro cell culture.
The HPV vaccine has become a research hotspot, including therapeutic vaccination for individuals infected with HPV virus and those with genital tract disease caused by HPV. Compared with prophylactic vaccines, therapeutic vaccines are mainly targeted against antigenic determinants of HPV early proteins, E6 or E7. E6 and E7 are the main transforming genes: cervical cancer cells cannot evade the immune response because of the absence of antigens, so E6 and E7 proteins can be used as targets for therapeutic vaccines against cervical cancer. However, HPV virus vaccine is highly specific; there are differences in immunogenicity among different HPV16 types, thus increasing the difficulty in the development of therapeutic vaccines. Therefore, the investigation and analysis of different HPV16 variants in the population of specific regions has significance in uncovering the carcinogenic mechanism of HPV16 and in developing preventive and therapeutic vaccines against HPV.
Samples of HPV16 in Xinjiang were mainly of the European variant, followed by the Asian variant; no African 1, 2 or Asia–America variant types were found.
Wu LX, Yun LI, Mu LI, Li FX, Cui MY, Yang Y, et al. The expression and clinical significance of Sox 11 in cervical carcinogenesis. J Xi’an Jiaotong Univ (Medical Sciences). 2018;39(1):73–7.
Iancu IV, Anton G, Botezatu A, Huica I, Nastase A, Socolov DG, et al. LINC01101 and LINC00277 expression levels as novel factors in HPV-induced cervical neoplasia. J Cell Mol Med. 2017;21(12):3787–94.
Song D, Shi Q, Hou X, Zhenghai MA. The effect of mutations in the upstream regulatory region of HPV-16 on the activity of virus early promoter. China Oncol. 2017;27(2):109–14.
Jackson R, Rosa BA, Lameiras S, Cuninghame S, Bernard J, Floriano WB, et al. Functional variants of human papillomavirus type 16 demonstrate host genome integration and transcriptional alterations corresponding to their unique cancer epidemiology. BMC Genom. 2016;17:851.
Zhao JW, Fang F, Guo Y, et al. HPV16 integration probably contributes to cervical oncogenesis through interrupting tumor suppressor genes and inducing chromosome instability. J Exp Clin Cancer Res. 2016;35(1):180.
Tzenov YR, Andrews PG, Voisey K, Popadiuk P, Xiong J, Popadiuk C, et al. Human papilloma virus (HPV) E7-mediated attenuation of retinoblastoma (Rb) induces hPygopus2 expression via Elf-1 in cervical cancer. Mol Cancer Res. 2013;11(1):19–30.
Lopera EA, Baena A, Florez V, Montiel J, Duque C, Ramirez T, et al. Unexpected inverse correlation between Native American ancestry and Asian American variants of HPVl6 in admixed colombian cervical cancer cases. Infect Genet Evol. 2014;28:339–48.
Mosmann JP, Monetti MS, Frutos MC, Kiguen AX, Venezuela RF, Cuffini CG. Mutation detection of E6 and LCR genes from HPV 16 associated with carcinogenesis. Asian Pac J Cancer Prev. 2015;16(3):1151–7.
Yang Y, Ren J, Zhang Q. Distribution of human papilloma virus type 16 E6/E7 gene mutation in cervical precancer or cancer: a case control study in Guizhou Province, China. J Med Virol. 2016;88(2):345–50.
Cornet I, Gheit T, Franceschi S, Vignat J, Burk RD, Sylla BS, et al. Human papillomavirus type 16 genetic variants: phylogeny and classification based on E6 and LCR. J Virol. 2012;86(12):6855–61.
Sun Z, Lu Z, Liu J, Wang G, Zhou W, Yang L, et al. Genetic variations of E6 and long control region of human papillomavirus type16 from patients with cervical lesion in Liaoning, China. BMC Cancer. 2013;13:459.
Gheit T, Cornet I, Clifford GM, Iftner T, Munk C, Tommasino M, et al. Risks for persistence and progression by human papillomavirus type 16 variant lineages among a population-based sample of Danish women. Cancer Epidemiol Prev Biomark. 2011;20(7):1315–21.
Smith B, Chen Z, Reimers L, van Doorslaer K, Schiffman M, DeSalle R, et al. Sequence imputation of HPV16 genomes for genetic association studies. PLoS ONE. 2011;6(6):e21375.
Berumen J, Ordoñez RM, Lazcano E, Salmeron J, Galvan SC, Estrada RA, et al. Asian–American variants of human papillomavirus 16 and risk for cervical cancer: a case-control study. J Natl Cancer Inst. 2001;93(17):1325–30.
Seedorf K, Krammer G, Durst M, Suhai S, Rowekamp WG. Human papillomavirus type 16 DNA sequence. Virology. 1985;145:181–5.
Tamura K, Stecher G, Peterson D, Filipski A, Kumar S. MEGA6: molecular evolutionary genetics analysis version 6.0. Mol Biol Evol. 2013;30:2725–9.
Ren J, Zhang Q, Yang Y, et al. Association between D32E/M28V + L94P variant of Human papilloma virus type 16 E6/E7 gene with CIN and the cervical cancer. Chin J Exp Clin Virol. 2015;29(5):417–20.
Cornet I, Gheit T, Iannacone MR, Vignat J, Sylla BS, Del Mistro A, et al. HPV16 genetic variation and the development of cervical cancer worldwide. Br J Cancer. 2013;108(1):240–4.
Cai HB, Chen CC, Ding XH. Human papillomavirus type 16 E6 gene variations in Chinese population. Eur J Surg Oncol. 2010;36(2):160–3.
Qmichou Z, Khyatti M, Berraho M, Ennaji MM, Benbacer L, Nejjari C, Benjaafar N, Benider A, Attaleb M, El Mzibri M. Analysis of mutations in the E6 oncogene of human papillomavirus 16 in cervical cancer isolates from Moroccan women. BMC Infect Dis. 2013;13:378.
Shang Q, Wang Y, Fang Y, Wei L, Chen S, Sun Y, et al. Human papillomavirus type 16 variant analysis of E6, E7 and L1 genes and long control region in cervical carcinomas in patients in northeast China. J Clin Microbiol. 2011;49(7):2656–63.
Mi J, Rhee JE, Jang DH, Kim SS. Gene expression profiles are altered in human papillomavirus-16 E6 D25E-expressing cell lines. Virol J. 2011;8(1):1–6.
Chansaenroj J, Theamboonlers A, Junyangdikul P, Swangvaree S, Karalak A, Poovorawan Y. Whole genome analysis of human papillomavirus type 16 multiple infection in cervical cancer patients. Asian Pacific J Cancer Prev. 2012;13(2):599–606.
Yang L, Yang H, Kun W, Shi X, Ma S, Sun Q. Prevalence of HPV and variation of HPV16/HPV18 E6/E7 genes in cervical cancer in women in South West China. J Med Virol. 2014;86(11):1926–36.
Jenison SA, Yu XP, Valentine JM, Galloway DA. Characterization of human antibody-reactive epitopes encoded by human papillomavirus types 16 and 18. J Virol. 1991;65(3):1208.
Youk EG, Ku JL, Park JG. Detection and typing of human papillomavirus in anal epidermoid carcinomas: sequence variation in the E7 gene of human papillomavirus type 16. Dis Colon Rectum. 2001;44(2):236–42.
Fujinaga Y, Okazawa K, Nishikawa A, Yamakawa Y, Fukushima M, Kato I, et al. Sequence variation of human papillomavirus type 16 E7 in preinvasive and invasive cervical neoplasias. Virus Genes. 1994;9(1):85–92.
Tornesello ML, Duraturo ML, Salatiello I, Buonaguro L, Losito S, Botti G, Stellato G, Greggi S, Piccoli R, Pilotti S, Stefanon B. Analysis of human papillomavirus type-16 variants in Italian women with cervical intraepithelial neoplasia and cervical cancer. J Med Virol. 2004;74(1):117–26.
Ding T, Wang X, Ye F, Cheng X, Lu W, Xie X. Distribution of human papillomavirus 16 E6/E7 variants in cervical cancer and intraepithelial neoplasia in Chinese women. Int J Gynecol Cancer. 2010;20(8):1391–8.
Anacker DC, Moody CA. Modulation of the DNA damage response during the life cycle of human papillomaviruses. Virus Res. 2017;231:41–9.
Xi J, Chen J, Xu M, Yang H, Luo J, Pan Y, et al. Genetic variability and functional implication of the long control region in HPV-16 variants in Southwest China. PLoS ONE. 2017;12(8):e0182388.
ZP designed the study. XZ and HX carried out the molecular genetic studies, participated in the gene analysis and wrote the manuscript. ZP and FJ carried out specimen collection and DNA extraction. CZ and SF carried out HPV16 identification. DC and YL performed the statistical analysis. WZ, HL, DL and RS participated in experimental design and revised the manuscript. All authors read and approved the final manuscript.
We acknowledge funding support from the National Natural Science Foundation of China (Grant Numbers U1503125, 81660462), the International Science and Technology Collaboration Projector of Xinjiang Production and Construction Corps (Grant Number 2013BC003), the Importance of Scientific Research and Innovation Project of Shihezi University (Grant Number gxjs2013-zdgg05), and the High Level Talent of Scientific Research for a Starting Project of Special Grant of Shihezi University (Grant Number RCZX201534).
The authors declare that they have no competing interests.
Availability of data and materials
The datasets generated and/or analyzed during the current study are not publicly available because (The data belong to our group) but are available from the corresponding author on reasonable request.
Consent for publication
Ethics approval and consent to participate
The National Natural Science Foundation of China (Grant Numbers U1503125, 81660462), the International Science and Technology Collaboration Projector of Xinjiang Production and Construction Corps (Grant Number 2013BC003), the Importance of Scientific Research and Innovation Project of Shihezi University (Grant Number gxjs2013-zdgg05), and the High Level Talent of Scientific Research for a Starting Project of Special Grant of Shihezi University (Grant Number RCZX201534).
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.