Changes in the N-glycosylation of porcine immune globulin G during postnatal development

N-glycosylation influences the effectiveness of immune globulin G (IgG) and thus the immunological downstream responses of immune cells. This impact arises from the presence of N-glycans within the Fc region, which not only alters the conformation of IgG but also influences its steric hindrance. Consequently, these modifications affect the interaction between IgG and its binding partners within the immune system. Moreover, this posttranslational modification vary according to the physiological condition of each individual. In this study, we examined the N-glycosylation of IgG in pigs from birth to five months of age. Our analysis identified a total of 48 distinct N-glycan structures. Remarkably, we observed defined changes in the composition of these N-glycans during postnatal development. The presence of agalactosylated and sialylated structures increases in relation to the number of N-glycans terminated by galactose residues during the first months of life. This shift may indicate a transition from passively transferred antibodies from the colostrum of the sow to the active production of endogenous IgG by the pig's own immune system.

Bile acid profiles and mRNA abundance of bile acid-related genes in adipose tissue of dairy cows with high versus normal body condition

Besides their lipid-digestive role, bile acids (BA) influence overall energy homeostasis, such as glucose and lipid metabolism. We hypothesized that BA along with their receptors, regulatory enzymes, and transporters are present in subcutaneous adipose tissue (scAT). In addition, we hypothesized that their mRNA abundance varies with the body condition of dairy cows around calving. Therefore, we analyzed BA in serum and scAT as well as the mRNA abundance of BA -related enzymes, transporters, and receptors in scAT during the transition period in cows with different body conditions around calving. In a previously established animal model, 38 German Holstein cows were divided into either a high (HBCS; n = 19) or normal BCS (NBCS; n = 19) group based on their body condition score (BCS) and back fat thickness (BFT). Cows were fed different diets to achieve the targeted differences in BCS and BFT (NBCS: BCS <3.5, BFT <1.2 cm; HBCS: BCS >3.75, BFT >1.4 cm) until dry-off at 7 wk ante partum. During the dry period and subsequent lactation, both groups were fed the same diets regarding their demands. Using a targeted metabolomics approach via LC-ESI-MS /MS, BA were analyzed in serum and scAT at wk -7, 1, 3, and 12 relative to parturition. In serum, 15 BA (cholic acid (CA), chenodeoxycholic acid (CDCA), glycocholic acid (GCA), taurocholic acid (TCA), glycochenodeoxycholic acid (GCDCA), taurochenodeoxycholic acid (TCDCA), deoxycholic acid (DCA), lithocholic acid (LCA), glycodeoxycholic acid (GDCA), glycolithocholic acid (GLCA), taurodeoxycholic acid (TDCA), taurolithocholic acid (TLCA), β-muricholic acid (MCA(b)), tauromuricholic acid (sum of α and β) (TMCA (a+b)), glycoursodeoxycholic acid (GUDCA)) were observed, whereas in scAT 7 BA (CA, GCA, TCA, GCDCA, TCDCA, GDCA, TDCA) were detected. In serum and scAT samples, the primary BA CA and its conjugate GCA were predominantly detected. Increasing serum concentrations of CA, CDCA, TCA, GCA, GCDCA, DCA, and MCA(b) with the onset of lactation might be related to the increasing DMI after parturition. Furthermore, serum concentrations of CA, CDCA, GCA, DCA, GCDCA, TCA, LCA, and GDCA were lower in HBCS cows compared with NBCS cows, concomitant with increased lipolysis in HBCS cows. The correlation between CA in serum and scAT may point to the transport of CA across cell membranes. Overall, the findings of the present study suggest a potential role of BA in lipid metabolism depending on the body condition of periparturient dairy cows.

RNA-Seq-based discovery of genetic variants and allele-specific expression of two layer lines and broiler chicken

Recent advances in the selective breeding of broilers and layers have made poultry production one of the fastest-growing industries. In this study, a transcriptome variant calling approach from RNA-seq data was used to determine population diversity between broilers and layers. In total, 200 individuals were analyzed from three different chicken populations (Lohmann Brown (LB), n = 90), Lohmann Selected Leghorn (LSL, n = 89), and Broiler (BR, n = 21). The raw RNA-sequencing reads were pre-processed, quality control checked, mapped to the reference genome, and made compatible with Genome Analysis ToolKit for variant detection. Subsequently, pairwise fixation index (F ST) analysis was performed between broilers and layers. Numerous candidate genes were identified, that were associated with growth, development, metabolism, immunity, and other economically significant traits. Finally, allele-specific expression (ASE) analysis was performed in the gut mucosa of LB and LSL strains at 10, 16, 24, 30, and 60 weeks of age. At different ages, the two-layer strains showed significantly different allele-specific expressions in the gut mucosa, and changes in allelic imbalance were observed across the entire lifespan. Most ASE genes are involved in energy metabolism, including sirtuin signaling pathways, oxidative phosphorylation, and mitochondrial dysfunction. A high number of ASE genes were found during the peak of laying, which were particularly enriched in cholesterol biosynthesis. These findings indicate that genetic architecture as well as biological processes driving particular demands relate to metabolic and nutritional requirements during the laying period shape allelic heterogeneity. These processes are considerably affected by breeding and management, whereby elucidating allele-specific gene regulation is an essential step towards deciphering the genotype to phenotype map or functional diversity between the chicken populations. Additionally, we observed that several genes showing significant allelic imbalance also colocalized with the top 1% of genes identified by the FST approach, suggesting a fixation of genes in cis-regulatory elements.

Exploration of genotype-by-environment interactions affecting gene expression responses in porcine immune cells

As one of the keys to healthy performance, robustness of farm animals is gaining importance, and with this comes increasing interest in genetic dissection of genotype-by-environment interactions (G×E). Changes in gene expression are among the most sensitive responses conveying adaptation to environmental stimuli. Environmentally responsive regulatory variation thus likely plays a central role in G×E. In the present study, we set out to detect action of environmentally responsive cis-regulatory variation by the analysis of condition-dependent allele specific expression (cd-ASE) in porcine immune cells. For this, we harnessed mRNA-sequencing data of peripheral blood mononuclear cells (PBMCs) stimulated in vitro with lipopolysaccharide, dexamethasone, or their combination. These treatments mimic common challenges such as bacterial infection or stress, and induce vast transcriptome changes. About two thirds of the examined loci showed significant ASE in at least one treatment, and out of those about ten percent exhibited cd-ASE. Most of the ASE variants were not yet reported in the PigGTEx Atlas. Genes showing cd-ASE were enriched in cytokine signaling in immune system and include several key candidates for animal health. In contrast, genes showing no ASE featured cell-cycle related functions. We confirmed LPS-dependent ASE for one of the top candidates, SOD2, which ranks among the major response genes in LPS-stimulated monocytes. The results of the present study demonstrate the potential of in vitro cell models coupled with cd-ASE analysis for the investigation of G×E in farm animals. The identified loci may benefit efforts to unravel the genetic basis of robustness and improvement of health and welfare in pigs.

Glucocorticoid receptor hypersensitivity enhances inflammatory signaling and inhibits cell cycle progression in porcine PBMCs

The consequences of glucocorticoid receptor (GR) hypersensitivity during infection have so far received little attention. We previously discovered that a natural gain-of-function Ala610Val substitution in the porcine GR aggravates response of pigs to lipopolysaccharide (LPS)-induced endotoxemia, which can be alleviated by dexamethasone (DEX) pretreatment. In this work, we investigated the relevant molecular basis of these phenotypes by transcriptomic profiling of porcine peripheral blood mononuclear cells (PBMCs) carrying different GR genotypes, in unstimulated conditions or in response to DEX and/or LPS in vitro. The Val allele differentially regulated abunda+nt genes in an additive-genetic manner. A subset of more than 200 genes was consistently affected by the substitution across treatments. This was associated with upregulation of genes related i.a. to endo-lysosomal system, lipid and protein catabolism, and immune terms including platelet activation, and antigen presentation, while downregulated genes were mainly involved in cell cycle regulation. Most importantly, the set of genes constitutively upregulated by Val includes members of the TLR4/LPS signaling pathway, such as LY96. Consequently, when exposing PBMCs to LPS treatment, the Val variant upregulated a panel of additional genes related to TLR4 and several other pattern recognition receptors, as well as cell death and lymphocyte signaling, ultimately amplifying the inflammatory responses. In contrast, when stimulated by DEX treatment, the Val allele orchestrated several genes involved in anti-inflammatory responses during infection. This study provides novel insights into the impact of GR hypersensitivity on the fate and function of immune cells, which may be useful for endotoxemia therapy.

Regulation of Porcine Oviduct Epithelium Functions via Progesterone and Estradiol Is Influenced by Cortisol

