Isthmus progenitor cells contribute to homeostatic cellular turnover and support regeneration following intestinal injury.

The currently accepted intestinal epithelial cell organization model proposes that Lgr5+ crypt-base columnar (CBC) cells represent the sole intestinal stem cell (ISC) compartment. However, previous studies have indicated that Lgr5+ cells are dispensable for intestinal regeneration, leading to two major hypotheses: one favoring the presence of a quiescent reserve ISC and the other calling for differentiated cell plasticity. To investigate these possibilities, we studied crypt epithelial cells in an unbiased fashion via high-resolution single-cell profiling. These studies, combined with in vivo lineage tracing, show that Lgr5 is not a specific ISC marker and that stemness potential exists beyond the crypt base and resides in the isthmus region, where undifferentiated cells participate in intestinal homeostasis and regeneration following irradiation (IR) injury. Our results provide an alternative model of intestinal epithelial cell organization, suggesting that stemness potential is not restricted to CBC cells, and neither de-differentiation nor reserve ISC are drivers of intestinal regeneration.

A modular master regulator landscape controls cancer transcriptional identity.

Despite considerable efforts, the mechanisms linking genomic alterations to the transcriptional identity of cancer cells remain elusive. Integrative genomic analysis, using a network-based approach, identified 407 master regulator (MR) proteins responsible for canalizing the genetics of individual samples from 20 cohorts in The Cancer Genome Atlas (TCGA) into 112 transcriptionally distinct tumor subtypes. MR proteins could be further organized into 24 pan-cancer, master regulator block modules (MRBs), each regulating key cancer hallmarks and predictive of patient outcome in multiple cohorts. Of all somatic alterations detected in each individual sample, >50% were predicted to induce aberrant MR activity, yielding insight into mechanisms linking tumor genetics and transcriptional identity and establishing non-oncogene dependencies. Genetic and pharmacological validation assays confirmed the predicted effect of upstream mutations and MR activity on downstream cellular identity and phenotype. Thus, co-analysis of mutational and gene expression profiles identified elusive subtypes and provided testable hypothesis for mechanisms mediating the effect of genetic alterations.

Transcriptional Architecture of Synaptic Communication Delineates GABAergic Neuron Identity.

Understanding the organizational logic of neural circuits requires deciphering the biological basis of neuronal diversity and identity, but there is no consensus on how neuron types should be defined. We analyzed single-cell transcriptomes of a set of anatomically and physiologically characterized cortical GABAergic neurons and conducted a computational genomic screen for transcriptional profiles that distinguish them from one another. We discovered that cardinal GABAergic neuron types are delineated by a transcriptional architecture that encodes their synaptic communication patterns. This architecture comprises 6 categories of ∼40 gene families, including cell-adhesion molecules, transmitter-modulator receptors, ion channels, signaling proteins, neuropeptides and vesicular release components, and transcription factors. Combinatorial expression of select members across families shapes a multi-layered molecular scaffold along the cell membrane that may customize synaptic connectivity patterns and input-output signaling properties. This molecular genetic framework of neuronal identity integrates cell phenotypes along multiple axes and provides a foundation for discovering and classifying neuron types.

Comparative analysis of gene regulatory networks: from network reconstruction to evolution.

Regulation of gene expression is central to many biological processes. Although reconstruction of regulatory circuits from genomic data alone is therefore desirable, this remains a major computational challenge. Comparative approaches that examine the conservation and divergence of circuits and their components across strains and species can help reconstruct circuits as well as provide insights into the evolution of gene regulatory processes and their adaptive contribution. In recent years, advances in genomic and computational tools have led to a wealth of methods for such analysis at the sequence, expression, pathway, module, and entire network level. Here, we review computational methods developed to study transcriptional regulatory networks using comparative genomics, from sequence to functional data. We highlight how these methods use evolutionary conservation and divergence to reliably detect regulatory components as well as estimate the extent and rate of divergence. Finally, we discuss the promise and open challenges in linking regulatory divergence to phenotypic divergence and adaptation.

The Immune Landscape of Cancer.

