vSNP: a SNP pipeline for the generation of transparent SNP matrices and phylogenetic trees from whole genome sequencing data sets.

BACKGROUND: Several single nucleotide polymorphism (SNP) pipelines exist, each offering its own advantages. Among them and described here is vSNP that has been developed over the past decade and is specifically tailored to meet the needs of diagnostic laboratories. Laboratories that aim to provide rapid whole genome sequencing results during outbreak investigations face unique challenges. vSNP addresses these challenges by enabling users to verify and validate sequence accuracy with ease- having utility across various pathogens, being fully auditable, and presenting results that are easy to interpret and can be comprehended by individuals with diverse backgrounds. RESULTS: vSNP has proven effective for real-time phylogenetic analysis of disease outbreaks and eradication efforts, including bovine tuberculosis, brucellosis, virulent Newcastle disease, SARS-CoV-2, African swine fever, and highly pathogenic avian influenza. The pipeline produces easy-to-read SNP matrices, sorted for convenience, as well as corresponding phylogenetic trees, making the output easily understandable. Essential data for verifying SNPs is included in the output, and the process has been divided into two steps for ease of use and faster processing times. vSNP requires minimal computational resources to run and can be run in a wide range of environments. Several utilities have been developed to make analysis more accessible for subject matter experts who may not have computational expertise. CONCLUSION: The vSNP pipeline integrates seamlessly into a diagnostic workflow and meets the criteria for quality control accreditation programs, such as 17025 by the International Organization for Standardization. Its versatility and robustness make it suitable for use with a diverse range of organisms, providing detailed, reproducible, and transparent results, making it a valuable tool in various applications, including phylogenetic analysis performed in real time.

Concurrent colonization of rodent kidneys with multiple species and serogroups of pathogenic Leptospira.

Rodents are important reservoir hosts of pathogenic leptospires in the US Virgin Islands. Our previous work determined that trapped rodents were colonized with Leptospira borgpetersenii serogroup Ballum (n = 48) and/or Leptospira kirschneri serogroup Icterohaemorrhagiae (n = 3). In addition, nine rodents appeared to be colonized with a mixed population comprising more than one species/serogroup. The aim of this study was to validate this finding by characterizing clonal isolates derived from cultures of mixed species. Cultures of presumptive mixed species (designated LR1, LR5, LR37, LR57, LR60, LR61, LR68, LR70, and LR72) were propagated in different media including Hornsby-Alt-Nally (HAN) media, incubated at both 29℃ and 37℃, and T80/40/LH incubated at 29℃. Polyclonal reference antisera specific for serogroup Ballum and Icterohaemorrhagiae were used to enrich for different serogroups followed by subculture on agar plates. Individual colonies were then selected for genotyping and serotyping. Of the nine cultures of mixed species/serogroups, a single clonal isolate was separated in five of them: L. borgpetersenii serogroup Ballum in LR1, LR5, and LR37, and L. kirschneri serogroup Icterohaemorrhagiae in LR60 and LR72. In four of the cultures with mixed species (LR57, LR61, LR68, and LR70), clonal isolates of both L. borgpetersenii serogroup Ballum and L. kirschneri serogroup Icterohaemorrhagiae were recovered. Our results definitively establish that rodents can be colonized with more than one species/serogroup of Leptospira concurrently. The identification and characterization of multiple species/serogroups of Leptospira from individual reservoir hosts of infection are essential to understand the epidemiology and transmission of disease to both human and domestic animal populations.IMPORTANCEPathogenic Leptospira, the causative agent of human and animal leptospirosis, comprise a diverse genus of species/serogroups which are inherently difficult to isolate from mammalian hosts due to fastidious growth requirements. Molecular evidence has indicated that reservoir hosts of Leptospira may shed multiple species concurrently. However, evidence of this phenomena by culture has been lacking. Culture is definitive and is essential for comprehensive characterization of recovered isolates by high-resolution genome sequencing and serotyping. In this work, a protocol using recently developed novel media formulations, in conjunction with reference antisera, was developed and validated to demonstrate the recovery of multiple species/serogroups of pathogenic Leptospira from the same host. The identification and characterization of multiple species/serogroups of Leptospira from individual reservoir hosts of infection are essential to understand the epidemiology and transmission of disease to both human and domestic animal populations.

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.

