- Volume 7, Issue 7, 2021

Natural products (NPs) are synthesized by biosynthetic gene clusters (BGCs), whose genes are involved in producing one or a family of chemically related metabolites. Advances in comparative genomics have been favourable for exploiting huge amounts of data and discovering previously unknown BGCs. Nonetheless, studying distribution patterns of novel BGCs and elucidating the biosynthesis of orphan metabolites remains a challenge. To fill this knowledge gap, our study developed a pipeline for high-quality comparative genomics for the actinomycete genus   Rhodococcus  , which is metabolically versatile, yet understudied in terms of NPs, leading to a total of 110 genomes, 1891 BGCs and 717 non-ribosomal peptide synthetases (NRPSs). Phylogenomic inferences showed four major clades retrieved from strains of several ecological habitats. BiG-SCAPE sequence similarity BGC networking revealed 44 unidentified gene cluster families (GCFs) for NRPS, which presented a phylogenomic-dependent evolution pattern, supporting the hypothesis of vertical gene transfer. As a proof of concept, we analysed in-depth one of our marine strains,   Rhodococcus   sp. H-CA8f, which revealed a unique BGC distribution within its phylogenomic clade, involved in producing a chloramphenicol-related compound. While this BGC is part of the most abundant and widely distributed NRPS GCF, corason analysis unveiled major differences regarding its genetic context, co-occurrence patterns and modularity. This BGC is composed of three sections, two well-conserved right/left arms flanking a very variable middle section, composed of nrps genes. The presence of two non-canonical domains in H-CA8f’s BGC may contribute to adding chemical diversity to this family of NPs. Liquid chromatography-high resolution MS and dereplication efforts retrieved a set of related orphan metabolites, the corynecins, which to our knowledge are reported here for the first time in   Rhodococcus  . Overall, our data provide insights to connect BGC uniqueness with orphan metabolites, by revealing key comparative genomic features supported by models of BGC distribution along phylogeny.

Hadal trenches are the deepest but underexplored ecosystems on the Earth. Inhabiting the trench bottom is a group of micro-organisms termed obligate piezophiles that grow exclusively under high hydrostatic pressures (HHP). To reveal the genetic and physiological characteristics of their peculiar lifestyles and microbial adaptation to extreme high pressures, we sequenced the complete genome of the obligately piezophilic bacterium Moritella yayanosii DB21MT-5 isolated from the deepest oceanic sediment at the Challenger Deep, Mariana Trench. Through comparative analysis against pressure sensitive and deep-sea piezophilic Moritella strains, we identified over a hundred genes that present exclusively in hadal strain DB21MT-5. The hadal strain encodes fewer signal transduction proteins and secreted polysaccharases, but has more abundant metal ion transporters and the potential to utilize plant-derived saccharides. Instead of producing osmolyte betaine from choline as other Moritella strains, strain DB21MT-5 ferments on choline within a dedicated bacterial microcompartment organelle. Furthermore, the defence systems possessed by DB21MT-5 are distinct from other Moritella strains but resemble those in obligate piezophiles obtained from the same geographical setting. Collectively, the intensive comparative genomic analysis of an obligately piezophilic strain Moritella yayanosii DB21MT-5 demonstrates a depth-dependent distribution of energy metabolic pathways, compartmentalization of important metabolism and use of distinct defence systems, which likely contribute to microbial adaptation to the bottom of hadal trench.

