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

Multilocus Sequence Analysis and Detection of Copper Ion Resistance of Xanthomonas phaseoli pv. manihotis Causing Bacterial Blight in Cassava

Cassava (Manihot esculenta Crantz) is an important tropical tuber crop around the world. Cassava bacterial blight, caused by Xanthomonas phaseoli pv. manihotis, is a key disease that influences cassava production worldwide. Between 2008 and 2020, 50 X. phaseoli pv. manihotis strains were isolated from diseased plant samples or acquired from China, Uganda, Cambodia, Colombia, Malaysia, and Micronesia. Using multilocus sequence analysis, the genetic diversity of X. phaseoli pv. manihotis strains was evaluated. A neighbor-joining phylogenetic dendrogram was constructed based on partial sequences of five housekeeping genes (atpD-dnaK-gyrB-efp-rpoD). The strains clustered into three groups whose clusters were consistent with atpD and RpoD gene sequences. Group I contained 46 strains from China, Uganda, Cambodia, and Micronesia, and the other two groups were comprised of strains from Colombia and Malaysia, respectively. The resistance of all these strains to copper ion (Cu²⁺) was determined, the minimal inhibitory concentration was between 1.3 and 1.7 mM, and there was no significant difference between strains from different geographic region. During genome annotation of the X. phaseoli pv. manihotis strain CHN01, homologous gene clusters of copLAB and xmeRSA were identified. The predicted amino acid sequences of two gene clusters were highly homologous with the copper-resistant protein from Xanthomonas strains. CopLAB and xmeRSA were amplified from all these strains, suggesting that the regulation of copper resistance is associated with two distinct metabolic pathways. CopLAB and xmeRSA were highly conserved among strains from different geographic regions, possibly associated with other conserved function.

Polymorphism of E6 and E7 Genes in Human Papillomavirus Types 31 and 33 in Northeast China

Persistent infection with human papillomavirus (HPV) types 31 and 33 is an important causative factor for cervical cancer. The E6/E7 genes are key oncogenes involved in the immortalization and transformation of human epithelial cells. Genetic polymorphism may lead to differences in the virus' carcinogenic potential, the immune reaction of the host, and the potencies of vaccines. Few studies on HPV31/33 E6/E7 genetic polymorphism have been carried out. To study the genetic polymorphism of HPV31 and HPV33 E6/E7 genes in northeast China, these genes (HPV31 E6/E7, n = 151; HPV33 E6/E7, n = 136) were sequenced and compared to reference sequences (J04353.1, M12732.1) using BioEdit. Phylogenetic trees were constructed by the neighbor-joining method using MegaX. The diversity of the secondary structure was estimated using the PSIPred server. The positively selected sites were analyzed using PAML4.9. The major histocompatibility complex (MHC) class I and MHCII epitopes were predicted using the ProPred-I server and ProPredserver. B-cell epitopes were predicted using the ABCpred server. In the 151 HPV31E6 sequences, 25 (25/450) single-nucleotide mutations were found, 14 of which were synonymous mutations and 11 were nonsynonymous. In the 151 HPV31E7 sequences, 8 (8/297) nucleotide mutations were found, 3 of which were synonymous mutations and 5 were nonsynonymous. In the 136 HPV33E6 sequences, 17 (17/450) nucleotide mutations were observed, 7 of which were synonymous mutations and 10 were nonsynonymous. C14T/G (T5I/S) was a triallelic mutation. Finally, in the 136 HPV33E7 sequences, 9 (9/294) nucleotide mutations were observed, 3 of which were synonymous mutations and 6 were nonsynonymous. C134T/A (A45V/E) and C278G/A (T93S/N) were triallelic mutations. Lineage A was the most common lineage in both HPV31 and HPV33. In all of the sequences, we only identified one positively selected site, HPV33 E6 (K93N). Most nonsynonymous mutations were localized at sites belonging to MHC and/or B-cell predicted epitopes. Data obtained in this study should contribute to the development and application of detection probes, targeted drugs, and vaccines.

A Pan-Global Study of Bacterial Leaf Spot of Chilli Caused by Xanthomonas spp

Bacterial Leaf Spot (BLS) is a serious bacterial disease of chilli (Capsicum spp.) caused by at least four different Xanthomonas biotypes: X. euvesicatoria pv. euvesicatoria, X. euvesicatoria pv. perforans, X. hortorum pv. gardneri, and X. vesicatoria. Symptoms include black lesions and yellow halos on the leaves and fruits, resulting in reports of up to 66% losses due to unsalable and damaged fruits. BLS pathogens are widely distributed in tropical and subtropical regions. Xanthomonas is able to survive in seeds and crop residues for short periods, leading to the infections in subsequent crops. The pathogen can be detected using several techniques, but largely via a combination of traditional and molecular approaches. Conventional detection is based on microscopic and culture observations, while a suite of Polymerase Chain Reaction (PCR) and Loop-Mediated Isothermal Amplification (LAMP) assays are available. Management of BLS is challenging due to the broad genetic diversity of the pathogens, a lack of resilient host resistance, and poor efficacy of chemical control. Some biological control agents have been reported, including bacteriophage deployment. Incorporating stable host resistance is a critical component in ongoing integrated management for BLS. This paper reviews the current status of BLS of chilli, including its distribution, pathogen profiles, diagnostic options, disease management, and the pursuit of plant resistance.

Xanthomonas hortorum - beyond gardens: Current taxonomy, genomics, and virulence repertoires

TAXONOMY

Bacteria; Phylum Proteobacteria; Class Gammaproteobacteria; Order Lysobacterales (earlier synonym of Xanthomonadales); Family Lysobacteraceae (earlier synonym of Xanthomonadaceae); Genus Xanthomonas; Species X. hortorum; Pathovars: pv. carotae, pv. vitians, pv. hederae, pv. pelargonii, pv. taraxaci, pv. cynarae, and pv. gardneri.

HOST RANGE

Xanthomonas hortorum affects agricultural crops, and horticultural and wild plants. Tomato, carrot, artichoke, lettuce, pelargonium, ivy, and dandelion were originally described as the main natural hosts of the seven separate pathovars. Artificial inoculation experiments also revealed other hosts. The natural and experimental host ranges are expected to be broader than initially assumed. Additionally, several strains, yet to be assigned to a pathovar within X. hortorum, cause diseases on several other plant species such as peony, sweet wormwood, lavender, and oak-leaf hydrangea.

