Single-Cell Co-expression Analysis Reveals Distinct Functional Modules, Co-regulation Mechanisms and Clinical Outcomes

Tracing the Shared Foundations of Gene Expression and Chromatin Structure

The three-dimensional organization of chromatin into topologically associating domains (TADs) may impact gene regulation by bringing distant genes into contact. However, many questions about TADs' function and their influence on transcription remain unresolved due to technical limitations in defining TAD boundaries and measuring the direct effect that TADs have on gene expression. Here, we develop consensus TAD maps for human and mouse with a novel "bag-of-genes" approach for defining the gene composition within TADs. This approach enables new functional interpretations of TADs by providing a way to capture species-level differences in chromatin organization. We also leverage a generative AI foundation model computed from 33 million transcriptomes to define contextual similarity, an embedding-based metric that is more powerful than co-expression at representing functional gene relationships. Our analytical framework directly leads to testable hypotheses about chromatin organization across cellular states. We find that TADs play an active role in facilitating gene co-regulation, possibly through a mechanism involving transcriptional condensates. We also discover that the TAD-linked enhancement of transcriptional context is strongest in early developmental stages and systematically declines with aging. Investigation of cancer cells show distinct patterns of TAD usage that shift with chemotherapy treatment, suggesting specific roles for TAD-mediated regulation in cellular development and plasticity. Finally, we develop "TAD signatures" to improve statistical analysis of single-cell transcriptomic data sets in predicting cancer cell-line drug response. These findings reshape our understanding of cellular plasticity in development and disease, indicating that chromatin organization acts through probabilistic mechanisms rather than deterministic rules.

Software availability

https://singhlab.net/tadmap.

Leveraging cell type-specificity for gene set analysis of single cell transcriptomics

Although single cell RNA-sequencing (scRNA-seq) provides unprecedented insights into the biology of complex tissues, analyzing such data on a gene-by-gene basis is challenging due to the large number of tested hypotheses and consequent low statistical power and difficult interpretation. These issues are magnified by the increased noise, significant sparsity and multi-modal distributions characteristic of single cell data. One promising approach for addressing these challenges is gene set testing, or pathway analysis. Unfortunately, statistical and biological differences between single cell and bulk transcriptomic data make it challenging to use existing gene set collections, which were developed for bulk tissue analysis, on scRNA-seq data. In this paper, we describe a procedure for customizing gene set collections originally created for bulk tissue analysis to reflect the structure of gene activity within specific cell types. Our approach leverages information about mean gene expression in the 81 human cell types profiled via scRNA-seq by the Human Protein Atlas (HPA) Single Cell Type Atlas. This HPA information is used to compute cell type-specific gene and gene set weights that can be used to filter or weight gene set collections. As demonstrated through the analysis of immune cell scRNA-seq data using gene sets from the Molecular Signatures Database (MSigDB), accounting for cell type-specificity can significantly improve gene set testing power and interpretability. An example vignette along with gene and gene set weights for the 81 HPA SCTA cell types and the MSigDB collections are available at https://hrfrost.host.dartmouth.edu/SCGeneSetOpt/.

Comparative proteomic analysis across the developmental stages of the Eimeria tenella

Eimeria tenella is the main pathogen responsible for coccidiosis in chickens. The life cycle of E. tenella is, arguably, the least complex of all Coccidia, with only one host. However, it presents different developmental stages, either in the environment or in the host and either intracellular or extracellular. Its signaling and metabolic pathways change with its different developmental stages. Until now, little is known about the developmental regulation and transformation mechanisms of its life cycle. In this study, protein profiles from the five developmental stages, including unsporulated oocysts (USO), partially sporulated (7 h) oocysts (SO7h), sporulated oocysts (SO), sporozoites (S) and second-generation merozoites (M2), were harvested using the label-free quantitative proteomics approach. Then the differentially expressed proteins (DEPs) for these stages were identified. A total of 314, 432, 689, and 665 DEPs were identified from the comparison of SO7h vs USO, SO vs SO7h, S vs SO, and M2 vs S, respectively. By conducting weighted gene coexpression network analysis (WGCNA), six modules were dissected. Proteins in blue and brown modules were calculated to be significantly positively correlated with the E. tenella developmental stages of sporozoites (S) and second-generation merozoites (M2), respectively. In addition, hub proteins with high intra-module degree were identified. Gene Ontology (GO) and Kyoto Encyclopedia of Gene and Genomes (KEGG) pathway enrichment analyses revealed that hub proteins in blue modules were involved in electron transport chain and oxidative phosphorylation. Hub proteins in the brown module were involved in RNA splicing. These findings provide new clues and ideas to enhance our fundamental understanding of the molecular mechanisms underlying parasite development.

