Bård Ove Karlsen

The genome sequence of Atlantic cod reveals a unique immune system

Atlantic cod (Gadus morhua) is a large, cold-adapted teleost that sustains long-standing commercial fisheries and incipient aquaculture. Here we present the genome sequence of Atlantic cod, showing evidence for complex thermal adaptations in its haemoglobin gene cluster and an unusual immune architecture compared to other sequenced vertebrates. The genome assembly was obtained exclusively by 454 sequencing of shotgun and paired-end libraries, and automated annotation identified 22,154 genes. The major histocompatibility complex (MHC) II is a conserved feature of the adaptive immune system of jawed vertebrates, but we show that Atlantic cod has lost the genes for MHC II, CD4 and invariant chain (Ii) that are essential for the function of this pathway. Nevertheless, Atlantic cod is not exceptionally susceptible to disease under natural conditions. We find a highly expanded number of MHC I genes and a unique composition of its Toll-like receptor (TLR) families. This indicates how the Atlantic cod immune system has evolved compensatory mechanisms in both adaptive and innate immunity in the absence of MHC II. These observations affect fundamental assumptions about the evolution of the adaptive immune system and its components in vertebrates.

Large-scale sequence analyses of Atlantic cod.

The Atlantic cod (Gadus morhua) is a key species in the North Atlantic ecosystem and commercial fisheries, with increasing aquacultural production in several countries. A Norwegian effort to sequence the complete 0.9Gbp genome by the 454 pyrosequencing technology has been initiated and is in progress. Here we review recent progress in large-scale sequence analyses of the nuclear genome, the mitochondrial genome and genome-wide microRNA identification in the Atlantic cod. The nuclear genome will be de novo sequenced with 25 times oversampling. A total of 120 mitochondrial genomes, sampled from several locations in the North Atlantic, are being completely sequenced by Sanger technology in a high-throughput pipeline. These sequences will be included in a new database for maternal marker reference of Atlantic cod diversity. High-throughput 454 sequencing, as well as Evolutionary Image Array (EvoArray) informatics, is used to investigate the complete set of expressed microRNAs and corresponding mRNA targets in various developmental stages and tissues. Information about microRNA profiles will be essential in the understanding of transcriptome complexity and regulation. Finally, developments and perspectives of Atlantic cod aquaculture are discussed in the light of next-generation high-throughput sequence technologies.

Halibut mitochondrial genomes contain extensive heteroplasmic tandem repeat arrays involved in DNA recombination

BackgroundHalibuts are commercially important flatfish species confined to the North Pacific and North Atlantic Oceans. We have determined the complete mitochondrial genome sequences of four specimens each of Atlantic halibut (Hippoglossus hippoglossus), Pacific halibut (Hippoglossus stenolepis) and Greenland halibut (Reinhardtius hippoglossoides), and assessed the nucleotide variability within and between species.ResultsAbout 100 variable positions were identified within the four specimens in each halibut species, with the control regions as the most variable parts of the genomes (10 times that of the mitochondrial ribosomal DNA). Due to tandem repeat arrays, the control regions have unusually large sizes compared to most vertebrate mtDNAs. The arrays are highly heteroplasmic in size and consist mainly of different variants of a 61-bp motif. Halibut mitochondrial genomes lacking arrays were also detected.ConclusionThe complexity, distribution, and biological role of the heteroplasmic tandem repeat arrays in halibut mitochondrial control regions are discussed. We conclude that the most plausible explanation for array maintenance includes both the slipped-strand mispairing and DNA recombination mechanisms.

Genomic divergence between the migratory and stationary ecotypes of Atlantic cod

