Genome-wide association of white blood cell counts in Hispanic/Latino Americans: the Hispanic Community Health Study/Study of Latinos

Differential white blood cell count and epigenetic clocks: a bidirectional Mendelian randomization study

BACKGROUND

Human aging and white blood cell (WBC) count are complex traits influenced by multiple genetic factors. Predictors of chronological age have been developed using epigenetic clocks. However, the bidirectional causal effects between epigenetic clocks and WBC count have not been fully examined.

METHODS

This study employed Mendelian randomization (MR) to analyze summary statistics from four epigenetic clocks involving 34,710 participants, alongside data from the Blood Cell Consortium encompassing 563,946 individuals. We primarily explored bidirectional causal relationships using the random-effects inverse-variance weighted method, supplemented by additional MR methods for comprehensive analysis. Additionally, multivariate MR was applied to investigate independent effects of WBC count on epigenetic age acceleration.

RESULTS

In the two-sample univariate MR (UVMR) analysis, we observed that a decrease in lymphocyte count markedly accelerated aging according to the PhenoAge, GrimAge, and HannumAge metrics (all P < 0.01, β < 0), though it did not affect Intrinsic Epigenetic Age Acceleration (IEAA). Conversely, an increase in neutrophil count significantly elevated PhenoAge levels (β: 0.38; 95% CI 0.14, 0.61; P = 1.65E-03 < 0.01). Reverse MR revealed no significant causal impacts of epigenetic clocks on overall WBC counts. Furthermore, in multivariate MR, the impact of lymphocyte counts on epigenetic aging metrics remained statistically significant. We also identified a marked causal association between neutrophil counts and PhenoAge, GrimAge, and HannumAge, with respective results showing strong associations (PhenoAge β: 0.78; 95% CI 0.47, 1.09; P = 8.26E-07; GrimAge β: 0.55; 95% CI 0.31, 0.79; P = 5.50E-06; HannumAge β: 0.42; 95% CI 0.18, 0.67; P = 6.30E-04). Likewise, eosinophil cell count demonstrated significant association with HannumAge (β: 0.33; 95% CI 0.13, 0.53; P = 1.43E-03 < 0.01).

CONCLUSION

These findings demonstrated that within WBCs, lymphocyte and neutrophil counts exert irreversible and independent causal effects on the acceleration of PhenoAge, GrimAge, and HannumAge. Our findings highlight the critical role of WBCs in influencing epigenetic clocks and underscore the importance of considering immune parameters when interpreting epigenetic age.

Bidirectional two-sample Mendelian randomization study of differential white blood cell counts and schizophrenia

BACKGROUND

Schizophrenia and white blood cell counts (WBC) are both complex and polygenic traits. Previous evidence suggests that increased WBC are associated with higher all-cause mortality, and other studies have found elevated WBC in first-episode psychosis and chronic schizophrenia. However, these observational findings may be confounded by antipsychotic exposures and their effects on WBC. Mendelian randomization (MR) is a useful method for examining the directions of genetically-predicted relationships between schizophrenia and WBC.

METHODS

We performed a two-sample MR using summary statistics from genome-wide association studies (GWAS) conducted by the Psychiatric Genomics Consortium Schizophrenia Workgroup (N = 130,644) and the Blood Cell Consortium (N = 563,946). The MR methods included inverse variance weighted (IVW), MR Egger, weighted median, MR-PRESSO, contamination mixture, and a novel approach called mixture model reciprocal causal inference (MRCI). False discovery rate was employed to correct for multiple testing.

RESULTS

Multiple MR methods supported bidirectional genetically-predicted relationships between lymphocyte count and schizophrenia: IVW (b = 0.026; FDR p-value = 0.008), MR Egger (b = 0.026; FDR p-value = 0.008), weighted median (b = 0.013; FDR p-value = 0.049), and MR-PRESSO (b = 0.014; FDR p-value = 0.010) in the forward direction, and IVW (OR = 1.100; FDR p-value = 0.021), MR Egger (OR = 1.231; FDR p-value < 0.001), weighted median (OR = 1.136; FDR p-value = 0.006) and MRCI (OR = 1.260; FDR p-value = 0.026) in the reverse direction. MR Egger (OR = 1.171; FDR p-value < 0.001) and MRCI (OR = 1.154; FDR p-value = 0.026) both suggested genetically-predicted eosinophil count is associated with schizophrenia, but MR Egger (b = 0.060; FDR p-value = 0.010) and contamination mixture (b = -0.013; FDR p-value = 0.045) gave ambiguous results on whether genetically predicted liability to schizophrenia would be associated with eosinophil count. MR Egger (b = 0.044; FDR p-value = 0.010) and MR-PRESSO (b = 0.009; FDR p-value = 0.045) supported genetically predicted liability to schizophrenia is associated with elevated monocyte count, and the opposite direction was also indicated by MR Egger (OR = 1.231; FDR p-value = 0.045). Lastly, unidirectional genetic liability from schizophrenia to neutrophil count were proposed by MR-PRESSO (b = 0.011; FDR p-value = 0.028) and contamination mixture (b = 0.011; FDR p-value = 0.045) method.

