Genetic and heat-stress related environmental influences on pig whole-blood gene expression levels

This study investigates the relative impacts of genetics and heat-related environmental factors on pig whole-blood gene expression, identifying thousands of differentially expressed genes and eQTLs while highlighting specific genetic mechanisms and candidate genes involved in thermoregulation and production traits.

The Gist

Imagine a massive, high-stakes cooking competition where the goal isn't to make the best dish, but to see which pigs can survive a scorching kitchen without melting. This is essentially what scientists did in this study, but instead of ovens, they used two different pig farms: one in the humid, tropical heat of Guadeloupe (the "Tropical Kitchen") and one in the cooler, temperate climate of France (the "Temperate Kitchen").

Here is the story of what they discovered, broken down into simple concepts.

1. The Setup: Nature vs. Nurture (and Heat)

The researchers took a group of pigs that were all closely related cousins (a mix of hardy tropical "Créole" pigs and fast-growing temperate "Large White" pigs). They split them up:

• Group A lived in the hot, tropical farm, getting used to the heat day after day.
• Group B lived in the cool farm, but at the end of the experiment, they were forced into a "heat stress challenge" (a 3-week sauna) to see how they reacted.

The scientists wanted to know: Is it the pigs' DNA (Nature) or the weather (Nurture) that changes how their bodies work? To find out, they looked at the "instruction manuals" inside the pigs' blood cells (gene expression).

2. The Big Discovery: The "Tropical Team" vs. The "Temperate Team"

When they compared the blood of the two groups, they found 1,967 genes that were behaving differently.

• The Tropical Pigs: Their bodies were in "survival mode." Their genes were busy fighting off stress, managing their immune systems, and dealing with the constant heat. Think of them as soldiers constantly on high alert, checking their armor and repairing damage.
• The Temperate Pigs: Their genes were focused on "growth mode." They were busy processing food and building muscle. Think of them as construction workers calmly building a skyscraper.

The Analogy: Imagine two cars. One is driving through a dusty, rocky desert (Tropical), so its engine is running hot and the suspension is working overtime. The other is driving on a smooth highway (Temperate), so it's focused on speed and fuel efficiency. The study showed that the "desert car" had to change its engine settings just to keep running.

3. The "Heat Shock" Experiment

When the cool-weather pigs were suddenly put into the heat for 3 weeks, their bodies went into panic mode.

• Day 3: Their bodies started scrambling to fix damaged proteins (like a mechanic trying to patch a leaky hose).
• Day 18: After three weeks, their immune systems were exhausted. They started showing signs of inflammation and were less able to fight off potential infections. It was like a marathon runner who started strong but was completely worn out by the finish line.

4. The Genetic "Remote Controls" (eQTLs)

This is the most technical part, but here's the simple version:
The scientists found 6,014 "remote controls" (called eQTLs) in the pigs' DNA that could turn specific genes on or off.

• Cis-eQTLs (The Local Remote): These are remotes that control the light switch right next to them. Most of the controls found were local, meaning a small change in DNA right next to a gene changed how that gene worked.
• Trans-eQTLs (The Master Remote): These are rare, powerful remotes that can control lights in other rooms. The scientists found a "Master Remote" on a specific chromosome (SSC12) controlled by a gene called GPATCH8. This one remote seemed to influence dozens of other genes related to blood cell production and immunity. It's like finding a single switch in the basement that controls the lights in the entire house.

5. The "Who's Who" of Pig Traits

The study also played a game of "connect the dots" to see which genes were responsible for specific pig traits, like how much fat they have (backfat) or how hot they get (thermoregulation).

• They found that two specific genes, TMCO1 and ZNF184, were the "architects" behind how much fat a pig carries.
• This is huge for farmers because if you know which genes control fat, you can breed pigs that are healthier and more efficient, even in hot weather.

6. Why Does This Matter?

Pigs are sensitive to heat. When it gets too hot, they stop eating, stop growing, and get sick. This costs the farming industry millions of dollars every year.

The Takeaway:
This study is like a user manual for the pig's body in the heat. It tells us:

1. Heat changes the body: It shifts the body from "growing" to "surviving."
2. Genetics matter: Some pigs are naturally better at handling the heat because of their DNA "remote controls."
3. We can breed better pigs: By identifying the specific genes (like GPATCH8, TMCO1, and ZNF184), scientists can help farmers breed pigs that stay cool, eat well, and grow fast, even when the temperature rises.

In short, the scientists didn't just watch pigs sweat; they decoded the secret language of their blood to figure out how to make them tougher against the heat.

Technical Summary

1. Problem Statement

Heat stress is a major constraint in the global pig industry, causing significant economic losses (estimated at US$300 million annually in the US) and compromising animal welfare, growth, and health. While phenotypic indicators of heat resilience (e.g., rectal temperature, feed intake) are well-documented, the molecular mechanisms underlying genetic variation in heat adaptation remain poorly understood. Specifically, there is a lack of comprehensive data quantifying the relative influence of genetics versus environmental factors (chronic heat vs. acute heat stress) on the pig transcriptome, and a scarcity of identified expression Quantitative Trait Loci (eQTLs) associated with heat stress responses.

2. Methodology

The study utilized a unique backcross population (1/4 Créole, 3/4 Large White) designed to study heat adaptation. The experimental design involved two distinct environments and multiple analytical layers:

• Animal Model & Experimental Design:
  • Population: 359 backcross pigs (181 from a tropical facility in Guadeloupe; 180 from a temperate facility in France).
  • Chronic Environment: Tropical pigs were raised in a naturally heat-stressed environment. Temperate pigs were raised in a standard environment but subjected to an experimental acute heat stress challenge (30°C) for 3 weeks at the end of their growth period.
  • Sampling: Whole blood was collected at 23 weeks (d0). Temperate pigs were sampled again at day 3 (d3) and day 18 (d18) of the heat stress challenge.
