Bioinformatic analysis of bacteria and host cell dual RNA-sequencing experiments.

Bacterial pathogens subvert host cells by manipulating cellular pathways for survival and replication; in turn, host cells respond to the invading pathogen through cascading changes in gene expression. Deciphering these complex temporal and spatial dynamics to identify novel bacterial virulence factors or host response pathways is crucial for improved diagnostics and therapeutics. Dual RNA sequencing (dRNA-Seq) has recently been developed to simultaneously capture host and bacterial transcriptomes from an infected cell. This approach builds on the high sensitivity and resolution of RNA sequencing technology and is applicable to any bacteria that interact with eukaryotic cells, encompassing parasitic, commensal or mutualistic lifestyles. Several laboratory protocols have been presented that outline the collection, extraction and sequencing of total RNA for dRNA-Seq experiments, but there is relatively little guidance available for the detailed bioinformatic analyses required. This protocol outlines a typical dRNA-Seq experiment, based on a Chlamydia trachomatis-infected host cell, with a detailed description of the necessary bioinformatic analyses with currently available software tools.

Comparative genomic analysis of human Chlamydia pneumoniae isolates from respiratory, brain and cardiac tissues.

Chlamydia pneumoniae is an obligate intracellular bacterium implicated in a wide range of human diseases including atherosclerosis and Alzheimer's disease. Efforts to understand the relationships between C. pneumoniae detected in these diseases have been hindered by the availability of sequence data for non-respiratory strains. In this study, we sequenced the whole genomes for C. pneumoniae isolates from atherosclerosis and Alzheimer's disease, and compared these to previously published C. pneumoniae genomes. Phylogenetic analyses of these new C. pneumoniae strains indicate two sub-groups within human C. pneumoniae, and suggest that both recombination and mutation events have driven the evolution of human C. pneumoniae. Further fine-detailed analyses of these new C. pneumoniae sequences show several genetically variable loci. This suggests that similar strains of C. pneumoniae are found in the brain, lungs and cardiovascular system and that only minor genetic differences may contribute to the adaptation of particular strains in human disease.

Genetic diversity in the plasticity zone and the presence of the chlamydial plasmid differentiates Chlamydia pecorum strains from pigs, sheep, cattle, and koalas.

BACKGROUND: Chlamydia pecorum is a globally recognised pathogen of livestock and koalas. To date, comparative genomics of C. pecorum strains from sheep, cattle and koalas has revealed that only single nucleotide polymorphisms (SNPs) and a limited number of pseudogenes appear to contribute to the genetic diversity of this pathogen. No chlamydial plasmid has been detected in these strains despite its ubiquitous presence in almost all other chlamydial species. Genomic analyses have not previously included C. pecorum from porcine hosts. We sequenced the genome of three C. pecorum isolates from pigs with differing pathologies in order to re-evaluate the genetic differences and to update the phylogenetic relationships between C. pecorum from each of the hosts. METHODS: Whole genome sequences for the three porcine C. pecorum isolates (L1, L17 and L71) were acquired using C. pecorum-specific sequence capture probes with culture-independent methods, and assembled in CLC Genomics Workbench. The pairwise comparative genomic analyses of 16 pig, sheep, cattle and koala C. pecorum genomes were performed using several bioinformatics platforms, while the phylogenetic analyses of the core C. pecorum genomes were performed with predicted recombination regions removed. Following the detection of a C. pecorum plasmid, a newly developed C. pecorum-specific plasmid PCR screening assay was used to evaluate the plasmid distribution in 227 C. pecorum samples from pig, sheep, cattle and koala hosts. RESULTS: Three porcine C. pecorum genomes were sequenced using C. pecorum-specific sequence capture probes with culture-independent methods. Comparative genomics of the newly sequenced porcine C. pecorum genomes revealed an increased average number of SNP differences (~11 500) between porcine and sheep, cattle, and koala C. pecorum strains, compared to previous C. pecorum genome analyses. We also identified a third copy of the chlamydial cytotoxin gene, found only in porcine C. pecorum isolates. Phylogenetic analyses clustered porcine isolates into a distinct clade, highlighting the polyphyletic origin of C. pecorum in livestock. Most surprising, we also discovered a plasmid in the porcine C. pecorum genome. Using this novel C. pecorum plasmid (pCpec) sequence, a) we developed a pCpec screening assay to evaluate the plasmid distribution in C. pecorum from different hosts; and b) to characterise the pCpec sequences from available previously sequenced C. pecorum genome data. pCpec screening showed that the pCpec is common in all hosts of C. pecorum, however not all C. pecorum strains carry pCpec. CONCLUSIONS: This study provides further insight into the complexity of C. pecorum epidemiology and novel genomic regions that may be linked to host specificity. C. pecorum plasmid characterisation may aid in improving our understanding of C. pecorum pathogenesis across the variety of host species this animal pathogen infects.

