Analyses of potential causal contributors to increased waist/hip ratio-associated cardiometabolic disease: A combined and sex-stratified Mendelian randomization study.

BACKGROUND: Increased waist/hip ratio (WHR) contributes to type 2 diabetes, fatty liver, dyslipidaemia, hypertension and coronary artery disease, with potential sex-differential effects. Postulated mediators include increased lipid flux, branched-chain amino acids, glycine and glycoprotein acetyl, but their relative contributions and sex-specific impact on WHR-associated cardiometabolic disease (CMD) are not established. METHODS: We therefore undertook combined and sex-stratified Mendelian randomization (MR) to assess the relative causal contributions of these mediators to WHR-associated CMD using summary statistics from the largest genome-wide association studies in European ancestries. RESULTS: In sex-combined MR analyses, increased WHR significantly reduces high-density lipoprotein (beta = -0.416, SE = 0.029, p = 2.87E-47), increases triglyceride (beta = 0.431, SE = 0.029, p = 1.87E-50), type 2 diabetes (odds ratio = 2.747, SE = 0.09, p = 26E-23), coronary artery disease (odds ratio = 1.478, SE = 0.045, p = 6.96E-18), alanine transaminase (beta = 0.062, SE = 0.004, p = 6.88E-67), and systolic (beta = 0.134, SE = 0.022, p = 7.81E-10) and diastolic blood pressure (beta = 0.162, SE = 0.026, p = 5.38E-10). Adjustment for the mediators attenuated WHR's effects, but the associations remained significant with concordant results in females. In males, a similar pattern was seen, except after adjusting for the effect of the ratio of monounsaturated fatty acid to total free fatty acid, the potential causal effect of WHR was no longer significant: high-density lipoprotein (beta = -0.117, SE = 0.069, p = .09) and triglyceride (beta = 0.051, SE = 0.068, p = .459). CONCLUSIONS: MR suggests WHR increases the risk of CMD independent of these mediators, with the exception of dyslipidaemia in males, which is largely driven by the monounsaturated fatty acid to total free fatty acid ratio.

Comprehensive whole-genome analyses of the UK Biobank reveal significant sex differences in both genotype missingness and allele frequency on the X chromosome.

The UK Biobank is the most used dataset for genome-wide association studies (GWAS). GWAS of sex, essentially sex differences in minor allele frequencies (sdMAF), has identified autosomal SNPs with significant sdMAF, including in the UK Biobank, but the X chromosome was excluded. Our recent report identified multiple regions on the X chromosome with significant sdMAF, using short-read sequencing of other datasets. We performed a whole genome sdMAF analysis, with ~410 k white British individuals from the UK Biobank, using array genotyped, imputed or exome sequencing data. We observed marked sdMAF on the X chromosome, particularly at the boundaries between the pseudo-autosomal regions (PAR) and the non-PAR (NPR), as well as throughout the NPR, consistent with our earlier report. A small fraction of autosomal SNPs also showed significant sdMAF. Using the centrally imputed data, which relied mostly on low-coverage whole genome sequence, resulted in 2.1% of NPR SNPs with significant sdMAF. The whole exome sequencing also displays sdMAF on the X chromosome, including some NPR SNPs with heterozygous genotype calls in males. Genotyping, sequencing and imputation of X chromosomal SNPs requires further attention to ensure the integrity for downstream association analysis.

Beyond apples and pears: sex-specific genetics of body fat percentage.

