Development of genomic markers associated to growth-related traits and sex determination in lumpfish (Cyclopterus lumpus).

Cleaner fish species have gained great importance in the control of sea lice, among them, lumpfish (Cyclopterus lumpus) has become one of the most popular. Lumpfish life cycle has been closed, and hatchery reproduction is now possible, however, current production is reliant on wild caught broodstock to meet the increasing demand. Selective breeding practices are called to play an important role in the successful breeding of most aquaculture species, including lumpfish. In this study we analysed a lumpfish population for the identification of genomic markers linked to production traits. Sequencing of RAD libraries allowed us to identify, 7193 informative markers within the sampled individuals. Genome wide association analysis for sex, weight, condition factor and standard length was performed. One single major QTL region was identified for sex, while nine QTL regions were detected for weight, and three QTL regions for standard length. A total of 177 SNP markers of interest (from QTL regions) and 399 high Fst SNP markers were combined in a low-density panel, useful to obtain relevant genetic information from lumpfish populations. Moreover, a robust combined subset of 29 SNP markers (10 associated to sex, 14 to weight and 18 to standard length) provided over 90% accuracy in predicting the animal's phenotype by machine learning. Overall, our findings provide significant insights into the genetic control of important traits in lumpfish and deliver important genomic resources that will facilitate the establishment of selective breeding programmes in lumpfish.

Tracking emissions in the US electricity system.

Understanding electricity consumption and production patterns is a necessary first step toward reducing the health and climate impacts of associated emissions. In this work, the economic input-output model is adapted to track emissions flows through electric grids and quantify the pollution embodied in electricity production, exchanges, and, ultimately, consumption for the 66 continental US Balancing Authorities (BAs). The hourly and BA-level dataset we generate and release leverages multiple publicly available datasets for the year 2016. Our analysis demonstrates the importance of considering location and temporal effects as well as electricity exchanges in estimating emissions footprints. While increasing electricity exchanges makes the integration of renewable electricity easier, importing electricity may also run counter to climate-change goals, and citizens in regions exporting electricity from high-emission-generating sources bear a disproportionate air-pollution burden. For example, 40% of the carbon emissions related to electricity consumption in California's main BA were produced in a different region. From 30 to 50% of the sulfur dioxide and nitrogen oxides released in some of the coal-heavy Rocky Mountain regions were related to electricity produced that was then exported. Whether for policymakers designing energy efficiency and renewable programs, regulators enforcing emissions standards, or large electricity consumers greening their supply, greater resolution is needed for electric-sector emissions indices to evaluate progress against current and future goals.

Considering adaptive genetic variation in climate change vulnerability assessment reduces species range loss projections.

Local adaptations can determine the potential of populations to respond to environmental changes, yet adaptive genetic variation is commonly ignored in models forecasting species vulnerability and biogeographical shifts under future climate change. Here we integrate genomic and ecological modeling approaches to identify genetic adaptations associated with climate in two cryptic forest bats. We then incorporate this information directly into forecasts of range changes under future climate change and assessment of population persistence through the spread of climate-adaptive genetic variation (evolutionary rescue potential). Considering climate-adaptive potential reduced range loss projections, suggesting that failure to account for intraspecific variability can result in overestimation of future losses. On the other hand, range overlap between species was projected to increase, indicating that interspecific competition is likely to play an important role in limiting species' future ranges. We show that although evolutionary rescue is possible, it depends on a population's adaptive capacity and connectivity. Hence, we stress the importance of incorporating genomic data and landscape connectivity in climate change vulnerability assessments and conservation management.

The nedd-8 activating enzyme gene underlies genetic resistance to infectious pancreatic necrosis virus in Atlantic salmon.

