The parallel evolution of dwarfism in Arctic charr is accompanied by adaptive divergence in mTOR-pathway gene expression

Nutritional status affects Igf1 regulation of skeletal muscle myogenesis, myostatin, and myofibrillar protein degradation pathways in gopher rockfish (Sebastes carnatus)

Insulin-like growth factor-1 (Igf1) regulates skeletal muscle growth in fishes by increasing protein synthesis and promoting muscle hypertrophy. When fish experience periods of insufficient food intake, they undergo slower muscle growth or even muscle wasting, and those changes emerge in part from nutritional modulation of Igf1 signaling. Here, we examined how food deprivation (fasting) modulates Igf1 regulation of liver and skeletal muscle gene expression in gopher rockfish (Sebastes carnatus), a nearshore rockfish of importance for commercial and recreational fisheries in the northeastern Pacific Ocean, to understand how food limitation impacts Igf regulation of muscle growth pathways. Rockfish were either fed or fasted for 14 d, after which a subset of fish from each group was treated with recombinant Igf1 from sea bream (Sparus aurata). Fish that were fasted lost body mass and had lower body condition, reduced hepatosomatic index, and lower plasma Igf1 concentrations, as well as a decreased abundance of igf1 gene transcripts in the liver, increased hepatic mRNAs for Igf binding proteins igfbp1a, igfbp1b, and igfbp3a, and decreased mRNA abundances for igfbp2b and a putative Igf acid labile subunit (igfals) gene. In skeletal muscle, fasted fish showed a reduced abundance of intramuscular igf1 mRNAs but elevated gene transcripts encoding Igf1 receptors A (igf1ra) and B (igf1rb), which also showed downregulation by Igf1. Fasting increased skeletal muscle mRNAs for myogenin and myostatin1, as well as ubiquitin ligase F-box only protein 32 (fbxo32) and muscle RING-finger protein-1 (murf1) genes involved in muscle atrophy, while concurrently downregulating mRNAs for myoblast determination protein 2 (myod2), myostatin2, and myogenic factors 5 (myf5) and 6 (myf6 encoding Mrf4). Treatment with Igf1 downregulated muscle myostatin1 and fbxo32 under both feeding conditions, but showed feeding-dependent effects on murf1, myf5, and myf6/Mrf4 gene expression indicating that Igf1 effects on muscle growth and atrophy pathways is contingent on recent food consumption experience.

Effects of different acute stressors on the regulation of appetite genes in the carp (Cyprinus carpio L.) brain

Our understanding of the timing of stress responses and specific roles of different regulatory pathways that drive stress responses is incomplete. In particular, the regulation of appetite genes as a consequence of exposure to different stressors has not been studied in sufficient detail in fish. Therefore, a stress trial was conducted with koi carp, aiming at identifying typical effects of stress on regulation of appetite genes. The stressors tank manipulation, air exposure and feed rewarding were chosen. The responses to these stressors were evaluated 10, 30 and 60 min after the stressors were applied. Orexigenic and anorexigenic genes were investigated in four different brain regions (telencephalon, hypothalamus, optic tectum and rhombencephalon). The results show that, apart from the typical appetite regulation in the hypothalamus, the different brain regions also display pronounced responses of appetite genes to the different stressors. In addition, several genes in the serotonergic, dopaminergic and gaba-related pathways were investigated. These genes revealed that rearing in pairs of two and opening of the tank lid affected anorexigenic genes, such as cart and cck, which were not changed by air exposure or feed rewarding. Moreover, distress and eustress led to limited, but distinguishable gene expression pattern changes in the investigated brain regions.

Genome resequencing clarifies phylogeny and reveals patterns of selection in the toxicogenomics model Pimephales promelas

Background

The fathead minnow (Pimephales promelas) is a model species for toxicological research. A high-quality genome reference sequence is available, and genomic methods are increasingly used in toxicological studies of the species. However, phylogenetic relationships within the genus remain incompletely known and little population-genomic data are available for fathead minnow despite the potential effects of genetic background on toxicological responses. On the other hand, a wealth of extant samples is stored in museum collections that in principle allow fine-scale analysis of contemporary and historical genetic variation.

Methods

Here we use short-read shotgun resequencing to investigate sequence variation among and within Pimephales species. At the genus level, our objectives were to resolve phylogenetic relationships and identify genes with signatures of positive diversifying selection. At the species level, our objective was to evaluate the utility of archived-sample resequencing for detecting selective sweeps within fathead minnow, applied to a population introduced to the San Juan River of the southwestern United States sometime prior to 1950.

Results

We recovered well-supported but discordant phylogenetic topologies for nuclear and mitochondrial sequences that we hypothesize arose from mitochondrial transfer among species. The nuclear tree supported bluntnose minnow (P. notatus) as sister to fathead minnow, with the slim minnow (P. tenellus) and bullhead minnow (P. vigilax) more closely related to each other. Using multiple methods, we identified 11 genes that have diversified under positive selection within the genus. Within the San Juan River population, we identified selective-sweep regions overlapping several sets of related genes, including both genes that encode the giant sarcomere protein titin and the two genes encoding the MTORC1 complex, a key metabolic regulator. We also observed elevated polymorphism and reduced differentation among populations (F_ST) in genomic regions containing certain immune-gene clusters, similar to what has been reported in other taxa. Collectively, our data clarify evolutionary relationships and selective pressures within the genus and establish museum archives as a fruitful resource for characterizing genomic variation. We anticipate that large-scale resequencing will enable the detection of genetic variants associated with environmental toxicants such as heavy metals, high salinity, estrogens, and agrichemicals, which could be exploited as efficient biomarkers of exposure in natural populations.

Genome-Wide Investigation of the Multiple Origins Hypothesis for Deep-Spawning Kokanee Salmon (Oncorhynchus nerka) across its Pan-Pacific Distribution