Culture-independent genome sequencing of clinical samples reveals an unexpected heterogeneity of infections by Chlamydia pecorum.

Chlamydia pecorum is an important global pathogen of livestock, and it is also a significant threat to the long-term survival of Australia's koala populations. This study employed a culture-independent DNA capture approach to sequence C. pecorum genomes directly from clinical swab samples collected from koalas with chlamydial disease as well as from sheep with arthritis and conjunctivitis. Investigations into single-nucleotide polymorphisms within each of the swab samples revealed that a portion of the reads in each sample belonged to separate C. pecorum strains, suggesting that all of the clinical samples analyzed contained mixed populations of genetically distinct C. pecorum isolates. This observation was independent of the anatomical site sampled and the host species. Using the genomes of strains identified in each of these samples, whole-genome phylogenetic analysis revealed that a clade containing a bovine and a koala isolate is distinct from other clades comprised of livestock or koala C. pecorum strains. Providing additional evidence to support exposure of koalas to Australian livestock strains, two minor strains assembled from the koala swab samples clustered with livestock strains rather than koala strains. Culture-independent probe-based genome capture and sequencing of clinical samples provides the strongest evidence yet to suggest that naturally occurring chlamydial infections are comprised of multiple genetically distinct strains.

Comparative genomics of koala, cattle and sheep strains of Chlamydia pecorum.

BACKGROUND: Chlamydia pecorum is an important pathogen of domesticated livestock including sheep, cattle and pigs. This pathogen is also a key factor in the decline of the koala in Australia. We sequenced the genomes of three koala C. pecorum strains, isolated from the urogenital tracts and conjunctiva of diseased koalas. The genome of the C. pecorum VR629 (IPA) strain, isolated from a sheep with polyarthritis, was also sequenced. RESULTS: Comparisons of the draft C. pecorum genomes against the complete genomes of livestock C. pecorum isolates revealed that these strains have a conserved gene content and order, sharing a nucleotide sequence similarity > 98%. Single nucleotide polymorphisms (SNPs) appear to be key factors in understanding the adaptive process. Two regions of the chromosome were found to be accumulating a large number of SNPs within the koala strains. These regions include the Chlamydia plasticity zone, which contains two cytotoxin genes (toxA and toxB), and a 77 kbp region that codes for putative type III effector proteins. In one koala strain (MC/MarsBar), the toxB gene was truncated by a premature stop codon but is full-length in IPTaLE and DBDeUG. Another five pseudogenes were also identified, two unique to the urogenital strains C. pecorum MC/MarsBar and C. pecorum DBDeUG, respectively, while three were unique to the koala C. pecorum conjunctival isolate IPTaLE. An examination of the distribution of these pseudogenes in C. pecorum strains from a variety of koala populations, alongside a number of sheep and cattle C. pecorum positive samples from Australian livestock, confirmed the presence of four predicted pseudogenes in koala C. pecorum clinical samples. Consistent with our genomics analyses, none of these pseudogenes were observed in the livestock C. pecorum samples examined. Interestingly, three SNPs resulting in pseudogenes identified in the IPTaLE isolate were not found in any other C. pecorum strain analysed, raising questions over the origin of these point mutations. CONCLUSIONS: The genomic data revealed that variation between C. pecorum strains were mainly due to the accumulation of SNPs, some of which cause gene inactivation. The identification of these genetic differences will provide the basis for further studies to understand the biology and evolution of this important animal pathogen.

