Genetic effects on the skin methylome in healthy older twins.

Whole-skin DNA methylation variation has been implicated in several diseases, including melanoma, but its genetic basis has not yet been fully characterized. Using bulk skin tissue samples from 414 healthy female UK twins, we performed twin-based heritability and methylation quantitative trait loci (meQTL) analyses for >400,000 DNA methylation sites. We find that the human skin DNA methylome is on average less heritable than previously estimated in blood and other tissues (mean heritability: 10.02%). meQTL analysis identified local genetic effects influencing DNA methylation at 18.8% (76,442) of tested CpG sites, as well as 1,775 CpG sites associated with at least one distal genetic variant. As a functional follow-up, we performed skin expression QTL (eQTL) analyses in a partially overlapping sample of 604 female twins. Colocalization analysis identified over 3,500 shared genetic effects affecting thousands of CpG sites (10,067) and genes (4,475). Mediation analysis of putative colocalized gene-CpG pairs identified 114 genes with evidence for eQTL effects being mediated by DNA methylation in skin, including in genes implicating skin disease such as ALOX12 and CSPG4. We further explored the relevance of skin meQTLs to skin disease and found that skin meQTLs and CpGs under genetic influence were enriched for multiple skin-related genome-wide and epigenome-wide association signals, including for melanoma and psoriasis. Our findings give insights into the regulatory landscape of epigenomic variation in skin.

Thyroid function, pernicious anemia and erythropoiesis: a two-sample Mendelian randomization study.

Autoimmune thyroid disease (AITD) and pernicious anemia (PA) often coexist, but the directionality is unknown. In a two-sample Mendelian randomization (MR) analysis, using summary statistics from large genome-wide association studies (GWASs) in Europeans (N = 49 269-755 406), we examined the genetic associations between thyroid function, PA and markers of erythropoiesis. We performed inverse variance weighted random-effects MR, several sensitivity MR analyses, and bidirectional MR and MR Steiger for directionality. AITD and PA were associated bidirectionally (P ≤ 8 × 10-6). Neither euthyroid thyroid stimulating hormone (TSH) nor free thyroxine (FT4) were causally associated with PA. One standard deviation (SD) increase in euthyroid FT4 regulated by genetic variants in deiodinases 1 and 2 genes (DIO1/DIO2), corresponding to low-normal free triiodothyronine (FT3) levels, was causally associated with a pernicious/macrocytic anemia pattern, i.e. decreased erythrocyte counts (rank-based inverse normal transformed ß = -0,064 [95% confidence interval: -0,085, -0,044], P = 8 × 10-10) and hemoglobin (-0.028 [-0.051, -0.005], P = 0.02) and increased mean corpuscular hemoglobin (0.058 [0.025, 0.091], P = 5 × 10-4) and mean corpuscular volume levels (0.075 [0.052, 0.098], P = 1 × 10-8). Meanwhile, subclinical hyperthyroidism mirrored that pattern. AITD was causally associated with increased erythrocyte distribution width (P = 0.007) and decreased reticulocyte counts (P ≤ 0.02), whereas high-normal FT4 regulated by DIO1/DIO2 variants was causally associated with decreased bilirubin (-0.039 (-0.064, -0.013), P = 0.003). In conclusion, the bidirectional association between AITD and PA suggests a shared heritability for these two autoimmune diseases. AITD was causally associated with impaired erythropoiesis and not autoimmune hemolysis. Additionally, in euthyroid individuals, local regulation of thyroid hormones by deiodinases likely plays a role in erythropoiesis.

Association of Factor V Leiden With Subsequent Atherothrombotic Events: A GENIUS-CHD Study of Individual Participant Data.

