CRISPR-Cas13d screens identify KILR, a breast cancer risk-associated lncRNA that regulates DNA replication and repair.

BACKGROUND: Long noncoding RNAs (lncRNAs) have surpassed the number of protein-coding genes, yet the majority have no known function. We previously discovered 844 lncRNAs that were genetically linked to breast cancer through genome-wide association studies (GWAS). Here, we show that a subset of these lncRNAs alter breast cancer risk by modulating cell proliferation, and provide evidence that a reduced expression on one lncRNA increases breast cancer risk through aberrant DNA replication and repair. METHODS: We performed pooled CRISPR-Cas13d-based knockdown screens in breast cells to identify which of the 844 breast cancer-associated lncRNAs alter cell proliferation. We selected one of the lncRNAs that increased cell proliferation, KILR, for follow-up functional studies. KILR pull-down followed by mass spectrometry was used to identify binding proteins. Knockdown and overexpression studies were performed to assess the mechanism by which KILR regulates proliferation. RESULTS: We show that KILR functions as a tumor suppressor, safeguarding breast cells against uncontrolled proliferation. The half-life of KILR is significantly reduced by the risk haplotype, revealing an alternative mechanism by which variants alter cancer risk. Mechanistically, KILR sequesters RPA1, a subunit of the RPA complex required for DNA replication and repair. Reduced KILR expression promotes breast cancer cell proliferation by increasing the available pool of RPA1 and speed of DNA replication. Conversely, KILR overexpression promotes apoptosis in breast cancer cells, but not normal breast cells. CONCLUSIONS: Our results confirm lncRNAs as mediators of breast cancer risk, emphasize the need to annotate noncoding transcripts in relevant cell types when investigating GWAS variants and provide a scalable platform for mapping phenotypes associated with lncRNAs.

Preventing recurrence in Sonic Hedgehog Subgroup Medulloblastoma using the OLIG2 inhibitor CT-179.

Recurrence is the primary life-threatening complication for medulloblastoma (MB). In Sonic Hedgehog (SHH)-subgroup MB, OLIG2-expressing tumor stem cells drive recurrence. We investigated the anti-tumor potential of the small-molecule OLIG2 inhibitor CT-179, using SHH-MB patient-derived organoids, patient-derived xenograft (PDX) tumors and mice genetically-engineered to develop SHH-MB. CT-179 disrupted OLIG2 dimerization, DNA binding and phosphorylation and altered tumor cell cycle kinetics in vitro and in vivo, increasing differentiation and apoptosis. CT-179 increased survival time in GEMM and PDX models of SHH-MB, and potentiated radiotherapy in both organoid and mouse models, delaying post-radiation recurrence. Single cell transcriptomic studies (scRNA-seq) confirmed that CT-179 increased differentiation and showed that tumors up-regulated Cdk4 post-treatment. Consistent with increased CDK4 mediating CT-179 resistance, CT-179 combined with CDK4/6 inhibitor palbociclib delayed recurrence compared to either single-agent. These data show that targeting treatment-resistant MB stem cell populations by adding the OLIG2 inhibitor CT-179 to initial MB treatment can reduce recurrence.

Redefining normal breast cell populations using long noncoding RNAs.

Single-cell RNAseq has allowed unprecedented insight into gene expression across different cell populations in normal tissue and disease states. However, almost all studies rely on annotated gene sets to capture gene expression levels and sequencing reads that do not align to known genes are discarded. Here, we discover thousands of long noncoding RNAs (lncRNAs) expressed in human mammary epithelial cells and analyze their expression in individual cells of the normal breast. We show that lncRNA expression alone can discriminate between luminal and basal cell types and define subpopulations of both compartments. Clustering cells based on lncRNA expression identified additional basal subpopulations, compared to clustering based on annotated gene expression, suggesting that lncRNAs can provide an additional layer of information to better distinguish breast cell subpopulations. In contrast, these breast-specific lncRNAs poorly distinguish brain cell populations, highlighting the need to annotate tissue-specific lncRNAs prior to expression analyses. We also identified a panel of 100 breast lncRNAs that could discern breast cancer subtypes better than protein-coding markers. Overall, our results suggest that lncRNAs are an unexplored resource for new biomarker and therapeutic target discovery in the normal breast and breast cancer subtypes.

CRISPR screens identify gene targets at breast cancer risk loci.

