Search

Triple negative breast cancer (TNBC) is the most aggressive breast cancer subtype and has the highest rate of recurrence. The predominant standard of care for advanced TNBC is systemic chemotherapy with or without immunotherapy, however responses are typically short-lived. Thus, there is an urgent need to develop more effective treatments. PI3K pathway components represent plausible therapeutic targets, as approximately 40% of TNBCs have PIK3CA/AKT1/PTEN alterations. However, unlike hormone receptor-positive tumors, it is still unclear if or how PI3K pathway inhibitors will be effective in triple-negative disease. Here we identify a promising AKT inhibitor-based therapeutic combination for TNBC. Specifically, we show that AKT inhibitors potently synergize with agents that suppress the histone methyltransferase, EZH2, and promote robust tumor regression in multiple TNBC models in vivo. AKT and EZH2 inhibitors exert these effects by first cooperatively driving basal-like TNBC cells into a more differentiated, luminal-like state, which cannot be effectively induced by either agent alone. More importantly, once differentiated, these agents kill TNBCs by hijacking signals that normally drive mammary gland involution. Together these findings identify a promising therapeutic strategy for this highly aggressive tumor type and illustrate how deregulated epigenetic enzymes can insulate tumors from oncogenic vulnerabilities. These studies also reveal how developmental tissue-specific cell death pathways may be co-opted for therapeutic benefit. Citation Format: Amy Schade, Naiara Perurena, Marina Watanabe, Carrie L. Rodriguez, Patrick Loi, Natalie Pilla, Rachel A. Davis, Kaia Mattioli, Dongxi Xiang, Jason J. Zoeller, Zhe Li, Ana C. Garrido-Castro, Sara Tolaney, Karen Cichowski. AKT and EZH2 inhibitors kill TNBCs by hijacking mechanisms of involution [abstract]. In: Proceedings of the 2022 San Antonio Breast Cancer Symposium; 2022 Dec 6-10; San Antonio, TX. Philadelphia (PA): AACR; Cancer Res 2023;83(5 Suppl):Abstract nr P4-08-11.

High blood pressure (BP) is the major risk factor for cardiovascular disease. Genome-wide association studies have identified genetic variants for BP, but functional insights into causality and related molecular mechanisms lag behind. We functionally characterize 4,608 genetic variants in linkage with 135 BP loci in vascular smooth muscle cells and cardiomyocytes by massively parallel reporter assays. High densities of regulatory variants at BP loci (i.e., ULK4, MAP4, CFDP1, PDE5A) indicate that multiple variants drive genetic association. Regulatory variants are enriched in repeats, alter cardiovascular-related transcription factor motifs, and spatially converge with genes controlling specific cardiovascular pathways. Using heuristic scoring, we define likely causal variants, and CRISPR prime editing finally determines causal variants for KCNK9, SFXN2, and PCGF6, which are candidates for developing high BP. Our systems-level approach provides a catalog of functionally relevant variants and their genomic architecture in two trait-relevant cell lines for a better understanding of BP gene regulation. We thank the Melé and Maass labs for intellectual input, Dr. Steven Erwood from Dr. Ronald Cohn’s lab for guidance in applying CRISPR prime editing, The Centre for Applied Genomics, The Structural & Biophysical Core Facility, and The Imaging Facility, The Hospital for Sick Children, Toronto, Canada, for assistance with high-throughput sequencing, luminescence detection, and imaging. We thank Dovetail Genomics, LLC, 100 Enterprise Way, Suite A101, Scotts Valley, CA 95066, USA, for generating Omni-C libraries and for the collaborative support throughout the project. W.O. was supported by a Fundació la Marató grant (ref. 321/C/2019), K.K. was supported by an OGS fellowship, J.W.L.B. was supported by a CGS-D fellowship, and D.F.L. was supported by an Ontario Genomics-CANSSI Ontario Postdoctoral Fellowship in Genome Data Science. This project was supported by Canada’s New Frontiers in Research Fund (NFRFE-2018-01305), the Canadian Institutes of Health Research (CIHR PJT 173542 [P.G.M.] and PJT 175034 [S.M., J.E.]), CIHR ENP 161429 under the frame of ERA PerMed (S.M.), the Ted Rogers Centre for Heart Research (S.M., J.E.), and the Heart and Stroke Foundation of Canada. J.E. holds a Canada Research Chair Tier 1 in Stem Cell Models of Childhood Disease, S.M. holds the Heart and Stroke Foundation of Canada & Robert M. Freedom Chair in Cardiovascular Science, M. Melé was supported by a Ramon y Cajal fellowship (RYC-2017-22249), and P.G.M. holds a Canada Research Chair Tier 2 in Non-coding Disease Mechanisms. Peer Reviewed "Article signat per 14 autors/es: Winona Oliveros, Kate Delfosse, Daniella F. Lato , Katerina Kiriakopulos, Milad Mokhtaridoost, Abdelrahman Said, Brandon J. McMurray, Jared W.L. Browning, Kaia Mattioli, Guoliang Meng, James Ellis, Seema Mital, Marta Melé, Philipp G. Maass"

