Publications

Population-based association studies have identified many genetic risk loci for coronary artery disease (CAD), but it is often unclear how genes within these loci are linked to CAD. Here, we perform interaction proteomics for 11 CAD-risk genes to map their protein-protein interactions (PPIs) in human vascular cells and elucidate their roles in CAD. The resulting PPI networks contain interactions that are outside of known biology in the vasculature and are enriched for genes involved in immunity-related and arterial-wall-specific mechanisms. Several PPI networks derived from smooth muscle cells are significantly enriched for genetic variants associated with CAD and related vascular phenotypes. Furthermore, the networks identify 61 genes that are found in genetic loci associated with risk of CAD, prioritizing them as the causal candidates within these loci. These findings indicate that the PPI networks we have generated are a rich resource for guiding future research into the molecular pathogenesis of CAD.

Hydrogen sulfide (H2S) played a pivotal role in the early evolution of life on Earth before the predominance of atmospheric oxygen. The legacy of a persistent role for H2S in life's processes recently emerged through its discovery in modern biochemistry as an endogenous cellular signalling modulator involved in numerous biological processes. One major mechanism through which H2S signals is protein cysteine persulfidation, an oxidative post-translational modification. In recent years, chemoproteomic technologies have been developed to allow the global scanning of protein persulfidation targets in mammalian cells and tissues, providing a powerful tool to elucidate the broader impact of altered H2S in organismal physiological health and human disease states. While hundreds of proteins were confirmed to be persulfidated by global persulfidome methodologies, the targeting of specific proteins of interest and the investigation of further mechanistic studies are still underdeveloped due to a lack of stringent specificity of the methods and the inherent instability of persulfides. This review provides an overview of the processes of endogenous H2S production, oxidation, and signalling and highlights the application and limitations of current persulfidation labelling approaches for investigation of this important evolutionarily conserved biological switch for protein function.

We outline a robust, simple, and cost-effective method to simultaneously visualize all mouse lower hindlimb skeletal muscles. We describe procedures for orientating the whole lower hindlimb in gum tragacanth prior to freezing, simplifying the proceeding experimental steps, and enhancing the clarity and comprehensiveness of characterizations. We then detail steps for quantifying muscle fiber size and fiber type characteristics in a single cryosection using immunofluorescence and histochemistry. This protocol can be applicable for commonly used histological and (immuno) histochemical evaluations such muscle degeneration/regeneration, fibrosis, immune cell infiltration, enzymatic activity and glycogen content.

Complement proteins facilitate synaptic elimination during neurodevelopmental pruning, but neural complement regulation is not well understood. CUB and Sushi Multiple Domains 1 (CSMD1) can regulate complement activity in vitro, is expressed in the brain, and is associated with increased schizophrenia risk. Beyond this, little is known about CSMD1 including whether it regulates complement activity in the brain or otherwise plays a role in neurodevelopment. We used biochemical, immunohistochemical, and proteomic techniques to examine the regional, cellular, and subcellular distribution as well as protein interactions of CSMD1 in the brain. To evaluate whether CSMD1 is involved in complement-mediated synapse elimination, we examined Csmd1-knockout mice and CSMD1-knockout human stem cell-derived neurons. We interrogated synapse and circuit development of the mouse visual thalamus, a process that involves complement pathway activity. We also quantified complement deposition on synapses in mouse visual thalamus and on cultured human neurons. Finally, we assessed uptake of synaptosomes by cultured microglia. We found that CSMD1 is present at synapses and interacts with complement proteins in the brain. Mice lacking Csmd1 displayed increased levels of complement component C3, an increased colocalization of C3 with presynaptic terminals, fewer retinogeniculate synapses, and aberrant segregation of eye-specific retinal inputs to the visual thalamus during the critical period of complement-dependent refinement of this circuit. Loss of CSMD1 in vivo enhanced synaptosome engulfment by microglia in vitro, and this effect was dependent on activity of the microglial complement receptor, CR3. Finally, human stem cell-derived neurons lacking CSMD1 were more vulnerable to complement deposition. These data suggest that CSMD1 can function as a regulator of complement-mediated synapse elimination in the brain during development.