Disease management at the wildlife-livestock interface: Using whole-genome sequencing to study the role of elk in Mycobacterium bovis transmission in Michigan, USA.

The role of wildlife in the persistence and spread of livestock diseases is difficult to quantify and control. These difficulties are exacerbated when several wildlife species are potentially involved. Bovine tuberculosis (bTB), caused by Mycobacterium bovis, has experienced an ecological shift in Michigan, with spillover from cattle leading to an endemically infected white-tailed deer (deer) population. It has potentially substantial implications for the health and well-being of both wildlife and livestock and incurs a significant economic cost to industry and government. Deer are known to act as a reservoir of infection, with evidence of M. bovis transmission to sympatric elk and cattle populations. However, the role of elk in the circulation of M. bovis is uncertain; they are few in number, but range further than deer, so may enable long distance spread. Combining Whole Genome Sequences (WGS) for M. bovis isolates from exceptionally well-observed populations of elk, deer and cattle with spatiotemporal locations, we use spatial and Bayesian phylogenetic analyses to show strong spatiotemporal admixture of M. bovis isolates. Clustering of bTB in elk and cattle suggests either intraspecies transmission within the two populations, or exposure to a common source. However, there is no support for significant pathogen transfer amongst elk and cattle, and our data are in accordance with existing evidence that interspecies transmission in Michigan is likely only maintained by deer. This study demonstrates the value of whole genome population studies of M. bovis transmission at the wildlife-livestock interface, providing insights into bTB management in an endemic system.

Identification of Source of Brucella suis Infection in Human by Using Whole-Genome Sequencing, United States and Tonga.

Brucella suis infection was diagnosed in a man from Tonga, Polynesia, who had butchered swine in Oregon, USA. Although the US commercial swine herd is designated brucellosis-free, exposure history suggested infection from commercial pigs. We used whole-genome sequencing to determine that the man was infected in Tonga, averting a field investigation.

Identification of equine mares as reservoir hosts for pathogenic species of Leptospira.

Equine leptospirosis can result in abortion, stillbirth, neonatal death, placentitis, and uveitis. Horses can also act as subclinical reservoir hosts of infection, which are characterized as asymptomatic carriers that persistently excrete leptospires and transmit disease. In this study, PCR and culture were used to assess urinary shedding of pathogenic Leptospira from 37 asymptomatic mares. Three asymptomatic mares, designated as H2, H8, and H9, were PCR-positive for lipL32, a gene specific for pathogenic species of Leptospira. One asymptomatic mare, H9, was culture-positive, and the recovered isolate was classified as L. kirschneri serogroup Australis serovar Rushan. DNA capture and enrichment of Leptospira genomic DNA from PCR-positive, culture-negative samples determined that asymptomatic mare H8 was also shedding L. kirschneri serogroup Australis, whereas asymptomatic mare H2 was shedding L. interrogans serogroup Icterohaemorrhagiae. Sera from all asymptomatic mares were tested by the microscopic agglutination test (MAT) and 35 of 37 (94.6%) were seropositive with titers ranging from 1:100 to 1:3200. In contrast to asymptomatic mares, mare H44 presented with acute spontaneous abortion and a serum MAT titer of 1:102,400 to L. interrogans serogroup Pomona serovar Pomona. Comparison of L. kirschneri serogroup Australis strain H9 with that of L. interrogans serogroup Pomona strain H44 in the hamster model of leptospirosis corroborated differences in virulence of strains. Since lipopolysaccharide (LPS) is a protective antigen in bacterin vaccines, the LPS of strain H9 (associated with subclinical carriage) was compared with strain H44 (associated with spontaneous abortion). This revealed different LPS profiles and immunoreactivity with reference antisera. It is essential to know what species and serovars of Leptospira are circulating in equine populations to design efficacious vaccines and diagnostic tests. Our results demonstrate that horses in the US can act as reservoir hosts of leptospirosis and shed diverse pathogenic Leptospira species via urine. This report also details the detection of L. kirschneri serogroup Australis serovar Rushan, a species and serotype of Leptospira, not previously reported in the US.

An accurate and interpretable model for antimicrobial resistance in pathogenic Escherichia coli from livestock and companion animal species.