BACKGROUND: Studies examining the role of factor V Leiden among patients at higher risk of atherothrombotic events, such as those with established coronary heart disease (CHD), are lacking. Given that coagulation is involved in the thrombus formation stage on atherosclerotic plaque rupture, we hypothesized that factor V Leiden may be a stronger risk factor for atherothrombotic events in patients with established CHD. METHODS: We performed an individual-level meta-analysis including 25 prospective studies (18 cohorts, 3 case-cohorts, 4 randomized trials) from the GENIUS-CHD (Genetics of Subsequent Coronary Heart Disease) consortium involving patients with established CHD at baseline. Participating studies genotyped factor V Leiden status and shared risk estimates for the outcomes of interest using a centrally developed statistical code with harmonized definitions across studies. Cox proportional hazards regression models were used to obtain age- and sex-adjusted estimates. The obtained estimates were pooled using fixed-effect meta-analysis. The primary outcome was composite of myocardial infarction and CHD death. Secondary outcomes included any stroke, ischemic stroke, coronary revascularization, cardiovascular mortality, and all-cause mortality. RESULTS: The studies included 69 681 individuals of whom 3190 (4.6%) were either heterozygous or homozygous (n=47) carriers of factor V Leiden. Median follow-up per study ranged from 1.0 to 10.6 years. A total of 20 studies with 61 147 participants and 6849 events contributed to analyses of the primary outcome. Factor V Leiden was not associated with the combined outcome of myocardial infarction and CHD death (hazard ratio, 1.03 [95% CI, 0.92-1.16]; I2=28%; P-heterogeneity=0.12). Subgroup analysis according to baseline characteristics or strata of traditional cardiovascular risk factors did not show relevant differences. Similarly, risk estimates for the secondary outcomes including stroke, coronary revascularization, cardiovascular mortality, and all-cause mortality were also close to identity. CONCLUSIONS: Factor V Leiden was not associated with increased risk of subsequent atherothrombotic events and mortality in high-risk participants with established and treated CHD. Routine assessment of factor V Leiden status is unlikely to improve atherothrombotic events risk stratification in this population.

Thyroid function, sex hormones and sexual function: a Mendelian randomization study.

Hypothyroidism and hyperthyroidism are observationally associated with sex hormone concentrations and sexual dysfunction, but causality is unclear. We investigated whether TSH, fT4, hypo- and hyperthyroidism are causally associated with sex hormones and sexual function. We used publicly available summary statistics from genome-wide association studies on TSH and fT4 and hypo- and hyperthyroidism from the ThyroidOmics Consortium (N ≤ 54,288). Outcomes from UK Biobank (women ≤ 194,174/men ≤ 167,020) and ReproGen (women ≤ 252,514) were sex hormones (sex hormone binding globulin [SHBG], testosterone, estradiol, free androgen index [FAI]) and sexual function (ovulatory function in women: duration of menstrual period, age at menarche and menopause, reproductive lifespan, and erectile dysfunction in men). We performed two-sample Mendelian randomization (MR) analyses on summary level, and unweighted genetic risk score (GRS) analysis on individual level data. One SD increase in TSH was associated with a 1.332 nmol/L lower (95% CI: - 0.717,- 1.946; p = 2 × 10-5) SHBG and a 0.103 nmol/l lower (- 0.051,V0.154; p = 9 × 10-5) testosterone in two-sample MR, supported by the GRS approach. Genetic predisposition to hypothyroidism was associated with decreased and genetic predisposition to hyperthyroidism with increased SHBG and testosterone in both approaches. The GRS for fT4 was associated with increased testosterone and estradiol in women only. The GRS for TSH and hypothyroidism were associated with increased and the GRS for hyperthyroidism with decreased FAI in men only. While genetically predicted thyroid function was associated with sex hormones, we found no association with sexual function.

Genetically modulated educational attainment and coronary disease risk.

