ClonalFrameML: efficient inference of recombination in whole bacterial genomes

For Someone, You Are the Whole World: Host-Specificity of Salmonella enterica

Salmonella enterica is a bacterial pathogen known to cause gastrointestinal infections in diverse hosts, including humans and animals. Despite extensive knowledge of virulence mechanisms, understanding the factors driving host specificity remains limited. In this study, we performed a comprehensive pangenome-wide analysis of S. enterica to identify potential loci determining preference towards certain hosts. We used a dataset of high-quality genome assemblies grouped into 300 reference clusters with a special focus on four host groups: humans, pigs, cattle, and birds. The reconstructed pangenome was shown to be open and enriched with the accessory component implying high genetic diversity. Notably, phylogenetic inferences did not correspond to the distribution of affected hosts, as large compact phylogenetic groups were absent. By performing a pangenome-wide association study, we identified potential host specificity determinants. These included multiple genes encoding proteins involved in distinct infection stages, e.g., secretion systems, surface structures, transporters, transcription regulators, etc. We also identified antibiotic resistance loci in host-adapted strains. Functional annotation corroborated the results obtained with significant enrichments related to stress response, antibiotic resistance, ion transport, and surface or extracellular localization. We suggested categorizing the revealed specificity factors into three main groups: pathogenesis, resistance to antibiotics, and propagation of mobile genetic elements (MGEs).

Whole-genome sequencing analysis of molecular epidemiology and silent transmissions causing meticillin-resistant Staphylococcus aureus bloodstream infections in a university hospital

BACKGROUND

The emergence of novel genomic-type clones, such as community-associated meticillin-resistant Staphylococcus aureus (MRSA) and livestock-associated MRSA, and their invasion into hospitals have become major concerns worldwide; however, little information is available regarding the prevalence of MRSA in Japan. Whole-genome sequencing (WGS) has been conducted to analyse various pathogens worldwide. Therefore, it is important to establish a genome database of clinical MRSA isolates available in Japan.

AIM

A molecular epidemiological analysis of MRSA strains isolated from bloodstream-infected patients in a Japanese university hospital was conducted using WGS and single-nucleotide polymorphism (SNP) analysis. Additionally, through a review of patients' clinical characteristics, the effectiveness of SNP analysis as a tool for detecting silent nosocomial transmission that may be missed by other methods was evaluated in diverse settings and various time points of detection.

METHODS

Polymerase-chain-reaction-based staphylococcal cassette chromosome mec (SCCmec) typing was performed using 135 isolates obtained between 2014 and 2018, and WGS was performed using 88 isolates obtained between 2015 and 2017.

FINDINGS

SCCmec type II strains, prevalent in 2014, became rare in 2018, whereas the prevalence of SCCmec type IV strains increased from 18.75% to 83.87% of the population, and became the dominant clones. Clonal complex (CC) 5 CC8 and CC1 were detected between 2015 and 2017, with CC1 being dominant. In 88 cases, SNP analyses revealed nosocomial transmissions among 20 patients which involved highly homologous strains.

CONCLUSIONS

Routine monitoring of MRSA by whole-genome analysis is effective not only for gaining knowledge regarding molecular epidemiology, but also for detecting silent nosocomial transmission.

Ecological differences among hydrothermal vent symbioses may drive contrasting patterns of symbiont population differentiation

The intra-host composition of horizontally transmitted microbial symbionts can vary across host populations due to interactive effects of host genetics, environmental, and geographic factors. While adaptation to local habitat conditions can drive geographic subdivision of symbiont strains, it is unknown how differences in ecological characteristics among host-symbiont associations influence the genomic structure of symbiont populations. To address this question, we sequenced metagenomes of different populations of the deep-sea mussel Bathymodiolus septemdierum, which are common at Western Pacific deep-sea hydrothermal vents and show characteristic patterns of niche partitioning with sympatric gastropod symbioses. Bathymodiolus septemdierum lives in close symbiotic relationship with sulfur-oxidizing chemosynthetic bacteria but supplements its symbiotrophic diet through filter-feeding, enabling it to occupy ecological niches with little exposure to geochemical reductants. Our analyses indicate that symbiont populations associated with B. septemdierum show structuring by geographic location, but that the dominant symbiont strain is uncorrelated with vent site. These patterns are in contrast to co-occurring Alviniconcha and Ifremeria gastropod symbioses that exhibit greater symbiont nutritional dependence and occupy habitats with higher spatial variability in environmental conditions. Our results suggest that relative habitat homogeneity combined with sufficient symbiont dispersal and genomic mixing might promote persistence of similar symbiont strains across geographic locations, while mixotrophy might decrease selective pressures on the host to affiliate with locally adapted symbiont strains. Overall, these data contribute to our understanding of the potential mechanisms influencing symbiont population structure across a spectrum of marine microbial symbioses that occupy contrasting ecological niches. IMPORTANCE Beneficial relationships between animals and microbial organisms (symbionts) are ubiquitous in nature. In the ocean, microbial symbionts are typically acquired from the environment and their composition across geographic locations is often shaped by adaptation to local habitat conditions. However, it is currently unknown how generalizable these patterns are across symbiotic systems that have contrasting ecological characteristics. To address this question, we compared symbiont population structure between deep-sea hydrothermal vent mussels and co-occurring but ecologically distinct snail species. Our analyses show that mussel symbiont populations are less partitioned by geography and do not demonstrate evidence for environmental adaptation. We posit that the mussel's mixotrophic feeding mode may lower its need to affiliate with locally adapted symbiont strains, while microhabitat stability and symbiont genomic mixing likely favors persistence of symbiont strains across geographic locations. Altogether, these findings further our understanding of the mechanisms shaping symbiont population structure in marine environmentally transmitted symbioses.

Multi-host infection and phylogenetically diverse lineages shape the recombination and gene pool dynamics of Staphylococcus aureus

BACKGROUND

Staphylococcus aureus can infect and adapt to multiple host species. However, our understanding of the genetic and evolutionary drivers of its generalist lifestyle remains inadequate. This is particularly important when considering local populations of S. aureus, where close physical proximity between bacterial lineages and between host species may facilitate frequent and repeated interactions between them. Here, we aim to elucidate the genomic differences between human- and animal-derived S. aureus from 437 isolates sampled from disease cases in the northeast region of the United States.

RESULTS

Multi-locus sequence typing revealed the existence of 75 previously recognized sequence types (ST). Our population genomic analyses revealed heterogeneity in the accessory genome content of three dominant S. aureus lineages (ST5, ST8, ST30). Genes related to antimicrobial resistance, virulence, and plasmid types were differentially distributed among isolates according to host (human versus non-human) and among the three major STs. Across the entire population, we identified a total of 1,912 recombination events that occurred in 765 genes. The frequency and impact of homologous recombination were comparable between human- and animal-derived isolates. Low-frequency STs were major donors of recombined DNA, regardless of the identity of their host. The most frequently recombined genes (clfB, aroA, sraP) function in host infection and virulence, which were also frequently shared between the rare lineages.

CONCLUSIONS

Taken together, these results show that frequent but variable patterns of recombination among co-circulating S. aureus lineages, including the low-frequency lineages, that traverse host barriers shape the structure of local gene pool and the reservoir of host-associated genetic variants. Our study provides important insights to the genetic and evolutionary factors that contribute to the ability of S. aureus to colonize and cause disease in multiple host species. Our study highlights the importance of continuous surveillance of S. aureus circulating in different ecological host niches and the need to systematically sample from them. These findings will inform development of effective measures to control S. aureus colonization, infection, and transmission across the One Health continuum.

Genotypic Comparison of Pasteurella multocida from Healthy Animals at Entry to the Feedlots with That and from Bovine Respiratory Disease-Affected Animals during the Fattening Period

The aim of this study was to investigate the possible genotypic differences between commensal Pasteurella multocida isolates from apparently healthy animals (AHA) at the time of entry to feedlots and those from BRD-affected animals (BRD-AA). A total of 20 batches of beef calves in seven feedlots were followed-up during the fattening period. P. multocida was isolated from 28.1% of AHA and 22.9% of BRD-AA. All isolates belonged to the A: L3 genotype. Most isolates from clinical cases (81.0%) grouped into a PFGE cluster were significantly associated with BRD cases (OR, 24.9; 95% CI, 6.4-96.2). The whole genomes of 14 isolates representative of the pulsotypes most frequently detected in BRD-AA and AHA were sequenced and compared with 53 bovine genomes belonging to the identified ST13, ST79, and ST80 genotypes for a global comparison. No differences were found in the virulence-associated gene content between sequence types (STs) globally or between BRD-AA and AHA isolates in this study. Significantly, ST79 isolates harbored ARGs, conferring resistance to different antimicrobials, including macrolides and tetracyclines, which are commonly used for the treatment of BRD. Two Spanish ST79 isolates carried an ICE highly similar to ICE Tn7407, which was recently detected in Germany, suggesting that ST79 P. multocida isolates in Europe and North America may be associated with different ICEs.

Pangenome analysis provides insights into the genetic diversity, metabolic versatility, and evolution of the genus Flavobacterium

