Cerebellar Neurodynamics Predict Decision Timing and Outcome on the Single-Trial Level.

Goal-directed behavior requires the interaction of multiple brain regions. How these regions and their interactions with brain-wide activity drive action selection is less understood. We have investigated this question by combining whole-brain volumetric calcium imaging using light-field microscopy and an operant-conditioning task in larval zebrafish. We find global, recurring dynamics of brain states to exhibit pre-motor bifurcations toward mutually exclusive decision outcomes. These dynamics arise from a distributed network displaying trial-by-trial functional connectivity changes, especially between cerebellum and habenula, which correlate with decision outcome. Within this network the cerebellum shows particularly strong and predictive pre-motor activity (>10 s before movement initiation), mainly within the granule cells. Turn directions are determined by the difference neuroactivity between the ipsilateral and contralateral hemispheres, while the rate of bi-hemispheric population ramping quantitatively predicts decision time on the trial-by-trial level. Our results highlight a cognitive role of the cerebellum and its importance in motor planning.

A scalable and robust variance components method reveals insights into the architecture of gene-environment interactions underlying complex traits.

Understanding the contribution of gene-environment interactions (GxE) to complex trait variation can provide insights into disease mechanisms, explain sources of heritability, and improve genetic risk prediction. While large biobanks with genetic and deep phenotypic data hold promise for obtaining novel insights into GxE, our understanding of GxE architecture in complex traits remains limited. We introduce a method to estimate the proportion of trait variance explained by GxE (GxE heritability) and additive genetic effects (additive heritability) across the genome and within specific genomic annotations. We show that our method is accurate in simulations and computationally efficient for biobank-scale datasets. We applied our method to common array SNPs (MAF ≥1%), fifty quantitative traits, and four environmental variables (smoking, sex, age, and statin usage) in unrelated white British individuals in the UK Biobank. We found 68 trait-E pairs with significant genome-wide GxE heritability (p<0.05/200) with a ratio of GxE to additive heritability of ≈6.8% on average. Analyzing ≈8 million imputed SNPs (MAF ≥0.1%), we documented an approximate 28% increase in genome-wide GxE heritability compared to array SNPs. We partitioned GxE heritability across minor allele frequency (MAF) and local linkage disequilibrium (LD) values, revealing that, like additive allelic effects, GxE allelic effects tend to increase with decreasing MAF and LD. Analyzing GxE heritability near genes highly expressed in specific tissues, we find significant brain-specific enrichment for body mass index (BMI) and basal metabolic rate in the context of smoking and adipose-specific enrichment for waist-hip ratio (WHR) in the context of sex.

The training process of many deep networks explores the same low-dimensional manifold.

We develop information-geometric techniques to analyze the trajectories of the predictions of deep networks during training. By examining the underlying high-dimensional probabilistic models, we reveal that the training process explores an effectively low-dimensional manifold. Networks with a wide range of architectures, sizes, trained using different optimization methods, regularization techniques, data augmentation techniques, and weight initializations lie on the same manifold in the prediction space. We study the details of this manifold to find that networks with different architectures follow distinguishable trajectories, but other factors have a minimal influence; larger networks train along a similar manifold as that of smaller networks, just faster; and networks initialized at very different parts of the prediction space converge to the solution along a similar manifold.

Specialized proresolving mediator resolvin E1 corrects the altered cystic fibrosis nasal epithelium cilia beating dynamics.

In cystic fibrosis (CF), impaired mucociliary clearance leads to chronic infection and inflammation. However, cilia beating features in a CF altered environment, consisting of dehydrated airway surface liquid layer and abnormal mucus, have not been fully characterized. Furthermore, acute inflammation is normally followed by an active resolution phase requiring specialized proresolving lipid mediators (SPMs) and allowing return to homeostasis. However, altered SPMs biosynthesis has been reported in CF. Here, we explored cilia beating dynamics in CF airways primary cultures and its response to the SPMs, resolvin E1 (RvE1) and lipoxin B4 (LXB4). Human nasal epithelial cells (hNECs) from CF and non-CF donors were grown at air-liquid interface. The ciliary beat frequency, synchronization, orientation, and density were analyzed from high-speed video microscopy using a multiscale Differential Dynamic Microscopy algorithm and an in-house developed method. Mucins and ASL layer height were studied by qRT-PCR and confocal microscopy. Principal component analysis showed that CF and non-CF hNEC had distinct cilia beating phenotypes, which was mostly explained by differences in cilia beat organization rather than frequency. Exposure to RvE1 (10 nM) and to LXB4 (10 nM) restored a non-CF-like cilia beating phenotype. Furthermore, RvE1 increased the airway surface liquid (ASL) layer height and reduced the mucin MUC5AC thickness. The calcium-activated chloride channel, TMEM16A, was involved in the RvE1 effect on cilia beating, hydration, and mucus. Altogether, our results provide evidence for defective cilia beating in CF airway epithelium and a role of RvE1 and LXB4 to restore the main epithelial functions involved in the mucociliary clearance.

