Shape-IT: new rapid and accurate algorithm for haplotype inference

Identifying individuals with rare disease variants by inferring shared ancestral haplotypes from SNP array data

We describe FoundHaplo, an identity-by-descent algorithm that can be used to screen untyped disease-causing variants using single nucleotide polymorphism (SNP) array data. FoundHaplo leverages knowledge of shared disease haplotypes for inherited variants to identify those who share the disease haplotype and are, therefore, likely to carry the rare [minor allele frequency (MAF) ≤ 0.01%] variant. We performed a simulation study to evaluate the performance of FoundHaplo across 33 disease-harbouring loci. FoundHaplo was used to infer the presence of two rare (MAF ≤ 0.01%) pathogenic variants, SCN1B c.363C>G (p.Cys121Trp) and WWOX c.49G>A (p.E17K), which can cause mild dominant and severe recessive epilepsy, respectively, in the Epi25 cohort and the UK Biobank. FoundHaplo demonstrated substantially better sensitivity at inferring the presence of these rare variants than existing genome-wide imputation. FoundHaplo is a valuable screening tool for searching disease-causing variants with known founder effects using only SNP genotyping data. It is also applicable to nonhuman applications and nondisease-causing traits, including rare-variant drivers of quantitative traits. The FoundHaplo algorithm is available at https://github.com/bahlolab/FoundHaplo (DOI:10.5281/zenodo.8058286).

Characterization of NAT, GST, and CYP2E1 Genetic Variation in Sub-Saharan African Populations: Implications for Treatment of Tuberculosis and Other Diseases

Tuberculosis (TB) is a major health burden in Africa. Although TB is treatable, anti-TB drugs are associated with adverse drug reactions (ADRs), which are partly attributed to pharmacogenetic variation. The distribution of star alleles (haplotypes) influencing anti-TB drug metabolism is unknown in many African populations. This presents challenges in implementing genotype-guided therapy in Africa to decrease the occurrence of ADRs and enhance the efficacy of anti-TB drugs. In this study, we used StellarPGx to call variants and star alleles in NAT1, NAT2, GSTM1, GSTT1, GSTP1, and CYP2E1, from 1079 high-depth African whole genomes. We present the distribution of common, rare, and potential novel star alleles across various Sub-Saharan African (SSA) populations, in comparison with other global populations. NAT1*10 (53.6%), GSTT1*0 (65%), GSTM1*0 (48%), and NAT2*5 (17.5%) were among the predominant functionally relevant star alleles. Additionally, we predicted varying phenotype distributions for NAT1 and NAT2 (acetylation) and the glutathione-S-transferase (GST) enzymes (detoxification activity) between SSA and other global populations. Forty-seven potentially novel haplotypes were identified computationally across the genes. This study provides insight into the distribution of key variants and star alleles potentially relevant to anti-TB drug metabolism and other drugs prescribed across various African populations. The high number of potentially novel star alleles exemplifies the need for pharmacogenomics studies in the African context. Overall, our study provides a foundation for functional pharmacogenetic studies and potential implementation of pharmacogenetic testing in Africa to reduce the risk of ADRs related to treatment of TB and other diseases.

Integrative Computational Analysis of Common EXO5 Haplotypes: Impact on Protein Dynamics, Genome Stability, and Cancer Progression

Understanding the impact of common germline variants on protein structure, function, and disease progression is crucial in cancer research. This study presents a comprehensive analysis of the EXO5 gene, which encodes a DNA exonuclease involved in DNA repair that was previously associated with cancer susceptibility. We employed an integrated approach combining genomic and clinical data analysis, deep learning variant effect prediction, and molecular dynamics (MD) simulations to investigate the effects of common EXO5 haplotypes on protein structure, dynamics, and cancer outcomes. We characterized the haplotype structure of EXO5 across diverse human populations, identifying five common haplotypes, and studied their impact on the EXO5 protein. Extensive, all-atom MD simulations revealed significant structural and dynamic differences among the EXO5 protein variants, particularly in their catalytic region. The L151P EXO5 protein variant exhibited the most substantial conformational changes, potentially disruptive for EXO5's function and nuclear localization. Analysis of The Cancer Genome Atlas data showed that cancer patients carrying L151P EXO5 had significantly shorter progression-free survival in prostate and pancreatic cancers and exhibited increased genomic instability. This study highlights the strength of our methodology in uncovering the effects of common genetic variants on protein function and their implications for disease outcomes.

A structural haplotype in the 17q21.31 MAPT region is associated with increased risk for chronic traumatic encephalopathy endophenotypes

Chronic traumatic encephalopathy (CTE) is a neurodegenerative tauopathy associated with repetitive head impact (RHI) exposure. Genetic variation in the 17q21.31 region, containing microtubule-associated protein tau (MAPT), has been implicated in tauopathies but has not been investigated in CTE. The region includes a megabase-long inversion (H1/H2) and copy-number variations, including α, β, and γ segments, which can be characterized as nine segregating structural haplotypes. We leveraged array SNP data and a reference panel across the 17q21.31 region to impute structural haplotypes and test their association with CTE endophenotypes in 447 European ancestry brain donors with RHI exposure. The H1β1γ1 haplotype was significantly associated with dementia and semi-quantitative tau burden in multiple cortical and medial temporal regions commonly affected in CTE. H1β1γ1 differential expression analyses in dorsolateral frontal cortex implicated cis-acting genes and inflammatory pathways. Taken together, the H1β1γ1 haplotype may help explain CTE heterogeneity among those with similar RHI exposure.

The impact of the COVID-19 pandemic on alcohol use disorder symptoms: Testing interactions with polygenic risk

Objective: The purpose of this study was to test whether COVID impact interacts with genetic risk (polygenic risk score/PRS) to predict alcohol use disorder (AUD) symptoms. Method: Participants were n = 455 college students (79.6% female, 51% European Ancestry/EA, 24% African Ancestry/AFR, 25% Americas Ancestry/AMER) from a longitudinal study during the initial stage (March-May 2020) of the pandemic. Path models allowed for the examination of PRS and previously identified COVID-19 impact constructs. Results: There was a main effect of the AUD PRS on AUD symptoms within the EA group (β: .165, p < .01). Additionally, food/housing insecurity was predictive in the AMER group (β.295, p < .05), and greater increases in substance use were associated with AUD symptoms for EA (β:.459, p < .001) and AMER groups (β:.468, p < .001). Conclusions: Greater food/housing instability and increases in substance use, as well higher scores on PRS are associated with more AUD symptoms for some ancestral groups within this college sample.

Red blood cell urate levels are linked to hemolysis in vitro and post-transfusion as a function of donor sex, population and genetic polymorphisms in SLC2A9 and ABCG2

BACKGROUND

Storage of packed red blood cells (RBCs) for transfusion leads to biochemical and morphological changes, increasing hemolysis risk. Urate levels in blood bags at donation contribute to the molecular heterogeneity and hemolytic propensity of stored RBCs. However, studies to date have been underpowered to investigate at scale the contribution of donor demographics and genetics to the heterogeneity in urate levels across donations.

STUDY DESIGN AND METHODS

Urate levels were measured in 13,091 RBC units from the REDS study. Characteristics tested included hemolysis parameters (spontaneous, osmotic, oxidative) at storage end and post-transfusion hemoglobin (Hb) increments in recipients. Donor demographics, urate levels, and genetic variants were analyzed for associations with these outcomes.

RESULTS

Elevated urate levels were linked to male sex, older age, high BMI, and Asian descent. Units with high urate levels exhibited increased spontaneous and osmotic hemolysis, while oxidative hemolysis was unaffected. Genetic variants in SLC2A9 (V282I) and ABCG2 (Q141K) were strongly associated with elevated urate, particularly in Asian donors. Post-transfusion analyses revealed that units from female donors carrying these variants were associated with reduced Hb increments, with up to a 31% reduction in efficacy. This effect was not observed in male donors.

DISCUSSION

RBC urate levels and genetic traits significantly impact storage quality and transfusion outcomes. These findings highlight the importance of donor molecular characteristics for optimizing transfusion strategies. Moreover, genetic and metabolic insights may inform donor recruitment efforts, providing health feedback to volunteers while ensuring effective transfusion products.

A meta-analysis of genome-wide studies of resilience in the German population

Resilience is the capacity to adapt to stressful life events. As such, this trait is associated with physical and mental functions and conditions. Here, we aimed to identify the genetic factors contributing to shape resilience. We performed variant- and gene-based meta-analyses of genome-wide association studies from six German cohorts (N = 15822) using the 11-item version of the Resilience Scale (RS-11) as outcome measure. Variant- and gene-level results were combined to explore the biological context using network analysis. In addition, we conducted tests of correlation between RS-11 and the polygenic scores (PGSs) for 12 personality and mental health traits in one of these cohorts (PROCAM-2, N = 3879). The variant-based analysis found no signals associated with resilience at the genome-wide level (p < 5 × 10-8), but suggested five genomic loci (p < 1 × 10-5). The gene-based analysis identified three genes (ROBO1, CIB3 and LYPD4) associated with resilience at genome-wide level (p < 2.48 × 10-6) and 32 potential candidates (p < 1 × 10-4). Network analysis revealed enrichment of biological pathways related to neuronal proliferation and differentiation, synaptic organization, immune responses and vascular homeostasis. We also found significant correlations (FDR < 0.05) between RS-11 and the PGSs for neuroticism and general happiness. Overall, our observations suggest low heritability of resilience. Large, international efforts will be required to uncover the genetic factors that contribute to shape trait resilience. Nevertheless, as the largest investigation of the genetics of resilience in general population to date, our study already offers valuable insights into the biology potentially underlying resilience and resilience's relationship with other personality traits and mental health.

