Cis-eQTL analysis and functional validation of candidate susceptibility genes for high-grade serous ovarian cancer

Integrative multi-omics analyses to identify the genetic and functional mechanisms underlying ovarian cancer risk regions

To identify credible causal risk variants (CCVs) associated with different histotypes of epithelial ovarian cancer (EOC), we performed genome-wide association analysis for 470,825 genotyped and 10,163,797 imputed SNPs in 25,981 EOC cases and 105,724 controls of European origin. We identified five histotype-specific EOC risk regions (p value <5 × 10-8) and confirmed previously reported associations for 27 risk regions. Conditional analyses identified an additional 11 signals independent of the primary signal at six risk regions (p value <10-5). Fine mapping identified 4,008 CCVs in these regions, of which 1,452 CCVs were located in ovarian cancer-related chromatin marks with significant enrichment in active enhancers, active promoters, and active regions for CCVs from each EOC histotype. Transcriptome-wide association and colocalization analyses across histotypes using tissue-specific and cross-tissue datasets identified 86 candidate susceptibility genes in known EOC risk regions and 32 genes in 23 additional genomic regions that may represent novel EOC risk loci (false discovery rate <0.05). Finally, by integrating genome-wide HiChIP interactome analysis with transcriptome-wide association study (TWAS), variant effect predictor, transcription factor ChIP-seq, and motifbreakR data, we identified candidate gene-CCV interactions at each locus. This included risk loci where TWAS identified one or more candidate susceptibility genes (e.g., HOXD-AS2, HOXD8, and HOXD3 at 2q31) and other loci where no candidate gene was identified (e.g., MYC and PVT1 at 8q24) by TWAS. In summary, this study describes a functional framework and provides a greater understanding of the biological significance of risk alleles and candidate gene targets at EOC susceptibility loci identified by a genome-wide association study.

Systems genetics analysis reveals the common genetic basis for pain sensitivity and cognitive function

BACKGROUND

There is growing evidence of a strong correlation between pain sensitivity and cognitive function under both physiological and pathological conditions. However, the detailed mechanisms remain largely unknown. In the current study, we sought to explore candidate genes and common molecular mechanisms underlying pain sensitivity and cognitive function with a transcriptome-wide association study using recombinant inbred mice from the BXD family.

METHODS

The pain sensitivity determined by Hargreaves' paw withdrawal test and cognition-related phenotypes were systematically analyzed in 60 strains of BXD mice and correlated with hippocampus transcriptomes, followed by quantitative trait locus (QTL) mapping and systems genetics analysis.

RESULTS

The pain sensitivity showed significant variability across the BXD strains and co-varies with cognitive traits. Pain sensitivity correlated hippocampual genes showed a significant involvement in cognition-related pathways, including glutamatergic synapse, and PI3K-Akt signaling pathway. Moreover, QTL mapping identified a genomic region on chromosome 4, potentially regulating the variation of pain sensitivity. Integrative analysis of expression QTL mapping, correlation analysis, and Bayesian network modeling identified Ring finger protein 20 (Rnf20) as the best candidate. Further pathway analysis indicated that Rnf20 may regulate the expression of pain sensitivity and cognitive function through the PI3K-Akt signaling pathway, particularly through interactions with genes Ppp2r2b, Ppp2r5c, Col9a3, Met, Rps6, Tnc, and Kras.

CONCLUSIONS

Our study demonstrated that pain sensitivity is associated with genetic background and Rnf20-mediated PI3K-Akt signaling may involve in the regulation of pain sensitivity and cognitive functions.

Immunosuppression causes dynamic changes in expression QTLs in psoriatic skin

Psoriasis is a chronic, systemic inflammatory condition primarily affecting skin. While the role of the immune compartment (e.g., T cells) is well established, the changes in the skin compartment are more poorly understood. Using longitudinal skin biopsies (n = 375) from the "Psoriasis Treatment with Abatacept and Ustekinumab: A Study of Efficacy"(PAUSE) clinical trial (n = 101), we report 953 expression quantitative trait loci (eQTLs). Of those, 116 eQTLs have effect sizes that were modulated by local skin inflammation (eQTL interactions). By examining these eQTL genes (eGenes), we find that most are expressed in the skin tissue compartment, and a subset overlap with the NRF2 pathway. Indeed, the strongest eQTL interaction signal - rs1491377616-LCE3C - links a psoriasis risk locus with a gene specifically expressed in the epidermis. This eQTL study highlights the potential to use biospecimens from clinical trials to discover in vivo eQTL interactions with therapeutically relevant environmental variables.

Mendelian randomization uncovers a protective effect of interleukin-1 receptor antagonist on kidney function

Interleukins (ILs), key cytokine family of inflammatory response, are closely associated with kidney function. However, the causal effect of various ILs on kidney function needs further investigation. Here we show two-sample summary-level Mendelian randomization (MR) analysis that examined the causality between serum IL levels and kidney function. Genetic variants with strong association with serum IL levels were obtained from a previous genome-wide association study meta-analysis. Summary-level data for estimated glomerular filtration rate (eGFR) were obtained from CKDGen database. As a main MR analysis, multiplicative random-effects inverse-variance weighted method was performed. Pleiotropy-robust MR analysis, including MR-Egger with bootstrapped error and weighted median methods, were also implemented. We tested the causal estimates from nine ILs on eGFR traits. Among the results, higher genetically predicted serum IL-1 receptor antagonist level was significantly associated with higher eGFR values in the meta-analysis of CKDGen and the UK Biobank data. In addition, the result was consistent towards eGFR decline phenotype of the outcome database. Otherwise, nonsignificant association was identified between other genetically predicted ILs and eGFR outcome. These findings support the clinical importance of IL-1 receptor antagonist-associated pathway in relation to kidney function in the general individuals, particularly highlighting the importance of IL-1 receptor antagonist.

InferLoop: leveraging single-cell chromatin accessibility for the signal of chromatin loop

Deciphering cell-type-specific 3D structures of chromatin is challenging. Here, we present InferLoop, a novel method for inferring the strength of chromatin interaction using single-cell chromatin accessibility data. The workflow of InferLoop is, first, to conduct signal enhancement by grouping nearby cells into bins, and then, for each bin, leverage accessibility signals for loop signals using a newly constructed metric that is similar to the perturbation of the Pearson correlation coefficient. In this study, we have described three application scenarios of InferLoop, including the inference of cell-type-specific loop signals, the prediction of gene expression levels and the interpretation of intergenic loci. The effectiveness and superiority of InferLoop over other methods in those three scenarios are rigorously validated by using the single-cell 3D genome structure data of human brain cortex and human blood, the single-cell multi-omics data of human blood and mouse brain cortex, and the intergenic loci in the GWAS Catalog database as well as the GTEx database, respectively. In addition, InferLoop can be applied to predict loop signals of individual spots using the spatial chromatin accessibility data of mouse embryo. InferLoop is available at https://github.com/jumphone/inferloop.