Multilevel superposition for deciphering the conformational variability of protein ensembles.

The dynamics and variability of protein conformations are directly linked to their functions. Many comparative studies of X-ray protein structures have been conducted to elucidate the relevant conformational changes, dynamics and heterogeneity. The rapid increase in the number of experimentally determined structures has made comparison an effective tool for investigating protein structures. For example, it is now possible to compare structural ensembles formed by enzyme species, variants or the type of ligands bound to them. In this study, the author developed a multilevel model for estimating two covariance matrices that represent inter- and intra-ensemble variability in the Cartesian coordinate space. Principal component analysis using the two estimated covariance matrices identified the inter-/intra-enzyme variabilities, which seemed to be important for the enzyme functions, with the illustrative examples of cytochrome P450 family 2 enzymes and class A $\beta$-lactamases. In P450, in which each enzyme has its own active site of a distinct size, an active-site motion shared universally between the enzymes was captured as the first principal mode of the intra-enzyme covariance matrix. In this case, the method was useful for understanding the conformational variability after adjusting for the differences between enzyme sizes. The developed method is advantageous in small ensemble-size problems and hence promising for use in comparative studies on experimentally determined structures where ensemble sizes are smaller than those generated, for example, by molecular dynamics simulations.

Projection-based two-sample inference for sparsely observed multivariate functional data.

Modern longitudinal studies collect multiple outcomes as the primary endpoints to understand the complex dynamics of the diseases. Oftentimes, especially in clinical trials, the joint variation among the multidimensional responses plays a significant role in assessing the differential characteristics between two or more groups, rather than drawing inferences based on a single outcome. We develop a projection-based two-sample significance test to identify the population-level difference between the multivariate profiles observed under a sparse longitudinal design. The methodology is built upon widely adopted multivariate functional principal component analysis to reduce the dimension of the infinite-dimensional multi-modal functions while preserving the dynamic correlation between the components. The test applies to a wide class of (non-stationary) covariance structures of the response, and it detects a significant group difference based on a single p-value, thereby overcoming the issue of adjusting for multiple p-values that arise due to comparing the means in each of components separately. Finite-sample numerical studies demonstrate that the test maintains the type-I error, and is powerful to detect significant group differences, compared to the state-of-the-art testing procedures. The test is carried out on two significant longitudinal studies for Alzheimer's disease and Parkinson's disease (PD) patients, namely, TOMMORROW study of individuals at high risk of mild cognitive impairment to detect differences in the cognitive test scores between the pioglitazone and the placebo groups, and Azillect study to assess the efficacy of rasagiline as a potential treatment to slow down the progression of PD.

Functional support vector machine.

Linear and generalized linear scalar-on-function modeling have been commonly used to understand the relationship between a scalar response variable (e.g. continuous, binary outcomes) and functional predictors. Such techniques are sensitive to model misspecification when the relationship between the response variable and the functional predictors is complex. On the other hand, support vector machines (SVMs) are among the most robust prediction models but do not take account of the high correlations between repeated measurements and cannot be used for irregular data. In this work, we propose a novel method to integrate functional principal component analysis with SVM techniques for classification and regression to account for the continuous nature of functional data and the nonlinear relationship between the scalar response variable and the functional predictors. We demonstrate the performance of our method through extensive simulation experiments and two real data applications: the classification of alcoholics using electroencephalography signals and the prediction of glucobrassicin concentration using near-infrared reflectance spectroscopy. Our methods especially have more advantages when the measurement errors in functional predictors are relatively large.

A comparison between similarity matrices for principal component analysis to assess population stratification in sequenced genetic data sets.

