Comprehensive Multi-omics Analysis Reveals Mitochondrial Stress as a Central Biological Hub for Spaceflight Impact.

Spaceflight is known to impose changes on human physiology with unknown molecular etiologies. To reveal these causes, we used a multi-omics, systems biology analytical approach using biomedical profiles from fifty-nine astronauts and data from NASA's GeneLab derived from hundreds of samples flown in space to determine transcriptomic, proteomic, metabolomic, and epigenetic responses to spaceflight. Overall pathway analyses on the multi-omics datasets showed significant enrichment for mitochondrial processes, as well as innate immunity, chronic inflammation, cell cycle, circadian rhythm, and olfactory functions. Importantly, NASA's Twin Study provided a platform to confirm several of our principal findings. Evidence of altered mitochondrial function and DNA damage was also found in the urine and blood metabolic data compiled from the astronaut cohort and NASA Twin Study data, indicating mitochondrial stress as a consistent phenotype of spaceflight.

A Multiscale Map of the Stem Cell State in Pancreatic Adenocarcinoma.

Drug resistance and relapse remain key challenges in pancreatic cancer. Here, we have used RNA sequencing (RNA-seq), chromatin immunoprecipitation (ChIP)-seq, and genome-wide CRISPR analysis to map the molecular dependencies of pancreatic cancer stem cells, highly therapy-resistant cells that preferentially drive tumorigenesis and progression. This integrated genomic approach revealed an unexpected utilization of immuno-regulatory signals by pancreatic cancer epithelial cells. In particular, the nuclear hormone receptor retinoic-acid-receptor-related orphan receptor gamma (RORγ), known to drive inflammation and T cell differentiation, was upregulated during pancreatic cancer progression, and its genetic or pharmacologic inhibition led to a striking defect in pancreatic cancer growth and a marked improvement in survival. Further, a large-scale retrospective analysis in patients revealed that RORγ expression may predict pancreatic cancer aggressiveness, as it positively correlated with advanced disease and metastasis. Collectively, these data identify an orthogonal co-option of immuno-regulatory signals by pancreatic cancer stem cells, suggesting that autoimmune drugs should be evaluated as novel treatment strategies for pancreatic cancer patients.

Functional landscape of SARS-CoV-2 cellular restriction.

A deficient interferon (IFN) response to severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) infection has been implicated as a determinant of severe coronavirus disease 2019 (COVID-19). To identify the molecular effectors that govern IFN control of SARS-CoV-2 infection, we conducted a large-scale gain-of-function analysis that evaluated the impact of human IFN-stimulated genes (ISGs) on viral replication. A limited subset of ISGs were found to control viral infection, including endosomal factors inhibiting viral entry, RNA binding proteins suppressing viral RNA synthesis, and a highly enriched cluster of endoplasmic reticulum (ER)/Golgi-resident ISGs inhibiting viral assembly/egress. These included broad-acting antiviral ISGs and eight ISGs that specifically inhibited SARS-CoV-2 and SARS-CoV-1 replication. Among the broad-acting ISGs was BST2/tetherin, which impeded viral release and is antagonized by SARS-CoV-2 Orf7a protein. Overall, these data illuminate a set of ISGs that underlie innate immune control of SARS-CoV-2/SARS-CoV-1 infection, which will facilitate the understanding of host determinants that impact disease severity and offer potential therapeutic strategies for COVID-19.

Collagenolysis-dependent DDR1 signalling dictates pancreatic cancer outcome.

