A clinically applicable integrative molecular classification of meningiomas.

Meningiomas are the most common primary intracranial tumour in adults1. Patients with symptoms are generally treated with surgery as there are no effective medical therapies. The World Health Organization histopathological grade of the tumour and the extent of resection at surgery (Simpson grade) are associated with the recurrence of disease; however, they do not accurately reflect the clinical behaviour of all meningiomas2. Molecular classifications of meningioma that reliably reflect tumour behaviour and inform on therapies are required. Here we introduce four consensus molecular groups of meningioma by combining DNA somatic copy-number aberrations, DNA somatic point mutations, DNA methylation and messenger RNA abundance in a unified analysis. These molecular groups more accurately predicted clinical outcomes compared with existing classification schemes. Each molecular group showed distinctive and prototypical biology (immunogenic, benign NF2 wild-type, hypermetabolic and proliferative) that informed therapeutic options. Proteogenomic characterization reinforced the robustness of the newly defined molecular groups and uncovered highly abundant and group-specific protein targets that we validated using immunohistochemistry. Single-cell RNA sequencing revealed inter-individual variations in meningioma as well as variations in intrinsic expression programs in neoplastic cells that mirrored the biology of the molecular groups identified.

Metabologenomic characterization uncovers a clinically aggressive IDH mutant glioma subtype.

Mutations in the pivotal metabolic isocitrate dehydrogenase (IDH) enzymes are recognized to drive the molecular footprint of diffuse gliomas, and patients with IDH mutant gliomas have overall favorable outcomes compared to patients with IDH wild-type tumors. However, survival still varies widely among patients with IDH mutated tumors. Here, we aimed to characterize molecular signatures that explain the range of IDH mutant gliomas. By integrating matched epigenome-wide methylome, transcriptome, and global metabolome data in 154 patients with gliomas, we identified a group of IDH mutant gliomas with globally altered metabolism that resembled IDH wild-type tumors. IDH-mutant gliomas with altered metabolism have significantly shorter overall survival from their IDH mutant counterparts that is not fully accounted for by recognized molecular prognostic markers of CDKN2A/B loss and glioma CpG Island Methylator Phenotype (GCIMP) status. IDH-mutant tumors with dysregulated metabolism harbored distinct epigenetic alterations that converged to drive proliferative and stem-like transcriptional profiles, providing a window to target novel dependencies in gliomas.

Increased mRNA expression of CDKN2A is a transcriptomic marker of clinically aggressive meningiomas.

Homozygous deletion of CDKN2A/B was recently incorporated into the World Health Organization classification for grade 3 meningiomas. While this marker is overall rare in meningiomas, its relationship to other CDKN2A alterations on a transcriptomic, epigenomic, and copy number level has not yet been determined. We therefore utilized multidimensional molecular data of 1577 meningioma samples from 6 independent cohorts enriched for clinically aggressive meningiomas to comprehensively interrogate the spectrum of CDKN2A alterations through DNA methylation, copy number variation, transcriptomics, and proteomics using an integrated molecular approach. Homozygous CDKN2A/B deletions were identified in only 7.1% of cases but were associated with significantly poorer outcomes compared to tumors without these deletions. Heterozygous CDKN2A/B deletions were identified in 2.6% of cases and had similarly poor outcomes as those with homozygous deletions. Among tumors with intact CDKN2A/B (without a homozygous or heterozygous deletion), we found a distinct difference in outcome based on mRNA expression of CDKN2A, with meningiomas that had elevated mRNA expression (CDKN2Ahigh) having a significantly shorter time to recurrence. The expression of CDKN2A was independently prognostic after accounting for copy number loss and consistently increased with WHO grade and more aggressive molecular and methylation groups irrespective of cohort. Despite the discordant and mutually exclusive status of the CDKN2A gene in these groups, both CDKN2Ahigh meningiomas and meningiomas with CDKN2A deletions were enriched for similar cell cycle pathways but at different checkpoints. High mRNA expression of CDKN2A was also associated with gene hypermethylation, Rb-deficiency, and lack of response to CDK inhibition. p16 immunohistochemistry could not reliably differentiate between meningiomas with and without CDKN2A deletions but appeared to correlate better with mRNA expression. These findings support the role of CDKN2A mRNA expression as a biomarker of clinically aggressive meningiomas with potential therapeutic implications.

