Genomic islands of divergence and their consequences for the resolution of spatial structure in an exploited marine fish.

As populations diverge, genomic regions associated with adaptation display elevated differentiation. These genomic islands of adaptive divergence can inform conservation efforts in exploited species, by refining the delineation of management units, and providing genomic tools for more precise and effective population monitoring and the successful assignment of individuals and products. We explored heterogeneity in genomic divergence and its impact on the resolution of spatial population structure in exploited populations of Atlantic cod, Gadus morhua, using genome wide expressed sequence derived single nucleotide polymorphisms in 466 individuals sampled across the range. Outlier tests identified elevated divergence at 5.2% of SNPs, consistent with directional selection in one-third of linkage groups. Genomic regions of elevated divergence ranged in size from a single position to several cM. Structuring at neutral loci was associated with geographic features, whereas outlier SNPs revealed genetic discontinuities in both the eastern and western Atlantic. This fine-scale geographic differentiation enhanced assignment to region of origin, and through the identification of adaptive diversity, fundamentally changes how these populations should be conserved. This work demonstrates the utility of genome scans for adaptive divergence in the delineation of stock structure, the traceability of individuals and products, and ultimately a role for population genomics in fisheries conservation.

Effects of body position, imaging plane, and observer on computed tomographic measurements of the lumbosacral intervertebral foraminal area in dogs.

OBJECTIVE: To evaluate effects of imaging plane, flexion and extension, patient weight, and observer on computed tomographic (CT) image measurements of the area of the lumbosacral (L7-S1) intervertebral foramen (LSIF) in dogs. SAMPLE: 12 dog cadavers (2 were excluded because of foraminal stenosis). PROCEDURES: In each cadaver, sagittal, sagittal oblique, transverse oblique, and double oblique CT images were obtained at 3 zones (entrance, middle, and exit zones) of the region of the lateral lumbar spinal canal that comprises the LSIF while the lumbosacral junction (LSJ) was positioned in flexion or extension. Barium-impregnated polymethylmethacrylate was used to fill the intervertebral foramina to aid boundary detection. Measurements of interest were obtained. RESULTS: Among the dog cadavers, there was large variability in LSIF cross-sectional areas (range, 0.12 to 0.44 cm2; SD, 0.1 cm2) and in foraminal angles required to obtain a double oblique plane in LSJ extension (SD, 8 ° to 9 °). For LSIF area measurements in standard sagittal CT images, interobserver variability was 23% to 44% and intraobserver variability was 4% to 5%. Sagittal oblique images obtained during LSJ extension yielded smaller mean LSIF areas (0.30 cm2), compared with findings in sagittal images (0.37 to 0.52 cm2). The exit and middle zone areas were smaller than the entrance zone area in sagittal images obtained during LSJ extension. CONCLUSIONS AND CLINICAL RELEVANCE: Repeated measurements of the LSIF area in images obtained during LSJ extension may be unreliable as a result of interobserver variability and the effects of dog positioning and CT slice orientation.

Evaluating SNP ascertainment bias and its impact on population assignment in Atlantic cod, Gadus morhua.

The increasing use of single nucleotide polymorphisms (SNPs) in studies of nonmodel organisms accentuates the need to evaluate the influence of ascertainment bias on accurate ecological or evolutionary inference. Using a panel of 1641 expressed sequence tag-derived SNPs developed for northwest Atlantic cod (Gadus morhua), we examined the influence of ascertainment bias and its potential impact on assignment of individuals to populations ranging widely in origin. We hypothesized that reductions in assignment success would be associated with lower diversity in geographical regions outside the location of ascertainment. Individuals were genotyped from 13 locations spanning much of the contemporary range of Atlantic cod. Diversity, measured as average sample heterozygosity and number of polymorphic loci, declined (c. 30%) from the western (H(e) = 0.36) to eastern (H(e) = 0.25) Atlantic, consistent with a signal of ascertainment bias. Assignment success was examined separately for pools of loci representing differing degrees of reductions in diversity. SNPs displaying the largest declines in diversity produced the most accurate assignment in the ascertainment region (c. 83%) and the lowest levels of correct assignment outside the ascertainment region (c. 31%). Interestingly, several isolated locations showed no effect of assignment bias and consistently displayed 100% correct assignment. Contrary to expectations, estimates of accurate assignment range-wide using all loci displayed remarkable similarity despite reductions in diversity. Our results support the use of large SNP panels in assignment studies of high geneflow marine species. However, our evidence of significant reductions in assignment success using some pools of loci suggests that ascertainment bias may influence assignment results and should be evaluated in large-scale assignment studies.

