Search

AimsCalcific aortic valve disease (CAVD) is the most common valve disease, which consists of a chronic interplay of inflammation, fibrosis, and calcification. In this study, sortilin (SORT1) was identified as a novel key player in the pathophysiology of CAVD, and its role in the transformation of valvular interstitial cells (VICs) into pathological phenotypes is explored.Methods and resultsAn aortic valve (AV) wire injury (AVWI) mouse model with sortilin deficiency was used to determine the effects of sortilin on AV stenosis, fibrosis, and calcification. In vitro experiments employed human primary VICs cultured in osteogenic conditions for 7, 14, and 21 days; and processed for imaging, proteomics, and transcriptomics including single-cell RNA-sequencing (scRNA-seq). The AVWI mouse model showed reduced AV fibrosis, calcification, and stenosis in sortilin-deficient mice vs. littermate controls. Protein studies identified the transition of human VICs into a myofibroblast-like phenotype mediated by sortilin. Sortilin loss-of-function decreased in vitro VIC calcification. ScRNA-seq identified 12 differentially expressed cell clusters in human VIC samples, where a novel combined inflammatory myofibroblastic-osteogenic VIC (IMO-VIC) phenotype was detected with increased expression of SORT1, COL1A1, WNT5A, IL-6, and serum amyloid A1. VICs sequenced with sortilin deficiency showed decreased IMO-VIC phenotype.ConclusionSortilin promotes CAVD by mediating valvular fibrosis and calcification, and a newly identified phenotype (IMO-VIC). This is the first study to examine the role of sortilin in valvular calcification and it may render it a therapeutic target to inhibit IMO-VIC emergence by simultaneously reducing inflammation, fibrosis, and calcification, the three key pathological processes underlying CAVD.

