A framework for conducting time-varying genome-wide association studies: An application to body mass index across childhood in six multiethnic cohorts.

Genetic effects on changes in human traits over time are understudied and may have important pathophysiological impact. We propose a framework that enables data quality control, implements mixed models to evaluate trajectories of change in traits, and estimates phenotypes to identify age-varying genetic effects in genome-wide association studies (GWASs). Using childhood body mass index (BMI) as an example, we included 71,336 participants from six cohorts and estimated the slope and area under the BMI curve within four time periods (infancy, early childhood, late childhood and adolescence) for each participant, in addition to the age and BMI at the adiposity peak and the adiposity rebound. GWAS on each of the estimated phenotypes identified 28 genome-wide significant variants at 13 loci across the 12 estimated phenotypes, one of which was novel (in DAOA) and had not been previously associated with childhood or adult BMI. Genetic studies of changes in human traits over time could uncover novel biological mechanisms influencing quantitative traits.

Pathogenic, Total Loss-of-Function DYRK1B Variants Cause Monogenic Obesity Associated With Type 2 Diabetes.

OBJECTIVE: Rare variants in DYRK1B have been described in some patients with central obesity, type 2 diabetes, and early-onset coronary disease. Owing to the limited number of conducted studies, the broader impact of DYRK1B variants on a larger scale has yet to be investigated. RESEARCH DESIGN AND METHODS: DYRK1B was sequenced in 9,353 participants from a case-control study for obesity and type 2 diabetes. Each DYRK1B variant was functionally assessed in vitro. Variant pathogenicity was determined using criteria from the American College of Medical Genetics and Genomics (ACMG). The effect of pathogenic or likely pathogenic (P/LP) variants on metabolic traits was assessed using adjusted mixed-effects score tests. RESULTS: Sixty-five rare, heterozygous DYRK1B variants were identified and were not associated with obesity or type 2 diabetes. Following functional analyses, 20 P/LP variants were pinpointed, including 6 variants that exhibited a fully inhibitory effect (P/LP-null) on DYRK1B activity. P/LP and P/LP-null DYRK1B variants were associated with increased BMI and obesity risk; however, the impact was notably more pronounced for the P/LP-null variants (effect of 8.0 ± 3.2 and odds ratio of 7.9 [95% CI 1.2-155]). Furthermore, P/LP-null variants were associated with higher fasting glucose and type 2 diabetes risk (effect of 2.9 ± 1.0 and odds ratio of 4.8 [95% CI 0.85-37]), while P/LP variants had no effect on glucose homeostasis. CONCLUSIONS: P/LP, total loss-of-function DYRK1B variants cause monogenic obesity associated with type 2 diabetes. This study underscores the significance of conducting functional assessments in order to accurately ascertain the tangible effects of P/LP DYRK1B variants.

High morbidity and mortality in children with untreated congenital deficiency of leptin or its receptor.

The long-term clinical outcomes of severe obesity due to leptin signaling deficiency are unknown. We carry out a retrospective cross-sectional investigation of a large cohort of children with leptin (LEP), LEP receptor (LEPR), or melanocortin 4 receptor (MC4R) deficiency (n = 145) to evaluate the progression of the disease. The affected individuals undergo physical, clinical, and metabolic evaluations. We report a very high mortality in children with LEP (26%) or LEPR deficiency (9%), mainly due to severe pulmonary and gastrointestinal infections. In addition, 40% of surviving children with LEP or LEPR deficiency experience life-threatening episodes of lung or gastrointestinal infections. Although precision drugs are currently available for LEP and LEPR deficiencies, as yet, they are not accessible in Pakistan. An appreciation of the severe impact of LEP or LEPR deficiency on morbidity and early mortality, educational attainment, and the attendant stigmatization should spur efforts to deliver the available life-saving drugs to these children as a matter of urgency.

Identification of myokines susceptible to improve glucose homeostasis after bariatric surgery.

