Antagonistic maternal and direct effects of the leptin receptor gene on body weight in pigs

Evaluating regional heritability mapping methods for identifying QTLs in a wild population of Soay sheep

The study of complex traits and their genetic underpinnings is crucial for understanding the evolutionary processes and mechanisms that shape natural populations. Regional heritability mapping (RHM) is a method for estimating the heritability of genomic segments that may contain both common and rare variants affecting a complex trait. This research is important because it advances our ability to detect genetic loci that contribute to phenotypic variation, even those that might be missed by traditional methods such as genome-wide association studies (GWAS). Here, we compare three RHM methods: SNP-RHM, which uses genomic relationship matrices (GRMs) based on SNP genotypes; Hap-RHM, which utilizes GRMs based on haplotypes; and SNHap-RHM, which integrates both SNP-based and haplotype-based GRMs jointly. These methods were applied to data from a wild population of sheep, focusing on the analysis of eleven polygenic traits. The results were compared with findings from previous GWAS to assess how RHM performed at identifying both known and novel associated loci. We found that while the inclusion of the regional matrix did not account for significant variation in all regions associated with trait variation as identified by GWAS, it did uncover several regions that were not previously linked to trait variation. This suggests that RHM methods can provide additional insights into the genetic architecture of complex traits, highlighting regions of the genome that may be overlooked by GWAS alone. This study underscores the importance of using complementary approaches to fully understand the genetic basis of complex traits in natural populations.

Association between mitochondrial DNA copy number and production traits in pigs

Mitochondria are essential organelles in the regulation of cellular energetic metabolism. Mitochondrial DNA copy number (mtDNA_CN) can be used as a proxy for mitochondria number, size, and activity. The aims of our study are to evaluate the effect of mtDNA_CN and mitochondrial haploblocks on production traits in pigs, and to identify the genetic background of this cellular phenotype. We collected performance data of 234 pigs and extracted DNA from skeletal muscle. Whole-genome sequencing data was used to determine mtDNA_CN. We found positive correlations of muscle mtDNA_CN with backfat thickness at 207 d (+0.14; p-value = 0.07) and negative correlations with carcase loin thickness (-0.14; p-value = 0.03). Pigs with mtDNA_CN values below the lower quartile had greater loin thickness (+4.1 mm; p-value = 0.01) and lower backfat thickness (-1.1 mm; p-value = 0.08), which resulted in greater carcase lean percentage (+2.4%; p-value = 0.04), than pigs with mtDNA_CN values above the upper quartile. These results support the hypothesis that a reduction of mitochondrial activity is associated with greater feed efficiency. Higher mtDNA_CN was also positively correlated with higher meat ultimate pH (+0.19; p-value <0.01) but we did not observe significant difference for meat ultimate pH between the two groups with extreme mtDNA_CN. We found no association of the most frequent mitochondrial haploblocks with mtDNA_CN or the production traits, but several genomic regions that harbour potential candidate genes with functions related to mitochondrial biogenesis and homeostasis were associated with mtDNA_CN. These regions provide new insights into the genetic background of this cellular phenotype but it is still uncertain if such associations translate into noticeable effects on the production traits.

Maternal effects and its importance in the genetic evaluations of preweaning live weight traits of beef cattle. A review

Maternal effects in cattle genetics are defined as the causal influence of the phenotype or maternal genotype on the offspring's phenotype by effects occurring when the genetic and environmental characteristics of the mother influence the phenotype of the offspring beyond the direct inheritance of genes. Its relevance has been strongly described in genetic models focused on the genetic improvement of preweaning traits in cow-calf beef cattle production systems. Here, basic concepts and the importance of maternal effects when using linear and animal model procedures for genetic evaluations of growth and live-weight traits in beef cattle are reviewed and discussed. A brief history of estimation methods from classical studies to recent studies used for the development of animal models for studying maternal effects is also provided. Some important biometric concepts for maternal effect estimation are described, and the antagonism between direct genetic effects and maternal effects, its biological basis, and sources of error in the estimation of direct genetic and maternal covariance are discussed. Finally, some genomic perspectives are presented.

Candidate markers for enhanced host response to PRRS have scarce adverse effects on pigs' growth and production

