<p>Antibiotic resistance is often linked to human use of antibiotics. However, antibiotics and antibiotic biosynthetic pathways have been evolving for millions of years suggesting that antibiotic resistance is an ancient phenomenon. As of now, there has been no systematic survey of environmental microbes proven to exist in the absence of human influence and Lechuguilla cave offers such environment<em>. </em> Resistance diversity in strains isolated from this cave was analyzed by a phenotypic screen against a panel of 26 different antibiotics. Resistant strains were further investigated through determination of minimal inhibitory concentration (MIC) and inactivation studies. Of particular interest was strain LC044 (<em>Brachybacterium paraconglomeratum</em>), observed to inactivate macrolide antibiotics by phosphorylation. Genome sequencing and bioinformatics (BLAST analysis) identified a putative macrolide phosphotransferase (MPH) in strain LC044 and biochemical characterization of the purified recombinant protein confirmed its macrolide inactivating properties. To investigate if characterized MPH was unique to cave isolate, available terrestrial <em>Brachybacterium faecium</em> DSM 4810 genome was mined for presence of MPH-like protein. The top hit to the MPH from LC044 (a protein with 282 amino acids and 72% identity) was heterologously expressed and purified. Complete biochemical analysis of this enzyme revealed (i) MPH-activity, despite its annotation as aminoglycoside phosphotransferase (APH), and (ii) no significant differences in substrate specificities or kinetic parameters between these two enzymes suggesting that these two enzymes were equally effective resistance enzymes. This work highlights the prevalence of antibiotic resistance in a pristine, cave ecosystem and provides further support for the theory that antibiotic resistance is everywhere. Furthermore, the <em>mph</em> resistance determinant found in cave isolate and closely related terrestrial isolate show homology to clinical<em> mph</em> genes, suggesting that environmental <em>mph</em> genes could have served as reservoir of clinical determinants.</p> / Master of Science (MSc)

Enterococcus faecalis is a Gram-positive bacterial species that is typically a member of the human gastrointestinal tract microbiota. However, E. faecalis is also a nosocomial pathogen, which is involved in urinary tract infections, soft tissue infections and endocarditis. In recent times, the occurrence of antibiotic resistance has complicated the treatment of these infections. One of the major differences between commensal and pathogenic strains of E. faecalis is that pathogens contain multiple mobile elements such as plasmids, transposons and integrative conjugative elements (ICE). These elements allow for the acquisition and transfer of virulence factors and resistance genes. Conjugative plasmids are a class of plasmids present in E. faecalis whose transfer to host cells is induced by a small pheromone peptide, cCF10 (LVTLVFV). This peptide is initially encoded as a 22-amino acid precursor (pre-cCF10) from the signal sequence of the chromosomal ccfA gene and is then proteolytically cleaved by signal peptidase II and Eep. Once pCF10 has been transferred a host E. faecalis cell, it is exceptionally stable. A replicon clone is maintained in greater than 85% of host cells over 100 generations in the absence of selection, suggesting the stability of pCF10 is intrinsic to the replicon. Three unique features of the replication initiation protein PrgW may contribute to this stability: (a) the interaction of PrgW with pre-cCF10, (b) disulfide bond formation at three conserved cysteines (C78, C275, and C307) in PrgW, and (c) processing of the nascent PrgW protein. Replication initiation proteins associated with theta replicons, such as pCF10, are often self-contained units. To initiate plasmid replication, the replication initiation protein (PrgW) binds to direct repeats (oriV) in its own coding sequence (prgW). In silico analysis of PrgW suggests the existence of three distinct domains within the protein. The first 122 amino acids are homologous to a conserved domain present in related replication initiation proteins, which includes a Helix-Turn-Helix (HTH) DNA binding domain. This suggests that this domain of PrgW has a DNA-binding function and binds to the oriV site in prgW. The following 61 amino acids are not similar to any known sequence, and are encoded by the DNA sequence containing the direct repeats in the oriV site. This domain may or may not have a distinct function. The remaining sequence forms a domain that contains cysteines C275 and C307, and is also similar to no known structure. It is hypothesized that this domain is related to the stability of pCF10. C307 appears to be critical, as previous experiments indicate that its mutation alone affects plasmid stability. Secondary structure analysis of this domain revealed the presence of multiple alpha-helices that contain distinct hydrophobic domains that possibly contribute to pre-cCF10 binding and PrgW tertiary structure. The positions of the conserved cysteines within these alpha-helices may stabilize a hydrophobic binding pocket that could potentially facilitate interaction with pre-cCF10. PrgW has a predicted molecular weight of 38.6 KDa and can be detected in Western blots as a band with an apparent approximate molecular weight (mw) of 36,000. Previous data from our lab indicates that, when overexpressed in E. faecalis, four bands of PrgW are present with observed molecular weights of 40,000, 36,000, 24,000 and 18,000. Time course experiments revealed that the 40,000 mw form is converted to a 36,000 mw form independent. The 40,000 mw form is unstable (with a complete turnover in 30 minutes) while the 36,000 mw form has a half-life of greater than 4 hours. The 24,000 mw band does not have a DNA binding motif and is likely a turnover product. When the three conserved cysteines (and only cysteines) in PrgW are replaced with alanine, the 40,000 mw form is still processed to the 36,000 mw form. However, the cysteine to alanine mutants accumulate the 36,000 mw form. / Microbiology and Immunology

