Identification of SNPs associated with muscle yield and quality traits using allelic-imbalance analyses of pooled RNA-Seq samples in rainbow trout

Application of genomic tools to study and potentially improve the upper thermal tolerance of farmed Atlantic salmon (Salmo salar)

BACKGROUND

The Atlantic salmon (Salmo salar) aquaculture industry must mitigate the impacts of rising ocean temperatures and the increased prevalence/severity of marine heat waves. Therefore, we investigated the genetic architecture and gene expression (transcriptomics) responsible for determining a salmon's upper thermal tolerance.

RESULTS

A genome-wide association study (GWAS) was conducted using fin clips of salmon from a previous incremental thermal maximum (ITMax) challenge (n = 251) and the North American 50 K SNP chip. ITMax was a highly polygenic trait with low/moderate heritability (mean SNP-based h2 = 0.20 and pedigree-based h2 = 0.25). Using data from the same fish, a separate GWAS assessed thermal-unit growth coefficient (TGC). Five significant SNPs were detected on chromosomes three and five, and high heritability estimates were calculated for TGC measured as fish grew from 12 to 20 °C (mean SNP-based h2 = 0.62 and pedigree-based h2 = 0.64). RNA-seq analyses of liver samples (n = 5-6 family-1 temperature-1) collected from the four most and four least tolerant families at 10 and 20 °C were also used to provide insights into potential mechanisms modulating this species' thermal tolerance. Between the top and bottom families, 347 and 175 differentially expressed transcripts (FDR-adjusted p < 0.01; fold-change ≥|2.0|) were identified at 10 and 20 °C, respectively. GO term enrichment analysis revealed unique responses to elevated temperature between family rankings (e.g., 'blood coagulation', 'sterol metabolic process' and 'synaptic growth at neuromuscular junction'). qPCR analyses further confirmed differences pertaining to cholesterol metabolism (lpl), inflammation (epx, elf3, ccl20), apoptosis (htra1b, htra2, anxa5b), angiogenesis (angl4, pdgfa), nervous system processes (insyn2a, kcnj11l) and heat stress (serpinh1b-1, serpinh1b-2). Three differentially expressed transcripts (i.e., ppp1r9a, gal3st1a, f5) were located in close proximity (± 120 kbp) to near-significant SNPs from the GWAS. Interestingly, ppp1r9a and gal3st1a have putative neurological functions, while f5 regulates blood coagulation.

CONCLUSIONS

These analyses provide several putative biomarkers of upper thermal tolerance in salmon that could prove valuable in helping the industry develop more temperature-tolerant fish. Further, our study supports previous reports that ITMax has low/moderate heritability in this species, and suggests that TGC at elevated temperatures is highly heritable.

Identification of Single-Nucleotide Polymorphisms in Differentially Expressed Genes Favoring Soybean Meal Tolerance in Higher-Growth Zebrafish (Danio rerio)

Genetic variability within the same fish species could confer soybean meal (SBM) tolerance in some individuals, thus favoring growth. This study investigates the single-nucleotide polymorphisms (SNPs) in differentially expressed genes (DEGs) favoring SBM tolerance in higher-growth zebrafish (Danio rerio). In a previous work, nineteen families of zebrafish were fed a fish meal diet (100FM control diet) or SBM-based diets supplemented with saponin (50SBM + 2SPN-experimental diet), from juvenile to adult stages. Individuals were selected from families with a genotype-by-environment interaction higher (170 ± 18 mg) or lower (76 ± 10 mg) weight gain on 50SBM + 2SPN in relation to 100FM. Intestinal transcriptomic analysis using RNA-seq revealed six hundred and sixty-five differentially expressed genes in higher-growth fish fed 50SBM + 2SPN diet. In this work, using these results, 47 SNPs in DEGs were selected. These SNPs were genotyped by Sequenom in 340 zebrafish that were fed with a 50SBM + 2SPN diet or with 100FM diet. Marker-trait analysis revealed 4 SNPs associated with growth in 3 immunity-related genes (aif1l, arid3c, and cst14b.2) in response to the 50SBM + 2SPN diet (p-value < 0.05). Two SNPs belonging to aif1l y arid3c produce a positive (+19 mg) and negative (-26 mg) effect on fish growth, respectively. These SNPs can be used as markers to improve the early selection of tolerant fish to SBM diet or other plant-based diets. These genes can be used as biomarkers to identify SNPs in commercial fish, thus contributing to the aquaculture sustainability.

