BISCUIT: an efficient, standards-compliant tool suite for simultaneous genetic and epigenetic inference in bulk and single-cell studies.

Data from both bulk and single-cell whole-genome DNA methylation experiments are under-utilized in many ways. This is attributable to inefficient mapping of methylation sequencing reads, routinely discarded genetic information, and neglected read-level epigenetic and genetic linkage information. We introduce the BISulfite-seq Command line User Interface Toolkit (BISCUIT) and its companion R/Bioconductor package, biscuiteer, for simultaneous extraction of genetic and epigenetic information from bulk and single-cell DNA methylation sequencing. BISCUIT's performance, flexibility and standards-compliant output allow large, complex experimental designs to be characterized on clinical timescales. BISCUIT is particularly suited for processing data from single-cell DNA methylation assays, with its excellent scalability, efficiency, and ability to greatly enhance mappability, a key challenge for single-cell studies. We also introduce the epiBED format for single-molecule analysis of coupled epigenetic and genetic information, facilitating the study of cellular and tissue heterogeneity from DNA methylation sequencing.

Formation of templated inclusions in a forebrain α-synuclein mouse model is independent of LRRK2.

Leucine-rich repeat kinase 2 (LRRK2) and α-synuclein share enigmatic roles in the pathobiology of Parkinson's disease (PD). LRRK2 mutations are a common genetic cause of PD which, in addition to neurodegeneration, often present with abnormal deposits of α-synuclein in the form of Lewy-related pathology. As Lewy-related pathology is a prominent neuropathologic finding in sporadic PD, the relationship between LRRK2 and α-synuclein has garnered considerable interest. However, whether and how LRRK2 might influence the accumulation of Lewy-related pathology remains poorly understood. Through stereotactic injection of mouse α-synuclein pre-formed fibrils (PFF), we modeled the spread of Lewy-related pathology within forebrain regions where LRRK2 is most highly expressed. The impact of LRRK2 genotype on the formation of α-synuclein inclusions was evaluated at 1-month post-injection. Neither deletion of LRRK2 nor G2019S LRRK2 knockin appreciably altered the burden of α-synuclein pathology at this early timepoint. These observations fail to provide support for a robust pathophysiologic interaction between LRRK2 and α-synuclein in the forebrain in vivo. There was, however, a modest reduction in microglial activation induced by PFF delivery in the hippocampus of LRRK2 knockout mice, suggesting that LRRK2 may contribute to α-synuclein-induced neuroinflammation. Collectively, our data indicate that the pathological accumulation of α-synuclein in the mouse forebrain is largely independent of LRRK2.

NF1 deficiency drives metabolic reprogramming in ER+ breast cancer.

OBJECTIVE: NF1 is a tumor suppressor gene and its protein product, neurofibromin, is a negative regulator of the RAS pathway. NF1 is one of the top driver mutations in sporadic breast cancer such that 27 % of breast cancers exhibit damaging NF1 alterations. NF1 loss-of-function is a frequent event in the genomic evolution of estrogen receptor (ER)+ breast cancer metastasis and endocrine resistance. Individuals with Neurofibromatosis type 1 (NF) - a disorder caused by germline NF1 mutations - have an increased risk of dying from breast cancer [1-4]. NF-related breast cancers are associated with decreased overall survival compared to sporadic breast cancer. Despite numerous studies interrogating the role of RAS mutations in tumor metabolism, no study has comprehensively profiled the NF1-deficient breast cancer metabolome to define patterns of energetic and metabolic reprogramming. The goals of this investigation were (1) to define the role of NF1 deficiency in estrogen receptor-positive (ER+) breast cancer metabolic reprogramming and (2) to identify potential targeted pathway and metabolic inhibitor combination therapies for NF1-deficient ER + breast cancer. METHODS: We employed two ER+ NF1-deficient breast cancer models: (1) an NF1-deficient MCF7 breast cancer cell line to model sporadic breast cancer, and (2) three distinct, Nf1-deficient rat models to model NF-related breast cancer [1]. IncuCyte proliferation analysis was used to measure the effect of NF1 deficiency on cell proliferation and drug response. Protein quantity was assessed by Western Blot analysis. We then used RNAseq to investigate the transcriptional effect of NF1 deficiency on global and metabolism-related transcription. We measured cellular energetics using Agilent Seahorse XF-96 Glyco Stress Test and Mito Stress Test assays. We performed stable isotope labeling and measured [U-13C]-glucose and [U-13C]-glutamine metabolite incorporation and measured total metabolite pools using mass spectrometry. Lastly, we used a Bliss synergy model to investigate NF1-driven changes in targeted and metabolic inhibitor synergy. RESULTS: Our results revealed that NF1 deficiency enhanced cell proliferation, altered neurofibromin expression, and increased RAS and PI3K/AKT pathway signaling while constraining oxidative ATP production and restricting energetic flexibility. Neurofibromin deficiency also increased glutamine influx into TCA intermediates and dramatically increased lipid pools, especially triglycerides (TG). Lastly, NF1 deficiency alters the synergy between metabolic inhibitors and traditional targeted inhibitors. This includes increased synergy with inhibitors targeting glycolysis, glutamine metabolism, mitochondrial fatty acid transport, and TG synthesis. CONCLUSIONS: NF1 deficiency drives metabolic reprogramming in ER+ breast cancer. This reprogramming is characterized by oxidative ATP constraints, glutamine TCA influx, and lipid pool expansion, and these metabolic changes introduce novel metabolic-to-targeted inhibitor synergies.

