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

Abstract

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.

1. Introduction

The aquaculture industry produces food fish to satisfy a growing US and worldwide demand. Rainbow trout is the most cultivated, cool, freshwater fish in the United States [1]. Aquaculture supplies protein with low saturated fat and cholesterol content and high omega-3 fatty acids [2,3]. Rainbow trout are reared to produce fillets, and high production efficiency is needed to meet the ever-increasing demand for quality products. A significant constraint is the lack of genetically improved fish strains with high fillet yields and good-quality fillets. The industry has worked to remedy the situation by introducing breeding programs to select the best animals as parents for the next generations. Many of these breeding programs are traditional, using phenotypic information from breeding candidates and their pedigree to make selection decisions [4]. However, traditional breeding programs can be time-consuming and inefficient, especially for lethal traits such as fillet yield and color that cannot be accurately measured on live fish [5].

Use of genomic information in breeding programs offers a faster and more accurate method of achieving genetic progress. Achieving this goal requires understanding the genetic architecture underlying the variability in these traits. Genome-wide association (GWA) studies can identify genome regions associated with desired traits. GWA studies take advantage of linkage disequilibrium between SNP markers and genetic loci controlling a trait of interest. GWA studies have been conducted in the rainbow trout breeding program at the National Center for Cool- and Cold-Water Aquaculture (NCCCWA) for growth [6], muscle yield [7], intramuscular fat [8], fillet firmness [9], and disease resistance [10]. Another essential trait that requires attention is the fillet color—targeted in this study.

Fillet color is an important quality trait, usually influencing consumers’ satisfaction and point-of-purchase decisions. There are two markets for rainbow trout fillets—one for a bright reddish/pink fillet and one for bright white fillets. Consumers usually reject or downgrade pale yellowish fillets. Factors affecting fillet color range from genetics to environmental factors and to harvest, handling, and storage conditions [11]. Postmortem fillet color stability also depends on the rate of myoglobin oxidation, which is influenced by oxidation of intramuscular lipids and mitochondrial activity [12]. Salmonids’ characteristic pink/bright reddish fillet color results from the deposition of naturally occurring carotenoids or synthetic pigments added to the diets [13]. Carotenoids supplied in the diet are transported through the intestinal wall, metabolized within the cells of the intestinal linings or in the liver. The unmetabolized portion is deposited in the muscle by binding to muscle α-actin [14,15]. Rainbow trout flesh will typically be less reddish or whitish when the diet is not supplemented with carotenoids, as salmonids cannot synthesize carotenoids de novo [16]. Atlantic salmon and rainbow trout fish fed an unpigmented diet yield fillets with higher L* (lightness) and lower a* (redness) and b* (yellowness) values in comparison to fish on astaxanthin-supplemented diets [17,18]. Brown et al. [19] reported a significant difference in the color retention indices (ECI, hue, and chroma) of fillets from rainbow trout fish that never received dietary astaxanthin compared to fillets from fish that received an astaxanthin-supplemented diet. Even when fed a non-pigmented diet, Crouse et al. [20] observed significant differences in fillet redness (a*) between different rainbow trout strains fed the same diet. Red/pink pigmented rainbow trout fillets are deemed more desirable and marketed at a higher price than white fillets [19,21,22], but some consumers, especially in the US, may prefer a whiter fish.

Astaxanthin, a carotenoid added to the salmonid feed to improve the reddish color of the fillet, is an expensive feed ingredient, accounting for up to 30% of the feed cost. Therefore, development of genetically improved rainbow trout strains that more efficiently incorporate carotenoids into the muscle will benefit the aquaculture industry by improving profitability and consumer satisfaction. When fed unpigmented diets, genetically improved strains will also use naturally occurring carotenoids in the feed ingredients.

Studies in humans [21,22], chicken [23], mice [24], and Atlantic salmon [25] identified the β-carotene 15,15′-oxygenase (BCO1) enzyme as responsible for variation in the ability to metabolize carotenoids. Additionally, a recent study on a rainbow trout line used for commercial production in France identified Bcmo1 (β,β-carotene 15,15-dioxygenase), dkk3a (dickkopf WNT signaling pathway inhibitor 3a), and bola3 (bolA family member 3) as possible genes whose functions regulate the color of rainbow trout fillets [26]. Sae-Lim et al. [27] used the multi-trait GWAS approach to account for the relationship between body weight and fillet color and identified BCO1 and ppa1b (inorganic pyrophosphatase) within QTL regions influencing fillet color in Atlantic salmon. Other studies identified ATP-binding cassette subfamily G member 2 (abcg2-1a) in Atlantic salmon [28], PyBCO-1 in Scallop [29], and BCO2 in Chinook salmon [30] as candidate genes for fillet color. However, there is still much to learn about the genetic architecture of fillet color before it can be incorporated into breeding programs through genomic selection.

This study aims to use GWA analysis to identify genomic regions associated with fillet color traits (redness, yellowness, lightness, and whiteness) in a population of rainbow trout developed at the NCCCWA that had undergone five generations of selection for growth rate. Fish were fed an unpigmented commercial fishmeal-based diet.

2. Materials and Methods

2.1. Fish Population and Phenotype Used for GWA in This Study

The rainbow trout fish population used in this study was from a growth-selected line from NCCCWA, as described by Leeds et al. [31]. Fish from the third (hatch-year 2010) and fourth (hatch-year 2012) generations belonging to 197 families were included in this study. The breeding, selection, feeding, rearing, and harvesting procedures are as described by Salem et al. [7]. The fish used in this study were fed an unpigmented commercial fishmeal-based diet (42% protein, 16% fat; Ziegler Bros Inc., Gardners, PA, USA) using automatic feeders (Arvotec, Huutokoski, Finland). Initially, young fish were fed at a daily rate of ∼2.5% of body weight (BW), gradually reduced to approximately 0.75% of BW.

Fillet color parameters, L*, a*, and b*, which represent lightness, redness, and yellowness, respectively, were obtained from the fresh fillet surface using the Minolta Chroma Meter CR-200 (Minolta, Model CR-300; Minolta Camera Co., Osaka, Japan). The parameters were recorded a day after harvest at three locations above the lateral line of the right-side fillet, as described by Al-Tobasei et al. [32]. In addition to the standard color parameters, L*, a*, and b*, the fillet whiteness index was calculated using the equation: Whiteness = 100 − [(100 − L*)² + a*² + b*²]^1/2 [33]. The data were obtained from 878 fish from 2 harvest years: 406 from hatch-year 2010 and 472 from hatch-year 2012.

2.2. Genotyping and Quality Control

The 878 fish were genotyped with the 50k transcribed SNP-chip developed and described before [7]. PREGSF90 [34] was used to perform quality control using the following criteria: call rate for SNP and samples > 0.90, MAF > 0.05, monomorphic = 1, and HWE < 0.15. In total, 32,868 SNPs passed the QC and were used for subsequent analysis.

2.3. Descriptive Statistics

The mean and standard deviation values of each fillet color phenotype were calculated. Heritability was estimated as the ratio of additive genetic variance to total phenotypic variance. Variance components were estimated using the restricted maximum likelihood method found in AIREML in BLUPF90 software [34] using the following linear mixed model:

y = Xb + Z₁a + Z₂w + e

where y is the vector of phenotypes, b is the vector of fixed effects (age, harvest group, and hatch-year), a is the vector of additive genetic effect, w is the vector of random family effect, and e is the residual effect. X, Z₁, and Z₂ are incidence matrices for the effects contained in b, a, and w, respectively.

2.4. Genome-Wide Association Analysis

The weighted single-step GBLUP (wssGBLUP) approach proposed by Wang et al. [35] was used to perform genome-wide association analysis using the BLUPF90 family programs [34]. This method allows the use of genotyped and ungenotyped animals while integrating phenotype, genotype, and pedigree information in a mixed model for single-trait analysis.

The four fillet color parameters (L*, a*, b*, and whiteness) were analyzed using the single-trait animal model in wssGBLUP according to the model below:

y = Xb + Z₁a + Z₂w + e

where y is the vector of phenotypes, b is the vector of fixed effects, a is the vector of additive genetic effect, w is the vector of random family effect, and e is the residual effect. X, Z₁, and Z₂ are incidence matrices for the effects contained in b, a, and w, respectively. The fixed effects used in this study are fish age, harvest group, and hatch-year. The assumptions are that a~N(0, H${\mathsf{\sigma }}_a^2$) and e~N(0, I${\mathsf{\sigma }}_e^2$), where ${\mathsf{\sigma }}_a^2$ and ${\mathsf{\sigma }}_e^2$ are the additive genetic variance and residual variance, respectively. The H is a blend of pedigree and SNP-derived matrix [36], while I denotes the Identity matrix. The inverse of H is used in the wssGBLUP mixed model analysis [37].

${\mathrm{H}}^{-1}=A^{-1}+\left[\begin{array}{cc}0& 0\\ 0& G^{-1}{A}_{22}^{-1}\end{array}\right]$

where A⁻¹ is the inverse of the pedigree relationship matrix for all animals, $A_{22}^{-1}$is the inverse of the pedigree relationship matrix of genotyped animals, and G⁻¹ is the inverse of the genomic relationship matrix. The random family effect is uncorrelated and only accounts for the fact that the animals within the same family were raised in a common environment, and the covariance structure is given by I${\mathsf{\sigma }}_w^2$, where I is an identity matrix and ${\mathsf{\sigma }}_w^2$ is the family variance.

AIREMLF90 was used to estimate variance components supplied to BLUPF90 to predict genomic estimated breeding values (GEBV). The inbreeding coefficient was calculated from the pedigree data of 1420 fish by RENUMF90 using the method of Meuwissen and Luo [38].

BLUPF90 was used to predict breeding values using a weighted genomic relationship matrix (G). The SNP marker effect and new weights were then computed with POSTGSF90 [34] using 50 adjacent SNP sliding windows. All SNPs were initially assumed to be equally weighted (i.e., given an equal weight of 1.0). The final SNP weights and SNP effects were estimated using the option “non-linear A”, which allows for stable SNP weights after some iterations. Non-linear prediction assumes prior non-normal distribution of the marker effect and that markers do not contribute equally to genetic variance [39]. The non-linear approach resulted in greater reliability in the genomic prediction breeding value for bulls [39].

The percentage of additive genetic variance explained by each SNP window was calculated as:

$\frac{var\left(a_i\right)}{{\mathsf{\sigma }}_a^2}\times 100\%=\frac{var\left({\mathsf{\Sigma }}_{j=i}^{50}{\mathrm{Z}}_j{\mu }_j\right)}{{\mathsf{\sigma }}_a^2}\times 100\%$

where a_i is the genetic value of the i-th window consisting of 50 adjacent SNPs, σ² is the total genetic variance, z_j is a vector genotype of the j-th SNP for all animals, and ${\mu }_j$ is the SNP effect of the j-th SNP within the i-th window.

The qqman package [40] was used to obtain Manhattan plots for the proportion of additive genetic variance explained by each SNP window.

2.5. Identification of Candidate Genes

Genomic windows explaining at least 1% of the genetic variance were selected as possible genetic regions associated with the fillet color traits. The 1% threshold was set based on the literature obtained [41,42,43]. The SNPs were annotated using the NCBI rainbow trout genome assembly (GCF_013265735.2) to identify SNP-harboring genes. We used a literature search to identify relevant gene pathways and functions to understand the possible mechanisms by which the candidate genes regulate the traits. Genes previously identified in the literature or found to be related to color traits were further discussed.

