The genetic regulation of the gastric transcriptome is associated with metabolic and obesity-related traits and diseases.

Tissue-specific gene expression and gene regulation lead to a better understanding of tissue-specific physiology and pathophysiology. We analyzed the transcriptome and genetic regulatory profiles of two distinct gastric sites, corpus and antrum, to identify tissue-specific gene expression and its regulation. Gastric corpus and antrum mucosa biopsies were collected during routine gastroscopies from up to 431 healthy individuals. We obtained genotype and transcriptome data and performed transcriptome profiling and expression quantitative trait locus (eQTL) studies. We further used data from genome-wide association studies (GWAS) of various diseases and traits to partition their heritability and to perform transcriptome-wide association studies (TWAS). The transcriptome data from corpus and antral mucosa highlights the heterogeneity of gene expression in the stomach. We identified enriched pathways revealing distinct and common physiological processes in gastric corpus and antrum. Furthermore, we found an enrichment of the single nucleotide polymorphism (SNP)-based heritability of metabolic, obesity-related, and cardiovascular traits and diseases by considering corpus- and antrum-specifically expressed genes. Particularly, we could prioritize gastric-specific candidate genes for multiple metabolic traits, like NQO1 which is involved in glucose metabolism, MUC1 which contributes to purine and protein metabolism or RAB27B being a regulator of weight and body composition. Our findings show that gastric corpus and antrum vary in their transcriptome and genetic regulatory profiles indicating physiological differences which are mostly related to digestion and epithelial protection. Moreover, our findings demonstrate that the genetic regulation of the gastric transcriptome is linked to biological mechanisms associated with metabolic, obesity-related, and cardiovascular traits and diseases. NEW & NOTEWORTHY We analyzed the transcriptomes and genetic regulatory profiles of gastric corpus and for the first time also of antrum mucosa in 431 healthy individuals. Through tissue-specific gene expression and eQTL analyses, we uncovered unique and common physiological processes across both primary gastric sites. Notably, our findings reveal that stomach-specific eQTLs are enriched in loci associated with metabolic traits and diseases, highlighting the pivotal role of gene expression regulation in gastric physiology and potential pathophysiology.

Pathogen-specific innate immune response patterns are distinctly affected by genetic diversity.

Innate immune responses vary by pathogen and host genetics. We analyze quantitative trait loci (eQTLs) and transcriptomes of monocytes from 215 individuals stimulated by fungal, Gram-negative or Gram-positive bacterial pathogens. We identify conserved monocyte responses to bacterial pathogens and a distinct antifungal response. These include 745 response eQTLs (reQTLs) and corresponding genes with pathogen-specific effects, which we find first in samples of male donors and subsequently confirm for selected reQTLs in females. reQTLs affect predominantly upregulated genes that regulate immune response via e.g., NOD-like, C-type lectin, Toll-like and complement receptor-signaling pathways. Hence, reQTLs provide a functional explanation for individual differences in innate response patterns. Our identified reQTLs are also associated with cancer, autoimmunity, inflammatory and infectious diseases as shown by external genome-wide association studies. Thus, reQTLs help to explain interindividual variation in immune response to infection and provide candidate genes for variants associated with a range of diseases.

Dissecting the genetic heterogeneity of gastric cancer.