CONCLUSION

This MR study utilised multiple MR methods to obtain results suggesting bidirectional genetic genetically-predicted relationships for elevated lymphocyte counts and schizophrenia risk. In addition, moderate evidence also showed bidirectional genetically-predicted relationships between schizophrenia and monocyte counts, and unidirectional effect from genetic liability for eosinophil count to schizophrenia and from genetic liability for schizophrenia to neutrophil count. The influence of schizophrenia to eosinophil count is less certain. Our findings support the role of WBC in schizophrenia and concur with the hypothesis of neuroinflammation in schizophrenia.

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.

Robust identification of regulatory variants (eQTLs) using a differential expression framework developed for RNA-sequencing

BACKGROUND

A gap currently exists between genetic variants and the underlying cell and tissue biology of a trait, and expression quantitative trait loci (eQTL) studies provide important information to help close that gap. However, two concerns that arise with eQTL analyses using RNA-sequencing data are normalization of data across samples and the data not following a normal distribution. Multiple pipelines have been suggested to address this. For instance, the most recent analysis of the human and farm Genotype-Tissue Expression (GTEx) project proposes using trimmed means of M-values (TMM) to normalize the data followed by an inverse normal transformation.

RESULTS

In this study, we reasoned that eQTL analysis could be carried out using the same framework used for differential gene expression (DGE), which uses a negative binomial model, a statistical test feasible for count data. Using the GTEx framework, we identified 35 significant eQTLs (P < 5 × 10^-8) following the ANOVA model and 39 significant eQTLs (P < 5 × 10^-8) following the additive model. Using a differential gene expression framework, we identified 930 and six significant eQTLs (P < 5 × 10^-8) following an analytical framework equivalent to the ANOVA and additive model, respectively. When we compared the two approaches, there was no overlap of significant eQTLs between the two frameworks. Because we defined specific contrasts, we identified trans eQTLs that more closely resembled what we expect from genetic variants showing complete dominance between alleles. Yet, these were not identified by the GTEx framework.

CONCLUSIONS

Our results show that transforming RNA-sequencing data to fit a normal distribution prior to eQTL analysis is not required when the DGE framework is employed. Our proposed approach detected biologically relevant variants that otherwise would not have been identified due to data transformation to fit a normal distribution.

Quantifying portable genetic effects and improving cross-ancestry genetic prediction with GWAS summary statistics

Polygenic risk scores (PRS) calculated from genome-wide association studies (GWAS) of Europeans are known to have substantially reduced predictive accuracy in non-European populations, limiting their clinical utility and raising concerns about health disparities across ancestral populations. Here, we introduce a statistical framework named X-Wing to improve predictive performance in ancestrally diverse populations. X-Wing quantifies local genetic correlations for complex traits between populations, employs an annotation-dependent estimation procedure to amplify correlated genetic effects between populations, and combines multiple population-specific PRS into a unified score with GWAS summary statistics alone as input. Through extensive benchmarking, we demonstrate that X-Wing pinpoints portable genetic effects and substantially improves PRS performance in non-European populations, showing 14.1%-119.1% relative gain in predictive R2 compared to state-of-the-art methods based on GWAS summary statistics. Overall, X-Wing addresses critical limitations in existing approaches and may have broad applications in cross-population polygenic risk prediction.

AI-assisted discovery of an ethnicity-influenced driver of cell transformation in esophageal and gastroesophageal junction adenocarcinomas

