Search

AimsIndividuals with type 1 diabetes (T1D) carry a markedly increased risk of stroke, with distinct clinical and neuroimaging characteristics as compared to those without diabetes. Using whole-genome sequencing (WGS) and whole-exome sequencing (WES), we aimed to find rare and low-frequency genomic variants associated with stroke in T1D. The lead findings were followed up in various datasets to replicate the findings and to assess their specificity to diabetes.Methods and ResultsWe studied stroke genetics in 1,051 individuals with T1D using WGS or WES. We analysed the genome with single-variant analyses, gene aggregate analyses, and aggregate analyses on genomic windows, enhancers and promoters. Furthermore, we attempted replication in T1D using a genome-wide association study (N=3,945) and direct genotyping (N=3,600), and in the general population from the FinnGen project and UK Biobank summary statistics. We identified a rare missense mutation onSREBF1associated with hemorrhagic stroke (rs114001633, p.Pro227Leu,p-value=8.96×10-9), which further replicated in T1D. Using gene aggregate analysis with protein altering or protein truncating variants, we identified exome-wide significant genes:ANK1andLRRN1displayed replication evidence in T1D, whileLRRN1,HAS1andUACAreplicated in the general population (UK Biobank). Furthermore, we performed sliding-window analyses and identified 14 genome-wide significant windows for stroke on 4q33-34.1, of which two replicated in T1D, and a suggestive genomic window onLINC01500, which replicated in T1D. Finally, with the regulatory region aggregate analyses, we identified a stroke-associatedTRPM2-ASpromoter (p-value=5.78×10-6), which we validated with an in vitro cell-based assay.TRPM2has been previously linked to ischemic stroke.ConclusionsHere, we report the first genome-wide analysis on stroke in individuals with diabetes. We identified multiple stroke risk loci with evidence of replication: 4q33-34.1,SREBF1, andANK1for stroke in T1D; andHAS1,UACA,LRRN1,LINC01500, andTRPM2-ASpromoter for stroke potentially generalizable to the non-diabetic population.

Background Type 1 diabetes (T1D) is an autoimmune disease affecting individuals in the early years of life. Although previous studies have identified genetic loci influencing T1D diagnosis age, these studies did not investigate the genome with high resolution. Objective and methods We performed a genome-wide meta-analysis for age at diagnosis with cohorts from Finland (Finnish Diabetic Nephropathy Study), the United Kingdom (UK Genetic Resource Investigating Diabetes) and Sardinia. Through SNP associations, transcriptome-wide association analysis linked T1D diagnosis age and gene expression. Results We identified two chromosomal regions associated with T1D diagnosis age: multiple independent variants in the HLA region on chromosome 6 and a locus on chromosome 17q12. We performed gene-level association tests with transcriptome prediction models from two whole blood datasets, lymphocyte cell line, spleen, pancreas and small intestine tissues. Of the non-HLA genes, lower PNMT expression in whole blood, and higher IKZF3 and ZPBP2, and lower ORMDL3 and GSDMB transcription levels in multiple tissues were associated with lower T1D diagnosis age (FDR = 0.05). These genes lie on chr17q12 which is associated with T1D, other autoimmune diseases, and childhood asthma. Additionally, higher expression of PHF20L1, a gene not previously implicated in T1D, was associated with lower diagnosis age in lymphocytes, pancreas, and spleen. Altogether, the non-HLA associations were enriched in open chromatin in various blood cells, blood vessel tissues and foetal thymus tissue. Conclusion Multiple genes on chr17q12 and PHF20L1 on chr8 were associated with T1D diagnosis age and only further studies may elucidate the role of these genes for immunity and T1D onset.

Background Rare variants in gene coding regions likely have a greater impact on disease-related phenotypes than common variants through disruption of their encoded protein. We searched for rare variants associated with onset of end stage kidney disease (ESKD) in individuals with type 1 diabetes at advanced kidney disease stage. Methods Gene-based exome array analysis of 15,449 genes in 5 large incidence cohorts of individuals with type 1 diabetes and proteinuria were analyzed for survival time-to-ESKD, testing the top gene in a 6th cohort (N=2,372/1,115 events all cohorts) and replicating in two retrospective case-control studies (N=1,072 cases, 752 controls). Deep resequencing of the top associated gene in 5 cohorts confirmed the findings. We performed immunohistochemistry and gene expression experiments in human control and diseased cells, and in mouse ischemia reperfusion and aristolochic acid nephropathy models. Results Protein coding variants in the hydroxysteroid 17-beta dehydrogenase 14 gene (HSD17B14), predicted to affect protein structure, had a net protective effect against development of ESKD at exome-wide significance (N=4,196; p-value=3.3x10-7). The HSD17B14 gene and encoded enzyme were robustly expressed in healthy human kidney, maximally in proximal tubular cells. Paradoxically, gene and protein expression were attenuated in human diabetic proximal tubules and in mouse kidney injury models. Expressed HSD17B14 gene and protein levels remained low without recovery after 21 days in a murine ischemic reperfusion injury model. Decreased gene expression was found in other chronic kidney disease-associated renal pathologies. Conclusions HSD17B14 gene is mechanistically involved in diabetic kidney disease. The encoded sex steroid enzyme is a druggable target, potentially opening a new avenue for therapeutic development.

