Evolution of duplicated growth hormone genes in autotetraploid salmonid fishes

GHRH-SST-GH-IGF axis regulates crosstalk between growth and immunity in rainbow trout (Oncorhynchus mykiss) infected with Vibrio anguillarum

An energy trade-off is existed between immunological competence and growth. The axis of growth hormone releasing hormone, somatostatin, growth hormone, insulin-like growth factor (GHRH-SST-GH-IGF axis) regulates growth performances and immune competences in rainbow trout (Oncorhynchus mykiss). The salmonid-specific whole genome duplication event is known to result in duplicated copies of several key genes in GHRH-SST-GH-IGF axis. In this study, we evaluated the physiological functions of GHRH-SST-GH-IGF axis in regulating crosstalk between growth and immunity. Based on principal components analysis (PCA), we observed the overall expression profiles of GHRH-SST-GH-IGF axis were significantly altered by Vibrio anguillarum infection. Trout challenged with Vibrio anguillarum showed down-regulated igf1s subtypes and up-regulated igfbp1a1. The brain sst genes (sst1a, sst1b, sst3b and sst5) and igfpbs genes (igfbp4s and igfbp5b2) were significantly affected by V. anguillarum infection, while the igfbp4s, igfbp5s, igfbp6s and igf2bps genes showed significant changes in peripheral immune tissues in response to V. anguillarum infection. Gene enrichment analyses showed functional and signaling pathways associated with apoptosis (such as p53, HIF-1 or FoxO signaling) were activated. We further proposed a possible model that describes the IGF and IGFBPs-regulated interaction between cell growth and programmed death. Our study provided new insights into the physiological functions and potentially regulatory mechanisms of the GHRH-SST-GH-IGF axis, indicating the pleiotropic effects of GHRH-SST-GH-IGF axis in regulating crosstalk between growth and immunity in trout.

Effect of dietary threonine on growth performance and muscle growth, protein synthesis and antioxidant-related signalling pathways of hybrid catfish Pelteobagrus vachelli♀ × Leiocassis longirostris♂

The experiment was conducted to investigate the effects of dietary threonine (Thr) on growth performance and muscle growth, protein synthesis and antioxidant-related signalling pathways of hybrid catfish Pelteobagrus vachelli♀ × Leiocassis longirostris♂. A total of 1200 fish (14·19 (se 0·13) g) were randomly distributed into six groups with four replicates each, fed six diets with graded level of Thr (9·5, 11·5, 13·5, 15·4, 17·4 and 19·3 g/kg diets) for 56 d. Results showed (P < 0·05) that dietary Thr (1) increased percentage weight gain, specific growth rate, feed efficiency and protein efficiency ratio; (2) up-regulated growth hormone (GH), insulin-like growth factor 1 (IGF-1), proliferating cell nuclear antigen, myogenic regulation factors (MyoD, Myf5, MyoG and Mrf4) and myosin heavy chain (MyHC) mRNA levels; (3) increased muscle protein content via regulating the protein kinase B/target of rapamycin signalling pathway and (4) decreased malondialdehyde and protein carbonyl contents, increased catalase, glutathione-S-transferase, glutathione reductase and GSH activities, up-regulated mRNA levels of antioxidant enzymes related to NFE2-related factor 2 and γ-glutamylcysteine ligase catalytic subunit. These results suggest that Thr has a potential role to improve muscle growth and protein synthesis, which might be due to the regulation of GH-IGF system, muscle growth-related gene, antioxidative capacity and protein synthesis-related signalling pathways. Based on the quadratic regression analysis of specific growth rate, the Thr requirement of hybrid catfish (14·19-25·77 g) was estimated to be 13·77 g/kg of the diet (33·40 g/kg of dietary protein).

Characterization and expression analyses of somatolactin-α and -β genes in rare minnows (Gobiocypris rarus) following waterborne cadmium exposure

