GWAS of blood cell traits identifies novel associated loci and epistatic interactions in Caucasian and African-American children

SPLENDID incorporates continuous genetic ancestry in biobank-scale data to improve polygenic risk prediction across diverse populations

Polygenic risk scores are widely used in disease risk stratification, but their accuracy varies across diverse populations. Recent methods large-scale leverage multi-ancestry data to improve accuracy in under-represented populations but require labelling individuals by ancestry for prediction. This poses challenges for practical use, as clinical practices are typically not based on ancestry. We propose SPLENDID, a novel penalized regression framework for diverse biobank-scale data. Our method utilizes ancestry principal component interactions to model genetic ancestry as a continuum within a single prediction model for all ancestries, eliminating the need for discrete labels. In extensive simulations and analyses of 9 traits from the All of Us Research Program (N=224,364) and UK Biobank (N=340,140), SPLENDID significantly outperformed existing methods in prediction accuracy and model sparsity. By directly incorporating continuous genetic ancestry in model training, SPLENDID stands as a valuable tool for robust risk prediction across diverse populations and fairer clinical implementation.

Ribosomal DNA copy number variation associates with hematological profiles and renal function in the UK Biobank

The phenotypic impact of genetic variation of repetitive features in the human genome is currently understudied. One such feature is the multi-copy 47S ribosomal DNA (rDNA) that codes for rRNA components of the ribosome. Here, we present an analysis of rDNA copy number (CN) variation in the UK Biobank (UKB). From the first release of UKB whole-genome sequencing (WGS) data, a discovery analysis in White British individuals reveals that rDNA CN associates with altered counts of specific blood cell subtypes, such as neutrophils, and with the estimated glomerular filtration rate, a marker of kidney function. Similar trends are observed in other ancestries. A range of analyses argue against reverse causality or common confounder effects, and all core results replicate in the second UKB WGS release. Our work demonstrates that rDNA CN is a genetic influence on trait variance in humans.

The association of TMPRSS6 gene polymorphism with iron status in Egyptian children (a pilot study)

Several studies have shown association of single nucleotide polymorphisms (SNPs) of hepcidin regulatory pathways genes with impaired iron status. The most common is in the TMPRSS6 gene. In Africa, very few studies have been reported. We aimed to investigate the correlation between the common SNPs in the transmembrane protease, serine 6 (TMPRSS6) gene and iron indicators in a sample of Egyptian children for identifying the suitable candidate for iron supplementation.Patients and methods One hundred and sixty children aged 5-13 years were included & classified into iron deficient, iron deficient anemia and normal healthy controls. All were subjected to assessment of serum iron, serum ferritin, total iron binding capacity, complete blood count, reticulocyte count, serum soluble transferrin receptor and serum hepcidin. Molecular study of TMPRSS6 genotyping polymorphisms (rs4820268, rs855791 and rs11704654) were also evaluated.Results There was an association of iron deficiency with AG of rs855791 SNP, (P = 0.01). The minor allele frequency for included children were 0.43, 0.45 & 0.17 for rs4820268, rs855791 & rs11704654 respectively. Genotype GG of rs4820268 expressed the highest hepcidin gene expression fold, the lowest serum ferroportin & iron store compared to AA and AG genotypes (p = 0.05, p = 0.05, p = 0.03 respectively). GG of rs855791 had lower serum ferritin than AA (p = 0.04), lowest iron store & highest serum hepcidin compared to AA and AG genotypes (p = 0.04, p = 0.01 respectively). Children having CC of rs11704654 had lower level of hemoglobin, serum ferritin and serum hepcidin compared with CT genotype (p = 0.01, p = 0.01, p = 0.02) respectively.Conclusion Possible contribution of SNPs (rs855791, rs4820268 and rs11704654) to low iron status.

Comparison of the loci associated with HbA1c and blood glucose levels identified by a genome-wide association study in the Japanese population

Aims

Hemoglobin A1c (HbA1c) levels are widely employed to diagnose diabetes. However, estimates of the heritability of HbA1c and glucose levels are different. Therefore, we explored HbA1c- and blood glucose-associated loci in a non-diabetic Japanese population.

Methods

We conducted a two-stage genome-wide association study (GWAS) on variants associated with HbA1c and blood glucose levels in a Japanese population. In the initial stage, data of 4911 participants of the Japan Multi-Institutional Collaborative Cohort (J-MICC) were subjected to discovery analysis. In the second stage, two datasets from the Tohoku Medical Megabank project, with 8175 and 40,519 participants, were used for the replication study. Association of the imputed variants with HbA1c and blood glucose levels was determined via linear regression analyses adjusted for age, sex, body mass index (BMI), smoking, and genetic principal components (PC1-PC10). Moreover, we performed a BMI-stratified GWAS on HbA1c levels in the J-MICC. The discovery analysis and BMI-stratified GWAS results were validated with re-analyses of normalized HbA1c levels adjusted for site in addition to the above, and blood glucose adjusted for fasting time as an additional covariate.

Results

Genetic variants associated with HbA1c levels were identified in KCNQ1 and TMC6. None of the genetic variants associated with blood glucose levels in the discovery analysis were replicated. Association of rs2299620 in KCNQ1 with HbA1c levels showed heterogeneity between individuals with BMI ≥ 25 kg/m2 and BMI < 25 kg/m2.

Conclusions

The variant rs2299620 in KCNQ1 might affect HbA1c levels differentially based on BMI grouping in the Japanese population.

Supplementary Information

The online version contains supplementary material available at 10.1007/s13340-023-00618-0.

A pleiotropic hypoxia-sensitive EPAS1 enhancer is disrupted by adaptive alleles in Tibetans

In Tibetans, noncoding alleles in EPAS1-whose protein product hypoxia-inducible factor 2α (HIF-2α) drives the response to hypoxia-carry strong signatures of positive selection; however, their functional mechanism has not been systematically examined. Here, we report that high-altitude alleles disrupt the activity of four EPAS1 enhancers in one or more cell types. We further characterize one enhancer (ENH5) whose activity is both allele specific and hypoxia dependent. Deletion of ENH5 results in down-regulation of EPAS1 and HIF-2α targets in acute hypoxia and in a blunting of the transcriptional response to sustained hypoxia. Deletion of ENH5 in mice results in dysregulation of gene expression across multiple tissues. We propose that pleiotropic adaptive effects of the Tibetan alleles in EPAS1 underlie the strong selective signal at this gene.

Genome-Wide Association Study on the Hematological Phenotypic Characteristics of the Han Population from Northwest China

Background

Hematological characteristics have positive reference value as clinical indicators in the evaluation of various diseases. The purpose of this study was to determine the gene loci associated with 20 hematological phenotypes in the Han population from northwest China.

Methods

A genome-wide association study (GWAS) was conducted on hematological indicators of 1005 Han people from northwest China. Genotyping was performed with a GeneTitan multichannel instrument and Axiom Analysis Suite 6.0. Using the 1000 Genomes Project (phase 3) as a reference, haplotype imputation was performed with IMPUTE2. SNVs (single nucleotide variants) significantly associated with hematological phenotypes were identified. The top SNV (p < 5E-7) was then selected for replication detection.

Results

Ninety genetic variations identified in the GWAS were significantly associated with hematological indicators. Among them, only rs35289401 (CCDC157) was significantly associated (genome-wide) with red blood cell distribution width (RDW) (p = 4.21E-08). The fourteen top SNVs were selected for replication verification and were significantly associated with hematological phenotypes. However, only HBS1 L-MYB rs1331309 was significantly associated with the mean hemoglobin content (p = 6.42E-07). We also found that the mean corpuscular hemoglobin (MCH) level in the rs1331309 GG/GT genotype was significantly higher than that in the TT genotype (p = 0.023).

Conclusion

The GWAS identified a total of 90 genetic variants significantly associated with hematological phenotypic indicators. In particular, rs1331309 (HBS1 L-MYB) is expected to be a biomarker for monitoring the dynamics of MCH levels. This study provides a reference for related studies on the genetic structure of hematological characteristics. It provides a valuable reference for the clinical diagnosis or prediction of a variety of diseases.

Uncertainty quantification in variable selection for genetic fine-mapping using bayesian neural networks

In this paper, we propose a new approach for variable selection using a collection of Bayesian neural networks with a focus on quantifying uncertainty over which variables are selected. Motivated by fine-mapping applications in statistical genetics, we refer to our framework as an "ensemble of single-effect neural networks" (ESNN) which generalizes the "sum of single effects" regression framework by both accounting for nonlinear structure in genotypic data (e.g., dominance effects) and having the capability to model discrete phenotypes (e.g., case-control studies). Through extensive simulations, we demonstrate our method's ability to produce calibrated posterior summaries such as credible sets and posterior inclusion probabilities, particularly for traits with genetic architectures that have significant proportions of non-additive variation driven by correlated variants. Lastly, we use real data to demonstrate that the ESNN framework improves upon the state of the art for identifying true effect variables underlying various complex traits.

