Transkingdom Network Analysis (TkNA): a systems framework for inferring causal factors underlying host-microbiota and other multi-omic interactions.

We present Transkingdom Network Analysis (TkNA), a unique causal-inference analytical framework that offers a holistic view of biological systems by integrating data from multiple cohorts and diverse omics types. TkNA helps to decipher key players and mechanisms governing host-microbiota (or any multi-omic data) interactions in specific conditions or diseases. TkNA reconstructs a network that represents a statistical model capturing the complex relationships between different omics in the biological system. It identifies robust and reproducible patterns of fold change direction and correlation sign across several cohorts to select differential features and their per-group correlations. The framework then uses causality-sensitive metrics, statistical thresholds and topological criteria to determine the final edges forming the transkingdom network. With the subsequent network's topological features, TkNA identifies nodes controlling a given subnetwork or governing communication between kingdoms and/or subnetworks. The computational time for the millions of correlations necessary for network reconstruction in TkNA typically takes only a few minutes, varying with the study design. Unlike most other multi-omics approaches that find only associations, TkNA focuses on establishing causality while accounting for the complex structure of multi-omic data. It achieves this without requiring huge sample sizes. Moreover, the TkNA protocol is user friendly, requiring minimal installation and basic familiarity with Unix. Researchers can access the TkNA software at https://github.com/CAnBioNet/TkNA/ .

Intestinal IL-22RA1 signaling regulates intrinsic and systemic lipid and glucose metabolism to alleviate obesity-associated disorders.

IL-22 is critical for ameliorating obesity-induced metabolic disorders. However, it is unknown where IL-22 acts to mediate these outcomes. Here we examine the importance of tissue-specific IL-22RA1 signaling in mediating long-term high fat diet (HFD) driven metabolic disorders. To do so, we generated intestinal epithelium-, liver-, and white adipose tissue (WAT)-specific Il22ra1 knockout and littermate control mice. Intestinal epithelium- and liver-specific IL-22RA1 signaling upregulated systemic glucose metabolism. Intestinal IL-22RA1 signaling also mediated liver and WAT metabolism in a microbiota-dependent manner. We identified an association between Oscillibacter and elevated WAT inflammation, likely induced by Mmp12 expressing macrophages. Mechanistically, transcription of intestinal lipid metabolism genes is regulated by IL-22 and potentially IL-22-induced IL-18. Lastly, we show that Paneth cell-specific IL-22RA1 signaling, in part, mediates systemic glucose metabolism after HFD. Overall, these results elucidate a key role of intestinal epithelium-specific IL-22RA1 signaling in regulating intestinal metabolism and alleviating systemic obesity-associated disorders.

Multi-omic network analysis identified betacellulin as a novel target of omega-3 fatty acid attenuation of western diet-induced nonalcoholic steatohepatitis.

Clinical and preclinical studies established that supplementing diets with ω3 polyunsaturated fatty acids (PUFA) can reduce hepatic dysfunction in nonalcoholic steatohepatitis (NASH) but molecular underpinnings of this action were elusive. Herein, we used multi-omic network analysis that unveiled critical molecular pathways involved in ω3 PUFA effects in a preclinical mouse model of western diet induced NASH. Since NASH is a precursor of liver cancer, we also performed meta-analysis of human liver cancer transcriptomes that uncovered betacellulin as a key EGFR-binding protein upregulated in liver cancer and downregulated by ω3 PUFAs in animals and humans with NASH. We then confirmed that betacellulin acts by promoting proliferation of quiescent hepatic stellate cells, inducing transforming growth factor-ß2 and increasing collagen production. When used in combination with TLR2/4 agonists, betacellulin upregulated integrins in macrophages thereby potentiating inflammation and fibrosis. Taken together, our results suggest that suppression of betacellulin is one of the key mechanisms associated with anti-inflammatory and anti-fibrotic effects of ω3 PUFA on NASH.

Early transcriptome changes associated with western diet induced NASH in Ldlr-/- mice points to activation of hepatic macrophages and an acute phase response.

