Long distance linkage disequilibrium and limited hybridization suggest cryptic speciation in atlantic cod.

Hybrid zones provide unprecedented opportunity for the study of the evolution of reproductive isolation, and the extent of hybridization across individuals and genomes can illuminate the degree of isolation. We examine patterns of interchromosomal linkage disequilibrium (ILD) and the presence of hybridization in Atlantic cod, Gadus morhua, in previously identified hybrid zones in the North Atlantic. Here, previously identified clinal loci were mapped to the cod genome with most (∼70%) occurring in or associated with (<5 kb) coding regions representing a diverse array of possible functions and pathways. Despite the observation that clinal loci were distributed across three linkage groups, elevated ILD was observed among all groups of clinal loci and strongest in comparisons involving a region of low recombination along linkage group 7. Evidence of ILD supports a hypothesis of divergence hitchhiking transitioning to genome hitchhiking consistent with reproductive isolation. This hypothesis is supported by Bayesian characterization of hybrid classes present and we find evidence of common F1 hybrids in several regions consistent with frequent interbreeding, yet little evidence of F2 or backcrossed individuals. This work suggests that significant barriers to hybridization and introgression exist among these co-occurring groups of cod either through strong selection against hybrid individuals, or genetic incompatibility and intrinsic barriers to hybridization. In either case, the presence of strong clinal trends, and little gene flow despite extensive hybridization supports a hypothesis of reproductive isolation and cryptic speciation in Atlantic cod. Further work is required to test the degree and nature of reproductive isolation in this species.

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.

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.

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.

Impact of asymptomatic nodavirus carrier state and intraperitoneal viral mimic injection on brain transcript expression in Atlantic cod (Gadus morhua).

Nodaviruses and other RNA viruses have a profoundly negative impact on the global aquaculture industry. Nodaviruses target nervous tissue causing viral nervous necrosis, a disease characterized by neurological damage, swimming abnormalities, and morbidity. This study used functional genomic techniques to study the Atlantic cod (Gadus morhua) brain transcript expression responses to asymptomatic high nodavirus carrier state and intraperitoneal injection of polyriboinosinic polyribocytidylic acid (pIC). Reciprocal suppression subtractive hybridization (SSH) cDNA libraries enriched for virus-responsive brain transcripts were constructed and characterized. We generated 1,938 expressed sequence tags (ESTs) from a forward brain SSH library (enriched for transcripts upregulated by nodavirus and/or pIC) and 1,980 ESTs from a reverse brain SSH library (enriched for transcripts downregulated by nodavirus and/or pIC). To examine the effect of nodavirus carrier state on individual brain gene expression in asymptomatic cod, 27 transcripts of interest were selected for quantitative reverse transcription-polymerase chain reaction (QPCR) studies. Transcripts found to be >10-fold upregulated in individuals with a high nodavirus carrier state relative to those in a no/low nodavirus carrier state were identified as ISG15, IL8, DHX58 (alias LGP2), ZNFX1, RSAD2 (alias viperin), and SACS (sacsin, alias spastic ataxia of Charlevoix-Saguenay). These and other SSH-identified transcripts were also found by QPCR to be significantly (P < 0.05) upregulated by pIC compared with saline-injected controls within 72 h of injection. Several transcripts identified in the reverse SSH library, including two putative ubiquitination pathway members (HERC4 and SUMO2), were found to be significantly (P < 0.05) downregulated in individuals with a high nodavirus carrier state. Our data shows that Atlantic cod brains have a strong interferon pathway response to asymptomatic high nodavirus carrier state and that many interferon pathway and other immune relevant transcripts are significantly induced in brain by both nodavirus and pIC.

Atlantic cod (Gadus morhua) CC chemokines: Diversity and expression analysis.

