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.

Escape from X-inactivation in twins exhibits intra- and inter-individual variability across tissues and is heritable.

X-chromosome inactivation (XCI) silences one X in female cells to balance sex-differences in X-dosage. A subset of X-linked genes escape XCI, but the extent to which this phenomenon occurs and how it varies across tissues and in a population is as yet unclear. To characterize incidence and variability of escape across individuals and tissues, we conducted a transcriptomic study of escape in adipose, skin, lymphoblastoid cell lines and immune cells in 248 healthy individuals exhibiting skewed XCI. We quantify XCI escape from a linear model of genes' allelic fold-change and XIST-based degree of XCI skewing. We identify 62 genes, including 19 lncRNAs, with previously unknown patterns of escape. We find a range of tissue-specificity, with 11% of genes escaping XCI constitutively across tissues and 23% demonstrating tissue-restricted escape, including cell type-specific escape across immune cells of the same individual. We also detect substantial inter-individual variability in escape. Monozygotic twins share more similar escape than dizygotic twins, indicating that genetic factors may underlie inter-individual differences in escape. However, discordant escape also occurs within monozygotic co-twins, suggesting environmental factors also influence escape. Altogether, these data indicate that XCI escape is an under-appreciated source of transcriptional differences, and an intricate phenotype impacting variable trait expressivity in females.

Cell-Type Heterogeneity in Adipose Tissue Is Associated with Complex Traits and Reveals Disease-Relevant Cell-Specific eQTLs.

Adipose tissue is an important endocrine organ with a role in many cardiometabolic diseases. It is comprised of a heterogeneous collection of cell types that can differentially impact disease phenotypes. Cellular heterogeneity can also confound -omic analyses but is rarely taken into account in analysis of solid-tissue transcriptomes. Here, we investigate cell-type heterogeneity in two population-level subcutaneous adipose-tissue RNA-seq datasets (TwinsUK, n = 766 and the Genotype-Tissue Expression project [GTEx], n = 326) by estimating the relative proportions of four distinct cell types (adipocytes, macrophages, CD4+ T cells, and micro-vascular endothelial cells). We find significant cellular heterogeneity within and between the TwinsUK and GTEx adipose datasets. We find that adipose cell-type composition is heritable and confirm the positive association between adipose-resident macrophage proportion and obesity (high BMI), but we find a stronger BMI-independent association with dual-energy X-ray absorptiometry (DXA) derived body-fat distribution traits. We benchmark the impact of adipose-tissue cell composition on a range of standard analyses, including phenotype-gene expression association, co-expression networks, and cis-eQTL discovery. Our results indicate that it is critical to account for cell-type composition when combining adipose transcriptome datasets in co-expression analysis and in differential expression analysis with obesity-related traits. We applied gene expression by cell-type proportion interaction models (G × Cell) to identify 26 cell-type-specific expression quantitative trait loci (eQTLs) in 20 genes, including four autoimmune disease genome-wide association study (GWAS) loci. These results identify cell-specific eQTLs and demonstrate the potential of in silico deconvolution of bulk tissue to identify cell-type-restricted regulatory variants.

ACE2 expression in adipose tissue is associated with cardio-metabolic risk factors and cell type composition-implications for COVID-19.

BACKGROUND: COVID-19 severity varies widely. Although some demographic and cardio-metabolic factors, including age and obesity, are associated with increasing risk of severe illness, the underlying mechanism(s) are uncertain. SUBJECTS/METHODS: In a meta-analysis of three independent studies of 1471 participants in total, we investigated phenotypic and genetic factors associated with subcutaneous adipose tissue expression of Angiotensin I Converting Enzyme 2 (ACE2), measured by RNA-Seq, which acts as a receptor for SARS-CoV-2 cellular entry. RESULTS: Lower adipose tissue ACE2 expression was associated with multiple adverse cardio-metabolic health indices, including type 2 diabetes (T2D) (P = 9.14 × 10-6), obesity status (P = 4.81 × 10-5), higher serum fasting insulin (P = 5.32 × 10-4), BMI (P = 3.94 × 10-4), and lower serum HDL levels (P = 1.92 × 10-7). ACE2 expression was also associated with estimated proportions of cell types in adipose tissue: lower expression was associated with a lower proportion of microvascular endothelial cells (P = 4.25 × 10-4) and higher proportion of macrophages (P = 2.74 × 10-5). Despite an estimated heritability of 32%, we did not identify any proximal or distal expression quantitative trait loci (eQTLs) associated with adipose tissue ACE2 expression. CONCLUSIONS: Our results demonstrate that individuals with cardio-metabolic features known to increase risk of severe COVID-19 have lower background ACE2 levels in this highly relevant tissue. Reduced adipose tissue ACE2 expression may contribute to the pathophysiology of cardio-metabolic diseases, as well as the associated increased risk of severe COVID-19.

