Expression Quantitative Trait Loci (eQTL) Mapping in Korean Patients With Crohn's Disease and Identification of Potential Causal Genes Through Integration With Disease Associations

Newcastle disease virus harboring the PTEN gene inhibits pancreatic cancer growth by inhibiting PI3K/AKT/mTOR signaling and activating apoptosis

Pancreatic ductal adenocarcinoma (PDAC) is a highly aggressive and intractable cancer that requires more effective therapies that can improve early detection, enhance treatment efficacy, and provide better patient outcomes. Kirsten rat sarcoma viral oncogene homolog (KRAS) mutations and reduced phosphatase and tensin homolog (PTEN) protein expression are key factors driving the proliferation and severity of PDAC. To address this, a recombinant Newcastle disease virus (rNDV) containing the PTEN gene (rNDV-PTEN) was created to investigate its PDAC cell-killing and tumor-suppression effects in PDAC cells transplanted into mice. PTEN expression induced by rNDV-PTEN virus infection in KRAS-mutated PDAC cells lowered phosphatidylinositol 3-kinases (PI3K)/protein kinase B (AKT)/mammalian target of rapamycin (mTOR) signaling, promoted PDAC cell death, and suppressed tumor growth. PTEN overexpression promotes apoptotic signaling pathways in both PANC-1 cells and orthotopic xenograft mice. Additionally, during virotherapy, rNDV-PTEN-injected mice exhibited a mild immune response and no abnormal responses in blood parameters such as glucose, triglyceride, and total cholesterol levels. These findings support the potential of rNDV-PTEN as a safe and effective therapy for PDAC with highly activated PI3K/AKT/mTOR signaling caused by KRAS and PTEN gene mutations. Thus, PTEN gene-containing rNDV may be a promising candidate for pancreatic cancer treatment.

Transcriptome-wide association studies associated with Crohn's disease: challenges and perspectives

Crohn's disease (CD) is regarded as a lifelong progressive disease affecting all segments of the intestinal tract and multiple organs. Based on genome-wide association studies (GWAS) and gene expression data, transcriptome-wide association studies (TWAS) can help identify susceptibility genes associated with pathogenesis and disease behavior. In this review, we overview seven reported TWASs of CD, summarize their study designs, and discuss the key methods and steps used in TWAS, which affect the prioritization of susceptibility genes. This article summarized the screening of tissue-specific susceptibility genes for CD, and discussed the reported potential pathological mechanisms of overlapping susceptibility genes related to CD in a certain tissue type. We observed that ileal lipid-related metabolism and colonic extracellular vesicles may be involved in the pathogenesis of CD by performing GO pathway enrichment analysis for susceptibility genes. We further pointed the low reproducibility of TWAS associated with CD and discussed the reasons for these issues, strategies for solving them. In the future, more TWAS are needed to be designed into large-scale, unified cohorts, unified analysis pipelines, and fully classified databases of expression trait loci.

Transcriptomics : Approaches to Quantifying Gene Expression and Their Application to Studying the Human Brain

To date, the field of transcriptomics has been characterized by rapid methods development and technological advancement, with new technologies continuously rendering older ones obsolete.This chapter traces the evolution of approaches to quantifying gene expression and provides an overall view of the current state of the field of transcriptomics, its applications to the study of the human brain, and its place in the broader emerging multiomics landscape.

Identification of asthma-related genes using asthmatic blood eQTLs of Korean patients

BACKGROUND

More than 200 asthma-associated genetic variants have been identified in genome-wide association studies (GWASs). Expression quantitative trait loci (eQTL) data resources can help identify causal genes of the GWAS signals, but it can be difficult to find an eQTL that reflects the disease state because most eQTL data are obtained from normal healthy subjects.

METHODS

We performed a blood eQTL analysis using transcriptomic and genotypic data from 433 Korean asthma patients. To identify asthma-related genes, we carried out colocalization, Summary-based Mendelian Randomization (SMR) analysis, and Transcriptome-Wide Association Study (TWAS) using the results of asthma GWASs and eQTL data. In addition, we compared the results of disease eQTL data and asthma-related genes with two normal blood eQTL data from Genotype-Tissue Expression (GTEx) project and a Japanese study.

RESULTS

We identified 340,274 cis-eQTL and 2,875 eGenes from asthmatic eQTL analysis. We compared the disease eQTL results with GTEx and a Japanese study and found that 64.1% of the 2,875 eGenes overlapped with the GTEx eGenes and 39.0% with the Japanese eGenes. Following the integrated analysis of the asthmatic eQTL data with asthma GWASs, using colocalization and SMR methods, we identified 15 asthma-related genes specific to the Korean asthmatic eQTL data.

CONCLUSIONS