Chemokines are a large, diverse group of small cytokines that can be classified into several families, including the CC chemokines that are characterized by two adjacent cysteines near their amino terminus. CC chemokines play a pivotal role in host defense mechanisms by inducing leukocyte chemotaxis under physiological and inflammatory conditions. Analysis of CC chemokines from teleost fishes indicates that the number of CC chemokine genes and their tissue expression patterns vary largely in this group of vertebrates. Here we describe 32 distinct CC chemokine sequences from Atlantic cod (Gadus morhua) identified by analysis of approximately 206,000 ESTs. Phylogenetic analysis of Atlantic cod CC chemokines placed these sequences in seven clusters, most likely resulting from species-specific gene duplications, and two unique sequences; 12 of these CC chemokines, including at least one member of each cluster, were analyzed by QPCR using four immune-related tissues (head kidney, liver, spleen and blood) obtained from unstimulated, polyriboinosinic polyribocytidylic acid (pIC)-stimulated and formalin-killed atypical Aeromonas salmonicida-stimulated individuals. EST abundance and QPCR analysis indicate that the expression of closely related CC chemokines GmSCYA101 and GmSCYA102, GmSCYA108 and GmSCYA109 or GmSCYA122 and GmSCYA124 can be highly tissue-specific despite substantial sequence identity. Stimulation with the viral mimic pIC or formalin-killed atypical A. salmonicida resulted in increased expression of most of the CC chemokines, indicating that they can be regarded as either inducible (inflammatory) or dual-function rather than constitutive (homeostatic). Tissue specificity, and the level of induction, varied broadly; for example, GmSCYA123 was at least 4-fold up-regulated by both inducers in all tissues analyzed, whereas pIC increased the expression of GmSCYA124 in liver over 1500 times.

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.

Heat-shock responsive genes identified and validated in Atlantic cod (Gadus morhua) liver, head kidney and skeletal muscle using genomic techniques.

BACKGROUND: Daily and seasonal changes in temperature are challenges that fish within aquaculture settings cannot completely avoid, and are known to elicit complex organismal and cellular stress responses. We conducted a large-scale gene discovery and transcript expression study in order to better understand the genes that are potentially involved in the physiological and cellular aspects of stress caused by heat-shock. We used suppression subtractive hybridization (SSH) cDNA library construction and characterization to identify transcripts that were dysregulated by heat-shock in liver, skeletal muscle and head kidney of Atlantic cod. These tissues were selected due to their roles in metabolic regulation, locomotion and growth, and immune function, respectively. Fish were exposed for 3 hours to an 8 degrees C elevation in temperature, and then allowed to recover for 24 hours at the original temperature (i.e. 10 degrees C). Tissue samples obtained before heat-shock (BHS), at the cessation of heat-shock (CS), and 3, 12, and 24 hours after the cessation of heat-shock (ACS), were used for reciprocal SSH library construction and quantitative reverse transcription - polymerase chain reaction (QPCR) analysis of gene expression using samples from a group that was transferred but not heat-shocked (CT) as controls. RESULTS: We sequenced and characterized 4394 ESTs (1524 from liver, 1451 from head kidney and 1419 from skeletal muscle) from three "forward subtracted" libraries (enriched for genes up-regulated by heat-shock) and 1586 from the liver "reverse subtracted" library (enriched for genes down-regulated by heat-shock), for a total of 5980 ESTs. Several cDNAs encoding putative chaperones belonging to the heat-shock protein (HSP) family were found in these libraries, and "protein folding" was among the gene ontology (GO) terms with the highest proportion in the libraries. QPCR analysis of HSP90alpha and HSP70-1 (synonym: HSPA1A) mRNA expression showed significant up-regulation in all three tissues studied. These transcripts were more than 100-fold up-regulated in liver following heat-shock. We also identified HSP47, GRP78 and GRP94-like transcripts, which were significantly up-regulated in all 3 tissues studied. Toll-like receptor 22 (TLR22) transcript, found in the liver reverse SSH library, was shown by QPCR to be significantly down-regulated in the head kidney after heat-shock. CONCLUSION: Chaperones are an important part of the cellular response to stress, and genes identified in this work may play important roles in resistance to thermal-stress. Moreover, the transcript for one key immune response gene (TLR22) was down-regulated by heat-shock, and this down-regulation may be a component of heat-induced immunosuppression.

Atlantic cod (Gadus morhua) hemoglobin genes: multiplicity and polymorphism.