BACKGROUND: Gastric cancer (GC) is clinically heterogenous according to location (cardia/non-cardia) and histopathology (diffuse/intestinal). We aimed to characterize the genetic risk architecture of GC according to its subtypes. Another aim was to examine whether cardia GC and oesophageal adenocarcinoma (OAC) and its precursor lesion Barrett's oesophagus (BO), which are all located at the gastro-oesophageal junction (GOJ), share polygenic risk architecture. METHODS: We did a meta-analysis of ten European genome-wide association studies (GWAS) of GC and its subtypes. All patients had a histopathologically confirmed diagnosis of gastric adenocarcinoma. For the identification of risk genes among GWAS loci we did a transcriptome-wide association study (TWAS) and expression quantitative trait locus (eQTL) study from gastric corpus and antrum mucosa. To test whether cardia GC and OAC/BO share genetic aetiology we also used a European GWAS sample with OAC/BO. FINDINGS: Our GWAS consisting of 5816 patients and 10,999 controls highlights the genetic heterogeneity of GC according to its subtypes. We newly identified two and replicated five GC risk loci, all of them with subtype-specific association. The gastric transcriptome data consisting of 361 corpus and 342 antrum mucosa samples revealed that an upregulated expression of MUC1, ANKRD50, PTGER4, and PSCA are plausible GC-pathomechanisms at four GWAS loci. At another risk locus, we found that the blood-group 0 exerts protective effects for non-cardia and diffuse GC, while blood-group A increases risk for both GC subtypes. Furthermore, our GWAS on cardia GC and OAC/BO (10,279 patients, 16,527 controls) showed that both cancer entities share genetic aetiology at the polygenic level and identified two new risk loci on the single-marker level. INTERPRETATION: Our findings show that the pathophysiology of GC is genetically heterogenous according to location and histopathology. Moreover, our findings point to common molecular mechanisms underlying cardia GC and OAC/BO. FUNDING: German Research Foundation (DFG).

GWAS meta-analysis of 16 790 patients with Barrett's oesophagus and oesophageal adenocarcinoma identifies 16 novel genetic risk loci and provides insights into disease aetiology beyond the single marker level.

OBJECTIVE: Oesophageal cancer (EC) is the sixth leading cause of cancer-related deaths. Oesophageal adenocarcinoma (EA), with Barrett's oesophagus (BE) as a precursor lesion, is the most prevalent EC subtype in the Western world. This study aims to contribute to better understand the genetic causes of BE/EA by leveraging genome wide association studies (GWAS), genetic correlation analyses and polygenic risk modelling. DESIGN: We combined data from previous GWAS with new cohorts, increasing the sample size to 16 790 BE/EA cases and 32 476 controls. We also carried out a transcriptome wide association study (TWAS) using expression data from disease-relevant tissues to identify BE/EA candidate genes. To investigate the relationship with reported BE/EA risk factors, a linkage disequilibrium score regression (LDSR) analysis was performed. BE/EA risk models were developed combining clinical/lifestyle risk factors with polygenic risk scores (PRS) derived from the GWAS meta-analysis. RESULTS: The GWAS meta-analysis identified 27 BE and/or EA risk loci, 11 of which were novel. The TWAS identified promising BE/EA candidate genes at seven GWAS loci and at five additional risk loci. The LDSR analysis led to the identification of novel genetic correlations and pointed to differences in BE and EA aetiology. Gastro-oesophageal reflux disease appeared to contribute stronger to the metaplastic BE transformation than to EA development. Finally, combining PRS with BE/EA risk factors improved the performance of the risk models. CONCLUSION: Our findings provide further insights into BE/EA aetiology and its relationship to risk factors. The results lay the foundation for future follow-up studies to identify underlying disease mechanisms and improving risk prediction.

Elucidation of the genetic causes of bicuspid aortic valve disease.

AIMS: The present study aims to characterize the genetic risk architecture of bicuspid aortic valve (BAV) disease, the most common congenital heart defect. METHODS AND RESULTS: We carried out a genome-wide association study (GWAS) including 2236 BAV patients and 11 604 controls. This led to the identification of a new risk locus for BAV on chromosome 3q29. The single nucleotide polymorphism rs2550262 was genome-wide significant BAV associated (P = 3.49 × 10-08) and was replicated in an independent case-control sample. The risk locus encodes a deleterious missense variant in MUC4 (p.Ala4821Ser), a gene that is involved in epithelial-to-mesenchymal transformation. Mechanistical studies in zebrafish revealed that loss of Muc4 led to a delay in cardiac valvular development suggesting that loss of MUC4 may also play a role in aortic valve malformation. The GWAS also confirmed previously reported BAV risk loci at PALMD (P = 3.97 × 10-16), GATA4 (P = 1.61 × 10-09), and TEX41 (P = 7.68 × 10-04). In addition, the genetic BAV architecture was examined beyond the single-marker level revealing that a substantial fraction of BAV heritability is polygenic and ∼20% of the observed heritability can be explained by our GWAS data. Furthermore, we used the largest human single-cell atlas for foetal gene expression and show that the transcriptome profile in endothelial cells is a major source contributing to BAV pathology. CONCLUSION: Our study provides a deeper understanding of the genetic risk architecture of BAV formation on the single marker and polygenic level.