EPIDEMIOLOGY AND CONTROL

X. hortorum pathovars are mainly disseminated by infected seeds (e.g., X. hortorum pvs carotae and vitians) or cuttings (e.g., X. hortorum pv. pelargonii) and can be further dispersed by wind and rain, or mechanically transferred during planting and cultivation. Global trade of plants, seeds, and other propagating material constitutes a major pathway for their introduction and spread into new geographical areas. The propagules of some pathovars (e.g., X. horturum pv. pelargonii) are spread by insect vectors, while those of others can survive in crop residues and soils, and overwinter until the following growing season (e.g., X. hortorum pvs vitians and carotae). Control measures against X. hortorum pathovars are varied and include exclusion strategies (i.e., by using certification programmes and quarantine regulations) to multiple agricultural practices such as the application of phytosanitary products. Copper-based compounds against X. hortorum are used, but the emergence of copper-tolerant strains represents a major threat for their effective management. With the current lack of efficient chemical or biological disease management strategies, host resistance appears promising, but is not without challenges. The intrastrain genetic variability within the same pathovar poses a challenge for breeding cultivars with durable resistance.

USEFUL WEBSITES

https://gd.eppo.int/taxon/XANTGA, https://gd.eppo.int/taxon/XANTCR, https://gd.eppo.int/taxon/XANTPE, https://www.euroxanth.eu, http://www.xanthomonas.org, http://www.xanthomonas.org/dokuwiki.

Bacterial Spot of Tomato and Pepper in Africa: Diversity, Emergence of T5 Race, and Management

Bacterial spot disease was first reported from South Africa by Ethel M. Doidge in 1920. In the ensuing century after the initial discovery, the pathogen has gained global attention in plant pathology research, providing insights into host-pathogen interactions, pathogen evolution, and effector discovery, such as the first discovery of transcription activation-like effectors, among many others. Four distinct genetic groups, including Xanthomonas euvesicatoria (proposed name: X. euvesicatoria pv. euvesicatoria), Xanthomonas perforans (proposed name: X. euvesicatoria pv. perforans), Xanthomonas gardneri (proposed name: Xanthomonas hortorum pv. gardneri), and Xanthomonas vesicatoria, are known to cause bacterial spot disease. Recently, a new race of a bacterial spot pathogen, race T5, which is a product of recombination between at least two Xanthomonas species, was reported in Nigeria. In this review, our focus is on the progress made on the African continent, vis-à-vis progress made in the global bacterial spot research community to provide a body of information useful for researchers in understanding the diversity, evolutionary changes, and management of the disease in Africa.

A centenary for bacterial spot of tomato and pepper

DISEASE SYMPTOMS

Symptoms include water-soaked areas surrounded by chlorosis turning into necrotic spots on all aerial parts of plants. On tomato fruits, small, water-soaked, or slightly raised pale-green spots with greenish-white halos are formed, ultimately becoming dark brown and slightly sunken with a scabby or wart-like surface.

HOST RANGE

Main and economically important hosts include different types of tomatoes and peppers. Alternative solanaceous and nonsolanaceous hosts include Datura spp., Hyoscyamus spp., Lycium spp., Nicotiana rustica, Physalis spp., Solanum spp., Amaranthus lividus, Emilia fosbergii, Euphorbia heterophylla, Nicandra physaloides, Physalis pubescens, Sida glomerata, and Solanum americanum.

TAXONOMIC STATUS OF THE PATHOGEN

Domain, Bacteria; phylum, Proteobacteria; class, Gammaproteobacteria; order, Xanthomonadales; family, Xanthomonadaceae; genus, Xanthomonas; species, X. euvesicatoria, X. hortorum, X. vesicatoria.

SYNONYMS (NONPREFERRED SCIENTIFIC NAMES)

Bacterium exitiosum, Bacterium vesicatorium, Phytomonas exitiosa, Phytomonas vesicatoria, Pseudomonas exitiosa, Pseudomonas gardneri, Pseudomonas vesicatoria, Xanthomonas axonopodis pv. vesicatoria, Xanthomonas campestris pv. vesicatoria, Xanthomonas cynarae pv. gardneri, Xanthomonas gardneri, Xanthomonas perforans.

MICROBIOLOGICAL PROPERTIES

Colonies are gram-negative, oxidase-negative, and catalase-positive and have oxidative metabolism. Pale-yellow domed circular colonies of 1-2 mm in diameter grow on general culture media.

DISTRIBUTION

The bacteria are widespread in Africa, Brazil, Canada and the USA, Australia, eastern Europe, and south-east Asia. Occurrence in western Europe is restricted.

PHYTOSANITARY CATEGORIZATION

A2 no. 157, EU Annex designation II/A2.

EPPO CODES

XANTEU, XANTGA, XANTPF, XANTVE.

Trends in Molecular Diagnosis and Diversity Studies for Phytosanitary Regulated Xanthomonas

Bacteria in the genus Xanthomonas infect a wide range of crops and wild plants, with most species responsible for plant diseases that have a global economic and environmental impact on the seed, plant, and food trade. Infections by Xanthomonas spp. cause a wide variety of non-specific symptoms, making their identification difficult. The coexistence of phylogenetically close strains, but drastically different in their phenotype, poses an added challenge to diagnosis. Data on future climate change scenarios predict an increase in the severity of epidemics and a geographical expansion of pathogens, increasing pressure on plant health services. In this context, the effectiveness of integrated disease management strategies strongly depends on the availability of rapid, sensitive, and specific diagnostic methods. The accumulation of genomic information in recent years has facilitated the identification of new DNA markers, a cornerstone for the development of more sensitive and specific methods. Nevertheless, the challenges that the taxonomic complexity of this genus represents in terms of diagnosis together with the fact that within the same bacterial species, groups of strains may interact with distinct host species demonstrate that there is still a long way to go. In this review, we describe and discuss the current molecular-based methods for the diagnosis and detection of regulated Xanthomonas, taxonomic and diversity studies in Xanthomonas and genomic approaches for molecular diagnosis.

