Anthropogenic habitat alteration leads to rapid loss of adaptive variation and restoration potential in wild salmon populations.

Phenotypic variation is critical for the long-term persistence of species and populations. Anthropogenic activities have caused substantial shifts and reductions in phenotypic variation across diverse taxa, but the underlying mechanism(s) (i.e., phenotypic plasticity and/or genetic evolution) and long-term consequences (e.g., ability to recover phenotypic variation) are unclear. Here we investigate the widespread and dramatic changes in adult migration characteristics of wild Chinook salmon caused by dam construction and other anthropogenic activities. Strikingly, we find an extremely robust association between migration phenotype (i.e., spring-run or fall-run) and a single locus, and that the rapid phenotypic shift observed after a recent dam construction is explained by dramatic allele frequency change at this locus. Furthermore, modeling demonstrates that continued selection against the spring-run phenotype could rapidly lead to complete loss of the spring-run allele, and an empirical analysis of populations that have already lost the spring-run phenotype reveals they are not acting as sustainable reservoirs of the allele. Finally, ancient DNA analysis suggests the spring-run allele was abundant in historical habitat that will soon become accessible through a large-scale restoration (i.e., dam removal) project, but our findings suggest that widespread declines and extirpation of the spring-run phenotype and allele will challenge reestablishment of the spring-run phenotype in this and future restoration projects. These results reveal the mechanisms and consequences of human-induced phenotypic change and highlight the need to conserve and restore critical adaptive variation before the potential for recovery is lost.

Temporal dynamics of migration-linked genetic variation are driven by streamflows and riverscape permeability.

Landscape permeability is often explored spatially, but may also vary temporally. Landscape permeability, including partial barriers, influences migratory animals that move across the landscape. Partial barriers are common in rivers where barrier passage varies with streamflow. We explore the influence of partial barriers on the spatial and temporal distribution of migration-linked genotypes of Oncorhynchus mykiss, a salmonid fish with co-occurring resident and migratory forms, in tributaries to the South Fork Eel River, California, USA, Elder and Fox Creeks. We genotyped >4,000 individuals using RAD-capture and classified individuals as resident, heterozygous or migratory genotypes using life history-associated loci. Across four years of study (2014-2017), the permeability of partial barriers varied across dry and wet years. In Elder Creek, the largest waterfall was passable for adults migrating up-river 4-39 days each year. In this stream, the overall spatial pattern, with fewer migratory genotypes above the waterfall, remained true across dry and wet years (67%-76% of migratory alleles were downstream of the waterfall). We also observed a strong relationship between distance upstream and proportion of migratory alleles. In Fox Creek, the primary barrier is at the mouth, and we found that the migratory allele frequency varied with the annual timing of high flow events. In years when rain events occurred during the peak breeding season, migratory allele frequency was high (60%-68%), but otherwise it was low (30% in two years). We highlight that partial barriers and landscape permeability can be temporally dynamic, and this effect can be observed through changing genotype frequencies in migratory animals.

Wildfire reveals transient changes to individual traits and population responses of a native bumble bee Bombus vosnesenskii.

Fire-induced changes in the abundance and distribution of organisms, especially plants, can alter resource landscapes for mobile consumers driving bottom-up effects on their population sizes, morphologies and reproductive potential. We expect these impacts to be most striking for obligate visitors of plants, like bees and other pollinators, but these impacts can be difficult to interpret due to the limited information provided by forager counts in the absence of survival or fitness proxies. Increased bumble bee worker abundance is often coincident with the pulses of flowers that follow recent fire. However, it is unknown if observed postfire activity is due to underlying population growth or a stable pool of colonies recruiting more foragers to abundant resource patches. This distinction is necessary for determining the net impact of disturbance on bumble bees: are there population-wide responses or do just a few colonies reap the rewards? We estimated colony abundance before and after fire in burned and unburned areas using a genetic mark-recapture framework. We paired colony abundance estimates with measures of body size, counts of queens, and estimates of foraging and dispersal to assess changes in worker size, reproductive output, and landscape-scale movements. Higher floral abundance following fire not only increased forager abundance but also the number of colonies from which those foragers came. Importantly, despite a larger population size, we also observed increased mean worker size. Two years following fire, queen abundance was higher in both burned and unburned sites, potentially due to the dispersal of queens from burned into unburned areas. The effects of fire were transient; within two growing seasons, worker abundance was substantially reduced across the entire sampling area and body sizes were similar between burned and unburned sites. Our results reveal how disturbance can temporarily release populations from resource limitation, boosting the genetic diversity, body size, and reproductive output of populations. Given that the effects of fire on bumble bees acted indirectly through pulsed resource availability, it is likely our results are generalizable to other situations, such as habitat restorations, where resource density is enhanced within the landscape.