BACKGROUND: Hemoglobin (Hb) polymorphism, assessed by protein gel electrophoresis, has been used almost exclusively to characterize the genetic structure of Atlantic cod (Gadus morhua) populations and to establish correlations with phenotypic traits such as Hb oxygen binding capacity, temperature tolerance and growth characteristics. The genetic system used to explain the results of gel electrophoresis entails the presence of one polymorphic locus with two major alleles (HbI-1; HbI-2). However, vertebrates have more than one gene encoding Hbs and recent studies have reported that more than one Hb gene is present in Atlantic cod. These observations prompted us to re-evaluate the number of Hb genes expressed in Atlantic cod, and to perform an in depth search for polymorphisms that might produce relevant phenotypes for breeding programs. RESULTS: Analysis of Expressed Sequence Tags (ESTs) led to the identification of nine distinct Hb transcripts; four corresponding to the alpha Hb gene family and five to the beta Hb gene family. To gain insights about the Hb genes encoding these transcripts, genomic sequence data was generated from heterozygous (HbI-1/2) parents and fifteen progeny; five of each HbI type, i.e., HbI-1/1, HbI-1/2 and HbI-2/2. beta Hb genes displayed more polymorphism than alpha Hb genes. Two major allele types (beta1A and beta1B) that differ by two linked non-synonymous substitutions (Met55Val and Lys62Ala) were found in the beta1 Hb gene, and the distribution of these beta1A and beta1B alleles among individuals was congruent with that of the HbI-1 and HbI-2 alleles determined by protein gel electrophoresis. RT-PCR and Q-PCR analysis of the nine Hb genes indicates that all genes are expressed in adult fish, but their level of expression varies greatly; higher expression of almost all Hb genes was found in individuals displaying the HbI-2/2 electrophoretic type. CONCLUSION: This study indicates that more Hb genes are present and expressed in adult Atlantic cod than previously documented. Our finding that nine Hb genes are expressed simultaneously in adult fish suggests that Atlantic cod, similarly to fish such as rainbow trout, carp, and goldfish, might be able to respond to environmental challenges such as chronic hypoxia or long-term changes in temperature by altering the level of expression of these genes. In this context, the role of the non-conservative substitution Lys62Ala found in the beta1 Hb gene, which appears to explain the occurrence of the HbI-1 and HbI-2 alleles described by gel electrophoresis, and which was found to be present in other fish such as eel, emerald rockcod, rainbow trout and moray, requires further investigation.

Identification and analysis of differentially expressed genes in immune tissues of Atlantic cod stimulated with formalin-killed, atypical Aeromonas salmonicida.

Physiological changes, elicited in animal immune tissues by exposure to pathogens, may be studied using functional genomics approaches. We created and characterized reciprocal suppression subtractive hybridization (SSH) cDNA libraries to identify differentially expressed genes in spleen and head kidney tissues of Atlantic cod (Gadus morhua) challenged with intraperitoneal injections of formalin-killed, atypical Aeromonas salmonicida. Of 4,154 ESTs from four cDNA libraries, 10 genes with immune-relevant functional annotations were selected for QPCR studies using individual fish templates to assess biological variability. Genes confirmed by QPCR as upregulated by A. salmonicida included interleukin-1 beta, interleukin-8, a small inducible cytokine, interferon regulatory factor 1 (IRF1), ferritin heavy subunit, cathelicidin, and hepcidin. This study is the first large-scale discovery of bacteria-responsive genes in cod and the first to demonstrate upregulation of IRF1 in fish immune tissues as a result of bacterial antigen stimulation. Given the importance of IRF1 in vertebrate immune responses to viral and bacterial pathogens, the full-length cDNA sequence of Atlantic cod IRF1 was obtained and compared with putative orthologous sequences from other organisms. Functional annotations of assembled SSH library ESTs showed that bacterial antigen stimulation caused changes in many biological processes including chemotaxis, regulation of apoptosis, antimicrobial peptide production, and iron homeostasis. Moreover, differences in spleen and head kidney gene expression responses to the bacterial antigens pointed to a potential role for the cod spleen in blood-borne pathogen clearance. Our data show that Atlantic cod immune tissue responses to bacterial antigens are similar to those seen in other fish species and higher vertebrates.

The genome of Thermosipho africanus TCF52B: lateral genetic connections to the Firmicutes and Archaea.

Lateral gene transfers (LGT) (also called horizontal gene transfers) have been a major force shaping the Thermosipho africanus TCF52B genome, whose sequence we describe here. Firmicutes emerge as the principal LGT partner. Twenty-six percent of phylogenetic trees suggest LGT with this group, while 13% of the open reading frames indicate LGT with Archaea.

Complete sequence and analysis of the mitochondrial genome of Hemiselmis andersenii CCMP644 (Cryptophyceae).