We performed an extensive immunogenomic analysis of more than 10,000 tumors comprising 33 diverse cancer types by utilizing data compiled by TCGA. Across cancer types, we identified six immune subtypes-wound healing, IFN-γ dominant, inflammatory, lymphocyte depleted, immunologically quiet, and TGF-ß dominant-characterized by differences in macrophage or lymphocyte signatures, Th1:Th2 cell ratio, extent of intratumoral heterogeneity, aneuploidy, extent of neoantigen load, overall cell proliferation, expression of immunomodulatory genes, and prognosis. Specific driver mutations correlated with lower (CTNNB1, NRAS, or IDH1) or higher (BRAF, TP53, or CASP8) leukocyte levels across all cancers. Multiple control modalities of the intracellular and extracellular networks (transcription, microRNAs, copy number, and epigenetic processes) were involved in tumor-immune cell interactions, both across and within immune subtypes. Our immunogenomics pipeline to characterize these heterogeneous tumors and the resulting data are intended to serve as a resource for future targeted studies to further advance the field.

Bayesian functional analysis for untargeted metabolomics data with matching uncertainty and small sample sizes.

Untargeted metabolomics based on liquid chromatography-mass spectrometry technology is quickly gaining widespread application, given its ability to depict the global metabolic pattern in biological samples. However, the data are noisy and plagued by the lack of clear identity of data features measured from samples. Multiple potential matchings exist between data features and known metabolites, while the truth can only be one-to-one matches. Some existing methods attempt to reduce the matching uncertainty, but are far from being able to remove the uncertainty for most features. The existence of the uncertainty causes major difficulty in downstream functional analysis. To address these issues, we develop a novel approach for Bayesian Analysis of Untargeted Metabolomics data (BAUM) to integrate previously separate tasks into a single framework, including matching uncertainty inference, metabolite selection and functional analysis. By incorporating the knowledge graph between variables and using relatively simple assumptions, BAUM can analyze datasets with small sample sizes. By allowing different confidence levels of feature-metabolite matching, the method is applicable to datasets in which feature identities are partially known. Simulation studies demonstrate that, compared with other existing methods, BAUM achieves better accuracy in selecting important metabolites that tend to be functionally consistent and assigning confidence scores to feature-metabolite matches. We analyze a COVID-19 metabolomics dataset and a mouse brain metabolomics dataset using BAUM. Even with a very small sample size of 16 mice per group, BAUM is robust and stable. It finds pathways that conform to existing knowledge, as well as novel pathways that are biologically plausible.

IDMIR: identification of dysregulated miRNAs associated with disease based on a miRNA-miRNA interaction network constructed through gene expression data.

Micro ribonucleic acids (miRNAs) play a pivotal role in governing the human transcriptome in various biological phenomena. Hence, the accumulation of miRNA expression dysregulation frequently assumes a noteworthy role in the initiation and progression of complex diseases. However, accurate identification of dysregulated miRNAs still faces challenges at the current stage. Several bioinformatics tools have recently emerged for forecasting the associations between miRNAs and diseases. Nonetheless, the existing reference tools mainly identify the miRNA-disease associations in a general state and fall short of pinpointing dysregulated miRNAs within a specific disease state. Additionally, no studies adequately consider miRNA-miRNA interactions (MMIs) when analyzing the miRNA-disease associations. Here, we introduced a systematic approach, called IDMIR, which enabled the identification of expression dysregulated miRNAs through an MMI network under the gene expression context, where the network's architecture was designed to implicitly connect miRNAs based on their shared biological functions within a particular disease context. The advantage of IDMIR is that it uses gene expression data for the identification of dysregulated miRNAs by analyzing variations in MMIs. We illustrated the excellent predictive power for dysregulated miRNAs of the IDMIR approach through data analysis on breast cancer and bladder urothelial cancer. IDMIR could surpass several existing miRNA-disease association prediction approaches through comparison. We believe the approach complements the deficiencies in predicting miRNA-disease association and may provide new insights and possibilities for diagnosing and treating diseases. The IDMIR approach is now available as a free R package on CRAN (https://CRAN.R-project.org/package=IDMIR).

Veterinary systems biology for bridging the phenotype-genotype gap via computational modeling for disease epidemiology and animal welfare.

Veterinary systems biology is an innovative approach that integrates biological data at the molecular and cellular levels, allowing for a more extensive understanding of the interactions and functions of complex biological systems in livestock and veterinary science. It has tremendous potential to integrate multi-omics data with the support of vetinformatics resources for bridging the phenotype-genotype gap via computational modeling. To understand the dynamic behaviors of complex systems, computational models are frequently used. It facilitates a comprehensive understanding of how a host system defends itself against a pathogen attack or operates when the pathogen compromises the host's immune system. In this context, various approaches, such as systems immunology, network pharmacology, vaccinology and immunoinformatics, can be employed to effectively investigate vaccines and drugs. By utilizing this approach, we can ensure the health of livestock. This is beneficial not only for animal welfare but also for human health and environmental well-being. Therefore, the current review offers a detailed summary of systems biology advancements utilized in veterinary sciences, demonstrating the potential of the holistic approach in disease epidemiology, animal welfare and productivity.

