Circulating metabolites, genetics and way of life components in relation to future threat of kind 2 diabetes

This web page was created programmatically, to learn the article in its authentic location you may go to the hyperlink bellow:
https://www.nature.com/articles/s41591-025-04105-8
and if you wish to take away this text from our web site please contact us


Study contributors and ethics approval

Our MWAS for incident T2D includes the usage of knowledge from ten potential cohorts, together with the Nurses’ Health Study (NHS; initiated in 1976 with 121,701 feminine nurses aged 30–55 years9,50), NHS2 (began in 1989 with 116,429 feminine nurses aged 25–42 years9,50), Health Professionals Follow-Up Study (HPFS; began in 1986 with 51,529 male well being professions aged 40–75 years9), Hispanic Community Health Study/Study of Latinos (SOL; enrolled 16,415 Hispanic/Latino adults aged 18–74 years throughout 2008–201151,52), Women’s Health Initiative (WHI; initiated in 1993 enrolling 68,132 girls aged 50–79 years to considered one of three scientific trials or an observational research53), Atherosclerosis Risk in Communities (ARIC) research (enrolled 15,792 largely Black and white US adults aged 45–64 years throughout 1987–198954), Framingham Heart Study Offspring cohort (FHS; enrolled 5,124 adults; we targeted on these attended the fifth examination throughout 1991–1995), Multi-Ethnic Study of Atherosclerosis (MESA; initiated in 2000 with 6,814 adults aged 45–84 years55,56), the Boston Puerto Rican Health Study (BPRHS; enrolled 1,500 self-identified Puerto Rican adults aged 45–75 years) and the Prevención con Dieta Mediterránea Study (PREDIMED; a 5-year dietary trial with 7,447 adults aged 55–80 years57). In every cohort, complete knowledge on demographics, medical and household historical past, food regimen, way of life and different well being info have been collected at baseline and have been up to date throughout longitudinal follow-ups. Blood samples have been collected at baseline and/or throughout follow-ups. Our MWAS for incident T2D included contributors with certified metabolomics knowledge, and have been freed from diabetes, heart problems and most cancers at research baseline. The remaining evaluation included 6,890 contributors from NHS; 3,692 from NHS2 and a couple of,529 from HPFS; 2,821 from SOL; 1,392 from WHI; 1,288 white and 1,433 Black contributors from ARIC; 1,424 from FHS; 902 from MESA; 378 from BPRHS and 885 from PREDIMED (Extended Data Table 1). Each research was permitted by Institutional Review Boards at respective establishments or research facilities, and all contributors offered knowledgeable consent. Our GWAS for metabolites included contributors from eight cohorts comprising NHS, NHS2, HPFS, SOL, WHI, ARIC, FHS and, as well as, the Cardiovascular Health Study (CHS; enrolled 5,201 adults throughout 1989–1990 and 678 predominantly Black contributors in 1992–199358,59) (Supplementary Table 7). The detailed descriptions of the design, knowledge assortment, moral overview of every cohort, and our inclusion and exclusion standards are offered in Supplementary Methods.

Ascertainment of T2D

In all cohorts, incident T2D was outlined when a participant was freed from diabetes at baseline however was recognized as having T2D throughout longitudinal follow-up. Detailed info on analysis standards in every cohort is included in Supplementary Methods, and follow-up years and numbers of incident circumstances are listed in Extended Data Table 1. Briefly, in NHS/HPFS, T2D have been recognized by follow-up questionnaires, and confirmed by way of a supplementary questionnaire primarily based on diagnostic standards from the National Diabetes Data Group earlier than 199860 and the American Diabetes Association (ADA) standards after 199861,62. In SOL, T2D was outlined if a participant had fasting glucose ≥7.0 mmol l−1, fasting ≤8 h and nonfasting glucose ≥11.1 mmol l−1, put up oral glucose tolerance check glucose ≥11.1 mmol l−1, HbA1c ≥ 6.5%, present use of antidiabetic medicines or self-reported physician-diagnosed diabetes63. In WHI, T2D was decided primarily based on self-reported historical past of diabetes or utilizing antidiabetic medicines (capsules or photographs) in any visits/interviews. In ARIC and FHS, T2D was recognized if an individual had fasting glucose ≥7.0 mmol l−1, fasting ≤8 h and nonfasting glucose ≥11.1 mmol l−1, or present use of antidiabetic medicines with ARIC additional contemplating self-reported physician-diagnosed diabetes64,65. T2D circumstances in MESA and BPRHS have been decided in accordance with the ADA standards66, which included fasting plasma glucose degree ≥7.0 mmol l−1 or the usage of antidiabetic medicines or insulin56,67. In PREDIMED, T2D was adjudicated by way of blind evaluation by a Clinical Endpoint and Adjudication of Events Committee, primarily based on the ADA standards68.

