Coordinated immune dysregulation in Juvenile Dermatomyositis revealed by single-cell genomics.

Juvenile Dermatomyositis (JDM) is one of several childhood-onset autoimmune disorders characterized by a type I interferon response and autoantibodies. Treatment options are limited due to incomplete understanding of how the disease emerges from dysregulated cell states across the immune system. We therefore investigated the blood of JDM patients at different stages of disease activity using single-cell transcriptomics paired with surface protein expression. By immunophenotyping peripheral blood mononuclear cells, we observed skewing of the B cell compartment towards an immature naive state as a hallmark of JDM at diagnosis. Furthermore, we find that these changes in B cells are paralleled by T cell signatures suggestive of Th2-mediated inflammation that persist despite disease quiescence. We applied network analysis to reveal that hyperactivation of the type I interferon response in all immune populations is coordinated with previously masked cell states including dysfunctional protein processing in CD4+ T cells and regulation of cell death programming in NK, CD8+ T cells and gdT cells. Together, these findings unveil the coordinated immune dysregulation underpinning JDM and provide insight into strategies for restoring balance in immune function.

Exposure to autoimmune disorders increases Alzheimer's disease risk in a multi-site electronic health record analysis.

Molecular studies of Alzheimer's disease (AD) implicate potential links between autoimmunity and AD, but the underlying clinical relationships between these conditions remain poorly understood. Electronic health records (EHRs) provide an opportunity to determine the clinical risk relationship between autoimmune disorders and AD and understand whether specific disorders and disorder subtypes affect AD risk at the phenotypic level in human populations. We evaluated relationships between 26 autoimmune disorders and AD across retrospective observational case-control and cohort study designs in the EHR systems at UCSF and Stanford. We quantified overall and sex-specific AD risk effects that these autoimmune disorders confer. We identified significantly increased AD risk in autoimmune disorder patients in both study designs at UCSF and at Stanford. This pattern was driven by specific autoimmunity subtypes including endocrine, gastrointestinal, dermatologic, and musculoskeletal disorders. We also observed increased AD risk from autoimmunity in both women and men, but women with autoimmune disorders continued to have a higher AD prevalence than men, indicating persistent sex-specificity. This study identifies autoimmune disorders as strong risk factors for AD that validate across several study designs and EHR databases. It sets the foundation for exploring how underlying autoimmune mechanisms increase AD risk and contribute to AD pathogenesis.

Data-driven insights can transform women's reproductive health.

This perspective explores the transformative potential of data-driven insights to understand and address women's reproductive health conditions. Historically, clinical studies often excluded women, hindering comprehensive research into conditions such as adverse pregnancy outcomes and endometriosis. Recent advances in technology (e.g., next-generation sequencing techniques, electronic medical records (EMRs), computational power) provide unprecedented opportunities for research in women's reproductive health. Studies of molecular data, including large-scale meta-analyses, provide valuable insights into conditions like preterm birth and preeclampsia. Moreover, EMRs and other clinical data sources enable researchers to study populations of individuals, uncovering trends and associations in women's reproductive health conditions. Despite these advancements, challenges such as data completeness, accuracy, and representation persist. We emphasize the importance of holistic approaches, greater inclusion, and refining and expanding on how we leverage data and computational integrative approaches for discoveries so that we can benefit not only women's reproductive health but overall human health.

Identifying therapeutic candidates for endometriosis through a transcriptomics-based drug repositioning approach.

Existing medical treatments for endometriosis-related pain are often ineffective, underscoring the need for new therapeutic strategies. In this study, we applied a computational drug repurposing pipeline to stratified and unstratified disease signatures based on endometrial gene expression data to identify potential therapeutics from existing drugs, based on expression reversal. Of 3,131 unique genes differentially expressed by at least one of six endometriosis signatures, only 308 (9.8%) were in common; however, 221 out of 299 drugs identified, (73.9%) were shared. We selected fenoprofen, an uncommonly prescribed NSAID that was the top therapeutic candidate for further investigation. When testing fenoprofen in an established rat model of endometriosis, fenoprofen successfully alleviated endometriosis-associated vaginal hyperalgesia, a surrogate marker for endometriosis-related pain. These findings validate fenoprofen as a therapeutic that could be utilized more frequently for endometriosis and suggest the utility of the aforementioned computational drug repurposing approach for endometriosis.