Genetic resistance to infectious pancreatic necrosis virus (IPNV) in Atlantic salmon is a rare example of a trait where a single locus (QTL) explains almost all of the genetic variation. Genetic marker tests based on this QTL on salmon chromosome 26 have been widely applied in selective breeding to markedly reduce the incidence of the disease. In the current study, whole genome sequencing and functional annotation approaches were applied to characterise genes and variants in the QTL region. This was complemented by an analysis of differential expression between salmon fry of homozygous resistant and homozygous susceptible genotypes challenged with IPNV. These analyses pointed to the NEDD-8 activating enzyme 1 (nae1) gene as a putative functional candidate underlying the QTL effect. The role of nae1 in IPN resistance was further assessed via CRISPR-Cas9 knockout of the nae1 gene and chemical inhibition of the nae1 protein activity in Atlantic salmon cell lines, both of which resulted in highly significant reduction in productive IPNV replication. In contrast, CRISPR-Cas9 knockout of a candidate gene previously purported to be a cellular receptor for the virus (cdh1) did not have a major impact on productive IPNV replication. These results suggest that nae1 is the causative gene underlying the major QTL affecting resistance to IPNV in salmon, provide further evidence for the critical role of neddylation in host-pathogen interactions, and highlight the value in combining high-throughput genomics approaches with targeted genome editing to understand the genetic basis of disease resistance.

DNA polymorphism underlying allozyme variation at a malic enzyme locus (mMEP-2*) in Atlantic salmon (Salmo salar L.).

A non-synonymous single nucleotide polymorphism (SNP) underlies a diallelic allozyme polymorphism at the mitochondrial NADP-dependent mMEP-2* locus in Atlantic salmon (Salmo salar L.). The resultant amino acid substitution, which alters the charge of the allelic products, matches the differential mobility of the two allozyme alleles, whereas allozyme and SNP assays revealed genotyping concordance in 257 of 258 individuals. A single mismatch, homozygous allozyme vs. heterozygote SNP, suggests the presence of a second, less common null allele.

A high-density linkage map and sex-linked markers for the Amazon Tambaqui Colossoma macropomum.

BACKGROUND: Tambaqui (Colossoma macropomum, Cuvier, 1818) is the most economically important native freshwater fish species in Brazil. It can reach a total length of over 1 m and a weight of over 40 kg. The species displays a clear sex dimorphism in growth performance, with females reaching larger sizes at harvest. In aquaculture, the production of monosex populations in selective breeding programmes has been therefore identified as a key priority. RESULTS: In the present study, a genetic linkage map was generated by double digest restriction-site associated DNA (ddRAD) sequencing from 248 individuals sampled from two F1 families. The map was constructed using 14,805 informative SNPs and spanned 27 linkage groups. From this, the tambaqui draft genome was improved, by ordering the scaffolds into chromosomes, and sex-linked markers were identified. A total of 235 markers on linkage group 26 showed a significant association with the phenotypic sex, supporting an XX/XY sex determination system in the species. The four most informative sex-linked markers were validated on another 206 sexed individuals, demonstrating an accuracy in predicting sex ranging from 90.0 to 96.7%. CONCLUSIONS: The genetic mapping and novel sex-linked DNA markers identified and validated offer new tools for rapid progeny sexing, thus supporting the development of monosex female production in the industry while also supporting breeding programmes of the species.

Transcriptomic comparison of communally reared wild, domesticated and hybrid Atlantic salmon fry under stress and control conditions.

BACKGROUND: Domestication is the process by which organisms become adapted to the human-controlled environment. Since the selection pressures that act upon cultured and natural populations differ, adaptations that favour life in the domesticated environment are unlikely to be advantageous in the wild. Elucidation of the differences between wild and domesticated Atlantic salmon may provide insights into some of the genomic changes occurring during domestication, and, help to predict the evolutionary consequences of farmed salmon escapees interbreeding with wild conspecifics. In this study the transcriptome of the offspring of wild and domesticated Atlantic salmon were compared using a common-garden experiment under standard hatchery conditions and in response to an applied crowding stressor. RESULTS: Transcriptomic differences between wild and domesticated crosses were largely consistent between the control and stress conditions, and included down-regulation of environmental information processing, immune and nervous system pathways and up-regulation of genetic information processing, carbohydrate metabolism, lipid metabolism and digestive and endocrine system pathways in the domesticated fish relative to their wild counterparts, likely reflective of different selection pressures acting in wild and cultured populations. Many stress responsive functions were also shared between crosses and included down-regulation of cellular processes and genetic information processing and up-regulation of some metabolic pathways, lipid and energy in particular. The latter may be indicative of mobilization and reallocation of energy resources in response to stress. However, functional analysis indicated that a number of pathways behave differently between domesticated and wild salmon in response to stress. Reciprocal F1 hybrids permitted investigation of inheritance patterns that govern transcriptomic differences between these genetically divergent crosses. Additivity and maternal dominance accounted for approximately 42 and 25% of all differences under control conditions for both hybrids respectively. However, the inheritance of genes differentially expressed between crosses under stress was less consistent between reciprocal hybrids, potentially reflecting maternal environmental effects. CONCLUSION: We conclude that there are transcriptomic differences between the domesticated and wild salmon strains studied here, reflecting the different selection pressures operating on them. Our results indicate that stress may affect certain biological functions differently in wild, domesticated and hybrid crosses and these should be further investigated.

