Exploring transcriptomic diversity in muscle revealed that cellular signaling pathways mainly differentiate five Western porcine breeds
© SanCristobal et al. 2015
Received: 26 May 2015
Accepted: 30 November 2015
Published: 12 December 2015
Among transcriptomic studies, those comparing species or populations can increase our understanding of the impact of the evolutionary forces on the differentiation of populations. A particular situation is the one of short evolution time with breeds of a domesticated species that underwent strong selective pressures. In this study, the gene expression diversity across five pig breeds has been explored in muscle. Samples came from: 24 Duroc, 33 Landrace, 41 Large White dam line, 10 Large White sire line and 39 Piétrain. From these animals, 147 muscle samples obtained at slaughter were analyzed using the porcine Agilent 44 K v1 microarray.
A total of 12,358 genes were identified as expressed in muscle after normalization and 1,703 genes were declared differential for at least one breed (FDR < 0.001). The functional analysis highlighted that gene expression diversity is mainly linked to cellular signaling pathways such as the PI3K (phosphoinositide 3-kinase) pathway. The PI3K pathway is known to be involved in the control of development of the skeletal muscle mass by affecting extracellular matrix - receptor interactions, regulation of actin cytoskeleton pathways and some metabolic functions. This study also highlighted 228 spots (171 unique genes) that differentiate the breeds from each other. A common subgroup of 15 genes selected by three statistical methods was able to differentiate Duroc, Large White and Piétrain breeds.
This study on transcriptomic differentiation across Western pig breeds highlighted a global picture: mainly signaling pathways were affected. This result is consistent with the selection objective of increasing muscle mass. These transcriptional changes may indicate selection pressure or simply breed differences which may be driven by human selection. Further work aiming at comparing genetic and transcriptomic diversities would further increase our understanding of the consequences of human impact on livestock species.
After the study of genome evolution, much effort has been recently put on the evolution of gene regulation and gene expression (see Romero et al.  for a review). Several recent studies report the existence of differentially expressed genes across species. These differences in gene expression could be due to various evolutionary processes, neutral or not ([2–13]). All these studies focused on comparisons between species, with large evolutionary divergences.
For shorter evolutionary times, Hufford et al.  observed that candidate genes for domestication in maize do not display any specific expression profile, contrary to the genomic patterns (DNA sequence). For the domesticated period of sorghum, Jiang et al.  observed that gene expression divergence between two lines was mainly determined by DNA sequence divergence. Nätt et al.  on the contrary observed numerous gene expression and methylation changes between wild and domesticated chickens, with an overrepresentation in selective sweeps. For a shorter time scale, Yang et al.  used gene expression in addition to genomic polymorphism (SNPs) to assign a human individual to its ethnic population. Muller et al.  suggest that gene expression changes could be related to the out-of-Africa adaptation in Drosophila with a clear sex-specificity. In budding yeast, Fraser et al.  observed that entire pathways can be affected by adaptation of gene expression. Transcriptomic differences, as well as proteomics and metabolomics were shown between two genetically diverse dry bean germplasm by Mensack et al. .
As briefly shown above, some work has been conducted on the differences in gene expression across populations, these populations being mostly at species level. The following question is however less studied. Are there any differentially expressed genes across breeds of the same species? If yes, what are their biological functions? We propose to explore this problem in the particular context of a domesticated species with a long history of selection pressure from human: the pig. Perez-Enciso et al.  compared gene expression among pig breeds in several tissues and provided interesting results on gene expression divergence. However, only 16 animals were used in their study. The aim of our study was to see if genes were differentially expressed among the main Western pig breeds, at a larger sampling scale, and what might be the biological implications. Secondly, we aim to answer the questions what kind of gene expression characterizes a particular breed, and what differentiates a breed on the basis of gene expression. For that purpose, we will focus on one tissue of interest in pig breeding: a post-mortem muscle. Indeed pig meat production has placed a strong selective pressure on various characteristics of meat (muscle).
The sampled animals came from 5 pig breeds: 24 Duroc (DU), 33 Landrace (LR), 41 Large White dam line (LWF), 10 Large White sire line (LWM) and 39 Piétrain (PI). The two Large White lines derived from a common and recent ancestor, and were specialized on “male” traits (like conformation) for LWM or on “female” traits (like maternal behaviour) for LWF. All animals were males (castrates) except Piétrain (females). They were reared in 10 different contemporary groups, and slaughtered around 100 kg as in . From these animals, 147 post-mortem Longissimus dorsi muscle samples were analyzed using the porcine Agilent 44 K v1 microarray. A total of 12,358 probes were detected and considered as expressed in muscle in the conditions of this experiment and after normalization. From these 12,358 probes, 9,055 (73 %) were mapped on the porcine genome (Sscrofa10.2 assembly) of which 8,758 (71 %) were localized on the autosomes.
Global differential analysis
Number of differentially expressed probes linked to the breed effect
FDR 0.1 %
Breed effect (global)
Global or pairwise
Global and pairwise
Number of differentially expressed probes between two breeds from the pairwise analysis
FDR 0.1 %
The overall intersect for comparisons with Duroc included the v-kit Hardy-Zuckerman 4 feline sarcoma viral oncogene homolog (KIT)  and oculocutaneous albinism II (OCA2) genes (among 19 probes in total corresponding to 17 unique genes and two unknown). KIT and OCA2 genes are both mainly known to be involved in the determination of skin colour . As expected, both genes are under-expressed in Duroc compared to the other pigs, white pigs or white spotted pigs (Fig. 3). This under-expression was observed in muscle tissue where KIT and OCA2 encode membrane transporters and KIT is involved in the activation of the PI3K pathway.
Twenty-two genes were differentially expressed between the Large White dam line and the other breeds excluding the Large White sire line. The most differentially expressed gene is RAB18, a gene coding a small G proteins belonging to the Ras superfamily. This superfamily is known to coordinate vesicular trafficking in the cell . Mutations in RAB18 have been detected in humans with progressive neurological deterioration, and a severe hypotonia. Moreover, the variance of RAB18 expression in skeletal muscle may suggest that this gene is especially genetically regulated in Large White dam line (Fig. 3).
Among the 37 genes differentially expressed between Piétrain and all the other breeds, 11 are localized on chromosomes X or Y. Two examples of gene expression are presented on Fig. 3 with the genes USP9X and USP9Y, respectively up and down regulated in Piétrain. This result corresponds to an over-representation of expression in females compared to males while a more widely held view is dosage compensation between XX and XY gene expression. In this work, as expected the USP9 gene localized on Y chromosome was not expressed in female; while the expression of the USP9 gene localized on X chromosome was overexpressed in Piétrain compared to its expression in the other breeds. Here where only Piétrains are female, we observed a doubled expression of the USP9X and no expression of USP9Y. The result would be an almost equal expression of USP9 proteins with probably similar function (assuming that both genes have exactly the same spatio-temporal expression profiles).
Differential analysis restricted to the genes localized on autosomes
According to our preceding results and the description of how chromosome X is subject to selective pressure with more highly sex-biased gene expression , another differential analysis was conducted on the restricted list of 8,758 genes located on autosomes in order to lessen the sex effect observed with the first analysis. In these conditions, 1,228 DE probes were identified globally across breeds (FDR 0.1 %) and 1,232 DE probes for the pairwise comparisons (FDR 0.1 %). The union of both lists led to a total of 1,374 DE probes with a FDR for the global analysis of less than 0.8 %. These analyses (all probes and the ones restricted to autosomes) were done in parallel. That would allow the identification of the molecular basis of the difference between breeds including the Piétrain. However some genes (e.g. IGF2) were not localized because they are absent from the current assembly, even if their locations are well known (e.g. IGF2 is on SSC2). Figure 1 (c and d) shows how avoiding X localized genes placed Piétrain pigs closer to the other breeds without changing the overall projection of the different breeds. The number of DE probes can be found in Tables 1 and 2. See Additional file 3 for details about probes on autosomes. No localized DE probes were identified between Large White sire and dam lines even with a FDR of 5 %.
Restricting the analysis to the autosomes is likely to be only a partial solution to the confounding effect of gender and breed. We chose to maintain Piétrain in the analysis because this breed is of prime importance for the breeders. Indeed, it has the highest muscularity and up to 80 % of the semen used in French production comes from Piétrain breed. Moreover, the gender effect may be reduced by the fact that all males were castrated in this work. We did the analysis without the Piétrain breed and the lists of DE probes were similar to those obtained with the Piétrain breed (not shown). The functional enrichment is hence similar for pairwise comparisons involving all breeds except the Piétrain.
Tables 1 and 2 underline the large number of gene lists to be analysed for biological interpretation. The functional analysis was undertaken in a sequential manner. First, the objective was to evaluate the impact of the sex effect of the Piétrain breed. While it would have been possible to restrict the functional analysis to autosomal genes, excluding important genes such as IGF2 (absent from the current assembly) may give unreliable results. Therefore, the list of the 1,374 DE probes from the union of the global and pairwise comparisons restricted to autosomes was compared to the same list plus the 61 DE probes localised on chromosomes X and Y (Additional file 2). The free GeneCodis software  was used and the comparison was restricted to KEGG pathways, which are a collection of manually drawn pathway maps representing our knowledge on the molecular interaction and reaction networks. The top 10 significant pathways were exactly the same between both lists with or without genes on sex chromosomes (data not shown). According to the weak effect of the analysis restricted to autosomes on the PCA projection and on functional analysis, we hypothesize a low impact of the sex effect on our breed comparison, even if this effect could not be absolutely excluded.
The lists of DE genes from the union of global and pairwise analysis restricted (1,374 DE spots) or not (1,858 DE spots) to autosomes were compared through GeneCodis.
Relevant and significantly enriched biological functions (KEGG pathways) for diversity of muscle expression for five pig breeds
Number of genes
ITGA1, MAPK1, comP, ITGA6, ERBB2, LAMA2, CCND2, BCAR1, DIAPH1, PTEN, ITGA2, SPP1, ACTB, LAMC1, CAPN2, COL4A1, EGFR, VAV2, THBS4, ITGA9, THBS2, PTK2, PIK3CG, VEGFA, fn1, PPP1CB, VAV3
Regulation of actin cytoskeleton
PIP4K2A, ITGA1, MAPK1, ARPC2, ITGA6, ITGB2, ARPC1B, F2R, BCAR1, DIAPH1, GNA12, ITGA2, NCKAP1L, GSN, ACTB, TMSB4X, EGFR, VAV2, ITGA9, CD14, PTK2, Kras, BAIAP2, PIK3CG, fn1, PPP1CB, VAV3
ITGA1, comP, ITGA6, LAMA2, ITGA2, GP1BB, SPP1, LAMC1, COL4A1, THBS4, ITGA9, THBS2, AGRN, fn1
MAPK1, PARD3, CTNNA3, ACP1, ERBB2, PTPRJ, ACTB, EGFR, BAIAP2, CREBBP, SORBS1, PVRL3
Valine, leucine and isoleucine degradation
ALDH3A2, BCAT1, OXCT1, EHHADH, IVD, BCKDHB, ECHS1, ACADSB, HADHA, PCCB, AUH, ALDH2
MLL5, ALDH3A2, EHHADH, NSD1, WHSC1L1, ECHS1, HADHA, setmar, ALDH2
NUDT2, Pgm2, GMPR2, POLR2G, ENPP1, RRM2B, FHIT, GUCY1A3, ADA, NT5C2, NME6, PGM1, PRPS2, ATIC, NT5E, POLR3E
ALDH3A2, EHHADH, ECHS1, CARNS1, HADHA, ALDH2
Arginine and proline metabolism
ALDH3A2, ACY1, P4HA1, GLUD1, CARNS1, AMD1, ALDH4A1, ALDH2
Fatty acid metabolism
ALDH3A2, EHHADH, ECHS1, ACADSB, ACSL4, HADHA, ALDH2
comP, CTSS, ATP6V1C1, ITGB2, ITGA2, PIK3C3, ACTB, CYBB, CYBA, TFRC, THBS4, CD14, THBS2, NCF4, NCF2, DYNC1I2, EEA1, FCGR2B, MRC1, TUBA8, TUBA4A
Protein processing in endoplasmic reticulum
nsfl1c, BAG2, UBE4B, PARK2, SAR1B, HSPH1, SEC24D, DNAJC1, CAPN2, DDOST, BCAP31, HSPA2, PDIA3, CUL1, UBE2D4, HERPUD1, EDEM2, DNAJA1, DERL1, CKAP4, TUSC3
mRNA surveillance pathway
SAP18, PPP2R5A, CPSF6, CLP1, PPP2R2A, PAPOLA, PAPOLB, PPP2CB, Fip1l1, ETF1, HBS1L, CPSF3, CSTF2T
RAB11FIP5, PARD3, KIT, F2R, DNM2, EHD1, RABEP1, ARRB2, GRK5, DNM1L, FOLR1, EGFR, TFRC, HSPA2, WWP1, ZFYVE16, SMAP2, EEA1, epn1, FAM125B
A differential expression analysis aims at giving a list of genes whose transcript abundances differ significantly among classes (here breeds). A discriminant analysis aims at identifying transcripts whose abundance differences help clustering (separating) each breed from the others. It helps predicting the breed to which a RNA sample belongs to, based on a list of “discriminant” genes. Although related, these 2 notions are different. In general a discriminant gene is a differential gene, but the reverse is not always true. Among a vast choice of discriminant methods, we performed a Random Forest (RF) analysis [28, 29] which is robust and non-linear and a sparse Partial Least Square – Discriminant Analysis (sPLS-DA)  that is pertaining to linear discriminators.
Top 15 discriminant genes for breeds
v-kit Hardy-zuckerman 4 fe line sarcoma viral oncogene homolog
axis 2 (LW>DU,LR)
Endocytosis, cell growth and/or maintenance
ER membrane protein complex sub unit 2
axis 2 (LW>DU,LR)
Component of the ER membrane protein complex (EMC)
tight junction protein 2 (zona occludens 2)
axis 3 (DU>LR)
proteasome (prosome, macropain)sub unit, beta type, 4
axis 2/3 (DU>LR,LW,PI)
vitamin D (1,25-dihydroxyvitamin D3) receptor
axis 2/3 (DU<LR, LW, PI)
oculocutaneous albinism II (pink-eye dilution homolog, mouse)
axis 2/3 (DU<LR, LW, PI)
Transport of tyrosine
RAB18, member RAS oncogene family
axis 2 (LW<DU,LR)
Two-component signal transduction system (phosphorelay)
protein kinase, cAMP-dependent, catalytic, beta
axis 2 (LW<DU,LR)
Protein amino acid phosphorylation; signal transduction; Insulin signaling pathway, Gap
tRNA-yW synthesizing protein 3 homolog (S. cerevisiae)
axis 2 (LW>DU,LR)
RNA translation, tRNA stabilization
granzyme B (granzyme 2, cytotoxic T-lymphocyte-associated serine esterase 1)
axis 2 (LW<DU,LR)
Apoptosis; cytolysis; proteolysis and peptidolysis
N(alpha)-acetyltransferase 20, NatBcatalytic subunit
axis 1 (PI<DU,LR,LW)
CotranslationalN(alpha)-terminal acetylation of methionine residues
MOB kinase activation3C
axis 1 (PI<DU,LR,LW)
Protein kinase essential for spindle pole body duplication and mitotic checkpoint regulation
phosphoinositide-3-kinase, catalytic, gamma polypeptide
axis 1 (PI>DU,LR,LW)
G-protein coupled receptorprotein signaling pathway; Regulation of actin cytoskeleton
RNA guanylyltransferase and 5’-phosphatase
axis 1 (PI<DU,LR,LW)
mRNA capping; protein amino acid dephosphorylation
diphthamide biosynthesis 5
axis 1 (PI<DU,LR,LW)
Diphthamide synthesis pathway
The present work highlights the variability of the expression of about 1,800 transcripts between five pig breeds. These genes were identified for global differences (using a linear mixed model, 1,703 transcripts), and a pairwise analysis (1,655 transcripts). A subset of 228 transcripts was identified to be able to discriminate the five breeds in one way or another. A short list of 15, as the intersection of all lists, was extracted. Altogether, the union of all these highlighted transcripts corresponded to 1,230 unique genes.
Variability of expression between breeds mainly concerns PI3Kinase signaling pathways involved in the regulation of muscle mass
Muscle mass is determined by both the number of muscle fibres and the size of these fibres. The development of muscle is temporally regulated and the total fibre number is fixed before birth (around day 90 of gestation, birth is around day 114) in pigs . Afterwards, increase in muscle mass may be the consequence of mechanisms regulating muscle hypertrophy like the size of the fibres. In our study, muscle samples were collected at slaughter. Animals were about six months old and weighted around 110 kg. Therefore the expression analysis corresponds probably more to the hypertrophic process if existing. This hypothetical statement is in accordance with the results of the enrichment analysis that highlighted signaling pathways most regulated by PI3 kinase and known to be involved in the regulation of muscle mass. Then transcriptomic diversity between the five breeds underlined how these pathways regulating muscle mass may have been targeted by selection to increase lean meat production.
Upstream to the PI3K signaling pathway, regulated DE genes referred to interactions between cells or with the extracellular matrix (ECM) together with focal adhesion. These pathways are the three among the four top signaling pathways presented in Table 3 and in Fig. 7. DE genes products are presented at the ECM and cell membrane localization. This result is consistent with the importance of the ECM highlighted in muscle cattle in development . As reviewed by , the extracellular matrix (ECM) is established to be a key player in muscle growth. ECM was most often described as an inactive component of cells. But many studies across species have now highlighted ECM functions, e.g. filtering, activating/inhibiting enzymatic activities, binding hormones, enzymes, and regulating the interaction of several ligands with their receptors. Most of the constituents of the ECM are mainly produced locally by the adjacent cells and most often by the fibroblasts present in the muscle tissue.
Also, PI3 kinase pathways regulate some energy metabolic pathways such as glycolysis and glycogenesis. Then selection on muscle mass may have also affected the regulation of these metabolic pathways. For example it has been suggested that decades of selection for more lean meat (more muscle mass) and larger litter size, may have increased piglet neonatal mortality . As piglets have no brown adipose tissue, piglet thermoregulation at birth is essentially carried out by the skeletal muscle (shivering; ). It may be hypothesized that impaired glycolytic metabolism at birth, which is essential to ensure body thermoregulation, has been affected by genetic selection, affecting gene expression between breeds [37–39].
Identification of some breeds specificities
15 genes that discriminate duroc, Pietrain and large white
A total of 171 unique annotated genes were identified as able to differentiate the breeds. Among these 171 genes (from 228 probes), 15 genes were common in the two discriminant methods (Random Forest and sPLS-DA) and the pairwise differential analysis (Fig. 6). These 15 genes may be very interesting to better characterize the differences between the breeds. None of these genes were located on the sexual chromosome, which may imply they are less affected by the sex effect (confounding with the Piétrain breed) observed in this work.
From these 15 genes, six genes differentiate the Duroc from the other breeds. Among the four genes under-expressed in Duroc, two genes are specifically known to be involved in skin coloration: KIT (v-kit Hardy-Zuckerman 4 feline sarcoma viral oncogene homolog or Dominant white locus in pigs; ) and OCA2 (oculocutaneous albinism II). Both genes are coding transporter proteins and are involved in biological processes as transmembrane receptor protein tyrosine kinase signaling pathway for KIT and tyrosine transport for OCA2. OCA2 gene was found over-expressed in white muscle (longissimus dorsi) compared to red muscle (soleus) in Meishan pig . It is interesting to observe how genes related to skin colour and maybe submitted to selection may affect another tissue important for production traits. The term “melanocyte differentiation” was the first enriched gene ontology (q-value = 0.0043) with KIT and OCA2 genes; Duroc pigs are brown-red.
Two genes are over-expressed in Duroc, TJP2 (tight junction protein 2 or ZO2, zonula occludens-2; Fig. 3) and PSMB4 (proteasome (prosome, macropain) subunit, beta type, 4). TJP2 is a membrane-associated guanylate kinase and encoded protein functions as a component of the tight junction barrier for cellular permeability involved in intercellular communication . Moreover there is some evidence that at least another zonula occludens protein (ZO1) is involved in vascular remodelling processes . These biological functions, cellular permeability and vascular remodelling, are interesting functions related to muscle meat traits, especially water holding capacity . Similarly, the proteasomic proteins (as PSMB4) has a major role in degrading proteins in muscle cells  and the proteasome might be one of the endogenous proteolytic system contributing to meat texture development . PSMB4 with VDR, KIT and OCA2 genes are involved in “reproduction process” (q-value = 0.0098). For example, the VDR (Vitamin D receptor) gene is known to be involved in male reproduction . In Humans, it was observed to be down-regulated with leaner patient  and allelic variations in the VDR gene were identified to be associated with lean body mass and height in Human . In French pig production the Duroc is used to obtain terminal boars.
In this study, five genes were identified to differentiate the Piétrain line from the others. Four genes were under-expressed in Piétrain compared to the others (NAA20, MOB3C, RNGTT, DPH5). The NAA20 gene (N(alpha)-acetyltransferase 20, NatB catalytic subunit or NAT5 gene) encodes a protein involved in normal cell proliferation and it functions in posttranslational protein N-terminal acetylation process . Some specific proteins targeted by this N-terminal acetyltransferase were identified e.g. actin and tropomyosin; the human NatB complex depletion perturbs actin cytoskeleton and focal adhesion organization . The MOB3C gene (MOB kinase activator 3C), member of the MOB protein family, has been shown to regulate mitosis, cell proliferation, apoptosis, centrosome biology and morphological changes . The RNGTT gene (RNA guanylyltransferase and 5′-phosphatase) encodes a RNA guanylyltransferase involved in the regulation of gene expression with capping the 5′end of the mRNA (from NCBI/BioSystems). RGNTT was identified to be a potential candidate gene for average daily feed intake but not in Piétrain pigs . The DPH5 (diphthamide biosynthesis 5) gene is a component of the diphthamide synthesis pathway (from NCBI/Gene). This pathway regulates post-translational modifications . None of these four genes are directly linked with the characteristic phenotypes of Piétrain (e.g. conformation). The fifth and up-regulated gene in Piétrain is PIK3CG (phosphatidylinositol-4,5-bisphosphate 3-kinase, catalytic subunit gamma; Fig. 3). PIK3CG is the only gene with a possible link with conformation as is one of the genes involved in the regulation of muscle mass (see first paragraph of discussion).
Large white sire and dam lines
The male and female Large White lines were difficult to differentiate from each other. Six probes corresponding to five unique annotated genes were found significantly different with a relatively low FDR (5 %). Among these five genes, four (CPNE2, DHX33, NLRC5, RAD1) were overexpressed in the male line and are coding for proteins with nuclear localizations. One of these genes was NLRC5 (NLR family, CARD domain containing 5) and is the largest member of the NLR protein family (a NOD-like receptor). NLRC5 is an intracellular receptor involved in innate immune sensing; it is induced by interferons in case of pathogen infection [54, 55]. It is also involved in the regulation of kinase activity and NF-kappaB transcription factor activity (Biological Process Ontology from GeneCodis). The three other genes have mechanistic nuclear roles. RAD1 codes a cell cycle checkpoint protein required for cell cycle arrest and DNA damage repair. DHX33 (DEAH (Asp-Glu-Ala-His) box polypeptide 33) is identified to code a protein with an important role in rRNA transcription and cell proliferation. CPNE2 (copine II) codes a calcium-dependent membrane-binding protein that may regulate molecular events at the interface of the cell membrane and cytoplasm. Only one gene, PRKRIP1 (PRKR interacting protein 1 (IL11 inducible)) gene was over-expressed in Large White dam line. Little functional information is available for this gene except a negative regulation of protein kinase activity (Biological Process Ontology from GenoCodis) but it may play a role in cytokine-mediated biological functions .
IGF2 was among the few genes able to discriminate the Landrace from the others with random forest and the pairwise analysis (see boxplot on Fig. 3). One surprising point is that the sPLS-DA didn’t identify IGF2 as a discriminant gene. The over-expression of IGF2 in the other breeds is probably the consequence of the mutation in the IGF2 intron3 g.3072G > A  which leads to an increased muscle mass and a reduced backfat deposition. This mutation has detrimental effects in prolificacy described as a result of an excess of leanness diminishing reproductive performance of the sow [58, 59]. In French breeds, the mutation favorable (allelic frequency of A > 96 %) to increase muscle mass seems fixed in Large White, Piétrain and Duroc, while the allelic frequency is always segregating in Landrace (about 70 % of allele A and 30 % of B; ).
The second gene identified as Landrace specific was LBR (lamin B receptor). LBR encodes an integral inner nuclear membrane protein. This protein is involved in the sequestration of heterochromatin near the periphery and the nucleoli in mammalian nuclei  and also in sterol metabolism .
Transcriptomic diversity among breeds
Genetic diversity of Western pig breeds and lines has been studied over the past years using a range of genetic markers [63–67]. A clear separation of breeds was observed at markers that were specifically chosen to vary across breeds. On the transcriptome level, we observed here that breeds were also clearly separated on differential probes. The differential aspect of probes is equivalent to the pre-selection of genetic markers. In this study, after eliminating the confounding effect of sex and Piétrain, Duroc animals are the most distant to other breeds; that corresponds to the genetic basis and historical knowledge (Duroc came from America while Piétrain, Landrace and Large White originated from Europe). Then Piétrain and Large White are the most distant, again as with the genetic markers. So the overall picture of breed differentiation is most probably the same at the genomic and transcriptomic level.
Specific statistical tools are now able to detect signatures of selection at the transcriptomic level and compare them to the ones at the genomic level. This will be the focus of future work.
However, transcriptomic differences across breeds may reveal the impact of selection and help understanding genetic and phenotypic changes. Perez-Enciso et al.  also reported interesting patterns of differential gene expression across breeds in a study involving 16 animals per breed (Large White, Duroc, Iberian, and a cross with a Sino-European hybrid line). They observed hierarchical clusterings of breeds differing with tissues. Even if the five tissues were chosen from the endocrine axis (hypothalamus, adenohypophysis, thyroid gland, gonads and fat tissue), it is interesting to note that muscle differentiation was highlighted in their breeds’ transcriptomic differentiation.
More studies are needed to understand the effects of the various evolutionary forces (amongst which is selection) on the transcriptome of pig breeds.
We found numerous differences in gene expression across the main pig Western breeds. Their functional analysis highlighted which biological function differed between these breeds and more precisely which genes from these functions differed between breeds. We hypothesized that these genes and related biological functions may have been targeted by human management (among which artificial selection) on production traits. The ideal situation would be that only terminal mechanisms (muscle mass, lipid metabolism…) for production traits vary across breeds. However our results showed that a multi-functional signaling pathway (PI3K) was affected. This pathway is well-known to be involved in the regulation of muscle mass and may impact the production traits that are selected for, but also “functional” traits around the regulation of metabolic pathway such as glycolysis. Moreover, most of the top discriminant genes between breeds were found to affect fundamental transcriptional and post-transcriptional mechanisms. These results could suggest how animal management and/or genetic selection may impact not only the terminal mechanisms of production traits but also the fundamental regulation of cellular processes such as regulation of gene expression and signaling pathways. Further analysis will help to decipher if robustness in farm animals was also impacted.
All procedures and facilities were approved by French veterinary services (ethics committee: Direction Départementale de la Cohésion Sociale et de la Protection des Populations in Rennes, France; agreement number A35-240-7). All animals were raised at the French central test Station in Le Rheu (France) in 2007 and 2008, and slaughtered in the same commercial slaughterhouse (Cooperl-Hunaudaye, Montfort-sur-Meu, France). A total of 150 individuals were sampled, and came from 5 breeds: 24 Duroc (DU), 33 Landrace (LR), 41 Large White dam line (LWF), 10 Large White sire line (LWM) and 39 Piétrain (PI). These animals were castrates in all breeds except females in PI, were reared in 10 different contempory groups, slaughtered in 29 different series, each containing several breeds. Further details can be found in . Muscle samples were biopsied from Longissimus dorsi (LD) muscle 20’ after stunning and exsanguination. Samples were immediately frozen in liquid N2 and kept at −80 °C until analysis.
Total RNA extraction
The total RNA extraction was previously described in . Total RNA was isolated from each of the 150 muscle samples. Briefly, the muscle samples were disrupted, homogenized and ground to a fine powder by rapid agitation for 1 min in a liquid-nitrogen cooled grinder with stainless steel beads. An aliquot of 250–300 mg of the fine powder was then processed for total RNA isolation and purification using RNeasy Fibrous Tissue Midi kit according to the manufacturer’s instructions (Qiagen SA France, Courtaboeuf, France). The method included a proteinase K digestion step to remove proteins and a DNase digestion step to remove contaminating DNA. The extracted total RNA was eluted in 300 μl of RNase-free water and stored at −80 °C. RNA quality and concentration were controlled using an AGILENT 2100 bioanalyzer (RNA solutions and RNA 6000 Nano Lab- Chip Kit, Agilent Technologies France, Massy, France).
Microarray data: hybridization
The 4 x 44 K Porcine Gene Expression Microarray (G2514F, V1: 020109 Agilent Technologies; GEO accession number GPL10162) used in this work was previously described by the manufacturer. The 43,803 porcine probes sourced from UniGene (Release 33, Feb 2008), RefSeq (Release 27, Jan 2008), TIGR (Release 12, Jun 2006) using Agilent 60-mer SurePrint technology.
RNA concentrations were determined using a NanoDrop® ND-1000 spectrophotometer (NanoDrop Technologies, Wilmington, DE). cRNA labeled with fluorescent Cyanine 3-CTP was used for hybridization at 65 °C for 17 h onto porcine oligo microarray slides, according to the manufacturer’s recommendations (Low Input Quick Amp Labeling with Low Input QuickAmp Cy3 labeling kit One-color and One-Color Microarray-Based Gene Expression Analysis; Agilent Technologies). Hybridized microarray slides were scanned with scanner Genepix 4000B (AXON INSTRUMENTS) at 5-μm resolution. The scanned images were analyzed numerically using Agilent Feature Extraction Software version 188.8.131.52. (Agilent Technologies).
Microarray data: normalization
The whole set of 45,220 spots were analysed on 181 chips (one sample per chip, a limited number of animals giving 2 to 3 samples). One chip was removed because of overall bad quality. Each spot on each chip was allocated a weight of 0 or 1 depending on various criteria of spot’s quality and signal to background level. On average, half of the spots per chip had a 0 weight. Only spots with a weight of 1 for at least 70 % of the chips were kept for further analysis, leading to 12,358 spots. All intensity data were log transformed, and the term “data” or “intensity” referred to log (natural logarithm) transformed data throughout the manuscript. The effect of the hybridization day was removed by subtracting the mean intensity of the hybridization day. Then the mean of each chip was subtracted to make all chips comparable.
A good linearity of intensity signals were observed between all samples. One animal was hybridized 3 times, and 32 individuals were sampled and hybridized twice. A very good adequacy was observed between samples from the same individual, with an average correlation of 0.98. Then, for each individual, only the sample with the best mean quality was kept. Finally, 147 chips representing 147 individuals were available for deeper statistical analysis. The normalized microarray dataset have been deposited in the NCBI Gene Expression Omnibus GSE56011.
A linear mixed model was fitted for each spot, with breed as fixed effect and contemporary groups as random effect, with the lme function of the R package nlme. The significance of the breed effect on gene expression was evaluated with an F-type test (R function anova). A False Discovery Rate (FDR, ) correction was performed on raw p-values (R function multtest). In the same linear mixed model, comparisons of each pair of breed were tested, with a FDR correction thereafter. In all cases, a threshold of 0.1 % has been applied except for the comparison between LWM and LWF (5 % threshold). This stringent threshold was chosen to ease the functional analysis; no enrichment can be highlighted with too many genes.
The fold change (FC) for any pair of breeds was calculated as the exponential of the difference between mean expression level in breed 1 and mean expression level in breed 2. This FC corresponds to the intensity ratio of mean raw signals on the microarray between breeds.
Hierarchical clusterings were built with the Euclidian distance and the Ward aggregation criterion for each clustering.
A Random Forest (RF) analysis was performed in order to predict the breed of each animal on the basis of its transcriptomic profile, using the random Forest package of the R software . A preliminary RF analysis was used to eliminate the 70 % less important genes to avoid too much noise. On the basis of the remaining 30 %, i.e. 3,336 spots, a second RF analysis was performed with 5,000 trees (the other parameters were the default ones).
The Agilent probes were annotated searching sequence homologies against following databases: SwissProt, TIGR Pig SsGI 12, UniGene Pig, Ensembl Human Transcripts NCBI36 (annotation from SIGENAE, http://www.sigenae.org/). The annotation of the DEG is summarized in Additional file 2. This file included the localization of the probes on Sus scrofa genome (Sscrofa10.2.69) and the human, bovine and mice orthologs were added when available (ortholog_one2one). Finally, some differentially expressed genes were manually annotated with blastn against the Refseq_RNA library (NCBI) or with blat against the pig genome (Ensembl, Sscrofa10.2 version).
From the Database for Annotation, Visualization and Integrated Discovery (DAVID), the software EASE (an Expression Analysis Systematic Explorer; https://david.ncifcrf.gov/ease/ease.jsp) was used to obtain functional Gene Ontology (GO) terms for each gene, the associated KEGG pathway, and the function summary. The systematic ontological annotation is given in Additional file 6.
In a second step, the GeneCodis website was used to identify co-occurrence in the functional annotation to highlight functions specifically enriched (http://genecodis.cnb.csic.es/; ). The Human genome was used as reference. We focused on KEGG pathways to avoid some redundant functional information. GeneCodis calculates an adjusted p-value. Most of the results for this work were obtained with GeneCodis with top pathways presented in Table 3 and details in Additional file 4.
Availability of supporting data
The data set supporting the results of this article is available in the Gene Expression Omnibus (GEO) repository, with an accession number GSE56011, at http://www.ncbi.nlm.nih.gov/geo/query/acc.cgi?acc=GSE56011.
database for annotation, visualization and integrated discovery
differentially expressed genes
DEAH (Asp-Glu-Ala-His) box polypeptide 33
GenBank accession number for a non-annotated gene
diphthamide biosynthesis 5
expression analysis systematic explorer
false discovery rate
insulin-like growth factor 2
ingenuity pathways analysis
kyoto encyclopedia of genes and genomes
V-kit hardy-zuckerman 4 feline sarcoma viral oncogene homolog
large white: large white dam line
large white sire line
MOB kinase activator 3C
N (Alpha)-acetyltransferase 20, NatB catalytic subunit
NOD-like receptor family CARD domain containing 5
oculocutaneous albinism II
principal component analysis
phosphatidylinositol 3-kinase, catalytic subunit type 3
phosphatidylinositol-4,5-bisphosphate 3-kinase, catalytic subunit gamma
partial least squares discriminant analysis
PRKR interacting protein 1 (IL11 Inducible)
proteasome (prosome, macropain) subunit, beta type, 4
RAB18, member RAS oncogene family
checkpoint DNA exonuclease
RNA guanylyltransferase and 5′-Phosphatase
sus scrofa chromosome
tight junction protein 2 or ZO2, zonula occludens-2
- USP9X or USP9Y:
ubiquitin specific peptidase 9, x-linked or y-linked
vitamin D receptor
This work is part of the DéLiSus project funded by BIOPORC (ADN, Choice Genetics France, Gene + and Nucléus breeding organizations and IFIP) and French ANR (ANR-07-GANI-001). We thank the animal and DNA providers (BIOPORC). F.R. acknowledges financial support from Région Midi-Pyrénées and Département de Génétique Animale of INRA. The microarray hybridizations and image analysis have been realized at the Biochip Platform http://biopuce.insa-toulouse.fr/. The microarray annotation was performed by SIGENAE computer group (Système d’information du projet d’analyse des génomes des animaux d’élevage, http://www.sigenae.org). We are grateful to Zhilei Li, Marine Dumont, Agnès Solomiac, Laurence Souchu (students of the Institut National des Sciences Appliquées de Toulouse, France) for an earlier contribution in the statistical analysis, and David Robelin, Bertrand Servin and Nathalie Villa-Vialaneix for their help. Thanks to the referees for their valuable comments.
Open AccessThis article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.
- Romero IG, Ruvinsky I, Gilad Y. Comparative studies of gene expression and the evolution of gene regulation. Nat Rev Genet. 2012;13.
- Whitehead A, Crawford D. Variation within and among species in gene expression: raw material for evolution. Mol Biol Evol. 2006;15(5):1197–211.Google Scholar
- Oleksiak M, Churchill G, Crawford D. Variation in gene expression within and among natural populations. Nat Genet. 2002;32(2):261–6.View ArticlePubMedGoogle Scholar
- Rifkin S, Kim J, White K. Evolution of gene expression in the Drosophila melanogaster subgroup. Nat Genet. 2003;33(2):138–44.View ArticlePubMedGoogle Scholar
- Jimenez-Guri E, Huerta-Cepas J, Cozzuto L, Wotton KR, Kang H, Himmelbauer H, et al. Comparative transcriptomics of early dipteran development. BMC Genomics. 2013;14:123.PubMed CentralView ArticlePubMedGoogle Scholar
- Enard W, Khaitovich P, Klose J, Zöllner S, Heissig F, Giavalisco P, et al. Intra- and interspecific variation in primate gene expression patterns. Science. 2002;296(5566):340–3.View ArticlePubMedGoogle Scholar
- Khaitovich P, Weiss G, Lachmann M, Hellmann I, Enard W, Muetzel B, et al. A neutral model of transcriptome evolution. PLoS Biol. 2004;2(5), E132.PubMed CentralView ArticlePubMedGoogle Scholar
- Busby MA, Gray JM, Costa AM, Stewart C, Stromberg MP, Barnett D, et al. Expression divergence measured by transcriptome sequencing of four yeast species. BMC Genomics. 2011;12:635.PubMed CentralView ArticlePubMedGoogle Scholar
- Le Gall T, Darlu P, Escobar-Páramo P, Picard B, Denamur E. Selection-driven transcriptome polymorphism in Escherichia coli/Shigella species. Genome Res. 2005;15(2):260–8.PubMed CentralView ArticlePubMedGoogle Scholar
- Giger T, Excoffier L, Amstutz U, Day P, Champigneulle A, Hansen M, et al. Population transcriptomics of life-history variation in the genus Salmo. Mol Ecol. 2008;17(13):3095–108.View ArticlePubMedGoogle Scholar
- Kristiansson E, Osterlund T, Gunnarsson L, Arne G, Larsson DG, Nerman O. A novel method for cross-species gene expression analysis. BMC Bioinformatics. 2013;14(1):70.PubMed CentralView ArticlePubMedGoogle Scholar
- Yang L, Zou M, Fu B, He S. Genome-wide identification, characterization, and expression analysis of lineage-specific genes within zebrafish. BMC Genomics. 2013;14:65.PubMed CentralView ArticlePubMedGoogle Scholar
- Atanasova L, Crom SL, Gruber S, Coulpier F, Seidl-Seiboth V, Kubicek CP, et al. Comparative transcriptomics reveals different strategies of Trichoderma mycoparasitism. BMC Genomics. 2013;14:121.PubMed CentralView ArticlePubMedGoogle Scholar
- Hufford MB, Xu X, Van Heerwaarden J, Pyhajarvi T, Chia JM, Cartwright RA, et al. Comparative population genomics of maize domestication and improvement. Nat Genet. 2012;44(7):808–11.View ArticlePubMedGoogle Scholar
- Jiang SY, Ma Z, Vanitha J, Ramachandran S. Genetic variation and expression diversity between grain and sweet sorghum lines. BMC Genomics. 2013;14:18.PubMed CentralView ArticlePubMedGoogle Scholar
- Natt D, Rubin CJ, Wright D, Johnsson M, Belteky J, Andersson L, et al. Heritable genome-wide variation of gene expression and promoter methylation between wild and domesticated chickens. BMC Genomics. 2012;13:59.PubMed CentralView ArticlePubMedGoogle Scholar
- Yang HC, Wang PL, Lin CW, Chen CH. Integrative analysis of single nucleotide polymorphisms and gene expression efficiently distinguishes samples from closely related ethnic populations. BMC Genomics. 2012;13(1):346.PubMed CentralView ArticlePubMedGoogle Scholar
- Müller L, Hutter S, Stamboliyska R, Saminadin-Peter S, Stephan W, Parsch J. Population transcriptomics of Drosophila melanogaster females. BMC Genomics. 2011;12:81.PubMed CentralView ArticlePubMedGoogle Scholar
- Fraser H, Moses A, Schadt E. Evidence for widespread adaptive evolution of gene expression in budding yeast. Proc Natl Acad Sci U S A. 2010;107(7):2977–82.PubMed CentralView ArticlePubMedGoogle Scholar
- Mensack M, Fitzgerald V, Ryan E, Lewis M, Thompson H, Brick M. Evaluation of diversity among common beans (Phaseolus vulgaris L.) from two centers of domestication using ‘omics’ technologies. BMC Genomics. 2010;11:686–720.PubMed CentralView ArticlePubMedGoogle Scholar
- Pérez-Enciso M, Ferraz A, Ojeda A, López-Béjar M. Impact of breed and sex on porcine endocrine transcriptome: a bayesian biometrical analysis. BMC Genomics. 2009;10:89.PubMed CentralView ArticlePubMedGoogle Scholar
- Rohart F, Paris A, Laurent B, Canlet C, Molina J, Mercat MJ, et al. Phenotypic prediction based on metabolomic data for growing pigs from three main European breeds. J Anim Sci. 2012;90(13):4729–40.View ArticlePubMedGoogle Scholar
- Rubin CJ, Megens HJ, Martinez Barrio A, Maqbool K, Sayyab S, Schwochow D, et al. Strong signatures of selection in the domestic pig genome. Proc Natl Acad Sci U S A. 2012;109(48):19529–36.PubMed CentralView ArticlePubMedGoogle Scholar
- Andersson L, Plastow GS. Molecular genetics of coat colour variation. In: Rothschild MF, Ruvinsky A, editors. The Genetics of the Pig. Wallingford, UK: CAB Internationa; 2011. p. 38–50.View ArticleGoogle Scholar
- Bem D, Yoshimura S, Nunes-Bastos R, Bond FC, Kurian MA, Rahman F, et al. Loss-of-function mutations in RAB18 cause Warburg micro syndrome. Am J Hum Genet. 2011;88(4):499–507.PubMed CentralView ArticlePubMedGoogle Scholar
- Vicoso B, Charlesworth B. Evolution on the X chromosome: unusual patterns and processes. Nat Rev Genet. 2006;7(8):645–53.View ArticlePubMedGoogle Scholar
- Tabas-Madrid D, Nogales-Cadenas R, Pascual-Montano A. GeneCodis3: a non-redundant and modular enrichment analysis tool for functional genomics. Nucleic Acids Res. 2012;40:W478–83. Web Server issue.PubMed CentralView ArticlePubMedGoogle Scholar
- Breiman L. Random forests. Mach Learn. 2001;45(1):5–32.View ArticleGoogle Scholar
- Bonnet A, Le Cao KA, Sancristobal M, Benne F, Robert-Granie C, Law-So G, et al. In vivo gene expression in granulosa cells during pig terminal follicular development. Reproduction. 2008;136(2):211–24.View ArticlePubMedGoogle Scholar
- Lê Cao K-A, Boitard S, Besse P. Sparse PLS discriminant analysis: biologically relevant feature selection and graphical displays for multiclass problems. BMC Bioinformatics. 2011;12:253.PubMed CentralView ArticlePubMedGoogle Scholar
- Glass DJ. Molecular mechanisms modulating muscle mass. Trends Mol Med. 2003;9(8):344–50.View ArticlePubMedGoogle Scholar
- Egerman MA, Glass DJ. Signaling pathways controlling skeletal muscle mass. Crit Rev Biochem Mol Biol. 2014;49(1):59–68.PubMed CentralView ArticlePubMedGoogle Scholar
- Foxcroft GR, Dixon WT, Novak S, Putman CT, Town SC, Vinsky MD. The biological basis for prenatal programming of postnatal performance in pigs. J Anim Sci. 2006;84(Suppl):E105–12.PubMedGoogle Scholar
- Guo B, Greenwood PL, Cafe LM, Zhou G, Zhang W, Dalrymple BP. Transcriptome analysis of cattle muscle identifies potential markers for skeletal muscle growth rate and major cell types. BMC Genomics. 2015;16:177.PubMed CentralView ArticlePubMedGoogle Scholar
- Thorsteinsdottir S, Deries M, Cachaco AS, Bajanca F. The extracellular matrix dimension of skeletal muscle development. Dev Biol. 2011;354(2):191–207.View ArticlePubMedGoogle Scholar
- Canario L, Pere MC, Tribout T, Thomas F, David C, Gogue J, et al. Estimation of genetic trends from 1977 to 1998 of body composition and physiological state of Large White pigs at birth. Int J Anim Biosci. 2007;1(10):1409–13.Google Scholar
- Herpin P, Damon M, Le Dividich J. Development of thermoregulation and neonatal survival in pigs. Livest Prod Sci. 2002;78(1):25–45.View ArticleGoogle Scholar
- van der Lende T, Knol EF, Leenhouwers JI. Prenatal development as a predisposing factor for perinatal losses in pigs. Reprod Suppl. 2001;58:247–61.PubMedGoogle Scholar
- Voillet V, SanCristobal M, Lippi Y, Martin PG, Iannuccelli N, Lascor C, et al. Muscle transcriptomic investigation of late fetal development identifies candidate genes for piglet maturity. BMC Genomics. 2014;15(1):797.PubMed CentralView ArticlePubMedGoogle Scholar
- Fontanesi L, D’Alessandro E, Scotti E, Liotta L, Crovetti A, Chiofalo V, et al. Genetic heterogeneity and selection signature at the KIT gene in pigs showing different coat colours and patterns. Anim Genet. 2010;41(5):478–92.View ArticlePubMedGoogle Scholar
- Li Y, Xu Z, Li H, Xiong Y, Zuo B. Differential transcriptional analysis between red and white skeletal muscle of Chinese Meishan pigs. Int J Biol Sci. 2010;6(4):350–60.PubMed CentralView ArticlePubMedGoogle Scholar
- Tkachuk N, Tkachuk S, Patecki M, Kusch A, Korenbaum E, Haller H, et al. The tight junction protein ZO-2 and Janus kinase 1 mediate intercellular communications in vascular smooth muscle cells. Biochem Biophys Res Commun. 2011;410(3):531–6.View ArticlePubMedGoogle Scholar
- Laing JG, Saffitz JE, Steinberg TH, Yamada KA. Diminished zonula occludens-1 expression in the failing human heart. Cardiovasc Pathol. 2007;16(3):159–64.View ArticlePubMedGoogle Scholar
- Huff-Lonergan E, Lonergan SM. Mechanisms of water-holding capacity of meat: The role of postmortem biochemical and structural changes. Meat Sci. 2005;71(1):194–204.View ArticlePubMedGoogle Scholar
- Lee SH, Joo ST, Ryu YC. Skeletal muscle fiber type and myofibrillar proteins in relation to meat quality. Meat Sci. 2010;86(1):166–70.View ArticlePubMedGoogle Scholar
- Ouali A, Gagaoua M, Boudida Y, Becila S, Boudjellal A, Herrera-Mendez CH, et al. Biomarkers of meat tenderness: present knowledge and perspectives in regards to our current understanding of the mechanisms involved. Meat Sci. 2013;95(4):854–70.View ArticlePubMedGoogle Scholar
- Blomberg Jensen M. Vitamin D and male reproduction. Nat Rev Endocrinol. 2014;10(3):175–86.View ArticlePubMedGoogle Scholar
- Patel HP, Al-Shanti N, Davies LC, Barton SJ, Grounds MD, Tellam RL, et al. Lean mass, muscle strength and gene expression in community dwelling older men: findings from the Hertfordshire Sarcopenia Study (HSS). Calcif Tissue Int. 2014;95(4):308–16.View ArticlePubMedGoogle Scholar
- Ferrarezi DA, Bellili-Munoz N, Nicolau C, Cheurfa N, Guazzelli IC, Frazzatto E, et al. Allelic variations in the vitamin D receptor gene, insulin secretion and parents’ heights are independently associated with height in obese children and adolescents. Metab Clin Exp. 2012;61(10):1413–21.View ArticlePubMedGoogle Scholar
- Starheim KK, Arnesen T, Gromyko D, Ryningen A, Varhaug JE, Lillehaug JR. Identification of the human N (alpha)-acetyltransferase complex B (hNatB): a complex important for cell-cycle progression. Biochem J. 2008;415(2):325–31.View ArticlePubMedGoogle Scholar
- Van Damme P, Lasa M, Polevoda B, Gazquez C, Elosegui-Artola A, Kim DS, et al. N-terminal acetylome analyses and functional insights of the N-terminal acetyltransferase NatB. Proc Natl Acad Sci U S A. 2012;109(31):12449–54.PubMed CentralView ArticlePubMedGoogle Scholar
- Chow A, Hao Y, Yang X. Molecular characterization of human homologs of yeast MOB1. Int J Cancer J Int Du Cancer. 2010;126(9):2079–89.Google Scholar
- Su X, Lin Z, Lin H. The biosynthesis and biological function of diphthamide. Crit Rev Biochem Mol Biol. 2013;48(6):515–21.PubMed CentralView ArticlePubMedGoogle Scholar
- Yao Y, Qian Y. Expression regulation and function of NLRC5. Protein Cell. 2013;4(3):168–75.View ArticlePubMedGoogle Scholar
- Kobayashi KS, van den Elsen PJ. NLRC5: a key regulator of MHC class I-dependent immune responses. Nat Rev Immunol. 2012;12(12):813–20.View ArticlePubMedGoogle Scholar
- Yin Z, Haynie J, Williams BR, Yang YC. C114 is a novel IL-11-inducible nuclear double-stranded RNA-binding protein that inhibits protein kinase R. J Biol Chem. 2003;278(25):22838–45.View ArticlePubMedGoogle Scholar
- Van Laere AS, Nguyen M, Braunschweig M, Nezer C, Collette C, Moreau L, et al. A regulatory mutation in IGF2 causes a major QTL effect on muscle growth in the pig. Nature. 2003;425(6960):832–6.View ArticlePubMedGoogle Scholar
- Buys N, van den Abeele G, Stinckens A, Deley J, Georges M. Effect of the IGF2-intron3-G3072A mutation on prolificacy in sows. In: 8th World Congress on Genetics Applied to Livestock Production: 13–18 August, 2006. Minas Gerais, Brazil: Belo Horizonte; 2006. p. 6–22.Google Scholar
- Stinckens A, Mathur P, Janssens S, Bruggeman V, Onagbesan OM, Schroyen M, et al. Indirect effect of IGF2 intron3 g.3072G > A mutation on prolificacy in sows. Anim Genet. 2010;41(5):493–8.View ArticlePubMedGoogle Scholar
- Mercat MJ, Feve K, Mûller N, Swob S, Le Roy P, Bidanel JP, et al. Estimation, dans un dispositif familial issu des populations porcines françaises en sélection, de l’effet quantitatif de mutations dans des gènes majeurs et des gènes candidats. Journées Recherche Porcine. 2012;44:1–6.Google Scholar
- Politz JC, Ragoczy T, Groudine M. When untethered, something silent inside comes. Nucleus. 2013;4(3):153–5.PubMed CentralView ArticlePubMedGoogle Scholar
- Clayton P, Fischer B, Mann A, Mansour S, Rossier E, Veen M, et al. Mutations causing Greenberg dysplasia but not Pelger anomaly uncouple enzymatic from structural functions of a nuclear membrane protein. Nucleus. 2010;1(4):354–66.PubMed CentralView ArticlePubMedGoogle Scholar
- Laval G, Iannuccelli N, Legault C, Milan D, Groenen MA, Giuffra E, et al. Genetic diversity of eleven European pig breeds. Genet Sel Evol. 2000;32(2):187–203.PubMed CentralView ArticlePubMedGoogle Scholar
- SanCristobal M, Chevalet C, Haley CS, Joosten R, Rattink AP, Harlizius B, et al. Genetic diversity within and between European pig breeds using microsatellite markers. Anim Genet. 2006;37(3):189–98.View ArticlePubMedGoogle Scholar
- SanCristobal M, Chevalet C, Peleman J, Heuven H, Brugmans B, Van Schriek M, et al. Genetic diversity in European pigs utilizing amplified fragment length polymorphism markers. Anim Genet. 2006;37(3):232–8.View ArticlePubMedGoogle Scholar
- Megens HJ, Crooijmans RPMA, Cristobal MS, Hui X, Li N, Groenen MAM. Biodiversity of pig breeds from China and Europe estimated from pooled DNA samples: differences in microsatellite variation between two areas of domestication. Genet Sel Evol. 2008;40(1):103–28.PubMed CentralPubMedGoogle Scholar
- Amaral AJ, Ferretti L, Megens HJ, Crooijmans RP, Nie H, Ramos-Onsins SE, et al. Genome-wide footprints of pig domestication and selection revealed through massive parallel sequencing of pooled DNA. PLoS One. 2011;6(4):e14782.PubMed CentralView ArticlePubMedGoogle Scholar
- Ferre PJ, Liaubet L, Concordet D, SanCristobal M, Uro-Coste E, Tosser-Klopp G, et al. Longitudinal analysis of gene expression in porcine skeletal muscle after post-injection local injury. Pharm Res. 2007;24(8):1480–9.View ArticlePubMedGoogle Scholar
- Benjamini Y, Hochberg Y. Controlling the False Discovery Rate - a Practical and Powerful Approach to Multiple Testing. J Roy Stat Soc B Met. 1995;57(1):289–300.Google Scholar
- Liaw A, Wiener M. Classification and Regression by randomForest. R News. 2002;2(3):18–22.Google Scholar