Integrative Functional Genomics Research Articles

Genome-wide functional integration identified MAZ-controlled RPS14 dysregulation in hepatocellular carcinoma.

Chronic infection with Hepatitis B virus (HBV) significantly increases the risk of hepatocellular carcinoma (HCC), particularly in Eastern Asia. However, only a subset of individuals with chronic HBV infection develop HCC, suggesting the role for genetic factors in HCC etiology. Despite genome-wide association studies (GWASs) identifying multiple single nucleotide polymorphisms (SNPs) associated with HBV-related HCC susceptibility, the underlying mechanisms and causal genetic polymorphisms remain largely unclear. To address this, we developed The Updated Integrative Functional Genomics Approach (TUIFGA), an methodology that combines data from transcription factor (TF) cistromics, ATAC-seq, DNAase-seq, and the 1000 Genomes Project to identify cancer susceptibility SNPs within TF-binding sites across human genome. Using TUIFGA, we discovered SNP rs13170300 which located in the TF MAZ binding motif of RPS14. The RPS14 rs13170300 was significantly associated with HCC risk in two case-control sets, with the T allele as the protective allele (Shandong discovery set: TT OR = 0.60, 95% CI = 0.49-0.74, P = 1.0 × 10-6; CT OR = 0.69, 95% CI = 0.55-0.86, P = 0.001; Jiangsu validation set: TT OR = 0.70, 95% CI = 0.56-0.87, P = 0.001; CT OR = 0.65, 95% CI = 0.53-0.82, P = 1.6 × 10-4). SNP rs13170300 affected MAZ binding in the RPS14 promoter, resulting in allele-specific changes in gene expression. RPS14 functions as a novel oncogene in HCC, specifically via activating the AKT signaling. Our findings present important insights into the functional genetics underlying HBV-related HCC development and may contribute to personalized approaches for cancer prevention and novel therapeutics.

DEK Regulates B-Cell Proliferative Capacity and Is Associated with Aggressive Disease in Low-Grade B-Cell Lymphoma

Introduction: Non-Follicular Low-Grade B-Cell Lymphomas (LGBCL) are a subset of lymphomas that share overlapping clinical and histologic features. While most cases of LGBCL have an indolent course, some experience aggressive disease with no definitive biomarkers to identify aggressive cases at diagnosis. Significant advances have been made in the molecular classification of other lymphomas, yet molecular profiling and identification of mechanisms driving aggressive disease in LGBCL remain an unmet need. We previously used a multi-omic approach to characterize the molecular and immune landscapes of LGBCL tumors and identified a transcriptomic signature associated with cell proliferation and inferior overall survival (OS) in two independent cohorts (HR 7.82; 95% CI 2.4-25.4; p < 0.001 & HR 10.07; 95% CI 2.00-50.61; p = 0.005) and cases that transformed to DLBCL (Hopper et al., 2023). Here, we identify the oncoprotein DEK as a novel transcriptional regulator of this signature and use a functional genomic approach to further interrogate the role of DEK in B-cell lymphoma. Methods: Tumors from 64 newly-diagnosed/untreated LGBCL cases enrolled in the University of Iowa/Mayo Clinic Molecular Epidemiology Resource were included in this study (SMZL = 48, NMZL = 6, LPL = 5, B-NOS = 3, EMZL = 2). Tumor RNA-seq and paired tumor-normal WES were used to assess DEK expression and identify copy number alterations (CNAs), respectively. DEK deficient (KO) cell line models were generated using CRIPSR/Cas9. DEK expression was assessed by western blot (WB) and immunofluorescence (IF) using standard techniques. Proliferation was assessed by MTT assay. Results: To identify transcriptional regulators of our previously identified signature, we used the tool MAGIC to mine publicly available ChIP-seq data and identified DEK, an oncoprotein associated with tumorigenesis and poor prognosis in several cancers, as a key regulator of genes in our signature (p = 0.008). RNA-seq indicated that high levels of DEK expression associated with inferior OS in LGBCL (HR 2.99; 95% CI 1.21-7.36; p = 0.02). Additionally, we identified DEK as a target of copy number gains/amps in LGBCL (7% of cases), which corresponded to significantly higher DEK expression (p = 0.04). The association of DEK with aggressive lymphoma was confirmed when we assessed its expression in public cell line RNA-seq data and observed the highest expression in lines derived from aggressive lymphomas, which we further validated at the protein level in a panel of B-cell lymphoma cell lines by WB and IF confocal microscopy. Additionally, we detected nuclear localization of DEK, consistent with its role as a transcriptional regulator. As a model of DEK overexpression for functional characterization, we used the SUDHL6 lymphoma cell line, which has an amplification of DEK at 6p22.3 and expresses high levels of DEK protein. Using CRISPR technology, we generated DEK deficient SUDHL6 and tested the impact of loss of DEK expression. DEK KO cells proliferated significantly slower than WT SUDHL6 cells in vitro (p < 0.0001)(Panel A), indicating DEK plays a role in regulating proliferative capacity in lymphoma. To identify mechanisms responsible for reduced proliferation in DEK KOs, we assessed baseline protein expression levels of apoptosis proteins in WT versus KO cells. We found near complete abrogation of anti-apoptotic protein Bcl-2 and reduction of Bcl-xL expression in DEK KO cells (Panel B). These results indicate DEK may contribute to tumor proliferative capacity by regulating susceptibility to apoptosis. Conclusion: DEK has been implicated in regulation of gene expression, but its precise role has remained elusive, with little known about its role in lymphoma. Using an integrative molecular and functional genomics approach, we identified DEK as a novel oncoprotein in LGBCL that is upregulated in cases with inferior clinical outcome and is regulated by copy number gains/amps. DEK plays a role in regulating tumor proliferative capacity by altering susceptibility to apoptosis through regulation of Bcl-2 and Bcl-xL. Collectively, this data suggests high-risk LGBCL tumors may enrich for cells with DEK overexpression or CNAs as a mechanism driving oncogene addiction, highlighting DEK, and its downstream targets, as potential therapeutic targets in a subset of high-risk LGBCL.