BACKGROUND

Porcine Reproductive and Respiratory Syndrome (PRRS) is one of the most challenging viral diseases that cause substantial economic losses in the pig industry worldwide. The clinical signs of PRRS depend on, among others, the immunomodulatory properties of the PRRS virus strain, farm health status, herd immunity, and host genetics. The high virulence and mutation rate of PRRS virus limit the efficacy of vaccination programs. In recent years, several candidate genetic markers associated with PRRS resilience have been identified, and selective breeding was suggested as an additional approach to control PRRS under field conditions. Even so, it is essential to investigate the effects of these genetic markers on pigs' productivity. Our study aimed to assess the association between seven previously reported candidate genetic markers for host response to PRRS (rs80800372 in GBP1, rs340943904 in GBP5, rs322187731 in GBP6, rs1107556229 in CD163, rs338508371 in SGK1, rs80928141 in TAP1, and a 275-bp insertion in the promoter of MX1) and production traits in pigs under non-challenging conditions.

RESULTS

About 600 high-health Duroc pigs were genotyped for the selected genetic markers and their effects on production traits (live body weight, carcass weight, backfat thickness, intramuscular fat content and composition) were assessed using a linear model. The genetic markers GBP5_rs340943904, GBP6_rs322187731, CD163_rs1107556229, and the 275-bp insertion at the promoter of MX1 showed no relevant associations with growth and carcass traits at slaughter. Regarding GBP1_rs80800372 (WUR1000125), the favourable G allele for PRRS resilience displayed significant additive effects on backfat thickness (+ 1.18 ± 0.42 mm; p = 0.005) and lean content (-1.72 ± 0.56%; p ≤ 0.01) at slaughter. In addition, the genetic markers SGK1_rs338508371 and TAP1_rs8092814 were associated with the palmitoleic content in gluteus medius, without affecting the total of the monounsaturated fatty acids.

CONCLUSIONS

Our results indicate that genetic markers for PRRS resilience have no relevant effects on growth and carcass traits in pigs reared under non-challenging conditions, except for GBP1_rs80800372 where the favourable allele for PRRS response has a negative impact on lean content. Therefore, since the effects of GBP1_rs80800372 were attributed to the causal variant GBP5_rs340943904, it seems beneficial to select pigs for the genetic marker at GBP5 instead of GBP1. Overall, pigs might be selected for enhanced PRRS resilience without compromising their overall productivity.

Impact of the leptin receptor gene on pig performance and quality traits

The recessive T allele of the missense polymorphism rs709596309 C > T of the leptin receptor gene is associated with intramuscular fat. However, its overall impact on pork production is still partial. In this work, we investigated the all-round effects of the TT genotype on lean growth efficiency and carcass, meat and fat quality using data from an experiment that compared the performance of 48 TT and 48 C- (24 CT and 24 CC) Duroc barrows. The TT pigs were less efficient for lean growth than the C- pigs. Although heavier, their carcasses had less lean content, were shorter and had lighter loins. Apart from increasing marbling and saturated fatty acid content, changes caused by the TT genotype in meat and fat quality are likely not enough to be perceived by consumers. The effect on visual marbling score exceeded that on intramuscular fat content, which suggests a direct influence of the T allele on the pattern of fat distribution in muscle. With current low-protein diets, the T allele is expected to be cost-effective only in niche markets where a very high level of marbling is critical.

Identification of breed-specific genomic variants in Colombian Creole pig breeds by whole-genome sequencing