Using reverse transcription-polymerase chain reaction (RT-PCR) and RACE (rapid amplification of cDNA ends), somatolactin-α (rmSLα) and -β (rmSLβ) were identified from the pituitary gland of rare minnows (Gobiocypris rarus). The full-length cDNAs of these two genes were 1288 and 801 bp, encoding prepeptides of 250 and 228 amino acids residues, respectively. rmSLβ can be detected in the brain (including the pituitary), ovary, testis, and gill, while rmSLα was mainly expressed in the brain. On the other hand, rmSLα was expressed in all the fetal developmental stages; however, rmSLβ can just be detected in the stages since from 14 h post-fertilization (hpf). After exposure to acute waterborne cadmium (Cd), rmSLα was distinctly upregulated in juvenile rare minnows at all detected time points, from 24 to 96 h and 10 days, while rmSLβ was significantly altered only in 96 h or 10-day treatment groups. As for adults, acute Cd exposure caused alterations of both rmSLα and rmSLβ in the brain (containing the pituitary) at the 24 h; subchronic waterborne Cd treatment led to upregulation of rmSLα, while decrease of mSLβ in the brain. Alteration of rmSL transcripts following waterborne Cd exposure further confirmed the endocrine disruption of this heavy metal. Besides, exposure to as low as 5 μg/L Cd caused alteration of rmSLα, which suggested that rmSLα might be a potential biomarker for risk assessment of aquatic Cd.

Growth hormone transgenesis in coho salmon disrupts muscle immune function impacting cross-talk with growth systems

Suppression of growth during infection may aid resource allocation towards effective immune function. Past work supporting this hypothesis in salmonid fish revealed an immune-responsive regulation of the insulin-like growth factor (IGF) system, an endocrine pathway downstream of growth hormone (GH). Skeletal muscle is the main target for growth and energetic storage in fish, yet little is known about how its growth is regulated during an immune response. We addressed this knowledge gap by characterizing muscle immune responses in size-matched coho salmon (Oncorhynchus kisutch) achieving different growth rates. We compared a wild-type strain with two GH transgenic groups from the same genetic background achieving either maximal or suppressed growth, a design separating GH's direct effects from its influence on growth rate and nutritional state. Fish were sampled 30h post-injection with PBS (control) or mimics of bacterial or viral infection. We quantified mRNA expression levels for genes from the GH, GH receptor, IGF hormone, IGF1 receptor and IGF-binding protein families, along with immune genes involved in inflammatory or antiviral responses and muscle growth status marker genes. We demonstrate dampened immune function in GH transgenics compared to wild-type. The muscle of GH transgenics achieving rapid growth showed no detectable antiviral response, coupled with evidence of a constitutive inflammatory state. GH and IGF system gene expression was strongly altered by GH transgenesis and fast growth, both for baseline expression and responses to immune stimulation. Thus, GH transgenesis strongly disrupts muscle immune status and normal GH and IGF system expression responses to immune stimulation.

Linkage maps of the Atlantic salmon (Salmo salar) genome derived from RAD sequencing

BACKGROUND

Genetic linkage maps are useful tools for mapping quantitative trait loci (QTL) influencing variation in traits of interest in a population. Genotyping-by-sequencing approaches such as Restriction-site Associated DNA sequencing (RAD-Seq) now enable the rapid discovery and genotyping of genome-wide SNP markers suitable for the development of dense SNP linkage maps, including in non-model organisms such as Atlantic salmon (Salmo salar). This paper describes the development and characterisation of a high density SNP linkage map based on SbfI RAD-Seq SNP markers from two Atlantic salmon reference families.

RESULTS

Approximately 6,000 SNPs were assigned to 29 linkage groups, utilising markers from known genomic locations as anchors. Linkage maps were then constructed for the four mapping parents separately. Overall map lengths were comparable between male and female parents, but the distribution of the SNPs showed sex-specific patterns with a greater degree of clustering of sire-segregating SNPs to single chromosome regions. The maps were integrated with the Atlantic salmon draft reference genome contigs, allowing the unique assignment of ~4,000 contigs to a linkage group. 112 genome contigs mapped to two or more linkage groups, highlighting regions of putative homeology within the salmon genome. A comparative genomics analysis with the stickleback reference genome identified putative genes closely linked to approximately half of the ordered SNPs and demonstrated blocks of orthology between the Atlantic salmon and stickleback genomes. A subset of 47 RAD-Seq SNPs were successfully validated using a high-throughput genotyping assay, with a correspondence of 97% between the two assays.

