Cell-type Specificity Research Articles

BiSpec Pairwise AI: guiding the selection of bispecific antibody target combinations with pairwise learning and GPT augmentation

PurposeBispecific antibodies (BsAbs), capable of targeting two antigens simultaneously, represent a significant advancement by employing dual mechanisms of action for tumor suppression. However, how to pair targets to develop effective and safe bispecific drugs is a major challenge for pharmaceutical companies.MethodsUsing machine learning models, we refined the biological characteristics of currently approved or in clinical development BsAbs and analyzed hundreds of membrane proteins as bispecific targets to predict the likelihood of successful drug development for various target combinations. Moreover, to enhance the interpretability of prediction results in bispecific target combination, we combined machine learning models with Large Language Models (LLMs). Through a Retrieval-Augmented Generation (RAG) approach, we supplement each pair of bispecific targets’ machine learning prediction with important features and rationales, generating interpretable analytical reports.ResultsIn this study, the XGBoost model with pairwise learning was employed to predict the druggability of BsAbs. By analyzing extensive data on BsAbs and designing features from perspectives such as target activity, safety, cell type specificity, pathway mechanism, and gene embedding representation, our model is able to predict target combinations of BsAbs with high market potential. Specifically, we integrated XGBoost with the GPT model to discuss the efficacy of each bispecific target pair, thereby aiding the decision-making for drug developers.ConclusionThe novelty of this study lies in the integration of machine learning and GPT techniques to provide a novel framework for the design of BsAbs drugs. This holistic approach not only improves prediction accuracy, but also enhances the interpretability and innovativeness of drug design.

Investigating the role of local acidosis in hypothalamic vasopressin cell excitation during hypoxia

The supraoptic nucleus (SON) of the hypothalamus plays a key role in the regulation of fluid homeostasis through the synthesis and release of vasopressin (VP). In this population, VP release can occur at both dendritic and axonal sites, allowing for both localized and distal effects such as renal water retention and vasoconstriction. Recent work from our laboratory demonstrated an inverse neurovascular coupling (iNVC) response in which VP cell activation by acute salt loading resulted in dendritic release of VP, leading to localized vasoconstriction and hypoxia. Additionally, we demonstrated that the hypoxic milieu in turn potentiates VP cell activity and increases firing, suggesting that these mechanism constitutes a positive feedback loop that optimized VP activation in order to cope with the homeostatic challenge (Roy et al., Cell Reports 2021). Still, the biophysical mechanism by which the iNVC-hypoxia elicits increased firing in VP neurons has not yet been explored. We hypothesized that hypoxia-induced acidosis, via activation of acid sensitive ion channels (ASICs) (which are present in SON VP neurons; Ohbuchi et al in 2010) is a key mechanism contributing to VP cell excitation in the iNVC positive feedback loop. We investigated this with whole cell patch clamp electrophysiology using eGFP-VP transgenic rats. We recorded from identified SON VP cells and measured deflections in voltage as bath pH was changed from 7.25 to 6.3 (10 mins). Bath solution was either HCL-buffered artificial cerebrospinal fluid (aCSF) aerated with 95% O2/5% CO2 or HEPES-buffered solution (HBS). We observed an excitatory response to the decreases in pH in most SON VP cells recorded (68.75%, n = 11/16). These responses were reversible and, in many cases, repeatable following a second stimulation. To examine cell-type specificity, ongoing studies are focusing on comparing these responses to those obtained from oxytocin cells. Further, we plan to repeat these experiments in the presence of the ASIC-blocker amiloride and determine whether the hypoxia-induced excitatory effect during iNVC is indeed dependent on ASIC-channels. Together, our results support our hypothesis that tissue acidosis (likely via activation of ASIC-channels) contributes to the excitatory feedback loop evoked during salt-induced iNVC in VP neurons. This abstract was supported by the National Heart Lung and Blood Institute of the National Institutes of Health [R01 HL162575-01 (J.E.S.)]. This is the full abstract presented at the American Physiology Summit 2024 meeting and is only available in HTML format. There are no additional versions or additional content available for this abstract. Physiology was not involved in the peer review process.

Prediction of cell-type-specific cohesin-mediated chromatin loops based on chromatin state

