Cell-type Specificity Research Articles

ScGWAS: landscape of trait-cell type associations by integrating single-cell transcriptomics-wide and genome-wide association studies

BackgroundThe rapid accumulation of single-cell RNA sequencing (scRNA-seq) data presents unique opportunities to decode the genetically mediated cell-type specificity in complex diseases. Here, we develop a new method, scGWAS, which effectively leverages scRNA-seq data to achieve two goals: (1) to infer the cell types in which the disease-associated genes manifest and (2) to construct cellular modules which imply disease-specific activation of different processes.ResultsscGWAS only utilizes the average gene expression for each cell type followed by virtual search processes to construct the null distributions of module scores, making it scalable to large scRNA-seq datasets. We demonstrated scGWAS in 40 genome-wide association studies (GWAS) datasets (average sample size N ≈ 154,000) using 18 scRNA-seq datasets from nine major human/mouse tissues (totaling 1.08 million cells) and identified 2533 trait and cell-type associations, each with significant modules for further investigation. The module genes were validated using disease or clinically annotated references from ClinVar, OMIM, and pLI variants.ConclusionsWe showed that the trait-cell type associations identified by scGWAS, while generally constrained to trait-tissue associations, could recapitulate many well-studied relationships and also reveal novel relationships, providing insights into the unsolved trait-tissue associations. Moreover, in each specific cell type, the associations with different traits were often mediated by different sets of risk genes, implying disease-specific activation of driving processes. In summary, scGWAS is a powerful tool for exploring the genetic basis of complex diseases at the cell type level using single-cell expression data.

Dynamic proteomic and phosphoproteomic atlas of corticostriatal axons in neurodevelopment.

Mammalian axonal development begins in embryonic stages and continues postnatally. After birth, axonal proteomic landscape changes rapidly, coordinated by transcription, protein turnover, and post-translational modifications. Comprehensive profiling of axonal proteomes across neurodevelopment is limited, with most studies lacking cell-type and neural circuit specificity, resulting in substantial information loss. We create a Cre-dependent APEX2 reporter mouse line and map cell-type-specific proteome of corticostriatal projections across postnatal development. We synthesize analysis frameworks to define temporal patterns of axonal proteome and phosphoproteome, identifying co-regulated proteins and phosphorylations associated with genetic risk for human brain disorders. We discover proline-directed kinases as major developmental regulators. APEX2 transgenic reporter proximity labeling offers flexible strategies for subcellular proteomics with cell type specificity in early neurodevelopment, a critical period for neuropsychiatric disease.

Epigenetics of neural differentiation: Spotlight on enhancers.

Neural induction, both in vivo and in vitro, includes cellular and molecular changes that result in phenotypic specialization related to specific transcriptional patterns. These changes are achieved through the implementation of complex gene regulatory networks. Furthermore, these regulatory networks are influenced by epigenetic mechanisms that drive cell heterogeneity and cell-type specificity, in a controlled and complex manner. Epigenetic marks, such as DNA methylation and histone residue modifications, are highly dynamic and stage-specific during neurogenesis. Genome-wide assessment of these modifications has allowed the identification of distinct non-coding regulatory regions involved in neural cell differentiation, maturation, and plasticity. Enhancers are short DNA regulatory regions that bind transcription factors (TFs) and interact with gene promoters to increase transcriptional activity. They are of special interest in neuroscience because they are enriched in neurons and underlie the cell-type-specificity and dynamic gene expression profiles. Classification of the full epigenomic landscape of neural subtypes is important to better understand gene regulation in brain health and during diseases. Advances in novel next-generation high-throughput sequencing technologies, genome editing, Genome-wide association studies (GWAS), stem cell differentiation, and brain organoids are allowing researchers to study brain development and neurodegenerative diseases with an unprecedented resolution. Herein, we describe important epigenetic mechanisms related to neurogenesis in mammals. We focus on the potential roles of neural enhancers in neurogenesis, cell-fate commitment, and neuronal plasticity. We review recent findings on epigenetic regulatory mechanisms involved in neurogenesis and discuss how sequence variations within enhancers may be associated with genetic risk for neurological and psychiatric disorders.

Postnatal developmental trajectory of sex-biased gene expression in the mouse pituitary gland