SAR11 bacteria dominate ocean surface bacterioplankton communities, and play an important role in marine carbon and nutrient cycling. The biology and ecology of SAR11 are impacted by SAR11 phages (pelagiphages) that are highly diverse and abundant in the ocean. Among the currently known pelagiphages, HTVC010P represents an extremely abundant but under-studied phage group in the ocean. In this study, we have isolated seven new HTVC010P-type pelagiphages, and recovered 77 nearly full-length HTVC010P-type metagenomic viral genomes from marine metagenomes. Comparative genomic and phylogenomic analyses showed that HTVC010P-type pelagiphages display genome synteny and can be clustered into two major subgroups, with subgroup I consisting of strictly lytic phages and subgroup II mostly consisting of phages with potential lysogenic life cycles. All but one member of the subgroup II contain an integrase gene. Site-specific integration of subgroup II HTVC010P-type pelagiphage was either verified experimentally or identified by in silico genomic sequence analyses, which revealed that various SAR11 tRNA genes can serve as the integration sites of HTVC010P-type pelagiphages. Moreover, HTVC010P-type pelagiphage integration was confirmed by the detection of several Global Ocean Survey (GOS) fragments that contain hybrid phage–host integration sites. Metagenomic recruitment analysis revealed that these HTVC010P-type phages were globally distributed and most lytic subgroup I members exhibited higher relative abundance. Altogether, this study significantly expands our knowledge about the genetic diversity, life strategies and ecology of HTVC010P-type pelagiphages.

  Stenotrophomonas maltophilia   is a typical example of an environmental originated opportunistic human pathogen, which can thrive at different habitats including the human body and can cause a wide range of infections. It must cope with heat stress during transition from the environment to the human body as the physiological temperature of the human body (37 °C) is higher than environmental niches (22–30 °C). Interestingly,   S. rhizophila   a phylogenetic neighbour of   S. maltophilia   within genus   Stenotrophomonas   is unable to grow at 37 °C. Thus, it is crucial to understand how   S. maltophilia   is adapted to human body temperature, which could suggest its evolution as an opportunistic human pathogen. In this study, we have performed comparative transcriptome analysis of   S. maltophilia   grown at 28 and 37 °C as temperature representative for environmental niches and the human body, respectively. RNA-Seq analysis revealed several interesting findings showing alterations in gene-expression levels at 28 and 37 °C, which can play an important role during infection. We have observed downregulation of genes involved in cellular motility, energy production and metabolism, replication and repair whereas upregulation of VirB/D4 type IV secretion system, aerotaxis, cation diffusion facilitator family transporter and LacI family transcriptional regulators at 37 °C. Microscopy and plate assays corroborated altered expression of genes involved in motility. The results obtained enhance our understanding of the strategies employed by   S. maltophilia   during adaptation towards the human body.

Antibiotic resistance is regarded as one of the most serious threats to human health worldwide. The rapid increase in resistance rates has been attributed to the extensive use of antibiotics since they became commercially available. The use of antibiotics as growth promotors has been banned in numerous regions for this reason. Mannan-rich fraction (MRF) has been reported to show similar growth-promoting effects to antibiotics. We investigated the effect of MRF on the microbial community, resistome and metabolic pathways within the caecum of commercial broilers at two different timepoints within the growth of the broiler, day 27 and day 34. The data indicated an overall increase in health and economic gain for the producer with the addition of MRF to the diet of the broilers. The only significant difference across the microbial composition of the samples was in the richness of the microbial communities across all samples. While all samples harboured resistance genes conferring resistance to the same classes of antibiotics, there was significant variation in the antimicrobial resistance gene richness across time and treatment and across combinations of time and treatment. The taxa with positive correlation comprised Bacilli and Clostridia. The negative correlation taxa were also dominated by Bacilli, specifically the Streptococcus genera. The KEGG-pathway analysis identified an age-related change in the metabolism pathway abundances of the caecal microflora. We suggest that the MRF-related increases in health and weight gain in the broilers may be associated with changes in the metabolism of the microbiomes rather than the microbial composition. The resistome variations across samples were correlated with specific genera. These data may be used to further enhance the development of feed supplements to reduce the presence of antibiotic resistance genes (ARGs) within poultry. While the ARGs of greatest concern to human or animal health were not detected in this study, it has identified the potential to reduce the presence of ARGs by the increase in specific genera.