Members of the genus Flavobacterium are widely distributed and produce various polysaccharide-degrading enzymes. Many species in the genus have been isolated and characterized. However, few studies have focused on marine isolates or fish pathogens, and in-depth genomic analyses, particularly comparative analyses of isolates from different habitat types, are lacking. Here, we isolated 20 strains of the genus from various environments in South Korea and sequenced their full-length genomes. Combined with published sequence data, we examined genomic traits, evolution, environmental adaptation, and putative metabolic functions in total 187 genomes of isolated species in Flavobacterium categorized as marine, host-associated, and terrestrial including freshwater. A pangenome analysis revealed a correlation between genome size and coding or noncoding density. Flavobacterium spp. had high levels of diversity, allowing for novel gene repertories via recombination events. Defense-related genes only accounted for approximately 3% of predicted genes in all Flavobacterium genomes. While genes involved in metabolic pathways did not differ with respect to isolation source, there was substantial variation in genomic traits; in particular, the abundances of tRNAs and rRNAs were higher in the host-associdated group than in other groups. One genome in the host-associated group contained a Microviridae prophage closely related to an enterobacteria phage. The proteorhodopsin gene was only identified in four terrestrial strains isolated for this study. Furthermore, recombination events clearly influenced genomic diversity and may contribute to the response to environmental stress. These findings shed light on the high genetic variation in Flavobacterium and functional roles in diverse ecosystems as a result of their metabolic versatility. IMPORTANCE The genus Flavobacterium is a diverse group of bacteria that are found in a variety of environments. While most species of this genus are harmless and utilize organic substrates such as proteins and polysaccharides, some members may play a significant role in the cycling for organic substances within their environments. Nevertheless, little is known about the genomic dynamics and/or metabolic capacity of Flavobacterium. Here, we found that Flavobacterium species may have an open pangenome, containing a variety of diverse and novel gene repertoires. Intriguingly, we discovered that one genome (classified into host-associated group) contained a Microviridae prophage closely related to that of enterobacteria. Proteorhodopsin may be expressed under conditions of light or oxygen pressure in some strains isolated for this study. Our findings significantly contribute to the understanding of the members of the genus Flavobacterium diversity exploration and will provide a framework for the way for future ecological characterizations.

Portrait of a generalist bacterium: pathoadaptation, metabolic specialization and extreme environments shape diversity of Staphylococcus saprophyticus

Staphylococcus saprophyticus is a Gram-positive, coagulase-negative staphylococcus found in diverse environments including soil and freshwater, meat, and dairy foods. S. saprophyticus is also an important cause of urinary tract infections (UTIs) in humans, and mastitis in cattle. However, the genetic determinants of virulence have not yet been identified, and it remains unclear whether there are distinct sub-populations adapted to human and animal hosts. Using a diverse sample of S. saprophyticus isolates from food, animals, environmental sources, and human infections, we characterized the population structure and diversity of global populations of S. saprophyticus . We found that divergence of the two major clades of S. saprophyticus is likely facilitated by barriers to horizontal gene transfer (HGT) and differences in metabolism. Using genome-wide association study (GWAS) tools we identified the first Type VII secretion system (T7SS) described in S. saprophyticus and its association with bovine mastitis. Finally, we found that in general, strains of S. saprophyticus from different niches are genetically similar with the exception of built environments, which function as a 'sink' for S. saprophyticus populations. This work increases our understanding of the ecology of S. saprophyticus and of the genomics of bacterial generalists.

Data summary

Raw sequencing data for newly sequenced S. saprophyticus isolates have been deposited to the NCBI SRA under the project accession PRJNA928770. A list of all genomes used in this work and their associated metadata are available in the supplementary material. Custom scripts used in the comparative genomics and GWAS analyses are available here: https://github.com/myoungblom/sapro_genomics .

Impact statement

It is not known whether human and cattle diseases caused by S. saprophyticus represent spillover events from a generalist adapted to survive in a range of environments, or whether the capacity to cause disease represents a specific adaptation. Seasonal cycles of S. saprophyticus UTIs and molecular epidemiological evidence suggest that these infections may be environmentally-acquired rather than via transmission from person to person. Using comparative genomics and genome wide association study tools, we found that S. saprophyticus appears adapted to inhabit a wide range of environments (generalist), with isolates from animals, food, natural environments and human infections being closely related. Bacteria that routinely switch environments, particularly between humans and animals, are of particular concern when it comes to the spread of antibiotic resistance from farm environments into human populations. This work provides a framework for comparative genomic analyses of bacterial generalists and furthers our understanding of how bacterial populations move between humans, animals, and the environment.

Dynamics of the MRSA Population in a Chilean Hospital: a Phylogenomic Analysis (2000-2016)

The global dissemination of methicillin-resistant Staphylococcus aureus (MRSA) is associated with the emergence and establishment of clones in specific geographic areas. The Chilean-Cordobes clone (ChC) (ST5-SCCmecI) has been the predominant MRSA clone in Chile since its first description in 1998, despite the report of other emerging MRSA clones in recent years. Here, we characterize the evolutionary history of MRSA from 2000 to 2016 in a Chilean tertiary health care center using phylogenomic analyses. We sequenced 469 MRSA isolates collected between 2000 and 2016. We evaluated the temporal trends of the circulating clones and performed a phylogenomic reconstruction to characterize the clonal dynamics. We found a significant increase in the diversity and richness of sequence types (STs; Spearman r = 0.8748, P < 0.0001) with a Shannon diversity index increasing from 0.221 in the year 2000 to 1.33 in 2016, and an effective diversity (Hill number; q = 2) increasing from 1.12 to 2.71. The temporal trend analysis revealed that in the period 2000 to 2003 most of the isolates (94.2%; n = 98) belonged to the ChC clone. However, since then, the frequency of the ChC clone has decreased over time, accounting for 52% of the collection in the 2013 to 2016 period. This decline was accompanied by the rise of two emerging MRSA lineages, ST105-SCCmecII and ST72-SCCmecVI. In conclusion, the ChC clone remains the most frequent MRSA lineage, but this lineage is gradually being replaced by several emerging clones, the most important of which is clone ST105-SCCmecII. To the best of our knowledge, this is the largest study of MRSA clonal dynamics performed in South America. IMPORTANCE Methicillin-resistant Staphylococcus aureus (MRSA) is a major public health pathogen that disseminates through the emergence of successful dominant clones in specific geographic regions. Knowledge of the dissemination and molecular epidemiology of MRSA in Latin America is scarce and is largely based on small studies or more limited typing techniques that lack the resolution to represent an accurate description of the genomic landscape. We used whole-genome sequencing to study 469 MRSA isolates collected between 2000 and 2016 in Chile providing the largest and most detailed study of clonal dynamics of MRSA in South America to date. We found a significant increase in the diversity of MRSA clones circulating over the 17-year study period. Additionally, we describe the emergence of two novel clones (ST105-SCCmecII and ST72-SCCmecVI), which have been gradually increasing in frequency over time. Our results drastically improve our understanding of the dissemination and update our knowledge about MRSA in Latin America.

Characterization of the two tandem repeats for the KPC-2 core structures on a plasmid from hospital-derived Klebsiella pneumoniae

Today, Klebsiella pneumoniae strains are sophisticatedly associated with the transmission of KPC, and ST11 clones carrying KPC-2 are an important target for anti-infective clinical therapy, posing a very high threat to patients. To present the detailed genetic features of two KPC-2 core structures of F94_plasmid pA, the whole genome of K. pneumoniae strain F94 was sequenced by nanopore and illumina platform, and mobile genetic elements associated with antibiotic-resistance genes were analyzed with a series of bioinformatics methods. K. pneumoniae strain F94, identified as a class A carbapenemase-resistant Enterobacteriaceae, was resistant to most tested antibiotics, especially to low-levels of ceftazidime/avibactam (avibactam ≤ 4 mg/L), owing to overexpression of the two KPC-2 in F94_plasmid pA. However, strain F94 was sensitive to high-levels of ceftazidime/avibactam (avibactam ≥ 8 mg/L), which correlated with further inhibition of ceftazidime hydrolysis by the KPC-2 enzyme due to the multiplication of avibactam. Collinearity analysis indicated that multi-drug resistance (MDR) regions of plasmids with the tandam repeats of two or more KPC-2 core structures share highly similar structures. This study characterized the MDR region of the F94_ plasmid pA as homologous to plasmids pKPC2_090050, pKPC2_090374, plasmid unnamed 2, pC2414-2-KPC, pKPC2-020037, pBS1014-KPC2, pKPC-J5501, and pKPC2-020002, which contained the tandem repeats of one, two, or more KPC-2 core structures, providing insight into the evolution of multidrug resistance in K. pneumoniae. An alternative theoretical basis for exploring the tandem repeats of two or more KPC-2 core structures was developed by analyzing and constructing the homologous sequence of F94_ plasmid pA.

Herbarium specimen sequencing allows precise dating of Xanthomonas citri pv. citri diversification history

Herbarium collections are an important source of dated, identified and preserved DNA, whose use in comparative genomics and phylogeography can shed light on the emergence and evolutionary history of plant pathogens. Here, we reconstruct 13 historical genomes of the bacterial crop pathogen Xanthomonas citri pv. citri (Xci) from infected Citrus herbarium specimens. Following authentication based on ancient DNA damage patterns, we compare them with a large set of modern genomes to estimate their phylogenetic relationships, pathogenicity-associated gene content and several evolutionary parameters. Our results indicate that Xci originated in Southern Asia ~11,500 years ago (perhaps in relation to Neolithic climate change and the development of agriculture) and diversified during the beginning of the 13th century, after Citrus diversification and before spreading to the rest of the world (probably via human-driven expansion of citriculture through early East-West trade and colonization).

Acquisition, co-option, and duplication of the rtx toxin system and the emergence of virulence in Kingella

The bacterial genus Kingella includes two pathogenic species, namely Kingella kingae and Kingella negevensis, as well as strictly commensal species. Both K. kingae and K. negevensis secrete a toxin called RtxA that is absent in the commensal species. Here we present a phylogenomic study of the genus Kingella, including new genomic sequences for 88 clinical isolates, genotyping of another 131 global isolates, and analysis of 52 available genomes. The phylogenetic evidence supports that the toxin-encoding operon rtxCA was acquired by a common ancestor of the pathogenic Kingella species, and that a preexisting type-I secretion system was co-opted for toxin export. Subsequent genomic reorganization distributed the toxin machinery across two loci, with 30-35% of K. kingae strains containing two copies of the rtxA toxin gene. The rtxA duplication is largely clonal and is associated with invasive disease. Assays with isogenic strains show that a single copy of rtxA is associated with reduced cytotoxicity in vitro. Thus, our study identifies key steps in the evolutionary transition from commensal to pathogen, including horizontal gene transfer, co-option of an existing secretion system, and gene duplication.

Antibiotic Consumption During the Coronavirus Disease 2019 Pandemic and Emergence of Carbapenemase-Producing Klebsiella pneumoniae Lineages Among Inpatients in a Chilean Hospital: A Time-Series Study and Phylogenomic Analysis

