Impact of breed and sex on porcine endocrine transcriptome: a bayesian biometrical analysis

Changes in Biceps femoris Transcriptome along Growth in Iberian Pigs Fed Different Energy Sources and Comparative Analysis with Duroc Breed

Simple Summary

The genetic mechanisms that regulate biological processes, such as skeletal muscle development and growth, or intramuscular fat deposition, have attracted great interest, given their impact on production traits and meat quality. In this sense, a comparison of the transcriptome of skeletal muscle between phenotypically different pig breeds, or along growth, could be useful to improve the understanding of the molecular processes underlying the differences in muscle metabolism and phenotypic traits, potentially driving the identification of causal genes, regulators and metabolic pathways involved in their variability.

Abstract

This experiment was conducted to investigate the effects of developmental stage, breed, and diet energy source on the genome-wide expression, meat quality traits, and tissue composition of biceps femoris muscle in growing pure Iberian and Duroc pigs. The study comprised 59 Iberian (IB) and 19 Duroc (DU) animals, who started the treatment at an average live weight (LW) of 19.9 kg. The animals were kept under identical management conditions and fed two diets with different energy sources (6% high oleic sunflower oil or carbohydrates). Twenty-nine IB animals were slaughtered after seven days of treatment at an average LW of 24.1 kg, and 30 IB animals plus all the DU animals were slaughtered after 47 days at an average LW of 50.7 kg. The main factors affecting the muscle transcriptome were age, with 1832 differentially expressed genes (DEGs), and breed (1055 DEGs), while the effect of diet on the transcriptome was very small. The results indicated transcriptome changes along time in Iberian animals, being especially related to growth and tissue development, extracellular matrix (ECM) composition, and cytoskeleton organization, with DEGs affecting relevant functions and biological pathways, such as myogenesis. The breed also affected functions related to muscle development and cytoskeleton organization, as well as functions related to solute transport and lipid and carbohydrate metabolism. Taking into account the results of the two main comparisons (age and breed effects), we can postulate that the Iberian breed is more precocious than the Duroc breed, regarding myogenesis and muscle development, in the studied growing stage.

Conservation Genomic Analysis of the Croatian Indigenous Black Slavonian and Turopolje Pig Breeds

The majority of the nearly 400 existing local pig breeds are adapted to specific environments and human needs. The demand for large production quantities and the industrialized pig production have caused a rapid decline of many local pig breeds in recent decades. Black Slavonian pig and Turopolje pig, the latter highly threatened, are the two Croatian local indigenous breeds typically grown in extensive or semi-intensive systems. In order to guide a long-term breeding program to prevent the disappearance of these breeds, we analyzed their genetic diversity, inbreeding level and relationship with other local breeds across the world, as well as modern breeds and several wild populations, using high throughput genomic data obtained using the Illumina Infinium PorcineSNP60 v2 BeadChip. Multidimensional scaling analysis positioned Black Slavonian pigs close to the UK/North American breeds, while the Turopolje pig clustered within the Mediterranean breeds. Turopolje pig showed a very high inbreeding level (F_ROH>4Mb = 0.400 and F_ROH>8Mb = 0.332) that considerably exceeded the level of full-sib mating, while Black Slavonian pig showed much lower inbreeding (F_ROH>4Mb = 0.098 and F_ROH>8Mb = 0.074), indicating a planned mating strategy. In Croatian local breeds we identified several genome regions showing adaptive selection signals that were not present in commercial breeds. The results obtained in this study reflect the current genetic status and breeding management of the two Croatian indigenous local breeds. Given the small populations of both breeds, a controlled management activity has been implemented in Black Slavonian pigs since their commercial value has been recognized. In contrast, the extremely high inbreeding level observed in Turopolje pig argues for an urgent conservation plan with a long-term, diversity-oriented breeding program.

Breed, Diet, and Interaction Effects on Adipose Tissue Transcriptome in Iberian and Duroc Pigs Fed Different Energy Sources