Understanding the microbial genomic contributors to antimicrobial resistance (AMR) is essential for early detection of emerging AMR infections, a pressing global health threat in human and veterinary medicine. Here we used whole genome sequencing and antibiotic susceptibility test data from 980 disease causing Escherichia coli isolated from companion and farm animals to model AMR genotypes and phenotypes for 24 antibiotics. We determined the strength of genotype-to-phenotype relationships for 197 AMR genes with elastic net logistic regression. Model predictors were designed to evaluate different potential modes of AMR genotype translation into resistance phenotypes. Our results show a model that considers the presence of individual AMR genes and total number of AMR genes present from a set of genes known to confer resistance was able to accurately predict isolate resistance on average (mean F1 score = 98.0%, SD = 2.3%, mean accuracy = 98.2%, SD = 2.7%). However, fitted models sometimes varied for antibiotics in the same class and for the same antibiotic across animal hosts, suggesting heterogeneity in the genetic determinants of AMR resistance. We conclude that an interpretable AMR prediction model can be used to accurately predict resistance phenotypes across multiple host species and reveal testable hypotheses about how the mechanism of resistance may vary across antibiotics within the same class and across animal hosts for the same antibiotic.

Advancing animal tuberculosis surveillance using culture-independent long-read whole-genome sequencing.

Animal tuberculosis is a significant infectious disease affecting both livestock and wildlife populations worldwide. Effective disease surveillance and characterization of Mycobacterium bovis (M. bovis) strains are essential for understanding transmission dynamics and implementing control measures. Currently, sequencing of genomic information has relied on culture-based methods, which are time-consuming, resource-demanding, and concerning in terms of biosafety. This study explores the use of culture-independent long-read whole-genome sequencing (WGS) for a better understanding of M. bovis epidemiology in African buffaloes (Syncerus caffer). By comparing two sequencing approaches, we evaluated the efficacy of Illumina WGS performed on culture extracts and culture-independent Oxford Nanopore adaptive sampling (NAS). Our objective was to assess the potential of NAS to detect genomic variants without sample culture. In addition, culture-independent amplicon sequencing, targeting mycobacterial-specific housekeeping and full-length 16S rRNA genes, was applied to investigate the presence of microorganisms, including nontuberculous mycobacteria. The sequencing quality obtained from DNA extracted directly from tissues using NAS is comparable to the sequencing quality of reads generated from culture-derived DNA using both NAS and Illumina technologies. We present a new approach that provides complete and accurate genome sequence reconstruction, culture independently, and using an economically affordable technique.

Whole genome sequencing improves the discrimination between Mycobacterium bovis strains on the southern border of Kruger National Park, South Africa.

Background: Mycobacterium bovis forms part of the Mycobacterium tuberculosis complex and has an extensive host range and zoonotic potential. Various genotyping methods (e.g., spoligotyping) have been used to describe the molecular epidemiology of M. bovis. Advances in whole genome sequencing (WGS) have increased resolution to enable detection of genomic variants to the level of single nucleotide polymorphisms. This is especially relevant to One Health research on tuberculosis which benefits by being able to use WGS to identify epidemiologically linked cases, especially recent transmission. The use of WGS in molecular epidemiology has been extensively used in humans and cattle but is limited in wildlife. This approach appears to overcome the limitations of conventional genotyping methods due to lack of genetic diversity in M. bovis. Methods: This pilot study investigated the spoligotype and WGS of M. bovis isolates (n = 7) from wildlife in Marloth Park (MP) and compared these with WGS data from other South African M. bovis isolates. In addition, the greater resolution of WGS was used to explore the phylogenetic relatedness of M. bovis isolates in neighbouring wildlife populations. Results: The phylogenetic analyses showed the closest relatives to the seven isolates from MP were isolates from wildlife in Kruger National Park (KNP), which shares a border with MP. However, WGS data indicated that the KNP and MP isolates formed two distinct clades, even though they had similar spoligotypes and identical in silico genetic regions of difference profiles. Conclusions: Mycobacterium bovis isolates from MP were hypothesized to be directly linked to KNP wildlife, based on spoligotyping. However, WGS indicated more complex epidemiology. The presence of two distinct clades which were genetically distinct (SNP distance of 19-47) and suggested multiple transmission events. Therefore, WGS provided new insight into the molecular epidemiology of the M. bovis isolates from MP and their relationship to isolates from KNP. This approach will facilitate greater understanding of M. bovis transmission at wildlife-livestock-human interfaces and advances One Health research on tuberculosis, especially across different host species.