Current smoking and COVID-19 risk: results from a population symptom app in over 2.4 million people.

BACKGROUND: The association between current tobacco smoking, the risk of developing symptomatic COVID-19 and the severity of illness is an important information gap. METHODS: UK users of the Zoe COVID-19 Symptom Study app provided baseline data including demographics, anthropometrics, smoking status and medical conditions, and were asked to log their condition daily. Participants who reported that they did not feel physically normal were then asked by the app to complete a series of questions, including 14 potential COVID-19 symptoms and about hospital attendance. The main study outcome was the development of 'classic' symptoms of COVID-19 during the pandemic defined as fever, new persistent cough and breathlessness and their association with current smoking. The number of concurrent COVID-19 symptoms was used as a proxy for severity and the pattern of association between symptoms was also compared between smokers and non-smokers. RESULTS: Between 24 March 2020 and 23 April 2020, data were available on 2 401 982 participants, mean (SD) age 43.6 (15.1) years, 63.3% female, overall smoking prevalence 11.0%. 834 437 (35%) participants reported being unwell and entered one or more symptoms. Current smokers were more likely to report symptoms suggesting a diagnosis of COVID-19; classic symptoms adjusted OR (95% CI) 1.14 (1.10 to 1.18); >5 symptoms 1.29 (1.26 to 1.31); >10 symptoms 1.50 (1.42 to 1.58). The pattern of association between reported symptoms did not vary between smokers and non-smokers. INTERPRETATION: These data are consistent with people who smoke being at an increased risk of developing symptomatic COVID-19.

Probable delirium is a presenting symptom of COVID-19 in frail, older adults: a cohort study of 322 hospitalised and 535 community-based older adults.

BACKGROUND: Frailty, increased vulnerability to physiological stressors, is associated with adverse outcomes. COVID-19 exhibits a more severe disease course in older, comorbid adults. Awareness of atypical presentations is critical to facilitate early identification. OBJECTIVE: To assess how frailty affects presenting COVID-19 symptoms in older adults. DESIGN: Observational cohort study of hospitalised older patients and self-report data for community-based older adults. SETTING: Admissions to St Thomas' Hospital, London with laboratory-confirmed COVID-19. Community-based data for older adults using the COVID Symptom Study mobile application. SUBJECTS: Hospital cohort: patients aged 65 and over (n = 322); unscheduled hospital admission between 1 March 2020 and 5 May 2020; COVID-19 confirmed by RT-PCR of nasopharyngeal swab. Community-based cohort: participants aged 65 and over enrolled in the COVID Symptom Study (n = 535); reported test-positive for COVID-19 from 24 March (application launch) to 8 May 2020. METHODS: Multivariable logistic regression analysis performed on age-matched samples from hospital and community-based cohorts to ascertain association of frailty with symptoms of confirmed COVID-19. RESULTS: Hospital cohort: significantly higher prevalence of probable delirium in the frail sample, with no difference in fever or cough. Community-based cohort: significantly higher prevalence of possible delirium in frailer, older adults and fatigue and shortness of breath. CONCLUSIONS: This is the first study demonstrating higher prevalence of probable delirium as a COVID-19 symptom in older adults with frailty compared to other older adults. This emphasises need for systematic frailty assessment and screening for delirium in acutely ill older patients in hospital and community settings. Clinicians should suspect COVID-19 in frail adults with delirium.

Fasting and time of day independently modulate circadian rhythm relevant gene expression in adipose and skin tissue.

