Identification of genetic loci associated with higher resistance to pancreas disease (PD) in Atlantic salmon (Salmo salar L.).

BACKGROUND: Pancreas disease (PD) is a contagious disease caused by salmonid alphavirus (SAV) with significant economic and welfare impacts on salmon farming. Previous work has shown that higher resistance against PD has underlying additive genetic components and can potentially be improved through selective breeding. To better understand the genetic basis of PD resistance in Atlantic salmon, we challenged 4506 smolts from 296 families of the SalmoBreed strain. Fish were challenged through intraperitoneal injection with the most virulent form of the virus found in Norway (i.e., SAV3). Mortalities were recorded, and more than 900 fish were further genotyped on a 55 K SNP array. RESULTS: The estimated heritability for PD resistance was 0.41 ± 0.017. The genetic markers on two chromosomes, ssa03 and ssa07, showed significant associations with higher disease resistance. Collectively, markers on these two QTL regions explained about 60% of the additive genetic variance. We also sequenced and compared the cardiac transcriptomics of moribund fish and animals that survived the challenge with a focus on candidate genes within the chromosomal segments harbouring QTL. Approximately 200 genes, within the QTL regions, were found to be differentially expressed. Of particular interest, we identified various components of immunoglobulin-heavy-chain locus B (IGH-B) on ssa03 and immunoglobulin-light-chain on ssa07 with markedly higher levels of transcription in the resistant animals. These genes are closely linked to the most strongly QTL associated SNPs, making them likely candidates for further investigation. CONCLUSIONS: The findings presented here provide supporting evidence that breeding is an efficient tool for increasing PD resistance in Atlantic salmon populations. The estimated heritability is one of the largest reported for any disease resistance in this species, where the majority of the genetic variation is explained by two major QTL. The transcriptomic analysis has revealed the activation of essential components of the innate and the adaptive immune responses following infection with SAV3. Furthermore, the complementation of the genomic with the transcriptomic data has highlighted the possible critical role of the immunoglobulin loci in combating PD virus.

Genome-Wide Association Study of Piscine Myocarditis Virus (PMCV) Resistance in Atlantic Salmon (Salmo salar).

Cardiomyopathy syndrome is a severe, viral disease of Atlantic salmon that mostly affects farmed animals during their late production stage at sea. Caused by piscine myocarditis virus (PMCV), over the past few years outbreaks due to this disease have resulted in significant losses to the aquaculture industry. However, there is currently no vaccine that has proven effective against this virus. In this study, using a challenge model, we investigated the genetic variation for resistance to PMCV, by screening a large number of animals using a 55 K SNP array. In particular, we aimed to identify genetic markers that are tightly linked to higher disease resistance and can potentially be used in breeding programs. Using genomic information, we estimated a heritability of 0.51 ± 0.06, suggesting that resistance against this virus, to a great extent, is controlled by genetic factors. Through association analysis, we identified a significant quantitative trait locus (QTL) on chromosome 27, explaining approximately 57% of the total additive genetic variation. The region harboring this QTL contains various immune-related candidate genes, many of which have previously been shown to have a different expression profile between the naïve and infected animals. We also identified a suggestive association on chromosome 12, with the QTL linked markers located in 2 putatively immune-related genes. These results are of particular interest, as they can readily be implemented into breeding programs, can further assist in fine-mapping the causative mutations, and help in better understanding the biology of the disease and the immunological mechanisms underlying resistance against PMCV.

Lack of Truncated IFITM3 Transcripts in Cells Homozygous for the rs12252-C Variant That is Associated With Severe Influenza Infection.

Interferon-induced transmembrane 3 (IFITM3) is known to restrict the entry of a range of enveloped viruses. The single nucleotide polymorphism rs12252-C within IFITM3 has been shown to be associated with severe influenza A virus infection. It has been suggested that rs12252-C results in expression of a truncated IFITM3 protein lacking the first 21 amino acids. By performing high-throughput RNA sequencing on primary dendritic cells and peripheral blood mononuclear cells isolated from pandemic H1N1 influenza and human immunodeficiency virus-1 (HIV-1) infected patients we show that full-length IFITM3 mRNA is dominantly expressed (>99%) across all rs12252 genotypes. Full-length IFITM3 protein can be detected in all genotypes.