BACKGROUND

The impact of coronavirus disease 2019 (COVID-19) on antimicrobial use (AU) and resistance has not been well evaluated in South America. These data are critical to inform national policies and clinical care.

METHODS

At a tertiary hospital in Santiago, Chile, between 2018 and 2022, subdivided into pre- (3/2018-2/2020) and post-COVID-19 onset (3/2020-2/2022), we evaluated intravenous AU and frequency of carbapenem-resistant Enterobacterales (CRE). We grouped monthly AU (defined daily doses [DDD]/1000 patient-days) into broad-spectrum β-lactams, carbapenems, and colistin and used interrupted time-series analysis to compare AU during pre- and post-pandemic onset. We studied the frequency of carbapenemase-producing (CP) CRE and performed whole-genome sequencing analyses of all carbapenem-resistant (CR) Klebsiella pneumoniae (CRKpn) isolates collected during the study period.

RESULTS

Compared with pre-pandemic, AU (DDD/1000 patient-days) significantly increased after the pandemic onset, from 78.1 to 142.5 (P < .001), 50.9 to 110.1 (P < .001), and 4.1 to 13.3 (P < .001) for broad-spectrum β-lactams, carbapenems, and colistin, respectively. The frequency of CP-CRE increased from 12.8% pre-COVID-19 to 51.9% after pandemic onset (P < .001). The most frequent CRE species in both periods was CRKpn (79.5% and 76.5%, respectively). The expansion of CP-CRE harboring blaNDM was particularly noticeable, increasing from 40% (n = 4/10) before to 73.6% (n = 39/53) after pandemic onset (P < .001). Our phylogenomic analyses revealed the emergence of two distinct genomic lineages of CP-CRKpn: ST45, harboring blaNDM, and ST1161, which carried blaKPC.

CONCLUSIONS

AU and the frequency of CP-CRE increased after COVID-19 onset. The increase in CP-CRKpn was driven by the emergence of novel genomic lineages. Our observations highlight the need to strengthen infection prevention and control and antimicrobial stewardship efforts.

Emergence of linezolid-resistant Enterococcus faecium in a tertiary hospital in Copenhagen

Linezolid is used as first-line treatment of infections caused by vancomycin-resistant Enterococcus faecium. However, resistance to linezolid is increasingly detected. The aim of the present study was to elucidate the causes and mechanisms for the increase in linezolid-resistant E. faecium at Copenhagen University Hospital - Rigshospitalet. We therefore combined patient information on linezolid treatment with whole-genome sequencing data for vancomycin- or linezolid-resistant E. faecium isolates that had been systematically collected since 2014 (n=458). Whole-genome sequencing was performed for multilocus sequence typing (MLST), identification of linezolid resistance-conferring genes/mutations and determination of phylogenetically closely related strains. The collection of E. faecium isolates belonged to prevalent vancomycin-resistant MLST types. Among these, we identified clusters of closely related linezolid-resistant strains compatible with nosocomial transmission. We also identified linezolid-resistant enterococcus isolates not genetically closely related to other isolates compatible with de novo generation of linezolid resistance. Patients with the latter isolates were significantly more frequently exposed to linezolid treatment than patients with related linezolid-resistant enterococcus isolates. We also identified six patients who initially carried a vancomycin-resistant, linezolid-sensitive enterococcus, but from whom vancomycin-resistant, linezolid-resistant enterococci (LVRE) closely related to their initial isolate were recovered after linezolid treatment. Our data illustrate that linezolid resistance may develop in the individual patient subsequent to linezolid exposure and can be transmitted between patients in a hospital setting.

Phylogeny, Virulence, and Antimicrobial Resistance Gene Profiles of Enterococcus faecium Isolated from Australian Feedlot Cattle and Their Significance to Public and Environmental Health

The extent of similarity between E. faecium strains found in healthy feedlot beef cattle and those causing extraintestinal infections in humans is not yet fully understood. This study used whole-genome sequencing to analyse the antimicrobial resistance profile of E. faecium isolated from beef cattle (n = 59) at a single feedlot and compared them to previously reported Australian isolates obtained from pig (n = 60) and meat chicken caecal samples (n = 8), as well as human sepsis cases (n = 302). The E. faecium isolated from beef cattle and other food animal sources neither carried vanA/vanB responsible for vancomycin nor possessed gyrA/parC and liaR/liaS gene mutations associated with high-level fluoroquinolone and daptomycin resistance, respectively. A small proportion (7.6%) of human isolates clustered with beef cattle and pig isolates, including a few isolates belonging to the same sequence types ST22 (one beef cattle, one pig, and two human isolates), ST32 (eight beef cattle and one human isolate), and ST327 (two beef cattle and one human isolate), suggesting common origins. This provides further evidence that these clonal lineages may have broader host range but are unrelated to the typical hospital-adapted human strains belonging to clonal complex 17, significant proportions of which contain vanA/vanB and liaR/liaS. Additionally, none of the human isolates belonging to these STs contained resistance genes to WHO critically important antimicrobials. The results confirm that most E. faecium isolated from beef cattle in this study do not pose a significant risk for resistance to critically important antimicrobials and are not associated with current human septic infections.

Human Anelloviruses: Influence of Demographic Factors, Recombination, and Worldwide Diversity

Anelloviruses represent the major and most diverse component of the healthy human virome, referred to as the anellome. In this study, we determined the anellome of 50 blood donors, forming two sex- and age-matched groups. Anelloviruses were detected in 86% of the donors. The number of detected anelloviruses increased with age and was approximately twice as high in men as in women. A total of 349 complete or nearly complete genomes were classified as belonging to torque teno virus (TTV), torque teno mini virus (TTMV), and torque teno midi virus (TTMDV) anellovirus genera (197, 88, and 64 sequences, respectively). Most donors had intergenus (69.8%) or intragenus (72.1%) coinfections. Despite the limited number of sequences, intradonor recombination analysis showed 6 intragenus recombination events in ORF1. As thousands of anellovirus sequences have been described recently, we finally analyzed the global diversity of human anelloviruses. Species richness and diversity were close to saturation in each anellovirus genus. Recombination was found to be the main factor promoting diversity, although its effect was significantly lower in TTV than in TTMV and TTMDV. Overall, our results suggest that differences in diversity between genera may be caused by variations in the relative contribution of recombination. IMPORTANCE Anelloviruses are the most common human infectious viruses and are considered essentially harmless. Compared to other human viruses, they are characterized by enormous diversity, and recombination is suggested to play an important role in their diversification and evolution. Here, by analyzing the composition of the plasma anellome of 50 blood donors, we find that recombination is also a determinant of viral evolution at the intradonor level. On a larger scale, analysis of anellovirus sequences currently available in databases shows that their diversity is close to saturation and differs among the three human anellovirus genera and that recombination is the main factor explaining this intergenus variability. Global characterization of anellovirus diversity could provide clues about possible associations between certain virus variants and pathologies, as well as facilitate the implementation of unbiased PCR-based detection protocols, which may be relevant for using anelloviruses as endogenous markers of immune status.

Stenotrophomonas maltophilia natural history and evolution in the airways of adults with cystic fibrosis

Introduction

Stenotrophomonas maltophilia is an opportunistic pathogen infecting persons with cystic fibrosis (pwCF) and portends a worse prognosis. Studies of S. maltophilia infection dynamics have been limited by cohort size and follow-up. We investigated the natural history, transmission potential, and evolution of S. maltophilia in a large Canadian cohort of 321 pwCF over a 37-year period.

Methods

One-hundred sixty-two isolates from 74 pwCF (23%) were typed by pulsed-field gel electrophoresis, and shared pulsotypes underwent whole-genome sequencing.

Results

S. maltophilia was recovered at least once in 82 pwCF (25.5%). Sixty-four pwCF were infected by unique pulsotypes, but shared pulsotypes were observed between 10 pwCF. In chronic carriage, longer time periods between positive sputum cultures increased the likelihood that subsequent isolates were unrelated. Isolates from individual pwCF were largely clonal, with differences in gene content being the primary source of genetic diversity objectified by gene content differences. Disproportionate progression of CF lung disease was not observed amongst those infected with multiple strains over time (versus a single) or amongst those with shared clones (versus strains only infecting one patient). We did not observe evidence of patient-to-patient transmission despite relatedness between isolates. Twenty-four genes with ≥ 2 mutations accumulated over time were identified across 42 sequenced isolates from all 11 pwCF with ≥ 2 sequenced isolates, suggesting a potential role for these genes in adaptation of S. maltophilia to the CF lung.

Discussion

Genomic analyses suggested common, indirect sources as the origins of S. maltophilia infections in the clinic population. The information derived from a genomics-based understanding of the natural history of S. maltophilia infection within CF provides unique insight into its potential for in-host evolution.

Hybridization breaks species barriers in long-term coevolution of a cyanobacterial population

Bacterial species often undergo rampant recombination yet maintain cohesive genomic identity. Ecological differences can generate recombination barriers between species and sustain genomic clusters in the short term. But can these forces prevent genomic mixing during long-term coevolution? Cyanobacteria in Yellowstone hot springs comprise several diverse species that have coevolved for hundreds of thousands of years, providing a rare natural experiment. By analyzing more than 300 single-cell genomes, we show that despite each species forming a distinct genomic cluster, much of the diversity within species is the result of hybridization driven by selection, which has mixed their ancestral genotypes. This widespread mixing is contrary to the prevailing view that ecological barriers can maintain cohesive bacterial species and highlights the importance of hybridization as a source of genomic diversity.

Parallel host shifts in a bacterial plant pathogen suggest independent genetic solutions

