Genome-wide association study identifies host genetic variants influencing oral microbiota diversity and metabolic health.

The microbial communities of the oral cavity are important elements of oral and systemic health. With emerging evidence highlighting the heritability of oral bacterial microbiota, this study aimed to identify host genome variants that influence oral microbial traits. Using data from 16S rRNA gene amplicon sequencing, we performed genome-wide association studies with univariate and multivariate traits of the salivary microbiota from 610 unrelated adults from the Danish ADDITION-PRO cohort. We identified six single nucleotide polymorphisms (SNPs) in human genomes that showed associations with abundance of bacterial taxa at different taxonomical tiers (P < 5 × 10-8). Notably, SNP rs17793860 surpassed our study-wide significance threshold (P < 1.19 × 10-9). Additionally, rs4530093 was linked to bacterial beta diversity (P < 5 × 10-8). Out of these seven SNPs identified, six exerted effects on metabolic traits, including glycated hemoglobin A1c, triglyceride and high-density lipoprotein cholesterol levels, the risk of type 2 diabetes and stroke. Our findings highlight the impact of specific host SNPs on the composition and diversity of the oral bacterial community. Importantly, our results indicate an intricate interplay between host genetics, the oral microbiota, and metabolic health. We emphasize the need for integrative approaches considering genetic, microbial, and metabolic factors.

Circulating metabolomic markers in association with overall burden of microvascular complications in type 1 diabetes.

INTRODUCTION: Diabetic retinopathy (DR), diabetic kidney disease (DKD) and distal symmetric polyneuropathy (DSPN) share common pathophysiology and pose an additive risk of early mortality. RESEARCH DESIGN AND METHODS: In adults with type 1 diabetes, 49 metabolites previously associated with either DR or DKD were assessed in relation to presence of DSPN. Metabolites overlapping in significance with presence of all three complications were assessed in relation to microvascular burden severity (additive number of complications-ie, presence of DKD±DR±DSPN) using linear regression models. Subsequently, the same metabolites were assessed with progression to endpoints: soft microvascular events (progression in albuminuria grade, ≥30% estimated glomerular filtration rate (eGFR) decline, or any progression in DR grade), hard microvascular events (progression to proliferative DR, chronic kidney failure, or ≥40% eGFR decline), and hard microvascular or macrovascular events (hard microvascular events, cardiovascular events (myocardial infarction, stroke, or arterial interventions), or cardiovascular mortality), using Cox models. All models were adjusted for sex, baseline age, diabetes duration, systolic blood pressure, HbA1c, body mass index, total cholesterol, smoking, and statin treatment. RESULTS: The full cohort investigated consisted of 487 participants. Mean (SD) follow-up was 4.8 (2.9, 5.7) years. Baseline biothesiometry was available in 202 participants, comprising the cross-sectional cohort. Eight metabolites were significantly associated with presence of DR, DKD, and DSPN, and six with additive microvascular burden severity. In the full cohort longitudinal analysis, higher levels of 3,4-dihydroxybutanoic acid (DHBA), 2,4-DHBA, ribonic acid, glycine, and ribitol were associated with development of events in both crude and adjusted models. Adding 3,4-DHBA, ribonic acid, and glycine to a traditional risk factor model improved the discrimination of hard microvascular events. CONCLUSIONS: While prospective studies directly assessing the predictive ability of these markers are needed, our results strengthen the role of clinical metabolomics in relation to risk assessment of diabetic complications in chronic type 1 diabetes.

Randomized Trial of SGLT2 Inhibitor Identifies Target Proteins in Diabetic Kidney Disease.

