The mutual information: detecting and evaluating dependencies between variables

Unraveling the evolutionary patterns of construction accidents: a risk assessment framework based on average mutual information theory

Studies of the evolutionary law of construction accidents and the formulation of effective risk assessments are crucial for ensuring construction safety and reducing accidents. However, existing research often focuses exclusively on one of the four critical perspectives: human, facilities, environment, and management, with limited systematic theoretical analysis. Current risk assessments typically emphasized the impact of risks at specific points in time, neglecting the cumulative effects of risks over extended periods. This study has explored how hazardous and harmful factors affect the evolution of accidents, then has established a risk assessment model utilizing average mutual information theory. The findings indicate that "Command violations", "Command errors", and "Illegal operations" are the primary direct causes, while "Inadequate occupational safety and health management structure and staffing" and "Inadequate or unimplemented occupational health management system" are key management factors. In addition, it also provides practical guidance for construction enterprises to strengthen safety management from aspects such as safety production investment and safety management system integration. These contributions are expected to significantly improve the safety level of construction projects and reduce the occurrence of accidents.

Joint analysis of chromatin accessibility and gene expression in the same single cells reveals cancer-specific regulatory programs

Biological analyses conducted at the single-cell scale have revealed profound impacts of heterogeneity and plasticity of chromatin states and gene expression on physiology and cancer. Here, we developed Parallel-seq, a technology for simultaneously measuring chromatin accessibility and gene expression in the same single cells. By combining combinatorial cell indexing and droplet overloading, Parallel-seq generates high-quality data in an ultra-high-throughput fashion and at a cost two orders of magnitude lower than alternative technologies (10× Multiome and ISSAAC-seq). We applied Parallel-seq to 40 lung tumor and tumor-adjacent clinical samples and obtained over 200,000 high-quality joint scATAC-and-scRNA profiles. Leveraging this large dataset, we characterized copy-number variations (CNVs) and extrachromosomal circular DNA (eccDNA) heterogeneity in tumor cells, predicted hundreds of thousands of cell-type-specific regulatory events, and identified enhancer mutations affecting tumor progression. Our analyses highlight Parallel-seq's power in investigating epigenetic and genetic factors driving cancer development at the cell-type-specific level and its utility for revealing vulnerable therapeutic targets.

ICRL: independent causality representation learning for domain generalization

Domain generalization (DG) addresses the challenge of out-of-distribution (OOD) data; however, the reliance on statistical correlations during model development often introduces shortcut learning problems. Current approaches to mitigating these issues commonly involve the integration of causal inference, formalizing DG problems through a general structural causal model. Nevertheless, ensuring the independence of features when incorporating causal models is often overlooked, leading to spurious causal relationships. In this work, we design three independent feature modules using GAN variants (GAN, WGAN, and WGAN-GP) and select the best-performing WGAN module to integrate into the existing causal model framework, constructing an independent causal relationship learning (ICRL) model. Extensive experiments on widely used datasets demonstrate that our model, with independent causal representations, outperforms the original model in both performance and efficiency, thereby validating the effectiveness of our proposed approach. The code for ICRL can be accessed at: https://github.com/22Shao/ICRL.git .

Differential Transcriptional Programs Reveal Modular Network Rearrangements Associated with Late-Onset Alzheimer's Disease

Alzheimer's disease (AD) is a complex, genetically heterogeneous disorder. The diverse phenotypes associated with AD result from interactions between genetic and environmental factors, influencing multiple biological pathways throughout disease progression. Network-based approaches offer a way to assess phenotype-specific states. In this study, we calculated key network metrics to characterize the network transcriptional structure and organization in LOAD, focusing on genes and pathways implicated in AD pathology within the dorsolateral prefrontal cortex (DLPFC). Our findings revealed disease-specific coexpression markers associated with diverse metabolic functions. Additionally, significant differences were observed at both the mesoscopic and local levels between AD and control networks, along with a restructuring of gene coexpression and biological functions into distinct transcriptional modules. These results show the molecular reorganization of the transcriptional program occurring in LOAD, highlighting specific adaptations that may contribute to or result from cellular responses to pathological stressors. Our findings may support the development of a unified model for the causal mechanisms of AD, suggesting that its diverse manifestations arise from multiple pathways working together to produce the disease's complex clinical patho-phenotype.

Optimizing functional brain network analysis by incorporating nonlinear factors and frequency band selection with machine learning models

The accurate assessment of the brain's functional network is seen as crucial for the understanding of complex relationships between different brain regions. Hidden information within different frequency bands, which is often overlooked by traditional linear correlation-based methods such as Pearson correlation (PC) and partial correlation, fails to be revealed, leading to the neglect of more intricate nonlinear factors. These limitations were aimed to be overcome in this study by the combination of fast continuous wavelet transform and normalized mutual information (NMI) to develop a novel approach. Original time-domain signals from resting-state functional magnetic resonance imaging were decomposed into different frequency domains using fast continuous wavelet transform, and adjacency matrices were constructed to enhance feature separation across brain regions. Both linear and nonlinear aspects between brain regions were comprehensively considered through the integration of complex correlation coefficient and NMI. The construction of functional brain networks was enabled by the adaptive selection of optimal frequency band combinations. The construction of the model was facilitated by feature extraction using tree models with extreme gradient boosting. It was demonstrated through comparative analysis that the method outperformed baseline methods such as PC and NMI, achieving an area under the curve of 0.9054. The introduction of nonlinear factors was found to increase precision by 14.25% and recall by 17.14%. Importantly, the approach optimized the original data without significantly altering the feature topology. Overall, this innovation advances the understanding of brain function, offering more accurate potential for future research and clinical applications.

ISCAZIM: Integrated statistical correlation analysis for zero-inflated microbiome data

Microbiome-metabolome association analysis is critical to reveal the key pairs of gut microbiota and metabolites for discovery of the microbial biomarkers in chronic diseases. However, the characteristics of microbiome data, such as zero inflation, over dispersion, may impair the confidence of association analysis between microbiome and metabolome data. The objectives of this study are to evaluate the strengths and weaknesses of existing statistical methods and to develop a computational framework tailored to the unique characteristics of microbiome data. We designed a computational framework called Integrated Statistical Correlation Analysis for Zero-Inflated Microbiome data (ISCAZIM) that takes account of complicated microbiome data characteristics, including zero inflation rates (ZIRs), dispersion and correlation patterns. ISCAZIM first benchmarked prevalent statistical correlation methods, Pearson, Spearman, zero inflated negative binomial (ZINB) model, mutual information and Maximal Information Coefficient. ISCAZIM then classifies the correlation pattern to linear or non-linear and applies the correlation method according to the ZIRs status. Applying to multiple real-world microbiome-metabolomics data, ISCAZIM is overall more accurate than using a single method with more truly significant association pairs included. Therefore, ISCAZIM will significantly facilitate the association analysis using zero-inflated microbiome data for multi-omics integration.