Cell type evolution reconstruction across species through cell phylogenies of single-cell RNA sequencing data

The origin and evolution of cell types has emerged as a key topic in evolutionary biology. Driven by rapidly accumulating single-cell datasets, recent attempts to infer cell type evolution have largely been limited to pairwise comparisons because we lack approaches to build cell phylogenies using model-based approaches. Here we approach the challenges of applying explicit phylogenetic methods to single-cell data by using principal components as phylogenetic characters. We infer a cell phylogeny from a large, comparative single-cell dataset of eye cells from five distantly related mammals. Robust cell type clades enable us to provide a phylogenetic, rather than phenetic, definition of cell type, allowing us to forgo marker genes and phylogenetically classify cells by topology. We further observe evolutionary relationships between diverse vessel endothelia and identify the myelinating and non-myelinating Schwann cells as sister cell types. Finally, we examine principal component loadings and describe the gene expression dynamics underlying the function and identity of cell type clades that have been conserved across the five species. A cell phylogeny provides a rigorous framework towards investigating the evolutionary history of cells and will be critical to interpret comparative single-cell datasets that aim to ask fundamental evolutionary questions.

Identification of genetic variants that impact gene co-expression relationships using large-scale single-cell data

BACKGROUND

Expression quantitative trait loci (eQTL) studies show how genetic variants affect downstream gene expression. Single-cell data allows reconstruction of personalized co-expression networks and therefore the identification of SNPs altering co-expression patterns (co-expression QTLs, co-eQTLs) and the affected upstream regulatory processes using a limited number of individuals.

RESULTS

We conduct a co-eQTL meta-analysis across four scRNA-seq peripheral blood mononuclear cell datasets using a novel filtering strategy followed by a permutation-based multiple testing approach. Before the analysis, we evaluate the co-expression patterns required for co-eQTL identification using different external resources. We identify a robust set of cell-type-specific co-eQTLs for 72 independent SNPs affecting 946 gene pairs. These co-eQTLs are replicated in a large bulk cohort and provide novel insights into how disease-associated variants alter regulatory networks. One co-eQTL SNP, rs1131017, that is associated with several autoimmune diseases, affects the co-expression of RPS26 with other ribosomal genes. Interestingly, specifically in T cells, the SNP additionally affects co-expression of RPS26 and a group of genes associated with T cell activation and autoimmune disease. Among these genes, we identify enrichment for targets of five T-cell-activation-related transcription factors whose binding sites harbor rs1131017. This reveals a previously overlooked process and pinpoints potential regulators that could explain the association of rs1131017 with autoimmune diseases.

CONCLUSION

Our co-eQTL results highlight the importance of studying context-specific gene regulation to understand the biological implications of genetic variation. With the expected growth of sc-eQTL datasets, our strategy and technical guidelines will facilitate future co-eQTL identification, further elucidating unknown disease mechanisms.

siVAE: interpretable deep generative models for single-cell transcriptomes

Neural networks such as variational autoencoders (VAE) perform dimensionality reduction for the visualization and analysis of genomic data, but are limited in their interpretability: it is unknown which data features are represented by each embedding dimension. We present siVAE, a VAE that is interpretable by design, thereby enhancing downstream analysis tasks. Through interpretation, siVAE also identifies gene modules and hubs without explicit gene network inference. We use siVAE to identify gene modules whose connectivity is associated with diverse phenotypes such as iPSC neuronal differentiation efficiency and dementia, showcasing the wide applicability of interpretable generative models for genomic data analysis.

A multimodal deep learning model to infer cell-type-specific functional gene networks

BACKGROUND

Functional gene networks (FGNs) capture functional relationships among genes that vary across tissues and cell types. Construction of cell-type-specific FGNs enables the understanding of cell-type-specific functional gene relationships and insights into genetic mechanisms of human diseases in disease-relevant cell types. However, most existing FGNs were developed without consideration of specific cell types within tissues.

RESULTS

In this study, we created a multimodal deep learning model (MDLCN) to predict cell-type-specific FGNs in the human brain by integrating single-nuclei gene expression data with global protein interaction networks. We systematically evaluated the prediction performance of the MDLCN and showed its superior performance compared to two baseline models (boosting tree and convolutional neural network). Based on the predicted cell-type-specific FGNs, we observed that cell-type marker genes had a higher level of hubness than non-marker genes in their corresponding cell type. Furthermore, we showed that risk genes underlying autism and Alzheimer's disease were more strongly connected in disease-relevant cell types, supporting the cellular context of predicted cell-type-specific FGNs.

CONCLUSIONS