An integrated approach to gene discovery and marker development in Atlantic cod (Gadus morhua).

Atlantic cod is a species that has been overexploited by the capture fishery. Programs to domesticate this species are underway in several countries, including Canada, to provide an alternative route for production. Selective breeding programs have been successfully applied in the domestication of other species, with genomics-based approaches used to augment conventional methods of animal production in recent years. Genomics tools, such as gene sequences and sets of variable markers, also have the potential to enhance and accelerate selective breeding programs in aquaculture, and to provide better monitoring tools to ensure that wild cod populations are well managed. We describe the generation of significant genomics resources for Atlantic cod through an integrated genomics/selective breeding approach. These include 158,877 expressed sequence tags (ESTs), a set of annotated putative transcripts and several thousand single nucleotide polymorphism markers that were developed from, and have been shown to be highly variable in, fish enrolled in two selective breeding programs. Our EST collection was generated from various tissues and life cycle stages. In some cases, tissues from which libraries were generated were isolated from fish exposed to stressors, including elevated temperature, or antigen stimulation (bacterial and viral) to enrich for transcripts that are involved in these response pathways. The genomics resources described here support the developing aquaculture industry, enabling the application of molecular markers within selective breeding programs. Marker sets should also find widespread application in fisheries management.

Development and experimental validation of a 20K Atlantic cod (Gadus morhua) oligonucleotide microarray based on a collection of over 150,000 ESTs.

The collapse of Atlantic cod (Gadus morhua) wild populations strongly impacted the Atlantic cod fishery and led to the development of cod aquaculture. In order to improve aquaculture and broodstock quality, we need to gain knowledge of genes and pathways involved in Atlantic cod responses to pathogens and other stressors. The Atlantic Cod Genomics and Broodstock Development Project has generated over 150,000 expressed sequence tags from 42 cDNA libraries representing various tissues, developmental stages, and stimuli. We used this resource to develop an Atlantic cod oligonucleotide microarray containing 20,000 unique probes. Selection of sequences from the full range of cDNA libraries enables application of the microarray for a broad spectrum of Atlantic cod functional genomics studies. We included sequences that were highly abundant in suppression subtractive hybridization (SSH) libraries, which were enriched for transcripts responsive to pathogens or other stressors. These sequences represent genes that potentially play an important role in stress and/or immune responses, making the microarray particularly useful for studies of Atlantic cod gene expression responses to immune stimuli and other stressors. To demonstrate its value, we used the microarray to analyze the Atlantic cod spleen response to stimulation with formalin-killed, atypical Aeromonas salmonicida, resulting in a gene expression profile that indicates a strong innate immune response. These results were further validated by quantitative PCR analysis and comparison to results from previous analysis of an SSH library. This study shows that the Atlantic cod 20K oligonucleotide microarray is a valuable new tool for Atlantic cod functional genomics research.

Integrating the markers Pan I and haemoglobin with the genetic linkage map of Atlantic cod (Gadus morhua).