Genetic similarity matrices are commonly used to assess population substructure (PS) in genetic studies. Through simulation studies and by the application to whole-genome sequencing (WGS) data, we evaluate the performance of three genetic similarity matrices: the unweighted and weighted Jaccard similarity matrices and the genetic relationship matrix. We describe different scenarios that can create numerical pitfalls and lead to incorrect conclusions in some instances. We consider scenarios in which PS is assessed based on loci that are located across the genome ('globally') and based on loci from a specific genomic region ('locally'). We also compare scenarios in which PS is evaluated based on loci from different minor allele frequency bins: common (>5%), low-frequency (5-0.5%) and rare (<0.5%) single-nucleotide variations (SNVs). Overall, we observe that all approaches provide the best clustering performance when computed based on rare SNVs. The performance of the similarity matrices is very similar for common and low-frequency variants, but for rare variants, the unweighted Jaccard matrix provides preferable clustering features. Based on visual inspection and in terms of standard clustering metrics, its clusters are the densest and the best separated in the principal component analysis of variants with rare SNVs compared with the other methods and different allele frequency cutoffs. In an application, we assessed the role of rare variants on local and global PS, using WGS data from multiethnic Alzheimer's disease data sets and European or East Asian populations from the 1000 Genome Project.

Identifying potential small molecule-miRNA associations via Robust PCA based on γ-norm regularization.

Dysregulation of microRNAs (miRNAs) is closely associated with refractory human diseases, and the identification of potential associations between small molecule (SM) drugs and miRNAs can provide valuable insights for clinical treatment. Existing computational techniques for inferring potential associations suffer from limitations in terms of accuracy and efficiency. To address these challenges, we devise a novel predictive model called RPCA$\Gamma $NR, in which we propose a new Robust principal component analysis (PCA) framework based on $\gamma $-norm and $l_{2,1}$-norm regularization and design an Augmented Lagrange Multiplier method to optimize it, thereby deriving the association scores. The Gaussian Interaction Profile Kernel Similarity is calculated to capture the similarity information of SMs and miRNAs in known associations. Through extensive evaluation, including Cross Validation Experiments, Independent Validation Experiment, Efficiency Analysis, Ablation Experiment, Matrix Sparsity Analysis, and Case Studies, RPCA$\Gamma $NR outperforms state-of-the-art models concerning accuracy, efficiency and robustness. In conclusion, RPCA$\Gamma $NR can significantly streamline the process of determining SM-miRNA associations, thus contributing to advancements in drug development and disease treatment.

Integrating multiple traits for improving polygenic risk prediction in disease and pharmacogenomics GWAS.

Polygenic risk score (PRS) has been recently developed for predicting complex traits and drug responses. It remains unknown whether multi-trait PRS (mtPRS) methods, by integrating information from multiple genetically correlated traits, can improve prediction accuracy and power for PRS analysis compared with single-trait PRS (stPRS) methods. In this paper, we first review commonly used mtPRS methods and find that they do not directly model the underlying genetic correlations among traits, which has been shown to be useful in guiding multi-trait association analysis in the literature. To overcome this limitation, we propose a mtPRS-PCA method to combine PRSs from multiple traits with weights obtained from performing principal component analysis (PCA) on the genetic correlation matrix. To accommodate various genetic architectures covering different effect directions, signal sparseness and across-trait correlation structures, we further propose an omnibus mtPRS method (mtPRS-O) by combining P values from mtPRS-PCA, mtPRS-ML (mtPRS based on machine learning) and stPRSs using Cauchy Combination Test. Our extensive simulation studies show that mtPRS-PCA outperforms other mtPRS methods in both disease and pharmacogenomics (PGx) genome-wide association studies (GWAS) contexts when traits are similarly correlated, with dense signal effects and in similar effect directions, and mtPRS-O is consistently superior to most other methods due to its robustness under various genetic architectures. We further apply mtPRS-PCA, mtPRS-O and other methods to PGx GWAS data from a randomized clinical trial in the cardiovascular domain and demonstrate performance improvement of mtPRS-PCA in both prediction accuracy and patient stratification as well as the robustness of mtPRS-O in PRS association test.

Cooperation of structural motifs controls drug selectivity in cyclin-dependent kinases: an advanced theoretical analysis.

Understanding drug selectivity mechanism is a long-standing issue for helping design drugs with high specificity. Designing drugs targeting cyclin-dependent kinases (CDKs) with high selectivity is challenging because of their highly conserved binding pockets. To reveal the underlying general selectivity mechanism, we carried out comprehensive analyses from both the thermodynamics and kinetics points of view on a representative CDK12 inhibitor. To fully capture the binding features of the drug-target recognition process, we proposed to use kinetic residue energy analysis (KREA) in conjunction with the community network analysis (CNA) to reveal the underlying cooperation effect between individual residues/protein motifs to the binding/dissociating process of the ligand. The general mechanism of drug selectivity in CDKs can be summarized as that the difference of structural cooperation between the ligand and the protein motifs leads to the difference of the energetic contribution of the key residues to the ligand. The proposed mechanisms may be prevalent in drug selectivity issues, and the insights may help design new strategies to overcome/attenuate the drug selectivity associated problems.