While there are documented host shifts in many bacterial plant pathogens, the genetic foundation of host shifts is largely unknown. Xylella fastidiosa is a bacterial pathogen found in over 600 host plant species. Two parallel host shifts occurred-in Brazil and Italy-in which X. fastidiosa adapted to infect olive trees, whereas related strains infected coffee. Using 10 novel whole-genome sequences from an olive-infecting population in Brazil, we investigated whether these olive-infecting strains diverged from closely related coffee-infecting strains. Several single-nucleotide polymorphisms, many derived from recombination events, and gene gain and loss events separated olive-infecting strains from coffee-infecting strains in this clade. The olive-specific variation suggests that this event was a host jump with genetic isolation between coffee- and olive-infecting X. fastidiosa populations. Next, we investigated the hypothesis of genetic convergence in the host shift from coffee to olive in both populations (Brazil and Italy). Each clade had multiple mutations and gene gain and loss events unique to olive, yet no overlap between clades. Using a genome-wide association study technique, we did not find any plausible candidates for convergence. Overall, this work suggests that the two populations adapted to infect olive trees through independent genetic solutions.

Phylodynamics Uncovers the Transmission of Antibiotic-Resistant Escherichia coli between Canines and Humans in an Urban Environment

The role of canines in transmitting antibiotic resistant bacteria to humans in the urban environment is poorly understood. To elucidate this role, we utilized genomic sequencing and phylogenetics to characterize the burden and transmission dynamics of antibiotic resistant Escherichia coli (ABR-Ec) cultured from canine and human feces present on urban sidewalks in San Francisco, California. We collected a total of fifty-nine ABR-Ec from human (n=12) and canine (n=47) fecal samples from the Tenderloin and South of Market (SoMa) neighborhoods of San Francisco. We then analyzed phenotypic and genotypic antibiotic resistance (ABR) of the isolates, as well as clonal relationships based on cgMLST and single nucleotide polymorphisms (SNPs) of the core genomes. Using Bayesian inference, we reconstructed the transmission dynamics between humans and canines from multiple local outbreak clusters using the marginal structured coalescent approximation (MASCOT). Overall, we found human and canine samples to carry similar amounts and profiles of ABR genes. Our results provide evidence for multiple transmission events of ABR-Ec between humans and canines. In particular, we found one instance of likely transmission from canines to humans as well as an additional local outbreak cluster consisting of one canine and one human sample. Based on this analysis, it appears that canine feces act as an important reservoir of clinically relevant ABR-Ec within the urban environment. Our findings support that public health measures should continue to emphasize proper canine feces disposal practices, access to public toilets and sidewalk and street cleaning. Importance: Antibiotic resistance in E. coli is a growing public health concern with global attributable deaths projected to reach millions annually. Current research has focused heavily on clinical routes of antibiotic resistance transmission to design interventions while the role of alternative reservoirs such as domesticated animals remain less well understood. Our results suggest canines are part of the transmission network that disseminates high-risk multidrug resistance in E. coli within the urban San Francisco community. As such, this study highlights the need to consider canines, and potentially domesticated animals more broadly, when designing interventions to reduce the prevalence of antibiotic resistance in the community. Additionally, it showcases the utility of genomic epidemiology to reconstruct the pathways by which antimicrobial resistance spreads.

Estimating the fitness cost and benefit of antimicrobial resistance from pathogen genomic data

Increasing levels of antibiotic resistance in many bacterial pathogen populations are a major threat to public health. Resistance to an antibiotic provides a fitness benefit when the bacteria are exposed to this antibiotic, but resistance also often comes at a cost to the resistant pathogen relative to susceptible counterparts. We lack a good understanding of these benefits and costs of resistance for many bacterial pathogens and antibiotics, but estimating them could lead to better use of antibiotics in a way that reduces or prevents the spread of resistance. Here, we propose a new model for the joint epidemiology of susceptible and resistant variants, which includes explicit parameters for the cost and benefit of resistance. We show how Bayesian inference can be performed under this model using phylogenetic data from susceptible and resistant lineages and that by combining data from both we are able to disentangle and estimate the resistance cost and benefit parameters separately. We applied our inferential methodology to several simulated datasets to demonstrate good scalability and accuracy. We analysed a dataset of Neisseria gonorrhoeae genomes collected between 2000 and 2013 in the USA. We found that two unrelated lineages resistant to fluoroquinolones shared similar epidemic dynamics and resistance parameters. Fluoroquinolones were abandoned for the treatment of gonorrhoea due to increasing levels of resistance, but our results suggest that they could be used to treat a minority of around 10% of cases without causing resistance to grow again.

The Multifaceted Role of Homologous Recombination in a Fastidious Bacterial Plant Pathogen

Homologous recombination plays a key function in the evolution of bacterial genomes. Within Xylella fastidiosa, an emerging plant pathogen with increasing host and geographic ranges, it has been suggested that homologous recombination facilitates host switching, speciation, and the development of virulence. We used 340 whole-genome sequences to study the relationship between inter- and intrasubspecific homologous recombination, random mutation, and natural selection across individual X. fastidiosa genes. Individual gene orthologs were identified and aligned, and a maximum likelihood (ML) gene tree was generated. Each gene alignment and tree pair were then used to calculate gene-wide and branch-specific r/m values (relative effect of recombination to mutation), gene-wide and branch-site nonsynonymous over synonymous substitution rates (dN/dS values; episodic selection), and branch length (as a proxy for mutation rate). The relationships between these variables were evaluated at the global level (i.e., for all genes among and within a subspecies), among specific functional classes (i.e., COGs), and between pangenome components (i.e., accessory versus core genes). Our analysis showed that r/m varied widely among genes as well as across X. fastidiosa subspecies. While r/m and dN/dS values were positively correlated in some instances (e.g., core genes in X. fastidiosa subsp. fastidiosa and both core and accessory genes in X. fastidiosa subsp. multiplex), low correlation coefficients suggested no clear biological significance. Overall, our results indicate that, in addition to its adaptive role in certain genes, homologous recombination acts as a homogenizing and a neutral force across phylogenetic clades, gene functional groups, and pangenome components. IMPORTANCE There is ample evidence that homologous recombination occurs frequently in the economically important plant pathogen Xylella fastidiosa. Homologous recombination has been known to occur among sympatric subspecies and is associated with host-switching events and virulence-linked genes. As a consequence, is it generally assumed that recombinant events in X. fastidiosa are adaptive. This mindset influences expectations of how homologous recombination acts as an evolutionary force as well as how management strategies for X. fastidiosa diseases are determined. Yet, homologous recombination plays roles beyond that of a source for diversification and adaptation. Homologous recombination can act as a DNA repair mechanism, as a means to facilitate nucleotide compositional change, as a homogenization mechanism within populations, or even as a neutral force. Here, we provide a first assessment of long-held beliefs regarding the general role of recombination in adaptation for X. fastidiosa. We evaluate gene-specific variations in homologous recombination rate across three X. fastidiosa subspecies and its relationship to other evolutionary forces (e.g., natural selection, mutation, etc.). These data were used to assess the role of homologous recombination in X. fastidiosa evolution.

Mutation rates and adaptive variation among the clinically dominant clusters of Mycobacterium abscessus

Mycobacterium abscessus (Mab) is a multidrug-resistant pathogen increasingly responsible for severe pulmonary infections. Analysis of whole-genome sequences (WGS) of Mab demonstrates dense genetic clustering of clinical isolates collected from disparate geographic locations. This has been interpreted as supporting patient-to-patient transmission, but epidemiological studies have contradicted this interpretation. Here, we present evidence for a slowing of the Mab molecular clock rate coincident with the emergence of phylogenetic clusters. We performed phylogenetic inference using publicly available WGS from 483 Mab patient isolates. We implement a subsampling approach in combination with coalescent analysis to estimate the molecular clock rate along the long internal branches of the tree, indicating a faster long-term molecular clock rate compared to branches within phylogenetic clusters. We used ancestry simulation to predict the effects of clock rate variation on phylogenetic clustering and found that the degree of clustering in the observed phylogeny is more easily explained by a clock rate slowdown than by transmission. We also find that phylogenetic clusters are enriched in mutations affecting DNA repair machinery and report that clustered isolates have lower spontaneous mutation rates in vitro. We propose that Mab adaptation to the host environment through variation in DNA repair genes affects the organism's mutation rate and that this manifests as phylogenetic clustering. These results challenge the model that phylogenetic clustering in Mab is explained by person-to-person transmission and inform our understanding of transmission inference in emerging, facultative pathogens.

Erythromycin resistance of clinical Campylobacter jejuni and Campylobacter coli in Shanghai, China

Campylobacter species are zoonotic pathogens, as well as the prevalent cause of foodborne bacterial gastroenteritis. The spread of antimicrobial-resistant strains poses a serious threat to global public health and attracts attention worldwide, but information about clinical Campylobacter is relatively limited compared to isolates from food and animals. The current study illustrated the prevalence and antimicrobial resistance profiles of Campylobacter jejuni and Campylobacter coli isolates collected from a consecutive surveillance program between 2012 and 2019 in Shanghai, China, using antimicrobial susceptibility testing and whole-genome sequencing. Among the 891 Campylobacter strains (761 C. jejuni and 130 C. coli) isolates collected, high portions above 90% of resistance to ciprofloxacin, nalidixic acid, and tetracycline were observed for both C. jejuni and C. coli. The most common MDR profiles represented by C. jejuni and C. coli were combination of ciprofloxacin, tetracycline, florfenicol and nalidixic acid (5.39%), and azithromycin, ciprofloxacin, erythromycin, gentamicin, tetracycline, clindamycin, nalidixic acid (28.46%), respectively. The erythromycin resistance of C. coli (59.23%) is higher than C. jejuni (2.50%). A total of 76 erythromycin resistant isolates (16 C. jejuni and 60 C. coli) were sequenced using Illumina platform for determining the genotypes, antimicrobial resistance patterns and phylogeny analysis. Multilocus sequence typing (MLST) analysis showed a high genetic diversity with 47 sequence types (STs), including 4 novel alleles and 12 new STs. The most abundant clonal complexes (CCs) were CC-403 (31.25%) and CC-828 (88.33%) for C. jejuni and C. coli, respectively. Among the 76 erythromycin-resistant isolates, mutation A2075G in 23S rRNA and erm(B) gene were detected in 53.95 and 39.47%, respectively. The erm(B) gene was identified exclusively in 30 C. coli isolates. All these erm(B) positive isolates were multi-drug resistant. Furthermore, comparison of the erm(B)-carrying isolates of multiple sources worldwide demonstrated the possibility of zoonotic transmission of erm(B) in Campylobacter. These findings highlight the importance of continuous surveillance of erythromycin resistance dissemination in Campylobacter which may compromise the effectiveness of antimicrobial therapy.

