Whole-genome Phylogenetic Tree Research Articles

Bioinformatic identification of monkeypox virus phylogenetic gene trees that are representative of its whole-genome phylogenetic tree.

Phylogenetic analysis based on whole-genome sequences is the gold standard for monkeypox virus (MPXV) phylogeny. However, genomic epidemiology capability and capacity are lacking or limited in resource poor countries of sub-Saharan Africa. Therefore, these make real-time genome surveillance of MPXV virtually impossible. We hypothesized that phylogenetic analysis based on single, conserved genes will produce phylogenetic tree topology consistent with MPXV whole-genome phylogeny, thus serving as a reliable proxy to phylogenomic analysis. In this study, we analyzed 62 conserved MPXV genes and showed that Bayesian phylogenetic analysis based on five genes (OPG 066/E4L, OPG068/E6R, OPG079/I3L, OPG145/A18R, and OPG150/A23R) generated phylogenetic trees with 72.2-96.3% topology similarity index to the reference phylogenomic tree topology. Our results showed that phylogenetic analysis of the identified five genes singly or in combination can serve as surrogate for whole-genome phylogenetic analysis, and thus obviates the need for whole-genome sequencing and phylogenomic analysis in regions where genomic epidemiology competence and capacity are lacking or unavailable. This study is relevant to evolution and genome surveillance of MPXV in resource limited countries.

Molecular characterization of recombinant LSDV isolates from 2022 outbreak in Indonesia through phylogenetic networks and whole-genome SNP-based analysis

Lumpy skin disease (LSD) is a transboundary viral disease of cattle and water buffaloes caused by the LSD virus, leading to high morbidity, low mortality, and a significant economic impact. Initially endemic to Africa only, LSD has spread to the Middle East, Europe, and Asia in the past decade. The most effective control strategy for LSD is the vaccination of cattle with live-attenuated LSDV vaccines. Consequently, the emergence of two groups of LSDV strains in Asian countries, one closely related to the ancient Kenyan LSDV isolates and the second made of recombinant viruses with a backbone of Neethling-vaccine and field isolates, emphasized the need for constant molecular surveillance. This current study investigated the first outbreak of LSD in Indonesia in 2022. Molecular characterization of the isolate circulating in the country based on selected LSDV-marker genes: RPO30, GPCR, EEV glycoprotein gene, and B22R, as well as whole genome analysis using several analytical tools, indicated the Indonesia LSDV isolate as a recombinant of LSDV_Neethling_vaccine_LW_1959 and LSDV_NI-2490. The analysis clustered the Indonesia_LSDV with the previously reported LSDV recombinants circulating in East and Southeast Asia, but different from the recombinant viruses in Russia and the field isolates in South-Asian countries. Additionally, this study has demonstrated alternative accurate ways of LSDV whole genome analysis and clustering of isolates, including the recombinants, instead of whole-genome phylogenetic tree analysis. These data will strengthen our understanding of the pathogens’ origin, the extent of their spread, and determination of suitable control measures required.

Genomic characterization of two metagenome-assembled genomes of Tropheryma whipplei from China.

Whipple’s disease is a rare chronic systemic disease that affects almost any organ system of the body caused by the intracellular bacterium Tropheryma whipplei, which is found ubiquitously in the environment. Sequencing of the T. whipplei genome has revealed that it has a reduced genome (0.93 Mbp), a characteristic shared with other intracellular bacteria. Until our research started, 19 T. whipplei strains had been sequenced from cultures originated in France, Canada, and Germany. The genome of T. whipplei bacterium has not been studied in Asia yet. Here, two metagenome-assembled genomes (MAGs) of T. whipplei from China were reconstructed through metagenomic next-generation sequencing (mNGS) and genome binning. We also provided genomic insights into the geographical role and genomic features by analyzing the whole genome. The whole-genome phylogenetic tree was constructed based on single-nucleotide polymorphism (SNP) distance calculations and then grouped by distance similarity. The phylogenetic tree shows inconsistencies with geographic origins, thus suggesting that the variations in geographical origins cannot explain the phylogenetic relationships among the 21 T. whipplei strains. The two Chinese strains were closely related to each other, and also found to be related to strains from Germany (T. whipplei TW08/27) and France (T. whipplei Bcu26 and T. whipplei Neuro1). Furthermore, the Average Nucleotide Identity (ANI) matrix also showed no association between geographic origins and genomic similarities. The pan-genome analysis revealed that T. whipplei has a closed pan-genome composed of big core-genomes and small accessory genomes, like other intracellular bacteria. By examining the genotypes of the sequenced strains, all 21 T. whipplei strains were found to be resistant to fluoroquinolones, due to the genetic mutations in genes gyrA, gyrB, parC, and parE. The 21 T. Whipplei strains shared the same virulence factors, except for the alpC gene, which existed in 7 out of the 21 T. whipplei strains. When comparing 21 entire T. whipplei pan-genomes from various nations, it was discovered that the bacterium also possessed a closed genome, which was a trait shared by intracellular pathogens.