Development and application of a 1K functional liquid chip for lactation performance in Bactrian camels

Introduction

The advancement of high-throughput, high-quality, flexible, and cost-effective genotyping platforms is crucial for the progress of dairy breeding in Bactrian camels. This study focuses on developing and evaluating a 1K functional liquid single nucleotide polymorphism (SNP) array specifically designed for milk performance in Bactrian camels.

Methods

We utilized RNA sequencing data from 125 lactating camels to identify and select 1,002 loci associated with milk production traits for inclusion in the SNP array. The array's performance was then assessed using 24 randomly selected camels. Additionally, the array was employed to genotype 398 individuals, which allowed for population validation to assess the polymorphism of SNP sites.

Results

The SNP array demonstrated high overall SNP call rates (> 99%) and a remarkable 100% consistency in genotyping. Population validation results indicate that camels from six breeding areas in Northwest China share a similar genetic background regarding lactation functionality.

Discussion

This study highlights the potential of the SNP array to accelerate the breeding process of lactating Bactrian camels and provides a robust technical foundation for improving lactation performance.

FishSNP: a high quality cross-species SNP database of fishes

The progress of aquaculture heavily depends on the efficient utilization of diverse genetic resources to enhance production efficiency and maximize profitability. Single nucleotide polymorphisms (SNPs) have been widely used in the study of aquaculture genomics, genetics, and breeding research since they are the most prevalent molecular markers on the genome. Currently, a large number of SNP markers from cultured fish species are scattered in individual studies, making querying complicated and data reuse problematic. We compiled relevant SNP data from literature and public databases to create a fish SNP database, FishSNP ( http://bioinfo.ihb.ac.cn/fishsnp ), and also used a unified analysis pipeline to process raw data that the author of the literature did not perform SNP calling on to obtain SNPs with high reliability. This database presently contains 45,690,243 (45 million) nonredundant SNP data for 13 fish species, with 30,288,958 (30 million) of those being high-quality SNPs. The main function of FishSNP is to search, browse, annotate and download SNPs, which provide researchers various and comprehensive associated information.

Functional SNPs in SYISL promoter significantly affect muscle fiber density and muscle traits in pigs

Our previous studies showed that SYISL is a negative regulator of muscle growth and regeneration in mice, pigs and humans. SYISL knockout resulted in an increase in the density of muscle fibers and muscle growth. However, it is unclear whether there are natural mutations in pig SYNPO2 intron sense-overlapping lncRNA (pSYISL) that affect the expression of pSYISL and muscle growth traits. In this study, three SNPs in exons and six SNPs within the promoter of pSYISL were identified. Association analysis showed that the two SNPs in exons are significantly associated with loin muscle area (p < 0.05); the six SNPs in the promoter that show complete linkage are significantly associated with live backfat thickness and live loin muscle area in American Large White pigs. Bioinformatics and luciferase reporter assays as well as in vitro binding experiments indicated that the mutation of SNP rs702045770 (g.539G>A) leads to the loss of YY1 binding to the promoter, thus affecting the expression level of pSYISL, and we found that Jiangshan Black pigs with genotype GG have a higher expression level of pSYISL than genotype AA individuals, but the muscle fiber density was significantly lower than in genotype AA individuals. Furthermore, the association analysis showed that the carcass backfat thickness of genotype GG of SNP rs702045770 was significantly higher than that of other genotypes in (Pietrain × Duroc) × (Landrace × Yorkshire) crossbred pigs (p < 0.05). The glycolytic potential of genotype GG was significantly higher than that of other genotypes (p < 0.05). These results provide novel insight into the identification of functional SNPs in non-coding genomic regions.

Association of Genetic Polymorphisms of TERT with Telomere Length in Coke Oven Emissions-Exposed Workers

We explored the association between variations in the telomere maintenance genes and change in telomere length (TL) in workers. The TL of peripheral blood leukocytes from 544 coke oven workers and 238 controls were detected using the Real-time PCR method. Variations in four genes were then detected using the PCR based restriction fragment length polymorphism. The effects of environmental and genetic factors on TL were subsequently analyzed through covariance analysis and a generalized linear model .The TL of subjects with GG genotypes were longer than those with AG genotype in the TERT rs2736098 locus amongst the controls (P = .032). The combined effect of COEs exposure and AG+AA genotypes had a significant effect on TL (P < .001). The interaction between the COEs exposure factor and the rs2736098AG+AA genotypes had a significant effect on the TL (P < .05). The TL in coke oven workers is associated with the interactions between TERT rs2736098 AG+AA and COEs exposure.

Identification and Functional Analysis of Transcriptome Profiles, Long Non-Coding RNAs, Single-Nucleotide Polymorphisms, and Alternative Splicing from the Oocyte to the Preimplantation Stage of Sheep by Single-Cell RNA Sequencing

Numerous dynamic and complicated processes characterize development from the oocyte to the embryo. However, given the importance of functional transcriptome profiles, long non-coding RNAs, single-nucleotide polymorphisms, and alternative splicing during embryonic development, the effect that these features have on the blastomeres of 2-, 4-, 8-, 16-cell, and morula stages of development has not been studied. Here, we carried out experiments to identify and functionally analyze the transcriptome profiles, long non-coding RNAs, single-nucleotide polymorphisms (SNPs), and alternative splicing (AS) of cells from sheep from the oocyte to the blastocyst developmental stages. We found between the oocyte and zygote groups significantly down-regulated genes and the second-largest change in gene expression occurred between the 8- and 16-cell stages. We used various methods to construct a profile to characterize cellular and molecular features and systematically analyze the related GO and KEGG profile of cells of all stages from the oocyte to the blastocyst. This large-scale, single-cell atlas provides key cellular information and will likely assist clinical studies in improving preimplantation genetic diagnosis.

Weighted Single-Step GWAS Identifies Genes Influencing Fillet Color in Rainbow Trout

The visual appearance of the fish fillet is a significant determinant of consumers' purchase decisions. Depending on the rainbow trout diet, a uniform bright white or reddish/pink fillet color is desirable. Factors affecting fillet color are complex, ranging from the ability of live fish to accumulate carotenoids in the muscle to preharvest environmental conditions, early postmortem muscle metabolism, and storage conditions. Identifying genetic markers of fillet color is a desirable goal but a challenging task for the aquaculture industry. This study used weighted, single-step GWAS to explore the genetic basis of fillet color variation in rainbow trout. We identified several SNP windows explaining up to 3.5%, 2.5%, and 1.6% of the additive genetic variance for fillet redness, yellowness, and whiteness, respectively. SNPs are located within genes implicated in carotenoid metabolism (β,β-carotene 15,15'-dioxygenase, retinol dehydrogenase) and myoglobin homeostasis (ATP synthase subunit β, mitochondrial (ATP5F1B)). These genes are involved in processes that influence muscle pigmentation and postmortem flesh coloration. Other identified genes are involved in the maintenance of muscle structural integrity (kelch protein 41b (klh41b), collagen α-1(XXVIII) chain (COL28A1), and cathepsin K (CTSK)) and protection against lipid oxidation (peroxiredoxin, superoxide dismutase 2 (SOD2), sestrin-1, Ubiquitin carboxyl-terminal hydrolase-10 (USP10)). A-to-G single-nucleotide polymorphism in β,β-carotene 15,15'-dioxygenase, and USP10 result in isoleucine-to-valine and proline-to-leucine non-synonymous amino acid substitutions, respectively. Our observation confirms that fillet color is a complex trait regulated by many genes involved in carotenoid metabolism, myoglobin homeostasis, protection against lipid oxidation, and maintenance of muscle structural integrity. The significant SNPs identified in this study could be prioritized via genomic selection in breeding programs to improve fillet color in rainbow trout.