Metagenomics has become a reliable tool for the analysis of the microbial diversity and the molecular mechanisms carried out by microbial communities. By the use of next generation sequencing, metagenomic studies can generate millions of short sequencing reads that are processed by computational tools. However, with the rapid adoption of metagenomics a large amount of data has been generated. This situation requires the development of computational tools and pipelines to manage the data scalability, accessibility, and performance. In this thesis, several strategies varying from command line, web-based platforms to machine learning have been developed to address these computational challenges. Interpretation of specific information from metagenomic data is especially a challenge for environmental samples as current annotation systems only offer broad classification of microbial diversity and function. Therefore, I developed MetaStorm, a public web-service that facilitates customization of computational analysis for metagenomic data. The identification of antibiotic resistance genes (ARGs) from metagenomic data is carried out by searches against curated databases producing a high rate of false negatives. Thus, I developed DeepARG, a deep learning approach that uses the distribution of sequence alignments to predict over 30 antibiotic resistance categories with a high accuracy. Curation of ARGs is a labor intensive process where errors can be easily propagated. Thus, I developed ARGminer, a web platform dedicated to the annotation and inspection of ARGs by using crowdsourcing. Effective environmental monitoring tools should ideally capture not only ARGs, but also mobile genetic elements and indicators of co-selective forces, such as metal resistance genes. Here, I introduce NanoARG, an online computational resource that takes advantage of the long reads produced by nanopore sequencing technology to provide insights into mobility, co-selection, and pathogenicity. Sequence alignment has been one of the preferred methods for analyzing metagenomic data. However, it is slow and requires high computing resources. Therefore, I developed MetaMLP, a machine learning approach that uses a novel representation of protein sequences to perform classifications over protein functions. The method is accurate, is able to identify a larger number of hits compared to sequence alignments, and is >50 times faster than sequence alignment techniques. / Doctor of Philosophy / Antimicrobial resistance (AMR) is one of the biggest threats to human public health. It has been estimated that the number of deaths caused by AMR will surpass the ones caused by cancer on 2050. The seriousness of these projections requires urgent actions to understand and control the spread of AMR. In the last few years, metagenomics has stand out as a reliable tool for the analysis of the microbial diversity and the AMR. By the use of next generation sequencing, metagenomic studies can generate millions of short sequencing reads that are processed by computational tools. However, with the rapid adoption of metagenomics, a large amount of data has been generated. This situation requires the development of computational tools and pipelines to manage the data scalability, accessibility, and performance. In this thesis, several strategies varying from command line, web-based platforms to machine learning have been developed to address these computational challenges. In particular, by the development of computational pipelines to process metagenomics data in the cloud and distributed systems, the development of machine learning and deep learning tools to ease the computational cost of detecting antibiotic resistance genes in metagenomic data, and the integration of crowdsourcing as a way to curate and validate antibiotic resistance genes.