CONCLUSIONS

This Atlantic salmon RAD-Seq linkage map is a resource for salmonid genomics research as genotyping-by-sequencing becomes increasingly common. This is aided by the integration of the SbfI RAD-Seq SNPs with existing reference maps and the draft reference genome, as well as the identification of putative genes proximal to the SNPs. Differences in the distribution of recombination events between the sexes is evident, and regions of homeology have been identified which are reflective of the recent salmonid whole genome duplication.

The role of hybridization in plant speciation

The importance of hybridization in plant speciation and evolution has been debated for decades, with opposing views of hybridization as either a creative evolutionary force or evolutionary noise. Hybrid speciation may occur at either the homoploid (i.e., between two species of the same ploidy) or the polyploid level, each with its attendant genetic and evolutionary consequences. Whereas allopolyploidy (i.e., resulting from hybridization and genome doubling) has long been recognized as an important mode of plant speciation, the implications of genome duplication have typically not been taken into account in most fields of plant biology. Recent developments in genomics are revolutionizing our views of angiosperm genomes, demonstrating that perhaps all angiosperms have likely undergone at least one round of polyploidization and that hybridization has been an important force in generating angiosperm species diversity. Hybridization and polyploid formation continue to generate species diversity, with several new allopolyploids having originated just within the past century or so. The origins of polyploid species-whether via hybridization between species or between genetically differentiated populations of a single species-and the immediate genetic consequences of polyploid formation are therefore receiving enthusiastic attention. The time is therefore right for a review of the role of hybridization in plant speciation.

A salmonid EST genomic study: genes, duplications, phylogeny and microarrays

Background

Salmonids are of interest because of their relatively recent genome duplication, and their extensive use in wild fisheries and aquaculture. A comprehensive gene list and a comparison of genes in some of the different species provide valuable genomic information for one of the most widely studied groups of fish.

Results

298,304 expressed sequence tags (ESTs) from Atlantic salmon (69% of the total), 11,664 chinook, 10,813 sockeye, 10,051 brook trout, 10,975 grayling, 8,630 lake whitefish, and 3,624 northern pike ESTs were obtained in this study and have been deposited into the public databases. Contigs were built and putative full-length Atlantic salmon clones have been identified. A database containing ESTs, assemblies, consensus sequences, open reading frames, gene predictions and putative annotation is available. The overall similarity between Atlantic salmon ESTs and those of rainbow trout, chinook, sockeye, brook trout, grayling, lake whitefish, northern pike and rainbow smelt is 93.4, 94.2, 94.6, 94.4, 92.5, 91.7, 89.6, and 86.2% respectively. An analysis of 78 transcript sets show Salmo as a sister group to Oncorhynchus and Salvelinus within Salmoninae, and Thymallinae as a sister group to Salmoninae and Coregoninae within Salmonidae. Extensive gene duplication is consistent with a genome duplication in the common ancestor of salmonids. Using all of the available EST data, a new expanded salmonid cDNA microarray of 32,000 features was created. Cross-species hybridizations to this cDNA microarray indicate that this resource will be useful for studies of all 68 salmonid species.

Conclusion

An extensive collection and analysis of salmonid RNA putative transcripts indicate that Pacific salmon, Atlantic salmon and charr are 94–96% similar while the more distant whitefish, grayling, pike and smelt are 93, 92, 89 and 86% similar to salmon. The salmonid transcriptome reveals a complex history of gene duplication that is consistent with an ancestral salmonid genome duplication hypothesis. Genome resources, including a new 32 K microarray, provide valuable new tools to study salmonids.

Isolation, characterization and comparison of Atlantic and Chinook salmon growth hormone 1 and 2

Background

