Genome-wide association study of white blood cell count in 16,388 African Americans: the continental origins and genetic epidemiology network (COGENT)

The C11orf24 Gene as a Useful Biomarker for Predicting Severe Neutropenia in Modified FOLFIRINOX for Pancreatic Cancer

Pancreatic cancer (PC) is an aggressive and lethal tumor with a poor prognosis. FOLFIRINOX improves the prognosis of patients with PC; however, despite UGT1A1 screening, adverse events, such as severe neutropenia, occur frequently. This study aimed to identify the novel biomarkers of severe neutropenia in patients treated with modified FOLFIRINOX (mFFX) for PC. In this study, patients with PC treated with mFFX (n = 71) and gemcitabine plus nab-paclitaxel (GnP) (n = 92) and patients with colorectal cancer treated with FOLFOXIRI (n = 50) were included. Genome-wide screening using whole-exome sequencing was performed during the screening phase. Validation analysis was performed using polymerase chain reaction genotyping, the Cochran-Armitage trend test, and multivariate analysis. The diagnostic performance of combined risk factors for severe neutropenia was examined using logistic regression with leave-one-out cross-validation. Three gene polymorphisms were selected from the screening phase and subjected to the validation phase. In the validation phase, a single nucleotide polymorphism in C11orf24 (c.448C>T, rs901827) was significantly correlated with ≥ Grade 3 neutropenia in mFFX and FOLFOXIRI but not in GnP. Multivariate analysis showed C11orf24 and baseline neutrophil count as independent risk factors for ≥ Grade 3 neutropenia. The diagnostic performance of the neutropenia prediction model showed areas under the curve of 0.754 (sensitivity = 0.605, specificity = 0.848) and 0.856 (sensitivity = 0.800, specificity = 0.893) for ≥ Grade 3 and 4 neutropenia, respectively. The C11orf24 gene and baseline neutrophil count may be useful biomarkers for predicting severe neutropenia following irinotecan-containing triplet chemotherapy.

Characterization of the phenotypic consequences of the Duffy-null genotype

ABSTRACT

A wealth of research focused on African American populations has connected rs2814778-CC ("Duffy-null") to decreased neutrophil (neutropenia) and leukocyte counts (leukopenia). Although it has been proposed that this variant is benign, prior studies have shown that the misinterpretation of Duffy-null-associated neutropenia and leukopenia can lead to unnecessary bone marrow biopsies, inequities in cytotoxic and chemotherapeutic treatment courses, underenrollment in clinical trials, and other disparities. To investigate the phenotypic correlates of Duffy-null status, we conducted a phenome-wide association study across >1400 clinical conditions in All of Us, the Vanderbilt University Medical Center's Biobank, and the Million Veteran Program. This reveals that Duffy-null status is only reproducibly associated with changes in white blood cell count and not with any disease outcomes. Moreover, we find that Duffy-null-associated neutropenia is on average less severe than other neutropenia cases in All of Us. We also show that this genotype is present in considerable frequencies in All of Us populations that are genetically similar to African (68%) and Middle Eastern (14%) 1000 Genomes/Human Genome Diversity Project reference populations as well as those who identify with >1 race (12%), as Pacific Islander (7%), and as Hispanic (5%). Furthermore, we find that race is not an accurate predictor of Duffy-null status or associated benign neutropenia. Our research suggests that broad genetic screening of rs2814778 across all populations could provide a more robust and accurate understanding of white blood cell count and mitigate resulting health disparities.

Cross-cohort analysis of expression and splicing quantitative trait loci in TOPMed

Most genetic variants associated with complex traits and diseases occur in non-coding genomic regions and are hypothesized to regulate gene expression. To understand the genetics underlying gene expression variability, we characterize 14,324 ancestrally diverse RNA-sequencing samples from the NHLBI Trans-Omics for Precision Medicine (TOPMed) program and integrate whole genome sequencing data to perform cis and trans expression and splicing quantitative trait locus (cis-/trans-e/sQTL) analyses in six tissues and cell types, most notably whole blood (N=6,454) and lung (N=1,291). We show this dataset enables greater detection of secondary cis-e/sQTL signals than was achieved in previous studies, and that secondary cis-eQTL and primary trans-eQTL signal discovery is not saturated even though eGene discovery is. Most TOPMed trans-eQTL signals colocalize with cis-e/sQTL signals, suggesting many trans signals are mediated by cis signals. We fine-map European UK BioBank GWAS signals from 164 traits and colocalize the resulting 34,107 fine-mapped GWAS signals with TOPMed e/sQTL signals, finding that of 10,611 GWAS signals with a colocalization, 7,096 GWAS signals colocalize with at least one secondary e/sQTL signal. These results demonstrate that larger e/sQTL analyses will continue to uncover secondary e/sQTL signals, and that these new signals will benefit GWAS interpretation.

Genetic variants associated with white blood cell count amongst individuals with sickle cell disease

BACKGROUND

Sickle cell disease (SCD) is a Mendelian disorder characterized by a point mutation in the β-globin gene that leads to sickling of erythrocytes. Several studies have shown that absolute neutrophil count is strongly associated with clinical severity of SCD, suggesting an apparent role of white blood cells (WBC) in SCD pathology. However, the mechanism by which genetic variants lead to WBC count differences in SCD patients remains unclear.

METHODS

Genome-wide association (GWA) analyses were carried out amongst a cohort of 2409 Brazil SCD participants. Association of WBC count and genetic markers were investigated in homozygous sickle cell anaemia participants and compound heterozygous sickle cell haemoglobin C participants.

RESULTS

GWA analysis showed that variants in genes TERT, ACKR1, and FAM3C are associated with WBC count variation. The well-studied association between WBC count and Duffy null phenotype (variant in ACKR1) in healthy populations was replicated, reinforcing the influence of the SNP rs2814778 (T>C) in WBC count.

CONCLUSION

Genetics plays an important role in regulating WBC count in patients with SCD. Our results point to possible mechanisms involved in WBC count variation and as increased WBC count is associated with more severe SCD, these results could suggest potential therapeutic targets for individuals with SCD.

Lower Neutrophil Count Without Clinical Consequence Among Children of African Ancestry Living With HIV in Canada

OBJECTIVE

To investigate the association between African ancestry and neutrophil counts among children living with HIV (CLWH). We also examined whether medications, clinical conditions, hospitalization, or HIV virologic control were associated with low neutrophil counts or African ancestry.

DESIGN

We conducted a secondary analysis of the Early Pediatric Initiation Canada Child Cure Cohort (EPIC4) Study, a multicenter prospective cohort study of CLWH across 8 Canadian pediatric HIV care centers.

METHODS

We classified CLWH according to African ancestry, defined as "African," "Caribbean," or "Black" maternal race. Longitudinal laboratory data (white blood cells, neutrophils, lymphocytes, viral load, and CD4 count) and clinical data (hospitalizations, AIDS-defining conditions, and treatments) were abstracted from medical records.

RESULTS

