Trans-ethnic and Ancestry-Specific Blood-Cell Genetics in 746,667 Individuals from 5 Global Populations

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.

The GWAS SNP rs80207740 modulates erythrocyte traits via allele-specific binding of IKZF1 and targeting XPO7 gene

Genome-wide association studies have identified many single nucleotide polymorphisms (SNPs) associated with erythrocyte traits. However, the functional variants and their working mechanisms remain largely unknown. Here, we reported that the SNP of rs80207740, which was associated with red blood cell (RBC) volume and hemoglobin content across populations, conferred enhancer activity to XPO7 gene via allele-differentially binding to Ikaros family zinc finger 1 (IKZF1). We showed that the region around rs80207740 was an erythroid-specific enhancer using reporter assays, and that the G-allele further enhanced activity. 3D genome evidence showed that the enhancer interacted with the XPO7 promoter, and eQTL analysis suggested that the G-allele upregulated expression of XPO7. We further showed that the rs80207740-G allele facilitated the binding of transcription factor IKZF1 in EMSA and ChIP analyses. Knockdown of IKZF1 and GATA1 resulted in decreased expression of Xpo7 in both human and mouse erythroid cells. Finally, we constructed Xpo7 knockout mouse by CRISPR/Cas9 and observed anemic phenotype with reduced volume and hemoglobin content of RBC, consistent to the effect of rs80207740 on erythrocyte traits. Overall, our study demonstrated that rs80207740 modulated erythroid indices by regulating IKZF1 binding and Xpo7 expression.

Cross-ancestry genome-wide association studies of brain imaging phenotypes

Genome-wide association studies of brain imaging phenotypes are mainly performed in European populations, but other populations are severely under-represented. Here, we conducted Chinese-alone and cross-ancestry genome-wide association studies of 3,414 brain imaging phenotypes in 7,058 Chinese Han and 33,224 white British participants. We identified 38 new associations in Chinese-alone analyses and 486 additional new associations in cross-ancestry meta-analyses at P < 1.46 × 10-11 for discovery and P < 0.05 for replication. We pooled significant autosomal associations identified by single- or cross-ancestry analyses into 6,443 independent associations, which showed uneven distribution in the genome and the phenotype subgroups. We further divided them into 44 associations with different effect sizes and 3,557 associations with similar effect sizes between ancestries. Loci of these associations were shared with 15 brain-related non-imaging traits including cognition and neuropsychiatric disorders. Our results provide a valuable catalog of genetic associations for brain imaging phenotypes in more diverse populations.

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.

Inferring causal cell types of human diseases and risk variants from candidate regulatory elements

The heritability of human diseases is extremely enriched in candidate regulatory elements (cRE) from disease-relevant cell types. Critical next steps are to infer which and how many cell types are truly causal for a disease (after accounting for co-regulation across cell types), and to understand how individual variants impact disease risk through single or multiple causal cell types. Here, we propose CT-FM and CT-FM-SNP, two methods that leverage cell-type-specific cREs to fine-map causal cell types for a trait and for its candidate causal variants, respectively. We applied CT-FM to 63 GWAS summary statistics (average N = 417K) using nearly one thousand cRE annotations, primarily coming from ENCODE4. CT-FM inferred 81 causal cell types with corresponding SNP-annotations explaining a high fraction of trait SNP-heritability (∼2/3 of the SNP-heritability explained by existing cREs), identified 16 traits with multiple causal cell types, highlighted cell-disease relationships consistent with known biology, and uncovered previously unexplored cellular mechanisms in psychiatric and immune-related diseases. Finally, we applied CT-FM-SNP to 39 UK Biobank traits and predicted high confidence causal cell types for 2,798 candidate causal non-coding SNPs. Our results suggest that most SNPs impact a phenotype through a single cell type, and that pleiotropic SNPs target different cell types depending on the phenotype context. Altogether, CT-FM and CT-FM-SNP shed light on how genetic variants act collectively and individually at the cellular level to impact disease risk.

Genotype-based "virtual" metabolomics in a clinical biobank identifies novel metabolite-disease associations

Introduction: Circulating metabolites act as biomarkers of dysregulated metabolism and may inform disease pathophysiology. A portion of the inter-individual variability in circulating metabolites is influenced by common genetic variation. We evaluated whether a genetics-based "virtual" metabolomics approach can identify novel metabolite-disease associations. Methods: We examined the association between polygenic scores for 724 metabolites with 1,247 clinical phenotypes in the BioVU DNA biobank, comprising 57,735 European ancestry and 15,754 African ancestry participants. We applied Mendelian randomization (MR) to probe significant relationships and validated significant MR associations using independent GWAS of candidate phenotypes. Results and Discussion: We found significant associations between 336 metabolites and 168 phenotypes in European ancestry and 107 metabolites and 56 phenotypes in African ancestry. Of these metabolite-disease pairs, MR analyses confirmed associations between 73 metabolites and 53 phenotypes in European ancestry. Of 22 metabolitephenotype pairs evaluated for replication in independent GWAS, 16 were significant (false discovery rate p < 0.05). These included associations between bilirubin and X-21796 with cholelithiasis, phosphatidylcholine (16:0/22:5n3,18:1/20:4) and arachidonate with inflammatory bowel disease and Crohn's disease, and campesterol with coronary artery disease and myocardial infarction. These associations may represent biomarkers or potentially targetable mediators of disease risk.

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.

Causal association of blood cell traits with inflammatory bowel diseases: a Mendelian randomization study

Background

Observational studies have found associations between blood cell traits and inflammatory bowel diseases (IBDs), whereas the causality and dose-effect relationships are still undetermined.

Methods

Two-sample Mendelian randomization (MR) analyses using linear regression approaches, as well as Bayesian model averaging (MR-BMA), were conducted to identify and prioritize the causal blood cell traits for Crohn's disease (CD) and ulcerative colitis (UC). An observational study was also performed using restricted cubic spline (RCS) to explore the relationship between important blood cell traits and IBDs.

Results

Our uvMR analysis using the random effects inverse variance weighted (IVW) method identified eosinophil (EOS) as a causal factor for UC (OR = 1.36; 95% CI: 1.13, 1.63). Our MR-BMA analysis further prioritized that high level of lymphocyte (LYM) decreased CD risk (MIP = 0.307; θ ^ MACE = -0.059; PP = 0.189; θ ^ λ = -0.173), whereas high level of EOS increased UC risk (MIP = 0.824; θ ^ MACE = 0.198; PP = 0.627; θ ^ λ = 0.239). Furthermore, the observational study clearly depicts the nonlinear relationship between important blood cell traits and the risk of IBDs.

