http://orcid.org/0000-0002-3332-8454
Figures
Abstract
Introduction
Variability in red blood cell volumes (distribution width, RDW) increases with age and is strongly predictive of mortality, incident coronary heart disease and cancer. We investigated inherited genetic variation associated with RDW in 116,666 UK Biobank human volunteers.
Results
A large proportion RDW is explained by genetic variants (29%), especially in the older group (60+ year olds, 33.8%, <50 year olds, 28.4%). RDW was associated with 194 independent genetic signals; 71 are known for conditions including autoimmune disease, certain cancers, BMI, Alzheimer’s disease, longevity, age at menopause, bone density, myositis, Parkinson’s disease, and age-related macular degeneration. Exclusion of anemic participants did not affect the overall findings. Pathways analysis showed enrichment for telomere maintenance, ribosomal RNA, and apoptosis. The majority of RDW-associated signals were intronic (119 of 194), including SNP rs6602909 located in an intron of oncogene GAS6, an eQTL in whole blood.
Conclusions
Although increased RDW is predictive of cardiovascular outcomes, this was not explained by known CVD or related lipid genetic risks, and a RDW genetic score was not predictive of incident disease. The predictive value of RDW for a range of negative health outcomes may in part be due to variants influencing fundamental pathways of aging.
Citation: Pilling LC, Atkins JL, Duff MO, Beaumont RN, Jones SE, Tyrrell J, et al. (2017) Red blood cell distribution width: Genetic evidence for aging pathways in 116,666 volunteers. PLoS ONE 12(9): e0185083. https://doi.org/10.1371/journal.pone.0185083
Editor: Dan E. Arking, Johns Hopkins University, UNITED STATES
Received: May 22, 2017; Accepted: September 6, 2017; Published: September 28, 2017
Copyright: © 2017 Pilling et al. This is an open access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.
Data Availability: The data is available from figshare at https://doi.org/10.6084/m9.figshare.5395504.v1.
Funding: This work was supported by an award to DM, TF, AM and LH by the UK Medical Research Council (grant number MR/M023095/1). SEJ is funded by the Medical Research Council (grant: MR/M005070/1). JT is funded by a Diabetes Research and Wellness Foundation Fellowship. RB is funded by the Wellcome Trust and Royal Society grant: 104150/Z/14/Z. MAT, MNW and AM are supported by the Wellcome Trust Institutional Strategic Support Award (WT097835MF). ARW, HY, and TF are supported by the European Research Council grant: 323195:GLUCOSEGENES-FP7-IDEAS-ERC. LF is supported by the Intramural Research Program of the National Institute on Aging, U.S. National Institutes of Health. Input from MD, CLK and GK was supported by the University of Connecticut Health Center. This research has been conducted using the UK Biobank Resource under Application Number 14631. The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript.
Competing interests: The authors have declared that no competing interests exist.
Introduction
Increased variation in a person’s Red Blood Cell (RBC) volumes (RBC distribution width (RDW), also termed anisocytosis) is strongly predictive of a range of incident cardiovascular conditions, cancers and mortality [1–3]. Although RDW is routinely measured in clinical hematology reporting–it is calculated by dividing the standard deviation of mean cell volume (MCV) by the MCV and multiplying by 100, to yield a RDW percentage [4]–it is only used clinically for diagnosis of anemia subtypes. Understanding the mechanisms involved in the links between increased RDW and negative health outcomes could provide clues to potential interventions to improve prognosis in those with high RDW who are not anemic, particularly in older people.
Established clinical causes of increased RDW include anemia and other iron or folate deficiencies [5], dyslipidemia [6] and other metabolic abnormalities, and inflammation [7]. Proposed mechanisms for increased RDW also include impaired erythropoiesis (the generation of new RBC) perhaps due to effects of inflammation or senescence of erythropoietic cells in the bone marrow, plus variations in RBC survival [8]. A previous analysis of 36 blood cell traits identified genetic variants [9], but the genetic signals for RDW were not investigated in depth in relation to the mechanisms that might explain the predictive value of RDW for negative health outcomes in people with normal hemoglobin levels.
We aimed to investigate RDW (overall and excluding anemia) using genetic analysis in a large population cohort, to identify underlying mechanisms. This involved genome-wide analysis of associations to find independent signals, investigations of biological pathways implicated by the results, and overlap with known risk alleles. We also examined associations between RDW and known variant genetic risk score analysis for conditions predicted by RDW, including cardiovascular disease. For this analysis, we used the exceptionally large UK Biobank volunteer sample with standardized measures of RDW across the cohort.
Results
We included 116,666 UK Biobank participants of white/British descent with complete hematology measures, covariate data, and genetics data from the interim data release (May 2015) in our analyses. The mean age was 57 years (SD: 7.9, min = 40, max = 70), with 52,541 aged 60 to 70 years old, and the majority (52.6%) were female: Table 1.
Variance explained by genotypes
We estimated that 29.3% (SE = 0.5%) of the variance in RDW was accounted for by 457,643 directly genotyped variants with MAF>0.1%. In a secondary analysis we estimated the proportion explained in two groups: those aged ≥60 years (n = 52,541), and those aged <50 years (n = 24,988). The proportion of variance in RDW explained by the genetic variants was greater in the older group (33.8%, SE = 1.0%) compared to the younger group (28.4%, SE = 2.0%), and this difference is statistically significant (p = 0.012).
In sensitivity analyses including only those directly genotyped SNPs associated with RDW (at 2 separate p-value cut-offs: p<5x10-5 and p<5x10-8; n SNPs = 4,110 and 2,095 respectively) we observed a reduction in the variance in RDW explained, as expected, but the conclusions remained consistent: in all participants the variance in RDW explained by the included SNPs was 18% (SE = 0.6%) and 14.2% (SE = 0.7%) for the 2 p-value thresholds respectively. There was still a marked difference in the variance explained between the older group (19.6%, SE = 0.8%; and 15.4%, SE = 0.8%) and the younger group (17.3%, SE = 0.9%; and 13.5%, SE = 1%).
We estimated the proportion of variance in CHD (10,280 cases in 116,666 participants) accounted for by the variants to be 5.95% (SE = 0.45%). The proportion of the variance shared between RDW and CHD attributable to genetics is 6.62% (SE = 2.69%: 95% CIs = 1.35 to 11.9%).
Genome-wide association study
Of the 16,832,071 genetic variants included in this GWAS, 30,988 were significantly (p<5x10-8) associated with RDW (Fig 1; full results available to download here: https://doi.org/10.6084/m9.figshare.5395504.v1) after adjustment for age, sex, assessment center and array type (genetic relatedness is accounted for in the linear mixed models approach so no PCs are included–see methods). A quantile-quantile plot shows clear inflation (S1 Fig; Lambda = 1.15), consistent with polygenic inheritance of many causal variants of small effects which can now be discovered using large sample sizes [10]. The 30,988 variants were mapped to 141 loci (runs of variants separated by <2Mb) on the genome, and included 194 independent signals after conditional analysis (“conditional SNPS”) (S1 Table).
The variants are grouped into 194 independent signals, colored blue if a variant in the signal is associated with any trait in the NHGRI-EBI GWAS catalog of known associations, otherwise colored in red. The y-axis (–log10 p-values) is limited to 30 for clarity, as the max value is 200. See S1 Table for RDW associations for each signal, and S4 Table for mapping to the catalog. Horizontal line p = 5x10-8.
After excluding 17,795 participants with anemia the GWAS of RDW was repeated on the 98,871 remaining participants; all 194 conditional SNPs remained nominally associated with RDW (p<0.001), but 24 were no longer genome-wide significant (p<5x10-8), possibly due to reduced power in the smaller sample size. All subsequent analyses are based on the results from the analysis of all participants.
Functional implications of RDW-associated genetic variants
We utilized the UCSC Variant Annotation Integrator (http://www.genome.ucsc.edu/cgi-bin/hgVai) and the Ensembl Variant Effect Predictor (http://grch37.ensembl.org/Homo_sapiens/Tools/VEP) to interrogate a number of genomic annotation databases for the conditional SNPs. The majority (119 out of 194 total) were intronic, and 12 were located in the 3’ or 5’ un-translated regions (see S2 Table for complete variant annotation output). We found that 37 of 194 independent RDW-associated signals are known eQTLs (i.e. affect the expression of a gene) in whole blood (http://genenetwork.nl/bloodeqtlbrowser): 15 of these affect the gene predicted by the variant annotation integrator, including rs6602909, located in an intron of oncogene GAS6 [11] (see S3 Table for complete matching of RDW-associated signals to eQTL data). SNP rs7775698 is an trans-eQTL for more than 30 genes.
Fifteen of the RDW-associated signals were exonic: 11 missense, 3 synonymous, and a 17 base-pair exonic deletion in gene SMIM1. The missense and deletion variants are shown in Table 2. PolyPhen-2 predicted that three missense variants are “probably damaging” to the protein function of genes TRIM58, PLD1 and PNPLA3. Variant rs2075995 was predicted to be “possibly damaging” to the protein function of gene E2F2. TRIM58 is a ubiquitin ligase induced during late erythropoiesis [11]; PLD1 is a phospholipase implicated in processes including membrane trafficking [11]; PNPLA3 is a triacylglycerol lipase in adipocytes and variant rs738409 is associated with susceptibility to Non-alcoholic Fatty Liver Disease [12]. Gene SMIM1 is involved in RBC formation [11] and a 17 base-pair deletion (rs566629828) causing an exonic frameshift is strongly associated with increased RDW.
GWAS catalogue of known genetic associations
Of the 194 conditional SNPs associated with RDW, 77 mapped to at least one trait in the catalogue of published GWAS (downloaded 13th March 2017). This was arrived at by; filtering the 33,005 SNP-trait associations to those with p<5x10-8 (leaving 14,148 SNP-trait pairs for analysis); matching the positions to the UK Biobank results (13,146); filtering to those with significant RDW association in our analysis (p<5x10-8), leaving 923 SNP-trait pairs (420 unique SNPs; some SNPs are associated with multiple traits). These 420 unique SNPs mapped to 77 of the 194 conditional SNPs associated with RDW. These are shown in Fig 1; see S4 Table for further detail.
Traits also associated with the individual RDW variants included iron metabolism and several other red cell measures. Variants present were also associated with BMI, several lipids, hemoglobin A1C and metabolic syndrome, as well as height. Autoimmune associated conditions included autoimmune thyroid disease, type 1 diabetes, Crohn’s disease, inflammatory bowel disease, rheumatoid arthritis, systemic lupus erythematosus and ulcerative colitis. Variants linked to Lung, ovary and nasopharyngeal cancers were present. In addition, conditions associated with aging were represented, including variants linked to Alzheimer’s disease, age at menopause, bone density, Myositis, Parkinson’s disease, macular degeneration, C-reactive protein levels and longevity. For Alzheimer’s and longevity these were known SNPs in the APOE gene region.
Gene ontology enrichment
MAGENTA software [13] identified pathways enriched in the genes mapped to variants significantly associated with RDW, including telomere maintenance, ribosomal RNA transcription and histone modifications (Table 3), plus apoptosis. In addition, pathways related to lipid metabolism (in particular chylomicrons) were also enriched. Additional information, such as the expected/observed number of associated genes for each pathway is included in S5 Table.
Genetic risk score associations
We tested 20 Genetic Risk Scores (GRS) for potentially explanatory traits for the predictive value of RDW for negative health outcomes; 7 were significant after adjustment for multiple testing (p<0.0025). Three were associated with raised RDW (HDL, type-1 diabetes, and BMI); four were associated with lower RDW (triglycerides, LDL, systolic blood pressure, and Alzheimer’s disease) (Fig 2; S6 Table). After exclusion of the ApoE locus the HDL GRS remained significant (p = 0.003), but the AD GRS was no longer significant (p = 0.84).
* FDR = false-discovery rate adjusted significant association. Genetic Risk Scores (GRS) were z-transformed prior to analysis. Linear regression model against RDW (z-transformed) including 116,666 participants, adjusted for age, sex, assessment center and population structure (genetic PCs 1–5). LDL (low-density lipoprotein), HDL (high-density lipoprotein), TG (triglycerides), SBP (systolic blood pressure), CAD (coronary artery disease), T1D (type-1 diabetes), T2D (type-2 diabetes), AMD (age-related macular degeneration), AD (Alzheimer’s disease), FVC (forced vital capacity), BMI (body mass index), IBD (inflammatory bowel disease). Full results in S6 Table.
Genetic risks for Crohn’s disease and inflammatory bowel disease were nominally associated with increased RDW (Beta = 0.009: 95% CIs = 0.003 to 0.015; Beta = 0.009: 95% CIs = 0.003 to 0.014; respectively). Risk scores for coronary artery disease were not associated with RDW, including after removal of the lipid related variants (Beta = 0.00: 95% CIs = -0.006 to 0.005). A genetic risk score for telomere length was not significantly associated with RDW.
We also created a GRS for RDW using the 194 conditional SNPs identified in this analysis. Although the GRS was associated with anemia (OR = 1.10: 95% CIs = 1.08 to 1.12) there was no association with mortality, coronary artery disease, or cancer (p = 0.24, p = 0.24, and p = 0.44, respectively). In sensitivity analyses excluding RDW-SNPs associated with CAD p<5x10-8 in a recent meta-analysis [14], in case results were biased due to invalid genetic instruments, or using the MR-Egger method to check for pleiotropy [15], the conclusions remained consistent.
Discussion
Variation in RBC size (RBC Distribution Width, RDW) increases markedly with age [16] and high RDW values are strongly predictive of increased mortality, plus incident cardiovascular disease and certain cancers [3,17]. However, RDW is not generally considered as being a clinically useful measure outside the assessment of anemia sub-type, perhaps because the mechanisms explaining its prognostic value in people without anemia are unclear. In this study we investigated RDW using genetic analysis to understand the molecular mechanisms underpinning variation in RBC size.
A large proportion of RDW variation (29.3%) was attributable to common genetic variants in this analysis, and the variation explained by genetic variants appeared to increase with age, contrary to common assumptions that genetic effects decrease with advancing age. This could be due to effects accumulating over a time. Many of the “hallmarks of aging” also have this property [18]. Additionally, many of the RDW-associated genetic variants (in 71 of 194 conditionally independent signals) have previously been associated with other traits including metabolic syndrome, certain cancers, and autoimmune disease as well as aging related conditions including menopausal age.
Our analysis of genetic risk scores (GRS) showed that participants with genetically-increased risk for coronary artery disease or cancer did not have significantly higher RDW. Similarly, we found that a GRS for RDW (using the 194 variants) had was not associated with mortality, incident coronary heart disease (CHD), or cancer. Consistent with this, the proportion of the variance shared between RDW and CHD attributable to genetics was only 6.6%; therefore the majority of genetically influenced red cell variation is independent of CHD.
As higher RDW is associated with CVD epidemiologically, we might have hypothesized that genetic risks for adverse lipid levels or blood pressures affect RDW, but instead we found associations in the opposite directions: GRS analysis showed that participants with genetically lower LDL levels, triglyceride levels, or systolic blood pressure, had higher RDW, and genetically higher HDL was associated with greater RDW. Published observational associations between lipids and RDW are only partially consistent with our findings here: Lippi et al found that LDL and HDL cholesterol were both negatively related to RDW, and triglyceride levels were positively related [6]. It is known that RBC have a role in cholesterol homeostasis by transporting cholesterol in the plasma membranes, with significant inter-individual differences not entirely explained by age or cholesterol levels [19]. The relationship between lipids and RDW is complex and requires further investigation. We also observed that genetically increased risk of type-1 diabetes was associated with increased RDW–further evidence for autoimmune involvement, in addition to the overlap in significant SNPs in the GWAS catalogue–and that genetically increased BMI is associated with increased RDW.
Genetic variants associated with RDW were enriched in expected pathways, including iron homeostasis, but we also found evidence for telomere maintenance, ribosomal RNA production, and a number of nucleosome and histone pathways. Short telomere length is a hallmark of cellular aging to senescence [18] and is associated with many risk factors of disease, however the causal direction is still uncertain [20], and longer telomeres have been linked to risk of cancer [21]. Kozlitina et al reported in 2012 that increased RDW is associated with shorter telomeres in leukocytes [22]. We created a genetic risk score for telomere length but it was not associated with RDW (Beta = 0.004: 95% CIs = -0.002 to 0.010); further work is required to clarify this association.
We also found enrichment of RNA polymerase I (which transcribes ribosomal RNA) and RNA polymerase III (which transcribes transfer RNA); both are required for protein synthesis, including hemoglobin, and can even function as regulators of gene expression in their own right [23], suggesting these are key factors for consistent production of RBC. Deregulation of transcription and proteostasis are hallmarks of aging [18], and we have previously reported deregulation of gene expression of the transcriptional machinery with advancing age [24].
Four of the conditionally independent genetic variants associated with RDW are exonic and affect the amino-acid sequence of the protein products. Most others are intronic or intergenic, and may be regulatory; this is supported by published eQTL data [25], in which 37 of the RDW-associated signals have been reported to affect the expression of a gene in whole blood. Of the 194 independent variants identified, 119 are intronic; this does not appear to be an uncommonly high proportion, as 23 of 39 novel BMI variants identified by Locke et al. 2015 were intronic [26] and 6 of 12 novel blood pressure variants identified by Kato et al. 2015 were intronic [27]. Intronic variants may have direct effects on protein expression by affecting the splicing or processing of the pre-mRNA, a complex process to produce the correct mRNA and subsequent protein [28], or by altering expression of a regulatory non-coding RNA such as a micro RNA. These will be useful targets for future research to determine how these variants ultimately affect consistency in red blood cell size.
Limitations
The UK Biobank is a volunteer study which achieved only a 5% response rate, so at assessment the participants were healthier than the general population. However, there was substantial variation in RDW within the participants so the results can still be generalized to the wider population [29].
More work is needed to replicate these findings in independent cohorts, and to establish the effects in other populations, particularly non-European ancestries. No data have been released regarding the UK Biobank participant’s lipid and other blood assays; once this is available, further investigations into the complex relationship between RDW and lipids can be performed.
The GWAS catalogue does not contain every published GWAS, especially the most recent studies, but contains many of the largest meta-analyses for traits such as cardiovascular disease and cancer. It is therefore likely we have missed some studies using this method, therefore our results present an approximation of the overlap between RDW signals and other traits.
Conclusions
Variation in RDW has a substantial genetic component, and this increases with increasing age. Although increased RDW is predictive of cardiovascular outcomes, this was not explained by known CVD or related lipid genetic risks. The predictive value of RDW for a range of negative health outcomes may in part be due to variants influencing fundamental pathways of aging, but it may also be reflective of exposures or underlying conditions.
Materials and methods
The UK Biobank study recruited 503,325 volunteers aged 40–70 who were seen between 2006 and 2010. Data includes RBC distribution width (RDW) and other clinical hematology measures, extensive questionnaires including smoking behavior and education history, and follow-up using electronic medical records. Currently one third of the participants have available genotype information. We utilized data from 116,666 participants of white/British descent with all available data. UK Biobank has approval from the North West Multi-centre Research Ethics Committee (MREC), which covers the UK. UK Biobank data is available to all bona fide researchers, in the UK and internationally, without preferential access. Applications are reviewed, with ethical advice from from both the University of Oxford’s Ethox Centre and the Ethics and Governance Council, where appropriate. Following approval and signing of a materials transfer agreement, de-identified data is transferred to researchers [30]. Information and access to data can be found online (www.ukbiobank.ac.uk/scientists).
Phenotypes
RDW is a measure of the variability in the mean size of the RBC in each participant (in % units). It was measured using four Beckman Coulter LH750 instruments within 24 hours of blood draw, with extensive quality control performed by UK Biobank [31]. RDW is a continuous, highly skewed trait, therefore we used quantile normalization of the continuous measure to create a Gaussian distribution, so that the normality assumption of the linear regression models were not violated.
Anemia was determined both using self-reported diagnosis, electronic medical records (ICD10: D64* and D5* categories), or by low hemoglobin levels at the baseline assessment (<120g/L in females, <130g/L in males: from WHO definition [32]).
Coronary heart disease (CHD) was defined using self-reported myocardial infarction or angina, or diagnosis in the electronic medical records (ICD10: I20-I25).
Variance explained and genetic correlation
Variance-components analysis determines the heritable phenotypic variation in complex traits explained by genetic variants. We used BOLT-REML to determine the variance in RDW explained by the common, genotyped variants (n = 457,643 directly genotyped variants with minor allele frequency (MAF) >0.1%, HWE p>1x10-6 and missingness <1.5%) using restricted maximum likelihood estimation [33]. The method performs multiple iterations to approximate BOLT-REML was also used to estimate the genetic correlation between RDW and CHD.
To determine whether differences in heritability between the young/old groups was statistically significant we computed the t-statistic using Eq (1) from which a two-tailed p-value is derived.
Genome-wide association study
We used data from the May 2015 release of the UK Biobank genetics data, where 152,248 participants were genotyped. Two custom Affymetrix genotyping arrays were used to directly genotype >800,000 genetic variants; the UK BiLEVE array (n = 49,922) and the UK Biobank Axiom array (n = 102,326), which share >95% of genomic markers. Both were enriched for known disease variants and include >600,000 variants specifically chosen to optimize imputation in populations of European descent. Quality control steps performed centrally by UK Biobank [34] included filtering less reliable genotyping results (e.g. setting genotypes to missing in batches where the SNP deviated from Hardy-Weinberg Equilibrium) and identifying poor quality samples (e.g. extreme heterozygosity and/or low call rate) according to Affymetrix recommendations.
Genotype imputation, also performed centrally by UK Biobank [35], was performed on a subset of genotype data, where SNPs with MAF <1% were excluded, leaving 641,018 variants. Phasing and imputation were performed separately, using a custom SHAPEIT2 algorithm to handle large sample size and IMPUTE2 with a combined UK10K and 1000 Genomes Phase 3 reference panel. This resulted in 73,355,667 genetic variations available in the dataset.
We performed a GWAS in 116,666 white/British participants (by self-report & genetic ancestry) with complete genetic data to determine genetic variants associated with RDW. After quality control and filtering (we included autosomal variants with MAF ≥0.1%, missingness <1.5%, imputation quality >0.4 and with Hardy-Weinberg equilibrium (HWE) p>1x10-6 within the white/British participants) 16,889,199 imputed genetic variants were available for GWAS analysis: methods described in detail previously [34–36]. We also utilized data directly from the microarrays for variants on the X (n = 19,381) and Y (n = 284) chromosomes, and on the mitochondrial genome (n = 135), which were unavailable in the imputed dataset.
GWAS was performed using BOLT-LMM, a software that uses linear mixed-effect models to determine associations between each variant and the outcome, incorporating genetic relatedness [33]. We provided all directly genotyped SNPs that were present in the imputed data to BOLT-LMM to construct the model. For X, Y and mitochondrial variants we used Plink (v1.9) [37] to determine the associations between genotype and phenotype, with additional adjustment for the first 5 principal components from FlashPCA [38] based on 95,535 independent SNPs (pairwise r2<0.1) with MAF >2.5%, missingness <1.5%, and HWE p>1x10-6.
The main outcome of the GWAS was RDW residuals from a linear regression model adjusted for age, sex, and assessment center, and were quantile-normalized prior to analysis. Models were adjusted for array type (two different Affymetrix arrays were used, which are >95% identical) at run time. Variants were classed as significant if the p-value for the association with RDW was less than 5x10-8.
Identifying conditionally independent GWAS signals
Many of the identified variants are correlated and may therefore not be independent; we used conditional analysis to determine independent signals by adjusting each variant in a locus for the most significant variant in that locus (loci defined as runs of SNPs separated by <2Mb on a chromosome). This process was repeated until only “conditional SNPs” remained that were significantly associated with RDW independent of one another.
Conditional SNPs were checked for their consistency of association with RDW in two sensitivity analyses: once excluding participants with prevalent anemia (either clinical diagnosis, self-report, or raised hemoglobin), and in a separate analysis excluding participants <60 years of age.
GWAS-significant SNPs follow-up
The determine the gene-locations and possible functional consequences of the conditionally independent SNPs we submitted them to the UCSC Variant Annotation Integrator (https://genome.ucsc.edu/cgi-bin/hgVai), which combines information from several sources to determine the probable effect of a genetic variant, including on specific genes (for example intronic, missense, splice site, intergenic etc.) and PolyPhen-2 (a tool for predicting the impact of amino acid substitutions on the protein product). We used data from the “Blood eQTL browser” (http://genenetwork.nl/bloodeqtlbrowser) published by Westra et al. which reports associations between genetic variants and gene expression in whole blood [25] to determine whether the independent SNPs affect gene expression in human whole blood.
The GWAS catalogue of published variant-trait associations was searched for all SNPs (not just conditional SNPs) associated with RDW to determine which loci had previously been associated with another trait in GWAS analyses (p<5x10-8), and which were novel [39]. We used the UCSC `liftOver`tool (https://genome.ucsc.edu/cgi-bin/hgLiftOver) to match the genomic coordinates between the GWAS catalogue (GRCh38) and the UK Biobank genetics data (GRCh37).
We used the software package MAGENTA [13] to determine whether any biological pathways were enriched in the GWAS results, using the pathways database file GO_PANTHER_INGENUITY_KEGG_REACTOME_BIOCARTA. These methods involve adjustment for multiple testing and include information on the expected and observed number of significant genes in each pathway.
Genetic risk scores
Twenty Genetic Risk Scores (GRS) were computed for each participant based on the number of trait-raising alleles they have for a particular phenotype, such as LDL cholesterol or type-2 diabetes. The derivation of each individual score is described in the S1 Methods document. To compute the GRS we used Plink function `scores`[37] which computes the number of trait-raising alleles (0, 1 or 2) in each participant, multiplied (weighted) by the effect size (coefficient or odds ratio) from the previously published study. Each of the 20 GRS computed was tested for its association with RDW using linear regression models, adjusted for age, sex, assessment center, genotype array, and population stratification (using the first 5 principal components (PCs)).
Supporting information
S1 Table. Genetics variants (signals) independently associated with RDW.
https://doi.org/10.1371/journal.pone.0185083.s001
(XLSX)
S2 Table. Variant annotation information for the independent signals associated with RDW.
https://doi.org/10.1371/journal.pone.0185083.s002
(XLSX)
S3 Table. Whole blood eQTL information for the RDW signals.
https://doi.org/10.1371/journal.pone.0185083.s003
(XLSX)
S4 Table. Known associations from the GWAS catalogue for the RDW signals.
https://doi.org/10.1371/journal.pone.0185083.s004
(XLSX)
S5 Table. Pathways significantly enriched in the genetic variants associated with RDW.
This analysis using the MAGENTA software included all genetic variants used in the analysis, not just the significant independent signals.
https://doi.org/10.1371/journal.pone.0185083.s005
(XLSX)
S6 Table. Genetic risk score associations with RDW.
https://doi.org/10.1371/journal.pone.0185083.s006
(XLSX)
Acknowledgments
This research has been conducted using the UK Biobank Resource under Application Number 14631. We thank UK Biobank participants and coordinators for this dataset.
References
- 1. Patel K V., Ferrucci L, Ershler WB, Longo DL, Guralnik JM. Red blood cell distribution width and the risk of death in middle-aged and older adults. Arch Intern Med. 2009;169: 515–23. pmid:19273783
- 2. Patel K V, Semba RD, Ferrucci L, Newman AB, Fried LP, Wallace RB, et al. Red cell distribution width and mortality in older adults: a meta-analysis. J Gerontol A Biol Sci Med Sci. 2010;65: 258–65. pmid:19880817
- 3. Turcato G, Serafini V, Dilda A, Bovo C, Caruso B, Ricci G, et al. Red blood cell distribution width independently predicts medium-term mortality and major adverse cardiac events after an acute coronary syndrome. Ann Transl Med. 2016;4: 254–254. pmid:27500155
- 4. Danese E, Lippi G, Montagnana M. Red blood cell distribution width and cardiovascular diseases. J Thorac Dis. 2015;7: E402–11. pmid:26623117
- 5. Inuzuka R, Abe J. Red Blood Cell Distribution Width as a Link Between Ineffective Erythropoiesis and Chronic Inflammation in Heart Failure. Circ J. 2015;79: 3–4. pmid:25808125
- 6. Lippi G, Sanchis-Gomar F, Danese E, Montagnana M. Association of red blood cell distribution width with plasma lipids in a general population of unselected outpatients. Kardiol Pol. 2013;71: 931–6. pmid:24065422
- 7. Sharma R, Agrawal V V. The Relationship Between Red Blood Cell Distribution Width (RDW CV) And C Reactive Protein (CRP) With The Clinical Outcomes In Non-St Elevation Myocardial Infarction And Unstable Angina Pectoris: A 6 Months Follow Up Study. Int Cardiovasc Forum J. 2015;2: 27–31.
- 8. Salvagno GL, Sanchis-Gomar F, Picanza A, Lippi G. Red blood cell distribution width: A simple parameter with multiple clinical applications. Crit Rev Clin Lab Sci. 2015;52: 86–105. pmid:25535770
- 9. Astle WJ, Elding H, Jiang T, Allen D, Ruklisa D, Mann AL, et al. The Allelic Landscape of Human Blood Cell Trait Variation and Links to Common Complex Disease. Cell. 2016;167: 1415–1429.e19. pmid:27863252
- 10. Yang J, Weedon MN, Purcell S, Lettre G, Estrada K, Willer CJ, et al. Genomic inflation factors under polygenic inheritance. Eur J Hum Genet. 2011;19: 807–812. pmid:21407268
- 11. Rebhan M, Chalifa-Caspi V, Prilusky J, Lancet D. GeneCards: integrating information about genes, proteins and diseases. Trends Genet. 1997;13: 163. Available: http://www.ncbi.nlm.nih.gov/pubmed/9097728 pmid:9097728
- 12.
NCBI ClinVar. ClinVar record for rs738409 [Internet]. [cited 15 Feb 2017]. Available: www.ncbi.nlm.nih.gov/clinvar/RCV000299548
- 13. Segrè A V, DIAGRAM Consortium, MAGIC investigators, Groop L, Mootha VK, Daly MJ, et al. Common inherited variation in mitochondrial genes is not enriched for associations with type 2 diabetes or related glycemic traits. PLoS Genet. 2010;6. pmid:20714348
- 14. Nikpay M, Goel A, Won H-H, Hall LM, Willenborg C, Kanoni S, et al. A comprehensive 1000 Genomes–based genome-wide association meta-analysis of coronary artery disease. Nat Genet. 2015;47: 1121–1130. pmid:26343387
- 15. Bowden J, Davey Smith G, Burgess S. Mendelian randomization with invalid instruments: effect estimation and bias detection through Egger regression. Int J Epidemiol. 2015;44: 512–25. pmid:26050253
- 16. Cheng CK-W, Chan J, Cembrowski GS, van Assendelft OW. Complete blood count reference interval diagrams derived from NHANES III: stratification by age, sex, and race. Lab Hematol. 2004;10: 42–53. Available: http://www.ncbi.nlm.nih.gov/pubmed/15070217 pmid:15070217
- 17. Lippi G, Cervellin G, Sanchis-Gomar F. Red blood cell distribution width and cardiovascular disorders. Does it really matter which comes first, the chicken or the egg? Int J Cardiol. Elsevier Ireland Ltd; 2016;206: 129–30. pmid:26788688
- 18. López-Otín C, Blasco MA, Partridge L, Serrano M, Kroemer G. The hallmarks of aging. Cell. 2013;153: 1194–217. pmid:23746838
- 19. Chakrabarti RS, Ingham SA, Kozlitina J, Gay A, Cohen JC, Radhakrishnan A, et al. Variability of cholesterol accessibility in human red blood cells measured using a bacterial cholesterol-binding toxin. Elife. 2017;6. pmid:28169829
- 20.
Fyhrquist FY, Saijonmaa OJ. Modifiable Factors Influencing Telomere Length and Aging. In: Bondy SC, Campbell A, editors. Inflammation, Aging, and Oxidative Stress. Cham: Springer International Publishing; 2016. pp. 67–80. https://doi.org/10.1007/978-3-319-33486-8_4
- 21. Telomeres Mendelian Randomization Collaboration, Haycock PC, Burgess S, Nounu A, Zheng J, Okoli GN, et al. Association Between Telomere Length and Risk of Cancer and Non-Neoplastic Diseases: A Mendelian Randomization Study. JAMA Oncol. 2017; pmid:28241208
- 22. Kozlitina J, Garcia CK. Red blood cell size is inversely associated with leukocyte telomere length in a large multi-ethnic population. PLoS One. 2012;7: e51046. pmid:23226558
- 23. Megel C, Morelle G, Lalande S, Duchêne A-M, Small I, Maréchal-Drouard L. Surveillance and cleavage of eukaryotic tRNAs. Int J Mol Sci. 2015;16: 1873–93. pmid:25599528
- 24. Harries LW, Hernandez D, Henley W, Wood AR, Holly AC, Bradley-Smith RM, et al. Human aging is characterized by focused changes in gene expression and deregulation of alternative splicing. Aging Cell. 2011;10: 868–78. pmid:21668623
- 25. Westra H-J, Peters MJ, Esko T, Yaghootkar H, Schurmann C, Kettunen J, et al. Systematic identification of trans eQTLs as putative drivers of known disease associations. Nat Genet. Nature Publishing Group; 2013;45: 1238–43. pmid:24013639
- 26. Locke AE, Kahali B, Berndt SI, Justice AE, Pers TH, Day FR, et al. Genetic studies of body mass index yield new insights for obesity biology. Nature. 2015;518: 197–206. pmid:25673413
- 27. Kato N, Loh M, Takeuchi F, Verweij N, Wang X, Zhang W, et al. Trans-ancestry genome-wide association study identifies 12 genetic loci influencing blood pressure and implicates a role for DNA methylation. Nat Genet. 2015;47: 1282–93. pmid:26390057
- 28. Wang G-S, Cooper T a. Splicing in disease: disruption of the splicing code and the decoding machinery. Nat Rev Genet. 2007;8: 749–761. pmid:17726481
- 29. Manolio TA, Collins R. Enhancing the feasibility of large cohort studies. JAMA. 2010;304: 2290–1. pmid:21098774
- 30. Sudlow C, Gallacher J, Allen N, Beral V, Burton P, Danesh J, et al. UK biobank: an open access resource for identifying the causes of a wide range of complex diseases of middle and old age. PLoS Med. 2015;12: e1001779. pmid:25826379
- 31.
Sheard S, Nicholls R, Froggat J. UK Biobank Haematology Data Companion Document [Internet]. 2016 [cited 25 May 2016]. Available: https://biobank.ctsu.ox.ac.uk/crystal/docs/haematology.pdf
- 32.
WHO. WHO definition of anemia by hemoglobin concentration [Internet]. [cited 20 Oct 2016]. Available: http://www.who.int/vmnis/indicators/haemoglobin.pdf
- 33. Loh P-R, Tucker G, Bulik-Sullivan BK, Vilhjálmsson BJ, Finucane HK, Salem RM, et al. Efficient Bayesian mixed-model analysis increases association power in large cohorts. Nat Genet. Nature Publishing Group; 2015;47: 284–90. pmid:25642633
- 34.
UK Biobank. UK Biobank genotype collections and QC on 150k participants, April 2014 [Internet]. 2014. Available: http://www.ukbiobank.ac.uk/wp-content/uploads/2014/04/UKBiobank_genotyping_QC_documentation-web.pdf
- 35.
UK Biobank. UK Biobank imputation on 150k participants, May 2015 [Internet]. 2015. Available: http://www.ukbiobank.ac.uk/wp-content/uploads/2014/04/imputation_documentation_May2015.pdf
- 36. Pilling LC, Atkins JL, Bowman K, Jones SE, Tyrrell J, Beaumont RN, et al. Human longevity is influenced by many genetic variants: evidence from 75,000 UK Biobank participants. Aging (Albany NY). 2016;8: 547–63.
- 37. Purcell S, Neale B, Todd-Brown K, Thomas L, Ferreira MAR, Bender D, et al. PLINK: a tool set for whole-genome association and population-based linkage analyses. Am J Hum Genet. 2007;81: 559–75. pmid:17701901
- 38. Abraham G, Inouye M. Fast principal component analysis of large-scale genome-wide data. PLoS One. 2014;9: e93766. pmid:24718290
- 39.
Burdett T (EBI), Hall P (NHGRI), Hastings E (EBI), Hindorff L (NHGRI), Junkins H (NHGRI), Klemm A (NHGRI), et al. The NHGRI-EBI Catalog of published genome-wide association studies [Internet]. 2016 [cited 23 Oct 2016]. Available: www.ebi.ac.uk/gwas