Host species and habitat shape fish-associated bacterial communities: phylosymbiosis between fish and their microbiome.

BACKGROUND: While many studies have reported that the structure of the gut and skin microbiota is driven by both species-specific and habitat-specific factors, the relative importance of host-specific versus environmental factors in wild vertebrates remains poorly understood. The aim of this study was to determine the diversity and composition of fish skin, gut, and surrounding water bacterial communities (hereafter referred to as microbiota) and assess the extent to which host habitat and phylogeny predict microbiota similarity. Skin swabs and gut samples from 334 fish belonging to 17 species were sampled in three Laurentian Great Lakes (LGLs) habitats (Detroit River, Lake Erie, Lake Ontario). We also collected and filtered water samples at the time of fish collection. We analyzed bacterial community composition using 16S metabarcoding and tested for community variation. RESULTS: We found that the water microbiota was distinct from the fish microbiota, although the skin microbiota more closely resembled the water microbiota. We also found that environmental (sample location), habitat, fish diet, and host species factors shape and promote divergence or convergence of the fish microbiota. Since host species significantly affected both gut and skin microbiota (separately from host species effects), we tested for phylosymbiosis using pairwise host species phylogenetic distance versus bacterial community dissimilarity. We found significant phylogenetic effects on bacterial community dissimilarity, consistent with phylosymbiosis for both the fish skin and gut microbiota, perhaps reflecting the longstanding co-evolutionary relationship between the host species and their microbiomes. CONCLUSIONS: Analyzing the gut and skin mucus microbiota across diverse fish species in complex natural ecosystems such as the LGLs provides insights into the potential for habitat and species-specific effects on the microbiome, and ultimately the health, of the host. Video Abstract.

Regulation of host gene expression by gastrointestinal tract microbiota in Chinook Salmon (Oncorhynchus tshawytscha).

Differences in gut microbiome composition are linked with health, disease and ultimately host fitness; however, the molecular mechanisms underlying that relationship are not well characterized. Here, we modified the fish gut microbiota using antibiotic and probiotic feed treatments to address the effect of host microbiome on gene expression patterns. Chinook salmon (Oncorhynchus tshawytscha) gut gene expression was evaluated using whole transcriptome sequencing (RNA-Seq) on hindgut mucosa samples from individuals treated with antibiotic, probiotic and control diets to determine differentially expressed (DE) host genes. Fifty DE host genes were selected for further characterization using nanofluidic qPCR chips. We used 16S rRNA gene metabarcoding to characterize the rearing water and host gut microbiome (bacterial) communities. Daily administration of antibiotics and probiotics resulted in significant changes in fish gut and aquatic microbiota as well as more than 100 DE genes in the antibiotic and probiotic treatment fish, relative to healthy controls. Normal microbiota depletion by antibiotics mostly led to downregulation of different aspects of immunity and upregulation of apoptotic process. In the probiotic treatment, genes related to post-translation modification and inflammatory responses were up-regulated relative to controls. Our qPCR results revealed significant effects of treatment (antibiotic and probiotic) on rabep2, aifm3, manf, prmt3 gene transcription. Moreover, we found significant associations between members of Lactobacillaceae and Bifidobacteriaceae with host gene expression patterns. Overall, our analysis showed that the microbiota had significant impacts on many host signalling pathways, specifically targeting immune, developmental and metabolic processes. Our characterization of some of the molecular mechanisms involved in microbiome-host interactions will help develop new strategies for preventing/ treating microbiome disruption-related diseases.

The effects of host quantitative genetic architecture on the gut microbiota composition of Chinook salmon (Oncorhynchus tshawytscha).

