Evolutionary genomic analyses of canine E. coli infections identify a relic capsular locus associated with resistance to multiple classes of antimicrobials.

Infections caused by antimicrobial-resistant Escherichia coli are the leading cause of death attributed to antimicrobial resistance (AMR) worldwide, and the known AMR mechanisms involve a range of functional proteins. Here, we employed a pan-genome wide association study (GWAS) approach on over 1,000 E. coli isolates from sick dogs collected across the US and Canada and identified a strong statistical association (empirical P < 0.01) of AMR, involving a range of antibiotics to a group 1 capsular (CPS) gene cluster. This cluster included genes under relaxed selection pressure, had several loci missing, and had pseudogenes for other key loci. Furthermore, this cluster is widespread in E. coli and Klebsiella clinical isolates across multiple host species. Earlier studies demonstrated that the octameric CPS polysaccharide export protein Wza can transmit macrolide antibiotics into the E. coli periplasm. We suggest that the CPS in question, and its highly divergent Wza, functions as an antibiotic trap, preventing antimicrobial penetration. We also highlight the high diversity of lineages circulating in dogs across all regions studied, the overlap with human lineages, and regional prevalence of resistance to multiple antimicrobial classes. IMPORTANCE: Much of the human genomic epidemiology data available for E. coli mechanism discovery studies has been heavily biased toward shiga-toxin producing strains from humans and livestock. E. coli occupies many niches and produces a wide variety of other significant pathotypes, including some implicated in chronic disease. We hypothesized that since dogs tend to share similar strains with their owners and are treated with similar antibiotics, their pathogenic isolates will harbor unexplored AMR mechanisms of importance to humans as well as animals. By comparing over 1,000 genomes with in vitro antimicrobial susceptibility data from sick dogs across the US and Canada, we identified a strong multidrug resistance association with an operon that appears to have once conferred a type 1 capsule production system.

Evolutionary genomic and bacteria GWAS analysis of Mycobacterium avium subsp. paratuberculosis and dairy cattle Johne's disease phenotypes.

Mycobacterium avium subsp. paratuberculosis (MAP) is the causative agent of Johne's disease in ruminants, which has important health consequences for dairy cattle. The Regional Dairy Quality Management Alliance (RDQMA) project is a multistate research program involving MAP isolates taken from three intensively studied commercial dairy farms in the northeastern United States, which emphasized longitudinal data collection of both MAP isolates and animal health in three regional dairy herds for a period of about 7 years. This paper reports the results of a pan-GWAS analysis involving 318 MAP isolates and dairy cow Johne's disease phenotypes, taken from these three farms. Based on our highly curated accessory gene count the pan-GWAS analysis identified several MAP genes associated with bovine Johne's disease phenotypes scored from these three farms, with some of the genes having functions suggestive of possible cause/effect relationships to these phenotypes. This paper reports a pan-genomic comparative analysis between MAP and Mycobacterium tuberculosis, assessing functional Gene Ontology category enrichments between these taxa. Finally, we also provide a population genomic perspective on the effectiveness of herd isolation, involving closed dairy farms, in preventing MAP inter-farm cross infection on a micro-geographic scale.IMPORTANCE Mycobacterium avium subsp. paratuberculosis (MAP) is the causative agent of Johne's disease in ruminants, which has important health consequences for dairy cattle, and enormous economic consequences for the dairy industry. Understanding which genes in this bacterium are correlated with key disease phenotypes can lead to functional experiments targeting these genes and ultimately lead to improved control strategies. This study represents a rare example of a prolonged longitudinal study of dairy cattle where the disease was measured and the bacteria were isolated from the same cows. The genome sequences of over 300 MAP isolates were analyzed for genes that were correlated with a wide range of Johne's disease phenotypes. A number of genes were identified that were significantly associated with several aspects of the disease and suggestive of further experimental follow-up.

Design Elements of Listeria Environmental Monitoring Programs in Food Processing Facilities: A Scoping Review of Research and Guidance Materials.

