A mutation associated with resistance to synthetic pyrethroids is widespread in US populations of the tropical lineage of Rhipicephalus sanguineus s.l

The brown dog tick, Rhipicephalus sanguineus sensu lato (s.l.), is an important vector for Rickettsia rickettsii, causative agent of Rocky Mountain spotted fever. Current public health prevention and control efforts to protect people involve preventing tick infestations on domestic animals and in and around houses. Primary prevention tools rely on acaricides, often synthetic pyrethroids (SPs); resistance to this chemical class is widespread in ticks and other arthropods. Rhipicephalus sanguineus s.l. is a complex that likely contains multiple unique species and although the distribution of this complex is global, there are differences in morphology, ecology, and perhaps vector competence among these major lineages. Two major lineages within Rh. sanguineus s.l., commonly referred to as temperate and tropical, have been documented from multiple locations in North America, but are thought to occupy different ecological niches. To evaluate potential acaricide resistance and better define the distributions of the tropical and temperate lineages throughout the US and in northern Mexico, we employed a highly multiplexed amplicon sequencing approach to characterize sequence diversity at: 1) three loci within the voltage-gated sodium channel (VGSC) gene, which contains numerous genetic mutations associated with resistance to SPs; 2) a region of the gamma-aminobutyric acid-gated chloride channel gene (GABA-Cl) containing several mutations associated with dieldrin/fipronil resistance in other species; and 3) three mitochondrial genes (COI, 12S, and 16S). We utilized a geographically diverse set of Rh sanguineus s.l. collected from domestic pets in the US in 2013 and a smaller set of ticks collected from canines in Baja California, Mexico in 2021. We determined that a single nucleotide polymorphism (T2134C) in domain III segment 6 of the VGSC, which has previously been associated with SP resistance in Rh. sanguineus s.l., was widespread and abundant in tropical lineage ticks (>50 %) but absent from the temperate lineage, suggesting that resistance to SPs may be common in the tropical lineage. We found evidence of multiple copies of GABA-Cl in ticks from both lineages, with some copies containing mutations associated with fipronil resistance in other species, but the effects of these patterns on fipronil resistance in Rh. sanguineus s.l. are currently unknown. The tropical lineage was abundant and geographically widespread, accounting for 79 % of analyzed ticks and present at 13/14 collection sites. The temperate and tropical lineages co-occurred in four US states, and as far north as New York. None of the ticks we examined were positive for Rickettsia rickettsii or Rickettsia massiliae.

Knockdown resistance mutations are common and widely distributed in Xenopsylla cheopis fleas that transmit plague in Madagascar

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.

Quantification of mechanical properties in long-term in vivo used silicone catheter lines according to DIN 10555-3

BACKGROUND

Totally implantable central venous access port devices are crucial for intravenous application of chemotherapeutics and long-term therapy for chronic disease. Common complications include thrombosis and device fracture related to altered material properties through exposure in situ. This study exhibits whether uniaxial tensile properties (DIN 10555-3) of in vivo used catheters prove inferior to unused catheters.

MATERIAL AND METHODS

5 unused, originally packed silicone catheters were cut into 6 segments of 50 mm: 3 segments each were cleaned via cleaning solution (n = 15) while 3 segments were left unattended (n = 15). Distal segments (50 mm) of long-term in vivo used silicone catheters were cleaned before testing (n = 33). Overall mechanical behavior was tested in a custom-made self-centering, torsion free carrier. Maximum force stress at failure, strain at failure and Young's modulus were determined and statistically analyzed.

FINDINGS

Unused catheters showed no significant difference in testing. in vivo used catheters exhibited 20% lower maximal force than unused catheters (p < 0.001), strain at break (p 〈0,001), and 7% higher elastic modulus (p = 0.004; power: 0.845). Due to a constant cross section area, stress at failure was proportional to maximum force (p < 0.001). Relation between determined parameters and dwell times was non-significant.

INTERPRETATION

In vivo long-term used silicone catheters showed significantly lower ultimate strength than unused ones. It is likely that in situ altering changes the mechanical properties of catheters and may potentially lead to failure.

DNA Capture and Enrichment: A Culture-Independent Approach for Characterizing the Genomic Diversity of Pathogenic Leptospira Species

Because they are difficult to culture, obtaining genomic information from Leptospira spp. is challenging, hindering the overall understanding of leptospirosis. We designed and validated a culture-independent DNA capture and enrichment system for obtaining Leptospira genomic information from complex human and animal samples. It can be utilized with a variety of complex sample types and diverse species as it was designed using the pan-genome of all known pathogenic Leptospira spp. This system significantly increases the proportion of Leptospira DNA contained within DNA extracts obtained from complex samples, oftentimes reaching >95% even when some estimated starting proportions were <1%. Sequencing enriched extracts results in genomic coverage similar to sequenced isolates, thereby enabling enriched complex extracts to be analyzed together with whole genome sequences from isolates, which facilitates robust species identification and high-resolution genotyping. The system is flexible and can be readily updated when new genomic information becomes available. Implementation of this DNA capture and enrichment system will improve efforts to obtain genomic data from unculturable Leptospira-positive human and animal samples. This, in turn, will lead to a better understanding of the overall genomic diversity and gene content of Leptospira spp. that cause leptospirosis, aiding epidemiology and the development of improved diagnostics and vaccines.

Burkholderia thailandensis Isolated from the Environment, United States

Burkholderia thailandensis, an opportunistic pathogen found in the environment, is a bacterium closely related to B. pseudomallei, the cause of melioidosis. Human B. thailandensis infections are uncommon. We isolated B. thailandensis from water in Texas and Puerto Rico and soil in Mississippi in the United States, demonstrating a potential public health risk.

Low risk of acquiring melioidosis from the environment in the continental United States

Melioidosis is an underreported human disease of tropical and sub-tropical regions caused by the saprophyte Burkholderia pseudomallei. Although most global melioidosis cases are reported from tropical regions in Southeast Asia and northern Australia, there are multiple occurrences from sub-tropical regions, including the United States (U.S.). Most melioidosis cases reported from the continental U.S. are the result of acquiring the disease during travel to endemic regions or from contaminated imported materials. Only two human melioidosis cases from the continental U.S. have likely acquired B. pseudomallei directly from local environments and these cases lived only ~7 km from each other in rural Texas. In this study, we assessed the risk of acquiring melioidosis from the environment within the continental U.S. by surveying for B. pseudomallei in the environment in Texas where these two human melioidosis cases likely acquired their infections. We sampled the environment near the homes of the two cases and at additional sampling locations in surrounding counties in Texas that were selected based on ecological niche modeling. B. pseudomallei was not detected at the residences of these two cases or in the surrounding region. These negative data are important to demonstrate that B. pseudomallei is rare in the environment in the U.S. even at locations where locally acquired human cases likely have occurred, documenting the low risk of acquiring B. pseudomallei infection from the environment in the continental U.S.

Expanding the Burkholderia pseudomallei Complex with the Addition of Two Novel Species: Burkholderia mayonis sp. nov. and Burkholderia savannae sp. nov

