Background Human papillomavirus (HPV) E6 and E7 oncoproteins play a crucial role in HPV-related diseases, such as cervical cancer, and can be used as ideal targets for therapeutic vaccines. Human leukocyte antigen (HLA) participates in the immune response to block HPV infection and invasion by its target/recognition function. HPV-33 and HPV-58 are highly prevalent among Chinese women. Therefore, it is of great significance to study the E6 and E7 region-specific gene polymorphisms of HPV-33 and HPV-58 in Southwest China and to identify ideal epitopes for vaccine design. Both HPV-33 and HPV-58 belong to α-9 genus HPV and are highly homologous, so their correlations are included in our research. Methods To study the E6 and E7 variations and polymorphisms of HPV-33 and HPV-58 in Southwest China, we collected samples, extracted and sequenced DNA, and identified variants. Nucleotide sequences were translated into amino acids by Mega 6.0 software. The physical/chemical properties, amino acid-conserved sequences and secondary structure of protein sequences were analysed by the Protparam server, ConSurf server and PSIPRED software. The T and B cell epitopes of the E6/E7 reference and variant sequences in HPV-33 and HPV-58 were predicted by the Immune Epitope Database (IEDB) analysis server and the ABCpred server, respectively. Results Five and seven optimal HLA-I restricted T cell epitopes were selected from HPV-33 and HPV-58 E6, respectively, and these optimal epitopes are mainly located in 41-58 EVYDFAFADLTVVYREGN of HPV-33 E6 and 40-60 SEVYDFVFADLRIVYRDGNPF of HPV-58 E6. Six optimal HLA-I-restricted T cell epitopes were selected from HPV-33 and HPV-58 E7, and these epitopes are mainly located in 77-90 RTIQQLLMGTVNIV of HPV-33 E7 and 78-91 RTLQQLLMGTCTIV of HPV-58 E7. Conclusions HPV-33/HPV-58 E6/E7 gene polymorphisms and T/B cell epitopes of their reference and variant sequences were studied, and candidate epitopes were selected by bioinformatics techniques for therapeutic vaccine design for people in Southwest China. This study was the first to investigate the correlation of epitopes between HPV-33 and HPV-58. After experimental validation, these selected epitopes will be employed to induce a wide range of immune responses in heterogeneous HLA populations. Electronic supplementary material The online version of this article (10.1186/s12985-019-1168-y) contains supplementary material, which is available to authorized users.

The aim of this study was to determine the prevalence of human papillomavirus (HPV) infections in women with different gynaecological diseases (GDs) and ages in southwest China. With the assay of reverse dot blot, a total of 5544 cervical samples, including 1008 normal and 4536 abnormal specimens from women with GDs, were assayed for HPV detection. For the normal group, 20.2% (204/1008) tested positive for HPV, of which 40.2% (82/204) were high-risk (HR-HPV) genotypes. In the 4536 abnormal samples, 1569 women (34.6%) tested positive for HPV: 73.5% (1153/1569) had a single infection, 17.7% (277/1569) dual infection and 8.8% (138/1569) multiple HPV infection. Of 1569 HPV-positive specimens, 58.7% (921/1569) were infected with only HR-HPV genotypes. Significant differences in HPV infections were found among women of different ages (P < 0.01), number of pregnancies (P < 0.01), GDs (P < 0.01) and age at first sex (P < 0.01). In the present study, we found a high prevalence of HPV infection in women with GDs in southwest China. In addition to HPV types 16 and 18, a significant proportion of other HR-HPV genotypes were detected in this population.

BackgroundAfter the first epidemic wave from February through May 2013, the influenza A (H7N9) virus emerged and has followed a second epidemic wave since June 2013. As of June 27, 2014, the outbreak of H7N9 had caused 450 confirmed cases of human infection, with 165 deaths included. The case-fatality rate of all confirmed cases is about 36%, making the H7N9 virus a significant threat to people’s health. At present, neuraminidase inhibitors are the only licensed antiviral medications available to treat H7N9 infections in humans. Oseltamivir is the most commonly used inhibitor, and it is also a front-line drug for the threatening H7N9. Unfortunately, it has been reported that patients treated with oseltamivir can induce R294K (Arg294Lys) substitution in the H7N9 virus, which is a rare mutation and can reduce the antiviral efficacy of inhibitors. Even worse, deaths caused by such mutation after oseltamivir treatment have already been reported, indicating that the need to find substitutive neuraminidase inhibitors for currently available drugs to treat drug-resistant H7N9 is really pressing.Materials and methodsFirst, the structure of H7N9 containing the R294K substitution was downloaded from the Protein Data Bank, and structural information of approved drugs was downloaded from the ZINC (ZINC Is Not Commercial) database. Taking oseltamivir carboxylate as a reference drug, we then filtered these molecules through virtual screening to find out potential inhibitors targeting the mutated H7N9 virus. For further evaluation, we carried out a 14 ns molecular dynamic simulation for each H7N9–drug complex and calculated the binding energy for each candidate drug.ResultsWe found five inhibitors that could be candidate drugs for treating the mutated H7N9 virus. Docking poses showed these drugs could bind to the virus effectively, with the contribution of hydrogen bonds and hydrophobic interactions. With regard to the molecular dynamic simulations, receptor–ligand complexes formed by these candidate drugs were more stable than the one formed by oseltamivir carboxylate. The binding energy of oseltamivir carboxylate was −122.4 kJ/mol, while those for these potential inhibitors were −417.5, −404.7, −372.2, −304.3, and −289.9 kJ/mol, much better than the reference drug.ConclusionGiven the current and future threat of the mutated H7N9 virus, it is urgent that potent drugs and effective antiviral therapeutics be found. Our study therefore is able to complement currently available drugs for influenza A infectors and helps to prevent the ongoing threat of H7N9 virus.

