Background:Maternal smoking has been linked to adverse health outcomes in newborns but the extent to which it impacts newborn health has not been quantified through an aggregated cord blood DNA methylation (DNAm) score. Here, we examine the feasibility of using cord blood DNAm scores leveraging large external studies as discovery samples to capture the epigenetic signature of maternal smoking and its influence on newborns in White European and South Asian populations.Methods:We first examined the association between individual CpGs and cigarette smoking during pregnancy, and smoking exposure in two White European birth cohorts (n=744). Leveraging established CpGs for maternal smoking, we constructed a cord blood epigenetic score of maternal smoking that was validated in one of the European-origin cohorts (n=347). This score was then tested for association with smoking status, secondary smoking exposure during pregnancy, and health outcomes in offspring measured after birth in an independent White European (n=397) and a South Asian birth cohort (n=504).Results:Several previously reported genes for maternal smoking were supported, with the strongest and most consistent association signal from the GFI1 gene (6 CpGs with p<5 × 10-5). The epigenetic maternal smoking score was strongly associated with smoking status during pregnancy (OR = 1.09 [1.07, 1.10], p=5.5 × 10-33) and more hours of self-reported smoking exposure per week (1.93 [1.27, 2.58], p=7.8 × 10-9) in White Europeans. However, it was not associated with self-reported exposure (p>0.05) among South Asians, likely due to a lack of smoking in this group. The same score was consistently associated with a smaller birth size (–0.37±0.12 cm, p=0.0023) in the South Asian cohort and a lower birth weight (–0.043±0.013 kg, p=0.0011) in the combined cohorts.Conclusions:This cord blood epigenetic score can help identify babies exposed to maternal smoking and assess its long-term impact on growth. Notably, these results indicate a consistent association between the DNAm signature of maternal smoking and a small body size and low birth weight in newborns, in both White European mothers who exhibited some amount of smoking and in South Asian mothers who themselves were not active smokers.Funding:This study was funded by the Canadian Institutes of Health Research Metabolomics Team Grant: MWG-146332.

Epigenetic modifications, particularly DNA methylation (DNAm) in cord blood, are an important biological marker of how external exposures during gestation can influence the in-utero environment and subsequent offspring development. Despite the recognized importance of DNAm during gestation, comparative studies to determine the consistency of these epigenetic signals across different ethnic groups are largely absent. To address this gap, we first performed epigenome-wide association studies (EWAS) of gestational age (GA) using newborn cord blood DNAm comparatively in a white European (n = 342) and a South Asian (n = 490) birth cohort living in Canada. Then, we capitalized on established cord blood epigenetic GA clocks to examine the associations between maternal exposures, offspring characteristics and epigenetic GA, as well as GA acceleration, defined as the residual difference between epigenetic and chronological GA at birth. Individual EWASs confirmed 1,211 and 1,543 differentially methylated CpGs previously reported to be associated with GA, in white European and South Asian cohorts, respectively, with a similar distribution of effects. We confirmed that Bohlin’s cord blood GA clock was robustly correlated with GA in white Europeans (r = 0.71; p = 6.0 × 10–54) and South Asians (r = 0.66; p = 6.9 × 10–64). In both cohorts, Bohlin’s clock was positively associated with newborn weight and length and negatively associated with parity, newborn female sex, and gestational diabetes. Exclusive to South Asians, the GA clock was positively associated with the newborn ponderal index, while pre-pregnancy weight and gestational weight gain were strongly predictive of increased epigenetic GA in white Europeans. Important predictors of GA acceleration included gestational diabetes mellitus, newborn sex, and parity in both cohorts. These results demonstrate the consistent DNAm signatures of GA and the utility of Bohlin’s GA clock across the two populations. Although the overall pattern of DNAm is similar, its connections with the mother's environment and the baby's anthropometrics can differ between the two groups. Further research is needed to understand these unique relationships.