Genome sequencing and comparative analysis of three Chlamydia pecorum strains associated with different pathogenic outcomes.

BACKGROUND: Chlamydia pecorum is the causative agent of a number of acute diseases, but most often causes persistent, subclinical infection in ruminants, swine and birds. In this study, the genome sequences of three C. pecorum strains isolated from the faeces of a sheep with inapparent enteric infection (strain W73), from the synovial fluid of a sheep with polyarthritis (strain P787) and from a cervical swab taken from a cow with metritis (strain PV3056/3) were determined using Illumina/Solexa and Roche 454 genome sequencing. RESULTS: Gene order and synteny was almost identical between C. pecorum strains and C. psittaci. Differences between C. pecorum and other chlamydiae occurred at a number of loci, including the plasticity zone, which contained a MAC/perforin domain protein, two copies of a >3400 amino acid putative cytotoxin gene and four (PV3056/3) or five (P787 and W73) genes encoding phospholipase D. Chlamydia pecorum contains an almost intact tryptophan biosynthesis operon encoding trpABCDFR and has the ability to sequester kynurenine from its host, however it lacks the genes folA, folKP and folB required for folate metabolism found in other chlamydiae. A total of 15 polymorphic membrane proteins were identified, belonging to six pmp families. Strains possess an intact type III secretion system composed of 18 structural genes and accessory proteins, however a number of putative inc effector proteins widely distributed in chlamydiae are absent from C. pecorum. Two genes encoding the hypothetical protein ORF663 and IncA contain variable numbers of repeat sequences that could be associated with persistence of infection. CONCLUSIONS: Genome sequencing of three C. pecorum strains, originating from animals with different disease manifestations, has identified differences in ORF663 and pseudogene content between strains and has identified genes and metabolic traits that may influence intracellular survival, pathogenicity and evasion of the host immune system.

Genome organization and genomics in Chlamydia: whole genome sequencing increases understanding of chlamydial virulence, evolution, and phylogeny.

The genus Chlamydia contains important obligate intracellular bacterial pathogens to humans and animals, including C. trachomatis and C. pneumoniae. Since 1998, when the first Chlamydia genome was published, our understanding of how these microbes interact, evolved and adapted to different intracellular host environments has been transformed due to the expansion of chlamydial genomes. This review explores the current state of knowledge in Chlamydia genomics and how whole genome sequencing has revolutionised our understanding of Chlamydia virulence, evolution, and phylogeny over the past two and a half decades. This review will also highlight developments in multi-omics and other approaches that have complemented whole genome sequencing to advance knowledge of Chlamydia pathogenesis and future directions for chlamydial genomics.

Phylogenomic analysis of a global collection of Escherichia coli ST38: evidence of interspecies and environmental transmission?

IMPORTANCE: Extraintestinal pathogenic Escherichia coli (ExPEC) sequence type (ST) 38 is one of the top 10 human pandemic lineages. Although a major cause of urinary tract and blood stream infections, ST38 has been poorly characterized from a global phylogenomic perspective. A comprehensive genome-scale analysis of 925 ST38 isolate genomes identified two broad ancestral clades and linkage of discrete ST38 clusters with specific bla CTX-M variants. In addition, the clades and clusters carry important virulence genes, with diverse but poorly characterized plasmids. Numerous putative interhost and environment transmission events were identified here by the presence of ST38 clones (defined as isolates with ≤35 SNPs) within humans, companion animals, food sources, urban birds, wildlife, and the environment. A small cluster of international ST38 clones from diverse sources, likely representing progenitors of a hospital outbreak that occurred in Brisbane, Australia, in 2017, was also identified. Our study emphasizes the importance of characterizing isolate genomes derived from nonhuman sources and geographical locations, without any selection bias.