We provided Korean asthmatic cis-eQTL data and identified asthma-related genes by integrating them with GWAS data. In addition, we suggested these asthma-related genes as therapeutic targets for asthma. We envisage that our findings will contribute to understanding the etiological mechanisms of asthma and provide novel therapeutic targets.

Multi-Omics Studies in Historically Excluded Populations: The Road to Equity

Over the past few decades, genomewide association studies (GWASs) have identified the specific genetics variants contributing to many complex diseases by testing millions of genetic variations across the human genome against a variety of phenotypes. However, GWASs are limited in their ability to uncover mechanistic insight given that most significant associations are found in non-coding region of the genome. Furthermore, the lack of diversity in studies has stymied the advance of precision medicine for many historically excluded populations. In this review, we summarize most popular multi-omics approaches (genomics, transcriptomics, proteomics, and metabolomics) related to precision medicine and highlight if diverse populations have been included and how their findings have advance biological understanding of disease and drug response. New methods that incorporate local ancestry have been to improve the power of GWASs for admixed populations (such as African Americans and Latinx). Because most signals from GWAS are in the non-coding region, other machine learning and omics approaches have been developed to identify the potential causative single-nucleotide polymorphisms and genes that explain these phenotypes. These include polygenic risk scores, expression quantitative trait locus mapping, and transcriptome-wide association studies. Analogous protein methods, such as proteins quantitative trait locus mapping, proteome-wide association studies, and metabolomic approaches provide insight into the consequences of genetic variation on protein abundance. Whereas, integrated multi-omics studies have improved our understanding of the mechanisms for genetic association, we still lack the datasets and cohorts for historically excluded populations to provide equity in precision medicine and pharmacogenomics.

Case-case genome-wide association analysis identifying genetic loci with divergent effects on Crohn's disease and ulcerative colitis

Crohn's disease (CD) and ulcerative colitis (UC), two major subtypes of inflammatory bowel disease, show substantial differences in their clinical course and treatment response. To identify the genetic factors underlying the distinct characteristics of these two diseases, we performed a genome-wide association study (GWAS) between CD (n = 2359) and UC (n = 2175) in a Korean population, followed by replication in an independent sample of 772 CD and 619 UC cases. Two novel loci were identified with divergent effects on CD and UC: rs9842650 in CD200 and rs885026 in NCOR2. In addition, the seven established susceptibility loci [major histocompatibility complex (MHC), TNFSF15, OTUD3, USP12, IL23R, FCHSD2 and RIPK2] reached genome-wide significance. Of the nine loci, six (MHC, TNFSF15, OTUD3, USP12, IL23R and CD200) were replicated in the case-case GWAS of European populations. The proportion of variance explained in CD-UC status by polygenic risk score analysis was up to 22.6%. The area under the receiver-operating characteristic curve value was 0.74, suggesting acceptable discrimination between CD and UC. This CD-UC GWAS provides new insights into genetic differences between the two diseases with similar symptoms and might be useful in improving their diagnosis and treatment.

Аpplication of massive parallel reporter analysis in biotechnology and medicine

The development and functioning of an organism relies on tissue-specific gene programs. Genome regulatory elements play a key role in the regulation of such programs, and disruptions in their function can lead to the development of various pathologies, including cancers, malformations and autoimmune diseases. The emergence of high-throughput genomic studies has led to massively parallel reporter analysis (MPRA) methods, which allow the functional verification and identification of regulatory elements on a genome-wide scale. Initially MPRA was used as a tool to investigate fundamental aspects of epigenetics, but the approach also has great potential for clinical and practical biotechnology. Currently, MPRA is used for validation of clinically significant mutations, identification of tissue-specific regulatory elements, search for the most promising loci for transgene integration, and is an indispensable tool for creating highly efficient expression systems, the range of application of which extends from approaches for protein development and design of next-generation therapeutic antibody superproducers to gene therapy. In this review, the main principles and areas of practical application of high-throughput reporter assays will be discussed.

High-resolution sequencing of nine elite upland cotton cultivars uncovers genic variations and breeding improvement targets

Structural variations (SVs) are critical factors affecting genome evolution and important traits. However, identification results and functional analyses of SVs in upland cotton are rare. Here, based on the genetic relationships, breeding history and cumulative planting area of upland cotton in China, nine predominant cultivars from the past 60 years (1950s-2010s) were selected for long read sequencing to uncover genic variations and breeding improvement targets for this crop. Based on the ZM24 reference genome, 0.88-1.47 × 10⁴ SVs per cultivar were identified, and an SV set was constructed. SVs affected the expression of a large number of genes during fiber elongation, and a transposable element insertion resulted in the glandless phenotype in upland cotton. Six widespread inversions were identified based on nine draft genomes and high-throughput chromosome conformation capture data. Multiple haplotype blocks that were always associated with aggregated SVs were demonstrated to play a pivotal role in the agronomic traits of upland cotton and drove its adaptation to the northern planting region. Exotic introgression was the source of these haplotype blocks and increased the genetic diversity of upland cotton. Our results enrich the genome resources of upland cotton, and the identified SVs will promote genetic and breeding research in cotton.

