The genome and transcriptome of the snail Biomphalaria sudanica s.l.: immune gene diversification and highly polymorphic genomic regions in an important African vector of Schistosoma mansoni.

BACKGROUND: Control and elimination of schistosomiasis is an arduous task, with current strategies proving inadequate to break transmission. Exploration of genetic approaches to interrupt Schistosoma mansoni transmission, the causative agent for human intestinal schistosomiasis in sub-Saharan Africa and South America, has led to genomic research of the snail vector hosts of the genus Biomphalaria. Few complete genomic resources exist, with African Biomphalaria species being particularly underrepresented despite this being where the majority of S. mansoni infections occur. Here we generate and annotate the first genome assembly of Biomphalaria sudanica sensu lato, a species responsible for S. mansoni transmission in lake and marsh habitats of the African Rift Valley. Supported by whole-genome diversity data among five inbred lines, we describe orthologs of immune-relevant gene regions in the South American vector B. glabrata and present a bioinformatic pipeline to identify candidate novel pathogen recognition receptors (PRRs). RESULTS: De novo genome and transcriptome assembly of inbred B. sudanica originating from the shoreline of Lake Victoria (Kisumu, Kenya) resulted in a haploid genome size of ~ 944.2 Mb (6,728 fragments, N50 = 1.067 Mb), comprising 23,598 genes (BUSCO = 93.6% complete). The B. sudanica genome contains orthologues to all described immune genes/regions tied to protection against S. mansoni in B. glabrata, including the polymorphic transmembrane clusters (PTC1 and PTC2), RADres, and other loci. The B. sudanica PTC2 candidate immune genomic region contained many PRR-like genes across a much wider genomic region than has been shown in B. glabrata, as well as a large inversion between species. High levels of intra-species nucleotide diversity were seen in PTC2, as well as in regions linked to PTC1 and RADres orthologues. Immune related and putative PRR gene families were significantly over-represented in the sub-set of B. sudanica genes determined as hyperdiverse, including high extracellular diversity in transmembrane genes, which could be under pathogen-mediated balancing selection. However, no overall expansion in immunity related genes was seen in African compared to South American lineages. CONCLUSIONS: The B. sudanica genome and analyses presented here will facilitate future research in vector immune defense mechanisms against pathogens. This genomic/transcriptomic resource provides necessary data for the future development of molecular snail vector control/surveillance tools, facilitating schistosome transmission interruption mechanisms in Africa.

The genome and transcriptome of the snail Biomphalaria sudanica s.l.: Immune gene diversification and highly polymorphic genomic regions in an important African vector of Schistosoma mansoni.

Background: Control and elimination of schistosomiasis is an arduous task, with current strategies proving inadequate to break transmission. Exploration of genetic approaches to interrupt Schistosoma mansoni transmission, the causative agent for human intestinal schistosomiasis in sub-Saharan Africa and South America, has led to genomic research of the snail vector hosts of the genus Biomphalaria. Few complete genomic resources exist, with African Biomphalaria species being particularly underrepresented despite this being where the majority of S. mansoni infections occur. Here we generate and annotate the first genome assembly of Biomphalaria sudanica sensu lato, a species responsible for S. mansoni transmission in lake and marsh habitats of the African Rift Valley. Supported by whole-genome diversity data among five inbred lines, we describe orthologs of immune-relevant gene regions in the South American vector B. glabrata and present a bioinformatic pipeline to identify candidate novel pathogen recognition receptors (PRRs). Results: De novo genome and transcriptome assembly of inbred B. sudanica originating from the shoreline of Lake Victoria (Kisumu, Kenya) resulted in a haploid genome size of ~944.2 Mb (6732 fragments, N50=1.067 Mb), comprising 23,598 genes (BUSCO=93.6% complete). The B. sudanica genome contains orthologues to all described immune genes/regions tied to protection against S. mansoni in B. glabrata. The B. sudanica PTC2 candidate immune genomic region contained many PRR-like genes across a much wider genomic region than has been shown in B. glabrata, as well as a large inversion between species. High levels of intra-species nucleotide diversity were seen in PTC2, as well as in regions linked to PTC1 and RADres orthologues. Immune related and putative PRR gene families were significantly over-represented in the sub-set of B. sudanica genes determined as hyperdiverse, including high extracellular diversity in transmembrane genes, which could be under pathogen-mediated balancing selection. However, no overall expansion in immunity related genes were seen in African compared to South American lineages. Conclusions: The B. sudanica genome and analyses presented here will facilitate future research in vector immune defense mechanisms against pathogens. This genomic/transcriptomic resource provides necessary data for the future development of molecular snail vector control/surveillance tools, facilitating schistosome transmission interruption mechanisms in Africa.