Introduction: Sodium-glucose cotransporter 2 inhibitors (SGLT2i) have emerged as novel therapeutics to treat diabetic kidney disease (DKD). Although the beneficial effects of SGLT2i have been demonstrated, their target mechanisms on kidney function are unknown. The current study aimed to elucidate these mechanisms by studying SGLT2i-induced changes in the urinary proteome of persons with type 2 diabetes (T2D) and DKD. Methods: A total of 40 participants with T2D were enrolled in a double-blinded randomized cross-over trial at the Steno Diabetes Center Copenhagen, Denmark. They were treated with 10 mg of dapagliflozin for 12 weeks. Thirty-two participants with complete urinary proteomics measures before and after the trial were included. All participants received renin-angiotensin system blockade and had albuminuria, (urine albumin-to-creatinine ratio [UACR] ≥30 mg/g). A type 1 diabetes (T1D) cohort consisting of healthy controls and persons with DKD was included for validation. Urinary proteome changes were analyzed using Wilcoxon signed-rank test. Functional enrichment analysis was conducted to discover affected biological processes. Results: Dapagliflozin treatment significantly (Padjusted < 0.05) affected 36 urinary peptide fragments derived from 19 proteins. Eighteen proteins were correspondingly reflected in the validation cohort. A multifold change in peptide abundance was observed in many proteins (A1BG, urinary albumin [ALB], Caldesmon 1, COLCRNN, heat shock protein 90-ß [HSP90AB1], IGLL5, peptidase inhibitor 16 [PI16], prostaglandin-H2-D-isomerase [PTGDS], SERPINA1). These also included urinary biomarkers of kidney fibrosis and function (type I and III collagens and albumin). Biological processes relating to inflammation, wound healing, and kidney fibrosis were enriched. Conclusion: The current study discovers the urinary proteome impacted by the SGLT2i, thereby providing new potential target sites and pathways, especially relating to wound healing and inflammation.

Genome-Wide Analyses of Vocabulary Size in Infancy and Toddlerhood: Associations With Attention-Deficit/Hyperactivity Disorder, Literacy, and Cognition-Related Traits.

BACKGROUND: The number of words children produce (expressive vocabulary) and understand (receptive vocabulary) changes rapidly during early development, partially due to genetic factors. Here, we performed a meta-genome-wide association study of vocabulary acquisition and investigated polygenic overlap with literacy, cognition, developmental phenotypes, and neurodevelopmental conditions, including attention-deficit/hyperactivity disorder (ADHD). METHODS: We studied 37,913 parent-reported vocabulary size measures (English, Dutch, Danish) for 17,298 children of European descent. Meta-analyses were performed for early-phase expressive (infancy, 15-18 months), late-phase expressive (toddlerhood, 24-38 months), and late-phase receptive (toddlerhood, 24-38 months) vocabulary. Subsequently, we estimated single nucleotide polymorphism-based heritability (SNP-h2) and genetic correlations (rg) and modeled underlying factor structures with multivariate models. RESULTS: Early-life vocabulary size was modestly heritable (SNP-h2 = 0.08-0.24). Genetic overlap between infant expressive and toddler receptive vocabulary was negligible (rg = 0.07), although each measure was moderately related to toddler expressive vocabulary (rg = 0.69 and rg = 0.67, respectively), suggesting a multifactorial genetic architecture. Both infant and toddler expressive vocabulary were genetically linked to literacy (e.g., spelling: rg = 0.58 and rg = 0.79, respectively), underlining genetic similarity. However, a genetic association of early-life vocabulary with educational attainment and intelligence emerged only during toddlerhood (e.g., receptive vocabulary and intelligence: rg = 0.36). Increased ADHD risk was genetically associated with larger infant expressive vocabulary (rg = 0.23). Multivariate genetic models in the ALSPAC (Avon Longitudinal Study of Parents and Children) cohort confirmed this finding for ADHD symptoms (e.g., at age 13; rg = 0.54) but showed that the association effect reversed for toddler receptive vocabulary (rg = -0.74), highlighting developmental heterogeneity. CONCLUSIONS: The genetic architecture of early-life vocabulary changes during development, shaping polygenic association patterns with later-life ADHD, literacy, and cognition-related traits.

Bone mineral density loci specific to the skull portray potential pleiotropic effects on craniosynostosis.

Skull bone mineral density (SK-BMD) provides a suitable trait for the discovery of key genes in bone biology, particularly to intramembranous ossification, not captured at other skeletal sites. We perform a genome-wide association meta-analysis (n ~ 43,800) of SK-BMD, identifying 59 loci, collectively explaining 12.5% of the trait variance. Association signals cluster within gene-sets involved in skeletal development and osteoporosis. Among the four novel loci (ZIC1, PRKAR1A, AZIN1/ATP6V1C1, GLRX3), there are factors implicated in intramembranous ossification and as we show, inherent to craniosynostosis processes. Functional follow-up in zebrafish confirms the importance of ZIC1 on cranial suture patterning. Likewise, we observe abnormal cranial bone initiation that culminates in ectopic sutures and reduced BMD in mosaic atp6v1c1 knockouts. Mosaic prkar1a knockouts present asymmetric bone growth and, conversely, elevated BMD. In light of this evidence linking SK-BMD loci to craniofacial abnormalities, our study provides new insight into the pathophysiology, diagnosis and treatment of skeletal diseases.