The use of antibiotics in dairy cattle has the potential to stimulate the development and subsequent fecal dissemination of antibiotic resistance genes (ARGs) in bacteria. The objectives were to use metagenomic techniques to evaluate the effect of antibiotic treatment on ARG prevalence in the fecal microbiome of the dairy cow and to determine the temporal excretion pattern of ARGs. Twelve Holstein cows were assigned to one of four antibiotic treatments: control, pirlimycin, ceftiofur, or cephapirin. Fecal samples were collected on d -1, 1, 3, 5, 7, 14, 21, and 28. Samples were freeze-dried and subjected to DNA extraction followed by Illumina paired-end HiSeq sequencing and quantitative polymerase chain reaction (qPCR). Illumina sequences were analyzed using MG-RAST and the Antibiotic Resistance Gene Database (ARDB) via BLAST. Abundance of ampC, ermB, tetO, tetW, and 16S rRNA genes were determined using qPCR. All data were statistically analyzed with PROC GLIMMIX in SAS. Antibiotic treatment resulted in a shift in bacterial cell functions. Sequences associated with 'resistance to antibiotics and toxic compounds' were higher in ceftiofur-treated cows than control cows. Ceftiofur-treated cows had a higher abundance of 𝛽-lactam and multidrug resistance sequences than control cows. There was no effect of treatment or day on fecal tetO and ermB excretion. The relative abundances of tetW and ampC were higher on d 3 post-treatment than d 5 and d 28. In conclusion, antibiotic use in dairy cattle shifted bacterial cell functions and temporarily increased antibiotic resistance in the fecal microbiome. / Master of Science

Exposure to antibiotics has the potential to increase the incidence and proliferation of antibiotic resistance genes (ARG) in the gut and fecal microbiome. Non-saleable, antibiotic-containing milk from cows treated with antibiotics (waste milk) is commonly fed to dairy calves but the effects of ingestion of antibiotics at an early age on the gut microbiome and the development of ARG in the naive gut are not well understood. Pirlimycin, a lincosamide antibiotic acting against Gram positive bacteria through inhibiting protein synthesis by binding to the 50S ribosome, is commonly used as mastitis therapy. Lincosamides are also considered highly important in human medicine, often used against Staphylococcus aureus and Clostridium difficile infections. Emerging microbial resistance to pirlimycin is of concern for both animal and human health. The objective of this study was to determine the effect of early lincosamide antibiotic exposure on the abundance of ARG in feces of milk-fed calves. Eight female Holstein calves were blocked by age, paired by block, and randomly assigned to pasteurized whole milk (control; n = 4) or milk containing 0.2 mg/L of pirlimycin (treatment; n = 4). Calves were enrolled after receiving two colostrum feedings and were fed 5.68 L of pasteurized whole milk, treatment, or control, divided into two daily feedings, from d 1 to d 50 of age. After weaning calves were fed non-medicated starter grain ad libitum. Fecal samples were collected weekly until 85 d of age and freeze-dried. DNA was extracted using QiaAmp® Fast DNA Stool Mini Kit and qPCR was used to quantify the absolute abundance (gene copies/g of wet feces) and relative abundance (gene copies/copies of 16S rRNA genes) of erm(B), tet(O), tet(W) and 16S rRNA genes. Data was analyzed using PROC GLIMMIX in SAS. Abundance of 16S rRNA genes, tet(O) and tet(W) were not different between control and pirlimycin-fed calves nor were the relative abundance of tet(O) (mean = 0.050 tet(O) copies/16S rRNA genes) or tet(W) (0.561 tet(W) copies/16S rRNA genes). While abundance of erm(B) was higher in pirlimycin-fed calves compared to control calves (6.46 and 6.04 log gene copies/g wet feces; P = 0.04) the relative abundance of erm(B) (0.273 gene copies/16S rRNA genes) in feces of calves was not influenced by treatment. There was an effect of day (P < 0.10) for absolute abundance of tet(O), tet(W), and erm(B) indicating that the levels change with time as the fecal microbiome develops. This study suggests that feeding pirlimycin-containing non-saleable milk to growing calves may increase environmental loading of erm(B), which codes for resistance to highly important macrolide and lincosamide antibiotics. Additional research is needed on effects of feeding waste milk to calves on other fecal ARG and on the post-excretion and post-application fate of these genes. / Master of Science

