BACKGROUND
psoriasis vulgaris is one of the most prevalent chronic inflammatory skin diseases affecting approximately 2% of individuals in western societies, and found worldwide in all populations.
psoriasis is a complex disease affecting cellular, gene and protein levels and presented as skin lesions.
the skin lesions are characterized by abnormal keratinocyte differentiation, hyperproliferation of keratinocytes, and infiltration of inflammatory cells  <cit> .
the mechanisms of psoriasis pathology are complex and involve genetic and environmental factors.
as we gain more knowledge about molecular pathways implicated in the disease, novel therapies emerge .
in recent years, microarray mrna expression profiling  <cit>  of lesional psoriatic skin revealed over  <dig>  differentially expressed genes.
enrichment analysis  showed that these genes encode proteins involved in regeneration, hyperkeratosis, metabolic function, immune response, and inflammation and revealed a number of modulating signaling pathways.
these efforts may help to develop new-generation drugs.
however, enrichment analysis limits our understanding of altered molecular interactions in psoriasis as it provides a relative ranking based on ontology terms resulting in the representation of fragmented and disconnected perturbed pathways.
furthermore, analysis of gene expression alone is not sufficient for understanding the whole variety of pathological changes at different levels of cellular organization.
indeed, new methodologies have been applied to the analysis of omics data in complex diseases that include algorithm-based biological network analysis  <cit>  and meta-analysis of multiple datasets of different types  <cit> .
here, we applied several techniques of network and meta-analysis to reveal the similarities and differences between transcriptomics- and proteomics-level perturbations in psoriasis lesions.
we particularly focused on revealing novel regulatory pathways playing a role in psoriasis development and progression.
methods
skin biopsies
acquisition of the human tissue was approved by the vavilov institute of general genetics of russian academy of sciences review board and the study was conducted after patient's consent and according to the declaration of helsinki principles.
a total of  <dig> paired nonlesional and lesional  skin biopsies from  <dig> psoriatic patients were profiled using 2d electrophoresis.
all the donors who gave biopsy tissue  provided a written informed consent for the tissue to be taken and used in this study.
clinical data for all patients are listed in additional file  <dig>
full-thickness punch biopsies were taken from uninvolved skin  and from the involved margin of a psoriatic plaque  from every patient.
sample preparation
skin biopsies from lesional  and non-lesional  skin  were mechanically homogenized on ice with mortar and pestle in solubilization buffer comprised of  <dig> m urea ,  <dig> m thiourea , 4%  chaps ,  <dig> %  triton x- <dig> ,  <dig> % ampholines 3/ <dig> ,  <dig> mm tris base ,  <dig> mm mgcl <dig> , and  <dig>  mm pmsf .
the amount of solubilization buffer taken per sample was seven times tissue wet weight.
samples were then carefully sonicated at 4°c for  <dig> sec  with bandelin sonicator at 50% power output.
following homogenization and sonication dtt  was added to reach the final concentration of  <dig> mm, and samples left at 4°c for  <dig> min.
then na2edta was added  and mixture was incubated for additional  <dig> h and then centrifuged at  <dig>  × g for  <dig> min to remove insoluble debris.
relative protein concentration was determined with gs- <dig> calibrated densitometer  after 1d-sds-page of  <dig> μl sample solution loaded on  <dig> ×  <dig> cm mini-gel stained with coomassie r- <dig> after electrophoresis.
two-dimensional electrophoresis
for the first dimension ca-ief method was used  <cit> .
glass tubes  <dig> ×  <dig>  mm  was filled with gel mixture containing  <dig>  m urea , 4%  acrylamide monomers , 2% chaps ,  <dig> %  5/ <dig> and  <dig> %  3/ <dig> ampholines , loaded with  <dig> μg protein and run  <dig> min at  <dig> v,  <dig> min at  <dig> v,  <dig> hrs at  <dig> v, and  <dig> h at  <dig> v. isoelectrofocusing was performed in protean ii xi cell .
before second dimension gels were extruded from tubes and equilibrated in  <dig> mm tris-buffer ph  <dig> , containing 2%  sds, 10%  glycerol, and 1%  dtt for  <dig> min.
for protein separation in the second dimension, 9-16% linear gradient slab gels prepared according to the standard protocol  <cit>  were used.
gel image analysis
protein spots on 2-de gels were visualized by silver staining  <cit>  and scanned with resolution  <dig> dpi.
images were analyzed using melanie iii software .
the conventional analysis involved  protein spot relative volume  determination, which was expressed as the sum of pixel intensities in the certain spot divided to the sum of pixel intensities in all spots on the gel;  gel alignment; and  spot matching.
further, sets of %vol values for every spot were processed by student test, thereby testing whether there was a significant variation of the certain protein level between two specified groups.
protein identification by maldi-tof mass-spectrometry
protein spots were cut out  from 2-de gels, destained, and in-gel digested with trypsin.
mass-spectrometry of trypsin digested proteins  was performed using a microflex maldi-tof mass-spectrometer .
peptide samples  were mixed with an equal volume of  <dig> -dihydroxybenzoic acid solution  in 20% acetonitrile and  <dig> % trifluoroacetic acid, and the resulting droplets were dried in air.
mass-spectra were obtained for mass range from  <dig> to  <dig> daltons in reflection mode and calibrated using internal standards .
peptide peak lists were formed by the snap algorithm .
proteins were identified using the mascot database search engine.
the search parameters were as follows: mass tolerance  <dig> ppm, ncbi protein sequence database, homo sapiens taxon, one missed cleavage, variable modifications by propionamide for cysteines and oxidation for methionines.
low molecular weight proteins  were identified using nanolc-ms/ms mass-spectrometry for higher convenience, regarding lack of cleavage sites in such proteins.
analysis of trypsin digest was performed on electrospray ion trap  equipped with chip-cube head.
one μl of each sample was subjected  onto reverse-phase in-chip column  for  <dig> minutes under isocratic buffer .
following sample application peptides were separated by linear gradient  of solution a  and solution b  for  <dig> minutes.
spectrum scanning was repeated three times for each sample of protein gel spots and tissue hydrolysate.
mass spectra of eluted peptide were simultaneously obtained under positive polarity for 425- <dig> m/z range both in ms and ms/ms mode,  <dig>  kv applied accumulation of  <dig> ions for  <dig> milliseconds, averages on  <dig> spectra.
mass spectra were processed with spectrum mill ms proteomics workbench software .
proteins were identified using swissprot human database in the following parameters: score ≥  <dig> for peptide and ≥  <dig> for protein, minimum s/n ratio  <dig>  maximum peptide ion charge + <dig>  precursor mass tolerance ±  <dig>  da, product mass tolerance ±  <dig>  da, proteins identification was accomplished with detection of minimum  <dig> peaks of the same peptide ion with maximum of  <dig> missed cleavages.
microarray data analysis
we used recently published data set  <cit>  from geo data base .
we compared  <dig> pairs of samples .
values for each sample were normalized by sample median value in order to unify distributions of expression signals.
for assessment of differential expression we used paired welch t-test with fdr correction  <cit> .
probe set was considered as differentially expressed if its average fold change exceeded  <dig>  and fdr corrected p-value was less than  <dig> .
overconnection analysis
all network-based analyses were conducted with metacore software suite http://www.genego.com.
this software employs a dense and manually curated database of interactions between biological objects and variety of tools for functional analysis of high-throughput data.
we defined a gene as overconnected with the gene set of interest if the corresponding node had more direct interactions with the nodes of interest than it would be expected by chance.
significance of overconnection was estimated using hypergeometric distribution with parameters r - number of interactions between examined node and the list of interest; r - degree of examined node, n -sum of interactions involving genes of interest and n - total number of interactions in the database:
hidden nodes analysis
in addition to direct interacting objects, we also used objects that may not interact directly with objects of interest but are important upstream regulators of those <cit> .
the approach is generally the same as described above, but the shortest paths instead of direct links are taken into account.
as we were interested in transcriptional regulation, we defined a transcriptional activation shortest path as the preferred shortest path from any object in the metacore database to the transcription factor target object from the data set.
we added an additional condition to include the uneven number of inhibiting interactions in the path .
if the number of such paths containing examined gene and leading to one of objects of interest were higher than expected by chance, this gene was considered as significant hidden regulator.
the significance of a node's importance was estimated using hypergeometric distribution with parameters r - number of shortest paths between containing currently examined gene; r - total number of shortest paths leading to a gene of interest through transcriptional factor, n -total number of transcription activation shortest paths containing examined gene and n - total number of transcription activation shortest paths in the database.
rank aggregation
both topology significance approaches produced lists of genes significantly linked to a gene or protein set of interest, ranked by corresponding p-values.
to combine results of these two approaches, we used a weighted rank aggregation method described in  <cit> .
weighted spearman distance was used as distance measure and the genetic algorithm was employed to select the optimal aggregated list of size  <dig>  this part of work was accomplished in r  <dig> .
<dig> http://www.r-project.org.
network analysis
in addition to topology analysis, we examined overexpressed genes and proteins using various algorithms for selecting connected biologically meaningful subnetworks enriched with objects of interest.
significance of enrichment is estimated using hypergeometric distribution.
we first used an algorithm intended to find regulatory pathways that are presumably activated under pathological conditions.
it defines a set of transcription factors that are directly regulating genes of interest and a set of receptors whose ligands are in the list of interest and then constructs series of networks; one for each receptor.
each network contains all shortest paths from a receptor to the selected transcriptional factors and their targets.
this approach allows us to reveal the most important areas of regulatory machinery affected under the investigated pathological condition.
networks are sorted by enrichment p-value.
the second applied algorithm used was aimed to define the most influential transcription factors.
it considers a transcriptional factor from the data base and gradually expands the subnetwork around it until it reaches a predefined threshold size .
networks are sorted by enrichment p-value.
RESULTS
differentially abundant proteins
protein abundance was determined by densitometric quantification of the protein spots on 2d-electophoresis gel  followed by maldi-tof mass spectrometry.
total of  <dig> proteins were over-abundant at least 2-fold in lesional skin compared with uninvolved skin: keratin  <dig>  keratin  <dig>  keratin  <dig>  squamous cell carcinoma antigen, squamous cell carcinoma antigen- <dig>  enolase  <dig>  superoxide dismutase , galectin- <dig>  s <dig> calcium-binding protein a <dig> and s <dig> calcium-binding protein a <dig> .
several of these proteins were previously reported to be over-abundant in psoriatic plaques  <cit> .
the proteins belonged to a diverse set of pathways and processes.
thus, keratin  <dig>  keratin  <dig>  and keratin  <dig> are a member of the keratin gene family.
the keratins are intermediate filament proteins responsible for the structural integrity of epithelial cells.
serpinb <dig> and serpinb <dig> are serine protease inhibitor to modulate the host immune response against tumor cells.
enolase  <dig>  more commonly known as alpha-enolase, is a glycolytic enzyme expressed in most tissues, one of the isozymes of enolase.
superoxide dismutase  <dig> protein  transforms toxic superoxide, a byproduct of the mitochondrial electron transport chain, into hydrogen peroxide and diatomic oxygen.
galectins are a family of beta-galactoside-binding proteins implicated in modulating cell-cell and cell-matrix interactions.
differential and in situ hybridizations indicate that this lectin is specifically expressed in keratinocytes.
the cellular localization and its striking down-regulation in cultured keratinocytes imply a role in cell-cell and/or cell-matrix interactions necessary for normal growth control.
s100a <dig> and s100a <dig> proteins are localized in the cytoplasm and/or nucleus of a wide range of cells, and involved in the regulation of a number of cellular processes such as cell cycle progression and differentiation.
s100a <dig> is markedly over-expressed in the skin lesions of psoriatic patients.
we attempted to connect the proteins into a network using a collection of over  <dig>  manually curated protein interactions and several variants of "shortest path" algorithms applied in metacore suite  <cit>  .
the genes encoding overabundant proteins were found to be regulated by several common transcription factors  including members of the nf-kb and ap- <dig> complexes, stat <dig>  stat <dig>  c-myc and sp <dig>  moreover, the upstream pathways activating these tfs were initiated by the overabundant s100a <dig> through its receptor rage  <cit>  and signal transduction kinases .
this network also included a positive feedback loop as s100a <dig> expression was determined to be controlled by nf-kb <cit> .
the topology of this proteomics-derived network was confirmed by several transcriptomics studies  <cit>  which showed overexpression of these tfs in psoriasis lesions.
transiently expressed tfs normally have low protein level and, therefore, usually fail to be detected by proteomics methods.
rage receptor is clearly the key regulator on this network and plays the major role in orchestrating observed changes of protein abundance.
this protein is abundant in both keratinocytes and leukocytes, though normally its expression is low  <cit> .
rage participates in a range of processes in these cell types, including inflammation.
it is being investigated as a drug target for treatment of various inflammatory disorders  <cit> .
thus, we may propose that rage can also play significant role in psoriasis.
differentially expressed genes
we used affymetrix gene expression data set from the recent study  <cit>  involving  <dig> psoriasis patients.
originally, more than  <dig> probe sets were found to be upregulated in lesions as compared with unlesional skin of the same people.
we identified  <dig> genes overexpressed in lesional skin under more stringent statistical criteria .
the list of overexpressed genes can be found in additional file  <dig>  the genes encoding  <dig> out of  <dig> proteomic markers were overexpressed, well consistent with proteomics data.
expression of enolase  <dig>  keratin  <dig> and galectin  <dig> was not altered.
common transcription regulation for overexpressed genes and differentially abundant proteins
despite good consistency between the proteomics and expression datasets, the two orders of magnitude difference in list size make direct correlation analysis difficult.
therefore, we applied interactome methods for the analysis of common upstream regulation of the two datasets at the level of transcription factors.
first, we defined the sets of the most influential transcription factors using two recently developed methods of interactome analysis  <cit>  and the "hidden nodes" algorithm  <cit> .
the former method ranks tfs based on their one-step overconnectivity with the dataset of interest compared to randomly expected number of interactions.
the latter approach takes into account direct and more distant regulation, calculating the p-values for local subnetworks by an aggregation algorithm  <cit> .
we calculated and ranked the top  <dig> tfs for each data type and added several tfs identified by network analysis approaches .
the tfs common for both data types were taken as set of 'important pathological signal transducers' .
noticeably, they closely resemble the set of tfs regulating the protein network on figure  <dig>
some genes were considered significant only by one of two topological approaches .
missing p-value means that correspondent gene has not passed  <dig>  significance threshold and has been listed among top factors only due to low p-value for other topological approach.
only one p-value was determined in the case of network analysis .
transcriptional factors common for proteomics and transcriptomics level are in bold text.
identification of influential receptors
in the next step, we applied "hidden nodes" algorithm to identify the most influential receptors that could trigger maximal possible transcriptional response.
in total, we found  <dig> membrane receptors significantly involved into regulation of  <dig> differentially expressed genes .
the complete list of receptors can be found in additional file  <dig>  assuming that topological significance alone does not necessarily prove that all receptors are involved in real signaling or are even expressed in the sample; we filtered this list by expression performance.
the receptors used were those whose encoding genes or corresponding ligands were overexpressed greater than  <dig>  fold.
we assumed that the pathways initiated by over-expressed receptors and ligands are more likely to be activated in psoriasis.
here we assumed that expression alterations and protein abundance are at least collinear.
an additional criterion was that the candidate receptors had to participate in the same signaling pathways with at least one of the common tfs.
no receptor was rejected based on this criterion.
in total,  <dig> receptors passed the transcription cut-off.
of these  <dig> receptor genes were overexpressed;  <dig> had overexpressed ligands and  <dig> cases had overexpression of both ligands and receptors .
interestingly, for several receptors, more than one ligand was over-expressed .
several receptors are composed of several subunits, only one of which was upregulated .
out of  <dig> receptors we identified by topology analysis,  <dig> were previously reported as psoriasis markers .
the other  <dig> receptors were not reported to be linked to psoriasis or known to be implicated in other inflammatory diseases.
these receptors belong to different cellular processes  .
'possible' term was used if protein name co-occurred with psoriasis in articles, but no clear evidence of its implication was shown.
in some cases, ligands are associated with psoriasis .
discussion
meta-analysis of multiple omics data types and studies is becoming an important research tool in understanding complex diseases.
several methods were developed for correlation analysis between the datasets of different type, such as mrna and proteomics  <cit> .
however, there are many technological challenges to resolve, including mismatching protein ids and mrna probes, fundamental differences in omics technologies, differences in experimental set-ups in studies done by different groups etc  <cit> .
moreover, biological reasons such as differences in rna and protein degradation processes also contribute to variability of different data types.
as a result, transcriptome and proteome datasets usually show only weak positive correlation although were considered as complimentary.
more recent studies focused on functional similarities and differences observed for different levels of cellular organization and reflected in different types of omics data  <cit> .
for example, common interacting objects were found for distinct altered transcripts and proteins in type  <dig> diabetes  <cit> .
in one leukemia study  <cit>  authors found that distinct alterations at transcriptomics and proteomic levels reflect different sides of the same deregulated cellular processes.
the overall concordance between mrna and protein expression landscapes was addressed in earlier studies, although the data types were compared mostly at the gene/protein level with limited functional analysis  <cit> .
later, ontology enrichment co-examination of transcriptomics and proteomic data has shown that the two data types affect similar biological processes and are complimentary  <cit> .
however, the key issue of biological causality and functional consequences of distinct regulation events at both mrna and protein levels of cellular organization were not yet specifically addressed.
these issues cannot be resolved by low resolution functional methods like enrichment analysis.
instead, one has to apply more precise computational methods such as topology and biological networks, which take into consideration directed binary interactions and multi-step pathways connecting objects between the datasets of different types regardless of their direct overlap at gene/protein level  <cit> .
for example, topology methods such as "hidden nodes"  <cit>  can identify and rank the upstream regulatory genes responsible for expression and protein level alterations while network tools help to uncover functional modules most affected in the datasets, identify the most influential genes/proteins within the modules and suggest how specific modules contribution to clinical phenotype  <cit> .
in this study, we observed substantial direct overlap between transcriptomics and proteomics data, as  <dig> out of  <dig> over-abundant proteins in psoriasis lesions were encoded by differentially over-expressed genes.
however, the two orders of magnitude difference in dataset size  made the standard correlation methods inapplicable.
besides, proteomics datasets display a systematic bias in function of abundant proteins, favoring "effector" proteins such as structural, inflammatory, core metabolism proteins but not the transiently expressed and fast degradable signaling proteins.
therefore, we applied topological network methods to identify common regulators for two datasets such as the most influential transcription factors and receptors.
we have identified some key regulators of the "proteomics" set among differentially expressed genes, including transcription factors, membrane receptors and extracellular ligands, thus reconstructing upstream signaling pathways in psoriasis.
in particular, we identified  <dig> receptors previously not linked to psoriasis.
importantly, many ligands and receptors defined as putative starts of signaling pathways were activated by transcription factors at the same pathways, clearly indicating on positive regulatory loops activated in psoriasis.
the versatility and the variety of signaling pathways activated in psoriasis is also impressive, which is evident from differentially overexpression of  <dig> membrane receptors and ligands in skin lesions.
this complexity and redundancy of psoriasis signaling likely contributes to the inefficiency of current treatments, even novel therapies such as monoclonal antibodies against tnf-α and il- <dig>  thus, the key regulator, rage receptor, triggers multiple signaling pathways which stay activated even when certain immunological pathways are blocked.
our study suggests that combination therapy targeting multiple pathways may be more efficient for psoriasis .
in addition, the  <dig> receptors we identified by topology analysis and previously not linked with psoriasis can be tested as potential novel targets for disease therapy.
the functional machinery of psoriasis is still not complete and additional studies can be helpful in "filling the gaps" of our understanding of its molecular mechanisms.
for instance, kinase activity is still unaccounted for, as signaling kinases are activated only transiently and are often missed in gene expression studies.
topological analysis methods such as "hidden nodes"  <cit>  may help to reconstruct regulatory events missing in the data.
also, the emerging phosphoproteomics methodology may prove to become a helpful and complimentary omics technology.
the network analysis methodology is not dependent on the type of data analyzed and or any gene/protein content overlap between the studies and is well applicable for functional integration of multiple data types.
CONCLUSIONS
we have successfully applied network-based methods to integrate and explore two distinct high-throughput disease data sets of different origin and size.
through identification of common regulatory machinery that is likely to cause overexpression of genes and proteins, we came to the signaling pathways that might contribute to the altered state of regulatory network in psoriatic lesion.
our approach allows easy integrative investigation of different data types and produces biologically meaningful results, leading to new potential therapy targets.
we have demonstrated that pathology can be caused and maintained by a great amount of various cascades, many previously not described as implicated in psoriasis; therefore, combined therapies targeting multiple pathways might be effective in treatment.
authors' contributions
ep and tn conceived the study and participated in its design and coordination; sb, sm and stm performed the proteomics analysis of psoriatic tissues, ai performed statistical and network-based analyses and wrote the manuscript, ra helped with bioinformatical analysis of the data, yn edited the manuscript.
all authors have read and approved the manuscript.
supplementary material
additional file 1
table s <dig>  patient description
click here for file
additional file 2
table s <dig>  list of genes significantly upregulated in lesions at the mrna level
click here for file
additional file 3
table s <dig>  list of receptors significantly involved in regulation of genes overexpressed in psoriatic plaque
click here for file
additional file 4
legend.
network layout legend explaining meaning of object icons and interaction types
click here for file
acknowledgements
this work was supported by grant  <dig> .
<dig>  from ministry of science and technologies of russian federation, and grant from russian academy of sciences

 psoriasis vulgaris is one of the most prevalent chronic inflammatory skin diseases affecting approximately 2% of individuals in western societies, and found worldwide in all populations. psoriasis is a complex disease affecting cellular, gene and protein levels and presented as skin lesions. the skin lesions are characterized by abnormal keratinocyte differentiation, hyperproliferation of keratinocytes, and infiltration of inflammatory cells. the mechanisms of psoriasis pathology are complex and involve genetic and environmental factors. as we gain more knowledge about molecular pathways implicated in the disease, novel therapies emerge. in recent years, microarray mrna expression profiling of lesional psoriatic skin revealed over differentially expressed genes. enrichment analysis showed that these genes encode proteins involved in regeneration, hyperkeratosis, metabolic function, immune response, and inflammation and revealed a number of modulating signaling pathways. these efforts may help to develop new-generation drugs. however, enrichment analysis limits our understanding of altered molecular interactions in psoriasis as it provides a relative ranking based on ontology terms resulting in the representation of fragmented and disconnected perturbed pathways. furthermore, analysis of gene expression alone is not sufficient for understanding the whole variety of pathological changes at different levels of cellular organization. indeed, new methodologies have been applied to the analysis of omics data in complex diseases that include algorithm-based biological network analysis and meta-analysis of multiple datasets of different types. here, we applied several techniques of network and meta-analysis to reveal the similarities and differences between transcriptomics- and proteomics-level perturbations in psoriasis lesions. we particularly focused on revealing novel regulatory pathways playing a role in psoriasis development and progression. methods skin biopsies acquisition of the human tissue was approved by the vavilov institute of general genetics of russian academy of sciences review board and the study was conducted after patient's consent and according to the declaration of helsinki principles. a total of paired nonlesional and lesional skin biopsies from psoriatic patients were profiled using 2d electrophoresis. all the donors who gave biopsy tissue provided a written informed consent for the tissue to be taken and used in this study. clinical data for all patients are listed in additional file full-thickness punch biopsies were taken from uninvolved skin and from the involved margin of a psoriatic plaque from every patient. sample preparation skin biopsies from lesional and non-lesional skin were mechanically homogenized on ice with mortar and pestle in solubilization buffer comprised of m urea, m thiourea, 4% chaps, % triton x-, % ampholines 3/, mm tris base, mm mgcl, and mm pmsf. the amount of solubilization buffer taken per sample was seven times tissue wet weight. samples were then carefully sonicated at 4°c for sec with bandelin sonicator at 50% power output. following homogenization and sonication dtt was added to reach the final concentration of mm, and samples left at 4°c for min. then na2edta was added and mixture was incubated for additional h and then centrifuged at × g for min to remove insoluble debris. relative protein concentration was determined with gs- calibrated densitometer after 1d-sds-page of μl sample solution loaded on × cm mini-gel stained with coomassie r- after electrophoresis. two-dimensional electrophoresis for the first dimension ca-ief method was used. glass tubes × mm was filled with gel mixture containing m urea, 4% acrylamide monomers, 2% chaps, % 5/ and % 3/ ampholines, loaded with μg protein and run min at v, min at v, hrs at v, and h at v. isoelectrofocusing was performed in protean ii xi cell. before second dimension gels were extruded from tubes and equilibrated in mm tris-buffer ph, containing 2% sds, 10% glycerol, and 1% dtt for min. for protein separation in the second dimension, 9-16% linear gradient slab gels prepared according to the standard protocol were used. gel image analysis protein spots on 2-de gels were visualized by silver staining and scanned with resolution dpi. images were analyzed using melanie iii software. the conventional analysis involved protein spot relative volume determination, which was expressed as the sum of pixel intensities in the certain spot divided to the sum of pixel intensities in all spots on the gel; gel alignment; and spot matching. further, sets of %vol values for every spot were processed by student test, thereby testing whether there was a significant variation of the certain protein level between two specified groups. protein identification by maldi-tof mass-spectrometry protein spots were cut out from 2-de gels, destained, and in-gel digested with trypsin. mass-spectrometry of trypsin digested proteins was performed using a microflex maldi-tof mass-spectrometer. peptide samples were mixed with an equal volume of -dihydroxybenzoic acid solution in 20% acetonitrile and % trifluoroacetic acid, and the resulting droplets were dried in air. mass-spectra were obtained for mass range from to daltons in reflection mode and calibrated using internal standards. peptide peak lists were formed by the snap algorithm. proteins were identified using the mascot database search engine. the search parameters were as follows: mass tolerance ppm, ncbi protein sequence database, homo sapiens taxon, one missed cleavage, variable modifications by propionamide for cysteines and oxidation for methionines. low molecular weight proteins were identified using nanolc-ms/ms mass-spectrometry for higher convenience, regarding lack of cleavage sites in such proteins. analysis of trypsin digest was performed on electrospray ion trap equipped with chip-cube head. one μl of each sample was subjected onto reverse-phase in-chip column for minutes under isocratic buffer. following sample application peptides were separated by linear gradient of solution a and solution b for minutes. spectrum scanning was repeated three times for each sample of protein gel spots and tissue hydrolysate. mass spectra of eluted peptide were simultaneously obtained under positive polarity for 425- m/z range both in ms and ms/ms mode, kv applied accumulation of ions for milliseconds, averages on spectra. mass spectra were processed with spectrum mill ms proteomics workbench software. proteins were identified using swissprot human database in the following parameters: score ≥ for peptide and ≥ for protein, minimum s/n ratio maximum peptide ion charge + precursor mass tolerance ± da, product mass tolerance ± da, proteins identification was accomplished with detection of minimum peaks of the same peptide ion with maximum of missed cleavages. microarray data analysis we used recently published data set from geo data base. we compared pairs of samples. values for each sample were normalized by sample median value in order to unify distributions of expression signals. for assessment of differential expression we used paired welch t-test with fdr correction. probe set was considered as differentially expressed if its average fold change exceeded and fdr corrected p-value was less than. overconnection analysis all network-based analyses were conducted with metacore software suite http://www.genego.com. this software employs a dense and manually curated database of interactions between biological objects and variety of tools for functional analysis of high-throughput data. we defined a gene as overconnected with the gene set of interest if the corresponding node had more direct interactions with the nodes of interest than it would be expected by chance. significance of overconnection was estimated using hypergeometric distribution with parameters r - number of interactions between examined node and the list of interest; r - degree of examined node, n -sum of interactions involving genes of interest and n - total number of interactions in the database: hidden nodes analysis in addition to direct interacting objects, we also used objects that may not interact directly with objects of interest but are important upstream regulators of those. the approach is generally the same as described above, but the shortest paths instead of direct links are taken into account. as we were interested in transcriptional regulation, we defined a transcriptional activation shortest path as the preferred shortest path from any object in the metacore database to the transcription factor target object from the data set. we added an additional condition to include the uneven number of inhibiting interactions in the path. if the number of such paths containing examined gene and leading to one of objects of interest were higher than expected by chance, this gene was considered as significant hidden regulator. the significance of a node's importance was estimated using hypergeometric distribution with parameters r - number of shortest paths between containing currently examined gene; r - total number of shortest paths leading to a gene of interest through transcriptional factor, n -total number of transcription activation shortest paths containing examined gene and n - total number of transcription activation shortest paths in the database. rank aggregation both topology significance approaches produced lists of genes significantly linked to a gene or protein set of interest, ranked by corresponding p-values. to combine results of these two approaches, we used a weighted rank aggregation method described in. weighted spearman distance was used as distance measure and the genetic algorithm was employed to select the optimal aggregated list of size this part of work was accomplished in r. http://www.r-project.org. network analysis in addition to topology analysis, we examined overexpressed genes and proteins using various algorithms for selecting connected biologically meaningful subnetworks enriched with objects of interest. significance of enrichment is estimated using hypergeometric distribution. we first used an algorithm intended to find regulatory pathways that are presumably activated under pathological conditions. it defines a set of transcription factors that are directly regulating genes of interest and a set of receptors whose ligands are in the list of interest and then constructs series of networks; one for each receptor. each network contains all shortest paths from a receptor to the selected transcriptional factors and their targets. this approach allows us to reveal the most important areas of regulatory machinery affected under the investigated pathological condition. networks are sorted by enrichment p-value. the second applied algorithm used was aimed to define the most influential transcription factors. it considers a transcriptional factor from the data base and gradually expands the subnetwork around it until it reaches a predefined threshold size. networks are sorted by enrichment p-value.  differentially abundant proteins protein abundance was determined by densitometric quantification of the protein spots on 2d-electophoresis gel followed by maldi-tof mass spectrometry. total of proteins were over-abundant at least 2-fold in lesional skin compared with uninvolved skin: keratin keratin keratin squamous cell carcinoma antigen, squamous cell carcinoma antigen- enolase superoxide dismutase, galectin- s calcium-binding protein a and s calcium-binding protein a. several of these proteins were previously reported to be over-abundant in psoriatic plaques. the proteins belonged to a diverse set of pathways and processes. thus, keratin keratin and keratin are a member of the keratin gene family. the keratins are intermediate filament proteins responsible for the structural integrity of epithelial cells. serpinb and serpinb are serine protease inhibitor to modulate the host immune response against tumor cells. enolase more commonly known as alpha-enolase, is a glycolytic enzyme expressed in most tissues, one of the isozymes of enolase. superoxide dismutase protein transforms toxic superoxide, a byproduct of the mitochondrial electron transport chain, into hydrogen peroxide and diatomic oxygen. galectins are a family of beta-galactoside-binding proteins implicated in modulating cell-cell and cell-matrix interactions. differential and in situ hybridizations indicate that this lectin is specifically expressed in keratinocytes. the cellular localization and its striking down-regulation in cultured keratinocytes imply a role in cell-cell and/or cell-matrix interactions necessary for normal growth control. s100a and s100a proteins are localized in the cytoplasm and/or nucleus of a wide range of cells, and involved in the regulation of a number of cellular processes such as cell cycle progression and differentiation. s100a is markedly over-expressed in the skin lesions of psoriatic patients. we attempted to connect the proteins into a network using a collection of over manually curated protein interactions and several variants of "shortest path" algorithms applied in metacore suite. the genes encoding overabundant proteins were found to be regulated by several common transcription factors including members of the nf-kb and ap- complexes, stat stat c-myc and sp moreover, the upstream pathways activating these tfs were initiated by the overabundant s100a through its receptor rage and signal transduction kinases. this network also included a positive feedback loop as s100a expression was determined to be controlled by nf-kb. the topology of this proteomics-derived network was confirmed by several transcriptomics studies which showed overexpression of these tfs in psoriasis lesions. transiently expressed tfs normally have low protein level and, therefore, usually fail to be detected by proteomics methods. rage receptor is clearly the key regulator on this network and plays the major role in orchestrating observed changes of protein abundance. this protein is abundant in both keratinocytes and leukocytes, though normally its expression is low. rage participates in a range of processes in these cell types, including inflammation. it is being investigated as a drug target for treatment of various inflammatory disorders. thus, we may propose that rage can also play significant role in psoriasis. differentially expressed genes we used affymetrix gene expression data set from the recent study involving psoriasis patients. originally, more than probe sets were found to be upregulated in lesions as compared with unlesional skin of the same people. we identified genes overexpressed in lesional skin under more stringent statistical criteria. the list of overexpressed genes can be found in additional file the genes encoding out of proteomic markers were overexpressed, well consistent with proteomics data. expression of enolase keratin and galectin was not altered. common transcription regulation for overexpressed genes and differentially abundant proteins despite good consistency between the proteomics and expression datasets, the two orders of magnitude difference in list size make direct correlation analysis difficult. therefore, we applied interactome methods for the analysis of common upstream regulation of the two datasets at the level of transcription factors. first, we defined the sets of the most influential transcription factors using two recently developed methods of interactome analysis and the "hidden nodes" algorithm. the former method ranks tfs based on their one-step overconnectivity with the dataset of interest compared to randomly expected number of interactions. the latter approach takes into account direct and more distant regulation, calculating the p-values for local subnetworks by an aggregation algorithm. we calculated and ranked the top tfs for each data type and added several tfs identified by network analysis approaches. the tfs common for both data types were taken as set of'important pathological signal transducers'. noticeably, they closely resemble the set of tfs regulating the protein network on figure some genes were considered significant only by one of two topological approaches. missing p-value means that correspondent gene has not passed significance threshold and has been listed among top factors only due to low p-value for other topological approach. only one p-value was determined in the case of network analysis. transcriptional factors common for proteomics and transcriptomics level are in bold text. identification of influential receptors in the next step, we applied "hidden nodes" algorithm to identify the most influential receptors that could trigger maximal possible transcriptional response. in total, we found membrane receptors significantly involved into regulation of differentially expressed genes. the complete list of receptors can be found in additional file assuming that topological significance alone does not necessarily prove that all receptors are involved in real signaling or are even expressed in the sample; we filtered this list by expression performance. the receptors used were those whose encoding genes or corresponding ligands were overexpressed greater than fold. we assumed that the pathways initiated by over-expressed receptors and ligands are more likely to be activated in psoriasis. here we assumed that expression alterations and protein abundance are at least collinear. an additional criterion was that the candidate receptors had to participate in the same signaling pathways with at least one of the common tfs. no receptor was rejected based on this criterion. in total, receptors passed the transcription cut-off. of these receptor genes were overexpressed; had overexpressed ligands and cases had overexpression of both ligands and receptors. interestingly, for several receptors, more than one ligand was over-expressed. several receptors are composed of several subunits, only one of which was upregulated. out of receptors we identified by topology analysis, were previously reported as psoriasis markers. the other receptors were not reported to be linked to psoriasis or known to be implicated in other inflammatory diseases. these receptors belong to different cellular processes.'possible' term was used if protein name co-occurred with psoriasis in articles, but no clear evidence of its implication was shown. in some cases, ligands are associated with psoriasis. discussion meta-analysis of multiple omics data types and studies is becoming an important research tool in understanding complex diseases. several methods were developed for correlation analysis between the datasets of different type, such as mrna and proteomics. however, there are many technological challenges to resolve, including mismatching protein ids and mrna probes, fundamental differences in omics technologies, differences in experimental set-ups in studies done by different groups etc. moreover, biological reasons such as differences in rna and protein degradation processes also contribute to variability of different data types. as a result, transcriptome and proteome datasets usually show only weak positive correlation although were considered as complimentary. more recent studies focused on functional similarities and differences observed for different levels of cellular organization and reflected in different types of omics data. for example, common interacting objects were found for distinct altered transcripts and proteins in type diabetes. in one leukemia study authors found that distinct alterations at transcriptomics and proteomic levels reflect different sides of the same deregulated cellular processes. the overall concordance between mrna and protein expression landscapes was addressed in earlier studies, although the data types were compared mostly at the gene/protein level with limited functional analysis. later, ontology enrichment co-examination of transcriptomics and proteomic data has shown that the two data types affect similar biological processes and are complimentary. however, the key issue of biological causality and functional consequences of distinct regulation events at both mrna and protein levels of cellular organization were not yet specifically addressed. these issues cannot be resolved by low resolution functional methods like enrichment analysis. instead, one has to apply more precise computational methods such as topology and biological networks, which take into consideration directed binary interactions and multi-step pathways connecting objects between the datasets of different types regardless of their direct overlap at gene/protein level. for example, topology methods such as "hidden nodes" can identify and rank the upstream regulatory genes responsible for expression and protein level alterations while network tools help to uncover functional modules most affected in the datasets, identify the most influential genes/proteins within the modules and suggest how specific modules contribution to clinical phenotype. in this study, we observed substantial direct overlap between transcriptomics and proteomics data, as out of over-abundant proteins in psoriasis lesions were encoded by differentially over-expressed genes. however, the two orders of magnitude difference in dataset size made the standard correlation methods inapplicable. besides, proteomics datasets display a systematic bias in function of abundant proteins, favoring "effector" proteins such as structural, inflammatory, core metabolism proteins but not the transiently expressed and fast degradable signaling proteins. therefore, we applied topological network methods to identify common regulators for two datasets such as the most influential transcription factors and receptors. we have identified some key regulators of the "proteomics" set among differentially expressed genes, including transcription factors, membrane receptors and extracellular ligands, thus reconstructing upstream signaling pathways in psoriasis. in particular, we identified receptors previously not linked to psoriasis. importantly, many ligands and receptors defined as putative starts of signaling pathways were activated by transcription factors at the same pathways, clearly indicating on positive regulatory loops activated in psoriasis. the versatility and the variety of signaling pathways activated in psoriasis is also impressive, which is evident from differentially overexpression of membrane receptors and ligands in skin lesions. this complexity and redundancy of psoriasis signaling likely contributes to the inefficiency of current treatments, even novel therapies such as monoclonal antibodies against tnf-α and il- thus, the key regulator, rage receptor, triggers multiple signaling pathways which stay activated even when certain immunological pathways are blocked. our study suggests that combination therapy targeting multiple pathways may be more efficient for psoriasis. in addition, the receptors we identified by topology analysis and previously not linked with psoriasis can be tested as potential novel targets for disease therapy. the functional machinery of psoriasis is still not complete and additional studies can be helpful in "filling the gaps" of our understanding of its molecular mechanisms. for instance, kinase activity is still unaccounted for, as signaling kinases are activated only transiently and are often missed in gene expression studies. topological analysis methods such as "hidden nodes" may help to reconstruct regulatory events missing in the data. also, the emerging phosphoproteomics methodology may prove to become a helpful and complimentary omics technology. the network analysis methodology is not dependent on the type of data analyzed and or any gene/protein content overlap between the studies and is well applicable for functional integration of multiple data types.  we have successfully applied network-based methods to integrate and explore two distinct high-throughput disease data sets of different origin and size. through identification of common regulatory machinery that is likely to cause overexpression of genes and proteins, we came to the signaling pathways that might contribute to the altered state of regulatory network in psoriatic lesion. our approach allows easy integrative investigation of different data types and produces biologically meaningful results, leading to new potential therapy targets. we have demonstrated that pathology can be caused and maintained by a great amount of various cascades, many previously not described as implicated in psoriasis; therefore, combined therapies targeting multiple pathways might be effective in treatment. authors' contributions ep and tn conceived the study and participated in its design and coordination; sb, sm and stm performed the proteomics analysis of psoriatic tissues, ai performed statistical and network-based analyses and wrote the manuscript, ra helped with bioinformatical analysis of the data, yn edited the manuscript. all authors have read and approved the manuscript. supplementary material additional file 1 table s patient description click here for file additional file 2 table s list of genes significantly upregulated in lesions at the mrna level click here for file additional file 3 table s list of receptors significantly involved in regulation of genes overexpressed in psoriatic plaque click here for file additional file 4 legend. network layout legend explaining meaning of object icons and interaction types click here for file acknowledgements this work was supported by grant. from ministry of science and technologies of russian federation, and grant from russian academy of sciences 