In this study, we analyzed the effects of breed, diet energy source, and their interaction on adipose tissue transcriptome in growing Iberian and Duroc pigs. The study comprised 29 Iberian and 19 Duroc males, which were kept under identical management conditions except the nutritional treatment. Two isoenergetic diets were used with 6% high oleic sunflower oil (HO) or carbohydrates (CH) as energy sources. All animals were slaughtered after 47 days of treatment at an average live weight of 51.2 kg. Twelve animals from each breed (six fed each diet) were employed for ham subcutaneous adipose tissue RNA-Seq analysis. The data analysis was performed using two different bioinformatic pipelines. We detected 837 and 1456 differentially expressed genes (DEGs) according to breed, depending on the pipeline. Due to the strong effect of breed on transcriptome, the effect of the diet was separately evaluated in the two breeds. We identified 207 and 57 DEGs depending on diet in Iberian and Duroc pigs, respectively. A joint analysis of both effects allowed the detection of some breed–diet interactions on transcriptome, which were inferred from RNA-Seq and quantitative PCR data. The functional analysis showed the enrichment of functions related to growth and tissue development, inflammatory response, immune cell trafficking, and carbohydrate and lipid metabolism, and allowed the identification of potential regulators. The results indicate different effects of diet on adipose tissue gene expression between breeds, affecting relevant biological pathways.

Modulatory Effects of Breed, Feeding Status, and Diet on Adipogenic, Lipogenic, and Lipolytic Gene Expression in Growing Iberian and Duroc Pigs

Meat quality depends on tissue composition which is in turn influenced by different factors, such as diet, genotype, age, or sex. We evaluated the effects of breed, 24 h fasting, and dietary energy source (HO: oleic acid versus CH: carbohydrates) on the expression of candidate genes involved in adipogenesis, lipogenesis, and lipolysis in the adipose tissue from Iberian and Duroc growing pigs. The Iberian pigs showed greater feed intake, backfat thickness, and saturated fatty acids (SFA) content in the subcutaneous fat, whereas the Duroc pigs had greater ham weight and polyunsaturated fatty acids (PUFA) content. In both breeds, the diet induced changes in the fatty acid (FA) composition of subcutaneous fat samples. The HO group had higher monounsaturated fatty acids (MUFA) and oleic acid, and lower SFA than the CH group. Regarding gene expression, breed and feeding status (fasting versus postprandial) had significant effects on gene expression, with quantitative interactions between them, while diet showed negligible effects. In general, adipogenic and lipogenic genes were upregulated in the Iberian pigs and in postprandial samples. In contrast, the expression of lipolytic genes showed complex interaction effects. Our results agree with the phenotypic differences between the Iberian and Duroc breeds and with the inhibition of lipogenesis by fasting. Quantitative interactions between breed and feeding status effects were observed, which indicates a different response to fasting of the two breeds, with the obese Iberian breed showing a more stable expression of lipogenic genes. These results highlight the complexity of lipid metabolism regulation, especially in relation to lipolysis processes.

Developmental Stage, Muscle and Genetic Type Modify Muscle Transcriptome in Pigs: Effects on Gene Expression and Regulatory Factors Involved in Growth and Metabolism

Iberian pig production includes purebred (IB) and Duroc-crossbred (IBxDU) pigs, which show important differences in growth, fattening and tissue composition. This experiment was conducted to investigate the effects of genetic type and muscle (Longissimus dorsi (LD) vs Biceps femoris (BF)) on gene expression and transcriptional regulation at two developmental stages. Nine IB and 10 IBxDU piglets were slaughtered at birth, and seven IB and 10 IBxDU at four months of age (growing period). Carcass traits and LD intramuscular fat (IMF) content were measured. Muscle transcriptome was analyzed on LD samples with RNA-Seq technology. Carcasses were smaller in IB than in IBxDU neonates (p < 0.001), while growing IB pigs showed greater IMF content (p < 0.05). Gene expression was affected (p < 0.01 and Fold change > 1.5) by the developmental stage (5,812 genes), muscle type (135 genes), and genetic type (261 genes at birth and 113 at growth). Newborns transcriptome reflected a highly proliferative developmental stage, while older pigs showed upregulation of catabolic and muscle functioning processes. Regarding the genetic type effect, IBxDU newborns showed enrichment of gene pathways involved in muscle growth, in agreement with the higher prenatal growth observed in these pigs. However, IB growing pigs showed enrichment of pathways involved in protein deposition and cellular growth, supporting the compensatory gain experienced by IB pigs during this period. Moreover, newborn and growing IB pigs showed more active glucose and lipid metabolism than IBxDU pigs. Moreover, LD muscle seems to have more active muscular and cell growth, while BF points towards lipid metabolism and fat deposition. Several regulators controlling transcriptome changes in both genotypes were identified across muscles and ages (SIM1, PVALB, MEFs, TCF7L2 or FOXO1), being strong candidate genes to drive expression and thus, phenotypic differences between IB and IBxDU pigs. Many of the identified regulators were known to be involved in muscle and adipose tissues development, but others not previously associated with pig muscle growth were also identified, as PVALB, KLF1 or IRF2. The present study discloses potential molecular mechanisms underlying phenotypic differences observed between IB and IBxDU pigs and highlights candidate genes implicated in these molecular mechanisms.