Background Variations in human papillomavirus (HPV) E6 and E7 have been shown to be closely related to the persistence of the virus and the occurrence and development of cervical cancer. Long control region (LCR) of HPV has been shown multiple functions on regulating viral transcription. In recent years, there have been reports on E6/E7/LCR of HPV-16 and HPV-58, but there are few studies on HPV-52, especially for LCR. In this study, we focused on gene polymorphism of the HPV-52 E6/E7/LCR sequences, assessed the effects of variations on the immune recognition of viral E6 and E7 antigens, predicted the effect of LCR variations on transcription factor binding sites and provided more basic date for further study of E6/E7/LCR in Chengdu, China. Methods LCR/E6/E7 of the HPV-52 were amplified and sequenced to do polymorphic and phylogenetic analysis. Sequences were aligned with the reference sequence by MEGA 7.0 to identify SNP. A neighbor-joining phylogenetic tree was constructed by MEGA 7.0, followed by the secondary structure prediction of the related proteins using PSIPRED 4.0. The selection pressure of E6 and E7 coding regions were estimated by Bayes empirical Bayes analysis of PAML 4.9. The HLA class-I and II binding peptides were predicted by the Immune Epitope Database server. The B cell epitopes were predicted by ABCpred server. Transcription factor binding sites in LCR were predicted by JASPAR database. Results 50 SNP sites (6 in E6, 10 in E7, 34 in LCR) were found. From the most variable to the least variable, the nucleotide variations were LCR > E7 > E6. Two deletions were found between the nucleotide sites 7387–7391 (TTATG) and 7698–7700 (CTT) in all samples. A deletion was found between the nucleotide sites 7287–7288 (TG) in 97.56% (40/41) of the samples. The combinations of all the SNP sites and deletions resulted in 12 unique sequences. As shown in the neighbor-joining phylogenetic tree, except for one belonging to sub-lineage C2, others sequences clustered into sub-lineage B2. No positive selection was observed in E6 and E7. 8 non-synonymous amino acid substitutions (including E3Q and K93R in the E6, and T37I, S52D, Y59D, H61Y, D64N and L99R in the E7) were potential affecting multiple putative epitopes for both CD4+ and CD8+ T-cells and B-cells. A7168G was the most variable site (100%) and the binding sites for transcription factor VAX1 in LCR. In addition, the prediction results showed that LCR had the high probability binding sites for transcription factors SOX9, FOS, RAX, HOXA5, VAX1 and SRY. Conclusion This study provides basic data for understanding the relation among E6/E7/LCR mutations, lineages and carcinogenesis. Furthermore, it provides an insight into the intrinsic geographical relatedness and biological differences of the HPV-52 variants, and contributes to further research on the HPV-52 therapeutic vaccine development.

Artocarpus nanchuanensis (Moraceae), which is naturally distributed in China, is a representative and extremely endangered tree species. In this study, we obtained a high-quality chromosome-scale genome assembly and annotation information for A. nanchuanensis using integrated approaches, including Illumina, Nanopore sequencing platform, and Hi-C. A total of 128.71 Gb of raw Nanopore reads were generated from 20-kb libraries, and 123.38 Gb of clean reads were obtained after filtration with 160.34× coverage depth and a 17.48-kb average read length. The final assembled A. nanchuanensis genome was 769.44 Mb with a 2.09 Mb contig N50, and 99.62% (766.50 Mb) of the assembled data was assigned to 28 pseudochromosomes. In total, 39,596 genes (95.10%, 39,596/41,636) were successfully annotated, and 129 metabolic pathways were detected. Plants disease resistance/insect resistance genes, plant–pathogen interaction metabolic pathways, and abundant biosynthesis pathways of vitamins, flavonoid, and gingerol were detected. Unigene reveals the basis of species-specific functions, and gene family in contraction and expansion generally implies strong functional differences in the evolution. Compared with other related species, a total of 512 unigenes, 309 gene families in contraction, and 559 gene families in expansion were detected in A. nanchuanensis. This A. nanchuanensis genome information provides an important resource to expand our understanding of the unique biological processes, nutritional and medicinal benefits, and evolutionary relationship of this species. The study of gene function and metabolic pathway in A. nanchuanensis may reveal the theoretical basis of a special trait in A. nanchuanensis and promote the study and utilization of its rare medicinal value.