The contribution of DNA repair pathways to Staphylococcus aureus fitness and fidelity during nitric oxide stress.

Staphylococcus aureus is a major human pathogen that causes a variety of illnesses, ranging from minor skin and soft tissue infections to more severe systemic infections. Although the primary host immune response can typically clear bacterial infections, S. aureus is uniquely resistant to inflammation. For instance, our laboratory has determined that S. aureus is highly resistant to nitric oxide (NO⋅), an important component of the innate immune response that plays a role in both immunomodulatory and antibacterial processes. Additionally, NO⋅ and its derivatives can cause damage to S. aureus DNA, more specifically, deamination and/or oxidation of DNA bases; however, regulation and repair mechanisms of DNA in S. aureus are understudied. Thus, we hypothesize that several DNA repair mechanisms may account for the replication fidelity of S. aureus and may contribute to fitness in the presence of NO⋅. Here, we show the role of several DNA repair mechanisms in S. aureus. More specifically, we found that recombinational repair genes recJ, recG, and polA may play a role in the repair of NO⋅-induced replication fork collapses. We also show the role of the base excision repair pathway protein, MutY, in reducing NO⋅-mediated mutagenesis. Overall, our results suggest that NO⋅ leads to DNA damage, which subsequently induces the activity of several DNA repair pathways, contributing to the replication fidelity and fitness of S. aureus.IMPORTANCEPathogenic bacteria must evolve various mechanisms in order to evade the host immune response that they are infecting. One aspect of the primary host immune response to an infection is the production of an inflammatory effector component, nitric oxide (NO⋅). Staphylococcus aureus has uniquely evolved a diverse array of strategies to circumvent the inhibitory activity of nitric oxide. One such mechanism by which S. aureus has evolved allows the pathogen to survive and maintain its genomic integrity in this environment. For instance, here, our results suggest that S. aureus employs several DNA repair pathways to ensure replicative fitness and fidelity under NO⋅ stress. Thus, our study presents evidence of an additional strategy that allows S. aureus to evade the cytotoxic effects of host NO⋅.

TheiaEuk: a species-agnostic bioinformatics workflow for fungal genomic characterization.

Introduction: The clinical incidence of antimicrobial-resistant fungal infections has dramatically increased in recent years. Certain fungal pathogens colonize various body cavities, leading to life-threatening bloodstream infections. However, the identification and characterization of fungal isolates in laboratories remain a significant diagnostic challenge in medicine and public health. Whole-genome sequencing provides an unbiased and uniform identification pipeline for fungal pathogens but most bioinformatic analysis pipelines focus on prokaryotic species. To this end, TheiaEuk_Illumina_PE_PHB (TheiaEuk) was designed to focus on genomic analysis specialized to fungal pathogens. Methods: TheiaEuk was designed using containerized components and written in the workflow description language (WDL) to facilitate deployment on the cloud-based open bioinformatics platform Terra. This species-agnostic workflow enables the analysis of fungal genomes without requiring coding, thereby reducing the entry barrier for laboratory scientists. To demonstrate the usefulness of this pipeline, an ongoing outbreak of C. auris in southern Nevada was investigated. We performed whole-genome sequence analysis of 752 new C. auris isolates from this outbreak. Furthermore, TheiaEuk was utilized to observe the accumulation of mutations in the FKS1 gene over the course of the outbreak, highlighting the utility of TheiaEuk as a monitor of emerging public health threats when combined with whole-genome sequencing surveillance of fungal pathogens. Results: A primary result of this work is a curated fungal database containing 5,667 unique genomes representing 245 species. TheiaEuk also incorporates taxon-specific submodules for specific species, including clade-typing for Candida auris (C. auris). In addition, for several fungal species, it performs dynamic reference genome selection and variant calling, reporting mutations found in genes currently associated with antifungal resistance (FKS1, ERG11, FUR1). Using genome assemblies from the ATCC Mycology collection, the taxonomic identification module used by TheiaEuk correctly assigned genomes to the species level in 126/135 (93.3%) instances and to the genus level in 131/135 (97%) of instances, and provided zero false calls. Application of TheiaEuk to actual specimens obtained in the course of work at a local public health laboratory resulted in 13/15 (86.7%) correct calls at the species level, with 2/15 called at the genus level. It made zero incorrect calls. TheiaEuk accurately assessed clade type of Candida auris in 297/302 (98.3%) of instances. Discussion: TheiaEuk demonstrated effectiveness in identifying fungal species from whole genome sequence. It further showed accuracy in both clade-typing of C. auris and in the identification of mutations known to associate with drug resistance in that organism.

