Search

It is not possible to systematically screen the environment for rabies virus (RABV) using current approaches. We sought to determine under what conditions RABV is detectable from feces and other accessible samples from infected wildlife to broaden the number of biological samples that could be used to test for RABV. We employed a recently-developed quantitative RT-PCR assay called the "LN34 panlyssavirus real-time RT-PCR assay", which is highly sensitive and specific for all variants of RABV. We harvested and tested brain tissue, fecal, and/or mouth swab samples from 25 confirmed RABV positive bats of six species. To determine if rabies RNA lasts in feces sufficiently long post-defecation to use it as a surveillance tool, we tested fecal samples from 10 bats at the time of sample collection and after 24 hours of exposure to ambient conditions, with an additional test on six bats out to 72 hours. To assess whether we could pool fecal pellets and still detect a positive, we generated dilutions of known positives at 1:1, 1:10, 1:50, and 1:200. For six individuals for which matched brain, mouth swab, and fecal samples were tested, results were positive for 100%, 67%, and 67%, respectively. For the first time test to 24 hours, 63% of feces that were positive at time 0 were still positive after 24 hours, and 50% of samples at 72 hours were positive across all three replicates. Pooling tests revealed that fecal positives were detected at 1:10 dilution, but not at 1:50 or 1:200. Our preliminary results suggest that fecal samples hold promise for a rapid and non-invasive environmental screening system. [ABSTRACT FROM AUTHOR]

Background: Plague, caused by the bacterium Yersinia pestis, remains an important disease in Madagascar, where the oriental rat flea, Xenopsylla cheopis, is a primary vector. To control fleas, synthetic pyrethroids (SPs) have been used for >20 years, resulting in resistance in many X. cheopis populations. The most common mechanisms of SP resistance are target site mutations in the voltage-gated sodium channel (VGSC) gene. Methodology/Principal findings: We obtained 25 collections of X. cheopis from 22 locations across Madagascar and performed phenotypic tests to determine resistance to deltamethrin, permethrin, and/or dichlorodiphenyltrichloroethane (DDT). Most populations were resistant to all these insecticides. We sequenced a 535 bp segment of the VGSC gene and identified two different mutations encoding distinct substitutions at amino acid position 1014, which is associated with knockdown resistance (kdr) to SPs in insects. Kdr mutation L1014F occurred in all 25 collections; a rarer mutation, L1014H, was found in 12 collections. There was a significant positive relationship between the frequency of kdr alleles and the proportion of individuals surviving exposure to deltamethrin. Phylogenetic comparisons of 12 VGSC alleles in Madagascar suggested resistant alleles arose from susceptible lineages at least three times. Because genotype can reasonably predict resistance phenotype, we developed a TaqMan PCR assay for the rapid detection of kdr resistance alleles. Conclusions/Significance: Our study provides new insights into VGSC mutations in Malagasy populations of X. cheopis and is the first to report a positive correlation between VGSC genotypes and SP resistance phenotypes in fleas. Widespread occurrence of these two SP resistance mutations in X. cheopis populations in Madagascar reduces the viability of these insecticides for flea control. However, the TaqMan assay described here facilitates rapid detection of kdr mutations to inform when use of these insecticides is still warranted to reduce transmission of plague. Author summary: In Madagascar, the Oriental rat flea Xenopsylla cheopis is the primary vector of Yersinia pestis, which causes plague. This flea has been reported to be resistant to insecticides used for vector control. The objective of this study was to examine the target site mutations in the voltage-gated sodium channel (VGSC) gene in X. cheopis populations sampled throughout the island and previously tested for synthetic pyrethroid (SPs) susceptibility. Two widespread mutations associated with knock-down resistance (kdr) to SPs were identified, and a significant relationship between the frequency of kdr alleles and the SP resistance phenotype was revealed. The results raise concerns about the overall efficacy of SPs for controlling Yersinia pestis transmission by X. cheopis in Madagascar. [ABSTRACT FROM AUTHOR]

