Epigenetic variation impacts individual differences in the transcriptional response to influenza infection.

Humans display remarkable interindividual variation in their immune response to identical challenges. Yet, our understanding of the genetic and epigenetic factors contributing to such variation remains limited. Here we performed in-depth genetic, epigenetic and transcriptional profiling on primary macrophages derived from individuals of European and African ancestry before and after infection with influenza A virus. We show that baseline epigenetic profiles are strongly predictive of the transcriptional response to influenza A virus across individuals. Quantitative trait locus (QTL) mapping revealed highly coordinated genetic effects on gene regulation, with many cis-acting genetic variants impacting concomitantly gene expression and multiple epigenetic marks. These data reveal that ancestry-associated differences in the epigenetic landscape can be genetically controlled, even more than gene expression. Lastly, among QTL variants that colocalized with immune-disease loci, only 7% were gene expression QTL, while the remaining genetic variants impact epigenetic marks, stressing the importance of considering molecular phenotypes beyond gene expression in disease-focused studies.

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.

Genome graphs detect human polymorphisms in active epigenomic state during influenza infection.

Genetic variants, including mobile element insertions (MEIs), are known to impact the epigenome. We hypothesized that genome graphs, which encapsulate genetic diversity, could reveal missing epigenomic signals. To test this, we sequenced the epigenome of monocyte-derived macrophages from 35 ancestrally diverse individuals before and after influenza infection, allowing us to investigate the role of MEIs in immunity. We characterized genetic variants and MEIs using linked reads and built a genome graph. Mapping epigenetic data revealed 2.3%-3% novel peaks for H3K4me1, H3K27ac chromatin immunoprecipitation sequencing (ChIP-seq), and ATAC-seq. Additionally, the use of a genome graph modified some quantitative trait loci estimates and revealed 375 polymorphic MEIs in an active epigenomic state. Among these is an AluYh3 polymorphism whose chromatin state changed after infection and was associated with the expression of TRIM25, a gene that restricts influenza RNA synthesis. Our results demonstrate that graph genomes can reveal regulatory regions that would have been overlooked by other approaches.

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.

Large-scale analysis of circulating glutamate and adipose gene expression in relation to abdominal obesity.

Circulating levels of the amino acid glutamate are associated with central fat accumulation, yet the pathophysiology of this relationship remains unknown. We aimed to (i) refine and validate the association between circulating glutamate and abdominal obesity in a large twin cohort, and (ii) investigate whether transcriptomic profiles in adipose tissue could provide insight into the biological mechanisms underlying the association. First, in a cohort of 4665 individuals from the TwinsUK resource, we identified individuals with abdominal obesity and compared prevalence of the latter across circulating glutamate quintiles. Second, we used transcriptomic signatures generated from adipose tissue, both subcutaneous and visceral, to investigate associations with circulating glutamate levels. Individuals in the top circulating glutamate quintile had a sevenfold higher prevalence of abdominal obesity compared to those in the bottom quintile. The adipose tissue transcriptomic analyses identified GLUL, encoding Glutamate-Ammonia Ligase, as being associated with circulating glutamate and abdominal obesity, with pronounced signatures in the visceral depot. In conclusion, circulating glutamate is positively associated with the prevalence of abdominal obesity which relates to dysregulated GLUL expression specifically in visceral adipose tissue.

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.

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.

Single-cell analysis of human adipose tissue identifies depot and disease specific cell types.

The complex relationship between metabolic disease risk and body fat distribution in humans involves cellular characteristics which are specific to body fat compartments. Here we show depot-specific differences in the stromal vascual fraction of visceral and subcutaneous adipose tissue by performing single-cell RNA sequencing of tissue specimen from obese individuals. We characterize multiple immune cells, endothelial cells, fibroblasts, adipose and hematopoietic stem cell progenitors. Subpopulations of adipose-resident immune cells are metabolically active and associated with metabolic disease status and those include a population of potential dysfunctional CD8+ T cells expressing metallothioneins. We identify multiple types of adipocyte progenitors that are common across depots, including a subtype enriched in individuals with type 2 diabetes. Depot-specific analysis reveals a class of adipocyte progenitors unique to visceral adipose tissue, which shares common features with beige preadipocytes. Our human single-cell transcriptome atlas across fat depots provides a resource to dissect functional genomics of metabolic disease.

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.