Molecular Management of High-Grade Serous Ovarian Carcinoma

High-grade serous ovarian carcinoma (HGSOC) represents the most common form of epithelial ovarian carcinoma. The absence of specific symptoms leads to late-stage diagnosis, making HGSOC one of the gynecological cancers with the worst prognosis. The cellular origin of HGSOC and the role of reproductive hormones, genetic traits (such as alterations in P53 and DNA-repair mechanisms), chromosomal instability, or dysregulation of crucial signaling pathways have been considered when evaluating prognosis and response to therapy in HGSOC patients. However, the detection of HGSOC is still based on traditional methods such as carbohydrate antigen 125 (CA125) detection and ultrasound, and the combined use of these methods has yet to support significant reductions in overall mortality rates. The current paradigm for HGSOC management has moved towards early diagnosis via the non-invasive detection of molecular markers through liquid biopsies. This review presents an integrated view of the relevant cellular and molecular aspects involved in the etiopathogenesis of HGSOC and brings together studies that consider new horizons for the possible early detection of this gynecological cancer.

Long Noncoding RNA LINC00909 Induces Epithelial-Mesenchymal Transition and Contributes to Osteosarcoma Tumorigenesis and Metastasis

Background

Osteosarcoma (OS) is a malignant tumor that is highly metastatic with a high mortality rate. Although mounting evidence suggests that LINC00909 is strongly associated with the malignant progression of various tumors, the exact role of LINC00909 in OS remains unknown. Therefore, the current study was designed to investigate the relationship between LINC00909 and the malignant progression of OS.

Methods

LINC00909 expression was measured in OS cell lines and clinical specimens using RT-qPCR assays. The effects of LINC00909 on OS proliferation, invasion, and migration were calculated both in vitro and in vivo. Apart from that, bioinformatics analyses, FISH, RIP, and luciferase reporter assays were carried out to investigate the downstream target of LINC00909. Rescue experiments were also conducted to investigate the potential mechanisms of action of competitive endogenous RNAs (ceRNAs).

Results

In this study, we found that LINC00909 was highly expressed in OS cell lines and clinical specimens. In vivo and in vitro experiments demonstrated that LINC00909 induces epithelial-to-mesenchymal transition (EMT) and contributes to OS tumorigenesis and metastasis. FISH, RIP, and luciferase assays indicated that miR-875-5p is a direct target of LINC00909. Moreover, HOXD9 was validated as the downstream target of miR-875-5p in a luciferase reporter assay and western blotting experiments. Rescue experiments revealed that HOXD9 reversed the effect of LV-sh-LINC00909 on OS cells by positively regulating the PI3K/AKT/mTOR signaling pathway.

Conclusion

Collectively, LINC00909 induces EMT and contributes to OS tumorigenesis and metastasis through the PI3K/AKT/mTOR pathway by binding to miR-875-5p to upregulate HOXD9 expression. Targeting the LINC00909/miR-875-5p/HOXD9 axis may have potential in treating OS.

chromMAGMA: regulatory element-centric interrogation of risk variants

Candidate causal risk variants from genome-wide association studies reside almost exclusively in noncoding regions of the genome and innovative approaches are necessary to understand their biological function. Multi-marker analysis of genomic annotation (MAGMA) is a widely used program that nominates candidate risk genes by mapping single-nucleotide polymorphism summary statistics from genome-wide association studies to gene bodies. We augmented MAGMA to create chromatin-MAGMA (chromMAGMA), a method to nominate candidate risk genes based on the presence of risk variants within noncoding regulatory elements (REs). We applied chromMAGMA to a genetic susceptibility dataset for epithelial ovarian cancer (EOC), a rare gynecologic malignancy characterized by high mortality. This identified 155 unique candidate EOC risk genes across five EOC histotypes; 83% (105/127) of high-grade serous ovarian cancer risk genes had not previously been implicated in this EOC histotype. Risk genes nominated by chromMAGMA converged on mRNA splicing and transcriptional dysregulation pathways. chromMAGMA is a pipeline that nominates candidate risk genes through a gene regulation-focused approach and helps interpret the biological mechanism of noncoding risk variants for complex diseases.

Identification of UAP1L1 as a critical factor for prostate cancer and underlying molecular mechanism in tumorigenicity

Background

Prostate cancer is the second most common cancer in men, and some new target genes are needed to predict the risk of prostate cancer progression and the treatment.

Methods

In this study, the effects of UAP1L1 (UAP1-like-1) on prostate cancer were investigated by detecting the proliferation, migration, invasion and apoptosis of prostate cancer cells in vitro using MTT, wound healing, Transwell and flow cytometry assay, and the tumor growth in vivo. The downstream genes and pathways of UAP1L1 were explored using Ingenuity Pathway Analysis (IPA), and screened by qRT-PCR and western blot. The effects of CDCA8 on prostate cancer cells were also verified in vitro, which was through detecting the change of proliferation, migration, invasion and apoptosis of prostate cancer cells after CDCA8 knockdown.

Results

The results indicated that UAP1L1 promoted the proliferation, migration and invasion of prostate cancer cells, which was inhibited by downregulating CDCA8. Furthermore, the promotion of CDCA8 knockdown on cell apoptosis was reduced when UAP1L1 was simultaneously overexpressed.

Conclusions

In conclusion, the results in this study revealed that UAP1L1 promoted the progression of prostate cancer through the downstream gene CDCA8.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12967-022-03291-0.

Genome and transcriptome profiling of spontaneous preterm birth phenotypes

Preterm birth (PTB) occurs before 37 weeks of gestation. Risk factors include genetics and infection/inflammation. Different mechanisms have been reported for spontaneous preterm birth (SPTB) and preterm birth following preterm premature rupture of membranes (PPROM). This study aimed to identify early pregnancy biomarkers of SPTB and PPROM from the maternal genome and transcriptome. Pregnant women were recruited at the Liverpool Women’s Hospital. Pregnancy outcomes were categorised as SPTB, PPROM (≤ 34 weeks gestation, n = 53), high-risk term (HTERM, ≥ 37 weeks, n = 126) or low-risk (no history of SPTB/PPROM) term (LTERM, ≥ 39 weeks, n = 188). Blood samples were collected at 16 and 20 weeks gestation from which, genome (UK Biobank Axiom array) and transcriptome (Clariom D Human assay) data were acquired. PLINK and R were used to perform genetic association and differential expression analyses and expression quantitative trait loci (eQTL) mapping. Several significant molecular signatures were identified across the analyses in preterm cases. Genome-wide significant SNP rs14675645 (ASTN1) was associated with SPTB whereas microRNA-142 transcript and PPARG1-FOXP3 gene set were associated with PPROM at week 20 of gestation and is related to inflammation and immune response. This study has determined genomic and transcriptomic candidate biomarkers of SPTB and PPROM that require validation in diverse populations.

Functional analysis of the 1p34.3 risk locus implicates GNL2 in high-grade serous ovarian cancer