On the Ecology and Distribution of Steelhead (Oncorhynchus mykiss) in California's Eel River.

The preservation of life history and other phenotypic complexity is central to the resilience of Pacific salmon stocks. Steelhead (Oncorhynchus mykiss) express a diversity of life-history strategies such as the propensity to migrate (anadromy/residency) and the timing and state of maturation upon return to freshwater (run-timing), providing an opportunity to study adaptive phenotypic complexity. Historically, the Eel River supported upwards of 1 million salmon and steelhead, but the past century has seen dramatic declines of all salmonids in the watershed. Here we investigate life-history variation in Eel River steelhead by using Rapture sequencing, on thousands of individuals, to genotype the region diagnostic for run-timing (GREB1L) and the region strongly associated with residency/anadromy (OMY5) in the Eel River and other locations, as well as determine patterns of overall genetic differentiation. Our results provide insight into many conservation-related issues. For example, we found that distinct segregation between winter and summer-run steelhead correlated with flow-dependent barriers in major forks of the Eel, that summer-run steelhead inhabited the upper Eel prior to construction of an impassable dam, and that both life history and overall genetic diversity have been maintained in the resident trout population above; and we found no evidence of the summer-run allele in the South Fork Eel, indicating that summer run-timing cannot be expected to arise from standing genetic variation in this and other populations that lack the summer-run phenotype. The results presented in this study provide valuable information for designing future restoration and management strategies for O. mykiss in Northern California and beyond.

Hybridization between two parapatric ranid frog species in the northern Sierra Nevada, California, USA.

Contact zones between species provide a unique opportunity to test whether taxa can hybridize or not. Cross-breeding or hybridization between closely related taxa can promote gene flow (introgression) between species, adaptation, or even speciation. Though hybridization events may be short-lived and difficult to detect in the field, genetic data can provide information about the level of introgression between closely related taxa. Hybridization can promote introgression between species, which may be an important evolutionary mechanism for either homogenization (reversing initial divergence between species) or reproductive isolation (potentially leading to speciation). Here, we used thousands of genetic markers from nuclear DNA to detect hybridization between two parapatric frog species (Rana boylii and Rana sierrae) in the Sierra Nevada of California. Based on principal components analysis, admixture, and analysis of heterozygosity at species diagnostic SNPs, we detected two F1 hybrid individuals in the Feather River basin, as well as a weak signal of introgression and gene flow between the frog species compared with frog populations from two other adjacent watersheds. This study provides the first documentation of hybridization and introgression between these two species, which are of conservation concern.

Best available science still supports an ancient common origin of Devils Hole and Devils Hole pupfish.

The age of DHP and how pupfish colonized Devils Hole have always been a topic of interest. Recently, two different publications (Martin, Crawford, Turner, & Simons, & Saglam et al., ) tackled this issue using genomic data sets and demographic models but came to widely different conclusions. In their comment, Martin and Höhne () argue that our results (Saglam et al., ) were misleading because we used inappropriate calibration information and biased a priori assumptions. They then re-analysed our data using a "biologically informed" mutation rate prior and concluded that our data support a much younger age of DHP (12.6 kya) as opposed to 60 kya reported in our study. Below we will summarize why their arguments do not hold up and explore some of the inconsistencies between their claims and what was actually presented in our study. Furthermore, we will demonstrate their re-analyses provide no new information compared to what was presented in our original manuscript and reinforce our estimate of a 60 kya divergence of DHP as outweighing competing hypotheses.