Aquaculture omics: An update on the current status of research and data analysis

Aquaculture is the fast-growing agricultural sector and has the ability to meet the growing demand for protein nutritional security for future population. In future aquaculture is going to be the major source of fish proteins as capture fisheries reached at its maximum. However, several challenges need to overcome such as lack of genetically improved strains/varieties, lack of species-specific feed/functional feed, round the year availability of quality fish seed, pollution of ecosystems and increased frequencies of disease occurrence etc. In recent years, the continuous development of high throughput sequencing technology has revolutionized the biological sciences and provided necessary tools. Application of 'omics' in aquaculture research have been successfully used to resolve several productive and reproductive issues and thus ensure its sustainability and profitability. To date, high quality draft genomes of over fifty fish species have been generated and successfully used to develop large number of single nucleotide polymorphism markers (SNPs), marker panels and other genomic resources etc in several aquaculture species. Similarly, transcriptome profiling and miRNAs analysis have been used in aquaculture research to identify key transcripts and expression analysis of candidate genes/miRNAs involved in reproduction, immunity, growth, development, stress toxicology and disease. Metagenome analysis emerged as a promising scientific tool to analyze the complex genomes contained within microbial communities. Metagenomics has been successfully used in the aquaculture sector to identify novel and potential pathogens, antibiotic resistance genes, microbial roles in microcosms, microbial communities forming biofloc, probiotics etc. In the current review, we discussed application of high-throughput technologies (NGS) in the aquaculture sector.