Among 217 CLWH (median age 14, 55% female), 145 were of African ancestry and 72 were of non-African ancestry. African ancestry was associated with lower neutrophil counts, white blood cell counts, and neutrophil-lymphocyte ratios. Neutrophil count <1.5 × 109/L was detected in 60% of CLWH of African ancestry, compared with 31% of CLWH of non-African ancestry (P < 0.0001), representing a 2.0-fold higher relative frequency (95% CI: 1.4-2.9). Neutrophil count was on average 0.74 × 109/L (95% CI: 0.45 to 1.0) lower in CLWH of African ancestry (P < 0.0001). Neither neutrophil count<1.5 × 109/L nor African ancestry was associated with medications, hospitalizations, AIDS-defining conditions, or markers of virologic control (viral load, sustained viral suppression, and lifetime nadir CD4).

CONCLUSIONS

In CLWH, African ancestry is associated with lower neutrophil counts, without clinical consequences. A flexible evaluation of neutrophil counts in CLWH of African ancestry may avoid unnecessary interventions.

Discerning clinicopathological features of congenital neutropenia syndromes: an approach to diagnostically challenging differential diagnoses

The congenital neutropenia syndromes are rare haematological conditions defined by impaired myeloid precursor differentiation or function. Patients are prone to severe infections with high mortality rates in early life. While some patients benefit from granulocyte colony-stimulating factor treatment, they may still face an increased risk of bone marrow failure, myelodysplastic syndrome and acute leukaemia. Accurate diagnosis is crucial for improved outcomes; however, diagnosis depends on familiarity with a heterogeneous group of rare disorders that remain incompletely characterised. The clinical and pathological overlap between reactive conditions, primary and congenital neutropenias, bone marrow failure, and myelodysplastic syndromes further clouds diagnostic clarity.We review the diagnostically useful clinicopathological and morphological features of reactive causes of neutropenia and the most common primary neutropenia disorders: constitutional/benign ethnic neutropenia, chronic idiopathic neutropenia, cyclic neutropenia, severe congenital neutropenia (due to mutations in ELANE, GFI1, HAX1, G6PC3, VPS45, JAGN1, CSF3R, SRP54, CLPB and WAS), GATA2 deficiency, Warts, hypogammaglobulinaemia, infections and myelokathexis syndrome, Shwachman-Diamond Syndrome, the lysosomal storage disorders with neutropenia: Chediak-Higashi, Hermansky-Pudlak, and Griscelli syndromes, Cohen, and Barth syndromes. We also detail characteristic cytogenetic and molecular factors at diagnosis and in progression to myelodysplastic syndrome/leukaemia.

Severe neutropenia unrelated to clozapine in patients receiving clozapine

BACKGROUND

Clozapine is known to cause agranulocytosis. Mandatory monitoring schemes are aimed at reducing the risk of agranulocytosis and of the consequences of agranulocytosis. All cases of agranulocytosis occurring in people prescribed clozapine are assumed to be caused by clozapine.

METHODS

In a previous study, we examined a cohort of patients listed on our hospital database as having had clozapine-induced agranulocytosis and applied specific criteria to identify those with confirmed clozapine-related, life-threatening agranulocytosis. In this study, we examine the cases not meeting these specific criteria.

RESULTS

In the original study, 9 of 23 cases met the criteria for clozapine-induced, life-threatening agranulocytosis. Of the 13 remaining cases for whom data were available, 5 were probably caused by clozapine but were not life-threatening. Three cases were the result of concomitant cancer chemotherapy. Three were anomalous results probably related to measurement error. For the remaining two cases, the cause was not identified.

CONCLUSION

Not all cases of agranulocytosis occurring in people taking clozapine are caused by clozapine. The widely used threshold criterion-based diagnosis overestimates the risk of agranulocytosis. True clozapine-related agranulocytosis is best identified by pattern-based criteria: rapid fall in neutrophil counts over around 2 weeks to below 0.5 × 109/L for two consecutive days (unless clozapine is stopped very early or granulocyte colony stimulating factor is given) where other possible causes (benign ethnic neutropenia, cancer chemotherapy) can be ruled out.

Clinical associations with a polygenic predisposition to benign lower white blood cell counts

Polygenic variation unrelated to disease contributes to interindividual variation in baseline white blood cell (WBC) counts, but its clinical significance is uncharacterized. We investigated the clinical consequences of a genetic predisposition toward lower WBC counts among 89,559 biobank participants from tertiary care centers using a polygenic score for WBC count (PGSWBC) comprising single nucleotide polymorphisms not associated with disease. A predisposition to lower WBC counts was associated with a decreased risk of identifying pathology on a bone marrow biopsy performed for a low WBC count (odds-ratio = 0.55 per standard deviation increase in PGSWBC [95%CI, 0.30-0.94], p = 0.04), an increased risk of leukopenia (a low WBC count) when treated with a chemotherapeutic (n = 1724, hazard ratio [HR] = 0.78 [0.69-0.88], p = 4.0 × 10-5) or immunosuppressant (n = 354, HR = 0.61 [0.38-0.99], p = 0.04). A predisposition to benign lower WBC counts was associated with an increased risk of discontinuing azathioprine treatment (n = 1,466, HR = 0.62 [0.44-0.87], p = 0.006). Collectively, these findings suggest that there are genetically predisposed individuals who are susceptible to escalations or alterations in clinical care that may be harmful or of little benefit.

A genome-wide association study of neutrophil count in individuals associated to an African continental ancestry group facilitates studies of malaria pathogenesis

BACKGROUND

'Benign ethnic neutropenia' (BEN) is a heritable condition characterized by lower neutrophil counts, predominantly observed in individuals of African ancestry, and the genetic basis of BEN remains a subject of extensive research. In this study, we aimed to dissect the genetic architecture underlying neutrophil count variation through a linear-mixed model genome-wide association study (GWAS) in a population of African ancestry (N = 5976). Malaria caused by P. falciparum imposes a tremendous public health burden on people living in sub-Saharan Africa. Individuals living in malaria endemic regions often have a reduced circulating neutrophil count due to BEN, raising the possibility that reduced neutrophil counts modulate severity of malaria in susceptible populations. As a follow-up, we tested this hypothesis by conducting a Mendelian randomization (MR) analysis of neutrophil counts on severe malaria (MalariaGEN, N = 17,056).

RESULTS

We carried out a GWAS of neutrophil count in individuals associated to an African continental ancestry group within UK Biobank, identifying 73 loci (r2 = 0.1) and 10 index SNPs (GCTA-COJO loci) associated with neutrophil count, including previously unknown rare loci regulating neutrophil count in a non-European population. BOLT-LMM was reliable when conducted in a non-European population, and additional covariates added to the model did not largely alter the results of the top loci or index SNPs. The two-sample bi-directional MR analysis between neutrophil count and severe malaria showed the greatest evidence for an effect between neutrophil count and severe anaemia, although the confidence intervals crossed the null.

CONCLUSION

Our GWAS of neutrophil count revealed unique loci present in individuals of African ancestry. We note that a small sample-size reduced our power to identify variants with low allele frequencies and/or low effect sizes in our GWAS. Our work highlights the need for conducting large-scale biobank studies in Africa and for further exploring the link between neutrophils and severe malaria.

Human 'knockouts' of CSF3 display severe congenital neutropenia