BackgroundThe pituitary gland regulates essential physiological processes such as growth, pubertal onset, stress response, metabolism, reproduction, and lactation. While sex biases in these functions and hormone production have been described, the underlying identity, temporal deployment, and cell-type specificity of sex-biased pituitary gene regulatory networks are not fully understood.MethodsTo capture sex differences in pituitary gene regulation dynamics during postnatal development, we performed 3’ untranslated region sequencing and small RNA sequencing to ascertain gene and microRNA expression, respectively, across five postnatal ages (postnatal days 12, 22, 27, 32, 37) that span the pubertal transition in female and male C57BL/6J mouse pituitaries (n = 5–6 biological replicates for each sex at each age).ResultsWe observed over 900 instances of sex-biased gene expression and 17 sex-biased microRNAs, with the majority of sex differences occurring with puberty. Using miRNA–gene target interaction databases, we identified 18 sex-biased genes that were putative targets of 5 sex-biased microRNAs. In addition, by combining our bulk RNA-seq with publicly available male and female mouse pituitary single-nuclei RNA-seq data, we obtained evidence that cell-type proportion sex differences exist prior to puberty and persist post-puberty for three major hormone-producing cell types: somatotropes, lactotropes, and gonadotropes. Finally, we identified sex-biased genes in these three pituitary cell types after accounting for cell-type proportion differences between sexes.ConclusionOur study reveals the identity and postnatal developmental trajectory of sex-biased gene expression in the mouse pituitary. This work also highlights the importance of considering sex biases in cell-type composition when understanding sex differences in the processes regulated by the pituitary gland.

Characterization of two novel neutral sphingomyelinase 2 inhibitors in endosomal sorting and Extracellular Vesicle biogenesis

Sphingomyelinase hydrolyzes the phosphodiester bond of the sphingomyelin to ceramide and phosphorylcholine and have been involved in extracellular vesicle (EV) biogenesis and more recently in membrane repair. Here we describe an initial testing of two recently discovered neutral sphingomyelinase 2 (nSMase2) inhibitors ((R)-(1-(3-(3,4-dimethoxyphenyl)-2,6-dimethylimidazo[1,2-b]pyridazin-8-yl)pyrrolidin-3-yl)-carbamate (PDDC) and 2,6-dimethoxy-4-[4-phenyl-5-(2-thienyl)-1H-imidazol-2-yl]phenol (DPTIP)). PDDC and DPTIP show differential effects on cell viability, and EV marker secretion, indicating that side effects of these inhibitors on lysosomal and autophagic degradation pathways need to be considered. Moreover, similar to commonly used nSMase2 inhibitor GW4869, cell type specificity seems to play a role in the endosomal trafficking routes that can be explored to unravel mechanisms of specific EV biogenesis and secretion pathways.

CellCallEXT: Analysis of Ligand-Receptor and Transcription Factor Activities in Cell-Cell Communication of Tumor Immune Microenvironment.

Simple SummaryCellCall is an R package tool that is used to analyze cell–cell communication based on transcription factor (TF) activities calculated by cell-type specificity of target genes and thus cannot directly handle two-condition comparisons. We developed CellCallEXT to complement CellCall. CellCallEXT can directly identify ligand–receptor (L–R) interactions that alter the expression profiles of downstream genes between two conditions, such as tumor and healthy tissue. Scoring in CellCallEXT quantitatively integrates expression of ligands, receptors, TFs, and target genes (TGs). The pathway enrichment analysis and visualization modules allow biologists to investigate how disease alters cell–cell communication. Furthermore, Reactome pathways were added into CellCallEXT to expand the L–R–TF database.(1) Background: Single-cell RNA sequencing (scRNA-seq) data are useful for decoding cell–cell communication. CellCall is a tool that is used to infer inter- and intracellular communication pathways by integrating paired ligand–receptor (L–R) and transcription factor (TF) activities from steady-state data and thus cannot directly handle two-condition comparisons. For tumor and healthy status, it can only individually analyze cells from tumor or healthy tissue and examine L–R pairs only identified in either tumor or healthy controls, but not both together. Furthermore, CellCall is highly affected by gene expression specificity in tissues. (2) Methods: CellCallEXT is an extension of CellCall that deconvolutes intercellular communication and related internal regulatory signals based on scRNA-seq. Information on Reactome was retrieved and integrated with prior knowledge of L–R–TF signaling and gene regulation datasets of CellCall. (3) Results: CellCallEXT was successfully applied to examine tumors and immune cell microenvironments and to identify the altered L–R pairs and downstream gene regulatory networks among immune cells. Application of CellCallEXT to scRNA-seq data from patients with deficiency of adenosine deaminase 2 demonstrated its ability to impute dysfunctional intercellular communication and related transcriptional factor activities. (4) Conclusions: CellCallEXT provides a practical tool to examine intercellular communication in disease based on scRNA-seq data.