The interictal suppression hypothesis is the dominant differentiator of seizure onset zones in focal epilepsy.

Successful surgical treatment of drug-resistant epilepsy traditionally relies on the identification of seizure onset zones (SOZs). Connectome-based analyses of electrographic data from stereo electroencephalography (SEEG) may empower improved detection of SOZs. Specifically, connectome-based analyses based on the interictal suppression hypothesis posit that when the patient is not having a seizure, SOZs are inhibited by non-SOZs through high inward connectivity and low outward connectivity. However, it is not clear whether there are other motifs that can better identify potential SOZs. Thus, we sought to use unsupervised machine learning to identify network motifs that elucidate SOZs and investigate if there is another motif that outperforms the ISH. Resting-state SEEG data from 81 patients with drug-resistant epilepsy undergoing a pre-surgical evaluation at Vanderbilt University Medical Center were collected. Directed connectivity matrices were computed using the alpha band (8-13 Hz). Principal component analysis (PCA) was performed on each patient's connectivity matrix. Each patient's components were analysed qualitatively to identify common patterns across patients. A quantitative definition was then used to identify the component that most closely matched the observed pattern in each patient. A motif characteristic of the interictal suppression hypothesis (high-inward and low-outward connectivity) was present in all individuals and found to be the most robust motif for identification of SOZs in 64/81 (79%) patients. This principal component demonstrated significant differences in SOZs compared to non-SOZs. While other motifs for identifying SOZs were present in other patients, they differed for each patient, suggesting that seizure networks are patient specific, but the ISH is present in nearly all networks. We discovered that a potentially suppressive motif based on the interictal suppression hypothesis was present in all patients, and it was the most robust motif for SOZs in 79% of patients. Each patient had additional motifs that further characterized SOZs, but these motifs were not common across all patients. This work has the potential to augment clinical identification of SOZs to improve epilepsy treatment.

VCF2PCACluster: a simple, fast and memory-efficient tool for principal component analysis of tens of millions of SNPs.

Principal component analysis (PCA) is an important and widely used unsupervised learning method that determines population structure based on genetic variation. Genome sequencing of thousands of individuals usually generate tens of millions of SNPs, making it challenging for PCA analysis and interpretation. Here we present VCF2PCACluster, a simple, fast and memory-efficient tool for Kinship estimation, PCA and clustering analysis, and visualization based on VCF formatted SNPs. We implemented five Kinship estimation methods and three clustering methods for its users to choose from. Moreover, unlike other PCA tools, VCF2PCACluster possesses a clustering function based on PCA result, which enabling users to automatically and clearly know about population structure. We demonstrated the same accuracy but a higher performance of this tool in performing PCA analysis on tens of millions of SNPs compared to another popular PLINK2 software, especially in peak memory usage that is independent of the number of SNPs in VCF2PCACluster.

Exploring combinations of dimensionality reduction, transfer learning, and regularization methods for predicting binary phenotypes with transcriptomic data.

BACKGROUND: Numerous transcriptomic-based models have been developed to predict or understand the fundamental mechanisms driving biological phenotypes. However, few models have successfully transitioned into clinical practice due to challenges associated with generalizability and interpretability. To address these issues, researchers have turned to dimensionality reduction methods and have begun implementing transfer learning approaches. METHODS: In this study, we aimed to determine the optimal combination of dimensionality reduction and regularization methods for predictive modeling. We applied seven dimensionality reduction methods to various datasets, including two supervised methods (linear optimal low-rank projection and low-rank canonical correlation analysis), two unsupervised methods [principal component analysis and consensus independent component analysis (c-ICA)], and three methods [autoencoder (AE), adversarial variational autoencoder, and c-ICA] within a transfer learning framework, trained on > 140,000 transcriptomic profiles. To assess the performance of the different combinations, we used a cross-validation setup encapsulated within a permutation testing framework, analyzing 30 different transcriptomic datasets with binary phenotypes. Furthermore, we included datasets with small sample sizes and phenotypes of varying degrees of predictability, and we employed independent datasets for validation. RESULTS: Our findings revealed that regularized models without dimensionality reduction achieved the highest predictive performance, challenging the necessity of dimensionality reduction when the primary goal is to achieve optimal predictive performance. However, models using AE and c-ICA with transfer learning for dimensionality reduction showed comparable performance, with enhanced interpretability and robustness of predictors, compared to models using non-dimensionality-reduced data. CONCLUSION: These findings offer valuable insights into the optimal combination of strategies for enhancing the predictive performance, interpretability, and generalizability of transcriptomic-based models.