Distinct Burkholderia strains were isolated from soil samples collected in tropical northern Australia (Northern Territory and the Torres Strait Islands, Queensland). Phylogenetic analysis of 16S rRNA and whole genome sequences revealed these strains were distinct from previously described Burkholderia species and assigned them to two novel clades within the B. pseudomallei complex (Bpc). Because average nucleotide identity and digital DNA-DNA hybridization calculations are consistent with these clades representing distinct species, we propose the names Burkholderia mayonis sp. nov. and Burkholderia savannae sp. nov. Strains assigned to B. mayonis sp. nov. include type strain BDU6T (=TSD-80; LMG 29941; ASM152374v2) and BDU8. Strains assigned to B. savannae sp. nov. include type strain MSMB266T (=TSD-82; LMG 29940; ASM152444v2), MSMB852, BDU18, and BDU19. Comparative genomics revealed unique coding regions for both putative species, including clusters of orthologous genes associated with phage. Type strains of both B. mayonis sp. nov. and B. savannae sp. nov. yielded biochemical profiles distinct from each other and from other species in the Bpc, and profiles also varied among strains within B. mayonis sp. nov. and B. savannae sp. nov. Matrix-assisted laser desorption ionization time-of-flight (MLST) analysis revealed a B. savannae sp. nov. cluster separate from other species, whereas B. mayonis sp. nov. strains did not form a distinct cluster. Neither B. mayonis sp. nov. nor B. savannae sp. nov. caused mortality in mice when delivered via the subcutaneous route. The addition of B. mayonis sp. nov. and B. savannae sp. nov. results in a total of eight species currently within the Bpc. IMPORTANCEBurkholderia species can be important sources of novel natural products, and new species are of interest to diverse scientific disciplines. Although many Burkholderia species are saprophytic, Burkholderia pseudomallei is the causative agent of the disease melioidosis. Understanding the genomics and virulence of the closest relatives to B. pseudomallei, i.e., the other species within the B. pseudomallei complex (Bpc), is important for identifying robust diagnostic targets specific to B. pseudomallei and for understanding the evolution of virulence in B. pseudomallei. Two proposed novel species, B. mayonis sp. nov. and B. savannae sp. nov., were isolated from soil samples collected from multiple locations in northern Australia. The two proposed species belong to the Bpc but are phylogenetically distinct from all other members of this complex. The addition of B. mayonis sp. nov. and B. savannae sp. nov. results in a total of eight species within this significant complex of bacteria that are available for future studies.

No evidence for enzootic plague within black-tailed prairie dog (Cynomys ludovicianus) populations

Yersinia pestis, causative agent of plague, occurs throughout the western United States in rodent populations and periodically causes epizootics in susceptible species, including black‐tailed prairie dogs (Cynomys ludovicianus). How Y. pestis persists long‐term in the environment between these epizootics is poorly understood but multiple mechanisms have been proposed, including, among others, a separate enzootic transmission cycle that maintains Y. pestis without involvement of epizootic hosts and persistence of Y. pestis within epizootic host populations without causing high mortality within those populations. We live‐trapped and collected fleas from black‐tailed prairie dogs and other mammal species from sites with and without black‐tailed prairie dogs in 2004 and 2005 and tested all fleas for presence of Y. pestis. Y. pestis was not detected in 2126 fleas collected in 2004 but was detected in 294 fleas collected from multiple sites in 2005, before and during a widespread epizootic that drastically reduced black‐tailed prairie dog populations in the affected colonies. Temporal and spatial patterns of Y. pestis occurrence in fleas and genotyping of Y. pestis present in some infected fleas suggest Y. pestis was introduced multiple times from sources outside the study area and once introduced, was dispersed between several sites. We conclude Y. pestis likely was not present in these black‐tailed prairie dog colonies prior to epizootic activity in these colonies. Although we did not identify likely enzootic hosts, we found evidence that deer mice (Peromyscus maniculatus) may serve as bridging hosts for Y. pestis between unknown enzootic hosts and black‐tailed prairie dogs.

Long-Term Survival of Virulent Tularemia Pathogens outside a Host in Conditions That Mimic Natural Aquatic Environments

Francisella tularensis, the causative agent of the zoonotic disease tularemia, can cause seasonal outbreaks of acute febrile illness in humans with disease peaks in late summer to autumn. Interestingly, its mechanisms for environmental persistence between outbreaks are poorly understood. One hypothesis is that F. tularensis forms biofilms in aquatic environments. We utilized two fully virulent wild-type strains: FSC200 (Francisella tularensis subsp. holarctica) and Schu S4 (Francisella tularensis subsp. tularensis) and three control strains, the attenuated live vaccine strain (LVS; F. tularensis subsp. holarctica), a Schu S4 ΔwbtI mutant that is documented to form biofilms, and the low-virulence strain U112 of the closely related species Francisella novicida Strains were incubated in saline solution (0.9% NaCl) microcosms for 24 weeks at both 4°C and 20°C, whereupon viability and biofilm formation were measured. These temperatures were selected to approximate winter and summer temperatures of fresh water in Scandinavia, respectively. U112 and Schu S4 ΔwbtI formed biofilms, but F. tularensis strains FSC200 and Schu S4 and the LVS did not. All strains exhibited prolonged viability at 4°C compared to 20°C. U112 and FSC200 displayed remarkable long-term persistence at 4°C, with only 1- and 2-fold log reductions, respectively, of viable cells after 24 weeks. Schu S4 exhibited lower survival, yielding no viable cells by week 20. At 24 weeks, cells from FSC200, but not from Schu S4, were still fully virulent in mice. Taken together, these results demonstrate biofilm-independent, long-term survival of pathogenic F. tularensis subsp. holarctica in conditions that mimic overwinter survival in aquatic environments.IMPORTANCE Tularemia, a disease caused by the environmental bacterium Francisella tularensis, is characterized by acute febrile illness. F. tularensis is highly infectious: as few as 10 organisms can cause human disease. Tularemia is not known to be spread from person to person. Rather, all human infections are independently acquired from the environment via the bite of blood-feeding arthropods, ingestion of infected food or water, or inhalation of aerosolized bacteria. Despite the environmental origins of human disease events, the ecological factors governing the long-term persistence of F. tularensis in nature between seasonal human outbreaks are poorly understood. The significance of our research is in identifying conditions that promote long-term survival of fully virulent F. tularensis outside a mammalian host or insect vector. These conditions are similar to those found in natural aquatic environments in winter and provide important new insights on how F. tularensis may persist long-term in the environment.

Burkholderia pseudomallei in Soil, US Virgin Islands, 2019

The distribution of Burkholderia pseudomallei in the Caribbean is poorly understood. We isolated B. pseudomallei from US Virgin Islands soil. The soil isolate was genetically similar to other isolates from the Caribbean, suggesting that B. pseudomallei might have been introduced to the islands multiple times through severe weather events.

COYOTES (CANIS LATRANS) IN ARIZONA, USA, EXHIBIT IMMUNE AND GENETIC EVIDENCE OF RICKETTSIAL INFECTIONS

Rocky Mountain spotted fever (RMSF), caused by the bacterium Rickettsia rickettsii, was recognized as endemic in Arizona, US after a 2002 outbreak and has since been a public health concern. The brown dog tick (Rhipicephalus sanguineus sensu lato) is the principal vector of this pathogen in Arizona. Domesticated dogs (Canis lupus familiaris) are the tick's main host, so free-roaming dogs in peridomestic areas have been named the primary risk factor for human cases of RMSF. However, the sudden emergence and long-distance dispersal of the pathogen have not been adequately explained, and one possible mechanism could include wildlife. Coyotes (Canis latrans) are wide ranging in Arizona and closely related to dogs, so it is possible that brown dog ticks parasitize coyotes and infect them. Although R. rickettsii is the most severe spotted fever group (SFG) rickettsial pathogen in humans, others occur in Arizona, and antibodies raised against them are cross-reactive, so we more-broadly hypothesized that coyotes in Arizona are exposed to SFG rickettsiae. We collected coyote tissues in spring 2016 and 2017. We tested sera for antibodies to R. rickettsii and found 9% (8/94) of samples were antibody-positive with titers of ≥256. Subsequent quantitative PCR analyses of skin showed evidence for Rickettsia spp. in 2.9% (4/138) of samples. These data suggest that coyotes have a role in the maintenance of SFG rickettsiae in Arizona. Further investigation is warranted to reveal which specific pathogen-vector complexes act on coyotes in the region and whether they represent a risk to human health.

