The relationship between daily behavior, hormones, and a color dimorphism in a seabird under natural continuous light.

The predictable oscillation between the light of day and the dark of night across the diel cycle is a powerful selective force that has resulted in anticipatory mechanisms in nearly all taxa. At polar latitude, however, this oscillation becomes highly attenuated during the continuous light of polar day during summer. A general understanding of how animals keep time under these conditions is poorly understood. We tested the hypothesis that the common murre (a seabird, Uria aalge) can use melatonin and corticosterone, hormones associated with timekeeping, to track the diel cycle despite continuous light. We also tested the assumption that common murres breeding during polar summer schedule their colony attendance by time of day and sex, as they do at subpolar latitude. In the Atlantic population, common murres have a plumage color dimorphism associated with fitness-related traits, and we investigated the relationship of this dimorphism with colony attendance, melatonin, and corticosterone. The common murres did not schedule their attendance behavior by time of day or sex, yet they had higher concentrations of melatonin and, to a more limited extent, corticosterone during "night" than "day". Melatonin also linked to behavioral state. The two color morphs tended to have different colony-attendance behavior and melatonin concentrations, lending support for balancing selection maintaining the plumage dimorphism. In common murres, melatonin can signal time of day despite continuous light, and the limited diel variation of corticosterone contributes to the mounting evidence that polar-adapted birds and mammals require little or no diel variation in circulating glucocorticoids during polar day.

Complete loss of the MHC II pathway in an anglerfish, Lophius piscatorius.

Genome studies in fish provide evidence for the adaptability of the vertebrate immune system, revealing alternative immune strategies. The reported absence of the major compatibility complex (MHC) class II pathway components in certain species of pipefish (genus Syngnathus) and cod-like fishes (order Gadiformes) is of particular interest. The MHC II pathway is responsible for immunization and defence against extracellular threats through the presentation of exogenous peptides to T helper cells. Here, we demonstrate the absence of all genes encoding MHC II components (CD4, CD74 A/B, and both classical and non-classical MHC II α/ß) in the genome of an anglerfish, Lophius piscatorius, indicating loss of the MHC II pathway. By contrast, it has previously been reported that another anglerfish, Antennarius striatus, retains all MHC II genes, placing the loss of MHC II in the Lophius clade to their most recent common ancestor. In the three taxa where MHC II loss has occurred, the gene loss has been restricted to four or five core MHC II components, suggesting that, in teleosts, only these genes have functions that are restricted to the MHC II pathway.

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.

An evolutionary preserved intergenic spacer in gadiform mitogenomes generates a long noncoding RNA.

BACKGROUND: Vertebrate mitogenomes are economically organized and usually lack intergenic sequences other than the control region. Intergenic spacers located between the tRNA(Thr) and tRNA(Pro) genes ("T-P spacers") have been observed in several taxa, including gadiform species, but information about their biological roles and putative functions is still lacking. RESULTS: Sequence characterization of the complete European hake Merluccius merluccius mitogenome identified a complex T-P spacer ranging in size from 223-532 bp. Further analyses of 32 gadiform species, representing 8 families and 28 genera, revealed the evolutionary preserved presence of T-P spacers across all taxa. Molecular complexity of the T-P spacers was found to be coherent with the phylogenetic relationships, supporting a common ancestral origin and gain of function during codfish evolution. Intraspecific variation of T-P spacer sequences was assessed in 225 Atlantic cod specimens and revealed 26 haplotypes. Pyrosequencing data representing the mito-transcriptome poly (A) fraction in Atlantic cod identified an abundant H-strand specific long noncoding RNA of about 375 nt. The T-P spacer corresponded to the 5' part of this transcript, which terminated within the control region in a tail-to-tail configuration with the L-strand specific transcript (the 7S RNA). CONCLUSIONS: The T-P spacer is inferred to be evolutionary preserved in gadiform mitogenomes due to gain of function through a long noncoding RNA. We suggest that the T-P spacer adds stability to the H-strand specific long noncoding RNA by forming stable hairpin structures and additional protein binding sites.

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.

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.

Blood plasma clinical-chemical parameters as biomarker endpoints for organohalogen contaminant exposure in Norwegian raptor nestlings.