Manure is a commonly used soil fertilizer, but there are concerns that this practice could affect the spread of antibiotic resistance genes (ARGs) from farm to fork. A microcosm-scale study evaluated the effect of prior antibiotic use (manure-based soil amendments generated from dairy and beef cattle with or without antibiotic administration), composting, and soil type on the quantity of ARGs and the microbial community composition of dairy and beef manure applied soil. ARGs were analyzed through novel metagenomic techniques and quantitative polymerase chain reaction of sul1, tet(W), and 16S rRNA gene, while the microbial community composition was determined via 16S rRNA amplicon sequencing. The results indicated that while prior antibiotic administration elevated the relative abundance of ARGs and changed the microbial community of raw manure applied soils, composting reduced this effect. However, compost applied soils still had a higher relative abundance of ARGs than the unamended soils and occasionally soil applied with raw manure of untreated cattle. Soil type may be a mediating factor as there were differences observed between the three soil types (sandy loam, silty clay loam, and silty loam) with sandy loam amended soils often having the least attenuation of ARGs. As the relative abundance of ARGs was still elevated and the microbial community composition still significantly different from the unamended soils after 120 days, these results suggest that 120 days is not a long enough waiting period between biological soil amendments and crop harvest for ARG dissipation. / Master of Science

Wastewater treatment plants (WWTPs) receive a confluence of sewage containing antibiotics, antibiotic resistant bacteria, antibiotic resistance genes (ARGs), and pathogens, thus serving as key point of interest for the surveillance of antibiotic resistance (AR) dissemination. This thesis advances knowledge about the fate of AR indicators throughout treatment and reuse. The field study informs approaches for monitoring AR at a WWTP by characterizing the resistome (i.e., full profile of ARGs) and microbiome across eight sampling events via metagenomic sequencing, complemented by antibiotic data. The WWTP significantly reduced the total load of ARGs and antibiotics, although correlations between ARGs and antibiotics were generally weak. Quantitative polymerase chain reaction was applied to validate the quantitative capacity of metagenomics, whereby we found strong correlations. The influent and effluent to the WWTP were remarkably stable with time, providing further insight into the sampling frequency necessary for adequate surveillance. The laboratory study examined the effects of commonly applied disinfection processes (chlorination, chloramination, and ultraviolet irradiation [UV]) on the inactivation of antibiotic resistant pathogens and corresponding susceptible pathogens in recycled and potable water. Further, we evaluated their regrowth following disinfection by simulating distribution. Acinetobacter baumannii, an environmental opportunistic pathogen, regrew especially well following UV disinfection, although not when a disinfectant residual was present. Enterococcus faecium, a fecal pathogen, did not regrow following any disinfection process. There were no significant differences between water types. The findings of this study emphasize a need to move beyond the framework of assessing treatment efficacy based on the attenuation of fecal pathogens. / Master of Science / Wastewater treatment plants (WWTPs) have traditionally been designed and further enhanced to minimize environmental contamination caused by solid waste, fecal pathogens, nutrients (e.g., nitrogen), and organic matter. However, treatment has not been optimized to remove the contaminants of emerging concern (CECs) investigated in this thesis: antibiotic resistant bacteria (ARB), antibiotic resistance genes (ARGs), and antibiotics. WWTPs are key point of interest for local and global surveillance of antibiotic resistance as they can receive the aforementioned CECs (via human excretion or improper disposal) from various sources (e.g., residences, hospitals). Antibiotic resistant bacteria have caused 2.8 million infections and subsequently 35,000 deaths in the United States each year. Considering treated wastewater can serve as a route of exposure for humans, potential spread of antibiotic resistance by WWTPs is of high priority to mitigate from a public health perspective. In the first study utilizing a technology to assess the full complement of ARGs in a given sample, we observed that the total load of ARGs was removed by approximately 50% across wastewater treatment, on average; total antibiotic load exhibited a similar reduction. The second study demonstrated that antibiotic resistant environmental opportunistic pathogen (i.e., pathogens which take advantage of the "opportunity" to infect an immunocompromised host, especially thriving in low nutrient engineered systems), Acinetobacter baumannii, possesses the ability to regrow following disinfection in the absence of a disinfectant residual. In contrast, antibiotic resistant Enterococcus faecium, an opportunistic pathogen of fecal origin, was successfully inactivated and unable to regrow. The findings of this study emphasize a need to move beyond the framework of assessing treatment efficacy based on the attenuation of fecal pathogens.

