LINE1 and PRC2 control nucleolar organization and repression of the 8C state in human ESCs.

The mechanisms that ensure developmental progression in the early human embryo remain largely unknown. Here, we show that the family of long interspersed nuclear element 1 (LINE1) transposons prevents the reversion of naive human embryonic stem cells (hESCs) to 8-cell-like cells (8CLCs). LINE1 RNA contributes to maintenance of H3K27me3 levels, particularly at chromosome 19 (Chr19). Chr19 is enriched for key 8C regulators, H3K27me3, and genes derepressed upon LINE1 knockdown or PRC2 inhibition. Moreover, Chr19 is strongly associated with the nucleolus in hESCs but less in 8CLCs. Direct inhibition of PRC2 activity induces the 8C program and leads to a relocalization of Chr19 away from the nucleolus. LINE1 KD or PRC2 inhibition induces nucleolar stress, and disruption of nucleolar architecture is sufficient to de-repress the 8C program. These results indicate that LINE1 RNA and PRC2 maintain H3K27me3-mediated gene repression and 3D nuclear organization to prevent developmental reversion of hESCs.

Autism-associated CHD8 controls reactive gliosis and neuroinflammation via remodeling chromatin in astrocytes.

Reactive changes of glial cells during neuroinflammation impact brain disorders and disease progression. Elucidating the mechanisms that control reactive gliosis may help us to understand brain pathophysiology and improve outcomes. Here, we report that adult ablation of autism spectrum disorder (ASD)-associated CHD8 in astrocytes attenuates reactive gliosis via remodeling chromatin accessibility, changing gene expression. Conditional Chd8 deletion in astrocytes, but not microglia, suppresses reactive gliosis by impeding astrocyte proliferation and morphological elaboration. Astrocyte Chd8 ablation alleviates lipopolysaccharide-induced neuroinflammation and septic-associated hypothermia in mice. Astrocytic CHD8 plays an important role in neuroinflammation by altering the chromatin landscape, regulating metabolic and lipid-associated pathways, and astrocyte-microglia crosstalk. Moreover, we show that reactive gliosis can be directly mitigated in vivo using an adeno-associated virus (AAV)-mediated Chd8 gene editing strategy. These findings uncover a role of ASD-associated CHD8 in the adult brain, which may warrant future exploration of targeting chromatin remodelers in reactive gliosis and neuroinflammation in injury and neurological diseases.

Transposable elements are associated with the variable response to influenza infection.

Influenza A virus (IAV) infections are frequent every year and result in a range of disease severity. Here, we wanted to explore the potential contribution of transposable elements (TEs) to the variable human immune response. Transcriptome profiling in monocyte-derived macrophages from 39 individuals following IAV infection revealed significant inter-individual variation in viral load post-infection. Using transposase-accessible chromatin using sequencing (ATAC-seq), we identified a set of TE families with either enhanced or reduced accessibility upon infection. Of the enhanced families, 15 showed high variability between individuals and had distinct epigenetic profiles. Motif analysis showed an association with known immune regulators (e.g., BATFs, FOSs/JUNs, IRFs, STATs, NFkBs, NFYs, and RELs) in stably enriched families and with other factors in variable families, including KRAB-ZNFs. We showed that TEs and host factors regulating TEs were predictive of viral load post-infection. Our findings shed light on the role TEs and KRAB-ZNFs may play in inter-individual variation in immunity.

Capturing sex-specific and hypofertility-linked effects of assisted reproductive technologies on the cord blood DNA methylome.

