Next Article in Journal
Genetic and Phenotypic Characteristics of Belted Pig Breeds: A Review
Next Article in Special Issue
Estimation of (Co) Variance Components and Genetic Parameters for Pre- and Post-Weaning Growth Traits in Dağlıç Sheep
Previous Article in Journal
Multi-Trait Single-Step Genomic Prediction for Milk Yield and Milk Components for Polish Holstein Population
Previous Article in Special Issue
Perspectives in Genome-Editing Techniques for Livestock
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Genome-Wide Association Study of Maternal Genetic Effects on Intramuscular Fat and Fatty Acid Composition in Rabbits

by
Ayman G. EL Nagar
1,2,
Imen Heddi
3,
Bolívar Samuel Sosa-Madrid
1,
Agustín Blasco
1,
Pilar Hernández
1 and
Noelia Ibáñez-Escriche
1,*
1
Institute for Animal Science and Technology, Universitat Politècnica de València, 46022 Valencia, Spain
2
Department of Animal Production, Faculty of Agriculture at Moshtohor, Benha University, Benha 13736, Egypt
3
Centro Regional de Selección y Reproducción Animal (CERSYRA), Av. del Vino, 10, 13300 Valdepeñas, Spain
*
Author to whom correspondence should be addressed.
Animals 2023, 13(19), 3071; https://doi.org/10.3390/ani13193071
Submission received: 31 August 2023 / Revised: 25 September 2023 / Accepted: 28 September 2023 / Published: 30 September 2023
(This article belongs to the Collection Applications of Quantitative Genetics in Livestock Production)

Abstract

:

Simple Summary

Ten generations of divergent selection for intramuscular fat (IMF) content was performed in rabbits. Through the generations, the IMF content has improved. Different Bayesian animal models were used to assess the importance of maternal genetic effects on the genetic control of IMF and fatty acid composition of the two rabbit lines divergently selected for IMF content. Additionally, maternal genome-wide association studies were carried out in these lines to identify genomic regions affecting the IMF and intramuscular fatty acid composition of the offspring. The maternal genetic variance explained an important part of the phenotypic variance of the IMF and intramuscular fatty acid composition. Potential candidate genes were identified on the associated genomic regions related to the traits under study.

Abstract

Maternal genetic effects (MGE) could affect meat quality traits such as intramuscular fat (IMF) and its fatty acid composition. However, it has been scarcely studied, especially in rabbits. The objectives of the present study were, first, to assess the importance of MGE on intramuscular fat and fatty acid composition by applying a Bayesian maternal animal model in two rabbit lines divergently selected for IMF. The second objective was to identify genomic regions and candidate genes of MGE that are associated with the traits of these offspring, using Bayesian methods in a Genome Wide Association Study (GWAS). Quantitative analyses were performed using data from 1982 rabbits, and 349 animals from the 9th generation and 76 dams of the 8th generation with 88,512 SNPs were used for the GWAS. The studied traits were IMF, saturated fatty acids (total SFA, C14:0; myristic acid, C16:0; palmitic acid and C18:0; stearic acid), monounsaturated fatty acids (total MUFA, C16:1n-7; palmitoleic acid and C18:1n-9; oleic acid), polyunsaturated fatty acids (total PUFA, C18:2n-6; linoleic acid, C18:3n-3; α-linolenic acid and C20:4n-6; arachidonic acid), MUFA/SFA and PUFA/SFA. The proportion of phenotypic variance explained by the maternal genetic effect ranged from 8 to 22% for IMF, depending on the model. For fatty acid composition, the proportion of phenotypic variance explained by maternal genetic effects varied from 10% (C18:0) to 46% (MUFA) in a model including both direct and additive maternal genetic effects, together with the common litter effect as a random variable. In particular, there were significant direct maternal genetic correlations for C16:0, C18:1n9, C18:2n6, SFA, MUFA, and PUFA with values ranging from −0.53 to −0.89. Relevant associated genomic regions were located on the rabbit chromosomes (OCU) OCU1, OCU5 and OCU19 containing some relevant candidates (TANC2, ACE, MAP3K3, TEX2, PRKCA, SH3GL2, CNTLN, RPGRIP1L and FTO) related to lipid metabolism, binding, and obesity. These regions explained about 1.2 to 13.9% of the total genomic variance of the traits studied. Our results showed an important maternal genetic effect on IMF and its fatty acid composition in rabbits and identified promising candidate genes associated with these traits.

1. Introduction

Many traits of interest in livestock production are affected by the mother of the individual in whom the trait is observed [1]. Maternal effects are defined by [2] as the causal influence of the maternal genotype (maternal genetic effects; MGE) or phenotype (maternal environmental effects) on the offspring phenotype. For instance, the intrauterine environment, dam’s milk production and the mothering ability represent maternal effects, which may be a combination of genetic and environmental factors that influence the offspring [3]. Maternal effects constitute a sizeable source of variation in traits of young mammals [3]. In rabbits and pigs, the importance of maternal effects for traits recorded in the fattening period is a consequence of the short interval of time between weaning and slaughtering. Additionally, the maternal effect might also influence post-weaning growth as a carry-over effect from weaning weight [4]. Similarly, the importance of the litter effect on intramuscular fat (IMF) content and its fatty acid composition has been confirmed in fat pigs such as Iberian [5] and Duroc [6]. However, these studies did not distinguish between the MGE and maternal environment or litter effect. Intramuscular fat content and its fatty acid composition are of interest in meat-producing animals since they are directly related to the quality attributes of meats, including organoleptic properties, fat consistency, shelf life, nutritional value and health issues [7,8,9]. The direct genetic control of IMF and its fatty acid composition has been determined in rabbits [10,11] and pigs [12,13,14] by both classical quantitative genetics and genome-wide association studies. However, to our knowledge, studies concerning the maternal genetic effects on IMF and fatty acid composition in livestock, particularly rabbits, are scarce. Thus, the main objectives of the present study were: (1) to quantify the importance of maternal genetic effects in the estimation of variance components for IMF and fatty acid composition using different Bayesian animal models and to determine the most appropriate model, (2) to identify genomic regions and potential candidate genes for maternal genetic effects associated with the traits studied in two rabbit lines divergently selected for IMF.

2. Materials and Methods

2.1. Ethical Statement

Animal management and the experimental procedures were approved by the Ethical Committee of the Universitat Politècnica de València, according to Council Directives 98/58/EC and 2010/63/EU, (reference number 2017/VSC/PEA/00212).

2.2. Animals and Phenotypes

The data used in this study were obtained from two rabbit lines involved in a divergent selection experiment aimed at modifying the intramuscular fat (IMF) content in the Longissimus thoracis et lumborum (LTL) muscle. This experiment was carried out over 10 generations (during the period from 2011 to 2019) in a rabbit selection nucleus located at the farm of the Institute of Animal Science and Technology of the Universitat Politècnica de València, Valencia, Spain. The base population consisted of 13 sires and 83 dams. Up to the 7th generation, each line consisted of 8 sires and 40 dams per generation. From the 8th generation on, the number increased to 10 sires and 60 dams per line. The selection was based on the offspring phenotypic value of IMF. For each dam, two full sibs (one male and one female) of first-parity offspring were sacrificed at nine weeks of age, and their IMF values were measured. The dams were then ranked according to these measurements. The top twenty percent of dams, as determined by this ranking, were selected to provide all the females for the next generation. Each sire was mated with five dams, and, in order to minimize the inbreeding, only one male offspring of the sire from the highest-ranking mating was selected for breeding the next generation. Further information on the divergent selection experiment can be found in previous publications [15,16]. The present study was carried out with data from 1982 rabbits, 166 of which were from the base population, 874 from the High-IMF line (H) and 942 from the Low-IMF line (L). The pedigree used had 3677 animals. Both lines (H and L) were contemporarily reared in the same rabbitry under the same environmental and management conditions, and fed ad libitum with a standard commercial diet, containing 15.1% crude protein, 14.5% crude fiber and 2.8% fat.
Rabbits were slaughtered at 9 weeks of age by exsanguination after electrical stunning. Carcasses were stored at 4 °C for 24 h, and then LTL muscle was excised from the carcass. LTL muscle minced and lyophilized. The IMF content of 1st to 8th generation samples were determined by near-infrared reflectance spectroscopy (NIRS) using the calibration equation developed by [17]. Twenty percent of 1st to 8th generation samples were chemically analysed by ether extraction after acid hydrolysis to validate the NIRS results. Chemical analysis was performed on all 8th and 9th generation LTL samples. IMF content was quantified as grams of IMF per 100 g of fresh muscle. Fatty acid methyl esters (Fame) were prepared as described by [18] and analyzed in a Focus Gas Chromatograph. More details of the experiment can be found in a previous publication [19]. IMF and fatty acid content were measured in grams per 100 g of fresh muscle, with fatty acids expressed as a percentage of total fatty acids. The fatty acids analyzed included the major individual fatty acids C14:0, C16:0, C18:0, C16:1n-7, C18:1n-9, C18:2n-6, C18:3n-3 and C20:4n-6 and the categories total SFA (saturated fatty acids), MUFA (monounsaturated fatty acids) and PUFA and the ratios MUFA/SFA and PUFA/SFA.

2.3. Quantifying the Importance of Maternal Genetic Effects

