Combined meta-analysis of systemic effects of allogeneic stem cell transplantation and systemic sclerosis
© Grigoryev et al.; licensee BioMed Central Ltd. 2014
Received: 8 July 2013
Accepted: 6 March 2014
Published: 22 March 2014
Chronic graft-versus-host disease (cGVHD) is a major factor of morbidity and mortality for allogeneic stem cell transplantation (aSCT). The skin and internal organ involvement is the most common systemic complication of cGVHD and closely resembles systemic sclerosis (SSc). Circulating lymphocytes characterize the autoimmune nature of both conditions. Therefore we hypothesized that the common clinical manifestation (systemic organ and skin injury) and the common underlying players (lymphocytes) justify the combined meta-analysis of these diseases.
The aSCT and SSc datasets were uploaded from Gene Expression Omnibus (GEO), a public functional genomics data repository. The available microarray studies of peripheral blood mononuclear cells (PBMCs) and isolated lymphocytes were limited to well established microarray platforms (Affymetrix, Agilent, Canvac, and Illumina) and experimental settings with ≥10 patients per group. The resulting pools of data were merged by unique gene identifier and analyzed by the expression genome-wide association studies (eGWAS) coupled with the subtraction of the cGVHD+ and cGVHD− molecular signatures. The eGWAS was applied to 47 and 50 lymphocyte profiles from aSCT and SSc patients, respectively. The identified 35 candidates were represented by 8 known cGVHD genes (including CXCR4, LTBR and PML) and 28 new candidate genes (including SEPX1 and DNJGB1). The further mutual subtraction of cGVHD+ and cGVHD− candidates and pathway analysis identified a list of 25 genes. Seven of these genes belong to the fibroblast development and function pathway, consisting of the well known cGVHD genes CCND1, JUN, and FOS, and the new molecular targets MMP2, FOSB, TNFAIP8, and DUSP1. These genes become primary candidates for a potential link of systemic effects of cGVHD and SSc.
We designed a new approach for meta-analysis by combining data from different diseases using common clinical manifestation as a linker. This allowed us to power up the insufficient standalone meta-analysis of aSCT microarray studies, by adding SSc samples to the data pool. This new method has successfully identified novel molecular targets for systemic effects of both aSCT and SSc. We believe that this approach is generalizable and can be applied to an array of diseases with common clinical manifestations.
KeywordsMeta-analysis Allogeneic stem cell transplantation Systemic sclerosis Microarray Molecular signature Gene expression Public data repository Peripheral mononuclear cells Circulating lymphocytes
Allogeneic stem cell transplantation (aSCT) has become a curative therapy for increasing numbers of diseases. To date, it is the only successful cellular immunotherapy for high-risk malignancies such as leukemia. In pediatrics, it also offers curative therapy for nonmalignant blood disorders such as thalassemia, immune dysregulation, congenital bone marrow failure syndromes, inborn errors of metabolism and autoimmune conditions . However, 30-70% of aSCT cases develop cGVHD, which occurs 100 days after transplantation with median time to onset of 4–6 months . The precise mechanisms of cGVHD are unknown and evaluation of global transcriptional changes of reconstituted immune cells from aSCT patients is a feasible approach for studying cGVHD . For instance, the transcript profiling of CD4+ and CD8+ lymphocytes from donors for aSCT demonstrated that predictive biomarkers for cGVHD can be detected .
Due to the paucity of microarray data for standalone meta-analysis of cGVHD, we used similarities in systemic responses between cGVHD and systemic sclerosis (SSc) to link their microarray datasets [5, 6]. The similarity in systemic responses between cGVHD and SSc was demonstrated in the experimental settings and was attributed to the excessive activation of T and B cells in both conditions .Therefore, we hypothesized that the combined meta-analysis of both diseases would lead to better understanding of systemic effects of circulating lymphocytes in aSCT and SSc.
In this study we combined publicly available datasets for aSCT (with or without cGVHD) and SSc (diffuse and limited) and analyzed them using a meta-analysis technique.
Results and discussion
Eight genes among these candidates were previously linked to cGVHD (Figure 1), which demonstrated that despite the amplification of cGVHD data with unrelated but similar clinical conditions, we generated a feasible candidate list for cGVHD. The top candidate was lymphotoxin beta receptor or TNFR superfamily, member 3 (LTBR) with P = 1.73 × 10−13 (Figure 1). Concordantly, it was recently reported that the expression of LTBR is elevated in patients with cGVHD+ as compared to cGVHD−, and that the application of lymphotoxin beta receptor-immunoglobulin fusion protein in mouse model of cGVHD is beneficial [8, 9].
Given that eGWAS confirmed the feasibility of combined meta-analysis of two different diseases with common clinical manifestations, we proceeded with the mutual subtraction of molecular signatures of cGVHD and SSc.
The comparison of molecular changes of CD4+ and CD8+ T-lymphocytes 1 year after aSCT was conducted by chi-square test and demonstrated a significant difference (P < 0.0001) in ratios of up- and down-regulated genes between these cell types. In CD4+ T-cells the molecular signature was shifted towards suppression of gene expression (75 up-regulated genes vs 135 down-regulated genes), while in CD8+ T-cells the molecular signature was shifted towards activation (95 up-regulated genes vs 51 down-regulated genes). The overall transcriptional changes were more pronounced in CD4+ cells (210 genes) versus CD8+ cells (146 genes). These findings suggested that CD4+ plays a bigger role in aSCT than CD8+.
The significance of co-expression of these 25 genes was evaluated by Student’s T-test. The contribution of CD4+ lymphocytes to both cGVHD and SSc was 2.55 times higher (P = 0.033) than the contribution of CD8+ lymphocytes, which confirmed our observation of more pronounced involvement of CD4+. We next tested which type of SSc the molecular signature of cGVHD lymphocytes most resembles. While we were able to demonstrate that the molecular signature of cGVHD lymphocytes resembles diffuse SSc by 2.2 fold more than limited SSc, the t-test of this comparison did not reach significance (p = 0.124).
The greater contribution of CD4+ to cGVHD molecular signature was further confirmed by the pathway analysis. We evaluated the relevance of cGVHD transcripts to known biological processes using the Ingenuity Pathway Analysis (IPA) tool (http://www.ingenuity.com). 163 out of 210 transcripts of CD4+ lymphocytes and 127 out of 146 transcripts of CD8+ lymphocytes were mapped by the software to known genes, while the unmapped transcripts consisted of hypothetical proteins and open reading frames.
These genes were unreported by the original studies, most likely due to different goals and focuses of the authors. This is an example of how the re-analysis of publicly available data using new approaches will bring to light molecular targets unsuspected by the original studies.
Gene candidates for skin manifestation of cGVHD
Stem cell transplantation
Graft versus host disease
Lymphotoxin beta receptor
FBJ murine osteosarcoma viral oncogene homolog
CDC-like kinase 1
Matrix metallopeptidase 2 (gelatinase A)
DEAD (Asp-Glu-Ala-Asp) box polypeptide 3, X-linked
GRB2-associated binding protein 1
FBJ murine osteosarcoma viral oncogene homolog B
Tumor necrosis factor, alpha-induced protein 8
SATB homeobox 1
Dual specificity phosphatase 1
Ubiquitin specific peptidase 18
Dnaj (Hsp40) homolog, subfamily B, member 1
Elongation factor, RNA polymerase II, 2
Serum/glucocorticoid regulated kinase
Retinoid X receptor, alpha
Zinc finger protein 331
Protein kinase, AMP-activated, alpha 1 catalytic subunit
Methionine sulfoxide reductase B1
ARP6 actin-related protein 6 homolog (yeast)
Ubiquitin specific peptidase 47
Of the remaining genes considered new candidates for cGVHD, MMP2, FOSB, TNFAIP8, and DUSP1 were from the fibroblast proliferation pathway and all, but TNFAIP8, were assigned to the same gene-gene interaction network (Figure 6). Moreover, amongst these genes, DUSP1 was detected by eGWAS studies as well (Figure 1), therefore DUSP1 became our primary candidate.
The expression of DUSP1 in lymphocytes had been detected by many groups (18 citations), however, this gene has not yet been implicated in aSCT or SSc (Table 1). Despite this, the recent studies of circulating cells in peripheral blood of kidney transplant patients identified DUSP1 as a potential biomarker for monitoring renal graft status . We believe that results of our meta-analysis and the reported allograft sensitivity of DUSP1 warrants further studies of this gene in the cGVHD setting.
The other new candidate, MMP2, had the strongest link to aSCT (30 citations) and was detected in scleroderma (2 citations). However it has not yet been linked to cGVHD (Table 1 and Figure 5). In addition to its ability to degrade extracellular matrix proteins, MMP2 can also act on several nonmatrix proteins such as calcitonin gene-related peptide, thus promoting vasoconstriction. Moreover, its C-terminal non-catalytic fragment possesses an anti-angiogenic property (http://www.uniprot.org). One can speculate that the combination of vasoconstriction and anti-angiogenisity can be contributory to hypoxia and systemic sclerosis. Furthermore, the potential immunogenetic role of MMP2 was observed in a mouse model of heart transplant. It was demonstrated that alloreactivity of T-cells was significantly lower in MMP2-deficient mice compared to their wild type littermates . All these findings make MMP2 an extremely appealing target for cGVHD.
Our final candidate, FOSB, is less known in aSCT (1 citation), but was reported in supplemental materials of original studies submitted to GEO . The gene-gene interaction analysis put this gene into the central cross-linker nodule of two major nodes of cGVHD candidates (Figure 6). The main function of FOSB is believed to enhance binding activity of JUN proteins to DNA. Furthermore, the finding that during induced clonal anergy of CD4+ T-cells, JUN and FOS proteins, but not FOSB, reduce their DNA-binding ability  makes FOSB an interesting candidate that may prolong the active state of T lymphocytes under suppressive conditions, thus promoting the aSCT complications.
Presented here data indicate that: 1) for the first time in the field of meta-analysis of microarray data we conducted cross-disease analysis based on common clinical manifestations; 2) CD4+ T-cells undergo more pronounced changes during cGVHD than their CD8+ counterparts and these changes are shifted towards down-regulation; 3) the molecular signature of cGVHD does not demonstrate preferable similarity with either type of systemic sclerosis; 4) the pathway analysis linked expressional changes in lymphocytes from cGVHD patients to fibroblast proliferation; 5) we identified 25 gene candidates common to both cGVHD and SSc, of which primary candidates MMP2, FOSB, and DUSP1are the most appealing for further studies. We believe that our approach of in silico amplification of sparse microarray data by linking molecular signatures of different diseases with the same systemic complications and conducting meta-analysis on such combined dataset is generalizable and capable of detecting new molecular targets.
The aSCT and SSc datasets were uploaded from gene expression omnibus (GEO, NCBI, http://www.ncbi.nlm.nih.gov/geo). The term “bone marrow transplant” was submitted to the database query windows and this search returned 81 entries (as of December 14, 2012), 23 of which were data series, while the remaining data were either annotation of microarray platforms or individual samples.
The inclusion criteria for retrieved data were the array samples, which must: 1) represent genome wide studies of human peripheral blood mononuclear cells (PBMCs) or isolated lymphocytes; 2) be generated with established microarray platforms and the number of interrogated sequences should be >5000, thus excluding custom platforms that are pathway oriented and will introduce bias towards given biological process; and 3) the experimental settings must have ≥10 patients per condition.
Summary of microarray data used for meta-analysis
Baron et al .
aSCT + cGVHD
Baron et al. 
Cheadle et al. 
Risbano et al. 
Pendergrass et al. 
Baron et al. : GSM103566-67, GSM103610-11, GSM103616-17, GSM103642-43, GSM103680-81, GSM103684-85, GSM103700-01, GSM103712-13, GSM103716-17, GSM103720-21, GSM103724-25, and GSM103728-29. Their CD4+/CD8+ cells were from 12 donors, whose cells did not cause cGVHD in a recipient 1 year after aSCT; GSM103626-27, GSM103630-31, GSM103634-35, GSM103638-39, GSM103650-51, GSM103668-69, GSM103678-79, GSM103686-87, GSM103689, GSM103696-97, GSM103568-69, GSM103570-71, GSM103573, GSM103584-85, GSM103592-93, GSM103594-95, GSM103604-05, GSM103704-05, and GSM103707. Their CD4+/CD8+ cells were from 18 donors, whose cells did cause cGVHD in a recipient 1 year after aSCT (one CD8+ sample was not reported, see Table 2); GSM103620-21, GSM103640-41, GSM103710-11, GSM103714-15, GSM103718-19, GSM103722-23, GSM103726-27, GSM103652-53, GSM103682-83, and GSM103698-99. Their CD4+/CD8+ cells were from 10 aSCT patients, who remain cGVHD-free 1 year after aSCT; GSM103622-23, GSM103624-25, GSM103628-29, GSM103632-33, GSM103636-37, GSM103572, GSM103582-83, GSM103590-91, GSM103618-19, GSM103648-49, GSM103666-67, GSM103688, GSM103694-95, GSM103702-03, and GSM103706. Their CD4+/CD8+ cells were from 14 cGVHD patients, who developed cGVHD 1 year after aSCT (one CD8+ sample is not reported, see Table 2).
Cheadle et al. : GSM827665-705 – PBMCs were from 41 healthy controls, and GSM827778-96, PBMCs from 19 patients with diffuse SSc.
Risbano et al. : GSM556441-50, PBMCs were from 10 healthy controls, and GSM556413-22, PBMCs from 10 patients with diffuse SSc.
Pendergrass et al.: GSM489236-41, GSM489243, GSM489245-47, PBMCs were from 10 healthy controls, and GSM489194-98, GSM489200-01, GSM489204, GSM489207-218, PBMCs from 21 patients with limited SSc.
Linking different microarray platforms
Obtained GSMs were cross-linked based on the GEO provided Minimum Information About a Microarray Experiment (MIAME) data as described previously [16, 17]. Briefly, the probe IDs across different microarray platforms were linked using the Array Information Library Universal Navigator (AILUN) tool (http://ailun.stanford.edu) and the probes that remained unmatched by AILUN were added to the main AILUN created dataset by matching HUGO-approved gene symbols.
Gene candidate selection
To estimate differences between groups of samples from aSCT and aSCT + cGVHD patients, and healthy and SSc patients, raw post-quantitation microarray data was median trimmed using minimal trim of 1 highest and 1 lowest expressors. Given the general trend of microarray data series to have low number of biological replicates, the more vigorous trimming was not applied, due to the great reduction of statistical power. The probes without a positive signal throughout all the experiments were then removed, and the remaining data were reanalyzed using SAM 2.0 software , applying default settings without application of arbitrary restrictions , as described previously . For each gene in every microarray experiment, the d score, which denotes the standardized change in gene expression of diseased lymphocytes versus corresponding healthy lymphocytes, was calculated. Given that the total of 32671 probes were analyzed, the Bonferroni threshold (P = 1.53 × 10−6) was used for detecting significance.
For eGWAS analysis, 32671 probes from all experiments were merged into one dataset and significant and not significant changes counted and  evaluated by chi-square test as described previously . Genes with an absolute value of d score ≥2  or an absolute value of fold change ≥1.2  between healthy and diseased groups were considered significantly dysregulated.
Molecular signature subtraction
The gene lists for subtraction analysis were generated using “ranking order” approach. SAM output was ranked according to the fold change and genes from the top quartile were considered significantly affected by a corresponding disease. The subtraction of transcript lists was conducted using VENNY tool (http://bioinfogp.cnb.csic.es/tools/venny/index.html). The significance of association of lymphocyte types (CD4+ or CD8+) in cGVHD samples with the PBMC samples from scleroderma types (diffuse or limited) was assessed using two tailed T-test. Given that scleroderma is represented in the GEO collection by PBMC data only, but knowing that PBMCs are comprised of up to 60% of CD4+ and up to 30% of CD8+, we considered it acceptable to directly compare CD4+ and CD8+ fractions of cGVHD samples with PBMC samples of SSc patients. Values submitted to the T-test were first normalized to the total number of lymphocyte type or scleroderma type, respectively. Association between molecular signatures with P value <0.05 was considered significant.
Automated literature search
PubMatrix (multiplex literature mining tool) analysis  was conducted as described previously . We restricted our search to human symbols approved by HUGO Gene Nomenclature Committee (HGNC), which were enriched by all aliases and former (discontinued) symbols for selected candidate genes (http://www.genenames.org). The symbols with the hyphen and double names separated by a slash were filtered out. The biomedical literature reference count for a given gene was represented by the highest number among HUGO symbol or its aliases. The genes with number of references N ≥ 5 were considered to have established association with cGVHD, 5 > N > 0 were considered novel (have been noticed during cGVHD studies, but were not linked to the disease), and N = 0 were considered new.
The pathway analysis, which links the most relevant biological processes to a provided list of candidate genes, was conducted using the Ingenuity Pathways Knowledge Base tool (Ingenuity Systems, Inc., Redwood City, CA.) as described previously [24, 25]. Gene symbols from CD4+ and CD8+ significant gene lists were submitted to the IPA and analyzed using “Core” tool. Significance of the identified pathways was tested by the IPA imbedded Fisher Exact test and expressed as p-value . Pathways, that were built of 5 or more genes and had p < 0.05 were considered significant.
Array Information Library Universal Navigator
Allogeneic stem cell transplantation
Chronic graft-versus-host disease
Expression genome-wide association studies
Gene expression omnibus
GEO sample IDs or gene expression series matrix
HUGO Gene Nomenclature Committee
Ingenuity pathways analysis
Lymphotoxin beta receptor or TNFR superfamily, member 3
Minimum information about a microarray experiment
National Center for Biotechnology Information
Peripheral blood mononuclear cells
We thank Dr. Daniel Heruth for his critical reading and editing of our manuscript. This work is in part supported by the start-up fund and William R. Brown/Missouri Endowment of Children’s Mercy Hospitals and Clinics, University of Missouri at Kansas City (Ye, S.Q.).
- Rozmus J, Schultz KR: Biomarkers in chronic graft-versus-host disease. Expert Rev Hematol. 2011, 4 (3): 329-342. 10.1586/ehm.11.27.View ArticlePubMedPubMed CentralGoogle Scholar
- Lee SJ, Flowers ME: Recognizing and managing chronic graft-versus-host disease. Hematology Am Soc Hematol Educ Program. 2008, 1: 134-141.Google Scholar
- Buzzeo MP, Yang J, Casella G, Reddy V: A preliminary gene expression profile of acute graft-versus-host disease. Cell Transplant. 2008, 17 (5): 489-494. 10.3727/096368908785096042.View ArticlePubMedGoogle Scholar
- Baron C, Somogyi R, Greller LD, Rineau V, Wilkinson P, Cho CR, Cameron MJ, Kelvin DJ, Chagnon P, Roy DC, Busque L, Sekaly RP, Perreault C: Prediction of graft-versus-host disease in humans by donor gene-expression profiling. PLoS medicine. 2007, 4 (1): e23-10.1371/journal.pmed.0040023.View ArticlePubMedPubMed CentralGoogle Scholar
- Kitko CL, White ES, Baird K: Fibrotic and sclerotic manifestations of chronic graft-versus-host disease. Biol Blood Marrow Transplant. 2012, 18 (1 Suppl): S46-52.View ArticlePubMedPubMed CentralGoogle Scholar
- White JM, Creamer D, du Vivier AW, Pagliuca A, Ho AY, Devereux S, Salisbury JR, Mufti GJ: Sclerodermatous graft-versus-host disease: clinical spectrum and therapeutic challenges. Br J Dermatol. 2007, 156 (5): 1032-1038. 10.1111/j.1365-2133.2007.07827.x.View ArticlePubMedGoogle Scholar
- Zhang C, Todorov I, Zhang Z, Liu Y, Kandeel F, Forman S, Strober S, Zeng D: Donor CD4+ T and B cells in transplants induce chronic graft-versus-host disease with autoimmune manifestations. Blood. 2006, 107 (7): 2993-3001. 10.1182/blood-2005-09-3623.View ArticlePubMedGoogle Scholar
- Poloni A, Sartini D, Emanuelli M, Trappolini S, Mancini S, Pozzi V, Costantini B, Serrani F, Berardinelli E, Renzi E, Olivieri A, Leoni P: Gene expression profile of cytokines in patients with chronic graft-versus-host disease after allogeneic hematopoietic stem cell transplantation with reduced conditioning. Cytokine. 2011, 53 (3): 376-383. 10.1016/j.cyto.2010.12.008.View ArticlePubMedGoogle Scholar
- Srinivasan M, Flynn R, Price A, Ranger A, Browning JL, Taylor PA, Ritz J, Antin JH, Murphy WJ, Luznik L, Shlomchik MJ, Panoskaltsis-Mortari A, Blazar BR: Donor B-cell alloantibody deposition and germinal center formation are required for the development of murine chronic GVHD and bronchiolitis obliterans. Blood. 2012, 119 (6): 1570-1580. 10.1182/blood-2011-07-364414.View ArticlePubMedPubMed CentralGoogle Scholar
- Li L, Khatri P, Sigdel TK, Tran T, Ying L, Vitalone MJ, Chen A, Hsieh S, Dai H, Zhang M, Naesens M, Zarkhin V, Sansanwal P, Chen R, Mindrinos M, Xiao W, Benfield M, Ettenger RB, Dharnidharka V, Mathias R, Portale A, McDonald R, Harmon W, Kershaw D, Vehaskari VM, Kamil E, Baluarte HJ, Warady B, Davis R, Butte AJ, Salvatierra O, Sarwal MM: A peripheral blood diagnostic test for acute rejection in renal transplantation. Am J Transplant. 2012, 12 (10): 2710-2718. 10.1111/j.1600-6143.2012.04253.x.View ArticlePubMedPubMed CentralGoogle Scholar
- Campbell LG, Ramachandran S, Liu W, Shipley JM, Itohara S, Rogers JG, Moazami N, Senior RM, Jaramillo A: Different roles for matrix metalloproteinase-2 and matrix metalloproteinase-9 in the pathogenesis of cardiac allograft rejection. Am J Transplant. 2005, 5 (3): 517-528. 10.1111/j.1600-6143.2005.00744.x.View ArticlePubMedGoogle Scholar
- Heisel O, Keown P: Alterations in transcription factor binding at the IL-2 promoter region in anergized human CD4+ T lymphocytes. Transplantation. 2001, 72 (8): 1416-1422. 10.1097/00007890-200110270-00015.View ArticlePubMedGoogle Scholar
- Cheadle C, Berger AE, Mathai SC, Grigoryev DN, Watkins TN, Sugawara Y, Barkataki S, Fan J, Boorgula M, Hummers L, Zaiman AL, Girgis R, McDevitt MA, Johns RA, Wigley F, Barnes KC, Hassoun PM: Erythroid-specific transcriptional changes in PBMCs from pulmonary hypertension patients. PloS One. 2012, 7 (4): e34951-10.1371/journal.pone.0034951.View ArticlePubMedPubMed CentralGoogle Scholar
- Risbano MG, Meadows CA, Coldren CD, Jenkins TJ, Edwards MG, Collier D, Huber W, Mack DG, Fontenot AP, Geraci MW, Bull TM: Altered immune phenotype in peripheral blood cells of patients with scleroderma-associated pulmonary hypertension. Clin Transl Sci. 2010, 3 (5): 210-218. 10.1111/j.1752-8062.2010.00218.x.View ArticlePubMedPubMed CentralGoogle Scholar
- Pendergrass SA, Hayes E, Farina G, Lemaire R, Farber HW, Whitfield ML, Lafyatis R: Limited systemic sclerosis patients with pulmonary arterial hypertension show biomarkers of inflammation and vascular injury. PloS One. 2010, 5 (8): e12106-10.1371/journal.pone.0012106.View ArticlePubMedPubMed CentralGoogle Scholar
- Chen R, Li L, Butte AJ: AILUN: reannotating gene expression data automatically. Nat Methods. 2007, 4 (11): 879-10.1038/nmeth1107-879.View ArticlePubMedPubMed CentralGoogle Scholar
- Grigoryev DN, Cheranova DI, Heruth DP, Huang P, Zhang LQ, Rabb H, Ye SQ: Meta-analysis of molecular response of kidney to ischemia reperfusion injury for the identification of new candidate genes. BMC Nephrol. 2013, 14: 231-10.1186/1471-2369-14-231.View ArticlePubMedPubMed CentralGoogle Scholar
- Tusher VG, Tibshirani R, Chu G: Significance analysis of microarrays applied to the ionizing radiation response. Proc Natl Acad Sci U S A. 2001, 98 (9): 5116-5121. 10.1073/pnas.091062498.View ArticlePubMedPubMed CentralGoogle Scholar
- Larsson O, Wahlestedt C, Timmons JA: Considerations when using the significance analysis of microarrays (SAM) algorithm. BMC Bioinform. 2005, 6: 129-10.1186/1471-2105-6-129.View ArticleGoogle Scholar
- Grigoryev DN, Mathai SC, Fisher MR, Girgis RE, Zaiman AL, Housten-Harris T, Cheadle C, Gao L, Hummers LK, Champion HC, Garcia JG, Wigley FM, Tuder RM, Barnes KC, Hassoun PM: Identification of candidate genes in scleroderma-related pulmonary arterial hypertension. Transl Res. 2008, 151 (4): 197-207.View ArticlePubMedPubMed CentralGoogle Scholar
- Kodama K, Horikoshi M, Toda K, Yamada S, Hara K, Irie J, Sirota M, Morgan AA, Chen R, Ohtsu H, Maeda S, Kadowaki T, Butte AJ: Expression-based genome-wide association study links the receptor CD44 in adipose tissue with type 2 diabetes. Proc Natl Acad Sci U S A. 2012, 109 (18): 7049-7054. 10.1073/pnas.1114513109.View ArticlePubMedPubMed CentralGoogle Scholar
- Grigoryev DN, Ma SF, Irizarry RA, Ye SQ, Quackenbush J, Garcia JG: Orthologous gene-expression profiling in multi-species models: search for candidate genes. Genome Biol. 2004, 5 (5): R34-10.1186/gb-2004-5-5-r34.View ArticlePubMedPubMed CentralGoogle Scholar
- Becker KG, Hosack DA, Dennis G, Lempicki RA, Bright TJ, Cheadle C, Engel J: PubMatrix: a tool for multiplex literature mining. BMC Bioinform. 2003, 4: 61-10.1186/1471-2105-4-61.View ArticleGoogle Scholar
- Grigoryev DN, Liu M, Hassoun HT, Cheadle C, Barnes KC, Rabb H: The local and systemic inflammatory transcriptome after acute kidney injury. J Am Soc Nephrol. 2008, 19 (3): 547-558. 10.1681/ASN.2007040469.View ArticlePubMedPubMed CentralGoogle Scholar
- Hassoun HT, Lie ML, Grigoryev DN, Liu M, Tuder RM, Rabb H: Kidney ischemia-reperfusion injury induces caspase-dependent pulmonary apoptosis. Am J Physiol Renal Physiol. 2009, 297 (1): F125-137. 10.1152/ajprenal.90666.2008.View ArticlePubMedPubMed CentralGoogle Scholar
- Jimenez-Marin A, Collado-Romero M, Ramirez-Boo M, Arce C, Garrido JJ: Biological pathway analysis by ArrayUnlock and Ingenuity Pathway Analysis. BMC Proc. 2009, 3 (Suppl 4): S6-10.1186/1753-6561-3-s4-s6.View ArticlePubMedPubMed CentralGoogle Scholar
- The pre-publication history for this paper can be accessed here:http://www.biomedcentral.com/2052-1839/14/7/prepub
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly credited.