Conclusion

Using MR approaches, several blood cell traits were identified as risk factors of CD and UC, which could be used as potential targets for the management of IBDs. Stratified genome-wide association studies (GWASs) based on the concentration of traits would be helpful owing to the nonlinear relationships between blood cell traits and IBDs, as demonstrated in our clinical observational study. Together, these findings could shed light on the clinical strategies applied to the management of CD and UC.

Analysis of blood methylation quantitative trait loci in East Asians reveals ancestry-specific impacts on complex traits

Methylation quantitative trait loci (mQTLs) are essential for understanding the role of DNA methylation changes in genetic predisposition, yet they have not been fully characterized in East Asians (EAs). Here we identified mQTLs in whole blood from 3,523 Chinese individuals and replicated them in additional 1,858 Chinese individuals from two cohorts. Over 9% of mQTLs displayed specificity to EAs, facilitating the fine-mapping of EA-specific genetic associations, as shown for variants associated with height. Trans-mQTL hotspots revealed biological pathways contributing to EA-specific genetic associations, including an ERG-mediated 233 trans-mCpG network, implicated in hematopoietic cell differentiation, which likely reflects binding efficiency modulation of the ERG protein complex. More than 90% of mQTLs were shared between different blood cell lineages, with a smaller fraction of lineage-specific mQTLs displaying preferential hypomethylation in the respective lineages. Our study provides new insights into the mQTL landscape across genetic ancestries and their downstream effects on cellular processes and diseases/traits.

Leveraging functional genomic annotations and genome coverage to improve polygenic prediction of complex traits within and between ancestries

We develop a method, SBayesRC, that integrates genome-wide association study (GWAS) summary statistics with functional genomic annotations to improve polygenic prediction of complex traits. Our method is scalable to whole-genome variant analysis and refines signals from functional annotations by allowing them to affect both causal variant probability and causal effect distribution. We analyze 50 complex traits and diseases using ∼7 million common single-nucleotide polymorphisms (SNPs) and 96 annotations. SBayesRC improves prediction accuracy by 14% in European ancestry and up to 34% in cross-ancestry prediction compared to the baseline method SBayesR, which does not use annotations, and outperforms other methods, including LDpred2, LDpred-funct, MegaPRS, PolyPred-S and PRS-CSx. Investigation of factors affecting prediction accuracy identifies a significant interaction between SNP density and annotation information, suggesting whole-genome sequence variants with annotations may further improve prediction. Functional partitioning analysis highlights a major contribution of evolutionary constrained regions to prediction accuracy and the largest per-SNP contribution from nonsynonymous SNPs.

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.

Unraveling COVID-19 relationship with anxiety disorders and symptoms using genome-wide data

BACKGROUND

There is still a limited understanding of the dynamics contributing to the comorbidity of COVID-19 and anxiety outcomes.

METHODS

To dissect the pleiotropic mechanisms contributing to COVID-19/anxiety comorbidity, we used genome-wide data from UK Biobank (up to 420,531 participants), FinnGen Project (up to 329,077 participants), Million Veteran Program (175,163 participants), and COVID-19 Host Genetics Initiative (up to 122,616 cases and 2,475,240 controls). Specifically, we assessed global and local genetic correlation and genetically inferred effects linking COVID-19 outcomes (infection, hospitalization, and severe respiratory symptoms) to anxiety disorders and symptoms.

RESULTS

We observed a strong genetic correlation of anxiety disorder with COVID-19 positive status (rg = 0.35, p = 2×10-4) and COVID-19 hospitalization (rg = 0.31, p = 7.2×10-4). Among anxiety symptoms, "Tense, sore, or aching muscles during worst period of anxiety" was genetically correlated with COVID-19 positive status (rg = 0.33, p = 0.001), while "Frequent trouble falling or staying asleep during worst period of anxiety" was genetically correlated with COVID-19 hospitalization (rg = 0.24, p = 0.004). Through a latent causal variable analysis, we observed that COVID-19 outcomes have statistically significant genetic causality proportion (gcp) on anxiety symptoms (e.g., COVID-19 positive status→"Recent easy annoyance or irritability" │gcp│ = 0.18, p = 6.72×10-17). Conversely, anxiety disorders appear to have a possible causal effect on COVID-19 (│gcp│ = 0.38, p = 3.17×10-9). Additionally, we also identified multiple loci with evidence of local genetic correlation between anxiety and COVID-19. These appear to be related to genetic effects shared with lung function, brain morphology, alcohol and tobacco use, and hematologic parameters.

CONCLUSIONS

This study provided insights into the pleiotropic mechanisms linking COVID-19 and anxiety outcomes, suggesting differences between dynamics related to anxiety disorders and those related to anxiety symptoms.

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.

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

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

Genomic and phenotypic correlates of mosaic loss of chromosome Y in blood

Mosaic loss of Y (mLOY) is the most common somatic chromosomal alteration detected in human blood. The presence of mLOY is associated with altered blood cell counts and increased risk of Alzheimer's disease, solid tumors, and other age-related diseases. We sought to gain a better understanding of genetic drivers and associated phenotypes of mLOY through analyses of whole genome sequencing of a large set of genetically diverse males from the Trans-Omics for Precision Medicine (TOPMed) program. This approach enabled us to identify differences in mLOY frequencies across populations defined by genetic similarity, revealing a higher frequency of mLOY in the European American (EA) ancestry group compared to those of Hispanic American (HA), African American (AA), and East Asian (EAS) ancestry. Further, we identified two genes ( CFHR1 and LRP6 ) that harbor multiple rare, putatively deleterious variants associated with mLOY susceptibility, show that subsets of human hematopoietic stem cells are enriched for activity of mLOY susceptibility variants, and that certain alleles on chromosome Y are more likely to be lost than others.

Multi-trait GWAS for diverse ancestries: mapping the knowledge gap

BACKGROUND

Approximately 95% of samples analyzed in univariate genome-wide association studies (GWAS) are of European ancestry. This bias toward European ancestry populations in association screening also exists for other analyses and methods that are often developed and tested on European ancestry only. However, existing data in non-European populations, which are often of modest sample size, could benefit from innovative approaches as recently illustrated in the context of polygenic risk scores.

METHODS