Onchocerca lupi is a filarial nematode that causes ocular onchocercosis in canines globally including North America and areas of Europe, North Africa, and the Middle East. Reported incidence of this parasite in canines has continued to steadily escalate since the early 21st century and was more recently documented in humans. Whole genome sequencing (WGS) of this parasite can provide insight into gene content, provide novel surveillance targets, and elucidate the origin and range expansion. However, past attempts of whole genome sequencing of other Onchocerca species reported a substantial portion of their data unusable due to the variable over-abundance of host DNA in samples. Here, we have developed a method to determine the host-to-parasite DNA ratio using a quantitative PCR (qPCR) approach that relies on two standard plasmids each of which contains a single copy gene specific to the parasite genus Onchocerca (major body wall myosin gene, myosin) or a single copy gene specific to the canine host (polycystin-1 precursor, pkd1). These plasmid standards were used to determine the copy number of the myosin and pkd1 genes within a sample to calculate the ratio of parasite and host DNA. Furthermore, whole genome sequence (WGS) data for three O. lupi isolates were consistent with our host-to-parasite DNA ratio results. Our study demonstrates, despite unified DNA extraction methods, variable quantities of host DNA within any one sample which will likely affect downstream WGS applications. Our quantification assay of host-to-parasite genome copy number provides a robust and accurate method of assessing canine host DNA load in an O. lupi specimen that will allow informed sample selection for WGS. This study has also provided the first whole genome draft sequence for this species. This approach is also useful for future focused WGS studies of other parasites. [ABSTRACT FROM AUTHOR]

Genomic surveillance and wastewater tracking strategies were used to strengthen the public health response to an outbreak of the SARS-CoV-2 Delta AY.25 lineage associated with a university campus in Arizona. Epidemiologic and clinical data routinely gathered through contact tracing were matched to SARS-CoV-2 genomes belonging to an outbreak of AY.25 identified through ongoing phylogenomic analyses. Continued phylogenetic analyses were conducted to further describe the AY.25 outbreak. Wastewater collected twice weekly from sites across campus was tested for SARS-CoV-2 by RT-qPCR, and subsequently sequenced to identify variants. The AY.25 outbreak was defined by a single mutation (C18804T) and comprised 379 genomes from SARS-CoV-2 positive cases associated with the university and community. Several undergraduate student gatherings and congregate living settings on campus likely contributed to the rapid spread of COVID-19 across the university with secondary transmission into the community. The clade defining mutation was also found in wastewater samples collected from around student dormitories a week before the semester began, and 9 days before cases were identified. Genomic, epidemiologic, and wastewater surveillance provided evidence that an AY.25 clone was likely imported into the university setting just prior to the onset of the Fall 2021 semester, rapidly spread through a subset of the student population, and then subsequent spillover occurred in the surrounding community. The university and local public health department worked closely together to facilitate timely reporting of cases, identification of close contacts, and other necessary response and mitigation strategies. The emergence of new SARS-CoV-2 variants and potential threat of other infectious disease outbreaks on university campuses presents an opportunity for future comprehensive One Health genomic data driven, targeted interventions. [ABSTRACT FROM AUTHOR]

In this paper, we describe a population of mothers who are opioid dependent at the time of giving birth and neonates exposed to opioids in utero who experience withdrawal following birth. While there have been studies of national trends in this population, there remains a gap in studies of regional trends. Using data from the Arizona Department of Health Services Hospital Discharge Database, this study aimed to characterize the population of neonates with neonatal opioid withdrawal syndrome (NOWS) and mothers who were opioid dependent at the time of giving birth, in Arizona. We analyzed approximately 1.2 million electronic medical records from the Arizona Department of Health Services Hospital Discharge Database to identify patterns and disparities across socioeconomic, ethnic, racial, and/or geographic groupings. In addition, we identified comorbid conditions that are differentially associated with NOWS in neonates or opioid dependence in mothers. Our analysis was designed to assess whether indicators such as race/ethnicity, insurance payer, marital status, and comorbidities are related to the use of opioids while pregnant. Our findings suggest that women and neonates who are non-Hispanic White and economically disadvantaged, tend be part of our populations of interest more frequently than expected. Additionally, women who are opioid dependent at the time of giving birth are unmarried more often than expected, and we suggest that marital status could be a proxy for support. Finally, we identified comorbidities associated with neonates who have NOWS and mothers who are opioid dependent not previously reported. [ABSTRACT FROM AUTHOR]