Unity in variety--the pan-genome of the Chlamydiae.

Chlamydiae are evolutionarily well-separated bacteria that live exclusively within eukaryotic host cells. They include important human pathogens such as Chlamydia trachomatis as well as symbionts of protozoa. As these bacteria are experimentally challenging and genetically intractable, our knowledge about them is still limited. In this study, we obtained the genome sequences of Simkania negevensis Z, Waddlia chondrophila 2032/99, and Parachlamydia acanthamoebae UV-7. This enabled us to perform the first comprehensive comparative and phylogenomic analysis of representative members of four major families of the Chlamydiae, including the Chlamydiaceae. We identified a surprisingly large core gene set present in all genomes and a high number of diverse accessory genes in those Chlamydiae that do not primarily infect humans or animals, including a chemosensory system in P. acanthamoebae and a type IV secretion system. In S. negevensis, the type IV secretion system is encoded on a large conjugative plasmid (pSn, 132 kb). Phylogenetic analyses suggested that a plasmid similar to the S. negevensis plasmid was originally acquired by the last common ancestor of all four families and that it was subsequently reduced, integrated into the chromosome, or lost during diversification, ultimately giving rise to the extant virulence-associated plasmid of pathogenic chlamydiae. Other virulence factors, including a type III secretion system, are conserved among the Chlamydiae to variable degrees and together with differences in the composition of the cell wall reflect adaptation to different host cells including convergent evolution among the four chlamydial families. Phylogenomic analysis focusing on chlamydial proteins with homology to plant proteins provided evidence for the acquisition of 53 chlamydial genes by a plant progenitor, lending further support for the hypothesis of an early interaction between a chlamydial ancestor and the primary photosynthetic eukaryote.

Chlamydia pneumoniae is genetically diverse in animals and appears to have crossed the host barrier to humans on (at least) two occasions.

Chlamydia pneumoniae is a common human and animal pathogen associated with a wide range of diseases. Since the first isolation of C. pneumoniae TWAR in 1965, all human isolates have been essentially clonal, providing little evolutionary insight. To address this gap, we investigated the genetic diversity of 30 isolates from diverse geographical locations, from both human and animal origin (amphibian, reptilian, equine and marsupial). Based on the level of variation that we observed at 23 discreet gene loci, it was clearly evident that the animal isolates were more diverse than the isolates of human origin. Furthermore, we show that C. pneumoniae isolates could be grouped into five major genotypes, A-E, with A, B, D and E genotypes linked by geographical location, whereas genotype C was found across multiple continents. Our evidence strongly supports two separate animal-to-human cross species transfer events in the evolutionary history of this pathogen. The C. pneumoniae human genotype identified in the USA, Canada, Taiwan, Iran, Japan, Korea and Australia (non-Indigenous) most likely originated from a single amphibian or reptilian lineage, which appears to have been previously geographically widespread. We identified a separate human lineage present in two Australian Indigenous isolates (independent geographical locations). This lineage is distinct and is present in Australian amphibians as well as a range of Australian marsupials.

Genome sequence of the obligate intracellular animal pathogen Chlamydia pecorum E58.

Chlamydia pecorum is an obligate intracellular bacterial pathogen that causes diverse disease in a wide variety of economically important mammals. We report the finished complete genome sequence of C. pecorum E58, the type strain for the species.

Genome sequences of the zoonotic pathogens Chlamydia psittaci 6BC and Cal10.