W chromosome expression responds to female-specific selection.

The W chromosome is predicted to be subject to strong female-specific selection stemming from its female-limited inheritance and therefore should play an important role in female fitness traits. However, the overall importance of directional selection in shaping the W chromosome is unknown because of the powerful degradative forces that act to decay the nonrecombining sections of the genome. Here we greatly expand the number of known W-linked genes and assess the expression of the W chromosome after >100 generations of different female-specific selection regimens in different breeds of chicken and in the wild ancestor, the Red Jungle Fowl. Our results indicate that female-specific selection has a significant effect on W chromosome gene-expression patterns, with a strong convergent pattern of up-regulation associated with increased female-specific selection. Many of the transcriptional changes in the female-selected breeds are the product of positive selection, suggesting that selection is an important force in shaping the evolution of gene expression on the W chromosome, a finding consistent with both the importance of the W chromosome in female fertility and the haploid nature of the W. Taken together, these data provide evidence for the importance of the sex-limited chromosome in a female heterogametic species and show that sex-specific selection can act to preserve sex-limited chromosomes from degrading forces.

Identification and Genomic Localization of Autosomal sdY Locus in a Population of Atlantic Salmon (Salmo salar).

The determination of sex in salmonid fishes is controlled by genetic mechanisms, with males being the heterogametic sex. The master sex-determining gene, the sexually dimorphic gene on the Y chromosome (sdY), is a conserved gene across various salmonid species. Nevertheless, variations in the genomic location of sdY have been observed both within and between species. Furthermore, different studies have reported discordances in the association between the sdY and the phenotypic gender. While some males seem to lack this locus, there have been reports of females carrying sdY. Although the exact reasons behind this discordance remain under investigation, some recent studies have proposed the existence of an autosomal, non-functional copy of sdY as a potential cause. In this study, we confirmed the presence of this autosomal sdY in the SalmoBreed strain of Atlantic salmon using a genotyping platform through a novel approach that allows for high-throughput screening of a large number of individuals. We further characterized the segregation profile of this locus across families and found the ratio of genetically assigned female-to-male progeny to be in accordance with the expected profile of a single autosomal sdY locus. Additionally, our mapping efforts localized this locus to chromosome 3 and suggested a putative copy on chromosome 6.

Characterization and 454 pyrosequencing of major histocompatibility complex class I genes in the great tit reveal complexity in a passerine system.

BACKGROUND: The critical role of Major Histocompatibility Complex (Mhc) genes in disease resistance and their highly polymorphic nature make them exceptional candidates for studies investigating genetic effects on survival, mate choice and conservation. Species that harbor many Mhc loci and high allelic diversity are particularly intriguing as they are potentially under strong selection and studies of such species provide valuable information as to the mechanisms maintaining Mhc diversity. However comprehensive genotyping of complex multilocus systems has been a major challenge to date with the result that little is known about the consequences of this complexity in terms of fitness effects and disease resistance. RESULTS: In this study, we genotyped the Mhc class I exon 3 of the great tit (Parus major) from two nest-box breeding populations near Oxford, UK that have been monitored for decades. Characterization of Mhc class I exon 3 was adopted and bidirectional sequencing was carried using the 454 sequencing platform. Full analysis of sequences through a stepwise variant validation procedure allowed reliable typing of more than 800 great tits based on 214,357 reads; from duplicates we estimated the repeatability of typing as 0.94. A total of 862 alleles were detected, and the presence of at least 16 functional loci was shown - the highest number characterized in a wild bird species. Finally, the functional alleles were grouped into 17 supertypes based on their antigen binding affinities. CONCLUSIONS: We found extreme complexity at the Mhc class I of the great tit both in terms of allelic diversity and gene number. The presence of many functional loci was shown, together with a pseudogene family and putatively non-functional alleles; there was clear evidence that functional alleles were under strong balancing selection. This study is the first step towards an in-depth analysis of this gene complex in this species, which will help understanding how parasite-mediated and sexual selection shape and maintain host genetic variation in nature. We believe that study systems like ours can make important contributions to the field of evolutionary biology and emphasize the necessity of integrating long-term field-based studies with detailed genetic analysis to unravel complex evolutionary processes.