The high-grade serous ovarian cancer (HGSOC) risk locus at chromosome 1p34.3 resides within a frequently amplified genomic region signifying the presence of an oncogene. Here, we integrate in silico variant-to-function analysis with functional studies to characterize the oncogenic potential of candidate genes in the 1p34.3 locus. Fine mapping of genome-wide association statistics identified candidate causal SNPs local to H3K27ac-demarcated enhancer regions that exhibit allele-specific binding for CTCF in HGSOC and normal fallopian tube secretory epithelium cells (FTSECs). SNP risk associations colocalized with eQTL for six genes (DNALI1, GNL2, RSPO1, SNIP1, MEAF6, and LINC01137) that are more highly expressed in carriers of the risk allele, and three (DNALI1, GNL2, and RSPO1) were upregulated in HGSOC compared to normal ovarian surface epithelium cells and/or FTSECs. Increased expression of GNL2 and MEAF6 was associated with shorter survival in HGSOC with 1p34.3 amplifications. Despite its activation of β-catenin signaling, RSPO1 overexpression exerted no effects on proliferation or colony formation in our study of ovarian cancer and FTSECs. Instead, GNL2, MEAF6, and SNIP1 silencing impaired in vitro ovarian cancer cell growth. Additionally, GNL2 silencing diminished xenograft tumor formation, whereas overexpression stimulated proliferation and colony formation in FTSECs. GNL2 influences 60S ribosomal subunit maturation and global protein synthesis in ovarian cancer and FTSECs, providing a potential mechanism of how GNL2 upregulation might promote ovarian cancer development and mediate genetic susceptibility of HGSOC.

Scratching Below the Ovarian Cancer GWAS Surface

Despite recent notable treatment advancements, ovarian cancer survival rates remain poor, with about half of women surviving five years after diagnosis. Uncovering novel prognostic factors is critical to better understand and reduce mortality from this deadly disease. While genome-wide association studies have identified numerous loci associated with risk of epithelial ovarian cancer, the investigation of genetic factors associated with outcomes among women with ovarian cancer has been limited due to several challenges summarized in the present commentary. Using data from the Ovarian Cancer Association Consortium, Quinn and colleagues conducted a genome-wide association study of patients with ovarian cancer receiving debulking surgery and standard chemotherapy as first-line treatment, revealing a locus at 12q24.33 associated with progression-free survival. Experimental evidence suggests that ULK1, a gene coding for a serine/threonine kinase implicated in autophagy, is the target of the association. We discuss the novelty of these findings, unanswered questions, and next steps for the road ahead in translating the work of Quinn and colleagues into clinical practice.See related article by Quinn et al., p. 1669.

Supervised Dimension Reduction for Large-Scale "Omics" Data With Censored Survival Outcomes Under Possible Non-Proportional Hazards

The past two decades have witnessed significant advances in high-throughput "omics" technologies such as genomics, proteomics, metabolomics, transcriptomics and radiomics. These technologies have enabled simultaneous measurement of the expression levels of tens of thousands of features from individual patient samples and have generated enormous amounts of data that require analysis and interpretation. One specific area of interest has been in studying the relationship between these features and patient outcomes, such as overall and recurrence-free survival, with the goal of developing a predictive "omics" profile. Large-scale studies often suffer from the presence of a large fraction of censored observations and potential time-varying effects of features, and methods for handling them have been lacking. In this paper, we propose supervised methods for feature selection and survival prediction that simultaneously deal with both issues. Our approach utilizes continuum power regression (CPR) - a framework that includes a variety of regression methods - in conjunction with the parametric or semi-parametric accelerated failure time (AFT) model. Both CPR and AFT fall within the linear models framework and, unlike black-box models, the proposed prognostic index has a simple yet useful interpretation. We demonstrate the utility of our methods using simulated and publicly available cancer genomics data.

Effects of long-term norepinephrine treatment on normal immortalized ovarian and fallopian tube cells

Sustained adrenergic stimulation by norepinephrine (NE) contributes to ovarian carcinoma metastasis and impairment of chemotherapy response. Although the effect of sustained NE stimulation in cancer progression is well established, less is known about its role in cancer initiation. To determine the extent to which stress hormones influence ovarian cancer initiation, we conducted a long-term (> 3 months; > 40 population doublings) experiment in which normal immortalized fallopian tube secretory (iFTSEC283) and ovarian surface epithelial (iOSE11) cell lines and their isogenic pairs containing a p53 mutation (iFTSEC283p53R175H; iOSE11p53R175H), were continuously exposed to NE (100 nM, 1 μM, 10 μM). Fallopian tube cells displayed a p53-independent increase in proliferation and colony-forming ability in response to NE, while ovarian surface epithelial cells displayed a p53-independent decrease in both assays. Fallopian tube cells with mutant p53 showed a mild loss of chromosomes and TP53 status was also a defining factor in transcriptional response of fallopian tube cells to long-term NE treatment.

Transcription factors in epithelial ovarian cancer: histotype-specific drivers and novel therapeutic targets

Transcription factors (TFs) are major contributors to cancer risk and somatic development. In preclinical and clinical studies, direct or indirect inhibition of TF-mediated oncogenic gene expression profiles have proven to be effective in many tumor types, highlighting this group of proteins as valuable therapeutic targets. In spite of this, our understanding of TFs in epithelial ovarian cancer (EOC) is relatively limited. EOC is a heterogeneous disease composed of five major histologic subtypes; high-grade serous, low-grade serous, endometrioid, clear cell and mucinous. Each histology is associated with unique clinical etiologies, sensitivity to therapies, and molecular signatures - including diverse transcriptional regulatory programs. While some TFs are shared across EOC subtypes, a set of TFs are expressed in a histotype-specific manner and likely explain part of the histologic diversity of EOC subtypes. Targeting TFs present with unique opportunities for development of novel precision medicine strategies for ovarian cancer. This article reviews the critical TFs in EOC subtypes and highlights the potential of exploiting TFs as biomarkers and therapeutic targets.

Cis- and Trans-Acting Expression Quantitative Trait Loci of Long Non-Coding RNA in 2,549 Cancers With Potential Clinical and Therapeutic Implications

Many cancer risk loci act as expression quantitative trait loci (eQTLs) of transcripts including non-coding RNA. Long non-coding RNAs (lncRNAs) are implicated in various human cancers. However, the pathological and clinical impacts of the genetic determinants of lncRNAs in cancers remain largely unknown. In this study, we performed eQTL mapping of lncRNA expression (elncRNA) in 11 TCGA cancer types and characterized the biological processes of elncRNAs in the setting of genomic location, cancer treatment responses, and immune microenvironment. As a result, 10.86% of the cis-eQTLs and 1.67% of the trans-eQTLs of lncRNA were related to known genome-wide association studies (GWAS) cancer risk loci. The elncRNAs are significantly enriched for those which are previously annotated as predictive of drug sensitivities in cancer cell lines. We further revealed the downstream transcriptomic effectors of eQTL-elncRNA pairs. Our data specifically suggested that the genes affected by eQTL-elncRNA associations are enriched in the immune system processes and eQTL-elncRNA associations influence the constitution of tumor infiltrating lymphocytes. In ovarian cancer, the "rs34631313-AC092580.4" pair was associated with increased fraction of CD8+ T cells and M1 Macrophage; whereas in KIRC, the "rs9546285-LINC00426" pair was associated with increased fraction of CD8+ T cells and a decreased fraction of M2 macrophages. Our findings provide a systematic view of the transcriptomic impacts of the eQTL landscape of lncRNA in human cancers and suggest its strong potential relevance to cancer immunity and treatment.

