Exome sequence analysis identifies rare coding variants associated with a machine learning-based marker for coronary artery disease.

Coronary artery disease (CAD) exists on a spectrum of disease represented by a combination of risk factors and pathogenic processes. An in silico score for CAD built using machine learning and clinical data in electronic health records captures disease progression, severity and underdiagnosis on this spectrum and could enhance genetic discovery efforts for CAD. Here we tested associations of rare and ultrarare coding variants with the in silico score for CAD in the UK Biobank, All of Us Research Program and BioMe Biobank. We identified associations in 17 genes; of these, 14 show at least moderate levels of prior genetic, biological and/or clinical support for CAD. We also observed an excess of ultrarare coding variants in 321 aggregated CAD genes, suggesting more ultrarare variant associations await discovery. These results expand our understanding of the genetic etiology of CAD and illustrate how digital markers can enhance genetic association investigations for complex diseases.

Lipid-siRNA conjugate accesses a perivascular transport mechanism and achieves widespread and durable knockdown in the central nervous system.

Short-interfering RNA (siRNA) has gained significant interest for treatment of neurological diseases by providing the capacity to achieve sustained inhibition of nearly any gene target. Yet, achieving efficacious drug delivery throughout deep brain structures of the CNS remains a considerable hurdle. We herein describe a lipid-siRNA conjugate that, following delivery into the cerebrospinal fluid (CSF), is transported effectively through perivascular spaces, enabling broad dispersion within CSF compartments and through the CNS parenchyma. We provide a detailed examination of the temporal kinetics of gene silencing, highlighting potent knockdown for up to five months from a single injection without detectable toxicity. Single-cell RNA sequencing further demonstrates gene silencing activity across diverse cell populations in the parenchyma and at brain borders, which may provide new avenues for neurological disease-modifying therapies.

Genotype-guided prescribing predictors in CYP2C19 intermediate metabolizers receiving percutaneous coronary intervention.

Background: Previous differences in guideline recommendation strength for CYP2C19 intermediate metabolizers may have limited genotype (PGx)-optimal post-percutaneous coronary intervention antiplatelet prescribing. Results: In this single-center retrospective observational cohort study of CYP2C19 intermediate metabolizers, patients prescribed PGx-optimal therapy were younger and less likely on anticoagulation (2 vs 12%; p = 0.006). More patients prescribed PGx-optimal therapy possessed commercial insurance (36 vs 7%; p < 0.001), which was a predictor for PGx-optimal selection (OR: 6.464; 95% CI: 2.386-17.516; p < 0.001). Conclusion: Anticoagulation use was significantly associated with clopidogrel use (OR: 0.138; 95% CI: 0.026-0.730; p = 0.020). No statistical difference in composite major adverse cardiovascular events (5 vs 14%; p = 0.173) or bleeding (8 vs 6%; Not significant) was observed between PGx-optimal and PGx-suboptimal therapy.

Not all CYP2C19 intermediate metabolizers undergoing PCI are prescribed genotype-optimal P2Y12 antiplatelet therapy. Commercial insurance and no anticoagulant were found to be associated with ticagrelor and prasugrel prescribing in this population.

Anandamide is an Early Blood Biomarker of Hermansky-Pudlak Syndrome Pulmonary Fibrosis.

Hermansky-Pudlak syndrome (HPS) is a group of rare genetic disorders, with several subtypes leading to fatal adult-onset pulmonary fibrosis (PF) and no effective treatment. Circulating biomarkers detecting early PF have not been identified. We investigated whether endocannabinoids could serve as blood biomarkers of PF in HPS. We measured endocannabinoids in the serum of HPS, IPF, and healthy human subjects and in a mouse model of HPSPF. Pulmonary function tests (PFT) were correlated with endocannabinoid measurements. In a pale ear mouse model of bleomycin-induced HPSPF, serum endocannabinoid levels were measured with and without treatment with zevaquenabant (MRI-1867), a peripheral CB1R and iNOS antagonist. In three separate cohorts, circulating anandamide levels were increased in HPS-1 patients with or without PF, compared to healthy volunteers. This increase was not observed in IPF patients or in HPS-3 patients, who do not have PF. Circulating anandamide (AEA) levels were negatively correlated with PFT. Furthermore, a longitudinal study over the course of 5-14 years with HPS-1 patients indicated that circulating AEA levels begin to increase with the fibrotic lung process even at the subclinical stages of HPSPF. In pale ear mice with bleomycin-induced HpsPF, serum AEA levels were significantly increased in the earliest stages of PF and remained elevated at a later fibrotic stage. Zevaquenabant treatment reduced the increased AEA levels and attenuated progression in bleomycin-induced HpsPF. Circulating AEA may be a prognostic blood biomarker for PF in HPS-1 patients. Further studies are indicated to evaluate endocannabinoids as potential surrogate biomarkers in progressive fibrotic lung diseases.