Introduction: Biological sex influences both overall adiposity and fat distribution. Further, testosterone and sex hormone binding globulin (SHBG) influence adiposity and metabolic function, with differential effects of testosterone in men and women. Here, we aimed to perform sex-stratified genome-wide association studies (GWAS) of body fat percentage (BFPAdj) (adjusting for testosterone and sex hormone binding globulin (SHBG)) to increase statistical power. Methods: GWAS were performed in white British individuals from the UK Biobank (157,937 males and 154,337 females). To avoid collider bias, loci associated with SHBG or testosterone were excluded. We investigated association of BFPAdj loci with high density cholesterol (HDL), triglyceride (TG), type 2 diabetes (T2D), coronary artery disease (CAD), and MRI-derived abdominal subcutaneous adipose tissue (ASAT), visceral adipose tissue (VAT) and gluteofemoral adipose tissue (GFAT) using publicly available data from large GWAS. We also performed 2-sample Mendelian Randomization (MR) using identified BFPAdj variants as instruments to investigate causal effect of BFPAdj on HDL, TG, T2D and CAD in males and females separately. Results: We identified 195 and 174 loci explaining 3.35% and 2.60% of the variation in BFPAdj in males and females, respectively at genome-wide significance (GWS, p<5x10-8). Although the direction of effect at these loci was generally concordant in males and females, only 38 loci were common to both sexes at GWS. Seven loci in males and ten loci in females have not been associated with any adiposity/cardiometabolic traits previously. BFPAdj loci generally did not associate with cardiometabolic traits; several had paradoxically beneficial cardiometabolic effects with favourable fat distribution. MR analyses did not find convincing supportive evidence that increased BFPAdj has deleterious cardiometabolic effects in either sex with highly significant heterogeneity. Conclusions: There was limited genetic overlap between BFPAdj in males and females at GWS. BFPAdj loci generally did not have adverse cardiometabolic effects which may reflect the effects of favourable fat distribution and cardiometabolic risk modulation by testosterone and SHBG.

SLCO5A1 and synaptic assembly genes contribute to impulsivity in juvenile myoclonic epilepsy.

Elevated impulsivity is a key component of attention-deficit hyperactivity disorder (ADHD), bipolar disorder and juvenile myoclonic epilepsy (JME). We performed a genome-wide association, colocalization, polygenic risk score, and pathway analysis of impulsivity in JME (n = 381). Results were followed up with functional characterisation using a drosophila model. We identified genome-wide associated SNPs at 8q13.3 (P = 7.5 × 10-9) and 10p11.21 (P = 3.6 × 10-8). The 8q13.3 locus colocalizes with SLCO5A1 expression quantitative trait loci in cerebral cortex (P = 9.5 × 10-3). SLCO5A1 codes for an organic anion transporter and upregulates synapse assembly/organisation genes. Pathway analysis demonstrates 12.7-fold enrichment for presynaptic membrane assembly genes (P = 0.0005) and 14.3-fold enrichment for presynaptic organisation genes (P = 0.0005) including NLGN1 and PTPRD. RNAi knockdown of Oatp30B, the Drosophila polypeptide with the highest homology to SLCO5A1, causes over-reactive startling behaviour (P = 8.7 × 10-3) and increased seizure-like events (P = 6.8 × 10-7). Polygenic risk score for ADHD genetically correlates with impulsivity scores in JME (P = 1.60 × 10-3). SLCO5A1 loss-of-function represents an impulsivity and seizure mechanism. Synaptic assembly genes may inform the aetiology of impulsivity in health and disease.

Investigating the association between fasting insulin, erythrocytosis and HbA1c through Mendelian randomization and observational analyses.

Background: Insulin resistance (IR) with associated compensatory hyperinsulinemia (HI) are early abnormalities in the etiology of prediabetes (preT2D) and type 2 diabetes (T2D). IR and HI also associate with increased erythrocytosis. Hemoglobin A1c (HbA1c) is commonly used to diagnose and monitor preT2D and T2D, but can be influenced by erythrocytosis independent of glycemia. Methods: We undertook bidirectional Mendelian randomization (MR) in individuals of European ancestry to investigate potential causal associations between increased fasting insulin adjusted for BMI (FI), erythrocytosis and its non-glycemic impact on HbA1c. We investigated the association between the triglyceride-glucose index (TGI), a surrogate measure of IR and HI, and glycation gap (difference between measured HbA1c and predicted HbA1c derived from linear regression of fasting glucose) in people with normoglycemia and preT2D. Results: Inverse variance weighted MR (IVWMR) suggested that increased FI increases hemoglobin (Hb, b=0.54 ± 0.09, p=2.7 x 10-10), red cell count (RCC, b=0.54 ± 0.12, p=5.38x10-6) and reticulocyte (RETIC, b=0.70 ± 0.15, p=2.18x10-6). Multivariable MR indicated that increased FI did not impact HbA1c (b=0.23 ± 0.16, p=0.162) but reduced HbA1c after adjustment for T2D (b=0.31 ± 0.13, p=0.016). Increased Hb (b=0.03 ± 0.01, p=0.02), RCC (b=0.02 ± 0.01, p=0.04) and RETIC (b=0.03 ± 0.01, p=0.002) might modestly increase FI. In the observational cohort, increased TGI associated with decreased glycation gap, (i.e., measured HbA1c was lower than expected based on fasting glucose, (b=-0.09 ± 0.009, p<0.0001)) in people with preT2D but not in those with normoglycemia (b=0.02 ± 0.007, p<0.0001). Conclusions: MR suggests increased FI increases erythrocytosis and might potentially decrease HbA1c by non-glycemic effects. Increased TGI, a surrogate measure of increased FI, associates with lower-than-expected HbA1c in people with preT2D. These findings merit confirmatory studies to evaluate their clinical significance.

