Identification of independent association signals and putative functional variants for breast cancer risk through fine-scale mapping of the 12p11 locus

Enhancing Disease Risk Gene Discovery by Integrating Transcription Factor-Linked Trans-located Variants into Transcriptome-Wide Association Analyses

Transcriptome-wide association studies (TWAS) have been successful in identifying disease susceptibility genes by integrating cis-variants predicted gene expression with genome-wide association studies (GWAS) data. However, trans-located variants for predicting gene expression remain largely unexplored. Here, we introduce transTF-TWAS, which incorporates transcription factor (TF)-linked trans-located variants to enhance model building. Using data from the Genotype-Tissue Expression project, we predict gene expression and alternative splicing and applied these models to large GWAS datasets for breast, prostate, and lung cancers. We demonstrate that transTF-TWAS outperforms other existing TWAS approaches in both constructing gene prediction models and identifying disease-associated genes, as evidenced by simulations and real data analysis. Our transTF-TWAS approach significantly contributes to the discovery of disease risk genes. Findings from this study have shed new light on several genetically driven key regulators and their associated regulatory networks underlying disease susceptibility.

Regulatory network and targeted interventions for CCDC family in tumor pathogenesis

CCDC (coiled-coil domain-containing) is a coiled helix domain that exists in natural proteins. There are about 180 CCDC family genes, encoding proteins that are involved in intercellular transmembrane signal transduction and genetic signal transcription, among other functions. Alterations in expression, mutation, and DNA promoter methylation of CCDC family genes have been shown to be associated with the pathogenesis of many diseases, including primary ciliary dyskinesia, infertility, and tumors. In recent studies, CCDC family genes have been found to be involved in regulation of growth, invasion, metastasis, chemosensitivity, and other biological behaviors of malignant tumor cells in various cancer types, including nasopharyngeal carcinoma, lung cancer, colorectal cancer, and thyroid cancer. In this review, we summarize the involvement of CCDC family genes in tumor pathogenesis and the relevant upstream and downstream molecular mechanisms. In addition, we summarize the potential of CCDC family genes as tumor therapy targets. The findings discussed here help us to further understand the role and the therapeutic applications of CCDC family genes in tumors.

Integrating transcription factor occupancy with transcriptome-wide association analysis identifies susceptibility genes in human cancers

Transcriptome-wide association studies (TWAS) have successfully discovered many putative disease susceptibility genes. However, TWAS may suffer from inaccuracy of gene expression predictions due to inclusion of non-regulatory variants. By integrating prior knowledge of susceptible transcription factor occupied elements, we develop sTF-TWAS and demonstrate that it outperforms existing TWAS approaches in both simulation and real data analyses. Under the sTF-TWAS framework, we build genetic models to predict alternative splicing and gene expression in normal breast, prostate and lung tissues from the Genotype-Tissue Expression project and apply these models to data from large genome-wide association studies (GWAS) conducted among European-ancestry populations. At Bonferroni-corrected P < 0.05, we identify 354 putative susceptibility genes for these cancers, including 189 previously unreported in GWAS loci and 45 in loci unreported by GWAS. These findings provide additional insight into the genetic susceptibility of human cancers. Additionally, we show the generalizability of the sTF-TWAS on non-cancer diseases.

PTHrP induces STAT5 activation, secretory differentiation and accelerates mammary tumor development

BACKGROUND

Parathyroid hormone-related protein (PTHrP) is required for embryonic breast development and has important functions during lactation, when it is produced by alveolar epithelial cells and secreted into the maternal circulation to mobilize skeletal calcium used for milk production. PTHrP is also produced by breast cancers, and GWAS studies suggest that it influences breast cancer risk. However, the exact functions of PTHrP in breast cancer biology remain unsettled.

METHODS

We developed a tetracycline-regulated, MMTV (mouse mammary tumor virus)-driven model of PTHrP overexpression in mammary epithelial cells (Tet-PTHrP mice) and bred these mice with the MMTV-PyMT (polyoma middle tumor-antigen) breast cancer model to analyze the impact of PTHrP overexpression on normal mammary gland biology and in breast cancer progression.

RESULTS

