Functional annotation of colon cancer risk SNPs

An oncogenic enhancer promotes melanoma progression via regulating ETV4 expression

BACKGROUND

Enhancers are important gene regulatory elements that promote the expression of critical genes in development and disease. Aberrant enhancer can modulate cancer risk and activate oncogenes that lead to the occurrence of various cancers. However, the underlying mechanism of most enhancers in cancer remains unclear. Here, we aim to explore the function and mechanism of a crucial enhancer in melanoma.

METHODS

Multi-omics data were applied to identify an enhancer (enh17) involved in melanoma progression. To evaluate the function of enh17, CRISPR/Cas9 technology were applied to knockout enh17 in melanoma cell line A375. RNA-seq, ChIP-seq and Hi-C data analysis integrated with luciferase reporter assay were performed to identify the potential target gene of enh17. Functional experiments were conducted to further validate the function of the target gene ETV4. Multi-omics data integrated with CUT&Tag sequencing were performed to validate the binding profile of the inferred transcription factor STAT3.

RESULTS

An enhancer, named enh17 here, was found to be aberrantly activated and involved in melanoma progression. CRISPR/Cas9-mediated deletion of enh17 inhibited cell proliferation, migration, and tumor growth of melanoma both in vitro and in vivo. Mechanistically, we identified ETV4 as a target gene regulated by enh17, and functional experiments further support ETV4 as a target gene that is involved in cancer-associated phenotypes. In addition, STAT3 acts as a transcription factor binding with enh17 to regulate the transcription of ETV4.

CONCLUSIONS

Our findings revealed that enh17 plays an oncogenic role and promotes tumor progression in melanoma, and its transcriptional regulatory mechanisms were fully elucidated, which may open a promising window for melanoma prevention and treatment.

Genome-wide association studies and Mendelian randomization analyses provide insights into the causes of early-onset colorectal cancer

BACKGROUND

The incidence of early-onset colorectal cancer (EOCRC; diagnosed <50 years of age) is rising globally; however, the causes underlying this trend are largely unknown. CRC has strong genetic and environmental determinants, yet common genetic variants and causal modifiable risk factors underlying EOCRC are unknown. We conducted the first EOCRC-specific genome-wide association study (GWAS) and Mendelian randomization (MR) analyses to explore germline genetic and causal modifiable risk factors associated with EOCRC.

PATIENTS AND METHODS

We conducted a GWAS meta-analysis of 6176 EOCRC cases and 65 829 controls from the Genetics and Epidemiology of Colorectal Cancer Consortium (GECCO), the Colorectal Transdisciplinary Study (CORECT), the Colon Cancer Family Registry (CCFR), and the UK Biobank. We then used the EOCRC GWAS to investigate 28 modifiable risk factors using two-sample MR.

RESULTS

We found two novel risk loci for EOCRC at 1p34.1 and 4p15.33, which were not previously associated with CRC risk. We identified a deleterious coding variant (rs36053993, G396D) at polyposis-associated DNA repair gene MUTYH (odds ratio 1.80, 95% confidence interval 1.47-2.22) but show that most of the common genetic susceptibility was from noncoding signals enriched in epigenetic markers present in gastrointestinal tract cells. We identified new EOCRC-susceptibility genes, and in addition to pathways such as transforming growth factor (TGF) β, suppressor of Mothers Against Decapentaplegic (SMAD), bone morphogenetic protein (BMP) and phosphatidylinositol kinase (PI3K) signaling, our study highlights a role for insulin signaling and immune/infection-related pathways in EOCRC. In our MR analyses, we found novel evidence of probable causal associations for higher levels of body size and metabolic factors-such as body fat percentage, waist circumference, waist-to-hip ratio, basal metabolic rate, and fasting insulin-higher alcohol drinking, and lower education attainment with increased EOCRC risk.

CONCLUSIONS

Our novel findings indicate inherited susceptibility to EOCRC and suggest modifiable lifestyle and metabolic targets that could also be used to risk-stratify individuals for personalized screening strategies or other interventions.

Integrated meta-analysis of colorectal cancer public proteomic datasets for biomarker discovery and validation

The cancer biomarker field has been an object of thorough investigation in the last decades. Despite this, colorectal cancer (CRC) heterogeneity makes it challenging to identify and validate effective prognostic biomarkers for patient classification according to outcome and treatment response. Although a massive amount of proteomics data has been deposited in public data repositories, this rich source of information is vastly underused. Here, we attempted to reuse public proteomics datasets with two main objectives: i) to generate hypotheses (detection of biomarkers) for their posterior/downstream validation, and (ii) to validate, using an orthogonal approach, a previously described biomarker panel. Twelve CRC public proteomics datasets (mostly from the PRIDE database) were re-analysed and integrated to create a landscape of protein expression. Samples from both solid and liquid biopsies were included in the reanalysis. Integrating this data with survival annotation data, we have validated in silico a six-gene signature for CRC classification at the protein level, and identified five new blood-detectable biomarkers (CD14, PPIA, MRC2, PRDX1, and TXNDC5) associated with CRC prognosis. The prognostic value of these blood-derived proteins was confirmed using additional public datasets, supporting their potential clinical value. As a conclusion, this proof-of-the-concept study demonstrates the value of re-using public proteomics datasets as the basis to create a useful resource for biomarker discovery and validation. The protein expression data has been made available in the public resource Expression Atlas.

The Clinical Significance of MicroRNAs in Colorectal Cancer Signaling Pathways: A Review

Colorectal carcinoma (colon and rectum) is currently considered among the most prevalent malignancies of Western societies. The pathogenesis and etiological mechanisms underlying colorectal cancer (CRC) development remain complex and heterogeneous. The homeostasis and function of normal human intestinal cells is highly regulated by microRNAs. Therefore, it is not surprising that mutations and inactivation of these molecules appear to be linked with progression of colorectal tumors. Recent studies have reported significant alterations of microRNA expression in adenomas and CRCs compared with adjacent normal tissues. This observed deviation has been proposed to correlate with the progression and survival of disease as well as with choice of optimal treatment and drug resistance. MicroRNAs can adopt either oncogenic or tumor-suppressive roles during regulation of pathways that drive carcinogenesis. Typically, oncogenic microRNAs termed oncomirs, target and silence endogenous tumor-suppressor genes. On the other hand, tumor-suppressive microRNAs are critical in downregulating genes associated with cell growth and malignant capabilities. By extensively evaluating robust studies, we have emphasized and distinguished a discrete set of microRNAs that can modulate tumor progression by silencing specific driver genes crucial in signaling pathways including Wnt/b-catenin, epidermal growth factor receptor, P53, mismatch repair DNA repair, and transforming-growth factor beta.

Emerging insights into enhancer biology and function

Cell type-specific gene expression is coordinated by DNA-encoded enhancers and the transcription factors (TFs) that bind to them in a sequence-specific manner. As such, these enhancers and TFs are critical mediators of normal development and altered enhancer or TF function is associated with the development of diseases such as cancer. While initially defined by their ability to activate gene transcription in reporter assays, putative enhancer elements are now frequently defined by their unique chromatin features including DNase hypersensitivity and transposase accessibility, bidirectional enhancer RNA (eRNA) transcription, CpG hypomethylation, high H3K27ac and H3K4me1, sequence-specific transcription factor binding, and co-factor recruitment. Identification of these chromatin features through sequencing-based assays has revolutionized our ability to identify enhancer elements on a genome-wide scale, and genome-wide functional assays are now capitalizing on this information to greatly expand our understanding of how enhancers function to provide spatiotemporal coordination of gene expression programs. Here, we highlight recent technological advances that are providing new insights into the molecular mechanisms by which these critical cis-regulatory elements function in gene control. We pay particular attention to advances in our understanding of enhancer transcription, enhancer-promoter syntax, 3D organization and biomolecular condensates, transcription factor and co-factor dependencies, and the development of genome-wide functional enhancer screens.