AIMS: Genetic disposition and lifestyle factors are understood as independent components underlying the risk of multiple diseases. In this study, we aim to investigate the interplay between genetics, educational attainment-an important denominator of lifestyle-and coronary artery disease (CAD) risk. METHODS AND RESULTS: Based on the effect sizes of 74 genetic variants associated with educational attainment, we calculated a 'genetic education score' in 13 080 cases and 14 471 controls and observed an inverse correlation between the score and risk of CAD [P = 1.52 × 10-8; odds ratio (OR) 0.79, 95% confidence interval (CI) 0.73-0.85 for the higher compared with the lowest score quintile]. We replicated in 146 514 individuals from UK Biobank (P = 1.85 × 10-6) and also found strong associations between the 'genetic education score' with 'modifiable' risk factors including smoking (P = 5.36 × 10-23), body mass index (BMI) (P = 1.66 × 10-30), and hypertension (P = 3.86 × 10-8). Interestingly, these associations were only modestly attenuated by adjustment for years spent in school. In contrast, a model adjusting for BMI and smoking abolished the association signal between the 'genetic education score' and CAD risk suggesting an intermediary role of these two risk factors. Mendelian randomization analyses performed with summary statistics from large genome-wide meta-analyses and sensitivity analysis using 1271 variants affecting educational attainment (OR 0.68 for the higher compared with the lowest score quintile; 95% CI 0.63-0.74; P = 3.99 × 10-21) further strengthened these findings. CONCLUSION: Genetic variants known to affect educational attainment may have implications for a health-conscious lifestyle later in life and subsequently affect the risk of CAD.

Identification of Genetic Variants Associated With Myocardial Infarction in Saudi Arabia.

The genetic variants associated with various genetic disorders have not been identified decisively in Saudi Arabia. Among these variants, six known for their association with coronary artery disease or myocardial infarction (MI) were studied on Saudi patients. Reference single nucleotide polymorphisms (SNPs) of these variants are rs5174, rs11591147, rs2259816, rs111245230, rs3782886 and rs2259820, referring to genes LRP8, PCSK9, HNF1A, SVEP1, BRAP and HNF1A, respectively. The analysis employed polymerase chain reaction panel coupled with mini-sequencing (SNapShot multiplex system) in order to identify these variants. A total of 100 MI patients and 103 healthy control individuals participated in this study. The six variants (SNPs) were evaluated for the risk of developing MI in the Saudi patients. Analysis of allele frequencies indicated that A allele of rs11591147 variant can be a protective allele, thus, is associated with the decreased risk of MI in Saudi individuals. Rare allele of rs111245230 variant (e.g., C allele) was extremely reduced, while rare allele of rs3782886 variant (e.g., G allele) does not exist in the ethnic signature of the Saudi population. This study elucidates the possible prediction of risk factors associated with severe diseases in Saudi population utilizing SNapShot multiplex system.

Polymorphism in a lincRNA Associates with a Doubled Risk of Pneumococcal Bacteremia in Kenyan Children.

Bacteremia (bacterial bloodstream infection) is a major cause of illness and death in sub-Saharan Africa but little is known about the role of human genetics in susceptibility. We conducted a genome-wide association study of bacteremia susceptibility in more than 5,000 Kenyan children as part of the Wellcome Trust Case Control Consortium 2 (WTCCC2). Both the blood-culture-proven bacteremia case subjects and healthy infants as controls were recruited from Kilifi, on the east coast of Kenya. Streptococcus pneumoniae is the most common cause of bacteremia in Kilifi and was thus the focus of this study. We identified an association between polymorphisms in a long intergenic non-coding RNA (lincRNA) gene (AC011288.2) and pneumococcal bacteremia and replicated the results in the same population (p combined = 1.69 × 10(-9); OR = 2.47, 95% CI = 1.84-3.31). The susceptibility allele is African specific, derived rather than ancestral, and occurs at low frequency (2.7% in control subjects and 6.4% in case subjects). Our further studies showed AC011288.2 expression only in neutrophils, a cell type that is known to play a major role in pneumococcal clearance. Identification of this novel association will further focus research on the role of lincRNAs in human infectious disease.

Platelet-Related Variants Identified by Exomechip Meta-analysis in 157,293 Individuals.