An outbreak of acute respiratory disease caused by HAdV-55 in Beijing, China, 2020.

Human adenoviruses (HAdVs) can cause acute respiratory diseases (ARDs) worldwide, and HAdV-55 is a reemergent pathogen in recent years. In the study, we investigated an outbreak of ARD at a school due to HAdV-55 in Beijing, China, during the early outbreak of coronavirus disease 2019 (COVID-19). The epidemic prevention team was dispatched to the school to collect epidemiologic data and nasopharyngeal samples. Then, real-time reverse transcription polymerase chain reaction (PCR) and multiplex PCR assays were used to detect severe acute respiratory syndrome coronavirus 2 and other respiratory pathogens, respectively. One representative HAdV-55 isolate was selected and submitted for whole-genome sequencing using a MiSeq system and the whole-genome phylogenetic tree was conducted based on the maximum likelihood method. The outbreak lasted from January 27 to February 6, 2020, and 108 students developed fever, among whom 60 (55.56%) cases were diagnosed with HAdV-55 infection in the laboratory using real-time PCR and 56 cases were hospitalized. All the confirmed cases hada fever and 11 cases (18.33%) presented with a fever above 39°C. Other main clinical symptoms included sore throat (43.33%) and headache (43.33%). We obtained and assembled the full genome of one isolate, BJ-446, with 34 761 nucleotides in length. HAdV-55 isolate BJ-446 was 99.85% identical to strain QS-DLL, which was the first HAdV-55 strain in China isolated from an ARD outbreak in Shanxi in 2006. One and four amino acid mutations were observed in the hexon gene and the coding region of L2 pV 40.1 kDa protein, respectively. We identified the first HAdV-55 infection associated with the ARD outbreak in Beijing since the emergence of COVID-19. The study suggests that improved surveillance of HAdV is needed, althoughCOVID-19 is still prevalent in the world.

Insights into complex infection by two Pectobacterium species causing potato blackleg and soft rot

Pectobacterium spp. are causative agents of blackleg and soft rot of potato. However, little is known about the relationship between the pathogenicity of mixed infections of different Pectobacterium spp. at different temperatures. In this study, two pectinolytic strains of Pectobacterium spp. were isolated from the same potato plant with typical symptoms of blackleg and identified as P. brasiliense and P. carotovorum by multilocus sequence analysis (MLSA), whole-genome phylogenetic tree construction, average nucleotide identity (ANI) analysis and digital DNA–DNA hybridization (dDDH). Plant cell wall degrading enzyme, including pectinases, cellulases and proteases, as the most important virulence factors, as well as pathogenicity toward potato tuber, were compared between the strains P. brasiliense BL-2 and P. carotovorum BL-4 at 28 ℃. The results showed that P. carotovorum had higher cell wall-degrading enzyme activities and brought more severe disease symptoms to potato tubers than P. brasiliense. Moreover, the pathogenicity of P. carotovorum and P. brasiliense increased with increasing temperature (20, 25, 28, 32 ℃). The pathogenicity was more severe when P. carotovorum strain BL-4 was co-inoculated with P. brasiliense strain BL-2, especially when the former exhibited an advantage in bacterial number at the initial time. The results of this study provide new insight for understanding the pathogenicity caused by mixed infections with different species of Pectobacterium spp., and they may provide some guidance for controlling potato blackleg and soft rot.

The success of particular Acinetobacter baumannii clones: accumulating resistance and virulence inside a sugary shield.

