A Conservation Hatchery Population of Delta Smelt Shows Evidence of Genetic Adaptation to Captivity After 9 Generations.

Genetic adaptation to captivity is a concern for threatened and endangered species held in conservation hatcheries. Here, we present evidence of genetic adaptation to captivity in a conservation hatchery for the endangered delta smelt (Fish Conservation and Culture Laboratory, University of California Davis; FCCL). The FCCL population is genetically managed with parentage analysis and the addition of wild fish each year. Molecular monitoring indicates little loss of genetic variation and low differentiation between the wild and conservation populations. Yet, we found an increase in offspring survival to reproductive maturity during the subsequent spawning season (recovery rate) in crosses that included one or both cultured parents. Crosses with higher levels of hatchery ancestry tend to produce a greater number of offspring that are recovered the following year. The recovery rate of a cross decreases when offspring are raised in a tank with fish of high levels of hatchery ancestry. We suggest changes in fish rearing practices at the FCCL to reduce genetic adaptation to captivity, as delta smelt numbers in the wild continue to decline and the use of FCCL fish for reintroduction becomes more likely.

Migration-related phenotypic divergence is associated with epigenetic modifications in rainbow trout.

Migration is essential for the reproduction and survival of many animals, yet little is understood about its underlying molecular mechanisms. We used the salmonid Oncorhynchus mykiss to gain mechanistic insight into smoltification, which is a morphological, physiological and behavioural transition undertaken by juveniles in preparation for seaward migration. O. mykiss is experimentally tractable and displays intra- and interpopulation variation in migration propensity. Migratory individuals can produce nonmigratory progeny and vice versa, indicating a high degree of phenotypic plasticity. One potential way that phenotypic plasticity might be linked to variation in migration-related life history tactics is through epigenetic regulation of gene expression. To explore this, we quantitatively measured genome-scale DNA methylation in fin tissue using reduced representation bisulphite sequencing of F2 siblings produced from a cross between steelhead (migratory) and rainbow trout (nonmigratory) lines. We identified 57 differentially methylated regions (DMRs) between smolt and resident O. mykiss juveniles. DMRs were high in magnitude, with up to 62% differential methylation between life history types, and over half of the gene-associated DMRs were in transcriptional regulatory regions. Many of the DMRs encode proteins with activity relevant to migration-related transitions (e.g. circadian rhythm pathway, nervous system development, protein kinase activity). This study provides the first evidence of a relationship between epigenetic variation and life history divergence associated with migration-related traits in any species.

Evaluating the performance of captive breeding techniques for conservation hatcheries: a case study of the delta smelt captive breeding program.

The delta smelt, an endangered fish species endemic to the San Francisco Bay-Delta, California, United States, was recently brought into captivity for species preservation. This study retrospectively evaluates the implementation of a genetic management plan for the captive delta smelt population. The captive genetic management plan entails tagging fish, molecular data collection, pedigree reconstruction, relatedness estimation, and recommending fish crosses annually in an effort to minimize the average coancestry in the population and limit inbreeding. We employed 12 microsatellite DNA markers to examine temporal genetic diversity in consecutive, discrete generations to determine the effects of intensive genetic management on the population and to quantify the amount of wild genetic diversity present within each captive generation. Wild fish are incorporated into the captive population each generation to minimize genetic drift, and 91% of the original founders are still represented in the F(3) generation. The average mean kinship in the third generation in captivity was 0.0035. There was no evidence of significant genetic divergence of the captive population from the wild population. The results of this study yield management insights into the practical application of genetic management plans for captive populations and conservation hatcheries, in an attempt to preserve the genetic integrity of endangered species.

Detecting drug diversion in health-system data using machine learning and advanced analytics.

PURPOSE: The theft of drugs from healthcare facilities, also known as drug diversion, occurs frequently but is often undetected. This paper describes a research study to develop and test novel drug diversion detection methods. Improved diversion detection and reduction in diversion improves patient safety, limits harm to the person diverting, reduces the public health impact of substance use disorder, and mitigates significant liability risk to pharmacists and their organizations. METHODS: Ten acute care inpatient hospitals across 4 independent health systems extracted 2 datasets from various health information technology systems. Both datasets were consolidated, normalized, classified, and sampled to provide a harmonious dataset for analysis. Supervised machine learning methods were iteratively used on the initial sample dataset to train algorithms to classify medication movement transactions as involving a low or high risk of diversion. Thereafter, the resulting machine learning model classified the risk of diversion in a historical dataset capturing 8 to 24 months of history that included 27.9 million medication movement transactions by 19,037 nursing, 1,047 pharmacy, and 712 anesthesia clinicians and that included 22 known, blinded diversion cases to measure when the model would have detected the diversion compared to when the diversion was actually detected by existing methods. RESULTS: The machine learning model had 96.3% accuracy, 95.9% specificity, and 96.6% sensitivity in detecting transactions involving a high risk of diversion using the initial sample dataset. In subsequent testing using the much larger historical dataset, the analytics detected known diversion cases (n = 22) in blinded data faster than existing detection methods (a mean of 160 days and a median of 74 days faster; range, 7-579 days faster). CONCLUSION: The study showed that (1) consolidated datasets and (2) supervised machine learning can detect known diversion cases faster than existing detection methods. Users of the technology also noted improved investigation efficiency.

