MHC class I expression dependent on bacterial infection and parental factors in whitefish embryos (Salmonidae)

Ecological conditions can influence not only the expression of a phenotype, but also the heritability of a trait. As such, heritable variation for a trait needs to be studied across environments. We have investigated how pathogen challenge affects the expression of MHC genes in embryos of the lake whitefish Coregonus palaea. In order to experimentally separate paternal (i.e. genetic) from maternal and environmental effects, and determine whether and how stress affects the heritable variation for MHC expression, embryos were produced in full‐factorial in vitro fertilizations, reared singly, and exposed at 208 degree days (late‐eyed stage) to either one of two strains of Pseudomonas fluorescens that differ in their virulence characteristics (one increased mortality, while both delayed hatching time). Gene expression was assessed 48 h postinoculation, and virulence effects of the bacterial infection were monitored until hatching. We found no evidence of MHC class II expression at this stage of development. MHC class I expression was markedly down‐regulated in reaction to both pseudomonads. While MHC expression could not be linked to embryo survival, the less the gene was expressed, the earlier the embryos hatched within each treatment group, possibly due to trade‐offs between immune function and developmental rate or further factors that affect both hatching timing and MHC expression. We found significant additive genetic variance for MHC class I expression in some treatments. That is, changes in pathogen pressures could induce rapid evolution in MHC class I expression. However, we found no additive genetic variance in reaction norms in our study population.     

Genetic variation underlies temperature tolerance of embryos in the sea urchin Heliocidaris erythrogramma armigera

Ocean warming can alter natural selection on marine systems, and in many cases, the long‐term persistence of affected populations will depend on genetic adaptation. In this study, we assess the potential for adaptation in the sea urchin Heliocidaris erythrogramma armigera, an Australian endemic, that is experiencing unprecedented increases in ocean temperatures. We used a factorial breeding design to assess the level of heritable variation in larval hatching success at two temperatures. Fertilized eggs from each full‐sibling family were tested at 22 °C (current spawning temperature) and 25 °C (upper limit of predicted warming this century). Hatching success was significantly lower at higher temperatures, confirming that ocean warming is likely to exert selection on this life‐history stage. Our analyses revealed significant additive genetic variance and genotype‐by‐environment interactions underlying hatching success. Consistent with prior work, we detected significant nonadditive (sire‐by‐dam) variance in hatching success, but additionally found that these interactions were modified by temperature. Although these findings suggest the potential for genetic adaptation, any evolutionary responses are likely to be influenced (and possibly constrained) by complex genotype‐by‐environment and sire‐by‐dam interactions and will additionally depend on patterns of genetic covariation with other fitness traits.     

Pedigree‐free quantitative genetic approach provides evidence for heritability of movement tactics in wild roe deer

Assessing the evolutionary potential of animal populations in the wild is crucial to understanding how they may respond to selection mediated by rapid environmental change (e.g. habitat loss and fragmentation). A growing number of studies have investigated the adaptive role of behaviour, but assessments of its genetic basis in a natural setting remain scarce. We combined intensive biologging technology with genome‐wide data and a pedigree‐free quantitative genetic approach to quantify repeatability, heritability and evolvability for a suite of behaviours related to the risk avoidance‐resource acquisition trade‐off in a wild roe deer (Capreolus capreolus) population inhabiting a heterogeneous, human‐dominated landscape. These traits, linked to the stress response, movement and space‐use behaviour, were all moderately to highly repeatable. Furthermore, the repeatable among‐individual component of variation in these traits was partly due to additive genetic variance, with heritability estimates ranging from 0.21 ± 0.08 to 0.70 ± 0.11 and evolvability ranging from 1.1% to 4.3%. Changes in the trait mean can therefore occur under hypothetical directional selection over just a few generations. To the best of our knowledge, this is the first empirical demonstration of additive genetic variation in space‐use behaviour in a free‐ranging population based on genomic relatedness data. We conclude that wild animal populations may have the potential to adjust their spatial behaviour to human‐driven environmental modifications through microevolutionary change.     

Hybridization,polyploidy and clonality influence geographic patterns of diversity and salt tolerance in the model halophyte seashore paspalum (Paspalum vaginatum)