Although Barrett's metaplasia of the esophagus (BE) is the only known precursor lesion to esophageal adenocarcinomas (EACs), drivers of cellular transformation in BE remain incompletely understood. We use an artificial intelligence-guided network approach to study EAC initiation and progression. Key predictions are subsequently validated in a human organoid model, in patient-derived biopsy specimens of BE, a case-control study of genomics of BE progression, and in a cross-sectional study of 113 patients with BE and EACs. Our model classified healthy esophagus from BE and BE from EACs in several publicly available gene expression data sets (n = 932 samples). The model confirmed that all EACs must originate from BE and pinpointed a CXCL8/IL8↔neutrophil immune microenvironment as a driver of cellular transformation in EACs and gastroesophageal junction adenocarcinomas. This driver is prominent in White individuals but is notably absent in African Americans (AAs). Network-derived gene signatures, independent signatures of neutrophil processes, CXCL8/IL8 expression, and an absolute neutrophil count (ANC) are associated with risk of progression. SNPs associated with changes in ANC by ethnicity (e.g., benign ethnic neutropenia [BEN]) modify that risk. Findings define a racially influenced immunological basis for cell transformation and suggest that BEN in AAs may be a deterrent to BE→EAC progression.

Metabolomic Associations of Asthma in the Hispanic Community Health Study/Study of Latinos

Asthma disproportionally affects Hispanic and/or Latino backgrounds; however, the relation between circulating metabolites and asthma remains unclear. We conducted a cross-sectional study associating 640 individual serum metabolites, as well as twelve metabolite modules, with asthma in 3347 Hispanic/Latino background participants (514 asthmatics, 15.36%) from the Hispanic/Latino Community Health Study/Study of Latinos. Using survey logistic regression, per standard deviation (SD) increase in 1-arachidonoyl-GPA (20:4) was significantly associated with 32% high odds of asthma after accounting for clinical risk factors (p = 6.27 × 10−5), and per SD of the green module, constructed using weighted gene co-expression network, was suggestively associated with 25% high odds of asthma (p = 0.006). In the stratified analyses by sex and Hispanic and/or Latino backgrounds, the effect of 1-arachidonoyl-GPA (20:4) and the green module was predominantly observed in women (OR = 1.24 and 1.37, p < 0.001) and people of Cuban and Puerto-Rican backgrounds (OR = 1.25 and 1.27, p < 0.01). Mutations in Fatty Acid Desaturase 2 (FADS2) affected the levels of 1-arachidonoyl-GPA (20:4), and Mendelian Randomization analyses revealed that high genetically regulated 1-arachidonoyl-GPA (20:4) levels were associated with increased odds of asthma (p < 0.001). The findings reinforce a molecular basis for asthma etiology, and the potential causal effect of 1-arachidonoyl-GPA (20:4) on asthma provides an opportunity for future intervention.

Immunological and psychosocial functioning in parents of children with cancer

PURPOSE

Research has shown that parents of children with cancer exhibit an altered immune profile compared to parents of healthy children, reflective of increased susceptibility to illness. These parents are also at risk for poorer psychosocial outcomes and quality of life. The current study compares peripheral blood cell analyses and psychosocial self-reports from parents of children being treated for cancer (n = 21) to parents of healthy children (n = 30).

METHODS

A blood sample was drawn from parents to analyze immune profiles. Parents also completed the Perceived Stress Scale (PSS), Medical Outcomes Study Short Form-36 (MOS), and Patient-Reported Outcomes Measurement Information System Short Form v1.0 Emotional Distress-Anxiety 8a, and Emotional Distress-Depression 8a (PROMIS). Mann-Whitney U tests and independent samples t-tests were conducted to examine differences in outcomes between parent groups.

RESULTS

Parents of children with cancer exhibited higher monocyte percentages in their peripheral blood compared to peers with healthy children. Parents of children with cancer also reported poorer psychosocial outcomes: higher perceived stress, higher anxiety and depression symptoms, more role disability resulting from emotional problems, poorer general and mental health, and poorer social functioning.

CONCLUSION

These findings support research that has shown a direct effect of chronic stress on the immune system. Symptoms reported by parents of children with cancer indicate unmet psychosocial needs that could potentially affect long-term health. Given the central role of parents in their children's cancer care, it is compelling to address and work to improve parent immunological and psychosocial well-being.

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.

Genetic architecture of gene regulation in Indonesian populations identifies QTLs associated with global and local ancestries

