Meta-analysis of cell- specific transcriptomic data using fuzzy c-means clustering discovers versatile viral responsive genes

Type 2 Diabetes Induced by Changes in Proteomic Profiling of Zebrafish Chronically Exposed to a Mixture of Organochlorine Pesticides at Low Concentrations

Effect of organochlorine pesticides (OCPs) mixtures on development of type 2 diabetes mellitus (T2DM) and the underlying mechanism, especially at protein levels, are largely unknown. We exposed a mixture of five OCPs to zebrafish at concentrations of 0, 0.05, 0.25, 2.5, and 25 μg/L for 12 weeks. Differentially expressed proteins (DEPs) were quantitatively identified in female zebrafish livers, and its functional study was conducted. The significantly high glucose and low insulin levels were observed only at 0.05 μg/L, linking to the different pattern of DEPs than other concentrations. A total of 1082 proteins was quantified, of which 321 proteins formed 6 clusters in protein dynamics analysis. The enriched pathways in cluster 3 showing distinct pattern of DEPs could explain the nonlinear response at 0.05 μg/L, indicating that OCP mixtures adversely affected proteins associated with mitochondrial function and energy metabolism. We proposed a feasible mechanism that decrease in expression of aldehyde dehydrogenase led to abnormal accumulation of aldehydes, reducing expression of glyceraldehyde 3-phosphate dehydrogenase, and resulting in disruption of glucose homeostasis. Our findings help to better understand the causality of T2DM by exposure to OCP mixtures and to identify biomarkers in the protein expression level.

Optimizing weighted gene co-expression network analysis with a multi-threaded calculation of the topological overlap matrix

Biomolecular networks are often assumed to be scale-free hierarchical networks. The weighted gene co-expression network analysis (WGCNA) treats gene co-expression networks as undirected scale-free hierarchical weighted networks. The WGCNA R software package uses an Adjacency Matrix to store a network, next calculates the topological overlap matrix (TOM), and then identifies the modules (sub-networks), where each module is assumed to be associated with a certain biological function. The most time-consuming step of WGCNA is to calculate TOM from the Adjacency Matrix in a single thread. In this paper, the single-threaded algorithm of the TOM has been changed into a multi-threaded algorithm (the parameters are the default values of WGCNA). In the multi-threaded algorithm, Rcpp was used to make R call a C++ function, and then C++ used OpenMP to start multiple threads to calculate TOM from the Adjacency Matrix. On shared-memory MultiProcessor systems, the calculation time decreases as the number of CPU cores increases. The algorithm of this paper can promote the application of WGCNA on large data sets, and help other research fields to identify sub-networks in undirected scale-free hierarchical weighted networks. The source codes and usage are available at https://github.com/do-somethings-haha/multi-threaded_calculate_unsigned_TOM_from_unsigned_or_signed_Adjacency_Matrix_of_WGCNA.

Unbiased analysis of peripheral blood mononuclear cells reveals CD4 T cell response to RSV matrix protein

Respiratory syncytial virus (RSV) is the most important cause of respiratory tract illness especially in young infants that develop severe disease requiring hospitalization, and accounting for 74,000–126,000 admissions in the United States (Rezaee et al., 2017; Resch, 2017). Observations of neonatal and infant T cells suggest that they may express different immune markers compared to T-cells from older children. Flow cytometry analysis of cellular responses using “conventional” anti-viral markers (IL2, IFN-γ, TNF, IL10 and IL4) upon RSV-peptide stimulation detected an overall low RSV response in peripheral blood. Therefore we sought an unbiased approach to identify RSV-specific immune markers using RNA-sequencing upon stimulation of infant PBMCs with overlapping peptides representing RSV antigens. To understand the cellular response using transcriptional signatures, transcription factors and cell-type specific signatures were used to investigate breadth of response across peptides. Unexpected from the ICS data, M peptide induced a response equivalent to the F-peptide and was characterized by activation of GATA2, 3, STAT3 and IRF1. This along with upregulation of several unconventional T cell signatures was only observed upon M-peptide stimulation. Moreover, signatures of natural RSV infections were identified from the data available in the public domain to investigate similarities between transcriptional signatures from PBMCs and upon peptide stimulation. This analysis also suggested activation of T cell response upon M-peptide stimulation. Hence, based on transcriptional response, markers were chosen to validate the role of M-peptide in activation of T cells. Indeed, CD4+CXCL9+ cells were identified upon M-peptide stimulation by flow cytometry. Future work using additional markers identified in this study could reveal additional unconventional T cells responding to RSV infections in infants. In conclusion, T cell responses to RSV in infants may not follow the canonical Th1/Th2 patterns of effector responses but include additional functions that may be unique to the neonatal period and correlate with clinical outcomes.

Network Analysis of Large-Scale Data and Its Application to Immunology

Diseases and infections elicit a multilayered immune response which consists of molecular and cellular interaction cascades. Recent advances in high-throughput technologies have facilitated multiparameter investigation of immune cells involved in human immune responses. These multiparameter investigations generate large-scale datasets and advanced computational techniques are required to gain useful information from them. Networks or graphs offer a practical way to represent complex information and develop advanced algorithms to unveil the underlying mechanisms. Here we discuss ways to assemble and analyze networks using genome-wide transcriptional profiles. Additionally, we discuss ways to integrate information available in primary literature and databases with the networks assembled using large-scale datasets. Finally, we describe ways in which network analysis offers insights into human immune responses.

Multi-Objective Optimized Fuzzy Clustering for Detecting Cell Clusters from Single-Cell Expression Profiles

Rapid advance in single-cell RNA sequencing (scRNA-seq) allows measurement of the expression of genes at single-cell resolution in complex disease or tissue. While many methods have been developed to detect cell clusters from the scRNA-seq data, this task currently remains a main challenge. We proposed a multi-objective optimization-based fuzzy clustering approach for detecting cell clusters from scRNA-seq data. First, we conducted initial filtering and SCnorm normalization. We considered various case studies by selecting different cluster numbers ( c l = 2 to a user-defined number), and applied fuzzy c-means clustering algorithm individually. From each case, we evaluated the scores of four cluster validity index measures, Partition Entropy ( P E ), Partition Coefficient ( P C ), Modified Partition Coefficient ( M P C ), and Fuzzy Silhouette Index ( F S I ). Next, we set the first measure as minimization objective (↓) and the remaining three as maximization objectives (↑), and then applied a multi-objective decision-making technique, TOPSIS, to identify the best optimal solution. The best optimal solution (case study) that had the highest TOPSIS score was selected as the final optimal clustering. Finally, we obtained differentially expressed genes (DEGs) using Limma through the comparison of expression of the samples between each resultant cluster and the remaining clusters. We applied our approach to a scRNA-seq dataset for the rare intestinal cell type in mice [GEO ID: GSE62270, 23,630 features (genes) and 288 cells]. The optimal cluster result (TOPSIS optimal score= 0.858) comprised two clusters, one with 115 cells and the other 91 cells. The evaluated scores of the four cluster validity indices, F S I , P E , P C , and M P C for the optimized fuzzy clustering were 0.482, 0.578, 0.607, and 0.215, respectively. The Limma analysis identified 1240 DEGs (cluster 1 vs. cluster 2). The top ten gene markers were Rps21, Slc5a1, Crip1, Rpl15, Rpl3, Rpl27a, Khk, Rps3a1, Aldob and Rps17. In this list, Khk (encoding ketohexokinase) is a novel marker for the rare intestinal cell type. In summary, this method is useful to detect cell clusters from scRNA-seq data.