A genome-wide association study with tissue transcriptomics identifies genetic drivers for classic bladder exstrophy.

Classic bladder exstrophy represents the most severe end of all human congenital anomalies of the kidney and urinary tract and is associated with bladder cancer susceptibility. Previous genetic studies identified one locus to be involved in classic bladder exstrophy, but were limited to a restrict number of cohort. Here we show the largest classic bladder exstrophy genome-wide association analysis to date where we identify eight genome-wide significant loci, seven of which are novel. In these regions reside ten coding and four non-coding genes. Among the coding genes is EFNA1, strongly expressed in mouse embryonic genital tubercle, urethra, and primitive bladder. Re-sequence of EFNA1 in the investigated classic bladder exstrophy cohort of our study displays an enrichment of rare protein altering variants. We show that all coding genes are expressed and/or significantly regulated in both mouse and human embryonic developmental bladder stages. Furthermore, nine of the coding genes residing in the regions of genome-wide significance are differentially expressed in bladder cancers. Our data suggest genetic drivers for classic bladder exstrophy, as well as a possible role for these drivers to relevant bladder cancer susceptibility.

Genome-wide association study in patients with posterior urethral valves.

Congenital lower urinary tract obstructions (LUTO) are most often caused by posterior urethral valves (PUV), a male limited anatomical obstruction of the urethra affecting 1 in 4,000 male live births. Little is known about the genetic background of PUV. Here, we report the largest genome-wide association study (GWAS) for PUV in 4 cohorts of patients and controls. The final meta-analysis included 756 patients and 4,823 ethnicity matched controls and comprised 5,754,208 variants that were genotyped or imputed and passed quality control in all 4 cohorts. No genome-wide significant locus was identified, but 33 variants showed suggestive significance (P < 1 × 10-5). When considering only loci with multiple variants residing within < 10 kB of each other showing suggestive significance and with the same effect direction in all 4 cohorts, 3 loci comprising a total of 9 variants remained. These loci resided on chromosomes 13, 16, and 20. The present GWAS and meta-analysis is the largest genetic study on PUV performed to date. The fact that no genome-wide significant locus was identified, can be explained by lack of power or may indicate that common variants do not play a major role in the etiology of PUV. Nevertheless, future studies are warranted to replicate and validate the 3 loci that yielded suggestive associations.

First genome-wide association study of esophageal atresia identifies three genetic risk loci at CTNNA3, FOXF1/FOXC2/FOXL1, and HNF1B.

Esophageal atresia with or without tracheoesophageal fistula (EA/TEF) is the most common congenital malformation of the upper digestive tract. This study represents the first genome-wide association study (GWAS) to identify risk loci for EA/TEF. We used a European case-control sample comprising 764 EA/TEF patients and 5,778 controls and observed genome-wide significant associations at three loci. On chromosome 10q21 within the gene CTNNA3 (p = 2.11 × 10-8; odds ratio [OR] = 3.94; 95% confidence interval [CI], 3.10-5.00), on chromosome 16q24 next to the FOX gene cluster (p = 2.25 × 10-10; OR = 1.47; 95% CI, 1.38-1.55) and on chromosome 17q12 next to the gene HNF1B (p = 3.35 × 10-16; OR = 1.75; 95% CI, 1.64-1.87). We next carried out an esophageal/tracheal transcriptome profiling in rat embryos at four selected embryonic time points. Based on these data and on already published data, the implicated genes at all three GWAS loci are promising candidates for EA/TEF development. We also analyzed the genetic EA/TEF architecture beyond the single marker level, which revealed an estimated single-nucleotide polymorphism (SNP)-based heritability of around 37% ± 14% standard deviation. In addition, we examined the polygenicity of EA/TEF and found that EA/TEF is less polygenic than other complex genetic diseases. In conclusion, the results of our study contribute to a better understanding on the underlying genetic architecture of ET/TEF with the identification of three risk loci and candidate genes.