The multiomic landscape of meningiomas: a review and update.

PURPOSE: Meningiomas are the most common primary brain tumor in adults. Traditionally they have been understudied compared to other central nervous system (CNS) tumors. However over the last decade, there has been renewed interest in uncovering the molecular topography of these tumors, with landmark studies identifying key driver alterations contributing to meningioma development and progression. Recent work from several independent research groups have integrated different genomic and epigenomic platforms to develop a molecular-based classification scheme for meningiomas that could supersede histopathological grading in terms of diagnostic accuracy, biological relevance, and outcome prediction, keeping pace with contemporary grading schemes for other CNS tumors including gliomas and medulloblastomas. METHODS: Here we summarize the studies that have uncovered key alterations in meningiomas which builds towards the discovery of consensus molecular groups in meningiomas by integrating these findings. These groups supersede WHO grade and other clinical factors in being able to accurately predict tumor biology and clinical outcomes following surgery. RESULTS: Despite differences in the nomenclature of recently uncovered molecular groups across different studies, the biological similarities between these groups enables us to likely reconciliate these groups into four consensus molecular groups: two benign groups largely dichotomized by NF2-status, and two clinically aggressive groups defined by their hypermetabolic transcriptome, and by their preponderance of proliferative, cell-cycling pathways respectively. CONCLUSION: Future work, including by our group and others are underway to validate these molecular groups and harmonize the nomenclature for routine clinical use.

Cortical-Bone Fragility--Insights from sFRP4 Deficiency in Pyle's Disease.

BACKGROUND: Cortical-bone fragility is a common feature in osteoporosis that is linked to nonvertebral fractures. Regulation of cortical-bone homeostasis has proved elusive. The study of genetic disorders of the skeleton can yield insights that fuel experimental therapeutic approaches to the treatment of rare disorders and common skeletal ailments. METHODS: We evaluated four patients with Pyle's disease, a genetic disorder that is characterized by cortical-bone thinning, limb deformity, and fractures; two patients were examined by means of exome sequencing, and two were examined by means of Sanger sequencing. After a candidate gene was identified, we generated a knockout mouse model that manifested the phenotype and studied the mechanisms responsible for altered bone architecture. RESULTS: In all affected patients, we found biallelic truncating mutations in SFRP4, the gene encoding secreted frizzled-related protein 4, a soluble Wnt inhibitor. Mice deficient in Sfrp4, like persons with Pyle's disease, have increased amounts of trabecular bone and unusually thin cortical bone, as a result of differential regulation of Wnt and bone morphogenetic protein (BMP) signaling in these two bone compartments. Treatment of Sfrp4-deficient mice with a soluble Bmp2 receptor (RAP-661) or with antibodies to sclerostin corrected the cortical-bone defect. CONCLUSIONS: Our study showed that Pyle's disease was caused by a deficiency of sFRP4, that cortical-bone and trabecular-bone homeostasis were governed by different mechanisms, and that sFRP4-mediated cross-regulation between Wnt and BMP signaling was critical for achieving proper cortical-bone thickness and stability. (Funded by the Swiss National Foundation and the National Institutes of Health.).

ChemiRs: a web application for microRNAs and chemicals.

BACKGROUND: MicroRNAs (miRNAs) are about 22 nucleotides, non-coding RNAs that affect various cellular functions, and play a regulatory role in different organisms including human. Until now, more than 2500 mature miRNAs in human have been discovered and registered, but still lack of information or algorithms to reveal the relations among miRNAs, environmental chemicals and human health. Chemicals in environment affect our health and daily life, and some of them can lead to diseases by inferring biological pathways. RESULTS: We develop a creditable online web server, ChemiRs, for predicting interactions and relations among miRNAs, chemicals and pathways. The database not only compares gene lists affected by chemicals and miRNAs, but also incorporates curated pathways to identify possible interactions. CONCLUSIONS: Here, we manually retrieved associations of miRNAs and chemicals from biomedical literature. We developed an online system, ChemiRs, which contains miRNAs, diseases, Medical Subject Heading (MeSH) terms, chemicals, genes, pathways and PubMed IDs. We connected each miRNA to miRBase, and every current gene symbol to HUGO Gene Nomenclature Committee (HGNC) for genome annotation. Human pathway information is also provided from KEGG and REACTOME databases. Information about Gene Ontology (GO) is queried from GO Online SQL Environment (GOOSE). With a user-friendly interface, the web application is easy to use. Multiple query results can be easily integrated and exported as report documents in PDF format. Association analysis of miRNAs and chemicals can help us understand the pathogenesis of chemical components. ChemiRs is freely available for public use at http://omics.biol.ntnu.edu.tw/ChemiRs .