Overexpression of PTHrP in luminal epithelial cells caused alveolar hyperplasia and secretory differentiation of the mammary epithelium with milk production. This was accompanied by activation of Stat5 and increased expression of E74-like factor-5 (Elf5) as well as a delay in post-lactation involution. In MMTV-PyMT mice, overexpression of PTHrP (Tet-PTHrP;PyMT mice) shortened tumor latency and accelerated tumor growth, ultimately reducing overall survival. Tumors overproducing PTHrP also displayed increased expression of nuclear pSTAT5 and Elf5, increased expression of markers of secretory differentiation and milk constituents, and histologically resembled secretory carcinomas of the breast. Overexpression of PTHrP within cells isolated from tumors, but not PTHrP exogenously added to cell culture media, led to activation of STAT5 and milk protein gene expression. In addition, neither ablating the Type 1 PTH/PTHrP receptor (PTH1R) in epithelial cells nor treating Tet-PTHrP;PyMT mice with an anti-PTH1R antibody prevented secretory differentiation or altered tumor latency. These data suggest that PTHrP acts in a cell-autonomous, intracrine manner. Finally, expression of PTHrP in human breast cancers is associated with expression of genes involved in milk production and STAT5 signaling.

CONCLUSIONS

Our study suggests that PTHrP promotes pathways leading to secretory differentiation and proliferation in both normal mammary epithelial cells and in breast tumor cells.

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.

Genome-wide and transcriptome-wide association studies of mammographic density phenotypes reveal novel loci

BACKGROUND

Mammographic density (MD) phenotypes, including percent density (PMD), area of dense tissue (DA), and area of non-dense tissue (NDA), are associated with breast cancer risk. Twin studies suggest that MD phenotypes are highly heritable. However, only a small proportion of their variance is explained by identified genetic variants.

METHODS

We conducted a genome-wide association study, as well as a transcriptome-wide association study (TWAS), of age- and BMI-adjusted DA, NDA, and PMD in up to 27,900 European-ancestry women from the MODE/BCAC consortia.

RESULTS

We identified 28 genome-wide significant loci for MD phenotypes, including nine novel signals (5q11.2, 5q14.1, 5q31.1, 5q33.3, 5q35.1, 7p11.2, 8q24.13, 12p11.2, 16q12.2). Further, 45% of all known breast cancer SNPs were associated with at least one MD phenotype at p < 0.05. TWAS further identified two novel genes (SHOX2 and CRISPLD2) whose genetically predicted expression was significantly associated with MD phenotypes.

CONCLUSIONS

Our findings provided novel insight into the genetic background of MD phenotypes, and further demonstrated their shared genetic basis with breast cancer.

Alu insertion variants alter gene transcript levels

Alu are high copy number interspersed repeats that have accumulated near genes during primate and human evolution. They are a pervasive source of structural variation in modern humans. Impacts that Alu insertions may have on gene expression are not well understood, although some have been associated with expression quantitative trait loci (eQTLs). Here, we directly test regulatory effects of polymorphic Alu insertions in isolation of other variants on the same haplotype. To screen insertion variants for those with such effects, we used ectopic luciferase reporter assays and evaluated 110 Alu insertion variants, including more than 40 with a potential role in disease risk. We observed a continuum of effects with significant outliers that up- or down-regulate luciferase activity. Using a series of reporter constructs, which included genomic context surrounding the Alu, we can distinguish between instances in which the Alu disrupts another regulator and those in which the Alu introduces new regulatory sequence. We next focused on three polymorphic Alu loci associated with breast cancer that display significant effects in the reporter assay. We used CRISPR to modify the endogenous sequences, establishing cell lines varying in the Alu genotype. Our findings indicate that Alu genotype can alter expression of genes implicated in cancer risk, including PTHLH, RANBP9, and MYC These data show that commonly occurring polymorphic Alu elements can alter transcript levels and potentially contribute to disease risk.

Circular RNA hsa_circ_0043278 inhibits breast cancer progression via the miR-455-3p/EI24 signalling pathway

BACKGROUND

Breast cancer (BC) is one of the major malignancies worldwide. Circular ribonucleic acids (circRNAs) are a class of conserved ribonucleic acid (RNA) molecules that play important roles in various diseases. Recently, circRNAs have been suggested to have diagnostic value and may function as potential diagnostic biomarkers for BC. Previously, hsa_circ_0043278 was found to be downregulated in human BC. However, its role in human BC has not yet been identified.

METHODS