Lack of diversity in human genomics limits our understanding of the genetic underpinnings of complex traits, hinders precision medicine, and contributes to health disparities. To map genetic effects on gene regulation in the underrepresented Indonesian population, we have integrated genotype, gene expression, and CpG methylation data from 115 participants across three island populations that capture the major sources of genomic diversity in the region. In a comparison with European datasets, we identify eQTLs shared between Indonesia and Europe as well as population-specific eQTLs that exhibit differences in allele frequencies and/or overall expression levels between populations. By combining local ancestry and archaic introgression inference with eQTLs and methylQTLs, we identify regulatory loci driven by modern Papuan ancestry as well as introgressed Denisovan and Neanderthal variation. GWAS colocalization connects QTLs detected here to hematological traits, and further comparison with European datasets reflects the poor overall transferability of GWAS statistics across diverse populations. Our findings illustrate how population-specific genetic architecture, local ancestry, and archaic introgression drive variation in gene regulation across genetically distinct and in admixed populations and highlight the need for performing association studies on non-European populations.

Racial/ethnic differences in eligibility for asthma biologics among pediatric populations

BACKGROUND

Asthma is a heterogeneous disease. Clinical blood parameters differ by race/ethnicity and are used to distinguish asthma subtypes and inform therapies. Differences in subtypes may explain population-specific trends in asthma outcomes. However, these differences in racial/ethnic minority pediatric populations are unclear.

OBJECTIVE

We investigated the association of blood parameters and asthma subtypes with asthma outcomes and examined population-specific eligibility for biologic therapies in minority pediatric populations.

METHODS

Using data from 2 asthma case-control studies of pediatric minority populations, we performed case-control (N = 3738) and case-only (N = 2743) logistic regressions to quantify the association of blood parameters and asthma subtypes with asthma outcomes. Heterogeneity of these associations was tested using an interaction term between race/ethnicity and each exposure. Differences in therapeutic eligibility were investigated using chi-square tests.

RESULTS

Race/ethnicity modified the association between total IgE and asthma exacerbations. Elevated IgE level was associated with worse asthma outcomes in Puerto Ricans. Allergic asthma was associated with worse outcomes in Mexican Americans, whereas eosinophilic asthma was associated with worse outcomes in Puerto Ricans. A lower proportion of Puerto Ricans met dosing criteria for allergic asthma-directed biologic therapy than other groups. A higher proportion of Puerto Ricans qualified for eosinophilic asthma-directed biologic therapy than African Americans.

CONCLUSIONS

We found population-specific associations between blood parameters and asthma subtypes with asthma outcomes. Our findings suggest that eligibility for asthma biologic therapies differs across pediatric racial/ethnic populations. These findings call for more studies in diverse populations for equitable treatment of minority patients with asthma.

Association between systemic immune-inflammation index and neoadjuvant chemotherapy efficacy as well as prognosis in triple-negative breast cancer

OBJECTIVES

Triple-negative breast cancer (TNBC) is a highly aggressive subtype of breast cancer with a relatively poor prognosis. Neoadjuvant chemotherapy (NAC) is the main treatment method. Due to the heterogeneity of the tumor, the chemotherapy response of TNBC patients is significantly different. Inflammation is closely related to the occurrence and development of cancer. The systemic immune-inflammation index (SII) is an indicator that can comprehensively reflect the state of systemic inflammation. This study aims to explore the association between SII and the NAC efficacy as well as the prognosis in TNBC.

METHODS

The data of TNBC patients who underwent NAC and systemic treatment in Xiangya Hospital of Central South University from January 2015 to June 2019 were collected. According to the inclusion and exclusion criteria, 231 TNBC patients were finally included. The pre-NAC SII was calculated according to the blood routine results of the patients at 1 week before chemotherapy, and the patients were divided into a pre-NAC low SII group (SII<412, 115 cases) and a pre-NAC high SII group (SII≥412, 116 cases). The SII after chemotherapy was calculated according to the blood routine results of the patients at 2 to 3 months after the end of chemotherapy, and the patients were divided into a low SII group after chemotherapy (SII<474, 115 cases) and a high SII group after chemotherapy (SII≥474, 116 cases). Pearson's chi-square test was used to analyze the relationship between SII and other clinical characteristics of TNBC patients, and the relationship between the NAC efficacy and clinical characteristics of TNBC patients. Binary logistic regression analysis was used to find independent factors that affect the efficacy of NAC in TNBC patients. Kaplan-Meier curve analysis was used to analyze factors affecting the prognosis of TNBC patients. Cox regression model was used to find independent factors affecting the prognosis of TNBC patients.