Genetic polymorphisms of CYP24A1 gene and cancer susceptibility: a meta-analysis including 40640 subjects

BACKGROUND

Whether cytochrome P450 24A1 (CYP24A1) polymorphism is associated with cancer susceptibility, the individual study results are still controversial. Therefore, we performed a comprehensive study to identify the association of CYP24A1 polymorphisms (rs4809960, rs6068816, rs2296241, rs4809957, rs2762939) with cancer susceptibility.

METHODS

Electronic databases including Cochrane Library, PubMed, and Embase were systematically retrieved for relevant publications. Fixed or random-effect model was selected to calculate odds ratios (ORs) with their 95% confidence intervals (95%CI).

RESULTS

Eighteen published articles were identified. The results indicated that rs4809960 polymorphism was associated with a decreased cancer risk in Caucasian (TT vs. TC+CC: P=0.035; C vs. T: P=0.016) and Asian population (CC vs. TC+TT: OR P=0.044; TT vs. TC+CC: P=0.021; CC vs. TT: P=0.020; C vs. T: P=0.008) and breast cancer risk (TT vs. TC+CC: P = 0.007; TC vs. TT: P=0.004; C vs. T: P=0.033). A significant association was found between rs2296241 polymorphism and esophageal squamous cell carcinoma risk (AA vs. GG+AG: P = 0.023) and prostate cancer susceptibility (A vs. G: P=0.022). Furthermore, rs4809957 polymorphism was associated with prostate cancer susceptibility in Caucasian (GG vs. GA+AA: P=0.029; GA vs. GG: P=0.022) and breast cancer susceptibility (AA vs. GG+GA: P=0.012; AA vs. GG, P=0.010; A vs. G: P=0.024). Additionally, rs6068816 polymorphism significantly decreased the lung cancer (CC vs. CT+TT: P = 0.016; TT vs. CC: P = 0.044; CT vs. CC: P = 0.036; T vs. C: P = 0.016) and breast cancer risk (TT vs. CC+CT: P = 0.043; TT vs. CC: P = 0.039). No association was found for rs2762939 polymorphism with overall cancer risk. However, for rs2296241, rs4809957, and rs6068816 polymorphisms, there were no significant differences after the Bonferroni correction.

CONCLUSION

The meta-analysis suggested that rs4809960 was associated with cancer risk and might be a genetic marker for predicting cancer risk. More large-scale and large-sample studies are necessary to further confirm these results.

Exome-based new allele-specific PCR markers and transferability for sodicity tolerance in bread wheat (Triticum aestivum L.)