The levels of hsa_circ_0043278 in BC cell lines were verified by quantitative reverse transcription-polymerase chain reaction (qRT-PCR). The overexpression vector and short hairpin RNA (shRNA) of hsa_circ_0043278 were transfected into MDA-MB-231 and MCF-7 cells, respectively. The effects of hsa_circ_0043278 on tumour cell growth, migration and invasion were measured by 3-(4,5-Dimethylthiazol-2-yl)-2,5-diphenyl-tetrazolium bromide (MTT), colony formation, wound healing and Transwell assays in vitro. A xenograft experiment was conducted to validate the inhibitory effect of hsa_circ_0043278 on tumour growth. The interaction between hsa_circ_0043278 and miR-455-3p was confirmed by a dual-luciferase reporter assay. Mimics and inhibitors of miR-455-3p were designed to confirm the influence of hsa_circ_0043278 on the hsa_circ_0043278/miR-455-3p/etoposide-induced gene 24 (EI24) axis.

RESULTS

Hsa_circ_0043278 was downregulated in BC cell lines. Furthermore, overexpression of hsa_circ_0043278 notably decreased BC cell viability and inhibited BC cell migration and invasion in vitro and suppressed tumour growth in vivo. Downregulation of hsa_circ_0043278 led to the opposite results. Hsa_circ_0043278 expression was negatively correlated with that of miR-455-3p. In addition, mechanistic investigation proved that hsa_circ_0043278 directly bound to miR-455-3p and regulated EI24 and NF-κB expression in BC cells.

CONCLUSION

Hsa_circ_0043278 acts as a tumour suppressor gene in BC through the hsa_circ_0043278/miR-455-3p/EI24 axis and may be regarded as a new prognostic predictor or potential therapeutic target in BC.

Genetic variations of DNA bindings of FOXA1 and co-factors in breast cancer susceptibility

Identifying transcription factors (TFs) whose DNA bindings are altered by genetic variants that regulate susceptibility genes is imperative to understand transcriptional dysregulation in disease etiology. Here, we develop a statistical framework to analyze extensive ChIP-seq and GWAS data and identify 22 breast cancer risk-associated TFs. We find that, by analyzing genetic variations of TF-DNA bindings, the interaction of FOXA1 with co-factors such as ESR1 and E2F1, and the interaction of TFs with chromatin features (i.e., enhancers) play a key role in breast cancer susceptibility. Using genetic variants occupied by the 22 TFs, transcriptome-wide association analyses identify 52 previously unreported breast cancer susceptibility genes, including seven with evidence of essentiality from functional screens in breast relevant cell lines. We show that FOXA1 and co-factors form a core TF-transcriptional network regulating the susceptibility genes. Our findings provide additional insights into genetic variations of TF-DNA bindings (particularly for FOXA1) underlying breast cancer susceptibility.

The identification of transcription factors (TFs) whose binding sites are affected by risk genetic variants remains crucial. Here, the authors develop a statistical framework to analyse ChIP-seq and GWAS data, identify 22 breast cancer risk-associated TFs and a core TF-transcriptional network for FOXA1 and co-factors.

Functional annotation of the 2q35 breast cancer risk locus implicates a structural variant in influencing activity of a long-range enhancer element

A combination of genetic and functional approaches has identified three independent breast cancer risk loci at 2q35. A recent fine-scale mapping analysis to refine these associations resulted in 1 (signal 1), 5 (signal 2), and 42 (signal 3) credible causal variants at these loci. We used publicly available in silico DNase I and ChIP-seq data with in vitro reporter gene and CRISPR assays to annotate signals 2 and 3. We identified putative regulatory elements that enhanced cell-type-specific transcription from the IGFBP5 promoter at both signals (30- to 40-fold increased expression by the putative regulatory element at signal 2, 2- to 3-fold by the putative regulatory element at signal 3). We further identified one of the five credible causal variants at signal 2, a 1.4 kb deletion (esv3594306), as the likely causal variant; the deletion allele of this variant was associated with an average additional increase in IGFBP5 expression of 1.3-fold (MCF-7) and 2.2-fold (T-47D). We propose a model in which the deletion allele of esv3594306 juxtaposes two transcription factor binding regions (annotated by estrogen receptor alpha ChIP-seq peaks) to generate a single extended regulatory element. This regulatory element increases cell-type-specific expression of the tumor suppressor gene IGFBP5 and, thereby, reduces risk of estrogen receptor-positive breast cancer (odds ratio = 0.77, 95% CI 0.74-0.81, p = 3.1 × 10-31).

Long Noncoding RNA RUNXOR Promotes Myeloid-Derived Suppressor Cell Expansion and Functions via Enhancing Immunosuppressive Molecule Expressions during Latent HIV Infection