Background: Nonalcoholic fatty liver disease (NAFLD) is a global health problem. Identifying early gene indicators contributing to the onset and progression of NAFLD has the potential to develop novel targets for early therapeutic intervention. We report on the early and late transcriptomic signatures of western diet (WD)-induced nonalcoholic steatohepatitis (NASH) in female and male Ldlr-/- mice, with time-points at 1 week and 40 weeks on the WD. Control Ldlr-/- mice were maintained on a low-fat diet (LFD) for 1 and 40 weeks. Methods: The approach included quantitation of anthropometric and hepatic histology markers of disease as well as the hepatic transcriptome. Results: Only mice fed the WD for 40 weeks revealed evidence of NASH, i.e., hepatic steatosis and fibrosis. RNASeq transcriptome analysis, however, revealed multiple cell-specific changes in gene expression after 1 week that persisted to 40 weeks on the WD. These early markers of disease include induction of acute phase response (Saa1-2, Orm2), fibrosis (Col1A1, Col1A2, TGFß) and NASH associated macrophage (NAM, i.e., Trem2 high, Mmp12 low). We also noted the induction of transcripts associated with metabolic syndrome, including Mmp12, Trem2, Gpnmb, Lgals3 and Lpl. Finally, 1 week of WD feeding was sufficient to significantly induce TNFα, a cytokine involved in both hepatic and systemic inflammation. Conclusion: This study revealed early onset changes in the hepatic transcriptome that develop well before any anthropometric or histological evidence of NALFD or NASH and pointed to cell-specific targeting for the prevention of disease progression.

Transkingdom Network Analysis (TkNA): a systems framework for inferring causal factors underlying host-microbiota and other multi-omic interactions.

Technological advances have generated tremendous amounts of high-throughput omics data. Integrating data from multiple cohorts and diverse omics types from new and previously published studies can offer a holistic view of a biological system and aid in deciphering its critical players and key mechanisms. In this protocol, we describe how to use Transkingdom Network Analysis (TkNA), a unique causal-inference analytical framework that can perform meta-analysis of cohorts and detect master regulators among measured parameters that govern pathological or physiological responses of host-microbiota (or any multi-omic data) interactions in a particular condition or disease. TkNA first reconstructs the network that represents a statistical model capturing the complex relationships between the different omics of the biological system. Here, it selects differential features and their per-group correlations by identifying robust and reproducible patterns of fold change direction and sign of correlation across several cohorts. Next, a causality-sensitive metric, statistical thresholds, and a set of topological criteria are used to select the final edges that form the transkingdom network. The second part of the analysis involves interrogating the network. Using the network's local and global topology metrics, it detects nodes that are responsible for control of given subnetwork or control of communication between kingdoms and/or subnetworks. The underlying basis of the TkNA approach involves fundamental principles including laws of causality, graph theory and information theory. Hence, TkNA can be used for causal inference via network analysis of any host and/or microbiota multi-omics data. This quick and easy-to-run protocol requires very basic familiarity with the Unix command-line environment.

Microbiota and adipocyte mitochondrial damage in type 2 diabetes are linked by Mmp12+ macrophages.

Microbiota contribute to the induction of type 2 diabetes by high-fat/high-sugar (HFHS) diet, but which organs/pathways are impacted by microbiota remain unknown. Using multiorgan network and transkingdom analyses, we found that microbiota-dependent impairment of OXPHOS/mitochondria in white adipose tissue (WAT) plays a primary role in regulating systemic glucose metabolism. The follow-up analysis established that Mmp12+ macrophages link microbiota-dependent inflammation and OXPHOS damage in WAT. Moreover, the molecular signature of Mmp12+ macrophages in WAT was associated with insulin resistance in obese patients. Next, we tested the functional effects of MMP12 and found that Mmp12 genetic deficiency or MMP12 inhibition improved glucose metabolism in conventional, but not in germ-free mice. MMP12 treatment induced insulin resistance in adipocytes. TLR2-ligands present in Oscillibacter valericigenes bacteria, which are expanded by HFHS, induce Mmp12 in WAT macrophages in a MYD88-ATF3-dependent manner. Thus, HFHS induces Mmp12+ macrophages and MMP12, representing a microbiota-dependent bridge between inflammation and mitochondrial damage in WAT and causing insulin resistance.

Identifying and Tackling Emergent Vulnerability in Drug-Resistant Mycobacteria.

The global mechanisms and associated molecular alterations that occur in drug-resistant mycobacteria are poorly understood. To address this, we obtain genomics data and then construct a genome-scale response network in isoniazid-resistant Mycobacterium smegmatis and apply a network-mining algorithm. Through this, we decipher global alterations in an unbiased manner and identify emergent vulnerabilities in resistant bacilli, of which redox response was prominent. Using phenotypic profiling, we find that resistant bacilli exhibit collateral sensitivity to several compounds that block antioxidant responses. We find that nanogram/milliliter concentrations of ebselen, vancomycin, and phenylarsine oxide, in combination with isoniazid, are highly effective against Mycobacterium tuberculosis H37Rv and three clinical drug-resistant strains. Dynamic measurements of cytoplasmic redox potential revealed a surprisingly diminished capacity of clinical drug-resistant strains to counteract oxidative stress, providing a mechanistic basis for efficient and synergistic mycobactericidal activity of the drug combinations. Ebselen and vancomycin appear to be promising repurposable drugs.