BACKGROUND: Haemoglobin (Hb) and pantophysin (Pan I) markers have been used intensively in population studies of Atlantic cod (Gadus morhua) and in the analysis of traits such as temperature tolerance, growth characteristics and sexual maturation. We used an Illumina GoldenGate panel and the KASPar SNP genotyping system to analyse SNPs in three Atlantic cod families, one of which was polymorphic at the Hb ß1 locus, and to generate a genetic linkage map integrating Pan I and multiple Hb loci. FINDINGS: Data generated allowed the mapping of nine Hb loci, the Pan I locus, and other 122 SNPs onto an existing linkage genetic map for Atlantic cod. Four Hb genes (i.e. α1, α4, ß1 and ß5) have been mapped on linkage group (LG) 2 while the other five (i.e. α2, α3, ß2, ß3 and ß4) were placed on LG18. Pan I was mapped on LG 1 using a newly developed KASPar assay for a SNP variable only in Pan IA allelic variants. The new linkage genetic map presented here comprises 1046 SNPs distributed between 23 linkage groups, with a length of 1145.6 cM. A map produced by forcing additional loci, resulting in a reduced goodness-of-fit for mapped markers, allowed the mapping of a total of 1300 SNPs. Finally, we compared our genetic linkage map data with the genetic linkage map data produced by a different group and identified 29 shared SNPs distributed on 10 different linkage groups. CONCLUSIONS: The genetic linkage map presented here incorporates the marker Pan I, together with multiple Hb loci, and integrates genetic linkage data produced by two different research groups. This represents a useful resource to further explore if Pan I and Hbs or other genes underlie quantitative trait loci (QTL) for temperature sensitivity/tolerance or other phenotypes.

Barium toxicosis in a dog.

CASE DESCRIPTION: A 2-year-old 14.9-kg (32.8-lb) neutered female Shetland Sheepdog was admitted to the University of Liverpool Small Animal Teaching Hospital for evaluation of acute collapse. CLINICAL FINDINGS: At admission, the dog was tachypneic and had reduced limb reflexes and muscle tone in all limbs consistent with diffuse lower motor neuron dysfunction. The dog was severely hypokalemic (1.7 mEq/L; reference range, 3.5 to 5.8 mEq/L). Clinical status of the dog deteriorated; there was muscle twitching, flaccid paralysis, and respiratory failure, which was considered a result of respiratory muscle weakness. Ventricular arrhythmias and severe acidemia (pH, 7.18; reference range, 7.35 to 7.45) developed. Intoxication was suspected, and plasma and urine samples submitted for barium analysis had barium concentrations comparable with those reported in humans with barium toxicosis. Analysis of barium concentrations in 5 control dogs supported the diagnosis of barium toxicosis in the dog. TREATMENT AND OUTCOME: Fluids and potassium supplementation were administered IV. The dog recovered rapidly. Electrolyte concentrations measured after recovery were consistently unremarkable. Quantification of plasma barium concentration 56 days after the presumed episode of intoxication revealed a large decrease; however, the plasma barium concentration remained elevated, compared with that in control dogs. CLINICAL RELEVANCE: To our knowledge, this case represented the first description of barium toxicosis in the veterinary literature. Barium toxicosis can cause life-threatening hypokalemia; however, prompt supportive treatment can yield excellent outcomes. Barium toxicosis is a rare but important differential diagnosis in animals with hypokalemia and appropriate clinical signs.

Parallel adaptive evolution of Atlantic cod on both sides of the Atlantic Ocean in response to temperature.

