CHARACTERIZATION OF SURVIVAL ASSOCIATED GENE INTERACTIONS AND LYMPHOCYTE HETEROGENEITY IN CANCER

Cancer is the second leading cause of death globally. Tumors form intricate ecosystems in which malignant and immune cells interact to shape disease progression. Yet, the molecular underpinnings of tumorigenesis and immunological responses to tumors are poorly understood, limiting their manipulation to elicit favorable clinical outcomes. This thesis lays conceptual frameworks for investigating the molecular interactions taking place in tumors as well as the diversity of the immune response to cancer. In the molecular level of individual cancer cells, the phenotypic effect of perturbing a gene’s activity depends on the activity level of other genes, reflecting the notion that phenotypes are emergent properties of a network of functionally interacting genes. In the context of cancer, contemporary investigations have primarily focused on just one type of functional genetic interaction (GI) – synthetic lethality (SL). However, there may be additional types of GIs whose systematic identification would enrich the molecular and functional characterization of cancer. This thesis describes a novel data-driven approach called EnGIne, that applied to large-scale cancer data identifies 71,946 GIs spanning 12 distinct types, only a small minority of which are SLs. The detected GIs explain cancer driver genes’ tissue- specificity and differences in patients’ response to drugs, and stratify breast cancer tumors into refined subtypes. These results expand the scope of cancer GIs and lay a conceptual and computational basis for future studies of additional types of GIs and their translational applications. Furthermore, tumor growth is continuously shaped by the immune response. However, T cells typically adopt a dysfunctional phenotype may be reversed using immunotherapy strategies. Most current tumor immunotherapies leverage cytotoxic CD8+ T cells to elicit an effective anti-tumor response. Despite evidence for clinical potential of CD4+ tumor-infiltrating lymphocytes (TILs), their functional diversity has limited our ability to harness their anti-tumor activity. To address this issue, we have used single-cell mRNA sequencing (scRNAseq) to analyze the response of CD4+ T cells specific for a defined recombinant tumor antigen, both in the tumor microenvironment and draining lymph nodes (dLN). New computational approaches to characterize subpopulations identified TIL transcriptomic patterns strikingly distinct from those elicited by responses to infection, and dominated by diversity among T-bet-expressing T helper type 1 (Th1)-like cells. In contrast, the dLN response includes Follicular helper (Tfh)-like cells but lacks Th1 cells. We identify an interferon-driven signature in Th1-like TILs, and show that it is found in human liver cancer and melanoma, in which it is negatively associated with response to checkpoint therapy. Our study unveils unsuspected differences between tumor and virus CD4+ T cell responses, and provides a proof-of-concept methodology to characterize tumor- control CD4+ T cell effector programs. Targeting these programs should help improve immunotherapy strategies

Single cell dissection of plasma cell heterogeneity in symptomatic and asymptomatic myeloma

Multiple myeloma, a plasma cell malignancy, is the second most common blood cancer. Despite extensive research, disease heterogeneity is poorly characterized, hampering efforts for early diagnosis and improved treatments. Here, we apply single cell RNA sequencing to study the heterogeneity of 40 individuals along the multiple myeloma progression spectrum, including 11 healthy controls, demonstrating high interindividual variability that can be explained by expression of known multiple myeloma drivers and additional putative factors. We identify extensive subclonal structures for 10 of 29 individuals with multiple myeloma. In asymptomatic individuals with early disease and in those with minimal residual disease post-treatment, we detect rare tumor plasma cells with molecular characteristics similar to those of active myeloma, with possible implications for personalized therapies. Single cell analysis of rare circulating tumor cells allows for accurate liquid biopsy and detection of malignant plasma cells, which reflect bone marrow disease. Our work establishes single cell RNA sequencing for dissecting blood malignancies and devising detailed molecular characterization of tumor cells in symptomatic and asymptomatic patients

asmagen/SPAGEfinder 0.0.1

Computational approach to identify Survival associated Pairwise Gene Expression states

A network diffusion approach to inferring sample-specific function reveals functional changes associated with breast cancer

<div><p><i>Guilt-by-association</i> codifies the empirical observation that a gene’s function is informed by its neighborhood in a biological network. This would imply that when a gene’s network context is altered, for instance in disease condition, so could be the gene’s function. Although context-specific changes in biological networks have been explored, the potential changes they may induce on the functional roles of genes are yet to be characterized. Here we analyze, for the first time, the network-induced potential functional changes in breast cancer. Using transcriptomic samples for 1047 breast tumors and 110 healthy breast tissues from TCGA, we derive sample-specific protein interaction networks and assign sample-specific functions to genes via a diffusion strategy. Testing for significant changes in the inferred functions between normal and cancer samples, we find several functions to have significantly gained or lost genes in cancer, not due to differential expression of genes known to perform the function, but rather due to changes in the network topology. Our predicted functional changes are supported by mutational and copy number profiles in breast cancers. Our diffusion-based functional assignment provides a novel characterization of a tumor that is complementary to the standard approach based on functional annotation alone. Importantly, this characterization is effective in predicting patient survival, as well as in predicting several known histopathological subtypes of breast cancer.</p></div

Contingency table.

<p>The following table is generated to determine if elevated missense (respectively nonsense) mutation frequencies are enriched among functions with net gain (respectively net loss).</p

Functions ranked based on functional variability of genes.

<p>Top 10 gained (green) and lost (red) functions are shown, along with <i>Δ</i>_<i>f</i>, <i>Δ</i>_<i>f</i> divided (normalized) by the number of genes annotated by the function, followed by the sample shuffling, and the log fold change, which is the log ratio of the average number of expressed genes annotated by <i>f</i> in cancer and normal samples.</p

Diffusion based functional heterogeneity across clinical subtypes.

<p>The following figure displays the log ratio between the average numbers of genes assigned to each function by diffusion (represented by columns) across samples annotated with a subtype (represented by rows) versus the rest of the samples.</p

Overall approach.

<p>The reference gene is depicted by black circle. The initial static global PIN is projected onto normal and cancer samples based on gene expression, and each function (red and green) are diffused through each PIN. In this case, the reference gene is assigned green function in normal and red function in cancer, i.e., the gene gained red and lost the green function in cancer.</p