Chromatin loop is of crucial importance for the regulation of gene transcription. Cohesin is a type of chromatin-associated protein that mediates the interaction of chromatin through the loop extrusion. Cohesin-mediated chromatin interactions have strong cell-type specificity, posing a challenge for predicting chromatin loops. Existing computational methods perform poorly in predicting cell-type-specific chromatin loops. To address this issue, we propose a random forest model to predict cell-type-specific cohesin-mediated chromatin loops based on chromatin states identified by ChromHMM and the occupancy of related factors. Our results show that chromatin state is responsible for cell-type-specificity of loops. Using only chromatin states as features, the model achieved high accuracy in predicting cell-type-specific loops between two cell types and can be applied to different cell types. Furthermore, when chromatin states are combined with the occurrence frequency of CTCF, RAD21, YY1, and H3K27ac ChIP-seq peaks, more accurate prediction can be achieved. Our feature extraction method provides novel insights into predicting cell-type-specific chromatin loops and reveals the relationship between chromatin state and chromatin loop formation.

Expression characteristics and potential function of non-coding RNA in mouse cortical cells.

Non-coding RNAs (ncRNAs) play essential regulatory functions in various physiological and pathological processes in the brain. To systematically characterize the ncRNA profile in cortical cells, we downloaded single-cell SMART-Seq v4 data of mouse cerebral cortex. Our results revealed that the ncRNAs alone are sufficient to define the identity of most cortical cell types. We identified 1,600 ncRNAs that exhibited cell type specificity, even yielding to distinguish microglia from perivascular macrophages with ncRNA. Moreover, we characterized cortical layer and region specific ncRNAs, in line with the results by spatial transcriptome (ST) data. By constructing a co-expression network of ncRNAs and protein-coding genes, we predicted the function of ncRNAs. By integrating with genome-wide association studies data, we established associations between cell type-specific ncRNAs and traits related to neurological disorders. Collectively, our study identified differentially expressed ncRNAs at multiple levels and provided the valuable resource to explore the functions and dysfunctions of ncRNAs in cortical cells.

Single-cell long-read sequencing-based mapping reveals specialized splicing patterns in developing and adult mouse and human brain

RNA isoforms influence cell identity and function. However, a comprehensive brain isoform map was lacking. We analyze single-cell RNA isoforms across brain regions, cell subtypes, developmental time points and species. For 72% of genes, full-length isoform expression varies along one or more axes. Splicing, transcription start and polyadenylation sites vary strongly between cell types, influence protein architecture and associate with disease-linked variation. Additionally, neurotransmitter transport and synapse turnover genes harbor cell-type variability across anatomical regions. Regulation of cell-type-specific splicing is pronounced in the postnatal day 21-to-postnatal day 28 adolescent transition. Developmental isoform regulation is stronger than regional regulation for the same cell type. Cell-type-specific isoform regulation in mice is mostly maintained in the human hippocampus, allowing extrapolation to the human brain. Conversely, the human brain harbors additional cell-type specificity, suggesting gain-of-function isoforms. Together, this detailed single-cell atlas of full-length isoform regulation across development, anatomical regions and species reveals an unappreciated degree of isoform variability across multiple axes.

Global spatiotemporal synchronizing structures of spontaneous neural activities in different cell types

Increasing evidence has revealed the large-scale nonstationary synchronizations as traveling waves in spontaneous neural activity. However, the interplay of various cell types in fine-tuning these spatiotemporal patters remains unclear. Here, we performed comprehensive exploration of spatiotemporal synchronizing structures across different cell types, states (awake, anesthesia, motion) and developmental axis in male mice. We found traveling waves in glutamatergic neurons exhibited greater variety than those in GABAergic neurons. Moreover, the synchronizing structures of GABAergic neurons converged toward those of glutamatergic neurons during development, but the evolution of waves exhibited varying timelines for different sub-type interneurons. Functional connectivity arises from both standing and traveling waves, and negative connections can be elucidated by the spatial propagation of waves. In addition, some traveling waves were correlated with the spatial distribution of gene expression. Our findings offer further insights into the neural underpinnings of traveling waves, functional connectivity, and resting-state networks, with cell-type specificity and developmental perspectives.

Morphological correlates of pyramidal cell axonal myelination in mouse and human neocortex.