Chlamydia psittaci is a highly prevalent avian pathogen and the cause of a potentially lethal zoonosis, causing life-threatening pneumonia in humans. We report the genome sequences of C. psittaci 6BC, the prototype strain of the species, and C. psittaci Cal10, a widely used laboratory strain.

Two carbapenem-resistant ST1:ST231:KL1:OCL1 Acinetobacter baumannii strains recovered in Tehran, Iran, carry AbaR31 in the chromosome and AbaR4 and TnaphA6 in a RepAci6 plasmid.

OBJECTIVES: To analyse the context of genes conferring antibiotic resistance in two carbapenem-resistant Acinetobacter baumannii isolates recovered in Tehran, Iran. METHODS: The antibiotic resistance phenotype for 28 antibiotics was determined using disc diffusion. The whole genome sequences of ABH008 and ABS200 were determined using the Illumina HiSeq X Ten platform. Resistance genes were identified using ResFinder and multilocus sequence types were determined using the Oxford and Institut Pasteur schemes. RESULTS: Isolates ABH008 and ABS200, recovered in 2012 and 2013, respectively, in two different Tehran hospitals, belong to the common global clone 1 lineage, ST1IP and ST231OX. They are resistant to sulfamethoxazole, tetracycline, gentamicin, amikacin, third-generation cephalosporins and carbapenems. Despite being isolated in different hospitals, phylogenetic analysis indicated they are closely related. Consistent with this, both isolates carry catA1, sul1, aacC1 and aadA1 in a novel variant of the AbaR3-type resistance island, named AbaR31. Both isolates are resistant to amikacin and carbapenems owing to aphA6 and oxa23, respectively. The oxa23 gene is located in the AbaR4 resistance island, and aphA6 in TnaphA6, and both mobile elements are in an ∼90 kbp plasmid encoding the putative RepAci6 replication initiation protein. Resistance to third-generation cephalosporins is due to the acquisition by homologous recombination of a 5 kb DNA segment that contains ISAba1-ampC from a ST623 strain. CONCLUSIONS: The resistance gene complements of ABH008 and ABS200 were found in AbaR31 and a plasmid that encodes RepAci6. The close genetic relationship of ABH008 and ABS200, despite each being recovered from different hospitals, indicates transmission between the two hospitals.

Dual RNA-seq analysis of in vitro infection multiplicity and RNA depletion methods in Chlamydia-infected epithelial cells.

Dual RNA-seq experiments examining viral and bacterial pathogens are increasing, but vary considerably in their experimental designs, such as infection rates and RNA depletion methods. Here, we have applied dual RNA-seq to Chlamydia trachomatis infected epithelial cells to examine transcriptomic responses from both organisms. We compared two time points post infection (1 and 24 h), three multiplicity of infection (MOI) ratios (0.1, 1 and 10) and two RNA depletion methods (rRNA and polyA). Capture of bacterial-specific RNA were greatest when combining rRNA and polyA depletion, and when using a higher MOI. However, under these conditions, host RNA capture was negatively impacted. Although it is tempting to use high infection rates, the implications on host cell survival, the potential reduced length of infection cycles and real world applicability should be considered. This data highlights the delicate nature of balancing host-pathogen RNA capture and will assist future transcriptomic-based studies to achieve more specific and relevant infection-related biological insights.

Phylogenomics of two ST1 antibiotic-susceptible non-clinical Acinetobacter baumannii strains reveals multiple lineages and complex evolutionary history in global clone 1.

