A single generation of domestication heritably alters the expression of hundreds of genes

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.

Next-generation data filtering in the genomics era

Genomic data are ubiquitous across disciplines, from agriculture to biodiversity, ecology, evolution and human health. However, these datasets often contain noise or errors and are missing information that can affect the accuracy and reliability of subsequent computational analyses and conclusions. A key step in genomic data analysis is filtering - removing sequencing bases, reads, genetic variants and/or individuals from a dataset - to improve data quality for downstream analyses. Researchers are confronted with a multitude of choices when filtering genomic data; they must choose which filters to apply and select appropriate thresholds. To help usher in the next generation of genomic data filtering, we review and suggest best practices to improve the implementation, reproducibility and reporting standards for filter types and thresholds commonly applied to genomic datasets. We focus mainly on filters for minor allele frequency, missing data per individual or per locus, linkage disequilibrium and Hardy-Weinberg deviations. Using simulated and empirical datasets, we illustrate the large effects of different filtering thresholds on common population genetics statistics, such as Tajima's D value, population differentiation (FST), nucleotide diversity (π) and effective population size (Ne).

Selection among critically endangered landlocked salmon (Salmo salar m. sebago) families in survival and growth traits across early life stages and in different environments

Endangered wild fish populations are commonly supported by hatchery propagation. However, hatchery-reared fish experience very different selective pressures compared to their wild counterparts, potentially causing genotype-by-environment interactions (G × E) in essential fitness traits. We experimentally studied early selection in a critically endangered landlocked Atlantic salmon population, first from fertilization to the swim-up stage in a common hatchery setting, and thereafter until the age of 5 months in two contrasting rearing environments. Swim-up progeny were moved either to standard indoor hatchery tanks involving conventional husbandry or to seminatural outdoor channels providing only natural food. After the first summer, sampled survivors were assigned to their families by genotyping. Early survival until the swim-up stage was mostly determined by maternal effects, but also involved significant variation due to sires and full-sib families (potential genetic effects). High on-growing survival in hatchery tanks (88.7%) maintained a more even distribution among families (relative share 1.5%-4.2%) than the seminatural environment (0.0%-5.4%). This heterogeneity was mostly maternal, whereas no independent paternal effect occurred. Heritability estimates were high for body size traits in both environments (0.62-0.69). Genetic correlations between the environments were significantly positive for body size traits (0.67-0.69), and high body condition in hatchery was also genetically linked to rapid growth in the seminatural environment (0.54). Additive and phenotypic growth variation increased in the seminatural environment, but scaling effects probably played a less significant role for G × E, compared to re-ranking of genotypes. Our results suggest that not only maternal effects, but also genetic effects, direct selection according to the environmental conditions experienced. Consistently high genetic variation in growth implies that, despite its low overall genetic diversity and long history in captive rearing (>50 years), this landlocked Atlantic salmon population still possesses adaptive potential for response to change from hatchery rearing back to more natural conditions.

A single generation in the wild increases fitness for descendants of hatchery-origin Chinook salmon (Oncorhynchus tshawytscha)

Reintroduction is an important tool for the recovery of imperiled species. For threatened Pacific salmonids (Oncorhynchus spp.) species, hatchery-origin (HOR) individuals from a nearby source are often used to reestablish populations in vacant, historically occupied habitat. However, this approach is challenged by the relatively low reproductive success that HOR Pacific salmonids experience when they spawn in the wild, relative to their natural-origin (NOR) counterparts. In this study, we used genetic parentage analysis to compare the reproductive success of three groups of adult Chinook salmon (Oncorhynchus tshawytscha) reintroduced above Cougar Dam on the South Fork McKenzie River, Oregon: HOR Chinook salmon from an integrated stock; first-generation, wild-born descendants (hereafter F 1s) of Chinook salmon produced at the same hatchery; and NOR Chinook salmon that are presumed to have been produced below the dam, on the mainstem McKenzie River, or elsewhere and volitionally entered a trap below Cougar Dam. We found that F 1s produced nearly as many adult offspring as NORs, and 1.8-fold more adult offspring than HORs. This result suggests that, for the South Fork McKenzie reintroduction program, a single generation in the wild increases fitness for the descendants of HOR Chinook salmon. Although these results are encouraging, care must be taken before extrapolating our results to other systems.

Physiological differences between wild and captive animals: a century-old dilemma

Laboratory-based research dominates the fields of comparative physiology and biomechanics. The power of lab work has long been recognized by experimental biologists. For example, in 1932, Georgy Gause published an influential paper in Journal of Experimental Biology describing a series of clever lab experiments that provided the first empirical test of competitive exclusion theory, laying the foundation for a field that remains active today. At the time, Gause wrestled with the dilemma of conducting experiments in the lab or the field, ultimately deciding that progress could be best achieved by taking advantage of the high level of control offered by lab experiments. However, physiological experiments often yield different, and even contradictory, results when conducted in lab versus field settings. This is especially concerning in the Anthropocene, as standard laboratory techniques are increasingly relied upon to predict how wild animals will respond to environmental disturbances to inform decisions in conservation and management. In this Commentary, we discuss several hypothesized mechanisms that could explain disparities between experimental biology in the lab and in the field. We propose strategies for understanding why these differences occur and how we can use these results to improve our understanding of the physiology of wild animals. Nearly a century beyond Gause's work, we still know remarkably little about what makes captive animals different from wild ones. Discovering these mechanisms should be an important goal for experimental biologists in the future.

DNA hydroxymethylation differences underlie phenotypic divergence of somatic growth in Nile tilapia reared in common garden

Phenotypic plasticity of metabolism and growth are essential for adaptation to new environmental conditions, such as those experienced during domestication. Epigenetic regulation plays a key role in this process but the underlying mechanisms are poorly understood, especially in the case of hydroxymethylation. Using reduced representation 5-hydroxymethylcytosine profiling, we compared the liver hydroxymethylomes in full-sib Nile tilapia with distinct growth rates (3.8-fold difference) and demonstrated that DNA hydroxymethylation is strongly associated with phenotypic divergence of somatic growth during the early stages of domestication. The 2677 differentially hydroxymethylated cytosines between fast- and slow-growing fish were enriched within gene bodies (79%), indicating a pertinent role in transcriptional regulation. Moreover, they were found in genes involved in biological processes related to skeletal system and muscle structure development, and there was a positive association between somatic growth and 5hmC levels in genes coding for growth factors, kinases and receptors linked to myogenesis. Single nucleotide polymorphism analysis revealed no genetic differentiation between fast- and slow-growing fish. In addition to unveiling a new link between DNA hydroxymethylation and epigenetic regulation of growth in fish during the initial stages of domestication, this study suggests that epimarkers may be applied in selective breeding programmes for superior phenotypes.

Environmental and genetic effects of captivity - are there lessons for monarch butterfly conservation?

Rearing monarch butterflies in captivity for later release is a popular but contentious activity due to concerns about its potential negative effects on the wild population. In this review, I discuss how captive rearing and breeding could impact monarch fitness in the wild, the current evidence for such impacts in monarchs and other captive-reared/released organisms, and how this should inform our efforts to conserve monarchs and other species.

Multi-level assessment of the origin, feeding area and organohalogen contamination on salmon from the Baltic Sea