Genome-wide analysis of expression QTL (eQTL) and allele-specific expression (ASE) in pig muscle identifies candidate genes for meat quality traits

BACKGROUND

Genetic analysis of gene expression level is a promising approach for characterizing candidate genes that are involved in complex economic traits such as meat quality. In the present study, we conducted expression quantitative trait loci (eQTL) and allele-specific expression (ASE) analyses based on RNA-sequencing (RNAseq) data from the longissimus muscle of 189 Duroc × Luchuan crossed pigs in order to identify some candidate genes for meat quality traits.

RESULTS

Using a genome-wide association study based on a mixed linear model, we identified 7192 cis-eQTL corresponding to 2098 cis-genes (p ≤ 1.33e-3, FDR ≤ 0.05) and 6400 trans-eQTL corresponding to 863 trans-genes (p ≤ 1.13e-6, FDR ≤ 0.05). ASE analysis using RNAseq SNPs identified 9815 significant ASE-SNPs in 2253 unique genes. Integrative analysis between the cis-eQTL and ASE target genes identified 540 common genes, including 33 genes with expression levels that were correlated with at least one meat quality trait. Among these 540 common genes, 63 have been reported previously as candidate genes for meat quality traits, such as PHKG1 (q-value = 1.67e-6 for the leading SNP in the cis-eQTL analysis), NUDT7 (q-value = 5.67e-13), FADS2 (q-value = 8.44e-5), and DGAT2 (q-value = 1.24e-3).

CONCLUSIONS

The present study confirmed several previously published candidate genes and identified some novel candidate genes for meat quality traits via eQTL and ASE analyses, which will be useful to prioritize candidate genes in further studies.

Ovarian Cancer Risk Variants Are Enriched in Histotype-Specific Enhancers and Disrupt Transcription Factor Binding Sites

Quantifying the functional effects of complex disease risk variants can provide insights into mechanisms underlying disease biology. Genome-wide association studies have identified 39 regions associated with risk of epithelial ovarian cancer (EOC). The vast majority of these variants lie in the non-coding genome, where they likely function through interaction with gene regulatory elements. In this study we first estimated the heritability explained by known common low penetrance risk alleles for EOC. The narrow sense heritability (hg2) of EOC overall and high-grade serous ovarian cancer (HGSOCs) were estimated to be 5%-6%. Partitioned SNP heritability across broad functional categories indicated a significant contribution of regulatory elements to EOC heritability. We collated epigenomic profiling data for 77 cell and tissue types from Roadmap Epigenomics and ENCODE, and from H3K27Ac ChIP-seq data generated in 26 ovarian cancer and precursor-related cell and tissue types. We identified significant enrichment of risk single-nucleotide polymorphisms (SNPs) in active regulatory elements marked by H3K27Ac in HGSOCs. To further investigate how risk SNPs in active regulatory elements influence predisposition to ovarian cancer, we used motifbreakR to predict the disruption of transcription factor binding sites. We identified 469 candidate causal risk variants in H3K27Ac peaks that are predicted to significantly break transcription factor (TF) motifs. The most frequently broken motif was REST (p value = 0.0028), which has been reported as both a tumor suppressor and an oncogene. Overall, these systematic functional annotations with epigenomic data improve interpretation of EOC risk variants and shed light on likely cells of origin.

Integration of CRISPR-engineering and hiPSC-based models of psychiatric genomics

Neuropsychiatric disorders are highly heritable polygenic disorders arising from the complex interplay of highly penetrant rare variants and common variants of small effect. There is a large index of comorbidity and shared genetic risk between disorders, reflecting the pleiotropy of individual variants as well as predicted downstream pathway-level convergence. Importantly, the mechanism(s) through which psychiatric disease-associated variants interact to contribute to disease risk remains unknown. Human induced pluripotent stem cell (hiPSC)-based models are increasingly useful for the systematic study of the complex genetics associated with brain diseases, particularly when combined with CRISPR-mediated genomic engineering, which together facilitate isogenic comparisons of defined neuronal cell types. In this review, we discuss the latest CRISPR technologies and consider how they can be successfully applied to the functional characterization of the growing list genetic variants linked to psychiatric disease.

Non-coding somatic mutations converge on the PAX8 pathway in ovarian cancer

The functional consequences of somatic non-coding mutations in ovarian cancer (OC) are unknown. To identify regulatory elements (RE) and genes perturbed by acquired non-coding variants, here we establish epigenomic and transcriptomic landscapes of primary OCs using H3K27ac ChIP-seq and RNA-seq, and then integrate these with whole genome sequencing data from 232 OCs. We identify 25 frequently mutated regulatory elements, including an enhancer at 6p22.1 which associates with differential expression of ZSCAN16 (P = 6.6 × 10-4) and ZSCAN12 (P = 0.02). CRISPR/Cas9 knockout of this enhancer induces downregulation of both genes. Globally, there is an enrichment of single nucleotide variants in active binding sites for TEAD4 (P = 6 × 10-11) and its binding partner PAX8 (P = 2×10-10), a known lineage-specific transcription factor in OC. In addition, the collection of cis REs associated with PAX8 comprise the most frequently mutated set of enhancers in OC (P = 0.003). These data indicate that non-coding somatic mutations disrupt the PAX8 transcriptional network during OC development.

Identification of Candidate Genes Associated with Susceptibility to Ovarian Clear Cell Adenocarcinoma Using cis-eQTL Analysis

Ovarian clear cell adenocarcinoma (Ov-CCA) has a higher prevalence in the Japanese ancestry than other populations. The ancestral disparities in Ov-CCA prevalence suggests the presence of Ov-CCA-specific genetic alterations and may provide an opportunity to identify the novel genes associated with Ov-CCA tumorigenesis. Using 94 previously reported genes as the phenotypic trait, we conducted multistep expression quantitative trait loci (eQTL) analysis with the HapMap3 project datasets. Four single-nucleotide polymorphisms (SNPs) (rs4873815, rs12976454, rs11136002, and rs13259097) that had different allele frequencies in the Japanese ancestry and seven genes associated in cis (APBA3, C8orf58, KIAA1967, NAPRT1, RHOBTB2, TNFRSF10B, and ZNF707) were identified. In silico functional annotation analysis and in vitro promoter assay validated the regulatory effect of rs4873815-TT on ZNF707 and rs11136002-TT on TNFRSF10B. Furthermore, ZNF707 was highly expressed in Ov-CCA and had a negative prognostic value in disease recurrence in our sample cohort. This prognostic power was consistently observed in The Cancer Genome Atlas (TCGA) clear cell renal cell carcinoma dataset, suggesting that ZNF707 may have prognostic value in clear cell histology regardless of tissue origin. In conclusion, rs4873815-TT/ZNF707 may have clinical significance in the prognosis and tumorigenesis of Ov-CCA, which may be more relevant to clear cell histology. Besides, this study may underpin the evidence that cis-eQTL analysis based on ancestral disparities can facilitate the discovery of causal genetic alterations in complex diseases, such as cancer.