BACKGROUND: Intermittent fasting and time-restricted diets are associated with lower risk biomarkers for cardio-metabolic disease. The shared mechanisms underpinning the similar physiological response to these events is not established, but circadian rhythm could be involved. Here we investigated the transcriptional response to fasting in a large cross-sectional study of adipose and skin tissue from healthy volunteers (N = 625) controlling for confounders of circadian rhythm: time of day and season. RESULTS: We identified 367 genes in adipose and 79 in skin whose expression levels were associated (FDR < 5%) with hours of fasting conditionally independent of time of day and season, with 19 genes common to both tissues. Among these genes, we replicated 38 in human, 157 in non-human studies, and 178 are novel associations. Fasting-responsive genes were enriched for regulation of and response to circadian rhythm. We identified 99 genes in adipose and 54 genes in skin whose expression was associated to time of day; these genes were also enriched for circadian rhythm processes. In genes associated to both exposures the effect of time of day was stronger and in an opposite direction to that of hours fasted. We also investigated the relationship between fasting and genetic regulation of gene expression, including GxE eQTL analysis to identify personal responses to fasting. CONCLUSION: This study robustly implicates circadian rhythm genes in the response to hours fasting independently of time of day, seasonality, age and BMI. We identified tissue-shared and tissue-specific differences in the transcriptional response to fasting in a large sample of healthy volunteers.

Beneficial effect of a high number of copies of salivary amylase AMY1 gene on obesity risk in Mexican children.

AIMS/HYPOTHESIS: Childhood obesity is a major public health problem in Mexico, affecting one in every three children. Genome-wide association studies identified genetic variants associated with childhood obesity, but a large missing heritability remains to be elucidated. We have recently shown a strong association between a highly polymorphic copy number variant encompassing the salivary amylase gene (AMY1 also known as AMY1A) and obesity in European and Asian adults. In the present study, we aimed to evaluate the association between AMY1 copy number and obesity in Mexican children. METHODS: We evaluated the number of AMY1 copies in 597 Mexican children (293 obese children and 304 normal weight controls) through highly sensitive digital PCR. The effect of AMY1 copy number on obesity status was assessed using a logistic regression model adjusted for age and sex. RESULTS: We identified a marked effect of AMY1 copy number on reduced risk of obesity (OR per estimated copy 0.84, with the number of copies ranging from one to 16 in this population; p = 4.25 × 10(-6)). The global association between AMY1 copy number and reduced risk of obesity seemed to be mostly driven by the contribution of the highest AMY1 copy number. Strikingly, all children with >10 AMY1 copies were normal weight controls. CONCLUSIONS/INTERPRETATION: Salivary amylase initiates the digestion of dietary starch, which is highly consumed in Mexico. Our current study suggests putative benefits of high number of AMY1 copies (and related production of salivary amylase) on energy metabolism in Mexican children.

Novel association approach for variable number tandem repeats (VNTRs) identifies DOCK5 as a susceptibility gene for severe obesity.

Variable number tandem repeats (VNTRs) constitute a relatively under-examined class of genomic variants in the context of complex disease because of their sequence complexity and the challenges in assaying them. Recent large-scale genome-wide copy number variant mapping and association efforts have highlighted the need for improved methodology for association studies using these complex polymorphisms. Here we describe the in-depth investigation of a complex region on chromosome 8p21.2 encompassing the dedicator of cytokinesis 5 (DOCK5) gene. The region includes two VNTRs of complex sequence composition which flank a common 3975 bp deletion, all three of which were genotyped by polymerase chain reaction and fragment analysis in a total of 2744 subjects. We have developed a novel VNTR association method named VNTRtest, suitable for association analysis of multi-allelic loci with binary and quantitative outcomes, and have used this approach to show significant association of the DOCK5 VNTRs with childhood and adult severe obesity (P(empirical)= 8.9 × 10(-8) and P= 3.1 × 10(-3), respectively) which we estimate explains ~0.8% of the phenotypic variance. We also identified an independent association between the 3975 base pair (bp) deletion and obesity, explaining a further 0.46% of the variance (P(combined)= 1.6 × 10(-3)). Evidence for association between DOCK5 transcript levels and the 3975 bp deletion (P= 0.027) and both VNTRs (P(empirical)= 0.015) was also identified in adipose tissue from a Swedish family sample, providing support for a functional effect of the DOCK5 deletion and VNTRs. These findings highlight the potential role of DOCK5 in human obesity and illustrate a novel approach for analysis of the contribution of VNTRs to disease susceptibility through association studies.

Insights into the genetic susceptibility to type 2 diabetes from genome-wide association studies of glycaemic traits.