PPM1D mutations are oncogenic drivers of de novo diffuse midline glioma formation

The role of PPM1D mutations in de novo gliomagenesis has not been systematically explored. Here we analyze whole genome sequences of 170 pediatric high-grade gliomas and find that truncating mutations in PPM1D that increase the stability of its phosphatase are clonal driver events in 11% of Diffuse Midline Gliomas (DMGs) and are enriched in primary pontine tumors. Through the development of DMG mouse models, we show that PPM1D mutations potentiate gliomagenesis and that PPM1D phosphatase activity is required for in vivo oncogenesis. Finally, we apply integrative phosphoproteomic and functional genomics assays and find that oncogenic effects of PPM1D truncation converge on regulators of cell cycle, DNA damage response, and p53 pathways, revealing therapeutic vulnerabilities including MDM2 inhibition.

Identification of genetic variants influencing methylation in brain with pleiotropic effects on psychiatric disorders

Psychiatric disorders affect 29% of the global population at least once in the lifespan, and genetic studies have proved a shared genetic basis among them, although the underlying molecular mechanisms remain largely unknown. DNA methylation plays an important role in complex disorders and, remarkably, enrichment of common genetic variants influencing allele-specific methylation (ASM) has been reported among variants associated with specific psychiatric disorders.In the present study we assessed the contribution of ASM to a set of eight psychiatric disorders by combining genetic, epigenetic and expression data. We interrogated a list of 3896 ASM tagSNPs in the brain in the summary statistics of a cross-disorder GWAS meta-analysis of eight psychiatric disorders from the Psychiatric Genomics Consortium, including more than 162,000 cases and 276,000 controls. We identified 80 SNPs with pleiotropic effects on psychiatric disorders that show an opposite directional effect on methylation and gene expression. These SNPs converge on eight candidate genes: ZSCAN29, ZSCAN31, BTN3A2, DDAH2, HAPLN4, ARTN, FAM109B and NAGA. ZSCAN29 shows the broadest pleiotropic effects, showing associations with five out of eight psychiatric disorders considered, followed by ZSCAN31 and BTN3A2, associated with three disorders. All these genes overlap with CNVs related to cognitive phenotypes and psychiatric traits, they are expressed in the brain, and seven of them have previously been associated with specific psychiatric disorders, supporting our results. To sum up, our integrative functional genomics analysis identified eight psychiatric disease risk genes that impact a broad list of disorders and highlight an etiologic role of SNPs that influence DNA methylation and gene expression in the brain.