Accelerating bioinformatics implementation in public health.

We have adopted an open bioinformatics ecosystem to address the challenges of bioinformatics implementation in public health laboratories (PHLs). Bioinformatics implementation for public health requires practitioners to undertake standardized bioinformatic analyses and generate reproducible, validated and auditable results. It is essential that data storage and analysis are scalable, portable and secure, and that implementation of bioinformatics fits within the operational constraints of the laboratory. We address these requirements using Terra, a web-based data analysis platform with a graphical user interface connecting users to bioinformatics analyses without the use of code. We have developed bioinformatics workflows for use with Terra that specifically meet the needs of public health practitioners. These Theiagen workflows perform genome assembly, quality control, and characterization, as well as construction of phylogeny for insights into genomic epidemiology. Additonally, these workflows use open-source containerized software and the WDL workflow language to ensure standardization and interoperability with other bioinformatics solutions, whilst being adaptable by the user. They are all open source and publicly available in Dockstore with the version-controlled code available in public GitHub repositories. They have been written to generate outputs in standardized file formats to allow for further downstream analysis and visualization with separate genomic epidemiology software. Testament to this solution meeting the requirements for bioinformatic implementation in public health, Theiagen workflows have collectively been used for over 5 million sample analyses in the last 2 years by over 90 public health laboratories in at least 40 different countries. Continued adoption of technological innovations and development of further workflows will ensure that this ecosystem continues to benefit PHLs.

The use of whole-genome sequencing and development of bioinformatics to monitor overlapping outbreaks of Candida auris in southern Nevada.

A Candida auris outbreak has been ongoing in Southern Nevada since August 2021. In this manuscript we describe the sequencing of over 200 C. auris isolates from patients at several facilities. Genetically distinct subgroups of C. auris were detected from Clade I (3 distinct lineages) and III (1 lineage). Open-source bioinformatic tools were developed and implemented to aid in the epidemiological investigation. The work herein compares three methods for C. auris whole genome analysis: Nullarbor, MycoSNP and a new pipeline TheiaEuk. We also describe a novel analysis method focused on elucidating phylogenetic linkages between isolates within an ongoing outbreak. Moreover, this study places the ongoing outbreaks in a global context utilizing existing sequences provided worldwide. Lastly, we describe how the generated results were communicated to the epidemiologists and infection control to generate public health interventions.

Pathogen genomics in public health laboratories: successes, challenges, and lessons learned from California's SARS-CoV-2 Whole-Genome Sequencing Initiative, California COVIDNet.

The capacity for pathogen genomics in public health expanded rapidly during the coronavirus disease 2019 (COVID-19) pandemic, but many public health laboratories did not have the infrastructure in place to handle the vast amount of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) sequence data generated. The California Department of Public Health, in partnership with Theiagen Genomics, was an early adopter of cloud-based resources for bioinformatics and genomic epidemiology, resulting in the creation of a SARS-CoV-2 genomic surveillance system that combined the efforts of more than 40 sequencing laboratories across government, academia and industry to form California COVIDNet, California's SARS-CoV-2 Whole-Genome Sequencing Initiative. Open-source bioinformatics workflows, ongoing training sessions for the public health workforce, and automated data transfer to visualization tools all contributed to the success of California COVIDNet. While challenges remain for public health genomic surveillance worldwide, California COVIDNet serves as a framework for a scaled and successful bioinformatics infrastructure that has expanded beyond SARS-CoV-2 to other pathogens of public health importance.

The Nutritional Environment Is Sufficient To Select Coexisting Biofilm and Quorum Sensing Mutants of Pseudomonas aeruginosa.