The microbiota consists of microbes living in or on an organism and has been implicated in host health and function. Environmental and host-related factors were shown to shape host microbiota composition and diversity in many fish species, but the role of host quantitative architecture across populations and among families within a population is not fully characterized. Here, Chinook salmon were used to determine if inter-population differences and additive genetic variation within populations influenced the gut microbiota diversity and composition. Specifically, hybrid stocks of Chinook salmon were created by crossing males from eight populations with eggs from an inbred line created from self-fertilized hermaphrodite salmon. Based on high-throughput sequencing of the 16S rRNA gene, significant gut microbial community diversity and composition differences were found among the hybrid stocks. Furthermore, additive genetic variance components varied among hybrid stocks, indicative of population-specific heritability patterns, suggesting the potential to select for specific gut microbiota composition for aquaculture purposes. Determining the role of host genetics in shaping their gut microbiota has important implications for predicting population responses to environmental changes and will thus impact conservation efforts for declining populations of Chinook salmon.

Variation in the diversity of bacterial communities and potential sources of fecal contamination of beaches in the Huron to Erie corridor.

Understanding the diversity of bacteria and E.coli levels at beaches is important for managing health risks. This study compared temporal changes of the bacterial communities of Belle Isle Beach (Detroit, MI) and Sand Point Beach (Windsor, ONT), both located near the Lake St. Clair origin of the Detroit River. Water samples collected 4 days/week for 12 weeks in summer, were subjected to 16S rRNA analysis of amplicon sequencing and E. coli enumeration. Bacterial communities changed over time, as determined by cluster dendrogram analysis, exhibiting different communities in July and August than in June and different communities at the two beaches. After June, alpha diversity decreased and relative abundance of Enterobacter (Gammaproteobacteria) increased at Sand Point; whereas, Belle Isle maintained its alpha diversity and dominance by Betaproteobacteria and Actinobacteria. Contamination at both beaches is dominated by birds (23% to 50% of samples), while only ∼10% had evidence of human-associated bacteria. High E. coli at both beaches was often associated with precipitation. Nearshore sampling counts were higher than waist-deep sampling counts. Despite the dynamic changes in bacterial communities between the two beaches, this analysis based on 16S rRNA amplicon sequencing is able to provide information about bacterial types associated with high E. coli levels and to use bacterial sequences to more precisely determine sources and health relevance of contaminants.

Acute thermal stress elicits interactions between gene expression and alternative splicing in a fish of conservation concern.

Transcriptomic research provides a mechanistic understanding of an organism's response to environmental challenges such as increasing temperatures, which can provide key insights into the threats posed by thermal challenges associated with urbanization and climate change. Differential gene expression and alternative splicing are two elements of the transcriptomic stress response that may work in tandem, but relatively few studies have investigated these interactions in fishes of conservation concern. We studied the imperilled redside dace (Clinostomus elongatus) as thermal stress is hypothesized to be an important cause of population declines. We tested the hypothesis that gene expression-splicing interactions contribute to the thermal stress response. Wild fish exposed to acute thermal stress were compared with both handling controls and fish sampled directly from a river. Liver tissue was sampled to study the transcriptomic stress response. With a gene set enrichment analysis, we found that thermally stressed fish showed a transcriptional response related to transcription regulation and responses to unfolded proteins, and alternatively spliced genes related to gene expression regulation and metabolism. One splicing factor, prpf38b, was upregulated in the thermally stressed group compared with the other treatments. This splicing factor may have a role in the Jun/AP-1 cellular stress response, a pathway with wide-ranging and context-dependent effects. Given large gene interaction networks and the context-dependent nature of transcriptional responses, our results highlight the importance of understanding interactions between gene expression and splicing for understanding transcriptomic responses to thermal stress. Our results also reveal transcriptional pathways that can inform conservation breeding, translocation and reintroduction programs for redside dace and other imperilled species by identifying appropriate source populations.

Conservation implications of revised genetic structure resulting from new population discovery: the threatened eastern sand darter (Ammocrypta pellucida) in Canada.