In plant species, variation in levels of clonality, ploidy and interspecific hybridization can interact to influence geographic patterns of genetic diversity. These factors commonly vary in plants that specialize on saline habitats (halophytes) and may play a role in how they adapt to salinity variation across their range. One such halophyte is the turfgrass and emerging genomic model system seashore paspalum (Paspalum vaginatum Swartz). To investigate how clonal propagation, ploidy variation, and interspecific hybridization vary across ecotypes and local salinity levels in wild P. vaginatum, we employed genotyping‐by‐sequencing, cpDNA sequencing and flow cytometry in 218 accessions representing > 170 wild collections from throughout the coastal southern United States plus USDA germplasm. We found that the two morphologically distinct ecotypes of P. vaginatum differ in their adaptive strategies. The fine‐textured ecotype is diploid and appears to reproduce in the wild both sexually and by clonal propagation; in contrast, the coarse‐textured ecotype consists largely of clonally‐propagating triploid and diploid genotypes. The coarse‐textured ecotype appears to be derived from hybridization between fine‐textured P. vaginatum and an unidentified Paspalum species. These clonally propagating hybrid genotypes are more broadly distributed than clonal fine‐textured genotypes and may represent a transition to a more generalist adaptive strategy. Additionally, the triploid genotypes vary in whether they carry one or two copies of the P. vaginatum subgenome, indicating multiple evolutionary origins. This variation in subgenome composition shows associations with local ocean salinity levels across the sampled populations and may play a role in local adaptation.     

The genetic basis of female multiple mating in a polyandrous livebearing fish

The widespread occurrence of female multiple mating (FMM) demands evolutionary explanation, particularly in the light of the costs of mating. One explanation encapsulated by “good sperm” and “sexy‐sperm” (GS‐SS) theoretical models is that FMM facilitates sperm competition, thus ensuring paternity by males that pass on genes for elevated sperm competitiveness to their male offspring. While support for this component of GS‐SS theory is accumulating, a second but poorly tested assumption of these models is that there should be corresponding heritable genetic variation in FMM – the proposed mechanism of postcopulatory preferences underlying GS‐SS models. Here, we conduct quantitative genetic analyses on paternal half‐siblings to test this component of GS‐SS theory in the guppy (Poecilia reticulata), a freshwater fish with some of the highest known rates of FMM in vertebrates. As with most previous quantitative genetic analyses of FMM in other species, our results reveal high levels of phenotypic variation in this trait and a correspondingly low narrow‐sense heritability (h² = 0.11). Furthermore, although our analysis of additive genetic variance in FMM was not statistically significant (probably owing to limited statistical power), the ensuing estimate of mean‐standardized additive genetic variance (I_A = 0.7) was nevertheless relatively low compared with estimates published for life‐history traits across a broad range of taxa. Our results therefore add to a growing body of evidence that FMM is characterized by relatively low additive genetic variation, thus apparently contradicting GS‐SS theory. However, we qualify this conclusion by drawing attention to potential deficiencies in most designs (including ours) that have tested for genetic variation in FMM, particularly those that fail to account for intersexual interactions that underlie FMM in many systems.     

Thermal adaptation and phenotypic plasticity in a warming world: Insights from common garden experiments on Alaskan sockeye salmon

An important unresolved question is how populations of coldwater‐dependent fishes will respond to rapidly warming water temperatures. For example, the culturally and economically important group, Pacific salmon (Oncorhynchus spp.), experience site‐specific thermal regimes during early development that could be disrupted by warming. To test for thermal local adaptation and heritable phenotypic plasticity in Pacific salmon embryos, we measured the developmental rate, survival, and body size at hatching in two populations of sockeye salmon (Oncorhynchus nerka) that overlap in timing of spawning but incubate in contrasting natural thermal regimes. Using a split half‐sibling design, we exposed embryos of 10 families from each of two populations to variable and constant thermal regimes. These represented both experienced temperatures by each population, and predicted temperatures under plausible future conditions based on a warming scenario from the downscaled global climate model (MIROC A1B scenario). We did not find evidence of thermal local adaptation during the embryonic stage for developmental rate or survival. Within treatments, populations hatched within 1 day of each other, on average, and among treatments, did not differ in survival in response to temperature. We did detect plasticity to temperature; embryos developed 2.5 times longer (189 days) in the coolest regime compared to the warmest regime (74 days). We also detected variation in developmental rates among families within and among temperature regimes, indicating heritable plasticity. Families exhibited a strong positive relationship between thermal variability and phenotypic variability in developmental rate but body length and mass at hatching were largely insensitive to temperature. Overall, our results indicated a lack of thermal local adaptation, but a presence of plasticity in populations experiencing contrasting conditions, as well as family‐specific heritable plasticity that could facilitate adaptive change.     