Comparative Analysis of Muscle Transcriptome between Pig Genotypes Identifies Genes and Regulatory Mechanisms Associated to Growth, Fatness and Metabolism

Iberian ham production includes both purebred (IB) and Duroc-crossbred (IBxDU) Iberian pigs, which show important differences in meat quality and production traits, such as muscle growth and fatness. This experiment was conducted to investigate gene expression differences, transcriptional regulation and genetic polymorphisms that could be associated with the observed phenotypic differences between IB and IBxDU pigs. Nine IB and 10 IBxDU pigs were slaughtered at birth. Morphometric measures and blood samples were obtained and samples from Biceps femoris muscle were employed for compositional and transcriptome analysis by RNA-Seq technology. Phenotypic differences were evident at this early age, including greater body size and weight in IBxDU and greater Biceps femoris intramuscular fat and plasma cholesterol content in IB newborns. We detected 149 differentially expressed genes between IB and IBxDU neonates (p < 0.01 and Fold-Change > 1. 5). Several were related to adipose and muscle tissues development (DLK1, FGF21 or UBC). The functional interpretation of the transcriptomic differences revealed enrichment of functions and pathways related to lipid metabolism in IB and to cellular and muscle growth in IBxDU pigs. Protein catabolism, cholesterol biosynthesis and immune system were functions enriched in both genotypes. We identified transcription factors potentially affecting the observed gene expression differences. Some of them have known functions on adipogenesis (CEBPA, EGRs), lipid metabolism (PPARGC1B) and myogenesis (FOXOs, MEF2D, MYOD1), which suggest a key role in the meat quality differences existing between IB and IBxDU hams. We also identified several polymorphisms showing differential segregation between IB and IBxDU pigs. Among them, non-synonymous variants were detected in several transcription factors as PPARGC1B and TRIM63 genes, which could be associated to altered gene function. Taken together, these results provide information about candidate genes, metabolic pathways and genetic polymorphisms potentially involved in phenotypic differences between IB and IBxDU pigs associated to meat quality and production traits.

Exploring transcriptomic diversity in muscle revealed that cellular signaling pathways mainly differentiate five Western porcine breeds

Background

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.

Results

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.

Conclusions

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.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-015-2259-9) contains supplementary material, which is available to authorized users.

Longissimus dorsi transcriptome analysis of purebred and crossbred Iberian pigs differing in muscle characteristics

BACKGROUND

The two main genetic types in Iberian pig production show important phenotypic differences in growth, fattening and tissue composition since early developmental stages. The objective of this work was the evaluation of muscle transcriptome profile in piglets of both genetic types, in order to identify genes, pathways and regulatory factors responsible for their phenotypic differences. Contemporary families coming from pure Iberian pigs (IB) or from crossing with Duroc boars (DU×IB) were generated. Piglets (14 from each genetic type) were slaughtered at weaning (28 days) and longissimus dorsi was sampled for composition and gene expression studies. RNA was obtained and hybridized to Affymetrix Porcine Genechip expression arrays.

RESULTS

