Single generation epigenetic change in captivity and reinforcement in subsequent generations in a delta smelt (Hypomesus transpacificus) conservation hatchery.

A refugial population of the endangered delta smelt (Hypomesus transpacificus) has been maintained at the Fish Conservation and Culture Laboratory (FCCL) at UC Davis since 2008. Despite intense genetic management, fitness differences between wild and cultured fish have been observed at the FCCL. To investigate the molecular underpinnings of hatchery domestication, we used whole-genome bisulfite sequencing to quantify epigenetic differences between wild and hatchery-origin delta smelt. Differentially methylated regions (DMRs) were identified from 104 individuals by comparing the methylation patterns in different generations of hatchery fish (G1, G2, G3) with their wild parents (G0). We discovered a total of 132 significant DMRs (p < .05) between G0 and G1, 132 significant DMRs between G0 and G2, and 201 significant DMRs between G0 and G3. Our results demonstrate substantial differences in methylation patterns emerged between the wild and hatchery-reared fish in the early generations in the hatchery, with a higher proportion of hypermethylated DMRs in hatchery-reared fish. The rearing environment was found to be a stronger predictor of individual clustering based on methylation patterns than family, sex or generation. Our study indicates a reinforcement of the epigenetic status with successive generations in the hatchery environment, as evidenced by an increase in methylation in hypermethylated DMRs and a decrease in methylation in hypomethylated DMRs over time. Lastly, our results demonstrated heterogeneity in inherited methylation pattern in families across generations. These insights highlight the long-term consequences of hatchery practices on the epigenetic landscape, potentially impacting wild fish populations.

Evaluating environmental DNA detection of a rare fish in turbid water using field and experimental approaches.

Detection sensitivity of aquatic species using environmental DNA (eDNA) generally decreases in turbid water but is poorly characterized. In this study, eDNA detection targeted delta smelt (Hypomesus transpacificus), a critically endangered estuarine fish associated with turbid water. eDNA sampling in the field was first paired with a trawl survey. Species-specific detection using a Taqman qPCR assay showed concordance between the methods, but a weak eDNA signal. Informed by the results of field sampling, an experiment was designed to assess how turbidity and filtration methods influence detection of a rare target. Water from non-turbid (5 NTU) and turbid (50 NTU) estuarine sites was spiked with small volumes (0.5 and 1 mL) of water from a delta smelt tank to generate low eDNA concentrations. Samples were filtered using four filter types: cartridge filters (pore size 0.45 µm) and 47 mm filters (glass fiber, pore size 1.6 µm and polycarbonate, pore sizes 5 and 10 µm). Prefiltration was also tested as an addition to the filtration protocol for turbid water samples. eDNA copy numbers were analyzed using a censored data method for qPCR data. The assay limits and lack of PCR inhibition indicated an optimized assay. Glass fiber filters yielded the highest detection rates and eDNA copies in non-turbid and turbid water. Prefiltration improved detection in turbid water only when used with cartridge and polycarbonate filters. Statistical analysis identified turbidity as a significant effect on detection probability and eDNA copies detected; filter type and an interaction between filter type and prefilter were significant effects on eDNA copies detected, suggesting that particulate-filter interactions can affect detection sensitivity. Pilot experiments and transparent criteria for positive detection could improve eDNA surveys of rare species in turbid environments.

Loss of plasticity in maturation timing after ten years of captive spawning in a delta smelt conservation hatchery.

Adaptation to captivity in spawning programs can lead to unintentional consequences, such as domestication that results in reduced fitness in the wild. The timing of sexual maturation has been shown to be a trait under domestication selection in fish hatcheries, which affects a fish's access to mating opportunities and aligning their offspring's development with favorable environmental conditions. Earlier maturing fish may be favored in hatchery settings where managers provide artificially optimal growing conditions, but early maturation may reduce fitness in the wild if, for example, there is a mismatch between timing of reproduction and availability of resources that support recruitment. We investigated patterns of maturation timing in a delta smelt (Hypomesus transpacificus) conservation hatchery by quantifying changes to the median age at maturity since the captive spawning program was initiated in 2008. Over the span of a decade, we observed a small, but significant increase in age at maturity among broodstock by 2.2 weeks. This trait had low heritability and was largely controlled by phenotypic plasticity that was dependent on the time of year fish were born. Fish that were born later in the year matured faster, potentially a carryover from selection favoring synchronous spawning in the wild. However, higher DI (domestication index) fish showed a loss of plasticity, we argue, as a result of hatchery practices that breed individuals past peak periods of female ripeness. Our findings suggest that the hatchery setting has relaxed selection pressures for fish to mature quickly at the end of the year and, consequently, has led to a loss of plasticity in age at maturity. Hatchery fish that are re-introduced in the wild may not be able to align maturation with population peaks if their maturation rates are too slow with reduced plasticity, potentially resulting in lower fitness.