2.6. MicroRNA Target Prediction

SNPs located in the 3’UTR of genes associated with fillet color in this study were investigated if their 3’UTR served as a target site for rainbow trout microRNAs. MicroRNA targets were predicted using three algorithms (PITA, miRanda, and TargetSpy) from the sRNAToolbox (http://bioinfo5.ugr.es/srnatoolbox). The rainbow trout microRNA repertoire was obtained from Juanchich et al. [44].

3. Results

3.1. Descriptive Statistics and Heritability Estimates for the Color Traits

There was more variation in the redness (54%) and yellowness (30%) in comparison to lightness and whiteness (6%) (Table 1). The heritability estimates of the traits in this population were moderate (0.16–0.39). The phenotypic correlation between lightness (L*) and whiteness color indices was 0.99 (R² = 0.98), whiteness and redness (a*) was 0.28, whiteness and yellowness (b*) was 0.67, and redness (a*) and yellowness (b*) was 0.45 (Table 1).

3.2. Genome-Wide Association Study and QTL Identification

A weighted single-step GBLUP approach was implemented in the BLUPF90 family of programs [34] to identify SNPs associated with fillet color traits. The GWAS results for whiteness and lightness color indices follow the same pattern, as expected, because of the high phenotypic correlation. Subsequently, only the whiteness trait will be discussed further. We identified 244, 161, and 115 SNPs in genomic windows, explaining at least 1% of the genetic variation in fillet redness, yellowness, and whiteness, respectively (Table 2 and Table 3, and Supplementary Table S1). The SNPs were identified within a genomic sliding window of 50 SNPs.

For redness (a*), chromosome 7 harbors the majority (33%) of the SNPs (80), followed by chromosome 9 (67 SNPs) (Figure 1, Table 2). Forty-five percent of the SNPs are in untranslated regions of genes, forty-two percent are in the coding regions. The highest peak corresponds to a SNP window on chromosome 7 that explains ~3.5% of the genetic variance.

For the yellowness trait (b*), most of the SNPs (66) are resident in chromosome 6 (41%), followed by 46 SNPs on chromosome 4 (29%). The peak SNP window, resident on chromosome 6, explains up to ~2.5% of the genetic variance for this trait (Figure 2, Table 2). Forty percent of the SNPs are in untranslated regions (UTR), while forty-seven percent are in coding regions.

Lightness (L*) and whiteness are similar in their genetic architecture, with peak SNP on chromosome 8 explaining only 1.6% of the genetic variance for this trait (Figure 3, Table 3). Forty-three percent of the SNPs are found within gene-coding regions, while forty-five percent are located in untranslated regions.

3.3. MicroRNA Target Prediction

Our results revealed that the 3’UTR region of ANKH (ANKH inorganic pyrophosphate transport regulator), RETRIG1 (reticulophagy regulator 1), and HSPB1 (heat-shock protein, α-crystallin-related, 1) genes are target sites for the omy-mir-1388-3p, omy-mir-219-5p, and omy-miR-724-5p microRNAs, respectively. An A-to-T single-nucleotide substitution at the target site of omy-mir-1388-3p causes a loss of its miRNA target site. Likewise, a C-to-T transition at the 3’UTR of HSPB1 resulted in a loss of the target site for the omy-miR-724-5p miRNA. Single-nucleotide substitution at the target site of omy-mir-219-5p does not lead to a loss of the target site.

4. Discussion

Fillet color is an important quality trait in salmonids influencing consumers’ purchasing decisions. Therefore, the industry is interested in selecting rainbow trout with superior genetic merit in their ability to produce a bright red or white fillet. Understanding the trait’s genetic architecture is required to determine the best genetic improvement approach. In this study, a genome-wide association investigation identifies regions of the genome influencing variability in fillet color traits in rainbow trout.

4.1. Descriptive Statistics and Heritability Estimates for the Color Traits

There is more variation in redness (a*) and yellowness (b*) compared to whiteness. The estimated heritability for fish in this population is low to moderate, similar to an estimate of 0.27 obtained for a rainbow trout fillet color score by Gjerde and Schaeffer [45]. A heritability estimate of 0.30 was recorded for fillet redness by Haffray et al. [46]. Blay et al. [26] reported higher heritability estimates of 0.46, 0.45, and 0.28 for rainbow trout fillet lightness (L*), redness (a*), and yellowness (b*), respectively. Overall, these studies demonstrate the possibility of achieving genetic improvement for fillet color traits through selection.

4.2. Summary of wssGWAS for Fillet Color Traits

The SNP windows explaining the highest genetic variance are found on chromosomes 7, 6, and 8 for fillet redness, yellowness, and whiteness, respectively. The SNP-harboring genes were classified according to their function and relevance to fillet color into the following categories.

4.3. Genes Involved in Carotenoid Metabolism

β, betacarotene 15,15-dioxygenase and retinol dehydrogenase are involved in carotenoid metabolism [47,48,49]. Fish species such as Atlantic salmon and rainbow trout deposit carotenoids in their muscle that enhance the reddish coloration of the fillet [15] and variation in carotenoid metabolism is associated with β-carotene oxygenase-1 function. Similar to the findings of this study, β, betacarotene 15,15-dioxygenase was implicated in its association with rainbow trout fillet yellowness [26]. Helgeland et al. [25] identified β-carotene oxygenase-1 (BCO1) and its paralogue β-carotene oxygenase-1 (BCO1L) as two probable causal genes influencing flesh color in Atlantic salmon. They supported their findings with functional studies of mRNA and protein expression, which pointed to BCO1L as the most likely of the two genes to influence flesh color variation. Several studies have identified single-nucleotide polymorphism within BCO1 that was associated with breast meat color in chicken [25,48,49]. In the mollusk, Yesso scallop, GWAS, and gene expression studies were used to confirm that PyBCO (a homolog of BCO1 in fish) was responsible for carotenoid metabolism and subsequent muscle coloration [29].

β,β-carotene 15,15′-dioxygenase on chromosome 6 explains 2.2% of the phenotypic variance for yellowness, while retinol dehydrogenase-7 found on chromosome 7 explains 2.9% of the variation in the redness trait (Table 2). A-G SNP in β,β-carotene 15,15′-dioxygenase causes isoleucine-to-valine non-synonymous amino acid substitution. Transcriptome analysis identifies retinol dehydrogenase-12 as a candidate gene regulating body-color formation in ornamental shrimp [50]. Carotenoids can serve as exogenous antioxidants to prevent cell oxidative damage, and these pigments inhibit lipid peroxidation and hemoglobin oxidation in human erythrocytes [51].

4.4. Genes Involved in Myoglobin Homeostasis and Protection against Lipid Oxidation

ATP5F1B, methylmalonyl-CoA mutase, ABC11, calsequestrin, cytochrome b5 (CYB5), ubiquitin carboxyl-terminal hydrolase 10 (USP10), peroxiredoxin, superoxide dismutase 2 (SOD2), sestrin-1, myosin X, and protein PRRC2C are genes that were found to affect fillet color in the study (Table 2 and Table 3). They are known to play a role in either myoglobin homeostasis or regulation of lipid peroxidation.

ATP synthase subunit β, mitochondrial (ATP5F1B) on chromosome 7 is another identified gene explaining over 3.5% of the genetic variability of fillet redness in rainbow trout (Table 2). It generates ATP from ADP through the electron transport system of the respiratory chain in the mitochondria [52]. Myoglobin is a muscle protein that binds oxygen and is responsible for muscle coloration [53]. Myoglobin exists in three forms: deoxymyoglobin, oxymyoglobin, and metmyoglobin. Although with low concentration of heme in the muscle, studies in salmonids have indicated that flesh color is, to some extent, dependent on the status of myoglobin [54,55]. The oxymyoglobin form promotes bright-reddish coloration in beef and salmon fillets, while metmyoglobin promotes fillet lightness (L*) [53,54]. The mitochondria function influences conversion between the three myoglobin forms [56]. Mitochondrial function can remain in postmortem muscle, influencing the conversion between myoglobin forms and the meat’s color [56,57]. Ramanathan et al. [53] suggested that understanding factors that influence mitochondrial function is key to unraveling the regulation of beef color appearance. Similarly, the gene ATP5F1B may influence rainbow trout fillet color by regulating mitochondrial integrity and function. Methylmalonyl-CoA mutase regulates mitochondria function by catalyzing the isomerization of methylmalonyl-CoA to succinyl-CoA [58]. It explains 1.29% of the genetic variability for fillet whiteness (Table 3).

Bile salt export pump (ABCB11) on chromosome 7 explains ~2.9% of the genetic variation in fillet redness (Table 2). This gene participates in bile acid homeostasis in an ATP-dependent manner [59,60]. It affects lipid metabolism and oxidation by regulating biliary tract lipid acid secretion through its action on bile salts’ excretion [61,62,63]. The influence of lipid oxidation on myoglobin, and thus meat color, is essential in meat color research. Postmortem meat color stability is affected by the muscle’s lipid oxidation rate [12]. The lipid auto-oxidation process generates free radicals and secondary products such as aldehydes and ketones that accelerate myoglobin oxidation [12] and, consequently, meat color deterioration [64,65,66]. Lipid peroxidation in the bile may generate pro-inflammatory agents by converting free fatty acids into lipid peroxides and aldehydes [67,68]. Chen et al. [66] discovered that aldehydes, a lipid oxidation product, accelerate the rate of myoglobin oxidation and promote permeability of the mitochondrial membrane. This process inhibits electron-transport chain-mediated metmyoglobin reduction and could profoundly affect fillet color stability, as discussed above with the ATP synthase subunit β, mitochondrial (ATP5F1B) gene. Blay et al. [26] identified two genes, dkk3 and bola3, known to be involved in adipogenesis, as genes harboring regulatory regions associated with fillet color. This work supports a relationship between fillet color and intramuscular fat content.

Conversion between the three myoglobin forms is influenced by mitochondrial function [56]. Cytochrome B5, a metmyoglobin reductase, reduces ferric myoglobin (methemoglobin) to ferrous myoglobin within muscle mitochondria [69,70]. In this study, cytochrome b5 (CYB5), on chromosome 6, explained up to 2.3% of the genetic variability for fillet yellowness (Table 2). This gene may play a role in the interconversion of three myoglobin forms, thereby influencing fillet coloration. The cytochrome c oxidase subunit II gene was a differentially expressed gene between red and chocolate ornamental shrimp [50].

Various studies have implicated ubiquitination as one of the regulatory mechanisms that determine meat quality in pork [71], lamb [72,73], and broiler chicken [74]. Ubiquitin carboxyl-terminal hydrolase-10 (USP10) is a member of the deubiquitinating enzyme family known as deubiquitinases, which include ubiquitin C-terminal hydroxylase-1 (UCH-L1) [75]. USP10 and ubiquitin carboxyl-terminal hydrolase 47 (both on chromosome 6), respectively, explain 2.3% and ~2.0% of genetic variance associated with the yellowness phenotype in this study (Table 2). UCH-L1 was implicated as influencing meat quality traits in pigs [76] and sheep [73]. UCH-L1 reportedly regulates oxidative activity in skeletal muscle [77] and plays a role in myogenesis [78]. Polymorphism (A/G) in USP10 causes a non-synonymous change in the amino acid from proline to leucine.

The peroxiredoxin family is a group of proteins capable of detoxifying peroxides and protecting cells against oxidation [79]. Peroxiredoxin-6 (PRDX6), on chromosome 4, explains 2.1% of genetic variance in fillet yellowness in the present fish population (Table 2). Proteome analysis of beef longissimus muscle revealed that peroxiredoxin-1 accounted for up to 70% of variances in color traits (L* a* b*) of muscle [80], and Wu et al. [79] identified peroxiredoxin as a possible marker for beef color. Peroxiredoxin-6 enzyme protects oxymyoglobin from peroxide attacks, thereby improving postmortem color stability [81]. Activator protein (AP-1) transcription factor on chromosome 6, which explains 1.7% of the genetic variance in fillet yellowness (Table 2), has been identified as a regulator of oxidative stress [82,83]. It protects the cell against reactive oxygen species. Other studies have identified a relationship between peroxiredoxin and meat quality or color traits in beef [84,85,86,87] and chevon [88]. Activation of the AP-1 transcription factor induces the expression of many antioxidants, including peroxiredoxin and glutathione reductase [89,90]. It is possible that these genes (PRDX6 and AP-1) function in the homeostatic regulation of the myoglobin redox state, protecting oxymyoglobin against oxidation and thereby enhancing the reddish coloration of the fillet.

Superoxide dismutase 2 (SOD2) encodes for muscle antioxidant enzyme. This enzyme reduces the damage caused by superoxide anion radicals [91]. Nohl et al. [92] identified superoxide dismutase as one of the agents protecting the mitochondria against lipid peroxidation and damage. Lipid oxidation and mitochondrial damage inhibit metmyoglobin reduction, and this causes muscle color deterioration [66]. SOD2 on chromosome 8 explained 1.5% of the genetic variance in fillet whiteness (Table 3) in this study. A proteomics study on color stability in lamb identified SOD2 as one of the proteins protecting the muscle against postmortem discoloration [93]. Superoxide dismutase was also a possible predictor of meat color stability in cattle [86] and chicken [94].

Sestrin-1 (SESN-1) on chromosome 8 explains ~1.4% of the genetic variance in fillet whiteness (Table 3). SESN1 is known to confer resistance to oxidative stress through regenerating peroxiredoxins [95,96,97].

Hanan and Shaklai [98] reported a peroxidative interaction between myoglobin and myosin that regulates myoglobin homeostasis when attacked by a peroxide. In vitro oxidation of oxymyoglobin was significantly greater (p < 0.05) when in the presence of myosin compared to when myosin is absent in Tuna fish and Sardine [99]. Myosin X (MYO10) on chromosome 8 explains 1.3% of the genetic variance for fillet whiteness (Table 3). MYO10 encodes for a myosin protein belonging to the myosin superfamily [100]. Myosin X may play a role in determining fillet color through its effect on oxymyoglobin oxidation.

Protein PRRC2C explains 2% of the genetic variability for fillet redness in this study (Table 2). Protein PRRC2C is required to efficiently form stress granules [101]. It is involved in the aggregation, arrangement, and bonding of proteins and RNA molecules to form a stress granule [101]. Stress granules are critical for facilitating responses against oxidative and cellular stress [102,103,104].

4.5. Genes Involved in Maintenance of Muscle Structural Integrity

The kelch protein 41b (KLH41B), collagen α-1(XXVIII) chain (COL28A1), myocilin (MYOC), F-actin-methionine sulfoxide oxidase (MICAL2), and cathepsin K (CTSK) are genes that are found to affect fillet color in the study. They are known to be involved in the maintenance of muscle structural integrity. KLH4LB on chromosome 7 explains up to 3.4% of this study’s variance in the redness trait (Table 2). KLH4LB is involved in skeletal muscle cell differentiation, muscle fiber development, and sarcomere organization [105]. Functional studies of the role of the KLH41B gene in zebrafish revealed that its knockout resulted in myofibrillar disorganization and muscle weakness [106]. The relationship between structure and fillet color has been reported in the literature. Kiessling et al. [107] reported that fillets with higher L* (lightness) values were softer than those with low lightness values. Gagaoua et al. [108] identified protein biomarkers (α-actin and connectin) for beef color traits that are also structural proteins. The structural attributes of the muscle could influence the extent of light scattering for meat [109]. In their study on mice, Ramirez-Martinez et al. [110] showed that KLH41B maintains muscle function by preferentially helping stabilize nebulin, a protein needed to maintain muscle sarcomere integrity. They revealed that proteins involved in sarcomere organization and muscle contraction regulation were downregulated in KLH41B knockout mice. Loss of nebulin causes nemaline myopathy in humans, a condition associated with severe muscle weakness [111].

Collagen α-1(XXVIII) chain (COL28A1) harbors the SNP marker for muscle color in broiler chicken [112]. The same gene explained ~3.5% of the variance in fillet redness in this study (Table 2). Collagen is a connective tissue protein. The muscle extracellular matrix is mainly composed of collagen family proteins [113]. The relative amount and distribution of collagen fibers in the muscle can influence muscle quality [114].

Cathepsin K activity influenced skeletal muscle repair in mice [115]. Cathepsin K (CTSK) explains up to 2.8% of this study’s genetic variance in fillet redness (Table 2).

Myocilin (MYOC) encodes the protein myocilin, which is involved in regulating the actin cytoskeleton [116]. It explained 1.95% of the genetic variance in fillet yellowness in this study (Table 2).

F-actin-methionine sulfoxide oxidase (MICAL2) encodes methionine monooxygenase, which promotes depolymerization of F-actin by mediating the oxidation of residues on actin to form methionine-sulfoxide, resulting in actin filament disassembly and preventing repolymerization [117,118]. The gene is also involved in cytoskeleton organization [118]. It explains 1.9% of the genetic variance in fillet yellowness in this study (Table 2).

The cysteine-rich secretory protein LCCL domain-containing 2 (CRISPLD2) gene encodes a protein binding heparin and glycosaminoglycans and is involved in regulating the innate immune system [119]. It was downregulated as part of broiler chickens’ regulatory mechanisms for muscle pigmentation [120]. It explains 2.45% of the genetic variance in fillet yellowness in this study (Table 2).

4.6. SNP Variants Alter MicroRNA Binding Sites

MicroRNAs (miRNAs) are short non-coding RNAs between 20 and 24 nucleotides in length and can regulate gene expression post-transcriptionally by binding to the 3’UTR of its target mRNA [121,122]. This binding process can form the RNA-induced silencing complex (RISC) and subsequently repression of translation [123]. Mutation/polymorphism in miRNA and/or the target 3’UTR sequence have been associated with phenotypic variation in economically important traits. A G-to-A SNP substitution in the myostatin 3’UTR changes the miRNA target site and affects muscularity in sheep [124]. A C/G polymorphism in the precursor region of microRNA affects body weight, pelvis breadth, and chest depth in chickens [125,126].

The 3’UTR region of ANKH (ANKH inorganic pyrophosphate transport regulator), RETRIG1 (reticulophagy regulator 1), and HSPB1 (heat-shock protein, α-crystallin-related, 1) genes are target sites for omy-mir-1388-3p, omy-mir-219-5p, and omy-miR-724-5p microRNAs, respectively (Table 2 and Table 3). An A-to-T single nucleotide substitution at the target site of omy-mir-1388-3p causes a loss of its miRNA target site. Likewise, a C-to-T transition at the 3’UTR of HSPB1 resulted in a loss of the target site for the omy-miR-724-5p miRNA. Single-nucleotide substitution at the target site of omy-mir-219-5p does not lead to loss of the target site.

Heat-shock proteins are common effectors of the cellular stress response. Thermal, environmental, or oxidative stress can trigger the transcription of genes encoding heat-shock proteins [127,128,129]. Protection of oxymyoglobin against oxidative stress is required to preserve bright-reddish meat coloration [66]. HSPB1 encodes for a heat-shock protein that can protect against oxidative stress. Diet supplementation with antioxidant vitamins resulted in a significant drop in HSPB1 expression in athletes after an exercise period compared with athletes fed an un-supplemented diet [130]. Over-expression of HSPB1 has been shown to improve stress resistance (including oxidative stress) [130,131,132,133,134]. MicroRNA can repress the translation of its target mRNA. It is possible that the loss of the HSPB1 target site facilitates the translation of the gene and induces resistance against oxidative stress.

Another gene that exhibits a loss of the miRNA target site under single-nucleotide polymorphism is the ANKH (ANKH inorganic pyrophosphate transport regulator). The gene encodes a protein that controls the extracellular level of pyrophosphate [135]. Inorganic pyrophosphate also plays an active role in oxidative stress resistance in several organisms [136,137,138].

5. Conclusions

We used weighted single-step GWAS to identify genetic variants associated with variability in fillet color traits in rainbow trout. Our result confirms that fillet color is a complex trait with no major gene but many SNP variants contributing to its regulation. We established that regulatory genes are involved in maintaining muscle structural integrity, carotenoid metabolism, or protection against myoglobin and lipid oxidation. An isoleucine-to-valine non-synonymous amino acid substitution mutation in β,β-carotene 15,15′-dioxygenase explained 2.2% of the phenotypic variance for yellowness, while SNP variants in retinol dehydrogenase-7 explained 2.9% of the variance in the muscle redness.

Author Contributions

M.S., T.L. and B.K. conceived, designed, and performed the experiment. R.O.A., R.A.-T. and A.A. analyzed the data. R.O.A., A.A., R.A.-T. and M.S. wrote the paper. All authors have read and agreed to the published version of the manuscript.

Institutional Review Board Statement

The Institutional Animal Care and Use Committee of the United States Department of Agriculture, National Center for Cool- and Cold-Water Aquaculture (Leetown, WV) specifically reviewed and approved all husbandry practices used in this study (IACUC protocol #056).

Informed Consent Statement

Not applicable.

Data Availability Statement

All datasets generated for this study are included in the manuscript and/or the Supplementary Files. The genotypes (ped and map files) and phenotypes are available in our previous publication [9].

Conflicts of Interest

The authors declare no conflict of interest.

Abbreviations

GWA: genome-wide association; HWE: Hardy–Weinberg equilibrium; LD: linkage disequilibrium; MAF: minor allele frequency; NCCCWA: USDA National Center of Cool- and Cold-Water Aquaculture; QC: quality control; QTL: quantitative trait loci; SNP: single-nucleotide polymorphism; UTR: untranslated region; WssGBLUP: weighted single-step GBLUP; YC: year class, L, a, b, w.

Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Figures and Tables

Figure 1. Manhattan plot of percent of genetic variance explained by 50 adjacent SNP windows for fillet redness (a*).

Figure 2. Manhattan plot of percent of genetic variance explained by 50 adjacent SNP windows for fillet yellowness (b*).

Figure 3. Manhattan plot of percent of genetic variance explained by 50 adjacent SNP windows for the fillet whiteness.

Table 1

Descriptive statistics of the observed phenotypes.

Trait	N	Mean	SD	Min	Max	CV (%)	${\mathsf{\sigma }}_{\mathit{a}}^2$	${\mathsf{\sigma }}_{\mathit{w}}^2$	${\mathsf{\sigma }}_{\mathit{e}}^2$	h2 (SE)
Redness	878	1.98	1.06	−0.17	5.833	0.54	0.08	0.04	0.38	0.16 ± 0.06
Yellowness	878	4.41	1.31	−0.79	8.123	0.30	0.52	0.16	0.67	0.39 ± 0.07
Lightness	878	44.54	2.74	38.17	54.81	0.06	1.23	0.33	4.13	0.22 ± 0.07
Whiteness	878	44.3	2.64	38.11	54.22	0.06	1.13	0.31	3.92	0.21 ± 0.06

Where ${\mathsf{\sigma }}_a^2$, ${\mathsf{\sigma }}_w^2$, and ${\mathsf{\sigma }}_e^2$ are the additive genetic variance, family variance, and residual variance, respectively, and h² is the heritability estimate.

Table 2

Selected SNP markers within 50 SNPs’ genomic sliding windows, explaining at least 1% of the additive genetic variance for fillet redness and yellowness traits.

Redness
Chr	POS	%Var	Gene ID	Gene Annotation	Region/Effect
7	10,996,914	2.43	LOC110527401	Radixin	CDS/syn
7	11,138,396	2.49	LOC110527405	Calsequestrin-2	CDS/Syn
7	11,312,252	2.73	LOC110527407	Zinc finger protein Dzip1	CDS/syn
7	11,399,310	3.45	LOC110527414	Kelch protein 41b	CDS/syn
7	11,402,881	3.47	LOC110527413	Collagen α-1(XXVIII) chain	3’UTR
7	11,438,574	3.29	LOC100136600	ATP synthase subunit β, mitochondrial	CDS/syn
7	11,444,638	3.02	LOC110527417	Retinol dehydrogenase 7	CDS/syn
7	11,459,018	2.98	abcb11	Bile salt export pump	CDS/syn
7	11,477,215	2.88	LOC100136260	Cathepsin K	CDS/syn
9	52,063,734	2.27	LOC110532529	Tyrosine-protein phosphatase non-receptor type 1	3’UTR
9	52,106,708	2.26	LOC110532530	Ubiquitin-conjugating enzyme E2 variant 1	3’UTR
9	52,291,239	2.28	LOC110532539	Partner of Y14 and mago A	CDS/syn
12	53,800,425	1.1	hspb1	Heat-shock protein, α-crystallin-related-1	3’UTR/miRNA target
Yellowness
4	22,957,625	2.09	prdx6	Peroxiredoxin 6	CDS/syn
4	22,973,619	2.11	plpp6	Phospholipid phosphatase 6	5’UTR
4	23,074,540	2	LOC110521622	Protein PRRC2C	3’UTR
4	23,103,208	1.92	vamp4	Vesicle-associated membrane protein 4	3’UTR
4	23,115,313	1.95	LOC110521624	Myocilin	CDS/Syn
6	61,578,946	1.9	LOC110526379	F-actin-methionine Sulfoxide oxidase MICAL2	5’UTR
6	61,592,297	1.99	LOC110526380	Ubiquitin carboxyl-terminal hydrolase 47	CDS/syn
6	61,666,093	2.25	LOC110526946	β,β-carotene 15,15′-dioxygenase-l	CDS/Non-syn
6	61,805,211	2.11	LOC110526388	Nuclear factor of activated T-cells 5	CDS/syn
6	61,837,913	2.39	LOC110526389	Lysine-tRNA ligase	3’UTR
6	61,847,413	2.38	LOC110526390	60S ribosomal protein L13	CDS/syn
6	61,998,041	2.36	LOC110526393	Cytochrome b5	3’UTR
6	62,768,347	2.45	LOC110526402	Cysteine-rich Secretory protein LCCL domain-containing 2	3’UTR
6	62,812,905	2.32	LOC110526403	Ubiquitin carboxyl-terminal hydrolase 10	CDS/Non-syn
6	62,896,859	2.24	LOC110526405	AP-1 complex subunit γ-1	3’UTR
6	62,961,238	2.26	LOC110526408	Myotubularin-related Protein 10	3’UTR
6	63,056,828	2.39	LOC100136691	Cyclin B2	CDS/syn

Chr = chromosome, POS = SNP position %Var = % variance explained, Syn = synonymous amino acid substitution, Non-Syn = non-synonymous amino acid substitution. Color intensities (green, yellow, and red) reflect changes in additive genetic variance explained by the SNP genomic sliding window for the fillet trait. A color gradient indicates differences in additive genetic variance explained by windows containing the representative SNP marker (green is the highest and red is the lowest).

Table 3

Selected SNP markers within 50 SNPs’ genomic sliding windows, explaining at least 1% of the additive genetic variance for the fillet whiteness trait.

		Whiteness
Chr	POS	%Var	Gene ID	Gene Annotation	Region/Effect
8	34,097,292	1.17	LOC110529884	Peptidyl-prolyl cis-trans isomerase FKBP1B	CDS/Syn
8	34,136,112	1.29	mut	Methylmalonyl-CoA mutase	3’UTR
8	34,495,040	1.49	sod2	Superoxide dismutase 2	3’UTR
8	34,936,875	1.57	LOC110529892	cGMP-dependent protein kinase 1	3’UTR
8	36,538,411	1.42	LOC110529899	SAM and SH3 domain-containing protein 1	3’UTR
8	37,290,793	1.54	LOC110529911	Sialomucin core protein 24	3’UTR
8	37,412,186	1.38	LOC110529910	Sestrin-1	3’UTR
8	37,829,107	1.38	ostm1	Osteopetrosis-associated transmembrane protein 1	3’UTR
8	38,254,068	1.3	LOC110529920	Poly(U)-binding-splicing factor PUF60	CDS/Syn
8	39,295,098	1.37	ankh	ANKH inorganic pyrophosphate transport regulator	3UTR/miRNA target
8	40,954,559	1.3	myo10	Myosin X	3’UTR
8	40,978,990	1.36	znf622	Zinc finger protein 622	3’UTR
8	41,002,542	1.26	retreg1	Reticulophagy regulator 1	3’UTR/miRNA target
19	41,952,271	1.18	LOC110497982	Uncharacterized protein C15orf52	3’UTR
27	1,675,710	1.19	LOC110507317	Protein IWS1 homolog	3’UTR
27	3,976,684	1.18	LOC110507360	Serine/threonine-protein phosphatase 2A 65 kDa regulatory subunit A β isoform	CDS/Syn

Chr = chromosome, POS = SNP position, %Var = % variance explained, Syn = synonymous amino acid substitution. Color intensities (green, yellow, and red) reflect changes in additive genetic variance explained by the SNP genomic sliding window for the fillet trait. A color gradient indicates differences in additive genetic variance explained by windows containing the representative SNP marker (green is the highest and red is the lowest).

Supplementary Materials

The following supporting information can be downloaded at: https://www.mdpi.com/article/10.3390/genes13081331/s1, Table S1: SNPs and genes identified in this study explaining at least 1% of the genetic variation in fillet redness, yellowness, and whiteness.

References

1. Thorgaard, G.H.; Bailey, G.S.; Williams, D.; Buhler, D.R.; Kaattari, S.L.; Ristow, S.S.; Hansen, J.D.; Winton, J.R.; Bartholomew, J.L.; Nagler, J.J. et al. Status and opportunities for genomics research with rainbow trout. Comp. Biochem. Physiol. B Biochem. Mol. Biol.; 2002; 133, pp. 609-646. [DOI: https://dx.doi.org/10.1016/S1096-4959(02)00167-7]

2. Turchini, G.M.; Francis, D.S.; Keast, R.S.; Sinclair, A.J. Transforming salmonid aquaculture from a consumer to a producer of long chain omega-3 fatty acids. Food Chem.; 2011; 124, pp. 609-614. [DOI: https://dx.doi.org/10.1016/j.foodchem.2010.06.083]

3. Harlioglu, A.G. Fatty acid composition, fat soluble vitamins and cholesterol content of farmed rainbow trout (Oncorhynchus mykiss). Pak. J. Zool.; 2012; 44, pp. 1013-1019.

4. Zenger, K.; Khatkar, M.; Jerry, D.; Raadsma, H. The next wave in selective breeding: Implementing genomic selection in aquaculture. Proc. Proc. Assoc. Advmt. Anim. Breed. Genet.; 2017; 22, pp. 105-112.

5. Ibtisham, F.; Zhang, L.; Xiao, M.; An, L.; Ramzan, M.B.; Nawab, A.; Zhao, Y.; Li, G.; Xu, Y. Genomic selection and its application in animal breeding. Thai J. Vet. Med.; 2017; 47, 301.

6. Ali, A.; Al-Tobasei, R.; Lourenco, D.; Leeds, T.; Kenney, B.; Salem, M. Genome-wide identification of loci associated with growth in rainbow trout. BMC Genom.; 2020; 21, 209. [DOI: https://dx.doi.org/10.1186/s12864-020-6617-x]

7. Salem, M.; Al-Tobasei, R.; Ali, A.; Lourenco, D.; Gao, G.; Palti, Y.; Kenney, B.; Leeds, T.D. Genome-Wide Association Analysis With a 50K Transcribed Gene SNP-Chip Identifies QTL Affecting Muscle Yield in Rainbow Trout. Front. Genet.; 2018; 9, 387. [DOI: https://dx.doi.org/10.3389/fgene.2018.00387]

8. Ali, A.; Al-Tobasei, R.; Lourenco, D.; Leeds, T.; Kenney, B.; Salem, M. Genome-wide scan for common variants associated with intramuscular fat and moisture content in rainbow trout. BMC Genom.; 2020; 21, 529. [DOI: https://dx.doi.org/10.1186/s12864-020-06932-0]

9. Ali, A.; Al-Tobasei, R.; Lourenco, D.; Leeds, T.; Kenney, B.; Salem, M. Genome-Wide Association Study Identifies Genomic Loci Affecting Filet Firmness and Protein Content in Rainbow Trout. Front. Genet.; 2019; 10, 386. [DOI: https://dx.doi.org/10.3389/fgene.2019.00386]

10. Vallejo, R.L.; Silva, R.M.O.; Evenhuis, J.P.; Gao, G.; Liu, S.; Parsons, J.E.; Martin, K.E.; Wiens, G.D.; Lourenco, D.A.L.; Leeds, T.D. et al. Accurate genomic predictions for BCWD resistance in rainbow trout are achieved using low-density SNP panels: Evidence that long-range LD is a major contributing factor. J. Anim. Breed. Genet.; 2018; 135, pp. 263-274. [DOI: https://dx.doi.org/10.1111/jbg.12335]

11. Colihueque, N. Genetics of salmonid skin pigmentation: Clues and prospects for improving the external appearance of farmed salmonids. Rev. Fish Biol. Fish.; 2010; 20, pp. 71-86. [DOI: https://dx.doi.org/10.1007/s11160-009-9121-6]

12. Joseph, P.; Suman, S.P.; Li, S.; Beach, C.M.; Steinke, L.; Fontaine, M. Characterization of bison (Bison bison) myoglobin. Meat Sci.; 2010; 84, pp. 71-78. [DOI: https://dx.doi.org/10.1016/j.meatsci.2009.08.014]

13. Buttle, L.; Crampton, V.; Williams, P. The effect of feed pigment type on flesh pigment deposition and colour in farmed Atlantic salmon, Salmo salar L. Aquac. Res.; 2001; 32, pp. 103-111. [DOI: https://dx.doi.org/10.1046/j.1365-2109.2001.00536.x]

14. Castenmiller, J.J.; West, C.E. Bioavailability and bioconversion of carotenoids. Annu. Rev. Nutr.; 1998; 18, pp. 19-38. [DOI: https://dx.doi.org/10.1146/annurev.nutr.18.1.19]

15. Matthews, S.J.; Ross, N.W.; Lall, S.P.; Gill, T.A. Astaxanthin binding protein in Atlantic salmon. Comp. Biochem. Physiol. B Biochem. Mol. Biol.; 2006; 144, pp. 206-214. [DOI: https://dx.doi.org/10.1016/j.cbpb.2006.02.007]

16. Storebakken, T.; No, H.K. Pigmentation of rainbow trout. Aquaculture; 1992; 100, pp. 209-229. [DOI: https://dx.doi.org/10.1016/0044-8486(92)90372-R]

17. Thomas, A.C. Astaxanthin in Juvenile Farmed Chinook salmon (Oncorhynchus tshawytscha): Effective Dietary Levels for Flesh Pigmentation and Influence on Fatty Acid Profile during Cold Temperature Storage of Fillets. Ph.D. Thesis; University of British Columbia: Vancouver, BC, Canada, 1999.

18. Arredondo-Figueroa, J.L.; Mora, G.I.d.l.; Ponce-Palafox, J.T.; Barriga-Sosa, I.d.l.A.; Vernon-Carter, E.J. Color of raw, frozen, and smoked fillets of rainbow trout (Oncorhynchus mykiss) fed diets supplemented with astaxanthin and saponified red chilli (Capsicum annuum) extracts. J. Aquat. Food Prod. Technol.; 2007; 16, pp. 35-50. [DOI: https://dx.doi.org/10.1300/J030v16n01_04]

19. Brown, K.R.; Barnes, M.; Parker, T.; Fletcher, B. Retention of fillet coloration in rainbow trout after dietary astaxanthin cessation. Fish. Aquac. J.; 2016; 7, 163. [DOI: https://dx.doi.org/10.4172/2150-3508.1000163]

20. Crouse, C.C.; Davidson, J.W.; Good, C.M.; May, T.C.; Summerfelt, S.T.; Kenney, P.B.; Leeds, T.D.; Cleveland, B.M. Growth and fillet quality attributes of five genetic strains of rainbow trout (Oncorhynchus mykiss) reared in a partial water reuse system and harvested at different sizes. Aquac. Res.; 2018; 49, pp. 1672-1681. [DOI: https://dx.doi.org/10.1111/are.13623]

21. Yabuta, S.; Masaki, M.; Shidoji, Y. Associations of Buccal Cell Telomere Length with Daily Intake of beta-Carotene or alpha-Tocopherol Are Dependent on Carotenoid Metabolism-related Gene Polymorphisms in Healthy Japanese Adults. J. Nutr. Health Aging; 2016; 20, pp. 267-274. [DOI: https://dx.doi.org/10.1007/s12603-015-0577-x]

22. Lietz, G.; Oxley, A.; Leung, W.; Hesketh, J. Single nucleotide polymorphisms upstream from the beta-carotene 15,15′-monoxygenase gene influence provitamin A conversion efficiency in female volunteers. J. Nutr.; 2012; 142, pp. 161S-165S. [DOI: https://dx.doi.org/10.3945/jn.111.140756] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/22113863]

23. Le Bihan-Duval, E.; Nadaf, J.; Berri, C.; Pitel, F.; Graulet, B.; Godet, E.; Leroux, S.Y.; Demeure, O.; Lagarrigue, S.; Duby, C. Detection of a Cis eQTL controlling BMCO1 gene expression leads to the identification of a QTG for chicken breast meat color. PLoS ONE; 2011; 6, e14825.

24. Hessel, S.; Eichinger, A.; Isken, A.; Amengual, J.; Hunzelmann, S.; Hoeller, U.; Elste, V.; Hunziker, W.; Goralczyk, R.; Oberhauser, V. et al. CMO1 deficiency abolishes vitamin A production from beta-carotene and alters lipid metabolism in mice. J. Biol. Chem.; 2007; 282, pp. 33553-33561. [DOI: https://dx.doi.org/10.1074/jbc.M706763200]

25. Helgeland, H.; Sodeland, M.; Zoric, N.; Torgersen, J.S.; Grammes, F.; von Lintig, J.; Moen, T.; Kjoglum, S.; Lien, S.; Vage, D.I. Genomic and functional gene studies suggest a key role of beta-carotene oxygenase 1 like (bco1l) gene in salmon flesh color. Sci. Rep.; 2019; 9, 20061. [DOI: https://dx.doi.org/10.1038/s41598-019-56438-3] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/31882713]

26. Blay, C.; Haffray, P.; Bugeon, J.; D’ambrosio, J.; Dechamp, N.; Collewet, G.; Enez, F.; Petit, V.; Cousin, X.; Corraze, G. Genetic parameters and genome-wide association studies of quality traits characterised using imaging technologies in Rainbow trout, Oncorhynchus mykiss. Front. Genet.; 2021; 12, 219. [DOI: https://dx.doi.org/10.3389/fgene.2021.639223] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/33692832]

27. Sae-Lim, P.; Boison, S.; Baranski, M. Multi-Trait Genome Wide Association Study in Correlated Traits: Fillet Colour and Body Weight in Atlantic Salmon; Wageningen Academic Publishers: Wageningen, The Netherlands, Available online: https://www.wageningenacademic.com/pb-assets/wagen/WCGALP2022/12_002.pdf (accessed on 22 July 2022).

28. Zoric, N. Characterization of genes and gene products influencing carotenoid metabolism in Atlantic salmon. Ph.D. Thesis; Norwegian University of Life Sciences: As, Norway, 2017.

29. Li, X.; Wang, S.; Xun, X.; Zhang, M.; Wang, S.; Li, H.; Zhao, L.; Fu, Q.; Wang, H.; Li, T. A carotenoid oxygenase is responsible for muscle coloration in scallop. Biochim. Biophys. Acta (BBA) Mol. Cell Biol. Lipids; 2019; 1864, pp. 966-975. [DOI: https://dx.doi.org/10.1016/j.bbalip.2019.03.003]

30. Lehnert, S.; Christensen, K.; Vandersteen, W.; Sakhrani, D.; Pitcher, T.; Heath, J.; Koop, B.; Heath, D.; Devlin, R. Carotenoid pigmentation in salmon: Variation in expression at BCO2-l locus controls a key fitness trait affecting red coloration. Proc. R. Soc. B; 2019; 286, 20191588. [DOI: https://dx.doi.org/10.1098/rspb.2019.1588]

31. Leeds, T.D.; Vallejo, R.L.; Weber, G.M.; Gonzalez-Pena, D.; Silverstein, J.T. Response to five generations of selection for growth performance traits in rainbow trout (Oncorhynchus mykiss). Aquaculture; 2016; 465, pp. 341-351. [DOI: https://dx.doi.org/10.1016/j.aquaculture.2016.08.036]

32. Al-Tobasei, R.; Ali, A.; Leeds, T.D.; Liu, S.; Palti, Y.; Kenney, B.; Salem, M. Identification of SNPs associated with muscle yield and quality traits using allelic-imbalance analyses of pooled RNA-Seq samples in rainbow trout. BMC Genom.; 2017; 18, 582. [DOI: https://dx.doi.org/10.1186/s12864-017-3992-z]

33. Park, J.W. Functional protein additives in surimi gels. J. Food Sci.; 1994; 59, pp. 525-527. [DOI: https://dx.doi.org/10.1111/j.1365-2621.1994.tb05554.x]

34. Misztal, I.; Tsuruta, S.; Strabel, T.; Auvray, B.; Druet, T.; Lee, D. BLUPF90 and related programs (BGF90). Proceedings of the 7th World Congress on Genetics Applied to Livestock Production; Montpellier, France, 19–23 August 2002.

35. Wang, H.; Misztal, I.; Aguilar, I.; Legarra, A.; Muir, W. Genome-wide association mapping including phenotypes from relatives without genotypes. Genet. Res.; 2012; 94, pp. 73-83. [DOI: https://dx.doi.org/10.1017/S0016672312000274] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/22624567]

36. Legarra, A.; Aguilar, I.; Misztal, I. A relationship matrix including full pedigree and genomic information. J. Dairy Sci.; 2009; 92, pp. 4656-4663. [DOI: https://dx.doi.org/10.3168/jds.2009-2061] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/19700729]

37. Aguilar, I.; Misztal, I.; Johnson, D.L.; Legarra, A.; Tsuruta, S.; Lawlor, T.J. Hot topic: A unified approach to utilize phenotypic, full pedigree, and genomic information for genetic evaluation of Holstein final score. J. Dairy Sci.; 2010; 93, pp. 743-752. [DOI: https://dx.doi.org/10.3168/jds.2009-2730]

38. Luo, Z. Computing inbreeding coefficients in large populations. Genet. Sel. Evol.; 1992; 24, pp. 305-313.

39. VanRaden, P.M. Efficient methods to compute genomic predictions. J. Dairy Sci.; 2008; 91, pp. 4414-4423. [DOI: https://dx.doi.org/10.3168/jds.2007-0980]

40. Turner, S.D. qqman: An R package for visualizing GWAS results using QQ and manhattan plots. Biorxiv; 2014; 005165.

41. Aguilar, I.; Misztal, I.; Legarra, A.; Tsuruta, S. Efficient computation of the genomic relationship matrix and other matrices used in single-step evaluation. J. Anim. Breed. Genet.; 2011; 128, pp. 422-428. [DOI: https://dx.doi.org/10.1111/j.1439-0388.2010.00912.x]

42. Gonzalez-Pena, D.; Gao, G.; Baranski, M.; Moen, T.; Cleveland, B.M.; Kenney, P.B.; Vallejo, R.L.; Palti, Y.; Leeds, T.D. Genome-wide association study for identifying loci that affect fillet yield, carcass, and body weight traits in rainbow trout (Oncorhynchus mykiss). Front. Genet.; 2016; 7, 203. [DOI: https://dx.doi.org/10.3389/fgene.2016.00203]

43. Gao, N.; Chen, Y.; Liu, X.; Zhao, Y.; Zhu, L.; Liu, A.; Jiang, W.; Peng, X.; Zhang, C.; Tang, Z. et al. Weighted single-step GWAS identified candidate genes associated with semen traits in a Duroc boar population. BMC Genom.; 2019; 20, 797. [DOI: https://dx.doi.org/10.1186/s12864-019-6164-5]

44. Juanchich, A.; Bardou, P.; Rue, O.; Gabillard, J.C.; Gaspin, C.; Bobe, J.; Guiguen, Y. Characterization of an extensive rainbow trout miRNA transcriptome by next generation sequencing. BMC Genom.; 2016; 17, 164. [DOI: https://dx.doi.org/10.1186/s12864-016-2505-9]

45. Gjerde, B.a.; Schaeffer, L. Body traits in rainbow trout: II. Estimates of heritabilities and of phenotypic and genetic correlations. Aquaculture; 1989; 80, pp. 25-44. [DOI: https://dx.doi.org/10.1016/0044-8486(89)90271-8]

46. Haffray, P.; Enez, F.; Bugeon, J.; Chapuis, H.; Dupont-Nivet, M.; Chatain, B.; Vandeputte, M. Accuracy of BLUP breeding values in a factorial mating design with mixed families and marker-based parentage assignment in rainbow trout Oncorhynchus mykiss. Aquaculture; 2018; 490, pp. 350-354. [DOI: https://dx.doi.org/10.1016/j.aquaculture.2018.03.003]

47. Yan, W.; Jang, G.-F.; Haeseleer, F.; Esumi, N.; Chang, J.; Kerrigan, M.; Campochiaro, M.; Campochiaro, P.; Palczewski, K.; Zack, D.J. Cloning and characterization of a human β,β-carotene-15, 15′-dioxygenase that is highly expressed in the retinal pigment epithelium. Genomics; 2001; 72, pp. 193-202. [DOI: https://dx.doi.org/10.1006/geno.2000.6476] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/11401432]

48. Lei, C.; Li, J.; Zheng, Z.; Du, X.; Deng, Y. Molecular cloning, expression pattern of beta-carotene 15,15-dioxygenase gene and association analysis with total carotenoid content in pearl oyster Pinctada fucata martensii. Comp. Biochem. Physiol. B Biochem. Mol. Biol.; 2019; 229, pp. 34-41. [DOI: https://dx.doi.org/10.1016/j.cbpb.2018.11.006] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/30502389]

49. dela Seña, C.; Riedl, K.M.; Narayanasamy, S.; Curley, R.W.; Schwartz, S.J.; Harrison, E.H. The human enzyme that converts dietary provitamin A carotenoids to vitamin A is a dioxygenase. J. Biol. Chem.; 2014; 289, pp. 13661-13666. [DOI: https://dx.doi.org/10.1074/jbc.M114.557710] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/24668807]

50. Huang, Y.; Zhang, L.; Wang, G.; Huang, S. De novo assembly transcriptome analysis reveals the genes associated with body color formation in the freshwater ornamental shrimps Neocaridina denticulate sinensis. Gene; 2022; 806, 145929. [DOI: https://dx.doi.org/10.1016/j.gene.2021.145929] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/34461150]

51. Chisté, R.C.; Freitas, M.; Mercadante, A.Z.; Fernandes, E. Carotenoids inhibit lipid peroxidation and hemoglobin oxidation, but not the depletion of glutathione induced by ROS in human erythrocytes. Life Sci.; 2014; 99, pp. 52-60. [DOI: https://dx.doi.org/10.1016/j.lfs.2014.01.059]

52. Grauso, M.; Lan, A.; Andriamihaja, M.; Bouillaud, F.; Blachier, F. Hyperosmolar environment and intestinal epithelial cells: Impact on mitochondrial oxygen consumption, proliferation, and barrier function in vitro. Sci. Rep.; 2019; 9, 11360. [DOI: https://dx.doi.org/10.1038/s41598-019-47851-9]

53. Ramanathan, R.; Nair, M.; Hunt, M.; Suman, S. Mitochondrial functionality and beef colour: A review of recent research. South Afr. J. Anim. Sci.; 2019; 49, pp. 9-19. [DOI: https://dx.doi.org/10.4314/sajas.v49i1.2]

54. Ottestad, S.; Sørheim, O.; Heia, K.; Skaret, J.; Wold, J.P. Effects of storage atmosphere and heme state on the color and visible reflectance spectra of salmon (Salmo salar) fillets. J. Agric. Food Chem.; 2011; 59, pp. 7825-7831. [DOI: https://dx.doi.org/10.1021/jf201150x]

55. Tintchev, F.; Kuhlmann, U.; Wackerbarth, H.; Töpfl, S.; Heinz, V.; Knorr, D.; Hildebrandt, P. Redox processes in pressurised smoked salmon studied by resonance Raman spectroscopy. Food Chem.; 2009; 112, pp. 482-486. [DOI: https://dx.doi.org/10.1016/j.foodchem.2008.05.088]

56. Mancini, R.A.; Ramanathan, R. Effects of postmortem storage time on color and mitochondria in beef. Meat Sci.; 2014; 98, pp. 65-70. [DOI: https://dx.doi.org/10.1016/j.meatsci.2014.04.007] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/24862957]

57. Tang, J.; Faustman, C.; Hoagland, T.A.; Mancini, R.A.; Seyfert, M.; Hunt, M.C. Postmortem oxygen consumption by mitochondria and its effects on myoglobin form and stability. J. Agric. Food Chem.; 2005; 53, pp. 1223-1230. [DOI: https://dx.doi.org/10.1021/jf048646o] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/15713045]

58. Banerjee, R.; Vlasie, M. Controlling the reactivity of radical intermediates by coenzyme B12-dependent methylmalonyl-CoA mutase. Biochem. Soc. Trans.; 2002; 30, pp. 621-624. [DOI: https://dx.doi.org/10.1042/bst0300621] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/12196149]

59. Yabuuchi, H.; Tanaka, K.; Maeda, M.; Takemura, M.; Oka, M.; Ohashi, R.; Tamai, I. Cloning of the dog bile salt export pump (BSEP.; ABCB11) and functional comparison with the human and rat proteins. Biopharm. Drug Dispos.; 2008; 29, pp. 441-448. [DOI: https://dx.doi.org/10.1002/bdd.629] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/18985798]

60. Paulusma, C.C.; de Waart, D.R.; Kunne, C.; Mok, K.S.; Elferink, R.P.O. Activity of the bile salt export pump (ABCB11) is critically dependent on canalicular membrane cholesterol content. J. Biol. Chem.; 2009; 284, pp. 9947-9954. [DOI: https://dx.doi.org/10.1074/jbc.M808667200] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/19228692]

61. Hayashi, H.; Inamura, K.; Aida, K.; Naoi, S.; Horikawa, R.; Nagasaka, H.; Takatani, T.; Fukushima, T.; Hattori, A.; Yabuki, T. et al. AP2 adaptor complex mediates bile salt export pump internalization and modulates its hepatocanalicular expression and transport function. Hepatology; 2012; 55, pp. 1889-1900. [DOI: https://dx.doi.org/10.1002/hep.25591]

62. Bellanti, F.; Villani, R.; Facciorusso, A.; Vendemiale, G.; Serviddio, G. Lipid oxidation products in the pathogenesis of non-alcoholic steatohepatitis. Free Radic. Biol. Med.; 2017; 111, pp. 173-185. [DOI: https://dx.doi.org/10.1016/j.freeradbiomed.2017.01.023]

63. Stofan, M.; Guo, G.L. Bile Acids and FXR: Novel Targets for Liver Diseases. Front. Med.; 2020; 7, 544. [DOI: https://dx.doi.org/10.3389/fmed.2020.00544]

64. Faustman, C.; Sun, Q.; Mancini, R.; Suman, S.P. Myoglobin and lipid oxidation interactions: Mechanistic bases and control. Meat Sci.; 2010; 86, pp. 86-94. [DOI: https://dx.doi.org/10.1016/j.meatsci.2010.04.025]

65. Liu, F.; Dai, R.; Zhu, J.; Li, X. Optimizing color and lipid stability of beef patties with a mixture design incorporating with tea catechins, carnosine, and α-tocopherol. J. Food Eng.; 2010; 98, pp. 170-177. [DOI: https://dx.doi.org/10.1016/j.jfoodeng.2009.12.023]

66. Chen, C.; Yu, Q.; Han, L.; Zhang, J.; Guo, Z. Effects of aldehyde products of lipid oxidation on the color stability and metmyoglobin reducing ability of bovine Longissimus muscle. Anim. Sci. J.; 2018; 89, pp. 810-816. [DOI: https://dx.doi.org/10.1111/asj.12993] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/29460317]

67. Esterbauer, H.; Schaur, R.J.; Zollner, H. Chemistry and biochemistry of 4-hydroxynonenal, malonaldehyde and related aldehydes. Free Radic. Biol. Med.; 1991; 11, pp. 81-128. [DOI: https://dx.doi.org/10.1016/0891-5849(91)90192-6]

68. Grattagliano, I.; Ciampi, S.A.; Portincasa, P. Gallbladder disease: Relevance of oxidative stress. Gastrointestinal Tissue; Elsevier: Amsterdam, The Netherlands, 2017; pp. 187-194.

69. Arihara, K.; Cassens, R.G.; Greaser, M.L.; Luchansky, J.B.; Mozdziak, P.E. Localization of metmyoglobin-reducing enzyme (NADH-cytochrome b5 reductase) system components in bovine skeletal muscle. Meat Sci.; 1995; 39, pp. 205-213. [DOI: https://dx.doi.org/10.1016/0309-1740(94)P1821-C]

70. Seyfert, M.; Mancini, R.A.; Hunt, M.C.; Tang, J.; Faustman, C.; Garcia, M. Color stability, reducing activity, and cytochrome c oxidase activity of five bovine muscles. J. Agric. Food Chem.; 2006; 54, pp. 8919-8925. [DOI: https://dx.doi.org/10.1021/jf061657s]

71. Mora, L.; Gallego, M.; Aristoy, M.C.; Fraser, P.D.; Toldrá, F. Peptides naturally generated from ubiquitin-60S ribosomal protein as potential biomarkers of dry-cured ham processing time. Food Control; 2015; 48, pp. 102-107. [DOI: https://dx.doi.org/10.1016/j.foodcont.2013.12.029]

72. Liu, Y.; Du, M.; Li, X.; Chen, L.; Shen, Q.; Tian, J.; Zhang, D. Role of the ubiquitin-proteasome pathway on proteolytic activity in postmortem proteolysis and tenderisation of sheep skeletal muscle. Int. J. Food Sci. Technol.; 2016; 51, pp. 2353-2359. [DOI: https://dx.doi.org/10.1111/ijfs.13214]

73. Fayemi, P.; Muchenje, V. Expression of ovine ubiquitin C-terminal hydroxylase 1, pH and colour of variety meats from head-stunned Dohne Merino sheep. S. Afr. J. Anim. Sci.; 2018; 48, pp. 88-97. [DOI: https://dx.doi.org/10.4314/sajas.v48i1.11]

74. Zheng, Q.; Zhang, Y.; Chen, Y.; Yang, N.; Wang, X.J.; Zhu, D. Systematic identification of genes involved in divergent skeletal muscle growth rates of broiler and layer chickens. BMC Genom.; 2009; 10, 87. [DOI: https://dx.doi.org/10.1186/1471-2164-10-87]

75. Harrigan, J.A.; Jacq, X.; Martin, N.M.; Jackson, S.P. Deubiquitylating enzymes and drug discovery: Emerging opportunities. Nat. Rev. Drug Discov.; 2018; 17, pp. 57-78. [DOI: https://dx.doi.org/10.1038/nrd.2017.152]

76. Damon, M.; Wyszynska-Koko, J.; Vincent, A.; Herault, F.; Lebret, B. Comparison of muscle transcriptome between pigs with divergent meat quality phenotypes identifies genes related to muscle metabolism and structure. PLoS ONE; 2012; 7, e33763. [DOI: https://dx.doi.org/10.1371/journal.pone.0033763] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/22470472]

77. Gao, H.; Antony, R.; Srinivasan, R.; Wu, P.; Wang, X.; Li, Y. UCHL1 regulates oxidative activity in skeletal muscle. PLoS ONE; 2020; 15, e0241716. [DOI: https://dx.doi.org/10.1371/journal.pone.0241716] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/33137160]

78. Gao, H.; Hartnett, S.; Li, Y. Ubiquitin C-Terminal Hydrolase L1 regulates myoblast proliferation and differentiation. Biochem. Biophys. Res. Commun.; 2017; 492, pp. 96-102. [DOI: https://dx.doi.org/10.1016/j.bbrc.2017.08.027] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/28803986]

79. Wu, W.; Gao, X.-G.; Dai, Y.; Fu, Y.; Li, X.-M.; Dai, R.-T. Post-mortem changes in sarcoplasmic proteome and its relationship to meat color traits in M. semitendinosus of Chinese Luxi yellow cattle. Food Res. Int.; 2015; 72, pp. 98-105. [DOI: https://dx.doi.org/10.1016/j.foodres.2015.03.030]

80. Zhang, Y.-M.; Zhang, X.-Z.; Wang, T.-T.; Hopkins, D.L.; Mao, Y.-W.; Liang, R.-R.; Yang, G.-F.; Xin, L.; Zhu, L.-X. Implications of step-chilling on meat color investigated using proteome analysis of the sarcoplasmic protein fraction of beef longissimus lumborum muscle. J. Integr. Agric.; 2018; 17, pp. 2118-2125. [DOI: https://dx.doi.org/10.1016/S2095-3119(18)62028-3]

81. Joseph, P.; Suman, S.P.; Rentfrow, G.; Li, S.; Beach, C.M. Proteomics of muscle-specific beef color stability. J. Agric. Food Chem.; 2012; 60, pp. 3196-3203. [DOI: https://dx.doi.org/10.1021/jf204188v]

82. Cartwright, G.M.; Scott, B. Redox regulation of an AP-1-like transcription factor, YapA, in the fungal symbiont Epichloe festucae. Eukaryot. Cell; 2013; 12, pp. 1335-1348. [DOI: https://dx.doi.org/10.1128/EC.00129-13]

83. Takahashi, M.; Yamashita, K.; Shiozawa, A.; Ichiishi, A.; Fukumori, F.; Fujimura, M. An AP-1-like transcription factor, NAP-1, regulates expression of the glutathione S-transferase and NADH: Flavin oxidoreductase genes in Neurospora crassa. Biosci. Biotechnol. Biochem.; 2010; 74, pp. 746-752. [DOI: https://dx.doi.org/10.1271/bbb.90790]

84. Jia, X.; Veiseth-Kent, E.; Grove, H.; Kuziora, P.; Aass, L.; Hildrum, K.; Hollung, K. Peroxiredoxin-6—A potential protein marker for meat tenderness in bovine longissimus thoracis muscle. J. Anim. Sci.; 2009; 87, pp. 2391-2399. [DOI: https://dx.doi.org/10.2527/jas.2009-1792]

85. Subramaniyan, S.A.; Kang, D.R.; Jung, Y.C.; Jung, J.H.; Choi, Y.I.; Lee, M.J.; Choe, H.S.; Shim, K.S. Comparative studies of meat quality traits and the proteome profile between low-and high-pH muscles in longissimus dorsi of Berkshire. Can. J. Anim. Sci.; 2017; 97, pp. 640-649. [DOI: https://dx.doi.org/10.1139/CJAS-2016-0230]

86. Wu, W.; Yu, Q.-Q.; Fu, Y.; Tian, X.-J.; Jia, F.; Li, X.-M.; Dai, R.-T. Towards muscle-specific meat color stability of Chinese Luxi yellow cattle: A proteomic insight into post-mortem storage. J. Proteom.; 2016; 147, pp. 108-118. [DOI: https://dx.doi.org/10.1016/j.jprot.2015.10.027] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/26546560]

87. Yang, X.; Wu, S.; Hopkins, D.L.; Liang, R.; Zhu, L.; Zhang, Y.; Luo, X. Proteomic analysis to investigate color changes of chilled beef longissimus steaks held under carbon monoxide and high oxygen packaging. Meat Sci.; 2018; 142, pp. 23-31. [DOI: https://dx.doi.org/10.1016/j.meatsci.2018.04.001] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/29635219]

88. Jia, W.; Zhang, R.; Liu, L.; Zhu, Z.; Xu, M.; Shi, L. Molecular mechanism of protein dynamic change for Hengshan goat meat during freezing storage based on high-throughput proteomics. Food Res. Int.; 2021; 143, 110289. [DOI: https://dx.doi.org/10.1016/j.foodres.2021.110289] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/33992389]

89. Grant, C.M.; Collinson, L.P.; Roe, J.H.; Dawes, I.W. Yeast glutathione reductase is required for protection against oxidative stress and is a target gene for yAP-1 transcriptional regulation. Mol. Microbiol.; 1996; 21, pp. 171-179. [DOI: https://dx.doi.org/10.1046/j.1365-2958.1996.6351340.x] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/8843443]

90. Lee, J.; Godon, C.; Lagniel, G.; Spector, D.; Garin, J.; Labarre, J.; Toledano, M.B. Yap1 and Skn7 control two specialized oxidative stress response regulons in yeast. J. Biol. Chem.; 1999; 274, pp. 16040-16046. [DOI: https://dx.doi.org/10.1074/jbc.274.23.16040] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/10347154]

91. Chan, K.M.; Decker, E.A. Endogenous skeletal muscle antioxidants. Crit. Rev. Food Sci. Nutr.; 1994; 34, pp. 403-426. [DOI: https://dx.doi.org/10.1080/10408399409527669]

92. Nohl, H.; Breuninger, V.; Hegner, D. Influence of mitochondrial radical formation on energy-linked respiration. Eur. J. Biochem.; 1978; 90, pp. 385-390. [DOI: https://dx.doi.org/10.1111/j.1432-1033.1978.tb12615.x]

93. Gao, X.; Wu, W.; Ma, C.; Li, X.; Dai, R. Postmortem changes in sarcoplasmic proteins associated with color stability in lamb muscle analyzed by proteomics. Eur. Food Res. Technol.; 2016; 242, pp. 527-535. [DOI: https://dx.doi.org/10.1007/s00217-015-2563-2]

94. Lee, T.T.; Ciou, J.Y.; Chen, C.L.; Yu, B. Effect of Echinacea purpurea L. on oxidative status and meat quality in Arbor Acres broilers. J. Sci. Food Agric.; 2013; 93, pp. 166-172. [DOI: https://dx.doi.org/10.1002/jsfa.5745]

95. Budanov, A.V.; Shoshani, T.; Faerman, A.; Zelin, E.; Kamer, I.; Kalinski, H.; Gorodin, S.; Fishman, A.; Chajut, A.; Einat, P. et al. Identification of a novel stress-responsive gene Hi95 involved in regulation of cell viability. Oncogene; 2002; 21, pp. 6017-6031. [DOI: https://dx.doi.org/10.1038/sj.onc.1205877]

96. Sablina, A.A.; Budanov, A.V.; Ilyinskaya, G.V.; Agapova, L.S.; Kravchenko, J.E.; Chumakov, P.M. The antioxidant function of the p53 tumor suppressor. Nat. Med.; 2005; 11, pp. 1306-1313. [DOI: https://dx.doi.org/10.1038/nm1320] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/16286925]

97. Budanov, A.V.; Sablina, A.A.; Feinstein, E.; Koonin, E.V.; Chumakov, P.M. Regeneration of peroxiredoxins by p53-regulated sestrins, homologs of bacterial AhpD. Science; 2004; 304, pp. 596-600. [DOI: https://dx.doi.org/10.1126/science.1095569] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/15105503]

98. Hanan, T.; Shaklai, N. Peroxidative interaction of myoglobin and myosin. Eur. J. Biochem.; 1995; 233, pp. 930-936. [DOI: https://dx.doi.org/10.1111/j.1432-1033.1995.930_3.x]

99. Chaijan, M.; Benjakul, S.; Visessanguan, W.; Faustman, C. Interaction between fish myoglobin and myosin in vitro. Food Chem.; 2007; 103, pp. 1168-1175. [DOI: https://dx.doi.org/10.1016/j.foodchem.2006.10.019]

100. Sellers, J.R. Myosins: A diverse superfamily. Biochim. Biophys. Acta; 2000; 1496, pp. 3-22. [DOI: https://dx.doi.org/10.1016/S0167-4889(00)00005-7]

101. Nonhoff, U.; Ralser, M.; Welzel, F.; Piccini, I.; Balzereit, D.; Yaspo, M.-L.; Lehrach, H.; Krobitsch, S. Ataxin-2 interacts with the DEAD/H-box RNA helicase DDX6 and interferes with P-bodies and stress granules. Mol. Biol. Cell; 2007; 18, pp. 1385-1396. [DOI: https://dx.doi.org/10.1091/mbc.e06-12-1120]

102. Kobayashi, T.; Winslow, S.; Sunesson, L.; Hellman, U.; Larsson, C. PKCα binds G3BP2 and regulates stress granule formation following cellular stress. PLoS ONE; 2012; 7, e35820. [DOI: https://dx.doi.org/10.1371/journal.pone.0035820]

103. Byrd, A.K.; Zybailov, B.L.; Maddukuri, L.; Gao, J.; Marecki, J.C.; Jaiswal, M.; Bell, M.R.; Griffin, W.C.; Reed, M.R.; Chib, S. Evidence that G-quadruplex DNA accumulates in the cytoplasm and participates in stress granule assembly in response to oxidative stress. J. Biol. Chem.; 2016; 291, pp. 18041-18057. [DOI: https://dx.doi.org/10.1074/jbc.M116.718478]

104. Chen, L.; Liu, B. Relationships between stress granules, oxidative stress, and neurodegenerative diseases. Oxidative Med. Cell. Longev.; 2017; 2017, 1809592. [DOI: https://dx.doi.org/10.1155/2017/1809592]

105. Gupta, V.A.; Beggs, A.H. Kelch proteins: Emerging roles in skeletal muscle development and diseases. Skelet Muscle; 2014; 4, 11. [DOI: https://dx.doi.org/10.1186/2044-5040-4-11]

106. Gupta, V.A.; Ravenscroft, G.; Shaheen, R.; Todd, E.J.; Swanson, L.C.; Shiina, M.; Ogata, K.; Hsu, C.; Clarke, N.F.; Darras, B.T. Identification of KLHL41 mutations implicates BTB-Kelch-mediated ubiquitination as an alternate pathway to myofibrillar disruption in nemaline myopathy. Am. J. Hum. Genet.; 2013; 93, pp. 1108-1117. [DOI: https://dx.doi.org/10.1016/j.ajhg.2013.10.020]

107. Kiessling, A.; Espe, M.; Ruohonen, K.; Mørkøre, T. Texture, gaping and colour of fresh and frozen Atlantic salmon flesh as affected by pre-slaughter iso-eugenol or CO₂ anaesthesia. Aquaculture; 2004; 236, pp. 645-657. [DOI: https://dx.doi.org/10.1016/j.aquaculture.2004.02.030]

108. Gagaoua, M.; Bonnet, M.; De Koning, L.; Picard, B. Reverse Phase Protein array for the quantification and validation of protein biomarkers of beef qualities: The case of meat color from Charolais breed. Meat Sci.; 2018; 145, pp. 308-319. [DOI: https://dx.doi.org/10.1016/j.meatsci.2018.06.039] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/30015160]

109. Hughes, J.; Oiseth, S.; Purslow, P.; Warner, R. A structural approach to understanding the interactions between colour, water-holding capacity and tenderness. Meat Sci.; 2014; 98, pp. 520-532. [DOI: https://dx.doi.org/10.1016/j.meatsci.2014.05.022] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/25034451]

110. Ramirez-Martinez, A.; Cenik, B.K.; Bezprozvannaya, S.; Chen, B.; Bassel-Duby, R.; Liu, N.; Olson, E.N. KLHL41 stabilizes skeletal muscle sarcomeres by nonproteolytic ubiquitination. Elife; 2017; 6, e26439. [DOI: https://dx.doi.org/10.7554/eLife.26439] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/28826497]

111. Ravenscroft, G.; Miyatake, S.; Lehtokari, V.L.; Todd, E.J.; Vornanen, P.; Yau, K.S.; Hayashi, Y.K.; Miyake, N.; Tsurusaki, Y.; Doi, H. et al. Mutations in KLHL40 are a frequent cause of severe autosomal-recessive nemaline myopathy. Am. J. Hum. Genet.; 2013; 93, pp. 6-18. [DOI: https://dx.doi.org/10.1016/j.ajhg.2013.05.004]

112. Kong, H.R.; Anthony, N.B.; Rowland, K.C.; Khatri, B.; Kong, B.C. Genome re-sequencing to identify single nucleotide polymorphism markers for muscle color traits in broiler chickens. Asian-Australas J. Anim Sci.; 2018; 31, pp. 13-18. [DOI: https://dx.doi.org/10.5713/ajas.17.0479]

113. Bailey, A.J.; Light, N.D. Connective Tissue in Meat and Meat Products; Elsevier Applied Science: Amsterdam, The Netherlands, 1989.

114. McCormick, R.J. Extracellular modifications to muscle collagen: Implications for meat quality. Poult. Sci.; 1999; 78, pp. 785-791. [DOI: https://dx.doi.org/10.1093/ps/78.5.785]

115. Ogasawara, S.; Cheng, X.W.; Inoue, A.; Hu, L.; Piao, L.; Yu, C.; Goto, H.; Xu, W.; Zhao, G.; Lei, Y. et al. Cathepsin K activity controls cardiotoxin-induced skeletal muscle repair in mice. J. Cachexia Sarcopenia Muscle; 2018; 9, pp. 160-175. [DOI: https://dx.doi.org/10.1002/jcsm.12248]

116. Kwon, H.-S.; Lee, H.-S.; Rubin, J.; Tomarev, S. Myocilin May Regulate Actin Cytoskeleton Through Components of Wnt Signaling Pathway. Investig. Ophthalmol. Vis. Sci.; 2008; 49, 5111.

117. Wu, H.; Yesilyurt, H.G.; Yoon, J.; Terman, J.R. The MICALs are a family of F-actin dismantling oxidoreductases conserved from Drosophila to humans. Sci. Rep.; 2018; 8, 937. [DOI: https://dx.doi.org/10.1038/s41598-017-17943-5] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/29343822]

118. Lundquist, M.R.; Storaska, A.J.; Liu, T.-C.; Larsen, S.D.; Evans, T.; Neubig, R.R.; Jaffrey, S.R. Redox modification of nuclear actin by MICAL-2 regulates SRF signaling. Cell; 2014; 156, pp. 563-576. [DOI: https://dx.doi.org/10.1016/j.cell.2013.12.035] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/24440334]

119. Ługowska, A.; Hetmańczyk-Sawicka, K.; Iwanicka-Nowicka, R.; Fogtman, A.; Cieśla, J.; Purzycka-Olewiecka, J.K.; Sitarska, D.; Płoski, R.; Filocamo, M.; Lualdi, S. Gene expression profile in patients with Gaucher disease indicates activation of inflammatory processes. Sci. Rep.; 2019; 9, 6060. [DOI: https://dx.doi.org/10.1038/s41598-019-42584-1] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/30988500]

120. Tarique, T.; Yang, S.; Mohsina, Z.; Qiu, J.; Yan, Z.; Chen, G.; Chen, A. Identification of genes involved in regulatory mechanism of pigments in broiler chickens. Genet. Mol. Res.; 2014; 13, pp. 7201-7216. [DOI: https://dx.doi.org/10.4238/2014.September.5.6]

121. Bartel, D.P. MicroRNAs: Genomics, biogenesis, mechanism, and function. Cell; 2004; 116, pp. 281-297. [DOI: https://dx.doi.org/10.1016/S0092-8674(04)00045-5]

122. Kutter, C.; Svoboda, P. miRNA, siRNA, piRNA: Knowns of the Unknown; Taylor & Francis: Oxfordshire, UK, 2008.

123. He, L.; Hannon, G.J. MicroRNAs: Small RNAs with a big role in gene regulation. Nat. Rev. Genet.; 2004; 5, pp. 522-531. [DOI: https://dx.doi.org/10.1038/nrg1379]

124. Clop, A.; Marcq, F.; Takeda, H.; Pirottin, D.; Tordoir, X.; Bibe, B.; Bouix, J.; Caiment, F.; Elsen, J.M.; Eychenne, F. et al. A mutation creating a potential illegitimate microRNA target site in the myostatin gene affects muscularity in sheep. Nat. Genet.; 2006; 38, pp. 813-818. [DOI: https://dx.doi.org/10.1038/ng1810]

125. Shi, J.; Sun, G. Effect of pre-miRNA-1658 gene polymorphism on chicken growth and carcass traits. Asian-Australas J. Anim. Sci.; 2017; 30, pp. 455-461. [DOI: https://dx.doi.org/10.5713/ajas.16.0305] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/27503275]

126. Li, H.; Jiang, K.; Wang, S.; Liu, X.; Kang, X.; Jiang, R.; Li, Z.; Sun, G. Assessment of correlation between pre-miRNA-1757 polymorphism and chicken performance traits. Genet. Mol. Res.; 2015; 14, pp. 12184-12195. [DOI: https://dx.doi.org/10.4238/2015.October.9.7]

127. Morimoto, R.I. Regulation of the heat shock transcriptional response: Cross talk between a family of heat shock factors, molecular chaperones, and negative regulators. Genes Dev.; 1998; 12, pp. 3788-3796. [DOI: https://dx.doi.org/10.1101/gad.12.24.3788]

128. Arya, R.; Mallik, M.; Lakhotia, S.C. Heat shock genes—Integrating cell survival and death. J. Biosci.; 2007; 32, pp. 595-610. [DOI: https://dx.doi.org/10.1007/s12038-007-0059-3] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/17536179]

129. Kaźmierczuk, A.; Kiliańska, Z.M. Plejotropowa aktywność białek szoku cieplnego The pleiotropic activity of heat-shock proteins. Postepy Hig. Med. Dosw.; 2009; 63, pp. 502-521.

130. Zychowska, M.; Jastrzebski, Z.; Chruscinski, G.; Michalowska-Sawczyn, M.; Nowak-Zaleska, A. Vitamin C, A and E supplementation decreases the expression of HSPA1A and HSPB1 genes in the leukocytes of young polish figure skaters during a 10-day training camp. J. Int. Soc. Sports Nutr.; 2015; 12, 9. [DOI: https://dx.doi.org/10.1186/s12970-015-0069-8]

131. Alexander, C.C.; Munkascy, E.; Tillmon, H.; Fraker, T.; Scheirer, J.; Holstein, D.; Lozano, D.; Khan, M.; Gidalevitz, T.; Lechleiter, J.D. et al. HspB1 Overexpression Improves Life Span and Stress Resistance in an Invertebrate Model. J. Gerontol. A Biol. Sci. Med. Sci.; 2022; 77, pp. 268-275. [DOI: https://dx.doi.org/10.1093/gerona/glab296] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/34610126]

132. Terra, L.F.; Wailemann, R.A.; Dos Santos, A.F.; Gomes, V.M.; Silva, R.P.; Laporte, A.; Meotti, F.C.; Terra, W.R.; Palmisano, G.; Lortz, S. Heat shock protein B1 is a key mediator of prolactin-induced beta-cell cytoprotection against oxidative stress. Free Radic. Biol. Med.; 2019; 134, pp. 394-405. [DOI: https://dx.doi.org/10.1016/j.freeradbiomed.2019.01.023] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/30699366]

133. Liu, X.; Xiao, W.; Jiang, Y.; Zou, L.; Chen, F.; Xiao, W.; Zhang, X.; Cao, Y.; Xu, L.; Zhu, Y. Bmal1 Regulates the Redox Rhythm of HSPB1, and Homooxidized HSPB1 Attenuates the Oxidative Stress Injury of Cardiomyocytes. Oxidative Med. Cell. Longev.; 2021; 2021, 5542815. [DOI: https://dx.doi.org/10.1155/2021/5542815]

134. Liu, X.; Liu, Y.; Jiang, Y.; Zou, L.; Liu, K.; Liu, M.; Zou, J.; Tan, S.; Wang, K.; Chen, F. Homo-oxidized HSPB1 protects cardiomyocytes against oxidative stress via targeting Keap1/Nrf-2 signaling pathway. J. Mol. Cell. Cardiol.; 2020; 140, pp. 37-38. [DOI: https://dx.doi.org/10.1016/j.yjmcc.2019.11.089]

135. Zhang, Y.; Johnson, K.; Russell, R.G.G.; Wordsworth, B.P.; Carr, A.J.; Terkeltaub, R.A.; Brown, M.A. Association of sporadic chondrocalcinosis with a− 4-basepair G-to-A transition in the 5′-untranslated region of ANKH that promotes enhanced expression of ANKH protein and excess generation of extracellular inorganic pyrophosphate. Arthritis Rheum.; 2005; 52, pp. 1110-1117. [DOI: https://dx.doi.org/10.1002/art.20978]

136. Wu, H.J.; Seib, K.L.; Srikhanta, Y.N.; Edwards, J.; Kidd, S.P.; Maguire, T.L.; Hamilton, A.; Pan, K.T.; Hsiao, H.H.; Yao, C.W. et al. Manganese regulation of virulence factors and oxidative stress resistance in Neisseria gonorrhoeae. J. Proteom.; 2010; 73, pp. 899-916. [DOI: https://dx.doi.org/10.1016/j.jprot.2009.12.001]

137. Yuan, F.L.; Xu, R.S.; Ye, J.X.; Zhao, M.D.; Ren, L.J.; Li, X. Apoptotic bodies from endplate chondrocytes enhance the oxidative stress-induced mineralization by regulating PPi metabolism. J. Cell. Mol. Med.; 2019; 23, pp. 3665-3675. [DOI: https://dx.doi.org/10.1111/jcmm.14268]

138. Gutiérrez-Luna, F.M.; Hernández-Domínguez, E.E.; Valencia-Turcotte, L.G.; Rodríguez-Sotres, R. Pyrophosphate and pyrophosphatases in plants, their involvement in stress responses and their possible relationship to secondary metabolism. Plant Sci.; 2018; 267, pp. 11-19. [DOI: https://dx.doi.org/10.1016/j.plantsci.2017.10.016] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/29362089]