In Portugal, carbapenem-resistant Acinetobacter baumannii (CRAB) has been associated with ST98, ST103 and ST208 (Oxford Scheme, Oxf) and a clone has usually been associated with a particular period of time. These clonal shifts were primarily explained by an increased antimicrobial resistance profile. Here we explore genomic and biochemical differences among these and more recent clones, which could further explain the diversity and evolution of this species. A total of 116 CRAB isolates (2010-15), together with representatives of a previously described CRAB collection (4 isolates, 2001-06) were characterized by attenuated total reflection Fourier transform infrared spectroscopy (FTIR-ATR) and MLST. Representatives of different FTIR-ATR/MLST clusters were selected for WGS (n = 13), which allowed the in silico extraction of resistance and virulence genes, capsule locus and SNP analysis. A. baumannii clonal shifts of OXA-58-producing ST103Oxf (2001-04), OXA-40-producing ST98Oxf (2002-06), OXA-23-producing ST208Oxf (2006-10) and OXA-23-producing ST218Oxf (2010-15) were accompanied by an increase in AMR genes and virulence factors. FTIR-ATR clustering was congruent with sugar composition predicted from the capsular locus: a fucosamine cluster comprising ST98Oxf, ST103Oxf and a single ST218Oxf isolate; a pseudaminic acid cluster of ST208Oxf and ST1557Oxf isolates; and legionaminic acid, resembling the sialic acid from mammalian cells, in a cluster comprising ST218Oxf isolates. The whole-genome phylogenetic tree was congruent with MLST, with isolates presenting 5-28 938 SNPs. ST208Oxf and ST218Oxf presented ∼1900 SNPs while ST103Oxf and ST1557Oxf showed a greater number of SNPs (∼28 000). Clonal shifts of CRAB were promoted, in our country, by consecutive virulence and AMR gene pool enlargement, together with features increasing pathogen-host adaptation. Worldwide dominance of ST218Oxf is supported by the combination of high AMR and virulence levels.

Genomic analyses of multidrug-resistant Salmonella Indiana, Typhimurium, and Enteritidis isolates using MinION and MiSeq sequencing technologies.

We sequenced 25 isolates of phenotypically multidrug-resistant Salmonella Indiana (n = 11), Typhimurium (n = 8), and Enteritidis (n = 6) using both MinION long-read [SQK-LSK109 and flow cell (R9.4.1)] and MiSeq short-read (Nextera XT and MiSeq Reagent Kit v2) sequencing technologies to determine the advantages of each approach in terms of the characteristics of genome structure, antimicrobial resistance (AMR), virulence potential, whole-genome phylogeny, and pan-genome. The MinION reads were base-called in real-time using MinKnow 3.4.8 integrated with Guppy 3.0.7. The long-read-only assembly, Illumina-only assembly, and hybrid assembly pipelines of Unicycler 0.4.8 were used to generate the MinION, MiSeq, and hybrid assemblies, respectively. The MinION assemblies were highly contiguous compared to the MiSeq assemblies but lacked accuracy, a deficiency that was mitigated by adding the MiSeq short reads through the Unicycler hybrid assembly which corrected erroneous single nucleotide polymorphisms (SNPs). The MinION assemblies provided similar predictions of AMR and virulence potential compared to the MiSeq and hybrid assemblies, although they produced more total false negatives of AMR genotypes, primarily due to failure in identifying tetracycline resistance genes in 11 of the 19 MinION assemblies of tetracycline-resistant isolates. The MinION assemblies displayed a large genetic distance from their corresponding MiSeq and hybrid assemblies on the whole-genome phylogenetic tree, indicating that the lower read accuracy of MinION sequencing caused incorrect clustering. The pan-genome of the MinION assemblies contained significantly more accessory genes and less core genes compared to the MiSeq and hybrid assemblies, suggesting that although these assemblies were more contiguous, their sequencing errors reduced accurate genome annotations. Our research demonstrates that MinION sequencing by itself provides an efficient assessment of the genome structure, antimicrobial resistance, and virulence potential of Salmonella; however, it is not sufficient for whole-genome phylogenetic and pan-genome analyses. MinION in combination with MiSeq facilitated the most accurate genomic analyses.

A30 Avian influenza viruses in wild birds: Virus evolution in a multi-host ecosystem