Identification of ATP2B4 Regulatory Element Containing Functional Genetic Variants Associated with Severe Malaria

Genome-wide association studies for severe malaria (SM) have identified 30 genetic variants mostly located in non-coding regions. Here, we aimed to identify potential causal genetic variants located in these loci and demonstrate their functional activity. We systematically investigated the regulatory effect of the SNPs in linkage disequilibrium (LD) with the malaria-associated genetic variants. Annotating and prioritizing genetic variants led to the identification of a regulatory region containing five ATP2B4 SNPs in LD with rs10900585. We found significant associations between SM and rs10900585 and our candidate SNPs (rs11240734, rs1541252, rs1541253, rs1541254, and rs1541255) in a Senegalese population. Then, we demonstrated that both individual SNPs and the combination of SNPs had regulatory effects. Moreover, CRISPR/Cas9-mediated deletion of this region decreased ATP2B4 transcript and protein levels and increased Ca²⁺ intracellular concentration in the K562 cell line. Our data demonstrate that severe malaria-associated genetic variants alter the expression of ATP2B4 encoding a plasma membrane calcium-transporting ATPase 4 (PMCA4) expressed on red blood cells. Altering the activity of this regulatory element affects the risk of SM, likely through calcium concentration effect on parasitaemia.

Whole genome sequence analysis of platelet traits in the NHLBI Trans-Omics for Precision Medicine (TOPMed) initiative

Platelets play a key role in thrombosis and hemostasis. Platelet count (PLT) and mean platelet volume (MPV) are highly heritable quantitative traits, with hundreds of genetic signals previously identified, mostly in European ancestry populations. We here utilize whole genome sequencing (WGS) from NHLBI's Trans-Omics for Precision Medicine initiative (TOPMed) in a large multi-ethnic sample to further explore common and rare variation contributing to PLT (n = 61 200) and MPV (n = 23 485). We identified and replicated secondary signals at MPL (rs532784633) and PECAM1 (rs73345162), both more common in African ancestry populations. We also observed rare variation in Mendelian platelet-related disorder genes influencing variation in platelet traits in TOPMed cohorts (not enriched for blood disorders). For example, association of GP9 with lower PLT and higher MPV was partly driven by a pathogenic Bernard-Soulier syndrome variant (rs5030764, p.Asn61Ser), and the signals at TUBB1 and CD36 were partly driven by loss of function variants not annotated as pathogenic in ClinVar (rs199948010 and rs571975065). However, residual signal remained for these gene-based signals after adjusting for lead variants, suggesting that additional variants in Mendelian genes with impacts in general population cohorts remain to be identified. Gene-based signals were also identified at several genome-wide association study identified loci for genes not annotated for Mendelian platelet disorders (PTPRH, TET2, CHEK2), with somatic variation driving the result at TET2. These results highlight the value of WGS in populations of diverse genetic ancestry to identify novel regulatory and coding signals, even for well-studied traits like platelet traits.

Secondary analyses for genome-wide association studies using expression quantitative trait loci

Genome-wide association studies (GWAS) have successfully identified thousands of single nucleotide polymorphisms (SNPs) associated with complex traits; however, the identified SNPs account for a fraction of trait heritability, and identifying the functional elements through which genetic variants exert their effects remains a challenge. Recent evidence suggests that SNPs associated with complex traits are more likely to be expression quantitative trait loci (eQTL). Thus, incorporating eQTL information can potentially improve power to detect causal variants missed by traditional GWAS approaches. Using genomic, transcriptomic, and platelet phenotype data from the Genetic Study of Atherosclerosis Risk family-based study, we investigated the potential to detect novel genomic risk loci by incorporating information from eQTL in the relevant target tissues (i.e., platelets and megakaryocytes) using established statistical principles in a novel way. Permutation analyses were performed to obtain family-wise error rates for eQTL associations, substantially lowering the genome-wide significance threshold for SNP-phenotype associations. In addition to confirming the well known association between PEAR1 and platelet aggregation, our eQTL-focused approach identified a novel locus (rs1354034) and gene (ARHGEF3) not previously identified in a GWAS of platelet aggregation phenotypes. A colocalization analysis showed strong evidence for a functional role of this eQTL.

Genome-Wide Association and Mendelian Randomization Analysis Reveal the Causal Relationship Between White Blood Cell Subtypes and Asthma in Africans

Background: White blood cell (WBC) traits and their subtypes such as basophil count (Bas), eosinophil count (Eos), lymphocyte count (Lym), monocyte count (Mon), and neutrophil counts (Neu) are known to be associated with diseases such as stroke, peripheral arterial disease, and coronary heart disease. Methods: We meta-analyze summary statistics from genome-wide association studies in 17,802 participants from the African Partnership for Chronic Disease Research (APCDR) and African ancestry individuals from the Blood Cell Consortium (BCX2) using GWAMA. We further carried out a Bayesian fine mapping to identify causal variants driving the association with WBC subtypes. To access the causal relationship between WBC subtypes and asthma, we conducted a two-sample Mendelian randomization (MR) analysis using summary statistics of the Consortium on Asthma among African Ancestry Populations (CAAPA: n cases = 7,009, n control = 7,645) as our outcome phenotype. Results: Our metanalysis identified 269 loci at a genome-wide significant value of (p = 5 × 10-9) in a composite of the WBC subtypes while the Bayesian fine-mapping analysis identified genetic variants that are more causal than the sentinel single-nucleotide polymorphism (SNP). We found for the first time five novel genes (LOC126987/MTCO3P14, LINC01525, GAPDHP32/HSD3BP3, FLG-AS1/HMGN3P1, and TRK-CTT13-1/MGST3) not previously reported to be associated with any WBC subtype. Our MR analysis showed that Mon (IVW estimate = 0.38, CI: 0.221, 0.539, p < 0.001), Neu (IVW estimate = 0.189, CI: 0.133, 0.245, p < 0.001), and WBCc (IVW estimate = 0.185, CI: 0.108, 0.262, p < 0.001) are associated with increased risk of asthma. However, there was no evidence of causal relationship between Lym and asthma risk. Conclusion: This study provides insight into the relationship between some WBC subtypes and asthma and potential route in the treatment of asthma and may further inform a new therapeutic approach.

Genetic variants associated with platelet count are predictive of human disease and physiological markers

Platelets play an important role in hemostasis and other aspects of vascular biology. We conducted a meta-analysis of platelet count GWAS using data on 536,974 Europeans and identified 577 independent associations. To search for mechanisms through which these variants affect platelets, we applied cis-expression quantitative trait locus, DEPICT and IPA analyses and assessed genetic sharing between platelet count and various traits using polygenic risk scoring. We found genetic sharing between platelet count and counts of other blood cells (except red blood cells), in addition to several other quantitative traits, including markers of cardiovascular, liver and kidney functions, height, and weight. Platelet count polygenic risk score was predictive of myeloproliferative neoplasms, rheumatoid arthritis, ankylosing spondylitis, hypertension, and benign prostate hyperplasia. Taken together, these results advance understanding of diverse aspects of platelet biology and how they affect biological processes in health and disease.

Evgenia Mikaelsdottir et al. report a study of variants associated with platelet count among European individuals where they identify 577 associations. They also report a genetic overlap between platelet count and human diseases, including myeloproliferative neoplasms, rheumatoid arthritis, and hypertension, as well as a genetic overlap between platelet count and various physiological markers.

Common host variation drives malaria parasite fitness in healthy human red cells

The replication of Plasmodium falciparum parasites within red blood cells (RBCs) causes severe disease in humans, especially in Africa. Deleterious alleles like hemoglobin S are well-known to confer strong resistance to malaria, but the effects of common RBC variation are largely undetermined. Here, we collected fresh blood samples from 121 healthy donors, most with African ancestry, and performed exome sequencing, detailed RBC phenotyping, and parasite fitness assays. Over one-third of healthy donors unknowingly carried alleles for G6PD deficiency or hemoglobinopathies, which were associated with characteristic RBC phenotypes. Among non-carriers alone, variation in RBC hydration, membrane deformability, and volume was strongly associated with P. falciparum growth rate. Common genetic variants in PIEZO1, SPTA1/SPTB, and several P. falciparum invasion receptors were also associated with parasite growth rate. Interestingly, we observed little or negative evidence for divergent selection on non-pathogenic RBC variation between Africans and Europeans. These findings suggest a model in which globally widespread variation in a moderate number of genes and phenotypes modulates P. falciparum fitness in RBCs.