Targeted exome-based genotype by sequencing (t-GBS), a sequencing technology that tags SNPs and haplotypes in gene-rich regions was used in previous genome-wide association studies (GWAS) for sodicity tolerance in bread wheat. Thirty-nine novel SNPs including 18 haplotypes for yield and yield-components were identified. The present study aimed at developing SNP-derived markers by precisely locating new SNPs on ~180 bp allelic sequence of t-GBS, marker validation, and SNP functional characterization based on its exonic location. We identified unknown locations of significant SNPs/haplotypes by aligning allelic sequences on to IWGSC RefSeqv1.0 on respective chromosomes. Eighteen out of the target 39 SNP locations fulfilled the criteria for producing PCR markers, among which only eight produced polymorphic signals. These eight markers associated with yield, plants m-2, heads m-2, and harvest index, including a pleiotropic marker for yield, harvest index, and grains/head were validated for its amplification efficiency and phenotypic effects in focused identification germplasm strategy (FIGS) wheat set and a doubled haploid (DH) population (Scepter/IG107116). The phenotypic variation explained by these markers are in the range of 4.1-37.6 in the FIGS population. High throughput PCR-based genotyping using new markers and association with phenotypes in FIGS wheat set and DH population validated the effect of functional SNP on closely associated genes-calcineurin B-like- and dirigent protein, basic helix-loop-helix (BHLH-), plant homeodomain (PHD-) and helix-turn-helix myeloblastosis (HTH myb) type -transcription factor. Further, genome-wide SNP annotation using SnpEff tool confirmed that these SNPs are in gene regulatory regions (upstream, 3'-UTR, and intron) modifying gene expression and protein-coding. This integrated approach of marker design for t-GBS alleles, SNP functional annotation, and high-throughput genotyping of functional SNP offers translation solutions across crops and complex traits in crop improvement programs.

The EN-TEx resource of multi-tissue personal epigenomes & variant-impact models

Understanding how genetic variants impact molecular phenotypes is a key goal of functional genomics, currently hindered by reliance on a single haploid reference genome. Here, we present the EN-TEx resource of 1,635 open-access datasets from four donors (∼30 tissues × ∼15 assays). The datasets are mapped to matched, diploid genomes with long-read phasing and structural variants, instantiating a catalog of >1 million allele-specific loci. These loci exhibit coordinated activity along haplotypes and are less conserved than corresponding, non-allele-specific ones. Surprisingly, a deep-learning transformer model can predict the allele-specific activity based only on local nucleotide-sequence context, highlighting the importance of transcription-factor-binding motifs particularly sensitive to variants. Furthermore, combining EN-TEx with existing genome annotations reveals strong associations between allele-specific and GWAS loci. It also enables models for transferring known eQTLs to difficult-to-profile tissues (e.g., from skin to heart). Overall, EN-TEx provides rich data and generalizable models for more accurate personal functional genomics.

Application of CRISPR-Cas9 gene editing technology in basic research, diagnosis and treatment of colon cancer

Colon cancer is the fourth leading cause of cancer death worldwide, and its progression is accompanied by a complex array of genetic variations. CRISPR/Cas9 can identify new drug-resistant or sensitive mutations in colon cancer, and can use gene editing technology to develop new therapeutic targets and provide personalized treatments, thereby significantly improving the treatment of colon cancer patients. CRISPR/Cas9 systems are driving advances in biotechnology. RNA-directed Cas enzymes have accelerated the pace of basic research and led to clinical breakthroughs. This article reviews the rapid development of CRISPR/Cas in colon cancer, from gene editing to transcription regulation, gene knockout, genome-wide CRISPR tools, therapeutic targets, stem cell genomics, immunotherapy, metabolism-related genes and inflammatory bowel disease. In addition, the limitations and future development of CRISPR/Cas9 in colon cancer studies are reviewed. In conclusion, this article reviews the application of CRISPR-Cas9 gene editing technology in basic research, diagnosis and treatment of colon cancer.

Role of Non-Coding RNAs in Colorectal Cancer: Focus on Long Non-Coding RNAs

Colorectal cancer is one of the most common causes of cancer-related deaths worldwide. Despite the advances in the knowledge of pathogenetic molecular mechanisms and the implementation of more effective drug treatments in recent years, the overall survival rate of patients remains unsatisfactory. The high death rate is mainly due to metastasis of cancer in about half of the cancer patients and the emergence of drug-resistant populations of cancer cells. Improved understanding of cancer molecular biology has highlighted the role of non-coding RNAs (ncRNAs) in colorectal cancer development and evolution. ncRNAs regulate gene expression through various mechanisms, including epigenetic modifications and interactions of long non-coding RNAs (lncRNAs) with both microRNAs (miRNAs) and proteins, and through the action of lncRNAs as miRNA precursors or pseudogenes. LncRNAs can also be detected in the blood and circulating ncRNAs have become a new source of non-invasive cancer biomarkers for the diagnosis and prognosis of colorectal cancer, as well as for predicting the response to drug therapy. In this review, we focus on the role of lncRNAs in colorectal cancer development, progression, and chemoresistance, and as possible therapeutic targets.

An enhancer variant at 16q22.1 predisposes to hepatocellular carcinoma via regulating PRMT7 expression

Most cancer causal variants are found in gene regulatory elements, e.g., enhancers. However, enhancer variants predisposing to hepatocellular carcinoma (HCC) remain unreported. Here we conduct a genome-wide survey of HCC-susceptible enhancer variants through a three-stage association study in 11,958 individuals and identify rs73613962 (T > G) within the intronic region of PRMT7 at 16q22.1 as a susceptibility locus of HCC (OR = 1.41, P = 6.02 × 10^-10). An enhancer dual-luciferase assay indicates that the rs73613962-harboring region has allele-specific enhancer activity. CRISPR-Cas9/dCas9 experiments further support the enhancer activity of this region to regulate PRMT7 expression. Mechanistically, transcription factor HNF4A binds to this enhancer region, with preference to the risk allele G, to promote PRMT7 expression. PRMT7 upregulation contributes to in vitro, in vivo, and clinical HCC-associated phenotypes, possibly by affecting the p53 signaling pathway. This concept of HCC pathogenesis may open a promising window for HCC prevention/treatment.

Integrative epigenomic and high-throughput functional enhancer profiling reveals determinants of enhancer heterogeneity in gastric cancer

BACKGROUND

Enhancers are distal cis-regulatory elements required for cell-specific gene expression and cell fate determination. In cancer, enhancer variation has been proposed as a major cause of inter-patient heterogeneity-however, most predicted enhancer regions remain to be functionally tested.

METHODS

We analyzed 132 epigenomic histone modification profiles of 18 primary gastric cancer (GC) samples, 18 normal gastric tissues, and 28 GC cell lines using Nano-ChIP-seq technology. We applied Capture-based Self-Transcribing Active Regulatory Region sequencing (CapSTARR-seq) to assess functional enhancer activity. An Activity-by-contact (ABC) model was employed to explore the effects of histone acetylation and CapSTARR-seq levels on enhancer-promoter interactions.

RESULTS

We report a comprehensive catalog of 75,730 recurrent predicted enhancers, the majority of which are GC-associated in vivo (> 50,000) and associated with lower somatic mutation rates inferred by whole-genome sequencing. Applying CapSTARR-seq to the enhancer catalog, we observed significant correlations between CapSTARR-seq functional activity and H3K27ac/H3K4me1 levels. Super-enhancer regions exhibited increased CapSTARR-seq signals compared to regular enhancers, even when decoupled from native chromatin contexture. We show that combining histone modification and CapSTARR-seq functional enhancer data improves the prediction of enhancer-promoter interactions and pinpointing of germline single nucleotide polymorphisms (SNPs), somatic copy number alterations (SCNAs), and trans-acting TFs involved in GC expression. We identified cancer-relevant genes (ING1, ARL4C) whose expression between patients is influenced by enhancer differences in genomic copy number and germline SNPs, and HNF4α as a master trans-acting factor associated with GC enhancer heterogeneity.

CONCLUSIONS

Our results indicate that combining histone modification and functional assay data may provide a more accurate metric to assess enhancer activity than either platform individually, providing insights into the relative contribution of genetic (cis) and regulatory (trans) mechanisms to GC enhancer functional heterogeneity.

Novel insights into the molecular mechanisms underlying risk of colorectal cancer from smoking and red/processed meat carcinogens by modeling exposure in normal colon organoids

Tobacco smoke and red/processed meats are well-known risk factors for colorectal cancer (CRC). Most research has focused on studies of normal colon biopsies in epidemiologic studies or treatment of CRC cell lines in vitro. These studies are often constrained by challenges with accuracy of self-report data or, in the case of CRC cell lines, small sample sizes and lack of relationship to normal tissue at risk. In an attempt to address some of these limitations, we performed a 24-hour treatment of a representative carcinogens cocktail in 37 independent organoid lines derived from normal colon biopsies. Machine learning algorithms were applied to bulk RNA-sequencing and revealed cellular composition changes in colon organoids. We identified 738 differentially expressed genes in response to carcinogens exposure. Network analysis identified significantly different modules of co-expression, that included genes related to MSI-H tumor biology, and genes previously implicated in CRC through genome-wide association studies. Our study helps to better define the molecular effects of representative carcinogens from smoking and red/processed meat in normal colon epithelial cells and in the etiology of the MSI-H subtype of CRC, and suggests an overlap between molecular mechanisms involved in inherited and environmental CRC risk.

CD8+ T-cell Immune Surveillance against a Tumor Antigen Encoded by the Oncogenic Long Noncoding RNA PVT1

CD8⁺ T cells recognize peptides displayed by HLA class I molecules on cell surfaces, monitoring pathologic conditions such as cancer. Advances in proteogenomic analysis of HLA ligandomes have demonstrated that cells present a subset of cryptic peptides derived from noncoding regions of the genome; however, the roles of cryptic HLA ligands in tumor immunity remain unknown. In the current study, we comprehensively and quantitatively investigated the HLA class I ligandome of a set of human colorectal cancer and matched normal tissues, showing that cryptic translation products accounted for approximately 5% of the HLA class I ligandome. We also found that a peptide encoded by the long noncoding RNA (lncRNA) PVT1 was predominantly enriched in multiple colorectal cancer tissues. The PVT1 gene is located downstream of the MYC gene in the genome and is aberrantly overexpressed across a variety of cancers, reflecting its oncogenic property. The PVT1 peptide was recognized by patient CD8⁺ tumor-infiltrating lymphocytes, as well as peripheral blood mononuclear cells, suggesting the presence of patient immune surveillance. Our findings show that peptides can be translated from lncRNAs and presented by HLA class I and that cancer patient T cells are capable of sensing aberrations in noncoding regions of the genome.

Epigenetic alteration contributes to the transcriptional reprogramming in T-cell prolymphocytic leukemia

T cell prolymphocytic leukemia (T-PLL) is a rare disease with aggressive clinical course. Cytogenetic analysis, whole-exome and whole-genome sequencing have identified primary structural alterations in T-PLL, including inversion, translocation and copy number variation. Recurrent somatic mutations were also identified in genes encoding chromatin regulators and those in the JAK-STAT signaling pathway. Epigenetic alterations are the hallmark of many cancers. However, genome-wide epigenomic profiles have not been reported in T-PLL, limiting the mechanistic study of its carcinogenesis. We hypothesize epigenetic mechanisms also play a key role in T-PLL pathogenesis. To systematically test this hypothesis, we generated genome-wide maps of regulatory regions using H3K4me3 and H3K27ac ChIP-seq, as well as RNA-seq data in both T-PLL patients and healthy individuals. We found that genes down-regulated in T-PLL are mainly associated with defense response, immune system or adaptive immune response, while up-regulated genes are enriched in developmental process, as well as WNT signaling pathway with crucial roles in cell fate decision. In particular, our analysis revealed a global alteration of regulatory landscape in T-PLL, with differential peaks highly enriched for binding motifs of immune related transcription factors, supporting the epigenetic regulation of oncogenes and genes involved in DNA damage response and T-cell activation. Together, our work reveals a causal role of epigenetic dysregulation in T-PLL.

Validation of prostate cancer risk variants rs10993994 and rs7098889 by CRISPR/Cas9 mediated genome editing

GWAS have identified numerous SNPs associated with prostate cancer risk. One such SNP is rs10993994. It is located in the β-microseminoprotein (MSMB) promoter region, mediates MSMB prostate secretion levels, and is linked to mRNA expression changes in both MSMB and the adjacent gene NCOA4. In addition, our previous work showed a second SNP, rs7098889, is in positive linkage disequilibrium with rs10993994 and associated with MSMB expression independent of rs10993994. Here, we generate a series of clones with single alleles removed by double guide RNA (gRNA) mediated CRISPR/Cas9 deletions, through which we demonstrate that each of these SNPs independently and greatly alters MSMB expression in an allele-specific manner. We further show that these SNPs have no substantial effect on the expression of NCOA4. These data demonstrate that a single SNP can have a large effect on gene expression and illustrate the importance of functional validation studies to deconvolute observed correlations. The method we have developed is generally applicable to test any SNP for which a relevant heterozygous cell line is available. AUTHOR SUMMARY: In pursuing the underlying biological mechanism of prostate cancer pathogenesis, scientists utilized the existence of common single nucleotide polymorphisms (SNPs) in the human genome as genetic markers to perform large scale genome wide association studies (GWAS) and have so far identified more than a hundred prostate cancer risk variants. Such variants provide an unbiased and systematic new venue to study the disease mechanism, and the next big challenge is to translate these genetic associations to the causal role of altered gene function in oncogenesis. The majority of these variants are waiting to be studied and lots of them may act in oncogenesis through gene expression regulation. To prove the concept, we took rs10993994 and its linked rs7098889 as an example and engineered single cell clones by allelic-specific CRISPR/Cas9 deletion to separate the effect of each allele. We observed that a single nucleotide difference would lead to surprisingly high level of MSMB gene expression change in a gene specific and cell-type specific manner. Our study strongly supports the notion that differential level of gene expression caused by risk variants and their associated genetic locus play a major role in oncogenesis and also highlights the importance of studying the function of MSMB encoded β-MSP in prostate cancer pathogenesis.

Computational Prediction of Probable Single Nucleotide Polymorphism-Cancer Relationships

Genetic variations such as single nucleotide polymorphisms (SNPs) can cause susceptibility to cancer. Although thousands of genetic variants have been identified to be associated with different cancers, the molecular mechanisms of cancer remain unknown. There is not a particular dataset of relationships between cancer and SNPs, as a bipartite network, for computational analysis and prediction. Link prediction as a computational graph analysis method can help us to gain new insight into the network. In this article, after creating a network between cancer and SNPs using SNPedia and Cancer Research UK databases, we evaluated the computational link prediction methods to foresee new SNP-Cancer relationships. Results show that among the popular scoring methods based on network topology, for relation prediction, the preferential attachment (PA) algorithm is the most robust method according to computational and experimental evidence, and some of its computational predictions are corroborated in recent publications. According to the PA predictions, rs1801394-Non-small cell lung cancer, rs4880-Non-small cell lung cancer, and rs1805794-Colorectal cancer are some of the best probable SNP-Cancer associations that have not yet been mentioned in any published article, and they are the most probable candidates for additional laboratory and validation studies. Also, it is feasible to improve the predicting algorithms to produce new predictions in the future.

PWAS: proteome-wide association study-linking genes and phenotypes by functional variation in proteins

We introduce Proteome-Wide Association Study (PWAS), a new method for detecting gene-phenotype associations mediated by protein function alterations. PWAS aggregates the signal of all variants jointly affecting a protein-coding gene and assesses their overall impact on the protein's function using machine learning and probabilistic models. Subsequently, it tests whether the gene exhibits functional variability between individuals that correlates with the phenotype of interest. PWAS can capture complex modes of heritability, including recessive inheritance. A comparison with GWAS and other existing methods proves its capacity to recover causal protein-coding genes and highlight new associations. PWAS is available as a command-line tool.

Integrative Functional Annotation of 52 Genetic Loci Influencing Myocardial Mass Identifies Candidate Regulatory Variants and Target Genes

BACKGROUND

Regulatory elements may be involved in the mechanisms by which 52 loci influence myocardial mass, reflected by abnormal amplitude and duration of the QRS complex on the ECG. Functional annotation thus far did not take into account how these elements are affected in disease context.

METHODS

We generated maps of regulatory elements on hypertrophic cardiomyopathy patients (ChIP-seq N=14 and RNA-seq N=11) and nondiseased hearts (ChIP-seq N=4 and RNA-seq N=11). We tested enrichment of QRS-associated loci on elements differentially acetylated and directly regulating differentially expressed genes between hypertrophic cardiomyopathy patients and controls. We further performed functional annotation on QRS-associated loci using these maps of differentially active regulatory elements.

RESULTS

Regions differentially affected in disease showed a stronger enrichment ( P=8.6×10^-5) for QRS-associated variants than those not showing differential activity ( P=0.01). Promoters of genes differentially regulated between hypertrophic cardiomyopathy patients and controls showed more enrichment ( P=0.001) than differentially acetylated enhancers ( P=0.8) and super-enhancers ( P=0.025). We also identified 74 potential causal variants overlapping these differential regulatory elements. Eighteen of the genes mapped confirmed previous findings, now also pinpointing the potentially affected regulatory elements and candidate causal variants. Fourteen new genes were also mapped.

CONCLUSIONS

Our results suggest differentially active regulatory elements between hypertrophic cardiomyopathy patients and controls can offer more insights into the mechanisms of QRS-associated loci than elements not affected by disease.

A functional variant in TNXB promoter associates with the risk of esophageal squamous-cell carcinoma

Our previous study identified a tag single-nucleotide polymorphism (SNP) rs204900 in TNXB associated with risk of esophageal squamous-cell carcinoma (ESCC) in the Chinese population. However, the functional role of TNXB and causal variants had not been interrogated in that study. In the present study, we explored the effects of TNXB expression in the development of ESCC and searched for functional variants in this gene. We found TNXB was downregulated in ESCC tumors. Using small interfering RNAs and CRISPR-Cas9 methods, we identified that both knockdown and knockout of TNXB significantly promoted ESCC cell growth in vitro, suggesting a tumor suppressor role of this gene in ESCC. Through further fine-mapping analysis, we identified that a noncoding variant in the promoter of TNXB, rs411337, predisposed to ESCC risk (odds ratio = 1.36, 95% confidence interval: 1.22-1.51, P = 9.10 × 10^-9 ). These findings revealed the functional mechanism of TNXB in the development of ESCC and may contribute to the prevention and treatment of this disease in the future.

PVT1 Long Non-coding RNA in Gastrointestinal Cancer

Whole genome and transcriptome sequencing technologies have led to the identification of many long non-coding RNAs (lncRNAs) and stimulated the research of their role in health and disease. LncRNAs participate in the regulation of critical signaling pathways including cell growth, motility, apoptosis, and differentiation; and their expression has been found dysregulated in human tumors. Thus, lncRNAs have emerged as new players in the initiation, maintenance and progression of tumorigenesis. PVT1 (plasmacytoma variant translocation 1) lncRNA is located on chromosomal 8q24.21, a large locus frequently amplified in human cancers and predictive of increased cancer risk in genome-wide association studies (GWAS). Combined, colorectal and gastric adenocarcinomas are the most frequent tumor malignancies and also the leading cause of cancer-related deaths worldwide. PVT1 expression is elevated in gastrointestinal tumors and correlates with poor patient prognosis. In this review, we discuss the mechanisms of action underlying PVT1 oncogenic role in colorectal and gastric cancer such as MYC upregulation, miRNA production, competitive endogenous RNA (ceRNA) function, protein stabilization, and epigenetic regulation. We also illustrate the potential role of PVT1 as prognostic biomarker and its relationship with resistance to current chemotherapeutic treatments.

Interactions between the gut microbiome and host gene regulation in cystic fibrosis

BACKGROUND

Cystic fibrosis is the most common autosomal recessive genetic disease in Caucasians. It is caused by mutations in the CFTR gene, leading to poor hydration of mucus and impairment of the respiratory, digestive, and reproductive organ functions. Advancements in medical care have led to markedly increased longevity of patients with cystic fibrosis, but new complications have emerged, such as early onset of colorectal cancer. Although the pathogenesis of colorectal cancer in cystic fibrosis remains unclear, altered host-microbe interactions might play a critical role. To investigate this, we characterized changes in the microbiome and host gene expression in the colonic mucosa of cystic fibrosis patients relative to healthy controls, and identified host gene-microbiome interactions in the colon of cystic fibrosis patients.

METHODS

We performed RNA-seq on colonic mucosa samples from cystic fibrosis patients and healthy controls to determine differentially expressed host genes. We also performed 16S rRNA sequencing to characterize the colonic mucosal microbiome and identify gut microbes that are differentially abundant between patients and healthy controls. Lastly, we modeled associations between relative abundances of specific bacterial taxa in the gut mucosa and host gene expression.

RESULTS

We find that 1543 genes, including CFTR, show differential expression in the colon of cystic fibrosis patients compared to healthy controls. These genes are enriched with functions related to gastrointestinal and colorectal cancer, such as metastasis of colorectal cancer, tumor suppression, p53, and mTOR signaling pathways. In addition, patients with cystic fibrosis show decreased gut microbial diversity, decreased abundance of butyrate producing bacteria, such as Ruminococcaceae and Butyricimonas, and increased abundance of other taxa, such as Actinobacteria and Clostridium. An integrative analysis identified colorectal cancer-related genes, including LCN2 and DUOX2, for which gene expression is correlated with the abundance of colorectal cancer-associated bacteria, such as Ruminococcaceae and Veillonella.

CONCLUSIONS

In addition to characterizing host gene expression and mucosal microbiome in cystic fibrosis patients, our study explored the potential role of host-microbe interactions in the etiology of colorectal cancer in cystic fibrosis. Our results provide biomarkers that may potentially serve as targets for stratifying risk of colorectal cancer in patients with cystic fibrosis.

A functional variant in CHK1 contributes to increased risk of nasopharyngeal carcinoma in a Han Chinese population

DNA damage checkpoints act as a supervisor by preventing the course of cell cycle upon DNA damage and keeping the steadiness of genome. Checkpoint kinase 1 (CHK1) cannot be ignore in the etiology of numerous human cancers including nasopharyngeal cancer (NPC). To discuss genetic polymorphisms of CHK1 rs492510 in the occurrence of NPC was our objective. Rs492510 polymorphism of CHK1 was genotyped in 684 patients with NPC and 823 cancer-free controls. We utilize logistic regression models to appraise the correlation of rs492510 and susceptibility of NPC. Comparative expression level about CHK1 in nasopharyngeal carcinoma tissues were determined by real-time polymerase chain reaction. And we made use of Dual-Luciferase Reporter Assay to assess the transcriptional ability of CHK1 with different rs492510 allele. Adjusting multivariate logistic regression based on age, sex, body mass index, smoking, and drinking status showed that CHK1 rs492510 GA + GG genotype carriers presented prominent higher risk in NPC (odds ratio = 1.376, 95% confidence interval: 1.087-1.742; P = .008). As a consequence, we revealed that CHK1 relative expression levels in NPC tissues was higher than rhinitis tissues. Besides, the expressions of CHK1 in rs492510 GA genotype carriers were higher compared with people in AA genotype. The G allele of rs492510 generated remarkable higher transcription activity of CHK1 vs A allele by luciferase reporter assay. Our study considered that single nucleotide polymorphism rs492510 could increase transcription activity of CHK1 with the functionality, contributing to the susceptibility of NPC.

Analyzing Metabolomics Data for Environmental Health and Exposome Research

The exposome is the cumulative measure of environmental influences and associated biological responses across the life span, with critical relevance for understanding how exposures can impact human health. Metabolomics analysis of biological samples offers unique advantages for examining the exposome. Simultaneous analysis of external exposures, biological responses, and host susceptibility at a systems level can help establish links between external exposures and health outcomes. As metabolomics technologies continue to evolve for the study of the exposome, metabolomics ultimately will help provide valuable insights for exposure risk assessment, and disease prevention and management. Here, we discuss recent advances in metabolomics, and describe data processing protocols that can enable analysis of the exposome. This chapter focuses on using liquid chromatography-mass spectrometry (LC-MS)-based untargeted metabolomics for analysis of the exposome, including (1) preprocessing of untargeted metabolomics data, (2) identification of exposure chemicals and their metabolites, and (3) methods to establish associations between exposures and diseases.

Genetic association of polymorphisms at the intergenic region between PRDM1 and ATG5 with hepatitis B virus infection in Han Chinese patients

Chronic hepatitis B virus (HBV) infection is related to chronic hepatitis, cirrhosis, and hepatocellular carcinoma (HCC), and the interplay between the virus and host immune response leads to different outcomes of the infection. PR domain zinc finger protein 1 (PRDM1) and autophagy-related protein 5 (ATG5) are involved in immune response and HBV infection. An intergenic region between PRDM1 and ATG5 (PRDM1-ATG5 region) has been identified, and single-nucleotide polymorphisms (SNPs) in this region were shown to be involved in immune regulation. This study investigated the functionally relevant rs548234, rs6937876, and rs6568431 polymorphisms at the PRDM1-ATG5 region in a Han Chinese population (403 patients with chronic HBV infection [171 chronic hepatitis, 119 cirrhosis, and 113 HCC], 70 infection resolvers, and 196 healthy controls). The frequencies of the rs6568431 allele A in HBV patients (P = .005) and genotype CA in infection resolvers (P = .005) were significantly higher than in healthy controls. In the dominant model, HCC patients had significantly higher frequencies of rs548234 genotypes CC + TC than cirrhosis patients (P = .009). Rs548234 was an independent factor for HCC in comparison with either cirrhosis (P = .005) or all chronic HBV infection without HCC (P = .018). Functional annotation showed evidence of the role of the SNPs in gene regulation. In conclusion, through this study it is revealed for the first time that rs6568431 may be associated with susceptibility to HBV infection and that rs548234 may be associated with HCC risk in chronic HBV infection, supporting the presence of HBV-related disease-causing regulatory polymorphisms in the PRDM1-ATG5 intergenic region.

The prostate cancer risk variant rs55958994 regulates multiple gene expression through extreme long-range chromatin interaction to control tumor progression

Genome-wide association studies identified single-nucleotide polymorphism (SNP) rs55958994 as a significant variant associated with increased susceptibility to prostate cancer. However, the mechanisms by which this SNP mediates increased risk to cancer are still unknown. In this study, we show that this variant is located in an enhancer active in prostate cancer cells. Deletion of this enhancer from prostate tumor cells resulted in decreased tumor initiation, tumor growth, and invasive migration, as well as a loss of stem-like cells. Using a combination of capture chromosome conformation capture (Capture-C) and RNA sequencing, we identified genes on the same and different chromosomes as targets regulated by the enhancer. Furthermore, we show that expression of individual candidate target genes in an enhancer-deleted cell line rescued different aspects of tumorigenesis. Our data suggest that the rs55958994-associated enhancer affects prostate cancer progression by influencing expression of multiple genes via long-range chromatin interactions.

Integrated Analysis of the Gene Expression Changes During Colorectal Cancer Progression by Bioinformatic Methods

We attempted to analyze the aberrant pathways and genes underlying the successive stages of colorectal cancer (CRC). The CRC related microarray data (GSE77953) were retrieved from Gene Expression Omnibus database, which included 17 colonic adenoma, 17 carcinoma, 11 CRC metastases, and 13 normal colonic epithelium samples. The differential expression patterns in colonic adenoma, carcinoma, and metastases were analyzed. Gene functional interaction (FI) and coexpressed network were constructed. Pathway enrichment analysis was performed to investigate the perturbed pathways, and disease-related genes were explored based on the Comparative Toxicogenomics Database. Total 438 genes were identified to be differentially expressed in colonic adenoma, 885 in carcinoma and 736 in metastases. The upregulated genes in adenoma were significantly related with ribosome, oxidative phosphorylation, and protein export related pathways. The downregulated genes in carcinoma and metastases were enriched in the same pathways, such as nitrogen metabolism, mineral absorption, and steroid hormone biosynthesis. FI network was constructed with 219 and 3914 edges, which were further divided to 12 modules. The genes in module 0 were closely related with ribosome, protein export, and RNA transport. Coexpressed genes were enriched in ribosome, protein export, and mineral absorption pathways. Total eight common upregulated genes were found to be the CRC-related genes such as RNF43 (ring finger protein 43), EIF3H (eukaryotic translation initiation factor 3 subunit H), and STRAP (serine/threonine kinase receptor associated protein). The common downregulated genes included ABCG2 (ATP binding cassette subfamily G member 2), GCG (glucagon), and SULT1A1 (sulfotransferase family 1A member 1). Oxidative phosphorylation, nitrogen metabolism, mineral absorption, and protein synthesis may significantly be perturbed in the progression of CRC. The overexpression of EIF3H may be the predictor for CRC formation.

Genome-wide enhancer annotations differ significantly in genomic distribution, evolution, and function

Background

Non-coding gene regulatory enhancers are essential to transcription in mammalian cells. As a result, a large variety of experimental and computational strategies have been developed to identify cis-regulatory enhancer sequences. Given the differences in the biological signals assayed, some variation in the enhancers identified by different methods is expected; however, the concordance of enhancers identified by different methods has not been comprehensively evaluated. This is critically needed, since in practice, most studies consider enhancers identified by only a single method. Here, we compare enhancer sets from eleven representative strategies in four biological contexts.

Results

All sets we evaluated overlap significantly more than expected by chance; however, there is significant dissimilarity in their genomic, evolutionary, and functional characteristics, both at the element and base-pair level, within each context. The disagreement is sufficient to influence interpretation of candidate SNPs from GWAS studies, and to lead to disparate conclusions about enhancer and disease mechanisms. Most regions identified as enhancers are supported by only one method, and we find limited evidence that regions identified by multiple methods are better candidates than those identified by a single method. As a result, we cannot recommend the use of any single enhancer identification strategy in all settings.

Conclusions

Our results highlight the inherent complexity of enhancer biology and identify an important challenge to mapping the genetic architecture of complex disease. Greater appreciation of how the diverse enhancer identification strategies in use today relate to the dynamic activity of gene regulatory regions is needed to enable robust and reproducible results.

Electronic supplementary material

The online version of this article (10.1186/s12864-019-5779-x) contains supplementary material, which is available to authorized users.

Comprehensive functional annotation of susceptibility SNPs prioritized 10 genes for schizophrenia

Nearly 95% of susceptibility SNPs identified by genome-wide association studies (GWASs) are located in non-coding regions, which causes a lot of difficulty in deciphering their biological functions on disease pathogenesis. Here, we aimed to conduct a comprehensive functional annotation for all the schizophrenia susceptibility loci obtained from GWASs. Considering varieties of epigenomic regulatory elements, we annotated all 22,688 acquired susceptibility SNPs according to their genomic positions to obtain functional SNPs. The comprehensive annotation indicated that these functional SNPs are broadly involved in diverse biological processes. Histone modification enrichment showed that H3K27ac, H3K36me3, H3K4me1, and H3K4me3 were related to the development of schizophrenia. Transcription factors (TFs) prediction, methylation quantitative trait loci (meQTL) analyses, expression quantitative trait loci (eQTL) analyses, and proteomic quantitative trait loci analyses (pQTL) identified 447 target protein-coding genes. Subsequently, differential expression analyses between schizophrenia cases and controls, nervous system phenotypes from mouse models, and protein-protein interaction with known schizophrenia-related pathways and genes were carried out with our target genes. We finaly prioritized 10 target genes for schizophrenia (CACNA1C, CLU, CSNK2B, GABBR1, GRIN2A, MAPK3, NOTCH4, SRR, TNF, and SYNGAP1). Our results may serve as an encyclopedia of schizophrenia susceptibility SNPs and offer holistic guides for post-GWAS functional experiments.

Functional microRNA binding site variants

Germline single nucleotide polymorphisms are one of the most common genetic variations. Polymorphisms that cause nonsynonymous mutations in gene coding regions are known to cause serious deleterious downstream effects. However, even polymorphisms in noncoding regions can have profound functional consequences by disrupting essential regulatory sites. Specifically, polymorphisms that alter microRNA binding sites can disrupt the regulation of hallmark biological pathways implicated in tumorigenesis and tumor progression. Many of these microRNA-associated polymorphisms (miR-SNPs) have recently been shown to be important biomarkers of cancer risk, prognosis, and treatment outcomes. This review will summarize the functional impact of key miR-SNPs and define a subset of miR-SNPs that may be clinically useful prognostic or predictive biomarkers.

Using 3D epigenomic maps of primary olfactory neuronal cells from living individuals to understand gene regulation

As part of PsychENCODE, we developed a three-dimensional (3D) epigenomic map of primary cultured neuronal cells derived from olfactory neuroepithelium (CNON). We mapped topologically associating domains and high-resolution chromatin interactions using Hi-C and identified regulatory elements using chromatin immunoprecipitation and nucleosome positioning assays. Using epigenomic datasets from biopsies of 63 living individuals, we found that epigenetic marks at distal regulatory elements are more variable than marks at proximal regulatory elements. By integrating genotype and metadata, we identified enhancers that have different levels corresponding to differences in genetic variation, gender, smoking, and schizophrenia. Motif searches revealed that many CNON enhancers are bound by neuronal-related transcription factors. Last, we combined 3D epigenomic maps and gene expression profiles to predict enhancer-target gene interactions on a genome-wide scale. This study not only provides a framework for understanding individual epigenetic variation using a primary cell model system but also contributes valuable data resources for epigenomic studies of neuronal epithelium.

Regulatory SNPs and their widespread effects on the transcriptome

Currently, it is generally accepted that the cis-acting effects of noncoding variants on gene expression are a major factor for phenotypic variation in complex traits and disease susceptibility. Meanwhile, the protein products of many target genes for the identified cis-regulatory variants (rSNPs) are regulatory molecules themselves (transcription factors, effectors, components of signal transduction pathways, etc.), which implies dramatic downstream effects of these variations on complex gene networks. Here, we brief the results of recent most comprehensive studies on the role of rSNPs in transcriptional regulation across the genome.

The MYC Enhancer-ome: Long-Range Transcriptional Regulation of MYC in Cancer

MYC is one of the most important oncogenes in cancer. Indeed, MYC is upregulated in 50-60% of all tumors. MYC overexpression can be achieved through a variety of mechanisms, including gene duplications, chromosomal translocations, or somatic mutations leading to increased MYC stability. However, recent studies have identified numerous tissue-specific noncoding enhancers of MYC that play major roles in cancer, highlighting long-range transcriptional regulation of MYC as a critical novel mechanism leading to MYC hyperactivation and as a potential target for new therapeutic strategies in the near future. Here we summarize the regions and mechanisms involved in the long-range transcriptional regulation of MYC, underscoring the relevance of MYC enhancers both in normal physiological development and in MYC-driven cancer initiation and progression.

CRISPR-mediated deletion of prostate cancer risk-associated CTCF loop anchors identifies repressive chromatin loops

BACKGROUND

Recent genome-wide association studies (GWAS) have identified more than 100 loci associated with increased risk of prostate cancer, most of which are in non-coding regions of the genome. Understanding the function of these non-coding risk loci is critical to elucidate the genetic susceptibility to prostate cancer.

RESULTS

We generate genome-wide regulatory element maps and performed genome-wide chromosome confirmation capture assays (in situ Hi-C) in normal and tumorigenic prostate cells. Using this information, we annotate the regulatory potential of 2,181 fine-mapped prostate cancer risk-associated SNPs and predict a set of target genes that are regulated by prostate cancer risk-related H3K27Ac-mediated loops. We next identify prostate cancer risk-associated CTCF sites involved in long-range chromatin loops. We use CRISPR-mediated deletion to remove prostate cancer risk-associated CTCF anchor regions and the CTCF anchor regions looped to the prostate cancer risk-associated CTCF sites, and we observe up to 100-fold increases in expression of genes within the loops when the prostate cancer risk-associated CTCF anchor regions are deleted.

CONCLUSIONS

We identify GWAS risk loci involved in long-range loops that function to repress gene expression within chromatin loops. Our studies provide new insights into the genetic susceptibility to prostate cancer.

Positional integration of lung adenocarcinoma susceptibility loci with primary human alveolar epithelial cell epigenomes

AIM

To identify functional lung adenocarcinoma (LUAD) risk SNPs.

MATERIALS & METHODS

Eighteen validated LUAD risk SNPs (p ≤ 5 × 10^-8) and 930 SNPs in high linkage disequilibrium (r² > 0.5) were integrated with epigenomic information from primary human alveolar epithelial cells. Enhancer-associated SNPs likely affecting transcription factor-binding sites were predicted. Three SNPs were functionally investigated using luciferase assays, expression quantitative trait loci and cancer-specific expression.

RESULTS

Forty-seven SNPs mapped to putative enhancers; 11 located to open chromatin. Of these, seven altered predicted transcription factor-binding motifs. Rs6942067 showed allele-specific luciferase expression and expression quantitative trait loci analysis indicates that it influences expression of DCBLD1, a gene that encodes an unknown membrane protein and is overexpressed in LUAD.

CONCLUSION

Integration of candidate LUAD risk SNPS with epigenomic marks from normal alveolar epithelium identified numerous candidate functional LUAD risk SNPs including rs6942067, which appears to affect DCBLD1 expression. Data deposition: Data are provided in GEO record GSE84273.

Comprehensive review and annotation of susceptibility SNPs associated with obesity-related traits

We aimed to summarize the results of genetic association studies for obesity and provide a comprehensive annotation of all susceptibility single nucleotide polymorphisms (SNPs). A total of 72 studies were summarized, resulting in 90,361 susceptibility SNPs (738 index SNPs and 89,623 linkage disequilibrium SNPs). Over 90% of the susceptibility SNPs are located in non-coding regions, and it is challenging to understand their functional significance. Therefore, we annotated these SNPs by using various functional databases. We identified 24,623 functional SNPs, including 4 nonsense SNPs, 479 missense SNPs, 399 untranslated region SNPs which might affect microRNA binding, 262 promoter and 5,492 enhancer SNPs which might affect transcription factor binding, 7 splicing sites, 76 SNPs which might affect gene methylation levels, 1,839 SNPs under natural selection and 17,351 SNPs which might modify histone binding. Expression quantitative trait loci analyses for functional SNPs identified 98 target genes, including 69 protein coding genes, 27 long non-coding RNAs and 3 processed transcripts. The percentage of protein coding genes that could be correlated with obesity-related pathways directly or through gene-gene interaction is 75.36 (52/69). Our results may serve as an encyclopaedia of obesity susceptibility SNPs and offer guide for functional experiments.

A catalog of potential putative functional variants in psoriasis genome-wide association regions

Psoriasis is a common inflammatory skin disease, with considerable genetic contribution. Genome-wide association studies have successfully identified a number of genomic regions for the risk of psoriasis. However, it is challenging to pinpoint the functional causal variants and then further decipher the genetic mechanisms underlying each region. In order to prioritize potential functional causal variants within psoriasis susceptibility regions, we integrated the genetic association findings and functional genomic data publicly available, i.e. histone modifications in relevant immune cells. We characterized a pervasive enrichment pattern of psoriasis variants in five core histone marks across immune cells/tissues. We discovered that genetic alleles within psoriasis association regions might influence gene expression levels through significantly affecting the binding affinities of 17 transcription factors. We established a catalog of 654 potential functional causal variants for psoriasis and suggested that they significantly overlapped with causal variants for autoimmune diseases. We identified potential causal variant rs79824801 overlay with the peaks of five histone marks in primary CD4+ T cells. Its alternative allele affected the binding affinity of transcription factor IKZF1. This study highlights the complex genetic architecture and complicated mechanisms for psoriasis. The findings will inform the functional experiment design for psoriasis.

Integrative omics for health and disease

Advances in omics technologies - such as genomics, transcriptomics, proteomics and metabolomics - have begun to enable personalized medicine at an extraordinarily detailed molecular level. Individually, these technologies have contributed medical advances that have begun to enter clinical practice. However, each technology individually cannot capture the entire biological complexity of most human diseases. Integration of multiple technologies has emerged as an approach to provide a more comprehensive view of biology and disease. In this Review, we discuss the potential for combining diverse types of data and the utility of this approach in human health and disease. We provide examples of data integration to understand, diagnose and inform treatment of diseases, including rare and common diseases as well as cancer and transplant biology. Finally, we discuss technical and other challenges to clinical implementation of integrative omics.

CRISPR-based strategies for studying regulatory elements and chromatin structure in mammalian gene control

The development of high-throughput methods has enabled the genome-wide identification of putative regulatory elements in a wide variety of mammalian cells at an unprecedented resolution. Extensive genomic studies have revealed the important role of regulatory elements and genetic variation therein in disease formation and risk. In most cases, there is only correlative evidence for the roles of these elements and non-coding changes within these elements in pathogenesis. With the advent of genome- and epigenome-editing tools based on the CRISPR technology, it is now possible to test the functional relevance of the regulatory elements and alterations on a genomic scale. Here, we review the various CRISPR-based strategies that have been developed to functionally validate the candidate regulatory elements in mammals as well as the non-coding genetic variants found to be associated with human disease. We also discuss how these synthetic biology tools have helped to elucidate the role of three-dimensional nuclear architecture and higher-order chromatin organization in shaping functional genome and controlling gene expression.

A novel variant associated with HDL-C levels by modifying DAGLB expression levels: An annotation-based genome-wide association study

Although numbers of genome-wide association studies (GWAS) have been performed for serum lipid levels, limited heritability has been explained. Studies showed that combining data from GWAS and expression quantitative trait loci (eQTLs) signals can both enhance the discovery of trait-associated SNPs and gain a better understanding of the mechanism. We performed an annotation-based, multistage genome-wide screening for serum-lipid-level-associated loci in totally 6863 Han Chinese. A serum high-density lipoprotein cholesterol (HDL-C) associated variant rs1880118 (hg19 chr7:g. 6435220G>C) was replicated (P_combined = 1.4E-10). rs1880118 was associated with DAGLB (diacylglycerol lipase, beta) expression levels in subcutaneous adipose tissue (P = 5.9E-42) and explained 47.7% of the expression variance. After the replication, an active segment covering variants tagged by rs1880118 near 5' of DAGLB was annotated using histone modification and transcription factor binding signals. The luciferase report assay revealed that the segment containing the minor alleles showed increased transcriptional activity compared with segment contains the major alleles, which was consistent with the eQTL analyses. The expression-trait association tests indicated the association between the DAGLB and serum HDL-C levels using gene-based approaches called "TWAS" (P = 3.0E-8), "SMR" (P = 1.1E-4), and "Sherlock" (P = 1.6E-6). To summarize, we identified a novel HDL-C-associated variant which explained nearly half of the expression variance of DAGLB. Integrated analyses established a genotype-gene-phenotype three-way association and expanded our knowledge of DAGLB in lipid metabolism.

Association between the HOTAIR polymorphisms and cancer risk: an updated meta-analysis

Purpose

LncRNA HOTAIR plays an important role in many cancer. Several studies have shown that some HOTAIR SNPs might be associated with tumor risk in case-control studies, but the results are inconsistent and inconclusive. Therefore, it is necessary to better evaluate association between the HOTAIR SNPs and the risk of cancer.

Results

rs920778, rs7958904 and rs874945 but not rs4759314 and rs1899663 loci were significantly related to cancer risk, among of which rs920778 and rs874945 increased and rs7958904 decreased cancer risk, respectively. Moreover, rs920778 is significantly susceptible in both Asian population and digestive cancer risks.

Materials and Methods

Data were collected from PubMed, Embase and Web of Science. A total of 11 case-control studies were selected for the quantitative analysis. Software Stata (Version 12) was used to calculate Odds ratios (ORs) and 95% confidence intervals (CIs) to evaluate the strength of the associations. Subgroup analysis, sensitivity analysis, and publication bias were also performed. Five HOTAIR SNPs were finally enrolled in the study.

Conclusions

HOTAIR SNP rs920778, rs7958904 and rs874945 are susceptible to cancer risk. SNP rs920778 is also a useful risk factor in evaluation of Asian population and digestive cancer. In addition, the cancer risk SNP rs874945 is first reported in the meta-analysis.

Integrative expression quantitative trait locus-based analysis of colorectal cancer identified a functional polymorphism regulating SLC22A5 expression

Multiple single nucleotide polymorphisms (SNPs) have been found to be highly correlated with colorectal cancer (CRC) risk. However, the variants identified thus far only explain a small part of the cases, suggesting the existence of many uncharacterised genetic determinants. In this study, using the multilevel 'omics' data provided in The Cancer Genome Atlas, we systematically performed expression quantitative trait locus (eQTL) analysis for CRC and identified nine SNPs with significant effects on mRNA expression (correlation |r| > 0.3 and FDR < 0.01). Then we conducted a two-stage case-control study consisting of 1528 cases and 1528 controls to examine the associations between candidate SNPs and CRC risk. We found that rs27437 in SLC22A5 was significantly correlated with CRC risk in both stages and the combined study (additive model, OR = 1.31, 95%CI = 1.17-1.47, P = 1.97 × 10^-6). eQTL analysis showed that rs27437 GG and GA genotypes were associated with lower expression of SLC22A5 compared with the AA genotype. Dual-luciferase reporter assays confirmed that the G risk allele could decrease the expression of luciferase. SLC22A5 was significantly decreased in CRC tumour tissues compared with adjacent normal tissues, indicating that SLC22A5 may play important roles in CRC, and rs27437 in SLC22A5 might serve as a novel biomarker for early detection and prevention of CRC.

Assessing sufficiency and necessity of enhancer activities for gene expression and the mechanisms of transcription activation

Enhancers are important genomic regulatory elements directing cell type-specific transcription. They assume a key role during development and disease, and their identification and functional characterization have long been the focus of scientific interest. The advent of next-generation sequencing and clustered regularly interspaced short palindromic repeat (CRISPR)/Cas9-based genome editing has revolutionized the means by which we study enhancer biology. In this review, we cover recent developments in the prediction of enhancers based on chromatin characteristics and their identification by functional reporter assays and endogenous DNA perturbations. We discuss that the two latter approaches provide different and complementary insights, especially in assessing enhancer sufficiency and necessity for transcription activation. Furthermore, we discuss recent insights into mechanistic aspects of enhancer function, including findings about cofactor requirements and the role of post-translational histone modifications such as monomethylation of histone H3 Lys4 (H3K4me1). Finally, we survey how these approaches advance our understanding of transcription regulation with respect to promoter specificity and transcriptional bursting and provide an outlook covering open questions and promising developments.

Identification of a functional variant for colorectal cancer risk mapping to chromosome 5q31.1

Genome-wide association studies (GWASs) have established chromosome 5q31.1 as a risk locus for colorectal cancer (CRC). We previously identified a potentially regulatory single nucleotide polymorphism (SNP) rs17716310 within 5q31.1. Now, we extended our study with another independent Chinese population, functional assays and analyses of TCGA (The Cancer Genome Atlas) data. Significant associations between rs17716310 and CRC risk were found in Present Study including 1075 CRC cases and 1999 controls (additive model: OR = 1.149, 95% CI = 1.027–1.286, P = 0.016), and in Combined Study including 1766 cases and 2708 controls (additive model: OR = 1.145, 95% CI = 1.045–1.254, P = 0.004). Dual luciferase reporter gene assays indicated that the variant C allele obviously increased transcriptional activity. Using TCGA datasets, we indicated rs17716310 as a cis expression quantitative trait locus (eQTL) for the gene SMAD5, whose expression was significantly higher in CRC tissues. These findings suggested that the functional polymorphism rs17716310 A > C might be a genetic modifier for CRC, promoting the expression of SMAD5 that belonged to the transforming growth factor beta (TGF-β) signaling pathway.

Systematic identification of regulatory variants associated with cancer risk

Background

Most cancer risk-associated single nucleotide polymorphisms (SNPs) identified by genome-wide association studies (GWAS) are noncoding and it is challenging to assess their functional impacts. To systematically identify the SNPs that affect gene expression by modulating activities of distal regulatory elements, we adapt the self-transcribing active regulatory region sequencing (STARR-seq) strategy, a high-throughput technique to functionally quantify enhancer activities.

Results

From 10,673 SNPs linked with 996 cancer risk-associated SNPs identified in previous GWAS studies, we identify 575 SNPs in the fragments that positively regulate gene expression, and 758 SNPs in the fragments with negative regulatory activities. Among them, 70 variants are regulatory variants for which the two alleles confer different regulatory activities. We analyze in depth two regulatory variants—breast cancer risk SNP rs11055880 and leukemia risk-associated SNP rs12142375—and demonstrate their endogenous regulatory activities on expression of ATF7IP and PDE4B genes, respectively, using a CRISPR-Cas9 approach.

Conclusions

By identifying regulatory variants associated with cancer susceptibility and studying their molecular functions, we hope to help the interpretation of GWAS results and provide improved information for cancer risk assessment.

Electronic supplementary material

The online version of this article (doi:10.1186/s13059-017-1322-z) contains supplementary material, which is available to authorized users.