Five bivariate Bayesian animal models [20] of increasing complexity were implemented to estimate variance components and consequently assess the relevance of maternal genetic effects on IMF and their fatty acids. The variance components of IMF were estimated with univariate models, whereas bivariate models, always including the IMF trait to account for the selection effect, were fitted for fatty acids. All models and traits had the same fixed effects: month of slaughtering (51 levels), sex (2 levels: male and female), parity order (2 levels: 1st parity, ≥2nd parity), and method of fatty acid measurement (2 levels: NIRS, chemistry analysis). The five following bivariate models included different random effects:
Model 1. Animal model with direct additive genetic effects:
y 1 y 2 = X 1 0 0 X 2 b 1 b 2 + Z d 1 0 0 Z d 2 a d 1 a d 2 + e 1 e 2
Model 2. Animal model with direct additive genetic effects and maternal genetic effects:
y 1 y 2 = X 1 0 0 X 2 b 1 b 2 + Z d 1 0 0 Z d 2 a d 1 a d 2 + Z m 1 0 0 Z m 2 a m 1 a m 2 + e 1 e 2
Model 3. Animal model with direct additive genetic effects and common litter effects:
y 1 y 2 = X 1 0 0 X 2 b 1 b 2 + Z d 1 0 0 Z d 2 a d 1 a d 2 + W c 1 0 0 W c 2 c 1 c 2 + e 1 e 2
Model 4. Animal model with direct and maternal genetic effects as well as the common litter effects:
y 1 y 2 = X 1 0 0 X 2 b 1 b 2 + Z d 1 0 0 Z d 2 a d 1 a d 2 + Z m 1 0 0 Z m 2 a m 1 a m 2 + W c 1 0 0 W c 2 c 1 c 2 + e 1 e 2
Model 5. Animal model with direct and maternal genetic effects, and the maternal environmental effects:
y 1 y 2 = X 1 0 0 X 2 b 1 b 2 + Z d 1 0 0 Z d 2 a d 1 a d 2 + Z m 1 0 0 Z m 2 a m 1 a m 2 + W m e 1 0 0 W m e 2 m e 1 m e 2 + e 1 e 2
In all models, data were assumed to be conditionally distributed as a multinormal distribution. For instance, in model 5:
y 1 y 2 b , a d , a m , m e 1 , R 0 ~   N X 1 0 0 X 2 b 1 b 2 + Z d 1 0 0 Z d 2 a d 1 a d 2 + Z m 1 0 0 Z m 2 a m 1 a m 2 + W m e 1 0 0 W m e 2 m e 1 m e 2 ,   I n R 0
where   y 1   and y 2   are the vectors of the observations for the first (IMF) and second (fatty acid composition) traits sorted by individual and trait within an individual; b 1 and b 2   are the fixed effects vectors; a d 1 and a d 2 are the vectors of direct additive genetic effects; a m 1  and a m 2 are the maternal genetic effects vectors; m e 1 and  m e 2 are the maternal environmental effects vectors; R 0 is the residual co(variance) matrix between the two traits; and X , Z d , Z m , W c , W m e  are the incidence matrices relating the observations to the corresponding effects.
The direct additive genetic effects and the maternal genetic effects had a multinormal distribution:
a d a m ~ N 0 0 , A G 0 .
The maternal environmental effects and residuals ( e 1  and e 2 ) were assumed to be distributed as:
m e ~ N 0 , I m e C m e ,
e ~ N 0 , I n R 0 ,
where G 0 is the 4 × 4 co(variance) matrix of the direct additive and maternal genetic effects, and C m e and R 0 are the 2 × 2 co(variance) matrices of the maternal environmental and residual effects, respectively. A is the genetic relationship matrix; Ime and I n   are identity matrices of the same order as the number of levels of maternal environmental effects ( m e ) and the number of observations (n). Likewise, for model 3 the common litter effects ( c 1 and  c 2 ) were assumed to be distributed as c ~ N 0 , I c C c , where C c is the the 2 × 2 co(variance) matrix and I c is the identity matrix of the same order as the number of litters.
The univariate models for IMF included the same effects as the bivariate models, but the G 0 was a 2 × 2 co(variance) matrix of the direct additive genetic and maternal genetic effects, Cme, Cc and R 0 , where 1 × 1 co(variance) matrices, and corresponded to the variance of the maternal environmental, litter and residual effects, respectively. Flat priors were assumed for the fixed effects ( b ) and R 0 CmeCc and G 0   m a t r i c e s .
Marginal posterior distributions were estimated using Gibbs sampling with the GIBBS2F90 software [21]. The results were derived from Markov Chain Monte Carlo (MCMC) chains comprising 300,000 iterations. A burn-in period of 20,000 iterations was applied and only 1 in 200 samples was retained for subsequent inference. Subsequently, the POSTGIBBSF90 software [21] was used to estimate the following parameters: posterior means, median, mode, posterior standard deviation, and the highest posterior density region at 95% probability (HPD95%).

2.4. Genotype Data

After slaughtering the animals, about 50 g of the obliquus abdominis muscle specimens were used to extract the genomic DNA following a standard protocol [22]. The genotypes used in this study correspond to 76 8th generation dams and 350 of their 9th generation progeny, using the Affymetrix Axiom OrcunSNP Array (Affymetrix Inc., Santa Clara, CA, USA). The SNP array contains 199,692 genetic molecular markers. After quality control by removing SNPs with minor allele frequency (MAF) lower than 0.05, missing genotype rate per marker >5%, missing genotype rate per individual >3% and mapped to sex chromosomes, a total of 349 animals of the 9th generation and 76 dams of the 8th generation with 88,512 SNPs were utilized for GWAS. The missing genotypes were imputed using the software BEAGLE v4.0 [23], saving only the SNPs with imputation quality score R 2 > 0.75 .  Quality control was carried out using Axiom Analysis Suite v3.0.1.4 and ZANARDI [24].

2.5. Maternal Genome-Wide Association and Gene Annotation

The following model named MGWA [20] was fitted to identify and verify the associated maternal genomic regions. MGWA was a Bayesian multiple marker regression model used with a Bayes B model, in which the genotype of the dam is assigned to the phenotypes of her offspring:
y = 1 μ + X b + j = 1 k z j α j δ j + e
where y is the vector of offspring’s phenotypes, 1 is a vector of ones, μ is the trait mean, X is a design matrix for the systematic effects, b is a vector of fixed effects (Sex; 2 levels, Month of slaughtering; 5 levels, Parity order; 3 levels), z j is a vector of genotypic covariates for each SNP or locus j (0 = heterozygote Aa, 1 = homozygote AA or 2 = homozygote aa reference alleles), k represents the total number of SNPs remaining post-quality control (88,512), α j denotes the random allele substitution effect for SNP j, and its distribution follows a t ν ( 0 ,   σ 2 α ) , with σ 2 α obtained as described in [20], using the genetic variance of each trait estimated across all generations of the divergent selection experiment for IMF. δ j is a random binary variable (0 or 1) denoting the presence ( δ j = 1 with probability 1-π) or the absence ( δ j = 0 with probability π) of SNPs in the model for a given iteration, e is a vector of the residual values distributed as N ( 0 , I σ 2 e ) . Therefore, π probability is the proportion of SNPs that have no effect or variance in the model and is a constant calculated as a function of the number of records and the total number of SNPs as 0.9987 [11]. The MGWA Bayes B model was implemented with the GenSel software v4.9 [25]. The posterior distributions of the unknown effects comprised 500,000 iterations, discarding the first 100,000. Afterward, one in 10 samples was saved to avoid the high correlation between consecutive samples and model parameters of the marginal posterior distribution. Inference of the associations between SNPs and phenotypes was based on 1-Mb genomic windows. A total of 1973 genomic windows were analyzed across the 21 autosomes. The relevance of the association was assessed by calculating the posterior mean of the percentage of the genomic variance explained by each window and by calculating the Bayes Factor (BF) for each SNP. A genomic region was considered as associated when a window explained at least 1% of the total genomic variance and had at least one SNP with BF greater than 10 (relevant SNP). Genomic windows accounting for more than 0.5% of the trait’s genomic variance and having SNPs with a BF greater than 10 were also thought to be associated with the trait. The 1% and 0.5% thresholds, respectively, are 20 and 10 times the expected percentage of the genomic variance explained by each genomic window [26].
Candidate genes in the genomic regions associated with the traits of interest were identified using the Ensembl release 99 databases [27]. The reference genome used was the “Oryctolagus cuniculus” genome (OryCun 2.0.97). The functional annotation to uncover the biological functions of the annotated genes was performed using the DAVID v6.8 online annotation database [28,29]. Those genes reported to be involved in the physiological processes related to fat deposition, regulation of fat cell differentiation, and fatty acid and lipid metabolism were considered as potential candidate genes.

3. Results and Discussion

3.1. Quantifying the Importance of Maternal Genetic Effects

3.1.1. Descriptive Statistics

Table 1 shows the main statistical features of IMF and fatty acids. On average, IMF was 1.08 g/100 g of muscle. The percentage of PUFA had the highest average with a value of 39.29%, closely followed by the SFA with a value of 37.07% and finally the MUFA with a value of 23.77%. The most abundant individual fatty acids were linoleic acid (C18:2n6), palmitic acid (C16:0) and oleic acid (C18:1n9) with 26.95%, 26.70% and 20.45%, respectively, followed by stearic acid (C18:0) and arachidonic acid (C20:4n6) with 9.05 and 6.75%, respectively. Myristic acid (C14:0), palmitoleic acid (C16:1), and linolenic acid (C18:3n3) were found in lower percentages of less than 2% each. These results agree with other studies in rabbits [30]. The present results showed a greater percentage of PUFA, C18:2n6, C18:3n3 and C20:4n6 than in pig, beef and sheep, and lower percentages of C14:0 and C18:0 than in beef and veal [31] which explain the high quality of rabbit meat in terms of human nutrition.

3.1.2. Estimate of Genetic, Maternal, and Common Environment Parameters for IMF

Table 2 shows the posterior mean and the highest posterior density region at 95% of the heritability for direct and maternal effects, and the ratio between the litter and environmental effects and the total phenotypic variance for the IMF using five models with increasing complexity. The direct additive heritability was 0.84 when model 1 was applied (only included direct genetic effect). This heritability decreased to 0.54 when direct additive genetic and litter effects were included in model 3 and reached the lowest values for models 2, 4, and 5, which included direct and maternal genetic effects with heritability estimates of 0.44, 0.45 and 0.46, respectively. Similar results were found in Duroc, where estimates of direct heritability for IMF were reduced by approximately 33% when maternal effects were included in the model [32]. The high heritability value found with model 1 can be explained by the overestimation of the direct effect due to the absence of the maternal and the common environmental effects in the model.
A previous study [33] using the same experimental lines and model 3 reported a lower direct heritability value (0.38) compared to our study. It is worth noting that this previous study [33] only analysed data from rabbits within the first three generations. However, in a subsequent study [15], which included data from seven generations, a heritability estimate (0.53) similar to our results was obtained. To our knowledge, no other estimates of IMF heritability have been published in other lines of rabbits. In pigs, with the same model as model 3 of the present study, similar results were found in Duroc [6,13], whereas lower h 2 d   was reported for Large White, at 0.38, Pietrain, at 0.42 [34] and Iberian, at 0.26 [5], and a higher h 2 d for Landrace, at 0.67 [34]. In cattle, lower heritability (0.29) was found for Red Angus [35] and also for Hereford (0.27) with a model similar to model 5 of the present study [36]. The posterior mean of the heritability of the maternal genetic effects ( h 2 m ) was lower for models 4 (0.09) and 5 (0.08) than for model 2 (0.22), and the genetic correlations between direct and maternal genetic effects were nearly null in all models (models 2, 4, and 5). This disparity in heritability estimates could be explained by the absence in model 2 of the common litter effect and the maternal environment effect added in models 4 and 5. In Duroc pigs, [32] using a model similar to model 5 of the present study, slightly higher estimates of maternal heritability were reported for IMF (0.13). In the Rotes Höhenvieh cattle breed, a similar h 2 m (0.25) was found by applying a model equal to model 2 [37]. Nevertheless, in this study a negative direct-maternal genetic correlation ranging from −0.18 to −0.49 was found using different models [37].
Regarding the litter effect (C2), a posterior mean value of 0.18 was found with model 3, which decreased to 0.14 when model 4 including maternal genetic effects was used. This result suggests that the estimate of the litter effect found with model 3 partly includes the maternal genetic effects. The estimates for the litter effect and maternal environmental effect were similar in our study due to the reduced number of litters per dam. Similar results of C2 were reported in some breeds of pigs, with estimates of 0.14 for Large White, 0.16 for Pietrain, and 0.18 for Duroc [34,38], whereas an important mean value (0.35) was found for Iberian [5]. In contrast, a weak value (0.03) was reported for Landrace [34].
Ignoring the maternal effects leads to an overestimation of the direct heritability (models 1 and 3). Furthermore, the exclusion of the common environmental effects results in an overestimation of maternal and direct heritability (model 2). As stated by [39], in genetic evaluation models, maternal genetic effects should be taken into account; otherwise, direct heritability may be overestimated. Therefore, models 4 and 5 seem to be the most appropriate models. In fact, there is no difference between them due to the reduced number of litters per dam. Hence, we considered model 4 for the estimation of the variance components of the fatty acids.