Genotype and diet affect resistance,survival, and fecundity but not fecundity tolerance

Insects are exposed to a variety of potential pathogens in their environment, many of which can severely impact fitness and health. Consequently, hosts have evolved resistance and tolerance strategies to suppress or cope with infections. Hosts utilizing resistance improve fitness by clearing or reducing pathogen loads, and hosts utilizing tolerance reduce harmful fitness effects per pathogen load. To understand variation in, and selective pressures on, resistance and tolerance, we asked to what degree they are shaped by host genetic background, whether plasticity in these responses depends upon dietary environment, and whether there are interactions between these two factors. Females from ten wild‐type Drosophila melanogaster genotypes were kept on high‐ or low‐protein (yeast) diets and infected with one of two opportunistic bacterial pathogens, Lactococcus lactis or Pseudomonas entomophila. We measured host resistance as the inverse of bacterial load in the early infection phase. The relationship (slope) between fly fecundity and individual‐level bacteria load provided our fecundity tolerance measure. Genotype and dietary yeast determined host fecundity and strongly affected survival after infection with pathogenic P. entomophila. There was considerable genetic variation in host resistance, a commonly found phenomenon resulting from for example varying resistance costs or frequency‐dependent selection. Despite this variation and the reproductive cost of higher P. entomophila loads, fecundity tolerance did not vary across genotypes. The absence of genetic variation in tolerance may suggest that at this early infection stage, fecundity tolerance is fixed or that any evolved tolerance mechanisms are not expressed under these infection conditions.     

Evolutionary potential of morphological traits across different life‐history stages

Despite accumulating examples of selection acting on heritable traits in the wild, predicted evolutionary responses are often different from observed phenotypic trends. Various explanations have been suggested for these mismatches. These include within‐individual changes across lifespan that can create important variation in genetic architecture of traits and selection acting on them, but also potential problems with the methodological approach used to predict evolutionary responses of traits. Here, we used an 8‐year data set on tree swallow (Tachycineta bicolor) to first assess the effects of differences among three nestling life‐history stages on the genetic (co)variances of two morphological traits (body mass and primary feather length) and the selection acting on them over three generations. We then estimated the evolutionary potential of these traits by predicting their evolutionary responses using the breeder's equation and the secondary theorem of selection approaches. Our results showed variation in strength and direction of selection and slight changes in trait variance across ages. Predicted evolutionary responses differed importantly between both approaches for half of the trait–age combinations we studied, suggesting the presence of environmentally induced correlations between focal traits and fitness possibly biasing breeder's equation predictions. Our results emphasize that predictions of evolutionary potential for morphological traits are likely to be highly variable, both in strength and direction, depending on the life stage and method used, thus mitigating our capacity to predict adaptation and persistence of wild populations.     

Genetic variation in variability: Phenotypic variability of fledging weight and its evolution in a songbird population

Variation in traits is essential for natural selection to operate and genetic and environmental effects can contribute to this phenotypic variation. From domesticated populations, we know that families can differ in their level of within‐family variance, which leads to the intriguing situation that within‐family variance can be heritable. For offspring traits, such as birth weight, this implies that within‐family variance in traits can vary among families and can thus be shaped by natural selection. Empirical evidence for this in wild populations is however lacking. We investigated whether within‐family variance in fledging weight is heritable in a wild great tit (Parus major) population and whether these differences are associated with fitness. We found significant evidence for genetic variance in within‐family variance. The genetic coefficient of variation (GCV) was 0.18 and 0.25, when considering fledging weight a parental or offspring trait, respectively. We found a significant quadratic relationship between within‐family variance and fitness: families with low or high within‐family variance had lower fitness than families with intermediate within‐family variance. Our results show that within‐family variance can respond to selection and provides evidence for stabilizing selection on within‐family variance.     

Opportunities and challenges of next‐generation sequencing applications in ecological epigenetics