Massively parallel in vivo Perturb-seq reveals cell-type-specific transcriptional networks in cortical development.

Leveraging AAVs' versatile tropism and labeling capacity, we expanded the scale of in vivo CRISPR screening with single-cell transcriptomic phenotyping across embryonic to adult brains and peripheral nervous systems. Through extensive tests of 86 vectors across AAV serotypes combined with a transposon system, we substantially amplified labeling efficacy and accelerated in vivo gene delivery from weeks to days. Our proof-of-principle in utero screen identified the pleiotropic effects of Foxg1, highlighting its tight regulation of distinct networks essential for cell fate specification of Layer 6 corticothalamic neurons. Notably, our platform can label >6% of cerebral cells, surpassing the current state-of-the-art efficacy at <0.1% by lentivirus, to achieve analysis of over 30,000 cells in one experiment and enable massively parallel in vivo Perturb-seq. Compatible with various phenotypic measurements (single-cell or spatial multi-omics), it presents a flexible approach to interrogate gene function across cell types in vivo, translating gene variants to their causal function.

Muesli Intake May Protect Against Coronary Artery Disease: Mendelian Randomization on 13 Dietary Traits.

BACKGROUND: Diet is a key modifiable risk factor of coronary artery disease (CAD). However, the causal effects of specific dietary traits on CAD risk remain unclear. With the expansion of dietary data in population biobanks, Mendelian randomization (MR) could help enable the efficient estimation of causality in diet-disease associations. OBJECTIVES: The primary goal was to test causality for 13 common dietary traits on CAD risk using a systematic 2-sample MR framework. A secondary goal was to identify plasma metabolites mediating diet-CAD associations suspected to be causal. METHODS: Cross-sectional genetic and dietary data on up to 420,531 UK Biobank and 184,305 CARDIoGRAMplusC4D individuals of European ancestry were used in 2-sample MR. The primary analysis used fixed effect inverse-variance weighted regression, while sensitivity analyses used weighted median estimation, MR-Egger regression, and MR-Pleiotropy Residual Sum and Outlier. RESULTS: Genetic variants serving as proxies for muesli intake were negatively associated with CAD risk (OR: 0.74; 95% CI: 0.65-0.84; P = 5.385 × 10-4). Sensitivity analyses using weighted median estimation supported this with a significant association in the same direction. Additionally, we identified higher plasma acetate levels as a potential mediator (OR: 0.03; 95% CI: 0.01-0.12; P = 1.15 × 10-4). CONCLUSIONS: Muesli, a mixture of oats, seeds, nuts, dried fruit, and milk, may causally reduce CAD risk. Circulating levels of acetate, a gut microbiota-derived short-chain fatty acid, could be mediating its cardioprotective effects. These findings highlight the role of gut flora in cardiovascular health and help prioritize randomized trials on dietary interventions for CAD.

Genome-first evaluation with exome sequence and clinical data uncovers underdiagnosed genetic disorders in a large healthcare system.