Integrative Functional Genomics Implicated the Key T-/B-Cell Deficiency Regulator RAG1 in Transarterial Chemoembolization of Hepatocellular Carcinoma.

Transarterial chemoembolization (TACE) has significantly prolonged overall survival (OS) of unresectable hepatocellular carcinoma (HCC) patients. Unfortunately, there are still a portion of patients without therapeutic responses to TACE. Although genome-wide association studies identified multiple HCC susceptibility SNPs, it is still largely unclear how genome-wide identified functional SNPs impacting gene expression contribute to the prognosis of TACE-treated HCC patients. In this study, we developed an integrative functional genomics methodology to identify gene expression-related SNPs significantly contributing to prognosis of TACE-treated HCC patients across the whole genome. Employing integration of data from expression quantitative trait locus (eQTLs) analyses of The Cancer Genome Atlas (TCGA) liver hepatocellular carcinoma (LIHC) as well as the 1000 Genomes project, we successfully annotated 60 gene expression-related SNPs which are associated with OS of the TCGA patients. After genotyping these 60 SNPs in our TACE cohort, we identified four SNPs (rs12574873, rs12513391, rs34597395, and rs35624901) which are significantly associated with OS of HCC patients treated with TACE. For instance, multivariate Cox proportional hazards model indicated that the rs35624901 Deletion.Deletion (Del.Del) genotype carriers had markedly prolonged OS and a 55% decreased death risk compared with individuals with the GG genotype after TACE therapy (p = 8.3 × 10–5). In support of this, the rs35624901 Del.Del genotype is correlated to higher expression of RAG1, a key T-/B-cell deficiency regulator. Our findings reported the first evidence supporting the prognostic value of four eQTL SNPs in TACE-treated HCC patients. Importantly, our data implicated that antitumor immunity might contribute to TACE efficiency for unresectable HCC patients.

In vivo Functional Genomics for Undiagnosed Patients: The Impact of Small GTPases Signaling Dysregulation at Pan-Embryo Developmental Scale.

While individually rare, disorders affecting development collectively represent a substantial clinical, psychological, and socioeconomic burden to patients, families, and society. Insights into the molecular mechanisms underlying these disorders are required to speed up diagnosis, improve counseling, and optimize management toward targeted therapies. Genome sequencing is now unveiling previously unexplored genetic variations in undiagnosed patients, which require functional validation and mechanistic understanding, particularly when dealing with novel nosologic entities. Functional perturbations of key regulators acting on signals’ intersections of evolutionarily conserved pathways in these pathological conditions hinder the fine balance between various developmental inputs governing morphogenesis and homeostasis. However, the distinct mechanisms by which these hubs orchestrate pathways to ensure the developmental coordinates are poorly understood. Integrative functional genomics implementing quantitative in vivo models of embryogenesis with subcellular precision in whole organisms contribute to answering these questions. Here, we review the current knowledge on genes and mechanisms critically involved in developmental syndromes and pediatric cancers, revealed by genomic sequencing and in vivo models such as insects, worms and fish. We focus on the monomeric GTPases of the RAS superfamily and their influence on crucial developmental signals and processes. We next discuss the effectiveness of exponentially growing functional assays employing tractable models to identify regulatory crossroads. Unprecedented sophistications are now possible in zebrafish, i.e., genome editing with single-nucleotide precision, nanoimaging, highly resolved recording of multiple small molecules activity, and simultaneous monitoring of brain circuits and complex behavioral response. These assets permit accurate real-time reporting of dynamic small GTPases-controlled processes in entire organisms, owning the potential to tackle rare disease mechanisms.

Targeting CDK4/6 in mantle cell lymphoma.

