Ragas: integration and enhanced visualization for single cell subcluster analysis.

SUMMARY: Subcluster analysis is a powerful means to improve clustering and characterization of single cell RNA-Seq data. However, there are no existing tools to systematically integrate results from multiple subclusters, which creates hurdles for accurate data quantification, visualization, and interpretation in downstream analysis. To address this issue, we developed Ragas, an R package that integrates multi-level subclustering objects for streamlined analysis and visualization. A new data structure was implemented to seamlessly connect and assemble miscellaneous single cell analyses from different levels of subclustering, along with several new or enhanced visualization functions. Moreover, a re-projection algorithm was developed to integrate nearest-neighbor graphs from multiple subclusters in order to maximize their separability on the combined cell embeddings, which significantly improved the presentation of rare and homogeneous subpopulations. AVAILABILITY AND IMPLEMENTATION: The Ragas package and its documentation can be accessed through https://github.com/jig4003/Ragas and its source code is also available at https://zenodo.org/records/11244921.

Human Rhinovirus Infection of the Respiratory Tract Affects Sphingolipid Synthesis.

The 17q21 asthma susceptibility locus includes asthma risk alleles associated with decreased sphingolipid synthesis, likely resulting from increased expression of ORMDL3. ORMDL3 inhibits serine-palmitoyl transferase (SPT), the rate-limiting enzyme of de novo sphingolipid synthesis. There is evidence that decreased sphingolipid synthesis is critical to asthma pathogenesis. Children with asthma and 17q21 asthma risk alleles display decreased sphingolipid synthesis in blood cells. Reduced SPT activity results in airway hyperreactivity, a hallmark feature of asthma. 17q21 asthma risk alleles are also linked to childhood infections with human rhinovirus (RV). This study evaluates the interaction of RV with the de novo sphingolipid synthesis pathway, and the alterative effects of concurrent SPT inhibition in SPT-deficient mice and human airway epithelial cells. In mice, RV infection shifted lung sphingolipid synthesis gene expression to a pattern that resembles genetic SPT deficiency, including decreased expression of Sptssa, a small SPT subunit. This pattern was pronounced in lung epithelial cellular adhesion molecule (EpCAM+) cells and reproduced in human bronchial epithelial cells. RV did not affect Sptssa expression in lung CD45+ immune cells. RV increased sphingolipids unique to the de novo synthesis pathway in mouse lung and human airway epithelial cells. Interestingly, these de novo sphingolipid species were reduced in the blood of RV-infected wild-type mice. RV exacerbated SPT deficiency-associated airway hyperreactivity. Airway inflammation was similar in RV-infected wild-type and SPT-deficient mice. This study reveals the effects of RV infection on the de novo sphingolipid synthesis pathway, elucidating a potential mechanistic link between 17q21 asthma risk alleles and rhinoviral infection.

A human Barrett's esophagus organoid system reveals epithelial-mesenchymal plasticity induced by acid and bile salts.

The pathogenesis of subsquamous intestinal metaplasia (SSIM), in which glands of Barrett's esophagus (BE) are buried under esophageal squamous epithelium, is unknown. In a rat model of reflux esophagitis, we found that columnar-lined esophagus developed via a wound-healing process involving epithelial-mesenchymal plasticity (EMP) that buried glands under ulcerated squamous epithelium. To explore a role for reflux-induced EMP in BE, we established and characterized human Barrett's organoids and sought evidence of EMP after treatment with acidic bile salts (AB). We optimized media to grow human BE organoids from immortalized human Barrett's cells and from BE biopsies from seven patients, and we characterized histological, morphological, and molecular features of organoid development. Features and markers of EMP were explored following organoid exposure to AB, with and without a collagen I (COL1) matrix to simulate a wound-healing environment. All media successfully initiated organoid growth, but advanced DMEM/F12 (aDMEM) was best at sustaining organoid viability. Using aDMEM, organoids comprising nongoblet and goblet columnar cells that expressed gastric and intestinal cell markers were generated from BE biopsies of all seven patients. After AB treatment, early-stage Barrett's organoids exhibited EMP with loss of membranous E-cadherin and increased protrusive cell migration, events significantly enhanced by COL1. Using human BE biopsies, we have established Barrett's organoids that recapitulate key histological and molecular features of BE to serve as high-fidelity BE models. Our findings suggest that reflux can induce EMP in human BE, potentially enabling Barrett's cells to migrate under adjacent squamous epithelium to form SSIM.NEW & NOTEWORTHY Using Barrett's esophagus (BE) biopsies, we established organoids recapitulating key BE features. During early stages of organoid development, a GERD-like wound environment-induced features of epithelial-mesenchymal plasticity (EMP) in Barrett's progenitor cells, suggesting that reflux-induced EMP can enable Barrett's cells to migrate underneath squamous epithelium to form subsquamous intestinal metaplasia, a condition that may underlie Barrett's cancers that escape detection by endoscopic surveillance, and recurrences of Barrett's metaplasia following endoscopic eradication therapy.