Population-based genomic screening may help diagnose individuals with disease-risk variants. Here, we perform a genome-first evaluation for nine disorders in 29,039 participants with linked exome sequences and electronic health records (EHRs). We identify 614 individuals with 303 pathogenic/likely pathogenic or predicted loss-of-function (P/LP/LoF) variants, yielding 644 observations; 487 observations (76%) lack a corresponding clinical diagnosis in the EHR. Upon further investigation, 75 clinically undiagnosed observations (15%) have evidence of symptomatic untreated disease, including familial hypercholesterolemia (3 of 6 [50%] undiagnosed observations with disease evidence) and breast cancer (23 of 106 [22%]). These genetic findings enable targeted phenotyping that reveals new diagnoses in previously undiagnosed individuals. Disease yield is greater with variants in penetrant genes for which disease is observed in carriers in an independent cohort. The prevalence of P/LP/LoF variants exceeds that of clinical diagnoses, and some clinically undiagnosed carriers are discovered to have disease. These results highlight the potential of population-based genomic screening.

Disrupting CD38-driven T cell dysfunction restores sensitivity to cancer immunotherapy.

A central problem in cancer immunotherapy with immune checkpoint blockade (ICB) is the development of resistance, which affects 50% of patients with metastatic melanoma1,2. T cell exhaustion, resulting from chronic antigen exposure in the tumour microenvironment, is a major driver of ICB resistance3. Here, we show that CD38, an ecto-enzyme involved in nicotinamide adenine dinucleotide (NAD+) catabolism, is highly expressed in exhausted CD8+ T cells in melanoma and is associated with ICB resistance. Tumour-derived CD38hiCD8+ T cells are dysfunctional, characterised by impaired proliferative capacity, effector function, and dysregulated mitochondrial bioenergetics. Genetic and pharmacological blockade of CD38 in murine and patient-derived organotypic tumour models (MDOTS/PDOTS) enhanced tumour immunity and overcame ICB resistance. Mechanistically, disrupting CD38 activity in T cells restored cellular NAD+ pools, improved mitochondrial function, increased proliferation, augmented effector function, and restored ICB sensitivity. Taken together, these data demonstrate a role for the CD38-NAD+ axis in promoting T cell exhaustion and ICB resistance, and establish the efficacy of CD38 directed therapeutic strategies to overcome ICB resistance using clinically relevant, patient-derived 3D tumour models.

Development of a human genetics-guided priority score for 19,365 genes and 399 drug indications.

Studies have shown that drug targets with human genetic support are more likely to succeed in clinical trials. Hence, a tool integrating genetic evidence to prioritize drug target genes is beneficial for drug discovery. We built a genetic priority score (GPS) by integrating eight genetic features with drug indications from the Open Targets and SIDER databases. The top 0.83%, 0.28% and 0.19% of the GPS conferred a 5.3-, 9.9- and 11.0-fold increased effect of having an indication, respectively. In addition, we observed that targets in the top 0.28% of the score were 1.7-, 3.7- and 8.8-fold more likely to advance from phase I to phases II, III and IV, respectively. Complementary to the GPS, we incorporated the direction of genetic effect and drug mechanism into a directional version of the score called the GPS with direction of effect. We applied our method to 19,365 protein-coding genes and 399 drug indications and made all results available through a web portal.

Massively parallel in vivo Perturb-seq reveals cell type-specific transcriptional networks in cortical development.

Systematic analysis of gene function across diverse cell types in vivo is hindered by two challenges: obtaining sufficient cells from live tissues and accurately identifying each cell's perturbation in high-throughput single-cell assays. Leveraging AAV's versatile cell type tropism and high labeling capacity, we expanded the resolution and scale of in vivo CRISPR screens: allowing phenotypic analysis at single-cell resolution across a multitude of cell types in the embryonic brain, adult brain, and peripheral nervous system. We undertook extensive tests of 86 AAV serotypes, combined with a transposon system, to substantially amplify labeling and accelerate in vivo gene delivery from weeks to days. Using this platform, we performed an in utero genetic screen as proof-of-principle and identified pleiotropic regulatory networks of Foxg1 in cortical development, including Layer 6 corticothalamic neurons where it tightly controls distinct networks essential for cell fate specification. Notably, our platform can label >6% of cerebral cells, surpassing the current state-of-the-art efficacy at <0.1% (mediated by lentivirus), and achieve analysis of over 30,000 cells in one experiment, thus enabling massively parallel in vivo Perturb-seq. Compatible with various perturbation techniques (CRISPRa/i) and phenotypic measurements (single-cell or spatial multi-omics), our platform presents a flexible, modular approach to interrogate gene function across diverse cell types in vivo, connecting gene variants to their causal functions.