Our study proposes a powerful deep learning approach (MDLCN) to predict FGNs underlying a diverse set of cell types in human brain. The MDLCN model enhances prediction accuracy of cell-type-specific FGNs compared to single modality convolutional neural network (CNN) and boosting tree models, as shown by higher areas under both receiver operating characteristic (ROC) and precision-recall curves for different levels of independent test datasets. The predicted FGNs also show evidence for the cellular context and distinct topological features (i.e. higher hubness and topological score) of cell-type marker genes. Moreover, we observed stronger modularity among disease-associated risk genes in FGNs of disease-relevant cell types. For example, the strength of connectivity among autism risk genes was stronger in neurons, but risk genes underlying Alzheimer's disease were more connected in microglia.

Weighted Gene Co-Expression Network Analysis (WGCNA) Discovered Novel Long Non-Coding RNAs for Polycystic Ovary Syndrome

Polycystic ovary syndrome (PCOS) affects reproductive-age women. This condition causes infertility, insulin resistance, obesity, and heart difficulties. The molecular basis and mechanism of PCOS might potentially generate effective treatments. Long non-coding RNAs (lncRNAs) show control over multifactorial disorders' growth and incidence. Numerous studies have emphasized its significance and alterations in PCOS. We used bioinformatic methods to find novel dysregulated lncRNAs in PCOS. To achieve this objective, the gene expression profile of GSE48301, comprising PCOS patients and normal control tissue samples, was evaluated using the R limma package with the following cut-off criterion: p-value < 0.05. Firstly, weighted gene co-expression network analysis (WGCNA) was used to determine the co-expression genes of lncRNAs; subsequently, hub gene identification and pathway enrichment analysis were used. With the defined criteria, nine novel dysregulated lncRNAs were identified. In WGCNA, different colors represent different modules. In the current study, WGCNA resulted in turquoise, gray, blue, and black co-expression modules with dysregulated lncRNAs. The pathway enrichment analysis of these co-expressed modules revealed enrichment in PCOS-associated pathways, including gene expression, signal transduction, metabolism, and apoptosis. In addition, CCT7, EFTUD2, ESR1, JUN, NDUFAB1, CTTNB1, GRB2, and CTNNB1 were identified as hub genes, and some of them have been investigated in PCOS. This study uncovered nine novel PCOS-related lncRNAs. To confirm how these lncRNAs control translational modification in PCOS, functional studies are required.

Gut Microbiota Alterations in Trace Amine-Associated Receptor 9 (TAAR9) Knockout Rats

Trace amine-associated receptors (TAAR1-TAAR9) are a family of G-protein-coupled monoaminergic receptors which might have great pharmacological potential. It has now been well established that TAAR1 plays an important role in the central nervous system. Interestingly, deletion of TAAR9 in rats leads to alterations in the periphery. Previously, we found that knockout of TAAR9 in rats (TAAR9-KO rats) decreased low-density lipoprotein cholesterol levels in the blood. TAAR9 was also identified in intestinal tissues, and it is known that it responds to polyamines. To elucidate the role of TAAR9 in the intestinal epithelium, we analyzed TAAR9-co-expressed gene clusters in public data for cecum samples. As identified by gene ontology enrichment analysis, in the intestine, TAAR9 is co-expressed with genes involved in intestinal mucosa homeostasis and function, including cell organization, differentiation, and death. Additionally, TAAR9 was co-expressed with genes implicated in dopamine signaling, which may suggest a role for this receptor in the regulation of peripheral dopaminergic transmission. To further investigate how TAAR9 might be involved in colonic mucosal homeostasis, we analyzed the fecal microbiome composition in TAAR9-KO rats and their wild-type littermates. We identified a significant difference in the number of observed taxa between the microbiome of TAAR9-KO and wild-type rats. In TAAR9-KO rats, the gut microbial community became more variable compared with the wild-type rats. Furthermore, it was found that the family Saccharimonadaceae, which is one of the top 10 most abundant families in TAAR9-KO rat feces, is almost completely absent in wild-type animal fecal samples. Taken together, these data indicate a role of TAAR9 in intestinal function.

SpaceX: gene co-expression network estimation for spatial transcriptomics

MOTIVATION

The analysis of spatially resolved transcriptome enables the understanding of the spatial interactions between the cellular environment and transcriptional regulation. In particular, the characterization of the gene-gene co-expression at distinct spatial locations or cell types in the tissue enables delineation of spatial co-regulatory patterns as opposed to standard differential single gene analyses. To enhance the ability and potential of spatial transcriptomics technologies to drive biological discovery, we develop a statistical framework to detect gene co-expression patterns in a spatially structured tissue consisting of different clusters in the form of cell classes or tissue domains.

RESULTS