Early Transcriptional Divergence Marks Virus-Specific Primary Human CD8+ T Cells in Chronic versus Acute Infection.

Distinct molecular pathways govern the differentiation of CD8+ effector T cells into memory or exhausted T cells during acute and chronic viral infection, but these are not well studied in humans. Here, we employed an integrative systems immunology approach to identify transcriptional commonalities and differences between virus-specific CD8+ T cells from patients with persistent and spontaneously resolving hepatitis C virus (HCV) infection during the acute phase. We observed dysregulation of metabolic processes during early persistent infection that was linked to changes in expression of genes related to nucleosomal regulation of transcription, T cell differentiation, and the inflammatory response and correlated with subject age, sex, and the presence of HCV-specific CD4+ T cell populations. These early changes in HCV-specific CD8+ T cell transcription preceded the overt establishment of T cell exhaustion, making this signature a prime target in the search for the regulatory origins of T cell dysfunction in chronic viral infection.

Comprehensive analysis of the proximity-dependent nuclear interactome for the oncoprotein NOTCH1 in live cells.

Notch signaling plays a critical role in cell fate decisions in all cell types. Furthermore, gain-of-function mutations in NOTCH1 have been uncovered in many human cancers. Disruption of Notch signaling has recently emerged as an attractive disease treatment strategy. However, the nuclear interaction landscape of the oncoprotein NOTCH1 remains largely unexplored. We therefore employed here a proximity-dependent biotin identification approach to identify in vivo protein associations with the nuclear Notch1 intracellular domain in live cells. We identified a large set of previously reported and unreported proteins that associate with NOTCH1, including general transcription and elongation factors, DNA repair and replication factors, coactivators, corepressors, and components of the NuRD and SWI/SNF chromatin remodeling complexes. We also found that Notch1 intracellular domain associates with protein modifiers and components of other signaling pathways that may influence Notch signal transduction and protein stability such as USP7. We further validated the interaction of NOTCH1 with histone deacetylase 1 or GATAD2B using protein network analysis, proximity-based ligation, in vivo cross-linking and coimmunoprecipitation assays in several Notch-addicted cancer cell lines. Through data mining, we also revealed potential drug targets for the inhibition of Notch signaling. Collectively, these results provide a valuable resource to uncover the mechanisms that fine-tune Notch signaling in tumorigenesis and inform therapeutic targets for Notch-addicted tumors.

How network-based approaches can complement gene identification studies in frontotemporal dementia.

Frontotemporal dementia (FTD) is a primary cause of dementia encompassing a broad range of clinical phenotypes and cellular pathologies. Genetic discoveries in FTD have largely been driven by linkage studies in well-documented extended families, explaining most of the patients with a known pathogenic mutation. In the context of complex diseases, it is hypothesized that mutations with reduced penetrance or a combination of low-effect size variants with environmental factors drive disease. Furthermore, these genes are likely to be part of the interaction networks of known FTD genes, contributing to converging cellular processes. In this review, we examine gene discovery approaches in FTD and introduce network biology concepts as tools to assist gene identification studies in genetically complex disease.

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.

Protein-protein interaction network-based integration of GWAS and functional data for blood pressure regulation analysis.

