In vivo CRISPR screening reveals nutrient signaling processes underpinning CD8+ T cell fate decisions.

How early events in effector T cell (TEFF) subsets tune memory T cell (TMEM) responses remains incompletely understood. Here, we systematically investigated metabolic factors in fate determination of TEFF and TMEM cells using in vivo pooled CRISPR screening, focusing on negative regulators of TMEM responses. We found that amino acid transporters Slc7a1 and Slc38a2 dampened the magnitude of TMEM differentiation, in part through modulating mTORC1 signaling. By integrating genetic and systems approaches, we identified cellular and metabolic heterogeneity among TEFF cells, with terminal effector differentiation associated with establishment of metabolic quiescence and exit from the cell cycle. Importantly, Pofut1 (protein-O-fucosyltransferase-1) linked GDP-fucose availability to downstream Notch-Rbpj signaling, and perturbation of this nutrient signaling axis blocked terminal effector differentiation but drove context-dependent TEFF proliferation and TMEM development. Our study establishes that nutrient uptake and signaling are key determinants of T cell fate and shape the quantity and quality of TMEM responses.

CXCR6 orchestrates brain CD8+ T cell residency and limits mouse Alzheimer's disease pathology.

Neurodegenerative diseases, including Alzheimer's disease (AD), are characterized by innate immune-mediated inflammation, but functional and mechanistic effects of the adaptive immune system remain unclear. Here we identify brain-resident CD8+ T cells that coexpress CXCR6 and PD-1 and are in proximity to plaque-associated microglia in human and mouse AD brains. We also establish that CD8+ T cells restrict AD pathologies, including ß-amyloid deposition and cognitive decline. Ligand-receptor interaction analysis identifies CXCL16-CXCR6 intercellular communication between microglia and CD8+ T cells. Further, Cxcr6 deficiency impairs accumulation, tissue residency programming and clonal expansion of brain PD-1+CD8+ T cells. Ablation of Cxcr6 or CD8+ T cells ultimately increases proinflammatory cytokine production from microglia, with CXCR6 orchestrating brain CD8+ T cell-microglia colocalization. Collectively, our study reveals protective roles for brain CD8+ T cells and CXCR6 in mouse AD pathogenesis and highlights that microenvironment-specific, intercellular communication orchestrates tissue homeostasis and protection from neuroinflammation.

G3BP1 Is a Tunable Switch that Triggers Phase Separation to Assemble Stress Granules.

The mechanisms underlying ribonucleoprotein (RNP) granule assembly, including the basis for establishing and maintaining RNP granules with distinct composition, are unknown. One prominent type of RNP granule is the stress granule (SG), a dynamic and reversible cytoplasmic assembly formed in eukaryotic cells in response to stress. Here, we show that SGs assemble through liquid-liquid phase separation (LLPS) arising from interactions distributed unevenly across a core protein-RNA interaction network. The central node of this network is G3BP1, which functions as a molecular switch that triggers RNA-dependent LLPS in response to a rise in intracellular free RNA concentrations. Moreover, we show that interplay between three distinct intrinsically disordered regions (IDRs) in G3BP1 regulates its intrinsic propensity for LLPS, and this is fine-tuned by phosphorylation within the IDRs. Further regulation of SG assembly arises through positive or negative cooperativity by extrinsic G3BP1-binding factors that strengthen or weaken, respectively, the core SG network.

Metabolic control of TFH cells and humoral immunity by phosphatidylethanolamine.