Discovery of genes implicated in whirling disease infection and resistance in rainbow trout using genome-wide expression profiling.

BACKGROUND: Whirling disease, caused by the pathogen Myxobolus cerebralis, afflicts several salmonid species. Rainbow trout are particularly susceptible and may suffer high mortality rates. The disease is persistent and spreading in hatcheries and natural waters of several countries, including the U.S.A., and the economic losses attributed to whirling disease are substantial. In this study, genome-wide expression profiling using cDNA microarrays was conducted for resistant Hofer and susceptible Trout Lodge rainbow trout strains following pathogen exposure with the primary objective of identifying specific genes implicated in whirling disease resistance. RESULTS: Several genes were significantly up-regulated in skin following pathogen exposure for both the resistant and susceptible rainbow trout strains. For both strains, response to infection appears to be linked with the interferon system. Expression profiles for three genes identified with microarrays were confirmed with qRT-PCR. Ubiquitin-like protein 1 was up-regulated over 100 fold and interferon regulating factor 1 was up-regulated over 15 fold following pathogen exposure for both strains. Expression of metallothionein B, which has known roles in inflammation and immune response, was up-regulated over 5 fold in the resistant Hofer strain but was unchanged in the susceptible Trout Lodge strain following pathogen exposure. CONCLUSION: The present study has provided an initial view into the genetic basis underlying immune response and resistance of rainbow trout to the whirling disease parasite. The identified genes have allowed us to gain insight into the molecular mechanisms implicated in salmonid immune response and resistance to whirling disease infection.

Characterization of a chromosomal rearrangement responsible for producing "apparent" XY-female fall-run Chinook salmon in California.

Fluorescence in situ hybridization (FISH) was used to identify the X and Y chromosomes of offspring produced by normal and "apparent" XY-female fall-run Chinook salmon (Oncorhynchus tshawytscha) from California. FISH experiments were performed using probes to 2 sex-linked loci, growth hormone pseudogene (GH-Psi), and OtY1, as well as a probe to a sex-linked microsatellite (Omy7INRA). Comparison of FISH staining patterns between the offspring produced by normal and apparent XY-females revealed that the apparent XY-female examined transmitted a "Y-like" chromosome with an attenuated OtY1 and GH-Psi signal to half of its offspring. Segregation analysis of microsatellites derived from rainbow trout (Oncorhynchus mykiss) with respect to phenotypic sex was carried out for 2 normal and 2 apparent XY-female Chinook salmon families. Inheritance patterns of Omy7INRA were consistent with this locus being closely linked to GH-Psi in males and in apparent XY-females carrying the Y-like chromosome. Another microsatellite locus (Omm1077) was closely linked to the primary sex-determining locus (SEX) in males but not to GH-Psi/OtY1 in apparent XY-females. The FISH analyses suggest that apparent XY-female fall-run Chinook salmon in California are not the product of a Y chromosome to autosome translocation. Despite the combined FISH and inheritance analyses, we were unable to differentiate between 2 alternative explanations for apparent XY-females, namely, recombination of markers between the sex chromosomes, or a Y chromosome with a dysfunctional or missing sex-determining region.

An induced mass spawn of the hermaphroditic lion-paw scallop, Nodipecten subnodosus: genetic assignment of maternal and paternal parentage.