Molecular variability of the Ancylostoma secreted Protein-2 (Aca-asp-2) gene from Ancylostoma caninum contributes to expand information on population genetic studies of hookworms.

Hookworm infection is a major public health problem in many regions of the world. Given the high levels of host morbidity and even mortality of the host caused by these infections, it is crucial to understand the genetic structure of hookworm populations. This understanding can provide insights into the ecology, transmission patterns, mechanisms of drug resistance, and the development of vaccines and immunotherapeutic strategies. Previously, we examined presumably neutral molecular markers, such as microsatellites and COI (Cytochrome C oxidase subunit 1) in Brazilian populations of Ancylostoma caninum. Here we analyze the molecular variability of a genomic fragment of the Aca-asp-2 (Ancylostoma secreted protein-2) gene from Ancylostoma caninum. This gene is a highly expressed and activated following the infection of the L3 larvae in the host. We obtained individuals of A. caninum from five different geographic locations in Brazil, sequenced and analyzed parts of the gene. The results revealed extensive polymorphism at this fragment, especially in the intronic region, indicating low selective pressure acting on these sequences. However, we also observed irregular distributions of nucleotides and polymorphisms in the coding region of this gene, resulting in the identification of 27 alleles. The data presented here contribute to expanding the understanding of population genetic studies of hookworms.

PTC2 region genotypes counteract Biomphalaria glabrata population differences between M-line and BS90 in resistance to infection by Schistosoma mansoni.

Background: Biomphalaria glabrata is a snail intermediate host for Schistosoma mansoni, a trematode responsible for human schistosomiasis. BS90 is one of the most well studied strains of B. glabrata owing to its high resistance to infection by most strains of S. mansoni. An F2 mapping study from 1999 identified two RAPD markers that associated with what appeared to be single-locus, dominant resistance by the BS90 population relative to the susceptible M-line population. One marker cannot be mapped, but the other, OPM-04, maps to within 5 Mb of PTC2, a region we recently showed has a very large effect on resistance within another snail population challenged by the same strain of parasite (PR1). Here we tested the hypothesis that the PTC2 region contains the causal gene/s that explain the iconic resistance of BS90 snails. Methods: We used marker-assisted backcrossing to drive the BS90 version of the PTC2 region (+/-~1 Mb on either side) into an M-line (susceptible strain) genetic background, and the M-line version into a BS90 genetic background. We challenged the offspring with PR1-strain schistosomes and tested for effects of allelic variation in the PTC2 region in a common genetic background. Results: Relative to M-line haplotypes, the BS90 haplotype actually confers enhanced susceptibility. So we reject our original hypothesis. One possible explanation for our result was that the causal gene linked to OPM-04 is near, but not in the PTC2 block that we introgressed into each line. So we used an F2 cross to independently test the effects of the PTC2 and OPM-04 regions in a randomized genetic background. We confirmed that the BS90 haplotype confers increased susceptibility, and we see a similar, although non-significant effect at OPM-04. We discuss possible reasons why our results differed so dramatically from those of the 1999 study. We also present Pacbio assemblies of the PTC2 and flanking region in BS90 and M-line, compare with previously published PTC2 haplotypes, and discuss candidate genes that might be behind the enhanced susceptibility of the BS90 haplotype.

Offspring of first-generation hatchery steelhead trout (Oncorhynchus mykiss) grow faster in the hatchery than offspring of wild fish, but survive worse in the wild: Possible mechanisms for inadvertent domestication and fitness loss in hatchery salmon.