Raptors are exposed to biomagnifying and toxic organohalogenated compounds (OHCs) such as organochlorines, brominated flame retardants and perfluorinated compounds. To investigate how OHC exposure may affect biochemical pathways we collected blood plasma from Norwegian northern goshawk (n=56), golden eagle (n=12) and white-tailed eagle (n=36) nestlings during three consecutive breeding seasons. We found that blood plasma concentrations of calcium, sodium, creatinine, cholesterol, albumin, total protein, urea, inorganic phosphate, protein:creatinine, urea:creatinine and uric acid:creatinine ratios and liver enzymes ALKP and ALAT were positively correlated to PCBs, chlordanes, p,p'-DDE, HCB, PFCs and/or PBDEs. Total bilirubin and glucose were negatively correlated to PCBs while magnesium and potassium were negatively correlated to HCB and p,p'-DDE. In addition, protein:creatinine and ALAT were also negatively correlated to PCBs and PFCs, respectively. The most significant relationships were found for the highly contaminated northern goshawks and white-tailed eagles. The statistical relationships between OHCs and BCCPs indicate that biochemical pathways could be influenced while it is uncertain if such changes have any health effects. The OHC concentrations were below concentrations causing reproductive toxicity in adults of other raptor species but similar to those of concern for endocrine disruption of thyroid hormones in e.g., bald eagles.

Exploring the potential of mRNA for taxonomic delineation of marine benthic eukaryotes.