Circulating Metabolomic and Lipidomic Signatures Identify a Type 2 Diabetes Risk Profile in Low-Birth-Weight Men with Non-Alcoholic Fatty Liver Disease.

The extent to which increased liver fat content influences differences in circulating metabolites and/or lipids between low-birth-weight (LBW) individuals, at increased risk of type 2 diabetes (T2D), and normal-birth-weight (NBW) controls is unknown. The objective of the study was to perform untargeted serum metabolomics and lipidomics analyses in 26 healthy, non-obese early-middle-aged LBW men, including five men with screen-detected and previously unrecognized non-alcoholic fatty liver disease (NAFLD), compared with 22 age- and BMI-matched NBW men (controls). While four metabolites (out of 65) and fifteen lipids (out of 279) differentiated the 26 LBW men from the 22 NBW controls (p ≤ 0.05), subgroup analyses of the LBW men with and without NAFLD revealed more pronounced differences, with 11 metabolites and 56 lipids differentiating (p ≤ 0.05) the groups. The differences in the LBW men with NAFLD included increased levels of ornithine and tyrosine (PFDR ≤ 0.1), as well as of triglycerides and phosphatidylcholines with shorter carbon-chain lengths and fewer double bonds. Pathway and network analyses demonstrated downregulation of transfer RNA (tRNA) charging, altered urea cycling, insulin resistance, and an increased risk of T2D in the LBW men with NAFLD. Our findings highlight the importance of increased liver fat in the pathogenesis of T2D in LBW individuals.

Dissecting the causal association between inflammation and post-traumatic stress disorder: A bidirectional Mendelian randomization study.

BACKGROUND: Accumulating evidence showed a bidirectional association between post-traumatic stress disorder and inflammation. However, whether the association is causal remains unclear. We aimed to evaluate the causal relationships between inflammatory cytokines and post-traumatic stress disorder using two-sample bi-directional Mendelian randomization analysis. METHODS: Single nucleotide polymorphism from genome-wide association studies of inflammatory cytokines, C-reactive protein, and post-traumatic stress disorder (23,212 patients and 151,447 controls) was selected as instrumental variables. The causal associations were estimated by inverse variance weighting with sensitivity analyses using weighted median, MR-Egger, and MR-PRESSO methods. RESULTS: We observed suggestive associations of genetically predicted interleukin-17 (IL-17) and RANTES with post-traumatic stress disorder. One standard deviation (SD) increase in genetically predicted level of IL-17 lowered the risk of post-traumatic stress disorder with an odds ratio (OR) of 0.902 (95 % CI = 0.828, 0.984, P = 0.02). One SD higher genetically predicted RANTES (CCL5) concentration increased post-traumatic stress disorder risk (OR = 1.067, 95 % CI = 1.005, 1.133, P = 0.032). However, we found no evidence of causal associations of post-traumatic stress disorder with the selected inflammatory cytokines and biomarkers. We observed no evidence supporting the presence of pleiotropy. The results of sensitivity analyses demonstrated the same directions and similar effect sizes as the primary findings. LIMITATIONS: Potential pleiotropy, possible weak instruments, and low statistical power limited our findings. CONCLUSION: Inflammation was suggestively causally associated with the risk of post-traumatic stress disorder, and inflammatory cytokines had no downstream effect on post-traumatic stress disorder. Further studies are needed to explain the mechanisms of systemic inflammation and neuroinflammation in post-traumatic stress disorder.

Four missense genetic variants in CUBN are associated with higher levels of eGFR in non-diabetes but not in diabetes mellitus or its subtypes: A genetic association study in Europeans.

