THE INNATE IMMUNE RESPONSE MAY BE IMPORTANT FOR SURVIVING PLAGUE IN WILD GUNNISON'S PRAIRIE DOGS

Polygenic plague resistance in the great gerbil uncovered by population sequencing

Pathogens can elicit high selective pressure on hosts, potentially altering genetic diversity over short evolutionary timescales. Intraspecific variation in immune response is observable as variable survivability from specific infections. The great gerbil (Rhombomys opimus) is a rodent plague host with a heterogenic but highly resistant phenotype. Here, we investigate the genomic basis for plague-resistant phenotypes by exposing wild-caught great gerbils to plague (Yersinia pestis). Whole genome sequencing of 10 survivors and 10 moribund individuals revealed a subset of genomic regions showing elevated differentiation. Gene ontology analysis of candidate genes in these regions demonstrated enrichment of genes directly involved in immune functions, cellular metabolism and the regulation of apoptosis as well as pathways involved in transcription, translation, and gene regulation. Transcriptomic analysis revealed that the early activated great gerbil immune response to plague consisted of classical components of the innate immune system. Our approach combining challenge experiments with transcriptomics and population level sequencing, provides new insight into the genetic background of plague-resistance and confirms its complex nature, most likely involving multiple genes and pathways of both the immune system and regulation of basic cellular functions.

Reevaluation of the Role of Blocked Oropsylla hirsuta Prairie Dog Fleas (Siphonaptera: Ceratophyllidae) in Yersinia pestis (Enterobacterales: Enterobacteriaceae) Transmission

Prairie dogs in the western United States experience periodic epizootics of plague, caused by the flea-borne bacterial pathogen Yersinia pestis. An early study indicated that Oropsylla hirsuta (Baker), often the most abundant prairie dog flea vector of plague, seldom transmits Y. pestis by the classic blocked flea mechanism. More recently, an alternative early-phase mode of transmission has been proposed as the driving force behind prairie dog epizootics. In this study, using the same flea infection protocol used previously to evaluate early-phase transmission, we assessed the vector competence of O. hirsuta for both modes of transmission. Proventricular blockage was evident during the first two weeks after infection and transmission during this time was at least as efficient as early-phase transmission 2 d after infection. Thus, both modes of transmission likely contribute to plague epizootics in prairie dogs.

Plague reservoir species throughout the world

Plague has been known since ancient times as a re-emerging infectious disease, causing considerable socioeconomic burden in regional hotspots. To better understand the epidemiological cycle of the causative agent of the plague, its potential occurrence, and possible future dispersion, one must carefully consider the taxonomy, distribution, and ecological requirements of reservoir-species in relation either to natural or human-driven changes (e.g. climate change or urbanization). In recent years, the depth of knowledge on species taxonomy and species composition in different landscapes has undergone a dramatic expansion, driven by modern taxonomic methods such as synthetic surveys that take into consideration morphology, genetics, and the ecological setting of captured animals to establish their species identities. Here, we consider the recent taxonomic changes of the rodent species in known plague reservoirs and detail their distribution across the world, with a particular focus on those rodents considered to be keystone host species. A complete checklist of all known plague-infectable vertebrates living in plague foci is provided as a Supporting Information table.

First Genome Sequence of the Gunnison's Prairie Dog (Cynomys gunnisoni), a Keystone Species and Player in the Transmission of Sylvatic Plague

Prairie dogs (genus Cynomys) are a charismatic symbol of the American West. Their large social aggregations and complex vocalizations have been the subject of scientific and popular interest for decades. A large body of literature has documented their role as keystone species of western North America's grasslands: They generate habitat for other vertebrates, increase nutrient availability for plants, and act as a food source for mammalian, squamate, and avian predators. An additional keystone role lies in their extreme susceptibility to sylvatic plague (caused by Yersinia pestis), which results in periodic population extinctions, thereby generating spatiotemporal heterogeneity in both biotic communities and ecological processes. Here, we report the first Cynomys genome for a Gunnison's prairie dog (C. gunnisoni gunnisoni) from Telluride, Colorado (USA). The genome was constructed using a hybrid assembly of PacBio and Illumina reads and assembled with MaSuRCA and PBJelly, which resulted in a scaffold N50 of 824 kb. Total genome size was 2.67 Gb, with 32.46% of the bases occurring in repeat regions. We recovered 94.9% (91% complete) of the single copy orthologs using the mammalian Benchmarking Universal Single-Copy Orthologs database and detected 49,377 gene models (332,141 coding regions). Pairwise Sequentially Markovian Coalescent showed support for long-term stable population size followed by a steady decline beginning near the end of the Pleistocene, as well as a recent population reduction. The genome will aid in studies of mammalian evolution, disease resistance, and the genomic basis of life history traits in ground squirrels.

The Genome of the Great Gerbil Reveals Species-Specific Duplication of an MHCII Gene