BACKGROUND: Children conceived through assisted reproduction are at an increased risk for growth and genomic imprinting disorders, often linked to DNA methylation defects. It has been suggested that assisted reproductive technology (ART) and underlying parental infertility can induce epigenetic instability, specifically interfering with DNA methylation reprogramming events during germ cell and preimplantation development. To date, human studies exploring the association between ART and DNA methylation defects have reported inconsistent or inconclusive results, likely due to population heterogeneity and the use of technologies with limited coverage of the epigenome. In our study, we explored the epigenetic risk of ART by comprehensively profiling the DNA methylome of 73 human cord blood samples of singleton pregnancies (n = 36 control group, n = 37 ART/hypofertile group) from a human prospective longitudinal birth cohort, the 3D (Design, Develop, Discover) Study, using a high-resolution sequencing-based custom capture panel that examines over 2.4 million autosomal CpGs in the genome. RESULTS: We identified evidence of sex-specific effects of ART/hypofertility on cord blood DNA methylation patterns. Our genome-wide analyses identified ~ 46% more CpGs affected by ART/hypofertility in female than in male infant cord blood. We performed a detailed analysis of three imprinted genes which have been associated with altered DNA methylation following ART (KCNQ1OT1, H19/IGF2 and GNAS) and found that female infant cord blood was associated with DNA hypomethylation. When compared to less invasive procedures such as intrauterine insemination, more invasive ARTs (in vitro fertilization, intracytoplasmic sperm injection, embryo culture) resulted in more marked and distinct effects on the cord blood DNA methylome. In the in vitro group, we found a close to fourfold higher proportion of significantly enriched Gene Ontology terms involved in development than in the in vivo group. CONCLUSIONS: Our study highlights the ability of a sensitive, targeted, sequencing-based approach to uncover DNA methylation perturbations in cord blood associated with hypofertility and ART and influenced by offspring sex and ART technique invasiveness.

Targeting inflammatory monocytes by immune-modifying nanoparticles prevents acute kidney allograft rejection.

Inflammatory monocytes are a major component of the cellular infiltrate in acutely rejecting human kidney allografts. Since immune-modifying nanoparticles (IMPs) bind to circulating inflammatory monocytes via the specific scavenger receptor MARCO, causing diversion to the spleen and subsequent apoptosis, we investigated the therapeutic potential of negatively charged, 500-nm diameter polystyrene IMPs to prevent kidney allograft rejection. Kidney transplants were performed from BALB/c (H2d) to C57BL/6 (H2b) mice in two groups: controls (allo) and allo mice infused with IMPs. Groups were studied for 14 (acute rejection) or 100 (chronic rejection) days. Allo mice receiving IMPs exhibited superior survival and markedly less acute rejection, with better kidney function, less tubulitis, and diminished inflammatory cell density, cytokine and cytotoxic molecule expression in the allograft and lower titers of donor-specific IgG2c antibody in serum at day 14, as compared to allo mice. Cells isolated from kidneys from allo mice receiving IMPs showed reduced Ly6Chi monocytes, CD11b+ cells and NKT+ cells compared to allo mice. IMPs predominantly bound CD11b+ cells in the bloodstream and CD11b+ and CD11c-B220+ marginal zone B cells in the spleen. In the spleen, IMPs were found predominantly in red pulp, colocalized with MARCO and expression of cleaved caspase-3. At day 100, allo mice receiving IMPs exhibited reduced macrophage M1 responses but were not protected from chronic rejection. IMPs afforded significant protection from acute rejection, inhibiting both innate and adaptive alloimmunity. Thus, our current experimental findings, coupled with our earlier demonstration of IMP-induced protection in kidney ischemia-reperfusion injury, identify IMPs as a potential induction agent in kidney transplantation.

Differentially methylated CpGs in response to growth hormone administration in children with idiopathic short stature.