Occurrence of Listeria monocytogenes (Lm), the causative agent of listeriosis, in food processing facilities presents considerable challenges to food producers and food safety authorities. Design of an effective, risk-based environmental monitoring (EM) program is essential for finding and eliminating Lm from the processing environment to prevent product contamination. A scoping review was conducted to collate and synthesize available research and guidance materials on Listeria EM in food processing facilities. An exhaustive search was performed to identify all available research, industry and regulatory documents, and search results were screened for relevance based on eligibility criteria. After screening, 198 references were subjected to an in-depth review and categorized according to objectives for conducting Listeria sampling in food processing facilities and food sector. Mapping of the literature revealed research and guidance gaps by food sector, as fresh produce was the focus in only 10 references, compared to 72 on meat, 52 on fish and seafood, and 50 on dairy. Review of reported practices and guidance highlighted key design elements of EM, including the number, location, timing and frequency of sampling, as well as methods of detection and confirmation, and record-keeping. While utilization of molecular subtyping methods is a trend that will continue to advance understanding of Listeria contamination risks, improved study design and reporting standards by researchers will be essential to assist the food industry optimize their EM design and decision-making. The comprehensive collection of documents identified and synthesized in this review aids continued efforts to minimize the risk of Lm contaminated foods.

Natural selection differences detected in key protein domains between non-pathogenic and pathogenic Feline Coronavirus phenotypes.

Feline Coronaviruses (FCoVs) commonly cause mild enteric infections in felines worldwide (termed Feline Enteric Coronavirus [FECV]), with around 12% developing into deadly Feline Infectious Peritonitis (FIP; Feline Infectious Peritonitis Virus [FIPV]). Genomic differences between FECV and FIPV have been reported, yet the putative genotypic basis of the highly pathogenic phenotype remains unclear. Here, we used state-of-the-art molecular evolutionary genetic statistical techniques to identify and compare differences in natural selection pressure between FECV and FIPV sequences, as well as to identify FIPV and FECV specific signals of positive selection. We analyzed full length FCoV protein coding genes thought to contain mutations associated with FIPV (Spike, ORF3abc, and ORF7ab). We identified two sites exhibiting differences in natural selection pressure between FECV and FIPV: one within the S1/S2 furin cleavage site, and the other within the fusion domain of Spike. We also found 15 sites subject to positive selection associated with FIPV within Spike, 11 of which have not previously been suggested as possibly relevant to FIP development. These sites fall within Spike protein subdomains that participate in host cell receptor interaction, immune evasion, tropism shifts, host cellular entry, and viral escape. There were 14 sites (12 novel) within Spike under positive selection associated with the FECV phenotype, almost exclusively within the S1/S2 furin cleavage site and adjacent C domain, along with a signal of relaxed selection in FIPV relative to FECV, suggesting that furin cleavage functionality may not be needed for FIPV. Positive selection inferred in ORF7b was associated with the FECV phenotype, and included 24 positively selected sites, while ORF7b had signals of relaxed selection in FIPV. We found evidence of positive selection in ORF3c in FCoV wide analyses, but no specific association with the FIPV or FECV phenotype. We hypothesize that some combination of mutations in FECV may contribute to FIP development, and that is unlikely to be one singular "switch" mutational event. This work expands our understanding of the complexities of FIP development and provides insights into how evolutionary forces may alter pathogenesis in coronavirus genomes.

Genomes of endangered great hammerhead and shortfin mako sharks reveal historic population declines and high levels of inbreeding in great hammerhead.

Despite increasing threats of extinction to Elasmobranchii (sharks and rays), whole genome-based conservation insights are lacking. Here, we present chromosome-level genome assemblies for the Critically Endangered great hammerhead (Sphyrna mokarran) and the Endangered shortfin mako (Isurus oxyrinchus) sharks, with genetic diversity and historical demographic comparisons to other shark species. The great hammerhead exhibited low genetic variation, with 8.7% of the 2.77 Gbp genome in runs of homozygosity (ROH) > 1 Mbp and 74.4% in ROH >100 kbp. The 4.98 Gbp shortfin mako genome had considerably greater diversity and <1% in ROH > 1 Mbp. Both these sharks experienced precipitous declines in effective population size (Ne) over the last 250 thousand years. While shortfin mako exhibited a large historical Ne that may have enabled the retention of higher genetic variation, the genomic data suggest a possibly more concerning picture for the great hammerhead, and a need for evaluation with additional individuals.

