Toby Andrew

This study investigated the potential genetic mechanisms which underlie adipose tissue mitochondrial dysfunction in Type 2 diabetes (T2D), by systematically identifying nuclear-encoded mitochondrial genes (NEMGs) among the genes regulated by T2D-associated genetic loci. The target genes of these ‘disease loci’ were identified by mapping genetic loci associated with both disease and gene expression levels (expression quantitative trait loci, eQTL) using high resolution genetic maps, with independent estimates co-locating to within a small genetic distance. These co-locating signals were defined as T2D-eQTL and the target genes as T2D cis-genes. In total, 763 cis-genes were associated with T2D-eQTL, of which 50 were NEMGs. Independent gene expression datasets for T2D and insulin resistant cases and controls confirmed that the cis-genes and cis-NEMGs were enriched for differential expression in cases, providing independent validation that genetic maps can identify informative functional genes. Two additional results were consistent with a potential role of T2D-eQTL in regulating the 50 identified cis-NEMGs in the context of T2D risk: (1) the 50 cis-NEMGs showed greater differential expression compared to other NEMGs and (2) other NEMGs showed a trend towards significantly decreased expression if their expression levels correlated more highly with the subset of 50 cis-NEMGs. These 50 cis-NEMGs, which are differentially expressed and associated with mapped T2D disease loci, encode proteins acting within key mitochondrial pathways, including some of current therapeutic interest such as the metabolism of branched-chain amino acids, GABA and biotin.

Background: Adipogenesis, the process whereby preadipocytes differentiate into mature adipocytes, is crucial for maintaining metabolic homeostasis. Cholesterol lowering statins increase type 2 diabetes (T2D) risk possibly by affecting adipogenesis and insulin resistance but the (epi)genetic mechanisms involved are unknown. Here, we characterised the effects of statin treatment on adipocyte differentiation using in vitro human preadipocytes cell model to identify putative effective genes.Results: Statin treatment during adipocyte differentiation caused a reduction in key genes involved in adipogenesis, such as ADIPOQ, GLUT4 and ABCG1. Using Illumina’s Infinium ‘850K’ Methylation EPIC array, we found a significant hypomethylation of cg14566882, located in the promoter of the histone deacetylase 9 (HDAC9) gene, in response to two types of statins (atorvastatin and mevastatin), which correlates with an increased HDAC9 mRNA expression. HDAC9 is a transcriptional repressor of the cholesterol efflux ABCG1 gene expression, which is epigenetically modified in obesity and prediabetic states. Thus, we assessed the putative impact of ABCG1 knockdown in mimicking the effect of statin in adipogenesis. ABCG1 KD reduced the expression of key genes involved in adipocyte differentiation and decreased insulin signalling and glucose uptake. In human blood cells from two cohorts, ABCG1 expression was impaired in response to statins, confirming that ABCG1 is in vivo targeted by these drugs.Conclusions: We identified an epigenetic link between adipogenesis and adipose tissue insulin resistance in the context of T2D risk associated with statin use, which has important implications as HDAC9 and ABCG1 are considered potential therapeutic targets for obesity and metabolic diseases.

Numerous functional studies have implicated PARL in relation to type 2 diabetes (T2D). We hypothesised that conflicting human association studies may be due to neighbouring causal variants being in linkage disequilibrium (LD) with PARL. We conducted a comprehensive candidate gene study of the extended LD genomic region that includes PARL and transporter ABCC5 using three data sets (two European and one African American), in relation to healthy glycaemic variation, visceral fat accumulation and T2D disease. We observed no evidence for previously reported T2D association with Val262Leu or PARL using array and fine-map genomic and expression data. By contrast, we observed strong evidence of T2D association with ABCC5 (intron 26) for European and African American samples (P = 3E−07) and with ABCC5 adipose expression in Europeans [odds ratio (OR) = 3.8, P = 2E−04]. The genomic location estimate for the ABCC5 functional variant, associated with all phenotypes and expression data (P = 1E−11), was identical for all samples (at Chr3q 185,136 kb B36), indicating that the risk variant is an expression quantitative trait locus (eQTL) with increased expression conferring risk of disease. That the association with T2D is observed in populations of disparate ancestry suggests the variant is a ubiquitous risk factor for T2D.

Observational studies have consistently shown associations between mild deficiencies in folate and vitamin B₁₂ with increased risk of a myriad of common diseases. These findings have invariably translated into null outcomes in intervention trials due in part to our ignorance of the specific genomic and environmental factors that underpin population variability in requirements to these B-vitamins. Although genome-wide association studies have shed initial light on the genetic architecture of variability in status of these vitamins, particularly vitamin B₁₂, the causal mechanisms remain uncharacterised. A recent study by Grarup et al. (PLoS Genet 9(6):e1003530, 2013) used next-generation whole-genome sequencing to gain further insight into the genetic architecture of vitamin B₁₂ and folate status in the general population. Their study represents the analysis of approximately ten times greater number of genetic variants and nearly four times the number of individuals compared to the largest previous GWAS study of these B-vitamins. In light of this, we purport that although the study may be viewed as the state of the art in the roadmap to personalised or precision nutrition, the lack of insight provided by the study serves as a cautionary reminder of the importance of study design, particularly when leveraging large-scale data, such as those from whole-genome sequences. We believe that the precedent set by such large-scale “proof of principle” type projects will wrongly enforce a negative outlook for nutrigenomics research and present alternative study designs, which although less opportunistic are far more likely to be informative and yield novel results.