Platelet production, maintenance, and clearance are tightly controlled processes indicative of platelets' important roles in hemostasis and thrombosis. Platelets are common targets for primary and secondary prevention of several conditions. They are monitored clinically by complete blood counts, specifically with measurements of platelet count (PLT) and mean platelet volume (MPV). Identifying genetic effects on PLT and MPV can provide mechanistic insights into platelet biology and their role in disease. Therefore, we formed the Blood Cell Consortium (BCX) to perform a large-scale meta-analysis of Exomechip association results for PLT and MPV in 157,293 and 57,617 individuals, respectively. Using the low-frequency/rare coding variant-enriched Exomechip genotyping array, we sought to identify genetic variants associated with PLT and MPV. In addition to confirming 47 known PLT and 20 known MPV associations, we identified 32 PLT and 18 MPV associations not previously observed in the literature across the allele frequency spectrum, including rare large effect (FCER1A), low-frequency (IQGAP2, MAP1A, LY75), and common (ZMIZ2, SMG6, PEAR1, ARFGAP3/PACSIN2) variants. Several variants associated with PLT/MPV (PEAR1, MRVI1, PTGES3) were also associated with platelet reactivity. In concurrent BCX analyses, there was overlap of platelet-associated variants with red (MAP1A, TMPRSS6, ZMIZ2) and white (PEAR1, ZMIZ2, LY75) blood cell traits, suggesting common regulatory pathways with shared genetic architecture among these hematopoietic lineages. Our large-scale Exomechip analyses identified previously undocumented associations with platelet traits and further indicate that several complex quantitative hematological, lipid, and cardiovascular traits share genetic factors.

Emerging applications of genome-editing technology to examine functionality of GWAS-associated variants for complex traits.

Over the last decade, genome-wide association studies (GWAS) have propelled the discovery of thousands of loci associated with complex diseases. The focus is now turning toward the function of these association signals, determining the causal variant(s) among those in strong linkage disequilibrium, and identifying their underlying mechanisms, such as long-range gene regulation. Genome-editing techniques utilizing zinc-finger nucleases (ZFN), transcription activator-like effector nucleases (TALENs), and clustered regularly-interspaced short palindromic repeats with Cas9 nuclease (CRISPR-Cas9) are becoming the tools of choice to establish functionality for these variants, due to the ability to assess effects of single variants in vivo. This review will discuss examples of how these technologies have begun to aid functional analysis of GWAS loci for complex traits such as cardiovascular disease, Type 2 diabetes, cancer, obesity, and autoimmune disease. We focus on analysis of variants occurring within noncoding genomic regions, as these comprise the majority of GWAS variants, providing the greatest challenges to determining functionality, and compare editing strategies that provide different levels of evidence for variant functionality. The review describes molecular insights into some of these potentially causal variants and how these may relate to the pathology of the trait and look toward future directions for these technologies in post-GWAS analysis, such as base-editing.

Meta-analysis of 65,734 individuals identifies TSPAN15 and SLC44A2 as two susceptibility loci for venous thromboembolism.

Venous thromboembolism (VTE), the third leading cause of cardiovascular mortality, is a complex thrombotic disorder with environmental and genetic determinants. Although several genetic variants have been found associated with VTE, they explain a minor proportion of VTE risk in cases. We undertook a meta-analysis of genome-wide association studies (GWASs) to identify additional VTE susceptibility genes. Twelve GWASs totaling 7,507 VTE case subjects and 52,632 control subjects formed our discovery stage where 6,751,884 SNPs were tested for association with VTE. Nine loci reached the genome-wide significance level of 5 × 10(-8) including six already known to associate with VTE (ABO, F2, F5, F11, FGG, and PROCR) and three unsuspected loci. SNPs mapping to these latter were selected for replication in three independent case-control studies totaling 3,009 VTE-affected individuals and 2,586 control subjects. This strategy led to the identification and replication of two VTE-associated loci, TSPAN15 and SLC44A2, with lead risk alleles associated with odds ratio for disease of 1.31 (p = 1.67 × 10(-16)) and 1.21 (p = 2.75 × 10(-15)), respectively. The lead SNP at the TSPAN15 locus is the intronic rs78707713 and the lead SLC44A2 SNP is the non-synonymous rs2288904 previously shown to associate with transfusion-related acute lung injury. We further showed that these two variants did not associate with known hemostatic plasma markers. TSPAN15 and SLC44A2 do not belong to conventional pathways for thrombosis and have not been associated to other cardiovascular diseases nor related quantitative biomarkers. Our findings uncovered unexpected actors of VTE etiology and pave the way for novel mechanistic concepts of VTE pathophysiology.