The emergence of antibiotic resistance in riverine systems has become a growing issue worldwide. The use of watershed-scale models is popular with many other water quality issues but not in the case of antibiotic resistance. In this study, we introduce a watershed-scale bacteria fate and transport model to simulate antibiotic resistance in E. coli. This model was developed through amendments to an existing watershed-scale physically based hydrological model (SWAT), and the newly modified model was called SWAT-ARB. The SWAT-ARB model was employed in the receiving environment of a WWTP in the Adyar River basin, India. The SWAT-ARB model simulations of resistant fractions (resistant E. coli concentration/E. coli concentration) in stream water were analyzed by the flow levels with the application of a range of parameter values. It is concluded that the model can be used to test prevailing hypotheses and evaluate the current state of knowledge. For instance, model simulations suggest that the influx of ARB can be a primary driver of antibiotic resistance in rivers compared to ambient antibiotic concentrations. We used the SWAT-ARB model to quantify the impact of climate change on antibiotic resistance. Six climate models were used to obtain the future climates in two distinct scenarios. The model was applied to three watersheds as Adyar basin- India, Crab Creek basin- USA, and upper Viskan basin- Sweden. It was concluded that temperature increase may greatly affect the colder climates (Crab Creek and Viskan) with higher simulated resistant fractions. In case of Adyar basin, resistant fractions are alleviated in high flow conditions, while aggravated in low flow conditions. / Doctor of Philosophy / The antibiotic resistance occurs when bacteria no longer responds to antibiotics. Hence, the diseases that caused by resistant bacteria are harder to treat. These antibiotic resistant bacteria end up in our rivers because of our heavy use of antibiotics in human and animal treatments. Thus, the spread of antibiotic resistance has become a water quality issue in the rivers worldwide. Scientists generally use computer models to understand water quality issues in rivers. These computer models are important because of high cost of monitoring and their use in finding how environment works. Up to the date of this publication, there is no sophisticated enough model to simulate antibiotic resistance in rivers. Hence, we created a river basin scale model to simulate antibiotic resistance. We found that the influx of ARB can be a primary driver of antibiotic resistance in rivers compared to ambient antibiotic concentrations. The model was applied to three watersheds as Adyar basin- India, Crab Creek basin- USA, and upper Viskan basin- Sweden. It was concluded that temperature increase may greatly affect the colder climates (Crab Creek and Viskan) with higher antibiotic resistant bacteria compared to susceptible bacteria.