Pancreatic ductal adenocarcinoma (PDAC) is a highly desmoplastic, aggressive cancer that frequently progresses and spreads by metastasis to the liver1. Cancer-associated fibroblasts, the extracellular matrix and type I collagen (Col I) support2,3 or restrain the progression of PDAC and may impede blood supply and nutrient availability4. The dichotomous role of the stroma in PDAC, and the mechanisms through which it influences patient survival and enables desmoplastic cancers to escape nutrient limitation, remain poorly understood. Here we show that matrix-metalloprotease-cleaved Col I (cCol I) and intact Col I (iCol I) exert opposing effects on PDAC bioenergetics, macropinocytosis, tumour growth and metastasis. Whereas cCol I activates discoidin domain receptor 1 (DDR1)-NF-κB-p62-NRF2 signalling to promote the growth of PDAC, iCol I triggers the degradation of DDR1 and restrains the growth of PDAC. Patients whose tumours are enriched for iCol I and express low levels of DDR1 and NRF2 have improved median survival compared to those whose tumours have high levels of cCol I, DDR1 and NRF2. Inhibition of the DDR1-stimulated expression of NF-κB or mitochondrial biogenesis blocks tumorigenesis in wild-type mice, but not in mice that express MMP-resistant Col I. The diverse effects of the tumour stroma on the growth and metastasis of PDAC and on the survival of patients are mediated through the Col I-DDR1-NF-κB-NRF2 mitochondrial biogenesis pathway, and targeting components of this pathway could provide therapeutic opportunities.

Role of Hepatic Stellate and Liver Sinusoidal Endothelial Cells in a Human Primary Cell Three-Dimensional Model of Nonalcoholic Steatohepatitis.

Nonalcoholic steatohepatitis (NASH) is an inflammatory and fibrotic liver disease that has reached epidemic proportions and has no approved pharmacologic therapies. Research and drug development efforts are hampered by inadequate preclinical models. This research describes a three-dimensional bioprinted liver tissue model of NASH built using primary human hepatocytes and nonparenchymal liver cells (hepatic stellate cells, liver sinusoidal endothelial cells, and Kupffer cells) from either healthy or NASH donors. Three-dimensional tissues bioprinted with cells sourced from diseased patients showed a NASH phenotype, including fibrosis. More importantly, this NASH phenotype occurred without the addition of disease-inducing agents. Bioprinted tissues composed entirely of healthy cells exhibited significantly less evidence of disease. The role of individual cell types in driving the NASH phenotype was examined by producing chimeric bioprinted tissues composed of healthy cells together with the addition of one or more diseased nonparenchymal cell types. These experiments reveal a role for both hepatic stellate and liver sinusoidal endothelial cells in the disease process. This model represents a fully human system with potential to detect clinically active targets and eventually therapies.

Axon-derived PACSIN1 binds to the Schwann cell survival receptor, LRP1, and transactivates TrkC to promote gliatrophic activities.

Schwann cells (SCs) undergo phenotypic transformation and then orchestrate nerve repair following PNS injury. The ligands and receptors that activate and sustain SC transformation remain incompletely understood. Proteins released by injured axons represent important candidates for activating the SC Repair Program. The low-density lipoprotein receptor-related protein-1 (LRP1) is acutely up-regulated in SCs in response to injury, activating c-Jun, and promoting SC survival. To identify novel LRP1 ligands released in PNS injury, we applied a discovery-based approach in which extracellular proteins in the injured nerve were captured using Fc-fusion proteins containing the ligand-binding motifs of LRP1 (CCR2 and CCR4). An intracellular neuron-specific protein, Protein Kinase C and Casein Kinase Substrate in Neurons (PACSIN1) was identified and validated as an LRP1 ligand. Recombinant PACSIN1 activated c-Jun and ERK1/2 in cultured SCs. Silencing Lrp1 or inhibiting the LRP1 cell-signaling co-receptor, the NMDA-R, blocked the effects of PACSIN1 on c-Jun and ERK1/2 phosphorylation. Intraneural injection of PACSIN1 into crush-injured sciatic nerves activated c-Jun in wild-type mice, but not in mice in which Lrp1 is conditionally deleted in SCs. Transcriptome profiling of SCs revealed that PACSIN1 mediates gene expression events consistent with transformation to the repair phenotype. PACSIN1 promoted SC migration and viability following the TNFα challenge. When Src family kinases were pharmacologically inhibited or the receptor tyrosine kinase, TrkC, was genetically silenced or pharmacologically inhibited, PACSIN1 failed to induce cell signaling and prevent SC death. Collectively, these studies demonstrate that PACSIN1 is a novel axon-derived LRP1 ligand that activates SC repair signaling by transactivating TrkC.

The genomic landscape of Ménière's disease: a path to endolymphatic hydrops.

