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.

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.

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.

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.

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.

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.

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.

A Common Variant at the 14q32 Endometrial Cancer Risk Locus Activates AKT1 through YY1 Binding.

A recent meta-analysis of multiple genome-wide association and follow-up endometrial cancer case-control datasets identified a novel genetic risk locus for this disease at chromosome 14q32.33. To prioritize the functional SNP(s) and target gene(s) at this locus, we employed an in silico fine-mapping approach using genotyped and imputed SNP data for 6,608 endometrial cancer cases and 37,925 controls of European ancestry. Association and functional analyses provide evidence that the best candidate causal SNP is rs2494737. Multiple experimental analyses show that SNP rs2494737 maps to a silencer element located within AKT1, a member of the PI3K/AKT/MTOR intracellular signaling pathway activated in endometrial tumors. The rs2494737 risk A allele creates a YY1 transcription factor-binding site and abrogates the silencer activity in luciferase assays, an effect mimicked by transfection of YY1 siRNA. Our findings suggest YY1 is a positive regulator of AKT1, mediating the stimulatory effects of rs2494737 increasing endometrial cancer risk. Identification of an endometrial cancer risk allele within a member of the PI3K/AKT signaling pathway, more commonly activated in tumors by somatic alterations, raises the possibility that well tolerated inhibitors targeting this pathway could be candidates for evaluation as chemopreventive agents in individuals at high risk of developing endometrial cancer.

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.

Long-range regulators of the lncRNA HOTAIR enhance its prognostic potential in breast cancer.

Predicting response to endocrine therapy and survival in oestrogen receptor positive breast cancer is a significant clinical challenge and novel prognostic biomarkers are needed. Long-range regulators of gene expression are emerging as promising biomarkers and therapeutic targets for human diseases, so we have explored the potential of distal enhancer elements of non-coding RNAs in the prognostication of breast cancer survival. HOTAIR is a long non-coding RNA that is overexpressed, promotes metastasis and is predictive of decreased survival. Here, we describe a long-range transcriptional enhancer of the HOTAIR gene that binds several hormone receptors and associated transcription factors, interacts with the HOTAIR promoter and augments transcription. This enhancer is dependent on Forkhead-Box transcription factors and functionally interacts with a novel alternate HOTAIR promoter. HOTAIR expression is negatively regulated by oestrogen, positively regulated by FOXA1 and FOXM1, and is inversely correlated with oestrogen receptor and directly correlated with FOXM1 in breast tumours. The combination of HOTAIR and FOXM1 enables greater discrimination of endocrine therapy responders and non-responders in patients with oestrogen receptor positive breast cancer. Consistent with this, HOTAIR expression is increased in cell-line models of endocrine resistance. Analysis of breast cancer gene expression data indicates that HOTAIR is co-expressed with FOXA1 and FOXM1 in HER2-enriched tumours, and these factors enhance the prognostic power of HOTAIR in aggressive HER2+ breast tumours. Our study elucidates the transcriptional regulation of HOTAIR, identifies HOTAIR and its regulators as novel biomarkers of patient response to endocrine therapy and corroborates the importance of transcriptional enhancers in cancer.

Endometriosis risk alleles at 1p36.12 act through inverse regulation of CDC42 and LINC00339.

Genome-wide association studies (GWAS) have identified markers within the WNT4 region on chromosome 1p36.12 showing consistent and strong association with increasing endometriosis risk. Fine mapping using sequence and imputed genotype data has revealed strong candidates for the causal SNPs within these critical regions; however, the molecular pathogenesis of these SNPs is currently unknown. We used gene expression data collected from whole blood from 862 individuals and endometrial tissue from 136 individuals from independent populations of European descent to examine the mechanism underlying endometriosis susceptibility. Association mapping results from 7,090 individuals (2,594 cases and 4,496 controls) supported rs3820282 as the SNP with the strongest association for endometriosis risk (P = 1.84 × 10−5, OR = 1.244 (1.126-1.375)). SNP rs3820282 is a significant eQTL in whole blood decreasing expression of LINC00339 (also known as HSPC157) and increasing expression of CDC42 (P = 2.0 ×10−54 and 4.5x10−4 respectively). The largest effects were for two LINC00339 probes (P = 2.0 ×10−54; 1.0 × 10−34). The eQTL for LINC00339 was also observed in endometrial tissue (P = 2.4 ×10−8) with the same direction of effect for both whole blood and endometrial tissue. There was no evidence for eQTL effects for WNT4. Chromatin conformation capture provides evidence for risk SNPs interacting with the promoters of both LINC00339 and CDC4 and luciferase reporter assays suggest the risk SNP rs12038474 is located in a transcriptional silencer for CDC42 and the risk allele increases expression of CDC42. However, no effect of rs3820282 was observed in the LINC00339 expression in Ishikawa cells. Taken together, our results suggest that SNPs increasing endometriosis risk in this region act through CDC42, but further functional studies are required to rule out inverse regulation of both LINC00339 and CDC42.