Increased exercise tolerance in humanized G6PD-deficient mice

ABSTRACT

Glucose-6-phosphate dehydrogenase (G6PD) deficiency affects 500 million people globally, affecting red blood cell (RBC) antioxidant pathways and increasing susceptibility to hemolysis under oxidative stress. Despite the systemic generation of reactive oxygen species during exercise, the effects of exercise on individuals with G6PD deficiency remain poorly understood This study used humanized mouse models expressing the G6PD Mediterranean variant (S188F, with 10% enzymatic activity) to investigate exercise performance and molecular outcomes. Surprisingly, despite decreased enzyme activity, G6PD-deficient mice have faster critical speed than mice expressing human canonical G6PD. After exercise, deficient mice did not exhibit differences in RBC morphology or hemolysis, but had improved cardiac function, including cardiac output, stroke volume, sarcomere length, and mitochondrial content. Proteomics analyses of cardiac and skeletal muscles (gastrocnemius, soleus) from G6PD-deficient compared with sufficient mice revealed improvements in mitochondrial function and increased protein turnover via ubiquitination, especially for mitochondrial and structural myofibrillar proteins. Mass spectrometry-based metabolomics revealed alterations in energy metabolism and fatty acid oxidation. These findings challenge the traditional assumptions regarding hemolytic risk during exercise in G6PD deficiency, suggesting a potential metabolic advantage in exercise performance for individuals carrying noncanonical G6PD variants.

Neuropathology-based approach reveals novel Alzheimer's Disease genes and highlights female-specific pathways and causal links to disrupted lipid metabolism: insights into a vicious cycle

Dementia refers to an umbrella phenotype of many different underlying pathologies with Alzheimer's disease (AD) being the most common type. Neuropathological examination remains the gold standard for accurate AD diagnosis, however, most that we know about AD genetics is based on Genome-Wide Association Studies (GWAS) of clinically defined AD. Such studies have identified multiple AD susceptibility variants with a significant portion of the heritability unexplained and highlighting the phenotypic and genetic heterogeneity of the clinically defined entity. Furthermore, despite women's increased susceptibility to dementia, there is a lack of sex-specific genetic studies and understanding of sex-specific background for the disorder. Here, we aim to tackle the heterogeneity of AD by specifically concentrating on neuropathological features and pursuing sex-specific analysis. We bring together 14 different genomic and neuropathology datasets (6960 individuals) and we integrate our GWAS findings with transcriptomic and phenotypic data aiming to also identify biomarkers for AD progression. We uncover novel genetic associations to AD neuropathology, including BIN1 and OPCML. Our sex-specific analysis points to a role for BIN1 specifically in women as well as novel AD loci including QRFPR and SGCZ. Post-GWAS analyses illuminate the functional and biological mechanisms underlying AD and reveal sex-specific differences. Finally, through PheWAS and Mendelian Randomization analysis, we identify causal links with AD neuropathology pointing to disrupted lipid metabolism, as well as impaired peripheral immune response and liver dysfunction as part of a vicious cycle that fuels neurodegeneration.

Genomic signatures of adaptation in native lizards exposed to human-introduced fire ants

Understanding the process of genetic adaptation in response to human-mediated ecological change will help elucidate the eco-evolutionary impacts of human activity. In the 1930s red imported fire ants (Solenopsis invicta) were accidently introduced to the Southeastern USA, where today they are both venomous predators and toxic prey to native eastern fence lizards (Sceloporus undulatus). Here, we investigate potential lizard adaptation to invasive fire ants by generating whole-genome sequences from 420 lizards across three populations: one with long exposure to fire ants, and two unexposed populations. Signatures of positive selection exclusive to the exposed population overlap immune system, growth factor pathway, and morphological development genes. Among invaded lizards, longer limbs (used to remove stinging ants) are associated with increased survival. We identify alleles associated with longer limbs that are highly differentiated from the unexposed populations, a pattern counter to the pre-invasion latitudinal cline for limb lengths based on museum specimens. While we cannot rule out other environmental differences between populations driving these patterns, these results do constitute plausible genetic adaptations in lizards invaded by fire ants.

Polygenic Risk Scores and Twin Concordance for Schizophrenia and Bipolar Disorder

Importance

Schizophrenia and bipolar disorder are highly heritable psychiatric disorders with strong genetic and phenotypic overlap. Twin and molecular methods can be leveraged to predict the shared genetic liability to these disorders.

Objective

To investigate whether twin concordance for psychosis depends on the level of polygenic risk score (PRS) for psychosis and zygosity and compare PRS from cases and controls from several large samples and estimate the twin heritability of psychosis.

Design, Setting, and Participants

In this case-control study, psychosis PRS were generated from a genome-wide association study (GWAS) combining schizophrenia and bipolar disorder into a single psychosis phenotype and compared between cases and controls from the Schizophrenia and Bipolar Twin Study in Sweden (STAR) project. Further tests were conducted to ascertain if twin concordance for psychosis depended on the mean PRS for psychosis. Structural equation modeling was used to estimate heritability. This study constituted an analysis of existing clinical and population datasets with genotype and/or twin data. Included were twins from the STAR cohort and from the Swedish Twin Registry. Data were collected during the 2006 to 2013 period and analyzed from March 2023 to June 2024.

Exposures

PRS for psychosis based on the most recent GWAS of combined schizophrenia/bipolar disorder.

Main Outcomes and Measures

Psychosis case status was assessed by clinical interviews and/or Swedish National Register data.

Results

The final cohort comprised 87 pairs of twins with 1 or both affected and 59 unaffected pairs from the STAR project (for a total of 292 twins) as well as 443 pairs with 1 or both affected and 20 913 unaffected pairs from the Swedish Twin Registry. Among the 292 twins (mean [SD] birth year, 1960 [10.8] years; 158 female [54.1%]; 134 male [45.9%]), 134 were monozygotic twins, and 158 were dyzygotic twins. PRS for psychosis was higher in cases than in controls and associated with twin concordance for psychosis (1-SD increase in PRS, odds ratio [OR], 2.12; 95% CI, 1.23-3.87 on case status in monozygotic twins and OR, 2.74; 95% CI, 1.56-5.30 in dizygotic twins). The association between PRS for psychosis and concordance was not modified by zygosity. The twin heritability was estimated at 0.73 (95% CI, 0.30-1.00), which overlapped with the estimate in the full Swedish Twin Registry (0.69; 95% CI, 0.43-0.85).

Conclusions and Relevance

In this case-control study, using the natural experiment of twins, results suggest that twins with greater inherited liability for psychosis were more likely to have an affected co-twin. Results from twin and molecular designs largely aligned. Even as illness vulnerability is not solely genetic, PRS carried predictive power for psychosis even in a modest sample size.

Biological and genetic determinants of glycolysis: Phosphofructokinase isoforms boost energy status of stored red blood cells and transfusion outcomes

Mature red blood cells (RBCs) lack mitochondria and thus exclusively rely on glycolysis to generate adenosine triphosphate (ATP) during aging in vivo or storage in blood banks. Here, we leveraged 13,029 volunteers from the Recipient Epidemiology and Donor Evaluation Study to identify associations between end-of-storage levels of glycolytic metabolites and donor age, sex, and ancestry-specific genetic polymorphisms in regions encoding phosphofructokinase 1, platelet (detected in mature RBCs); hexokinase 1 (HK1); and ADP-ribosyl cyclase 1 and 2 (CD38/BST1). Gene-metabolite associations were validated in fresh and stored RBCs from 525 Diversity Outbred mice and via multi-omics characterization of 1,929 samples from 643 human RBC units during storage. ATP and hypoxanthine (HYPX) levels-and the genetic traits linked to them-were associated with hemolysis in vitro and in vivo, both in healthy autologous transfusion recipients and in 5,816 critically ill patients receiving heterologous transfusions, suggesting their potential as markers to improve transfusion outcomes.

Haplotype analysis identifies functional elements in monoclonal gammopathy of unknown significance

Genome-wide association studies (GWASs) based on common single nucleotide polymorphisms (SNPs) have identified several loci associated with the risk of monoclonal gammopathy of unknown significance (MGUS), a precursor condition for multiple myeloma (MM). We hypothesized that analyzing haplotypes might be more useful than analyzing individual SNPs, as it could identify functional chromosomal units that collectively contribute to MGUS risk. To test this hypothesis, we used data from our previous GWAS on 992 MGUS cases and 2910 controls from three European populations. We identified 23 haplotypes that were associated with the risk of MGUS at the genome-wide significance level (p < 5 × 10-8) and showed consistent results among all three populations. In 10 genomic regions, strong promoter, enhancer and regulatory element-related histone marks and their connections to target genes as well as genome segmentation data supported the importance of these regions in MGUS susceptibility. Several associated haplotypes affected pathways important for MM cell survival such as ubiquitin-proteasome system (RNF186, OTUD3), PI3K/AKT/mTOR (HINT3), innate immunity (SEC14L1, ZBP1), cell death regulation (BID) and NOTCH signaling (RBPJ). These pathways are important current therapeutic targets for MM, which may highlight the advantage of the haplotype approach homing to functional units.