3.1.3. Estimates of Genetic, Maternal, and Common Environment Parameters for Fatty Acid Composition

Table 3 shows the posterior mean and highest posterior density region at 95% heritability for direct and maternal effects, and the ratio of litter effect to total phenotypic variance for fatty acid composition estimated by model 4. In line with IMF, fatty acids have high to moderate direct heritability. PUFA had higher posterior mean heritability (0.43) than SFA and MUFA (0.35). The largest h 2 d for individual fatty acids were found for C16:0, C14:0 and C18:3n3, at 0.54, 0.50, and 0.50, respectively, followed by C18:2n6, C16:1, and C18:0 with direct heritability of 0.46, 0.45 and 0.40, respectively. C20:4n6 and C18:1n9 showed the lowest h2d (0.38 and 0.33, respectively). Another study with data of seven generations of these same experimental lines, using a model similar to model 3 of our study, reported a similar result of h 2 d for PUFA (0.59), lower for SFA (0.12), and higher for MUFA (0.56) [10].
Regarding individual fatty acids, [10] found similar results for most fatty acids, except for C16:0 and C18:3n3, which reported lower values (0.16 and 0.18, respectively), and for C18:1n9, which reported higher values (0.53). The size of the data and the different models used for fatty acids in our study may explain some of these observed differences. In pigs, in line with our study, moderate to high heritability were found for most fatty acid compositions, but using a model with common litter effect and without maternal genetic effects [6,40,41]. The posterior mean of the heritability of maternal genetic effects ( h 2 m ) ranged from 0.10 to 0.46. MUFA had the highest posterior mean of h 2 m (0.46), followed by PUFA (0.17) and SFA with the lowest value of 0.14. For the individual fatty acids, C16:1, C18:3n3 and C18:2n6 showed the greatest h 2 m estimates (0.46, 0.33 and 0.31, respectively). Meanwhile, the other individual fatty acids had heritabilities ranging from 0.10 to 0.20. To our knowledge, no estimates of maternal genetic effects have been previously reported for fatty acids. The posterior mean C2 did not exceed 0.09 for C18:1n9, MUFA and PUFA, and ranged from 0.11 to 0.18 for the other fatty acids, except for SFA, the only one with C2 (0.21) higher than for the maternal genetic effect. In Iberian pigs, the common environmental litter effect had an important impact on the fatty acid composition, ranging from 0.06 to 0.53 [5]. However, in Duroc pigs, the C2 values ranged from 8% (for C18:2n6) to 15% (for C20:4n6) [6]. Nevertheless, the model used in these studies did not include the maternal genetic effect.

3.1.4. Estimates of Direct Additive and Maternal Genetic Correlations between IMF and Fatty Acid Composition, and the Direct-Maternal Genetic Correlations of Fatty Acid Composition

The relationships between IMF and the fatty acid composition are of interest because both affect organoleptic and nutritional properties of meat [7]. Table 4 shows the posterior mean and the highest posterior density region at 95% of direct additive ( ρ d) and maternal genetic ( ρ m) correlations between IMF and fatty acid composition and the correlation between direct and maternal genetic effects ( ρ dm) for fatty acid composition. The posterior mean of direct genetic correlations was strong and positive between IMF and C14:0, C16:1, and SFA, ranging from 0.92 to 0.97, whereas it was strong and negative for C18:0 with an estimate of −0.80. For C18:1n9, C18:3n3 and MUFA, the ρ d was moderate and positive ranging between 0.36 and 0.52. The IMF was negatively correlated with C20:4n6 (−0.63) and with PUFA (−0.51). A previous study reported similar results for C14:0, C16:0, C18:0, C16:1, and C18:3n3, but different estimates for other fatty acid compositions, showing a lower genetic correlation for SFA (0.30), and higher correlations for C18:1, C12:2n6, C20:4n6, MUFA and PUFA compared to the present findings [10]. They used a model similar to our model 3, but their study was carried out with data from eight generations instead of ten. This led to different estimates and wider confidence intervals. Overall, our results showed more robust correlations between intramuscular fat (IMF) and meat fatty acid composition than those observed in previous studies in pigs [41] and cattle [42,43]. The posterior mean of the maternal genetic correlation ( ρ m) between IMF and its fatty acids was of the same sign as ρ d (except for C16:0 and SFA). C14:0, C16:1, C18:1n9 and MUFA had a relevant ρ m with IMF, ranging from 0.52 to 0.74, whereas C20:4n6 was strongly and negatively correlated with IMF (−0.81). For the other fatty acids, ρ m was weak with a wide HPD95%, providing little information on the true value of the parameter. The genetic correlation between the direct and maternal additive genetic effect ( ρ dm) within the trait was almost negative, except for the C20:4n6, which showed a wide HPD95%, including the zero. The direct-maternal genetic correlations were important for C16:0, C18:1n9, C18:2n6, SFA, MUFA and PUFA, giving estimates between −0.53 and −0.89. A considerable negative correlation suggests that both direct and maternal genetic effects are influenced by common genetic factors but have opposing influences on the expression of offspring traits. These results are in line with the conclusions of other studies that were reviewed by [44], who reported that direct-maternal genetic correlations showed high and negative estimates in cattle. As stated by [45], a definitive biological rationale for the conflicting interactions between maternal and direct genetic effects on meat quality traits remains elusive. Based on the results from models that included maternal genetic effects (models 2, 4 and 5), it can be concluded that maternal genetic effects play a significant role in shaping the genetic basis of the traits studied and should not be ignored in the model analyses. These findings therefore support performing GWAS analyses to identify genomic regions and potential candidate genes associated with maternal genetic effects related to the traits studied in two rabbit lines under divergent selection for IMF.

3.2. Maternal GWAS for IMF and Fatty Acid Composition