Molecular Epidemiology of Xanthomonas euvesicatoria Strains from the Balkan Peninsula Revealed by a New Multiple-Locus Variable-Number Tandem-Repeat Analysis Scheme

Bacterial spot of pepper and tomato is caused by at least three species of Xanthomonas, among them two pathovars of Xanthomonas euvesicatoria, which are responsible for significant yield losses on all continents. In order to trace back the spread of bacterial spot pathogens within and among countries, we developed the first multilocus variable number of tandem repeat analyses (MLVA) scheme for pepper- and tomato-pathogenic strains of X. euvesicatoria. In this work, we assessed the repeat numbers by DNA sequencing of 16 tandem repeat loci and applied this new tool to analyse a representative set of 88 X. euvesicatoria pepper strains from Bulgaria and North Macedonia. The MLVA-16 scheme resulted in a Hunter–Gaston Discriminatory Index (HGDI) score of 0.944 and allowed to resolve 36 MLVA haplotypes (MTs), thus demonstrating its suitability for high-resolution molecular typing. Strains from the different regions of Bulgaria and North Macedonia were found to be widespread in genetically distant clonal complexes or singletons. Sequence types of the variable number of tandem repeats (VNTR) amplicons revealed cases of size homoplasy and suggested the coexistence of different populations and different introduction events. The large geographical distribution of MTs and the existence of epidemiologically closely related strains in different regions and countries suggest long dispersal of strains on pepper in this area.

Clarifying the taxonomy of the causal agent of bacterial leaf spot of lettuce through a polyphasic approach reveals that Xanthomonas cynarae Trébaol et al. 2000 emend. Timilsina et al. 2019 is a later heterotypic synonym of Xanthomonas hortorum Vauterin et al. 1995

Assessment of the taxonomy and diversity of Xanthomonas strains causing bacterial leaf spot of lettuce (BLSL), commonly referred to as Xanthomonas campestris pv. vitians, has been a long-lasting issue which held back the global efforts made to understand this pathogen. In order to provide a sound basis essential to its study, we conducted a polyphasic approach on strains obtained through sampling campaigns or acquired from collections. Results of a multilocus sequence analysis crossed with phenotypic assays revealed that the pathotype strain does not match the description of the nomenspecies provided by Brown in 1918. However, strain LMG 938=CFBP 8686 does fit this description. Therefore, we propose that it replaces LMG 937=CFBP 2538 as pathotype strain of X. campestris pv. vitians. Then, whole-genome based phylogenies and overall genome relatedness indices calculated on taxonomically relevant strains exhibited the intermediate position of X. campestris pv. vitians between closely related species Xanthomonas hortorum and Xanthomonas cynarae. Phenotypic profiles characterized using Biolog microplates did not reveal stable diagnostic traits legitimizing their distinction. Therefore, we propose that X. cynarae Trébaol et al. 2000 emend. Timilsina et al. 2019 is a later heterotypic synonym of X. hortorum, to reclassify X. campestris pv. vitians as X. hortorum pv. vitians comb. nov. and to transfer X. cynarae pathovars in X. hortorum as X. hortorum pv. cynarae comb. nov. and X. hortorum pv. gardneri comb. nov. An emended description of X. hortorum is provided, making this extended species a promising model for the study of Xanthomonas quick adaptation to different hosts.

Phylogenetic Analyses of Xanthomonads Causing Bacterial Leaf Spot of Tomato and Pepper: Xanthomonas euvesicatoria Revealed Homologous Populations Despite Distant Geographical Distribution

Bacterial leaf spot of tomato and pepper (BLS), an economically important bacterial disease caused by four species of Xanthomonas (X. euvesicatoria (Xe), X. vesicatoria (Xv), X. gardneri (Xg), and X. perforans (Xp)), is a global problem and can cause over 50% crop loss under unfavorable conditions. Among the four species, Xe and Xv are prevalent worldwide. Characterization of the pathogens is crucial for disease management and regulatory purposes. In this study, we performed a multilocus sequence analysis (MLSA) with six genes (hrcN, dnaA gyrB, gapA, pdg, and hmbs) on BLS strains. Other Xanthomonas species were included to determine phylogenetic relationships within and among the tested strains. Four BLS species comprising 76 strains from different serological groups and diverse geographical locations were resolved into three major clades. BLS xanthomonads formed distinct clusters in the phylogenetic analyses. Three other xanthomonads, including X. albilineans, X. sacchari, and X. translucens pv. undolusa revealed less than 85%, 88%, and 89% average nucleotide identity (ANI), respectively, with the other species of Xanthomonas included in this study. Both antibody and MLSA data showed that Xv was clearly separated from Xe and that the latter strains were remarkably clonal, even though they originated from distant geographical locations. The Xe strains formed two separate phylogenetic groups; Xe group A1 consisted only of tomato strains, whereas Xe group A2 included strains from pepper and tomato. In contrast, the Xv group showed greater heterogeneity. Some Xv strains from South America were closely related to strains from California, while others grouped closer to a strain from Indiana and more distantly to a strain from Hawaii. Using this information molecular tests can now be devised to track distribution of clonal populations that may be introduced into new geographic areas through seeds and other infected plant materials.

Independent Evolution with the Gene Flux Originating from Multiple Xanthomonas Species Explains Genomic Heterogeneity in Xanthomonas perforans