Burkholderia pseudomallei, the causative agent of melioidosis, is rare but ecologically established and widely dispersed in the environment in Puerto Rico

BACKGROUND

Burkholderia pseudomallei is a soil-dwelling bacterium and the causative agent of melioidosis. The global burden and distribution of melioidosis is poorly understood, including in the Caribbean. B. pseudomallei was previously isolated from humans and soil in eastern Puerto Rico but the abundance and distribution of B. pseudomallei in Puerto Rico as a whole has not been thoroughly investigated.

METHODOLOGY/PRINCIPAL FINDINGS

We collected 600 environmental samples (500 soil and 100 water) from 60 sites around Puerto Rico. We identified B. pseudomallei by isolating it via culturing and/or using PCR to detect its DNA within complex DNA extracts. Only three adjacent soil samples from one site were positive for B. pseudomallei with PCR; we obtained 55 isolates from two of these samples. The 55 B. pseudomallei isolates exhibited fine-scale variation in the core genome and contained four novel genomic islands. Phylogenetic analyses grouped Puerto Rico B. pseudomallei isolates into a monophyletic clade containing other Caribbean isolates, which was nested inside a larger clade containing all isolates from Central/South America. Other Burkholderia species were commonly observed in Puerto Rico; we cultured 129 isolates from multiple soil and water samples collected at numerous sites around Puerto Rico, including representatives of B. anthina, B. cenocepacia, B. cepacia, B. contaminans, B. glumae, B. seminalis, B. stagnalis, B. ubonensis, and several unidentified novel Burkholderia spp.

CONCLUSIONS/SIGNIFICANCE

B. pseudomallei was only detected in three soil samples collected at one site in north central Puerto Rico with only two of those samples yielding isolates. All previous human and environmental B. pseudomallei isolates were obtained from eastern Puerto Rico. These findings suggest B. pseudomallei is ecologically established and widely dispersed in the environment in Puerto Rico but rare. Phylogeographic patterns suggest the source of B. pseudomallei populations in Puerto Rico and elsewhere in the Caribbean may have been Central or South America.

Climate-driven reduction of genetic variation in plant phenology alters soil communities and nutrient pools

We examined the hypothesis that climate-driven evolution of plant traits will influence associated soil microbiomes and ecosystem function across the landscape. Using a foundation tree species, Populus angustifolia, observational and common garden approaches, and a base population genetic collection that spans 17 river systems in the western United States, from AZ to MT, we show that (a) as mean annual temperature (MAT) increases, genetic and phenotypic variation for bud break phenology decline; (b) soil microbiomes, soil nitrogen (N), and soil carbon (C) vary in response to MAT and conditioning by trees; and (c) with losses of genetic variation due to warming, population-level regulation of community and ecosystem functions strengthen. These results demonstrate a relationship between the potential evolutionary response of populations and subsequent shifts in ecosystem function along a large temperature gradient.

Domestic canines do not display evidence of gut microbial dysbiosis in the presence of Clostridioides (Clostridium) difficile, despite cellular susceptibility to its toxins

Clostridioides difficile infection (CDI) is an emerging public health threat and C. difficile is the most common cause of antimicrobial-associated diarrhea worldwide and the leading cause of hospital-associated infections in the US, yet the burden of community-acquired infections (CAI) is poorly understood. Characterizing C. difficile isolated from canines is important for understanding the role that canines may play in CAI. In addition, several studies have suggested that canines carry toxigenic C. difficile asymptomatically, which may imply that there are mechanisms responsible for resistance to CDI in canines that could be exploited to help combat human CDI. To assess the virulence potential of canine-derived C. difficile, we tested whether toxins TcdA and TcdB (hereafter toxins) derived from a canine isolate were capable of causing tight junction disruptions to colonic epithelial cells. Additionally, we addressed whether major differences exist between human and canine cells regarding C. difficile pathogenicity by exposing them to identical toxins. We then examined the canine gut microbiome associated with C. difficile carriage using 16S rRNA gene sequencing and searched for deviations from homeostasis as an indicator of CDI. Finally, we queried 16S rRNA gene sequences for bacterial taxa that may be associated with resistance to CDI in canines. Clostridioides difficile isolated from a canine produced toxins that reduced tight junction integrity in both human and canine cells in vitro. However, canine guts were not dysbiotic in the presence of C. difficile. These findings support asymptomatic carriage in canines and, furthermore, suggest that there are features of the gut microbiome and/or a canine-specific immune response that may protect canines against CDI. We identified two biologically relevant bacteria that may aid in CDI resistance in canines: 1) Clostridium hiranonis, which synthesizes secondary bile acids that have been shown to provide resistance to CDI in mice; and 2) Sphingobacterium faecium, which produces sphingophospholipids that may be associated with regulating homeostasis in the canine gut. Our findings suggest that canines may be cryptic reservoirs for C. difficile and, furthermore, that mechanisms of CDI resistance in the canine gut could provide insights into targeted therapeutics for human CDI.

Caprine humoral response to Burkholderia pseudomallei antigens during acute melioidosis from aerosol exposure

Burkholderia pseudomallei causes melioidosis, a common source of pneumonia and sepsis in Southeast Asia and Northern Australia that results in high mortality rates. A caprine melioidosis model of aerosol infection that leads to a systemic infection has the potential to characterize the humoral immune response. This could help identify immunogenic proteins for new diagnostics and vaccine candidates. Outbred goats may more accurately mimic human infection, in contrast to the inbred mouse models used to date. B. pseudomallei infection was delivered as an intratracheal aerosol. Antigenic protein profiling was generated from the infecting strain MSHR511. Humoral immune responses were analyzed by ELISA and western blot, and the antigenic proteins were identified by mass spectrometry. Throughout the course of the infection the assay results demonstrated a much greater humoral response with IgG antibodies, in both breadth and quantity, compared to IgM antibodies. Pre-infection sera showed multiple immunogenic proteins already reactive for IgG (7-20) and IgM (0-12) in most of the goats despite no previous exposure to B. pseudomallei. After infection, the number of IgG reactive proteins showed a marked increase as the disease progressed. Early stage infection (day 7) showed immune reaction to chaperone proteins (GroEL, EF-Tu, and DnaK). These three proteins were detected in all serum samples after infection, with GroEL immunogenically dominant. Seven common reactive antigens were selected for further analysis using ELISA. The heat shock protein GroEL1 elicited the strongest goat antibody immune response compared to the other six antigens. Most of the six antigens showed the peak IgM reactivity at day 14, whereas the IgG reactivity increased further as the disease progressed. An overall MSHR511 proteomic comparison between the goat model and human sera showed that many immune reactive proteins are common between humans and goats with melioidosis.

Phylogenetic analysis of West Nile Virus in Maricopa County, Arizona: Evidence for dynamic behavior of strains in two major lineages in the American Southwest

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.

Range-wide genetic analysis of Dermacentor variabilis and its Francisella-like endosymbionts demonstrates phylogeographic concordance between both taxa

BACKGROUND

The American dog tick, Dermacentor variabilis, is an important vector of pathogens to humans, wildlife and domestic animals in North America. Although this tick species is widely distributed in the USA and Canada, knowledge of its range-wide phylogeographic patterns remains incomplete.

METHODS