Integration of eQTL Analysis and GWAS Highlights Regulation Networks in Cotton under Stress Condition

The genus Gossypium is one of the most economically important crops in the world. Here, we used RNA-seq to quantify gene expression in a collection of G. arboreum seedlings and performed eGWAS on 28,382 expressed genes. We identified a total of 30,089 eQTLs in 10,485 genes, of which >90% were trans-regulate target genes. Using luciferase assays, we confirmed that different cis-eQTL haplotypes could affect promoter activity. We found ~6600 genes associated with ~1300 eQTL hotspots. Moreover, hotspot 309 regulates the expression of 325 genes with roles in stem length, fresh weight, seed germination rate, and genes related to cell wall biosynthesis and salt stress. Transcriptome-wide association study (TWAS) identified 19 candidate genes associated with the cotton growth and salt stress response. The variation in gene expression across the population played an essential role in population differentiation. Only a small number of the differentially expressed genes between South China, the Yangtze River region, and the Yellow River region sites were located in different chromosomal regions. The eQTLs found across the duplicated gene pairs showed conservative cis- or trans- regulation and that the expression levels of gene pairs were correlated. This study provides new insights into the evolution of gene expression regulation in cotton, and identifies eQTLs in stress-related genes for use in breeding improved cotton varieties.

Development of a Machine Learning Model to Predict Non-Durable Response to Anti-TNF Therapy in Crohn’s Disease Using Transcriptome Imputed from Genotypes

Almost half of patients show no primary or secondary response to monoclonal anti-tumor necrosis factor α (anti-TNF) antibody treatment for inflammatory bowel disease (IBD). Thus, the exact mechanisms of a non-durable response (NDR) remain inadequately defined. We used our genome-wide genotype data to impute expression values as features in training machine learning models to predict a NDR. Blood samples from various IBD cohorts were used for genotyping with the Korea Biobank Array. A total of 234 patients with Crohn’s disease (CD) who received their first anti-TNF therapy were enrolled. The expression profiles of 6294 genes in whole-blood tissue imputed from the genotype data were combined with clinical parameters to train a logistic model to predict the NDR. The top two and three most significant features were genetic features (DPY19L3, GSTT1, and NUCB1), not clinical features. The logistic regression of the NDR vs. DR status in our cohort by the imputed expression levels showed that the β coefficients were positive for DPY19L3 and GSTT1, and negative for NUCB1, concordant with the known eQTL information. Machine learning models using imputed gene expression features effectively predicted NDR to anti-TNF agents in patients with CD.

Lupus Susceptibility Region Containing CDKN1B rs34330 Mechanistically Influences Expression and Function of Multiple Target Genes, Also Linked to Proliferation and Apoptosis

OBJECTIVE

In a recent genome-wide association study, a significant genetic association between rs34330 of CDKN1B and risk of systemic lupus erythematosus (SLE) in Han Chinese was identified. This study was undertaken to validate the reported association and elucidate the biochemical mechanisms underlying the effect of the variant.

METHODS

We performed an allelic association analysis in patients with SLE, followed by a meta-analysis assessing genome-wide association data across 11 independent cohorts (n = 28,872). In silico bioinformatics analysis and experimental validation in SLE-relevant cell lines were applied to determine the functional consequences of rs34330.

RESULTS

We replicated a genetic association between SLE and rs34330 (meta-analysis P = 5.29 × 10-22 , odds ratio 0.84 [95% confidence interval 0.81-0.87]). Follow-up bioinformatics and expression quantitative trait locus analysis suggested that rs34330 is located in active chromatin and potentially regulates several target genes. Using luciferase and chromatin immunoprecipitation-real-time quantitative polymerase chain reaction, we demonstrated substantial allele-specific promoter and enhancer activity, and allele-specific binding of 3 histone marks (H3K27ac, H3K4me3, and H3K4me1), RNA polymerase II (Pol II), CCCTC-binding factor, and a critical immune transcription factor (interferon regulatory factor 1 [IRF-1]). Chromosome conformation capture revealed long-range chromatin interactions between rs34330 and the promoters of neighboring genes APOLD1 and DDX47, and effects on CDKN1B and the other target genes were directly validated by clustered regularly interspaced short palindromic repeat (CRISPR)-based genome editing. Finally, CRISPR/dead CRISPR-associated protein 9-based epigenetic activation/silencing confirmed these results. Gene-edited cell lines also showed higher levels of proliferation and apoptosis.

