Evaluation of low-density SNP panels and imputation for cost-effective genomic selection in four aquaculture species.

Genomic selection can accelerate genetic progress in aquaculture breeding programmes, particularly for traits measured on siblings of selection candidates. However, it is not widely implemented in most aquaculture species, and remains expensive due to high genotyping costs. Genotype imputation is a promising strategy that can reduce genotyping costs and facilitate the broader uptake of genomic selection in aquaculture breeding programmes. Genotype imputation can predict ungenotyped SNPs in populations genotyped at a low-density (LD), using a reference population genotyped at a high-density (HD). In this study, we used datasets of four aquaculture species (Atlantic salmon, turbot, common carp and Pacific oyster), phenotyped for different traits, to investigate the efficacy of genotype imputation for cost-effective genomic selection. The four datasets had been genotyped at HD, and eight LD panels (300-6,000 SNPs) were generated in silico. SNPs were selected to be: i) evenly distributed according to physical position ii) selected to minimise the linkage disequilibrium between adjacent SNPs or iii) randomly selected. Imputation was performed with three different software packages (AlphaImpute2, FImpute v.3 and findhap v.4). The results revealed that FImpute v.3 was faster and achieved higher imputation accuracies. Imputation accuracy increased with increasing panel density for both SNP selection methods, reaching correlations greater than 0.95 in the three fish species and 0.80 in Pacific oyster. In terms of genomic prediction accuracy, the LD and the imputed panels performed similarly, reaching values very close to the HD panels, except in the pacific oyster dataset, where the LD panel performed better than the imputed panel. In the fish species, when LD panels were used for genomic prediction without imputation, selection of markers based on either physical or genetic distance (instead of randomly) resulted in a high prediction accuracy, whereas imputation achieved near maximal prediction accuracy independently of the LD panel, showing higher reliability. Our results suggests that, in fish species, well-selected LD panels may achieve near maximal genomic selection prediction accuracy, and that the addition of imputation will result in maximal accuracy independently of the LD panel. These strategies represent effective and affordable methods to incorporate genomic selection into most aquaculture settings.

Temperature-associated selection linked to putative chromosomal inversions in king scallop (Pecten maximus).

The genomic landscape of divergence-the distribution of differences among populations or species across the genome-is increasingly characterized to understand the role that microevolutionary forces such as natural selection and recombination play in causing and maintaining genetic divergence. This line of inquiry has also revealed chromosome structure variation to be an important factor shaping the landscape of adaptive genetic variation. Owing to a high prevalence of chromosome structure variation and the strong pressure for local adaptation necessitated by their sessile nature, bivalve molluscs are an ideal taxon for exploring the relationship between chromosome structure variation and local adaptation. Here, we report a population genomic survey of king scallop (Pecten maximus) across its natural range in the northeastern Atlantic Ocean, using a recent chromosome-level genome assembly. We report the presence of at least three large (12-22 Mb), putative chromosomal inversions associated with sea surface temperature and whose frequencies are in contrast to neutral population structure. These results highlight a potentially large role for recombination-suppressing chromosomal inversions in local adaptation and suggest a hypothesis to explain the maintenance of differences in reproductive timing found at relatively small spatial scales across king scallop populations.

Untangling Mucosal Drug Delivery: Engineering, Designing, and Testing Nanoparticles to Overcome the Mucus Barrier.

Mucus is a complex viscoelastic gel and acts as a barrier covering much of the soft tissue in the human body. High vascularization and accessibility have motivated drug delivery to various mucosal surfaces; however, these benefits are hindered by the mucus layer. To overcome the mucus barrier, many nanomedicines have been developed, with the goal of improving the efficacy and bioavailability of drug payloads. Two major nanoparticle-based strategies have emerged to facilitate mucosal drug delivery, namely, mucoadhesion and mucopenetration. Generally, mucoadhesive nanoparticles promote interactions with mucus for immobilization and sustained drug release, whereas mucopenetrating nanoparticles diffuse through the mucus and enhance drug uptake. The choice of strategy depends on many factors pertaining to the structural and compositional characteristics of the target mucus and mucosa. While there have been promising results in preclinical studies, mucus-nanoparticle interactions remain poorly understood, thus limiting effective clinical translation. This article reviews nanomedicines designed with mucoadhesive or mucopenetrating properties for mucosal delivery, explores the influence of site-dependent physiological variation among mucosal surfaces on efficacy, transport, and bioavailability, and discusses the techniques and models used to investigate mucus-nanoparticle interactions. The effects of non-homeostatic perturbations on protein corona formation, mucus composition, and nanoparticle performance are discussed in the context of mucosal delivery. The complexity of the mucosal barrier necessitates consideration of the interplay between nanoparticle design, tissue-specific differences in mucus structure and composition, and homeostatic or disease-related changes to the mucus barrier to develop effective nanomedicines for mucosal delivery.