There are few studies focusing on conducting GWAS of maternal genetic effects and their relationship to animal phenotype in livestock. To our knowledge, two GWAS studies in dairy cattle investigated the maternal genetic effects on stillbirth and dystocia [46] and body weight at different ages [47]. In the present study, we assigned the genotype of the dam to her offspring and thus performed a GWAS with the phenotype of the offspring and the genotype of the dam using the MGWA model.
Table 5 shows the 1-Mb SNP window regions that accounted for more than 1% of the genomic variance, with the requirement of having SNPs with a BF greater than 10. These particular regions were considered for the identification of potential candidate genes (PCGs). The maternal genomic regions associated with the traits studied were distributed over four rabbit chromosomes (OCU1, OCU5, OCU19 and OCU20). On OCU1, an extended genomic region (35.9–37.3 Mb; three consecutive windows) was associated with PUFA, C16:1n-7, C16:0 and PUFA/SFA ratio, accounting for 3.2%, 9.9%, 12.2% and 13.9% of the genomic variance, respectively (Table 5).
This associated region harbors the potential candidate genes ADAMTSL1, CNTLN, BNC2 and SH3GL2 (Table 6). SH3GL2 and CNTLN (centlein, centrosomal protein; a neighbouring upstream gene of SH3GL2) genes are located on porcine chromosome 1 and appear to be related to synaptic vesicle endocytosis and are associated with lipid binding [48]. A previous GWAS study was performed to investigate the relationship between these two genes and obesity in pigs [49].
The genomic regions identified on OCU5, represented by a single genomic window, included genes with major biological functions directly relevant to the traits studied. Consistently, a maternal genomic region of 10.1 to 10.6 Mb on OCU5 was found to be associated with certain fatty acid composition traits. This region explained 3.5%, 5.1% and 4.3% of the genomic variance for C16:0, SFA and PUFA/SFA ratio, respectively. This result would be expected, as the main component of SFA is the palmitic fatty acid (C16:0). Two important putative candidate genes (RPGRIP1L and FTO) were found in this region; these genes have biological functions directly related to fat metabolism and deposition. The FTO (Fat mass and obesity-associated protein) gene plays an important role in adipocyte proliferation, regulation adipocyte differentiation, regulation of white and brown adipocyte differentiation, key regulators of feed intake/appetite and energy homeostasis. The association between the FTO gene and obesity in humans has been demonstrated in many populations, and this gene has been identified as an important candidate gene for fat mass, fat-deposition-related traits and obesity [50]. In addition, the FTO gene plays an essential role in postnatal growth by influencing milk fat composition and regulating lipid storage [51,52,53]. This result could partially explain the findings of [6], who found important MGE effects on the fatty acid composition of the meat of Duroc pigs. This may be due to their effects on the milk composition of the dams [54]. Other studies in lambs and calves have confirmed that meat fat composition, fat deposition and fatty acid composition of suckling animals are strongly related to the fat composition of the milk they consumed during the lactation period [55,56,57]. Based on the present results, further research is needed to test this hypothesis by evaluating the relationship between the fatty acid composition of the maternal milk fat and that of the offspring’s meat, comparing the two rabbit lines with high and low IMF.
Several studies confirmed that the FTO gene affected significantly muscle development and obesity in pigs. A polymorphism in intron 4 of the FTO gene (AM931150:g.276T>G) was found to be linked to fat-related traits such as marbling, back fat thickness, back fat between 3rd and 4th last ribs, and IMF content [58,59,60]. Our results are in agreement with those of [52] who reported that, using the candidate gene approach, polymorphism of the FTO gene located on OCU5 was significantly associated with growth and meat quality traits in rabbits. In cattle, the FTO gene has been associated with growth, carcass, marbling score and traits related to meat quality [61,62,63] as well as fat and protein yield in milk [64]. In addition, a number of human GWAS have identified the FTO gene as a primary candidate among the genomic regions that are associated with obesity [65].
The neighboring gene within this genomic region is RPGRIP1L (retinitis pigmentosa GTPase regulator-interacting protein-1-like), which is located upstream of the FTO gene. RPGRIP1L is involved in adipogenesis and insulin-regulated adipocyte metabolism. It encodes a protein with a conserved C2 domain commonly found in calcium-dependent membrane proteins. This protein binds to phospholipids, inositol polyphosphates, and intracellular proteins. Studies have shown that RPGRIP1L is associated with hyperphagic obesity, reduced leptin sensitivity and is expressed in adipose tissue. It is also associated with anthropometric and metabolic parameters, including measures of adiposity and affecting milk fat composition [64,66]. Several studies in humans and mice have found that FTO and RPGRIP1L genes have a co-expression pattern and that the effects of obesity associated with the FTO region could be modulated by changes in RPGRIP1L expression, leading to downstream effects on leptin signalling [67]. On OCU19, a genomic region (three consecutive windows) at 48.6–50.1 Mb was relevantly associated with IMF and C20:4n-6 fatty acid, accounting for about 8.0% and 1.2% of the genomic variance of the two traits, respectively. This result would be consistent with the high genetic correlation found between IMF and C20:4n-6 fatty acid in rabbits (−0.89) [10]. Of the 41 positionally annotated genes, 30 are known genes and 11 are novel genes identified by their Ensembl gene ID. Some of the existing genes in this region showed functional annotations related to lipid metabolism. In a previous GWAS, [68] found a strong association of the TANC2 gene with fat content in the dorsal muscle of common carp fish (Cyprinus carpio). Another GWAS revealed that TANC2 was significantly associated with C14:0 fatty acid in Japanese Black cattle [69]. Furthermore, this gene was differentially expressed in pigs with different IMF levels [70]. ACE (angiotensin-I-converting enzyme) gene plays a direct role in fatty acid transport, lipid localization and binding, the process of cell proliferation, differentiation, apoptosis and angiogenesis [71]. The association between polymorphism in the ACE gene with overweight, obesity and fat distribution in humans has been previously reported [72,73,74]. The MAP3K3 gene is involved in fatty acid metabolism and regulates the expression of leptin, which is secreted by adipocytes to regulate food intake and energy expenditure and plays a crucial role in maintaining body weight homeostasis [75]. Previous studies have demonstrated that elevated levels of leptin in the bloodstream can suppress fatty acid esterification in mice [76] and skeletal muscle in pigs [77]. Higher leptin expression is expected when a greater number and higher fat content of adipocytes are present [78,79]. Two divergent selection experiments in pigs have revealed significant associations between leptin and the high accumulation of IMF in muscle tissue [75,80]. A GWAS showed that the TEX2 gene was associated with C18:1n-9 fatty acid in Japanese Black cattle [81]. PITPNC1 (phosphatidylinositol transfer protein, cytoplasmic 1), also known as the M-rdgB beta gene, is incorporated in lipid metabolism through encoding for a protein of the Nir/rdgB family that is implicated in wide-ranging cellular functions such as regulation of lipid transport, metabolism, and signalling, among others [82]. A study in pigs found a relevant association between the PITPNC1 gene and the C18:3n-3 content of IMF [83]. The PRKCA (protein kinase C, alpha) gene was found to be associated with IMF accumulation in bovine, as it regulates the insulin receptor signaling pathway [84] and is involved in the regulation of IMF deposition and adipose tissue development in domestic yak [85]. In rabbits, prior to this study, GWAS analyses were conducted to investigate direct genetic effects related to IMF [11] and its fatty acid composition [19], using the same rabbit High-IMF and Low-IMF lines. These studies revealed significant genomic regions at OCU1, OCU8 and OCU13 in relation to IMF content. In particular, the most prominent region at 24.9–26.95 Mb on OCU8 accounted for 7.34% of the genomic variance [11]. In terms of fatty acid composition, a genomic region on OCU1 was identified, showing associations with C14:0, C16:0, SFA, and C18:2n6, explaining 3.5%, 11.2%, 11.3%, and 3.2% of the genomic variance, respectively. Furthermore, an additional region on OCU18 was identified, responsible for up to 8% of the genomic variance in the MUFA/SFA ratio. These studies revealed that the identified regions harbored numerous genes associated with IMF and its fatty acid composition (SCD, PLIN2, ERLIN1, and MTMR2). Manhattan plots of the IMF, C16:0, SFA, C16:1n-7, PUFA, and PUFA/SFA are shown in Figure 1, Figure 2, Figure 3, Figure 4, Figure 5 and Figure 6, respectively.

4. Conclusions

To our knowledge, this is the first study to assess the importance of maternal genetic effects and identify maternal genomic regions for intramuscular fat and intramuscular fatty acid composition traits in rabbits. Our results revealed that maternal genetic effects are not negligible and should be included in the genetic evaluation models for these traits. Ignoring maternal genetic effects led to an overestimation of direct heritability. Furthermore, excluding common environmental effects led to an overestimation of maternal and direct heritability. Therefore, models 4 and 5 seem to be the most appropriate models to estimate the genetic parameters for IMF and its fatty acid composition. The proportion of phenotypic variance explained by maternal genetic effects for IMF ranged from 8 to 22%, depending on the model. Using model 4, important maternal genetic effects for fatty acids ranged from 10% to 46%. The genetic correlation between the direct and maternal genetic effects within the trait was almost negative. It was important for C16:0, C18:1n9, C18:2n6, SFA, MUFA and PUFA, with estimates between −0.53 and −0.89. This suggests that direct and maternal genetic effects have opposing influences on the expression of offspring traits, although they share a common genetic basis. Maternal genetic effects had important influences on IMF and fatty acid composition and are recommended to be included in genetic evaluation models for these traits. A maternal genomic region associated with IMF was detected on OCU19, which contains some relevant candidate genes (TANC2, ACE, MAP3K3, TEX2 and PRKCA). Two other annotated genes (SH3GL2 and CNTLN) were detected in a genomic region on OCU1 associated with C16:0, SFA, C16:1n-7, PUFA, and PUFA/SFA fatty acids. A genomic region on OCU5 was associated with C16:0, SFA, PUFA, and PUFA/SFA fatty acids and included two important candidate genes directly related to lipid metabolism, binding, and obesity (FTO and RPGRIP1L). Our results confirm the polygenic nature of IMF and fatty acid composition traits, which are controlled by many genes distributed on many chromosomes. Further studies and fine mapping analyses are needed to validate the associations detected by GWAS. In addition, it is important to determine the relationship between the milk fatty acid composition of the dams and the fatty acid composition of the meat of their offspring, comparing the two High-IMF and Low-IMF rabbit lines.

Author Contributions

Conceptualization, N.I.-E. and P.H.; methodology, N.I.-E.; software, B.S.S.-M.; validation, A.G.E.N., I.H. and B.S.S.-M.; formal analysis, A.G.E.N. and I.H.; investigation, P.H.; resources, A.B. and P.H.; data curation, P.H.; writing—original draft preparation, A.G.E.N. and I.H.; writing—review and editing, A.B., P.H., B.S.S.-M. and N.I.-E.; visualization, N.I.-E. and A.G.E.N.; supervision, N.I.-E.; project administration, A.B.; funding acquisition, A.B. All authors have read and agreed to the published version of the manuscript.

Funding

The work was funded by the Spanish Ministry of Science and Innovation, project number PID2020-115558GB-C21, and by the Conselleria for Innovation, Universities, Science and Digital Society, project number AICO/2020/349. Ayman EL Nagar was supported by a Short-term Post-Doctoral Scholarship-Ministry of Higher Education-Cultural Affairs & Missions Sector-Egypt.

Institutional Review Board Statement

The study was conducted in accordance with Council Directives 98/58/EC and 2010/63/EU and approved by the Ethical Committee of the Universitat Politècnica de València (reference number 2017/VSC/PEA/00212).

Informed Consent Statement

Not applicable.

Data Availability Statement