We develop SpaceX (spatially dependent gene co-expression network), a Bayesian methodology to identify both shared and cluster-specific co-expression network across genes. SpaceX uses an over-dispersed spatial Poisson model coupled with a high-dimensional factor model which is based on a dimension reduction technique for computational efficiency. We show via simulations, accuracy gains in co-expression network estimation and structure by accounting for (increasing) spatial correlation and appropriate noise distributions. In-depth analysis of two spatial transcriptomics datasets in mouse hypothalamus and human breast cancer using SpaceX, detected multiple hub genes which are related to cognitive abilities for the hypothalamus data and multiple cancer genes (e.g. collagen family) from the tumor region for the breast cancer data.

AVAILABILITY AND IMPLEMENTATION

The SpaceX R-package is available at github.com/bayesrx/SpaceX.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

scGENA: A Single-Cell Gene Coexpression Network Analysis Framework for Clustering Cell Types and Revealing Biological Mechanisms

Single-cell RNA-sequencing (scRNA-seq) is a recent high-throughput technique that can measure gene expression, reveal cell heterogeneity, rare and complex cell populations, and discover cell types and their relationships. The analysis of scRNA-seq data is challenging because of transcripts sparsity, replication noise, and outlier cell populations. A gene coexpression network (GCN) analysis effectively deciphers phenotypic differences in specific states by describing gene–gene pairwise relationships. The underlying gene modules with different coexpression patterns partially bridge the gap between genotype and phenotype. This study presents a new framework called scGENA (single-cell gene coexpression network analysis) for GCN analysis based on scRNA-seq data. Although there are several methods for scRNA-seq data analysis, we aim to build an integrative pipeline for several purposes that cover primary data preprocessing, including data exploration, quality control, normalization, imputation, and dimensionality reduction of clustering as downstream of GCN analysis. To demonstrate this integrated workflow, an scRNA-seq dataset of the human diabetic pancreas with 1600 cells and 39,851 genes was implemented to perform all these processes in practice. As a result, scGENA is demonstrated to uncover interesting gene modules behind complex diseases, which reveal biological mechanisms. scGENA provides a state-of-the-art method for gene coexpression analysis for scRNA-seq data.

Deeper insights into long-term survival heterogeneity of pancreatic ductal adenocarcinoma (PDAC) patients using integrative individual- and group-level transcriptome network analyses

Pancreatic ductal adenocarcinoma (PDAC) is categorized as the leading cause of cancer mortality worldwide. However, its predictive markers for long-term survival are not well known. It is interesting to delineate individual-specific perturbed genes when comparing long-term (LT) and short-term (ST) PDAC survivors and integrate individual- and group-based transcriptome profiling. Using a discovery cohort of 19 PDAC patients from CHU-Liège (Belgium), we first performed differential gene expression analysis comparing LT to ST survivor. Second, we adopted systems biology approaches to obtain clinically relevant gene modules. Third, we created individual-specific perturbation profiles. Furthermore, we used Degree-Aware disease gene prioritizing (DADA) method to develop PDAC disease modules; Network-based Integration of Multi-omics Data (NetICS) to integrate group-based and individual-specific perturbed genes in relation to PDAC LT survival. We identified 173 differentially expressed genes (DEGs) in ST and LT survivors and five modules (including 38 DEGs) showing associations to clinical traits. Validation of DEGs in the molecular lab suggested a role of REG4 and TSPAN8 in PDAC survival. Via NetICS and DADA, we identified various known oncogenes such as CUL1 and TGFB1. Our proposed analytic workflow shows the advantages of combining clinical and omics data as well as individual- and group-level transcriptome profiling.

Risk subtyping and prognostic assessment of prostate cancer based on consensus genes

Prostate cancer (PCa) is the most frequent malignancy in male urogenital system around worldwide. We performed molecular subtyping and prognostic assessment based on consensus genes in patients with PCa. Five cohorts containing 1,046 PCa patients with RNA expression profiles and recorded clinical follow-up information were included. Univariate, multivariate Cox regression analysis and least absolute shrinkage and selection operator (LASSO) Cox regression were used to select prognostic genes and establish the signature. Immunohistochemistry staining, cell proliferation, migration and invasion assays were used to assess the biological functions of key genes. Thirty-nine intersecting consensus prognostic genes from five independent cohorts were identified. Subsequently, an eleven-consensus-gene classifier was established. In addition, multivariate Cox regression analyses showed that the classifier served as an independent indicator of recurrence-free survival in three of the five cohorts. Combined receiver operating characteristic (ROC) analysis achieved synthesized effects by combining the classifier with clinicopathological features in four of five cohorts. SRD5A2 inhibits cell proliferation, while ITGA11 promotes cell migration and invasion, possibly through the PI3K/AKT signaling pathway. To conclude, we established and validated an eleven-consensus-gene classifier, which may add prognostic value to the currently available staging system.