Association Between Obesity and Chronic Kidney Disease: Multivariable Mendelian Randomization Analysis and Observational Data From a Bariatric Surgery Cohort.

Obesity is postulated to independently increase chronic kidney disease (CKD), even after adjusting for type 2 diabetes (T2D) and hypertension. Dysglycemia below T2D thresholds, frequently seen with obesity, also increases CKD risk. Whether obesity increases CKD independent of dysglycemia and hypertension is unknown and likely influences the optimal weight loss (WL) needed to reduce CKD. T2D remission rates plateau with 20-25% WL after bariatric surgery (BS), but further WL increases normoglycemia and normotension. We undertook bidirectional inverse variance weighted Mendelian randomization (IVWMR) to investigate potential independent causal associations between increased BMI and estimated glomerular filtration rate (eGFR) in CKD (CKDeGFR) (<60 mL/min/1.73 m2) and microalbuminuria (MA). In 5,337 BS patients, we assessed whether WL influences >50% decline in eGFR (primary outcome) or CKD hospitalization (secondary outcome), using <20% WL as a comparator. IVWMR results suggest that increased BMI increases CKDeGFR (b = 0.13, P = 1.64 × 10-4; odds ratio [OR] 1.14 [95% CI 1.07, 1.23]) and MA (b = 0.25; P = 2.14 × 10-4; OR 1.29 [1.13, 1.48]). After adjusting for hypertension and fasting glucose, increased BMI did not significantly increase CKDeGFR (b = -0.02; P = 0.72; OR 0.98 [0.87, 1.1]) or MA (b = 0.19; P = 0.08; OR 1.21 [0.98, 1.51]). Post-BS WL significantly reduced the primary outcome with 30 to <40% WL (hazard ratio [HR] 0.53 [95% CI 0.32, 0.87]) but not 20 to <30% WL (HR 0.72 [0.44, 1.2]) and ≥40% WL (HR 0.73 [0.41, 1.30]). For CKD hospitalization, progressive reduction was seen with increased WL, which was significant for 30 to <40% WL (HR 0.37 [0.17, 0.82]) and ≥40% WL (HR 0.24 [0.07, 0.89]) but not 20 to <30% WL (HR 0.60 [0.29, 1.23]). The data suggest that obesity is likely not an independent cause of CKD. WL thresholds previously associated with normotension and normoglycemia, likely causal mediators, may reduce CKD after BS.

Role of 19 SNPs in 10 genes with type 2 diabetes in the Pakistani population.

Meta-GWAS report numerous variants associated with type 2 diabetes (T2D), however, for diabetic retinopathy (DR) no loci achieved genome-wide significance. There are limited candidate gene analyses for T2D and/or DR reported from the Pakistani population. Therefore, the current study was designed to evaluate the genetic association of 10 loci with T2D, non-proliferative DR (NPDR), and proliferative DR (PDR). In total 375 T2D cases and 205 controls were collected. The T2D cases included diabetic no retinopathy (n = 196), NPDR (n = 95), and PDR (n = 84). Genomic DNA was isolated, and 19 SNPs were genotyped. To determine association of SNPs with T2D, logistic regression analyses were performed adjusting for age and sex. Moreover, for association of SNPs with NPDR and PDR logistic regression analyses adjusting for diabetes duration and age of T2D onset were performed. In multivariate analysis, the minor alleles of rs1043618 [G > C, odds ratio (OR) 95% confidence interval (CI) = 1.45 (1.13-1.87), p = 4.00E-3], rs3807987 [G > A, 1.87 (1.22-2.94), p = 0.01], rs12672038 [G > A, 1.53 (1.04-2.30), p = 0.03] and rs2055858 [G > C, 1.70 (1.20-2.43), p = 3.00E-3] were associated with higher risk while rs1801133 (C > T, 0.59 (0.42-0.83), p = 2.28E-3) was associated with a lower risk of T2D. Moreover, minor alleles of rs2055858 [G > C, 1.77 (1.17-2.68), p = 0.02], and rs3759890 [C > G, 2.17 (1.39-3.39), p = 4.00E-3] showed an association with PDR when compared with DNR. However, only the association of rs1801133 survived multiple test correction. Hence, we report that rs1801133 is associated with T2D in the Pakistani population. In addition, out of studied 10 genes 8 proteins had higher interactions among themselves that are predicted to be partially biologically connected, as a group.