Salmonid fish raised in hatcheries often have lower fitness (number of returning adult offspring) than wild fish when both spawn in the wild. Body size at release from hatcheries is positively correlated with survival at sea. So one explanation for reduced fitness is that hatcheries inadvertently select for trait values that enhance growth rate under the unnatural environment of a hatchery, but that are maladaptive in the wild environment. A simple prediction of this hypothesis is that juveniles of hatchery origin should grow more quickly than fish of wild origin under hatchery conditions, but should have lower survival under wild conditions. We tested that hypothesis using multiple full sibling families of steelhead (Oncorhynchus mykiss) that were spawned using either two wild parents (WxW) or two first-generation hatchery (HxH) parents. Offspring from all the families were grown together under hatchery conditions and under semi-natural conditions in artificial streams. HxH families grew significantly faster in the hatchery, but had significantly lower survival in the streams. That we see this tradeoff after only a single generation of selection suggests that the traits involved are under very strong selection. We also considered one possible alteration to the hatchery environment that might reduce the intensity of selection among families in size at release. Here we tested whether reducing the fat content of hatchery feed would reduce the variance among families in body size. Although fish raised under a low-fat diet were slightly smaller, the variation among families in final size was unchanged. Thus, there is no evidence that reducing the fat content of hatchery feed would reduce the opportunity for selection among families on size at release.

Clusters of polymorphic transmembrane genes control resistance to schistosomes in snail vectors.

Schistosomiasis is a debilitating parasitic disease infecting hundreds of millions of people. Schistosomes use aquatic snails as intermediate hosts. A promising avenue for disease control involves leveraging innate host mechanisms to reduce snail vectorial capacity. In a genome-wide association study of Biomphalaria glabrata snails, we identify genomic region PTC2 which exhibits the largest known correlation with susceptibility to parasite infection (>15 fold effect). Using new genome assemblies with substantially higher contiguity than the Biomphalaria reference genome, we show that PTC2 haplotypes are exceptionally divergent in structure and sequence. This variation includes multi-kilobase indels containing entire genes, and orthologs for which most amino acid residues are polymorphic. RNA-Seq annotation reveals that most of these genes encode single-pass transmembrane proteins, as seen in another resistance region in the same species. Such groups of hyperdiverse snail proteins may mediate host-parasite interaction at the cell surface, offering promising targets for blocking the transmission of schistosomiasis.

Schistosomiasis is a widespread parasitic disease, affecting over 200 million people in tropical countries. It is caused by schistosome worms, which are carried by freshwater snails. These snails release worm larvae into the water, where they can infect humans ­ for example, after bathing or swimming. Treatment options for schistosomiasis are limited. Eliminating the freshwater snails is one way to control the disease, but this is not always effective in the long term and the chemicals used can also harm other animals in the water. Another way to manage schistosomiasis could be to stop the worms from infecting their snail host by breaking the parasites' life cycle without killing the snails. It is already known that some snails are naturally resistant to infection by some strains of schistosomes. Since this immunity is also inherited by the offspring of resistant snails, there is likely a genetic mechanism behind it. However, very little else is known about any genes that might be involved. Tennessen et al. therefore set out to identify what genes were responsible for schistosome resistance and how they worked. The experiments used a large laboratory colony of snails, whose susceptibility to schistosome infection varied among individual animals. To determine the genes behind this variation, Tennessen et al. first searched for areas of DNA that also differed between the immune and infected snails. Comparing genetic sequences across over 1,000 snails revealed a distinct region of DNA that had a large effect on how likely they were to be infected. This section of DNA turned out to be highly diverse, with different snails carrying varying numbers and different forms of the genes within this region. Many of these genes appear to encode proteins found on the surface of snail cells, which could affect whether snails and worms can recognize each other when they come into contact. This in turn could determine whether or not the worms can infect their hosts. These results shed new light on how the snails that carry schistosomes may be able to resist infections. In the future, this knowledge could be key to controlling schistosomiasis, either by releasing genetically engineered, immune snails into the wild (thus making it harder for the parasites to reproduce) or by using the snails' mechanism of resistance to design better drug therapies.

Heat shock increases hydrogen peroxide release from circulating hemocytes of the snail Biomphalaria glabrata.

Planorbid freshwater snails are important intermediate hosts for parasitic diseases caused by parasitic worms, most notably schistosomiasis. There are numerous reports of snails, specifically Biomphalaria glabrata, having compromised defences against schistosomes after being exposed to thermal stress. Environmental modifications to the defenses of schistosome transmitting snails could have negative ramifications for human disease risk in the context of climate change. Here the effects of heat shock on the production of hydrogen peroxide, a primary anti-microbial effector in many molluscs, were examined. The present findings show that heat shock increases NADPH oxidase 2 mRNA levels and hydrogen peroxide produced by snail hemocytes, and that both of these phenotypes could be reversed by an HSP-90 inhibitor. These findings indicate that snail defense systems are altered by heat shock at a molecular level in B. glabrata, and that snail immunity to many pathogens may be altered by the rapid variations in temperature that are associated with global climate change.