Due to the large number of negative tests, individually screening large populations for rare pathogens can be wasteful and expensive. Sample pooling methods improve the efficiency of large-scale pathogen screening campaigns by reducing the number of tests and reagents required to accurately categorize positive and negative individuals. Such methods rely on group testing theory which mainly focuses on minimizing the total number of tests; however, many other practical concerns and tradeoffs must be considered when choosing an appropriate method for a given set of circumstances. Here we use computational simulations to determine how several theoretical approaches compare in terms of (a) the number of tests, to minimize costs and save reagents, (b) the number of sequential steps, to reduce the time it takes to complete the assay, (c) the number of samples per pool, to avoid the limits of detection, (d) simplicity, to reduce the risk of human error, and (e) robustness, to poor estimates of the number of positive samples. We found that established methods often perform very well in one area but very poorly in others. Therefore, we introduce and validate a new method which performs fairly well across each of the above criteria making it a good general use approach. [ABSTRACT FROM AUTHOR]

Eastern tarsiers (Tarsius tarsier complex) are small nocturnal primates endemic to Sulawesi Island and small adjacent islands of Indonesia. In 2004, the hybrid biogeography hypothesis predicted this species complex might contain 16 or more taxa, each corresponding to a region of endemism, based on: 1) geological evidence of the development of the archipelago, 2) biological evidence in the form of concordant distributions of monkeys and toads, and 3) the distribution of tarsier acoustic groups. Since then, 11 tarsier species have been recognized, potentially leaving more to be described. Efforts to identify these cryptic species are urgently needed so that habitat conversion, pet trade, and cultural activities will not render some species extinct before they are recognized. We gathered data to test the hypothesis of cryptic tarsier species on three volcanic islands in Bunaken National Park, North Sulawesi, namely Bunaken, Manadotua, and Mantehage, during May-August 2018. We sequenced individuals at 5 nuclear genes (ABCA1, ADORA3, AXIN1, RAG, and TTR) and made comparisons to existing genotypes at 14 mainland sites. Bayesian phylogenetic analyses revealed that island populations are genetically identical in all 5 genes, and formed a clade separated from the mainland ones. The eastern tarsiers first diverged from the western tarsiers approximately 2.5 MYA. The three island populations diverged from mainland tarsiers approximately 2,000–150,000 YA, due to either human activities or natural rafting. This study provides information for tarsier conservation, advances the understanding of biogeography of Sulawesi, and contributes to Indonesian awareness of biodiversity. Further quantitative genetics research on tarsiers, especially the island populations, will offer significant insights to establish more efficient and strategic tarsier conservation actions. [ABSTRACT FROM AUTHOR]

Although acute melioidosis is the most common outcome of Burkholderia pseudomallei infection, we have documented a case, P314, where disease severity lessened with time, and the pathogen evolved towards a commensal relationship with the host. In the current study, we used whole-genome sequencing to monitor this long-term symbiotic relationship to better understand B. pseudomallei persistence in P314's sputum despite intensive initial therapeutic regimens. We collected and sequenced 118 B. pseudomallei isolates from P314's airways over a >16-year period, and also sampled the patient's home environment, recovering six closely related B. pseudomallei isolates from the household water system. Using comparative genomics, we identified 126 SNPs in the core genome of the 124 isolates or 162 SNPs/indels when the accessory genome was included. The core SNPs were used to construct a phylogenetic tree, which demonstrated a close relationship between environmental and clinical isolates and detailed within-host evolutionary patterns. The phylogeny had little homoplasy, consistent with a strictly clonal mode of genetic inheritance. Repeated sampling revealed evidence of genetic diversification, but frequent extinctions left only one successful lineage through the first four years and two lineages after that. Overall, the evolution of this population is nonadaptive and best explained by genetic drift. However, some genetic and phenotypic changes are consistent with in situ adaptation. Using a mouse model, P314 isolates caused greatly reduced morbidity and mortality compared to the environmental isolates. Additionally, potentially adaptive phenotypes emerged and included differences in the O-antigen, capsular polysaccharide, motility, and colony morphology. The >13-year co-existence of two long-lived lineages presents interesting hypotheses that can be tested in future studies to provide additional insights into selective pressures, niche differentiation, and microbial adaptation. This unusual melioidosis case presents a rare example of the evolutionary progression towards commensalism by a highly virulent pathogen within a single human host. Author summary: Pathogens frequently jump between different hosts, and associated adaptation may lead to the emergence of new infectious agents. Such host-jumping evolution is witnessed through endpoint analyses but these cannot capture genetic changes in lineages that have gone extinct. In this study, we have identified and monitored an example of the adaptive evolution of a bacterium often deadly to its mammalian host, in an unprecedented case whereby disease lessened through time and the pathogen adapted to become a part of the commensal human flora. We used genomic analyses to characterize more than 16 years of this evolutionary process and the stepwise mutations that control pathogen interactions with the patient. Soon after infection, mutational changes occurred that allowed the bacterium to remain in the airways without causing disease. This shift towards avirulence was determined based on clinical data and virulence testing in an animal model. In addition, mutations occurred that contributed to the persistence of the bacteria in the patient's lungs. Finally, we found evidence for the evolutionary emergence and persistence of two distinct lineages of the bacterium over the last 13 years, presenting interesting questions about niche utilization. Bacteria are ubiquitous in the human body and almost all are beneficial or benign. In this study, we document the evolutionary conversion of a normally deadly bacterium towards a commensal. [ABSTRACT FROM AUTHOR]