Aim: Rare genetic variants in the CUBN gene encoding the main albumin-transporter in the proximal tubule of the kidneys have previously been associated with microalbuminuria and higher urine albumin levels, also in diabetes. Sequencing studies in isolated proteinuria suggest that these variants might not affect kidney function, despite proteinuria. However, the relation of these CUBN missense variants to the estimated glomerular filtration rate (eGFR) is largely unexplored. We hereby broadly examine the associations between four CUBN missense variants and eGFRcreatinine in Europeans with Type 1 (T1D) and Type 2 Diabetes (T2D). Furthermore, we sought to deepen our understanding of these variants in a range of single- and aggregate- variant analyses of other kidney-related traits in individuals with and without diabetes mellitus. Methods: We carried out a genetic association-based linear regression analysis between four CUBN missense variants (rs141640975, rs144360241, rs45551835, rs1801239) and eGFRcreatinine (ml/min/1.73 m2, CKD-EPIcreatinine(2012), natural log-transformed) in populations with T1D (n ~ 3,588) or T2D (n ~ 31,155) from multiple European studies and in individuals without diabetes from UK Biobank (UKBB, n ~ 370,061) with replication in deCODE (n = 127,090). Summary results of the diabetes-group were meta-analyzed using the fixed-effect inverse-variance method. Results: Albeit we did not observe associations between eGFRcreatinine and CUBN in the diabetes-group, we found significant positive associations between the minor alleles of all four variants and eGFRcreatinine in the UKBB individuals without diabetes with rs141640975 being the strongest (Effect=0.02, PeGFR_creatinine=2.2 × 10-9). We replicated the findings for rs141640975 in the Icelandic non-diabetes population (Effect=0.026, PeGFR_creatinine=7.7 × 10-4). For rs141640975, the eGFRcreatinine-association showed significant interaction with albuminuria levels (normo-, micro-, and macroalbuminuria; p = 0.03). An aggregated genetic risk score (GRS) was associated with higher urine albumin levels and eGFRcreatinine. The rs141640975 variant was also associated with higher levels of eGFRcreatinine-cystatin C (ml/min/1.73 m2, CKD-EPI2021, natural log-transformed) and lower circulating cystatin C levels. Conclusions: The positive associations between the four CUBN missense variants and eGFR in a large population without diabetes suggests a pleiotropic role of CUBN as a novel eGFR-locus in addition to it being a known albuminuria-locus. Additional associations with diverse renal function measures (lower cystatin C and higher eGFRcreatinine-cystatin C levels) and a CUBN-focused GRS further suggests an important role of CUBN in the future personalization of chronic kidney disease management in people without diabetes.

Multiomics signatures of type 1 diabetes with and without albuminuria.

Aims/hypothesis: To identify novel pathophysiological signatures of longstanding type 1 diabetes (T1D) with and without albuminuria we investigated the gut microbiome and blood metabolome in individuals with T1D and healthy controls (HC). We also mapped the functional underpinnings of the microbiome in relation to its metabolic role. Methods: One hundred and sixty-one individuals with T1D and 50 HC were recruited at the Steno Diabetes Center Copenhagen, Denmark. T1D cases were stratified based on levels of albuminuria into normoalbuminuria, moderate and severely increased albuminuria. Shotgun sequencing of bacterial and viral microbiome in stool samples and circulating metabolites and lipids profiling using mass spectroscopy in plasma of all participants were performed. Functional mapping of microbiome into Gut Metabolic Modules (GMMs) was done using EggNog and KEGG databases. Multiomics integration was performed using MOFA tool. Results: Measures of the gut bacterial beta diversity differed significantly between T1D and HC, either with moderately or severely increased albuminuria. Taxonomic analyses of the bacterial microbiota identified 51 species that differed in absolute abundance between T1D and HC (17 higher, 34 lower). Stratified on levels of albuminuria, 10 species were differentially abundant for the moderately increased albuminuria group, 63 for the severely increased albuminuria group while 25 were common and differentially abundant both for moderately and severely increased albuminuria groups, when compared to HC. Functional characterization of the bacteriome identified 23 differentially enriched GMMs between T1D and HC, mostly involved in sugar and amino acid metabolism. No differences in relation to albuminuria stratification was observed. Twenty-five phages were differentially abundant between T1D and HC groups. Six of these varied with albuminuria status. Plasma metabolomics indicated differences in the steroidogenesis and sugar metabolism and circulating sphingolipids in T1D individuals. We identified association between sphingolipid levels and Bacteroides sp. abundances. MOFA revealed reduced interactions between gut microbiome and plasma metabolome profiles albeit polar metabolite, lipids and bacteriome compositions contributed to the variance in albuminuria levels among T1D individuals. Conclusions: Individuals with T1D and progressive kidney disease stratified on levels of albuminuria show distinct signatures in their gut microbiome and blood metabolome.

