Transcriptional responses of the black-chinned tilapia Sarotherodon melanotheron to salinity extremes

Coping with salinity extremes: Gill transcriptome profiling in the black-chinned tilapia (Sarotherodon melanotheron)

Steeper and sometimes extreme salinity gradients increasingly affect aquatic organisms because of climate change. Hypersalinity habitats demand powerful physiological adaptive strategies. Few teleost species have the capacity to spend their whole life cycle in salinities way over seawater levels. Focusing on the multifunctional gill, we unraveled the tilapia S. melanotheron key strategies to cope with different environmental conditions, ranging from freshwater up to hypersaline habitats. De novo transcriptome assembly based on RNAseq allowed for the analysis of 40,967 annotated transcripts among samples collected in three wild populations at 0, 40 and 80 ‰. A trend analysis of the expression patterns revealed responses across the salinity gradient with different gene pathways involved. Genes linked to ion transport, pH regulation and cell surface receptor signaling were mainly upregulated in the high salinity habitat. We identified tight junction proteins that were critical in high salinity habitats and that were different from the well-known tightening junctional proteins identified and expressed in fresh water. Expression profiles also suggest a change in the vascular tone that could be linked to an osmorespiratory compromise not only in fresh water, but also in high salinity environments. A striking downregulation of genes linked to the immune system and to the heat shock response was observed suggesting an energetic trade-off between immunity and acclimation/adaptation in the hypersaline habitat. The high expression of transcripts coding for immune and heat shock response in the freshwater habitat suggests the establishment of powerful mechanisms to protect gills from environmental threats and to maintain protein integrity. Non-directional expression trends were also detected with an upregulation of genes only in the hypersaline habitat (80 ‰) or only in the marine habitat (40 ‰). Unravel physiological strategies in S. melanotheron populations will help to better understand the molecular basis of fish euryhalinity in salinity-contrasted environments.

Multidisciplinary long-term survey of Manila clam grown in farming sites subjected to different environmental conditions

In recent years recurrent bivalve mass mortalities considerably increased around the world, causing the collapse of natural and farmed populations. Venice Lagoon has historically represented one of the major production areas of the Manila clam Ruditapes philippinarum in Europe. However, in the last 20 years a 75 % decrease in the annual production has been experienced. While climate change and anthropogenic interventions may have played a key role in natural and farmed stocks reductions, no studies investigated at multiple levels the environmental stressors affecting farmed Manila clam to date. In this work we carried out a long-term monitoring campaign on Manila clam reared in four farming sites located at different distances from the southern Venice Lagoon inlet, integrating (meta)genomic approaches (i.e. RNA-seq; microbiota characterization), biometric measurements and chemical-physical parameters. Our study allowed to characterize the molecular mechanisms adopted by this species to cope with the different environmental conditions characterizing farming sites and to propose hypotheses to explain mortality events observed in recent years. Among the most important findings, the disruption of clam's immune response, the spread of Vibrio spp., and the up-regulation of molecular pathways involved in xenobiotic metabolism suggested major environmental stressors affecting clams farmed in sites placed close to Chioggia's inlet, where highest mortality was also observed. Overall, our study provides knowledge-based tools for managing Manila clam farming on-growing areas. In addition, the collected data is a snapshot of the time immediately before the commissioning of MoSE, a system of mobile barriers aimed at protecting Venice from high tides, and will represent a baseline for future studies on the effects of MoSE on clams farming and more in general on the ecology of the Venice Lagoon.

Effects of seawater and freshwater challenges on the Gh/Igf system in the saline-tolerant blackchin tilapia (Sarotherodon melanotheron)

Prolactin (Prl) and growth hormone (Gh) as well as insulin-like growth factor 1 (Igf1) are involved in the physiological adaptation of fish to varying salinities. The Igfs have been also ascribed other physiological roles during development, growth, reproduction and immune regulation. However, the main emphasis in the investigation of osmoregulatory responses has been the endocrine, liver-derived Igf1 route and local regulation within the liver and osmoregulatory organs. Few studies have focused on the impact of salinity alterations on the Gh/Igf-system within the neuroendocrine and immune systems and particularly in a salinity-tolerant species, such as the blackchin tilapia Sarotherodon melanotheron. This species is tolerant to hypersalinity and saline variations, but it is confronted by severe climate changes in the Saloum inverse estuary. Here we investigated bidirectional effects of increased salinity followed by its decrease on the gene regulation of prl, gh, igf1, igf2, Gh receptor and the tumor-necrosis factor a. A mixed population of sexually mature 14-month old blackchin tilapia adapted to freshwater were first exposed to seawater for one week and then to fresh water for another week. Brain, pituitary, head kidney and spleen were excised at 4 h, 1, 2, 3 and 7 days after both exposures and revealed differential expression patterns. This investigation should give us a better understanding of the role of the Gh/Igf system within the neuroendocrine and immune organs and the impact of bidirectional saline challenges on fish osmoregulation in non-osmoregulatory organs, notably the complex orchestration of growth factors and cytokines.