Preimplantation maternal stress, characterized by elevated glucocorticoids (GCs), has been linked to reproductive failures caused by impaired oviduct functionality, which is known to be predominantly regulated by the sex steroids, progesterone (P4) and (17)estradiol (E2). Although steroid receptors share analogous structures and binding preferences, the interaction between GCs and E2/P4 in the oviduct has attracted little attention. Using an air-liquid interface culture model, porcine oviduct epithelial cells were stimulated with single (cortisol, E2, P4) or hormone mixtures (cortisol/E2, cortisol/P4) for 12 hours and 72 hours. Cultures were subsequently assessed for epithelial morphometry, bioelectrical properties, and gene expression responses (steroid hormone signaling, oviductal function, immune response, and apoptosis). Results confirmed the suppressive role of P4 in regulating oviduct epithelium characteristics, which was partially opposed by E2. Besides increasing the ratio of ciliated cells, cortisol antagonized the effect of P4 on epithelial polarity and modified sex steroid-induced changes in transepithelial electrical properties. Both sex steroids affected the glucocorticoid receptor expression, while cortisol downregulated the expression of progesterone receptor. The overall gene expression pattern suggests that sex steroid dominates the cotreatment, but cortisol contributes by altering the gene responses to sex steroids. We conclude that besides its individual action, maternal cortisol interplays with sex steroids at phenotypic and molecular levels in the oviduct epithelium, thereby influencing the microenvironment of gametes and early embryos.

SNP discovery and association study for growth, fatness and meat quality traits in Iberian crossbred pigs

Iberian pigs and its crosses are produced to obtain high-quality meat products. The objective of this work was to evaluate a wide panel of DNA markers, selected by biological and functional criteria, for association with traits related to muscle growth, fatness, meat quality and metabolism. We used 18 crossbred Iberian pigs with divergent postnatal growth patterns for whole genome sequencing and SNP discovery, with over 13 million variants being detected. We selected 1023 missense SNPs located on annotated genes and showing different allele frequencies between pigs with makerdly different growth patterns. We complemented this panel with 192 candidate SNPs obtained from literature mining and from muscle RNAseq data. The selected markers were genotyped in 480 Iberian × Duroc pigs from a commercial population, in which phenotypes were obtained, and an association study was performed for the 1005 successfully genotyped SNPs showing segregation. The results confirmed the effects of several known SNPs in candidate genes (such as LEPR, ACACA, FTO, LIPE or SCD on fatness, growth and fatty acid composition) and also disclosed interesting effects of new SNPs in less known genes such as LRIG3, DENND1B, SOWAHB, EPHX1 or NFE2L2 affecting body weight, average daily gain and adiposity at different ages, or KRT10, NLE1, KCNH2 or AHNAK affecting fatness and FA composition. The results provide a valuable basis for future implementation of marker-assisted selection strategies in swine and contribute to a better understanding of the genetic architecture of relevant traits.

DNA methylation landscapes from pig's limbic structures underline regulatory mechanisms relevant for brain plasticity

Epigenetic dynamics are essential for reconciling stress-induced responses in neuro-endocrine routes between the limbic brain and adrenal gland. CpG methylation associates with the initiation and end of regulatory mechanisms underlying responses critical for survival, and learning. Using Reduced Representation Bisulfite Sequencing, we identified methylation changes of functional relevance for mediating tissue-specific responses in the hippocampus, amygdala, hypothalamus, and adrenal gland in pigs. We identified 4186 differentially methylated CpGs across all tissues, remarkably, enriched for promoters of transcription factors (TFs) of the homeo domain and zinc finger classes. We also detected 5190 differentially methylated regions (DMRs, 748 Mb), with about half unique to a single pairwise. Two structures, the hypothalamus and the hippocampus, displayed 860 unique brain-DMRs, with many linked to regulation of chromatin, nervous development, neurogenesis, and cell-to-cell communication. TF binding motifs for TFAP2A and TFAP2C are enriched amount DMRs on promoters of other TFs, suggesting their role as master regulators, especially for pathways essential in long-term brain plasticity, memory, and stress responses. Our results reveal sets of TF that, together with CpG methylation, may serve as regulatory switches to modulate limbic brain plasticity and brain-specific molecular genetics in pigs.

Prenatal transcript levels and metabolomics analyses reveal metabolic changes associated with intrauterine growth restriction and sex

The metabolic changes associated with intrauterine growth restriction (IUGR) particularly affect the liver, which is a central metabolic organ and contributes significantly to the provision of energy and specific nutrients and metabolites. Therefore, our aim was to decipher and elucidate the molecular pathways of developmental processes mediated by miRNAs and mRNAs, as well as the metabolome in fetal liver tissue in IUGR compared to appropriate for gestational age groups (AGA). Discordant siblings representing the extremes in fetal weight at day 63 post conception (dpc) were selected from F2 fetuses of a cross of German Landrace and Pietrain. Most of the changes in the liver of IUGR at midgestation involved various lipid metabolic pathways, both on transcript and metabolite levels, especially in the category of sphingolipids and phospholipids. Differentially expressed miRNAs, such as miR-34a, and their differentially expressed mRNA targets were identified. Sex-specific phenomena were observed at both the transcript and metabolite levels, particularly in male. This suggests that sex-specific adaptations in the metabolic system occur in the liver during midgestation (63 dpc). Our multi-omics network analysis reveals interactions and changes in the metabolic system associated with IUGR and identified an important biosignature that differs between IUGR and AGA piglets.

Tissue-Wide Expression of Genes Related to Vitamin D Metabolism and FGF23 Signaling following Variable Phosphorus Intake in Pigs

Calcium (Ca) and phosphorus (P) homeostasis is maintained by several regulators, including vitamin D and fibroblast growth factor 23 (FGF23), and their tissue-specific activation and signaling cascades. In this study, the tissue-wide expression of key genes linked to vitamin D metabolism (CYP2R1, CYP27A1, CYP27B1, CYP24A1, GC, VDR) and FGF23 signaling (FGF23, FGFR1-4, KL) were investigated in pigs fed conventional (trial 1) and divergent P diets (trial 2). The tissue set comprised kidney, liver, bone, lung, aorta, and gastrointestinal tract sections. Expression patterns revealed that non-renal tissues and cells (NRTC) express genes to form active vitamin D [1,25(OH)₂D₃] according to site-specific requirements. A low P diet resulted in higher serum calcitriol and increased CYP24A1 expression in the small intestine, indicating local suppression of vitamin D signaling. A high P diet prompted increased mRNA abundances of CYP27B1 for local vitamin D synthesis, specifically in bone. For FGF23 signaling, analyses revealed ubiquitous expression of FGFR1-4, whereas KL was expressed in a tissue-specific manner. Dietary P supply did not affect skeletal FGF23; however, FGFR4 and KL showed increased expression in bone at high P supply, suggesting regulation to balance mineralization. Specific NRTC responses influence vitamin D metabolism and P homeostasis, which should be considered for a thrifty but healthy P supply.

Brain Transcriptome Responses to Dexamethasone Depending on Dose and Sex Reveal Factors Contributing to Sex-Specific Vulnerability to Stress-Induced Disorders

BACKGROUND

Glucocorticoid (GC) receptor (GR) signaling in the hypothalamus (Hyp) and in the superordinate limbic structures, such as the hippocampus (Hip), conveys feedback regulation of the neuroendocrine stress response and acts upon other neurobiological functions that ultimately influence mental health. These responses are strongly influenced by sex, but the molecular causes are still largely unexplored.

METHODS

To investigate GR targets and their GC sensitivity in the Hyp and Hip, we treated juvenile male and female piglets with 10 (D10) or 60 (D60) µg/kg dexamethasone (DEX), a selective GR agonist, and analyzed transcriptome responses compared to a saline control group using RNA sequencing.

RESULTS

Both doses influenced similar biological functions, including cellular response to lipid and immune cell-related functions, but the transcriptional response to D10 was considerably weaker, particularly in the Hip. Weighted Gene Co-expression Network Analysis revealed a network of genes coordinately regulated by DEX in both structures, among which the alpha-arrestin ARRDC2 takes a central position. Distinct functional groups of genes were differentially regulated by DEX between sexes depending on the dose; at D10, these included particularly mitochondrial genes, whereas at D60 interferon signaling and lipid homeostasis genes were enriched. The general and sex-specific transcriptional responses to DEX highlight microglia as the prominent target. Several key marker genes of disease-associated microglia were regulated by DEX depending on sex, such as TREM2 and LPL.

CONCLUSION

The discovered expression signatures suggest that DEX induced a dysfunctional state of microglia in males, while in females microglia were primed, which could entail predisposition for different mental disorders.

Multi-Transcript Level Profiling Revealed Distinct mRNA, miRNA, and tRNA-Derived Fragment Bio-Signatures for Coping Behavior Linked Haplotypes in HPA Axis and Limbic System