DLoopCaller: A deep learning approach for predicting genome-wide chromatin loops by integrating accessible chromatin landscapes.

In recent years, major advances have been made in various chromosome conformation capture technologies to further satisfy the needs of researchers for high-quality, high-resolution contact interactions. Discriminating the loops from genome-wide contact interactions is crucial for dissecting three-dimensional(3D) genome structure and function. Here, we present a deep learning method to predict genome-wide chromatin loops, called DLoopCaller, by combining accessible chromatin landscapes and raw Hi-C contact maps. Some available orthogonal data ChIA-PET/HiChIP and Capture Hi-C were used to generate positive samples with a wider contact matrix which provides the possibility to find more potential genome-wide chromatin loops. The experimental results demonstrate that DLoopCaller effectively improves the accuracy of predicting genome-wide chromatin loops compared to the state-of-the-art method Peakachu. Moreover, compared to two of most popular loop callers, such as HiCCUPS and Fit-Hi-C, DLoopCaller identifies some unique interactions. We conclude that a combination of chromatin landscapes on the one-dimensional genome contributes to understanding the 3D genome organization, and the identified chromatin loops reveal cell-type specificity and transcription factor motif co-enrichment across different cell lines and species.

Human cardiomyocytes are more susceptible to irreversible electroporation by pulsed electric field than human esophageal cells.

Pulse electric field‐based (PEF) ablation is a technique whereby short high‐intensity electric fields inducing irreversible electroporation (IRE) are applied to various tissues. Here, we implemented a standardized in vitro model to compare the effects of biphasic symmetrical pulses (100 pulses, 1–10 μs phase duration (d), 10–1000 Hz pulse repetition rate (f)) using two different human cellular models: human‐induced pluripotent stem cell‐derived cardiomyocytes (hiPSC‐CMs) and human esophageal smooth muscle cells (hESMCs) cultured in monolayer format. We report the PEF‐induced irreversibly electroporated cell monolayer areas and the corresponding electric field thresholds (EFTs) for both cardiac and esophageal cultures. Our results suggest marked cell type specificity with EFT estimated to be 2–2.5 times lower in hiPSC‐CMs than in hESMCs when subjected to identical PEF treatments (e.g., 0.90 vs 1.85 kV/cm for the treatment of 100 pulses with d = 5 μs, f = 10 Hz, and 0.65 vs 1.67 kV/cm for the treatment of 100 pulses with d = 10 μs, f = 10 Hz). PEF treatment can result in increased temperature around the stimulating electrodes and lead to unanticipated thermal tissue damage that is proportional to the peak temperature rise and to the duration of the PEF‐induced elevated temperatures. In our study, temperature increases ranged from less than 1°C to as high as 30°C, however, all temperature changes were transient and quickly returned to baseline and the highest observed ∆T returned to 50% of its maximum recorded temperature in tens of seconds.

Nucleosome proteostasis and histone turnover.

Maintenance of protein folding homeostasis, or proteostasis is critical for cell survival as well as for execution of cell type specific biological processes such as muscle cell contractility, neuronal synapse and memory formation, and cell transition from a mitotic to post-mitotic cell type. Cell type specification is driven largely by chromatin organization, which dictates which genes are turned off or on, depending on cell needs and function. Loss of chromatin organization can have catastrophic consequences either on cell survival or cell type specific function. Chromatin organization is highly dependent on organization of nucleosomes, spatiotemporal nucleosome assembly and disassembly, and histone turnover. In this review our goal is to highlight why nucleosome proteostasis is critical for chromatin organization, how this process is mediated by histone chaperones and ATP-dependent chromatin remodelers and outline potential and established mechanisms of disrupted nucleosome proteostasis during disease. Finally, we highlight how these mechanisms of histone turnover and nucleosome proteostasis may conspire with unfolded protein response programs to drive histone turnover in cell growth and development.