BACKGROUND: It is valuable to analyze the genome-wide association studies (GWAS) data for a complex disease phenotype in the context of the protein-protein interaction (PPI) network, as the related pathophysiology results from the function of interacting polyprotein pathways. The analysis may include the design and curation of a phenotype-specific GWAS meta-database incorporating genotypic and eQTL data linking to PPI and other biological datasets, and the development of systematic workflows for PPI network-based data integration toward protein and pathway prioritization. Here, we pursued this analysis for blood pressure (BP) regulation. METHODS: The relational scheme of the implemented in Microsoft SQL Server BP-GWAS meta-database enabled the combined storage of: GWAS data and attributes mined from GWAS Catalog and the literature, Ensembl-defined SNP-transcript associations, and GTEx eQTL data. The BP-protein interactome was reconstructed from the PICKLE PPI meta-database, extending the GWAS-deduced network with the shortest paths connecting all GWAS-proteins into one component. The shortest-path intermediates were considered as BP-related. For protein prioritization, we combined a new integrated GWAS-based scoring scheme with two network-based criteria: one considering the protein role in the reconstructed by shortest-path (RbSP) interactome and one novel promoting the common neighbors of GWAS-prioritized proteins. Prioritized proteins were ranked by the number of satisfied criteria. RESULTS: The meta-database includes 6687 variants linked with 1167 BP-associated protein-coding genes. The GWAS-deduced PPI network includes 1065 proteins, with 672 forming a connected component. The RbSP interactome contains 1443 additional, network-deduced proteins and indicated that essentially all BP-GWAS proteins are at most second neighbors. The prioritized BP-protein set was derived from the union of the most BP-significant by any of the GWAS-based or the network-based criteria. It included 335 proteins, with ~ 2/3 deduced from the BP PPI network extension and 126 prioritized by at least two criteria. ESR1 was the only protein satisfying all three criteria, followed in the top-10 by INSR, PTN11, CDK6, CSK, NOS3, SH2B3, ATP2B1, FES and FINC, satisfying two. Pathway analysis of the RbSP interactome revealed numerous bioprocesses, which are indeed functionally supported as BP-associated, extending our understanding about BP regulation. CONCLUSIONS: The implemented workflow could be used for other multifactorial diseases.

The dynamic causality brain network reflects whether the working memory is solidified.

Working memory, which is foundational to higher cognitive function, is the "sketchpad of volitional control." Successful working memory is the inevitable outcome of the individual's active control and manipulation of thoughts and turning them into internal goals during which the causal brain processes information in real time. However, little is known about the dynamic causality among distributed brain regions behind thought control that underpins successful working memory. In our present study, given that correct responses and incorrect ones did not differ in either contralateral delay activity or alpha suppression, further rooting on the high-temporal-resolution EEG time-varying directed network analysis, we revealed that successful working memory depended on both much stronger top-down connections from the frontal to the temporal lobe and bottom-up linkages from the occipital to the temporal lobe, during the early maintenance period, as well as top-down flows from the frontal lobe to the central areas as the delay behavior approached. Additionally, the correlation between behavioral performance and casual interactions increased over time, especially as memory-guided delayed behavior approached. Notably, when using the network metrics as features, time-resolved multiple linear regression of overall behavioral accuracy was exactly achieved as delayed behavior approached. These results indicate that accurate memory depends on dynamic switching of causal network connections and shifting to more task-related patterns during which the appropriate intervention may help enhance memory.

Uncovering Protein Networks in Cardiovascular Proteomics.

Biological networks have been widely used in many different diseases to identify potential biomarkers and design drug targets. In the present review, we describe the main computational techniques for reconstructing and analyzing different types of protein networks and summarize the previous applications of such techniques in cardiovascular diseases. Existing tools are critically compared, discussing when each method is preferred such as the use of co-expression networks for functional annotation of protein clusters and the use of directed networks for inferring regulatory associations. Finally, we are presenting examples of reconstructing protein networks of different types (regulatory, co-expression, and protein-protein interaction networks). We demonstrate the necessity to reconstruct networks separately for each cardiovascular tissue type and disease entity and provide illustrative examples of the importance of taking into consideration relevant post-translational modifications. Finally, we demonstrate and discuss how the findings of protein networks could be interpreted using single-cell RNA-sequencing data.

Calculation of centralities in protein kinase A.

Topological analysis of protein residue networks (PRNs) is a common method that can help to understand the roles of individual residues. Here, we used protein kinase A as a study object and asked what already known functionally important residues can be detected by network analysis. Along several traditional approaches to weight edges in PRNs we used local spatial pattern (LSP) alignment that assigns high weights to edges only if CαCß vectors for the corresponding residues retain their mutual positions and orientation. Our results show that even short molecular dynamic simulations of 10 to 20 ns can give convergent values for betweenness and degree centralities calculated from the LSP-based PRNs. Using these centralities, we were able to clearly distinguish a group of residues that are highly conserved in protein kinases and play important functional and regulatory roles. In comparison, traditional methods based on cross-correlation and linear mutual information were much less efficient for this particular task. These results call for reevaluation of the current methods to generate PRNs.

Disentangling direct from indirect relationships in association networks.