Background: Mitochondrial DNA copy number (mtDNA-CN) is an accessible blood-based measurement believed to capture underlying mitochondrial (MT) function. The specific biological processes underpinning its regulation, and whether those processes are causative for disease, is an area of active investigation. Methods: We developed a novel method for array-based mtDNA-CN estimation suitable for biobank-scale studies, called ‘automatic mitochondrial copy (AutoMitoC).’ We applied AutoMitoC to 395,781 UKBiobank study participants and performed genome- and exome-wide association studies, identifying novel common and rare genetic determinants. Finally, we performed two-sample Mendelian randomization to assess whether genetically low mtDNA-CN influenced select MT phenotypes. Results: Overall, genetic analyses identified 71 loci for mtDNA-CN, which implicated several genes involved in rare mtDNA depletion disorders, deoxynucleoside triphosphate (dNTP) metabolism, and the MT central dogma. Rare variant analysis identified SAMHD1 mutation carriers as having higher mtDNA-CN (beta = 0.23 SDs; 95% CI, 0.18–0.29; p=2.6 × 10-19), a potential therapeutic target for patients with mtDNA depletion disorders, but at increased risk of breast cancer (OR = 1.91; 95% CI, 1.52–2.40; p=2.7 × 10-8). Finally, Mendelian randomization analyses suggest a causal effect of low mtDNA-CN on dementia risk (OR = 1.94 per 1 SD decrease in mtDNA-CN; 95% CI, 1.55–2.32; p=7.5 × 10-4). Conclusions: Altogether, our genetic findings indicate that mtDNA-CN is a complex biomarker reflecting specific MT processes related to mtDNA regulation, and that these processes are causally related to human diseases. Funding: No funds supported this specific investigation. Awards and positions supporting authors include: Canadian Institutes of Health Research (CIHR) Frederick Banting and Charles Best Canada Graduate Scholarships Doctoral Award (MC, PM); CIHR Post-Doctoral Fellowship Award (RM); Wellcome Trust Grant number: 099313/B/12/A; Crasnow Travel Scholarship; Bongani Mayosi UCT-PHRI Scholarship 2019/2020 (TM); Wellcome Trust Health Research Board Irish Clinical Academic Training (ICAT) Programme Grant Number: 203930/B/16/Z (CJ); European Research Council COSIP Grant Number: 640580 (MO); E.J. Moran Campbell Internal Career Research Award (MP); CISCO Professorship in Integrated Health Systems and Canada Research Chair in Genetic and Molecular Epidemiology (GP)

Background and Objectives Low buffy coat mitochondrial DNA copy number (mtDNA-CN) is associated with incident risk of stroke and poststroke mortality; however, its prognostic utility has not been extensively explored. Our goal was to investigate whether low buffy coat mtDNA-CN is a marker and causal determinant of poststroke outcomes using epidemiologic and genetic studies. Methods First, we performed association testing between baseline buffy coat mtDNA-CN measurements and 1-month poststroke outcomes in 3,498 cases of acute, first stroke from 25 countries from the international, multicenter case-control study Importance of Conventional and Emerging Risk Factors of Stroke in Different Regions and Ethnic Groups of the World (INTERSTROKE). Then, we performed 2-sample mendelian randomization analyses to evaluate potential causative effects of low mtDNA-CN on 3-month modified Rankin Scale (mRS) score. Genetic variants associated with mtDNA-CN levels were derived from the UK Biobank study (N = 383,476), and corresponding effects on 3-month mRS score were ascertained from the Genetics of Ischemic Stroke Functional Outcome (GISCOME; N = 6,021) study. Results A 1-SD lower mtDNA-CN at baseline was associated with stroke severity (baseline mRS score: odds ratio [OR] 1.27, 95% confidence interval [CI] 1.19–1.36; p = 4.7 × 10−12). Independently of baseline stroke severity, lower mtDNA-CN was associated with increased odds of greater 1-month disability (ordinal mRS score: OR 1.16, 95% CI 1.08–1.24; p = 4.4 × 10−5), poor functional outcome status (mRS score 3–6 vs 0–2: OR 1.21, 95% CI 1.08–1.34; p = 6.9 × 10−4), and mortality (OR 1.35, 95% CI 1.14–1.59; p = 3.9 × 10−4). Subgroup analyses demonstrated consistent effects across stroke type, sex, age, country income level, and education level. In addition, mtDNA-CN significantly improved reclassification of poor functional outcome status (net reclassification index [NRI] score 0.16, 95% CI 0.08–0.23; p = 3.6 × 10−5) and mortality (NRI score 0.31, 95% CI 0.19–0.43; p = 1.7 × 10−7) beyond known prognosticators. With the use of independent datasets, mendelian randomization revealed that a 1-SD decrease in genetically determined mtDNA-CN was associated with increased odds of greater 3-month disability quantified by ordinal mRS score (OR 2.35, 95% CI 1.13–4.90; p = 0.02) and poor functional outcome status (OR 2.68, 95% CI 1.05–6.86; p = 0.04). Discussion Buffy coat mtDNA-CN is a novel and robust marker of poststroke prognosis that may also be a causal determinant of poststroke outcomes. Classification of Evidence This study provides Class II evidence that low buffy coat mtDNA-CN (>1 SD) was associated with worse baseline severity and 1-month outcomes in patients with ischemic or hemorrhagic stroke.