Salmonids have emerged as important study systems for investigating molecular processes underlying parallel evolution given their tremendous life history variation. Kokanee, the resident form of anadromous sockeye salmon (Oncorhynchus nerka), have evolved multiple times across the species' pan-Pacific distribution, exhibiting multiple reproductive ecotypes including those that spawn in streams, on lake-shores, and at lake depths >50 m. The latter has only been detected in 5 locations in Japan and British Columbia, Canada. Here, we investigated the multiple origins hypothesis for deep-spawning kokanee, using 9721 single nucleotide polymorphisms distributed across the genome analyzed for the vast majority of known populations in Japan (Saiko Lake) and Canada (Anderson, Seton, East Barrière Lakes) relative to stream-spawning populations in both regions. We detected 397 outlier loci, none of which were robustly identified in paired-ecotype comparisons in Japan and Canada independently. Bayesian clustering and principal components analyses based on neutral loci revealed 6 distinct clusters, largely associated with geography or translocation history, rather than ecotype. Moreover, a high level of divergence between Canadian and Japanese populations, and between deep- and stream-spawning populations regionally, suggests the deep-spawning ecotype independently evolved on the 2 continents. On a finer level, Japanese kokanee populations exhibited low estimates of heterozygosity, significant levels of inbreeding, and reduced effective population sizes relative to Canadian populations, likely associated with transplantation history. Along with preliminary evidence for hybridization between deep- and stream-spawning ecotypes in Saiko Lake, these findings should be considered within the context of on-going kokanee fisheries management in Japan.

Comparative transcriptome profiles of large and small bodied large-scale loaches cultivated in paddy fields

Fish culture in paddy fields is a traditional aquaculture mode, which has a long history in East Asia. Large-scale loach (Paramisgurnus dabryanus) fast growth is suitable for paddy fields aquaculture in China. The objective of this study was to identify differential expression genes (DEGs) in the brain, liver and muscle tissues between large (LG, top 5% of maximum total length) and small (SG, top 5% of minimum total length) groups using RNA-seq. In total, 150 fish were collected each week and 450 fish were collected at twelfth week from three paddy fields for all the experimental. Histological observation found that the muscle fibre diameter of LG loaches was greater than that of SG loaches. Transcriptome results revealed that the high expression genes (HEGs) in LG loaches (fold change ≥ 2, p < 0.05) were mainly concentrated in metabolic pathways, such as "Thyroid hormone signalling pathway", "Citrate cycle (TCA cycle)", "Carbon metabolism", "Fatty acid metabolism", and "Cholesterol metabolism", and the HEGs in SG loaches were enriched in the pathways related to environmental information processing such as "Cell adhesion molecules (CAMs)", "ECM- receptor interaction" and "Rap1 signalling pathway"; cellular processes such as "Tight junction", "Focal adhesion", "Phagosome" and "Adherens junction". Furthermore, IGFs gene family may play an important role in loach growth for their different expression pattern between the two groups. These findings can enhance our understanding about the molecular mechanism of different growth and development levels of loaches in paddy fields.

Parallelism in gene expression between foothill and alpine ecotypes in Arabidopsis arenosa

Parallel adaptation results from the independent evolution of similar traits between closely related lineages and allows us to test to what extent evolution is repeatable. Similar gene expression changes are often detected but the identity of genes shaped by parallel selection and the causes of expression parallelism remain largely unknown. By comparing genomes and transcriptomes of four distinct foothill-alpine population pairs across four treatments, we addressed the genetic underpinnings, plasticity and functional consequences of gene expression parallelism in alpine adaptation. Seeds of eight populations of Arabidopsis arenosa were raised under four treatments that differed in temperature and irradiance, factors varying strongly with elevation. Parallelism in differential gene expression between the foothill and alpine ecotypes was quantified by RNA-seq in leaves of young plants. By manipulating temperature and irradiance, we also tested for parallelism in plasticity (i.e., gene-environment interaction, GEI). In spite of global non-parallel patterns transcriptome wide, we found significant parallelism in gene expression at the level of individual loci with an over-representation of genes involved in biotic stress response. In addition, we demonstrated significant parallelism in GEI, indicating a shared differential response of the originally foothill versus alpine populations to environmental variation across mountain regions. A fraction of genes showing expression parallelism also encompassed parallel outliers for genomic differentiation, with greater enrichment of such variants in cis-regulatory elements in some mountain regions. In summary, our results suggest frequent evolutionary repeatability in gene expression changes associated with the colonization of a challenging environment that combines constitutive expression differences and plastic interaction with the surrounding environment.

Effects of food deprivation on plasma insulin-like growth factor-1 (Igf1) and Igf binding protein (Igfbp) gene transcription in juvenile cabezon (Scorpaenichthys marmoratus)

The growth hormone (GH)/insulin-like growth factor (Igf) endocrine axis regulates somatic growth in the face of changing environmental conditions. In actinopterygian fishes, food availability is a key modulator of the somatotropic axis, with lower food intake generally depressing liver Igf1 release to diminish growth. Igf1 signaling, however, also involves several distinct IGF binding proteins (Igfbps), and the functional roles of many of these Igfbps in affecting growth during shifting food availability remain uncertain. Here, we tested how complete food deprivation (fasting) affected gene transcription for paralogs of all six types of Igfbps in the liver and fast-twitch skeletal muscle of cabezon (Scorpaenichthys marmoratus), a nearshore marine fish important for recreational fisheries in the eastern North Pacific Ocean. Juvenile cabezon were maintained as either fed (6% mass food⋅g fish wet mass^-1⋅d^-1) or fasted for 14 d. Fasted fish exhibited a lower body condition (K), a depressed mass-specific growth rate (SGR), and reduced plasma concentrations of Igf1. In the liver, fasting reduced the relative abundance of gene transcripts encoding Igfbps igfbp2a and igfbp2b, while significantly elevating mRNA levels for igfbp1a, igfbp1b, igfbp3b, and igfbp4. Fasting also reduced hepatic mRNA levels of GH receptor-1 (ghr1) - but not GH receptor-2 (ghr2) - supporting the idea that changes in liver sensitivity to GH may underlie the decline in plasma Igf1 during food deprivation. In skeletal muscle, fasting downregulated gene transcripts encoding igf1, igfbp2b, igfbp5b, and igfbp6b, while also upregulating mRNAs for igf2 and ghr2. These data demonstrate isoform-specific regulation of Igfbps in liver and skeletal muscle in cabezon experiencing food deprivation and reinforce the idea that the repertoire of duplicated Igfbp genes that evolved in actinopterygian fishes supports a diverse scope of endocrine and paracrine functions.

Extensive genetic differentiation between recently evolved sympatric Arctic charr morphs