Analyzing postprandial metabolomics data using multiway models: a simulation study.

BACKGROUND: Analysis of time-resolved postprandial metabolomics data can improve the understanding of metabolic mechanisms, potentially revealing biomarkers for early diagnosis of metabolic diseases and advancing precision nutrition and medicine. Postprandial metabolomics measurements at several time points from multiple subjects can be arranged as a subjects by metabolites by time points array. Traditional analysis methods are limited in terms of revealing subject groups, related metabolites, and temporal patterns simultaneously from such three-way data. RESULTS: We introduce an unsupervised multiway analysis approach based on the CANDECOMP/PARAFAC (CP) model for improved analysis of postprandial metabolomics data guided by a simulation study. Because of the lack of ground truth in real data, we generate simulated data using a comprehensive human metabolic model. This allows us to assess the performance of CP models in terms of revealing subject groups and underlying metabolic processes. We study three analysis approaches: analysis of fasting-state data using principal component analysis, T0-corrected data (i.e., data corrected by subtracting fasting-state data) using a CP model and full-dynamic (i.e., full postprandial) data using CP. Through extensive simulations, we demonstrate that CP models capture meaningful and stable patterns from simulated meal challenge data, revealing underlying mechanisms and differences between diseased versus healthy groups. CONCLUSIONS: Our experiments show that it is crucial to analyze both fasting-state and T0-corrected data for understanding metabolic differences among subject groups. Depending on the nature of the subject group structure, the best group separation may be achieved by CP models of T0-corrected or full-dynamic data. This study introduces an improved analysis approach for postprandial metabolomics data while also shedding light on the debate about correcting baseline values in longitudinal data analysis.

A novel method for cell deconvolution using DNA methylation in PCA space.

BACKGROUND: In this study, we present a novel method for reference-based cell deconvolution using data from DNA methylation arrays. Different from existing methods like IDOL-Ext, which operate on probe-level data, our approach represents features in the principal component analysis (PCA) space for cell type deconvolution. RESULTS: Our method's accuracy in estimating cell compositions is validated across various public datasets, including blood samples from glioma patients. It demonstrates precision comparable to IDOL-Ext, with R2 values ranging from 0.73 to 0.99 for most cell types, while offering improved discrimination between similar cell types, particularly T cell subtypes in glioma patient samples (R2 0.42-0.75 vs. 0.36-0.66 for IDOL-Ext). However, both methods showed lower accuracy for certain cell types, such as memory CD8 T cells in glioma patients (R2 0.42 vs. 0.36 for IDOL-Ext), highlighting the challenges in distinguishing closely related cell populations. We have made this method available as an R package "BloodCellDecon" on GitHub. CONCLUSIONS: Our study confirms the efficacy of cell type deconvolution in PCA space. The results indicate wide-ranging applicability and potential for adaptation to other forms of genomic data.

Altered static and dynamic functional brain network in knee osteoarthritis: A resting-state functional magnetic resonance imaging study: Static and dynamic FNC in KOA.

This study aimed to investigate altered static and dynamic functional network connectivity (FNC) and its correlation with clinical symptoms in patients with knee osteoarthritis (KOA). One hundred and fifty-nine patients with KOA and 73 age- and gender-matched healthy subjects (HS) underwent resting-state functional magnetic resonance imaging (rs-fMRI) and clinical evaluations. Group independent component analysis (GICA) was applied, and seven resting-state networks were identified. Patients with KOA had decreased static FNC within the default mode network (DM), visual network (VS), and cerebellar network (CB) and increased static FNC between the subcortical network (SC) and VS (p < 0.05, FDR corrected). Four reoccurring FNC states were identified using k-means clustering analysis. Although abnormalities in dynamic FNCs of KOA patients have been found using the common window size (22 TR, 44 s), but the results of the clustering analysis were inconsistent when using different window sizes, suggesting dynamic FNCs might be an unstable method to compare brain function between KOA patients and HS. These recent findings illustrate that patients with KOA have a wide range of abnormalities in the static and dynamic FNCs, which provided a reference for the identification of potential central nervous therapeutic targets for KOA treatment and might shed light on the other musculoskeletal pain neuroimaging studies.