Swordtail fish hybrids reveal that genome evolution is surprisingly predictable after initial hybridization

Over the past 2 decades, biologists have come to appreciate that hybridization, or genetic exchange between distinct lineages, is remarkably common-not just in particular lineages but in taxonomic groups across the tree of life. As a result, the genomes of many modern species harbor regions inherited from related species. This observation has raised fundamental questions about the degree to which the genomic outcomes of hybridization are repeatable and the degree to which natural selection drives such repeatability. However, a lack of appropriate systems to answer these questions has limited empirical progress in this area. Here, we leverage independently formed hybrid populations between the swordtail fish Xiphophorus birchmanni and X. cortezi to address this fundamental question. We find that local ancestry in one hybrid population is remarkably predictive of local ancestry in another, demographically independent hybrid population. Applying newly developed methods, we can attribute much of this repeatability to strong selection in the earliest generations after initial hybridization. We complement these analyses with time-series data that demonstrates that ancestry at regions under selection has remained stable over the past approximately 40 generations of evolution. Finally, we compare our results to the well-studied X. birchmanni × X. malinche hybrid populations and conclude that deeper evolutionary divergence has resulted in stronger selection and higher repeatability in patterns of local ancestry in hybrids between X. birchmanni and X. cortezi.

Validation of the Utility of the Genetically Shared Regions of Chromosomes (GD-ICS) Measuring Method in Identifying Complicated Genetic Relatedness

BACKGROUND

Relatives share more genomic regions than unrelated individuals, with closer relatives sharing more regions. This concept, paired with the increased availability of high-throughput single nucleotide polymorphism (SNP) genotyping technologies, has made it feasible to measure the shared chromosomal regions between individuals to assess their level of relation to each other. However, such techniques have remained in the conceptual rather than practical stages in terms of applying measures or indices. Recently, we developed an index called "genetic distance-based index of chromosomal sharing (GD-ICS)" utilizing large-scale SNP data from Korean family samples and demonstrated its potential for practical applications in kinship determination. In the current study, we present validation results from various real cases demonstrating the utility of this method in resolving complex familial relationships where information obtained from traditional short tandem repeats (STRs) or lineage markers is inconclusive.

METHODS

We obtained large-scale SNP data through microarray analysis from Korean individuals involving 13 kinship cases and calculated GD-ICS values using the method described in our previous study. Based on the GD-ICS reference constructed for Korean families, each disputed kinship was evaluated and validated using a combination of traditional STRs and lineage markers.

RESULTS

The cases comprised those A) that were found to be inconclusive using the traditional approach, B) for which it was difficult to apply traditional testing methods, and C) that were more conclusively resolved using the GD-ICS method. This method has overcome the limitations faced by traditional STRs in kinship testing, particularly in a paternity case with STR mutational events and in confirming distant kinship where the individual of interest is unavailable for testing. It has also been demonstrated to be effective in identifying various relationships without specific presumptions and in confirming a lack of genetic relatedness between individuals.

CONCLUSION

This method has been proven effective in identifying familial relationships across diverse complex and practical scenarios. It is not only useful when traditional testing methods fail to provide conclusive results, but it also enhances the resolution of challenging kinship cases, which suggests its applicability in various types of practical casework.

Key subphenotypes of bipolar disorder are differentially associated with polygenic liabilities for bipolar disorder, schizophrenia, and major depressive disorder

Bipolar disorder (BD) features heterogenous clinical presentation and course of illness. It remains unclear how subphenotypes associate with genetic loadings of BD and related psychiatric disorders. We investigated associations between the subphenotypes and polygenic risk scores (PRS) for BD, schizophrenia, and major depressive disorder (MDD) in two BD cohorts from Sweden (N = 5180) and the UK (N = 2577). Participants were assessed through interviews and medical records for inter-episode remission, psychotic features during mood episodes, global assessment of functioning (GAF, function and symptom burden dimensions), and comorbid anxiety disorders. Meta-analyses based on both cohorts showed that inter-episode remission and GAF-function were positively correlated with BD-PRS but negatively correlated with schizophrenia-PRS (SCZ-PRS) and MDD-PRS. Moreover, BD-PRS was negatively, and MDD-PRS positively, associated with the risk of comorbid anxiety disorders. Finally, SCZ-PRS was positively associated with psychotic symptoms during mood episodes. Assuming a higher PRS of certain psychiatric disorders in cases with a positive family history, we further tested the associations between subphenotypes in index BD people and occurrence of BD, schizophrenia, or MDD in their relatives using Swedish national registries. BD patients with a relative diagnosed with BD had: (1) higher GAF and lower risk of comorbid anxiety than those with a relative diagnosed with schizophrenia or MDD, (2) lower risk of psychotic symptoms than those with a relative diagnosed with schizophrenia. Our findings shed light on the genetic underpinnings of the heterogeneity in clinical manifestations and course of illness in BD, which ultimately provide insights for developing personalized approaches to the diagnosis and treatment.

Genetic regulation of carnitine metabolism controls lipid damage repair and aging RBC hemolysis in vivo and in vitro

ABSTRACT

Recent large-scale multiomics studies suggest that genetic factors influence the chemical individuality of donated blood. To examine this concept, we performed metabolomics analyses of 643 blood units from volunteers who donated units of packed red blood cells (RBCs) on 2 separate occasions. These analyses identified carnitine metabolism as the most reproducible pathway across multiple donations from the same donor. We also measured l-carnitine and acyl-carnitines in 13 091 packed RBC units from donors in the Recipient Epidemiology and Donor Evaluation study. Genome-wide association studies against 879 000 polymorphisms identified critical genetic factors contributing to interdonor heterogeneity in end-of-storage carnitine levels, including common nonsynonymous polymorphisms in genes encoding carnitine transporters (SLC22A16, SLC22A5, and SLC16A9); carnitine synthesis (FLVCR1 and MTDH) and metabolism (CPT1A, CPT2, CRAT, and ACSS2), and carnitine-dependent repair of lipids oxidized by ALOX5. Significant associations between genetic polymorphisms on SLC22 transporters and carnitine pools in stored RBCs were validated in 525 Diversity Outbred mice. Donors carrying 2 alleles of the rs12210538 SLC22A16 single-nucleotide polymorphism exhibited the lowest l-carnitine levels, significant elevations of in vitro hemolysis, and the highest degree of vesiculation, accompanied by increases in lipid peroxidation markers. Separation of RBCs by age, via in vivo biotinylation in mice, and Percoll density gradients of human RBCs, showed age-dependent depletions of l-carnitine and acyl-carnitine pools, accompanied by progressive failure of the reacylation process after chemically induced membrane lipid damage. Supplementation of stored murine RBCs with l-carnitine boosted posttransfusion recovery, suggesting this could represent a viable strategy to improve RBC storage quality.

Biological and Genetic Determinants of Glycolysis: Phosphofructokinase Isoforms Boost Energy Status of Stored Red Blood Cells and Transfusion Outcomes

Mature red blood cells (RBCs) lack mitochondria, and thus exclusively rely on glycolysis to generate adenosine triphosphate (ATP) during aging in vivo or storage in the blood bank. Here we leveraged 13,029 volunteers from the Recipient Epidemiology and Donor Evaluation Study to identify an association between end-of-storage levels of glycolytic metabolites and donor age, sex, and ancestry-specific genetic polymorphisms in regions encoding phosphofructokinase 1, platelet (detected in mature RBCs), hexokinase 1, ADP-ribosyl cyclase 1 and 2 (CD38/BST1). Gene-metabolite associations were validated in fresh and stored RBCs from 525 Diversity Outbred mice, and via multi-omics characterization of 1,929 samples from 643 human RBC units during storage. ATP and hypoxanthine levels - and the genetic traits linked to them - were associated with hemolysis in vitro and in vivo, both in healthy autologous transfusion recipients and in 5,816 critically ill patients receiving heterologous transfusions, suggesting their potential as markers to improve transfusion outcomes.

eTOC and Highlights

Highlights

Blood donor age and sex affect glycolysis in stored RBCs from 13,029 volunteers;Ancestry, genetic polymorphisms in PFKP, HK1, CD38/BST1 influence RBC glycolysis;Modeled PFKP effects relate to preventing loss of the total AXP pool in stored RBCs;ATP and hypoxanthine are biomarkers of hemolysis in vitro and in vivo.

Ferroptosis regulates hemolysis in stored murine and human red blood cells

Red blood cell (RBC) metabolism regulates hemolysis during aging in vivo and in the blood bank. Here, we leveraged a diversity outbred mouse population to map the genetic drivers of fresh/stored RBC metabolism and extravascular hemolysis upon storage and transfusion in 350 mice. We identify the ferrireductase Steap3 as a critical regulator of a ferroptosis-like process of lipid peroxidation. Steap3 polymorphisms were associated with RBC iron content, in vitro hemolysis, and in vivo extravascular hemolysis both in mice and 13,091 blood donors from the Recipient Epidemiology and Donor evaluation Study. Using metabolite Quantitative Trait Loci analyses, we identified a network of gene products (FADS1/2, EPHX2 and LPCAT3) - enriched in donors of African descent - associated with oxylipin metabolism in stored human RBCs and related to Steap3 or its transcriptional regulator, the tumor protein TP53. Genetic variants were associated with lower in vivo hemolysis in thousands of single-unit transfusion recipients.

Highlights