BACKGROUND: Genome-wide association studies (GWAS) have identified > 200 loci associated with breast cancer risk. The majority of candidate causal variants are in non-coding regions and likely modulate cancer risk by regulating gene expression. However, pinpointing the exact target of the association, and identifying the phenotype it mediates, is a major challenge in the interpretation and translation of GWAS. RESULTS: Here, we show that pooled CRISPR screens are highly effective at identifying GWAS target genes and defining the cancer phenotypes they mediate. Following CRISPR mediated gene activation or suppression, we measure proliferation in 2D, 3D, and in immune-deficient mice, as well as the effect on DNA repair. We perform 60 CRISPR screens and identify 20 genes predicted with high confidence to be GWAS targets that promote cancer by driving proliferation or modulating the DNA damage response in breast cells. We validate the regulation of a subset of these genes by breast cancer risk variants. CONCLUSIONS: We demonstrate that phenotypic CRISPR screens can accurately pinpoint the gene target of a risk locus. In addition to defining gene targets of risk loci associated with increased breast cancer risk, we provide a platform for identifying gene targets and phenotypes mediated by risk variants.

C16orf72/HAPSTR1 is a molecular rheostat in an integrated network of stress response pathways.

All cells contain specialized signaling pathways that enable adaptation to specific molecular stressors. Yet, whether these pathways are centrally regulated in complex physiological stress states remains unclear. Using genome-scale fitness screening data, we quantified the stress phenotype of 739 cancer cell lines, each representing a unique combination of intrinsic tumor stresses. Integrating dependency and stress perturbation transcriptomic data, we illuminated a network of genes with vital functions spanning diverse stress contexts. Analyses for central regulators of this network nominated C16orf72/HAPSTR1, an evolutionarily ancient gene critical for the fitness of cells reliant on multiple stress response pathways. We found that HAPSTR1 plays a pleiotropic role in cellular stress signaling, functioning to titrate various specialized cell-autonomous and paracrine stress response programs. This function, while dispensable to unstressed cells and nematodes, is essential for resilience in the presence of stressors ranging from DNA damage to starvation and proteotoxicity. Mechanistically, diverse stresses induce HAPSTR1, which encodes a protein expressed as two equally abundant isoforms. Perfectly conserved residues in a domain shared between HAPSTR1 isoforms mediate oligomerization and binding to the ubiquitin ligase HUWE1. We show that HUWE1 is a required cofactor for HAPSTR1 to control stress signaling and that, in turn, HUWE1 feeds back to ubiquitinate and destabilize HAPSTR1. Altogether, we propose that HAPSTR1 is a central rheostat in a network of pathways responsible for cellular adaptability, the modulation of which may have broad utility in human disease.

Auxin-Inducible Degron System Reveals Temporal-Spatial Roles of HSF-1 and Its Transcriptional Program in Lifespan Assurance.

HSF-1 is a key regulator of cellular proteotoxic stress response and is required for animal lifespan. In C. elegans, HSF-1 mediated heat shock response (HSR) declines sharply on the first day of adulthood, and HSF-1 was proposed to function primarily during larval stages for lifespan assurance based on studies using RNAi. The tissue requirement for HSF-1 in lifespan, however, is not well understood. Using the auxin-inducible degron (AID) system, we manage to uncouple the roles of HSF-1 in development and longevity. In wild-type animals, we find HSF-1 is required during the whole self-reproductive period for lifespan. This period is extended in long-lived animals that have arrested germline stem cells (GSC) or reduced insulin/IGF-1 signaling (IIS). While depletion of HSF-1 from any major somatic tissues during development results in severe defects, HSF-1 primarily functions in the intestine and likely neural system of adults to support lifespan. Finally, by combining AID and genome-wide transcriptional analyses, we find HSF-1 directly activates the transcription of constitutively-expressed chaperone and co-chaperone genes among others in early adulthood, which underlies its roles in longevity assurance.

Insulin/IGF-1 signaling and heat stress differentially regulate HSF1 activities in germline development.

Germline development is sensitive to nutrient availability and environmental perturbation. Heat shock transcription factor 1 (HSF1), a key transcription factor driving the cellular heat shock response (HSR), is also involved in gametogenesis. The precise function of HSF1 (HSF-1 in C. elegans) and its regulation in germline development are poorly understood. Using the auxin-inducible degron system in C. elegans, we uncovered a role of HSF-1 in progenitor cell proliferation and early meiosis and identified a compact but important transcriptional program of HSF-1 in germline development. Interestingly, heat stress only induces the canonical HSR in a subset of germ cells but impairs HSF-1 binding at its developmental targets. Conversely, insulin/insulin growth factor 1 (IGF-1) signaling dictates the requirement for HSF-1 in germline development and functions through repressing FOXO/DAF-16 in the soma to activate HSF-1 in germ cells. We propose that this non-cell-autonomous mechanism couples nutrient-sensing insulin/IGF-1 signaling to HSF-1 activation to support homeostasis in rapid germline growth.