RESULTS

Before NAC, the differences in SII between groups with different ages and tumor sizes were significant (P=0.007 and P=0.002, respectively); after chemotherapy, there were no significant differences in SII between different ages, tumor sizes, histological grades, lymph node staging, and Ki-67 groups (all P>0.05). There were 115 patients with low SII before NAC, with a pathological complete response (pCR) rate of 15.7%; there were 116 patients with high SII before NAC, with a pCR rate of 6.0%. Patients with low SII before NAC had a higher pCR rate than patients with high SII before NAC, and the difference was statistically significant (P=0.019).There were 156 patients with lymph node staging pN0, with a pCR rate of 14.7%; and there were 75 patients with lymph node staging pN1-pN2, with a pCR rate of 2.7%. Patients with lymph node staging pN0 had a higher pCR rate than those with lymph node staging pN1-pN2, and the difference was significant (P=0.006). During the follow-up, 34 patients had local recurrence or distant metastasis. The Kaplan-Meier survival curve showed that the 3-year disease-free survival (DFS) rates for patients with low SII before NAC and high SII before NAC were 87.8% and 82.8%, respectively, and the former was significantly higher than the latter (P=0.005); the 3-year DFS rates for patients with tumor sizes of T1-T2 and T3 were 89.0% and 67.5%, respectively, and the former was significantly higher than the latter (P=0.001); the 3-year DFS rates for patients with lymph node staging of pN0 and pN1-pN2 were 87.8% and 82.8%, respectively, and the former was significantly higher than the latter (P=0.009). Cox analysis showed that SII before NAC and tumor size were independent influencing factors of patients' DFS (P=0.038, P=0.010, respectively).

CONCLUSIONS

SII has important clinical significance in predicting the efficacy and prognosis of NAC in TNBC patients, and it has the potential to be a biomarker.

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.

Genetic risk scores for cardiometabolic traits in sub-Saharan African populations

BACKGROUND

There is growing support for the use of genetic risk scores (GRS) in routine clinical settings. Due to the limited diversity of current genomic discovery samples, there are concerns that the predictive power of GRS will be limited in non-European ancestry populations. GRS for cardiometabolic traits were evaluated in sub-Saharan Africans in comparison with African Americans and European Americans.

METHODS

We evaluated the predictive utility of GRS for 12 cardiometabolic traits in sub-Saharan Africans (AF; n = 5200), African Americans (AA; n = 9139) and European Americans (EUR; n = 9594). GRS were constructed as weighted sums of the number of risk alleles. Predictive utility was assessed using the additional phenotypic variance explained and the increase in discriminatory ability over traditional risk factors [age, sex and body mass index (BMI)], with adjustment for ancestry-derived principal components.

RESULTS

Across all traits, GRS showed up to a 5-fold and 20-fold greater predictive utility in EUR relative to AA and AF, respectively. Predictive utility was most consistent for lipid traits, with percentage increase in explained variation attributable to GRS ranging from 10.6% to 127.1% among EUR, 26.6% to 65.8% among AA and 2.4% to 37.5% among AF. These differences were recapitulated in the discriminatory power, whereby the predictive utility of GRS was 4-fold greater in EUR relative to AA and up to 44-fold greater in EUR relative to AF. Obesity and blood pressure traits showed a similar pattern of greater predictive utility among EUR.

CONCLUSIONS

This work demonstrates the poorer performance of GRS in AF and highlights the need to improve representation of multiple ethnic populations in genomic studies to ensure equitable clinical translation of GRS.

Synaptic processes and immune-related pathways implicated in Tourette syndrome