The axons of neocortical pyramidal neurons are frequently myelinated. Heterogeneity in the topography of axonal myelination in the cerebral cortex has been attributed to a combination of electrophysiological activity, axonal morphology, and neuronal-glial interactions. Previously, we showed that axonal segment length and caliber are critical local determinants of fast-spiking interneuron myelination. However, the factors that determine the myelination of individual axonal segments along neocortical pyramidal neurons remain largely unexplored. Here, we used structured illumination microscopy to examine the extent to which axonal morphology is predictive of the topography of myelination along neocortical pyramidal neurons. We identified critical thresholds for axonal caliber and interbranch distance that are necessary, but not sufficient, for myelination of pyramidal cell axons in mouse primary somatosensory cortex (S1). Specifically, we found that pyramidal neuron axonal segments with a caliber < 0.24μm or interbranch distance < 18.10μm are rarely myelinated. Moreover, we further confirmed that these findings in mice are similar for human neocortical pyramidal cell myelination (caliber < 0.25μm, interbranch distance < 19.00μm), suggesting that this mechanism is evolutionarily conserved. Taken together, our findings suggest that axonal morphology is a critical correlate of the topography and cell-type specificity of neocortical myelination.

Integrative multi-omics analysis identifies genetically supported druggable targets and immune cell specificity for myasthenia gravis

BackgroundMyasthenia gravis (MG) is a chronic autoimmune disorder characterized by fluctuating muscle weakness. Despite the availability of established therapies, the management of MG symptoms remains suboptimal, partially attributed to lack of efficacy or intolerable side-effects. Therefore, new effective drugs are warranted for treatment of MG.MethodsBy employing an analytical framework that combines Mendelian randomization (MR) and colocalization analysis, we estimate the causal effects of blood druggable expression quantitative trait loci (eQTLs) and protein quantitative trait loci (pQTLs) on the susceptibility of MG. We subsequently investigated whether potential genetic effects exhibit cell-type specificity by utilizing genetic colocalization analysis to assess the interplay between immune-cell-specific eQTLs and MG risk.ResultsWe identified significant MR results for four genes (CDC42BPB, CD226, PRSS36, and TNFSF12) using cis-eQTL genetic instruments and three proteins (CTSH, PRSS8, and CPN2) using cis-pQTL genetic instruments. Six of these loci demonstrated evidence of colocalization with MG susceptibility (posterior probability > 0.80). We next undertook genetic colocalization to investigate cell-type-specific effects at these loci. Notably, we identified robust evidence of colocalization, with a posterior probability of 0.854, linking CTSH expression in TH2 cells and MG risk.ConclusionsThis study provides crucial insights into the genetic and molecular factors associated with MG susceptibility, singling out CTSH as a potential candidate for in-depth investigation and clinical consideration. It additionally sheds light on the immune-cell regulatory mechanisms related to the disease. However, further research is imperative to validate these targets and evaluate their feasibility for drug development.

Abstract 3024: HIF1A and HIF2A differentially contribute to early and adaptive changes following in vivo Vhl inactivation in the kidney

