Sunil M. Kurian

BACKGROUND: Genomic 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. RESULTS: We 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. CONCLUSIONS: The 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.

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 clinicians. Lastly, we conducted bioinformatic analyses to identify biological pathways, mechanisms and medication targets. Overall, suicidality may be underlined, at least in part, by biological mechanisms related to stress, inflammation and apoptosis.