The structural variation landscape in 492 Atlantic salmon genomes.

Structural variants (SVs) are a major source of genetic and phenotypic variation, but remain challenging to accurately type and are hence poorly characterized in most species. We present an approach for reliable SV discovery in non-model species using whole genome sequencing and report 15,483 high-confidence SVs in 492 Atlantic salmon (Salmo salar L.) sampled from a broad phylogeographic distribution. These SVs recover population genetic structure with high resolution, include an active DNA transposon, widely affect functional features, and overlap more duplicated genes retained from an ancestral salmonid autotetraploidization event than expected. Changes in SV allele frequency between wild and farmed fish indicate polygenic selection on behavioural traits during domestication, targeting brain-expressed synaptic networks linked to neurological disorders in humans. This study offers novel insights into the role of SVs in genome evolution and the genetic architecture of domestication traits, along with resources supporting reliable SV discovery in non-model species.

Demonstration of the Use of Environmental DNA for the Non-Invasive Genotyping of a Bivalve Mollusk, the European Flat Oyster (Ostrea edulis).

Accurate SNP (single nucleotide polymorphism) genotype information is critical for a wide range of selective breeding applications in aquaculture, including parentage assignment, marker-assisted, and genomic selection. However, the sampling of tissue for genetic analysis can be invasive for juvenile animals or taxa where sampling tissue is difficult or may cause mortality (e.g. bivalve mollusks). Here, we demonstrate a novel, non-invasive technique for sampling DNA based on the collection of environmental DNA using European Flat Oysters (Ostrea edulis) as an example. The live animals are placed in individual containers until sufficient genetic material is released into the seawater which is then recovered by filtration. We compared the results of tissue and eDNA derived SNP genotype calls using a PCR based genotyping platform. We found that 100% accurate genotype calls from eDNA are possible, but depend on appropriate filtration and the dilution of the sample throughout the workflow. We also developed an additional low-cost DNA extraction technique which provided >99% correct SNP genotype calls in comparison to tissue. It was concluded that eDNA sampling can be used in hatchery and selective breeding programs applicable to any aquatic organism for which direct tissue sampling may result in animal welfare concerns or mortality.

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.

Duplication of a Single myhz1.1 Gene Facilitated the Ability of Goldfish (Carassius auratus) to Alter Fast Muscle Contractile Properties With Seasonal Temperature Change.