RUNX1 overlapping RNA (RUNXOR) is a long noncoding RNA and a key regulator of myeloid-derived suppressor cells (MDSCs) via targeting runt-related transcription factor 1 (RUNX1). We and others have previously reported MDSC expansion and inhibition of host immune responses during viral infections; however, the mechanisms regulating MDSC differentiation and suppressive functions, especially the role of RUNXOR-RUNX1 in the regulation of MDSCs in people living with HIV (PLHIV), remain unknown. In this study, we demonstrate that RUNXOR and RUNX1 expressions are upregulated in MDSCs that expand and accumulate in human PBMCs derived from PLHIV. We found that the upregulation of RUNXOR and RUNX1 is associated with the expressions of several key immunosuppressive molecules, including arginase 1, inducible NO synthase, STAT3, IL-6, and reactive oxygen species. RUNXOR and RUNX1 could positively regulate each other's expression and control the expressions of these suppressive mediators. Specifically, silencing RUNXOR or RUNX1 expression in MDSCs from PLHIV attenuated MDSC expansion and immunosuppressive mediator expressions, whereas overexpressing RUNXOR in CD33⁺ myeloid precursors from healthy subjects promoted their differentiation into MDSCs and enhanced the expression of these mediators. Moreover, loss of RUNXOR-RUNX1 function in MDSCs improved IFN-γ production from cocultured autologous CD4 T cells derived from PLHIV. These results suggest that the RUNXOR-RUNX1 axis promotes the differentiation and suppressive functions of MDSCs via regulating multiple immunosuppressive signaling molecules and may represent a potential target for immunotherapy in conjunction with antiviral therapy in PLHIV.

Identifying Novel Susceptibility Genes for Colorectal Cancer Risk From a Transcriptome-Wide Association Study of 125,478 Subjects

BACKGROUND AND AIMS

Susceptibility genes and the underlying mechanisms for the majority of risk loci identified by genome-wide association studies (GWAS) for colorectal cancer (CRC) risk remain largely unknown. We conducted a transcriptome-wide association study (TWAS) to identify putative susceptibility genes.

METHODS

Gene-expression prediction models were built using transcriptome and genetic data from the 284 normal transverse colon tissues of European descendants from the Genotype-Tissue Expression (GTEx), and model performance was evaluated using data from The Cancer Genome Atlas (n = 355). We applied the gene-expression prediction models and GWAS data to evaluate associations of genetically predicted gene-expression with CRC risk in 58,131 CRC cases and 67,347 controls of European ancestry. Dual-luciferase reporter assays and knockdown experiments in CRC cells and tumor xenografts were conducted.

RESULTS

We identified 25 genes associated with CRC risk at a Bonferroni-corrected threshold of P < 9.1 × 10-6, including genes in 4 novel loci, PYGL (14q22.1), RPL28 (19q13.42), CAPN12 (19q13.2), MYH7B (20q11.22), and MAP1L3CA (20q11.22). In 9 known GWAS-identified loci, we uncovered 9 genes that have not been reported previously, whereas 4 genes remained statistically significant after adjusting for the lead risk variant of the locus. Through colocalization analysis in GWAS loci, we additionally identified 12 putative susceptibility genes that were supported by TWAS analysis at P < .01. We showed that risk allele of the lead risk variant rs1741640 affected the promoter activity of CABLES2. Knockdown experiments confirmed that CABLES2 plays a vital role in colorectal carcinogenesis.

CONCLUSIONS

Our study reveals new putative susceptibility genes and provides new insight into the biological mechanisms underlying CRC development.

A case-only study to identify genetic modifiers of breast cancer risk for BRCA1/BRCA2 mutation carriers

Breast cancer (BC) risk for BRCA1 and BRCA2 mutation carriers varies by genetic and familial factors. About 50 common variants have been shown to modify BC risk for mutation carriers. All but three, were identified in general population studies. Other mutation carrier-specific susceptibility variants may exist but studies of mutation carriers have so far been underpowered. We conduct a novel case-only genome-wide association study comparing genotype frequencies between 60,212 general population BC cases and 13,007 cases with BRCA1 or BRCA2 mutations. We identify robust novel associations for 2 variants with BC for BRCA1 and 3 for BRCA2 mutation carriers, P < 10-8, at 5 loci, which are not associated with risk in the general population. They include rs60882887 at 11p11.2 where MADD, SP11 and EIF1, genes previously implicated in BC biology, are predicted as potential targets. These findings will contribute towards customising BC polygenic risk scores for BRCA1 and BRCA2 mutation carriers.

Overexpression of miR-331 Indicates Poor Prognosis and Promotes Progression of Breast Cancer