Assessing the impact of alcohol consumption on the genetic contribution to mean corpuscular volume

The relationship between the genetic loci that influence mean corpuscular volume (MCV) and those associated with excess alcohol drinking are unknown. We used white British participants from the UK Biobank (n = 362 595) to assess the association between alcohol consumption and MCV, and whether this was modulated by genetic factors. Multivariable regression was applied to identify predictors of MCV. GWAS, with and without stratification for alcohol consumption, determined how genetic variants influence MCV. SNPs in ADH1B, ADH1C and ALDH1B were used to construct a genetic score to test the assumption that acetaldehyde formation is an important determinant of MCV. Additional investigations using mendelian randomisation and phenome-wide association analysis were conducted. Increasing alcohol consumption by 40 g/week resulted in a 0.30% (95% CI: 0.30 to 0.31%) increase in MCV (P < 1.0x10-320). Unstratified (irrespective of alcohol intake) GWAS identified 212 loci associated with MCV, of which 108 were novel. There was no heterogeneity of allelic effects by drinking status. No association was found between MCV and the genetic score generated from alcohol metabolising genes. Mendelian randomisation demonstrated a causal effect for alcohol on MCV. Seventy-one SNP-outcome pairs reached statistical significance in phenome-wide association analysis, with evidence of shared genetic architecture for MCV and thyroid dysfunction, and mineral metabolism disorders. MCV increases linearly with alcohol intake in a causal manner. Many genetic loci influence MCV, with new loci identified in this analysis that provide novel biological insights. However, there was no interaction between alcohol consumption and the allelic variants associated with MCV.

Association of common TMPRSS6 and TF gene variants with hepcidin and iron status in healthy rural Gambians

Genome-wide association studies in Europeans and Asians have identified numerous variants in the transmembrane protease serine 6 (TMPRSS6) and transferrin (TF) genes that are associated with changes in iron status. We sought to investigate the effects of common TMPRSS6 and TF gene SNPs on iron status indicators in a cohort of healthy Africans from rural Gambia. We measured iron biomarkers and haematology traits on individuals participating in the Keneba Biobank with genotype data on TMPRSS6 (rs2235321, rs855791, rs4820268, rs2235324, rs2413450 and rs5756506) and TF (rs3811647 and rs1799852), n = 1316. After controlling for inflammation, age and sex, we analysed the effects of carrying either single or multiple iron-lowering alleles on iron status. TMPRSS6 rs2235321 significantly affected plasma hepcidin concentrations (AA genotypes having lower hepcidin levels; F ratio 3.7, P = 0.014) with greater impact in individuals with low haemoglobin or ferritin. No other TMPRSS6 variant affected hepcidin. None of the TMPRSS6 variants nor a TMPRSS6 allele risk score affected other iron biomarkers or haematological traits. TF rs3811647 AA carriers had 21% higher transferrin (F ratio 16.0, P < 0.0001), 24% higher unsaturated iron-binding capacity (F ratio 12.8, P < 0.0001) and 25% lower transferrin saturation (F ratio 4.3, P < 0.0001) compared to GG carriers. TF rs3811647 was strongly associated with transferrin, unsaturated iron-binding capacity (UIBC) and transferrin saturation (TSAT) with a single allele effect of 8-12%. There was no association between either TF SNP and any haematological traits or iron biomarkers. We identified meaningful associations between TMPRSS6 rs2235321 and hepcidin and replicated the previous findings on the effects of TF rs3811647 on transferrin and iron binding capacity. However, the effects are subtle and contribute little to population variance. Further genetic and functional studies, including polymorphisms frequent in Africa populations, are needed to identify markers for genetically stratified approaches to prevention or treatment of iron deficiency anaemia.

RCL1 copy number variants are associated with a range of neuropsychiatric phenotypes

Mendelian and early-onset severe psychiatric phenotypes often involve genetic variants having a large effect, offering opportunities for genetic discoveries and early therapeutic interventions. Here, the index case is an 18-year-old boy, who at 14 years of age had a decline in cognitive functioning over the course of a year and subsequently presented with catatonia, auditory and visual hallucinations, paranoia, aggression, mood dysregulation, and disorganized thoughts. Exome sequencing revealed a stop-gain mutation in RCL1 (NM_005772.4:c.370 C > T, p.Gln124Ter), encoding an RNA 3'-terminal phosphate cyclase-like protein that is highly conserved across eukaryotic species. Subsequent investigations across two academic medical centers identified eleven additional cases of RCL1 copy number variations (CNVs) with varying neurodevelopmental or psychiatric phenotypes. These findings suggest that dosage variation of RCL1 contributes to a range of neurological and clinical phenotypes.

miR-125a-5p regulates megakaryocyte proplatelet formation via the actin-bundling protein L-plastin

There is heritability to interindividual variation in platelet count, and better understanding of the regulating genetic factors may provide insights for thrombopoiesis. MicroRNAs (miRs) regulate gene expression in health and disease, and megakaryocytes (MKs) deficient in miRs have lower platelet counts, but information about the role of miRs in normal human MK and platelet production is limited. Using genome-wide miR profiling, we observed strong correlations among human bone marrow MKs, platelets, and differentiating cord blood-derived MK cultures, and identified MK miR-125a-5p as associated with human platelet number but not leukocyte or hemoglobin levels. Overexpression and knockdown studies showed that miR-125a-5p positively regulated human MK proplatelet (PP) formation in vitro. Inhibition of miR-125a-5p in vivo lowered murine platelet counts. Analyses of MK and platelet transcriptomes identified LCP1 as a miR-125a-5p target. LCP1 encodes the actin-bundling protein, L-plastin, not previously studied in MKs. We show that miR-125a-5p directly targets and reduces expression of MK L-plastin. Overexpression and knockdown studies show that L-plastin promotes MK progenitor migration, but negatively correlates with human platelet count and inhibits MK PP formation (PPF). This work provides the first evidence for the actin-bundling protein, L-plastin, as a regulator of human MK PPF via inhibition of the late-stage MK invagination system, podosome and PPF, and PP branching. We also provide resources of primary and differentiating MK transcriptomes and miRs associated with platelet counts. miR-125a-5p and L-plastin may be relevant targets for increasing in vitro platelet manufacturing and for managing quantitative platelet disorders.

Novel Approaches to Fine-Tune Therapeutic Targeting of Platelets in Atherosclerosis: A Critical Appraisal

The pathogenesis of atherosclerotic vascular disease is driven by a multitude of risk factors intertwining metabolic and inflammatory pathways. Increasing knowledge about platelet biology sheds light on how platelets take part in these processes from early to later stages of plaque development. Recent insights from experimental studies and mouse models substantiate platelets as initiators and amplifiers in atherogenic leukocyte recruitment. These studies are complemented by results from genetics studies shedding light on novel molecular mechanisms which provide an interesting prospect as novel targets. For instance, experimental studies provide further details how platelet-decorated von Willebrand factor tethered to activated endothelial cells plays a role in atherogenic monocyte recruitment. Novel aspects of platelets as atherogenic inductors of neutrophil extracellular traps and particularities in signaling pathways such as cyclic guanosine monophosphate and the inhibitory adaptor molecule SHB23/LNK associating platelets with atherogenesis are shared. In summary, it was our intention to balance insights from recent experimental data that support a plausible role for platelets in atherogenesis against a paucity of clinical evidence needed to validate this concept in humans.

Genomic heterogeneity in myeloproliferative neoplasms and applications to clinical practice

The myeloproliferative neoplasms (MPN) polycythaemia vera, essential thrombocythaemia and primary myelofibrosis are chronic myeloid disorders associated most often with mutations in JAK2, MPL and CALR, and in some patients with additional acquired genomic lesions. Whilst the molecular mechanisms downstream of these mutations are now clearer, it is apparent that clinical phenotype in MPN is a product of complex interactions, acting between individual mutations, between disease subclones, and between the tumour and background host factors. In this review we first discuss MPN phenotypic driver mutations and the factors that interact with them to influence phenotype. We consider the importance of ongoing studies of clonal haematopoiesis, which may inform a better understanding of why MPN develop in specific individuals. We then consider how best to deploy genomic testing in a clinical environment and the challenges as well as opportunities that may arise from more routine, comprehensive genomic analysis of patients with MPN.

Estimation of non-null SNP effect size distributions enables the detection of enriched genes underlying complex traits