Seasonal temperature changes markedly effect the swimming performance of some cyprinid fish acutely tested at different temperatures, involving a restructuring of skeletal muscle phenotype including changes in contractile properties and myosin heavy chain expression. We analyzed the transcriptome of fast myotomal muscle from goldfish (Carassius auratus L.) acclimated to either 8 or 25°C for 4 weeks (12 h light: 12 h dark) and identified 10 myosin heavy chains (myh) and 13 myosin light chain (myl) transcripts. Goldfish orthologs were classified based on zebrafish nomenclature as myhz1.1α, myhz1.1ß, myhz1.1γ, myha, myhb, embryo_myh1, myh9b, smyh2, symh3, and myh11 (myosin heavy chains) and myl1a, myl1b, myl2, myl9a, myl9b, myl3, myl13, myl6, myl12.1a, myl12.1b, myl12.2a, myl12.2b, and myl10 (myosin light chains). The most abundantly expressed transcripts myhz1.1α, myhz1.1ß, myhz1.1γ, myha, myl1a, myl1b, myl2, and myl3) were further investigated in fast skeletal muscle of goldfish acclimated to either 4, 8, 15, or 30°C for 12 weeks (12 h light:12 h dark). Total copy number for the myosin heavy chains showed a distinct optimum at 15°C (P < 0.01). Together myhz1.1α and myhz1.1ß comprised 90 to 97% of myhc transcripts below 15°C, but only 62% at 30°C. Whereas myhz1.1α and myhz1.1ß were equally abundant at 4 and 8°C, myhz1.1ß transcripts were 17 and 12 times higher than myhz1.1α at 15 and 30°C, respectively, (P < 0.01). Myhz1.1γ expression was at least nine-fold higher at 30°C than at cooler temperatures (P < 0.01). In contrast, the expression of myha and myosin light chains showed no consistent pattern with acclimation temperature. A phylogenetic analysis indicated that the previously reported ability of goldfish and common carp to alter contractile properties and myofibrillar ATPase activity with temperature acclimation was related to the duplication of a single myhz1.1 fast muscle myosin heavy chain found in basal cyprinids such as the zebrafish (Danio rerio).

Genomic Tools and Selective Breeding in Molluscs.

The production of most farmed molluscs, including mussels, oysters, scallops, abalone, and clams, is heavily dependent on natural seed from the plankton. Closing the lifecycle of species in hatcheries can secure independence from wild stocks and enables long-term genetic improvement of broodstock through selective breeding. Genomic techniques have the potential to revolutionize hatchery-based selective breeding by improving our understanding of the characteristics of mollusc genetics that can pose a challenge for intensive aquaculture and by providing a new suite of tools for genetic improvement. Here we review characteristics of the life history and genetics of molluscs including high fecundity, self-fertilization, high genetic diversity, genetic load, high incidence of deleterious mutations and segregation distortion, and critically assess their impact on the design and effectiveness of selective breeding strategies. A survey of the results of current breeding programs in the literature show that selective breeding with inbreeding control is likely the best strategy for genetic improvement of most molluscs, and on average growth rate can be improved by 10% per generation and disease resistance by 15% per generation across the major farmed species by implementing individual or family-based selection. Rapid advances in sequencing technology have resulted in a wealth of genomic resources for key species with the potential to greatly improve hatchery-based selective breeding of molluscs. In this review, we catalog the range of genomic resources currently available for molluscs of aquaculture interest and discuss the bottlenecks, including lack of high-quality reference genomes and the relatively high cost of genotyping, as well as opportunities for applying genomics-based selection.

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.

Comparison of the transcriptional responses of skeletal muscle and bone to a flooding dose of leucine in the gilthead sea bream (Sparus aurata).

Skeletal muscle, cartilage and bone must function in a co-ordinated fashion during locomotion and growth. In the present study on the gilthead sea bream (Sparus aurata) we tested the hypothesis that muscle and bone differ in their responsiveness to stimuli eliciting fast growth, providing a potential mechanism for generating the skeletal deformities observed in aquaculture. To investigate transcription regulation in skeletal muscle and bone we stimulated protein synthesis using a flooding dose of the branched chain amino acid leucine and compared the results with saline-injected controls. To increase the amount of available sequence information for gene expression analysis a de novo transcriptome was assembled using publicly available Next Generation Sequencing libraries from embryo, fast skeletal muscle, bone and cartilage. The resulting 5 million reads were assembled into 125,646 isotigs representing around 16,000 unique genes, including most components of the Pi3k/Akt/mTor signalling pathway. Principal components analysis was able to distinguish the transcriptional responses between leucine and saline injected controls in skeletal muscle, but not in the bone. General Linear Modelling revealed significant temporal changes in gene expression following leucine injection including the tissue-specific markers sparc, bglap (bone), mlc2 and myod2 (muscle) and gene transcripts associated with Pi3k/Akt/mTor signalling, p70sk6, akt2, ampka and mtor. Skeletal muscle showed more pronounced and rapid changes in transcript abundance than the bone to the same pro-growth signal. The observed differences in transcriptional response are consistent with the idea that fast growth results in a miss-match between muscle and bone development and may contribute to a higher incidence of skeletal deformities.