The evolution of bacterial populations during infections can be influenced by various factors including available nutrients, the immune system, and competing microbes, rendering it difficult to identify the specific forces that select on evolved traits. The genomes of Pseudomonas aeruginosa isolated from the airways of people with cystic fibrosis (CF), for example, have revealed commonly mutated genes, but which phenotypes led to their prevalence is often uncertain. Here, we focus on effects of nutritional components of the CF airway on genetic adaptations by P. aeruginosa grown in either well-mixed (planktonic) or biofilm-associated conditions. After only 80 generations of experimental evolution in a simple medium with glucose, lactate, and amino acids, all planktonic populations diversified into lineages with mutated genes common to CF infections: morA, encoding a regulator of biofilm formation, or lasR, encoding a quorum sensing regulator that modulates the expression of virulence factors. Although mutated quorum sensing is often thought to be selected in vivo due to altered virulence phenotypes or social cheating, isolates with lasR mutations demonstrated increased fitness when grown alone and outcompeted the ancestral PA14 strain. Nonsynonymous SNPs in morA increased fitness in a nutrient concentration-dependent manner during planktonic growth and surprisingly also increased biofilm production. Populations propagated in biofilm conditions also acquired mutations in loci associated with chronic infections, including lasR and cyclic di-GMP regulators roeA and wspF. These findings demonstrate that nutrient conditions and biofilm selection are sufficient to select mutants with problematic clinical phenotypes including increased biofilm and altered quorum sensing. IMPORTANCE Pseudomonas aeruginosa produces dangerous chronic infections that are known for their rapid diversification and recalcitrance to treatment. We performed evolution experiments to identify adaptations selected by two specific aspects of the CF respiratory environment: nutrient levels and surface attachment. Propagation of P. aeruginosa in nutrients present within the CF airway was sufficient to drive diversification into subpopulations with identical mutations in regulators of biofilm and quorum sensing to those arising during infection. Thus, the adaptation of opportunistic pathogens to nutrients found in the host may select mutants with phenotypes that complicate treatment and clearance of infection.

Quantitative mapping of mRNA 3' ends in Pseudomonas aeruginosa reveals a pervasive role for premature 3' end formation in response to azithromycin.

Pseudomonas aeruginosa produces serious chronic infections in hospitalized patients and immunocompromised individuals, including patients with cystic fibrosis. The molecular mechanisms by which P. aeruginosa responds to antibiotics and other stresses to promote persistent infections may provide new avenues for therapeutic intervention. Azithromycin (AZM), an antibiotic frequently used in cystic fibrosis treatment, is thought to improve clinical outcomes through a number of mechanisms including impaired biofilm growth and quorum sensing (QS). The mechanisms underlying the transcriptional response to AZM remain unclear. Here, we interrogated the P. aeruginosa transcriptional response to AZM using a fast, cost-effective genome-wide approach to quantitate RNA 3' ends (3pMap). We also identified hundreds of P. aeruginosa genes with high incidence of premature 3' end formation indicative of riboregulation in their transcript leaders using 3pMap. AZM treatment of planktonic and biofilm cultures alters the expression of hundreds of genes, including those involved in QS, biofilm formation, and virulence. Strikingly, most genes downregulated by AZM in biofilms had increased levels of intragenic 3' ends indicating premature transcription termination, transcriptional pausing, or accumulation of stable intermediates resulting from the action of nucleases. Reciprocally, AZM reduced premature intragenic 3' end termini in many upregulated genes. Most notably, reduced termination accompanied robust induction of obgE, a GTPase involved in persister formation in P. aeruginosa. Our results support a model in which AZM-induced changes in 3' end formation alter the expression of central regulators which in turn impairs the expression of QS, biofilm formation and stress response genes, while upregulating genes associated with persistence.

Parallel Evolution of Tobramycin Resistance across Species and Environments.