Evolutionary theory posits that adaptation can result when populations harbour heritable phenotypic variation for traits that increase tolerance to local conditions. However, the actual mechanisms that underlie heritable phenotypic variation are not completely understood (Keller 2014 ). Recently, the potential role of epigenetic mechanisms in the process of adaptive evolution has been the subject of much debate (Pigliucci & Finkelman 2014 ). Studies of variation in DNA methylation in particular have shown that natural populations harbour high amounts of epigenetic variation, which can be inherited across generations and can cause heritable trait variation independently of genetic variation (Kilvitis et al. 2014 ). While we have made some progress addressing the importance of epigenetics in ecology and evolution using methylation‐sensitive AFLP (MS‐AFLP), this approach provides relatively few anonymous and dominant markers per individual. MS‐AFLP are difficult to link to functional genomic elements or phenotype and are difficult to compare directly to genetic variation, which has limited the insights drawn from studies of epigenetic variation in natural nonmodel populations (Schrey et al. 2013 ). In this issue, Platt et al. provide an example of a promising approach to address this problem by applying a reduced representation bisulphite sequencing (RRBS) approach based on next‐generation sequencing methods in an ecological context.     

Genetic basis and fitness correlates of dynamic carotenoid‐based ornamental coloration in male and female common kestrels Falco tinnunculus

Knowledge of the genetic basis of sexual ornaments is essential to understand their evolution through sexual selection. Although carotenoid‐based ornaments have been instrumental in the study of sexual selection, given the inability of animals to synthesize carotenoids de novo, they are generally assumed to be influenced solely by environmental variation. However, very few studies have directly estimated the role of genes and the environment in shaping variation in carotenoid‐based traits. Using long‐term individual‐based data, we here explore the evolutionary potential of a dynamic, carotenoid‐based ornament (namely skin coloration), in male and female common kestrels. We first estimate the amount of genetic variation underlying variation in hue, chroma and brightness. After correcting for sex differences, the chroma of the orange‐yellow eye ring coloration was significantly heritable (h² ± SE = 0.40 ± 0.17), whereas neither hue (h² = 0) nor brightness (h² = 0.02) was heritable. Second, we estimate the strength and shape of selection acting upon chromatic (hue and chroma) and achromatic (brightness) variation and show positive and negative directional selection on female but not male chroma and hue, respectively, whereas brightness was unrelated to fitness in both sexes. This suggests that different components of carotenoid‐based signals traits may show different evolutionary dynamics. Overall, we show that carotenoid‐based coloration is a complex and multifaceted trait. If we are to gain a better understanding of the processes responsible for the generation and maintenance of variation in carotenoid‐based coloration, these complexities need to be taken into account.     

Risk‐induced hatching timing shows low heritability and evolves independently of spontaneous hatching in red‐eyed treefrogs

Plasticity in the timing of transitions between stages of complex life cycles allows organisms to adjust their growth and development to local environmental conditions. Genetic variation in such plasticity is common, but the evolution of context‐dependent transition timing may be constrained by information reliability, lag‐time and developmental constraints. We studied the genetic architecture of hatching plasticity in embryos of the red‐eyed treefrog (Agalychnis callidryas) in response to simulated predator attacks using a series of paternal and maternal half‐sibs from a captive breeding colony of wild‐collected animals. We compared the developmental timing of induced early hatching across sibships and estimated cross‐environment genetic correlations between induced and spontaneous hatching traits. Additive genetic variance for induced early hatching was very low, indicating a constraint on the short‐term evolution of earlier hatching timing. This constraint is likely related to the maturation of the hatching mechanism. The most plastic genotypes produced the most extreme spontaneous hatching phenotypes, indicating that developmental range, per se, is not constrained. Cross‐environment genetic correlation in hatching timing was negligible, so the evolution of spontaneous hatching in this species has not depended on the evolution of risk‐induced hatching and vice versa.     

Heritability and genetic correlations of personality traits in a wild population of yellow‐bellied marmots (Marmota flaviventris)

