Transcriptomic and genomic evolution under constant cold in Antarctic notothenioid fish

Seasonal changes in hepatic gene expression reveal modulation of multiple processes in rainbow smelt (Osmerus mordax)

Rainbow smelt (Osmerus mordax) are freeze-resistant fish that accumulate glycerol and produce an antifreeze protein during winter. Quantitative reverse transcription PCR (qPCR) and subtractive hybridization studies have previously revealed five genes in rainbow smelt liver to be differentially regulated in winter in comparison with the fall when water temperatures are warmer. In order to further define the suite of processes that are regulated seasonally, we undertook a large-scale analysis of gene expression by hybridization of smelt cDNA to the salmonid 16K cGRASP microarray. In total, 69 genes were identified as up-regulated and 14 genes as down-regulated under winter conditions. A subset of these genes was examined for differential regulation by qPCR in the individual cDNA samples that were pooled for microarray analysis. Ten of the 15 genes tested showed significant change in the same direction as microarray results, whereas one showed significant change in the opposite direction. Fructose-bisphosphate aldolase B and the cytosolic NAD-dependent glycerol-3-phosphate dehydrogenase were among the most highly up-regulated genes, a result supporting a metabolic focus on glycerol synthesis during winter. Modulation of other processes, including endoplasmic reticulum stress, lipid metabolism and transport, and protein synthesis, was also suggested by the qPCR analysis of array-identified genes. The 15 genes were subsequently examined by qPCR for seasonal variation in expression over five sampling times between October and March, and ten showed significant variation in expression over the sampling period. Taken together, these results provide new understanding of the biochemical adaptations of vertebrates to an extremely low seasonal temperature.

Cold-adapted bacteria and the globin case study in the Antarctic bacterium Pseudoalteromonas haloplanktis TAC125

Environmental oxygen availability may play an important role in the evolution of polar marine organisms, as suggested by the physiological and biochemical strategies adopted by these organisms to acquire, deliver and scavenge oxygen. Stress conditions such as extreme temperatures increase the production of reactive oxygen species (ROS) in cells. Thus, in order to prevent cellular damage, adjustments in antioxidant defences are needed to maintain the steady-state concentration of ROS. Cold-adapted bacteria are generally acknowledged to achieve their physiological and ecological success in cold environments through structural and functional properties developed in their genomes. A short overview on the molecular adaptations of polar bacteria and in particular on the biological function of oxygen-binding proteins in Pseudoalteromonas haloplanktis TAC125, selected as a model, will be provided together with the role of oxygen and oxidative/nitrosative stress in regulating adaptive responses at cellular and molecular levels.

Genome enablement of the notothenioidei: genome size estimates from 11 species and BAC libraries from 2 representative taxa

The perciform suborder Notothenoidei provides a compelling opportunity to study the adaptive radiation of a marine species flock in the cold Southern Ocean surrounding Antarctica. To enable genome-level studies of these psychrophilic fishes, we estimated the sizes of the genomes of 11 Antarctic species and generated high-quality BAC libraries for 2, the notothen Notothenia coriiceps and the icefish Chaenocephalus aceratus. Our results indicate that evolution of phylogenetically derived notothenioid families, [e.g., the icefishes (Channichthyidae)], was accompanied by genome expansion. Species (n=6) of the basal family Nototheniidae had C values that ranged between 0.98 and 1.20 pg, whereas those of the icefishes, the notothenioid crown group, were 1.66-1.83 pg (n=4 species). The BAC libraries VMRC-19 (N. coriiceps) and VMRC-21 (C. aceratus) comprised 12X and 10X coverage of the respective genomes and had average insert sizes of 138 and 168 kb. Greater than 60% of paired BAC ends sampled from each library ( approximately 0.1% of each genome) contained repetitive sequences, and the repetitive element landscapes of the 2 genomes (13.4% of the N. coriiceps genome and 14.5% for C. aceratus) were similar. The representation and depth of coverage of the libraries were verified by identification of multiple Hox gene contigs: six discrete Hox clusters were found in N. coriiceps and at least five Hox clusters were found in C. aceratus. Given the unusual anatomical and physiological adaptations of the notothenioids, the availability of these BAC libraries sets the stage for expanded analysis of the psychrophilic mode of life.