Genetic architecture of body weight, condition factor and age of sexual maturation in Icelandic Arctic charr (Salvelinus alpinus).

The high commercial value from the aquaculture of salmonid fishes has prompted many studies into the genetic architecture of complex traits and the need to identify genomic regions that have repeatable associations with trait variation both within and among species. We searched for quantitative trait loci (QTL) for body weight (BW), condition factor (CF) and age of sexual maturation (MAT) in families of Arctic charr (Salvelinus alpinus) from an Icelandic breeding program. QTL with genome-wide significance were detected for each trait on multiple Arctic charr (AC) linkage groups (BW: AC-4, AC-20; CF: AC-7, AC-20, AC-23, AC-36; MAT: AC-13/34, AC-39). In addition to the genome-wide significant QTL for both BW and CF on AC-20, linkage groups AC-4, AC-7, AC-8, and AC-16 contain QTL for both BW and CF with chromosome-wide significance. These regions had effects (albeit weaker) on MAT with the exception of the region on AC-8. Comparisons with a North American cultured strain of Arctic charr, as well as North American populations of Atlantic salmon (Salmo salar), and rainbow trout (Oncorhynchus mykiss), reveal some conservation in QTL location and structure, particularly with respect to the joint associations of QTL influencing BW and CF. The detection of some differences in genetic architecture between the two aquaculture strains of Arctic charr may be reflective of the differential evolutionary histories experienced by these fishes, and illustrates the importance of including different strains to investigate genetic variation in a species where the intent is to use that variation in selective breeding programs.

Genome evolution in the fish family salmonidae: generation of a brook charr genetic map and comparisons among charrs (Arctic charr and brook charr) with rainbow trout.

BACKGROUND: Salmonids are regarded as 4R derivative species, having experienced 4 whole genome duplication events in their ancestry. Many duplicated chromosome regions still share extensive homology with one another which is maintained primarily through male-based homeologous chromosome pairings during meiosis. The formation of quadrivalents during meiosis leads to pseudolinkage. This phenomenon is more prevalent within 5 of the 12 ancestral teleost linkage groups in salmonids. RESULTS: We constructed a genetic linkage map for brook charr and used this in combination with the genetic map from Arctic charr, to make comparisons with the genetic map of rainbow trout. Although not all chromosome arms are currently mapped, some homologous chromosome rearrangements were evident between Arctic charr and brook charr. Notably, 10 chromosome arms in brook charr representing 5 metacentric chromosomes in Arctic charr have undergone rearrangements. Three metacentrics have one arm translocated and fused with another chromosome arm in brook charr to a make a new metacentrics while two metacentrics are represented by 4 acrocentric pairs in brook charr. In two cases (i.e., BC-4 and BC-16), an apparent polymorphism was observed with the identification of both a putative metacentric structure (similar to metacentric AC-4 = BC-4 and a joining of acrocentric AC-16 + one arm of AC-28 = BC-16), as well as two separate acrocentric linkage groups evident in the mapping parents. Forty-six of the expected 50 karyotypic arms could be inter-generically assigned. SEX in brook charr (BC-4) was localized to the same homologous linkage group region as in Arctic charr (AC-4). The homeologous affinities detected in the two charr species facilitated the identification of 20 (expected number = 25) shared syntenic regions with rainbow trout, although it is likely that some of these regions were partial or overlapping arm regions. CONCLUSIONS: Inter-generic comparisons among 2 species of charr (genus Salvelinus) and a trout (genus Oncorhynchus) have identified that linkage group arm arrangements are largely retained among these species. Previous studies have revealed that up to 7 regions of high duplicate marker retention occur between Salmo species (i.e., Atlantic salmon and brown trout) and rainbow trout, with 5 of these regions exhibiting higher levels of pseudolinkage. Pseudolinkage was detected in the charr species (i.e., BC-1/21, AC-12/27, AC-6/23, = RT-2p/29q, RT-12p/16p, and RT-27p/31p, respectively) consistent with three of the five 'salmonid-specific' pseudolinkage regions. Chromosome arms with the highest number of duplicated markers in rainbow trout are the linkage group arms with the highest retention of duplicated markers in both charr species.