The Atlantic salmon (Salmo salar) population in the Baltic Sea consists of wild and hatchery-reared fish that have been released into the sea to support salmon stocks. During feeding migration, salmon migrate to different parts of the Baltic Sea and are exposed to various biotic and abiotic stressors, such as organohalogen compounds (OHCs). The effects of salmon origin (wild or hatchery-reared), feeding area (Baltic Main Basin, Bothnian Sea, and Gulf of Finland), and OHC concentration on the differences in hepatic proteome of salmon were investigated. Multi-level analysis of the OHC concentration, transcriptome, proteome, and oxidative stress biomarkers measured from the same salmon individuals were performed to find the key variables (origin, feeding area, OHC concentrations, and oxidative stress) that best account for the differences in the transcriptome and proteome between the salmon groups. When comparing wild and hatchery-reared salmon, differences were found in xenobiotic and amino acid metabolism-related pathways. When comparing salmon from different feeding areas, the amino acid and carbohydrate metabolic pathways were notably different. Several proteins found in these pathways are correlated with the concentrations of polychlorinated biphenyls (PCBs). The multi-level analysis also revealed amino acid metabolic pathways in connection with PCBs and oxidative stress variables related to glutathione metabolism. Other pathways found in the multi-level analysis included genetic information processes related to ribosomes, signaling and cellular processes related to the cytoskeleton, and the immune system, which were connected mainly to the concentrations of Polychlorinated biphenyls and Dichlorodiphenyltrichloroethane and their metabolites. These results suggest that the hepatic proteome of salmon in the Baltic Sea, together with the transcriptome, is more affected by the OHC concentrations and oxidative stress of the feeding area than the origin of the salmon.

Gene expression differentials driven by mass rearing and artificial selection in black soldier fly colonies

The micro-evolutionary forces that shape genetic diversity during domestication have been assessed in many plant and animal systems. However, the impact of these processes on gene expression, and consequent functional adaptation to artificial environments, remains under-investigated. In this study, whole-transcriptome dynamics associated with the early stages of domestication of the black soldier fly (BSF), Hermetia illucens, were assessed. Differential gene expression (DGE) was evaluated in relation to (i) generational time within the cultured environment (F2 vs. F3), and (ii) two selection strategies [no artificial selective pressure (NS); and selection for greater larval mass (SEL)]. RNA-seq was conducted on 5th instar BSF larvae (n = 36), representing equal proportions of the NS (F2 = 9; F3 = 9) and SEL (F2 = 9; F3 = 9) groups. A multidimensional scaling plot revealed greater gene expression variability within the NS and F2 subgroups, while the SEL group clustered separately with lower levels of variation. Comparisons between generations revealed 898 differentially expressed genes (DEGs; FDR-corrected p < 0.05), while between selection strategies, 213 DEGs were observed (FDR-corrected p < 0.05). Enrichment analyses revealed that metabolic, developmental, and defence response processes were over-expressed in the comparison between F2 and F3 larvae, while metabolic processes were the main differentiating factor between NS and SEL lines. This illustrates the functional adaptations that occur in BSF colonies across generations due to mass rearing; as well as highlighting genic dynamics associated with artificial selection for production traits that might inform future selective breeding strategies.

Genome-wide hydroxymethylation profiles in liver of female Nile tilapia with distinct growth performance

The mechanisms underlying the fast genome evolution that occurs during animal domestication are poorly understood. Here, we present a genome-wide epigenetic dataset that quantifies DNA hydroxymethylation at single nucleotide resolution among full-sib Nile tilapia (Oreochromis niloticus) with distinct growth performance. In total, we obtained 355 million, 75 bp reads from 5 large- and 5 small-sized fish on an Illumina NextSeq500 platform. We identified several growth-related genes to be differentially hydroxymethylated, especially within gene bodies and promoters. Previously, we proposed that DNA hydroxymethylation greatly affects the earliest responses to adaptation and potentially drives genome evolution through its targeted enrichment and elevated nucleotide transversion rates. This dataset can be analysed in various contexts (e.g., epigenetics, evolution and growth) and compared to other epigenomic datasets in the future, namely DNA methylation and histone modifications. With forthcoming advancements in genome research, this hydroxymethylation dataset will also contribute to better understand the epigenetic regulation of key genomic features, such as cis-regulatory and transposable elements.

Captive rearing effects on the methylome of Atlantic salmon after oceanic migration: Sex-specificity and intergenerational stability

Captive rearing in salmon hatcheries can have considerable impacts on both fish phenotype and fitness within a single generation, even in the absence of genetic change. Evidence for hatchery-induced changes in DNA methylation is becoming abundant, though questions remain on the sex-specificity of these effects, their persistence until spawning and potential for transmission to future generations. Here we performed whole genome methylation sequencing of fin tissue for 16 hatchery and 16 wild Atlantic salmon (Salmo salar) returning to spawn in the Rimouski River, Québec, Canada. We identified two cohorts of hatchery-reared salmon through methylation analysis, one of which was epigenetically similar to wild fish, suggesting that supplementation efforts may be able to minimize the epigenetic effects of hatchery rearing. We found considerable sex-specific effects of hatchery rearing, with few genomic regions being affected in both males and females. We also analysed the methylome of 32 F1 offspring from four groups (pure wild, pure hatchery origin and reciprocal hybrids). We found that few epigenetic changes due to parental hatchery rearing persisted in the F1 offspring though the patterns of inheritance appear to be complex, involving nonadditive effects. Our results suggest that the epigenetic effects of hatchery rearing can be minimal in F0 . There may also be minimal epigenetic inheritance and rapid loss of epigenetic changes associated with hatchery rearing. However, due to sex-specificity and nonadditive patterns of inheritance, methylation changes due to captive rearing are rather complex and the field would benefit from further research on minimizing the epigenetic effects of captive rearing in conservation efforts.

Influence of RNA-Seq library construction, sampling methods, and tissue harvesting time on gene expression estimation

RNA sequencing (RNA-Seq) is popular for measuring gene expression in non-model organisms, including wild populations. While RNA-Seq can detect gene expression variation among wild-caught individuals and yield important insights into biological function, sampling methods can also affect gene expression estimates. We examined the influence of multiple technical variables on estimated gene expression in a non-model fish, the westslope cutthroat trout (Oncorhynchus clarkii lewisi), using two RNA-Seq library types: 3' RNA-Seq (QuantSeq) and whole mRNA-Seq (NEB). We evaluated effects of dip netting versus electrofishing, and of harvesting tissue immediately versus 5 min after euthanasia on estimated gene expression in blood, gill, and muscle. We found no significant differences in gene expression between sampling methods or tissue collection times with either library type. When library types were compared using the same blood samples, 58% of genes detected by both NEB and QuantSeq showed significantly different expression between library types, and NEB detected 31% more genes than QuantSeq. Although the two library types recovered different numbers of genes and expression levels, results with NEB and QuantSeq were consistent in that neither library type showed differences in gene expression between sampling methods and tissue harvesting times. Our study suggests that researchers can safely rely on different fish sampling strategies in the field. In addition, while QuantSeq is more cost effective, NEB detects more expressed genes. Therefore, when it is crucial to detect as many genes as possible (especially low expressed genes), when alternative splicing is of interest, or when working with an organism lacking good genomic resources, whole mRNA-Seq is more powerful.

A Model-Based Hierarchical Bayesian Approach to Sholl Analysis

