BackgroundGenomic and other high dimensional analyses often require one to summarize multiple related variables by a single representative. This task is also variously referred to as collapsing, combining, reducing, or aggregating variables. Examples include summarizing several probe measurements corresponding to a single gene, representing the expression profiles of a co-expression module by a single expression profile, and aggregating cell-type marker information to de-convolute expression data. Several standard statistical summary techniques can be used, but network methods also provide useful alternative methods to find representatives. Currently few collapsing functions are developed and widely applied.ResultsWe introduce the R function collapseRows that implements several collapsing methods and evaluate its performance in three applications. First, we study a crucial step of the meta-analysis of microarray data: the merging of independent gene expression data sets, which may have been measured on different platforms. Toward this end, we collapse multiple microarray probes for a single gene and then merge the data by gene identifier. We find that choosing the probe with the highest average expression leads to best between-study consistency. Second, we study methods for summarizing the gene expression profiles of a co-expression module. Several gene co-expression network analysis applications show that the optimal collapsing strategy depends on the analysis goal. Third, we study aggregating the information of cell type marker genes when the aim is to predict the abundance of cell types in a tissue sample based on gene expression data ("expression deconvolution"). We apply different collapsing methods to predict cell type abundances in peripheral human blood and in mixtures of blood cell lines. Interestingly, the most accurate prediction method involves choosing the most highly connected "hub" marker gene. Finally, to facilitate biological interpretation of collapsed gene lists, we introduce the function userListEnrichment, which assesses the enrichment of gene lists for known brain and blood cell type markers, and for other published biological pathways.ConclusionsThe R function collapseRows implements several standard and network-based collapsing methods. In various genomic applications we provide evidence that both types of methods are robust and biologically relevant tools.

A major challenge for kidney transplantation is balancing the need for immunosuppression to prevent rejection, while minimizing drug-induced toxicities.We used DNA microarrays (HG-U95Av2 GeneChips, Affymetrix) to determine gene expression profiles for kidney biopsies and peripheral blood lymphocytes (PBLs) in transplant patients including normal donor kidneys, well-functioning transplants without rejection, kidneys undergoing acute rejection, and transplants with renal dysfunction without rejection. We developed a data analysis schema based on expression signal determination, class comparison and prediction, hierarchical clustering, statistical power analysis and real-time quantitative PCR validation. We identified distinct gene expression signatures for both biopsies and PBLs that correlated significantly with each of the different classes of transplant patients. This is the most complete report to date using commercial arrays to identify unique expression signatures in transplant biopsies distinguishing acute rejection, acute dysfunction without rejection and well-functioning transplants with no rejection history. We demonstrate for the first time the successful application of high density DNA chip analysis of PBL as a diagnostic tool for transplantation. The significance of these results, if validated in a multicenter prospective trial, would be the establishment of a metric based on gene expression signatures for monitoring the immune status and immunosuppression of transplanted patients.

Our aim was to determine outcomes with transplanting kidneys from deceased donors with acute kidney injury, defined as a donor with terminal serum creatinine ≥2.0 mg/dL, or a donor requiring acute renal replacement therapy. We included all patients who received deceased donor kidney transplant from June 2004 to October 2013. There were 162 AKI donor transplant recipients (21% of deceased donor transplants): 139 in the standard criteria donor (SCD) and 23 in the expanded criteria donor (ECD) cohort. 71% of the AKI donors had stage 3 (severe AKI), based on acute kidney injury network (AKIN) staging. Protocol biopsies were done at 1, 4, and 12 months posttransplant. One and four month formalin‐fixed paraffin embedded (FFPE) biopsies from 48 patients (24 AKI donors, 24 non‐AKI) underwent global gene expression profiling using DNA microarrays (96 arrays). DGF was more common in the AKI group but eGFR, graft survival at 1 year and proportion with IF/TA>2 at 1 year were similar for the two groups. At 1 month, there were 898 differentially expressed genes in the AKI group (p‐value <0.005; FDR <10%), but by 4 months there were no differences. Transplanting selected kidneys from deceased donors with AKI is safe and has excellent outcomes.

We performed a randomized prospective trial comparing calcineurin inhibitor (CNI)-free to CNI-based imished prevalence of CAN and down-regulated expression of genes responsible for progression of CAN. All may provide for an alternative natural history with improved graft survival.

Suicides are a leading cause of death in psychiatric patients, and in society at large. Developing more quantitative and objective ways (biomarkers) for predicting and tracking suicidal states would have immediate practical applications and positive societal implications. We undertook such an endeavor. First, building on our previous blood biomarker work in mood disorders and psychosis, we decided to identify blood gene expression biomarkers for suicidality, looking at differential expression of genes in the blood of subjects with a major mood disorder (bipolar disorder), a high-risk population prone to suicidality. We compared no suicidal ideation (SI) states and high SI states using a powerful intrasubject design, as well as an intersubject case–case design, to generate a list of differentially expressed genes. Second, we used a comprehensive Convergent Functional Genomics (CFG) approach to identify and prioritize from the list of differentially expressed gene biomarkers of relevance to suicidality. CFG integrates multiple independent lines of evidence—genetic and functional genomic data—as a Bayesian strategy for identifying and prioritizing findings, reducing the false-positives and false-negatives inherent in each individual approach. Third, we examined whether expression levels of the blood biomarkers identified by us in the live bipolar subject cohort are actually altered in the blood in an age-matched cohort of suicide completers collected from the coroner's office, and report that 13 out of the 41 top CFG scoring biomarkers (32%) show step-wise significant change from no SI to high SI states, and then to the suicide completers group. Six out of them (15%) remained significant after strict Bonferroni correction for multiple comparisons. Fourth, we show that the blood levels of SAT1 (spermidine/spermine N1–acetyltransferase 1), the top biomarker identified by us, at the time of testing for this study, differentiated future as well as past hospitalizations with suicidality, in a live cohort of bipolar disorder subjects, and exhibited a similar but weaker pattern in a live cohort of psychosis (schizophrenia/schizoaffective disorder) subjects. Three other (phosphatase and tensin homolog (PTEN), myristoylated alanine-rich protein kinase C substrate (MARCKS), and mitogen-activated protein kinase kinase kinase 3 (MAP3K3)) of the six biomarkers that survived Bonferroni correction showed similar but weaker effects. Taken together, the prospective and retrospective hospitalization data suggests SAT1, PTEN, MARCKS and MAP3K3 might be not only state biomarkers but trait biomarkers as well. Fifth, we show how a multi-dimensional approach using SAT1 blood expression levels and two simple visual-analog scales for anxiety and mood enhances predictions of future hospitalizations for suicidality in the bipolar cohort (receiver-operating characteristic curve with area under the curve of 0.813). Of note, this simple approach does not directly ask about SI, which some individuals may deny or choose not to share with clinician...