RNAseq analysis of fast skeletal muscle in restriction-fed transgenic coho salmon (Oncorhynchus kisutch): an experimental model uncoupling the growth hormone and nutritional signals regulating growth.

BACKGROUND: Coho salmon (Oncorhynchus kisutch) transgenic for growth hormone (Gh) express Gh in multiple tissues which results in increased appetite and continuous high growth with satiation feeding. Restricting Gh-transgenics to the same lower ration (TR) as wild-type fish (WT) results in similar growth, but with the recruitment of fewer, larger diameter, muscle skeletal fibres to reach a given body size. In order to better understand the genetic mechanisms behind these different patterns of muscle growth and to investigate how the decoupling of Gh and nutritional signals affects gene regulation we used RNA-seq to compare the fast skeletal muscle transcriptome in TR and WT coho salmon. RESULTS: Illumina sequencing of individually barcoded libraries from 6 WT and 6 TR coho salmon yielded 704,550,985 paired end reads which were used to construct 323,115 contigs containing 19,093 unique genes of which >10,000 contained >90 % of the coding sequence. Transcripts coding for 31 genes required for myoblast fusion were identified with 22 significantly downregulated in TR relative to WT fish, including 10 (vaspa, cdh15, graf1, crk, crkl, dock1, trio, plekho1a, cdc42a and dock5) associated with signaling through the cell surface protein cadherin. Nineteen out of 44 (43 %) translation initiation factors and 14 of 47 (30 %) protein chaperones were upregulated in TR relative to WT fish. CONCLUSIONS: TR coho salmon showed increased growth hormone transcripts and gene expression associated with protein synthesis and folding than WT fish even though net rates of protein accretion were similar. The uncoupling of Gh and amino acid signals likely results in additional costs of transcription associated with protein turnover in TR fish. The predicted reduction in the ionic costs of homeostasis in TR fish associated with increased fibre size were shown to involve multiple pathways regulating myotube fusion, particularly cadherin signaling.

Characterisation and expression analysis of cathepsins and ubiquitin-proteasome genes in gilthead sea bream (Sparus aurata) skeletal muscle.

BACKGROUND: The proteolytic enzymes involved in normal protein turnover in fish muscle are also responsible for post-mortem softening of the flesh and are therefore potential determinants of product quality. The main enzyme systems involved are calpains, cathepsins, and the ubiquitin-proteasome (UbP). In this study on Sparus aurata (Sa), the coding sequences of cathepsins (SaCTSB and SaCTSDb) and UbP family members (SaN3 and SaUb) were cloned from fast skeletal muscle, and their expression patterns were examined during ontogeny and in a fasting/re-feeding experiment. RESULTS: The amino acid sequences identified shared 66-100% overall identity with their orthologues in other vertebrates, with well conserved characteristic functional domains and catalytic residues. SaCTSDb showed phylogenetic, sequence and tissue distribution differences with respect to its paralogue SaCTSDa, previously identified in the ovary. Expression of gilthead sea bream cathepsins (B, L, Da, Db) and UbP members (N3, Ub, MuRF1 and MAFbx) in fast skeletal muscle was determined at three different life-history stages and in response to fasting and re-feeding in juveniles. Most of the proteolytic genes analysed were significantly up-regulated during fasting, and down-regulated with re-feeding and, between the fingerling (15 g) and juvenile/adult stages (~50/500 g), consistent with a decrease in muscle proteolysis in both later contexts. In contrast, SaCTSDa and SaMuRF1 expression was relatively stable with ontogeny and SaUb had higher expression in fingerlings and adults than juveniles. CONCLUSIONS: The data obtained in the present study suggest that cathepsins and UbP genes in gilthead sea bream are co-ordinately regulated during ontogeny to control muscle growth, and indicate that feeding regimes can modulate their expression, providing a potential dietary method of influencing post-mortem fillet tenderisation, and hence, product quality.

Characterization of the transcriptome of fast and slow muscle myotomal fibres in the pacu (Piaractus mesopotamicus).