Colony-stimulating factor 3 (CSF3) is a key factor in neutrophil production and function, and recombinant forms have been used clinically for decades to treat congenital and acquired neutropenia. Although biallelic inactivation of its receptor CSF3R is a well-established cause of severe congenital neutropenia (SCN), no corresponding Mendelian disease has been ascribed to date to CSF3. Here, we describe three patients from two families each segregating a different biallelic inactivating variant in CSF3 with SCN. Complete deficiency of CSF3 as a result of nonsense-mediated decay (NMD) could be demonstrated on RT-PCR using skin fibroblasts-derived RNA. The phenotype observed in this cohort mirrors that documented in mouse and zebrafish models of CSF3 deficiency. Our results suggest that CSF3 deficiency in humans causes a novel autosomal recessive form of SCN.

Clinical consequences of a polygenic predisposition to benign lower white blood cell counts: Consequences of benign WBC count genetics

Polygenic variation unrelated to disease contributes to interindividual variation in baseline white blood cell (WBC) counts, but its clinical significance is undefined. We investigated the clinical consequences of a genetic predisposition toward lower WBC counts among 89,559 biobank participants from tertiary care centers using a polygenic score for WBC count (PGSWBC) comprising single nucleotide polymorphisms not associated with disease. A predisposition to lower WBC counts was associated with a decreased risk of identifying pathology on a bone marrow biopsy performed for a low WBC count (odds-ratio=0.55 per standard deviation increase in PGSWBC [95%CI, 0.30 - 0.94], p=0.04), an increased risk of leukopenia (a low WBC count) when treated with a chemotherapeutic (n=1,724, hazard ratio [HR]=0.78 [0.69 - 0.88], p=4.0×10-5) or immunosuppressant (n=354, HR=0.61 [0.38 - 0.99], p=0.04). A predisposition to benign lower WBC counts was associated with an increased risk of discontinuing azathioprine treatment (n=1,466, HR=0.62 [0.44 - 0.87], p=0.006). Collectively, these findings suggest that a WBC count polygenic score identifies individuals who are susceptible to escalations or alterations in clinical care that may be harmful or of little benefit.

Shared genetic architecture of blood eosinophil counts and asthma in UK Biobank

Rationale

Asthma is a complex, heterogeneous disease strongly associated with type 2 inflammation, and blood eosinophil counts guide therapeutic interventions in moderate and severe asthma. Eosinophils are leukocytes involved in type 2 immune responses. Despite these critical associations between asthma and blood eosinophil counts, the shared genetic architecture of these two traits remains unknown. The objective of the present study was to characterise the genetic architecture of blood eosinophil counts and asthma in the UK Biobank.

Methods

We performed genome-wide association studies (GWAS) of doctor-diagnosed asthma, blood eosinophil, neutrophil, lymphocyte and monocyte counts in the UK Biobank. Genetic correlation analysis was performed on GWAS results and validated in the Trans-National Asthma Genetic Consortium (TAGC) study of asthma.

Results

GWAS of doctor-diagnosed asthma and blood eosinophil counts in the UK Biobank identified 585 and 3429 significant variants, respectively. STAT6, a transcription factor involved in interleukin-4 signalling, was a key shared pathway between asthma and blood eosinophil counts. Genetic correlation analysis demonstrated a positive correlation between doctor-diagnosed asthma and blood eosinophil counts (r=0.38±0.10, correlation±se; p=4.7×10-11). As a validation of this association, we found a similar correlation between TAGC and blood eosinophil counts in the UK Biobank (0.37±0.08, correlation±se; p=1.2×10-6).

Conclusions

These findings define the shared genetic architecture between blood eosinophil counts and asthma risk in subjects of European ancestry and point to a genetic link to the STAT6 signalling pathway in these two traits.

A statistical framework to identify cell types whose genetically regulated proportions are associated with complex diseases

Finding disease-relevant tissues and cell types can facilitate the identification and investigation of functional genes and variants. In particular, cell type proportions can serve as potential disease predictive biomarkers. In this manuscript, we introduce a novel statistical framework, cell-type Wide Association Study (cWAS), that integrates genetic data with transcriptomics data to identify cell types whose genetically regulated proportions (GRPs) are disease/trait-associated. On simulated and real GWAS data, cWAS showed good statistical power with newly identified significant GRP associations in disease-associated tissues. More specifically, GRPs of endothelial and myofibroblasts in lung tissue were associated with Idiopathic Pulmonary Fibrosis and Chronic Obstructive Pulmonary Disease, respectively. For breast cancer, the GRP of blood CD8+ T cells was negatively associated with breast cancer (BC) risk as well as survival. Overall, cWAS is a powerful tool to reveal cell types associated with complex diseases mediated by GRPs.

Evaluating significance of European-associated index SNPs in the East Asian population for 31 complex phenotypes

BACKGROUND

Genome-wide association studies (GWASs) have identified many single-nucleotide polymorphisms (SNPs) associated with complex phenotypes in the European (EUR) population; however, the extent to which EUR-associated SNPs can be generalized to other populations such as East Asian (EAS) is not clear.

RESULTS

By leveraging summary statistics of 31 phenotypes in the EUR and EAS populations, we first evaluated the difference in heritability between the two populations and calculated the trans-ethnic genetic correlation. We observed the heritability estimates of some phenotypes varied substantially across populations and 53.3% of trans-ethnic genetic correlations were significantly smaller than one. Next, we examined whether EUR-associated SNPs of these phenotypes could be identified in EAS using the trans-ethnic false discovery rate method while accounting for winner's curse for SNP effect in EUR and difference of sample sizes in EAS. We found on average 54.5% of EUR-associated SNPs were also significant in EAS. Furthermore, we discovered non-significant SNPs had higher effect heterogeneity, and significant SNPs showed more consistent linkage disequilibrium and allele frequency patterns between the two populations. We also demonstrated non-significant SNPs were more likely to undergo natural selection.

CONCLUSIONS

Our study revealed the extent to which EUR-associated SNPs could be significant in the EAS population and offered deep insights into the similarity and diversity of genetic architectures underlying phenotypes in distinct ancestral groups.

A genome-wide association study identifies distinct variants associated with pulmonary function among European and African ancestries from the UK Biobank

Pulmonary function is an indicator of well-being, and pulmonary pathologies are the third major cause of death worldwide. We analysed the UK Biobank genome-wide association summary statistics of pulmonary function for Europeans and individuals of recent African descent to identify variants associated with the trait in the two ancestries. Here, we show 627 variants in Europeans and 3 in Africans associated with three pulmonary function parameters. In addition to the 110 variants in Europeans previously reported to be associated with phenotypes related to pulmonary function, we identify 279 novel loci, including an ISX intergenic variant rs369476290 on chromosome 22 in Africans. Remarkably, we find no shared variants among Africans and Europeans. Furthermore, enrichment analyses of variants separately for each ancestry background reveal significant enrichment for terms related to pulmonary phenotypes in Europeans but not Africans. Further analysis of studies of pulmonary phenotypes reveals that individuals of European background are disproportionally overrepresented in datasets compared to Africans, with the gap widening over the past five years. Our findings extend our understanding of the different variants that modify the pulmonary function in Africans and Europeans, a promising finding for future GWASs and medical studies.

Correlations between complex human phenotypes vary by genetic background, gender, and environment