Expression of periodontitis susceptibility genes in human gingiva using single-cell RNA sequencing.

Single-cell transcriptomics was used to determine the possible cell-type specificity of periodontitis susceptibility genes. The last decade has witnessed remarkable advances in the field of human genomics. Despite many advances, the genetic factors associated with or contributing to the periodontitis pathogenesis have only been identified to a limited extent and are often poorly validated. Confirming whether a given single nucleotide polymorphism has an association with periodontitis requires a robust mechanistic explanation on the functional consequences of a given genetic variant. We globally assessed the expression of 26 disease-associated genes identified by GWAS within the gingival mucosa. A total of 12 411 cells from 4 different donors were analysed. Differentially expressed genes were analysed using Seurat, a non-parametric Wilcoxon rank sum test. The minimum threshold for significance was defined as p < .05. This exploration at a cellular-level suggests diverse populations contributing to disease pathogenesis, with macrophages expressing a higher number of the analysed disease-associated genes. IL1B, PTGS2, FCGR2A, IL10 and IL1A specifically showed a more restricted expression in the myeloid lineages. This short report combines human genetics and single-cell genomics to better understand periodontitis by mapping variants to predict their cells of action and putative functions. These findings seem to suggest that innate cell dysfunction is linked to disease susceptibility.

Identifying pleiotropic genes for major psychiatric disorders with GWAS summary statistics using multivariate adaptive association tests

BackgroundGenome wide association studies (GWAS) have discovered a few of single nucleotide polymorphisms (SNPs) related to major psychiatric disorders. However, it is not completely clear which genes play a pleiotropic role in multiple disorders. The study aimed to identify the pleiotropic genes across five psychiatric disorders using multivariate adaptive association tests. MethodsSummary statistics of five psychiatric disorders were downloaded from Psychiatric Genomics Consortium. We applied linkage disequilibrium score regression (LDSC) to estimate genetic correlation and conducted tissue and cell type specificity analyses based on Multi-marker Analysis of GenoMic Annotation (MAGMA). Then, we identified the pleiotropic genes using MTaSPUsSet and aSPUs tests. We ultimately performed the functional analysis for pleiotropic genes. ResultsWe confirmed the significant genetic correlation and brain tissue and neuron specificity among five disorders. 100 pleiotropic genes were detected to be significantly associated with five psychiatric disorders, of which 55 were novel genes. These genes were functionally enriched in neuron differentiation and synaptic transmission. LimitationsThe effect direction of pleiotropic genes couldn't be distinguished due to without individual-level data. ConclusionWe identified pleiotropic genes using multivariate adaptive association tests and explored their biological function. The findings may provide novel insight into the development and implementation of prevention and treatment as well as targeted drug discovery in practice.

Analysis of low-level somatic mosaicism reveals stage and tissue-specific mutational features in human development.

Most somatic mutations that arise during normal development are present at low levels in single or multiple tissues depending on the developmental stage and affected organs. However, the effect of human developmental stages or mutations of different organs on the features of somatic mutations is still unclear. Here, we performed a systemic and comprehensive analysis of low-level somatic mutations using deep whole-exome sequencing (average read depth ~500×) of 498 multiple organ tissues with matched controls from 190 individuals. Our results showed that early clone-forming mutations shared between multiple organs were lower in number but showed higher allele frequencies than late clone-forming mutations [0.54 vs. 5.83 variants per individual; 6.17% vs. 1.5% variant allele frequency (VAF)] along with less nonsynonymous mutations and lower functional impacts. Additionally, early and late clone-forming mutations had unique mutational signatures that were distinct from mutations that originated from tumors. Compared with early clone-forming mutations that showed a clock-like signature across all organs or tissues studied, late clone-forming mutations showed organ, tissue, and cell-type specificity in the mutation counts, VAFs, and mutational signatures. In particular, analysis of brain somatic mutations showed a bimodal occurrence and temporal-lobe-specific signature. These findings provide new insights into the features of somatic mosaicism that are dependent on developmental stage and brain regions.

De novo annotation of lncRNA HOTAIR transcripts by long-read RNA capture-seq reveals a differentiation-driven isoform switch