BACKGROUND: Recombinant human growth hormone (rhGH) has shown a great growth-promoting potential in children with idiopathic short stature (ISS). However, the response to rhGH differs across individuals, largely due to genetic and epigenetic heterogeneity. Since epigenetic marks on the methylome can be dynamically influenced by GH, we performed a comprehensive pharmacoepigenomics analysis of DNA methylation changes associated with long-term rhGH administration in children with ISS. RESULTS: We measured DNA methylation profiles before and after GH treatment (with a duration of ~ 18 months in average) on 47 healthy children using customized methylC-seq capture sequencing. Their changes were compared and associated with changes in plasma IGF1 by adjusting sex, age, treatment duration and estimated blood proportions. We observed a considerable inter-individual heterogeneity of DNA methylation changes responding to GH treatment. We identified 267 response-associated differentially methylated cytosines (DMCs) that were enriched in promoter regions, CpG islands and blood cell-type-specific regulatory elements. Furthermore, the genes associated with these DMCs were enriched in the biology process of "cell development," "neuron differentiation" and "developmental growth," and in the TGF-beta signaling pathway, PPAR Alpha pathway, endoderm differentiation pathway, adipocytokine signaling pathway as well as PI3K-Akt signaling pathway, and cAMP signaling pathway. CONCLUSION: Our study provides a first insight in DNA methylation changes associated with rhGH administration, which may help understand mechanisms of epigenetic regulation on GH-responsive genes.

Application of ATAC-Seq for genome-wide analysis of the chromatin state at single myofiber resolution.

Myofibers are the main components of skeletal muscle, which is the largest tissue in the body. Myofibers are highly adaptive and can be altered under different biological and disease conditions. Therefore, transcriptional and epigenetic studies on myofibers are crucial to discover how chromatin alterations occur in the skeletal muscle under different conditions. However, due to the heterogenous nature of skeletal muscle, studying myofibers in isolation proves to be a challenging task. Single-cell sequencing has permitted the study of the epigenome of isolated myonuclei. While this provides sequencing with high dimensionality, the sequencing depth is lacking, which makes comparisons between different biological conditions difficult. Here, we report the first implementation of single myofiber ATAC-Seq, which allows for the sequencing of an individual myofiber at a depth sufficient for peak calling and for comparative analysis of chromatin accessibility under various physiological and disease conditions. Application of this technique revealed significant differences in chromatin accessibility between resting and regenerating myofibers, as well as between myofibers from a mouse model of Duchenne Muscular Dystrophy (mdx) and wild-type (WT) counterparts. This technique can lead to a wide application in the identification of chromatin regulatory elements and epigenetic mechanisms in muscle fibers during development and in muscle-wasting diseases.

Whole-genome sequencing of H3K4me3 and DNA methylation in human sperm reveals regions of overlap linked to fertility and development.

The paternal environment has been linked to infertility and negative outcomes. Such effects may be transmitted via sperm through histone modifications. To date, in-depth profiling of the sperm chromatin in men has been limited. Here, we use deep sequencing to characterize the sperm profiles of histone H3 lysine 4 tri-methylation (H3K4me3) and DNA methylation in a representative reference population of 37 men. Our analysis reveals that H3K4me3 is localized throughout the genome and at genes for fertility and development. Remarkably, enrichment is also found at regions that escape epigenetic reprogramming in primordial germ cells, embryonic enhancers, and short-interspersed nuclear elements (SINEs). There is significant overlap in H3K4me3 and DNA methylation throughout the genome, suggesting a potential interplay between these marks previously reported to be mutually exclusive in sperm. Comparisons made between H3K4me3 marked regions in sperm and the embryonic transcriptome suggest an influence of paternal chromatin on embryonic gene expression.

Non-CG methylation and multiple histone profiles associate child abuse with immune and small GTPase dysregulation.

Early-life adversity (ELA) is a major predictor of psychopathology, and is thought to increase lifetime risk by epigenetically regulating the genome. Here, focusing on the lateral amygdala, a major brain site for emotional homeostasis, we describe molecular cross-talk among multiple mechanisms of genomic regulation, including 6 histone marks and DNA methylation, and the transcriptome, in subjects with a history of ELA and controls. In the healthy brain tissue, we first uncover interactions between different histone marks and non-CG methylation in the CAC context. Additionally, we find that ELA associates with methylomic changes that are as frequent in the CAC as in the canonical CG context, while these two forms of plasticity occur in sharply distinct genomic regions, features, and chromatin states. Combining these multiple data indicates that immune-related and small GTPase signaling pathways are most consistently impaired in the amygdala of ELA individuals. Overall, this work provides insights into genomic brain regulation as a function of early-life experience.