The great gerbil (Rhombomys opimus) is a social rodent living in permanent, complex burrow systems distributed throughout Central Asia, where it serves as the main host of several important vector-borne infectious pathogens including the well-known plague bacterium (Yersinia pestis). Here, we present a continuous annotated genome assembly of the great gerbil, covering over 96% of the estimated 2.47-Gb genome. Taking advantage of the recent genome assemblies of the sand rat (Psammomys obesus) and the Mongolian gerbil (Meriones unguiculatus), comparative immunogenomic analyses reveal shared gene losses within TLR gene families (i.e., TLR8, TLR10, and the entire TLR11-subfamily) for Gerbillinae, accompanied with signs of diversifying selection of TLR7 and TLR9. Most notably, we find a great gerbil-specific duplication of the MHCII DRB locus. In silico analyses suggest that the duplicated gene provides high peptide binding affinity for Yersiniae epitopes as well as Leishmania and Leptospira epitopes, putatively leading to increased capability to withstand infections by these pathogens. Our study demonstrates the power of whole-genome sequencing combined with comparative genomic analyses to gain deeper insight into the immunogenomic landscape of the great gerbil and its close relatives.

Differential plague susceptibility in species and populations of prairie dogs

Laboratory trials conducted over the past decade at U.S. Geological Survey National Wildlife Health Center indicate that wild populations of prairie dogs (Cynomys spp.) display different degrees of susceptibility to experimental challenge with fully virulent Yersinia pestis, the causative agent of plague. We evaluated patterns in prairie dog susceptibility to plague to determine whether the historical occurrence of plague at location of capture was related to survival times of prairie dogs challenged with Y. pestis. We found that black-tailed prairie dogs (Cynomys ludovicianus) from South Dakota (captured prior to the detection of plague in the state), Gunnison's prairie dogs (Cynomys gunnisoni) from Colorado, and Utah prairie dogs (Cynomys parvidens) from Utah were most susceptible to plague. Though the susceptibility of black-tailed prairie dogs in South Dakota compared with western locations supports our hypothesis regarding historical exposure, both Colorado and Utah prairie dogs have a long history of exposure to plague. It is possible that for these populations, genetic isolation/bottle necks have made them more susceptible to plague outbreaks.

Immunology of Yersinia pestis Infection

As a pathogen of plague, Yersinia pestis caused three massive pandemics in history that killed hundreds of millions of people. Yersinia pestis is highly invasive, causing severe septicemia which, if untreated, is usually fatal to its host. To survive in the host and maintain a persistent infection, Yersinia pestis uses several stratagems to evade the innate and the adaptive immune responses. For example, infections with this organism are biphasic, involving an initial "noninflammatory" phase where bacterial replication occurs initially with little inflammation and following by extensive phagocyte influx, inflammatory cytokine production, and considerable tissue destruction, which is called "proinflammatory" phase. In contrast, the host also utilizes its immune system to eliminate the invading bacteria. Neutrophil and macrophage are the first defense against Yersinia pestis invading through phagocytosis and killing. Other innate immune cells also play different roles, such as dendritic cells which help to generate more T helper cells. After several days post infection, the adaptive immune response begins to provide organism-specific protection and has a long-lasting immunological memory. Thus, with the cooperation and collaboration of innate and acquired immunity, the bacterium may be eliminated from the host. The research of Yersinia pestis and host immune systems provides an important topic to understand pathogen-host interaction and consequently develop effective countermeasures.

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.

Evaluation of Yersinia pestis Transmission Pathways for Sylvatic Plague in Prairie Dog Populations in the Western U.S

Sylvatic plague, caused by the bacterium Yersinia pestis, is periodically responsible for large die-offs in rodent populations that can spillover and cause human mortalities. In the western US, prairie dog populations experience nearly 100% mortality during plague outbreaks, suggesting that multiple transmission pathways combine to amplify plague dynamics. Several alternate pathways in addition to flea vectors have been proposed, such as transmission via direct contact with bodily fluids or inhalation of infectious droplets, consumption of carcasses, and environmental sources of plague bacteria, such as contaminated soil. However, evidence supporting the ability of these proposed alternate pathways to trigger large-scale epizootics remains elusive. Here we present a short review of potential plague transmission pathways and use an ordinary differential equation model to assess the contribution of each pathway to resulting plague dynamics in black-tailed prairie dogs (Cynomys ludovicianus) and their fleas (Oropsylla hirsuta). Using our model, we found little evidence to suggest that soil contamination was capable of producing plague epizootics in prairie dogs. However, in the absence of flea transmission, direct transmission, i.e., contact with bodily fluids or inhalation of infectious droplets, could produce enzootic dynamics, and transmission via contact with or consumption of carcasses could produce epizootics. This suggests that these pathways warrant further investigation.

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.