Growth hormone (GH) is an important regulator of skeletal growth, as well as other adapted processes in salmonids. The GH gene (gh) in salmonids is represented by duplicated, non-allelic isoforms designated as gh1 and gh2. We have isolated and characterized gh-containing bacterial artificial chromosomes (BACs) of both Atlantic and Chinook salmon (Salmo salar and Oncorhynchus tshawytscha) in order to further elucidate our understanding of the conservation and regulation of these loci.

Results

BACs containing gh1 and gh2 from both Atlantic and Chinook salmon were assembled, annotated, and compared to each other in their coding, intronic, regulatory, and flanking regions. These BACs also contain the genes for skeletal muscle sodium channel oriented in the same direction. The sequences of the genes for interferon alpha-1, myosin alkali light chain and microtubule associated protein Tau were also identified, and found in opposite orientations relative to gh1 and gh2. Viability of each of these genes was examined by PCR. We show that transposon insertions have occurred differently in the promoters of gh, within and between each species. Other differences within the promoters and intronic and 3'-flanking regions of the four gh genes provide evidence that they have distinct regulatory modes and possibly act to function differently and/or during different times of salmonid development.

Conclusion

A core proximal promoter for transcription of both gh1 and gh2 is conserved between the two species of salmon. Nevertheless, transposon integration and regulatory element differences do exist between the promoters of gh1 and gh2. Additionally, organization of transposon families into the BACs containing gh1 and for the BACs containing gh2, are very similar within orthologous regions, but much less clear conservation is apparent in comparisons between the gh1- and gh2-containing paralogous BACs for the two fish species. This is consistent with the hypothesis that a burst of transposition activity occurred during the speciation events which led to Atlantic and Pacific salmon. The Chinook and other Oncorhynchus GH1s are strikingly different in comparison to the other GHs and this change is not apparent in the surrounding non-coding sequences.

Rainbow smelt (Osmerus mordax) genomic library and EST resources

Genomic resources in rainbow smelt (Osmerus mordax) enable us to examine the genome duplication process in salmonids and test hypotheses relating to the fate of duplicated genes. They further enable us to pursue physiological and ecological studies in smelt. A bacterial artificial chromosome library containing 52,410 clones with an average insert size of 146 kb was constructed. This library represents an 11-fold average coverage of the rainbow smelt (O. mordax) genome. In addition, several complementary deoxyribonucleic acid libraries were constructed, and 36,758 sequences were obtained and combined into 12,159 transcripts. Over half of these transcripts have been identified, several of which have been associated with cold adaptation. These basic resources show high levels of similarity (86%) to salmonid genes and provide initial support for genome duplication in the salmonid ancestor. They also facilitate identification of genes important to fish and direct us toward new technologies for other studies in fish biology.

Major quantitative trait loci affect resistance to infectious pancreatic necrosis in Atlantic salmon (Salmo salar)

Infectious pancreatic necrosis (IPN) is a viral disease currently presenting a major problem in the production of Atlantic salmon (Salmon salar). IPN can cause significant mortality to salmon fry within freshwater hatcheries and to smolts following transfer to seawater, although challenged populations show clear genetic variation in resistance. To determine whether this genetic variation includes loci of major effect, a genomewide quantitative trait loci (QTL) scan was performed within 10 full-sib families that had received a natural seawater IPN challenge. To utilize the large difference between Atlantic salmon male and female recombination rates, a two-stage mapping strategy was employed. Initially, a sire-based QTL analysis was used to detect linkage groups with significant effects on IPN resistance, using two to three microsatellite markers per linkage group. A dam-based analysis with additional markers was then used to confirm and position any detected QTL. Two genomewide significant QTL and one suggestive QTL were detected in the genome scan. The most significant QTL was mapped to linkage group 21 and was significant at the genomewide level in both the sire and the dam-based analyses. The identified QTL can be applied in marker-assisted selection programs to improve the resistance of salmon to IPN and reduce disease-related mortality.

Phylogenetic analysis of the insulin-like growth factor binding protein (IGFBP) and IGFBP-related protein gene families