Antarctic notothenioid fishes: genomic resources and strategies for analyzing an adaptive radiation

The perciform suborder Notothenoidei provides a compelling opportunity to study the adaptive radiation of a marine species-flock in the cold Southern Ocean that surrounds Antarctica. To facilitate genome-level studies of the diversification of these fishes, we present estimates of the genome sizes of 11 Antarctic species and describe the production of high-quality bacterial artificial chromosome (BAC) libraries for two, the red-blooded notothen Notothenia coriiceps and the white-blooded icefish Chaenocephalus aceratus. Our results indicate that evolution of phylogenetically derived notothenioid families (e.g., the crown group Channichthyidae [icefishes]), was accompanied by genome expansion. Six species from the basal family Nototheniidae had C-values between 0.98 and 1.20 pg, a range that is consistent with the genome sizes of proposed outgroups (e.g., percids) of the notothenioid suborder. In contrast, four icefishes had C-values in the range 1.66-1.83 pg. The BAC libraries VMRC-19 (N. coriiceps) and VMRC-21 (C. aceratus) comprise 12× and 10× coverage of the respective genomes and have average insert sizes of 138 and 168 kb. Paired BAC-end reads representing ∼0.1% of each genome showed that the repetitive element landscapes of the two genomes (13.4% of the N. coriiceps genome and 14.5% for C. aceratus) were similar. The availability of these high-quality and well-characterized BAC libraries sets the stage for targeted genomic analyses of the unusual anatomical and physiological adaptations of the notothenioids, some of which mimic human diseases. Here we consider the evolution of secondary pelagicism by various taxa of the group and illustrate the utility of Antarctic icefishes as an evolutionary-mutant model of human osteopenia (low-mineral density of bones).

A gene family-based method for interspecies comparisons of sequencing-based transcriptomes and its use in environmental adaptation analysis