Conformational dynamics and multi-modal interaction of Paxillin with the Focal Adhesion Targeting Domain

Paxillin (PXN) and focal adhesion kinase (FAK) are two major components of the focal adhesion complex, a multiprotein structure linking the intracellular cytoskeleton to the cell exterior. PXN interacts directly with the C-terminal targeting domain of FAK (FAT) via its intrinsically disordered N-terminal domain. This interaction is necessary and sufficient for localizing FAK to focal adhesions. Furthermore, PXN serves as a platform for recruiting other proteins that together control the dynamic changes needed for cell migration and survival. Here, we show that the PXN disordered region undergoes large-scale conformational restriction upon binding to FAT, forming a 48-kDa multi-modal complex consisting of four major interconverting states. Although the complex is flexible, each state has unique sets of contacts involving disordered regions that are both highly represented in ensembles and conserved. Moreover, conserved intramolecular contacts from glutamine-rich regions in PXN contribute to high entropy and thus stability of the FAT bound complex. As PXN is a hub protein, the results provide a structural basis for understanding how perturbations that lead to cellular network rewiring, such as ligand binding and phosphorylation, may lead to shifts in the multi-state equilibrium and phenotypic switching.

PyNetCor: a high-performance Python package for large-scale correlation analysis

The development of multi-omics technologies has generated an abundance of biological datasets, providing valuable resources for investigating potential relationships within complex biological systems. However, most correlation analysis tools face computational challenges when dealing with these high-dimensional datasets containing millions of features. Here, we introduce pyNetCor, a fast and scalable tool for constructing correlation networks on large-scale and high-dimensional data. PyNetCor features optimized algorithms for both full correlation coefficient matrix computation and top-k correlation search, outperforming other tools in the field in terms of runtime and memory consumption. It utilizes a linear interpolation strategy to rapidly estimate P-values and achieve false discovery rate control, demonstrating a speedup of over 110 times compared to existing methods. Overall, pyNetCor supports large-scale correlation analysis, a crucial foundational step for various bioinformatics workflows, and can be easily integrated into downstream applications to accelerate the process of extracting biological insights from data.

Multiscale detrended cross-correlation coefficient: estimating coupling in non-stationary neurophysiological signals

The brain consists of a vastly interconnected network of regions, the connectome. By estimating the statistical interdependence of neurophysiological time series, we can measure the functional connectivity (FC) of this connectome. Pearson's correlation (r P) is a common metric of coupling in FC studies. Yet r P does not account properly for the non-stationarity of the signals recorded in neuroimaging. In this study, we introduced a novel estimator of coupled dynamics termed multiscale detrended cross-correlation coefficient (MDC3). Firstly, we showed that MDC3 had higher accuracy compared to r P and lagged covariance using simulated time series with known coupling, as well as simulated functional magnetic resonance imaging (fMRI) signals with known underlying structural connectivity. Next, we computed functional brain networks based on empirical magnetoencephalography (MEG) and fMRI. We found that by using MDC3 we could construct networks of healthy populations with significantly different properties compared to r P networks. Based on our results, we believe that MDC3 is a valid alternative to r P that should be incorporated in future FC studies.

A cross-omics data analysis strategy for metabolite-microbe pair identification

Given the pivotal roles of metabolomics and microbiomics, numerous data mining approaches aim to uncover their intricate connections. However, the complex many-to-many associations between metabolome-microbiome profiles yield numerous statistically significant but biologically unvalidated candidates. To address these challenges, we introduce BiOFI, a strategic framework for identifying metabolome-microbiome correlation pairs (Bi-Omics). BiOFI employs a comprehensive scoring system, incorporating intergroup differences, effects on feature correlation networks, and organism abundance. Meanwhile, it establishes a built-in database of metabolite-microbe-KEGG functional pathway linking relationships. Furthermore, BiOFI can rank related feature pairs by combining importance scores and correlation strength. Validation on a dataset of cesarean-section infants confirms the strategy's validity and interpretability. The BiOFI R package is freely accessible at https://github.com/chentianlu/BiOFI.

NetSci: A Library for High Performance Biomolecular Simulation Network Analysis Computation

We present the NetSci program-an open-source scientific software package designed for estimating mutual information (MI) between data sets using GPU acceleration and a k-nearest-neighbor algorithm. This approach significantly enhances calculation speed, achieving improvements of several orders of magnitude over traditional CPU-based methods, with data set size limits dictated only by available hardware. To validate NetSci, we accurately compute MI for an analytically verifiable two-dimensional Gaussian distribution and replicate the generalized correlation (GC) analysis previously conducted on the B1 domain of protein G. We also apply NetSci to molecular dynamics simulations of the Sarcoendoplasmic Reticulum Calcium-ATPase (SERCA) pump, exploring the allosteric mechanisms and pathways influenced by ATP and 2'-deoxy-ATP (dATP) binding. Our analysis reveals distinct allosteric effects induced by ATP compared to dATP, with predicted information pathways from the bound nucleotide to the calcium-binding domain differing based on the nucleotide involved. NetSci proves to be a valuable tool for estimating MI and GC in various data sets and is particularly effective for analyzing intraprotein communication and information transfer.

Relevance of the regulation of the brain-placental axis to the nocturnal bottleneck of mammals

INTRODUCTION

Evolutionary theory suggests that the ancestors of all placental animals were nocturnal. Visual perceptive function of mammalian brain has evolved extensively, but nearly 70 % of today's mammals are still nocturnal. While placental influence on brain development is known, if placenta plays a role in the visual perceptive function of mammalian brain remains untested. The present study aims to test this hypothesis.

METHODS

In this study, single-nuclei RNA sequencing was performed to identify genes expressed in the pig placenta and fetal brain, and then compared with the orthologous genes expressed in the placenta and fetal brain cells of mouse. Differential gene expression analysis was performed to identify placental genes regulated differentially between nocturnal and diurnal animals. Phylogenetic modeling was performed to test correlated evolution between placenta type, and the nocturnal or diurnal activity among different mammals.