Over the past 8 years, the genetics of complex traits have benefited from an unprecedented advancement in the identification of common variant loci for diseases such as type 2 diabetes (T2D). The ability to undertake genome-wide association studies in large population-based samples for quantitative glycaemic traits has permitted us to explore the hypothesis that models arising from studies in non-diabetic individuals may reflect mechanisms involved in the pathogenesis of diabetes. Amongst 88 T2D risk and 72 glycaemic trait loci, only 29 are shared and show disproportionate magnitudes of phenotypic effects. Important mechanistic insights have been gained regarding the physiological role of T2D loci in disease predisposition through the elucidation of their contribution to glycaemic trait variability. Further investigation is warranted to define causal variants within these loci, including functional characterisation of associated variants, to dissect their role in disease mechanisms and to enable clinical translation.

Adipose tissue eQTL meta-analysis reveals the contribution of allelic heterogeneity to gene expression regulation and cardiometabolic traits.

Complete characterization of the genetic effects on gene expression is needed to elucidate tissue biology and the etiology of complex traits. Here, we analyzed 2,344 subcutaneous adipose tissue samples and identified 34K conditionally distinct expression quantitative trait locus (eQTL) signals in 18K genes. Over half of eQTL genes exhibited at least two eQTL signals. Compared to primary signals, non-primary signals had lower effect sizes, lower minor allele frequencies, and less promoter enrichment; they corresponded to genes with higher heritability and higher tolerance for loss of function. Colocalization of eQTL with conditionally distinct genome-wide association study signals for 28 cardiometabolic traits identified 3,605 eQTL signals for 1,861 genes. Inclusion of non-primary eQTL signals increased colocalized signals by 46%. Among 30 genes with ≥2 pairs of colocalized signals, 21 showed a mediating gene dosage effect on the trait. Thus, expanded eQTL identification reveals more mechanisms underlying complex traits and improves understanding of the complexity of gene expression regulation.

famCNV: copy number variant association for quantitative traits in families.

UNLABELLED: A program package to enable genome-wide association of copy number variants (CNVs) with quantitative phenotypes in families of arbitrary size and complexity. Intensity signals that act as proxies for the number of copies are modeled in a variance component framework and association with traits is assessed through formal likelihood testing. AVAILABILITY AND IMPLEMENTATION: The Java package is made available at www.imperial.ac.uk/medicine/people/m.falchi/. CONTACT: m.falchi@imperial.ac.uk.

Age acquired skewed X chromosome inactivation is associated with adverse health outcomes in humans.

Background: Ageing is a heterogenous process characterised by cellular and molecular hallmarks, including changes to haematopoietic stem cells and is a primary risk factor for chronic diseases. X chromosome inactivation (XCI) randomly transcriptionally silences either the maternal or paternal X in each cell of 46, XX females to balance the gene expression with 46, XY males. Age acquired XCI-skew describes the preferential selection of cells across a tissue resulting in an imbalance of XCI, which is particularly prevalent in blood tissues of ageing females, and yet its clinical consequences are unknown. Methods: We assayed XCI in 1575 females from the TwinsUK population cohort using DNA extracted from whole blood. We employed prospective, cross-sectional, and intra-twin study designs to characterise the relationship of XCI-skew with molecular and cellular measures of ageing, cardiovascular disease risk, and cancer diagnosis. Results: We demonstrate that XCI-skew is independent of traditional markers of biological ageing and is associated with a haematopoietic bias towards the myeloid lineage. Using an atherosclerotic cardiovascular disease risk score, which captures traditional risk factors, XCI-skew is associated with an increased cardiovascular disease risk both cross-sectionally and within XCI-skew discordant twin pairs. In a prospective 10 year follow-up study, XCI-skew is predictive of future cancer incidence. Conclusions: Our study demonstrates that age acquired XCI-skew captures changes to the haematopoietic stem cell population and has clinical potential as a unique biomarker of chronic disease risk. Funding: KSS acknowledges funding from the Medical Research Council [MR/M004422/1 and MR/R023131/1]. JTB acknowledges funding from the ESRC [ES/N000404/1]. MM acknowledges funding from the National Institute for Health Research (NIHR)-funded BioResource, Clinical Research Facility and Biomedical Research Centre based at Guy's and St Thomas' NHS Foundation Trust in partnership with King's College London. TwinsUK is funded by the Wellcome Trust, Medical Research Council, European Union, Chronic Disease Research Foundation (CDRF), Zoe Global Ltd and the National Institute for Health Research (NIHR)-funded BioResource, Clinical Research Facility and Biomedical Research Centre based at Guy's and St Thomas' NHS Foundation Trust in partnership with King's College London.