By analysis of gene expression profiles of prostate cancer patients from multiple platforms, an eleven-consensus-gene classifier is constructed to provide a robust tool for the prediction of recurrence-free survival.

JEBIN: analyzing gene co-expressions across multiple datasets by joint network embedding

The inference of gene co-expression associations is one of the fundamental tasks for large-scale transcriptomic data analysis. Due to the high dimensionality and high noises in transcriptomic data, it is difficult to infer stable gene co-expression associations from single dataset. Meta-analysis of multisource data can effectively tackle this problem. We proposed Joint Embedding of multiple BIpartite Networks (JEBIN) to learn the low-dimensional consensus representation for genes by integrating multiple expression datasets. JEBIN infers gene co-expression associations in a nonlinear and global similarity manner and can integrate datasets with different distributions in linear time complexity with the gene and total sample size. The effectiveness and scalability of JEBIN were verified by simulation experiments, and its superiority over the commonly used integration methods was proved by three indexes on real biological datasets. Then, JEBIN was applied to study the gene co-expression patterns of hepatocellular carcinoma (HCC) based on multiple expression datasets of HCC and adjacent normal tissues, and further on latest HCC single-cell RNA-seq data. Results show that gene co-expressions are highly different between bulk and single-cell datasets. Finally, many differentially co-expressed ligand-receptor pairs were discovered by comparing HCC with adjacent normal data, providing candidate HCC targets for abnormal cell-cell communications.

Linking co-expression modules with phenotypes

The method for quantifying the association between co-expression module and clinical trait of interest requires application of dimensionality reduction to summaries modules as one dimensional (1D) vector. However, these methods are often linked with information loss. The amount of information lost depends upon the percentage of variance captured by the reduced 1D vector. Therefore, it is of interest to describe a method using analysis of rank (AOR) to assess the association between module and clinical trait of interest. This method works with clinical traits represented as binary class labels and can be adopted for clinical traits measured in continuous scale by dividing samples in two groups around median value. Application of the AOR method on test data for muscle gene expression profiles identifies modules significantly associated with diabetes status.

Distinct characteristics of correlation analysis at the single-cell and the population level

Correlation analysis is widely used in biological studies to infer molecular relationships within biological networks. Recently, single-cell analysis has drawn tremendous interests, for its ability to obtain high-resolution molecular phenotypes. It turns out that there is little overlap of co-expressed genes identified in single-cell level investigations with that of population level investigations. However, the nature of the relationship of correlations between single-cell and population levels remains unclear. In this manuscript, we aimed to unveil the origin of the differences between the correlation coefficients at the single-cell level and that at the population level, and bridge the gap between them. Through developing formulations to link correlations at the single-cell and the population level, we illustrated that aggregated correlations could be stronger, weaker or equal to the corresponding individual correlations, depending on the variations and the correlations within the population. When the correlation within the population is weaker than the individual correlation, the aggregated correlation is stronger than the corresponding individual correlation. Besides, our data indicated that aggregated correlation is more likely to be stronger than the corresponding individual correlation, and it was rare to find gene-pairs exclusively strongly correlated at the single-cell level. Through a bottom-up approach to model interactions between molecules in a signaling cascade or a multi-regulator-controlled gene expression, we surprisingly found that the existence of interaction between two components could not be excluded simply based on their low correlation coefficients, suggesting a reconsideration of connectivity within biological networks which was derived solely from correlation analysis. We also investigated the impact of technical random measurement errors on the correlation coefficients for the single-cell level and the population level. The results indicate that the aggregated correlation is relatively robust and less affected. Because of the heterogeneity among single cells, correlation coefficients calculated based on data of the single-cell level might be different from that of the population level. Depending on the specific question we are asking, proper sampling and normalization procedure should be done before we draw any conclusions.

Single-cell transcriptome atlas of human mesenchymal stem cells exploring cellular heterogeneity

BACKGROUND

The heterogeneity of mesenchymal stem cells (MSCs) is poorly understood, thus limiting clinical application and basic research reproducibility. Advanced single-cell RNA sequencing (scRNA-seq) is a robust tool used to analyse for dissecting cellular heterogeneity. However, the comprehensive single-cell atlas for human MSCs has not been achieved.

METHODS

This study used massive parallel multiplexing scRNA-seq to construct an atlas of > 130 000 single-MSC transcriptomes across multiple tissues and donors to assess their heterogeneity. The most widely clinically utilised tissue resources for MSCs were collected, including normal bone marrow (n = 3), adipose (n = 3), umbilical cord (n = 2), and dermis (n = 3).

RESULTS