Steap3 regulates lipid peroxidation and extravascular hemolysis in 350 diversity outbred miceSteap3 SNPs are linked to RBC iron, hemolysis, vesiculation in 13,091 blood donorsmQTL analyses of oxylipins identified ferroptosis-related gene products FADS1/2, EPHX2, LPCAT3Ferroptosis markers are linked to hemoglobin increments in transfusion recipients.

Graphical abstract

Integrative multi-omics analyses to identify the genetic and functional mechanisms underlying ovarian cancer risk regions

To identify credible causal risk variants (CCVs) associated with different histotypes of epithelial ovarian cancer (EOC), we performed genome-wide association analysis for 470,825 genotyped and 10,163,797 imputed SNPs in 25,981 EOC cases and 105,724 controls of European origin. We identified five histotype-specific EOC risk regions (p value <5 × 10-8) and confirmed previously reported associations for 27 risk regions. Conditional analyses identified an additional 11 signals independent of the primary signal at six risk regions (p value <10-5). Fine mapping identified 4,008 CCVs in these regions, of which 1,452 CCVs were located in ovarian cancer-related chromatin marks with significant enrichment in active enhancers, active promoters, and active regions for CCVs from each EOC histotype. Transcriptome-wide association and colocalization analyses across histotypes using tissue-specific and cross-tissue datasets identified 86 candidate susceptibility genes in known EOC risk regions and 32 genes in 23 additional genomic regions that may represent novel EOC risk loci (false discovery rate <0.05). Finally, by integrating genome-wide HiChIP interactome analysis with transcriptome-wide association study (TWAS), variant effect predictor, transcription factor ChIP-seq, and motifbreakR data, we identified candidate gene-CCV interactions at each locus. This included risk loci where TWAS identified one or more candidate susceptibility genes (e.g., HOXD-AS2, HOXD8, and HOXD3 at 2q31) and other loci where no candidate gene was identified (e.g., MYC and PVT1 at 8q24) by TWAS. In summary, this study describes a functional framework and provides a greater understanding of the biological significance of risk alleles and candidate gene targets at EOC susceptibility loci identified by a genome-wide association study.

Genomic Signatures of Microgeographic Adaptation in Anopheles coluzzii Along an Anthropogenic Gradient in Gabon

Species distributed across heterogeneous environments often evolve locally adapted populations, but understanding how these persist in the presence of homogenizing gene flow remains puzzling. In Gabon, Anopheles coluzzii, a major African malaria mosquito is found along an ecological gradient, including a sylvatic population, away of any human presence. This study identifies into the genomic signatures of local adaptation in populations from distinct environments including the urban area of Libreville, and two proximate sites 10km apart in the La Lopé National Park (LLP), a village and its sylvatic neighborhood. Whole genome re-sequencing of 96 mosquitoes unveiled ∼ 5.7millions high-quality single nucleotide polymorphisms. Coalescent-based demographic analyses suggest an ∼ 8,000-year-old divergence between Libreville and La Lopé populations, followed by a secondary contact ( ∼ 4,000 ybp) resulting in asymmetric effective gene flow. The urban population displayed reduced effective size, evidence of inbreeding, and strong selection pressures for adaptation to urban settings, as suggested by the hard selective sweeps associated with genes involved in detoxification and insecticide resistance. In contrast, the two geographically proximate LLP populations showed larger effective sizes, and distinctive genomic differences in selective signals, notably soft-selective sweeps on the standing genetic variation. Although neutral loci and chromosomal inversions failed to discriminate between LLP populations, our findings support that microgeographic adaptation can swiftly emerge through selection on standing genetic variation despite high gene flow. This study contributes to the growing understanding of evolution of populations in heterogeneous environments amid ongoing gene flow and how major malaria mosquitoes adapt to human.

Significance

Anopheles coluzzii , a major African malaria vector, thrives from humid rainforests to dry savannahs and coastal areas. This ecological success is linked to its close association with domestic settings, with human playing significant roles in driving the recent urban evolution of this mosquito. Our research explores the assumption that these mosquitoes are strictly dependent on human habitats, by conducting whole-genome sequencing on An. coluzzii specimens from urban, rural, and sylvatic sites in Gabon. We found that urban mosquitoes show de novo genetic signatures of human-driven vector control, while rural and sylvatic mosquitoes exhibit distinctive genetic evidence of local adaptations derived from standing genetic variation. Understanding adaptation mechanisms of this mosquito is therefore crucial to predict evolution of vector control strategies.

Phenotypic Associations With the HMOX1 GT(n) Repeat in European Populations

Heme oxygenase 1 is a key enzyme in the management of heme in humans. A GT(n) repeat length in the heme oxygenase 1 gene (HMOX1) has been widely associated with a variety of phenotypes, including susceptibility to and outcomes in diabetes, cancer, infections, and neonatal jaundice. However, studies have generally been small and results inconsistent. In this study, we imputed the GT(n) repeat length in participants from 2 UK cohort studies (the UK Biobank study (n = 463,005; recruited in 2006-2010) and the Avon Longitudinal Study of Parents and Children (ALSPAC; n = 937; recruited in 1990-1991)), with the reliability of imputation tested in other cohorts (1000 Genomes Project, Human Genome Diversity Project, and Personal Genome Project UK). Subsequently, we measured the relationship between repeat length and previously identified associations (diabetes, chronic obstructive pulmonary disease, pneumonia, and infection-related mortality in the UK Biobank; neonatal jaundice in ALSPAC) and performed a phenomewide association study in the UK Biobank. Despite high-quality imputation (correlation between true repeat length and imputed repeat length > 0.9 in test cohorts), clinical associations were not identified in either the phenomewide association study or specific association studies. These findings were robust to definitions of repeat length and sensitivity analyses. Despite multiple smaller studies identifying associations across a variety of clinical settings, we could not replicate or identify any relevant phenotypic associations with the HMOX1 GT(n) repeat.

A novel approach of kinship determination based on the physical length of genetically shared regions of chromosomes

BACKGROUND

Determination of genetic relatedness between individuals plays a crucial role in resolving numerous civil cases involving familial relationships and in forensic investigation concerning missing persons. Short tandem repeats (STRs), known for their high degree of DNA polymorphism, have traditionally been the primary choice of DNA markers in genetic testing, but their application for kinships testing is limited to cases involving close kinship. SNPs have emerged as promising supplementary markers for kinship determination. Nevertheless, the challenging remains in discriminating between third-degree or more distant relatives, such as first cousins, using SNPs.

OBJECTIVE

To investigate a kinship analysis method for distant degree of familial relationships using high-density SNP data.

METHODS

A high-density SNP data from 337 individuals of Korean families using Affymetrix Axiom KORV1.0-96 Array was obtained for this study. SNPs were aligned by chromosomal positions, and identity-by-state (IBS) was determined, and then shared regions as consecutive SNPs with IBS of 1 or 2 were investigated. The physical lengths of these IBS segments were measured and summed them to create an Index, as a measure of kinship.

RESULTS

The kinship was determined by the physical length of shared chromosomal regions that are distinguished by each kinship. Using this method, the relationship was able be distinguished up to the fourth degree of kinship, and non-relatives were clearly distinguished from true relatives. We also found a potential for this approach to be used universally, regardless of microarray platforms for SNP genotyping and populations.

CONCLUSION

This method has a potential to determine the different degree of kinship between individuals and to distinguish non-relatives from true relatives, which can be of great help for practical applications in kinship determination.

Epigenetic variation impacts individual differences in the transcriptional response to influenza infection

Humans display remarkable interindividual variation in their immune response to identical challenges. Yet, our understanding of the genetic and epigenetic factors contributing to such variation remains limited. Here we performed in-depth genetic, epigenetic and transcriptional profiling on primary macrophages derived from individuals of European and African ancestry before and after infection with influenza A virus. We show that baseline epigenetic profiles are strongly predictive of the transcriptional response to influenza A virus across individuals. Quantitative trait locus (QTL) mapping revealed highly coordinated genetic effects on gene regulation, with many cis-acting genetic variants impacting concomitantly gene expression and multiple epigenetic marks. These data reveal that ancestry-associated differences in the epigenetic landscape can be genetically controlled, even more than gene expression. Lastly, among QTL variants that colocalized with immune-disease loci, only 7% were gene expression QTL, while the remaining genetic variants impact epigenetic marks, stressing the importance of considering molecular phenotypes beyond gene expression in disease-focused studies.

Recombination, admixture and genome instability shape the genomic landscape of Saccharomyces cerevisiae derived from spontaneous grape ferments

Cultural exchange of fermentation techniques has driven the spread of Saccharomyces cerevisiae across the globe, establishing natural populations in many countries. Despite this, Oceania is thought to lack native populations of S. cerevisiae, only being introduced after colonisation. Here we investigate the genomic landscape of 411 S. cerevisiae isolated from spontaneous grape fermentations in Australia across multiple locations, years, and grape cultivars. Spontaneous fermentations contained highly recombined mosaic strains that exhibited high levels of genome instability. Assigning genomic windows to putative ancestral origin revealed that few closely related starter lineages have come to dominate the genetic landscape, contributing most of the genetic variation. Fine-scale phylogenetic analysis of loci not observed in strains of commercial wine origin identified widespread admixture with European derived beer yeast along with three independent admixture events from potentially endemic Oceanic lineages that was associated with genome instability. Finally, we investigated Australian ecological niches for basal isolates, identifying phylogenetically distinct S. cerevisiae of non-European, non-domesticated origin associated with admixture loci. Our results illustrate the effect commercial use of microbes may have on local microorganism genetic diversity and demonstrates the presence of non-domesticated, potentially endemic lineages of S. cerevisiae in Australian niches that are actively admixing.