A Multi-Omics Perspective of Quantitative Trait Loci in Precision Medicine

Quantitative trait loci (QTL) analysis is an important approach to investigate the effects of genetic variants identified through an increasing number of large-scale, multidimensional 'omics data sets. In this 'big data' era, the research community has identified a significant number of molecular QTLs (molQTLs) and increased our understanding of their effects. Herein, we review multiple categories of molQTLs, including those associated with transcriptome, post-transcriptional regulation, epigenetics, proteomics, metabolomics, and the microbiome. We summarize approaches to identify molQTLs and to infer their causal effects. We further discuss the integrative analysis of molQTLs through a multi-omics perspective. Our review highlights future opportunities to better understand the functional significance of genetic variants and to utilize the discovery of molQTLs in precision medicine.

cis-regulatory variation modulates susceptibility to enteric infection in the Drosophila genetic reference panel

BACKGROUND

Resistance to enteric pathogens is a complex trait at the crossroads of multiple biological processes. We have previously shown in the Drosophila Genetic Reference Panel (DGRP) that resistance to infection is highly heritable, but our understanding of how the effects of genetic variants affect different molecular mechanisms to determine gut immunocompetence is still limited.

RESULTS

To address this, we perform a systems genetics analysis of the gut transcriptomes from 38 DGRP lines that were orally infected with Pseudomonas entomophila. We identify a large number of condition-specific, expression quantitative trait loci (local-eQTLs) with infection-specific ones located in regions enriched for FOX transcription factor motifs. By assessing the allelic imbalance in the transcriptomes of 19 F1 hybrid lines from a large round robin design, we independently attribute a robust cis-regulatory effect to only 10% of these detected local-eQTLs. However, additional analyses indicate that many local-eQTLs may act in trans instead. Comparison of the transcriptomes of DGRP lines that were either susceptible or resistant to Pseudomonas entomophila infection reveals nutcracker as the only differentially expressed gene. Interestingly, we find that nutcracker is linked to infection-specific eQTLs that correlate with its expression level and to enteric infection susceptibility. Further regulatory analysis reveals one particular eQTL that significantly decreases the binding affinity for the repressor Broad, driving differential allele-specific nutcracker expression.

CONCLUSIONS

Our collective findings point to a large number of infection-specific cis- and trans-acting eQTLs in the DGRP, including one common non-coding variant that lowers enteric infection susceptibility.

Candidate Causal Variants at the 8p12 Breast Cancer Risk Locus Regulate DUSP4

Genome-wide association studies have revealed a locus at 8p12 that is associated with breast cancer risk. Fine-mapping of this locus identified 16 candidate causal variants (CCVs). However, as these variants are intergenic, their function is unclear. To map chromatin looping from this risk locus to a previously identified candidate target gene, DUSP4, we performed chromatin conformation capture analyses in normal and tumoural breast cell lines. We identified putative regulatory elements, containing CCVs, which looped to the DUSP4 promoter region. Using reporter gene assays, we found that the risk allele of CCV rs7461885 reduced the activity of a DUSP4 enhancer element, consistent with the function of DUSP4 as a tumour suppressor gene. Furthermore, the risk allele of CCV rs12155535, located in another DUSP4 enhancer element, was negatively correlated with looping of this element to the DUSP4 promoter region, suggesting that this allele would be associated with reduced expression. These findings provide the first evidence that CCV risk alleles downregulate DUSP4 expression, suggesting that this gene is a regulatory target of the 8p12 breast cancer risk locus.

Characterization of genetic-origin-dependent monoallelic expression in mouse embryonic stem cells

Monoallelic gene expression occurs in various mammalian cells and can be regulated genetically, epigenetically and/or stochastically. We identified 145 monoallelically expressed genes (MoEGs), including seven known imprinted genes, in mouse embryonic stem cells (ESCs) derived from reciprocal F1 hybrid blastocysts and cultured in 2i/LIF. As all MoEGs except for the imprinted genes were expressed in a genetic-origin-dependent manner, we focused on this class of MoEGs for mechanistic studies. We showed that a majority of the genetic-origin-dependent MoEGs identified in 2i/LIF ESCs remain monoallelically expressed in serum/LIF ESCs, but become more relaxed or even biallelically expressed upon differentiation. These MoEGs and their regulatory regions were highly enriched for single nucleotide polymorphisms. In addition, some MoEGs were associated with retrotransposon insertions/deletions, consistent with the fact that certain retrotransposons act as regulatory elements in pluripotent stem cells. Interestingly, most MoEGs showed allelic differences in enrichment of histone H3K27me and H3K4me marks, linking allelic epigenetic differences and monoallelic expression. In contrast, there was little or no allelic difference in CpG methylation or H3K9me. Taken together, our study highlights the impact of genetic variation including single nucleotide polymorphisms and retrotransposon insertions/deletions on monoallelic epigenetic marks and expression in ESCs.

SNP eQTL status and eQTL density in the adjacent region of the SNP are associated with its statistical significance in GWA studies

Background

Over the relatively short history of Genome Wide Association Studies (GWASs), hundreds of GWASs have been published and thousands of disease risk-associated SNPs have been identified. Summary statistics from the conducted GWASs are often available and can be used to identify SNP features associated with the level of GWAS statistical significance. Those features could be used to select SNPs from gray zones (SNPs that are nominally significant but do not reach the genome-wide level of significance) for targeted analyses.

Methods

We used summary statistics from recently published breast and lung cancer and scleroderma GWASs to explore the association between the level of the GWAS statistical significance and the expression quantitative trait loci (eQTL) status of the SNP. Data from the Genotype-Tissue Expression Project (GTEx) were used to identify eQTL SNPs.

Results

We found that SNPs reported as eQTLs were more significant in GWAS (higher -log₁₀p) regardless of the tissue specificity of the eQTL. Pan-tissue eQTLs (those reported as eQTLs in multiple tissues) tended to be more significant in the GWAS compared to those reported as eQTL in only one tissue type. eQTL density in the ±5 kb adjacent region of a given SNP was also positively associated with the level of GWAS statistical significance regardless of the eQTL status of the SNP. We found that SNPs located in the regions of high eQTL density were more likely to be located in regulatory elements (transcription factor or miRNA binding sites).