Despite the enormous economic and ecological importance of marine organisms, the spatial scales of adaptation and biocomplexity remain largely unknown. Yet, the preservation of local stocks that possess adaptive diversity is critical to the long-term maintenance of productive stable fisheries and ecosystems. Here, we document genomic evidence of range-wide adaptive differentiation in a broadcast spawning marine fish, Atlantic cod (Gadus morhua), using a genome survey of single nucleotide polymorphisms. Of 1641 gene-associated polymorphisms examined, 70 (4.2%) tested positive for signatures of selection using a Bayesian approach. We identify a subset of these loci (n=40) for which allele frequencies show parallel temperature-associated clines (p<0.001, r2=0.89) in the eastern and western north Atlantic. Temperature associations were robust to the statistical removal of geographic distance or latitude effects, and contrasted 'neutral' loci, which displayed no temperature association. Allele frequencies at temperature-associated loci were significantly correlated, spanned three linkage groups and several were successfully annotated supporting the involvement of multiple independent genes. Our results are consistent with the evolution and/or selective sweep of multiple genes in response to ocean temperature, and support the possibility of a new conservation paradigm for non-model marine organisms based on genomic approaches to resolving functional and adaptive diversity.

Development of a SNP resource and a genetic linkage map for Atlantic cod (Gadus morhua).

BACKGROUND: Atlantic cod (Gadus morhua) is a species with increasing economic significance for the aquaculture industry. The genetic improvement of cod will play a critical role in achieving successful large-scale aquaculture. While many microsatellite markers have been developed in cod, the number of single nucleotide polymorphisms (SNPs) is currently limited. Here we report the identification of SNPs from sequence data generated by a large-scale expressed sequence tag (EST) program, focusing on fish originating from Canadian waters. RESULTS: A total of 97976 ESTs were assembled to generate 13448 contigs. We detected 4753 SNPs that met our selection criteria (depth of coverage > or = 4 reads; minor allele frequency > 25%). 3072 SNPs were selected for testing. The percentage of successful assays was 75%, with 2291 SNPs amplifying correctly. Of these, 607 (26%) SNPs were monomorphic for all populations tested. In total, 64 (4%) of SNPs are likely to represent duplicated genes or highly similar members of gene families, rather than alternative alleles of the same gene, since they showed a high frequency of heterozygosity. The remaining polymorphic SNPs (1620) were categorised as validated SNPs. The mean minor allele frequency of the validated loci was 0.258 (+/- 0.141). Of the 1514 contigs from which validated SNPs were selected, 31% have a significant blast hit. For the SNPs predicted to occur in coding regions (141), we determined that 36% (51) are non-synonymous. Many loci (1033 SNPs; 64%) are polymorphic in all populations tested. However a small number of SNPs (184) that are polymorphic in the Western Atlantic were monomorphic in fish tested from three European populations. A preliminary linkage map has been constructed with 23 major linkage groups and 924 mapped SNPs. CONCLUSIONS: These SNPs represent powerful tools to accelerate the genetic improvement of cod aquaculture. They have been used to build a genetic linkage map that can be applied to quantitative trait locus (QTL) discovery. Since these SNPs were generated from ESTs, they are linked to specific genes. Genes that map within QTL intervals can be prioritized for testing to determine whether they contribute to observed phenotypes.

Mutations in the UBIAD1 gene, encoding a potential prenyltransferase, are causal for Schnyder crystalline corneal dystrophy.

Schnyder crystalline corneal dystrophy (SCCD, MIM 121800) is a rare autosomal dominant disease characterized by progressive opacification of the cornea resulting from the local accumulation of lipids, and associated in some cases with systemic dyslipidemia. Although previous studies of the genetics of SCCD have localized the defective gene to a 1.58 Mbp interval on chromosome 1p, exhaustive sequencing of positional candidate genes has thus far failed to reveal causal mutations. We have ascertained a large multigenerational family in Nova Scotia affected with SCCD in which we have confirmed linkage to the same general area of chromosome 1. Intensive fine mapping in our family revealed a 1.3 Mbp candidate interval overlapping that previously reported. Sequencing of genes in our interval led to the identification of five putative causal mutations in gene UBIAD1, in our family as well as in four other small families of various geographic origins. UBIAD1 encodes a potential prenyltransferase, and is reported to interact physically with apolipoprotein E. UBIAD1 may play a direct role in intracellular cholesterol biochemistry, or may prenylate other proteins regulating cholesterol transport and storage.