Traditional univariate genome-wide association studies generate false positives and negatives due to difficulties distinguishing associated variants from variants with spurious nonzero effects that do not directly influence the trait. Recent efforts have been directed at identifying genes or signaling pathways enriched for mutations in quantitative traits or case-control studies, but these can be computationally costly and hampered by strict model assumptions. Here, we present gene-ε, a new approach for identifying statistical associations between sets of variants and quantitative traits. Our key insight is that enrichment studies on the gene-level are improved when we reformulate the genome-wide SNP-level null hypothesis to identify spurious small-to-intermediate SNP effects and classify them as non-causal. gene-ε efficiently identifies enriched genes under a variety of simulated genetic architectures, achieving greater than a 90% true positive rate at 1% false positive rate for polygenic traits. Lastly, we apply gene-ε to summary statistics derived from six quantitative traits using European-ancestry individuals in the UK Biobank, and identify enriched genes that are in biologically relevant pathways.

Ancestry-specific associations identified in genome-wide combined-phenotype study of red blood cell traits emphasize benefits of diversity in genomics

BACKGROUND

Quantitative red blood cell (RBC) traits are highly polygenic clinically relevant traits, with approximately 500 reported GWAS loci. The majority of RBC trait GWAS have been performed in European- or East Asian-ancestry populations, despite evidence that rare or ancestry-specific variation contributes substantially to RBC trait heritability. Recently developed combined-phenotype methods which leverage genetic trait correlation to improve statistical power have not yet been applied to these traits. Here we leveraged correlation of seven quantitative RBC traits in performing a combined-phenotype analysis in a multi-ethnic study population.

RESULTS

We used the adaptive sum of powered scores (aSPU) test to assess combined-phenotype associations between ~ 21 million SNPs and seven RBC traits in a multi-ethnic population (maximum n = 67,885 participants; 24% African American, 30% Hispanic/Latino, and 43% European American; 76% female). Thirty-nine loci in our multi-ethnic population contained at least one significant association signal (p < 5E-9), with lead SNPs at nine loci significantly associated with three or more RBC traits. A majority of the lead SNPs were common (MAF > 5%) across all ancestral populations. Nineteen additional independent association signals were identified at seven known loci (HFE, KIT, HBS1L/MYB, CITED2/FILNC1, ABO, HBA1/2, and PLIN4/5). For example, the HBA1/2 locus contained 14 conditionally independent association signals, 11 of which were previously unreported and are specific to African and Amerindian ancestries. One variant in this region was common in all ancestries, but exhibited a narrower LD block in African Americans than European Americans or Hispanics/Latinos. GTEx eQTL analysis of all independent lead SNPs yielded 31 significant associations in relevant tissues, over half of which were not at the gene immediately proximal to the lead SNP.

CONCLUSION

This work identified seven loci containing multiple independent association signals for RBC traits using a combined-phenotype approach, which may improve discovery in genetically correlated traits. Highly complex genetic architecture at the HBA1/2 locus was only revealed by the inclusion of African Americans and Hispanics/Latinos, underscoring the continued importance of expanding large GWAS to include ancestrally diverse populations.

Mosaic chromosome Y loss is associated with alterations in blood cell counts in UK Biobank men

Mosaic loss of Y chromosome (mLOY) is the most frequently detected somatic copy number alteration in leukocytes of men. In this study, we investigate blood cell counts as a potential mechanism linking mLOY to disease risk in 206,353 UK males. Associations between mLOY, detected by genotyping arrays, and blood cell counts were assessed by multivariable linear models adjusted for relevant risk factors. Among the participants, mLOY was detected in 39,809 men. We observed associations between mLOY and reduced erythrocyte count (−0.009 [−0.014, −0.005] × 10¹² cells/L, p = 2.75 × 10⁻⁵) and elevated thrombocyte count (5.523 [4.862, 6.183] × 10⁹ cells/L, p = 2.32 × 10⁻⁶⁰) and leukocyte count (0.218 [0.198, 0.239] × 10⁹ cells/L, p = 9.22 × 10⁻⁹⁵), particularly for neutrophil count (0.174 × [0.158, 0.190]10⁹ cells/L, p = 1.24 × 10⁻⁹⁹) and monocyte count (0.021 [0.018 to 0.024] × 10⁹ cells/L, p = 6.93 × 10⁻⁵⁷), but lymphocyte count was less consistent (0.016 [0.007, 0.025] × 10⁹ cells/L, p = 8.52 × 10⁻⁴). Stratified analyses indicate these associations are independent of the effects of aging and smoking. Our findings provide population-based evidence for associations between mLOY and blood cell counts that should stimulate investigation of the underlying biological mechanisms linking mLOY to cancer and chronic disease risk.

Use of >100,000 NHLBI Trans-Omics for Precision Medicine (TOPMed) Consortium whole genome sequences improves imputation quality and detection of rare variant associations in admixed African and Hispanic/Latino populations

Most genome-wide association and fine-mapping studies to date have been conducted in individuals of European descent, and genetic studies of populations of Hispanic/Latino and African ancestry are limited. In addition, these populations have more complex linkage disequilibrium structure. In order to better define the genetic architecture of these understudied populations, we leveraged >100,000 phased sequences available from deep-coverage whole genome sequencing through the multi-ethnic NHLBI Trans-Omics for Precision Medicine (TOPMed) program to impute genotypes into admixed African and Hispanic/Latino samples with genome-wide genotyping array data. We demonstrated that using TOPMed sequencing data as the imputation reference panel improves genotype imputation quality in these populations, which subsequently enhanced gene-mapping power for complex traits. For rare variants with minor allele frequency (MAF) < 0.5%, we observed a 2.3- to 6.1-fold increase in the number of well-imputed variants, with 11-34% improvement in average imputation quality, compared to the state-of-the-art 1000 Genomes Project Phase 3 and Haplotype Reference Consortium reference panels. Impressively, even for extremely rare variants with minor allele count <10 (including singletons) in the imputation target samples, average information content rescued was >86%. Subsequent association analyses of TOPMed reference panel-imputed genotype data with hematological traits (hemoglobin (HGB), hematocrit (HCT), and white blood cell count (WBC)) in ~21,600 African-ancestry and ~21,700 Hispanic/Latino individuals identified associations with two rare variants in the HBB gene (rs33930165 with higher WBC [p = 8.8x10-15] in African populations, rs11549407 with lower HGB [p = 1.5x10-12] and HCT [p = 8.8x10-10] in Hispanics/Latinos). By comparison, neither variant would have been genome-wide significant if either 1000 Genomes Project Phase 3 or Haplotype Reference Consortium reference panels had been used for imputation. Our findings highlight the utility of the TOPMed imputation reference panel for identification of novel rare variant associations not previously detected in similarly sized genome-wide studies of under-represented African and Hispanic/Latino populations.

Genetic Study of Severe Prolonged Lymphopenia in Multiple Sclerosis Patients Treated With Dimethyl Fumarate

In delayed-release dimethyl fumarate (DMF)-treated patients, absolute lymphocyte count (ALC) often declines in the first year and stabilizes thereafter; early declines have been associated with development of severe prolonged lymphopenia (SPL). Prolonged moderate or severe lymphopenia is a known risk factor for progressive multifocal leukoencephalopathy (PML); DMF-associated PML is very rare. It is unknown whether genetic predictors of SPL secondary to DMF treatment exist. We aimed to identify genetic predictors of reduced white blood cell (WBC) counts in DMF-treated multiple sclerosis (MS) patients. Genotyping (N = 1,258) and blood transcriptional profiling (N = 1,133) were performed on MS patients from DEFINE/CONFIRM. ALCs were categorized as: SPL, < 500 cells/µL for ≥6 months; moderate prolonged lymphopenia (MPL), < 800 cells/µL for ≥6 months, excluding SPL; mildly reduced lymphocytes, < 910 cells/µL at any point, excluding SPL and MPL; no lymphopenia, ≥910 cells/µL. Genome-wide association, HLA, and cross-sectional gene expression studies were performed. No common variants, HLA alleles, or expression profiles clinically useful for predicting SPL or MPL were identified. There was no overlap between genetic peaks and genetic loci known to be associated with WBC. Gene expression profiles were not associated with lymphopenia status. A classification model including gene expression features was not more predictive of lymphopenia status than standard covariates. There were no genetic predictors of SPL (or MPL) secondary to DMF treatment. Our results support ALC monitoring during DMF treatment as the most effective way to identify patients at risk of SPL.

Evolutionary history of disease-susceptibility loci identified in longitudinal exome-wide association studies

Background

Our longitudinal exome‐wide association studies previously detected various genetic determinants of complex disorders using ~26,000 single‐nucleotide polymorphisms (SNPs) that passed quality control and longitudinal medical examination data (mean follow‐up period, 5 years) in 4884–6022 Japanese subjects. We found that allele frequencies of several identified SNPs were remarkably different among four ethnic groups. Elucidating the evolutionary history of disease‐susceptibility loci may help us uncover the pathogenesis of the related complex disorders.

Methods