RESULTS

The results showed that genes differentially regulated in the fetal brain were related to visual perception whereas the placental genes were related to the nocturnal or diurnal activity in placental animals. Phylogenetic modeling of these genes in thirty-four diverse mammalian species showed evidence for evolutionary link between placenta and the nocturnal/diurnal activity in animals.

DISCUSSION

The findings of this study suggest that the placenta plays a role in the evolution of visual perceptive function of brain to shape the nocturnal or diurnal activity of placental animals.

Inference of single-cell network using mutual information for scRNA-seq data analysis

BACKGROUND

With the advance in single-cell RNA sequencing (scRNA-seq) technology, deriving inherent biological system information from expression profiles at a single-cell resolution has become possible. It has been known that network modeling by estimating the associations between genes could better reveal dynamic changes in biological systems. However, accurately constructing a single-cell network (SCN) to capture the network architecture of each cell and further explore cell-to-cell heterogeneity remains challenging.

RESULTS

We introduce SINUM, a method for constructing the SIngle-cell Network Using Mutual information, which estimates mutual information between any two genes from scRNA-seq data to determine whether they are dependent or independent in a specific cell. Experiments on various scRNA-seq datasets with different cell numbers based on eight performance indexes (e.g., adjusted rand index and F-measure index) validated the accuracy and robustness of SINUM in cell type identification, superior to the state-of-the-art SCN inference method. Additionally, the SINUM SCNs exhibit high overlap with the human interactome and possess the scale-free property.

CONCLUSIONS

SINUM presents a view of biological systems at the network level to detect cell-type marker genes/gene pairs and investigate time-dependent changes in gene associations during embryo development. Codes for SINUM are freely available at https://github.com/SysMednet/SINUM .

Beyond Pearson's Correlation: Modern Nonparametric Independence Tests for Psychological Research

When examining whether two continuous variables are associated, tests based on Pearson's, Kendall's, and Spearman's correlation coefficients are typically used. This paper explores modern nonparametric independence tests as an alternative, which, unlike traditional tests, have the ability to potentially detect any type of relationship. In addition to existing modern nonparametric independence tests, we developed and considered two novel variants of existing tests, most notably the Heller-Heller-Gorfine-Pearson (HHG-Pearson) test. We conducted a simulation study to compare traditional independence tests, such as Pearson's correlation, and the modern nonparametric independence tests in situations commonly encountered in psychological research. As expected, no test had the highest power across all relationships. However, the distance correlation and the HHG-Pearson tests were found to have substantially greater power than all traditional tests for many relationships and only slightly less power in the worst case. A similar pattern was found in favor of the HHG-Pearson test compared to the distance correlation test. However, given that distance correlation performed better for linear relationships and is more widely accepted, we suggest considering its use in place or additional to traditional methods when there is no prior knowledge of the relationship type, as is often the case in psychological research.

Microbial reaction rate estimation using proteins and proteomes

Microbes transform their environments using diverse enzymatic reactions. However, it remains challenging to measure microbial reaction rates in natural environments. Despite advances in global quantification of enzyme abundances, the individual relationships between enzyme abundances and their reaction rates have not been systematically examined. Using matched proteomic and reaction rate data from microbial cultures, we show that enzyme abundance is often insufficient to predict its corresponding reaction rate. However, we discovered that global proteomic measurements can be used to make accurate rate predictions of individual reaction rates (median R 2 = 0.78). Accurate rate predictions required only a small number of proteins and they did not need explicit prior mechanistic knowledge or environmental context. These results indicate that proteomes are encoders of cellular reaction rates, potentially enabling proteomic measurements in situ to estimate the rates of microbially mediated reactions in natural systems.

A comparative analysis of mutual information methods for pairwise relationship detection in metagenomic data

BACKGROUND

Construction of co-occurrence networks in metagenomic data often employs correlation to infer pairwise relationships between microbes. However, biological systems are complex and often display qualities non-linear in nature. Therefore, the reliance on correlation alone may overlook important relationships and fail to capture the full breadth of intricacies presented in underlying interaction networks. It is of interest to incorporate metrics that are not only robust in detecting linear relationships, but non-linear ones as well.

RESULTS

In this paper, we explore the use of various mutual information (MI) estimation approaches for quantifying pairwise relationships in biological data and compare their performances against two traditional measures-Pearson's correlation coefficient, r, and Spearman's rank correlation coefficient, ρ. Metrics are tested on both simulated data designed to mimic pairwise relationships that may be found in ecological systems and real data from a previous study on C. diff infection. The results demonstrate that, in the case of asymmetric relationships, mutual information estimators can provide better detection ability than Pearson's or Spearman's correlation coefficients. Specifically, we find that these estimators have elevated performances in the detection of exploitative relationships, demonstrating the potential benefit of including them in future metagenomic studies.

CONCLUSIONS

Mutual information (MI) can uncover complex pairwise relationships in biological data that may be missed by traditional measures of association. The inclusion of such relationships when constructing co-occurrence networks can result in a more comprehensive analysis than the use of correlation alone.

Assessment of Fractal Synchronization during an Epileptic Seizure

In this paper, we define fractal synchronization (FS) based on the idea of stochastic synchronization and propose a mathematical apparatus for estimating FS. One major advantage of our proposed approach is that fractal synchronization makes it possible to estimate the aggregate strength of the connection on multiple time scales between two projections of the attractor, which are time series with a fractal structure. We believe that one of the promising uses of FS is the assessment of the interdependence of encephalograms. To demonstrate this approach in evaluating the cross-dependence between channels in a network of electroencephalograms, we evaluated the FS of encephalograms during an epileptic seizure. Fractal synchronization demonstrates the presence of desynchronization during an epileptic seizure.

Targeted co-expression networks for the study of traits

Weighted Gene Co-expression Network Analysis (WGCNA) is a widely used approach for the generation of gene co-expression networks. However, networks generated with this tool usually create large modules with a large set of functional annotations hard to decipher. We have developed TGCN, a new method to create Targeted Gene Co-expression Networks. This method identifies the transcripts that best predict the trait of interest based on gene expression using a refinement of the LASSO regression. Then, it builds the co-expression modules around those transcripts. Algorithm properties were characterized using the expression of 13 brain regions from the Genotype-Tissue Expression project. When comparing our method with WGCNA, TGCN networks lead to more precise modules that have more specific and yet rich biological meaning. Then, we illustrate its applicability by creating an APP-TGCN on The Religious Orders Study and Memory and Aging Project dataset, aiming to identify the molecular pathways specifically associated with APP role in Alzheimer's disease. Main biological findings were further validated in two independent cohorts. In conclusion, we provide a new framework that serves to create targeted networks that are smaller, biologically relevant and useful in high throughput hypothesis driven research. The TGCN R package is available on Github: https://github.com/aliciagp/TGCN .