Abstract The von Hippel Lindau (VHL) tumor suppressor gene is a classical tumor suppressor whose biallelic inactivation is a truncal event in clear cell renal cell carcinoma (ccRCC) evolution. VHL is best characterized as the recognition component of an E3 ubiquitin ligase complex that targets hypoxia inducible transcription factors (HIF1A and HIF2A) for conditional oxygen-dependent degradation. HIF activation following VHL inactivation is an early event in ccRCC progression. While the functions of HIF1A and HIF2A in ccRCC biology and progression have been studied, those in VHL-null non-neoplastic cells of the renal tubular epithelium (RTE) and their contribution to tumorigenesis have not been fully explored. Here, using a combination of lineage marking and single-cell RNA sequencing (scRNA-seq), we provide the first description of HIF1A- and/or HIF2A-dependent and -independent consequences of in vivo Vhl inactivation in the RTE with single-cell resolution. We have previously reported on the development of a novel mouse model designed to mark and track cells that have undergone Vhl inactivation in vivo. Using this model, we have shown that Vhl inactivation has extensive, cell-type specific effects on gene expression in cells of the RTE within 3 weeks of recombination. Furthermore, we have shown that over a course of 4-8 months after Vhl inactivation, Vhl-null cells undergo specific adaptive changes in gene expression that encompass genes not regulated early after Vhl inactivation. In the work reported here, we have coupled this mouse model to a tamoxifen-inducible RTE-specific Pax8-CreERT2 and conditional alleles for Hif1a and/or Epas1 (codes for HIF2A). We have analyzed ‘marked’ recombined cells at early (1-3 weeks) or late (4-8 months) timepoints after recombination and compared these cells to time-matched Vhl-null cells using immunohistochemistry, in situ RNA hybridization and scRNA-seq. We found that though both HIF isoforms were present in Vhl-null cells at the early timepoint, expression of genes regulated early after Vhl inactivation was affected more by Hif1a inactivation than Epas1 inactivation. In contrast, expression of genes that are regulated in Vhl-null cells over time was drastically impacted by Epas1 deletion but was unaffected by Hif1a deletion. Interestingly, most of these Epas1-dependent, time-dependent changes in Vhl-null cells involved genes that have not been identified as direct transcriptional targets of HIF. Collectively, our work suggests that HIF1A and HIF2A have temporally distinct contributions to gene expression in Vhl-null cells in the RTE. We are currently analyzing the cell-type specificity of these contributions and ascertaining their relevance to the adoption of ccRCC-associated gene expression profiles. Citation Format: Samvid Kurlekar, Joanna Lima, Ayslan B. Barros, Norma Masson, Sophia Zhai, Christopher W. Pugh, Julie Adam, Peter J. Ratcliffe. HIF1A and HIF2A differentially contribute to early and adaptive changes following in vivo Vhl inactivation in the kidney [abstract]. In: Proceedings of the American Association for Cancer Research Annual Meeting 2024; Part 1 (Regular Abstracts); 2024 Apr 5-10; San Diego, CA. Philadelphia (PA): AACR; Cancer Res 2024;84(6_Suppl):Abstract nr 3024.

Abstract 7355: Integration analysis of lncRNAs and R-loops in 3D genome organization and gene regulation in the cancer genome

Abstract Existing evidence shows that some long non-coding RNAs (lncRNAs) can recruit chromatin regulators and form R-loops to mediate CCTCC-binding factor (CTCF) binding at topologically associated domain (TAD) boundaries to activate transcription of target genes. However, it remains unclear to what extent such lncRNA-dependent chromatin organization events occur genome-wide in normal and cancer cells. By curated public data of lncRNA, CTCF, RAD21 profiles and integrative analysis of multi-omics data using innovative computational approaches, we investigated global associations of lncRNAs and R-loops with activities of transcription factors and chromatin regulators in gene expression regulations across multiple human cell types. We developed a comprehensive computational pipeline for cross-cell-type analysis of lncRNA and protein-coding gene expression profiles and transcription factor binding profiles. We identified a total of 542 candidate genomic regions where lncRNAs play a pivotal role in regulating target genes by facilitating CTCF/cohesin binding. R-loops are enriched at CTCF/cohesin binding sites, suggesting their potential contribution to CTCF binding. Additionally, our findings underscore the significant cell-type specificity of CTCF binding patterns across the genome. This work adds knowledge to cancer genome organization and gene regulation, and builds foundations for further functional and mechanistic studies of lncRNAs and R-loops in the human genome. Citation Format: Shengyuan Wang, Zhenjia Wang, Yuh-Hwa Wang, Hui Li, Peter Todd Stukenburg, Chongzhi Zang. Integration analysis of lncRNAs and R-loops in 3D genome organization and gene regulation in the cancer genome [abstract]. In: Proceedings of the American Association for Cancer Research Annual Meeting 2024; Part 1 (Regular Abstracts); 2024 Apr 5-10; San Diego, CA. Philadelphia (PA): AACR; Cancer Res 2024;84(6_Suppl):Abstract nr 7355.

An Automated and Fast Sample Preparation Workflow for Laser Microdissection Guided Ultrasensitive Proteomics

Spatial tissue proteomics integrating whole-slide imaging, laser microdissection, and ultrasensitive mass spectrometry is a powerful approach to link cellular phenotypes to functional proteome states in (patho)physiology. To be applicable to large patient cohorts and low sample input amounts, including single-cell applications, loss-minimized and streamlined end-to-end workflows are key. We here introduce an automated sample preparation protocol for laser microdissected samples utilizing the cellenONE robotic system, which has the capacity to process 192 samples in 3 h. Following laser microdissection collection directly into the proteoCHIP LF 48 or EVO 96 chip, our optimized protocol facilitates lysis, formalin de-crosslinking, and tryptic digest of low-input archival tissue samples. The seamless integration with the Evosep ONE LC system by centrifugation allows ‘on-the-fly’ sample clean-up, particularly pertinent for laser microdissection workflows. We validate our method in human tonsil archival tissue, where we profile proteomes of spatially-defined B-cell, T-cell, and epithelial microregions of 4000 μm2 to a depth of ∼2000 proteins and with high cell type specificity. We finally provide detailed equipment templates and experimental guidelines for broad accessibility.

