GLP1R Gene Expression and Kidney Disease Progression.

Importance: Glucagon-like peptide 1 receptor agonists (GLP-1RAs) may have nephroprotective properties beyond those related to weight loss and glycemic control. Objective: To investigate the association of genetically proxied GLP-1RAs with kidney disease progression. Design, Setting, and Participants: This genetic association study assembled a national retrospective cohort of veterans aged 18 years or older from the US Department of Veterans Affairs Million Veteran Program between January 10, 2011, and December 31, 2021. Data were analyzed from November 2023 to February 2024. Exposures: Genetic risk score for systemic GLP1R gene expression that was calculated for each study participant based on genetic variants associated with GLP1R mRNA levels across all tissue samples within the Genotype-Tissue Expression project. Main Outcomes and Measures: The primary composite outcome was incident end-stage kidney disease or a 40% decline in estimated glomerular filtration rate. Cox proportional hazards regression survival analysis assessed the association between genetically proxied GLP-1RAs and kidney disease progression. Results: Among 353 153 individuals (92.5% men), median age was 66 years (IQR, 58.0-72.0 years) and median follow-up was 5.1 years (IQR, 3.1-7.2 years). Overall, 25.7% had diabetes, and 45.0% had obesity. A total of 4.6% experienced kidney disease progression. Overall, higher genetic GLP1R gene expression was associated with a lower risk of kidney disease progression in the unadjusted model (hazard ratio [HR], 0.96; 95% CI, 0.92-0.99; P = .02) and in the fully adjusted model accounting for baseline patient characteristics, body mass index, and the presence or absence of diabetes (HR, 0.96; 95% CI, 0.92-1.00; P = .04). The results were similar in sensitivity analyses stratified by diabetes or obesity status. Conclusions and Relevance: In this genetic association study, higher GLP1R gene expression was associated with a small reduction in risk of kidney disease progression. These findings support pleiotropic nephroprotective mechanisms of GLP-1RAs independent of their effects on body weight and glycemic control.

Exome wide association study for blood lipids in 1,158,017 individuals from diverse populations.

Rare coding alleles play crucial roles in the molecular diagnosis of genetic diseases. However, the systemic identification of these alleles has been challenging due to their scarcity in the general population. Here, we discovered and characterized rare coding alleles contributing to genetic dyslipidemia, a principal risk for coronary artery disease, among over a million individuals combining three large contemporary genetic datasets (the Million Veteran Program, n = 634,535, UK Biobank, n = 431,178, and the All of Us Research Program, n = 92,304) totaling 1,158,017 multi-ancestral individuals. Unlike previous rare variant studies in lipids, this study included 238,243 individuals (20.6%) from non-European-like populations. Testing 2,997,401 rare coding variants from diverse backgrounds, we identified 800 exome-wide significant associations across 209 genes including 176 predicted loss of function and 624 missense variants. Among these exome-wide associations, 130 associations were driven by non-European-like populations. Associated alleles are highly enriched in functional variant classes, showed significant additive and recessive associations, exhibited similar effects across populations, and resolved pathogenicity for variants enriched in African or South-Asian populations. Furthermore, we identified 5 lipid-related genes associated with coronary artery disease (RORC, CFAP65, GTF2E2, PLCB3, and ZNF117). Among them, RORC is a potentially novel therapeutic target through the down regulation of LDLC by its silencing. This study provides resources and insights for understanding causal mechanisms, quantifying the expressivity of rare coding alleles, and identifying novel drug targets across diverse populations.

Military and occupational exposures among Veterans in the Million Veteran Program by survey self-report: a descriptive study.