Acinetobacter baumannii is an opportunistic pathogen that is difficult to treat due to its resistance to extreme conditions, including desiccation and antibiotics. Most strains causing outbreaks around the world belong to two main global lineages, namely global clones 1 and 2 (GC1 and GC2). Here, we used a combination of Illumina short read and MinION (Oxford Nanopore) long-read sequence data with a hybrid assembly approach to complete the genome sequence of two antibiotic-sensitive GC1 strains, Ex003 and Ax270, recovered in Lebanon from water and a rectal swab of a cat, respectively. Phylogenetic analysis of Ax270 and Ex003 with 186 publicly available GC1 genomes revealed two major clades, including five main lineages (L1-L5), and four single-isolate lineages outside of the two clades. Ax270 and Ex003, along with AB307-0294 and MRSN7213 (both predicted antibiotic-susceptible isolates) represent these individual lineages. Antibiotic resistance islands and transposons interrupting the comM gene remain important features in L1-L5, with L1 associated with the AbaR-type resistance islands, L2 with AbaR4, L3 strains containing either AbaR4 or its variants as well as Tn6022::ISAba42, and L4 and L5 associated with Tn6022 or its variants. Analysis of the capsule (KL) and outer core (OCL) polysaccharide loci further revealed a complex evolutionary history probably involving many recombination events. As more genomes become available, more GC1 lineages continue to emerge. However, genome sequence data from more diverse geographical regions are needed to draw a more accurate population structure of this globally distributed clone.

Comparison of koala LPCoLN and human strains of Chlamydia pneumoniae highlights extended genetic diversity in the species.

BACKGROUND: Chlamydia pneumoniae is a widespread pathogen causing upper and lower respiratory tract infections in addition to a range of other diseases in humans and animals. Previous whole genome analyses have focused on four essentially clonal (> 99% identity) C. pneumoniae human genomes (AR39, CWL029, J138 and TW183), providing relatively little insight into strain diversity and evolution of this species. RESULTS: We performed individual gene-by-gene comparisons of the recently sequenced C. pneumoniae koala genome and four C. pneumoniae human genomes to identify species-specific genes, and more importantly, to gain an insight into the genetic diversity and evolution of the species. We selected genes dispersed throughout the chromosome, representing genes that were specific to C. pneumoniae, genes with a demonstrated role in chlamydial biology and/or pathogenicity (n = 49), genes encoding nucleotide salvage or amino acid biosynthesis proteins (n = 6), and extrachromosomal elements (9 plasmid and 2 bacteriophage genes). CONCLUSIONS: We have identified strain-specific differences and targets for detection of C. pneumoniae isolates from both human and animal origin. Such characterisation is necessary for an improved understanding of disease transmission and intervention.

Genome sequence and identification of candidate vaccine antigens from the animal pathogen Dichelobacter nodosus.

Dichelobacter nodosus causes ovine footrot, a disease that leads to severe economic losses in the wool and meat industries. We sequenced its 1.4-Mb genome, the smallest known genome of an anaerobe. It differs markedly from small genomes of intracellular bacteria, retaining greater biosynthetic capabilities and lacking any evidence of extensive ongoing genome reduction. Comparative genomic microarray studies and bioinformatic analysis suggested that, despite its small size, almost 20% of the genome is derived from lateral gene transfer. Most of these regions seem to be associated with virulence. Metabolic reconstruction indicated unsuspected capabilities, including carbohydrate utilization, electron transfer and several aerobic pathways. Global transcriptional profiling and bioinformatic analysis enabled the prediction of virulence factors and cell surface proteins. Screening of these proteins against ovine antisera identified eight immunogenic proteins that are candidate antigens for a cross-protective vaccine.

Analysis of Complete Genome Sequence of Acinetobacter baumannii Strain ATCC 19606 Reveals Novel Mobile Genetic Elements and Novel Prophage.

Acinetobacter baumannii isolate ATCC 19606 was recovered in the US prior to 1948. It has been used as a reference and model organism in many studies involving antibiotic resistance and pathogenesis of A. baumannii, while, until recently, a complete genome of this strain was not available. Here, we present an analysis of the complete 3.91-Mbp genome sequence, generated via a combination of short-read sequencing (Illumina) and long-read sequencing (MinION), and show it contains two small cryptic plasmids and a novel complete prophage of size 41.2 kb. We also characterised several regions of the ATCC 19606 genome, leading to the identification of a novel cadmium/mercury transposon, which was named Tn6551. ATCC 19606 is an antibiotic-sensitive strain, but a comparative analysis of all publicly available ST52 strains predicts a resistance to modern antibiotics by the accumulation of antibiotic-resistance genes via plasmids in recent isolates that belong to this sequence type.