Gene expression meta-analysis reveals aging and cellular senescence signatures in scleroderma-associated interstitial lung disease.

Aging and cellular senescence are increasingly recognized as key contributors to pulmonary fibrosis. However, our understanding in the context of scleroderma-associated interstitial lung disease (SSc-ILD) is limited. To investigate, we leveraged previously established lung aging- and cell-specific senescence signatures to determine their presence and potential relevance to SSc-ILD. We performed a gene expression meta-analysis of lung tissues from 38 SSc-ILD and 18 healthy controls and found that markers (GDF15, COMP, and CDKN2A) and pathways (p53) of senescence were significantly increased in SSc-ILD. When probing the established aging and cellular senescence signatures, we found that epithelial and fibroblast senescence signatures had a 3.6- and 3.7-fold enrichment, respectively, in the lung tissue of SSc-ILD and that lung aging genes (CDKN2A, FRZB, PDE1A, and NAPI12) were increased in SSc-ILD. These signatures were also enriched in SSc skin and associated with degree of skin involvement (limited vs. diffuse cutaneous). To further support these findings, we examined telomere length (TL), a surrogate for aging, in the lung tissue and found that, independent of age, SSc-ILD had significantly shorter telomeres than controls in type II alveolar cells in the lung. TL in SSc-ILD was comparable to idiopathic pulmonary fibrosis, a disease of known aberrant aging. Taken together, this study provides novel insight into the possible mechanistic effects of accelerated aging and aberrant cellular senescence in SSc-ILD pathogenesis.

Miscarriage risk assessment: a bioinformatic approach to identifying candidate lethal genes and variants.

PURPOSE: Miscarriage, often resulting from a variety of genetic factors, is a common pregnancy outcome. Preconception genetic carrier screening (PGCS) identifies at-risk partners for newborn genetic disorders; however, PGCS panels currently lack miscarriage-related genes. In this study, we evaluated the potential impact of both known and candidate genes on prenatal lethality and the effectiveness of PGCS in diverse populations. METHODS: We analyzed 125,748 human exome sequences and mouse and human gene function databases. Our goals were to identify genes crucial for human fetal survival (lethal genes), to find variants not present in a homozygous state in healthy humans, and to estimate carrier rates of known and candidate lethal genes in various populations and ethnic groups. RESULTS: This study identified 138 genes in which heterozygous lethal variants are present in the general population with a frequency of 0.5% or greater. Screening for these 138 genes could identify 4.6% (in the Finnish population) to 39.8% (in the East Asian population) of couples at risk of miscarriage. This explains the cause of pregnancy loss in approximately 1.1-10% of cases affected by biallelic lethal variants. CONCLUSION: This study has identified a set of genes and variants potentially associated with lethality across different ethnic backgrounds. The variation of these genes across ethnic groups underscores the need for a comprehensive, pan-ethnic PGCS panel that includes genes related to miscarriage.

Leveraging electronic health records and knowledge networks for Alzheimer's disease prediction and sex-specific biological insights.

Identification of Alzheimer's disease (AD) onset risk can facilitate interventions before irreversible disease progression. We demonstrate that electronic health records from the University of California, San Francisco, followed by knowledge networks (for example, SPOKE) allow for (1) prediction of AD onset and (2) prioritization of biological hypotheses, and (3) contextualization of sex dimorphism. We trained random forest models and predicted AD onset on a cohort of 749 individuals with AD and 250,545 controls with a mean area under the receiver operating characteristic of 0.72 (7 years prior) to 0.81 (1 day prior). We further harnessed matched cohort models to identify conditions with predictive power before AD onset. Knowledge networks highlight shared genes between multiple top predictors and AD (for example, APOE, ACTB, IL6 and INS). Genetic colocalization analysis supports AD association with hyperlipidemia at the APOE locus, as well as a stronger female AD association with osteoporosis at a locus near MS4A6A. We therefore show how clinical data can be utilized for early AD prediction and identification of personalized biological hypotheses.

Discovery of sparse, reliable omic biomarkers with Stabl.

