Subcutaneous fat (SCF) and intramuscular fat (IMF) deposition is relevant to health in humans, as well as meat production and quality in pigs. In this study, we generated RNA sequence data for 122 ...SCF, 120 IMF, and 87 longissimus dorsi muscle (LDM) samples using 155 F
pigs from a specially designed heterogeneous population generated by intercrossing four highly selected European commercial breeds and four indigenous Chinese pig breeds. The phenotypes including waist back fat thickness and intramuscular fat content were also measured in the 155 F
pigs. We found that the genes in SCF and IMF differed largely in both expression levels and network connectivity, and highlighted network modules that exhibited strongest gain of connectivity in SCF and IMF, containing genes that were associated with the immune process and DNA double-strand repair, respectively. We identified 215 SCF genes related to kinase inhibitor activity, mitochondrial fission, and angiogenesis, and 90 IMF genes related to lipolysis and fat cell differentiation, displayed a tissue-specific association with back fat thickness and IMF content, respectively. We found that cis-expression QTL for trait-associated genes in the two adipose tissues tended to have tissue-dependent predictability for the two adipose traits. Alternative splicing of genes was also found to be associated with SCF or IMF deposition, but the association was much less extensive than that based on expression levels. This study provides a better understanding of SCF and IMF gene transcription and network organization and identified critical genes and network modules that displayed tissue-specific associations with subcutaneous and intramuscular fat deposition. These features are helpful for designing breeding programs to genetically improve the two adipose traits in a balanced way.
Full text
Available for:
EMUNI, FIS, FZAB, GEOZS, GIS, IJS, IMTLJ, KILJ, KISLJ, MFDPS, NLZOH, NUK, OILJ, PNG, SAZU, SBCE, SBJE, SBMB, SBNM, UKNU, UL, UM, UPUK, VKSCE, ZAGLJ
DnaJ heat shock protein family (Hsp40) member A3 (DNAJA3) plays an important role in viral infections. However, the role of DNAJA3 in replication of foot-and-mouth-disease virus (FMDV) remains ...unknown. In this study, DNAJA3, a novel binding partner of VP1, was identified using yeast two-hybrid screening. The DNAJA3-VP1 interaction was further confirmed by coimmunoprecipitation and colocalization in FMDV-infected cells. The J domain of DNAJA3 (amino acids 1 to 168) and the lysine at position 208 (K208) of VP1 were shown to be critical for the DNAJA3-VP1 interaction. Overexpression of DNAJA3 dramatically dampened FMDV replication, whereas loss of function of DNAJA3 elicited opposing effects against FMDV replication. Mechanistical study demonstrated that K208 of VP1 was critical for reducing virus titer caused by DNAJA3 using K208A mutant virus. DNAJA3 induced lysosomal degradation of VP1 by interacting with LC3 to enhance the activation of lysosomal pathway. Meanwhile, we discovered that VP1 suppressed the beta interferon (IFN-β) signaling pathway by inhibiting the phosphorylation, dimerization, and nuclear translocation of IRF3. This inhibitory effect was considerably boosted in DNAJA3-knockout cells. In contrast, overexpression of DNAJA3 markedly attenuated VP1-mediated suppression on the IFN-β signaling pathway. Poly(I⋅C)-induced phosphorylation of IRF3 was also decreased in DNAJA3-knockout cells compared to that in the DNAJA3-WT cells. In conclusion, our study described a novel role for DNAJA3 in the host's antiviral response by inducing the lysosomal degradation of VP1 and attenuating the VP1-induced suppressive effect on the IFN-β signaling pathway.
This study pioneeringly determined the antiviral role of DNAJA3 in FMDV. DNAJA3 was found to interact with FMDV VP1 and trigger its degradation via the lysosomal pathway. In addition, this study is also the first to clarify the mechanism by which VP1 suppressed IFN-β signaling pathway by inhibiting the phosphorylation, dimerization, and nuclear translocation of IRF3. Moreover, DNAJA3 significantly abrogated VP1-induced inhibitive effect on the IFN-β signaling pathway. These data suggested that DNAJA3 plays an important antiviral role against FMDV by both degrading VP1 and restoring of IFN-β signaling pathway.
Glycolytic potential (GP) in skeletal muscle is economically important in the pig industry because of its effect on pork processing yield. We have previously mapped a major quantitative trait loci ...(QTL) for GP on chromosome 3 in a White Duroc × Erhualian F2 intercross. We herein performed a systems genetic analysis to identify the causal variant underlying the phenotype QTL (pQTL). We first conducted genome-wide association analyses in the F2 intercross and an F19 Sutai pig population. The QTL was then refined to an 180-kb interval based on the 2-LOD drop method. We then performed expression QTL (eQTL) mapping using muscle transcriptome data from 497 F2 animals. Within the QTL interval, only one gene (PHKG1) has a cis-eQTL that was colocolizated with pQTL peaked at the same SNP. The PHKG1 gene encodes a catalytic subunit of the phosphorylase kinase (PhK), which functions in the cascade activation of glycogen breakdown. Deep sequencing of PHKG1 revealed a point mutation (C>A) in a splice acceptor site of intron 9, resulting in a 32-bp deletion in the open reading frame and generating a premature stop codon. The aberrant transcript induces nonsense-mediated decay, leading to lower protein level and weaker enzymatic activity in affected animals. The mutation causes an increase of 43% in GP and a decrease of>20% in water-holding capacity of pork. These effects were consistent across the F2 and Sutai populations, as well as Duroc × (Landrace × Yorkshire) hybrid pigs. The unfavorable allele exists predominantly in Duroc-derived pigs. The findings provide new insights into understanding risk factors affecting glucose metabolism, and would greatly contribute to the genetic improvement of meat quality in Duroc related pigs.
Full text
Available for:
DOBA, IZUM, KILJ, NUK, PILJ, PNG, SAZU, SIK, UILJ, UKNU, UL, UM, UPUK
Carcass length is very important for body size and meat production for swine, thus understanding the genetic mechanisms that underly this trait is of great significance in genetic improvement ...programs for pigs. Although many quantitative trait loci (QTL) have been detected in pigs, very few have been fine-mapped to the level of the causal mutations. The aim of this study was to identify potential causal single nucleotide polymorphisms (SNPs) for carcass length by integrating a genome-wide association study (GWAS) and functional assays.
Here, we present a GWAS in a commercial Duroc × (Landrace × Yorkshire) (DLY) population that reveals a prominent association signal (P = 4.49E-07) on pig chromosome 17 for carcass length, which was further validated in two other DLY populations. Within the detected 1 Mb region, the BMP2 gene stood out as the most likely causal candidate because of its functions in bone growth and development. Whole-genome gene expression studies showed that the BMP2 gene was differentially expressed in the cartilage tissues of pigs with extreme carcass length. Then, we genotyped an additional 267 SNPs in 500 selected DLY pigs, followed by further whole-genome SNP imputation, combined with deep genome resequencing data on multiple pig breeds. Reassociation analyses using genotyped and imputed SNP data revealed that the rs320706814 SNP, located approximately 123 kb upstream of the BMP2 gene, was the strongest candidate causal mutation, with a large association with carcass length, with a ~ 4.2 cm difference in length across all three DLY populations (N = 1501; P = 3.66E-29). This SNP segregated in all parental lines of the DLY (Duroc, Large White and Landrace) and was also associated with a significant effect on body length in 299 pure Yorkshire pigs (P = 9.2E-4), which indicates that it has a major value for commercial breeding. Functional assays showed that this SNP is likely located within an enhancer and may affect the binding affinity of transcription factors, thereby regulating BMP2 gene expression.
Taken together, these results suggest that the rs320706814 SNP on pig chromosome 17 is a putative causal mutation for carcass length in the widely used DLY pigs and has great value in breeding for body size in pigs.
The negative regulation of antiviral immune responses is essential for the host to maintain homeostasis. Jumonji domain-containing protein 6 (JMJD6) was previously identified with a number of ...functions during RNA virus infection. Upon viral RNA recognition, retinoic acid-inducible gene-I-like receptors (RLRs) physically interact with the mitochondrial antiviral signaling protein (MAVS) and activate TANK-binding kinase 1 (TBK1) to induce type-I interferon (IFN-I) production. Here, JMJD6 was demonstrated to reduce type-I interferon (IFN-I) production in response to cytosolic poly (I:C) and RNA virus infections, including Sendai virus (SeV) and Vesicular stomatitis virus (VSV). Genetic inactivation of JMJD6 enhanced IFN-I production and impaired viral replication. Our unbiased proteomic screen demonstrated JMJD6 contributes to IRF3 K48 ubiquitination degradation in an RNF5-dependent manner. Mice with gene deletion of JMJD6 through piggyBac transposon-mediated gene transfer showed increased VSV-triggered IFN-I production and reduced susceptibility to the virus. These findings classify JMJD6 as a negative regulator of the host's innate immune responses to cytosolic viral RNA.
Full text
Available for:
DOBA, IZUM, KILJ, NUK, PILJ, PNG, SAZU, SIK, UILJ, UKNU, UL, UM, UPUK
Copy number variation (CNV) is a major source of structural variants and has been commonly identified in mammalian genome. It is associated with gene expression and may present a major genetic ...component of phenotypic diversity. Unlike many other mammalian genomes where CNVs have been well annotated, studies of porcine CNV in diverse breeds are still limited.
Here we used Porcine SNP60 BeadChip and PennCNV algorithm to identify 1,315 putative CNVs belonging to 565 CNV regions (CNVRs) in 1,693 pigs from 18 diverse populations. Total 538 out of 683 CNVs identified in a White Duroc × Erhualian F2 population fit Mendelian transmission and 6 out of 7 randomly selected CNVRs were confirmed by quantitative real time PCR. CNVRs were non-randomly distributed in the pig genome. Several CNV hotspots were found on pig chromosomes 6, 11, 13, 14 and 17. CNV numbers differ greatly among different pig populations. The Duroc pigs were identified to have the most number of CNVs per individual. Among 1,765 transcripts located within the CNVRs, 634 genes have been reported to be copy number variable genes in the human genome. By integrating analysis of QTL mapping, CNVRs and the description of phenotypes in knockout mice, we identified 7 copy number variable genes as candidate genes for phenotypes related to carcass length, backfat thickness, abdominal fat weight, length of scapular, intermuscle fat content of logissimus muscle, body weight at 240 day, glycolytic potential of logissimus muscle, mean corpuscular hemoglobin, mean corpuscular volume and humerus diameter.
We revealed the distribution of the unprecedented number of 565 CNVRs in pig genome and investigated copy number variable genes as the possible candidate genes for phenotypic traits. These findings give novel insights into porcine CNVs and provide resources to facilitate the identification of trait-related CNVs.
Full text
Available for:
DOBA, IZUM, KILJ, NUK, PILJ, PNG, SAZU, SIK, UILJ, UKNU, UL, UM, UPUK
Sequencing-based genome-wide association studies (GWAS) have facilitated the identification of causal associations between genetic variants and traits in diverse species. However, it is ...cost-prohibitive for the majority of research groups to sequence a large number of samples. Here, we carried out genotype imputation to increase the density of single nucleotide polymorphisms in a large-scale Swine F
population using a reference panel including 117 individuals, followed by a series of GWAS analyses. The imputation accuracies reached 0.89 and 0.86 for allelic concordance and correlation, respectively. A quantitative trait nucleotide (QTN) affecting the chest vertebrate was detected directly, while the investigation of another QTN affecting the residual glucose failed due to the presence of similar haplotypes carrying wild-type and mutant allelesin the reference panel used in this study. A high imputation accuracy was confirmed by Sanger sequencing technology for the most significant loci. Two candidate genes, CPNE5 and MYH3, affecting meat-related traits were proposed. Collectively, we illustrated four scenarios in imputation-based GWAS that may be encountered by researchers, and our results will provide an extensive reference for future genotype imputation-based GWAS analyses in the future.
Full text
Available for:
EMUNI, FIS, FZAB, GEOZS, GIS, IJS, IMTLJ, KILJ, KISLJ, MFDPS, NLZOH, NUK, OILJ, PNG, SAZU, SBCE, SBJE, SBMB, SBNM, UKNU, UL, UM, UPUK, VKSCE, ZAGLJ
Thousands of QTLs for meat quality traits have been identified by linkage mapping studies, but most of them lack precise position or replication between populations, which hinder their application in ...pig breeding programs. To localize QTLs for meat quality traits to precise genomic regions, we performed a genome-wide association (GWA) study using the Illumina PorcineSNP60K Beadchip in two swine populations: 434 Sutai pigs and 933 F2 pigs from a White Duroc×Erhualian intercross. Meat quality traits, including pH, color, drip loss, moisture content, protein content and intramuscular fat content (IMF), marbling and firmness scores in the M. longissimus (LM) and M. semimembranosus (SM) muscles, were recorded on the two populations. In total, 127 chromosome-wide significant SNPs for these traits were identified. Among them, 11 SNPs reached genome-wise significance level, including 1 on SSC3 for pH, 1 on SSC3 and 3 on SSC15 for drip loss, 3 (unmapped) for color a*, and 2 for IMF each on SSC9 and SSCX. Except for 11 unmapped SNPs, 116 significant SNPs fell into 28 genomic regions of approximately 10 Mb or less. Most of these regions corresponded to previously reported QTL regions and spanned smaller intervals than before. The loci on SSC3 and SSC7 appeared to have pleiotropic effects on several related traits. Besides them, a few QTL signals were replicated between the two populations. Further, we identified thirteen new candidate genes for IMF, marbling and firmness, on the basis of their positions, functional annotations and reported expression patterns. The findings will contribute to further identification of the causal mutation underlying these QTLs and future marker-assisted selection in pigs.
Full text
Available for:
DOBA, IZUM, KILJ, NUK, PILJ, PNG, SAZU, SIK, UILJ, UKNU, UL, UM, UPUK
Muscle glycolytic potential (GP) is a key factor affecting multiple meat quality traits. It is calculated based on the contents of residual glycogen and glucose (RG), glucose-6-phosphate (G6P), and ...lactate (LAT) contents in muscle. However, the genetic mechanism of glycolytic metabolism in skeletal muscle of pigs remains poorly understood. With a history of more than 400 years and some unique characteristics, the Erhualian pig is called the "giant panda" (very precious) in the world's pig species by Chinese animal husbandry.
Here, we performed a genome-wide association study (GWAS) using 1.4M single nucleotide polymorphisms (SNPs) chips for longissimus RG, G6P, LAT, and GP levels in 301 purebred Erhualian pigs.
We found that the average GP value of Erhualian was unusually low (68.09 μmol/g), but the variation was large (10.4-112.7 μmol/g). The SNP-based heritability estimates for the four traits ranged from 0.16-0.32. In total, our GWAS revealed 31 quantitative trait loci (QTLs), including eight for RG, nine for G6P, nine for LAT, five for GP. Of these loci, eight were genome-wide significant (
< 3.8 × 10
), and six loci were common to two or three traits. Multiple promising candidate genes such as
,
,
,
,
and
were identified. The genotype combinations of the five GP-associated SNPs also showed significant effect on other meat quality traits.
These results not only provide insights into the genetic architecture of GP related traits in Erhualian, but also are useful for pig breeding programs involving this breed.
Elucidation of the pig transcriptome is essential for interpreting functional elements of the genome and understanding the genetic architecture of complex traits such as fat deposition, metabolism ...and growth.
Here we used massive parallel high-throughput RNA sequencing to generate a high-resolution map of the porcine mRNA and miRNA transcriptome in liver, longissimus dorsi and abdominal fat from two full-sib F2 hybrid pigs with segregated phenotypes on growth, blood physiological and biochemical parameters, and fat deposition. We obtained 8,508,418-10,219,332 uniquely mapped reads that covered 78.0% of the current annotated transcripts and identified 48,045-122,931 novel transcript fragments, which constituted 17,085-29,499 novel transcriptional active regions in six tested samples. We found that about 18.8% of the annotated genes showed alternative splicing patterns, and alternative 3' splicing is the most common type of alternative splicing events in pigs. Cross-tissue comparison revealed that many transcriptional events are tissue-differential and related to important biological functions in their corresponding tissues. We also detected a total of 164 potential novel miRNAs, most of which were tissue-specifically identified. Integrated analysis of genome-wide association study and differential gene expression revealed interesting candidate genes for complex traits, such as IGF2, CYP1A1, CKM and CES1 for heart weight, hemoglobin, pork pH value and serum cholesterol, respectively.
This study provides a global view of the complexity of the pig transcriptome, and gives an extensive new knowledge about alternative splicing, gene boundaries and miRNAs in pigs. Integrated analysis of genome wide association study and differential gene expression allows us to find important candidate genes for porcine complex traits.
Full text
Available for:
DOBA, IZUM, KILJ, NUK, PILJ, PNG, SAZU, SIK, UILJ, UKNU, UL, UM, UPUK