We develop a closed-form Haseman-Elston estimator for genetic and environmental correlation coefficients between complex phenotypes, which we term HEc, that is as precise as GCTA yet ∼20× faster. We estimate genetic and environmental correlations between over 7,000 phenotype pairs in subgroups from the Trans-Omics in Precision Medicine (TOPMed) program. We demonstrate substantial differences in both heritabilities and genetic correlations for multiple phenotypes and phenotype pairs between individuals of self-reported Black, Hispanic/Latino, and White backgrounds. We similarly observe differences in many of the genetic and environmental correlations between genders. To estimate the contribution of genetics to the observed phenotypic correlation, we introduce "fractional genetic correlation" as the fraction of phenotypic correlation explained by genetics. Finally, we quantify the enrichment of correlations between phenotypic domains, each of which is comprised of multiple phenotypes. Altogether, we demonstrate that the observed correlations between complex human phenotypes depend on the genetic background of the individuals, their gender, and their environment.

The causal relationship between white blood cell counts and hepatocellular carcinoma: a Mendelian randomization study

BACKGROUND

Most of hepatocellular carcinoma (HCC) arises on the background of chronic inflammation. The presence of infiltrating inflammatory cells is associated with tumour initiation, progression and clinical response to treatment. The influence of white blood cell (WBC) subtype counts on HCC progression remains unclear.

METHODS

In this study, we performed a Mendelian randomization (MR) study with the validation of two datasets. The summary data for WBC counts were extracted from a recent large GWAS of individuals of European ancestry. The GWAS data related to HCC were obtained from the UK Biobank (UKB). Univariable and multivariable MR analyses were used to identify risk factors genetically associated with HCC risks.

RESULTS

In the discovery dataset, multivariable MR analysis revealed that sum basophil neutrophil counts had an independent causal effect on the occurrence of HCC, with the sum basophil neutrophil counts as follows: (OR = 0.437, P = 0.003, CI 0.252-0.757). Similarly, in the validation dataset, total basophil neutrophil counts were also been identified as an independent risk factor for HCC. The sum basophil neutrophil counts were as follows: (OR = 0.574, P = 0.021, CI 0.358-0.920).

CONCLUSION

In the European population, genetically predicted lower total basophil neutrophil counts might be an independent risk factor for HCC.

Sickle cell cerebrovascular reactivity to a CO2 stimulus: Too little, too slow

Background: Despite increased cerebral blood flow (CBF), cerebral infarcts occur in patients with sickle cell disease (SCD). This suggests increased CBF does not meet metabolic demand possibly due to compromised cerebral vasodilatory response. Hypothesis: In adult SCD patients, cerebrovascular reactivity (CVR) and speed of vasodilatory response (tau) to a standardized vasodilatory stimulus, are reduced compared to normal subjects.

Methods: Functional brain imaging performed as part of routine care in adult SCD patients without known large vessel cerebral vasculopathy was reviewed retrospectively. CVR was calculated as the change in CBF measured as the blood-oxygenation-level-dependent (BOLD)-magnetic resonance imaging signal, in response to a standard vasoactive stimulus of carbon dioxide (CO₂). The tau corresponding to the best fit between the convolved end-tidal partial pressures of CO₂ and BOLD signal was defined as the speed of vascular response. CVR and tau were normalized using a previously generated atlas of 42 healthy controls.

Results: Fifteen patients were included. CVR was reduced in grey and white matter (mean Z-score for CVR −0.5 [−1.8 to 0.3] and −0.6 [−2.3 to 0.7], respectively). Tau Z-scores were lengthened in grey and white matter (+0.9 [−0.5 to 3.3] and +0.8 [−0.7 to 2.7], respectively). Hematocrit was the only significant independent predictor of CVR on multivariable regression.

Conclusion: Both measures of cerebrovascular health (CVR and tau) in SCD patients were attenuated compared to normal controls. These findings show that CVR represents a promising tool to assess disease state, stroke risk, and therapeutic efficacy of treatments in SCD and merits further investigation.

A high-density genome-wide association with absolute blood monocyte count in domestic sheep identifies novel loci

Monocytes are a core component of the immune system that arise from bone marrow and differentiate into cells responsible for phagocytosis and antigen presentation. Their derivatives are often responsible for the initiation of the adaptive immune response. Monocytes and macrophages are central in both controlling and propagating infectious diseases such as infection by Coxiella burnetii and small ruminant lentivirus in sheep. Genotypes from 513 Rambouillet, Polypay, and Columbia sheep (Ovis aries) were generated using the Ovine SNP50 BeadChip. Of these sheep, 222 animals were subsequently genotyped with the Ovine Infinium^® HD SNP BeadChip to increase SNP coverage. Data from the 222 HD genotyped sheep were combined with the data from an additional 258 unique sheep to form a 480-sheep reference panel; this panel was used to impute the low-density genotypes to the HD genotyping density. Then, a genome-wide association analysis was conducted to identify loci associated with absolute monocyte counts from blood. The analysis used a single-locus mixed linear model implementing EMMAX with age and ten principal components as fixed effects. Two genome-wide significant peaks (p < 5x10^-7) were identified on chromosomes 9 and 1, and ten genome-wide suggestive peaks (p < 1x10^-5) were identified on chromosomes 1, 2, 3, 4, 9, 10, 15, and 16. The identified loci were within or near genes including KCNK9, involved into cytokine production, LY6D, a member of a superfamily of genes, some of which subset monocyte lineages, and HMGN1, which encodes a chromatin regulator associated with myeloid cell differentiation. Further investigation of these loci is being conducted to understand their contributions to monocyte counts. Investigating the genetic basis of monocyte lineages and numbers may in turn provide information about pathogens of veterinary importance and elucidate fundamental immunology.

GMEPS: a fast and efficient likelihood approach for genome-wide mediation analysis under extreme phenotype sequencing

Due to many advantages such as higher statistical power of detecting the association of genetic variants in human disorders and cost saving, extreme phenotype sequencing (EPS) is a rapidly emerging study design in epidemiological and clinical studies investigating how genetic variations associate with complex phenotypes. However, the investigation of the mediation effect of genetic variants on phenotypes is strictly restrictive under the EPS design because existing methods cannot well accommodate the non-random extreme tails sampling process incurred by the EPS design. In this paper, we propose a likelihood approach for testing the mediation effect of genetic variants through continuous and binary mediators on a continuous phenotype under the EPS design (GMEPS). Besides implementing in EPS design, it can also be utilized as a general mediation analysis procedure. Extensive simulations and two real data applications of a genome-wide association study of benign ethnic neutropenia under EPS design and a candidate-gene study of neurocognitive performance in patients with sickle cell disease under random sampling design demonstrate the superiority of GMEPS under the EPS design over widely used mediation analysis procedures, while demonstrating compatible capabilities under the general random sampling framework.

Genetic Influence on Frequencies of Myeloid-Derived Cell Subpopulations in Mouse