Here, we extend and assess the potential limitations and gains of our multi-trait GWAS pipeline, JASS (Joint Analysis of Summary Statistics), for the analysis of non-European ancestries. To this end, we conducted the joint GWAS of 19 hematological traits and glycemic traits across five ancestries (European (EUR), admixed American (AMR), African (AFR), East Asian (EAS), and South-East Asian (SAS)).

RESULTS

We detected 367 new genome-wide significant associations in non-European populations (15 in Admixed American (AMR), 72 in African (AFR) and 280 in East Asian (EAS)). New associations detected represent 5%, 17% and 13% of associations in the AFR, AMR and EAS populations, respectively. Overall, multi-trait testing increases the replication of European associated loci in non-European ancestry by 15%. Pleiotropic effects were highly similar at significant loci across ancestries (e.g. the mean correlation between multi-trait genetic effects of EUR and EAS ancestries was 0.88). For hematological traits, strong discrepancies in multi-trait genetic effects are tied to known evolutionary divergences: the ARKC1 loci, which is adaptive to overcome p.vivax induced malaria.

CONCLUSIONS

Multi-trait GWAS can be a valuable tool to narrow the genetic knowledge gap between European and non-European populations.

Genetic evidence for involvement of β2-adrenergic receptor in brown adipose tissue thermogenesis in humans

BACKGROUND

Sympathetic activation of brown adipose tissue (BAT) thermogenesis can ameliorate obesity and related metabolic abnormalities. However, crucial subtypes of the β-adrenergic receptor (AR), as well as effects of its genetic variants on functions of BAT, remains unclear in humans. We conducted association analyses of genes encoding β-ARs and BAT activity in human adults.

METHODS

Single nucleotide polymorphisms (SNPs) in β1-, β2-, and β3-AR genes (ADRB1, ADRB2, and ADRB3) were tested for the association with BAT activity under mild cold exposure (19 °C, 2 h) in 399 healthy Japanese adults. BAT activity was measured using fluorodeoxyglucose-positron emission tomography and computed tomography (FDG-PET/CT). To validate the results, we assessed the effects of SNPs in the two independent populations comprising 277 healthy East Asian adults using near-infrared time-resolved spectroscopy (NIRTRS) or infrared thermography (IRT). Effects of SNPs on physiological responses to intensive cold exposure were tested in 42 healthy Japanese adult males using an artificial climate chamber.

RESULTS

We found a significant association between a functional SNP (rs1042718) in ADRB2 and BAT activity assessed with FDG-PET/CT (p < 0.001). This SNP also showed an association with cold-induced thermogenesis in the population subset. Furthermore, the association was replicated in the two other independent populations; BAT activity was evaluated by NIRTRS or IRT (p < 0.05). This SNP did not show associations with oxygen consumption and cold-induced thermogenesis under intensive cold exposure, suggesting the irrelevance of shivering thermogenesis. The SNPs of ADRB1 and ADRB3 were not associated with these BAT-related traits.

CONCLUSIONS

The present study supports the importance of β2-AR in the sympathetic regulation of BAT thermogenesis in humans. The present collection of DNA samples is the largest to which information on the donor's BAT activity has been assigned and can serve as a reference for further in-depth understanding of human BAT function.

Improved multi-ancestry fine-mapping identifies cis-regulatory variants underlying molecular traits and disease risk

Multi-ancestry statistical fine-mapping of cis-molecular quantitative trait loci (cis-molQTL) aims to improve the precision of distinguishing causal cis-molQTLs from tagging variants. However, existing approaches fail to reflect shared genetic architectures. To solve this limitation, we present the Sum of Shared Single Effects (SuShiE) model, which leverages LD heterogeneity to improve fine-mapping precision, infer cross-ancestry effect size correlations, and estimate ancestry-specific expression prediction weights. We apply SuShiE to mRNA expression measured in PBMCs (n=956) and LCLs (n=814) together with plasma protein levels (n=854) from individuals of diverse ancestries in the TOPMed MESA and GENOA studies. We find SuShiE fine-maps cis-molQTLs for 16% more genes compared with baselines while prioritizing fewer variants with greater functional enrichment. SuShiE infers highly consistent cis-molQTL architectures across ancestries on average; however, we also find evidence of heterogeneity at genes with predicted loss-of-function intolerance, suggesting that environmental interactions may partially explain differences in cis-molQTL effect sizes across ancestries. Lastly, we leverage estimated cis-molQTL effect-sizes to perform individual-level TWAS and PWAS on six white blood cell-related traits in AOU Biobank individuals (n=86k), and identify 44 more genes compared with baselines, further highlighting its benefits in identifying genes relevant for complex disease risk. Overall, SuShiE provides new insights into the cis-genetic architecture of molecular traits.

DNA binding analysis of rare variants in homeodomains reveals homeodomain specificity-determining residues

Homeodomains (HDs) are the second largest class of DNA binding domains (DBDs) among eukaryotic sequence-specific transcription factors (TFs) and are the TF structural class with the largest number of disease-associated mutations in the Human Gene Mutation Database (HGMD). Despite numerous structural studies and large-scale analyses of HD DNA binding specificity, HD-DNA recognition is still not fully understood. Here, we analyze 92 human HD mutants, including disease-associated variants and variants of uncertain significance (VUS), for their effects on DNA binding activity. Many of the variants alter DNA binding affinity and/or specificity. Detailed biochemical analysis and structural modeling identifies 14 previously unknown specificity-determining positions, 5 of which do not contact DNA. The same missense substitution at analogous positions within different HDs often exhibits different effects on DNA binding activity. Variant effect prediction tools perform moderately well in distinguishing variants with altered DNA binding affinity, but poorly in identifying those with altered binding specificity. Our results highlight the need for biochemical assays of TF coding variants and prioritize dozens of variants for further investigations into their pathogenicity and the development of clinical diagnostics and precision therapies.

A genome-wide association study provides insights into the genetic etiology of 57 essential and non-essential trace elements in humans

Trace elements are important for human health but may exert toxic or adverse effects. Mechanisms of uptake, distribution, metabolism, and excretion are partly under genetic control but have not yet been extensively mapped. Here we report a comprehensive multi-element genome-wide association study of 57 essential and non-essential trace elements. We perform genome-wide association meta-analyses of 14 trace elements in up to 6564 Scandinavian whole blood samples, and genome-wide association studies of 43 trace elements in up to 2819 samples measured only in the Trøndelag Health Study (HUNT). We identify 11 novel genetic loci associated with blood concentrations of arsenic, cadmium, manganese, selenium, and zinc in genome-wide association meta-analyses. In HUNT, several genome-wide significant loci are also indicated for other trace elements. Using two-sample Mendelian randomization, we find several indications of weak to moderate effects on health outcomes, the most precise being a weak harmful effect of increased zinc on prostate cancer. However, independent validation is needed. Our current understanding of trace element-associated genetic variants may help establish consequences of trace elements on human health.