The availability of diverse ecological niches can promote adaptation of trophic specializations and related traits, as has been repeatedly observed in evolutionary radiations of freshwater fish. The role of genetics, environment, and history in ecologically driven divergence and adaptation, can be studied on adaptive radiations or populations showing ecological polymorphism. Salmonids, especially the Salvelinus genus, are renowned for both phenotypic diversity and polymorphism. Arctic charr (Salvelinus alpinus) invaded Icelandic streams during the glacial retreat (about 10,000 years ago) and exhibits many instances of sympatric polymorphism. Particularly, well studied are the four morphs in Lake Þingvallavatn in Iceland. The small benthic (SB), large benthic (LB), planktivorous (PL), and piscivorous (PI) charr differ in many regards, including size, form, and life history traits. To investigate relatedness and genomic differentiation between morphs, we identified variable sites from RNA-sequencing data from three of those morphs and verified 22 variants in population samples. The data reveal genetic differences between the morphs, with the two benthic morphs being more similar and the PL-charr more genetically different. The markers with high differentiation map to all linkage groups, suggesting ancient and pervasive genetic separation of these three morphs. Furthermore, GO analyses suggest differences in collagen metabolism, odontogenesis, and sensory systems between PL-charr and the benthic morphs. Genotyping in population samples from all four morphs confirms the genetic separation and indicates that the PI-charr are less genetically distinct than the other three morphs. The genetic separation of the other three morphs indicates certain degree of reproductive isolation. The extent of gene flow between the morphs and the nature of reproductive barriers between them remain to be elucidated.

A way forward with eco evo devo: an extended theory of resource polymorphism with postglacial fishes as model systems

A major goal of evolutionary science is to understand how biological diversity is generated and altered. Despite considerable advances, we still have limited insight into how phenotypic variation arises and is sorted by natural selection. Here we argue that an integrated view, which merges ecology, evolution and developmental biology (eco evo devo) on an equal footing, is needed to understand the multifaceted role of the environment in simultaneously determining the development of the phenotype and the nature of the selective environment, and how organisms in turn affect the environment through eco evo and eco devo feedbacks. To illustrate the usefulness of an integrated eco evo devo perspective, we connect it with the theory of resource polymorphism (i.e. the phenotypic and genetic diversification that occurs in response to variation in available resources). In so doing, we highlight fishes from recently glaciated freshwater systems as exceptionally well-suited model systems for testing predictions of an eco evo devo framework in studies of diversification. Studies on these fishes show that intraspecific diversity can evolve rapidly, and that this process is jointly facilitated by (i) the availability of diverse environments promoting divergent natural selection; (ii) dynamic developmental processes sensitive to environmental and genetic signals; and (iii) eco evo and eco devo feedbacks influencing the selective and developmental environments of the phenotype. We highlight empirical examples and present a conceptual model for the generation of resource polymorphism - emphasizing eco evo devo, and identify current gaps in knowledge.

Influence of feed ration size on somatic and muscle growth in landlocked dwarf and farmed Atlantic salmon Salmo salar

We examined the possible adaptation of the dwarf Bleke population of Atlantic salmon Salmo salar from Lake Byglandsfjord in southern Norway to limited food resources. The growth performance and muscle development in juvenile Bleke and farmed S. salar under satiated or restricted (50%) feeding were examined for 10 months, starting 3 weeks after first-feeding stage. Four-thousand fish were divided into four replicated groups and random samples of 16-40 fish per group were measured six times during the experiment. The two strains showed no significant difference in mean body mass when fed restricted ration, but the individual variation was considerably higher in the farmed fish. Both Bleke and farmed S. salar grew significantly faster when fed to satiation, but the farmed S. salar showed much higher gain in mass and were three times heavier (201.5 g vs 66.7 g) and possessed twice as many fast muscle fibres (179,682 vs 84,779) compared with landlocked S. salar after 10 months. Farmed fish fed full ration displayed both hypertrophic and hyperplasic muscle growth, while the increased growth in Bleke S. salar was entirely associated with a larger fibre diameter. The landlocked Bleke strain has apparently adapted to low food availability by minimising the metabolic costs of maintenance and growth through reduced dominance hierarchies and by an increase in average muscle fibre diameter relative to the ancestral condition.

Differential gene expression during early development in recently evolved and sympatric Arctic charr morphs

Phenotypic differences between closely related taxa or populations can arise through genetic variation or be environmentally induced, leading to altered transcription of genes during development. Comparative developmental studies of closely related species or variable populations within species can help to elucidate the molecular mechanisms related to evolutionary divergence and speciation. Studies of Arctic charr (Salvelinus alpinus) and related salmonids have revealed considerable phenotypic variation among populations and in Arctic charr many cases of extensive variation within lakes (resource polymorphism) have been recorded. One example is the four Arctic charr morphs in the ∼10,000 year old Lake Thingvallavatn, which differ in numerous morphological and life history traits. We set out to investigate the molecular and developmental roots of this polymorphism by studying gene expression in embryos of three of the morphs reared in a common garden set-up. We performed RNA-sequencing, de-novo transcriptome assembly and compared gene expression among morphs during an important timeframe in early development, i.e., preceding the formation of key trophic structures. Expectedly, developmental time was the predominant explanatory variable. As the data were affected by some form of RNA-degradation even though all samples passed quality control testing, an estimate of 3'-bias was the second most common explanatory variable. Importantly, morph, both as an independent variable and as interaction with developmental time, affected the expression of numerous transcripts. Transcripts with morph effect, separated the three morphs at the expression level, with the two benthic morphs being more similar. However, Gene Ontology analyses did not reveal clear functional enrichment of transcripts between groups. Verification via qPCR confirmed differential expression of several genes between the morphs, including regulatory genes such as AT-Rich Interaction Domain 4A (arid4a) and translin (tsn). The data are consistent with a scenario where genetic divergence has contributed to differential expression of multiple genes and systems during early development of these sympatric Arctic charr morphs.

Insulin-Like Growth Factor-Binding Proteins of Teleost Fishes

The insulin-like growth factor (Igf) binding protein (Igfbp) family has a broad range of physiological functions and a fascinating evolutionary history. This review focuses on the Igfbps of teleost fishes, where genome duplication events have diversified gene repertoire, function, and physiological regulation-with six core Igfbps expanded into a family of over twenty genes in some lineages. In addition to briefly summarizing the current state of knowledge on teleost Igfbp evolution, function, and expression-level regulation, we highlight gaps in our understanding and promising areas for future work.