Xanthomonas perforans is the predominant pathogen responsible for bacterial leaf spot of tomato and X. euvesicatoria for that of pepper in the southeast United States. Previous studies have indicated significant changes in the X. perforans population collected from Florida tomato fields over the span of 2 decades, including a shift in race and diversification into three phylogenetic groups driven by genome-wide homologous-recombination events derived from X. euvesicatoria In our sampling of Xanthomonas strains associated with bacterial spot disease in Alabama, we were readily able to isolate X. perforans from symptomatic pepper plants grown in several Alabama counties, indicating a recent shift in the host range of the pathogen. To investigate the diversity of these pepper-pathogenic strains and their relation to populations associated with tomatoes grown in the southeast United States, we sequenced the genomes of eight X. perforans strains isolated from tomatoes and peppers grown in Alabama and compared them with previously published genome data available from GenBank. Surprisingly, reconstruction of the X. perforans core genome revealed the presence of two novel genetic groups in Alabama that each harbored a different transcription activation-like effector (TALE). While one TALE, AvrHah1, was associated with an emergent lineage pathogenic to both tomato and pepper, the other was identified as a new class within the AvrBs3 family, here designated PthXp1, and was associated with enhanced symptom development on tomato. Examination of patterns of homologous recombination across the larger X. euvesicatoria species complex revealed a dynamic pattern of gene flow, with multiple donors of Xanthomonas spp. associated with diverse hosts of isolation.IMPORTANCE Bacterial leaf spot of tomato and pepper is an endemic plant disease with a global distribution. In this study, we investigated the evolutionary processes leading to the emergence of novel X. perforans lineages identified in Alabama. While one lineage was isolated from symptomatic tomato and pepper plants, confirming the host range expansion of X. perforans, the other lineage was isolated from tomato and acquired a novel transcription activation-like effector, here designated PthXp1. Functional analysis of PthXp1 indicated that it does not induce Bs4-mediated resistance in tomato and contributes to virulence, providing an adaptive advantage to strains on tomato. Our findings also show that different phylogenetic groups of the pathogen have experienced independent recombination events originating from multiple Xanthomonas species. This suggests a continuous gene flux between related xanthomonads associated with diverse plant hosts that results in the emergence of novel pathogen lineages and associated phenotypes, including host range.

E6 and E7 gene polymorphisms in human papillomavirus Type-6 identified in Southwest China

BACKGROUND

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.

METHODS

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.

RESULTS

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.

CONCLUSIONS

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.

AFLP protocol comparison for microbial diversity fingerprinting

Over the last decade, several methods based on genomic DNA have been developed for the identification and genotyping of prokaryotic and eukaryotic organisms. These genomic methods differ regarding taxonomic range, discriminatory power, reproducibility, and ease of interpretation and standardization. The amplified fragment length polymorphism (AFLP) technique is a very powerful DNA fingerprinting technique for DNA of any source or complexity, varying in both size and base composition. In addition, this method shows high discriminatory power and good reproducibility allowing it to be efficient in discriminating at both the species and strain levels. The development and application of AFLP have allowed significant progress in the study of biodiversity and taxonomy of microorganisms. In the last years, the Applied Biosystems AFLP Microbial Fingerprinting Kit, now out of production, was widely used in various studies to perform AFLP characterization of selected bacteria strains (described by Vos et al. (Nucleic Acids Res 23(21):4407-4414, 1995)). Its replacement gives the possibility for laboratories to continue the use of the previous AFLP data as a reference for bacteria genetic fingerprinting analysis in biodiversity studies. To overcome this issue a result comparison, by using an improved AFLP protocol and the AFLP commercial kit, was performed. In particular, previous results on different species (Listeria monocytogenes, Lactobacillus plantarum, and Streptococcus thermophilus) obtained with the commercial kit were compared with the improved AFLP procedure to validate the protocol. When compared with the AFLP Microbial Fingerprinting Kit, the improved protocol shows high reproducibility, resolution, and overall, is a faster method with lower costs.

A duplex quantitative real-time PCR assay for the detection and quantification of Xanthomonas phaseoli pv. dieffenbachiae from diseased and latently infected anthurium tissue

Anthurium bacterial blight caused by Xanthomonas phaseoli pv. dieffenbachiae (formerly Xanthomonas axonopodis pv. dieffenbachiae) is the major phytosanitary threat in many anthurium growing areas worldwide. Reliable and sensitive diagnostic tools are required for surveillance and certification programs. A duplex real-time quantitative PCR assay was developed for the detection and quantification of X. phaseoli pv. dieffenbachiae from anthurium tissue. This PCR assay targeted a X. phaseoli pv. dieffenbachiae-specific gene encoding an ABC transporter and an internal control encoding for chalcone synthase in Anthurium andreanum. A cycle threshold (Ct), using a receiver-operating characteristic approach (ROC), was implemented to ensure that the declaration of a positive sample was reliable. The duplex real-time assay displayed very high performance with regards to analytical specificity (100% inclusivity, 98.9% exclusivity), analytical sensitivity (LOD_95% = 894 bacteria/ml corresponding to 18 bacteria per reaction) and repeatability. We demonstrated the pertinence of this real-time quantitative PCR assay for detecting X. phaseoli pv. dieffenbachiae from diseased leaf tissue (collected from outbreaks on anthurium) and from asymptomatic, latently infected anthurium plants. This assay could be useful for surveillance, as well as for indexing propagative plant material for the presence of X. phaseoli pv. dieffenbachiae.

Genetic variability in E5, E6, E7 and L1 genes of human papillomavirus type 31

Human papillomavirus (HPV) type 31 is an important pathogenic subtype associated with cervical cancer. The aims of the present study were to analyze E5, E6, E7 and L1 gene mutations of HPV‑31 among females, and to elucidate the evolutionary associations between them. In total, 87 positive samples were collected. The E5, E6, E7 and L1 genes were amplified by polymerase chain reaction and sequenced. Subsequently, two phylogenetic trees were constructed from the nucleotide sequences of the E5, E6 and E7 and the L1 variants of HPV‑31. In total, 31 mutation sites of E5, E6 and E7 genes were identified, of which 16 were non‑synonymous. T4053A (F80I), C285T (H60Y), C520T (A138V) and A743G (K62E) were the most common non‑synonymous mutations. A total of 30 mutation sites of L1 genes were identified, of which four were non‑synonymous. The most common non‑synonymous mutations of L1 genes were A6350G (T29A) and C6372A (T36N). By phylogenetic analysis, A and C variants were most frequently detected, while B variants were less frequently detected in this population. The sequence variation data obtained in the present study provides a foundation for future research regarding HPV‑induced oncogenesis, and may prove valuable for developing diagnostic probes and in the design of HPV vaccines for targeted populations.

Genetic Diversity of the Ralstonia solanacearum Species Complex in the Southwest Indian Ocean Islands