Loin muscle chemical composition showed significant differences between genetic types in intramuscular fat content (6.1% vs. 4.3% in IB and DUxIB animals, respectively, P = 0.009) and in saturated (P = 0.019) and monounsaturated fatty acid proportions (P = 0.044). The statistical analysis of gene expression data allowed the identification of 256 differentially expressed (DE) genes between genetic types (FDR < 0.10), 102 upregulated in IB and 154 upregulated in DU×IB. Transcript differences were validated for a subset of DE genes by qPCR. We observed alteration in biological functions related to extracellular matrix function and organization, cellular adhesion, muscle growth, lipid metabolism and proteolysis. Candidate genes with known effects on muscle growth were found among the DE genes upregulated in DU×IB. Genes related to lipid metabolism and proteolysis were found among those upregulated in IB. Regulatory factors (RF) potentially involved in the expression differences were identified by calculating the regulatory impact factors. Twenty-nine RF were found, some of them with known relationship with tissue development (MSTN, SIX4, IRX3), adipogenesis (CEBPD, PPARGC1B), or extracellular matrix processes (MAX, MXI1). Correlation among the expression of these RF and DE genes show relevant differences between genetic types.

CONCLUSION

These results provide valuable information about genetic mechanisms determining the phenotypic differences on growth and meat quality between the genetic types studied, mainly related to the development and function of the extracellular matrix and also to some metabolic processes as proteolysis and lipid metabolism. Transcription factors and regulatory mechanisms are proposed for these altered biological functions.

Identification of genes regulating growth and fatness traits in pig through hypothalamic transcriptome analysis

Previous studies on Iberian × Landrace (IBMAP) pig intercrosses have enabled the identification of several quantitative trait locus (QTL) regions related to growth and fatness traits; however, the genetic variation underlying those QTLs are still unknown. These traits are not only relevant because of their impact on economically important production traits, but also because pig constitutes a widely studied animal model for human obesity and obesity-related diseases. The hypothalamus is the main gland regulating growth, food intake, and fat accumulation. Therefore, the aim of this work was to identify genes and/or gene transcripts involved in the determination of growth and fatness in pig by a comparison of the whole hypothalamic transcriptome (RNA-Seq) in two groups of phenotypically divergent IBMAP pigs. Around 16,000 of the ∼25.010 annotated genes were expressed in these hypothalamic samples, with most of them showing intermediate expression levels. Functional analyses supported the key role of the hypothalamus in the regulation of growth, fat accumulation, and energy expenditure. Moreover, 58,927 potentially new isoforms were detected. More than 250 differentially expressed genes and novel transcript isoforms were identified between the two groups of pigs. Twenty-one DE genes/transcripts that colocalized in previously identified QTL regions and/or whose biological functions are related to the traits of interest were explored in more detail. Additionally, the transcription factors potentially regulating these genes and the subjacent networks and pathways were also analyzed. This study allows us to propose strong candidate genes for growth and fatness based on expression patterns, genomic location, and network interactions.

First survey and functional annotation of prohormone and convertase genes in the pig

Background

The pig is a biomedical model to study human and livestock traits. Many of these traits are controlled by neuropeptides that result from the cleavage of prohormones by prohormone convertases. Only 45 prohormones have been confirmed in the pig. Sequence homology can be ineffective to annotate prohormone genes in sequenced species like the pig due to the multifactorial nature of the prohormone processing. The goal of this study is to undertake the first complete survey of prohormone and prohormone convertases genes in the pig genome. These genes were functionally annotated based on 35 gene expression microarray experiments. The cleavage sites of prohormone sequences into potentially active neuropeptides were predicted.

Results

We identified 95 unique prohormone genes, 2 alternative calcitonin-related sequences, 8 prohormone convertases and 1 cleavage facilitator in the pig genome 10.2 assembly and trace archives. Of these, 11 pig prohormone genes have not been reported in the UniProt, UniGene or Gene databases. These genes are intermedin, cortistatin, insulin-like 5, orexigenic neuropeptide QRFP, prokineticin 2, prolactin-releasing peptide, parathyroid hormone 2, urocortin, urocortin 2, urocortin 3, and urotensin 2-related peptide. In addition, a novel neuropeptide S was identified in the pig genome correcting the previously reported pig sequence that is identical to the rabbit sequence. Most differentially expressed prohormone genes were under-expressed in pigs experiencing immune challenge relative to the un-challenged controls, in non-pregnant relative to pregnant sows, in old relative to young embryos, and in non-neural relative to neural tissues. The cleavage prediction based on human sequences had the best performance with a correct classification rate of cleaved and non-cleaved sites of 92% suggesting that the processing of prohormones in pigs is similar to humans. The cleavage prediction models did not find conclusive evidence supporting the production of the bioactive neuropeptides urocortin 2, urocortin 3, torsin family 2 member A, tachykinin 4, islet amyloid polypeptide, and calcitonin receptor-stimulating peptide 2 in the pig.