The data presented in this study are available upon request to the corresponding author.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Bijma, P. Estimating Maternal Genetic Effects in Livestock. J. Anim. Sci. 2006, 84, 800–806. [Google Scholar] [CrossRef] [PubMed]
  2. Wolf, J.B.; Wade, M.J. What Are Maternal Effects (and What Are They Not)? Philos. Trans. R. Soc. B Biol. Sci. 2009, 364, 1107–1115. [Google Scholar] [CrossRef] [PubMed]
  3. Maniatis, N.; Pollott, G.E. Maternal Effects on Weight and Ultrasonically Measured Traits of Lambs in a Small Closed Suffolk Flock. Small Rumin. Res. 2002, 45, 235–246. [Google Scholar] [CrossRef]
  4. Meyer, K.; Carrick, M.J.; Donnelly, B.J.P. Genetic Parameters for Growth Traits of Australian Beef Cattle from a Multibreed Selection Experiment. J. Anim. Sci. 1993, 71, 2614–2622. [Google Scholar] [CrossRef]
  5. Ibáñez-Escriche, N.; Magallón, E.; Gonzalez, E.; Tejeda, J.F.; Noguera, J.L. Genetic Parameters and Crossbreeding Effects of Fat Deposition and Fatty Acid Profiles in Iberian Pig Lines. J. Anim. Sci. 2016, 94, 28–37. [Google Scholar] [CrossRef]
  6. Gol, S.; González-Prendes, R.; Bosch, L.; Tor, M.; Reixach, J.; Pena, R.N.; Estany, J. Linoleic Acid Metabolic Pathway Allows for an Efficient Increase of Intramuscular Fat Content in Pigs. J. Anim. Sci. Biotechnol. 2019, 10, 33. [Google Scholar] [CrossRef]
  7. Wood, J.D.; Nute, G.R.; Richardson, R.I.; Whittington, F.M.; Southwood, O.; Plastow, G.; Mansbridge, R.; Da Costa, N.; Chang, K.C. Effects of Breed, Diet and Muscle on Fat Deposition and Eating Quality in Pigs. Meat Sci. 2004, 67, 651–667. [Google Scholar] [CrossRef]
  8. Hocquette, J.F.; Gondret, F.; Baza, E.; Mdale, F.; Jurie, C.; Pethick, D.W. Intramuscular Fat Content in Meat-Producing Animals: Development, Genetic and Nutritional Control, and Identification of Putative Markers. Animal 2010, 4, 303–319. [Google Scholar] [CrossRef]
  9. Aboagye, G.; Zappaterra, M.; Laghi, L.; Dall’Olio, S.; Petracci, M.; Nanni Costa, L. Water Status in Meat from Pig Breeds Strongly Differing in Growth Performances. Food Chem. 2020, 305, 125445. [Google Scholar] [CrossRef]
  10. Martinez-Alvaro, M.; Blasco, A.; Hernandez, P. Effect of Selection for Intramuscular Fat on the Fatty Acid Composition of Rabbit Meat. Animal 2018, 12, 2002–2008. [Google Scholar] [CrossRef]
  11. Sosa-Madrid, B.S.; Hernández, P.; Blasco, A.; Haley, C.S.; Fontanesi, L.; Santacreu, M.A.; Pena, R.N.; Navarro, P.; Ibáñez-Escriche, N. Genomic Regions Influencing Intramuscular Fat in Divergently Selected Rabbit Lines. Anim. Genet. 2020, 51, 58–69. [Google Scholar] [CrossRef]
  12. Gao, Y.; Zhang, R.; Hu, X.; Li, N. Application of Genomic Technologies to the Improvement of Meat Quality of Farm Animals. Meat Sci. 2007, 77, 36–45. [Google Scholar] [CrossRef] [PubMed]
  13. Ros-Freixedes, R.; Gol, S.; Pena, R.N.; Tor, M.; Ibáñez-Escriche, N.; Dekkers, J.C.M.; Estany, J. Genome-Wide Association Study Singles out SCD and LEPR as the Two Main Loci Influencing Intramuscular Fat Content and Fatty Acid Composition in Duroc Pigs. PLoS ONE 2016, 11, e0152496. [Google Scholar] [CrossRef] [PubMed]
  14. Pena, R.N.; Noguera, J.L.; García-Santana, M.J.; González, E.; Tejeda, J.F.; Ros-Freixedes, R.; Ibáñez-Escriche, N. Five Genomic Regions Have a Major Impact on Fat Composition in Iberian Pigs. Sci. Rep. 2019, 9, 2031. [Google Scholar] [CrossRef] [PubMed]
  15. Martínez-Álvaro, M.; Hernández, P.; Blasco, A. Divergent Selection on Intramuscular Fat in Rabbits: Responses to Selection and Genetic Parameters. J. Anim. Sci. 2016, 94, 4993–5003. [Google Scholar] [CrossRef]
  16. Zubiri-Gaitán, A.; Blasco, A.; Ccalta, R.; Satué, K.; Hernández, P. Intramuscular Fat Selection in Rabbits Modifies the Fatty Acid Composition of Muscle and Liver Tissues. Animals 2022, 12, 893. [Google Scholar] [CrossRef] [PubMed]
  17. Zomeño, C.; Juste, V.; Hernández, P. Application of NIRS for Predicting Fatty Acids in Intramuscular Fat of Rabbit. Meat Sci. 2012, 91, 155–159. [Google Scholar] [CrossRef]
  18. O’Fallon, J.V.; Busboom, J.R.; Nelson, M.L.; Gaskins, C.T. A Direct Method for Fatty Acid Methyl Ester Synthesis: Application to Wet Meat Tissues, Oils, and Feedstuffs. J. Anim. Sci. 2007, 85, 1511–1521. [Google Scholar] [CrossRef]
  19. Laghouaouta, H.; Sosa-Madrid, B.S.; Zubiri-Gaitán, A.; Hernández, P.; Blasco, A. Novel Genomic Regions Associated with Intramuscular Fatty Acid Composition in Rabbits. Animals 2020, 10, 2090. [Google Scholar] [CrossRef]
  20. Blasco, A. Bayesian Data Analysis for Animal Scientists; Springer: Berlin/Heidelberg, Germany, 2017. [Google Scholar]
  21. Misztal, I.; Tsuruta, S.; Strabel, T.; Auvray, B.; Druet, T.; Lee, D.H. Blupf90 and Related Programs (Bgf90). In Proceedings of the 7th World Congress on Genetics Applied to Livestock Production, Montpellier, France, 19–23 August 2002; pp. 743–744. [Google Scholar]
  22. Green, M.R.; Sambrook, J. Quantification of DNA and RNA by Real-Time Polymerase Chain Reaction. In Molecular Cloning: Laboratory Manual, 4th ed.; Cold Spring Harb. Lab. Press: New York, NY, USA, 2012; pp. 631–679. [Google Scholar]
  23. Browning, B.L.; Browning, S.R. Genotype Imputation with Millions of Reference Samples. Am. J. Hum. Genet. 2016, 98, 116–126. [Google Scholar] [CrossRef]
  24. Marras, G.; Rossoni, A.; Schwarzenbacher, H.; Biffani, S.; Biscarini, F.; Nicolazzi, E.L. Zanardi: An Open-Source Pipeline for Multiple-Species Genomic Analysis of SNP Array Data. Anim. Genet. 2017, 48, 121. [Google Scholar] [CrossRef] [PubMed]
  25. Garrick, D.J.; Fernando, R.L. Implementing a QTL Detection Study (GWAS) Using Genomic Prediction Methodology. In Genome-Wide Association Studies and Genomic Prediction; Springer Nature: Berlin/Heidelberg, Germany, 2013; pp. 275–298. [Google Scholar]
  26. Kass, R.E.; Raftery, A.E. Bayes Factors. J. Am. Stat. Assoc. 1995, 90, 773–795. [Google Scholar] [CrossRef]
  27. Cunningham, F.; Achuthan, P.; Akanni, W.; Allen, J.; Amode, M.R.; Armean, I.M.; Bennett, R.; Bhai, J.; Billis, K.; Boddu, S.; et al. Ensembl 2019. Nucleic Acids Res. 2019, 47, D745–D751. [Google Scholar] [CrossRef]
  28. Huang, D.W.; Sherman, B.T.; Lempicki, R.A. Bioinformatics Enrichment Tools: Paths toward the Comprehensive Functional Analysis of Large Gene Lists. Nucleic Acids Res. 2009, 37, 1–13. [Google Scholar] [CrossRef] [PubMed]
  29. Jiao, X.; Sherman, B.T.; Huang, D.W.; Stephens, R.; Baseler, M.W.; Lane, H.C.; Lempicki, R.A. DAVID-WS: A Stateful Web Service to Facilitate Gene/Protein List Analysis. Bioinformatics 2012, 28, 1805–1806. [Google Scholar] [CrossRef]
  30. Hernández, P.; Zotte, A.D. Influence of Diet on Rabbit Meat Quality. In Nutrition of the Rabbit; CABI: Wallingford, UK, 2010; pp. 163–178. [Google Scholar]
  31. Dalle Zotte, A. Perception of Rabbit Meat Quality and Major Factors Influencing the Rabbit Carcass and Meat Quality. Livest. Prod. Sci. 2002, 75, 11–32. [Google Scholar] [CrossRef]
  32. Hoque, M.A.; Kadowaki, H.; Shibata, T.; Suzuki, K. Maternal and Direct Genetic Parameters for Production Traits and Maternal Correlations among Production and Feed Efficiency Traits in Duroc Pigs. Asian-Australas. J. Anim. Sci. 2008, 21, 961–966. [Google Scholar] [CrossRef]
  33. Zomeño, C.; Hernández, P.; Blasco, A. Divergent Selection for Intramuscular Fat Content in Rabbits. I. Direct Response to Selection. J. Anim. Sci. 2013, 91, 4526–4531. [Google Scholar] [CrossRef]
  34. Knapp, P.; Willam, A.; Sölkner, J. Genetic Parameters for Lean Meat Content and Meat Quality Traits in Different Pig Breeds. Livest. Prod. Sci. 1997, 52, 69–73. [Google Scholar] [CrossRef]
  35. Mcallister, C.M.; Speidel, S.E.; Crews, D.H., Jr.; Enns, R.M. Genetic Parameters for Intramuscular Fat Percentage, Marbling Score, Scrotal Circumference, and Heifer Pregnancy in Red Angus Cattle 1. J. Anim. Sci. 2011, 89, 2068–2072. [Google Scholar] [CrossRef]
  36. Torres-vazquez, J.A. Genetic Parameters for Docility, Weaning Weight, Yearling Weight and Intramuscular Fat Percentage in Hereford Cattle. J. Anim. Sci. 2016, 94, 21–27. [Google Scholar] [CrossRef] [PubMed]
  37. Halli, K.; Bohlouli, M.; Schulz, L.; Sundrum, A.; König, S. Estimation of Direct and Maternal Genetic Effects and Annotation of Potential Candidate Genes for Weight and Meat Quality Traits in a Genotyped Outdoor Dual-Purpose Cattle Breed. Transl. Anim. Sci. 2022, 6, txac022. [Google Scholar] [CrossRef]
  38. Suzuki, K.; Irie, M.; Kadowaki, H.; Shibata, T.; Kumagai, M.; Nishida, A. Genetic Parameter Estimates of Meat Quality Traits in Duroc Pigs Selected for Average Daily Gain, Longissimus Muscle Area, Backfat Thickness, and Intramuscular Fat Content. J. Anim. Sci. 2005, 83, 2058–2065. [Google Scholar] [CrossRef] [PubMed]
  39. Meyer, K.; Carrick, M.J.; Donnelly, B.J.P. Genetic Parameters for Milk Production of Australian Beef Cows and Weaning Weight of Their Calves. J. Anim. Sci. 1994, 72, 1155–1165. [Google Scholar] [CrossRef] [PubMed]
  40. Sellier, P.; Maignel, L.; Bidanel, J.P. Genetic Parameters for Tissue and Fatty Acid Composition of Backfat, Perirenal Fat and Longissimus Muscle in Large White and Landrace Pigs. Animal 2010, 4, 497–504. [Google Scholar] [CrossRef] [PubMed]
  41. Ros-Freixedes, R.; Reixach, J.; Bosch, L.; Tor, M.; Estany, J. Genetic Correlations of Intramuscular Fat Content and Fatty Acid Composition among Muscles and with Subcutaneous Fat in Duroc Pigs. J. Anim. Sci. 2014, 92, 5417–5425. [Google Scholar] [CrossRef]
  42. Nogi, T.; Honda, T.; Mukai, F.; Okagaki, T.; Oyama, K. Heritabilities and Genetic Correlations of Fatty Acid Compositions in Longissimus Muscle Lipid with Carcass Traits in Japanese Black Cattle. J. Anim. Sci. 2011, 89, 615–621. [Google Scholar] [CrossRef]
  43. Buchanan, J.W.; Reecy, J.M.; Garrick, D.J.; Duan, Q.; Beitz, D.C.; Mateescu, R.G. Genetic Parameters and Genetic Correlations among Triacylglycerol and Phospholipid Fractions in Angus Cattle. J. Anim. Sci. 2015, 93, 522–528. [Google Scholar] [CrossRef]
  44. Quintanilla, R.; Piedrafita, J. Efectos Maternos En El Peso Al Destete Del Ganado Vacuno de Carne: Revisión. Itea 2000, 96, 7–39. [Google Scholar]
  45. Lee, C. On the Negative Estimates of Direct and Maternal Genetic Correlation—A Review. Asian-Australas. J. Anim. Sci. 2002, 15, 1222–1226. [Google Scholar] [CrossRef]
  46. Olsen, H.G.; Hayes, B.J.; Kent, M.P.; Nome, T.; Svendsen, M.; Lien, S. A Genome Wide Association Study for QTL Affecting Direct and Maternal Effects of Stillbirth and Dystocia in Cattle. Anim. Genet. 2010, 41, 273–280. [Google Scholar] [CrossRef] [PubMed]
  47. Yin, T.; König, S. Genome-Wide Associations and Detection of Potential Candidate Genes for Direct Genetic and Maternal Genetic Effects Influencing Dairy Cattle Body Weight at Different Ages. Genet. Sel. Evol. 2019, 51, 4. [Google Scholar] [CrossRef] [PubMed]
  48. Huttner, W.B.; Schmidt, A. Lipids, Lipid Modification and Lipid-Protein Interaction in Membrane Budding and Fission—Insights from the Roles of Endophilin A1 and Synaptophysin in Synaptic Vesicle Endocytosis. Curr. Opin. Neurobiol. 2000, 10, 543–551. [Google Scholar] [CrossRef] [PubMed]
  49. Kogelman, L.J.A.; Pant, S.D.; Fredholm, M.; Kadarmideen, H.N. Systems Genetics of Obesity in an F2 Pig Model by Genome-Wide Association, Genetic Network, and Pathway Analyses. Front. Genet. 2014, 5, 214. [Google Scholar] [CrossRef]
  50. Hess, M.E.; Brüning, J.C. The Fat Mass and Obesity-Associated (FTO) Gene: Obesity and Beyond? Biochim. Biophys. Acta Mol. Basis Dis. 2014, 1842, 2039–2047. [Google Scholar] [CrossRef]
  51. Fan, B.; Du, Z.-Q.; Rothschild, M.F. The Fat Mass and Obesity-Associated (FTO) Gene Is Associated with Intramuscular Fat Content and Growth Rate in the Pig. Anim. Biotechnol. 2009, 20, 58–70. [Google Scholar] [CrossRef]
  52. Zhang, G.W.; Gao, L.; Chen, S.Y.; Zhao, X.B.; Tian, Y.F.; Wang, X.; Deng, X.S.; Lai, S.J. Single Nucleotide Polymorphisms in the FTO Gene and Their Association with Growth and Meat Quality Traits in Rabbits. Gene 2013, 527, 553–557. [Google Scholar] [CrossRef]
  53. Reyer, H.; Varley, P.F.; Murani, E.; Ponsuksili, S.; Wimmers, K. Genetics of Body Fat Mass and Related Traits in a Pig Population Selected for Leanness. Sci. Rep. 2017, 7, 9118. [Google Scholar] [CrossRef]
  54. Jin, C.; Fang, Z.; Lin, Y.; Che, L.; Wu, C.; Xu, S.; Feng, B.; Li, J.; Wu, D. Influence of Dietary Fat Source on Sow and Litter Performance, Colostrum and Milk Fatty Acid Profile in Late Gestation and Lactation. Anim. Sci. J. 2017, 88, 1768–1778. [Google Scholar] [CrossRef]
  55. Velasco, S.; Cañeque, V.; Lauzurica, S.; Pérez, C.; Huidobro, F. Effect of Different Feeds on Meat Quality and Fatty Acid Composition of Lambs Fattened at Pasture. Meat Sci. 2004, 66, 457–465. [Google Scholar] [CrossRef]
  56. Moreno, T.; Varela, A.; Oliete, B.; Carballo, J.A.; Sánchez, L.; Montserrat, L. Nutritional Characteristics of Veal from Weaned and Unweaned Calves: Discriminatory Ability of the Fat Profile. Meat Sci. 2006, 73, 209–217. [Google Scholar] [CrossRef] [PubMed]
  57. Alfaia, C.M.M.; Castro, M.L.F.; Martins, S.I.V.; Portugal, A.P.V.; Alves, S.P.A.; Fontes, C.M.G.A.; Bessa, R.J.B.; Prates, J.A.M. Effect of Slaughter Season on Fatty Acid Composition, Conjugated Linoleic Acid Isomers and Nutritional Value of Intramuscular Fat in Barrosã-PDO Veal. Meat Sci. 2007, 75, 44–52. [Google Scholar] [CrossRef]
  58. Fontanesi, L.; Scotti, E.; Buttazzoni, L.; Davoli, R.; Russo, V. The Porcine Fat Mass and Obesity Associated (FTO) Gene Is Associated with Fat Deposition in Italian Duroc Pigs. Anim. Genet. 2009, 40, 90–93. [Google Scholar] [CrossRef] [PubMed]
  59. Fontanesi, L.; Scotti, E.; Buttazzoni, L.; Dall’Olio, S.; Davoli, R.; Russo, V. A Single Nucleotide Polymorphism in the Porcine Cathepsin K (CTSK) Gene Is Associated with Back Fat Thickness and Production Traits in Italian Duroc Pigs. Mol. Biol. Rep. 2010, 37, 491–495. [Google Scholar] [CrossRef] [PubMed]
  60. Fontanesi, L.; Schiavo, G.; Galimberti, G.; Calò, D.G.; Russo, V. A Genomewide Association Study for Average Daily Gain in Italian Large White Pigs. J. Anim. Sci. 2014, 92, 1385–1394. [Google Scholar] [CrossRef]
  61. Rempel, L.A.; Casas, E.; Shackelford, S.D.; Wheeler, T.L. Relationship of Polymorphisms within Metabolic Genes and Carcass Traits in Crossbred Beef Cattle. J. Anim. Sci. 2012, 90, 1311–1316. [Google Scholar] [CrossRef] [PubMed]
  62. Chung, E.-R. Novel SNP in the Coding Region of the FTO Gene Is Associated with Marbling Score in Hanwoo (Korean Cattle). J. Anim. Sci. Technol. 2014, 56, 27. [Google Scholar] [CrossRef] [PubMed]
  63. Jevsinek Skok, D.; Kunej, T.; Kovac, M.; Malovrh, S.; Potocnik, K.; Petric, N.; Zgur, S.; Dovc, P.; Horvat, S. FTO Gene Variants Are Associated with Growth and Carcass Traits in Cattle. Anim. Genet. 2016, 47, 219–222. [Google Scholar] [CrossRef] [PubMed]
  64. Zielke, L.G.; Bortfeldt, R.H.; Reissmann, M.; Tetens, J.; Thaller, G.; Brockmann, G.A. Impact of Variation at the FTO Locus on Milk Fat Yield in Holstein Dairy Cattle. PLoS ONE 2013, 8, e63406. [Google Scholar] [CrossRef]
  65. Babenko, V.; Babenko, R.; Gamieldien, J.; Markel, A. FTO Haplotyping Underlines High Obesity Risk for European Populations. BMC Med. Genom. 2019, 12, 46. [Google Scholar] [CrossRef]
  66. Tews, D.; Fischer-Posovszky, P.; Wabitsch, M. Regulation of FTO and FTM Expression during Human Preadipocyte Differentiation. Horm. Metab. Res. 2011, 43, 17–21. [Google Scholar] [CrossRef] [PubMed]
  67. Vogan, K. RPGRIP1L, FTO and Obesity. Nat. Genet. 2014, 46, 532. [Google Scholar] [CrossRef]
  68. Zheng, X.; Kuang, Y.; Lv, W.; Cao, D.; Sun, Z.; Sun, X. Genome-Wide Association Study for Muscle Fat Content and Abdominal Fat Traits in Common Carp (Cyprinus Carpio). PLoS ONE 2016, 11, e0169127. [Google Scholar] [CrossRef] [PubMed]
  69. Ishii, A.; Yamaji, K.; Uemoto, Y.; Sasago, N.; Kobayashi, E.; Kobayashi, N.; Matsuhashi, T.; Maruyama, S.; Matsumoto, H.; Sasazaki, S.; et al. Genome-Wide Association Study for Fatty Acid Composition in Japanese Black Cattle. Anim. Sci. J. 2013, 84, 675–682. [Google Scholar] [CrossRef]
  70. Liu, J.; Damon, M.; Guitton, N.; Guisle, I.; Ecolan, P.; Vincent, A.; Cherel, P.; Gondret, F. Differentially-Expressed Genes in Pig Longissimus Muscles with Contrasting Levels of Fat, as Identified by Combined Transcriptomic, Reverse Transcription PCR, and Proteomic Analyses. J. Agric. Food Chem. 2009, 57, 3808–3817. [Google Scholar] [CrossRef]
  71. Leung, P.S. The Peptide Hormone Angiotensin II: Its New Functions in Tissues and Organs. Curr. Protein Pept. Sci. 2004, 5, 267–273. [Google Scholar] [CrossRef]
  72. Riera-Fortuny, C.; Real, J.T.; Chaves, F.J.; Morales-Suárez-Varela, M.; Martínez-Triguero, M.L.; Morillas-Arino, C.; Hernández-Mijares, A. The Relation between Obesity, Abdominal Fat Deposit and the Angiotensin-Converting Enzyme Gene I/D Polymorphism and Its Association with Coronary Heart Disease. Int. J. Obes. 2005, 29, 78–84. [Google Scholar] [CrossRef]
  73. Lagou, V.; Manios, Y.; Moran, C.N.; Bailey, M.E.S.; Grammatikaki, E.; Oikonomou, E.; Ioannou, E.; Moschonis, G.; Wilson, R.H.; Pitsiladis, Y.P. Developmental Changes in Adiposity in Toddlers and Preschoolers in the GENESIS Study and Associations with the ACE I/D Polymorphism. Int. J. Obes. 2007, 31, 1052–1060. [Google Scholar] [CrossRef]
  74. Bienertova-Vasku, J.; Bienert, P.; Sablikova, L.; Slovackova, L.; Forejt, M.; Piskackova, Z.; Kucerova, L.; Heczkova, K.; Brazdova, Z.; Vasku, A. Effect of ID ACE Gene Polymorphism on Dietary Composition and Obesity-Related Anthropometric Parameters in the Czech Adult Population. Genes Nutr. 2009, 4, 207–213. [Google Scholar] [CrossRef]
  75. Hamill, R.M.; Aslan, O.; Mullen, A.M.; O’Doherty, J.V.; McBryan, J.; Morris, D.G.; Sweeney, T. Transcriptome Analysis of Porcine, M. Semimembranosus Divergent in Intramuscular Fat as a Consequence of Dietary Protein Restriction. BMC Genom. 2013, 14, 453. [Google Scholar] [CrossRef]
  76. Steinberg, G.R.; Dyck, D.J.; Calles-Escandon, J.; Tandon, N.N.; Luiken, J.J.F.P.; Glatz, J.F.C.; Bonen, A. Chronic Leptin Administration Decreases Fatty Acid Uptake and Fatty Acid Transporters in Rat Skeletal Muscle. J. Biol. Chem. 2002, 277, 8854–8860. [Google Scholar] [CrossRef] [PubMed]
  77. Ramsay, T.G. Porcine Leptin Alters Isolated Adipocyte Glucose and Fatty Acid Metabolism. Domest. Anim. Endocrinol. 2004, 26, 11–21. [Google Scholar] [CrossRef]
  78. Steiner, A.A.; Romanovsky, A.A. Leptin: At the Crossroads of Energy Balance and Systemic Inflammation. Prog. Lipid Res. 2007, 46, 89–107. [Google Scholar] [CrossRef] [PubMed]
  79. Yau, M.H.; Wang, Y.; Lam, K.S.L.; Zhang, J.; Wu, D.; Xu, A. A Highly Conserved Motif within the NH2-Terminal Coiled-Coil Domain of Angiopoietin-like Protein 4 Confers Its Inhibitory Effects on Lipoprotein Lipase by Disrupting the Enzyme Dimerization. J. Biol. Chem. 2009, 284, 11942–11952. [Google Scholar] [CrossRef]
  80. Cánovas, A.; Quintanilla, R.; Amills, M.; Pena, R.N. Muscle Transcriptomic Profiles in Pigs with Divergent Phenotypes for Fatness Traits. BMC Genom. 2010, 11, 372. [Google Scholar] [CrossRef] [PubMed]
  81. Uemoto, Y.; Abe, T.; Tameoka, N.; Hasebe, H.; Inoue, K.; Nakajima, H.; Shoji, N.; Kobayashi, M.; Kobayashi, E. Whole-Genome Association Study for Fatty Acid Composition of Oleic Acid in Japanese Black Cattle. Anim. Genet. 2011, 42, 141–148. [Google Scholar] [CrossRef] [PubMed]
  82. Lev, S. The Role of the Nir/RdgB Protein Family in Membrane Trafficking and Cytoskeleton Remodeling. Exp. Cell Res. 2004, 297, 1–10. [Google Scholar] [CrossRef]
  83. Muñoz, M.; Alves, E.; Corominas, J.; Folch, J.M.; Casellas, J.; Noguera, J.L.; Silió, L.; Fernández, A.I. Survey of SSC12 Regions Affecting Fatty Acid Composition of Intramuscular Fat Using High-Density SNP Data. Front. Genet. 2012, 2, 101. [Google Scholar] [CrossRef]
  84. Lim, D.; Kim, N.K.; Park, H.S.; Lee, S.H.; Cho, Y.M.; Oh, S.J.; Kim, T.H.; Kim, H. Identification of Candidate Genes Related to Bovine Marbling Using Pro-Tein-Protein Interaction Networks. Int. J. Biol. Sci. 2011, 7, 992–1002. [Google Scholar] [CrossRef]
  85. Wang, H.; Zhong, J.; Zhang, C.; Chai, Z.; Cao, H.; Wang, J.; Zhu, J.; Wang, J.; Ji, Q. The Whole-Transcriptome Landscape of Muscle and Adipose Tissues Reveals the CeRNA Regulation Network Related to Intramuscular Fat Deposition in Yak. BMC Genom. 2020, 21, 347. [Google Scholar] [CrossRef]
