DNA methylation plays a fundamental role in the regulation of the genome, but the optimal strategy for analysis of genome-wide DNA methylation data remains to be determined. We developed a ...comprehensive analysis pipeline for epigenome-wide association studies (EWAS) using the Illumina Infinium HumanMethylation450 BeadChip, based on 2,687 individuals, with 36 samples measured in duplicate. We propose new approaches to quality control, data normalisation and batch correction through control-probe adjustment and establish a null hypothesis for EWAS using permutation testing. Our analysis pipeline outperforms existing approaches, enabling accurate identification of methylation quantitative trait loci for hypothesis driven follow-up experiments.
Abstract
Background
Disparities in COVID-19 outcomes exist on the basis of ethnicity and comorbidities. Minority ethnic groups in the UK are known to have poorer COVID-19 outcomes, but also an ...increased prevelance of certain comorbidities associated with severe outcomes. Additionally, despite the prevalence of certain psychiatric disorders there is a lack of research establishing their relationship with COVID-19 outcomes.
Methods
We used UK Biobank data, involving 472,182 participants, to test for an association between comorbidities and COVID-19 diagnosis (
n
= 30,901); and to test for an association between comorbidities and severe COVID-19 (
n
= 3182). This was done by performing univariable and multivariable logistic regression analysis, estimating odds ratios (ORs) and their 95% confidence intervals (95% CIs). The comorbidities studied were coronary heart disease (CHD), hypertension, type II diabetes mellitus (T2DM), obesity, chronic kidney disease (CKD), depression and anxiety. Multivariable models were adjusted for various socioeconomic, demographic and health-related confounders. We then performed sub-group analysis by common UK ethnic groups (White, South Asian, and Black).
Results
Increased prevalence of all studied comorbidities was seen in both outcomes, compared to the rest of the cohort. All studied comorbidities were associated with an increased risk of COVID-19 infection and severity across all models. For example, the adjusted ORs (95% CI) for depression were 1.112 (1.083 – 1.161) for COVID-19 diagnosis and 2.398 (2.163 – 2.658) for severe COVID-19. Sub-group analysis revealed stronger associations of COVID-19 diagnosis and severe COVID-19 for South-Asian participants for CHD (OR 1.585 95% CI 1.194–2.105 for COVID-19 diagnosis and 3.021 1.683–5.390 for severe COVID-19), hypertension (1.488 1.231–1.799; 3.399 1.862–6.206) and T2DM (1.671 1.346–2.076; 5.412 3.130–9.357) compared to White participants (1.264 1.195–1.336 and 1.627 1.441–1.837 for CHD; 1.131 1.097–1.116 and 2.075 1.885–2.284 for hypertension; 1.402 1.331–1.476 and 2.890 2.596–3.216 for T2DM). Similar results were seen for Black participants with CKD and hypertension.
Conclusion
Specific comorbidities are risk factors for poorer COVID-19 outcomes, supporting targeted interventions and policy aimed at individuals with these comorbidities. Although further research is required, there’s also a need for targeted policies for ethnic minorities assessing the unique reasons they are at greater risk of poor COVID-19 outcomes.
Insulin resistance and related metabolic disturbances are major risk factors for the higher T2D risk and associated morbidity and mortality amongst South Asians. The contribution of physical activity ...to the increased prevalence of insulin resistance and related disturbances amongst South Asians is unknown.
We recruited 902 South Asian and European men and women, aged 35-85 years from the ongoing LOLIPOP study. Clinical characterisation comprised standardised questionnaire and measurement of height, weight, waist and hip circumference and blood pressure. Fasting bloods were taken for assessment of glucose, insulin, lipids and HbA1c. Physical activity was quantified using a validated accelerometer, Actigraph GT3X+, worn for 7 days. Univariate and multivariate approaches were used to investigate the relationship between ethnicity, physical activity, insulin resistance and related metabolic disturbances.
Total physical activity was ~31% (P = 0.01) lower amongst South Asians compared to Europeans (Mean MET.minutes SD: 1505.2 52 vs. 2050.9 86.6, P<0.001). After adjusting for age and sex, total physical activity had a negative association with HOMA-IR (B SE: -0.18 0.08, P = 0.04) and fasting glucose levels (BSE: -0.11 0.04, P = 0.02). There was no association between physical activity and other glycemic and lipid parameters. Total physical activity per week contributed towards the differences in insulin resistance and associated metabolic disturbances between South Asians and Europeans.
Lower levels of physical activity may contribute to the increased insulin resistance in South Asians compared to Europeans. Our results suggest that lifestyle modification through increased physical activity may help to improve glucose metabolism and reduce the burden of excess T2D and related complications amongst South Asians.
COVID-19 manifests with huge heterogeneity in susceptibility and severity outcomes. UK Black Asian and Minority Ethnic (BAME) groups have demonstrated disproportionate burdens. Some variability ...remains unexplained, suggesting potential genetic contribution. Polygenic Risk Scores (PRS) can determine genetic predisposition to disease based on Single Nucleotide Polymorphisms (SNPs) within the genome. COVID-19 PRS analyses within non-European samples are extremely limited. We applied a multi-ethnic PRS to a UK-based cohort to understand genetic contribution to COVID-19 variability.
We constructed two PRS for susceptibility and severity outcomes based on leading risk-variants from the COVID-19 Host Genetics Initiative. Scores were applied to 447,382 participants from the UK-Biobank. Associations with COVID-19 outcomes were assessed using binary logistic regression and discriminative power was validated using incremental area under receiver operating curve (ΔAUC). Variance explained was compared between ethnic groups via incremental pseudo-R
(ΔR
).
Compared to those at low genetic risk, those at high risk had a significantly greater risk of severe COVID-19 for White (odds ratio OR 1.57, 95% confidence interval CI 1.42-1.74), Asian (OR 2.88, 95% CI 1.63-5.09) and Black (OR 1.98, 95% CI 1.11-3.53) ethnic groups. Severity PRS performed best within Asian (ΔAUC 0.9%, ΔR
0.98%) and Black (ΔAUC 0.6%, ΔR
0.61%) cohorts. For susceptibility, higher genetic risk was significantly associated with COVID-19 infection risk for the White cohort (OR 1.31, 95% CI 1.26-1.36), but not for Black or Asian groups.
Significant associations between PRS and COVID-19 outcomes were elicited, establishing a genetic basis for variability in COVID-19. PRS showed utility in identifying high-risk individuals. The multi-ethnic approach allowed applicability of PRS to diverse populations, with the severity model performing well within Black and Asian cohorts. Further studies with larger sample sizes of non-White samples are required to increase statistical power and better assess impacts within BAME populations.
The mechanisms that predispose to hypertension, coronary artery disease (CAD), and type 2 diabetes (T2D) in individuals of normal weight are poorly understood. In contrast, in monogenic primary ...lipodystrophy-a reduction in subcutaneous adipose tissue-it is clear that it is adipose dysfunction that causes severe insulin resistance (IR), hypertension, CAD, and T2D. We aimed to test the hypothesis that common alleles associated with IR also influence the wider clinical and biochemical profile of monogenic IR. We selected 19 common genetic variants associated with fasting insulin-based measures of IR. We used hierarchical clustering and results from genome-wide association studies of eight nondisease outcomes of monogenic IR to group these variants. We analyzed genetic risk scores against disease outcomes, including 12,171 T2D cases, 40,365 CAD cases, and 69,828 individuals with blood pressure measurements. Hierarchical clustering identified 11 variants associated with a metabolic profile consistent with a common, subtle form of lipodystrophy. A genetic risk score consisting of these 11 IR risk alleles was associated with higher triglycerides (β = 0.018; P = 4 × 10(-29)), lower HDL cholesterol (β = -0.020; P = 7 × 10(-37)), greater hepatic steatosis (β = 0.021; P = 3 × 10(-4)), higher alanine transaminase (β = 0.002; P = 3 × 10(-5)), lower sex-hormone-binding globulin (β = -0.010; P = 9 × 10(-13)), and lower adiponectin (β = -0.015; P = 2 × 10(-26)). The same risk alleles were associated with lower BMI (per-allele β = -0.008; P = 7 × 10(-8)) and increased visceral-to-subcutaneous adipose tissue ratio (β = -0.015; P = 6 × 10(-7)). Individuals carrying ≥17 fasting insulin-raising alleles (5.5% population) were slimmer (0.30 kg/m(2)) but at increased risk of T2D (odds ratio OR 1.46; per-allele P = 5 × 10(-13)), CAD (OR 1.12; per-allele P = 1 × 10(-5)), and increased blood pressure (systolic and diastolic blood pressure of 1.21 mmHg per-allele P = 2 × 10(-5) and 0.67 mmHg per-allele P = 2 × 10(-4), respectively) compared with individuals carrying ≤9 risk alleles (5.5% population). Our results provide genetic evidence for a link between the three diseases of the "metabolic syndrome" and point to reduced subcutaneous adiposity as a central mechanism.
Genome-wide association studies (GWAS) in coronary artery disease (CAD) had identified 66 loci at 'genome-wide significance' (P < 5 × 10
) at the time of this analysis, but a much larger number of ...putative loci at a false discovery rate (FDR) of 5% (refs. 1,2,3,4). Here we leverage an interim release of UK Biobank (UKBB) data to evaluate the validity of the FDR approach. We tested a CAD phenotype inclusive of angina (SOFT; n
= 10,801) as well as a stricter definition without angina (HARD; n
= 6,482) and selected cases with the former phenotype to conduct a meta-analysis using the two most recent CAD GWAS. This approach identified 13 new loci at genome-wide significance, 12 of which were on our previous list of loci meeting the 5% FDR threshold, thus providing strong support that the remaining loci identified by FDR represent genuine signals. The 304 independent variants associated at 5% FDR in this study explain 21.2% of CAD heritability and identify 243 loci that implicate pathways in blood vessel morphogenesis as well as lipid metabolism, nitric oxide signaling and inflammation.
Objectives The purpose of this study was to define age-, sex-, and ethnicity-specific reference values for 3-dimensional echocardiographic (3DE) left ventricular (LV) volumes and LV ejection fraction ...(LVEF) in a large cohort of European white and Indian Asian subjects. Background Transthoracic 3DE imaging is recommended for the routine evaluation of LV volumes and function. However, there remains a lack of population-based reference values for 3DE LV volumes and LVEF, hindering adoption of this technique into routine clinical practice. Methods We identified subjects from the LOLIPOP (London Life Sciences Prospective Population) study who were free of clinical cardiovascular disease, hypertension, and type 2 diabetes. All subjects underwent transthoracic 2-dimensional and 3D echocardiography for quantification of LV end-systolic volume index, LV end-diastolic volume index, and LVEF. Results 3DE image quality was satisfactory in 978 subjects (89%) for the purposes of LV volumetric analysis. Indexed 3DE LV volumes were significantly smaller in female compared with male subjects and in Indian Asians compared with European whites. Upper limit of normal (mean ± 2 SD) reference values for the LV end-systolic volume index and LV end-diastolic volume index for the 4 ethnicity-sex subgroups were, respectively, as follows: European white men, 29 ml/m2 and 67 ml/m2 ; Indian Asian men, 26 ml/m2 and 59 ml/m2 ; European white women, 24 ml/m2 and 58 ml/m2 ; Indian Asian women, 23 ml/m2 and 55 ml/m2 , respectively. Compared with 3DE studies, 2-dimensional echocardiography underestimated the LV end-systolic volume index and LV end-diastolic volume index by an average of 2.0 ml/m2 and 4.7 ml/m2 , respectively. LVEF was similar between in all 4 groups and between 2- and 3-dimensional techniques, with a lower cutoff of 52% for the whole cohort. Conclusions These reference values are based on the largest 3DE study performed to date that should facilitate the standardization of the technique and encourage its adoption for the routine assessment of LV volumes and LVEF in the clinical echocardiography laboratory. This study supports the application of ethnicity-specific reference values for indexed LV volumes.
We carried out a genome-wide association study (318,237 SNPs) for insulin resistance and related phenotypes in 2,684 Indian Asians, with further testing in 11,955 individuals of Indian Asian or ...European ancestry. We found associations of rs12970134 near MC4R with waist circumference (P = 1.7 × 10−9) and, independently, with insulin resistance. Homozygotes for the risk allele of rs12970134 have ∼2 cm increased waist circumference. Common genetic variation near MC4R is associated with risk of adiposity and insulin resistance.
We determined the relationships between DNA sequence variation and DNA methylation using blood samples from 3,799 Europeans and 3,195 South Asians. We identify 11,165,559 SNP-CpG associations ...(methylation quantitative trait loci (meQTL), P < 10
), including 467,915 meQTL that operate in trans. The meQTL are enriched for functionally relevant characteristics, including shared chromatin state, High-throuhgput chromosome conformation interaction, and association with gene expression, metabolic variation and clinical traits. We use molecular interaction and colocalization analyses to identify multiple nuclear regulatory pathways linking meQTL loci to phenotypic variation, including UBASH3B (body mass index), NFKBIE (rheumatoid arthritis), MGA (blood pressure) and COMMD7 (white cell counts). For rs6511961 , chromatin immunoprecipitation followed by sequencing (ChIP-seq) validates zinc finger protein (ZNF)333 as the likely trans acting effector protein. Finally, we used interaction analyses to identify population- and lineage-specific meQTL, including rs174548 in FADS1, with the strongest effect in CD8
T cells, thus linking fatty acid metabolism with immune dysregulation and asthma. Our study advances understanding of the potential pathways linking genetic variation to human phenotype.
While increased age is an established risk factor for COVID-19, there is great heterogeneity in outcomes within age groups. This is because chronological age does not reflect health, unlike ...biological age. We intend to investigate the association between accelerated ageing and COVID-19 outcomes through the lens of three measures, namely phenotypic age acceleration (PhenoAgeAccel), telomere length (Adjusted T/S Ratio) and facial ageing, and to examine whether there are differences across ethnic groups.
Taking participants from the UK Biobank, we associated accelerated ageing with severe COVID-19 outcomes, defined as COVID-related hospitalisation or death. Separate logistic regressions models were created for age and the three accelerated ageing-related variables, adjusting for a variety of covariates in each model. Multivariable logistic regression models were also created within White, Black, Asian and Other ethnic groups to assess for potential differing associations. Forward likelihood ratio logistic regression models were created to evaluate importance of the variables and to assess for patterns of association across the total population and ethnic groups.
After adjusting for all covariates, the odds ratio (OR) and 95% confidence interval (95% CI) of COVID-19 severe outcomes for age was 1.080 (1.074-1.086). After further adjusting age for the accelerated ageing variables, the ORs were 1.029 (1.020-1.039) for PhenoAgeAccel and 0.847 (0.772-0.929) for Facial Ageing's "Younger Than You Are" while Adjusted T/S ratio and "Older Than You Are" were statistically insignificant. The OR for age remained similar across ethnic groups. Both PhenoAgeAccel and younger facial ages in the White population and PhenoAgeAccel in the Black population had ORs of 1.031 (1.021-1.042), 0.853 (0.774-0.939), and 1.049 (1.001-1.100), respectively. Both Adjusted T/S Ratio and older facial ages showed statistical insignificance in all ethnicities. In forward logistic regression, age and PhenoAgeAccel were the age-related variables selected most frequently in all models.
Accelerated ageing is associated with increased COVID-19 severity. The mechanisms at work here are likely immunosenescence and inflamaging. This association indicates that anti-ageing treatment may improve COVID-19 outcome. The results within ethnic groups and that of telomere length were inconclusive, but point to a need for future, more focused research on the topic.