Wild ducks and gulls are the major reservoirs for avian influenza A viruses (AIVs). The mechanisms that drive AIV evolution are complex at sites where various duck and gull species from multiple flyways breed, winter, or stage. The Republic of Georgia is located at the intersection of three migratory flyways: the Central Asian Flyway, East Asian/East African Flyway, and Black Sea/Mediterranean Flyway. For six consecutive years (2010–6), we collected AIV samples from various duck and gull species that breed, migrate, and overwinter in Georgia. We found substantial subtype diversity of viruses that varied in prevalence from year to year. Low pathogenic (LP)AIV subtypes included H1N1, H2N3, H2N5, H2N7, H3N8, H4N2, H6N2, H7N3, H7N7, H9N1, H9N3, H10N4, H10N7, H11N1, H13N2, H13N6, H13N8, and H16N3, plus two H5N5 and H5N8 highly pathogenic (HP)AIVs belonging to clade 2.3.4.4. Whole-genome phylogenetic trees showed significant host species lineage restriction for nearly all gene segments and significant differences for LPAIVs among different host species in observed reassortment rates, as defined by quantification of phylogenetic incongruence, and in nucleotide diversity. Hemagglutinin clade 2.3.4.4 H5N8 viruses, circulated in Eurasia during 2014–5 did not reassort, but analysis after its subsequent dissemination during 2016–7 revealed reassortment in all gene segments except NP and NS. Some virus lineages appeared to be unrelated to AIVs in wild bird populations in other regions with maintenance of local AIV viruses in Georgia, whereas other lineages showed considerable genetic inter-relationship with viruses circulating in other parts of Eurasia and Africa, despite relative under-sampling in the area.

Clustering of Vibrio parahaemolyticus Isolates Using MLST and Whole-Genome Phylogenetics and Protein Motif Fingerprinting.

Vibrio parahaemolyticus is a ubiquitous and abundant member of native microbial assemblages in coastal waters and shellfish. Though V. parahaemolyticus is predominantly environmental, some strains have infected human hosts and caused outbreaks of seafood-related gastroenteritis. In order to understand differences among clinical and environmental V. parahaemolyticus strains, we used high quality DNA sequencing data to compare the genomes of V. parahaemolyticus isolates (n = 43) from a variety of geographic locations and clinical and environmental sample matrices. We used phylogenetic trees inferred from multilocus sequence typing (MLST) and whole-genome (WG) alignments, as well as a novel classification and genome clustering approach that relies on protein motif fingerprints (MFs), to assess relationships between V. parahaemolyticus strains and identify novel molecular targets associated with virulence. Differences in strain clustering at more than one position were observed between the MLST and WG phylogenetic trees. The WG phylogeny had higher support values and strain resolution since isolates of the same sequence type could be differentiated. The MF analysis revealed groups of protein motifs that were associated with the pathogenic MLST type ST36 and a large group of clinical strains isolated from human stool. A subset of the stool and ST36-associated protein motifs were selected for further analysis and the motif sequences were found in genes with a variety of functions, including transposases, secretion system components and effectors, and hypothetical proteins. DNA sequences associated with these protein motifs are candidate targets for future molecular assays in order to improve surveys of pathogenic V. parahaemolyticus in the environment and seafood.

Streptococcus chenjunshii sp. nov. isolated from feces of Tibetan antelopes.

Three Gram-stain-positive, catalase-negative, α-haemolytic, chain-forming and coccus-shaped microorganisms (strains Z15T, Z1 and Z2) were isolated from feces of Tibetan antelopes collected from the Qinghai-Tibet plateau, PR China. The results of 16S rRNA gene sequence studies indicated that Z15T shared 94.5, 93.1 and 92.2 % similarity with Streptococcus pantholopis DSM 102135T, Streptococcus ursoris NUM 1615T and Streptococcus dentapri NUM 1529T, respectively. rpoB and groEL-based sequence analysis of our three novel isolates revealed interspecies divergence of 16.7 and 14.3 % from Streptococcus pantholopis DSM 102135T. The genomic DNA G+C content of Z15T is 42.48 mol%. Z15T has an average nucleotide identity (ANI) value of 81.19 % with S. pantholopis DSM 102135T and a DNA-DNA relatedness value of less than 70 % in the in-silico DNA-DNA hybridization (isDDH) with other species of genus Streptococcus deposited in the GenBank database. A whole-genome phylogenetic tree based on 246 core genes of 78 genomes of members of the genus Streptococcusindicated that Z15T represents a member of genus Streptococcus but one well separated from the currently recognized species. Z15T contains C18 : 1ω7c (25.5 %), C18 : 1ω9c (19.6 %), C16 : 0 (17.5 %) and C16 : 1ω9c (13.3 %) as its major cellular fatty acids. According to the morphological, biochemical and molecular phylogenetic features of the three novel isolates, they represent a novel species of the genus Streptococcus, and Streptococcus chenjunshii sp. nov. is thus proposed. The type strain is Z15T (=CGMCC 1.16529=DSM 106182).