Captive-reared Delta Smelt (Hypomesus transpacificus) exhibit high survival in natural conditions using in situ enclosures.

Conservation of endangered fishes commonly includes captive breeding, applied research, and management. Since 1996, a captive breeding program has existed for the federally threatened and California endangered Delta Smelt Hypomesus transpacificus, an osmerid fish endemic to the upper San Francisco Estuary. Although this program serves as a captive refuge population, with experimental releases being initiated to supplement the wild population, it was uncertain how individuals would survive, feed, and maintain condition outside hatchery conditions. We evaluated this and the effects of three enclosure designs (41% open, 63% open, and 63% open with partial outer mesh wrap) on growth, survival, and feeding efficacy of cultured Delta Smelt at two locations (Sacramento River near Rio Vista, CA and in Sacramento River Deepwater Ship Channel) in the wild. Enclosures exposed fish to semi-natural conditions (ambient environmental fluctuations and wild food resources) but prevented escape and predation. After four weeks, survival was high for all enclosure types (94-100%) at both locations. The change in condition and weight was variable between sites, increasing at the first location but decreasing at the second location. Gut content analysis showed that fish consumed wild zooplankton that came into the enclosures. Cumulatively, results show that captive-reared Delta Smelt can survive and forage successfully when housed in enclosures under semi-natural conditions in the wild. When comparing enclosure types, we observed no significant difference in fish weight changes (p = 0.58-0.81 across sites). The success of housing captive-reared Delta Smelt in enclosures in the wild provides preliminary evidence that these fish may be suitable to supplement the wild population in the San Francisco Estuary. Furthermore, these enclosures are a new tool to test the efficacy of habitat management actions or to acclimate fish to wild conditions as a soft release strategy for recently initiated supplementation efforts.

Molecular Systematics of Redband Trout from Genome-Wide DNA Sequencing Substantiates the Description of a New Taxon (Salmonidae: Oncorhynchus mykiss calisulat) from the McCloud River.

Rainbow Trout, Oncorhynchus mykiss, exhibit high levels of phenotypic diversity leading to the recognition of numerous subspecies. A major distinction among Rainbow Trout subspecies exists between Coastal Rainbow Trout (O. m. irideus), which occurs west of the Cascade and Sierra Nevada mountain ranges, and interior Redband Trout (O. mykiss sspp.), largely distributed to the east. Interior Redband Trout are composed of three primary lineages and can share various outward, anatomical or physiological characteristics that are often symplesiomorphies or examples of convergence. We examine high-throughput DNA sequence data from Sacramento Redband Trout O. m. stonei from the Upper Pit and Upper McCloud Rivers along with representatives of Rainbow Trout and Golden Trout lineages to clarify the composition and relationships of the Sacramento Redband Trout. We find O. m. stonei to be polyphyletic, divided between populations in the Pit River and the Upper McCloud River. Redband Trout obtained from the Pit River are most-closely related to Great Basin Redband Trout O. m. newberrii and to fish of the Warner Lakes Basin and Surprise Valley within the O. m. newberrii lineage. The type specimen of O. m. stonei, collected from the Lower McCloud River, is phenotypically similar to Great Basin Redband Trout. We find as well that the isolated populations of trout in the Upper McCloud River Basin represent a lineage of Rainbow Trout now restricted to that region, are monophyletic and are not most closely related to the interior Redband Trout lineages of O. m. newberrii or O. m. gairdnerii. Furthermore, they are not represented by the type specimens of O. m. stonei or O. m. shasta. Consequently, we formally describe the McCloud River Redband Trout O. mykiss calisulat, new subspecies.

A draft reference genome of the Vernal Pool Fairy Shrimp, Branchinecta lynchi.

We present the reference genome of the Vernal Pool Fairy Shrimp Branchinecta lynchi. This branchiopod crustacean is endemic to California's freshwater ephemeral ponds. It faces enormous habitat loss and fragmentation as urbanization and agriculture have fundamentally changed the vernal pool landscape over the past 3 centuries. The assembled genome consists of 22 chromosome-length scaffolds that account for 96.85% of the total sequence. One hundred and ninety-five unscaffolded contigs comprise the rest of the genome's 575.6 Mb length. The genome is substantially complete with a BUSCO score of 90.0%. There is no immediately identifiable sex chromosome, typical for this class of organism. This new resource will permit researchers to better understand the adaptive capacity of this imperiled species, as well as answer lingering questions about anostracan physiology, sex determination, and development.

A chromosome-level reference genome for the Versatile Fairy Shrimp, Branchinecta lindahli.