Targeting the cell cycle represents a rational approach to mantle cell lymphoma (MCL) therapy, as aberrant expression of cyclin D1 and dysregulation of CDK4 underlie cell cycle progression and proliferation of MCL cells. Although cell cycle cancer therapy was historically ineffective due to a lack of selective and effective drugs, this landscape changed with the advent of selective and potent small-molecule oral CDK4/6 inhibitors. Here, we review the anti-tumor activities and clinical data of selective CDK4/6 inhibitors in MCL. We summarize the known mechanism of action of palbociclib, the most specific CDK4/6 inhibitor to date, and the strategy to leverage this specificity to reprogram MCL for a deeper and more durable clinical response to partner drugs. We also discuss integrative longitudinal functional genomics as a strategy to discover tumor-intrinsic genomic biomarkers and tumor-immune interactions that potentially contribute to the clinical response to palbociclib in combination therapy for MCL. Understanding the genomic basis for targeting CDK4/6 and the mechanisms of action and resistance in MCL may advance personalized therapy for MCL and shed light on drug resistance in other cancers.

Cisplatin-resistant triple-negative breast cancer subtypes: multiple mechanisms of resistance

BackgroundUnderstanding mechanisms underlying specific chemotherapeutic responses in subtypes of cancer may improve identification of treatment strategies most likely to benefit particular patients. For example, triple-negative breast cancer (TNBC) patients have variable response to the chemotherapeutic agent cisplatin. Understanding the basis of treatment response in cancer subtypes will lead to more informed decisions about selection of treatment strategies.MethodsIn this study we used an integrative functional genomics approach to investigate the molecular mechanisms underlying known cisplatin-response differences among subtypes of TNBC. To identify changes in gene expression that could explain mechanisms of resistance, we examined 102 evolutionarily conserved cisplatin-associated genes, evaluating their differential expression in the cisplatin-sensitive, basal-like 1 (BL1) and basal-like 2 (BL2) subtypes, and the two cisplatin-resistant, luminal androgen receptor (LAR) and mesenchymal (M) subtypes of TNBC.ResultsWe found 20 genes that were differentially expressed in at least one subtype. Fifteen of the 20 genes are associated with cell death and are distributed among all TNBC subtypes. The less cisplatin-responsive LAR and M TNBC subtypes show different regulation of 13 genes compared to the more sensitive BL1 and BL2 subtypes. These 13 genes identify a variety of cisplatin-resistance mechanisms including increased transport and detoxification of cisplatin, and mis-regulation of the epithelial to mesenchymal transition.ConclusionsWe identified gene signatures in resistant TNBC subtypes indicative of mechanisms of cisplatin. Our results indicate that response to cisplatin in TNBC has a complex foundation based on impact of treatment on distinct cellular pathways. We find that examination of expression data in the context of heterogeneous data such as drug-gene interactions leads to a better understanding of mechanisms at work in cancer therapy response.

Epithelial Function in the Drosophila Malpighian Tubule: An In Vivo Renal Model.

The insect renal (Malpighian) tubule has long been a model system for the study of fluid secretion and its neurohormonal control, as well as studies on ion transport mechanisms. To extend these studies beyond the boundaries of classical physiology, a molecular genetic approach together with the 'omics technologies is required. To achieve this in any vertebrate transporting epithelium remains a daunting task, as the genetic tools available are still relatively unsophisticated. Drosophila melanogaster, however, is an outstanding model organism for molecular genetics. Here we describe a technique for fluid secretion assays in the D. melanogaster equivalent of the kidney nephron. The development of this first physiological assay for a Drosophila epithelium, allowing combined approaches of integrative physiology and functional genomics, has now provided biologists with an entirely new model system, the Drosophila Malpighian tubule, which is utilized in multiple fields as diverse as kidney disease research and development of new modes of pest insect control.

Integrative functional genomics identifies regulatory genetic variant modulating RAB31 expression and altering susceptibility to breast cancer.

Despite the successes of genome-wide association study (GWAS) in identifying breast cancer (BC) risk-associated variants, only a small fraction of the heritability can be explained. The greatest challenge in the post-GWAS is to identify causal variants and underlying mechanisms responsible for BC susceptibility. In this study, we integrated functional genomic data from ENCODE ChIP-seq, ANNOVAR, and the TRANSFAC matrix to identify potentially regulatory variants with modulating FOXA1-binding affinity across the whole genome, and then conducted a two-stage case-control study including 2164 cases and 2382 controls to investigate the associations between candidate SNPs and BC susceptibility. We identified a BC susceptibility SNP, rs6506689 G>T, with an odds ratio (OR) of 1.23 (95% confidence interval = 1.07-1.40, P = 0.003) under a dominant model in the combined study. Biological assays indicated that the germline G>T variation at rs6506689 creates a FOXA1-binding site and up-regulates the expression of RAB31, thus playing an important role in the development of BC. Our results highlight the importance of regulatory genetic variants in the development of BC by influencing TF-DNA interaction and provide critical insights to pinpoint causal genetic variants.