Physiological and metabolic responses to hypersalinity reveal interpopulation tolerance in the green macroalga Ulva compressa with different pollution histories

There is scarce investigation addressing interpopulation tolerance responses to address the influence of a history of chronic stress exposure, as that occurring in polluted environments, in photoautotrophs. We evaluated ecophysiological (photosynthetic activity) and metabolic (oxidative stress and damage) responses of two populations of green macroalga Ulva compressa from polluted (Ventanas) and non-polluted (Cachagua) localions of central Chile, and exposed to controlled hypersalinity conditions of 32 (control), 42, 62 and 82 psu (practical salinity units) for 6 h, 48 h and 6 d. Both primary production (ETR_max) and photosynthetic efficiency (α_ETR) were generally higher in the population from Cachagua compared to Ventanas at all times and salinities. Moreover, at most experimental times and salinities the population from Ventanas had greater levels of H₂O₂ and lipid peroxidation that individuals from Cachagua. Total ascorbate was higher in the population of Cachagua than Ventanas at 42 and 82 psu after 6 and 48 h, respectively, while at 6 d concentrations were similar between both populations at all salinities. Total glutathione was greater in both populations after 6 h at all salinities, but at 48 h its concentrations were higher only in the population from Cachagua, a trend that was maintained at 6 d under 82 psu only. Reduced and oxidized ascorbate (ASC and DHA, respectively) and glutathione (GSH and GSSG, respectively) demonstrated similar patterns between U. compressa populations, with an increase oxidation with greater salinities but efficient recycling to maintain sufficient batch of ASC and GSH. When assessing the expression of antioxidant enzymes catalase (CAT), superoxide dismutase (SOD) and dehydroascorbate reductase (DHAR), while the population of Ventanas displayed a general trend of upregulation with increasing salinities along the experiments, U. compressa from Cachagua revealed patterns of downregulation. Results demonstrated that although both populations were still viable after the applied hypersalinities during all experimental times, biological performance was usually more affected in the population from the Ventanas than Cachagua, likely due to a depressed baseline metabolism after a long history of exposition to environmental pollution.

Comparative analysis of the liver transcriptome in the red-eared slider Trachemys scripta elegans under chronic salinity stress

The red-eared slider (Trachemys scripta elegans), identified as one of the 100 most invasive species in the world, is a freshwater turtle originally from the eastern United States and northeastern Mexico. Field investigations have shown that T. s. elegans can survive and lay eggs in saline habitats. In order to understand the molecular mechanisms of salinity adaptation, high-throughput RNA-Seq was utilized to identify the changes in gene expression profiles in the liver of T. s. elegans in response to elevated salinity. We exposed individuals to 0, 5, or 15 psu (practical salinity units) for 30 days. A total of 157.21 million reads were obtained and assembled into 205138 unigenes with an average length of 620 bp and N50 of 964 bp. Of these, 1019 DEGs (differentially expressed genes) were found in the comparison of 0 vs. 5 psu, 1194 DEGs in 0 vs. 15 psu and 1180 DEGs in 5 vs. 15 psu, which are mainly related to macromolecule metabolic process, ion transport, oxidoreductase activity and generation of precursor metabolites and energy by GO (Gene Ontology) enrichment analyses. T. s. elegans can adapt itself into salinity by balancing the entry of sodium and chloride ions via the up-regulation expression genes of ion transport (potassium voltage-gated channel subfamily H member 5, KCNH5; erine/threonine-protein kinase 32, STK32; salt-inducible kinase 1, SIK1; adiponectin, ACDC), and by accumulating plasma urea and free amino acid via the up-regulation expression genes of amino acid metabolism (ornithine decarboxylase antizyme 3, OAZ3; glutamine synthetase, GLUL; asparaginase-like protein 1b, ASRGL; L-amino-acid oxidase-like, LAAO; sodium-dependent neutral amino acid transporter B, SLC6A15s; amino acid permease, SLC7A9) in response to osmotic regulation. An investment of energy to maintain their homeostatic balance is required to salinity adaptation, therefore, the genes related to energy production and conversion (F-ATPase protein 6, ATP6; cytochrome c oxidase subunit I, COX1; cytochrome c oxidase subunit III, COX3; cytochrome b, CYTb; cytochrome P450 17A1, CYP17A1) were up-regulated with the increase of gene expression associated with lipid metabolism (apolipoprotein E precursor, APoE; coenzyme Q-binding protein, CoQ10; high-density lipoprotein particle, SAA) and carbohydrate metabolism (HK, MIP). These findings improve our understanding of the underlying molecular mechanisms involved in salinity adaptation and provide general guidance to illuminate the invasion potential of T. s. elegans into saline environments.