Dissecting genetic variation of local breeds is important for the success of conservation. In this research, we investigated the genomic variation of Colombian Creole (CR) pigs, with a focus on the breed-specific variants in the exonic region of 34 genes with reported effects on adaptive and economic traits. Seven individuals of each of the three CR breeds (CM, Casco de Mula; SP, San Pedreño; and ZU, Zungo) were whole-genome sequenced along with 7 Iberian (IB) pigs and 7 pigs of each of the four most used cosmopolitan (CP) breeds (Duroc, Landrace × Large White, and Pietrain). Molecular variability in CR (6,451,218 variants; from 3,919,242, in SP, to 4,648,069, in CM) was comparable to that in CP, but higher than in IB. For the investigated genes, SP pigs displayed less exonic variants (178) than ZU (254), CM (263), IB (200), and the individual CP genetic types (201 to 335). Sequence variation in these genes confirmed the resemblance of CR to IB and indicates that CR pigs, particularly ZU and CM, are not exempt from selective introgression of other breeds. A total of 50 exonic variants were identified as being potentially specific to CR, including a high-impact deletion in the intron between exons 15 and 16 of the leptin receptor gene, which was only found in CM and ZU. The identification of breed-specific variants in genes related to adaptive and economical traits can bolster the understanding of the role of gene-environment interactions on local adaptation and points the way for effective breeding and conservation of CR pigs.

The leptin receptor gene affects piglet behavior and growth

Piglets with low birth weight present low vitality after farrowing, often leading to impaired weight gain during lactation. A recessive missense variant (C > T) for increased appetite and fatness in the porcine leptin receptor gene (rs709596309) causes a negative maternal effect on the weight of piglets at weaning. However, it is not known whether this variant already exerts an effect on the birth weight and vitality of newborn piglets and on their growing capacity during lactation. An experiment was conducted using 668 purebred Duroc piglets (131 CC, 311 CT, and 226 TT) from 74 multiparous sows (9 CC, 43 CT, and 22 TT) and 14 boars (1 CC, 10 CT, and 3 TT). All piglets were individually weighed at birth and tested for vitality, which was assessed on a scale from 1 (low vitality) to 3 (high vitality) based on behavioral observations, including the status of the piglet immediately before the test. Only non-adopted piglets were considered for piglet performance at weaning. Inferences on the effect of the genotype on birth and weaning traits were done on a Bayesian setting based on 2-trait bivariate models including the effects of the piglet and the litter, as well as the genotype of the sow and the piglet, the sex of the piglet, and the parity number. Vitality and the status of the piglet before the test were analyzed using a liability threshold (probit) model. As compared to other genotypes, TT newborn piglets were 28 g heavier, were more vital (the probability of being scored as highly vital was 6.5% higher) and were more often found suckling before the test (the probability of being suckling at test was 6.5% higher). As a result, they grew more during lactation (153 g) and were heavier at weaning (169 g) than littermates of the two other genotypes, thus partly compensating for the limited maternal capacity of TT sows. Our findings provide evidence that appetite-influencing genes, such as the leptin receptor gene, have developmental implications from very early life stages.

Genome-wide analyses reveal a strong association between LEPR gene variants and body fat reserves in ewes

BACKGROUND

Among the adaptive capacities of animals, the management of energetic body reserves (BR) through the BR mobilization and accretion processes (BR dynamics, BRD) has become an increasingly valuable attribute for livestock sustainability, allowing animals to cope with more variable environments. BRD has previously been reported to be heritable in ruminants. In the present work, we conducted genome-wide studies (GWAS) in sheep to determine genetic variants associated with BRD. BR (i.e. levels) and BRD (i.e. changes over time) were obtained through body condition score measurements at eight physiological stages throughout each productive cycle in Romane ewes (n = 1034) and were used as phenotypes for GWAS. After quality controls and imputation, 48,593 single nucleotide polymorphisms (SNP) were included in the GWAS.

RESULTS

Among the 23 QTL regions identified, a major QTL associated with BR during pregnancy and lactation was identified on chromosome 1. In this region, several significant SNPs mapped to the leptin receptor gene (LEPR), among which one SNP mapped to the coding sequence. The point mutation induces the p.P1019S substitution in the cytoplasmic domain, close to tyrosine phosphorylation sites. The frequency of the SNP associated with increased BR levels was 32%, and the LEPR genotype explained up to 5% of the variance of the trait. Higher fatness levels in ewes carrying the LEPR p.P1019S mutation were observed all along the productive cycle.

CONCLUSIONS

These results provide strong evidences for involvement of LEPR in the regulation of BR in sheep and highlight it as a major candidate gene for improving adaptive capacities by genetic selection.

Heritability and genetic correlation estimates of semen production traits with litter traits and pork production traits in purebred Duroc pigs

