Draft Genome Sequence of an Extremely Drug-Resistant KPC-Producing Klebsiella pneumoniae ST258 Epidemic Strain

Extremely drug-resistant (XDR) Klebsiella pneumoniae carbapenemase-producing clone ST258 has rapidly disseminated worldwide. We report here the draft genome sequence of the K. pneumoniae ST258 XDR clinical strain from Israel.     

Complex Population Structure and Virulence Differences among Serotype 2 Streptococcus suis Strains Belonging to Sequence Type 28

Streptococcus suis is a major swine pathogen and a zoonotic agent. Serotype 2 strains are the most frequently associated with disease. However, not all serotype 2 lineages are considered virulent. Indeed, sequence type (ST) 28 serotype 2 S. suis strains have been described as a homogeneous group of low virulence. However, ST28 strains are often isolated from diseased swine in some countries, and at least four human ST28 cases have been reported. Here, we used whole-genome sequencing and animal infection models to test the hypothesis that the ST28 lineage comprises strains of different genetic backgrounds and different virulence. We used 50 S. suis ST28 strains isolated in Canada, the United States and Japan from diseased pigs, and one ST28 strain from a human case isolated in Thailand. We report a complex population structure among the 51 ST28 strains. Diversity resulted from variable gene content, recombination events and numerous genome-wide polymorphisms not attributable to recombination. Phylogenetic analysis using core genome single-nucleotide polymorphisms revealed four discrete clades with strong geographic structure, and a fifth clade formed by US, Thai and Japanese strains. When tested in experimental animal models, strains from this latter clade were significantly more virulent than a Canadian ST28 reference strain, and a closely related Canadian strain. Our results highlight the limitations of MLST for both phylogenetic analysis and virulence prediction and raise concerns about the possible emergence of ST28 strains in human clinical cases.     

Population Structure and Antimicrobial Resistance Profiles of Streptococcus suis Serotype 2 Sequence Type 25 Strains

Strains of serotype 2 Streptococcus suis are responsible for swine and human infections. Different serotype 2 genetic backgrounds have been defined using multilocus sequence typing (MLST). However, little is known about the genetic diversity within each MLST sequence type (ST). Here, we used whole-genome sequencing to test the hypothesis that S. suis serotype 2 strains of the ST25 lineage are genetically heterogeneous. We evaluated 51 serotype 2 ST25 S. suis strains isolated from diseased pigs and humans in Canada, the United States of America, and Thailand. Whole-genome sequencing revealed numerous large-scale rearrangements in the ST25 genome, compared to the genomes of ST1 and ST28 S. suis strains, which result, among other changes, in disruption of a pilus island locus. We report that recombination and lateral gene transfer contribute to ST25 genetic diversity. Phylogenetic analysis identified two main and distinct Thai and North American clades grouping most strains investigated. These clades also possessed distinct patterns of antimicrobial resistance genes, which correlated with acquisition of different integrative and conjugative elements (ICEs). Some of these ICEs were found to be integrated at a recombination hot spot, previously identified as the site of integration of the 89K pathogenicity island in serotype 2 ST7 S. suis strains. Our results highlight the limitations of MLST for phylogenetic analysis of S. suis, and the importance of lateral gene transfer and recombination as drivers of diversity in this swine pathogen and zoonotic agent.     

Genomic insights into evolution of extensive drug resistance in Stenotrophomonas maltophilia complex