Phylogenetic Analysis of Escherichia coli Isolated from Australian Feedlot Cattle in Comparison to Pig Faecal and Poultry/Human Extraintestinal Isolates

The similarity of commensal Escherichia coli isolated from healthy cattle to antimicrobial-resistant bacteria causing extraintestinal infections in humans is not fully understood. In this study, we used a bioinformatics approach based on whole genome sequencing data to determine the genetic characteristics and phylogenetic relationships among faecal Escherichia coli isolates from beef cattle (n = 37) from a single feedlot in comparison to previously analysed pig faecal (n = 45), poultry extraintestinal (n = 19), and human extraintestinal E. coli isolates (n = 40) from three previous Australian studies. Most beef cattle and pig isolates belonged to E. coli phylogroups A and B1, whereas most avian and human isolates belonged to B2 and D, although a single human extraintestinal isolate belonged to phylogenetic group A and sequence type (ST) 10. The most common E. coli sequence types (STs) included ST10 for beef cattle, ST361 for pig, ST117 for poultry, and ST73 for human isolates. Extended-spectrum and AmpC β-lactamase genes were identified in seven out of thirty-seven (18.9%) beef cattle isolates. The most common plasmid replicons identified were IncFIB (AP001918), followed by IncFII, Col156, and IncX1. The results confirm that feedlot cattle isolates examined in this study represent a reduced risk to human and environmental health with regard to being a source of antimicrobial-resistant E. coli of clinical importance.

Host Association and Spatial Proximity Shape but Do Not Constrain Population Structure in the Mutualistic Symbiont Xenorhabdus bovienii

To what extent are generalist species cohesive evolutionary units rather than a compilation of recently diverged lineages? We examine this question in the context of host specificity and geographic structure in the insect pathogen and nematode mutualist Xenorhabdus bovienii. This bacterial species partners with multiple nematode species across two clades in the genus Steinernema. We sequenced the genomes of 42 X. bovienii strains isolated from four different nematode species and three field sites within a 240-km² region and compared them to globally available reference genomes. We hypothesized that X. bovienii would comprise several host-specific lineages, such that bacterial and nematode phylogenies would be largely congruent. Alternatively, we hypothesized that spatial proximity might be a dominant signal, as increasing geographic distance might lower shared selective pressures and opportunities for gene flow. We found partial support for both hypotheses. Isolates clustered largely by nematode host species but did not strictly match the nematode phylogeny, indicating that shifts in symbiont associations across nematode species and clades have occurred. Furthermore, both genetic similarity and gene flow decreased with geographic distance across nematode species, suggesting differentiation and constraints on gene flow across both factors, although no absolute barriers to gene flow were observed across the regional isolates. Several genes associated with biotic interactions were found to be undergoing selective sweeps within this regional population. The interactions included several insect toxins and genes implicated in microbial competition. Thus, gene flow maintains cohesiveness across host associations in this symbiont and may facilitate adaptive responses to a multipartite selective environment. IMPORTANCE Microbial populations and species are notoriously hard to delineate. We used a population genomics approach to examine the population structure and the spatial scale of gene flow in Xenorhabdus bovienii, an intriguing species that is both a specialized mutualistic symbiont of nematodes and a broadly virulent insect pathogen. We found a strong signature of nematode host association, as well as evidence for gene flow connecting isolates associated with different nematode host species and collected from distinct study sites. Furthermore, we saw signatures of selective sweeps for genes involved with nematode host associations, insect pathogenicity, and microbial competition. Thus, X. bovienii exemplifies the growing consensus that recombination not only maintains cohesion but can also allow the spread of niche-beneficial alleles.

Host phylogeny and ecological associations best explain Wolbachia host shifts in scale insects

Wolbachia are among the most prevalent and widespread endosymbiotic bacteria on Earth. Wolbachia's success in infecting an enormous number of arthropod species is attributed to two features: the range of phenotypes they induce in their hosts, and their ability to switch between host species. Whilst much progress has been made in elucidating their induced phenotypes, our understanding of Wolbachia host-shifting is still very limited: we lack answers to even fundamental questions concerning Wolbachia's routes of transfer and the importance of factors influencing host shifts. Here, we investigate the diversity and host-shift patterns of Wolbachia in scale insects, a group of arthropods with intimate associations with other insects that make them well suited to studying host shifts. Using Illumina multitarget amplicon sequencing of Wolbachia-infected scale insects and their direct associates we determined the identity of all Wolbachia strains. We then fitted a generalized additive mixed model to our data to estimate the influence of host phylogeny and the geographical distribution on Wolbachia strain sharing among scale insect species. The model predicts no significant contribution of host geography but strong effects of host phylogeny, with high rates of Wolbachia sharing among closely related species and a sudden drop-off in sharing with increasing phylogenetic distance. We also detected the same Wolbachia strain in scale insects and several intimately associated species (ants, wasps and flies). This indicates putative host shifts and potential routes of transfers via these associates and highlights the importance of ecological connectivity in Wolbachia host-shifting.

Genome-wide association reveals host-specific genomic traits in Escherichia coli

BACKGROUND

Escherichia coli is an opportunistic pathogen which colonizes various host species. However, to what extent genetic lineages of E. coli are adapted or restricted to specific hosts and the genomic determinants of such adaptation or restriction is poorly understood.

RESULTS

We randomly sampled E. coli isolates from four countries (Germany, UK, Spain, and Vietnam), obtained from five host species (human, pig, cattle, chicken, and wild boar) over 16 years, from both healthy and diseased hosts, to construct a collection of 1198 whole-genome sequenced E. coli isolates. We identified associations between specific E. coli lineages and the host from which they were isolated. A genome-wide association study (GWAS) identified several E. coli genes that were associated with human, cattle, or chicken hosts, whereas no genes associated with the pig host could be found. In silico characterization of nine contiguous genes (collectively designated as nan-9) associated with the human host indicated that these genes are involved in the metabolism of sialic acids (Sia). In contrast, the previously described sialic acid regulon known as sialoregulon (i.e. nanRATEK-yhcH, nanXY, and nanCMS) was not associated with any host species. In vitro growth experiments with a Δnan-9 E. coli mutant strain, using the sialic acids 5-N-acetylneuraminic acid (Neu5Ac) and N-glycolylneuraminic acid (Neu5Gc) as sole carbon source, showed impaired growth behaviour compared to the wild-type.

CONCLUSIONS

This study provides an extensive analysis of genetic determinants which may contribute to host specificity in E. coli. Our findings should inform risk analysis and epidemiological monitoring of (antimicrobial resistant) E. coli.

Genomic surveillance uncovers a pandemic clonal lineage of the wheat blast fungus

Wheat, one of the most important food crops, is threatened by a blast disease pandemic. Here, we show that a clonal lineage of the wheat blast fungus recently spread to Asia and Africa following two independent introductions from South America. Through a combination of genome analyses and laboratory experiments, we show that the decade-old blast pandemic lineage can be controlled by the Rmg8 disease resistance gene and is sensitive to strobilurin fungicides. However, we also highlight the potential of the pandemic clone to evolve fungicide-insensitive variants and sexually recombine with African lineages. This underscores the urgent need for genomic surveillance to track and mitigate the spread of wheat blast outside of South America and to guide preemptive wheat breeding for blast resistance.

Genomic surveillance uncovers a pandemic clonal lineage of the wheat blast fungus

Wheat, one of the most important food crops, is threatened by a blast disease pandemic. Here, we show that a clonal lineage of the wheat blast fungus recently spread to Asia and Africa following two independent introductions from South America. Through a combination of genome analyses and laboratory experiments, we show that the decade-old blast pandemic lineage can be controlled by the Rmg8 disease resistance gene and is sensitive to strobilurin fungicides. However, we also highlight the potential of the pandemic clone to evolve fungicide-insensitive variants and sexually recombine with African lineages. This underscores the urgent need for genomic surveillance to track and mitigate the spread of wheat blast outside of South America and to guide preemptive wheat breeding for blast resistance.

An integrated eco-evolutionary framework to predict population-level responses of climate-sensitive pathogens

It is critical to gain insight into how climate change impacts evolutionary responses within climate-sensitive pathogen populations, such as increased resilience, opportunistic responses and the emergence of dominant variants from highly variable genomic backgrounds and subsequent global dispersal. This review proposes a framework to support such analysis, by combining genomic evolutionary analysis with climate time-series data in a novel spatiotemporal dataframe for use within machine learning applications, to understand past and future evolutionary pathogen responses to climate change. Recommendations are presented to increase the feasibility of interdisciplinary applications, including the importance of robust spatiotemporal metadata accompanying genome submission to databases. Such workflows will inform accessible public health tools and early-warning systems, to aid decision-making and mitigate future human health threats.

Unlocking the microbial studies through computational approaches: how far have we reached?

The metagenomics approach accelerated the study of genetic information from uncultured microbes and complex microbial communities. In silico research also facilitated an understanding of protein-DNA interactions, protein-protein interactions, docking between proteins and phyto/biochemicals for drug design, and modeling of the 3D structure of proteins. These in silico approaches provided insight into analyzing pathogenic and nonpathogenic strains that helped in the identification of probable genes for vaccines and antimicrobial agents and comparing whole-genome sequences to microbial evolution. Artificial intelligence, more precisely machine learning (ML) and deep learning (DL), has proven to be a promising approach in the field of microbiology to handle, analyze, and utilize large data that are generated through nucleic acid sequencing and proteomics. This enabled the understanding of the functional and taxonomic diversity of microorganisms. ML and DL have been used in the prediction and forecasting of diseases and applied to trace environmental contaminants and environmental quality. This review presents an in-depth analysis of the recent application of silico approaches in microbial genomics, proteomics, functional diversity, vaccine development, and drug design.

Rhometa: Population recombination rate estimation from metagenomic read datasets