Cardiomyopathy (CMP) is a heritable disorder. Over 50% of cases are gene-elusive on clinical gene panel testing. The contribution of variants in non-coding DNA elements that result in cryptic splicing and regulate gene expression has not been explored. We analyzed whole-genome sequencing (WGS) data in a discovery cohort of 209 pediatric CMP patients and 1953 independent replication genomes and exomes. We searched for protein-coding variants, and non-coding variants predicted to affect the function or expression of genes. Thirty-nine percent of cases harbored pathogenic coding variants in known CMP genes, and 5% harbored high-risk loss-of-function (LoF) variants in additional candidate CMP genes. Fifteen percent harbored high-risk regulatory variants in promoters and enhancers of CMP genes (odds ratio 2.25, p = 6.70 × 10−7 versus controls). Genes involved in α-dystroglycan glycosylation (FKTN, DTNA) and desmosomal signaling (DSC2, DSG2) were most highly enriched for regulatory variants (odds ratio 6.7–58.1). Functional effects were confirmed in patient myocardium and reporter assays in human cardiomyocytes, and in zebrafish CRISPR knockouts. We provide strong evidence for the genomic contribution of functionally active variants in new genes and in regulatory elements of known CMP genes to early onset CMP. This project was supported by the Ted Rogers Centre for Heart Research (SM, JE), the Canadian Institutes of Health Research (PJT 175034) (SM, JE) and by the Canadian Institutes of Health Research (ENP 161429), under the frame of ERA PerMed (SM). SM holds the Heart and Stroke Foundation of Canada & Robert M Freedom Chair in Cardiovascular Science. SWS holds the GlaxoSmithKline Endowed Chair in Genome Sciences at the Hospital for Sick Children and the University of Toronto. PGM holds a Canada Research Chair Tier 2 in Non-coding Disease Mechanisms. PGM acknowledges the support of the Government of Canada’s New Frontiers in Research Fund (NFRF), [NFRFE-2018-01305]. EO holds the Bitove Family Professorship of Adult Congenital Heart Disease. MM holds a Ramon y Cajal grant from the Spanish Ministry of Science and Innovation (RYC-2017-22249). WO is supported by funding from Fundació La Marató (321/C/2019). JB is funded by a Frans Van de Werf fellowship for clinical cardiovascular research, and by a senior clinical investigator fellowship of the FWO Flanders. KM was a National Science Foundation Graduate Research Fellow under grant no. DGE1144152 during the majority of the project. CS is the recipient of a National Health and Medical Research Council (NHMRC) Practitioner Fellowship (1154992). JI is the recipient of an NHMRC Career Development Fellowship (1162929). RDB is the recipient of a New South Wales Health Cardiovascular Disease Senior Scientist Grant. PSD is supported by the DBT/Wellcome Trust- Indian Alliance. We acknowledge the Labatt Family Heart Centre Biobank at the Hospital for Sick Children for access to DNA samples, and The Centre for Applied Genomics at the Hospital for Sick Children for performing WGS. We thank Xiucheng Cui and Emanuela Pannia for performing the zebrafish experiments at the SickKids Zebrafish Genetics and Disease Models Core (CRISPR-Cas9 and gRNA syntheses, zebrafish embryo microinjections, gRNA PCR validation, qRT-PCR, cardiac imaging). This research was made possible through access to the data and findings generated by the 100,000 Genomes Project. The 100,000 Genomes Project is managed by Genomics England Limited (a wholly owned company of the Department of Health and Social Care). The 100,000 Genomes Project is funded by the National Institute for Health Research and NHS England. The Wellcome Trust, Cancer Research UK and the Medical Research Council have also funded research infrastructure. The 100,000 Genomes Project uses data provided by patients and collected by the National Health Service as part of their care and support. We thank members of the ICGC/PCAWG working groups for generating the variant calls used in our case-control burden analyses. Peer Reviewed "Article signat per 38 autors/es: Robert Lesurf, Abdelrahman Said, Oyediran Akinrinade, Jeroen Breckpot, Kathleen Delfosse, Ting Liu, Roderick Yao, Gabrielle Persad, Fintan McKenna, Ramil R. Noche, Winona Oliveros, Kaia Mattioli, Shreya Shah, Anastasia Miron, Qian Yang, Guoliang Meng, Michelle Chan Seng Yue, Wilson W. L. Sung, Bhooma Thiruvahindrapuram, Jane Lougheed, Erwin Oechslin, Tapas Mondal, Lynn Bergin, John Smythe, Shashank Jayappa, Vinay J. Rao, Jayaprakash Shenthar, Perundurai S. Dhandapany, Christopher Semsarian, Robert G. Weintraub, Richard D. Bagnall, Jodie Ingles, Genomics England Research Consortium, Marta Melé, Philipp G. Maass, James Ellis, Stephen W. Scherer & Seema Mital"

