Background Polar environments are characterized by extreme seasonal changes in day length, light intensity and spectrum, the extent of sea ice during the winter, and food availability. A key species of the Southern Ocean ecosystem, the Antarctic krill (Euphausia superba) has evolved rhythmic physiological and behavioral mechanisms to adapt to daily and seasonal changes. The molecular organization of the clockwork underlying these biological rhythms is, nevertheless, still only partially understood. Methodology/Principal Findings The genome sequence of the Antarctic krill is not yet available. A normalized cDNA library was produced and pyrosequenced in the attempt to identify large numbers of transcripts. All available E. superba sequences were then assembled to create the most complete existing oligonucleotide microarray platform with a total of 32,217 probes. Gene expression signatures of specimens collected in the Ross Sea at five different time points over a 24-hour cycle were defined, and 1,308 genes differentially expressed were identified. Of the corresponding transcripts, 609 showed a significant sinusoidal expression pattern; about 40% of these exibithed a 24-hour periodicity while the other 60% was characterized by a shorter (about 12-hour) rhythm. We assigned the differentially expressed genes to functional categories and noticed that those concerning translation, proteolysis, energy and metabolic process, redox regulation, visual transduction and stress response, which are most likely related to daily environmental changes, were significantly enriched. Two transcripts of peroxiredoxin, thought to represent the ancestral timekeeping system that evolved about 2.5 billion years ago, were also identified as were two isoforms of the EsRh1 opsin and two novel arrestin1 sequences involved in the visual transduction cascade. Conclusions Our work represents the first characterization of the krill diurnal transcriptome under natural conditions and provides a first insight into the genetic regulation of physiological changes, which occur around the clock during an Antarctic summer day.
References
[1]
Atkinson A, Siegel V, Pakhomov EA, Jessopp MJ, Loeb V (2009) A re-appraisal of the total biomass and annual production of antarctic krill. Deep Sea Research Part I: Oceanographic Research Papers 56: 727–740.
[2]
Atkinson A, Siegel V, Pakhomov EA, Rothery P, Loeb V, et al. (2008) Oceanic circumpolar habitats of antarctic krill. Mar Ecol Prog Ser 362: 1–23.
[3]
Clarke A, Tyler PA (2008) Adult antarctic krill feeding at abyssal depths. Curr Biol 18: 282–285.
[4]
Flores H, van Franeker JA, Siegel V, Haraldsson M, Strass V, et al. (2012) The association of antarctic krill Euphausia superba with the under-ice habitat. PLoS One 7: e31775.
[5]
Gigliotti JC, Jaczynski J, Tou JC (2008) Determination of the nutritional value, protein quality and safety of krill protein concentrate isolated using an isoelectric solubilization/precipitation technique. Food Chem 111: 209–214.
[6]
Atkinson A, Siegel V, Pakhomov E, Rothery P (2004) Long-term decline in krill stock and increase in salps within the southern ocean. Nature 432: 100–103.
[7]
Forcada J, Trathan PN, Reid K, Murphy EJ (2005) The effects of global climate variability in pup production of antartic fur seals. Ecology 86: 2408–2417.
[8]
Trathan PN, Forcada J, Murphy EJ (2007) Environmental forcing and southern ocean marine predator populations: Effects of climate change and variability. Philosophical Transactions of the Royal Society B: Biological Sciences 362: 2351–2365.
[9]
Kawaguchi S, Nicol S, Press AJ (2009) Direct effects of climate change on the antarctic krill fishery. Fish Manage Ecol 16: 424–427.
[10]
Teschke M, Kawaguchi S, Meyer B (2007) Simulated light regimes affect feeding and metabolism of antarctic krill, Euphausia superba. Limnol Oceanogr 52: 1046–1054.
[11]
Meyer B, Auerswald L, Siegel V, Spahi？ S, Pape C, et al. (2010) Seasonal variation in body composition, metabolic activity, feeding, and growth of adult krill Euphausia superba in the lazarev sea. Mar Ecol Prog Ser 398: 1–18.
[12]
Seear PJ, Goodall-Copestake WP, Fleming AH, Rosato E, Tarling GA (2012) Seasonal and spatial influences on gene expression in antarctic krill Euphausia superba. Mar Ecol Prog Ser 467: 61–75.
[13]
Strauss J, Dircksen H (2010) Circadian clocks in crustaceans: identified neuronal and cellular systems. Front Biosci 15: 1040–1074.
[14]
Teschke M, Wendt S, Kawaguchi S, Kramer A, Meyer B (2011) A circadian clock in antarctic krill: an endogenous timing system governs metabolic output rhythms in the euphausid species Euphausia superba. PLoS One 6: e26090.
[15]
Gaten E, Tarling G, Dowse H, Kyriacou C, Rosato E (2008) Is vertical migration in antarctic krill (Euphausia superba) influenced by an underlying circadian rhythm? J Genet 87: 473–483.
[16]
Godlewska M (1996) Vertical migrations of krill (Euphausia superba DANA). Polskie Archiwum Hydrobiologii 43: 9–63.
[17]
Tarling GA, Johnson ML (2006) Satiation gives krill that sinking feeling. Curr Biol 16: R83–4.
[18]
Aguzzi J (2010) Company JB (2010) Chronobiology of deep-water decapod crustaceans on continental margins. Adv Mar Biol 58: 155–225.
[19]
Naylor E (2006) Orientation and navigation in coastal and estuarine zooplankton. Mar Freshwat Behav Physiol 39: 13–24.
[20]
Tessmar-Raible K, Raible F, Arboleda E (2011) Another place, another timer: marine species and the rhythms of life. Bioessays 33: 165–172.
[21]
Mazzotta GM, De Pittà C, Benna C, Tosatto SC, Lanfranchi G, et al. (2010) A cry from the krill. Chronobiol Int 27: 425–445.
[22]
De Pittà C, Bertolucci C, Mazzotta GM, Bernante F, Rizzo G, et al. (2008) Systematic sequencing of mRNA from the antarctic krill (Euphausia superba) and first tissue specific transcriptional signature. BMC Genomics 9: 45.
[23]
Seear PJ, Tarling GA, Burns G, Goodall-Copestake WP, Gaten E, et al. (2010) Differential gene expression during the moult cycle of antarctic krill (Euphausia superba). BMC Genomics 11: 582.
[24]
Clark MS, Thorne MA, Toullec JY, Meng Y, Guan LL, et al. (2011) Antarctic krill 454 pyrosequencing reveals chaperone and stress transcriptome. PLoS One 6: e15919.
[25]
Mundry M, Bornberg-Bauer E, Sammeth M, Feulner PG (2012) Evaluating characteristics of de novo assembly software on 454 transcriptome data: A simulation approach. PLoS One 7: e31410.
[26]
Jung H, Lyons RE, Dinh H, Hurwood DA, McWilliam S, et al. (2011) Transcriptomics of a giant freshwater prawn (Macrobrachium rosenbergii): De novo assembly, annotation and marker discovery. PLoS One 6: e27938.
[27]
Ma K, Qiu G, Feng J, Li J (2012) Transcriptome analysis of the oriental river prawn, Macrobrachium nipponense using 454 pyrosequencing for discovery of genes and markers. PLoS One 7: e39727.
[28]
Craft JA, Gilbert JA, Temperton B, Dempsey KE, Ashelford K, et al. (2010) Pyrosequencing of Mytilus galloprovincialis cDNAs: Tissue-specific expression patterns. PLoS One 5: e8875.
[29]
Clark MS, Thorne MA, Vieira FA, Cardoso JC, Power DM, et al. (2010) Insights into shell deposition in the antarctic bivalve Laternula elliptica: Gene discovery in the mantle transcriptome using 454 pyrosequencing. BMC Genomics 11: 362.
[30]
Meyer E, Aglyamova GV, Wang S, Buchanan-Carter J, Abrego D, et al. (2009) Sequencing and de novo analysis of a coral larval transcriptome using 454 GSFlx. BMC Genomics 10: 219.
[31]
Milan M, Coppe A, Reinhardt R, Cancela LM, Leite RB, et al. (2011) Transcriptome sequencing and microarray development for the manila clam, Ruditapes philippinarum: Genomic tools for environmental monitoring. BMC Genomics 12: 234.
[32]
Storey JD, Tibshirani R (2003) Statistical significance for genomewide studies. Proc Natl Acad Sci U S A 100: 9440–9445.
[33]
Doherty CJ, Kay SA (2010) Circadian control of global gene expression patterns. Annu Rev Genet 44: 419–444.
[34]
Yang R, Su Z (2010) Analyzing circadian expression data by harmonic regression based on autoregressive spectral estimation. Bioinformatics 26: i168–74.
[35]
Chudova D, Ihler A, Lin KK, Andersen B, Smyth P (2009) Bayesian detection of non-sinusoidal periodic patterns in circadian expression data. Bioinformatics 25: 3114–3120.
[36]
Ceriani MF, Hogenesch JB, Yanovsky M, Panda S, Straume M, et al. (2002) Genome-wide expression analysis in Drosophila reveals genes controlling circadian behavior. J Neurosci 22: 9305–9319.
[37]
Claridge-Chang A, Wijnen H, Naef F, Boothroyd C, Rajewsky N, et al. (2001) Circadian regulation of gene expression systems in the Drosophila head. Neuron 32: 657–671.
[38]
Panda S, Antoch MP, Miller BH, Su AI, Schook AB, et al. (2002) Coordinated transcription of key pathways in the mouse by the circadian clock. Cell 109: 307–320.
[39]
Weger BD, Sahinbas M, Otto GW, Mracek P, Armant O, et al. (2011) The light responsive transcriptome of the zebrafish: Function and regulation. PLoS One 6: e17080.
[40]
Fukuyama T, Doi M, Matsuo M, Nishinaga H, Miyake S, et al. (2008) Circadian expression of 86- and 84-kDa heat shock proteins in the mouse suprachiasmatic nucleus. Biomed Res 29: 93–98.
[41]
Doi M, Okano T, Yujnovsky I, Sassone-Corsi P, Fukada Y (2004) Negative control of circadian clock regulator E4BP4 by casein kinase iepsilon-mediated phosphorylation. Curr Biol 14: 975–980.
[42]
Henriques R, Jang IC, Chua NH (2009) Regulated proteolysis in light-related signaling pathways. Curr Opin Plant Biol 12: 49–56.
[43]
Masri S, Zocchi L, Katada S, Mora E, Sassone-Corsi P (2012) The circadian clock transcriptional complex: metabolic feedback intersects with epigenetic control. Ann N Y Acad Sci 1264: 103–109.
[44]
Lin JD, Liu C, Li S (2008) Integration of energy metabolism and the mammalian clock. Cell Cycle 7: 453–457.
[45]
Fanjul-Moles ML, Escamilla-Chimal EG, Salceda R, Giulianini PG, Sanchez-Chavez G (2010) Circadian modulation of crustacean hyperglycemic hormone in crayfish eyestalk and retina. Chronobiol Int 27: 34–51.
[46]
Mandilaras K, Missirlis F (2012) Genes for iron metabolism influence circadian rhythms in Drosophila melanogaster. Metallomics 4: 928–936.
[47]
Edgar RS, Green EW, Zhao Y, van Ooijen G, Olmedo M, et al. (2012) Peroxiredoxins are conserved markers of circadian rhythms. Nature 485: 459–464.
[48]
Golombek DA, Rosenstein RE (2010) Physiology of circadian entrainment. Physiol Rev 90: 1063–1102.
[49]
Nieoullon A, Canolle B, Masmejean F, Guillet B, Pisano P, et al. (2006) The neuronal excitatory amino acid transporter EAAC1/EAAT3: Does it represent a major actor at the brain excitatory synapse? J Neurochem 98: 1007–1018.
[50]
Gotow T, Nishi T (2008) Simple photoreceptors in some invertebrates: physiological properties of a new photosensory modality. Brain Res 1225: 3–16.
[51]
Golombek DA, Agostino PV, Plano SA, Ferreyra GA (2004) Signaling in the mammalian circadian clock: The NO/cGMP pathway. Neurochem Int 45: 929–936.
[52]
Im SH, Taghert PH (2010) PDF receptor expression reveals direct interactions between circadian oscillators in drosophila. J Comp Neurol 518: 1925–1945.
[53]
Rajkumar P, Rollmann SM, Cook TA, Layne JE (2010) Molecular evidence for color discrimination in the atlantic sand fiddler crab, Uca pugilator. J Exp Biol 213: 4240–4248.
[54]
McGinnis JF, Whelan JP, Donoso LA (1992) Transient, cyclic changes in mouse visual cell gene products during the light-dark cycle. J Neurosci Res 31: 584–590.
[55]
Craft CM, Whitmore DH, Donoso LA (1990) Differential expression of mRNA and protein encoding retinal and pineal S-antigen during the light/dark cycle. J Neurochem 55: 1461–1473.
[56]
Sasagawa H, Narita R, Kitagawa Y, Kadowaki T (2003) The expression of genes encoding visual components is regulated by a circadian clock, light environment and age in the honeybee (Apis mellifera). Eur J Neurosci 17: 963–970.
[57]
Battelle BA, Williams CD, Schremser-Berlin JL, Cacciatore C (2000) Regulation of arrestin mRNA levels in Limulus lateral eye: Separate and combined influences of circadian efferent input and light. Vis Neurosci 17: 217–227.
[58]
Jeffery N (2012) The first genome size estimates for six species of krill (Malacostraca, Euphausiidae): Large genomes at the north and south poles. Polar Biol 35: 959–962.
[59]
Ward B, Heuermann K (2010) Transplex？ whole transcriptome amplification of RNA from low cell-number samples. Journal of Biomolecular Techniques 21: S30.
[60]
Zhulidov PA, Bogdanova EA, Shcheglov AS, Vagner LL, Khaspekov GL, et al. (2004) Simple cDNA normalization using kamchatka crab duplex-specific nuclease. Nucleic Acids Res 32: e37.
[61]
Zhulidov PA, Bogdanova EA, Shcheglov AS, Shagina IA, Vagner LL, et al. (2005) A method for the preparation of normalized cDNA libraries enriched with full-length sequences. Bioorg Khim 31: 186–194.
[62]
Bolstad BM, Irizarry RA, Astrand M, Speed TP (2003) A comparison of normalization methods for high density oligonucleotide array data based on variance and bias. Bioinformatics 19: 185–193.
Oster H, Damerow S, Hut RA, Eichele G (2006) Transcriptional profiling in the adrenal gland reveals circadian regulation of hormone biosynthesis genes and nucleosome assembly genes. J Biol Rhythms 21: 350–361.
[65]
Huang da W, Sherman BT, Lempicki RA (2009) Systematic and integrative analysis of large gene lists using DAVID bioinformatics resources. Nat Protoc 4: 44–57.
[66]
Thompson JD, Higgins DG, Gibson TJ (1994) CLUSTAL W: Improving the sensitivity of progressive multiple sequence alignment through sequence weighting, position-specific gap penalties and weight matrix choice. Nucleic Acids Res 22: 4673–4680.
[67]
Tamura K, Dudley J, Nei M, Kumar S (2007) MEGA4: Molecular evolutionary genetics analysis (MEGA) software version 4.0. Mol Biol Evol 24: 1596–1599.
[68]
Untergasser A, Cutcutache I, Koressaar T, Ye J, Faircloth BC, et al. (2012) Primer3 – new capabilities and interfaces. Nucleic Acids Res 40: e115.
[69]
Livak KJ, Schmittgen TD (2001) Analysis of relative gene expression data using real-time quantitative PCR and the 2(-delta delta C(T)) method. Methods 25: 402–408.