The Pacific lion-paw scallop is commonly propagated for aquaculture by induced mass spawns of few individuals. Parentage of a mass spawn of this species has not been evaluated nor has the maternal and paternal contribution of each of these functional hermaphrodites to the progeny. Genotypes of 6 spawners and 374 resulting progeny at 6 microsatellite loci were coupled with mitochondrial DNA sequencing to assign maternal and paternal parentage. After the identification of a high proportion of null alleles (9.7%), microsatellite data revealed that 51.7% of the progenies were full siblings, with a significant, unequal contribution of the 6 spawners to the progeny. Three progenies were the result of self-fertilization. All spawners contributed paternally (though unequally); however, 2 spawners were the maternal parents of all but 7 progenies resulting in a variance effective population size of 3.52. DNA sequencing confirmed 4 microsatellite mutations within 4476 alleles scored, all in the paternal germ line. With minor exception, the loci conformed to Mendelian rules of segregation when null alleles were accounted for, and 2 loci were found to be linked. These results lend insight to the genetic composition of induced mass spawns and provide a basis for the development of more effective spawning techniques.

Range-wide phylogeographic structure of the vernal pool fairy shrimp (Branchinecta lynchi).

Wetland habitats across the world are experiencing rapid modification and loss due to accelerating habitat conversion. Impacts to wetland habitats are particularly acute in California where up to 90% of wetland habitats have been modified or lost. Vernal pool ecosystems have therefore undergone a dramatic loss in habitat and along with them an entire endemic fauna is under threat of extinction. Recent efforts to conserve vernal pool habitat and associated species have involved restoration and creation of vernal pools as well as translocations of threatened species. The vernal pool fairy shrimp, Branchinecta lynchi, is one of several endemic and federally listed species being targeted for translocations. To guide reintroduction and conservation, detailed information on range-wide population structure and diversity is needed. We collected genetic data from two mitochondrial genes throughout the known extant range of B. lynchi to elucidate population structure and diversity of the species. We found support for phylogeographic structure throughout the range of B. lynch associated with isolated watersheds and vernal pool regions previously identified in the recovery plan for the species. The underlying mechanisms responsible for this broad pattern of genetic structure have yet to be identified. However, the evidence of only a few haplotypes being shared across the species range and patterns of isolation by distance within vernal pool regions suggests dispersal limitation may play a role. These results stress that conservation programs, at a minimum, should consider using individuals from regional populations as sources for reintroductions to maintain historical patterns of genetic differentiation. Additionally, because genetic structure is associated with vernal pool regions which are based on local hydrology and geology, translocations should proceed considering the distance between donor and recipient sites.

Sequencing improves our ability to study threatened migratory species: Genetic population assignment in California's Central Valley Chinook salmon.

Effective conservation and management of migratory species requires accurate identification of unique populations, even as they mix along their migratory corridors. While telemetry has historically been used to study migratory animal movement and habitat use patterns, genomic tools are emerging as a superior alternative in many ways, allowing large-scale application at reduced costs. Here, we demonstrate the usefulness of genomic resources for identifying single-nucleotide polymorphisms (SNPs) that allow fast and accurate identification of the imperiled Chinook salmon in the Great Central Valley of California. We show that 80 well-chosen loci, drawn from a pool of over 11,500 SNPs developed from restriction site-associated DNA sequencing, can accurately identify Chinook salmon runs and select populations within run. No other SNP panel for Central Valley Chinook salmon has been able to achieve the high accuracy of assignment we show here. This panel will greatly improve our ability to study and manage this ecologically, economically, and socially important species and demonstrates the great utility of using genomics to study migratory species.

Genetic structure of the tick Ornithodoros coriaceus (Acari: Argasidae) in California, Nevada, and Oregon.

The argasid tick Ornithodoros coriaceus (Koch) is the only confirmed vector of epizootic bovine abortion (EBA) in the United States. The disease and its tick vector have historically been reported in the foothills of the Sierra Nevada Mountains and coast ranges of California. In the past two decades, the range of EBA has apparently expanded into southern Oregon and northern Nevada. Possible explanations for this expansion include 1) increased recognition and reporting of EBA in these regions; 2) widespread movement of tick-infested and EBA-infected hosts with subsequent colonization of these regions by infected ticks; and 3) widespread movement of the EBA agent, independent of tick movements, into extant tick populations in these new regions. The current study was performed to evaluate these hypotheses by examining patterns of variability in a 420-bp segment of the 16S mitochondrial rDNA gene sequence among 210 O. coriaceus individuals from 14 sites in California, Oregon, and Nevada. Sixty-three unique haplotypes were identified in the ticks tested, with 84% of the sequence variation attributable to among-population variation and 16% to within-population variation. A majority of the haplotypes were unique to their particular collection site, whereas only four collection sites shared haplotypes. Overall, very little evidence of gene flow among tick populations was detected, making it unlikely that widespread tick movement had introduced O. coriaceus and the EBA agent into new regions.

Ten real-time PCR assays for detection of fish predation at the community level in the San Francisco Estuary-Delta.