In vitro and in vivo pharmacokinetic characterization, chiral conversion and PBPK scaling towards human PK simulation of S-MRI-1867, a drug candidate for Hermansky-Pudlak syndrome pulmonary fibrosis.

Hermansky-Pudlak syndrome (HPS) is a rare autosomal recessive disorder that affects lysosome-related organelles, often leading to fatal pulmonary fibrosis (PF). The search for a treatment for HPS pulmonary fibrosis (HPSPF) is ongoing. S-MRI-1867, a dual cannabinoid receptor 1 (CB1R)/inducible nitric oxide synthase (iNOS) inhibitor, has shown great promise for the treatment of several fibrotic diseases, including HPSPF. In this study, we investigated the in vitro ADME characteristics of S-MRI-1867, as well as its pharmacokinetic (PK) properties in mice, rats, dogs, and monkeys. S-MRI-1867 showed low aqueous solubility (< 1 µg/mL), high plasma protein binding (>99%), and moderate to high metabolic stability. In its preclinical PK studies, S-MRI-1867 exhibited moderate to low plasma clearance (CLp) and high steady-state volume of distribution (Vdss) across all species. Despite the low solubility and P-gp efflux, S-MRI-1867 showed great permeability and metabolic stability leading to a moderate bioavailability (21-60%) across mouse, rat, dog, and monkey. Since the R form of MRI-1867 is CB1R-inactive, we investigated the potential conversion of S-MRI-1867 to R-MRI-1867 in mice and found that the chiral conversion was negligible. Furthermore, we developed and validated a PBPK model that adequately fits the PK profiles of S-MRI-1867 in mice, rats, dogs, and monkeys using various dosing regimens. We employed this PBPK model to simulate the human PK profiles of S-MRI-1867, enabling us to inform human dose selection and support the advancement of this promising drug candidate in the treatment of HPSPF.

Characterization of Pik1 function in fission yeast reveals its conserved role in lipid synthesis and not cytokinesis.

Phosphatidylinositol (PI)-4-phosphate (PI4P) is a lipid found at the plasma membrane (PM) and Golgi in cells from yeast to humans. PI4P is generated from PI by PI4-kinases and can be converted into PI-4,5-bisphosphate [PI(4,5)P2]. Schizosaccharomyces pombe have two essential PI4-kinases - Stt4 and Pik1. Stt4 localizes to the PM, and its loss from the PM results in a decrease of PM PI4P and PI(4,5)P2. As a result, cells divide non-medially due to disrupted cytokinetic ring-PM anchoring. However, the localization and function of S. pombe Pik1 has not been thoroughly examined. Here, we found that Pik1 localizes exclusively to the trans-Golgi and is required for Golgi PI4P production. We determined that Ncs1 regulates Pik1, but unlike in other organisms, it is not required for Pik1 Golgi localization. When Pik1 function was disrupted, PM PI4P but not PI(4,5)P2 levels were reduced, a major difference compared with Stt4. We conclude that Stt4 is the chief enzyme responsible for producing the PI4P that generates PI(4,5)P2. Also, that cells with disrupted Pik1 do not divide asymmetrically highlights the specific importance of PM PI(4,5)P2 for cytokinetic ring-PM anchoring.

Characterization of Pik1 function in fission yeast reveals its conserved role in lipid synthesis and not cytokinesis.

Phosphatidylinositol (PI)-4-phosphate (PI4P) is a lipid found at the plasma membrane (PM) and Golgi in cells from yeast to humans. PI4P is generated from PI by PI4-kinases and can be converted to PI-4,5-bisphosphate [PI(4,5)P 2 ]. Schizosaccharomyces pombe have 2 essential PI4-kinases: Stt4 and Pik1. Stt4 localizes to the PM and its loss from the PM results in a decrease of PM PI4P and PI(4,5)P 2 . As a result, cells divide non-medially due to disrupted cytokinetic ring-PM anchoring. However, the localization and function of S. pombe Pik1 has not been thoroughly examined. Here, we found that Pik1 localizes exclusively to the trans-Golgi and is required for Golgi PI4P production. We determined that Ncs1 regulates Pik1, but unlike in other organisms, it is not required for Pik1 Golgi localization. When Pik1 function was disrupted, PM PI4P but not PI(4,5)P 2 levels were reduced, a major difference with Stt4. We conclude that Stt4 is the chief enzyme responsible for producing the PI4P that generates PI(4,5)P 2 . Also, that cells with disrupted Pik1 do not divide asymmetrically highlights the specific importance of PM PI(4,5)P 2 for cytokinetic ring-PM anchoring. Summary statement: Fission yeast Pik1 localizes exclusively to the trans-Golgi independently of Ncs1, where it contributes to PI4P but not PI(4,5)P 2 synthesis. Pik1 does not affect cytokinesis.