Avian Influenza Viruses in Wild Birds: Virus Evolution in a Multihost Ecosystem.

ABSTRACTWild ducks and gulls are the major reservoirs for avian influenza A viruses (AIVs). The mechanisms that drive AIV evolution are complex at sites where various duck and gull species from multiple flyways breed, winter, or stage. The Republic of Georgia is located at the intersection of three migratory flyways: the Central Asian flyway, the East Africa/West Asia flyway, and the Black Sea/Mediterranean flyway. For six complete study years (2010 to 2016), we collected AIV samples from various duck and gull species that breed, migrate, and overwinter in Georgia. We found a substantial subtype diversity of viruses that varied in prevalence from year to year. Low-pathogenic AIV (LPAIV) subtypes included H1N1, H2N3, H2N5, H2N7, H3N8, H4N2, H6N2, H7N3, H7N7, H9N1, H9N3, H10N4, H10N7, H11N1, H13N2, H13N6, H13N8, and H16N3, and two highly pathogenic AIVs (HPAIVs) belonging to clade 2.3.4.4, H5N5 and H5N8, were found. Whole-genome phylogenetic trees showed significant host species lineage restriction for nearly all gene segments and significant differences in observed reassortment rates, as defined by quantification of phylogenetic incongruence, and in nucleotide sequence diversity for LPAIVs among different host species. Hemagglutinin clade 2.3.4.4 H5N8 viruses, which circulated in Eurasia during 2014 and 2015, did not reassort, but analysis after their subsequent dissemination during 2016 and 2017 revealed reassortment in all gene segments except NP and NS. Some virus lineages appeared to be unrelated to AIVs in wild bird populations in other regions, with maintenance of local AIVs in Georgia, whereas other lineages showed considerable genetic interrelationships with viruses circulating in other parts of Eurasia and Africa, despite relative undersampling in the area.IMPORTANCE Waterbirds (e.g., gulls and ducks) are natural reservoirs of avian influenza viruses (AIVs) and have been shown to mediate the dispersal of AIVs at intercontinental scales during seasonal migration. The segmented genome of influenza viruses enables viral RNA from different lineages to mix or reassort when two viruses infect the same host. Such reassortant viruses have been identified in most major human influenza pandemics and several poultry outbreaks. Despite their importance, we have only recently begun to understand AIV evolution and reassortment in their natural host reservoirs. This comprehensive study illustrates AIV evolutionary dynamics within a multihost ecosystem at a stopover site where three major migratory flyways intersect. Our analysis of this ecosystem over a 6-year period provides a snapshot of how these viruses are linked to global AIV populations. Understanding the evolution of AIVs in the natural host is imperative to mitigating both the risk of incursion into domestic poultry and the potential risk to mammalian hosts, including humans.

Isolation and characterization of Streptococcus respiraculi sp. nov. from Marmota himalayana (Himalayan marmot) respiratory tract.

Two bacterial strains were individually isolated from Marmota himalayana respiratory tracts; the animals were from the Tibet-Qinghai Plateau, PR China. The isolates were Gram-stain-positive, catalase-negative, coccus-shaped, chain-forming organisms. Analysis of 16S rRNA gene sequences indicated that the type strain HTS25T shared 98.0, 97.4, 97.2 and 97.1 % similarity with Streptococcus cuniculi, Streptococcus acidominimus, Streptococcus marmotae and Streptococcus himalayensis respectively. Sequence analysis of the sodA and rpoB genes indicated that HTS25T was closely related to S. marmotae (similarities of 94.7 and 91.4 % respectively). Analysis of groEL sequences showed interspecies similarity of 84.8 % between HTS25T and S. himalayensis. A whole-genome phylogenetic tree reconstructed from 81 core genes from the genomes of 17 members of the genus Streptococcus was used to validate that HTS25T forms a distinct subline from other recognized species of the genus Streptococcus. DNA-DNA hybridization of HTS25T showed a maximum estimated DNA reassociation value of 32.1 % to Streptococcus cuniculi CCUG 65085T. On the basis of the results of phenotypic and phylogenetic analyses, we propose that the two isolates be classified as representing a novel species of the genus Streptococcus, named Streptococcus respiraculi sp. nov. The type strain is HTS25T (=DSM 101998T=CGMCC 1.15531T). The genome of Streptococcus respiraculi sp. nov. strain HTS25T (2 067 971 bp) contains 2001 genes with an average DNA G+C content of 42.7 mol%.

