pathfindR: An R Package for Comprehensive Identification of Enriched Pathways in Omics Data Through Active Subnetworks

The type of DNA damage response after decitabine treatment depends on the level of DNMT activity

Decitabine and azacytidine are considered as epigenetic drugs that induce DNA methyltransferase (DNMT)-DNA crosslinks, resulting in DNA hypomethylation and damage. Although they are already applied against myeloid cancers, important aspects of their mode of action remain unknown, highly limiting their clinical potential. Using a combinatorial approach, we reveal that the efficacy profile of both compounds primarily depends on the level of induced DNA damage. Under low DNMT activity, only decitabine has a substantial impact. Conversely, when DNMT activity is high, toxicity and cellular response to both compounds are dramatically increased, but do not primarily depend on DNA hypomethylation or RNA-associated processes. By investigating proteome dynamics on chromatin, we show that decitabine induces a strictly DNMT-dependent multifaceted DNA damage response based on chromatin recruitment, but not expression-level changes of repair-associated proteins. The choice of DNA repair pathway hereby depends on the severity of decitabine-induced DNA lesions. Although under moderate DNMT activity, mismatch (MMR), base excision (BER), and Fanconi anaemia-dependent DNA repair combined with homologous recombination are activated in response to decitabine, high DNMT activity and therefore immense replication stress induce activation of MMR and BER followed by non-homologous end joining.

Increased [18F]FDG uptake in the infarcted myocardial area displayed by combined PET/CMR correlates with snRNA-seq-detected inflammatory cell invasion

Combined [18F]FDG PET-cardiac MRI imaging (PET/CMR) is a useful tool to assess myocardial viability and cardiac function in patients with acute myocardial infarction (AMI). Here, we evaluated the prognostic value of PET/CMR in a porcine closed-chest reperfused AMI (rAMI) model. Late gadolinium enhancement by PET/CMR imaging displayed tracer uptake defect at the infarction site by 3 days after the rAMI in the majority of the animals (group Match, n = 28). Increased [18F]FDG uptake at the infarcted area (metabolism/contractility mismatch) with reduced tracer uptake in the remote viable myocardium (group Mismatch, n = 12) 3 days after rAMI was observed in the animals with larger infarct size and worse left ventricular ejection fraction (LVEF) (34 ± 8.7 vs 42.0 ± 5.2%), with lower LVEF also at the 1-month follow-up (35.8 ± 9.5 vs 43.0 ± 6.3%). Transcriptome analyses by bulk and single-nuclei RNA sequencing of the infarcted myocardium and border zones (n = 3 of each group, and 3 sham-operated controls) revealed a strong inflammatory response with infiltration of monocytes and macrophages in the infarcted and border areas in Mismatch animals. Our data indicate a high prognostic relevance of combined PET/MRI in the subacute phase of rAMI for subsequent impairment of heart function and underline the adverse effects of an excessive activation of the innate immune system in the initial phase after rAMI.

The clock gene BHLHE40 and atypical CCNG2 control androgen-induced cellular senescence as a novel tumor suppressive pathway in prostate cancer

BACKGROUND

The androgen receptor (AR) is a drug target used to inhibit AR and prostate cancer (PCa) growth. Surprisingly, treatment with supraphysiological androgen level (SAL), used in bipolar androgen therapy, inhibits growth of PCa suggesting a tumor-suppressive activity by SAL. SAL was shown to induce cellular senescence in PCa.

METHODS

RNA-seq and transcriptome analysis, ChIP-seq, human 3D PCa spheroids, mouse xenografted castration-resistant PCa, knockdown and overexpression, Co-immunoprecipitation (Co-IP), translocation analysis, immune detection, qRT-PCR, protein-protein interaction modelling.

RESULTS

Here, mice xenografts with castration-resistant PCa tumors show that SAL inhibits cancer growth in vivo suggesting that SAL activates a tumor-suppressive mechanism. RNA-seq and ChIP-seq revealed the clock gene BHLHE40 is a novel direct AR target. Compared to adjacent human prostate tissues, the expression of BHLHE40 is reduced in PCa tumors and associated with reduced survival. Knockdown suggests that BHLHE40 mediates SAL-induced cellular senescence including tumor spheroids. Interestingly, a large overlap of differentially expressed gene sets was identified between BHLHE40 and SAL leading to the identification of four classes of SAL-BHLHE40 transcriptome landscapes. Co-IP and modelling suggest binding of BHLHE40 to AR and their co-translocation into nucleus by SAL treatment. Further, RNA-seq and ChIP-seq analysis indicate that the atypical tumor suppressive cyclin G2 emerged as a novel downstream target of BHLHE40 and a mediator of SAL-induced cellular senescence.

CONCLUSIONS

The data provide evidence of the tumor suppressive activity of SAL and a novel signaling by the AR-BHLHE40-CCNG2 axis for androgen-induced cellular senescence, linking circadian rhythm factor to androgen signaling as a novel tumor suppressive pathway.

Transcriptome Study of 2 Black Cohorts Reveals cis Long Noncoding RNAs Associated With Hypertension-Related mRNAs

BACKGROUND

Long noncoding RNAs (lncRNAs) have emerged as critical regulators of the expression of genes involved in cardiovascular diseases. This project aims to identify circulating lncRNAs associated with protein-coding mRNAs differentially expressed between hypertensive and normotensive individuals and establish their link with hypertension.

METHODS AND RESULTS

The analyses were conducted in 3 main steps: (1) an unbiased whole blood transcriptome-wide analysis was conducted to identify and replicate protein-coding genes differentially expressed by hypertension status in 497 and 179 Black individuals from the GENE-FORECAST (Genomics, Environmental Factors and the Social Determinants of Cardiovascular Disease in African-Americans Study) and MH-GRID (Minority Health Genomics and Translational Research Bio-Repository Database) studies, respectively. Subsequently, (2) proximal lncRNAs, termed cis lncRNA quantitative trait loci, associated with each mRNA were identified in the GENE-FORECAST study and replicated in the MH-GRID study. Finally, (3) the lncRNA quantitative trait loci were used as predictors in a random forest model to predict hypertension in both data sets. A total of 129 mRNAs were significantly differentially expressed between normotensive and hypertensive individuals in both data sets. The lncRNA-mRNA association analysis revealed 249 cis lncRNA quantitative trait loci associated with 102 mRNAs, including VAMP2 (vesicle-associated membrane protein 2), mitogen-activated protein kinase kinase 3, CCAAT enhancer binding protein beta, and lymphocyte antigen 6 complex, locus E. The 249 lncRNA quantitative trait loci predicted hypertension with an area under the curve of 0.79 and 0.71 in GENE-FORECAST and MH-GRID studies, respectively.

CONCLUSIONS

This study leveraged a significant sample of Black individuals, a population facing a disproportionate burden of hypertension. The analyses unveiled a total of 271 lncRNA-mRNA relationships involving mRNAs that play critical roles in vascular pathways relevant to blood pressure regulation. The compelling findings, consistent across 2 independent data sets, establish a reliable foundation for designing in vitro/in vivo experiments.

MYCT1 controls environmental sensing in human haematopoietic stem cells

The processes that govern human haematopoietic stem cell (HSC) self-renewal and engraftment are poorly understood and challenging to recapitulate in culture to reliably expand functional HSCs1-3. Here we identify MYC target 1 (MYCT1; also known as MTLC) as a crucial human HSC regulator that moderates endocytosis and environmental sensing in HSCs. MYCT1 is selectively expressed in undifferentiated human haematopoietic stem and progenitor cells (HSPCs) and endothelial cells but becomes markedly downregulated during HSC culture. Lentivirus-mediated knockdown of MYCT1 prevented human fetal liver and cord blood (CB) HSPC expansion and engraftment. By contrast, restoring MYCT1 expression improved the expansion and engraftment of cultured CB HSPCs. Single-cell RNA sequencing of human CB HSPCs in which MYCT1 was knocked down or overexpressed revealed that MYCT1 governs important regulatory programmes and cellular properties essential for HSC stemness, such as ETS factor expression and low mitochondrial activity. MYCT1 is localized in the endosomal membrane in HSPCs and interacts with vesicle trafficking regulators and signalling machinery. MYCT1 loss in HSPCs led to excessive endocytosis and hyperactive signalling responses, whereas restoring MYCT1 expression balanced culture-induced endocytosis and dysregulated signalling. Moreover, sorting cultured CB HSPCs on the basis of lowest endocytosis rate identified HSPCs with preserved MYCT1 expression and MYCT1-regulated HSC stemness programmes. Our work identifies MYCT1-moderated endocytosis and environmental sensing as essential regulatory mechanisms required to preserve human HSC stemness. Our data also pinpoint silencing of MYCT1 as a cell-culture-induced vulnerability that compromises human HSC expansion.

Rewiring Saccharomyces cerevisiae metabolism for optimised Taxol® precursors production

Saccharomyces cerevisiae has been conveniently used to produce Taxol® anticancer drug early precursors. However, the harmful impact of oxidative stress by the first cytochrome P450-reductase enzymes (CYP725A4-POR) of Taxol® pathway has hampered sufficient progress in yeast. Here, we evolved an oxidative stress-resistant yeast strain with three-fold higher titre of their substrate, taxadiene. The performance of the evolved and parent strains were then evaluated in galactose-limited chemostats before and under the oxidative stress by an oxidising agent. The interaction of evolution and oxidative stress was comprehensively evaluated through transcriptomics and metabolite profiles integration in yeast enzyme-constrained genome scale model. Overall, the evolved strain showed improved respiration, reduced overflow metabolites production and oxidative stress re-induction tolerance. The cross-protection mechanism also potentially contributed to better heme, flavin and NADPH availability, essential for CYP725A4 and POR optimal activity in yeast. The results imply that the evolved strain is a robust cell factory for future efforts towards Taxol© production.