High-resolution analyses of human sperm dynamic methylome reveal thousands of novel age-related epigenetic alterations.

BACKGROUND: Children of aged fathers are at a higher risk of developing mental disorders. Alterations in sperm DNA methylation have been implicated as a potential cause. However, age-dependent modifications of the germ cells' epigenome remain poorly understood. Our objective was to assess the DNA methylation profile of human spermatozoa during aging. RESULTS: We used a high throughput, customized methylC-capture sequencing (MCC-seq) approach to characterize the dynamic DNA methylation in spermatozoa from 94 fertile and infertile men, who were categorized as young, 48 men between 18-38 years or old 46 men between 46-71 years. We identified more than 150,000 age-related CpG sites that are significantly differentially methylated among 2.65 million CpG sites covered. We conducted machine learning using our dataset to predict the methylation age of subjects; the age prediction accuracy based on our assay provided a more accurate prediction than that using the 450 K chip approach. In addition, we found that there are more hypermethylated (62%) than hypomethylated (38%) CpG sites in sperm of aged men, corresponding to 798 of total differential methylated regions (DMRs), of which 483 are hypermethylated regions (HyperDMR), and 315 hypomethylated regions (HypoDMR). Moreover, the distribution of age-related hyper- and hypomethylated CpGs in sperm is not random; the CpG sites that were hypermethylated with advanced age were frequently located in the distal region to genes, whereas hypomethylated sites were near to gene transcription start sites (TSS). We identified a high density of age-associated CpG changes in chromosomes 4 and 16, particularly HyperDMRs with localized clusters, the chr4 DMR cluster overlaps PGC1α locus, a protein involved in metabolic aging and the chr16 DMR cluster overlaps RBFOX1 locus, a gene implicated in neurodevelopmental disease. Gene ontology analysis revealed that the most affected genes by age were associated with development, neuron projection, differentiation and recognition, and behaviour, suggesting a potential link to the higher risk of neurodevelopmental disorders in children of aged fathers. CONCLUSION: We identified thousands of age-related and sperm-specific epigenetic alterations. These findings provide novel insight in understanding human sperm DNA methylation dynamics during paternal aging, and the subsequently affected genes potentially related to diseases in offspring.

Eosinophil microRNAs Play a Regulatory Role in Allergic Diseases Included in the Atopic March.

(1) Background: The atopic march is defined by the increased prevalence of allergic diseases after atopic dermatitis onset. In fact, atopic dermatitis is believed to play an important role in allergen sensitization via the damaged skin barrier, leading to allergic diseases such as allergic asthma and allergic rhinitis. The eosinophil, a pro-inflammatory cell that contributes to epithelial damage, is one of the various cells recruited in the inflammatory reactions characterizing these diseases. Few studies were conducted on the transcriptome of this cell type and even less on their specific microRNA (miRNA) profile, which could modulate pathogenesis of allergic diseases and clinical manifestations post-transcriptionally. Actually, their implication in allergic diseases is not fully understood, but they are believed to play a role in inflammation-related patterns and epithelial cell proliferation. (2) Methods: Next-generation sequencing was performed on RNA samples from eosinophils of individuals with atopic dermatitis, atopy, allergic rhinitis and asthma to obtain differential counts of primary miRNA (pri-miRNA); these were also analyzed for asthma-related phenotypes such as forced expiratory volume in one second (FEV1), immunoglobulin E (IgE) and provocative concentration of methacholine inducing a 20% fall in forced expiratory volume in 1 s (PC20) levels, as well as FEV1 to forced vital capacity (FEV1/FVC) ratio. (3) Results: Eighteen miRNAs from eosinophils were identified to be significantly different between affected individuals and unaffected ones. Based on counts from these miRNAs, individuals were then clustered into groups using Ward's method on Euclidian distances. Groups were found to be explained by asthma diagnosis, familial history of respiratory diseases and allergic rhinitis as well as neutrophil counts. (4) Conclusions: The 18 differential miRNA counts for the studying phenotypes allow a better understanding of the epigenetic mechanisms underlying the development of the allergic diseases included in the atopic march.