We report first complete genomic investigation of extensive drug resistance (XDR) in a nosocomial Stenotrophomonas maltophilia complex strain that is resistant to mainstream drugs (trimethoprim/sulfamethoxazole and levofloxacin). Comprehensive genomic investigation revealed its exclusive fourteen dynamic regions and highly enriched resistome comprising of two sulfonamide resistance genes on two diverse super-integrons of chromosomal origin. In addition, both these integrons harbour array of antibiotic resistance and commonly used disinfectant's resistance genes linked to ISCR elements. Isolation of a novel XDR strain from Indian tertiary care unit belonging to novel ST with diverse array of resistance genes on ISCR linked super-integrons indicates extent and nature of selection pressure in hospitals. Since, repetitive elements have major role in their spread and due to limitations of draft genomes, there is an urgent need to employ complete genome-based investigation for tracking the emergence of XDR at global level and designing strategies of antimicrobial stewardship and disinfection.ImportanceHospital settings in India have one of the highest usages of antimicrobials and a heavy patient load. We hereby report a novel clinical isolate of S. maltophilia complex with two super-integrons that harbour array of antimicrobial resistance genes along with biocide and heavy metal resistance genes. Further, the presence of ISCR type of transposable elements on both the integrons indicates their propensity to transfer resistome while their chromosomal origin suggests possibilities for further genomic/phenotypic complexities according to selection pressure. Such complex mobile cassettes in a novel strain is a potential threat to global health care. Hence, to understand the evolution of opportunistic nosocomial pathogen, there is an urgent need to employ cost-effective long read technologies to keep vigilance on novel and XDR pathogens in populous countries. There is also need for surveillance of the usage of disinfectants and other antimicrobials for environmental hygiene and linked/rapid co-evolution of XDR in nosocomial pathogens.Repositories: Complete genome sequence of Stenotrophomonas maltophilia SM866: CP031058.     

Comparative genomics inferred two distinct populations of piscine pathogenic Streptococcus agalactiae,serotype Ia ST7 and serotype III ST283, in Thailand and Vietnam

The genomes of Streptococcus agalactiae (group B streptococcus; GBS) collected from diseased fish in Thailand and Vietnam over a nine-year period (2008–2016) were sequenced and compared (n = 21). Based on capsular serotype and multilocus sequence typing (MLST), GBS isolates are divided into 2 groups comprised of i) serotype Ia; sequence type (ST)7 and ii) serotype III; ST283. Population structure inferred by core genome (cg)MLST and Bayesian clustering analysis also strongly indicated distribution of two GBS populations in both Thailand and Vietnam. Deep phylogenetic analysis implied by CRISPR array's spacer diversity was able to cluster GBS isolates according to their temporal and geographic origins, though ST7 has varying CRISPR1-spacer profiles when compared to ST283 strains. Based on overall genotypic features, Thai ST283 strains were closely related to the Singaporean ST283 strain causing foodborne illness in humans in 2015, thus, signifying zoonotic potential of this GBS population in the country.     

Discrimination among <Emphasis Type="Italic">Bacillus thuringiensis</Emphasis> H serotypes,serovars and strains based on 16S rRNA, <Emphasis Type="Italic">gyrB</Emphasis> and <Emphasis Type="Italic">aroE</Emphasis> gene sequence analyses

Our aim was to investigate the capability of each of three genes, 16S rRNA, gyrB and aroE, to discriminate, first, among Bacillus thuringiensis H serotypes; second, among B. thuringiensis serovars from the same H serotype; and third, among B. thuringiensis strains from the same serovar. The 16S rRNA, gyrB and aroE genes were amplified from 21 B. thuringiensis H serotypes and their nucleotide sequences determined. Additional strains from four B. cereus sensu lato species were included for comparison purposes. These sequences were pair-wise compared and phylogenetic relationships were revealed. Each of the three genes under study could discriminate among B. thuringiensis H serotypes. The gyrB and aroE genes showed a discriminatory power among B. thuringiensis H serotypes up to nine fold greater than that of the 16S rRNA gene. The gyrB gene was retained for subsequent analyses to discriminate B. thuringiensis serovars from the same H serotype and to discriminate strains from same serovar. A total of 42 B. thuringiensis strains, which encompassed 25 serovars from 12 H serotypes, were analyzed. The gyrB gene nucleotide sequences were different enough as to be sufficient to discriminate among B. thuringiensis serovars from the same H serotype and among B. thuringiensis strains from the same serovar. Electronic supplementary material  The online version of this article (doi:) contains supplementary material, which is available to authorized users.     

Chromosomal Rearrangements between Serotype A and D Strains in Cryptococcus neoformans