Assessment of food regimen, way of life components and covariates

Detailed info on knowledge assortment in every cohort is in Supplementary Methods. Briefly, demographic components (for instance, self-reported intercourse, and race and ethnicity), socioeconomic standing, well being info (for instance, medical circumstances and household historical past) and way of life (for instance, smoking historical past and PAs), anthropometrics and blood strain, have been collected at baseline and follow-up visits, by way of self-administrated questionnaires, or in-person or telephone-based interviews by skilled employees. PA was quantified as metabolic equal (MET) in hours per week. We calculated BMI primarily based on baseline weight and peak, and WHR primarily based on waist and hip circumferences. Blood scientific biomarkers have been measured utilizing commonplace assays. Among contributors with serum creatinine knowledge, eGFR was estimated utilizing the CKD-EPI (Chronic Kidney Disease Epidemiology Collaboration) formulation, primarily based on age, intercourse and race in NHS/HPFS, WHI, ARIC and PREDIMED69, and commonplace reference equations for Hispanics adjusting for age and intercourse in SOL. In PREDIMED, two propensity scores have been estimated to account for the likelihood of project to intervention teams57.

In NHS/HPFS, food regimen was assessed utilizing a semi-quantitative meals frequency questionnaire (FFQ) each 4 years; in our evaluation we averaged the intakes from the 2 FFQs closest to the time of blood draw (NHS: 1986 and 1990; NHS2: 1995 and 1999; HPFS: 1994 and 1998). In WHI, ARIC, FHS, MESA and BPRHS, food regimen was equally assessed by FFQs designed and validated for software to their focused populations (for instance, multiethnic and geographically various populations in WHI70,71,72 and Puerto Rican inhabitants in BPRHS73). In SOL, food regimen was assessed utilizing two 24-h dietary recollects and a meals propensity questionnaire74. The general dietary high quality was assessed by the Alternate Healthy Eating Index-2010 (AHEI-2010)75 in all cohorts aside from the PREDIMED trial, through which it was assessed by a 14-item Mediterranean Diet Adherence Screener rating57. In NHS/HPFS, SOL and WHI, we additionally calculated baseline consumptions of 15 most important meals teams within the unit of servings per day.

Metabolomic profiling, high quality management and knowledge harmonization

Metabolomic profiling in NHS/HPFS, WHI, MESA, PREDIMED, FHS and CHS was carried out with the Metabolomics Platforms on the Broad Institute of MIT and Harvard University, utilizing three to 4 complementary LC–MS strategies9,65,76. Metabolomic profiling in SOL and ARIC (serum samples) and BPRHS (plasma samples) was carried out utilizing LC–MS primarily based strategies by the Metabolon DiscoveryHD4 Panel on the Metabolon Inc.63,77,78. Detailed protocols for each platforms have been described beforehand53,79.

Data processing was carried out inside every research and, if relevant, individually inside every batch (or substudy) if a number of batches/substudies have been carried out inside a cohort. Samples have been eliminated if their metabolite detection charge was <80%, or have been recognized as outliers by multidimensional scaling evaluation inside a particular race/ethnic group. Metabolites have been filtered if their detection charge throughout samples was <80% and, if relevant, had a coefficient of variation >20% for high quality management (QC) samples. After high quality filtering, missingness of every metabolite have been imputed utilizing the half minimal worth, and the info have been then standardized for evaluation. Across all cohorts, we matched metabolites by their HMDB ID and/or PubChem ID, offered by the corresponding metabolomic laboratories. A complete of 1,273 named metabolites have been initially certified for evaluation in at the very least one cohort. To cut back single-study bias, we restricted our analyses to 469 metabolites that have been out there in at the very least 4 unbiased cohorts, or out there in at the very least three unbiased cohorts if the three cohorts coated each Metabolomic platforms. Finally, 407 metabolites from NHS, 363 from NHS2, 291 from HPFS, 364 from WHI, 327 from MESA, 274 from PREDIMED, 188 from FHS, 283 from SOL, 139 from ARIC and 231 from BPRHS have been harmonized for our evaluation (Extended Data Table 1). In CHS, 411 metabolites have been included in genetic analyses (Supplementary Table 7). Details of the metabolomic profiling, QC and knowledge processing are within the Supplementary Methods.