Due to the link between microglial morphology and function, morphological changes in microglia are frequently used to identify pathological immune responses in the central nervous system. In the absence of pathology, microglia are responsible for maintaining homeostasis, and their morphology can be indicative of how the healthy brain behaves in the presence of external stimuli and genetic differences. Despite recent interest in high throughput methods for morphological analysis, Sholl analysis is still the gold standard for quantifying microglia morphology via imaging data. Often, the raw data are naturally hierarchical, minimally including many cells per image and many images per animal. However, existing methods for performing downstream inference on Sholl data rely on truncating this hierarchy so rudimentary statistical testing procedures can be used. To fill this longstanding gap, we introduce a fully parametric model-based approach for analyzing Sholl data. We generalize our model to a hierarchical Bayesian framework so that inference can be performed without aggressive reduction of otherwise very rich data. We apply our model to three real data examples and perform simulation studies comparing the proposed method with a popular alternative.

Environment and genotype predict the genomic nature of domestication of salmonids as revealed by gene expression

Billions of salmonids are produced annually by artificial reproduction for harvest and conservation. Morphologically, behaviourally and physiologically these fish differ from wild-born fish, including in ways consistent with domestication. Unlike most studied domesticates, which diverged from wild ancestors millennia ago, salmonids offer a tractable model for early-stage domestication. Here, we review a fundamental mechanism for domestication-driven differences in early-stage domestication, differentially expressed genes (DEGs), in salmonids. We found 34 publications examining DEGs under domestication driven by environment and genotype, covering six species, over a range of life-history stages and tissues. Three trends emerged. First, domesticated genotypes have increased expression of growth hormone and related metabolic genes, with differences magnified under artificial environments with increased food. Regulatory consequences of these DEGs potentially drive overall DEG patterns. Second, immune genes are often DEGs under domestication and not simply owing to release from growth-immune trade-offs under increased food. Third, domesticated genotypes exhibit reduced gene expression plasticity, with plasticity further reduced in low-complexity environments typical of production systems. Recommendations for experimental design improvements, coupled with tissue-specific expression and emerging analytical approaches for DEGs present tractable avenues to understand the evolution of domestication in salmonids and other species.

Captivity induces large and population-dependent brain transcriptomic changes in wild-caught cane toads (Rhinella marina)

Gene expression levels are key molecular phenotypes at the interplay between genotype and environment. Mounting evidence suggests that short-term changes in environmental conditions, such as those encountered in captivity, can substantially affect gene expression levels. Yet, the exact magnitude of this effect, how general it is, and whether it results in parallel changes across populations are not well understood. Here, we take advantage of the well-studied cane toad, Rhinella marina, to examine the effect of short-term captivity on brain gene expression levels, and determine whether effects of captivity differ between long-colonized and vanguard populations of the cane toad's Australian invasion range. We compared the transcriptomes of wild-caught toads immediately assayed with those from toads captured from the same populations but maintained in captivity for seven months. We found large differences in gene expression levels between captive and wild-caught toads from the same population, with an over-representation of processes related to behaviour and the response to stress. Captivity had a much larger effect on both gene expression levels and gene expression variability in toads from vanguard populations compared to toads from long-colonized areas, potentially indicating an increased plasticity in toads at the leading edge of the invasion. Overall, our findings indicate that short-term captivity can induce large and population-specific transcriptomic changes, which has significant implications for studies comparing phenotypic traits of wild-caught organisms from different populations that have been held in captivity.

Longitudinal monitoring of neutral and adaptive genomic diversity in a reintroduction

Restoration programs in the form of ex-situ breeding combined with reintroductions are becoming critical to counteract demographic declines and species losses. Such programs are increasingly using genetic management to improve conservation outcomes. However, the lack of long-term monitoring of genetic indicators following reintroduction prevents assessments of the trajectory and persistence of reintroduced populations. We carried out an extensive monitoring program in the wild for a threatened small-bodied fish (southern pygmy perch, Nannoperca australis) to assess the long-term genomic effects of its captive breeding and reintroduction. The species was rescued prior to its extirpation from the terminal lakes of Australia's Murray-Darling Basin, and then used for genetically informed captive breeding and reintroductions. Subsequent annual or biannual monitoring of abundance, fitness, and occupancy over a period of 11 years, combined with postreintroduction genetic sampling, revealed survival and recruitment of reintroduced fish. Genomic analyses based on data from the original wild rescued, captive born, and reintroduced cohorts revealed low inbreeding and strong maintenance of neutral and candidate adaptive genomic diversity across multiple generations. An increasing trend in the effective population size of the reintroduced population was consistent with field monitoring data in demonstrating successful re-establishment of the species. This provides a rare empirical example that the adaptive potential of a locally extinct population can be maintained during genetically informed ex-situ conservation breeding and reintroduction into the wild. Strategies to improve biodiversity restoration via ex-situ conservation should include genetic-based captive breeding and longitudinal monitoring of standing genomic variation in reintroduced populations.

Dual resistance to Flavobacterium psychrophilum and Myxobolus cerebralis in rainbow trout (Oncorhynchus mykiss, Walbaum)

Aquatic pathogens are a major concern for fish hatchery production, fisheries management, and conservation, and disease control needs to be addressed. Two important salmonid pathogens are Myxobolus cerebralis and Flavobacterium psychrophilum that cause whirling disease and bacterial coldwater disease (BCWD), respectively. Innate disease resistance is a potential option for reducing disease-related mortality in hatchery-reared rainbow trout (Oncorhynchus mykiss, Walbaum). Two experiments were conducted to assess pathogen resistance of first-generation (F1) rainbow trout created by crossing M. cerebralis- and F. psychrophilum-resistant strains. In the first experiment, we exposed two rainbow trout strains and one F1 cross to six treatments: control (no exposure), mock injection, F. psychrophilum only, M. cerebralis only, F. psychrophilum then M. cerebralis, and M. cerebralis then F. psychrophilum. Results indicated that the F1 cross was not resistant to either pathogen. In the second experiment, we exposed five rainbow trout strains and four rainbow trout crosses to F. psychrophilum. The second experiment indicated that at least one rainbow trout cross was F. psychrophilum-resistant. Achieving dual resistance may be possible using selective breeding but only some multigenerational strains are suitable candidates for further evaluation.

Effects of stocking at the parr stage on the reproductive fitness and genetic diversity of a wild population of Atlantic salmon (Salmo salar L.)

Captive-breeding programs are among the most adopted conservation practices to mitigate the loss of biodiversity, including genetic diversity. However, both genetic and nongenetic changes occurring in captivity can reduce the fitness of supplemented individuals, which complicate rehabilitation efforts. In the case of Atlantic salmon, the intensity of changes that occur in captivity and their impact on fitness will vary with the stocking practice adopted. In this study, we test whether salmon stocked at the parr stage have reduced reproductive success compared with their wild conspecifics and whether they contribute to increase genetic diversity in the targeted population. To do so, we use high-throughput microsatellite sequencing of 38 loci to accurately assign 2381 offspring to a comprehensive set of possible parents from a supplemented Atlantic salmon population in Québec, Canada. Captive-bred salmon stocked at the parr stage had fewer mates than their wild conspecifics, as well as a reduced relative reproductive success (RSS) compared with their wild counterparts. Nonetheless, in comparison with previous studies, stocking at the parr stage significantly improved RSS compared with salmon stocked as smolts and they displayed a reduction in reproductive success similar to salmon stocked as fry, which spend less time in captivity than parr. Moreover, supplementation of captive-bred salmon significantly contributed to increasing genetic diversity. These results should contribute to informing resource managers in determining the best stocking practice to enhance Atlantic salmon populations.