Genetic Parameters and Genome-Wide Association Studies of Quality Traits Characterised Using Imaging Technologies in Rainbow Trout, Oncorhynchus mykiss

One of the top priorities of the aquaculture industry is the genetic improvement of economically important traits in fish, such as those related to processing and quality. However, the accuracy of genetic evaluations has been hindered by a lack of data on such traits from a sufficiently large population of animals. The objectives of this study were thus threefold: (i) to estimate genetic parameters of growth-, yield-, and quality-related traits in rainbow trout (Oncorhynchus mykiss) using three different phenotyping technologies [invasive and non-invasive: microwave-based, digital image analysis, and magnetic resonance imaging (MRI)], (ii) to detect quantitative trait loci (QTLs) associated with these traits, and (iii) to identify candidate genes present within these QTL regions. Our study collected data from 1,379 fish on growth, yield-related traits (body weight, condition coefficient, head yield, carcass yield, headless gutted carcass yield), and quality-related traits (total fat, percentage of fat in subcutaneous adipose tissue, percentage of fat in flesh, flesh colour); genotypic data were then obtained for all fish using the 57K SNP Axiom® Trout Genotyping array. Heritability estimates for most of the 14 traits examined were moderate to strong, varying from 0.12 to 0.67. Most traits were clearly polygenic, but our genome-wide association studies (GWASs) identified two genomic regions on chromosome 8 that explained up to 10% of the genetic variance (cumulative effects of two QTLs) for several traits (weight, condition coefficient, subcutaneous and total fat content, carcass and headless gutted carcass yields). For flesh colour traits, six QTLs explained 1-4% of the genetic variance. Within these regions, we identified several genes (htr1, gnpat, ephx1, bcmo1, and cyp2x) that have been implicated in adipogenesis or carotenoid metabolism, and thus represent good candidates for further functional validation. Finally, of the three techniques used for phenotyping, MRI demonstrated particular promise for measurements of fat content and distribution, while the digital image analysis-based approach was very useful in quantifying colour-related traits. This work provides new insights that may aid the development of commercial breeding programmes in rainbow trout, specifically with regard to the genetic improvement of yield and flesh-quality traits as well as the use of invasive and/or non-invasive technologies to predict such traits.

Genomic predictions for fillet yield and firmness in rainbow trout using reduced-density SNP panels

Background

One of the most important goals for the rainbow trout aquaculture industry is to improve fillet yield and fillet quality. Previously, we showed that a 50 K transcribed-SNP chip can be used to detect quantitative trait loci (QTL) associated with fillet yield and fillet firmness. In this study, data from 1568 fish genotyped for the 50 K transcribed-SNP chip and ~ 774 fish phenotyped for fillet yield and fillet firmness were used in a single-step genomic BLUP (ssGBLUP) model to compute the genomic estimated breeding values (GEBV). In addition, pedigree-based best linear unbiased prediction (PBLUP) was used to calculate traditional, family-based estimated breeding values (EBV).

Results

The genomic predictions outperformed the traditional EBV by 35% for fillet yield and 42% for fillet firmness. The predictive ability for fillet yield and fillet firmness was 0.19–0.20 with PBLUP, and 0.27 with ssGBLUP. Additionally, reducing SNP panel densities indicated that using 500–800 SNPs in genomic predictions still provides predictive abilities higher than PBLUP.