BACKGROUND

With the increasing number of cases of breast cancer every year, the exploration of novel biomarkers has drawn attention. miR-331 has been demonstrated to play a role in various cancers, but its role in breast cancer is still unknown.

METHODS

In this study, we included 121 patients with breast cancer treated at Affiliated Hospital of Weifang Medical University. Breast cancer tissues and adjacent normal tissues were collected during the surgery. The expression of miR-331 in breast cancer tissues and cell lines was detected by qRT-PCR assay. Then, with the help of Kaplan-Meier survival and Cox regression analyses, the role of miR-331 in the prognosis of breast cancer was analyzed. Finally, the effect of miR-331 on cell proliferation, migration, and invasion was investigated with CCK-8 assay and transwell assay.

RESULTS

miR-331 was significantly upregulated in breast cancer tissues compared with normal tissues. The overexpression of miR-331 was associated with lymph node metastasis, TNM stage, and poor prognosis. From the results of Cox regression analyses, it was found that miR-331 served as an independent indicator in the prognosis of breast cancer. In addition, miR-331 was also found to be upregulated in breast cancer cells, which promoted cell proliferation, migration, and invasion of breast cancer.

CONCLUSION

As shown from our data, miR-331 may be a potential prognostic biomarker in breast cancer. Moreover, the development and progression of breast cancer may involve miR-331. These findings suggest a novel therapeutic target for the treatment of breast cancer.

Characterization of hepatitis B virus infection and viral DNA integration in non-Hodgkin lymphoma

Hepatitis B virus (HBV) infection has been reported to be associated with non-Hodgkin lymphoma (NHL). However, the evidence is limited to the seroepidemiological study. There is a lack of evidence showing the HBV infection and integration in NHL cells. Here, we reported that in the Shanghai area, the positive rates of serum HBsAg (OR: 3.11; 95% CI: 2.20-4.41) and HBeAg (OR: 3.99; 95% CI: 1.73-9.91) were significantly higher in patients with NHL. HBsAg, HBcAg and HBV DNA were detected in 34.4%, 45.2% and 47.0% of the NHL tissues, respectively. Furthermore, by using a high-throughput viral integration detection approach (HIVID), integrated HBV DNA was identified from 50% (6/12) HBV-related NHL tissues. There were a total of 313 HBV integration sites isolated from the NHL tissues, among which four protein-coding genes (FAT2, SETX, ITGA10 and CD63) were interrupted by HBV DNA in their exons. Seven HBV preferential target genes (ANKS1B, HDAC4, EGFLAM, MAN1C1, XKR6, ZBTB38 and CCDC91) showed significantly altered expression levels in NHL, suggesting a potential role of these genes in NHL development. Taken together, HBV integration is a common phenomenon in NHL. This finding opens up a new direction of research into the mechanistic link between HBV infection and NHL.

Fine-mapping of a novel premenopausal breast cancer susceptibility locus at Chr4q31.22 in Caucasian women and validation in African and Chinese women

We previously identified a novel breast cancer susceptibility variant on chromosome 4q31.22 locus (rs1429142) conferring risk among women of European ancestry. Here, we report replication of findings, validation of the variant in diverse populations and fine-mapping of the associated locus in Caucasian population. The SNP rs1429142 (C/T, minor allele frequency 18%) showed association for the overall breast cancer risk in Stages 1-4 (n = 4,331 cases/4271 controls; p = 4.35 × 10-8 ; odds ratio, ORC-allele ,1.25), and an elevated risk among premenopausal women (n = 1,503 cases/4271 controls; p = 5.81 × 10-10 ; ORC-allele 1.40) in European populations. SNP rs1429142 was associated with premenopausal breast cancer risk in women of African (T/C; p-value 1.45 × 10-02 ; ORC-allele 1.2) but not from Chinese ancestry. Fine-mapping of the locus revealed several potential causal variants which are present within a single association signal, revealed from the conditional regression analysis. Functional annotation of the potential causal variants revealed three putative SNPs rs1366691, rs1429139 and rs7667633 with active enhancer functions inferred based on histone marks, DNase hypersensitive sites in breast cell line data. These putative variants were bound by transcription factors (C-FOS, STAT1/3 and POL2/3) with known roles in inflammatory pathways. Furthermore, Hi-C data revealed several short-range interactions in the fine-mapped locus harboring the putative variants. The fine mapped locus was predicted to be within a single topologically associated domain, potentially facilitating enhancer-promoter interactions possibly leading to the regulation of nearby genes.