Identification of a New Infectious Pancreatic Necrosis Virus (IPNV) Variant in Atlantic Salmon (Salmo salar L.) that can Cause High Mortality Even in Genetically Resistant Fish.

Infectious pancreatic necrosis (IPN) is an important viral disease of salmonids that can affect fish during various life cycles. In Atlantic salmon, selecting for genetically resistant fish against IPN has been one of the most highly praised success stories in the history of fish breeding. During the late 2000s, the findings that resistance against this disease has a significant genetic component, which is mainly controlled by variations in a single gene, have helped to reduce the IPN outbreaks to a great extent. In this paper, we present the identification of a new variant of the IPN virus from a field outbreak in Western Norway that had caused mortality, even in genetically resistant salmon. We recovered and assembled the full-length genome of this virus, following the deep-sequencing of the head-kidney transcriptome. The comparative sequence analysis revealed that for the critical amino acid motifs, previously found to be associated with the degree of virulence, the newly identified variant is similar to the virus's avirulent form. However, we detected a set of deduced amino acid residues, particularly in the hypervariable domain of the VP2, that collectively are unique to this variant compared to all other reference sequences assessed in this study. We suggest that these mutations have likely equipped the virus with the capacity to escape the host defence mechanism more efficiently, even in the genetically deemed IPN resistant fish.

Growth-related quantitative trait loci in domestic and wild rainbow trout (Oncorhynchus mykiss).

BACKGROUND: Somatic growth is a complex process that involves the action and interaction of genes and environment. A number of quantitative trait loci (QTL) previously identified for body weight and condition factor in rainbow trout (Oncorhynchus mykiss), and two other salmonid species, were used to further investigate the genetic architecture of growth-influencing genes in this species. Relationships among previously mapped candidate genes for growth and their co-localization to identified QTL regions are reported. Furthermore, using a comparative genomic analysis of syntenic rainbow trout linkage group clusters to their homologous regions within model teleost species such as zebrafish, stickleback and medaka, inferences were made regarding additional possible candidate genes underlying identified QTL regions. RESULTS: Body weight (BW) QTL were detected on the majority of rainbow trout linkage groups across 10 parents from 3 strains. However, only 10 linkage groups (i.e., RT-3, -6, -8, -9, -10, -12, -13, -22, -24, -27) possessed QTL regions with chromosome-wide or genome-wide effects across multiple parents. Fewer QTL for condition factor (K) were identified and only six instances of co-localization across families were detected (i.e. RT-9, -15, -16, -23, -27, -31 and RT-2/9 homeologs). Of note, both BW and K QTL co-localize on RT-9 and RT-27. The incidence of epistatic interaction across genomic regions within different female backgrounds was also examined, and although evidence for interaction effects within certain QTL regions were evident, these interactions were few in number and statistically weak. Of interest, however, was the fact that these predominantly occurred within K QTL regions. Currently mapped growth candidate genes are largely congruent with the identified QTL regions. More QTL were detected in male, compared to female parents, with the greatest number evident in an F1 male parent derived from an intercross between domesticated and wild strain of rainbow trout which differed strongly in growth rate. CONCLUSIONS: Strain background influences the degree to which QTL effects are evident for growth-related genes. The process of domestication (which primarily selects faster growing fish) may largely reduce the genetic influences on growth-specific phenotypic variation. Although heritabilities have been reported to be relatively high for both BW and K growth traits, the genetic architecture of K phenotypic variation appears less defined (i.e., fewer major contributing QTL regions were identified compared with BW QTL regions).