Broadening gene therapy applications requires manufacturable vectors that efficiently transduce target cells in humans and preclinical models. Conventional selections of adeno-associated virus (AAV) capsid libraries are inefficient at searching the vast sequence space for the small fraction of vectors possessing multiple traits essential for clinical translation. Here, we present Fit4Function, a generalizable machine learning (ML) approach for systematically engineering multi-trait AAV capsids. By leveraging a capsid library that uniformly samples the manufacturable sequence space, reproducible screening data are generated to train accurate sequence-to-function models. Combining six models, we designed a multi-trait (liver-targeted, manufacturable) capsid library and validated 88% of library variants on all six predetermined criteria. Furthermore, the models, trained only on mouse in vivo and human in vitro Fit4Function data, accurately predicted AAV capsid variant biodistribution in macaque. Top candidates exhibited production yields comparable to AAV9, efficient murine liver transduction, up to 1000-fold greater human hepatocyte transduction, and increased enrichment relative to AAV9 in a screen for liver transduction in macaques. The Fit4Function strategy ultimately makes it possible to predict cross-species traits of peptide-modified AAV capsids and is a critical step toward assembling an ML atlas that predicts AAV capsid performance across dozens of traits.

ABSTRACT: Cancer progression and response to therapy are inextricably reliant on the coevolution of a supportive tissue microenvironment. This is particularly evident in pancreatic ductal adenocarcinoma, a tumor type characterized by expansive and heterogeneous stroma. Herein, we employed single-cell RNA sequencing and spatial transcriptomics of normal, inflamed, and malignant pancreatic tissues to contextualize stromal dynamics associated with disease and treatment status, identifying temporal and spatial trajectories of fibroblast differentiation. Using analytical tools to infer cellular communication, together with a newly developed assay to annotate genomic alterations in cancer cells, we additionally explored the complex intercellular networks underlying tissue circuitry, highlighting a fibroblast-centric interactome that grows in strength and complexity in the context of malignant transformation. Our study yields new insights on the stromal remodeling events favoring the development of a tumor-supportive microenvironment and provides a powerful resource for the exploration of novel points of therapeutic intervention in pancreatic ductal adenocarcinoma.

SIGNIFICANCE: Pancreatic cancer remains a high unmet medical need. Understanding the interactions between stroma and cancer cells in this disease may unveil new opportunities for therapeutic intervention.

IMPORTANCE: Immune checkpoint inhibitors (ICIs) have revolutionized cancer care; however, accompanying immune-related adverse events (irAEs) confer substantial morbidity and occasional mortality. Life-threatening irAEs may require permanent cessation of ICI, even in patients with positive tumor response. Therefore, it is imperative to comprehensively define the spectrum of irAEs to aid individualized decision-making around the initiation of ICI therapy.

OBJECTIVE: To define incidence, risk factors, and clinical spectrum of an irreversible and life-threatening irAE: ICI-induced diabetes.

DESIGN, SETTING, AND PARTICIPANTS: This cohort study, conducted at an academic integrated health care system examined 14 328 adult patients treated with ICIs, including 64 patients who developed ICI-induced diabetes, from July 2010 to January 2022. The data were analyzed from 2022 to 2023. Cases of ICI-induced diabetes were manually confirmed; detailed clinical phenotyping was performed at diagnosis and 1-year follow-up. For 862 patients, genotyping data were available, and polygenic risk for type 1 diabetes was determined.

MAIN OUTCOMES AND MEASURES: For ICI-induced diabetes cases and controls, demographic characteristics, comorbidities, tumor category, and ICI category were compared. Among ICI-induced diabetes cases, markers of glycemic physiology were examined at diagnosis and 1-year follow-up. For patients with available genotyping, a published type 1 diabetes polygenic score (T1D GRS2) was calculated.

RESULTS: Of 14 328 participants, 6571 (45.9%) were women, and the median (range) age was 66 (8-106) years. The prevalence of ICI-induced diabetes among ICI-treated patients was 0.45% (64 of 14 328), with an incidence of 124.8 per 100 000 person-years. Preexisting type 2 diabetes (odds ratio [OR], 5.91; 95% CI, 3.34-10.45) and treatment with combination ICI (OR, 2.57; 95% CI, 1.44-4.59) were significant clinical risk factors of ICI-induced diabetes. T1D GRS2 was associated with ICI-induced diabetes risk, with an OR of 4.4 (95% CI, 1.8-10.5) for patients in the top decile of T1D GRS2, demonstrating a genetic association between spontaneous autoimmunity and irAEs. Patients with ICI-induced diabetes were in 3 distinct phenotypic categories based on autoantibodies and residual pancreatic function, with varying severity of initial presentation.