Dyslipidemia is a major risk factor for cardiovascular disease. While common genetic variants are known to modestly affect the serum lipid concentrations, rare genetic mutations can cause monogenic forms of hypercholesteremia and other genetic disorders of lipid metabolism. Aiming to identify low-frequency protein-altering variants (PAVs) affecting lipoprotein and lipid traits, we analyzed whole-exome and whole-genome sequencing data of 481 and 573 individuals with type 1 diabetes, respectively. The phenotypic data consisted of 97 serum lipid, apolipoprotein, or other metabolic phenotypes obtained with clinical laboratory measurements and nuclear magnetic resonance (NMR) technology. Single variant analysis identified a novel association between LIPC p.Thr405Met (rs113298164) and serum apolipoprotein-A1 levels (p=7.8×10−8). In the APOB gene, we identified novel associations at two protein-truncating variants (PTVs) resulting in lower serum apolipoprotein B levels (p=5.6×10−4). The burden of PAVs was significantly associated with lipid phenotypes in LIPC, RBM47, TRMT5, and GTF3C5 (p−6). The RBM47 gene is required for apolipoprotein-B post-translational modifications, and in our data, the association between RBM47 and apolipoprotein C-III levels was led by a rare 21 base pair Ala496-Ala502 deletion; as replication, the burden of rare deleterious variants in RBM47 was associated with TG-to-HDLC ratio in WES of 20,917 individuals (p=0.0093). Two PAVs in GTF3C5 were highly Finnish-enriched and associated with cardiovascular phenotypes in external data, whereby the TRMT5 p.Ser185Cys lead variant was associated with stroke phenotypes. Altogether, we identified both novel variant associations in known lipid genes, as well as novel genes implicated in lipoprotein metabolism.

Background Chronic kidney disease (CKD) shows different clinical features in Types1 (T1D) and 2 diabetes (T2D). Metabolomics have recently provided useful contribution to the identification of biomarkers of CKD progression in either form of the disease. However, no studies have so far compared plasma metabolomics between T1D and T2D in order to identify differential signatures of progression of estimated glomerular filtration rate (eGFR) decline. Methods We used two large cohorts of T1D (from Finland) and T2D (from Italy) patients followed up to 7 and 3 years, respectively. In both groups, progression was defined as the top quartile of yearly decline in eGFR. Pooled data from the two groups were analysed by univariate and bivariate random forest (RF), and confirmed by bivariate partial least squares (PLS) analysis, the response variables being type of diabetes and eGFR progression. Results In progressors, yearly eGFR loss was significantly larger in T2D [−5.3 (3.0), median (interquartile range)mL/min/1.73 m2/year] than T1D [−3.7 (3.1) mL/min/1.73 m2/year ; P = 0.018]. Out of several hundreds, bivariate RF extracted 22 metabolites associated with diabetes type (all higher in T1D than T2D except for 5-methylthioadenosine, pyruvate and β-hydroxypyruvate) and 13 molecules associated with eGFR progression (all higher in progressors than non-progressors except for sphyngomyelin). Three of the selected metabolites (histidylphenylalanine, leucylphenylalanine, tryptophylasparagine) showed a significant interaction between disease type and progression. Only eight metabolites were common to both bivariate RF and PLS. Conclusions Identification of metabolomic signatures of CKD progression is partially dependent on the statistical model. Dual analysis identified molecules specifically associated with progressive renal impairment in both T1D and T2D.