Genome-wide association study identifies novel loci associated with skin autofluorescence in individuals without diabetes.

BACKGROUND: Skin autofluorescence (SAF) is a non-invasive measure reflecting accumulation of advanced glycation endproducts (AGEs) in the skin. Higher SAF levels are associated with an increased risk of developing type 2 diabetes and cardiovascular disease. An earlier genome-wide association study (GWAS) revealed a strong association between NAT2 variants and SAF. The aim of this study was to calculate SAF heritability and to identify additional genetic variants associated with SAF through genome-wide association studies (GWAS). RESULTS: In 27,534 participants without diabetes the heritability estimate of lnSAF was 33% ± 2.0% (SE) in a model adjusted for covariates. In meta-GWAS for lnSAF five SNPs, on chromosomes 8, 11, 15 and 16 were associated with lnSAF (P < 5 × 10-8): 1. rs2846707 (Chr11:102,576,358,C > T), which results in a Met30Val missense variant in MMP27 exon 1 (NM_022122.3); 2. rs2470893 (Chr15:75,019,449,C > T), in intergenic region between CYP1A1 and CYP1A2; with attenuation of the SNP-effect when coffee consumption was included as a covariate; 3. rs12931267 (Chr16:89,818,732,C > G) in intron 30 of FANCA and near MC1R; and following conditional analysis 4. rs3764257 (Chr16:89,800,887,C > G) an intronic variant in ZNF276, 17.8 kb upstream from rs12931267; finally, 30 kb downstream from NAT2 5. rs576201050 (Chr8:18,288,053,G > A). CONCLUSIONS: This large meta-GWAS revealed five SNPs at four loci associated with SAF in the non-diabetes population. Further unravelling of the genetic architecture of SAF will help in improving its utility as a tool for screening and early detection of diseases and disease complications.

Genomic architecture of autism from comprehensive whole-genome sequence annotation.

Fully understanding autism spectrum disorder (ASD) genetics requires whole-genome sequencing (WGS). We present the latest release of the Autism Speaks MSSNG resource, which includes WGS data from 5,100 individuals with ASD and 6,212 non-ASD parents and siblings (total n = 11,312). Examining a wide variety of genetic variants in MSSNG and the Simons Simplex Collection (SSC; n = 9,205), we identified ASD-associated rare variants in 718/5,100 individuals with ASD from MSSNG (14.1%) and 350/2,419 from SSC (14.5%). Considering genomic architecture, 52% were nuclear sequence-level variants, 46% were nuclear structural variants (including copy-number variants, inversions, large insertions, uniparental isodisomies, and tandem repeat expansions), and 2% were mitochondrial variants. Our study provides a guidebook for exploring genotype-phenotype correlations in families who carry ASD-associated rare variants and serves as an entry point to the expanded studies required to dissect the etiology in the ∼85% of the ASD population that remain idiopathic.

Insulin Response to Oral Glucose and Cardiometabolic Disease: A Mendelian Randomization Study to Assess Potential Causality.