A machine learning model identifies patients in need of autoimmune disease testing using electronic health records.

Systemic autoimmune rheumatic diseases (SARDs) can lead to irreversible damage if left untreated, yet these patients often endure long diagnostic journeys before being diagnosed and treated. Machine learning may help overcome the challenges of diagnosing SARDs and inform clinical decision-making. Here, we developed and tested a machine learning model to identify patients who should receive rheumatological evaluation for SARDs using longitudinal electronic health records of 161,584 individuals from two institutions. The model demonstrated high performance for predicting cases of autoantibody-tested individuals in a validation set, an external test set, and an independent cohort with a broader case definition. This approach identified more individuals for autoantibody testing compared with current clinical standards and a greater proportion of autoantibody carriers among those tested. Diagnoses of SARDs and other autoimmune conditions increased with higher model probabilities. The model detected a need for autoantibody testing and rheumatology encounters up to five years before the test date and assessment date, respectively. Altogether, these findings illustrate that the clinical manifestations of a diverse array of autoimmune conditions are detectable in electronic health records using machine learning, which may help systematize and accelerate autoimmune testing.

An Integrative Multiomics Framework for Identification of Therapeutic Targets in Pulmonary Fibrosis.

Pulmonary fibrosis (PF) is a heterogeneous disease with a poor prognosis. Therefore, identifying additional therapeutic modalities is required to improve outcome. However, the lack of biomarkers of disease progression hampers the preclinical to clinical translational process. Here, this work assesses and identifies progressive alterations in pulmonary function, transcriptomics, and metabolomics in the mouse lung at 7, 14, 21, and 28 days after a single dose of oropharyngeal bleomycin. By integrating multi-omics data, this work identifies two central gene subnetworks associated with multiple critical pathological changes in transcriptomics and metabolomics as well as pulmonary function. This work presents a multi-omics-based framework to establish a translational link between the bleomycin-induced PF model in mice and human idiopathic pulmonary fibrosis to identify druggable targets and test therapeutic candidates. This work also indicates peripheral cannabinoid receptor 1 (CB1 R) antagonism as a rational therapeutic target for clinical translation in PF. Mouse Lung Fibrosis Atlas can be accessed freely at https://niaaa.nih.gov/mouselungfibrosisatlas.

Customized gaming system engages young children in reaching and balance training.

Purpose: Trunk stability, an important prerequisite for many activities of daily living, can be impaired in children with movement disorders. Current treatment options can be costly and fail to fully engage young participants. We developed an affordable, smart screen-based intervention and tested if it engages young children in physical therapy goal driven exercises. Methods: Here we describe the ADAPT system, Aiding Distanced and Accessible Physical Therapy, which is a large touch-interactive device with customizable games. One such game, "Bubble Popper," encourages high repetitions of weight shifts, reaching, and balance training as the participant pops bubbles in sitting, kneeling, or standing positions. Results: Sixteen participants aged 2-18 years were tested during physical therapy sessions. The number of screen touches and length of game play indicate high participant engagement. In trials lasting less than 3 min, on average, older participants (12-18 years) made 159 screen touches per trial while the younger participants (2-7 years) made 97. In a 30-min session, on average, older participants actively played the game for 12.49 min while younger participants played for 11.22 min. Conclusion: The ADAPT system is a feasible means to engage young participants in reaching and balance training during physical therapy.

Phenome-wide Mendelian randomization study of plasma triglyceride levels and 2600 disease traits.