Seven tissue-specific and five conserved MSC subpopulations with distinct gene-expression signatures were identified from multiple tissue origins based on the high-quality data, which has not been achieved previously. This study showed that extracellular matrix (ECM) highly contributes to MSC heterogeneity. Notably, tissue-specific MSC subpopulations were substantially heterogeneous on ECM-associated immune regulation, antigen processing/presentation, and senescence, thus promoting inter-donor and intra-tissue heterogeneity. The variable dynamics of ECM-associated genes had discrete trajectory patterns across multiple tissues. Additionally, the conserved and tissue-specific transcriptomic-regulons and protein-protein interactions were identified, potentially representing common or tissue-specific MSC functional roles. Furthermore, the umbilical-cord-specific subpopulation possessed advantages in immunosuppressive properties.

CONCLUSION

In summary, this work provides timely and great insights into MSC heterogeneity at multiple levels. This MSC atlas taxonomy also provides a comprehensive understanding of cellular heterogeneity, thus revealing the potential improvements in MSC-based therapeutic efficacy.

Evolving Transcriptomic Profiles From Single-Cell RNA-Seq Data Using Nature-Inspired Multiobjective Optimization

Transcriptomic profiling plays an important role in post-genomic analysis. Especially, the single-cell RNA-seq technology has advanced our understanding of gene expression from cell population level into individual cell level. Many computational methods have been proposed to decipher transcriptomic profiles from those RNA-seq data. However, most of the related algorithms suffer from realistic restrictions such as high dimensionality and premature convergence. In this paper, we propose and formulate an evolutionary multiobjective blind compressed sensing (EMOBCS) to address those problems for evolving transcriptomic profiles from single-cell RNA-seq data. In the proposed framework, to characterize various gene expression profile models, two objective functions including chi-squared kernel score and euclidean distance of different gene expression profiles are formulated. After that, multiobjective blind compressed sensing based on artificial bee colony is designed to optimize the two objective functions on single-cell RNA-seq data by proposing a rank probability model and two new search strategies into the cooperative convolution framework in an unbiased manner. To demonstrate its effectiveness, extensive experiments have been conducted, comparing the proposed algorithm with 14 algorithms including eight state-of-the-art algorithms and six different EMOBCS algorithms under different search strategies on 10 single-cell RNA-seq datasets and one case study. The experimental results reveal that the proposed algorithm is better than or comparable with those compared algorithms. Furthermore, we also conduct the time complexity analysis, convergence analysis, and parameter analysis to demonstrate various properties of EMOBCS.

Comprehensive Characterization of Multitissue Expression Landscape, Co-Expression Networks and Positive Selection in Pikeperch

Promising efforts are ongoing to extend genomics resources for pikeperch (Sander lucioperca), a species of high interest for the sustainable European aquaculture sector. Although previous work, including reference genome assembly, transcriptome sequence, and single-nucleotide polymorphism genotyping, added a great wealth of genomic tools, a comprehensive characterization of gene expression across major tissues in pikeperch still remains an unmet research need. Here, we used deep RNA-Sequencing of ten vital tissues collected in eight animals to build a high-confident and annotated trancriptome atlas, to detect the tissue-specificity of gene expression and co-expression network modules, and to investigate genome-wide selective signatures in the Percidae fish family. Pathway enrichment and protein–protein interaction network analyses were performed to characterize the unique biological functions of tissue-specific genes and co-expression modules. We detected strong functional correlations and similarities of tissues with respect to their expression patterns—but also significant differences in the complexity and composition of their transcriptomes. Moreover, functional analyses revealed that tissue-specific genes essentially play key roles in the specific physiological functions of the respective tissues. Identified network modules were also functionally coherent with tissues’ main physiological functions. Although tissue specificity was not associated with positive selection, several genes under selection were found to be involved in hypoxia, immunity, and gene regulation processes, that are crucial for fish adaption and welfare. Overall, these new resources and insights will not only enhance the understanding of mechanisms of organ biology in pikeperch, but also complement the amount of genomic resources for this commercial species.

Mutant IDH1 promotes phagocytic function of microglia/macrophages in gliomas by downregulating ICAM1