Atlantic cod displays a range of phenotypic and genotypic variations, which includes the differentiation into coastal stationary and offshore migratory types of cod that co‐occur in several parts of its distribution range and are often sympatric on the spawning grounds. Differentiation of these ecotypes may involve both historical separation and adaptation to ecologically distinct environments, the genetic basis of which is now beginning to be unravelled. Genomic analyses based on recent sequencing advances are able to document genomic divergence in more detail and may facilitate the exploration of causes and consequences of genome‐wide patterns. We examined genomic divergence between the stationary and migratory types of cod in the Northeast Atlantic, using next‐generation sequencing of pooled DNA from each of two population samples. Sequence data was mapped to the published cod genome sequence, arranged in more than 6000 scaffolds (611 Mb). We identified 25 divergent scaffolds (26 Mb) with a higher than average gene density, against a backdrop of overall moderate genomic differentiation. Previous findings of localized genomic divergence in three linkage groups were confirmed, including a large (15 Mb) genomic region, which seems to be uniquely involved in the divergence of migratory and stationary cod. The results of the pooled sequencing approach support and extend recent findings based on single‐nucleotide polymorphism markers and suggest a high degree of reproductive isolation between stationary and migratory cod in the North‐east Atlantic.

RNA deep sequencing of the Atlantic cod transcriptome

The Atlantic cod (Gadus morhua) is an emerging aquaculture species. Efforts to develop and characterize its genomic recourses, including draft-grade genome sequencing, have been initiated by the research community. The transcriptome represents the whole complement of RNA transcripts in cells and tissues and reflects the expressed genes at various life stages, tissue types, physiological states, and environmental conditions. We are investigating the Atlantic cod transcriptome by Roche 454, Illumina GA, and ABI SOLiD deep sequencing platforms and corresponding bioinformatics. Both embryonic developmental stages and adult tissues are studied. Here we summarize our recent progress in the analyses of nuclear and mitochondrial polyA mRNAs, non-protein-coding intermediate RNAs, and regulatory microRNAs.

A novel beta-defensin antimicrobial peptide in Atlantic cod with stimulatory effect on phagocytic activity.

A novel defensin antimicrobial peptide gene was identified in Atlantic cod, Gadus morhua. This three exon/two intron defensin gene codes for a peptide precursor consisting of two domains: a signal peptide of 26 amino acids and a mature peptide of 40 residues. The mature cod defensin has six conserved cysteine residues that form 1–5, 2–4 and 3–6 disulphide bridges. This pattern is typical of beta-defensins and this gene was therefore named cod beta-defensin (defb). The tertiary structure of Defb exhibits an α/β fold with one α helix and β1β2β3 sheets. RT-PCR analysis indicated that defb transcripts were present mainly in the swim bladder and peritoneum wall but could also be detected at moderate to low levels in skin, head- and excretory kidneys. In situ hybridisation revealed that defb was specifically expressed by cells located in the swim bladder submucosa and the oocytes. During embryonic development, defb gene transcripts were detectable from the golden eye stage onwards and their expression was restricted to the swim bladder and retina. Defb was differentially expressed in several tissues following antigenic challenge with Vibrio anguillarum, being up-regulated up to 25-fold in head kidney. Recombinant Defb displayed antibacterial activity, with a minimal inhibitory concentration of 0.4–0.8 µM and 25–50 µM against the Gram-(+) bacteria Planococcus citreus and Micrococcus luteus, respectively. In addition, Defb stimulated phagocytic activity of cod head kidney leucocytes in vitro. These findings imply that beta-defensins may play an important role in the innate immune response of Atlantic cod.

Digital Marine Bioprospecting : Mining New Neurotoxin Drug Candidates from the Transcriptomes of Cold-Water Sea Anemones

Marine bioprospecting is the search for new marine bioactive compounds and large-scale screening in extracts represents the traditional approach. Here, we report an alternative complementary protocol, called digital marine bioprospecting, based on deep sequencing of transcriptomes. We sequenced the transcriptomes from the adult polyp stage of two cold-water sea anemones, Bolocera tuediae and Hormathia digitata. We generated approximately 1.1 million quality-filtered sequencing reads by 454 pyrosequencing, which were assembled into approximately 120,000 contigs and 220,000 single reads. Based on annotation and gene ontology analysis we profiled the expressed mRNA transcripts according to known biological processes. As a proof-of-concept we identified polypeptide toxins with a potential blocking activity on sodium and potassium voltage-gated channels from digital transcriptome libraries.