Fine-mapping of 150 breast cancer risk regions identifies 191 likely target genes

Genome-wide association studies have identified breast cancer risk variants in over 150 genomic regions, but the mechanisms underlying risk remain largely unknown. These regions were explored by combining association analysis with in silico genomic feature annotations. We defined 205 independent risk-associated signals with the set of credible causal variants in each one. In parallel, we used a Bayesian approach (PAINTOR) that combines genetic association, linkage disequilibrium and enriched genomic features to determine variants with high posterior probabilities of being causal. Potentially causal variants were significantly over-represented in active gene regulatory regions and transcription factor binding sites. We applied our INQUSIT pipeline for prioritizing genes as targets of those potentially causal variants, using gene expression (expression quantitative trait loci), chromatin interaction and functional annotations. Known cancer drivers, transcription factors and genes in the developmental, apoptosis, immune system and DNA integrity checkpoint gene ontology pathways were over-represented among the highest-confidence target genes.

Identification of BRCA1/2 p.Ser1613Gly, p.Pro871Leu, p.Lys1183Arg, p.Glu1038Gly, p.Ser1140Gly, p.Ala2466Val, p.His2440Arg variants in women under 45 years old with breast nodules suspected of having breast cancer in Burkina Faso

Breast cancer is the top cause of cancer mortality among women in the world and the second in Africa. The aims of this study were to: i) identify women with breast nodules suspected of having breast cancer ii) sequence the BRCA1 and BRCA2 genes and iii) screen mutations. From 2015 to 2016, 112 women aged from 35 to 44 years, who had come for consultation in the gynecology/obstetrics and the oncology department of the University Hospital Yalgado Ouedraogo, voluntarily agreed to participate to this study. Whole blood was collected from those with mammary nodules. The genomic DNA was extracted using Qiagen kit. FAST KAPA was used for genomic DNA amplification and the purified PCR products were analyzed by direct sequencing using Big Dye v1.1 and ABI 3730 automated sequencer. Nucleotides substitutions were determined. We identified BRCA1 SNPs rs1799966, rs799917, rs16942, rs16941, rs2227945, and BRCA2 SNPs rs169547, rs4986860. These identified variants are found mostly in cases of benign tumors of breast or ovarian cancer with familial history of breast cancer. This study in Burkina-Faso, is the basis for improved and more specific genetic testing, and suggests that additional genes contributing to an increased risk of breast cancer should be analyzed.

The SuprMam1 breast cancer susceptibility locus disrupts the vitamin D/ calcium/ parathyroid hormone pathway and alters bone structure in congenic mice

Breast cancer is a complex disease, and approximately 30% of cases are considered to be hereditary or familial, with a large fraction of this being polygenic. However, it is difficult to demonstrate the functional importance of genes of small effect in population studies, and these genes are not always easily targeted for prevention. The SuprMam (suppressor of mammary tumour) breast cancer susceptibility alleles were previously identified as contributors to spontaneous mammary tumour development in Trp53^+/- mice. In this study, we have generated and characterised congenic mice that contain the BALB/c SuprMam1 (susceptibility) locus on a C57BL/6 (resistant) background and discovered a subtle impairment in the vitamin D/ calcium/ parathyroid hormone (PTH) pathway. This was evident as altered gene expression in the mammary glands of key players in this pathway. Further functional analysis of the mice revealed elevated PTH levels, reduced Cyp27b1 expression in kidneys, and reduced trabecular bone volume/ tissue volume percentage. Plasma 25(OH)D and serum calcium were unchanged. This impairment was a result of genetic differences and occurred only in females, but the elevated PTH levels could be overcome with either calcium or vitamin D dietary supplementation. Either low levels of active vitamin D (1,25(OH)₂D) or chronically elevated PTH levels may contribute to increased breast cancer susceptibility. These indicators are not easily measured in human population studies, but either mechanism may be preventable with dietary calcium or vitamin D supplements. Therefore, SuprMam congenic mice could serve as a valuable model for studying the role of gene-hormone-environment interactions of the vitamin D/ calcium/ PTH pathway in cancer and other diseases and for testing preventive interventions.

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.

Early Diagnosis of Sepsis: Is an Integrated Omics Approach the Way Forward?