Gut Microbial Composition of Pacific Salmonids Differs across Oregon River Basins and Hatchery Ancestry

The gut microbiome may represent a relatively untapped resource in the effort to manage and conserve threatened or endangered fish populations, including wild and hatchery-reared Pacific salmonids. To clarify this potential, we defined how steelhead trout gut microbiome composition varies across watersheds and as a function of ancestry. First, we measured this variation across watersheds using wild steelhead trout sampled from nine locations spanning three river basins. While gut microbial composition differs across basins, there exist bacterial clades that are ubiquitous across all populations. Correlating the phylogenetic composition of clades with geographic distance reveals 395 clades of bacteria whose ecological distribution implicates their co-diversification with steelheads. Second, we quantified how microbiome composition varies between first generation hatchery-reared steelhead and traditional hatchery-reared steelhead. Despite being subject to the same hatchery management strategies, fish bred from wild parents carry distinct microbiomes from those bred from hatchery broodstock, implicating the role of genotype on microbiome composition. Finally, we integrated all data from both studies to reveal two distinct, yet robust clusters of community composition. Collectively, our study documents for the first time how the steelhead gut microbiome varies by geography or broodstock and uncovers microbial taxa that may indicate the watershed or hatchery from which an individual was sourced.

Strong parallel differential gene expression induced by hatchery rearing weakly associated with methylation signals in adult Coho Salmon (O. kisutch)

Human activities and resource exploitation led to a massive decline of wild salmonid populations, consequently, numerous conservation programs have been developed to supplement wild populations. However, many studies documented reduced fitness of hatchery-born relative to wild fish. Here, by using both RNA sequencing and Whole Genome Bisulfite Sequencing of hatchery and wild-born adult Coho salmon (Oncorhynchus kisutch) originating from two previously studied river systems, we show that early-life hatchery-rearing environment-induced significant and parallel gene expression differentiation is maintained until Coho come back to their natal river for reproduction. A total of 3,643 genes differentially expressed and 859 coexpressed genes were downregulated in parallel in hatchery-born fish from both rivers relative to their wild congeners. Among those genes, 26 displayed a significant relationship between gene expression and the median gene body methylation and 669 single CpGs displayed a significant correlation between methylation level and the associated gene expression. The link between methylation and gene expression was weak suggesting that DNA methylation is not the only player in mediating hatchery-related expression differences. Yet, significant gene expression differentiation was observed despite 18 months spent in a common environment (i.e., the sea). Finally, the differentiation is observed in parallel in two different river systems, highlighting the fact that early-life environment may account for at least some of the reduced fitness of the hatchery salmon in the wild. These results illustrate the relevance and importance of considering both epigenome and transcriptome to evaluate the costs and benefits of large-scale supplementation programs.

Distinct early life stage gene expression effects of hybridization among European and North American farmed and wild Atlantic salmon populations

Due to multi-generation domestication selection, farmed and wild Atlantic salmon diverge genetically, which raises concerns about potential genetic interactions among escaped farmed and wild populations and disruption of local adaptation through introgression. When farmed strains of distant geographic origin are used, it is unknown whether the genetic consequences posed by escaped farmed fish will be greater than if more locally derived strains are used. Quantifying gene transcript expression differences among divergent farmed, wild and F₁ hybrids under controlled conditions is one of the ways to explore the consequences of hybridization. We compared the transcriptomes of fry at the end of yolk sac absorption of a European (EO) farmed ("StofnFiskur", Norwegian strain), a North American (NA) farmed (Saint John River, NB strain), a Newfoundland (NF) wild population with EO ancestry, and related F₁ hybrids using 44K microarrays. Our findings indicate that the wild population showed greater transcriptome differences from the EO farmed strain than that of the NA farmed strain. We also found the largest differences in global gene expression between the two farmed strains. We detected the fewest differentially expressed transcripts between F₁ hybrids and domesticated/wild maternal strains. We also found that the differentially expressed genes between cross types over-represented GO terms associated with metabolism, development, growth, immune response, and redox homeostasis processes. These findings suggest that the interbreeding of escaped EO/NA farmed and NF wild population would alter gene transcription, and the consequences of hybridization would be greater from escaped EO farmed than NA farmed salmon, resulting in potential effects on the wild populations.

Early fish domestication affects methylation of key genes involved in the rapid onset of the farmed phenotype

Animal domestication is a process of environmental modulation and artificial selection leading to permanent phenotypic modifications. Recent studies showed that phenotypic changes occur very early in domestication, i.e., within the first generation in captivity, which raises the hypothesis that epigenetic mechanisms may play a critical role on the early onset of the domestic phenotype. In this context, we applied reduced representation bisulphite sequencing to compare methylation profiles between wild Nile tilapia females and their offspring reared under farmed conditions. Approximately 700 differentially methylated CpG sites were found, many of them associated not only with genes involved in muscle growth, immunity, autophagy and diet response but also related to epigenetic mechanisms, such as RNA methylation and histone modifications. This bottom-up approach showed that the phenotypic traits often related to domestic animals (e.g., higher growth rate and different immune status) may be regulated epigenetically and prior to artificial selection on gene sequences. Moreover, it revealed the importance of diet in this process, as reflected by differential methylation patterns in genes critical to fat metabolism. Finally, our study highlighted that the TGF-β1 signalling pathway may regulate and be regulated by several differentially methylated CpG-associated genes. This could be an important and multifunctional component in promoting adaptation of fish to a domestic environment while modulating growth and immunity-related traits.

Environment-driven reprogramming of gamete DNA methylation occurs during maturation and is transmitted intergenerationally in Atlantic Salmon

An epigenetic basis for transgenerational plasticity in animals is widely theorized, but convincing empirical support is limited by taxa-specific differences in the presence and role of epigenetic mechanisms. In teleost fishes, DNA methylation generally does not undergo extensive reprogramming and has been linked with environmentally induced intergenerational effects, but solely in the context of early life environmental differences. Using whole-genome bisulfite sequencing, we demonstrate that differential methylation of sperm occurs in response to captivity during the maturation of Atlantic Salmon (Salmo salar), a species of major economic and conservation significance. We show that adult captive exposure further induces differential methylation in an F1 generation that is associated with fitness-related phenotypic differences. Some genes targeted with differential methylation were consistent with genes differential methylated in other salmonid fishes experiencing early-life hatchery rearing, as well as genes under selection in domesticated species. Our results support a mechanism of transgenerational plasticity mediated by intergenerational inheritance of DNA methylation acquired late in life for salmon. To our knowledge, this is the first-time environmental variation experienced later in life has been directly demonstrated to influence gamete DNA methylation in fish.

Captive breeding and the conservation of the threatened houbara bustards