AICAR confers prophylactic cardioprotection in doxorubicin-induced heart failure in rats

Doxorubicin (DOX) is a widely used chemotherapeutic agent that can cause serious cardiotoxic side effects, leading to heart failure (HF). Impaired mitochondrial function is thought to be key factor driving progression into HF. We have previously shown in a rat model of DOX-HF that heart failure with reduced ejection fraction correlates with mitochondrial loss and dysfunction. Adenosine monophosphate-dependent kinase (AMPK) is a cellular energy sensor, regulating mitochondrial biogenesis and energy metabolism, including fatty acid oxidation. We hypothesised that AMPK activation could restore mitochondrial function and therefore be a novel cardioprotective strategy for the prevention of DOX-HF. Consequently, we set out to assess whether 5-aminoimidazole-4-carboxamide 1-β-D-ribofuranoside (AICAR), an activator of AMPK, could prevent cardiac functional decline in this chronic intravenous rat model of DOX-HF. In line with our hypothesis, AICAR improved cardiac systolic function. AICAR furthermore improved cardiac mitochondrial fatty acid oxidation, independent of mitochondrial number, and in the absence of observable AMPK-activation. In addition, we found that AICAR prevented loss of myocardial mass. RNAseq analysis showed that this may be driven by normalisation of pathways associated with ribosome function and protein synthesis, which are impaired in DOX-treated rat hearts. AICAR furthermore prevented dyslipidemia and excessive body-weight loss in DOX-treated rats, which may contribute to preservation of myocardial mass. Though it is unclear whether AICAR exerted its cardioprotective effect through cardiac or extra-cardiac AMPK-activation or via an AMPK-independent effect, these results show promise for the use of AICAR as a cardioprotective agent in DOX-HF to both preserve cardiac function and mass.

Novel Multiplexed Plasma Biomarker Panel Has Diagnostic and Prognostic Potential in Children With Hypertrophic Cardiomyopathy

BACKGROUND

Hypertrophic cardiomyopathy (HCM) is defined clinically by pathological left ventricular hypertrophy. We have previously developed a plasma proteomics biomarker panel that correlates with clinical markers of disease severity and sudden cardiac death risk in adult patients with HCM. The aim of this study was to investigate the utility of adult biomarkers and perform new discoveries in proteomics for childhood-onset HCM.

METHODS

Fifty-nine protein biomarkers were identified from an exploratory plasma proteomics screen in children with HCM and augmented into our existing multiplexed targeted liquid chromatography-tandem/mass spectrometry-based assay. The association of these biomarkers with clinical phenotypes and outcomes was prospectively tested in plasma collected from 148 children with HCM and 50 healthy controls. Machine learning techniques were used to develop novel pediatric plasma proteomic biomarker panels.

RESULTS

Four previously identified adult HCM markers (aldolase fructose-bisphosphate A, complement C3a, talin-1, and thrombospondin 1) and 3 new markers (glycogen phosphorylase B, lipoprotein a and profilin 1) were elevated in pediatric HCM. Using supervised machine learning applied to training (n=137) and validation cohorts (n=61), this 7-biomarker panel differentiated HCM from healthy controls with an area under the curve of 1.0 in the training data set (sensitivity 100% [95% CI, 95-100]; specificity 100% [95% CI, 96-100]) and 0.82 in the validation data set (sensitivity 75% [95% CI, 59-86]; specificity 88% [95% CI, 75-94]). Reduced circulating levels of 4 other peptides (apolipoprotein L1, complement 5b, immunoglobulin heavy constant epsilon, and serum amyloid A4) found in children with high sudden cardiac death risk provided complete separation from the low and intermediate risk groups and predicted mortality and adverse arrhythmic outcomes (hazard ratio, 2.04 [95% CI, 1.0-4.2]; P=0.044).

CONCLUSIONS

In children, a 7-biomarker proteomics panel can distinguish HCM from controls with high sensitivity and specificity, and another 4-biomarker panel identifies those at high risk of adverse arrhythmic outcomes, including sudden cardiac death.

Vaginal microbes alter epithelial transcriptome and induce epigenomic modifications providing insight into mechanisms for susceptibility to adverse reproductive outcomes

The cervicovaginal microbiome is highly associated with women's health, with microbial communities dominated by Lactobacillus species considered optimal. Conversely, a lack of lactobacilli and a high abundance of strict and facultative anaerobes, including Gardnerella vaginalis, have been associated with adverse reproductive outcomes. However, how host-microbial interactions alter specific molecular pathways and impact cervical and vaginal epithelial function remains unclear. Using RNA-sequencing, we characterized the in vitro cervicovaginal epithelial transcriptional response to different vaginal bacteria and their culture supernatants. We showed that G. vaginalis upregulates genes associated with an activated innate immune response. Unexpectedly, G. vaginalis specifically induced inflammasome pathways through activation of NLRP3-mediated increases in caspase-1, IL-1β and cell death, while live L. crispatus had minimal transcriptomic changes on epithelial cells. L. crispatus culture supernatants resulted in a shift in the epigenomic landscape of cervical epithelial cells that was confirmed by ATAC-sequencing showing reduced chromatin accessibility. This study reveals new insights into host-microbe interactions in the lower reproductive tract and suggests potential therapeutic strategies leveraging the vaginal microbiome to improve reproductive health.

Targeting neddylation and sumoylation in chemoresistant triple negative breast cancer

Triple negative breast cancer (TNBC) accounts for 15-20% of breast cancer cases in the United States. Systemic neoadjuvant chemotherapy (NACT), with or without immunotherapy, is the current standard of care for patients with early-stage TNBC. However, up to 70% of TNBC patients have significant residual disease once NACT is completed, which is associated with a high risk of developing recurrence within two to three years of surgical resection. To identify targetable vulnerabilities in chemoresistant TNBC, we generated longitudinal patient-derived xenograft (PDX) models from TNBC tumors before and after patients received NACT. We then compiled transcriptomes and drug response profiles for all models. Transcriptomic analysis identified the enrichment of aberrant protein homeostasis pathways in models from post-NACT tumors relative to pre-NACT tumors. This observation correlated with increased sensitivity in vitro to inhibitors targeting the proteasome, heat shock proteins, and neddylation pathways. Pevonedistat, a drug annotated as a NEDD8-activating enzyme (NAE) inhibitor, was prioritized for validation in vivo and demonstrated efficacy as a single agent in multiple PDX models of TNBC. Pharmacotranscriptomic analysis identified a pathway-level correlation between pevonedistat activity and post-translational modification (PTM) machinery, particularly involving neddylation and sumoylation targets. Elevated levels of both NEDD8 and SUMO1 were observed in models exhibiting a favorable response to pevonedistat compared to those with a less favorable response in vivo. Moreover, a correlation emerged between the expression of neddylation-regulated pathways and tumor response to pevonedistat, indicating that targeting these PTM pathways may prove effective in combating chemoresistant TNBC.

Neocarzilin Inhibits Cancer Cell Proliferation via BST-2 Degradation, Resulting in Lipid Raft-Trapped EGFR

Neocarzilin (NCA) is a natural product exhibiting potent antimigratory as well as antiproliferative effects. While vesicle amine transport protein 1 (VAT-1) was previously shown to inhibit migration upon NCA binding, the molecular mechanisms responsible for impaired proliferation remained elusive. We here introduce a chemical probe closely resembling the structural and stereochemical features of NCA and unravel bone marrow stromal antigen 2 (BST-2) as one of the targets responsible for the antiproliferative effect of NCA in cancer cells. The antiproliferative mechanism of NCA was confirmed in corresponding BST-2 knockout (KO) HeLa cells, which were less sensitive to compound treatment. Vice versa, reconstitution of BST-2 in the KO cells again reduced proliferation upon NCA addition, comparable to that of wild-type (wt) HeLa cells. Whole proteome mass spectrometric (MS) analysis of NCA-treated wt and KO cancer cells revealed regulated pathways and showed reduced levels of BST-2 upon NCA treatment. In-depth analysis of BST-2 levels in response to proteasome and lysosome inhibitors unraveled a lysosomal degradation path upon NCA treatment. As BST-2 mediates the release of epidermal growth factor receptor (EGFR) from lipid rafts to turn on proliferation signaling pathways, reduced BST-2 levels led to attenuated phosphorylation of STAT3. Furthermore, fluorescence microscopy confirmed increased colocalization of EGFR and lipid rafts in the presence of NCA. Overall, NCA represents a versatile anticancer natural product with a unique dual mode of action and unconventional inhibition of proliferation via BST-2 degradation.

Multiregional transcriptomics identifies congruent consensus subtypes with prognostic value beyond tumor heterogeneity of colorectal cancer

Intra-tumor heterogeneity compromises the clinical value of transcriptomic classifications of colorectal cancer. We investigated the prognostic effect of transcriptomic heterogeneity and the potential for classifications less vulnerable to heterogeneity in a single-hospital series of 1093 tumor samples from 692 patients, including multiregional samples from 98 primary tumors and 35 primary-metastasis sets. We show that intra-tumor heterogeneity of the consensus molecular subtypes (CMS) is frequent and has poor-prognostic associations independently of tumor microenvironment markers. Multiregional transcriptomics uncover cancer cell-intrinsic and low-heterogeneity signals that recapitulate the intrinsic CMSs proposed by single-cell sequencing. Further subclassification identifies congruent CMSs that explain a larger proportion of variation in patient survival than intra-tumor heterogeneity. Plasticity is indicated by discordant intrinsic phenotypes of matched primary and metastatic tumors. We conclude that multiregional sampling reconciles the prognostic power of tumor classifications from single-cell and bulk transcriptomics in the context of intra-tumor heterogeneity, and phenotypic plasticity challenges the reconciliation of primary and metastatic subtypes.