Differences in frequencies of blood cell subpopulations were reported to influence the course of infections, atopic and autoimmune diseases, and cancer. We have discovered a unique mouse strain B10.O20 containing extremely high frequency of myeloid-derived cells (MDC) in spleen. B10.O20 carries 3.6% of genes of the strain O20 on the C57BL/10 genetic background. It contains much higher frequency of CD11b⁺Gr1⁺ cells in spleen than both its parents. B10.O20 carries O20-derived segments on chromosomes 1, 15, 17, and 18. Their linkage with frequencies of blood cell subpopulations in spleen was tested in F₂ hybrids between B10.O20 and C57BL/10. We found 3 novel loci controlling MDC frequencies: Mydc1, 2, and 3 on chromosomes 1, 15, and 17, respectively, and a locus controlling relative spleen weight (Rsw1) that co-localizes with Mydc3 and also influences proportion of white and red pulp in spleen. Mydc1 controls numbers of CD11b⁺Gr1⁺ cells. Interaction of Mydc2 and Mydc3 regulates frequency of CD11b⁺Gr1⁺ cells and neutrophils (Gr1⁺Siglec-F^- cells from CD11b⁺ cells). Interestingly, Mydc3/Rsw1 is orthologous with human segment 6q21 that was shown previously to determine counts of white blood cells. Bioinformatics analysis of genomic sequence of the chromosomal segments bearing these loci revealed polymorphisms between O20 and C57BL/10 that change RNA stability and genes’ functions, and we examined expression of relevant genes. This identified potential candidate genes Smap1, Vps52, Tnxb, and Rab44. Definition of genetic control of MDC can help to personalize therapy of diseases influenced by these cells.

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.

Insights into the genetic architecture of haematological traits from deep phenotyping and whole-genome sequencing for two Mediterranean isolated populations

Haematological traits are linked to cardiovascular, metabolic, infectious and immune disorders, as well as cancer. Here, we examine the role of genetic variation in shaping haematological traits in two isolated Mediterranean populations. Using whole-genome sequencing data at 22× depth for 1457 individuals from Crete (MANOLIS) and 1617 from the Pomak villages in Greece, we carry out a genome-wide association scan for haematological traits using linear mixed models. We discover novel associations (p < 5 × 10^–9) of five rare non-coding variants with alleles conferring effects of 1.44–2.63 units of standard deviation on red and white blood cell count, platelet and red cell distribution width. Moreover, 10.0% of individuals in the Pomak population and 6.8% in MANOLIS carry a pathogenic mutation in the Haemoglobin Subunit Beta (HBB) gene. The mutational spectrum is highly diverse (10 different mutations). The most frequent mutation in MANOLIS is the common Mediterranean variant IVS-I-110 (G>A) (rs35004220). In the Pomak population, c.364C>A (“HbO-Arab”, rs33946267) is most frequent (4.4% allele frequency). We demonstrate effects on haematological and other traits, including bilirubin, cholesterol, and, in MANOLIS, height and gestation age. We find less severe effects on red blood cell traits for HbS, HbO, and IVS-I-6 (T>C) compared to other b+ mutations. Overall, we uncover allelic diversity of HBB in Greek isolated populations and find an important role for additional rare variants outside of HBB.

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.

Circulating white blood cells and lung function impairment: the observational studies and Mendelian randomization analysis

BACKGROUND

Circulating white blood cell (WBC) counts have been related to lung function impairment, but causal relationship was not established. We aimed to evaluate independent effects and causal relationships of WBC subtypes with lung function.

METHODS

The 19,159 participants from NHANES 2011-2012 (n = 3570), coke-oven workers (COW, n = 1762) and Dongfeng-Tongji (DFTJ, n = 13,827) cohorts were included in the observational studies. The associations between circulating counts of WBC subtypes and prebronchodilator lung function were evaluated by linear regression models and LASSO regression was used to select effective WBC subtypes. Summary statistics for WBC-associated SNPs were extracted from literature, and Mendelian randomization (MR) analysis with inverse-variance weighted (IVW) method was applied to estimate the causal effects of total WBC and subtypes on lung function among 4012 subjects from COW (n = 1126) and DFTJ cohorts (n = 2886).

RESULTS

Total WBC counts were negatively associated with lung function among three populations and their pooled analysis indicated that per 1 × 109 cells/L increase in total WBC was associated with 36.13 (95% CI: 30.35, 41.91) mL and 25.23 (95% CI: 19.97, 30.50) mL decrease in FVC and FEV1, respectively. Independent associations with lung function were found for neutrophils, monocytes, eosinophils and basophils (all p < .05), except lymphocytes. Besides, IVW MR analysis showed that genetically predicted total WBC and neutrophil counts were associated with reduced FVC (p = .017 and .021, respectively) and FEV1 (p = .048 and .043, respectively).

CONCLUSIONS

WBC subtypes were independently associated with lower lung function except lymphocytes. Our findings suggest that circulating neutrophils may be causal factors in lung function impairment.KEY MESSAGESWhite blood cell (WBC) subtypes were negatively associated with lung function level except lymphocytes in the observational studies.Associations of WBC subtypes with lung function may be modified by sex and smoking.Mendelian randomization analysis shows that neutrophils may be causal factors in lung function impairment.

Whole-genome sequencing in diverse subjects identifies genetic correlates of leukocyte traits: The NHLBI TOPMed program

Many common and rare variants associated with hematologic traits have been discovered through imputation on large-scale reference panels. However, the majority of genome-wide association studies (GWASs) have been conducted in Europeans, and determining causal variants has proved challenging. We performed a GWAS of total leukocyte, neutrophil, lymphocyte, monocyte, eosinophil, and basophil counts generated from 109,563,748 variants in the autosomes and the X chromosome in the Trans-Omics for Precision Medicine (TOPMed) program, which included data from 61,802 individuals of diverse ancestry. We discovered and replicated 7 leukocyte trait associations, including (1) the association between a chromosome X, pseudo-autosomal region (PAR), noncoding variant located between cytokine receptor genes (CSF2RA and CLRF2) and lower eosinophil count; and (2) associations between single variants found predominantly among African Americans at the S1PR3 (9q22.1) and HBB (11p15.4) loci and monocyte and lymphocyte counts, respectively. We further provide evidence indicating that the newly discovered eosinophil-lowering chromosome X PAR variant might be associated with reduced susceptibility to common allergic diseases such as atopic dermatitis and asthma. Additionally, we found a burden of very rare FLT3 (13q12.2) variants associated with monocyte counts. Together, these results emphasize the utility of whole-genome sequencing in diverse samples in identifying associations missed by European-ancestry-driven GWASs.

Association between Polygenetic Risk Scores of Low Immunity and Interactions between These Scores and Moderate Fat Intake in a Large Cohort

White blood cell (WBC) counts represent overall immunity. However, a few studies have been conducted to explore the genetic impacts of immunity and their interaction with lifestyles. We aimed to identify genetic variants associated with a low-WBC risk and document interactions between polygenetic risk scores (PRS), lifestyle factors, and nutrient intakes that influence low-WBC risk in a large hospital-based cohort. Single nucleotide polymorphisms (SNPs) were selected by genome-wide association study of participants with a low-WBC count (<4 × 109/L, n = 4176; low-WBC group) or with a normal WBC count (≥4 × 109/L, n = 36,551; control group). The best model for gene-gene interactions was selected by generalized multifactor dimensionality reduction. PRS was generated by summing selected SNP risk alleles of the best genetic model. Adjusted odds ratio (ORs) of the low-WBC group were 1.467 (1.219-1.765) for cancer incidence risk and 0.458 (0.385-0.545) for metabolic syndrome risk. Vitamin D intake, plant-based diet, and regular exercise were positively related to the low-WBC group, but smoking and alcohol intake showed an inverse association. The 7 SNPs included in the best genetic model were PSMD3_rs9898547, LCT_rs80157389, HLA-DRB1_rs532162239 and rs3097649, HLA-C rs2308575, CDKN1A_rs3176337 and THRA_rs7502539. PRS with 7 SNP model were positively associated with the low-WBC risk by 2.123-fold (1.741 to 2.589). PRS interacted with fat intake and regular exercise but not with other nutrient intakes or lifestyles. The proportion with the low WBC in the participants with high-PRS was lower among those with moderate-fat intake and regular exercise than those with low-fat intake and no exercise. In conclusion, adults with high-PRS had a higher risk of a low WBC count, and they needed to be advised to have moderate fat intake (20-25 energy percent) and regular exercise.