Translocation of captive-bred individuals to reinforce wild populations may be an important conservation approach for some species, but can be detrimental when employed to boost exploited wild populations, particularly where repeated long-term reinforcement aims to compensate for repeated unregulated offtake. We review evidence that captive breeding alters multiple physiological, life-history and temperamental traits through founder effects, genetic drift and unintended adaption to captivity; degrades learnt behaviours; and compromises biogeography, population structure and viability through introgression. We highlight these risks for the globally threatened African houbara Chlamydotis undulata and Asian houbara C. macqueenii, 2 bustard species hunted throughout much of their ranges and now subject to multiple large-scale captive-breeding programmes and translocations. In eastern Morocco, annual releases of captive-bred African houbara are 2‒3 times higher than original wild numbers, but no investigation of their potentially deleterious effects has, to our knowledge, been published, although most wild populations may now have been replaced by captive-bred domestic stock, which are reportedly not self-sustaining. Despite multiple decades of reinforcement, we are not aware of any analysis of the contribution of captive breeding to African houbara population dynamics, or of the genomic consequences. Asian houbara release programmes may also be promoting rather than preventing declines, and need to contextualise themselves through rigorous analyses of wild population numbers, demographic rates and threats, maintenance of phylogeographic concordance of released with supplemented populations, profiling of traits crucial to survival and the measurement and modelling of the impacts of reinforcement on physiological and behavioural fitness of wild populations.

Epigenomic modifications induced by hatchery rearing persist in germ line cells of adult salmon after their oceanic migration

Human activities induce direct or indirect selection pressure on natural population and may ultimately affect population's integrity. While numerous conservation programs aimed to minimize human-induced genomic variation, human-induced environmental variation may generate epigenomic variation potentially affecting fitness through phenotypic modifications. Major questions remain pertaining to how much epigenomic variation arises from environmental heterogeneity, whether this variation can persist throughout life, and whether it can be transmitted across generations. We performed whole genome bisulfite sequencing (WGBS) on the sperm of genetically indistinguishable hatchery and wild-born migrating adults of Coho salmon (Oncorhynchus kisutch) from two geographically distant rivers at different epigenome scales. Our results showed that coupling WGBS with fine-scale analyses (local and chromosomal) allows the detection of parallel early-life hatchery-induced epimarks that differentiate wild from hatchery-reared salmon. Four chromosomes and 183 differentially methylated regions (DMRs) displayed a significant signal of methylation differentiation between hatchery and wild-born Coho salmon. Moreover, those early-life epimarks persisted in germ line cells despite about 1.5 year spent in the ocean following release from hatchery, opening the possibility for transgenerational inheritance. Our results strengthen the hypothesis that epigenomic modifications environmentally induced during early-life development persist in germ cells of adults until reproduction, which could potentially impact their fitness.

Redox-related gene expression and sugar accumulation patterns are altered in the edible inflorescence produced by the cultivated form of pacaya palm (Chamaedorea tepejilote)

BACKGROUND AND AIMS

The pacaya palm is a dioecious neotropical palm species that is exploited in Latin America for its male inflorescence, which is edible when immature. It is cultivated, in a non-intensive manner, in Guatemala, where a morphotype occurs that produces much larger, more highly branched inflorescences compared with wild palms. We sought to identify molecular factors underlying this phenotypic divergence, which is likely to be a product of domestication.

METHODS

We performed RNA-seq-based studies on immature pacaya palm male inflorescences in order to identify genes that might be directly or indirectly affected in their expression in relation to domestication. We also measured the accumulation of a range of soluble sugar molecules to provide information on the biochemical status of the two different types of material.

KEY RESULTS

A total of 408 genes were found to display significantly different expression levels between the wild and cultivated morphotypes. Three different functional categories were found to be enriched in the gene set that was upregulated in the cultivated morphotype: redox balance; secondary metabolism; and transport. Several sugars were found to accumulate at higher levels in inflorescences of the cultivated morphotype, in particular myo-inositol, fructose and glucose.

CONCLUSIONS

The observed upregulation of redox-related genes in the cultivated morphotype is corroborated by the observation of higher myo-inositol accumulation, which has been shown to be associated with enhanced scavenging of reactive oxygen species in other plants and which may affect meristem activity.

Population differences in Chinook salmon (Oncorhynchus tshawytscha) DNA methylation: Genetic drift and environmental factors

Local adaptation and phenotypic differences among populations have been reported in many species, though most studies focus on either neutral or adaptive genetic differentiation. With the discovery of DNA methylation, questions have arisen about its contribution to individual variation in and among natural populations. Previous studies have identified differences in methylation among populations of organisms, although most to date have been in plants and model animal species. Here we obtained eyed eggs from eight populations of Chinook salmon (Oncorhynchus tshawytscha) and assayed DNA methylation at 23 genes involved in development, immune function, stress response, and metabolism using a gene-targeted PCR-based assay for next-generation sequencing. Evidence for population differences in methylation was found at eight out of 23 gene loci after controlling for developmental timing in each individual. However, we found no correlation between freshwater environmental parameters and methylation variation among populations at those eight genes. A weak correlation was identified between pairwise DNA methylation dissimilarity among populations and pairwise F ST based on 15 microsatellite loci, indicating weak effects of genetic drift or geographic distance on methylation. The weak correlation was primarily driven by two genes, GTIIBS and Nkef. However, single-gene Mantel tests comparing methylation and pairwise F ST were not significant after Bonferroni correction. Thus, population differences in DNA methylation are more likely related to unmeasured oceanic environmental conditions, local adaptation, and/or genetic drift. DNA methylation is an additional mechanism that contributes to among population variation, with potential influences on organism phenotype, adaptive potential, and population resilience.

The use of non-lethal sampling for transcriptomics to assess the physiological status of wild fishes

Fishes respond to different abiotic and biotic stressors through changes in gene expression as a part of an integrated physiological response. Transcriptomics approaches have been used to quantify gene expression patterns as a reductionist approach to understand responses to environmental stressors in animal physiology and have become more commonly used to study wild fishes. We argue that non-lethal sampling for transcriptomics should become the norm for assessing the physiological status of wild fishes, especially when there are conservation implications. Processes at the level of the transcriptome provide a "snapshot" of the cellular conditions at a given time; however, by using a non-lethal sampling protocol, researchers can connect the transcriptome profile with fitness-relevant ecological endpoints such as reproduction, movement patterns and survival. Furthermore, telemetry is a widely used approach in fisheries to understand movement patterns in the wild, and when combined with transcriptional profiling, provides arguably the most powerful use of non-lethal sampling for transcriptomics in wild fishes. In this review, we discuss the different tissues that can be successfully incorporated into non-lethal sampling strategies, which is particularly useful in the context of the emerging field of conservation transcriptomics. We briefly describe different methods for transcriptional profiling in fishes from high-throughput qPCR to whole transcriptome approaches. Further, we discuss strategies and the limitations of using transcriptomics for non-lethally studying fishes. Lastly, as 'omics' technology continues to advance, transcriptomics paired with different omics approaches to study wild fishes will provide insight into the factors that regulate phenotypic variation and the physiological responses to changing environmental conditions in the future.

Offspring survival changes over generations of captive breeding

Conservation breeding programs such as zoos play a major role in preventing extinction, but their sustainability may be impeded by neutral and adaptive population genetic change. These changes are difficult to detect for a single species or context, and impact global conservation efforts. We analyse pedigree data from 15 vertebrate species – over 30,000 individuals – to examine offspring survival over generations of captive breeding. Even accounting for inbreeding, we find that the impacts of increasing generations in captivity are highly variable across species, with some showing substantial increases or decreases in offspring survival over generations. We find further differences between dam and sire effects in first- versus multi-generational analysis. Crucially, our multispecies analysis reveals that responses to captivity could not be predicted from species’ evolutionary (phylogenetic) relationships. Even under best-practice captive management, generational fitness changes that cannot be explained by known processes (such as inbreeding depression), are occurring.

Captive breeding could prevent species extinctions, but selection for captivity may decrease fitness. Here the authors analyse pedigree data on 15 long-running vertebrate breeding programs and find generational fitness changes that processes such as inbreeding depression cannot explain.

