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.

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.

Combining developmental and sleep health measures for autism spectrum disorder screening: an ECHO study.

BACKGROUND: Sleep problems are reported for up to 80% of autistic individuals. We examined whether parsimonious sets of items derived from the Modified Checklist for Autism in Toddlers, Revised (M-CHAT-R) and the Brief Infant Sleep Questionnaire (BISQ) are superior to the standard M-CHAT-R in predicting subsequent autism spectrum disorder (ASD) diagnoses. METHODS: Participants from 11 Environmental influences on Child Health Outcomes (ECHO) cohorts were included. We performed logistic LASSO regression models with 10-fold cross-validation to identify whether a combination of items derived from the M-CHAT-R and BISQ are superior to the standard M-CHAT-R in predicting ASD diagnoses. RESULTS: The final sample comprised 1552 children. The standard M-CHAT-R had a sensitivity of 44% (95% CI: 34, 55), specificity of 92% (95% CI: 91, 94), and AUROC of 0.726 (95% CI: 0.663, 0.790). A higher proportion of children with ASD had difficulty falling asleep or resisted bedtime during infancy/toddlerhood. However, LASSO models revealed parental reports of sleep problems did not improve the accuracy of the M-CHAT-R in predicting ASD diagnosis. CONCLUSION: While children with ASD had higher rates of sleep problems during infancy/toddlerhood, there was no improvement in ASD developmental screening through the incorporation of parent-report sleep metrics. IMPACT: Parental-reported sleep problems are common in autism spectrum disorder (ASD). We investigated whether the inclusion of parental-reports of infant/toddler sleep patterns enhanced the effectiveness of developmental screening for autism. We reported higher rates of difficulty falling asleep and resisting bedtime during infancy and toddlerhood among children later diagnosed with ASD; however, we did not find an improvement in ASD developmental screening through the incorporation of parent-report sleep metrics. In our sample, the standard M-CHAT-R had a sensitivity of 39% among children of mothers with government insurance compared with a sensitivity of 53% among children of mothers with employer-based insurance.

Smoking-informed methylation and expression QTLs in human brain and colocalization with smoking-associated genetic loci.

Smoking is a leading cause of preventable morbidity and mortality. Smoking is heritable, and genome-wide association studies (GWASs) of smoking behaviors have identified hundreds of significant loci. Most GWAS-identified variants are noncoding with unknown neurobiological effects. We used genome-wide genotype, DNA methylation, and RNA sequencing data in postmortem human nucleus accumbens (NAc) to identify cis-methylation/expression quantitative trait loci (meQTLs/eQTLs), investigate variant-by-cigarette smoking interactions across the genome, and overlay QTL evidence at smoking GWAS-identified loci to evaluate their regulatory potential. Active smokers (N = 52) and nonsmokers (N = 171) were defined based on cotinine biomarker levels and next-of-kin reporting. We simultaneously tested variant and variant-by-smoking interaction effects on methylation and expression, separately, adjusting for biological and technical covariates and correcting for multiple testing using a two-stage procedure. We found >2 million significant meQTL variants (padj < 0.05) corresponding to 41,695 unique CpGs. Results were largely driven by main effects, and five meQTLs, mapping to NUDT12, FAM53B, RNF39, and ADRA1B, showed a significant interaction with smoking. We found 57,683 significant eQTL variants for 958 unique eGenes (padj < 0.05) and no smoking interactions. Colocalization analyses identified loci with smoking-associated GWAS variants that overlapped meQTLs/eQTLs, suggesting that these heritable factors may influence smoking behaviors through functional effects on methylation/expression. One locus containing MUSTN1 and ITIH4 colocalized across all data types (GWAS, meQTL, and eQTL). In this first genome-wide meQTL map in the human NAc, the enriched overlap with smoking GWAS-identified genetic loci provides evidence that gene regulation in the brain helps explain the neurobiology of smoking behaviors.

Comparison of maternal venous blood metabolomics collected as dried blood spots, dried blood microsamplers, and plasma for integrative environmental health research.