Rare variants in gene coding regions likely have a greater impact on disease-related phenotypes than common variants through disruption of their encoded protein. We searched for rare variants associated with onset of ESKD in individuals with type 1 diabetes at advanced kidney disease stage.Gene-based exome array analyses of 15,449 genes in five large incidence cohorts of individuals with type 1 diabetes and proteinuria were analyzed for survival time to ESKD, testing the top gene in a sixth cohort (Protein coding variants in the hydroxysteroid 17

Identification of biomarkers associated with protection from developing diabetic complications is a prerequisite for an effective prevention and treatment. The aim of the present study was to identify clinical and plasma metabolite markers associated with freedom from vascular complications in people with very long duration of type 1 diabetes (T1D). Individuals with T1D, who despite having longer than 30 years of diabetes duration never developed major macro- or microvascular complications (non-progressors; NP) were compared with those who developed vascular complications within 25 years from diabetes onset (rapid progressors; RP) in the Scandinavian PROLONG (n = 385) and DIALONG (n = 71) cohorts. The DIALONG study also included 75 healthy controls. Plasma metabolites were measured using gas and/or liquid chromatography coupled to mass spectrometry. Lower hepatic fatty liver indices were significant common feature characterized NPs in both studies. Higher insulin sensitivity and residual ß-cell function (C-peptide) were also associated with NPs in PROLONG. Protection from diabetic complications was associated with lower levels of the glycolytic metabolite pyruvate and APOCIII in PROLONG, and with lower levels of thiamine monophosphate and erythritol, a cofactor and intermediate product in the pentose phosphate pathway as well as higher phenylalanine, glycine and serine in DIALONG. Furthermore, T1D individuals showed elevated levels of picolinic acid as compared to the healthy individuals. The present findings suggest a potential beneficial shunting of glycolytic substrates towards the pentose phosphate and one carbon metabolism pathways to promote nucleotide biosynthesis in the liver. These processes might be linked to higher insulin sensitivity and lower liver fat content, and might represent a mechanism for protection from vascular complications in individuals with long-term T1D.

BACKGROUND: Although diabetic kidney disease demonstrates both familial clustering and single nucleotide polymorphism heritability, the specific genetic factors influencing risk remain largely unknown.METHODS: To identify genetic variants predisposing to diabetic kidney disease, we performed genome-wide association study (GWAS) analyses. Through collaboration with the Diabetes Nephropathy Collaborative Research Initiative, we assembled a large collection of type 1 diabetes cohorts with harmonized diabetic kidney disease phenotypes. We used a spectrum of ten diabetic kidney disease definitions based on albuminuria and renal function.RESULTS: Our GWAS meta-analysis included association results for up to 19,406 individuals of European descent with type 1 diabetes. We identified 16 genome-wide significant risk loci. The variant with the strongest association (rs55703767) is a common missense mutation in the collagen type IV alpha 3 chain (COL4A3) gene, which encodes a major structural component of the glomerular basement membrane (GBM). Mutations in COL4A3 are implicated in heritable nephropathies, including the progressive inherited nephropathy Alport syndrome. The rs55703767 minor allele (Asp326Tyr) is protective against several definitions of diabetic kidney disease, including albuminuria and ESKD, and demonstrated a significant association with GBM width; protective allele carriers had thinner GBM before any signs of kidney disease, and its effect was dependent on glycemia. Three other loci are in or near genes with known or suggestive involvement in this condition (BMP7) or renal biology (COLEC11 and DDR1).CONCLUSIONS: The 16 diabetic kidney disease-associated loci may provide novel insights into the pathogenesis of this condition and help identify potential biologic targets for prevention and treatment.

Diabetic kidney disease (DKD) is a heritable but poorly understood complication of diabetes. To identify genetic variants predisposing to DKD, we performed genome-wide association analyses in 19,406 individuals with type 1 diabetes (T1D) using a spectrum of DKD definitions basedon albuminuria and renal function. We identified 16 genome-wide significant loci. The variant with the strongest association (rs55703767) is a common missense mutation in the collagen type IV alpha 3 chain(COL4A3)gene, which encodes a major structural component of the glomerular basement membrane (GBM) implicated in heritable nephropathies. The rs55703767 minor allele (Asp326Tyr) is protective against several definitions of DKD, including albuminuria and end-stage renal disease. Three other loci are in or near genes with known or suggestive involvement in DKD(BMP7)or renal biology (COLEC11andDDR1). The 16 DKD-associated loci provide novel insights into the pathogenesis of DKD, identifying potential biological targets for prevention and treatment.

Elevated urinary albumin excretion (microalbuminuria) is an early marker of diabetic nephropathy, but there is an unmet need for better biomarkers that capture the individuals at risk with higher accuracy and earlier than the current markers do. We performed an untargeted metabolomic study to assess baseline differences between individuals with type 1 diabetes who either developed microalbuminuria or remained normoalbuminuric. A total of 102 individuals progressed to microalbuminuria during a median follow-up of 3.2 years, whereas 98 sex-, age- and body mass index (BMI) matched non-progressors remained normoalbuminuric during a median follow-up of 7.1 years. Metabolomic screening identified 1,242 metabolites, out of which 111 differed significantly between progressors and non-progressors after adjustment for age of diabetes onset, baseline glycosylated hemoglobin A1c (HbA1c), and albumin excretion rate (AER). The metabolites that predicted development of microalbumiuria included several uremic toxins and carnitine metabolism related molecules. Iterative variable selection indicated erythritol, 3-phenylpropionate, and N-trimethyl-5-aminovalerate as the best set of variables to predict development of microalbuminuria. A metabolomic index based on these metabolites improved the prediction of incident microalbuminuria on top of the clinical variables age of diabetes onset, baseline HbA1c and AER (ROCAUC = 0.842 vs 0.797), highlighting their ability to predict early-phase diabetic nephropathy.

Made available in DSpace on 2015-10-21T20:30:40Z (GMT). No. of bitstreams: 0 Previous issue date: 2015-04-01. Added 1 bitstream(s) on 2015-10-22T09:48:02Z : No. of bitstreams: 1 WOS000355148800002.pdf: 2850380 bytes, checksum: 09d6a629a3febf17eadca613619ec995 (MD5) National Institutes of Health Marie Curie International Incoming Fellowship from the European Union Fundação de Amparo à Pesquisa do Estado de São Paulo (FAPESP) University of Tampere Finland Distinguished Professor Programme of the Academy of Finland The large number of complete mitochondrial DNA (mtDNA) sequences available for metazoan species makes it a good system for studying genome diversity, although little is known about the mechanisms that promote and/or are correlated with the evolution of this organellar genome. By investigating the molecular evolutionary history of the catalytic and accessory subunits of the mtDNA polymerase, pol gamma, we sought to develop mechanistic insight into its function that might impact genome structure by exploring the relationships between DNA replication and animal mitochondrial genome diversity. We identified three evolutionary patterns among metazoan pol gamma s. First, a trend toward stabilization of both sequence and structure occurred in vertebrates, with both subunits evolving distinctly from those of other animal groups, and acquiring at least four novel structural elements, the most important of which is the HLH-3 beta (helix-loop-helix, 3 beta-sheets) domain that allows the accessory subunit to homodimerize. Second, both subunits of arthropods and tunicates have become shorter and evolved approximately twice as rapidly as their vertebrate homologs. And third, nematodes have lost the gene for the accessory subunit, which was accompanied by the loss of its interacting domain in the catalytic subunit of pol gamma, and they show the highest rate of molecular evolution among all animal taxa. These findings correlate well with the mtDNA genomic features of each group described above, and with their modes of DNA replication, although a substantive amount of biochemical work is needed to draw conclusive links regarding the latter. Describing the parallels between evolution of pol gamma and metazoan mtDNA architecture may also help in understanding the processes that lead to mitochondrial dysfunction and to human disease-related phenotypes. Institute of Biosciences and Medical Technology, University of Tampere, Finland Department of Biochemistry and Molecular Biology and Center for Mitochondrial Science and Medicine, Michigan State University Departamento de Tecnologia, Faculdade de Ciências Agrárias e Veterinárias, Universidade Estadual Paulista “Júlio de Mesquita Filho,” Jaboticabal, SP, Brazil National Institutes of Health: 45295 Marie Curie International Incoming Fellowship from the European Union: GA328988 FAPESP: 2014/02253-6

Previous studies have reported that chromosome synteny in Lepidoptera has been well conserved, yet the number of haploid chromosomes varies widely from 5 to 223. Here we report the genome (393 Mb) of the Glanville fritillary butterfly (Melitaea cinxia; Nymphalidae), a widely recognized model species in metapopulation biology and eco-evolutionary research, which has the putative ancestral karyotype of n=31. Using a phylogenetic analyses of Nymphalidae and of other Lepidoptera, combined with orthologue-level comparisons of chromosomes, we conclude that the ancestral lepidopteran karyotype has been n=31 for at least 140 My. We show that fusion chromosomes have retained the ancestral chromosome segments and very few rearrangements have occurred across the fusion sites. The same, shortest ancestral chromosomes have independently participated in fusion events in species with smaller karyotypes. The short chromosomes have higher rearrangement rate than long ones. These characteristics highlight distinctive features of the evolutionary dynamics of butterflies and moths., Marie Curie International Fellowship (PIOF-GA-2011-303312)