Evidence of the Complexity of Gene Expression Analysis in Fish Wild Populations

The present work examines the induction of the band 3 anion transport protein, mitogen-activated protein kinase, and lactate dehydrogenase, respectively related to osmolyte transport, cell volume regulation, and energy production in the gills of two tilapia strains exposed to either freshwater or hypersaline water. Overall, genes showed similar expression patterns between strains. However, a wild population survey across a range of natural habitats and salinities did not reveal the expected patterns. Although significant, the correlations between gene expression and salinity were slightly ambiguous and did not show any link with phenotypic differences in life history traits previously reported between the same populations. The differential expression was also not associated with the population genetic structure inferred from neutral markers. The results suggest that the differential expression observed is not the result of evolutionary forces such as genetic drift or adaptation by natural selection. Instead, it can be speculated that genes responded to various abiotic and biotic stressors, including factors intrinsic to animals. This study provides clear evidence of the complexity of gene expression analysis in wild populations and shows that more attention needs to be paid when selecting candidates as potential biomarkers for monitoring adaptive responses to a specific environmental perturbation.

Evaluation of potential candidate genes involved in salinity tolerance in striped catfish (Pangasianodon hypophthalmus) using an RNA-Seq approach

Increasing salinity levels in freshwater and coastal environments caused by sea level rise linked to climate change is now recognized to be a major factor that can impact fish growth negatively, especially for freshwater teleost species. Striped catfish (Pangasianodon hypophthalmus) is an important freshwater teleost that is now widely farmed across the Mekong River Delta in Vietnam. Understanding the basis for tolerance and adaptation to raised environmental salinity conditions can assist the regional culture industry to mitigate predicted impacts of climate change across this region. Attempt of next generation sequencing using the ion proton platform results in more than 174 million raw reads from three tissue libraries (gill, kidney and intestine). Reads were filtered and de novo assembled using a variety of assemblers and then clustered together to generate a combined reference transcriptome. Downstream analysis resulted in a final reference transcriptome that contained 60,585 transcripts with an N50 of 683 bp. This resource was further annotated using a variety of bioinformatics databases, followed by differential gene expression analysis that resulted in 3062 transcripts that were differentially expressed in catfish samples raised under two experimental conditions (0 and 15 ppt). A number of transcripts with a potential role in salinity tolerance were then classified into six different functional gene categories based on their gene ontology assignments. These included; energy metabolism, ion transportation, detoxification, signal transduction, structural organization and detoxification. Finally, we combined the data on functional salinity tolerance genes into a hypothetical schematic model that attempted to describe potential relationships and interactions among target genes to explain the molecular pathways that control adaptive salinity responses in P. hypophthalmus. Our results indicate that P. hypophthalmus exhibit predictable plastic regulatory responses to elevated salinity by means of characteristic gene expression patterns, providing numerous candidate genes for future investigations.

Transcriptomic Analysis of Trout Gill Ionocytes in Fresh Water and Sea Water Using Laser Capture Microdissection Combined with Microarray Analysis