Use of capillary blood devices for exposome research can deepen our understanding of the intricate relationship between environment and health, and open up new avenues for preventive and personalized medicine, particularly for vulnerable populations. While the potential of these whole blood devices to accurately measure chemicals and metabolites has been demonstrated, how untargeted metabolomics data from these samplers can be integrated with previous and ongoing environmental health studies that have used conventional blood collection approaches is not yet clear. Therefore, we performed a comprehensive comparison between relative-quantitative metabolite profiles measured in venous blood collected with dried whole blood microsamplers (DBM), dried whole blood spots (DBS), and plasma from 54 mothers in an ethnically diverse population. We determined that a majority of the 309 chemicals and metabolites showed similar median intensity rank, moderate correlation, and moderate agreement between participant-quantiled intraclass correlation coefficients (ICCs) for pair-wise comparisons among the three biomatrices. In particular, whole blood sample types, DBM and DBS, were in highest agreement across metabolite comparison metrics, followed by metabolites measured in DBM and plasma, and then metabolites measured in DBS and plasma. We provide descriptive characteristics and measurement summaries as a reference database. This includes unique metabolites that were particularly concordant or discordant in pairwise comparisons. Our results demonstrate that the range of metabolites from untargeted metabolomics data collected with DBM, DBS, and plasma provides biologically relevant information for use in independent exposome investigations. However, before meta-analysis with combined datasets are performed, robust statistical approaches that integrate untargeted metabolomics data collected on different blood matrices need to be developed.

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.

Functional effects of an African glucose-6-phosphate dehydrogenase (G6PD) polymorphism (Val68Met) on red blood cell hemolytic propensity and post-transfusion recovery.

BACKGROUND: Donor genetic variation is associated with red blood cell (RBC) storage integrity and post-transfusion recovery. Our previous large-scale genome-wide association study demonstrated that the African G6PD deficient A- variant (rs1050828, Val68Met) is associated with higher oxidative hemolysis after cold storage. Despite a high prevalence of X-linked G6PD mutation in African American population (>10%), blood donors are not routinely screened for G6PD status and its importance in transfusion medicine is relatively understudied. STUDY DESIGN AND METHODS: To further evaluate the functional effects of the G6PD A- mutation, we created a novel mouse model carrying this genetic variant using CRISPR-Cas9. We hypothesize that this humanized G6PD A- variant is associated with reduced G6PD activity with a consequent effect on RBC hemolytic propensity and post-transfusion recovery. RESULTS: G6PD A- RBCs had reduced G6PD protein with ~5% residual enzymatic activity. Significantly increased in vitro hemolysis induced by oxidative stressors was observed in fresh and stored G6PD A- RBCs, along with a lower GSH:GSSG ratio. However, no differences were observed in storage hemolysis, osmotic fragility, mechanical fragility, reticulocytes, and post-transfusion recovery. Interestingly, a 14% reduction of 24-h survival following irradiation was observed in G6PD A- RBCs compared to WT RBCs. Metabolomic assessment of stored G6PD A- RBCs revealed an impaired pentose phosphate pathway (PPP) with increased glycolytic flux, decreasing cellular antioxidant capacity. DISCUSSION: This novel mouse model of the common G6PD A- variant has impaired antioxidant capacity like humans and low G6PD activity may reduce survival of transfused RBCs when irradiation is performed.

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. 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.

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.

Genetic polymorphisms and expression of Rhesus blood group RHCE are associated with 2,3-bisphosphoglycerate in humans at high altitude.

Red blood cell (RBC) metabolic reprogramming upon exposure to high altitude contributes to physiological human adaptations to hypoxia, a multifaceted process critical to health and disease. To delve into the molecular underpinnings of this phenomenon, first, we performed a multi-omics analysis of RBCs from six lowlanders after exposure to high-altitude hypoxia, with longitudinal sampling at baseline, upon ascent to 5,100 m and descent to sea level. Results highlighted an association between erythrocyte levels of 2,3-bisphosphoglycerate (BPG), an allosteric regulator of hemoglobin that favors oxygen off-loading in the face of hypoxia, and expression levels of the Rhesus blood group RHCE protein. We then expanded on these findings by measuring BPG in RBCs from 13,091 blood donors from the Recipient Epidemiology and Donor Evaluation Study. These data informed a genome-wide association study using BPG levels as a quantitative trait, which identified genetic polymorphisms in the region coding for the Rhesus blood group RHCE as critical determinants of BPG levels in erythrocytes from healthy human volunteers. Mechanistically, we suggest that the Rh group complex, which participates in the exchange of ammonium with the extracellular compartment, may contribute to intracellular alkalinization, thus favoring BPG mutase activity.