CoVar: A generalizable machine learning approach to identify the coordinated regulators driving variational gene expression

Network inference is used to model transcriptional, signaling, and metabolic interactions among genes, proteins, and metabolites that identify biological pathways influencing disease pathogenesis. Advances in machine learning (ML)-based inference models exhibit the predictive capabilities of capturing latent patterns in genomic data. Such models are emerging as an alternative to the statistical models identifying causative factors driving complex diseases. We present CoVar, an ML-based framework that builds upon the properties of existing inference models, to find the central genes driving perturbed gene expression across biological states. Unlike differentially expressed genes (DEGs) that capture changes in individual gene expression across conditions, CoVar focuses on identifying variational genes that undergo changes in their expression network interaction profiles, providing insights into changes in the regulatory dynamics, such as in disease pathogenesis. Subsequently, it finds core genes from among the nearest neighbors of these variational genes, which are central to the variational activity and influence the coordinated regulatory processes underlying the observed changes in gene expression. Through the analysis of simulated as well as yeast expression data perturbed by the deletion of the mitochondrial genome, we show that CoVar captures the intrinsic variationality and modularity in the expression data, identifying key driver genes not found through existing differential analysis methodologies.

From local to global gene co-expression estimation using single-cell RNA-seq data

In genomics studies, the investigation of gene relationships often brings important biological insights. Currently, the large heterogeneous datasets impose new challenges for statisticians because gene relationships are often local. They change from one sample point to another, may only exist in a subset of the sample, and can be nonlinear or even nonmonotone. Most previous dependence measures do not specifically target local dependence relationships, and the ones that do are computationally costly. In this paper, we explore a state-of-the-art network estimation technique that characterizes gene relationships at the single cell level, under the name of cell-specific gene networks. We first show that averaging the cell-specific gene relationship over a population gives a novel univariate dependence measure, the averaged Local Density Gap (aLDG), that accumulates local dependence and can detect any nonlinear, nonmonotone relationship. Together with a consistent nonparametric estimator, we establish its robustness on both the population and empirical levels. Then, we show that averaging the cell-specific gene relationship over mini-batches determined by some external structure information (eg, spatial or temporal factor) better highlights meaningful local structure change points. We explore the application of aLDG and its minibatch variant in many scenarios, including pairwise gene relationship estimation, bifurcating point detection in cell trajectory, and spatial transcriptomics structure visualization. Both simulations and real data analysis show that aLDG outperforms existing ones.

Inference of Dynamic Growth Regulatory Network in Cancer Using High-Throughput Transcriptomic Data

Growth is regulated by gene expression variation at different developmental stages of biological processes such as cell differentiation, disease progression, or drug response. In cancer, a stage-specific regulatory model constructed to infer the dynamic expression changes in genes contributing to tissue growth or proliferation is referred as a dynamic growth regulatory network (dGRN). Over the past decade, gene expression data has been widely used for reconstructing dGRN by computing correlations between the differentially expressed genes (DEGs). A wide variety of pipelines are available to construct the GRNs using DEGs and the choice of a particular method or tool depends on the nature of the study. In this protocol, we have outlined a step-by-step guide for the analysis of DEGs using RNA-Seq data, beginning from data acquisition, pre-processing, mapping to reference genome, and construction of a correlation-based co-expression network to further downstream analysis. We have also outlined the steps for the inclusion of publicly available interaction/regulation information into the dGRN followed by relevant topological inferences. This tutorial has been designed in a way that early researchers can refer to for an easy and comprehensive glimpse of methodologies used in the inference of dGRN using transcriptomics data.

Does DNA methylation in the fetal brain leave an epigenetic memory in the blood?

Epigenetic memory is an emerging concept that refers to the process in which epigenetic changes occurring early-in life can lead to long-term programs of gene regulation in time and space. By leveraging neural network regression modeling of DNA methylation data in pigs, we show that specific methylations in the adult blood can reliably predict methylation changes that occurred in the fetal brain. Genes associated with these methylations represented known markers of specific cell types of blood including bone marrow hematopoietic progenitor cells, and ependymal and oligodendrocyte cells of brain. This suggested that methylation changes that occurred in the developing brain were maintained as an epigenetic memory in the blood through the adult life.

The non-specific lethal complex regulates genes and pathways genetically linked to Parkinson's disease

Genetic variants conferring risks for Parkinson's disease have been highlighted through genome-wide association studies, yet exploration of their specific disease mechanisms is lacking. Two Parkinson's disease candidate genes, KAT8 and KANSL1, identified through genome-wide studies and a PINK1-mitophagy screen, encode part of the histone acetylating non-specific lethal complex. This complex localizes to the nucleus, where it plays a role in transcriptional activation, and to mitochondria, where it has been suggested to have a role in mitochondrial transcription. In this study, we sought to identify whether the non-specific lethal complex has potential regulatory relationships with other genes associated with Parkinson's disease in human brain. Correlation in the expression of non-specific lethal genes and Parkinson's disease-associated genes was investigated in primary gene co-expression networks using publicly-available transcriptomic data from multiple brain regions (provided by the Genotype-Tissue Expression Consortium and UK Brain Expression Consortium), whilst secondary networks were used to examine cell type specificity. Reverse engineering of gene regulatory networks generated regulons of the complex, which were tested for heritability using stratified linkage disequilibrium score regression. Prioritized gene targets were then validated in vitro using a QuantiGene multiplex assay and publicly-available chromatin immunoprecipitation-sequencing data. Significant clustering of non-specific lethal genes was revealed alongside Parkinson's disease-associated genes in frontal cortex primary co-expression modules, amongst other brain regions. Both primary and secondary co-expression modules containing these genes were enriched for mainly neuronal cell types. Regulons of the complex contained Parkinson's disease-associated genes and were enriched for biological pathways genetically linked to disease. When examined in a neuroblastoma cell line, 41% of prioritized gene targets showed significant changes in mRNA expression following KANSL1 or KAT8 perturbation. KANSL1 and H4K8 chromatin immunoprecipitation-sequencing data demonstrated non-specific lethal complex activity at many of these genes. In conclusion, genes encoding the non-specific lethal complex are highly correlated with and regulate genes associated with Parkinson's disease. Overall, these findings reveal a potentially wider role for this protein complex in regulating genes and pathways implicated in Parkinson's disease.