Genes Differentially Expressed Across Major Arteries Are Enriched in Endothelial Dysfunction-Related Gene Sets: Implications for Relative Inter-artery Atherosclerosis Risk

Atherosclerosis differs across major arteries. Although the biological basis is not fully understood, limited evidence of genetic differences has been documented. This study, therefore, was aimed to identify differentially expressed genes between clinically relevant major arteries and investigate their enrichment in endothelial dysfunction-related gene sets. A bioinformatic analysis of publicly available gene-level read counts for coronary, aortic, and tibial arteries was performed. Differential gene expression was conducted with DeSeq2 at a false discovery rate of 0.05. Differentially expressed genes were then subjected to over-representation analysis and active-subnetwork-oriented enrichment analysis, both at a false discovery rate of 0.005. Enriched terms common to both analyses were categorized for each contrast into immunity/inflammation-, membrane biology-, lipid metabolism-, and coagulation-related terms, and the top differentially expressed genes validated against Swiss Institute of Bioinformatics' Bgee database. There was mostly upregulation of differentially expressed genes for the coronary/tibial and aorta/tibial contrasts, but milder changes for the coronary/aorta contrast. Transcriptomic differences between coronary or aortic versus tibial samples largely involved immunity/inflammation-, membrane biology-, lipid metabolism-, and coagulation-related genes, suggesting potential to modulate endothelial dysfunction and atherosclerosis. These results imply atheroprone coronary and aortic environments compared with tibial artery tissue, which may explain observed relative inter-artery atherosclerosis risk.

Genetic and transcriptomic landscape of colonic diverticulosis

OBJECTIVE

Colonic diverticulosis is a prevalent condition among older adults, marked by the presence of thin-walled pockets in the colon wall that can become inflamed, infected, haemorrhage or rupture. We present a case-control genetic and transcriptomic study aimed at identifying the genetic and cellular determinants underlying this condition and the relationship with other gastrointestinal disorders.

DESIGN

We conducted DNA and RNA sequencing on colonic tissue from 404 patients with (N=172) and without (N=232) diverticulosis. We investigated variation in the transcriptome associated with diverticulosis and further integrated this variation with single-cell RNA-seq data from the human intestine. We also integrated our expression quantitative trait loci with genome-wide association study using Mendelian randomisation (MR). Furthermore, a Polygenic Risk Score analysis gauged associations between diverticulosis severity and other gastrointestinal disorders.

RESULTS

We discerned 38 genes with differential expression and 17 with varied transcript usage linked to diverticulosis, indicating tissue remodelling as a primary diverticula formation mechanism. Diverticula formation was primarily linked to stromal and epithelial cells in the colon including endothelial cells, myofibroblasts, fibroblasts, goblet, tuft, enterocytes, neurons and glia. MR highlighted five genes including CCN3, CRISPLD2, ENTPD7, PHGR1 and TNFSF13, with potential causal effects on diverticulosis. Notably, ENTPD7 upregulation was confirmed in diverticulosis cases. Additionally, diverticulosis severity was positively correlated with genetic predisposition to diverticulitis.

CONCLUSION

Our results suggest that tissue remodelling is a primary mechanism for diverticula formation. Individuals with an increased genetic proclivity to diverticulitis exhibit a larger numbers of diverticula on colonoscopy.

The Long Non-Coding RNA Obesity-Related (Obr) Contributes To Lipid Metabolism Through Epigenetic Regulation

Obesity is a multifactorial disease that is part of today's epidemic and also increases the risk of other metabolic diseases. Long noncoding RNAs (lncRNAs) provide one tier of regulatory mechanisms to maintain metabolic homeostasis. Although lncRNAs are a significant constituent of the mammalian genome, studies aimed at their metabolic significance, including obesity, are only beginning to be addressed. Here, a developmentally regulated lncRNA, termed as obesity related (Obr), whose expression in metabolically relevant tissues such as skeletal muscle, liver, and pancreas is altered in diet-induced obesity, is identified. The Clone 9 cell line and high-fat diet-induced obese Wistar rats are used as a model system to verify the function of Obr. By using stable expression and antisense oligonucleotide-mediated downregulation of the expression of Obr followed by different molecular biology experiments, its role in lipid metabolism is verified. It is shown that Obr associates with the cAMP response element-binding protein (Creb) and activates different transcription factors involved in lipid metabolism. Its association with the Creb histone acetyltransferase complex, which includes the cAMP response element-binding protein (CBP) and p300, positively regulates the transcription of genes involved in lipid metabolism. In addition, Obr is regulated by Pparγ in response to lipid accumulation.

Osteopontin drives neuroinflammation and cell loss in MAPT-N279K frontotemporal dementia patient neurons

Frontotemporal dementia (FTD) is an incurable group of early-onset dementias that can be caused by the deposition of hyperphosphorylated tau in patient brains. However, the mechanisms leading to neurodegeneration remain largely unknown. Here, we combined single-cell analyses of FTD patient brains with a stem cell culture and transplantation model of FTD. We identified disease phenotypes in FTD neurons carrying the MAPT-N279K mutation, which were related to oxidative stress, oxidative phosphorylation, and neuroinflammation with an upregulation of the inflammation-associated protein osteopontin (OPN). Human FTD neurons survived less and elicited an increased microglial response after transplantation into the mouse forebrain, which we further characterized by single nucleus RNA sequencing of microdissected grafts. Notably, downregulation of OPN in engrafted FTD neurons resulted in improved engraftment and reduced microglial infiltration, indicating an immune-modulatory role of OPN in patient neurons, which may represent a potential therapeutic target in FTD.

Inhibition of Crif1 protects fatty acid-induced POMC neuron-like cell-line damage by increasing CPT-1 function

Hypothalamic proopiomelanocortin (POMC) neurons are sensors of signals that reflect the energy stored in the body. Inducing mild stress in proopiomelanocortin neurons protects them from the damage promoted by the consumption of a high-fat diet, mitigating the development of obesity; however, the cellular mechanisms behind these effects are unknown. Here, we induced mild stress in a proopiomelanocortin neuron cell line by inhibiting Crif1. In proopiomelanocortin neurons exposed to high levels of palmitate, the partial inhibition of Crif1 reverted the defects in mitochondrial respiration and ATP production; this was accompanied by improved mitochondrial fusion/fission cycling. Furthermore, the partial inhibition of Crif1 resulted in increased reactive oxygen species production, increased fatty acid oxidation, and reduced dependency on glucose for mitochondrial respiration. These changes were dependent on the activity of CPT-1. Thus, we identified a CPT-1-dependent metabolic shift toward greater utilization of fatty acids as substrates for respiration as the mechanism behind the protective effect of mild stress against palmitate-induced damage of proopiomelanocortin neurons.NEW & NOTEWORTHY Saturated fats can damage hypothalamic neurons resulting in positive energy balance, and this is mitigated by mild cellular stress; however, the mechanisms behind this protective effect are unknown. Using a proopiomelanocortin cell line, we show that under exposure to a high concentration of palmitate, the partial inhibition of the mitochondrial protein Crif1 results in protection due to a metabolic shift warranted by the increased expression and activity of the mitochondrial fatty acid transporter CPT-1.

Transcriptomic crosstalk between viral and host factors drives aberrant homeostasis of T-cell proliferation and cell death in HIV-infected immunological non-responders

BACKGROUND

Immunological non-responders (INRs) among people living with HIV have inherently higher mortality and morbidity rates. The underlying immunological mechanisms whereby failure of immune reconstitution occurs in INRs require elucidation.

METHOD

HIV-1 DNA and HIV-1 cell-associated RNA (CA-HIV RNA) quantifications were conducted via RT-qPCR. Transcriptome sequencing (RNA-seq), bioinformatics, and biological verifications were performed to discern the crosstalk between host and viral factors. Flow cytometry was employed to analyze cellular activation, proliferation, and death.

RESULTS

HIV-1 DNA and CA-HIV RNA levels were observed to be significantly higher in INRs compared to immunological responders (IRs). Evaluation of CD4/CD8 ratios showed a significantly negative correlation with HIV-1 DNA in IRs, but not in INRs. Bioinformatics analyses and biological verifications showed IRF7/INF-α regulated antiviral response was intensified in INRs. PBMCs of INRs expressed significantly more HIV integrase-mRNA (p31) than IRs. Resting (CD4+CD69- T-cells) and activated (CD4+CD69+ T-cells) HIV-1 reservoir harboring cells were significantly higher in INRs, with the co-occurrence of significantly higher cellular proliferation and cell death in CD4+ T-cells of INRs.

CONCLUSION

In INRs, the systematic crosstalk between the HIV-1 reservoir and host cells tends to maintain a persistent antiviral response-associated inflammatory environment, which drives aberrant cellular activation, proliferation, and death of CD4+ T-cells.

Data-Driven Characterization of Genetic Variability in Disease Pathways and Pesticide-Induced Nervous System Disease in the United States Population

BACKGROUND

Genetic susceptibility to chemicals is incompletely characterized. However, nervous system disease development following pesticide exposure can vary in a population, implying some individuals may have higher genetic susceptibility to pesticide-induced nervous system disease.

OBJECTIVES

We aimed to build a computational approach to characterize single-nucleotide polymorphisms (SNPs) implicated in chemically induced adverse outcomes and used this framework to assess the link between differential population susceptibility to pesticides and human nervous system disease.

METHODS

We integrated publicly available datasets of Chemical-Gene, Gene-Pathway, and SNP-Disease associations to build Chemical-Pathway-Gene-SNP-Disease linkages for humans. As a case study, we integrated these linkages with spatialized pesticide application data for the US from 1992 to 2018 and spatialized nervous system disease rates for 2018. Through this, we characterized SNPs that may be important in states with high disease occurrence based on the pesticides used there.