Probing the Druggability Limits for Enzymes of the NAD Biosynthetic Network in Glioma.

The biosynthesis of NAD constitutes an important metabolic module in the cell, since NAD is an essential cofactor involved in several metabolic reactions. NAD concentrations are known to be significantly increased in several cancers, particularly in glioma, consistent with the observation of up-regulation of several enzymes of the network. Modulating NAD biosynthesis in glioma is therefore an attractive therapeutic strategy. Here we report reconstruction of a biochemical network of NAD biosynthesis consisting of 22 proteins, 36 metabolites, and 86 parameters, tuned to mimic the conditions in glioma. Kinetic simulations of the network provide comprehensive insights about the role of individual enzymes. Further, quantitative changes in the same network between different states of health and disease enable identification of drug targets, based on specific alterations in the given disease. Through simulations of enzyme inhibition titrations, we identify NMPRTase as a potential drug target, while eliminating other possible candidates NMNAT, NAPRTase, and NRK. We have also simulated titrations of both binding affinities as well as inhibitor concentrations, which provide insights into the druggability limits of the target, a novel aspect that can provide useful guidelines for designing inhibitors with optimal affinities. Our simulations suggest that an inhibitor affinity of 10 nM used in a concentration range of 0.1 to 10 µM achieves a near maximal inhibition response for NMPRTase and that increasing the affinity any further is not likely to have a significant advantage. Thus, the quantitative appreciation defines a maximal extent of inhibition possible for a chosen enzyme in the context of its network. Knowledge of this type enables an upper affinity threshold to be defined as a goal in lead screening and refinement stages in drug discovery.

Gene expression profiles of wild-type and isoniazid-resistant strains of Mycobacterium smegmatis.

The global variations in the gene expression pattern of drug treated (½X) and laboratory evolved drug-resistant strains (2XR and 4XR) of Mycobacterium smegmatis were obtained and compared with the M. smegmatis mc(2) 155 (WT) strain. The genes exhibiting two-fold change and p-value <0.05 under the treated conditions have been considered as differentially expressed genes (DEGs). Overall, the numbers of DEGs observed are 1529 in ½X (596-up, 933-down), 1381 (899-up, 482-down) in 2XR and 716 in 4XR (267-up, 449-down) conditions. The data is publicly available through the GEO database with accession number GSE64132.

Complete Genome Sequences of a Mycobacterium smegmatis Laboratory Strain (MC2 155) and Isoniazid-Resistant (4XR1/R2) Mutant Strains.

We report the whole genome sequences of a Mycobacterium smegmatis laboratory wild-type strain (MC(2) 155) and mutants (4XR1, 4XR2) resistant to isoniazid. Compared to Mycobacterium smegmatis MC(2) 155 (NC_008596), a widely used strain in laboratory experiments, the MC(2) 155, 4XR1, and 4XR2 strains are 60, 128 and 93 bp longer, respectively.

Identifying feasible metabolic routes in Mycobacterium smegmatis and possible alterations under diverse nutrient conditions.

BACKGROUND: Many studies on M. tuberculosis have emerged from using M. smegmatis MC(2)155 (Msm), since they share significant similarities and yet Msm is non-pathogenic and faster growing. Although several individual molecules have been studied from Msm, many questions remain open about its metabolism as a whole and its capability to be versatile. Adaptability and versatility are emergent properties of a system, warranting a molecular systems perspective to understand them. RESULTS: We identify feasible metabolic pathways in Msm in reference condition with transcriptome, phenotypic microarray, along with functional annotation of the genome. Together with transcriptome data, specific genes from a set of alternatives have been mapped onto different pathways. About 257 metabolic pathways can be considered to be feasible in Msm. Next, we probe cellular metabolism with an array of alternative carbon and nitrogen sources and identify those that are utilized and favour growth as well as those that do not support growth. In all, about 135 points in the entire metabolic map are probed. Analyzing growth patterns under these conditions, lead us to hypothesize different pathways that can become active in various conditions and possible alternate routes that may be induced, thus explaining the observed physiological adaptations. CONCLUSIONS: The study provides the first detailed analysis of feasible pathways towards adaptability. We obtain mechanistic insights that explain observed phenotypic behaviour by studying gene-expression profiles and pathways inferred from the genome sequence. Comparison of transcriptome and phenome analysis of Msm and Mtb provides a rationale for understanding commonalities in metabolic adaptability.