Fish gills represent a complex organ composed of several cell types that perform multiple physiological functions. Among these cells, ionocytes are implicated in the maintenance of ion homeostasis. However, because the ionocyte represents only a small percent of whole gill tissue, its specific transcriptome can be overlooked among the numerous cell types included in the gill. The objective of this study is to better understand ionocyte functions by comparing the RNA expression of this cell type in freshwater and seawater acclimated rainbow trout. To realize this objective, ionocytes were captured from gill cryosections using laser capture microdissection after immunohistochemistry. Then, transcriptome analyses were performed on an Agilent trout oligonucleotide microarray. Gene expression analysis identified 108 unique annotated genes differentially expressed between freshwater and seawater ionocytes, with a fold change higher than 3. Most of these genes were up-regulated in freshwater cells. Interestingly, several genes implicated in ion transport, extracellular matrix and structural cellular proteins appeared up-regulated in freshwater ionocytes. Among them, several ion transporters, such as CIC2, SLC26A6, and NBC, were validated by qPCR and/or in situ hybridization. The latter technique allowed us to localize the transcripts of these ion transporters in only ionocytes and more particularly in the freshwater cells. Genes involved in metabolism and also several genes implicated in transcriptional regulation, cell signaling and the cell cycle were also enhanced in freshwater ionocytes. In conclusion, laser capture microdissection combined with microarray analysis allowed for the determination of the transcriptional signature of scarce cells in fish gills, such as ionocytes, and aided characterization of the transcriptome of these cells in freshwater and seawater acclimated trout.

Plasticity of gene expression according to salinity in the testis of broodstock and F1 black-chinned tilapia, Sarotherodon melanotheron heudelotii

The black-chinned tilapia Sarotherodon melanotheron heudelotii Rüppell 1852 (Teleostei, Cichlidae) displays remarkable acclimation capacities. When exposed to drastic changes of salinity, which can be the case in its natural habitat, it develops quick physiological responses and keeps reproducing. The present study focused on the physiological impact of salinity on male reproductive capacities, using gene expression as a proxy of acclimation process. Two series of experimental fish were investigated: the first one was composed of fish maintained in freshwater for several generations and newly acclimated to salinities of 35 and 70, whereas the second one consisted of the descendants of the latter born and were raised under their native salinity. Expression patterns of 43 candidate genes previously identified from the testes of wild males was investigated in the three salinities and two generations. Twenty of them showed significant expression differences between salinities, and their predicted function revealed that most of them are involved in the osmotic tolerance of sperm cells and/or in the maintenance of sperm motility. A high level of expression variation was evidenced, especially for fish maintained in freshwater. In spite of this, gene expression patterns allowed the differentiation between fish raised in freshwater and those maintained in hypersaline water in both generations. Altogether, the results presented here suggest that this high variability of expression is likely to ensure the reproductive success of this species under varying salinities.

A transcriptomic analysis of striped catfish (Pangasianodon hypophthalmus) in response to salinity adaptation: De novo assembly, gene annotation and marker discovery

The striped catfish (Pangasianodon hypophthalmus) culture industry in the Mekong Delta in Vietnam has developed rapidly over the past decade. The culture industry now however, faces some significant challenges, especially related to climate change impacts notably from predicted extensive saltwater intrusion into many low topographical coastal provinces across the Mekong Delta. This problem highlights a need for development of culture stocks that can tolerate more saline culture environments as a response to expansion of saline water-intruded land. While a traditional artificial selection program can potentially address this need, understanding the genomic basis of salinity tolerance can assist development of more productive culture lines. The current study applied a transcriptomic approach using Ion PGM technology to generate expressed sequence tag (EST) resources from the intestine and swim bladder from striped catfish reared at a salinity level of 9ppt which showed best growth performance. Total sequence data generated was 467.8Mbp, consisting of 4,116,424 reads with an average length of 112bp. De novo assembly was employed that generated 51,188 contigs, and allowed identification of 16,116 putative genes based on the GenBank non-redundant database. GO annotation, KEGG pathway mapping, and functional annotation of the EST sequences recovered with a wide diversity of biological functions and processes. In addition, more than 11,600 simple sequence repeats were also detected. This is the first comprehensive analysis of a striped catfish transcriptome, and provides a valuable genomic resource for future selective breeding programs and functional or evolutionary studies of genes that influence salinity tolerance in this important culture species.

An integrated transcriptomic and comparative genomic analysis of differential gene expression in Arctic charr (Salvelinus alpinus) following seawater exposure