Human activity can put non-game fishes at higher risk of extinction because of inappropriate management action. Eastern sand darter (Ammocrypta pellucida), a small benthic fish classified as threatened across much of its northern range, inhabits increasingly fragmented sandy habitats and, as a non-game fish, may be easily overlooked in conservation efforts. In this study, the authors use genotype data from nine microsatellite loci and cytochrome oxidase I (COI) sequencing data across its northern native range to re-assess genetic structure and to characterize a newly discovered, geographically disjunct population. Previous microsatellite marker analyses had identified seven distinct population genetic clusters across the region sampled; the analysis of this study showed that the newly discovered population (West Lake, Ontario) exhibits a divergent structure. COI haplotype analysis suggests that a single haplotype recolonized the Great Lakes and surrounding water bodies after the Wisconsinan glacial period, and subsequent fluctuation in water levels and habitat fragmentation resulted in divergence of genetic clusters. Although the novel West Lake population has a common ancestral source with other populations in the broader region, its divergent genetic signature merits its consideration as a separate conservation unit. The analyses of this study highlight the potential conservation implications of the discovery of new populations, particularly those of at-risk species, even within a region that has been genetically well characterized.

Averting an Outbreak of SARS-CoV-2 in a University Residence Hall through Wastewater Surveillance.

A wastewater surveillance program targeting a university residence hall was implemented during the spring semester 2021 as a proactive measure to avoid an outbreak of COVID-19 on campus. Over a period of 7 weeks from early February through late March 2021, wastewater originating from the residence hall was collected as grab samples 3 times per week. During this time, there was no detection of SARS-CoV-2 by reverse transcriptase quantitative PCR (RT-qPCR) in the residence hall wastewater stream. Aiming to obtain a sample more representative of the residence hall community, a decision was made to use passive samplers beginning in late March onwards. Adopting a Moore swab approach, SARS-CoV-2 was detected in wastewater samples just 2 days after passive samplers were deployed. These samples also tested positive for the B.1.1.7 (Alpha) variant of concern (VOC) using RT-qPCR. The positive result triggered a public health case-finding response, including a mobile testing unit deployed to the residence hall the following day, with testing of nearly 200 students and staff, which identified two laboratory-confirmed cases of Alpha variant COVID-19. These individuals were relocated to a separate quarantine facility, averting an outbreak on campus. Aggregating wastewater and clinical data, the campus wastewater surveillance program has yielded the first estimates of fecal shedding rates of the Alpha VOC of SARS-CoV-2 in individuals from a nonclinical setting. IMPORTANCE Among early adopters of wastewater monitoring for SARS-CoV-2 have been colleges and universities throughout North America, many of whom are using this approach to monitor congregate living facilities for early evidence of COVID-19 infection as an integral component of campus screening programs. Yet, while there have been numerous examples where wastewater monitoring on a university campus has detected evidence for infection among community members, there are few examples where this monitoring triggered a public health response that may have averted an actual outbreak. This report details a wastewater-testing program targeting a residence hall on a university campus during spring 2021, when there was mounting concern globally over the emergence of SARS-CoV-2 variants of concern, reported to be more transmissible than the wild-type Wuhan strain. In this communication, we present a clear example of how wastewater monitoring resulted in actionable responses by university administration and public health, which averted an outbreak of COVID-19 on a university campus.

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.

Spatial and temporal genetic variation in an exploited reef fish: The effects of exploitation on cohort genetic structure.