Chronic stress inhibits growth and induces proteolytic mechanisms through two different nonoverlapping pathways in the skeletal muscle of a teleost fish

Chronic stress detrimentally affects animal health and homeostasis, with somatic growth, and thus skeletal muscle, being particularly affected. A detailed understanding of the underlying endocrine and molecular mechanisms of how chronic stress affects skeletal muscle growth remains lacking. To address this issue, the present study assessed primary (plasma cortisol), secondary (key components of the GH/IGF system, muscular proteolytic pathways, and apoptosis), and tertiary (growth performance) stress responses in fine flounder ( Paralichthys adspersus) exposed to crowding chronic stress. Levels of plasma cortisol, glucocorticoid receptor 2 ( gr2), and its target genes ( klf15 and redd1) mRNA increased significantly only at 4 wk of crowding ( P < 0.05). The components of the GH/IGF system, including ligands, receptors, and their signaling pathways, were significantly downregulated at 7 wk of crowding ( P < 0.05). Interestingly, chronic stress upregulated the ubiquitin-proteasome pathway and the intrinsic apoptosis pathways at 4wk ( P < 0.01), whereas autophagy was only significantly activated at 7 wk ( P < 0.05), and meanwhile the ubiquitin-proteasome and the apoptosis pathways returned to control levels. Overall growth was inhibited in fish in the 7-wk chronic stress trial ( P < 0.05). In conclusion, chronic stress directly affects muscle growth and downregulates the GH/IGF system, an action through which muscular catabolic mechanisms are promoted by two different and nonoverlapping proteolytic pathways. These findings provide new information on molecular mechanisms involved in the negative effects that chronic stress has on muscle anabolic/catabolic signaling balance.

Upwelling-derived oceanographic conditions impact growth performance and growth-related gene expression in intertidal fish

Growth is one of the main biological processes in aquatic organisms that is affected by environmental fluctuations such as upwelling (characterized by food-rich waters). In fish, growth is directly related with skeletal muscle increase; which represents the largest tissue of body mass. However, the effects of upwelling on growth, at the physiological and molecular level, are unknown. This study used Girella laevifrons (one of the most abundant intertidal fish in Eastern South Pacific) as a biological model, considering animals from upwelling (U) and non-upwelling (NU) areas. Here, we evaluated the effect of nutritional composition and food availability on growth performance and expression of key growth-related genes (insulin-kike growth factor 1 (igf1) and myosin heavy-chain (myhc)) and atrophy-related genes (muscle ring-finger 1 (murf1), F-box only protein 32 (atrogin-1) and BCL2/adenovirus E1B 19kDa-interacting protein 3 (bnip3)). We reported that, among zones, U fish displayed higher growth performance in response to nutritional composition, specifically between protein- and fiber-rich diets (~1g). We also found in NU fish that atrophy-related genes were upregulated with fiber-rich diet and during fasting (~2-fold at minimum respect U). In conclusion, our results suggest that the growth potential of upwelling fish may be a consequence of differential muscle gene expression. Our data provide a preliminary approach contributing on how upwelling influence fish growth at the physiological and molecular levels. Future studies are required to gain further knowledge about molecular differences between U and NU animals, as well as the possible applications of this knowledge in the aquaculture industry.

Divergent regulation of insulin-like growth factor binding protein genes in cultured Atlantic salmon myotubes under different models of catabolism and anabolism

Much attention has been given to insulin-like growth factor (Igf) pathways that regulate the balance of skeletal muscle protein synthesis and breakdown in response to a range of extrinsic and intrinsic signals. However, we have a less complete understanding of how the same signals modulate muscle mass upstream of such signalling, through a family of functionally-diverse Igf-binding proteins (Igfbps) that modify the availability of Igfs to the cell receptor Igf1r. We exposed cultured myotubes from Atlantic salmon (Salmo salar L.) to treatments recapturing three catabolic signals: inflammation (interleukin-1β), stress (dexamethasone) and fasting (amino acid deprivation), plus one anabolic signal: recovery of muscle mass post-fasting (supplementation of fasted myotubes with Igf-I and amino acids). The intended phenotype of treatments was confirmed by significant changes in myotube diameter and immunofluorescent staining of structural proteins. We quantified the mRNA-level regulation of the full expressed Igf and Igfbp gene complement across a post-treatment time course, along with marker genes for muscle structural protein synthesis, as well as muscle breakdown, via the ubiquitin-proteasome and autophagy systems. Our results highlight complex, non-overlapping responses of Igfbp family members to the different treatments, suggesting that the profile of expressed Igfbps is differentially regulated by distinct signals promoting similar muscle remodelling phenotypes. We also demonstrate divergent regulation of salmonid-specific gene duplicates of igfbp5b1 and igfbp5b2 under distinct catabolic and anabolic conditions. Overall, this study increases our understanding of the regulation of Igfbp genes in response to signals that promote remodelling of skeletal muscle.

Targeted sequencing for high-resolution evolutionary analyses following genome duplication in salmonid fish: Proof of concept for key components of the insulin-like growth factor axis

High-throughput sequencing has revolutionised comparative and evolutionary genome biology. It has now become relatively commonplace to generate multiple genomes and/or transcriptomes to characterize the evolution of large taxonomic groups of interest. Nevertheless, such efforts may be unsuited to some research questions or remain beyond the scope of some research groups. Here we show that targeted high-throughput sequencing offers a viable alternative to study genome evolution across a vertebrate family of great scientific interest. Specifically, we exploited sequence capture and Illumina sequencing to characterize the evolution of key components from the insulin-like growth (IGF) signalling axis of salmonid fish at unprecedented phylogenetic resolution. The IGF axis represents a central governor of vertebrate growth and its core components were expanded by whole genome duplication in the salmonid ancestor ~95Ma. Using RNA baits synthesised to genes encoding the complete family of IGF binding proteins (IGFBP) and an IGF hormone (IGF2), we captured, sequenced and assembled orthologous and paralogous exons from species representing all ten salmonid genera. This approach generated 299 novel sequences, most as complete or near-complete protein-coding sequences. Phylogenetic analyses confirmed congruent evolutionary histories for all nineteen recognized salmonid IGFBP family members and identified novel salmonid-specific IGF2 paralogues. Moreover, we reconstructed the evolution of duplicated IGF axis paralogues across a replete salmonid phylogeny, revealing complex historic selection regimes - both ancestral to salmonids and lineage-restricted - that frequently involved asymmetric paralogue divergence under positive and/or relaxed purifying selection. Our findings add to an emerging literature highlighting diverse applications for targeted sequencing in comparative-evolutionary genomics. We also set out a viable approach to obtain large sets of nuclear genes for any member of the salmonid family, which should enable insights into the evolutionary role of whole genome duplication before additional nuclear genome sequences become available.