Long-Range Modulation of PAG1 Expression by 8q21 Allergy Risk Variants.

The gene(s) whose expression is regulated by allergy risk variants is unknown for many loci identified through genome-wide association studies. Addressing this knowledge gap might point to new therapeutic targets for allergic disease. The aim of this study was to identify the target gene(s) and the functional variant(s) underlying the association between rs7009110 on chromosome 8q21 and allergies. Eight genes are located within 1 Mb of rs7009110. Multivariate association analysis of publicly available exon expression levels from lymphoblastoid cell lines (LCLs) identified a significant association between rs7009110 and the expression of a single gene, PAG1 (p = 0.0017), 732 kb away. Analysis of histone modifications and DNase I hypersensitive sites in LCLs identified four putative regulatory elements (PREs) in the region. Chromosome conformation capture confirmed that two PREs interacted with the PAG1 promoter, one in allele-specific fashion. To determine whether these PREs were functional, LCLs were transfected with PAG1 promoter-driven luciferase reporter constructs. PRE3 acted as a transcriptional enhancer for PAG1 exclusively when it carried the rs2370615:C allergy predisposing allele, a variant in complete linkage disequilibrium with rs7009110. As such, rs2370615, which overlaps RelA transcription factor (TF) binding in LCLs and was found to disrupt Foxo3a binding to PRE3, represents the putative functional variant in this locus. Our studies suggest that the risk-associated allele of rs2370615 predisposes to allergic disease by increasing PAG1 expression, which might promote B cell activation and have a pro-inflammatory effect. Inhibition of PAG1 expression or function might have therapeutic potential for allergic diseases.

Beyond GWASs: illuminating the dark road from association to function.

Genome-wide association studies (GWASs) have enabled the discovery of common genetic variation contributing to normal and pathological traits and clinical drug responses, but recognizing the precise targets of these associations is now the major challenge. Here, we review recent approaches to the functional follow-up of GWAS loci, including fine mapping of GWAS signal(s), prioritization of putative functional SNPs by the integration of genetic epidemiological and bioinformatic methods, and in vitro and in vivo experimental verification of predicted molecular mechanisms for identifying the targeted genes. The majority of GWAS-identified variants fall in noncoding regions of the genome. Therefore, this review focuses on strategies for assessing likely mechanisms affected by noncoding variants; such mechanisms include transcriptional regulation, noncoding RNA function, and epigenetic regulation. These approaches have already accelerated progress from genetic studies to biological knowledge and might ultimately guide the development of prognostic, preventive, and therapeutic measures.

Functional variants at the 11q13 risk locus for breast cancer regulate cyclin D1 expression through long-range enhancers.

Analysis of 4,405 variants in 89,050 European subjects from 41 case-control studies identified three independent association signals for estrogen-receptor-positive tumors at 11q13. The strongest signal maps to a transcriptional enhancer element in which the G allele of the best candidate causative variant rs554219 increases risk of breast cancer, reduces both binding of ELK4 transcription factor and luciferase activity in reporter assays, and may be associated with low cyclin D1 protein levels in tumors. Another candidate variant, rs78540526, lies in the same enhancer element. Risk association signal 2, rs75915166, creates a GATA3 binding site within a silencer element. Chromatin conformation studies demonstrate that these enhancer and silencer elements interact with each other and with their likely target gene, CCND1.

Fine-scale mapping of the FGFR2 breast cancer risk locus: putative functional variants differentially bind FOXA1 and E2F1.

The 10q26 locus in the second intron of FGFR2 is the locus most strongly associated with estrogen-receptor-positive breast cancer in genome-wide association studies. We conducted fine-scale mapping in case-control studies genotyped with a custom chip (iCOGS), comprising 41 studies (n = 89,050) of European ancestry, 9 Asian ancestry studies (n = 13,983), and 2 African ancestry studies (n = 2,028) from the Breast Cancer Association Consortium. We identified three statistically independent risk signals within the locus. Within risk signals 1 and 3, genetic analysis identified five and two variants, respectively, highly correlated with the most strongly associated SNPs. By using a combination of genetic fine mapping, data on DNase hypersensitivity, and electrophoretic mobility shift assays to study protein-DNA binding, we identified rs35054928, rs2981578, and rs45631563 as putative functional SNPs. Chromatin immunoprecipitation showed that FOXA1 preferentially bound to the risk-associated allele (C) of rs2981578 and was able to recruit ERα to this site in an allele-specific manner, whereas E2F1 preferentially bound the risk variant of rs35054928. The risk alleles were preferentially found in open chromatin and bound by Ser5 phosphorylated RNA polymerase II, suggesting that the risk alleles are associated with changes in transcription. Chromatin conformation capture demonstrated that the risk region was able to interact with the promoter of FGFR2, the likely target gene of this risk region. A role for FOXA1 in mediating breast cancer susceptibility at this locus is consistent with the finding that the FGFR2 risk locus primarily predisposes to estrogen-receptor-positive disease.

Resistance to therapy caused by intragenic deletion in BRCA2.

Cells with loss of BRCA2 function are defective in homologous recombination (HR) and are highly sensitive to inhibitors of poly(ADP-ribose) polymerase (PARP), which provides the basis for a new therapeutic approach. Here we show that resistance to PARP inhibition can be acquired by deletion of a mutation in BRCA2. We derived PARP-inhibitor-resistant (PIR) clones from the human CAPAN1 pancreatic cancer cell line, which carries the protein-truncating c.6174delT frameshift mutation. PIR clones could form DNA-damage-induced RAD51 nuclear foci and were able to limit genotoxin-induced genomic instability, both hallmarks of a competent HR pathway. New BRCA2 isoforms were expressed in the resistant lines as a result of intragenic deletion of the c.6174delT mutation and restoration of the open reading frame (ORF). Reconstitution of BRCA2-deficient cells with these revertant BRCA2 alleles rescued PARP inhibitor sensitivity and HR deficiency. Most of the deletions in BRCA2 were associated with small tracts of homology, and possibly arose from error-prone repair caused by BRCA2 deficiency. Similar ORF-restoring mutations were present in carboplatin-resistant ovarian tumours from c.6174delT mutation carriers. These observations have implications for understanding drug resistance in BRCA mutation carriers as well as in defining functionally important domains within BRCA2.

Long-range transcriptional regulation of breast cancer genes.

Breast cancer is a major health problem and understanding the genetic basis of this disease is crucial for predicting risk and developing effective targeted therapeutics. Several breast cancer predisposing genes have been identified, but mutations in the coding regions of these genes only accounts for a small proportion of risk. Research now suggests that combinations of multiple non-coding changes in breast cancer susceptibility genes, which cause moderate alterations in gene expression, will be responsible for the remaining inherited risk. These non-coding changes will include variants in proximal and distal transcriptional and post-transcriptional regulatory elements and may affect the levels and function of trans-acting factors, including proteins and RNAs, which act on these elements. Somatic changes in such elements and factors have also been associated with breast cancer progression. With the recent advent of techniques allowing the detection of long-range DNA interactions spanning the human genome, it has become increasingly clear that long-range regulatory elements constitute an important mechanism for gene regulation. Recent studies have identified several such elements that are important for regulating genes involved in breast cancer, raising the possibility that defects in these sequences may contribute to breast cancer predisposition and progression. In this review, we discuss the emerging functions of cis-regulatory elements and a subset of trans-acting factors in breast tumorigenesis. We also discuss some recent progress in our understanding of how dysregulation in these transcriptional components may contribute to breast cancer, and the potential implications for molecular diagnosis, prognosis prediction, and the treatment of this disease.

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.

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.

Identification of fifteen novel germline variants in the BRCA1 3'UTR reveals a variant in a breast cancer case that introduces a functional miR-103 target site.

Mutations in the BRCA1 gene confer a substantial increase in breast cancer risk, yet routine clinical genetic screening is limited to the coding regions and intron-exon boundaries, precluding the identification of mutations in noncoding and untranslated regions (UTR). As 3'UTR mutations can influence cancer susceptibility by altering protein and microRNA (miRNA) binding regions, we screened the BRCA1 3'UTR for mutations in a large series of BRCA-mutation negative, population and clinic-based breast cancer cases, and controls. Fifteen novel BRCA1 3'UTR variants were identified, the majority of which were unique to either cases or controls. Using luciferase reporter assays, three variants found in cases, c.* 528G>C, c.* 718A>G, and c.* 1271T>C and four found in controls, c.* 309T>C, c.* 379G>A, c.* 823C>T, and c.* 264C>T, reduced 3'UTR activity (P < 0.02), whereas two variants found in cases, c.* 291C>T and c.* 1139G>T, increased 3'UTR activity (P < 0.01). Three case variants, c.* 718A>G, c.* 800T>C, and c.* 1340_1342delTGT, were predicted to create new miRNA binding sites and c.* 1340_1342delTGT caused a reduction (25%, P = 0.0007) in 3'UTR reporter activity when coexpressed with the predicted targeting miRNA, miR-103. This is the most comprehensive identification and analysis of BRCA1 3'UTR variants published to date.