Population Genomics Training for the Next Generation of Conservation Geneticists: ConGen 2018 Workshop

The increasing availability and complexity of next-generation sequencing (NGS) data sets make ongoing training an essential component of conservation and population genetics research. A workshop entitled “ConGen 2018” was recently held to train researchers in conceptual and practical aspects of NGS data production and analysis for conservation and ecological applications. Sixteen instructors provided helpful lectures, discussions, and hands-on exercises regarding how to plan, produce, and analyze data for many important research questions. Lecture topics ranged from understanding probabilistic (e.g., Bayesian) genotype calling to the detection of local adaptation signatures from genomic, transcriptomic, and epigenomic data. We report on progress in addressing central questions of conservation genomics, advances in NGS data analysis, the potential for genomic tools to assess adaptive capacity, and strategies for training the next generation of conservation genomicists.

Transcriptome Analyses Throughout Chili Pepper Fruit Development Reveal Novel Insights into the Domestication Process

Chili pepper (Capsicum spp.) is an important crop, as well as a model for fruit development studies and domestication. Here, we performed a time-course experiment to estimate standardized gene expression profiles with respect to fruit development for six domesticated and four wild chili pepper ancestors. We sampled the transcriptomes every 10 days from flowering to fruit maturity, and found that the mean standardized expression profiles for domesticated and wild accessions significantly differed. The mean standardized expression was higher and peaked earlier for domesticated vs. wild genotypes, particularly for genes involved in the cell cycle that ultimately control fruit size. We postulate that these gene expression changes are driven by selection pressures during domestication and show a robust network of cell cycle genes with a time shift in expression, which explains some of the differences between domesticated and wild phenotypes.

Farm Animals Are Long Away from Natural Behavior: Open Questions and Operative Consequences on Animal Welfare

Simple Summary

Animal welfare is a very important issue. One of the tasks of researchers is to provide explanations and possible solutions to questions arising from non-experts. This work analyzes part of the extensive literature on relationships between selection and domestic, mainly farm, animals’ behavior and deals with some very important themes, such as the role of regulations, domestication, and selection.

Abstract

The concept of welfare applied to farm animals has undergone a remarkable evolution. The growing awareness of citizens pushes farmers to guarantee the highest possible level of welfare to their animals. New perspectives could be opened for animal welfare reasoning around the concept of domestic, especially farm, animals as partial human artifacts. Therefore, it is important to understand how much a particular behavior of a farm animal is far from the natural one of its ancestors. This paper is a contribution to better understand the role of genetics of the farm animals on their behavior. This means that the naïve approach to animal welfare regarding returning animals to their natural state should be challenged and that welfare assessment should be considered.

Incipient resistance to an effective pesticide results from genetic adaptation and the canalization of gene expression

The resistance of pest species to chemical controls has vast ecological, economic, and societal costs. In most cases, resistance is only detected after spreading throughout an entire population. Detecting resistance in its incipient stages, by comparison, provides time to implement preventative strategies. Incipient resistance can be detected by coupling standard toxicology assays with large-scale gene expression experiments. We apply this approach to a system where an invasive parasite, sea lamprey (Petromyzon marinus), has been treated with the highly effective pesticide 3-trifluoromethyl-4-nitrophenol (TFM) for 60 years. Toxicological experiments revealed that lamprey from treated populations did not have higher survival to TFM exposure than lamprey from untreated populations, demonstrating that full-fledged resistance has not yet evolved. In contrast, we find hundreds of genes differentially expressed in response to TFM in the population with the longest history of exposure, many of which relate to TFM's primary mode of action, the uncoupling of oxidative phosphorylation, and subsequent depletion of ATP. Three genes critical to oxidative phosphorylation, ATP5PB, PLCB1, and NDUFA9, were nearly fixed for alternative alleles in comparisons of SNPs between treated and untreated populations (FST  > 5 SD from the mean). ATP5PB encodes subunit b of ATP synthase and an additional subunit, ATP5F1B, was canalized for high expression in treated populations, but remained plastic in response to TFM treatment in individuals from the untreated population. These combined genomic and transcriptomic results demonstrate that an adaptive, genetic response to TFM is likely driving incipient resistance in a damaging pest species.

Does parental angling selection affect the behavior or metabolism of brown trout parr?

The behavior of organisms can be subject to human-induced selection such as that arising from fishing. Angling is expected to induce mortality on fish with bold and explorative behavior, which are behaviors commonly linked to a high standard metabolic rate. We studied the transgenerational response of brown trout (Salmo trutta) to angling-induced selection by examining the behavior and metabolism of 1-year-old parr between parents that were or were not captured by experimental fly fishing. We performed the angling selection experiment on both a wild and a captive population, and compared the offspring for standard metabolic rate and behavior under predation risk in common garden conditions. Angling had population-specific effects on risk taking and exploration tendency, but no effects on standard metabolic rate. Our study adds to the evidence that angling can induce transgenerational responses on fish personality. However, understanding the mechanisms of divergent responses between the populations requires further study on the selectivity of angling in various conditions.

A Bioinformatics Model of Human Diseases on the Basis of Differentially Expressed Genes (of Domestic Versus Wild Animals) That Are Orthologs of Human Genes Associated with Reproductive-Potential Changes

Earlier, after our bioinformatic analysis of single-nucleotide polymorphisms of TATA-binding protein-binding sites within gene promoters on the human Y chromosome, we suggested that human reproductive potential diminishes during self-domestication. Here, we implemented bioinformatics models of human diseases using animal in vivo genome-wide RNA-Seq data to compare the effect of co-directed changes in the expression of orthologous genes on human reproductive potential and during the divergence of domestic and wild animals from their nearest common ancestor (NCA). For example, serotonin receptor 3A (HTR3A) deficiency contributes to sudden death in pregnancy, consistently with Htr3a underexpression in guinea pigs (Cavia porcellus) during their divergence from their NCA with cavy (C. aperea). Overall, 25 and three differentially expressed genes (hereinafter, DEGs) in domestic animals versus 11 and 17 DEGs in wild animals show the direction consistent with human orthologous gene-markers of reduced and increased reproductive potential. This indicates a reliable association between DEGs in domestic animals and human orthologous genes reducing reproductive potential (Pearson’s χ² test p < 0.001, Fisher’s exact test p < 0.05, binomial distribution p < 0.0001), whereas DEGs in wild animals uniformly match human orthologous genes decreasing and increasing human reproductive potential (p > 0.1; binomial distribution), thus enforcing the norm (wild type).

Correspondence between a recreational fishery index and ecological condition for U.S.A. streams and rivers

Sport fishing is an important recreational and economic activity, especially in Australia, Europe and North America, and the condition of sport fish populations is a key ecological indicator of water body condition for millions of anglers and the public. Despite its importance as an ecological indicator representing the status of sport fish populations, an index for measuring this ecosystem service has not been quantified by analyzing actual fish taxa, size and abundance data across the U.S.A. Therefore, we used game fish data collected from 1,561 stream and river sites located throughout the conterminous U.S.A. combined with specific fish species and size dollar weights to calculate site-specific recreational fishery index (RFI) scores. We then regressed those scores against 38 potential site-specific environmental predictor variables, as well as site-specific fish assemblage condition (multimetric index; MMI) scores based on entire fish assemblages, to determine the factors most associated with the RFI scores. We found weak correlations between RFI and MMI scores and weak to moderate correlations with environmental variables, which varied in importance with each of 9 ecoregions. We conclude that the RFI is a useful indicator of a stream ecosystem service, which should be of greater interest to the USA public and traditional fishery management agencies than are MMIs, which tend to be more useful for ecologists, environmentalists and environmental quality agencies.