Epigenome-wide association study (EWAS) of BMI, BMI change and waist circumference in African American adults identifies multiple replicated loci.

Obesity is an important component of the pathophysiology of chronic diseases. Identifying epigenetic modifications associated with elevated adiposity, including DNA methylation variation, may point to genomic pathways that are dysregulated in numerous conditions. The Illumina 450K Bead Chip array was used to assay DNA methylation in leukocyte DNA obtained from 2097 African American adults in the Atherosclerosis Risk in Communities (ARIC) study. Mixed-effects regression models were used to test the association of methylation beta value with concurrent body mass index (BMI) and waist circumference (WC), and BMI change, adjusting for batch effects and potential confounders. Replication using whole-blood DNA from 2377 White adults in the Framingham Heart Study and CD4+ T cell DNA from 991 Whites in the Genetics of Lipid Lowering Drugs and Diet Network Study was followed by testing using adipose tissue DNA from 648 women in the Multiple Tissue Human Expression Resource cohort. Seventy-six BMI-related probes, 164 WC-related probes and 8 BMI change-related probes passed the threshold for significance in ARIC (P < 1 × 10(-7); Bonferroni), including probes in the recently reported HIF3A, CPT1A and ABCG1 regions. Replication using blood DNA was achieved for 37 BMI probes and 1 additional WC probe. Sixteen of these also replicated in adipose tissue, including 15 novel methylation findings near genes involved in lipid metabolism, immune response/cytokine signaling and other diverse pathways, including LGALS3BP, KDM2B, PBX1 and BBS2, among others. Adiposity traits are associated with DNA methylation at numerous CpG sites that replicate across studies despite variation in tissue type, ethnicity and analytic approaches.

Gene × dietary pattern interactions in obesity: analysis of up to 68 317 adults of European ancestry.

Obesity is highly heritable. Genetic variants showing robust associations with obesity traits have been identified through genome-wide association studies. We investigated whether a composite score representing healthy diet modifies associations of these variants with obesity traits. Totally, 32 body mass index (BMI)- and 14 waist-hip ratio (WHR)-associated single nucleotide polymorphisms were genotyped, and genetic risk scores (GRS) were calculated in 18 cohorts of European ancestry (n = 68 317). Diet score was calculated based on self-reported intakes of whole grains, fish, fruits, vegetables, nuts/seeds (favorable) and red/processed meats, sweets, sugar-sweetened beverages and fried potatoes (unfavorable). Multivariable adjusted, linear regression within each cohort followed by inverse variance-weighted, fixed-effects meta-analysis was used to characterize: (a) associations of each GRS with BMI and BMI-adjusted WHR and (b) diet score modification of genetic associations with BMI and BMI-adjusted WHR. Nominally significant interactions (P = 0.006-0.04) were observed between the diet score and WHR-GRS (but not BMI-GRS), two WHR loci (GRB14 rs10195252; LYPLAL1 rs4846567) and two BMI loci (LRRN6C rs10968576; MTIF3 rs4771122), for the respective BMI-adjusted WHR or BMI outcomes. Although the magnitudes of these select interactions were small, our data indicated that associations between genetic predisposition and obesity traits were stronger with a healthier diet. Our findings generate interesting hypotheses; however, experimental and functional studies are needed to determine their clinical relevance.

Genome-wide association study of nevirapine hypersensitivity in a sub-Saharan African HIV-infected population.