T follicular helper (TFH) cells are crucial for B cell-mediated humoral immunity1. Although transcription factors such as BCL6 drive the differentiation of TFH cells2,3, it is unclear whether and how post-transcriptional and metabolic programs enforce TFH cell programming. Here we show that the cytidine diphosphate (CDP)-ethanolamine pathway co-ordinates the expression and localization of CXCR5 with the responses of TFH cells and humoral immunity. Using in vivo CRISPR-Cas9 screening and functional validation in mice, we identify ETNK1, PCYT2, and SELENOI-enzymes in the CDP-ethanolamine pathway for de novo synthesis of phosphatidylethanolamine (PE)-as selective post-transcriptional regulators of TFH cell differentiation that act by promoting the surface expression and functional effects of CXCR5. TFH cells exhibit unique lipid metabolic programs and PE is distributed to the outer layer of the plasma membrane, where it colocalizes with CXCR5. De novo synthesis of PE through the CDP-ethanolamine pathway co-ordinates these events to prevent the internalization and degradation of CXCR5. Genetic deletion of Pcyt2, but not of Pcyt1a (which mediates the CDP-choline pathway), in activated T cells impairs the differentiation of TFH cells, and this is associated with reduced humoral immune responses. Surface levels of PE and CXCR5 expression on B cells also depend on Pcyt2. Our results reveal that phospholipid metabolism orchestrates post-transcriptional mechanisms for TFH cell differentiation and humoral immunity, highlighting the metabolic control of context-dependent immune signalling and effector programs.

CRISPR screens unveil signal hubs for nutrient licensing of T cell immunity.

Nutrients are emerging regulators of adaptive immunity1. Selective nutrients interplay with immunological signals to activate mechanistic target of rapamycin complex 1 (mTORC1), a key driver of cell metabolism2-4, but how these environmental signals are integrated for immune regulation remains unclear. Here we use genome-wide CRISPR screening combined with protein-protein interaction networks to identify regulatory modules that mediate immune receptor- and nutrient-dependent signalling to mTORC1 in mouse regulatory T (Treg) cells. SEC31A is identified to promote mTORC1 activation by interacting with the GATOR2 component SEC13 to protect it from SKP1-dependent proteasomal degradation. Accordingly, loss of SEC31A impairs T cell priming and Treg suppressive function in mice. In addition, the SWI/SNF complex restricts expression of the amino acid sensor CASTOR1, thereby enhancing mTORC1 activation. Moreover, we reveal that the CCDC101-associated SAGA complex is a potent inhibitor of mTORC1, which limits the expression of glucose and amino acid transporters and maintains T cell quiescence in vivo. Specific deletion of Ccdc101 in mouse Treg cells results in uncontrolled inflammation but improved antitumour immunity. Collectively, our results establish epigenetic and post-translational mechanisms that underpin how nutrient transporters, sensors and transducers interplay with immune signals for three-tiered regulation of mTORC1 activity and identify their pivotal roles in licensing T cell immunity and immune tolerance.

Genomic landscape of Down syndrome-associated acute lymphoblastic leukemia.

Trisomy 21, the genetic cause of Down syndrome (DS), is the most common congenital chromosomal anomaly. It is associated with a 20-fold increased risk of acute lymphoblastic leukemia (ALL) during childhood and results in distinctive leukemia biology. To comprehensively define the genomic landscape of DS-ALL, we performed whole-genome sequencing and whole-transcriptome sequencing (RNA-Seq) on 295 cases. Our integrated genomic analyses identified 15 molecular subtypes of DS-ALL, with marked enrichment of CRLF2-r, IGH::IGF2BP1, and C/EBP altered (C/EBPalt) subtypes compared with 2257 non-DS-ALL cases. We observed abnormal activation of the CEBPD, CEBPA, and CEBPE genes in 10.5% of DS-ALL cases via a variety of genomic mechanisms, including chromosomal rearrangements and noncoding mutations leading to enhancer hijacking. A total of 42.3% of C/EBP-activated DS-ALL also have concomitant FLT3 point mutations or insertions/deletions, compared with 4.1% in other subtypes. CEBPD overexpression enhanced the differentiation of mouse hematopoietic progenitor cells into pro-B cells in vitro, particularly in a DS genetic background. Notably, recombination-activating gene-mediated somatic genomic abnormalities were common in DS-ALL, accounting for a median of 27.5% of structural alterations, compared with 7.7% in non-DS-ALL. Unsupervised hierarchical clustering analyses of CRLF2-rearranged DS-ALL identified substantial heterogeneity within this group, with the BCR::ABL1-like subset linked to an inferior event-free survival, even after adjusting for known clinical risk factors. These results provide important insights into the biology of DS-ALL and point to opportunities for targeted therapy and treatment individualization.