Mendelian randomization (MR) suggests that postprandial hyperinsulinemia (unadjusted for plasma glucose) increases BMI, but its impact on cardiometabolic disease, a leading cause for mortality and morbidity in people with obesity, is not established. Fat distribution i.e., increased centripetal and/or reduced femoro-gluteal adiposity, is causally associated with and better predicts cardiometabolic disease than BMI. We therefore undertook bidirectional MR to assess the effect of corrected insulin response (CIR) (insulin 30 min after a glucose challenge adjusted for plasma glucose) on BMI, waist-to-hip ratio (WHR), leg fat, type 2 diabetes (T2D), triglyceride (TG), HDL, liver fat, hypertension (HTN), and coronary artery disease (CAD) in people of European descent. Inverse variance-weighted MR suggests a potential causal association between increased CIR and increased BMI (b = 0.048 ± 0.02, P = 0.03), increased leg fat (b = 0.029 ± 0.012, P = 0.01), reduced T2D (b = -0.73 ± 0.15, P = 6 × 10-7, odds ratio [OR] 0.48 [95% CI 0.36-0.64]), reduced TG (b = -0.07 ± 0.02, P = 0.003), and increased HDL (b = 0.04 ± 0.01, P = 0.006) with some evidence of horizontal pleiotropy. CIR had neutral effects on WHR (b = 0.009 ± 0.02, P = 0.69), liver fat (b = -0.08 ± 0.04, P = 0.06), HTN (b = -0.001 ± 0.004, P = 0.7, OR 1.00 [95% CI 0.99-1.01]), and CAD (b = -0.002 ± 0.002, P = 0.48, OR 0.99 [95% CI 0.81-1.21]). T2D decreased CIR (b -0.22 ± 0.04, P = 1.3 × 10-7), with no evidence that BMI, TG, HDL, liver fat, HTN, and CAD modulate CIR. In conclusion, we did not find evidence that increased CIR increases cardiometabolic disease. It might increase BMI with favorable fat distribution, reduce T2D, and improve lipids.

Cohort profile: genomic data for 26 622 individuals from the Canadian Longitudinal Study on Aging (CLSA).

PURPOSE: The Canadian Longitudinal Study on Aging (CLSA) Comprehensive cohort was established to provide unique opportunities to study the genetic and environmental contributions to human disease as well as ageing process. The aim of this report was to describe the genomic data included in CLSA. PARTICIPANTS: A total of 26 622 individuals from the CLSA Comprehensive cohort of men and women aged 45-85 recruited between 2010 and 2015 underwent genome-wide genotyping of DNA samples collected from blood. Comprehensive quality control metrics were measured for genetic markers and samples, respectively. The genotypes were imputed to the TOPMed reference panel. Sex chromosome abnormalities were identified by copy number profiling. Classical human leukocyte antigen gene haplotypes were imputed at two-field (four-digit). FINDINGS TO DATE: Of the 26 622 genotyped participants, 24 655 (92.6%) were identified as having European ancestry. These genomic data were linked to physical, lifestyle, medical, economic, environmental and psychosocial factors collected longitudinally in CLSA. The combined analysis, including CLSA genomic data, uncovered over 100 novel loci associated with key parameters to define glaucoma. The CLSA genomic dataset validated the contribution of a polygenic risk score to screen individuals with high fracture risk. It is also a valuable resource to directly identify common genetic variations associated with conditions related to complex traits. Taking advantage of the comprehensive interview and physical information collected in CLSA, this genomic dataset has been linked to psychosocial factors to investigate both the independent and interactive effects on cardiovascular disease. FUTURE PLANS: The CLSA overall is ongoing. Follow-up data will continue to be collected from participants in the current genomic subcohort, including the DNA methylation and metabolomic data. Ongoing studies focus on elucidating the role of genetic factors in cognitive decline and cardiovascular diseases. This genomic data resource is available on request through the CLSA data access application process.

A recurrent SHANK3 frameshift variant in Autism Spectrum Disorder.

Autism Spectrum Disorder (ASD) is genetically complex with ~100 copy number variants and genes involved. To try to establish more definitive genotype and phenotype correlations in ASD, we searched genome sequence data, and the literature, for recurrent predicted damaging sequence-level variants affecting single genes. We identified 18 individuals from 16 unrelated families carrying a heterozygous guanine duplication (c.3679dup; p.Ala1227Glyfs*69) occurring within a string of 8 guanines (genomic location [hg38]g.50,721,512dup) affecting SHANK3, a prototypical ASD gene (0.08% of ASD-affected individuals carried the predicted p.Ala1227Glyfs*69 frameshift variant). Most probands carried de novo mutations, but five individuals in three families inherited it through somatic mosaicism. We scrutinized the phenotype of p.Ala1227Glyfs*69 carriers, and while everyone (17/17) formally tested for ASD carried a diagnosis, there was the variable expression of core ASD features both within and between families. Defining such recurrent mutational mechanisms underlying an ASD outcome is important for genetic counseling and early intervention.