Transcriptional dynamics of immune, growth and stress related genes in skeletal muscle of the fine flounder (Paralichthys adpersus) during different nutritional statuses

The effects of stress on immune activity and growth in early vertebrates have not been studied in detail. The present study used fine flounder (Paralichthys adspersus) skeletal muscle as a model to evaluate molecules involved in the stress response, including the glucocorticoid receptors, foxo1/3, and the target genes of these. Additionally, immune markers (il-1β and tnfα) and effector molecules of atrophy (bnip3, caspase-3, and lc3) were assessed. These molecules were analyzed during periods of long-term fasting and refeeding. During fasting, gene expression related to the stress response and atrophy increased; whereas immune markers were down-regulated. During refeeding, atrophy- and stress-related gene expression significantly decreased. In contrast, immune markers were up-regulated. These results provide novel insight on the control of growth in the skeletal muscle of a non-mammalian species under a stressful condition, suggesting that growth, stress, and immune activity in muscle are closely related and coordinated by orchestrated transcriptional dynamics.

The developmental transcriptome of contrasting Arctic charr (Salvelinus alpinus) morphs

Species and populations with parallel evolution of specific traits can help illuminate how predictable adaptations and divergence are at the molecular and developmental level. Following the last glacial period, dwarfism and specialized bottom feeding morphology evolved rapidly in several landlocked Arctic charrSalvelinus alpinuspopulations in Iceland.  To study the genetic divergence between small benthic morphs and limnetic morphs, we conducted RNA-sequencing charr embryos at four stages in early development. We studied two stocks with contrasting morphologies: the small benthic (SB) charr from Lake Thingvallavatn and Holar aquaculture (AC) charr.The data reveal significant differences in expression of several biological pathways during charr development. There was also an expression difference between SB- and AC-charr in genes involved in energy metabolism and blood coagulation genes. We confirmed differing expression of five genes in whole embryos with qPCR, includinglysozymeandnatterin-likewhich was previously identified as a fish-toxin of a lectin family that may be a putative immunopeptide. We also verified differential expression of 7 genes in the developing head that associated consistently with benthic v.s.limnetic morphology (studied in 4 morphs). Comparison of single nucleotide polymorphism (SNP) frequencies reveals extensive genetic differentiation between the SB and AC-charr (~1300 with more than 50% frequency difference). Curiously, three derived alleles in the otherwise conserved 12s and 16s mitochondrial ribosomal RNA genes are found in benthic charr.The data implicate multiple genes and molecular pathways in divergence of small benthic charr and/or the response of aquaculture charr to domestication. Functional, genetic and population genetic studies on more freshwater and anadromous populations are needed to confirm the specific loci and mutations relating to specific ecological traits in Arctic charr.

The developmental transcriptome of contrasting Arctic charr ( Salvelinus alpinus) morphs

Species and populations with parallel evolution of specific traits can help illuminate how predictable adaptations and divergence are at the molecular and developmental level. Following the last glacial period, dwarfism and specialized bottom feeding morphology evolved rapidly in several landlocked Arctic charr Salvelinus alpinus populations in Iceland.  

To study the genetic divergence between small benthic morphs and limnetic morphs, we conducted RNA-sequencing charr embryos at four stages in early development. We studied two stocks with contrasting morphologies: the small benthic (SB) charr from Lake Thingvallavatn and Holar aquaculture (AC) charr.

The data reveal significant differences in expression of several biological pathways during charr development. There was also an expression difference between SB- and AC-charr in genes involved in energy metabolism and blood coagulation genes. We confirmed differing expression of five genes in whole embryos with qPCR, including lysozyme and natterin-like which was previously identified as a fish-toxin of a lectin family that may be a putative immunopeptide. We also verified differential expression of 7 genes in the developing head that associated consistently with benthic v.s.limnetic morphology (studied in 4 morphs). Comparison of single nucleotide polymorphism (SNP) frequencies reveals extensive genetic differentiation between the SB and AC-charr (~1300 with more than 50% frequency difference). Curiously, three derived alleles in the otherwise conserved 12s and 16s mitochondrial ribosomal RNA genes are found in benthic charr.

The data implicate multiple genes and molecular pathways in divergence of small benthic charr and/or the response of aquaculture charr to domestication. Functional, genetic and population genetic studies on more freshwater and anadromous populations are needed to confirm the specific loci and mutations relating to specific ecological traits in Arctic charr.

The TORC1/P70S6K and TORC1/4EBP1 signaling pathways have a stronger contribution on skeletal muscle growth than MAPK/ERK in an early vertebrate: Differential involvement of the IGF system and atrogenes

Knowledge about the underlying mechanisms, particularly the signaling pathways that account for muscle growth in vivo in early vertebrates is still scarce. Fish (Paralichthys adspersus) were fasted for 3weeks to induce a catabolic period of strong muscle atrophy. Subsequently, fish were refed for 2weeks to induce compensatory muscle hypertrophy. During refeeding, the fish were treated daily with either rapamycin (TORC blocker), PD98059 (MEK blocker), or PBS (V; vehicle), or were untreated (C; control). Rapamycin and PD98059 differentially impaired muscle cellularity in vivo, growth performance, and the expression of growth-related genes, and the inhibition of TORC1 had a greater impact on fish muscle growth than the inhibition of MAPK. Blocking TORC1 inhibited the phosphorylation of P70S6K and 4EBP1, two downstream components activated by TORC1, thus affecting protein contents in muscle. Concomitantly, the gene expression in muscle of igf-1, 2 and igfbp-4, 5 was down-regulated while the expression of atrogin-1, murf-1, and igfbp-2, 3 was up-regulated. Muscle hypertrophy was abolished and muscle atrophy was promoted, which finally affected body weight. TORC2 complex was not affected by rapamycin. On the other hand, the PD98059 treatment triggered ERK inactivation, a downstream component activated by MEK. mRNA contents of igf-1 in muscle were down-regulated, and muscle hypertrophy was partially impaired. The present study provides the first direct data on the in vivo contribution of TORC1/P70S6K, TORC1/4EBP1, and MAPK/ERK signaling pathways in the skeletal muscle of an earlier vertebrate, and highlights the transcendental role of TORC1 in growth from the cellular to organism level.