The molecular basis of porcine coping behavior (CB) relies on a sophisticated interplay of genetic and epigenetic features. Deep sequencing technologies allowed the identification of a plethora of new regulatory small non-coding RNA (sncRNA). We characterized mRNA and sncRNA profiles of central parts of the physiological stress response system including amygdala, hippocampus, hypothalamus and adrenal gland using systems biology for integration. Therefore, ten each of high- (HR) and low- (LR) reactive pigs (n = 20) carrying a CB associated haplotype in a prominent QTL-region on SSC12 were selected for mRNA and sncRNA expression profiling. The molecular markers related to the LR group included ATP1B2, MPDU1, miR-19b-5p, let-7g-5p, and 5′-tiRNA^Leu in the adrenal gland, miR-194a-5p, miR-125a-5p, miR-7-1-5p, and miR-107-5p in the hippocampus and CBL and PVRL1 in the hypothalamus. Interestingly, amygdalae of the LR group showed 5′-tiRNA and 5′-tRF (5′-tRF^Lys, 5′-tiRNA^Lys, 5′-tiRNA^Cys, and 5′-tiRNA^Gln) enrichment. Contrarily, molecular markers associated with the HR group encompassed miR-26b-5p, tRNA^Arg, tRNA^GlyiF in the adrenal gland, IGF1 and APOD in the amygdala and PBX1, TOB1, and C18orf1 in the hippocampus and miR-24 in the hypothalamus. In addition, hypothalami of the HR group were characterized by 3′-tiRNA enrichment (3′-tiRNAGln, 3′-tiRNA^Asn, 3′-tiRNA^Val, 3′-tRF^Pro, 3′-tiRNA^Cys, and 3′-tiRNA^Ala) and 3′-tRFs enrichment (3′-tRF^Asn, 3′-tRF^Glu, and 3′-tRF^Val). These evidence suggest that tRNA-derived fragments and their cleavage activity are a specific marker for coping behavior. Data integration revealed new bio-signatures of important molecular interactions on a multi-transcript level in HPA axis and limbic system of pigs carrying a CB-associated haplotype.

In Utero Fetal Weight in Pigs Is Regulated by microRNAs and Their Target Genes

Impaired skeletal muscle growth in utero can result in reduced birth weight and poor carcass quality in pigs. Recently, we showed the role of microRNAs (miRNAs) and their target genes in prenatal skeletal muscle development and pathogenesis of intrauterine growth restriction (IUGR). In this study, we performed an integrative miRNA-mRNA transcriptomic analysis in longissimus dorsi muscle (LDM) of pig fetuses at 63 days post conception (dpc) to identify miRNAs and genes correlated to fetal weight. We found 13 miRNAs in LDM significantly correlated to fetal weight, including miR-140, miR-186, miR-101, miR-15, miR-24, miR-29, miR-449, miR-27, miR-142, miR-99, miR-181, miR-199, and miR-210. The expression of these miRNAs decreased with an increase in fetal weight. We also identified 1315 genes significantly correlated to fetal weight at 63 dpc, of which 135 genes were negatively correlated as well as identified as potential targets of the above-listed 13 miRNAs. These miRNAs and their target genes enriched pathways and biological processes important for fetal growth, development, and metabolism. These results indicate that the transcriptomic profile of skeletal muscle can be used to predict fetal weight, and miRNAs correlated to fetal weight can serve as potential biomarkers of prenatal fetal health and growth.

A 192 bp ERV fragment insertion in the first intron of porcine TLR6 may act as an enhancer associated with the increased expressions of TLR6 and TLR1

Background

Toll-like receptors (TLRs) play important roles in building innate immune and inducing adaptive immune responses. Associations of the TLR genes polymorphisms with disease susceptibility, which are the basis of molecular breeding for disease resistant animals, have been reported extensively. Retrotransposon insertion polymorphisms (RIPs), as a new type of molecular markers developed recently, have great potential in population genetics and quantitative trait locus mapping. In this study, bioinformatic prediction combined with PCR-based amplification was employed to screen for RIPs in porcine TLR genes. Their population distribution was examined, and for one RIP the impact on gene activity and phenotype was further evaluated.

Results

Five RIPs, located at the 3' flank of TLR3, 5' flank of TLR5, intron 1 of TLR6, intron 1 of TLR7, and 3' flank of TLR8 respectively, were identified. These RIPs were detected in different breeds with an uneven distribution among them. By using the dual luciferase activity assay a 192 bp endogenous retrovirus (ERV) in the intron 1 of TLR6 was shown to act as an enhancer increasing the activities of TLR6 putative promoter and two mini-promoters. Furthermore, real-time quantitative polymerase chain reaction (qPCR) analysis revealed significant association (p < 0.05) of the ERV insertion with increased mRNA expression of TLR6, the neighboring gene TLR1, and genes downstream in the TLR signaling pathway such as MyD88 (Myeloid differentiation factor 88), Rac1 (Rac family small GTPase 1), TIRAP (TIR domain containing adaptor protein), Tollip (Toll interacting protein) as well as the inflammatory factors IL6 (Interleukin 6), IL8 (Interleukin 8), and TNFα (Tumor necrosis factor alpha) in tissues of 30 day-old piglet. In addition, serum IL6 and TNFα concentrations were also significantly upregulated by the ERV insertion (p < 0.05).

Conclusions

A total of five RIPs were identified in five different TLR loci. The 192 bp ERV insertion in the first intron of TLR6 was associated with higher expression of TLR6, TLR1, and several genes downstream in the signaling cascade. Thus, the ERV insertion may act as an enhancer affecting regulation of the TLR signaling pathways, and can be potentially applied in breeding of disease resistant animals.

Supplementary Information

The online version contains supplementary material available at 10.1186/s13100-021-00248-w.

Genome-Wide Analysis for Early Growth-Related Traits of the Locally Adapted Egyptian Barki Sheep

Sheep play a critical role in the agricultural and livestock sector in Egypt. For sheep meat production, growth traits such as birth and weaning weights are very important and determine the supply and income of local farmers. The Barki sheep originates from the northeastern coastal zone of Africa, and due to its good adaptation to the harsh environmental conditions, it contributes significantly to the meat production in these semi-arid regions. This study aimed to use a genome-wide SNP panel to identify genomic regions that are diversified between groups of individuals of Egyptian Barki sheep with high and low growth performance traits. In this context, from a phenotyped population of 140 lambs of Barki sheep, 69 lambs were considered for a genome-wide scan with the Illumina OvineSNP50 V2 BeadChip. The selected lambs were grouped into divergent subsets with significantly different performance for birth weight and weaning weight. After quality control, 63 animals and 40,383 SNPs were used for analysis. The fixation index (F_ST) for each SNP was calculated between the groups. The results verified genomic regions harboring some previously proposed candidate genes for traits related to body growth, i.e., EYA2, GDF2, GDF10, MEF2B, SLC16A7, TBX15, TFAP2B, and TNNC2. Moreover, novel candidate genes were proposed with known functional implications on growth processes such as CPXM2 and LRIG3. Subsequent association analysis showed significant effects of the considered SNPs on birth and weaning weights. Results highlight the genetic diversity associated with performance traits and thus the potential to improve growth traits in the Barki sheep breed.

Transcriptome analysis of porcine PBMCs reveals lipopolysaccharide-induced immunomodulatory responses and crosstalk of immune and glucocorticoid receptor signaling

The current level of knowledge on transcriptome responses triggered by endotoxins and glucocorticoids in immune cells in pigs is limited. Therefore, in the present study, we treated porcine peripheral blood mononuclear cells (PBMCs) with lipopolysaccharide (LPS) and dexamethasone (DEX) separately or combined for 2 hours. The resultant transcriptional responses were examined by mRNA sequencing. We found that the LPS treatment triggered pronounced inflammatory responses as evidenced by upregulation of pro-inflammatory cytokines, chemokines, and related signaling pathways like NF-κB. Concurrently, a series of downregulated pro-inflammatory and upregulated anti-inflammatory molecules were identified. These are involved in the inhibition of TLR, NF-κB, and MAPK cascades and activation of signaling mediated by Tregs and STAT3, respectively. These findings suggested that LPS initiated also an anti-inflammatory process to prevent an overwhelming inflammatory response. The transcriptome responses further revealed substantial crosstalk of immune responses and glucocorticoid receptor (GR) signaling. This was apparent in four aspects: constitutive inhibition of T cell signaling by DEX through a subset of genes showing no response to LPS; inhibition of LPS-induced inflammatory genes by DEX; attenuation of DEX action by LPS paralleled by the regulation of genes implicated in cytokine and calcium signaling; and DEX-induced changes in genes associated with the activation of pro-inflammatory TLR, NF-κB, iNOS, and IL-1 signaling. Consequently, our study provides novel insights into inflammatory and GR signaling in pigs, as well as an understanding of the application of glucocorticoid drugs for the treatment of inflammatory disorders.

SINE jumping contributes to large-scale polymorphisms in the pig genomes