Large-effect loci affect survival in Tasmanian devils (Sarcophilus harrisii) infected with a transmissible cancer.

Identifying the genetic architecture of complex phenotypes is a central goal of modern biology, particularly for disease-related traits. Genome-wide association methods are a classical approach for identifying the genomic basis of variation in disease phenotypes, but such analyses are particularly challenging in natural populations due to sample size difficulties. Extensive mark-recapture data, strong linkage disequilibrium and a lethal transmissible cancer make the Tasmanian devil (Sarcophilus harrisii) an ideal model for such an association study. We used a RAD-capture approach to genotype 624 devils at ~16,000 loci and then used association analyses to assess the heritability of three cancer-related phenotypes: infection case-control (where cases were infected devils and controls were devils that were never infected), age of first infection and survival following infection. The SNP array explained much of the phenotypic variance for female survival (>80%) and female case-control (>61%). We found that a few large-effect SNPs explained much of the variance for female survival (~5 SNPs explained >61% of the total variance), whereas more SNPs (~56) of smaller effect explained less of the variance for female case-control (~23% of the total variance). By contrast, these same SNPs did not account for a significant proportion of phenotypic variance in males, suggesting that the genetic bases of these traits and/or selection differ across sexes. Loci involved with cell adhesion and cell-cycle regulation underlay trait variation, suggesting that the devil immune system is rapidly evolving to recognize and potentially suppress cancer growth through these pathways. Overall, our study provided necessary data for genomics-based conservation and management in Tasmanian devils.

Phylogenetics support an ancient common origin of two scientific icons: Devils Hole and Devils Hole pupfish.