Direct sequencing of mRNA isolated from environmental samples has been commonly used to analyze the functional activity of ambient communities and, occasionally used for taxonomic identification. Here we assess the viability of using mRNA for investigating the species composition of marine benthic eukaryotes. Total RNA was extracted from sediments sampled close to fish farms and mRNA was sequenced after poly-A enrichment on an Illumina MiSeq sequencer. We investigated the origin of both raw reads and assembled contigs by aligning them to the NCBI non-redundant (nr/nt) nucleotide database. Although sequences were predominantly of eukaryotic origin, the analyses were complicated both by experimental and database artefacts. These issues were addressed by applying filtering procedures that removed the majority of ambiguous sequences from downstream analyses. These processes resulted in a set of 436 high-confidence contigs, the vast majority of which mapped to benthic organisms. Our alignments were dominated by annelids, consistent with burrowing groups found in a parallel morphological analysis. This study shows that it is possible to obtain adequate taxonomic information from the RNA of an eukaryotic community from a limited sample at a moderate cost, demonstrates how both laboratory and in silico artefacts can be overcome through appropriate bioinformatic procedures, and finally highlights some of the drawbacks and caveats of using NCBI as a reference database for such a dataset.

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 tRNA(Ser 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.

Relationships between organohalogen contaminants and blood plasma clinical-chemical parameters in chicks of three raptor species from Northern Norway.

Organohalogen contaminants (OHCs) may affect various physiological parameters in birds including blood chemistry. We therefore examined blood plasma clinical-chemical parameters and OHCs in golden eagle, white-tailed eagle and goshawk chicks from Northern Norway. Correlation analyses on pooled data showed that alkaline phosphatase (ALKP), glucose and creatinine were significantly negatively correlated to various OHCs (all: p<0.05; r: -0.43 to -0.55; n=23), while alanine aminotransferase (ALAT), total protein, cholesterol, uric acid, total bilirubin, ratios protein:creatinine and uric acid:creatinine were significantly positively correlated to various OHCs (all: p<0.05; r: 0.43-0.96). Based on these relationships, we suggest that the OHC concentrations found in certain raptor chicks of Northern Scandinavia may impact blood plasma biochemistry in a way that indicates impacts on liver, kidney, bone, endocrinology and metabolism. In order to elaborate further on these relationships and mechanisms, we recommend that a larger study should take place in the near future.

Correction: On the genus Crossaster (Echinodermata: Asteroidea) and its distribution.

[This corrects the article DOI: 10.1371/journal.pone.0227223.].

On the genus Crossaster (Echinodermata: Asteroidea) and its distribution.

Several starfish (Echinodermata, Asteroidea) are keystone species of marine ecosystems, but some of the species are difficult to identify using morphological criteria only. The common sunstar, Crossaster papposus (Linnaeus, 1767), is a conspicuous species with a wide circumboreal distribution. In 1900, a closely similar species, C. squamatus (Döderlein, 1900) was described from the NE Atlantic Ocean, but subsequent authors have differed in their views on whether this is a valid taxon or rather an ecotype associated with temperature variations. We assessed the differentiating morphological characters of specimens from Norwegian and Greenland waters identified as C. papposus and C. squamatus and compared their distributions in the NE Atlantic as inferred from research cruises. The field data show that C. papposus is found mainly in temperate and shallow waters, whereas C. squamatus resides on the shelf-break in colder, mixed water masses. Intraspecific diversity and interspecific genetic differentiation of the two putative species, and their phylogenetic relationships to several Crossaster congeners worldwide, were explored using mitochondrial and nuclear DNA sequences. The molecular evidence suggests that C. papposus is the more diverse and geographically structured taxon, in line with its wide distribution. C. papposus and C. squamatus are closely related, yet clearly distinct taxa, while C. papposus and C. multispinus H.L. Clark, 1916, the latter from the South Pacific Ocean, are closely related, possibly sister taxa.

Succession of embryonic and the intestinal bacterial communities of Atlantic salmon (Salmo salar) reveals stage-specific microbial signatures.

Host-associated microbiota undergoes a continuous transition, from the birth to adulthood of the host. These developmental stage-related transitions could lead to specific microbial signatures that could impact the host biological processes. In this study, the succession of early-life and intestinal bacterial communities of Atlantic salmon (starting from embryonic stages to 80-week post hatch; wph) was studied using amplicon sequencing of 16S rRNA. Stage-specific bacterial community compositions and the progressive transitions of the communities were evident in both the early life and the intestine. The embryonic communities showed lower richness and diversity (Shannon and PD whole tree) compared to the hatchlings. A marked transition of the intestinal communities also occurred during the development; Proteobacteria were dominant in the early stages (both embryonic and intestinal), though the abundant genera under this phylum were stage-specific. Firmicutes were the most abundant group in the intestine of late freshwater; Weissella being the dominant genus at 20 wph and Anaerofilum at 62 wph. Proteobacteria regained its dominance after the fish entered seawater. Furthermore, LEfSe analysis identified genera under the above - mentioned phyla that are significant features of specific stages. The environmental (water) bacterial community was significantly different from that of the fish, indicating that the host is a determinant of microbial assemblage. Overall the study demonstrated the community dynamics during the development of Atlantic salmon.

Deep-water sea anemone with a two-chromosome mitochondrial genome.

Mitochondrial genome organization of sea anemones appears conserved among species and families, and is represented by a single circular DNA molecule of 17 to 21 kb. The mitochondrial gene content corresponds to the same 13 protein components of the oxidative phosphorylation (OxPhos) system as in vertebrates. Hallmarks, however, include a highly reduced tRNA gene repertoire and the presence of autocatalytic group I introns. Here we demonstrate that the mitochondrial genome of the deep-water sea anemone Protanthea simplex deviates significantly from that of other known sea anemones. The P. simplex mitochondrial genome contains a heavily scrambled order of genes that are coded on both DNA strands and organized along two circular mito-chromosomes, MCh-I and MCh-II. We found MCh-I to be representative of the prototypic sea anemone mitochondrial genome, encoding 12 OxPhos proteins, two ribosomal RNAs, two transfer RNAs, and a group I intron. In contrast, MCh-II was found to be a laterally transferred plasmid-like DNA carrying the conserved cytochrome oxidase II gene and a second allele of the small subunit ribosomal RNA gene.

A mitochondrial long noncoding RNA in atlantic cod harbors complex heteroplasmic tandem repeat motifs.

A heteroplasmic tandem repeat (HTR) array occupies 100 to 300 bp of the mitochondrial DNA control region in the Atlantic cod, and recently we noted that the repeat appeared integrated in a polyadenylated mitochondrial long noncoding RNA. Here we provide a more detailed analysis of the mitochondrial HTR in the mitochondrial genome of 134 Atlantic cod specimens. We report all specimens to harbor mitochondrial HTRs in the control region, and identified 26 distinct variants among the 402 repeat motifs assessed. Whereas most specimens contained HTR profiles of 2-5 copies consisting of the same 40-bp motif, 22 specimens showed compound HTR arrays of at least two types of motifs present in the same mitochondrial DNA molecule. We found HTR profiles to be highly conserved between different tissue types of a single individual, and strictly maternally inherited in a mating experiment between parental Atlantic cod expressing different HTR profiles and array motifs. We conclude that mitochondrial heteroplasmy in the control region is very common in Atlantic cod, and results in length heterogenity of the long noncoding RNA lncCR-H.

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

BACKGROUND: Halibuts 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. RESULTS: About 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. CONCLUSION: The 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.

Mitochondrial genome variation of Atlantic cod.

OBJECTIVE: The objective of this study was to analyse intraspecific sequence variation of Atlantic cod mitochondrial DNA, based on a comprehensive collection of completely sequenced mitochondrial genomes. RESULTS: We determined the complete mitochondrial DNA sequence of 124 cod specimens from the eastern and western part of the species' distribution range in the North Atlantic Ocean. All specimens harboured a unique mitochondrial DNA haplotype. Nine hundred and fifty-two polymorphic sites were identified, including 109 non-synonymous sites within protein coding regions. Eighteen variable sites were identified as indels, exclusively distributed in structural RNA genes and non-coding regions. Phylogeographic analyses based on 156 available cod mitochondrial genomes did not reveal a clear structure. There was a lack of mitochondrial genetic differentiation between two ecotypes of cod in the eastern North Atlantic, but eastern and western cod were differentiated and mitochondrial genome diversity was higher in the eastern than the western Atlantic, suggesting deviating population histories. The geographic distribution of mitochondrial genome variation seems to be governed by demographic processes and gene flow among ecotypes that are otherwise characterized by localized genomic divergence associated with chromosomal inversions.

A paretic condition in an Anaplasma phagocytophilum infected roe deer calf.

This paper describes a case of Anaplasma phagocytophilum infection in a roe deer (Capreolus capreolus) calf in Norway. The calf was found deserted, paretic, and heavily infested with Ixodes ricinus ticks. It was euthanized and investigated postmortem. Anaplasma phagocytophilum was detected in several tissues by polymerase chain reaction (PCR) and 16S rRNA sequence analyses. Analyses for Borrelia burgdorferi sensu lato and tick-borne encephalitis (TBE) virus infections were negative. This is the first report of a possible paretic condition in A. phagocytophilum infected roe deer.

The stress hormone corticosterone in a marine top predator reflects short-term changes in food availability.

In many seabird studies, single annual proxies of prey abundance have been used to explain variability in breeding performance, but much more important is probably the timing of prey availability relative to the breeding season when energy demand is at a maximum. Until now, intraseasonal variation in prey availability has been difficult to quantify in seabirds. Using a state-of-the-art ocean drift model of larval cod Gadus morhua, an important constituent of the diet of common guillemots Uria aalge in the southwestern Barents Sea, we were able to show clear, short-term correlations between food availability and measurements of the stress hormone corticosterone (CORT) in parental guillemots over a 3-year period (2009-2011). The model allowed the extraction of abundance and size of cod larvae with very high spatial (4 km) and temporal resolutions (1 day) and showed that cod larvae from adjacent northern spawning grounds in Norway were always available near the guillemot breeding colony while those from more distant southerly spawning grounds were less frequent, but larger. The latter arrived in waves whose magnitude and timing, and thus overlap with the guillemot breeding season, varied between years. CORT levels in adult guillemots were lower in birds caught after a week with high frequencies of southern cod larvae. This pattern was restricted to the two years (2009 and 2010) in which southern larvae arrived before the end of the guillemot breeding season. Any such pattern was masked in 2011 by already exceptionally high numbers of cod larvae in the region throughout chick-rearing period. The findings suggest that CORT levels in breeding birds increase when the arrival of southern sizable larvae does not match the period of peak energy requirements during breeding.

Mitogenome sequence variation in migratory and stationary ecotypes of North-east Atlantic cod.

Sequencing of mitochondrial gene fragments from specimens representing a wide range of geographical locations has indicated limited population structuring in Atlantic cod (Gadus morhua). We recently performed whole genome analysis based on next-generation sequencing of two pooled ecotype samples representing offshore migratory and inshore stationary cod from the North-east Atlantic Ocean. Here we report molecular features and variability of the 16.7kb mitogenome component that was collected from the datasets. These sequences represented more than 25 times coverage of each individual and more than 1100 times coverage of each ecotype sample. We estimated the mitogenome to have evolved 14 times more rapidly than the nuclear genome. Among the 365 single nucleotide polymorphism (SNP) sites identified, 121 were shared between ecotypes, and 151 and 93 were private within the migratory and stationary cod, respectively. We found 323 SNPs to be located in protein coding genes, of which 29 were non-synonymous. One synonymous site in ND2 was likely to be under positive selection. FST measurements indicated weak differentiation in ND1 and ND2 between ecotypes. We conclude that the Atlantic cod mitogenome and the nuclear genome apparently evolved by distinct evolutionary constraints, and that the reproductive isolation observed from whole genome analysis was not visible in the mtDNA sequences.