Background: The antiretroviral nevirapine is associated with hypersensitivity reactions in 6%-10% of patients, including hepatotoxicity, maculopapular exanthema, Stevens-Johnson syndrome (SJS) and toxic epidermal necrolysis (TEN). Objectives: To undertake a genome-wide association study (GWAS) to identify genetic predisposing factors for the different clinical phenotypes associated with nevirapine hypersensitivity. Methods: A GWAS was undertaken in a discovery cohort of 151 nevirapine-hypersensitive and 182 tolerant, HIV-infected Malawian adults. Replication of signals was determined in a cohort of 116 cases and 68 controls obtained from Malawi, Uganda and Mozambique. Interaction with ERAP genes was determined in patients positive for HLA-C*04:01 . In silico docking studies were also performed for HLA-C*04:01 . Results: Fifteen SNPs demonstrated nominal significance ( P < 1 × 10 -5 ) with one or more of the hypersensitivity phenotypes. The most promising signal was seen in SJS/TEN, where rs5010528 ( HLA-C locus) approached genome-wide significance ( P < 8.5 × 10 -8 ) and was below HLA -wide significance ( P < 2.5 × 10 -4 ) in the meta-analysis of discovery and replication cohorts [OR 4.84 (95% CI 2.71-8.61)]. rs5010528 is a strong proxy for HLA-C*04:01 carriage: in silico docking showed that two residues (33 and 123) in the B pocket were the most likely nevirapine interactors. There was no interaction between HLA-C*04:01 and ERAP1 , but there is a potential protective effect with ERAP2 [ P = 0.019, OR 0.43 (95% CI 0.21-0.87)]. Conclusions: HLA-C*04:01 predisposes to nevirapine-induced SJS/TEN in sub-Saharan Africans, but not to other hypersensitivity phenotypes. This is likely to be mediated via binding to the B pocket of the HLA-C peptide. Whether this risk is modulated by ERAP2 variants requires further study.

Human metabolic individuality in biomedical and pharmaceutical research.

Genome-wide association studies (GWAS) have identified many risk loci for complex diseases, but effect sizes are typically small and information on the underlying biological processes is often lacking. Associations with metabolic traits as functional intermediates can overcome these problems and potentially inform individualized therapy. Here we report a comprehensive analysis of genotype-dependent metabolic phenotypes using a GWAS with non-targeted metabolomics. We identified 37 genetic loci associated with blood metabolite concentrations, of which 25 show effect sizes that are unusually high for GWAS and account for 10-60% differences in metabolite levels per allele copy. Our associations provide new functional insights for many disease-related associations that have been reported in previous studies, including those for cardiovascular and kidney disorders, type 2 diabetes, cancer, gout, venous thromboembolism and Crohn's disease. The study advances our knowledge of the genetic basis of metabolic individuality in humans and generates many new hypotheses for biomedical and pharmaceutical research.

Cis and trans effects of human genomic variants on gene expression.

Gene expression is a heritable cellular phenotype that defines the function of a cell and can lead to diseases in case of misregulation. In order to detect genetic variations affecting gene expression, we performed association analysis of single nucleotide polymorphisms (SNPs) and copy number variants (CNVs) with gene expression measured in 869 lymphoblastoid cell lines of the Avon Longitudinal Study of Parents and Children (ALSPAC) cohort in cis and in trans. We discovered that 3,534 genes (false discovery rate (FDR) = 5%) are affected by an expression quantitative trait locus (eQTL) in cis and 48 genes are affected in trans. We observed that CNVs are more likely to be eQTLs than SNPs. In addition, we found that variants associated to complex traits and diseases are enriched for trans-eQTLs and that trans-eQTLs are enriched for cis-eQTLs. As a variant affecting both a gene in cis and in trans suggests that the cis gene is functionally linked to the trans gene expression, we looked specifically for trans effects of cis-eQTLs. We discovered that 26 cis-eQTLs are associated to 92 genes in trans with the cis-eQTLs of the transcriptions factors BATF3 and HMX2 affecting the most genes. We then explored if the variation of the level of expression of the cis genes were causally affecting the level of expression of the trans genes and discovered several causal relationships between variation in the level of expression of the cis gene and variation of the level of expression of the trans gene. This analysis shows that a large sample size allows the discovery of secondary effects of human variations on gene expression that can be used to construct short directed gene regulatory networks.