Targeted Pten deletion plus p53-R270H mutation in mouse mammary epithelium induces aggressive claudin-low and basal-like breast cancer.

BACKGROUND: Triple-negative breast cancer (TNBC), an aggressive disease comprising several subtypes including basal-like and claudin-low, involves frequent deletions or point mutations in TP53, as well as loss of PTEN. We previously showed that combined deletion of both tumor suppressors in the mouse mammary epithelium invariably induced claudin-low-like TNBC. The effect of p53 mutation plus Pten deletion on mammary tumorigenesis and whether this combination can induce basal-like TNBC in the mouse are unknown. METHODS: WAP-Cre:Pten(f/f):p53(lox.stop.lox_R270H) composite mice were generated in which Pten is deleted and a p53-R270H mutation in the DNA-binding domain is induced upon expression of Cre-recombinase in pregnancy-identified alveolar progenitors. Tumors were characterized by histology, marker analysis, transcriptional profiling [GEO-GSE75989], bioinformatics, high-throughput (HTP) FDA drug screen as well as orthotopic injection to quantify tumor-initiating cells (TICs) and tail vein injection to identify lung metastasis. RESULTS: Combined Pten deletion plus induction of p53-R270H mutation accelerated formation of four distinct mammary tumors including poorly differentiated adenocarcinoma (PDA) and spindle/mesenchymal-like lesions. Transplantation assays revealed highest frequency of TICs in PDA and spindle tumors compared with other subtypes. Hierarchical clustering demonstrated that the PDA and spindle tumors grouped closely with human as well as mouse models of basal and claudin-low subtypes, respectively. HTP screens of primary Pten(∆):p53(∆) vs. Pten(∆):p53(R270H) spindle tumor cells with 1120 FDA-approved drugs identified 8-azaguanine as most potent for both tumor types, but found no allele-specific inhibitor. A gene set enrichment analysis revealed increased expression of a metastasis pathway in Pten(∆):p53(R270H) vs. Pten(∆):p53(∆) spindle tumors. Accordingly, following tail vein injection, both Pten(∆):p53(R270H) spindle and PDA tumor cells induced lung metastases and morbidity significantly faster than Pten(∆):p53(∆) double-deletion cells, and this was associated with the ability of Pten(∆):p53(R270H) tumor cells to upregulate E-cadherin expression in lung metastases. CONCLUSIONS: Our results demonstrate that WAP-Cre:Pten(f/f):p53(lox.stop.lox_R270H) mice represent a tractable model to study basal-like breast cancer because unlike p53 deletion, p53(R270H) mutation in the mouse does not skew tumors toward the claudin-low subtype. The WAP-Cre:Pten(f/f):p53(lox.stop.lox_R270H) mice develop basal-like breast cancer that is enriched in TICs, can readily form lung metastasis, and provides a preclinical model to study both basal-like and claudin-low TNBC in immune-competent mice.

Seventeen-gene signature from enriched Her2/Neu mammary tumor-initiating cells predicts clinical outcome for human HER2+:ERα- breast cancer.