Effects of genetically predicted posttraumatic stress disorder on autoimmune phenotypes

Observational studies suggest that posttraumatic stress disorder (PTSD) increases risk for various autoimmune diseases. Insights into shared biology and causal relationships between these diseases may inform intervention approaches to PTSD and co-morbid autoimmune conditions. We investigated the shared genetic contributions and causal relationships between PTSD, 18 autoimmune diseases, and 3 immune/inflammatory biomarkers. Univariate MiXeR was used to contrast the genetic architectures of phenotypes. Genetic correlations were estimated using linkage disequilibrium score regression. Bi-directional, two-sample Mendelian randomization (MR) was performed using independent, genome-wide significant single nucleotide polymorphisms; inverse variance weighted and weighted median MR estimates were evaluated. Sensitivity analyses for uncorrelated (MR PRESSO) and correlated horizontal pleiotropy (CAUSE) were also performed. PTSD was considerably more polygenic (10,863 influential variants) than autoimmune diseases (median 255 influential variants). However, PTSD evidenced significant genetic correlation with nine autoimmune diseases and three inflammatory biomarkers. PTSD had putative causal effects on autoimmune thyroid disease (p = 0.00009) and C-reactive protein (CRP) (p = 4.3 × 10-7). Inferences were not substantially altered by sensitivity analyses. Additionally, the PTSD-autoimmune thyroid disease association remained significant in multivariable MR analysis adjusted for genetically predicted inflammatory biomarkers as potential mechanistic pathway variables. No autoimmune disease had a significant causal effect on PTSD (all p values > 0.05). Although causal effect models were supported for associations of PTSD with CRP, shared pleiotropy was adequate to explain a putative causal effect of CRP on PTSD (p = 0.18). In summary, our results suggest a significant genetic overlap between PTSD, autoimmune diseases, and biomarkers of inflammation. PTSD has a putative causal effect on autoimmune thyroid disease, consistent with existing epidemiologic evidence. A previously reported causal effect of CRP on PTSD is potentially confounded by shared genetics. Together, results highlight the nuanced links between PTSD, autoimmune disorders, and associated inflammatory signatures, and suggest the importance of targeting related pathways to protect against disease and disability.

Tissue-specific enhancer-gene maps from multimodal single-cell data identify causal disease alleles

Translating genome-wide association study (GWAS) loci into causal variants and genes requires accurate cell-type-specific enhancer-gene maps from disease-relevant tissues. Building enhancer-gene maps is essential but challenging with current experimental methods in primary human tissues. Here we developed a nonparametric statistical method, SCENT (single-cell enhancer target gene mapping), that models association between enhancer chromatin accessibility and gene expression in single-cell or nucleus multimodal RNA sequencing and ATAC sequencing data. We applied SCENT to 9 multimodal datasets including >120,000 single cells or nuclei and created 23 cell-type-specific enhancer-gene maps. These maps were highly enriched for causal variants in expression quantitative loci and GWAS for 1,143 diseases and traits. We identified likely causal genes for both common and rare diseases and linked somatic mutation hotspots to target genes. We demonstrate that application of SCENT to multimodal data from disease-relevant human tissue enables the scalable construction of accurate cell-type-specific enhancer-gene maps, essential for defining noncoding variant function.

Genetic associations of protein-coding variants in venous thromboembolism

Previous genetic studies of venous thromboembolism (VTE) have been largely limited to common variants, leaving the genetic determinants relatively incomplete. We performed an exome-wide association study of VTE among 14,723 cases and 334,315 controls. Fourteen known and four novel genes (SRSF6, PHPT1, CGN, and MAP3K2) were identified through protein-coding variants, with broad replication in the FinnGen cohort. Most genes we discovered exhibited the potential to predict future VTE events in longitudinal analysis. Notably, we provide evidence for the additive contribution of rare coding variants to known genome-wide polygenic risk in shaping VTE risk. The identified genes were enriched in pathways affecting coagulation and platelet activation, along with liver-specific expression. The pleiotropic effects of these genes indicated the potential involvement of coagulation factors, blood cell traits, liver function, and immunometabolic processes in VTE pathogenesis. In conclusion, our study unveils the valuable contribution of protein-coding variants in VTE etiology and sheds new light on its risk stratification.

Haemoglobin thresholds to define anaemia from age 6 months to 65 years: estimates from international data sources

BACKGROUND

Detection of anaemia is crucial for clinical medicine and public health. Current WHO anaemia definitions are based on statistical thresholds (fifth centiles) set more than 50 years ago. We sought to establish evidence for the statistical haemoglobin thresholds for anaemia that can be applied globally and inform WHO and clinical guidelines.

METHODS

In this analysis we identified international data sources from populations in the USA, England, Australia, China, the Netherlands, Canada, Ecuador, and Bangladesh with sufficient clinical and laboratory information collected between 1998 and 2020 to obtain a healthy reference sample. Individuals with clinical or biochemical evidence of a condition that could reduce haemoglobin concentrations were excluded. We estimated haemoglobin thresholds (ie, 5th centiles) for children aged 6-23 months, 24-59 months, 5-11 years, and 12-17 years, and adults aged 18-65 years (including during pregnancy) for individual datasets and pooled across data sources. We also collated findings from three large-scale genetic studies to summarise genetic variants affecting haemoglobin concentrations in different ancestral populations.

FINDINGS

We identified eight data sources comprising 18 individual datasets that were eligible for inclusion in the analysis. In pooled analyses, the haemoglobin fifth centile was 104·4 g/L (90% CI 103·5-105·3) in 924 children aged 6-23 months, 110·2 g/L (109·5-110·9) in 1874 children aged 24-59 months, and 114·4 g/L (113·6-115·2) in 1839 children aged 5-11 years. Values diverged by sex in adolescents and adults. In pooled analyses, the fifth centile was 122·2 g/L (90% CI 121·3-123·1) in 1741 female adolescents aged 12-17 years and 128·2 g/L (126·4-130·0) in 1103 male adolescents aged 12-17 years. In pooled analyses of adults aged 18-65 years, the fifth centile was 119·7 g/L (90% CI 119·1-120·3) in 3640 non-pregnant females and 134·9 g/L (134·2-135·6) in 2377 males. Fifth centiles in pregnancy were 110·3 g/L (90% CI 109·5-111·0) in the first trimester (n=772) and 105·9 g/L (104·0-107·7) in the second trimester (n=111), with insufficient data for analysis in the third trimester. There were insufficient data for adults older than 65 years. We did not identify ancestry-specific high prevalence of non-clinically relevant genetic variants that influence haemoglobin concentrations.