Adoption of high-content omic technologies in clinical studies, coupled with computational methods, has yielded an abundance of candidate biomarkers. However, translating such findings into bona fide clinical biomarkers remains challenging. To facilitate this process, we introduce Stabl, a general machine learning method that identifies a sparse, reliable set of biomarkers by integrating noise injection and a data-driven signal-to-noise threshold into multivariable predictive modeling. Evaluation of Stabl on synthetic datasets and five independent clinical studies demonstrates improved biomarker sparsity and reliability compared to commonly used sparsity-promoting regularization methods while maintaining predictive performance; it distills datasets containing 1,400-35,000 features down to 4-34 candidate biomarkers. Stabl extends to multi-omic integration tasks, enabling biological interpretation of complex predictive models, as it hones in on a shortlist of proteomic, metabolomic and cytometric events predicting labor onset, microbial biomarkers of pre-term birth and a pre-operative immune signature of post-surgical infections. Stabl is available at https://github.com/gregbellan/Stabl .

Osteosarcoma PDX-Derived Cell Line Models for Preclinical Drug Evaluation Demonstrate Metastasis Inhibition by Dinaciclib through a Genome-Targeted Approach.

PURPOSE: Models to study metastatic disease in rare cancers are needed to advance preclinical therapeutics and to gain insight into disease biology. Osteosarcoma is a rare cancer with a complex genomic landscape in which outcomes for patients with metastatic disease are poor. As osteosarcoma genomes are highly heterogeneous, multiple models are needed to fully elucidate key aspects of disease biology and to recapitulate clinically relevant phenotypes. EXPERIMENTAL DESIGN: Matched patient samples, patient-derived xenografts (PDX), and PDX-derived cell lines were comprehensively evaluated using whole-genome sequencing and RNA sequencing. The in vivo metastatic phenotype of the PDX-derived cell lines was characterized in both an intravenous and an orthotopic murine model. As a proof-of-concept study, we tested the preclinical effectiveness of a cyclin-dependent kinase inhibitor on the growth of metastatic tumors in an orthotopic amputation model. RESULTS: PDXs and PDX-derived cell lines largely maintained the expression profiles of the patient from which they were derived despite the emergence of whole-genome duplication in a subset of cell lines. The cell lines were heterogeneous in their metastatic capacity, and heterogeneous tissue tropism was observed in both intravenous and orthotopic models. Single-agent dinaciclib was effective at dramatically reducing the metastatic burden. CONCLUSIONS: The variation in metastasis predilection sites between osteosarcoma PDX-derived cell lines demonstrates their ability to recapitulate the spectrum of the disease observed in patients. We describe here a panel of new osteosarcoma PDX-derived cell lines that we believe will be of wide use to the osteosarcoma research community.

Exploring applications of non-targeted analysis in the characterization of the prenatal exposome.

Capturing the breadth of chemical exposures in utero is critical in understanding their long-term health effects for mother and child. We explored methodological adaptations in a Non-Targeted Analysis (NTA) pipeline and evaluated the effects on chemical annotation and discovery for maternal and infant exposure. We focus on lesser-known/underreported chemicals in maternal and umbilical cord serum analyzed with liquid chromatography-quadrupole time-of-flight mass spectrometry (LC-QTOF/MS). The samples were collected from a demographically diverse cohort of 296 maternal-cord pairs (n = 592) recruited in San Francisco Bay area. We developed and evaluated two data processing pipelines, primarily differing by detection frequency cut-off, to extract chemical features from non-targeted analysis (NTA). We annotated the detected chemical features by matching with EPA CompTox Chemicals Dashboard (n = 860,000 chemicals) and Human Metabolome Database (n = 3140 chemicals) and applied a Kendrick Mass Defect filter to detect homologous series. We collected fragmentation spectra (MS/MS) on a subset of serum samples and matched to an experimental MS/MS database within the MS-Dial website and other experimental MS/MS spectra collected from standards in our lab. We annotated ~72 % of the features (total features = 32,197, levels 1-4). We confirmed 22 compounds with analytical standards, tentatively identified 88 compounds with MS/MS spectra, and annotated 4862 exogenous chemicals with an in-house developed annotation algorithm. We detected 36 chemicals that appear to not have been previously reported in human blood and 9 chemicals that were reported in less than five studies. Our findings underline the importance of NTA in the discovery of lesser-known/unreported chemicals important to characterize human exposures.