Genome comparison of different Zymomonas mobilis strains provides insights on conservation of the evolution

Zymomonas mobilis has the special Entner-Doudoroff (ED) pathway and it has excellent industrial characteristics, including low cell mass formation, high-specific productivity,ethanol yield, notable ethanol tolerance and wide pH range, a relatively small genome size. In this study, the genome sequences of NRRL B-14023 and NRRL B-12526 were sequenced and compared with other strains to explore their evolutionary relationships and the genetic basis of Z. mobilis. The comparative genomic analyses revealed that the 8 strains share a conserved core chromosomal backbone. ZM4, NRRL B-12526, NRRL B-14023, NCIMB 11163 and NRRL B-1960 share 98% sequence identity across the whole genome sequences. Highly similar plasmids and CRISPR repeats were detected in these strains. A whole-genome phylogenetic tree of the 8 strains indicated that NRRL B-12526, NRRL B-14023 and ATCC 10988 had a close evolutionary relationship with the strain ZM4. Furthermore, strains ATCC29191 and ATCC29192 had distinctive CRISPR with a far distant relationship. The size of the pan-genome was 1945 genes, including 1428 core genes and 517 accessory genes. The genomes of Z. mobilis were highly conserved; particularly strains ZM4, NRRL B-12526, NRRL B-14023, NCIMB 11163 and NRRL B-1960 had a close genomic relationship. This comparative study of Z. mobilis presents a foundation for future functional analyses and applications.

Enterococcus wangshanyuanii sp. nov., isolated from faeces of yaks (Bos grunniens)

The taxonomic position of four phenotypically closely related strains isolated from faecal samples of yaks (Bos grunniens) collected from the Qinghai-Tibetan plateau, China, was determined by a polyphasic approach. The strains were non-spore-forming, non-motile Gram-stain-positive, ovoid cocci, occurring predominantly in pairs and short chains or in irregular clusters. The 16S rRNA gene of strain MN05T was related phylogenetically to those of Enterococcushaemoperoxidus, Enterococcusrotai, Enterococcusquebecensis, Enterococcusplantarum, Enterococcuscrotali, Enterococcusmoraviensis, Enterococcussilesiacus, Enterococcuscaccae, Enterococcustermitis, Enterococcusureasiticus and Enterococcusureilyticus, all belonging to the Enterococcusfaecalis species group. The sequence similarities of three selected genes of MN05T to those of the type strains of phylogenetically related species were measured, with values within the range of 99.2-99.5 % (16S rRNA gene), 90.0-97.3 % (rpoA) and 80.0-85.3 % (pheS), respectively. The genome of MN05T (3 842 361 bp) contained 4299 genes with a DNA G+C content of 37.5 mol%. A whole-genome phylogenetic tree based on 808 core genes confirmed that MN05T belongs to a distinct lineage, well separated from all recognized species of the Enterococcusfaecalis species group. DNA-DNA hybridization in silico showed that MN05T displayed less than 70 % DNA-DNA relatedness with the other 13 species of the Enterococcusfaecalis species group. Moreover, their phenotypic features distinguished the four strains from the other species of the Enterococcusfaecalis species group. Based upon these data obtained from the polyphasic characterization performed in the present study, a novel species of the genus Enterococcus, Enterococcus wangshanyuanii sp. nov., is proposed, with the type strain MN05T (=DSM 104047T=CGMCC 1.15942T).

Streptococcus himalayensis sp. nov., isolated from the respiratory tract of Marmota himalayana.