BACKGROUND: Cryptophytes are an enigmatic group of unicellular eukaryotes with plastids derived by secondary (i.e., eukaryote-eukaryote) endosymbiosis. Cryptophytes are unusual in that they possess four genomes-a host cell-derived nuclear and mitochondrial genome and an endosymbiont-derived plastid and 'nucleomorph' genome. The evolutionary origins of the host and endosymbiont components of cryptophyte algae are at present poorly understood. Thus far, a single complete mitochondrial genome sequence has been determined for the cryptophyte Rhodomonas salina. Here, the second complete mitochondrial genome of the cryptophyte alga Hemiselmis andersenii CCMP644 is presented. RESULTS: The H. andersenii mtDNA is 60,553 bp in size and encodes 30 structural RNAs and 36 protein-coding genes, all located on the same strand. A prominent feature of the genome is the presence of a approximately 20 Kbp long intergenic region comprised of numerous tandem and dispersed repeat units of between 22-336 bp. Adjacent to these repeats are 27 copies of palindromic sequences predicted to form stable DNA stem-loop structures. One such stem-loop is located near a GC-rich and GC-poor region and may have a regulatory function in replication or transcription. The H. andersenii mtDNA shares a number of features in common with the genome of the cryptophyte Rhodomonas salina, including general architecture, gene content, and the presence of a large repeat region. However, the H. andersenii mtDNA is devoid of inverted repeats and introns, which are present in R. salina. Comparative analyses of the suite of tRNAs encoded in the two genomes reveal that the H. andersenii mtDNA has lost or converted its original trnK(uuu) gene and possesses a trnS-derived 'trnK(uuu)', which appears unable to produce a functional tRNA. Mitochondrial protein coding gene phylogenies strongly support a variety of previously established eukaryotic groups, but fail to resolve the relationships among higher-order eukaryotic lineages. CONCLUSION: Comparison of the H. andersenii and R. salina mitochondrial genomes reveals a number of cryptophyte-specific genomic features, most notably the presence of a large repeat-rich intergenic region. However, unlike R. salina, the H. andersenii mtDNA does not possess introns and lacks a Lys-tRNA, which is presumably imported from the cytosol.

Functional genomic analysis of the response of Atlantic cod (Gadus morhua) spleen to the viral mimic polyriboinosinic polyribocytidylic acid (pIC).

In order to improve our understanding of how Atlantic cod (Gadus morhua) respond to viruses, we characterized immune-related gene expression in spleen tissues following stimulation with a synthetic double-stranded RNA polyriboinosinic polyribocytidylic acid (pIC). We used reciprocal suppression subtractive hybridization (SSH) cDNA libraries and quantitative RTPCR (QPCR) to identify and quantify pIC-responsive transcripts. A total of 3874 expressed sequence tags (ESTs) were generated from SSH libraries enriched for genes responsive to pIC. Thirteen immune-relevant genes from the libraries were subjected to QPCR. Genes confirmed as up-regulated by pIC included interferon stimulated gene 15, a small inducible cytokine, interferon regulatory factors (1, 7, and 10), MHC class I, viperin, and ATP-dependent helicase LGP2. Alpha-1-microglobulin (bikunin) was down-regulated, suggesting that pIC may suppress the acute phase response. Since the SSH libraries built for this study identified genes involved in the antiviral response, they are important resources for studying the responses of Atlantic cod to viruses. Evidence is provided for the existence of a RIG-I-like RNA helicase viral recognition pathway in Atlantic cod. Taken together, our data show that Atlantic cod can recognize double-stranded RNA and mount a rapid and potent interferon pathway response that is similar to that observed in other fish species and higher vertebrates.

Nucleomorph genome of Hemiselmis andersenii reveals complete intron loss and compaction as a driver of protein structure and function.