Cryptococcus neoformans is a major human pathogenic fungus that can cause meningoencephalitis in immunocompromised hosts. It contains two divergent varieties, var. grubii (serotype A) and var. neoformans (serotype D), as well as hybrids (serotype AD) between these two varieties. In this study, we investigated the extent of chromosomal rearrangements between the two varieties, estimated the effects of chromosomal rearrangements on recombination frequencies, and surveyed the potential polymorphisms of the rearrangements among natural strains of the three serotypes. Through the analyses of two sequenced genomes from strains H99 (representing var. grubii) and JEC21 (representing var. neoformans), we revealed a total of 32 unambiguous chromosome rearrangements, including five translocations, nine simple inversions, and 18 complex rearrangements. Our analyses identified that overall, rearranged regions had recombination frequencies about half of those around syntenic regions. Using a direct PCR screening strategy, we examined the potential polymorphisms of 11 rearrangements among 64 natural C. neoformans strains from five countries. We found no polymorphism within var. neoformans and very limited polymorphism within var. grubii. However, strains of serotype AD showed significant polymorphism, consistent with their hybrid origins coupled with differential loss of heterozygosity. We discuss the implications of these results on the genome structure, ecology, and evolution of C. neoformans.     

<Emphasis Type="Italic">Mycobacterium tuberculosis</Emphasis> mutants with multidrug resistance: History of origin,genetic and molecular mechanisms of resistance,and emerging challenges

The review summarizes the data on the Mycobacterium tuberculosis mutations that lead to multidrug resistance (MDR) to various antibiotics. MDR strains arose over the past 30 years as a variety of antituberculosis drugs were introduced in medicine, and they largely discount the results of chemotherapy for tuberculosis. The most dangerous of them are strains with extensive drug resistance (XDR), which are resistant to four or five different drugs on average. The molecular mechanisms that make a strain resistant are considered. XDR and MDR strains result from successive and usually independent resistance mutations, which arise in various regions of the mycobacterial genome. In addition, the formation of resistant strains is affected by the phenomenon of tolerance and mycobacterial latency in infected tissues.     

Homologous recombination mediated by two palindromic repeated sequences in the mitochondrial genome of Oryza

Palindromic repeated sequences (PRSs) are distributed in at least ten regions of the mitochondrial (mt) genome of rice and are, apparently, mobile. In the present study, we examined the possibility of homologous recombination via some PRSs during the course of evolution of Oryza. We first performed Southern hybridization of the DNA from 11 species (18 strains) of Oryza in order to identify the distribution of PRSs in the mitochondrial genome of Oryza. The hybridization patterns revealed genome type-specific and/or species-specific variations. We speculated that homologous recombination via some PRSs might have made a contribution to such variations. After subsequent polymerase chain reaction, Southern hybridization and sequencing, we concluded that homologous recombination mediated by two PRSs occurred in the mtDNA of Oryza after divergence of the BB genome type and the other genome types of Oryza. Evidence was obtained that some PRSs were involved in both insertion and recombination events during the evolution of Oryza. Our results indicate, therefore, that PRSs have contributed considerably to the polymorphism of Oryza mtDNAs.     

Characterization of a putative colonization factor (PCFO166) of enterotoxigenic Escherichia coli of serogroup O166

Enterotoxigenic Escherichia coli (ETEC) of serogroup O166 gave mannose-resistant haemagglutination (MRHA) with bovine and human erythrocytes. The strains did not react with antisera prepared against the known colonization factors CFA/I, CFA/II, CFA/III, CFA/IV and PCFO159:H4. Strain E7476 of serotype O166:H27, which produced heat-stable enterotoxin (ST), was examined initially. It produced fimbriae about 7 nm in diameter. On SDS-PAGE two possible fimbrial polypeptides of molecular mass 15.5 and 17.0 kDa were seen. When variants of strain E7476 were isolated, loss of ST and MRHA together was associated with loss of a 98 MDa plasmid, while loss of ST alone correlated with plasmid deletion. An absorbed anti-strain E7476 antiserum reacted specifically with the 15.5 and 17.0 kDa polypeptides in Western immunoblotting and bound to the intact fimbriae by immuno-electron microscopy. When this antiserum was used in an ELISA to examine other strains of serogroup O166, a positive reaction was obtained with all the ST- and MRHA-positive strains. One strain of serotype O71:H27 and two strains of serotype O98:H- also reacted with the absorbed anti-strain E7476 antiserum. The antiserum did not react with ETEC carrying known colonization factors. E. coli K12 and a number of E. coli of different serotypes carrying a plasmid coding for ST transferred from strain E7476, all gave MRHA and reacted with the absorbed anti-strain E7476 antiserum. The term putative colonization factor O166 (PCFO166) is proposed to describe the adhesive factor(s) on ETEC of serogroup O166 because of the similarity of properties with those of known colonization factors.     

