Identification of the Molecular Components of Enhancer-Mediated Gene Expression Variation in Multiple Tissues Regulating Blood Pressure.

BACKGROUND: Inter-individual variation in blood pressure (BP) arises in part from sequence variants within enhancers modulating the expression of causal genes. We propose that these genes, active in tissues relevant to BP physiology, can be identified from tissue-level epigenomic data and genotypes of BP-phenotyped individuals. METHODS: We used chromatin accessibility data from the heart, adrenal, kidney, and artery to identify cis-regulatory elements (CREs) in these tissues and estimate the impact of common human single-nucleotide variants within these CREs on gene expression, using machine learning methods. To identify causal genes, we performed a gene-wise association test. We conducted analyses in 2 separate large-scale cohorts: 77 822 individuals from the Genetic Epidemiology Research on Adult Health and Aging and 315 270 individuals from the UK Biobank. RESULTS: We identified 309, 259, 331, and 367 genes (false discovery rate <0.05) for diastolic BP and 191, 184, 204, and 204 genes for systolic BP in the artery, kidney, heart, and adrenal, respectively, in Genetic Epidemiology Research on Adult Health and Aging; 50% to 70% of these genes were replicated in the UK Biobank, significantly higher than the 12% to 15% expected by chance (P<0.0001). These results enabled tissue expression prediction of these 988 to 2875 putative BP genes in individuals of both cohorts to construct an expression polygenic score. This score explained ≈27% of the reported single-nucleotide variant heritability, substantially higher than expected from prior studies. CONCLUSIONS: Our work demonstrates the power of tissue-restricted comprehensive CRE analysis, followed by CRE-based expression prediction, for understanding BP regulation in relevant tissues and provides dual-modality supporting evidence, CRE and expression, for the causality genes.

RET enhancer haplotype-dependent remodeling of the human fetal gut development program.

Hirschsprung disease (HSCR) is associated with deficiency of the receptor tyrosine kinase RET, resulting in loss of cells of the enteric nervous system (ENS) during fetal gut development. The major contribution to HSCR risk is from common sequence variants in RET enhancers with additional risk from rare coding variants in many genes. Here, we demonstrate that these RET enhancer variants specifically alter the human fetal gut development program through significant decreases in gene expression of RET, members of the RET-EDNRB gene regulatory network (GRN), other HSCR genes, with an altered transcriptome of 2,382 differentially expressed genes across diverse neuronal and mesenchymal functions. A parsimonious hypothesis for these results is that beyond RET's direct effect on its GRN, it also has a major role in enteric neural crest-derived cell (ENCDC) precursor proliferation, its deficiency reducing ENCDCs with relative expansion of non-ENCDC cells. Thus, genes reducing RET proliferative activity can potentially cause HSCR. One such class is the 23 RET-dependent transcription factors enriched in early gut development. We show that their knockdown in human neuroblastoma SK-N-SH cells reduces RET and/or EDNRB gene expression, expanding the RET-EDNRB GRN. The human embryos we studied had major remodeling of the gut transcriptome but were unlikely to have had HSCR: thus, genetic or epigenetic changes in addition to those in RET are required for aganglionosis.

Tissue-specific and tissue-agnostic effects of genome sequence variation modulating blood pressure.

Genome-wide association studies (GWASs) have identified numerous variants associated with polygenic traits and diseases. However, with few exceptions, a mechanistic understanding of which variants affect which genes in which tissues to modulate trait variation is lacking. Here, we present genomic analyses to explain trait heritability of blood pressure (BP) through the genetics of transcriptional regulation using GWASs, multiomics data from different tissues, and machine learning approaches. Approximately 500,000 predicted regulatory variants across four tissues explain 33.4% of variant heritability: 2.5%, 5.3%, 7.7%, and 11.8% for kidney-, adrenal-, heart-, and artery-specific variants, respectively. Variation in the enhancers involved shows greater tissue specificity than in the genes they regulate, suggesting that gene regulatory networks perturbed by enhancer variants in a tissue relevant to a phenotype are the major source of interindividual variation in BP. Thus, our study provides an approach to scan human tissue and cell types for their physiological contribution to any trait.

Wnt signaling enhances macrophage responses to IL-4 and promotes resolution of atherosclerosis.

Atherosclerosis is a disease of chronic inflammation. We investigated the roles of the cytokines IL-4 and IL-13, the classical activators of STAT6, in the resolution of atherosclerosis inflammation. Using Il4-/-Il13-/- mice, resolution was impaired, and in control mice, in both progressing and resolving plaques, levels of IL-4 were stably low and IL-13 was undetectable. This suggested that IL-4 is required for atherosclerosis resolution, but collaborates with other factors. We had observed increased Wnt signaling in macrophages in resolving plaques, and human genetic data from others showed that a loss-of-function Wnt mutation was associated with premature atherosclerosis. We now find an inverse association between activation of Wnt signaling and disease severity in mice and humans. Wnt enhanced the expression of inflammation resolving factors after treatment with plaque-relevant low concentrations of IL-4. Mechanistically, activation of the Wnt pathway following lipid lowering potentiates IL-4 responsiveness in macrophages via a PGE2/STAT3 axis.

In search of a genetic explanation for LDLc variability in an FH family: common SNPs and a rare mutation in MTTP explain only part of LDL variability in an FH family.

We previously identified a highly consanguineous familial hypercholesterolemia (FH) family demonstrating segregation of the JD Bari mutation in the LDL receptor as well as a putative cholesterol-lowering trait. We aimed to identify genes related to the latter effect. LDL cholesterol (LDLc) values were normalized for FH affectation status, age, and gender. Using genome-wide SNP data, we examined whether known SNPs gleaned from a genome-wide association study could explain the variation observed in LDLc. Four individuals with markedly reduced LDL levels underwent whole exome sequencing. After prioritizing all potential mutations, we identified the most promising candidate genes and tested them for segregation with the lowering trait. We transfected a plasmid carrying the top candidate mutation, microsomal triglyceride transfer protein (MTTP) R634C, into COS-7 cells to test enzymatic activity. The SNP score explained 3% of the observed variability. MTTP R634C showed reduced activity (49.1 nmol/ml) compared with the WT allele (185.8 nmol/ml) (P = 0.0012) and was marginally associated with reduced LDLc in FH patients (P = 0.05). Phenotypic variability in a FH pedigree can only partially be explained by a combination of common SNPs and a rare mutation and a rare variant in the MTTP gene. LDLc variability in FH patients may have nongenetic causes.