SNPs associated with colorectal cancer at 15q13.3 affect risk enhancers that modulate GREM1 gene expression

Prioritization of risk genes in colorectal cancer by integrative analysis of multi-omics data and gene networks

Genome-wide association studies (GWASs) have identified over 140 colorectal cancer (CRC)-associated loci; however, target genes at the majority of loci and underlying molecular mechanisms are poorly understood. Here, we utilized a Bayesian approach, integrative risk gene selector (iRIGS), to prioritize risk genes at CRC GWAS loci by integrating multi-omics data. As a result, a total of 105 high-confidence risk genes (HRGs) were identified, which exhibited strong gene dependencies for CRC and enrichment in the biological processes implicated in CRC. Among the 105 HRGs, CEBPB, located at the 20q13.13 locus, acted as a transcription factor playing critical roles in cancer. Our subsequent assays indicated the tumor promoter function of CEBPB that facilitated CRC cell proliferation by regulating multiple oncogenic pathways such as MAPK, PI3K-Akt, and Ras signaling. Next, by integrating a fine-mapping analysis and three independent case-control studies in Chinese populations consisting of 8,039 cases and 12,775 controls, we elucidated that rs1810503, a putative functional variant regulating CEBPB, was associated with CRC risk (OR=0.90, 95%CI=0.86-0.93, P=1.07×10-7). The association between rs1810503 and CRC risk was further validated in three additional multi-ancestry populations consisting of 24,254 cases and 58,741 controls. Mechanistically, the rs1810503 A to T allele change weakened the enhancer activity in an allele-specific manner to decrease CEBPB expression via long-range promoter-enhancer interactions, mediated by the transcription factor, REST, and thus decreased CRC risk. In summary, our study provides a genetic resource and a generalizable strategy for CRC etiology investigation, and highlights the biological implications of CEBPB in CRC tumorigenesis, shedding new light on the etiology of CRC.

Assessment of Colorectal Cancer Risk Factors through the Application of Network-Based Approaches in a Racially Diverse Cohort of Colon Organoid Stem Cells

Numerous demographic factors have been associated with colorectal cancer (CRC) risk. To better define biological mechanisms underlying these associations, we performed RNA sequencing of stem-cell-enriched organoids derived from the healthy colons of seven European Americans and eight African Americans. A weighted gene co-expression network analysis was performed following RNA sequencing. Module-trait relationships were determined through the association testing of each module and five CRC risk factors (age, body mass index, sex, smoking history, and race). Only modules that displayed a significantly positive correlation for gene significance and module membership were considered for further investigation. In total, 16 modules were associated with known CRC risk factors (p < 0.05). To contextualize the role of risk modules in CRC, publicly available RNA-sequencing data from TCGA-COAD were downloaded and re-analyzed. Differentially expressed genes identified between tumors and matched normal-adjacent tissue were overlaid across each module. Loci derived from CRC genome-wide association studies were additionally overlaid across modules to identify robust putative targets of risk. Among them, MYBL2 and RXRA represented strong plausible drivers through which cigarette smoking and BMI potentially modulated CRC risk, respectively. In summary, our findings highlight the potential of the colon organoid system in identifying novel CRC risk mechanisms in an ancestrally diverse and cellularly relevant population.

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.

A risk variant for Barrett's esophagus and esophageal adenocarcinoma at chr8p23.1 affects enhancer activity and implicates multiple gene targets

Nineteen genetic susceptibility loci for esophageal adenocarcinoma (EAC) and its precursor Barrett's esophagus (BE) have been identified through genome-wide association studies (GWAS). Clinical translation of such discoveries, however, has been hindered by the slow pace of discovery of functional/causal variants and gene targets at these loci. We previously developed a systematic informatics pipeline to prioritize candidate functional variants using functional potential scores, applied the pipeline to select high-scoring BE/EAC risk loci, and validated a functional variant at chr19p13.11 (rs10423674). Here, we selected two additional prioritized loci for experimental interrogation: chr3p13/rs1522552 and chr8p23.1/rs55896564. Candidate enhancer regions encompassing these variants were evaluated using luciferase reporter assays in two EAC cell lines. One of the two regions tested exhibited allele-specific enhancer activity - 8p23.1/rs55896564. CRISPR-mediated deletion of the putative enhancer in EAC cell lines correlated with reduced expression of three candidate gene targets: B lymphocyte kinase (BLK), nei like DNA glycosylase 2 (NEIL2), and cathepsin B (CTSB). Expression quantitative trait locus (eQTL) mapping in normal esophagus and stomach revealed strong associations between the BE/EAC risk allele at rs55896564 (G) and lower expression of CTSB, a protease gene implicated in epithelial wound repair. These results further support the utility of functional potential scores for GWAS variant prioritization, and provide the first experimental evidence of a functional variant and risk enhancer at the 8p23.1 GWAS locus. Identification of CTSB, BLK, and NEIL2 as candidate gene targets suggests that altered expression of these genes may underlie the genetic risk association at 8p23.1 with BE/EAC.

The landscape of GWAS validation; systematic review identifying 309 validated non-coding variants across 130 human diseases

Background

The remarkable growth of genome-wide association studies (GWAS) has created a critical need to experimentally validate the disease-associated variants, 90% of which involve non-coding variants.