We carried out a phylogenetic analysis of D. variabilis using samples collected from 26 USA states and five Canadian provinces. Tick samples (n = 1053 in total) originated from two main sources: existing archives (2000-2011), and new collections made from 2012 to 2013. We sequenced a 691 bp fragment of the cox1 gene from a subset (n = 332) of geographically diverse D. variabilis. DNA extracted from individual ticks (n = 1053) was also screened for a Francisella-like endosymbiont, using a targeted 16S rRNA sequencing approach, and important pathogens (Rickettsia spp. and Coxiella burnetii), using species-specific quantitative PCR assays.

RESULTS

Maximum parsimony analysis of cox1 sequences revealed two major groups within D. variabilis with distinct geographical distributions: one from the eastern USA/Canada (Group 1) and one from the west coast states of the USA (California and Washington; Group 2). However, genetic subdivisions within both of these two major groups were weak to moderate and not tightly correlated with geography. We found molecular signatures consistent with Francisella-like endosymbionts in 257 of the DNA extracts from the 1053 individual ticks, as well as Rickettsia spp. and Coxiella burnetii in a small number of ticks (n = 29 and 2, respectively). Phylogenetic patterns for Francisella-like endosymbionts, constructed using sequence data from the bacterial 16S rRNA locus, were similar to those for D. variabilis, with two major groups that had a nearly perfect one-to-one correlation with the two major groups within D. variabilis.

CONCLUSIONS

Our findings reveal a distinct phylogenetic split between the two major D. variabilis populations. However, high levels of genetic mixture among widely separated geographical localities occur within each of these two major groups. Furthermore, our phylogenetic analyses provide evidence of long-term tick-symbiont co-evolution. This work has implications for understanding the dispersal and evolutionary ecology of D. variabilis and associated vector-borne diseases.

Test method for mechanical properties of implantable catheters according to DIN 10555-3

To enable causal analysis of port catheter failure, this study aimed to develop an experimental setup for uniaxial tensile tests that addresses the specific requirements of highly elastic medical catheters; and to quantify parameters of the catheters' mechanical competence with respect to effects of artificial aging. Segments of 6F-polyurethane catheters were tested in their native status, after chemical and after mechanical aging. Tension experiments were performed with a rate of 220 mm/min until catheter failure. Material behavior was analyzed based on load cell measurements of the universal test system and an additional optical distance registration. The Young's modulus, the ultimate stress and the ultimate strain were determined. Chemical aging significantly decreased Young's modulus (84%; p = 0.001) and ultimate stress (83%; p < 0.001), whereas mechanical aged samples demonstrated similar results for the Young's modulus (p = 0.772) and a non-significant rise of ultimate stress (13%; p = 0.128). Ultimate strain did not differ significantly regardless of the pretreatment. The results proof reliability, reproducibility and sensitivity to quantify artificial aging induced variations and also promise to detect deviations in material features caused by long-term clinical usage of catheters.

Phylogeographic, genomic, and meropenem susceptibility analysis of Burkholderia ubonensis

The bacterium Burkholderia ubonensis is commonly co-isolated from environmental specimens harbouring the melioidosis pathogen, Burkholderia pseudomallei. B. ubonensis has been reported in northern Australia and Thailand but not North America, suggesting similar geographic distribution to B. pseudomallei. Unlike most other Burkholderia cepacia complex (Bcc) species, B. ubonensis is considered non-pathogenic, although its virulence potential has not been tested. Antibiotic resistance in B. ubonensis, particularly towards drugs used to treat the most severe B. pseudomallei infections, has also been poorly characterised. This study examined the population biology of B. ubonensis, and includes the first reported isolates from the Caribbean. Phylogenomic analysis of 264 B. ubonensis genomes identified distinct clades that corresponded with geographic origin, similar to B. pseudomallei. A small proportion (4%) of strains lacked the 920kb chromosome III replicon, with discordance of presence/absence amongst genetically highly related strains, demonstrating that the third chromosome of B. ubonensis, like other Bcc species, probably encodes for a nonessential pC3 megaplasmid. Multilocus sequence typing using the B. pseudomallei scheme revealed that one-third of strains lack the "housekeeping" narK locus. In comparison, all strains could be genotyped using the Bcc scheme. Several strains possessed high-level meropenem resistance (≥32 μg/mL), a concern due to potential transmission of this phenotype to B. pseudomallei. In silico analysis uncovered a high degree of heterogeneity among the lipopolysaccharide O-antigen cluster loci, with at least 35 different variants identified. Finally, we show that Asian B. ubonensis isolate RF23-BP41 is avirulent in the BALB/c mouse model via a subcutaneous route of infection. Our results provide several new insights into the biology of this understudied species.

Burkholderia humptydooensis sp. nov., a New Species Related to Burkholderia thailandensis and the Fifth Member of the Burkholderia pseudomallei Complex

During routine screening for Burkholderia pseudomallei from water wells in northern Australia in areas where it is endemic, Gram-negative bacteria (strains MSMB43T, MSMB121, and MSMB122) with a similar morphology and biochemical pattern to B. pseudomallei and B. thailandensis were coisolated with B. pseudomallei on Ashdown's selective agar. To determine the exact taxonomic position of these strains and to distinguish them from B. pseudomallei and B. thailandensis, they were subjected to a series of phenotypic and molecular analyses. Biochemical and fatty acid methyl ester analysis was unable to distinguish B. humptydooensis sp. nov. from closely related species. With matrix-assisted laser desorption ionization-time of flight analysis, all isolates grouped together in a cluster separate from other Burkholderia spp. 16S rRNA and recA sequence analyses demonstrated phylogenetic placement for B. humptydooensis sp. nov. in a novel clade within the B. pseudomallei group. Multilocus sequence typing (MLST) analysis of the three isolates in comparison with MLST data from 3,340 B. pseudomallei strains and related taxa revealed a new sequence type (ST318). Genome-to-genome distance calculations and the average nucleotide identity of all isolates to both B. thailandensis and B. pseudomallei, based on whole-genome sequences, also confirmed B. humptydooensis sp. nov. as a novel Burkholderia species within the B. pseudomallei complex. Molecular analyses clearly demonstrated that strains MSMB43T, MSMB121, and MSMB122 belong to a novel Burkholderia species for which the name Burkholderia humptydooensis sp. nov. is proposed, with the type strain MSMB43T (American Type Culture Collection BAA-2767; Belgian Co-ordinated Collections of Microorganisms LMG 29471; DDBJ accession numbers CP013380 to CP013382).IMPORTANCEBurkholderia pseudomallei is a soil-dwelling bacterium and the causative agent of melioidosis. The genus Burkholderia consists of a diverse group of species, with the closest relatives of B. pseudomallei referred to as the B. pseudomallei complex. A proposed novel species, B. humptydooensis sp. nov., was isolated from a bore water sample from the Northern Territory in Australia. B. humptydooensis sp. nov. is phylogenetically distinct from B. pseudomallei and other members of the B. pseudomallei complex, making it the fifth member of this important group of bacteria.

Radiologically inserted gastrostomy: differences of maintenance of balloon- vs. loop-retained devices

PURPOSE

To compare outcome and associated complications of ballon- vs. loop-retained devices for radiologically inserted gastrostomy (RIG).

METHODS

From 2007 to 2011 233 patients (age 63.7 ± 10.6 years) were referred for a RIG because of pharyngeal stricture Intervention was performed with four different devices: balloon-retained - Freka^® GastroTube, Fresenius Kabi (n = 121); MIC^® Gastrostomy Feeding Tube, Kimberly-Clark (n = 34); Russell^® Gastrostomy Tray, Cook Medical Inc. (n = 17); and loop-retained - Tilma^® Gastrostomy Set, Cook Medical Inc. (n = 50). Follow-up was performed with regard to RIG-related complications, cause of removal and fatalities. Revision-free survival times after RIG were evaluated using Kaplan-Meier analysis and group differences by log-rank tests. For analysis of demographic and methodical variables multivariate Cox regression models were used.