Metabolome-wide affiliation evaluation for incident T2D

Details of analytical approaches and fashions are offered in Supplementary Methods and Supplementary Table 1. Briefly, all affiliation analyses have been carried out individually for every cohort, stratified by main racial/ethnic teams when pattern sizes permitted. Metabolites have been inversely regular remodeled by every substudy and racial/ethnic group (if relevant) in every cohort. To analyze the affiliation between every metabolite and T2D threat, we utilized Cox regression for research of longitudinal cohort design (NHS excluding the T2D nested case–management substudy, NHS2, HPFS, SOL, ARIC, WHI, FHS, MESA and BPRHS); logistic regression for the NHS T2D nested case–management substudy; and Cox regression with Barlow weights80 and strong estimators for the PREDIMED T2D nested case–cohort research. The primary multivariate mannequin (mannequin 1) was adjusted for age, intercourse, smoking standing, alcohol consumption and, if relevant, schooling, household revenue, fasting standing, lipid-lowering medicines, anti-hypertensive medicines, household historical past of diabetes, self-reported physician-diagnosed hypertension, self-reported physician-diagnosed dyslipidemia and study-specific covariates. The most important mannequin was additional adjusted for BMI and WHR (mannequin 2). In sensitivity analyses, mannequin 1 was additional adjusted for PA and dietary high quality index (mannequin 3); high-density lipoprotein (HDL)-cholesterol, low-density lipoprotein (LDL)-cholesterol and triglycerides (mannequin 4), or systolic and diastolic blood pressures (mannequin 5). In one other sensitivity evaluation, mannequin 2 was additional adjusted for eGFR in NHS, NHS2, HPFS, SOL, ARIC, WHI and PREDIMED. For every metabolite, affiliation outcomes from all out there cohorts and racial/ethnic teams have been mixed utilizing a fixed-effect, inverse-variance-weighted (IVW) meta-analysis, and a meta-analyzed FDR < 0.05 was thought of statistically vital. In secondary analyses, meta-analysis was carried out combining outcomes from the identical racial/ethnic teams, or cohorts utilizing the identical platforms.

To annotate the novelty of the recognized associations, we reviewed earlier potential cohort research linking circulating metabolites to T2D threat. We used a literature-review-based meta-analysis4 that included all research printed earlier than 6 March 2021 as an anchor, and searched for added research printed from 2021 to 202421,81,82,83,84,85,86,87,88,89,90,91,92,93,94. We thought of an affiliation as ‘previously reported,’ if the affiliation was statistically vital in a printed research after a number of testing correction primarily based on the research’s prespecified evaluation plan.

GWAS of metabolites

Detailed info on genotyping arrays, imputation strategies, pattern dimension and GWAS and meta-analysis strategies, is offered in Supplementary Methods and Supplementary Table 7. Briefly, genotyping have been carried out utilizing a number of sorts of array by earlier research in NHS/HPFS95, SOL96, ARIC7, WHI97, CHS98 and FHS43. Imputation was carried out primarily based on the HRC reference panel in NHS/HPFS and CHS; 1000 Genomes Project section 3 worldwide reference panel in SOL, 1000 Genomes Project section 3 v.5 in WHI and HapMap CEU inhabitants launch v.22 in FHS with complete pre- and postimputation QC. GWAS of metabolites have been carried out beforehand within the NHS/HPFS (median n = 6,610, vary 971–8,054) and WHI (n = 1,256) utilizing the RVTESTS device6,42,99, in SOL (n = 3,933) utilizing a linear mixed-effect mannequin in GMMAT7 and in ARIC (n = 1,772 and n = 1509 for African American and non-Hispanic white contributors, respectively)7, CHS (n = 263) and FHS (n = 1,802)43, with detailed evaluation procedures described in earlier publications7,42,43.