In the present study, we conducted evolutionary analyses such as extended haplotype homozygosity, focusing on genomic regions containing disease‐susceptibility loci and based on genotyping data of our previous studies and datasets from the 1000 Genomes Project.

Results

Our evolutionary analyses suggest that derived alleles of rs78338345 of GGA3, rs7656604 at 4q13.3, rs34902660 of SLC17A3, and six SNPs closely located at 12q24.1 associated with type 2 diabetes mellitus, obesity, dyslipidemia, and three complex disorders (hypertension, hyperuricemia, and dyslipidemia), respectively, rapidly expanded after the human dispersion from Africa (Out‐of‐Africa). Allele frequencies of GGA3 and six SNPs at 12q24.1 appeared to have remarkably changed in East Asians, whereas the derived alleles of rs34902660 of SLC17A3 and rs7656604 at 4q13.3 might have spread across Japanese and non‐Africans, respectively, although we cannot completely exclude the possibility that allele frequencies of disease‐associated loci may be affected by demographic events.

Conclusion

Our findings indicate that derived allele frequencies of nine disease‐associated SNPs (rs78338345 of GGA3, rs7656604 at 4q13.3, rs34902660 of SLC17A3, and six SNPs at 12q24.1) identified in the longitudinal exome‐wide association studies largely increased in non‐Africans after Out‐of‐Africa.

Genome-wide association studies of severe P. falciparum malaria susceptibility: progress, pitfalls and prospects

BACKGROUND

P. falciparum malaria has been recognized as one of the prominent evolutionary selective forces of human genome that led to the emergence of multiple host protective alleles. A comprehensive understanding of the genetic bases of severe malaria susceptibility and resistance can potentially pave ways to the development of new therapeutics and vaccines. Genome-wide association studies (GWASs) have recently been implemented in malaria endemic areas and identified a number of novel association genetic variants. However, there are several open questions around heritability, epistatic interactions, genetic correlations and associated molecular pathways among others. Here, we assess the progress and pitfalls of severe malaria susceptibility GWASs and discuss the biology of the novel variants.

RESULTS

We obtained all severe malaria susceptibility GWASs published thus far and accessed GWAS dataset of Gambian populations from European Phenome Genome Archive (EGA) through the MalariaGen consortium standard data access protocols. We noticed that, while some of the well-known variants including HbS and ABO blood group were replicated across endemic populations, only few novel variants were convincingly identified and their biological functions remain to be understood. We estimated SNP-heritability of severe malaria at 20.1% in Gambian populations and showed how advanced statistical genetic analytic methods can potentially be implemented in malaria susceptibility studies to provide useful functional insights.

CONCLUSIONS

The ultimate goal of malaria susceptibility study is to discover a novel causal biological pathway that provide protections against severe malaria; a fundamental step towards translational medicine such as development of vaccine and new therapeutics. Beyond singe locus analysis, the future direction of malaria susceptibility requires a paradigm shift from single -omics to multi-stage and multi-dimensional integrative functional studies that combines multiple data types from the human host, the parasite, the mosquitoes and the environment. The current biotechnological and statistical advances may eventually lead to the feasibility of systems biology studies and revolutionize malaria research.

Analysis of Hematological Traits in Polled Yak by Genome-Wide Association Studies Using Individual SNPs and Haplotypes

Yak (Bos grunniens) is an important domestic animal living in high-altitude plateaus. Due to inadequate disease prevention, each year, the yak industry suffers significant economic losses. The identification of causal genes that affect blood- and immunity-related cells could provide preliminary reference guidelines for the prevention of diseases in the population of yaks. The genome-wide association studies (GWASs) utilizing a single-marker or haplotype method were employed to analyze 15 hematological traits in the genome of 315 unrelated yaks. Single-marker GWASs identified a total of 43 significant SNPs, including 35 suggestive and eight genome-wide significant SNPs, associated with nine traits. Haplotype analysis detected nine significant haplotype blocks, including two genome-wide and seven suggestive blocks, associated with seven traits. The study provides data on the genetic variability of hematological traits in the yak. Five essential genes (GPLD1, EDNRA,APOB, HIST1H1E, and HIST1H2BI) were identified, which affect the HCT, HGB, RBC, PDW, PLT, and RDWSD traits and can serve as candidate genes for regulating hematological traits. The results provide a valuable reference to be used in the analysis of blood properties and immune diseases in the yak.

GWAS and PheWAS of red blood cell components in a Northern Nevadan cohort

In this study, we perform a full genome-wide association study (GWAS) to identify statistically significantly associated single nucleotide polymorphisms (SNPs) with three red blood cell (RBC) components and follow it with two independent PheWASs to examine associations between phenotypic data (case-control status of diagnoses or disease), significant SNPs, and RBC component levels. We first identified associations between the three RBC components: mean platelet volume (MPV), mean corpuscular volume (MCV), and platelet counts (PC), and the genotypes of approximately 500,000 SNPs on the Illumina Infimum DNA Human OmniExpress-24 BeadChip using a single cohort of 4,673 Northern Nevadans. Twenty-one SNPs in five major genomic regions were found to be statistically significantly associated with MPV, two regions with MCV, and one region with PC, with p<5x10^-8. Twenty-nine SNPs and nine chromosomal regions were identified in 30 previous GWASs, with effect sizes of similar magnitude and direction as found in our cohort. The two strongest associations were SNP rs1354034 with MPV (p = 2.4x10^-13) and rs855791 with MCV (p = 5.2x10^-12). We then examined possible associations between these significant SNPs and incidence of 1,488 phenotype groups mapped from International Classification of Disease version 9 and 10 (ICD9 and ICD10) codes collected in the extensive electronic health record (EHR) database associated with Healthy Nevada Project consented participants. Further leveraging data collected in the EHR, we performed an additional PheWAS to identify associations between continuous red blood cell (RBC) component measures and incidence of specific diagnoses. The first PheWAS illuminated whether SNPs associated with RBC components in our cohort were linked with other hematologic phenotypic diagnoses or diagnoses of other nature. Although no SNPs from our GWAS were identified as strongly associated to other phenotypic components, a number of associations were identified with p-values ranging between 1x10^-3 and 1x10^-4 with traits such as respiratory failure, sleep disorders, hypoglycemia, hyperglyceridemia, GERD and IBS. The second PheWAS examined possible phenotypic predictors of abnormal RBC component measures: a number of hematologic phenotypes such as thrombocytopenia, anemias, hemoglobinopathies and pancytopenia were found to be strongly associated to RBC component measures; additional phenotypes such as (morbid) obesity, malaise and fatigue, alcoholism, and cirrhosis were also identified to be possible predictors of RBC component measures.

Natural history of benign ethnic neutropenia in individuals of African ancestry

BACKGROUND

Benign ethnic neutropenia (BEN), defined by neutrophil count <1.5 k/μL in the absence of other causes, is an asymptomatic condition more commonly observed in individuals of African ancestry. However, the natural history of this condition has been less well described.

METHODS

Individuals with BEN were retrospectively identified by chart review or referral to hematology clinics. They were then invited to enroll in a prospective natural history study. Retrospective and prospective clinical and laboratory data were combined for descriptive analyses.

FINDINGS

46 participants, younger and older adults from 2 institutions, had BEN. Hypertension was reported in 30%, musculoskeletal disorders in 15%, and upper respiratory infection in 33% of these adults. Their leukopenia resulted from isolated neutropenia, ranging from 1000 and 1500 cells/μL. The severity of infections was mild and the frequency was similar to other healthy individuals in the ambulatory clinic.

INTERPRETATION

In this group of BEN participants, their leukopenia was stable over time, and they had low rates of infections or common medical disorders, confirming the benign nature of this condition. The presence of BEN in children, younger adults, and older adults suggest a hereditary pattern for BEN.

A decade in psychiatric GWAS research

After more than 10 years of accumulated efforts, genome-wide association studies (GWAS) have led to many findings, most of which have been deposited into the GWAS Catalog. Between GWAS's inception and March 2017, the GWAS Catalog has collected 2429 studies, 1818 phenotypes, and 28,462 associated SNPs. We reclassified the psychology-related phenotypes into 217 reclassified phenotypes, which accounted for 514 studies and 7052 SNPs. In total, 1223 of the SNPs reached genome-wide significance. Of these, 147 were replicated for the same psychological trait in different studies. Another 305 SNPs were replicated within one original study. The SNPs rs2075650 and rs4420638 were linked to the most replications within a single reclassified phenotype or very similar reclassified phenotypes; both were associated with Alzheimer's disease (AD). Schizophrenia was associated with 74 within-phenotype SNPs reported in independents studies. Alzheimer's disease and schizophrenia were both linked to some physical phenotypes, including cholesterol and body mass index, through common GWAS signals. Alzheimer's disease also shared risk SNPs with age-related phenotypes such as age-related macular degeneration and longevity. Smoking-related SNPs were linked to lung cancer and respiratory function. Alcohol-related SNPs were associated with cardiovascular and digestive system phenotypes and disorders. Two separate studies also identified a shared risk SNP for bipolar disorder and educational attainment. This review revealed a list of reproducible SNPs worthy of future functional investigation. Additionally, by identifying SNPs associated with multiple phenotypes, we illustrated the importance of studying the relationships among phenotypes to resolve the nature of their causal links. The insights within this review will hopefully pave the way for future evidence-based genetic studies.