Abstract 2358: Uncovering novel mechanisms of resistance in AML using integrative functional genomics

Abstract Treatment options that effectively cure patients diagnosed with acute myeloid leukemia (AML) continue to represent an area of unmet need in oncology clinical care. While remission rates in AML patients can reach upwards of 90% under the current frontline therapy paradigm, nearly all patients relapse with treatment refractory disease less than 5 years after diagnosis. Relapse driven by therapy resistant cells that persist in the body after treatment (defined as minimal residual disease, or MRD) is the principal source of fatality in AML patients. Therefore, understanding how and where these leukemic cells survive treatment in vivo may help advance the rational development of highly synergistic combination therapies for the treatment of AML. Using a functional genomic approach (in vivo RNAi) and a new mouse model of AML chemoresistance (ChemoR) generated in our labs, we have identified several putative mediators of therapy resistance. Transcriptional profiling of the ChemoR model allowed us to generate a chemoresistance gene signature that we overlapped with the results of the shRNA screen to identify high-confidence genes of interest. The top genes from a ranked list of the most highly overexpressed (OE) genes in ChemoR cells and the top depleted genes from the shRNA screen in the context of therapy treatment were selected as high interest hits and subsequently tagged for individual follow-up experiments. Early validation studies suggest that some of our hits protect leukemic cells resident only at specific anatomical niches from therapy, while other proteins appear to be general chemosensitizers both in vivo and in vitro, and in both murine and human AML. Subsequent work will focus on elucidating the mechanisms by which these proteins promote resistance to frontline therapies in AML. Citation Format: Azucena Ramos, Luis R. Millan Barea, Alexandre Puissant, Nina Fenouille, Gabriela Alexe, Kimberly Stegmaier, Michael T. Hemann. Uncovering novel mechanisms of resistance in AML using integrative functional genomics [abstract]. In: Proceedings of the American Association for Cancer Research Annual Meeting 2018; 2018 Apr 14-18; Chicago, IL. Philadelphia (PA): AACR; Cancer Res 2018;78(13 Suppl):Abstract nr 2358.

A Mixed-Effects Model for Powerful Association Tests in Integrative Functional Genomics

Genome-wide association studies (GWASs) have successfully identified thousands of genetic variants for many complex diseases; however, these variants explain only a small fraction of the heritability. Recently, genetic association studies that leverage external transcriptome data have received much attention and shown promise for discovering novel variants. One such approach, PrediXcan, is to use predicted gene expression through genetic regulation. However, there are limitations in this approach. The predicted gene expression may be biased, resulting from regularized regression applied to moderately sample-sized reference studies. Further, some variants can individually influence disease risk through alternative functional mechanisms besides expression. Thus, testing only the association of predicted gene expression as proposed in PrediXcan will potentially lose power. To tackle these challenges, we consider a unified mixed effects model that formulates the association of intermediate phenotypes such as imputed gene expression through fixed effects, while allowing residual effects of individual variants to be random. We consider a set-based score testing framework, MiST (mixed effects score test), and propose two data-driven combination approaches to jointly test for the fixed and random effects. We establish the asymptotic distributions, which enable rapid calculation of p values for genome-wide analyses, and provide p values for fixed and random effects separately to enhance interpretability over GWASs. Extensive simulations demonstrate that our approaches are more powerful than existing ones. We apply our approach to a large-scale GWAS of colorectal cancer and identify two genes, POU5F1B and ATF1, which would have otherwise been missed by PrediXcan, after adjusting for all known loci.

An integrative functional genomics framework for effective identification of novel regulatory variants in genome\u2013phenome studies