Networks are vital tools for understanding and modeling interactions in complex systems in science and engineering, and direct and indirect interactions are pervasive in all types of networks. However, quantitatively disentangling direct and indirect relationships in networks remains a formidable task. Here, we present a framework, called iDIRECT (Inference of Direct and Indirect Relationships with Effective Copula-based Transitivity), for quantitatively inferring direct dependencies in association networks. Using copula-based transitivity, iDIRECT eliminates/ameliorates several challenging mathematical problems, including ill-conditioning, self-looping, and interaction strength overflow. With simulation data as benchmark examples, iDIRECT showed high prediction accuracies. Application of iDIRECT to reconstruct gene regulatory networks in Escherichia coli also revealed considerably higher prediction power than the best-performing approaches in the DREAM5 (Dialogue on Reverse Engineering Assessment and Methods project, #5) Network Inference Challenge. In addition, applying iDIRECT to highly diverse grassland soil microbial communities in response to climate warming showed that the iDIRECT-processed networks were significantly different from the original networks, with considerably fewer nodes, links, and connectivity, but higher relative modularity. Further analysis revealed that the iDIRECT-processed network was more complex under warming than the control and more robust to both random and target species removal (P < 0.001). As a general approach, iDIRECT has great advantages for network inference, and it should be widely applicable to infer direct relationships in association networks across diverse disciplines in science and engineering.

The generality of cryptic dietary niche differences in diverse large-herbivore assemblages.

Ecological niche differences are necessary for stable species coexistence but are often difficult to discern. Models of dietary niche differentiation in large mammalian herbivores invoke the quality, quantity, and spatiotemporal distribution of plant tissues and growth forms but are agnostic toward food plant species identity. Empirical support for these models is variable, suggesting that additional mechanisms of resource partitioning may be important in sustaining large-herbivore diversity in African savannas. We used DNA metabarcoding to conduct a taxonomically explicit analysis of large-herbivore diets across southeastern Africa, analyzing ∼4,000 fecal samples of 30 species from 10 sites in seven countries over 6 y. We detected 893 food plant taxa from 124 families, but just two families-grasses and legumes-accounted for the majority of herbivore diets. Nonetheless, herbivore species almost invariably partitioned food plant taxa; diet composition differed significantly in 97% of pairwise comparisons between sympatric species, and dissimilarity was pronounced even between the strictest grazers (grass eaters), strictest browsers (nongrass eaters), and closest relatives at each site. Niche differentiation was weakest in an ecosystem recovering from catastrophic defaunation, indicating that food plant partitioning is driven by species interactions, and was stronger at low rainfall, as expected if interspecific competition is a predominant driver. Diets differed more between browsers than grazers, which predictably shaped community organization: Grazer-dominated trophic networks had higher nestedness and lower modularity. That dietary differentiation is structured along taxonomic lines complements prior work on how herbivores partition plant parts and patches and suggests that common mechanisms govern herbivore coexistence and community assembly in savannas.

Comparative proteomic analysis across the developmental stages of the Eimeria tenella.

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

Molecular pathways underlying lung-brain axis signaling in asthma: Relevance for psychopathology and neuroinflammation.

BACKGROUND: Accumulating evidence indicates that asthma has systemic effects and affects brain function. Although airway inflammation is proposed to initiate afferent communications with the brain, the signaling pathways have not been established. OBJECTIVE: We sought to identify the cellular and molecular pathways involved in afferent lung-brain communication during airway inflammation in asthma. METHODS: In 23 adults with mild asthma, segmental bronchial provocation with allergen (SBP-Ag) was used to provoke airway inflammation and retrieve bronchoalveolar lavage fluid for targeted protein analysis and RNA sequencing to determine gene expression profiles. Neural responses to emotional cues in nodes of the salience network were assessed with functional magnetic resonance imaging at baseline and 48 hours after SBP-Ag. RESULTS: Cell deconvolution and gene coexpression network analysis identified 11 cell-associated gene modules that changed in response to SBP-Ag. SBP-Ag increased bronchoalveolar lavage eosinophils and expression of an eosinophil-associated module enriched for genes related to TH17-type inflammation (eg, IL17A), as well as cell proliferation in lung and brain (eg, NOTCH1, VEGFA, and LIF). Increased expression of genes in this module, as well as several TH17-type inflammation-related proteins, was associated with an increase from baseline in salience network reactivity. CONCLUSIONS: Our results identify a specific inflammatory pathway linking asthma-related airway inflammation and emotion-related neural function. Systemically, TH17-type inflammation has been implicated in both depression and neuroinflammation, with impacts on long-term brain health. Thus, our data emphasize that inflammation in the lung in asthma may have profound effects outside of the lung that may be targetable with novel therapeutic approaches.