- Open Access
E6 and E7 gene polymorphisms in human papillomavirus Type-6 identified in Southwest China
Virology Journal volume 16, Article number: 114 (2019)
Human papillomavirus type-6 (HPV6) is the major etiological agent of anogenital warts both men and women. The present study aimed to characterize the genetic diversity among HPV6 in Southwest China, and to investigate the origin of, selective pressure experienced by, and impact of the resultantly identified genetic variants on the HPV6 secondary structure.
Phylogenetic trees were constructed by Maximum-likelihood and the Kimura 2-parameters methods by Molecular Evolutionary Genetics Analysis version 6.0. The diversity of secondary structure was analyzed by PSIPred software. The selection pressures acting on the E6/E7 genes were estimated by Phylogenetic Analyses by Maximum Likelihood version 4.8 software.
HPV6 was the most prevalent low risk HPV type in southwest China. In total, 143 E6 and E7 gene sequences of HPV6 isolated from patients were sequenced and compared to GenBank HPV6 reference sequence X00203. The results of these analyses revealed that both the HPV6 E6 and E7 were highly conserved within the analyzed patient samples, and comprised only 3 types of variant sequence, respectively. Furthermore, the analysis of HPV6 E6 and E7 sequences revealed seven/five single-nucleotide mutations, two/four and five/one of which were non-synonymous and synonymous, respectively. The phylogenetic analyses of the E6 and E7 sequences indicated that they belonged to sub-lineage A1 and sub-lineage B1, whereas the selective pressure analyses showed that only the E7 mutation sites 4R, 34E, and 52F were positive selection.
HPV6 (detection rate = 13.10%) was very prevalent in southwest China, both the HPV6 E6 and E7 sequences were highly conserved within the analyzed patient samples in southwest China, indicating that the low risk HPV6 can adapt to the environment well without much evolution.
Human papillomaviruses (HPVs) are associated with a variety of epithelial lesions, including benign genital warts and cervical intraepithelial neoplasia . To date, more than 250 HPV types have been identified and each of these genotypes are associated with infection at particular anatomical sites. HPV6 may be the most prevalent low risk alpha-papillomavirus type and is commonly associated with genital warts . For example, anogenital warts are primarily caused by HPV6 (family Papillomaviridae, genus Alphapapillomavirus, species 10) , which brings a significant burden to both the healthcare system and patients. Similarly, one third of Dutch primary school children have cutaneous warts, of which approximately 20% seek medical treatment each year . Generally speaking, these infections are classified as “not carcinogenic” or “low risk”, they often attract negative attention, and thereby cause significant psychological distress . However, some HPV6 variants are classified as “carcinogenic”, because they cause infections that lead to potentially fatal conditions, such as tonsillar, and malignant laryngeal carcinoma and/or malignant laryngeal papilloma [6,7,8,9].
To date, extensive research has been conducted to investigate sequence variation among carcinogenic HPV types; nevertheless, only limited data is available regarding HPV6 variants, despite their significant impact on human health. Structurally, the HPV is a double-stranded, circular DNA virus that encodes E1, E2, E4, E5, E6, E7, L1 and L2 proteins . HPVs infect cells via the basal layer of the stratified epithelium, and viral gene expression is closely linked to the endogenous differentiation program of the host cells . Of the HPV-encoded proteins, E6 and E7 have been shown to be the most important pathogenic HPV proteins. They have been previously shown to function as oncoproteins that critically regulate HPV-induced tumorigenesis . Furthermore, they have also been demonstrated to be essential to maintain the extrachromosomal forms of HPV in undifferentiated basal cells .
Genetic variability analyses have proven essential to facilitate an improved understanding of the evolution of the papillomavirus. A number of carcinogenic variants have been identified in HPV variants isolated from populations in Southwest China; however, only limited research has been conducted to identify low risk HPV variants. Thus, the present study aimed to analyze E6 and E7 sequence variability among HPV6 isolated from cervical papilloma samples collected from patients in Southwest China. Phylogenetic analyses were conducted to compare the identified nucleotide sequences with those previously described in other ethnic populations. In addition, the secondary structure of the identified sequences were predicted to assess the probably impact of the low risk variants on overall viral function. The results of the study could provide important data for the research on HPV6 prevention, diagnostic, therapeutic and even the design of therapeutic vaccines based on proteins E6 and E7 in Southwest China.
Methods & Materials
Clinical samples and HPV typing
From May 8, 2013 to June 1, 2016, cervical swabs were obtained from patients (with informed consent and ethical approval) at the Affiliated Hospital of Zunyi Medical University, Angel Women’s and Children’s Hospital, Sichuan Reproductive Health Research Centre Affiliated Hospital, and the Chengdu Western Hospital Maternity Unit. Women over 18 years old with visible cervical lesions and/or HPV-related diseases (e.g. cervical papilloma) were eligible for inclusion. Specimens were stored at − 20 °C until DNA extraction and HPV typing. Specimens DNA were extracted and tested using the Human Papillomavirus Genotyping Kit For 23 Types (PCR-RDB, reverse Dot Blot) according to the manufacturer’s instructions (Yaneng Bio, Shenzhen, China). This kit enabled the classification of the 23 HPV types (HPV 16, 18, 31, 33, 35, 39, 45, 51, 52, 53, 56, 58, 59, 66, 68, 73, 83, MM4, 6, 11, 42, 43, and 44).
In total, 216 samples were subjected to PCR amplification with E6 and E7 gene primers (see Additional file 1) that were designed using Primer Premier 5.0 software (Premier Biosoft, California, USA) and the HPV6 reference sequences (X00203) listed in the GenBank database (https://www.ncbi.nlm.nih.gov/genbank/). The reaction mixture comprised 2.5 mM dNTPs, 2 U Taq DNA polymerase, 90 pmol of each primer, and sufficient 10× Mg2+ PCR buffer to reach a final volume of 25 μL. Amplification of the sample sequences was achieved using thermal cycling conditions comprising an initial denaturation step (5 min at 95 °C), followed by 39 amplification cycles (each consisting of a 45 s denaturation step at 94 °C, a 60 s annealing step at 53 °C, and a 60 s elongation step at 72 °C). All amplified HPV6 E6 and E7 DNA products were detected via electrophoresis using an agarose gel that was impregnated with ethidium bromide, and sequenced using Sanger sequencing (Sango Biotech, Shanghai, China), and all the data were confirmed by repeating the PCR amplification and sequence analysis at least twice.
The secondary structures of the identified sequences were predicted with PSIPRED online software v.3.2 (http://bioinf.cs.ucl.ac.uk/psipred/), using default parameters. This software uses a very stringent cross-validation technique to achieve an average Q3 score of 81.6% . Any mutation observed to occur with a frequency of ≥10% was considered to be a major mutation. The sequences and variants were subsequently analyzed with various other software programs, including NCBI BLAST (https://blast.ncbi.nlm.nih.gov/Blast.cgi), Primer Premier 5 (Premier Biosoft), and DNAMAN version 5.2.2 (Lynnon Biosoft, California, USA). HPV6 nucleotide positions were numbered in accordance with the GenBank HPV6 reference sequence (X00203).
Phylogenetic analysis of the identified HPV6 variants
“Maximum-likelihood” phylogenetic trees were constructed for the identified E6 and E7 variants using MEGA v.6 software (http://www.megasoftware.net/home) to apply Kimura’s two-parameter model. Tree topology was evaluated by bootstrap resampling 1000 times . The reference viral sequences used to construct the distinct phylogenetic branches were collected from the GenBank database. Numbers above the branches indicate the bootstrap values that are greater than 60%.
Selective pressure analysis
Sites within the E6 and E7 HPV6 gene sequences that were likely to be subject to positive selection were predicted by using PAML 4.8 software (codeml tool; http://abacus.gene.ucl.ac.uk/software/paml.html) to perform likelihood ratio tests (LRTs) to infer nonsynonymous/synonymous nucleotide divergence, according to the method described by Nei and Gojobor [16,17,18]. The E6 and E7 protein sequences were aligned using ClustalX v. 2.1 software (ftp://ftp.ebi.ac.uk/pub/software/clustalw2/) .
The sequence results were compared to reference sequences in NCBI (https://www.ncbi.nlm.nih.gov/nuccore/). Seven reference sequences were obtained (KU298876; HE599232; FR751328; FR751337; AF092932; L41216; JN252318).
Characteristics of HPV6 prevalence in Southwest China
Nine thousand, three hundred and forty nine specimens (mean age 31.56 ± 8.42 years old) were collected. The overall positive rate of detectable 23 HPV types was 33.33% (3116/9349). 22.01% (2058/9349) specimens were high risk HPV types and 11.32% (1058/9349) were low risk HPV types. In all HPV-positive women involved in the study, the most common detected HPV types were HPV16 (detection rate = 13.12%) and HPV6 (detection rate = 13.10%) (Fig. 1). HPV6 was the most prevalent low risk HPV, followed by HPV11 (265, 8.50%), HPV-43 (234, 7.50%) and HPV42 (152, 4.87%). No HPV59 and HPV44 types were detected in the patients tested.
E6 sequence variations
In total, 143 HPV6 E6 gene sequences were obtained from the analyzed patient specimens using Sanger sequencing. A comparison of these with the GenBank HPV6 reference sequence (X00203) identified seven polymorphic sites, two of which comprised non-synonymous (D14Y, H50Q), and five of which comprised synonymous mutations. One of the identified non-synonymous mutations, D14Y, occurred within a region that encodes the alpha helix, whereas three of the synonymous mutations occurred in regions encoding the α-helix or β-strand. Notably, the mutation A120T and G372A were discovered in all specimens (Table 1). The table is designed with the reference to previous study [20, 21]. No insertion or deletion mutations, nor any evidence of introduced premature stop codons, were detected within the analyzed E6 HPV6 sequences. The conducted alignment of the analyzed (453 bp) E6 sequences with reference sequences belonged to sub-lineage A1 (70.63%) and sub-lineage B1 (29.37%) (Fig. 2).
E7 sequence variations
Compared with the HPV6 reference sequence (X00203), five single-nucleotide changes were identified within the 297 bp E7 open reading frame (ORF), of which four substitutions were non-synonymous, and one substitutions were synonymous. The most common non-synonymous mutation, T155A, occurred with a frequency of 34.27%, resulted in the amino acid (AA) change of F52Y (Phe to Tyr). The C294A was the most common synonymous mutation, and occurred with a frequency of 34.26%. One non-synonymous mutation (F52Y) occurred in the E7 sequence encoding β-strand (Table 2). No insertions, deletions, or premature stop codons were identified within the analyzed E7 variants. The (297 bp) E7 variants mainly belonged to sublineages A1 (65.73%) and sublineages B1 (31.47%) (Fig. 3). HPV6E701 was newly found.
None of the identified variant sequences were predicted to introduce changes to the reference HPV6 E6 and E7 secondary structures (see Additional files 2, 3, 4 and 5). Overall, the secondary structure of the E6 protein was predicted to consist of 30.0% helix, 16.7% strand, and 53.3% coil elements, whereas that of the E7 protein was predicted to comprise 12.3% helix, 18.4% strand, and 69.3% coil elements. There was no obvious change in the predicted HPV6 E6 and E7 secondary structures.
The non-synonymous/synonymous rate ratio (dN/dS), can measures selective pressure. When one non-synonymous mutation has advantage in fitness, it will rise to a higher rate than synonymous mutation, resulting in dN/dS > 1. The variable dN/dS ratios were tested among various lineages using the PAML 4.8 software. This program uses the Nei and Gojobori method of estimating synonymous substitutions, which is an unweighted pathway method. The first step in the procedure requires enumeration of the number of synonymous and non-synonymous sites present at each codon, where each site may be both partially synonymous and non-synonymous. The next stage is to determine the number of synonymous and non-synonymous changes between each pair of aligned sequences, codon-by-codon. The results of the conducted selective pressure (i.e, variable dN/dS ratio) analysis of the identified HPV6 E6 and E7 gene sequences are summarized in Tables 3 and 4. No positively selected sites were identified within the E6 sequence; however, in contrast, the E7 variant sites R4E (arginine to glutamic acid), E34K (glutamic acid to lysine), and F52Y (phenylalanine to tyrosine), were predicted to be subjected to positive selection.
The analysis of sequence results and reference sequence
Comparation between the results of our sequence with the reference sequence has been listed in the Table 5. For E6, the mutations A120T, C150G, A264T, and C291T were discovered in all sequence listed. The mutation G40 T was only found in HPV6vc (Table 5). Point mutations altered the E6 AA sequence of non-prototypic isolates, which contained specific mutations at nucleotide positions 150 (H50Q). AA exchanges identified at nucleotide positions 40 (D14Y) were additional characteristic features. The mutation at nucleotide positions 55 (T19S) was not found in this study. For E7, A10G, G11A, and G100A were proved as new mutations in this study. The mutations C156T and A262G were not discovered in the research. All sequence listed contain the mutation C294A. Point mutations altered the E7 AA sequence of non-prototypic isolates, which contained specific mutations at nucleotide positions 155 (F52K). AA exchanges identified at nucleotide positions 10, 11 (R4E) and 100(E34K) were additional characteristic features. The mutation at nucleotide positions 262 (D88N) was not found in this study.
E6 and E7 are essential HPV E-gene products, they target P53 and retinoblastoma (Rb) tumuor-suppressor proteins, respectively. The degradation of PRb can initiate abnormal cell replication; the inhibition of p53 can cause abnormal replication cell to lose control. The HPV genome can only be replicated along with the replication of the host genome. Therefore, E6 and E7 are especially important in HPV lifecycle. The present study sequenced the E6 and E7 ORFs in 143 HPV6 isolates from patients in Southwest China, to assess the HPV6 genetic diversity and evolution characteristics within this population, and help further identify specific HPV6 variants.
HPV6 was the second prevalent HPV type and the most prevalent low risk HPV type in southwest China, indicating that HPV6 is highly adaptable to the environment in southwest China compare with most HPV types that are not so prevalent.
Seven E6 sequence variants were identified, two of which induced corresponding E6 amino acid changes. One of these, H50Q (His to Gln), was identified in 30.08% of the analyzed HPV6 sequences, the amino acid at position 50 of prototype HPV6 E6 may be located centrally between the two internal Zn-binding motifs, which is important for the E6 protein stability . These variants have not been discovered previously in Southwest China, and their functional impact requires further analysis. Overall, the HPV6 E6 sequence was shown to be highly conserved within the analyzed patient samples, which exhibited only two main genotypes, HPV6E601 and HPV6E603, that represented 67.13% (96/143) and 29.37% (42/143) of the samples, respectively. Likewise, five E7 variants were identified, four of which exhibited a corresponding E7 amino-acid change. The most common E7 nucleotide mutations were T155A and C294A. Notably, T155A mutation affected 34.27% (49/143) patients, and caused a F52Y (Phe to Tyr) amino acid conversion that was predicted impact the protein’s secondary structure. Furthermore, the E7 sequence also was proved to be highly conserved within the patient specimens, comprising only two main genotypes HPV6E702 and HPV6E703 that were exhibited by 31.47% (45/143) and 65.73% (94/143) of the analyzed patients, respectively, and the sequence of HPV6E703 is in conformity with the reference (X00203). The fact that both sequences were so highly conserved supports the hypothesis that they play vital roles in HPV6 structure and function [23, 24], consistent with previous studies by Dartmann . Thus, they are likely also promising targets for HPV6 primer design and diagnostic detection. Compared to deadly cancerogenic high risk HPV, low risk HPV is able to get along better with the host; low risk HPV is more likely to be ignored in clinical treatment and prevention. HPV6 E6 and E7 were much more conserved than main high risk HPV types (like 16, 33, 53 and 58), may indicating low risk HPV has lower evolutionary pressure in gene level [20, 21, 26].
Structurally, some of the detected mutations affected amino acids at critical positions related to known biological functions. In E6, one (I97) and three (D14Y, H50Q, and I74) mutations were found to affect the β-strand and α-helix-encoding regions, respectively, that are critical for structural stability. In E7, only one mutation (F52Y) was detected may affect the β-strand-encoding region. The nucleotides as positions 58, 61, 91, and 94 in the E7 protein sequence have been previously shown to act as zinc binding sites , but no mutations were identified by the present study to occur in these positions, nor in the consensus LXCXE RB1-binding site (positions 22–26), or zinc-binding motif.
Notably, the present study is first to conduct an analysis to assess whether the E6 and E7 sequences are subject to positive selection in Southwest China. The main characteristic of positive selection is that it causes an unusually rapid rise in allele frequency, thereby enabling a species to adapt rapidly to environmental changes . The selective pressure analysis showed that all the sites that evolved under positive selection were common non-synonymous mutations, indicating that the positively selected variations beneficial for HPV6 to accommodate their environments are wide-spread. Moreover, remarkably, the positive sites R4E, E34K, and F52Y were observed in HPV6 E7, they may have evolutionary significance in making HPV6 adaptive to their environments.
Previous studies revealed the existence of two variant lineages (lineage A, lineage B) and five variant sub-lineages (sub-lineage B1, B2, B3, B4, B5)  among HPV6 variants. The conducted phylogenetic analyses showed that 97.56% of the E6 and E7 sequences within the patient specimens belonged to sub-lineage A1 (Reference, HPV6b) and sub-lineage B1 (HE599232, LP243) . Lineage A has been shown to predominate in Asia; in contrast, lineage B is distributed globally [30, 31], and consists mainly of variants from sub-lineage B1. Given its global distribution, lineage B1 may represent the oldest HPV6 sub-lineage, and likely disseminated during early human evolution. This may have enabled it to migrate to different regions of the world prior to the emergence of other HPV6 sub-lineages , and for E7, a new sub-lineage was found, however, only three samples were detected, and therefore, it requires further study.
Certain strengths of this study were that all clinical samples were collected from the Southwest region of China which had a strong regional representation and several new mutations were discovered, which will provide real and valid data for development of therapeutic vaccines for affected people in the Southwest China. Thus, the results of the present study significantly expand the current knowledge of HPV6 genetic diversity in Southwest China, and also provide a valuable resource for future studies of HPV6 epidemiology, evolution, function, pathogenesis, and use as a therapeutic target. However, the cellular level research is required, and a series of cellular experiments about the mutations should be designed for future studies.
HPV6 was very prevalent in southwest China, both the HPV6 E6 and E7 sequences were highly conserved within the analysed patient samples in southwest China, may indicating that the low risk HPV6 can adapt to the environment well without much evolution.
Availability of data and materials
All data generated or analyzed during this study are included in this published article and GeneBank.
Bayes empirical Bayes
Molecular evolutionary genetics analysis version
Open reading frame
Phylogenetic analyses by maximun likelihood
Doorbar J, Griffin H. Intrabody strategies for the treatment of human papillomavirus-associated disease [J]. Expert Opin Biol Ther. 2007;7(5):677–89.
Cheah PL, Looi LM. Biology and pathological associations of the human papillomaviruses: a review [J]. Malays J Pathol. 1998;20(1):1–10.
Singh M, Thakral D, Kar HK, et al. Distinct clinico-immunological profile of patients infected with human papilloma virus genotypes 6 and 11[J]. Virus Dis. 2017;28(2):200–4.
De Koning MNC, Quint KD, Bruggink SC, et al. High prevalence of cutaneous warts in elementary school children and the ubiquitous presence of wart-associated human papillomavirus on clinically normal skin [J]. Br J Dermatol. 2015;172(1):196–201.
Stone A, Potton A. Emotional responses to disfigured faces and disgust sensitivity: an eye-tracking study [J]. J Health Psychol. 2017;1359105317692856.
DiLorenzo TP, Tamsen A, Abramson AL, et al. Human papillomavirus type 6a DNA in the lung carcinoma of a patient with recurrent laryngeal papillomatosis is characterized by a partial duplication [J]. J Gen Virol. 1992;73(2):423–8.
Bercovich JA, Centeno CR, Aguilar OG, et al. Presence and integration of human papillomavirus type 6 in a tonsillar carcinoma [J]. J Gen Virol. 1991;72(10):2569–72.
Zarod AP, Rutherford JD, Corbitt G. Malignant progression of laryngeal papilloma associated with human papilloma virus type 6 (HPV-6) DNA [J]. J Clin Pathol. 1988;41(3):280–3.
Venuti A, Manni V, Morello R, et al. Physical state and expression of human papillomavirus in laryngeal carcinoma and surrounding normal mucosa [J]. J Med Virol. 2000;60(4):396–402.
Garcã-A-Vallvã S, Alonso A, Bravo IG. Papillomaviruses: different genes have different histories [J]. Trends Microbiol. 2005;13(11):514–21.
Mcmurray HR, Nguyen D, Westbrook TF, et al. Biology of human papillomaviruses.[J]. Int J Exp Pathol. 2010;82(1):15–33.
Raiol WT, Amorim PSD, et al. Genetic variability and phylogeny of the high-risk HPV-31, −33, −35, −52, and −58 in Central Brazil [J]. J Med Virol. 2010;81(4):685–92.
Thomas JT, Hubert WG, Ruesch MN, et al. Human papillomavirus type 31 oncoproteins E6 and E7 are required for the maintenance of episomes during the viral life cycle in normal human keratinocytes [J]. Proc Natl Acad Sci. 1999;96(15):8449–54.
Buchan DWA, Minneci F, Nugent TCO, et al. Scalable web services for the PSIPRED protein analysis workbench [J]. Nucleic Acids Res. 2013;41(W1):W349–57.
Tamura K, Stecher G, Peterson D, et al. MEGA6: molecular evolutionary genetics analysis version 6.0[J]. Mol Biol Evol. 2013;30(12):2725–9.
Nei M, Gojobori T. Simple methods for estimating the numbers of synonymous and nonsynonymous nucleotide substitutions [J]. Mol Biol Evol. 1986;3(5):418–26.
Yang Z. PAML 4: phylogenetic analysis by maximum likelihood [J]. Mol Biol Evol. 2007;24(8):1586–91.
Hamza AA, Robene-Soustrade I, Jouen E, et al. MultiLocus sequence analysis- and amplified fragment length polymorphism-based characterization of xanthomonads associated with bacterial spot of tomato and pepper and their relatedness to Xanthomonas species [J]. Syst Appl Microbiol. 2012;35(3):183–90.
Hwang T. Detection and typing of human papillomavirus DNA by PCR using consensus primers in various cervical lesions of Korean women [J]. J Korean Med Sci. 1999;14(6).
Chen Z, Jing Y, Wen Q, et al. E6 and E7 gene polymorphisms in human papillomavirus Types-58 and 33 identified in Southwest China [J]. PLoS One. 2017;12(1):e0171140.
Chen Z, Jing Y, Wen Q, et al. L1andL2gene polymorphisms in HPV-58 and HPV-33: implications for vaccine design and diagnosis [J]. Virol J. 2016;13(1):167.
Barbosa MS, Lowy DR, Schiller JT. Papillomavirus polypeptides E6 and E7 are zinc-binding proteins [J]. J Virol. 1989;63(3):1404–7.
Donne AJ, Hampson L, Homer JJ, et al. The role of HPV type in recurrent respiratory papillomatosis [J]. Int J Pediatr Otorhinolaryngol. 2010;74(1):0–14.
Hummel M, Hudson JB, Laimins LA. Differentiation-induced and constitutive transcription of human papillomavirus type 31b in cell lines containing viral episomes [J]. J Virol. 1992;66(10):6070–80.
Dartmann K, Schwarz E, Gissmann L, et al. The nucleotide sequence and genome organization of human papilloma virus type 11[J]. Virology. 1986;151(1):124–30.
Cao M, Chenzhang Y, Ding X, et al. Genetic variability and lineage phylogeny of human papillomavirus Type-16 and -53 based on the E6, E7, and L1 genes in Southwest China.[J]. Gene. 2016;592(1):49–59.
Chemes LB, Camporeale G, Sánchez IE, et al. Cysteine-rich positions outside the structural zinc motif of human papillomavirus E7 provide conformational modulation and suggest functional redox roles [J]. Biochemistry. 2014;53(10):1680–96.
Iyengar VK, Reeve HK, Eisner T. Paternal inheritance of a female moth\"s mating preference [J]. Nature, 2002, 419(6909):830–832.
Jelen M M, Chen Z, Kocjan B J, et al. Global genomic diversity of human papillomavirus 6 based on 724 isolates and 190 complete genome sequences[J]. Journal of virology. 2014;88(13):7307–7316.
Danielewski JA, Garland SM, Jenny MC, et al. Human papillomavirus type 6 and 11 genetic variants found in 71 Oral and Anogenital epithelial samples from Australia [J]. PLoS One. 2013;8(5):e63892.
Kocjan BJ, Poljak M, Cimerman M, et al. Prevaccination genomic diversity of human papillomavirus genotype 6 (HPV 6)[J]. Virology. 2009;391(2):274–83.
We thank for Affiliated Hospital of Zunyi Medical University, the Angel Women’s and Children’s Hospital, Sichuan Reproductive Health Research Center Affiliated Hospital, and Chengdu Western Hospital Maternity Unit Hospital, and some additional hospitals.
This reseach was funded by Affiliated Hospital of Zunyi Medical University (number 201808).
Ethics approval and consent to participate
The study was approved by the Ethics Committee of Affiliated Hospital of Zunyi Medical University, the approval number was ZYFYLS2018(81). Before sample collection, a written informed consent was obtained from all the patients or their guardians, and patient/study subject privacy was carefully protected.
Consent for publication
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.