Author Correction: An atlas of genetic influences on osteoporosis in humans and mice.

In the version of this article initially published, in Fig. 5a, the data in the right column of 'DAAM2 gRNA1' were incorrectly plotted as circles indicating 'untreated' rather than as squares indicating 'treated'. The error has been corrected in the HTML and PDF versions of the article.

Dissecting features of epigenetic variants underlying cardiometabolic risk using full-resolution epigenome profiling in regulatory elements.

Sparse profiling of CpG methylation in blood by microarrays has identified epigenetic links to common diseases. Here we apply methylC-capture sequencing (MCC-Seq) in a clinical population of ~200 adipose tissue and matched blood samples (Ntotal~400), providing high-resolution methylation profiling (>1.3 M CpGs) at regulatory elements. We link methylation to cardiometabolic risk through associations to circulating plasma lipid levels and identify lipid-associated CpGs with unique localization patterns in regulatory elements. We show distinct features of tissue-specific versus tissue-independent lipid-linked regulatory regions by contrasting with parallel assessments in ~800 independent adipose tissue and blood samples from the general population. We follow-up on adipose-specific regulatory regions under (1) genetic and (2) epigenetic (environmental) regulation via integrational studies. Overall, the comprehensive sequencing of regulatory element methylomes reveals a rich landscape of functional variants linked genetically as well as epigenetically to plasma lipid traits.

An atlas of genetic influences on osteoporosis in humans and mice.

Osteoporosis is a common aging-related disease diagnosed primarily using bone mineral density (BMD). We assessed genetic determinants of BMD as estimated by heel quantitative ultrasound in 426,824 individuals, identifying 518 genome-wide significant loci (301 novel), explaining 20% of its variance. We identified 13 bone fracture loci, all associated with estimated BMD (eBMD), in ~1.2 million individuals. We then identified target genes enriched for genes known to influence bone density and strength (maximum odds ratio (OR) = 58, P = 1 × 10-75) from cell-specific features, including chromatin conformation and accessible chromatin sites. We next performed rapid-throughput skeletal phenotyping of 126 knockout mice with disruptions in predicted target genes and found an increased abnormal skeletal phenotype frequency compared to 526 unselected lines (P < 0.0001). In-depth analysis of one gene, DAAM2, showed a disproportionate decrease in bone strength relative to mineralization. This genetic atlas provides evidence linking associated SNPs to causal genes, offers new insight into osteoporosis pathophysiology, and highlights opportunities for drug development.

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

BACKGROUND: The functional impact of genetic variation has been extensively surveyed, revealing that genetic changes correlated to phenotypes lie mostly in non-coding genomic regions. Studies have linked allele-specific genetic changes to gene expression, DNA methylation, and histone marks but these investigations have only been carried out in a limited set of samples. RESULTS: We describe a large-scale coordinated study of allelic and non-allelic effects on DNA methylation, histone mark deposition, and gene expression, detecting the interrelations between epigenetic and functional features at unprecedented resolution. We use information from whole genome and targeted bisulfite sequencing from 910 samples to perform genotype-dependent analyses of allele-specific methylation (ASM) and non-allelic methylation (mQTL). In addition, we introduce a novel genotype-independent test to detect methylation imbalance between chromosomes. Of the ~2.2 million CpGs tested for ASM, mQTL, and genotype-independent effects, we identify ~32% as being genetically regulated (ASM or mQTL) and ~14% as being putatively epigenetically regulated. We also show that epigenetically driven effects are strongly enriched in repressed regions and near transcription start sites, whereas the genetically regulated CpGs are enriched in enhancers. Known imprinted regions are enriched among epigenetically regulated loci, but we also observe several novel genomic regions (e.g., HOX genes) as being epigenetically regulated. Finally, we use our ASM datasets for functional interpretation of disease-associated loci and show the advantage of utilizing naïve T cells for understanding autoimmune diseases. CONCLUSIONS: Our rich catalogue of haploid methylomes across multiple tissues will allow validation of epigenome association studies and exploration of new biological models for allelic exclusion in the human genome.