Figure 1. Manhattan plot of the genome-wide association study for intramuscular fat of Longissimus thoracis et lumborum muscle in rabbits using the Bayes factors by each SNP along the rabbit chromosomes under the MGWA model. The Bayes factor threshold of 10 is indicated by the black dashed line.
Figure 1. Manhattan plot of the genome-wide association study for intramuscular fat of Longissimus thoracis et lumborum muscle in rabbits using the Bayes factors by each SNP along the rabbit chromosomes under the MGWA model. The Bayes factor threshold of 10 is indicated by the black dashed line.
Animals 13 03071 g001
Figure 2. Manhattan plot of the genome-wide association study for palmitic fatty acid (C16:0) of Longissimus thoracis et lumborum muscle in rabbits using the Bayes factors by each SNP along the rabbit chromosomes under the MGWA model. The Bayes factor threshold of 10 is indicated by the black dashed line.
Figure 2. Manhattan plot of the genome-wide association study for palmitic fatty acid (C16:0) of Longissimus thoracis et lumborum muscle in rabbits using the Bayes factors by each SNP along the rabbit chromosomes under the MGWA model. The Bayes factor threshold of 10 is indicated by the black dashed line.
Animals 13 03071 g002
Figure 3. Manhattan plot of the genome-wide association study for saturated fatty acids (SFA) of Longissimus thoracis et lumborum muscle in rabbits using the Bayes factors by each SNP along the rabbit chromosomes under the MGWA model. The Bayes factor threshold of 10 is indicated by the black dashed line.
Figure 3. Manhattan plot of the genome-wide association study for saturated fatty acids (SFA) of Longissimus thoracis et lumborum muscle in rabbits using the Bayes factors by each SNP along the rabbit chromosomes under the MGWA model. The Bayes factor threshold of 10 is indicated by the black dashed line.
Animals 13 03071 g003
Figure 4. Manhattan plot of the genome-wide association study for palmitoleic fatty acid (C16:1n7) of Longissimus thoracis et lumborum muscle in rabbits using the Bayes factors by each SNP along the rabbit chromosomes under the MGWA model. The Bayes factor threshold of 10 is indicated by the black dashed line.
Figure 4. Manhattan plot of the genome-wide association study for palmitoleic fatty acid (C16:1n7) of Longissimus thoracis et lumborum muscle in rabbits using the Bayes factors by each SNP along the rabbit chromosomes under the MGWA model. The Bayes factor threshold of 10 is indicated by the black dashed line.
Animals 13 03071 g004
Figure 5. Manhattan plot of the genome-wide association study for polyunsaturated fatty acids (PUFA) of Longissimus thoracis et lumborum muscle in rabbits using the Bayes factors by each SNP along the rabbit chromosomes under the MGWA model. The Bayes factor threshold of 10 is indicated by the black dashed line.
Figure 5. Manhattan plot of the genome-wide association study for polyunsaturated fatty acids (PUFA) of Longissimus thoracis et lumborum muscle in rabbits using the Bayes factors by each SNP along the rabbit chromosomes under the MGWA model. The Bayes factor threshold of 10 is indicated by the black dashed line.
Animals 13 03071 g005
Figure 6. Manhattan plot of the genome-wide association study for the ratio between polyunsaturated and saturated fatty acids (PUFA/SFA) of Longissimus thoracis et lumborum muscle in rabbits using the Bayes factors by each SNP along the rabbit chromosomes under the MGWA model. The Bayes factor threshold of 10 is indicated by the black dashed line.
Figure 6. Manhattan plot of the genome-wide association study for the ratio between polyunsaturated and saturated fatty acids (PUFA/SFA) of Longissimus thoracis et lumborum muscle in rabbits using the Bayes factors by each SNP along the rabbit chromosomes under the MGWA model. The Bayes factor threshold of 10 is indicated by the black dashed line.
Animals 13 03071 g006
Table 1. Mean, median, standard deviation (SD), and coefficient of variation (CV %) for intramuscular fat (g/100 g of muscle) and fatty acid composition of Longissimus thoracis et lumborum muscle (% of total fatty acids).
Table 1. Mean, median, standard deviation (SD), and coefficient of variation (CV %) for intramuscular fat (g/100 g of muscle) and fatty acid composition of Longissimus thoracis et lumborum muscle (% of total fatty acids).
TraitMeanMedianSDCV (%)
IMF1.081.060.2321.30
C14:01.251.310.5241.60
C16:026.7026.981.575.88
C18:09.059.021.2613.92
SFA37.0737.011.674.50
C16:1n71.401.350.7855.71
C18:1n920.4520.923.2615.94
MUFA23.7724.524.5519.14
C18:2n626.9527.062.9711.02
C18:3n31.641.690.4326.22
C20:4n66.576.121.9930.29
PUFA39.2938.654.8612.37
MUFA/SFA0.640.660.034.69
PUFA/SFA1.061.040.065.66
IMF: intramuscular fat; SFA: saturated fatty acids = C14:0 + C15:0 + C16:0 + C17:0 + C18:0; MUFA: monounsaturated fatty acids = C16:1n7 + C18:1n7 + C18:1n9; PUFA: polyunsaturated fatty acids = C18:2n6 + C18:3n3 + C20:2n6 + C20:3n6 + C20:4n6 + C20:5n3 + C22:4n6 + C22:5n3 + C22:6n3; C14:0: myristic acid; C16:0: palmitic acid; C18:0: stearic acid; C16:1n7: palmitoleic acid; C18:1n9: oleic acid; C18:2n-6: linoleic acid; C18:3n3: α-linolenic acid; C20:4n-6: arachidonic acid. The total number of slaughtered rabbits is 1982.
Table 2. Marginal posterior mean and highest posterior density region at 95% (between brackets) of the heritability for direct ( h 2 d ) and maternal ( h 2 m ) genetic effects, the correlation between direct and maternal genetic effects ( ρ dm ) and the ratio between the litter (C2) and maternal environmental (Me2) effects and the total phenotypic variance for the Intramuscular Fat of Longissimus thoracis et lumborum muscle using five models.
Table 2. Marginal posterior mean and highest posterior density region at 95% (between brackets) of the heritability for direct ( h 2 d ) and maternal ( h 2 m ) genetic effects, the correlation between direct and maternal genetic effects ( ρ dm ) and the ratio between the litter (C2) and maternal environmental (Me2) effects and the total phenotypic variance for the Intramuscular Fat of Longissimus thoracis et lumborum muscle using five models.
Trait h 2 d h 2 m ρ dmC2Me2
Model 10.84
[0.69, 0.99]
Model 20.44
[0.24, 0.64]
0.22
[0.08, 0.36]
−0.02
[−0.54, 0.51]
Model 30.54
[0.39, 0.69]
0.18
[0.12, 0.24]
Model 40.45
[0.28, 0.65]
0.09
[0.01, 0.19]
−0.02
[−0.63, 0.59]
0.14
[0.07, 0.21]
Model 50.46
[0.28, 0.63]
0.08
[0.01, 0.19]
0
[−0.64, 0.64]
0.13
[0.06, 0.21]
Model 1: a model fitting the direct additive genetic effects; Model 2 = Model 1 + maternal genetic effects; Model 3 = Model 1 + common litter effects; Model 4 = Model 1 + maternal genetic effects + common litter effects; Model 5 = Model 1 + maternal genetic effects + maternal environmental effects.
Table 3. Marginal posterior mean and highest posterior density region at 95% (between brackets) of the heritability for direct ( h 2 d ) and maternal genetic ( h 2 m ) effects and the ratio between the litter effects and the total phenotypic variance (C2) for the fatty acid composition of Longissimus thoracis et lumborum muscle.
Table 3. Marginal posterior mean and highest posterior density region at 95% (between brackets) of the heritability for direct ( h 2 d ) and maternal genetic ( h 2 m ) effects and the ratio between the litter effects and the total phenotypic variance (C2) for the fatty acid composition of Longissimus thoracis et lumborum muscle.
Trait h 2 d h 2 m C2
C14:00.50
[0.34, 0.67]
0.16
[0.07, 0.26]
0.16
[0.09, 0.23]
C16:00.54
[0.38, 0.71]
0.15
[0.06, 0.24]
0.14
[0.09, 0.20]
C18:00.40
[0.24, 0.56]
0.10
[0.01, 0.19]
0.14
[0.07, 0.20]
SFA0.35
[0.22, 0.47]
0.14
[0.04, 0.24]
0.21
[0.15, 0.27]
C16:1n70.45
[0.27, 0.62]
0.46
[0.38, 0.54]
0.15
[0.08, 0.22]
C18:1n90.33
[0.23, 0.43]
0.17
[0.07, 0.26]
0.05
[−0.01, 0.10]
MUFA0.35
[0.24, 0.46]
0.46
[0.33, 0.60]
0.03
[0, 0.07]
C18:2n60.46
[0.35, 0.57]
0.31
[0.18, 0.44]
0.12
[0.05, 0.18]
C18:3n30.50
[0.32, 0.67]
0.33
[0.19, 0.47]
0.18
[0.10, 0.26]
C20:4n60.38
[0.20, 0.56]
0.20
[0.05, 0.36]
0.11
[0.05, 0.18]
PUFA0.43
[0.30, 0.55]
0.17
[0.08, 0.27]
0.09
[0.02, 0.15]
Model 4 = direct additive genetic effects + maternal genetic effects + common litter effects; SFA: saturated fatty acids = C14:0 + C15:0 + C16:0 + C17:0 + C18:0; MUFA: monounsaturated fatty acids= C16:1n7 + C18:1n7 + C18:1n9; PUFA: polyunsaturated fatty acids = C18:2n6 + C18:3n3 + C20:2n6 + C20:3n6 + C20:4n6 + C20:5n3 + C22:4n6 + C22:5n3 + C22:6n3; C14:0: myristic acid; C16:0: palmitic acid; C18:0: stearic acid; C16:1n7: palmitoleic acid; C18:1n9: oleic acid; C18:2n-6: linoleic acid; C18:3n3: α-linolenic acid; C20:4n-6: arachidonic acid.
Table 4. Marginal posterior mean and highest posterior density region at 95% (between brackets) of direct and maternal genetic correlations between intramuscular fat and fatty acid composition and the genetic correlation between additive and maternal genetic effects for fatty acid composition of Longissimus thoracis et lumborum muscle.
Table 4. Marginal posterior mean and highest posterior density region at 95% (between brackets) of direct and maternal genetic correlations between intramuscular fat and fatty acid composition and the genetic correlation between additive and maternal genetic effects for fatty acid composition of Longissimus thoracis et lumborum muscle.
Trait ρ d 1 ρ m 2 ρ dm 3
C14:00.96
[0.90, 1.00]
0.68
[0.27, 1.00]
−0.21
[−0.60, 0.18]
C16:00.69
[0.46, 0.91]
−0.02
[−0.56, 0.52]
−0.66
[−0.94, −0.39]
C18:0−0.81
[−0.98, −0.64]
−0.33
[−0.81, 0.16]
−0.15
[−0.61, 0.32]
SFA0.97
[0.88, 1.00]
−0.03
[−0.55, 0.49]
−0.89
[−1.00, −0.74]
C16:1n70.92
[0.81, 1.00]
0.74
[0.50, 0.97]
−0.39
[−0.73, −0.05]
C18:1n90.36
[0.15, 0.57]
0.68
[0.44, 0.91]
−0.66
[−0.83, −0.49]
MUFA0.40
[0.19, 0.61]
0.52
[0.17, 0.88]
−0.53
[−0.69, −0.37]
C18:2n6−0.27
[−0.51, −0.03]
−0.13
[−0.57, 0.31]
−0.54
[−0.77, −0.31]
C18:3n30.52
[0.20, 0.84]
0.33
[−0.39, 1.00]
−0.38
[−0.82, 0.06]
C20:4n6−0.63
[−0.84, −0.42]
−0.81
[−1.00, −0.61]
0.15
[−0.30, 0.60]
PUFA−0.51
[−0.69, −0.32]
−0.19
[−0.58, 0.20]
−0.56
[−0.72, −0.40]
1 ρ d: genetic correlation of direct additive genetic effects between IMF and fatty acid composition; 2  ρ m: genetic correlation of maternal genetic effects between IMF and fatty acid composition; 3  ρ dm: genetic correlation between direct and maternal additive genetic effects for fatty acid composition; IMF: intramuscular fat; SFA: saturated fatty acids = C14:0 + C15:0 + C16:0 + C17:0 + C18:0; MUFA: monounsaturated fatty acids = C16:1n7 + C18:1n7 + C18:1n9; PUFA: polyunsaturated fatty acids = C18:2n6 + C18:3n3 + C20:2n6 + C20:3n6 + C20:4n6 + C20:5n3 + C22:4n6 + C22:5n3 + C22:6n3; C14:0: myristic acid; C16:0: palmitic acid; C18:0: stearic acid; C16:1n7: palmitoleic acid; C18:1n9: oleic acid; C18:2n-6: linoleic acid; C18:3n3: α-linolenic acid; C20:4n-6: arachidonic acid.
Table 5. The percentage of genomic variance explained by each relevant genomic region 1 associated with intramuscular fat and fatty acid composition of Longissimus thoracis et lumborum muscle in rabbits.
Table 5. The percentage of genomic variance explained by each relevant genomic region 1 associated with intramuscular fat and fatty acid composition of Longissimus thoracis et lumborum muscle in rabbits.
OCUMap Position (bp)Trait (% of Explained Variance 2)No. of SNPs 3
135,902,181–37,321,725PUFA/SFA (13.9%)50
135,902,181–37,337,488C16:0 (12.2%)48
135,902,181–37,337,488C16:1n7 (9.9%)43
135,902,181–37,337,488PUFA (3.2%)44
510,141,256–10,568,012C16:0 (3.5%)12
510,141,256–10,568,012SFA (5.1%)12
510,141,256–10,568,012PUFA/SFA (4.3%)12
510,161,032–10,568,012C18:2n6 (3.8%)8
1948,610,554–50,061,214IMF (8.0%)36
1949,573,782–49,993,977C20:4n6 (1.2%)21
2016,022,958–16,482,079C20:4n6 (1.6%)19
1: Genomic windows explaining ≥ 1% of the genetic variance and containing SNPs with a Bayes Factor ≥ 10; OCU: rabbit chromosome; 2: percentage of genomic variance explained by the window; IMF: intramuscular fat; C16:0: palmitic acid; SFA: saturated fatty acids; C16:1n7: palmitoleic acid; C18:2n-6: linoleic acid; C20:4n-6: arachidonic acid; PUFA: polyunsaturated fatty acids; 3: number of SNPs with a Bayes Factor ≥ 10.
Table 6. Genomic regions associated with intramuscular fat and fatty acid composition of Longissimus thoracis et lumborum muscle in rabbits.
Table 6. Genomic regions associated with intramuscular fat and fatty acid composition of Longissimus thoracis et lumborum muscle in rabbits.
TraitOCUMap Position (bp)No. of GenesPCG 1
StartEnd
IMF1948,000,58850,994,79141TANC2, CYB561, ACE, ENSOCUG00000034369, KCNH6, DCAF7, MAP3K3, LIMD2, STRADA, CCDC47, DDX42, FTSJ3, ENSOCUG00000004329, SMARCD2, ENSOCUG00000026494, GHB1, CD79B, SCN4A, ENSOCUG00000035912, ICAM2,
ENSOCUG00000029995, ENSOCUG00000004362, TEX2, ENSOCUG00000034748, PECAM1, ENSOCUG00000007896, POLG2, DDX5, CEP95, SMURF2, ENSOCUG00000023204, ENSOCUG00000015913, BPTF, NOL11, PITPNC1, PSMD12, HELZ, ENSOCUG00000008850, CACNG4, CACNG5, PRKCA
C16:0135,072,64936,992,5664ADAMTSL1, SH3GL2, CNTLN, BNC2
510,022,60710,999,5402RPGRIP1L, FTO
SFA136,002,87036,992,5662SH3GL2, CNTLN
510,022,60710,999,5402RPGRIP1L, FTO
C16:1n7135,072,64937,977,7244ADAMTSL1, SH3GL2, CNTLN, BNC2
PUFA135,072,64937,977,7244ADAMTSL1, SH3GL2, CNTLN, BNC2
510,022,60710,999,5402RPGRIP1L, FTO
PUFA/SFA135,017,58137,620,9964ADAMTSL1, SH3GL2, CNTLN, BNC2
59,928,94610,444,6692RPGRIP1L, FTO
OCU, rabbit chromosome; bp, base pair; N, number of genes; PCG, positional candidate gene (only protein-coding genes); 1 Novel genes are named according to their Ensembl gene ID and retrieved from DAVID database; IMF, intramuscular fat; C16:0, palmitic acid; SFA, saturated fatty acids; C16:1n-7 acid, palmitoleic; PUFA, polyunsaturated fatty acids.
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