Molecular epidemiology of Brucella abortus isolates from cattle, elk, and bison in the United States, 1998 to 2011.

A variable-number tandem repeat (VNTR) protocol targeting 10 loci in the Brucella abortus genome was used to assess genetic diversity among 366 field isolates recovered from cattle, bison, and elk in the Greater Yellowstone Area (GYA) and Texas during 1998 to 2011. Minimum spanning tree (MST) and unweighted-pair group method with arithmetic mean (UPGMA) analyses of VNTR data identified 237 different VNTR types, among which 14 prominent clusters of isolates could be identified. Cattle isolates from Texas segregated into three clusters: one comprised of field isolates from 1998 to 2005, one comprised of vaccination-associated infections, and one associated with an outbreak in Starr County in January 2011. An isolate obtained from a feral sow trapped on property adjacent to the Starr County herd in May 2011 clustered with the cattle isolates, suggesting a role for feral swine as B. abortus reservoirs in Starr County. Isolates from a 2005 cattle outbreak in Wyoming displayed VNTR-10 profiles matching those of strains recovered from Wyoming and Idaho elk. Additionally, isolates associated with cattle outbreaks in Idaho in 2002, Montana in 2008 and 2011, and Wyoming in 2010 primarily clustered with isolates recovered from GYA elk. This study indicates that elk play a predominant role in the transmission of B. abortus to cattle located in the GYA.

Isolation of Leptospira kirschneri serovar Grippotyphosa from a red panda (Ailurus fulgens) after antimicrobial therapy: Case report.

A 1-year-old female red panda started showing symptoms of illness, including lethargy, anorexia, abdominal discomfort, and vomiting, shortly after transfer to a new zoo. Serum was tested for leptospirosis using the microscopic agglutination test, and a titer of 1:25,600 to serogroup Grippotyphosa was detected. Antimicrobial treatment with doxycycline was initiated. After completion of treatment and resolution of clinical symptoms, a urine sample was collected to ensure clearance of leptospires and cessation of urinary shedding prior to co-housing with other red pandas. A repeat serum sample taken 13 days later had a lower titer of 1:6,400 to serogroup Grippotyphosa. A sample of the animal's urine was cultured in HAN media and was culture positive for Leptospira. The recovered isolate was completely characterized by whole genome sequencing and serotyping with reference antisera, and the isolate was classified as Leptospira kirschneri serogroup Grippotyphosa serovar Grippotyphosa strain RedPanda1.

Use of Whole-Genome Sequencing to Unravel the Genetic Diversity of a Prevalent Mycobacterium bovis Spoligotype in a Multi-Host Scenario in Spain.

Despite the efforts invested in the eradication of bovine tuberculosis in Spain, herd prevalence has remained constant in the country during the last 15 years (~1.5-1.9%) due to a combination of epidemiological factors impairing disease control, including between-species transmission. Here, our aim was to investigate the molecular diversity of Mycobacterium bovis isolates belonging to the highly prevalent SB0339 spoligotype in the cattle-wildlife interface in different regions of Spain using whole-genome sequencing (WGS). Genomic data of 136 M. bovis isolates recovered from different animal species (cattle, wild boar, fallow deer, and red deer) and locations between 2005 and 2018 were analyzed to investigate between- and within-species transmission, as well as within-herds. All sequenced isolates differed by 49-88 single nucleotide polymorphisms from their most recent common ancestor. Genetic heterogeneity was geographic rather than host species-specific, as isolates recovered from both cattle and wildlife from a given region were more closely related compared to isolates from the same species but geographically distant. In fact, a strong association between the geographic and the genetic distances separating pairs of M. bovis isolates was found, with a significantly stronger effect when cattle isolates were compared with wildlife or cattle-wildlife isolates in Spain. The same results were obtained in Madrid, the region with the largest number of sequenced isolates, but no differences depending on the host were observed. Within-herd genetic diversity was limited despite the considerable time elapsed between isolations. The detection of closely related strains in different hosts demonstrates the complex between-host transmission dynamics present in endemic areas in Spain. In conclusion, WGS results a valuable tool to track bTB infection at a high resolution and may contribute to achieve its eradication in Spain.