Natural selection differences detected in key protein domains between non-pathogenic and pathogenic feline coronavirus phenotypes.

Feline coronaviruses (FCoVs) commonly cause mild enteric infections in felines worldwide (termed feline enteric coronavirus [FECV]), with around 12 per cent developing into deadly feline infectious peritonitis (FIP; feline infectious peritonitis virus [FIPV]). Genomic differences between FECV and FIPV have been reported, yet the putative genotypic basis of the highly pathogenic phenotype remains unclear. Here, we used state-of-the-art molecular evolutionary genetic statistical techniques to identify and compare differences in natural selection pressure between FECV and FIPV sequences, as well as to identify FIPV- and FECV-specific signals of positive selection. We analyzed full-length FCoV protein coding genes thought to contain mutations associated with FIPV (Spike, ORF3abc, and ORF7ab). We identified two sites exhibiting differences in natural selection pressure between FECV and FIPV: one within the S1/S2 furin cleavage site (FCS) and the other within the fusion domain of Spike. We also found fifteen sites subject to positive selection associated with FIPV within Spike, eleven of which have not previously been suggested as possibly relevant to FIP development. These sites fall within Spike protein subdomains that participate in host cell receptor interaction, immune evasion, tropism shifts, host cellular entry, and viral escape. There were fourteen sites (twelve novel sites) within Spike under positive selection associated with the FECV phenotype, almost exclusively within the S1/S2 FCS and adjacent to C domain, along with a signal of relaxed selection in FIPV relative to FECV, suggesting that furin cleavage functionality may not be needed for FIPV. Positive selection inferred in ORF7b was associated with the FECV phenotype and included twenty-four positively selected sites, while ORF7b had signals of relaxed selection in FIPV. We found evidence of positive selection in ORF3c in FCoV-wide analyses, but no specific association with the FIPV or FECV phenotype. We hypothesize that some combination of mutations in FECV may contribute to FIP development, and that it is unlikely to be one singular 'switch' mutational event. This work expands our understanding of the complexities of FIP development and provides insights into how evolutionary forces may alter pathogenesis in coronavirus genomes.

Recent Zoonotic Spillover and Tropism Shift of a Canine Coronavirus Is Associated with Relaxed Selection and Putative Loss of Function in NTD Subdomain of Spike Protein.

A canine coronavirus (CCoV) has now been reported from two independent human samples from Malaysia (respiratory, collected in 2017-2018; CCoV-HuPn-2018) and Haiti (urine, collected in 2017); these two viruses were nearly genetically identical. In an effort to identify any novel adaptations associated with this apparent shift in tropism we carried out detailed evolutionary analyses of the spike gene of this virus in the context of related Alphacoronavirus 1 species. The spike 0-domain retains homology to CCoV2b (enteric infections) and Transmissible Gastroenteritis Virus (TGEV; enteric and respiratory). This domain is subject to relaxed selection pressure and an increased rate of molecular evolution. It contains unique amino acid substitutions, including within a region important for sialic acid binding and pathogenesis in TGEV. Overall, the spike gene is extensively recombinant, with a feline coronavirus type II strain serving a prominent role in the recombinant history of the virus. Molecular divergence time for a segment of the gene where temporal signal could be determined, was estimated at around 60 years ago. We hypothesize that the virus had an enteric origin, but that it may be losing that particular tropism, possibly because of mutations in the sialic acid binding region of the spike 0-domain.

Characterizing infectious disease progression through discrete states using hidden Markov models.

Infectious disease management relies on accurate characterization of disease progression so that transmission can be prevented. Slowly progressing infectious diseases can be difficult to characterize because of a latency period between the time an individual is infected and when they show clinical signs of disease. The introduction of Mycobacterium avium ssp. paratuberculosis (MAP), the cause of Johne's disease, onto a dairy farm could be undetected by farmers for years before any animal shows clinical signs of disease. In this time period infected animals may shed thousands of colony forming units. Parameterizing trajectories through disease states from infection to clinical disease can help farmers to develop control programs based on targeting individual disease state, potentially reducing both transmission and production losses due to disease. We suspect that there are two distinct progression pathways; one where animals progress to a high-shedding disease state, and another where animals maintain a low-level of shedding without clinical disease. We fit continuous-time hidden Markov models to multi-year longitudinal fecal sampling data from three US dairy farms, and estimated model parameters using a modified Baum-Welch expectation maximization algorithm. Using posterior decoding, we observed two distinct shedding patterns: cows that had observations associated with a high-shedding disease state, and cows that did not. This model framework can be employed prospectively to determine which cows are likely to progress to clinical disease and may be applied to characterize disease progression of other slowly progressing infectious diseases.

