Search

The sequencing of the genome of Plasmodium falciparum promises to revolutionize the way in which malaria research will be carried out. Beyond simple gene discovery, the genome sequence will facilitate the comprehensive determination of the parasite?s gene expression during its developmental phases, pathology, and in response to environmental variables, such as drug treatment and host genetic background. This article reviews the current status of the P. falciparum genome sequencing project and the unique insights it has generated. We also summarize the application of bioinformatics and analytical tools that have been developed for functional genomics. The aim of these activities is the rational, information-based identification of new therapeutic strategies and targets, based on a thorough insight into the biology of Plasmodium spp.

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.

Author(s): Malcolm J. Gardner (corresponding author) [1]; Neil Hall [2]; Eula Fung [3]; Owen White [1]; Matthew Berriman [2]; Richard W. Hyman [3]; Jane M. Carlton [1]; Arnab Pain [2]; [...]

Author(s): Hernando A. del Portillo (corresponding author) [1]; Carmen Fernandez-Becerra [1]; Sharen Bowman [2]; Karen Oliver [2]; Martin Preuss [3]; Cecilia P. Sanchez [3]; Nick K. Schneider [3]; Juan M. [...]

Abstract 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.

Abstract 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 ≥ 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.

Abstract 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°C elevation in temperature, and then allowed to recover for 24 hours at the original temperature (i.e. 10°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 HSP90α 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.

Abstract 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 α Hb gene family and five to the β 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. β Hb genes displayed more polymorphism than α Hb genes. Two major allele types (β1A and β1B) that differ by two linked non-synonymous substitutions (Met55Val and Lys62Ala) were found in the β1 Hb gene, and the distribution of these β1A and β1B 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 β1 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.

Abstract 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 ~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.

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. [ABSTRACT FROM AUTHOR]

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 onethird 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. [ABSTRACT FROM AUTHOR]

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 = 0.36) to eastern ( H = 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. [ABSTRACT FROM AUTHOR]

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 ≥ 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. [ABSTRACT FROM AUTHOR]

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°C elevation in temperature, and then allowed to recover for 24 hours at the original temperature (i.e. 10°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 HSP90a and HSP70-1 (synonym: HSPA1A) mRNA expression showed significant upregulation 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. [ABSTRACT FROM AUTHOR]

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. [ABSTRACT FROM AUTHOR]

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 ~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 GCrich 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. [ABSTRACT FROM AUTHOR]

Trypanosoma brucei evades the immune system by switching between Variant Surface Glycoprotein (VSG) genes. The active VSG gene is transcribed in one of approximately 20 telomeric expression sites (ESs). It has been postulated that ES polymorphism plays a role in host adaptation. To gain more insight into ES architecture, we have determined the complete sequence of Bacterial Artificial Chromosomes (BACs) containing DNA from three ESs and their flanking regions. There was variation in the order and number of ES-associated genes (ESAGs). ESAGs 6 and 7, encoding transferrin receptor subunits, are the only ESAGs with functional copies in every ES that has been sequenced until now. A BAC clone containing the VO2 ES sequences comprised approximately half of a 330 kb ‘intermediate’ chromosome. The extensive similarity between this intermediate chromosome and the left telomere of T. brucei 927 chromosome I, suggests that this previously uncharacterised intermediate size class of chromosomes could have arisen from breakage of megabase chromosomes. Unexpected conservation of sequences, including pseudogenes, indicates that the multiple ESs could have arisen through a relatively recent amplification of a single ES. [Copyright &y& Elsevier]

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.

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.

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.

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.

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.

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.