Sex determination in the GIFT strain of tilapia is controlled by a locus in linkage group 23.

BACKGROUND: Tilapias (Family Cichlidae) are the second most important group of aquaculture species in the world. They have been the subject of much research on sex determination due to problems caused by early maturation in culture and their complex sex-determining systems. Different sex-determining loci (linkage group 1, 20 and 23) have been detected in various tilapia stocks. The 'genetically improved farmed tilapia' (GIFT) stock, founded from multiple Nile tilapia (Oreochromis niloticus) populations, with some likely to have been introgressed with O. mossambicus, is a key resource for tilapia aquaculture. The sex-determining mechanism in the GIFT stock was unknown, but potentially complicated due to its multiple origins. RESULTS: A bulk segregant analysis (BSA) version of double-digest restriction-site associated DNA sequencing (BSA-ddRADseq) was developed and used to detect and position sex-linked single nucleotide polymorphism (SNP) markers in 19 families from the GIFT strain breeding nucleus and two Stirling families as controls (a single XY locus had been previously mapped to LG1 in the latter). About 1500 SNPs per family were detected across the genome. Phenotypic sex in Stirling families showed strong association with LG1, whereas only SNPs located in LG23 showed clear association with sex in the majority of the GIFT families. No other genomic regions linked to sex determination were apparent. This region was validated using a series of LG23-specific DNA markers (five SNPs with highest association to sex from this study, the LG23 sex-associated microsatellite UNH898 and ARO172, and the recently isolated amhy marker for individual fish (n = 284). CONCLUSIONS: Perhaps surprisingly given its multiple origins, sex determination in the GIFT strain breeding nucleus was associated only with a locus in LG23. BSA-ddRADseq allowed cost-effective analysis of multiple families, strengthening this conclusion. This technique has potential to be applied to other complex traits. The sex-linked SNP markers identified will be useful for potential marker-assisted selection (MAS) to control sex-ratio in GIFT tilapia to suppress unwanted reproduction during growout.

Genetic diversity and structure in Arapaima gigas populations from Amazon and Araguaia-Tocantins river basins.

BACKGROUND: Arapaima gigas (Schinz, 1822) is the largest freshwater scaled fish in the world, and an emerging species for tropical aquaculture development. Conservation of the species, and the expansion of aquaculture requires the development of genetic tools to study polymorphism, differentiation, and stock structure. This study aimed to investigate genomic polymorphism through ddRAD sequencing, in order to identify a panel of single nucleotide polymorphisms (SNPs) and to simultaneously assess genetic diversity and structure in wild (from rivers Amazon, Solimões, Tocantins and Araguaia) and captive populations. RESULTS: Compared to many other teleosts, the degree of polymorphism in A. gigas was low with only 2.3% of identified RAD-tags (135 bases long) containing SNPs. A panel of 393 informative SNPs was identified and screened across the five populations. Higher genetic diversity indices (number of polymorphic loci and private alleles, Shannon's Index and HO) were found in populations from the Amazon and Solimões, intermediate levels in Tocantins and Captive, and very low levels in the Araguaia population. These results likely reflect larger population sizes from less urbanized environments in the Amazon basin compared to Araguaia. Populations were significantly differentiated with pairwise FST values ranging from 0.086 (Amazon × Solimões) to 0.556 (Amazon × Araguaia). Mean pairwise relatedness among individuals was significant in all populations (P < 0.01), reflecting a degree of inbreeding possibly due to severe depletion of natural stocks, the species sedentary behaviour and possible sampling biases. Although Mantel test was not significant (P = 0.104; R2 = 0.65), Bayesian analysis in STRUCTURE and discriminant analysis of principal components (DAPC) showed populations of Amazon and Solimões to be genetically differentiated from Araguaia, with Tocantins comprising individuals from both identified stocks. CONCLUSIONS: This relatively rapid genotyping by sequencing approach proved to be successful in delineating arapaima stocks. The approach and / or SNP panels identified should prove valuable for more detailed genetic studies of arapaima populations, including the elucidation of the genetic status of described discrete morphotypes and aid in delivery of conservation programs to maintain genetic diversity in reservoirs across the Amazon region.

Molecular epidemiological study on Infectious Pancreatic Necrosis Virus isolates from aquafarms in Scotland over three decades.

In order to obtain an insight into genomic changes and associated evolution and adaptation of Infectious Pancreatic Necrosis Virus (IPNV), the complete coding genomes of 57 IPNV isolates collected from Scottish aquafarms from 1982 to 2014 were sequenced and analysed. Phylogenetic analysis of the sequenced IPNV strains showed separate clustering of genogroups I, II, III and V. IPNV isolates with genetic reassortment of segment A/B of genogroup III/II were determined. About 59 % of the IPNV isolates belonged to the persistent type and 32 % to the low-virulent type, and only one highly pathogenic strain (1.79 %) was identified. Codon adaptation index calculations indicated that the IPNV major capsid protein VP2 has adapted to its salmonid host. Under-representation of CpG dinucleotides in the IPNV genome to minimize detection by the innate immunity receptors, and observed positive selection in the virulence determination sites of VP2 embedded in the variable region of the main antigenic region, suggest an immune escape mechanism driving virulence evolution. The prevalence of mostly persistent genotypes, together with the assumption of adaptation and immune escape, indicates that IPNV is evolving with the host.

Whole genome duplication and transposable element proliferation drive genome expansion in Corydoradinae catfishes.

Genome size varies significantly across eukaryotic taxa and the largest changes are typically driven by macro-mutations such as whole genome duplications (WGDs) and proliferation of repetitive elements. These two processes may affect the evolutionary potential of lineages by increasing genetic variation and changing gene expression. Here, we elucidate the evolutionary history and mechanisms underpinning genome size variation in a species-rich group of Neotropical catfishes (Corydoradinae) with extreme variation in genome size-0.6 to 4.4 pg per haploid cell. First, genome size was quantified in 65 species and mapped onto a novel fossil-calibrated phylogeny. Two evolutionary shifts in genome size were identified across the tree-the first between 43 and 49 Ma (95% highest posterior density (HPD) 36.2-68.1 Ma) and the second at approximately 19 Ma (95% HPD 15.3-30.14 Ma). Second, restriction-site-associated DNA (RAD) sequencing was used to identify potential WGD events and quantify transposable element (TE) abundance in different lineages. Evidence of two lineage-scale WGDs was identified across the phylogeny, the first event occurring between 54 and 66 Ma (95% HPD 42.56-99.5 Ma) and the second at 20-30 Ma (95% HPD 15.3-45 Ma) based on haplotype numbers per contig and between 35 and 44 Ma (95% HPD 30.29-64.51 Ma) and 20-30 Ma (95% HPD 15.3-45 Ma) based on SNP read ratios. TE abundance increased considerably in parallel with genome size, with a single TE-family (TC1-IS630-Pogo) showing several increases across the Corydoradinae, with the most recent at 20-30 Ma (95% HPD 15.3-45 Ma) and an older event at 35-44 Ma (95% HPD 30.29-64.51 Ma). We identified signals congruent with two WGD duplication events, as well as an increase in TE abundance across different lineages, making the Corydoradinae an excellent model system to study the effects of WGD and TEs on genome and organismal evolution.

Gene-centromere mapping in meiotic gynogenetic European seabass.

BACKGROUND: Fully isogenic lines in fish can be developed using "mitotic" gynogenesis (suppression of first zygotic mitosis following inactivation of the sperm genome). However, genome-wide verification of the steps in this process has seldom been applied. We used ddRADseq to generate SNP markers in a meiotic gynogenetic family of European seabass (Dicentrarchus labrax): (i) to verify the lack of paternal contribution in a meiotic gynogenetic family; (ii) to generate a gene-centromere map from this family; (iii) to identify telomeric markers that could distinguish mitotic gynogenetics from meiotic gynogenetics, which sometimes arise spontaneously in mitotic gynogenetic families. RESULTS: From a single meiotic gynogenetic family consisting of 79 progeny, 42 million sequencing reads (Illumina, trimmed to 148 bases) resolved 6866 unique RAD-tags. The 340 male-informative SNP markers that were identified confirmed the lack of paternal contribution. A gene-centromere map was constructed based on 804 female-informative SNPs in 24 linkage groups (2n = 48) with a total length of 1251.02 cM (initial LG assignment was based on the seabass genome assembly, dicLab v1). Chromosome arm structure could be clearly discerned from the pattern of heterozygosity in each linkage group in 18 out of 24 LGs: the other six showed anomalies that appeared to be related to issues in the genome assembly. CONCLUSION: Genome-wide screening enabled substantive verification of the production of the gynogenetic family used in this study. The large number of telomeric and subtelomeric markers with high heterozygosity values in the meiotic gynogenetic family indicate that such markers could be used to clearly distinguish between meiotic and mitotic gynogenetics.

Mapping the sex determination locus in the hapuku (Polyprion oxygeneios) using ddRAD sequencing.

BACKGROUND: Hapuku (Polyprion oxygeneios) is a member of the wreckfish family (Polyprionidae) and is highly regarded as a food fish. Although adults grow relatively slowly, juveniles exhibit low feed conversion ratios and can reach market size in 1-2 years, making P. oxygeneios a strong candidate for aquaculture. However, they can take over 5 years to reach sexual maturity in captivity and are not externally sexually dimorphic, complicating many aspects of broodstock management. Understanding the sex determination system of P. oxygeneios and developing accurate assays to assign genetic sex will contribute significantly towards its full-scale commercialisation. RESULTS: DNA from parents and sexed offspring (n = 57) from a single family of captive bred P. oxygeneios was used as a template for double digestion Restriction-site Associated DNA (ddRAD) sequencing. Two libraries were constructed using SbfI - SphI and SbfI - NcoI restriction enzyme combinations, respectively. Two runs on an Illumina MiSeq platform generated 70,266,464 raw reads, identifying 19,669 RAD loci. A combined sex linkage map (1367 cM) was constructed based on 1575 Single Nucleotide Polymorphism (SNP) markers that resolved into 35 linkage groups. Sex-specific linkage maps were of similar size (1132 and 1168 cM for male and female maps respectively). A single major sex-determining locus, found to be heterogametic in males, was mapped to linkage group 14. Several markers were found to be in strong linkage disequilibrium with the sex-determining locus. Allele-specific PCR assays were developed for two of these markers, SphI6331 and SphI8298, and demonstrated to accurately differentiate sex in progeny within the same pedigree. Comparative genomic analyses indicated that many of the linkage groups within the P. oxygeneios map share a relatively high degree of homology with those published for the European seabass (Dicentrarchus labrax). CONCLUSION: P. oxygeneios has an XX/XY sex determination system. Evaluation of allele-specific PCR assays, based on the two SNP markers most closely associated with phenotypic sex, indicates that a simple molecular assay for sexing P. oxygeneios should be readily attainable. The high degree of synteny observed with D. labrax should aid further molecular genetic study and exploitation of hapuku as a food fish.

Gene expression comparison of resistant and susceptible Atlantic salmon fry challenged with Infectious Pancreatic Necrosis virus reveals a marked contrast in immune response.

BACKGROUND: Infectious Pancreatic Necrosis (IPN) is a highly contagious birnavirus disease of farmed salmonid fish, which often causes high levels of morbidity and mortality. A large host genetic component to resistance has been previously described for Atlantic salmon (Salmo salar L.), which mediates high mortality rates in some families and zero mortality in others. However, the molecular and immunological basis for this resistance is not yet fully known. This manuscript describes a global comparison of the gene expression profiles of resistant and susceptible Atlantic salmon fry following challenge with the IPN virus. RESULTS: Salmon fry from two IPNV-resistant and two IPNV-susceptible full sibling families were challenged with the virus and sampled at 1 day, 7 days and 20 days post-challenge. Significant viral titre was observed in both resistant and susceptible fish at all timepoints, although generally at higher levels in susceptible fish. Gene expression profiles combined with gene ontology and pathway analyses demonstrated that while a clear immune response was observed in both resistant and susceptible fish, there were striking differences between the two phenotypes. The susceptible fish showed marked up-regulation of genes related to cytokine activity and inflammatory response that evidently failed to protect against the virus. In contrast, the resistant fish demonstrated a less pronounced immune response including up-regulation of genes relating to the M2 macrophage system. CONCLUSIONS: While only the susceptible phenotype shows appreciable mortality levels, both resistant and susceptible fish can become infected with IPNV. Susceptible fish are characterized by a much larger, yet ineffective, immune response, largely related to cytokine and inflammatory systems. Resistant fish demonstrate a more moderate, putative macrophage-mediated inflammatory response, which may contribute to their survival.

Comparing the transcriptomes of embryos from domesticated and wild Atlantic salmon (Salmo salar L.) stocks and examining factors that influence heritability of gene expression.

BACKGROUND: Due to selective breeding, domesticated and wild Atlantic salmon are genetically diverged, which raises concerns about farmed escapees having the potential to alter the genetic composition of wild populations and thereby disrupting local adaptation. Documenting transcriptional differences between wild and domesticated stocks under controlled conditions is one way to explore the consequences of domestication and selection. We compared the transcriptomes of wild and domesticated Atlantic salmon embryos, by using a custom 44k oligonucleotide microarray to identify perturbed gene pathways between the two stocks, and to document the inheritance patterns of differentially-expressed genes by examining gene expression in their reciprocal hybrids. RESULTS: Data from 24 array interrogations were analysed: four reciprocal cross types (W♀ × W♂, D♀ × W♂; W♀ × D♂, D♀ × D♂) × six biological replicates. A common set of 31,491 features on the microarrays passed quality control, of which about 62 % were assigned a KEGG Orthology number. A total of 6037 distinct genes were identified for gene-set enrichment/pathway analysis. The most highly enriched functional groups that were perturbed between the two stocks were cellular signalling and immune system, ribosome and RNA transport, and focal adhesion and gap junction pathways, relating to cell communication and cell adhesion molecules. Most transcripts that were differentially expressed between the stocks were governed by additive gene interaction (33 to 42 %). Maternal dominance and over-dominance were also prevalent modes of inheritance, with no convincing evidence for a stock effect. CONCLUSIONS: Our data indicate that even at this relatively early developmental stage, transcriptional differences exist between the two stocks and affect pathways that are relevant to wild versus domesticated environments. Many of the identified differentially perturbed pathways are involved in organogenesis, which is expected to be an active process at the eyed egg stage. The dominant effects are more largely due to the maternal line than to the origin of the stock. This finding is particularly relevant in the context of potential introgression between farmed and wild fish, since female escapees tend to have a higher spawning success rate compared to males.

Genomic prediction of host resistance to sea lice in farmed Atlantic salmon populations.

BACKGROUND: Sea lice have significant negative economic and welfare impacts on marine Atlantic salmon farming. Since host resistance to sea lice has a substantial genetic component, selective breeding can contribute to control of lice. Genomic selection uses genome-wide marker information to predict breeding values, and can achieve markedly higher accuracy than pedigree-based methods. Our aim was to assess the genetic architecture of host resistance to sea lice, and test the utility of genomic prediction of breeding values. Individual lice counts were measured in challenge experiments using two large Atlantic salmon post-smolt populations from a commercial breeding programme, which had genotypes for ~33 K single nucleotide polymorphisms (SNPs). The specific objectives were to: (i) estimate the heritability of host resistance; (ii) assess its genetic architecture by performing a genome-wide association study (GWAS); (iii) assess the accuracy of predicted breeding values using varying SNP densities (0.5 to 33 K) and compare it to that of pedigree-based prediction; and (iv) evaluate the accuracy of prediction in closely and distantly related animals. RESULTS: Heritability of host resistance was significant (0.22 to 0.33) in both populations using either pedigree or genomic relationship matrices. The GWAS suggested that lice resistance is a polygenic trait, and no genome-wide significant quantitative trait loci were identified. Based on cross-validation analysis, genomic predictions were more accurate than pedigree-based predictions for both populations. Although prediction accuracies were highest when closely-related animals were used in the training and validation sets, the benefit of having genomic-versus pedigree-based predictions within a population increased as the relationships between training and validation sets decreased. Prediction accuracy reached an asymptote with a SNP density of ~5 K within populations, although higher SNP density was advantageous for cross-population prediction. CONCLUSIONS: Host resistance to sea lice in farmed Atlantic salmon has a significant genetic component. Phenotypes relating to host resistance can be predicted with moderate to high accuracy within populations, with a major advantage of genomic over pedigree-based methods, even at relatively sparse SNP densities. Prediction accuracies across populations were low, but improved with higher marker densities. Genomic selection can contribute to lice control in salmon farming.

Development and validation of a mixed-tissue oligonucleotide DNA microarray for Atlantic bluefin tuna, Thunnus thynnus (Linnaeus, 1758).

BACKGROUND: The largest of the tuna species, Atlantic bluefin tuna (Thunnus thynnus), inhabits the North Atlantic Ocean and the Mediterranean Sea and is considered to be an endangered species, largely a consequence of overfishing. T. thynnus aquaculture, referred to as fattening or farming, is a capture based activity dependent on yearly renewal from the wild. Thus, the development of aquaculture practices independent of wild resources can provide an important contribution towards ensuring security and sustainability of this species in the longer-term. The development of such practices is today greatly assisted by large scale transcriptomic studies. RESULTS: We have used pyrosequencing technology to sequence a mixed-tissue normalised cDNA library, derived from adult T. thynnus. A total of 976,904 raw sequence reads were assembled into 33,105 unique transcripts having a mean length of 893 bases and an N50 of 870. Of these, 33.4% showed similarity to known proteins or gene transcripts and 86.6% of them were matched to the congeneric Pacific bluefin tuna (Thunnus orientalis) genome, compared to 70.3% for the more distantly related Nile tilapia (Oreochromis niloticus) genome. Transcript sequences were used to develop a novel 15 K Agilent oligonucleotide DNA microarray for T. thynnus and comparative tissue gene expression profiles were inferred for gill, heart, liver, ovaries and testes. Functional contrasts were strongest between gills and ovaries. Gills were particularly associated with immune system, signal transduction and cell communication, while ovaries displayed signatures of glycan biosynthesis, nucleotide metabolism, transcription, translation, replication and repair. CONCLUSIONS: Sequence data generated from a novel mixed-tissue T. thynnus cDNA library provide an important transcriptomic resource that can be further employed for study of various aspects of T. thynnus ecology and genomics, with strong applications in aquaculture. Tissue-specific gene expression profiles inferred through the use of novel oligo-microarray can serve in the design of new and more focused transcriptomic studies for future research of tuna physiology and assessment of the welfare in a production environment.

A comparative analysis of the response of the hepatic transcriptome to dietary docosahexaenoic acid in Atlantic salmon (Salmo salar) post-smolts.

BACKGROUND: The present study aimed to explore the impact of dietary docosahexaenoic acid (DHA) on aspects of the metabolism of Atlantic salmon (Salmo salar). The effects of diets containing increasing levels of DHA (1 g kg(-1), 3 g kg(-1), 6 g kg(-1), 10 g kg(-1) and 13 g kg(-1)) on the liver transcriptome of post-smolt salmon was examined to elucidate patterns of gene expression and responses of specific metabolic pathways. Total RNA was isolated from the liver of individual fish and analyzed using a custom gene expression 44K feature Atlantic salmon oligo-microarray. RESULTS: The expression of up to 911 unique annotated genes was significantly affected by dietary DHA inclusion relative to a low DHA reference diet. Analysis of a total of 797 unique genes were found with a significant linear correlation between expression level and dietary DHA. Gene-Set Enrichment Analysis (GSEA) identified a range of pathways that were significantly affected by dietary DHA content. CONCLUSIONS: Pathways that showed a significant response to dietary DHA level included those for long-chain polyunsaturated fatty acid biosynthesis, fatty acid elongation, steroid biosynthesis, glycan biosynthesis, protein export and protein processing in the endoplasmic reticulum. These findings suggest that in addition to clear roles in influencing lipid metabolic pathways, DHA might also have key functional roles in other pathways distinct from lipid metabolism.

A novel sex-determining QTL in Nile tilapia (Oreochromis niloticus).

BACKGROUND: Fish species often exhibit significant sexual dimorphism for commercially important traits. Accordingly, the control of phenotypic sex, and in particular the production of monosex cultures, is of particular interest to the aquaculture industry. Sex determination in the widely farmed Nile tilapia (Oreochromis niloticus) is complex, involving genomic regions on at least three chromosomes (chromosomes 1, 3 and 23) and interacting in certain cases with elevated early rearing temperature as well. Thus, sex ratios may vary substantially from 50%. RESULTS: This study focused on mapping sex-determining quantitative trait loci (QTL) in families with skewed sex ratios. These included four families that showed an excess of males (male ratio varied between 64% and 93%) when reared at standard temperature (28°C) and a fifth family in which an excess of males (96%) was observed when fry were reared at 36°C for ten days from first feeding. All the samples used in the current study were genotyped for two single-nucleotide polymorphisms (rs397507167 and rs397507165) located in the expected major sex-determining region in linkage group 1 (LG 1). The only misassigned individuals were phenotypic males with the expected female genotype, suggesting that those offspring had undergone sex-reversal with respect to the major sex-determining locus. We mapped SNPs identified from double digest Restriction-site Associated DNA (ddRAD) sequencing in these five families. Three genetic maps were constructed consisting of 641, 175 and 1,155 SNPs from the three largest families. QTL analyses provided evidence for a novel genome-wide significant QTL in LG 20. Evidence was also found for another sex-determining QTL in the fifth family, in the proximal region of LG 1. CONCLUSIONS: Overall, the results from this study suggest that these previously undetected QTLs are involved in sex determination in the Nile tilapia, causing sex reversal (masculinisation) with respect to the XX genotype at the major sex-determining locus in LG 1.

The effects of feeding ß-glucan to Pangasianodon hypophthalmus on immune gene expression and resistance to Edwardsiella ictaluri.

Pangasianodon hypophthalmus (striped catfish) is an important aquaculture species and intensification of farming has increased disease problems, particularly Edwardsiella ictaluri. The effects of feeding ß-glucans on immune gene expression and resistance to E. ictaluri in P. hypophthalmus were explored. Fish were fed 0.1% fungal-derived ß-glucan or 0.1% commercial yeast-derived ß-glucan or a basal control diet without glucan. After 14 days of feeding, the mRNA expression of immune genes (transferrin, C-reactive protein, precerebellin-like protein, Complement C3 and factor B, 2a MHC class II and interleukin-1 beta) in liver, kidney and spleen were determined. Following this fish from each of the three diet treatment groups were infected with E. ictaluri and further gene expression measured 24 h post-infection (h.p.i.), while the remaining fish were monitored over 2 weeks for mortalities. Cumulative percentage mortality at 14 days post-infection (d.p.i.) was less in ß-glucan fed fish compared to controls. There was no difference in gene expression between dietary groups after feeding for 14 days, but there was a clear difference between infected and uninfected fish at 24 h.p.i., and based on principal component analysis ß-glucans stimulated the overall expression of immune genes in the liver, kidney and spleen at 24 h.p.i.