Describing and quantifying animal personality is now an integral part of behavioural studies because individually distinctive behaviours have ecological and evolutionary consequences. Yet, to fully understand how personality traits may respond to selection, one must understand the underlying heritability and genetic correlations between traits. Previous studies have reported a moderate degree of heritability of personality traits, but few of these studies have either been conducted in the wild or estimated the genetic correlations between personality traits. Estimating the additive genetic variance and covariance in the wild is crucial to understand the evolutionary potential of behavioural traits. Enhanced environmental variation could reduce heritability and genetic correlations, thus leading to different evolutionary predictions. We estimated the additive genetic variance and covariance of docility in the trap, sociability (mirror image stimulation), and exploration and activity in two different contexts (open‐field and mirror image simulation experiments) in a wild population of yellow‐bellied marmots (Marmota flaviventris). We estimated both heritability of behaviours and of personality traits and found nonzero additive genetic variance in these traits. We also found nonzero maternal, permanent environment and year effects. Finally, we found four phenotypic correlations between traits, and one positive genetic correlation between activity in the open‐field test and sociability. We also found permanent environment correlations between activity in both tests and docility and exploration in the MIS test. This is one of a handful of studies to adopt a quantitative genetic approach to explain variation in personality traits in the wild and, thus, provides important insights into the potential variance available for selection.     

Potential for adaptation in response to thermal stress in an intertidal macroalga

Understanding responses of marine algae to changing ocean temperatures requires knowledge of the impacts of elevated temperatures and the likelihood of adaptation to thermal stress. The potential for rapid evolution of thermal tolerance is dependent on the levels of heritable genetic variation in response to thermal stress within a population. Here, we use a quantitative genetic breeding design to establish whether there is a heritable variation in thermal sensitivity in two populations of a habitat‐forming intertidal macroalga, Hormosira banksii (Turner) Descaisne. Gametes from multiple parents were mixed and growth and photosynthetic performance were measured in the resulting embryos, which were incubated under control and elevated temperature (20°C and 28°C). Embryo growth was reduced at 28°C, but significant interactions between male genotype and temperature in one population indicated the presence of genetic variation in thermal sensitivity. Selection for more tolerant genotypes thus has the ability to result in the evolution of increased thermal tolerance. Furthermore, genetic correlations between embryos grown in the two temperatures were positive, indicating that those genotypes that performed well in elevated temperature also performed well in control temperature. Chlorophyll a fluorescence measurements showed a marked decrease in maximum quantum yield of photosystem II (PSII) under elevated temperature. There was an increase in the proportion of energy directed to photoinhibition (nonregulated nonphotochemical quenching) and a concomitant decrease in energy used to drive photochemistry and xanthophyll cycling (regulated nonphotochemical quenching). However, PSII performance between genotypes was similar, suggesting that thermal sensitivity is related to processes other than photosynthesis.     

Heritability and cross-sex genetic correlations of early-life circulating testosterone levels in a wild mammal

Testosterone is an important hormone that has been shown to have sex-specific links to fitness in numerous species. Although testosterone concentrations vary substantially between individuals in a population, little is known about its heritable genetic basis or between-sex genetic correlations that determine its evolutionary potential. We found circulating neonatal testosterone levels to be both heritable (0.160 ± 0.064 s.e.) and correlated between the sexes (0.942 ± 0.648 s.e.) in wild red deer calves (Cervus elaphus). This may have important evolutionary implications if, as in adults, the sexes have divergent optima for circulating testosterone levels.     

Uncovering symbiont‐driven genetic diversity across North American pea aphids

Heritable genetic variation is required for evolution, and while typically encoded within nuclear and organellar genomes, several groups of invertebrates harbour heritable microbes serving as additional sources of genetic variation. Hailing from the symbiont‐rich insect order Hemiptera, pea aphids (Acyrthosiphon pisum) possess several heritable symbionts with roles in host plant utilization, thermotolerance and protection against natural enemies. As pea aphids vary in the numbers and types of harboured symbionts, these bacteria provide heritable and functionally important variation within field populations. In this study, we quantified the cytoplasmically inherited genetic variation contributed by symbionts within North American pea aphids. Through the use of Denaturing Gradient Gel Electrophoresis (DGGE) and 454 amplicon pyrosequencing of 16S rRNA genes, we explored the diversity of bacteria harboured by pea aphids from five populations, spanning three locations and three host plants. We also characterized strain variation by analysing 16S rRNA, housekeeping and symbiont‐associated bacteriophage genes. Our results identified eight species of facultative symbionts, which often varied in frequency between locations and host plants. We detected 28 cytoplasmic genotypes across 318 surveyed aphids, considering only the various combinations of secondary symbiont species infecting single hosts. Yet the detection of multiple Regiella insecticola, Hamiltonella defensa and Rickettsia strains, and diverse bacteriophage genotypes from H. defensa, suggest even greater diversity. Combined, these findings reveal that heritable bacteria contribute substantially to genetic variation in A. pisum. Given the costs and benefits of these symbionts, it is likely that fluctuating selective forces play a role in the maintenance of this diversity.     