INTERPRETATION

Our results enable global harmonisation of clinical and public health haemoglobin thresholds for diagnosis of anaemia. Haemoglobin thresholds are similar between sexes until adolescence, after which males have higher thresholds than females. We did not find any evidence that thresholds should differ between people of differering ancestries.

FUNDING

World Health Organization and the Bill & Melinda Gates Foundation.

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.

Causal associations of genetically predicted gut microbiota and blood metabolites with inflammatory states and risk of infections: a Mendelian randomization analysis

Background

Inflammation serves as a key pathologic mediator in the progression of infections and various diseases, involving significant alterations in the gut microbiome and metabolism. This study aims to probe into the potential causal relationships between gut microbial taxa and human blood metabolites with various serum inflammatory markers (CRP, SAA1, IL-6, TNF-α, WBC, and GlycA) and the risks of seven common infections (gastrointestinal infections, dysentery, pneumonia, bacterial pneumonia, bronchopneumonia and lung abscess, pneumococcal pneumonia, and urinary tract infections).

Methods

Two-sample Mendelian randomization (MR) analysis was performed using inverse variance weighted (IVW), maximum likelihood, MR-Egger, weighted median, and MR-PRESSO.

Results

After adding other MR models and sensitivity analyses, genus Roseburia was simultaneously associated adversely with CRP (Beta IVW = -0.040) and SAA1 (Beta IVW = -0.280), and family Bifidobacteriaceae was negatively associated with both CRP (Beta IVW = -0.034) and pneumonia risk (Beta IVW = -0.391). After correction by FDR, only glutaroyl carnitine remained significantly associated with elevated CRP levels (Beta IVW = 0.112). Additionally, threonine (Beta IVW = 0.200) and 1-heptadecanoylglycerophosphocholine (Beta IVW = -0.246) were found to be significantly associated with WBC levels. Three metabolites showed similar causal effects on different inflammatory markers or infectious phenotypes, stearidonate (18:4n3) was negatively related to SAA1 and urinary tract infections, and 5-oxoproline contributed to elevated IL-6 and SAA1 levels. In addition, 7-methylguanine showed a positive correlation with dysentery and bacterial pneumonia.

Conclusion

This study provides novel evidence confirming the causal effects of the gut microbiome and the plasma metabolite profile on inflammation and the risk of infection. These potential molecular alterations may aid in the development of new targets for the intervention and management of disorders associated with inflammation and infections.

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.

SCGB1D2 inhibits growth of Borrelia burgdorferi and affects susceptibility to Lyme disease

Lyme disease is a tick-borne disease caused by bacteria of the genus Borrelia. The host factors that modulate susceptibility for Lyme disease have remained mostly unknown. Using epidemiological and genetic data from FinnGen and Estonian Biobank, we identify two previously known variants and an unknown common missense variant at the gene encoding for Secretoglobin family 1D member 2 (SCGB1D2) protein that increases the susceptibility for Lyme disease. Using live Borrelia burgdorferi (Bb) we find that recombinant reference SCGB1D2 protein inhibits the growth of Bb in vitro more efficiently than the recombinant protein with SCGB1D2 P53L deleterious missense variant. Finally, using an in vivo murine infection model we show that recombinant SCGB1D2 prevents infection by Borrelia in vivo. Together, these data suggest that SCGB1D2 is a host defense factor present in the skin, sweat, and other secretions which protects against Bb infection and opens an exciting therapeutic avenue for Lyme disease.

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.

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.

Genetic Analysis of Perceived Youthfulness Reveals Differences in How Men's and Women's Age Is Assessed

Skin aging is a natural process that occurs over time but can be accelerated by sun exposure. Measuring skin age in a large population can provide insight into the extent of skin damage from sun exposure and skin cancer risk. Understanding the genetics of skin aging, within and across sexes (males and females), could improve our understanding of the genetic drivers of both skin aging and skin cancer. We used UK Biobank data to examine the genetic overlap between perceived youthfulness and traits relevant to actinic photoaging. Our GWAS identified 22 genome-wide significant loci for women and 43 for men. The genetic correlation (rg) between perceived youthfulness in men and women was significantly less than unity (rg = 0.75, 95% confidence interval = 0.69-0.80), suggesting a gene-by-sex interaction. In women, perceived youthfulness was modestly correlated with keratinocyte cancer (rg = -0.19) and skin tanning (rg = 0.18). In men, perceived youthfulness was correlated with male-pattern baldness (rg = -0.23). This suggests that the genetic architecture of perceived youthfulness may differ between sexes, with genes influencing skin tanning and skin cancer susceptibility driving the difference in women, whereas genes influencing male-pattern baldness and other puberty-related traits drive the difference in men. We recommend that future genetic analysis of skin aging include a sex-stratified component.

Downregulation of the kidney glucagon receptor, essential for renal function and systemic homeostasis, contributes to chronic kidney disease

The glucagon receptor (GCGR) in the kidney is expressed in nephron tubules. In humans and animal models with chronic kidney disease, renal GCGR expression is reduced. However, the role of kidney GCGR in normal renal function and in disease development has not been addressed. Here, we examined its role by analyzing mice with constitutive or conditional kidney-specific loss of the Gcgr. Adult renal Gcgr knockout mice exhibit metabolic dysregulation and a functional impairment of the kidneys. These mice exhibit hyperaminoacidemia associated with reduced kidney glucose output, oxidative stress, enhanced inflammasome activity, and excess lipid accumulation in the kidney. Upon a lipid challenge, they display maladaptive responses with acute hypertriglyceridemia and chronic proinflammatory and profibrotic activation. In aged mice, kidney Gcgr ablation elicits widespread renal deposition of collagen and fibronectin, indicative of fibrosis. Taken together, our findings demonstrate an essential role of the renal GCGR in normal kidney metabolic and homeostatic functions. Importantly, mice deficient for kidney Gcgr recapitulate some of the key pathophysiological features of chronic kidney disease.