Combined genetic deletion of GDF15 and FGF21 has modest effects on body weight, hepatic steatosis and insulin resistance in high fat fed mice.

OBJECTIVES: Obesity in humans and mice is associated with elevated levels of two hormones responsive to cellular stress, namely GDF15 and FGF21. Over-expression of each of these is associated with weight loss and beneficial metabolic changes but where they are secreted from and what they are required for physiologically in the context of overfeeding remains unclear. METHODS: Here we used tissue selective knockout mouse models and human transcriptomics to determine the source of circulating GDF15 in obesity. We then generated and characterized the metabolic phenotypes of GDF15/FGF21 double knockout mice. RESULTS: Circulating GDF15 and FGF21 are both largely derived from the liver, rather than adipose tissue or skeletal muscle, in obese states. Combined whole body deletion of FGF21 and GDF15 does not result in any additional weight gain in response to high fat feeding but it does result in significantly greater hepatic steatosis and insulin resistance than that seen in GDF15 single knockout mice. CONCLUSIONS: Collectively the data suggest that overfeeding activates a stress response in the liver which is the major source of systemic rises in GDF15 and FGF21. These hormones then activate pathways which reduce this metabolic stress.

Symptom clusters in COVID-19: A potential clinical prediction tool from the COVID Symptom Study app.

As no one symptom can predict disease severity or the need for dedicated medical support in coronavirus disease 2019 (COVID-19), we asked whether documenting symptom time series over the first few days informs outcome. Unsupervised time series clustering over symptom presentation was performed on data collected from a training dataset of completed cases enlisted early from the COVID Symptom Study Smartphone application, yielding six distinct symptom presentations. Clustering was validated on an independent replication dataset between 1 and 28 May 2020. Using the first 5 days of symptom logging, the ROC-AUC (receiver operating characteristic - area under the curve) of need for respiratory support was 78.8%, substantially outperforming personal characteristics alone (ROC-AUC 69.5%). Such an approach could be used to monitor at-risk patients and predict medical resource requirements days before they are required.

Real-time tracking of self-reported symptoms to predict potential COVID-19.

A total of 2,618,862 participants reported their potential symptoms of COVID-19 on a smartphone-based app. Among the 18,401 who had undergone a SARS-CoV-2 test, the proportion of participants who reported loss of smell and taste was higher in those with a positive test result (4,668 of 7,178 individuals; 65.03%) than in those with a negative test result (2,436 of 11,223 participants; 21.71%) (odds ratio = 6.74; 95% confidence interval = 6.31-7.21). A model combining symptoms to predict probable infection was applied to the data from all app users who reported symptoms (805,753) and predicted that 140,312 (17.42%) participants are likely to have COVID-19.

X chromosome dosage of histone demethylase KDM5C determines sex differences in adiposity.

Males and females differ in body composition and fat distribution. Using a mouse model that segregates gonadal sex (ovaries and testes) from chromosomal sex (XX and XY), we showed that XX chromosome complement in combination with a high-fat diet led to enhanced weight gain in the presence of male or female gonads. We identified the genomic dosage of Kdm5c, an X chromosome gene that escapes X chromosome inactivation, as a determinant of the X chromosome effect on adiposity. Modulating Kdm5c gene dosage in XX female mice to levels that are normally present in males resulted in reduced body weight, fat content, and food intake to a degree similar to that seen with altering the entire X chromosome dosage. In cultured preadipocytes, the levels of KDM5C histone demethylase influenced chromatin accessibility (ATAC-Seq), gene expression (RNA-Seq), and adipocyte differentiation. Both in vitro and in vivo, Kdm5c dosage influenced gene expression involved in extracellular matrix remodeling, which is critical for adipocyte differentiation and adipose tissue expansion. In humans, adipose tissue KDM5C mRNA levels and KDM5C genetic variants were associated with body mass. These studies demonstrate that the sex-dependent dosage of Kdm5c contributes to male/female differences in adipocyte biology and highlight X-escape genes as a critical component of female physiology.

ACE2 expression in adipose tissue is associated with COVID-19 cardio-metabolic risk factors and cell type composition.