Oncogenic Cell Tagging and Single-Cell Transcriptomics Reveal Cell Type-Specific and Time-Resolved Responses to Vhl Inactivation in the Kidney.

Single-cell analysis of heterogeneous and dynamic responses to Vhl inactivation in the kidney suggests that early events shape the cell type specificity of oncogenesis, providing a focus for mechanistic understanding and therapeutic targeting.

Distinct cellular microenvironment with cytotypic effects regulates orderly regeneration of vascular tissues

Regeneration of the architecturally complex blood vascular system requires precise temporal and spatial control of cell behaviours. Additional components must be integrated into the structure to achieve clinical success for in situ tissue engineering. Consequently, this study proposed a universal method for including any substrate type in vascular cell extracellular matrices (VCEM) via regulating selective adhesion to promote vascular tissue regeneration. The results uncovered that the VCEM worked as cell adhesion substrates, exhibited cell type specificity, and functioned as an address signal for recognition by vascular cells, which resulted in matching with the determined cells. The qPCR and immunofluorescence results revealed that a cell type-specific VCEM could be designed to promote or inhibit cell adhesion, consistenting with the expression patterns of eyes absent 3 (Eya3). In addition, a 3D vascular graft combined with VCEM which could recapitulate the vascular cell-like microenvironment was fabricated. The vascular graft revealed a prospective role for cellular microenvironment in the establishment of vascular cell distribution and tissue architecture, and potentiated the orderly regeneration and functional recovery of vascular tissues in vivo. The findings demonstrate that differential adhesion between cell types due to the cellular microenvironment is sufficient to drive the complex assembly of engineered blood vessel functional units, and underlies hierarchical organization during vascular regeneration.

Performance Assessment of Sysmex DI-60: Is Digital Morphology Analyzer Reliable for White Blood Cell Differentials in Body Fluids?

Few studies have evaluated digital morphology (DM) analyzers on body fluids (BF). We evaluated the performance of a DM analyzer, Sysmex DI-60 (Sysmex, Kobe, Japan) for white blood cell (WBC) differentials in BF samples. In five BF samples (two pleural fluids and three ascites) containing a single, dominant cell type (>80%, neutrophils, lymphocytes, macrophages, abnormal lymphocytes, and malignant cells in each sample), we evaluated the precision of the DI-60 and compared the WBC differentials and turnaround times (TAT) between DI-60 and manual counting. The precision of the DI-60 pre-classification and verification was excellent (%CV, 0.01-3.16%). After verification, the DI-60 showed high sensitivity, specificity, and efficiency (ranges: 90.8-98.1%, 96.8-97.9%, and 92.5-98.0%, respectively) for the dominant cell types in neutrophil- and lymphocyte-dominant samples. For all samples, the DI-60 and manual counting showed high correlations for major cell types (neutrophils, lymphocytes, macrophages, and others, r = 0.72 to 0.94) after verification. The agreement between the pre-classification and verification of the DI-60 was strong in the neutrophil-dominant sample (κ = 0.81). The DI-60 showed a significantly longer TAT (min: s) than manual counting for all samples (median TAT/slide: 6:28 vs. 1:53, p < 0.0001), with remarkable differences in abnormal lymphocyte- and malignant cell-dominant samples (21:05 vs. 2:06; 12:34 vs. 2:25). The DI-60 may provide reliable data in neutrophil- and lymphocyte-dominant BF samples. However, it may require longer times and higher workloads for WBC differentials especially in BF samples containing atypical cells. Further improvement would be needed before applying DM analyzers for routine clinical practice in BF analysis.

Regulatory activity is the default DNA state in eukaryotes.