BACKGROUND AND OBJECTIVES:Low buffy coat mitochondrial DNA copy number (mtDNA-CN) is associated with incident risk of stroke and poststroke mortality; however, its prognostic utility has not been extensively explored. Our goal was to investigate whether low buffy coat mtDNA-CN is a marker and causal determinant of poststroke outcomes using epidemiologic and genetic studies. METHODS:First, we performed association testing between baseline buffy coat mtDNA-CN measurements and 1-month poststroke outcomes in 3,498 cases of acute, first stroke from 25 countries from the international, multicenter case-control study Importance of Conventional and Emerging Risk Factors of Stroke in Different Regions and Ethnic Groups of the World (INTERSTROKE). Then, we performed 2-sample mendelian randomization analyses to evaluate potential causative effects of low mtDNA-CN on 3-month modified Rankin Scale (mRS) score. Genetic variants associated with mtDNA-CN levels were derived from the UK Biobank study (N = 383,476), and corresponding effects on 3-month mRS score were ascertained from the Genetics of Ischemic Stroke Functional Outcome (GISCOME; N = 6,021) study. RESULTS:A 1-SD lower mtDNA-CN at baseline was associated with stroke severity (baseline mRS score: odds ratio [OR] 1.27, 95% confidence interval [CI] 1.19-1.36; p = 4.7 × 10-12). Independently of baseline stroke severity, lower mtDNA-CN was associated with increased odds of greater 1-month disability (ordinal mRS score: OR 1.16, 95% CI 1.08-1.24; p = 4.4 × 10-5), poor functional outcome status (mRS score 3-6 vs 0-2: OR 1.21, 95% CI 1.08-1.34; p = 6.9 × 10-4), and mortality (OR 1.35, 95% CI 1.14-1.59; p = 3.9 × 10-4). Subgroup analyses demonstrated consistent effects across stroke type, sex, age, country income level, and education level. In addition, mtDNA-CN significantly improved reclassification of poor functional outcome status (net reclassification index [NRI] score 0.16, 95% CI 0.08-0.23; p = 3.6 × 10-5) and mortality (NRI score 0.31, 95% CI 0.19-0.43; p = 1.7 × 10-7) beyond known prognosticators. With the use of independent datasets, mendelian randomization revealed that a 1-SD decrease in genetically determined mtDNA-CN was associated with increased odds of greater 3-month disability quantified by ordinal mRS score (OR 2.35, 95% CI 1.13-4.90; p = 0.02) and poor functional outcome status (OR 2.68, 95% CI 1.05-6.86; p = 0.04). DISCUSSION:Buffy coat mtDNA-CN is a novel and robust marker of poststroke prognosis that may also be a causal determinant of poststroke outcomes. CLASSIFICATION OF EVIDENCE:This study provides Class II evidence that low buffy coat mtDNA-CN (>1 SD) was associated with worse baseline severity and 1-month outcomes in patients with ischemic or hemorrhagic stroke.

AbstractMitochondrial DNA copy number (mtDNA-CN) is an accessible blood-based measurement believed to capture underlying mitochondrial function. The specific biological processes underpinning its regulation, and whether those processes are causative for disease, is an area of active investigation. We developed a novel method for array-based mtDNA-CN estimation suitable for biobank-scale studies, called “AutoMitoC”. We applied AutoMitoC to 395,781 UKBiobank study participants and performed genome and exome-wide association studies, identifying novel common and rare genetic determinants. Overall, genetic analyses identified 73 loci for mtDNA-CN, which implicated several genes involved in rare mtDNA depletion disorders, dNTP metabolism, and the mitochondrial central dogma. Rare variant analysis identified SAMHD1 mutation carriers as having higher mtDNA-CN (beta=0.23 SDs; 95% CI, 0.18-0.29; P=2.6×10−19), a potential therapeutic target for patients with mtDNA depletion disorders, but at increased risk of breast cancer (OR=1.91; 95% CI, 1.52-2.40; P=2.7×10−8). Finally, Mendelian randomization analyses suggest a causal effect of low mtDNA-CN on dementia risk (OR=1.94 per 1 SD decrease in mtDNA-CN; 95% CI, 1.55-2.32; P=7.5×10−4). Altogether, our genetic findings indicate that mtDNA-CN is a complex biomarker reflecting specific mitochondrial processes related to mtDNA regulation, and that these processes are causally related to human diseases.