Paired rRNA-depleted and polyA-selected RNA sequencing data and supporting multi-omics data from human T cells.

Both poly(A) enrichment and ribosomal RNA depletion are commonly used for RNA sequencing. Either has its advantages and disadvantages that may lead to biases in the downstream analyses. To better access these effects, we carried out both ribosomal RNA-depleted and poly(A)-selected RNA-seq for CD4+ T naive cells isolated from 40 healthy individuals from the Blueprint Project. For these 40 individuals, the genomic and epigenetic data were also available. This dataset offers a unique opportunity to understand how library construction influences differential gene expression, alternative splicing and molecular QTL (quantitative loci) analyses for human primary cells.

Gut Microbial Metabolites Induce Donor-Specific Tolerance of Kidney Allografts through Induction of T Regulatory Cells by Short-Chain Fatty Acids.

BACKGROUND: Short-chain fatty acids derived from gut microbial fermentation of dietary fiber have been shown to suppress autoimmunity through mechanisms that include enhanced regulation by T regulatory cells (Tregs). METHODS: Using a murine kidney transplantation model, we examined the effects on alloimmunity of a high-fiber diet or supplementation with the short-chain fatty acid acetate. Kidney transplants were performed from BALB/c(H2d) to B6(H2b) mice as allografts in wild-type and recipient mice lacking the G protein-coupled receptor GPR43 (the metabolite-sensing receptor of acetate). Allograft mice received normal chow, a high-fiber diet, or normal chow supplemented with sodium acetate. We assessed rejection at days 14 (acute) and 100 (chronic), and used 16S rRNA sequencing to determine gut microbiota composition pretransplantation and post-transplantation. RESULTS: Wild-type mice fed normal chow exhibited dysbiosis after receiving a kidney allograft but not an isograft, despite the avoidance of antibiotics and immunosuppression for the latter. A high-fiber diet prevented dysbiosis in allograft recipients, who demonstrated prolonged survival and reduced evidence of rejection compared with mice fed normal chow. Allograft mice receiving supplemental sodium acetate exhibited similar protection from rejection, and subsequently demonstrated donor-specific tolerance. Depletion of CD25+ Tregs or absence of the short-chain fatty acid receptor GPR43 abolished this survival advantage. CONCLUSIONS: Manipulation of the microbiome by a high-fiber diet or supplementation with sodium acetate modified alloimmunity in a kidney transplant model, generating tolerance dependent on Tregs and GPR43. Diet-based therapy to induce changes in the gut microbiome can alter systemic alloimmunity in mice, in part through the production of short-chain fatty acids leading to Treg cell development, and merits study as a potential clinical strategy to facilitate transplant acceptance.

Personalized and graph genomes reveal missing signal in epigenomic data.

BACKGROUND: Epigenomic studies that use next generation sequencing experiments typically rely on the alignment of reads to a reference sequence. However, because of genetic diversity and the diploid nature of the human genome, we hypothesize that using a generic reference could lead to incorrectly mapped reads and bias downstream results. RESULTS: We show that accounting for genetic variation using a modified reference genome or a de novo assembled genome can alter histone H3K4me1 and H3K27ac ChIP-seq peak calls either by creating new personal peaks or by the loss of reference peaks. Using permissive cutoffs, modified reference genomes are found to alter approximately 1% of peak calls while de novo assembled genomes alter up to 5% of peaks. We also show statistically significant differences in the amount of reads observed in regions associated with the new, altered, and unchanged peaks. We report that short insertions and deletions (indels), followed by single nucleotide variants (SNVs), have the highest probability of modifying peak calls. We show that using a graph personalized genome represents a reasonable compromise between modified reference genomes and de novo assembled genomes. We demonstrate that altered peaks have a genomic distribution typical of other peaks. CONCLUSIONS: Analyzing epigenomic datasets with personalized and graph genomes allows the recovery of new peaks enriched for indels and SNVs. These altered peaks are more likely to differ between individuals and, as such, could be relevant in the study of various human phenotypes.

