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.

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.

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.

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.

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.

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.

Genotypic variation in host response to infection affects parasite reproductive rate.

Parasite fitness is largely influenced by a variation in host response due to the host's genetic background. Here we investigated the impact of host genotype on pathogen success in the snail vector of its castrating parasite, Schistosoma mansoni. We infected five inbred lines of Biomphalaria glabrata with two infection doses and followed their growth, reproductive output and parasite production throughout the course of infection. There was no difference in resistance to infection among inbred lines, but lines varied in their responses to infection and the numbers of parasites produced. Snails did not compensate for castration by increasing their fecundity during the early phase of infection (fecundity compensation). However, some lines were able to delay parasite shedding for up to 30 weeks, thus prolonging reproduction before the onset of castration. Here we propose this strategy as a novel defense against castrating pathogens in snails. Gigantism, a predicted outcome of castration due to energy reallocation, occurred early in infection (<15 weeks) and was not universal among the snail lines. Lines that did not show gigantism were also characterised by a high parasite production rate and low survivorship, perhaps indicating energy reallocation into parasite production and costly immune defense. We observed no differences in total parasite production among lines throughout the entire course of infection, although lines differed in their parasite reproductive rate. The average rate of parasite production varied among lines from 1300 to 2450 cercariae within a single 2h shedding period, resulting in a total production of 6981-29,509 cercariae over the lifetime of a single snail. Regardless of genetic background, snail size was a strong predictor of parasite reproduction: each millimetre increase in snail size at the time of the first shed resulted in up to 3500 more cercariae over the lifetime of the snail. The results of this study provide a detailed picture of variation in hosts' responses to infection and the resulting impacts on parasite fitness, further defining the intricacies of snail-schistosome compatibility.

Genome-Wide Scan and Test of Candidate Genes in the Snail Biomphalaria glabrata Reveal New Locus Influencing Resistance to Schistosoma mansoni.

BACKGROUND: New strategies to combat the global scourge of schistosomiasis may be revealed by increased understanding of the mechanisms by which the obligate snail host can resist the schistosome parasite. However, few molecular markers linked to resistance have been identified and characterized in snails. METHODOLOGY/PRINCIPAL FINDINGS: Here we test six independent genetic loci for their influence on resistance to Schistosoma mansoni strain PR1 in the 13-16-R1 strain of the snail Biomphalaria glabrata. We first identify a genomic region, RADres, showing the highest differentiation between susceptible and resistant inbred lines among 1611 informative restriction-site associated DNA (RAD) markers, and show that it significantly influences resistance in an independent set of 439 outbred snails. The additive effect of each RADres resistance allele is 2-fold, similar to that of the previously identified resistance gene sod1. The data fit a model in which both loci contribute independently and additively to resistance, such that the odds of infection in homozygotes for the resistance alleles at both loci (13% infected) is 16-fold lower than the odds of infection in snails without any resistance alleles (70% infected). Genome-wide linkage disequilibrium is high, with both sod1 and RADres residing on haplotype blocks >2 Mb, and with other markers in each block also showing significant effects on resistance; thus the causal genes within these blocks remain to be demonstrated. Other candidate loci had no effect on resistance, including the Guadeloupe Resistance Complex and three genes (aif, infPhox, and prx1) with immunological roles and expression patterns tied to resistance, which must therefore be trans-regulated. CONCLUSIONS/SIGNIFICANCE: The loci RADres and sod1 both have strong effects on resistance to S. mansoni. Future approaches to control schistosomiasis may benefit from further efforts to characterize and harness this natural genetic variation.

Hyperdiverse gene cluster in snail host conveys resistance to human schistosome parasites.

Schistosomiasis, a neglected global pandemic, may be curtailed by blocking transmission of the parasite via its intermediate hosts, aquatic snails. Elucidating the genetic basis of snail-schistosome interaction is a key to this strategy. Here we map a natural parasite-resistance polymorphism from a Caribbean population of the snail Biomphalaria glabrata. In independent experimental evolution lines, RAD genotyping shows that the same genomic region responds to selection for resistance to the parasite Schistosoma mansoni. A dominant allele in this region conveys an 8-fold decrease in the odds of infection. Fine-mapping and RNA-Seq characterization reveal a <1Mb region, the Guadeloupe Resistance Complex (GRC), with 15 coding genes. Seven genes are single-pass transmembrane proteins with putative immunological roles, most of which show strikingly high nonsynonymous divergence (5-10%) among alleles. High linkage disequilibrium among three intermediate-frequency (>25%) haplotypes across the GRC, a significantly non-neutral pattern, suggests that balancing selection maintains diversity at the GRC. Thus, the GRC resembles immune gene complexes seen in other taxa and is likely involved in parasite recognition. The GRC is a potential target for controlling transmission of schistosomiasis, including via genetic manipulation of snails.