Rearing environment affects the genetic architecture and plasticity of DNA methylation in Chinook salmon

Genetic architecture and phenotypic plasticity are important considerations when studying trait variation within and among populations. Since environmental change can induce shifts in the genetic architecture and plasticity of traits, it is important to consider both genetic and environmental sources of phenotypic variation. While there is overwhelming evidence for environmental effects on phenotype, the underlying mechanisms are less clear. Variation in DNA methylation is a potential mechanism mediating environmental effects on phenotype due to its sensitivity to environmental stimuli, transgenerational inheritance, and influences on transcription. To characterize the effect of environment on methylation, we created two 6 × 6 (North Carolina II) Chinook salmon breeding crosses and reared the offspring in two environments: uniform hatchery tanks and seminatural stream channels. We sampled the fish twice during development, at the alevin (larval) and fry (juvenile) stages. We measured DNA methylation at 13 genes using a PCR-based bisulfite sequencing protocol. The genetic architecture of DNA methylation differed between rearing environments, with greater additive and nonadditive genetic variance in hatchery fish and greater maternal effects in seminatural channel fish, though gene-specific variation was evident. We observed plasticity in methylation across all assayed genes, as well as gene-specific effects at two genes in alevin and six genes in fry, indicating developmental stage-specific effects of rearing environment on methylation. Characterizing genetic and environmental influences on methylation is critical for future studies on DNA methylation as a potential mechanism for acclimation and adaptation.

Differential DNA methylation in somatic and sperm cells of hatchery vs wild (natural-origin) steelhead trout populations

Environmental factors such as nutrition, stress, and toxicants can influence epigenetic programming and phenotypes of a wide variety of species from plants to humans. The current study was designed to investigate the impacts of hatchery spawning and rearing on steelhead trout (Oncorhynchus mykiss) vs the wild fish on a molecular level. Additionally, epigenetic differences between feeding practices that allow slow growth (2 years) and fast growth (1 year) hatchery trout were investigated. The sperm and red blood cells (RBC) from adult male slow growth/maturation hatchery steelhead, fast growth/maturation hatchery steelhead, and wild (natural-origin) steelhead were collected for DNA preparation to investigate potential alterations in differential DNA methylation regions (DMRs) and genetic mutations, involving copy number variations (CNVs). The sperm and RBC DNA both had a large number of DMRs when comparing the hatchery vs wild steelhead trout populations. The DMRs were cell type specific with negligible overlap. Slow growth/maturation compared to fast growth/maturation steelhead also had a larger number of DMRs in the RBC samples. A number of the DMRs had associated genes that were correlated to various biological processes and pathologies. Observations demonstrate a major epigenetic programming difference between the hatchery and wild natural-origin fish populations, but negligible genetic differences. Therefore, hatchery conditions and growth/maturation rate can alter the epigenetic developmental programming of the steelhead trout. Interestingly, epigenetic alterations in the sperm allow for potential epigenetic transgenerational inheritance of phenotypic variation to future generations. The impacts of hatchery exposures are not only important to consider on the fish exposed, but also on future generations and evolutionary trajectory of fish in the river populations.

Effects on the Liver Transcriptome in Baltic Salmon: Contributions of Contamination with Organohalogen Compounds and Origin of Salmon

Hatchery-reared Atlantic salmon (Salmo salar) has been released to support the wild salmon stocks in the Baltic Sea for decades. During their feeding migration, salmon are exposed to organohalogen compounds (OHCs). Here, we investigated the OHC levels and transcriptome profiles in the liver of wild and hatchery-reared salmon collected from the Baltic main basin (BMB), the Bothnian Sea (BS), and the Gulf of Finland (GoF) and examined whether salmon origin and OHC levels contributed to the hepatic transcriptome profiles. There were no differences in the OHC concentrations between wild and reared fish but larger differences between areas. Several transcript levels were associated with non-dioxin-like polychlorinated biphenyls, polybrominated diphenylethers, chlordanes, and dichlorodiphenyltrichloroethane in a concentration-dependent manner. Between wild and reared salmon, lipid metabolism and related signaling pathways were enriched within the BMB and BS, while amino acid metabolism was altered within the GoF. When comparing the different areas, lipid metabolism, environmental stress and cell growth, and death-related pathways were enriched. Class coinertia analysis showed that the covariation in the OHC levels and the transcriptome were significantly similar. These results suggest that the hepatic transcriptomes in wild and hatchery-reared salmon are more affected by the OHC levels rather than the origin of salmon.

Major gene expression changes and epigenetic remodelling in Nile tilapia muscle after just one generation of domestication

The historically recent domestication of fishes has been essential to meet the protein demands of a growing human population. Selection for traits of interest during domestication is a complex process whose epigenetic basis is poorly understood. Cytosine hydroxymethylation is increasingly recognized as an important DNA modification involved in epigenetic regulation. In the present study, we investigated if hydroxymethylation plays a role in fish domestication and demonstrated for the first time at a genome-wide level and single nucleotide resolution that the muscle hydroxymethylome changes after a single generation of Nile tilapia (Oreochromis niloticus, Linnaeus) domestication. The overall decrease in hydroxymethylcytosine levels was accompanied by the downregulation of 2015 genes in fish reared in captivity compared to their wild progenitors. In contrast, several myogenic and metabolic genes that can affect growth potential were upregulated. There were 126 differentially hydroxymethylated cytosines between groups, which were not due to genetic variation; they were associated with genes involved in immune-, growth- and neuronal-related pathways. Taken together, our data unveil a new role for DNA hydroxymethylation in epigenetic regulation of fish domestication with impact in aquaculture and implications in artificial selection, environmental adaptation and genome evolution.

The Effects of Quantitative Trait Architecture on Detection Power in Short-Term Artificial Selection Experiments

Evolve and resequence (E&R) experiments, in which artificial selection is imposed on organisms in a controlled environment, are becoming an increasingly accessible tool for studying the genetic basis of adaptation. Previous work has assessed how different experimental design parameters affect the power to detect the quantitative trait loci (QTL) that underlie adaptive responses in such experiments, but so far there has been little exploration of how this power varies with the genetic architecture of the evolving traits. In this study, we use forward simulation to build a more realistic model of an E&R experiment in which a quantitative polygenic trait experiences a short, but strong, episode of truncation selection. We study the expected power for QTL detection in such an experiment and how this power is influenced by different aspects of trait architecture, including the number of QTL affecting the trait, their starting frequencies, effect sizes, clustering along a chromosome, dominance, and epistasis patterns. We show that all of these parameters can affect allele frequency dynamics at the QTL and linked loci in complex and often unintuitive ways, and thus influence our power to detect them. One consequence of this is that existing detection methods based on models of independent selective sweeps at individual QTL often have lower detection power than a simple measurement of allele frequency differences before and after selection. Our findings highlight the importance of taking trait architecture into account when designing and interpreting studies of molecular adaptation with temporal data. We provide a customizable modeling framework that will enable researchers to easily simulate E&R experiments with different trait architectures and parameters tuned to their specific study system, allowing for assessment of expected detection power and optimization of experimental design.