Genomes encode for genes and non-coding DNA, both capable of transcriptional activity. However, unlike canonical genes, many transcripts from non-coding DNA have limited evidence of conservation or function. Here, to determine how much biological noise is expected from non-genic sequences, we quantify the regulatory activity of evolutionarily naive DNA using RNA-seq in yeast and computational predictions in humans. In yeast, more than 99% of naive DNA bases were transcribed. Unlike the evolved transcriptome, naive transcripts frequently overlapped with opposite sense transcripts, suggesting selection favored coherent gene structures in the yeast genome. In humans, regulation-associated chromatin activity is predicted to be common in naive dinucleotide-content-matched randomized DNA. Here, naive and evolved DNA have similar co-occurrence and cell-type specificity of chromatin marks, challenging these as indicators of selection. However, in both yeast and humans, extreme high activities were rare in naive DNA, suggesting they result from selection. Overall, basal regulatory activity seems to be the default, which selection can hone to evolve a function or, if detrimental, repress.

Anterograde trans-neuronal labeling of striatal interneurons in relation to dopamine neurons in the substantia nigra pars compacta.

Recent advances in neural tracing have unveiled numerous neural circuits characterized by brain region and cell type specificity, illuminating the underpinnings of specific functions and behaviors. Dopaminergic (DA) neurons in the midbrain are highly heterogeneous in terms of gene and protein expression and axonal projections. Different cell types within the substantia nigra pars compacta (SNc) tend to project to the striatum in a cell-type-dependent manner characterized by specific topography. Given the wide and dense distribution of DA axons, coupled with a combination of synaptic and volume transmission, it remains unclear how DA release is spatially and temporally regulated, to appropriately achieve specific behaviors and functions. Our hypothesis posits that hidden rules governing synapse formation between pre-synaptic DA neuron types and striatal neuron types may modulate the effect of DA at a single-cell level. To address this conjecture, we employed adeno-associated virus serotype 1 (AAV1) to visualize the neural circuitry of DA neurons. AAV1 has emerged as a potent anatomical instrument capable of labeling and visualizing pre- and post-synaptic neurons simultaneously through anterograde trans-synaptic labeling. First, AAV1-Cre was injected into the SNc, resulting in Cre expression in both medium spiny neurons and interneurons in the striatum. Due to the potential occurrence of the retrograde transfer of AAV1, only striatal interneurons were considered for trans-synaptic or trans-neuronal labeling. Interneuron types expressing parvalbumin, choline acetyltransferase, somatostatin, or nitrogen oxide synthase exhibited Cre expression. Using a combination of AAV1-Cre and Cre-driven fluorophore expressing AAVs, striatal interneurons and the axons originating from the SNc were visualized in distinct colors. Using immunofluorescence against neurotransmitter transporters, almost all axons in the striatum visualized using this approach were confirmed to be dopaminergic. Moreover, individual DA axons established multiple appositions on the somata and proximal dendrites of interneurons. This finding suggests that irrespective of the extensive and widespread axonal arborization of DA neurons, a particular DA neuron may exert a significant influence on specific interneurons. Thus, AAV1-based labeling of the DA system can be a valuable tool to uncover the concealed rules governing these intricate relationships.

Tetherless Optical Neuromodulation: Wavelength from Orange-red to Mid-infrared.

Optogenetics, a technique that employs light for neuromodulation, has revolutionized the study of neural mechanisms and the treatment of neurological disorders due to its high spatiotemporal resolution and cell-type specificity. However, visible light, particularly blue and green light, commonly used in conventional optogenetics, has limited penetration in biological tissue. This limitation necessitates the implantation of optical fibers for light delivery, especially in deep brain regions, leading to tissue damage and experimental constraints. To overcome these challenges, the use of orange-red and infrared light with greater tissue penetration has emerged as a promising approach for tetherless optical neuromodulation. In this review, we provide an overview of the development and applications of tetherless optical neuromodulation methods with long wavelengths. We first discuss the exploration of orange-red wavelength-responsive rhodopsins and their performance in tetherless optical neuromodulation. Then, we summarize two novel tetherless neuromodulation methods using near-infrared light: upconversion nanoparticle-mediated optogenetics and photothermal neuromodulation. In addition, we discuss recent advances in mid-infrared optical neuromodulation.

Identifying novel proteins for suicide attempt by integrating proteomes from brain and blood with genome-wide association data