Sepsis remains one of the leading causes of death in the USA and it is expected to get worse as the population ages. Moreover, the standard of care, which recommends aggressive treatment with appropriate antibiotics, has led to an increase in multiple drug-resistant organisms. There is a dire need for the development of new antibiotics, improved antibiotic stewardship, and therapies that treat the host response. Development of new sepsis therapeutics has been a disappointment as no drugs are currently approved to treat the various complications from sepsis. Much of the failure has been blamed on animal models that do not accurately reflect the course of the disease. However, recent improvements in metabolomic, transcriptomic, genomic, and proteomic platforms have allowed for a broad-spectrum look at molecular changes in the host response using clinical samples. Integration of these multi-omic datasets allows researchers to perform systems biology approaches to identify novel pathophysiology of the disease. In this review, we highlight what is currently known about sepsis and how integrative omics has identified new diagnostic and predictive models of sepsis as well as novel mechanisms. These changes may improve patient care as well as guide future preclinical analysis of sepsis.

Long non-coding RNA RUNXOR accelerates MDSC-mediated immunosuppression in lung cancer

Background

RUNX1 overlapping RNA (RUNXOR) is a long non-coding RNA that has been indicated as a key regulator in the development of myeloid cells by targeting runt-related transcription factor 1 (RUNX1). Myeloid-derived suppressor cells (MDSCs) are a heterogeneous population of cells consisting of immature granulocytes and monocytes with immunosuppression. However, the impact of lncRNA RUNXOR on the development of MDSCs remains unknown.

Methods

Both the expressions of RUNXOR and RUNX1 in the peripheral blood were measured by qRT-PCR. Human MDSCs used in this study were isolated from tumor tissue of patients with lung cancer by FCM or induced from PBMCs of healthy donors with IL-1β + GM-CSF. Specific siRNA was used to knockdown the expression of RUNXOR in MDSCs.

Results

In this study, we found that the lncRNA RUNXOR was upregulated in the peripheral blood of lung cancer patients. In addition, as a target gene of RUNXOR, the expression of RUNX1 was downregulated in lung cancer patients. Finally, the expression of RUNXOR was higher in MDSCs isolated from the tumor tissues of lung cancer patients compared with cells from adjacent tissue. In addition, RUNXOR knockdown decreased Arg1 expression in MDSCs.

Conclusions

Based on our findings, it is illustrated that RUNXOR is significantly associated with the immunosuppression induced by MDSCs in lung cancer patients and may be a target of anti-tumor therapy.

A Comprehensive cis-eQTL Analysis Revealed Target Genes in Breast Cancer Susceptibility Loci Identified in Genome-wide Association Studies

Genome-wide association studies (GWASs) have identified more than 150 common genetic loci for breast cancer risk. However, the target genes and underlying mechanisms remain largely unknown. We conducted a cis-expression quantitative trait loci (cis-eQTL) analysis using normal or tumor breast transcriptome data from the Molecular Taxonomy of Breast Cancer International Consortium (METABRIC), The Cancer Genome Atlas (TCGA), and the Genotype-Tissue Expression (GTEx) project. We identified a total of 101 genes for 51 lead variants after combing the results of a meta-analysis of METABRIC and TCGA, and the results from GTEx at a Benjamini-Hochberg (BH)-adjusted p < 0.05. Using luciferase reporter assays in both estrogen-receptor positive (ER⁺) and negative (ER^-) cell lines, we showed that alternative alleles of potential functional single-nucleotide polymorphisms (SNPs), rs11552449 (DCLRE1B), rs7257932 (SSBP4), rs3747479 (MRPS30), rs2236007 (PAX9), and rs73134739 (ATG10), could significantly change promoter activities of their target genes compared to reference alleles. Furthermore, we performed in vitro assays in breast cancer cell lines, and our results indicated that DCLRE1B, MRPS30, and ATG10 played a vital role in breast tumorigenesis via certain disruption of cell behaviors. Our findings revealed potential target genes for associations of genetic susceptibility risk loci and provided underlying mechanisms for a better understanding of the pathogenesis of breast cancer.

Identification of circular RNAs as a promising new class of diagnostic biomarkers for human breast cancer