IMPORTANCE AND OBJECTIVE: The identification of myokines susceptible to improve glucose homeostasis following bariatric surgery could lead to new therapeutic approaches for type 2 diabetes. METHODS: Changes in the homeostasis model assessment (HOMA) test were assessed in patients before and 3 months after bariatric surgery. Changes in myokines expression and circulating levels were assessed using real-time quantitative polymerase chain reaction (RT-qPCR) and enzyme-linked immunosorbent assay (ELISA). Myokines known to regulate glucose homeostasis were identified using literature (targeted study) and putative myokines using RNA-sequencing (untargeted study). A linear regression analysis adjusted for age and sex was used to search for associations between changes in the HOMA test and changes in myokines. RESULTS: In the targeted study, brain-derived neurotrophic factor (BDNF) expression was upregulated (+30%, P = .006) while BDNF circulating levels were decreased (-12%, P = .001). Upregulated BDNF expression was associated with decreased HOMA of insulin resistance (HOMA-IR) (adjusted estimate [95% confidence interval {CI}]: -0.51 [-0.88 to -0.13], P = .010). Decreased BDNF serum levels were associated with decreased HOMA of beta-cell function (HOMA-B) (adjusted estimate [95% CI] = 0.002 [0.00002-0.0031], P = .046). In the untargeted study, upregulated putative myokines included XYLT1 (+64%, P < .001), LGR5 (+57, P< .001), and SPINK5 (+46%, P < .001). Upregulated LGR5 was associated with decreased HOMA-IR (adjusted estimate [95% CI] = -0.50 [-0.86 to -0.13], P = .009). Upregulated XYLT1 and SPINK5 were associated with increased HOMA of insulin sensitivity (HOMA-S) (respectively, adjusted estimate [95% CI] = 109.1 [28.5-189.8], P = .009 and 16.5 [0.87-32.19], P = .039). CONCLUSIONS: Improved glucose homeostasis following bariatric surgery is associated with changes in myokines expression and circulating levels. In particular, upregulation of BDNF, XYLT1, SPINK5, and LGR5 is associated with improved insulin sensitivity. These results suggest that these myokines could contribute to improved glucose homeostasis following bariatric surgery. STUDY REGISTRATION: NCT03341793 on ClinicalTrials.gov (https://clinicaltrials.gov/).

Biallelic Mutations in P4HTM Cause Syndromic Obesity.

We previously demonstrated that 50% of children with obesity from consanguineous families from Pakistan carry pathogenic variants in known monogenic obesity genes. Here, we have discovered a novel monogenetic recessive form of severe childhood obesity using an in-house computational staged approach. The analysis included whole-exome sequencing data of 366 children with severe obesity, 1,000 individuals of the Pakistan Risk of Myocardial Infarction Study (PROMIS) study, and 200,000 participants of the UK Biobank to prioritize genes harboring rare homozygous variants with putative effect on human obesity. We identified five rare or novel homozygous missense mutations predicted deleterious in five consanguineous families in P4HTM encoding prolyl 4-hydroxylase transmembrane (P4H-TM). We further found two additional homozygous missense mutations in children with severe obesity of Indian and Moroccan origin. Molecular dynamics simulation suggested that these mutations destabilized the active conformation of the substrate binding domain. Most carriers also presented with hypotonia, cognitive impairment, and/or developmental delay. Three of the five probands died of pneumonia during the first 2 years of the follow-up. P4HTM deficiency is a novel form of syndromic obesity, affecting 1.5% of our children with obesity associated with high mortality. P4H-TM is a hypoxia-inducible factor that is necessary for survival and adaptation under oxygen deprivation, but the role of this pathway in energy homeostasis and obesity pathophysiology remains to be elucidated.

Epigenetic changes associated with hyperglycaemia exposure in the longitudinal D.E.S.I.R. cohort.

AIM: - Understanding DNA methylation dynamics associated with progressive hyperglycaemia exposure could provide early diagnostic biomarkers and an avenue for delaying type 2 diabetes mellitus (T2DM). We aimed to identify DNA methylation changes during a 6-year period associated with early hyperglycaemia exposure using the longitudinal D.E.S.I.R. METHODS: - We selected individuals with progressive hyperglycaemia exposure based on T2DM diagnostic criteria: 27 with long-term exposure, 34 with short-term exposure and 34 normoglycaemic controls. DNA from blood at inclusion and at the 6-year visit was subjected to methylation analysis using 850K methylation-EPIC arrays. A linear mixed model was used to perform an epigenome-wide association study (EWAS) and identify methylated changes associated with hyperglycaemia exposure during a 6-year time-period. RESULTS: - We did not identify differentially methylated sites that reached false discovery rate (FDR)-significance in our cohort. Based on EWAS, we focused our analysis on methylation sites that had a constant effect during the 6 years across the hyperglycaemia groups compared to controls and found the most statistically significant site was the reported cg19693031 probe (TXNIP). We also performed an EWAS with HbA1c, using the inclusion and the 6-year methylation data and did not identify any FDR-significant CpGs. CONCLUSIONS: - Our study reveals that DNA methylation changes are not robustly associated with hyperglycaemia exposure or HbA1c during a short-term period, however, our top loci indicate potential interest and should be replicated in larger cohorts.

Rare Variant Analysis of Obesity-Associated Genes in Young Adults With Severe Obesity From a Consanguineous Population of Pakistan.

Recent advances in genetic analysis have significantly helped in progressively attenuating the heritability gap of obesity and have brought into focus monogenic variants that disrupt the melanocortin signaling. In a previous study, next-generation sequencing revealed a monogenic etiology in ∼50% of the children with severe obesity from a consanguineous population in Pakistan. Here we assess rare variants in obesity-causing genes in young adults with severe obesity from the same region. Genomic DNA from 126 randomly selected young adult obese subjects (BMI 37.2 ± 0.3 kg/m2; age 18.4 ± 0.3 years) was screened by conventional or augmented whole-exome analysis for point mutations and copy number variants (CNVs). Leptin, insulin, and cortisol levels were measured by ELISA. We identified 13 subjects carrying 13 different pathogenic or likely pathogenic variants in LEPR, PCSK1, MC4R, NTRK2, POMC, SH2B1, and SIM1. We also identified for the first time in the human, two homozygous stop-gain mutations in ASNSD1 and IFI16 genes. Inactivation of these genes in mouse models has been shown to result in obesity. Additionally, we describe nine homozygous mutations (seven missense, one stop-gain, and one stop-loss) and four copy-loss CNVs in genes or genomic regions previously linked to obesity-associated traits by genome-wide association studies. Unexpectedly, in contrast to obese children, pathogenic mutations in LEP and LEPR were either absent or rare in this cohort of young adults. High morbidity and mortality risks and social disadvantage of children with LEP or LEPR deficiency may in part explain this difference between the two cohorts.