When SNPs were stratified by the level of statistical significance, the proportion of eQTLs was positively associated with the mean level of statistical significance in the group. The association curve reaches a plateau around -log₁₀p ≈ 5. The observed associations suggest that quasi-significant SNPs (10^− 5 < p < 5 × 10^− 8) and SNPs at the genome wide level of statistical significance (p < 5 × 10^− 8) may have a similar proportions of risk associated SNPs.

Conclusions

The results of this study indicate that the SNP’s eQTL status, as well as eQTL density in the adjacent region are positively associated with the level of statistical significance of the SNP in GWAS.

The prospect of discovering new biomarkers for ovarian cancer based on current knowledge of susceptibility loci and genetic variation (Review)

Ovarian cancer is the most lethal gynaecological malignancy. The cancer initially presents with non‑specific symptoms; thus, it is typically not discovered until the patient has reached the late, considerably more lethal, stages of the disease. Research focus is currently on finding novel biomarkers, especially for early detection and stratification of the disease. One promising approach has been to focus on mutations or variations in the genetic code that are associated with the risk of developing ovarian cancer. A certain heritable component is already known regarding genes such as BRCA1/2, TP53, MSH6, BRIP1 and RAD51C, yet these are estimated to only account for ~3.1% of the total risk. Recent advances in sequencing technologies have enabled the investigation of hundreds of thousands of genetic variants in genome‑wide association studies in tens of thousands of patients, which has led to the discovery of 108 (39 loci with P<5.0x10‑8) novel susceptibility loci for ovarian cancer, presented in this review. Using the published variants in a patient cohort screening, together with variants identified in our ongoing whole exome sequencing project, future aims are to ascertain whether certain of the novel variants could be used as biomarkers for early diagnosis and/or treatment decisions.

A transcriptome-wide association study of high-grade serous epithelial ovarian cancer identifies new susceptibility genes and splice variants

We sought to identify susceptibility genes for high-grade serous ovarian cancer (HGSOC) by performing a transcriptome-wide association study of gene expression and splice junction usage in HGSOC-relevant tissue types (N = 2,169) and the largest genome-wide association study available for HGSOC (N = 13,037 cases and 40,941 controls). We identified 25 transcriptome-wide association study significant genes, 7 at the junction level only, including LRRC46 at 19q21.32, (P = 1 × 10-9), CHMP4C at 8q21 (P = 2 × 10-11) and a PRC1 junction at 15q26 (P = 7 × 10-9). In vitro assays for CHMP4C showed that the associated variant induces allele-specific exon inclusion (P = 0.0024). Functional screens in HGSOC cell lines found evidence of essentiality for three of the new genes we identified: HAUS6, KANSL1 and PRC1, with the latter comparable to MYC. Our study implicates at least one target gene for 6 out of 13 distinct genome-wide association study regions, identifying 23 new candidate susceptibility genes for HGSOC.

Redefining the IBDs using genome-scale molecular phenotyping

The IBDs, Crohn's disease and ulcerative colitis, are chronic inflammatory conditions of the gastrointestinal tract resulting from an aberrant immune response to enteric microbiota in genetically susceptible individuals. Disease presentation and progression within and across IBDs, especially Crohn's disease, are highly heterogeneous in location, severity of inflammation and other phenotypes. Current clinical classifications fail to accurately predict disease course and response to therapies. Genome-wide association studies have identified >240 loci that confer risk of IBD, but the clinical utility of these findings remains unclear, and mechanisms by which the genetic variants contribute to disease are largely unknown. In the past 5 years, the profiling of genome-wide gene expression, epigenomic features and gut microbiota composition in intestinal tissue and faecal samples has uncovered distinct molecular signatures that define IBD subtypes, including within Crohn's disease and ulcerative colitis. In this Review, we summarize studies in both adult and paediatric patients that have identified different IBD subtypes, which in some cases have been associated with distinct clinical phenotypes. We posit that genome-scale molecular phenotyping in large cohorts holds great promise not only to further our understanding of the diverse molecular causes of IBD but also for improving clinical trial design to develop more personalized disease management and treatment.

The Role of HOX Transcription Factors in Cancer Predisposition and Progression

Homeobox (HOX) transcription factors, encoded by a subset of homeodomain superfamily genes, play pivotal roles in many aspects of cellular physiology, embryonic development, and tissue homeostasis. Findings over the past decade have revealed that mutations in HOX genes can lead to increased cancer predisposition, and HOX genes might mediate the effect of many other cancer susceptibility factors by recognizing or executing altered genetic information. Remarkably, several lines of evidence highlight the interplays between HOX transcription factors and cancer risk loci discovered by genome-wide association studies, thereby gaining molecular and biological insight into cancer etiology. In addition, deregulated HOX gene expression impacts various aspects of cancer progression, including tumor angiogenesis, cell autophagy, proliferation, apoptosis, tumor cell migration, and metabolism. In this review, we will discuss the fundamental roles of HOX genes in cancer susceptibility and progression, highlighting multiple molecular mechanisms of HOX involved gene misregulation, as well as their potential implications in clinical practice.

eQTL analysis from co-localization of 2739 GWAS loci detects associated genes across 14 human cancers

Genetic variants can predict other "linked" diseases because alterations in one or more genes in vivo may affect relevant phenotype properties. Our study systematically explored the pan-cancer common gene and cancer type-specific genes based on GWAS loci and TCGA data of multiple cancers. It was found that there were 17 SNPs were significantly associated with the expression of 18 genes. Associations between the 18 cis-regulatory genes and the pathologic stage of each cancer showed that MYL2 and PTGFR in HNSC, 4 genes (F8, SATB2, G6PD and UGT1A6) in KIRP, 3 genes (CHMP4C, MAP3K1 and MECP2) in LUAD were all strongly associated with cancer stage levels. Additionally, the survival association analysis showed that SATB2 was correlated with HNSC survival, and MPP1 was strongly associated with the survival of SARC. This study will shed light on the biological pathways involved in cancer-genetic associations, and has the potential to be applied to the predictions of the risk of cancers developing in healthy individuals.

Maximizing the Utility of Cancer Transcriptomic Data

Transcriptomic profiling has been applied to large numbers of cancer samples, by large-scale consortia, including The Cancer Genome Atlas, International Cancer Genome Consortium, and Cancer Cell Line Encyclopedia. Advances in mining cancer transcriptomic data enable us to understand the endless complexity of the cancer transcriptome and thereby to discover new biomarkers and therapeutic targets. In this paper, we review computational resources for deep mining of transcriptomic data to identify, quantify, and determine the functional effects and clinical utility of transcriptomic events, including noncoding RNAs, post-transcriptional regulation, exogenous RNAs, and transcribed genetic variants. These approaches can be applied to other complex diseases, thereby greatly leveraging the impact of this work.

Exploring the Role of Fallopian Ciliated Cells in the Pathogenesis of High-Grade Serous Ovarian Cancer