Background: Causality between plasma triglyceride (TG) levels and atherosclerotic cardiovascular disease (ASCVD) risk remains controversial despite more than four decades of study and two recent landmark trials, STRENGTH, and REDUCE-IT. Further unclear is the association between TG levels and non-atherosclerotic diseases across organ systems. Methods: Here, we conducted a phenome-wide, two-sample Mendelian randomization (MR) analysis using inverse-variance weighted (IVW) regression to systematically infer the causal effects of plasma TG levels on 2600 disease traits in the European ancestry population of UK Biobank. For replication, we externally tested 221 nominally significant associations (p<0.05) in an independent cohort from FinnGen. To account for potential horizontal pleiotropy and the influence of invalid instrumental variables, we performed sensitivity analyses using MR-Egger regression, weighted median estimator, and MR-PRESSO. Finally, we used multivariable MR (MVMR) controlling for correlated lipid fractions to distinguish the independent effect of plasma TG levels. Results: Our results identified seven disease traits reaching Bonferroni-corrected significance in both the discovery (p<1.92 × 10-5) and replication analyses (p<2.26 × 10-4), suggesting a causal relationship between plasma TG levels and ASCVDs, including coronary artery disease (OR 1.33, 95% CI 1.24-1.43, p=2.47 × 10-13). We also identified 12 disease traits that were Bonferroni-significant in the discovery or replication analysis and at least nominally significant in the other analysis (p<0.05), identifying plasma TG levels as a novel potential risk factor for nine non-ASCVD diseases, including uterine leiomyoma (OR 1.19, 95% CI 1.10-1.29, p=1.17 × 10-5). Conclusions: Taking a phenome-wide, two-sample MR approach, we identified causal associations between plasma TG levels and 19 disease traits across organ systems. Our findings suggest unrealized drug repurposing opportunities or adverse effects related to approved and emerging TG-lowering agents, as well as mechanistic insights for future studies. Funding: RD is supported by the National Institute of General Medical Sciences of the National Institutes of Health (NIH) (R35-GM124836) and the National Heart, Lung, and Blood Institute of the NIH (R01-HL139865 and R01-HL155915).

A TREM2-activating antibody with a blood-brain barrier transport vehicle enhances microglial metabolism in Alzheimer's disease models.

Loss-of-function variants of TREM2 are associated with increased risk of Alzheimer's disease (AD), suggesting that activation of this innate immune receptor may be a useful therapeutic strategy. Here we describe a high-affinity human TREM2-activating antibody engineered with a monovalent transferrin receptor (TfR) binding site, termed antibody transport vehicle (ATV), to facilitate blood-brain barrier transcytosis. Upon peripheral delivery in mice, ATV:TREM2 showed improved brain biodistribution and enhanced signaling compared to a standard anti-TREM2 antibody. In human induced pluripotent stem cell (iPSC)-derived microglia, ATV:TREM2 induced proliferation and improved mitochondrial metabolism. Single-cell RNA sequencing and morphometry revealed that ATV:TREM2 shifted microglia to metabolically responsive states, which were distinct from those induced by amyloid pathology. In an AD mouse model, ATV:TREM2 boosted brain microglial activity and glucose metabolism. Thus, ATV:TREM2 represents a promising approach to improve microglial function and treat brain hypometabolism found in patients with AD.

Mental health of undergraduates one year after the start of the COVID-19 pandemic: Findings from the national college health assessment III.

OBJECTIVE: To compare the mental health of undergraduates before the COVID-19 pandemic lockdown to their mental health one year later. PARTICIPANTS: Data from the American College Health Association (ACHA)'s National College Health Assessment III (ACHA-NCHAIII) were used, averaging a sample size of 54,844 undergraduate students and 106 schools nationwide per time point of assessment in the study. METHODS: Secondary analyses of the ACHA-NCHAIII compared undergraduates' scores on five measures of mental health measures (loneliness, psychological distress, suicidality, flourishing, and resilience) from Spring 2020 to Spring 2021. RESULTS: Undergraduates' responses showed an increase in loneliness, psychological distress, and suicidality as well as a decrease in flourishing and resilience. CONCLUSIONS: The worsening of undergraduates' mental health calls for greater action by schools to alleviate students' distress and improve their wellbeing.