Dietary Fiber Protects against Diabetic Nephropathy through Short-Chain Fatty Acid-Mediated Activation of G Protein-Coupled Receptors GPR43 and GPR109A.

BACKGROUND: Studies have reported "dysbiotic" changes to gut microbiota, such as depletion of gut bacteria that produce short-chain fatty acids (SCFAs) through gut fermentation of fiber, in CKD and diabetes. Dietary fiber is associated with decreased inflammation and mortality in CKD, and SCFAs have been proposed to mediate this effect. METHODS: To explore dietary fiber's effect on development of experimental diabetic nephropathy, we used streptozotocin to induce diabetes in wild-type C57BL/6 and knockout mice lacking the genes encoding G protein-coupled receptors GPR43 or GPR109A. Diabetic mice were randomized to high-fiber, normal chow, or zero-fiber diets, or SCFAs in drinking water. We used proton nuclear magnetic resonance spectroscopy for metabolic profiling and 16S ribosomal RNA sequencing to assess the gut microbiome. RESULTS: Diabetic mice fed a high-fiber diet were significantly less likely to develop diabetic nephropathy, exhibiting less albuminuria, glomerular hypertrophy, podocyte injury, and interstitial fibrosis compared with diabetic controls fed normal chow or a zero-fiber diet. Fiber beneficially reshaped gut microbial ecology and improved dysbiosis, promoting expansion of SCFA-producing bacteria of the genera Prevotella and Bifidobacterium, which increased fecal and systemic SCFA concentrations. Fiber reduced expression of genes encoding inflammatory cytokines, chemokines, and fibrosis-promoting proteins in diabetic kidneys. SCFA-treated diabetic mice were protected from nephropathy, but not in the absence of GPR43 or GPR109A. In vitro, SCFAs modulated inflammation in renal tubular cells and podocytes under hyperglycemic conditions. CONCLUSIONS: Dietary fiber protects against diabetic nephropathy through modulation of the gut microbiota, enrichment of SCFA-producing bacteria, and increased SCFA production. GPR43 and GPR109A are critical to SCFA-mediated protection against this condition. Interventions targeting the gut microbiota warrant further investigation as a novel renoprotective therapy in diabetic nephropathy.

Computational Analysis of HLA-presentation of Non-synonymous Recipient Mismatches Indicates Effect on the Risk of Chronic Graft-vs.-Host Disease After Allogeneic HSCT.

Genetic mismatches in protein coding genes between allogeneic hematopoietic stem cell transplantation (allo-HSCT) recipient and donor can elicit an alloimmunity response via peptides presented by the recipient HLA receptors as minor histocompatibility antigens (mHAs). While the impact of individual mHAs on allo-HSCT outcome such as graft-vs.-host and graft-vs.-leukemia effects has been demonstrated, it is likely that established mHAs constitute only a small fraction of all immunogenic non-synonymous variants. In the present study, we have analyzed the genetic mismatching in 157 exome-sequenced sibling allo-HSCT pairs to evaluate the significance of polymorphic HLA class I associated peptides on clinical outcome. We applied computational mismatch estimation approaches based on experimentally verified HLA ligands available in public repositories, published mHAs, and predicted HLA-peptide affinites, and analyzed their associations with chronic graft-vs.-host disease (cGvHD) grades. We found that higher estimated recipient mismatching consistently increased the risk of severe cGvHD, suggesting that HLA-presented mismatching influences the likelihood of long-term complications in the patient. Furthermore, computational approaches focusing on estimation of HLA-presentation instead of all non-synonymous mismatches indiscriminately may be beneficial for analysis sensitivity and could help identify novel mHAs.