On the reproductive success of early-generation hatchery fish in the wild.

Large numbers of hatchery salmon spawn in wild populations each year. Hatchery fish with multiple generations of hatchery ancestry often have heritably lower reproductive success than wild fish and may reduce the fitness of an entire population. Whether this reduced fitness also occurs for hatchery fish created with local- and predominantly wild-origin parents remains controversial. Here, we review recent studies on the reproductive success of such 'early-generation' hatchery fish that spawn in the wild. Combining 51 estimates from six studies on four salmon species, we found that (i) early-generation hatchery fish averaged only half the reproductive success of their wild-origin counterparts when spawning in the wild, (ii) the reduction in reproductive success was more severe for males than for females, and (iii) all species showed reduced fitness due to hatchery rearing. We review commonalities among studies that point to possible mechanisms (e.g., environmental versus genetic effects). Furthermore, we illustrate that sample sizes typical of these studies result in low statistical power to detect fitness differences unless the differences are substantial. This review demonstrates that reduced fitness of early-generation hatchery fish may be a general phenomenon. Future research should focus on determining the causes of those fitness reductions and whether they lead to long-term reductions in the fitness of wild populations.

Sequencing and characterization of the anadromous steelhead (Oncorhynchus mykiss) transcriptome.

Identifying the traits that differ between hatchery and wild fish may allow for pragmatic changes to hatchery practice. To meet those ends, we sequenced, assembled, and characterized the anadromous steelhead (Oncorhynchus mykiss) transcriptome. Using the Illumina sequencing platform, we sequenced nearly 41million 76-mer reads representing 3.1 Gbp of steelhead transcriptome. Upon final assembly, this sequence data yielded 86,402 transcript scaffolds, of which, 66,530 (77%) displayed homology to proteins of the non-redundant NCBI database. Gene descriptions and gene ontology terms were used to annotate the transcriptome resulting in 4030 unique gene ontology (GO) annotations attributed to the assembled sequences. We also conducted a comparative analysis that identified homologous genes within four other fish species including zebrafish (Danio rerio), stickleback (Gasterosteus aculeatus), and two pufferfish species (Tetraodon nigroviridis and Takifugu rubripes). Comparing our steelhead reference assembly directly to the transcriptome for rainbow trout (the fresh water life-history variant of the same species) revealed that while the steelhead and rainbow trout transcriptomes are complementary, the steelhead data will be useful for investigating questions related to anadromous (ocean-going) fishes. These sequence data and web tools provide a useful set of resources for salmonid researchers and the broader genomics community (available at http://salmon.cgrb.oregonstate.edu).

How much does inbreeding contribute to the reduced fitness of hatchery-born steelhead (Oncorhynchus mykiss) in the wild?

Many declining populations are supplemented with captive-born individuals that are released directly into the wild. Because captive-born individuals can have lower fitness in the wild than their wild-born counterparts, a comprehensive understanding of the mechanisms responsible for the reduced fitness of these individuals is required for appropriate conservation and management decisions. Inbreeding among captive-born individuals is one plausible mechanism because captive breeding programs frequently use small numbers of breeders to create large numbers of siblings that are subsequently released together into the wild. We tested this hypothesis in a supplementation program for steelhead (Oncorhynchus mykiss) from the Hood River, Oregon, for which first-generation hatchery fish were demonstrated to have lower fitness in the wild than their wild-born counterparts. To determine the contribution of inbreeding to this fitness decline, we first assigned 11 run-years of hatchery steelhead (3005 fish) back to their broodstock parents (462 fish) using 8 polymorphic microsatellite loci. By combining pedigree analyses with species-specific estimates of genetic load, we found that inbreeding could at most account for a 1-4% reduction in the fitness of hatchery fish relative to wild fish. Thus, inbreeding alone cannot adequately explain the 15% average fitness decline observed in first-generation hatchery fish from this population.