Machine Learning Prediction of Treatment Response to Inhaled Corticosteroids in Asthma

BACKGROUND

Although inhaled corticosteroids (ICS) are the first-line therapy for patients with persistent asthma, many patients continue to have exacerbations. We developed machine learning models to predict the ICS response in patients with asthma.

METHODS

The subjects included asthma patients of European ancestry (n = 1371; 448 children; 916 adults). A genome-wide association study was performed to identify the SNPs associated with ICS response. Using the SNPs identified, two machine learning models were developed to predict ICS response: (1) least absolute shrinkage and selection operator (LASSO) regression and (2) random forest.

RESULTS

The LASSO regression model achieved an AUC of 0.71 (95% CI 0.67-0.76; sensitivity: 0.57; specificity: 0.75) in an independent test cohort, and the random forest model achieved an AUC of 0.74 (95% CI 0.70-0.78; sensitivity: 0.70; specificity: 0.68). The genes contributing to the prediction of ICS response included those associated with ICS responses in asthma (TPSAB1, FBXL16), asthma symptoms and severity (ABCA7, CNN2, PTRN3, and BSG/CD147), airway remodeling (ELANE, FSTL3), mucin production (GAL3ST), leukotriene synthesis (GPX4), allergic asthma (ZFPM1, SBNO2), and others.

CONCLUSIONS

An accurate risk prediction of ICS response can be obtained using machine learning methods, with the potential to inform personalized treatment decisions. Further studies are needed to examine if the integration of richer phenotype data could improve risk prediction.

Regulation of kynurenine metabolism by blood donor genetics and biology impacts red cell hemolysis in vitro and in vivo

ABSTRACT

In the field of transfusion medicine, the clinical relevance of the metabolic markers of the red blood cell (RBC) storage lesion is incompletely understood. Here, we performed metabolomics of RBC units from 643 donors enrolled in the Recipient Epidemiology and Donor Evaluation Study, REDS RBC Omics. These units were tested on storage days 10, 23, and 42 for a total of 1929 samples and also characterized for end-of-storage hemolytic propensity after oxidative and osmotic insults. Our results indicate that the metabolic markers of the storage lesion poorly correlated with hemolytic propensity. In contrast, kynurenine was not affected by storage duration and was identified as the top predictor of osmotic fragility. RBC kynurenine levels were affected by donor age and body mass index and were reproducible within the same donor across multiple donations from 2 to 12 months apart. To delve into the genetic underpinnings of kynurenine levels in stored RBCs, we thus tested kynurenine levels in stored RBCs on day 42 from 13 091 donors from the REDS RBC Omics study, a population that was also genotyped for 879 000 single nucleotide polymorphisms. Through a metabolite quantitative trait loci analysis, we identified polymorphisms in SLC7A5, ATXN2, and a series of rate-limiting enzymes (eg, kynurenine monooxygenase, indoleamine 2,3-dioxygenase, and tryptophan dioxygenase) in the kynurenine pathway as critical factors affecting RBC kynurenine levels. By interrogating a donor-recipient linkage vein-to-vein database, we then report that SLC7A5 polymorphisms are also associated with changes in hemoglobin and bilirubin levels, suggestive of in vivo hemolysis in 4470 individuals who were critically ill and receiving single-unit transfusions.

Examining interactions between polygenic scores and interpersonal trauma exposure on alcohol consumption and use disorder in an ancestrally diverse college cohort

Introduction: Genetic factors impact alcohol consumption and use disorder (AUD), with large-scale genome-wide association studies (GWAS) identifying numerous associated variants. Aggregate genetic methods in combination with important environmental factors (e.g., interpersonal trauma [IPT]) can be applied to expand our understanding of the ways by which genetic and environmental variables work together to influence alcohol consumption and disordered use. The present study aimed to detail the relationships between genome-wide polygenic scores (PGS) for alcohol phenotypes (i.e., alcohol consumption and AUD status) and IPT exposure as well as the interaction between them across ancestry. Methods: Data were drawn from the Spit for Science (S4S) study, a US college student population, where participants reported on IPT exposure prior to college and alcohol consumption and problems during college (N = 9,006; ancestry: 21.3% African [AFR], 12.5% Admixed Americas [AMR], 9.6% East Asian [EAS], 48.1% European [EUR], 8.6% South Asian [SAS]). Two trans-ancestry PGS were constructed, one for alcohol consumption and another for AUD, using large-scale GWAS summary statistics from multiple ancestries weighted using PRS-CSx. Regression models were applied to test for the presence of associations between alcohol-PGS and IPT main and interaction effects. Results: In the meta-analysis across ancestry groups, IPT exposure and PGS were significantly associated with alcohol consumption (βIPT = 0.31, P IPT = 0.0002; βPGS = 0.09, P PGS = 0.004) and AUD (ORIPT = 1.12, P IPT = 3.5 × 10-8; ORPGS = 1.02, P PGS = 0.002). No statistically significant interactions were detected between IPT and sex nor between IPT and PGS. When inspecting ancestry specific results, the alcohol consumption-PGS and AUD-PGS were only statistically significant in the EUR ancestry group (βPGS = 0.09, P PGS = 0.04; ORPGS = 1.02, P PGS = 0.022, respectively). Discussion: IPT exposure prior to college was strongly associated with alcohol outcomes in this college-age sample, which could be used as a preventative measure to identify students at high risk for problematic alcohol use. Additionally, results add to developing evidence of polygenic score association in meta-analyzed samples, highlighting the importance of continued efforts to increase ancestral representation in genetic studies and inclusive analytic approaches to increase the generalizability of results from genetic association studies.

Molecular analysis and transfusion management in a rare case of cis-AB blood group: A report from India

Molecular characterization of a rare cis-AB blood group has not been done in the Indian subcontinent. Herein, we report a case of A2B3 blood group in an Indian patient which was subsequently confirmed to be a case of cis-AB phenotype. Blood grouping was performed by the column agglutination technique (CAT), conventional tube technique (CTT) and subsequently, whole exome sequencing for molecular analysis. The patient was initially typed as AB, RhD positive in forward grouping. However, serum grouping showed agglutination (2+) with the B red cells in CAT. In CTT, an extra reaction was observed with A1 red cells and a strong agglutination was seen with Anti-H lectin. Thus, the blood group was identified serologically as A2B3. During the next-generation sequencing, a total of 10 exonic variants in the ABO gene were filtered, of which 2 (rs8176747 and rs7853989) were found to be non-synonymous and occurring on the same allele. The other allele was found to be ABO*A1.01. The sample analyzed in the study was found to carry two previously reported nucleotide changes of cis-AB (c.803G > C and c.526C > G) on the same allele which had not been reported before. Transfusion requirement was managed with type O red cells and type AB plasma.

Association study of human leukocyte antigen variants and idiopathic pulmonary fibrosis

Introduction

Idiopathic pulmonary fibrosis (IPF) is a chronic interstitial pneumonia marked by progressive lung fibrosis and a poor prognosis. Recent studies have highlighted the potential role of infection in the pathogenesis of IPF, and a prior association of the HLA-DQB1 gene with idiopathic fibrotic interstitial pneumonia (including IPF) has been reported. Owing to the important role that the human leukocyte antigen (HLA) region plays in the immune response, here we evaluated if HLA genetic variation was associated specifically with IPF risk.

Methods

We performed a meta-analysis of associations of the HLA region with IPF risk in individuals of European ancestry from seven independent case-control studies of IPF (comprising 5159 cases and 27 459 controls, including a prior study of fibrotic interstitial pneumonia). Single nucleotide polymorphisms, classical HLA alleles and amino acids were analysed and signals meeting a region-wide association threshold of p<4.5×10-4 and a posterior probability of replication >90% were considered significant. We sought to replicate the previously reported HLA-DQB1 association in the subset of studies independent of the original report.

Results

The meta-analysis of all seven studies identified four significant independent single nucleotide polymorphisms associated with IPF risk. However, none met the posterior probability for replication criterion. The HLA-DQB1 association was not replicated in the independent IPF studies.

Conclusion

Variation in the HLA region was not consistently associated with risk in studies of IPF. However, this does not preclude the possibility that other genomic regions linked to the immune response may be involved in the aetiology of IPF.

Genome evolution is surprisingly predictable after initial hybridization

Over the past two decades, evolutionary biologists have come to appreciate that hybridization, or genetic exchange between distinct lineages, is remarkably common - not just in particular lineages but in taxonomic groups across the tree of life. As a result, the genomes of many modern species harbor regions inherited from related species. This observation has raised fundamental questions about the degree to which the genomic outcomes of hybridization are repeatable and the degree to which natural selection drives such repeatability. However, a lack of appropriate systems to answer these questions has limited empirical progress in this area. Here, we leverage independently formed hybrid populations between the swordtail fish Xiphophorus birchmanni and X. cortezi to address this fundamental question. We find that local ancestry in one hybrid population is remarkably predictive of local ancestry in another, demographically independent hybrid population. Applying newly developed methods, we can attribute much of this repeatability to strong selection in the earliest generations after initial hybridization. We complement these analyses with time-series data that demonstrates that ancestry at regions under selection has remained stable over the past ~40 generations of evolution. Finally, we compare our results to the well-studied X. birchmanni×X. malinche hybrid populations and conclude that deeper evolutionary divergence has resulted in stronger selection and higher repeatability in patterns of local ancestry in hybrids between X. birchmanni and X. cortezi.