Allele-specific transcription factor binding in a cellular model of orofacial clefting.

Non-syndromic cleft lip with/without cleft palate (nsCL/P) is a frequent congenital malformation with multifactorial etiology. While recent genome-wide association studies (GWAS) have identified several nsCL/P risk loci, the functional effects of the associated non-coding variants are largely unknown. Furthermore, additional risk loci remain undetected due to lack of power. As genetic variants might alter binding of transcription factors (TF), we here hypothesized that the integration of data from TF binding sites, expression analyses and nsCL/P GWAS might help to (i) identify functionally relevant variants at GWAS loci, and (ii) highlight novel risk variants that have been previously undetected. Analysing the craniofacial TF TFAP2A in human embryonic palatal mesenchyme (HEPM) cells, we identified 2845 TFAP2A ChIP-seq peaks, several of which were located near nsCL/P candidate genes (e.g. MSX1 and SPRY2). Comparison with independent data suggest that 802 of them might be specific to craniofacial development, and genes near these peaks are enriched in processes relevant to nsCL/P. Integration with nsCL/P GWAS data, however, did not show robust evidence for co-localization of common nsCL/P risk variants with TFAP2A ChIP-seq peaks. This data set represents a new resource for the analyses of craniofacial processes, and similar approaches with additional cell lines and TFs could be applied to generate further insights into nsCL/P etiology.

Human exome and mouse embryonic expression data implicate ZFHX3, TRPS1, and CHD7 in human esophageal atresia.

INTRODUCTION: Esophageal atresia with or without tracheoesophageal fistula (EA/TEF) occurs approximately 1 in 3.500 live births representing the most common malformation of the upper digestive tract. Only half a century ago, EA/TEF was fatal among affected newborns suggesting that the steady birth prevalence might in parts be due to mutational de novo events in genes involved in foregut development. METHODS: To identify mutational de novo events in EA/TEF patients, we surveyed the exome of 30 case-parent trios. Identified and confirmed de novo variants were prioritized using in silico prediction tools. To investigate the embryonic role of genes harboring prioritized de novo variants we performed targeted analysis of mouse transcriptome data of esophageal tissue obtained at the embryonic day (E) E8.5, E12.5, and postnatal. RESULTS: In total we prioritized 14 novel de novo variants in 14 different genes (APOL2, EEF1D, CHD7, FANCB, GGT6, KIAA0556, NFX1, NPR2, PIGC, SLC5A2, TANC2, TRPS1, UBA3, and ZFHX3) and eight rare de novo variants in eight additional genes (CELSR1, CLP1, GPR133, HPS3, MTA3, PLEC, STAB1, and PPIP5K2). Through personal communication during the project, we identified an additional EA/TEF case-parent trio with a rare de novo variant in ZFHX3. In silico prediction analysis of the identified variants and comparative analysis of mouse transcriptome data of esophageal tissue obtained at E8.5, E12.5, and postnatal prioritized CHD7, TRPS1, and ZFHX3 as EA/TEF candidate genes. Re-sequencing of ZFHX3 in additional 192 EA/TEF patients did not identify further putative EA/TEF-associated variants. CONCLUSION: Our study suggests that rare mutational de novo events in genes involved in foregut development contribute to the development of EA/TEF.

Toll-like Receptor Signaling Rewires Macrophage Metabolism and Promotes Histone Acetylation via ATP-Citrate Lyase.