Many coral reef fishes are fished, often resulting in detrimental genetic effects; however, reef fishes often show unpredictable patterns of genetic variation, which potentially mask the effects of fishing. Our goals were to characterize spatial and temporal genetic variation and determine the effects of fishing on an exploited reef fish, Plectropomus leopardus, Lacepède (the common coral trout). To determine population structure, we genotyped 417 Great Barrier Reef coral trout from four populations sampled in 2 years (1996 and 2004) at nine microsatellite loci. To test for exploitation effects, we additionally genotyped 869 individuals from a single cohort (ages 3-5) across eight different reefs, including fished and control populations. Genetic structure differed substantially in the two sampled years, with only 1 year exhibiting isolation by distance. Thus, genetic drift likely plays a role in shaping population genetic structure in this species. Although we found no loss of genetic diversity associated with exploitation, our relatedness patterns show that pulse fishing likely affects population genetics. Additionally, genetic structure in the cohort samples likely reflected spatial variation in recruitment contributing to genetic structure at the population level. Overall, we show that fishing does impact coral reef fishes, highlighting the importance of repeated widespread sampling to accurately characterize the genetic structure of reef fishes, as well as the power of analysing cohorts to avoid the impacts of recruitment-related genetic swamping. The high temporal and spatial variability in genetic structure, combined with possible selection effects, will make conservation/management of reef fish species complex.

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.

Microbial community and abiotic effects on aquatic bacterial communities in north temperate lakes.

The aquatic bacterial community (BC) plays a vital role in determining the nature and rate of ecosystem function. However, the biotic and abiotic factors influencing BC structure and function are largely unknown. Hence, the current study characterizes the impact of biotic and abiotic factors on aquatic bacterial biodiversity to determine whether the dominant effects are biotic or abiotic by partitioning their relative effects across temperate Canadian lakes. We collected water samples from sixty southern Ontario lakes and characterized their BC and microbial eukaryotic community (MEC) compositions using high throughput metabarcode sequencing of 16S and 18S rRNA gene fragments. The diversity and richness of aquatic BCs differed considerably among our study lakes, and those differences were explained by environmental, spatial, and biotic (MEC) factors (31%, 23%, and 23% of variance explained, respectively). The relatively large contribution from biotic and abiotic factors (54%), relative to spatial effects, shows deterministic processes prevail in shaping BC assembly in freshwater lakes. However, spatial effects also contributed significantly, highlighting the role of stochastic processes (ecological drift and coupled with limited dispersal) in shaping BC structure. Furthermore, our co-occurrence network analysis showed strong positive and negative interactions within and between the BCs and MECs, indicating mutualistic or antagonistic co-occurrence patterns relationships play important roles in driving the variation in BC composition among our sampled lakes. Considered together, our community analyses show that deterministic and stochastic processes combined contribute to determining the aquatic BC composition, and hence likely function as well, across a broad array of temperate freshwater lakes.

Diel Dynamics of Freshwater Bacterial Communities at Beaches in Lake Erie and Lake St. Clair, Windsor, Ontario.

Bacteria play a key role in freshwater biogeochemical cycling as well as water safety, but short-term trends in freshwater bacterial community composition and dynamics are not yet well characterized. We sampled four public beaches in southern Ontario, Canada; in June, July, and August (2016) over a 24-h (diel) cycle at 2-h intervals. Using high-throughput sequencing of 16S rRNA gene, we found substantial bi-hourly and day/night variation in the bacterial communities with considerable fluctuation in the relative abundance of Actinobacteria and Proteobacteria phyla. Moreover, relative abundance of Enterobacteriaceae (associated with potential health risk) was significantly high at night in some dial cycles. Diversity was significantly high at night across most of the diel sampling events. qPCR assays showed a substantial bi-hourly variation of Escherichia coli levels with a significant high level of E. coli at night hours in comparison with day hours and the lowest levels at noon and during the afternoon hours. Taken together, these findings highlighted a considerable short-term temporal variation of bacterial communities which helps better understanding of freshwater bacterial dynamics and their ecology. E. coli monitoring showed that multiple samples in different hours will provide more accurate picture of freshwater safety and human health risk. Graphical abstract.

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.

Spatial and environmental effects on Coho Salmon life history trait variation.

