Modeling biological age using blood biomarkers and physical measurements in Chinese adults

Background This study aimed to: 1) assess the associations of biological age acceleration based on Klemera and Doubal's method (KDM-AA) with long-term risk of all-cause mortality; and 2) compare the association of KDM-AA with all-cause mortality among participants potentially at different stages of the cardiovascular disease (CVD) continuum. Methods The present study was based on a subpopulation of the China Kadoorie Biobank, with baseline survey during 2004–08. A total of 12,377 participants free of ischemic heart disease, stroke, or cancer at baseline were included, in which 8180 participants were identified to develop major coronary event (MCE), ischemic stroke (IS), intracerebral hemorrhage (ICH) or subarachnoid hemorrhage (SAH), and 4197 remained free of these cardiovascular diseases before 1 January 2014. These participants were followed up until 1 Jan 2018. KDM-AA was calculated by regressing biological age measurement, which was constructed based on baseline 16 physical and 9 biochemical markers using Klemera and Doubal's method, on chronological age. We estimated the associations of KDM-AA with the mortality risk using the hazard ratio (HR) and 95% confidence interval (CI) from Cox proportional hazard models. We assessed discrimination performance by Harrell's C-index and net reclassification index (NRI). Findings The participants who developed MCE (mean KDM-AA = 0.1 year, standard deviation [SD] = 1.6 years) or ICH/SAH (0.3 ± 1.5 years) during subsequent follow-up showed accelerated aging at baseline compared to those of IS (0.0 ± 1.2 years) and control (−0.3 ± 1.3 years) groups. The KDM-AA was positively associated with long-term risk of all-cause mortality (HR = 1.20; 95% CI: 1.17, 1.23), and the association was robust for participants potentially at different stages of the CVD continuum. Adding KDM-AA improved mortality prediction compared to the model only with sociodemographic and lifestyle factors in whole participants, with the Harrell's C-index increasing from 0.813 (0.807, 0.819) to 0.821 (0.815, 0.826) (NRI = 0.011; 95% CI: 0.003, 0.019). Interpretation In this middle-aged and elderly Chinese population, the KDM-AA is a promising measurement for biological age, and can capture the difference in cardiovascular health and predict the risk of all-cause mortality over a decade. Funding This work was supported by National Natural Science Foundation of China (82192904, 82192901, 82192900, 81941018). The CKB baseline survey and the first re-survey were supported by a grant from the Kadoorie Charitable Foundation Hong Kong. The long-term follow-up is supported by grants from the UK Wellcome Trust (212946/Z/18/Z, 202922/Z/16/Z, 104085/Z/14/Z, 088158/Z/09/Z), grants (2016YFC0900500) from the National Key R&D Program of China, National Natural Science Foundation of China (81390540, 91846303), and Chinese Ministry of Science and Technology (2011BAI09B01)

Healthy lifestyle, DNA methylation age acceleration, and incident risk of coronary heart disease