The novel inflammatory biomarker GlycA and triglyceride-rich lipoproteins are associated with the presence of subclinical myocardial dysfunction in subjects with type 1 diabetes mellitus.

BACKGROUND: Subjects with Type 1 diabetes mellitus (T1DM) have an increased incidence of heart failure (HF). Several pathophysiological mechanisms have been involved in its development. The aim of this study was to analyze the potential contribution of the advanced lipoprotein profile and plasma glycosylation (GlycA) to the presence of subclinical myocardial dysfunction in subjects with T1DM. METHODS: We included subjects from a Danish cohort of T1DM subjects (Thousand & 1 study) with either diastolic and/or systolic subclinical myocardial dysfunction, and a control group without myocardial dysfunction, matched by age, sex and HbA1c. All underwent a transthoracic echocardiogram and an advanced lipoprotein profile obtained by using the NMR-based Liposcale® test. GlycA NMR signal was also analyzed. Systolic dysfunction was defined as left ventricular ejection fraction ≤ 45% and diastolic dysfunction was considered as E/e'≥12 or E/e' 8-12 + volume of the left atrium > 34 ml/m2. To identify a metabolic profile associated with the presence of subclinical myocardial dysfunction, a multivariate supervised model of classification based on least squares regression (PLS-DA regression) was performed. RESULTS: One-hundred forty-six subjects had diastolic dysfunction and 18 systolic dysfunction. Compared to the control group, patients with myocardial dysfunction had longer duration of diabetes (p = 0.005), and higher BMI (p = 0.013), serum NTproBNP concentration (p = 0.001), systolic blood pressure (p < 0.001), albuminuria (p < 0.001), and incidence of advanced retinopathy (p < 0.001). The supervised classification model identified a specific pattern associated with myocardial dysfunction, with a capacity to discriminate patients with myocardial dysfunction from controls. PLS-DA showed that triglyceride-rich lipoproteins (TGRLs), such as VLDL (total VLDL particles, large VLDL subclass and VLDL-TG content) and IDL (IDL cholesterol content), as well as the plasma concentration of GlycA, were associated with the presence of subclinical myocardial dysfunction. CONCLUSION: Proatherogenic TGRLs and the proinflammatory biomarker Glyc A are strongly associated to myocardial dysfunction in T1DM. These findings suggest a pivotal role of TGRLs and systemic inflammation in the development of subclinical myocardial dysfunction in T1DM.

A saturated map of common genetic variants associated with human height.

Common single-nucleotide polymorphisms (SNPs) are predicted to collectively explain 40-50% of phenotypic variation in human height, but identifying the specific variants and associated regions requires huge sample sizes1. Here, using data from a genome-wide association study of 5.4 million individuals of diverse ancestries, we show that 12,111 independent SNPs that are significantly associated with height account for nearly all of the common SNP-based heritability. These SNPs are clustered within 7,209 non-overlapping genomic segments with a mean size of around 90 kb, covering about 21% of the genome. The density of independent associations varies across the genome and the regions of increased density are enriched for biologically relevant genes. In out-of-sample estimation and prediction, the 12,111 SNPs (or all SNPs in the HapMap 3 panel2) account for 40% (45%) of phenotypic variance in populations of European ancestry but only around 10-20% (14-24%) in populations of other ancestries. Effect sizes, associated regions and gene prioritization are similar across ancestries, indicating that reduced prediction accuracy is likely to be explained by linkage disequilibrium and differences in allele frequency within associated regions. Finally, we show that the relevant biological pathways are detectable with smaller sample sizes than are needed to implicate causal genes and variants. Overall, this study provides a comprehensive map of specific genomic regions that contain the vast majority of common height-associated variants. Although this map is saturated for populations of European ancestry, further research is needed to achieve equivalent saturation in other ancestries.