BackgroundLncRNAs are tissue-specific and emerge as important regulators of various biological processes and as disease biomarkers. HOTAIR is a well-established pro-oncogenic lncRNA which has been attributed a variety of functions in cancer and native contexts. However, a lack of an exhaustive, cell type-specific annotation questions whether HOTAIR functions are supported by the expression of multiple isoforms.ResultsUsing a capture long-read sequencing approach, we characterize HOTAIR isoforms expressed in human primary adipose stem cells. We find HOTAIR isoforms population displays varied splicing patterns, frequently leading to the exclusion or truncation of canonical LSD1 and PRC2 binding domains. We identify a highly cell type-specific HOTAIR isoform pool regulated by distinct promoter usage, and uncover a shift in the HOTAIR TSS usage that modulates the balance of HOTAIR isoforms at differentiation onset.ConclusionOur results highlight the complexity and cell type-specificity of HOTAIR isoforms and open perspectives on functional implications of these variants and their balance to key cellular processes.

WNT5A-RHOA Signaling Is a Driver of Tumorigenesis and Represents a Therapeutically Actionable Vulnerability in Small Cell Lung Cancer.

The p130-WNT5A-RHOA pathway drives SCLC progression and is a potential target for the development of therapeutic interventions and biomarkers to improve patient treatment.

Jiyuan oridonin A induces differentiation of acute myeloid leukemia cells including leukemic stem-like cells.

Acute myeloid leukemia (AML) is an aggressive form of hematological neoplasia characterized by failure of myeloid differentiation. AML is a leading cause of death from leukemia. Cytarabine chemotherapy resistance is a major source of refractory/relapsed AML. A major obstacle to the successful treatment of AML results from residual disease maintained by leukemic stem cells (LSCs), which are mostly resistant to conventional chemotherapy. Here, we determined the effect of a natural compound, Jiyuan oridonin A (JOA), on the differentiation blockade in the M2 subtype [particularly t (8;21)] of AML cells, M3 subtype of AML cells (APL cells), and leukemic stem-like cells both in vitro and in vivo. We found that JOA induced cell differentiation and suppressed the colony formation capacity in various AML cell lines (Kasumi-1, KG-1, MUTZ-8, NB4, and HL-60) without eliciting apoptosis. The mechanism of JOA-induced cell differentiation depends on the specificity of cell type. JOA mediated the differentiation of Kasumi-1 cells by activating the hematopoietic cell lineage signaling pathway, while inhibition of c-MYC was involved in the JOA-induced differentiation of NB4 cells. Moreover, JOA was identified to target leukemic stem-like cells by induced cell differentiation in vivo. These findings demonstrated that JOA could inhibit the proliferation of M2 and M3 subtypes of AML cells and leukemic stem-like cells by overcoming the differentiation blockade, which may offer a novel therapeutic strategy for AML to overcome relapse and drug resistance in patients with AML. Our findings highlight the possibility of using compounds like JOA as a promising differentiation-induced agent for the treatment of AML.

Transition dynamics and optogenetic control of epileptiform activity in a modified mean filed model of human cortex

Physiological experiments have shown that inhibitory interneurons can induce and maintain epileptiform activities, and different interneuron subtypes may be responsible for different types of seizures. Here we aim to link these electrophysiological experimental phenomena with theoretical dynamic mechanisms based on an improved Liley mean field model. Fascinatingly, the number of synapses between inhibitory neural populations can induce a rich state transition. Typically, the system will experience normal rhythm discharges, multispike discharges, spike wave discharges (SWD), generalized periodic discharges (GPD), the beta band oscillation, and eventually return to the normal state. Interestingly, the transition process can also be reversed. Meanwhile, disinhibition circuits can cause more epileptiform activities after taking into account the delay of synaptic information transmission. Furthermore, we are committed to designing different control strategies for epileptiform activities. As we expected, both deep brain stimulation and optogenetic technology can destroy or even eliminate pathological waves. It is need to emphasis that the cell type specificity of optogenetic regulation allows it to precisely target inhibitory neural populations, which is agree with experiment. Particularly, three different optogenetic regulatory strategies targeting inhibitory neural populations are modeled and proposed, of which the intermittent is designed to save energy. These modeling results reproduce the experimental phenomena and more significantly, help understand the mechanism of epilepsy to guide the clinical practice.