Targeting REGNASE-1 programs long-lived effector T cells for cancer therapy.

Adoptive cell therapy represents a new paradigm in cancer immunotherapy, but it can be limited by the poor persistence and function of transferred T cells1. Here we use an in vivo pooled CRISPR-Cas9 mutagenesis screening approach to demonstrate that, by targeting REGNASE-1, CD8+ T cells are reprogrammed to long-lived effector cells with extensive accumulation, better persistence and robust effector function in tumours. REGNASE-1-deficient CD8+ T cells show markedly improved therapeutic efficacy against mouse models of melanoma and leukaemia. By using a secondary genome-scale CRISPR-Cas9 screening, we identify BATF as the key target of REGNASE-1 and as a rheostat that shapes antitumour responses. Loss of BATF suppresses the increased accumulation and mitochondrial fitness of REGNASE-1-deficient CD8+ T cells. By contrast, the targeting of additional signalling factors-including PTPN2 and SOCS1-improves the therapeutic efficacy of REGNASE-1-deficient CD8+ T cells. Our findings suggest that T cell persistence and effector function can be coordinated in tumour immunity and point to avenues for improving the efficacy of adoptive cell therapy for cancer.

miR-424(322)/503 is a breast cancer tumor suppressor whose loss promotes resistance to chemotherapy.

The female mammary gland is a very dynamic organ that undergoes continuous tissue remodeling during adulthood. Although it is well established that the number of menstrual cycles and pregnancy (in this case transiently) increase the risk of breast cancer, the reasons are unclear. Growing clinical and experimental evidence indicates that improper involution plays a role in the development of this malignancy. Recently, we described the miR-424(322)/503 cluster as an important regulator of mammary epithelial involution after pregnancy. Here, through the analysis of ∼3000 primary tumors, we show that miR-424(322)/503 is commonly lost in a subset of aggressive breast cancers and describe the genetic aberrations that inactivate its expression. Furthermore, through the use of a knockout mouse model, we demonstrate for the first time that loss of miR-424(322)/503 promotes breast tumorigenesis in vivo. Remarkably, we found that loss of miR-424(322)/503 promotes chemoresistance due to the up-regulation of two of its targets: BCL-2 and insulin-like growth factor-1 receptor (IGF1R). Importantly, targeted therapies blocking the aberrant activity of these targets restore sensitivity to chemotherapy. Overall, our studies reveal miR-424(322)/503 as a tumor suppressor in breast cancer and provide a link between mammary epithelial involution, tumorigenesis, and the phenomenon of chemoresistance.

Single-cell analysis reveals the Comma-1D cell line as a unique model for mammary gland development and breast cancer.

The mammary gland epithelial tree contains two distinct cell populations, luminal and basal. The investigation of how this heterogeneity is developed and how it influences tumorigenesis has been hampered by the need to perform studies on these populations using animal models. Comma-1D is an immortalized mouse mammary epithelial cell line that has unique morphogenetic properties. By performing single-cell RNA-seq studies, we found that Comma-1D cultures consist of two main populations with luminal and basal features, and a smaller population with mixed lineage and bipotent characteristics. We demonstrated that multiple transcription factors associated with the differentiation of the mammary epithelium in vivo also modulate this process in Comma-1D cultures. Additionally, we found that only cells with luminal features were able to acquire transformed characteristics after an oncogenic HER2 (also known as ERBB2) mutant was introduced in their genomes. Overall, our studies characterize, at a single-cell level, the heterogeneity of the Comma-1D cell line and illustrate how Comma-1D cells can be used as an experimental model to study both the differentiation and the transformation processes in vitro.

Impact of T-cell immunity on chemotherapy response in childhood acute lymphoblastic leukemia.