High-throughput RNA sequencing was employed to compare expression profiles in two Arctic charr (Salvelinus alpinus) families post seawater exposure to identify genes and biological processes involved in hypo-osmoregulation and regulation of salinity tolerance. To further understand the genetic architecture of hypo-osmoregulation, the genomic organization of differentially expressed (DE) genes was also analysed. Using a de novo gill transcriptome assembly we found over 2300 contigs to be DE. Major transporters from the seawater mitochondrion-rich cell (MRC) complex were up-regulated in seawater. Expression ratios for 257 differentially expressed contigs were highly correlated between families, suggesting they are strictly regulated. Based on expression profiles and known molecular pathways we inferred that seawater exposure induced changes in methylation states and elevated peroxynitrite formation in gill. We hypothesized that concomitance between DE immune genes and the transition to a hypo-osmoregulatory state could be related to Cl- sequestration by antimicrobial defence mechanisms. Gene Ontology analysis revealed that cell division genes were up-regulated, which could reflect the proliferation of ATP1α1b-type seawater MRCs. Comparative genomics analyses suggest that hypo-osmoregulation is influenced by the relative proximities among a contingent of genes on Arctic charr linkage groups AC-4 and AC-12 that exhibit homologous affinities with a region on stickleback chromosome Ga-I. This supports the hypothesis that relative gene location along a chromosome is a property of the genetic architecture of hypo-osmoregulation. Evidence of non-random structure between hypo-osmoregulation candidate genes was found on AC-1/11 and AC-28, suggesting that interchromosomal rearrangements played a role in the evolution of hypo-osmoregulation in Arctic charr.

Transcriptomics of salinity tolerance capacity in Arctic charr (Salvelinus alpinus): a comparison of gene expression profiles between divergent QTL genotypes

Osmoregulatory capabilities have played an important role in the evolution, dispersal, and diversification of vertebrates. To better understand the genetic architecture of hypo-osmoregulation in fishes and to determine which genes and biological processes affect intraspecific variation in salinity tolerance, we used mRNA sequence libraries from Arctic charr gill tissue to compare gene expression profiles in fish exhibiting divergent salinity tolerance quantitative trait locus (QTL) genotypes. We compared differentially expressed genes with QTL positions to gain insight about the nature of the underlying polymorphisms and examined gene expression within the context of genome organization to gain insight about the evolution of hypo-osmoregulation in fishes. mRNA sequencing of 18 gill tissue libraries produced 417 million reads, and the final reduced de novo transcriptome assembly consisted of 92,543 contigs. Families contained a similar number of differentially expressed contigs between high and low salinity tolerance capacity groups, and log2 expression ratios ranged from 10.4 to -8.6. We found that intraspecific variation in salinity tolerance capacity correlated with differential expression of immune response genes. Some differentially expressed genes formed clusters along linkage groups. Most clusters comprised gene pairs, though clusters of three, four, and eight genes were also observed. We postulated that conserved synteny of gene clusters on multiple ancestral and teleost chromosomes may have been preserved via purifying selection. Colocalization of QTL with differentially expressed genes suggests that polymorphisms in cis-regulatory elements are part of a majority of QTL.

Analysis of the black-chinned tilapia Sarotherodon melanotheron heudelotii reproducing under a wide range of salinities: from RNA-seq to candidate genes

The black-chinned tilapia Sarotherodon melanotheron heudelotii is an ecologically appealing model as it shows exceptional adaptive capacities, especially with regard to salinity. In spite of this, this species is devoid of genomic resources, which impedes the understanding of such remarkable features. De novo assembly of transcript sequences produced by next-generation sequencing technologies offers a rapid approach to obtain expressed gene sequences for non-model organisms. It also facilitates the development of quantitative real-time PCR (qPCR) assays for analysing gene expression under different environmental conditions. Nevertheless, obtaining accurate and reliable qPCR results from such data requires a number of validations prior to interpretation. The transcriptome of S. melanotheron was sequenced to discover transcripts potentially involved in the plasticity of male reproduction in response to salinity variations. A set of 54 candidate and reference genes was selected through a digital gene expression (DGE) approach, and a de novo qPCR assay using these genes was validated for further detailed expression analyses. A user-friendly web interface was created for easy handling of the sequence data. This sequence collection represents a major transcriptomic resource for S. melanotheron and will provide a useful tool for functional genomics and genetics studies.

A novel primer set for multilocus phylogenetic inference in East African cichlid fishes