Cell-Specific Transposable Element Gene Expression Analysis Identifies Associations with Systemic Lupus Erythematosus Phenotypes.

There is an established yet unexplained link between interferon (IFN) and systemic lupus erythematosus (SLE). The expression of sequences derived from transposable elements (TEs) may contribute to production of type I IFNs and generation of autoantibodies. We profiled cell-sorted RNA-seq data (CD4+ T cells, CD14+ monocytes, CD19+ B cells, and NK cells) from PBMCs of 120 SLE patients and quantified TE expression identifying 27,135 TEs. We tested for differential TE expression across 10 SLE phenotypes including autoantibody production and disease activity and discovered 731 differentially expressed (DE) TEs whose effects were mostly cell-specific and phenotype-specific. DE TEs were enriched for specific families and viral genes encoded in TE sequences. Increased expression of DE TEs was associated with genes involved in antiviral activity such as LY6E, ISG15, TRIM22 and pathways such as interferon signaling. These findings suggest that expression of TEs contributes to activation of SLE-related mechanisms in a cell-specific manner, which can impact disease diagnostics and therapeutics.

Associations with spontaneous and indicated preterm birth in a densely phenotyped EHR cohort.

Background: Preterm birth (PTB) is the leading cause of infant mortality and follows multiple biological pathways, many of which are poorly understood. Some PTBs result from medically indicated labor following complications from hypertension and/or diabetes, while many others are spontaneous with unknown causes. Previously, investigation of potential risk factors has been limited by lack of data on maternal medical history and the difficulty of classifying PTBs as indicated or spontaneous. Here, we leverage electronic health record (EHR) data (patient health information including demographics, diagnoses, and medications) and a supplemental curated pregnancy database to overcome these limitations. Novel associations may provide new insight into the pathophysiology of PTB as well as help identify individuals who would be at risk of PTB. Methods: We quantified associations between maternal diagnoses and preterm birth using logistic regression controlling for maternal age and socioeconomic factors within a University of California, San Francisco (UCSF), EHR cohort with 10,643 births ( nterm = 9692, nspontaneous_preterm = 449, nindicated_preterm = 418) and maternal pre-conception diagnosis phenotypes derived from International Classification of Diseases (ICD) 9 and 10 codes. Results: Eighteen conditions significantly and robustly (False Discovery Rate (FDR)<0.05) associated with PTBs compared to term. We discovered known (hypertension, diabetes, and chronic kidney disease) and less established (blood, cardiac, gynecological, and liver conditions) associations. Type 1 diabetes was the most significant overall association (adjusted p = 1.6×10 -14 , adjusted OR = 7 (95% CI 5, 12)), and the odds ratios for the significant phenotypes ranged from 3 to 13. We further carried out analysis stratified by spontaneous vs. indicated PTB. No phenotypes significantly associated with spontaneous PTB; however, the results for indicated PTB largely recapitulated the phenotype associations with all PTBs. Conclusions: Our study underscores the limitations of approaches that combine indicated and spontaneous births together. When combined, significant associations were almost entirely driven by indicated PTBs, although our spontaneous and indicated groups were of a similar size. Investigating the spontaneous population has the potential to reveal new pathways and understanding of the heterogeneity of PTB.

MaLiAmPi enables generalizable and taxonomy-independent microbiome features from technically diverse 16S-based microbiome studies.

For studies using microbiome data, the ability to robustly combine data from technically and biologically distinct microbiome studies is a crucial means of supporting more robust and clinically relevant inferences. Formidable technical challenges arise when attempting to combine data from technically diverse 16S rRNA gene variable region amplicon sequencing (16S) studies. Closed operational taxonomic units and taxonomy are criticized as being heavily dependent upon reference sets and with limited precision relative to the underlying biology. Phylogenetic placement has been demonstrated to be a promising taxonomy-free manner of harmonizing microbiome data, but it has lacked a validated count-based feature suitable for use in machine learning and association studies. Here we introduce a phylogenetic-placement-based, taxonomy-independent, compositional feature of microbiota: phylotypes. Phylotypes were predictive of clinical outcomes such as obesity or pre-term birth on technically diverse independent validation sets harmonized post hoc. Thus, phylotypes enable the rigorous cross-validation of 16S-based clinical prognostic models and associative microbiome studies.