Gene regulation plays a fundamental role in shaping tissue identity, function, and response to perturbation. Regulatory processes are controlled by complex networks of interacting elements, including transcription factors, miRNAs and their target genes. The structure of these networks helps to determine phenotypes and can ultimately influence the development of disease or response to therapy. We developed GRAND (https://grand.networkmedicine.org) as a database for computationally-inferred, context-specific gene regulatory network models that can be compared between biological states, or used to predict which drugs produce changes in regulatory network structure. The database includes 12 468 genome-scale networks covering 36 human tissues, 28 cancers, 1378 unperturbed cell lines, as well as 173 013 TF and gene targeting scores for 2858 small molecule-induced cell line perturbation paired with phenotypic information. GRAND allows the networks to be queried using phenotypic information and visualized using a variety of interactive tools. In addition, it includes a web application that matches disease states to potentially therapeutic small molecule drugs using regulatory network properties., Graphical Abstract Graphical AbstractModeling gene regulation across human conditions integrates cancer tissues and cell lines, small molecules and normal tissue networks.

ObjectiveTo use network science to model complex diet relationships a decade before onset of dementia in a large French cohort, the 3-City Bordeaux study.MethodsWe identified cases of dementia incident to the baseline food frequency questionnaire over 12 years of follow-up. For each case, we randomly selected 2 controls among individuals at risk at the age at case diagnosis and matched for age at diet assessment, sex, education, and season of the survey. We inferred food networks in both cases and controls using mutual information, a measure to detect nonlinear associations, and compared food consumption patterns between groups.ResultsIn the nested case-control study, the mean (SD) duration of follow-up and number of visits were 5.0 (2.5) vs 4.9 (2.6) years and 4.1 (1.0) vs 4.4 (0.9) for cases (n = 209) vs controls (n = 418), respectively. While there were few differences in simple, average food intakes, food networks differed substantially between cases and controls. The network in cases was focused and characterized by charcuterie as the main hub, with connections to foods typical of French southwestern diet and snack foods. In contrast, the network of controls included several disconnected subnetworks reflecting diverse and healthier food choices.ConclusionHow foods are consumed (and not only the quantity consumed) may be important for dementia prevention. Differences in predementia diet networks, suggesting worse eating habits toward charcuterie and snacking, were evident years before diagnosis in this cohort. Network methods, which are designed to model complex systems, may advance our understanding of risk factors for dementia.

Mass-spectrometry-enabled ADP-ribosylation workflows are developing rapidly, providing researchers a variety of ADP-ribosylome enrichment strategies and mass spectrometric acquisition options. Despite the growth spurt in upstream technologies, systematic ADP-ribosyl (ADPr) peptide mass spectral annotation methods are lacking. HCD-dependent ADP-ribosylome studies are common, but the resulting MS2 spectra are complex, owing to a mixture of b/y-ions and the m/p-ion peaks representing one or more dissociation events of the ADPr moiety (m-ion) and peptide (p-ion). In particular, p-ions that dissociate further into one or more fragment ions can dominate HCD spectra but are not recognized by standard spectral annotation workflows. As a result, annotation strategies that are solely reliant upon the b/y-ions result in lower spectral scores that in turn reduce the number of reportable ADPr peptides. To improve the confidence of spectral assignments, we implemented an ADPr peptide annotation and scoring strategy. All MS2 spectra are scored for the ADPr m-ions, but once spectra are assigned as an ADPr peptide, they are further annotated and scored for the p-ions. We implemented this novel workflow to ADPr peptides enriched from the liver and spleen isolated from mice post 4 h exposure to systemic IFN-γ. HCD collision energy experiments were first performed on the Orbitrap Fusion Lumos and the Q Exactive, with notable ADPr peptide dissociation properties verified with CID (Lumos). The m-ion and p-ion series score distributions revealed that ADPr peptide dissociation properties vary markedly between instruments and within instrument collision energy settings, with consequences on ADPr peptide reporting and amino acid localization. Consequentially, we increased the number of reportable ADPr peptides by 25% (liver) and 17% (spleen) by validation and the inclusion of lower confidence ADPr peptide spectra. This systematic annotation strategy will streamline future reporting of ADPr peptides that have been sequenced using any HCD/CID-based method.

Respiratory syncytial virus (RSV) is a major cause of lower respiratory tract infections and hospital visits during infancy and childhood. Although risk factors for RSV infection have been identified, the role of microbial species in the respiratory tract is only partially known. We aimed to understand the impact of interactions between the nasal microbiome and host transcriptome on the severity and clinical outcomes of RSV infection. We used 16 S rRNA sequencing to characterize the nasal microbiome of infants with RSV infection. We used RNA sequencing to interrogate the transcriptome of CD4+ T cells obtained from the same set of infants. After dimension reduction through principal component (PC) analysis, we performed an integrative analysis to identify significant co-variation between microbial clade and gene expression PCs. We then employed LIONESS (Linear Interpolation to Obtain Network Estimates for Single Samples) to estimate the clade-gene association patterns for each infant. Our network-based integrative analysis identified several clade-gene associations significantly related to the severity of RSV infection. The microbial taxa with the highest loadings in the implicated clade PCs included Moraxella, Corynebacterium, Streptococcus, Haemophilus influenzae, and Staphylococcus. Interestingly, many of the genes with the highest loadings in the implicated gene PCs are encoded in mitochondrial DNA, while others are involved in the host immune response. This study on microbiome-transcriptome interactions provides insights into how the host immune system mounts a response against RSV and specific infectious agents in nasal microbiota.

The biological processes that drive cellular function can be represented by a complex network of interactions between regulators (transcription factors) and their targets (genes). A cell’s epigenetic state plays an important role in mediating these interactions, primarily by influencing chromatin accessibility. However, effectively leveraging epigenetic information when constructing regulatory networks remains a challenge. We developed SPIDER, which incorporates epigenetic information (DNase-Seq) into a message passing framework in order to estimate gene regulatory networks. We validated SPIDER’s predictions using ChlP-Seq data from ENCODE and found that SPIDER networks were more accurate than other publicly available, epigenetically informed regulatory networks as well as networks based on methods that leverage epigenetic data to predict transcription factor binding sites. SPIDER was also able to improve the detection of cell line specific regulatory interactions. Notably, SPIDER can recover ChlP-seq verified transcription factor binding events in the regulatory regions of genes that do not have a corresponding sequence motif. Constructing biologically interpretable, epigenetically informed networks using SPIDER will allow us to better understand gene regulation as well as aid in the identification of cell-specific drivers and biomarkers of cellular phenotypes.

Background The International Agency for Research on Cancer classified radon and its decay-products as Group-1-human-carcinogens, and with the current knowledge they are linked specifically to lung cancer. Biokinetic models predict that radon could deliver a carcinogenic dose to breast tissue. Our previous work suggested that low-dose radon was associated with estrogen-receptor (ER)-negative breast cancer risk. However, there is limited research to examine the role of radon in breast cancer biology at the tissue level. We aim to understand molecular pathways linking radon exposure with breast cancer biology using transcriptome-wide-gene-expression from breast tumor and normal-adjacent tissues. Methods Our study included 943 women diagnosed with breast cancer from the Nurses’ Health Study (NHS) and NHSII. We estimated cumulative radon concentration for each participant up-to the year of breast cancer diagnosis by linking residential addresses with a radon exposure model. Transcriptome-wide-gene-expression was measured with the Affymetrix-Glue-Human-Transcriptome-Array-3.0 and Human-Transcriptome-Array-2.0. We performed covariate-adjusted linear-regression for individual genes and further employed pathway-analysis. All analyses were conducted separately for tumor and normal-adjacent samples and by ER-status. Results No individual gene was associated with cumulative radon exposure in ER-positive tumor, ER-negative tumor, or ER-negative normal-adjacent tissues at FDR PLCH2—reached transcriptome-wide-significance (FDR Conclusion Our findings provide insights into the molecular pathways of radon exposure that may influence breast cancer etiology.

Numerous studies have investigated individual biomarkers in relation to risk of type 2 diabetes. However, few have considered the interconnectivity of these biomarkers in the etiology of diabetes as well as the potential changes in the biomarker correlation network during diabetes development. We conducted a secondary analysis of 27 plasma biomarkers representing glucose metabolism, inflammation, adipokines, endothelial dysfunction, IGF axis, and iron store plus age and BMI at blood collection from an existing case-control study nested in the Nurses’ Health Study (NHS), including 1,303 incident diabetes case subjects and 1,627 healthy women. A correlation network was constructed based on pairwise Spearman correlations of the above factors that were statistically different between case and noncase subjects using permutation tests (P < 0.0005). We further evaluated the network structure separately among diabetes case subjects diagnosed 10 years after blood collection versus noncase subjects. Although pairwise biomarker correlations tended to have similar directions comparing diabetes case subjects to noncase subjects, most correlations were stronger in noncase than in case subjects, with the largest differences observed for the insulin/HbA1c and leptin/adiponectin correlations. Leptin and soluble leptin receptor were two hubs of the network, with large numbers of different correlations with other biomarkers in case versus noncase subjects. When examining the correlation network by timing of diabetes onset, there were more perturbations in the network for case subjects diagnosed >10 years versus

Sex differences manifest in many diseases and may drive sex-specific therapeutic responses. To understand the molecular basis of sex differences, we evaluated sex-biased gene regulation by constructing sample-specific gene regulatory networks in 29 human healthy tissues using 8,279 whole-genome expression profiles from the Genotype-Tissue Expression (GTEx) project. We find sex-biased regulatory network structures in each tissue. Even though most transcription factors (TFs) are not differentially expressed between males and females, many have sex-biased regulatory targeting patterns. In each tissue, genes that are differentially targeted by TFs between the sexes are enriched for tissue-related functions and diseases. In brain tissue, for example, genes associated with Parkinson’s disease and Alzheimer’s disease are targeted by different sets of TFs in each sex. Our systems-based analysis identifies a repertoire of TFs that play important roles in sex-specific architecture of gene regulatory networks, and it underlines sex-specific regulatory processes in both health and disease.

Epigenetic mechanisms integrate both genetic variability and environmental exposures. However, comprehensive epigenome-wide analysis has not been performed across major childhood allergic phenotypes. We examined the association of epigenome-wide DNA methylation in mid-childhood peripheral blood (Illumina HumanMethyl450K) with mid-childhood atopic sensitization, environmental/inhalant and food allergen sensitization in 739 children in two birth cohorts (Project Viva–Boston, and the Generation R Study–Rotterdam). We performed covariate-adjusted epigenome-wide association meta-analysis and employed pathway and regional analyses of results. Seven-hundred and five methylation sites (505 genes) were significantly cross-sectionally associated with mid-childhood atopic sensitization, 1411 (905 genes) for environmental and 45 (36 genes) for food allergen sensitization (FDR

Characterizing the collective regulatory impact of genetic variants on complex phenotypes is a major challenge in developing a genotype to phenotype map. Using expression quantitative trait locus (eQTL) analyses, we constructed bipartite networks in which edges represent significant associations between genetic variants and gene expression levels and found that the network structure informs regulatory function. We show, in 13 tissues, that these eQTL networks are organized into dense, highly modular communities grouping genes often involved in coherent biological processes. We find communities representing shared processes across tissues, as well as communities associated with tissue-specific processes that coalesce around variants in tissue-specific active chromatin regions. Node centrality is also highly informative, with the global and community hubs differing in regulatory potential and likelihood of being disease associated.

Summary Although all human tissues carry out common processes, tissues are distinguished by gene expression patterns, implying that distinct regulatory programs control tissue specificity. In this study, we investigate gene expression and regulation across 38 tissues profiled in the Genotype-Tissue Expression project. We find that network edges (transcription factor to target gene connections) have higher tissue specificity than network nodes (genes) and that regulating nodes (transcription factors) are less likely to be expressed in a tissue-specific manner as compared to their targets (genes). Gene set enrichment analysis of network targeting also indicates that the regulation of tissue-specific function is largely independent of transcription factor expression. In addition, tissue-specific genes are not highly targeted in their corresponding tissue network. However, they do assume bottleneck positions due to variability in transcription factor targeting and the influence of non-canonical regulatory interactions. These results suggest that tissue specificity is driven by context-dependent regulatory paths, providing transcriptional control of tissue-specific processes., Graphical abstract

SummaryExpression quantitative trait locus (eQTL) analysis associates genotype with gene expression, but most eQTL studies only includecis-acting variants and generally examine a single tissue. We used data from 13 tissues obtained by the Genotype-Tissue Expression (GTEx) project v6.0 and, in each tissue, identified bothcis- andtrans-eQTLs. For each tissue, we represented significant associations between single nucleotide polymorphisms (SNPs) and genes as edges in a bipartite network. These networks are organized into dense, highly modular communities often representing coherent biological processes. Global network hubs are enriched in distal gene regulatory regions such as enhancers, but are devoid of disease-associated SNPs from genome wide association studies. In contrast, local, community-specific network hubs (core SNPs) are preferentially located in regulatory regions such as promoters and enhancers and highly enriched for trait and disease associations. These results provide help explain how many weak-effect SNPs might together influence cellular function and phenotype.

SummaryCell lines are an indispensable tool in biomedical research and often used as surrogates for tissues. An important question is how well a cell line’s transcriptional and regulatory processes reflect those of its tissue of origin. We analyzed RNA-Seq data from GTEx for 127 paired Epstein-Barr virus transformed lymphoblastoid cell lines and whole blood samples; and 244 paired fibroblast cell lines and skin biopsies. A combination of gene expression and network analyses shows that while cell lines carry the expression signatures of their primary tissues, albeit at reduced levels, they also exhibit changes in their patterns of transcription factor regulation. Cell cycle genes are over-expressed in cell lines compared to primary tissue, and they have a reduction of repressive transcription factor targeting. Our results provide insight into the expression and regulatory alterations observed in cell lines and suggest that these changes should be considered when using cell lines as models.HighlightsCell lines differ from their source tissues in gene expression and regulationDistinct cell lines share altered patterns of cell cycle regulationCell cycle genes are less strongly targeted by repressive TFs in cell linesCell lines share expression with their source tissue, but at reduced levels

SummarySexual dimorphism manifests in many diseases and may drive sex-specific therapeutic responses. To understand the molecular basis of sexual dimorphism, we conducted a comprehensive assessment of gene expression and regulatory network modeling in 31 tissues using 8716 human transcriptomes from GTEx. We observed sexually dimorphic patterns of gene expression involving as many as 60% of autosomal genes, depending on the tissue. Interestingly, sex hormone receptors do not exhibit sexually dimorphic expression in most tissues; however, differential network targeting by hormone receptors and other transcription factors (TFs) captures their downstream sexually dimorphic gene expression. Furthermore, differential network wiring was found extensively in several tissues, particularly in brain, in which not all regions exhibit strong differential expression. This systems-based analysis provides a new perspective on the drivers of sexual dimorphism, one in which a repertoire of TFs plays important roles in sex-specific rewiring of gene regulatory networks.HighlightsSexual dimorphism manifests in both gene expression and gene regulatory networksSubstantial sexual dimorphism in regulatory networks was found in several tissuesMany differentially regulated genes are not differentially expressedSex hormone receptors do not exhibit sexually dimorphic expression in most tissues

Although ultrahigh-throughput RNA-Sequencing has become the dominant technology for genome-wide transcriptional profiling, the vast majority of RNA-Seq studies typically profile only tens of samples, and most analytical pipelines are optimized for these smaller studies. However, projects are generating ever-larger data sets comprising RNA-Seq data from hundreds or thousands of samples, often collected at multiple centers and from diverse tissues. These complex data sets present significant analytical challenges due to batch and tissue effects, but provide the opportunity to revisit the assumptions and methods that we use to preprocess, normalize, and filter RNA-Seq data – critical first steps for any subsequent analysis. We find analysis of large RNA-Seq data sets requires both careful quality control and that one account for sparsity due to the heterogeneity intrinsic in multi-group studies. An R package instantiating our method for large-scale RNA-Seq normalization and preprocessing, YARN, is available at bioconductor.org/packages/yarn.HighlightsOverview of assumptions used in preprocessing and normalizationPipeline for preprocessing, quality control, and normalization of large heterogeneous dataA Bioconductor package for the YARN pipeline and easy manipulation of count dataPreprocessed GTEx data set using the YARN pipeline available as a resource

Extreme events taking place on networks are not uncommon. We show that it is possible to manipulate the extreme events occurrence probabilities and its distribution over the nodes on scale-free networks by tuning the nodal capacity. This can be used to reduce the number of extreme events occurrences on a network. However monotonic nodal capacity enhancements, beyond a point, do not lead to any substantial reduction in the number of extreme events. We point out the practical implication of this result for network design in the context of reducing extreme events occurrences., 5 pages, 4 figures

Networks with bimodal degree distribution are most robust to targeted and random attacks. We present a model for constructing a network with bimodal degree distribution. The procedure adopted is to add nodes to the network with a probability p and delete the links between nodes with probability (1 − p). We introduce an additional constraint in the process through an immunity score, which controls the dynamics of the growth process based on the feedback value of the last few time steps. This results in bimodal nature for the degree distribution. We study the standard quantities which characterize the networks, like average path length and clustering coefficient in the context of our growth process and show that the resultant network is in the small world family. It is interesting to note that bimodality in degree distribution is an emergent phenomenon.

We study a lattice model where the coupling stochastically switches between repulsive (subtractive) and attractive (additive) at each site with probability p at every time instant. We observe that such a kind of coupling stabilizes the local fixed point of a circle map, with the resultant globally stable attractor providing a unique absorbing state. Interestingly, a continuous phase transition is observed from the absorbing state to spatiotemporal chaos via spatiotemporal intermittency for a range of values of p . It is interesting to note that the transition falls in class of directed percolation. Static and spreading exponents along with relevant scaling laws are found to be obeyed confirming the directed percolation universality class in spatiotemporal intermittency regime.

We study the spatiotemporal dynamics of a network of coupled chaotic maps modelling neuronal activity, under variation of coupling strength $ϵ$ and degree of randomness in coupling $p$. We find that at high coupling strengths $(ϵg{ϵ}_{\text{fixed}})$ the unstable saddle point solution of the local chaotic maps gets stabilized. The range of coupling where this spatiotemporal fixed point gains stability is unchanged in the presence of randomness in the connections, namely ${ϵ}_{\text{fixed}}$ is invariant under changes in $p$. As coupling gets weaker $(ϵl{ϵ}_{\text{fixed}})$, the spatiotemporal fixed point loses stability, and one obtains chaos. In this regime, when the coupling connections are completely regular $(p=0)$, the network becomes spatiotemporally chaotic. Interestingly however, in the presence of random links $(pg0)$ one obtains spatial synchronization in the network. We find that this range of synchronized chaos increases exponentially with the fraction of random links in the network. Further, in the space of fixed coupling strengths, the synchronization transition occurs at a finite value of $p$, a scenario quite distinct from the many examples of synchronization transitions at $p\ensuremath{\rightarrow}0$. Further we show that the synchronization here is robust in the presence of parametric noise, namely in a network of nonidentical neuronal maps. Finally we check the generality of our observations in networks of neurons displaying both spiking and bursting dynamics.

We investigate coupled circle maps in the presence of feedback and explore various dynamical phases observed in this system of coupled high dimensional maps. We observe an interesting transition from localized chaos to spatiotemporal chaos. We study this transition as a dynamic phase transition. We observe that persistence acts as an excellent quantifier to describe this transition. Taking the location of the fixed point of circle map (which does not change with feedback) as a reference point, we compute a number of sites which have been greater than (less than) the fixed point until time t. Though local dynamics is high dimensional in this case, this definition of persistence which tracks a single variable is an excellent quantifier for this transition. In most cases, we also obtain a well defined persistence exponent at the critical point and observe conventional scaling as seen in second order phase transitions. This indicates that persistence could work as a good order parameter for transitions from fully or partially arrested phase. We also give an explanation of gaps in eigenvalue spectrum of the Jacobian of localized state.