RESULTS

We found that the number of SNP hits per pesticide in US states positively correlated with disease incidence and prevalence for Alzheimer's disease, Parkinson disease, and multiple sclerosis. We performed frequent itemset mining to differentiate pesticides used over time in states with high and low disease occurrence and found that only 19% of pesticide sets overlapped between 10 states with high disease occurrence and 10 states with low disease occurrence rates, and more SNPs were implicated in pathways in high disease occurrence states. Through a cross-validation of subsets of five high and low disease occurrence states, we characterized SNPs, genes, pathways, and pesticides more frequently implicated in high disease occurrence states.

DISCUSSION

Our findings support that pesticides contribute to nervous system disease, and we developed priority lists of SNPs, pesticides, and pathways for further study. This data-driven approach can be adapted to other chemicals, diseases, and locations to characterize differential population susceptibility to chemical exposures. https://doi.org/10.1289/EHP14108.

Alterations in genes associated with cytosolic RNA sensing in whole blood are associated with coronary microvascular disease in SLE

Systemic lupus erythematosus (SLE) patients are 90% women and over three times more likely to die of cardiovascular disease than women in the general population. Chest pain with no obstructive cardiac disease is associated with coronary microvascular disease (CMD), where narrowing of the small blood vessels can lead to ischemia, and frequently reported by SLE patients. Using whole blood RNA samples, we asked whether gene signatures discriminate SLE patients with coronary microvascular dysfunction (CMD) on cardiac MRI (n=4) from those without (n=7) and whether any signaling pathway is linked to the underlying pathobiology of SLE CMD. RNA-seq analysis revealed 143 differentially expressed (DE) genes between the SLE and healthy control (HC) groups, with virus defense and interferon (IFN) signaling being the key pathways identified as enriched in SLE as expected. We next conducted a comparative analysis of genes differentially expressed in SLE-CMD and SLE-non-CMD relative to HC samples. Our analysis highlighted differences in IFN signaling, RNA sensing and ADP-ribosylation pathways between SLE-CMD and SLE-non-CMD. This is the first study to investigate possible gene signatures associating with CMD in SLE, and our data strongly suggests that distinct molecular mechanisms underly vascular changes in CMD and non-CMD involvement in SLE.

Pten regulates endocytic trafficking of cell adhesion and Wnt signaling molecules to pattern the retina

The retina is exquisitely patterned, with neuronal somata positioned at regular intervals to completely sample the visual field. Here, we show that phosphatase and tensin homolog (Pten) controls starburst amacrine cell spacing by modulating vesicular trafficking of cell adhesion molecules and Wnt proteins. Single-cell transcriptomics and double-mutant analyses revealed that Pten and Down syndrome cell adhesion molecule Dscam) are co-expressed and function additively to pattern starburst amacrine cell mosaics. Mechanistically, Pten loss accelerates the endocytic trafficking of DSCAM, FAT3, and MEGF10 off the cell membrane and into endocytic vesicles in amacrine cells. Accordingly, the vesicular proteome, a molecular signature of the cell of origin, is enriched in exocytosis, vesicle-mediated transport, and receptor internalization proteins in Pten conditional knockout (PtencKO) retinas. Wnt signaling molecules are also enriched in PtencKO retinal vesicles, and the genetic or pharmacological disruption of Wnt signaling phenocopies amacrine cell patterning defects. Pten thus controls vesicular trafficking of cell adhesion and signaling molecules to establish retinal amacrine cell mosaics.

The molecular signature of BCR::ABLP210 and BCR::ABLT315I in a Drosophila melanogaster chronic myeloid leukemia model

Chronic myeloid leukemia (CML) is a clonal hematopoietic stem cell disorder resulting from a balanced translocation leading to BCR::ABL1 oncogene with increased tyrosine kinase activity. Despite the advancements in the development of tyrosine kinase inhibitors (TKIs), the T315I gatekeeper point mutation in the BCR::ABL1 gene remains a challenge. We have previously reported in a Drosophila CML model an increased hemocyte count and disruption in sessile hemocyte patterns upon expression of BCR::ABL1p210 and BCR::ABL1T315I in the hemolymph. In this study, we performed RNA sequencing to determine if there is a distinct gene expression that distinguishes BCR::ABL1p210 and BCR::ABL1T315I. We identified six genes that were consistently upregulated in the fly CML model and validated in adult and pediatric CML patients and in a mouse cell line expressing BCR::ABL1T315I. This study provides a comprehensive analysis of gene signatures in BCR::ABL1p210 and BCR::ABL1T315I, laying the groundwork for targeted investigations into the role of these genes in CML pathogenesis.

Establishment, characterization, and biobanking of 36 pancreatic cancer organoids: prediction of metastasis in resectable pancreatic cancer

PURPOSE

Early dissemination of primary pancreatic ductal adenocarcinoma (PDAC) is the main cause of dismal prognosis as it highly limits possible treatment options. A number of PDAC patients experience distant metastasis even after treatment due to the metastatic clones. We aimed to demonstrate the molecular architecture of borderline resectable PDAC manifests cancer dissemination of PDAC.

METHODS

Here, 36 organoids isolated from primary tumor masses of PDAC patients with diverse metastatic statues are presented. Whole-exome sequencing and RNA sequencing were performed and drug responses to clinically relevant 18 compounds were assessed.

RESULTS

Our results revealed that borderline resectable PDAC organoids exhibited distinct patterns according to their metastatic potency highlighted by multiple genetic and transcriptional factors and strong variances in drug responses.

CONCLUSIONS

These data suggest that the presence of metastatic PDAC can be identified by integrating molecular compositions and drug responses of borderline resectable PDAC.

Plasma Proteomics Elucidated a Protein Signature in COVID-19 Patients with Comorbidities and Early-Diagnosis Biomarkers

Despite great scientific efforts, deep understanding of coronavirus-19 disease (COVID-19) immunopathology and clinical biomarkers remains a challenge. Pre-existing comorbidities increase the mortality rate and aggravate the exacerbated immune response against the severe acute respiratory syndrome coronavirus-2 (SARS-CoV-2) infection, which can result in more severe symptoms as well as long-COVID and post-COVID complications. In this study, we applied proteomics analysis of plasma samples from 28 patients with SARS-CoV-2, with and without pre-existing comorbidities, as well as their corresponding controls to determine the systemic protein changes caused by the SARS-CoV-2 infection. As a result, the protein signature shared amongst COVID-19 patients with comorbidities was revealed to be characterized by alterations in the coagulation and complement pathways, acute-phase response proteins, tissue damage and remodeling, as well as cholesterol metabolism. These altered proteins may play a relevant role in COVID-19 pathophysiology. Moreover, several novel potential biomarkers for early diagnosis of the SARS-CoV-2 infection were detected, such as increased levels of keratin K22E, extracellular matrix protein-1 (ECM1), and acute-phase response protein α-2-antiplasmin (A2AP). Importantly, elevated A2AP may contribute to persistent clotting complications associated with the long-COVID syndrome in patients with comorbidities. This study provides new insights into COVID-19 pathogenesis and proposes novel potential biomarkers for early diagnosis that could be facilitated for clinical application by further validation studies.

The protein composition of exosomes released by prostate cancer cells is distinctly regulated by androgen receptor-antagonists and -agonist to stimulate growth of target cells

BACKGROUND

Prostate cancer (PCa) is a prevalent malignancy in men worldwide, ranking as the second leading cause of cancer-related death in Western countries. Various PCa hormone therapies, such as androgen receptor (AR)-antagonists or supraphysiological androgen level (SAL) reduce cancer cell proliferation. However, treated cells may influence the growth of neighboring cells through secreted exosomes in the tumor microenvironment (TME). Here, the change of protein content of exosomes secreted from PCa cells through treatment with different AR-antagonists or SAL has been analyzed.

METHODS

Isolation of exosomes via ultracentrifugation of treated human PCa LNCaP cells with AR-agonist and various AR-antagonists; analysis of cellular senescence by detection of senescence associated beta galactosidase activity (SA β-Gal); Western blotting and immunofluorescence staining; Mass spectrometry (MS-spec) of exosomes and bioinformatic analyses to identify ligand-specific exosomal proteins. Growth assays to analyze influence of exosomes on non-treated cells.

RESULTS

MS-spec analysis identified ligand-specific proteins in exosomes. One thousand seventy proteins were up- and 52 proteins downregulated by SAL whereas enzalutamide upregulated 151 proteins and downregulated 42 exosomal proteins. The bioinformatic prediction indicates an up-regulation of pro-proliferative pathways. AR ligands augment hub factors in exosomes that include AKT1, CALM1, PAK2 and CTNND1. Accordingly, functional assays confirmed that the isolated exosomes from AR-ligand treated cells promote growth of untreated PCa cells.

CONCLUSION

The data suggest that the cargo of exosomes is controlled by AR-agonist and -antagonists and distinct among the AR-antagonists. Further, exosomes promote growth that might influence the TME. This finding sheds light into the complex interplay between AR signaling and exosome-mediated communication between PCa cells.

Longitudinal Transcriptomic Analysis of Human Cortical Spheroids Identifies Axonal Dysregulation in the Prenatal Brain as a Mediator of Genetic Risk for Schizophrenia

BACKGROUND

Schizophrenia (SCZ) has a known neurodevelopmental etiology, but limited access to human prenatal brain tissue hampers the investigation of basic disease mechanisms in early brain development. Here, we elucidate the molecular mechanisms contributing to SCZ risk in a disease-relevant model of the prenatal human brain.

METHODS

We generated induced pluripotent stem cell-derived organoids, termed human cortical spheroids (hCSs), from a large, genetically stratified sample of 14 SCZ cases and 14 age- and sex-matched controls. The hCSs were differentiated for 150 days, and comprehensive molecular characterization across 4 time points was carried out.