Transcriptional dynamics of a conserved gene expression network associated with craniofacial divergence in Arctic charr

Background

Understanding the molecular basis of craniofacial variation can provide insights into key developmental mechanisms of adaptive changes and their role in trophic divergence and speciation. Arctic charr (Salvelinus alpinus) is a polymorphic fish species, and, in Lake Thingvallavatn in Iceland, four sympatric morphs have evolved distinct craniofacial structures. We conducted a gene expression study on candidates from a conserved gene coexpression network, focusing on the development of craniofacial elements in embryos of two contrasting Arctic charr morphotypes (benthic and limnetic).

Results

Four Arctic charr morphs were studied: one limnetic and two benthic morphs from Lake Thingvallavatn and a limnetic reference aquaculture morph. The presence of morphological differences at developmental stages before the onset of feeding was verified by morphometric analysis. Following up on our previous findings that Mmp2 and Sparc were differentially expressed between morphotypes, we identified a network of genes with conserved coexpression across diverse vertebrate species. A comparative expression study of candidates from this network in developing heads of the four Arctic charr morphs verified the coexpression relationship of these genes and revealed distinct transcriptional dynamics strongly correlated with contrasting craniofacial morphologies (benthic versus limnetic). A literature review and Gene Ontology analysis indicated that a significant proportion of the network genes play a role in extracellular matrix organization and skeletogenesis, and motif enrichment analysis of conserved noncoding regions of network candidates predicted a handful of transcription factors, including Ap1 and Ets2, as potential regulators of the gene network. The expression of Ets2 itself was also found to associate with network gene expression. Genes linked to glucocorticoid signalling were also studied, as both Mmp2 and Sparc are responsive to this pathway. Among those, several transcriptional targets and upstream regulators showed differential expression between the contrasting morphotypes. Interestingly, although selected network genes showed overlapping expression patterns in situ and no morph differences, Timp2 expression patterns differed between morphs.

Conclusion

Our comparative study of transcriptional dynamics in divergent craniofacial morphologies of Arctic charr revealed a conserved network of coexpressed genes sharing functional roles in structural morphogenesis. We also implicate transcriptional regulators of the network as targets for future functional studies.

Electronic supplementary material

The online version of this article (doi:10.1186/2041-9139-5-40) contains supplementary material, which is available to authorized users.

Differential gene expression in small and large rainbow trout derived from two seasonal spawning groups

Background

Growth in fishes is regulated via many environmental and physiological factors and is shaped by the genetic background of each individual. Previous microarray studies of salmonid growth have examined fish experiencing either muscle wastage or accelerated growth patterns following refeeding, or the influence of growth hormone and transgenesis. This study determines the gene expression profiles of genetically unmanipulated large and small fish from a domesticated salmonid strain reared on a typical feeding regime. Gene expression profiles of white muscle and liver from rainbow trout (Oncorhynchus mykiss) from two seasonal spawning groups (September and December lots) within a single strain were examined when the fish were 15 months of age to assess the influence of season (late fall vs. onset of spring) and body size (large vs. small).

Results

Although IGFBP1 gene expression was up-regulated in the livers of small fish in both seasonal lots, few expression differences were detected in the liver overall. Faster growing Dec. fish showed a greater number of differences in white muscle expression compared to Sept. fish. Significant differences in the GO Generic Level 3 categories ‘response to external stimulus’, ‘establishment of localization’, and ‘response to stress’ were detected in white muscle tissue between large and small fish. Larger fish showed up-regulation of cytoskeletal component genes while many genes related to myofibril components of muscle tissue were up-regulated in small fish. Most of the genes up-regulated in large fish within the ‘response to stress’ category are involved in immunity while in small fish most of these gene functions are related to apoptosis.

Conclusions

A higher proportion of genes in white muscle compared to liver showed similar patterns of up- or down-regulation within the same size class across seasons supporting their utility as biomarkers for growth in rainbow trout. Differences between large and small Sept. fish in the ‘response to stress’ and ‘response to external stimulus’ categories for white muscle tissue, suggests that smaller fish have a greater inability to handle stress compared to the large fish. Sampling season had a significant impact on the expression of genes related to the growth process in rainbow trout.

Regulation of skeletal muscle growth in fish by the growth hormone--insulin-like growth factor system

The growth hormone (GH)-insulin-like growth factor (IGF) system is the key promoter of growth in vertebrates; however, how this system modulates muscle mass in fish is just recently becoming elucidated. In fish, the GH induces muscle growth by modulating the expression of several genes belonging to the myostatin (MSTN), atrophy, GH, and IGF systems as well as myogenic regulatory factors (MRFs). The GH controls the expression of igf1 via Janus kinase 2 (JAK2)/signal transducers and activators of the transcription 5 (STAT5) signaling pathway, but it seems that it is not the major regulator. These mild effects of the GH on igf1 expression in fish muscle seem to be related with the presence of higher contents of truncated GH receptor1 (tGHR1) than full length GHR (flGHR1). IGFs in fish stimulate myogenic cell proliferation, differentiation, and protein synthesis through the MAPK/ERK and PI3K/AKT/TOR signaling pathways, concomitant with abolishing protein degradation and atrophy via the PI3K/AKT/FOXO signaling pathway. Besides these signaling pathways control the expression of several genes belonging to the atrophy and IGF systems. Particularly, IGFs and amino acid control the expression of igf1, thus, suggesting other of alternative signaling pathways regulating the transcription of this growth factor. The possible role of IGF binding proteins (IGFBPs) and the contribution of muscle-derived versus hepatic-produced IGF1 on fish muscle growth is also addressed. Thus, a comprehensive overview on the GH-IGF system regulating fish skeletal muscle growth is presented, as well as perspectives for future research in this field.

Isolation and selection of suitable reference genes for real-time PCR analyses in the skeletal muscle of the fine flounder in response to nutritional status: assessment and normalization of gene expression of growth-related genes