CONCLUSIONS AND RELEVANCE: The results of this analysis of 14 328 ICI-treated patients followed up from ICI initiation determined the incidence, risk factors and clinical spectrum of ICI-induced diabetes. Widespread implementation of this approach across organ-specific irAEs may enhance diagnosis and management of these conditions, and this becomes especially pertinent as ICI treatment rapidly expands to treat a wide spectrum of cancers and is used at earlier stages of treatment.

The prevalence of youth-onset type 2 diabetes (T2D) and childhood obesity has been rising steadily1, producing a growing public health concern1 that disproportionately affects minority groups2. The genetic basis of youth-onset T2D and its relationship to other forms of diabetes are unclear3. Here we report a detailed genetic characterization of youth-onset T2D by analysing exome sequences and common variant associations for 3,005 individuals with youth-onset T2D and 9,777 adult control participants matched for ancestry, including both males and females. We identify monogenic diabetes variants in 2.4% of individuals and three exome-wide significant (P < 2.6 × 10-6) gene-level associations (HNF1A, MC4R, ATXN2L). Furthermore, we report rare variant association enrichments within 25 gene sets related to obesity, monogenic diabetes and β-cell function. Many youth-onset T2D associations are shared with adult-onset T2D, but genetic risk factors of all frequencies-and rare variants in particular-are enriched within youth-onset T2D cases (5.0-fold increase in the rare variant and 3.4-fold increase in common variant genetic liability relative to adult-onset cases). The clinical presentation of participants with youth-onset T2D is influenced in part by the frequency of genetic risk factors within each individual. These findings portray youth-onset T2D as a heterogeneous disease situated on a spectrum between monogenic diabetes and adult-onset T2D.

Reduced insulin sensitivity (insulin resistance) is a hallmark of normal physiology in late pregnancy and also underlies gestational diabetes mellitus (GDM). We conducted transcriptomic profiling of 434 human placentas and identified a positive association between insulin-like growth factor binding protein 1 gene (IGFBP1) expression in the placenta and insulin sensitivity at  26 weeks gestation. Circulating IGFBP1 protein levels rose over the course of pregnancy and declined postpartum, which, together with high gene expression levels in our placenta samples, suggests a placental or decidual source. Higher circulating IGFBP1 levels were associated with greater insulin sensitivity (lesser insulin resistance) at  26 weeks gestation in the same cohort and in two additional pregnancy cohorts. In addition, low circulating IGFBP1 levels in early pregnancy predicted subsequent GDM diagnosis in two cohorts of pregnant women. These results implicate IGFBP1 in the glycemic physiology of pregnancy and suggest a role for placental IGFBP1 deficiency in GDM pathogenesis.

Partitioned polygenic scores (pPS) have been developed to capture pathophysiologic processes underlying type 2 diabetes (T2D). We investigated the association of T2D pPS with diabetes-related traits and T2D incidence in the Diabetes Prevention Program. We generated five T2D pPS (β-cell, proinsulin, liver/lipid, obesity, lipodystrophy) in 2,647 participants randomized to intensive lifestyle, metformin, or placebo arms. Associations were tested with general linear models and Cox regression with adjustment for age, sex, and principal components. Sensitivity analyses included adjustment for BMI. Higher β-cell pPS was associated with lower insulinogenic index and corrected insulin response at 1-year follow-up with adjustment for baseline measures (effect per pPS SD -0.04, P = 9.6 × 10-7, and -8.45 μU/mg, P = 5.6 × 10-6, respectively) and with increased diabetes incidence with adjustment for BMI at nominal significance (hazard ratio 1.10 per SD, P = 0.035). The liver/lipid pPS was associated with reduced 1-year baseline-adjusted triglyceride levels (effect per SD -4.37, P = 0.001). There was no significant interaction between T2D pPS and randomized groups. The remaining pPS were associated with baseline measures only. We conclude that despite interventions for diabetes prevention, participants with a high genetic burden of the β-cell cluster pPS had worsening in measures of β-cell function.