Cobalamin (Cbl) is an essential B vitamin involved in the normal functioning of the nervous system, the formation of key components of blood, DNA synthesis and methylation, and energy production. Physiological levels of Cbl vary greatly within populations, although the basis for this variability remains largely unknown. We conducted a twin study to characterise the basis of variation in plasma Cbl levels and to test whether common genetic polymorphisms in genes known to cause defects in inborn errors of Cbl metabolism and transport are also associated with mean plasma Cbl levels in the general population. The present results showed that plasma levels of Cbl were heritable, with genetic and phenotypic variance increasing with age, and levels significantly correlated with age, BMI, exercise, alcohol consumption, smoking status, social class and folate levels, which collectively accounted for up to 15 % of Cbl variation. Of eight genes responsible for the defects of the Cbl metabolic pathway (cblA–G and mut), MMAA, MMACHC, MTRR and MUT harboured polymorphisms that showed evidence of association with Cbl levels. Characterisation of the heritable component of variation in Cbl levels can facilitate the early diagnosis and prognosis of Cbl insufficiency/deficiency in individuals at a higher risk of associated diseases.

Background: Excess accumulation of visceral fat is a prominent risk factor for cardiovascular and metabolic morbidity. While computed tomography (CT) is the gold standard to measure visceral adiposity, this is often not possible for large studies - thus valid, but less expensive and intrusive proxy measures of visceral fat are required such as dual-energy X-ray absorptiometry (DXA). Study aims were to a) identify a valid DXA-based measure of visceral adipose tissue (VAT), b) estimate VAT heritability and c) assess visceral fat association with morbidity in relation to body fat distribution. Methods: A validation sample of 54 females measured for detailed body fat composition - assessed using CT, DXA and anthropometry – was used to evaluate previously published predictive models of CT-measured visceral fat. Based upon a validated model, we realised an out-of-sample estimate of abdominal VAT area for a study sample of 3457 female volunteer twins and estimated VAT area heritability using a classical twin study design. Regression and residuals analyses were used to assess the relationship between adiposity and morbidity. Results: Published models applied to the validation sample explained >80% of the variance in CT-measured visceral fat. While CT visceral fat was best estimated using a linear regression for waist circumference, CT body cavity area and total abdominal fat (R ² = 0.91), anthropometric measures alone predicted VAT almost equally well (CT body cavity area and waist circumference, R ² = 0.86). Narrow sense VAT area heritability for the study sample was estimated to be 58% (95% CI: 51-66%) with a shared familial component of 24% (17-30%). VAT area is strongly associated with type 2 diabetes (T2D), hypertension (HT), subclinical atherosclerosis and liver function tests. In particular, VAT area is associated with T2D, HT and liver function (alanine transaminase) independent of DXA total abdominal fat and body mass index (BMI). Conclusions: DXA and anthropometric measures can be utilised to derive estimates of visceral fat as a reliable alternative to CT. Visceral fat is heritable and appears to mediate the association between body adiposity and morbidity. This observation is consistent with hypotheses that suggest excess visceral adiposity is causally related to cardiovascular and metabolic disease.

Purpose.: Astigmatism is a common refractive error that reduces vision, where the curvature and refractive power of the cornea in one meridian are less than those of the perpendicular axis. It is a complex trait likely to be influenced by both genetic and environmental factors. Twin studies of astigmatism have found approximately 60% of phenotypic variance is explained by genetic factors. This study aimed to identify susceptibility loci for astigmatism. Methods.: We performed a meta-analysis of seven genome-wide association studies that included 22,100 individuals of European descent, where astigmatism was defined as the number of diopters of cylinder prescription, using fixed effect inverse variance-weighted methods. Results.: A susceptibility locus was identified with lead single nucleotide polymorphism rs3771395 on chromosome 2p13.3 (meta-analysis, P = 1.97 × 10⁻⁷) in the VAX2 gene. VAX2 plays an important role in the development of the dorsoventral axis of the eye. Animal studies have shown a gradient in astigmatism along the vertical plane, with corresponding changes in refraction, particularly in the ventral field. Conclusions.: This finding advances the understanding of refractive error, and provides new potential pathways to be evaluated with regard to the development of astigmatism.

The mitochondrial theory of ageing proposes that damage to mitochondria and diminished mitochondrial DNA (mtDNA) repair are major contributors to cellular dysfunction and age-related diseases. We investigate the prevalence of heteroplasmy in the mtDNA control region in buccal swab and blood derived samples for 178 women from the TwinsUK cohort (41 DZ pair 39 MZ pairs, 18 singletons, mean age 57.5 range 28–82) and its relationship to age, BMI and fasting insulin and glucose serum levels. The overall estimated prevalence of heteroplasmy for both tissues in the control region measured for 37 sites was 17%. The prevalence of heteroplasmy was higher among the older half of the study subjects than in the younger half (23% vs 10% p−11). The presence of the C allele at position 16093 in blood was associated with the presence of heteroplasmy in buccal swabs at this position (p = 3.5×10⁻¹⁴) and also at VRII (p = 2×10⁻⁴) suggesting a possible predisposing role for this site in the accumulation of heteroplasmy. Our data indicate that BMI is potentially associated with control region heteroplasmy.