Epigenetic aging differences between Wichí and Criollos from Argentina: Insights from genomic history and ecology

Background and objectives

Epigenetic estimators based on DNA methylation levels have emerged as promising biomarkers of human aging. These estimators exhibit natural variations across human groups, but data about indigenous populations remain underrepresented in research. This study aims to investigate differences in epigenetic estimators between two distinct human populations, both residing in the Gran Chaco region of Argentina, the Native-American Wichí, and admixed Criollos who are descendants of intermarriages between Native Americans and the first European colonizers, using a population genetic approach.

Methodology

We analyzed 24 Wichí (mean age: 39.2 ± 12.9 yo) and 24 Criollos (mean age: 41.1 ± 14.0 yo) for DNA methylation levels using the Infinium MethylationEPIC (Illumina) to calculate 16 epigenetic estimators. Additionally, we examined genome-wide genetic variation using the HumanOmniExpress BeadChip (Illumina) to gain insights into the genetic history of these populations.

Results

Our results indicate that Native-American Wichí are epigenetically older compared to Criollos according to five epigenetic estimators. Analyses within the Criollos population reveal that global ancestry does not influence the differences observed, while local (chromosomal) ancestry shows positive associations between specific SNPs located in genomic regions over-represented by Native-American ancestry and measures of epigenetic age acceleration (AgeAccelHannum). Furthermore, we demonstrate that differences in population ecologies also contribute to observed epigenetic differences.

Conclusions and implications

Overall, our study suggests that while the genomic history may partially account for the observed epigenetic differences, non-genetic factors, such as lifestyle and ecological factors, play a substantial role in the variability of epigenetic estimators, thereby contributing to variations in human epigenetic aging.

The effects of SARS-CoV-2 on susceptible human cells

The biological consequences of viral infection result from biochemical, physiological, structural, morphological and genetic changes in infected cells. In productive infections, virus-induced biological changes in cells may be closely related to the efficiency of viral replication or to the recognition of these cells by the immune system. These changes are usually associated with cytocidal viruses, as in the case of the pandemic coronavirus SARS-CoV-2, which causes COVID-19. Many of these changes are required for effective viral replication. The physiological state of living cells has a significant impact on the outcome of viral infection, as the host cell provides the synthetic machinery, key regulatory molecules and precursors for newly synthesised viral proteins and nucleic acids. This review focuses on novel target cell types for SARS-CoV-2 exposure outside the respiratory tract. Findings and examples are collected that provide information on virus-cell interactions. The identification of unusual target cells for SARS-CoV-2 may help to explain the diverse symptoms in COVID-19 patients and the long-lasting effects after infection. In particular, the discovery of previously undescribed target cells for SARS-CoV-2 action needs to be considered to improve treatment of patients and prevention of infection.

Association study of human leukocyte antigen (HLA) variants and idiopathic pulmonary fibrosis

Introduction

Idiopathic pulmonary fibrosis (IPF) is a chronic interstitial pneumonia marked by progressive lung fibrosis and a poor prognosis. Recent studies have highlighted the potential role of infection in the pathogenesis of IPF and a prior association of the HLA-DQB1 gene with idiopathic fibrotic interstitial pneumonia (including IPF) has been reported. Due to the important role that the Human Leukocyte Antigen (HLA) region plays in the immune response, here we evaluated if HLA genetic variation was associated specifically with IPF risk.

Methods

We performed a meta-analysis of associations of the HLA region with IPF risk in individuals of European ancestry from seven independent case-control studies of IPF (comprising a total of 5,159 cases and 27,459 controls, including the prior study of fibrotic interstitial pneumonia). Single nucleotide polymorphisms, classical HLA alleles and amino acids were analysed and signals meeting a region-wide association threshold p<4.5×10-4 and a posterior probability of replication >90% were considered significant. We sought to replicate the previously reported HLA-DQB1 association in the subset of studies independent of the original report.

Results

The meta-analysis of all seven studies identified four significant independent single nucleotide polymorphisms associated with IPF risk. However, none met the posterior probability for replication criterion. The HLA-DQB1 association was not replicated in the independent IPF studies.

Conclusion

Variation in the HLA region was not consistently associated with risk in studies of IPF. However, this does not preclude the possibility that other genomic regions linked to the immune response may be involved in the aetiology of IPF.

High-Altitude Andean H194R HIF2A Allele Is a Hypomorphic Allele

For over 10,000 years, Andeans have resided at high altitude where the partial pressure of oxygen challenges human survival. Recent studies have provided evidence for positive selection acting in Andeans on the HIF2A (also known as EPAS1) locus, which encodes for a central transcription factor of the hypoxia-inducible factor pathway. However, the precise mechanism by which this allele might lead to altitude-adaptive phenotypes, if any, is unknown. By analyzing whole genome sequencing data from 46 high-coverage Peruvian Andean genomes, we confirm evidence for positive selection acting on HIF2A and a unique pattern of variation surrounding the Andean-specific single nucleotide variant (SNV), rs570553380, which encodes for an H194R amino acid substitution in HIF-2α. Genotyping the Andean-associated SNV rs570553380 in a group of 299 Peruvian Andeans from Cerro de Pasco, Peru (4,338 m), reveals a positive association with increased fraction of exhaled nitric oxide, a marker of nitric oxide biosynthesis. In vitro assays show that the H194R mutation impairs binding of HIF-2α to its heterodimeric partner, aryl hydrocarbon receptor nuclear translocator. A knockin mouse model bearing the H194R mutation in the Hif2a gene displays decreased levels of hypoxia-induced pulmonary Endothelin-1 transcripts and protection against hypoxia-induced pulmonary hypertension. We conclude the Andean H194R HIF2A allele is a hypomorphic (partial loss of function) allele.

A biologically informed polygenic score of neuronal plasticity moderates the association between cognitive aptitudes and cortical thickness in adolescents

Although many studies of the adolescent brain identified positive associations between cognitive abilities and cortical thickness, little is known about mechanisms underlying such brain-behavior relationships. With experience-induced plasticity playing an important role in shaping the cerebral cortex throughout life, it is likely that some of the inter-individual variations in cortical thickness could be explained by genetic variations in relevant molecular processes, as indexed by a polygenic score of neuronal plasticity (PGS-NP). Here, we studied associations between PGS-NP, cognitive abilities, and thickness of the cerebral cortex, estimated from magnetic resonance images, in the Saguenay Youth Study (SYS, 533 females, 496 males: age=15.0 ± 1.8 years of age; cross-sectional), and the IMAGEN Study (566 females, 556 males; between 14 and 19 years; longitudinal). Using Gene Ontology, we first identified 199 genes implicated in neuronal plasticity, which mapped to 155,600 single nucleotide polymorphisms (SNPs). Second, we estimated their effect sizes from an educational attainment meta-GWAS to build a PGS-NP. Third, we examined a possible moderating role of PGS-NP in the relationship between performance intelligence quotient (PIQ), and its subtests, and the thickness of 34 cortical regions. In SYS, we observed a significant interaction between PGS-NP and object assembly vis-à-vis thickness in male adolescents (p = 0.026). A median-split analysis showed that, in males with a 'high' PGS-NP, stronger associations between object assembly and thickness were found in regions with larger age-related changes in thickness (r = 0.55, p = 0.00075). Although the interaction between PIQ and PGS-NP was non-significant (p = 0.064), we performed a similar median-split analysis. Again, in the high PGS-NP males, positive associations between PIQ and thickness were observed in regions with larger age-related changes in thickness (r = 0.40, p = 0.018). In the IMAGEN cohort, we did not replicate the first set of results (interaction between PGS-NP and cognitive abilities via-a-vis cortical thickness) while we did observe the same relationship between the brain-behaviour relationship and (longitudinal) changes in cortical thickness (Matrix reasoning: r = 0.63, p = 6.5e-05). No statistically significant results were observed in female adolescents in either cohort. Overall, these cross-sectional and longitudinal results suggest that molecular mechanisms involved in neuronal plasticity may contribute to inter-individual variations of cortical thickness related to cognitive abilities during adolescence in a sex-specific manner.

PharmVar GeneFocus: SLCO1B1

The Pharmacogene Variation Consortium (PharmVar) is now providing star (*) allele nomenclature for the highly polymorphic human SLCO1B1 gene encoding the organic anion transporting polypeptide 1B1 (OATP1B1) drug transporter. Genetic variation within the SLCO1B1 gene locus impacts drug transport, which can lead to altered pharmacokinetic profiles of several commonly prescribed drugs. Variable OATP1B1 function is of particular importance regarding hepatic uptake of statins and the risk of statin-associated musculoskeletal symptoms. To introduce this important drug transporter gene into the PharmVar database and serve as a unified reference of haplotype variation moving forward, an international group of gene experts has performed an extensive review of all published SLCO1B1 star alleles. Previously published star alleles were self-assigned by authors and only loosely followed the star nomenclature system that was first developed for cytochrome P450 genes. This nomenclature system has been standardized by PharmVar and is now applied to other important pharmacogenes such as SLCO1B1. In addition, data from the 1000 Genomes Project and investigator-submitted data were utilized to confirm existing haplotypes, fill knowledge gaps, and/or define novel star alleles. The PharmVar-developed SLCO1B1 nomenclature has been incorporated by the Clinical Pharmacogenetics Implementation Consortium (CPIC) 2022 guideline on statin-associated musculoskeletal symptoms.