BackgroundGenome–phenome studies have identified thousands of variants that are statistically associated with disease or traits; however, their functional roles are largely unclear. A comprehensive investigation of regulatory mechanisms and the gene regulatory networks between phenome-wide association study (PheWAS) and genome-wide association study (GWAS) is needed to identify novel regulatory variants contributing to risk for human diseases.MethodsIn this study, we developed an integrative functional genomics framework that maps 215,107 significant single nucleotide polymorphism (SNP) traits generated from the PheWAS Catalog and 28,870 genome-wide significant SNP traits collected from the GWAS Catalog into a global human genome regulatory map via incorporating various functional annotation data, including transcription factor (TF)-based motifs, promoters, enhancers, and expression quantitative trait loci (eQTLs) generated from four major functional genomics databases: FANTOM5, ENCODE, NIH Roadmap, and Genotype-Tissue Expression (GTEx). In addition, we performed a tissue-specific regulatory circuit analysis through the integration of the identified regulatory variants and tissue-specific gene expression profiles in 7051 samples across 32 tissues from GTEx.ResultsWe found that the disease-associated loci in both the PheWAS and GWAS Catalogs were significantly enriched with functional SNPs. The integration of functional annotations significantly improved the power of detecting novel associations in PheWAS, through which we found a number of functional associations with strong regulatory evidence in the PheWAS Catalog. Finally, we constructed tissue-specific regulatory circuits for several complex traits: mental diseases, autoimmune diseases, and cancer, via exploring tissue-specific TF-promoter/enhancer-target gene interaction networks. We uncovered several promising tissue-specific regulatory TFs or genes for Alzheimer’s disease (e.g. ZIC1 and STX1B) and asthma (e.g. CSF3 and IL1RL1).ConclusionsThis study offers powerful tools for exploring the functional consequences of variants generated from genome–phenome association studies in terms of their mechanisms on affecting multiple complex diseases and traits.

Cellular microRNA networks regulate host dependency of hepatitis C virus infection

Cellular microRNAs (miRNAs) have been shown to regulate hepatitis C virus (HCV) replication, yet a systematic interrogation of the repertoire of miRNAs impacting HCV life cycle is lacking. Here we apply integrative functional genomics strategies to elucidate global HCV–miRNA interactions. Through genome-wide miRNA mimic and hairpin inhibitor phenotypic screens, and miRNA–mRNA transcriptomics analyses, we identify three proviral and nine antiviral miRNAs that interact with HCV. These miRNAs are functionally linked to particular steps of HCV life cycle and related viral host dependencies. Further mechanistic studies demonstrate that miR-25, let-7, and miR-130 families repress essential HCV co-factors, thus restricting viral infection at multiple stages. HCV subverts the antiviral actions of these miRNAs by dampening their expression in cell culture models and HCV-infected human livers. This comprehensive HCV–miRNA interaction map provides fundamental insights into HCV-mediated pathogenesis and unveils molecular pathways linking RNA biology to viral infections.

Drawing Links from Transcriptome to Metabolites: The Evolution of Aroma in the Ripening Berry of Moscato Bianco (Vitis vinifera L.).

Monoterpenes confer typical floral notes to “Muscat” grapevine varieties and, to a lesser extent, to other aromatic non-Muscat varieties. Previous studies have led to the identification and functional characterization of some enzymes and genes in this pathway. However, the underlying genetic map is still far from being complete. For example, the specific steps of monoterpene metabolism and its regulation are largely unknown. With the aim of identifying new candidates for the missing links, we applied an integrative functional genomics approach based on the targeted metabolic and genome-wide transcript profiling of Moscato Bianco ripening berries. In particular, gas chromatography-mass spectrometry analysis of free and bound terpenoid compounds was combined with microarray analysis in the skins of berries collected at five developmental stages from pre-veraison to over-ripening. Differentially expressed metabolites and probes were identified in the pairwise comparison between time points by using the early stage as a reference. Metabolic and transcriptomic data were integrated through pairwise correlation and clustering approaches to discover genes linked with particular metabolites or groups of metabolites. These candidate transcripts were further checked for co-localization with quantitative trait loci (QTLs) affecting aromatic compounds. Our findings provide insights into the biological networks of grapevine secondary metabolism, both at the catalytic and regulatory levels. Examples include a nudix hydrolase as component of a terpene synthase-independent pathway for monoterpene biosynthesis, genes potentially involved in monoterpene metabolism (cytochrome P450 hydroxylases, epoxide hydrolases, glucosyltransferases), transport (vesicle-associated proteins, ABCG transporters, glutathione S-transferases, amino acid permeases), and transcriptional control (transcription factors of the ERF, MYB and NAC families, intermediates in light- and circadian cycle-mediated regulation with supporting evidence from the literature and additional regulatory genes with a previously unreported association to monoterpene accumulation).