Conclusions

The present genomic and functional characterization supports the use of the pig as an effective animal model to gain a deeper understanding of prohormones, prohormone convertases and neuropeptides in biomedical and agricultural research.

Porcine tissue-specific regulatory networks derived from meta-analysis of the transcriptome

The processes that drive tissue identity and differentiation remain unclear for most tissue types. So are the gene networks and transcription factors (TF) responsible for the differential structure and function of each particular tissue, and this is particularly true for non model species with incomplete genomic resources. To better understand the regulation of genes responsible for tissue identity in pigs, we have inferred regulatory networks from a meta-analysis of 20 gene expression studies spanning 480 Porcine Affymetrix chips for 134 experimental conditions on 27 distinct tissues. We developed a mixed-model normalization approach with a covariance structure that accommodated the disparity in the origin of the individual studies, and obtained the normalized expression of 12,320 genes across the 27 tissues. Using this resource, we constructed a network, based on the co-expression patterns of 1,072 TF and 1,232 tissue specific genes. The resulting network is consistent with the known biology of tissue development. Within the network, genes clustered by tissue and tissues clustered by site of embryonic origin. These clusters were significantly enriched for genes annotated in key relevant biological processes and confirm gene functions and interactions from the literature. We implemented a Regulatory Impact Factor (RIF) metric to identify the key regulators in skeletal muscle and tissues from the central nervous systems. The normalization of the meta-analysis, the inference of the gene co-expression network and the RIF metric, operated synergistically towards a successful search for tissue-specific regulators. Novel among these findings are evidence suggesting a novel key role of ERCC3 as a muscle regulator. Together, our results recapitulate the known biology behind tissue specificity and provide new valuable insights in a less studied but valuable model species.

Mitochondrial function, fatty acid metabolism, and immune system are relevant features of pig adipose tissue development

The molecular mechanisms underlying the genetic control of fat development in humans and livestock species still require characterization. To gain insights on gene expression patterns associated with genetic propensity for adiposity, we compared subcutaneous adipose tissue (SCAT) transcriptomics profiles from two contrasted pig breeds for body fatness. Samples were obtained from Large White (LW; lean phenotype) and Basque pigs (B; low growth and high fat content) at 35 kg (n = 5 per breed) or 145 kg body weight (n = 10 per breed). Using a custom adipose tissue microarray, we found 271 genes to be differentially expressed between the two breeds at both stages, out of which 123 were highly expressed in LW pigs and 148 genes were highly expressed in B pigs. Functional enrichment analysis based on gene ontology (GO) terms highlighted gene groups corresponding to the mitochondrial energy metabolism in LW pigs, whereas immune response was found significantly enriched in B pigs. Genes associated with lipid metabolism, such as ELOVL6, a gene involved in fatty acid elongation, had a lower expression in B compared with LW pigs. Furthermore, despite enlarged adipocyte diameters and higher plasma leptin concentration, B pigs displayed reduced lipogenic enzyme activities compared with LW pigs at 145 kg. Altogether, our results suggest that the development of adiposity was associated with a progressive worsening of the metabolic status, leading to a low-grade inflammatory state, and may thus be of significant interest for both livestock production and human health.

Exploring the gonad transcriptome of two extreme male pigs with RNA-seq

BACKGROUND

Although RNA-seq greatly advances our understanding of complex transcriptome landscapes, such as those found in mammals, complete RNA-seq studies in livestock and in particular in the pig are still lacking. Here, we used high-throughput RNA sequencing to gain insight into the characterization of the poly-A RNA fraction expressed in pig male gonads. An expression analysis comparing different mapping approaches and detection of allele specific expression is also discussed in this study.

RESULTS