Environmental factors influencing DDT-DDE spatial distribution in an agricultural drainage system determined by using machine learning techniques

The presence and persistence of pesticides in the environment are environmental problems of great concern due to the health implications for humans and wildlife. The persistence of DDT-DDE in a Mediterranean coastal plain where pesticides were widely used and were banned decades ago is the aim of this study. Different sources of analytical information from water and soil analysis and topography and geographical variables were combined with the purpose of analyzing which environmental factors are more likely to condition the spatial distribution of DDT-DDE in the drainage watercourses of the area. An approach combining machine learning techniques, such as Random Forest and Mutual Information (MI), for classifying DDT-DDE concentration levels based on other environmental predictive variables was applied. In addition, classification procedure was iteratively performed with different training/validation partitions in order to extract the most informative parameters denoted by the highest MI scores and larger accuracy assessment metrics. Distance to drain canals, soil electrical conductivity, and soil sand texture fraction were the most informative environmental variables for predicting DDT-DDE water concentration clusters.

Role of caveolin-1 in metabolic programming of fetal brain

Mice lacking caveolin-1 (Cav1), a key protein of plasma membrane, exhibit brain aging at an early adult stage. Here, integrative analyses of metabolomics, transcriptomics, epigenetics, and single-cell data were performed to test the hypothesis that metabolic deregulation of fetal brain due to the ablation of Cav1 is linked to brain aging in these mice. The results of this study show that lack of Cav1 caused deregulation in the lipid and amino acid metabolism in the fetal brain, and genes associated with these deregulated metabolites were significantly altered in the brain upon aging. Moreover, ablation of Cav1 deregulated several metabolic genes in specific cell types of the fetal brain and impacted DNA methylation of those genes in coordination with mouse epigenetic clock. The findings of this study suggest that the aging program of brain is confounded by metabolic abnormalities in the fetal stage due to the absence of Cav1.

Capturing synchronization with complexity measure of ordinal pattern transition network constructed by crossplot

To evaluate the synchronization of bivariate time series has been a hot topic, and a number of measures have been proposed. In this work, by introducing the ordinal pattern transition network into the crossplot, a new method for measuring the synchronization of bivariate time series is proposed. After the crossplot been partitioned and coded, the coded partitions are defined as network nodes and a directed weighted network is constructed based on the temporal adjacency of the nodes. The crossplot transition entropy of the network is proposed as an indicator of the synchronization between two time series. To test the characteristics and performance of the method, it is used to analyse the unidirectional coupled Lorentz model and compared it with existing methods. The results showed the new method had the advantages of easy parameter setting, efficiency, robustness, good consistency and suitability for short time series. Finally, electroencephalogram (EEG) data from auditory-evoked potential EEG-biometric dataset are investigated, and some useful and interesting results are obtained.

Normalized mutual information of fNIRS signals as a measure for accessing typical and atypical brain activity

Normalized mutual information (NMI) can be used to detect statistical correlations between time series. We showed possibility of using NMI to quantify synchronicity of information transmission in different brain regions, thus to characterize functional connections, and ultimately analyze differences in physiological states of brain. Resting-state brain signals were recorded from bilateral temporal lobes by functional near-infrared spectroscopy (fNIRS) in 19 young healthy (YH) adults, 25 children with autism spectrum disorder (ASD), and 22 children with typical development (TD). Using NMI of the fNIRS signals, common information volume was assessed for each of three groups. Results showed that mutual information of children with ASD was significantly smaller than that of TD children, while mutual information of YH adults was slightly larger than that of TD children. This study may suggest that NMI could be a measure for assessing brain activity with different development states.

Computational design of dynamic receptor-peptide signaling complexes applied to chemotaxis

Engineering protein biosensors that sensitively respond to specific biomolecules by triggering precise cellular responses is a major goal of diagnostics and synthetic cell biology. Previous biosensor designs have largely relied on binding structurally well-defined molecules. In contrast, approaches that couple the sensing of flexible compounds to intended cellular responses would greatly expand potential biosensor applications. Here, to address these challenges, we develop a computational strategy for designing signaling complexes between conformationally dynamic proteins and peptides. To demonstrate the power of the approach, we create ultrasensitive chemotactic receptor-peptide pairs capable of eliciting potent signaling responses and strong chemotaxis in primary human T cells. Unlike traditional approaches that engineer static binding complexes, our dynamic structure design strategy optimizes contacts with multiple binding and allosteric sites accessible through dynamic conformational ensembles to achieve strongly enhanced signaling efficacy and potency. Our study suggests that a conformationally adaptable binding interface coupled to a robust allosteric transmission region is a key evolutionary determinant of peptidergic GPCR signaling systems. The approach lays a foundation for designing peptide-sensing receptors and signaling peptide ligands for basic and therapeutic applications.

Gene regulation network inference using k-nearest neighbor-based mutual information estimation: revisiting an old DREAM

BACKGROUND

A cell exhibits a variety of responses to internal and external cues. These responses are possible, in part, due to the presence of an elaborate gene regulatory network (GRN) in every single cell. In the past 20 years, many groups worked on reconstructing the topological structure of GRNs from large-scale gene expression data using a variety of inference algorithms. Insights gained about participating players in GRNs may ultimately lead to therapeutic benefits. Mutual information (MI) is a widely used metric within this inference/reconstruction pipeline as it can detect any correlation (linear and non-linear) between any number of variables (n-dimensions). However, the use of MI with continuous data (for example, normalized fluorescence intensity measurement of gene expression levels) is sensitive to data size, correlation strength and underlying distributions, and often requires laborious and, at times, ad hoc optimization.

RESULTS

In this work, we first show that estimating MI of a bi- and tri-variate Gaussian distribution using k-nearest neighbor (kNN) MI estimation results in significant error reduction as compared to commonly used methods based on fixed binning. Second, we demonstrate that implementing the MI-based kNN Kraskov-Stoögbauer-Grassberger (KSG) algorithm leads to a significant improvement in GRN reconstruction for popular inference algorithms, such as Context Likelihood of Relatedness (CLR). Finally, through extensive in-silico benchmarking we show that a new inference algorithm CMIA (Conditional Mutual Information Augmentation), inspired by CLR, in combination with the KSG-MI estimator, outperforms commonly used methods.