OBJECTIVE: We aimed to characterise self-reported military and occupational exposures including Agent Orange, chemical/biological warfare agents, solvents, fuels, pesticides, metals and burn pits among Veterans in the Department of Veterans Affairs Million Veteran Program (MVP). METHODS: MVP is an ongoing longitudinal cohort and mega-biobank of over one million US Veterans. Over 500 000 MVP participants reported military exposures on the baseline survey, and over 300 000 reported occupational exposures on the lifestyle survey. We determined frequencies of selected self-reported occupational exposures by service era, specific deployment operation (1990-1991 Gulf War, Operation Enduring Freedom/Operation Iraqi Freedom (OEF/OIF)), service in a combat zone and occupational categories. We also explored differences in self-reported exposures by sex and race. RESULTS: Agent Orange exposure was mainly reported by Vietnam-era Veterans. Gulf War and OEF/OIF Veterans deployed to a combat zone were more likely to report exposures to burn pits, chemical/biological weapons, anthrax vaccination and pyridostigmine bromide pill intake as compared with non-combat deployers and those not deployed. Occupational categories related to combat (infantry, combat engineer and helicopter pilot) often had the highest percentages of self-reported exposures, whereas those in healthcare-related occupations (dentists, physicians and occupational therapists) tended to report exposures much less often. Self-reported exposures also varied by race and sex. CONCLUSIONS: Our results demonstrate that the distribution of self-reported exposures varied by service era, demographics, deployment, combat experience and military occupation in MVP. Overall, the pattern of findings was consistent with previous population-based studies of US military Veterans.

Recalibrating the Genetics and Epidemiology of Colorectal Cancer Consortium (GECCO) Environmental Risk Score for Use in US Veterans.

BACKGROUND: Risk for colorectal cancer (CRC) may accumulate through multiple environmental factors. Understanding their effects, along with genetics, age and family history, could allow improvements in clinical decisions for screening protocols. We aimed to extend previous work by recalibrating an environmental risk score (e-Score) for CRC among a sample of US Veteran participants of the Million Veteran Program (MVP). METHODS: Demographic, lifestyle, and CRC data from 2011-2022 were abstracted from survey responses and health records of 227,504 male MVP participants. Weighting for each environmental factor's effect size was recalculated using VA training data to create a recalibrated e-Score. This recalibrated score was compared with the original weighted e-Score in a validation sample of 113,752 (n cases=590). Nested multiple logistic regression models tested associations between quintiles for recalibrated and original e-Scores. Likelihood Ratio Tests were used to compare model performance. RESULTS: Age (p<0.0001), education (p<0.0001), diabetes (p<0.0001), physical activity (p<0.0001), smoking (p<0.0001), NSAID use (p<0.0001), calcium (p=0.015), folate (p=0.020), and fruit consumption (p=0.019) were significantly different between CRC case and control groups. In the validation sample, the recalibrated e-Score model significantly improved the base model performance (p<0.001), but the original e-Score model did not (p=0.07). The recalibrated e-Score model quintile 5 was associated with significantly higher odds for CRC compared with quintile 1 (Q5 vs Q1: 1.79, 95% CI: 1.38-2.33). CONCLUSIONS: Multiple environmental factors, and the recalibrated e-Score quintiles were significantly associated with CRC cases. IMPACT: A recalibrated, Veteran-specific e-Score could be used to help personalize CRC screening and prevention strategies.

Comparison of methods for building polygenic scores for diverse populations.

Polygenic scores (PGSs) are a promising tool for estimating individual-level genetic risk of disease based on the results of genome-wide association studies (GWASs). However, their promise has yet to be fully realized because most currently available PGSs were built with genetic data from predominantly European-ancestry populations, and PGS performance declines when scores are applied to target populations different from the populations from which they were derived. Thus, there is a great need to improve PGS performance in currently under-studied populations. In this work we leverage data from two large and diverse cohorts the Million Veterans Program (MVP) and All of Us (AoU), providing us the unique opportunity to compare methods for building PGSs for multi-ancestry populations across multiple traits. We build PGSs for five continuous traits and five binary traits using both multi-ancestry and single-ancestry approaches with popular Bayesian PGS methods and both MVP META GWAS results and population-specific GWAS results from the respective African, European, and Hispanic MVP populations. We evaluate these scores in three AoU populations genetically similar to the respective African, Admixed American, and European 1000 Genomes Project superpopulations. Using correlation-based tests, we make formal comparisons of the PGS performance across the multiple AoU populations. We conclude that approaches that combine GWAS data from multiple populations produce PGSs that perform better than approaches that utilize smaller single-population GWAS results matched to the target population, and specifically that multi-ancestry scores built with PRS-CSx outperform the other approaches in the three AoU populations.

Measuring Frailty Using Self-Report or Automated Tools to Identify Risk of Cardiovascular Events and Mortality: The Million Veteran Program.