The Devils Hole pupfish (Cyprinodon diabolis; DHP) is an icon of conservation biology. Isolated in a 50 m(2) pool (Devils Hole), DHP is one of the rarest vertebrate species known and an evolutionary anomaly, having survived in complete isolation for thousands of years. However, recent findings suggest DHP might be younger than commonly thought, potentially introduced to Devils Hole by humans in the past thousand years. As a result, the significance of DHP from an evolutionary and conservation perspective has been questioned. Here we present a high-resolution genomic analysis of DHP and two closely related species, with the goal of thoroughly examining the temporal divergence of DHP. To this end, we inferred the evolutionary history of DHP from multiple random genomic subsets and evaluated four historical scenarios using the multispecies coalescent. Our results provide substantial information regarding the evolutionary history of DHP. Genomic patterns of secondary contact present strong evidence that DHP were isolated in Devils Hole prior to 20-10 ka and the model best supported by geological history and known mutation rates predicts DHP diverged around 60 ka, approximately the same time Devils Hole opened to the surface. We make the novel prediction that DHP colonized and have survived in Devils Hole since the cavern opened, and the two events (colonization and collapse of the cavern's roof) were caused by a common geologic event. Our results emphasize the power of evolutionary theory as a predictive framework and reaffirm DHP as an important evolutionary novelty, worthy of continued conservation and exploration.

The conserved protein DCN-1/Dcn1p is required for cullin neddylation in C. elegans and S. cerevisiae.

SCF-type E3 ubiquitin ligases are multi-protein complexes required for polyubiquitination and subsequent degradation of target proteins by the 26S proteasome. Cullins, together with the RING-finger protein Rbx1, form the catalytic core of the ligase, and recruit the substrate-recognition module. Cycles of covalent modification of cullins by the ubiquitin-like molecule Nedd8 (neddylation) and removal of Nedd8 by the COP9 signalosome (deneddylation) positively regulate E3 ligase activity. Here we report the identification and analysis of a widely conserved protein that is required for cullin neddylation in the nematode Caenorhabditis elegans and the yeast Saccharomyces cerevisiae. C. elegans DCN-1 and S. cerevisiae Dcn1p (defective in cullin neddylation) are characterized by a novel UBA-like ubiquitin-binding domain and a DUF298 domain of unknown function. Consistent with their requirements for neddylation, DCN-1 and Dcn1p directly bind Nedd8 and physically associate with cullins in both species. Moreover, overexpression of Dcn1p in yeast results in the accumulation of Nedd8-modified cullin Cdc53p. Both in vivo and in vitro experiments indicate that Dcn1p does not inhibit deneddylation of Cdc53p by the COP9 signalosome, but greatly increases the kinetics of the neddylation reaction.

Dynein modifiers in C. elegans: light chains suppress conditional heavy chain mutants.

Cytoplasmic dynein is a microtubule-dependent motor protein that functions in mitotic cells during centrosome separation, metaphase chromosome congression, anaphase spindle elongation, and chromosome segregation. Dynein is also utilized during interphase for vesicle transport and organelle positioning. While numerous cellular processes require cytoplasmic dynein, the mechanisms that target and regulate this microtubule motor remain largely unknown. By screening a conditional Caenorhabditis elegans cytoplasmic dynein heavy chain mutant at a semipermissive temperature with a genome-wide RNA interference library to reduce gene functions, we have isolated and characterized twenty dynein-specific suppressor genes. When reduced in function, these genes suppress dynein mutants but not other conditionally mutant loci, and twelve of the 20 specific suppressors do not exhibit sterile or lethal phenotypes when their function is reduced in wild-type worms. Many of the suppressor proteins, including two dynein light chains, localize to subcellular sites that overlap with those reported by others for the dynein heavy chain. Furthermore, knocking down any one of four putative dynein accessory chains suppresses the conditional heavy chain mutants, suggesting that some accessory chains negatively regulate heavy chain function. We also identified 29 additional genes that, when reduced in function, suppress conditional mutations not only in dynein but also in loci required for unrelated essential processes. In conclusion, we have identified twenty genes that in many cases are not essential themselves but are conserved and when reduced in function can suppress conditionally lethal C. elegans cytoplasmic dynein heavy chain mutants. We conclude that conserved but nonessential genes contribute to dynein function during the essential process of mitosis.

Rapture facilitates inexpensive and high-throughput parent-based tagging in salmonids.

Accurate methods for tracking individuals are crucial to the success of fisheries and aquaculture management. Management of migratory salmonid populations, which are important for the health of many economies, ecosystems, and indigenous cultures, is particularly dependent on data gathered from tagged fish. However, the physical tagging methods currently used have many challenges including cost, variable marker retention, and information limited to tagged individuals. Genetic tracking methods combat many of the problems associated with physical tags, but have their own challenges including high cost, potentially difficult marker design, and incompatibility of markers across species. Here we show the feasibility of a new genotyping method for parent-based tagging (PBT), where individuals are tracked through the inherent genetic relationships with their parents. We found that Rapture sequencing, a combination of restriction-site associated DNA and capture sequencing, provides sufficient data for parentage assignment. Additionally, the same capture bait set, which targets specific restriction-site associated DNA loci, can be used for both Rainbow Trout Oncorhynchus mykiss and Chinook Salmon Oncorhynchus tshawytscha. We input 248 single nucleotide polymorphisms from 1,121 samples to parentage assignment software and compared parent-offspring relationships of the spawning pairs recorded in a hatchery. Interestingly, our results suggest sperm contamination during hatchery spawning occurred in the production of 14% of offspring, further confirming the need for genetic tagging in accurately tracking individuals. PBT with Rapture successfully assigned progeny to parents with a 98.86% accuracy with sufficient genetic data. Cost for this pilot study was approximately $3 USD per sample. As costs vary based on the number of markers used and individuals sequenced, we expect that when implemented at a large-scale, per sample costs could be further decreased. We conclude that Rapture PBT provides a cost-effective and accurate alternative to the physical coded wire tags, and other genetic-based methods.

A versatile Rapture (RAD-Capture) platform for genotyping marine turtles.

Advances in high-throughput sequencing (HTS) technologies coupled with increased interdisciplinary collaboration are rapidly expanding capacity in the scope and scale of wildlife genetic studies. While existing HTS methods can be directly applied to address some evolutionary and ecological questions, certain research goals necessitate tailoring methods to specific study organisms, such as high-throughput genotyping of the same loci that are comparable over large spatial and temporal scales. These needs are particularly common for studies of highly mobile species of conservation concern like marine turtles, where life history traits, limited financial resources and other constraints require affordable, adaptable methods for HTS genotyping to meet a variety of study goals. Here, we present a versatile marine turtle HTS targeted enrichment platform adapted from the recently developed Rapture (RAD-Capture) method specifically designed to meet these research needs. Our results demonstrate consistent enrichment of targeted regions throughout the genome and discovery of candidate variants in all species examined for use in various conservation genetics applications. Accurate species identification confirmed the ability of our platform to genotype over 1,000 multiplexed samples and identified areas for future methodological improvement such as optimization for low initial concentration samples. Finally, analyses within green turtles supported the ability of this platform to identify informative SNPs for stock structure, population assignment and other applications over a broad geographic range of interest to management. This platform provides an additional tool for marine turtle genetic studies and broadens capacity for future large-scale initiatives such as collaborative global marine turtle genetic databases.

The Local South American Chicken Populations Are a Melting-Pot of Genomic Diversity.

Chicken have a considerable impact in South American rural household economy as a source of animal protein (eggs and meat) and a major role in cultural traditions (e.g., cockfighting, religious ceremonies, folklore). A large number of phenotypes and its heterogeneity are due to the multitude of environments (from arid to tropical rain forest and high altitude) and agricultural systems (highly industrialized to subsistence agriculture). This heterogeneity also represents the successive introduction of domestic chicken into this continent, which some consider predating Columbus' arrival to South America. In this study, we have used next-generation restriction site-associated DNA sequencing to scan for genome-wide variation across 145 South American chickens representing local populations from six countries of South America (Colombia, Brazil, Ecuador, Peru, Bolivia, and Chile). After quality control, the genotypes of 122,801 single nucleotide polymorphisms (SNPs) were used to assess the genomic diversity and interpopulation genetic relationship between those populations and their potential sources. The estimated population genetic diversity displayed that the gamefowl has the least diverse population (θπ = 0.86; θS = 0.70). This population is also the most divergent (F ST = 0.11) among the South American populations. The allele-sharing analysis and the admixture analysis revealed that the current diversity displayed by these populations resulted from multiple admixture events with a strong influence of the modern commercial egg-layer chicken (ranging between 44% and 79%). It also revealed an unknown genetic component that is mostly present in the Easter Island population that is also present in local chicken populations from the South American Pacific fringe.

Yeast go the whole HOG for the hyperosmotic response.

An evolutionarily conserved mitogen-activated protein kinase pathway--the high osmolarity glycerol (HOG) pathway--mediates the hyperosmotic response in Saccharomyces cerevisiae. A variety of powerful approaches has generated a comprehensive picture of how cells respond to this stress condition. Several presumptive osmosensors on the cell surface recruit and activate downstream signaling components, which regulate the activity of transcription factors to control gene expression.

A third osmosensing branch in Saccharomyces cerevisiae requires the Msb2 protein and functions in parallel with the Sho1 branch.

Two Saccharomyces cerevisiae plasma membrane-spanning proteins, Sho1 and Sln1, function during increased osmolarity to activate a mitogen-activated protein (MAP) kinase cascade. One of these proteins, Sho1, utilizes the MAP kinase kinase kinase Ste11 to activate Pbs2. We previously used the FUS1 gene of the pheromone response pathway as a reporter to monitor cross talk in hog1 mutants. Cross talk requires the Sho1-Ste11 branch of the HOG pathway, but some residual signaling, which is STE11 dependent, still occurs in the absence of Sho1. These observations led us to propose the existence of another osmosensor upstream of Ste11. To identify such an osmosensor, we screened for mutants in which the residual signaling in a hog1 sho1 mutant was further reduced. We identified the MSB2 gene, which encodes a protein with a single membrane-spanning domain and a large presumptive extracellular domain. Assay of the FUS1-lacZ reporter (in a hog1 mutant background) showed that sho1 and msb2 mutations both reduced the expression of the reporter partially and that the hog1 sho1 msb2 mutant was severely defective in the expression of the reporter. The use of DNA microarrays to monitor gene expression revealed that Sho1 and Msb2 regulate identical gene sets in hog1 mutants. A role for MSB2 in HOG1 strains was also seen in strains defective in the two known branches that activate Pbs2: an ssk1 sho1 msb2 strain was more osmosensitive than an ssk1 sho1 MSB2 strain. These observations indicate that Msb2 is partially redundant with the Sho1 osmosensing branch for the activation of Ste11.

Unique and redundant roles for HOG MAPK pathway components as revealed by whole-genome expression analysis.

The Saccharomyces cerevisiae high osmolarity glycerol (HOG) mitogen-activated protein kinase pathway is required for osmoadaptation and contains two branches that activate a mitogen-activated protein kinase (Hog1) via a mitogen-activated protein kinase kinase (Pbs2). We have characterized the roles of common pathway components (Hog1 and Pbs2) and components in the two upstream branches (Ste11, Sho1, and Ssk1) in response to elevated osmolarity by using whole-genome expression profiling. Several new features of the HOG pathway were revealed. First, Hog1 functions during gene induction and repression, cross talk inhibition, and in governing the regulatory period. Second, the phenotypes of pbs2 and hog1 mutants are identical, indicating that the sole role of Pbs2 is to activate Hog1. Third, the existence of genes whose induction is dependent on Hog1 and Pbs2 but not on Ste11 and Ssk1 suggests that there are additional inputs into Pbs2 under our inducing conditions. Fourth, the two upstream pathway branches are not redundant: the Sln1-Ssk1 branch has a much more prominent role than the Sho1-Ste11 branch for activation of Pbs2 by modest osmolarity. Finally, the general stress response pathway and both branches of the HOG pathway all function at high osmolarity. These studies demonstrate that cells respond to increased osmolarity by using different signal transduction machinery under different conditions.

Genomics clarifies taxonomic boundaries in a difficult species complex.

Efforts to taxonomically delineate species are often confounded with conflicting information and subjective interpretation. Advances in genomic methods have resulted in a new approach to taxonomic identification that stands to greatly reduce much of this conflict. This approach is ideal for species complexes, where divergence times are recent (evolutionarily) and lineages less well defined. The California Roach/Hitch fish species complex is an excellent example, experiencing a convoluted geologic history, diverse habitats, conflicting species designations and potential admixture between species. Here we use this fish complex to illustrate how genomics can be used to better clarify and assign taxonomic categories. We performed restriction-site associated DNA (RAD) sequencing on 255 Roach and Hitch samples collected throughout California to discover and genotype thousands of single nucleotide polymorphism (SNPs). Data were then used in hierarchical principal component, admixture, and FST analyses to provide results that consistently resolved a number of ambiguities and provided novel insights across a range of taxonomic levels. At the highest level, our results show that the CA Roach/Hitch complex should be considered five species split into two genera (4 + 1) as opposed to two species from distinct genera (1 +1). Subsequent levels revealed multiple subspecies and distinct population segments within identified species. At the lowest level, our results indicate Roach from a large coastal river are not native but instead introduced from a nearby river. Overall, this study provides a clear demonstration of the power of genomic methods for informing taxonomy and serves as a model for future studies wishing to decipher difficult species questions. By allowing for systematic identification across multiple scales, taxonomic structure can then be tied to historical and contemporary ecological, geographic or anthropogenic factors.

The evolutionary basis of premature migration in Pacific salmon highlights the utility of genomics for informing conservation.

The delineation of conservation units (CUs) is a challenging issue that has profound implications for minimizing the loss of biodiversity and ecosystem services. CU delineation typically seeks to prioritize evolutionary significance, and genetic methods play a pivotal role in the delineation process by quantifying overall differentiation between populations. Although CUs that primarily reflect overall genetic differentiation do protect adaptive differences between distant populations, they do not necessarily protect adaptive variation within highly connected populations. Advances in genomic methodology facilitate the characterization of adaptive genetic variation, but the potential utility of this information for CU delineation is unclear. We use genomic methods to investigate the evolutionary basis of premature migration in Pacific salmon, a complex behavioral and physiological phenotype that exists within highly connected populations and has experienced severe declines. Strikingly, we find that premature migration is associated with the same single locus across multiple populations in each of two different species. Patterns of variation at this locus suggest that the premature migration alleles arose from a single evolutionary event within each species and were subsequently spread to distant populations through straying and positive selection. Our results reveal that complex adaptive variation can depend on rare mutational events at a single locus, demonstrate that CUs reflecting overall genetic differentiation can fail to protect evolutionarily significant variation that has substantial ecological and societal benefits, and suggest that a supplemental framework for protecting specific adaptive variation will sometimes be necessary to prevent the loss of significant biodiversity and ecosystem services.

Genome Report: Identification and Validation of Antigenic Proteins from Pajaroellobacter abortibovis Using De Novo Genome Sequence Assembly and Reverse Vaccinology.

Epizootic bovine abortion (EBA), or "foothill abortion," is the leading cause of beef cattle abortion in California and has also been reported in Nevada and Oregon. In the 1970s, the soft-shelled tick Ornithodoros coriaceus, or "pajaroello tick," was confirmed as the disease-transmitting vector. In 2005, a novel Deltaproteobacterium was discovered as the etiologic agent of EBA (aoEBA), recently named Pajaroellobacter abortibovis This organism cannot be grown in culture using traditional microbiological techniques; it can only be grown in experimentally-infected severe combined immunodeficient (SCID) mice. The objectives of this study were to perform a de novo genome assembly for P. abortibovis and identify and validate potential antigenic proteins as candidates for future recombinant vaccine development. DNA and RNA were extracted from spleen tissue collected from experimentally-infected SCID mice following exposure to P. abortibovis This combination of mouse and bacterial DNA was sequenced and aligned to the mouse genome. Mouse sequences were subtracted from the sequence pool and the remaining sequences were de novo assembled at 50x coverage into a 1.82 Mbp complete closed circular Deltaproteobacterial genome containing 2250 putative protein-coding sequences. Phylogenetic analysis of P. abortibovis predicts that this bacterium is most closely related to the organisms of the order Myxococcales, referred to as Myxobacteria. In silico prediction of vaccine candidates was performed using a reverse vaccinology approach resulting in the identification and ranking of the top 10 candidate proteins that are likely to be antigenic. Immunologic testing of these candidate proteins confirmed antigenicity of seven of the nine expressed protein candidates using serum from P. abortibovis immunized mice.

Whole-Genome Sequencing and Concordance Between Antimicrobial Susceptibility Genotypes and Phenotypes of Bacterial Isolates Associated with Bovine Respiratory Disease.

Extended laboratory culture and antimicrobial susceptibility testing timelines hinder rapid species identification and susceptibility profiling of bacterial pathogens associated with bovine respiratory disease, the most prevalent cause of cattle mortality in the United States. Whole-genome sequencing offers a culture-independent alternative to current bacterial identification methods, but requires a library of bacterial reference genomes for comparison. To contribute new bacterial genome assemblies and evaluate genetic diversity and variation in antimicrobial resistance genotypes, whole-genome sequencing was performed on bovine respiratory disease-associated bacterial isolates (Histophilus somni, Mycoplasma bovis, Mannheimia haemolytica, and Pasteurella multocida) from dairy and beef cattle. One hundred genomically distinct assemblies were added to the NCBI database, doubling the available genomic sequences for these four species. Computer-based methods identified 11 predicted antimicrobial resistance genes in three species, with none being detected in M. bovis While computer-based analysis can identify antibiotic resistance genes within whole-genome sequences (genotype), it may not predict the actual antimicrobial resistance observed in a living organism (phenotype). Antimicrobial susceptibility testing on 64 H. somni, M. haemolytica, and P. multocida isolates had an overall concordance rate between genotype and phenotypic resistance to the associated class of antimicrobials of 72.7% (P < 0.001), showing substantial discordance. Concordance rates varied greatly among different antimicrobial, antibiotic resistance gene, and bacterial species combinations. This suggests that antimicrobial susceptibility phenotypes are needed to complement genomically predicted antibiotic resistance gene genotypes to better understand how the presence of antibiotic resistance genes within a given bacterial species could potentially impact optimal bovine respiratory disease treatment and morbidity/mortality outcomes.