High-grade serous epithelial ovarian cancer (HGSOC) is the fifth leading cause of cancer death in women and the first among gynecological malignancies. Despite an initial response to standard chemotherapy, most HGSOC patients relapse. To improve treatment options, we must continue investigating tumor biology. Tumor characteristics (e.g., risk factors and epidemiology) are valuable clues to accomplish this task. The two most frequent risk factors for HGSOC are the lifetime number of ovulations, which is associated with increased oxidative stress in the pelvic area caused by ovulation fluid, and a positive family history due to genetic factors. In the attempt to identify novel genetic factors (i.e., genes) associated with HGSOC, we observed that several genes in linkage with HGSOC are expressed in the ciliated cells of the fallopian tube. This finding made us hypothesize that ciliated cells, despite not being the cell of origin for HGSOC, may take part in HGSOC tumor initiation. Specifically, malfunction of the ciliary beat impairs the laminar fluid flow above the fallopian tube epithelia, thus likely reducing the clearance of oxidative stress caused by follicular fluid. Herein, we review the up-to-date findings dealing with HGSOC predisposition with the hypothesis that fallopian ciliated cells take part in HGSOC onset. Finally, we review the up-to-date literature concerning genes that are located in genomic loci associated with epithelial ovarian cancer (EOC) predisposition that are expressed by the fallopian ciliated cells.

Identification of Genes and Pathways Involved in Ovarian Epithelial Cancer by Bioinformatics Analysis

Ovarian epithelial cancer (OEC) is an often fatal disease with poor prognosis in women with high-stage disease. In contrast, ovarian low malignant potential (LMP) tumors with favorable prognosis behaves as a disease between benign and malignant tumors. The involved genes and pathways between benign-like LMP and aggressive OEC are largely unknown. This study integrated two cohorts profile datasets to investigate the potential key candidate genes and pathways associated with OEC. Gene expression in two datasets (GSE9891 and GSE12172), including 327 OECs and 48 LMP tumors, were analyzed. 559 differentially expressed genes were found to overlap, 251 up-regulated and 308 down-regulated. Subsequently, analysis of gene ontology, signaling pathway enrichment and protein-protein interaction (PPI) network was performed. Gene ontology analysis clustered the up-regulated and down-regulated genes based on significant enrichment. 282 nodes/ differentially expressed genes (DEGs) were identified from DEGs PPI network complex, and two most significant k-clique modules were identified from PPI. In a summary, using integrated bioinformatics analysis, we are able to identify biomarkers potentially significant in the pathogenesis of OEC, which can improve our understanding of the cause and molecular events. These candidate genes and pathways could be used for further confirmation, and lead to better disease diagnose and therapy.

Early transcriptional response of human ovarian and fallopian tube surface epithelial cells to norepinephrine

Evidence from human and animal studies suggests that chronic behavioral stress and resulting activation of the sympathetic nervous system may influence initiation and progression of tumors. However, the underlying mechanisms for these observations are poorly understood. The purpose of this study is to explore the effects of adrenergic signaling on cell line models derived from normal cells presumed to originate epithelial ovarian cancers. Here we explored the effects of the stress-related hormone, norepinephrine, on the transcriptional program of normal immortalized ovarian (iOSE) and fallopian tube (iFTSEC) surface epithelial cells. Analysis of RNA-Seq data of treated and untreated cells revealed a significant overlap between the responses in iOSE and iFTSEC cells. Most genes modulated by norepinephrine in ovarian and fallopian tube epithelial cells are already expressed in normal ovarian and fallopian tissue and cells. For several genes, expression changes were reflected at the protein level. Genes in immune-related and developmental pathways were enriched in the set of genes modulated by norepinephrine. We identified HOXA5, SPIB, REL, SRF, SP1, NFKB1, MEF2A, E2F1, and EGR1 transcription factor binding sites to be highly enriched in our dataset. These data represent the early transcriptional response to norepinephrine in cells postulated to originate epithelial ovarian cancer.

Common Genetic Variation and Susceptibility to Ovarian Cancer: Current Insights and Future Directions

In this review, we summarize current progress in the genetic epidemiology of epithelial ovarian cancer (EOC), focusing exclusively on elucidating the role of common germline genetic variation in conferring susceptibility to EOC. We provide an overview of the more than 30 EOC risk loci identified to date by genome-wide association studies (GWAS) and describe the contribution of large-scale, cross-cancer type, custom genotyping projects, such as the OncoArray and the Collaborative Oncological Gene-Environment Study, to locus discovery and replication. We discuss the histotype-specific nature of these EOC risk loci, pleiotropy, or overlapping genetic effects between EOC and other hormone-related cancer types, and the application of findings to polygenic risk prediction for EOC. The second part of the article offers a concise review of primarily laboratory-based studies that have led to the identification of several putative EOC susceptibility genes using common variants at the known EOC risk loci as starting points. More global biological insights emerging from network- and pathway-based analyses of GWAS for EOC susceptibility are also highlighted. Finally, we delve into potential future directions, including the need to identify EOC risk loci in non-European populations and the next generation of GWAS functional studies that are likely to involve genome editing to establish the cell type-specific carcinogenic effects of EOC risk variants Cancer Epidemiol Biomarkers Prev; 27(4); 395-404. ©2018 AACRSee all articles in this CEBP Focus section, "Genome-Wide Association Studies in Cancer."

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

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

Genetic epidemiology of ovarian cancer and prospects for polygenic risk prediction

Epithelial ovarian cancer (EOC) is a heterogeneous disease with a major heritable component. The different histotypes of invasive disease - high grade serous, clear cell, endometrioid and mucinous - are associated with different underlying genetic susceptibility and epidemiological and lifestyle risk factors, all of which contribute to the different biology and clinical characteristics of each histotype. A combination of familial and population based sequencing studies, and genome wide association studies (GWAS) have identified a range of genetic susceptibility alleles for EOC comprising rare but highly penetrant genes (e.g. BRCA1, BRCA2) that are responsible for familial clustering of ovarian cancer cases; more moderate penetrance susceptibility genes (e.g. BRIP1, RAD51C/D, MSH6); and multiple common but low penetrance susceptibility alleles identified by GWAS. Identifying genetic risk alleles for ovarian cancer has had a significant impact on disease prevention strategies; for example it is now routine clinical practice for individuals with germline BRCA1 and BRCA2 mutations to undergo risk reducing salpingo-oophorectomy. Because ovarian cancers are commonly diagnosed at a late clinical stage when the prognosis is poor, the continued development of genetic risk prediction and prevention strategies will represent an important approach to reduce mortality due to ovarian cancer. Advances in genomics technologies that enable more high-throughput genetic testing, combined with research studies that identify additional EOC risk alleles will likely provide further opportunities to establish polygenic risk prediction approaches, based on combinations of rare high/moderate penetrance susceptibility genes and common, low penetrance susceptibility alleles. This article reviews the current literature describing the genetic and epidemiological components of ovarian cancer risk, and discusses both the opportunities and challenges in using this information for clinical risk prediction and prevention.