BACKGROUND: The Pacu (Piaractus mesopotamicus) is a member of the Characiform family native to the Prata Basin (South America) and a target for the aquaculture industry. A limitation for the development of a selective breeding program for this species is a lack of available genetic information. The primary objectives of the present study were 1) to increase the genetic resources available for the species, 2) to exploit the anatomical separation of myotomal fibres types to compare the transcriptomes of slow and fast muscle phenotypes and 3) to systematically investigate the expression of Ubiquitin Specific Protease (USP) family members in fast and slow muscle in response to fasting and refeeding. RESULTS: We generated 0.6 Tb of pair-end reads from slow and fast skeletal muscle libraries. Over 665 million reads were assembled into 504,065 contigs with an average length of 1,334 bp and N50 = 2,772 bp. We successfully annotated nearly 47% of the transcriptome and identified around 15,000 unique genes and over 8000 complete coding sequences. 319 KEGG metabolic pathways were also annotated and 380 putative microsatellites were identified. 956 and 604 genes were differentially expressed between slow and fast skeletal muscle, respectively. 442 paralogues pairs arising from the teleost-specific whole genome duplication were identified, with the majority showing different expression patterns between fibres types (301 in slow and 245 in fast skeletal muscle). 45 members of the USP family were identified in the transcriptome. Transcript levels were quantified by qPCR in a separate fasting and refeeding experiment. USP genes in fast muscle showed a similar transient increase in expression with fasting as the better characterized E3 ubiquitin ligases. CONCLUSION: We have generated a 53-fold coverage transcriptome for fast and slow myotomal muscle in the pacu (Piaractus mesopotamicus) significantly increasing the genetic resources available for this important aquaculture species. We describe significant differences in gene expression between muscle fibre types for fundamental components of general metabolism, the Pi3k/Akt/mTor network and myogenesis, including detailed analysis of paralogue expression. We also provide a comprehensive description of USP family member expression between muscle fibre types and with changing nutritional status.

Profiling of the embryonic Atlantic halibut (Hippoglossus hippoglossus L.) transcriptome reveals maternal transcripts as potential markers of embryo quality.

BACKGROUND: Commercial Atlantic halibut (Hippoglossus hippoglossus) farming is restricted by variable oocyte quality, slow growth, and early maturation of male fish. Maternally transferred components regulate early developmental processes; therefore, they have an effect on the future viability of the embryo. Using a newly developed Agilent 10 k custom-made oligonucleotide array, we profiled components of the transcriptome involved in immune defence as well as germline and muscle development during early developmental stages: 8-cell embryos (8CS), germ ring stage (GR), 10-somite stage (10SS), and hatched embryos (HT). In addition, we identified differentially expressed transcripts in low (≤9 ± 3% hatching) and high (≥86 ± 3°% hatching) quality eggs at 8CS to identify potential maternal markers for embryo quality. RESULTS: Out of 2066 differentially expressed transcripts, 160 were identified as maternal transcripts being specifically expressed at 8CS only. Twenty transcripts were differentially expressed in 8-cell embryos between low and high quality egg groups. Several immune-related transcripts were identified as promising molecular markers of hatching success including interferon regulatory factor 7 and mhc class 2A chain. Differential expression was positively validated with quantitative real-time PCR. CONCLUSIONS: We have demonstrated maternal transfer of innate and adaptive immune system transcripts into Atlantic halibut embryos and their relation with future embryo developmental potential. We identified several transcripts as potential molecular markers of embryo quality. The developed microarray represents a useful resource for improving the commercial production of Atlantic halibut.

Muscle fibre size optimisation provides flexibility for energy budgeting in calorie-restricted coho salmon transgenic for growth hormone.

Coho salmon (Oncorhynchus kisutch) transgenic for growth hormone (GH) show substantially faster growth than wild-type (WT) fish. We fed GH-transgenic salmon either to satiation (1 year; TF) or the same smaller ration of wild-type fish (2 years; TR), resulting in groups matched for body size to WT salmon. The myotomes of TF and WT fish had the same number and size distribution of muscle fibres, indicating a twofold higher rate of fibre recruitment in the GH transgenics. Unexpectedly, calorie restriction was found to decrease the rate of fibre production in transgenics, resulting in a 20% increase in average fibre size and reduced costs of ionic homeostasis. Genes for myotube formation were downregulated in TR relative to TF and WT fish. We suggest that muscle fibre size optimisation allows the reallocation of energy from maintenance to locomotion, explaining the observation that calorie-restricted transgenics grow at the same rate as WT fish whilst exhibiting markedly higher foraging activity.