Insulin-like growth factor (IGF) activity is regulated by six high affinity binding proteins (IGFBPs) and possibly by some of the nine IGFBP-related proteins (IGFBP-rPs). To determine the phylogenetic relationship of this proposed gene superfamily, we conducted maximum likelihood (ML) and Bayesian inference analyses on a matrix of amino acid sequences from a diversity of vertebrate species. A single most likely phylogram, ML bootstrap, and Bayesian consensus tree of 10,000,000 generations revealed a monophyletic IGFBP lineage independent of the IGFBP-rPs. The IGFBPs segregated into three distinct clades: IGFBP-1, -3, and -6. Subsequent gene duplication events within the IGFBP-1 and -3 clades resulted in the production and divergence of IGFBP-2 and -4 within the IGFBP-1 clade and IGFBP-5 in the IGFBP-3 clade. By contrast, the IGFBP-rPs were distributed paraphyletically into two clades: IGFBP-rP1, 5, and 6 in one clade and the CCN family (IGFBP-rP2-4,7-9) in another. A recently identified IGFBP-3 homolog in rainbow trout localized to the IGFBP-2 subclade. Subsequence analysis identified a RGD motif common to IGFBP-2 orthologs, but did not identify the nuclear localization sequence present in IGFBP-3 and -5 homologs. The putative trout IGFBP-3 was 36-55% identical to different IGFBP-2 proteins, but only 24-27% identical to IGFBP-3 proteins. These results suggest that the IGFBPs and IGFBP-rPs are at best distantly related and that the limited similarities likely resulted from exon shuffling. They also suggest that rainbow trout, and possibly other salmonids, possess two IGFBP-2 paralogs as the putative trout IGFBP-3 is misannotated.

Bursts and horizontal evolution of DNA transposons in the speciation of pseudotetraploid salmonids

Background

Several genome duplications have occurred in the evolutionary history of teleost fish. In returning to a stable diploid state, the polyploid genome reorganized, and large portions are lost, while the fish lines evolved to numerous species. Large scale transposon movement has been postulated to play an important role in the genome reorganization process. We analyzed the DNA sequence of several large loci in Salmo salar and other species for the presence of DNA transposon families.

Results

We have identified bursts of activity of 14 families of DNA transposons (12 Tc1-like and 2 piggyBac-like families, including 11 novel ones) in genome sequences of Salmo salar. Several of these families have similar sequences in a number of closely and distantly related fish, lamprey, and frog species as well as in the parasite Schistosoma japonicum. Analysis of sequence similarities between copies within the families of these bursts demonstrates several waves of transposition activities coinciding with salmonid species divergence. Tc1-like families show a master gene-like copying process, illustrated by extensive but short burst of copying activity, while the piggyBac-like families show a more random copying pattern. Recent families may include copies with an open reading frame for an active transposase enzyme.

Conclusion

We have identified defined bursts of transposon activity that make use of master-slave and random mechanisms. The bursts occur well after hypothesized polyploidy events and coincide with speciation events. Parasite-mediated lateral transfer of transposons are implicated.

The role of aquaporin 3 in teleost fish

The aquaporin isoform, AQP3 has now been identified in a number of different teleost fish species, with additional DNA sequence information on AQP3 genes in further fish species available in genome databases. In zebrafish (Danio rerio), the AQP3 gene is present as two duplicate isoforms resulting from a teleostean fish genome-wide duplication. A further splicoform/isoform has also been identified in rainbow trout (Oncorhynchus mykiss). The identification of these AQP3 isoforms in other fish species is consequently explored. The role of AQP3 in physiological/osmoregulatory processes, in various teleost organs is then described. In teleost gill, AQP3 is expressed in 'chloride' cells, and in some species, in other epithelial cell types, where it may have a number of different functions including the prevention of dehydration. In eel esophagus, immunohistochemistry shows that AQP3 is expressed in surface epithelial cells in the anterior esophagus, but in mucus cells within the epithelium of the posterior esophagus. In eel intestine, AQP3 is found in macrophage-like cells and probably plays no part in osmoregulatory processes. In the rectum, as in the posterior esophagus AQP3 is expressed in mucus cells. In eel kidney, AQP3 is expressed in a subset of renal tubules, and localizes to the apical pole of tubule cells. There is no apparent change in the location or protein abundance of renal AQP3 following the acclimation of eels from freshwater to seawater.