Tourette syndrome (TS) is a neuropsychiatric disorder of complex genetic architecture involving multiple interacting genes. Here, we sought to elucidate the pathways that underlie the neurobiology of the disorder through genome-wide analysis. We analyzed genome-wide genotypic data of 3581 individuals with TS and 7682 ancestry-matched controls and investigated associations of TS with sets of genes that are expressed in particular cell types and operate in specific neuronal and glial functions. We employed a self-contained, set-based association method (SBA) as well as a competitive gene set method (MAGMA) using individual-level genotype data to perform a comprehensive investigation of the biological background of TS. Our SBA analysis identified three significant gene sets after Bonferroni correction, implicating ligand-gated ion channel signaling, lymphocytic, and cell adhesion and transsynaptic signaling processes. MAGMA analysis further supported the involvement of the cell adhesion and trans-synaptic signaling gene set. The lymphocytic gene set was driven by variants in FLT3, raising an intriguing hypothesis for the involvement of a neuroinflammatory element in TS pathogenesis. The indications of involvement of ligand-gated ion channel signaling reinforce the role of GABA in TS, while the association of cell adhesion and trans-synaptic signaling gene set provides additional support for the role of adhesion molecules in neuropsychiatric disorders. This study reinforces previous findings but also provides new insights into the neurobiology of TS.

Functional annotation of genetic associations by transcriptome-wide association analysis provides insights into neutrophil development regulation

Genome-wide association studies (GWAS) have identified multiple genomic loci linked to blood cell traits, however understanding the biological relevance of these genetic loci has proven to be challenging. Here, we performed a transcriptome-wide association study (TWAS) integrating gene expression and splice junction usage in neutrophils (N = 196) with a neutrophil count GWAS (N = 173,480 individuals). We identified a total of 174 TWAS-significant genes enriched in target genes of master transcription factors governing neutrophil specification. Knockout of a TWAS candidate at chromosome 5q13.2, TAF9, in CD34⁺ hematopoietic and progenitor cells (HSPCs) using CRISPR/Cas9 technology showed a significant effect on neutrophil production in vitro. In addition, we identified 89 unique genes significant only for splice junction usage, thus emphasizing the importance of alternative splicing beyond gene expression underlying granulopoiesis. Our results highlight the advantages of TWAS, followed by gene editing, to determine the functions of GWAS loci implicated in hematopoiesis.

The Absolute Basophil Count

The absolute basophil count (cells/L) can be determined by manual counting of peripheral blood smears or using cell counting chambers as well as by automated hematology analyzers and fluorescence flow cytometry. Manual basophil counting of peripheral blood smears is currently regarded as the reference method, although the limitations of this method (distribution, observer, and statistical errors) are widely recognized. Automated hematology analyzers offer an advantage of larger numbers of counted cells and high throughput but are characterized by inconsistent analytical performance for basophil enumeration. Flow cytometric enumeration of circulating basophils using panels of monoclonal antibodies is being developed as novel candidate reference method for the absolute basophil count in peripheral blood. Basophil counting using fluorescence flow cytometry is characterized by high precision and statistical superiority. Emerging innovative technologies for absolute cell counts include imaging flow cytometry, mass cytometry, and on-chip blood counting, but their analytical performance for absolute basophil counts is yet to be established. Here, we describe various techniques for absolute basophil counting in peripheral blood including manual basophil counts in smears and hemocytometers and flow cytometric methodologies using double-platform, bead-based, and volumetric approaches.

FLT3 stop mutation increases FLT3 ligand level and risk of autoimmune thyroid disease

Autoimmune thyroid disease is the most common autoimmune disease and is highly heritable¹. Here, by using a genome-wide association study of 30,234 cases and 725,172 controls from Iceland and the UK Biobank, we find 99 sequence variants at 93 loci, of which 84 variants are previously unreported^2-7. A low-frequency (1.36%) intronic variant in FLT3 (rs76428106-C) has the largest effect on risk of autoimmune thyroid disease (odds ratio (OR) = 1.46, P = 2.37 × 10^-24). rs76428106-C is also associated with systemic lupus erythematosus (OR = 1.90, P = 6.46 × 10^-4), rheumatoid factor and/or anti-CCP-positive rheumatoid arthritis (OR = 1.41, P = 4.31 × 10^-4) and coeliac disease (OR = 1.62, P = 1.20 × 10^-4). FLT3 encodes fms-related tyrosine kinase 3, a receptor that regulates haematopoietic progenitor and dendritic cells. RNA sequencing revealed that rs76428106-C generates a cryptic splice site, which introduces a stop codon in 30% of transcripts that are predicted to encode a truncated protein, which lacks its tyrosine kinase domains. Each copy of rs76428106-C doubles the plasma levels of the FTL3 ligand. Activating somatic mutations in FLT3 are associated with acute myeloid leukaemia⁸ with a poor prognosis and rs76428106-C also predisposes individuals to acute myeloid leukaemia (OR = 1.90, P = 5.40 × 10^-3). Thus, a predicted loss-of-function germline mutation in FLT3 causes a reduction in full-length FLT3, with a compensatory increase in the levels of its ligand and an increased disease risk, similar to that of a gain-of-function mutation.