GWAS abstract statistics from every cohort have been lifted over to Genome Build v.37 and filtered, retaining single nucleotide polymorphisms with a minor allele frequency ≥ 0.01 and imputation ratio ≥0.3. For every metabolite, an IVW fixed-effect meta-analysis, applied in METAL100, was used to mix GWAS outcomes from the cohorts through which the metabolite was out there. Genomic management was applied earlier than and after meta-analysis100. The remaining GWAS have been out there for 458 out of 469 harmonized metabolites, with the entire pattern dimension starting from 1,074 to 18,590 (median n = 8,611). We in contrast vital mQTLs recognized at P < 5 × 10−8 and 1.09 × 10−10 (that’s, 5 × 10−8 additional correcting for 458 metabolites) ranges. Manhattan plots have been derived utilizing R package deal CMplot and regional plots have been draw with LocusZoom101. In a secondary evaluation, we in contrast genetic impact heterogeneity between racial/ethnic teams on the recognized mQTLs for T2D-associated metabolites (Supplementary Methods).

We annotate the novelty of our vital mQTLs for the 165 T2D-associated metabolites at P < 1.09 × 10−10, by evaluating our outcomes to eight earlier research (with N ≥ 4,000 and used LC–MS primarily based metabolomic platforms)8,24,25,26,102,103,104,105. We thought of a locus for a particular metabolite as ‘previously reported’ if the reported lead genetic variant was the identical lead variant, or not the identical lead variant however was vital in our research; or not in our research however throughout the clumping vary of our recognized locus. We thought of a locus for a metabolite as doubtlessly new if our locus was not beforehand reported for this metabolite, or this metabolite was not beforehand reported in these research.

Lead variants for metabolites, pathway evaluation and proportion of variance defined

We used the PLINK clumping operate (P < 5 × 10−8 and r2 < 0.01 in a 1,000-kb window) to establish unbiased genetic variants related to every metabolite. For metabolite with no variant at P < 5 × 10−8, a single lead variant with the smallest P was chosen. Gene annotation for prime variants was carried out utilizing the SNPNexus net device106. Canonical pathway enrichment analyses was carried out utilizing the MetaCore software program with the default background107; and we in contrast prime enriched pathways for genes annotated to mQTLs of T2D-related metabolites versus these of non-associated metabolites. We calculated the R2 of every metabolite defined by unbiased lead genetic variants utilizing the formulation ({sum }_{i=1}^{okay}beta instances beta instances 2times {rm{MAF}}instances (1-{rm{MAF}})), through which okay is the variety of unbiased lead variants, and β is the affiliation coefficient between the variant and the metabolite. We in contrast the R2 distribution for the T2D-associated versus non-associated metabolites utilizing Wilcoxon check.

Genetic correlation r
g between metabolites and T2D-related traits

We acquired publicly out there GWAS abstract statistics from giant consortium research for T2D (180,834 circumstances and 1,159,055 controls)27, fasting insulin (N = 98,210)108, proinsulin (N = 45,861)109, HOMA-IR and HOMA-B (N = 51,750)110, BMI-adjusted insulin sensitivity index (ISI, N = 53,657) and insulin fold-change (IFC; N = 55,124)111, BMI and WHR (N = 700,000)112 and lipids (N = 1,500,000)113. We carried out GWAS for HBA1c (N = 390,982), subcutaneous fats quantity (N = 37,912), visceral fats quantity (N = 37,912), liver proton density fats fraction (PDFF; N = 29,512), pancreas PDFF (N = 28,624) and liver enzymes (N = 390,000) within the UK Biobank utilizing BOLT-LMM (Supplementary Methods). We calculated rg between every metabolite and every scientific trait utilizing linkage disequilibrium rating regression, primarily based on their GWAS abstract knowledge overlapping with the 1.2 M HapMap3 variants after excluding the main histocompatibility complicated area within the European inhabitants114. For every trait, we in contrast the distribution of its rg with T2D-associated versus non-associated metabolites, utilizing chi-squared check, and regarded FDR < 0.05 (correcting for numbers of comparisons examined) as statistically vital.

Genetic colocalization