Antibiotic resistance is currently a global threat spanning clinical, environmental, and geopolitical research domains. The environment is increasingly recognized as a key node in the spread of antibiotic resistance genes (ARGs), which confer antibiotic resistance to bacteria. Detecting ARGs in the environment is the first step in monitoring and controlling antibiotic resistance. In recent years, next-generation sequencing of environmental samples (metagenomic sequencing data) has become a prolific tool for the field of surveillance. Metagenomic data are nucleic acid sequences, or nucleotides, of environmental samples. Metagenomic sequencing data has been used over the years to detect and analyze ARGs. An intriguing instance of ARGs is the multi-functional ARG, where one ARG encodes two or more different antibiotic resistance functions. Multi-functional ARGs provide resistance to two or more antibiotics, thus should have evolutionary advantage over ARGs with resistance to single antibiotic. However, there is no tool readily available to detect these multi-functional ARGs in metagenomic data. In this study, we develop a computational pipeline to detect multi-functional ARGs in metagenomic data. The pipeline takes raw metagenomic data as the input and generates a list of potential multi-functional ARGs. A plot for each potential multi-functional ARG is also created, showing the location of the multi-functionalities in the sequence and the sequencing coverage level. We collected samples from three different sources: influent samples of a wastewater treatment plant, hospital wastewater samples, and reclaimed water samples, ran the pipeline, and identified 19, 57, and 8 potentially bi-functional ARGs in each source, respectively. Manual inspection of the results identified three most likely bi-functional ARGs. Interestingly, one bi-functional ARG, encoding both aminoglycoside and tetracycline resistance, appeared in all three data sets, indicating its prevalence in different environments. As the amount of antibiotics keeps increasing in the environment, multi-functional ARGs might become more and more common. The pipeline will be a useful computational tool for initial screening and identification of multi-functional ARGs in metagenomic data. / Master of Science / Antibiotics are the drug to fight against the infection of bacteria. Since the first antibiotic was discovered in 1928, many antibiotic drugs have been developed. At the same time, scientists discovered many genes responsible for the resistance of antibiotic drugs. Nowadays, antibiotic resistance is a global threat. Detecting antibiotic resistance genes in the environment is the first step toward monitoring and controlling antibiotic resistance. In recent years, next-generation sequencing has been widely used to get the DNA sequence from the environment. Metagenomics analysis has been used over the years to detect and analyze ARGs. In the literature, it has been reported that a single gene could carry two parts of sequences corresponding to two different ARGs, thus conferring resistance to two different antibiotics. This fusion might have some evolutionary advantages. In this study, we developed a novel computational tool to detect multi-functional ARGs. We collected data from three sources: the treatment plant water, the hospital wastewater, and the reclaimed water, and identified 19, 57, and 8 potential bi-functional ARGs in each source, respectively. After we manually inspected the result, we found three most likely bi-functional ARGs. We also found one bi-functional ARG that appears in all three datasets. The gene is responsible for aminoglycoside and tetracycline resistance. The tool will serve as the initial screening step to detect multi-functional ARGs.

Mycobacterium avium isolates resistant to clarithromycin and azithromycin have been recovered from patients undergoing antibiotic therapy. Comparison of DNA fingerprints of sensitive and resistant isolates showed that resistance resulted from mutation of the original, sensitive isolate in five of seven patients. In the other two patients, the clarithromycin-resistant isolates were unrelated to the sensitive isolate, suggesting that the resistant isolate resulted from either superinfection or selection of a resistant strain from a polyclonal population. Investigation of the mechanisms of clarithromycin and azithromycin resistance in M. avium showed that high-level resistance resulted from a point mutation at position A-2058 in the 23S rRNA. Based on this finding, a rapid screen for clarithromycin-resistance in M. avium was developed based on PCR. Twenty-three clinical isolates were analyzed, seven of which were clarithromycin-resistant. The target product was amplified only in clarithromycin-resistant strains, all of which had mutations at position 2058. A polyuridylic acid (poly U)-dependent in vitro translation system from M. avium was developed to investigate the effect of antibiotics on protein synthesis. Clarithromycin was an effective inhibitor of protein synthesis in cell-free extracts of a susceptible M. avium strain, whereas a high-level resistant strain was less susceptible to clarithromycin in vitro. Mixtures of extracts from sensitive and resistant strains showed a pattern of clarithromycin inhibition similar to the resistant strain, suggesting that resistance may be dominant in partial diploids. Three M. avium strains exhibiting step-wise, intermediate resistance to azithromycin were characterized in comparison to the sensitive parent. All strains were similar in hydrophobicity, growth medium requirements, and growth response to temperature. The azithromycin-resistant strains were resistant to several unrelated agents, including ciprofloxacin, rifabutin, and ethidium bromide. Addition of carbonyl cyanide m-chlorophenylhydrazone (CCCP) did not lower minimal inhibitory concentrations (MICs) for ciprofloxacin or ethidium bromide. Cell-free extracts of the strains were as sensitive to azithromycin in vitro as the parent strain. The results rule out inactivation, efflux, and mutations in the target as resistance mechanisms, and suggest intermediate resistance may be due to altered permeability of the cell wall or membrane. / Ph. D.