RESULTS

With a primary technical success rate of 95.3% (222/233) a total of 92 instances of revisions were necessary in 66 patients (66/233, 28.3%) during follow-up (mean 182.8 ± 86.6 days). The most common complication was tube dislodgement (14.3%). There were no significant differences between the distinct devices (p = 0.098), but analyzing the data in subgroups of balloon-compared to loop-retained gastrostomy tubes we observed a significantly higher probability of minor complications for the latter (p = 0.023).

CONCLUSION

As it is significantly less prone to minor complications we recommend the use of balloon-retained gastrostomy tubes to improve the practicability and maintenance of RIG.

VNTR diversity in Yersinia pestis isolates from an animal challenge study reveals the potential for in vitro mutations during laboratory cultivation

Underlying mutation rates and other evolutionary forces shape the population structure of bacteria in nature. Although easily overlooked, similar forces are at work in the laboratory and may influence observed mutations. Here, we investigated tissue samples and Yersinia pestis isolates from a rodent laboratory challenge with strain CO92 using whole genome sequencing and multi-locus variable-number tandem repeat (VNTR) analysis (MLVA). We identified six VNTR mutations that were found to have occurred in vitro during laboratory cultivation rather than in vivo during the rodent challenge. In contrast, no single nucleotide polymorphism (SNP) mutations were observed, either in vivo or in vitro. These results were consistent with previously published mutation rates and the calculated number of Y. pestis generations that occurred during the in vitro versus the in vivo portions of the experiment. When genotyping disease outbreaks, the potential for in vitro mutations should be considered, particularly when highly variable genetic markers such as VNTRs are used.

More than 50% of Clostridium difficile Isolates from Pet Dogs in Flagstaff, USA, Carry Toxigenic Genotypes

Nosocomial acquisition of Clostridium difficile is well documented, yet recent studies have highlighted the importance of community acquired infections and identified community associated reservoirs for this pathogen. Multiple studies have implicated companion pets and farm animals as possible sources of community acquired C. difficile infections in humans. To explore the potential role of pet dogs in human C. difficile infections we systematically collected canine fecal samples (n = 197) in Flagstaff, AZ. Additionally, nineteen fecal samples were collected at a local veterinary clinic from diarrheic dogs. We used these combined samples to investigate important questions regarding C. difficile colonization in pet canines: 1) What is the prevalence and diversity of C. difficile in this companion pet population, and 2) Do C. difficile isolates collected from canines genetically overlap with isolates that cause disease in humans? We used a two-step sequence typing approach, including multilocus sequence typing to determine the overall genetic diversity of C. difficile present in Flagstaff canines, and whole-genome sequencing to assess the fine-scale diversity patterns within identical multilocus sequence types from isolates obtained within and among multiple canine hosts. We detected C. difficile in 17% of the canine fecal samples with 10% containing toxigenic strains that are known to cause human disease. Sequencing analyses revealed similar genotypes in dogs and humans. These findings suggest that companion pets are a potential source of community acquired C. difficile infections in humans.

Genetic variation at the MHC DRB1 locus is similar across Gunnison's prairie dog (Cynomys gunnisoni) colonies regardless of plague history

Yersinia pestis was introduced to North America around 1900 and leads to nearly 100% mortality in prairie dog (Cynomys spp.) colonies during epizootic events, which suggests this pathogen may exert a strong selective force. We characterized genetic diversity at an MHC class II locus (DRB1) in Gunnison's prairie dog (C. gunnisoni) and quantified population genetic structure at the DRB1 versus 12 microsatellite loci in three large Arizona colonies. Two colonies, Seligman (SE) and Espee Ranch (ES), have experienced multiple plague‐related die‐offs in recent years, whereas plague has never been documented at Aubrey Valley (AV). We found fairly low allelic diversity at the DRB1 locus, with one allele (DRB1*01) at high frequency (0.67–0.87) in all colonies. Two other DRB1 alleles appear to be trans‐species polymorphisms shared with the black‐tailed prairie dog (C. ludovicianus), indicating that these alleles have been maintained across evolutionary time frames. Estimates of genetic differentiation were generally lower at the MHC locus (F_ST = 0.033) than at microsatellite markers (F_ST = 0.098). The reduced differentiation at DRB1 may indicate that selection has been important for shaping variation at MHC loci, regardless of the presence or absence of plague in recent decades. However, genetic drift has probably also influenced the DRB1 locus because its level of differentiation was not different from that of microsatellites in an F_ST outlier analysis. We then compared specific MHC alleles to plague survivorship in 60 C. gunnisoni that had been experimentally infected with Y. pestis. We found that survival was greater in individuals that carried at least one copy of the most common allele (DRB1*01) compared to those that did not (60% vs. 20%). Although the sample sizes of these two groups were unbalanced, this result suggests the possibility that this MHC class II locus, or a nearby linked gene, could play a role in plague survival.

Age-related differences in the activity of arterial mineral deposition and regional bone metabolism: a 18F-sodium fluoride positron emission tomography study

Functional (18)F-fluoride PET demonstrated an inverse relationship between the activity of arterial mineral deposition and regional bone metabolism. While bone metabolism decreases with age, the activity of arterial mineral deposition increases.

INTRODUCTION

The extent of arterial calcification increases with age, whereas bone mineral density decreases, evidencing a well-known inverse correlation on morphological basis. The aim of this study was to evaluate the functional relationship between the activity of arterial mineral deposition and regional bone metabolism as assessed by (18)F-sodium fluoride (NaF) PET/CT.

METHODS

Three hundred four subjects were examined by (18)F-NaF PET/CT. Tracer accumulation in the femoral arteries was analyzed both qualitatively and semiquantitatively by measuring the blood-pool-corrected standardized uptake value (target-to-background ratio). Uptake was compared with cardiovascular risk factors (RFs), calcified plaque burden, and regional bone metabolism as assessed by PET/CT.

RESULTS

The activity of arterial mineral deposition significantly increased with age (p < 0.001), whereas regional bone metabolism significantly decreased (p < 0.001). There was a significant inverse correlation between bone metabolism and arterial mineral deposition (unadjusted, p < 0.001); that association was not significant (p = 0.79) when controlled for age and other RFs. Both high activity of arterial mineral deposition and low bone metabolism were significantly associated with cardiovascular events and other RFs.

CONCLUSION

(18)F-NaF PET/CT provides a tool to visualize and quantify the activity of arterial mineral deposition and regional bone metabolism. In this study, we observed an inverse correlation between the activity of arterial mineral deposition and regional bone metabolism. While the activity of arterial mineral deposition significantly increases with age, regional bone metabolism decreases.

Multiple mutations in the para-sodium channel gene are associated with pyrethroid resistance in Rhipicephalus microplus from the United States and Mexico

BACKGROUND

Acaricide resistant Rhipicephalus microplus populations have become a major problem for many cattle producing areas of the world. Pyrethroid resistance in arthropods is typically associated with mutations in domains I, II, III, and IV of voltage-gated sodium channel genes. In R. microplus, known resistance mutations include a domain II change (C190A) in populations from Australia, Africa, and South America and a domain III mutation (T2134A) that only occurs in Mexico and the U.S.

METHODS

We investigated pyrethroid resistance in cattle fever ticks from Texas and Mexico by estimating resistance levels in field-collected ticks using larval packet discriminating dose (DD) assays and identifying single nucleotide polymorphisms (SNPs) in the para-sodium channel gene that associated with resistance. We then developed qPCR assays for three SNPs and screened a larger set of 1,488 R. microplus ticks, representing 77 field collections and four laboratory strains, for SNP frequency.

RESULTS