Cardiac myoglobin deficit has evolved repeatedly in teleost fishes.

Myoglobin (Mb) is the classic vertebrate oxygen-binding protein present in aerobic striated muscles. It functions principally in oxygen delivery and provides muscle with its characteristic red colour. Members of the Antarctic icefish family (Channichthyidae) are widely thought to be extraordinary for lacking cardiac Mb expression, a fact that has been attributed to their low metabolic rate and unusual evolutionary history. Here, we report that cardiac Mb deficit, associated with pale heart colour, has evolved repeatedly during teleost evolution. This trait affects both gill- and air-breathing species from temperate to tropical habitats across a full range of salinities. Cardiac Mb deficit results from total pseudogenization in three-spined stickleback and is associated with a massive reduction in mRNA level in two species that evidently retain functional Mb. The results suggest that near or complete absence of Mb-assisted oxygen delivery to heart muscle is a common facet of teleost biodiversity, even affecting lineages with notable oxygen demands. We suggest that Mb deficit may affect how different teleost species deal with increased tissue oxygen demands arising under climate change.

Systematic variation in the pattern of gene paralog retention between the teleost superorders Ostariophysi and Acanthopterygii.

Teleost fish underwent whole-genome duplication around 450 Ma followed by diploidization and loss of 80-85% of the duplicated genes. To identify a deep signature of this teleost-specific whole-genome duplication (TSGD), we searched for duplicated genes that were systematically and uniquely retained in one or other of the superorders Ostariophysi and Acanthopterygii. TSGD paralogs comprised 17-21% of total gene content. Some 2.6% (510) of TSGD paralogs were present as pairs in the Ostariophysi genomes of Danio rerio (Cypriniformes) and Astyanax mexicanus (Characiformes) but not in species from four orders of Acanthopterygii (Gasterosteiformes, Gasterosteus aculeatus; Tetraodontiformes, Tetraodon nigroviridis; Perciformes, Oreochromis niloticus; and Beloniformes, Oryzias latipes) where a single copy was identified. Similarly, 1.3% (418) of total gene number represented cases where TSGD paralogs pairs were systematically retained in the Acanthopterygian but conserved as a single copy in Ostariophysi genomes. We confirmed the generality of these results by phylogenetic and synteny analysis of 40 randomly selected linage-specific paralogs (LSPs) from each superorder and completed with the transcriptomes of three additional Ostariophysi species (Ictalurus punctatus [Siluriformes], Sinocyclocheilus species [Cypriniformes], and Piaractus mesopotamicus [Characiformes]). No chromosome bias was detected in TSGD paralog retention. Gene ontology (GO) analysis revealed significant enrichment of GO terms relative to the human GO SLIM database for "growth," "Cell differentiation," and "Embryo development" in Ostariophysi and for "Transport," "Signal Transduction," and "Vesicle mediated transport" in Acanthopterygii. The observed patterns of paralog retention are consistent with different diploidization outcomes having contributed to the evolution/diversification of each superorder.

A well-constrained estimate for the timing of the salmonid whole genome duplication reveals major decoupling from species diversification.