The effect of predation on native fish by introduced species in the San Francisco Estuary-Delta (SFE) has not been thoroughly studied despite its potential to impact species abundances. Species-specific quantitative PCR (qPCR) is an accurate method for identifying species from exogenous DNA samples. Quantitative PCR assays can be used for detecting prey in gut contents or faeces, discriminating between cryptic species, or detecting rare aquatic species. We designed ten TaqMan qPCR assays for fish species from the SFE watershed most likely to be affected by non-native piscivores. The assays designed are highly specific, producing no signal from co-occurring or related species, and sensitive, with a limit of detection between 3.2 and 0.013 pg/µL of target DNA. These assays will be used in conjunction with a high-throughput qPCR platform to compare predation rates between native and non-native piscivores and assess the impacts of predation in the system.

Transcriptional Response to Acute Thermal Exposure in Juvenile Chinook Salmon Determined by RNAseq.

Thermal exposure is a serious and growing challenge facing fish species worldwide. Chinook salmon (Oncorhynchus tshawytscha) living in the southern portion of their native range are particularly likely to encounter warmer water due to a confluence of factors. River alterations have increased the likelihood that juveniles will be exposed to warm water temperatures during their freshwater life stage, which can negatively impact survival, growth, and development and pose a threat to dwindling salmon populations. To better understand how acute thermal exposure affects the biology of salmon, we performed a transcriptional analysis of gill tissue from Chinook salmon juveniles reared at 12° and exposed acutely to water temperatures ranging from ideal to potentially lethal (12° to 25°). Reverse-transcribed RNA libraries were sequenced on the Illumina HiSeq2000 platform and a de novo reference transcriptome was created. Differentially expressed transcripts were annotated using Blast2GO and relevant gene clusters were identified. In addition to a high degree of downregulation of a wide range of genes, we found upregulation of genes involved in protein folding/rescue, protein degradation, cell death, oxidative stress, metabolism, inflammation/immunity, transcription/translation, ion transport, cell cycle/growth, cell signaling, cellular trafficking, and structure/cytoskeleton. These results demonstrate the complex multi-modal cellular response to thermal stress in juvenile salmon.

Genetic structure of mountain lion (Puma concolor) populations in California.

Analysis of 12 microsatellite loci from431 mountain lions (Puma concolor)revealed distinct genetic subdivision that wasassociated with geographic barriers andisolation by distance in California. Levels ofgenetic variation differed among geographicregions, and mountain lions that inhabitedcoastal areas exhibited less heterozygositythan those sampled inland. The San FranciscoBay and Sacramento-San Joaquin River Delta, theCentral Valley, and the Los Angeles Basinappeared to be substantial barriers to geneflow, and allele frequencies of populationsseparated by those features differedsubstantially. A partial barrier to gene flowappeared to exist along the crest of the SierraNevada. Estimated gene flow was high amongmountain lions inhabiting the Modoc Plateau,the western Sierra Nevada, and northern sectionof the eastern Sierra Nevada. SouthernCalifornia mountain lion populations mayfunction as a metapopulation; however, humandevelopments threaten to eliminate habitat andmovement corridors. While north-south geneflow along the western Sierra Nevada wasestimated to be very high, projected loss andfragmentation of foothill habitat may reducegene flow and subdivide populations. Preservation of existing movement corridorsamong regions could prevent population declinesand loss of genetic variation. This studyshows that mountain lion management andconservation efforts should be individualizedaccording to region and incorporatelandscape-level considerations to protecthabitat connectivity.

Increased accuracy of species lists developed for alpine lakes using morphology and cytochrome oxidase I for identification of specimens.

The first step in many community ecology studies is to produce a species list from a sample of individuals. Community ecologists now have two viable ways of producing a species list: morphological and barcode identification. In this study, we compared the taxonomic resolution gained by a combined use of both methods and tested whether a change in taxonomic resolution significantly impacted richness estimates for benthic macroinvertebrates sampled from ten lakes in Sequoia National Park, USA. Across all lakes, 77 unique taxa were identified and 42% (32) were reliably identified to species using both barcode and morphological identification. Of the 32 identified to species, 63% (20) were identified solely by comparing the barcode sequence from cytochrome oxidase I to the Barcode of Life reference library. The increased resolution using a combined identification approach compared to identifications based solely on morphology resulted in a significant increase in estimated richness within a lake at the order, family, genus and species levels of taxonomy (P < 0.05). Additionally, young or damaged individuals that could not be identified using morphology were identified using their COI sequences to the genus or species level on average 75% of the time. Our results demonstrate that a combined identification approach improves accuracy of benthic macroinvertebrate species lists in alpine lakes and subsequent estimates of richness. We encourage the use of barcodes for identification purposes and specifically when morphology is insufficient, as in the case of damaged and early life stage specimens of benthic macroinvertebrates.