Sex-stratified phenotyping of comorbidities associated with an inpatient delirium diagnosis using real world data.

Delirium is a heterogeneous and detrimental mental condition often seen in older, hospitalized patients and is currently hard to predict. In this study, we leverage large-scale, real- world data using the electronic health records (EHR) to identify two cohorts comprised of 7,492 UCSF patients and 19,417 UC health system patients (excluding UCSF patients) with an inpatient delirium diagnosis and the same number of propensity score-matched control patients without delirium. We found significant associations between comorbidities or laboratory test values and an inpatient delirium diagnosis which were validated independently. Most of these associations were those previously-identified as risk factors for delirium, including metabolic abnormalities, mental health diagnoses, and infections. Some of the associations were sex- specific, including those related to dementia subtypes and infections. We further explored the diagnostic associations with anemia and bipolar disorder by conducting longitudinal analyses from the time of first diagnosis of the risk factor to development of delirium demonstrating a significant relationship across time. Finally, we show that an inpatient delirium diagnosis leads to dramatic increases in mortality outcome across both cohorts. These results demonstrate the powerful application of leveraging EHR data to shed insights into prior diagnoses and laboratory test values that could help predict development of inpatient delirium and emphasize the importance of considering patient demographic characteristics including documented sex when making these assessments. One Sentence Summary: Longitudinal analysis of electronic health record data reveals associations between inpatient delirium, comorbidities, and mortality.

Coordinated immune dysregulation in Juvenile Dermatomyositis revealed by single-cell genomics.

Juvenile Dermatomyositis (JDM) is one of several childhood-onset autoimmune disorders characterized by a type I interferon response and autoantibodies. Treatment options are limited due to incomplete understanding of how the disease emerges from dysregulated cell states across the immune system. We therefore investigated the blood of JDM patients at different stages of disease activity using single-cell transcriptomics paired with surface protein expression. By immunophenotyping peripheral blood mononuclear cells, we observed skewing of the B cell compartment towards an immature naive state as a hallmark of JDM. Furthermore, we find that these changes in B cells are paralleled by signatures of Th2-mediated inflammation. Additionally, our work identified SIGLEC-1 expression in monocytes as a composite measure of heterogeneous type I interferon activity in disease. We applied network analysis to reveal that hyperactivation of the type I interferon response in all immune populations is coordinated with dysfunctional protein processing and regulation of cell death programming. This analysis separated the ubiquitously expressed type I interferon response into a central hub and revealed previously masked cell states. Together, these findings reveal the coordinated immune dysregulation underpinning JDM and provide novel insight into strategies for restoring balance in immune function.

Gene Expression Meta-Analysis Reveals Aging and Cellular Senescence Signatures in Scleroderma-associated Interstitial Lung Disease.

Aging and cellular senescence are increasingly recognized as key contributors to pulmonary fibrosis. However, our understanding in the context of scleroderma associated interstitial lung disease (SSc-ILD) is limited. To investigate, we leveraged previously established lung aging and cell-specific senescence signatures to determine their presence and potential relevance to SSc-ILD. We performed a gene expression meta-analysis of lung tissue from 38 SSc-ILD and 18 healthy controls and found markers (GDF15, COMP, CDKN2A) and pathways (p53) of senescence were significantly increased in SSc-ILD. When probing the established aging and cellular senescence signatures, we found epithelial and fibroblast senescence signatures had a 3.6-fold and 3.7-fold enrichment respectively in the lung tissue of SSc-ILD and that lung aging genes ( CDKN2A, FRZB, PDE1A, NAPI12) were increased in SSc-ILD. These signatures were also enriched in SSc skin and associated with degree of skin involvement (limited vs. diffuse cutaneous). To further support these findings, we examined telomere length (TL), a surrogate for aging, in lung tissue and found independent of age, SSc-ILD had significantly shorter telomeres than controls in type II alveolar cells in the lung. TL in SSc-ILD was comparable to idiopathic pulmonary fibrosis, a disease of known aberrant aging. Taken together, this study provides novel insight into the possible mechanistic effects of accelerated aging and aberrant cellular senescence in SSc-ILD pathogenesis.