RESULTS

The transcriptional and cellular architecture of hCSs closely resembled that of fetal brain tissue at 10 to 24 postconception weeks, showing strongest spatial overlap with frontal regions of the cerebral cortex. A total of 3520 genes were differentially modulated between SCZ and control hCSs across organoid maturation, displaying a significant contribution of genetic loading, an overrepresentation of risk genes for autism spectrum disorder and SCZ, and the strongest enrichment for axonal processes in all hCS stages. The two axon guidance genes SEMA7A and SEMA5A, the first a promoter of synaptic functions and the second a repressor, were downregulated and upregulated, respectively, in SCZ hCSs. This expression pattern was confirmed at the protein level and replicated in a large postmortem sample.

CONCLUSIONS

Applying a disease-relevant model of the developing fetal brain, we identified consistent dysregulation of axonal genes as an early risk factor for SCZ, providing novel insights into the effects of genetic predisposition on the neurodevelopmental origins of the disorder.

Identification of Key Genes Associated with Tumor Microenvironment Infiltration and Survival in Gastric Adenocarcinoma via Bioinformatics Analysis

OBJECTIVE

Gastric carcinoma (GC) is the fifth most commonly diagnosed cancer and the third leading cause of cancer-related deaths globally. The tumor microenvironment plays a significant role in the pathogenesis, prognosis, and response to immunotherapy. However, the immune-related molecular mechanisms underlying GC remain elusive. Bioinformatics analysis of the gene expression of GC and paracancerous healthy tissues from the same patient was performed to identify the key genes and signaling pathways, as well as their correlation to the infiltration of the tumor microenvironment (TME) by various immune cells related to GC development.

METHODS

We employed GSE19826, a gene expression profile from the Gene Expression Omnibus (GEO), for our analysis. Functional enrichment analysis of Differentially Expressed Genes (DEGs) was conducted using the Gene Ontology and Kyoto Encyclopedia of Genes and Genomes database.

RESULTS

Cytoscape software facilitated the identification of nine hub DEGs, namely, FN1, COL1A1, COL1A2, THBS2, COL3A1, COL5A1, APOE, SPP1, and BGN. Various network analysis algorithms were applied to determine their high connectivity. Among these hub genes, FN1, COL1A2, THBS2, COL3A1, COL5A1, and BGN were found to be associated with a poor prognosis for GC patients. Subsequent analysis using the TIMER database revealed the infiltration status of the TME concerning the overexpression of these six genes. Specifically, the abovementioned genes demonstrated direct correlations with cancer-associated fibroblasts, M1 and M2 macrophages, myeloid-derived suppressor cells, and activated dendritic cells.

CONCLUSION

Our findings suggest that the identified hub genes, particularly BGN, FN1, COL1A2, THBS2, COL3A1, and COL5A1, play crucial roles in GC prognosis and TME cell infiltration. This comprehensive analysis enhances our understanding of the molecular mechanisms underlying GC development and may contribute to the identification of potential therapeutic targets and prognostic markers for GC patients.

A Porphyromonas gingivalis hypothetical protein controlled by the type I-C CRISPR-Cas system is a novel adhesin important in virulence

The ability of many human pathogens to infect requires their ability to adhere to the host surfaces as a first step in the process. Porphyromonas gingivalis, a keystone oral pathogen, uses adhesins to adhere to the surface of the gingival epithelium and other members of the oral microbiome. In a previous study, we identified several proteins potentially linked to virulence whose mRNA levels are regulated by CRISPR-Cas type I-C. Among those, PGN_1547 was highly upregulated in the CRISPR-Cas 3 mutant. PGN_1547 is annotated as a hypothetical protein. Employing homology searching, our data support that PGN_1547 resembles an auto-transporter adhesin of P. gingivalis based on containing the DUF2807 domain. To begin to characterize the function of PGN_1547, we found that a deletion mutant displayed a significant decrease in virulence using a Galleria mellonela model. Furthermore, this mutant was significantly impaired in forming biofilms and attaching to the macrophage-like cell THP-1. Luminex revealed that the PGN_1547 mutant elicited a less robust cytokine and chemokine response from THP-1 cells, and TLR2 predominantly sensed that recombinant PGN_1547. Taken together, these findings broaden our understanding of the toolbox of virulence factors possessed by P. gingivalis. Importantly, PGN_1547, a hypothetical protein, has homologs in another member of the order Bacteroidales whose function is unknown, and our results could shed light on the role of this family of proteins as auto-transport adhesins in this phylogenetic group.IMPORTANCEPeriodontal diseases are among humans' most common infections, and besides their effect on the oral cavity, they have been associated with systemic inflammatory conditions. Among members of the oral microbiome implicated in the development of periodontitis, Porphyromonas gingivalis is considered a keystone pathogen. We have identified a new adhesin that acts as a virulence factor, PGN_1547, which contains the DUF2807 domain, which belongs to the putative auto-transporter adhesin, head GIN domain family. Deletion of this gene lowers the virulence of P. gingivalis and impacts the ability of P. gingivalis to form biofilm and attach to host cells. Furthermore, the broad distribution of these receptors in the order Bacteroidales suggests their importance in colonization by this important group of organisms.

Investigation of the impact of bromodomain inhibition on cytoskeleton stability and contraction

BACKGROUND

Injury to contractile organs such as the heart, vasculature, urinary bladder and gut can stimulate a pathological response that results in loss of normal contractility. PDGF and TGFβ are among the most well studied initiators of the injury response and have been shown to induce aberrant contraction in mechanically active cells of hollow organs including smooth muscle cells (SMC) and fibroblasts. However, the mechanisms driving contractile alterations downstream of PDGF and TGFβ in SMC and fibroblasts are incompletely understood, limiting therapeutic interventions.

METHODS

To identify potential molecular targets, we have leveraged the analysis of publicly available data, comparing transcriptomic changes in mechanically active cells stimulated with PDGF and TGFβ. Additional Analysis of publicly available data sets were performed on SMC and fibroblasts treated in the presence or absence of the MYC inhibitor JQ1. Validation of in silico findings were performed with qPCR, immunoblots, and collagen gel contraction assays measure the effect of JQ1 on cytoskeleton associated genes, proteins and contractility in mechanically active cells. Likelihood ratio test and FDR adjusted p-values were used to determine significant differentially expressed genes. Student ttest were used to calculate statistical significance of qPCR and contractility analyses.

RESULTS

Comparing PDGF and TGFβ stimulated SMC and fibroblasts identified a shared molecular profile regulated by MYC and members of the AP-1 transcription factor complex. Additional in silico analysis revealed a unique set of cytoskeleton-associated genes that were sensitive to MYC inhibition with JQ1. In vitro validation demonstrated JQ1 was also able to attenuate TGFβ and PDGF induced changes to the cytoskeleton and contraction of smooth muscle cells and fibroblasts in vitro.

CONCLUSIONS

These findings identify MYC as a key driver of aberrant cytoskeletal and contractile changes in fibroblasts and SMC, and suggest that JQ1 could be used to restore normal contractile function in hollow organs.

Evaluating synergistic effects of metformin and simvastatin on ovarian cancer cells

BACKGROUND

Ovarian Cancer (OC) stands as the most lethal gynecological malignancy, presenting an urgent clinical challenge in the quest to improve response rates. One approach to address this challenge is through drug repurposing, exemplified by the investigation of metabolic-modulating drugs such as Metformin (MTF) and Simvastatin (SIM). This study aims to explore the molecular mechanisms contributing to the potential synergistic anti-cancer effects between MTF and SIM on ovarian cancer cells.

METHODS

We assessed the effects of the combination on the proliferation and viability of two cell lines OVCAR-3 and SKOV-3. IC50 concentrations of MTF and SIM were determined using a proliferation assay, followed by subtoxic concentrations to explore the potential synergistic effects on the viability of both cell lines. Transcriptomic analysis was conducted on OVCAR-3 treated cells, and the findings were validated by assessing the expression levels of differentially expressed genes (DEGs) through real-time PCR in both cell lines SK-OV-3 and OVCAR-3.

RESULTS

Cytotoxicity analysis guided the selection of treatment concentrations as such MTF 10 mM and SIM 5 μM. The combined treatment of MTF and SIM demonstrated a synergistic inhibition of proliferation and viability in both cell lines. In OVCAR-3, exclusive identification of 507 DEGs was seen in the combination arm. Upregulation of FOXO3, RhoA, and TNFα, along with downregulation of PIK3R1, SKP2, and ATP6V1D levels, was observed in OVCAR-3 treated cells. Real-time PCR validation confirmed the consistency of expression levels for the mentioned DEGs.

CONCLUSION

Our data strongly supports the presence of synergy between MTF and SIM in OC cells. The combination's effect is associated with the dysregulation of genes in the key regulators AMPK and mTOR alongside other interconnected pathways.

An integrative proteomics approach identifies tyrosine kinase KIT as a therapeutic target for SPINK1-positive prostate cancer

Elevated serine peptidase inhibitor, Kazal type 1 (SPINK1) levels in ∼10%-25% of prostate cancer (PCa) patients associate with aggressive phenotype, for which there are limited treatment choices and dismal clinical outcomes. Using an integrative proteomics approach involving label-free phosphoproteome and proteome profiling, we delineated the downstream signaling pathways involved in SPINK1-mediated tumorigenesis and identified tyrosine kinase KIT as highly enriched. Furthermore, high to moderate levels of KIT expression were detected in ∼85% of SPINK1-positive PCa specimens. We show KIT signaling orchestrates SPINK1-mediated oncogenesis, and treatment with KIT inhibitor reduces tumor growth and metastases in preclinical mice models. Mechanistically, KIT signaling modulates WNT/β-catenin pathway and confers stemness-related features in PCa. Notably, inhibiting KIT signaling led to restoration of AR/REST levels, forming a feedback loop enabling SPINK1 repression. Overall, we uncover the role of KIT signaling downstream of SPINK1 in maintaining lineage plasticity and provide distinct treatment modalities for advanced-stage SPINK1-positive patients.