Conclusion

These results suggest that genomic evaluation is a feasible strategy to identify and select fish with superior genetic merit within rainbow trout families, even with low-density SNP panels.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-021-07404-9.

Distinct microbial assemblages associated with genetic selection for high- and low- muscle yield in rainbow trout

Background

Fish gut microbial assemblages play a crucial role in the growth rate, metabolism, and immunity of the host. We hypothesized that the gut microbiota of rainbow trout was correlated with breeding program based genetic selection for muscle yield. To test this hypothesis, fecal samples from 19 fish representing an F2 high-muscle genetic line (ARS-FY-H) and 20 fish representing an F1 low-muscle yield genetic line (ARS-FY-L) were chosen for microbiota profiling using the 16S rRNA gene. Significant differences in microbial assemblages between these two genetic lines might represent the effect of host genetic selection in structuring the gut microbiota of the host.

Results

Tukey’s transformed inverse Simpson indices indicated that high muscle yield genetic line (ARS-FY-H) samples have higher microbial diversity compared to those of the low muscle yield genetic line (ARS-FY-L) (LMM, χ2(1) =14.11, p < 0.05). The fecal samples showed statistically distinct structure in microbial assemblages between the genetic lines (F_1,36 = 4.7, p < 0.05, R² = 11.9%). Functional profiling of bacterial operational taxonomic units predicted characteristic functional capabilities of the microbial communities in the high (ARS-FY-H) and low (ARS-FY-L) muscle yield genetic line samples.

Conclusion

The significant differences of the microbial assemblages between high (ARS-FY-H) and low (ARS-FY-L) muscle yield genetic lines indicate a possible effect of genetic selection on the microbial diversity of the host. The functional composition of taxa demonstrates a correlation between bacteria and improving the muscle accretion in the host, probably, by producing various metabolites and enzymes that might aid in digestion. Further research is required to elucidate the mechanisms involved in shaping the microbial community through host genetic selection.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-020-07204-7.

Genome-wide scan for common variants associated with intramuscular fat and moisture content in rainbow trout

BACKGROUND

Genetic improvement of fillet quality attributes is a priority of the aquaculture industry. Muscle composition impacts quality attributes such as flavor, appearance, texture, and juiciness. Fat and moisture make up about ~ 80% of the tissue weight. The genetic architecture underlying the fat and moisture content of the muscle is still to be fully explored in fish. A 50 K gene transcribed SNP chip was used for genotyping 789 fish with available phenotypic data for fat and moisture content. Genotyped fish were obtained from two consecutive generations produced in the National Center for Cool and Cold Water Aquaculture (NCCCWA) growth-selective breeding program. Estimates of SNP effects from weighted single-step GBLUP (WssGBLUP) were used to perform genome-wide association (GWA) analysis to identify quantitative trait loci (QTL) associated with the studied traits.

RESULTS

Using genomic sliding windows of 50 adjacent SNPs, 137 and 178 SNPs were identified as associated with fat and moisture content, respectively. Chromosomes 19 and 29 harbored the highest number of SNPs explaining at least 2% of the genetic variation in fat and moisture content. A total of 61 common SNPs on chromosomes 19 and 29 affected the aforementioned traits; this association suggests common mechanisms underlying intramuscular fat and moisture content. Additionally, based on single-marker GWA analyses, 8 and 24 SNPs were identified in association with fat and moisture content, respectively.

CONCLUSION

SNP-harboring genes were primarily involved in lipid metabolism, cytoskeleton remodeling, and protein turnover. This work provides putative SNP markers that could be prioritized and used for genomic selection in breeding programs.

Comparative Analysis of the Transcriptome and Distribution of Putative SNPs in Two Rainbow Trout (Oncorhynchus mykiss) Breeding Strains by Using Next-Generation Sequencing