Genome-Wide Association Study Confirms Previous Findings of Major Loci Affecting Resistance to Piscine myocarditis virus in Atlantic Salmon (Salmo salar L.).

Cardiomyopathy syndrome is a viral disease of Atlantic salmon, mostly affecting fish during the late stages of production, resulting in significant losses to the industry. It has been shown that resistance to this disease has a strong genetic component, with quantitative trait loci (QTL) on chromosomes 27 (Ssa27) and Ssa12 to explain most of the additive genetic variance. Here, by analysing animals from a different year-class and a different population, we further aimed to confirm and narrow down the locations of these QTL. The data support the existence of the two QTL and suggest that the causative mutation on Ssa27 is most likely within the 10-10.5 Mbp segment of this chromosome. This region contains a cluster of major histocompatibility complex class I (MHC I) genes with the most strongly associated marker mapped to one of these loci. On Ssa12, the data confirmed the previous finding that the location of the causative mutation is within the 61.3 to 61.7 Mbp region. This segment contains several immune-related genes, but of particular interest are genes related to MHC II. Together, these findings highlight the likely key role of MHC genes in Atlantic salmon following infection with Piscine myocarditis virus (PMCV) and their potential impact on influencing the trajectory of this disease.

Comparative genomics and evolution of conserved noncoding elements (CNE) in rainbow trout.

BACKGROUND: Recent advances in the accumulation of genetic mapping and DNA sequence information from several salmonid species support the long standing view of an autopolyploid origin of these fishes (i.e., 4R). However, the paralogy relationships of the chromosomal segments descendent from earlier polyploidization events (i.e., 2R/3R) largely remain unknown, mainly due to an unbalanced pseudogenization of paralogous genes that were once resident on the ancient duplicated segments. Inter-specific conserved noncoding elements (CNE) might hold the key in identifying these regions, if they are associated with arrays of genes that have been highly conserved in syntenic blocks through evolution. To test this hypothesis, we investigated the chromosomal positions of subset of CNE in the rainbow trout genome using a comparative genomic framework. RESULTS: Through a genome wide analysis, we selected 41 pairs of adjacent CNE located on various chromosomes in zebrafish and obtained their intervening, less conserved, sequence information from rainbow trout. We identified 56 distinct fragments corresponding to about 150 Kbp of sequence data that were localized to 67 different chromosomal regions in the rainbow trout genome. The genomic positions of many duplicated CNE provided additional support for some previously suggested homeologies in this species. Additionally, we now propose 40 new potential paralogous affinities by analyzing the variation in the segregation patterns of some multi-copy CNE along with the synteny association comparison using several model vertebrates. Some of these regions appear to carry signatures of the 1R, 2R or 3R duplications. A subset of these CNE markers also demonstrated high utility in identifying homologous chromosomal segments in the genomes of Atlantic salmon and Arctic charr. CONCLUSION: CNE seem to be more efficacious than coding sequences in providing insights into the ancient paralogous affinities within the vertebrate genomes. Such a feature makes these elements extremely attractive for comparative genomics studies, as they can be treated as 'anchor' markers to investigate the association of distally located candidate genes on the homologous genomic segments of closely or distantly related organisms.

Genomic and Transcriptomic Analysis of Amoebic Gill Disease Resistance in Atlantic Salmon (Salmo salar L.).