eQTL Colocalization Analyses Identify NTN4 as a Candidate Breast Cancer Risk Gene.

Breast cancer genome-wide association studies (GWASs) have identified 150 genomic risk regions containing more than 13,000 credible causal variants (CCVs). The CCVs are predominantly noncoding and enriched in regulatory elements. However, the genes underlying breast cancer risk associations are largely unknown. Here, we used genetic colocalization analysis to identify loci at which gene expression could potentially explain breast cancer risk phenotypes. Using data from the Breast Cancer Association Consortium (BCAC) and quantitative trait loci (QTL) from the Genotype-Tissue Expression (GTEx) project and The Cancer Genome Project (TCGA), we identify shared genetic relationships and reveal novel associations between cancer phenotypes and effector genes. Seventeen genes, including NTN4, were identified as potential mediators of breast cancer risk. For NTN4, we showed the rs61938093 CCV at this region was located within an enhancer element that physically interacts with the NTN4 promoter, and the risk allele reduced NTN4 promoter activity. Furthermore, knockdown of NTN4 in breast cells increased cell proliferation in vitro and tumor growth in vivo. These data provide evidence linking risk-associated variation to genes that may contribute to breast cancer predisposition.

Chromatin interactome mapping at 139 independent breast cancer risk signals.

BACKGROUND: Genome-wide association studies have identified 196 high confidence independent signals associated with breast cancer susceptibility. Variants within these signals frequently fall in distal regulatory DNA elements that control gene expression. RESULTS: We designed a Capture Hi-C array to enrich for chromatin interactions between the credible causal variants and target genes in six human mammary epithelial and breast cancer cell lines. We show that interacting regions are enriched for open chromatin, histone marks for active enhancers, and transcription factors relevant to breast biology. We exploit this comprehensive resource to identify candidate target genes at 139 independent breast cancer risk signals and explore the functional mechanism underlying altered risk at the 12q24 risk region. CONCLUSIONS: Our results demonstrate the power of combining genetics, computational genomics, and molecular studies to rationalize the identification of key variants and candidate target genes at breast cancer GWAS signals.

Non-coding RNAs underlie genetic predisposition to breast cancer.

BACKGROUND: Genetic variants identified through genome-wide association studies (GWAS) are predominantly non-coding and typically attributed to altered regulatory elements such as enhancers and promoters. However, the contribution of non-coding RNAs to complex traits is not clear. RESULTS: Using targeted RNA sequencing, we systematically annotated multi-exonic non-coding RNA (mencRNA) genes transcribed from 1.5-Mb intervals surrounding 139 breast cancer GWAS signals and assessed their contribution to breast cancer risk. We identify more than 4000 mencRNA genes and show their expression distinguishes normal breast tissue from tumors and different breast cancer subtypes. Importantly, breast cancer risk variants, identified through genetic fine-mapping, are significantly enriched in mencRNA exons, but not the promoters or introns. eQTL analyses identify mencRNAs whose expression is associated with risk variants. Furthermore, chromatin interaction data identify hundreds of mencRNA promoters that loop to regions that contain breast cancer risk variants. CONCLUSIONS: We have compiled the largest catalog of breast cancer-associated mencRNAs to date and provide evidence that modulation of mencRNAs by GWAS variants may provide an alternative mechanism underlying complex traits.

De Novo Variants in MAPK8IP3 Cause Intellectual Disability with Variable Brain Anomalies.

Using exome sequencing, we have identified de novo variants in MAPK8IP3 in 13 unrelated individuals presenting with an overlapping phenotype of mild to severe intellectual disability. The de novo variants comprise six missense variants, three of which are recurrent, and three truncating variants. Brain anomalies such as perisylvian polymicrogyria, cerebral or cerebellar atrophy, and hypoplasia of the corpus callosum were consistent among individuals harboring recurrent de novo missense variants. MAPK8IP3 has been shown to be involved in the retrograde axonal-transport machinery, but many of its specific functions are yet to be elucidated. Using the CRISPR-Cas9 system to target six conserved amino acid positions in Caenorhabditis elegans, we found that two of the six investigated human alterations led to a significantly elevated density of axonal lysosomes, and five variants were associated with adverse locomotion. Reverse-engineering normalized the observed adverse effects back to wild-type levels. Combining genetic, phenotypic, and functional findings, as well as the significant enrichment of de novo variants in MAPK8IP3 within our total cohort of 27,232 individuals who underwent exome sequencing, we implicate de novo variants in MAPK8IP3 as a cause of a neurodevelopmental disorder with intellectual disability and variable brain anomalies.