Evaluating 17 methods incorporating biological function with GWAS summary statistics to accelerate discovery demonstrates a tradeoff between high sensitivity and high positive predictive value.

Where sufficiently large genome-wide association study (GWAS) samples are not currently available or feasible, methods that leverage increasing knowledge of the biological function of variants may illuminate discoveries without increasing sample size. We comprehensively evaluated 17 functional weighting methods for identifying novel associations. We assessed the performance of these methods using published results from multiple GWAS waves across each of five complex traits. Although no method achieved both high sensitivity and positive predictive value (PPV) for any trait, a subset of methods utilizing pleiotropy and expression quantitative trait loci nominated variants with high PPV (>75%) for multiple traits. Application of functionally weighting methods to enhance GWAS power for locus discovery is unlikely to circumvent the need for larger sample sizes in truly underpowered GWAS, but these results suggest that applying functional weighting to GWAS can accurately nominate additional novel loci from available samples for follow-up studies.

Smoking-informed methylation and expression QTLs in human brain and colocalization with smoking-associated genetic loci.

Smoking is a leading cause of preventable morbidity and mortality. Smoking is heritable, and genome-wide association studies (GWAS) of smoking behaviors have identified hundreds of significant loci. Most GWAS-identified variants are noncoding with unknown neurobiological effects. We used genome-wide genotype, DNA methylation, and RNA sequencing data in postmortem human nucleus accumbens (NAc) to identify cis-methylation/expression quantitative trait loci (meQTLs/eQTLs), investigate variant-by-cigarette smoking interactions across the genome, and overlay QTL evidence at smoking GWAS-identified loci to evaluate their regulatory potential. Active smokers (N=52) and nonsmokers (N=171) were defined based on cotinine biomarker levels and next-of-kin reporting. We simultaneously tested variant and variant-by-smoking interaction effects on methylation and expression, separately, adjusting for biological and technical covariates and using a two-stage multiple testing approach with eigenMT and Bonferroni corrections. We found >2 million significant meQTL variants (padj<0.05) corresponding to 41,695 unique CpGs. Results were largely driven by main effects; five meQTLs, mapping to NUDT12, FAM53B, RNF39, and ADRA1B, showed a significant interaction with smoking. We found 57,683 significant eQTLs for 958 unique eGenes (padj<0.05) and no smoking interactions. Colocalization analyses identified loci with smoking-associated GWAS variants that overlapped meQTLs/eQTLs, suggesting that these heritable factors may influence smoking behaviors through functional effects on methylation/expression. One locus containing MUSTIN1 and ITIH4 colocalized across all data types (GWAS + meQTL + eQTL). In this first genome-wide meQTL map in the human NAc, the enriched overlap with smoking GWAS-identified genetic loci provides evidence that gene regulation in the brain helps explain the neurobiology of smoking behaviors.

Evaluation of pediatric epigenetic clocks across multiple tissues.

BACKGROUND: Epigenetic clocks are promising tools for assessing biological age. We assessed the accuracy of pediatric epigenetic clocks in gestational and chronological age determination. RESULTS: Our study used data from seven tissue types on three DNA methylation profiling microarrays and found that the Knight and Bohlin clocks performed similarly for blood cells, while the Lee clock was superior for placental samples. The pediatric-buccal-epigenetic clock performed the best for pediatric buccal samples, while the Horvath clock is recommended for children's blood cell samples. The NeoAge clock stands out for its unique ability to predict post-menstrual age with high correlation with the observed age in infant buccal cell samples. CONCLUSIONS: Our findings provide valuable guidance for future research and development of epigenetic clocks in pediatric samples, enabling more accurate assessments of biological age.