Establishment of Reference Intervals for Blood Cell Analysis of Adult Tibetan Farmers and Herdsmen Over 4100 Meters Above Sea Level in Tibet Based on a Health Survey

Yang, Wei, Song Zhao, Dan Liu, Guangming Su, Dongwei Zhang, Gui La, Yu Deng, Pian Ni, and Xiuru Guan. Establishment of reference intervals for blood cell analysis of adult Tibetan farmers and herdsmen over 4100 meters above sea level in Tibet based on a health survey. High Alt Med Biol. 21:223-231, 2020. Background: High altitude has substantial influence on reference intervals (RIs) for blood cell analysis. However, the currently used RIs for blood cell analysis in China have not been investigated in people from high-altitude areas, which may potentially cause confusion regarding disease diagnosis and treatment of Tibetan farmers and herdsmen from high-altitude areas. Therefore, it is necessary to establish blood cell analysis RIs that are suitable for people in high-altitude areas. Methods: In total, 1145 healthy Tibetan farmers and herdsmen were recruited who have lived at an altitude from 4100 to 5280 m in Kangma, Tibet, for a long time. The C28-A3 guideline from the Clinical and Laboratory Standards Institute was used to establish the RIs for hemoglobin (HGB) concentrations and red blood cell (RBC), white blood cell (WBC), and platelet (PLT) counts from the blood results of Tibetan participants. Results: There were significant gender differences in blood cell RIs for Tibetan farmers and herdsmen who lived above 4100 m above sea level in Kangma, Tibet. Compared with currently used RIs in China, RIs for HGB concentrations and RBCs in Kangma were significantly higher, whereas WBCs were similar and PLT count was lower. Conclusions: The currently used RIs in China may not be applicable to adult Tibetan farmers and herdsmen above 4100 m in Tibet. The newly obtained RIs can supplement the currently used RIs.

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.

Genome-Wide Association Study of White Blood Cell Counts in Patients With Ischemic Stroke

Background and Purpose- Immune cells play a key role in the first 24h poststroke (acute phase), being associated with stroke outcome. We aimed to find genetic risk factors associated with leukocyte counts during the acute phase of stroke. Methods- Ischemic stroke patients with leukocyte counts data during the first 24h were included. Genome-wide association study and gene expression studies were performed. Results- Our genome-wide association study, which included 2064 (Discovery) and 407 (Replication) patients, revealed a new locus (14q24.3) associated with leukocyte counts. After Joint analysis (n=2471) 5 more polymorphisms reached genome-wide significance (P<5×10^-8). The 14q24.3 locus was associated with acute stroke outcome (rs112809786, P=0.036) and with ACOT1 and PTGR2 gene expression. Previous polymorphisms associated with leukocyte counts in general-population did not show any significance in our study. Conclusions- We have found the first locus associated with leukocyte counts in ischemic stroke, also associated with acute outcome. Genetic analysis of acute endophenotypes could be useful to find the genetic factors associated with stroke outcome. Our findings suggested a different modulation of immune cells in stroke compared with healthy conditions.

Generalizing polygenic risk scores from Europeans to Hispanics/Latinos

Polygenic risk scores (PRSs) are weighted sums of risk allele counts of single-nucleotide polymorphisms (SNPs) associated with a disease or trait. PRSs are typically constructed based on published results from Genome-Wide Association Studies (GWASs), and the majority of which has been performed in large populations of European ancestry (EA) individuals. Although many genotype-trait associations have generalized across populations, the optimal choice of SNPs and weights for PRSs may differ between populations due to different linkage disequilibrium (LD) and allele frequency patterns. We compare various approaches for PRS construction, using GWAS results from both large EA studies and a smaller study in Hispanics/Latinos: The Hispanic Community Health Study/Study of Latinos (HCHS/SOL, n = 12 , 803 ). We consider multiple approaches for selecting SNPs and for computing SNP weights. We study the performance of the resulting PRSs in an independent study of Hispanics/Latinos from the Women's Health Initiative (WHI, n = 3 , 582 ). We support our investigation with simulation studies of potential genetic architectures in a single locus. We observed that selecting variants based on EA GWASs generally performs well, except for blood pressure trait. However, the use of EA GWASs for weight estimation was suboptimal. Using non-EA GWAS results to estimate weights improved results.