Characterisation and expression of the paired box protein 7 (Pax7) gene in polymorphic Arctic charr (Salvelinus alpinus)

Arctic charr (Salvelinus alpinus L.) from Lake Thingvallavatn, Iceland occur as four distinct morphs: large benthivorous (LB), dwarf benthivorous (DB), piscivorous (PI) and planktonivorous (PL). The morphs differ with respect to body size, head morphology, growth rate, and life history. The aim of this study was to investigate the paired box protein 7 (Pax7) gene as a candidate for such polymorphisms due to its importance in cranio-facial, skeletal muscle, and central nervous system development. No variation in coding and intronic sequences was found between morphs. We identified 10 alternate Pax7 isoforms with insertions/deletions: a four-residue (GNRT) deletion, a GEASS insertion truncated by the first serine residue (GEAS), and a thirteen-residue insertion (GQYA/TGPEYVYCGT). The latter insertion with a threonine (T) contains a putative casein kinase II (CK-2) phosphorylation site. Pax7 spatial expression patterns were identical in embryos of DB-, LB-, and PL-morphs, and were similar to those described for zebrafish Pax7c, but a difference in temporal expression for segmentation was observed between DB and LB morphs. At the end of segmentation, novel expression was observed in the mandibular region as two bilateral domains. The potential role of multiple alternative splicing of the Pax7 gene for the generation of different Arctic charr morphs is briefly discussed.

Varying signals of the effects of natural selection during teleost growth hormone gene evolution

The growth hormone (GH) gene of teleost fish exhibits a higher degree of variability compared with other vertebrate groups. However, the different selective constraints at the sequence level are not well understood. In this study, maximum-likelihood (ML) models of codon substitutions were used to investigate Darwinian adaptive evolution of the GH gene in teleost fishes. Complete GH gene sequences of 54 fish species were classified into 4 orders, and the variable nature of GH was examined by determining the dNand dSrate variation and the rates of molecular evolution for each teleost order. The results indicate that although the overall evolution rate for teleost GH is high ((1.15 ± 0.01) × 10–9substitutions/(aa site·y)) compared with the "slow phases" in mammals ((0.21 to 0.28 ± 0.05) × 10–9), the vital structure of this gene has been retained. While the majority of the amino acid changes appear to be due to relaxation of purifying selection, some positively selected sites were detected in regions with no specifically identified role in protein function. The positively selected regions observed in salmoniformes lineage suggests a possible role for positive selection driving functional divergence in paralogous forms of the GH gene after whole-genome duplication in this lineage.Key words: teleost fish, growth hormone, positive selection, synonymous substitution, non-synonymous substitution, molecular evolution.

Dietary rapeseed oil affects the expression of genes involved in hepatic lipid metabolism in Atlantic salmon (Salmo salar L.)

Supplies of marine fish oils (FO) are limited, and sustainable production in aquaculture dictates that alternatives that do not compromise fish health and product quality, such as vegetable oils, must be found. Nutrigenomics will increase our understanding of how nutrition influences metabolic pathways and homeostatic control, and may be used to measure and validate subtle changes in organ-specific, metabolic gene expression signatures. We compared 2 groups of Atlantic salmon fed diets containing 100% FO or 75% rapeseed oil (RO) for 42 wk. A small-scale cDNA microarray was constructed to screen for changes in the expression of lipid metabolism genes in the liver resulting from this partial substitution of RO for FO. Delta5 fatty acid desaturase gene expression was significantly greater in fish fed 75% RO than in fish fed the control diet; this was confirmed by quantitative real time PCR analysis. In addition, several genes, among these mitochondrial proteins, peroxisome proliferator-activated receptor gamma, as well as other transcription factors, coactivators, and signal transducers, showed significant differential regulation. This partially validated microarray may be used for further gene expression profiling using other dietary comparisons, and for further characterization of selected genes.