Although many long non-coding RNAs (lncRNAs) exhibit lineage-specific expression, the vast majority remain functionally uncharacterized in the context of development. Here, we report the first described human embryonic stem cell (hESC) lines to repress (CRISPRi) or activate (CRISPRa) transcription during differentiation into all three germ layers, facilitating the modulation of lncRNA expression during early development. We performed an unbiased, genome-wide CRISPRi screen targeting thousands of lncRNA loci expressed during endoderm differentiation. While dozens of lncRNA loci were required for proper differentiation, most differentially expressed lncRNAs were not, supporting the necessity for functional screening instead of relying solely on gene expression analyses. In parallel, we developed a clustering approach to infer mechanisms of action of lncRNA hits based on a variety of genomic features. We subsequently identified and validated FOXD3-AS1 as a functional lncRNA essential for pluripotency and differentiation. Taken together, the cell lines and methodology described herein can be adapted to discover and characterize novel regulators of differentiation into any lineage.

Although many long non-coding RNAs (lncRNAs) exhibit lineage-specific expression, the vast majority remain functionally uncharacterized in the context of development. Here, we report the first described human embryonic stem cell (hESC) lines to repress (CRISPRi) or activate (CRISPRa) transcription during differentiation into all three germ layers, facilitating the modulation of lncRNA expression during early development. We performed an unbiased, genome-wide CRISPRi screen targeting thousands of lncRNA loci expressed during endoderm differentiation. While dozens of lncRNA loci were required for proper differentiation, most differentially expressed lncRNAs were not, supporting the necessity for functional screening instead of relying solely on gene expression analyses. In parallel, we developed a clustering approach to infer mechanisms of action of lncRNA hits based on a variety of genomic features. We subsequently identified and validated FOXD3-AS1 as a functional lncRNA essential for pluripotency and differentiation. Taken together, the cell lines and methodology described herein can be adapted to discover and characterize novel regulators of differentiation into any lineage.

