Whole-exome sequencing identifies novel protein-altering variants associated with serum apolipoprotein and lipid concentrations.

BACKGROUND: Dyslipidemia is a major risk factor for cardiovascular disease, and diabetes impacts the lipid metabolism through multiple pathways. In addition to the standard lipid measurements, apolipoprotein concentrations provide added awareness of the burden of circulating lipoproteins. While common genetic variants modestly affect the serum lipid concentrations, rare genetic mutations can cause monogenic forms of hypercholesterolemia and other genetic disorders of lipid metabolism. We aimed to identify low-frequency protein-altering variants (PAVs) affecting lipoprotein and lipid traits. METHODS: We analyzed whole-exome (WES) and whole-genome sequencing (WGS) data of 481 and 474 individuals with type 1 diabetes, respectively. The phenotypic data consisted of 79 serum lipid and apolipoprotein phenotypes obtained with clinical laboratory measurements and nuclear magnetic resonance spectroscopy. RESULTS: The single-variant analysis identified an association between the LIPC p.Thr405Met (rs113298164) and serum apolipoprotein A1 concentrations (p=7.8×10-8). The burden of PAVs was significantly associated with lipid phenotypes in LIPC, RBM47, TRMT5, GTF3C5, MARCHF10, and RYR3 (p<2.9×10-6). The RBM47 gene is required for apolipoprotein B post-translational modifications, and in our data, the association between RBM47 and apolipoprotein C-III concentrations was due to a rare 21 base pair p.Ala496-Ala502 deletion; in replication, the burden of rare deleterious variants in RBM47 was associated with lower triglyceride concentrations in WES of >170,000 individuals from multiple ancestries (p=0.0013). Two PAVs in GTF3C5 were highly enriched in the Finnish population and associated with cardiovascular phenotypes in the general population. In the previously known APOB gene, we identified novel associations at two protein-truncating variants resulting in lower serum non-HDL cholesterol (p=4.8×10-4), apolipoprotein B (p=5.6×10-4), and LDL cholesterol (p=9.5×10-4) concentrations. CONCLUSIONS: We identified lipid and apolipoprotein-associated variants in the previously known LIPC and APOB genes, as well as PAVs in GTF3C5 associated with LDLC, and in RBM47 associated with apolipoprotein C-III concentrations, implicated as an independent CVD risk factor. Identification of rare loss-of-function variants has previously revealed genes that can be targeted to prevent CVD, such as the LDL cholesterol-lowering loss-of-function variants in the PCSK9 gene. Thus, this study suggests novel putative therapeutic targets for the prevention of CVD.

Association of Coding Variants in Hydroxysteroid 17-beta Dehydrogenase 14 (HSD17B14) with Reduced Progression to End Stage Kidney Disease in Type 1 Diabetes.

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 ESKD in individuals with type 1 diabetes at advanced kidney disease stage. METHODS: 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 (n=2372/1115 events all cohorts) and replicating in two retrospective case-control studies (n=1072 cases, 752 controls). Deep resequencing of the top associated gene in five 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-ß dehydrogenase 14 gene (HSD17B14), predicted to affect protein structure, had a net protective effect against development of ESKD at exome-wide significance (n=4196; P value=3.3 × 10-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 CKD-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.

Differential metabolomic signatures of declining renal function in Types 1 and 2 diabetes.

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.

Liver nucleotide biosynthesis is linked to protection from vascular complications in individuals with long-term type 1 diabetes.

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.

Genome-Wide Association Study of Diabetic Kidney Disease Highlights Biology Involved in Glomerular Basement Membrane Collagen.

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.

Metabolomic Profile Predicts Development of Microalbuminuria in Individuals with Type 1 Diabetes.

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.

Confirmation of GLRA3 as a susceptibility locus for albuminuria in Finnish patients with type 1 diabetes.

Urinary albumin excretion is an early sign of diabetic kidney disease, affecting every third individual with diabetes. Despite substantial estimated heritability, only variants in the GLRA3 gene have been genome-wide significantly associated (p-value < 5 × 10-8) with diabetic albuminuria, in Finnish individuals with type 1 diabetes; However, replication attempt in non-Finnish Europeans with type 1 diabetes showed nominally significant association in the opposite direction, suggesting a population-specific effect, but simultaneously leaving the finding controversial. In this study, the association between the common rs10011025 variant in the GLRA3 locus, and albuminuria, was confirmed in 1259 independent Finnish individuals with type 1 diabetes (p = 0.0013), and meta-analysis of all Finnish individuals yielded a genome-wide significant association. The association was particularly pronounced in subjects not reaching the treatment target for blood glucose levels (HbA1c > 7%; N = 2560, p = 1.7 × 10-9). Even though further studies are needed to pinpoint the causal variants, dissecting the association at the GLRA3 locus may uncover novel molecular mechanisms for diabetic albuminuria irrespective of population background.

The Glanville fritillary genome retains an ancient karyotype and reveals selective chromosomal fusions in Lepidoptera.

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.