In the present study, different reference genes were isolated, and their stability in the skeletal muscle of fine flounder subjected to different nutritional states was assessed using geNorm and NormFinder. The combinations between 18S and ActB; Fau and 18S; and Fau and Tubb were chosen as the most stable gene combinations in feeding, long-term fasting and refeeding, and short-term refeeding conditions, respectively. In all periods, ActB was identified as the single least stable gene. Subsequently, the expression of the myosin heavy chain (MYH) and the insulin-like growth factor-I receptor (IGF-IR) was assessed. A large variation in MYH and IGF-IR expression was found depending on the reference gene that was chosen for normalizing the expression of both genes. Using the most stable reference genes, mRNA levels of MYH decreased and IGF-IR increased during fasting, with both returning to basal levels during refeeding. However, the drop in mRNA levels for IGF-IR occurred during short-term refeeding, in contrast with the observed events in the expression of MYH, which occurred during long-term refeeding. The present study highlights the vast differences incurred when using unsuitable versus suitable reference genes for normalizing gene expression, pointing out that normalization without proper validation could result in a bias of gene expression.

Differentiation at the MHCIIα and Cath2 loci in sympatric Salvelinus alpinus resource morphs in Lake Thingvallavatn

Northern freshwater fish may be suitable for the genetic dissection of ecological traits because they invaded new habitats after the last ice age (∼10.000 years ago). Arctic charr (Salvelinus alpinus) colonizing streams and lakes in Iceland gave rise to multiple populations of small benthic morphotypes, often in sympatry with a pelagic morphotype. Earlier studies have revealed significant, but subtle, genetic differentiation between the three most common morphs in Lake Thingvallavatn. We conducted a population genetic screen on four immunological candidate genes Cathelicidin 2 (Cath2), Hepcidin (Hamp), Liver expressed antimicrobial peptide 2a (Leap-2a), and Major Histocompatibility Complex IIα (MHCIIα) and a mitochondrial marker (D-loop) among the three most common Lake Thingvallavatn charr morphs. Significant differences in allele frequencies were found between morphs at the Cath2 and MHCIIα loci. No such signal was detected in the D-loop nor in the other two immunological genes. In Cath2 the small benthic morph deviated from the other two (FST  = 0.13), one of the substitutions detected constituting an amino acid replacement polymorphism in the antimicrobial peptide. A more striking difference was found in the MHCIIα. Two haplotypes were very common in the lake, and their frequency differed greatly between the morphotypes (from 22% to 93.5%, FST  = 0.67). We then expanded our study by surveying the variation in Cath2 and MHCIIα in 9 Arctic charr populations from around Iceland. The populations varied greatly in terms of allele frequencies at Cath2, but the variation did not correlate with morphotype. At the MHCIIα locus, the variation was nearly identical to the variation in the two benthic morphs of Lake Thingvallavatn. The results are consistent with a scenario where parts of the immune systems have diverged substantially among Arctic charr populations in Iceland, after colonizing the island ∼10.000 years ago.

Validation of reference genes for expression studies during craniofacial development in arctic charr

Arctic charr (Salvelinus alpinus) is a highly polymorphic species and in Lake Thingvallavatn, Iceland, four phenotypic morphs have evolved. These differences in morphology, especially in craniofacial structures are already apparent during embryonic development, indicating that genes important in the formation of the craniofacial features are expressed differentially between the morphs. In order to generate tools to examine these expression differences in Arctic charr, the aim of the present study was to identify reference genes for quantitative real-time PCR (qPCR). The specific aim was to select reference genes which are able to detect very small expression differences among different morphs. We selected twelve candidate reference genes from the literature, identified corresponding charr sequences using data derived from transcriptome sequencing (RNA-seq) and examined their expression using qPCR. Many of the candidate reference genes were found to be stably expressed, yet their quality-rank as reference genes varied considerably depending on the type of analysis used. In addition to commonly used software for reference gene validation, we used classical statistics to evaluate expression profiles avoiding a bias for reference genes with similar expression patterns (co-regulation). Based on these analyses we chose three reference genes, ACTB, UB2L3 and IF5A1 for further evaluation. Their consistency was assessed in an expression study of three known craniofacially expressed genes, sparc (or osteonectin), matrix metalloprotease 2 (mmp2) and sox9 (sex-determining region Y box 9 protein) using qPCR in embryo heads derived from four charr groups at three developmental time points. The three reference genes were found to be very suitable for studying expression differences between the morphotypes, enabling robust detection of small relative expression changes during charr development. Further, the results showed that sparc and mmp2 are differentially expressed in embryos of different Arctic charr morphotypes.

Evolution of ancient functions in the vertebrate insulin-like growth factor system uncovered by study of duplicated salmonid fish genomes

Whole-genome duplication (WGD) was experienced twice by the vertebrate ancestor (2 rounds; 2R), again by the teleost fish ancestor (3R) and most recently in certain teleost lineages (4R). Consequently, vertebrate gene families are often expanded in 3R and 4R genomes. Arguably, many types of “functional divergence” present across 2R gene families will exceed that between 3R/4R paralogs of genes comprising 2R families. Accordingly, 4R offers a form of replication of 2R. Examining whether this concept has implications for molecular evolutionary research, we studied insulin-like growth factor (IGF) binding proteins (IGFBPs), whose six 2R family members carry IGF hormones and regulate interactions between IGFs and IGF1-receptors (IGF1Rs). Using phylogenomic approaches, we resolved the complete IGFBP repertoire of 4R-derived salmonid fishes (19 genes; 13 more than human) and established evolutionary relationships/nomenclature with respect to WGDs. Traits central to IGFBP action were determined for all genes, including atomic interactions in IGFBP–IGF1/IGF2 complexes regulating IGF–IGF1R binding. Using statistical methods, we demonstrate that attributes of these protein interfaces are overwhelming a product of 2R IGFBP family membership, explain 49–68% of variation in IGFBP mRNA concentration in several different tissues, and strongly predict the strength and direction of IGFBP transcriptional regulation under differing nutritional states. The results support a model where vertebrate IGFBP family members evolved divergent structural attributes to provide distinct competition for IGFs with IGF1Rs, predisposing different functions in the regulation of IGF signaling. Evolution of gene expression then acted to ensure the appropriate physiological production of IGFBPs according to their structural specializations, leading to optimal IGF-signaling according to nutritional-status and the endocrine/local mode of action. This study demonstrates that relatively recent gene family expansion can facilitate inference of functional evolution within ancient genetic systems.