We present the novel reference genome of the Versatile Fairy Shrimp, Branchinecta lindahli. The Versatile Fairy Shrimp is a freshwater anostracan crustacean found across the western United States from Iowa to Oregon and from Alberta to Baja California. It is an ephemeral pool specialist, living in prairie potholes, irrigation ditches, tire treads, vernal pools, and other temporary freshwater wetlands. Anostracan fairy shrimp are facing global declines with 3 species in California on the Endangered Species list. This species was included in the California Conservation Genomics Project to provide an easily accessible reference genome, and to provide whole-genome resources for a generalist species, which may lead to new insights into Anostracan resiliency in the face of climate change. The final gapped genome comprises 15 chromosome-length scaffolds covering 98.63% of the 384.8 Mb sequence length, and an additional 55 unscaffolded contigs.

Polygenic discrimination of migratory phenotypes in an estuarine forage fish.

Migration is a complex phenotypic trait with some species containing migratory and nonmigratory individuals. Such life history variation may be attributed in part to plasticity, epigenetics, or genetics. Although considered semianadromous, recent studies using otolith geochemistry have revealed life history variation within the critically endangered Delta Smelt. Broadly categorizable as migratory or freshwater residents, we examined Restriction site Associated DNA sequencing data to test for a relationship between genetic variation and migratory behaviors. As previously shown, we found no evidence for neutral population genetic structure within Delta Smelt; however, we found significant evidence for associations between genetic variants and life history phenotypes. Furthermore, discriminant analysis of principal components, hierarchical clustering, and machine learning resulted in accurate assignment of fish into the freshwater resident or migratory classes based on their genotypes. These results suggest the presence of adaptive genetic variants relating to life history variation within a panmictic population. Mechanisms that may lead to this observation are genotype dependent habitat choice and spatially variable selection, both of which could operate each generation and are not exclusive. Given that the population of cultured Delta Smelt are being used as a refugial population for conservation, as a supply for wild population supplementation, and currently represent the majority of all living individuals of this species, we recommend that the hatchery management strategy consider the frequencies of life history-associated alleles and how to maintain this important aspect of Delta Smelt biological variation while under captive propagation.

Environmental DNA sampling provides new management strategies for vernal pool branchiopods in California.

California's vernal pools are declining ecosystems that support valuable native plant and animal diversity. Vernal pool branchiopods are particularly at risk from vernal pool habitat loss and conservation efforts have targeted their long-term protection through the establishment of preserves and conservation banks. These conservation strategies require repeated, perpetual monitoring of preserved habitat, which is currently carried out through dip-net surveys and visual identification of specimens. Dip-netting may be destructive and frequently requires some sacrifice of protected species. Environmental DNA offers a new, modern method to monitor many protected freshwater organisms. We designed qPCR-based species-specific assays for four of California's vernal pool branchiopods: The Vernal Pool Fairy Shrimp Branchinecta lynchi (BRLY), the Midvalley Fairy Shrimp Branchinecta mesovallensis (BRME), and the Conservancy Fairy Shrimp Branchinecta conservatio (BRCO), and the Vernal Pool Tadpole Shrimp Lepidurus packardi (LEPA). We tested these assays using eDNA sampling protocols alongside traditional dip-net surveys to assess their viability as an alternative method to monitor vernal pool branchiopods. Based on occupancy modeling, each of our assays achieved a 95% or higher detection rate when using optimized sampling protocols.

Role of freshwater floodplain-tidal slough complex in the persistence of the endangered delta smelt.

Seasonal floodplain wetland is one of the most variable and diverse habitats found in coastal ecosystems, yet it is also one of the most highly altered by humans. The Yolo Bypass, the primary floodplain of the Sacramento River in California's Central Valley, USA, has been shown to provide various benefits to native fishes when inundated. However, the Yolo Bypass exists as a tidal dead-end slough during dry periods and its value to native fishes has been less studied in this state. During the recent drought (2012-2016), we found higher abundance of the endangered Delta Smelt (Hypomesus transpacificus), than the previous 14 years of fish monitoring within the Yolo Bypass. Meanwhile, Delta Smelt abundance elsewhere in the estuary was at record lows during this time. To determine the value of the Yolo Bypass as a nursery habitat for Delta Smelt, we compared growth, hatch dates, and diets of juvenile Delta Smelt collected within the Yolo Bypass with fish collected among other putative nursery habitats in the San Francisco Estuary between 2010 and 2016. Our results indicated that when compared to other areas of the estuary, fish in the Yolo Bypass spawned earlier, and offspring experienced both higher quality feeding conditions and growth rates. The occurrence of healthy juvenile Delta Smelt in the Yolo Bypass suggested that the region may have acted as a refuge for the species during the drought years of 2012-2016. However, our results also demonstrated that no single region provided the best rearing habitat for juvenile Delta Smelt. It will likely require a mosaic of habitats that incorporates floodplain-tidal sloughs in order to promote the resilience of this declining estuarine fish species.

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.

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.