New York State dairy farmers' perceptions of antibiotic use and resistance: A qualitative interview study.

Antibiotic resistance is a global problem affecting both human and animal health. Ensuring the strategic and effective use of antibiotics is paramount to combatting the emergence and spread of resistance. This study explored New York State (NYS) dairy farmers' perceptions regarding antibiotic use in dairy farming and antibiotic resistance. Dairy farmers' perceptions were assessed through semi-structured, in-person interviews. Twenty interviews with farm owners and/or managers of 15 conventional and five USDA certified organic dairy farms with 40 to 2,300 lactating cows were conducted. Thematic analysis was used to assess, compare and contrast transcripts for farmers' characterization of their beliefs, values, and concerns. Conventional dairy farmers had a low level of concern about the possible impacts of on-farm antibiotic resistance on human health and believed their antibiotic use was already judicious. Generally, they believed their cattle's health would suffer if antibiotic use were further curtailed. Conventional farmers expressed frustration over the possibility of more stringent governmental, milk cooperative, buyer, or marketer requirements for antibiotic use and associated animal welfare in the future. They attributed expanding regulations in part to misinformed consumer preferences, that farmers felt were influenced by the marketing of organic dairy products. Organic dairy farmers were generally more concerned about issues related to antibiotic resistance than conventional farmers. Both conventional and organic farmers placed emphasis on disease prevention through herd health management rather than treatment. In conclusion, the conventional NYS dairy farmers in this study were skeptical of the need for and benefits of reduced antibiotic use on their dairy farms. Interventions for farmers, delivered by a trusted source such as a veterinarian, that provide training about proper antibiotic use practices and information of possible financial benefits of refining antibiotic use may hold promise.

Immunohistochemical Labeling of Low-Density Lipoprotein Receptor and Scavenger Receptor Class B Type 1 Are Increased in Canine Lymphoma.

Altered lipid metabolism is a well-documented hallmark of neoplastic transformation and impacts disease progression. Two major lipoprotein receptors, the low-density lipoprotein receptor (LDL-R) and scavenger receptor class B, type 1 (SR-BI) are overexpressed in a number of cancer types in people. These receptors serve to deliver cholesterol to the tumor cells and have been used to target drug therapies. In this study, we performed a retrospective analysis of LDL-R and SR-B1 expression in canine lymphoma using archived formalin-fixed tissue samples. Cases were immunophenotyped and classified according to World Health Organization (WHO) standards prior to immunostaining for the LDL_R and SR-B1. A total of 45 cases were evaluated; 21 high grade B (HGB), 11 low grade B (LGB), 7 high grade T (HGT), and 6 low grade T (LGT) lymphomas. One sided Wilcoxon rank sum tests were used to compare staining intensity between neoplastic and hyperplastic lymphoid tissue. The relationships between histological score and tumor grade and score and stage at presentation were assessed using non-parametric Kruskal-Wallis tests. Neoplastic lymphoid tissue expressed higher levels of both receptors compared to reactive lymph nodes. Median LDL-R score was 85.0 (interquartile range = 101.7), Median SR-B1 score was 209.0 (interquartile range 105.2). No relationship between LDL-R or SR-B1 staining score and tumor grade or phenotype was found. Serum cholesterol concentration was compared between dogs with high and low grade tumors using a two sample T-test, and correlations between cholesterol concentration and histological score, and between the score for the two receptors were determined using a Spearman correlation. The high expression level of these lipoprotein receptors on most of the tumors could underlie the lack of relationship between score and tumor grade. The overexpression of LDL-R and SR-B1 in canine lymphoma holds therapeutic potential particularly in dogs that overexpress one or both of these receptors, and this warrants further investigation.