BACKGROUND

Molecular markers based on retrotransposon insertion polymorphisms (RIPs) have been developed and are widely used in plants and animals. Short interspersed nuclear elements (SINEs) exert wide impacts on gene activity and even on phenotypes. However, SINE RIP profiles in livestock remain largely unknown, and not be revealed in pigs.

RESULTS

Our data revealed that SINEA1 displayed the most polymorphic insertions (22.5 % intragenic and 26.5 % intergenic), followed by SINEA2 (10.5 % intragenic and 9 % intergenic) and SINEA3 (12.5 % intragenic and 5.0 % intergenic). We developed a genome-wide SINE RIP mining protocol and obtained a large number of SINE RIPs (36,284), with over 80 % accuracy and an even distribution in chromosomes (14.5/Mb), and 74.34 % of SINE RIPs generated by SINEA1 element. Over 65 % of pig SINE RIPs overlap with genes, most of them (> 95 %) are in introns. Overall, about one forth (23.09 %) of the total genes contain SINE RIPs. Significant biases of SINE RIPs in the transcripts of protein coding genes were observed. Nearly half of the RIPs are common in these pig breeds. Sixteen SINE RIPs were applied for population genetic analysis in 23 pig breeds, the phylogeny tree and cluster analysis were generally consistent with the geographical distributions of native pig breeds in China.

CONCLUSIONS

Our analysis revealed that SINEA1-3 elements, particularly SINEA1, are high polymorphic across different pig breeds, and generate large-scale structural variations in the pig genomes. And over 35,000 SINE RIP markers were obtained. These data indicate that young SINE elements play important roles in creating new genetic variations and shaping the evolution of pig genome, and also provide strong evidences to support the great potential of SINE RIPs as genetic markers, which can be used for population genetic analysis and quantitative trait locus (QTL) mapping in pig.

Prenatal Skeletal Muscle Transcriptome Analysis Reveals Novel MicroRNA-mRNA Networks Associated with Intrauterine Growth Restriction in Pigs

Intrauterine growth restriction (IUGR) occurs in 15–20% of pig neonates and poses huge economic losses to the pig industry. IUGR piglets have reduced skeletal muscle growth, which may persist after birth. Prenatal muscle growth is regulated by complex molecular pathways that are not well understood. MicroRNAs (miRNAs) have emerged as the main regulators of vital pathways and biological processes in the body. This study was designed to identify miRNA–mRNA networks regulating prenatal skeletal muscle development in pigs. We performed an integrative miRNA–mRNA transcriptomic analysis in longissimus dorsi muscle from IUGR fetuses and appropriate for gestational age (AGA) fetuses at 63 days post conception. Our data showed that 47 miRNAs and 3257 mRNAs were significantly upregulated, and six miRNAs and 477 mRNAs were significantly downregulated in IUGR compared to AGA fetuses. Moreover, 47 upregulated miRNAs were negatively correlated and can potentially target 326 downregulated genes, whereas six downregulated miRNAs were negatively correlated and can potentially target 1291 upregulated genes. These miRNA–mRNA networks showed enrichment in biological processes and pathways critical for fetal growth, development, and metabolism. The miRNA–mRNA networks identified in this study can potentially serve as indicators of prenatal fetal growth and development as well as postnatal carcass quality.

A natural Ala610Val substitution causing glucocorticoid receptor hypersensitivity aggravates consequences of endotoxemia

Despite the crucial role of glucocorticoid receptor (GR) in proper immune responses, the effect of GR hypersensitivity on inflammation is rarely reported. To fill this knowledge gap, we exploited the natural gain-of-function substitution in the porcine glucocorticoid receptor (GR_Ala610Val) and challenged pigs carrying normal or hypersensitive GR using 50 µg/kg lipopolysaccharide (LPS) following pretreatment with either saline or single bolus of 60 µg/kg dexamethasone (DEX). The GR_Ala610Val substitution reduced baseline cortisol, adrenocorticotropic hormone (ACTH), and triglyceride concentration and granulocyte proportion whereas baseline platelet counts were elevated. Val-carriers, i.e. AlaVal as well as ValVal pigs, showed less LPS-induced cortisol rise but the cortisol fold change was similar in all genotypes. Differently, ACTH response to LPS was most significant in GR_Ala610Val heterozygotes (AlaVal). LPS-induced disorders, including sickness behaviors, anorexia, thrombocytopenia, cytokine production, and metabolic alterations were more intense in Val-carriers. On the other hand, Val-carriers were more sensitive to DEX effect than wild types (AlaAla) during endotoxemia, but not under unchallenged conditions. This is the first report revealing aggravated responses to endotoxemia by GR gain-of-function. Together, these results imply that GR hypersensitivity is difficult to diagnose but may represent a risk factor for endotoxemia and sepsis.

Deep sequencing of small non-coding RNA highlights brain-specific expression patterns and RNA cleavage

With the advance of high-throughput sequencing technology numerous new regulatory small RNAs have been identified, that broaden the variety of processing mechanisms and functions of non-coding RNA. Here we explore small non-coding RNA (sncRNA) expression in central parts of the physiological stress and anxiety response system. Therefore, we characterize the sncRNA profile of tissue samples from Amygdala, Hippocampus, Hypothalamus and Adrenal Gland, obtained from 20 pigs. Our analysis reveals that all tissues but Amygdala and Hippocampus possess distinct, tissue-specific expression pattern of miRNA that are associated with Hypoxia, stress responses as well as memory and fear conditioning. In particular, we observe marked differences in the expression profile of limbic tissues compared to those associated to the HPA/stress axis, with a surprisingly high aggregation of 3´-tRNA halves in Amygdala and Hippocampus. Since regulation of sncRNA and RNA cleavage plays a pivotal role in the central nervous system, our work provides seminal insights in the role/involvement of sncRNA in the transcriptional and post-transcriptional regulation of negative emotion, stress and coping behaviour in pigs, and mammals in general.

Genetic Contribution to Variation in Blood Calcium, Phosphorus, and Alkaline Phosphatase Activity in Pigs

Blood values of calcium (Ca), inorganic phosphorus (IP), and alkaline phosphatase activity (ALP) are valuable indicators for mineral status and bone mineralization. The mineral homeostasis is maintained by absorption, retention, and excretion processes employing a number of known and unknown sensing and regulating factors with implications on immunity. Due to the high inter-individual variation of Ca and P levels in the blood of pigs and to clarify molecular contributions to this variation, the genetics of hematological traits related to the Ca and P balance were investigated in a German Landrace population, integrating both single-locus and multi-locus genome-wide association study (GWAS) approaches. Genomic heritability estimates suggest a moderate genetic contribution to the variation of hematological Ca (N = 456), IP (N = 1049), ALP (N = 439), and the Ca/P ratio (N = 455), with values ranging from 0.27 to 0.54. The genome-wide analysis of markers adds a number of genomic regions to the list of quantitative trait loci, some of which overlap with previous results. Despite the gaps in knowledge of genes involved in Ca and P metabolism, genes like THBS2, SHH, PTPRT, PTGS1, and FRAS1 with reported connections to bone metabolism were derived from the significantly associated genomic regions. Additionally, genomic regions included TRAFD1 and genes coding for phosphate transporters (SLC17A1-SLC17A4), which are linked to Ca and P homeostasis. The study calls for improved functional annotation of the proposed candidate genes to derive features involved in maintaining Ca and P balance. This gene information can be exploited to diagnose and predict characteristics of micronutrient utilization, bone development, and a well-functioning musculoskeletal system in pig husbandry and breeding.

Haplotypes of coping behavior associated QTL regions reveal distinct transcript profiles in amygdala and hippocampus

Stress response and coping behavior in pigs are largely shaped by hypothalamic-pituitary-adrenal axis and sympatho-adrenomedullary system action. However, the dynamic interaction between amygdala and hippocampus crucially modulates the behavioral response towards significant emotional events. While this functional relationship is well documented, the molecular underpinnings still remain insufficiently understood. Our study used transcriptome profiling of porcine amygdala and hippocampus to identify molecular pathways that are differentially activated depending on the haplotype of a significantly coping behavior-associated region on pig chromosome 12 (SSC12). The pigs were classified into two groups based on the haplotype information of this QTL-region discovered in our previous genome-wide association study. Ten each of high- (HR) and low- (LR) reactive pigs (n = 20) were selected for differential gene expression analysis and weighted gene co-expression analysis with subsequent pathway analysis. Differentially expressed genes identified in the amygdala include SELL, CXCR7 and NTS, while TRAF3, PTGS2 and CFI were detected in the hippocampus indicating a role of neuroinflammation and immunological processes. Pathway analysis revealed IL-8 signaling, NF-κB signaling, glutamate and GABA metabolism, glucocorticoid receptor signaling and chemokine signaling in the amygdala and ephrin receptor signaling, as well as NF-κB signaling in the hippocampus. We discovered candidate genes in regions detected by genome-wide association study including ARRB2, ADRBK2, THRB, NEK7 and ACVR2B, which relate to dopaminergic and other monoaminergic neurotransmitter systems, neuroimmunomodulation, neuroinflammation and GABA-ergic neurotransmission. These findings provide insights into the molecular underpinning of divergent coping behavior and associated haplotypes in limbic forebrain system in pig.