We estimated heritabilities of semen production traits and their genetic correlations with litter traits and pork production traits in purebred Duroc pigs. Semen production traits were semen volume, sperm concentration, proportion of morphologically normal sperms, total number of sperm, and total number of morphologically normal sperm. Litter traits at farrowing were total number born, number born alive, number stillborn, total litter weight at birth, mean litter weight at birth, and piglet survival rate at birth. Litter traits at weaning were litter size at weaning, total litter weight at weaning, mean litter weight at weaning, and piglet survival rate from birth to weaning. Pork production traits were average daily gain, backfat thickness, and loin muscle area. We analyzed 45,913 semen collection records of 896 boars, 6,950 farrowing performance records of 1,400 sows, 2,237 weaning performance records of 586 sows, and individual growth performance records of 9,550 animals measured at approximately 5 mo of age. Heritabilities were estimated using a single-trait animal model. Genetic correlations were estimated using a 2-trait animal model. Estimated heritabilities of semen production traits ranged from 0.20 for sperm concentration to 0.29 for semen volume and were equal to or higher than those of litter traits, ranging from 0.06 for number stillborn and piglet survival rate at birth to 0.25 for mean litter weight at birth, but lower than those of pork production traits, ranging from 0.50 for average daily gain to 0.63 for backfat thickness. In many cases, the absolute values of estimated genetic correlations between semen production traits and other traits were smaller than 0.3. These estimated genetic parameters provide useful information for establishing a comprehensive pig breeding scheme.

Genetic parameters of 5 semen production traits, 10 litter traits, and 3 pork production traits in purebred Duroc pigs was estimated. Heritabilities of semen production traits ranged from 0.20 for sperm concentration to 0.29 for semen volume and were equal to or higher than those of litter traits, ranging from 0.06 for number stillborn and piglet survival rate at birth to 0.25 for mean litter weight at birth, but lower than those of pork production traits, ranging from 0.50 for average daily gain to 0.63 for backfat thickness. In many cases, the absolute values of genetic correlations between semen production traits and other traits were smaller than 0.3. These estimated genetic parameters provide useful information for establishing a comprehensive pig breeding scheme.

Genetic architecture and major genes for backfat thickness in pig lines of diverse genetic backgrounds

Background

Backfat thickness is an important carcass composition trait for pork production and is commonly included in swine breeding programmes. In this paper, we report the results of a large genome-wide association study for backfat thickness using data from eight lines of diverse genetic backgrounds.

Methods

Data comprised 275,590 pigs from eight lines with diverse genetic backgrounds (breeds included Large White, Landrace, Pietrain, Hampshire, Duroc, and synthetic lines) genotyped and imputed for 71,324 single-nucleotide polymorphisms (SNPs). For each line, we estimated SNP associations using a univariate linear mixed model that accounted for genomic relationships. SNPs with significant associations were identified using a threshold of p < 10^–6 and used to define genomic regions of interest. The proportion of genetic variance explained by a genomic region was estimated using a ridge regression model.

Results

We found significant associations with backfat thickness for 264 SNPs across 27 genomic regions. Six genomic regions were detected in three or more lines. The average estimate of the SNP-based heritability was 0.48, with estimates by line ranging from 0.30 to 0.58. The genomic regions jointly explained from 3.2 to 19.5% of the additive genetic variance of backfat thickness within a line. Individual genomic regions explained up to 8.0% of the additive genetic variance of backfat thickness within a line. Some of these 27 genomic regions also explained up to 1.6% of the additive genetic variance in lines for which the genomic region was not statistically significant. We identified 64 candidate genes with annotated functions that can be related to fat metabolism, including well-studied genes such as MC4R, IGF2, and LEPR, and more novel candidate genes such as DHCR7, FGF23, MEDAG, DGKI, and PTN.

Conclusions

Our results confirm the polygenic architecture of backfat thickness and the role of genes involved in energy homeostasis, adipogenesis, fatty acid metabolism, and insulin signalling pathways for fat deposition in pigs. The results also suggest that several less well-understood metabolic pathways contribute to backfat development, such as those of phosphate, calcium, and vitamin D homeostasis.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12711-021-00671-w.