West Nile Virus (WNV) has been detected annually in Maricopa County, Arizona, since 2003. With this in mind, we sought to determine if contemporary strains are endemic to the county or are annually imported. As part of this effort, we developed a new protocol for tiled amplicon sequencing of WNV to efficiently attain greater than 99% coverage of 14 WNV genomes collected directly from positive mosquito pools distributed throughout Maricopa County between 2014 and 2017. Bayesian phylogenetic analyses revealed that contemporary genomes fall within two major lineages; NA/WN02 and SW/WN03. We found that all of the Arizona strains possessed an amino acid substitution known to be under positive selection, which has arisen independently at least four times in Arizona. The SW/WN03 strains exhibited transient behavior, with at least 10 separate introductions into Arizona when considering both historical and contemporary strains. However, NA/WN02 strains are geographically differentiated and appear to be endemic in Arizona, with two clades that have been circulating for four and seven years. This establishment in Maricopa County provides the first evidence of local overwintering by a WNV strain over the course of several years in Arizona. Within a national context, the placement of eleven contemporary Arizona strains in the NA/WN02 lineage indicates while WNV first entered the northeastern United States in 1999, the most ancestral extant strains of WNV are now circulating in the American southwest. [ABSTRACT FROM AUTHOR]

Burkholderia pseudomallei is the environmental bacillus that causes melioidosis; a disease clinically significant in Australia and Southeast Asia but emerging in tropical and sub-tropical regions around the globe. Previous studies have placed the ancestral population of the organism in Australia with a single lineage disseminated to Southeast Asia. We have previously characterized B. pseudomallei isolates from New Guinea and the Torres Strait archipelago; remote regions that share paleogeographic ties with Australia. These studies identified regional biogeographical boundaries. In this study, we utilize whole-genome sequencing to reconstruct ancient evolutionary relationships and ascertain correlations between paleogeography and present-day distribution of this bacterium in Australasia. Our results indicate that B. pseudomallei from New Guinea fall into a single clade within the Australian population. Furthermore, clades from New Guinea are region-specific; an observation possibly linked to limited recent anthropogenic influence in comparison to mainland Australia and Southeast Asia. Isolates from the Torres Strait archipelago were distinct yet scattered among those from mainland Australia. These results provide evidence that the New Guinean and Torres Strait lineages may be remnants of an ancient portion of the Australian population. Rising sea levels isolated New Guinea and the Torres Strait Islands from each other and the Australian mainland, and may have allowed long-term isolated evolution of these lineages, providing support for a theory of microbial biogeography congruent with that of macro flora and fauna. Moreover, these findings indicate that contemporary microbial biogeography theories should consider recent and ongoing impacts of globalisation and human activity. [ABSTRACT FROM AUTHOR]