Global endometrial DNA methylation analysis reveals insights into mQTL regulation and associated endometriosis disease risk and endometrial function.

Endometriosis is a leading cause of pain and infertility affecting millions of women globally. Herein, we characterize variation in DNA methylation (DNAm) and its association with menstrual cycle phase, endometriosis, and genetic variants through analysis of genotype data and methylation in endometrial samples from 984 deeply-phenotyped participants. We estimate that 15.4% of the variation in endometriosis is captured by DNAm and identify significant differences in DNAm profiles associated with stage III/IV endometriosis, endometriosis sub-phenotypes and menstrual cycle phase, including opening of the window for embryo implantation. Menstrual cycle phase was a major source of DNAm variation suggesting cellular and hormonally-driven changes across the cycle can regulate genes and pathways responsible for endometrial physiology and function. DNAm quantitative trait locus (mQTL) analysis identified 118,185 independent cis-mQTLs including 51 associated with risk of endometriosis, highlighting candidate genes contributing to disease risk. Our work provides functional evidence for epigenetic targets contributing to endometriosis risk and pathogenesis. Data generated serve as a valuable resource for understanding tissue-specific effects of methylation on endometrial biology in health and disease.

Endometriosis in the era of precision medicine and impact on sexual and reproductive health across the lifespan and in diverse populations.

Endometriosis is a common estrogen-dependent disorder wherein uterine lining tissue (endometrium) is found mainly in the pelvis where it causes inflammation, chronic pelvic pain, pain with intercourse and menses, and infertility. Recent evidence also supports a systemic inflammatory component that underlies associated co-morbidities, e.g., migraines and cardiovascular and autoimmune diseases. Genetics and environment contribute significantly to disease risk, and with the explosion of omics technologies, underlying mechanisms of symptoms are increasingly being elucidated, although novel and effective therapeutics for pain and infertility have lagged behind these advances. Moreover, there are stark disparities in diagnosis, access to care, and treatment among persons of color and transgender/nonbinary identity, socioeconomically disadvantaged populations, and adolescents, and a disturbing low awareness among health care providers, policymakers, and the lay public about endometriosis, which, if left undiagnosed and under-treated can lead to significant fibrosis, infertility, depression, and markedly diminished quality of life. This review summarizes endometriosis epidemiology, compelling evidence for its pathogenesis, mechanisms underlying its pathophysiology in the age of precision medicine, recent biomarker discovery, novel therapeutic approaches, and issues around reproductive justice for marginalized populations with this disorder spanning the past 100 years. As we enter the next revolution in health care and biomedical research, with rich molecular and clinical datasets, single-cell omics, and population-level data, endometriosis is well positioned to benefit from data-driven research leveraging computational and artificial intelligence approaches integrating data and predicting disease risk, diagnosis, response to medical and surgical therapies, and prognosis for recurrence.

Preconception Genetic Carrier Screening for Miscarriage Risk Assessment: A Bioinformatic Approach to Identifying Candidate Lethal Genes and Variants.

Purpose: Miscarriage, due to genetically heterogeneous etiology, is a common outcome of pregnancy. Preconception genetic carrier screening (PGCS) identifies at-risk partners for newborn genetic disorders; however, PGCS panels currently lack miscarriage-related genes. Here we assessed the theoretical impact of known and candidate genes on prenatal lethality and the PGCS among diverse populations. Methods: Human exome sequencing and mouse gene function databases were analyzed to define genes essential for human fetal survival (lethal genes), identify variants that are absent in a homozygous state in healthy human population, and to estimate carrier rates for known and candidate lethal genes. Results: Among 138 genes, potential lethal variants are present in the general population with a frequency of 0.5% or greater. Preconception screening for these 138 genes would identify from 4.6% (Finnish population) to 39.8% (East Asian population) of couples that are at-risk for miscarriage, explaining a cause for pregnancy loss for ∼1.1-10% of conceptions affected by biallelic lethal variants. Conclusion: This study identified a set of genes and variants potentially associated with lethality across different ethnic backgrounds. The diversity of these genes amongst the various ethnic groups highlights the importance of designing a pan-ethnic PGCS panel comprising miscarriage-related genes.