Genetic Predisposition to Adverse Neurodevelopmental Outcome of Extremely Low Birth Weight Infants.

OBJECTIVE: This study aimed to evaluate whether there are genetic variants associated with adverse neurodevelopmental outcomes in extremely low birth weight (ELBW) infants. STUDY DESIGN: We conducted a candidate gene association study in two well-defined cohorts of ELBW infants (<1,000 g). One cohort was for discovery and the other for replication. The discovery case-control analysis utilized anonymized DNA samples and evaluated 1,614 single-nucleotide polymorphisms (SNPs) in 145 genes concentrated in inflammation, angiogenesis, brain development, and oxidation pathways. Cases were children who died by age one or who were diagnosed with cerebral palsy (CP) or neurodevelopmental delay (Bayley II mental developmental index [MDI] or psychomotor developmental index [PDI] < 70) by 18 to 22 months. Controls were survivors with normal neurodevelopment. We assessed significant epidemiological variables and SNPs associated with the combined outcome of CP or death, CP, mental delay (MDI < 70) and motor delay (PDI < 70). Multivariable analyses adjusted for gestational age at birth, small for gestational age, sex, antenatal corticosteroids, multiple gestation, racial admixture, and multiple comparisons. SNPs associated with adverse neurodevelopmental outcomes with p < 0.01 were selected for validation in the replication cohort. Successful replication was defined as p < 0.05 in the replication cohort. RESULTS: Of 1,013 infants analyzed (452 cases, 561 controls) in the discovery cohort, 917 were successfully genotyped for >90% of SNPs and passed quality metrics. After adjusting for covariates, 26 SNPs with p < 0.01 for one or more outcomes were selected for replication cohort validation, which included 362 infants (170 cases and 192 controls). A variant in SERPINE1, which encodes plasminogen activator inhibitor (PAI1), was associated with the combined outcome of CP or death in the discovery analysis (p = 4.1 × 10-4) and was significantly associated with CP or death in the replication cohort (adjusted odd ratio: 0.4; 95% confidence interval: 0.2-1.0; p = 0.039). CONCLUSION: A genetic variant in SERPINE1, involved in inflammation and coagulation, is associated with CP or death among ELBW infants. KEY POINTS: · Early preterm and ELBW infants have dramatically increased risks of CP and developmental delay.. · A genetic variant in SERPINE1 is associated with CP or death among ELBW infants.. · The SERPINE1 gene encodes the serine protease inhibitor plasminogen activator inhibitor..

The Environmental Influences on Child Health Outcomes (ECHO)-Wide Cohort.

The Environmental Influences on Child Health Outcomes (ECHO)-Wide Cohort Study (EWC), a collaborative research design comprising 69 cohorts in 31 consortia, was funded by the National Institutes of Health (NIH) in 2016 to improve children's health in the United States. The EWC harmonizes extant data and collects new data using a standardized protocol, the ECHO-Wide Cohort Data Collection Protocol (EWCP). EWCP visits occur at least once per life stage, but the frequency and timing of the visits vary across cohorts. As of March 4, 2022, the EWC cohorts contributed data from 60,553 children and consented 29,622 children for new EWCP data and biospecimen collection. The median (interquartile range) age of EWCP-enrolled children was 7.5 years (3.7-11.1). Surveys, interviews, standardized examinations, laboratory analyses, and medical record abstraction are used to obtain information in 5 main outcome areas: pre-, peri-, and postnatal outcomes; neurodevelopment; obesity; airways; and positive health. Exposures include factors at the level of place (e.g., air pollution, neighborhood socioeconomic status), family (e.g., parental mental health), and individuals (e.g., diet, genomics).

Genome-wide association study of early ischaemic stroke risk in Brazilian individuals with sickle cell disease implicates ADAMTS2 and CDK18 and uncovers novel loci.