• Omics Data Generation:
  • Transcriptomics: Whole-blood RNA was analyzed using Agilent SurePrint G3 microarrays (8x60K). Data underwent rigorous filtering, normalization, and deduplication.
  • Genomics: Animals were genotyped using the Illumina PorcineSNP60 BeadChip. After quality control, 44,506 SNPs were retained.
• Statistical Analyses:
  • Differential Expression (DEG): Used limma to identify genes differentially expressed between facilities (Tropical vs. Temperate) and across timepoints (d0, d3, d18).
  • Expression QTL (eQTL) Mapping: Performed Genome-Wide Association Studies (GWAS) on gene expression levels using linear mixed models (GEMMA) to identify SNPs regulating gene expression. Defined cis- and trans-eQTLs based on genomic distance and linkage disequilibrium (LD).
  • Transcriptome-Wide Association Study (TWAS): Correlated gene expression levels with phenotypic traits (growth, backfat thickness, thermoregulation).
  • Colocalization: Integrated TWAS, eQTL, and previously published QTL data to pinpoint candidate genes.
  • Hotspot Analysis: Identified trans-eQTL hotspots (genomic regions regulating many distal genes) and candidate master regulators.

3. Key Contributions

• Quantification of Variance: Provided a precise estimation of the relative contribution of genetics vs. environment to gene expression variance in pigs.
• Comprehensive eQTL Catalog: Generated a large-scale catalog of 6,014 eQTLs in pig whole blood, distinguishing between cis and trans effects in a heat-adaptation context.
• Discovery of Regulatory Hotspots: Identified specific trans-eQTL hotspots and proposed candidate master regulator genes (e.g., GPATCH8) driving hematopoietic mechanisms.
• Multi-Omics Integration: Successfully integrated TWAS, eQTL, and QTL data to fine-map genomic regions associated with complex traits like backfat thickness and thermoregulation.
• Dataset Availability: Made a significant dataset of 718 transcriptomic samples and genotypic data publicly available, serving as a baseline for future farm animal eQTL studies.

4. Key Results

A. Differential Gene Expression (DGE)

• Chronic Environment: Identified 1,967 Differentially Expressed Genes (DEGs) between tropical and temperate facilities.
  • Tropical: Overexpression of immune-related pathways, prolactin signaling, and macromolecule metabolic processes.
  • Temperate: Overexpression of RNA/mRNA processing mechanisms.
• Acute Heat Stress: Identified 472 DEGs during the 3-week heat stress challenge.
  • Early Response (d3): Downregulation of chaperone-mediated autophagy (HSPA8); upregulation of methylation and syndecan interaction genes.
  • Late Response (d18): Upregulation of immune system pathways (neutrophil degranulation, NOD-like receptor signaling) and downregulation of Vitamin C metabolism.

B. Genetic Control of Expression (eQTLs)

• Variance Explained: Genetic factors explained an average of 36.3% of gene expression variance, acting as the primary source for 27.7% of transcripts.
• eQTL Landscape: Identified 6,014 eQTLs affecting 3,297 genes.
  • Cis-eQTLs: 4,222 (70.2% of total), affecting 2,887 genes. These had larger effect sizes.
  • Trans-eQTLs: 995 (16.5% of total), affecting 695 genes. These explained a higher percentage of variance (62% vs. 60% for cis) but had smaller effect sizes.
• Hotspots: Identified 38 trans-eQTL hotspots. Two major hotspots on Chromosome 12 (SSC12) were highlighted:
  • One linked to candidate genes GPATCH8 and CWC25, associated with platelet functions and early erythroid development.
  • Another linked to ANKRD40, CDK12, and others, associated with hemorrhage-related functions.

C. Genotype-by-Environment (G×E) Interactions

• Detected 73 SNPs with significant G×E interactions, primarily affecting 9 genes via cis-eQTLs. However, the study noted low statistical power to reliably detect these interactions given the sample size.

D. TWAS and Colocalization

• TWAS: Identified 139 genes associated with production and thermoregulation traits. Notable associations included FUNDC1, PLGRKT, and PSMA3 (negatively correlated with growth) and ING2, PER1 (positively correlated).
• Candidate Genes for Backfat Thickness: Integration of TWAS, eQTL, and QTL data pointed to TMCO1 and ZNF184 as strong candidates for regulating backfat thickness. ZNF184 was specifically linked to the activation of the FTO gene (obesity regulation).

5. Significance and Conclusion

This study provides a foundational resource for understanding the molecular basis of heat resilience in pigs.

• Biological Insight: It reveals that while chronic heat stress alters immune and metabolic pathways, acute stress triggers specific inflammatory and hematopoietic responses.
• Genetic Architecture: It demonstrates that genetic variation is a dominant driver of transcriptomic variance, often outweighing environmental effects, and that trans-regulatory networks (hotspots) play a crucial role in coordinating complex immune responses.
• Breeding Applications: The identification of candidate genes (TMCO1, ZNF184, GPATCH8) and the quantification of heritability for expression traits offer new targets for genomic selection to breed pigs with improved heat tolerance and production efficiency.
• Limitations: The study acknowledges limitations regarding the microarray technology (predefined probes) and the statistical power to detect subtle G×E interactions, suggesting future work should utilize RNA-seq and larger cohorts.

In summary, the paper successfully bridges the gap between genotype, environment, and phenotype in livestock, offering a roadmap for improving heat resilience in the face of climate change.