Human Epidermal Growth Factor Receptor 2-positive (HER2(+)) breast cancer (BC) is a highly aggressive disease commonly treated with chemotherapy and anti-HER2 drugs, including trastuzumab. There is currently no way to predict which HER2(+) BC patients will benefit from these treatments. Previous prognostic signatures for HER2(+) BC were developed irrespective of the subtype or the hierarchical organization of cancer in which only a fraction of cells, tumor-initiating cells (TICs), can sustain tumor growth. Here, we used serial dilution and single-cell transplantation assays to identify MMTV-Her2/Neu mouse mammary TICs as CD24(+):JAG1(-) at a frequency of 2-4.5%. A 17-gene Her2-TIC-enriched signature (HTICS), generated on the basis of differentially expressed genes in TIC versus non-TIC fractions and trained on one HER2(+) BC cohort, predicted clinical outcome on multiple independent HER2(+) cohorts. HTICS included up-regulated genes involved in S/G2/M transition and down-regulated genes involved in immune response. Its prognostic power was independent of other predictors, stratified lymph node(+) HER2(+) BC into low and high-risk subgroups, and was specific for HER2(+):estrogen receptor alpha-negative (ERα(-)) patients (10-y overall survival of 83.6% for HTICS(-) and 24.0% for HTICS(+) tumors; hazard ratio = 5.57; P = 0.002). Whereas HTICS was specific to HER2(+):ERα(-) tumors, a previously reported stroma-derived signature was predictive for HER2(+):ERα(+) BC. Retrospective analyses revealed that patients with HTICS(+) HER2(+):ERα(-) tumors resisted chemotherapy but responded to chemotherapy plus trastuzumab. HTICS is, therefore, a powerful prognostic signature for HER2(+):ERα(-) BC that can be used to identify high risk patients that would benefit from anti-HER2 therapy.

Meta-analysis of genome-wide scans for total body BMD in children and adults reveals allelic heterogeneity and age-specific effects at the WNT16 locus.

To identify genetic loci influencing bone accrual, we performed a genome-wide association scan for total-body bone mineral density (TB-BMD) variation in 2,660 children of different ethnicities. We discovered variants in 7q31.31 associated with BMD measurements, with the lowest P = 4.1 × 10(-11) observed for rs917727 with minor allele frequency of 0.37. We sought replication for all SNPs located ± 500 kb from rs917727 in 11,052 additional individuals from five independent studies including children and adults, together with de novo genotyping of rs3801387 (in perfect linkage disequilibrium (LD) with rs917727) in 1,014 mothers of children from the discovery cohort. The top signal mapping in the surroundings of WNT16 was replicated across studies with a meta-analysis P = 2.6 × 10(-31) and an effect size explaining between 0.6%-1.8% of TB-BMD variance. Conditional analyses on this signal revealed a secondary signal for total body BMD (P = 1.42 × 10(-10)) for rs4609139 and mapping to C7orf58. We also examined the genomic region for association with skull BMD to test if the associations were independent of skeletal loading. We identified two signals influencing skull BMD variation, including rs917727 (P = 1.9 × 10(-16)) and rs7801723 (P = 8.9 × 10(-28)), also mapping to C7orf58 (r(2) = 0.50 with rs4609139). Wnt16 knockout (KO) mice with reduced total body BMD and gene expression profiles in human bone biopsies support a role of C7orf58 and WNT16 on the BMD phenotypes observed at the human population level. In summary, we detected two independent signals influencing total body and skull BMD variation in children and adults, thus demonstrating the presence of allelic heterogeneity at the WNT16 locus. One of the skull BMD signals mapping to C7orf58 is mostly driven by children, suggesting temporal determination on peak bone mass acquisition. Our life-course approach postulates that these genetic effects influencing peak bone mass accrual may impact the risk of osteoporosis later in life.

WNT16 influences bone mineral density, cortical bone thickness, bone strength, and osteoporotic fracture risk.