Amoebic gill disease (AGD) is one of the most important parasitic diseases of farmed Atlantic salmon. It is a source of major economic loss to the industry and poses significant threats to animal welfare. Previous studies have shown that resistance against this disease has a moderate, heritable genetic component, although the genes and the genetic pathways that contribute to this process have yet to be elucidated. In this study, to identify the genetic mechanisms of AGD resistance, we first investigated the molecular signatures of AGD infection in Atlantic salmon through a challenge model, where we compared the transcriptome profiles of the naïve and infected animals. We then conducted a genome-wide association analysis with 1,333 challenged tested fish to map the AGD resistance genomic regions, supported by the results from the transcriptomic data. Further, we investigated the potential of incorporating gene expression analysis results in genomic prediction to improve prediction accuracy. Our data suggest thousands of genes have modified their expression following infection, with a significant increase in the transcription of genes with functional properties in cell adhesion and a sharp decline in the abundance of various components of the immune system genes. From the genome-wide association analysis, QTL regions on chromosomes ssa04, ssa09, and ssa13 were detected to be linked with AGD resistance. In particular, we found that QTL region on ssa04 harbors members of the cadherin gene family. These genes play a critical role in target recognition and cell adhesion. The QTL region on ssa09 also is associated with another member of the cadherin gene family, protocadherin Fat 4. The associated genetic markers on ssa13 span a large genomic region that includes interleukin-18-binding protein, a gene with function essential in inhibiting the proinflammatory effect of cytokine IL18. Incorporating gene expression information through a weighted genomic relationship matrix approach decreased genomic prediction accuracy and increased bias of prediction. Together, these findings help to improve our breeding programs and animal welfare against AGD and advance our knowledge of the genetic basis of host-pathogen interactions.

Distribution of ancestral proto-Actinopterygian chromosome arms within the genomes of 4R-derivative salmonid fishes (Rainbow trout and Atlantic salmon).

BACKGROUND: Comparative genomic studies suggest that the modern day assemblage of ray-finned fishes have descended from an ancestral grouping of fishes that possessed 12-13 linkage groups. All jawed vertebrates are postulated to have experienced two whole genome duplications (WGD) in their ancestry (2R duplication). Salmonids have experienced one additional WGD (4R duplication event) compared to most extant teleosts which underwent a further 3R WGD compared to other vertebrates. We describe the organization of the 4R chromosomal segments of the proto-ray-finned fish karyotype in Atlantic salmon and rainbow trout based upon their comparative syntenies with two model species of 3R ray-finned fishes. RESULTS: Evidence is presented for the retention of large whole-arm affinities between the ancestral linkage groups of the ray-finned fishes, and the 50 homeologous chromosomal segments in Atlantic salmon and rainbow trout. In the comparisons between the two salmonid species, there is also evidence for the retention of large whole-arm homeologous affinities that are associated with the retention of duplicated markers. Five of the 7 pairs of chromosomal arm regions expressing the highest level of duplicate gene expression in rainbow trout share homologous synteny to the 5 pairs of homeologs with the greatest duplicate gene expression in Atlantic salmon. These regions are derived from proto-Actinopterygian linkage groups B, C, E, J and K. CONCLUSION: Two chromosome arms in Danio rerio and Oryzias latipes (descendants of the 3R duplication) can, in most instances be related to at least 4 whole or partial chromosomal arms in the salmonid species. Multiple arm assignments in the two salmonid species do not clearly support a 13 proto-linkage group model, and suggest that a 12 proto-linkage group arrangement (i.e., a separate single chromosome duplication and ancestral fusion/fissions/recombination within the putative G/H/I groupings) may have occurred in the more basal soft-rayed fishes. We also found evidence supporting the model that ancestral linkage group M underwent a single chromosome duplication following the 3R duplication. In the salmonids, the M ancestral linkage groups are localized to 5 whole arm, and 3 partial arm regions (i.e., 6 whole arm regions expected). Thus, 3 distinct ancestral linkage groups are postulated to have existed in the G/H and M lineage chromosomes in the ancestor of the salmonids.

Determination of quantitative trait loci (QTL) for early maturation in rainbow trout (Oncorhynchus mykiss).