Neither heat pulse, nor multigenerational exposure to a modest increase in water temperature, alters the susceptibility of Guadeloupean Biomphalaria glabrata to Schistosoma mansoni infection.

There are increasing concerns regarding the role global climate change will have on many vector-borne diseases. Both mathematical models and laboratory experiments suggest that schistosomiasis risk may change as a result of the effects of increasing temperatures on the planorbid snails that host schistosomes. Heat pulse/heat shock of the BS90 strain of Biomphalaria glabrata was shown to increase the rate of infection by Schistosoma mansoni, but the result was not replicable in a follow up experiment by a different lab. We characterised the susceptibility and cercarial shedding of Guadeloupean B. glabrata after infection with S. mansoni under two temperature regimes: multigenerational exposure to small increases in temperature, and extreme heat pulse events. Neither long-term, multigenerational rearing at elevated temperatures, nor transient heat pulse modified the susceptibility of Guadeloupean B. glabrata to infection (prevalence) or shedding of schistosome cercaria (intensity of infection). These findings suggest that heat pulse-induced susceptibility in snail hosts may be dependent on the strain of the snail and/or schistosome, or on some as-yet unidentified environmental co-factor.

Allelic variation in a single genomic region alters the hemolymph proteome in the snail Biomphalaria glabrata.

Freshwater snails are obligate intermediate hosts for numerous parasitic trematodes, most notably schistosomes. Schistosomiasis is a devastating human and veterinary illness, which is primarily controlled by limiting the transmission of these parasites from their intermediate snail hosts. Understanding how this transmission occurs, as well as the basic immunobiology of these snails may be important for controlling this disease in the future. Allelic variation in the Guadeloupe resistance complex (GRC) of Biomphalaria glabrata partially determines their susceptibility to parasitic infection, and can influence the microbiome diversity and microbial defenses in the hemolymph of these snails. In the present study, we examine the most abundant proteins present in the hemolymph of snails that are resistant or susceptible to schistosomes, as determined by their GRC genotype. Using proteomic analysis, we found that snails with different GRC genotypes have differentially abundant hemolymph proteins that are not explained by differences in transcription. There are 13 revealed hemolymph proteins that differ significantly between resistant and susceptible genotypes, nearly 40% of which are involved in immune responses. These findings build on the mounting evidence that genes in the GRC region have multiple physiological roles, and likely contribute more extensively to the general immune response than previously believed. These data also raise the intriguing possibility that the GRC region controls resistance to schistosomes, not directly, but indirectly via its effects on the snail's proteome and potentially its microbiome.

Life history variation is maintained by fitness trade-offs and negative frequency-dependent selection.

The maintenance of diverse life history strategies within and among species remains a fundamental question in ecology and evolutionary biology. By using a near-complete 16-year pedigree of 12,579 winter-run steelhead (Oncorhynchus mykiss) from the Hood River, Oregon, we examined the continued maintenance of two life history traits: the number of lifetime spawning events (semelparous vs. iteroparous) and age at first spawning (2-5 years). We found that repeat-spawning fish had more than 2.5 times the lifetime reproductive success of single-spawning fish. However, first-time repeat-spawning fish had significantly lower reproductive success than single-spawning fish of the same age, suggesting that repeat-spawning fish forego early reproduction to devote additional energy to continued survival. For single-spawning fish, we also found evidence for a fitness trade-off for age at spawning: older, larger males had higher reproductive success than younger, smaller males. For females, in contrast, we found that 3-year-old fish had the highest mean lifetime reproductive success despite the observation that 4- and 5-year-old fish were both longer and heavier. This phenomenon was explained by negative frequency-dependent selection: as 4- and 5-year-old fish decreased in frequency on the spawning grounds, their lifetime reproductive success became greater than that of the 3-year-old fish. Using a combination of mathematical and individual-based models parameterized with our empirical estimates, we demonstrate that both fitness trade-offs and negative frequency-dependent selection observed in the empirical data can theoretically maintain the diverse life history strategies found in this population.

Allelic Variation in a Single Genomic Region Alters the Microbiome of the Snail Biomphalaria glabrata.