We describe a new method for sequencing-based cross-species transcriptome comparisons and define a new metric for evaluating gene expression across species using protein-coding families as units of comparison. Using this measure transcriptomes from different species were evaluated by mapping them to gene families and integrating the mapping results with expression data. Statistical tests were applied to the transcriptome evaluation results to identify differentially expressed families. A Perl program named Pro-Diff was compiled to implement this method. To evaluate the method and provide an example of its use, two liver EST transcriptomes from two closely related fish that live in different temperature zones were compared. One EST library was from a recent sequencing project of Dissosticus mawsoni, a fish that lives in cold Antarctic sea waters, while the other was newly sequenced data (available at: http://www.fishgenome.org/polarbank/) from Notothenia angustata, a species that lives in temperate near-shore water of southern New Zealand. Results from the comparison were consistent with results inferred from phenotype differences and also with our previously published Gene Ontology-based method. The Pro-Diff program and operation manual can be downloaded from: http://www.fishgenome.org/download/Prodiff.rar.

A critical assessment of cross-species detection of gene duplicates using comparative genomic hybridization

Background

Comparison of genomic DNA among closely related strains or species is a powerful approach for identifying variation in evolutionary processes. One potent source of genomic variation is gene duplication, which is prevalent among individuals and species. Array comparative genomic hybridization (aCGH) has been successfully utilized to detect this variation among lineages. Here, beyond the demonstration that gene duplicates among species can be quantified with aCGH, we consider the effect of sequence divergence on the ability to detect gene duplicates.

Results

Using the X chromosome genomic content difference between male D. melanogaster and female D. yakuba and D. simulans, we describe a decrease in the ability to accurately measure genomic content (copy number) for orthologs that are only 90% identical. We demonstrate that genome characteristics (e.g. chromatin environment and non-orthologous sequence similarity) can also affect the ability to accurately measure genomic content. We describe a normalization strategy and statistical criteria to be used for the identification of gene duplicates among any species group for which an array platform is available from a closely related species.

Conclusions

Array CGH can be used to effectively identify gene duplication and genome content; however, certain biases are present due to sequence divergence and other genome characteristics resulting from the divergence between lineages. Highly conserved gene duplicates will be more readily recovered by aCGH. Duplicates that have been retained for a selective advantage due to directional selection acting on many loci in one or both gene copies are likely to be under-represented. The results of this study should inform the interpretation of both previously published and future work that employs this powerful technique.

Did gene family expansions during the Eocene-Oligocene boundary climate cooling play a role in Pooideae adaptation to cool climates?

Adaptation to cool environments is a common feature in the core group of the grass subfamily Pooideae (Triticeae and Poeae). This suggest an ancient evolutionary origin of low temperature stress tolerance dating back prior to the initiation of taxonomic divergence of core Pooideae species. Viewing the Pooideae evolution in a palaeo-climatic perspective reveals that taxonomic divergence of the core Pooideae group initiated shortly after a global super-cooling period at the Eocene-Oligocene boundary (approximately 33.5-26 Ma). This global climate cooling altered distributions of plants and animals and must have imposed selection pressure for improved low temperature stress responses. Lineage-specific gene family expansions are known to be involved in adaptation to new environmental stresses. In Pooideae, two gene families involved in low temperature stress response, the C-repeat binding factor (CBF) and fructosyl transferase (FT) gene families, has undergone lineage-specific expansions. We investigated the timing of these gene family expansions by molecular dating and found that Pooideae-specific expansion events in CBF and FT gene families took place during Eocene-Oligocene super-cooling period. We hypothesize that the E-O super-cooling exerted selection pressure for improved low temperature stress response and frost tolerance in a core Pooideae ancestor, and that those individuals with multiple copies of CBF and FT genes were favoured.

Genomic approaches with natural fish populations

Natural populations v. inbred stocks provide a much richer resource for identifying the effects of nucleotide substitutions because natural populations have greater polymorphism. Additionally, natural populations offer an advantage over most common research organisms because they are subject to natural selection, and analyses of these adaptations can be used to identify biologically important changes. Among fishes, these analyses are enhanced by having a wide diversity of species (>28 000 species, more than any other group of vertebrates) living in a huge range of environments (from below freezing to > 46 degrees C, in fresh water to salinities >40 ppt.). Moreover, fishes exhibit many different life-history and reproductive strategies and have many different phenotypes and social structures. Although fishes provide numerous advantages over other vertebrate models, there is still a dearth of available genomic tools for fishes. Fishes make up approximately half of all known vertebrate species, yet <0.2% of fish species have significant genomic resources. Nonetheless, genomic approaches with fishes have provided some of the first measures of individual variation in gene expression and insights into environmental and ecological adaptations. Thus, genomic approaches with natural fish populations have the potential to revolutionize fundamental studies of diverse fish species that offer myriad ecological and evolutionary questions.

Transcription profiling of acute temperature stress in the Antarctic plunderfish Harpagifer antarcticus

Harpagifer antarcticus (the Antarctic plunderfish), a shallow-water benthic fish distributed around the Antarctic Peninsula, is a member of the notothenioid family, one of whose adaptations to the cold waters of Antarctica has been the loss of the classic heat shock response. In order to gain a more comprehensive understanding of the effects of temperature stress on H. antarcticus, we constructed a liver cDNA library and a 10,371 feature microarray. This was hybridized with material from a time course series of animals held at 6°C for 48h. The resulting expression profiles show that this fish displays the classical vertebrate acute inflammatory response. There was also a pronounced signal for increased energy requirements via up-regulation of genes involved in the β oxidation of fatty acids and also a strong signature of response to oxidative stress. Genes in the latter category did not include the "classic" antioxidants such as glutathione S-transferase, but genes involved in the production of reducing potential in the form of NADPH, peroxisome proliferation via peroxisomal acyl co-enzyme A oxidase 1 and genes known to be up-regulated by hypoxia-inducible factor 1 (HIF1). These identifications provide clear support for oxygen being the whole animal limiting factor at least in acute short-term temperature challenges. The classical heat shock proteins were not up-regulated during this trial, although numerous clones for each were present on the gene chip, confirming the lack of this response in this species. These data significantly increase our knowledge of the cellular stress response from animals in this unique environment.

Molecular adaptations in haemoglobins of notothenioid fishes

Since haemoglobins of all animal species have the same haem group, differences in their properties, including oxygen affinity, electrophoretic mobility and pH sensitivity, must result from the interaction of the prosthetic group with specific amino-acid residues in the primary structure. For this reason, fish globins have been the subject of extensive studies in recent years, not only for their structural characteristics, but also because they offer the possibility to investigate the evolutionary history of these ancient molecules in marine and freshwater species living in a great variety of environmental conditions. This review summarizes the current knowledge on the structure, function and phylogeny of haemoglobins of notothenioid fishes. On the basis of crystallographic analysis, the evolution of the Root effect is analysed. Adaptation of the oxygen transport system in notothenioids seems to be based on evolutionary changes, involving levels of biological organization higher than the structure of haemoglobin. These include changes in the rate of haemoglobin synthesis or in regulation by allosteric effectors, which affect the amount of oxygen transported in blood. These factors are thought to be more important for short-term response to environmental challenges than previously believed.

Physiological adaptations of stressed fish to polluted environments: role of heat shock proteins

Fish are subjected to a wide variety of environmental stressors. Stressors affect fish at all life stages and the stress-specific responses that occur at the biochemical and physiological levels affect the overall health and longevity of such animals. In this review, the organ-specific alterations in fish that inhabit polluted environments are addressed in detail. Fish,like other vertebrates, have evolved strategies to counteract stress-mediated effects. Among the key strategies that fish have developed is the induction of HSPs. The primary functions of HSPs are to promote the proper folding or refolding of proteins, to prevent potentially damaging interactions with proteins, and aiding in the disassembly of formations of protein aggregates. Stress, a state of unbalanced tissue oxidation, causes a general disturbance in the cellular antioxidant and redox balance and evokes HSP70 overexpression. Distinct families of HSPs have diverse physiological functions, and their induction, which is regulated at the transcriptional level, is mediated by the activation of heat shock factors. Interestingly, HSPs also interact with a wide variety of signaling molecules that modulate stress-mediated apoptotic effects. Hence, HSP induction is of major importance for maintenance of cell homeostasis. HSP-mediated adaptation processes are regarded as a fundamental protective mechanism that decreases cellular sensitivity to damaging events. Thus, the adaptive expression of HSPs is a protective response that helps combat stress-induced conformational damage to proteins. Additional research is needed to gain further information on the functional significance and role of individual HSPs and to enhance the understanding of the molecular mechanisms by which they act. In addition, field studies are needed to allow comprehensive evaluation of the potential use of HSPs as biomarkers for environmental monitoring. Furthermore, the expression of HSPs in fish fluctuates in response to seasonal variation. Because HSPs serves as a tool for assessing the stressed state of individuals and/or populations, the impact of seasonal influences on constitutive and inducible factors of these proteins should also be elucidated. Such research will lead to a fundamental improvement in the understanding of the functional role of HSPs in response to natural environmental changes and may allow correlation of the action of HSPs at the molecular level with the whole organismal stress response, which, so far, remains unexplained.

Comparison of fatty acid composition in nine organs of the sympatric Antarctic teleost fish species Notothenia coriiceps and Notothenia rossii (Perciformes: Nototheniidae)

Fatty acid (FA) composition of nine organs from two closely related Antarctic fish species, Notothenia coriiceps and Notothenia rossii, was determined through gas chromatography with flame ionization detection. A data set for each species was obtained using major FA profiles from specimens caught in the sea waters of Admiralty Bay during the summer season. The FA profiles for both species are overall similar, but organ peculiarities have been found, which could reflect metabolic specificities and feeding habits between species. With the exception of liver, the most abundant FA in organs was the n-3 polyunsaturated FA. The total n-6 polyunsaturated FAs were minor components in all evaluated organs. Palmitic acid was identified as the major saturated FA, whereas oleic acid was the most represented of the monounsaturated FA in almost all assessed organs of both species. The n-3/n-6 ratios of all organs were higher than 3.5. Differences in individual FA and FA metabolic profiles of some organs observed between N. coriiceps and N. rossii suggest specific requirements in the mobilization, transport, incorporation, and/or catabolism of lipids that were reinforced by differences on some FA ratios expressing the activity coefficient of enzymes implicated on the FA pathway flux.

Evolution of the Antarctic teleost immunoglobulin heavy chain gene

Notothenioid teleosts underwent major modifications of their genome to adapt to the cooling of the Antarctic environment. In order to identify specific features of the Antarctic teleost immunoglobulin, transcripts encoding the constant region of the IgM heavy chain from 13 Antarctic and non-Antarctic notothenioid species were sequenced. The primary mRNA splicing for the membrane form was found to be atypical in the majority of Antarctic species, because it led to exclusion of two entire constant exons, and to inclusion of 39-nucleotide exons encoding an unusually long Extracellular Membrane-Proximal Domain (EMPD). Genomic DNA analysis revealed that each 39-nucleotide exon fell within a long sequence that was the reverse complement of an upstream region. Deduced amino acid sequence analysis lead to the identification of cysteine encoding codons in the 39-nucleotide exons, but not in the respective sequence counterpart, suggesting that these residues might play an important role in the folding of the EMPD.

Genomic reaction norms: using integrative biology to understand molecular mechanisms of phenotypic plasticity

Phenotypic plasticity is the development of different phenotypes from a single genotype, depending on the environment. Such plasticity is a pervasive feature of life, is observed for various traits and is often argued to be the result of natural selection. A thorough study of phenotypic plasticity should thus include an ecological and an evolutionary perspective. Recent advances in large-scale gene expression technology make it possible to also study plasticity from a molecular perspective, and the addition of these data will help answer long-standing questions about this widespread phenomenon. In this review, we present examples of integrative studies that illustrate the molecular and cellular mechanisms underlying plastic traits, and show how new techniques will grow in importance in the study of these plastic molecular processes. These techniques include: (i) heterologous hybridization to DNA microarrays; (ii) next generation sequencing technologies applied to transcriptomics; (iii) techniques for studying the function of noncoding small RNAs; and (iv) proteomic tools. We also present recent studies on genetic model systems that uncover how environmental cues triggering different plastic responses are sensed and integrated by the organism. Finally, we describe recent work on changes in gene expression in response to an environmental cue that persist after the cue is removed. Such long-term responses are made possible by epigenetic molecular mechanisms, including DNA methylation. The results of these current studies help us outline future avenues for the study of plasticity.

Molecular cold-adaptation of protein function and gene regulation: The case for comparative genomic analyses in marine ciliated protozoa

Euplotes focardii is a marine ciliated protozoan discovered in the Ross Sea near Terra Nova Bay, Antarctica. This organism is strictly psychrophilic, survives and reproduces optimally at 4-5 °C, and has a genome rich in A/T base pairs. Like other ciliated protozoans, Euplotes spp. are characterized by nuclear dimorphism: 1) the germline micronucleus contains the entire genome as large chromosomes; and 2) the somatic macronucleus (∼50 megabases, or 5% of the micronuclear genome) contains small linear DNA nanochromosomes [1-12 kilobases], each of which constitutes a single genetic unit. These characteristics make E. focardii an ideal model for genome-level analysis to understand the evolutionary mechanisms that determine the adaptation of organisms to cold environments. Here we describe two examples that are controlled by phylogenetically appropriate comparison with mesophilic and psychrotolerant Euplotes species: 1) the genes and encoded proteins of the E. focardii tubulin superfamily, including α-, β-, and γ-tubulins; and 2) the genes of the heat-shock protein (Hsp) 70 family. The tubulins provide particular insight into protein-level structural changes that are likely to facilitate microtubule nucleation and polymerization in an energy poor environment. By contrast, the hsp70 genes of E. focardii and of its psychrotolerant relative E. nobilii reveal adaptive alterations in the regulation of gene expression in the cold. The unique characteristics of the E. focardii genome and the results that we present here argue strongly for a concerted effort to characterize the relatively low complexity macronuclear genome of this psychrophilic organism.