Genome-wide association studies of cancer: current insights and future perspectives

Genome-wide association studies (GWAS) provide an agnostic approach for investigating the genetic basis of complex diseases. In oncology, GWAS of nearly all common malignancies have been performed, and over 450 genetic variants associated with increased risks have been identified. As well as revealing novel pathways important in carcinogenesis, these studies have shown that common genetic variation contributes substantially to the heritable risk of many common cancers. The clinical application of GWAS is starting to provide opportunities for drug discovery and repositioning as well as for cancer prevention. However, deciphering the functional and biological basis of associations is challenging and is in part a barrier to fully unlocking the potential of GWAS.

Enrichment of putative PAX8 target genes at serous epithelial ovarian cancer susceptibility loci

BACKGROUND

Genome-wide association studies (GWAS) have identified 18 loci associated with serous ovarian cancer (SOC) susceptibility but the biological mechanisms driving these findings remain poorly characterised. Germline cancer risk loci may be enriched for target genes of transcription factors (TFs) critical to somatic tumorigenesis.

METHODS

All 615 TF-target sets from the Molecular Signatures Database were evaluated using gene set enrichment analysis (GSEA) and three GWAS for SOC risk: discovery (2196 cases/4396 controls), replication (7035 cases/21 693 controls; independent from discovery), and combined (9627 cases/30 845 controls; including additional individuals).

RESULTS

The PAX8-target gene set was ranked 1/615 in the discovery (PGSEA<0.001; FDR=0.21), 7/615 in the replication (PGSEA=0.004; FDR=0.37), and 1/615 in the combined (PGSEA<0.001; FDR=0.21) studies. Adding other genes reported to interact with PAX8 in the literature to the PAX8-target set and applying an alternative to GSEA, interval enrichment, further confirmed this association (P=0.006). Fifteen of the 157 genes from this expanded PAX8 pathway were near eight loci associated with SOC risk at P<10-5 (including six with P<5 × 10-8). The pathway was also associated with differential gene expression after shRNA-mediated silencing of PAX8 in HeyA8 (PGSEA=0.025) and IGROV1 (PGSEA=0.004) SOC cells and several PAX8 targets near SOC risk loci demonstrated in vitro transcriptomic perturbation.

CONCLUSIONS

Putative PAX8 target genes are enriched for common SOC risk variants. This finding from our agnostic evaluation is of particular interest given that PAX8 is well-established as a specific marker for the cell of origin of SOC.

The OncoArray Consortium: A Network for Understanding the Genetic Architecture of Common Cancers

BACKGROUND

Common cancers develop through a multistep process often including inherited susceptibility. Collaboration among multiple institutions, and funding from multiple sources, has allowed the development of an inexpensive genotyping microarray, the OncoArray. The array includes a genome-wide backbone, comprising 230,000 SNPs tagging most common genetic variants, together with dense mapping of known susceptibility regions, rare variants from sequencing experiments, pharmacogenetic markers, and cancer-related traits.

METHODS

The OncoArray can be genotyped using a novel technology developed by Illumina to facilitate efficient genotyping. The consortium developed standard approaches for selecting SNPs for study, for quality control of markers, and for ancestry analysis. The array was genotyped at selected sites and with prespecified replicate samples to permit evaluation of genotyping accuracy among centers and by ethnic background.

RESULTS

The OncoArray consortium genotyped 447,705 samples. A total of 494,763 SNPs passed quality control steps with a sample success rate of 97% of the samples. Participating sites performed ancestry analysis using a common set of markers and a scoring algorithm based on principal components analysis.

CONCLUSIONS

Results from these analyses will enable researchers to identify new susceptibility loci, perform fine-mapping of new or known loci associated with either single or multiple cancers, assess the degree of overlap in cancer causation and pleiotropic effects of loci that have been identified for disease-specific risk, and jointly model genetic, environmental, and lifestyle-related exposures.

IMPACT

Ongoing analyses will shed light on etiology and risk assessment for many types of cancer. Cancer Epidemiol Biomarkers Prev; 26(1); 126-35. ©2016 AACR.

Genetically regulated hepatic transcripts and pathways orchestrate haematological, biochemical and body composition traits

The liver is the central metabolic organ and exhibits fundamental functions in haematological traits. Hepatic expression, haematological, plasma biochemical, and body composition traits were assessed in a porcine model (n = 297) to establish tissue-specific genetic variations that influence the function of immune-metabolism-correlated expression networks. At FDR (false discovery rate) <1%, more than 3,600 transcripts were jointly correlated (r = |0.22-0.48|) with the traits. Functional enrichment analysis demonstrated common links of metabolic and immune traits. To understand how immune and metabolic traits are affected via genetic regulation of gene expression, eQTLs were assessed. 20517 significant (FDR < 5%) eQTLs for 1401 transcripts were identified, among which 443 transcripts were associated with at least one of the examined traits and had cis-eQTL (such as ACO1 (6.52 × 10^-7) and SOD1 (6.41 × 10^-30). The present study establishes a comprehensive view of hepatic gene activity which links together metabolic and immune traits in a porcine model for medical research.

Inherited variants affecting RNA editing may contribute to ovarian cancer susceptibility: results from a large-scale collaboration

RNA editing in mammals is a form of post-transcriptional modification in which adenosine is converted to inosine by the adenosine deaminases acting on RNA (ADAR) family of enzymes. Based on evidence of altered ADAR expression in epithelial ovarian cancers (EOC), we hypothesized that single nucleotide polymorphisms (SNPs) in ADAR genes modify EOC susceptibility, potentially by altering ovarian tissue gene expression. Using directly genotyped and imputed data from 10,891 invasive EOC cases and 21,693 controls, we evaluated the associations of 5,303 SNPs in ADAD1, ADAR, ADAR2, ADAR3, and SND1. Unconditional logistic regression was used to estimate odds ratios (OR) and 95% confidence intervals (CI), with adjustment for European ancestry. We conducted gene-level analyses using the Admixture Maximum Likelihood (AML) test and the Sequence-Kernel Association test for common and rare variants (SKAT-CR). Association analysis revealed top risk-associated SNP rs77027562 (OR (95% CI)= 1.39 (1.17-1.64), P=1.0x10-4) in ADAR3 and rs185455523 in SND1 (OR (95% CI)= 0.68 (0.56-0.83), P=2.0x10-4). When restricting to serous histology (n=6,500), the magnitude of association strengthened for rs185455523 (OR=0.60, P=1.0x10-4). Gene-level analyses revealed that variation in ADAR was associated (P<0.05) with EOC susceptibility, with PAML=0.022 and PSKAT-CR=0.020. Expression quantitative trait locus analysis in EOC tissue revealed significant associations (P<0.05) with ADAR expression for several SNPs in ADAR, including rs1127313 (G/A), a SNP in the 3' untranslated region. In summary, germline variation involving RNA editing genes may influence EOC susceptibility, warranting further investigation of inherited and acquired alterations affecting RNA editing.