Selective breeding can significantly improve the establishment of sustainable and profitable aquaculture fish farming. For rainbow trout (Oncorhynchus mykiss), one of the main aquaculture coldwater species in Europe, a variety of selected hatchery strains are commercially available. In this study, we investigated the genetic variation between the local Born strain, selected for survival, and the commercially available Silver Steelhead strain, selected for growth. We sequenced the transcriptome of six tissues (gills, head kidney, heart, liver, spleen, and white muscle) from eight healthy individuals per strain, using RNA-seq technology to identify strain-specific gene-expression patterns and single nucleotide polymorphisms (SNPs). In total, 1760 annotated genes were differentially expressed across all tissues. Pathway analysis assigned them to different gene networks. We also identified a set of SNPs, which are heterozygous for one of the two breeding strains: 1229 of which represent polymorphisms over all tissues and individuals. Our data indicate a strong genetic differentiation between Born and Silver Steelhead trout, despite the relatively short time of evolutionary separation of the two breeding strains. The results most likely reflect their specifically adapted genotypes and might contribute to the understanding of differences regarding their robustness toward high stress and pathogenic challenge described in former studies.

Whole genome re-sequencing reveals recent signatures of selection in three strains of farmed Nile tilapia (Oreochromis niloticus)

Nile tilapia belongs to the second most cultivated group of fish in the world, mainly because of its favorable characteristics for production. Genetic improvement programs and domestication process of Nile tilapia may have modified the genome through selective pressure, leaving signals that can be detected at the molecular level. In this work, signatures of selection were identified using genome-wide SNP data, by two haplotype-based (iHS and Rsb) and one F_ST based method. Whole-genome re-sequencing of 326 individuals from three strains (A, B and C) of farmed tilapia maintained in Brazil and Costa Rica was carried out using Illumina HiSeq 2500 technology. After applying conventional SNP-calling and quality-control filters, ~ 1.3 M high-quality SNPs were inferred and used as input for the iHS, Rsb and F_ST based methods. We detected several candidate genes putatively subjected to selection in each strain. A considerable number of these genes are associated with growth (e.g. NCAPG, KLF3, TBC1D1, TTN), early development (e.g. FGFR3, PFKFB3), and immunity traits (e.g. NLRC3, PIGR, MAP1S). These candidate genes represent putative genomic landmarks that could be associated to traits of biological and commercial interest in farmed Nile tilapia.

Genome-wide identification of loci associated with growth in rainbow trout

Background

Growth is a major economic production trait in aquaculture. Improvements in growth performance will reduce time and cost for fish to reach market size. However, genes underlying growth have not been fully explored in rainbow trout.

Results

A previously developed 50 K gene-transcribed SNP chip, containing ~ 21 K SNPs showing allelic imbalances potentially associated with important aquaculture production traits including body weight, muscle yield, was used for genotyping a total of 789 fish with available phenotypic data for bodyweight gain. Genotyped fish were obtained from two consecutive generations produced in the NCCCWA growth-selection breeding program. Weighted single-step GBLUP (WssGBLUP) was used to perform a genome-wide association (GWA) analysis to identify quantitative trait loci (QTL) associated with bodyweight gain. Using genomic sliding windows of 50 adjacent SNPs, 247 SNPs associated with bodyweight gain were identified. SNP-harboring genes were involved in cell growth, cell proliferation, cell cycle, lipid metabolism, proteolytic activities, chromatin modification, and developmental processes. Chromosome 14 harbored the highest number of SNPs (n = 50). An SNP window explaining the highest additive genetic variance for bodyweight gain (~ 6.4%) included a nonsynonymous SNP in a gene encoding inositol polyphosphate 5-phosphatase OCRL-1. Additionally, based on a single-marker GWA analysis, 33 SNPs were identified in association with bodyweight gain. The highest SNP explaining variation in bodyweight gain was identified in a gene coding for thrombospondin-1 (THBS1) (R² = 0.09).

Conclusion

The majority of SNP-harboring genes, including OCRL-1 and THBS1, were involved in developmental processes. Our results suggest that development-related genes are important determinants for growth and could be prioritized and used for genomic selection in breeding programs.

Re-thinking adaptive immunity in the beetles: Evolutionary and functional trajectories of lncRNAs