Tumor-associated microglia/macrophages (TAMs) are the main innate immune effector cells in malignant gliomas and have both pro- and anti-tumor functions. The plasticity of TAMs is partially dictated by oncogenic mutations in tumor cells. Heterozygous IDH1 mutation is a cancer driver gene prevalent in grade II/III gliomas, and IDH1 mutant gliomas have relatively favorable clinical outcomes. It is largely unknown how IDH mutation alters TAM phenotypes to influence glioma growth. Here we established clinically relevant isogenic glioma models carrying monoallelic IDH1 R132H mutation (IDH1^R132H/WT) and found that IDH1^R132H/WT significantly downregulated immune response-related pathways in glioma cells, indicating an immunomodulation role of mutant IDH1. Co-culturing IDH1^R132H/WT glioma cells with human macrophages promoted anti-tumor phenotypes of macrophages and increased macrophage migration and phagocytic capacity. In orthotopic xenografts, IDH1^R132H/WT decreased tumor growth and prolonged animal survival, accompanied by increased TAM recruitment and upregulated phagocytosis markers, suggesting the induction of anti-tumor TAM functions. Using human cytokine arrays that query 36 proteins, we identified significant downregulation of ICAM-1/CD54 in IDH1^R132H/WT gliomas, which was further confirmed by ELISA and immunoblotting analyses. ICAM1 gain-of-function studies revealed that ICAM1 downregulation in IDH1^R132H/WT cells played a mechanistic role to mediate the immunomodulation function of IDH1^R132H/WT. ICAM-1 silencing in IDH1 wild-type glioma cells decreased tumor growth and increased the anti-tumor function of TAMs. Together, our studies support a new TAM-mediated phagocytic function within IDH1 mutant gliomas, and improved understanding of this process may uncover novel approaches to targeting IDH1 wild type gliomas.

Comprehensive network modeling from single cell RNA sequencing of human and mouse reveals well conserved transcription regulation of hematopoiesis

Background

Presently, there is no comprehensive analysis of the transcription regulation network in hematopoiesis. Comparison of networks arising from gene co-expression across species can facilitate an understanding of the conservation of functional gene modules in hematopoiesis.

Results

We used single-cell RNA sequencing to profile bone marrow from human and mouse, and inferred transcription regulatory networks in each species in order to characterize transcriptional programs governing hematopoietic stem cell differentiation. We designed an algorithm for network reconstruction to conduct comparative transcriptomic analysis of hematopoietic gene co-expression and transcription regulation in human and mouse bone marrow cells. Co-expression network connectivity of hematopoiesis-related genes was found to be well conserved between mouse and human. The co-expression network showed “small-world” and “scale-free” architecture. The gene regulatory network formed a hierarchical structure, and hematopoiesis transcription factors localized to the hierarchy’s middle level.

Conclusions

Transcriptional regulatory networks are well conserved between human and mouse. The hierarchical organization of transcription factors may provide insights into hematopoietic cell lineage commitment, and to signal processing, cell survival and disease initiation.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-020-07241-2.

Nuclear gene proximity and protein interactions shape transcript covariations in mammalian single cells

Single-cell RNA sequencing studies on gene co-expression patterns could yield important regulatory and functional insights, but have so far been limited by the confounding effects of differentiation and cell cycle. We apply a tailored experimental design that eliminates these confounders, and report thousands of intrinsically covarying gene pairs in mouse embryonic stem cells. These covariations form a network with biological properties, outlining known and novel gene interactions. We provide the first evidence that miRNAs naturally induce transcriptome-wide covariations and compare the relative importance of nuclear organization, transcriptional and post-transcriptional regulation in defining covariations. We find that nuclear organization has the greatest impact, and that genes encoding for physically interacting proteins specifically tend to covary, suggesting importance for protein complex formation. Our results lend support to the concept of post-transcriptional RNA operons, but we further present evidence that nuclear proximity of genes may provide substantial functional regulation in mammalian single cells.

Gene expression covariation can be studied by single-cell RNA sequencing. Here the authors analyze intrinsically covarying gene pairs by eliminating the confounding effects in single-cell experiments and observe covariation of proximal genes and miRNA-induced covariation of target mRNAs.

A sparse Bayesian factor model for the construction of gene co-expression networks from single-cell RNA sequencing count data

BACKGROUND

Gene co-expression networks (GCNs) are powerful tools that enable biologists to examine associations between genes during different biological processes. With the advancement of new technologies, such as single-cell RNA sequencing (scRNA-seq), there is a need for developing novel network methods appropriate for new types of data.

RESULTS

We present a novel sparse Bayesian factor model to explore the network structure associated with genes in scRNA-seq data. Latent factors impact the gene expression values for each cell and provide flexibility to account for common features of scRNA-seq: high proportions of zero values, increased cell-to-cell variability, and overdispersion due to abnormally large expression counts. From our model, we construct a GCN by analyzing the positive and negative associations of the factors that are shared between each pair of genes.

CONCLUSIONS

Simulation studies demonstrate that our methodology has high power in identifying gene-gene associations while maintaining a nominal false discovery rate. In real data analyses, our model identifies more known and predicted protein-protein interactions than other competing network models.

A lncRNA landscape in breast cancer reveals a potential role for AC009283.1 in proliferation and apoptosis in HER2-enriched subtype