A Human Stem Cell-Derived Neurosensory-Epithelial Circuitry on a Chip to Model Herpes Simplex Virus Reactivation.

Both emerging viruses and well-known viral pathogens endowed with neurotropism can either directly impair neuronal functions or induce physio-pathological changes by diffusing from the periphery through neurosensory–epithelial connections. However, developing a reliable and reproducible in vitro system modeling the connectivity between the different human sensory neurons and peripheral tissues is still a challenge and precludes the deepest comprehension of viral latency and reactivation at the cellular and molecular levels. This study shows a stable topographic neurosensory–epithelial connection on a chip using human stem cell-derived dorsal root ganglia (DRG) organoids. Bulk and single-cell transcriptomics showed that different combinations of key receptors for herpes simplex virus 1 (HSV-1) are expressed by each sensory neuronal cell type. This neuronal–epithelial circuitry enabled a detailed analysis of HSV infectivity, faithfully modeling its dynamics and cell type specificity. The reconstitution of an organized connectivity between human sensory neurons and keratinocytes into microfluidic chips provides a powerful in vitro platform for modeling viral latency and reactivation of human viral pathogens.

Sustained deep-tissue voltage recording using a fast indicator evolved for two-photon microscopy

SUMMARYGenetically encoded voltage indicators are emerging tools for monitoring voltage dynamics with cell-type specificity. However, current indicators enable a narrow range of applications due to poor performance under two-photon microscopy, a method of choice for deep-tissue recording. To improve indicators, we developed a multiparameter high-throughput platform to optimize voltage indicators for two-photon microscopy. Using this system, we identified JEDI-2P, an indicator that is faster, brighter, and more sensitive and photostable than its predecessors. We demonstrate that JEDI-2P can report light-evoked responses in axonal termini of Drosophila interneurons and the dendrites and somata of amacrine cells of isolated mouse retina. JEDI-2P can also optically record the voltage dynamics of individual cortical neurons in awake behaving mice for more than 30 min using both resonant-scanning and ULoVE random-access microscopy. Finally, ULoVE recording of JEDI-2P can robustly detect spikes at depths exceeding 400 μm and report voltage correlations in pairs of neurons.

Modeling chromatin state from sequence across angiosperms using recurrent convolutional neural networks.

Accessible chromatin regions are critical components of gene regulation but modeling them directly from sequence remains challenging, especially within plants, whose mechanisms of chromatin remodeling are less understood than in animals. We trained an existing deep-learning architecture, DanQ, on data from 12 angiosperm species to predict the chromatin accessibility in leaf of sequence windows within and across species. We also trained DanQ on DNA methylation data from 10 angiosperms because unmethylated regions have been shown to overlap significantly with ACRs in some plants. The across-species models have comparable or even superior performance to a model trained within species, suggesting strong conservation of chromatin mechanisms across angiosperms. Testing a maize (Zea mays L.) held-out model on a multi-tissue chromatin accessibility panel revealed our models are best at predicting constitutively accessible chromatin regions, with diminishing performance as cell-type specificity increases. Using a combination of interpretation methods, we ranked JASPAR motifs by their importance to each model and saw that the TCP and AP2/ERF transcription factor (TF) families consistently ranked highly. We embedded the top three JASPAR motifs for each model at all possible positions on both strands in our sequence window and observed position- and strand-specific patterns in their importance to the model. With our publicly available across-species 'a2z' model it is now feasible to predict the chromatin accessibility and methylation landscape of any angiosperm genome.

All-optical inter-layers functional connectivity investigation in the mouse retina

SummaryWe developed a multi-unit microscope for all-optical inter-layers circuits interrogation. The system performs two-photon (2P) functional imaging and 2P multiplexed holographic optogenetics at axially distinct planes. We demonstrated the capability of the system to map, in the mouse retina, the functional connectivity between rod bipolar cells (RBCs) and ganglion cells (GCs) by activating single or defined groups of RBCs while recording the evoked response in the GC layer with cell-type specificity and single-cell resolution. We then used a logistic model to probe the functional connectivity between cell types by deriving the “cellular receptive field” describing how RBCs impact each GC type. With the capability to simultaneously image and control neuronal activity at axially distinct planes, the system enables a precise interrogation of multi-layered circuits. Understanding this information transfer is a promising avenue to dissect complex neural circuits and understand the neural basis of computations.