Freshwater snails are the intermediate hosts for numerous parasitic worms which can have negative consequences for human health and agriculture. Understanding the transmission of these diseases requires a more complete characterization of the immunobiology of snail hosts. This includes the characterization of its microbiome and genetic factors which may interact with this important commensal community. Allelic variation in the Guadeloupe resistance complex (GRC) genomic region of Guadeloupean Biomphalaria glabrata influences their susceptibility to schistosome infection and may have other roles in the snail immune response. In the present study, we examined whether a snail's GRC genotype has a role in shaping the bacterial diversity and composition present on or in whole snails. We show that the GRC haplotype, including the resistant genotype, has a significant effect on the diversity of bacterial species present in or on whole snails, including the relative abundances of Gemmatimonas aurantiaca and Micavibrio aeruginosavorus. These findings support the hypothesis that the GRC region is likely involved in pathways that can modify the microbial community of these snails and may have more immune roles in B. glabrata than originally believed. This is also one of few examples in which allelic variation at a particular locus has been shown to affect the microbiome in any species.

Clearance of schistosome parasites by resistant genotypes at a single genomic region in Biomphalaria glabrata snails involves cellular components of the hemolymph.

Schistosomiasis is one of the most detrimental neglected tropical diseases. Controlling the spread of this parasitic illness requires effective sanitation, access to chemotherapeutic drugs, and control over populations of the freshwater snails, such as Biomphalaria glabrata, that are essential intermediate hosts for schistosomes. Effectively controlling this disease, while minimising ecological implications of such control, will require an extensive understanding of the immunological interactions between schistosomes and their molluscan intermediate hosts. Here we histologically characterise the clearance of schistosome larvae by snails that exhibit allelic variation at a single genomic region, the Guadeloupe resistance complex. We show that snails with a resistant Guadeloupe resistance complex genotype clear schistosomes within the first 24-48 h, and that this resistance can be transferred to susceptible snails via whole hemolymph but not cell-free plasma. These findings imply that Guadeloupe resistance complex-coded proteins help to coordinate hemocyte-mediated immune responses to schistosome infections in Guadeloupean snails.

Allelic variation partially regulates galactose-dependent hydrogen peroxide release from circulating hemocytes of the snail Biomphalaria glabrata.

Freshwater snails are the intermediate hosts for numerous parasitic worms that are detrimental to human and agricultural health. Understanding the immune responses of these snails could be vital for finding ways to block transmission of those parasites. Allelic variation in a recently discovered genomic region in the snail, Biomphalaria glabrata, influences their susceptibility to schistosomes. Here we tested whether genes in that region, termed the Guadeloupe Resistance Complex (GRC), are involved in recognition of common pathogen-associated molecules that have been shown to be stimulants of the hydrogen peroxide defense pathway. We show that hemocytes extracted from individuals with one of the three GRC genotypes released less hydrogen peroxide than the other two genotypes, after stimulation with galactose. This difference was not observed after stimulation with several other microbial-associated carbohydrates, despite those ligands sharing the same putative pathway for hydrogen peroxide release. Therefore, we conclude that allelic variation in the GRC region may influence the recognition of galactose, rather than the conserved downstream steps in the hydrogen peroxide pathway. These results thus are consistent with the hypothesis that proteins produced by this region are involved in pathogen recognition.

Corrigendum: Whole genome analysis of a schistosomiasis-transmitting freshwater snail.

This corrects the article DOI: 10.1038/ncomms15451.

Whole genome analysis of a schistosomiasis-transmitting freshwater snail.

Biomphalaria snails are instrumental in transmission of the human blood fluke Schistosoma mansoni. With the World Health Organization's goal to eliminate schistosomiasis as a global health problem by 2025, there is now renewed emphasis on snail control. Here, we characterize the genome of Biomphalaria glabrata, a lophotrochozoan protostome, and provide timely and important information on snail biology. We describe aspects of phero-perception, stress responses, immune function and regulation of gene expression that support the persistence of B. glabrata in the field and may define this species as a suitable snail host for S. mansoni. We identify several potential targets for developing novel control measures aimed at reducing snail-mediated transmission of schistosomiasis.

A Targeted Capture Linkage Map Anchors the Genome of the Schistosomiasis Vector Snail, Biomphalaria glabrata.