Free energy profile of the substrate-induced occlusion of the human serotonin transporter.

The serotonin transporter (SERT) is a member of the Solute Carrier 6 (SLC6) family and is responsible for maintaining the appropriate level of serotonin in the brain. Dysfunction of SERT has been linked to several neuropsychiatric disorders, including depression, anxiety and obsessive-compulsive disorder. Therefore, an in-depth understanding of the mechanism on an atomistic level, coupled with a quantification of transporter dynamics and the associated free energies is required. Here, we constructed Markov state models (MSMs) from extensive unbiased molecular dynamics simulations to quantify the free energy profile of serotonin (5HT) triggered SERT occlusion and explored the driving forces of the mechanism of occlusion. Our results reveal that SERT occludes via multiple intermediate conformations and show that the motion of occlusion is energetically downhill for the 5HT-bound transporter. Force distribution analyses show that the interactions of 5HT with the bundle domain are crucial. During occlusion, attractive forces steadily increase and pull on the bundle domain, which leads to SERT occlusion. Some interactions become repulsive upon full occlusion, suggesting that SERT creates pressure on 5HT to promote its movement towards the cytosol.

Population-genetic and comparative interpopulation studies of the 15 autosomal STR markers in the population living in the Northwest of Iran.

INTRODUCTION: Iran, a country in the Middle East, has several ethnic and ethno-religious groups and needs its own ethnic-specific databases for the forensic statistical parameters and allele frequency of STR markers. METHODS: We have investigated 600 unrelated Turk individuals from four northwestern provinces of Iran using the Identifiler™ system (TPOX, FGA, vWA, TH01, CSF1PO, D2S1338, D3S1358, D5S818, D7S820, D8S1179, D13S317, D16S539, D18S51, D19S433, and D21S11). Furthermore, STR allelic frequencies were compared to previously population-based data. RESULTS AND CONCLUSION: After Bonferroni correction, deviation from Hardy-Weinberg equilibrium (HWE) was observed in the FGA, TPOX, VWA, and D19S433 loci (P value < 0.05). The combined power of discrimination (CPD) and exclusion (CPE) values for all 15 STR loci were 0.9999999999999999999984 and 0.9999999, respectively. In comparison with Azerbaijani and Turkish populations, there were no significant differences on all STR markers. However, in the Chinese Han population, differences at 13 STR loci were detected. Additionally, comparisons of Fischer genetic distance indices (FST) P-values did not reveal any statistically significant difference between Northwestern Iran, Azerbaijan and Iran (Fars) populations. PCA and PCoA analyses showed that our population was grouped with different populations in different quarters, showing a positive and negative correlation, respectively. In the NJ and UPGMA phylogenetic trees, Iranian populations were grouped together. These results demonstrated that the given set of STR markers can be confidently used for all identification tests in Northwestern Iran.

Alterations of resting-state functional network connectivity in patients with noise-induced hearing loss: A study based on independent component analysis.

Functional reorganization is a response to auditory deficits or deprivation, and less is known about the overall brain network alterations involving resting-state networks (RSNs) and multiple functional networks in patients with occupational noise-induced hearing loss (NIHL). So this study evaluated resting-state functional network connectivity (FNC) alterations in occupational NIHL using an independent component analysis (ICA). In total, 79 mild NIHL patients (MP), 32 relatively severe NIHL patients (RSP), and 84 age- and education- matched healthy controls (HC) were recruited. All subjects were tested using the Mini-mental State Examination scale, the tinnitus Handicap Inventory scale, the Hamilton Anxiety scale (HAMA) and scanned by T1-3DFSPGR, resting-state functional magnetic resonance imaging sequence in 3.0 T and analysed by the ICA. Seven RSNs were identified, compared with the HC, the MP showed increased FNC within the executive control network (ECN) and enhanced FNC within the default mode network (DMN) and the visual network (VN); compared with the HC, the RSP showed decreased FNC within the ECN and auditory network (AUN), DMN and VN; no significant changes in FNC were found in the MP compared with the RSP. Furthermore, the correlation analysis between the noise exposure time and hearing loss level, HAMA were both negative, and there were no significant correlations between the abnormal RSNs and the hearing level, noise exposure time and HAMA. These findings indicate that different degrees of NIHL involve different alterations in RSNs connectivity and may reveal the neural mechanisms related to emotion-related features and functional abnormalities following long-term NIHL.