Leptospira borgpetersenii serovar Hardjo and Leptospira santarosai serogroup Pyrogenes isolated from bovine dairy herds in Puerto Rico.

Leptospirosis is one of the most common zoonotic diseases in the world and endemic in the Caribbean Islands. Bovine leptospirosis is an important reproductive disease. Globally, cattle are recognized as a reservoir host for L. borgpetersenii serovar Hardjo, which is transmitted via urine, semen, and uterine discharges, and can result in abortion and poor reproductive performance. The dairy industry in Puerto Rico comprises up to 25% of agriculture-related income and is historically the most financially important agricultural commodity on the island. In this study, we report the isolation of two different pathogenic Leptospira species, from two different serogroups, from urine samples collected from dairy cows in Puerto Rico: L. borgpetersenii serogroup Sejroe serovar Hardjo and L. santarosai serogroup Pyrogenes. Recovered isolates were classified using whole-genome sequencing, serotyping with reference antisera and monoclonal antibodies, and immunoblotting. These results demonstrate that dairy herds in Puerto Rico can be concurrently infected with more than one species and serovar of Leptospira, and that bacterin vaccines and serologic diagnostics should account for this when applying intervention and diagnostic strategies.

Assessing rodents as carriers of pathogenic Leptospira species in the U.S. Virgin Islands and their risk to animal and public health.

Leptospirosis is a global zoonotic disease caused by pathogenic bacteria of the genus Leptospira. We sought to determine if rodents in U.S. Virgin Islands (USVI) are carriers of Leptospira. In total, 140 rodents were sampled, including 112 Mus musculus and 28 Rattus rattus. A positive carrier status was identified for 64/140 (45.7%); 49 (35.0%) were positive by dark-field microscopy, 60 (42.9%) by culture, 63 (45.0%) by fluorescent antibody testing, and 61 (43.6%) by real-time polymerase chain reaction (rtPCR). Molecular typing indicated that 48 isolates were L. borgpetersenii and 3 were L. kirschneri; the remaining nine comprised mixed species. In the single culture-negative sample that was rtPCR positive, genotyping directly from the kidney identified L. interrogans. Serotyping of L. borgpetersenii isolates identified serogroup Ballum and L. kirschneri isolates as serogroup Icterohaemorrhagiae. These results demonstrate that rodents are significant Leptospira carriers and adds to understanding the ecoepidemiology of leptospirosis in USVI.

Detection of SARS-CoV-2 Omicron variant (B.1.1.529) infection of white-tailed deer.

White-tailed deer ( Odocoileus virginianus ) are highly susceptible to infection by SARS-CoV-2, with multiple reports of widespread spillover of virus from humans to free-living deer. While the recently emerged SARS-CoV-2 B.1.1.529 Omicron variant of concern (VoC) has been shown to be notably more transmissible amongst humans, its ability to cause infection and spillover to non-human animals remains a challenge of concern. We found that 19 of the 131 (14.5%; 95% CI: 0.10-0.22) white-tailed deer opportunistically sampled on Staten Island, New York, between December 12, 2021, and January 31, 2022, were positive for SARS-CoV-2 specific serum antibodies using a surrogate virus neutralization assay, indicating prior exposure. The results also revealed strong evidence of age-dependence in antibody prevalence. A significantly (χ 2 , p < 0.001) greater proportion of yearling deer possessed neutralizing antibodies as compared with fawns (OR=12.7; 95% CI 4-37.5). Importantly, SARS-CoV-2 nucleic acid was detected in nasal swabs from seven of 68 (10.29%; 95% CI: 0.0-0.20) of the sampled deer, and whole-genome sequencing identified the SARS-CoV-2 Omicron VoC (B.1.1.529) is circulating amongst the white-tailed deer on Staten Island. Phylogenetic analyses revealed the deer Omicron sequences clustered closely with other, recently reported Omicron sequences recovered from infected humans in New York City and elsewhere, consistent with human to deer spillover. Interestingly, one individual deer was positive for viral RNA and had a high level of neutralizing antibodies, suggesting either rapid serological conversion during an ongoing infection or a "breakthrough" infection in a previously exposed animal. Together, our findings show that the SARS-CoV-2 B.1.1.529 Omicron VoC can infect white-tailed deer and highlights an urgent need for comprehensive surveillance of susceptible animal species to identify ecological transmission networks and better assess the potential risks of spillback to humans. KEY FINDINGS: These studies provide strong evidence of infection of free-living white-tailed deer with the SARS-CoV-2 B.1.1.529 Omicron variant of concern on Staten Island, New York, and highlight an urgent need for investigations on human-to-animal-to-human spillovers/spillbacks as well as on better defining the expanding host-range of SARS-CoV-2 in non-human animals and the environment.