Rationalization and prediction of drug resistant mutations in targets for clinical anti-tubercular drugs.

Resistance to therapy limits the effectiveness of drug treatment in many diseases. Drug resistance can be considered as a successful outcome of the bacterial struggle to survive in the hostile environment of a drug-exposed cell. An important mechanism by which bacteria acquire drug resistance is through mutations in the drug target. Drug resistant strains (multi-drug resistant and extensively drug resistant) of Mycobacterium tuberculosis are being identified at alarming rates, increasing the global burden of tuberculosis. An understanding of the nature of mutations in different drug targets and how they achieve resistance is therefore important. An objective of this study is to first decipher sequence as well as structural bases for the observed resistance in known drug resistant mutants and then to predict positions in each target that are more prone to acquiring drug resistant mutations. A curated database containing hundreds of mutations in the 38 drug targets of nine major clinical drugs, associated with resistance is studied here. Mutations have been classified into those that occur in the binding site itself, those that occur in residues interacting with the binding site and those that occur in outer zones. Structural models of the wild type and mutant forms of the target proteins have been analysed to seek explanations for reduction in drug binding. Stability analysis of an entire array of 19 mutations at each of the residues for each target has been computed using structural models. Conservation indices of individual residues, binding sites and whole proteins are computed based on sequence conservation analysis of the target proteins. The analyses lead to insights about which positions in the polypeptide chain have a higher propensity to acquire drug resistant mutations. Thus critical insights can be obtained about the effect of mutations on drug binding, in terms of which amino acid positions and therefore which interactions should not be heavily relied upon, which in turn can be translated into guidelines for modifying the existing drugs as well as for designing new drugs. The methodology can serve as a general framework to study drug resistant mutants in other micro-organisms as well.

Network approaches to drug discovery.

INTRODUCTION: Advances in genomics technologies are providing a very large amount of data on genome-wide gene expression profiles, protein molecules and their interactions with other macromolecules and metabolites. Molecular interaction networks provide a useful way to capture this complex data and comprehend it. Networks are beginning to be used in drug discovery, in many steps of the modern discovery pipeline, with large-scale molecular networks being particularly useful for the understanding of the molecular basis of the disease. AREAS COVERED: The authors discuss network approaches used for drug target discovery and lead identification in the drug discovery pipeline. By reconstructing networks of targets, drugs and drug candidates as well as gene expression profiles under normal and disease conditions, the paper illustrates how it is possible to find relationships between different diseases, find biomarkers, explore drug repurposing and study emergence of drug resistance. Furthermore, the authors also look at networks which address particular important aspects such as off-target effects, combination-targets, mechanism of drug action and drug safety. EXPERT OPINION: The network approach represents another paradigm shift in drug discovery science. A network approach provides a fresh perspective of understanding important proteins in the context of their cellular environments, providing a rational basis for deriving useful strategies in drug design. Besides drug target identification and inferring mechanism of action, networks enable us to address new ideas that could prove to be extremely useful for new drug discovery, such as drug repositioning, drug synergy, polypharmacology and personalized medicine.

Protein-protein interaction networks suggest different targets have different propensities for triggering drug resistance.

Emergence of drug resistance is a major problem in the treatment of many diseases including tuberculosis. To tackle the problem from a wholistic perspective, it is essential to understand the molecular mechanisms by which bacteria acquire drug resistance using a systems approach. Availability of genome-scale data of expression profiles under different drug exposed conditions and protein-protein interactions, makes it feasible to reconstruct and analyze systems-level models. A number of proteins involved in different resistance mechanisms, referred to as the resistome are identified from literature. The interaction of the drug directly with the resistome is unable to explain most resistance processes adequately, including that of increased mutations in the target's binding site. We recently hypothesized that some communication might exist from the drug environment to the resistome to trigger emergence of drug resistance. We report here a network based approach to identify most plausible paths of such communication in Mycobacterium tuberculosis. Networks capturing both structural and functional linkages among various proteins were weighted based on gene expression profiles upon exposure to specific drugs and betweenness centrality of the interactions. Our analysis suggests that different drug targets and hence different drugs could trigger the resistome to different extents and through different routes. The identified paths correlate well with the mechanisms known through experiment. Some examples of the top ranked hubs in multiple drug specific networks are PolA, FadD1, CydA, a monoxygenase and GltS, which could serve as co-targets, that could be inhibited in order to retard resistance related communication in the cell.