The cichlid fishes in the East African Great Lakes are a prime model system for the study of adaptive radiation. Therefore, the availability of an elaborate phylogenetic framework is an important prerequisite. Previous phylogenetic hypotheses on East African cichlids are mainly based on mitochondrial and/or fragment-based markers, and, to date, no taxon-rich phylogeny exists that is based on multilocus DNA sequence data. Here, we present the design of an extensive new primer set (24 nuclear makers) for East African cichlids that will be used for multilocus phylogenetic analyses in the future. The primers are designed to work for both Sanger sequencing and next-generation sequencing with the 454 technology. As a proof of principle, we validate these primers in a phylogenetically representative set of 16 cichlid species from Lake Tanganyika and main river systems in the area and provide a basic evaluation of the markers with respect to marker length and diversity indices.

A high-resolution map of the Nile tilapia genome: a resource for studying cichlids and other percomorphs

Background

The Nile tilapia (Oreochromis niloticus) is the second most farmed fish species worldwide. It is also an important model for studies of fish physiology, particularly because of its broad tolerance to an array of environments. It is a good model to study evolutionary mechanisms in vertebrates, because of its close relationship to haplochromine cichlids, which have undergone rapid speciation in East Africa. The existing genomic resources for Nile tilapia include a genetic map, BAC end sequences and ESTs, but comparative genome analysis and maps of quantitative trait loci (QTL) are still limited.

Results

We have constructed a high-resolution radiation hybrid (RH) panel for the Nile tilapia and genotyped 1358 markers consisting of 850 genes, 82 markers corresponding to BAC end sequences, 154 microsatellites and 272 single nucleotide polymorphisms (SNPs). From these, 1296 markers could be associated in 81 RH groups, while 62 were not linked. The total size of the RH map is 34,084 cR₃₅₀₀ and 937,310 kb. It covers 88% of the entire genome with an estimated inter-marker distance of 742 Kb. Mapping of microsatellites enabled integration to the genetic map. We have merged LG8 and LG24 into a single linkage group, and confirmed that LG16-LG21 are also merged. The orientation and association of RH groups to each chromosome and LG was confirmed by chromosomal in situ hybridizations (FISH) of 55 BACs. Fifty RH groups were localized on the 22 chromosomes while 31 remained small orphan groups. Synteny relationships were determined between Nile tilapia, stickleback, medaka and pufferfish.

Conclusion

The RH map and associated FISH map provide a valuable gene-ordered resource for gene mapping and QTL studies. All genetic linkage groups with their corresponding RH groups now have a corresponding chromosome which can be identified in the karyotype. Placement of conserved segments indicated that multiple inter-chromosomal rearrangements have occurred between Nile tilapia and the other model fishes. These maps represent a valuable resource for organizing the forthcoming genome sequence of Nile tilapia, and provide a foundation for evolutionary studies of East African cichlid fishes.

Variation in gene expression along a salinity gradient in wild populations of the euryhaline black-chinned tilapia Sarotherodon melanotheron

This study evaluated variation in expression of 11 genes within and among six wild populations of the black-chinned tilapia Sarotherodon melanotheron distributed along a salinity gradient from 0 to 100. Previous laboratory studies had shown that expression of these genes was sensitive to water salinity; the current study confirmed that a number of them also varied in expression in wild populations along the salinity gradient. Principal component analysis (PCA) first distinguished two, not mutually exclusive, sets of genes: trade-off genes that were highly expressed at one or other extreme of the salinity gradient and stress genes that were up-regulated at the two salinity extremes (i.e. a U-shaped expression pattern). The PCA clearly partitioned the populations into three groups based on their gene expression patterns and their position along the salinity gradient: a freshwater (GL; 0) population, four brackish and seawater (GB, HB, SM, SF; ranging from 20 to 50) populations and a hypersaline (SK, 100) population. Individual variation in gene expression was significantly greater within the populations at the extreme compared to intermediate salinities. These results reveal phenotypically plastic regulation of gene expression in S. melanotheron, and greater osmoregulatory and plasticity costs at extreme salinities, where fitness-related traits are known to be altered.

The genetic basis of salinity tolerance traits in Arctic charr (Salvelinus alpinus)

BACKGROUND

The capacity to maintain internal ion homeostasis amidst changing conditions is particularly important for teleost fishes whose reproductive cycle is dependent upon movement from freshwater to seawater. Although the physiology of seawater osmoregulation in mitochondria-rich cells of fish gill epithelium is well understood, less is known about the underlying causes of inter- and intraspecific variation in salinity tolerance. We used a genome-scan approach in Arctic charr (Salvelinus alpinus) to map quantitative trait loci (QTL) correlated with variation in four salinity tolerance performance traits and six body size traits. Comparative genomics approaches allowed us to infer whether allelic variation at candidate gene loci (e.g., ATP1α1b, NKCC1, CFTR, and cldn10e) could have underlain observed variation.