Characterization of mitochondrial mRNAs in codfish reveals unique features compared to mammals

Expression and processing of mitochondrial gene transcripts are fundamental to mitochondrial function, but information from early vertebrates like teleost fishes is essentially lacking. We have analyzed mitogenome sequences of ten codfishes (family Gadidae), and provide complete sequences from three new species (Saithe, Pollack and Blue whiting). Characterization of the mitochondrial mRNAs in Saithe and Atlantic cod identified a set of ten poly(A) transcripts, and six UAA stop codons are generated by posttranscriptional polyadenylation. Structural assessment of poly(A) sites is consistent with an RNaseP cleavage activity 5′ of tRNA acceptor-like stems. COI, ND5 and ND6 mRNAs were found to harbor 3′ UTRs with antisense potential extending into neighboring gene regions. While the 3′ UTR of COI mRNA is complementary to the tRNASer (UCN) and highly similar to that detected in human mitochondria, the ND5 and ND6 3′ UTRs appear more heterogenic. Deep sequencing confirms expression of all mitochondrial mRNAs and rRNAs, and provides information about the precise 5′ ends in mature transcripts. Our study supports an overall evolutionary conservation in mitochondrial RNA processing events among vertebrates, but reveals some unique 5′ and 3′ end characteristics in codfish mRNAs with implications to antisense regulation of gene expression.

Approaching marine bioprospecting in hexacorals by RNA deep sequencing

RNA deep sequencing represents a new complementary approach in marine bioprospecting. Next-generation sequencing platforms have recently been developed for de novo whole transcriptome analysis, small RNA discovery and gene expression profiling. Deep sequencing transcriptomics (sequencing the complete set of cellular transcripts at a specific stage or condition) leads to sequential identification of all expressed genes in a sample. When combined to high-throughput bioinformatics and protein synthesis, RNA deep sequencing represents a new powerful approach in gene product discovery and bioprospecting. Here we summarize recent progress in the analyses of hexacoral transcriptomes with the focus on cold-water sea anemones and related organisms.

Sea anemones possess dynamic mitogenome structures

A notable feature of hexacoral mitogenomes is the presence of complex self-catalytic group I introns. We investigated mitogenome structural variations and evolutionary mechanisms in actiniarian sea anemones based on the complete mitogenome sequence of the cold-water sea anemone species Urticina eques, Bolocera tuediae, Hormathia digitata and Metridium senile, and two isolates of the sub-tropical Aiptasia pulchella. Whole genome sequencing at 50 times coverage of B. tuediae and H. digitata indicated low mtDNA copy number of per haploid nuclear genome and presence of rare haplotypes. A group I intron inserted in ND5 was found to host essential mitochondrial protein genes in all species, and an additional truncated copy of ND5 in B. tuediae. A second group I intron (inserted in COI) that contained a homing endonuclease gene (HEG) was present in all mtDNA examined. Different variants of HEGs were observed, and included expressed elements fused in-frame with upstream exons and free-standing HEGs embedded within the intron. A notable hallmark of HEGs was a high extent of overlap with ribozyme structural elements; the U. eques HEG overlapped with the entire intron. We reconstructed the evolutionary history of the COI intron from insertion at unoccupied cognate sites, through HEG degradation, to intron loss. We also identified a novel insertion element in U. eques that contained two expressed protein-coding genes. An evolutionary analysis of the sea anemone mtDNA genes revealed higher substitution rates in the HEG and the insertion sequence as compared to the other loci, indicating relaxed selective pressures in these elements. We conclude that sea anemone mitogenomes are surprisingly dynamic in structure despite the economical organization and low sequence mutation rate.