Background: Yersinia pestis appears to be maintained in multiple, geographically separate, and phylogenetically distinct subpopulations within the highlands of Madagascar. However, the dynamics of these locally differentiated subpopulations through time are mostly unknown. To address that gap and further inform our understanding of plague epidemiology, we investigated the phylogeography of Y. pestis in Madagascar over an 18 year period. Methodology/Principal findings: We generated whole genome sequences for 31 strains and discovered new SNPs that we used in conjunction with previously identified SNPs and variable-number tandem repeats (VNTRs) to genotype 773 Malagasy Y. pestis samples from 1995 to 2012. We mapped the locations where samples were obtained on a fine geographic scale to examine phylogeographic patterns through time. We identified 18 geographically separate and phylogenetically distinct subpopulations that display spatial and temporal stability, persisting in the same locations over a period of almost two decades. We found that geographic areas with higher levels of topographical relief are associated with greater levels of phylogenetic diversity and that sampling frequency can vary considerably among subpopulations and from year to year. We also found evidence of various Y. pestis dispersal events, including over long distances, but no evidence that any dispersal events resulted in successful establishment of a transferred genotype in a new location during the examined time period. Conclusions/Significance: Our analysis suggests that persistent endemic cycles of Y. pestis transmission within local areas are responsible for the long term maintenance of plague in Madagascar, rather than repeated episodes of wide scale epidemic spread. Landscape likely plays a role in maintaining Y. pestis subpopulations in Madagascar, with increased topographical relief associated with increased levels of localized differentiation. Local ecological factors likely affect the dynamics of individual subpopulations and the associated likelihood of observing human plague cases in a given year in a particular location. [ABSTRACT FROM AUTHOR]

Background: Leptospirosis is a zoonotic disease responsible for high morbidity around the world, especially in tropical and low income countries. Rats are thought to be the main vector of human leptospirosis in urban settings. However, differences between urban and low-income rural communities provide additional insights into the epidemiology of the disease. Methodology/Principal findings: Our study was conducted in two low-income rural communities near the coast of Ecuador. We detected and characterized infectious leptospira DNA in a wide variety of samples using new real time quantitative PCR assays and amplicon sequencing. We detected infectious leptospira in a high percentage of febrile patients (14.7%). In contrast to previous studies on leptospirosis risk factors, higher positivity was not found in rats (3.0%) but rather in cows (35.8%) and pigs (21.1%). Six leptospira species were identified (L. borgpetersenii, L kirschnerii, L santarosai, L. interrogans, L noguchii, and an intermediate species within the L. licerasiae and L. wolffii clade) and no significant differences in the species of leptospira present in each animal species was detected (χ2 = 9.89, adj.p-value = 0.27). Conclusions/Significance: A large portion of the world’s human population lives in low-income, rural communities, however, there is limited information about leptospirosis transmission dynamics in these settings. In these areas, exposure to peridomestic livestock is particularly common and high prevalence of infectious leptospira in cows and pigs suggest that they may be the most important reservoir for human transmission. Genotyping clinical samples show that multiple species of leptospira are involved in human disease. As these genotypes were also detected in samples from a variety of animals, genotype data must be used in conjunction with epidemiological data to provide evidence of transmission and the importance of different potential leptospirosis reservoirs. [ABSTRACT FROM AUTHOR]

Staphylococcus aureus is an important clinical pathogen worldwide and understanding this organism's phylogeny and, in particular, the role of recombination, is important both to understand the overall spread of virulent lineages and to characterize outbreaks. To further elucidate the phylogeny of S. aureus, 35 diverse strains were sequenced using whole genome sequencing. In addition, 29 publicly available whole genome sequences were included to create a single nucleotide polymorphism (SNP)-based phylogenetic tree encompassing 11 distinct lineages. All strains of a particular sequence type fell into the same clade with clear groupings of the major clonal complexes of CC8, CC5, CC30, CC45 and CC1. Using a novel analysis method, we plotted the homoplasy density and SNP density across the whole genome and found evidence of recombination throughout the entire chromosome, but when we examined individual clonal lineages we found very little recombination. However, when we analyzed three branches of multiple lineages, we saw intermediate and differing levels of recombination between them. These data demonstrate that in S. aureus, recombination occurs across major lineages that subsequently expand in a clonal manner. Estimated mutation rates for the CC8 and CC5 lineages were different from each other. While the CC8 lineage rate was similar to previous studies, the CC5 lineage was 100-fold greater. Fifty known virulence genes were screened in all genomes in silico to determine their distribution across major clades. Thirty-three genes were present variably across clades, most of which were not constrained by ancestry, indicating horizontal gene transfer or gene loss. [ABSTRACT FROM AUTHOR]