Whole genome duplication (WGD) is often considered to be mechanistically associated with species diversification. Such ideas have been anecdotally attached to a WGD at the stem of the salmonid fish family, but remain untested. Here, we characterized an extensive set of gene paralogues retained from the salmonid WGD, in species covering the major lineages (subfamilies Salmoninae, Thymallinae and Coregoninae). By combining the data in calibrated relaxed molecular clock analyses, we provide the first well-constrained and direct estimate for the timing of the salmonid WGD. Our results suggest that the event occurred no later in time than 88 Ma and that 40-50 Myr passed subsequently until the subfamilies diverged. We also recovered a Thymallinae-Coregoninae sister relationship with maximal support. Comparative phylogenetic tests demonstrated that salmonid diversification patterns are closely allied in time with the continuous climatic cooling that followed the Eocene-Oligocene transition, with the highest diversification rates coinciding with recent ice ages. Further tests revealed considerably higher speciation rates in lineages that evolved anadromy--the physiological capacity to migrate between fresh and seawater--than in sister groups that retained the ancestral state of freshwater residency. Anadromy, which probably evolved in response to climatic cooling, is an established catalyst of genetic isolation, particularly during environmental perturbations (for example, glaciation cycles). We thus conclude that climate-linked ecophysiological factors, rather than WGD, were the primary drivers of salmonid diversification.

Characterisation and expression of myogenesis regulatory factors during in vitro myoblast development and in vivo fasting in the gilthead sea bream (Sparus aurata).

The aim of this study was to characterise a primary cell culture isolated from fast skeletal muscle of the gilthead sea bream. Gene expression profiles during culture maturation were compared with those obtained from a fasting-refeeding model which is widely used to modulate myogenesis in vivo. Myogenesis is controlled by numerous extracellular signals together with intracellular transcriptional factors whose coordinated expression is critical for the appropriate development of muscle fibres. Full-length cDNAs for the transcription factors Myf5, Mrf4, Pax7 and Sox8 were cloned and sequenced for gilthead sea bream. Pax7, sox8, myod2 and myf5 levels were up-regulated during the proliferating phase of the myogenic cultures coincident with the highest expression of proliferating cell nuclear antigen (PCNA). In contrast, myogenin and mrf4 transcript abundance was highest during the differentiation phase of the culture when myotubes were present, and was correlated with increased myosin heavy chain (mhc) and desmin expression. In vivo, 30days of fasting resulted in muscle fibre atrophy, a reduction in myod2, myf5 and igf1 expression, lower number of Myod-positive cells, and decreased PCNA protein expression, whereas myogenin expression was not significantly affected. Myostatin1 (mstn1) and pax7 expression were up-regulated in fasted relative to well-fed individuals, consistent with a role for Pax7 in the reduction of myogenic cell activity with fasting. The primary cell cultures and fasting-feeding experiments described provide a foundation for the future investigations on the regulation of muscle growth in gilthead sea bream.

Characterisation and expression of calpain family members in relation to nutritional status, diet composition and flesh texture in gilthead sea bream (Sparus aurata).

Calpains are non-lysosomal calcium-activated neutral proteases involved in a wide range of cellular processes including muscle proteolysis linked to post-mortem flesh softening. The aims of this study were (a) to characterise several members of the calpain system in gilthead sea bream and (b) to examine their expression in relation to nutritional status and muscle tenderisation. We identified the complete open reading frame of gilthead sea bream calpains1-3, sacapn1, sacapn2, sacapn3, and two paralogs of the calpain small subunit1, sacapns1a and sacapns1b. Proteins showed 63-90% sequence identity compared with sequences from mammals and other teleost fishes, and the characteristic domain structure of vertebrate calpains. Transcripts of sacapn1, sacapn2, sacapns1a and sacapns1b had a wide tissue distribution, whereas sacapn3 was almost exclusively detected in skeletal muscle. Next, we assessed transcript expression in skeletal muscle following alteration of nutritional status by (a) fasting and re-feeding or (b) feeding four experimental diets with different carbohydrate-to-protein ratios. Fasting significantly reduced plasma glucose and increased free fatty acids and triglycerides, together with a significant increase in sacapns1b expression. Following 7 days of re-feeding, plasma parameters returned to fed values and sacapn1, sacapn2, sacapns1a and sacapns1b expression was significantly reduced. Furthermore, an increase in dietary carbohydrate content (11 to 39%) diminished growth but increased muscle texture, which showed a significant correlation with decreased sacapn1 and sacapns1a expression, whilst the other calpains remained unaffected. This study has demonstrated that calpain expression is modulated by nutritional status and diet composition in gilthead sea bream, and that the expression of several calpain members is correlated with muscle texture, indicating their potential use as molecular markers for flesh quality in aquaculture production.