To identify quantitative trait loci (QTL) influencing early maturation (EM) in rainbow trout (Oncorhynchus mykiss), a genome scan was performed using 100 microsatellite loci across 29 linkage groups. Six inter-strain paternal half-sib families using three inter-strain F(1) brothers (approximately 50 progeny in each family) derived from two strains that differ in the propensity for EM were used in the study. Alleles derived from both parental sources were observed to contribute to the expression of EM in the progeny of the brothers. Four genome-wide significant QTL regions (i.e., RT-8, -17, -24, and -30) were observed. EM QTL detected on RT-8 and -24 demonstrated significant and suggestive QTL effects in both male and female progeny. Furthermore, within both male and female full-sib groupings, QTL on RT-8 and -24 were detected in two or more of the five parents used. Significant genome-wide and several strong chromosome-wide QTL for EM localized to different regions in males and females, suggesting some sex-specific control. Namely, QTL detected on RT-13, -15, -21, and -30 were associated with EM only in females, and those on RT-3, -17, and -19 were associated with EM only in males. Within the QTL regions identified, a comparison of syntenic EST markers from the rainbow trout linkage map with the zebrafish (Danio rerio) genome identified several putative candidate genes that may influence EM.

Impacts of Early Life Stress on the Methylome and Transcriptome of Atlantic Salmon.

Exposure to environmental stressors during early-life stages can change the rate and timing of various developmental processes. Epigenetic marks affecting transcriptional regulation can be altered by such environmental stimuli. To assess how stress might affect the methylome and transcriptome in salmon, fish were treated using cold-shock and air-exposure from the eye-stage until start-feeding. The fish were either stressed prior to hatching (E), post-hatching (PH), pre- and post-hatching (EPH) or not stressed (CO). Assessing transcriptional abundances just prior to start feeding, E and PH individuals were found to have modified the expression of thousands of genes, many with important functions in developmental processes. The EPH individuals however, showed expression similar to those of CO, suggesting an adaptive response to extended periods of stress. The methylome of stressed individuals differed from that of the CO, suggesting the importance of environment in shaping methylation signatures. Through integration of methylation with transcription, we identified bases with potential regulatory functions, some 10s of kb away from the targeted genes. We then followed fish growth for an additional year. Individuals in EPH showed superior growth compared to other treatment groups, highlighting how stress can potentially have long-lasting effects on an organism's ability to adapt to environmental perturbations.

The Atlantic Bonito (Sarda sarda, Bloch 1793) transcriptome and detection of differential expression during larvae development.

The Atlantic bonito (Sarda sarda, Bloch 1793) belongs to the important marine fish species with a wide geographical distribution covering the Atlantic Ocean, the Mediterranean and its bordering seas. Aquaculture practices for this species are still in their infancies and scientific studies are seldom undertaken, mainly because of difficulties in sampling. Thus for small tuna species like the Atlantic bonito only little is known about its biology and regarding the molecular background even less information is available. In the production of marine fish it is known that the most critical period is the larval stages, as high growth rates as well as significant developmental changes take place. In this study we have investigated the transcriptome of the Atlantic bonito of five larvae stages applying Illumina sequencing technology. For non-model species like aquaculture species, transcriptome analysis of RNA samples from individuals using Illumina sequencing technology is technically efficient and cost effective. In the present study a total number of 169,326,711 paired-end reads with a read length of 100 base pairs were generated resulting in a reference transcriptome of 68,220 contigs with an average length of 2054 base pairs. For differential expression analyses single end reads were obtained from different developmental stages and mapped to the constructed reference transcriptome. Differential expression analyses revealed in total 18,657 differentially expressed transcripts and were assigned to five distinguished groups. Each of the five clusters shows stage specific gene expression. We present for the first time in the Atlantic bonito an extensive RNA-Seq based characterization of its transcriptome as well as significant information on differential expression among five developmental larvae stages. The generated transcripts, including SNP and microsatellite information for candidate molecular markers and gene expression information will be a valuable resource for future genetic and molecular studies.

The plover neurotranscriptome assembly: transcriptomic analysis in an ecological model species without a reference genome.