Although acute lymphoblastic leukemia (ALL) is highly responsive to chemotherapy, it is unknown how or which host immune factors influence the long-term remission of this cancer. To this end, we systematically evaluated the effects of T-cell immunity on Ph+ ALL therapy outcomes. Using a murine Arf-/-BCR-ABL1 B-cell ALL model, we showed that loss of T cells in the host drastically increased leukemia relapse after dasatinib or cytotoxic chemotherapy. Although ABL1 mutations emerged early during dasatinib treatment in both immunocompetent and immunocompromised hosts, T-cell immunity was essential for suppressing the outgrowth of drug-resistant leukemia. Bulk and single-cell transcriptome profiling of T cells during therapy pointed to the activation of type 1 immunity-related cytokine signaling being linked to long-term leukemia remission in mice. Consistent with these observations, interferon γ and interleukin 12 directly modulated dasatinib antileukemia efficacy in vivo. Finally, we evaluated peripheral blood immune cell composition in 102 children with ALL during chemotherapy and observed a significant association of T-cell abundance with treatment outcomes. Together, these results suggest that T-cell immunity plays pivotal roles in maintaining long-term remission of ALL, highlighting that the interplay between host immunity and drug resistance can be harnessed to improve ALL chemotherapy outcomes.

Hippo/Mst signalling couples metabolic state and immune function of CD8α+ dendritic cells.

Dendritic cells orchestrate the crosstalk between innate and adaptive immunity. CD8α+ dendritic cells present antigens to CD8+ T cells and elicit cytotoxic T cell responses to viruses, bacteria and tumours 1 . Although lineage-specific transcriptional regulators of CD8α+ dendritic cell development have been identified 2 , the molecular pathways that selectively orchestrate CD8α+ dendritic cell function remain elusive. Moreover, metabolic reprogramming is important for dendritic cell development and activation3,4, but metabolic dependence and regulation of dendritic cell subsets are largely uncharacterized. Here we use a data-driven systems biology algorithm (NetBID) to identify a role of the Hippo pathway kinases Mst1 and Mst2 (Mst1/2) in selectively programming CD8α+ dendritic cell function and metabolism. Our NetBID analysis reveals a marked enrichment of the activities of Hippo pathway kinases in CD8α+ dendritic cells relative to CD8α- dendritic cells. Dendritic cell-specific deletion of Mst1/2-but not Lats1 and Lats2 (Lats1/2) or Yap and Taz (Yap/Taz), which mediate canonical Hippo signalling-disrupts homeostasis and function of CD8+ T cells and anti-tumour immunity. Mst1/2-deficient CD8α+ dendritic cells are impaired in presentation of extracellular proteins and cognate peptides to prime CD8+ T cells, while CD8α- dendritic cells that lack Mst1/2 have largely normal function. Mechanistically, compared to CD8α- dendritic cells, CD8α+ dendritic cells exhibit much stronger oxidative metabolism and critically depend on Mst1/2 signalling to maintain bioenergetic activities and mitochondrial dynamics for their functional capacities. Further, selective expression of IL-12 by CD8α+ dendritic cells depends on Mst1/2 and the crosstalk with non-canonical NF-κB signalling. Our findings identify Mst1/2 as selective drivers of CD8α+ dendritic cell function by integrating metabolic activity and cytokine signalling, and highlight that the interplay between immune signalling and metabolic reprogramming underlies the unique functions of dendritic cell subsets.

A developmentally prometastatic niche to hepatoblastoma in neonatal liver mediated by the Cxcl1/Cxcr2 axis.