Genome wide evolutionary analyses reveal serotype specific patterns of positive selection in selected Salmonella serotypes

Background  

The bacterium Salmonella enterica includes a diversity of serotypes that cause disease in humans and different animal species. Some Salmonella serotypes show a broad host range, some are host restricted and exclusively associated with one particular host, and some are associated with one particular host species, but able to cause disease in other host species and are thus considered "host adapted". Five Salmonella genome sequences, representing a broad host range serotype (Typhimurium), two host restricted serotypes (Typhi [two genomes] and Paratyphi) and one host adapted serotype (Choleraesuis) were used to identify core genome genes that show evidence for recombination and positive selection.     

Multilocus Sequence Typing of Listeria monocytogenes by Use of Hypervariable Genes Reveals Clonal and Recombination Histories of Three Lineages

In an attempt to develop a method to discriminate among isolates of Listeria monocytogenes, the sequences of all of the annotated genes from the fully sequenced strain L. monocytogenes EGD-e (serotype 1/2a) were compared by BLASTn to a file of the unfinished genomic sequence of L. monocytogenes ATCC 19115 (serotype 4b). Approximately 7% of the matching genes demonstrated 90% or lower identity between the two strains, and the lowest observed identity was 80%. Nine genes (hisJ, cbiE, truB, ribC, comEA, purM, aroE, hisC, and addB) in the 80 to 90% identity group and two genes (gyrB and rnhB) with approximately 97% identity were selected for multilocus sequence analysis in two sets of L. monocytogenes isolates (a 15-strain diversity set and a set of 19 isolates from a single food-processing plant). Based on concatenated sequences, a total of 33 allotypes were differentiated among the 34 isolates tested. Population genetics analyses revealed three lineages of L. monocytogenes that differed in their history of apparent recombination. Lineage I appeared to be completely clonal, whereas representatives of the other lineages demonstrated evidence of horizontal gene transfer and recombination. Although most of the gene sequences for lineage II strains were distinct from those of lineage I, a few strains with the majority of genes characteristic of lineage II had some that were characteristic of lineage I. Genes from lineage III organisms were mostly similar to lineage I genes, with instances of genes appearing to be mosaics with lineage II genes. Even though lineage I and lineage II generally demonstrated very distinct sequences, the sequences for the 11 selected genes demonstrated little discriminatory power within each lineage. In the L. monocytogenes isolate set obtained from one food-processing plant, lineage I and lineage II were found to be almost equally prevalent. While it appears that different lineages of L. monocytogenes can share habitats, they appear to differ in their histories of horizontal gene transfer.     

Comparative hybridization reveals extensive genome variation in the AIDS-associated pathogen Cryptococcus neoformans

Background

Genome variability can have a profound influence on the virulence of pathogenic microbes. The availability of genome sequences for two strains of the AIDS-associated fungal pathogen Cryptococcus neoformans presented an opportunity to use comparative genome hybridization (CGH) to examine genome variability between strains of different mating type, molecular subtype, and ploidy.

Results

Initially, CGH was used to compare the approximately 100 kilobase MAT a and MATα mating-type regions in serotype A and D strains to establish the relationship between the Log2 ratios of hybridization signals and sequence identity. Subsequently, we compared the genomes of the environmental isolate NIH433 (MAT a) and the clinical isolate NIH12 (MATα) with a tiling array of the genome of the laboratory strain JEC21 derived from these strains. In this case, CGH identified putative recombination sites and the origins of specific segments of the JEC21 genome. Similarly, CGH analysis revealed marked variability in the genomes of strains representing the VNI, VNII, and VNB molecular subtypes of the A serotype, including disomy for chromosome 13 in two strains. Additionally, CGH identified differences in chromosome content between three strains with the hybrid AD serotype and revealed that chromosome 1 from the serotype A genome is preferentially retained in all three strains.