By sequencing testicle mRNA of two phenotypically extreme pigs, one Iberian and one Large White, we identified hundreds of unannotated protein-coding genes (PcGs) in intergenic regions, some of them presenting orthology with closely related species. Interestingly, we also detected 2047 putative long non-coding RNA (lncRNA), including 469 with human homologues. Two methods, DEGseq and Cufflinks, were used for analyzing expression. DEGseq identified 15% less expressed genes than Cufflinks, because DEGseq utilizes only unambiguously mapped reads. Moreover, a large fraction of the transcriptome is made up of transposable elements (14500 elements encountered), as has been reported in previous studies. Gene expression results between microarray and RNA-seq technologies were relatively well correlated (r = 0.71 across individuals). Differentially expressed genes between Large White and Iberian showed a significant overrepresentation of gamete production and lipid metabolism gene ontology categories. Finally, allelic imbalance was detected in ~ 4% of heterozygous sites.

CONCLUSIONS

RNA-seq is a powerful tool to gain insight into complex transcriptomes. In addition to uncovering many unnanotated genes, our study allowed us to determine that a considerable fraction is made up of long non-coding transcripts and transposable elements. Their biological roles remain to be determined in future studies. In terms of differences in expression between Large White and Iberian pigs, these were largest for genes involved in spermatogenesis and lipid metabolism, which is consistent with phenotypic extreme differences in prolificacy and fat deposition between these two breeds.

Association between plasma metabolites and gene expression profiles in five porcine endocrine tissues

Background

Endocrine tissues play a fundamental role in maintaining homeostasis of plasma metabolites such as non-esterified fatty acids and glucose, the levels of which reflect the energy balance or the health status of animals. However, the relationship between the transcriptome of endocrine tissues and plasma metabolites has been poorly studied.

Methods

We determined the blood levels of 12 plasma metabolites in 27 pigs belonging to five breeds, each breed consisting of both females and males. The transcriptome of five endocrine tissues i.e. hypothalamus, adenohypophysis, thyroid gland, gonads and backfat tissues from 16 out of the 27 pigs was also determined. Sex and breed effects on the 12 plasma metabolites were investigated and associations between genes expressed in the five endocrine tissues and the 12 plasma metabolites measured were analyzed. A probeset was defined as a quantitative trait transcript (QTT) when its association with a particular metabolic trait achieved a nominal P value < 0.01.

Results

A larger than expected number of QTT was found for non-esterified fatty acids and alanine aminotransferase in at least two tissues. The associations were highly tissue-specific. The QTT within the tissues were divided into co-expression network modules enriched for genes in Kyoto Encyclopedia of Genes and Genomes or gene ontology categories that are related to the physiological functions of the corresponding tissues. We also explored a multi-tissue co-expression network using QTT for non-esterified fatty acids from the five tissues and found that a module, enriched in hypothalamus QTT, was positioned at the centre of the entire multi-tissue network.

Conclusions

These results emphasize the relationships between endocrine tissues and plasma metabolites in terms of gene expression. Highly tissue-specific association patterns suggest that candidate genes or gene pathways should be investigated in the context of specific tissues.

Regulatory impact factors: unraveling the transcriptional regulation of complex traits from expression data

MOTIVATION

Although transcription factors (TF) play a central regulatory role, their detection from expression data is limited due to their low, and often sparse, expression. In order to fill this gap, we propose a regulatory impact factor (RIF) metric to identify critical TF from gene expression data.

RESULTS

To substantiate the generality of RIF, we explore a set of experiments spanning a wide range of scenarios including breast cancer survival, fat, gonads and sex differentiation. We show that the strength of RIF lies in its ability to simultaneously integrate three sources of information into a single measure: (i) the change in correlation existing between the TF and the differentially expressed (DE) genes; (ii) the amount of differential expression of DE genes; and (iii) the abundance of DE genes. As a result, RIF analysis assigns an extreme score to those TF that are consistently most differentially co-expressed with the highly abundant and highly DE genes (RIF1), and to those TF with the most altered ability to predict the abundance of DE genes (RIF2). We show that RIF analysis alone recovers well-known experimentally validated TF for the processes studied. The TF identified confirm the importance of PPAR signaling in adipose development and the importance of transduction of estrogen signals in breast cancer survival and sexual differentiation. We argue that RIF has universal applicability, and advocate its use as a promising hypotheses generating tool for the systematic identification of novel TF not yet documented as critical.