In Esophageal Squamous Cells From Eosinophilic Esophagitis Patients, Th2 Cytokines Increase Eotaxin-3 Secretion Through Effects on Intracellular Calcium and a Non-Gastric Proton Pump.

BACKGROUND & AIMS: In upper airway cells, T helper 2 cytokines that signal through interleukin-4 (IL-4) receptor-α have been shown to stimulate eotaxin-3 secretion via a nongastric proton pump (ngH+,K+ATPase). To seek novel targets for eosinophilic esophagitis (EoE) treatments, we evaluated ngH+,K+ATPase expression in EoE squamous cells, and explored molecular pathways involved in eotaxin-3 secretion by IL-4 receptor-α signaling. METHODS: ngH+,K+ATPase expression in EoE cells was evaluated by quantitative real-time polymerase chain reaction and Western blotting. IL-4-stimulated eotaxin-3 secretion was measured by enzyme-linked immunosorbent assay after treatment with omeprazole, SCH 28080 (potassium-competitive acid blocker), ethylene glycol-bis(ß-aminoethyl)-N,N,N',N'-tetraacetoxymethyl ester (calcium chelator), 2-aminoethoxydiphenyl borate (inhibitor of endoplasmic reticulum calcium release), verapamil, and diltiazem (L-type calcium channel inhibitors). Intracellular calcium transients were measured by Fluo-4 fluorescence. Key experiments were confirmed in EoE primary cells and in RNA sequencing datasets from mucosal biopsies of patients with EoE and controls. RESULTS: EoE cells expressed ngH+,K+ATPase messenger RNA and protein. Omeprazole and SCH 28080 decreased IL-4-stimulated eotaxin-3 secretion. IL-4 increased intracellular calcium transients, and IL-4-stimulated eotaxin-3 secretion was blocked by ethylene glycol-bis(ß-aminoethyl)-N,N,N',N'-tetraacetoxymethyl ester, 2-aminoethoxydiphenyl borate, verapamil, and diltiazem. The combination of omeprazole and verapamil suppressed IL-4-stimulated eotaxin-3 secretion more than either agent alone. EoE biopsies expressed higher ngH+,K+ATPase and exhibited more calcium signaling than controls. CONCLUSIONS: EoE cells express a nongastric proton pump that mediates T helper 2 cytokine-stimulated eotaxin-3 secretion. IL-4 induces calcium release from the endoplasmic reticulum and calcium entry via L-type calcium channels, increasing intracellular calcium that contributes to eotaxin-3 secretion by EoE cells. L-type calcium channel inhibitors block T helper 2 cytokine-stimulated eotaxin-3 secretion, suggesting a potential role for these agents in EoE treatment.

Sensitive and specific post-call filtering of genetic variants in xenograft and primary tumors.

Motivation: Tumor genome sequencing offers great promise for guiding research and therapy, but spurious variant calls can arise from multiple sources. Mouse contamination can generate many spurious calls when sequencing patient-derived xenografts. Paralogous genome sequences can also generate spurious calls when sequencing any tumor. We developed a BLAST-based algorithm, Mouse And Paralog EXterminator (MAPEX), to identify and filter out spurious calls from both these sources. Results: When calling variants from xenografts, MAPEX has similar sensitivity and specificity to more complex algorithms. When applied to any tumor, MAPEX also automatically flags calls that potentially arise from paralogous sequences. Our implementation, mapexr, runs quickly and easily on a desktop computer. MAPEX is thus a useful addition to almost any pipeline for calling genetic variants in tumors. Availability and implementation: The mapexr package for R is available at https://github.com/bmannakee/mapexr under the MIT license. Contact: mannakee@email.arizona.edu or rgutenk@email.arizona.edu or eknudsen@email.arizona.edu. Supplementary information: Supplementary data are available at Bioinformatics online.

Pancreatic cancer cell lines as patient-derived avatars: genetic characterisation and functional utility.

OBJECTIVE: Pancreatic ductal adenocarcinoma (PDAC) is a therapy recalcitrant disease with the worst survival rate of common solid tumours. Preclinical models that accurately reflect the genetic and biological diversity of PDAC will be important for delineating features of tumour biology and therapeutic vulnerabilities. DESIGN: 27 primary PDAC tumours were employed for genetic analysis and development of tumour models. Tumour tissue was used for derivation of xenografts and cell lines. Exome sequencing was performed on the originating tumour and developed models. RNA sequencing, histological and functional analyses were employed to determine the relationship of the patient-derived models to clinical presentation of PDAC. RESULTS: The cohort employed captured the genetic diversity of PDAC. From most cases, both cell lines and xenograft models were developed. Exome sequencing confirmed preservation of the primary tumour mutations in developed cell lines, which remained stable with extended passaging. The level of genetic conservation in the cell lines was comparable to that observed with patient-derived xenograft (PDX) models. Unlike historically established PDAC cancer cell lines, patient-derived models recapitulated the histological architecture of the primary tumour and exhibited metastatic spread similar to that observed clinically. Detailed genetic analyses of tumours and derived models revealed features of ex vivo evolution and the clonal architecture of PDAC. Functional analysis was used to elucidate therapeutic vulnerabilities of relevance to treatment of PDAC. CONCLUSIONS: These data illustrate that with the appropriate methods it is possible to develop cell lines that maintain genetic features of PDAC. Such models serve as important substrates for analysing the significance of genetic variants and create a unique biorepository of annotated cell lines and xenografts that were established simultaneously from same primary tumour. These models can be used to infer genetic and empirically determined therapeutic sensitivities that would be germane to the patient.

Mechanistic insights into anticancer properties of oligomeric proanthocyanidins from grape seeds in colorectal cancer.

Although the anticancer properties of oligomeric proanthocyanidins (OPCs) from grape seeds have been well recognized, the molecular mechanisms by which they exert anticancer effects are poorly understood. In this study, through comprehensive RNA-sequencing-based gene expression profiling in multiple colorectal cancer cell lines, we for the first time illuminate the genome-wide effects of OPCs from grape seeds in colorectal cancer. Our data revealed that OPCs affect several key cancer-associated genes. In particular, genes involved in cell cycle and DNA replication were most significantly and consistently altered by OPCs across multiple cell lines. Intriguingly, our in vivo experiments showed that OPCs were significantly more potent at decreasing xenograft tumor growth compared with the unfractionated grape seed extract (GSE) that includes the larger polymers of proanthocyanidins. These findings were further confirmed in colorectal cancer patient-derived organoids, wherein OPCs more potently inhibited the formation of organoids compared with GSE. Furthermore, we validated alteration of cell cycle and DNA replication-associated genes in cancer cell lines, mice xenografts as well as patient-derived organoids. Overall, this study provides an unbiased and comprehensive look at the mechanisms by which OPCs exert anticancer properties in colorectal cancer.

Astrocyte elevated gene-1 and c-Myc cooperate to promote hepatocarcinogenesis in mice.

UNLABELLED: Astrocyte elevated gene-1 (AEG-1) and c-Myc are overexpressed in human hepatocellular carcinoma (HCC) functioning as oncogenes. AEG-1 is transcriptionally regulated by c-Myc, and AEG-1 itself induces c-Myc by activating the Wnt/ß-catenin-signaling pathway. We now document the cooperation of AEG-1 and c-Myc in promoting hepatocarcinogenesis by analyzing hepatocyte-specific transgenic mice expressing either AEG-1 (albumin [Alb]/AEG-1), c-Myc (Alb/c-Myc), or both (Alb/AEG-1/c-Myc). Wild-type and Alb/AEG-1 mice did not develop spontaneous HCC. Alb/c-Myc mice developed spontaneous HCC without distant metastasis, whereas Alb/AEG-1/c-Myc mice developed highly aggressive HCC with frank metastasis to the lungs. Induction of carcinogenesis by N-nitrosodiethylamine significantly accelerated the kinetics of tumor formation in all groups. However, in Alb/AEG-1/c-Myc, the effect was markedly pronounced with lung metastasis. In vitro analysis showed that Alb/AEG-1/c-Myc hepatocytes acquired increased proliferation and transformative potential with sustained activation of prosurvival and epithelial-mesenchymal transition-signaling pathways. RNA-sequencing analysis identified a unique gene signature in livers of Alb/AEG-1/c-Myc mice that was not observed when either AEG-1 or c-Myc was overexpressed. Specifically, Alb/AEG-1/c-Myc mice overexpressed maternally imprinted noncoding RNAs (ncRNAs), such as Rian, Meg-3, and Mirg, which are implicated in hepatocarcinogenesis. Knocking down these ncRNAs significantly inhibited proliferation and invasion by Alb/AEG-1/c-Myc hepatocytes. CONCLUSION: Our studies reveal a novel cooperative oncogenic effect of AEG-1 and c-Myc that might explain the mechanism of aggressive HCC. Alb/AEG-1/c-Myc mice provide a useful model to understand the molecular mechanism of cooperation between these two oncogenes and other molecules involved in hepatocarcinogenesis. This model might also be of use for evaluating novel therapeutic strategies targeting HCC.

Retinoblastoma protein potentiates the innate immune response in hepatocytes: significance for hepatocellular carcinoma.

UNLABELLED: Cancers mediated by viral etiology must exhibit deregulated cellular proliferation and evade immune recognition. The role of the retinoblastoma tumor suppressor (RB) pathway, which is lost at relatively high frequency in hepatocellular carcinoma (HCC), has recently been expanded to include the regulation of innate immune responsiveness. In this study we investigated the coordinate impact of RB-loss on cell cycle control and immune function in the liver. We found that RB depletion in hepatoma cells resulted in a compromised immunological response to multiple stimuli and reduced the potential of these cells to recruit myeloid cells. Viral-mediated liver-specific RB deletion in vivo led to the induction of genes associated with proliferation and cell cycle entry as well as the significant attenuation of genes associated with immune function, as evidenced by decreases in cytokine and chemokine expression, leukocyte recruitment, and hepatic inflammation. To determine if these changes in gene expression were instructive in human disease, we compared our liver-specific RB-loss gene signature to existing profiles of HCC and found that this signature was associated with disease progression and confers a worse prognosis. CONCLUSION: Our data confirm that RB participates in the regulation of innate immunity in liver parenchymal cells both in vitro and in vivo and to our knowledge describes the first gene signature associated with HCC that includes both immunoregulatory and proliferative genes and that can also be attributed to the alteration of a single gene in vitro.

Khdc3 Regulates Metabolism Across Generations in a DNA-Independent Manner.

Genetic variants can alter the profile of heritable molecules such as small RNAs in sperm and oocytes, and in this manner ancestral genetic variants can have a significant effect on offspring phenotypes even if they are not themselves inherited. Here we show that wild type female mice descended from ancestors with a mutation in the mammalian germ cell gene Khdc3 have hepatic metabolic defects that persist over multiple generations. We find that genetically wild type females descended from Khdc3 mutants have transcriptional dysregulation of critical hepatic metabolic genes, which persist over multiple generations and pass through both female and male lineages. This was associated with dysregulation of hepatically-metabolized molecules in the blood of these wild type mice with mutational ancestry. The oocytes of Khdc3-null females, as well as their wild type descendants, had dysregulation of multiple small RNAs, suggesting that these epigenetic changes in the gametes transmit the phenotype between generations. Our results demonstrate that ancestral mutation in Khdc3 can produce transgenerational inherited phenotypes, potentially indefinitely.

An unconventional mechanism of IL-1ß secretion that requires Type I IFN in lupus monocytes.

Systemic Lupus Erythematosus (SLE) is characterized by autoreactive B cell activation, upregulation of Type I Interferon (IFN) and widespread inflammation. Mitochondrial nucleic acids (NAs) are increasingly recognized as triggers of IFN 1 . Thus, defective removal of mitochondria from mature red blood cells (Mito + RBCs), a feature of SLE, contributes to IFN production by myeloid cells 2 . Here we identify blood monocytes (Mo) that have internalized RBCs and co-express IFN-stimulated genes (ISGs) and interleukin-1ß (IL-1ß) in SLE patients with active disease. We show that ISG expression requires the interaction between Mito + RBC-derived mitochondrial DNA (mtDNA) and cGAS, while IL-1ß production entails Mito + RBC-derived mitochondrial RNA (mtRNA) triggering of RIG-I-like receptors (RLRs). This leads to the cytosolic release of Mo-derived mtDNA that activates the NLRP3 inflammasome. Importantly, IL-1ß release depends on the IFN-inducible myxovirus resistant protein 1 (MxA), which enables the translocation of this cytokine into a trans-Golgi network (TGN)-mediated unconventional secretory pathway. Our study highlights a novel and synergistic pathway involving IFN and the NLRP3 inflammasome in SLE.

Dietary long-chain omega 3 fatty acids modify sphingolipid metabolism to facilitate airway hyperreactivity.

Omega-3 polyunsaturated fatty acids (n-3 PUFAs) are essential nutrients that can affect inflammatory responses. While n-3 PUFAs are generally considered beneficial for cardiovascular disease and obesity, the effects on asthma, the most common inflammatory lung disease are unclear. While prenatal dietary n-3 PUFAs decrease the risk for childhood wheezing, postnatal dietary n-3 PUFAs can worsen allergic airway inflammation. Sphingolipid metabolism is also affected by dietary n-3 PUFAs. Decreased sphingolipid synthesis leads to airway hyperreactivity, besides inflammation, a cardinal feature of asthma, and common genetic asthma risk alleles lead to lower sphingolipid synthesis. We investigated the effect of dietary n-3 PUFAs on sphingolipid metabolism and airway reactivity. Comparing a fish-oil diet with a high n-3 PUFA content (FO) to an isocaloric coconut oil-enriched diet (CO), we found an n-3 PUFA-dependent effect on increased airway reactivity, that was not accompanied by inflammation. Lung and whole blood content of dihydroceramides, ceramides, sphingomyelins, and glucosylceramides were lower in mice fed the n-3 PUFA enriched diet consistent with lower sphingolipid synthesis. In contrast, phosphorylated long chain bases such as sphingosine 1-phosphate were increased. These findings suggest that dietary n-3 PUFAs affect pulmonary sphingolipid composition to favor innate airway hyperreactivity, independent of inflammation, and point to an important role of n-3 PUFAs in sphingolipid metabolism.

Dysregulated DNA methylation in the pathogenesis of Fabry disease.

Fabry disease is an X-linked lysosomal storage disorder caused by a deficiency of α-galactosidase A and subsequent accumulation of glycosphingolipids with terminal α-D-galactosyl residues. The molecular process through which this abnormal metabolism of glycosphingolipids causes multisystem dysfunction in Fabry disease is not fully understood. We sought to determine whether dysregulated DNA methylation plays a role in the development of this disease. In the present study, using isogenic cellular models derived from Fabry patient endothelial cells, we tested whether manipulation of α-galactosidase A activity and glycosphingolipid metabolism affects DNA methylation. Bisulfite pyrosequencing revealed that changes in α-galactosidase A activity were associated with significantly altered DNA methylation in the androgen receptor promoter, and this effect was highly CpG loci-specific. Methylation array studies showed that α-galactosidase A activity and glycosphingolipid levels were associated with differential methylation of numerous CpG sites throughout the genome. We identified 15 signaling pathways that may be susceptible to methylation alterations in Fabry disease. By incorporating RNA sequencing data, we identified 21 genes that have both differential mRNA expression and methylation. Upregulated expression of collagen type IV alpha 1 and alpha 2 genes correlated with decreased methylation of these two genes. Methionine levels were elevated in Fabry patient cells and Fabry mouse tissues, suggesting that a perturbed methionine cycle contributes to the observed dysregulated methylation patterns. In conclusion, this study provides evidence that α-galactosidase A deficiency and glycosphingolipid storage may affect DNA methylation homeostasis and highlights the importance of epigenetics in the pathogenesis of Fabry disease and, possibly, of other lysosomal storage disorders.

Screening for SARS-CoV-2 via PCR and serological testing in asymptomatic healthcare workers.

The prevalence and seroconversion rate of SARS-CoV-2 infection among asymptomatic health care workers in the US is unclear. Our study utilized real-time polymerase chain reaction (RT-PCR) SARS-CoV-2 testing and serological evaluation to detect IgG antibodies specific to SARS-CoV-2 antigens in asymptomatic health care workers. A total of 197 subjects with a mean age of 35 years were recruited into the study. While most (67%) reported prolonged contact with known COVID-19 patients, only 8 (4.2%) tested positive on RT-PCR and 23 (11.7%) had detectable levels of IgG antibody to SARS-CoV-2. Out of 19 subjects with detectable IgG antibody at week 1, 11 (57.9%) lost their antibody response by week 3. No statistically significant difference was found in baseline characteristics or exposure status between subjects with positive and negative results on RT-PCR or antibody positivity. In conclusion, we found a low incidence of PCR positivity for SARS-CoV-2 in a high-risk group. This likely demonstrates the effectiveness of proper personal protective equipment use and low transmission risk in health care settings. The detectable IgG antibody titer was low, and a significant portion of subjects lost their antibody response on repeat testing. This may mean that antibody response in asymptomatic patients is categorically different than in symptomatic hospitalized patients with COVID-19.

Cerina: systematic circRNA functional annotation based on integrative analysis of ceRNA interactions.

Circular RNAs, a family of covalently circularized RNAs with tissue-specific expression, were recently demonstrated to play important roles in mammalian biology. Regardless of extensive research to predict, quantify, and annotate circRNAs, our understanding of their functions is still in its infancy. In this study, we developed a novel computational tool: Competing Endogenous RNA for INtegrative Annotations (Cerina), to predict biological functions of circRNAs based on the competing endogenous RNA model. Pareto Frontier Analysis was employed to integrate ENCODE mRNA/miRNA data with predicted microRNA response elements to prioritize tissue-specific ceRNA interactions. Using data from several circRNA-disease databases, we demonstrated that Cerina significantly improved the functional relevance of the prioritized ceRNA interactions by several folds, in terms of precision and recall. Proof-of-concept studies on human cancers and cardiovascular diseases further showcased the efficacy of Cerina on predicting potential circRNA functions in human diseases.

A combination of curcumin and oligomeric proanthocyanidins offer superior anti-tumorigenic properties in colorectal cancer.

Combining anti-cancer agents in cancer therapies is becoming increasingly popular due to improved efficacy, reduced toxicity and decreased emergence of resistance. Here, we test the hypothesis that dietary agents such as oligomeric proanthocyanidins (OPCs) and curcumin cooperatively modulate cancer-associated cellular mechanisms to inhibit carcinogenesis. By a series of in vitro assays in colorectal cancer cell lines, we showed that the anti-tumorigenic properties of the OPCs-curcumin combination were superior to the effects of individual compounds. By RNA-sequencing based gene-expression profiling in six colorectal cancer cell lines, we identified the cooperative modulation of key cancer-associated pathways such as DNA replication and cell cycle pathways. Moreover, several pathways, including protein export, glutathione metabolism and porphyrin metabolism were more effectively modulated by the combination of OPCs and curcumin. We validated genes belonging to these pathways, such as HSPA5, SEC61B, G6PD, HMOX1 and PDE3B to be cooperatively modulated by the OPCs-curcumin combination. We further confirmed that the OPCs-curcumin combination more potently suppresses colorectal carcinogenesis and modulated expression of genes identified by RNA-sequencing in mice xenografts and in colorectal cancer patient-derived organoids. Overall, by delineating the cooperative mechanisms of action of OPCs and curcumin, we make a case for the clinical co-administration of curcumin and OPCs as a treatment therapy for patients with colorectal cancer.

Stratification of Pancreatic Ductal Adenocarcinoma: Combinatorial Genetic, Stromal, and Immunologic Markers.

Purpose: Pancreatic ductal adenocarcinoma (PDAC) is associated with an immunosuppressive milieu that supports immune system evasion and disease progression. Here, we interrogated genetic, stromal, and immunologic features of PDAC to delineate impact on prognosis and means to more effectively employ immunotherapy.Experimental Design: A cohort of 109 PDAC cases annotated for overall survival was utilized as a primary discovery cohort. Gene expression analysis defined immunologic subtypes of PDAC that were confirmed in the Cancer Genome Atlas dataset. Stromal and metabolic characteristics of PDAC cases were evaluated by histologic analysis and immunostaining. Enumeration of lymphocytes, as well as staining for CD8, FOXP3, CD68, CD163, PDL1, and CTLA4 characterized immune infiltrate. Neoantigens were determined by analysis of whole-exome sequencing data. Random-forest clustering was employed to define multimarker subtypes, with univariate and multivariate analyses interrogating prognostic significance.Results: PDAC cases exhibited distinct stromal phenotypes that were associated with prognosis, glycolytic and hypoxic biomarkers, and immune infiltrate composition. Immune infiltrate was diverse among PDAC cases and enrichment for M2 macrophages and select immune checkpoints regulators were specifically associated with survival. Composite analysis with neoantigen burden, immunologic, and stromal features defined novel subtypes of PDAC that could have bearing on sensitivity to immunologic therapy approaches. In addition, a subtype with low levels of neoantigens and minimal lymphocyte infiltrate was associated with improved overall survival.Conclusions: The mutational burden of PDAC is associated with distinct immunosuppressive mechanisms that are conditioned by the tumor stromal environment. The defined subtypes have significance for utilizing immunotherapy in the treatment of PDAC. Clin Cancer Res; 23(15); 4429-40. ©2017 AACR.

Metabolic Reprogramming of Pancreatic Cancer Mediated by CDK4/6 Inhibition Elicits Unique Vulnerabilities.

Due to loss of p16ink4a in pancreatic ductal adenocarcinoma (PDA), pharmacological suppression of CDK4/6 could represent a potent target for treatment. In PDA models, CDK4/6 inhibition had a variable effect on cell cycle but yielded accumulation of ATP and mitochondria. Pharmacological CDK4/6 inhibitors induce cyclin D1 protein levels; however, RB activation was required and sufficient for mitochondrial accumulation. CDK4/6 inhibition stimulated glycolytic and oxidative metabolism and was associated with an increase in mTORC1 activity. MTOR and MEK inhibitors potently cooperate with CDK4/6 inhibition in eliciting cell-cycle exit. However, MTOR inhibition fully suppressed metabolism and yielded apoptosis and suppression of tumor growth in xenograft models. The metabolic state mediated by CDK4/6 inhibition increases mitochondrial number and reactive oxygen species (ROS). Concordantly, the suppression of ROS scavenging or BCL2 antagonists cooperated with CDK4/6 inhibition. Together, these data define the impact of therapeutics on PDA metabolism and provide strategies for converting cytostatic response to tumor cell killing.

Unique metabolic features of pancreatic cancer stroma: relevance to the tumor compartment, prognosis, and invasive potential.

Pancreatic ductal adenocarcinoma (PDAC) has a dismal prognosis. The aggressiveness and therapeutic recalcitrance of this malignancy has been attributed to multiple factors including the influence of an active desmoplastic stroma. How the stromal microenvironment of PDAC contributes to the fatal nature of this disease is not well defined. In the analysis of clinical specimens, we observed diverse expression of the hypoxic marker carbonic anhydrase IX and the lactate transporter MCT4 in the stromal compartment. These stromal features were associated with the epithelial to mesenchymal phenotype in PDAC tumor cells, and with shorter patient survival. Cultured cancer-associated fibroblasts (CAFs) derived from primary PDAC exhibited a high basal level of hypoxia inducible factor 1a (HIF1α) that was both required and sufficient to modulate the expression of MCT4. This event was associated with increased transcription and protein synthesis of HIF1α in CAFs relative to PDAC cell lines, while surprisingly the protein turnover rate was equivalent. CAFs utilized glucose predominantly for glycolytic intermediates, whereas glutamine was the preferred metabolite for the TCA cycle. Unlike PDAC cell lines, CAFs were resistant to glucose withdrawal but sensitive to glutamine depletion. Consistent with the lack of reliance on glucose, CAFs could survive the acute depletion of MCT4. In co-culture and xenograft studies CAFs stimulated the invasive potential and metastatic spread of PDAC cell lines through a mechanism dependent on HIF1α and MCT4. Together, these data indicate that stromal metabolic features influence PDAC tumor cells to promote invasiveness and metastatic potential and associate with poor outcome in patients with PDAC.

Immunologic and Metabolic Features of Pancreatic Ductal Adenocarcinoma Define Prognostic Subtypes of Disease.

PURPOSE: Pancreatic ductal adenocarcinoma (PDA) is associated with an immunosuppressive microenvironment that supports the growth of the malignancy as well as immune system evasion. Here we examine markers of immunosuppression in PDA within the context of the glycolytic tumor microenvironment, their interrelationship with tumor biology and association with overall survival. EXPERIMENTAL DESIGN: We utilized tissue microarrays consisting of 223 PDA patients annotated for clinical stage, tumor size, lymph node involvement, and survival. Expression of CD163, FoxP3, PD-L1, and MCT4 was assessed by IHC and statistical associations were evaluated by univariate and multivariate analysis. Multimarker subtypes were defined by random forest analysis. Mechanistic interactions were evaluated using PDA cell lines and models for myeloid differentiation. RESULTS: PDA exhibits discrete expression of CD163, FoxP3, and PD-L1 with modest individual significance. However, combined low expression of these markers was associated with improved prognosis (P = 0.02). PDA tumor cells altered macrophage phenotype and function, which supported enhanced invasiveness in cell-based models. Lactate efflux mediated by MCT4 was associated with, and required for, the selective conversion of myeloid cells. Correspondingly, MCT4 expression correlated with immune markers in PDA cases, and increased the significance of prognostic subtypes (P = 0.002). CONCLUSIONS: There exists a complex interplay between PDA tumor cells and the host immune system wherein immunosuppression is associated with negative outcome. MCT4 expression, representative of the glycolytic state of PDA, contributes to the phenotypic conversion of myeloid cells. Thus, metabolic status of PDA tumors is an important determinant of the immunosuppressive environment. Clin Cancer Res; 22(14); 3606-17. ©2016 AACR.