BACKGROUND AND AIMS: Hepatoblastoma (HB) is the most common pediatric liver cancer. Its predominant occurrence in very young children led us to investigate whether the neonatal liver provides a protumorigenic niche to HB development. APPROACH AND RESULTS: HB development was compared between orthotopic transplantation models established in postnatal day 5 (P5) and 60 (P60) mice (P5Tx and P60Tx models). Single-cell RNA-sequencing (sc-RNAseq) was performed using tumor and liver tissues from both models and the top candidate cell types and genes identified are investigated for their roles in HB cell growth, migration, and survival. CONCLUSIONS: We found that various HB cell lines including HepG2 cells were consistently and considerably more tumorigenic and metastatic in the P5Tx model than in the P60Tx models. Sc-RNAseq of the P5Tx and P60Tx HepG2 models revealed that the P5Tx tumor was more hypoxic and had a larger number of activated hepatic stellate cells (aHSCs) in the tumor-surrounding liver that express significantly higher levels of Cxcl1 than those from the P60Tx model. We found these differences were developmentally present in normal P5 and P60 liver. We showed that the Cxcl1/Cxcr2 axis mediated HB cell migration and was critical to HB cell survival under hypoxia. Treating HepG2 P60Tx model with recombinant CXCL1 protein induced intrahepatic and pulmonary metastasis and CXCR2 knockout (KO) in HepG2 cells abolished their metastatic potential in the P5Tx model. Lastly, we showed that in tumors from patients with metastatic HB, there was a similar larger population of aHSCs in the tumor-surrounding liver than in localized tumors, and tumor hypoxia was uniquely associated with prognosis of patients with HB among pediatric cancers. We demonstrated that the neonatal liver provides a prometastatic niche to HB development through the Cxcl1/Cxcr2 axis.

Modeling and targeting of erythroleukemia by hematopoietic genome editing.

Acute erythroid leukemia (AEL) is characterized by a distinct morphology, mutational spectrum, lack of preclinical models, and poor prognosis. Here, using multiplexed genome editing of mouse hematopoietic stem and progenitor cells and transplant assays, we developed preclinical models of AEL and non-erythroid acute leukemia and describe the central role of mutational cooperativity in determining leukemia lineage. Different combination of mutations in Trp53, Bcor, Dnmt3a, Rb1, and Nfix resulted in the development of leukemia with an erythroid phenotype, accompanied by the acquisition of alterations in signaling and transcription factor genes that recapitulate human AEL by cross-species genomic analysis. Clonal expansion during tumor evolution was driven by mutational cooccurrence, with clones harboring a higher number of founder and secondary lesions (eg, mutations in signaling genes) showing greater evolutionary fitness. Mouse and human AEL exhibited deregulation of genes regulating erythroid development, notably Gata1, Klf1, and Nfe2, driven by the interaction of mutations of the epigenetic modifiers Dnmt3a and Tet2 that perturbed methylation and thus expression of lineage-specific transcription factors. The established mouse leukemias were used as a platform for drug screening. Drug sensitivity was associated with the leukemia genotype, with the poly (ADP-ribose) polymerase inhibitor talazoparib and the demethylating agent decitabine efficacious in Trp53/Bcor-mutant AEL, CDK7/9 inhibitors in Trp53/Bcor/Dnmt3a-mutant AEL, and gemcitabine and bromodomain inhibitors in NUP98-KDM5A leukemia. In conclusion, combinatorial genome editing has shown the interplay of founding and secondary genetic alterations in phenotype and clonal evolution, epigenetic regulation of lineage-specific transcription factors, and therapeutic tractability in erythroid leukemogenesis.

miR-424/503 modulates Wnt/ß-catenin signaling in the mammary epithelium by targeting LRP6.

During the female lifetime, the expansion of the epithelium dictated by the ovarian cycles is supported by a transient increase in the mammary epithelial stem cell population (MaSCs). Notably, activation of Wnt/ß-catenin signaling is an important trigger for MaSC expansion. Here, we report that the miR-424/503 cluster is a modulator of canonical Wnt signaling in the mammary epithelium. We show that mammary tumors of miR-424(322)/503-depleted mice exhibit activated Wnt/ß-catenin signaling. Importantly, we show a strong association between miR-424/503 deletion and breast cancers with high levels of Wnt/ß-catenin signaling. Moreover, miR-424/503 cluster is required for Wnt-mediated MaSC expansion induced by the ovarian cycles. Lastly, we show that miR-424/503 exerts its function by targeting two binding sites at the 3'UTR of the LRP6 co-receptor and reducing its expression. These results unveil an unknown link between the miR-424/503, regulation of Wnt signaling, MaSC fate, and tumorigenesis.

Inhibition of the autocrine IL-6-JAK2-STAT3-calprotectin axis as targeted therapy for HR-/HER2+ breast cancers.