Background: Glanders, caused by the gram-negative bacterium Burkholderia mallei, is a highly infectious zoonotic disease of solipeds causing severe disease in animals and men. Although eradicated from many Western countries, it recently emerged in Asia, the Middle-East, Africa, and South America. Due to its rareness, little is known about outbreak dynamics of the disease and its epidemiology. Methodology/Principal Findings: We investigated a recent outbreak of glanders in Bahrain by applying high resolution genotyping (multiple locus variable number of tandem repeats, MLVA) and comparative whole genome sequencing to B. mallei isolated from infected horses and a camel. These results were compared to samples obtained from an outbreak in the United Arab Emirates in 2004, and further placed into a broader phylogeographic context based on previously published B. mallei data. The samples from the outbreak in Bahrain separated into two distinct clusters, suggesting a complex epidemiological background and evidence for the involvement of multiple B. mallei strains. Additionally, the samples from Bahrain were more closely related to B. mallei isolated from horses in the United Arab Emirates in 2004 than other B. mallei which is suggestive of repeated importation to the region from similar geographic sources. Conclusion/Significance: High-resolution genotyping and comparative whole genome analysis revealed the same phylogenetic patterns among our samples. The close relationship of the Dubai/UAE B. mallei populations to each other may be indicative of a similar geographic origin that has yet to be identified for the infecting strains. The recent emergence of glanders in combination with worldwide horse trading might pose a new risk for human infections. Author Summary: Glanders is a disease of antiquity, recognized as a malady of equines by Hippocrates and Aristotle. The causative agent, Burkholderia mallei, is currently feared as a potential biological weapon and has been used as such in the American Civil War and both World Wars to cripple equine military components. In the more economically developed countries, glanders has been eradicated through large scale culling. As a result, our understanding of transmission dynamics and networks is limited. However, regions of endemicity still exist in Asia, the Middle-East, Africa, and South America where it infects solipeds and camels. These areas provide reservoirs for re-introduction of glanders into countries previously listed as glanders-free. Here, we demonstrate the utility of high-resolution genotyping and whole genome sequence analysis in the investigation of a recent outbreak of glanders in horses and camels in Bahrain, a previously declared glanders-free country. Our analyses demonstrate that not one, but two strains likely caused this outbreak, and that these strains probably came from a similar geographic region via importation of infected animals. Even with careful monitoring, the global trade of animals from glanders-endemic regions can re-introduce and possibly re-establish this disease in animal populations of countries that have previously eradicated it. [ABSTRACT FROM AUTHOR]

The implementation of non-pharmaceutical public health interventions can have simultaneous impacts on pathogen transmission rates as well as host mobility rates. For instance, with SARS-CoV-2, masking can influence host-to-host transmission, while stay-at-home orders can influence mobility. Importantly, variations in transmission rates and mobility patterns can influence pathogen-induced hospitalization rates. This poses a significant challenge for the use of mathematical models of disease dynamics in forecasting the spread of a pathogen; to create accurate forecasts in spatial models of disease spread, we must simultaneously account for time-varying rates of transmission and host movement. In this study, we develop a statistical model-fitting algorithm to estimate dynamic rates of SARS-CoV-2 transmission and host movement from geo-referenced hospitalization data. Using simulated data sets, we then test whether our method can accurately estimate these time-varying rates simultaneously, and how this accuracy is influenced by the spatial population structure. Our model-fitting method relies on a highly parallelized process of grid search and a sliding window technique that allows us to estimate time-varying transmission rates with high accuracy and precision, as well as movement rates with somewhat lower precision. Estimated parameters also had lower precision in more rural data sets, due to lower hospitalization rates (i.e., these areas are less data-rich). This model-fitting routine could easily be generalized to any stochastic, spatially-explicit modeling framework, offering a flexible and efficient method to estimate time-varying parameters from geo-referenced data sets., Competing Interests: The authors have declared that no competing interests exist., (Copyright: © 2022 Ratnavale et al. This is an open access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.)