EL Nagar, A.G.; Heddi, I.; Sosa-Madrid, B.S.; Blasco, A.; Hernández, P.; Ibáñez-Escriche, N. Genome-Wide Association Study of Maternal Genetic Effects on Intramuscular Fat and Fatty Acid Composition in Rabbits. Animals 2023, 13, 3071. https://doi.org/10.3390/ani13193071

AMA Style

EL Nagar AG, Heddi I, Sosa-Madrid BS, Blasco A, Hernández P, Ibáñez-Escriche N. Genome-Wide Association Study of Maternal Genetic Effects on Intramuscular Fat and Fatty Acid Composition in Rabbits. Animals. 2023; 13(19):3071. https://doi.org/10.3390/ani13193071

Chicago/Turabian Style

EL Nagar, Ayman G., Imen Heddi, Bolívar Samuel Sosa-Madrid, Agustín Blasco, Pilar Hernández, and Noelia Ibáñez-Escriche. 2023. "Genome-Wide Association Study of Maternal Genetic Effects on Intramuscular Fat and Fatty Acid Composition in Rabbits" Animals 13, no. 19: 3071. https://doi.org/10.3390/ani13193071

APA Style

EL Nagar, A. G., Heddi, I., Sosa-Madrid, B. S., Blasco, A., Hernández, P., & Ibáñez-Escriche, N. (2023). Genome-Wide Association Study of Maternal Genetic Effects on Intramuscular Fat and Fatty Acid Composition in Rabbits. Animals, 13(19), 3071. https://doi.org/10.3390/ani13193071

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop