Cis-eQTL-based trans-ethnic meta-analysis reveals novel genes associated with breast cancer risk

Multi-tissue transcriptome-wide association studies identified 235 genes for intrinsic subtypes of breast cancer

BACKGROUND

Although genome-wide association studies (GWAS) of breast cancer (BC) identified common variants which differ between intrinsic subtypes, genes through which these variants act to impact BC risk have not been fully established. Transcriptome-wide association studies (TWAS) have identified genes associated with overall BC risk, but subtype-specific differences are largely unknown.

METHODS

We performed two multi-tissue TWAS for each BC intrinsic subtype, including an expression-based approach that collated TWAS signals from expression quantitative trait loci (eQTLs) across multiple tissues and a novel splicing-based approach that collated signals from splicing QTLs (sQTLs) across intron clusters and subsequently across tissues. We used summary statistics for five intrinsic subtypes including Luminal A-like, Luminal B-like, Luminal B/HER2-negative-like, HER2-enriched-like, and triple-negative BC, generated from 106 278 BC cases and 91 477 controls in the Breast Cancer Association Consortium.

RESULTS

Overall, we identified 235 genes in 88 loci that were associated with at least one of the five intrinsic subtypes. Most genes were subtype-specific, and many have not been reported in previous TWAS. We discovered common variants that modulate expression of CHEK2 confer increased risk to Luminal A-like BC, in contrast to the viewpoint that CHEK2 primarily harbors rare, penetrant mutations. Additionally, our splicing-based TWAS provided population-level support for MDM4 splice variants that increased the risk of triple-negative BC.

CONCLUSION

Our comprehensive, multi-tissue TWAS corroborated previous GWAS loci for overall BC risk and intrinsic subtypes, while underscoring how common variation that impacts expression and splicing of genes in multiple tissue types can be used to further elucidate the etiology of BC.

A multi-tissue, splicing-based joint transcriptome-wide association study identifies susceptibility genes for breast cancer

Splicing-based transcriptome-wide association studies (splicing-TWASs) of breast cancer have the potential to identify susceptibility genes. However, existing splicing-TWASs test the association of individual excised introns in breast tissue only and thus have limited power to detect susceptibility genes. In this study, we performed a multi-tissue joint splicing-TWAS that integrated splicing-TWAS signals of multiple excised introns in each gene across 11 tissues that are potentially relevant to breast cancer risk. We utilized summary statistics from a meta-analysis that combined genome-wide association study (GWAS) results of 424,650 women of European ancestry. Splicing-level prediction models were trained in GTEx (v.8) data. We identified 240 genes by the multi-tissue joint splicing-TWAS at the Bonferroni-corrected significance level; in the tissue-specific splicing-TWAS that combined TWAS signals of excised introns in genes in breast tissue only, we identified nine additional significant genes. Of these 249 genes, 88 genes in 62 loci have not been reported by previous TWASs, and 17 genes in seven loci are at least 1 Mb away from published GWAS index variants. By comparing the results of our splicing-TWASs with previous gene-expression-based TWASs that used the same summary statistics and expression prediction models trained in the same reference panel, we found that 110 genes in 70 loci that are identified only by the splicing-TWASs. Our results showed that for many genes, expression quantitative trait loci (eQTL) did not show a significant impact on breast cancer risk, whereas splicing quantitative trait loci (sQTL) showed a strong impact through intron excision events.

Cis- and trans-eQTL TWASs of breast and ovarian cancer identify more than 100 susceptibility genes in the BCAC and OCAC consortia

Transcriptome-wide association studies (TWASs) have investigated the role of genetically regulated transcriptional activity in the etiologies of breast and ovarian cancer. However, methods performed to date have focused on the regulatory effects of risk-associated SNPs thought to act in cis on a nearby target gene. With growing evidence for distal (trans) regulatory effects of variants on gene expression, we performed TWASs of breast and ovarian cancer using a Bayesian genome-wide TWAS method (BGW-TWAS) that considers effects of both cis- and trans-expression quantitative trait loci (eQTLs). We applied BGW-TWAS to whole-genome and RNA sequencing data in breast and ovarian tissues from the Genotype-Tissue Expression project to train expression imputation models. We applied these models to large-scale GWAS summary statistic data from the Breast Cancer and Ovarian Cancer Association Consortia to identify genes associated with risk of overall breast cancer, non-mucinous epithelial ovarian cancer, and 10 cancer subtypes. We identified 101 genes significantly associated with risk with breast cancer phenotypes and 8 with ovarian phenotypes. These loci include established risk genes and several novel candidate risk loci, such as ACAP3, whose associations are predominantly driven by trans-eQTLs. We replicated several associations using summary statistics from an independent GWAS of these cancer phenotypes. We further used genotype and expression data in normal and tumor breast tissue from the Cancer Genome Atlas to examine the performance of our trained expression imputation models. This work represents an in-depth look into the role of trans eQTLs in the complex molecular mechanisms underlying these diseases.

Using genome and transcriptome data from African-ancestry female participants to identify putative breast cancer susceptibility genes

African-ancestry (AA) participants are underrepresented in genetics research. Here, we conducted a transcriptome-wide association study (TWAS) in AA female participants to identify putative breast cancer susceptibility genes. We built genetic models to predict levels of gene expression, exon junction, and 3' UTR alternative polyadenylation using genomic and transcriptomic data generated in normal breast tissues from 150 AA participants and then used these models to perform association analyses using genomic data from 18,034 cases and 22,104 controls. At Bonferroni-corrected P < 0.05, we identified six genes associated with breast cancer risk, including four genes not previously reported (CTD-3080P12.3, EN1, LINC01956 and NUP210L). Most of these genes showed a stronger association with risk of estrogen-receptor (ER) negative or triple-negative than ER-positive breast cancer. We also replicated the associations with 29 genes reported in previous TWAS at P < 0.05 (one-sided), providing further support for an association of these genes with breast cancer risk. Our study sheds new light on the genetic basis of breast cancer and highlights the value of conducting research in AA populations.

Expression- and splicing-based multi-tissue transcriptome-wide association studies identified multiple genes for breast cancer by estrogen-receptor status

BACKGROUND

Although several transcriptome-wide association studies (TWASs) have been performed to identify genes associated with overall breast cancer (BC) risk, only a few TWAS have explored the differences in estrogen receptor-positive (ER+) and estrogen receptor-negative (ER-) breast cancer. Additionally, these studies were based on gene expression prediction models trained primarily in breast tissue, and they did not account for alternative splicing of genes.

METHODS

In this study, we utilized two approaches to perform multi-tissue TWASs of breast cancer by ER subtype: (1) an expression-based TWAS that combined TWAS signals for each gene across multiple tissues and (2) a splicing-based TWAS that combined TWAS signals of all excised introns for each gene across tissues. To perform this TWAS, we utilized summary statistics for ER + BC from the Breast Cancer Association Consortium (BCAC) and for ER- BC from a meta-analysis of BCAC and the Consortium of Investigators of Modifiers of BRCA1 and BRCA2 (CIMBA).

RESULTS

In total, we identified 230 genes in 86 loci that were associated with ER + BC and 66 genes in 29 loci that were associated with ER- BC at a Bonferroni threshold of significance. Of these genes, 2 genes associated with ER + BC at the 1q21.1 locus were located at least 1 Mb from published GWAS hits. For several well-studied tumor suppressor genes such as TP53 and CHEK2 which have historically been thought to impact BC risk through rare, penetrant mutations, we discovered that common variants, which modulate gene expression, may additionally contribute to ER + or ER- etiology.

CONCLUSIONS

Our study comprehensively examined how differences in common variation contribute to molecular differences between ER + and ER- BC and introduces a novel, splicing-based framework that can be used in future TWAS studies.

Multi-omics Mendelian randomization integrating GWAS, eQTL, and mQTL data identified genes associated with breast cancer

Breast cancer (BC) remains a major disease posing a threat to women's health, but the underlying biological interpretation remains largely unknown. Here, we aimed to identify genes associated with breast cancer and analyze their pathophysiological mechanisms based on multi-omics Mendelian randomization (MR). Summary-data-based MR (SMR) was performed to estimate the causal effects of blood and breast mammary tissue expression quantitative trait loci (eQTLs) on BC. External validation analysis was used to validate the identified genes. Integration analyses BC GWAS summaries with eQTLs and DNA methylation QTLs (mQTLs) from the blood were conducted using SMR to prioritize putative blood genes and their regulatory elements associated with BC risk. Finally, two prior genes (ATG10 and RCCD1) from blood tissue reached significant levels in both BCAC (ATG10: ORBRCR = 0.91, PBRCR = 1.29 × 10-11; RCCD1: ORBRCR = 0.90, PBRCR = 3.72 × 10-15) and FinnGen cohorts (ATG10: ORFinnGen = 0.89, PFinnGen = 8.55 × 10-5; RCCD1: ORFinnGen = 0.89, PFinnGen = 2.38 × 10-8). Additionally, those two genes from breast tissues also replicated in both BCAC (ATG10: ORBRCR = 0.95, PBRCR = 1.02 × 10-9; RCCD1: ORBRCR = 0.87, PBRCR = 4.70 × 10-10) and FinnGen cohorts (ATG10: ORFinnGen = 0.93, PFinnGen = 2.38 × 10-4; RCCD1: ORFinnGen = 0.85, PFinnGen = 3.81 × 10-6). Sensitive analysis and external validation analysis validated those two identified genes. Multi-omics MR analysis showed that the SNP signals associated with ATG10 and RCCD1 were significant across the data from BC Genome-wide association study (GWAS), eQTL, and mQTL studies. In conclusion, we identified two priority genes that are potentially associated with BC. These findings improve our limited understanding of the mechanism of BC and shed light on the development of therapeutic agents for treating BC.

LRRC25 expression during physiological aging and in mouse models of Alzheimer's disease and iPSC-derived neurons

The leucine-rich repeat-containing protein 25 (LRRC25) is relatively a novel protein with no information on its role in neuronal or brain function. A recent study suggested LRRC25 is a potential risk factor for Alzheimer's disease (AD). As a first step to understanding LRRC25's role in the brain and AD, we found LRRC25 is expressed in both cell membranes and cytoplasm in a punctuate appearance in astrocytes, microglia, and neurons in cell lines as well as mouse brain. We also found that LRRC25 expression is both age- and brain region-dependent and that 1-day-old (1D) pups expressed the least amount of LRRC25 protein compared to adult ages. In the APΔE9 mice, immunoblot quantified LRRC25 protein levels were increased by 166% (**p < 0.01) in the cortex (CX) and by 215% (***p < 0.001) in the hippocampus (HP) relative to wild-type (WT) controls. Both the brainstem (BS) and cerebellum (CB) showed no significant alterations. In the 3xTg mice, only CX showed an increase of LRRC25 protein by 91% (*p < 0.05) when compared to WT controls although the increased trend was noted in the other brain regions. In the AD patient brains also LRRC25 protein levels were increased by 153% (***p < 0.001) when compared to normal control (NC) subjects. Finally, LRRC25 expression in the iPSC-derived neurons quantified by immunofluorescence was increased by 181% (**p < 0.01) in AD-derived neurons when compared to NC-derived neurons. Thus increased LRRC25 protein in multiple models of AD suggests that LRRC25 may play a pathogenic role in either Aβ or tau pathology in AD. The mechanism for the increased levels of LRRC25 in AD is unknown at present, but a previous study showed that LRRC25 levels also increase during neonatal hypoxic-ischemia neuronal damage. Based on the evidence that autophagy is highly dysregulated in AD, the increased LRRC25 levels may be due to decreased autophagic degradation of LRRC25. Increased LRRC25 in turn may regulate the stability or activity of key enzymes involved in either Aβ or hyperphosphorylated tau generation and thus may contribute to increased plaques and neurofibrillary tangles.

Transcriptome-Wide Association Studies (TWAS): Methodologies, Applications, and Challenges

Transcriptome-wide association study (TWAS) methodologies aim to identify genetic effects on phenotypes through the mediation of gene transcription. In TWAS, in silico models of gene expression are trained as functions of genetic variants and then applied to genome-wide association study (GWAS) data. This post-GWAS analysis identifies gene-trait associations with high interpretability, enabling follow-up functional genomics studies and the development of genetics-anchored resources. We provide an overview of commonly used TWAS approaches, their advantages and limitations, and some widely used applications. © 2024 Wiley Periodicals LLC.

Cis- and trans-eQTL TWAS of breast and ovarian cancer identify more than 100 risk associated genes in the BCAC and OCAC consortia

Transcriptome-wide association studies (TWAS) have investigated the role of genetically regulated transcriptional activity in the etiologies of breast and ovarian cancer. However, methods performed to date have only considered regulatory effects of risk associated SNPs thought to act in cis on a nearby target gene. With growing evidence for distal (trans) regulatory effects of variants on gene expression, we performed TWAS of breast and ovarian cancer using a Bayesian genome-wide TWAS method (BGW-TWAS) that considers effects of both cis- and trans-expression quantitative trait loci (eQTLs). We applied BGW-TWAS to whole genome and RNA sequencing data in breast and ovarian tissues from the Genotype-Tissue Expression project to train expression imputation models. We applied these models to large-scale GWAS summary statistic data from the Breast Cancer and Ovarian Cancer Association Consortia to identify genes associated with risk of overall breast cancer, non-mucinous epithelial ovarian cancer, and 10 cancer subtypes. We identified 101 genes significantly associated with risk with breast cancer phenotypes and 8 with ovarian phenotypes. These loci include established risk genes and several novel candidate risk loci, such as ACAP3, whose associations are predominantly driven by trans-eQTLs. We replicated several associations using summary statistics from an independent GWAS of these cancer phenotypes. We further used genotype and expression data in normal and tumor breast tissue from the Cancer Genome Atlas to examine the performance of our trained expression imputation models. This work represents a first look into the role of trans-eQTLs in the complex molecular mechanisms underlying these diseases.

Polymorphic variants involved in methylation regulation: a strategy to discover risk loci for pancreatic ductal adenocarcinoma

INTRODUCTION

Only a small number of risk factors for pancreatic ductal adenocarcinoma (PDAC) has been established. Several studies identified a role of epigenetics and of deregulation of DNA methylation. DNA methylation is variable across a lifetime and in different tissues; nevertheless, its levels can be regulated by genetic variants like methylation quantitative trait loci (mQTLs), which can be used as a surrogate.

MATERIALS AND METHODS

We scanned the whole genome for mQTLs and performed an association study in 14 705 PDAC cases and 246 921 controls. The methylation data were obtained from whole blood and pancreatic cancer tissue through online databases. We used the Pancreatic Cancer Cohort Consortium and the Pancreatic Cancer Case-Control Consortium genome-wide association study (GWAS) data as discovery phase and the Pancreatic Disease Research consortium, the FinnGen project and the Japan Pancreatic Cancer Research consortium GWAS as replication phase.

RESULTS

The C allele of 15q26.1-rs12905855 showed an association with a decreased risk of PDAC (OR=0.90, 95% CI 0.87 to 0.94, p=4.93×10-8 in the overall meta-analysis), reaching genome-level statistical significance. 15q26.1-rs12905855 decreases the methylation of a 'C-phosphate-G' (CpG) site located in the promoter region of the RCCD1 antisense (RCCD1-AS1) gene which, when expressed, decreases the expression of the RCC1 domain-containing (RCCD1) gene (part of a histone demethylase complex). Thus, it is possible that the rs12905855 C-allele has a protective role in PDAC development through an increase of RCCD1 gene expression, made possible by the inactivity of RCCD1-AS1.

CONCLUSION

We identified a novel PDAC risk locus which modulates cancer risk by controlling gene expression through DNA methylation.

A joint transcriptome-wide association study across multiple tissues identifies candidate breast cancer susceptibility genes

Genome-wide association studies (GWASs) have identified more than 200 genomic loci for breast cancer risk, but specific causal genes in most of these loci have not been identified. In fact, transcriptome-wide association studies (TWASs) of breast cancer performed using gene expression prediction models trained in breast tissue have yet to clearly identify most target genes. To identify candidate genes, we performed a GWAS analysis in a breast cancer dataset from UK Biobank (UKB) and combined the results with the GWAS results of the Breast Cancer Association Consortium (BCAC) by a meta-analysis. Using the summary statistics from the meta-analysis, we performed a joint TWAS analysis that combined TWAS signals from multiple tissues. We used expression prediction models trained in 11 tissues that are potentially relevant to breast cancer from the Genotype-Tissue Expression (GTEx) data. In the GWAS analysis, we identified eight loci distinct from those reported previously. In the TWAS analysis, we identified 309 genes at 108 genomic loci to be significantly associated with breast cancer at the Bonferroni threshold. Of these, 17 genes were located in eight regions that were at least 1 Mb away from published GWAS hits. The remaining TWAS-significant genes were located in 100 known genomic loci from previous GWASs of breast cancer. We found that 21 genes located in known GWAS loci remained statistically significant after conditioning on previous GWAS index variants. Our study provides insights into breast cancer genetics through mapping candidate target genes in a large proportion of known GWAS loci and discovering multiple new loci.

OTTERS: a powerful TWAS framework leveraging summary-level reference data

Most existing TWAS tools require individual-level eQTL reference data and thus are not applicable to summary-level reference eQTL datasets. The development of TWAS methods that can harness summary-level reference data is valuable to enable TWAS in broader settings and enhance power due to increased reference sample size. Thus, we develop a TWAS framework called OTTERS (Omnibus Transcriptome Test using Expression Reference Summary data) that adapts multiple polygenic risk score (PRS) methods to estimate eQTL weights from summary-level eQTL reference data and conducts an omnibus TWAS. We show that OTTERS is a practical and powerful TWAS tool by both simulations and application studies.

ANKLE1 cleaves mitochondrial DNA and contributes to cancer risk by promoting apoptosis resistance and metabolic dysregulation

Alleles within the chr19p13.1 locus are associated with increased risk of both ovarian and breast cancer and increased expression of the ANKLE1 gene. ANKLE1 is molecularly characterized as an endonuclease that efficiently cuts branched DNA and shuttles between the nucleus and cytoplasm. However, the role of ANKLE1 in mammalian development and homeostasis remains unknown. In normal development ANKLE1 expression is limited to the erythroblast lineage and we found that ANKLE1's role is to cleave the mitochondrial genome during erythropoiesis. We show that ectopic expression of ANKLE1 in breast epithelial-derived cells leads to genome instability and mitochondrial DNA (mtDNA) cleavage. mtDNA degradation then leads to mitophagy and causes a shift from oxidative phosphorylation to glycolysis (Warburg effect). Moreover, mtDNA degradation activates STAT1 and expression of epithelial-mesenchymal transition (EMT) genes. Reduction in mitochondrial content contributes to apoptosis resistance, which may allow precancerous cells to avoid apoptotic checkpoints and proliferate. These findings provide evidence that ANKLE1 is the causal cancer susceptibility gene in the chr19p13.1 locus and describe mechanisms by which higher ANKLE1 expression promotes cancer risk.

MiXcan: a framework for cell-type-aware transcriptome-wide association studies with an application to breast cancer

Human bulk tissue samples comprise multiple cell types with diverse roles in disease etiology. Conventional transcriptome-wide association study approaches predict genetically regulated gene expression at the tissue level, without considering cell-type heterogeneity, and test associations of predicted tissue-level expression with disease. Here we develop MiXcan, a cell-type-aware transcriptome-wide association study approach that predicts cell-type-level expression, identifies disease-associated genes via combination of cell-type-level association signals for multiple cell types, and provides insight into the disease-critical cell type. As a proof of concept, we conducted cell-type-aware analyses of breast cancer in 58,648 women and identified 12 transcriptome-wide significant genes using MiXcan compared with only eight genes using conventional approaches. Importantly, MiXcan identified genes with distinct associations in mammary epithelial versus stromal cells, including three new breast cancer susceptibility genes. These findings demonstrate that cell-type-aware transcriptome-wide analyses can reveal new insights into the genetic and cellular etiology of breast cancer and other diseases.

Transforming Diagnosis and Therapeutics Using Cancer Genomics

In past quarter of the century, much has been understood about the genetic variation and abnormal genes that activate cancer in humans. All the cancers somehow possess alterations in the DNA sequence of cancer cell's genome. In present, we are heading toward the era where it is possible to obtain complete genome of the cancer cells for their better diagnosis, categorization and to explore treatment options.

Downregulation of LRRC19 Is Associated with Poor Prognosis in Colorectal Cancer

Objective. Colorectal cancer (CRC) is globally one of the most often diagnosed cancers with high mortality rates. This study aimed to explore novel biomarkers for the diagnosis and prognosis of CRC. Methods. We collected 4 datasets about CRC in GEO and sought differentially expressed genes (DEGs) with GEO2R. Leucine-rich repeat-containing protein 19 (LRRC19) expression was assessed through the Oncomine and TIMER database analyses, which was further confirmed by qRT-PCR of CRC samples. We used online survival analysis tools (GEPIA, PrognoScan, and Kaplan–Meier plotter) to examine the prognostic value of LRRC19 in CRC and other malignancies. GO and KEGG enrichment analyses were employed to explore the biological functions of LRRC19. Finally, we conducted network prediction by STRING and further validation on the GEPIA to discover other molecules that might interact with LRRC19. Results. A total of 21 upregulated and 46 downregulated DEGs were identified from the 4 datasets. The TIMER and Oncomine online analyses showed lower mRNA of LRRC19 in CRC tissues compared with adjacent normal tissues, which was validated by qRT-PCR in CRC patient samples. The survival analysis through the GEPIA and PrognoScan websites revealed that low LRRC19 expression was significantly correlated with poor prognosis in CRC patients. The Kaplan–Meier plotter survival analysis indicated that low LRRC19 expression was significantly associated with the disease progression of patients with ovarian cancer, gastric cancer, breast cancer, and lung cancer. The enrichment analysis suggested that low expression of LRRC19 could be involved in the retinol metabolism and the zymogen granule membrane. Through STRING and GEPIA, it was found that LRRC19 is clearly associated with ZCCHC10, MOB3B, IMMP2L, and TRMT11. Conclusion. LRRC19 mRNA was prominently decreased in human CRC tissues and was significantly associated with shorter survival in CRC patients. LRRC19 might serve as a possible target for early diagnosis and prognosis assessment in CRC.

‘Breast Cancer Resistance Likelihood and Personalized Treatment Through Integrated Multiomics’

In recent times, enormous progress has been made in improving the diagnosis and therapeutic strategies for breast carcinoma, yet it remains the most prevalent cancer and second highest contributor to cancer-related deaths in women. Breast cancer (BC) affects one in eight females globally. In 2018 alone, 1.4 million cases were identified worldwide in postmenopausal women and 645,000 cases in premenopausal females, and this burden is constantly increasing. This shows that still a lot of efforts are required to discover therapeutic remedies for this disease. One of the major clinical complications associated with the treatment of breast carcinoma is the development of therapeutic resistance. Multidrug resistance (MDR) and consequent relapse on therapy are prevalent issues related to breast carcinoma; it is due to our incomplete understanding of the molecular mechanisms of breast carcinoma disease. Therefore, elucidating the molecular mechanisms involved in drug resistance is critical. For management of breast carcinoma, the treatment decision not only depends on the assessment of prognosis factors but also on the evaluation of pathological and clinical factors. Integrated data assessments of these multiple factors of breast carcinoma through multiomics can provide significant insight and hope for making therapeutic decisions. This omics approach is particularly helpful since it identifies the biomarkers of disease progression and treatment progress by collective characterization and quantification of pools of biological molecules within and among the cancerous cells. The scrupulous understanding of cancer and its treatment at the molecular level led to the concept of a personalized approach, which is one of the most significant advancements in modern oncology. Likewise, there are certain genetic and non-genetic tests available for BC which can help in personalized therapy. Genetically inherited risks can be screened for personal predisposition to BC, and genetic changes or variations (mutations) can also be identified to decide on the best treatment. Ultimately, further understanding of BC at the molecular level (multiomics) will define more precise choices in personalized medicine. In this review, we have summarized therapeutic resistance associated with BC and the techniques used for its management.

Functional annotation of breast cancer risk loci: current progress and future directions

Genome-wide association studies coupled with large-scale replication and fine-scale mapping studies have identified more than 150 genomic regions that are associated with breast cancer risk. Here, we review efforts to translate these findings into a greater understanding of disease mechanism. Our review comes in the context of a recently published fine-scale mapping analysis of these regions, which reported 352 independent signals and a total of 13,367 credible causal variants. The vast majority of credible causal variants map to noncoding DNA, implicating regulation of gene expression as the mechanism by which functional variants influence risk. Accordingly, we review methods for defining candidate-regulatory sequences, methods for identifying putative target genes and methods for linking candidate-regulatory sequences to putative target genes. We provide a summary of available data resources and identify gaps in these resources. We conclude that while much work has been done, there is still much to do. There are, however, grounds for optimism; combining statistical data from fine-scale mapping with functional data that are more representative of the normal "at risk" breast, generated using new technologies, should lead to a greater understanding of the mechanisms that influence an individual woman's risk of breast cancer.

Aberrant epigenetic and transcriptional events associated with breast cancer risk

BACKGROUND

Genome-wide association studies have identified several breast cancer susceptibility loci. However, biomarkers for risk assessment are still missing. Here, we investigated cancer-related molecular changes detected in tissues from women at high risk for breast cancer prior to disease manifestation. Disease-free breast tissue cores donated by healthy women (N = 146, median age = 39 years) were processed for both methylome (MethylCap) and transcriptome (Illumina's HiSeq4000) sequencing. Analysis of tissue microarray and primary breast epithelial cells was used to confirm gene expression dysregulation.

RESULTS

Transcriptomic analysis identified 69 differentially expressed genes between women at high and those at average risk of breast cancer (Tyrer-Cuzick model) at FDR < 0.05 and fold change ≥ 2. Majority of the identified genes were involved in DNA damage checkpoint, cell cycle, and cell adhesion. Two genes, FAM83A and NEK2, were overexpressed in tissue sections (FDR < 0.01) and primary epithelial cells (p < 0.05) from high-risk breasts. Moreover, 1698 DNA methylation changes were identified in high-risk breast tissues (FDR < 0.05), partially overlapped with cancer-related signatures, and correlated with transcriptional changes (p < 0.05, r ≤ 0.5). Finally, among the participants, 35 women donated breast biopsies at two time points, and age-related molecular alterations enhanced in high-risk subjects were identified.

CONCLUSIONS

Normal breast tissue from women at high risk of breast cancer bears molecular aberrations that may contribute to breast cancer susceptibility. This study is the first molecular characterization of the true normal breast tissues, and provides an opportunity to investigate molecular markers of breast cancer risk, which may lead to new preventive approaches.

TIGAR-V2: Efficient TWAS tool with nonparametric Bayesian eQTL weights of 49 tissue types from GTEx V8

Standard transcriptome-wide association study (TWAS) methods first train gene expression prediction models using reference transcriptomic data and then test the association between the predicted genetically regulated gene expression and phenotype of interest. Most existing TWAS tools require cumbersome preparation of genotype input files and extra coding to enable parallel computation. To improve the efficiency of TWAS tools, we developed Transcriptome-Integrated Genetic Association Resource V2 (TIGAR-V2), which directly reads Variant Call Format (VCF) files, enables parallel computation, and reduces up to 90% of computation cost (mainly due to loading genotype data) compared to the original version. TIGAR-V2 can train gene expression imputation models using either nonparametric Bayesian Dirichlet process regression (DPR) or Elastic-Net (as used by PrediXcan), perform TWASs using either individual-level or summary-level genome-wide association study (GWAS) data, and implement both burden and variance-component statistics for gene-based association tests. We trained gene expression prediction models by DPR for 49 tissues using Genotype-Tissue Expression (GTEx) V8 by TIGAR-V2 and illustrated the usefulness of these Bayesian cis-expression quantitative trait locus (eQTL) weights through TWASs of breast and ovarian cancer utilizing public GWAS summary statistics. We identified 88 and 37 risk genes, respectively, for breast and ovarian cancer, most of which are either known or near previously identified GWAS (∼95%) or TWAS (∼40%) risk genes and three novel independent TWAS risk genes with known functions in carcinogenesis. These findings suggest that TWASs can provide biological insight into the transcriptional regulation of complex diseases. The TIGAR-V2 tool, trained Bayesian cis-eQTL weights, and linkage disequilibrium (LD) information from GTEx V8 are publicly available, providing a useful resource for mapping risk genes of complex diseases.

Joint Analysis of Genetic Correlation, Mendelian Randomization and Colocalization Highlights the Bi-Directional Causal Association Between Hypothyroidism and Primary Biliary Cirrhosis

Background: Hypothyroidism and primary biliary cirrhosis (PBC) are often co-existed in observational epidemiological studies. However, the causal relationship between them remains unclear.

Methods: Genetic correlation, Mendelian randomization (MR) and colocalization analysis were combined to assess the potential causal association between hypothyroidism and PBC by using summary statistics from large-scale genome-wide association studies. Various sensitivity analyses had been conducted to assess the robustness and the consistency of the findings.

Results: The linkage disequilibrium score regression demonstrated significant evidence of shared genetic architecture between hypothyroidism and PBC, with the genetic correlation estimated to be 0.117 (p = 0.006). The OR of hypothyroidism on PBC was 1.223 (95% CI, 1.072–1.396; p = 2.76 × 10⁻³) in MR analysis with inverse variance weighted (IVW) method. More importantly, the results from other 7MR methods with different model assumptions, were almost identical with that of IVW, suggesting the findings were robust and convincing. On the other hand, PBC was also causally associated with hypothyroidism (OR, 1.049; 95% CI, 1.010–1.089; p = 0.012), and, again, similar results can also be obtained from other MR methods. Various sensitivity analyses regarding the outlier detection and leave-one-out analysis were also performed. Besides, colocalization analysis suggested that there existed shared causal variants between hypothyroidism and PBC, further highlighting the robustness of the results.

Conclusion: Our results suggest evidence for the bi-directional causal association between hypothyroidism and PBC, which may provide insights into the etiology of hypothyroidism and PBC as well as inform prevention and intervention strategies directed toward both diseases.

Nuclear Receptor-Mediated Metabolic Reprogramming and the Impact on HR+ Breast Cancer

Simple Summary

Breast cancer is the most commonly diagnosed and second leading cause of cancer-related deaths in women in the United States, with hormone receptor positive (HR+) tumors representing more than two-thirds of new cases. Recent evidence has indicated that dysregulation of multiple metabolic programs, which can be driven through nuclear receptor activity, is essential for tumor genesis, progression, therapeutic resistance and metastasis. This study will review the current advances in our understanding of the impact and implication of altered metabolic processes driven by nuclear receptors, including hormone-dependent signaling, on HR+ breast cancer.

Abstract

Metabolic reprogramming enables cancer cells to adapt to the changing microenvironment in order to maintain metabolic energy and to provide the necessary biological macromolecules required for cell growth and tumor progression. While changes in tumor metabolism have been long recognized as a hallmark of cancer, recent advances have begun to delineate the mechanisms that modulate metabolic pathways and the consequence of altered signaling on tumorigenesis. This is particularly evident in hormone receptor positive (HR+) breast cancers which account for approximately 70% of breast cancer cases. Emerging evidence indicates that HR+ breast tumors are dependent on multiple metabolic processes for tumor progression, metastasis, and therapeutic resistance and that changes in metabolic programs are driven, in part, by a number of key nuclear receptors including hormone-dependent signaling. In this review, we discuss the mechanisms and impact of hormone receptor mediated metabolic reprogramming on HR+ breast cancer genesis and progression as well as the therapeutic implications of these metabolic processes in this disease.

Genome-Wide Identification of Cis-acting Expression QTLs in Large Yellow Croaker

The large yellow croaker is one of the largest marine economic fish in China with the farming yield about 30 tons per year. The genetic selection for growth and disease resistance has been performed in recent years. The identification of trait-associated molecular makers, causative variants, and causative genes is helpful for genetic selection in large yellow croaker. It has been discovered that most of polygenic traits are controlled with multiple genes via regulatory variant, and GWAS-identified loci are enriched in the regulatory variant. Cis-acting eQTL is a widespread regulatory variant that controls the expression of nearby gene. We herein take advantage of RNA-seq and whole genome sequencing technique to identify genome-wide eQTLs in liver tissue for large yellow croaker; a forward selection routine is applied for identification of multiple eQTLs. To fine map causative mutation for each eQTL, a credible set is built to confine causative variants. Totally, 2427 eQTLs have been identified, 69.7% (1,691/2,427) of them are primary eQTL signals, and the remaining are secondary signals, many functional important target genes have been discovered. We highlight several functional pivotal genes including SMC3, TUSC3, TITIN, MCPH1, and MDHC, in which the expression of MCPH1 is regulated by two eQTLs; the distance of these eQTLs from target genes is symmetrically distributed, 25.5% of them are within 1 Mb region from target genes, whereas 74.5% of them are between 1 and 2 Mb regions; most of the identified eQTL has been well resolved, and 19.3 (469/2427) of eQTL have the size of credible set (the number of variants in credible set) less than 50.

Multi-omics analysis to identify susceptibility genes for colorectal cancer

Most genetic variants for colorectal cancer (CRC) identified in genome-wide association studies (GWAS) are located in intergenic regions, implying pathogenic dysregulations of gene expression. However, comprehensive assessments of target genes in CRC remain to be explored. We conducted a multi-omics analysis using transcriptome and/or DNA methylation data from the Genotype-Tissue Expression, The Cancer Genome Atlas and the Colonomics projects. We identified 116 putative target genes for 45 GWAS-identified variants. Using summary-data-based Mendelian randomization approach (SMR), we demonstrated that the CRC susceptibility for 29 out of the 45 CRC variants may be mediated by cis-effects on gene regulation. At a cutoff of the Bonferroni-corrected PSMR < 0.05, we determined 66 putative susceptibility genes, including 39 genes that have not been previously reported. We further performed in vitro assays for two selected genes, DIP2B and SFMBT1, and provide functional evidence that they play a vital role in colorectal carcinogenesis via disrupting cell behavior, including migration, invasion and epithelial-mesenchymal transition. Our study reveals a large number of putative novel susceptibility genes and provides additional insight into the underlying mechanisms for CRC genetic risk loci.

Integrating GWAS and eQTL to predict genes and pathways for non-syndromic cleft lip with or without palate

OBJECTIVE

To explore susceptibility genes and pathways for non-syndromic cleft lip with or without cleft palate (NSCL/P).

MATERIALS AND METHODS

Two genome-wide association studies (GWAS) datasets, including 858 NSCL/P cases and 1,248 controls, were integrated with expression quantitative trait loci (eQTL) dataset identified by Genotype-Tissue Expression (GTEx) project in whole-blood samples. The expression of the candidate genes in mouse orofacial development was inquired from FaceBase. Protein-protein interaction (PPI) network was visualized to identify protein functions. Go and KEGG pathway analyses were performed to explore the underlying risk pathways.

RESULTS

A total of 233 eQTL single-nucleotide polymorphisms (SNPs) in 432 candidate genes were identified to be associated with the risk of NSCL/P. One hundred and eighty-three susceptible genes were expressed in mouse orofacial development according to FaceBase. PPI network analysis highlighted that these genes involved in ubiquitin-mediated proteolysis (KCTD7, ASB1, UBOX5, ANAPC4) and DNA synthesis (XRCC3, RFC3, KAT5, RHNO1) were associated with the risk of NSCL/P. GO and KEGG pathway analyses revealed that the fatty acid metabolism pathway (ACADL, HSD17B12, ACSL5, PPT1, MCAT) played an important role in the development of NSCL/P.

CONCLUSIONS

Our results identified novel susceptibility genes and pathways associated with the development of NSCL/P.

A Transcriptome-Wide Association Study Identifies Novel Candidate Susceptibility Genes for Pancreatic Cancer

BACKGROUND

Although 20 pancreatic cancer susceptibility loci have been identified through genome-wide association studies in individuals of European ancestry, much of its heritability remains unexplained and the genes responsible largely unknown.

METHODS

To discover novel pancreatic cancer risk loci and possible causal genes, we performed a pancreatic cancer transcriptome-wide association study in Europeans using three approaches: FUSION, MetaXcan, and Summary-MulTiXcan. We integrated genome-wide association studies summary statistics from 9040 pancreatic cancer cases and 12 496 controls, with gene expression prediction models built using transcriptome data from histologically normal pancreatic tissue samples (NCI Laboratory of Translational Genomics [n = 95] and Genotype-Tissue Expression v7 [n = 174] datasets) and data from 48 different tissues (Genotype-Tissue Expression v7, n = 74-421 samples).

RESULTS

We identified 25 genes whose genetically predicted expression was statistically significantly associated with pancreatic cancer risk (false discovery rate < .05), including 14 candidate genes at 11 novel loci (1p36.12: CELA3B; 9q31.1: SMC2, SMC2-AS1; 10q23.31: RP11-80H5.9; 12q13.13: SMUG1; 14q32.33: BTBD6; 15q23: HEXA; 15q26.1: RCCD1; 17q12: PNMT, CDK12, PGAP3; 17q22: SUPT4H1; 18q11.22: RP11-888D10.3; and 19p13.11: PGPEP1) and 11 at six known risk loci (5p15.33: TERT, CLPTM1L, ZDHHC11B; 7p14.1: INHBA; 9q34.2: ABO; 13q12.2: PDX1; 13q22.1: KLF5; and 16q23.1: WDR59, CFDP1, BCAR1, TMEM170A). The association for 12 of these genes (CELA3B, SMC2, and PNMT at novel risk loci and TERT, CLPTM1L, INHBA, ABO, PDX1, KLF5, WDR59, CFDP1, and BCAR1 at known loci) remained statistically significant after Bonferroni correction.

CONCLUSIONS

By integrating gene expression and genotype data, we identified novel pancreatic cancer risk loci and candidate functional genes that warrant further investigation.

The regulatory genome constrains protein sequence evolution: implications for the search for disease-associated genes

The development of explanatory models of protein sequence evolution has broad implications for our understanding of cellular biology, population history, and disease etiology. Here we analyze the GTEx transcriptome resource to quantify the effect of the transcriptome on protein sequence evolution in a multi-tissue framework. We find substantial variation among the central nervous system tissues in the effect of expression variance on evolutionary rate, with highly variable genes in the cortex showing significantly greater purifying selection than highly variable genes in subcortical regions (Mann-Whitney U p = 1.4 × 10^-4). The remaining tissues cluster in observed expression correlation with evolutionary rate, enabling evolutionary analysis of genes in diverse physiological systems, including digestive, reproductive, and immune systems. Importantly, the tissue in which a gene attains its maximum expression variance significantly varies (p = 5.55 × 10^-284) with evolutionary rate, suggesting a tissue-anchored model of protein sequence evolution. Using a large-scale reference resource, we show that the tissue-anchored model provides a transcriptome-based approach to predicting the primary affected tissue of developmental disorders. Using gradient boosted regression trees to model evolutionary rate under a range of model parameters, selected features explain up to 62% of the variation in evolutionary rate and provide additional support for the tissue model. Finally, we investigate several methodological implications, including the importance of evolutionary-rate-aware gene expression imputation models using genetic data for improved search for disease-associated genes in transcriptome-wide association studies. Collectively, this study presents a comprehensive transcriptome-based analysis of a range of factors that may constrain molecular evolution and proposes a novel framework for the study of gene function and disease mechanism.

Transcriptome-wide association study of breast cancer risk by estrogen-receptor status

Previous transcriptome-wide association studies (TWAS) have identified breast cancer risk genes by integrating data from expression quantitative loci and genome-wide association studies (GWAS), but analyses of breast cancer subtype-specific associations have been limited. In this study, we conducted a TWAS using gene expression data from GTEx and summary statistics from the hitherto largest GWAS meta-analysis conducted for breast cancer overall, and by estrogen receptor subtypes (ER+ and ER-). We further compared associations with ER+ and ER- subtypes, using a case-only TWAS approach. We also conducted multigene conditional analyses in regions with multiple TWAS associations. Two genes, STXBP4 and HIST2H2BA, were specifically associated with ER+ but not with ER- breast cancer. We further identified 30 TWAS-significant genes associated with overall breast cancer risk, including four that were not identified in previous studies. Conditional analyses identified single independent breast-cancer gene in three of six regions harboring multiple TWAS-significant genes. Our study provides new information on breast cancer genetics and biology, particularly about genomic differences between ER+ and ER- breast cancer.

ANKLE1 N6 -Methyladenosine-related variant is associated with colorectal cancer risk by maintaining the genomic stability

The N6 -Methyladenosine (m6 A) modification plays an important role in many biological processes, especially tumor development. However, little is still known about how it affects colorectal cancer (CRC) carcinogenesis. Here, we first systematically investigate the association of variants related to m6 A modification with the CRC risk in 1,062 CRC cases and 2,184 controls by using our exome-wide association data and followed by two replication sets including 7,341 CRC cases and 7,902 controls. The variant rs8100241 located in ANKLE1 was significantly associated with CRC risk (odds ratio = 0.88, 95% confidence interval = 0.84-0.92, p = 4.85 × 10-8 ) in 8,403 cases and 10,086 controls. This variant was previously identified to be associated with the susceptibility of breast cancer with BRCA1 mutation triple negative breast cancer. Further functional analysis indicated that overexpression of the rs8100241[A] allele significantly increased the ANKLE1 m6 A level and facilitated the ANKLE1 protein expression compared to that of rs8100241[G] allele. We further found the ANKLE1 m6 A modification was catalyzed by the "writer" complex (METTL3, METTL14, or WTAP) and recognized by the "reader" YTHDF1. Mechanistically, we found that the ANKLE1 functions as a potential tumor suppressor that inhibits cell proliferation and facilitates the genomic stability. An elevated frequency of micronucleated cells, increased cell proliferation, and colony formation ability were observed when ANKLE1 knockdown. Our study illustrated that the germline missense variant can increase CRC risk by influencing ANKLE1 m6 A level, highlighting a clinical potential of variants-associated m6 A modification as a risk marker for CRC prevention.

Prognostic Value of a Stemness Index-Associated Signature in Primary Lower-Grade Glioma

Objective

As a prevalent and infiltrative cancer type of the central nervous system, the prognosis of lower-grade glioma (LGG) in adults is highly heterogeneous. Recent evidence has demonstrated the prognostic value of the mRNA expression-based stemness index (mRNAsi) in LGG. Our aim was to develop a stemness index-based signature (SI-signature) for risk stratification and survival prediction.

Methods

Differentially expressed genes (DEGs) between LGG in the Cancer Genome Atlas (TCGA) and normal brain tissue samples from the Genotype-Tissue Expression (GTEx) project were screened out, and the weighted gene correlation network analysis (WGCNA) was employed to identify the mRNAsi-related gene sets. Meanwhile, the Gene Ontology and Kyoto Encyclopedia of Genes and Genomes enrichment analyses were performed for the functional annotation of the key genes. ESTIMATE was used to calculate tumor purity for acquiring the correct mRNAsi. Differences in overall survival (OS) between the high and low mRNAsi (corrected mRNAsi) groups were compared using the Kaplan Meier analysis. By combining the Lasso regression with univariate and multivariate Cox regression, the SI-signature was constructed and validated using the Chinese Glioma Genome Atlas (CGGA).

Results

There was a significant difference in OS between the high and low mRNAsi groups, which was also observed in the two corrected mRNAsi groups. Based on threshold limits, 86 DEGs were most significantly associated with mRNAsi via WGCNA. Seven genes (ADAP2, ALOX5AP, APOBEC3C, FCGRT, GNG5, LRRC25, and SP100) were selected to establish a risk signature for primary LGG. The ROC curves showed a fair performance in survival prediction in both the TCGA and the CGGA validation cohorts. Univariate and multivariate Cox regression revealed that the risk group was an independent prognostic factor in primary LGG. The nomogram was developed based on clinical parameters integrated with the risk signature, and its accuracy for predicting 3- and 5-years survival was assessed by the concordance index, the area under the curve of the time-dependent receiver operating characteristics curve, and calibration curves.

Conclusion

The SI-signature with seven genes could serve as an independent predictor, and suggests the importance of stemness features in risk stratification and survival prediction in primary LGG.

A framework for transcriptome-wide association studies in breast cancer in diverse study populations

BACKGROUND

The relationship between germline genetic variation and breast cancer survival is largely unknown, especially in understudied minority populations who often have poorer survival. Genome-wide association studies (GWAS) have interrogated breast cancer survival but often are underpowered due to subtype heterogeneity and clinical covariates and detect loci in non-coding regions that are difficult to interpret. Transcriptome-wide association studies (TWAS) show increased power in detecting functionally relevant loci by leveraging expression quantitative trait loci (eQTLs) from external reference panels in relevant tissues. However, ancestry- or race-specific reference panels may be needed to draw correct inference in ancestrally diverse cohorts. Such panels for breast cancer are lacking.

RESULTS

We provide a framework for TWAS for breast cancer in diverse populations, using data from the Carolina Breast Cancer Study (CBCS), a population-based cohort that oversampled black women. We perform eQTL analysis for 406 breast cancer-related genes to train race-stratified predictive models of tumor expression from germline genotypes. Using these models, we impute expression in independent data from CBCS and TCGA, accounting for sampling variability in assessing performance. These models are not applicable across race, and their predictive performance varies across tumor subtype. Within CBCS (N = 3,828), at a false discovery-adjusted significance of 0.10 and stratifying for race, we identify associations in black women near AURKA, CAPN13, PIK3CA, and SERPINB5 via TWAS that are underpowered in GWAS.

CONCLUSIONS

We show that carefully implemented and thoroughly validated TWAS is an efficient approach for understanding the genetics underpinning breast cancer outcomes in diverse populations.

Bioinformatic profiling of prognosis-related genes in the breast cancer immune microenvironment

In the microenvironment of breast cancer, immune cell infiltration is associated with an improved prognosis. To identify immune-related prognostic markers and therapeutic targets, we determined the lymphocyte-specific kinase (LCK) metagene scores of samples from breast cancer patients in The Cancer Genome Atlas. The LCK metagene score correlated highly with other immune-related scores, as well as with the clinical stage, prognosis and tumor suppressor gene mutation status (BRCA2, TP53, PTEN) of patients in the four breast cancer subtypes. A weighted gene co-expression network analysis was performed to detect representative genes from LCK metagene-related gene modules. In two of these modules, the levels of the co-expressed genes correlated highly with LCK metagene levels, so we conducted an enrichment analysis to discover their functions. We also identified differentially expressed genes in samples with high and low LCK metagene scores. By examining the overlapping results from these analyses, we obtained 115 genes, and found that 22 of them were independent predictors of overall survival in breast cancer patients. These genes were validated for their prognostic and diagnostic value with external data sets and paired tumor and non-tumor tissues. The genes identified herein could serve as diagnostic/prognostic markers and immune-related therapeutic targets in breast cancer.

LINC01419 promotes cell proliferation and metastasis in lung adenocarcinoma via sponging miR-519b-3p to up-regulate RCCD1

Lung adenocarcinoma (LUAD) is one of the most common type of lung cancer notorious for the high incidence and mortality around the world. Long non-coding RNAs (LncRNAs) are defined as a class of RNAs with length more than 200 nucleotides. Mounting studies have proved that lncRNAs are related closely to incidence of diseases and play crucial roles in cancer progression. Although LINC01419 has been studied in several cancers, the function and mechanism of LINC01419 in LUAD remains a mystery. Our findings showed that LINC01419 level was high in LUAD cells, and LINC01419 knockdown inhibited carcinogenesis via suppressing cell proliferation, migration as well as invasion. Moreover, bioinformatics prediction, luciferase reporter experiments and RIP assay were used to confirm miR-519-3p was sequestered and negatively regulated by LINC01419. Subsequently, RCCD1 was identified as a miR-519b-3p target and had inverse relationship with miR-519b-3p. LINC01419 was positively related to RCCD1. Furthermore, rescue assays confirmed that miR-519b-3p inhibitor or RCCD1 overexpression could neutralize the effect of LINC01419 silenced in cell proliferation, migration and invasion. Taken together, all the results indicated that LINC01419 exhibited oncogenic behaviors LUAD via binding to miR-519b-3p to enhance the expression of RCCD1, manifesting the underlying therapy values of LINC01419 in LUAD.

The Need for Multi-Omics Biomarker Signatures in Precision Medicine

Recent advances in omics technologies have led to unprecedented efforts characterizing the molecular changes that underlie the development and progression of a wide array of complex human diseases, including cancer. As a result, multi-omics analyses—which take advantage of these technologies in genomics, transcriptomics, epigenomics, proteomics, metabolomics, and other omics areas—have been proposed and heralded as the key to advancing precision medicine in the clinic. In the field of precision oncology, genomics approaches, and, more recently, other omics analyses have helped reveal several key mechanisms in cancer development, treatment resistance, and recurrence risk, and several of these findings have been implemented in clinical oncology to help guide treatment decisions. However, truly integrated multi-omics analyses have not been applied widely, preventing further advances in precision medicine. Additional efforts are needed to develop the analytical infrastructure necessary to generate, analyze, and annotate multi-omics data effectively to inform precision medicine-based decision-making.

Genome-wide association and transcriptome studies identify target genes and risk loci for breast cancer

Genome-wide association studies (GWAS) have identified more than 170 breast cancer susceptibility loci. Here we hypothesize that some risk-associated variants might act in non-breast tissues, specifically adipose tissue and immune cells from blood and spleen. Using expression quantitative trait loci (eQTL) reported in these tissues, we identify 26 previously unreported, likely target genes of overall breast cancer risk variants, and 17 for estrogen receptor (ER)-negative breast cancer, several with a known immune function. We determine the directional effect of gene expression on disease risk measured based on single and multiple eQTL. In addition, using a gene-based test of association that considers eQTL from multiple tissues, we identify seven (and four) regions with variants associated with overall (and ER-negative) breast cancer risk, which were not reported in previous GWAS. Further investigation of the function of the implicated genes in breast and immune cells may provide insights into the etiology of breast cancer.

Cross-Cancer Pleiotropic Associations with Lung Cancer Risk in African Americans

BACKGROUND

Identifying genetic variants with pleiotropic associations across multiple cancers can reveal shared biologic pathways. Prior pleiotropic studies have primarily focused on European-descent individuals. Yet population-specific genetic variation can occur, and potential pleiotropic associations among diverse racial/ethnic populations could be missed. We examined cross-cancer pleiotropic associations with lung cancer risk in African Americans.

METHODS

We conducted a pleiotropic analysis among 1,410 African American lung cancer cases and 2,843 controls. We examined 36,958 variants previously associated (or in linkage disequilibrium) with cancer in prior genome-wide association studies. Logistic regression analyses were conducted, adjusting for age, sex, global ancestry, study site, and smoking status.

RESULTS

We identified three novel genomic regions significantly associated (FDR-corrected P <0.10) with lung cancer risk (rs336958 on 5q14.3, rs7186207 on 16q22.2, and rs11658063 on 17q12). On chromosome16q22.2, rs7186207 was significantly associated with reduced risk [OR = 0.43; 95% confidence interval (CI), 0.73-0.89], and functional annotation using GTEx showed rs7186207 modifies DHODH gene expression. The minor allele at rs336958 on 5q14.3 was associated with increased lung cancer risk (OR = 1.47; 95% CI, 1.22-1.78), whereas the minor allele at rs11658063 on 17q12 was associated with reduced risk (OR = 0.80; 95% CI, 0.72-0.90).

CONCLUSIONS

We identified novel associations on chromosomes 5q14.3, 16q22.2, and 17q12, which contain HNF1B, DHODH, and HAPLN1 genes, respectively. SNPs within these regions have been previously associated with multiple cancers. This is the first study to examine cross-cancer pleiotropic associations for lung cancer in African Americans.

IMPACT

Our findings demonstrate novel cross-cancer pleiotropic associations with lung cancer risk in African Americans.

A statistical framework for cross-tissue transcriptome-wide association analysis

Transcriptome-wide association analysis is a powerful approach to studying the genetic architecture of complex traits. A key component of this approach is to build a model to impute gene expression levels from genotypes by using samples with matched genotypes and gene expression data in a given tissue. However, it is challenging to develop robust and accurate imputation models with a limited sample size for any single tissue. Here, we first introduce a multi-task learning method to jointly impute gene expression in 44 human tissues. Compared with single-tissue methods, our approach achieved an average of 39% improvement in imputation accuracy and generated effective imputation models for an average of 120% more genes. We describe a summary-statistic-based testing framework that combines multiple single-tissue associations into a powerful metric to quantify the overall gene-trait association. We applied our method, called UTMOST (unified test for molecular signatures), to multiple genome-wide-association results and demonstrate its advantages over single-tissue strategies.

Cancer expression quantitative trait loci (eQTLs) can be determined from heterogeneous tumor gene expression data by modeling variation in tumor purity

Expression quantitative trait loci (eQTLs) identified using tumor gene expression data could affect gene expression in cancer cells, tumor-associated normal cells, or both. Here, we have demonstrated a method to identify eQTLs affecting expression in cancer cells by modeling the statistical interaction between genotype and tumor purity. Only one third of breast cancer risk variants, identified as eQTLs from a conventional analysis, could be confidently attributed to cancer cells. The remaining variants could affect cells of the tumor microenvironment, such as immune cells and fibroblasts. Deconvolution of tumor eQTLs will help determine how inherited polymorphisms influence cancer risk, development, and treatment response.

Elucidating the Underlying Functional Mechanisms of Breast Cancer Susceptibility Through Post-GWAS Analyses

Genome-wide association studies (GWAS) have identified more than 170 single nucleotide polymorphisms (SNPs) associated with the susceptibility to breast cancer. Together, these SNPs explain 18% of the familial relative risk, which is estimated to be nearly half of the total familial breast cancer risk that is collectively explained by low-risk susceptibility alleles. An important aspect of this success has been the access to large sample sizes through collaborative efforts within the Breast Cancer Association Consortium (BCAC), but also collaborations between cancer association consortia. Despite these achievements, however, understanding of each variant's underlying mechanism and how these SNPs predispose women to breast cancer remains limited and represents a major challenge in the field, particularly since the vast majority of the GWAS-identified SNPs are located in non-coding regions of the genome and are merely tags for the causal variants. In recent years, fine-scale mapping studies followed by functional evaluation of putative causal variants have begun to elucidate the biological function of several GWAS-identified variants. In this review, we discuss the findings and lessons learned from these post-GWAS analyses of 22 risk loci. Identifying the true causal variants underlying breast cancer susceptibility and their function not only provides better estimates of the explained familial relative risk thereby improving polygenetic risk scores (PRSs), it also increases our understanding of the biological mechanisms responsible for causing susceptibility to breast cancer. This will facilitate the identification of further breast cancer risk alleles and the development of preventive medicine for those women at increased risk for developing the disease.

Complex polymorphisms in endocytosis genes suggest alpha-cyclodextrin as a treatment for breast cancer

Most breast cancer deaths are caused by metastasis and treatment options beyond radiation and cytotoxic drugs, which have severe side effects, and hormonal treatments, which are or become ineffective for many patients, are urgently needed. This study reanalyzed existing data from three genome-wide association studies (GWAS) using a novel computational biostatistics approach (muGWAS), which had been validated in studies of 600-2000 subjects in epilepsy and autism. MuGWAS jointly analyzes several neighboring single nucleotide polymorphisms while incorporating knowledge about genetics of heritable diseases into the statistical method and about GWAS into the rules for determining adaptive genome-wide significance. Results from three independent GWAS of 1000-2000 subjects each, which were made available under the National Institute of Health's "Up For A Challenge" (U4C) project, not only confirmed cell-cycle control and receptor/AKT signaling, but, for the first time in breast cancer GWAS, also consistently identified many genes involved in endo-/exocytosis (EEC), most of which had already been observed in functional and expression studies of breast cancer. In particular, the findings include genes that translocate (ATP8A1, ATP8B1, ANO4, ABCA1) and metabolize (AGPAT3, AGPAT4, DGKQ, LPPR1) phospholipids entering the phosphatidylinositol cycle, which controls EEC. These novel findings suggest scavenging phospholipids as a novel intervention to control local spread of cancer, packaging of exosomes (which prepare distant microenvironment for organ-specific metastases), and endocytosis of β1 integrins (which are required for spread of metastatic phenotype and mesenchymal migration of tumor cells). Beta-cyclodextrins (βCD) have already been shown to be effective in in vitro and animal studies of breast cancer, but exhibits cholesterol-related ototoxicity. The smaller alpha-cyclodextrins (αCD) also scavenges phospholipids, but cannot fit cholesterol. An in-vitro study presented here confirms hydroxypropyl (HP)-αCD to be twice as effective as HPβCD against migration of human cells of both receptor negative and estrogen-receptor positive breast cancer. If the previous successful animal studies with βCDs are replicated with the safer and more effective αCDs, clinical trials of adjuvant treatment with αCDs are warranted. Ultimately, all breast cancer are expected to benefit from treatment with HPαCD, but women with triple-negative breast cancer (TNBC) will benefit most, because they have fewer treatment options and their cancer advances more aggressively.

A transcriptome-wide association study of 229,000 women identifies new candidate susceptibility genes for breast cancer

The breast cancer risk variants identified in genome-wide association studies explain only a small fraction of the familial relative risk, and the genes responsible for these associations remain largely unknown. To identify novel risk loci and likely causal genes, we performed a transcriptome-wide association study evaluating associations of genetically predicted gene expression with breast cancer risk in 122,977 cases and 105,974 controls of European ancestry. We used data from the Genotype-Tissue Expression Project to establish genetic models to predict gene expression in breast tissue and evaluated model performance using data from The Cancer Genome Atlas. Of the 8,597 genes evaluated, significant associations were identified for 48 at a Bonferroni-corrected threshold of P < 5.82 × 10-6, including 14 genes at loci not yet reported for breast cancer. We silenced 13 genes and showed an effect for 11 on cell proliferation and/or colony-forming efficiency. Our study provides new insights into breast cancer genetics and biology.

Common Genetic Variation and Breast Cancer Risk-Past, Present, and Future

Breast cancer is the most common cancer among women in the United States, with up to 30% of those diagnosed displaying a family history of breast cancer. To date, 18% of the familial risk of breast cancer can be explained by SNPs. This review summarizes the discovery of risk-associated SNPs using candidate gene and genome-wide association studies (GWAS), including discovery and replication in large collaborative efforts such as The Collaborative Oncologic Gene-environment Study and OncoArray. We discuss the evolution of GWAS studies, efforts to discover additional SNPs, and methods for identifying causal variants. We summarize findings associated with overall breast cancer, pathologic subtypes, and mutation carriers (BRCA1, BRCA2, and CHEK2). In addition, we summarize the development of polygenic risk scores (PRS) using the risk-associated SNPs and show how PRS can contribute to estimation of individual risks for developing breast cancer. Cancer Epidemiol Biomarkers Prev; 27(4); 380-94. ©2018 AACRSee all articles in this CEBP Focus section, "Genome-Wide Association Studies in Cancer."

Exome sequencing and case-control analyses identify RCC1 as a candidate breast cancer susceptibility gene

Breast cancer is a genetic disease but the known genes explain a minority of cases. To elucidate the molecular basis of breast cancer in the Tunisian population, we performed exome sequencing on six BRCA1/BRCA2 mutation-negative patients with familial breast cancer and identified a novel frameshift mutation in RCC1, encoding the Regulator of Chromosome Condensation 1. Subsequent genotyping detected the 19-bp deletion in additional 5 out of 153 (3%) breast cancer patients but in none of 400 female controls (p = 0.0015). The deletion was enriched in patients with a positive family history (5%, p = 0.0009) and co-segregated with breast cancer in the initial pedigree. The mutant allele was lost in 4/6 breast tumors from mutation carriers which may be consistent with the hypothesis that RCC1 dysfunction provides a selective disadvantage at the stage of tumor progression. In summary, we propose RCC1 as a likely breast cancer susceptibility gene in the Tunisian population.