Variation in embryonic mortality and maternal transcript expression among Atlantic cod (Gadus morhua) broodstock: a functional genomics study.

Early life stage mortality is an important issue for Atlantic cod aquaculture, yet the impact of the cod maternal (egg) transcriptome on egg quality and mortality during embryonic development is poorly understood. In the present work, we studied embryonic mortality and maternal transcript expression using eggs from 15 females. Total mortality at 7days post-fertilization (7 dpf, segmentation stage) was used as an indice of egg quality. A 20,000 probe (20K) microarray experiment compared the 7hours post-fertilization (7 hpf, ~2-cell stage) egg transcriptome of the two lowest quality females (>90% mortality at 7 dpf) to that of the highest quality female (~16% mortality at 7 dpf). Forty-three microarray probes were consistently differentially expressed in both low versus high quality egg comparisons (25 higher expressed in low quality eggs, and 18 higher expressed in high quality eggs). The microarray experiment also identified many immune-relevant genes [e.g. interferon (IFN) pathway genes ifngr1 and ifrd1)] that were highly expressed in eggs of all 3 females regardless of quality. Twelve of the 43 candidate egg quality-associated genes, and ifngr1, ifrd1 and irf7, were included in a qPCR study with 7 hpf eggs from all 15 females. Then, the genes that were confirmed by qPCR to be greater than 2-fold differentially expressed between 7 hpf eggs from the lowest and highest quality females (dcbld1, ddc, and acy3 more highly expressed in the 2 lowest quality females; kpna7 and hacd1 more highly expressed in the highest quality female), and the 3 IFN pathway genes, were included in a second qPCR study with unfertilized eggs. While some maternal transcripts included in these qPCR studies were associated with extremes in egg quality, there was little correlation between egg quality and gene expression when all females were considered. Both dcbld1 and ddc showed greater than 100-fold differences in transcript expression between females and were potentially influenced by family. The Atlantic cod ddc (dopa decarboxylase) complete cDNA was characterized, and has a 1461bp open reading frame encoding a 486 amino acid protein that contains all eight residues of the conserved pyridoxal 5'-phosphate binding site including the catalytic lysine. This study provides valuable new information and resources related to the Atlantic cod egg transcriptome. Some of these microarray-identified, qPCR-confirmed, Atlantic cod egg transcripts (e.g. ddc, kpna7) play important roles during embryonic development of other vertebrate species, and may have similar functions in Atlantic cod.

Microevolution in time and space: SNP analysis of historical DNA reveals dynamic signatures of selection in Atlantic cod.

Little is known about how quickly natural populations adapt to changes in their environment and how temporal and spatial variation in selection pressures interact to shape patterns of genetic diversity. We here address these issues with a series of genome scans in four overfished populations of Atlantic cod (Gadus morhua) studied over an 80-year period. Screening of >1000 gene-associated single-nucleotide polymorphisms (SNPs) identified 77 loci that showed highly elevated levels of differentiation, likely as an effect of directional selection, in either time, space or both. Exploratory analysis suggested that temporal allele frequency shifts at certain loci may correlate with local temperature variation and with life history changes suggested to be fisheries induced. Interestingly, however, largely nonoverlapping sets of loci were temporal outliers in the different populations and outliers from the 1928 to 1960 period showed almost complete stability during later decades. The contrasting microevolutionary trajectories among populations resulted in sequential shifts in spatial outliers, with no locus maintaining elevated spatial differentiation throughout the study period. Simulations of migration coupled with observations of temporally stable spatial structure at neutral loci suggest that population replacement or gene flow alone could not explain all the observed allele frequency variation. Thus, the genetic changes are likely to at least partly be driven by highly dynamic temporally and spatially varying selection. These findings have important implications for our understanding of local adaptation and evolutionary potential in high gene flow organisms and underscore the need to carefully consider all dimensions of biocomplexity for evolutionarily sustainable management.