Toll-like receptor (TLR) activation induces inflammatory responses in macrophages by activating temporally defined transcriptional cascades. Whether concurrent changes in the cellular metabolism that occur upon TLR activation influence the quality of the transcriptional responses remains unknown. Here, we investigated how macrophages adopt their metabolism early after activation to regulate TLR-inducible gene induction. Shortly after TLR4 activation, macrophages increased glycolysis and tricarboxylic acid (TCA) cycle volume. Metabolic tracing studies revealed that TLR signaling redirected metabolic fluxes to generate acetyl-Coenzyme A (CoA) from glucose resulting in augmented histone acetylation. Signaling through the adaptor proteins MyD88 and TRIF resulted in activation of ATP-citrate lyase, which in turn facilitated the induction of distinct LPS-inducible gene sets. We postulate that metabolic licensing of histone acetylation provides another layer of control that serves to fine-tune transcriptional responses downstream of TLR activation. Our work highlights the potential of targeting the metabolic-epigenetic axis in inflammatory settings.

Evidence for PTGER4, PSCA, and MBOAT7 as risk genes for gastric cancer on the genome and transcriptome level.

Genetic associations between variants on chromosome 5p13 and 8q24 and gastric cancer (GC) have been previously reported in the Asian population. We aimed to replicate these findings and to characterize the associations at the genome and transcriptome level. We performed a fine-mapping association study in 1926 GC patients and 2012 controls of European descent using high dense SNP marker sets on both chromosomal regions. Next, we performed expression quantitative trait locus (eQTL) analyses using gastric transcriptome data from 143 individuals focusing on the GC associated variants. On chromosome 5p13 the strongest association was observed at rs6872282 (P = 2.53 × 10-04 ) and on chromosome 8q24 at rs2585176 (P = 1.09 × 10-09 ). On chromosome 5p13 we found cis-eQTL effects with an upregulation of PTGER4 expression in GC risk allele carrier (P = 9.27 × 10-11 ). On chromosome 8q24 we observed cis-eQTL effects with an upregulation of PSCA expression in GC risk allele carrier (P = 2.17 × 10-47 ). In addition, we found trans-eQTL effects for the same variants on 8q24 with a downregulation of MBOAT7 expression in GC risk allele carrier (P = 3.11 × 10-09 ). In summary, we confirmed and refined the previously reported GC associations at both chromosomal regions. Our data point to shared etiological factors between Asians and Europeans. Furthermore, our data imply an upregulated expression of PTGER4 and PSCA as well as a downregulated expression of MBOAT7 in gastric tissue as risk-conferring GC pathomechanisms.

Exome sequencing in syndromic brain malformations identifies novel mutations in ACTB, and SLC9A6, and suggests BAZ1A as a new candidate gene.

BACKGROUND: Syndromic brain malformations comprise a large group of anomalies with a birth prevalence of about 1 in 1,000 live births. Their etiological factors remain largely unknown. To identify causative mutations, we used whole-exome sequencing (WES) in aborted fetuses and children with syndromic brain malformations in which chromosomal microarray analysis was previously unremarkable. METHODS: WES analysis was applied in eight case-parent trios, six aborted fetuses, and two children. RESULTS: WES identified a novel de novo mutation (p.Gly268Arg) in ACTB (Baraitser-Winter syndrome-1), a homozygous stop mutation (p.R2442*) in ASPM (primary microcephaly type 5), and a novel hemizygous X-chromosomal mutation (p.I250V) in SLC9A6 (X-linked syndromic mentaly retardation, Christianson type). Furthermore, WES identified a de novo mutation (p.Arg1093Gln) in BAZ1A. This mutation was previously reported in only one allele in 121.362 alleles tested (dbSNP build 147). BAZ1A has been associated with neurodevelopmental impairment and dysregulation of several pathways including vitamin D metabolism. Here, serum vitamin-D (25-(OH)D) levels were insufficient and gene expression comparison between the child and her parents identified 27 differentially expressed genes. Of note, 10 out of these 27 genes are associated to cytoskeleton, integrin and synaptic related pathways, pinpointing to the relevance of BAZ1A in neural development. In situ hybridization in mouse embryos between E10.5 and E13.5 detected Baz1a expression in the central and peripheral nervous system. CONCLUSION: In syndromic brain malformations, WES is likely to identify causative mutations when chromosomal microarray analysis is unremarkable. Our findings suggest BAZ1A as a possible new candidate gene.