Adult size, egg size, fecundity, and mass of gonads are affected by trade-offs between reproductive investment and environmental conditions shaping the evolution of life history traits among populations for widely distributed species. Coho salmon Oncorhynchus kisutch have a large geographic distribution, and different environmental conditions are experienced by populations throughout their range. We examined the effect of environmental variables on female size, egg size, fecundity, and reproductive investment of populations of Coho Salmon from across British Columbia using an information theoretic approach. Female size increased with latitude and decreased with migration distance from the ocean to spawning locations. Egg size was lowest for intermediate intragravel temperature during incubation, decreased with migration distance, but increased in rivers below lakes. Fecundity increased with latitude, warmer temperature during the spawning period, and river size, but decreased in rivers below lakes compared with rivers with tributary sources. Relative gonad size increased with latitude and decreased with migration distance. Latitude of spawning grounds, migratory distance, and temperatures experienced by a population, but also hydrologic features-river size and headwater source-are influential in shaping patterns of reproductive investment, particularly egg size. Although, relative gonad size varied with latitude and migration distance, how gonadal mass was partitioned gives insight into the trade-off between egg size and fecundity. The lack of an effect of latitude on egg size suggests that local optima for egg size related to intragravel temperature may drive the variation in fecundity observed among years.

Metabarcoding of native and invasive species in stomach contents of Great Lakes fishes.

As aquatic invasive species (AIS) proliferate worldwide, a better understanding of their roles in invaded habitats is needed to inform management and introduction prevention strategies and priorities. Metabarcoding of stomach content DNA (scDNA) shows considerable promise in such regard. We thus metabarcoded scDNA from two non-native fish species (alewife (Alosa pseudoharengus) and rainbow smelt (Osmerus mordax)), and three native ones (bloater (Coregonus hoyi), ninespine stickleback (Pungitius pungitius), and slimy sculpin (Cottus cognatus)). Fishes (N = 376) were sampled in spring 2009 and 2010 from 73-128 m depths at three Lake Michigan sites. Four mitochondrial cytochrome oxidase 1 (CO1) primer sets designed to target five potential AIS prey, and a universal aquatic invertebrate CO1 primer set targeting both native and AIS prey were used. Quality controlled prey amplicons were matched to three AIS prey: Bythotrephes longimanus (mean percent frequency occurrence, all samples = 7%), Cercopagis pengoi (5%), and Dreissena rostriformis bugensis (11%). Neither invasive prey Dreissena polymorpha nor Hemimysis anomala were detected. Native prey Leptodiaptomus sicilis, Limnocalanus macrurus, and Mysis diluviana were relatively common in scDNA (respective mean percent occurrences, all samples: 48%, 25%, 42%). Analysis of variation in prey occurrences for sample site, predator species, sample year, sample depth, and predator total length (TL) indicated site and predator species were most important. However, B. longimanus occurrence in scDNA depended upon predator TL, perhaps indicative of its unique defensive spine limiting susceptibility to predation until fishes exceed species-specific gape-based limitations. Our analysis of native and invasive prey species indicated possible indirect AIS impacts such as native predators switching their diet due to AIS-driven losses of preferred native prey. Metabarcoding demonstrated that AIS are integrated components of the offshore Lake Michigan food web, with both native and non-native predators, and both invasive and native prey are affecting species interactions across multiple trophic levels.

Loci associated with variation in gene expression and growth in juvenile salmon are influenced by the presence of a growth hormone transgene.