Residual ß cell function in long-term type 1 diabetes associates with reduced incidence of hypoglycemia.

BACKGROUNDWe investigated residual ß cell function in Diabetes Control and Complications Trial/Epidemiology of Diabetes Interventions and Complications (DCCT/EDIC) study participants with an average 35-year duration of type 1 diabetes mellitus (T1DM).METHODSSerum C-peptide was measured during a 4-hour mixed-meal tolerance test. Associations with metabolic outcomes and complications were explored among nonresponders (all C-peptide values after meal <0.003 nmol/L) and 3 categories of responders, classified by peak C-peptide concentration (nmol/L) as high (>0.2), intermediate (>0.03 to ≤0.2), and low (≥ 0.003 to ≤0.03).RESULTSOf the 944 participants, 117 (12.4%) were classified as responders. Residual C-peptide concentrations were associated with higher DCCT baseline concentrations of stimulated C-peptide (P value for trend = 0.0001). Residual C-peptide secretion was not associated with current or mean HbA1c, HLA high-risk haplotypes for T1DM, or the current presence of T1DM autoantibodies. The proportion of subjects with a history of severe hypoglycemia was lower with high (27%) and intermediate (48%) residual C-peptide concentrations than with low (74%) and no (70%) residual C-peptide concentrations (P value for trend = 0.0001). Responders and nonresponders demonstrated similar rates of advanced microvascular complications.CONCLUSIONß Cell function can persist in long-duration T1DM. With a peak C-peptide concentration of >0.03 nmol/L, we observed clinically meaningful reductions in the prevalence of severe hypoglycemia.TRIAL REGISTRATIONClinicalTrials.gov NCT00360815 and NCT00360893.FUNDINGDivision of Diabetes Endocrinology and Metabolic Diseases of the National Institute of Diabetes and Digestive and Kidney Diseases (DP3-DK104438, U01 DK094176, and U01 DK094157).

DNA methylation mediates development of HbA1c-associated complications in type 1 diabetes.

Metabolic memory, the persistent benefits of early glycaemic control on preventing and/or delaying the development of diabetic complications, has been observed in the Diabetes Control and Complications Trial (DCCT) and in the Epidemiology of Diabetes Interventions and Complications (EDIC) follow-up study, but the underlying mechanisms remain unclear. Here, we show the involvement of epigenetic DNA methylation (DNAme) in metabolic memory by examining its associations with preceding glycaemic history, and with subsequent development of complications over an 18-yr period in the blood DNA of 499 randomly selected DCCT participants with type 1 diabetes who are also followed up in EDIC. We demonstrate the associations between DNAme near the closeout of DCCT and mean HbA1c during DCCT (mean-DCCT HbA1c) at 186 cytosine-guanine dinucleotides (CpGs) (FDR < 15%, including 43 at FDR < 5%), many of which were located in genes related to complications. Exploration studies into biological function reveal that these CpGs are enriched in binding sites for the C/EBP transcription factor, as well as enhancer/transcription regions in blood cells and haematopoietic stem cells, and open chromatin states in myeloid cells. Mediation analyses show that, remarkably, several CpGs in combination explain 68-97% of the association of mean-DCCT HbA1c with the risk of complications during EDIC. In summary, DNAme at key CpGs appears to mediate the association between hyperglycaemia and complications in metabolic memory, through modifying enhancer activity at myeloid and other cells.

DNA methylation age calculators reveal association with diabetic neuropathy in type 1 diabetes.