Obesity accelerates epigenetic aging of human liver.

Because of the dearth of biomarkers of aging, it has been difficult to test the hypothesis that obesity increases tissue age. Here we use a novel epigenetic biomarker of aging (referred to as an "epigenetic clock") to study the relationship between high body mass index (BMI) and the DNA methylation ages of human blood, liver, muscle, and adipose tissue. A significant correlation between BMI and epigenetic age acceleration could only be observed for liver (r = 0.42, P = 6.8 × 10(-4) in dataset 1 and r = 0.42, P = 1.2 × 10(-4) in dataset 2). On average, epigenetic age increased by 3.3 y for each 10 BMI units. The detected age acceleration in liver is not associated with the Nonalcoholic Fatty Liver Disease Activity Score or any of its component traits after adjustment for BMI. The 279 genes that are underexpressed in older liver samples are highly enriched (1.2 × 10(-9)) with nuclear mitochondrial genes that play a role in oxidative phosphorylation and electron transport. The epigenetic age acceleration, which is not reversible in the short term after rapid weight loss induced by bariatric surgery, may play a role in liver-related comorbidities of obesity, such as insulin resistance and liver cancer.

Higher chylomicron remnants and LDL particle numbers associate with CD36 SNPs and DNA methylation sites that reduce CD36.

Cluster of differentiation 36 (CD36) variants influence fasting lipids and risk of metabolic syndrome, but their impact on postprandial lipids, an independent risk factor for cardiovascular disease, is unclear. We determined the effects of SNPs within a ∼410 kb region encompassing CD36 and its proximal and distal promoters on chylomicron (CM) remnants and LDL particles at fasting and at 3.5 and 6 h following a high-fat meal (Genetics of Lipid Lowering Drugs and Diet Network study, n = 1,117). Five promoter variants associated with CMs, four with delayed TG clearance and five with LDL particle number. To assess mechanisms underlying the associations, we queried expression quantitative trait loci, DNA methylation, and ChIP-seq datasets for adipose and heart tissues that function in postprandial lipid clearance. Several SNPs that associated with higher serum lipids correlated with lower adipose and heart CD36 mRNA and aligned to active motifs for PPARγ, a major CD36 regulator. The SNPs also associated with DNA methylation sites that related to reduced CD36 mRNA and higher serum lipids, but mixed-model analyses indicated that the SNPs and methylation independently influence CD36 mRNA. The findings support contributions of CD36 SNPs that reduce adipose and heart CD36 RNA expression to inter-individual variability of postprandial lipid metabolism and document changes in CD36 DNA methylation that influence both CD36 expression and lipids.

Global analysis of DNA methylation variation in adipose tissue from twins reveals links to disease-associated variants in distal regulatory elements.

Epigenetic modifications such as DNA methylation play a key role in gene regulation and disease susceptibility. However, little is known about the genome-wide frequency, localization, and function of methylation variation and how it is regulated by genetic and environmental factors. We utilized the Multiple Tissue Human Expression Resource (MuTHER) and generated Illumina 450K adipose methylome data from 648 twins. We found that individual CpGs had low variance and that variability was suppressed in promoters. We noted that DNA methylation variation was highly heritable (h(2)median = 0.34) and that shared environmental effects correlated with metabolic phenotype-associated CpGs. Analysis of methylation quantitative-trait loci (metQTL) revealed that 28% of CpGs were associated with nearby SNPs, and when overlapping them with adipose expression quantitative-trait loci (eQTL) from the same individuals, we found that 6% of the loci played a role in regulating both gene expression and DNA methylation. These associations were bidirectional, but there were pronounced negative associations for promoter CpGs. Integration of metQTL with adipose reference epigenomes and disease associations revealed significant enrichment of metQTL overlapping metabolic-trait or disease loci in enhancers (the strongest effects were for high-density lipoprotein cholesterol and body mass index [BMI]). We followed up with the BMI SNP rs713586, a cg01884057 metQTL that overlaps an enhancer upstream of ADCY3, and used bisulphite sequencing to refine this region. Our results showed widespread population invariability yet sequence dependence on adipose DNA methylation but that incorporating maps of regulatory elements aid in linking CpG variation to gene regulation and disease risk in a tissue-dependent manner.