The aquatic planorbid snail Biomphalaria glabrata is one of the most intensively-studied mollusks due to its role in the transmission of schistosomiasis. Its 916 Mb genome has recently been sequenced and annotated, but it remains poorly assembled. Here, we used targeted capture markers to map over 10,000 B. glabrata scaffolds in a linkage cross of 94 F1 offspring, generating 24 linkage groups (LGs). We added additional scaffolds to these LGs based on linkage disequilibrium (LD) analysis of targeted capture and whole-genome sequences of 96 unrelated snails. Our final linkage map consists of 18,613 scaffolds comprising 515 Mb, representing 56% of the genome and 75% of genic and nonrepetitive regions. There are 18 large (> 10 Mb) LGs, likely representing the expected 18 haploid chromosomes, and > 50% of the genome has been assigned to LGs of at least 17 Mb. Comparisons with other gastropod genomes reveal patterns of synteny and chromosomal rearrangements. Linkage relationships of key immune-relevant genes may help clarify snail-schistosome interactions. By focusing on linkage among genic and nonrepetitive regions, we have generated a useful resource for associating snail phenotypes with causal genes, even in the absence of a complete genome assembly. A similar approach could potentially improve numerous poorly-assembled genomes in other taxa. This map will facilitate future work on this host of a serious human parasite.

Development of new microsatellites for the hookworm Ancylostoma caninum and analysis of genetic diversity in Brazilian populations.

Considering the great efforts towards formulating a vaccine against hookworms, and the concerns about the spread of drug resistance through hookworm populations, it is justified to study the molecular diversity and population genetic structure of these nematodes. This work had the aim to develop microsatellite markers to investigate the genetic structure and the molecular diversity of Brazilian populations of Ancylostoma caninum. Seven microsatellites markers were successfully used to characterize five Brazilian populations. These findings may contribute to a better comprehension of the ecology, patterns of transmission, drug resistances and development of immunotherapeutic strategies in hookworms.

Schistosome infectivity in the snail, Biomphalaria glabrata, is partially dependent on the expression of Grctm6, a Guadeloupe Resistance Complex protein.

Schistosomiasis is one of the most important neglected tropical diseases. Despite effective chemotherapeutic treatments, this disease continues to afflict hundreds of millions of people. Understanding the natural intermediate snail hosts of schistosome parasites is vital to the suppression of this disease. A recently identified genomic region in Caribbean Biomphalaria glabrata snails strongly influences their resistance to infection by Schistosoma mansoni. This region contains novel genes having structural similarity to known pathogen recognition proteins. Here we elaborate on the probable structure and role of one of these genes, grctm6. We characterised the expression of Grctm6 in a population of Caribbean snails, and performed a siRNA knockdown of Grctm6. We show that this protein is not only expressed in B. glabrata hemolymph, but that it also has a role in modulating the number of S. mansoni cercariae released by infected snails, making it a possible target for the biological control of schistosomiasis.

The behavioral effects of antibiotic treatment on the snail Biomphalaria glabrata.

Schistosomiasis is a detrimental neglected tropical disease that is transmitted by Planorbid snails. Understanding the transmission and control of this disease requires an extensive understanding of these intermediate hosts, which is only achieved by the effective rearing and study of species such as Biomphalaria glabrata. This species is the intermediate host for Schistosoma mansoni in the New World, and is also the main model for studying schistosomes in mollusks. Antibiotics are used routinely in B. glabrata tissue culture, and occasionally on live snails. Here we show that standard doses of three common antibiotics (penicillin, streptomycin and gentamicin) drastically diminish the activity of healthy B. glabrata, but that treated snails recover rapidly when placed in fresh water. Ampicillin treated snails did not show altered activity. We suggest that researchers keep these apparent toxicities in mind if a need for antibiotic treatment of live Planorbid snails arises.

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

The genetic underpinnings associated with the earliest stages of plant and animal domestication have remained elusive. Because a genome-wide response to selection can take many generations, the earliest detectable changes associated with domestication may first manifest as heritable changes to global patterns of gene expression. Here, to test this hypothesis, we measured differential gene expression in the offspring of wild and first-generation hatchery steelhead trout (Oncorhynchus mykiss) reared in a common environment. Remarkably, we find that there were 723 genes differentially expressed between the two groups of offspring. Reciprocal crosses reveal that the differentially expressed genes could not be explained by maternal effects or by chance differences in the background levels of gene expression among unrelated families. Gene-enrichment analyses reveal that adaptation to the novel hatchery environment involved responses in wound healing, immunity and metabolism. These findings suggest that the earliest stages of domestication may involve adaptation to highly crowded conditions.