Genomic data in the All of Us Research Program

Comprehensively mapping the genetic basis of human disease across diverse individuals is a long-standing goal for the field of human genetics1-4. The All of Us Research Program is a longitudinal cohort study aiming to enrol a diverse group of at least one million individuals across the USA to accelerate biomedical research and improve human health5,6. Here we describe the programme's genomics data release of 245,388 clinical-grade genome sequences. This resource is unique in its diversity as 77% of participants are from communities that are historically under-represented in biomedical research and 46% are individuals from under-represented racial and ethnic minorities. All of Us identified more than 1 billion genetic variants, including more than 275 million previously unreported genetic variants, more than 3.9 million of which had coding consequences. Leveraging linkage between genomic data and the longitudinal electronic health record, we evaluated 3,724 genetic variants associated with 117 diseases and found high replication rates across both participants of European ancestry and participants of African ancestry. Summary-level data are publicly available, and individual-level data can be accessed by researchers through the All of Us Researcher Workbench using a unique data passport model with a median time from initial researcher registration to data access of 29 hours. We anticipate that this diverse dataset will advance the promise of genomic medicine for all.

Human Leukocyte Antigen Signatures as Pathophysiological Discriminants of Microscopic Colitis Subtypes

BACKGROUND AND AIMS

Microscopic colitis [MC] is currently regarded as an inflammatory bowel disease that manifests as two subtypes: collagenous colitis [CC] and lymphocytic colitis [LC]. Whether these represent a clinical continuum or distinct entities is, however, an open question. Genetic investigations may contribute important insight into their respective pathophysiologies.

METHODS

We conducted a genome-wide association study [GWAS] meta-analysis in 1498 CC, 373 LC patients, and 13 487 controls from Europe and the USA, combined with publicly available MC GWAS data from UK Biobank and FinnGen [2599 MC cases and 552 343 controls in total]. Human leukocyte antigen [HLA] alleles and polymorphic residues were imputed and tested for association, including conditional analyses for the identification of key causative variants and residues. Genetic correlations with other traits and diagnoses were also studied.

RESULTS

We detected strong HLA association with CC, and conditional analyses highlighted the DRB1*03:01 allele and its residues Y26, N77, and R74 as key to this association (best p = 1.4 × 10-23, odds ratio [OR] = 1.96). Nominally significant genetic correlations were detected between CC and pneumonia [rg = 0.77; p = 0.048] and oesophageal diseases [rg = 0.45, p = 0.023]. An additional locus was identified in MC GWAS analyses near the CLEC16A and RMI2 genes on chromosome 16 [rs35099084, p = 2.0 × 10-8, OR = 1.31]. No significant association was detected for LC.

CONCLUSION

Our results suggest CC and LC have distinct pathophysiological underpinnings, characterised by an HLA predisposing role only in CC. This challenges existing classifications, eventually calling for a re-evaluation of the utility of MC umbrella definitions.

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.

Germline predisposition for clonal hematopoiesis

Clonal hematopoiesis (CH) is an entity hallmarked by skewed hematopoiesis with persistent overrepresentation of cells from a common stem/progenitor lineage harboring single-nucleotide variants and/or insertions/deletions. CH is a common and age-related phenomenon that is associated with an increased risk of hematological malignancies, cardiovascular disease, and all-cause mortality. While CH is a term of the hematological aspect, there exists a complex interaction with other organ systems, especially the cardiovascular system. The strongest factor in the development of CH is aging, however, other multiple factors also affect the development of CH including lifestyle-related factors and co-morbid diseases. In recent years, germline genetic factors have been linked to CH risk. In this review, we synthesize what is currently known about how genetic variation affects the risk of CH, how this genetic architecture intersects with myeloid neoplasms, and future prospects for CH.

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.

The Impact of Patterns in Linkage Disequilibrium and Sequencing Quality on the Imprint of Balancing Selection