Genome-wide association analyses of physical activity and sedentary behavior provide insights into underlying mechanisms and roles in disease prevention.

Although physical activity and sedentary behavior are moderately heritable, little is known about the mechanisms that influence these traits. Combining data for up to 703,901 individuals from 51 studies in a multi-ancestry meta-analysis of genome-wide association studies yields 99 loci that associate with self-reported moderate-to-vigorous intensity physical activity during leisure time (MVPA), leisure screen time (LST) and/or sedentary behavior at work. Loci associated with LST are enriched for genes whose expression in skeletal muscle is altered by resistance training. A missense variant in ACTN3 makes the alpha-actinin-3 filaments more flexible, resulting in lower maximal force in isolated type IIA muscle fibers, and possibly protection from exercise-induced muscle damage. Finally, Mendelian randomization analyses show that beneficial effects of lower LST and higher MVPA on several risk factors and diseases are mediated or confounded by body mass index (BMI). Our results provide insights into physical activity mechanisms and its role in disease prevention.

Non-dipping and higher nocturnal blood pressure are associated with risk of mortality and development of kidney disease in type 1 diabetes.

AIMS: People with type 1 diabetes have increased risk of cardiovascular (CV) and kidney disease. A 24-hour ambulatory blood pressure (BP) measurement (ABPM) examines diurnal variations in BP. We aimed to determine the prognostic significance of blunted decrease in nocturnal systolic BP of <10 % (non-dipping of nocturnal BP) for CV- and kidney disease and all-cause mortality in type 1 diabetes. METHODS: From 2009 to 2011, at Steno Diabetes Center Copenhagen, 654 participants with type 1 diabetes had 24-hour ABPM obtained with a tonometric wrist-watch device (BPro, HealthStats, Singapore). In 2017, outcomes (composite CV endpoint; all-cause mortality; decline in estimated glomerular filtration rate (eGFR) ≥30 %; end-stage kidney disease (ESKD); and a composite kidney endpoint including decline in eGFR ≥30 %, ESKD and all-cause mortality) were registered. Hazard Ratios (HR) were calculated using Cox regressions. RESULTS: Participants were mean ± SD 55 ± 13 years old and had median (IQR) 35 (24-44) years diabetes duration. Mean daytime and nocturnal systolic BP were 133 ± 16 and 121 ± 16 mmHg while 337 (52 %) participants demonstrated non-dipping. After CV risk factor adjustments, non-dipping was associated with all-cause mortality (HR 2.12 (1.09-4.11), p = 0.03) and the composite kidney endpoint (HR 1.92 (1.23-3.00), p = 0.004). CONCLUSIONS: Non-dipping entailed increased risk of all-cause mortality and kidney disease in type 1 diabetes.

A genome-wide association study of total child psychiatric problems scores.

Substantial genetic correlations have been reported across psychiatric disorders and numerous cross-disorder genetic variants have been detected. To identify the genetic variants underlying general psychopathology in childhood, we performed a genome-wide association study using a total psychiatric problem score. We analyzed 6,844,199 common SNPs in 38,418 school-aged children from 20 population-based cohorts participating in the EAGLE consortium. The SNP heritability of total psychiatric problems was 5.4% (SE = 0.01) and two loci reached genome-wide significance: rs10767094 and rs202005905. We also observed an association of SBF2, a gene associated with neuroticism in previous GWAS, with total psychiatric problems. The genetic effects underlying the total score were shared with common psychiatric disorders only (attention-deficit/hyperactivity disorder, anxiety, depression, insomnia) (rG > 0.49), but not with autism or the less common adult disorders (schizophrenia, bipolar disorder, or eating disorders) (rG < 0.01). Importantly, the total psychiatric problem score also showed at least a moderate genetic correlation with intelligence, educational attainment, wellbeing, smoking, and body fat (rG > 0.29). The results suggest that many common genetic variants are associated with childhood psychiatric symptoms and related phenotypes in general instead of with specific symptoms. Further research is needed to establish causality and pleiotropic mechanisms between related traits.

Genetic loci and prioritization of genes for kidney function decline derived from a meta-analysis of 62 longitudinal genome-wide association studies.