We detected resistance SNPs in 21 of 68 U.S. field collections and six of nine Mexico field collections. We expected to identify the domain III SNP (T2134A) at a high frequency; however, we only found it in three U.S. collections. A much more common SNP in the U.S. (detected in 19 of 21 field collections) was the C190A domain II mutation, which has never before been reported from North America. We also discovered a novel domain II SNP (T170C) in ten U.S. and two Mexico field collections. The T170C transition mutation has previously been associated with extreme levels of resistance (super-knockdown resistance) in insects. We found a significant correlation (r = 0.81) between the proportion of individuals in field collections that carried any two resistance SNPs and the percent survivorship of F1 larvae from these collections in DD assays. This relationship is accurately predicted by a simple linear regression model (R2 = 0.6635).

CONCLUSIONS

These findings demonstrate that multiple mutations in the para-sodium channel gene independently associate with pyrethroid resistance in R. microplus ticks, which is likely a consequence of human-induced selection.

Invasive potential of cattle fever ticks in the southern United States

Abstract'

Widespread movement of invasive cattle fever ticks (Rhipicephalus microplus) in southern Texas leads to shared local infestations on cattle and deer

Background

Rhipicephalus (Boophilus) microplus is a highly-invasive tick that transmits the cattle parasites (Babesia bovis and B. bigemina) that cause cattle fever. R. microplus and Babesia are endemic in Mexico and ticks persist in the United States inside a narrow tick eradication quarantine area (TEQA) along the Rio Grande. This containment area is threatened by unregulated movements of illegal cattle and wildlife like white-tailed deer (WTD; Odocoileus virginianus).

Methods

Using 11 microsatellite loci we genotyped 1,247 R. microplus from 63 Texas collections, including outbreak infestations from outside the TEQA. We used population genetic analyses to test hypotheses about ecological persistence, tick movement, and impacts of the eradication program in southern Texas. We tested acaricide resistance with larval packet tests (LPTs) on 47 collections.

Results

LPTs revealed acaricide resistance in 15/47 collections (32%); 11 were outside the TEQA and three were resistant to multiple acaricides. Some collections highly resistant to permethrin were found on cattle and WTD. Analysis of genetic differentiation over time at seven properties revealed local gene pools with very low levels of differentiation (F_ST 0.00-0.05), indicating persistence over timespans of up to 29 months. However, in one neighborhood differentiation varied greatly over a 12-month period (F_ST 0.03-0.13), suggesting recurring immigration from distinct sources as another persistence mechanism. Ticks collected from cattle and WTD at the same location are not differentiated (F_ST = 0), implicating ticks from WTD as a source of ticks on cattle (and vice versa) and emphasizing the importance of WTD to tick control strategies. We identified four major genetic groups (K = 4) using Bayesian population assignment, suggesting multiple introductions to Texas.

Conclusions

Two dispersal mechanisms give rise to new tick infestations: 1) frequent short-distance dispersal from the TEQA; and 2) rare long-distance, human-mediated dispersal from populations outside our study area, probably Mexico. The threat of cattle fever tick transport into Texas is increased by acaricide resistance and the ability of R. microplus to utilize WTD as an alternate host. Population genetic analyses may provide a powerful tool for tracking invasions in other parts of the world where these ticks are established.

Accurate and rapid identification of the Burkholderia pseudomallei near-neighbour, Burkholderia ubonensis, using real-time PCR

Burkholderia ubonensis is an environmental bacterium belonging to the Burkholderia cepacia complex (Bcc), a group of genetically related organisms that are associated with opportunistic but generally nonfatal infections in healthy individuals. In contrast, the near-neighbour species Burkholderia pseudomallei causes melioidosis, a disease that can be fatal in up to 95% of cases if left untreated. B. ubonensis is frequently misidentified as B. pseudomallei from soil samples using selective culturing on Ashdown’s medium, reflecting both the shared environmental niche and morphological similarities of these species. Additionally, B. ubonensis shows potential as an important biocontrol agent in B. pseudomallei-endemic regions as certain strains possess antagonistic properties towards B. pseudomallei. Current methods for characterising B. ubonensis are laborious, time-consuming and costly, and as such this bacterium remains poorly studied. The aim of our study was to develop a rapid and inexpensive real-time PCR-based assay specific for B. ubonensis. We demonstrate that a novel B. ubonensis-specific assay, Bu550, accurately differentiates B. ubonensis from B. pseudomallei and other species that grow on selective Ashdown’s agar. We anticipate that Bu550 will catalyse research on B. ubonensis by enabling rapid identification of this organism from Ashdown’s-positive colonies that are not B. pseudomallei.

Hybridization in natural sympatric populations of Dermacentor ticks in northwestern North America

Hybridization in ticks has been described in a handful of species and mostly as a result of laboratory experiments. We used 148 AFLP loci to describe putative hybridization events between D. andersoni and D. variabilis in sympatric populations from northwestern North America. Recently, D. variabilis has expanded its range westward into the natural range of D. andersoni. Using a sample of 235 D. andersoni and 62 D. variabilis, we identified 31 individuals as putative hybrids: four F2 individuals and 27 backcrosses to D. andersoni (as defined by NewHybrids). We found no evidence of hybrids backcrossing into D. variabilis. Furthermore, all hybrids presented 16S mtDNA signatures characteristic of D. andersoni, which indicates the directionality of the hybrid crosses: female D. andersoni × male D. variabilis. We also discovered 13 species-specific AFLP fragments for D. andersoni. These loci were found to have a decreased occurrence in the putative hybrids and were absent altogether in D. variabilis samples. AFLP profiles were also used to determine the levels of genetic population structure and gene flow among nine populations of D. andersoni and three of D. variabilis. Genetic structure exists in both species (D. andersoni, ΦST = 0.110; D. variabilis, ΦST = 0.304) as well as significant estimates of isolation by distance (D. andersoni, ρ = 0.066, P = 0.001; D. variabilis, ρ = 0.729, P = 0.001).

Genetic characterization of Theileria equi infecting horses in North America: evidence for a limited source of U.S. introductions

Background

Theileria equi is a tick-borne apicomplexan hemoparasite that causes equine piroplasmosis. This parasite has a worldwide distribution but the United States was considered to be free of this disease until recently.

Methods

We used samples from 37 horses to determine genetic relationships among North American T. equi using the 18S rRNA gene and microsatellites. We developed a DNA fingerprinting panel of 18 microsatellite markers using the first complete genome sequence of T. equi.

Results

A maximum parsimony analysis of 18S rRNA sequences grouped the samples into two major clades. The first clade (n = 36) revealed a high degree of nucleotide similarity in U.S. T. equi, with just 0–2 single nucleotide polymorphisms (SNPs) among samples. The remaining sample fell into a second clade that was genetically divergent (48 SNPs) from the other U.S. samples. This sample was collected at the Texas border, but may have originated in Mexico. We genotyped T. equi from the U.S. using microsatellite markers and found a moderate amount of genetic diversity (2–8 alleles per locus). The field samples were mostly from a 2009 Texas outbreak (n = 22) although samples from five other states were also included in this study. Using Weir and Cockerham’s F_ST estimator (θ) we found strong population differentiation of the Texas and Georgia subpopulations (θ = 0.414), which was supported by a neighbor-joining tree created with predominant single haplotypes. Single-clone infections were found in 27 of the 37 samples (73%), allowing us to identify 15 unique genotypes.

Conclusions

The placement of most T. equi into one monophyletic clade by 18S is suggestive of a limited source of introduction into the U.S. When applied to a broader cross section of worldwide samples, these molecular tools should improve source tracking of T. equi outbreaks and may help prevent the spread of this tick-borne parasite.