Chromatin accessibility dynamics of Chlamydia-infected epithelial cells.

Chlamydia are Gram-negative, obligate intracellular bacterial pathogens responsible for a broad spectrum of human and animal diseases. In humans, Chlamydia trachomatis is the most prevalent bacterial sexually transmitted infection worldwide and is the causative agent of trachoma (infectious blindness) in disadvantaged populations. Over the course of its developmental cycle, Chlamydia extensively remodels its intracellular niche and parasitises the host cell for nutrients, with substantial resulting changes to the host cell transcriptome and proteome. However, little information is available on the impact of chlamydial infection on the host cell epigenome and global gene regulation. Regions of open eukaryotic chromatin correspond to nucleosome-depleted regions, which in turn are associated with regulatory functions and transcription factor binding. We applied formaldehyde-assisted isolation of regulatory elements enrichment followed by sequencing (FAIRE-Seq) to generate temporal chromatin maps of C. trachomatis-infected human epithelial cells in vitro over the chlamydial developmental cycle. We detected both conserved and distinct temporal changes to genome-wide chromatin accessibility associated with C. trachomatis infection. The observed differentially accessible chromatin regions include temporally-enriched sets of transcription factors, which may help shape the host cell response to infection. These regions and motifs were linked to genomic features and genes associated with immune responses, re-direction of host cell nutrients, intracellular signalling, cell-cell adhesion, extracellular matrix, metabolism and apoptosis. This work provides another perspective to the complex response to chlamydial infection, and will inform further studies of transcriptional regulation and the epigenome in Chlamydia-infected human cells and tissues.

Genomic profiling of Escherichia coli isolates from bacteraemia patients: a 3-year cohort study of isolates collected at a Sydney teaching hospital.

This study sought to assess the genetic variability of Escherichia coli isolated from bloodstream infections (BSIs) presenting at Concord Hospital, Sydney during 2013-2016. Whole-genome sequencing was used to characterize 81 E. coli isolates sourced from community-onset (CO) and hospital-onset (HO) BSIs. The cohort comprised 64 CO and 17 HO isolates, including 35 multidrug-resistant (MDR) isolates exhibiting phenotypic resistance to three or more antibiotic classes. Phylogenetic analysis identified two major ancestral clades. One was genetically diverse with 25 isolates distributed in 16 different sequence types (STs) representing phylogroups A, B1, B2, C and F, while the other comprised phylogroup B2 isolates in subclades representing the ST131, ST73 and ST95 lineages. Forty-seven isolates contained a class 1 integron, of which 14 carried blaCTX -M-gene. Isolates with a class 1 integron carried more antibiotic resistance genes than isolates without an integron and, in most instances, resistance genes were localized within complex resistance loci (CRL). Resistance to fluoroquinolones could be attributed to point mutations in chromosomal parC and gyrB genes and, in addition, two isolates carried a plasmid-associated qnrB4 gene. Co-resistance to fluoroquinolone and broad-spectrum beta-lactam antibiotics was associated with ST131 (HO and CO), ST38 (HO), ST393 (CO), ST2003 (CO) and ST8196 (CO and HO), a novel ST identified in this study. Notably, 10/81 (12.3 %) isolates with ST95 (5 isolates), ST131 (2 isolates), ST88 (2 isolates) and a ST540 likely carry IncFII-IncFIB plasmid replicons with a full spectrum of virulence genes consistent with the carriage of ColV-like plasmids. Our data indicate that IncF plasmids play an important role in shaping virulence and resistance gene carriage in BSI E. coli in Australia.