Nucleomorphs are the remnant nuclei of algal endosymbionts that took up residence inside a nonphotosynthetic eukaryotic host. The nucleomorphs of cryptophytes and chlorarachniophytes are derived from red and green algal endosymbionts, respectively, and represent a stunning example of convergent evolution: their genomes have independently been reduced and compacted to <1 megabase pairs (Mbp) in size (the smallest nuclear genomes known) and to a similar three-chromosome architecture. The molecular processes underlying genome reduction and compaction in eukaryotes are largely unknown, as is the impact of reduction/compaction on protein structure and function. Here, we present the complete 0.572-Mbp nucleomorph genome of the cryptophyte Hemiselmis andersenii and show that it is completely devoid of spliceosomal introns and genes for splicing RNAs-a case of complete intron loss in a nuclear genome. Comparison of H. andersenii proteins to those encoded in the slightly smaller (0.551-Mbp) nucleomorph genome of another cryptophyte, Guillardia theta, and to their homologs in the unicellular red alga Cyanidioschyzon merolae reveal that (i) cryptophyte nucleomorph genomes encode proteins that are significantly smaller than those in their free-living algal ancestors, and (ii) the smaller, more compact G. theta nucleomorph genome encodes significantly smaller proteins than that of H. andersenii. These results indicate that genome compaction can eliminate both coding and noncoding DNA and, consequently, drive the evolution of protein structure and function. Nucleomorph proteins have the potential to reveal the minimal functional units required for basic eukaryotic cellular processes.

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.

Plastid genome sequence of the cryptophyte alga Rhodomonas salina CCMP1319: lateral transfer of putative DNA replication machinery and a test of chromist plastid phylogeny.

Cryptophytes are a group of unicellular algae with chlorophyll c-containing plastids derived from the uptake of a secondary (i.e., eukaryotic) endosymbiont. Biochemical and molecular data indicate that cryptophyte plastids are derived from red algae, yet the question of whether or not cryptophytes acquired their red algal plastids independent of those in heterokont, haptophyte, and dinoflagellate algae is of long-standing debate. To better understand the origin and evolution of the cryptophyte plastid, we have sequenced the plastid genome of Rhodomonas salina CCMP1319: at 135,854 bp, it is the largest secondary plastid genome characterized thus far. It also possesses interesting features not seen in the distantly related cryptophyte Guillardia theta or in other red secondary plastids, including pseudogenes, introns, and a bacterial-derived gene for the tau/gamma subunit of DNA polymerase III (dnaX), the first time putative DNA replication machinery has been found encoded in any plastid genome. Phylogenetic analyses indicate that dnaX was acquired by lateral gene transfer (LGT) in an ancestor of Rhodomonas, most likely from a firmicute bacterium. A phylogenomic survey revealed no additional cases of LGT, beyond a noncyanobacterial type rpl36 gene similar to that recently characterized in other cryptophytes and haptophytes. Rigorous concatenated analysis of 45 proteins encoded in 15 complete plastid genomes produced trees in which the heterokont, haptophyte, and cryptophyte (i.e., chromist) plastids were monophyletic, and heterokonts and haptophytes were each other's closest relatives. However, statistical support for chromist monophyly disappears when amino acids are recoded according to their chemical properties in order to minimize the impact of composition bias, and a significant fraction of the concatenate appears consistent with a sister-group relationship between cryptophyte and haptophyte plastids.

Retrotransposons and tandem repeat sequences in the nuclear genomes of cryptomonad algae.

The cryptomonads are an enigmatic group of unicellular eukaryotic algae that possess two nuclear genomes, having acquired photosynthesis by the uptake and retention of a eukaryotic algal endosymbiont. The endosymbiont nuclear genome, or nucleomorph, of the cryptomonad Guillardia theta has been completely sequenced: at only 551 kilobases (kb) and with a gene density of approximately 1 gene/kb, it is a model of compaction. In contrast, very little is known about the structure and composition of the cryptomonad host nuclear genome. Here we present the results of two small-scale sequencing surveys of fosmid clone libraries from two distantly related cryptomonads, Rhodomonas salina CCMP1319 and Cryptomonas paramecium CCAP977/2A, corresponding to approximately 150 and approximately 235 kb of sequence, respectively. Very few of the random end sequences determined in this study show similarity to known genes in other eukaryotes, underscoring the considerable evolutionary distance between the cryptomonads and other eukaryotes whose nuclear genomes have been completely sequenced. Using a combination of fosmid clone end-sequencing, Southern hybridizations, and PCR, we demonstrate that Ty3-gypsy long-terminal repeat (LTR) retrotransposons and tandem repeat sequences are a prominent feature of the nuclear genomes of both organisms. The complete sequence of a 30.9-kb genomic fragment from R. salina was found to contain a full-length Ty3-gypsy element with near-identical LTRs and a chromodomain, a protein module suggested to mediate the site-specific integration of the retrotransposon. The discovery of chromodomain-containing retroelements in cryptomonads further expands the known distribution of the so-called chromoviruses across the tree of eukaryotes.