Unlike vertebrate animals, invertebrates lack lymphocytes and therefore have historically been believed not to develop immune memory. A few studies have reported evidence of immune priming in insects; however, these studies lack the molecular mechanism and proposed it might be different among taxa. Since lncRNAs are known to regulate the immune response, we identified 10,120 lncRNAs in Tribolium castaneum genome-wide followed by transcriptome analysis of primed and unprimed larvae of different infectious status. A shift in lncRNA expression between Btt primed larvae and other treatment groups provides evidence of immune memory response. A few "priming" lncRNAs (n = 9) were uniquely regulated in Btt primed larvae. Evidence suggests these lncRNAs are likely controlling immune priming in Tribolium by regulating expression of genes involved in proteasomal machinery, Notch system, zinc metabolism, and methyltransferase activity, which are necessary to modulate phagocytosis. Our results support a conserved immune priming mechanism in a macrophage-dependent manner.

Genome-Wide Association Study Identifies Genomic Loci Affecting Filet Firmness and Protein Content in Rainbow Trout

Filet quality traits determine consumer satisfaction and affect profitability of the aquaculture industry. Soft flesh is a criterion for fish filet downgrades, resulting in loss of value. Filet firmness is influenced by many factors, including rate of protein turnover. A 50K transcribed gene SNP chip was used to genotype 789 rainbow trout, from two consecutive generations, produced in the USDA/NCCCWA selective breeding program. Weighted single-step GBLUP (WssGBLUP) was used to perform genome-wide association (GWA) analyses to identify quantitative trait loci affecting filet firmness and protein content. Applying genomic sliding windows of 50 adjacent SNPs, 212 and 225 SNPs were associated with genetic variation in filet shear force and protein content, respectively. Four common SNPs in the ryanodine receptor 3 gene (RYR3) affected the aforementioned filet traits; this association suggests common mechanisms underlying filet shear force and protein content. Genes harboring SNPs were mostly involved in calcium homeostasis, proteolytic activities, transcriptional regulation, chromatin remodeling, and apoptotic processes. RYR3 harbored the highest number of SNPs (n = 32) affecting genetic variation in shear force (2.29%) and protein content (4.97%). Additionally, based on single-marker analysis, a SNP in RYR3 ranked at the top of all SNPs associated with variation in shear force. Our data suggest a role for RYR3 in muscle firmness that may be considered for genomic- and marker-assisted selection in breeding programs of rainbow trout.

Genome-wide association analysis for body weight identifies candidate genes related to development and metabolism in rainbow trout (Oncorhynchus mykiss)

Growth is one of the most important traits from both a physiological and economic perspective in aquaculture species. Thus, identifying the genomic regions and genes underpinning genetic variation for this trait is of particular interest in several fish species, including rainbow trout. In this work, we perform a genome-wide association study (GWAS) to identify the genomic regions associated with body weight at tagging (BWT) and at 18 months (BW18M) using a dense SNP panel (57 k) and 4596 genotyped rainbow trout from 105 full-sib families belonging to a Chilean breeding population. Analysis was performed by means of single-step GBLUP approach. Genetic variance explained by 20 adjacent SNP windows across the whole genome is reported. To further explore candidate genes, we focused on windows that explained the highest proportion of genetic variance in the top 10 chromosomes for each trait. The main window from the top 10 chromosomes was explored by BLAST using the first and last SNP position of each window to determine the target nucleotide sequence. As expected, the percentage of genetic variance explained by windows was relatively low, due to the polygenic nature of body weight. The most important genomic region for BWT and BW18M were located on chromosomes 15 and 24 and they explained 2.14% and 3.02% of the genetic variance for each trait, respectively. Candidate genes including several growth factors, genes involved in development of skeletal muscle and bone tissue and nutrient metabolism were identified within the associated regions for both traits BWT and BW18M. These results indicate that body weight is polygenic in nature in rainbow trout, with the most important loci explaining as much as 3% of the genetic variance for the trait. The genes identified here represent good candidates for further functional validation to uncover biological mechanisms underlying variation for growth in rainbow trout.

Genome-Wide Association Analysis With a 50K Transcribed Gene SNP-Chip Identifies QTL Affecting Muscle Yield in Rainbow Trout