Climate change likely to facilitate the invasion of the non-native hydroid, Cordylophora caspia, in the San Francisco Estuary.

Climate change and invasive species can both have negative impacts on native species diversity. Additionally, climate change has the potential to favor invasive species over natives, dealing a double blow to native biodiversity. It is, therefore, vital to determine how changing climate conditions are directly linked to demographic rates and population growth of non-native species so we can quantitatively evaluate how invasive populations may be affected by changing conditions and, in turn, impact native species. Cordylophora caspia, a hydrozoan from the Ponto-Caspian region, has become established in the brackish water habitats of the San Francisco Estuary (SFE). We conducted laboratory experiments to study how temperature and salinity affect C. caspia population growth rates, in order to predict possible responses to climate change. C. Caspia population growth increased nonlinearly with temperature and leveled off at a maximum growth rate near the annual maximum temperature predicted under a conservative climate change scenario. Increasing salinity, however, did not influence growth rates. Our results indicate that C. caspia populations in the SFE will benefit from predicted regional warming trends and be little affected by changes in salinity. The population of C. caspia in the SFE has the potential to thrive under future climate conditions and may subsequently increase its negative impact on the food web.

TaqMan assays for the genetic identification of delta smelt (Hypomesus transpacificus) and wakasagi smelt (Hypomesus nipponensis).

We have developed species-specific TaqMan assays for two California fish species, the threatened delta smelt (Hypomesus transpacificus) and the introduced wakasagi smelt (Hypomesus nipponensis). The assays are capable of correctly identifying each species with 100% accuracy, with no cross-species amplification. We anticipate these assays will prove useful for future scientific studies requiring genetic species identification (e.g. predation of smelt) or monitoring (e.g. detection of delta smelt near water diversions).

Nuclear and mtDNA lineage diversity in wild and cultured Pacific lion-paw scallop, Nodipecten subnodosus (Baja California Peninsula, Mexico).

Pacific lion-paw scallops were collected from natural aggregations in Laguna Ojo de Liebre (Pacific Ocean), the Gulf of California, and from aquaculture facilities for genetic diversity analyses. Mitochondrial DNA sequencing uncovered two highly supported clades separated by 2.5% divergence. Data from ten microsatellite markers suggest individuals from these mitogroups are introgressed, raising questions about the mitotype origin. Some evidence suggests gene flow between La Paz and Ojo de Liebre; otherwise the Gulf of California and Ojo de Liebre are acting as two distinct populations. It is unclear whether translocations between sites have influenced the observed genetic structure or whether gene flow has been facilitated by past geologic events. Finally, scallops spawned for aquaculture are unique from the wild and have significantly less diversity. These results warrant the attention of managers and producers who should work to monitor and conserve genetic diversity in both wild and aquaculture populations.

Subspecies-informative SNP assays for evaluating introgression between native golden trout and introduced rainbow trout.

We characterize 20 single nucleotide polymorphism assays for evaluating hybridization between native golden trout subspecies (Oncorhynchus mykiss aguabonita and O. m. whitei) and introduced rainbow trout strains. These assays utilize the 5'-nuclease reaction, facilitating high-throughput genotyping of many individuals and making them useful in quantifying and monitoring introgression and potentially applicable to studies of other O. mykiss groups. Minor allele frequency differentials (δq) among native and introduced rainbow groups ranged from 0 to 1, with an average differential of 0.75 for both subspecies aguabonita and whitei relative to the hatchery rainbow trout strain.

Six diagnostic single nucleotide polymorphism markers for detecting introgression between cutthroat and rainbow trouts.

Ten primer pairs were screened to develop single nucleotide polymorphism (SNP) TaqMan assays that will distinguish California golden trout and some rainbow trouts (Oncorhynchus mykiss sspp., O. m. aguabonita) from the Paiute and Lahontan cutthroat trouts (Oncorhynchus clarkii seleniris, O. c. henshawi). From these 10 primer pairs, one mitochondrial and five nuclear fixed SNP differences were discovered and developed into TaqMan assays. These six assays will be useful for characterizing and monitoring hybridization between these groups. Additional Oncorhynchus clarkii sspp. and Oncorhynchus mykiss sspp. were assayed to determine if these assays are useful in closely related species.