Breast cancer is the most commonly diagnosed neoplasm in women worldwide with a well-recognized heterogeneous pathology, classified into four molecular subtypes: Luminal A, Luminal B, HER2-enriched and Basal-like, each one with different biological and clinical characteristics. Long non-coding RNAs (lncRNAs) represent 33% of the human transcriptome and play critical roles in breast carcinogenesis, but most of their functions are still unknown. Therefore, cancer research could benefit from continued exploration into the biology of lncRNAs in this neoplasm. We characterized lncRNA expression portraits in 74 breast tumors belonging to the four molecular subtypes using transcriptome microarrays. To infer the biological role of the deregulated lncRNAs in the molecular subtypes, we performed co-expression analysis of lncRNA-mRNA and gene ontology analysis. We identified 307 deregulated lncRNAs in tumor compared to normal tissue and 354 deregulated lncRNAs among the different molecular subtypes. Through co-expression analysis between lncRNAs and protein-coding genes, along with gene enrichment analysis, we inferred the potential function of the most deregulated lncRNAs in each molecular subtype, and independently validated our results taking advantage of TCGA data. Overexpression of the AC009283.1 was observed in the HER2-enriched subtype and it is localized in an amplification zone at chromosome 17q12, suggesting it to be a potential tumorigenic lncRNA. The functional role of lncRNA AC009283.1 was examined through loss of function assays in vitro and determining its impact on global gene expression. These studies revealed that AC009283.1 regulates genes involved in proliferation, cell cycle and apoptosis in a HER2 cellular model. We further confirmed these findings through ssGSEA and CEMITool analysis in an independent HER2-amplified breast cancer cohort. Our findings suggest a wide range of biological functions for lncRNAs in each breast cancer molecular subtype and provide a basis for their biological and functional study, as was conducted for AC009283.1, showing it to be a potential regulator of proliferation and apoptosis in the HER2-enriched subtype.

webCEMiTool: Co-expression Modular Analysis Made Easy

Co-expression analysis has been widely used to elucidate the functional architecture of genes under different biological processes. Such analysis, however, requires substantial knowledge about programming languages and/or bioinformatics skills. We present webCEMiTool,¹ a unique online tool that performs comprehensive modular analyses in a fully automated manner. The webCEMiTool not only identifies co-expression gene modules but also performs several functional analyses on them. In addition, webCEMiTool integrates transcriptomic data with interactome information (i.e., protein-protein interactions) and identifies potential hubs on each network. The tool generates user-friendly html reports that allow users to search for specific genes in each module, as well as check if a module contains genes overrepresented in specific pathways or altered in a specific sample phenotype. We used webCEMiTool to perform a modular analysis of single-cell RNA-seq data of human cells infected with either Zika virus or dengue virus.

scdNet: a computational tool for single-cell differential network analysis

Background

Single-cell RNA sequencing (scRNA-Seq) is an emerging technology that has revolutionized the research of the tumor heterogeneity. However, the highly sparse data matrices generated by the technology have posed an obstacle to the analysis of differential gene regulatory networks.

Results

Addressing the challenges, this study presents, as far as we know, the first bioinformatics tool for scRNA-Seq-based differential network analysis (scdNet). The tool features a sample size adjustment of gene-gene correlation, comparison of inter-state correlations, and construction of differential networks. A simulation analysis demonstrated the power of scdNet in the analyses of sparse scRNA-Seq data matrices, with low requirement on the sample size, high computation efficiency, and tolerance of sequencing noises. Applying the tool to analyze two datasets of single circulating tumor cells (CTCs) of prostate cancer and early mouse embryos, our data demonstrated that differential gene regulation plays crucial roles in anti-androgen resistance and early embryonic development.

Conclusions

Overall, the tool is widely applicable to datasets generated by the emerging technology to bring biological insights into tumor heterogeneity and other studies. MATLAB implementation of scdNet is available at https://github.com/ChenLabGCCRI/scdNet.

Brain Cell Type Specific Gene Expression and Co-expression Network Architectures

Elucidating brain cell type specific gene expression patterns is critical towards a better understanding of how cell-cell communications may influence brain functions and dysfunctions. We set out to compare and contrast five human and murine cell type-specific transcriptome-wide RNA expression data sets that were generated within the past several years. We defined three measures of brain cell type-relative expression including specificity, enrichment, and absolute expression and identified corresponding consensus brain cell “signatures,” which were well conserved across data sets. We validated that the relative expression of top cell type markers are associated with proxies for cell type proportions in bulk RNA expression data from postmortem human brain samples. We further validated novel marker genes using an orthogonal ATAC-seq dataset. We performed multiscale coexpression network analysis of the single cell data sets and identified robust cell-specific gene modules. To facilitate the use of the cell type-specific genes for cell type proportion estimation and deconvolution from bulk brain gene expression data, we developed an R package, BRETIGEA. In summary, we identified a set of novel brain cell consensus signatures and robust networks from the integration of multiple datasets and therefore transcend limitations related to technical issues characteristic of each individual study.