Prokaryotic evolution is influenced by the exchange of genetic information between species through a process referred to as recombination. The rate of recombination is a useful measure for the adaptive capacity of a prokaryotic population. We introduce Rhometa (https://github.com/sid-krish/Rhometa), a new software package to determine recombination rates from shotgun sequencing reads of metagenomes. It extends the composite likelihood approach for population recombination rate estimation and enables the analysis of modern short-read datasets. We evaluated Rhometa over a broad range of sequencing depths and complexities, using simulated and real experimental short-read data aligned to external reference genomes. Rhometa offers a comprehensive solution for determining population recombination rates from contemporary metagenomic read datasets. Rhometa extends the capabilities of conventional sequence-based composite likelihood population recombination rate estimators to include modern aligned metagenomic read datasets with diverse sequencing depths, thereby enabling the effective application of these techniques and their high accuracy rates to the field of metagenomics. Using simulated datasets, we show that our method performs well, with its accuracy improving with increasing numbers of genomes. Rhometa was validated on a real S. pneumoniae transformation experiment, where we show that it obtains plausible estimates of the rate of recombination. Finally, the program was also run on ocean surface water metagenomic datasets, through which we demonstrate that the program works on uncultured metagenomic datasets.

Phylogroup-specific variation shapes the clustering of antimicrobial resistance genes and defence systems across regions of genome plasticity in Pseudomonas aeruginosa

BACKGROUND

Pseudomonas aeruginosa is an opportunistic pathogen consisting of three phylogroups (hereafter named A, B, and C). Here, we assessed phylogroup-specific evolutionary dynamics across available and also new P. aeruginosa genomes.

METHODS

In this genomic analysis, we first generated new genome assemblies for 18 strains of the major P. aeruginosa clone type (mPact) panel, comprising a phylogenetically diverse collection of clinical and environmental isolates for this species. Thereafter, we combined these new genomes with 1991 publicly available P. aeruginosa genomes for a phylogenomic and comparative analysis. We specifically explored to what extent antimicrobial resistance (AMR) genes, defence systems, and virulence genes vary in their distribution across regions of genome plasticity (RGPs) and "masked" (RGP-free) genomes, and to what extent this variation differs among the phylogroups.

FINDINGS

We found that members of phylogroup B possess larger genomes, contribute a comparatively larger number of pangenome families, and show lower abundance of CRISPR-Cas systems. Furthermore, AMR and defence systems are pervasive in RGPs and integrative and conjugative/mobilizable elements (ICEs/IMEs) from phylogroups A and B, and the abundance of these cargo genes is often significantly correlated. Moreover, inter- and intra-phylogroup interactions occur at the accessory genome level, suggesting frequent recombination events. Finally, we provide here the mPact panel of diverse P. aeruginosa strains that may serve as a valuable reference for functional analyses.

INTERPRETATION

Altogether, our results highlight distinct pangenome characteristics of the P. aeruginosa phylogroups, which are possibly influenced by variation in the abundance of CRISPR-Cas systems and are shaped by the differential distribution of other defence systems and AMR genes.

FUNDING

German Science Foundation, Max-Planck Society, Leibniz ScienceCampus Evolutionary Medicine of the Lung, BMBF program Medical Infection Genomics, Kiel Life Science Postdoc Award.

Repeated Emergence of Variant TetR Family Regulator, FarR, and Increased Resistance to Antimicrobial Unsaturated Fatty Acid among Clonal Complex 5 Methicillin-Resistant Staphylococcus aureus

Resistance-nodulation-division (RND) superfamily efflux pumps promote antibiotic resistance in Gram-negative pathogens, but their role in Gram-positive pathogens, including methicillin-resistant Staphylococcus aureus (MRSA) is undocumented. However, recent in vitro selections for resistance of S. aureus to an antimicrobial fatty acid, linoleic acid, and an antibiotic, rhodomyrtone, identified H121Y and C116R substitution variants, respectively, in a TetR family regulator, FarR, promoting increased expression of the RND pump FarE. Hypothesizing that in vivo selection pressures have also promoted the emergence of FarR variants, we searched available genome data and found that strains with FarR^H121Y from human and bovine hosts have emerged sporadically in clonal complexes (CCs) CC1, CC30, CC8, CC22, and CC97, whereas multiple FarR variants have occurred within CC5 hospital-associated (HA)-MRSA. Of these, FarR^E160G and FarR^E93EE were exclusive to CC5, while FarR^C116Y, FarR^P165L, and FarR^G166D also occurred in nonrelated CCs, primarily from bovine hosts. Within CC5, FarR^C116Y and FarR^G166D strains were polyphyletic, each exhibiting two emergence events. FarR^C116Y and FarR^E160G were individually sufficient to confer increased expression of FarE and enhanced resistance to linoleic acid (LA). Isolates with FarR^E93EE were most closely related to S. aureus N315 MRSA and exhibited increased resistance independently of FarR^E93EE. Accumulation of pseudogenes and additional polymorphisms in FarR^E93EE strains contributed to a multiresistance phenotype which included fosfomycin and fusidic acid resistance in addition to increased linoleic acid resistance. These findings underscore the remarkable adaptive capacity of CC5 MRSA, which includes the polyphyletic USA100 lineage of HA-MRSA that is endemic in the Western hemisphere and known for the acquisition of multiple resistance phenotypes.

Pan-Chromosome and Comparative Analysis of Agrobacterium fabrum Reveal Important Traits Concerning the Genetic Diversity, Evolutionary Dynamics, and Niche Adaptation of the Species

Agrobacterium fabrum has been critical for the development of plant genetic engineering and agricultural biotechnology due to its ability to transform eukaryotic cells. However, the gene composition, evolutionary dynamics, and niche adaptation of this species is still unknown. Therefore, we established a comparative genomic analysis based on a pan-chromosome data set to evaluate the genetic diversity of A. fabrum. Here, 25 A. fabrum genomes were selected for analysis by core genome phylogeny combined with the average nucleotide identity (ANI), amino acid identity (AAI), and in silico DNA-DNA hybridization (DDH) values. An open pan-genome of A. fabrum exhibits genetic diversity with variable accessorial genes as evidenced by a consensus pan-genome of 12 representative genomes. The genomic plasticity of A. fabrum is apparent in its putative sequences for mobile genetic elements (MGEs), limited horizontal gene transfer barriers, and potentially horizontally transferred genes. The evolutionary constraints and functional enrichment in the pan-chromosome were measured by the Clusters of Orthologous Groups (COG) categories using eggNOG-mapper software, and the nonsynonymous/synonymous rate ratio (dN/dS) was determined using HYPHY software. Comparative analysis revealed significant differences in the functional enrichment and the degree of purifying selection between the core genome and non-core genome. We demonstrate that the core gene families undergo stronger purifying selection but have a significant bias to contain one or more positively selected sites. Furthermore, although they shared similar genetic diversity, we observed significant differences between chromosome 1 (Chr I) and the chromid in their functional features and evolutionary constraints. We demonstrate that putative genetic elements responsible for plant infection, ecological adaptation, and speciation represent the core genome, highlighting their importance in the adaptation of A. fabrum to plant-related niches. Our pan-chromosome analysis of A. fabrum provides comprehensive insights into the genetic properties, evolutionary patterns, and niche adaptation of the species. IMPORTANCE Agrobacterium spp. live in diverse plant-associated niches such as soil, the rhizosphere, and vegetation, which are challenged by multiple stressors such as diverse energy sources, plant defenses, and microbial competition. They have evolved the ability to utilize diverse resources, escape plant defenses, and defeat competitors. However, the underlying genetic diversity and evolutionary dynamics of Agrobacterium spp. remain unexplored. We examined the phylogeny and pan-genome of A. fabrum to define intraspecies evolutionary relationships. Our results indicate an open pan-genome and numerous MGEs and horizontally transferred genes among A. fabrum genomes, reflecting the flexibility of the chromosomes and the potential for genetic exchange. Furthermore, we observed significant differences in the functional features and evolutionary constraints between the core and accessory genomes and between Chr I and the chromid, respectively.

Helicobacter pylori Chronic Infection Selects for Effective Colonizers of Metaplastic Glands

Chronic gastric infection with Helicobacter pylori can lead to progressive tissue changes that culminate in cancer, but how H. pylori adapts to the changing tissue environment during disease development is not fully understood. In a transgenic mouse gastric metaplasia model, we found that strains from unrelated individuals differed in their ability to infect the stomach, to colonize metaplastic glands, and to alter the expression of the metaplasia-associated protein TFF3. H. pylori isolates from different stages of disease from a single individual had differential ability to colonize healthy and metaplastic gastric glands. Exposure to the metaplastic environment selected for high gastric colonization by one of these strains. Complete genome sequencing revealed a unique alteration in the frequency of a variant allele of the putative adhesin sabB, arising from a recombination event with the related sialic acid binding adhesin (SabA) gene. Mutation of sabB in multiple H. pylori strain backgrounds strongly reduced adherence to both normal and metaplastic gastric tissue, and highly attenuated stomach colonization in mice. Thus, the changing gastric environment during disease development promotes bacterial adhesin gene variation associated with enhanced gastric colonization. IMPORTANCE Chronic infection with Helicobacter pylori is the primary risk factor for developing stomach cancer. As disease progresses H. pylori must adapt to a changing host tissue environment that includes induction of new cell fates in the cells that line the stomach. We tested representative H. pylori isolates collected from the same patient during early and later stages of disease in a mouse model where we can rapidly induce disease-associated tissue changes. Only the later-stage H. pylori strains could robustly colonize the diseased stomach environment. We also found that the ability to colonize the diseased stomach was associated with genetic variation in a putative cell surface adhesin gene called sabB. Additional experiments revealed that SabB promotes binding to stomach tissue and is critical for stomach colonization by the late-stage strains. Thus, H. pylori diversifies its genome during disease progression and these genomic changes highlight critical factors for bacterial persistence.

Genomic diversity in Fructobacillus spp. isolated from fructose-rich niches

The Fructobacillus genus is a group of obligately fructophilic lactic acid bacteria (FLAB) that requires the use of fructose or another electron acceptor for their growth. In this work, we performed a comparative genomic analysis within the genus Fructobacillus by using 24 available genomes to evaluate genomic and metabolic differences among these organisms. In the genome of these strains, which varies between 1.15- and 1.75-Mbp, nineteen intact prophage regions, and seven complete CRISPR-Cas type II systems were found. Phylogenetic analyses located the studied genomes in two different clades. A pangenome analysis and a functional classification of their genes revealed that genomes of the first clade presented fewer genes involved in the synthesis of amino acids and other nitrogen compounds. Moreover, the presence of genes strictly related to the use of fructose and electron acceptors was variable within the genus, although these variations were not always related to the phylogeny.

Comparative Genomic and Metagenomic Investigations of the Corynebacterium tuberculostearicum Species Complex Reveals Potential Mechanisms Underlying Associations To Skin Health and Disease

Corynebacterium are a diverse genus and dominant member of the human skin microbiome. Recently, we reported that the most prevalent Corynebacterium species found on skin, including Corynebacterium tuberculostearicum and Corynebacterium kefirresidentii, comprise a narrow species complex despite the diversity of the genus. Here, we apply high-resolution phylogenomics and comparative genomics to describe the structure of the C. tuberculostearicum species complex and highlight genetic traits which are enriched or depleted in it relative to other Corynebacterium. Through metagenomic investigations, we also find that individual species within the complex can associate with specific body sites. Finally, we discover that one species from the complex, C. kefirresidentii, increases in relative abundance during atopic dermatitis flares, and show that most genomes of this species encode a colocalized set of putative virulence genes. IMPORTANCE Corynebacterium are commonly found bacteria on the human skin. In this study, we perform comparative genomics to gain insight into genetic traits which differentiate a phylogenetically related group of Corynebacterium, the Corynebacterium tuberculostearicum species complex, that includes the most prevalent species from the genus in skin microbiomes. After resolving the presence of distinct species within the complex, we applied metagenomic analysis to uncover biogeographic associations of individual species within the complex with specific body sites and discovered that one species, commonly found in the nares of individuals, increases in abundance across multiple body sites during atopic dermatitis flares.

Antimicrobial Resistance of and Genomic Insights into Pasteurella multocida Strains Isolated from Australian Pigs

Infection with Pasteurella multocida represents a significant economic threat to Australian pig producers, yet our knowledge of its antimicrobial susceptibilities is lagging, and genomic characterization of P. multocida strains associated with porcine lower respiratory disease is internationally scarce. This study utilized high-throughput robotics to phenotypically and genetically characterize an industry-wide collection of 252 clinical P. multocida isolates that were recovered between 2014 and 2019. Overall, antimicrobial resistance was found to be low, with clinical resistance below 1% for all tested antimicrobials except those from the tetracycline class. Five dominant sequence types, representing 64.8% of all isolates, were identified; they were disseminated across farms and had previously been detected in various animal hosts and countries. P. multocida in Australian farms remain controllable via current antimicrobial therapeutic protocols. The identification of highly dominant, interspecies-infecting strains provides insight into the epidemiology of the opportunistic pathogen, and it highlights a biosecurity threat to the Australian livestock industry. IMPORTANCE Pasteurellosis is rated by the World Animal Health Organisation (OIE) as a high-impact disease in livestock. Although it is well understood in many host-disease contexts, our understanding of the organism in porcine respiratory disease is limited. Given its high frequency of involvement in porcine respiratory disease complex (PRDC), it is important that we are aware of its antimicrobial susceptibilities so that we can respond quickly and appropriately with antimicrobial therapy. Genetic insights about the organism can help us to better understand its epidemiology and inform our biosecurity practices and prophylactic management.

Transcriptome and Proteome of Methicillin-Resistant Staphylococcus aureus Small-Colony Variants Reveal Changed Metabolism and Increased Immune Evasion

Methicillin-resistant Staphylococcus aureus (MRSA) infection has become a public health crisis. Recently, we isolated small-colony variants (SCVs) of MRSA, which are characterized by slow growth, decreased virulence, increased antibiotic resistance, and immune evasion. In the present study, we provided proteomic and transcriptomic profiles of clinical MRSA sequence type 239 (ST239) normal strains and SCVs and attempted to identify the key genes or pathways closely related to SCV formation and survival. RNAs and proteins were extracted and subjected to RNA sequencing and mass spectrometry, and the transcriptome and proteome were evaluated via bioinformatic analysis. The results were verified by functional assays. In total, 822 differentially expressed genes (DEGs) and 773 differentially expressed proteins (DEPs) were identified; of these, 286 DEGs and DEPs were correlated and subjected to Kyoto Encyclopedia Genes and Genomes analysis. Some pathways were significant, including ABC transporters, ribosome biogenesis, and metabolic pathways such as glycolysis/gluconeogenesis and the citrate cycle (tricarboxylic acid [TCA] cycle). Based on these results, we found that the downregulation of ABC transporters and the TCA cycle pathway resulted in electron transport chain deficiencies and reduced ATP production in SCVs, leading to a dependence on glycolysis and its upregulation. In addition, the upregulation of capsule polysaccharides and the downregulation of surface proteins prevented phagocytosis and reduced the adhesion of host cells, contributing to immune evasion by SCVs. These findings contribute to a better understanding of the mechanisms of SCV formation and survival. IMPORTANCE Small-colony variants (SCVs) of Staphylococcus aureus have drawn increasing research attention. Owing to their slow growth, atypical colony morphology, and unusual metabolic characteristics, SCVs often cause confusion in the laboratory. Furthermore, clinical treatment of SCVs is challenging owing to their antibiotic resistance and immune evasion, leading to persistent and recurrent infections. However, the mechanisms underlying their formation remain unclear. In this study, we isolated SCVs of methicillin-resistant S. aureus and provided transcriptomic and proteomic profiles of normal strains and SCVs. Based on our analysis, glycolysis upregulation and TCA cycle downregulation affected the electron transport chain and energy supply, leading to slower metabolism. Moreover, capsular biosynthesis was increased, while the number of surface proteins decreased, thus promoting immune evasion by SCVs.

Dynamics of the MRSA Population in A Chilean Hospital: A Phylogenomic Analysis (2000-2016)

The global dissemination of methicillin-resistant Staphylococcus aureus (MRSA) is associated with the emergence and establishment of clones in specific geographic areas. The Chilean-Cordobes clone (ChC) (ST5-SCC mec I) has been the predominant MRSA clone in Chile since its first description in 1998, despite the report of other emerging MRSA clones in the last years. Here, we characterize the evolutionary history of MRSA from 2000 to 2016 in a Chilean tertiary healthcare center using phylogenomic analyses. We sequenced 469 MRSA isolates collected between 2000-2016 in a tertiary healthcare center in Chile. We evaluated the temporal trends of the circulating clones and performed a phylogenomic reconstruction to characterize the clonal dynamics. We found a significant increase in the diversity and richness of sequence types (STs; Spearman r=0.8748, p<0.0001) with a Shannon diversity index increasing from 0.221 in the year 2000 to 1.33 in 2016. The temporal trend analysis revealed that in the period 2000-2003 most of the isolates (94.2%; n=98) belonged to the ChC clone. However, since then, the frequency of the ChC clone has decreased over time, accounting for 52% of the collection in the 2013-2016 period. This decline was accompanied by the rise of two emerging MRSA lineages, ST105-SCC mec II and ST72-SCC mec VI. In conclusion, the ChC clone remains the most frequent MRSA lineage in Chile. However, this lineage is gradually being replaced by several emerging clones, the most important of which is clone ST105-SCC mec II. To the best of our knowledge, this is the largest study of MRSA clonal dynamics performed in South America.

Importance

Methicillin-resistant Staphylococcus aureus (MRSA) is a major public health pathogen that disseminates through the emergence of successful dominant clones in specific geographic regions. Knowledge of the dissemination and molecular epidemiology of MRSA in Latin America is scarce and is largely based on small studies or classical typing techniques with several limitations to depict an accurate description of their genomic landscape. We used whole-genome sequencing to study 469 MRSA isolates collected between 2000-2016 in Chile to provide the largest and most detailed study of clonal dynamics of MRSA carried out in South America to date. We found a significant increase in the diversity of MRSA clones circulating over the 17-year study period. Additionally, we describe the emergence of two novel clones (ST105-SCCmecII and ST72-SCCmecVI), which have been gradually increasing their frequency over time. Our results drastically improve our understanding of the dissemination and update our knowledge about MRSA in Latin America.

Identification of evolutionary trajectories shared across human betacoronaviruses

Comparing the evolution of distantly related viruses can provide insights into common adaptive processes related to shared ecological niches. Phylogenetic approaches, coupled with other molecular evolution tools, can help identify mutations informative on adaptation, whilst the structural contextualization of these to functional sites of proteins may help gain insight into their biological properties. Two zoonotic betacoronaviruses capable of sustained human-to-human transmission have caused pandemics in recent times (SARS-CoV-1 and SARS-CoV-2), whilst a third virus (MERS-CoV) is responsible for sporadic outbreaks linked to animal infections. Moreover, two other betacoronaviruses have circulated endemically in humans for decades (HKU1 and OC43). To search for evidence of adaptive convergence between established and emerging betacoronaviruses capable of sustained human-to-human transmission (HKU1, OC43, SARS-CoV-1 and SARS-CoV-2), we developed a methodological pipeline to classify shared non-synonymous mutations as putatively denoting homoplasy (repeated mutations that do not share direct common ancestry) or stepwise evolution (sequential mutations leading towards a novel genotype). In parallel, we look for evidence of positive selection, and draw upon protein structure data to identify potential biological implications. We find 30 mutations, with four of these [codon sites 18121 (nsp14/residue 28), 21623 (spike/21), 21635 (spike/25) and 23948 (spike/796); SARS-CoV-2 genome numbering] displaying evolution under positive selection and proximity to functional protein regions. Our findings shed light on potential mechanisms underlying betacoronavirus adaptation to the human host and pinpoint common mutational pathways that may occur during establishment of human endemicity.

Postglacial adaptations enabled colonization and quasi-clonal dispersal of ammonia-oxidizing archaea in modern European large lakes

Ammonia-oxidizing archaea (AOA) play a key role in the aquatic nitrogen cycle. Their genetic diversity is viewed as the outcome of evolutionary processes that shaped ancestral transition from terrestrial to marine habitats. However, current genome-wide insights into AOA evolution rarely consider brackish and freshwater representatives or provide their divergence timeline in lacustrine systems. An unbiased global assessment of lacustrine AOA diversity is critical for understanding their origins, dispersal mechanisms, and ecosystem roles. Here, we leveraged continental-scale metagenomics to document that AOA species diversity in freshwater systems is remarkably low compared to marine environments. We show that the uncultured freshwater AOA, "Candidatus Nitrosopumilus limneticus," is ubiquitous and genotypically static in various large European lakes where it evolved 13 million years ago. We find that extensive proteome remodeling was a key innovation for freshwater colonization of AOA. These findings reveal the genetic diversity and adaptive mechanisms of a keystone species that has survived clonally in lakes for millennia.

Decoding geobiological evolution from microbiomes

Genomic records of genetic recombination and mutation rates indicate that freshwater ammonia-oxidizing archaea have evolved through paleoclimate and geohydrological history.

Genome-Wide Association with Uncertainty in the Genetic Similarity Matrix

Genome-wide association studies (GWASs) are often confounded by population stratification and structure. Linear mixed models (LMMs) are a powerful class of methods for uncovering genetic effects, while controlling for such confounding. LMMs include random effects for a genetic similarity matrix, and they assume that a true genetic similarity matrix is known. However, uncertainty about the phylogenetic structure of a study population may degrade the quality of LMM results. This may happen in bacterial studies in which the number of samples or loci is small, or in studies with low-quality genotyping. In this study, we develop methods for linear mixed models in which the genetic similarity matrix is unknown and is derived from Markov chain Monte Carlo estimates of the phylogeny. We apply our model to a GWAS of multidrug resistance in tuberculosis, and illustrate our methods on simulated data.

Genomic characterization and molecular dating of the novel bacterium Permianibacter aggregans HW001T, which originated from Permian ground water

The Permian Basin is a unique ecosystem located in the southwest of the USA. An unanswered question is whether the bacteria in the Permian Basin adapted to the changing paleomarine environment and survived in the remnants of Permian groundwater. In our previous study, a novel bacterial strain, Permianibacter aggregans HW001T, was isolated from microalgae cultures incubated with Permian Basin waters, and was shown to originate from the Permian Ocean. In this study, strain HW001T was shown to be the representative strain of a novel family, classified as 'Permianibacteraceae'. The results of molecular dating suggested that the strain HW001T diverged ~ 447 million years ago (mya), which is the early Permian period (~ 250 mya). Genome analysis was used to access its potential energy utilization and biosynthesis capacity. A large number of transporters, carbohydrate-active enzymes and protein-degradation related genes have been annotated in the genome of strain HW001T. In addition, a series of important metabolic pathways, such as peptidoglycan biosynthesis, osmotic stress response system and multifunctional quorum sensing were annotated, which may confer the ability to adapt to various unfavorable environmental conditions. Finally, the evolutionary history of strain HW001T was reconstructed and the horizontal transfer of genes was predicted, indicating that the adaptation of P. aggregans to a changing marine environment depends on the evolution of their metabolic capabilities, especially in signal transmission. In conclusion, the results of this study provide genomic information for revealing the adaptive mechanism of strain HW001T to the changing ancient oceans.

Supplementary Information

The online version contains supplementary material available at 10.1007/s42995-023-00164-3.

Streptococcus suis outbreak caused by an emerging zoonotic strain with acquired multi-drug resistance in Thailand

Streptococcus suis is an emerging zoonotic swine pathogen which can cause severe infections in humans. In March 2021, an outbreak of S. suis infections with 19 confirmed cases of septicemia and meningitis leading to two deaths, occurred in Nakhon Ratchasima province, Thailand. We characterized the outbreak through an epidemiological investigation combined with Illumina and Nanopore whole genome sequencing (WGS). The source of the outbreak was traced back to a raw pork dish prepared from a single pig during a Buddhist ceremony attended by 241 people. WGS analysis revealed that a single S. suis serotype 2 strain belonging to a novel sequence type (ST) of the emergent Thai zoonotic clade CC233/379, was responsible for the infections. The outbreak clone grouped together with other Thai zoonotic strains from CC233/379 and CC104 in a global S. suis phylogeny and capsule switching events between serotype 2 zoonotic strains and serotype 7 porcine strains were identified. The outbreak strain showed reduced susceptibility to penicillin corresponding with mutations in key residues in the penicillin binding proteins (PBPs). Furthermore, the outbreak strain was resistant to tetracycline, erythromycin, clindamycin, linezolid and chloramphenicol, having acquired an integrative and conjugative element (ICE) carrying resistance genes tetO and ermB, as well as a transposon from the IS1216 family carrying optrA and ermA. This investigation demonstrates that multi-drug resistant zoonotic lineages of S. suis which pose a threat to human health continue to emerge.

Unstable EBV latency drives inflammation in multiple sclerosis patient derived spontaneous B cells

Epidemiological studies have demonstrated that Epstein-Barr virus (EBV) is a known etiologic risk factor, and perhaps prerequisite, for the development of MS. EBV establishes life-long latent infection in a subpopulation of memory B cells. Although the role of memory B cells in the pathobiology of MS is well established, studies characterizing EBV-associated mechanisms of B cell inflammation and disease pathogenesis in EBV (+) B cells from MS patients are limited. Accordingly, we analyzed spontaneous lymphoblastoid cell lines (SLCLs) from multiple sclerosis patients and healthy controls to study host-virus interactions in B cells, in the context of an individual's endogenous EBV. We identify differences in EBV gene expression and regulation of both viral and cellular genes in SLCLs. Our data suggest that EBV latency is dysregulated in MS SLCLs with increased lytic gene expression observed in MS patient B cells, especially those generated from samples obtained during "active" disease. Moreover, we show increased inflammatory gene expression and cytokine production in MS patient SLCLs and demonstrate that tenofovir alafenamide, an antiviral that targets EBV replication, decreases EBV viral loads, EBV lytic gene expression, and EBV-mediated inflammation in both SLCLs and in a mixed lymphocyte assay. Collectively, these data suggest that dysregulation of EBV latency in MS drives a pro-inflammatory, pathogenic phenotype in memory B cells and that this response can be attenuated by suppressing EBV lytic activation. This study provides further support for the development of antiviral agents that target EBV-infection for use in MS.

Pre-epidemic evolution of the MRSA USA300 clade and a molecular key for classification

Introduction

USA300 has remained the dominant community and healthcare associated methicillin-resistant Staphylococcus aureus (MRSA) clone in the United States and in northern South America for at least the past 20 years. In this time, it has experienced epidemic spread in both of these locations. However, its pre-epidemic evolutionary history and origins are incompletely understood. Large sequencing databases, such as NCBI, PATRIC, and Staphopia, contain clues to the early evolution of USA300 in the form of sequenced genomes of USA300 isolates that are representative of lineages that diverged prior to the establishment of the South American epidemic (SAE) clade and North American epidemic (NAE) clade. In addition, historical isolates collected prior to the emergence of epidemics can help reconstruct early events in the history of this lineage.

Methods

Here, we take advantage of the accrued, publicly available data, as well as two newly sequenced pre-epidemic historical isolates from 1996, and a very early diverging ACME-negative NAE genome, to understand the pre-epidemic evolution of USA300. We use database mining techniques to emphasize genomes similar to pre-epidemic isolates, with the goal of reconstructing the early molecular evolution of the USA300 lineage.

Results

Phylogenetic analysis with these genomes confirms that the NAE and SAE USA300 lineages diverged from a most recent common ancestor around 1970 with high confidence, and it also pinpoints the independent acquisition events of the of the ACME and COMER loci with greater precision than in previous studies. We provide evidence for a North American origin of the USA300 lineage and identify multiple introductions of USA300 into South and North America. Notably, we describe a third major USA300 clade (the pre-epidemic branching clade; PEB1) consisting of both MSSA and MRSA isolates circulating around the world that diverged from the USA300 lineage prior to the establishment of the South and North American epidemics. We present a detailed analysis of specific sequence characteristics of each of the major clades, and present diagnostic positions that can be used to classify new genomes.

Population Structure and Genomic Characteristics of Australian Erysipelothrix rhusiopathiae Reveals Unobserved Diversity in the Australian Pig Industry

Erysipelothrix rhusiopathiae is a bacterial pathogen that is the causative agent of erysipelas in a variety of animals, including swine, emus, turkeys, muskox, caribou, moose, and humans. This study aims to investigate the population structure and genomic features of Australian isolates of E. rhusiopathiae in the Australian pig industry and compare them to the broader scope of isolates worldwide. A total of 178 isolates (154 Australian, seven vaccine isolates, six international isolates, and 11 of unknown origin) in this study were screened against an MLST scheme and publicly available reference isolates, identifying 59 new alleles, with isolates separating into two main single locus variant groups. Investigation with BLASTn revealed the presence of the spaA gene in 171 (96%) of the isolates, with three main groups of SpaA protein sequences observed amongst the isolates. Novel SpaA protein sequences, categorised here as group 3 sequences, consisted of two sequence types forming separate clades to groups 1 and 2, with amino acid variants at positions 195 (D/A), 303 (G/E) and 323(P/L). In addition to the newly identified groups, five new variant positions were identified, 124 (S/N), 307 (Q/R), 323 (P/L), 379 (M/I), and 400 (V/I). Resistance screening identified genes related to lincomycin, streptomycin, erythromycin, and tetracycline resistance. Of the 29 isolates carrying these resistance genes, 82% belonged to SpaA group 2-N101S (n = 22) or 2-N101S-I257L (n = 2). In addition, 79% (n = 23) of these 29 isolates belonged to MLST group ST 5. Our results illustrate that Australia appears to have a unique diversity of E. rhusiopathiae isolates in pig production industries within the wider global context of isolates.