SARS-CoV-2 Omicron (B.1.1.529) Infection of Wild White-Tailed Deer in New York City.

There is mounting evidence of SARS-CoV-2 spillover from humans into many domestic, companion, and wild animal species. Research indicates that humans have infected white-tailed deer, and that deer-to-deer transmission has occurred, indicating that deer could be a wildlife reservoir and a source of novel SARS-CoV-2 variants. We examined the hypothesis that the Omicron variant is actively and asymptomatically infecting the free-ranging deer of New York City. Between December 2021 and February 2022, 155 deer on Staten Island, New York, were anesthetized and examined for gross abnormalities and illnesses. Paired nasopharyngeal swabs and blood samples were collected and analyzed for the presence of SARS-CoV-2 RNA and antibodies. Of 135 serum samples, 19 (14.1%) indicated SARS-CoV-2 exposure, and 11 reacted most strongly to the wild-type B.1 lineage. Of the 71 swabs, 8 were positive for SARS-CoV-2 RNA (4 Omicron and 4 Delta). Two of the animals had active infections and robust neutralizing antibodies, revealing evidence of reinfection or early seroconversion in deer. Variants of concern continue to circulate among and may reinfect US deer populations, and establish enzootic transmission cycles in the wild: this warrants a coordinated One Health response, to proactively surveil, identify, and curtail variants of concern before they can spill back into humans.

Whole-Genome SNP Analysis Identifies Putative Mycobacterium bovis Transmission Clusters in Livestock and Wildlife in Catalonia, Spain.

The high-resolution WGS analyses of MTBC strains have provided useful insight for determining sources of infection for animal tuberculosis. In Spain, tuberculosis in livestock is caused by Mycobacterium bovis and Mycobacterium caprae, where wildlife reservoirs play an important role. We analyzed a set of 125 M. bovis isolates obtained from livestock and wildlife from Catalonia to investigate strain diversity and identify possible sources and/or causes of infection. Whole-genome SNP profiles were used for phylogenetic reconstruction and pairwise SNP distance analysis. Additionally, SNPs were investigated to identify virulence and antimicrobial resistance factors to investigate clade-specific associations. Putative transmission clusters (≤12 SNPs) were identified, and associated epidemiological metadata were used to determine possible explanatory factors for transmission. M. bovis distribution was heterogeneous, with 7 major clades and 21 putative transmission clusters. In order of importance, the explanatory factors associated were proximity and neighborhood, residual infection, livestock-wildlife interaction, shared pasture, and movement. Genes related to lipid transport and metabolism showed the highest number of SNPs. All isolates were pyrazinamide resistant, and five were additionally resistant to isoniazid, but no clade-specific associations could be determined. Our findings highlight the importance of high-resolution molecular surveillance to monitor bovine tuberculosis dynamics in a low-prevalence setting.

Who infects whom?-Reconstructing infection chains of Mycobacterium avium ssp. paratuberculosis in an endemically infected dairy herd by use of genomic data.