Integrative Functional Genomics for Systems Genetics in GeneWeaver.org.

The abundance of existing functional genomics studies permits an integrative approach to interpreting and resolving the results of diverse systems genetics studies. However, a major challenge lies in assembling and harmonizing heterogeneous data sets across species for facile comparison to the positional candidate genes and coexpression networks that come from systems genetic studies. GeneWeaver is an online database and suite of tools at www.geneweaver.org that allows for fast aggregation and analysis of gene set-centric data. GeneWeaver contains curated experimental data together with resource-level data such as GO annotations, MP annotations, and KEGG pathways, along with persistent stores of user entered data sets. These can be entered directly into GeneWeaver or transferred from widely used resources such as GeneNetwork.org. Data are analyzed using statistical tools and advanced graph algorithms to discover new relations, prioritize candidate genes, and generate function hypotheses. Here we use GeneWeaver to find genes common to multiple gene sets, prioritize candidate genes from a quantitative trait locus, and characterize a set of differentially expressed genes. Coupling a large multispecies repository curated and empirical functional genomics data to fast computational tools allows for the rapid integrative analysis of heterogeneous data for interpreting and extrapolating systems genetics results.

EPB41, a novel hepatoma susceptibility gene dysregulated by c-Myc: an integrative functional genomics study

BackgroundGenome-wide association studies (GWAS) have provided insight into cancer genetics. However, molecular mechanisms whereby many susceptibility single nucleotide polymorphisms (SNPs) identified by GWAS promote cancer heritability and risk are unknown. New research strategies to evaluate functionality are needed to systematically study causal genetic variants. MethodsIn this study, we developed an integrative functional genomics method to identify cancer susceptibility SNPs in transcription factor binding sites across the whole genome. By integration of functional genomic data from c-Myc cistromics, 1000 Genomes and the TRANSFAC matrix databases, we identified SNPs in the c-Myc cistrome that might modulate c-Myc binding affinity in hepatocellular carcinoma. We took blood samples from patients with hepatitis B-related hepatocellular carcinoma who were admitted to Shandong Cancer Hospital, Jinan, Shandong Province. China, or Huaian No. 2 Hospital, Huaian, Jiangsu Province, China. We extracted genomic DNA from these samples and did genotype analyses. Written informed consent was obtained from all participants and the study was ethically approved by the institutional review boards. FindingsWe identified 12 SNPs in 1806 people with hepatitis B-related hepatocellular carcinoma, we genotyped these SNPs and 1708 controls. We identified a novel hepatocellular carcinoma susceptibility SNP, EPB41 rs157224 G>T (OR of T allele 1·64 [95% CI 1·32–2·02]; p<0·00001), in Chinese patients with hepatocellular carcinoma. This SNP in EPB41 predisposes individuals to hepatocellular carcinoma through modification of c-Myc-mediated transcriptional regulation; the risk allele has decreased gene expression. EPB41 is a novel hepatocellular carcinoma susceptibility gene. Consistent with this notion, a EPB41 expression is decreased tissue samples from patients with hepatocellular carcinoma, especially portal vein metastasis or intrahepatic metastasis. InterpretationOur results show that SNP EPB41 rs157224 G>T could be associated with increased susceptibility to hepatocellular carcinoma and show that insight into malignancy pathology can be made using a genome-wide approach. FundingNational Natural Science Foundation of China (81201586 and 31271382), National High-Tech Research and Development Program of China (2015AA020950), State Key Laboratory of Molecular Oncology (SKL-KF-2015-05).

Integrative Functional Genomics Implicates EPB41 Dysregulation in Hepatocellular Carcinoma Risk