CONCLUSION

Collectively, these findings suggest a mechanism whereby the rs34330 risk allele (C) influences the presence of histone marks, RNA Pol II, and IRF-1 transcription factor to regulate expression of several target genes linked to proliferation and apoptosis. This process could potentially underlie the association of rs34330 with SLE.

Identification of Three Novel Susceptibility Loci for Inflammatory Bowel Disease in Koreans in an Extended Genome-Wide Association Study

BACKGROUND AND AIMS

Genome-wide association studies [GWAS] of inflammatory bowel disease [IBD] in multiple populations have identified over 240 susceptibility loci. We previously performed a largest-to-date Asian-specific IBD GWAS to identify two new IBD risk loci and confirm associations with 28 established loci. To identify additional susceptibility loci in Asians, we expanded our previous study design by doubling the case size with an additional dataset of 1726 cases and 378 controls.

METHODS

An inverse-variance fixed-effects meta-analysis was performed between the previous and the new GWAS dataset, comprising a total of 3195 cases and 4419 controls, followed by replication in an additional 1088 cases and 845 controls.

RESULTS

The meta-analysis of Korean GWAS identified one novel locus for ulcerative colitis at rs76227733 on 10q24 [pcombined = 6.56 × 10-9] and two novel loci for Crohn's disease [CD] at rs2240751 on 19p13 [pcombined = 3.03 × 10-8] and rs6936629 on 6q22 [pcombined = 3.63 × 10-8]. Pathway-based analysis of GWAS data using MAGMA showed that the MHC and antigenic stimulus-related pathways were more significant in Korean CD, whereas cytokine and transcription factor-related pathways were more significant in European CD. Phenotype variance explained by the polygenic risk scores derived from Korean data explained up to 14% of the variance of CD whereas those derived from European data explained 10%, emphasizing the need for large-scale genetic studies in this population.

CONCLUSIONS

The identification of novel loci not previously associated with IBD suggests the importance of studying IBD genetics in diverse populations.

GPR35 in Intestinal Diseases: From Risk Gene to Function

Diet and gut microbial metabolites mediate host immune responses and are central to the maintenance of intestinal health. The metabolite-sensing G-protein coupled receptors (GPCRs) bind metabolites and trigger signals that are important for the host cell function, survival, proliferation and expansion. On the contrary, inadequate signaling of these metabolite-sensing GPCRs most likely participate to the development of diseases including inflammatory bowel diseases (IBD). In the intestine, metabolite-sensing GPCRs are highly expressed by epithelial cells and by specific subsets of immune cells. Such receptors provide an important link between immune system, gut microbiota and metabolic system. Member of these receptors, GPR35, a class A rhodopsin-like GPCR, has been shown to be activated by the metabolites tryptophan-derived kynurenic acid (KYNA), the chemokine CXCL17 and phospholipid derivate lysophosphatidic acid (LPA) species. There have been studies on GPR35 in the context of intestinal diseases since its identification as a risk gene for IBD. In this review, we discuss the pharmacology of GPR35 including its proposed endogenous and synthetic ligands as well as its antagonists. We elaborate on the risk variants of GPR35 implicated in gut-related diseases and the mechanisms by which GPR35 contribute to intestinal homeostasis.

Therapeutic Opportunities and Challenges in Targeting the Orphan G Protein-Coupled Receptor GPR35

GPR35 is a class A, rhodopsin-like G protein-coupled receptor (GPCR) first identified more than 20 years ago. In the intervening period, identification of strong expression in the lower intestine and colon, in a variety of immune cells including monocytes and a variety of dendritic cells, and in dorsal root ganglia has suggested potential therapeutic opportunities in targeting this receptor in a range of conditions. GPR35 is, however, unusual in a variety of ways that challenge routes to translation. These include the following: (i) Although a substantial range and diversity of endogenous ligands have been suggested as agonist partners for this receptor, it officially remains defined as an "orphan" GPCR. (ii) Humans express two distinct protein isoform sequences, while rodents express only a single form. (iii) The pharmacologies of the human and rodent orthologues of GPR35 are very distinct, with variation between rat and mouse GPR35 being as marked as that between either of these species and the human forms. Herein we provide perspectives on each of the topics above as well as suggesting ways to overcome the challenges currently hindering potential translation. These include a better understanding of the extent and molecular basis for species selective GPR35 pharmacology and the production of novel mouse models in which both "on-target" and "off-target" effects of presumptive GPR35 ligands can be better defined, as well as a clear understanding of the human isoform expression profile and its significance at both tissue and individual cell levels.