Regions under balancing selection are characterized by dense polymorphisms and multiple persistent haplotypes, along with other sequence complexities. Successful identification of these patterns depends on both the statistical approach and the quality of sequencing. To address this challenge, at first, a new statistical method called LD-ABF was developed, employing efficient Bayesian techniques to effectively test for balancing selection. LD-ABF demonstrated the most robust detection of selection in a variety of simulation scenarios, compared against a range of existing tests/tools (Tajima's D, HKA, Dng, BetaScan, and BalLerMix). Furthermore, the impact of the quality of sequencing on detection of balancing selection was explored, as well, using: (i) SNP genotyping and exome data, (ii) targeted high-resolution HLA genotyping (IHIW), and (iii) whole-genome long-read sequencing data (Pangenome). In the analysis of SNP genotyping and exome data, we identified known targets and 38 new selection signatures in genes not previously linked to balancing selection. To further investigate the impact of sequencing quality on detection of balancing selection, a detailed investigation of the MHC was performed with high-resolution HLA typing data. Higher quality sequencing revealed the HLA-DQ genes consistently demonstrated strong selection signatures otherwise not observed from the sparser SNP array and exome data. The HLA-DQ selection signature was also replicated in the Pangenome samples using considerably less samples but, with high-quality long-read sequence data. The improved statistical method, coupled with higher quality sequencing, leads to more consistent identification of selection and enhanced localization of variants under selection, particularly in complex regions.

Uncovering the architecture of selection in two Bos taurus cattle breeds

Directional selection alters the genome via hard sweeps, soft sweeps, and polygenic selection. However, mapping polygenic selection is difficult because it does not leave clear signatures on the genome like a selective sweep. In populations with temporally stratified genotypes, the Generation Proxy Selection Mapping (GPSM) method identifies variants associated with generation number (or appropriate proxy) and thus variants undergoing directional allele frequency changes. Here, we use GPSM on two large datasets of beef cattle to detect associations between an animal's generation and 11 million imputed SNPs. Using these datasets with high power and dense mapping resolution, GPSM detected a total of 294 unique loci actively under selection in two cattle breeds. We observed that GPSM has a high power to detect selection in the very recent past (<10 years), even when allele frequency changes are small. Variants identified by GPSM reside in genomic regions associated with known breed-specific selection objectives, such as fertility and maternal ability in Red Angus, and carcass merit and coat color in Simmental. Over 60% of the selected loci reside in or near (<50 kb) annotated genes. Using haplotype-based and composite selective sweep statistics, we identify hundreds of putative selective sweeps that likely occurred earlier in the evolution of these breeds; however, these sweeps have little overlap with recent polygenic selection. This makes GPSM a complementary approach to sweep detection methods when temporal genotype data are available. The selected loci that we identify across methods demonstrate the complex architecture of selection in domesticated cattle.

A genome-wide association study reveals the relationship between human genetic variation and the nasal microbiome

The nasal cavity harbors diverse microbiota that contributes to human health and respiratory diseases. However, whether and to what extent the host genome shapes the nasal microbiome remains largely unknown. Here, by dissecting the human genome and nasal metagenome data from 1401 healthy individuals, we demonstrated that the top three host genetic principal components strongly correlated with the nasal microbiota diversity and composition. The genetic association analyses identified 63 genome-wide significant loci affecting the nasal microbial taxa and functions, of which 2 loci reached study-wide significance (p < 1.7 × 10-10): rs73268759 within CAMK2A associated with genus Actinomyces and family Actinomycetaceae; and rs35211877 near POM121L12 with Gemella asaccharolytica. In addition to respiratory-related diseases, the associated loci are mainly implicated in cardiometabolic or neuropsychiatric diseases. Functional analysis showed the associated genes were most significantly expressed in the nasal airway epithelium tissue and enriched in the calcium signaling and hippo signaling pathway. Further observational correlation and Mendelian randomization analyses consistently suggested the causal effects of Serratia grimesii and Yokenella regensburgei on cardiometabolic biomarkers (cystine, glutamic acid, and creatine). This study suggested that the host genome plays an important role in shaping the nasal microbiome.

The genetic etiology of body fluids on chronic obstructive airways disease

BACKGROUND

Numerous studies have documented significant alterations in the bodily fluids of Chronic Obstructive Pulmonary Disease (COPD) patients. However, existing literature lacks causal inference due to residual confounding and reverse causality.

METHODS

Summary-level data for COPD were obtained from two national biobanks: the UK Biobank, comprising 1,605 cases and 461,328 controls, and FinnGen, with 6,915 cases and 186,723 controls. We also validated our findings using clinical data from 2,690 COPD patients and 3,357 healthy controls from the First Affiliated Hospital of Guangzhou Medical University. A total of 44 bodily fluid biomarkers were selected as candidate risk factors. Mendelian randomization (MR) and meta-analyses were used to evaluate the causal effects of these bodily fluids on COPD and lung function (FEV1/FVC).

RESULTS

Mendelian randomization (MR) and meta-analyses, by integrating data from the UK Biobank and FinnGen cohort, found that 3 bodily fluids indicators (HDLC, EOS, and TP) were causally associated with the risk of COPD, two (EOS and TP) of which is consistent with our observational findings. Moreover, we noticed EOS and TP were causally associated with the risk of lung function (FEV1/FVC).

CONCLUSIONS

The MR findings and clinical data highlight the independent and significant roles of EOS and TP in the development of COPD and lung function (FEV1/FVC), which might provide a deeper insight into COPD risk factors and supply potential preventative strategies.

Causal relationship between the blood immune cells and intervertebral disc degeneration: univariable, bidirectional and multivariable Mendelian randomization

Background

Intervertebral disc degeneration (IVDD) is a prominent contributor to chronic low back pain, impacting millions of individuals annually. Current research on disc degeneration is placing a growing emphasis on the role of the immune system in this process. Nevertheless, the precise relationship between immunity and disc degeneration remains to be fully elucidated.

Method

We obtained GWAS data for immune cells from the latest summary-level GWAS, including 6,620 individuals from Sardinian and 746,667 individuals from five global populations. Summary results for IVDD were sourced from the FinnGen consortium, comprising 20,001 cases and 164,682 controls. We conducted a comprehensive univariable Mendelian randomization (MR) analysis to explore the potential causal relationship between immune cells and IVDD. Primary estimation was carried out using Inverse-Variance Weighting (IVW). To ensure robustness, we employed additional MR methods such as MR-Egger, Weighted Median, Weighted Mode, and Simple Mode. Various tests were employed to assess pleiotropy and heterogeneity, including the Cochran Q test, leave-one-out test, MR-Egger intercept analysis and MR-PRESSO test. To account for potential confounding factors among the immune cells, we conducted a multivariable MR analysis. Finally, we investigated the possibility of a reverse association between immune cells and IVDD through bidirectional MR.

Result

In total, our study identified 15 immune cells significantly associated with IVDD through univariable MR. Among these, 9 immune cell types were indicated as potential contributors to IVDD, while 6 were found to have protective effects. Importantly, we observed no evidence of heterogeneity or pleiotropy, signifying the robustness of our results. To mitigate confounding among immune cells, we utilized multivariable MR, leading to the discovery that only 9 immune cell types exerted independent effects on IVDD. These encompassed 7 as risk factors and 2 as protective factors. Additionally, our analysis revealed a bidirectional causal relationship between CD39+ CD4+ T cell %CD4+ T cell and IVDD.

Conclusion

Our findings suggest a connection between immune cells and the risk of IVDD, shedding light on potential therapeutic avenues for modulating immune cell function in individuals with IVDD. However, the specific underlying mechanisms warrant further investigation in future experiments.

COLOCdb: a comprehensive resource for multi-model colocalization of complex traits

Large-scale genome-wide association studies (GWAS) have provided profound insights into complex traits and diseases. Yet, deciphering the fine-scale molecular mechanisms of how genetic variants manifest to cause the phenotypes remains a daunting task. Here, we present COLOCdb (https://ngdc.cncb.ac.cn/colocdb), a comprehensive genetic colocalization database by integrating more than 3000 GWAS summary statistics and 13 types of xQTL to date. By employing two representative approaches for the colocalization analysis, COLOCdb deposits results from three key components: (i) GWAS-xQTL, pair-wise colocalization between GWAS loci and different types of xQTL, (ii) GWAS-GWAS, pair-wise colocalization between the trait-associated genetic loci from GWASs and (iii) xQTL-xQTL, pair-wise colocalization between the genetic loci associated with molecular phenotypes in xQTLs. These results together represent the most comprehensive colocalization analysis, which also greatly expands the list of shared variants with genetic pleiotropy. We expect that COLOCdb can serve as a unique and useful resource in advancing the discovery of new biological mechanisms and benefit future functional studies.

Genetic Overlap Between Global Cortical Brain Structure, C-Reactive Protein, and White Blood Cell Counts

BACKGROUND

For many brain disorders, a subset of patients jointly exhibit alterations in cortical brain structure and elevated levels of circulating immune markers. This may be driven in part by shared genetic architecture. Therefore, we investigated the phenotypic and genetic associations linking global cortical surface area and thickness with blood immune markers (i.e., white blood cell counts and plasma C-reactive protein levels).

METHODS

Linear regression was used to assess phenotypic associations in 30,823 UK Biobank participants. Genome-wide and local genetic correlations were assessed using linkage disequilibrium score regression and local analysis of covariance annotation. The number of shared trait-influencing genetic variants was estimated using MiXeR. Shared genetic architecture was assessed using a conjunctional false discovery rate framework, and mapped genes were included in gene-set enrichment analyses.

RESULTS

Cortical structure and blood immune markers exhibited predominantly inverse phenotypic associations. There were modest genome-wide genetic correlations, the strongest of which were for C-reactive protein levels (rg_surface_area = -0.13, false discovery rate-corrected p = 4.17 × 10-3; rg_thickness = -0.13, false discovery rate-corrected p = 4.00 × 10-2). Meanwhile, local genetic correlations showed a mosaic of positive and negative associations. White blood cells shared on average 46.24% and 38.64% of trait-influencing genetic variants with surface area and thickness, respectively. Additionally, surface area shared 55 unique loci with the blood immune markers while thickness shared 15. Overall, monocyte count exhibited the largest genetic overlap with cortical brain structure. A series of gene enrichment analyses implicated neuronal-, astrocytic-, and schizophrenia-associated genes.

CONCLUSIONS

The findings indicate shared genetic underpinnings for cortical brain structure and blood immune markers, with implications for neurodevelopment and understanding the etiology of brain-related disorders.

Perturbational phenotyping of human blood cells reveals genetically determined latent traits associated with subsets of common diseases

Although genome-wide association studies (GWAS) have successfully linked genetic risk loci to various disorders, identifying underlying cellular biological mechanisms remains challenging due to the complex nature of common diseases. We established a framework using human peripheral blood cells, physical, chemical and pharmacological perturbations, and flow cytometry-based functional readouts to reveal latent cellular processes and performed GWAS based on these evoked traits in up to 2,600 individuals. We identified 119 genomic loci implicating 96 genes associated with these cellular responses and discovered associations between evoked blood phenotypes and subsets of common diseases. We found a population of pro-inflammatory anti-apoptotic neutrophils prevalent in individuals with specific subsets of cardiometabolic disease. Multigenic models based on this trait predicted the risk of developing chronic kidney disease in type 2 diabetes patients. By expanding the phenotypic space for human genetic studies, we could identify variants associated with large effect response differences, stratify patients and efficiently characterize the underlying biology.

Lymphocyte Count Derived Polygenic Score and Interindividual Variability in CD4 T-cell Recovery in Response to Antiretroviral Therapy

Access to safe and effective antiretroviral therapy (ART) is a cornerstone in the global response to the HIV pandemic. Among people living with HIV, there is considerable interindividual variability in absolute CD4 T-cell recovery following initiation of virally suppressive ART. The contribution of host genetics to this variability is not well understood. We explored the contribution of a polygenic score which was derived from large, publicly available summary statistics for absolute lymphocyte count from individuals in the general population (PGSlymph) due to a lack of publicly available summary statistics for CD4 T-cell count. We explored associations with baseline CD4 T-cell count prior to ART initiation (n=4959) and change from baseline to week 48 on ART (n=3274) among treatment-naïve participants in prospective, randomized ART studies of the AIDS Clinical Trials Group. We separately examined an African-ancestry-derived and a European-ancestry-derived PGSlymph, and evaluated their performance across all participants, and also in the African and European ancestral groups separately. Multivariate models that included PGSlymph, baseline plasma HIV-1 RNA, age, sex, and 15 principal components (PCs) of genetic similarity explained ∼26-27% of variability in baseline CD4 T-cell count, but PGSlymph accounted for <1% of this variability. Models that also included baseline CD4 T-cell count explained ∼7-9% of variability in CD4 T-cell count increase on ART, but PGSlymph accounted for <1% of this variability. In univariate analyses, PGSlymph was not significantly associated with baseline or change in CD4 T-cell count. Among individuals of African ancestry, the African PGSlymph term in the multivariate model was significantly associated with change in CD4 T-cell count while not significant in the univariate model. When applied to lymphocyte count in a general medical biobank population (Penn Medicine BioBank), PGSlymph explained ∼6-10% of variability in multivariate models (including age, sex, and PCs) but only ∼1% in univariate models. In summary, a lymphocyte count PGS derived from the general population was not consistently associated with CD4 T-cell recovery on ART. Nonetheless, adjusting for clinical covariates is quite important when estimating such polygenic effects.

Circulating white blood cell traits and colorectal cancer risk: A Mendelian randomisation study

Observational studies have suggested a protective role for eosinophils in colorectal cancer (CRC) development and implicated neutrophils, but the causal relationships remain unclear. Here, we aimed to estimate the causal effect of circulating white blood cell (WBC) counts (N = ~550 000) for basophils, eosinophils, monocytes, lymphocytes and neutrophils on CRC risk (N = 52 775 cases and 45 940 controls) using Mendelian randomisation (MR). For comparison, we also examined this relationship using individual-level data from UK Biobank (4043 incident CRC cases and 332 773 controls) in a longitudinal cohort analysis. The inverse-variance weighted (IVW) MR analysis suggested a protective effect of increased basophil count and eosinophil count on CRC risk [OR per 1-SD increase: 0.88, 95% CI: 0.78-0.99, P = .04; OR: 0.93, 95% CI: 0.88-0.98, P = .01]. The protective effect of eosinophils remained [OR per 1-SD increase: 0.88, 95% CI: 0.80-0.97, P = .01] following adjustments for all other WBC subtypes, to account for genetic correlation between the traits, using multivariable MR. A protective effect of increased lymphocyte count on CRC risk was also found [OR: 0.84, 95% CI: 0.76-0.93, P = 6.70e-4] following adjustment. Consistent with MR results, a protective effect for eosinophils in the cohort analysis in the fully adjusted model [RR per 1-SD increase: 0.96, 95% CI: 0.93-0.99, P = .02] and following adjustment for the other WBC subtypes [RR: 0.96, 95% CI: 0.93-0.99, P = .001] was observed. Our study implicates peripheral blood immune cells, in particular eosinophils and lymphocytes, in CRC development, highlighting a need for mechanistic studies to interrogate these relationships.