Increased DNA methylation variability in type 1 diabetes across three immune effector cell types.

The incidence of type 1 diabetes (T1D) has substantially increased over the past decade, suggesting a role for non-genetic factors such as epigenetic mechanisms in disease development. Here we present an epigenome-wide association study across 406,365 CpGs in 52 monozygotic twin pairs discordant for T1D in three immune effector cell types. We observe a substantial enrichment of differentially variable CpG positions (DVPs) in T1D twins when compared with their healthy co-twins and when compared with healthy, unrelated individuals. These T1D-associated DVPs are found to be temporally stable and enriched at gene regulatory elements. Integration with cell type-specific gene regulatory circuits highlight pathways involved in immune cell metabolism and the cell cycle, including mTOR signalling. Evidence from cord blood of newborns who progress to overt T1D suggests that the DVPs likely emerge after birth. Our findings, based on 772 methylomes, implicate epigenetic changes that could contribute to disease pathogenesis in T1D.

Population whole-genome bisulfite sequencing across two tissues highlights the environment as the principal source of human methylome variation.

BACKGROUND: CpG methylation variation is involved in human trait formation and disease susceptibility. Analyses within populations have been biased towards CpG-dense regions through the application of targeted arrays. We generate whole-genome bisulfite sequencing data for approximately 30 adipose and blood samples from monozygotic and dizygotic twins for the characterization of non-genetic and genetic effects at single-site resolution. RESULTS: Purely invariable CpGs display a bimodal distribution with enrichment of unmethylated CpGs and depletion of fully methylated CpGs in promoter and enhancer regions. Population-variable CpGs account for approximately 15-20 % of total CpGs per tissue, are enriched in enhancer-associated regions and depleted in promoters, and single nucleotide polymorphisms at CpGs are a frequent confounder of extreme methylation variation. Differential methylation is primarily non-genetic in origin, with non-shared environment accounting for most of the variance. These non-genetic effects are mainly tissue-specific. Tobacco smoking is associated with differential methylation in blood with no evidence of this exposure impacting cell counts. Opposite to non-genetic effects, genetic effects of CpG methylation are shared across tissues and thus limit inter-tissue epigenetic drift. CpH methylation is rare, and shows similar characteristics of variation patterns as CpGs. CONCLUSIONS: Our study highlights the utility of low pass whole-genome bisulfite sequencing in identifying methylome variation beyond promoter regions, and suggests that targeting the population dynamic methylome of tissues requires assessment of understudied intergenic CpGs distal to gene promoters to reveal the full extent of inter-individual variation.

Characterization of functional methylomes by next-generation capture sequencing identifies novel disease-associated variants.

Most genome-wide methylation studies (EWAS) of multifactorial disease traits use targeted arrays or enrichment methodologies preferentially covering CpG-dense regions, to characterize sufficiently large samples. To overcome this limitation, we present here a new customizable, cost-effective approach, methylC-capture sequencing (MCC-Seq), for sequencing functional methylomes, while simultaneously providing genetic variation information. To illustrate MCC-Seq, we use whole-genome bisulfite sequencing on adipose tissue (AT) samples and public databases to design AT-specific panels. We establish its efficiency for high-density interrogation of methylome variability by systematic comparisons with other approaches and demonstrate its applicability by identifying novel methylation variation within enhancers strongly correlated to plasma triglyceride and HDL-cholesterol, including at CD36. Our more comprehensive AT panel assesses tissue methylation and genotypes in parallel at ∼4 and ∼3 M sites, respectively. Our study demonstrates that MCC-Seq provides comparable accuracy to alternative approaches but enables more efficient cataloguing of functional and disease-relevant epigenetic and genetic variants for large-scale EWAS.