BACKGROUND: Ménière's disease (MD) is a disorder of the inner ear that causes episodic bouts of severe dizziness, roaring tinnitus, and fluctuating hearing loss. To date, no targeted therapy exists. As such, we have undertaken a large whole genome sequencing study on carefully phenotyped unilateral MD patients with the goal of gene/pathway discovery and a move towards targeted intervention. This study was a retrospective review of patients with a history of Ménière's disease. Genomic DNA, acquired from saliva samples, was purified and subjected to whole genome sequencing. RESULTS: Stringent variant calling, performed on 511 samples passing quality checks, followed by gene-based filtering by recurrence and proximity in molecular interaction networks, led to 481 high priority MD genes. These high priority genes, including MPHOSPH8, MYO18A, TRIOBP, OTOGL, TNC, and MYO6, were previously implicated in hearing loss, balance, and cochlear function, and were significantly enriched in common variant studies of hearing loss. Validation in an independent MD cohort confirmed 82 recurrent genes. Pathway analysis pointed to cell-cell adhesion, extracellular matrix, and cellular energy maintenance as key mediators of MD. Furthermore, the MD-prioritized genes were highly expressed in human inner ear hair cells and dark/vestibular cells, and were differentially expressed in a mouse model of hearing loss. CONCLUSION: By enabling the development of model systems that may lead to targeted therapies and MD screening panels, the genes and variants identified in this study will inform diagnosis and treatment of MD.

Nonalcoholic fatty liver disease risk and histologic severity are associated with genetic polymorphisms in children.

BACKGROUND AND AIMS: NAFLD is the most common chronic liver disease in children. Large pediatric studies identifying single nucleotide polymorphisms (SNPs) associated with risk and histologic severity of NAFLD are limited. Study aims included investigating SNPs associated with risk for NAFLD using family trios and association of candidate alleles with histologic severity. APPROACH AND RESULTS: Children with biopsy-confirmed NAFLD were enrolled from the NASH Clinical Research Network. The Expert Pathology Committee reviewed liver histology. Genotyping was conducted with allele-specific primers for 60 candidate SNPs. Parents were enrolled for trio analysis. To assess risk for NAFLD, the transmission disequilibrium test was conducted in trios. Among cases, regression analysis assessed associations with histologic severity. A total of 822 children with NAFLD had mean age 13.2 years (SD 2.7) and mean ALT 101 U/L (SD 90). PNPLA3 (rs738409) demonstrated the strongest risk ( p = 2.24 × 10 -14 ) for NAFLD. Among children with NAFLD, stratifying by PNPLA3 s738409 genotype, the variant genotype associated with steatosis ( p = 0.005), lobular ( p = 0.03) and portal inflammation ( p = 0.002). Steatosis grade associated with TM6SF2 ( p = 0.0009), GCKR ( p = 0.0032), PNPLA3 rs738409 ( p = 0.0053), and MTTP ( p = 0.0051). Fibrosis stage associated with PARVB rs6006473 ( p = 0.0001), NR1I2 ( p = 0.0021), ADIPOR2 ( p = 0.0038), and OXTR ( p = 0.0065). PNPLA3 rs738409 ( p = 0.0002) associated with borderline zone 1 NASH. CONCLUSIONS: This study demonstrated disease-associated SNPs in children with NAFLD. In particular, rs6006473 was highly associated with severity of fibrosis. These hypothesis-generating results support future mechanistic studies of development of adverse outcomes such as fibrosis and generation of therapeutic targets for NAFLD in children.

Focal adhesion is associated with lithium response in bipolar disorder: evidence from a network-based multi-omics analysis.