CONCLUSIONS

Using three canonical datasets containing 15 synthetic networks, the newly developed method for GRN reconstruction-which combines CMIA, and the KSG-MI estimator-achieves an improvement of 20-35% in precision-recall measures over the current gold standard in the field. This new method will enable researchers to discover new gene interactions or better choose gene candidates for experimental validations.

Human and macaque pairs employ different coordination strategies in a transparent decision game

Many real-world decisions in social contexts are made while observing a partner's actions. To study dynamic interactions during such decisions, we developed a setup where two agents seated face-to-face to engage in game-theoretical tasks on a shared transparent touchscreen display ('transparent games'). We compared human and macaque pairs in a transparent version of the coordination game 'Bach-or-Stravinsky', which entails a conflict about which of two individually-preferred opposing options to choose to achieve coordination. Most human pairs developed coordinated behavior and adopted dynamic turn-taking to equalize the payoffs. All macaque pairs converged on simpler, static coordination. Remarkably, two animals learned to coordinate dynamically after training with a human confederate. This pair selected the faster agent's preferred option, exhibiting turn-taking behavior that was captured by modeling the visibility of the partner's action before one's own movement. Such competitive turn-taking was unlike the prosocial turn-taking in humans, who equally often initiated switches to and from their preferred option. Thus, the dynamic coordination is not restricted to humans but can occur on the background of different social attitudes and cognitive capacities in rhesus monkeys. Overall, our results illustrate how action visibility promotes the emergence and maintenance of coordination when agents can observe and time their mutual actions.

CoVar: A generalizable machine learning approach to identify the coordinated regulators driving variational gene expression

Network inference is used to model transcriptional, signaling, and metabolic interactions among genes, proteins, and metabolites that identify biological pathways influencing disease pathogenesis. Advances in machine learning (ML)-based inference models exhibit the predictive capabilities of capturing latent patterns in genomic data. Such models are emerging as an alternative to the statistical models identifying causative factors driving complex diseases. We present CoVar, an inference framework that builds upon the properties of existing inference models, to find the central genes driving perturbed gene expression across biological states. We leverage ML-based network inference to find networks that capture the strength of regulatory interactions. Our model first pinpoints a subset of genes, termed variational, whose expression variabilities typify the differences in network connectivity between the control and perturbed data. Variational genes, by being differentially expressed themselves or possessing differentially expressed neighbor genes, capture gene expression variability. CoVar then creates subnetworks comprising variational genes and their strongly connected neighbor genes and identifies core genes central to these subnetworks that influence the bulk of the variational activity. Through the analysis of yeast expression data perturbed by the deletion of the mitochondrial genome, we show that CoVar identifies key genes not found through independent differential expression analysis.

Human activity recognition from sensor data using spatial attention-aided CNN with genetic algorithm

Capturing time and frequency relationships of time series signals offers an inherent barrier for automatic human activity recognition (HAR) from wearable sensor data. Extracting spatiotemporal context from the feature space of the sensor reading sequence is challenging for the current recurrent, convolutional, or hybrid activity recognition models. The overall classification accuracy also gets affected by large size feature maps that these models generate. To this end, in this work, we have put forth a hybrid architecture for wearable sensor data-based HAR. We initially use Continuous Wavelet Transform to encode the time series of sensor data as multi-channel images. Then, we utilize a Spatial Attention-aided Convolutional Neural Network (CNN) to extract higher-dimensional features. To find the most essential features for recognizing human activities, we develop a novel feature selection (FS) method. In order to identify the fitness of the features for the FS, we first employ three filter-based methods: Mutual Information (MI), Relief-F, and minimum redundancy maximum relevance (mRMR). The best set of features is then chosen by removing the lower-ranked features using a modified version of the Genetic Algorithm (GA). The K-Nearest Neighbors (KNN) classifier is then used to categorize human activities. We conduct comprehensive experiments on five well-known, publicly accessible HAR datasets, namely UCI-HAR, WISDM, MHEALTH, PAMAP2, and HHAR. Our model significantly outperforms the state-of-the-art models in terms of classification performance. We also observe an improvement in overall recognition accuracy with the use of GA-based FS technique with a lower number of features. The source code of the paper is publicly available here https://github.com/apusarkar2195/HAR_WaveletTransform_SpatialAttention_FeatureSelection.

Association Measure and Compact Prediction for Chemical Process Data from an Information-Theoretic Perspective

Mutual information (MI) has been widely used for association mining in complex chemical processes, but how to precisely estimate MI between variables of different numerical types, discriminate their association relationships with targets and finally achieve compact and interpretable prediction has not been discussed in detail, which may limit MI in more complicated industrial applications. Therefore, this paper first reviews the existing information-based association measures and proposes a general framework, GIEF, to consistently detect associations and independence between different types of variables. Then, the study defines four mutually exclusive association relations of variables from an information-theoretic perspective to guide feature selection and compact prediction in high-dimensional processes. Based on GIEF and conditional mutual information maximization (CMIM), a new algorithm, CMIM-GIEF, is proposed and tested on a fluidized catalytic cracking (FCC) process with 217 variables, one which achieves significantly improved accuracies with fewer variables in predicting the yields of four crucial products. The compact variables identified are also consistent with the results of Shapley Additive exPlanations (SHAP) and industrial experience, proving good adaptivity of the method for chemical process data.

Bayesian SIR model with change points with application to the Omicron wave in Singapore

The Omicron variant has led to a new wave of the COVID-19 pandemic worldwide, with unprecedented numbers of daily confirmed new cases in many countries and areas. To analyze the impact of society or policy changes on the development of the Omicron wave, the stochastic susceptible-infected-removed (SIR) model with change points is proposed to accommodate the situations where the transmission rate and the removal rate may vary significantly at change points. Bayesian inference based on a Markov chain Monte Carlo algorithm is developed to estimate both the locations of change points as well as the transmission rate and removal rate within each stage. Experiments on simulated data reveal the effectiveness of the proposed method, and several stages are detected in analyzing the Omicron wave data in Singapore.

Vaginal microbiota networks as a mechanistic predictor of aerobic vaginitis