Methods

To determine how the field is addressing this urgent need, we performed a comprehensive literature review identifying 36,676 articles. These were reduced to 1454 articles through a set of filters using natural language processing and ontology-based text-mining. This was followed by manual curation and cross-referencing against the GWAS catalog, yielding a final set of 286 articles.

Results

We identified 309 experimentally validated non-coding GWAS variants, regulating 252 genes across 130 human disease traits. These variants covered a variety of regulatory mechanisms. Interestingly, 70% (215/309) acted through cis-regulatory elements, with the remaining through promoters (22%, 70/309) or non-coding RNAs (8%, 24/309). Several validation approaches were utilized in these studies, including gene expression (n = 272), transcription factor binding (n = 175), reporter assays (n = 171), in vivo models (n = 104), genome editing (n = 96) and chromatin interaction (n = 33).

Conclusions

This review of the literature is the first to systematically evaluate the status and the landscape of experimentation being used to validate non-coding GWAS-identified variants. Our results clearly underscore the multifaceted approach needed for experimental validation, have practical implications on variant prioritization and considerations of target gene nomination. While the field has a long way to go to validate the thousands of GWAS associations, we show that progress is being made and provide exemplars of validation studies covering a wide variety of mechanisms, target genes, and disease areas.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12920-022-01216-w.

Gene regulation by antisense transcription: A focus on neurological and cancer diseases

Advances in high-throughput sequencing over the past decades have led to the identification of thousands of non-coding RNAs (ncRNAs), which play a major role in regulating gene expression. One emerging class of ncRNAs is the natural antisense transcripts (NATs), the RNA molecules transcribed from the opposite strand of a protein-coding gene locus. NATs are known to concordantly and discordantly regulate gene expression in both cis and trans manners at the transcriptional, post-transcriptional, translational, and epigenetic levels. Aberrant expression of NATs can therefore cause dysregulation in many biological pathways and has been observed in many genetic diseases. This review outlines the involvements and mechanisms of NATs in the pathogenesis of various diseases, with a special emphasis on neurodegenerative diseases and cancer. We also summarize recent findings on NAT knockdown and/or overexpression experiments and discuss the potential of NATs as promising targets for future gene therapies.

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.

Prioritization and functional analysis of GWAS risk loci for Barrett's esophagus and esophageal adenocarcinoma

Genome-wide association studies (GWAS) have identified ~ 20 genetic susceptibility loci for esophageal adenocarcinoma (EAC), and its precursor, Barrett's esophagus (BE). Despite such advances, functional/causal variants and gene targets at these loci remain undefined, hindering clinical translation. A key challenge is that most causal variants map to non-coding regulatory regions such as enhancers, and typically, numerous potential candidate variants at GWAS loci require testing. We developed a systematic informatics pipeline for prioritizing candidate functional variants via integrative functional potential scores consolidated from multi-omics annotations, and used this pipeline to identify two high-scoring variants for experimental interrogation: chr9q22.32/rs11789015 and chr19p13.11/rs10423674. Minimal candidate enhancer regions spanning these variants were evaluated using luciferase reporter assays in two EAC cell lines. One of the two variants tested (rs10423674) exhibited allele-specific enhancer activity. CRISPR-mediated deletion of the putative enhancer region in EAC cell lines correlated with reduced expression of two genes-CREB-regulated transcription coactivator 1 (CRTC1) and Cartilage oligomeric matrix protein (COMP); expression of five other genes remained unchanged (CRLF1, KLHL26, TMEM59L, UBA52, RFXANK). Expression quantitative trait locus (eQTL) mapping indicated that rs10423674 genotype correlated with CRTC1 and COMP expression in normal esophagus. This study represents the first experimental effort to bridge GWAS associations to biology in BE/EAC, and supports the utility of functional potential scores to guide variant prioritization. Our findings reveal a functional variant and candidate risk enhancer at chr19p13.11, and implicate CRTC1 and COMP as putative gene targets, suggesting that altered expression of these genes may underlie the BE/EAC risk association.

SNPs associated with colorectal cancer at 15q13.3 affect risk enhancers that modulate GREM1 gene expression

Several genome wide association studies of colorectal cancer (CRC) have identified single nucleotide polymorphisms (SNPs) on chromosome 15q13.3 associated with CRC risk. To identify functional variant(s) underlying this association, we investigated SNPs in linkage disequilibrium with the risk‐associated SNP rs4779584 that overlapped regulatory regions/enhancer elements characterized in colon‐related tissues and cells. We identified several SNP‐containing regulatory regions that exhibited enhancer activity in vitro, including one SNP (rs1406389) that correlated with allele‐specific effects on enhancer activity. Deletion of either this enhancer or another enhancer that had previously been reported in this region correlated with decreased expression of GREM1 following CRISPR/Cas9 genome editing. That GREM1 is one target of these enhancers was further supported by an expression quantitative trait loci correlation between rs1406389 and GREM1 expression in the transverse but not sigmoid colon in the Genotype‐Tissue Expression dataset. Taken together, we conclude that the 15q13.3 region contains at least two functional variants that map to distinct enhancers and impact CRC risk through modulation of GREM1 expression.