Arsenic is highly toxic element to all forms of life and is a major environmental contaminant. Understanding acquisition, detoxification and adaptation mechanisms in bacteria that are associated with the host in arsenic-rich conditions can provide novel insights into the evolutionary dynamics of host–microbe–environment interactions. In the present study, we have investigated an arsenic-resistance mechanism acquired during the evolution of a particular lineage in the population of   Xanthomonas oryzae   pv. oryzae, which is a serious plant pathogen infecting rice. Our study revealed the horizontal acquisition of a novel chromosomal 12 kb ars cassette in   X. oryzae   pv. oryzae IXO1088 that confers high resistance to arsenate/arsenite. The ars cassette comprises several genes that constitute an operon induced in the presence of arsenate/arsenite. Transfer of the cloned ars cassette to   X. oryzae   pv. oryzae BXO512, which lacks the cassette, confers an arsenic-resistance phenotype. Furthermore, the transcriptional response of   X. oryzae   pv. oryzae IXO1088 under arsenate/arsenite exposure was analysed using RNA sequencing. Arsenic detoxification and efflux, oxidative stress, iron acquisition/storage, and damage repair are the main cellular responses to arsenic exposure. Our investigation has provided insights into the existence of a novel detoxification and adaptation mechanism within the   X. oryzae   pv. oryzae population to deal with high-arsenic conditions outside the rice plant.

Members of the   Bifidobacterium longum   species have been shown to possess adaptive abilities to allow colonization of different mammalian hosts, including humans, primates and domesticated mammalian species, such as dogs, horses, cattle and pigs. To date, three subspecies have formally been recognized to belong to this bifidobacterial taxon, i.e.   B. longum   subsp.   longum  ,   B. longum   subsp.   infantis   and   B. longum   subsp.   suis  . Although   B. longum   subsp.   longum   is widely distributed in the human gut irrespective of host age,   B. longum   subsp.   infantis   appears to play a significant role as a prominent member of the gut microbiota of breast-fed infants. Nevertheless, despite the considerable scientific relevance of these taxa and the vast body of genomic data now available, an accurate dissection of the genetic features that comprehensively characterize the   B. longum   species and its subspecies is still missing. In the current study, we employed 261 publicly available   B. longum   genome sequences, combined with those of 11 new isolates, to investigate genomic diversity of this taxon through comparative genomic and phylogenomic approaches. These analyses allowed us to highlight a remarkable intra-species genetic and physiological diversity. Notably, characterization of the genome content of members of   B. longum   subsp.   infantis   subspecies suggested that this taxon may have acquired genetic features for increased competitiveness in the gut environment of suckling hosts. Furthermore, specific   B. longum   subsp.   infantis   genomic features appear to be responsible for enhanced horizontal gene transfer (HGT) occurrences, underpinning an intriguing dedication toward acquisition of foreign DNA by HGT events.

Pathogens of the Mycobacterium tuberculosis complex (MTBC) are considered to be monomorphic, with little gene content variation between strains. Nevertheless, several genotypic and phenotypic factors separate strains of the different MTBC lineages (L), especially L5 and L6 (traditionally termed Mycobacterium africanum ) strains, from each other. However, this genome variability and gene content, especially of L5 strains, has not been fully explored and may be important for pathobiology and current approaches for genomic analysis of MTBC strains, including transmission studies. By comparing the genomes of 355 L5 clinical strains (including 3 complete genomes and 352 Illumina whole-genome sequenced isolates) to each other and to H37Rv, we identified multiple genes that were differentially present or absent between H37Rv and L5 strains. Additionally, considerable gene content variability was found across L5 strains, including a split in the L5.3 sub-lineage into L5.3.1 and L5.3.2. These gene content differences had a small knock-on effect on transmission cluster estimation, with clustering rates influenced by the selected reference genome, and with potential overestimation of recent transmission when using H37Rv as the reference genome. We conclude that full capture of the gene diversity, especially high-resolution outbreak analysis, requires a variation of the single H37Rv-centric reference genome mapping approach currently used in most whole-genome sequencing data analysis pipelines. Moreover, the high within-lineage gene content variability suggests that the pan-genome of M. tuberculosis is at least several kilobases larger than previously thought, implying that a concatenated or reference-free genome assembly (de novo) approach may be needed for particular questions.