Aerobic vaginitis (AV) is a complex vaginal dysbiosis that is thought to be caused by the micro-ecological change of the vaginal microbiota. While most studies have focused on how changes in the abundance of individual microbes are associated with the emergence of AV, we still do not have a complete mechanistic atlas of the microbe-AV link. Network modeling is central to understanding the structure and function of any microbial community assembly. By encapsulating the abundance of microbes as nodes and ecological interactions among microbes as edges, microbial networks can reveal how each microbe functions and how one microbe cooperate or compete with other microbes to mediate the dynamics of microbial communities. However, existing approaches can only estimate either the strength of microbe-microbe link or the direction of this link, failing to capture full topological characteristics of a network, especially from high-dimensional microbial data. We combine allometry scaling law and evolutionary game theory to derive a functional graph theory that can characterize bidirectional, signed, and weighted interaction networks from any data domain. We apply our theory to characterize the causal interdependence between microbial interactions and AV. From functional networks arising from different functional modules, we find that, as the only favorable genus from Firmicutes among all identified genera, the role of Lactobacillus in maintaining vaginal microbial symbiosis is enabled by upregulation from other microbes, rather than through any intrinsic capacity. Among Lactobacillus species, the proportion of L. crispatus to L. iners is positively associated with more healthy acid vaginal ecosystems. In a less healthy alkaline ecosystem, L. crispatus establishes a contradictory relationship with other microbes, leading to population decrease relative to L. iners. We identify topological changes of vaginal microbiota networks when the menstrual cycle of women changes from the follicular to luteal phases. Our network tool provides a mechanistic approach to disentangle the internal workings of the microbiota assembly and predict its causal relationships with human diseases including AV.

Local similarity of activity patterns during auditory and visual processing

Neuroimaging studies have shown that brain activity is variable and changes according to stimuli and the environmental context, reflecting brain coding or information representations at different processing levels. However, little is known about activity organization that reflects coding strategies. Here, we explored and compared two different coding approaches, spatial via cross-correlation and intensity-based coding using mutual information. Using two fMRI datasets and different seeds, we searched for the spatial and intensity-based similarities with the seeds in brain activity. Our results showed that, apart from the seed regions, significant regions detected by intensity-based similarity analysis differ completely from those found using cross-correlation. These findings may indicate that information shared through spatial coding differs from that transmitted via non-spatial coding processes. Our results suggest that brain coding is organized in several different ways to optimize information processing.

A study on the probabilistic quantification of heavy-truck crash risk under the influence of multi-factors

In order to manage and reduce a specific risk, its quantitative analysis is necessary. The key objective of this study is to illustrate the prevalence of multi-factors in fatal crashes involving heavy trucks and to quantify the crash risk under the influence of multi-factors. Data from a recent, nationally representative sample of fatal crashes was investigated to identify the risk factors contributing to crash causations and a novel risk index was obtained to develop a criterion for crash risk quantification. Then, based on the mutual information theory, the mutual dependence between risk factors was calculated to quantify the crash risk under different risk factor combinations. The results reveal that most heavy-truck fatal crashes are the result of co-occurring multi-factors rather than a single factor, and are mainly caused by simultaneous occurrence of two or three contributing factors. Moreover, crash risk increases with the increase of the number of risk factors influencing the driver. Furthermore, multi-factor interaction between certain risk factors, such as environmental and vehicular factors, makes incremental contribution to the crash risk by further increasing the crash probability. Specifically, when driver's aberrant behaviors (errors and/or violations) are exposed to both environmental and vehicular factors, driver's likelihood of being involved in a fatal crash increases significantly. This suggests that in addition to the number of risk factors, the crash risk also depends on the multi-factor interactions between different risk factors. Therefore, the effects of individual risk factors should be controlled at the outset to prevent the incremental effects of multi-factors on crash risk, in turn enabling risk minimization..

Data augmentation for depression detection using skeleton-based gait information

In recent years, the incidence of depression is rising rapidly worldwide, but large-scale depression screening is still challenging. Gait analysis provides a non-contact, low-cost, and efficient early screening method for depression. However, the early screening of depression based on gait analysis lacks sufficient effective sample data. In this paper, we propose a skeleton data augmentation method for assessing the risk of depression. First, we propose five techniques to augment skeleton data and apply them to depression and emotion datasets. Then, we divide augmentation methods into two types (non-noise augmentation and noise augmentation) based on the mutual information and the classification accuracy. Finally, we explore which augmentation strategies can capture the characteristics of human skeleton data more effectively. Experimental results show that the augmented training dataset that retains more of the raw skeleton data properties determines the performance of the detection model. Specifically, rotation augmentation and channel mask augmentation make the depression detection accuracy reach 92.15% and 91.34%, respectively.

Epigenetic regulation of fetal brain development in pig

How fetal brain development is regulated at the molecular level is not well understood. Due to ethical challenges associated with research on the human fetus, large animals particularly pigs are increasingly used to study development and disorders of fetal brain. The pig fetal brain grows rapidly during the last ∼ 50 days before birth which is around day 60 (d60) of pig gestation. But what regulates the onset of accelerated growth of the brain is unknown. The current study tests the hypothesis that epigenetic alteration around d60 is involved in the onset of rapid growth of fetal brain of pig. To test this hypothesis, DNA methylation changes of fetal brain was assessed in a genome-wide manner by Enzymatic Methyl-seq (EM-seq) during two gestational periods (GP): d45 vs. d60 (GP1) and d60 vs. d90 (GP2). The cytosine-guanine (CpG) methylation data was analyzed in an integrative manner with the RNA-seq data generated from the same brain samples from our earlier study. A neural network based modeling approach was implemented to learn changes in methylation patterns of the differentially expressed genes, and then predict methylations of the brain in a genome-wide manner during rapid growth. This approach identified specific methylations that changed in a mutually informative manner during rapid growth of the fetal brain. These methylations were significantly overrepresented in specific genic as well as intergenic features including CpG islands, introns, and untranslated regions. In addition, sex-bias methylations of known single nucleotide polymorphic sites were also identified in the fetal brain ide during rapid growth.

Detecting Unbiased Associations in Large Data Sets

Maximal information coefficient (MIC) explores the associations between pairwise variables in complex relationships. It approaches the correlation by optimized partition on the axis. However, when the relationships meet special noise, MIC may overestimate the correlated value, which leads to the misidentification of the relationship without noiseless. In this article, a novel method of weighted information coefficient mean (WICM) is proposed to detect unbiased associations in large data sets. First, we mathematically analyze the cause of giving an abnormal correlation value to a noisy relationship. Then, the WICM is presented in two core steps. One is to detect the potential overestimation from the relationships with high value, and the other is to rectify the overestimation by calculating information coefficient mean instead of just selecting the maximum element in the characteristic matrix. Finally, experiments in functional relationships and real-world data relationships show that the overestimation can be solved by WICM with both feasibility and effectiveness.