Background DNA methylation clocks emerged as a tool to determine biological aging and have been related to mortality and age-related diseases. Little is known about the association of DNA methylation age (DNAm age) with coronary heart disease (CHD), especially in the Asian population. Results Methylation level of baseline blood leukocyte DNA was measured by Infinium Methylation EPIC BeadChip for 491 incident CHD cases and 489 controls in the prospective China Kadoorie Biobank. We calculated the methylation age using a prediction model developed among Chinese. The correlation between chronological age and DNAm age was 0.90. DNA methylation age acceleration (Δage) was defined as the residual of regressing DNA methylation age on the chronological age. After adjustment for multiple risk factors of CHD and cell type proportion, compared with participants in the bottom quartile of Δage, the OR (95% CI) for CHD was 1.84 (1.17, 2.89) for participants in the top quartile. One SD increment in Δage was associated with 30% increased risk of CHD (OR = 1.30; 95% CI 1.09, 1.56; Ptrend = 0.003). The average number of cigarette equivalents consumed per day and waist-to-hip ratio were positively associated with Δage; red meat consumption was negatively associated with Δage, characterized by accelerated aging in those who never or rarely consumed red meat (all P  Conclusions We first identified the association between DNAm age acceleration and incident CHD in the Asian population, and provided evidence that unfavorable lifestyle-induced epigenetic aging may play an important part in the underlying pathway to CHD

Understanding the genetic complexity of puberty timing across the allele frequency spectrum

Pubertal timing varies considerably and is associated with later health outcomes. We performed multi-ancestry genetic analyses on ~800,000 women, identifying 1,080 signals for age at menarche. Collectively, these explained 11% of trait variance in an independent sample. Women at the top and bottom 1% of polygenic risk exhibited ~11 and ~14-fold higher risks of delayed and precocious puberty, respectively. We identified several genes harboring rare loss-of-function variants in ~200,000 women, including variants in ZNF483, which abolished the impact of polygenic risk. Variant-to-gene mapping approaches and mouse gonadotropin-releasing hormone neuron RNA sequencing implicated 665 genes, including an uncharacterized G-protein-coupled receptor, GPR83, which amplified the signaling of MC3R, a key nutritional sensor. Shared signals with menopause timing at genes involved in DNA damage response suggest that the ovarian reserve might signal centrally to trigger puberty. We also highlight body size-dependent and independent mechanisms that potentially link reproductive timing to later life disease

Associations of polygenic risk scores with risks of stroke and its subtypes in Chinese

Background and purpose Previous studies, mostly focusing on the European population, have reported polygenic risk scores (PRSs) might achieve risk stratification of stroke. We aimed to examine the association strengths of PRSs with risks of stroke and its subtypes in the Chinese population. Methods Participants with genome-wide genotypic data in China Kadoorie Biobank were split into a potential training set (n=22 191) and a population-based testing set (n=72 150). Four previously developed PRSs were included, and new PRSs for stroke and its subtypes were developed. The PRSs showing the strongest association with risks of stroke or its subtypes in the training set were further evaluated in the testing set. Cox proportional hazards regression models were used to estimate the association strengths of different PRSs with risks of stroke and its subtypes (ischaemic stroke (IS), intracerebral haemorrhage (ICH) and subarachnoid haemorrhage (SAH)). Results In the testing set, during 872 919 person-years of follow-up, 8514 incident stroke events were documented. The PRSs of any stroke (AS) and IS were both positively associated with risks of AS, IS and ICH (p<0.05). The HR for per SD increment (HRSD) of PRSAS was 1.10 (95% CI 1.07 to 1.12), 1.10 (95% CI 1.07 to 1.12) and 1.13 (95% CI 1.07 to 1.20) for AS, IS and ICH, respectively. The corresponding HRSD of PRSIS was 1.08 (95% CI 1.06 to 1.11), 1.08 (95% CI 1.06 to 1.11) and 1.09 (95% CI 1.03 to 1.15). PRSICH was positively associated with the risk of ICH (HRSD=1.07, 95% CI 1.01 to 1.14). PRSSAH was not associated with risks of stroke and its subtypes. The addition of current PRSs offered little to no improvement in stroke risk prediction and risk stratification. Conclusions In this Chinese population, the association strengths of current PRSs with risks of stroke and its subtypes were moderate, suggesting a limited value for improving risk prediction over traditional risk factors in the context of current genome-wide association study under-representing the East Asian population

Alcohol intake and cause-specific mortality: conventional and genetic evidence in a prospective cohort study of 512 000 adults in China

Background Genetic variants that affect alcohol use in East Asian populations could help assess the causal effects of alcohol consumption on cause-specific mortality. We aimed to investigate the associations between alcohol intake and cause-specific mortality using conventional and genetic epidemiological methods among more than 512 000 adults in China. Methods The prospective China Kadoorie Biobank cohort study enrolled 512 724 adults (210 205 men and 302 519 women) aged 30–79 years, during 2004–08. Residents with no major disabilities from ten diverse urban and rural areas of China were invited to participate, and alcohol use was self-reported. During 12 years of follow-up, 56 550 deaths were recorded through linkage to death registries, including 23 457 deaths among 168 050 participants genotyped for ALDH2-rs671 and ADH1B-rs1229984. Adjusted hazard ratios (HRs) for cause-specific mortality by self-reported and genotype-predicted alcohol intake were estimated using Cox regression. Findings 33% of men drank alcohol most weeks. In conventional observational analyses, ex-drinkers, non-drinkers, and heavy drinkers had higher risks of death from most major causes than moderate drinkers. Among current drinkers, each 100 g/week higher alcohol intake was associated with higher mortality risks from cancers (HR 1·18 [95% CI 1·14−1·22]), cardiovascular disease (CVD; HR 1·19 [1·15−1·24]), liver diseases (HR 1·51 [1·27−1·78]), non-medical causes (HR 1·15 [1·08−1·23]), and all causes (HR 1·18 [1·15−1·20]). In men, ALDH2-rs671 and ADH1B-rs1229984 genotypes predicted 60-fold differences in mean alcohol intake (4 g/week in the lowest group vs 255 g/week in the highest). Genotype-predicted alcohol intake was uniformly and positively associated with risks of death from all causes (n=12 939; HR 1·07 [95% CI 1·05−1·10]) and from pre-defined alcohol-related cancers (n=1274; 1·12 [1·04−1·21]), liver diseases (n=110; 1·31 [1·02−1·69]), and CVD (n=6109; 1·15 [1·10−1·19]), chiefly due to stroke (n=3285; 1·18 [1·12–1·24]) rather than ischaemic heart disease (n=2363; 1·06 [0·99–1·14]). Results were largely consistent using a polygenic score to predict alcohol intake, with higher intakes associated with higher risks of death from alcohol-related cancers, CVD, and all causes. Approximately 2% of women were current drinkers, and although power was low to assess observational associations of alcohol with mortality, the genetic evidence suggested that the excess risks in men were due to alcohol, not pleiotropy. Interpretation Higher alcohol intake increased the risks of death overall and from major diseases for men in China. There was no genetic evidence of protection from moderate drinking for all-cause and cause-specific mortality, including CVD. Funding Kadoorie Charitable Foundation, National Natural Science Foundation of China, British Heart Foundation, Cancer Research UK, GlaxoSmithKline, Wellcome Trust, Medical Research Council, and Chinese Ministry of Science and Technology

Genetic and healthy lifestyle factors in relation to the incidence and prognosis of severe liver disease in the Chinese population

Background: Severe liver disease (SLD), including cirrhosis and liver cancer, constitutes a major disease burden in China. We aimed to examine the association of genetic and healthy lifestyle factors with the incidence and prognosis of SLD. Methods: The study population included 504,009 participants from the prospective China Kadoorie Biobank aged 30–79 years. The individuals were from 10 diverse areas in China without a history of cancer or liver disease at baseline. Cox regression was used to estimate adjusted hazard ratios (HRs) for incident SLD and death after SLD diagnosis associated with healthy lifestyle factors (smoking, alcohol, physical activity, and central adiposity). Additionally, the contribution of genetic risk for hepatitis B virus (HBV, assessed by genetic variants in major histocompatibility complex, class II, DP/DQ [HLA-DP/DQ] genes) was also estimated. Results: Compared with those with 0–1 healthy lifestyle factor, participants with 2, 3, and 4 factors had 12% (HR 0.88 [95% confidence interval [CI] 0.85, 0.92]), 26% (HR 0.74 [95%CI: 0.69, 0.79]), and 44% (HR 0.56 [95%CI: 0.48, 0.65]) lower risks of SLD, respectively. Inverse associations were observed among participants with both low and high genetic risks (HR per 1-point increase 0.83 [95%CI: 0.74, 0.94] and 0.91 [95%CI: 0.82, 1.02], respectively; Pinteraction = 0.51), although with a non-significant trend among those with a high genetic risk. Inverse associations were also observed between healthy lifestyle factors and liver biomarkers regardless of the genetic risk. Despite the limited power, healthy lifestyle factors were associated with a lower risk of death after incident SLD among participants with a low genetic risk (HR 0.59 [95%CI: 0.37, 0.96]). Conclusions: Lifestyle modification may be beneficial in terms of lowering the risk of SLD regardless of the genetic risk. Moreover, it is also important for improving the prognosis of SLD in individuals with a low genetic risk. Future studies are warranted to examine the impact of healthy lifestyles on SLD prognosis, particularly among individuals with a high genetic risk

Associations of diabetes, circulating protein biomarkers, and risk of pancreatic cancer

Background: Type 2 diabetes (T2D) is associated with higher risk of pancreatic cancer (PC), but the underlying mechanisms are not fully understood. Methods: We conducted a case-subcohort study involving 610 PC cases and 623 subcohort participants with 92 protein biomarkers measured in baseline plasma samples. Genetically-instrumented T2D was derived using 86 single-nucleotide polymorphisms (SNPs), including insulin resistance (IR) SNPs. Results: In observational analyses of 623 subcohort participants (mean age, 52 years; 61% women), T2D was positively associated with 13 proteins (SD difference: IL6: 0.52 [0.23–0.81]; IL10: 0.41 [0.12–0.70]), of which 8 were nominally associated with incident PC. The 8 proteins potentially mediated 36.9% (18.7–75.0%) of the association between T2D and PC. In MR, no associations were observed for genetically-determined T2D with proteins, but there were positive associations of genetically-determined IR with IL6 and IL10 (SD difference: 1.23 [0.05–2.41] and 1.28 [0.31–2.24]). In two-sample MR, fasting insulin was associated with both IL6 and PC, but no association was observed between IL6 and PC. Conclusions: Proteomics were likely to explain the association between T2D and PC, but were not causal mediators. Elevated fasting insulin driven by insulin resistance might explain the associations of T2D, proteomics, and PC

Causal association between snoring and stroke: a Mendelian randomization study in a Chinese population

Background: Previous observational studies established a positive relationship between snoring and stroke. We aimed to investigate the causal effect of snoring on stroke. Methods: Based on 82,339 unrelated individuals with qualified genotyping data of Asian descent from the China Kadoorie Biobank (CKB), we conducted a Mendelian randomization (MR) analysis of snoring and stroke. Genetic variants identified in the genome-wide association analysis (GWAS) of snoring in CKB and UK Biobank (UKB) were selected for constructing genetic risk scores (GRS). A two-stage method was applied to estimate the associations of the genetically predicted snoring with stroke and its subtypes. Besides, MR analysis among the non-obese group (body mass index, BMI <24.0 kg/m2), as well as multivariable MR (MVMR), were performed to control for potential pleiotropy from BMI. In addition, the inverse-variance weighted (IVW) method was applied to estimate the causal association with genetic variants identified in CKB GWAS. Findings: Positive associations were found between snoring and total stroke, hemorrhagic stroke (HS), and ischemic stroke (IS). With GRS of CKB, the corresponding HRs (95% CIs) were 1.56 (1.15, 2.12), 1.50 (0.84, 2.69), 2.02 (1.36, 3.01), and the corresponding HRs (95% CIs) using GRS of UKB were 1.78 (1.30, 2.43), 1.94 (1.07, 3.52), and 1.74 (1.16, 2.61). The associations remained stable in the MR among the non-obese group, MVMR analysis, and MR analysis using the IVW method. Interpretation: This study suggests that, among Chinese adults, genetically predicted snoring could increase the risk of total stroke, IS, and HS, and the causal effect was independent of BMI

Association of vitamin D with risk of type 2 diabetes:A Mendelian randomisation study in European and Chinese adults

BACKGROUND:Observational studies have reported that higher plasma 25-hydroxyvitamin D (25[OH]D) concentrations are associated with lower risks of diabetes, but it is unclear if these associations are causal. The aim of this study was to test the relevance of 25(OH)D for type 2 diabetes using genetically instrumented differences in plasma 25(OH)D concentrations. METHODS AND FINDINGS:Data were available on four 25(OH)D single nucleotide polymorphisms (SNPs; n = 82,464), plasma 25(OH)D concentrations (n = 13,565), and cases with diabetes (n = 5,565) in the China Kadoorie Biobank (CKB). The effects on risk of diabetes were assessed by a genetic score using two 25(OH)D synthesis SNPs (DHCR7-rs12785878 and CYP2R1-rs10741657), with and without the addition of SNPs affecting the transport (GC/DBP-rs2282679) and catabolism (CYP24A1-rs6013897) of 25(OH)D. The CKB results were combined in a meta-analysis of 10 studies for the 2 synthesis SNPs (n = 58,312 cases) and 7 studies for all 4 SNPs (n = 32,796 cases). Mean (SD) 25(OH)D concentration was 62 (20) nmol/l in CKB, and the per allele effects of genetic scores on 25(OH)D were 2.87 (SE 0.39) for the synthesis SNPs and 3.54 (SE 0.32) for all SNPs. A 25-nmol/l higher biochemically measured 25(OH)D was associated with a 9% (95% CI: 0%-18%) lower risk of diabetes in CKB. In a meta-analysis of all studies, a 25-nmol/l higher genetically instrumented 25(OH)D concentration was associated with a 14% (95% CI: 3%-23%) lower risk of diabetes (p = 0.01) using the 2 synthesis SNPs. An equivalent difference in 25(OH)D using a genetic score with 4 SNPs was not significantly associated with diabetes (odds ratio 8%, 95% CI: -1% to 16%, lower risk, p = 0.07), but had some evidence of pleiotropy. A limitation of the meta-analysis was the access only to study level rather than individual level data. CONCLUSIONS:The concordant risks of diabetes for biochemically measured and genetically instrumented differences in 25(OH)D using synthesis SNPs provide evidence for a causal effect of higher 25(OH)D for prevention of diabetes