Stromal mediated DNA damage promotes high grade serous ovarian cancer initiation.

The fundamental steps in high-grade serous ovarian cancer (HGSOC) initiation are unclear, thus providing critical barriers to the development of prevention or early detection strategies for this deadly disease. Increasing evidence demonstrates most HGSOC starts in the fallopian tube epithelium (FTE). Current models propose HGSOC initiates when FTE cells acquire increasing numbers of mutations allowing cells to evolve into serous tubal intraepithelial carcinoma (STIC) precursors and then to full blown cancer. Here we report that epigenetically altered mesenchymal stem cells (termed high risk MSC-hrMSCs) can be detected prior to the formation of ovarian cancer precursor lesions. These hrMSCs drive DNA damage in the form of DNA double strand breaks in FTE cells while also promoting the survival of FTE cells in the face of DNA damage. Indicating the hrMSC may actually drive cancer initiation, we find hrMSCs induce full malignant transformation of otherwise healthy, primary FTE resulting in metastatic cancer in vivo . Further supporting a role for hrMSCs in cancer initiation in humans, we demonstrate that hrMSCs are highly enriched in BRCA1/2 mutation carriers and increase with age. Combined these findings indicate that hrMSCs may incite ovarian cancer initiation. These findings have important implications for ovarian cancer detection and prevention.

Nf1 deficiency modulates the stromal environment in the pretumorigenic rat mammary gland.

Background: Neurofibromin, coded by the NF1 tumor suppressor gene, is the main negative regulator of the RAS pathway and is frequently mutated in various cancers. Women with Neurofibromatosis Type I (NF1)-a tumor predisposition syndrome caused by a germline NF1 mutation-have an increased risk of developing aggressive breast cancer with poorer prognosis. The mechanism by which NF1 mutations lead to breast cancer tumorigenesis is not well understood. Therefore, the objective of this work was to identify stromal alterations before tumor formation that result in the increased risk and poorer outcome seen among NF1 patients with breast cancer. Approach: To accurately model the germline monoallelic NF1 mutations in NF1 patients, we utilized an Nf1-deficient rat model with accelerated mammary development before presenting with highly penetrant breast cancer. Results: We identified increased collagen content in Nf1-deficient rat mammary glands before tumor formation that correlated with age of tumor onset. Additionally, gene expression analysis revealed that Nf1-deficient mature adipocytes in the rat mammary gland have increased collagen expression and shifted to a fibroblast and preadipocyte expression profile. This alteration in lineage commitment was also observed with in vitro differentiation, however, flow cytometry analysis did not show a change in mammary adipose-derived mesenchymal stem cell abundance. Conclusion: Collectively, this study uncovered the previously undescribed role of Nf1 in mammary collagen deposition and regulating adipocyte differentiation. In addition to unraveling the mechanism of tumor formation, further investigation of adipocytes and collagen modifications in preneoplastic mammary glands will create a foundation for developing early detection strategies of breast cancer among NF1 patients.

p53 modulates kinase inhibitor resistance and lineage plasticity in NF1-related MPNSTs.

Malignant peripheral nerve sheath tumors (MPNSTs) are chemotherapy resistant sarcomas that are a leading cause of death in neurofibromatosis type 1 (NF1). Although NF1-related MPNSTs derive from neural crest cell origin, they also exhibit intratumoral heterogeneity. TP53 mutations are associated with significantly decreased survival in MPNSTs, however the mechanisms underlying TP53-mediated therapy responses are unclear in the context of NF1-deficiency. We evaluated the role of two commonly altered genes, MET and TP53, in kinome reprograming and cellular differentiation in preclinical MPNST mouse models. We previously showed that MET amplification occurs early in human MPNST progression and that Trp53 loss abrogated MET-addiction resulting in MET inhibitor resistance. Here we demonstrate a novel mechanism of therapy resistance whereby p53 alters MET stability, localization, and downstream signaling leading to kinome reprogramming and lineage plasticity. Trp53 loss also resulted in a shift from RAS/ERK to AKT signaling and enhanced sensitivity to MEK and mTOR inhibition. In response to MET, MEK and mTOR inhibition, we observed broad and heterogeneous activation of key differentiation genes in Trp53-deficient lines suggesting Trp53 loss also impacts lineage plasticity in MPNSTs. These results demonstrate the mechanisms by which p53 loss alters MET dependency and therapy resistance in MPNSTS through kinome reprogramming and phenotypic flexibility.