Lithium (Li) is one of the most effective drugs for treating bipolar disorder (BD), however, there is presently no way to predict response to guide treatment. The aim of this study is to identify functional genes and pathways that distinguish BD Li responders (LR) from BD Li non-responders (NR). An initial Pharmacogenomics of Bipolar Disorder study (PGBD) GWAS of lithium response did not provide any significant results. As a result, we then employed network-based integrative analysis of transcriptomic and genomic data. In transcriptomic study of iPSC-derived neurons, 41 significantly differentially expressed (DE) genes were identified in LR vs NR regardless of lithium exposure. In the PGBD, post-GWAS gene prioritization using the GWA-boosting (GWAB) approach identified 1119 candidate genes. Following DE-derived network propagation, there was a highly significant overlap of genes between the top 500- and top 2000-proximal gene networks and the GWAB gene list (Phypergeometric = 1.28E-09 and 4.10E-18, respectively). Functional enrichment analyses of the top 500 proximal network genes identified focal adhesion and the extracellular matrix (ECM) as the most significant functions. Our findings suggest that the difference between LR and NR was a much greater effect than that of lithium. The direct impact of dysregulation of focal adhesion on axon guidance and neuronal circuits could underpin mechanisms of response to lithium, as well as underlying BD. It also highlights the power of integrative multi-omics analysis of transcriptomic and genomic profiling to gain molecular insights into lithium response in BD.

Immunotherapy-based targeting of MSLN+ activated portal fibroblasts is a strategy for treatment of cholestatic liver fibrosis.

We investigated the role of mesothelin (Msln) and thymocyte differentiation antigen 1 (Thy1) in the activation of fibroblasts across multiple organs and demonstrated that Msln-/- mice are protected from cholestatic fibrosis caused by Mdr2 (multidrug resistance gene 2) deficiency, bleomycin-induced lung fibrosis, and UUO (unilateral urinary obstruction)-induced kidney fibrosis. On the contrary, Thy1-/- mice are more susceptible to fibrosis, suggesting that a Msln-Thy1 signaling complex is critical for tissue fibroblast activation. A similar mechanism was observed in human activated portal fibroblasts (aPFs). Targeting of human MSLN+ aPFs with two anti-MSLN immunotoxins killed fibroblasts engineered to express human mesothelin and reduced collagen deposition in livers of bile duct ligation (BDL)-injured mice. We provide evidence that antimesothelin-based therapy may be a strategy for treatment of parenchymal organ fibrosis.

Mapping the gene network landscape of Alzheimer's disease through integrating genomics and transcriptomics.

Integration of multi-omics data with molecular interaction networks enables elucidation of the pathophysiology of Alzheimer's disease (AD). Using the latest genome-wide association studies (GWAS) including proxy cases and the STRING interactome, we identified an AD network of 142 risk genes and 646 network-proximal genes, many of which were linked to synaptic functions annotated by mouse knockout data. The proximal genes were confirmed to be enriched in a replication GWAS of autopsy-documented cases. By integrating the AD gene network with transcriptomic data of AD and healthy temporal cortices, we identified 17 gene clusters of pathways, such as up-regulated complement activation and lipid metabolism, down-regulated cholinergic activity, and dysregulated RNA metabolism and proteostasis. The relationships among these pathways were further organized by a hierarchy of the AD network pinpointing major parent nodes in graph structure including endocytosis and immune reaction. Control analyses were performed using transcriptomics from cerebellum and a brain-specific interactome. Further integration with cell-specific RNA sequencing data demonstrated genes in our clusters of immunoregulation and complement activation were highly expressed in microglia.

Nondegradable Collagen Increases Liver Fibrosis but Not Hepatocellular Carcinoma in Mice.

Although hepatocellular cancer (HCC) usually occurs in the setting of liver fibrosis, the causal relationship between liver fibrosis and HCC is unclear. in vivo and in vitro models of HCC involving Colr/r mice (that produce a collagenase-resistant type I collagen) or wild-type (WT) mice were used to assess the relationship between type I collagen, liver fibrosis, and experimental HCC. HCC was either chemically induced in WT and Colr/r mice or Hepa 1-6 cells were engrafted into WT and Colr/r livers. The effect of hepatic stellate cells (HSCs) from WT and Colr/r mice on the growth of Hepa 1-6 cells was studied by using multicellular tumor spheroids and xenografts. Collagen type I deposition and fibrosis were increased in Colr/r mice, but they developed fewer and smaller tumors. Hepa 1-6 cells had reduced tumor growth in the livers of Colr/r mice. Although Colr/r HSCs exhibited a more activated phenotype, Hepa 1-6 growth and malignancy were suppressed in multicellular tumor spheroids and in xenografts containing Colr/r HSCs. Treatment with vitronectin, which mimics the presence of degraded collagen fragments, converted the Colr/r phenotype into a WT phenotype. Although Colr/r mice have increased liver fibrosis, they exhibited decreased HCC in several models. Thus, increased liver type I collagen does not produce increased experimental HCC.