Genetic signatures of gene flow and malaria-driven natural selection in sub-Saharan populations of the "endemic Burkitt Lymphoma belt"

Populations in sub-Saharan Africa have historically been exposed to intense selection from chronic infection with falciparum malaria. Interestingly, populations with the highest malaria intensity can be identified by the increased occurrence of endemic Burkitt Lymphoma (eBL), a pediatric cancer that affects populations with intense malaria exposure, in the so called "eBL belt" in sub-Saharan Africa. However, the effects of intense malaria exposure and sub-Saharan populations' genetic histories remain poorly explored. To determine if historical migrations and intense malaria exposure have shaped the genetic composition of the eBL belt populations, we genotyped ~4.3 million SNPs in 1,708 individuals from Ghana and Northern Uganda, located on opposite sides of eBL belt and with ≥ 7 months/year of intense malaria exposure and published evidence of high incidence of BL. Among 35 Ghanaian tribes, we showed a predominantly West-Central African ancestry and genomic footprints of gene flow from Gambian and East African populations. In Uganda, the North West population showed a predominantly Nilotic ancestry, and the North Central population was a mixture of Nilotic and Southern Bantu ancestry, while the Southwest Ugandan population showed a predominant Southern Bantu ancestry. Our results support the hypothesis of diverse ancestral origins of the Ugandan, Kenyan and Tanzanian Great Lakes African populations, reflecting a confluence of Nilotic, Cushitic and Bantu migrations in the last 3000 years. Natural selection analyses suggest, for the first time, a strong positive selection signal in the ATP2B4 gene (rs10900588) in Northern Ugandan populations. These findings provide important baseline genomic data to facilitate disease association studies, including of eBL, in eBL belt populations.

Genetic associations of hemoglobin in children with chronic kidney disease in the PediGFR Consortium

BACKGROUND

Genome-wide association studies (GWAS) in healthy populations have identified variants associated with erythrocyte traits, but genetic causes of hemoglobin variation in children with CKD are incompletely understood.

METHODS

The Pediatric Investigation of Genetic Factors Linked with Renal Progression (PediGFR) Consortium comprises three pediatric CKD cohorts: Chronic Kidney Disease in Children (CKiD), Effect of Strict Blood Pressure Control and ACE Inhibition on the Progression of CRF in Pediatric Patients (ESCAPE), and Cardiovascular Comorbidity in Children with CKD (4C). We performed cross-sectional and longitudinal association studies of single-nucleotide polymorphisms (SNPs) in 1125 patients.

RESULTS

Children of European (n = 725) or Turkish (n = 400) ancestry (EA or TA) were included. In cross-sectional analysis, two SNPs (rs10758658 and rs12718597) previously associated with RBC traits were significantly associated with hemoglobin levels in children of EA and TA. In longitudinal analysis, SNP rs2540917 was nominally associated with hemoglobin in EA and TA children.

CONCLUSIONS

SNPs associated with erythrocyte traits in healthy populations were marginally significant for an association with hemoglobin. Further analyses/replication studies are needed in larger CKD cohorts to investigate SNPs of unknown significance associated with hemoglobin. Functional studies will be required to confirm that the observed associations between SNPs and clinical phenotype are causal.

Genetics of platelet traits in ischaemic stroke: focus on mean platelet volume and platelet count

Purpose/Aim of the study: The aim of this review is to summarize the role of genetic variants affecting mean platelet volume (MPV) and platelet count (PLT) leading to higher platelet reactivity and in turn to thrombotic events like stroke and cardiovascular diseases.

MATERIALS AND METHODS

A search was conducted in PUBMED, MEDLINE, EMBASE, PROQUEST, Science Direct, Cochrane Library, and Google Scholar related to the studies focussing on genome-wide association studies (GWAS), whole exome sequencing (WES), whole genome sequencing (WGS), phenome-wide association studies (PheWAS) and multi-omic analysis that have been employed to identify the genetic variants influencing MPV and PLT.

RESULTS

Antiplatelet agents underscore the crucial role of platelets in the pathogenesis of stroke. Higher platelet reactivity in terms of mean platelet volume (MPV) and platelet count (PLT) contributes significantly to the interindividual variation in platelet reaction at the site of vessel wall injury. Some individuals encounter thrombotic events as platelets get occluded at the site of vessel wall injury whereas others heal the injury without occluding the circulation. Evidence suggests that MPV and PLT have a strong genetic component. High throughput techniques including genome-wide association studies (GWAS), whole exome sequencing (WES), whole genome sequencing (WGS), phenome-wide association studies (PheWAS) and multi-omic analysis have identified different genetic variants influencing MPV and PLT.

CONCLUSIONS

Identification of complex genetic cross talks affecting PLT and MPV might help to develop novel treatment strategies in treating neurovascular diseases like stroke.

Unfolding of hidden white blood cell count phenotypes for gene discovery using latent class mixed modeling

Resting-state white blood cell (WBC) count is a marker of inflammation and immune system health. There is evidence that WBC count is not fixed over time and there is heterogeneity in WBC trajectory that is associated with morbidity and mortality. Latent class mixed modeling (LCMM) is a method that can identify unobserved heterogeneity in longitudinal data and attempts to classify individuals into groups based on a linear model of repeated measurements. We applied LCMM to repeated WBC count measures derived from electronic medical records of participants of the National Human Genetics Research Institute (NHRGI) electronic MEdical Record and GEnomics (eMERGE) network study, revealing two WBC count trajectory phenotypes. Advancing these phenotypes to GWAS, we found genetic associations between trajectory class membership and regions on chromosome 1p34.3 and chromosome 11q13.4. The chromosome 1 region contains CSF3R, which encodes the granulocyte colony-stimulating factor receptor. This protein is a major factor in neutrophil stimulation and proliferation. The association on chromosome 11 contain genes RNF169 and XRRA1; both involved in the regulation of double-strand break DNA repair.

Identification of nine novel loci related to hematological traits in a Japanese population

Recent genome-wide association studies have identified various genetic variants associated with hematological traits. Although it is possible that quantitative data of hematological traits are varied among the years examined, conventional genome-wide association studies have been conducted in a cross-sectional manner that measures traits at a single point in time. To address this issue, we have traced blood profiles in 4,884 Japanese individuals who underwent annual health check-ups for several years. In the present study, longitudinal exome-wide association studies were conducted to identify genetic variants related to 13 hematological phenotypes. The generalized estimating equation model showed that a total of 67 single nucleotide polymorphisms (SNPs) were significantly [false discovery rate (FDR) of <0.01] associated with hematological phenotypes. Of the 67 SNPs, nine SNPs were identified as novel hematological markers: rs4686683 of SENP2 for red blood cell count (FDR = 0.008, P = 5.5 × 10⁻⁶); rs3917688 of SELP for mean corpuscular volume (FDR = 0.005, P = 2.4 × 10⁻⁶); rs3133745 of C8orf37-AS1 for white blood cell count (FDR = 0.003, P = 1.3 × 10⁻⁶); rs13121954 at 4q31.2 for basophil count (FDR = 0.007, P = 3.1 × 10⁻⁵); rs7584099 at 2q22.3 (FDR = 2.6 × 10⁻⁵, P = 8.8 × 10⁻⁸), rs1579219 of HCG17 (FDR = 0.003, P = 2.0 × 10⁻⁵), and rs10757049 of DENND4C (FDR = 0.008, P = 5.6 × 10⁻⁵) for eosinophil count; rs12338 of CTSB for neutrophil count (FDR = 0.007, P = 2.9 × 10⁻⁵); and rs395967 of OSMR-AS1 for monocyte count (FDR = 0.008, P = 3.2 × 10⁻⁵).

Generalization and fine mapping of red blood cell trait genetic associations to multi-ethnic populations: The PAGE Study