We aimed to identify genetic variants associated with cortical bone thickness (CBT) and bone mineral density (BMD) by performing two separate genome-wide association study (GWAS) meta-analyses for CBT in 3 cohorts comprising 5,878 European subjects and for BMD in 5 cohorts comprising 5,672 individuals. We then assessed selected single-nucleotide polymorphisms (SNPs) for osteoporotic fracture in 2,023 cases and 3,740 controls. Association with CBT and forearm BMD was tested for ∼2.5 million SNPs in each cohort separately, and results were meta-analyzed using fixed effect meta-analysis. We identified a missense SNP (Thr>Ile; rs2707466) located in the WNT16 gene (7q31), associated with CBT (effect size of -0.11 standard deviations [SD] per C allele, P = 6.2 × 10(-9)). This SNP, as well as another nonsynonymous SNP rs2908004 (Gly>Arg), also had genome-wide significant association with forearm BMD (-0.14 SD per C allele, P = 2.3 × 10(-12), and -0.16 SD per G allele, P = 1.2 × 10(-15), respectively). Four genome-wide significant SNPs arising from BMD meta-analysis were tested for association with forearm fracture. SNP rs7776725 in FAM3C, a gene adjacent to WNT16, was associated with a genome-wide significant increased risk of forearm fracture (OR = 1.33, P = 7.3 × 10(-9)), with genome-wide suggestive signals from the two missense variants in WNT16 (rs2908004: OR = 1.22, P = 4.9 × 10(-6) and rs2707466: OR = 1.22, P = 7.2 × 10(-6)). We next generated a homozygous mouse with targeted disruption of Wnt16. Female Wnt16(-/-) mice had 27% (P<0.001) thinner cortical bones at the femur midshaft, and bone strength measures were reduced between 43%-61% (6.5 × 10(-13)

Prototype development of an electrical impedance based simultaneous respiratory and cardiac monitoring system for gated radiotherapy.

BACKGROUND: In radiotherapy, temporary translocations of the internal organs and tumor induced by respiratory and cardiac activities can undesirably lead to significantly lower radiation dose on the targeted tumor but more harmful radiation on surrounding healthy tissues. Respiratory and cardiac gated radiotherapy offers a potential solution for the treatment of tumors located in the upper thorax. The present study focuses on the design and development of simultaneous acquisition of respiratory and cardiac signal using electrical impedance technology for use in dual gated radiotherapy. METHODS: An electronic circuitry was developed for monitoring the bio-impedance change due to respiratory and cardiac motions and extracting the cardiogenic ECG signal. The system was analyzed in terms of reliability of signal acquisition, time delay, and functionality in a high energy radiation environment. The resulting signal of the system developed was also compared with the output of the commercially available Real-time Position Management™ (RPM) system in both time and frequency domains. RESULTS: The results demonstrate that the bioimpedance-based method can potentially provide reliable tracking of respiratory and cardiac motion in humans, alternative to currently available methods. When compared with the RPM system, the impedance-based system developed in the present study shows similar output pattern but different sensitivities in monitoring different respiratory rates. The tracking of cardiac motion was more susceptible to interference from other sources than respiratory motion but also provided synchronous output compared with the ECG signal extracted. The proposed hardware-based implementation was observed to have a worst-case time delay of approximately 33 ms for respiratory monitoring and 45 ms for cardiac monitoring. No significant effect on the functionality of the system was observed when it was tested in a radiation environment with the electrode lead wires directly exposed to high-energy X-Rays. CONCLUSION: The developed system capable of rendering quality signals for tracking both respiratory and cardiac motions can potentially provide a solution for simultaneous dual-gated radiotherapy.

Molecular classification to refine surgical and radiotherapeutic decision-making in meningioma.

Treatment of the tumor and dural margin with surgery and sometimes radiation are cornerstones of therapy for meningioma. Molecular classifications have provided insights into the biology of disease; however, response to treatment remains heterogeneous. In this study, we used retrospective data on 2,824 meningiomas, including molecular data on 1,686 tumors and 100 prospective meningiomas, from the RTOG-0539 phase 2 trial to define molecular biomarkers of treatment response. Using propensity score matching, we found that gross tumor resection was associated with longer progression-free survival (PFS) across all molecular groups and longer overall survival in proliferative meningiomas. Dural margin treatment (Simpson grade 1/2) prolonged PFS compared to no treatment (Simpson grade 3). Molecular group classification predicted response to radiotherapy, including in the RTOG-0539 cohort. We subsequently developed a molecular model to predict response to radiotherapy that discriminates outcome better than standard-of-care classification. This study highlights the potential for molecular profiling to refine surgical and radiotherapy decision-making.

Skeletal phenotypes in secreted frizzled-related protein 4 gene knockout mice mimic skeletal architectural abnormalities in subjects with Pyle's disease from SFRP4 mutations.

Mutations in SFRP4 cause Pyle's bone disease with wide metaphyses and increased skeletal fragility. The WNT signaling pathway plays important roles in determining skeletal architecture and SFRP4 is a secreted Frizzled decoy receptor that inhibits WNT signaling. Seven cohorts of male and female Sfrp4 gene knockout mice, examined through 2 years of age, had a normal lifespan but showed cortical and trabecular bone phenotypes. Mimicking human Erlenmeyer flask deformities, bone cross-sectional areas were elevated 2-fold in the distal femur and proximal tibia but only 30% in femur and tibia shafts. Reduced cortical bone thickness was observed in the vertebral body, midshaft femur and distal tibia. Elevated trabecular bone mass and numbers were observed in the vertebral body, distal femur metaphysis and proximal tibia metaphysis. Midshaft femurs retained extensive trabecular bone through 2 years of age. Vertebral bodies had increased compressive strength, but femur shafts had reduced bending strength. Trabecular, but not cortical, bone parameters in heterozygous Sfrp4 mice were modestly affected. Ovariectomy resulted in similar declines in both cortical and trabecular bone mass in wild-type and Sfrp4 KO mice. SFRP4 is critical for metaphyseal bone modeling involved in determining bone width. Sfrp4 KO mice show similar skeletal architecture and bone fragility deficits observed in patients with Pyle's disease with SFRP4 mutations.

Thinking (Metastasis) outside the (Primary Tumor) Box.

The metastasis of tumor cells into vital organs is a major cause of death from diverse types of malignancies [...].

Distinct shared and compartment-enriched oncogenic networks drive primary versus metastatic breast cancer.

Metastatic breast-cancer is a major cause of death in women worldwide, yet the relationship between oncogenic drivers that promote metastatic versus primary cancer is still contentious. To elucidate this relationship in treatment-naive animals, we hereby describe mammary-specific transposon-mutagenesis screens in female mice together with loss-of-function Rb, which is frequently inactivated in breast-cancer. We report gene-centric common insertion-sites (gCIS) that are enriched in primary-tumors, in metastases or shared by both compartments. Shared-gCIS comprise a major MET-RAS network, whereas metastasis-gCIS form three additional hubs: Rho-signaling, Ubiquitination and RNA-processing. Pathway analysis of four clinical cohorts with paired primary-tumors and metastases reveals similar organization in human breast-cancer with subtype-specific shared-drivers (e.g. RB1-loss, TP53-loss, high MET, RAS, ER), primary-enriched (EGFR, TGFß and STAT3) and metastasis-enriched (RHO, PI3K) oncogenic signaling. Inhibitors of RB1-deficiency or MET plus RHO-signaling cooperate to block cell migration and drive tumor cell-death. Thus, targeting shared- and metastasis- but not primary-enriched derivers offers a rational avenue to prevent metastatic breast-cancer.

Oncolytic DNX-2401 virotherapy plus pembrolizumab in recurrent glioblastoma: a phase 1/2 trial.

Immune-mediated anti-tumoral responses, elicited by oncolytic viruses and augmented with checkpoint inhibition, may be an effective treatment approach for glioblastoma. Here in this multicenter phase 1/2 study we evaluated the combination of intratumoral delivery of oncolytic virus DNX-2401 followed by intravenous anti-PD-1 antibody pembrolizumab in recurrent glioblastoma, first in a dose-escalation and then in a dose-expansion phase, in 49 patients. The primary endpoints were overall safety and objective response rate. The primary safety endpoint was met, whereas the primary efficacy endpoint was not met. There were no dose-limiting toxicities, and full dose combined treatment was well tolerated. The objective response rate was 10.4% (90% confidence interval (CI) 4.2-20.7%), which was not statistically greater than the prespecified control rate of 5%. The secondary endpoint of overall survival at 12 months was 52.7% (95% CI 40.1-69.2%), which was statistically greater than the prespecified control rate of 20%. Median overall survival was 12.5 months (10.7-13.5 months). Objective responses led to longer survival (hazard ratio 0.20, 95% CI 0.05-0.87). A total of 56.2% (95% CI 41.1-70.5%) of patients had a clinical benefit defined as stable disease or better. Three patients completed treatment with durable responses and remain alive at 45, 48 and 60 months. Exploratory mutational, gene-expression and immunophenotypic analyses revealed that the balance between immune cell infiltration and expression of checkpoint inhibitors may potentially inform on response to treatment and mechanisms of resistance. Overall, the combination of intratumoral DNX-2401 followed by pembrolizumab was safe with notable survival benefit in select patients (ClinicalTrials.gov registration: NCT02798406).

Single-Cell, Single-Nucleus, and Spatial RNA Sequencing of the Human Liver Identifies Cholangiocyte and Mesenchymal Heterogeneity.

The critical functions of the human liver are coordinated through the interactions of hepatic parenchymal and non-parenchymal cells. Recent advances in single-cell transcriptional approaches have enabled an examination of the human liver with unprecedented resolution. However, dissociation-related cell perturbation can limit the ability to fully capture the human liver's parenchymal cell fraction, which limits the ability to comprehensively profile this organ. Here, we report the transcriptional landscape of 73,295 cells from the human liver using matched single-cell RNA sequencing (scRNA-seq) and single-nucleus RNA sequencing (snRNA-seq). The addition of snRNA-seq enabled the characterization of interzonal hepatocytes at a single-cell resolution, revealed the presence of rare subtypes of liver mesenchymal cells, and facilitated the detection of cholangiocyte progenitors that had only been observed during in vitro differentiation experiments. However, T and B lymphocytes and natural killer cells were only distinguishable using scRNA-seq, highlighting the importance of applying both technologies to obtain a complete map of tissue-resident cell types. We validated the distinct spatial distribution of the hepatocyte, cholangiocyte, and mesenchymal cell populations by an independent spatial transcriptomics data set and immunohistochemistry. Conclusion: Our study provides a systematic comparison of the transcriptomes captured by scRNA-seq and snRNA-seq and delivers a high-resolution map of the parenchymal cell populations in the healthy human liver.

Loss of Epigenetic Regulation Disrupts Lineage Integrity, Induces Aberrant Alveogenesis, and Promotes Breast Cancer.

Systematically investigating the scores of genes mutated in cancer and discerning disease drivers from inconsequential bystanders is a prerequisite for precision medicine but remains challenging. Here, we developed a somatic CRISPR/Cas9 mutagenesis screen to study 215 recurrent "long-tail" breast cancer genes, which revealed epigenetic regulation as a major tumor-suppressive mechanism. We report that components of the BAP1 and COMPASS-like complexes, including KMT2C/D, KDM6A, BAP1, and ASXL1/2 ("EpiDrivers"), cooperate with PIK3CAH1047R to transform mouse and human breast epithelial cells. Mechanistically, we find that activation of PIK3CAH1047R and concomitant EpiDriver loss triggered an alveolar-like lineage conversion of basal mammary epithelial cells and accelerated formation of luminal-like tumors, suggesting a basal origin for luminal tumors. EpiDriver mutations are found in ∼39% of human breast cancers, and ∼50% of ductal carcinoma in situ express casein, suggesting that lineage infidelity and alveogenic mimicry may significantly contribute to early steps of breast cancer etiology. SIGNIFICANCE: Infrequently mutated genes comprise most of the mutational burden in breast tumors but are poorly understood. In vivo CRISPR screening identified functional tumor suppressors that converged on epigenetic regulation. Loss of epigenetic regulators accelerated tumorigenesis and revealed lineage infidelity and aberrant expression of alveogenesis genes as potential early events in tumorigenesis. This article is highlighted in the In This Issue feature, p. 2711.

Ascorbate synthesis pathway: dual role of ascorbate in bone homeostasis.

Using mouse gene knock-out models, we identify aldehyde reductase (EC 1.1.1.2, Akr1a4 (GR)) and aldose reductase (EC 1.1.1.21, Akr1b3 (AR)) as the enzymes responsible for conversion of D-glucuronate to L-gulonate, a key step in the ascorbate (ASC) synthesis pathway in mice. The gene knock-out (KO) mice show that the two enzymes, GR and AR, provide approximately 85 and approximately 15% of L-gulonate, respectively. GRKO/ARKO double knock-out mice are unable to synthesize ASC (>95% ASC deficit) and develop scurvy. The GRKO mice ( approximately 85% ASC deficit) develop and grow normally when fed regular mouse chow (ASC content = 0) but suffer severe osteopenia and spontaneous fractures with stresses that increase ASC requirements, such as pregnancy or castration. Castration greatly increases osteoclast numbers and activity in GRKO mice and promotes increased bone loss as compared with wild-type controls and additionally induces proliferation of immature dysplastic osteoblasts likely because of an ASC-sensitive block(s) in early differentiation. ASC and the antioxidants pycnogenol and resveratrol block osteoclast proliferation and bone loss, but only ASC feeding restores osteoblast differentiation and prevents their dysplastic proliferation. This is the first in vivo demonstration of two independent roles for ASC as an antioxidant suppressing osteoclast activity and number as well as a cofactor promoting osteoblast differentiation. Although humans have lost the ability to synthesize ASC, our mouse models suggest the mechanisms by which suboptimal ASC availability facilitates the development of osteoporosis, which has important implications for human osteoporosis.