Racial and ethnic disparities in pediatric cancer incidence among children and young adults in the United States by single year of age

BACKGROUND

Incidence rates of pediatric cancers in the United States are typically reported in 5-year age groups, obscuring variation by single year of age. Additionally, racial and ethnic variation in incidence is typically presented in broad categories rather than by narrow age ranges.

METHODS

The Surveillance, Epidemiology, and End Results (SEER) 18 data (2000-2017) were examined to calculate frequencies and age-adjusted incidence rates among individuals aged birth to 39 years. Incidence rate ratios (IRRs) and 95% confidence intervals (95% CIs) were estimated as the measure of association for rate comparisons by race and Hispanic origin overall and by single year of age.

RESULTS

Several histologic types showed substantial variation in race/ethnicity-specific and overall rates by single year of age. Overall, Black children and young adults experienced substantially decreased incidence of acute lymphoid leukemia (IRR, 0.52; 95% CI, 0.49-0.55) compared to Whites, and this decreased incidence was strongest at ages 1 through 7 years and 16 through 20 years. Hispanic individuals experienced decreased overall incidence of Hodgkin lymphoma (IRR, 0.50; 95% CI, 0.48-0.52) and astrocytoma (IRR, 0.54; 95% CI, 0.52-0.56) and increased risk of acute lymphoblastic leukemia (IRR, 1.46; 95% CI, 1.42-1.51) compared to non-Hispanic Whites, and the increased risk was strongest at ages 10 through 23 years. Substantial decreased risk across many tumor types was also observed for Asian/Pacific Islanders and American Indian/Alaska Natives.

CONCLUSIONS

Examination of incidence rates for pediatric cancers by narrow age groups may provide insights regarding etiological differences in subgroups. Additionally, variation in age-specific incidence rates by race and ethnicity may enable hypothesis generation on drivers of disparities observed.

Multi-ethnic genome-wide association analyses of white blood cell and platelet traits in the Population Architecture using Genomics and Epidemiology (PAGE) study

BACKGROUND

Circulating white blood cell and platelet traits are clinically linked to various disease outcomes and differ across individuals and ancestry groups. Genetic factors play an important role in determining these traits and many loci have been identified. However, most of these findings were identified in populations of European ancestry (EA), with African Americans (AA), Hispanics/Latinos (HL), and other races/ethnicities being severely underrepresented.

RESULTS

We performed ancestry-combined and ancestry-specific genome-wide association studies (GWAS) for white blood cell and platelet traits in the ancestrally diverse Population Architecture using Genomics and Epidemiology (PAGE) Study, including 16,201 AA, 21,347 HL, and 27,236 EA participants. We identified six novel findings at suggestive significance (P < 5E-8), which need confirmation, and independent signals at six previously established regions at genome-wide significance (P < 2E-9). We confirmed multiple previously reported genome-wide significant variants in the single variant association analysis and multiple genes using PrediXcan. Evaluation of loci reported from a Euro-centric GWAS indicated attenuation of effect estimates in AA and HL compared to EA populations.

CONCLUSIONS

Our results highlighted the potential to identify ancestry-specific and ancestry-agnostic variants in participants with diverse backgrounds and advocate for continued efforts in improving inclusion of racially/ethnically diverse populations in genetic association studies for complex traits.

Whole-genome sequencing association analysis of quantitative red blood cell phenotypes: The NHLBI TOPMed program

Whole-genome sequencing (WGS), a powerful tool for detecting novel coding and non-coding disease-causing variants, has largely been applied to clinical diagnosis of inherited disorders. Here we leveraged WGS data in up to 62,653 ethnically diverse participants from the NHLBI Trans-Omics for Precision Medicine (TOPMed) program and assessed statistical association of variants with seven red blood cell (RBC) quantitative traits. We discovered 14 single variant-RBC trait associations at 12 genomic loci, which have not been reported previously. Several of the RBC trait-variant associations (RPN1, ELL2, MIDN, HBB, HBA1, PIEZO1, and G6PD) were replicated in independent GWAS datasets imputed to the TOPMed reference panel. Most of these discovered variants are rare/low frequency, and several are observed disproportionately among non-European Ancestry (African, Hispanic/Latino, or East Asian) populations. We identified a 3 bp indel p.Lys2169del (g.88717175_88717177TCT[4]) (common only in the Ashkenazi Jewish population) of PIEZO1, a gene responsible for the Mendelian red cell disorder hereditary xerocytosis (MIM: 194380), associated with higher mean corpuscular hemoglobin concentration (MCHC). In stepwise conditional analysis and in gene-based rare variant aggregated association analysis, we identified several of the variants in HBB, HBA1, TMPRSS6, and G6PD that represent the carrier state for known coding, promoter, or splice site loss-of-function variants that cause inherited RBC disorders. Finally, we applied base and nuclease editing to demonstrate that the sentinel variant rs112097551 (nearest gene RPN1) acts through a cis-regulatory element that exerts long-range control of the gene RUVBL1 which is essential for hematopoiesis. Together, these results demonstrate the utility of WGS in ethnically diverse population-based samples and gene editing for expanding knowledge of the genetic architecture of quantitative hematologic traits and suggest a continuum between complex trait and Mendelian red cell disorders.

When non-Whiteness becomes a condition

The term "benign ethnic neutropenia" describes the phenotype of having an absolute neutrophil count (ANC) <1500 cells/μL with no increased risk of infection. It is most commonly seen in those of African ancestry. In addition, ANC reference ranges from countries in Africa emphasize that ANC levels <1500 cells/μL are common and harmless. The lower ANC levels are driven by the Duffy null [Fy(a-b-)] phenotype, which is protective against malaria and seen in 80% to 100% of those of sub-Saharan African ancestry and <1% of those of European descent. Benign ethnic neutropenia is clinically insignificant, but the average ANC values differ from what are typically seen in those of European descent. Thus, the predominantly White American medical system has described this as a condition. This labeling implicitly indicates that common phenotypes in non-White populations are abnormal or wrong. We believe that it is important to examine and rectify practices in hematology that contribute to systemic racism.

An Effective Machine Learning Approach for Identifying Non-Severe and Severe Coronavirus Disease 2019 Patients in a Rural Chinese Population: The Wenzhou Retrospective Study