Conclusion

The genomes of serotypes A, D, and AD strains exhibit extensive variation that spans the range from small differences (such as regions of divergence, deletion, or amplification) to the unexpected disomy for chromosome 13 in haploid strains and preferential retention of specific chromosomes in naturally occurring diploids.     

Lateral gene transfer of streptococcal ICE element RD2 (region of difference 2) encoding secreted proteins

Background  

The genome of serotype M28 group A Streptococcus (GAS) strain MGAS6180 contains a novel genetic element named Region of Difference 2 (RD2) that encodes seven putative secreted extracellular proteins. RD2 is present in all serotype M28 strains and strains of several other GAS serotypes associated with female urogenital infections. We show here that the GAS RD2 element is present in strain MGAS6180 both as an integrative chromosomal form and a circular extrachromosomal element. RD2-like regions were identified in publicly available genome sequences of strains representing three of the five major group B streptococcal serotypes causing human disease. Ten RD2-encoded proteins have significant similarity to proteins involved in conjugative transfer of Streptococcus thermophilus integrative chromosomal elements (ICEs).     

Genome‐wide patterns of transposon proliferation in an evolutionary young hybrid fish

Hybridization can induce transposons to jump into new genomic positions, which may result in their accumulation across the genome. Alternatively, transposon copy numbers may increase through nonallelic (ectopic) homologous recombination in highly repetitive regions of the genome. The relative contribution of transposition bursts versus recombination‐based mechanisms to evolutionary processes remains unclear because studies on transposon dynamics in natural systems are rare. We assessed the genomewide distribution of transposon insertions in a young hybrid lineage (“invasive Cottus”, n = 11) and its parental species Cottus rhenanus (n = 17) and Cottus perifretum(n = 9) using a reference genome assembled from long single molecule pacbio reads. An inventory of transposable elements was reconstructed from the same data and annotated. Transposon copy numbers in the hybrid lineage increased in 120 (15.9%) out of 757 transposons studied here. The copy number increased on average by 69% (range: 10%–197%). Given the age of the hybrid lineage, this suggests that they have proliferated within a few hundred generations since admixture began. However, frequency spectra of transposon insertions revealed no increase in novel and rare insertions across assembled parts of the genome. This implies that transposons were added to repetitive regions of the genome that remain difficult to assemble. Future studies will need to evaluate whether recombination‐based mechanisms rather than genomewide transposition may explain the majority of the recent transposon proliferation in the hybrid lineage. Irrespectively of the underlying mechanism, the observed overabundance in repetitive parts of the genome suggests that gene‐rich regions are unlikely to be directly affected.     

Structural variability of Tvv1 grapevine retrotransposons can be caused by illegitimate recombination

Structural variability of Tvv1, a grapevine retrotransposon Ty1 copia-like family, was investigated within the grape genome and the canonical sequence of Tvv1 determined. Then, two remarkable elements, Tvv1-Δ3001 and Tvv1-Δ3640, which had suffered large deletions 3,001 bp and 3,460 bp in length of their coding sequences were compared to the canonical copy. In both deleted elements, the deletion breakpoint was characterized by a stretch 13 bp-long in Tvv1-Δ3001 and 11 bp-long in Tvv1-Δ3640 found duplicated in the canonical copy at each bound of the deleted regions. Tvv1-Δ3001 and Tvv1-Δ3460 were both shown to be unique copies fixed at a single locus in the grapevine genome. Their presence was very variable in a set of 58 varieties and wild vines. These elements have most likely been dispersed through natural intermixing after their initial insertion whose chronology was estimated. The model that we propose to explain the structure of Tvv1-Δ3001 and Tvv1-Δ3640, implies illegitimate recombination involving template switching between two RNA molecules co-packaged in the VLP prior to the integration of the deleted daughter copy into the host genome.     