Epigenome-wide skeletal muscle DNA methylation profiles at the background of distinct metabolic types and ryanodine receptor variation in pigs

BACKGROUND

Epigenetic variation may result from selection for complex traits related to metabolic processes or appear in the course of adaptation to mediate responses to exogenous stressors. Moreover epigenetic marks, in particular the DNA methylation state, of specific loci are driven by genetic variation. In this sense, polymorphism with major gene effects on metabolic and cell signaling processes, like the variation of the ryanodine receptors in skeletal muscle, may affect DNA methylation.

METHODS

DNA-Methylation profiles were generated applying Reduced Representation Bisulfite Sequencing (RRBS) on 17 Musculus longissimus dorsi samples. We examined DNA methylation in skeletal muscle of pig breeds differing in metabolic type, Duroc and Pietrain. We also included F2 crosses of these breeds to get a first clue to DNA methylation sites that may contribute to breed differences. Moreover, we compared DNA methylation in muscle tissue of Pietrain pigs differing in genotypes at the gene encoding the Ca2+ release channel (RYR1) that largely affects muscle physiology.

RESULTS

More than 2000 differently methylated sites were found between breeds including changes in methylation profiles of METRNL, IDH3B, COMMD6, and SLC22A18, genes involved in lipid metabolism. Depending on RYR1 genotype there were 1060 differently methylated sites including some functionally related genes, such as CABP2 and EHD, which play a role in buffering free cytosolic Ca²⁺ or interact with the Na⁺/Ca²⁺ exchanger.

CONCLUSIONS

The change in the level of methylation between the breeds is probably the result of the long-term selection process for quantitative traits involving an infinite number of genes, or it may be the result of a major gene mutation that plays an important role in muscle metabolism and triggers extensive compensatory processes.

Transcriptome Responses to Dexamethasone Depending on Dose and Glucocorticoid Receptor Sensitivity in the Liver

Tissue sensitivity to glucocorticoids is a key factor dictating outcome of their homeostatic and therapeutic action, whereby liver represents one of the major peripheral targets. Here, we used pigs carrying a natural gain-of-function glucocorticoid receptor (GR) variant Ala610Val (GR_Ala610Val) as a model to identify genes and pathways related to differential glucocorticoid sensitivity. Animals with different GR_Ala610Val genotypes were treated either with saline or two different doses of dexamethasone. Genome-wide transcriptional responses depending on treatment, genotype, and their interaction in the liver were investigated using mRNA sequencing. Dexamethasone induced vast transcriptional responses, with more than 30% of present genes being affected. Functional annotation of genes differentially expressed due to dexamethasone treatment suggested that genes related to inflammation respond more sensitively, despite absence of an immune stimulus. In contrast, genes involved in glucose metabolism and cancer appeared to be less sensitive. Analysis of genotype and genotype × treatment interaction revealed that clustered protocadherins, particularly PCDHB7, are most prominently affected by GR_Ala610Val, mainly depending on dose. GR_Ala610Val influenced also expression of a set of glucose metabolism related genes, including PPARGC1A and CEBPB, in the absence of dexamethasone though no differences in basal plasma glucose level were observed. This might represent an adaptive response, keeping balance between receptor sensitivity, and level of circulating endogenous glucocorticoids. Administration of low dexamethasone dose changed their expression pattern and induced higher glucose response in carriers of the hypersensitive Val receptor. Our findings suggest that GR_Ala610Val modulates tissue responses to glucocorticoids dynamically, depending on their circulating level.

Genetic Regulation of Liver Metabolites and Transcripts Linking to Biochemical-Clinical Parameters

Given the central metabolic role of the liver, hepatic metabolites and transcripts reflect the organismal physiological state. Biochemical-clinical plasma biomarkers, hepatic metabolites, transcripts, and single nucleotide polymorphism (SNP) genotypes of some 300 pigs were integrated by weighted correlation networks and genome-wide association analyses. Network-based approaches of transcriptomic and metabolomics data revealed linked of transcripts and metabolites of the pentose phosphate pathway (PPP). This finding was evidenced by using a NADP/NADPH assay and HDAC4 and G6PD transcript quantification with the latter coding for first limiting enzyme of this pathway and by RNAi knockdown experiments of HDAC4. Other transcripts including ARG2 and SLC22A7 showed link to amino acids and biomarkers. The amino acid metabolites were linked with transcripts of immune or acute phase response signaling, whereas the carbohydrate metabolites were highly enrich in cholesterol biosynthesis transcripts. Genome-wide association analyses revealed 180 metabolic quantitative trait loci (mQTL) (p < 10^-4). Trans-4-hydroxy-L-proline (p = 6 × 10^-9), being strongly correlated with plasma creatinine (CREA), showed strongest association with SNPs on chromosome 6 that had pleiotropic effects on PRODH2 expression as revealed by multivariate analysis. Consideration of shared marker association with biomarkers, metabolites, and transcripts revealed 144 SNPs associated with 44 metabolites and 69 transcripts that are correlated with each other, representing 176 mQTL and expression quantitative trait loci (eQTL). This is the first work to report genetic variants associated with liver metabolite and transcript levels as well as blood biochemical-clinical parameters in a healthy porcine model. The identified associations provide links between variation at the genome, transcriptome, and metabolome level molecules with clinically relevant phenotypes. This approach has the potential to detect novel biomarkers displaying individual variation and promoting predictive biology in medicine and animal breeding.

Kinetics of Physiological and Behavioural Responses in Endotoxemic Pigs with or without Dexamethasone Treatment

Although dexamethasone (DEX) is a widely used immunoregulatory agent, knowledge about its pharmacological properties in farm animals, especially pigs, is insufficient. Previous studies suggest that compared to other species, pigs are less sensitive to the immunosuppression conferred by DEX and more sensitive to the threat of bacterial endotoxins. However, there is a paucity of studies examining DEX immunomodulation in endotoxemia in this species. In this study, a porcine endotoxemia model was established by lipopolysaccharide (LPS) and the effect of DEX-pretreatment on the magnitude and kinetics of neuroendocrine, metabolic, hematologic, inflammatory, and behavioural responses were examined. DEX decreased cortisol, adrenocorticotropic hormone (ACTH), red blood cell, hemoglobin, hematocrit, and lymphocyte whereas glucose concentration was increased under both normal and endotoxemic conditions. By contrast, DEX decreased triglyceride, lactate, and IL-6 concentrations and increased platelet count only under an endotoxemic condition. DEX also reduced the frequency of sickness behaviour following LPS challenge. PCA showed that glucose and triglyceride metabolism together with red blood cell count mainly contributed to the separation of clusters during DEX treatment. Our study demonstrates that DEX protects pigs from inflammation and morbidity in endotoxemia, in spite of their less sensitivity to DEX. Moreover, its considerable role in the regulation of the metabolic and hematologic responses in endotoxemic pigs is revealed for the first time.

Genetic architecture and regulatory impact on hepatic microRNA expression linked to immune and metabolic traits

Regulation of microRNA (miRNA) expression contributes to a wide range of target gene expression and phenotypes. The miRNA expression in the liver, the central metabolic organ, was examined in 209 pigs, and integrated with haematological and clinical biomarkers of metabolic and overall health, mRNA-target expression levels and single-nucleotide polymorphism (SNP) genotypes. The expression levels of 426 miRNA species correlated with plasma haematological or biochemical traits (r² = |0.19–0.45|, false discovery rate < 5%). Pairs of these miRNAs and their predicted target mRNAs showing expressing levels associated with the identical traits were examined to understand how immune and metabolic traits are affected by miRNA-mediated regulatory networks derived by mapping miRNA abundance as an expression quantitative trait. In total, 221 miRNA-expression-QTL correspond to 164 SNPs and 108 miRNAs, including miR-34a, miR-30e, miR-148-3p, miR-204, miR-181-5p, miR-143-5p and let-7 g that also correlate with the biomarkers. Sixty-one SNPs were simultaneously associated with 29 miRNA and 41 mRNA species. The expression levels of 13 out of 29 miRNA were correlated with one of the biochemical or haematological traits. For example, the expression levels of miR-34a were correlated with serum phosphorus and cholesterin levels; miR-204, miR-15a and miR-16b were correlated with triglyceride. For haematological traits, the expression levels of miR-652 and miR-204 were correlated with the mean corpuscular haemoglobin concentration, and the expression of miR-143 was correlated with plateletcrit. Pleiotropic association analyses revealed genetic links between mRNA and miRNA on SSC6 for miR-34a, SSC9 for miR-708 and SSC14 for miR-652. Our analysis of miRNA and mRNA transcript profiles, their correlation with clinically important plasma parameters of hepatic functions as well as information on their genetic regulation provide novel regulatory networks and potential new biomarkers for immune and metabolic traits.