Mesenchyme instructs growth while epithelium directs branching in the mouse mammary gland

The mammary gland is a unique organ that undergoes dynamic alterations throughout a female's reproductive life, making it an ideal model for developmental, stem cell and cancer biology research. Mammary gland development begins in utero and proceeds via a quiescent bud stage before the initial outgrowth and subsequent branching morphogenesis. How mammary epithelial cells transit from quiescence to an actively proliferating and branching tissue during embryogenesis and, importantly, how the branch pattern is determined remain largely unknown. Here, we provide evidence indicating that epithelial cell proliferation and onset of branching are independent processes, yet partially coordinated by the Eda signaling pathway. Through heterotypic and heterochronic epithelial-mesenchymal recombination experiments between mouse mammary and salivary gland tissues and ex vivo live imaging, we demonstrate that unlike previously concluded, the mode of branching is an intrinsic property of the mammary epithelium whereas the pace of growth and the density of ductal tree are determined by the mesenchyme. Transcriptomic profiling and ex vivo and in vivo functional studies in mice disclose that mesenchymal Wnt/ß-catenin signaling, and in particular IGF-1 downstream of it critically regulate mammary gland growth. These results underscore the general need to carefully deconstruct the different developmental processes producing branched organs.

Oxidative stress governs mosquito innate immune signalling to reduce chikungunya virus infection in Aedes-derived cells

Arboviruses such as chikungunya, dengue and zika viruses cause debilitating diseases in humans. The principal vector species that transmits these viruses is the Aedes mosquito. Lack of substantial knowledge of the vector species hinders the advancement of strategies for controlling the spread of arboviruses. To supplement our information on mosquitoes' responses to virus infection, we utilized Aedes aegypti-derived Aag2 cells to study changes at the transcriptional level during infection with chikungunya virus (CHIKV). We observed that genes belonging to the redox pathway were significantly differentially regulated. Upon quantifying reactive oxygen species (ROS) in the cells during viral infection, we further discovered that ROS levels are considerably higher during the early hours of infection; however, as the infection progresses, an increase in antioxidant gene expression suppresses the oxidative stress in cells. Our study also suggests that ROS is a critical regulator of viral replication in cells and inhibits intracellular and extracellular viral replication by promoting the Rel2-mediated Imd immune signalling pathway. In conclusion, our study provides evidence for a regulatory role of oxidative stress in infected Aedes-derived cells.

Nicotinamide riboside activates renal metabolism and protects the kidney in a model of Alport syndrome

Chronic kidney disease (CKD) is associated with renal metabolic disturbances, including impaired fatty acid oxidation (FAO). Nicotinamide adenine dinucleotide (NAD + ) is a small molecule that participates in hundreds of metabolism-related reactions. NAD + levels are decreased in CKD, and NAD + supplementation is protective. However, both the mechanism of how NAD + supplementation protects from CKD, as well as the cell types most responsible, are poorly understood. Using a mouse model of Alport syndrome, we show that nicotinamide riboside (NR), an NAD + precursor, stimulates renal peroxisome proliferator-activated receptor α signaling and restores FAO in the proximal tubules, thereby protecting from CKD in both sexes. Bulk RNA-sequencing shows that renal metabolic pathways are impaired in Alport mice and dramatically activated by NR in both sexes. These transcriptional changes are confirmed by orthogonal imaging techniques and biochemical assays. Single nuclei RNA-sequencing and spatial transcriptomics, both the first of their kind from Alport mice, show that NAD + supplementation restores FAO in the proximal tubules with minimal effects on the podocytes. Finally, we also report, for the first time, sex differences at the transcriptional level in this Alport model. Male Alport mice had more severe inflammation and fibrosis than female mice at the transcriptional level. In summary, the data herein identify both the protective mechanism and location of NAD + supplementation in this model of CKD.

Genome-wide transcriptome profiling and development of age prediction models in the human brain

Aging-related transcriptome changes in various regions of the healthy human brain have been explored in previous works, however, a study to develop prediction models for age based on the expression levels of specific panels of transcripts is lacking. Moreover, studies that have assessed sexually dimorphic gene activities in the aging brain have reported discrepant results, suggesting that additional studies would be advantageous. The prefrontal cortex (PFC) region was previously shown to have a particularly large number of significant transcriptome alterations during healthy aging in a study that compared different regions in the human brain. We harmonized neuropathologically normal PFC transcriptome datasets obtained from the Gene Expression Omnibus (GEO) repository, ranging in age from 21 to 105 years, and found a large number of differentially regulated transcripts in the old and elderly, compared to young samples overall, and compared female and male-specific expression alterations. We assessed the genes that were associated with age by employing ontology, pathway, and network analyses. Furthermore, we applied various established (least absolute shrinkage and selection operator (Lasso) and Elastic Net (EN)) and recent (eXtreme Gradient Boosting (XGBoost) and Light Gradient Boosting Machine (LightGBM)) machine learning algorithms to develop accurate prediction models for chronological age and validated them. Studies to further validate these models in other large populations and molecular studies to elucidate the potential mechanisms by which the transcripts identified may be related to aging phenotypes would be advantageous.

Homeostatic status of thyroid hormones and brain water movement as determinant factors in biology of cerebral gliomas: a pilot study using a bioinformatics approach

Introduction

The expression and localization of the water channel transporters, aquaporins (AQPs), in the brain are substantially modified in gliomas during tumorigenesis, cell migration, edema formation, and resolution. We hypothesized that the molecular changes associated with AQP1 and AQP4 in the brain may potentially be anticancer therapeutic targets. To test this hypothesis, a bioinformatics analysis of publicly available data from international consortia was performed.

Methods

We used RNA-seq as an experimental strategy and identified the number of differential AQP1 and AQP4 transcript expressions in glioma tissue compared to normal brain tissue.

Results