Gene expression data analysis using Hellinger correlation in weighted gene co-expression networks (WGCNA)

Weighted gene co-expression network analysis (WGCNA) is used to detect clusters with highly correlated genes. Measurements of correlation most typically rely on linear relationships. However, a linear relationship does not always model pairwise functional-related dependence between genes. In this paper, we first compared 6 different correlation methods in their ability to capture complex dependence between genes in three different tissues. Next, we compared their gene-pairwise coefficient results and corresponding WGCNA results. Finally, we applied a recently proposed correlation method, Hellinger correlation, as a more sensitive correlation measurement in WGCNA. To test this method, we constructed gene networks containing co-expression gene modules from RNA-seq data of human frontal cortex from Alzheimer's disease patients. To test the generality, we also used a microarray data set from human frontal cortex, single cell RNA-seq data from human prefrontal cortex, RNA-seq data from human temporal cortex, and GTEx data from heart. The Hellinger correlation method captures essentially similar results as other linear correlations in WGCNA, but provides additional new functional relationships as exemplified by uncovering a link between inflammation and mitochondria function. We validated the network constructed with the microarray and single cell sequencing data sets and a RNA-seq dataset of temporal cortex. We observed that this new correlation method enables the detection of non-linear biologically meaningful relationships among genes robustly and provides a complementary new approach to WGCNA. Thus, the application of Hellinger correlation to WGCNA provides a more flexible correlation approach to modelling networks in gene expression analysis that uncovers novel network relationships.

Approaches in Gene Coexpression Analysis in Eukaryotes

Simple Summary

Genes whose expression levels rise and fall similarly in a large set of samples, may be considered coexpressed. Gene coexpression analysis refers to the en masse discovery of coexpressed genes from a large variety of transcriptomic experiments. The type of biological networks that studies gene coexpression, known as Gene Coexpression Networks, consist of an undirected graph depicting genes and their coexpression relationships. Coexpressed genes are clustered in smaller subnetworks, the predominant biological roles of which can be determined through enrichment analysis. By studying well-annotated gene partners, the attribution of new roles to genes of unknown function or assumption for participation in common metabolic pathways can be achieved, through a guilt-by-association approach. In this review, we present key issues in gene coexpression analysis, as well as the most popular tools that perform it.

Abstract

Gene coexpression analysis constitutes a widely used practice for gene partner identification and gene function prediction, consisting of many intricate procedures. The analysis begins with the collection of primary transcriptomic data and their preprocessing, continues with the calculation of the similarity between genes based on their expression values in the selected sample dataset and results in the construction and visualisation of a gene coexpression network (GCN) and its evaluation using biological term enrichment analysis. As gene coexpression analysis has been studied extensively, we present most parts of the methodology in a clear manner and the reasoning behind the selection of some of the techniques. In this review, we offer a comprehensive and comprehensible account of the steps required for performing a complete gene coexpression analysis in eukaryotic organisms. We comment on the use of RNA-Seq vs. microarrays, as well as the best practices for GCN construction. Furthermore, we recount the most popular webtools and standalone applications performing gene coexpression analysis, with details on their methods, features and outputs.

Enhancer methylation dynamics drive core transcriptional regulatory circuitry in pan-cancer

Accumulating evidence has demonstrated that enhancer methylation has strong and dynamic regulatory effects on gene expression. Some transcription factors (TFs) can auto- and cross-regulate in a feed-forward manner, and cooperate with their enhancers to form core transcriptional regulatory circuitries (CRCs). However, the elaborated regulatory mechanism between enhancer methylation and CRC remains the tip of the iceberg. Here, we revealed that DNA methylation could drive the tissue-specific enhancer basal transcription and target gene expression in human cancers. By integrating methylome, transcriptome, and 3D genomic data, we identified enhancer methylation triplets (enhancer methylation-enhancer transcription-target gene expression) and dissected potential regulatory patterns within them. Moreover, we observed that cancer-specific core TFs regulated by enhancers were able to shape their enhancer methylation forming the enhancer methylation-driven CRCs (emCRCs). Further parsing of clinical implications showed rewired emCRCs could serve as druggable targets and prognostic risk markers. In summary, the integrative analysis of enhancer methylation regulome would facilitate portraying the cancer epigenomics landscape and developing the epigenetic anti-cancer approaches.

Fast and accurate inference of gene regulatory networks through robust precision matrix estimation

MOTIVATION

Transcriptional regulation mechanisms allow cells to adapt and respond to external stimuli by altering gene expression. The possible cell transcriptional states are determined by the underlying gene regulatory network (GRN), and reliably inferring such network would be invaluable to understand biological processes and disease progression.

RESULTS

In this article, we present a novel method for the inference of GRNs, called PORTIA, which is based on robust precision matrix estimation, and we show that it positively compares with state-of-the-art methods while being orders of magnitude faster. We extensively validated PORTIA using the DREAM and MERLIN+P datasets as benchmarks. In addition, we propose a novel scoring metric that builds on graph-theoretical concepts.

AVAILABILITY AND IMPLEMENTATION

The code and instructions for data acquisition and full reproduction of our results are available at https://github.com/AntoinePassemiers/PORTIA-Manuscript. PORTIA is available on PyPI as a Python package (portia-grn).

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

SNMRS: An advanced measure for Co-expression network analysis

The challenge of identifying modules in a gene interaction network is important for a better understanding of the overall network architecture. In this work, we develop a novel similarity measure called Scaling-and-Shifting Normalized Mean Residue Similarity (SNMRS), based on the existing NMRS technique [1]. SNMRS yields correlation values in the range of 0 to +1 corresponding to negative and positive dependency. To study the performance of our measure, internal validation of extracted clusters resulting from different methods is carried out. Based on the performance, we choose hierarchical clustering and apply the same using the corresponding dissimilarity (distance) values of SNMRS scores, and utilize a dynamic tree cut method for extracting dense modules. The modules are validated using a literature search, KEGG pathway analysis, and gene-ontology analyses on the genes that make up the modules. Moreover, our measure can handle absolute, shifting, scaling, and shifting-and-scaling correlations and provides better performance than several other measures in terms of cluster-validity indices. Also, SNMRS based module detection method results in interesting biologically relevant patterns from gene microarray and RNA-seq dataset. A set of crucial genes having high relevance with the ESCC are also identified.