Exome-wide analysis of rare coding variation identifies novel associations with COPD and airflow limitation in MOCS3, IFIT3 and SERPINA12.

BACKGROUND: Several regions of the genome have shown to be associated with COPD in genome-wide association studies of common variants. OBJECTIVE: To determine rare and potentially functional single nucleotide polymorphisms (SNPs) associated with the risk of COPD and severity of airflow limitation. METHODS: 3226 current or former smokers of European ancestry with lung function measures indicative of Global Initiative for Chronic Obstructive Lung Disease (GOLD) 2 COPD or worse were genotyped using an exome array. An analysis of risk of COPD was carried out using ever smoking controls (n=4784). Associations with %predicted FEV1 were tested in cases. We followed-up signals of interest (p<10(-5)) in independent samples from a subset of the UK Biobank population and also undertook a more powerful discovery study by meta-analysing the exome array data and UK Biobank data for variants represented on both arrays. RESULTS: Among the associated variants were two in regions previously unreported for COPD; a low frequency non-synonymous SNP in MOCS3 (rs7269297, pdiscovery=3.08×10(-6), preplication=0.019) and a rare SNP in IFIT3, which emerged in the meta-analysis (rs140549288, pmeta=8.56×10(-6)). In the meta-analysis of % predicted FEV1 in cases, the strongest association was shown for a splice variant in a previously unreported region, SERPINA12 (rs140198372, pmeta=5.72×10(-6)). We also confirmed previously reported associations with COPD risk at MMP12, HHIP, GPR126 and CHRNA5. No associations in novel regions reached a stringent exome-wide significance threshold (p<3.7×10(-7)). CONCLUSIONS: This study identified several associations with the risk of COPD and severity of airflow limitation, including novel regions MOCS3, IFIT3 and SERPINA12, which warrant further study.

Maps of open chromatin highlight cell type-restricted patterns of regulatory sequence variation at hematological trait loci.

Nearly three-quarters of the 143 genetic signals associated with platelet and erythrocyte phenotypes identified by meta-analyses of genome-wide association (GWA) studies are located at non-protein-coding regions. Here, we assessed the role of candidate regulatory variants associated with cell type-restricted, closely related hematological quantitative traits in biologically relevant hematopoietic cell types. We used formaldehyde-assisted isolation of regulatory elements followed by next-generation sequencing (FAIRE-seq) to map regions of open chromatin in three primary human blood cells of the myeloid lineage. In the precursors of platelets and erythrocytes, as well as in monocytes, we found that open chromatin signatures reflect the corresponding hematopoietic lineages of the studied cell types and associate with the cell type-specific gene expression patterns. Dependent on their signal strength, open chromatin regions showed correlation with promoter and enhancer histone marks, distance to the transcription start site, and ontology classes of nearby genes. Cell type-restricted regions of open chromatin were enriched in sequence variants associated with hematological indices. The majority (63.6%) of such candidate functional variants at platelet quantitative trait loci (QTLs) coincided with binding sites of five transcription factors key in regulating megakaryopoiesis. We experimentally tested 13 candidate regulatory variants at 10 platelet QTLs and found that 10 (76.9%) affected protein binding, suggesting that this is a frequent mechanism by which regulatory variants influence quantitative trait levels. Our findings demonstrate that combining large-scale GWA data with open chromatin profiles of relevant cell types can be a powerful means of dissecting the genetic architecture of closely related quantitative traits.