The Polygenic and Monogenic Basis of Blood Traits and Diseases

Genetic association between gut microbiome and blood pressure and blood cell count as mediator: A two-step Mendelian randomization analysis

BACKGROUND

Previous studies have established a genetic link between gut microbiota and hypertension, but whether blood cell count plays a mediating role in this remains unknown. This study aims to explore genetic associations and causal factors involving the gut microbiome, peripheral blood cell count, and blood pressure.

METHODS

We utilized summary statistics derived from genome-wide association studies to conduct a two-sample mediation Mendelian randomization analysis (https://gwas.mrcieu.ac.uk/). We applied inverse variance weighted (IVW) estimation method as the primary method, along with MR Egger, Weighted median, Simple mode and Weighted mode as complementary methods. To ensure the robustness of the results, several sensitivity analyses were conducted.

RESULTS

Genetic variants significantly associated with the microbiome, blood pressure, or peripheral blood cell counts were selected as instrumental variables. Fourteen microbial taxa were found to have suggestive associations with diastolic blood pressure (DBP), while fifteen microbial taxa showed suggestive associations with systolic blood pressure (SBP). Meanwhile, red blood cell count, lymphocyte count, and platelet count were identified to mediate the influence of the gut microbiome on blood pressure. Specifically, red cell count was identified to mediate the effects of the phylum Cyanobacteria on DBP (mediated proportion: 8.262 %). Lymphocyte count was found mediate the effects of the genus Subdoligranulum (mediated proportion: 2.642 %) and genus Collinsella (mediated proportion: 2.749 %) on SBP. Additionally, platelet count was found to mediate the relationship between the genus Eubacterium ventriosum group and SBP, explaining 3.421 % of the mediated proportion.

CONCLUSIONS

Our findings highlighted that gut microbiota may have causal influence on the blood pressure by modulating blood cell counts, which sheds new light on the pathogenesis and potential clinical interventions through the intricate axis of gut microbiome, blood cell counts, and blood pressure.

Single nucleotide polymorphisms in the cannabinoid CB2 receptor: Molecular pharmacology and disease associations

Preclinical evidence implicating cannabinoid receptor 2 (CB2) in various diseases has led researchers to question whether CB2 genetics influence aetiology or progression. Associations between conditions and genetic loci are often studied via single nucleotide polymorphism (SNP) prevalence in case versus control populations. In the CNR2 coding exon, ~36 SNPs have high overall population prevalence (minor allele frequencies [MAF] ~37%), including non-synonymous SNP (ns-SNP) rs2501432 encoding CB2 63Q/R. Interspersed are ~27 lower frequency SNPs, four being ns-SNPs. CNR2 introns also harbour numerous SNPs. This review summarises CB2 ns-SNP molecular pharmacology and evaluates evidence from ~70 studies investigating CB2 genetic variants with proposed linkage to disease. Although CNR2 genetic variation has been associated with a wide variety of conditions, including osteoporosis, immune-related disorders, and mental illnesses, further work is required to robustly validate CNR2 disease links and clarify specific mechanisms linking CNR2 genetic variation to disease pathophysiology and potential drug responses.

Abdominal aortic aneurysm and cardiometabolic traits share strong genetic susceptibility to lipid metabolism and inflammation

Abdominal aortic aneurysm has a high heritability and often co-occurs with other cardiometabolic disorders, suggesting shared genetic susceptibility. We investigate this commonality leveraging recent GWAS studies of abdominal aortic aneurysm and 32 cardiometabolic traits. We find significant genetic correlations between abdominal aortic aneurysm and 21 of the cardiometabolic traits investigated, including causal relationships with coronary artery disease, hypertension, lipid traits, and blood pressure. For each trait pair, we identify shared causal variants, genes, and pathways, revealing that cholesterol metabolism and inflammation are shared most prominently. Additionally, we show the tissue and cell type specificity in the shared signals, with strong enrichment across traits in the liver, arteries, adipose tissues, macrophages, adipocytes, and fibroblasts. Finally, we leverage drug-gene databases to identify several lipid-lowering drugs and antioxidants with high potential to treat abdominal aortic aneurysm with comorbidities. Our study provides insight into the shared genetic mechanism between abdominal aortic aneurysm and cardiometabolic traits, and identifies potential targets for pharmacological intervention.

Genetic analyses point to alterations in immune-related pathways underpinning the association between psychiatric disorders and COVID-19

Current literature suggests that people with psychiatric disorders have a higher risk of COVID-19 infection and a worse prognosis of the disease. We aimed to study the genetic contribution to these associations across seven psychiatric disorders as well as a general psychopathology factor (P-factor) and determine whether these are unique or shared across psychiatric disorders using statistical genetic techniques. Using the largest available genome-wide association studies (GWAS), we found a significant genetic overlap between depression, ADHD, PTSD, and the P-factor with both COVID-19 infection and hospitalization, and between anxiety and COVID-19 hospitalization. We used pairwise GWAS to examine this overlap on a fine-grained scale and identified specific regions of the genome shared between several psychiatric disorders, the P-factor, and COVID-19. Gene-based analysis in these genomic regions suggested possible links with immune-related pathways such as thyroid homeostasis, inflammation, and stress response. Finally, we show preliminary evidence for causal associations between depression, ADHD, PTSD, and the P-factor, and higher COVID-19 infection and hospitalization using Mendelian Randomization and Latent Causal Variable methods. Our results support the hypothesis that the relationship between psychiatric disorders and COVID-19 risk is likely due to shared alterations in immune-related pathways and is not a result of environmental factors alone, shedding light on potentially viable therapeutic targets.

LNK/SH2B3 loss of function increases susceptibility to murine and human atrial fibrillation

AIMS

The lymphocyte adaptor protein (LNK) is a negative regulator of cytokine and growth factor signalling. The rs3184504 variant in SH2B3 reduces LNK function and is linked to cardiovascular, inflammatory, and haematologic disorders, including stroke. In mice, deletion of Lnk causes inflammation and oxidative stress. We hypothesized that Lnk-/- mice are susceptible to atrial fibrillation (AF) and that rs3184504 is associated with AF and AF-related stroke in humans. During inflammation, reactive lipid dicarbonyls are the major components of oxidative injury, and we further hypothesized that these mediators are critical drivers of the AF substrate in Lnk-/- mice.

METHODS AND RESULTS

Lnk-/- or wild-type (WT) mice were treated with vehicle or 2-hydroxybenzylamine (2-HOBA), a dicarbonyl scavenger, for 3 months. Compared with WT, Lnk-/- mice displayed increased AF duration that was prevented by 2-HOBA. In the Lnk-/- atria, action potentials were prolonged with reduced transient outward K+ current, increased late Na+ current, and reduced peak Na+ current, pro-arrhythmic effects that were inhibited by 2-HOBA. Mitochondrial dysfunction, especially for Complex I, was evident in Lnk-/- atria, while scavenging lipid dicarbonyls prevented this abnormality. Tumour necrosis factor-α (TNF-α) and interleukin-1 beta (IL-1β) were elevated in Lnk-/- plasma and atrial tissue, respectively, both of which caused electrical and bioenergetic remodelling in vitro. Inhibition of soluble TNF-α prevented electrical remodelling and AF susceptibility, while IL-1β inhibition improved mitochondrial respiration but had no effect on AF susceptibility. In a large database of genotyped patients, rs3184504 was associated with AF, as well as AF-related stroke.

CONCLUSION

These findings identify a novel role for LNK in the pathophysiology of AF in both experimental mice and humans. Moreover, reactive lipid dicarbonyls are critical to the inflammatory AF substrate in Lnk-/- mice and mediate the pro-arrhythmic effects of pro-inflammatory cytokines, primarily through electrical remodelling.

A Novel Detection Method for High-Order SNP Epistatic Interactions Based on Explicit-Encoding-Based Multitasking Harmony Search

To elucidate the genetic basis of complex diseases, it is crucial to discover the single-nucleotide polymorphisms (SNPs) contributing to disease susceptibility. This is particularly challenging for high-order SNP epistatic interactions (HEIs), which exhibit small individual effects but potentially large joint effects. These interactions are difficult to detect due to the vast search space, encompassing billions of possible combinations, and the computational complexity of evaluating them. This study proposes a novel explicit-encoding-based multitasking harmony search algorithm (MTHS-EE-DHEI) specifically designed to address this challenge. The algorithm operates in three stages. First, a harmony search algorithm is employed, utilizing four lightweight evaluation functions, such as Bayesian network and entropy, to efficiently explore potential SNP combinations related to disease status. Second, a G-test statistical method is applied to filter out insignificant SNP combinations. Finally, two machine learning-based methods, multifactor dimensionality reduction (MDR) as well as random forest (RF), are employed to validate the classification performance of the remaining significant SNP combinations. This research aims to demonstrate the effectiveness of MTHS-EE-DHEI in identifying HEIs compared to existing methods, potentially providing valuable insights into the genetic architecture of complex diseases. The performance of MTHS-EE-DHEI was evaluated on twenty simulated disease datasets and three real-world datasets encompassing age-related macular degeneration (AMD), rheumatoid arthritis (RA), and breast cancer (BC). The results demonstrably indicate that MTHS-EE-DHEI outperforms four state-of-the-art algorithms in terms of both detection power and computational efficiency. The source code is available at https://github.com/shouhengtuo/MTHS-EE-DHEI.git .

Critically-ill COVID-19 susceptibility gene CCR3 shows natural selection in sub-Saharan Africans

The prevalence of COVID-19 critical illness varies across ethnicities, with recent studies suggesting that genetic factors may contribute to this variation. The aim of this study was to investigate natural selection signals of genes associated with critically-ill COVID-19 in sub-Saharan Africans. Severe COVID-19 SNPs were obtained from the HGI website. Selection signals were assessed in 661 sub-Sahara Africans from 1000 Genomes Project using integrated haplotype score (iHS), cross-population extended haplotype homozygosity (XP-EHH), and fixation index (Fst). Allele frequency trajectory analysis of ancient DNA samples were used to validate the existing of selection in sub-Sahara Africans. We also used Mendelian randomization to decipher the correlation between natural selection and critically-ill COVID-19. We identified that CCR3 exhibited significant natural selection signals in sub-Sahara Africans. Within the CCR3 gene, rs17217831-A showed both high iHS (Standardized iHS = 2) and high XP-EHH (Standardized XP-EHH = 2.5) in sub-Sahara Africans. Allele frequency trajectory of CCR3 rs17217831-A revealed natural selection occurring in the recent 1,500 years. Natural selection resulted in increased CCR3 expression in sub-Sahara Africans. Mendelian Randomization provided evidence that increased blood CCR3 expression and eosinophil counts lowered the risk of critically ill COVID-19. Our findings suggest that sub-Saharan Africans are resistant to critically ill COVID-19 due to natural selection and identify CCR3 as a potential novel therapeutic target.

Genetic drivers and cellular selection of female mosaic X chromosome loss

Mosaic loss of the X chromosome (mLOX) is the most common clonal somatic alteration in leukocytes of female individuals1,2, but little is known about its genetic determinants or phenotypic consequences. Here, to address this, we used data from 883,574 female participants across 8 biobanks; 12% of participants exhibited detectable mLOX in approximately 2% of leukocytes. Female participants with mLOX had an increased risk of myeloid and lymphoid leukaemias. Genetic analyses identified 56 common variants associated with mLOX, implicating genes with roles in chromosomal missegregation, cancer predisposition and autoimmune diseases. Exome-sequence analyses identified rare missense variants in FBXO10 that confer a twofold increased risk of mLOX. Only a small fraction of associations was shared with mosaic Y chromosome loss, suggesting that distinct biological processes drive formation and clonal expansion of sex chromosome missegregation. Allelic shift analyses identified X chromosome alleles that are preferentially retained in mLOX, demonstrating variation at many loci under cellular selection. A polygenic score including 44 allelic shift loci correctly inferred the retained X chromosomes in 80.7% of mLOX cases in the top decile. Our results support a model in which germline variants predispose female individuals to acquiring mLOX, with the allelic content of the X chromosome possibly shaping the magnitude of clonal expansion.

Obstructive Sleep Apnea Syndrome Exacerbates NASH Progression via Selective Autophagy-Mediated Eepd1 Degradation

Obstructive sleep apnea syndrome (OSAS), characterized by chronic intermittent hypoxia (CIH), is an independent risk factor for aggravating non-alcoholic steatohepatitis (NASH). The prevailing mouse model employed in CIH research is inadequate for the comprehensive exploration of the impact of CIH on NASH development due to reduced food intake observed in CIH-exposed mice, which deviates from human responses. To address this issue, a pair-feeding investigation with CIH-exposed and normoxia-exposed mice is conducted. It is revealed that CIH exposure aggravates DNA damage, leading to hepatic fibrosis and inflammation. The analysis of genome-wide association study (GWAS) data also discloses the association between Eepd1, a DNA repair enzyme, and OSAS. Furthermore, it is revealed that CIH triggered selective autophagy, leading to the autophagic degradation of Eepd1, thereby exacerbating DNA damage in hepatocytes. Notably, Eepd1 liver-specific knockout mice exhibit aggravated hepatic DNA damage and further progression of NASH. To identify a therapeutic approach for CIH-induced NASH, a drug screening is conducted and it is found that Retigabine dihydrochloride suppresses CIH-mediated Eepd1 degradation, leading to alleviated DNA damage in hepatocytes. These findings imply that targeting CIH-mediated Eepd1 degradation can be an adjunctive approach in the treatment of NASH exacerbated by OSAS.

Twin-based Mendelian Randomization Analyses Highlight Smoking's Effects on Blood DNA Methylation, with Putative Reverse Causation

Epigenome-wide association studies (EWAS) aim to identify differentially methylated loci associated with complex traits and disorders. EWAS of cigarette smoking shows some of the most widespread DNA methylation (DNAm) associations in blood. However, traditional EWAS cannot differentiate between causation and confounding, leading to ambiguity in etiological interpretations. Here, we apply an integrated approach combining Mendelian Randomization and twin-based Direction-of-Causation analyses (MR-DoC) to examine causality underlying smoking-associated blood DNAm changes in the Netherlands Twin Register (N=2577). Evidence across models suggests that current smoking's causal effects on DNAm likely drive many of the previous EWAS findings, implicating functional pathways relevant to several adverse health outcomes of smoking, including hemopoiesis, cell- and neuro-development, and immune regulation. Additionally, we find evidence of potential reverse causal influences at some DNAm sites, with 17 of these sites enriched for gene regulatory functional elements in the brain. The top three sites with evidence of DNAm's effects on smoking annotate to genes involved in G protein-coupled receptor signaling (GNG7, RGS3) and innate immune response (SLC15A4), elucidating potential biological risk factors for smoking. This study highlights the utility of integrating genotypic and DNAm measures in twin cohorts to clarify the causal relationships between health behaviors and blood DNAm.

CARD9 in the pathogenesis of axial spondyloarthritis

Axial spondyloarthritis (axSpA) has been long classified as an autoimmune disease caused by a breakdown in the ability of the immune system to delineate self from foreign, resulting in self-reactive T cells. The strong genetic association of HLA-B27 supports this role for T cells. More recently, genetic and clinical studies indicate a prominent role of the environment in triggering axSpA, including an important role for microbes and the innate immune response. As an example, mutations in genes associated with innate immunity, including the anti-fungal signaling molecule Caspase recruitment domain-containing protein 9 (CARD9), have been linked to axSpA susceptibility. Thus, current thought classifies axSpA as a "mixed pattern condition" caused by both autoimmune and autoinflammatory mechanisms. The goal of this review is to convey.

Fast and flexible joint fine-mapping of multiple traits via the Sum of Single Effects model

We introduce mvSuSiE, a multi-trait fine-mapping method for identifying putative causal variants from genetic association data (individual-level or summary data). mvSuSiE learns patterns of shared genetic effects from data, and exploits these patterns to improve power to identify causal SNPs. Comparisons on simulated data show that mvSuSiE is competitive in speed, power and precision with existing multi-trait methods, and uniformly improves on single-trait fine-mapping (SuSiE) in each trait separately. We applied mvSuSiE to jointly fine-map 16 blood cell traits using data from the UK Biobank. By jointly analyzing the traits and modeling heterogeneous effect sharing patterns, we discovered a much larger number of causal SNPs (>3,000) compared with single-trait fine-mapping, and with narrower credible sets. mvSuSiE also more comprehensively characterized the ways in which the genetic variants affect one or more blood cell traits; 68% of causal SNPs showed significant effects in more than one blood cell type.

SMIM1 absence is associated with reduced energy expenditure and excess weight

BACKGROUND

Obesity rates have nearly tripled in the past 50 years, and by 2030 more than 1 billion individuals worldwide are projected to be obese. This creates a significant economic strain due to the associated non-communicable diseases. The root cause is an energy expenditure imbalance, owing to an interplay of lifestyle, environmental, and genetic factors. Obesity has a polygenic genetic architecture; however, single genetic variants with large effect size are etiological in a minority of cases. These variants allowed the discovery of novel genes and biology relevant to weight regulation and ultimately led to the development of novel specific treatments.

METHODS

We used a case-control approach to determine metabolic differences between individuals homozygous for a loss-of-function genetic variant in the small integral membrane protein 1 (SMIM1) and the general population, leveraging data from five cohorts. Metabolic characterization of SMIM1-/- individuals was performed using plasma biochemistry, calorimetric chamber, and DXA scan.

FINDINGS

We found that individuals homozygous for a loss-of-function genetic variant in SMIM1 gene, underlying the blood group Vel, display excess body weight, dyslipidemia, altered leptin to adiponectin ratio, increased liver enzymes, and lower thyroid hormone levels. This was accompanied by a reduction in resting energy expenditure.

CONCLUSION

This research identified a novel genetic predisposition to being overweight or obese. It highlights the need to investigate the genetic causes of obesity to select the most appropriate treatment given the large cost disparity between them.

FUNDING

This work was funded by the National Institute of Health Research, British Heart Foundation, and NHS Blood and Transplant.

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

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

Genome-wide association study of the common retinal disorder epiretinal membrane: Significant risk loci in each of three American populations

Epiretinal membrane (ERM) is a common retinal condition characterized by the presence of fibrocellular tissue on the retinal surface, often with visual distortion and loss of visual acuity. We studied European American (EUR), African American (AFR), and Latino (admixed American, AMR) ERM participants in the Million Veteran Program (MVP) for genome-wide association analysis-a total of 38,232 case individuals and 557,988 control individuals. We completed a genome-wide association study (GWAS) in each population separately, and then results were meta-analyzed. Genome-wide significant (GWS) associations were observed in all three populations studied: 31 risk loci in EUR subjects, 3 in AFR, and 2 in AMR, with 48 in trans-ancestry meta-analysis. Many results replicated in the FinnGen sample. Several GWS variants associate to alterations in gene expression in the macula. ERM showed significant genetic correlation to multiple traits. Pathway enrichment analyses implicated collagen and collagen-adjacent mechanisms, among others. This well-powered ERM GWAS identified novel genetic associations that point to biological mechanisms for ERM.

Causal association between circulating blood cell traits and pulmonary embolism: a mendelian randomization study

BACKGROUND

Pulmonary embolism (PE) is a life-threatening thromboembolic disease for which there is limited evidence for effective prevention and treatment. Our goal was to determine whether genetically predicted circulating blood cell traits could influence the incidence of PE.

METHODS

Using single variable Mendelian randomization (SVMR) and multivariate Mendelian randomization (MVMR) analyses, we identified genetic associations between circulating blood cell counts and lymphocyte subsets and PE. GWAS blood cell characterization summary statistics were compiled from the Blood Cell Consortium. The lymphocyte subpopulation counts were extracted from summary GWAS statistics for samples from 3757 individuals that had been analyzed by flow cytometry. GWAS data related to PE were obtained from the FinnGen study.

RESULTS

According to the SVMR and reverse MR, increased levels of circulating white blood cells (odds ratio [OR]: 0.88, 95% confidence interval [CI]: 0.81-0.95, p = 0.0079), lymphocytes (OR: 0.90, 95% CI: 0.84-0.97, p = 0.0115), and neutrophils (OR: 0.88, 95% CI: 0.81-0.96, p = 0.0108) were causally associated with PE susceptibility. MVMR analysis revealed that lower circulating lymphocyte counts (OR: 0.84, 95% CI: 0.75-0.94, p = 0.0139) were an independent predictor of PE. According to further MR results, this association may be primarily related to HLA-DR+ natural killer (NK) cells.

CONCLUSIONS

Among European populations, there is a causal association between genetically predicted low circulating lymphocyte counts, particularly low HLA-DR+ NK cells, and an increased risk of PE. This finding supports observational studies that link peripheral blood cells to PE and provides recommendations for predicting and preventing this condition.

SpRY-mediated screens facilitate functional dissection of non-coding sequences at single-base resolution

CRISPR mutagenesis screens conducted with SpCas9 and other nucleases have identified certain cis-regulatory elements and genetic variants but at a limited resolution due to the absence of protospacer adjacent motif (PAM) sequences. Here, leveraging the broad targeting scope of the near-PAMless SpRY variant, we have demonstrated that saturated SpRY mutagenesis and base editing screens can faithfully identify functional regulatory elements and essential genetic variants for target gene expression at single-base resolution. We further extended this methodology to investigate a genome-wide association study (GWAS) locus at 10q22.1 associated with a red blood cell trait, where we identified potential enhancers regulating HK1 gene expression, despite not all of these enhancers exhibiting typical chromatin signatures. More importantly, our saturated base editing screens pinpoint multiple causal variants within this locus that would otherwise be missed by Bayesian statistical fine-mapping. Our approach is generally applicable to functional interrogation of all non-coding genomic elements while complementing other high-coverage CRISPR screens.

Oral and non-oral lichen planus show genetic heterogeneity and differential risk for autoimmune disease and oral cancer

Lichen planus (LP) is a T-cell-mediated inflammatory disease affecting squamous epithelia in many parts of the body, most often the skin and oral mucosa. Cutaneous LP is usually transient and oral LP (OLP) is most often chronic, so we performed a large-scale genetic and epidemiological study of LP to address whether the oral and non-oral subgroups have shared or distinct underlying pathologies and their overlap with autoimmune disease. Using lifelong records covering diagnoses, procedures, and clinic identity from 473,580 individuals in the FinnGen study, genome-wide association analyses were conducted on carefully constructed subcategories of OLP (n = 3,323) and non-oral LP (n = 4,356) and on the combined group. We identified 15 genome-wide significant associations in FinnGen and an additional 12 when meta-analyzed with UKBB (27 independent associations at 25 distinct genomic locations), most of which are shared between oral and non-oral LP. Many associations coincide with known autoimmune disease loci, consistent with the epidemiologic enrichment of LP with hypothyroidism and other autoimmune diseases. Notably, a third of the FinnGen associations demonstrate significant differences between OLP and non-OLP. We also observed a 13.6-fold risk for tongue cancer and an elevated risk for other oral cancers in OLP, in agreement with earlier reports that connect LP with higher cancer incidence. In addition to a large-scale dissection of LP genetics and comorbidities, our study demonstrates the use of comprehensive, multidimensional health registry data to address outstanding clinical questions and reveal underlying biological mechanisms in common but understudied diseases.

Genetic Evidence for a Causal Relationship Between Innate Leukocytes and the Risk of Digestive System Cancers in East Asians and Europeans

Background

Peripheral traditional immune cell disorder plays an important role in cancer onset and development. The causal relationships between leukocytes prior to cancer and the risk of digestive system cancer remain unknown. This study assesses the causal correlations between leukocytes and digestive system cancer risk in East Asians and Europeans.

Methods

Summary-level data on leukocyte-related genetic variation were extracted from Biobank Japan (107,964 participants) and a recent large-scale meta-analysis (563,946 participants). Summary-level data for the cancers were obtained from Biobank Japan (212,978 individuals) and the FinnGen consortium (178,802 participants). Univariable and multivariable Mendelian randomization (MR) analyses were performed on East Asians and Europeans separately.

Results

Univariable MR analysis demonstrated the significant association between circulating eosinophil counts and risk of colorectal cancer (CRC) in East Asians (odds ratio (OR) = 0.80, 95% confidence interval (CI): 0.69 - 0.92, P = 0.002) and a suggestive relationship in the European population (OR = 0.86, 95% CI: 0.77 - 0.97, P = 0.013). An inverse suggestive association was observed between levels of basophils and the risk of gastric cancer (GC) in East Asians (OR = 0.83, 95% CI: 0.72 - 0.97, P = 0.019). The multivariable MR analysis showed the independent causal effect of eosinophil count on CRC risk in East Asians (OR = 0.72, 95% CI: 0.57 - 0.92, P = 0.009) and Europeans (OR = 0.80, 95% CI: 0.70 - 0.92, P = 0.002). Circulating basophils served as the negative causal factor in GC risk in East Asians (OR = 0.80, 95% CI: 0.67 - 0.94, P = 0.007).

Conclusions

Our MR analyses revealed a genetic causal relationship between reduced blood eosinophils and an increased CRC risk in both Europeans and East Asians. Furthermore, our results suggested a causal association between decreased basophils and an elevated GC risk specifically in East Asians.

No evidence linking sleep traits with white blood cell counts: Multivariable-adjusted and Mendelian randomization analyses

BACKGROUND

Disturbances in habitual sleep have been associated with multiple age-associated diseases. However, the biological mechanisms underpinning these associations remain largely unclear. We assessed the possible involvement of the circulating immune system by determining the associations between sleep traits and white blood cell counts using multivariable-adjusted linear regression and Mendelian randomization.

METHODS

Cross-sectional multivariable-adjusted linear regression analyses were done using participants within the normal range of total white blood cell counts (>4.5 × 109 and <11.0 × 109/μL) from UK Biobank. For the sleep traits, we examined (short and long) sleep duration, chronotype, insomnia symptoms and daytime dozing. Two-sample Mendelian randomization analyses were done using instruments for sleep traits derived from European-ancestry participants from UK Biobank (over 410,000 participants) and using SNP-outcome data derived from European-ancestry participants from the Blood Cell Consortium (N = 563,946) to which no data from UK Biobank contributed.

RESULTS

Using data from 357,656 participants (mean [standard deviation] age: 56.5 [8.1] years, and 44.4% men), we did not find evidence that disturbances in any of the studied sleep traits were associated with differences in blood cell counts (total, lymphocytes, neutrophiles, eosinophiles and basophiles). Also, we did not find associations between disturbances in any of the studied sleep traits and white blood cell counts using Mendelian Randomization.

CONCLUSION

Based on the results from two different methodologies, disturbances in habitual sleep are unlikely to cause changes in blood cell counts and thereby differences in blood cell counts are unlikely to be underlying the observed sleep-disease associations.

The pivotal role of the X-chromosome in the genetic architecture of the human brain

Genes on the X-chromosome are extensively expressed in the human brain. However, little is known for the X-chromosome's impact on the brain anatomy, microstructure, and functional network. We examined 1,045 complex brain imaging traits from 38,529 participants in the UK Biobank. We unveiled potential autosome-X-chromosome interactions, while proposing an atlas outlining dosage compensation (DC) for brain imaging traits. Through extensive association studies, we identified 72 genome-wide significant trait-locus pairs (including 29 new associations) that share genetic architectures with brain-related disorders, notably schizophrenia. Furthermore, we discovered unique sex-specific associations and assessed variations in genetic effects between sexes. Our research offers critical insights into the X-chromosome's role in the human brain, underscoring its contribution to the differences observed in brain structure and functionality between sexes.

Mendelian randomization analysis reveals causal relationships between circulating cell traits and renal disorders

Purpose

The aim of this study was to investigate the causal relationships between circulating cell traits and risk of renal disorders.

Methods

We applied a comprehensive two-sample Mendelian randomization (MR) analysis. Single nucleotide polymorphisms (SNPs) from publicly available genome-wide association studies (GWAS) databases were utilized. Genetically predicted instrumental variables of human blood cell traits were extracted from Blood Cell Consortium (BCX) while data on renal diseases was obtained from Finngen consortium. The primary MR analysis was conducted using the inverse variance weighted (IVW) method, with the weighted median (WM) and MR-Egger models used as additional methods. Sensitivity analyses, including MR-PRESSO, radial regression and MR-Egger intercept were conducted to detect outliers and assess horizontal pleiotropy. We further utilized the leave-one-out analysis to assess the robustness of the results. Causal associations were considered significant based on false rate correction (FDR), specifically when the IVW method provided a pFDR < 0.05.

Results

Our results demonstrated that both white blood cell (WBC) count (OR = 1.50, 95% CI = 1.10-2.06, pFDR = 0.033, pIVW = 0.011) and lymphocyte count (OR = 1.50, 95% CI = 1.13-1.98, pFDR = 0.027, pIVW = 0.005) were causally associated with a higher risk of IgA nephropathy. Furthermore, WBC count was identified as a significant genetic risk factor for renal malignant neoplasms (OR = 1.23, 95% CI = 1.06-1.43, pFDR = 0.041, pIVW = 0.007). Additionally, an increased level of genetically predicted eosinophils was found to be causally associated with a higher risk of diabetic nephropathy (OR = 1.21, 95% CI = 1.08-1.36, pFDR = 0.007, pIVW = 0.001). No evidence of pleiotropy was determined.

Conclusion

Our findings provide evidence of causal associations of circulating WBC count, lymphocyte count and IgA nephropathy, WBC count and renal malignant neoplasms, and eosinophil count and diabetic nephropathy. These results have the potential to contribute to the development of novel diagnostic options and therapeutic strategies for renal disorders.

Improving polygenic prediction from summary data by learning patterns of effect sharing across multiple phenotypes

Polygenic prediction of complex trait phenotypes has become important in human genetics, especially in the context of precision medicine. Recently, Morgante et al. introduced mr.mash, a flexible and computationally efficient method that models multiple phenotypes jointly and leverages sharing of effects across such phenotypes to improve prediction accuracy. However, a drawback of mr.mash is that it requires individual-level data, which are often not publicly available. In this work, we introduce mr.mash-rss, an extension of the mr.mash model that requires only summary statistics from Genome-Wide Association Studies (GWAS) and linkage disequilibrium (LD) estimates from a reference panel. By using summary data, we achieve the twin goal of increasing the applicability of the mr.mash model to data sets that are not publicly available and making it scalable to biobank-size data. Through simulations, we show that mr.mash-rss is competitive with, and often outperforms, current state-of-the-art methods for single- and multi-phenotype polygenic prediction in a variety of scenarios that differ in the pattern of effect sharing across phenotypes, the number of phenotypes, the number of causal variants, and the genomic heritability. We also present a real data analysis of 16 blood cell phenotypes in UK Biobank, showing that mr.mash-rss achieves higher prediction accuracy than competing methods for the majority of traits, especially when the data has smaller sample size.

L-plastin associated syndrome of immune deficiency and hematologic cytopenia

BACKGROUND

LCP1 encodes L-plastin, an actin-bundling protein primarily expressed in hematopoietic cells. In mouse and fish models, LCP1 deficiency has been shown to result in hematologic and immune defects.

OBJECTIVE

This study aimed to determine the nature of a human inborn error of immunity resulting from a novel genetic variant of LCP1.

METHODS

We performed genetic, protein, and cellular analysis of PBMCs from a kindred with apparent autosomal dominant immune deficiency. We identified a candidate causal mutation in LCP1, which we evaluated by engineering the orthologous mutation in mice and Jurkat cells.

RESULTS

A splice-site variant in LCP1 segregated with lymphopenia, neutropenia, and thrombocytopenia. The splicing defect resulted in at least 2 aberrant transcripts, producing an in-frame deletion of 24 nucleotides, and a frameshift deletion of exon 8. Cellular analysis of the kindred revealed a proportionate reduction of T and B cells and a mild expansion of transitional B cells. Similarly, mice carrying the orthologous genetic variant exhibited the same in-frame aberrant transcript, reduced expression Lcp1 and gene dose-dependent leukopenia, mild thrombocytopenia, and lymphopenia, with a significant reduction of T-cell populations. Functional analysis revealed that LCP1c740-1G>A confers a defect in platelet development and function with aberrant spreading on collagen. Immunologic analysis revealed defective actin organization in T cells, reduced migration of PBMCs from patients, splenocytes from mutant mice, and a mutant Jurkat cell line in response to CXCL12; impaired germinal center B-cell expansion after immunization; and reduced cytokinesis during T cell proliferation.

CONCLUSIONS

We describe a unique human hematopoietic defect affecting neutrophils, lymphocytes, and platelets arising from partial LCP1 deficiency.

Genetic association of inflammatory marker GlycA with lung function and respiratory diseases

Association of circulating glycoprotein acetyls (GlycA), a systemic inflammation biomarker, with lung function and respiratory diseases remain to be investigated. We examined the genetic correlation, shared genetics, and potential causality of GlycA (N = 115,078) with lung function and respiratory diseases (N = 497,000). GlycA showed significant genetic correlation with FEV1 (rg = -0.14), FVC (rg = -0.18), asthma (rg = 0.21) and COPD (rg = 0.31). We consistently identified ten shared loci (including chr3p21.31 and chr8p23.1) at both SNP and gene level revealing potential shared biological mechanisms involving ubiquitination, immune response, Wnt/β-catenin signaling, cell growth and differentiation in tissues or cells including blood, epithelium, fibroblast, fetal thymus, and fetal intestine. Genetically elevated GlycA was significantly correlated with lung function and asthma susceptibility (354.13 ml decrement of FEV1, 442.28 ml decrement of FVC, and 144% increased risk of asthma per SD increment of GlycA) from MR analyses. Our findings provide insights into biological mechanisms of GlycA in relating to lung function, asthma, and COPD.

Revisiting the association of sudden infant death syndrome (SIDS) with polymorphisms of NHE3 and IL13

OBJECTIVES

Disturbances of the central nervous system and immune system are thought to play a role in sudden infant death syndrome (SIDS). Dysregulated expression of sodium (Na+)/hydrogen (H+) exchanger 3 (NHE3) in the brainstem and of interleukin 13 (IL13) in the lungs has been observed in SIDS. An association of single-nucleotide polymorphisms (SNPs) in NHE3 and IL13 with SIDS has been proposed, but controversial results were reported. Therefore, there is a need to revisit the association of SNPs in NHE3 and IL13 with SIDS.

METHODS

Genotyping of rs71597645 (G1131A) and rs2247114 (C2405T) in NHE3 and rs20541 (+ 4464A/G) in IL13 was performed in 201 SIDS cases and 338 controls. A meta-analysis was performed after merging our data with previously published data (all from European populations).

RESULTS

Polymorphisms rs2247114 (NHE3) and rs20541 (IL13) were significantly associated with SIDS overall and in multiple subgroups, but no association was found for rs71597645 (NHE3). After combining our data with previously published data, a fixed-effect meta-analysis showed that rs2247114 in NHE3 retained a significant association with SIDS under a recessive model (OR 2.78, 95%CI 1.53 to 5.06; p = 0.0008).

CONCLUSION

Our findings suggest an association of NHE3 variant rs2247114 (C2405T), though not rs71597645 (NHE3), with SIDS. A potential role of rs20541 (IL13) still has to be elucidated. Especially NHE3 seems to be an interesting topic for future SIDS research.

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.

Germline genetic regulation of the colorectal tumor immune microenvironment

OBJECTIVE

To evaluate the contribution of germline genetics to regulating the briskness and diversity of T cell responses in CRC, we conducted a genome-wide association study to examine the associations between germline genetic variation and quantitative measures of T cell landscapes in 2,876 colorectal tumors from participants in the Molecular Epidemiology of Colorectal Cancer Study (MECC).

METHODS

Germline DNA samples were genotyped and imputed using genome-wide arrays. Tumor DNA samples were extracted from paraffin blocks, and T cell receptor clonality and abundance were quantified by immunoSEQ (Adaptive Biotechnologies, Seattle, WA). Tumor infiltrating lymphocytes per high powered field (TILs/hpf) were scored by a gastrointestinal pathologist. Regression models were used to evaluate the associations between each variant and the three T-cell features, adjusting for sex, age, genotyping platform, and global ancestry. Three independent datasets were used for replication.

RESULTS

We identified a SNP (rs4918567) near RBM20 associated with clonality at a genome-wide significant threshold of 5 × 10- 8, with a consistent direction of association in both discovery and replication datasets. Expression quantitative trait (eQTL) analyses and in silico functional annotation for these loci provided insights into potential functional roles, including a statistically significant eQTL between the T allele at rs4918567 and higher expression of ADRA2A (P = 0.012) in healthy colon mucosa.

CONCLUSIONS

Our study suggests that germline genetic variation is associated with the quantity and diversity of adaptive immune responses in CRC. Further studies are warranted to replicate these findings in additional samples and to investigate functional genomic mechanisms.

Examining the causal relationship between circulating immune cells and the susceptibility to osteomyelitis: A Mendelian randomization study

BACKGROUND

Osteomyelitis is considered as a deleterious inflammatory condition affecting the bone, primarily attributed to pathogenic infection. However, the underlying factors predisposing individuals to osteomyelitis remain incompletely elucidated. The immune system plays a multifaceted role in the progression of this condition, yet previous observational studies and randomized controlled trials investigating the association between circulating immune cell counts and osteomyelitis have been constrained. In order to address this knowledge gap, we conducted a Mendelian randomization (MR) analysis to evaluate the impact of diverse immune cell counts on the risk of developing osteomyelitis.

METHODS

In our study, we utilized single nucleotide polymorphisms (SNPs) that have been strongly linked to circulating immune cells or specific lymphocyte subtypes, as identified in large-scale genome-wide association studies (GWAS). These SNPs served as instrumental variables (IVs) for our MR analysis. We employed a more relaxed clumping threshold to conduct MR analysis on several related lymphocyte subtypes. To estimate causal effects, we utilized the Wald ratio, as well as the random-effects inverse variance weighted (IVW) and weighted median (WM) methods. To enhance the credibility of our results, we performed F-statistic calculations and a series of sensitivity analyses.

RESULTS

Our findings revealed a significant correlation between the absolute count of circulating lymphocytes and the risk of osteomyelitis [odds ratio(OR) 1.20；95 % confidence interval (CI), 1.08-1.32；P = 0.0005]. Furthermore, we identified a causal relationship between the absolute count of CD8+ T cells and susceptibility to osteomyelitis (OR 1.16; 95 % CI, 1.04-1.30; P = 0.0098). Importantly, these findings remained robust across a wide range of sensitivity analyses.

CONCLUSION

Through our MR analysis, we have provided evidence supporting a causal relationship between genetic predisposition to higher circulating immune cell counts and an increased risk of osteomyelitis. Specifically, our findings highlight the association between elevated CD8+ T cell counts and a heightened susceptibility to osteomyelitis. These results offer valuable insights for the future exploration of immunotherapy approaches in the management of osteomyelitis.

Genome-wide analyses of variance in blood cell phenotypes provide new insights into complex trait biology and prediction

Blood cell phenotypes are routinely tested in healthcare to inform clinical decisions. Genetic variants influencing mean blood cell phenotypes have been used to understand disease aetiology and improve prediction; however, additional information may be captured by genetic effects on observed variance. Here, we mapped variance quantitative trait loci (vQTL), i.e. genetic loci associated with trait variance, for 29 blood cell phenotypes from the UK Biobank (N~408,111). We discovered 176 independent blood cell vQTLs, of which 147 were not found by additive QTL mapping. vQTLs displayed on average 1.8-fold stronger negative selection than additive QTL, highlighting that selection acts to reduce extreme blood cell phenotypes. Variance polygenic scores (vPGSs) were constructed to stratify individuals in the INTERVAL cohort (N~40,466), where genetically less variable individuals (low vPGS) had increased conventional PGS accuracy (by ~19%) than genetically more variable individuals. Genetic prediction of blood cell traits improved by ~10% on average combining PGS with vPGS. Using Mendelian randomisation and vPGS association analyses, we found that alcohol consumption significantly increased blood cell trait variances highlighting the utility of blood cell vQTLs and vPGSs to provide novel insight into phenotype aetiology as well as improve prediction.

Genetic basis of pregnancy-associated decreased platelet counts and gestational thrombocytopenia

ABSTRACT

Platelet count reduction occurs throughout pregnancy, with 5% to 12% of pregnant women being diagnosed with gestational thrombocytopenia (GT), characterized by a more marked decrease in platelet count during pregnancy. However, the underlying biological mechanism behind these phenomena remains unclear. Here, we used sequencing data from noninvasive prenatal testing of 100 186 Chinese pregnant individuals and conducted, to our knowledge, the hitherto largest-scale genome-wide association studies on platelet counts during 5 periods of pregnancy (the first, second, and third trimesters, delivery, and the postpartum period) as well as 2 GT statuses (GT platelet count < 150 × 109/L and severe GT platelet count < 100 × 109/L). Our analysis revealed 138 genome-wide significant loci, explaining 10.4% to 12.1% of the observed variation. Interestingly, we identified previously unknown changes in genetic effects on platelet counts during pregnancy for variants present in PEAR1 and CBL, with PEAR1 variants specifically associated with a faster decline in platelet counts. Furthermore, we found that variants present in PEAR1 and TUBB1 increased susceptibility to GT and severe GT. Our study provides insight into the genetic basis of platelet counts and GT in pregnancy, highlighting the critical role of PEAR1 in decreasing platelet counts during pregnancy and the occurrence of GT. Those with pregnancies carrying specific variants associated with declining platelet counts may experience a more pronounced decrease, thereby elevating the risk of GT. These findings lay the groundwork for further investigation into the biological mechanisms and causal implications of GT.

Lipoprotein lipids and apolipoproteins in primary immune thrombocytopenia: Results from a clinical characteristics and causal relationship verification, potential drug target identification by Mendelian randomization analyses

Primary immune thrombocytopenia (ITP) is an acquired autoimmune disease. Cellular and systemic lipid metabolism plays a significant role in the regulation of immune cell activities. However, the role of lipoprotein lipids and apolipoproteins in ITP remains elusive. The automatic biochemistry analyser was used to measure the levels of serum total cholesterol (TC), triglyceride (TG), low-density lipoprotein cholesterol (LDL-C), high-density lipoprotein cholesterol (HDL-C), apolipoprotein A-I (apoA-I), apoB, apoE and lipoprotein a [LP(a)]. Genetic variants strongly associated with circulating lipoprotein lipids and apolipoproteins (LDL-C, apoB, TG, HDL-C and apoA-I) were extracted to perform Mendelian randomization (MR) analyses. Finally, drug-target MR and passive ITP mice model was used to investigate the potential druggable targets of ITP. Levels of HDL-C, apoA-I, decreased and LP(a) increased in ITP patients compared with healthy controls. Low HDL-C was causally associated with ITP susceptibility. Through drug-target MR and animal modelling, ABCA1 was identified as a potential target to design drugs for ITP. Our study found that lipid metabolism is related to ITP. The causative association between HDL-C and the risk of ITP was also established. The study provided new evidence of the aetiology of ITP. ABCA1 might be a potential drug target for ITP.

rs2736098, a synonymous polymorphism, is associated with carcinogenesis and cell count in multiple tissue types by regulating TERT expression

rs2736098 is a synonymous polymorphism in TERT (telomerase reverse transcriptase), an enzyme involved in tumor onset of multiple tissues, and should play no roles in carcinogenesis. However, a search in cancer somatic mutation database indicated that the mutation frequency at rs2736098 is much higher than the average one for TERT. Moreover, there are significant H3K4me1 and H3K27Ac signals, two universal histone modifications for active enhancers, surrounding rs2736098. Therefore, we hypothesized that rs2736098 might be within an enhancer region, regulate TERT expression and influence cancer risk. Through luciferase assay, it was verified that the enhancer activity of rs2736098C allele is significantly higher than that of T in multiple tissues. Transfection of plasmids containing TERT coding region with two different alleles indicated that rs2736098C allele can induce a significantly higher TERT expression than T. By chromatin immunoprecipitation, it was observed that the fragment spanning rs2736098 can interact with USF1 (upstream transcription factor 1). The two alleles of rs2736098 present evidently different binding affinity with nuclear proteins. Database and literature search indicated that rs2736098 is significantly associated with carcinogenesis in multiple tissues and count of multiple cell types. All these facts indicated that rs2736098 is also an oncogenic polymorphism and plays important role in cell proliferation.

Peripheral blood indicators and COVID-19: an observational and bidirectional Mendelian randomization study

Blood is critical for health, supporting key functions like immunity and oxygen transport. While studies have found links between common blood clinical indicators and COVID-19, they cannot provide causal inference due to residual confounding and reverse causality. To identify indicators affecting COVID-19, we analyzed clinical data (n = 2,293, aged 18-65 years) from Guangzhou Medical University's first affiliated hospital (2022-present), identifying 34 significant indicators differentiating COVID-19 patients from healthy controls. Utilizing bidirectional Mendelian randomization analyses, integrating data from over 2.46 million participants from various large-scale studies, we established causal links for six blood indicators with COVID-19 risk, five of which is consistent with our observational findings. Specifically, elevated Troponin I and Platelet Distribution Width levels are linked with increased COVID-19 susceptibility, whereas higher Hematocrit, Hemoglobin, and Neutrophil counts confer a protective effect. Reverse MR analysis confirmed four blood biomarkers influenced by COVID-19, aligning with our observational data for three of them. Notably, COVID-19 exhibited a positive causal relationship with Troponin I (Tnl) and Serum Amyloid Protein A, while a negative association was observed with Plateletcrit. These findings may help identify high-risk individuals and provide further direction on the management of COVID-19.

Population-enriched innate immune variants may identify candidate gene targets at the intersection of cancer and cardio-metabolic disease

Both cancer and cardio-metabolic disease disparities exist among specific populations in the US. For example, African Americans experience the highest rates of breast and prostate cancer mortality and the highest incidence of obesity. Native and Hispanic Americans experience the highest rates of liver cancer mortality. At the same time, Pacific Islanders have the highest death rate attributed to type 2 diabetes (T2D), and Asian Americans experience the highest incidence of non-alcoholic fatty liver disease (NAFLD) and cancers induced by infectious agents. Notably, the pathologic progression of both cancer and cardio-metabolic diseases involves innate immunity and mechanisms of inflammation. Innate immunity in individuals is established through genetic inheritance and external stimuli to respond to environmental threats and stresses such as pathogen exposure. Further, individual genomes contain characteristic genetic markers associated with one or more geographic ancestries (ethnic groups), including protective innate immune genetic programming optimized for survival in their corresponding ancestral environment(s). This perspective explores evidence related to our working hypothesis that genetic variations in innate immune genes, particularly those that are commonly found but unevenly distributed between populations, are associated with disparities between populations in both cancer and cardio-metabolic diseases. Identifying conventional and unconventional innate immune genes that fit this profile may provide critical insights into the underlying mechanisms that connect these two families of complex diseases and offer novel targets for precision-based treatment of cancer and/or cardio-metabolic disease.

Structural and genetic diversity in the secreted mucins, MUC5AC and MUC5B

The secreted mucins MUC5AC and MUC5B play critical defensive roles in airway pathogen entrapment and mucociliary clearance by encoding large glycoproteins with variable number tandem repeats (VNTRs). These polymorphic and degenerate protein coding VNTRs make the loci difficult to investigate with short reads. We characterize the structural diversity of MUC5AC and MUC5B by long-read sequencing and assembly of 206 human and 20 nonhuman primate (NHP) haplotypes. We find that human MUC5B is largely invariant (5761-5762aa); however, seven haplotypes have expanded VNTRs (6291-7019aa). In contrast, 30 allelic variants of MUC5AC encode 16 distinct proteins (5249-6325aa) with cysteine-rich domain and VNTR copy number variation. We grouped MUC5AC alleles into three phylogenetic clades: H1 (46%, ~5654aa), H2 (33%, ~5742aa), and H3 (7%, ~6325aa). The two most common human MUC5AC variants are smaller than NHP gene models, suggesting a reduction in protein length during recent human evolution. Linkage disequilibrium (LD) and Tajima's D analyses reveal that East Asians carry exceptionally large MUC5AC LD blocks with an excess of rare variation (p<0.05). To validate this result, we used Locityper for genotyping MUC5AC haplogroups in 2,600 unrelated samples from the 1000 Genomes Project. We observed signatures of positive selection in H1 and H2 among East Asians and a depletion of the likely ancestral haplogroup (H3). In Africans and Europeans, H3 alleles show an excess of common variation and deviate from Hardy-Weinberg equilibrium, consistent with heterozygote advantage and balancing selection. This study provides a generalizable strategy to characterize complex protein coding VNTRs for improved disease associations.

Rare copy number variation in autoimmune Addison's disease

Autoimmune Addison's disease (AAD) is a rare but life-threatening endocrine disorder caused by an autoimmune destruction of the adrenal cortex. A previous genome-wide association study (GWAS) has shown that common variants near immune-related genes, which mostly encode proteins participating in the immune response, affect the risk of developing this condition. However, little is known about the contribution of copy number variations (CNVs) to AAD susceptibility. We used the genome-wide genotyping data from Norwegian and Swedish individuals (1,182 cases and 3,810 controls) to investigate the putative role of CNVs in the AAD aetiology. Although the frequency of rare CNVs was similar between cases and controls, we observed that larger deletions (>1,000 kb) were more common among patients (OR = 4.23, 95% CI 1.85-9.66, p = 0.0002). Despite this, none of the large case-deletions were conclusively pathogenic, and the clinical presentation and an AAD-polygenic risk score were similar between cases with and without the large CNVs. Among deletions exclusive to individuals with AAD, we highlight two ultra-rare deletions in the genes LRBA and BCL2L11, which we speculate might have contributed to the polygenic risk in these carriers. In conclusion, rare CNVs do not appear to be a major cause of AAD but further studies are needed to ascertain the potential contribution of rare deletions to the polygenic load of AAD susceptibility.

Effect of the gut microbiome, plasma metabolome, peripheral cells, and inflammatory cytokines on obesity: a bidirectional two-sample Mendelian randomization study and mediation analysis

Background

Obesity is a metabolic and chronic inflammatory disease involving genetic and environmental factors. This study aimed to investigate the causal relationship among gut microbiota abundance, plasma metabolomics, peripheral cell (blood and immune cell) counts, inflammatory cytokines, and obesity.

Methods

Summary statistics of 191 gut microbiota traits (N = 18,340), 1,400 plasma metabolite traits (N = 8,299), 128 peripheral cell counts (blood cells, N = 408,112; immune cells, N = 3,757), 41 inflammatory cytokine traits (N = 8,293), and 6 obesity traits were obtained from publicly available genome-wide association studies. Two-sample Mendelian randomization (MR) analysis was applied to infer the causal links using inverse variance-weighted, maximum likelihood, MR-Egger, weighted median, weighted mode, and Wald ratio methods. Several sensitivity analyses were also utilized to ensure reliable MR results. Finally, we used mediation analysis to identify the pathway from gut microbiota to obesity mediated by plasma metabolites, peripheral cells, and inflammatory cytokines.

Results

MR revealed a causal effect of 44 gut microbiota taxa, 281 plasma metabolites, 27 peripheral cells, and 8 inflammatory cytokines on obesity. Among them, five shared causal gut microbiota taxa belonged to the phylum Actinobacteria, order Bifidobacteriales, family Bifidobacteriaceae, genus Lachnospiraceae UCG008, and species Eubacterium nodatum group. Furthermore, we screened 42 shared causal metabolites, 7 shared causal peripheral cells, and 1 shared causal inflammatory cytokine. Based on known causal metabolites, we observed that the metabolic pathways of D-arginine, D-ornithine, linoleic acid, and glycerophospholipid metabolism were closely related to obesity. Finally, mediation analysis revealed 20 mediation relationships, including the causal pathway from gut microbiota to obesity, mediated by 17 metabolites, 2 peripheral cells, and 1 inflammatory cytokine. Sensitivity analysis represented no heterogeneity or pleiotropy in this study.

Conclusion

Our findings support a causal relationship among gut microbiota, plasma metabolites, peripheral cells, inflammatory cytokines, and obesity. These biomarkers provide new insights into the mechanisms underlying obesity and contribute to its prevention, diagnosis, and treatment.

Epidemiologic and genetic associations of female reproductive disorders with depression or dysthymia: a Mendelian randomization study

Observational studies have previously reported an association between depression and certain female reproductive disorders. However, the causal relationships between depression and different types of female reproductive disorders remain unclear in terms of direction and magnitude. We conducted a comprehensive investigation using a two-sample bi-directional Mendelian randomization analysis, incorporating publicly available GWAS summary statistics. Our aim was to establish a causal relationship between genetically predicted depression and the risk of various female reproductive pathological conditions, such as ovarian dysfunction, polycystic ovary syndrome(PCOS), ovarian cysts, abnormal uterine and vaginal bleeding(AUB), endometriosis, leiomyoma of the uterus, female infertility, spontaneous abortion, eclampsia, pregnancy hypertension, gestational diabetes, excessive vomiting in pregnancy, cervical cancer, and uterine/endometrial cancer. We analyzed a substantial sample size, ranging from 111,831 to 210,870 individuals, and employed robust statistical methods, including inverse variance weighted, MR-Egger, weighted median, and MR-PRESSO, to estimate causal effects. Sensitivity analyses, such as Cochran's Q test, MR-Egger intercept test, MR-PRESSO, leave-one-out analysis, and funnel plots, were also conducted to ensure the validity of our results. Furthermore, risk factor analyses were performed to investigate potential mediators associated with these observed relationships. Our results demonstrated that genetic predisposition to depression or dysthymia was associated with an increased risk of developing PCOS (OR = 1.43, 95% CI 1.28-1.59; P = 6.66 × 10-11), ovarian cysts (OR = 1.36, 95% CI 1.20-1.55; P = 1.57 × 10-6), AUB (OR = 1.41, 95% CI 1.20-1.66; P = 3.01 × 10-5), and endometriosis (OR = 1.43, 95% CI 1.27-1.70; P = 2.21 × 10-7) after Bonferroni correction, but no evidence for reverse causality. Our study did not find any evidence supporting a causal or reverse causal relationship between depression/dysthymia and other types of female reproductive disorders. In summary, our study provides evidence for a causal relationship between genetically predicted depression and specific types of female reproductive disorders. Our findings emphasize the importance of depression management in the prevention and treatment of female reproductive disorders, notably including PCOS, ovarian cysts, AUB, and endometriosis.

A two-sample Mendelian randomization study explores metabolic profiling of different glycemic traits

We assessed the causal relation of four glycemic traits and type 2 diabetes liability with 167 metabolites using Mendelian randomization with various sensitivity analyses and a reverse Mendelian randomization analysis. We extracted instruments for fasting glucose, 2-h glucose, fasting insulin, and glycated hemoglobin from the Meta-Analyses of Glucose and Insulin-related traits Consortium (n = 200,622), and those for type 2 diabetes liability from a meta-analysis of multiple cohorts (148,726 cases, 965,732 controls) in Europeans. Outcome data were from summary statistics of 167 metabolites from the UK Biobank (n = 115,078). Fasting glucose and 2-h glucose were not associated with any metabolite. Higher glycated hemoglobin was associated with higher free cholesterol in small low-density lipoprotein. Type 2 diabetes liability and fasting insulin were inversely associated with apolipoprotein A1, total cholines, lipoprotein subfractions in high-density-lipoprotein and intermediate-density lipoproteins, and positively associated with aromatic amino acids. These findings indicate hyperglycemia-independent patterns and highlight the role of insulin in type 2 diabetes development. Further studies should evaluate these glycemic traits in type 2 diabetes diagnosis and clinical management.

The genome-wide association study of serum IgE levels demonstrated a shared genetic background in allergic diseases

Immunoglobulin E (IgE) synthessis is highly related to a variety of atopic diseases, and several genome-wide association studies (GWASs) have demonstrated the association between genes and IgE level. In this study, we conducted the largest genome-wide association study of IgE involving a Taiwanese Han population. Eight independent variants exhibited genome-wide significance. Among them, an intronic SNP of CD28, rs1181388, and an intergenic SNP, rs1002957030, on 11q23.2 were identified as novel signals for IgE. Seven of the loci were replicated successfully in a meta-analysis using data on Japanese population. Among all the human leukocyte antigen (HLA) regions, HLA-DQA1*03:02 - HLA-DQB1*03:03 was the most significant haplotype (OR = 1.25, SE = 0.02, FDR = 1.6 × 10-14), corresponding to HLA-DQA1 Asp160 and HLA-DQB1 Leu87 amino acid residues. The genetic correlation showed significance between IgE and allergic diseases including asthma, atopic dermatitis, and pollinosis. IgE PRS was significantly correlated with total IgE levels. Furthermore, the top decile IgE polygenic risk score (PRS) group had the highest risk of asthma for the Taiwan Biobank and Biobank Japan cohorts. IgE PRS may be used to aid in predicting the occurrence of allergic reactions before symptoms occur and biomarkers are detectable. Our study provided a more comprehensive understanding of the impact of genomic variants, including complex HLA alleles, on serum IgE levels.

Genome-wide analysis of spina bifida risk variants in a case-control study from Bangladesh

BACKGROUND

Human studies of genetic risk factors for neural tube defects, severe birth defects associated with long-term health consequences in surviving children, have predominantly been restricted to a subset of candidate genes in specific biological pathways including folate metabolism.

METHODS

In this study, we investigated the association of genetic variants spanning the genome with risk of spina bifida (i.e., myelomeningocele and meningocele) in a subset of families enrolled from December 2016 through December 2022 in a case-control study in Bangladesh, a population often underrepresented in genetic studies. Saliva DNA samples were analyzed using the Illumina Global Screening Array. We performed genetic association analyses to compare allele frequencies between 112 case and 121 control children, 272 mothers, and 128 trios.

RESULTS

In the transmission disequilibrium test analyses with trios only, we identified three novel exonic spina bifida risk loci, including rs140199800 (SULT1C2, p = 1.9 × 10-7), rs45580033 (ASB2, p = 4.2 × 10-10), and rs75426652 (LHPP, p = 7.2 × 10-14), after adjusting for multiple hypothesis testing. Association analyses comparing cases and controls, as well as models that included their mothers, did not identify genome-wide significant variants.

CONCLUSIONS

This study identified three novel single nucleotide polymorphisms involved in biological pathways not previously associated with neural tube defects. The study warrants replication in larger groups to validate findings and to inform targeted prevention strategies.

Causal effects between circulating immune cells and heart failure: evidence from a bidirectional Mendelian randomization study

BACKGROUND

Heart failure (HF) is a prevalent cardiac condition characterized by high mortality and morbidity rates. Immune cells play a pivotal role as crucial biomarkers in assessing the overall immune status of individuals. However, the causal relationship between circulating immune cells and the pathogenesis of HF remains an area requiring further investigation.

OBJECTIVES

The aim of this study was to investigate the genetic interactions between circulating immune cells and HF, and to further elucidate the genetic associations between different lymphocyte subsets and HF.

METHODS

We obtained genetic variants associated with circulating immune cells as instrumental variables (IVs) from the Blood Cell Consortium and publicly available HF summary data. We conducted additional subsets analyses on lymphocyte counts. Our study utilized two-sample and multivariate Mendelian randomization (MVMR) analysis to investigate the causal effect of immune cells on HF. The primary analysis employed inverse variance weighting (IVW) and was complemented by a series of sensitivity analyses.

RESULTS

The findings of the study showed that the IVW model demonstrated a significant correlation between an elevation in lymphocyte count and a decreased risk of HF (OR = 0.97, 95% CI, 0.94 - 1.00, P = 0.032). However, no such correlation was evident in the MVMR analysis for lymphocytes and HF. Furthermore, the examination of the lymphocyte subsets indicated that an increase in CD39+ CD4+ T-cell counts was notably linked to a reduced risk of HF (OR = 0.96, 95% CI, 0.95 - 0.98, P = 0.0002). The MVMR results confirmed that the association between CD39+ CD4+ T-cell counts and HF remained significant. There was no substantial evidence of reverse causality observed between circulating immune cells and HF.

CONCLUSION

Our MR research provided evidence for a causal relationship between lymphocyte cell and HF. Subsets analyses revealed a causal relationship between CD39+ CD4+ T lymphocytes and HF. These findings will facilitate a future understanding of the mechanisms underlying HF.

Causal association between the peripheral immunity and the risk and disease severity of multiple sclerosis

Background

Growing evidence links immunological responses to Multiple sclerosis (MS), but specific immune factors are still unclear.

Methods

Mendelian randomization (MR) was performed to investigate the association between peripheral hematological traits, MS risk, and its severity. Then, further subgroup analysis of immune counts and circulating cytokines and growth factors were performed.

Results

MR revealed higher white blood cell count (OR [95%CI] = 1.26 [1.10,1.44], P = 1.12E-03, P adjust = 3.35E-03) and lymphocyte count (OR [95%CI] = 1.31 [1.15,1.50], P = 5.37E-05, P adjust = 3.22E-04) increased the risk of MS. In further analysis, higher T cell absolute count (OR [95%CI] = 2.04 [1.36,3.08], P = 6.37E-04, P adjust = 2.19E-02) and CD4+ T cell absolute count (OR [95%CI] = 2.11 [1.37,3.24], P = 6.37E-04, P adjust = 2.19E-02), could increase MS risk. While increasing CD25++CD4+ T cell absolute count (OR [95%CI] = 0.75 [0.66,0.86], P = 2.12E-05, P adjust = 1.72E-03), CD25++CD4+ T cell in T cell (OR [95%CI] = 0.79[0.70,0.89], P = 8.54E-05, P adjust = 5.29E-03), CD25++CD4+ T cell in CD4+ T cell (OR [95%CI] = 0.80[0.72,0.89], P = 1.85E-05, P adjust = 1.72E-03), and CD25++CD8+ T cell in T cell (OR [95%CI] = 0.68[0.57,0.81], P = 2.22E-05, P adjust = 1.72E-03), were proved to be causally defensive for MS. For the disease severity, the suggestive association between some traits related to CD4+ T cell, Tregs and MS severity were demonstrated. Moreover, elevated levels of IL-2Ra had a detrimental effect on the risk of MS (OR [95%CI] = 1.22 [1.12,1.32], P = 3.20E-06, P adjust = 1.34E-04).

Conclusions

This study demonstrated a genetically predicted causal relationship between elevated peripheral immune cell counts and MS. Subgroup analysis revealed a specific contribution of peripheral immune cells, holding potential for further investigations into the underlying mechanisms of MS and its severity.

Improving polygenic risk prediction in admixed populations by explicitly modeling ancestral-differential effects via GAUDI

Polygenic risk scores (PRS) have shown successes in clinics, but most PRS methods focus only on participants with distinct primary continental ancestry without accommodating recently-admixed individuals with mosaic continental ancestry backgrounds for different segments of their genomes. Here, we develop GAUDI, a novel penalized-regression-based method specifically designed for admixed individuals. GAUDI explicitly models ancestry-differential effects while borrowing information across segments with shared ancestry in admixed genomes. We demonstrate marked advantages of GAUDI over other methods through comprehensive simulation and real data analyses for traits with associated variants exhibiting ancestral-differential effects. Leveraging data from the Women's Health Initiative study, we show that GAUDI improves PRS prediction of white blood cell count and C-reactive protein in African Americans by > 64% compared to alternative methods, and even outperforms PRS-CSx with large European GWAS for some scenarios. We believe GAUDI will be a valuable tool to mitigate disparities in PRS performance in admixed individuals.

Biological basis of extensive pleiotropy between blood traits and cancer risk

BACKGROUND

The immune system has a central role in preventing carcinogenesis. Alteration of systemic immune cell levels may increase cancer risk. However, the extent to which common genetic variation influences blood traits and cancer risk remains largely undetermined. Here, we identify pleiotropic variants and predict their underlying molecular and cellular alterations.

METHODS

Multivariate Cox regression was used to evaluate associations between blood traits and cancer diagnosis in cases in the UK Biobank. Shared genetic variants were identified from the summary statistics of the genome-wide association studies of 27 blood traits and 27 cancer types and subtypes, applying the conditional/conjunctional false-discovery rate approach. Analysis of genomic positions, expression quantitative trait loci, enhancers, regulatory marks, functionally defined gene sets, and bulk- and single-cell expression profiles predicted the biological impact of pleiotropic variants. Plasma small RNAs were sequenced to assess association with cancer diagnosis.

RESULTS

The study identified 4093 common genetic variants, involving 1248 gene loci, that contributed to blood-cancer pleiotropism. Genomic hotspots of pleiotropism include chromosomal regions 5p15-TERT and 6p21-HLA. Genes whose products are involved in regulating telomere length are found to be enriched in pleiotropic variants. Pleiotropic gene candidates are frequently linked to transcriptional programs that regulate hematopoiesis and define progenitor cell states of immune system development. Perturbation of the myeloid lineage is indicated by pleiotropic associations with defined master regulators and cell alterations. Eosinophil count is inversely associated with cancer risk. A high frequency of pleiotropic associations is also centered on the regulation of small noncoding Y-RNAs. Predicted pleiotropic Y-RNAs show specific regulatory marks and are overabundant in the normal tissue and blood of cancer patients. Analysis of plasma small RNAs in women who developed breast cancer indicates there is an overabundance of Y-RNA preceding neoplasm diagnosis.

CONCLUSIONS

This study reveals extensive pleiotropism between blood traits and cancer risk. Pleiotropism is linked to factors and processes involved in hematopoietic development and immune system function, including components of the major histocompatibility complexes, and regulators of telomere length and myeloid lineage. Deregulation of Y-RNAs is also associated with pleiotropism. Overexpression of these elements might indicate increased cancer risk.

Inherited polygenic effects on common hematological traits influence clonal selection on JAK2V617F and the development of myeloproliferative neoplasms

Myeloproliferative neoplasms (MPNs) are chronic cancers characterized by overproduction of mature blood cells. Their causative somatic mutations, for example, JAK2V617F, are common in the population, yet only a minority of carriers develop MPN. Here we show that the inherited polygenic loci that underlie common hematological traits influence JAK2V617F clonal expansion. We identify polygenic risk scores (PGSs) for monocyte count and plateletcrit as new risk factors for JAK2V617F positivity. PGSs for several hematological traits influenced the risk of different MPN subtypes, with low PGSs for two platelet traits also showing protective effects in JAK2V617F carriers, making them two to three times less likely to have essential thrombocythemia than carriers with high PGSs. We observed that extreme hematological PGSs may contribute to an MPN diagnosis in the absence of somatic driver mutations. Our study showcases how polygenic backgrounds underlying common hematological traits influence both clonal selection on somatic mutations and the subsequent phenotype of cancer.

Genetic relationships between high blood eosinophil count, asthma susceptibility, and asthma severity

OBJECTIVE

Genetic relationships between blood eosinophil count (BEC), asthma susceptibility, and severity are unclear. We sought to identify the genetic difference between type 2 (T2) and nontype 2 (non-T2) asthma (defined by BEC) and investigate genetic relationships between high BEC, asthma susceptibility, and severity.

METHODS

Genome-wide association studies (GWASs) were performed for T2 (n = 9,064; BEC ≥ 300 cells/μL) versus non-T2 asthma (n = 14,379; BEC < 150 cells/μL) and asthma susceptibility (37,227 asthmatics vs. 124,132 nonasthma controls) in the UK Biobank and asthma severity (moderate-to-severe asthma [n = 2,153] vs. mild asthma [n = 5165]) in the All of Us Research Program (AoURP). Genetic causality between BEC, asthma susceptibility, and severity were dissected using Mendelian randomization (MR).

RESULTS

High BEC was associated with asthma and decreased pulmonary function. GWASs revealed four sets of genetic variants (p < 5 × 10-8): genes associated with only BEC or asthma and genes associated with high BEC and asthma in the same or opposite direction. The C allele of rs653178 in ATXN2 was associated with high BEC, risk for autoimmune diseases, and protection for asthma. Genetic variants associated with BEC or asthma were not associated with asthma severity. MR indicated high BEC and asthma were in bidirectional causal relationship (p < .001); however, they were not causal for asthma severity.

CONCLUSIONS

Genetic variants associated with asthma or BEC and asthma severity are distinctive. High BEC is a risk factor for asthma; however, it is neither necessary nor sufficient for asthma susceptibility and severity.