Epidemiological surveillance of plant pathogens based on genotyping methods is mandatory to improve disease management strategies. In the Southwest Indian Ocean (SWIO) islands, bacterial wilt (BW) caused by the Ralstonia solanacearum species complex (RSSC) is hampering the production of many sustainable and cash crops. To thoroughly analyze the genetic diversity of the RSSC in the SWIO, we performed a wide sampling survey (in Comoros, Mauritius, Reunion, Rodrigues, and Seychelles) that yielded 1,704 isolates from 129 plots, mainly from solanaceous crops. Classification of the isolates to the four major RSSC phylogenetic groups, named phylotypes, showed that 87% were phylotype I, representing the most prevalent strain in each of the SWIO islands. Additionally, 9.7% were phylotype II, and 3.3% were phylotype III; however, these isolates were found only in Reunion. Phylotype IV (2 isolates), known to be restricted to Indonesia-Australia-Japan, was reported in Mauritius, representing the first report of this group in the SWIO. Partial endoglucanase (egl) sequencing, based on the selection of 145 isolates covering the geographic and host diversity in the SWIO (also including strains from Mayotte and Madagascar), revealed 14 sequevars with Reunion and Mauritius displaying the highest sequevar diversity. Through a multilocus sequence analysis (MLSA) scheme based on the partial sequencing of 6 housekeeping genes (gdhA, gyrB, rplB, leuS, adk, and mutS) and 1 virulence-associated gene (egl), we inferred the phylogenetic relationships between these 145 SWIO isolates and 90 worldwide RSSC reference strains. Phylotype I was the most recombinogenic, although recombination events were detected among all phylotypes. A multilocus sequence typing (MLST) scheme identified 29 sequence types (STs) with variable geographic distributions in the SWIO. The outstanding epidemiologic feature was STI-13 (sequevar I-31), which was overrepresented in the SWIO and obviously reflected a lineage strongly adapted to the SWIO environment. A goeBURST analysis identified eight clonal complexes (CCs) including SWIO isolates, four CCs being geographically restricted to the SWIO, and four CCs being widespread beyond the SWIO. This work, which highlights notable genetic links between African and SWIO strains, provides a basis for the epidemiological surveillance of RSSC and will contribute to BW management in the SWIO.

Pathogenic, phenotypic and molecular characterisation of Xanthomonas nasturtii sp. nov. and Xanthomonas floridensis sp. nov., new species of Xanthomonas associated with watercress production in Florida

We describe two new species of the genus Xanthomonas, represented by yellow mucoid bacterial strains isolated from diseased leaves of watercress (Nasturtium officinale) produced in Florida, USA. One strain was pathogenic on watercress, but not in other species including a range of brassicas; other strains were not pathogenic in any of the tested plants. Data from Biolog carbon source utilization tests and nucleotide sequence data from 16S and gyrB loci suggested that both pathogenic and non-pathogenic strains were related to, yet distinct from, previously described Xanthomonas species. Multilocus sequence analysis and whole genome-wide comparisons of the average nucleotide identity (ANI) of genomes of two strains from watercress showed that these are distinct and share less than 95 % ANI with all other known species; the non-pathogenic strain WHRI 8848 is close to Xanthomonascassavae (ANI of 93.72 %) whilst the pathogenic strain WHRI 8853 is close to a large clade of species that includes Xanthomonasvesicatoria (ANI ≤90.25 %). Based on these results, we propose that both strains represent new Xanthomonas species named Xanthomonas floridensis sp. nov. (type strain WHRI 8848=ATCC TSD-60=ICMP 21312=LMG 29665=NCPPB 4601) and Xanthomonas nasturtii sp. nov. (type strain WHRI 8853=ATCC TSD-61=ICMP 21313=LMG 29666=NCPPB 4600), respectively. The presence of non-pathogenic Xanthomonas strains in watercress and their interaction with pathogenic strains needs to be further investigated. Although the importance of the new pathogenic species is yet to be determined, the bacterial disease that it causes constitutes a threat to watercress production and its distribution should be monitored.

E6 and E7 Gene Polymorphisms in Human Papillomavirus Types-58 and 33 Identified in Southwest China

