|Identification of expression QTL (eQTL) of genes expressed in porcine M. longissimus dorsi and associated with meat quality traits.|
|Jump to Full Text|
|PMID: 20950486 Owner: NLM Status: MEDLINE|
|BACKGROUND: Genetic analysis of transcriptional profiles is a promising approach for identifying and dissecting the genetics of complex traits like meat performance. Accordingly, expression levels obtained by microarray analysis were taken as phenotypes in a linkage analysis to map eQTL. Moreover, expression levels were correlated with traits related to meat quality and principle components with high loadings of these traits. By using an up-to-date annotation and localization of the respective probe-sets, the integration of eQTL mapping data and information of trait correlated expression finally served to point to candidate genes for meat quality traits.
RESULTS: Genome-wide transcriptional profiles of M. longissimus dorsi RNAs samples of 74 F2 animals of a pig resource population revealed 11,457 probe-sets representing genes expressed in the muscle. Linkage analysis of expression levels of these probe-sets provided 9,180 eQTL at the suggestive significance threshold of LOD > 2. We mapped 653 eQTL on the same chromosome as the corresponding gene and these were designated as 'putative cis-eQTL'. In order to link eQTL to the traits of interest, probe-sets were addressed with relative transcript abundances that showed correlation with meat quality traits at p ≤ 0.05. Out of the 653 'putative cis-eQTL', 262 transcripts were correlated with at least one meat quality trait. Furthermore, association of expression levels with composite traits with high loadings for meat quality traits generated by principle component analysis were taken into account leading to a list of 85 genes exhibiting cis-eQTL and trait dependent expression.
CONCLUSION: Holistic expression profiling was integrated with QTL analysis for meat quality traits. Correlations between transcript abundance and meat quality traits, combined with genetic positional information of eQTL allowed us to prioritise candidate genes for further study.
|Siriluck Ponsuksili; Eduard Murani; Manfred Schwerin; Karl Schellander; Klaus Wimmers|
Related Documents :
|20069416 - A new semi-dwarfing gene identified by molecular mapping of quantitative trait loci in ...
21724876 - Analysis and manipulation of aspartate pathway genes for l-lysine overproduction from m...
19331686 - Robust physical methods that enrich genomic regions identical by descent for linkage st...
18563486 - The genomic determinants of alcohol preference in mice.
1827756 - Human d-type cyclin.
8318696 - Local and systemic signals in the wound response.
|Type: Journal Article; Research Support, Non-U.S. Gov't Date: 2010-10-16|
|Title: BMC genomics Volume: 11 ISSN: 1471-2164 ISO Abbreviation: BMC Genomics Publication Date: 2010|
|Created Date: 2010-11-03 Completed Date: 2010-12-27 Revised Date: 2013-07-03|
Medline Journal Info:
|Nlm Unique ID: 100965258 Medline TA: BMC Genomics Country: England|
|Languages: eng Pagination: 572 Citation Subset: IM|
|Leibniz Institute for Farm Animal Biology, Research Unit Molecular Biology, Wilhelm-Stahl-Allee 2, 18196 Dummerstorf, Germany.|
|APA/MLA Format Download EndNote Download BibTex|
Chromosomes, Mammalian / genetics
DNA Probes / metabolism
Gene Expression Regulation*
Meat / standards*
Muscle, Skeletal / metabolism*
Principal Component Analysis
Quantitative Trait Loci / genetics*
Quantitative Trait, Heritable*
RNA, Messenger / genetics, metabolism
Sus scrofa / genetics*
|0/DNA Probes; 0/RNA, Messenger|
Journal ID (nlm-ta): BMC Genomics
Publisher: BioMed Central
Copyright ©2010 Ponsuksili et al; licensee BioMed Central Ltd.
Received Day: 23 Month: 5 Year: 2010
Accepted Day: 16 Month: 10 Year: 2010
collection publication date: Year: 2010
Electronic publication date: Day: 16 Month: 10 Year: 2010
Volume: 11First Page: 572 Last Page: 572
Publisher Id: 1471-2164-11-572
PubMed Id: 20950486
|Identification of expression QTL (eQTL) of genes expressed in porcine M. longissimus dorsi and associated with meat quality traits|
|Siriluck Ponsuksili1||Email: email@example.com|
|Eduard Murani2||Email: firstname.lastname@example.org|
|Manfred Schwerin1||Email: email@example.com|
|Karl Schellander3||Email: firstname.lastname@example.org|
|Klaus Wimmers2||Email: email@example.com|
1Leibniz Institute for Farm Animal Biology, Research Group Functional Genome Analysis, Wilhelm-Stahl-Allee 2, 18196 Dummerstorf, Germany
2Leibniz Institute for Farm Animal Biology, Research Unit Molecular Biology, Wilhelm-Stahl-Allee 2, 18196 Dummerstorf, Germany
3Institute of Animal Science, Animal Breeding and Husbandry Group, University of Bonn, Endenicher Allee 15, 53115 Bonn, Germany
Genetical genomics as a new approach which combines gene-expression data and marker genotypes in a segregating population, offers great perspectives to make a major contribution to the dissection of complex traits [1,2]. Genetical genomics aims at detecting genomic loci that control variation in gene expression, so-called expression QTL (eQTL; to distinguish them from functional QTL that affect traits at the whole-organism level, subsequently termed pheneQTL(pQTL)). The detected eQTL can represent a locus that lies close to the gene that is being controlled (cis-acting) or one or more loci that are unlinked to the gene that is being controlled (trans-acting) . Expression-QTL for genes showing high correlation with the phenotype may provide the necessary information required for identifying genes that control quantitative phenotypes. Those cis-eQTLs resulting from the correlation of expression profiles with phenotypic measurements represent candidate genes for the genetic regulation underlying the variation of the physiological traits [3,4].
We have previously reported on the identification of eQTL of genes showing expression levels correlated with waterholding capacity of meat measured as drip loss . Here in this study, we expand and up-dated the analysis towards global eQTL mapping of transcripts showing variable abundance in porcine muscle and using annotation and localisation data of Affymetrix microarray probe-sets based on the current porcine sequence information . We focus on eQTL of genes whose expression at slaughter is significantly correlated to technological meat quality traits addressed in commercial pig breeding schemes. Taking into account own observations on the reliability of eQTL mapping we highlight eQTL located on the same chromosome as the corresponding genes themselves .
Expression data were obtained from M. longissimus dorsi samples of 74 F2 animals of a resource population using Affymetrix Porcine Genome Arrays containing in total 24,123 probe-sets of which 20,689 probe-sets were assigned to a known gene . MAS5 analysis revealed consistent 'present calls' for 11,457 probe-sets. This pre-selected set was further analyzed with the more sophisticated hybrid algorithm PLIER [8,9] and the expression levels were subjected to linkage analysis. Out of 11,457 probe-sets 6,117 showed at least one eQTL at the 5% chromosome-wide significance threshold (average F-value 4.92, corresponding to LOD score = 2.0 and nominal significance of p ≤ 10-3) corresponding approximately to the suggestive linkage threshold proposed by Lander & Kruglyak  (Table 1). In total, the 11,457 probe-sets revealed 9,180 eQTL at LOD score > 2 (Additional file 1), 1,058 eQTL at LOD score > 3, and 160 eQTL at LOD score > 4. Only 29 had eQTL with LOD score > 5. This yielded an average of 1.5 eQTL per transcript (ranging from 1 to 6 eQTLs). The eQTL were distributed over all autosomes with some regions harboring particularly many eQTL with high significance (Figure 1). Out of the 9,180 eQTL 8,168 eQTL could be assigned to the porcine genome sequence (Ensembl Sscrofa9 database, released April 2009). In total 653 eQTL were mapped on the same chromosome as the corresponding gene itself. These eQTL are putative cis-eQTL.
In order to link the eQTL to the genetic background of a classical phenotypic trait of interest, it is necessary to establish a relationship between the variation of that classical phenotypic trait, the expression levels of various transcripts, and the mapping position of the eQTL and the corresponding transcripts. Therefore, we selected transcripts whose expression level showed significant correlation with individual meat quality traits and with the composite traits, PC2 and PC3 (Tables 2, 3, 4) (p ≤ 0.05, corresponding to correlation coefficients ranging between |0.24-0.50|) and we focussed on those transcripts with eQTL on their own chromosome as the structurally fixed linkage group (collateral transcript and eQTL position). Out of 653 eQTL that were mapped on the chromosome of their corresponding genes 262 transcripts showed correlation with at least one meat quality trait. These are displayed in Figure 2, 3, 4, 5 and 6. The linkage map showed accurate match of microsatellites marker order with the marker arrangement in the current porcine genome sequence assembly (Sscrofa9 assembly, April 2009). A total of 73 out 115 microsatellites marker were found to align on the porcine genome sequence assembly map position by BLAST analysis. Reflecting this alignment, the position of 262 transcripts in the porcine physical map was assigned to the corresponding eQTL position in the genetic linkage maps. Most of the eQTL of transcripts whose expression level correlated with meat quality traits were found on SSC1, 2, 3, 6, 9 and 14. Table 2 summarizes the number of eQTL related to various meat quality traits.
Most of these meat quality parameters were correlated or dependent on each other. Such correlations were reported by many studies [11-13]. In a previous study we used principle components (PCs) to reduce the multi dimensional data sets into lower dimensions . The principal component (PC2) is a meat quality vector and has high loadings for pH24ld, pH24st and meat color which are inversely correlated with drip loss. The principal component (PC3) is also a meat quality vector with large positive contribution of conductivity (LF1ld, LF24ld, LF24sm) and negative contribution of pH1ld. In order to scale down the list of 262 putative cis-eQTL for meat quality, transcript levels of genes correlated with the two PCs with high loadings of meat quality traits were considered. The level of expression of 38 and 47 probe-sets, respectively, were correlated with the new composite traits that were generated by principal component analysis PC2 and PC3 (Tables 3 and 4). In order to unravel the composite traits, the correlation of these 38 and 47 transcripts with each single meat quality trait contributing to the PC was screened and traits with significant correlations with the respective genes are given in tables 3 and 4. Thus the probe-sets represent a subset of probe-sets correlated with individual meat quality traits. As previously shown, analyses of enrichment of functional annotation groups as defined in the Gene Ontology (GO), the Kyoto Encyclopedia of Genes and Genomes (KEGG) databases, or in the Ingenuity Pathways Analysis library, highlight ubiquitination, phosphorylation, mitochondrion dysfunction, actin-, integrin-, PDGF-, EGF-, VEGF-, and Ca-signalling pathways . Most of the eQTL related to PC2 were found on chromosomes 2, 3, and 6; for PC3 on chromosomes 1, 3, 7, and 9.
Complex traits are genetically controlled by many loci. The identification of candidate genes for these traits is widely based on the principle of 'collecting evidences' in order to prioritize genes for further analysis from the huge lists of functional and positional candidate genes. The detection of eQTL provides information about the positional aspects of gene expression regulation. Together with the growing knowledge of genome sequences and gene annotation, this gives insight into the architecture of regulatory networks. In order to relate eQTL analysis to the genetic background of any complex trait, either a positional link to QTL for the complex trait of interest (pQTL) and/or a functional link to the trait expression is required. Here we report on the detection of eQTL for some 11,000 transcripts found in porcine skeletal muscle by microarray analysis. We demonstrate that the global microarray eQTL analysis can serve for narrowing down the candidate genes for quantitative traits related to meat quality when it is integrated with the analysis of the correlation of expression levels with the traits of interest.
The genetic architecture of transcript level variation in the porcine F2 resource population was highly variable and complex. eQTL were detected for 53% of the 11,457 probe-sets, with many of them exhibiting more than one eQTL summing up to more than 9000 eQTL in total (Table 1). The expression levels of the majority of the transcripts showed quantitative variation. So the transcript levels are probably quantitatively controlled. The link between the linkage maps and porcine assembly sequence enable discriminating cis or trans-eQTL. All 23,934 probe-sets sequence represented on the porcine microarray and the microsatellites sequences were BLASTed against the Ensembl porcine BAC sequence (Sscrofa9, April 2009). Our microsatellite order showed high accuracy compared to the Ensembl porcine BAC sequence.
Correlations between transcript abundance and meat quality, combined with genetic positional information of eQTL allowed us to prioritize a small number of candidate genes. We considered transcripts, which exhibited expression levels correlated with meat quality traits, and which had eQTL on the same chromosome as the transcripts itself.
Our previous studies showed that cis regulation is a stable characteristic of individual transcripts. Consequently, a global microarray eQTL analysis of a limited number of samples can be used for exploring functional and regulatory gene networks and scanning for cis-eQTL. In particular, the assignment of eQTL to chromosomes is reliable; though some cis-eQTL change their position, they were consistently assigned to the same chromosome when comparing analyses based on 74 microarrays or 276 real time RT-PCRs . Based on this observation, we decided to highlight only eQTL which were located on the gene itself or the same chromosome. Moreover, this relaxed measure was chosen in order to avoid the exclusion of any true cis-acting eQTL that are not precisely mapped, because of the limited resolution of the genetic map. The observation that cis-eQTL were more consistently detected than trans-eQTL was also made by comparing the results from different studies of eQTL with different numbers of animals and different tissue types [14-17].
Expression QTL mapping, with its potential to categorize cis and trans-effects, provides the mean to discriminate between "effect" and "cause" with respect to trait-associated differential expression. Though a relaxed window of cis-eQTLs was used, a low proportion of 10% of putative cis-eQTL was found in this study, compared to other previous studies [18,19]. Cis-regulated genes are of interest, because the underlying genes are expected to harbour genetic variants that influence their own expression level, which may also influence the physiological traits of interest, if transcript abundance is correlated with the target phenotype . This provides a function-related evidence of the candidacy of a particular gene.
The usefulness of eQTL for identification of quantitative trait genes was demonstrated [3,5,20,21]. The causal link between sequence variation, gene expression, and phenotype arises because the polymorphism might be responsible for both a cis-eQTL and the QTL contributing to a quantitative phenotype, so called pQTL. QTL for traits related to meat quality were previously mapped in the population used here [12,22]. The matching of pQTL, eQTL and the localization of the corresponding genes with trait-correlated expression provides positional and functional evidence for the potential role of the respective genes and strongly promote them as candidate genes. Mapping of eQTL enables displaying regulatory networks and localizing genomic variation affecting the mRNA expression of a gene either within the genes itself (cis) or distant from the gene (trans). The key advantage of eQTL mapping is that it connects variation at the level of RNA expression to variation at the level of DNA. Only the latter provides versatile tools for breeding whereas the first reveals information on the biology of a trait and directs to new candidate genes. Moreover, integration of (1) information on QTL for a trait of interest in breeding (pheneQTL = pQTL) with analysis of (2) trait correlated expression and with (3) mapping of expression QTL (eQTL) for the corresponding trait-dependent-regulated genes facilitates the identification of genes and pathways with cumulative evidence of their involvement in the biology of the traits of interest and enable to built priority lists of candidate genes . However, there are also some issues that limit of the use of the genetical genomics approach, in particular the resolution of the genetic maps that is depending on the number of markers and animals used and the structure of the population used and artifacts caused by the limited sensitivity and specificity of microarray experiments .
In a previous study we used the principle components with high loadings of meat quality traits to identify functional networks of genes and to gain knowledge of biological and physiological processes taking place during the conversion of muscle to meat . Here in this study, we combined principle component analysis of traits related to meat quality with eQTL to scale down the list of candidate genes for the traits.
Holistic expression profiling was integrated with QTL analyses for meat quality traits. This is, to our knowledge, the first report of a comprehensive scan for cis-eQTL associated with meat quality traits in the pig. Correlations between transcript abundance and meat quality traits, combined with genetic positional information of eQTL allowed us to prioritise candidate genes for further study. Accordingly, a list of candidate genes for meat quality was set up. The further identification of the causative polymorphisms and the determination of their functional role are even more challenging, since there are several different molecular mechanisms through which mRNA levels in cells can be regulated.
This study was based on trait measurements, genotyping procedures, expression profiling, and linkage analysis performed in the three-generation resource family (DuPi) founded by crossbreeding Duroc and Pietrain and described in detail by Liu et al. [12,22]. All animals were free of the mutation at the ryanodine receptor locus, RYR1, which is responsible for malignant hyperthermia syndrome. A total of 572 F2 animals comprising 31 full-sib families were used for recording of meat performance traits and construction of a linkage map comprising genotypes of 115 microsatellite markers. Expression profiling and eQTL analysis were done with 74 F2 animals of this resource population that represented a subset of the population covering 25 full-sib families derived from all five F1 boars of the population and 18 out of 27 F1 sows. The experimental research on animals was done according to the German Animal Welfare Act and approved by the animal welfare board of the Leibniz Institute of Farm Animal Biology, FBN Dummerstorf.
Phenotypic data of F2 animals were collected following the guidelines of the German performance test (ZDS 2003) . After slaughter technological parameters of meat quality, i.e. pH-value, conductivity and colour, were measured by using Star-series equipment (Rudolf Matthaeus Company, Germany). Measures of pH and conductivity were at 45 min post mortem (pH1) and 24 h post mortem (pH24), respectively; both in M. longissimus dorsi between 13th/14th rib (pH1ld, pH24ld, LF1ld, LF24ld) and in the ham (M. semimembranosus) (pH24sm, LF24sm), respectively. Muscle colour was measured at 24 h post mortem by Opto-Star. Drip loss was scored based on a bag-method with a size-standardized sample from the M. longissimus dorsi collected at 24 hours post mortem that was weighted, suspended in a plastic bag, held at 4°C for 48 h, and thereafter re-weighed [26,27]. To obtain cooking loss, a loin cube was taken from the M. longissimus dorsi, weighed, placed in a polyethylene bag and incubated in water at 75°C for 50 minutes. The bag was then immersed in flowing water at room temperature for 30 minutes and the solid portion was re-weighed. Thawing loss was determined similarly after at least 24 h of freezing at -20°C. Drip loss, cooking loss, and thawing loss were calculated as a percentage of weight loss based on the start weight of a sample. Shear force was measured by the Instron-4310 equipment and the average values of four replicates were used for analyses. The procedure 'Factor' of the SAS software package (SAS version 9.1 SAS Institute, Cary, NC) was used to derive four principle components (PCs) based on 19 traits related to muscle and carcass properties of 572 DUPI animals; i.e. the dimensionality of the data was reduced from the original 19 traits to 4 principal components (PCs), where each PC is a linear combination of all traits without a significant loss of information. PC1 and PC4 represent variation in carcass traits, whereas PC2 and PC3 represent aspects of meat quality, consequently PC2 and PC3 were used to correlate with expression profiles. The principal component analysis of meat quality traits is discussed in Ponsuksili et al., 2009 .
Immediately post mortem tissue samples were collected and snap frozen that were taken between the 13th and 14th rib from the center of M. longissimus dorsi of 74 animals. Total RNA was isolated using TRI Reagent (Sigma, Taufkirchen, Germany) according to the manufacturer's protocol. After DNaseI treatment the RNA was cleaned up using the RNeasy Kit (Qiagen, Hilden, Germany). The quantity of RNA was established using the NanoDrop ND-1000 spectrophotometer (Peqlab, Erlangen, Germany) and the integrity was checked by running 1 μg of RNA on 1% agarose gel. In addition absence of DNA contamination was checked using the RNA as a template in standard PCR amplifying fragments of RPL32 and HPRT1. Muscle expression pattern were assessed using 74 Porcine Genome Array which contains 24,123 probe-sets that interrogate 20,689 transcripts that were assigned to known genes . Preparation of target products, hybridization and scanning using the GeneChip scanner 3000 were performed according to Affymetrix protocols using 5 μg of total RNA to prepare antisense biotinylated RNA. The quality of hybridization was assessed in all samples following the manufacturer's recommendations. Data were analysed with Affymetrix GCOS 1.1.1 software using global scaling to a target signal of 500. Data were then imported into Arrays Assist software (Stratagene Europe, Amsterdam, The Netherlands) for subsequent analysis. The data were processed with MAS5.0 to generate cell intensity files (using default settings with detection P-values of < 0.04 for 'present', ≥0.04 and ≤ 0.06 for 'marginal', < 0.06 for 'absent'; only 'present' calls were used). Quantitative expression levels of the present transcripts were estimated using PLIER (Probe Logarithmic Intensity Error) for normalization. The microarray data related to all samples have been deposited in the Gene Expression Omnibus (GEO, ) public repository (GEO accession number: GSE10204).
Phenotypic data, i.e. expression levels and meat quality data as well as PC traits, were adjusted for systematic effects by analysis of variance performed with the procedure 'Mixed' of the SAS software package (SAS System for Windows, Release 9.02) before analysing their correlation and performing eQTL analysis. Full-sib family and sex were used as fixed effects, carcass weight as a covariate and slaughter date as random effect. Based on observations from 74 animals Pearson correlation coefficients were calculated between the residuals of the log2 transformed expression intensities of all 11,453 probes and each individual meat quality trait as well as between the expression levels and both composite traits, PC2 and PC3. For each pair of transcript level and phenotype, Pearson correlation together with the P-value was computed (significance threshold p ≤ 0.05). The correlation coefficients at p ≤ 0.05 ranged between |0.24-0.50|. The observed false discovery rates were 0.08-0.34 which is reasonable for a microarray study, in particular, considering the relatively relaxed p-value (p ≤ 0.05).
In order to map eQTL adjusted expression values of 11,457 probe-sets were regressed onto the additive and dominance coefficients in intervals of 1 cM using the F2 option of QTL express . Chromosome-wide significance levels were estimated by permutation tests using 5000 permutations . The 5% chromosome-wide threshold corresponds approximately to the suggestive linkage threshold proposed by Lander & Kruglyak (1995) . Average significance thresholds were 4.92 corresponding to LOD = 2.0 with nominal p ≤ 10-3. Annotations and localization of probe-sets was based on assembly Sscrofa9 (April 2009) . In order to link the genetic map to Ensembl Sscrofa9, all 115 microsatellites sequences, which were used in the linkage map, were BLASTed against the Ensembl porcine BAC sequence (Sscrofa9, April 2009). Only 73 microsatellites sequences could be localized on Sscrofa 9 as shown on Figure 2, 3, 4, 5, 6.
SP analyzed the microarray data and drafted the manuscript; EM, MS, and KS aided in data analysis and helped in drafting the manuscript; KW conceived and designed the study, contributed to data interpretation and helped in drafting the manuscript. All authors read and approved the final manuscript.
eQTL detected for probe-sets representing transcripts that are expressed in porcine M. longissimus dorsi. All 9,180 eQTL (at LOD score > 2) estimated for probe-sets of the Affymetrix Porcine Genome Array, their position in cM on linkage groups, corresponding gene names and chromosomal positions according to the current porcine genome sequence assembly (Sscrofa9 assembly, April 2009) 
Click here for additional data file (1471-2164-11-572-S1.XLS)
The authors thank Annette Jugert and Joana Bittner for excellent technical help. This research was supported by German Research Foundation (Deutsche Forschungs-gemeinschaft, DFG; Forschergruppe 'DRIP', FOR 753)
|Jansen RC,Nap JP,Genetical genomics: the added value from segregationTrends GenetYear: 20011738839110.1016/S0168-9525(01)02310-111418218|
|Jansen RC,Studying complex biological systems using multifactorial perturbationNat Rev GenetYear: 2003414515110.1038/nrg99612560811|
|Schadt EE,Lamb J,Yang X,Zhu J,Edwards S,Guhathakurta D,Sieberts SK,Monks S,Reitman M,Zhang C,Lum PY,Leonardson A,Thieringer R,Metzger JM,Yang L,Castle J,Zhu H,Kash SF,Drake TA,Sachs A,Lusis AJ,An integrative genomics approach to infer causal associations between gene expression and diseaseNat GenetYear: 20053771071710.1038/ng158915965475|
|Doss S,Schadt EE,Drake TA,Lusis AJ,Cis-acting expression quantitative trait loci in miceGenome ResYear: 20051568169110.1101/gr.321690515837804|
|Ponsuksili S,Jonas E,Murani E,Phatsara C,Srikanchai T,Walz C,Schwerin M,Schellander K,Wimmers K,Trait correlated expression combined with expression QTL analysis reveals biological pathways and candidate genes affecting water holding capacity of muscleBMC GenomicsYear: 2008936710.1186/1471-2164-9-36718671879|
|Naraballobh W,Chomdej S,Murani E,Wimmers K,Ponsuksili S,Annotation and in silico localization of the Affymetrix GeneChip Porcine Genome ArrayArch TierzYear: 201053230238|
|Ponsuksili S,Murani E,Phatsara C,Schwerin M,Schellander K,Wimmers K,Expression quantitative trait loci analysis of genes in porcine muscle by quantitative real time RT-PCR compared to microarray dataHeredityYear: 2010 in press . 20145673|
|Affymetrix Technical ManualAffymetrix GeneChip Expression AnalysisYear: 2001Affymetrix, Santa Clara, CA|
|Affymetrix Technical NoteGuide to Probe Logarithmic Intensity Error (PLIER) EstimationYear: 2005Affymetrix, Santa Clara, CA|
|Lander E,Kruglyak L,Genetic dissection of complex traits: guidelines for interpreting and reporting linkage resultsNat GenetYear: 19951124124710.1038/ng1195-2417581446|
|Hodgson RR,Belk KE,Savell JW,Cross HR,Williams FL,Development of a quantitative quality grading system for mature cow carcassesJ Anim SciYear: 199270184018471634408|
|Liu G,Jennen DG,Tholen E,Juengst H,Kleinwächter T,Hölker M,et al. A genome scan reveals QTL for growth, fatness, leanness and meat quality in a Duroc-Pietrain resource populationAnim GenetYear: 20073824125210.1111/j.1365-2052.2007.01592.x17459017|
|Ponsuksili S,Murani E,Phatsara C,Schwerin M,Schellander K,Wimmers K,Porcine muscle sensory attributes associate with major changes in gene networks involving CAPZB, ANKRD1, and CTBP2Funct Integr GenomicsYear: 2009945547110.1007/s10142-009-0131-119597856|
|Brem RB,Yvert G,Clinton R,Kruglyak L,Genetic dissection of transcriptional regulation in budding yeastScienceYear: 200229675275510.1126/science.106951611923494|
|Yvert G,Brem RB,Whittle J,Akey JM,Foss E,Smith EN,et al. Trans-acting regulatory variation in Saccharomyces cerevisiae and the role of transcription factorsNat GenetYear: 200335576410.1038/ng122212897782|
|Morley M,Molony CM,Weber TM,Devlin JL,Ewens KG,Spielman RS,et al. Genetic analysis of genome-wide variation in human gene expressionNatureYear: 200443074374710.1038/nature0279715269782|
|Göring HH,Curran JE,Johnson MP,Dyer TD,Charlesworth J,Cole SA,et al. Discovery of expression QTLs using large-scale transcriptional profiling in human lymphocytesNat GenetYear: 2007391208121610.1038/ng211917873875|
|Schadt EE,Monks SA,Drake TA,Lusis AJ,Che N,Colinayo V,Ruff TG,Milligan SB,Lamb JR,Cavet G,Linsley PS,Mao M,Stoughton RB,Friend SH,Genetics of gene expression surveyed in maize, mouse and manNatureYear: 200342229730110.1038/nature0143412646919|
|Hubner N,Wallace CA,Zimdahl H,Petretto E,Schulz H,Maciver F,Mueller M,Hummel O,Monti J,Zidek V,Musilova A,Kren V,Causton H,Game L,Born G,Schmidt S,Müller A,Cook SA,Kurtz TW,Whittaker J,Pravenec M,Aitman TJ,Integrated transcriptional profiling and linkage analysis for identification of genes underlying diseaseNat GenetYear: 20053724325310.1038/ng152215711544|
|Bystrykh L,Weersing E,Dontje B,Sutton S,Pletcher MT,Wiltshire T,et al. Uncovering regulatory pathways that affect hematopoietic stem cell function using 'genetical genomics'Nat GenetYear: 20053722523210.1038/ng149715711547|
|Mehrabian M,Allayee H,Stockton J,Lum PY,Drake TA,Castellani LW,Suh M,Armour C,Edwards S,Lamb J,Lusis AJ,Schadt EE,Integrating genotypic and expression data in a segregating mouse population to identify 5-lipoxygenase as a susceptibility gene for obesity and bone traitsNat GenetYear: 2005371224123310.1038/ng161916200066|
|Liu G,Kim JJ,Jonas E,Wimmers K,Ponsuksili S,Murani E,Phatsara C,Tholen E,Jüngst H,Tesfaye D,Chen JL,Schellander K,Combined line-cross and half-sib QTL analysis in Duroc-Pietrain populationMamm GenomeYear: 2008642943810.1007/s00335-008-9132-y|
|Wimmers K,Murani E,Ponsuksili S,Functional genomics and genetical genomics approaches towards elucidating networks of genes affecting meat performance in pigsBrief Funct GenomicsYear: 2010925125820211968|
|Verdugo A,Farber CR,Warden CH,Medrano JF,Serious limitations of the QTL/microarray approach for QTL gene discoveryBMC BiologyYear: 201089610.1186/1741-7007-8-9620624276|
|ZDS, Zentral Verband der Deutschen Schweineproduktion e. VRichtlinie fuer die Stationspruefung auf Mastleistung, Schlachtkoerperwert und Fleischbeschaffenheit beim SchweinYear: 2003Bonn, Germany|
|Honikel KO,Wasserbindungsvermögen von FleischMitteilungsblatt der BAFFYear: 1986671507154|
|Kauffman RG,Eikelenboom G,Wal van der PG,Merkus G,Zaar M,The use of filter paper to estimate drip loss of porcine musculatureMeat SciYear: 19861819120010.1016/0309-1740(86)90033-1|
|Gene Expression Omnibushttp://www.ncbi.nlm.nih.gov/geo/|
|Seaton G,Haley CS,Knott SA,Kearsey M,Visscher PM,QTL Express: mapping quantitative trait loci in simple and complex pedigreesBioinformaticsYear: 20021833934010.1093/bioinformatics/18.2.33911847090|
|Churchill GA,Doerge RW,Empirical threshold values for quantitative trait mappingGeneticsYear: 19941389639717851788|
Global eQTL summary
|Total # of probe-sets||Lod score threshold||Nominal P||Total # of eQTL||# of probe-sets with at least one eQTL||# of eQTL mapped||# of eQTL with ipsi-chromosomal location*|
*, i.e. # of eQTL located on the same chromosome as the corresponding probe-set
The number of eQTL of probe-sets exhibiting expression levels that are correlated with meat quality traits
|No. of eQTL with trait correlated express gene|
|Meat quality parameters||Traits||All eQTL||eQTL located on the same chromosome as the corresponding probe-set|
|- conductivity||LF1h loin||253||16|
|technological meat quality parameters||QPTO||273||12|
pH1ld pH-value in Mld at 13th/14th rib 45 min post mortem
pH24ld pH-value in Mld at 13th/14th rib 24 h post mortem
pH24sm pH-value in M. semimembranosus (Msm) at 24 h post mortem
LF1ld conductivity in Mld at 13th/14th rib45 min post mortem
LF24ld conductivity in Mld at 13th/14th rib post mortem
LF24sm conductivity in Msm at 24 h post mortem
DL Drip loss: percentage of weight loss after 48 h of Mld samples collected at 24 h post mortem
TL Thaw loss: percentage of weight loss of Mld samples frozen at -20°C.
CL Cooking loss: percentage of weight loss of Mld samples incubated in water at 75°C for 50 min
OPTO meat colour 24 h post mortem in M. longissimus dorsi (Mld)at 13th/14th rib; OPTO star
SF Shear force was measured by the Instron-4310 equipment
Probe-sets with expression levels that are correlated with principal component 2 and information about their eQTL.
|Annotation||eQTL||individual trait correlation*|
|Ssc.25132.3.S1||RAD23B||1||130||5.11||DL, -pH1ld, LF1ld, LF24sm|
|Ssc.28435.1.A1||EHD1||2||56||8.21||CL, -pH24ld, -pH24sm, LF24sm|
|Ssc.6529.1.A1||AFF4||2||132||9.3||pH24ld, pH24sm, -LF24ld|
|Ssc.15224.2.S1||FOSL2||3||86||6.2||-DL, pH24ld, pH24sm|
|Ssc.15224.1.S1||FOSL2||3||86||6.3||-DL, -SF, -TL, pH24ld, pH24sm|
|Ssc.26780.1.S1||BCL2L12||6||0||6.9||DL, -pH24ld, -pH24sm, LF24sm,|
|Ssc.10297.3.S1||CAPZB||6||95||6.2||DL, CL, TL, -pH24ld, LF24ld, LF24sm|
|Ssc.4217.1.S1||ITIH4||13||47||7.3||-DL, pH1ld, -LF24sm|
|Ssc.10429.1.S1||ANKRD1||14||0||6||-DL, pH24ld, -LF24sm|
|Ssc.10911.1.A1||ADRBK2||14||0||6||-CL, pH24ld, -LF24ld|
|Ssc.6058.1.S1||LOC100129026||16||55||4.9||-DL, -SF, -LF24sm|
"*"the probe-sets listed here represent a subset of probe-sets correlated with individual meat quality traits as indicated in this column naming the traits with which the corresponding probe-sets show correlated expression in addition to the correlated expression with PC2; "-" indicates negative correlation
Trait abbreviations are given at table 2.
Probe-sets with expression levels that are correlated with principal component 3 and information about their eQTL
|Annotation||eQTL||individual trait correlation*|
|Ssc.29222.1.S1||UBE2CBP||1||31||5.72||pH1ld, -LF1ld, -LF24ld|
|Ssc.3853.1.S1||C9orf61||1||54||7.98||-DL, pH1ld, -LF24sm|
|Ssc.7558.1.A1||1||127||7.96||-SF, pH1ld, -LF24ld|
|Ssc.25132.3.S1||RAD23B||1||130||5.11||DL, -pH1ld, LF1ld, LF24sm|
|Ssc.5415.1.S1||DDX1||3||86||4.96||pH1ld, -pH24ld, -LF24ld|
|Ssc.10297.3.S1||CAPZB||6||95||6.17||DL, CL, TL, -pH24ld, LF24ld, LF24sm|
|Ssc.25358.1.S1||GRAMD2||7||5||5.41||-CL, pH1ld, -LF24ld|
|Ssc.4217.1.S1||ITIH4||13||47||7.31||-DL, pH1ld, -LF24sm|
"*"the probe-sets listed here represent a subset of probe-sets correlated with individual meat quality traits as indicated in this column naming the traits with which the corresponding probe-sets show correlated expression in addition to the correlated expression with PC3; "-" indicates negative correlation
Trait abbreviations are given at table 2.
Previous Document: Autoantibodies to muscarinic acetylcholine receptors found in patients with primary biliary cirrhosi...
Next Document: Biological changes in auditory function following training in children with autism spectrum disorder...