RESULTS

Combined parental analyses yielded genome-wide significant QTL on linkage groups 8, 14 and 20 for salinity tolerance performance traits, and on 1, 19, 20 and 28 for body size traits. Several QTL exhibited chromosome-wide significance. Among the salinity tolerance performance QTL, trait co-localizations occurred on chromosomes 1, 4, 7, 18 and 20, while the greatest experimental variation was explained by QTL on chromosomes 20 (19.9%), 19 (14.2%), 4 (14.1%) and 12 (13.1%). Several QTL localized to linkage groups exhibiting homeologous affinities, and multiple QTL mapped to regions homologous with the positions of candidate gene loci in other teleosts. There was no gene × environment interaction among body size QTL and ambient salinity.

CONCLUSIONS

Variation in salinity tolerance capacity can be mapped to a subset of Arctic charr genomic regions that significantly influence performance in a seawater environment. The detection of QTL on linkage group 12 was consistent with the hypothesis that variation in salinity tolerance may be affected by allelic variation at the ATP1α1b locus. IGF2 may also affect salinity tolerance capacity as suggested by a genome-wide QTL on linkage group 19. The detection of salinity tolerance QTL in homeologous regions suggests that candidate loci duplicated from the salmonid-specific whole-genome duplication may have retained their function on both sets of homeologous chromosomes. Homologous affinities suggest that loci affecting salinity tolerance in Arctic charr may coincide with QTL for smoltification and salinity tolerance traits in rainbow trout. The effects of body size QTL appear to be independent of changes in ambient salinity.

An EST resource for tilapia based on 17 normalized libraries and assembly of 116,899 sequence tags

BACKGROUND

Large collections of expressed sequence tags (ESTs) are a fundamental resource for analysis of gene expression and annotation of genome sequences. We generated 116,899 ESTs from 17 normalized and two non-normalized cDNA libraries representing 16 tissues from tilapia, a cichlid fish widely used in aquaculture and biological research.

RESULTS

The ESTs were assembled into 20,190 contigs and 36,028 singletons for a total of 56,218 unique sequences and a total assembled length of 35,168,415 bp. Over the whole project, a unique sequence was discovered for every 2.079 sequence reads. 17,722 (31.5%) of these unique sequences had significant BLAST hits (e-value < 10(-10)) to the UniProt database.

CONCLUSION

Normalization of the cDNA pools with double-stranded nuclease allowed us to efficiently sequence a large collection of ESTs. These sequences are an important resource for studies of gene expression, comparative mapping and annotation of the forthcoming tilapia genome sequence.

Differential expression of the heat shock protein Hsp70 in natural populations of the tilapia, Sarotherodon melanotheron, acclimatised to a range of environmental salinities

Background

The relationship between environmental variation and induction of heat shock proteins (Hsps) has been much documented under experimental conditions. However, very little is known about such induction in natural populations acclimatised to prevailing environmental conditions. Furthermore, while induction of stress proteins has been well documented in response to environmental contaminants and thermal stressors, little is known about whether factors, such as extreme salinity, are also potential inductors. The black-chinned tilapia Sarotherodon melanotheron is unusual for its ability to colonise estuarine environments in West Africa that are characterised by extremely high salinities. The relationships between mRNA levels of the 70 kDa heat shock protein (Hsp70) and Na⁺, K⁺-ATPase1α (Naka) in the gills, environmental salinity, and a life-history trait (condition factor) were investigated in wild populations of this species sampled from three locations in the Saloum estuary, at salinities ranging from 40 to 100 psu.

Results

The highest Hsp70 and Naka mRNA levels, and the poorest condition factors were recorded in the most saline sampling site (100 psu). The Hsp70 and Naka mRNA were correlated amongst themselves and showed a direct positive correlation with environmental salinity, and a negative correlation with fish condition factor. Thus, the Hsp70 is constitutively overexpressed by S. melanotheron acclimatised to extreme hypersalinity.

Conclusions

These results indicate that, although S. melanotheron can colonise extremely saline environments, the overexpression of Hsp70 combined with the higher Naka mRNA expression reveals that this represents a chronic stress. The induction of Hsp70 was, therefore, a biomarker of chronic hyper-osmotic stress which presumably can be linked to the impaired growth performance and precocious reproduction that have been demonstrated in the populations at the extremely saline sites.