Different species exposed to a common stress may adapt by mutations in shared pathways or in unique systems, depending on how past environments have molded their genomes. Understanding how diverse bacterial pathogens evolve in response to an antimicrobial treatment is a pressing example of this problem, where discovery of molecular parallelism could lead to clinically useful predictions. Evolution experiments with pathogens in environments containing antibiotics, combined with periodic whole-population genome sequencing, can be used to identify many contending routes to antimicrobial resistance. We separately propagated two clinically relevant Gram-negative pathogens, Pseudomonas aeruginosa and Acinetobacter baumannii, in increasing concentrations of tobramycin in two different environments each: planktonic and biofilm. Independently of the pathogen, the populations adapted to tobramycin selection by parallel evolution of mutations in fusA1, encoding elongation factor G, and ptsP, encoding phosphoenolpyruvate phosphotransferase. As neither gene is a direct target of this aminoglycoside, mutations to either are unexpected and underreported causes of resistance. Additionally, both species acquired antibiotic resistance-associated mutations that were more prevalent in the biofilm lifestyle than in the planktonic lifestyle; these mutations were in electron transport chain components in A. baumannii and lipopolysaccharide biosynthesis enzymes in P. aeruginosa populations. Using existing databases, we discovered site-specific parallelism of fusA1 mutations that extends across bacterial phyla and clinical isolates. This study suggests that strong selective pressures, such as antibiotic treatment, may result in high levels of predictability in molecular targets of evolution, despite differences between organisms' genetic backgrounds and environments.IMPORTANCE The rise of antimicrobial resistance is a leading medical threat, motivating efforts to forecast both its evolutionary dynamics and its genetic causes. Aminoglycosides are a major class of antibiotics that disrupt translation, but resistance may occur by a number of mechanisms. Here, we show the repeated evolution of resistance to the aminoglycoside tobramycin in both P. aeruginosa and A. baumannii via mutations in fusA1, encoding elongation factor G, and ptsP, encoding the nitrogen-specific phosphotransferase system. Laboratory evolution and whole-population genome sequencing were used to identify these targets, but mutations at identical amino acid positions were also found in published genomes of diverse bacterial species and clinical isolates. We also identified other resistance mechanisms associated with growth in biofilms that likely interfere with drug binding or uptake. Characterizing the evolution of multiple species in the presence of antibiotics can identify new, repeatable causes of resistance that may be predicted and counteracted by alternative treatment.

Evolutionary pathways to antibiotic resistance are dependent upon environmental structure and bacterial lifestyle.

Bacterial populations vary in their stress tolerance and population structure depending upon whether growth occurs in well-mixed or structured environments. We hypothesized that evolution in biofilms would generate greater genetic diversity than well-mixed environments and lead to different pathways of antibiotic resistance. We used experimental evolution and whole genome sequencing to test how the biofilm lifestyle influenced the rate, genetic mechanisms, and pleiotropic effects of resistance to ciprofloxacin in Acinetobacter baumannii populations. Both evolutionary dynamics and the identities of mutations differed between lifestyle. Planktonic populations experienced selective sweeps of mutations including the primary topoisomerase drug targets, whereas biofilm-adapted populations acquired mutations in regulators of efflux pumps. An overall trade-off between fitness and resistance level emerged, wherein biofilm-adapted clones were less resistant than planktonic but more fit in the absence of drug. However, biofilm populations developed collateral sensitivity to cephalosporins, demonstrating the clinical relevance of lifestyle on the evolution of resistance.

Coordination by Cdc42 of Actin, Contractility, and Adhesion for Melanoblast Movement in Mouse Skin.

The individual molecular pathways downstream of Cdc42, Rac, and Rho GTPases are well documented, but we know surprisingly little about how these pathways are coordinated when cells move in a complex environment in vivo. In the developing embryo, melanoblasts originating from the neural crest must traverse the dermis to reach the epidermis of the skin and hair follicles. We previously established that Rac1 signals via Scar/WAVE and Arp2/3 to effect pseudopod extension and migration of melanoblasts in skin. Here we show that RhoA is redundant in the melanocyte lineage but that Cdc42 coordinates multiple motility systems independent of Rac1. Similar to Rac1 knockouts, Cdc42 null mice displayed a severe loss of pigmentation, and melanoblasts showed cell-cycle progression, migration, and cytokinesis defects. However, unlike Rac1 knockouts, Cdc42 null melanoblasts were elongated and displayed large, bulky pseudopods with dynamic actin bursts. Despite assuming an elongated shape usually associated with fast mesenchymal motility, Cdc42 knockout melanoblasts migrated slowly and inefficiently in the epidermis, with nearly static pseudopods. Although much of the basic actin machinery was intact, Cdc42 null cells lacked the ability to polarize their Golgi and coordinate motility systems for efficient movement. Loss of Cdc42 de-coupled three main systems: actin assembly via the formin FMNL2 and Arp2/3, active myosin-II localization, and integrin-based adhesion dynamics.