Twelve Loci Associated With Bone Density in Middle-aged and Elderly Chinese: The Shanghai Changfeng Study

CONTEXT

Previous genome-wide association studies (GWASs) of bone mineral density (BMD) were mainly conducted in Europeans.

OBJECTIVE

To explore genetic variants that affect BMD and sex differences in a Chinese population.

METHODS

A total of 5428 middle-aged and elderly Chinese were included. Dual-energy X-ray absorptiometry was used to measure BMD at the lumbar spine, and total and specific sites of the hip. A mixed linear model was used to analyze the associations between BMD and autosomal genetic variants, adjusting for age, age squared, sex, and menopausal women (model 1); model 2 was further adjusted for height and weight. A GWAS of osteoporosis in the Biobank Japan (BBJ) project was used for replication. GWAMA software was used to detect the statistical significance of sex differences of estimated effects. Gene annotation and pathway enrichment analysis were performed.

RESULTS

Women lost BMD at earlier ages and faster than men. The 2 models identified a total of 12 loci that were associated with BMD at any site. Single nucleotide polymorphisms rs72354346, rs2024219, rs1463093, rs10037512, and rs5880932 were successfully replicated in the BBJ. Variations of rs79262027 G>A (VKORC1L1) and rs4795209 A>G (DDX52) were associated with BMD only in men, and rs1239055408 G>GA (KCNJ2) was associated with BMD only in women. Gene enrichment analysis showed that BMD in a Chinese elderly population was mainly related to ossification, bone resorption, sex hormones, and kidney physiology.

CONCLUSION

The present GWAS identified 12 loci that were significantly associated with BMD at any site in a Chinese population, and 3 of them showed sex differences in their effects.

Natural selection of immune and metabolic genes associated with health in two lowland Bolivian populations

A growing body of work has addressed human adaptations to diverse environments using genomic data, but few studies have connected putatively selected alleles to phenotypes, much less among underrepresented populations such as Amerindians. Studies of natural selection and genotype-phenotype relationships in underrepresented populations hold potential to uncover previously undescribed loci underlying evolutionarily and biomedically relevant traits. Here, we worked with the Tsimane and the Moseten, two Amerindian populations inhabiting the Bolivian lowlands. We focused most intensively on the Tsimane, because long-term anthropological work with this group has shown that they have a high burden of both macro and microparasites, as well as minimal cardiometabolic disease or dementia. We therefore generated genome-wide genotype data for Tsimane individuals to study natural selection, and paired this with blood mRNA-seq as well as cardiometabolic and immune biomarker data generated from a larger sample that included both populations. In the Tsimane, we identified 21 regions that are candidates for selective sweeps, as well as 5 immune traits that show evidence for polygenic selection (e.g., C-reactive protein levels and the response to coronaviruses). Genes overlapping candidate regions were strongly enriched for known involvement in immune-related traits, such as abundance of lymphocytes and eosinophils. Importantly, we were also able to draw on extensive phenotype information for the Tsimane and Moseten and link five regions (containing PSD4, MUC21 and MUC22, TOX2, ANXA6, and ABCA1) with biomarkers of immune and metabolic function. Together, our work highlights the utility of pairing evolutionary analyses with anthropological and biomedical data to gain insight into the genetic basis of health-related traits.

Polygenic liability, stressful life events and risk for secondary-treated depression in early life: a nationwide register-based case-cohort study

BACKGROUND

In this study, we examined the relationship between polygenic liability for depression and number of stressful life events (SLEs) as risk factors for early-onset depression treated in inpatient, outpatient or emergency room settings at psychiatric hospitals in Denmark.

METHODS

Data were drawn from the iPSYCH2012 case-cohort sample, a population-based sample of individuals born in Denmark between 1981 and 2005. The sample included 18 532 individuals who were diagnosed with depression by a psychiatrist by age 31 years, and a comparison group of 20 184 individuals. Information on SLEs was obtained from nationwide registers and operationalized as a time-varying count variable. Hazard ratios and cumulative incidence rates were estimated using Cox regressions.

RESULTS

Risk for depression increased by 35% with each standard deviation increase in polygenic liability (p < 0.0001), and 36% (p < 0.0001) with each additional SLE. There was a small interaction between polygenic liability and SLEs (β = -0.04, p = 0.0009). The probability of being diagnosed with depression in a hospital-based setting between ages 15 and 31 years ranged from 1.5% among males in the lowest quartile of polygenic liability with 0 events by age 15, to 18.8% among females in the highest quartile of polygenic liability with 4+ events by age 15.

CONCLUSIONS

These findings suggest that although there is minimal interaction between polygenic liability and SLEs as risk factors for hospital-treated depression, combining information on these two important risk factors could potentially be useful for identifying high-risk individuals.

A Pipeline for Phasing and Genotype Imputation on Mixed Human Data (Parents-Offspring Trios and Unrelated Subjects) by Reviewing Current Methods and Software

Genotype imputation has become an essential prerequisite when performing association analysis. It is a computational technique that allows us to infer genetic markers that have not been directly genotyped, thereby increasing statistical power in subsequent association studies, which consequently has a crucial impact on the identification of causal variants. Many features need to be considered when choosing the proper algorithm for imputation, including the target sample on which it is performed, i.e., related individuals, unrelated individuals, or both. Problems could arise when dealing with a target sample made up of mixed data, composed of both related and unrelated individuals, especially since the scientific literature on this topic is not sufficiently clear. To shed light on this issue, we examined existing algorithms and software for performing phasing and imputation on mixed human data from SNP arrays, specifically when related subjects belong to trios. By discussing the advantages and limitations of the current algorithms, we identified LD-based methods as being the most suitable for reconstruction of haplotypes in this specific context, and we proposed a feasible pipeline that can be used for imputing genotypes in both phased and unphased human data.

Genome-wide metabolite quantitative trait loci analysis (mQTL) in red blood cells from volunteer blood donors

The red blood cell (RBC)-Omics study, part of the larger NHLBI-funded Recipient Epidemiology and Donor Evaluation Study (REDS-III), aims to understand the genetic contribution to blood donor RBC characteristics. Previous work identified donor demographic, behavioral, genetic, and metabolic underpinnings to blood donation, storage, and (to a lesser extent) transfusion outcomes, but none have yet linked the genetic and metabolic bodies of work. We performed a genome-wide association (GWA) analysis using RBC-Omics study participants with generated untargeted metabolomics data to identify metabolite quantitative trait loci in RBCs. We performed GWA analyses of 382 metabolites in 243 individuals imputed using the 1000 Genomes Project phase 3 all-ancestry reference panel. Analyses were conducted using ProbABEL and adjusted for sex, age, donation center, number of whole blood donations in the past 2 years, and first 10 principal components of ancestry. Our results identified 423 independent genetic loci associated with 132 metabolites (p < 5×10^-8). Potentially novel locus-metabolite associations were identified for the region encoding heme transporter FLVCR1 and choline and for lysophosphatidylcholine acetyltransferase LPCAT3 and lysophosphatidylserine 16.0, 18.0, 18.1, and 18.2; these associations are supported by published rare disease and mouse studies. We also confirmed previous metabolite GWA results for associations, including N(6)-methyl-L-lysine and protein PYROXD2 and various carnitines and transporter SLC22A16. Association between pyruvate levels and G6PD polymorphisms was validated in an independent cohort and novel murine models of G6PD deficiency (African and Mediterranean variants). We demonstrate that it is possible to perform metabolomics-scale GWA analyses with a modest, trans-ancestry sample size.

Uganda Genome Resource: A rich research database for genomic studies of communicable and non-communicable diseases in Africa

The Uganda Genome Resource (UGR) is a well-characterized genomic database with a range of phenotypic communicable and non-communicable diseases and risk factors generated from the Uganda General Population Cohort (GPC), a population-based open cohort established in 1989. The UGR comprises genotype data on ∼5,000 and whole-genome sequence data on ∼2,000 Ugandan GPC individuals from 10 ethno-linguistic groups. Leveraging other platforms at MRC/UVRI and LSHTM Uganda Research Unit, there is opportunity for additional sample collection to expand the UGR to advance scientific discoveries. Here, we describe UGR and highlight how it is providing opportunities for discovery of novel disease susceptibility genetic loci, refining association signals at new and existing loci, developing and testing polygenic scores to determine disease risk, assessing causal relations in diseases, and developing capacity for genomics research in Africa. The UGR has the potential to develop to a comparable level of European and Asian large-scale genomic initiatives.

Maternal-fetal genetic interactions, imprinting, and risk of placental abruption

RESULTS

Abruption cases were more likely to experience preeclampsia, have shorter gestational age, and deliver infants with lower birthweight compared with controls. Models with MFGI effects provided improved fit than models with only maternal and fetal genotype main effects for SNP rs12530904 (p-value = 1.2e-04) in calcium/calmodulin-dependent protein kinase [CaM kinase] II beta (CAMK2B), and, SNP rs73136795 (p-value = 1.9e-04) in peroxisome proliferator-activated receptor-gamma (PPARG), both MB genes. We identified 320 SNPs in 45 maternally-imprinted genes (including potassium voltage-gated channel subfamily Q member 1 [KCNQ1], neurotrimin [NTM], and, ATPase phospholipid transporting 10 A [ATP10A]) associated with abruption. Top hits included rs2012323 (p-value = 1.6E-16) and rs12221520 (p-value1.3e-13) in KCNQ1, rs8036892 (p-value = 9.3E-17) and rs188497582 in ATP10A, rs12589854 (p-value = 2.9E-11) and rs80203467 (p-value = 4.6e-11) in maternally expressed 8, small nucleolar RNA host (MEG8), and rs138281088 in solute carrier family 22 member 2 (SLC22A2) (p-value = 6.8e-9).