BACKGROUND: Growth regulation is a complex process influenced by genetic and environmental factors. We examined differences between growth hormone (GH) transgenic (T) and non-transgenic (NT) coho salmon to elucidate whether the same loci were involved in controlling body size and gene expression phenotypes, and to assess whether physiological transformations occurring from GH transgenesis were under the influence of alternative pathways. The following genomic techniques were used to explore differences between size classes within and between transgenotypes (T vs. NT): RNA-Seq/Differentially Expressed Gene (DEG) analysis, quantitative PCR (qPCR) and OpenArray analysis, Genotyping-by-Sequencing, and Genome-Wide Association Study (GWAS). RESULTS: DEGs identified in comparisons between the large and small tails of the size distributions of T and NT salmon (NTLarge, NTSmall, TLarge and TSmall) spanned a broad range of biological processes, indicating wide-spread influence of the transgene on gene expression. Overexpression of growth hormone led to differences in regulatory loci between transgenotypes and size classes. Expression levels were significantly greater in T fish at 16 of 31 loci and in NT fish for 10 loci. Eleven genes exhibited different mRNA levels when the interaction of size and transgenotype was considered (IGF1, IGFBP1, GH, C3-4, FAS, FAD6, GLUT1, G6PASE1, GOGAT, MID1IP1). In the GWAS, 649 unique SNPs were significantly associated with at least one study trait, with most SNPs associated with one of the following traits: C3_4, ELA1, GLK, IGF1, IGFBP1, IGFII, or LEPTIN. Only 1 phenotype-associated SNP was found in common between T and NT fish, and there were no SNPs in common between transgenotypes when size was considered. CONCLUSIONS: Multiple regulatory loci affecting gene expression were shared between fast-growing and slow-growing fish within T or NT groups, but no such regulatory loci were found to be shared between NT and T groups. These data reveal how GH overexpression affects the regulatory responses of the genome resulting in differences in growth, physiological pathways, and gene expression in T fish compared with the wild type. Understanding the complexity of regulatory gene interactions to generate phenotypes has importance in multiple fields ranging from applications in selective breeding to quantifying influences on evolutionary processes.

DNA Methylation Profiles Suggest Intergenerational Transfer of Maternal Effects.

The view of maternal effects (nongenetic maternal environmental influence on offspring phenotype) has changed from one of distracting complications in evolutionary genetics to an important evolutionary mechanism for improving offspring fitness. Recent studies have shown that maternal effects act as an adaptive mechanism to prepare offspring for stressful environments. Although research into the magnitude of maternal effects is abundant, the molecular mechanisms of maternal influences on offspring phenotypic variation are not fully understood. Despite recent work identifying DNA methylation as a potential mechanism of nongenetic inheritance, currently proposed links between DNA methylation and parental effects are indirect and primarily involve genomic imprinting. We combined a factorial breeding design and gene-targeted sequencing methods to assess inheritance of methylation during early life stages at 14 genes involved in growth, development, metabolism, stress response, and immune function of Chinook salmon (Oncorhynchus tshawytscha). We found little evidence for additive or nonadditive genetic effects acting on methylation levels during early development; however, we detected significant maternal effects. Consistent with conventional maternal effect data, maternal effects on methylation declined through development and were replaced with nonadditive effects when offspring began exogenous feeding. We mapped methylation at individual CpG sites across the selected candidate genes to test for variation in site-specific methylation profiles and found significant maternal effects at selected CpG sites that also declined with development stage. While intergenerational inheritance of methylated DNA is controversial, we show that CpG-specific methylation may function as an underlying molecular mechanism for maternal effects, with important implications for offspring fitness.

Transcriptional Basis of Copper-Induced Olfactory Impairment in the Sea Lamprey, a Primitive Invasive Fish.