Sentryn and SAD Kinase Link the Guided Transport and Capture of Dense Core Vesicles in Caenorhabditis elegans.

Dense core vesicles (DCVs) can transmit signals by releasing neuropeptides from specialized synaptic regions called active zones. DCVs reach the active zone by motorized transport through a long axon. A reverse motor frequently interrupts progress by taking DCVs in the opposite direction. "Guided transport" refers to the mechanism by which outward movements ultimately dominate to bring DCVs to the synaptic region. After guided transport, DCVs alter their interactions with motors and enter a "captured" state. The mechanisms of guided transport and capture of DCVs are unknown. Here, we discovered two proteins that contribute to both processes in Caenorhabditis elegans SAD kinase and a novel conserved protein we named Sentryn are the first proteins found to promote DCV capture. By imaging DCVs moving in various regions of single identified neurons in living animals, we found that DCV guided transport and capture are linked through SAD kinase, Sentryn, and Liprin-α. These proteins act together to regulate DCV motorized transport in a region-specific manner. Between the cell body and the synaptic region, they promote forward transport. In the synaptic region, where all three proteins are highly enriched at active zones, they promote DCV pausing by inhibiting transport in both directions. These three proteins appear to be part of a special subset of active zone-enriched proteins because other active zone proteins do not share their unique functions.

Sentryn Acts with a Subset of Active Zone Proteins To Optimize the Localization of Synaptic Vesicles in Caenorhabditis elegans.

Synaptic vesicles (SVs) transmit signals by releasing neurotransmitters from specialized synaptic regions of neurons. In the synaptic region, SVs are tightly clustered around small structures called active zones. The motor KIF1A transports SVs outward through axons until they are captured in the synaptic region. This transport must be guided in the forward direction because it is opposed by the dynein motor, which causes SVs to reverse direction multiple times en route. The core synapse stability (CSS) system contributes to both guided transport and capture of SVs. We identified Sentryn as a CSS protein that contributes to the synaptic localization of SVs in Caenorhabditis elegans Like the CSS proteins SAD Kinase and SYD-2 (Liprin-α), Sentryn also prevents dynein-dependent accumulation of lysosomes in dendrites in strains lacking JIP3. Genetic analysis showed that Sentryn and SAD Kinase each have at least one nonoverlapping function for the stable accumulation of SVs at synapses that, when combined with their shared functions, enables most of the functions of SYD-2 (Liprin-α) for capturing SVs. Also like other CSS proteins, Sentryn appears enriched at active zones and contributes to active zone structure, suggesting that it is a novel, conserved active zone protein. Sentryn is recruited to active zones by a process dependent on the active zone-enriched CSS protein SYD-2 (Liprin-α). Our results define a specialized group of active zone enriched proteins that can affect motorized transport throughout the neuron and that have roles in both guided transport and capture of SVs.

MiR-29b-1-5p is altered in BRCA1 mutant tumours and is a biomarker in basal-like breast cancer.

Depletion of BRCA1 protein in mouse mammary glands results in defects in lactational development and increased susceptibility to mammary cancer. Extensive work has focussed on the role of BRCA1 in the normal breast and in the development of breast cancer, the cell of origin for BRCA1 tumours and the protein-coding genes altered in BRCA1 deficient cells. However, the role of non-coding RNAs in BRCA1-deficient cells is poorly understood. To evaluate miRNA expression in BRCA1 deficient mammary cells, RNA sequencing was performed on the mammary glands of Brca1 knockout mice. We identified 140 differentially expressed miRNAs, 9 of which were also differentially expressed in human BRCA1 breast tumours or familial non-BRCA1 patients and during normal gland development. We show that BRCA1 binds to putative cis-elements in promoter regions of the miRNAs with the potential to regulate their expression, and that four miRNAs (miR-29b-1-5p, miR-664, miR-16-2 and miR-744) significantly stratified the overall survival of basal-like tumours. Importantly the prognostic value of miR-29b-1-5p was higher in significance than several commonly used clinical biomarkers. These results emphasise the role of Brca1 in modulating expression of miRNAs and highlights the potential for BRCA1 regulated miRNAs to be informative biomarkers associated with BRCA1 loss and survival in breast cancer.

Local microtubule organization promotes cargo transport in C. elegans dendrites.