BACKGROUND: Many CpGs become hyper or hypo-methylated with age. Multiple methods have been developed by Horvath et al. to estimate DNA methylation (DNAm) age including Pan-tissue, Skin & Blood, PhenoAge, and GrimAge. Pan-tissue and Skin & Blood try to estimate chronological age in the normal population whereas PhenoAge and GrimAge use surrogate markers associated with mortality to estimate biological age and its departure from chronological age. Here, we applied Horvath's four methods to calculate and compare DNAm age in 499 subjects with type 1 diabetes (T1D) from the Diabetes Control and Complications Trial/Epidemiology of Diabetes Interventions and Complications (DCCT/EDIC) study using DNAm data measured by Illumina EPIC array in the whole blood. Association of the four DNAm ages with development of diabetic complications including cardiovascular diseases (CVD), nephropathy, retinopathy, and neuropathy, and their risk factors were investigated. RESULTS: Pan-tissue and GrimAge were higher whereas Skin & Blood and PhenoAge were lower than chronological age (p < 0.0001). DNAm age was not associated with the risk of CVD or retinopathy over 18-20 years after DNAm measurement. However, higher PhenoAge (ß = 0.023, p = 0.007) and GrimAge (ß = 0.029, p = 0.002) were associated with higher albumin excretion rate (AER), an indicator of diabetic renal disease, measured over time. GrimAge was also associated with development of both diabetic peripheral neuropathy (OR = 1.07, p = 9.24E-3) and cardiovascular autonomic neuropathy (OR = 1.06, p = 0.011). Both HbA1c (ß = 0.38, p = 0.026) and T1D duration (ß = 0.01, p = 0.043) were associated with higher PhenoAge. Employment (ß = - 1.99, p = 0.045) and leisure time (ß = - 0.81, p = 0.022) physical activity were associated with lower Pan-tissue and Skin & Blood, respectively. BMI (ß = 0.09, p = 0.048) and current smoking (ß = 7.13, p = 9.03E-50) were positively associated with Skin & Blood and GrimAge, respectively. Blood pressure, lipid levels, pulse rate, and alcohol consumption were not associated with DNAm age regardless of the method used. CONCLUSIONS: Various methods of measuring DNAm age are sub-optimal in detecting people at higher risk of developing diabetic complications although some work better than the others.

Association of IGF1 and VEGFA polymorphisms with diabetic retinopathy in Pakistani population.

AIMS: The incidence of microvascular complications, including diabetic retinopathy (DR), increases with duration of type 2 diabetes (T2D). Meta-GWAS have reported numerous single-nucleotide polymorphisms (SNPs) associated with T2D; however, no loci, achieving genome-wide significance has been reported for DR. Vascular endothelial growth factor A (VEGFA) and insulin-like growth factor 1 (IGF1) are considered as potential genetic candidates involved in T2D and DR progression. Moreover, the association of serum levels of these proteins with diabetes-related traits is controversial. Therefore, the current study was designed to evaluate the possible genetic predisposition and role of these circulating growth factors in serum in the pathophysiology of T2D and DR. METHODS: A cohort of 1126 individuals with T2D was collected including those without retinopathy (DNR = 573), non-progressive diabetic retinopathy (NPDR = 301) and progressive diabetic retinopathy (PDR = 252), and 348 healthy controls. Genomic DNA was isolated, and six SNPs: rs833061, rs13207351, rs1570360, rs2010963, rs5742632 and rs6214, were genotyped and results statistically analyzed. ELISA was performed on a subset of the samples to measure serum levels of IGF1 and VEGFA. RESULTS: The minor allele of rs6214 was associated with T2D [OR = 1.67 (95% CI 1.39-2.01, p = 4.9E-8)], rs13207351 was associated with NPDR [OR = 1.97 (95% CI 1.28-3.03, p = 9.0E-3)]when compared with DNR, and rs5742632 showed positive association with PDR [OR = 1.66 (95% CI 1.33-2.05, p = 1.0E-4)] compared to DNR. Lowered IGF1 serum levels were found to be associated with T2D, NPDR and PDR. CONCLUSIONS: IGF1 was found to increase the T2DM susceptibility as well as advanced DR, i.e., PDR, while VEGFA was found to be associated with early DR stage, i.e., NPDR.

gpart: human genome partitioning and visualization of high-density SNP data by identifying haplotype blocks.

SUMMARY: For the analysis of high-throughput genomic data produced by next-generation sequencing (NGS) technologies, researchers need to identify linkage disequilibrium (LD) structure in the genome. In this work, we developed an R package gpart which provides clustering algorithms to define LD blocks or analysis units consisting of SNPs. The visualization tool in gpart can display the LD structure and gene positions for up to 20 000 SNPs in one image. The gpart functions facilitate construction of LD blocks and SNP partitions for vast amounts of genome sequencing data within reasonable time and memory limits in personal computing environments. AVAILABILITY AND IMPLEMENTATION: The R package is available at https://bioconductor.org/packages/gpart. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