Small‐scale spatial variation in evolvability for life‐history traits in the storm petrel

The evolutionary potential in the timing of recruitment and reproduction may be crucial for the ability of populations to buffer against environmental changes, allowing them to avoid unfavourable breeding conditions. The evolution of a trait in a local population is determined by its heritability and selection. In the present study, we performed pedigree‐based quantitative genetic analyses for two life‐history traits (recruiting age and laying date) using population data of the storm petrel over an 18‐year period in two adjacent breeding colonies (only 150 m apart) that share the same environmental conditions. In both traits, natal colony effect was the main source of the phenotypic variation among individuals, and cohort variance for recruitment age and additive genetic variance for laying date were natal colony‐specific. We found significant heritability only in laying date and, more specifically, only in birds born in one of the colonies. The difference in genetic variance between the colonies was statistically significant. Interestingly, selection on earlier breeding birds was detected only in the colony in which heritable variation in laying date was found. Therefore, local evolvability for a life‐history trait may vary within a unexpectedly small spatial scale, through the diversifying natural selection and insulating gene flow. © 2012 The Linnean Society of London, Biological Journal of the Linnean Society, 2012, 106 , 439–446.     

A de novo chromosome‐level genome assembly of Coregonus sp. “Balchen”: One representative of the Swiss Alpine whitefish radiation

Salmonids are of particular interest to evolutionary biologists due to their incredible diversity of life‐history strategies and the speed at which many salmonid species have diversified. In Switzerland alone, over 30 species of Alpine whitefish from the subfamily Coregoninae have evolved since the last glacial maximum, with species exhibiting a diverse range of morphological and behavioural phenotypes. This, combined with the whole genome duplication which occurred in the ancestor of all salmonids, makes the Alpine whitefish radiation a particularly interesting system in which to study the genetic basis of adaptation and speciation and the impacts of ploidy changes and subsequent rediploidization on genome evolution. Although well‐curated genome assemblies exist for many species within Salmonidae, genomic resources for the subfamily Coregoninae are lacking. To assemble a whitefish reference genome, we carried out PacBio sequencing from one wild‐caught Coregonus sp. “Balchen” from Lake Thun to ~90× coverage. PacBio reads were assembled independently using three different assemblers, falcon , canu and wtdbg2 and subsequently scaffolded with additional Hi‐C data. All three assemblies were highly contiguous, had strong synteny to a previously published Coregonus linkage map, and when mapping additional short‐read data to each of the assemblies, coverage was fairly even across most chromosome‐scale scaffolds. Here, we present the first de novo genome assembly for the Salmonid subfamily Coregoninae. The final 2.2‐Gb wtdbg2 assembly included 40 scaffolds, an N50 of 51.9 Mb and was 93.3% complete for BUSCOs. The assembly consisted of ~52% transposable elements and contained 44,525 genes.     

Tolerance of whitefish embryos to Pseudomonas fluorescens linked to genetic and maternal effects,and reduced by previous exposure

Juvenile or adult fish can alter their behaviour and rely on an innate and adaptive immune system to avoid/counteract pathogens, while fish embryos have to depend on egg characteristics and may be partly protected by their developing immune system that is building up from a certain age on. We developed an infection protocol that allows testing the reaction of individual whitefish embryos (Coregonus palaea) to repeated exposures to Pseudomonas fluorescens, an opportunistic bacterial fish pathogen. We used a full-factorial in vitro breeding design to separately test the effects of paternal and maternal contributions to the embryos' susceptibility to different kinds of pathogen exposure. We found that a first non-lethal exposure had immunosuppressive effects: pre-exposed embryos were more susceptible to future challenges with the same pathogen. At intermediate and high levels of pathogen intensity, maternal effects turned out to be crucial for the embryos' tolerance to infection. Paternal (i.e. genetic) effects played a significant role at the strongest level of infection, i.e. the embryos' own genetics already explained some of the variation in embryo susceptibility. Our findings suggest that whitefish embryos are largely protected by maternally transmitted substances, but build up some own innate immunocompetence several days before hatching.