Cancer of the cervix is associated with infection by certain types of human papillomavirus (HPV). The gene variants differ in immune responses and oncogenic potential. The E6 and E7 proteins encoded by high-risk HPV play a key role in cellular transformation. HPV-33 and HPV-58 types are highly prevalent among Chinese women. To study the gene intratypic variations, polymorphisms and positive selections of HPV-33 and HPV-58 E6/E7 in southwest China, HPV-33 (E6, E7: n = 216) and HPV-58 (E6, E7: n = 405) E6 and E7 genes were sequenced and compared to others submitted to GenBank. Phylogenetic trees were constructed by Maximum-likelihood and the Kimura 2-parameters methods by MEGA 6 (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 PAML 4.8 (Phylogenetic Analyses by Maximun Likelihood version4.8) software. The positive sites of HPV-33 and HPV-58 E6/E7 were contrasted by ClustalX 2.1. Among 216 HPV-33 E6 sequences, 8 single nucleotide mutations were observed with 6/8 non-synonymous and 2/8 synonymous mutations. The 216 HPV-33 E7 sequences showed 3 single nucleotide mutations that were non-synonymous. The 405 HPV-58 E6 sequences revealed 8 single nucleotide mutations with 4/8 non-synonymous and 4/8 synonymous mutations. Among 405 HPV-58 E7 sequences, 13 single nucleotide mutations were observed with 10/13 non-synonymous mutations and 3/13 synonymous mutations. The selective pressure analysis showed that all HPV-33 and 4/6 HPV-58 E6/E7 major non-synonymous mutations were sites of positive selection. All variations were observed in sites belonging to major histocompatibility complex and/or B-cell predicted epitopes. K93N and R145 (I/N) were observed in both HPV-33 and HPV-58 E6.

L1 and L2 gene polymorphisms in HPV-58 and HPV-33: implications for vaccine design and diagnosis

Background

Cervical cancer is associated with infection by certain subtypes of human papillomavirus (HPV). The L1 protein comprising HPV vaccine formulations elicits high-titre neutralizing antibodies and confers protection against specific HPV subtypes. HPV L2 protein is an attractive candidate for cross-protective vaccines. HPV-33 and HPV-58 are very prevalent among Chinese women.

Methods

To study the gene intratypic variations and polymorphisms of HPV-33 and HPV-58 L1/L2 in Sichuan China, HPV-33 and HPV-58 L1 and L2 genes were sequenced and compared with other genes submitted to GenBank. Phylogenetic trees were constructed by maximum-likelihood and the Kimura 2-parameters methods (MEGA 6). The secondary structure was analyzed by PSIPred software, and HPV-33 and HPV-58 L1 homology models were created by SWISS-MODEL software. The selection pressures acting on the L1/L2 genes were estimated by PAML 4.8.

Results

Among 124 HPV-33 L1 sequences 20 single nucleotide mutations were observed included 8/20 non-synonymous and 12/20 synonymous mutations. The 101 HPV-33 L2 sequences included 12 single nucleotide mutations comprising 7/12 non-synonymous and 5/12 synonymous mutations. The 223 HPV-58 L1 sequences included 32 single nucleotide mutations comprising 9/32 non-synonymous and 23/32 synonymous mutations. The 201 HPV-58 L2 sequences comprised 26 single nucleotide mutations including 9/26 non-synonymous and 17/26 synonymous mutations. Selective pressure analysis showed that most of the common non-synonymous mutations showed a positive selection. HPV-33 and HPV-58 L2 were more stable than HPV-33 and HPV-58 L1.

Conclusions

HPV-33 and HPV-58 L2 were better candidates as clinical diagnostic targets compared with HPV-33 and HPV-58 L1. Clinical diagnostic probes and second-generation polyvalent vaccines should be designed on the basis of the unique sequence of HPV-33 and 58 L1/L2 variations in Sichuan, to improve the accuracy of clinical detection and the protective efficiency of vaccines.

Electronic supplementary material

The online version of this article (doi:10.1186/s12985-016-0629-9) contains supplementary material, which is available to authorized users.

Xanthomonas euvesicatoria Causes Bacterial Spot Disease on Pepper Plant in Korea

In 2004, bacterial spot-causing xanthomonads (BSX) were reclassified into 4 species-Xanthomonas euvesicatoria, X. vesicatoria, X. perforans, and X. gardneri. Bacterial spot disease on pepper plant in Korea is known to be caused by both X. axonopodis pv. vesicatoria and X. vesicatoria. Here, we reidentified the pathogen causing bacterial spots on pepper plant based on the new classification. Accordingly, 72 pathogenic isolates were obtained from the lesions on pepper plants at 42 different locations. All isolates were negative for pectolytic activity. Five isolates were positive for amylolytic activity. All of the Korean pepper isolates had a 32 kDa-protein unique to X. euvesicatoria and had the same band pattern of the rpoB gene as that of X. euvesicatoria and X. perforans as indicated by PCR-restriction fragment length polymorphism analysis. A phylogenetic tree of 16S rDNA sequences showed that all of the Korean pepper plant isolates fit into the same group as did all the reference strains of X. euvesicatoria and X. perforans. A phylogenetic tree of the nucleotide sequences of 3 housekeeping genes-gapA, gyrB, and lepA showed that all of the Korean pepper plant isolates fit into the same group as did all of the references strains of X. euvesicatoria. Based on the phenotypic and genotypic characteristics, we identified the pathogen as X. euvesicatoria. Neither X. vesicatoria, the known pathogen of pepper bacterial spot, nor X. perforans, the known pathogen of tomato plant, was isolated. Thus, we suggest that the pathogen causing bacterial spot disease of pepper plants in Korea is X. euvesicatoria.

Genetic variability in E6, E7, and L1 genes of human papillomavirus genotype 52 from Southwest China

Human papillomavirus (HPV) is the major causative agent of cervical cancer, which accounts for the second highest cancer burden in women worldwide. HPV-52, the prevalent subtype in Asia, especially in southwest China, was analyzed in this study. To analyze polymorphisms, intratypic variants, and genetic variability in the E6-E7 (n=26) and L1 (n=53) genes of HPV-52, these genes were sequenced and the sequences were submitted to GenBank. Phylogenetic trees were constructed using the neighbor-joining and Kimura 2-parameters methods, followed by analysis of the diversity of secondary structure. Finally, we estimated the selection pressures acting on the E6-E7 and L1 genes. Fifty-one novel variants of HPV-52 L1, and two novel variants of HPV-52 E6-E7 were identified in this study. Thirty single nucleotide changes were observed in HPV-52 E6-E7 sequences with 19/30 non-synonymous mutations and 11/30 synonymous mutations (five in the alpha helix and five in the beta sheet). Fifty-five single nucleotide changes were observed in HPV-52 L1 sequences with 17/55 non-synonymous mutations (seven in the alpha helix and fourteen in the beta sheet) and 38/55 synonymous mutations. Selective pressure analysis predicted that most of these mutations reflect positive selection. Identifying new variants in HPV-52 may inform the rational design of new vaccines specifically for women in southwest China. Knowledge of genetic variation in HPV may be useful as an epidemiologic correlate of cervical cancer risk, or may even provide critical information for developing diagnostic probes.

Bacterial spot of tomato and pepper: diverse Xanthomonas species with a wide variety of virulence factors posing a worldwide challenge

TAXONOMIC STATUS

Bacteria; Phylum Proteobacteria; Class Gammaproteobacteria; Order Xanthomonadales; Family Xanthomonadaceae; Genus Xanthomonas; Species Xanthomonas euvesicatoria, Xanthomonas vesicatoria, Xanthomonas perforans and Xanthomonas gardneri.

MICROBIOLOGICAL PROPERTIES

Gram-negative, rod-shaped bacterium, aerobic, motile, single polar flagellum.

HOST RANGE

Causes bacterial spot disease on plants belonging to the Solanaceae family, primarily tomato (Solanum lycopersicum), pepper (Capsicum annuum) and chilli peppers (Capsicum frutescens).

DISEASE SYMPTOMS

Necrotic lesions on all above-ground plant parts.

DISTRIBUTION

Worldwide distribution of X. euvesicatoria and X. vesicatoria on tomato and pepper; X. perforans and X. gardneri increasingly being isolated from the USA, Canada, South America, Africa and Europe. A wide diversity within the bacterial spot disease complex, with an ability to cause disease at different temperatures, makes this pathogen group a worldwide threat to tomato and pepper production. Recent advances in genome analyses have revealed the evolution of the pathogen with a plethora of novel virulence factors. Current management strategies rely on the use of various chemical control strategies and sanitary measures to minimize pathogen spread through contaminated seed. Chemical control strategies have been a challenge because of resistance by the pathogen. Breeding programmes have been successful in developing commercial lines with hypersensitive and quantitative resistance. However, durability of resistance has been elusive. Recently, a transgenic approach has resulted in the development of tomato genotypes with significant levels of resistance and improved yield that hold promise. In this article, we discuss the current taxonomic status, distribution of the four species, knowledge of virulence factors, detection methods and strategies for disease control with possible directions for future research.

Classification and Taxonomy of Vegetable Macergens

Macergens are bacteria capable of releasing pectic enzymes (pectolytic bacteria). These enzymatic actions result in the separation of plant tissues leading to total plant destruction. This can be attributed to soft rot diseases in vegetables. These macergens primarily belong to the genus Erwinia and to a range of opportunistic pathogens namely: the Xanthomonas spp., Pseudomonas spp., Clostridium spp., Cytophaga spp., and Bacillus spp. They consist of taxa that displayed considerable heterogeneity and intermingled with members of other genera belonging to the Enterobacteriaceae. They have been classified based on phenotypic, chemotaxonomic and genotypic which obviously not necessary in the taxonomy of all bacterial genera for defining bacterial species and describing new ones These taxonomic markers have been used traditionally as a simple technique for identification of bacterial isolates. The most important fields of taxonomy are supposed to be based on clear, reliable and worldwide applicable criteria. Hence, this review clarifies the taxonomy of the macergens to the species level and revealed that their taxonomy is beyond complete. For discovery of additional species, further research with the use modern molecular methods like phylogenomics need to be done. This can precisely define classification of macergens resulting in occasional, but significant changes in previous taxonomic schemes of these macergens.

Genetic diversity analysis of isolates belonging to the Photobacterium phosphoreum species group collected from salmon products using AFLP fingerprinting

An accurate amplified fragment length polymorphism (AFLP) method, including three primer sets for the selective amplification step, was developed to display the phylogenetic position of Photobacterium isolates collected from salmon products. This method was efficient for discriminating the three species Photobacterium phosphoreum, Photobacterium iliopiscarium and Photobacterium kishitanii, until now indistinctly gathered in the P. phosphoreum species group known to be strongly responsible for seafood spoilage. The AFLP fingerprints enabled the isolates to be separated into two main clusters that, according to the type strains, were assigned to the two species P. phosphoreum and P. iliopiscarium. P. kishitanii was not found in the collection. The accuracy of the method was validated by using gyrB-gene sequencing and luxA-gene PCR amplification, which confirmed the species delineation. Most of the isolates of each species were clonally distinct and even those that were isolated from the same source showed some diversity. Moreover, this AFLP method may be an excellent tool for genotyping isolates in bacterial communities and for clarifying our knowledge of the role of the different members of the Photobacterium species group in seafood spoilage.

Genome sequencing reveals a new lineage associated with lablab bean and genetic exchange between Xanthomonas axonopodis pv. phaseoli and Xanthomonas fuscans subsp. fuscans

Common bacterial blight is a devastating seed-borne disease of common beans that also occurs on other legume species including lablab and Lima beans. We sequenced and analyzed the genomes of 26 strains of Xanthomonas axonopodis pv. phaseoli and X. fuscans subsp. fuscans, the causative agents of this disease, collected over four decades and six continents. This revealed considerable genetic variation within both taxa, encompassing both single-nucleotide variants and differences in gene content, that could be exploited for tracking pathogen spread. The bacterial strain from Lima bean fell within the previously described Genetic Lineage 1, along with the pathovar type strain (NCPPB 3035). The strains from lablab represent a new, previously unknown genetic lineage closely related to strains of X. axonopodis pv. glycines. Finally, we identified more than 100 genes that appear to have been recently acquired by Xanthomonas axonopodis pv. phaseoli from X. fuscans subsp. fuscans.

Multilocus sequence analysis of xanthomonads causing bacterial spot of tomato and pepper plants reveals strains generated by recombination among species and recent global spread of Xanthomonas gardneri

Four Xanthomonas species are known to cause bacterial spot of tomato and pepper, but the global distribution and genetic diversity of these species are not well understood. A collection of bacterial spot-causing strains from the Americas, Africa, Southeast Asia, and New Zealand were characterized for genetic diversity and phylogenetic relationships using multilocus sequence analysis of six housekeeping genes. By examining strains from different continents, we found unexpected phylogeographic patterns, including the global distribution of a single multilocus haplotype of X. gardneri, possible regional differentiation in X. vesicatoria, and high species diversity on tomato in Africa. In addition, we found evidence of multiple recombination events between X. euvesicatoria and X. perforans. Our results indicate that there have been shifts in the species composition of bacterial spot pathogen populations due to the global spread of dominant genotypes and that recombination between species has generated genetic diversity in these populations.

Development of multilocus variable-number tandem repeat analysis (MLVA) for Xanthomonas arboricola pathovars

Xanthomonas arboricola is an important bacterial species, the pathovars of which are responsible for bacterial blight diseases on stone fruit, hazelnut, Persian walnut, poplar, strawberry, poinsettia and banana. In this study, we evaluated variable number tandem repeats (VNTR) as a molecular typing tool for assessing the genetic diversity within pathovars of X. arboricola. Screening of the X. arboricola pv. pruni genome sequence (CFBP5530 strain) predicted 51 candidate VNTR loci. Primer pairs for polymerase chain reaction (PCR) amplification of all 51 loci were designed, and their discriminatory power was initially evaluated with a core collection of 8 X. arboricola strains representative of the different pathovars. Next, the 26 polymorphic VNTR loci present in all strains were used for genotyping a collection of 61 strains. MLVA is a typing method that clearly differentiates X. arboricola strains. The MLVA scheme described in this study is a rapid and reliable molecular typing tool that can be used for further epidemiological studies of bacterial diseases caused by X. arboricola pathovars.

Role of blossom colonization in pepper seed infestation by Xanthomonas euvesicatoria

Colonization of Xanthomonas euvesicatoria was investigated in pepper blossoms and the relationship between inoculum concentrations and seed infestation was determined. Inoculation of blossoms resulted in asymptomatic pepper fruit. However, real-time polymerase chain reaction detected X. euvesicatoria in 39% of the seed lots assayed and viable colonies were recovered from 35% of them. Successful transmission occurred in 16% of the seed lots tested. In a separate experiment, X. euvesicatoria reached populations of up to 1 × 10(5) CFU/blossom on stigmas 96 h after inoculation. Bacteria colonized stylar and ovary tissues with populations ranging from 1 × 10(5) to 1 × 10(6) CFU/blossom 96 h after inoculation. A positive correlation existed between inoculum concentration and percentage of infested seedlots. Blossoms inoculated with Acidovorax citrulli also resulted in infested pepper seedlots. Furthermore, A. citrulli colonized pepper blossoms significantly better than X. euvesicatoria by 96 h postinoculation. It was concluded that pepper blossoms can be a potential site of ingress for X. euvesicatoria into seed, and blossom colonization may be involved in pepper seed infestation. Data also indicated that seed infestation via blossoms may be nonspecific because nonhost plants can be colonized by incompatible pathogens. Thus, host-pathogen interactions may not be important for bacterial ingress through blossoms.

Harnessing population genomics to understand how bacterial pathogens emerge, adapt to crop hosts, and disseminate

Crop diseases emerge without warning. In many cases, diseases cross borders, or even oceans, before plant pathologists have time to identify and characterize the causative agents. Genome sequencing, in combination with intensive sampling of pathogen populations and application of population genetic tools, is now providing the means to unravel how bacterial crop pathogens emerge from environmental reservoirs, how they evolve and adapt to crops, and what international and intercontinental routes they follow during dissemination. Here, we introduce the field of population genomics and review the population genomics research of bacterial plant pathogens over the past 10 years. We highlight the potential of population genomics for investigating plant pathogens, using examples of population genomics studies of human pathogens. We also describe the complementary nature of the fields of population genomics and molecular plant-microbe interactions and propose how to translate new insights into improved disease prevention and control.

Genetic variability of HPV-58 E6 and E7 genes in Southwest China

HPV accounts for most of incidence of cervical cancer. Genetic variations of E6 and E7 may be associated with the development of cervical cancer in specific geographic regions. HPV-58 has been found to be a relatively prevalent high-risk HPV among southwest Chinese women. To explore gene intratypic variations and polymorphisms of HPV-58 E6 and E7 genes originating in Southwest China, a total of 2000 scraped cell samples were collected for DNA extraction and HPV typing. Then, the E6 and E7 genes of HPV-58 (n=22) were sequenced and compared to others submitted to GenBank, followed by an analysis of the diversity of secondary structure by DNASTAR software. Phylogenetic trees were then constructed by Neighbor-Joining and the Kimura 2-parameters methods, followed by an analysis of selection pressures acting on the E6/E7 genes by PAML software. 22 were HPV-58 positive among 215 high-risk types' samples. The nucleotide variation rate of E6 was 86.36% (19/22) among the 22 HPV-58 E6 sequences studied. 4 single nucleotide changes were identified among the E6 sequences with 3/4 synonymous mutations (C187T, A260C, C307T) and 1/4 non-synonymous mutations (A388C, from Lys to Asn, in alpha helix). The most common mutations of E6 genes are the C307T and A388C. 8 single nucleotide changes were identified among the HPV-58 E7 sequences with 2/8 synonymous mutations (T726C, T744G) and 6/8 non-synonymous mutations (G599A, C632T, G694A, G760A, G761A, T803C). The nucleotide variation rate of E7 was 72.73% (17/22). The most common mutations of E7 genes are C632T, G694A, T744G, G760A (from Gly to Ser, in turn), G761A and T803C. The phylogenetic analyses demonstrate that all HPV-58 E6/E7 variants identified belonged to the Southeast Asia lineage. There was no evidence of positive selection in the sequence alignment of HPV-58 E6 and E7 genes.

Genetic variability in L1 and L2 genes of HPV-16 and HPV-58 in Southwest China

HPV account for most of the incidence of cervical cancer. Approximately 90% of anal cancers and a smaller subset (<50%) of other cancers (oropharyngeal, penile, vaginal, vulvar) are also attributed to HPV. The L1 protein comprising HPV vaccine formulations elicits high-titre neutralizing antibodies and confers type restricted protection. The L2 protein is a promising candidate for a broadly protective HPV vaccine. In our previous study, we found the most prevalent high-risk HPV infectious serotypes were HPV-16 and HPV-58 among women of Southwest China. To explore gene polymorphisms and intratypic variations of HPV-16 and HPV-58 L1/L2 genes originating in Southwest China, HPV-16 (L1: n = 31, L2: n = 28) and HPV-58 (L1: n = 21, L2: n = 21) L1/L2 genes were sequenced and compared to others described and submitted to GenBank. Phylogenetic trees were then constructed by Neighbor-Joining and the Kimura 2-parameters methods (MEGA software), followed by an analysis of the diversity of secondary structure. Then selection pressures acting on the L1/L2 genes were estimated by PAML software. Twenty-nine single nucleotide changes were observed in HPV-16 L1 sequences with 16/29 non-synonymous mutations and 13/29 synonymous mutations (six in alpha helix and two in beta turns). Seventeen single nucleotide changes were observed in HPV-16 L2 sequences with 8/17 non-synonymous mutations (one in beta turn) and 9/17 synonymous mutations. Twenty-four single nucleotide changes were observed in HPV-58 L1 sequences with 10/24 non-synonymous mutations and 14/24 synonymous mutations (eight in alpha helix and four in beta turn). Seven single nucleotide changes were observed in HPV-58 L2 sequences with 4/7 non-synonymous mutations and 3/7 synonymous mutations. The result of selective pressure analysis showed that most of these mutations were of positive selection. This study may help understand the intrinsic geographical relatedness and biological differences of HPV-16/HPV-58 and contributes further to research on their infectivity, pathogenicity, and vaccine strategy.