Heterogeneity of HSCs in a Mouse Model of NASH.

BACKGROUND AND AIMS: In clinical and experimental NASH, the origin of the scar-forming myofibroblast is the HSC. We used foz/foz mice on a Western diet to characterize in detail the phenotypic changes of HSCs in a NASH model. APPROACH AND RESULTS: We examined the single-cell expression profiles (scRNA sequencing) of HSCs purified from the normal livers of foz/foz mice on a chow diet, in NASH with fibrosis of foz/foz mice on a Western diet, and in livers during regression of NASH after switching back to a chow diet. Selected genes were analyzed using immunohistochemistry, quantitative real-time PCR, and short hairpin RNA knockdown in primary mouse HSCs. Our analysis of the normal liver identified two distinct clusters of quiescent HSCs that correspond to their acinar position of either pericentral vein or periportal vein. The NASH livers had four distinct HSC clusters, including one representing the classic fibrogenic myofibroblast. The three other HSC clusters consisted of a proliferating cluster, an intermediate activated cluster, and an immune and inflammatory cluster. The livers with NASH regression had one cluster of inactivated HSCs, which was similar to, but distinct from, the quiescent HSCs. CONCLUSIONS: Analysis of single-cell RNA sequencing in combination with an interrogation of previous studies revealed an unanticipated heterogeneity of HSC phenotypes under normal and injured states.

Identification of Lineage-Specific Transcription Factors That Prevent Activation of Hepatic Stellate Cells and Promote Fibrosis Resolution.

BACKGROUND & AIMS: Development of liver fibrosis is associated with activation of quiescent hepatic stellate cells (HSCs) into collagen type I-producing myofibroblasts (activated HSCs). Cessation of liver injury often results in fibrosis resolution and inactivation of activated HSCs/myofibroblasts into a quiescent-like state (inactivated HSCs). We aimed to identify molecular features of phenotypes of HSCs from mice and humans. METHODS: We performed studies with LratCre, Ets1-floxed, Nf1-floxed, Pparγ-floxed, Gata6-floxed, Rag2-/-γc-/-, and C57/Bl6 (control) mice. Some mice were given carbon tetrachloride (CCl4) to induce liver fibrosis, with or without a peroxisome proliferator-activated receptor-γ (PPARγ) agonist. Livers from mice were analyzed by immunohistochemistry. Quiescent, activated, and inactivated HSCs were isolated from livers of Col1α1YFP mice and analyzed by chromatin immunoprecipitation and sequencing. Human HSCs were isolated from livers denied for transplantation. We compared changes in gene expression patterns and epigenetic modifications (histone H3 lysine 4 dimethylation and histone H3 lysine 27 acetylation) in primary mouse and human HSCs. Transcription factors were knocked down with small hairpin RNAs in mouse HSCs. RESULTS: Motif enrichment identified E26 transcription-specific transcription factors (ETS) 1, ETS2, GATA4, GATA6, interferon regulatory factor (IRF) 1, and IRF2 transcription factors as regulators of the mouse and human HSC lineage. Small hairpin RNA-knockdown of these transcription factors resulted in increased expression of genes that promote fibrogenesis and inflammation, and loss of HSC phenotype. Disruption of Gata6 or Ets1, or Nf1 or Pparγ (which are regulated by ETS1), increased the severity of CCl4-induced liver fibrosis in mice compared to control mice. Only mice with disruption of Gata6 or Pparγ had defects in fibrosis resolution after CCl4 administration was stopped, associated with persistent activation of HSCs. Administration of a PPARγ agonist accelerated regression of liver fibrosis after CCl4 administration in control mice but not in mice with disruption of Pparγ. CONCLUSIONS: Phenotypes of HSCs from humans and mice are regulated by transcription factors, including ETS1, ETS2, GATA4, GATA6, IRF1, and IRF2. Activated mouse and human HSCs can revert to a quiescent-like, inactivated phenotype. We found GATA6 and PPARγ to be required for inactivation of human HSCs and regression of liver fibrosis in mice.