Olfaction mediates behaviors necessary for survival and reproduction in fishes. Anthropogenic inputs of contaminants into aquatic environments, specifically copper, are known to disrupt a broad range of olfactory-mediated behaviors and can cause long-lasting damage even at low concentrations that have profound impacts on the biology of aquatic organisms. The sea lamprey (Petromyzon marinus) is a primitive fish species invasive to the North American Great Lakes that relies on olfaction to navigate during natal homing and in mate choice during reproduction. To investigate effects of copper on sea lamprey olfaction and the potential for maintenance of olfactory function during copper exposure, we exposed juvenile sea lamprey to environmentally ecologically relevant copper concentrations (0, 5, 10 and 30 µg/L) for 24 hr and characterized gene transcription response in olfactory tissue (i.e., peripheral olfactory organ and olfactory bulb) and forebrain using whole transcriptome sequencing. Copper exposure induced a pattern of positive dose-dependent transcriptional response. Expression changes primarily reflected up-regulation of genes involved in apoptosis and wound healing. Unlike higher vertebrates, genes specifically related to the olfactory senses of the sea lamprey, e.g., olfactory receptors, exhibited little transcriptional response to copper exposure, suggesting the mechanism of copper-induced olfactory impairment is through necrosis of the olfactory bulb and not copper-selective inhibition of olfactory receptors. Fully two-thirds of the differentially expressed genes at higher doses of copper have no known function and thus represent important candidates for further study of the responses to copper-induced olfactory injury. Our results shed light on the evolution of vertebrate olfactory repair mechanisms and have important implications for the conservation and management of both invasive and native populations of lamprey.

Variation in juvenile Chinook salmon (Oncorhynchus tshawytscha) transcription profiles among and within eight population crosses from British Columbia, Canada.

Phenotypic differences among populations within a species have been reported for a variety of traits, ranging from life history to physiology to gene transcription. Population-level phenotypic variation has been attributed to genetic differences resulting from genetic drift and/or local adaptation as well as environmental differences resulting from plasticity. We studied population- and family-level variation in gene transcription for 22 fitness-related genes, comprising immune, growth, metabolic, and stress processes in Chinook salmon (Oncorhynchus tshawytscha). We created hybrid Chinook salmon families from eight populations and treated them with an immune stimulus, a handling stress challenge, and held some as a no-treatment control group. Population effects, sire effects, and narrow-sense heritability (h2 ) were calculated for each candidate gene within each treatment group. We expected population to have a significant effect on gene transcription for many of our genes; however, we found a population effect for transcription at only one immune gene at rest. The limited number of significant population effects on gene transcription, combined with significant additive genetic variance within each population does not support the expectation of past strong selection pressures acting on heritable transcription profiles among populations. Instead, our results indicate that Chinook salmon likely adapt to their local environment through transcriptional plasticity rather than fixed differences. The expectation for fixed population-level differences in gene transcription at fitness-related genes, reflecting accepted models of local adaptation is high; however, comparisons among multiple populations using half-sibling breeding designs are rare. Our work fills an important gap in our growing understanding of the process of among and within-population divergence.

Significant differences in maternal carotenoid provisioning and effects on offspring fitness in Chinook salmon colour morphs.

In oviparous species, maternal carotenoid provisioning can deliver diverse fitness benefits to offspring via increased survival, growth and immune function. Despite demonstrated advantages of carotenoids, large intra- and interspecific variation in carotenoid utilization exists, suggesting trade-offs associated with carotenoids. In Chinook salmon (Oncorhynchus tshawytscha), extreme variation in carotenoid utilization delineates two colour morphs (red and white) that differ genetically in their ability to deposit carotenoids into tissues. Here, we take advantage of this natural variation to examine how large differences in maternal carotenoid provisioning influence offspring fitness. Using a full factorial breeding design crossing morphs and common-garden rearing, we measured differences in a suite of fitness-related traits, including survival, growth, viral susceptibility and host response, in offspring of red (carotenoid-rich eggs) and white (carotenoid-poor eggs) females. Eggs of red females had significantly higher carotenoid content than those of white females (6× more); however, this did not translate into measurable differences in offspring fitness. Given that white Chinook salmon may have evolved to counteract their maternal carotenoid deficiency, we also examined the relationship between egg carotenoid content and offspring fitness within each morph separately. Egg carotenoids only had a positive effect within the red morph on survival to eyed-egg (earliest measured trait), but not within the white morph. Although previous work shows that white females benefit from reduced egg predation, our study also supports a hypothesis that white Chinook salmon have evolved additional mechanisms to improve egg survival despite low carotenoids, providing novel insight into evolutionary mechanisms that maintain this stable polymorphism.