This paper has proposed an effective intelligent prediction model that can well discriminate and specify the severity of Coronavirus Disease 2019 (COVID-19) infection in clinical diagnosis and provide a criterion for clinicians to weigh scientific and rational medical decision-making. With indicators as the age and gender of the patients and 26 blood routine indexes, a severity prediction framework for COVID-19 is proposed based on machine learning techniques. The framework consists mainly of a random forest and a support vector machine (SVM) model optimized by a slime mould algorithm (SMA). When the random forest was used to identify the key factors, SMA was employed to train an optimal SVM model. Based on the COVID-19 data, comparative experiments were conducted between RF-SMA-SVM and several well-known machine learning algorithms performed. The results indicate that the proposed RF-SMA-SVM not only achieves better classification performance and higher stability on four metrics, but also screens out the main factors that distinguish severe COVID-19 patients from non-severe ones. Therefore, there is a conclusion that the RF-SMA-SVM model can provide an effective auxiliary diagnosis scheme for the clinical diagnosis of COVID-19 infection.

CXC Chemokine Receptor Type 5 Gene Polymorphisms in a Cohort of Egyptian Patients with Diffuse Large B-Cell Lymphoma

BACKGROUND

The chemokine receptor CXCR5 is selectively expressed on B cells; it is involved in lymphocyte homing and the development of normal lymphoid tissue. Its principle ligand is CXCL13 or B lymphocyte chemoattractant. Three polymorphisms in the CXCR5 gene, rs148351692 C/G, rs6421571 C/T, and rs78440425 G/A, have been identified.

OBJECTIVE

To assess the genetic polymorphisms of CXCR5 and evaluate their possible contribution to the susceptibility and response to therapy of diffuse large B-cell lymphoma (DLBCL).

PATIENTS AND METHODS

Fifty DLBCL (not otherwise specified) patients and 50 control subjects were included in this study. CXCR5 genotypes were determined by PCR-RFLP.

RESULTS

Our study revealed that the CXCR5 rs148351692 C/G and rs6421571 C/T gene polymorphisms are associated with an increased risk of developing DLBCL (OR 28.57 [95% CI 8.96-96.56] and 3.45 [1.67-11.83] respectively), while CXCR5 rs78440425 G/A showed no association with the risk of lymphoma. Moreover, the double and triple combined gene polymorphisms are associated with an increased risk of developing DLBCL of approximately 120-fold and 105-fold, respectively. CXCR5 gene polymorphisms had no significant impact on disease outcome or response to therapy.

CONCLUSIONS

CXCR5 gene polymorphisms could be considered a potential risk factor for the development of DLBCL.

The Duffy-null genotype and risk of infection

Many medical treatments, from oncology to psychiatry, can lower white blood cell counts and thus access to these treatments can be restricted to individuals with normal levels of white blood cells, principally in order to minimize risk of serious infection. This adversely affects individuals of African or Middle Eastern ancestries who have on average a reduced number of circulating white blood cells, because of the Duffy-null (CC) genotype at rs2814778 in the ACKR1 gene. Here, we investigate whether the Duffy-null genotype is associated with the risk of infection using the UK Biobank sample and the iPSYCH Danish case-cohort study, two population-based samples from different countries and age ranges. We found that a high proportion of those with the Duffy-null genotype (21%) had a neutrophil count below the threshold often used as a cut-off for access to relevant treatments, compared with 1% of those with the TC/TT genotype. In addition we found that despite its strong association with lower average neutrophil counts, the Duffy-null genotype was not associated with an increased risk of infection, viral or bacterial. These results have widespread implications for the clinical treatment of individuals of African ancestry and indicate that neutrophil thresholds to access treatments could be lowered in individuals with the Duffy-null genotype without an increased risk of infection.

Combining clinical and candidate gene data into a risk score for azathioprine-associated leukopenia in routine clinical practice

Leukopenia is a serious, frequent side effect associated with azathioprine use. Currently, we use thiopurine methyltransferase (TPMT) testing to predict leukopenia in patients taking azathioprine. We hypothesized that a risk score incorporating additional clinical and genetic variables would improve the prediction of azathioprine-associated leukopenia. In the discovery phase, we developed four risk score models: (1) age, sex, and TPMT metabolizer status; (2) model 1 plus additional clinical variables; (3) sixty candidate single nucleotide polymorphisms; and (4) model 2 plus model 3. The area under the receiver-operating-characteristic curve (AUC) of the risk scores was 0.59 (95% CI: 0.54-0.64), 0.75 (0.71-0.80), 0.66 (0.61-0.71), and 0.78 (0.74-0.82) for models 1, 2, 3, and 4, respectively. During the replication phase, models 2 and 4 (AUC = 0.64, 95% CI: 0.59-0.70 and AUC = 0.63, 95% CI: 0.58-0.69, respectively) were significant in an independent group. Compared with TPMT testing alone, additional genetic and clinical variables improve the prediction of azathioprine-associated leukopenia.

In The Blood: Connecting Variant to Function In Human Hematopoiesis

Genome-wide association studies (GWAS) have identified thousands of genetic variants associated with a range of human diseases and traits. However, understanding the mechanisms by which these genetic variants have an impact on associated diseases and traits, often referred to as the variant-to-function (V2F) problem, remains a significant hurdle. Solving the V2F challenge requires us to identify causative genetic variants, relevant cell types/states, target genes, and mechanisms by which variants can cause diseases or alter phenotypic traits. We discuss emerging functional approaches that are being applied to tackle the V2F problem for blood cell traits, illuminating how human genetic variation can impact on key mechanisms in hematopoiesis, as well as highlighting future prospects for this nascent field.

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.

Evaluation of mitochondrial DNA copy number estimation techniques

Mitochondrial DNA copy number (mtDNA-CN), a measure of the number of mitochondrial genomes per cell, is a minimally invasive proxy measure for mitochondrial function and has been associated with several aging-related diseases. Although quantitative real-time PCR (qPCR) is the current gold standard method for measuring mtDNA-CN, mtDNA-CN can also be measured from genotyping microarray probe intensities and DNA sequencing read counts. To conduct a comprehensive examination on the performance of these methods, we use known mtDNA-CN correlates (age, sex, white blood cell count, Duffy locus genotype, incident cardiovascular disease) to evaluate mtDNA-CN calculated from qPCR, two microarray platforms, as well as whole genome (WGS) and whole exome sequence (WES) data across 1,085 participants from the Atherosclerosis Risk in Communities (ARIC) study and 3,489 participants from the Multi-Ethnic Study of Atherosclerosis (MESA). We observe mtDNA-CN derived from WGS data is significantly more associated with known correlates compared to all other methods (p < 0.001). Additionally, mtDNA-CN measured from WGS is on average more significantly associated with traits by 5.6 orders of magnitude and has effect size estimates 5.8 times more extreme than the current gold standard of qPCR. We further investigated the role of DNA extraction method on mtDNA-CN estimate reproducibility and found mtDNA-CN estimated from cell lysate is significantly less variable than traditional phenol-chloroform-isoamyl alcohol (p = 5.44x10-4) and silica-based column selection (p = 2.82x10-7). In conclusion, we recommend the field moves towards more accurate methods for mtDNA-CN, as well as re-analyze trait associations as more WGS data becomes available from larger initiatives such as TOPMed.

A phenome-wide association study (PheWAS) in the Population Architecture using Genomics and Epidemiology (PAGE) study reveals potential pleiotropy in African Americans