We obtained tissue-specific cis-eQTLs abstract statistics from the GTEx undertaking v.8115,116. The shared causal variants between every metabolite and tissue-specific transcriptome from 47 tissue sorts, have been examined utilizing colocalization evaluation applied within the coloc.abf() operate in R package deal ‘coloc’ v.5117. For every metabolite, we enter the GWAS abstract statistics for all variants inside ±500 kb of its unbiased lead variants (Supplementary Methods). A posterior likelihood of H4 (PPH4) > 0.8 was thought of as robust proof for genetic colocalization. Within every tissue kind, we used univariant logistic regression to check whether or not the proportions of mQTL–eQTL colocalizations are increased for the T2D-associated versus non-associated metabolites, and a one-sided FDR < 0.05 (correcting for 47 tissue sorts) was thought of as statistically vital. We utilized an identical coloc strategy to look at genetic colocalizations between circulating metabolites and T2D27. We then aligned mQTL–T2D colocalizations with tissue-specific eQTL–mQTL colocalizations by metabolites and shared causal variants, to interpret the potential performance of metabolites in T2D pathogenesis.

MR evaluation

To infer the potential causal relationships between 233 T2D-associated metabolites (with genetic knowledge) and T2D threat, we utilized 4 MR strategies applied within the MendelianRandomization R package deal118: we used mode-based estimate (MBE) as the primary methodology as it’s usually conservative and strong to outliers; we additional utilized weighted-median, IVW and MR-egger to point end result consistency119. When testing the course from metabolites to T2D, we used unbiased variants from clumping (P < 5 × 10−8 and r2 < 0.01 in a 1,000-kb window) excluding the HLA area as genetic instrumental variables. If fewer than three variants have been recognized, we lowered the clumping P threshold till at the very least three variants have been recognized. We thought of a possible causal relationship when MBE–FDR < 0.05 and at the very least two different MR strategies confirmed the identical impact instructions as these from MBE. Sensitivity analyses have been carried out, both to take away variants mapped to the highest 3 recurrent loci (GCKR, ZNF259, FADS cluster) from the instrumental variables, or to make use of solely unbiased variants clumped at P < 1.09 × 10−10 because the instrumental variables of metabolites, utilizing the IVW MR methodology (on account of fewer variants retained). When testing the course from T2D to metabolites, we used unbiased lead variants related to T2D at P < 5 × 10−8 because the instrumental variables. For the 148 metabolites which are potential mediators between BMI and T2D threat, we utilized MR evaluation to check the course from BMI to metabolites. Details are offered in Supplementary Methods.

MWASs for modifiable threat components

We fitted linear fashions to regress inversely regular remodeled metabolite ranges on age, intercourse (solely in SOL), present smoking standing, BMI, PA, intakes of 15 most important meals teams and fasting standing, concurrently along with cohort-specific covariates. Analyses have been carried out in NHS/HPFS, SOL and WHI, individually, additional stratified by substudies or racial teams (Supplementary Methods). Association coefficients between metabolites and every specific threat issue have been then mixed throughout analytical units utilizing a fixed-effect IVW meta-analysis. The R2 of every metabolite defined by particular threat components have been first calculated in every analytical set utilizing the formulation (beta instances beta instances {mathrm{variance}}left({mathrm{threat}}; {mathrm{issue}}proper)!{/mathrm{variance}}left({mathrm{metabolite}}proper)), with the β being the affiliation coefficients between the metabolite and the chance issue; after which averaged throughout all analytical units. We in contrast the distributions of R2 for T2D-associated versus non-associated metabolites utilizing the Wilcoxon check.

Mediation evaluation between threat components, metabolites and T2D threat

Details for mediation evaluation are described in Supplementary Methods. Briefly, our evaluation targeted on BMI, PA, espresso/tea consumption and purple/processed meat consumption. For every threat issue, metabolites (1) that have been related to each the chance issue and T2D threat and (2) whose affiliation instructions with the chance issue and T2D threat have been in keeping with the pre-assumed epidemiological relationships between the chance issue and T2D threat, have been thought of. We examined whether or not, and to what diploma, every metabolite mediated the affiliation between a threat issue and T2D threat utilizing the CMAverse R package deal120, adjusting age, intercourse, smoking, BMI and PA (if not the examined threat issue), calorie consumption and different cohort-specific covariates, individually in NHS/HPFS, SOL and WHI. We mixed whole, oblique and direct results, respectively, from every analytical set utilizing a fixed-effect meta-analysis. The mediated proportion was calculated by dividing oblique impact to whole impact. Metabolites with an oblique impact FDR < 0.05 and a constant impact course between the oblique and whole results, was thought of as a possible mediator between a threat issue and T2D threat.

A multimetabolite signature for incident T2D prediction