COVID-19 severity has varied widely, with demographic and cardio-metabolic factors increasing risk of severe reactions to SARS-CoV-2 infection, but the underlying mechanisms for this remain uncertain. We investigated phenotypic and genetic factors associated with subcutaneous adipose tissue expression of Angiotensin I Converting Enzyme 2 ( ACE2 ), which has been shown to act as a receptor for SARS-CoV-2 cellular entry. In a meta-analysis of three independent studies including up to 1,471 participants, lower adipose tissue ACE2 expression was associated with adverse cardio-metabolic health indices including type 2 diabetes (T2D) and obesity status, higher serum fasting insulin and BMI, and lower serum HDL levels (P<5.32x10 -4 ). ACE2 expression levels were also associated with estimated proportions of cell types in adipose tissue; lower ACE2 expression was associated with a lower proportion of microvascular endothelial cells (P=4.25x10 -4 ) and higher macrophage proportion (P=2.74x10 -5 ), suggesting a link to inflammation. Despite an estimated heritability of 32%, we did not identify any proximal or distal genetic variants (eQTLs) associated with adipose tissue ACE2 expression. Our results demonstrate that at-risk individuals have lower background ACE2 levels in this highly relevant tissue. Further studies will be required to establish how this may contribute to increased COVID-19 severity.

Epigenetic findings in periodontitis in UK twins: a cross-sectional study.

BACKGROUND: Genetic and environmental risk factors contribute to periodontal disease, but the underlying susceptibility pathways are not fully understood. Epigenetic mechanisms are malleable regulators of gene function that can change in response to genetic and environmental stimuli, thereby providing a potential mechanism for mediating risk effects in periodontitis. The aim of this study is to identify epigenetic changes across tissues that are associated with periodontal disease. METHODS: Self-reported gingival bleeding and history of gum disease, or tooth mobility, were used as indicators of periodontal disease. DNA methylation profiles were generated using the Infinium HumanMethylation450 BeadChip in whole blood, buccal, and adipose tissue samples from predominantly older female twins (mean age 58) from the TwinsUK cohort. Epigenome-wide association scans (EWAS) of gingival bleeding and tooth mobility were conducted in whole blood in 528 and 492 twins, respectively. Subsequently, targeted candidate gene analysis at 28 genomic regions was carried out testing for phenotype-methylation associations in 41 (tooth mobility) and 43 (gingival bleeding) buccal, and 501 (tooth mobility) and 556 (gingival bleeding) adipose DNA samples. RESULTS: Epigenome-wide analyses in blood identified one CpG-site (cg21245277 in ZNF804A) associated with gingival bleeding (FDR = 0.03, nominal p value = 7.17e-8) and 58 sites associated with tooth mobility (FDR < 0.05) with the top signals in IQCE and XKR6. Epigenetic variation at 28 candidate regions (247 CpG-sites) for chronic periodontitis showed an enrichment for association with periodontal traits, and signals in eight genes (VDR, IL6ST, TMCO6, IL1RN, CD44, IL1B, WHAMM, and CXCL1) were significant in both traits. The methylation-phenotype association signals validated in buccal samples, and a subset (25%) also validated in adipose tissue. CONCLUSIONS: Epigenome-wide analyses in adult female twins identified specific DNA methylation changes linked to self-reported periodontal disease. Future work will explore the environmental basis and functional impact of these results to infer potential for strategic personalized treatments and prevention of chronic periodontitis.

Regulatory variants at KLF14 influence type 2 diabetes risk via a female-specific effect on adipocyte size and body composition.

Individual risk of type 2 diabetes (T2D) is modified by perturbations to the mass, distribution and function of adipose tissue. To investigate the mechanisms underlying these associations, we explored the molecular, cellular and whole-body effects of T2D-associated alleles near KLF14. We show that KLF14 diabetes-risk alleles act in adipose tissue to reduce KLF14 expression and modulate, in trans, the expression of 385 genes. We demonstrate, in human cellular studies, that reduced KLF14 expression increases pre-adipocyte proliferation but disrupts lipogenesis, and in mice, that adipose tissue-specific deletion of Klf14 partially recapitulates the human phenotype of insulin resistance, dyslipidemia and T2D. We show that carriers of the KLF14 T2D risk allele shift body fat from gynoid stores to abdominal stores and display a marked increase in adipocyte cell size, and that these effects on fat distribution, and the T2D association, are female specific. The metabolic risk associated with variation at this imprinted locus depends on the sex both of the subject and of the parent from whom the risk allele derives.