BACKGROUND: Frailty, a syndrome of physiologic vulnerability, increases cardiovascular disease (CVD) risk. Whether in person or automated frailty tools are ideal for identifying CVD risk remains unclear. We calculated 3 distinct frailty scores and examined their associations with mortality and CVD events in the Million Veteran Program, a prospective cohort of nearly 1 million US veterans. METHODS AND RESULTS: Veterans aged ≥50 years and enrolled from 2011 to 2018 were included. Two frailty indices (FI) based on the deficit accumulation theory were calculated: the questionnaire-based 36-item Million Veteran Program-FI and 31-item Veterans Affairs-FI using claims data. We calculated Fried physical frailty using the self-reported, 3-item Study of Osteoporotic Fractures. Multivariable-adjusted Cox models examined the association of frailty by each score with primary (all-cause and CVD mortality) and secondary (myocardial infarction, stroke, and heart failure) outcomes. In 190 688 veterans (69±9 years, 94% male, 85% White), 33, 233 (17%) all-cause and 10 115 (5%) CVD deaths occurred. Using Million Veteran Program-FI, 29% were robust, 42% pre-frail, and 29% frail. Frailty prevalence increased by age group (27% in 50-59 to 42% in ≥90 years). Using the Million Veteran Program-FI, over 6±2 years, frail veterans had a higher hazard of all-cause (hazard ratio [HR], 3.05 [95% CI, 2.95-3.16]) and CVD mortality (HR, 3.65 [95% CI, 3.43-3.90]). Findings were concordant for the Veterans Affairs-FI and Study of Osteoporotic Fractures frailty definitions, and remained significant even among younger veterans aged 50-59 years. CONCLUSIONS: Irrespective of frailty measure, frailty is associated with a higher risk of all-cause mortality and adverse CVD events. Further study of frailty in veterans aged <60 years old is warranted.

Characterising the genetic architecture of changes in adiposity during adulthood using electronic health records.

Obesity is a heritable disease, characterised by excess adiposity that is measured by body mass index (BMI). While over 1,000 genetic loci are associated with BMI, less is known about the genetic contribution to adiposity trajectories over adulthood. We derive adiposity-change phenotypes from 24.5 million primary-care health records in over 740,000 individuals in the UK Biobank, Million Veteran Program USA, and Estonian Biobank, to discover and validate the genetic architecture of adiposity trajectories. Using multiple BMI measurements over time increases power to identify genetic factors affecting baseline BMI by 14%. In the largest reported genome-wide study of adiposity-change in adulthood, we identify novel associations with BMI-change at six independent loci, including rs429358 (APOE missense variant). The SNP-based heritability of BMI-change (1.98%) is 9-fold lower than that of BMI. The modest genetic correlation between BMI-change and BMI (45.2%) indicates that genetic studies of longitudinal trajectories could uncover novel biology of quantitative traits in adulthood.

Protein Identification for Stroke Progression via Mendelian Randomization in Million Veteran Program and UK Biobank.

BACKGROUND: Individuals who have experienced a stroke, or transient ischemic attack, face a heightened risk of future cardiovascular events. Identification of genetic and molecular risk factors for subsequent cardiovascular outcomes may identify effective therapeutic targets to improve prognosis after an incident stroke. METHODS: We performed genome-wide association studies for subsequent major adverse cardiovascular events (MACE; ncases=51 929; ncontrols=39 980) and subsequent arterial ischemic stroke (AIS; ncases=45 120; ncontrols=46 789) after the first incident stroke within the Million Veteran Program and UK Biobank. We then used genetic variants associated with proteins (protein quantitative trait loci) to determine the effect of 1463 plasma protein abundances on subsequent MACE using Mendelian randomization. RESULTS: Two variants were significantly associated with subsequent cardiovascular events: rs76472767 near gene RNF220 (odds ratio, 0.75 [95% CI, 0.64-0.85]; P=3.69×10-8) with subsequent AIS and rs13294166 near gene LINC01492 (odds ratio, 1.52 [95% CI, 1.37-1.67]; P=3.77×10-8) with subsequent MACE. Using Mendelian randomization, we identified 2 proteins with an effect on subsequent MACE after a stroke: CCL27 ([C-C motif chemokine 27], effect odds ratio, 0.77 [95% CI, 0.66-0.88]; adjusted P=0.05) and TNFRSF14 ([tumor necrosis factor receptor superfamily member 14], effect odds ratio, 1.42 [95% CI, 1.24-1.60]; adjusted P=0.006). These proteins are not associated with incident AIS and are implicated to have a role in inflammation. CONCLUSIONS: We found evidence that 2 proteins with little effect on incident stroke appear to influence subsequent MACE after incident AIS. These associations suggest that inflammation is a contributing factor to subsequent MACE outcomes after incident AIS and highlights potential novel targets.