FungiQuant: a broad-coverage fungal quantitative real-time PCR assay

BACKGROUND

Fungal load quantification is a critical component of fungal community analyses. Limitation of current approaches for quantifying the fungal component in the human microbiome suggests the need for new broad-coverage techniques.

METHODS

We analyzed 2,085 18S rRNA gene sequences from the SILVA database for assay design. We generated and quantified plasmid standards using a qPCR-based approach. We evaluated assay coverage against 4,968 sequences and performed assay validation following the Minimum Information for Publication of Quantitative Real-Time PCR Experiments (MIQE) guidelines.

RESULTS

We designed FungiQuant, a TaqMan® qPCR assay targeting a 351 bp region in the fungal 18S rRNA gene. Our in silico analysis showed that FungiQuant is a perfect sequence match to 90.0% of the 2,617 fungal species analyzed. We showed that FungiQuant's is 100% sensitive and its amplification efficiencies ranged from 76.3% to 114.5%, with r(2)-values of >0.99 against the 69 fungal species tested. Additionally, FungiQuant inter- and intra-run coefficients of variance ranged from <10% and <20%, respectively. We further showed that FungiQuant has a limit of quantification 25 copies and a limit of detection at 5 copies. Lastly, by comparing results from human-only background DNA with low-level fungal DNA, we showed that amplification in two or three of a FungiQuant performed in triplicate is statistically significant for true positive fungal detection.

CONCLUSIONS

FungiQuant has comprehensive coverage against diverse fungi and is a robust quantification and detection tool for delineating between true fungal detection and non-target human DNA.

Phylogeography of Francisella tularensis subsp. holarctica, Europe

Francisella tularensis subsp. holarctica isolates from Austria, Germany, Hungary, Italy, and Romania were placed into an existing phylogeographic framework. Isolates from Italy were assigned to phylogenetic group B.FTNF002–00; the other isolates, to group B.13. Most F. tularensis subsp. holarctica isolates from Europe belong to these 2 geographically segregated groups.

Epidemiological tracking and population assignment of the non-clonal bacterium, Burkholderia pseudomallei

Rapid assignment of bacterial pathogens into predefined populations is an important first step for epidemiological tracking. For clonal species, a single allele can theoretically define a population. For non-clonal species such as Burkholderia pseudomallei, however, shared allelic states between distantly related isolates make it more difficult to identify population defining characteristics. Two distinct B. pseudomallei populations have been previously identified using multilocus sequence typing (MLST). These populations correlate with the major foci of endemicity (Australia and Southeast Asia). Here, we use multiple Bayesian approaches to evaluate the compositional robustness of these populations, and provide assignment results for MLST sequence types (STs). Our goal was to provide a reference for assigning STs to an established population without the need for further computational analyses. We also provide allele frequency results for each population to enable estimation of population assignment even when novel STs are discovered. The ability for humans and potentially contaminated goods to move rapidly across the globe complicates the task of identifying the source of an infection or outbreak. Population genetic dynamics of B. pseudomallei are particularly complicated relative to other bacterial pathogens, but the work here provides the ability for broad scale population assignment. As there is currently no independent empirical measure of successful population assignment, we provide comprehensive analytical details of our comparisons to enable the reader to evaluate the robustness of population designations and assignments as they pertain to individual research questions. Finer scale subdivision and verification of current population compositions will likely be possible with genotyping data that more comprehensively samples the genome. The approach used here may be valuable for other non-clonal pathogens that lack simple group-defining genetic characteristics and provides a rapid reference for epidemiologists wishing to track the origin of infection without the need to compile population data and learn population assignment algorithms.

Burkholderia pseudomallei is a soil-dwelling bacterium that can infect a large range of hosts. In humans, B. pseudomallei causes melioidosis, and typical routes of entry include open wounds, inhalation, or ingestion. Clinical features are diverse, although pneumonia and abscess formation are common. High rates of recombination within the genome of this bacterium have confounded attempts to match clinical samples to geographically defined populations. Here we provide a reference that simplifies source attribution issues. We applied population assignment software to previously generated sequence data from seven B. pseudomallei genes to define the major geographic populations within this species. We evaluated the robustness of our results by comparison with two additional population assignment programs. We present the likelihood that each variant is assigned to a particular geographic population. This information can be used to assign novel B. pseudomallei isolates to a geographic population without needing to learn and run cumbersome population assignment applications. This method can also be used for other bacteria that are difficult to source-attribute due to high levels of genomic variation and recombination.

Resistance to plague among black-tailed prairie dog populations

In some rodent species frequently exposed to plague outbreaks caused by Yersinia pestis, resistance to the disease has evolved as a population trait. As a first step in determining if plague resistance has developed in black-tailed prairie dogs (Cynomys ludovicianus), animals captured from colonies in a plague-free region (South Dakota) and two plague-endemic regions (Colorado and Texas) were challenged with Y. pestis at one of three doses (2.5, 250, or 2500 mouse LD50s). South Dakota prairie dogs were far more susceptible to plague than Colorado and Texas prairie dogs (p<0.001), with a mortality rate of nearly 100% over all doses. Colorado and Texas prairie dogs were quite similar in their response, with overall survival rates of 50% and 60%, respectively. Prairie dogs from these states were heterogeneous in their response, with some animals dying at the lowest dose (37% and 20%, respectively) and some surviving even at the highest dose (29% and 40%, respectively). Microsatellite analysis revealed that all three groups were distinct genetically, but further studies are needed to establish a genetic basis for the observed differences in plague resistance.

Phylogeography of Francisella tularensis subspecies holarctica from the country of Georgia

Background

Francisella tularensis, the causative agent of tularemia, displays subspecies-specific differences in virulence, geographic distribution, and genetic diversity. F. tularensis subsp. holarctica is widely distributed throughout the Northern Hemisphere. In Europe, F. tularensis subsp. holarctica isolates have largely been assigned to two phylogenetic groups that have specific geographic distributions. Most isolates from Western Europe are assigned to the B.Br.FTNF002-00 group, whereas most isolates from Eastern Europe are assigned to numerous lineages within the B.Br.013 group. The eastern geographic extent of the B.Br.013 group is currently unknown due to a lack of phylogenetic knowledge about populations at the European/Asian juncture and in Asia. In this study, we address this knowledge gap by describing the phylogenetic structure of F. tularensis subsp. holarctica isolates from the country of Georgia, and by placing these isolates into a global phylogeographic context.

Results

We identified a new genetic lineage of F. tularensis subsp. holarctica from Georgia that belongs to the B.Br.013 group. This new lineage is genetically and geographically distinct from lineages previously described from the B.Br.013 group from Central-Eastern Europe. Importantly, this new lineage is basal within the B.Br.013 group, indicating the Georgian lineage diverged before the diversification of the other known B.Br.013 lineages. Although two isolates from the Georgian lineage were collected nearby in the Ukrainian region of Crimea, all other global isolates assigned to this lineage were collected in Georgia. This restricted geographic distribution, as well as the high levels of genetic diversity within the lineage, is consistent with a relatively older origin and localized differentiation.

Conclusions

We identified a new lineage of F. tularensis subsp. holarctica from Georgia that appears to have an older origin than any other diversified lineages previously described from the B.Br.013 group. This finding suggests that additional phylogenetic studies of F. tularensis subsp. holarctica populations in Eastern Europe and Asia have the potential to yield important new insights into the evolutionary history and phylogeography of this broadly dispersed F. tularensis subspecies.

Identifying experimental surrogates for Bacillus anthracis spores: a review