Meta-genome-wide association studies identify a locus on chromosome 1 and multiple variants in the MHC region for serum C-peptide in type 1 diabetes.

AIMS/HYPOTHESIS: The aim of this study was to identify genetic variants associated with beta cell function in type 1 diabetes, as measured by serum C-peptide levels, through meta-genome-wide association studies (meta-GWAS). METHODS: We performed a meta-GWAS to combine the results from five studies in type 1 diabetes with cross-sectionally measured stimulated, fasting or random C-peptide levels, including 3479 European participants. The p values across studies were combined, taking into account sample size and direction of effect. We also performed separate meta-GWAS for stimulated (n = 1303), fasting (n = 2019) and random (n = 1497) C-peptide levels. RESULTS: In the meta-GWAS for stimulated/fasting/random C-peptide levels, a SNP on chromosome 1, rs559047 (Chr1:238753916, T>A, minor allele frequency [MAF] 0.24-0.26), was associated with C-peptide (p = 4.13 × 10-8), meeting the genome-wide significance threshold (p < 5 × 10-8). In the same meta-GWAS, a locus in the MHC region (rs9260151) was close to the genome-wide significance threshold (Chr6:29911030, C>T, MAF 0.07-0.10, p = 8.43 × 10-8). In the stimulated C-peptide meta-GWAS, rs61211515 (Chr6:30100975, T/-, MAF 0.17-0.19) in the MHC region was associated with stimulated C-peptide (ß [SE] = - 0.39 [0.07], p = 9.72 × 10-8). rs61211515 was also associated with the rate of stimulated C-peptide decline over time in a subset of individuals (n = 258) with annual repeated measures for up to 6 years (p = 0.02). In the meta-GWAS of random C-peptide, another MHC region, SNP rs3135002 (Chr6:32668439, C>A, MAF 0.02-0.06), was associated with C-peptide (p = 3.49 × 10-8). Conditional analyses suggested that the three identified variants in the MHC region were independent of each other. rs9260151 and rs3135002 have been associated with type 1 diabetes, whereas rs559047 and rs61211515 have not been associated with a risk of developing type 1 diabetes. CONCLUSIONS/INTERPRETATION: We identified a locus on chromosome 1 and multiple variants in the MHC region, at least some of which were distinct from type 1 diabetes risk loci, that were associated with C-peptide, suggesting partly non-overlapping mechanisms for the development and progression of type 1 diabetes. These associations need to be validated in independent populations. Further investigations could provide insights into mechanisms of beta cell loss and opportunities to preserve beta cell function.

New Locus for Skin Intrinsic Fluorescence in Type 1 Diabetes Also Associated With Blood and Skin Glycated Proteins.

Skin fluorescence (SF) noninvasively measures advanced glycation end products (AGEs) in the skin and is a risk indicator for diabetes complications. N-acetyltransferase 2 (NAT2) is the only known locus influencing SF. We aimed to identify additional genetic loci influencing SF in type 1 diabetes (T1D) through a meta-analysis of genome-wide association studies (N = 1,359) including Diabetes Control and Complications Trial/Epidemiology of Diabetes Interventions and Complications (DCCT/EDIC) and Wisconsin Epidemiologic Study of Diabetic Retinopathy (WESDR). A locus on chromosome 1, rs7533564 (P = 1.9 × 10(-9)), was associated with skin intrinsic fluorescence measured by SCOUT DS (excitation 375 nm, emission 435-655 nm), which remained significant after adjustment for time-weighted HbA1c (P = 1.7 × 10(-8)). rs7533564 was associated with mean HbA1c in meta-analysis (P = 0.0225), mean glycated albumin (P = 0.0029), and glyoxal hydroimidazolones (P = 0.049), an AGE measured in skin biopsy collagen, in DCCT. rs7533564 was not associated with diabetes complications in DCCT/EDIC or with SF in subjects without diabetes (nondiabetic [ND]) (N = 8,721). In conclusion, we identified a new locus associated with SF in T1D subjects that did not show similar effect in ND subjects, suggesting a diabetes-specific effect. This association needs to be investigated in type 2 diabetes.