Genome-wide association studies (GWASs) have identified risk loci for suicide attempt (SA), but deciphering how they confer risk for SA remains largely unknown. This study aims to identify the key proteins and gain insights into SA pathogenesis. We integrated data from the brain proteome (N = 376) and blood proteome (N = 35,559) and combined it with the largest SA GWAS summary statistics to date (N = 518,612). A comprehensive set of methods was employed, including Mendelian randomization (MR), Steiger filtering, Bayesian colocalization, proteome‑wide association studies (PWAS), transcript-levels, cell-type specificity, correlation, and protein-protein interaction (PPI) network analysis. Validation was performed using other protein datasets and the SA dataset from FinnGen study. We identified ten proteins (GLRX5, GMPPB, B3GALTL, FUCA2, TTLL12, ADCK1, MMAA, HIBADH, ACP1, DOC2A) associated with SA in brain proteomics. GLRX5, GMPPB, and FUCA2 showed strong colocalization evidence and were supported by PWAS and transcript-level analysis, and were predominantly expressed in glutamatergic neuronal cells. In blood proteomics, one significant protein (PEAR1) and three near-significant proteins (NDE1, EVA1C, B4GALT2) were identified, but lacked colocalization evidence. Moreover, despite the limited correlation between the same protein in brain and blood, the PPI network analysis provided new insights into the interaction between brain and blood in SA. Furthermore, GLRX5 was associated with the GSTP1, the target of Clozapine. The comprehensive analysis provides strong evidence supporting a causal association between three genetically determined brain proteins (GLRX5, GMPPB, and FUCA2) with SA. These findings offer valuable insights into SA’s underlying mechanisms and potential therapeutic approaches.

Integrative genomic analyses reveal putative cell type-specific targets of the Drosophila ets transcription factor Pointed

The Ets domain transcription factors direct diverse biological processes throughout all metazoans and are implicated in development as well as in tumor initiation, progression and metastasis. The Drosophila Ets transcription factor Pointed (Pnt) is the downstream effector of the Epidermal growth factor receptor (Egfr) pathway and is required for cell cycle progression, specification, and differentiation of most cell types in the larval eye disc. Despite its critical role in development, very few targets of Pnt have been reported previously. Here, we employed an integrated approach by combining genome-wide single cell and bulk data to identify putative cell type-specific Pnt targets. First, we used chromatin immunoprecipitation with high-throughput sequencing (ChIP-seq) to determine the genome-wide occupancy of Pnt in late larval eye discs. We identified enriched regions that mapped to an average of 6,941 genes, the vast majority of which are novel putative Pnt targets. Next, we integrated ChIP-seq data with two other larval eye single cell genomics datasets (scRNA-seq and snATAC-seq) to reveal 157 putative cell type-specific Pnt targets that may help mediate unique cell type responses upon Egfr-induced differentiation. Finally, our integrated data also predicts cell type-specific functional enhancers that were not reported previously. Together, our study provides a greatly expanded list of putative cell type-specific Pnt targets in the eye and is a resource for future studies that will allow mechanistic insights into complex developmental processes regulated by Egfr signaling.

PancrESS – a meta-analysis resource for understanding cell-type specific expression in the human pancreas

BackgroundThe human pancreas is composed of specialized cell types producing hormones and enzymes critical to human health. These specialized functions are the result of cell type-specific transcriptional programs which manifest in cell-specific gene expression. Understanding these programs is essential to developing therapies for pancreatic disorders. Transcription in the human pancreas has been widely studied by single-cell RNA technologies, however the diversity of protocols and analysis methods hinders their interpretability in the aggregate.ResultsIn this work, we perform a meta-analysis of pancreatic single-cell RNA sequencing data. We present a database for reference transcriptome abundances and cell-type specificity metrics. This database facilitates the identification and definition of marker genes within the pancreas. Additionally, we introduce a versatile tool which is freely available as an R package, and should permit integration into existing workflows. Our tool accepts count data files generated by widely-used single-cell gene expression platforms in their original format, eliminating an additional pre-formatting step. Although we designed it to calculate expression specificity of pancreas cell types, our tool is agnostic to the biological source of count data, extending its applicability to other biological systems.ConclusionsOur findings enhance the current understanding of expression specificity within the pancreas, surpassing previous work in terms of scope and detail. Furthermore, our database and tool enable researchers to perform similar calculations in diverse biological systems, expanding the applicability of marker gene identification and facilitating comparative analyses.