AQPs genes are overexpressed in patients with glioma. Among the glioma subtypes, AQP1 and AQP4 were overexpressed in astrocytoma (low-grade glioma) and classical (high-grade glioma). Overall survival analysis demonstrated that both AQP genes can be used as prognostic factors for patients with low-grade glioma. Additionally, we observed a correlation between the expression of genes involved in the tyrosine and thyroid hormone pathways and AQPs, namely: PNMT, ALDH1A3, AOC2, HGDATP1B1, ADCY5, PLCB4, ITPR1, ATP1A3, LRP2, HDAC1, MED24, MTOR, and ACTB1 (Spearman's coefficient = geq 0.20 and p-value = ≤ 0.05).

Conclusion

Our findings indicate that the thyroid hormone pathways and AQPs 1 and 4 are potential targets for new anti-tumor drugs and therapeutic biomarkers for malignant gliomas.

C9orf72 gene networks in the human brain correlate with cortical thickness in C9-FTD and implicate vulnerable cell types

Introduction

A hexanucleotide repeat expansion (HRE) intronic to chromosome 9 open reading frame 72 (C9orf72) is recognized as the most common genetic cause of amyotrophic lateral sclerosis (ALS), frontotemporal dementia (FTD), and ALS-FTD. Identifying genes that show similar regional co-expression patterns to C9orf72 may help identify novel gene targets and biological mechanisms that mediate selective vulnerability to ALS and FTD pathogenesis.

Methods

We leveraged mRNA expression data in healthy brain from the Allen Human Brain Atlas to evaluate C9orf72 co-expression patterns. To do this, we correlated average C9orf72 expression values in 51 regions across different anatomical divisions (cortex, subcortex, and cerebellum) with average gene expression values for 15,633 protein-coding genes, including 54 genes known to be associated with ALS, FTD, or ALS-FTD. We then performed imaging transcriptomic analyses to evaluate whether the identified C9orf72 co-expressed genes correlated with patterns of cortical thickness in symptomatic C9orf72 pathogenic HRE carriers (n = 19) compared to controls (n = 23). Lastly, we explored whether genes with significant C9orf72 imaging transcriptomic correlations (i.e., "C9orf72 imaging transcriptomic network") were enriched in specific cell populations in the brain and enriched for specific biological and molecular pathways.

Results

A total of 2,120 genes showed an anatomical distribution of gene expression in the brain similar to C9orf72 and significantly correlated with patterns of cortical thickness in C9orf72 HRE carriers. This C9orf72 imaging transcriptomic network was differentially expressed in cell populations previously implicated in ALS and FTD, including layer 5b cells, cholinergic neurons in the spinal cord and brainstem and medium spiny neurons of the striatum, and was enriched for biological and molecular pathways associated with protein ubiquitination, autophagy, cellular response to DNA damage, endoplasmic reticulum to Golgi vesicle-mediated transport, among others.

Conclusion

Considered together, we identified a network of C9orf72 associated genes that may influence selective regional and cell-type-specific vulnerabilities in ALS/FTD.

Nicotinamide Mononucleotide (NMN) Works in Type 2 Diabetes through Unexpected Effects in Adipose Tissue, Not by Mitochondrial Biogenesis

Nicotinamide mononucleotide (NMN) has emerged as a promising therapeutic intervention for age-related disorders, including type 2 diabetes. In this study, we confirmed the previously observed effects of NMN treatment on glucose uptake and investigated its underlying mechanisms in various tissues and cell lines. Through the most comprehensive proteomic analysis to date, we discovered a series of novel organ-specific effects responsible for glucose uptake as measured by the IPGTT: adipose tissue growing (suggested by increased protein synthesis and degradation and mTOR proliferation signaling upregulation). Notably, we observed the upregulation of thermogenic UCP1, promoting enhanced glucose conversion to heat in intermuscular adipose tissue while showing a surprising repressive effect on mitochondrial biogenesis in muscle and the brain. Additionally, liver and muscle cells displayed a unique response, characterized by spliceosome downregulation and concurrent upregulation of chaperones, proteasomes, and ribosomes, leading to mildly impaired and energy-inefficient protein synthesis machinery. Furthermore, our findings revealed remarkable metabolic rewiring in the brain. This involved increased production of ketone bodies, downregulation of mitochondrial OXPHOS and TCA cycle components, as well as the induction of well-known fasting-associated effects. Collectively, our data elucidate the multifaceted nature of NMN action, highlighting its organ-specific effects and their role in improving glucose uptake. These findings deepen our understanding of NMN's therapeutic potential and pave the way for novel strategies in managing metabolic disorders.

The shifting lipidomic landscape of blood monocytes and neutrophils during pneumonia

The lipidome of immune cells during infection has remained unexplored, although evidence of the importance of lipids in the context of immunity is mounting. In this study, we performed untargeted lipidomic analysis of blood monocytes and neutrophils from patients hospitalized for pneumonia and age- and sex-matched noninfectious control volunteers. We annotated 521 and 706 lipids in monocytes and neutrophils, respectively, which were normalized to an extensive set of internal standards per lipid class. The cellular lipidomes were profoundly altered in patients, with both common and distinct changes between the cell types. Changes involved every level of the cellular lipidome: differential lipid species, class-wide shifts, and altered saturation patterns. Overall, differential lipids were mainly less abundant in monocytes and more abundant in neutrophils from patients. One month after hospital admission, lipidomic changes were fully resolved in monocytes and partially in neutrophils. Integration of lipidomic and concurrently collected transcriptomic data highlighted altered sphingolipid metabolism in both cell types. Inhibition of ceramide and sphingosine-1-phosphate synthesis in healthy monocytes and neutrophils resulted in blunted cytokine responses upon stimulation with lipopolysaccharide. These data reveal major lipidomic remodeling in immune cells during infection, and link the cellular lipidome to immune functionality.

Technical optimization of spatially resolved single-cell transcriptomic datasets to study clinical liver disease

Single cell and spatially resolved 'omic' techniques have enabled deep characterization of clinical pathologies that remain poorly understood, providing unprecedented insights into molecular mechanisms of disease. However, transcriptomic platforms are costly, limiting sample size, which increases the possibility of pre-analytical variables such as tissue processing and storage procedures impacting RNA quality and downstream analyses. Furthermore, spatial transcriptomics have not yet reached single cell resolution, leading to the development of multiple deconvolution methods to predict individual cell types within each transcriptome 'spot' on tissue sections. In this study, we performed spatial transcriptomics and single nucleus RNA sequencing (snRNAseq) on matched specimens from patients with either histologically normal or advanced fibrosis to establish important aspects of tissue handling, data processing, and downstream analyses of biobanked liver samples. We observed that tissue preservation technique impacts transcriptomic data, especially in fibrotic liver. Single cell mapping of the spatial transcriptome using paired snRNAseq data generated a spatially resolved, single cell dataset with 24 unique liver cell phenotypes. We determined that cell-cell interactions predicted using ligand-receptor analysis of snRNAseq data poorly correlated with cellular relationships identified using spatial transcriptomics. Our study provides a framework for generating spatially resolved, single cell datasets to study gene expression and cell-cell interactions in biobanked clinical samples with advanced liver disease.

Whole genome sequencing of a family with autosomal dominant features within the oculoauriculovertebral spectrum

Background

Oculoauriculovertebral Spectrum (OAVS) encompasses a wide variety of anomalies on derivatives from the first and second pharyngeal arches including macrostomia, hemifacial microsomia, micrognathia, preauricular tags, ocular and vertebral anomalies. We present the genetic findings of a large three-generation family with multiple members affected with macrostomia, preauricular tags and uni- or bilateral ptosis following an autosomal dominant segregation pattern.

Methods

We generated whole genome sequencing data for the proband, affected parent and unaffected paternal grandparent followed by Sanger sequencing on 23 family members for the top 10 candidate genes: KCND2, PDGFRA, CASP9, NCOA3, WNT10A, SIX1, MTF1, KDR/VEGFR2, LRRK1, and TRIM2. We performed parent and sibling-based transmission disequilibrium tests and burden analysis to explore segregation and burden of candidate gene mutations. Bioinformatic analyses investigated the biological connection between genes and the abnormal phenotypes.

Results

Overall, rare missense mutations in SIX1, KDR/VEGFR2, and PDGFRA showed the best evidence of segregation with the OAV phenotypes in this family. When considering affection with any of the 3 OAVS phenotypes as an outcome, parent-TDTs and sib-TDTs (unadjusted p-values) found that SIX1 (p=0.025, p=0.052), followed by PDGFRA (p=0.180, p=0.069) and KDR/VEGFR2 (p=0.180, p=0.069) have the strongest associations in this family. Burden analysis via a penalized linear mixed model identified SIX1 (RC=0.87) and PDGFRA (RC=0.98) as having the strongest association with OAVS severity. Using phenotype-specific ogfrautcomes, sib-TDTs identified associations between (1) SIX1 with uni- or bilateral ptosis (p=0.049) and ear tags (p=0.01), (2) PDGFRA and KDR/VEGFR2 with ear tags (both p<0.01).

Conclusion

Our study reports the genomic findings of a large family with multiple individuals affected with OAVS phenotypes with autosomal dominant inheritance. Our findings narrow down to three potential candidate genes, SIX1, PDGFRA, and KDR/VEGFR2. Among these, SIX1 has been previously associated with OAVS ear malformations and it is co-expressed with EYA1 during ear development. Attempts to strengthen the genotype-phenotype co-relation underlying the OAVS of phenotypes are essential to discover the etiological factors leading to this complex and burdensome condition as well as for family counseling and prevention efforts.

Ketone flux through BDH1 supports metabolic remodeling of skeletal and cardiac muscles in response to intermittent time-restricted feeding

Time-restricted feeding (TRF) has gained attention as a dietary regimen that promotes metabolic health. This study questioned if the health benefits of an intermittent TRF (iTRF) schedule require ketone flux specifically in skeletal and cardiac muscles. Notably, we found that the ketolytic enzyme beta-hydroxybutyrate dehydrogenase 1 (BDH1) is uniquely enriched in isolated mitochondria derived from heart and red/oxidative skeletal muscles, which also have high capacity for fatty acid oxidation (FAO). Using mice with BDH1 deficiency in striated muscles, we discover that this enzyme optimizes FAO efficiency and exercise tolerance during acute fasting. Additionally, iTRF leads to robust molecular remodeling of muscle tissues, and muscle BDH1 flux does indeed play an essential role in conferring the full adaptive benefits of this regimen, including increased lean mass, mitochondrial hormesis, and metabolic rerouting of pyruvate. In sum, ketone flux enhances mitochondrial bioenergetics and supports iTRF-induced remodeling of skeletal muscle and heart.

3D genomic features across >50 diverse cell types reveal insights into the genomic architecture of childhood obesity

The prevalence of childhood obesity is increasing worldwide, along with the associated common comorbidities of type 2 diabetes and cardiovascular disease in later life. Motivated by evidence for a strong genetic component, our prior genome-wide association study (GWAS) efforts for childhood obesity revealed 19 independent signals for the trait; however, the mechanism of action of these loci remains to be elucidated. To molecularly characterize these childhood obesity loci we sought to determine the underlying causal variants and the corresponding effector genes within diverse cellular contexts. Integrating childhood obesity GWAS summary statistics with our existing 3D genomic datasets for 57 human cell types, consisting of high-resolution promoter-focused Capture-C/Hi-C, ATAC-seq, and RNA-seq, we applied stratified LD score regression and calculated the proportion of genome-wide SNP heritability attributable to cell type-specific features, revealing pancreatic alpha cell enrichment as the most statistically significant. Subsequent chromatin contact-based fine-mapping was carried out for genome-wide significant childhood obesity loci and their linkage disequilibrium proxies to implicate effector genes, yielded the most abundant number of candidate variants and target genes at the BDNF, ADCY3, TMEM18 and FTO loci in skeletal muscle myotubes and the pancreatic beta-cell line, EndoC-BH1. One novel implicated effector gene, ALKAL2 - an inflammation-responsive gene in nerve nociceptors - was observed at the key TMEM18 locus across multiple immune cell types. Interestingly, this observation was also supported through colocalization analysis using expression quantitative trait loci (eQTL) derived from the Genotype-Tissue Expression (GTEx) dataset, supporting an inflammatory and neurologic component to the pathogenesis of childhood obesity. Our comprehensive appraisal of 3D genomic datasets generated in a myriad of different cell types provides genomic insights into pediatric obesity pathogenesis.

Analysis of brain and blood single-cell transcriptomics in acute and subacute phases after experimental stroke

Cerebral ischemia triggers a powerful inflammatory reaction involving peripheral leukocytes and brain resident cells that contribute to both tissue injury and repair. However, their dynamics and diversity remain poorly understood. To address these limitations, we performed a single-cell transcriptomic study of brain and blood cells 2 or 14 days after ischemic stroke in mice. We observed a strong divergence of post-ischemic microglia, monocyte-derived macrophages and neutrophils over time, while endothelial cells and brain-associated macrophages showed altered transcriptomic signatures at 2 days poststroke. Trajectory inference predicted the in situ trans-differentiation of macrophages from blood monocytes into day 2 and day 14 phenotypes, while neutrophils were projected to be continuously de novo recruited from the blood. Brain single-cell transcriptomes from both female and male aged mice were similar to that of young male mice, but aged and young brains differed in their immune cell composition. Although blood leukocyte analysis also revealed altered transcriptomes after stroke, brain-infiltrating leukocytes displayed higher transcriptomic divergence than their circulating counterparts, indicating that phenotypic diversification occurs within the brain in the early and recovery phases of ischemic stroke. A portal ( https://anratherlab.shinyapps.io/strokevis/ ) is provided to allow user-friendly access to our data.

CDKN2A-p16 Deletion and Activated KRASG12D Drive Barrett's-Like Gland Hyperplasia-Metaplasia and Synergize in the Development of Dysplasia Precancer Lesions

BACKGROUND & AIMS

Barrett's esophagus is the precursor of esophageal dysplasia and esophageal adenocarcinoma. CDKN2A-p16 deletions were reported in 34%-74% of patients with Barrett's esophagus who progressed to dysplasia and esophageal adenocarcinoma, suggesting that p16 loss may drive neoplastic progression. KRAS activation frequently occurs in esophageal adenocarcinoma and precancer lesions. LGR5+ stem cells in the squamocolumnar-junction (SCJ) of mouse stomach contribute as Barrett's esophagus progenitors. We aimed to determine the functional effects of p16 loss and KRAS activation in Barrett's-like metaplasia and dysplasia development.

METHODS

We established mouse models with conditional knockout of CDKN2A-p16 (p16KO) and/or activated KRASG12D expression targeting SCJ LGR5+ cells in interleukin 1b transgenic mice and characterized histologic alterations (mucous-gland hyperplasia/metaplasia, inflammation, and dysplasia) in mouse SCJ. Gene expression was determined by microarray, RNA sequencing, and immunohistochemistry of SCJ tissues and cultured 3-dimensional organoids.

RESULTS

p16KO mice exhibited increased mucous-gland hyperplasia/metaplasia versus control mice (P = .0051). Combined p16KO+KRASG12D resulted in more frequent dysplasia and higher dysplasia scores (P = .0036), with 82% of p16KO+KRASG12D mice developing high-grade dysplasia. SCJ transcriptome analysis showed several activated pathways in p16KO versus control mice (apoptosis, tumor necrosis factor-α/nuclear factor-kB, proteasome degradation, p53 signaling, MAPK, KRAS, and G1-to-S transition).

CONCLUSIONS

p16 deletion in LGR5+ cell precursors triggers increased SCJ mucous-gland hyperplasia/metaplasia. KRASG12D synergizes with p16 deletion resulting in higher grades of SCJ glandular dysplasia, mimicking Barrett's high-grade dysplasia. These genetically modified mouse models establish a functional role of p16 and activated KRAS in the progression of Barrett's-like lesions to dysplasia in mice, representing an in vivo model of esophageal adenocarcinoma precancer. Derived 3-dimensional organoid models further provide in vitro modeling opportunities of esophageal precancer stages.

Longitudinal plasma proteomics reveals biomarkers of alveolar-capillary barrier disruption in critically ill COVID-19 patients

The pathobiology of respiratory failure in COVID-19 consists of a complex interplay between viral cytopathic effects and a dysregulated host immune response. In critically ill patients, imatinib treatment demonstrated potential for reducing invasive ventilation duration and mortality. Here, we perform longitudinal profiling of 6385 plasma proteins in 318 hospitalised patients to investigate the biological processes involved in critical COVID-19, and assess the effects of imatinib treatment. Nine proteins measured at hospital admission accurately predict critical illness development. Next to dysregulation of inflammation, critical illness is characterised by pathways involving cellular adhesion, extracellular matrix turnover and tissue remodelling. Imatinib treatment attenuates protein perturbations associated with inflammation and extracellular matrix turnover. These proteomic alterations are contextualised using external pulmonary RNA-sequencing data of deceased COVID-19 patients and imatinib-treated Syrian hamsters. Together, we show that alveolar capillary barrier disruption in critical COVID-19 is reflected in the plasma proteome, and is attenuated with imatinib treatment. This study comprises a secondary analysis of both clinical data and plasma samples derived from a clinical trial that was registered with the EU Clinical Trials Register (EudraCT 2020-001236-10, https://www.clinicaltrialsregister.eu/ctr-search/trial/2020-001236-10/NL ) and Netherlands Trial Register (NL8491, https://www.trialregister.nl/trial/8491 ).

The chromatin landscape of healthy and injured cell types in the human kidney

There is a need to define regions of gene activation or repression that control human kidney cells in states of health, injury, and repair to understand the molecular pathogenesis of kidney disease and design therapeutic strategies. Comprehensive integration of gene expression with epigenetic features that define regulatory elements remains a significant challenge. We measure dual single nucleus RNA expression and chromatin accessibility, DNA methylation, and H3K27ac, H3K4me1, H3K4me3, and H3K27me3 histone modifications to decipher the chromatin landscape and gene regulation of the kidney in reference and adaptive injury states. We establish a spatially-anchored epigenomic atlas to define the kidney's active, silent, and regulatory accessible chromatin regions across the genome. Using this atlas, we note distinct control of adaptive injury in different epithelial cell types. A proximal tubule cell transcription factor network of ELF3, KLF6, and KLF10 regulates the transition between health and injury, while in thick ascending limb cells this transition is regulated by NR2F1. Further, combined perturbation of ELF3, KLF6, and KLF10 distinguishes two adaptive proximal tubular cell subtypes, one of which manifested a repair trajectory after knockout. This atlas will serve as a foundation to facilitate targeted cell-specific therapeutics by reprogramming gene regulatory networks.

Disruption of lysosomal proteolysis in astrocytes facilitates midbrain organoid proteostasis failure in an early-onset Parkinson's disease model

Accumulation of advanced glycation end products (AGEs) on biopolymers accompanies cellular aging and drives poorly understood disease processes. Here, we studied how AGEs contribute to development of early onset Parkinson's Disease (PD) caused by loss-of-function of DJ1, a protein deglycase. In induced pluripotent stem cell (iPSC)-derived midbrain organoid models deficient for DJ1 activity, we find that lysosomal proteolysis is impaired, causing AGEs to accumulate, α-synuclein (α-syn) phosphorylation to increase, and proteins to aggregate. We demonstrated these processes are at least partly driven by astrocytes, as DJ1 loss reduces their capacity to provide metabolic support and triggers acquisition of a pro-inflammatory phenotype. Consistently, in co-cultures, we find that DJ1-expressing astrocytes are able to reverse the proteolysis deficits of DJ1 knockout midbrain neurons. In conclusion, astrocytes' capacity to clear toxic damaged proteins is critical to preserve neuronal function and their dysfunction contributes to the neurodegeneration observed in a DJ1 loss-of-function PD model.

Repression of latent NF-κB enhancers by PDX1 regulates β cell functional heterogeneity

Interactions between lineage-determining and activity-dependent transcription factors determine single-cell identity and function within multicellular tissues through incompletely known mechanisms. By assembling a single-cell atlas of chromatin state within human islets, we identified β cell subtypes governed by either high or low activity of the lineage-determining factor pancreatic duodenal homeobox-1 (PDX1). β cells with reduced PDX1 activity displayed increased chromatin accessibility at latent nuclear factor κB (NF-κB) enhancers. Pdx1 hypomorphic mice exhibited de-repression of NF-κB and impaired glucose tolerance at night. Three-dimensional analyses in tandem with chromatin immunoprecipitation (ChIP) sequencing revealed that PDX1 silences NF-κB at circadian and inflammatory enhancers through long-range chromatin contacts involving SIN3A. Conversely, Bmal1 ablation in β cells disrupted genome-wide PDX1 and NF-κB DNA binding. Finally, antagonizing the interleukin (IL)-1β receptor, an NF-κB target, improved insulin secretion in Pdx1 hypomorphic islets. Our studies reveal functional subtypes of single β cells defined by a gradient in PDX1 activity and identify NF-κB as a target for insulinotropic therapy.

Potential Biomarkers for Noninfectious Scleritis Identified by Serum and Tear Fluid Proteomics

Purpose

Scleritis is an extremely painful and potentially blinding inflammation of the sclera with unknown pathogenesis and unpredictable course. To gain insight in its disease process and identify biomarker candidates, we performed extensive proteomics in serum and tear fluid.

Design

Prospective multicenter cohort study.

Participants

A total of 121 patients with noninfectious scleritis (of which 39 active cases), 30 healthy controls, and 23 disease controls (uveitis and rheumatoid arthritis) were enrolled in the Netherlands from 2020 to 2022.

Methods

Serum, tear fluid of both eyes, and clinical data were gathered. The level of 368 inflammatory proteins was measured using proximity extension assays. Results were validated in an independent cohort of 15 patients with scleritis, and using addressable laser bead immunoassay, or enzyme-linked immunoassays. In addition, we studied an extended panel of matrix metalloproteinases in tear fluid of necrotizing scleritis with addressable laser bead immunoassay.

Main Outcome Measures

Statistically significant differences in the level of inflammatory proteins between patients with scleritis and control groups.

Results

Proteomics revealed 18 significantly upregulated or downregulated serum proteins in active scleritis cases compared with all control groups in both the discovery cohort and the validation cohort. The most upregulated protein was nuclear migration protein nudC (NudC; P = 0.0032), a protein involved in neurogenesis. The other significant hits included proteins involved in T-cell activation, apoptosis, epithelial barrier maintenance, and angiogenesis. Our tear fluid analysis showed matrix metalloproteinase 9 (MMP9) to be upregulated in the tear fluid of patients with scleral necrosis.

Conclusions

The results of our proteomics analysis suggest a role for neurogenesis, T-cell activation, disruption of epithelial barrier, and angiogenesis in the pathogenesis of scleritis, and highlight MMP9 and NudC as biomarkers with potential clinical relevance.

Funding Disclosures

The authors have no proprietary or commercial interest in any materials discussed in this article.