The evolutionary ecology of dwarfism in three-spined sticklebacks

1. Body size is a defining phenotypic trait, but the ecological causes of its evolution are poorly understood. Most studies have considered only a single putative causal agent and have failed to recognise that different environmental agents are often correlated. 2. Darwin suggested that although trait variation across populations is often associated with abiotic variation, evolution is more likely to be driven by biotic factors correlated with the abiotic variation. This hypothesis has received little explicit attention. 3. We use structural equation modelling to quantify the relative importance of abiotic (pH, metal concentrations) and biotic (competition, predation) factors in the evolution of body size in three-spined sticklebacks Gasterosteus aculeatus on the island of North Uist, Scotland. We combine phenotypic data from multiple isolated populations, detailed characterisation of their environment and a common garden experiment that establishes the genetic basis of size differences. 4. Three-spined sticklebacks on North Uist show almost unprecedented intraspecific evolution of body size that has taken place rapidly (<16,000 years). The smallest fish mature at only 7% of the mass of ancestral, anadromous fish. Dwarfism is associated with reduced abundance of a smaller competitor species, the nine-spined stickleback Pungitius pungitius, and with low pH indicative of poor resource conditions. Dwarfism also tends to occur where an important predator, the brown trout Salmo trutta, is also small. The abundance of P. pungitius and the size of S. trutta are themselves related to underlying abiotic environmental variation. 5. Despite the close association between abiotic and biotic factors across populations, our results support Darwin's hypothesis that biotic factors, associated with variation in the abiotic environment, are more important in explaining evolution than is abiotic variation per se. This study demonstrates the importance of considering the relationships between environmental variables before conclusions can be drawn about the causes of (body size) evolution on islands.

Experimental selection for body size at age modifies early life-history traits and muscle gene expression in adult zebrafish

The short generation time of the zebrafish (Danio rerio) was exploited to investigate the effects of selection for body size at age on early life-history traits and on the transcriptional response to a growth stimulus in skeletal muscle of adult fish. Replicate populations were either unselected (U-lineage) or subjected to four generations of experimental selection for small (S-lineage) or large (L-lineage) body size at 90 days post-fertilization. Body mass was on average 16.3% and 41.0% higher in the L- than in the U- and S-lineages, respectively. Egg diameter was 6.4% lower with 13% less yolk in the S-lineage compared with the other lineages. Maternal transcripts for igf2r, bmpr1aa, igf1ar, igf2a, igfbp5a, ghra and igfbp3 in 2-4 cell stage embryos were higher in the L- than in the S-lineage. Larvae from the L-lineage were significantly larger, but survivorship at the end of the first month was similar between lineages. Gene expression was measured in the fast muscle of adult fish fasted for 7 days and then re-fed to satiation for 48 h. The expression of 11 insulin-like growth factor pathway genes and 12 other nutritionally responsive genes was similar for the S- and L-lineages as was gut fullness with feeding. Transcript abundance for four genes (igf1a, igf2r, igfbp1a and igfbp1b) showed either regulated or constitutive differences between the S- and L-lineages. For example, igf2 receptor transcript abundance was higher and igbp1a/b transcript abundance was lower in the L- than in the S-lineage, consistent with an effect of selection on insulin-like growth factor signalling.

Dynamic transcriptional regulation of autocrine/paracrine igfbp1, 2, 3, 4, 5, and 6 in the skeletal muscle of the fine flounder during different nutritional statuses

The IGF-binding proteins (IGFBPs) play a dual role in the regulation of the activity and bioavailability of IGFs in different tissues. Diverse evidence has shown that IGFBPs can inhibit and/or potentiate IGF actions. In this study, igfbp1, 2, 3, 4, 5, and 6 were isolated in the fine flounder, a flat fish species that shows slow growth and inherent Gh resistance in muscle. Subsequently, the expression of all igfbps was assessed in the skeletal muscle of flounder that underwent different nutritional statuses. igfbp1 was not expressed in muscle during any of the nutritional conditions, whereas igfbp3 and igfbp5 were the lowest and the highest igfbps expressed respectively. A dynamic expression pattern was found in all the igfbps expressed in skeletal muscle, which depended on the nutritional status and sampling period. During the fasting period, igfbp2, 4, and 5 were downregulated, whereas igfbp3 was upregulated during part of the fasting period. The restoration of food modulated the expression of the igfbps dynamically, showing significant changes during both the long- and short-term refeeding. igfbp3 and igfbp6 were downregulated during short-term refeeding, whereas igfbp5 was upregulated, and igfbp2 and igfbp4 remained stable. During long-term refeeding, the expression of igfbp2, 4, 5, and 6 increased, while igfbp3 remained unchanged. In conclusion, this study shows for the first time the isolation of all igfbps in a single fish species, in addition to describing a dynamic nutritional and time-dependent response in the expression of igfbps in the skeletal muscle of a nonmammalian species.

Universal scaling rules predict evolutionary patterns of myogenesis in species with indeterminate growth

Intraspecific phenotypic variation is ubiquitous and often associated with resource exploitation in emerging habitats. For example, reduced body size has evolved repeatedly in Arctic charr (Salvelinus alpinus L.) and threespine stickleback (Gasterosteus aculeatus L.) across post-glacial habitats of the Northern Hemisphere. Exploiting these models, we examined how body size and myogenesis evolve with respect to the 'optimum fibre size hypothesis', which predicts that selection acts to minimize energetic costs associated with ionic homeostasis by optimizing muscle fibre production during development. In eight dwarf Icelandic Arctic charr populations, the ultimate production of fast-twitch muscle fibres (FN(max)) was only 39.5 and 15.5 per cent of that in large-bodied natural and aquaculture populations, respectively. Consequently, average fibre diameter (FD) scaled with a mass exponent of 0.19, paralleling the relaxation of diffusional constraints associated with mass-specific metabolic rate scaling. Similar reductions in FN(max) were observed for stickleback, including a small-bodied Alaskan population derived from a larger-bodied oceanic stock over a decadal timescale. The results suggest that in species showing indeterminate growth, body size evolution is accompanied by strong selection for fibre size optimization, theoretically allowing resources saved from ionic homeostasis to be allocated to other traits affecting fitness, including reproduction. Gene flow between small- and large-bodied populations residing in sympatry may counteract the evolution of this trait.