Five strains of Gram-positive-staining, catalase-negative, coccus-shaped, chain-forming organisms isolated separately from the respiratory tracts of five Marmota himalayana animals in the Qinghai-Tibet Plateau of China were subjected to phenotypic and molecular taxonomic analyses. Comparative analysis of the 16S rRNA gene indicated that these singular organisms represent a new member of the genus Streptococcus, being phylogenetically closest to Streptococcus marmotae DSM 101995T (98.4 % similarity). The groEL, sodA and rpoB sequence analysis showed interspecies similarity values between HTS2T and Streptococcus. marmotae DSM 101995T, its closest phylogenetic relative based on 16S rRNA gene sequences, of 98.2, 78.8 and 93.7 %, respectively. A whole-genome phylogenetic tree built from 82 core genes of genomes from 16 species of the genus Streptococcus validated that HTS2T forms a distinct subline and exhibits specific phylogenetic affinity with S. marmotae. In silico DNA-DNA hybridization of HTS2T showed an estimated DNA reassociation value of 40.5 % with Streptococcus. marmotae DSM 101995T. On the basis of their phenotypic characteristics and phylogenetic findings, it is proposed that the five isolates be classified as representatives of a novel species of the genus Streptococcus, Streptococcus himalayensis sp. nov. The type strain is HTS2T (=DSM 101997T=CGMCC 1.15533T). The genome of Streptococcus himalayensis sp. nov. strain HTS2T contains 2195 genes with a size of 2 275 471 bp and a mean DNA G+C content of 41.3 mol%.

Acinetobacter seifertii Isolated from China: Genomic Sequence and Molecular Epidemiology Analyses.

Clinical infections caused by Acinetobacter spp. have increasing public health concerns because of their global occurrence and ability to acquire multidrug resistance. Acinetobacter calcoaceticus–Acinetobacter baumannii (ACB) complex encompasses A. calcoaceticus, A. baumannii, A. pittii (formerly genomic species 3), and A nosocomial (formerly genomic species 13TU), which are predominantly responsible for clinical pathogenesis in the Acinetobacter genus.In our previous study, a putative novel species isolated from 385 non-A. baumannii spp. strains based on the rpoB gene phylogenetic tree was reported. Here, the putative novel species was identified as A. seifertii based on the whole-genome phylogenetic tree. A. seifertii was recognized as a novel member of the ACB complex and close to A. baumannii and A. nosocomials. Furthermore, we studied the characteristics of 10 A. seifertii isolates, which were distributed widely in 6 provinces in China and mainly caused infections in the elderly or children. To define the taxonomic status and characteristics, the biochemical reactions, antimicrobial susceptibility testing, pulsed field gel electrophoresis (PFGE), multilocus sequence typing (MLST), and whole-genome sequence analysis were performed.The phenotypic characteristics failed to distinguish A. serfertii from other species in the ACB complex. Most of the A. seifertii isolates were susceptible to antibiotics commonly used for nosocomial Acinetobacter spp. infections, but one isolate (strain A362) was resistant to ampicillin/sulbactam, ceftazidime and amikacin. The different patterns of MLST and PFGE suggested that the 10 isolates were not identical and lacked clonal relatedness.Our study reported for the first time the molecular epidemiological and genomic features of widely disseminated A. seifertii in China. These observations could enrich the knowledge of infections caused by non-A. baumannii and may provide a scientific basis for future clinical treatment.

Origin and Possible Genetic Recombination of the Middle East Respiratory Syndrome Coronavirus from the First Imported Case in China: Phylogenetics and Coalescence Analysis.

ABSTRACTThe Middle East respiratory syndrome coronavirus (MERS-CoV) causes a severe acute respiratory tract infection with a high fatality rate in humans. Coronaviruses are capable of infecting multiple species and can evolve rapidly through recombination events. Here, we report the complete genomic sequence analysis of a MERS-CoV strain imported to China from South Korea. The imported virus, provisionally named ChinaGD01, belongs to group 3 in clade B in the whole-genome phylogenetic tree and also has a similar tree topology structure in the open reading frame 1a and -b (ORF1ab) gene segment but clusters with group 5 of clade B in the tree constructed using the S gene. Genetic recombination analysis and lineage-specific single-nucleotide polymorphism (SNP) comparison suggest that the imported virus is a recombinant comprising group 3 and group 5 elements. The time-resolved phylogenetic estimation indicates that the recombination event likely occurred in the second half of 2014. Genetic recombination events between group 3 and group 5 of clade B may have implications for the transmissibility of the virus.

Comparative genomics of Klebsiella pneumoniae strains with different antibiotic resistance profiles.