We used metabolites shared between the Broad Institute and the Metabolon platforms (excluding glucose) to develop the signature to extend its generalizability to future research. To keep away from overfitting in mannequin growth and testing, we employed a leave-one-cohort-out cross-validation strategy, through which we put aside one cohort because the testing set every time, and skilled a prediction mannequin for the set-aside cohort utilizing knowledge from all different cohorts (Extended Data Fig. 8). Given the heterogeneity of our cohorts, we didn’t pool individual-level knowledge for mannequin coaching. Instead, we utilized a two-step strategy to coach the prediction mannequin in a representable cohort (that’s, WHI, which assessed essentially the most shared metabolites for all its contributors) but in addition leveraged affiliation knowledge from a number of different cohorts. In every iteration (that’s, for every held-out testing cohort), we first carried out a metabolome-wide meta-analysis for T2D threat utilizing all cohorts besides WHI and the held-out cohort. Then, metabolites related to T2D threat at FDR < 0.05 in step one and shared between the 2 metabolomic platforms, have been used as enter in a Cox regression with elastic internet regularization, applied utilizing the glmnet R package deal121, to assemble a metabolomic signature mannequin for T2D prediction in WHI. The derived mannequin was additional utilized to the held-out cohort to calculate a metabolomic signature rating. Within WHI, a leave-one-out cross-validation strategy was used to amass the unbiased metabolomic signature rating. For particulars, please see Supplementary Methods.

The metabolomic signature scores, calculated in every held-out cohort, have been then standardized. To consider whether or not the signature improved the T2D threat prediction, we fitted three units of logistic (in SOL, and T2D nested case–management substudy in NHS) or Cox fashions (all different datasets): one mannequin together with solely the metabolomic signature; a traditional threat issue mannequin together with age, intercourse, smoking, lipid-lowering remedy use, anti-hypertensive medicines, household historical past of diabetes, hypertension, dyslipidemia and BMI; and a 3rd mannequin together with all typical threat components and the metabolomic signature. We in contrast the AUC between the traditional mannequin versus the traditional plus metabolomic signature mannequin. In a secondary evaluation, we additional included blood glucose (from metabolomic assays) within the typical mannequin to judge the added worth of the metabolomic signatures past blood glucose.

In every cohort, we calculated the crude incident charge of T2D throughout deciles of the signature rating. We fitted logistic or Cox fashions to investigate the relative threat of T2D, evaluating increased versus lowest deciles of the metabolomic signature, adjusting for a similar covariates in the primary evaluation mannequin 2. In NHS/HPFS, SOL and WHI, we examined associations between the metabolomic signature with baseline threat components, by regressing the signature rating on age, intercourse (if acceptable), present smoking standing, BMI, PA, intakes of 15 most important meals teams and fasting standing concurrently, along with cohort-specific covariates, utilizing linear regression. All evaluation was carried out individually in every cohort, and outcomes have been mixed utilizing a meta-analysis. FDR < 0.05 was thought of as statistically vital.

We carried out two sensitivity analyses throughout mannequin growth. One was to make use of SOL (measured essentially the most metabolites utilizing the Metabolon platform) because the consultant coaching cohort as an alternative of WHI, which confirmed an identical, albeit barely weaker, mannequin efficiency in held-out cohorts (Extended Data Fig. 8). The different was to match between elastic internet versus lasso regularizations121, which reaffirmed that elastic internet regression had suitable however a barely higher efficiency versus lasso regression (Supplementary Fig. 13). Separately from the leave-one-cohort-out cross-validation, we offered a remaining metabolomic signature mannequin for future research, developed utilizing knowledge from all research cohorts. For this mannequin, we first carried out a metabolome-wide meta-analysis for T2D threat in all cohorts besides WHI, after which used vital metabolites (FDR < 0.05) as enter in a Cox regression with elastic internet regularization for T2D prediction in WHI. The chosen metabolites and their coefficients of this remaining mannequin are extremely in keeping with these of fashions utilized to every held-out cohort (Supplementary Table 18a).

Reporting abstract

Further info on analysis design is offered within the Nature Portfolio Reporting Summary linked to this text.


This web page was created programmatically, to learn the article in its authentic location you may go to the hyperlink bellow:
https://www.nature.com/articles/s41591-025-04105-8
and if you wish to take away this text from our web site please contact us