several cell types are perturbed in this pathology, and underlying signaling events are complex and still poorly understood. investigation of proteomics and transcriptomics data sets on psoriasis revealed versatility in regulatory machinery underlying pathology and showed complementarities between two levels of cellular organization. we developed network-based methodology for integrative analysis of high throughput data sets of different types. we also evaluated functional synergy between transcriptomics and proteomics results. network-based analysis revealed similarities in regulation at both proteomics and transcriptomics level. we identified a group of transcription factors responsible for overexpression of psoriasis genes and a number of previously unknown signaling pathways that may play a role in this process.
 psoriasis is complex inflammatory skin pathology of autoimmune origin. several cell types are perturbed in this pathology, and underlying signaling events are complex and still poorly understood.  in order to gain insight into molecular machinery underlying the disease, we conducted a comprehensive meta-analysis of proteomics and transcriptomics of psoriatic lesions from independent studies. network-based analysis revealed similarities in regulation at both proteomics and transcriptomics level. we identified a group of transcription factors responsible for overexpression of psoriasis genes and a number of previously unknown signaling pathways that may play a role in this process. we also evaluated functional synergy between transcriptomics and proteomics results.  we developed network-based methodology for integrative analysis of high throughput data sets of different types. investigation of proteomics and transcriptomics data sets on psoriasis revealed versatility in regulatory machinery underlying pathology and showed complementarities between two levels of cellular organization. 