There is a global emergence of multidrug-resistant (MDR) strains of Klebsiella pneumoniae, a Gram-negative enteric bacterium that causes nosocomial and urinary tract infections. While the epidemiology of K. pneumoniae strains and occurrences of specific antibiotic resistance genes, such as plasmid-borne extended-spectrum β-lactamases (ESBLs), have been extensively studied, only four complete genomes of K. pneumoniae are available. To better understand the multidrug resistance factors in K. pneumoniae, we determined by pyrosequencing the nearly complete genome DNA sequences of two strains with disparate antibiotic resistance profiles, broadly drug-susceptible strain JH1 and strain 1162281, which is resistant to multiple clinically used antibiotics, including extended-spectrum β-lactams, fluoroquinolones, aminoglycosides, trimethoprim, and sulfamethoxazoles. Comparative genomic analysis of JH1, 1162281, and other published K. pneumoniae genomes revealed a core set of 3,631 conserved orthologous proteins, which were used for reconstruction of whole-genome phylogenetic trees. The close evolutionary relationship between JH1 and 1162281 relative to other K. pneumoniae strains suggests that a large component of the genetic and phenotypic diversity of clinical isolates is due to horizontal gene transfer. Using curated lists of over 400 antibiotic resistance genes, we identified all of the elements that differentiated the antibiotic profile of MDR strain 1162281 from that of susceptible strain JH1, such as the presence of additional efflux pumps, ESBLs, and multiple mechanisms of fluoroquinolone resistance. Our study adds new and significant DNA sequence data on K. pneumoniae strains and demonstrates the value of whole-genome sequencing in characterizing multidrug resistance in clinical isolates.

Genetic diversity of the ordinary strain of Potato virus Y (PVY) and origin of recombinant PVY strains.

The ordinary strain of Potato virus Y (PVY), PVY(O), causes mild mosaic in tobacco and induces necrosis and severe stunting in potato cultivars carrying the Ny gene. A novel substrain of PVY(O) was recently reported, PVY(O)-O5, which is spreading in the United States and is distinguished from other PVY(O) isolates serologically (i.e., reacting to the otherwise PVY(N)-specific monoclonal antibody 1F5). To characterize this new PVY(O)-O5 subgroup and address possible reasons for its continued spread, we conducted a molecular study of PVY(O) and PVY(O)-O5 isolates from a North American collection of PVY through whole-genome sequencing and phylogenetic analysis. In all, 44 PVY(O) isolates were sequenced, including 31 from the previously defined PVY(O)-O5 group, and subjected to whole-genome analysis. PVY(O)-O5 isolates formed a separate lineage within the PVY(O) genome cluster in the whole-genome phylogenetic tree and represented a novel evolutionary lineage of PVY from potato. On the other hand, the PVY(O) sequences separated into at least two distinct lineages on the whole-genome phylogenetic tree. To shed light on the origin of the three most common PVY recombinants, a more detailed phylogenetic analysis of a sequence fragment, nucleotides 2,406 to 5,821, that is present in all recombinant and nonrecombinant PVY(O) genomes was conducted. The analysis revealed that PVY(N:O) and PVY(N-Wi) recombinants acquired their PVY(O) segments from two separate PVY(O) lineages, whereas the PVY(NTN) recombinant acquired its PVY(O) segment from the same lineage as PVY(N:O). These data suggest that PVY(N:O) and PVY(N-Wi) recombinants originated from two separate recombination events involving two different PVY(O) parental genomes, whereas the PVY(NTN) recombinants likely originated from the PVY(N:O) genome via additional recombination events.

Gene content phylogeny of herpesviruses.

Clusters of orthologous groups [COGs; Tatusov, R. L., Koonin, E. V. & Lipman, D. J. (1997) Science 278, 631-637] were identified for a set of 13 completely sequenced herpesviruses. Each COG represented a family of gene products conserved across several herpes genomes. These families were defined without using an arbitrary threshold criterion based on sequence similarity. The COG technique was modified so that variable stringency in COG construction was possible. High stringencies identify a core set of highly conserved genes. Varying COG stringency reveals differences in the degree of conservation between functional classes of genes. The COG data were used to construct whole-genome phylogenetic trees based on gene content. These trees agree well with trees based on other methods and are robust when tested by bootstrap analysis. The COG data also were used to construct a reciprocal tree that clustered genes with similar phylogenetic profiles. This clustering may give clues to genes with related functions or with related histories of acquisition and loss during herpesvirus evolution.