HER2-positive (HER2(+)) breast adenocarcinomas are a heterogeneous group in which hormone receptor (HR) status influences therapeutic decisions and patient outcome. By combining genome-wide RNAi screens with regulatory network analysis, we identified STAT3 as a critically activated master regulator of HR(-)/HER2(+) tumors, eliciting tumor dependency in these cells. Mechanistically, HR(-)/HER2(+) cells secrete high levels of the interleukin-6 (IL-6) cytokine, inducing the activation of STAT3, which in turn promotes a second autocrine stimulus to increase S100A8/9 complex (calprotectin) production and secretion. Increased calprotectin levels activate signaling pathways involved in proliferation and resistance. Importantly, we demonstrated that inhibition of the IL-6-Janus kinase 2 (JAK2)-STAT3-calprotectin axis with FDA-approved drugs, alone and in combination with HER2 inhibitors, reduced the tumorigenicity of HR(-)/HER2(+) breast cancers, opening novel targeted therapeutic opportunities.

SJARACNe: a scalable software tool for gene network reverse engineering from big data.

SUMMARY: Over the last two decades, we have observed an exponential increase in the number of generated array or sequencing-based transcriptomic profiles. Reverse engineering of biological networks from high-throughput gene expression profiles has been one of the grand challenges in systems biology. The Algorithm for the Reconstruction of Accurate Cellular Networks (ARACNe) represents one of the most effective and widely-used tools to address this challenge. However, existing ARACNe implementations do not efficiently process big input data with thousands of samples. Here we present an improved implementation of the algorithm, SJARACNe, to solve this big data problem, based on sophisticated software engineering. The new scalable SJARACNe package achieves a dramatic improvement in computational performance in both time and memory usage and implements new features while preserving the network inference accuracy of the original algorithm. Given that large-sampled transcriptomic data is increasingly available and ARACNe is extremely demanding for network reconstruction, the scalable SJARACNe will allow even researchers with modest computational resources to efficiently construct complex regulatory and signaling networks from thousands of gene expression profiles. AVAILABILITY AND IMPLEMENTATION: SJARACNe is implemented in C++ (computational core) and Python (pipelining scripting wrapper, ≥3.6.1). It is freely available at https://github.com/jyyulab/SJARACNe. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

Silymarin's Inhibition and Treatment Effects for Alzheimer's Disease.

As a longstanding problem, Alzheimer's disease (AD) has stymied researchers in the medical field with its increasing incidence and enormous treatment difficulty. Silymarin has always been valued by researchers for its good efficacy and safety in treating liver disease. Recent studies have shown that silymarin also has good pharmacological activity in the nervous system, especially for the treatment of AD. Silymarin can control the production of Aß by inhibiting the precursor substance of Aß (ß-amyloid precursor protein), and it can inhibit the polymerization of Aß. Silymarin can also increase the acetylcholine content in the nervous system by inhibiting cholinesterase activity. At the same time, it also has the effect of resisting oxidative stress and the inflammatory response of the nervous system. These pharmacological activities contribute to the inhibition of the onset of AD. The good efficacy of silymarin on AD and its high safety and availability give it huge potential for the treatment of AD.

A novel AP2/ERF family transcription factor from Glycine soja, GsERF71, is a DNA binding protein that positively regulates alkaline stress tolerance in Arabidopsis.