The specific organization of the neuronal microtubule cytoskeleton in axons and dendrites is an evolutionarily conserved determinant of neuronal polarity that allows for selective cargo sorting. However, how dendritic microtubules are organized and whether local differences influence cargo transport remains largely unknown. Here, we use live-cell imaging to systematically probe the microtubule organization in Caenorhabditiselegans neurons, and demonstrate the contribution of distinct mechanisms in the organization of dendritic microtubules. We found that most non-ciliated neurons depend on unc-116 (kinesin-1), unc-33 (CRMP) and unc-44 (ankyrin) for correct microtubule organization and polarized cargo transport, as previously reported. Ciliated neurons and the URX neuron, however, use an additional pathway to nucleate microtubules at the tip of the dendrite, from the base of the cilium in ciliated neurons. Since inhibition of distal microtubule nucleation affects distal dendritic transport, we propose a model in which the presence of a microtubule-organizing center at the dendrite tip ensures correct dendritic cargo transport.

BRCA1 and BRCA2 5' noncoding region variants identified in breast cancer patients alter promoter activity and protein binding.

The widespread use of next generation sequencing for clinical testing is detecting an escalating number of variants in noncoding regions of the genome. The clinical significance of the majority of these variants is currently unknown, which presents a significant clinical challenge. We have screened over 6,000 early-onset and/or familial breast cancer (BC) cases collected by the ENIGMA consortium for sequence variants in the 5' noncoding regions of BC susceptibility genes BRCA1 and BRCA2, and identified 141 rare variants with global minor allele frequency < 0.01, 76 of which have not been reported previously. Bioinformatic analysis identified a set of 21 variants most likely to impact transcriptional regulation, and luciferase reporter assays detected altered promoter activity for four of these variants. Electrophoretic mobility shift assays demonstrated that three of these altered the binding of proteins to the respective BRCA1 or BRCA2 promoter regions, including NFYA binding to BRCA1:c.-287C>T and PAX5 binding to BRCA2:c.-296C>T. Clinical classification of variants affecting promoter activity, using existing prediction models, found no evidence to suggest that these variants confer a high risk of disease. Further studies are required to determine if such variation may be associated with a moderate or low risk of BC.

Publisher Correction: Evidence that breast cancer risk at the 2q35 locus is mediated through IGFBP5 regulation.

This corrects the article DOI: 10.1038/ncomms5999.

Functional dissection of breast cancer risk-associated TERT promoter variants.

The multi-cancer susceptibility locus at 5p15.33 includes TERT, encoding the telomerase catalytic subunit. Genome-wide association studies (GWAS) have identified six single nucleotide polymorphisms (SNPs) in the TERT promoter associated with decreased breast cancer risk, although the precise causal variants and their mechanisms of action have remained elusive. Luciferase reporter assays indicated that the protective haplotype reduced TERT promoter activity in human mammary epithelial and cancer cells in an estrogen-independent manner. Using single variant constructs, we identified rs3215401 and rs2853669 as likely functional variants. Silencing of MYC decreased TERT promoter activity but neither MYC nor ETS2 silencing conferred allele-specificity. In chromatin immunoprecipitation experiments, the ETS protein GABPA, but not ETS2 or ELF1, bound rs2853669 in an allele-specific manner in mammary epithelial cells. Investigation of open chromatin in mammoplasty samples suggested involvement of three additional variants, though not rs3215401 or rs2853669. Chromosome conformation capture revealed no interaction of the TERT promoter with regulatory elements in the locus, indicating limited local impact of candidate variants on the TERT promoter. Collectively, our functional studies of the TERT-CLPTM1L breast cancer susceptibility locus describe rs2853669 as a functional variant of this association signal among three other potentially causal variants and demonstrate the versatile mechanisms by which TERT promoter variants may affect breast cancer risk.

Long Noncoding RNAs CUPID1 and CUPID2 Mediate Breast Cancer Risk at 11q13 by Modulating the Response to DNA Damage.

Breast cancer risk is strongly associated with an intergenic region on 11q13. We have previously shown that the strongest risk-associated SNPs fall within a distal enhancer that regulates CCND1. Here, we report that, in addition to regulating CCND1, this enhancer regulates two estrogen-regulated long noncoding RNAs, CUPID1 and CUPID2. We provide evidence that the risk-associated SNPs are associated with reduced chromatin looping between the enhancer and the CUPID1 and CUPID2 bidirectional promoter. We further show that CUPID1 and CUPID2 are predominantly expressed in hormone-receptor-positive breast tumors and play a role in modulating pathway choice for the repair of double-strand breaks. These data reveal a mechanism for the involvement of this region in breast cancer.