We performed a hypothesis-generating phenome-wide association study (PheWAS) to identify and characterize cross-phenotype associations, where one SNP is associated with two or more phenotypes, between thousands of genetic variants assayed on the Metabochip and hundreds of phenotypes in 5,897 African Americans as part of the Population Architecture using Genomics and Epidemiology (PAGE) I study. The PAGE I study was a National Human Genome Research Institute-funded collaboration of four study sites accessing diverse epidemiologic studies genotyped on the Metabochip, a custom genotyping chip that has dense coverage of regions in the genome previously associated with cardio-metabolic traits and outcomes in mostly European-descent populations. Here we focus on identifying novel phenome-genome relationships, where SNPs are associated with more than one phenotype. To do this, we performed a PheWAS, testing each SNP on the Metabochip for an association with up to 273 phenotypes in the participating PAGE I study sites. We identified 133 putative pleiotropic variants, defined as SNPs associated at an empirically derived p-value threshold of p<0.01 in two or more PAGE study sites for two or more phenotype classes. We further annotated these PheWAS-identified variants using publicly available functional data and local genetic ancestry. Amongst our novel findings is SPARC rs4958487, associated with increased glucose levels and hypertension. SPARC has been implicated in the pathogenesis of diabetes and is also known to have a potential role in fibrosis, a common consequence of multiple conditions including hypertension. The SPARC example and others highlight the potential that PheWAS approaches have in improving our understanding of complex disease architecture by identifying novel relationships between genetic variants and an array of common human phenotypes.

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.

Neutrophil extracellular traps drive inflammatory pathogenesis in malaria

Neutrophils are essential innate immune cells that extrude chromatin in the form of neutrophil extracellular traps (NETs) when they die. This form of cell death has potent immunostimulatory activity. We show that heme-induced NETs are essential for malaria pathogenesis. Using patient samples and a mouse model, we define two mechanisms of NET-mediated inflammation of the vasculature: activation of emergency granulopoiesis via granulocyte colony-stimulating factor production and induction of the endothelial cytoadhesion receptor intercellular adhesion molecule-1. Soluble NET components facilitate parasite sequestration and mediate tissue destruction. We demonstrate that neutrophils have a key role in malaria immunopathology and propose inhibition of NETs as a treatment strategy in vascular infections.

Genetic association of molecular traits: A help to identify causative variants in complex diseases

In the past 15 years, major progresses have been made in the understanding of the genetic basis of regulation of gene expression. These new insights have revolutionized our approach to resolve the genetic variation underlying complex diseases. Gene transcript levels were the first expression phenotypes that were studied. They are heritable and therefore amenable to genome-wide association studies. The genetic variants that modulate them are called expression quantitative trait loci. Their study has been extended to other molecular quantitative trait loci (molQTLs) that regulate gene expression at the various levels, from chromatin state to cellular responses. Altogether, these studies have generated a wealth of basic information on the genome-wide patterns of gene expression and their inter-individual variation. Most importantly, molQTLs have become an invaluable asset in the genetic study of complex diseases. Although the identification of the disease-causing variants on the basis of their overlap with molQTLs requires caution, molQTLs can help to prioritize the relevant candidate gene(s) in the disease-associated regions and bring a functional interpretation of the associated variants, therefore, bridging the gap between genotypes and clinical phenotypes.

Inferring the Nature of Missing Heritability in Human Traits Using Data from the GWAS Catalog

Thousands of genes responsible for many diseases and other common traits in humans have been detected by Genome Wide Association Studies (GWAS) in the last decade. However, candidate causal variants found so far usually explain only a small fraction of the heritability estimated by family data. The most common explanation for this observation is that the missing heritability corresponds to variants, either rare or common, with very small effect, which pass undetected due to a lack of statistical power. We carried out a meta-analysis using data from the NHGRI-EBI GWAS Catalog in order to explore the observed distribution of locus effects for a set of 42 complex traits and to quantify their contribution to narrow-sense heritability. With the data at hand, we were able to predict the expected distribution of locus effects for 16 traits and diseases, their expected contribution to heritability, and the missing number of loci yet to be discovered to fully explain the familial heritability estimates. Our results indicate that, for 6 out of the 16 traits, the additive contribution of a great number of loci is unable to explain the familial (broad-sense) heritability, suggesting that the gap between GWAS and familial estimates of heritability may not ever be closed for these traits. In contrast, for the other 10 traits, the additive contribution of hundreds or thousands of loci yet to be found could potentially explain the familial heritability estimates, if this were the case. Computer simulations are used to illustrate the possible contribution from nonadditive genetic effects to the gap between GWAS and familial estimates of heritability.

Benign ethnic neutropenia

Benign ethnic neutropenia (BEN) is one of the most common causes of chronic neutropenia seen in individuals of African, Middle Eastern and West Indian descent, affecting many individuals worldwide. Despite its prevalence, many physicians are not familiar with this benign condition, resulting in unnecessary evaluation and testing for neutropenia in otherwise healthy individuals. Clinically, patients with BEN are at no increased risk of infection despite their neutropenia. Implications of this condition are highlighted in those patients receiving therapies that have a known side effect of neutropenia, most commonly chemotherapy agents. Studies have suggested that disparities in survival among those patients receiving chemotherapy between patients of European decent and African decent may be attributed to measured neutropenia in these populations, questioning whether BEN could be an influential factor. This review encompasses all aspects of benign ethnic neutropenia, providing information about this condition and helping to guide clinical decision-making as to when an aggressive work up and referral are indicated and when it is appropriate to monitor. Additionally, we review the role of genetic studies in identifying the genes related to BEN, summarize the theories that offer the most accepted mechanisms behind the condition, and address the importance of pursuing larger studies to assess the implication of BEN in oncology patients as well as patients taking neutropenia-causing medications.

Blood donation barriers and facilitators of Sub-Saharan African migrants and minorities in Western high-income countries: a systematic review of the literature

OBJECTIVES

The present study aimed to gain more insight into, and summarise, blood donation determinants among migrants or minorities of Sub-Saharan heritage by systematically reviewing the current literature.

BACKGROUND

Sub-Saharan Africans are under-represented in the blood donor population in Western high-income countries. This causes a lack of specific blood types for transfusions and prevention of alloimmunisation among Sub-Saharan African patients.

METHODS/MATERIALS

Medline, EMBASE, PsycINFO and BIOSIS were searched for relevant empirical studies that focused on barriers and facilitators of blood donation among Sub-Saharan Africans in Western countries until 22 June 2017. Of the 679 articles screened by title and abstract, 152 were subsequently screened by full text. Paired reviewers independently assessed the studies based on predefined eligibility and quality criteria.

RESULTS

Of the 31 included studies, 24 used quantitative and 7 used qualitative research methods. Target cohorts varied from Black African Americans and refugees from Sub-Sahara Africa to specific Sub-Saharan migrant groups such as Comorians or Ethiopians. Main recurring barriers for Sub-Saharan Africans were haemoglobin deferral, fear of needles and pain, social exclusion, lack of awareness, negative attitudes and accessibility problems. Important recurring facilitators for Sub-Saharan Africans were altruism, free health checks and specific recruitment and awareness-raising campaigns.

CONCLUSION

The findings of this review can be used as a starting point to develop recruitment and retention strategies for Sub-Saharan African persons. Further research is needed to gain more insight in the role of these determinants in specific contexts as socioeconomic features, personal histories and host country regulations may differ per country.

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.