IL-17 signaling in steatotic hepatocytes and macrophages promotes hepatocellular carcinoma in alcohol-related liver disease.

BACKGROUND & AIMS: Chronic alcohol consumption is a leading risk factor for the development of hepatocellular carcinoma (HCC), which is associated with a marked increase in hepatic expression of pro-inflammatory IL-17A and its receptor IL-17RA. METHODS: Genetic deletion and pharmacological blocking were used to characterize the role of IL-17A/IL-17RA signaling in the pathogenesis of HCC in mouse models and human specimens. RESULTS: We demonstrate that the global deletion of the Il-17ra gene suppressed HCC in alcohol-fed diethylnitrosamine-challenged Il-17ra-/- and major urinary protein-urokinase-type plasminogen activator/Il-17ra-/- mice compared with wild-type mice. When the cell-specific role of IL-17RA signaling was examined, the development of HCC was decreased in both alcohol-fed Il-17raΔMΦ and Il-17raΔHep mice devoid of IL-17RA in myeloid cells and hepatocytes, but not in Il-17raΔHSC mice (deficient in IL-17RA in hepatic stellate cells). Deletion of Il-17ra in myeloid cells ameliorated tumorigenesis via suppression of pro-tumorigenic/inflammatory and pro-fibrogenic responses in alcohol-fed Il-17raΔMΦ mice. Remarkably, despite a normal inflammatory response, alcohol-fed Il-17raΔHep mice developed the fewest tumors (compared with Il-17raΔMΦ mice), with reduced steatosis and fibrosis. Steatotic IL-17RA-deficient hepatocytes downregulated the expression of Cxcl1 and other chemokines, exhibited a striking defect in tumor necrosis factor (TNF)/TNF receptor 1-dependent caspase-2-SREBP1/2-DHCR7-mediated cholesterol synthesis, and upregulated the production of antioxidant vitamin D3. The pharmacological blocking of IL-17A/Th-17 cells using anti-IL-12/IL-23 antibodies suppressed the progression of HCC (by 70%) in alcohol-fed mice, indicating that targeting IL-17 signaling might provide novel strategies for the treatment of alcohol-induced HCC. CONCLUSIONS: Overall, IL-17A is a tumor-promoting cytokine, which critically regulates alcohol-induced hepatic steatosis, inflammation, fibrosis, and HCC. LAY SUMMARY: IL-17A is a tumor-promoting cytokine, which critically regulates inflammatory responses in macrophages (Kupffer cells and bone-marrow-derived monocytes) and cholesterol synthesis in steatotic hepatocytes in an experimental model of alcohol-induced HCC. Therefore, IL-17A may be a potential therapeutic target for patients with alcohol-induced HCC.

Activated hepatic stellate cells and portal fibroblasts contribute to cholestatic liver fibrosis in MDR2 knockout mice.