Genome-wide association studies (GWASs) have provided many insights into cancer genetics. However, the molecular mechanisms of many susceptibility SNPs defined by GWASs in cancer heritability and in promoting cancer risk remain elusive. New research strategies, including functional evaluations, are warranted to systematically explore truly causal genetic variants. In this study, we developed an integrative functional genomics methodology to identify cancer susceptibility SNPs in transcription factor-binding sites across the whole genome. Employing integration of functional genomic data from c-Myc cistromics, 1000 Genomes, and the TRANSFAC matrix, we successfully annotated 12 SNPs present in the c-Myc cistrome with properties consistent with modulating c-Myc binding affinity in hepatocellular carcinoma (HCC). After genotyping these 12 SNPs in 1,806 HBV-related HCC case subjects and 1,708 control subjects, we identified a HCC susceptibility SNP, rs157224G>T, in Chinese populations (T allele: odds ratio = 1.64, 95% confidence interval = 1.32-2.02; p = 5.2× 10(-6)). This polymorphism leads to HCC predisposition through modifying c-Myc-mediated transcriptional regulation of EPB41, with the risk rs157224T allele showing significantly decreased gene expression. Based on cell proliferation, wound healing, and transwell assays as well as the mouse xenograft model, we identify EPB41 as a HCC susceptibility gene invitro and invivo. Consistent with this notion, we note that EPB41 expression is significantly decreased in HCC tissue specimens, especially in portal vein metastasis or intrahepatic metastasis, compared to normal tissues. Our results highlight the involvement of regulatory genetic variants in HCC and provide pathogenic insights of this malignancy via a genome-wide approach.

Integrative functional genomics identifies regulatory mechanisms at coronary artery disease loci.

Coronary artery disease (CAD) is the leading cause of mortality and morbidity, driven by both genetic and environmental risk factors. Meta-analyses of genome-wide association studies have identified >150 loci associated with CAD and myocardial infarction susceptibility in humans. A majority of these variants reside in non-coding regions and are co-inherited with hundreds of candidate regulatory variants, presenting a challenge to elucidate their functions. Herein, we use integrative genomic, epigenomic and transcriptomic profiling of perturbed human coronary artery smooth muscle cells and tissues to begin to identify causal regulatory variation and mechanisms responsible for CAD associations. Using these genome-wide maps, we prioritize 64 candidate variants and perform allele-specific binding and expression analyses at seven top candidate loci: 9p21.3, SMAD3, PDGFD, IL6R, BMP1, CCDC97/TGFB1 and LMOD1. We validate our findings in expression quantitative trait loci cohorts, which together reveal new links between CAD associations and regulatory function in the appropriate disease context.

Hepatitis C virus depends on E-cadherin as an entry factor and regulates its expression in epithelial-to-mesenchymal transition

Hepatitis C virus (HCV) enters the host cell through interactions with a cascade of cellular factors. Although significant progress has been made in understanding HCV entry, the precise mechanisms by which HCV exploits the receptor complex and host machinery to enter the cell remain unclear. This intricate process of viral entry likely depends on additional yet-to-be-defined cellular molecules. Recently, by applying integrative functional genomics approaches, we identified and interrogated distinct sets of host dependencies in the complete HCV life cycle. Viral entry assays using HCV pseudoparticles (HCVpps) of various genotypes uncovered multiple previously unappreciated host factors, including E-cadherin, that mediate HCV entry. E-cadherin silencing significantly inhibited HCV infection in Huh7.5.1 cells, HepG2/miR122/CD81 cells, and primary human hepatocytes at a postbinding entry step. Knockdown of E-cadherin, however, had no effect on HCV RNA replication or internal ribosomal entry site (IRES)-mediated translation. In addition, an E-cadherin monoclonal antibody effectively blocked HCV entry and infection in hepatocytes. Mechanistic studies demonstrated that E-cadherin is closely associated with claudin-1 (CLDN1) and occludin (OCLN) on the cell membrane. Depletion of E-cadherin drastically diminished the cell-surface distribution of these two tight junction proteins in various hepatic cell lines, indicating that E-cadherin plays an important regulatory role in CLDN1/OCLN localization on the cell surface. Furthermore, loss of E-cadherin expression in hepatocytes is associated with HCV-induced epithelial-to-mesenchymal transition (EMT), providing an important link between HCV infection and liver cancer. Our data indicate that a dynamic interplay among E-cadherin, tight junctions, and EMT exists and mediates an important function in HCV entry.