Ischaemic stroke is a common complication of sickle cell disease (SCD) and without intervention can affect 11% of children with SCD before the age of 20. Within the Trans-Omics for Precision Medicine (TOPMed), a genome-wide association study (GWAS) of ischaemic stroke was performed on 1333 individuals with SCD from Brazil (178 cases, 1155 controls). Via a novel Cox proportional-hazards analysis, we searched for variants associated with ischaemic stroke occurring at younger ages. Variants at genome-wide significance (p < 5 × 10-8 ) include two near genes previously linked to non-SCD early-onset stroke (<65 years): ADAMTS2 (rs147625068, p = 3.70 × 10-9 ) and CDK18 (rs12144136, p = 2.38 × 10-9 ). Meta-analysis, which included the independent SCD cohorts Walk-PHaSST and PUSH, exhibited consistent association for variants rs1209987 near gene TBC1D32 (p = 3.36 × 10-10 ), rs188599171 near CUX1 (p = 5.89 × 10-11 ), rs77900855 near BTG1 (p = 4.66 × 10-8 ), and rs141674494 near VPS13C (1.68 × 10-9 ). Findings from this study support a multivariant model of early ischaemic stroke risk and possibly a shared genetic architecture between SCD individuals and non-SCD individuals younger than 65 years.

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.

Erythrocyte mitogen-activated protein kinases mediate hemolytic events under osmotic and oxidative stress and in hemolytic diseases.

p38 MAPKs are key regulators of cellular adaptation to various stress stimuli, however, their role in mediating erythrocyte cell death and hemolysis is largely unknown. We hypothesized that activation of erythrocyte p38 MAPK is a common event in the stimulation of hemolysis, and that inhibition of p38 MAPK pathways could mitigate hemolysis in hemoglobinopathies. We exposed human erythrocytes to diamide-induced oxidative stress or to hypoosmotic shock in the presence or absence of p38 MAPK inhibitors (SCIO469, SB203580, CMPD1) and used immunoblotting to determine MAPK activity and to identify possible downstream effectors of p38 MAPK. We also evaluated the impact of p38 MAPK inhibitors on stress-induced hemolysis or hypoxia-induced sickling in erythrocytes from mouse models of sickle cell disease. We found that human erythrocytes express conventional MAPKs (MKK3, p38 MAPK, MAPKAPK2) and identified differential MAPK activation pathways in each stress condition. Specifically, p38 MAPK inhibition in diamide-treated erythrocytes was associated with decreased phosphorylation of Src tyrosine kinases and Band 3 protein. Conversely, hypoosmotic shock induced MAPKAPK2 and RSK2 phosphorylation, which was inhibited by SCIO469 or CMPD1. Relevant to hemoglobinopathies, sickle cell disease was associated with increased erythrocyte MKK3, p38 MAPK, and MAPKAPK2 expression and phosphorylation as compared with erythrocytes from healthy individuals. Furthermore, p38 MAPK inhibition was associated with decreased hemolysis in response to diamide treatments or osmotic shock, and with decreased erythrocyte sickling under experimental hypoxia. These findings provided insights into MAPK-mediated signaling pathways that regulate erythrocyte function and hemolysis in response to extracellular stressors or human diseases.

GAWMerge expands GWAS sample size and diversity by combining array-based genotyping and whole-genome sequencing.

Genome-wide association studies (GWAS) have made impactful discoveries for complex diseases, often by amassing very large sample sizes. Yet, GWAS of many diseases remain underpowered, especially for non-European ancestries. One cost-effective approach to increase sample size is to combine existing cohorts, which may have limited sample size or be case-only, with public controls, but this approach is limited by the need for a large overlap in variants across genotyping arrays and the scarcity of non-European controls. We developed and validated a protocol, Genotyping Array-WGS Merge (GAWMerge), for combining genotypes from arrays and whole-genome sequencing, ensuring complete variant overlap, and allowing for diverse samples like Trans-Omics for Precision Medicine to be used. Our protocol involves phasing, imputation, and filtering. We illustrated its ability to control technology driven artifacts and type-I error, as well as recover known disease-associated signals across technologies, independent datasets, and ancestries in smoking-related cohorts. GAWMerge enables genetic studies to leverage existing cohorts to validly increase sample size and enhance discovery for understudied traits and ancestries.