We assembled a de novo transcriptome of short-read Illumina RNA-Seq data generated from telencephalon and diencephalon tissue samples from the Kentish plover, Charadrius alexandrinus. This is a species of considerable interest in behavioural ecology for its highly variable mating system and parental behaviour, but it lacks genomic resources and is evolutionarily distant from the few available avian draft genome sequences. We assembled and identified over 21,000 transcript contigs with significant expression in our samples, showing high homology to exonic sequences in avian draft genomes. From these, we identified >31,000 high-quality SNPs and > 2500 simple sequence repeats (SSRs). We also analysed expression patterns in our data to identify potential candidate genes related to differences in male and female behaviour, identifying over 200 nonoverlapping putative autosomal transcripts that show significant expression differences between males and females. Gene ontology analysis revealed that female-biased transcripts were significantly enriched for cerebral functions related to learning, cognition and memory, and male-biased transcripts were mostly enriched for terms related to neural function such as neuron projection and synapses. This data set provides one of the first de novo transcriptome assemblies from non-normalized short-read next-generation data and outlines an effective strategy for measuring sequence and expression variability simultaneously without the aid of a reference genome.

Trade-off between selection for dosage compensation and masculinization on the avian Z chromosome.

Following the suppression of recombination, gene expression levels decline on the sex-limited chromosome, and this can lead to selection for dosage compensation in the heterogametic sex to rebalance average expression from the X or Z chromosome with average autosomal expression. At the same time, due to their unequal pattern of inheritance in males and females, the sex chromosomes are subject to unbalanced sex-specific selection, which contributes to a nonrandom distribution of sex-biased genes compared to the remainder of the genome. These two forces act against each other, and the relative importance of each is currently unclear. The Gallus gallus Z chromosome provides a useful opportunity to study the importance and trade-offs between sex-specific selection and dosage compensation in shaping the evolution of the genome as it shows incomplete dosage compensation and is also present twice as often in males than females, and therefore predicted to be enriched for male-biased genes. Here, we refine our understanding of the evolution of the avian Z chromosome, and show that multiple strata formed across the chromosome over ∼130 million years. We then use this evolutionary history to examine the relative strength of selection for sex chromosome dosage compensation vs. the cumulative effects of masculinizing selection on gene expression. We find that male-biased expression increases over time, indicating that selection for dosage compensation is relatively less important than masculinizing selection in shaping Z chromosome gene expression.

Clock genes and their genomic distributions in three species of salmonid fishes: Associations with genes regulating sexual maturation and cell cycling.

BACKGROUND: Clock family genes encode transcription factors that regulate clock-controlled genes and thus regulate many physiological mechanisms/processes in a circadian fashion. Clock1 duplicates and copies of Clock3 and NPAS2-like genes were partially characterized (genomic sequencing) and mapped using family-based indels/SNPs in rainbow trout (RT)(Oncorhynchus mykiss), Arctic charr (AC)(Salvelinus alpinus), and Atlantic salmon (AS)(Salmo salar) mapping panels. RESULTS: Clock1 duplicates mapped to linkage groups RT-8/-24, AC-16/-13 and AS-2/-18. Clock3/NPAS2-like genes mapped to RT-9/-20, AC-20/-43, and AS-5. Most of these linkage group regions containing the Clock gene duplicates were derived from the most recent 4R whole genome duplication event specific to the salmonids. These linkage groups contain quantitative trait loci (QTL) for life history and growth traits (i.e., reproduction and cell cycling). Comparative synteny analyses with other model teleost species reveal a high degree of conservation for genes in these chromosomal regions suggesting that functionally related or co-regulated genes are clustered in syntenic blocks. For example, anti-müllerian hormone (amh), regulating sexual maturation, and ornithine decarboxylase antizymes (oaz1 and oaz2), regulating cell cycling, are contained within these syntenic blocks. CONCLUSIONS: Synteny analyses indicate that regions homologous to major life-history QTL regions in salmonids contain many candidate genes that are likely to influence reproduction and cell cycling. The order of these genes is highly conserved across the vertebrate species examined, and as such, these genes may make up a functional cluster of genes that are likely co-regulated. CLOCK, as a transcription factor, is found within this block and therefore has the potential to cis-regulate the processes influenced by these genes. Additionally, clock-controlled genes (CCGs) are located in other life-history QTL regions within salmonids suggesting that at least in part, trans-regulation of these QTL regions may also occur via Clock expression.