Bacillus anthracis, the causative agent of anthrax, is a proven biological weapon. In order to study this threat, a number of experimental surrogates have been used over the past 70 years. However, not all surrogates are appropriate for B. anthracis, especially when investigating transport, fate and survival. Although B. atrophaeus has been widely used as a B. anthracis surrogate, the two species do not always behave identically in transport and survival models. Therefore, we devised a scheme to identify a more appropriate surrogate for B. anthracis. Our selection criteria included risk of use (pathogenicity), phylogenetic relationship, morphology and comparative survivability when challenged with biocides. Although our knowledge of certain parameters remains incomplete, especially with regards to comparisons of spore longevity under natural conditions, we found that B. thuringiensis provided the best overall fit as a non-pathogenic surrogate for B. anthracis. Thus, we suggest focusing on this surrogate in future experiments of spore fate and transport modelling.

Recent demographic bottlenecks are not accompanied by a genetic signature in banner-tailed kangaroo rats (Dipodomys spectabilis)

Single-sample methods of bottleneck detection are now routine analyses in studies of wild populations and conservation genetics. Three common approaches to bottleneck detection are the heterozygosity excess, mode-shift, and M-ratio tests. Empirical groundtruthing of these methods is difficult, but their performances are critical for the accurate reconstruction of population demography. We use two banner-tailed kangaroo rat (Dipodomys spectabilis) populations from southeastern Arizona (USA) that are known to have experienced recent demographic reductions to search for genetic bottleneck signals with eight microsatellite loci. Over eight total sample-years, neither population showed a genetic bottleneck signature. M-ratios in both populations were large, stable, and never fell below a critical significance value (Mc). The mode shift test did not detect any distortion of allele frequencies, and tests of heterozygosity excess were not significant in postbottleneck samples when we used standard microsatellite mutation models. The genetic effects of bottlenecks like those experienced by our study populations should be strongly influenced by rates of mutation and migration. We used genetic parentage data to estimate a relatively high mutation rate in D. spectabilis (0.0081 mutants/generation/locus), but mutation alone is unlikely to explain the temporal distribution of rare alleles that we observed. Migration (gene flow) is a more likely explanation, despite prior mark-recapture analysis that estimated very low rates of interpopulation dispersal. We interpret our kangaroo rat data in light of the broader literature and conclude that in natural populations connected by dispersal, demographic bottlenecks may prove difficult to detect using molecular genetic data.

Strain-specific single-nucleotide polymorphism assays for the Bacillus anthracis Ames strain

Highly precise diagnostics and forensic assays can be developed through a combination of evolutionary analysis and the exhaustive examination of genomic sequences. In Bacillus anthracis, whole-genome sequencing efforts revealed ca. 3,500 single-nucleotide polymorphisms (SNPs) among eight different strains and evolutionary analysis provides the identification of canonical SNPs. We have previously shown that SNPs are highly evolutionarily stable, and the clonal nature of B. anthracis makes them ideal signatures for subtyping this pathogen. Here we identified SNPs that define the lineage of B. anthracis that contains the Ames strain, the strain used in the 2001 bioterrorist attacks in the United States. Sequencing and real-time PCR were used to validate these SNPs across B. anthracis strains, including (i) 88 globally and genetically diverse isolates; (ii) isolates that were shown to be genetic relatives of the Ames strain by multiple-locus variable number tandem repeat analysis (MLVA); and (iii) several different lab stocks of the Ames strain, including a clinical isolate from the 2001 letter attack. Six SNPs were found to be highly specific for the Ames strain; four on the chromosome, one on the pX01 plasmid, and one on the pX02 plasmid. All six SNPs differentiated the B. anthracis Ames strain from the 88 unique B. anthracis strains, while five of the six separated Ames from its close genetic relatives. The use of these SNPs coupled with real-time PCR allows specific and sensitive (<100 fg of template DNA) identification of the Ames strain. This evolutionary and genomics-based approach provides an effective means for the discovery of strain-specific SNPs in B. anthracis.

Parentage analysis detects cryptic precapture dispersal in a philopatric rodent

Locating birthplaces using genetic parentage determination can increase the precision and accuracy with which animal dispersal patterns are established. We re-analyse patterns of movement away from the birthplace as a function of time, sex and population density for a sample of 303 banner-tailed kangaroo rats, Dipodomys spectabilis. We located birth sites using a combination of likelihood-based parentage analysis with live-trapping of mothers during the breeding season. The results demonstrate that natal-breeding site distances are density dependent in this species; in particular, both sexes emigrate earlier in the year, and females disperse farther than males, at low population densities. Banner-tailed kangaroo rats were chosen as a study system because live-trapping easily detects maternal and offspring locations; nevertheless, parentage analysis reveals that some offspring evade early detection and move substantial distances before their first capture. In a few cases, the approach even detects dispersal out of the natal 'deme' prior to first capture. Parentage analysis confirms the extreme philopatry of both sexes but indicates that prior estimates of median dispersal distance were too low. For D. spectabilis, more accurate location of individual birthplaces clarifies patterns of sex bias and density dependence in dispersal, and may resolve apparent discrepancies between direct and indirect estimates of dispersal distance. For species in which mothers can be more reliably trapped than juveniles, using offspring genotypes to locate parents is a novel way that genetic techniques can contribute to the analysis of animal dispersal.

Mass spectrometry provides accurate characterization of two genetic marker types in Bacillus anthracis

Epidemiological and forensic analyses of bioterrorism events involving Bacillus anthracis could be improved if both variable number tandem repeats (VNTRs) and single nucleotide polymorphisms (SNPs) could be combined on a single analysis platform. Here we present the use of electrospray ionization Fourier transform ion cyclotron resonance mass spectrometry (ESI-FTICR-MS) to characterize 24 alleles from 6 VNTR loci and 11 alleles from 7 SNP loci in B. anthracis. The results obtained with ESI-FTICR-MS were consistent with independent results obtained from traditional approaches using electrophoretic detection of fluorescent products. However, ESI-FTICR-MS improves on the traditional approaches because it does not require fluorescent labeling of PCR products, minimizes post-PCR processing, obviates electrophoresis, and provides unambiguous base composition of both SNP and VNTR PCR products. In addition, ESI-FTICR-MS allows both marker types to be examined simultaneously and at a rate of approximately 1 sample per min. This technology represents a significant advance in our ability to rapidly characterize B. anthracis isolates using VNTR and SNP loci.

Phylogenetic discovery bias in Bacillus anthracis using single-nucleotide polymorphisms from whole-genome sequencing

Phylogenetic reconstruction using molecular data is often subject to homoplasy, leading to inaccurate conclusions about phylogenetic relationships among operational taxonomic units. Compared with other molecular markers, single-nucleotide polymorphisms (SNPs) exhibit extremely low mutation rates, making them rare in recently emerged pathogens, but they are less prone to homoplasy and thus extremely valuable for phylogenetic analyses. Despite their phylogenetic potential, ascertainment bias occurs when SNP characters are discovered through biased taxonomic sampling; by using whole-genome comparisons of five diverse strains of Bacillus anthracis to facilitate SNP discovery, we show that only polymorphisms lying along the evolutionary pathway between reference strains will be observed. We illustrate this in theoretical and simulated data sets in which complex phylogenetic topologies are reduced to linear evolutionary models. Using a set of 990 SNP markers, we also show how divergent branches in our topologies collapse to single points but provide accurate information on internodal distances and points of origin for ancestral clades. These data allowed us to determine the ancestral root of B. anthracis, showing that it lies closer to a newly described "C" branch than to either of two previously described "A" or "B" branches. In addition, subclade rooting of the C branch revealed unequal evolutionary rates that seem to be correlated with ecological parameters and strain attributes. Our use of nonhomoplastic whole-genome SNP characters allows branch points and clade membership to be estimated with great precision, providing greater insight into epidemiological, ecological, and forensic questions.