Detection of coding/functional SNPs that change the biological function of a gene may lead to identification of putative causative alleles within QTL regions and discovery of genetic markers with large effects on phenotypes. This study has two-fold objectives, first to develop, and validate a 50K transcribed gene SNP-chip using RNA-Seq data. To achieve this objective, two bioinformatics pipelines, GATK and SAMtools, were used to identify ~21K transcribed SNPs with allelic imbalances associated with important aquaculture production traits including body weight, muscle yield, muscle fat content, shear force, and whiteness in addition to resistance/susceptibility to bacterial cold-water disease (BCWD). SNPs ere identified from pooled RNA-Seq data collected from ~620 fish, representing 98 families from growth- and 54 families from BCWD-selected lines with divergent phenotypes. In addition, ~29K transcribed SNPs without allelic-imbalances were strategically added to build a 50K Affymetrix SNP-chip. SNPs selected included two SNPs per gene from 14K genes and ~5K non-synonymous SNPs. The SNP-chip was used to genotype 1728 fish. The average SNP calling-rate for samples passing quality control (QC; 1,641 fish) was ≥ 98.5%. The second objective of this study was to test the feasibility of using the new SNP-chip in GWA (Genome-wide association) analysis to identify QTL explaining muscle yield variance. GWA study on 878 fish (representing 197 families from 2 consecutive generations) with muscle yield phenotypes and genotyped for 35K polymorphic markers (passing QC) identified several QTL regions explaining together up to 28.40% of the additive genetic variance for muscle yield in this rainbow trout population. The most significant QTLs were on chromosomes 14 and 16 with 12.71 and 10.49% of the genetic variance, respectively. Many of the annotated genes in the QTL regions were previously reported as important regulators of muscle development and cell signaling. No major QTLs were identified in a previous GWA study using a 57K genomic SNP chip on the same fish population. These results indicate improved detection power of the transcribed gene SNP-chip in the target trait and population, allowing identification of large-effect QTLs for important traits in rainbow trout.

Integrated analysis of lncRNA and mRNA expression in rainbow trout families showing variation in muscle growth and fillet quality traits

Muscle yield and quality traits are important for the aquaculture industry and consumers. Genetic selection for these traits is difficult because they are polygenic and result from multifactorial interactions. To study the genetic architecture of these traits, phenotypic characterization of whole body weight (WBW), muscle yield, fat content, shear force and whiteness were measured in ~500 fish representing 98 families from a growth-selected line. RNA-Seq was used to sequence the muscle transcriptome of different families exhibiting divergent phenotypes for each trait. We have identified 240 and 1,280 differentially expressed (DE) protein-coding genes and long noncoding RNAs (lncRNAs), respectively, in fish families exhibiting contrasting phenotypes. Expression of many DE lncRNAs (n = 229) was positively correlated with overlapping, neighboring or distantly located protein-coding genes (n = 1,030), resulting in 3,392 interactions. Three DE antisense lncRNAs were co-expressed with sense genes known to impact muscle quality traits. Forty-four DE lncRNAs had potential sponge functions to miRNAs that affect muscle quality traits. This study (1) defines muscle quality associated protein-coding and noncoding genes and (2) provides insight into non-coding RNAs involvement in regulating growth and fillet quality traits in rainbow trout.

Crosstalk among lncRNAs, microRNAs and mRNAs in the muscle 'degradome' of rainbow trout

In fish, protein-coding and noncoding genes involved in muscle atrophy are not fully characterized. In this study, we characterized coding and noncoding genes involved in gonadogenesis-associated muscle atrophy, and investigated the potential functional interplay between these genes. Using RNA-Seq, we compared expression pattern of mRNAs, long noncoding RNAs (lncRNAs) and microRNAs of atrophying skeletal muscle from gravid females and control skeletal muscle from age-matched sterile individuals. A total of 852 mRNAs, 1,160 lncRNAs and 28 microRNAs were differentially expressed (DE) between the two groups. Muscle atrophy appears to be mediated by many genes encoding ubiquitin-proteasome system, autophagy related proteases, lysosomal proteases and transcription factors. Transcripts encoding atrogin-1 and mir-29 showed exceptional high expression in atrophying muscle, suggesting an important role in bulk muscle proteolysis. DE genes were co-localized in the genome with strong expression correlation, and they exhibited extensive 'lncRNA-mRNA', 'lncRNA-microRNA', 'mRNA-microRNA' and 'lncRNA-protein' physical interactions. DE genes exhibiting potential functional interactions comprised the highly correlated 'lncRNA-mRNA-microRNA' gene network described as 'degradome'. This study pinpoints extensive coding and noncoding RNA interactions during muscle atrophy in fish, and provides valuable resources for future mechanistic studies.