Lowered dietary phosphorus affects intestinal and renal gene expression to maintain mineral homeostasis with immunomodulatory implications in weaned piglets

BACKGROUND

In monogastric animals, phosphorus (P) homeostasis is maintained by regulating intestinal absorption, bone mobilization, and renal excretion. Since P is a non-renewable resource, a shortage is imminent due to widespread over-usage in the farming and animal husbandry industries. As a consequence, P efficiency should be improved in pig production. We sought to characterize the transcriptional response in re-/absorbing and excreting tissues in pigs to diets varying in calcium: phosphorus ratios. Weaned piglets were assigned to one of three groups fed diets varying in digestible P content for a period of five weeks. Gene expression profiles were analyzed in jejunum, colon, and kidney.

RESULTS

Transcriptome analysis revealed that reduced dietary P intake affects gene expression in jejunum and kidney, but not in colon. The regulation of mineral homeostasis was reflected via altered mRNA abundances of CYP24A1, CYP27A1, TRPM6, SPP1, and VDR in jejunum and kidney. Moreover, lowered abundances of transcripts associated with the classical complement system pathway were observed in the jejunum. In kidney, shifted transcripts were involved in phospholipase C, calcium signaling, and NFAT signaling, which may have immunomodulatory implications.

CONCLUSIONS

Our results revealed local transcriptional consequences of variable P intake in intestinal and renal tissues. The adaptive responses are the result of organismal efforts to maintain systemic mineral homeostasis while modulating immune features at local tissue sites. Therefore, the deviation from the currently recommended dietary P supply must be carefully considered, as the endogenous mechanisms that respond to low P diets may impact important adaptive immune responses.

Mitochondrial-nuclear crosstalk, haplotype and copy number variation distinct in muscle fiber type, mitochondrial respiratory and metabolic enzyme activities

Genes expressed in mitochondria work in concert with those expressed in the nucleus to mediate oxidative phosphorylation (OXPHOS), a process that is relevant for muscle metabolism and meat quality. Mitochondrial genome activity can be efficiently studied and compared in Duroc and Pietrain pigs, which harbor different mitochondrial haplotypes and distinct muscle fiber types, mitochondrial respiratory activities, and fat content. Pietrain pigs homozygous-positive for malignant hyperthermia susceptibility (PiPP) carried only haplotype 8 and showed the lowest absolute mtDNA copy number accompanied by a decrease transcript abundance of mitochondrial-encoded subunits ND1, ND6, and ATP6 and nuclear-encoded subunits NDUFA11 and NDUFB8. In contrast, we found that haplotype 4 of Duroc pigs had significantly higher mitochondrial DNA (mtDNA) copy numbers and an increase transcript abundance of mitochondrial-encoded subunits ND1, ND6, and ATP6. These results suggest that the variation in mitochondrial and nuclear genetic background among these animals has an effect on mitochondrial content and OXPHOS system subunit expression. We observed the co-expression pattern of mitochondrial and nuclear encoded OXPHOS subunits suggesting that the mitochondrial-nuclear crosstalk functionally involves in muscle metabolism. The findings provide valuable information for understanding muscle biology processes and energy metabolism, and may direct use for breeding strategies to improve meat quality and animal health.

Exploring the genetics of feed efficiency and feeding behaviour traits in a pig line highly selected for performance characteristics

The consideration of feed efficiency traits is highly relevant in animal breeding due to economic and ecologic impacts of the efficient usage and utilization of feed resources. In pigs, corresponding observations are recorded using automatic feeding stations and serve as one of the main criteria in most pig selection programmes. Simultaneously, feeding stations also generate feeding behaviour data which represent a nearly unused resource and provide a valuable proxy measure of health status, animal welfare, and management practices. In the current study, an integrated approach was applied to a feed efficiency tested and genome-wide genotyped terminal sire line population. Therefore, genetic analyses were performed combining a single-marker based approach and a Bayesian multi-marker algorithm. Major quantitative trait loci (QTL) for feeding behaviour traits comprising daily occupation time, daily feeder visit, and daily feeding rate were identified on chromosomes 1, 4, 6, 7, 8, and 14. Feed efficiency was represented by feed conversion ratio and daily feed intake revealing prominent genomic regions on chromosomes 1, 6, 9, and 11. The positional and functional candidate genes identified are involved in transport processes like AQP4, SLC22A23, and SLC6A14 as well as energy sensing, generation, and utilization as exemplified by PPP3CA, IQGAP3, ECI2, and DnaJC15. These molecular features provide the first step towards the dissection of the genetic connection between distinct feeding behaviour patterns, feed efficiency and performance, health, and welfare traits driving the implementation of these traits in breeding programmes and pig husbandry.

Genetically regulated hepatic transcripts and pathways orchestrate haematological, biochemical and body composition traits