KEY MESSAGE: Here we first found that GsERF71, an ERF factor from wild soybean could increase plant alkaline stress tolerance by up-regulating H+-ATPase and by modifing the accumulation of Auxin. Alkaline soils are widely distributed all over the world and greatly limit plant growth and development. In our previous transcriptome analyses, we have identified several ERF (ethylene-responsive factor) genes that responded strongly to bicarbonate stress in the roots of wild soybean G07256 (Glycine soja). In this study, we cloned and functionally characterized one of the genes, GsERF71. When expressed in epidermal cells of onion, GsERF71 localized to the nucleus. It can activate the reporters in yeast cells, and the C-terminus of 170 amino acids is essential for its transactivation activity. Yeast one-hybrid and EMSA assays indicated that GsERF71 specifically binds to the cis-acting elements of the GCC-box, suggesting that GsERF71 may participate in the regulation of transcription of the relevant biotic and abiotic stress-related genes. Furthermore, transgenic Arabidopsis plants overexpressing GsERF71 showed significantly higher tolerance to bicarbonate stress generated by NaHCO3 or KHCO3 than the wild type (WT) plants, i.e., the transgenic plants had greener leaves, longer roots, higher total chlorophyll contents and lower MDA contents. qRT-PCR and rhizosphere acidification assays indicated that the expression level and activity of H+-ATPase (AHA2) were enhanced in the transgenic plants under alkaline stress. Further analysis indicated that the expression of auxin biosynthetic genes and IAA contents were altered to a lower extent in the roots of transgenic plants than WT plants under alkaline stress in a short-term. Together, our data suggest that GsERF71 enhances the tolerance to alkaline stress by up-regulating the expression levels of H+-ATPase and by modifying auxin accumulation in transgenic plants.

ScreenBEAM: a novel meta-analysis algorithm for functional genomics screens via Bayesian hierarchical modeling.

MOTIVATION: Functional genomics (FG) screens, using RNAi or CRISPR technology, have become a standard tool for systematic, genome-wide loss-of-function studies for therapeutic target discovery. As in many large-scale assays, however, off-target effects, variable reagents' potency and experimental noise must be accounted for appropriately control for false positives. Indeed, rigorous statistical analysis of high-throughput FG screening data remains challenging, particularly when integrative analyses are used to combine multiple sh/sgRNAs targeting the same gene in the library. METHOD: We use large RNAi and CRISPR repositories that are publicly available to evaluate a novel meta-analysis approach for FG screens via Bayesian hierarchical modeling, Screening Bayesian Evaluation and Analysis Method (ScreenBEAM). RESULTS: Results from our analysis show that the proposed strategy, which seamlessly combines all available data, robustly outperforms classical algorithms developed for microarray data sets as well as recent approaches designed for next generation sequencing technologies. Remarkably, the ScreenBEAM algorithm works well even when the quality of FG screens is relatively low, which accounts for about 80-95% of the public datasets. AVAILABILITY AND IMPLEMENTATION: R package and source code are available at: https://github.com/jyyu/ScreenBEAM. CONTACT: ac2248@columbia.edu, jose.silva@mssm.edu, yujiyang@gmail.com SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

HDAC6 activity is a non-oncogene addiction hub for inflammatory breast cancers.

INTRODUCTION: Inflammatory breast cancer (IBC) is the most lethal form of breast cancers with a 5-year survival rate of only 40 %. Despite its lethality, IBC remains poorly understood which has greatly limited its therapeutic management. We thus decided to utilize an integrative functional genomic strategy to identify the Achilles' heel of IBC cells. METHODS: We have pioneered the development of genetic tools as well as experimental and analytical strategies to perform RNAi-based loss-of-function studies at a genome-wide level. Importantly, we and others have demonstrated that these functional screens are able to identify essential functions linked to certain cancer phenotypes. Thus, we decided to use this approach to identify IBC specific sensitivities. RESULTS: We identified and validated HDAC6 as a functionally necessary gene to maintain IBC cell viability, while being non-essential for other breast cancer subtypes. Importantly, small molecule inhibitors for HDAC6 already exist and are in clinical trials for other tumor types. We thus demonstrated that Ricolinostat (ACY1215), a leading HDAC6 inhibitor, efficiently controls IBC cell proliferation both in vitro and in vivo. Critically, functional HDAC6 dependency is not associated with genomic alterations at its locus and thus represents a non-oncogene addiction. Despite HDAC6 not being overexpressed, we found that its activity is significantly higher in IBC compared to non-IBC cells, suggesting a possible rationale supporting the observed dependency. CONCLUSION: Our finding that IBC cells are sensitive to HDAC6 inhibition provides a foundation to rapidly develop novel, efficient, and well-tolerated targeted therapy strategies for IBC patients.