Customized MethylC-Capture Sequencing to Evaluate Variation in the Human Sperm DNA Methylome Representative of Altered Folate Metabolism.

BACKGROUND: The sperm DNA methylation landscape is unique and critical for offspring health. If gamete-derived DNA methylation escapes reprograming in early embryos, epigenetic defects in sperm may be transmitted to the next generation. Current techniques to assess sperm DNA methylation show bias toward CpG-dense regions and do not target areas of dynamic methylation, those predicted to be environmentally sensitive and tunable regulatory elements. OBJECTIVES: Our goal was to assess variation in human sperm DNA methylation and design a targeted capture panel to interrogate the human sperm methylome. METHODS: To characterize variation in sperm DNA methylation, we performed whole genome bisulfite sequencing (WGBS) on an equimolar pool of sperm DNA from a wide cross section of 30 men varying in age, fertility status, methylenetetrahydrofolate reductase (MTHFR) genotype, and exposures. With our targeted capture panel, in individual samples, we examined the effect of MTHFR genotype ([Formula: see text] 677CC, [Formula: see text] 677TT), as well as high-dose folic acid supplementation ([Formula: see text], per genotype, before and after supplementation). RESULTS: Through WGBS we discovered nearly 1 million CpGs possessing intermediate methylation levels (20-80%), termed dynamic sperm CpGs. These dynamic CpGs, along with 2 million commonly assessed CpGs, were used to customize a capture panel for targeted interrogation of the human sperm methylome and test its ability to detect effects of altered folate metabolism. As compared with MTHFR 677CC men, those with the 677TT genotype (50% decreased MTHFR activity) had both hyper- and hypomethylation in their sperm. High-dose folic acid supplement treatment exacerbated hypomethylation in MTHFR 677TT men compared with 677CC. In both cases, [Formula: see text] of altered methylation was found in dynamic sperm CpGs, uniquely measured by our assay. DISCUSSION: Our sperm panel allowed the discovery of differential methylation following conditions affecting folate metabolism in novel dynamic sperm CpGs. Improved ability to examine variation in sperm DNA methylation can facilitate comprehensive studies of environment-epigenome interactions. https://doi.org/10.1289/EHP4812.

Integrative analysis of vascular endothelial cell genomic features identifies AIDA as a coronary artery disease candidate gene.

BACKGROUND: Genome-wide association studies (GWAS) have identified hundreds of loci associated with coronary artery disease (CAD) and blood pressure (BP) or hypertension. Many of these loci are not linked to traditional risk factors, nor do they include obvious candidate genes, complicating their functional characterization. We hypothesize that many GWAS loci associated with vascular diseases modulate endothelial functions. Endothelial cells play critical roles in regulating vascular homeostasis, such as roles in forming a selective barrier, inflammation, hemostasis, and vascular tone, and endothelial dysfunction is a hallmark of atherosclerosis and hypertension. To test this hypothesis, we generate an integrated map of gene expression, open chromatin region, and 3D interactions in resting and TNFα-treated human endothelial cells. RESULTS: We show that genetic variants associated with CAD and BP are enriched in open chromatin regions identified in endothelial cells. We identify physical loops by Hi-C and link open chromatin peaks that include CAD or BP SNPs with the promoters of genes expressed in endothelial cells. This analysis highlights 991 combinations of open chromatin regions and gene promoters that map to 38 CAD and 92 BP GWAS loci. We validate one CAD locus, by engineering a deletion of the TNFα-sensitive regulatory element using CRISPR/Cas9 and measure the effect on the expression of the novel CAD candidate gene AIDA. CONCLUSIONS: Our data support an important role played by genetic variants acting in the vascular endothelium to modulate inter-individual risk in CAD and hypertension.

Correction to: Functional variation in allelic methylomes underscores a strong genetic contribution and reveals novel epigenetic alterations in the human epigenome.

Following publication of the original article [1], the authors reported an error in Additional file 1.