Endogenous noncoding circular RNAs (circRNAs) have gained attention for their involvement in carcinogenesis, but their expression pattern in breast cancer has remained largely unknown. In this two-stage study, we first used an Arraystar Human circRNA Array to construct a genome-wide circRNA profile. We then selected candidate circRNAs for validation using a quantitative real-time polymerase chain reaction system. CircRNA/miRNA interactions were predicted and sequence analyses were performed. Among 1155 differentially expressed circRNAs, 715 were upregulated and 440 were downregulated in breast cancer tissues. The validation study demonstrated that hsa_circ_103110, hsa_circ_104689 and hsa_circ_104821 levels were elevated in breast cancer tissues, whereas hsa_circ_006054, hsa_circ_100219 and hsa_circ_406697 were downregulated. These circRNAs targeted complementary miRNA response elements. The area under the receiver operating characteristic curve for distinguishing breast cancer was 0.82 (95% CI: 0.73-0.90) when hsa_circ_006054, hsa_circ_100219 and hsa_circ_406697 were used in combination. This study provides evidence that circRNAs are differentially expressed in breast cancer and are important in carcinogenesis because they participate in cancer-related pathways and sequester miRNAs.

Capture Hi-C identifies putative target genes at 33 breast cancer risk loci

Genome-wide association studies (GWAS) have identified approximately 100 breast cancer risk loci. Translating these findings into a greater understanding of the mechanisms that influence disease risk requires identification of the genes or non-coding RNAs that mediate these associations. Here, we use Capture Hi-C (CHi-C) to annotate 63 loci; we identify 110 putative target genes at 33 loci. To assess the support for these target genes in other data sources we test for associations between levels of expression and SNP genotype (eQTLs), disease-specific survival (DSS), and compare them with somatically mutated cancer genes. 22 putative target genes are eQTLs, 32 are associated with DSS and 14 are somatically mutated in breast, or other, cancers. Identifying the target genes at GWAS risk loci will lead to a greater understanding of the mechanisms that influence breast cancer risk and prognosis.

Characterizing Genetic Susceptibility to Breast Cancer in Women of African Ancestry

Background: Genome-wide association studies have identified approximately 100 common genetic variants associated with breast cancer risk, the majority of which were discovered in women of European ancestry. Because of different patterns of linkage disequilibrium, many of these genetic markers may not represent signals in populations of African ancestry.Methods: We tested 74 breast cancer risk variants and conducted fine-mapping of these susceptibility regions in 6,522 breast cancer cases and 7,643 controls of African ancestry from three genetic consortia (AABC, AMBER, and ROOT).Results: Fifty-four of the 74 variants (73%) were found to have ORs that were directionally consistent with those previously reported, of which 12 were nominally statistically significant (P < 0.05). Through fine-mapping, in six regions (3p24, 12p11, 14q13, 16q12/FTO, 16q23, 19p13), we observed seven markers that better represent the underlying risk variant for overall breast cancer or breast cancer subtypes, whereas in another two regions (11q13, 16q12/TOX3), we identified suggestive evidence of signals that are independent of the reported index variant. Overlapping chromatin features and regulatory elements suggest that many of the risk alleles lie in regions with biological functionality.Conclusions: Through fine-mapping of known susceptibility regions, we have revealed alleles that better characterize breast cancer risk in women of African ancestry.Impact: The risk alleles identified represent genetic markers for modeling and stratifying breast cancer risk in women of African ancestry. Cancer Epidemiol Biomarkers Prev; 26(7); 1016-26. ©2017 AACR.

Expression Quantitative Trait loci (QTL) in tumor adjacent normal breast tissue and breast tumor tissue

We investigate 71 single nucleotide polymorphisms (SNPs) identified in meta-analytic studies of genome-wide association studies (GWAS) of breast cancer, the majority of which are located in intergenic or intronic regions. To explore regulatory impacts of these variants we conducted expression quantitative loci (eQTL) analyses on tissue samples from 376 invasive postmenopausal breast cancer cases in the Nurses' Health Study (NHS) diagnosed from 1990-2004. Expression analysis was conducted on all formalin-fixed paraffin-embedded (FFPE) tissue samples (and on 264 adjacent normal samples) using the Affymetrix Human Transcriptome Array. Significance and ranking of associations between tumor receptor status and expression variation was preserved between NHS FFPE and TCGA fresh-frozen sample sets (Spearman r = 0.85, p<10^-10 for 17 of the 21 Oncotype DX recurrence signature genes). At an FDR threshold of 10%, we identified 27 trans-eQTLs associated with expression variation in 217 distinct genes. SNP-gene associations can be explored using an open-source interactive browser distributed in a Bioconductor package. Using a new a procedure for testing hypotheses relating SNP content to expression patterns in gene sets, defined as molecular function pathways, we find that loci on 6q14 and 6q25 affect various gene sets and molecular pathways (FDR < 10%). Although the ultimate biological interpretation of the GWAS-identified variants remains to be uncovered, this study validates the utility of expression analysis of this FFPE expression set for more detailed integrative analyses.