Complete genome sequence of Streptococcus troglodytae TKU31 isolated from the oral cavity of a chimpanzee (Pan troglodytes)

Streptococcus troglodytae TKU31 was isolated from the oral cavity of a chimpanzee (Pan troglodytes) and was found to be the most closely related species of the mutans group streptococci to Streptococcus mutans. The complete sequence of TKU31 genome consists of a single circular chromosome that is 2,097,874 base pairs long and has a G + C content of 37.18%. It possesses 2082 coding sequences (CDSs), 65 tRNAs and five rRNA operons (15 rRNAs). Two clustered regularly interspaced short palindromic repeats, six insertion sequences and two predicted prophage elements were identified. The genome of TKU31 harbors some putative virulence associated genes, including gtfB, gtfC and gtfD genes encoding glucosyltransferase and gbpA, gbpB, gbpC and gbpD genes encoding glucan‐binding cell wall‐anchored protein. The deduced amino acid identity of the rhamnose‐glucose polysaccharide F gene (rgpF), which is one of the serotype determinants, is 91% identical with that of S. mutans LJ23 (serotype k) strain. However, two other virulence‐associated genes cnm and cbm, which encode the collagen‐binding proteins, were not found in the TKU31 genome. The complete genome sequence of S. troglodytae TKU31 has been deposited at DDBJ/European Nucleotide Archive/GenBank under the accession no. AP014612.     

Reductive genome evolution in chemoautotrophic intracellular symbionts of deep-sea <Emphasis Type="Italic">Calyptogena</Emphasis> clams

To understand reductive genome evolution (RGE), we comparatively analyzed the recently reported small genomes of two chemoautotrophic, intracellular symbionts of deep-sea clams, Calyptogena okutanii and C. magnifica. Both genomes lack most genes for DNA recombination and repair such as recA and mutY. Their genome architectures were highly conserved except one inversion. Many deletions from small (<100 bp) to large (1–11 kbp) sizes were detected and the deletion numbers decreased exponentially with size. Densities of deletions and short-repeats, as well as A+T content were higher in non-coding regions than in coding regions. Because Calyptogena symbiont genomes lack recA, we propose that deletions and the single inversion occurred by RecA-independent recombination (RIR) at short-repeats with simultaneous consumption of repeats, and that short-repeats were regenerated by accelerated mutations with enhanced A+T bias due to the absence of mutY. We further propose that extant Calyptogena symbiont genomes are in an actively reducing stage of RGE consisting of small and large deletions, and the deletions are caused by short-repeat dependent RIR along with regeneration of short-repeats. In future, the RGE rate will slowdown when the gene repertoires approach the minimum gene set necessary for intracellular symbiotic life. Electronic supplementary material The online version of this article (doi:) contains supplementary material, which is available to authorized users. H. Kuwahara and Y. Takaki equally contributed to this work.     

Comparative genetic maps reveal extreme crossover localization in the Aegilops speltoides chromosomes

A total of 137 loci were mapped in Aegilops speltoides, the closest extant relative of the wheat B genome, using two F₂ mapping populations and a set of wheat-Ae. speltoides disomic addition (DA) lines. Comparisons of Ae. speltoides genetic maps with those of Triticum monococcum indicated that Ae. speltoides conserved the gross chromosome structure observed across the tribe Triticeae. A putative inversion involving the short arm of chromosome 2 was detected in Ae. speltoides. A translocation between chromosomes 2 and 6, present in the wheat B genome, was absent. The ligustica/aucheri spike dimorphism behaved as allelic variation at a single locus, which was mapped in the centromeric region of chromosome 3. The genetic length of each chromosome arm was about 50 cM, irrespective of its physical length. Compared to T. monococcum genetic maps, recombination was virtually eliminated from the proximal 50–100 cM and was localized in short distal regions, which were often expanded compared to the T. monococcum maps. The wheat B genome and the genome of Ae. longissima, a close relative of Ae. speltoides, do not show the extreme localization of crossovers observed in Ae. speltoides.