Red blood cell (RBC) traits provide insight into a wide range of physiological states and exhibit moderate to high heritability, making them excellent candidates for genetic studies to inform underlying biologic mechanisms. Previous RBC trait genome-wide association studies were performed primarily in European- or Asian-ancestry populations, missing opportunities to inform understanding of RBC genetic architecture in diverse populations and reduce intervals surrounding putative functional SNPs through fine-mapping. Here, we report the first fine-mapping of six correlated (Pearson's r range: |0.04 - 0.92|) RBC traits in up to 19,036 African Americans and 19,562 Hispanic/Latinos participants of the Population Architecture using Genomics and Epidemiology (PAGE) consortium. Trans-ethnic meta-analysis of race/ethnic- and study-specific estimates for approximately 11,000 SNPs flanking 13 previously identified association signals as well as 150,000 additional array-wide SNPs was performed using inverse-variance meta-analysis after adjusting for study and clinical covariates. Approximately half of previously reported index SNP-RBC trait associations generalized to the trans-ethnic study population (p<1.7x10-4 ); previously unreported independent association signals within the ABO region reinforce the potential for multiple functional variants affecting the same locus. Trans-ethnic fine-mapping did not reveal additional signals at the HFE locus independent of the known functional variants. Finally, we identified a potential novel association in the Hispanic/Latino study population at the HECTD4/RPL6 locus for RBC count (p=1.9x10-7 ). The identification of a previously unknown association, generalization of a large proportion of known association signals, and refinement of known association signals all exemplify the benefits of genetic studies in diverse populations. This article is protected by copyright. All rights reserved.

PheWAS and Beyond: The Landscape of Associations with Medical Diagnoses and Clinical Measures across 38,662 Individuals from Geisinger

Most phenome-wide association studies (PheWASs) to date have used a small to moderate number of SNPs for association with phenotypic data. We performed a large-scale single-cohort PheWAS, using electronic health record (EHR)-derived case-control status for 541 diagnoses using International Classification of Disease version 9 (ICD-9) codes and 25 median clinical laboratory measures. We calculated associations between these diagnoses and traits with ∼630,000 common frequency SNPs with minor allele frequency > 0.01 for 38,662 individuals. In this landscape PheWAS, we explored results within diseases and traits, comparing results to those previously reported in genome-wide association studies (GWASs), as well as previously published PheWASs. We further leveraged the context of functional impact from protein-coding to regulatory regions, providing a deeper interpretation of these associations. The comprehensive nature of this PheWAS allows for novel hypothesis generation, the identification of phenotypes for further study for future phenotypic algorithm development, and identification of cross-phenotype associations.

A robust and powerful two-step testing procedure for local ancestry adjusted allelic association analysis in admixed populations

Genetic association studies in admixed populations allow us to gain deeper understanding of the genetic architecture of human diseases and traits. However, population stratification, complicated linkage disequilibrium (LD) patterns, and the complex interplay of allelic and ancestry effects on phenotypic traits pose challenges in such analyses. These issues may lead to detecting spurious associations and/or result in reduced statistical power. Fortunately, if handled appropriately, these same challenges provide unique opportunities for gene mapping. To address these challenges and to take these opportunities, we propose a robust and powerful two-step testing procedure Local Ancestry Adjusted Allelic (LAAA) association. In the first step, LAAA robustly captures associations due to allelic effect, ancestry effect, and interaction effect, allowing detection of effect heterogeneity across ancestral populations. In the second step, LAAA identifies the source of association, namely allelic, ancestry, or the combination. By jointly modeling allele, local ancestry, and ancestry-specific allelic effects, LAAA is highly powerful in capturing the presence of interaction between ancestry and allele effect. We evaluated the validity and statistical power of LAAA through simulations over a broad spectrum of scenarios. We further illustrated its usefulness by application to the Candidate Gene Association Resource (CARe) African American participants for association with hemoglobin levels. We were able to replicate independent groups' previously identified loci that would have been missed in CARe without joint testing. Moreover, the loci, for which LAAA detected potential effect heterogeneity, were replicated among African Americans from the Women's Health Initiative study. LAAA is freely available at https://yunliweb.its.unc.edu/LAAA.

Analyses of genome wide association data, cytokines, and gene expression in African-Americans with benign ethnic neutropenia

BACKGROUND

Benign ethnic neutropenia (BEN) is a hematologic condition associated with people of African ancestry and specific Middle Eastern ethnic groups. Prior genetic association studies in large population showed that rs2814778 in Duffy Antigen Receptor for Chemokines (DARC) gene, specifically DARC null red cell phenotype, was associated with BEN. However, the mechanism of this red cell phenotype leading to low white cell count remained elusive.

METHODS

We conducted an extreme phenotype design genome-wide association study (GWAS), analyzed ~16 million single nucleotide polymorphisms (SNP) in 1,178 African-Americans individuals from the Reasons for Geographic and Racial Differences in Stroke (REGARDS) study and replicated from 819 African-American participants in the Atherosclerosis Risk in Communities (ARIC) study. Conditional analyses on rs2814778 were performed to identify additional association signals on chromosome 1q22. In a separate cohort of healthy individuals with and without BEN, whole genome gene expression from peripheral blood neutrophils were analyzed for DARC.

RESULTS

We confirmed that rs2814778 in DARC was associated with BEN (p = 4.09×10-53). Conditioning on rs2814778 abolished other significant chromosome 1 associations. Inflammatory cytokines (IL-2, 6, and 10) in participants in the Howard University Family Study (HUFS) and Multi-Ethnic Study in Atherosclerosis (MESA) showed similar levels in individuals homozygous for the rs2814778 allele compared to others, indicating cytokine sink hypothesis played a minor role in leukocyte homeostasis. Gene expression in neutrophils of individuals with and without BEN was also similar except for low DARC expression in BEN, suggesting normal function. BEN neutrophils had slightly activated profiles in leukocyte migration and hematopoietic stem cell mobilization pathways (expression fold change <2).

CONCLUSIONS

These results in humans support the notion of DARC null erythroid progenitors preferentially differentiating to myeloid cells, leading to activated DARC null neutrophils egressing from circulation to the spleen, and causing relative neutropenia. Collectively, these human data sufficiently explained the mechanism DARC null red cell phenotype causing BEN and further provided a biologic basis that BEN is clinically benign.

Toward the Relevance of Platelet Subpopulations for Transfusion Medicine

Circulating platelets consist of subpopulations with different age, maturation state and size. In this review, we address the association between platelet size and platelet function and summarize the current knowledge on platelet subpopulations including reticulated platelets, procoagulant platelets and platelets exposing signals to mediate their clearance. Thereby, we emphasize the impact of platelet turnover as an important condition for platelet production in vivo. Understanding of the features that characterize platelet subpopulations is very relevant for the methods of platelet concentrate production, which may enrich or deplete particular platelet subpopulations. Moreover, the concept of platelet size being associated with platelet function may be attractive for transfusion medicine as it holds the perspective to separate platelet subpopulations with specific functional capabilities.

2SNP heritability and effects of genetic variants for neutrophil-to-lymphocyte and platelet-to-lymphocyte ratio

Neutrophil-to-lymphocyte ratio (NLR) and platelet-to-lymphocyte ratio (PLR) are important biomarkers for disease development and progression. To gain insight into the genetic causes of variance in NLR and PLR in the general population, we conducted genome-wide association (GWA) analyses and estimated SNP heritability in a sample of 5901 related healthy Dutch individuals. GWA analyses identified a new genome-wide significant locus on the HBS1L-MYB intergenic region for PLR, which replicated in a sample of 2538 British twins. For platelet count, we replicated three known genome-wide significant loci in our cohort (at CCDC71L-PIK3CG, BAK1 and ARHGEF3). For neutrophil count, we replicated the PSMD3 locus. For the identified top SNPs, we found significant cis and trans expression quantitative trait loci effects for several loci involved in hematological and immunological pathways. Linkage Disequilibrium score (LD) regression analyses for PLR and NLR confirmed that both traits are heritable, with a polygenetic SNP heritability for PLR of 14.1%, and for NLR of 2.4%. Genetic correlations were present between ratios and the constituent counts, with the genetic correlation (r=0.45) of PLR with platelet count reaching statistical significance. In conclusion, we established that two important biomarkers have a significant heritable SNP component, and identified the first genome-wide locus for PLR.

An erythroid-specific ATP2B4 enhancer mediates red blood cell hydration and malaria susceptibility

The lack of mechanistic explanations for many genotype-phenotype associations identified by GWAS precludes thorough assessment of their impact on human health. Here, we conducted an expression quantitative trait locus (eQTL) mapping analysis in erythroblasts and found erythroid-specific eQTLs for ATP2B4, the main calcium ATPase of red blood cells (rbc). The same SNPs were previously associated with mean corpuscular hemoglobin concentration (MCHC) and susceptibility to severe malaria infection. We showed that Atp2b4-/- mice demonstrate increased MCHC, confirming ATP2B4 as the causal gene at this GWAS locus. Using CRISPR-Cas9, we fine mapped the genetic signal to an erythroid-specific enhancer of ATP2B4. Erythroid cells with a deletion of the ATP2B4 enhancer had abnormally high intracellular calcium levels. These results illustrate the power of combined transcriptomic, epigenomic, and genome-editing approaches in characterizing noncoding regulatory elements in phenotype-relevant cells. Our study supports ATP2B4 as a potential target for modulating rbc hydration in erythroid disorders and malaria infection.

The genetics of platelet count and volume in humans

The last decade has witnessed an explosion in the depth, variety, and amount of human genetic data that can be generated. This revolution in technical and analytical capacities has enabled the genetic investigation of human traits and disease in thousands to now millions of participants. Investigators have taken advantage of these advancements to gain insight into platelet biology and the platelet's role in human disease. To do so, large human genetics studies have examined the association of genetic variation with two quantitative traits measured in many population and patient based cohorts: platelet count (PLT) and mean platelet volume (MPV). This article will review the many human genetic strategies-ranging from genome-wide association study (GWAS), Exomechip, whole exome sequencing (WES), to whole genome sequencing (WGS)-employed to identify genes and variants that contribute to platelet traits. Additionally, we will discuss how these investigations have examined and interpreted the functional implications of these newly identified genetic factors and whether they also impart risk to human disease. The depth and size of genetic, phenotypic, and other -omic data are primed to continue their growth in the coming years and provide unprecedented opportunities to gain critical insights into platelet biology and how platelets contribute to disease.

SNP in human ARHGEF3 promoter is associated with DNase hypersensitivity, transcript level and platelet function, and Arhgef3 KO mice have increased mean platelet volume

Genome-wide association studies have identified a genetic variant at 3p14.3 (SNP rs1354034) that strongly associates with platelet number and mean platelet volume in humans. While originally proposed to be intronic, analysis of mRNA expression in primary human hematopoietic subpopulations reveals that this SNP is located directly upstream of the predominantly expressed ARHGEF3 isoform in megakaryocytes (MK). We found that ARHGEF3, which encodes a Rho guanine exchange factor, is dramatically upregulated during both human and murine MK maturation. We show that the SNP (rs1354034) is located in a DNase I hypersensitive region in human MKs and is an expression quantitative locus (eQTL) associated with ARHGEF3 expression level in human platelets, suggesting that it may be the causal SNP that accounts for the variations observed in human platelet traits and ARHGEF3 expression. In vitro human platelet activation assays revealed that rs1354034 is highly correlated with human platelet activation by ADP. In order to test whether ARHGEF3 plays a role in MK development and/or platelet function, we developed an Arhgef3 KO/LacZ reporter mouse model. Reflecting changes in gene expression, LacZ expression increases during MK maturation in these mice. Although Arhgef3 KO mice have significantly larger platelets, loss of Arhgef3 does not affect baseline MK or platelets nor does it affect platelet function or platelet recovery in response to antibody-mediated platelet depletion compared to littermate controls. In summary, our data suggest that modulation of ARHGEF3 gene expression in humans with a promoter-localized SNP plays a role in human MKs and human platelet function—a finding resulting from the biological follow-up of human genetic studies. Arhgef3 KO mice partially recapitulate the human phenotype.

Genome-wide association study of red blood cell traits in Hispanics/Latinos: The Hispanic Community Health Study/Study of Latinos

Prior GWAS have identified loci associated with red blood cell (RBC) traits in populations of European, African, and Asian ancestry. These studies have not included individuals with an Amerindian ancestral background, such as Hispanics/Latinos, nor evaluated the full spectrum of genomic variation beyond single nucleotide variants. Using a custom genotyping array enriched for Amerindian ancestral content and 1000 Genomes imputation, we performed GWAS in 12,502 participants of Hispanic Community Health Study and Study of Latinos (HCHS/SOL) for hematocrit, hemoglobin, RBC count, RBC distribution width (RDW), and RBC indices. Approximately 60% of previously reported RBC trait loci generalized to HCHS/SOL Hispanics/Latinos, including African ancestral alpha- and beta-globin gene variants. In addition to the known 3.8kb alpha-globin copy number variant, we identified an Amerindian ancestral association in an alpha-globin regulatory region on chromosome 16p13.3 for mean corpuscular volume and mean corpuscular hemoglobin. We also discovered and replicated three genome-wide significant variants in previously unreported loci for RDW (SLC12A2 rs17764730, PSMB5 rs941718), and hematocrit (PROX1 rs3754140). Among the proxy variants at the SLC12A2 locus we identified rs3812049, located in a bi-directional promoter between SLC12A2 (which encodes a red cell membrane ion-transport protein) and an upstream anti-sense long-noncoding RNA, LINC01184, as the likely causal variant. We further demonstrate that disruption of the regulatory element harboring rs3812049 affects transcription of SLC12A2 and LINC01184 in human erythroid progenitor cells. Together, these results reinforce the importance of genetic study of diverse ancestral populations, in particular Hispanics/Latinos.

Integrated genome-wide analysis of expression quantitative trait loci aids interpretation of genomic association studies

Background

Identification of single nucleotide polymorphisms (SNPs) associated with gene expression levels, known as expression quantitative trait loci (eQTLs), may improve understanding of the functional role of phenotype-associated SNPs in genome-wide association studies (GWAS). The small sample sizes of some previous eQTL studies have limited their statistical power. We conducted an eQTL investigation of microarray-based gene and exon expression levels in whole blood in a cohort of 5257 individuals, exceeding the single cohort size of previous studies by more than a factor of 2.

Results

We detected over 19,000 independent lead cis-eQTLs and over 6000 independent lead trans-eQTLs, targeting over 10,000 gene targets (eGenes), with a false discovery rate (FDR) < 5%. Of previously published significant GWAS SNPs, 48% are identified to be significant eQTLs in our study. Some trans-eQTLs point toward novel mechanistic explanations for the association of the SNP with the GWAS-related phenotype. We also identify 59 distinct blocks or clusters of trans-eQTLs, each targeting the expression of sets of six to 229 distinct trans-eGenes. Ten of these sets of target genes are significantly enriched for microRNA targets (FDR < 5%). Many of these clusters are associated in GWAS with multiple phenotypes.

Conclusions

These findings provide insights into the molecular regulatory patterns involved in human physiology and pathophysiology. We illustrate the value of our eQTL database in the context of a recent GWAS meta-analysis of coronary artery disease and provide a list of targeted eGenes for 21 of 58 GWAS loci.

Electronic supplementary material

The online version of this article (doi:10.1186/s13059-016-1142-6) contains supplementary material, which is available to authorized users.

Platelet to lymphocyte ratio in the prediction of adverse outcomes after acute coronary syndrome: a meta-analysis

Recent studies have shown platelet to lymphocyte ratio (PLR) to be a potential inflammatory marker in cardiovascular diseases. We performed a meta-analysis to systematically evaluate the prognostic role of PLR in acute coronary syndrome (ACS). A comprehensive literature search up to May 18, 2016 was conducted from PUBMED, EMBASE and Web of science to identify related studies. The risk ratio (RR) with 95% confidence interval (CI) was extracted or calculated for effect estimates. Totally ten studies involving 8932 patients diagnosed with ACS were included in our research. We demonstrated that patients with higher PLR level had significantly higher risk of in-hospital adverse outcomes (RR = 2.24, 95%CI = 1.81–2.77) and long-term adverse outcomes (RR = 2.32, 95%CI = 1.64–3.28). Sensitivity analyses confirmed the stability of our results. We didn’t detect significant publication bias by Begg’s and Egger’s test (p > 0.05). In conclusion, our meta-analysis revealed that PLR is promising biomarker in predicting worse prognosis in ACS patients. The results should be validated by future large-scale, standard investigations.

Genome-wide Trans-ethnic Meta-analysis Identifies Seven Genetic Loci Influencing Erythrocyte Traits and a Role for RBPMS in Erythropoiesis

Genome-wide association studies (GWASs) have identified loci for erythrocyte traits in primarily European ancestry populations. We conducted GWAS meta-analyses of six erythrocyte traits in 71,638 individuals from European, East Asian, and African ancestries using a Bayesian approach to account for heterogeneity in allelic effects and variation in the structure of linkage disequilibrium between ethnicities. We identified seven loci for erythrocyte traits including a locus (RBPMS/GTF2E2) associated with mean corpuscular hemoglobin and mean corpuscular volume. Statistical fine-mapping at this locus pointed to RBPMS at this locus and excluded nearby GTF2E2. Using zebrafish morpholino to evaluate loss of function, we observed a strong in vivo erythropoietic effect for RBPMS but not for GTF2E2, supporting the statistical fine-mapping at this locus and demonstrating that RBPMS is a regulator of erythropoiesis. Our findings show the utility of trans-ethnic GWASs for discovery and characterization of genetic loci influencing hematologic traits.