CONCLUSIONS

We identified novel PA-related maternal-fetal MB gene interactions and imprinting effects that highlight the role of the fetus in PA risk development. Findings can inform mechanistic investigations to understand the pathogenesis of PA.

Identification of Novel Loci Shared by Juvenile Idiopathic Arthritis Subtypes Through Integrative Genetic Analysis

OBJECTIVE

Juvenile idiopathic arthritis (JIA) is the most common chronic immune-mediated joint disease among children and encompasses a heterogeneous group of immune-mediated joint disorders classified into 7 subtypes according to clinical presentation. However, phenotype overlap and biologic evidence suggest a shared mechanistic basis between subtypes. This study was undertaken to systematically investigate shared genetic underpinnings of JIA subtypes.

METHODS

We performed a heterogeneity-sensitive genome-wide association study encompassing a total of 1,245 JIA cases (classified into 7 subtypes) and 9,250 controls, followed by fine-mapping of candidate causal variants at each genome-wide significant locus, functional annotation, and pathway and network analysis. We further identified candidate drug targets and drug repurposing opportunities by in silico analyses.

RESULTS

In addition to the major histocompatibility complex locus, we identified 15 genome-wide significant loci shared between at least 2 JIA subtypes, including 10 novel loci. Functional annotation indicated that candidate genes at these loci were expressed in diverse immune cell types.

CONCLUSION

This study identified novel genetic loci shared by JIA subtypes. Our findings identified candidate mechanisms underlying JIA subtypes and candidate targets with drug repurposing opportunities for JIA treatment.

Association of Genome-Wide Polygenic Risk Score for Body Mass Index With Cardiometabolic Health From Childhood Through Midlife

BACKGROUND

Genetic information may help to identify individuals in childhood who are at increased risk for cardiometabolic disease.

METHODS

We included 1201 BHS (Bogalusa Heart Study) participants (832 White participants and 369 Black participants) who were followed up to 42.3 years, starting at a mean age of 9.8 years. A validated genome-wide polygenic risk score (PRS) was tested for association with midlife body mass index (BMI), fasting plasma glucose, and systolic blood pressure using multiple linear regression models. Cox proportional hazards models tested associations of the PRS with incident obesity, diabetes, and hypertension. All analyses were conducted according to race and adjusted for baseline age, sex, ancestry, and BMI.

RESULTS

The constructed PRS was significantly and modestly correlated with midlife BMI in both White and Black participants, with correlation coefficients of 0.27 (P=1.94×10-8) and 0.16 (P=5.50×10-3), respectively. In White participants, per SD increase of PRS was associated with an average 1.29 kg/m2 higher BMI (P=4.44×10-9), 2.82 mg/dL higher fasting plasma glucose (P=1.17×10-3), and 1.09 mm Hg higher systolic blood pressure (P=3.57×10-2) at midlife. The PRS also conferred a 26% higher increased risk of obesity (P=3.50×10-6) in White participants. In addition, the variance in midlife BMI explained increased from 0.1973 to 0.2293 when PRS was added to the model including age, sex, principal components, and baseline BMI (P<0.0001). No associations were observed in Black participants.

CONCLUSIONS

Adiposity-related genetic information independently predicted cardiometabolic health in White BHS participants. Null associations observed in Black BHS participants highlight the urgent need for PRS development in multi-ancestry populations.

Ancient Maltese genomes and the genetic geography of Neolithic Europe

Archaeological consideration of maritime connectivity has ranged from a biogeographical perspective that considers the sea as a barrier to a view of seaways as ancient highways that facilitate exchange. Our results illustrate the former. We report three Late Neolithic human genomes from the Mediterranean island of Malta that are markedly enriched for runs of homozygosity, indicating inbreeding in their ancestry and an effective population size of only hundreds, a striking illustration of maritime isolation in this agricultural society. In the Late Neolithic, communities across mainland Europe experienced a resurgence of hunter-gatherer ancestry, pointing toward the persistence of different ancestral strands that subsequently admixed. This is absent in the Maltese genomes, giving a further indication of their genomic insularity. Imputation of genome-wide genotypes in our new and 258 published ancient individuals allowed shared identity-by-descent segment analysis, giving a fine-grained genetic geography of Neolithic Europe. This highlights the differentiating effects of seafaring Mediterranean expansion and also island colonization, including that of Ireland, Britain, and Orkney. These maritime effects contrast profoundly with a lack of migratory barriers in the establishment of Central European farming populations from Anatolia and the Balkans.

Association between Breast Cancer Polygenic Risk Score and Chemotherapy-Induced Febrile Neutropenia: Null Results

BACKGROUND

The hypothesis that breast cancer (BC) susceptibility variants are linked to chemotherapy-induced toxicity has been previously explored. Here, we investigated the association between a validated 313-marker-based BC polygenic risk score (PRS) and chemotherapy-induced neutropenia without fever and febrile neutropenia (FNc) in Asian BC patients.

METHODS

This observational case-control study of Asian BC patients treated with chemotherapy included 161 FNc patients, 219 neutropenia patients, and 936 patients who did not develop neutropenia. A continuous PRS was calculated by summing weighted risk alleles associated with overall, estrogen receptor- (ER-) positive, and ER-negative BC risk. PRS distributions neutropenia or FNc cases were compared to controls who did not develop neutropenia using two-sample t-tests. Odds ratios (OR) and corresponding 95% confidence intervals were estimated for the associations between PRS (quartiles and per standard deviation (SD) increase) and neutropenia-related outcomes compared to controls.

RESULTS

PRS distributions were not significantly different in any of the comparisons. Higher PRS_overall quartiles were negatively correlated with neutropenia or FNc. However, the associations were not statistically significant (PRS per SD increase OR neutropenia: 0.91 [0.79-1.06]; FNc: 0.87 [0.73-1.03]). No dose-dependent trend was observed for the ER-positive weighted PRS (PRS_ER-pos) and ER-negative weighted PRS (PRS_ER-neg).

CONCLUSION

BC PRS was not strongly associated with chemotherapy-induced neutropenia or FNc.

HaploMaker: An improved algorithm for rapid haplotype assembly of genomic sequences

BACKGROUND

In diploid organisms, whole-genome haplotype assembly relies on the accurate identification and assignment of heterozygous single-nucleotide polymorphism alleles to the correct homologous chromosomes. This appropriate phasing of these alleles ensures that combinations of single-nucleotide polymorphisms on any chromosome, called haplotypes, can then be used in downstream genetic analysis approaches including determining their potential association with important phenotypic traits. A number of statistical algorithms and complementary computational software tools have been developed for whole-genome haplotype construction from genomic sequence data. However, many algorithms lack the ability to phase long haplotype blocks and simultaneously achieve a competitive accuracy.

RESULTS

In this research we present HaploMaker, a novel reference-based haplotype assembly algorithm capable of accurately and efficiently phasing long haplotypes using paired-end short reads and longer Pacific Biosciences reads from diploid genomic sequences. To achieve this we frame the problem as a directed acyclic graph with edges weighted on read evidence and use efficient path traversal and minimization techniques to optimally phase haplotypes. We compared the HaploMaker algorithm with 3 other common reference-based haplotype assembly tools using public haplotype data of human individuals from the Platinum Genome project. With short-read sequences, the HaploMaker algorithm maintained a competitively low switch error rate across all haplotype lengths and was superior in phasing longer genomic regions. For longer Pacific Biosciences reads, the phasing accuracy of HaploMaker remained competitive for all block lengths and generated substantially longer block lengths than the competing algorithms.

CONCLUSIONS

HaploMaker provides an improved haplotype assembly algorithm for diploid genomic sequences by accurately phasing longer haplotypes. The computationally efficient and portable nature of the Java implementation of the algorithm will ensure that it has maximal impact in reference-sequence-based haplotype assembly applications.

Development of New Inhibitors of HDAC1-3 Enzymes Aided by In Silico Design Strategies

Histone deacetylases (HDACs) are overexpressed in cancer, and their inhibition shows promising results in cancer therapy. In particular, selective class I HDAC inhibitors such as entinostat are proposed to be more beneficial in breast cancer treatment. Computational drug design is an inevitable part of today's drug discovery projects because of its unequivocal role in saving time and cost. Using three HDAC inhibitors trichostatin, vorinostat, and entinostat as template structures and a diverse fragment library, all synthetically accessible compounds thereof (∼3200) were generated virtually and filtered based on similarity against the templates and PAINS removal. The 298 selected structures were docked into the active site of HDAC I and ranked using a calculated binding affinity. Top-ranking structures were inspected manually, and, considering the ease of synthesis and drug-likeness, two new structures (3a and 3b) were proposed for synthesis and biological evaluation. The synthesized compounds were purified to a degree of more than 95% and structurally verified using various methods. The designed compounds 3a and 3b showed 65-80 and 5% inhibition on HDAC 1, 2, and 3 isoforms at a concentration of 10 μM, respectively. The novel compound 3a may be used as a lead structure for designing new HDAC inhibitors.