Recent evidence of circulation of multiple strains within herds and mixed infections of cows marks the beginning of a rethink of our knowledge on Mycobacterium avium ssp. paratuberculosis (MAP) epidemiology. Strain typing opens new ways to investigate MAP transmission. This work presents a method for reconstructing infection chains in a setting of endemic Johne's disease on a well-managed dairy farm. By linking genomic data with demographic field data, strain-specific differences in spreading patterns could be quantified for a densely sampled dairy herd. Mixed infections of dairy cows with MAP are common, and some strains spread more successfully. Infected cows remain susceptible for co-infections with other MAP genotypes. The model suggested that cows acquired infection from 1-4 other cows and spread infection to 0-17 individuals. Reconstructed infection chains supported the hypothesis that high shedding animals that started to shed at an early age and showed a progressive infection pattern represented a greater risk for spreading MAP. Transmission of more than one genotype between animals was recorded. In this farm with a good MAP control management program, adult-to-adult contact was proposed as the most important transmission route to explain the reconstructed networks. For each isolate, at least one more likely ancestor could be inferred. Our study results help to capture underlying transmission processes and to understand the challenges of tracing MAP spread within a herd. Only the combination of precise longitudinal field data and bacterial strain type information made it possible to trace infection in such detail.

Exposure and Carriage of Pathogenic Leptospira in Livestock in St. Croix, U.S. Virgin Islands.

From 2019-2020, the Virgin Islands Department of Health (VIDOH) investigated potential animal reservoirs of Leptospira spp., the pathogenic bacteria that cause leptospirosis. We examined Leptospira exposure and carriage in livestock on the island of St. Croix, United States Virgin Islands (USVI). We utilized the microscopic agglutination test (MAT) to evaluate the sera, and the fluorescent antibody test (FAT), real time polymerase chain reaction (rt-PCR), and bacterial culture to evaluate urine specimens from livestock (n = 126): 28 cattle, 19 goats, 46 pigs, and 33 sheep. Seropositivity was 37.6% (47/125) with agglutinating antibodies to the following serogroups identified: Australis, Djasiman, Icterohaemorrhagiae, Ballum, Sejroe, Cynopteri, Autumnalis, Hebdomadis, Pomona, Canicola, Grippotyphosa, and Pyrogenes. Urine from 4 animals (4.0%, 4/101) was positive by rt-PCR for lipL32: 2 sheep, 1 goat, and 1 bull. Sequencing of secY amplicons identified L. interrogans in 1 sheep and 1 bull. Livestock in USVI harbor pathogenic Leptospira bacteria and could play a role in the zoonotic cycle of leptospirosis.

Whole Genome Sequencing Links Mycobacterium bovis From Cattle, Cheese and Humans in Baja California, Mexico.

Mycobacterium bovis causes tuberculosis (TB) in cattle, which in turn can transmit the pathogen to humans. Tuberculosis in dairy cattle is of particular concern where the consumption of raw milk and dairy products is customary. Baja California (BCA), Mexico, presents high prevalence of TB in both cattle and humans, making it important to investigate the molecular epidemiology of the disease in the region. A long-term study was undertaken to fully characterize the diversity of M. bovis genotypes circulating in dairy cattle, cheese and humans in BCA by whole-genome sequencing (WGS). During a 2-year period, 412 granulomatous tissue samples were collected from local abattoirs and 314 cheese samples were purchased from local stores and vendors in BCA and sent to the laboratory for mycobacterial culture, histology, direct PCR and WGS. For tissue samples M. bovis was recovered from 86.8%, direct PCR detected 90% and histology confirmed 85.9% as mycobacteriosis-compatible. For cheese, M. bovis was recovered from 2.5% and direct PCR detected 6% of the samples. There was good agreement between diagnostic tests. Subsequently, a total of 345 whole-genome SNP sequences were obtained. Phylogenetic analysis grouped these isolates into 10 major clades. SNP analysis revealed putative transmission clusters where the pairwise SNP distance between isolates from different dairies was ≤3 SNP. Also, human and/or cheese isolates were within 8.45 (range 0-17) and 5.8 SNP (range 0-15), respectively, from cattle isolates. Finally, a comparison between the genotypes obtained in this study and those reported previously suggests that the genetic diversity of M. bovis in BCA is well-characterized, and can be used to determine if BCA is the likely source of M. bovis in humans and cattle in routine epidemiologic investigations and future studies. In conclusion, WGS provided evidence of ongoing local transmission of M. bovis among the dairies in this high-TB burden region of BCA, as well as show close relationships between isolates recovered from humans, cheese, and cattle. This confirms the need for a coordinated One Health approach in addressing the elimination of TB in animals and humans. Overall, the study contributes to the knowledge of the molecular epidemiology of M. bovis in BCA, providing insight into the pathogen's dynamics in a high prevalence setting.