The liver is the central metabolic organ and exhibits fundamental functions in haematological traits. Hepatic expression, haematological, plasma biochemical, and body composition traits were assessed in a porcine model (n = 297) to establish tissue-specific genetic variations that influence the function of immune-metabolism-correlated expression networks. At FDR (false discovery rate) <1%, more than 3,600 transcripts were jointly correlated (r = |0.22-0.48|) with the traits. Functional enrichment analysis demonstrated common links of metabolic and immune traits. To understand how immune and metabolic traits are affected via genetic regulation of gene expression, eQTLs were assessed. 20517 significant (FDR < 5%) eQTLs for 1401 transcripts were identified, among which 443 transcripts were associated with at least one of the examined traits and had cis-eQTL (such as ACO1 (6.52 × 10^-7) and SOD1 (6.41 × 10^-30). The present study establishes a comprehensive view of hepatic gene activity which links together metabolic and immune traits in a porcine model for medical research.

Single- and Bayesian Multi-Marker Genome-Wide Association for Haematological Parameters in Pigs

Haematological traits are important traits that show associations with immune and metabolic status, as well as diseases in humans and animals. Mapping genome regions that affect the blood cell traits can contribute to the identification of genomic features useable as biomarkers for immune, disease and metabolic status. A genome-wide association study (GWAS) was conducted using PorcineSNP60 BeadChips. Single-marker and Bayesian multi-marker approaches were integrated to identify genomic regions and corresponding genes overlapping for both methods. GWAS was performed for haematological traits of 591 German Landrace pig. Heritability estimates for haematological traits were medium to high. In total 252 single SNPs associated with 12 haematological traits were identified (NegLog10 of p-value > 5). The Bayesian multi-marker approach revealed 102 QTL regions across the genome, indicated by 1-Mb windows with contribution to additive genetic variance above 0.5%. The integration of both methods resulted in 24 overlapping QTL regions. This study identified overlapping QTL regions from single- and multi-marker approaches for haematological traits. Identifying candidate genes that affect blood cell traits provides the first step towards the understanding of the molecular basis of haematological phenotypes.

Integrated Genome-wide association and hypothalamus eQTL studies indicate a link between the circadian rhythm-related gene PER1 and coping behavior

Animal personality and coping styles are basic concepts for evaluating animal welfare. Struggling response of piglets in so-called backtests early in life reflects their coping strategy. Behavioral reactions of piglets in backtests have a moderate heritability, but their genetic basis largely remains unknown. Here, latency, duration and frequency of struggling attempts during one-minute backtests were repeatedly recorded of piglets at days 5, 12, 19, and 26. A genome-wide association study for backtest traits revealed 465 significant SNPs (FDR ≤ 0.05) mostly located in QTL (quantitative trait locus) regions on chromosome 3, 5, 12 and 16. In order to capture genes in these regions, 37 transcripts with significant SNPs were selected for expressionQTL analysis in the hypothalamus. Eight genes (ASGR1, CPAMD8, CTC1, FBXO39, IL19, LOC100511790, RAD51B, UBOX5) had cis- and five (RANGRF, PER1, PDZRN3, SH2D4B, LONP2) had trans-expressionQTL. In particular, for PER1, with known physiological implications for maintenance of circadian rhythms, a role in coping behavior was evidenced by confirmed association in an independent population. For CTC1 a cis-expression QTL and the consistent relationship of gene polymorphism, mRNA expression level and backtest traits promoted its link to coping style. GWAS and eQTL analyses uncovered positional and functional gene candidates for coping behavior.

PBMC transcriptomic responses to primary and secondary vaccination differ due to divergent lean growth and antibody titers in a pig model

The genetic relationship between immune responsiveness and performance is not well understood, but a major topic of the evolution of resource allocation and of relevance in human medicine and livestock breeding, for instance. This study aims to show differences of transcript abundance changes during the time intervals before and after two tetanus toxoid (TT) vaccinations in domestic pigs differing in lean growth (LG) and anti-TT-antibody titers (AB) parameters of performance and immunocompetence. During response to the first vaccination all animals had a general decrease in transcript abundances related to various functional pathways as measured by comparative Affymetrix microarray hybridization and Ingenuity Pathway analyses. Low-AB phenotypes had predominantly decreased immune response transcripts. Combined phenotypes high-AB/high-LG had decreased transcripts related to growth factor signaling pathways. However, during the interval after the booster vaccination, high-LG and high-AB animals responded exclusively with increased immune transcripts, such as B-cell receptor signaling and cellular immune response pathways. In addition, high-LG and low-AB animals had predominantly increased transcripts of several cellular immune response pathways. Conversely, low-LG and high-AB animals had few elevated immune transcripts and decreased transcripts related to only two nonimmune-specific pathways. In response to booster vaccination high-LG phenotypes revealed enriched transcripts related to several different immune response pathways, regardless of AB phenotype. Different from the expected impact of AB titers, divergent AB phenotypes did not reflect considerable differences between immune transcripts. However, highly significant differences observed among divergent LG phenotypes suggest the compatibility of high performance and high vaccine responses.

MicroRNAs Regulate Cellular ATP Levels by Targeting Mitochondrial Energy Metabolism Genes during C2C12 Myoblast Differentiation

In our previous study, we identified an miRNA regulatory network involved in energy metabolism in porcine muscle. To better understand the involvement of miRNAs in cellular ATP production and energy metabolism, here we used C2C12 myoblasts, in which ATP levels increase during differentiation, to identify miRNAs modulating these processes. ATP level, miRNA and mRNA microarray expression profiles during C2C12 differentiation into myotubes were assessed. The results suggest 14 miRNAs (miR-423-3p, miR-17, miR-130b, miR-301a/b, miR-345, miR-15a, miR-16a, miR-128, miR-615, miR-1968, miR-1a/b, and miR-194) as cellular ATP regulators targeting genes involved in mitochondrial energy metabolism (Cox4i2, Cox6a2, Ndufb7, Ndufs4, Ndufs5, and Ndufv1) during C2C12 differentiation. Among these, miR-423-3p showed a high inverse correlation with increasing ATP levels. Besides having implications in promoting cell growth and cell cycle progression, its function in cellular ATP regulation is yet unknown. Therefore, miR-423-3p was selected and validated for the function together with its potential target, Cox6a2. Overexpression of miR-423-3p in C2C12 myogenic differentiation lead to decreased cellular ATP level and decreased expression of Cox6a2 compared to the negative control. These results suggest miR-423-3p as a novel regulator of ATP/energy metabolism by targeting Cox6a2.

Identification of common regulators of genes in co-expression networks affecting muscle and meat properties

Understanding the genetic contributions behind skeletal muscle composition and metabolism is of great interest in medicine and agriculture. Attempts to dissect these complex traits combine genome-wide genotyping, expression data analyses and network analyses. Weighted gene co-expression network analysis (WGCNA) groups genes into modules based on patterns of co-expression, which can be linked to phenotypes by correlation analysis of trait values and the module eigengenes, i.e. the first principal component of a given module. Network hub genes and regulators of the genes in the modules are likely to play an important role in the emergence of respective traits. In order to detect common regulators of genes in modules showing association with meat quality traits, we identified eQTL for each of these genes, including the highly connected hub genes. Additionally, the module eigengene values were used for association analyses in order to derive a joint eQTL for the respective module. Thereby major sites of orchestrated regulation of genes within trait-associated modules were detected as hotspots of eQTL of many genes of a module and of its eigengene. These sites harbor likely common regulators of genes in the modules. We exemplarily showed the consistent impact of candidate common regulators on the expression of members of respective modules by RNAi knockdown experiments. In fact, Cxcr7 was identified and validated as a regulator of genes in a module, which is involved in the function of defense response in muscle cells. Zfp36l2 was confirmed as a regulator of genes of a module related to cell death or apoptosis pathways. The integration of eQTL in module networks enabled to interpret the differentially-regulated genes from a systems perspective. By integrating genome-wide genomic and transcriptomic data, employing co-expression and eQTL analyses, the study revealed likely regulators that are involved in the fine-tuning and synchronization of genes with trait-associated expression.

Pre- and post-natal muscle microRNA expression profiles of two pig breeds differing in muscularity

miRNAs regulate the expression of target genes in diverse cellular processes and hence play important roles in physiological processes including developmental timing, patterning, embryogenesis, organogenesis, cell lineage, myogenesis and growth control. A comparative expression analysis of miRNAs expressed in the longissimus dorsi muscle at two prenatal stages (63 and 91 days post-conception (dpc)), and one adult stage (180 days post-natum) in both German Landrace (DL) and Pietrain (Pi) pig breeds was performed using a custom-designed array. During the prenatal stages, miR-199 and the miR-17 families were significantly up-regulated at 63 dpc, whereas miR-1 and miR-133a were overexpressed at 91 dpc. The abundance of several miRNAs was increased in the adult stage compared to 91 dpc including miR-1, miR-133, miR-22(a/b) and miR-29a. Some miRNAs were breed-specific, such as miR-199 and the miR-17 families which were all up-regulated in Pi pigs, while miR-133, miR-181 and miR-214 were up-regulated in DL pigs. Several pathways related to muscle development were enriched with predicted targets for the differentially expressed miRNAs. The dynamic expression and breed-associated regulation of porcine muscle miRNAs suggests a functional role for miRNA-mediated gene regulation during muscle development and phenotypic variations of muscle traits.

biomvRhsmm: genomic segmentation with hidden semi-Markov model

High-throughput technologies like tiling array and next-generation sequencing (NGS) generate continuous homogeneous segments or signal peaks in the genome that represent transcripts and transcript variants (transcript mapping and quantification), regions of deletion and amplification (copy number variation), or regions characterized by particular common features like chromatin state or DNA methylation ratio (epigenetic modifications). However, the volume and output of data produced by these technologies present challenges in analysis. Here, a hidden semi-Markov model (HSMM) is implemented and tailored to handle multiple genomic profile, to better facilitate genome annotation by assisting in the detection of transcripts, regulatory regions, and copy number variation by holistic microarray or NGS. With support for various data distributions, instead of limiting itself to one specific application, the proposed hidden semi-Markov model is designed to allow modeling options to accommodate different types of genomic data and to serve as a general segmentation engine. By incorporating genomic positions into the sojourn distribution of HSMM, with optional prior learning using annotation or previous studies, the modeling output is more biologically sensible. The proposed model has been compared with several other state-of-the-art segmentation models through simulation benchmarking, which shows that our efficient implementation achieves comparable or better sensitivity and specificity in genomic segmentation.

Discovery of candidate genes for muscle traits based on GWAS supported by eQTL-analysis

Biochemical and biophysical processes that take place in muscle under relaxed and stressed conditions depend on the abundance and activity of gene products of metabolic and structural pathways. In livestock at post-mortem, these muscle properties determine aspects of meat quality and are measurable. The conversion of muscle to meat mimics pathological processes associated with muscle ischemia, injury or damage in humans and it is an economic factor in pork production. Linkage, association, and expression analyses independently contributed to the identification of trait-associated molecular pathways and genes. We aim at providing multiple evidences for the role of specific genes in meat quality by integrating a genome-wide association study (GWAS) for meat quality traits and the detection of eQTL based on trait-correlated expressed genes and trait-associated markers. The GWAS revealed 51 and 200 SNPs significantly associated with meat quality in a crossbred Pietrain×(German Landrace×Large White) (Pi×(GL×LW)) and a purebred German Landrace (GL) population, respectively. Most significant SNPs in Pi×(GL×LW) were located on chromosomes (SSC) 4 and 6. The data of 47,836 eQTLs at a significance level of p<10^-5 were used to scale down the number candidate genes located in these regions. These SNPs on SSC4 showed association with expression levels of ZNF704, IMPA1, and OXSR1; SSC6 SNPs were associated with expression of SIGLEC10 and PIH1D1. Most significant SNPs in GL were located on SSC6 and associated with expression levels of PIH1D1, SIGLEC10, TBCB, LOC100518735, KIF1B, LOC100514845, and two unknown genes. The abundance of transcripts of these genes in muscle, in turn, is significantly correlated with meat quality traits. We identified several genes with evidence for their candidacy for meat quality arising from the integrative approach of a genome-wide association study and eQTL analysis.

High- and low-protein gestation diets do not provoke common transcriptional responses representing universal target-pathways in muscle and liver of porcine progeny

AIM

Maternal diets introduce transcriptional changes in the offspring, highlighting the concept of genetic and physiological plasticity. Our previous analyses investigated stage-dependent transcriptional responses to either maternal high or low protein/carbohydrate ratios in either muscle or liver. Foetal programming is proposed to be mediated by a small number of gatekeeper processes, such as cytoskeleton remodelling and cell-cycle regulation. Here, we conducted an overall analysis of a three-dimensional data set aiming to elucidate, whether there are universally targeted pathways of adaptive transcriptional response to different protein/carbohydrate ratios.

METHODS

Microarray analyses were performed on liver and skeletal muscle tissue sampled at 94 days post-conception and 1, 28 and 188 days post-natum from offspring (n = 253) of German Landrace gilts that were fed isoenergetic diets containing low, high or adequate protein.

RESULTS

Cluster analyses revealed a hierarchical influence of tissue, ontogenetic stage and diet on transcript levels. Considering results cumulatively over stages, liver showed only marginal transcriptional differences between the dietary groups, whereas considerable differences appeared in muscle. Considering results cumulatively over tissues, nutrition-responsive transcriptions were observed along ontogenesis. Pathway analyses revealed transcript differences in genes related to tissue remodelling, cell-cycle regulation and mitochondrial function.

CONCLUSION

The factors tissue, stage and diet impact gene expression in a hierarchical order. Porcine liver appeared to be a tissue that was more resilient to nutritional modulation compared with skeletal muscle tissue. Differential modulation between tissues and dietary groups reveal that there are no universal target-pathways of adaptive transcriptional response to different protein diets.

PBMC transcription profiles of pigs with divergent humoral immune responses and lean growth performance

Background: The identification of key genes and regulatory networks in the transcriptomic responses of blood cells to antigen stimulation could facilitate the understanding of host defence and disease resistance. Moreover, genetic relationships between immunocompetence and the expression of other phenotypes, such as those of metabolic interest, are debated but incompletely understood in farm animals. Both positive and negative associations between immune responsiveness and performance traits such as weight gain or lean growth have been reported.

We designed an in vivo microarray study of transcriptional changes in porcine peripheral blood mononuclear cells (PBMCs) during the immune response to tetanus toxoid (TT) as a model antigen for combined cellular (Th1) and humoral (Th2) responses. The aim of the study was to investigate the responsiveness of PBMCs against the background of divergent lean growth (LG) performance and anti-TT antibody (AB) titers and to compare lean growth and humoral immune performance phenotypes.

Results: In general, high LG phenotypes had increased cellular immune response transcripts, while low AB phenotypes had increased transcripts for canonical pathways that represented processes of intracellular and second messenger signaling and immune responses. Comparison of lean growth phenotypes in the context of high AB titers revealed higher cellular immune response transcripts in high LG phenotypes. Similar comparisons in the context of low AB titers failed to identify any corresponding pathways. When high and low AB titer phenotypes were differentially compared, low AB phenotypes had higher cellular immune response transcripts on a low LG background and higher cell signaling, growth, and proliferation transcripts on a high LG background.

Conclusions: Divergent phenotypes of both lean growth performance and humoral immune response are affected by significant and functional transcript abundance changes throughout the immune response. The selected high-performance phenotypes demonstrated both high AB titers and increased transcript abundance of cellular immune response genes, which were possibly offset by lower expression of other cellular functions. Further, indications of compensatory effects were observed between cellular and humoral immune responses that became visible only in low-performance phenotypes.

Expression variation of the porcine ADRB2 has a complex genetic background

Porcine adrenergic receptor beta 2 (ADRB2) gene exhibits differential allelic expression in skeletal muscle, and its genetic variation has been associated with muscle pH. Exploring the molecular-genetic background of expression variation for porcine ADRB2 will provide insight into the mechanisms driving its regulatory divergence and may also contribute to unraveling the genetic basis of muscle-related traits in pigs. In the present study, we therefore examined haplotype effects on the expression of porcine ADRB2 in four tissues: longissimus dorsi muscle, liver, subcutaneous fat, and spleen. The diversity and structure of haplotypes of the proximal gene region segregating in German commercial breeds were characterized. Seven haplotypes falling into three clades were identified. Two clades including five haplotypes most likely originated from introgression of Asian genetics during formation of modern breeds. Expression analyses revealed that the Asian-derived haplotypes increase expression of the porcine ADRB2 compared to the major, wild-type haplotype independently of tissue type. In addition, several tissue-specific differences in the expression of the Asian-derived haplotypes were found. Inspection of haplotype sequences showed that differentially expressed haplotypes exhibit polymorphisms in a polyguanine tract located in the core promoter region. These findings demonstrate that expression variation of the porcine ADRB2 has a complex genetic basis and suggest that the promoter polyguanine tract is causally involved. This study highlights the challenges of finding causal genetic variants underlying complex traits.

Gene expression profiling of articular cartilage reveals functional pathways and networks of candidate genes for osteochondrosis in pigs

Osteochondrosis (OC) is a joint disorder that frequently causes leg weakness in growing pigs, resulting in welfare problems and economic losses. We aimed to detect molecular pathways relevant to the emergence of the disease and to identify candidate genes for the liability to the disorder. Therefore, we compared microarray-based expression patterns of articular cartilage with (n=11) and without (n=11) histologically diagnosed OC lesions obtained from discordant sib-pairs. A total of 1,564 genes were found with different transcript abundance [differentially expressed (DE) genes] at q≤0.05. To further identify candidate genes, we integrated data from quantitative trait loci (QTL) and genome-wide association (GWA) studies with the expression analysis. We detected 317 DE genes within the QTL confidence intervals, of which 26 DE genes also overlapped GWA regions. Ingenuity Pathway Analysis suggests a pathogenic role of immune response, angiogenesis, and synthesis of extracellular matrix pathways for OC. These processes could facilitate the emergence of defects. But they may also promote the degradation of articular cartilage and the worsening of the disease. A functional network was derived that comprised genes with functional and positional clues of their role in bone and cartilage metabolisms and development, including extracellular matrix genes (e.g., LOX, OGN, and ASPN), angiogenesis genes (e.g., ANGPTL4 and PDGFA), and immune response genes (e.g., ICAM1, AZGP1, C1QB, C1QC, PDE4B, and CDA). The study identified molecular processes linked to OC and several genes with positional, genetic-statistical, and functional evidence for their role in the emergence of articular cartilage lesions and the liability to OC.

Correlated mRNAs and miRNAs from co-expression and regulatory networks affect porcine muscle and finally meat properties

Background

Physiological processes aiding the conversion of muscle to meat involve many genes associated with muscle structure and metabolic processes. MicroRNAs regulate networks of genes to orchestrate cellular functions, in turn regulating phenotypes.

Results

We applied weighted gene co-expression network analysis to identify co-expression modules that correlated to meat quality phenotypes and were highly enriched for genes involved in glucose metabolism, response to wounding, mitochondrial ribosome, mitochondrion, and extracellular matrix. Negative correlation of miRNA with mRNA and target prediction were used to select transcripts out of the modules of trait-associated mRNAs to further identify those genes that are correlated with post mortem traits.

Conclusions

Porcine muscle co-expression transcript networks that correlated to post mortem traits were identified. The integration of miRNA and mRNA expression analyses, as well as network analysis, enabled us to interpret the differentially-regulated genes from a systems perspective. Linking co-expression networks of transcripts and hierarchically organized pairs of miRNAs and mRNAs to meat properties yields new insight into several biological pathways underlying phenotype differences. These pathways may also be diagnostic for many myopathies, which are accompanied by deficient nutrient and oxygen supply of muscle fibers.

Transcript variants of the porcine glucocorticoid receptor gene (NR3C1)

Glucocorticoid receptor (GR) is a transcription factor activated by circulating glucocorticoids and mediates their effects on various biological functions in the body. The final cellular activity of GR is modulated by alternative splicing and cell-type specific expression of its encoding gene NR3C1. To enhance the current knowledge of alternative processing of NR3C1 in mammalian species and to facilitate future studies of its regulation in the pig we explored here structural diversity, and tissue-specific distribution of transcript variants of the porcine NR3C1, and the correlation between usage of alternative promoters and alternative splicing. We experimentally identified ten alternatively used untranslated first exons (1A, 1B, 1C, 1D, 1E, 1F, 1G, 1H, 1H2 and 1J) and four transcript variants encoding different GR subtypes (GR-alpha, GR-beta, GR-P and GR-gamma). Expression profiling of nine most important target tissues of glucocorticoids revealed that the promoter of exon 1C drives constitutive expression of the predominant GR-alpha subtype. We found compelling evidence that the occurrence of exon 1D influences abundance of the GR-P splice variant, while both seems to play an important role in regulating GR activity in neuroendocrine tissues. Exons 1A and 1B in turn appear to be important for the regulation of the expression of the porcine NR3C1 in liver and spleen. Our results demonstrate that tissue-specific actions of GR depend on the usage of alternative promoter regions that favour the processing of certain GR subtypes.