Despite contributing to the large disease burden in West Africa, little is known about the genomic epidemiology of Streptococcus pneumoniae which cause meningitis among children under 5 years old in the region. We analysed whole-genome sequencing data from 185  S . pneumoniae isolates recovered from suspected paediatric meningitis cases as part of the World Health Organization (WHO) invasive bacterial diseases surveillance from 2010 to 2016. The phylogeny was reconstructed, accessory genome similarity was computed and antimicrobial-resistance patterns were inferred from the genome data and compared to phenotypic resistance from disc diffusion. We studied the changes in the distribution of serotypes pre- and post-pneumococcal conjugate vaccine (PCV) introduction in the Central and Western sub-regions separately. The overall distribution of non-vaccine, PCV7 (4, 6B, 9V, 14, 18C, 19F and 23F) and additional PCV13 serotypes (1, 3, 5, 6A, 19A and 7F) did not change significantly before and after PCV introduction in the Central region (Fisher's test P value 0.27) despite an increase in the proportion of non-vaccine serotypes to 40 % (n=6) in the post-PCV introduction period compared to 21.9 % (n=14). In the Western sub-region, PCV13 serotypes were more dominant among isolates from The Gambia following the introduction of PCV7, 81 % (n=17), compared to the pre-PCV period in neighbouring Senegal, 51 % (n=27). The phylogeny illustrated the diversity of strains associated with paediatric meningitis in West Africa and highlighted the existence of phylogeographical clustering, with isolates from the same sub-region clustering and sharing similar accessory genome content. Antibiotic-resistance genotypes known to confer resistance to penicillin, chloramphenicol, co-trimoxazole and tetracycline were detected across all sub-regions. However, there was no discernible trend linking the presence of resistance genotypes with the vaccine introduction period or whether the strain was a vaccine or non-vaccine serotype. Resistance genotypes appeared to be conserved within selected sub-clades of the phylogenetic tree, suggesting clonal inheritance. Our data underscore the need for continued surveillance on the emergence of non-vaccine serotypes as well as chloramphenicol and penicillin resistance, as these antibiotics are likely still being used for empirical treatment in low-resource settings. This article contains data hosted by Microreact.

Whole-genome sequencing (WGS) is becoming the de facto standard for bacterial typing and outbreak surveillance of resistant bacterial pathogens. However, interoperability for WGS of bacterial outbreaks is poorly understood. We hypothesized that harmonization of WGS for outbreak surveillance is achievable through the use of identical protocols for both data generation and data analysis. A set of 30 bacterial isolates, comprising of various species belonging to the Enterobacteriaceae family and Enterococcus genera, were selected and sequenced using the same protocol on the Illumina MiSeq platform in each individual centre. All generated sequencing data were analysed by one centre using BioNumerics (6.7.3) for (i) genotyping origin of replications and antimicrobial resistance genes, (ii) core-genome multi-locus sequence typing (cgMLST) for Escherichia coli and Klebsiella pneumoniae and whole-genome multi-locus sequencing typing (wgMLST) for all species. Additionally, a split k-mer analysis was performed to determine the number of SNPs between samples. A precision of 99.0% and an accuracy of 99.2% was achieved for genotyping. Based on cgMLST, a discrepant allele was called only in 2/27 and 3/15 comparisons between two genomes, for E. coli and K. pneumoniae, respectively. Based on wgMLST, the number of discrepant alleles ranged from 0 to 7 (average 1.6). For SNPs, this ranged from 0 to 11 SNPs (average 3.4). Furthermore, we demonstrate that using different de novo assemblers to analyse the same dataset introduces up to 150 SNPs, which surpasses most thresholds for bacterial outbreaks. This shows the importance of harmonization of data-processing surveillance of bacterial outbreaks. In summary, multi-centre WGS for bacterial surveillance is achievable, but only if protocols are harmonized.

Travel to tropical regions is associated with high risk of acquiring extended-spectrum beta-lactamase-producing Enterobacterales (ESBL-E) that are typically cleared in less than 3 months following return. The conditions leading to persistent carriage that exceeds 3 months in some travellers require investigation. Whole-genome sequencing (Illumina MiSeq) was performed on the 82 ESBL-E isolates detected upon return and 1, 2, 3, 6 and 12 months later from the stools of 11 long-term (>3 months) ESBL-E carriers following travel abroad. One to five different ESBL Escherichia coli strains were detected per traveller upon return, and this diminished to one after 3 months. Long-term carriage was due to the presence of the same ESBL E. coli strain, for more than 3 months, in 9 out of 11 travellers, belonging to epidemic sequence type complexes (STc 10, 14, 38, 69, 131 and 648). The mean carriage duration of strains belonging to phylogroups B2/D/F, associated with extra-intestinal virulence, was higher than that for commensal-associated A/B1/E phylogroups (3.5 vs 0.5 months, P=0.021). Genes encoding iron capture systems (fyuA, irp), toxins (senB, sat), adhesins (flu, daaF, afa/nfaE, pap, ecpA) and colicin (cjrA) were more often present in persistent strains than in transient ones. Single-nucleotide polymorphism (SNP) analysis in persistent strains showed a maximum divergence of eight SNPs over 12 months without signs of adaptation. Genomic plasticity was observed during the follow-up with the loss or gain of mobile genetic elements such as plasmids, integrons and/or transposons that may contain resistance genes at different points in the follow-up. Long-term colonization of ESBL-E following travel is primarily due to the acquisition of E. coli strains belonging to epidemic clones and harbouring ‘virulence genes’, allowing good adaptation to the intestinal microbiota.

Achromobacter spp. are emerging pathogens in patients with cystic fibrosis (CF) and Achromobacter spp. caused infections are associated with more severe disease outcomes and high intrinsic antibiotic resistance. While conventional CF pathogens are studied extensively, little is known about the genetic determinants leading to antibiotic resistance and the genetic adaptation in Achromobacter spp. infections. Here, we analysed 101 Achromobacter spp. genomes from 51 patients with CF isolated during the course of up to 20 years of infection to identify within-host adaptation, mutational signatures and genetic variation associated with increased antibiotic resistance. We found that the same regulatory and inorganic ion transport genes were frequently mutated in persisting clone types within and between Achromobacter species, indicating convergent genetic adaptation. Genome-wide association study of six antibiotic resistance phenotypes revealed the enrichment of associated genes involved in inorganic ion transport, transcription gene enrichment in β-lactams, and energy production and translation gene enrichment in the trimethoprim/sulfonamide group. Overall, we provide insights into the pathogenomics of Achromobacter spp. infections in patients with CF airways. Since emerging pathogens are increasingly recognized as an important healthcare issue, our findings on evolution of antibiotic resistance and genetic adaptation can facilitate better understanding of disease progression and how mutational changes have implications for patients with CF.

Pseudomonas is a highly diverse genus that includes species that cause disease in both plants and animals. Recently, pathogenic pseudomonads from the Pseudomonas syringae and Pseudomonas fluorescens species complexes have caused significant outbreaks in several agronomically important crops in Turkey, including tomato, citrus, artichoke and melon. We characterized 169 pathogenic Pseudomonas strains associated with recent outbreaks in Turkey via multilocus sequence analysis and whole-genome sequencing, then used comparative and evolutionary genomics to characterize putative virulence mechanisms. Most of the isolates are closely related to other plant pathogens distributed among the primary phylogroups of P. syringae , although there are significant numbers of P. fluorescens isolates, which is a species better known as a rhizosphere-inhabiting plant-growth promoter. We found that all 39 citrus blast pathogens cluster in P. syringae phylogroup 2, although strains isolated from the same host do not cluster monophyletically, with lemon, mandarin orange and sweet orange isolates all being intermixed throughout the phylogroup. In contrast, 20 tomato pith pathogens are found in two independent lineages: one in the P. syringae secondary phylogroups, and the other from the P. fluorescens species complex. These divergent pith necrosis strains lack characteristic virulence factors like the canonical tripartite type III secretion system, large effector repertoires and the ability to synthesize multiple bacterial phytotoxins, suggesting they have alternative molecular mechanisms to cause disease. These findings highlight the complex nature of host specificity among plant pathogenic pseudomonads.

A key characteristic of Staphylococcus aureus infections, and one that also varies phenotypically between clones, is that of biofilm formation, which aids in bacterial persistence through increased adherence and immune evasion. Though there is a general understanding of the process of biofilm formation – adhesion, proliferation, maturation and dispersal – the tightly orchestrated molecular events behind each stage, and what drives variation between S. aureus strains, has yet to be unravelled. Herein we measure biofilm progression and dispersal in real-time across the five major S. aureus CDC-types (USA100-USA500) revealing adherence patterns that differ markedly amongst strains. To gain insight into this, we performed transcriptomic profiling on these isolates at multiple timepoints, compared to planktonically growing counterparts. Our findings support a model in which eDNA release, followed by increased positive surface charge, perhaps drives initial abiotic attachment. This is seemingly followed by cooperative repression of autolysis and activation of poly-N-acetylglucosamine (PNAG) production, which may indicate a developmental shift in structuring the biofilm matrix. As biofilms mature, diminished translational capacity was apparent, with 53 % of all ribosomal proteins downregulated, followed by upregulation of anaerobic respiration enzymes. These findings are noteworthy because reduced cellular activity and an altered metabolic state have been previously shown to contribute to higher antibiotic tolerance and bacterial persistence. In sum, this work is, to our knowledge, the first study to investigate transcriptional regulation during the early, establishing phase of biofilm formation, and to compare global transcriptional regulation both temporally and across multiple clonal lineages.

This study investigated the evolution and epidemiology of the community-associated and multidrug-resistant Staphylococcus aureus clone European CC1-MRSA-IV. Whole-genome sequences were obtained for 194 European CC1-MRSA-IV isolates (189 of human and 5 of animal origin) from 12 countries, and 10 meticillin-susceptible precursors (from North-Eastern Romania; all of human origin) of the clone. Phylogenetic analysis was performed using a maximum-likelihood approach, a time-measured phylogeny was reconstructed using Bayesian analysis, and in silico microarray genotyping was performed to identify resistance, virulence-associated and SCCmec (staphylococcal cassette chromosome mec) genes. Isolates were typically sequence type 1 (190/204) and spa type t127 (183/204). Bayesian analysis indicated that European CC1-MRSA-IV emerged in approximately 1995 before undergoing rapid expansion in the late 1990s and 2000s, while spreading throughout Europe and into the Middle East. Phylogenetic analysis revealed an unstructured meticillin-resistant S. aureus (MRSA) population, lacking significant geographical or temporal clusters. The MRSA were genotypically multidrug-resistant, consistently encoded seh, and intermittently (34/194) encoded an undisrupted hlb gene with concomitant absence of the lysogenic phage-encoded genes sak and scn. All MRSA also harboured a characteristic ~5350 nt insertion in SCCmec adjacent to orfX. Detailed demographic data from Denmark showed that there, the clone is typically (25/35) found in the community, and often (10/35) among individuals with links to South-Eastern Europe. This study elucidated the evolution and epidemiology of European CC1-MRSA-IV, which emerged from a meticillin-susceptible lineage prevalent in North-Eastern Romania before disseminating rapidly throughout Europe.

Salmonella enterica serovar Enteritidis is a major cause of foodborne Salmonella infections and outbreaks in humans. Effective surveillance and timely outbreak detection are essential for public health control. Multilevel genome typing (MGT) with multiple levels of resolution has been previously demonstrated as a promising tool for this purpose. In this study, we developed MGT with nine levels for S. Enteritidis and characterised the genomic epidemiology of S. Enteritidis in detail. We examined 26 670 publicly available S. Enteritidis genome sequences from isolates spanning 101 years from 86 countries to reveal their spatial and temporal distributions. Using the lower resolution MGT levels, globally prevalent and regionally restricted sequence types (STs) were identified; avian associated MGT4-STs were found that were common in human cases in the USA; temporal trends were observed in the UK with MGT5-STs from 2014 to 2018 revealing both long lived endemic STs and the rapid expansion of new STs. Using MGT3 to MGT6, we identified multidrug resistance (MDR) associated STs at various MGT levels, which improves precision of detection and global tracking of MDR clones. We also found that the majority of the global S. Enteritidis population fell within two predominant lineages, which had significantly different propensity of causing large scale outbreaks. An online open MGT database has been established for unified international surveillance of S. Enteritidis. We demonstrated that MGT provides a flexible and high-resolution genome typing tool for S. Enteritidis surveillance and outbreak detection.

Ascochyta blight disease, caused by the necrotrophic fungus Ascochyta rabiei, is a major biotic constraint to chickpea production in Australia and worldwide. Detailed knowledge of the structure of the pathogen population and its potential to adapt to our farming practices is key to informing optimal management of the disease. This includes understanding the molecular diversity among isolates and the frequency and distribution of the isolates that have adapted to overcome host resistance across agroecologically distinct regions. Thanks to continuous monitoring efforts over the past 6 years, a comprehensive collection of A. rabiei isolates was collated from the major Australian chickpea production regions. To determine the molecular structure of the entire population, representative isolates from each collection year and growing region have been genetically characterized using a DArTseq genotyping-by-sequencing approach. The genotyped isolates were further phenotyped to determine their pathogenicity levels against a differential set of chickpea cultivars and genotype-phenotype associations were inferred. Overall, the Australian A. rabiei population displayed a far lower genetic diversity (average Nei’s gene diversity of 0.047) than detected in other populations worldwide. This may be explained by the presence of a single mating-type in Australia, MAT1-2, limiting its reproduction to a clonal mode. Despite the low detected molecular diversity, clonal selection appears to have given rise to a subset of adapted isolates that are highly pathogenic on commonly employed resistance sources, and that are occurring at an increasing frequency. Among these, a cluster of genetically similar isolates was identified, with a higher proportion of highly aggressive isolates than in the general population. The discovery of distinct genetic clusters associated with high and low isolate pathogenicity forms the foundation for the development of a molecular pathotyping tool for the Australian A. rabiei population. Application of such a tool, along with continuous monitoring of the genetic structure of the population will provide crucial information for the screening of breeding material and integrated disease management packages.

Achromobacter species are increasingly being detected in cystic fibrosis (CF) patients, where they can establish chronic infections by adapting to the lower airway environment. To better understand the mechanisms contributing to a successful colonization by Achromobacter species, we sequenced the whole genome of 54 isolates from 26 patients with occasional and early/late chronic lung infection. We performed a phylogenetic analysis and compared virulence and resistance genes, genetic variants and mutations, and hypermutability mechanisms between chronic and occasional isolates. We identified five Achromobacter species as well as two non-affiliated genogroups (NGs). Among them were the frequently isolated Achromobacter xylosoxidans and four other species whose clinical importance is not yet clear: Achromobacter insuavis, Achromobacter dolens, Achromobacter insolitus and Achromobacter aegrifaciens . While A. insuavis and A. dolens were isolated only from chronically infected patients and A. aegrifaciens only from occasionally infected patients, the other species were found in both groups. Most of the occasional isolates lacked functional genes involved in invasiveness, chemotaxis, type 3 secretion system and anaerobic growth, whereas the great majority (>60%) of chronic isolates had these genomic features. Interestingly, almost all (n=22/23) late chronic isolates lacked functional genes involved in lipopolysaccharide production. Regarding antibiotic resistance, we observed a species-specific distribution of blaOXA genes, confirming what has been reported in the literature and additionally identifying blaOXA-2 in some A. insolitus isolates and observing no blaOXA genes in A. aegrifaciens or NGs. No significant difference in resistance genes was found between chronic and occasional isolates. The results of the mutator genes analysis showed that no occasional isolate had hypermutator characteristics, while 60% of early chronic (<1 year from first colonization) and 78% of late chronic (>1 year from first colonization) isolates were classified as hypermutators. Although all A. dolens, A. insuavis and NG isolates presented two different mutS genes, these seem to have a complementary rather than compensatory function. In conclusion, our results show that Achromobacter species can exhibit different adaptive mechanisms and some of these mechanisms might be more useful than others in establishing a chronic infection in CF patients, highlighting their importance for the clinical setting and the need for further studies on the less clinically characterized Achromobacter species.