BACKGROUND & AIMS: Chronic liver injury often results in the activation of hepatic myofibroblasts and the development of liver fibrosis. Hepatic myofibroblasts may originate from 3 major sources: hepatic stellate cells (HSCs), portal fibroblasts (PFs), and fibrocytes, with varying contributions depending on the etiology of liver injury. Here, we assessed the composition of hepatic myofibroblasts in multidrug resistance gene 2 knockout (Mdr2-/-) mice, a genetic model that resembles primary sclerosing cholangitis in patients. METHODS: Mdr2-/- mice expressing a collagen-GFP reporter were analyzed at different ages. Hepatic non-parenchymal cells isolated from collagen-GFP Mdr2-/- mice were sorted based on collagen-GFP and vitamin A. An NADPH oxidase (NOX) 1/4 inhibitor was administrated to Mdr2-/- mice aged 12-16 weeks old to assess the therapeutic approach of targeting oxidative stress in cholestatic injury. RESULTS: Thy1+ activated PFs accounted for 26%, 51%, and 54% of collagen-GFP+ myofibroblasts in Mdr2-/- mice at 4, 8, and 16 weeks of age, respectively. The remaining collagen-GFP+ myofibroblasts were composed of activated HSCs, suggesting that PFs and HSCs are both activated in Mdr2-/- mice. Bone-marrow-derived fibrocytes minimally contributed to liver fibrosis in Mdr2-/- mice. The development of cholestatic liver fibrosis in Mdr2-/- mice was associated with early recruitment of Gr1+ myeloid cells and upregulation of pro-inflammatory cytokines (4 weeks). Administration of a NOX inhibitor to 12-week-old Mdr2-/- mice suppressed the activation of myofibroblasts and attenuated the development of cholestatic fibrosis. CONCLUSIONS: Activated PFs and activated HSCs contribute to cholestatic fibrosis in Mdr2-/- mice, and serve as targets for antifibrotic therapy. LAY SUMMARY: Activated portal fibroblasts and hepatic stellate cells, but not fibrocytes, contributed to the production of the fibrous scar in livers of Mdr2-/- mice, and these cells can serve as targets for antifibrotic therapy in cholestatic injury. Therapeutic inhibition of the enzyme NADPH oxidase (NOX) in Mdr2-/- mice reversed cholestatic fibrosis, suggesting that targeting NOXs may be an effective strategy for the treatment of cholestatic fibrosis.

Interactive network visualization in Jupyter notebooks: visJS2jupyter.

Motivation: Network biology is widely used to elucidate mechanisms of disease and biological processes. The ability to interact with biological networks is important for hypothesis generation and to give researchers an intuitive understanding of the data. We present visJS2jupyter, a tool designed to embed interactive networks in Jupyter notebooks to streamline network analysis and to promote reproducible research. Results: The tool provides functions for performing and visualizing useful network operations in biology, including network overlap, network propagation around a focal set of genes, and co-localization of two sets of seed genes. visJS2jupyter uses the JavaScript library vis.js to create interactive networks displayed within Jupyter notebook cells with features including drag, click, hover, and zoom. We demonstrate the functionality of visJS2jupyter applied to a biological question, by creating a network propagation visualization to prioritize risk-related genes in autism. Availability and implementation: The visJS2jupyter package is distributed under the MIT License. The source code, documentation and installation instructions are freely available on GitHub at https://github.com/ucsd-ccbb/visJS2jupyter. The package can be downloaded at https://pypi.python.org/pypi/visJS2jupyter. Contact: sbrosenthal@ucsd.edu. Supplementary information: Supplementary data are available at Bioinformatics online.

Revealing the hidden networks of interaction in mobile animal groups allows prediction of complex behavioral contagion.

Coordination among social animals requires rapid and efficient transfer of information among individuals, which may depend crucially on the underlying structure of the communication network. Establishing the decision-making circuits and networks that give rise to individual behavior has been a central goal of neuroscience. However, the analogous problem of determining the structure of the communication network among organisms that gives rise to coordinated collective behavior, such as is exhibited by schooling fish and flocking birds, has remained almost entirely neglected. Here, we study collective evasion maneuvers, manifested through rapid waves, or cascades, of behavioral change (a ubiquitous behavior among taxa) in schooling fish (Notemigonus crysoleucas). We automatically track the positions and body postures, calculate visual fields of all individuals in schools of ∼150 fish, and determine the functional mapping between socially generated sensory input and motor response during collective evasion. We find that individuals use simple, robust measures to assess behavioral changes in neighbors, and that the resulting networks by which behavior propagates throughout groups are complex, being weighted, directed, and heterogeneous. By studying these interaction networks, we reveal the (complex, fractional) nature of social contagion and establish that individuals with relatively few, but strongly connected, neighbors are both most socially influential and most susceptible to social influence. Furthermore, we demonstrate that we can predict complex cascades of behavioral change at their moment of initiation, before they actually occur. Consequently, despite the intrinsic stochasticity of individual behavior, establishing the hidden communication networks in large self-organized groups facilitates a quantitative understanding of behavioral contagion.