Estimated glomerular filtration rate (eGFR) reflects kidney function. Progressive eGFR-decline can lead to kidney failure, necessitating dialysis or transplantation. Hundreds of loci from genome-wide association studies (GWAS) for eGFR help explain population cross section variability. Since the contribution of these or other loci to eGFR-decline remains largely unknown, we derived GWAS for annual eGFR-decline and meta-analyzed 62 longitudinal studies with eGFR assessed twice over time in all 343,339 individuals and in high-risk groups. We also explored different covariate adjustment. Twelve genome-wide significant independent variants for eGFR-decline unadjusted or adjusted for eGFR-baseline (11 novel, one known for this phenotype), including nine variants robustly associated across models were identified. All loci for eGFR-decline were known for cross-sectional eGFR and thus distinguished a subgroup of eGFR loci. Seven of the nine variants showed variant-by-age interaction on eGFR cross section (further about 350,000 individuals), which linked genetic associations for eGFR-decline with age-dependency of genetic cross-section associations. Clinically important were two to four-fold greater genetic effects on eGFR-decline in high-risk subgroups. Five variants associated also with chronic kidney disease progression mapped to genes with functional in-silico evidence (UMOD, SPATA7, GALNTL5, TPPP). An unfavorable versus favorable nine-variant genetic profile showed increased risk odds ratios of 1.35 for kidney failure (95% confidence intervals 1.03-1.77) and 1.27 for acute kidney injury (95% confidence intervals 1.08-1.50) in over 2000 cases each, with matched controls). Thus, we provide a large data resource, genetic loci, and prioritized genes for kidney function decline, which help inform drug development pipelines revealing important insights into the age-dependency of kidney function genetics.

Cardiovascular Autonomic Neuropathy in Type 1 Diabetes Is Associated With Disturbances in TCA, Lipid, and Glucose Metabolism.

Introduction: Diabetic cardiovascular autonomic neuropathy (CAN) is associated with increased mortality and morbidity. To explore metabolic mechanisms associated with CAN we investigated associations between serum metabolites and CAN in persons with type 1 diabetes (T1D). Materials and Methods: Cardiovascular reflex tests (CARTs) (heart rate response to: deep breathing; lying-to-standing test; and the Valsalva maneuver) were used to diagnose CAN in 302 persons with T1D. More than one pathological CARTs defined the CAN diagnosis. Serum metabolomics and lipidomic profiles were analyzed with two complementary non-targeted mass-spectrometry methods. Cross-sectional associations between metabolites and CAN were assessed by linear regression models adjusted for relevant confounders. Results: Participants were median (IQR) aged 55(49, 63) years, 48% males with diabetes duration 39(32, 47) years, HbA1c 63(55,69) mmol/mol and 34% had CAN. A total of 75 metabolites and 106 lipids were analyzed. In crude models, the CAN diagnosis was associated with higher levels of hydroxy fatty acids (2,4- and 3,4-dihydroxybutanoic acids, 4-deoxytetronic acid), creatinine, sugar derivates (ribitol, ribonic acid, myo-inositol), citric acid, glycerol, phenols, phosphatidylcholines and lower levels of free fatty acids and the amino acid methionine (p<0.05). Upon adjustment, positive associations with the CAN diagnoses were retained for hydroxy fatty acids, tricarboxylic acid (TCA) cycle-based sugar derivates, citric acid, and phenols (P<0.05). Conclusion: Metabolic pathways, including the TCA cycle, hydroxy fatty acids, phosphatidylcholines and sugar derivatives are associated with the CAN diagnosis in T1D. These pathway may be part of the pathogeneses leading to CAN and may be modifiable risk factors for the complication.

Early-life respiratory tract infections and the risk of school-age lower lung function and asthma: a meta-analysis of 150 000 European children.