Estimating the Effect of Bariatric Surgery on Cardiovascular Events Using Observational Data?

BACKGROUND: Observational studies have reported strongly protective effects of bariatric surgery on cardiovascular disease, but with oversimplified definitions of the intervention, eligibility criteria, and follow-up, which deviate from those in a randomized trial. We describe an attempt to estimate the effect of bariatric surgery on cardiovascular disease without introducing these sources of bias, which may not be entirely possible with existing observational data. METHODS: We propose two target trials among persons with diabetes: (1) bariatric operation (vs. no operation) among individuals who have undergone preoperative preparation (lifestyle modifications and screening) and (2) preoperative preparation and a bariatric operation (vs. neither preoperative nor operative component). We emulated both target trials using observational data of US veterans. RESULTS: Comparing bariatric surgery with no surgery (target trial #1; 8,087 individuals), the 7-year cardiovascular risk was 18.0% (95% CI = 6.9, 32.7) in the surgery group and 18.9% (95% CI = 17.7, 20.1) in the no-surgery group (risk difference -0.9, 95% CI = -12.0, 14.0). Comparing preoperative components plus surgery vs. neither (target trial #2; 10,065 individuals), the 7-year cardiovascular risk was 17.4% (95% CI = 13.6, 22.0) in the surgery group and 18.8% (95% CI = 17.8, 19.9) in the no-surgery group (risk difference -1.4, 95% CI = -5.1, 3.2). Body mass index and hemoglobin A1c were reduced with bariatric interventions in both emulations. CONCLUSIONS: Within limitations of available observational data, our estimates do not provide evidence that bariatric surgery reduces cardiovascular disease and support equipoise for a randomized trial of bariatric surgery for cardiovascular disease prevention.

SMIM1 absence is associated with reduced energy expenditure and excess weight.

BACKGROUND: Obesity rates have nearly tripled in the past 50 years, and by 2030 more than 1 billion individuals worldwide are projected to be obese. This creates a significant economic strain due to the associated non-communicable diseases. The root cause is an energy expenditure imbalance, owing to an interplay of lifestyle, environmental, and genetic factors. Obesity has a polygenic genetic architecture; however, single genetic variants with large effect size are etiological in a minority of cases. These variants allowed the discovery of novel genes and biology relevant to weight regulation and ultimately led to the development of novel specific treatments. METHODS: We used a case-control approach to determine metabolic differences between individuals homozygous for a loss-of-function genetic variant in the small integral membrane protein 1 (SMIM1) and the general population, leveraging data from five cohorts. Metabolic characterization of SMIM1-/- individuals was performed using plasma biochemistry, calorimetric chamber, and DXA scan. FINDINGS: We found that individuals homozygous for a loss-of-function genetic variant in SMIM1 gene, underlying the blood group Vel, display excess body weight, dyslipidemia, altered leptin to adiponectin ratio, increased liver enzymes, and lower thyroid hormone levels. This was accompanied by a reduction in resting energy expenditure. CONCLUSION: This research identified a novel genetic predisposition to being overweight or obese. It highlights the need to investigate the genetic causes of obesity to select the most appropriate treatment given the large cost disparity between them. FUNDING: This work was funded by the National Institute of Health Research, British Heart Foundation, and NHS Blood and Transplant.

Accelerating Genome- and Phenome-Wide Association Studies using GPUs - A case study using data from the Million Veteran Program.

The expansion of biobanks has significantly propelled genomic discoveries yet the sheer scale of data within these repositories poses formidable computational hurdles, particularly in handling extensive matrix operations required by prevailing statistical frameworks. In this work, we introduce computational optimizations to the SAIGE (Scalable and Accurate Implementation of Generalized Mixed Model) algorithm, notably employing a GPU-based distributed computing approach to tackle these challenges. We applied these optimizations to conduct a large-scale genome-wide association study (GWAS) across 2,068 phenotypes derived from electronic health records of 635,969 diverse participants from the Veterans Affairs (VA) Million Veteran Program (MVP). Our strategies enabled scaling up the analysis to over 6,000 nodes on the Department of Energy (DOE) Oak Ridge Leadership Computing Facility (OLCF) Summit High-Performance Computer (HPC), resulting in a 20-fold acceleration compared to the baseline model. We also provide a Docker container with our optimizations that was successfully used on multiple cloud infrastructures on UK Biobank and All of Us datasets where we showed significant time and cost benefits over the baseline SAIGE model.

Characterization of Post-COVID-19 Definitions and Clinical Coding Practices: Longitudinal Study.

BACKGROUND: Post-COVID-19 condition (colloquially known as "long COVID-19") characterized as postacute sequelae of SARS-CoV-2 has no universal clinical case definition. Recent efforts have focused on understanding long COVID-19 symptoms, and electronic health record (EHR) data provide a unique resource for understanding this condition. The introduction of the International Classification of Diseases, Tenth Revision (ICD-10) code U09.9 for "Post COVID-19 condition, unspecified" to identify patients with long COVID-19 has provided a method of evaluating this condition in EHRs; however, the accuracy of this code is unclear. OBJECTIVE: This study aimed to characterize the utility and accuracy of the U09.9 code across 3 health care systems-the Veterans Health Administration, the Beth Israel Deaconess Medical Center, and the University of Pittsburgh Medical Center-against patients identified with long COVID-19 via a chart review by operationalizing the World Health Organization (WHO) and Centers for Disease Control and Prevention (CDC) definitions. METHODS: Patients who were COVID-19 positive with either a U07.1 ICD-10 code or positive polymerase chain reaction test within these health care systems were identified for chart review. Among this cohort, we sampled patients based on two approaches: (1) with a U09.9 code and (2) without a U09.9 code but with a new onset long COVID-19-related ICD-10 code, which allows us to assess the sensitivity of the U09.9 code. To operationalize the long COVID-19 definition based on health agency guidelines, symptoms were grouped into a "core" cluster of 11 commonly reported symptoms among patients with long COVID-19 and an extended cluster that captured all other symptoms by disease domain. Patients having ≥2 symptoms persisting for ≥60 days that were new onset after their COVID-19 infection, with ≥1 symptom in the core cluster, were labeled as having long COVID-19 per chart review. The code's performance was compared across 3 health care systems and across different time periods of the pandemic. RESULTS: Overall, 900 patient charts were reviewed across 3 health care systems. The prevalence of long COVID-19 among the cohort with the U09.9 ICD-10 code based on the operationalized WHO definition was between 23.2% and 62.4% across these health care systems. We also evaluated a less stringent version of the WHO definition and the CDC definition and observed an increase in the prevalence of long COVID-19 at all 3 health care systems. CONCLUSIONS: This is one of the first studies to evaluate the U09.9 code against a clinical case definition for long COVID-19, as well as the first to apply this definition to EHR data using a chart review approach on a nationwide cohort across multiple health care systems. This chart review approach can be implemented at other EHR systems to further evaluate the utility and performance of the U09.9 code.

Meta-Analysis of Genome-Wide Association Studies Reveals Genetic Mechanisms of Supraventricular Arrhythmias.

BACKGROUND: Substantial data support a heritable basis for supraventricular tachycardias, but the genetic determinants and molecular mechanisms of these arrhythmias are poorly understood. We sought to identify genetic loci associated with atrioventricular nodal reentrant tachycardia (AVNRT) and atrioventricular accessory pathways or atrioventricular reciprocating tachycardia (AVAPs/AVRT). METHODS: We performed multiancestry meta-analyses of genome-wide association studies to identify genetic loci for AVNRT (4 studies) and AVAP/AVRT (7 studies). We assessed evidence supporting the potential causal effects of candidate genes by analyzing relations between associated variants and cardiac gene expression, performing transcriptome-wide analyses, and examining prior genome-wide association studies. RESULTS: Analyses comprised 2384 AVNRT cases and 106 489 referents, and 2811 AVAP/AVRT cases and 1,483 093 referents. We identified 2 significant loci for AVNRT, which implicate NKX2-5 and TTN as disease susceptibility genes. A transcriptome-wide association analysis supported an association between reduced predicted cardiac expression of NKX2-5 and AVNRT. We identified 3 significant loci for AVAP/AVRT, which implicate SCN5A, SCN10A, and TTN/CCDC141. Variant associations at several loci have been previously reported for cardiac phenotypes, including atrial fibrillation, stroke, Brugada syndrome, and electrocardiographic intervals. CONCLUSIONS: Our findings highlight gene regions associated with ion channel function (AVAP/AVRT), as well as cardiac development and the sarcomere (AVAP/AVRT and AVNRT) as important potential effectors of supraventricular tachycardia susceptibility.

Efficacy of serial casting protocols in idiopathic toe-walking.

OBJECTIVE: Idiopathic toe-walking (ITW) is a diagnosis of exclusion. A relationship between ITW and decreased range of motion (ROM) is postulated. Treatments focus on increasing ankle dorsiflexion including serial casting. There is no consensus for duration of serial casting. This study aimed to determine ROM changes with cast change intervals of one vs. two weeks, and the rate of ITW recurrence. METHODS: This was a retrospective study of 86 patients, ages 0-9 years with ITW undergoing weekly casting (N = 29) and two-week casting (N = 57) at a children's hospital from 2014-2020. ROM at baseline, two weeks, four weeks, and final cast removal were collected. Statistical analyses included chi-squared tests, two-sample t-tests, and linear mixed regression. Outcome distributions were assessed for normality. P-values < 0.05 were considered statistically significant. RESULTS: After adjusting for baseline ROM, the mean change in ROM from baseline to two weeks was 10.6∘ vs 7.5∘ in the one-week vs. two-week casting interval, respectively (p < 0.001). The baseline to final measurement was 13.4∘ vs 9.8∘ in the one-week vs. two-week casting interval, respectively (p < 0.001). The rate of recurrence of ITW was similar between the two groups. CONCLUSION: This study suggests greater improvement in ROM in the one-week vs. two-week casting interval group.

Type 1 Diabetes Genetic Risk in 109,954 Veterans With Adult-Onset Diabetes: The Million Veteran Program (MVP).

OBJECTIVE: To characterize high type 1 diabetes (T1D) genetic risk in a population where type 2 diabetes (T2D) predominates. RESEARCH DESIGN AND METHODS: Characteristics typically associated with T1D were assessed in 109,594 Million Veteran Program participants with adult-onset diabetes, 2011-2021, who had T1D genetic risk scores (GRS) defined as low (0 to <45%), medium (45 to <90%), high (90 to <95%), or highest (≥95%). RESULTS: T1D characteristics increased progressively with higher genetic risk (P < 0.001 for trend). A GRS ≥90% was more common with diabetes diagnoses before age 40 years, but 95% of those participants were diagnosed at age ≥40 years, and their characteristics resembled those of individuals with T2D in mean age (64.3 years) and BMI (32.3 kg/m2). Compared with the low-risk group, the highest-risk group was more likely to have diabetic ketoacidosis (low GRS 0.9% vs. highest GRS 3.7%), hypoglycemia prompting emergency visits (3.7% vs. 5.8%), outpatient plasma glucose <50 mg/dL (7.5% vs. 13.4%), a shorter median time to start insulin (3.5 vs. 1.4 years), use of a T1D diagnostic code (16.3% vs. 28.1%), low C-peptide levels if tested (1.8% vs. 32.4%), and glutamic acid decarboxylase antibodies (6.9% vs. 45.2%), all P < 0.001. CONCLUSIONS: Characteristics associated with T1D were increased with higher genetic risk, and especially with the top 10% of risk. However, the age and BMI of those participants resemble those of people with T2D, and a substantial proportion did not have diagnostic testing or use of T1D diagnostic codes. T1D genetic screening could be used to aid identification of adult-onset T1D in settings in which T2D predominates.

Neurological diagnoses in hospitalized COVID-19 patients associated with adverse outcomes: A multinational cohort study.

Few studies examining the patient outcomes of concurrent neurological manifestations during acute COVID-19 leveraged multinational cohorts of adults and children or distinguished between central and peripheral nervous system (CNS vs. PNS) involvement. Using a federated multinational network in which local clinicians and informatics experts curated the electronic health records data, we evaluated the risk of prolonged hospitalization and mortality in hospitalized COVID-19 patients from 21 healthcare systems across 7 countries. For adults, we used a federated learning approach whereby we ran Cox proportional hazard models locally at each healthcare system and performed a meta-analysis on the aggregated results to estimate the overall risk of adverse outcomes across our geographically diverse populations. For children, we reported descriptive statistics separately due to their low frequency of neurological involvement and poor outcomes. Among the 106,229 hospitalized COVID-19 patients (104,031 patients ≥18 years; 2,198 patients <18 years, January 2020-October 2021), 15,101 (14%) had at least one CNS diagnosis, while 2,788 (3%) had at least one PNS diagnosis. After controlling for demographics and pre-existing conditions, adults with CNS involvement had longer hospital stay (11 versus 6 days) and greater risk of (Hazard Ratio = 1.78) and faster time to death (12 versus 24 days) than patients with no neurological condition (NNC) during acute COVID-19 hospitalization. Adults with PNS involvement also had longer hospital stay but lower risk of mortality than the NNC group. Although children had a low frequency of neurological involvement during COVID-19 hospitalization, a substantially higher proportion of children with CNS involvement died compared to those with NNC (6% vs 1%). Overall, patients with concurrent CNS manifestation during acute COVID-19 hospitalization faced greater risks for adverse clinical outcomes than patients without any neurological diagnosis. Our global informatics framework using a federated approach (versus a centralized data collection approach) has utility for clinical discovery beyond COVID-19.

Heterogeneous associations between interleukin-6 receptor variants and phenotypes across ancestries and implications for therapy.

The Phenome-Wide Association Study (PheWAS) is increasingly used to broadly screen for potential treatment effects, e.g., IL6R variant as a proxy for IL6R antagonists. This approach offers an opportunity to address the limited power in clinical trials to study differential treatment effects across patient subgroups. However, limited methods exist to efficiently test for differences across subgroups in the thousands of multiple comparisons generated as part of a PheWAS. In this study, we developed an approach that maximizes the power to test for heterogeneous genotype-phenotype associations and applied this approach to an IL6R PheWAS among individuals of African (AFR) and European (EUR) ancestries. We identified 29 traits with differences in IL6R variant-phenotype associations, including a lower risk of type 2 diabetes in AFR (OR 0.96) vs EUR (OR 1.0, p-value for heterogeneity = 8.5 × 10-3), and higher white blood cell count (p-value for heterogeneity = 8.5 × 10-131). These data suggest a more salutary effect of IL6R blockade for T2D among individuals of AFR vs EUR ancestry and provide data to inform ongoing clinical trials targeting IL6 for an expanding number of conditions. Moreover, the method to test for heterogeneity of associations can be applied broadly to other large-scale genotype-phenotype screens in diverse populations.

Lipid levels and risk of acute pancreatitis using bidirectional Mendelian randomization.

Previous studies found lipid levels, especially triglycerides (TG), are associated with acute pancreatitis, but their causalities and bi-directions were not fully examined. We determined whether abnormal levels of TG, high-density lipoprotein cholesterol (HDL-C), and low-density lipoprotein cholesterol (LDL-C) are precursors and/or consequences of acute pancreatitis using bidirectional two-sample Mendelian randomization (MR) with two non-overlapping genome-wide association study (GWAS) summary statistics for lipid levels and acute pancreatitis. We found phenotypic associations that both higher TG levels and lower HDL-C levels contributed to increased risk of acute pancreatitis. Our GWAS meta-analysis of acute pancreatitis identified seven independent signals. Genetically predicted TG was positively associated with acute pancreatitis when using the variants specifically associated with TG using univariable MR [Odds ratio (OR), 95% CI 2.02, 1.22-3.31], but the reversed direction from acute pancreatitis to TG was not observed (mean difference = 0.003, SE = 0.002, P-value = 0.138). However, a bidirectional relationship of HDL-C and acute pancreatitis was observed: A 1-SD increment of genetically predicted HDL-C was associated with lower risk of acute pancreatitis (OR, 95% CI 0.84, 0.76-0.92) and genetically predisposed individuals with acute pancreatitis have, on average, 0.005 SD lower HDL-C (mean difference = - 0.005, SE = 0.002, P-value = 0.004). Our MR analysis confirms the evidence of TG as a risk factor of acute pancreatitis but not a consequence. A potential bidirectional relationship of HDL-C and acute pancreatitis occurs and raises the prospect of HDL-C modulation in the acute pancreatitis prevention and treatment.

Long-Term Body Mass Index Variability and Adverse Cardiovascular Outcomes.

Importance: Body mass index (BMI; calculated as weight in kilograms divided by height in meters squared) is a commonly used estimate of obesity, which is a complex trait affected by genetic and lifestyle factors. Marked weight gain and loss could be associated with adverse biological processes. Objective: To evaluate the association between BMI variability and incident cardiovascular disease (CVD) events in 2 distinct cohorts. Design, Setting, and Participants: This cohort study used data from the Million Veteran Program (MVP) between 2011 and 2018 and participants in the UK Biobank (UKB) enrolled between 2006 and 2010. Participants were followed up for a median of 3.8 (5th-95th percentile, 3.5) years. Participants with baseline CVD or cancer were excluded. Data were analyzed from September 2022 and September 2023. Exposure: BMI variability was calculated by the retrospective SD and coefficient of variation (CV) using multiple clinical BMI measurements up to the baseline. Main Outcomes and Measures: The main outcome was incident composite CVD events (incident nonfatal myocardial infarction, acute ischemic stroke, and cardiovascular death), assessed using Cox proportional hazards modeling after adjustment for CVD risk factors, including age, sex, mean BMI, systolic blood pressure, total cholesterol, high-density lipoprotein cholesterol, smoking status, diabetes status, and statin use. Secondary analysis assessed whether associations were dependent on the polygenic score of BMI. Results: Among 92 363 US veterans in the MVP cohort (81 675 [88%] male; mean [SD] age, 56.7 [14.1] years), there were 9695 Hispanic participants, 22 488 non-Hispanic Black participants, and 60 180 non-Hispanic White participants. A total of 4811 composite CVD events were observed from 2011 to 2018. The CV of BMI was associated with 16% higher risk for composite CVD across all groups (hazard ratio [HR], 1.16; 95% CI, 1.13-1.19). These associations were unchanged among subgroups and after adjustment for the polygenic score of BMI. The UKB cohort included 65 047 individuals (mean [SD] age, 57.30 (7.77) years; 38 065 [59%] female) and had 6934 composite CVD events. Each 1-SD increase in BMI variability in the UKB cohort was associated with 8% increased risk of cardiovascular death (HR, 1.08; 95% CI, 1.04-1.11). Conclusions and Relevance: This cohort study found that among US veterans, higher BMI variability was a significant risk marker associated with adverse cardiovascular events independent of mean BMI across major racial and ethnic groups. Results were consistent in the UKB for the cardiovascular death end point. Further studies should investigate the phenotype of high BMI variability.

Centralized Interactive Phenomics Resource: an integrated online phenomics knowledgebase for health data users.

OBJECTIVE: Development of clinical phenotypes from electronic health records (EHRs) can be resource intensive. Several phenotype libraries have been created to facilitate reuse of definitions. However, these platforms vary in target audience and utility. We describe the development of the Centralized Interactive Phenomics Resource (CIPHER) knowledgebase, a comprehensive public-facing phenotype library, which aims to facilitate clinical and health services research. MATERIALS AND METHODS: The platform was designed to collect and catalog EHR-based computable phenotype algorithms from any healthcare system, scale metadata management, facilitate phenotype discovery, and allow for integration of tools and user workflows. Phenomics experts were engaged in the development and testing of the site. RESULTS: The knowledgebase stores phenotype metadata using the CIPHER standard, and definitions are accessible through complex searching. Phenotypes are contributed to the knowledgebase via webform, allowing metadata validation. Data visualization tools linking to the knowledgebase enhance user interaction with content and accelerate phenotype development. DISCUSSION: The CIPHER knowledgebase was developed in the largest healthcare system in the United States and piloted with external partners. The design of the CIPHER website supports a variety of front-end tools and features to facilitate phenotype development and reuse. Health data users are encouraged to contribute their algorithms to the knowledgebase for wider dissemination to the research community, and to use the platform as a springboard for phenotyping. CONCLUSION: CIPHER is a public resource for all health data users available at https://phenomics.va.ornl.gov/ which facilitates phenotype reuse, development, and dissemination of phenotyping knowledge.