Influence of tissue context on gene prioritization for predicted transcriptome-wide association studies

Transcriptome-wide association studies (TWAS) have recently gained great attention due to their ability to prioritize complex trait-associated genes and promote potential therapeutics development for complex human diseases. TWAS integrates genotypic data with expression quantitative trait loci (eQTLs) to predict genetically regulated gene expression components and associates predictions with a trait of interest. As such, TWAS can prioritize genes whose differential expressions contribute to the trait of interest and provide mechanistic explanation of complex trait(s). Tissue-specific eQTL information grants TWAS the ability to perform association analysis on tissues whose gene expression profiles are otherwise hard to obtain, such as liver and heart. However, as eQTLs are tissue context-dependent, whether and how the tissue-specificity of eQTLs influences TWAS gene prioritization has not been fully investigated. In this study, we addressed this question by adopting two distinct TWAS methods, PrediXcan and UTMOST, which assume single tissue and integrative tissue effects of eQTLs, respectively. Thirty-eight baseline laboratory traits in 4,360 antiretroviral treatment-naïve individuals from the AIDS Clinical Trials Group (ACTG) studies comprised the input dataset for TWAS. We performed TWAS in a tissue-specific manner and obtained a total of 430 significant gene-trait associations (q-value < 0.05) across multiple tissues. Single tissue-based analysis by PrediXcan contributed 116 of the 430 associations including 64 unique gene-trait pairs in 28 tissues. Integrative tissue-based analysis by UTMOST found the other 314 significant associations that include 50 unique gene-trait pairs across all 44 tissues. Both analyses were able to replicate some associations identified in past variant-based genome-wide association studies (GWAS), such as high-density lipoprotein (HDL) and CETP (PrediXcan, q-value = 3.2e-16). Both analyses also identified novel associations. Moreover, single tissue-based and integrative tissuebased analysis shared 11 of 103 unique gene-trait pairs, for example, PSRC1-low-density lipoprotein (PrediXcan's lowest q-value = 8.5e-06; UTMOST's lowest q-value = 1.8e-05). This study suggests that single tissue-based analysis may have performed better at discovering gene-trait associations when combining results from all tissues. Integrative tissue-based analysis was better at prioritizing genes in multiple tissues and in trait-related tissue. Additional exploration is needed to confirm this conclusion. Finally, although single tissue-based and integrative tissue-based analysis shared significant novel discoveries, tissue context-dependency of eQTLs impacted TWAS gene prioritization. This study provides preliminary data to support continued work on tissue contextdependency of eQTL studies and TWAS.

Off-target phenotypes in forensic DNA phenotyping and biogeographic ancestry inference: A resource

With recent advances in DNA sequencing technologies it has become feasible and cost effective to genotype larger marker sets for forensic purposes. Two technologies that make use of the larger marker sets have come into focus in forensic research and applications; inference of biogeographic ancestry (BGA) and forensic DNA phenotyping (FDP). These methods hold the promise to reveal information about a yet unknown perpetrator from a DNA sample. In contrast, DNA-profiling, that is a standard practice in case work, relies on matching DNA-profiles between crime scene material and suspects on a database of DNA-profiles. Markers for DNA-profiling were developed under the premise to reveal as little additional information about the human source of the profile as possible, the rationale being that personal privacy rights have to be balanced against the public interest in solving a crime. The same argument holds for markers used in BGA and FDP; these markers might also reveal information on off-target phenotypes (OTPs), that go beyond BGA and the phenotypes targeted in FDP. In particular, health related OTPs might shift the balance between privacy protection and public interest. However, to our knowledge, there is currently no convenient resource available to incorporate knowledge on OTPs in BGA and FDP assay design and application. In order to provide such a resource, we performed a systematic search for OTPs associated with a comprehensive set of markers (1766 SNPs) used or suggested to be used for BGA inference and FDP. In this set, we identified a relatively small number of 27 SNPs (1.53%) that convey information on diverse health related OTPs such as cancer risk, induced asthma, or risk of alcoholism. Some of these SNPs are commonly used for FDP and BGA across different marker sets. We conclude that the effects of SNP markers used in FDP and BGA on OTPs are currently limited, with few exceptions that should be considered in a balanced decision on assay design and application.