BACKGROUND: Early-life respiratory tract infections might affect chronic obstructive respiratory diseases, but conclusive studies from general populations are lacking. Our objective was to examine if children with early-life respiratory tract infections had increased risks of lower lung function and asthma at school age. METHODS: We used individual participant data of 150 090 children primarily from the EU Child Cohort Network to examine the associations of upper and lower respiratory tract infections from age 6 months to 5 years with forced expiratory volume in 1 s (FEV1), forced vital capacity (FVC), FEV1/FVC, forced expiratory flow at 75% of FVC (FEF75%) and asthma at a median (range) age of 7 (4-15) years. RESULTS: Children with early-life lower, not upper, respiratory tract infections had a lower school-age FEV1, FEV1/FVC and FEF75% (z-score range: -0.09 (95% CI -0.14- -0.04) to -0.30 (95% CI -0.36- -0.24)). Children with early-life lower respiratory tract infections had a higher increased risk of school-age asthma than those with upper respiratory tract infections (OR range: 2.10 (95% CI 1.98-2.22) to 6.30 (95% CI 5.64-7.04) and 1.25 (95% CI 1.18-1.32) to 1.55 (95% CI 1.47-1.65), respectively). Adjustment for preceding respiratory tract infections slightly decreased the strength of the effects. Observed associations were similar for those with and without early-life wheezing as a proxy for early-life asthma. CONCLUSIONS: Our findings suggest that early-life respiratory tract infections affect development of chronic obstructive respiratory diseases in later life, with the strongest effects for lower respiratory tract infections.

Genome-wide study of early and severe childhood asthma identifies interaction between CDHR3 and GSDMB.

BACKGROUND: Asthma with severe exacerbation is one of the most common causes of hospitalization among young children. Exacerbations are typically triggered by respiratory infections, but the host factors causing recurrent infections and exacerbations in some children are poorly understood. As a result, current treatment options and preventive measures are inadequate. OBJECTIVE: We sought to identify genetic interaction associated with the development of childhood asthma. METHODS: We performed an exhaustive search for pairwise interaction between genetic single nucleotide polymorphisms using 1204 cases of a specific phenotype of early childhood asthma with severe exacerbations in patients aged 2 to 6 years combined with 5328 nonasthmatic controls. Replication was attempted in 3 independent populations, and potential underlying immune mechanisms were investigated in the COPSAC2010 and COPSAC2000 birth cohorts. RESULTS: We found evidence of interaction, including replication in independent populations, between the known childhood asthma loci CDHR3 and GSDMB. The effect of CDHR3 was dependent on the GSDMB genotype, and this interaction was more pronounced for severe and early onset of disease. Blood immune analyses suggested a mechanism related to increased IL-17A production after viral stimulation. CONCLUSIONS: We found evidence of interaction between CDHR3 and GSDMB in development of early childhood asthma, possibly related to increased IL-17A response to viral infections. This study demonstrates the importance of focusing on specific disease subtypes for understanding the genetic mechanisms of asthma.

Urinary Proteomics Identifies Cathepsin D as a Biomarker of Rapid eGFR Decline in Type 1 Diabetes.

OBJECTIVE: Understanding mechanisms underlying rapid estimated glomerular filtration rate (eGFR) decline is important to predict and treat kidney disease in type 1 diabetes (T1D). RESEARCH DESIGN AND METHODS: We performed a case-control study nested within four T1D cohorts to identify urinary proteins associated with rapid eGFR decline. Case and control subjects were categorized based on eGFR decline ≥3 and <1 mL/min/1.73 m2/year, respectively. We used targeted liquid chromatography-tandem mass spectrometry to measure 38 peptides from 20 proteins implicated in diabetic kidney disease. Significant proteins were investigated in complementary human cohorts and in mouse proximal tubular epithelial cell cultures. RESULTS: The cohort study included 1,270 participants followed a median 8 years. In the discovery set, only cathepsin D peptide and protein were significant on full adjustment for clinical and laboratory variables. In the validation set, associations of cathepsin D with eGFR decline were replicated in minimally adjusted models but lost significance with adjustment for albuminuria. In a meta-analysis with combination of discovery and validation sets, the odds ratio for the association of cathepsin D with rapid eGFR decline was 1.29 per SD (95% CI 1.07-1.55). In complementary human cohorts, urine cathepsin D was associated with tubulointerstitial injury and tubulointerstitial cathepsin D expression was associated with increased cortical interstitial fractional volume. In mouse proximal tubular epithelial cell cultures, advanced glycation end product-BSA increased cathepsin D activity and inflammatory and tubular injury markers, which were further increased with cathepsin D siRNA. CONCLUSIONS: Urine cathepsin D is associated with rapid eGFR decline in T1D and reflects kidney tubulointerstitial injury.