Cardiomyopathy (CMP) is a heritable genetic disorder. Protein-coding variants account for 20-30% of cases. The contribution of variants in non-coding DNA elements that regulate gene expression has not been explored. We performed whole-genome sequencing (WGS) of 228 unrelated CMP families. Besides pathogenic protein-coding variants in known CMP genes, 5% cases harbored rare loss-of-function variants in novel cardiac genes, with NRAP and FHOD3 being strong candidates. WGS also revealed a high burden of high-risk variants in promoters and enhancers of CMP genes in an additional 20% cases (Odds ratio 2.14, 95% CI 1.60-2.86, p=5.26×10−7 vs 1326 controls) with genes involved in α-dystroglycan glycosylation (FKTN, DTNA) and desmosomal signaling (DSC2, DSG2) specifically enriched for regulatory variants (False discovery rate

Deriving mechanisms of immune‐mediated disease from GWAS data remains a formidable challenge, with attempts to identify causal variants being frequently hampered by strong linkage disequilibrium. To determine whether causal variants could be identified from their functional effects, we adapted a massively parallel reporter assay for use in primary CD4 T cells, the cell type whose regulatory DNA is most enriched for immune‐mediated disease SNPs. This enabled the effects of candidate SNPs to be examined in a relevant cellular context and generated testable hypotheses into disease mechanisms. To illustrate the power of this approach, we investigated a locus that has been linked to six immune‐mediated diseases but cannot be fine‐mapped. By studying the lead expression‐modulating SNP, we uncovered an NF‐κB‐driven regulatory circuit which constrains T‐cell activation through the dynamic formation of a super‐enhancer that upregulates TNFAIP3 (A20), a key NF‐κB inhibitor. In activated T cells, this feedback circuit is disrupted—and super‐enhancer formation prevented—by the risk variant at the lead SNP, leading to unrestrained T‐cell activation via a molecular mechanism that appears to broadly predispose to human autoimmunity., Little progress has been made in resolving causal SNPs, genes and disease mechanisms at GWAS loci. An adapted massively‐parallel reporter assay (MPRA) allows to study immune‐mediated disease loci in CD4 T cells, the cell‐type whose regulatory DNA is most highly enriched for disease‐associated SNPs.

Deriving mechanisms of immune-mediated disease from GWAS data remains a formidable challenge, with attempts to identify causal variants being frequently hampered by strong linkage disequilibrium. To determine whether causal variants could be identified from their functional effects, we adapted a massively parallel reporter assay for use in primary CD4 T cells, the cell type whose regulatory DNA is most enriched for immune-mediated disease SNPs. This enabled the effects of candidate SNPs to be examined in a relevant cellular context and generated testable hypotheses into disease mechanisms. To illustrate the power of this approach, we investigated a locus that has been linked to six immune-mediated diseases but cannot be fine-mapped. By studying the lead expression-modulating SNP, we uncovered an NF-κB-driven regulatory circuit which constrains T-cell activation through the dynamic formation of a super-enhancer that upregulates TNFAIP3 (A20), a key NF-κB inhibitor. In activated T cells, this feedback circuit is disrupted-and super-enhancer formation prevented-by the risk variant at the lead SNP, leading to unrestrained T-cell activation via a molecular mechanism that appears to broadly predispose to human autoimmunity.

Background Gene expression differences between species are driven by both cis and trans effects. Whereas cis effects are caused by genetic variants located on the same DNA molecule as the target gene, trans effects are due to genetic variants that affect diffusible elements. Previous studies have mostly assessed the impact of cis and trans effects at the gene level. However, how cis and trans effects differentially impact regulatory elements such as enhancers and promoters remains poorly understood. Here, we use massively parallel reporter assays to directly measure the transcriptional outputs of thousands of individual regulatory elements in embryonic stem cells and measure cis and trans effects between human and mouse. Results Our approach reveals that cis effects are widespread across transcribed regulatory elements, and the strongest cis effects are associated with the disruption of motifs recognized by strong transcriptional activators. Conversely, we find that trans effects are rare but stronger in enhancers than promoters and are associated with a subset of transcription factors that are differentially expressed between human and mouse. While we find that cis-trans compensation is common within promoters, we do not see evidence of widespread cis-trans compensation at enhancers. Cis-trans compensation is inversely correlated with enhancer redundancy, suggesting that such compensation may often occur across multiple enhancers. Conclusions Our results highlight differences in the mode of evolution between promoters and enhancers in complex mammalian genomes and indicate that studying the evolution of individual regulatory elements is pivotal to understand the tempo and mode of gene expression evolution.

Gene expression differences between species are driven by bothcisandtranseffects. Whereasciseffects are caused by genetic variants in close proximity to the target gene,transeffects are due to distal genetic variants that affect diffusible elements such as transcription factors. Previous studies have mostly assessed the impact ofcisandtranseffects at the gene level. However, howcisandtranseffects differentially impact regulatory elements such as enhancers and promoters remains poorly understood. Here, we used massively parallel reporter assays to directly measurecisandtranseffects between human and mouse embryonic stem cells at thousands of individual regulatory elements. Our approach revealed thatciseffects are widespread across regulatory elements, and the strongestciseffects are associated with the disruption of motifs recognized by strong transcriptional activators. Conversely, we found thattranseffects are rare but stronger in enhancers than promoters, and can be attributed to a subset of transcription factors that are differentially expressed between human and mouse. While previous studies have found extensive co-occurrence ofcisandtranseffects in opposite directions that stabilize gene expression throughout evolution, we find thatcis-transcompensation is uncommon within individual regulatory elements. Thus, our results are consistent with a model wherein compensatorycis-transeffects at the gene level are explained bycisandtranseffects that separately impact several regulatory elements rather thancis-transeffects that occur simultaneously within a single regulatory element. Together, these results indicate that studying the evolution of individual regulatory elements is pivotal to understand the tempo and mode of gene expression evolution.

Numerous genomic loci have been implicated in autoimmune disorders, but attempts to identify causal variants, and thereby disease mechanisms, have been hampered by strong linkage disequilibrium; leaving most loci unresolved and the potential of GWAS largely unfulfilled. To overcome this, we developed a massively-parallel reporter assay for use in primary CD4 T-cells, a key effector of many autoimmune diseases, and sought to resolve potential causal variants via their functional effects. This provided testable hypotheses into disease mechanisms, which we illustrate using a multi-disease-associated locus on chromosome 6. By investigating the lead expression-modulating SNP, we uncover an NF-κB-driven regulatory circuit which constrains T-cell activation through the dynamic formation of a super-enhancer that upregulates expression of TNFAIP3 (A20), a key NF-κB inhibitor. In activated T-cells, this circuit is disrupted – and super-enhancer formation prevented – by the risk variant at the lead SNP, thus revealing a molecular mechanism that predisposes to multiple autoimmune diseases.

Transcription initiates at both coding and non-coding genomic elements, including mRNA and long non-coding RNA (lncRNA) core promoters and enhancer RNAs (eRNAs). However, each class has different expression profiles with lncRNAs and eRNAs being the most tissue-specific. How these complex differences in expression profiles and tissue-specificities are encoded in a single DNA sequence, however, remains unresolved. Here, we address this question using computational approaches and massively parallel reporter assays (MPRA) surveying hundreds of promoters and enhancers. We find that both divergent lncRNA and mRNA core promoters have higher capacities to drive transcription than non-divergent lncRNA and mRNA core promoters, respectively. Conversely, lincRNAs and eRNAs have lower capacities to drive transcription and are more tissue-specific than divergent genes. This higher tissue-specificity is strongly associated with having less complex TF motif profiles at the core promoter. We experimentally validated these findings by testing both engineered single-nucleotide deletions and human single-nucleotide polymorphisms (SNPs) in MPRA. In both cases, we observe that single nucleotides associated with many motifs are important drivers of promoter activity. Thus, we suggest that high TF motif density serves as a robust mechanism to increase promoter activity at the expense of tissue-specificity. Moreover, we find that 22% of common SNPs in core promoter regions have significant regulatory effects. Collectively, our findings show that high TF motif density provides redundancy and increases promoter activity at the expense of tissue specificity, suggesting that specificity of expression may be regulated by simplicity of motif usage.

Transcription initiates at both coding and noncoding genomic elements, including mRNA and long noncoding RNA (lncRNA) core promoters and enhancer RNAs (eRNAs). However, each class has a different expression profile with lncRNAs and eRNAs being the most tissue specific. How these complex differences in expression profiles and tissue specificities are encoded in a single DNA sequence remains unresolved. Here, we address this question using computational approaches and massively parallel reporter assays (MPRA) surveying hundreds of promoters and enhancers. We find that both divergent lncRNA and mRNA core promoters have higher capacities to drive transcription than nondivergent lncRNA and mRNA core promoters, respectively. Conversely, intergenic lncRNAs (lincRNAs) and eRNAs have lower capacities to drive transcription and are more tissue specific than divergent genes. This higher tissue specificity is strongly associated with having less complex transcription factor (TF) motif profiles at the core promoter. We experimentally validated these findings by testing both engineered single-nucleotide deletions and human single-nucleotide polymorphisms (SNPs) in MPRA. In both cases, we observe that single nucleotides associated with many motifs are important drivers of promoter activity. Thus, we suggest that high TF motif density serves as a robust mechanism to increase promoter activity at the expense of tissue specificity. Moreover, we find that 22% of common SNPs in core promoter regions have significant regulatory effects. Collectively, our findings show that high TF motif density provides redundancy and increases promoter activity at the expense of tissue specificity, suggesting that specificity of expression may be regulated by simplicity of motif usage.