The power and limitations of gene expression pathway analyses toward predicting population response to environmental stressors

Rapid environmental changes impact the global distribution and abundance of species, highlighting the urgency to understand and predict how populations will respond. The analysis of differentially expressed genes has elucidated areas of the genome involved in adaptive divergence to past and present environmental change. Such studies however have been hampered by large numbers of differentially expressed genes and limited knowledge of how these genes work in conjunction with each other. Recent methods (broadly termed "pathway analyses") have emerged that aim to group genes that behave in a coordinated fashion to a factor of interest. These methods aid in functional annotation and uncovering biological pathways, thereby collapsing complex datasets into more manageable units, providing more nuanced understandings of both the organism-level effects of modified gene expression, and the targets of adaptive divergence. Here, we reanalyze a dataset that investigated temperature-induced changes in gene expression in marine-adapted and freshwater-adapted threespine stickleback (Gasterosteus aculeatus), using Weighted Gene Co-expression Network Analysis (WGCNA) with PANTHER Gene Ontology (GO)-Slim overrepresentation and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analysis. Six modules exhibited a conserved response and six a divergent response between marine and freshwater stickleback when acclimated to 7°C or 22°C. One divergent module showed freshwater-specific response to temperature, and the remaining divergent modules showed differences in height of reaction norms. PPARAa, a transcription factor that regulates fatty acid metabolism and has been implicated in adaptive divergence, was located in a module that had higher expression at 7°C and in freshwater stickleback. This updated methodology revealed patterns that were not found in the original publication. Although such methods hold promise toward predicting population response to environmental stressors, many limitations remain, particularly with regard to module expression representation, database resources, and cross-database integration.

A river runs through it: The causes, consequences, and management of intraspecific diversity in river networks

Rivers are fascinating ecosystems in which the eco-evolutionary dynamics of organisms are constrained by particular features, and biologists have developed a wealth of knowledge about freshwater biodiversity patterns. Over the last 10 years, our group used a holistic approach to contribute to this knowledge by focusing on the causes and consequences of intraspecific diversity in rivers. We conducted empirical works on temperate permanent rivers from southern France, and we broadened the scope of our findings using experiments, meta-analyses, and simulations. We demonstrated that intraspecific (genetic) diversity follows a spatial pattern (downstream increase in diversity) that is repeatable across taxa (from plants to vertebrates) and river systems. This pattern can result from interactive processes that we teased apart using appropriate simulation approaches. We further experimentally showed that intraspecific diversity matters for the functioning of river ecosystems. It indeed affects not only community dynamics, but also key ecosystem functions such as litter degradation. This means that losing intraspecific diversity in rivers can yield major ecological effects. Our work on the impact of multiple human stressors on intraspecific diversity revealed that-in the studied river systems-stocking of domestic (fish) strains strongly and consistently alters natural spatial patterns of diversity. It also highlighted the need for specific analytical tools to tease apart spurious from actual relationships in the wild. Finally, we developed original conservation strategies at the basin scale based on the systematic conservation planning framework that appeared pertinent for preserving intraspecific diversity in rivers. We identified several important research avenues that should further facilitate our understanding of patterns of local adaptation in rivers, the identification of processes sustaining intraspecific biodiversity-ecosystem function relationships, and the setting of reliable conservation plans.

Genomic signals found using RNA sequencing show signatures of selection and subtle population differentiation in walleye (Sander vitreus) in a large freshwater ecosystem

RNA sequencing is an effective approach for studying aquatic species yielding both physiological and genomic data. However, its population genetic applications are not well-characterized. We investigate this possible role for RNA sequencing for population genomics in Lake Winnipeg, Manitoba, Canada, walleye (Sander vitreus). Lake Winnipeg walleye represent the largest component of the second-largest freshwater fishery in Canada. In the present study, large female walleye were sampled via nonlethal gill biopsy over two years at three spawning sites representing a latitudinal gradient in the lake. Genetic variation from sequenced mRNA was analyzed for neutral and adaptive markers to investigate population structure and possible adaptive variation. We find low population divergence (F ST = 0.0095), possible northward gene flow, and outlier loci that vary latitudinally in transcripts associated with cell membrane proteins and cytoskeletal function. These results indicate that Lake Winnipeg walleye may be effectively managed as a single demographically connected metapopulation with contributing subpopulations and suggest genomic differences possibly underlying observed phenotypic differences. Despite its high cost relative to other genotyping methods, RNA sequencing data can yield physiological in addition to genetic information discussed here. We therefore argue that it is useful for addressing diverse molecular questions in the conservation of freshwater species.

Multiple decades of stocking has resulted in limited hatchery introgression in wild brook trout (Salvelinus fontinalis) populations of Nova Scotia

Many populations of freshwater fishes are threatened with losses, and increasingly, the release of hatchery individuals is one strategy being implemented to support wild populations. However, stocking of hatchery individuals may pose long-term threats to wild populations, particularly if genetic interactions occur between wild and hatchery individuals. One highly prized sport fish that has been heavily stocked throughout its range is the brook trout (Salvelinus fontinalis). In Nova Scotia, Canada, hatchery brook trout have been stocked since the early 1900s, and despite continued stocking efforts, populations have suffered declines in recent decades. Before this study, the genetic structure of brook trout populations in the province was unknown; however, given the potential negative consequences associated with hatchery stocking, it is possible that hatchery programs have adversely affected the genetic integrity of wild populations. To assess the influence of hatchery supplementation on wild populations, we genotyped wild brook trout from 12 river systems and hatchery brook trout from two major hatcheries using 100 microsatellite loci. Genetic analyses of wild trout revealed extensive population genetic structure among and within river systems and significant isolation-by-distance. Hatchery stocks were genetically distinct from wild populations, and most populations showed limited to no evidence of hatchery introgression (<5% hatchery ancestry). Only a single location had a substantial number of hatchery-derived trout and was located in the only river where a local strain is used for supplementation. The amount of hatchery stocking within a watershed did not influence the level of hatchery introgression. Neutral genetic structure of wild populations was influenced by geography with some influence of climate and stocking indices. Overall, our study suggests that long-term stocking has not significantly affected the genetic integrity of wild trout populations, highlighting the variable outcomes of stocking and the need to evaluate the consequences on a case-by-case basis.

Genetic consequences of pond production of a pikeperch (Sander lucioperca L.) stock with natural origin: the effects of changed selection pressure and reduced population size

The pikeperch (Sander lucioperca L.) possesses great potential for diversifying European aquaculture. However, studies on the genetic risk of stocking natural waters with farmed individuals of this species have been limited. Even the effect of pond culture on the genetic composition of stocks with natural-origin has not yet been determined. Our study aimed to compare the genetic variability of a wild living pikeperch population, a pond cultured broodstock (originating from the wild population) and its offspring generation. We also aimed to detect the potential signs of selection using three different methods. By analyzing the molecular data with 14 microsatellite markers, we illustrated that the impact of pond culture on the genetic diversity of fish stocks is similar to hatchery rearing due to its diversity reducing effect caused by using lower effective population sizes. Although the heterozygosity was similar in all populations (H_o = 0.68–0.71), the average number of alleles and allelic richness were significantly lower in the pond cultured stocks (NA = 7.5 and 6; AR = 7.5 and 5.9) compared to the wild population (NA = 11.00, AR = 10.47). Despite the semi-natural conditions of the present study, we detected changing selection pressure in one of the 14 microsatellite markers.