Search

Hematopoietic stem cells (HSCs) derived from human induced pluripotent stem cells (iPS cells) have important biomedical applications. We identified differentiation conditions that generate HSCs defined by robust long-term multilineage engraftment in immune-deficient NOD,B6.Prkdc scid Il2rg tm1Wjl/SzJ Kit W41/W41 mice. We guided differentiating iPS cells, as embryoid bodies in a defined culture medium supplemented with retinyl acetate, through HOXA-patterned mesoderm to hemogenic endothelium specified by bone morphogenetic protein 4 and vascular endothelial growth factor (VEGF). Removal of VEGF facilitated an efficient endothelial-to-hematopoietic transition, evidenced by release into the culture medium of CD34 + blood cells, which were cryopreserved. Intravenous transplantation of two million thawed CD34 + cells differentiated from four independent iPS cell lines produced multilineage bone marrow engraftment in 25-50% of immune-deficient recipient mice. These functionally defined, multipotent CD34 + hematopoietic cells, designated iPS cell-derived HSCs (iHSCs), produced levels of engraftment similar to those achieved following umbilical cord blood transplantation. Our study provides a step toward the goal of generating HSCs for clinical translation., (© 2024. The Author(s).)

Cardiac function requires appropriate proteins in each chamber. Atria requires slow myosin to act as reservoirs, while ventricles demand fast myosin for swift pumping. Myosins are thus under chamber-biased cis-regulation, with myosin gene expression imbalances leading to congenital heart dysfunction. To identify regulatory inputs leading to cardiac chamber-biased expression, we computationally and molecularly dissected the quail Slow Myosin Heavy Chain III (SMyHC III) promoter that drives preferential expression to the atria. We show that SMyHC III gene states are orchestrated by a complex Nuclear Receptor Element (cNRE) of 32 base pairs. Using transgenesis in zebrafish and mice, we demonstrate that preferential atrial expression is achieved by a combinatorial regulatory input composed of atrial activation motifs and ventricular repression motifs. Using comparative genomics, we show that the cNRE might have emerged from an endogenous viral element through infection of an ancestral host germline, revealing an evolutionary pathway to cardiac chamber-specific expression., (© 2024. The Author(s).)

Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) primarily infects the respiratory tract, but pulmonary and cardiac complications occur in severe coronavirus disease 2019 (COVID-19). To elucidate molecular mechanisms in the lung and heart, we conducted paired experiments in human stem cell-derived lung alveolar type II (AT2) epithelial cell and cardiac cultures infected with SARS-CoV-2. With CRISPR-Cas9-mediated knockout of ACE2, we demonstrated that angiotensin-converting enzyme 2 (ACE2) was essential for SARS-CoV-2 infection of both cell types but that further processing in lung cells required TMPRSS2, while cardiac cells required the endosomal pathway. Host responses were significantly different; transcriptome profiling and phosphoproteomics responses depended strongly on the cell type. We identified several antiviral compounds with distinct antiviral and toxicity profiles in lung AT2 and cardiac cells, highlighting the importance of using several relevant cell types for evaluation of antiviral drugs. Our data provide new insights into rational drug combinations for effective treatment of a virus that affects multiple organ systems., Competing Interests: Conflict of interests E.P. is a co-founder and scientific advisor and holds equity in, and N.C. has a paid position with, Dynomics, a biotechnology company focused on the development of heart failure therapeutics. J.M.P. is a co-founder and shareholder of Mogrify, Ltd., a cell therapy company., (Copyright © 2023 The Authors. Published by Elsevier Inc. All rights reserved.)

SARS-CoV-2 primarily infects the respiratory tract, but pulmonary and cardiac complications occur in severe COVID-19. To elucidate molecular mechanisms in the lung and heart, we conducted paired experiments in human stem cell-derived lung alveolar type II (AT2) epithelial cell and cardiac cultures infected with SARS-CoV-2. With CRISPR- Cas9 mediated knock-out of ACE2, we demonstrated that angiotensin converting enzyme 2 (ACE2) was essential for SARS-CoV-2 infection of both cell types but further processing in lung cells required TMPRSS2 while cardiac cells required the endosomal pathway. Host responses were significantly different; transcriptome profiling and phosphoproteomics responses depended strongly on the cell type. We identified several antiviral compounds with distinct antiviral and toxicity profiles in lung AT2 and cardiac cells, highlighting the importance of using several relevant cell types for evaluation of antiviral drugs. Our data provide new insights into rational drug combinations for effective treatment of a virus that affects multiple organ systems., One-Sentence Summary: Rational treatment strategies for SARS-CoV-2 derived from human PSC models.

Human macrophages are a natural host of many mycobacterium species, including Mycobacterium abscessus (M. abscessus), an emerging pathogen affecting immunocompromised and cystic fibrosis patients with few available treatments. The search for an effective treatment is hindered by the lack of a tractable in vitro intracellular infection model. Here, we established a reliable model for M. abscessus infection using human pluripotent stem cell-derived macrophages (hPSC-macrophages). hPSC differentiation permitted reproducible generation of functional macrophages that were highly susceptible to M. abscessus infection. Electron microscopy demonstrated that M. abscessus was present in the hPSC-macrophage vacuoles. RNA sequencing analysis revealed a time-dependent host cell response, with differing gene and protein expression patterns post-infection. Engineered tdTOMATO-expressing hPSC-macrophages with GFP-expressing mycobacteria enabled rapid image-based high-throughput analysis of intracellular infection and quantitative assessment of antibiotic efficacy. Our study describes the first to our knowledge hPSC-based model for M. abscessus infection, representing a novel and accessible system for studying pathogen-host interaction and drug discovery., Competing Interests: Conflicts of interest E.G.S. is member of the editorial board of Stem Cell Reports., (Copyright © 2022 The Author(s). Published by Elsevier Inc. All rights reserved.)

Thymic epithelium is critical for the structural integrity of the thymus and for T cell development. Within the fully formed thymus, large numbers of hematopoietic cells shape the thymic epithelium into a scaffold-like structure which bears little similarity to classical epithelial layers, such as those observed in the skin, intestine or pancreas. Here, we show that human thymic epithelial cells (TECs) possess an epithelial identity that also incorporates the expression of mesenchymal cell associated genes, whose expression levels vary between medullary and cortical TECs (m/cTECs). Using pluripotent stem cell (PSC) differentiation systems, we identified a unique population of cells that co-expressed the master TEC transcription factor FOXN1 , as well as the epithelial associated marker EPCAM and the mesenchymal associated gene CD90. Using the same serum free culture conditions, we also observed co-expression of EPCAM and CD90 on cultured TECs derived from neonatal human thymus in vitro . Single cell RNA-sequencing revealed these cultured TECs possessed an immature mTEC phenotype and expressed epithelial and mesenchymal associated genes, such as EPCAM , CLDN4 , CD90 and COL1A1 . Importantly, flow cytometry and single cell RNA-sequencing analysis further confirmed the presence of an EPCAM+CD90+ population in the CD45- fraction of neonatal human thymic stromal cells in vivo . Using the human thymus cell atlas, we found that cTECs displayed more pronounced mesenchymal characteristics than mTECs during embryonic development. Collectively, these results suggest human TECs possess a hybrid gene expression program comprising both epithelial and mesenchymal elements, and provide a basis for the further exploration of thymus development from primary tissues and from the in vitro differentiation of PSCs., Competing Interests: The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest., (Copyright © 2022 Sun, Li, Nim, Piers, Ramialison, Porrello, Konstantinov, Elefanty and Stanley.)

Organ fibroblasts are essential components of homeostatic and diseased tissues. They participate in sculpting the extracellular matrix, sensing the microenvironment, and communicating with other resident cells. Recent studies have revealed transcriptomic heterogeneity among fibroblasts within and between organs. To dissect the basis of interorgan heterogeneity, we compare the gene expression of murine fibroblasts from different tissues (tail, skin, lung, liver, heart, kidney, and gonads) and show that they display distinct positional and organ-specific transcriptome signatures that reflect their embryonic origins. We demonstrate that expression of genes typically attributed to the surrounding parenchyma by fibroblasts is established in embryonic development and largely maintained in culture, bioengineered tissues and ectopic transplants. Targeted knockdown of key organ-specific transcription factors affects fibroblast functions, in particular genes involved in the modulation of fibrosis and inflammation. In conclusion, our data reveal that adult fibroblasts maintain an embryonic gene expression signature inherited from their organ of origin, thereby increasing our understanding of adult fibroblast heterogeneity. The knowledge of this tissue-specific gene signature may assist in targeting fibrotic diseases in a more precise, organ-specific manner., Competing Interests: EF, MR, HN, MM, JL, RC, SD, SB, ES, AE, JH, MC, NR, MF No competing interests declared, (© 2022, Forte et al.)

Background: Gene ontology (GO) enrichment analysis is frequently undertaken during exploration of various -omics data sets. Despite the wide array of tools available to biologists to perform this analysis, meaningful visualisation of the overrepresented GO in a manner which is easy to interpret is still lacking., Results: Monash Gene Ontology (MonaGO) is a novel web-based visualisation system that provides an intuitive, interactive and responsive interface for performing GO enrichment analysis and visualising the results. MonaGO supports gene lists as well as GO terms as inputs. Visualisation results can be exported as high-resolution images or restored in new sessions, allowing reproducibility of the analysis. An extensive comparison between MonaGO and 11 state-of-the-art GO enrichment visualisation tools based on 9 features revealed that MonaGO is a unique platform that simultaneously allows interactive visualisation within one single output page, directly accessible through a web browser with customisable display options., Conclusion: MonaGO combines dynamic clustering and interactive visualisation as well as customisation options to assist biologists in obtaining meaningful representation of overrepresented GO terms, producing simplified outputs in an unbiased manner. MonaGO will facilitate the interpretation of GO analysis and will assist the biologists into the representation of the results., (© 2022. The Author(s).)

Understanding the spatial gene expression and regulation in the heart is key to uncovering its developmental and physiological processes, during homeostasis and disease. Numerous techniques exist to gain gene expression and regulation information in organs such as the heart, but few utilize intuitive true-to-life three-dimensional representations to analyze and visualise results. Here we combined transcriptomics with 3D-modelling to interrogate spatial gene expression in the mammalian heart. For this, we microdissected and sequenced transcriptome-wide 18 anatomical sections of the adult mouse heart. Our study has unveiled known and novel genes that display complex spatial expression in the heart sub-compartments. We have also created 3D-cardiomics, an interface for spatial transcriptome analysis and visualization that allows the easy exploration of these data in a 3D model of the heart. 3D-cardiomics is accessible from http://3d-cardiomics.erc.monash.edu/., (Copyright © 2021. Published by Elsevier Ltd.)

Spatially resolved transcriptomics is an emerging class of high-throughput technologies that enable biologists to systematically investigate the expression of genes along with spatial information. Upon data acquisition, one major hurdle is the subsequent interpretation and visualization of the datasets acquired. To address this challenge, VR-Cardiomics is presented, which is a novel data visualization system with interactive functionalities designed to help biologists interpret spatially resolved transcriptomic datasets. By implementing the system in two separate immersive environments, fish tank virtual reality (FTVR) and head-mounted display virtual reality (HMD-VR), biologists can interact with the data in novel ways not previously possible, such as visually exploring the gene expression patterns of an organ, and comparing genes based on their 3D expression profiles. Further, a biologist-driven use-case is presented, in which immersive environments facilitate biologists to explore and compare the heart expression profiles of different genes., (© 2022. The Author(s).)

Background: Congenital heart diseases are the major cause of death in newborns, but the genetic etiology of this developmental disorder is not fully known. The conventional approach to identify the disease-causing genes focuses on screening genes that display heart-specific expression during development. However, this approach would have discounted genes that are expressed widely in other tissues but may play critical roles in heart development., Results: We report an efficient pipeline of genome-wide gene discovery based on the identification of a cardiac-specific cis-regulatory element signature that points to candidate genes involved in heart development and congenital heart disease. With this pipeline, we retrieve 76% of the known cardiac developmental genes and predict 35 novel genes that previously had no known connectivity to heart development. Functional validation of these novel cardiac genes by RNAi-mediated knockdown of the conserved orthologs in Drosophila cardiac tissue reveals that disrupting the activity of 71% of these genes leads to adult mortality. Among these genes, RpL14, RpS24, and Rpn8 are associated with heart phenotypes., Conclusions: Our pipeline has enabled the discovery of novel genes with roles in heart development. This workflow, which relies on screening for non-coding cis-regulatory signatures, is amenable for identifying developmental and disease genes for an organ without constraining to genes that are expressed exclusively in the organ of interest., (© 2021. The Author(s).)

Background: The ability to identify lupus patients in high disease activity status (HDAS) without knowledge of the SLEDAI could have application in selection of patients for treatment escalation or enrolment in trials. We sought to generate an algorithm that could calculate via model fitting the presence of HDAS using simple demographic and laboratory values., Methods: We examined the association of high disease activity (HDA) with demographic and laboratory parameters using prospectively collected data. An HDA visit is recorded when SLEDAI-2K ≥10. We utilized the use of combinatorial search to find algorithms to build a mathematical model predictive of HDA. Performance of each algorithm was evaluated using multi-class area under the receiver operating characteristic curve and the final model was compared with the naïve Bayes classifier, and analysed using the confusion matrix for accuracy and misclassification rate., Results: Data on 286 patients, followed for a median of 5.1 years were studied for a total of 5680 visits. Sixteen laboratory parameters were found to be significantly associated with HDA. A total of 216 algorithms were evaluated and the final algorithm chosen was based on seven pathology measures and three demographic variables. It has an accuracy of 88.6% and misclassification rate of 11.4%. When compared with the naïve Bayes classifier [area under the curve (AUC) = 0.663], our algorithm has a better accuracy with AUC = 0.829., Conclusion: This study shows that building an accurate model to calculate HDA using routinely available clinical parameters is feasible. Future studies to independently validate the algorithm will be needed to confirm its predictive performance., (© The Author(s) 2021. Published by Oxford University Press on behalf of the British Society for Rheumatology. All rights reserved. For permissions, please email: journals.permissions@oup.com.)

Background: We sought to examine the disease course of High Disease Activity Status (HDAS) patients and their different disease patterns in a real-world longitudinal cohort. Disease resolution till Lupus Low Disease Activity State (LLDAS) has been a general treatment goal, but there is limited information on this subset of patients who achieve this., Methods: All consenting patients of the Monash Lupus Cohort who had at least 12 months of observation were included. HDAS was defined as SLEDAI-2K ≥ 10 ever, and HDAS episode as the period from the first HDAS clinic visit until attainment of LLDAS. We examined the associations of different HDAS patterns with the likelihood of damage accrual., Results: Of 342 SLE patients, 151 experienced HDAS at least once, accounting for 298 HDAS episodes. The majority of HDAS patients (76.2%) experienced Recurrent HDAS (> 1 HDAS visit), and a smaller subset (47.7%) had Persistent HDAS (consecutive HDAS visits for longer than 2 months). Recurrent or Persistent HDAS patients were younger at diagnosis and more likely to experience renal or serositis manifestations; persistent HDAS patients were also more likely to experience neurological manifestations. Baseline SLEDAI greater than 10 was associated with longer HDAS episodes. Recurrent and Persistent HDAS were both associated with an increased likelihood of damage accrual. The total duration of HDAS episode greater than 2 years and experiencing multiple HDAS episodes (≥4) was also associated with an increased likelihood of damage accrual (OR 1.80, 95% CI 1.08-2.97, p = 0.02, and OR 3.31, 95% CI 1.66-13.26, p = 0.01, respectively)., Conclusion: HDAS episodes have a highly variable course. Recurrent and Persistent HDAS, and longer duration of HDAS episodes, increased the risk of damage accrual. In addition to a major signifier of severity in SLE, its resolution to LLDAS can determine the subsequent outcome in SLE patients., (© 2021. The Author(s).)

Heart formation in the zebrafish involves a rapid, complex series of morphogenetic events in three-dimensional space that spans cardiac lineage specification through to chamber formation and maturation. This process is tightly orchestrated by a cardiac gene regulatory network (GRN), which ensures the precise spatio-temporal deployment of genes critical for heart formation. Alterations of the timing or spatial localisation of gene expression can have a significant impact in cardiac ontogeny and may lead to heart malformations. Hence, a better understanding of the cellular and molecular basis of congenital heart disease relies on understanding the behaviour of cardiac GRNs with precise spatiotemporal resolution. Here, we review the recent technical advances that have expanded our capacity to interrogate the cardiac GRN in zebrafish. In particular, we focus on studies utilising high-throughput technologies to systematically dissect gene expression patterns, both temporally and spatially during heart development., (© The Author(s) 2021. Published by Oxford University Press. All rights reserved. For Permissions, please email: journals.permissions@oup.com.)

Background: Glucocorticoids, used as a therapy in systemic lupus erythematosus (SLE), interact with the cytoplasmic glucocorticoid receptor to modulate gene transcription. Minimising the use of glucocorticoids is a goal in SLE; however, pharmacological measures to support clinical guidelines are scarce. We evaluated glucocorticoid-regulated genes for their potential use as biomarkers of glucocorticoid exposure in SLE. We examined interactions between changes in gene expression that are induced by glucocorticoids and type I interferon., Methods: Genes regulated by glucocorticoids and type I interferon were analysed in relation to glucocorticoid exposure in adult patients meeting the American College of Rheumatology criteria for SLE from three cross-sectional cohorts: a local cohort from a tertiary hospital in Melbourne, VIC, Australia, and two public datasets (GSE49454, Hospital de la Conception, Marseille, France, and GSE88884, patients enrolled in a large, multicentre clinical trial). RNA sequencing was done using RNA from healthy donor leucocytes treated with the glucocorticoid dexamethasone, or type I interferon, or both., Findings: Glucocorticoid-regulated genes were analysed in a local SLE cohort (n=18) and public dataset GSE49454 (n=62). Five genes correlated with glucocorticoid dose in both cohorts and were combined to make a glucocorticoid gene signature. Validity of the glucocorticoid gene signature was tested in the public dataset GSE88884 (n=1756). A dose-dependent association was observed with glucocorticoid dose (p<0·0001), and the glucocorticoid gene signature had moderate ability to identify patients taking high-dose glucocorticoid (area under the curve [AUC]=0·77) although was less discriminatory when including all doses (AUC=0·69). We saw no effect of glucocorticoid dose on type I interferon -regulated gene expression. Patients with a high type I interferon gene signature had reduced glucocorticoid gene signature expression compared with patients with a low type I interferon gene signature matched for glucocorticoid dose, suggesting type I interferon inhibits glucocorticoid-stimulated gene expression. In RNA sequencing experiments, type I interferon impaired the expression of glucocorticoid-induced genes, whereas dexamethasone had minimal effect on the expression of type I interferon-stimulated genes. We identified genes regulated by dexamethasone but not affected by type I interferon; combined signatures using these genes also showed moderate ability to distinguish patients taking glucocorticoids., Interpretation: A gene signature for glucocorticoid exposure was identified, but the substantial effect of type I interferon on glucocorticoid-induced genes might limit its application in SLE. These data confirm the insensitivity of type I interferon-regulated genes to glucocorticoids, and together support the concept that type I interferon has a role in glucocorticoid resistance in SLE., Funding: Lupus Research Alliance and Australian National Health and Medical Research Council., Competing Interests: Declaration of interests MN received funding from National Health and Medical Research Council (NHMRC), Australia. EFM received grants from the Lupus Research Alliance, New York, NY, USA, and NHMRC, Australia, for the funding of this work; grants and personal fees from AstraZeneca, Bristol Myers Squibb, Eli Lilly, Janssen, Wolf, GlaxoSmithKline, and EMD Serono; and personal fees from Biogen, AbbVie, Neovacs, UCB, Sanofi, Novartis, and Amgen, outside the submitted work. All other authors declare no competing interests., (Copyright © 2021 Elsevier Ltd. All rights reserved.)

Unravelling spatio-temporal patterns of gene expression is crucial to understanding core biological principles from embryogenesis to disease. Here we review emerging technologies, providing automated, high-throughput, spatially resolved quantitative gene expression data. Novel techniques expand on current benchmark protocols, expediting their incorporation into ongoing research. These approaches digitally reconstruct patterns of embryonic expression in three dimensions, and have successfully identified novel domains of expression, cell types, and tissue features. Such technologies pave the way for unbiased and exhaustive recapitulation of gene expression levels in spatial and quantitative terms, promoting understanding of the molecular origin of developmental defects, and improving medical diagnostics.

Objective: Disease severity in SLE is an important concept related to disease activity, treatment burden and prognosis. We set out to evaluate if high disease activity status (HDAS), based on ever attainment of a Systemic Lupus Erythematosus Disease Activity Index (SLEDAI) disease activity score of ≥10, is an indicator for disease severity in SLE., Methods: Using prospectively collected data, we assessed the association of HDAS with sociodemographic and disease characteristics and adverse clinical outcomes using logistic regression or generalised estimating equations., Results: Of 286 patients with SLE, who were observed for a median (range) of 5.1 years (1-10.8 years), 43.7% experienced HDAS at least once during the observational period. Autoantibody positivity, particularly anti-dsDNA and anti-Sm positivity, were associated with increased likelihood of HDAS. Age ≥45 years at diagnosis was associated with reduced likelihood of HDAS (p=0.002). Patients with HDAS had higher Physician Global Assessment score (>1: OR 8.1, p<0.001) and were more likely to meet criteria for flare (mild/moderate flare: OR 4.4, p<0.001; severe flare: OR 17.2, p<0.001) at the time of experiencing HDAS. They were also more likely to have overall higher disease activity, as defined by time-adjusted mean SLEDAI-2K score in the highest quartile (OR 11.7, 95% CI 5.1 to 26.6; p>0.001), higher corticosteroid exposure (corticosteroid dose in highest quartile: OR 7.7, 95% CI 3.9 to 15.3; p<0.001) and damage accrual (OR 2.3, 95% CI 1.3 to 3.9; p=0.003) when compared with non-HDAS patients., Conclusions: HDAS is associated with more severe disease, as measured by higher disease activity across time, corticosteroid exposure and damage accrual. The occurrence of HDAS may be a useful prognostic marker in the management of SLE., Competing Interests: Competing interests: None declared., (© Author(s) (or their employer(s)) 2020. Re-use permitted under CC BY-NC. No commercial re-use. See rights and permissions. Published by BMJ.)

Background: Treat-to-target strategies have improved outcomes in single-organ diseases with simple clinical or laboratory endpoints. A lack of validated endpoints has prevented adoption of treat to target for complex multiorgan conditions, such as systemic lupus erythematosus (SLE). We report the first prospective study undertaken to specifically validate a treat-to-target endpoint for SLE., Methods: In this prospective cohort study, patients aged 18 years or older with SLE were recruited from 13 centres in eight countries and followed prospectively. Patients with at least two visits over the study period no more than 6 months apart were included in the longitudinal analysis. Patients with no visits in the final year of the study were censored from their last visit. Attainment of the lupus low disease activity state (LLDAS) was assessed at each visit. The primary outcome measure was accrual of irreversible end-organ damage, defined as at least a 1-point increase in the Systemic Lupus International Collaborating Clinics/American College of Rheumatology Damage Index. We used time-dependent hazard regression models and generalised linear models to measure the association between LLDAS (attainment at any timepoint, cumulative time in LLDAS, and sustained LLDAS) with accrual of irreversible end-organ damage or flare (key secondary outcome). This study is registered with ClinicalTrials.gov, NCT03138941., Findings: Between May 1, 2013, and Dec 31, 2016, 1707 patients were recruited and followed for a mean of 2·2 years (SD 0·9), totalling 12 689 visits. Attainment of LLDAS at any timepoint was associated with reduction in damage accrual (hazard ratio 0·59, 0·45-0·76; p<0·0001) and subsequent flare (0·65, 95% CI 0·56-0·75; p<0·0001). Cumulative time in LLDAS was associated with improved outcomes: compared with patients with less than 50% of observed time in LLDAS, those with at least 50% of observed time in LLDAS had reduced risk of damage accrual (0·54, 0·42-0·70; p<0·0001) and flare (0·41, 0·35-0·48; p<0·0001). Similarly, increased durations of sustained LLDAS were associated with incremental reductions in the risk of damage accrual. The association of LLDAS with reduced damage accrual was observed regardless of pre-existing damage or disease activity at study entry., Interpretation: LLDAS attainment is associated with significant protection against flare and damage accrual in SLE. These findings validate LLDAS as an endpoint for clinical studies in SLE., Funding: The Asia Pacific Lupus Collaboration received project support grants from UCB Pharma, GlaxoSmithKline, Janssen, Bristol-Myers Squibb, and AstraZeneca., (Copyright © 2019 Elsevier Ltd. All rights reserved.)

Objective: Systemic lupus erythematosus (SLE) is a multisystem autoimmune disease. SLE is characterized by high inter-patient variability, including fluctuations over time, a factor which most biomarker studies omit from consideration. We investigated relationships between disease activity and biomarker expression in SLE, using novel methods to control for time-dependent variability, in a proof-of-concept study to evaluate whether doing so revealed additional information. Methods: We measured 4 serum biomarkers (MIF, CCL2, CCL19, and CXCL10) and 13 routine clinical laboratory parameters, alongside disease activity measured by the SLE disease activity index-2k (SLEDAI-2k), collected longitudinally. We analyzed these data with unsupervised learning methods via ensemble clustering, incorporating temporal relationships using dynamic time warping for distance metric calculation. Results: Data from 843 visits in 110 patients (median age 47, 83% female) demonstrated highly heterogeneous time-dependent relationships between disease activity and biomarkers. Unbiased magnitude-based hierarchical clustering of biomarker expression levels isolated a patient subset ( n = 9) with distinctively heterogeneous expression of the 17 biological parameters, and who had MIF, CCL2, CCL19, and CXCL10 levels that were higher and more strongly associated with disease activity, based on leave-one-out cross-validated regression analysis. In the remaining subgroup, a time-dependent regression model revealed significantly stronger predictive power of biomarkers for disease activity, compared to a time-agnostic regression model. Despite no significant difference in simple magnitude, using dynamic time warping analysis to align longitudinal profiles revealed a large subset ( n = 69) with significantly stronger associations between biological parameters and disease activity. This subgroup had significantly lower flare rates, disease activity and damage scores, suggesting this clustering is clinically meaningful. Conclusions: These results suggest associations between biological parameters and disease activity in SLE exist in a multi-dimensional time-dependent pattern, with implications for the analysis of biomarkers in SLE often used to identify therapeutic targets. Novel methods to analyse high-dimensional data and control for time-dependent variability may have broad utility in the study complex relationships between clinical and biological parameters.

Understanding the movement of animals is important for a wide range of scientific interests including migration, disease spread, collective movement behaviour and analysing motion in relation to dynamic changes of the environment such as wind and thermal lifts. Particularly, the three-dimensional (3D) spatial-temporal nature of bird movement data, which is widely available with high temporal and spatial resolution at large volumes, presents a natural option to explore the potential of immersive analytics (IA). We investigate the requirements and benefits of a wide range of immersive environments for explorative visualization and analytics of 3D movement data, in particular regarding design considerations for such 3D immersive environments, and present prototypes for IA solutions. Tailored to biologists studying bird movement data, the immersive solutions enable geo-locational time-series data to be investigated interactively, thus enabling experts to visually explore interesting angles of a flock and its behaviour in the context of the environment. The 3D virtual world presents the audience with engaging and interactive content, allowing users to 'fly with the flock', with the potential to ascertain an intuitive overview of often complex datasets, and to provide the opportunity thereby to formulate and at least qualitatively assess hypotheses. This work also contributes to ongoing research efforts to promote better understanding of bird migration and the associated environmental factors at the global scale, thereby providing a visual vehicle for driving public awareness of environmental issues and bird migration patterns.

Quantification of analytes such as cytokines in serum samples is intrinsic to translational research in immune diseases. Optimising pre-analytical conditions is critical for ensuring study quality, including evaluation of cytokine stability. We aimed to evaluate the effect on cytokine stability of storage duration prior to freezing of serum, and compare to plasma samples obtained from patients with systemic lupus erythematosus (SLE). Protein stability was analysed by simultaneously quantifying 18 analytes using a custom multi-analyte profile in SLE patient serum and plasma samples that had been prospectively stored at 4 °C for pre-determined periods between 0 and 30 days, prior to freezing. Six analytes were excluded from analysis, because most tested samples were above or below the limit of detection. Amongst the 12 analysed proteins, 11 did not show significant signal degradation. Significant signal degradation was observed from the fourth day of storage for a single analyte, CCL19. Proteins levels were more stable in unseparated serum compared to plasma for most analytes, with the exception of IL-37 which appears slightly more stable in plasma. Based on this, a maximum 3 days of storage at 4 °C for unseparated serum samples is recommended for biobanked samples intended for cytokine analysis in studies of human immune disease., (Copyright © 2018 Elsevier Ltd. All rights reserved.)

Spatially adjacent habitats on coral reefs can represent highly distinct environments, often harbouring different coral communities. Yet, certain coral species thrive across divergent environments. It is unknown whether the forces of selection are sufficiently strong to overcome the counteracting effects of the typically high gene flow over short distances, and for local adaptation to occur. We screened the coral genome (using restriction site-associated sequencing) and characterized both the dinoflagellate photosymbiont- and tissue-associated prokaryote microbiomes (using metabarcoding) of a reef flat and slope population of the reef-building coral, Pocillopora damicornis, at two locations on Heron Island in the southern Great Barrier Reef. Reef flat and slope populations were separated by <100 m horizontally and ~5 m vertically, and the two study locations were separated by ~1 km. For the coral host, genetic divergence between habitats was much greater than between locations, suggesting limited gene flow between the flat and slope populations. Consistent with environmental selection, outlier loci primarily belonged to the conserved, minimal cellular stress response, likely reflecting adaptation to the different temperature and irradiance regimes on the reef flat and slope. The prokaryote community differed across both habitat and, to a lesser extent, location, whereas the dinoflagellate photosymbionts differed by habitat but not location. The observed intraspecific diversity associated with divergent habitats supports that environmental adaptation involves multiple members of the coral holobiont. Adaptive alleles or microbial associations present in coral populations from the environmentally variable reef flat may provide a source of adaptive variation for assisted evolution approaches, through assisted gene flow, artificial cross-breeding or probiotic inoculations, with the aim to increase climate resilience in the slope populations., (© 2018 John Wiley & Sons Ltd.)

The growth and progression of solid tumors involves dynamic cross-talk between cancer epithelium and the surrounding microenvironment. To date, molecular profiling has largely been restricted to the epithelial component of tumors; therefore, features underpinning the persistent protumorigenic phenotype of the tumor microenvironment are unknown. Using whole-genome bisulfite sequencing, we show for the first time that cancer-associated fibroblasts (CAFs) from localized prostate cancer display remarkably distinct and enduring genome-wide changes in DNA methylation, significantly at enhancers and promoters, compared to nonmalignant prostate fibroblasts (NPFs). Differentially methylated regions associated with changes in gene expression have cancer-related functions and accurately distinguish CAFs from NPFs. Remarkably, a subset of changes is shared with prostate cancer epithelial cells, revealing the new concept of tumor-specific epigenome modifications in the tumor and its microenvironment. The distinct methylome of CAFs provides a novel epigenetic hallmark of the cancer microenvironment and promises new biomarkers to improve interpretation of diagnostic samples., (© 2018 Pidsley et al.; Published by Cold Spring Harbor Laboratory Press.)

Background: A strong focus of the post-genomic era is mining of the non-coding regulatory genome in order to unravel the function of regulatory elements that coordinate gene expression (Nat 489:57-74, 2012; Nat 507:462-70, 2014; Nat 507:455-61, 2014; Nat 518:317-30, 2015). Whole-genome approaches based on next-generation sequencing (NGS) have provided insight into the genomic location of regulatory elements throughout different cell types, organs and organisms. These technologies are now widespread and commonly used in laboratories from various fields of research. This highlights the need for fast and user-friendly software tools dedicated to extracting cis-regulatory information contained in these regulatory regions; for instance transcription factor binding site (TFBS) composition. Ideally, such tools should not require prior programming knowledge to ensure they are accessible for all users., Results: We present TrawlerWeb, a web-based version of the Trawler&#95;standalone tool (Nat Methods 4:563-5, 2007; Nat Protoc 5:323-34, 2010), to allow for the identification of enriched motifs in DNA sequences obtained from next-generation sequencing experiments in order to predict their TFBS composition. TrawlerWeb is designed for online queries with standard options common to web-based motif discovery tools. In addition, TrawlerWeb provides three unique new features: 1) TrawlerWeb allows the input of BED files directly generated from NGS experiments, 2) it automatically generates an input-matched biologically relevant background, and 3) it displays resulting conservation scores for each instance of the motif found in the input sequences, which assists the researcher in prioritising the motifs to validate experimentally. Finally, to date, this web-based version of Trawler&#95;standalone remains the fastest online de novo motif discovery tool compared to other popular web-based software, while generating predictions with high accuracy., Conclusions: TrawlerWeb provides users with a fast, simple and easy-to-use web interface for de novo motif discovery. This will assist in rapidly analysing NGS datasets that are now being routinely generated. TrawlerWeb is freely available and accessible at: http://trawler.erc.monash.edu.au .

Organ growth requires a careful balance between stem cell self-renewal and lineage commitment to ensure proper tissue expansion. The cellular and molecular mechanisms that mediate this balance are unresolved in most organs, including skeletal muscle. Here we identify a long-lived stem cell pool that mediates growth of the zebrafish myotome. This population exhibits extensive clonal drift, shifting from random deployment of stem cells during development to reliance on a small number of dominant clones to fuel the vast majority of muscle growth. This clonal drift requires Meox1, a homeobox protein that directly inhibits the cell-cycle checkpoint gene ccnb1. Meox1 initiates G 2 cell-cycle arrest within muscle stem cells, and disrupting this G 2 arrest causes premature lineage commitment and the resulting defects in muscle growth. These findings reveal that distinct regulatory mechanisms orchestrate stem cell dynamics during organ growth, beyond the G 0 /G 1 cell-cycle inhibition traditionally associated with maintaining tissue-resident stem cells., (Copyright © 2017. Published by Elsevier Inc.)

Mutations in the Nkx2-5 gene are a main cause of congenital heart disease. Several studies have addressed the phenotypic consequences of disrupting the Nkx2-5 gene locus, although animal models to date failed to recapitulate the full spectrum of the human disease. Here, we describe a new Nkx2-5 point mutation murine model, akin to its human counterpart disease-generating mutation. Our model fully reproduces the morphological and physiological clinical presentations of the disease and reveals an understudied aspect of Nkx2-5 -driven pathology, a primary right ventricular dysfunction. We further describe the molecular consequences of disrupting the transcriptional network regulated by Nkx2-5 in the heart and show that Nkx2-5 -dependent perturbation of the Wnt signaling pathway promotes heart dysfunction through alteration of cardiomyocyte metabolism. Our data provide mechanistic insights on how Nkx2-5 regulates heart function and metabolism, a link in the study of congenital heart disease, and confirms that our models are the first murine genetic models to our knowledge to present all spectra of clinically relevant adult congenital heart disease phenotypes generated by NKX2-5 mutations in patients., Competing Interests: Conflict of interest: The authors have declared that no conflict of interest exists.

Background: Quantitative high-throughput data deposited in consortia such as International Cancer Genome Consortium (ICGC) and The Cancer Genome Atlas (TCGA) present opportunities and challenges for computational analyses., Methods: We present a computational strategy to systematically rank and investigate a large number (2 10 -2 20 ) of clinically testable gene sets, using combinatorial gene subset generation and disease-free survival (DFS) analyses. This approach integrates protein-protein interaction networks, gene expression, DNA methylation, and copy number data, in association with DFS profiles from patient clinical records., Results: As a case study, we applied this pipeline to systematically analyze the role of ALDH1A2 in prostate cancer (PCa). We have previously found this gene to have multiple roles in disease and homeostasis, and here we investigate the role of the associated ALDH1A2 gene/protein networks in PCa, using our methodology in combination with PCa patient clinical profiles from ICGC and TCGA databases. Relationships between gene signatures and relapse were analyzed using Kaplan-Meier (KM) log-rank analysis and multivariable Cox regression. Relative expression versus pooled mean from diploid population was used for z -statistics calculation. Gene/protein interaction network analyses generated 11 core genes associated with ALDH1A2 ; combinatorial ranking of the power set of these core genes identified two gene sets (out of 2 11 - 1 = 2,047 combinations) with significant correlation with disease relapse (KM log rank p  < 0.05). For the more significant of these two sets, referred to as the optimal gene set (OGS), patients have median survival 62.7 months with OGS alterations compared to >150 months without OGS alterations ( p  = 0.0248, hazard ratio = 2.213, 95% confidence interval = 1.1-4.098). Two genes comprising OGS ( CYP26A1 and RDH10 ) are strongly associated with ALDH1A2 in the retinoic acid (RA) pathways, suggesting a major role of RA signaling in early PCa progression. Our pipeline complements human expertise in the search for prognostic biomarkers in large-scale datasets.

In the adult, tissue repair after injury is generally compromised by fibrosis, which maintains tissue integrity with scar formation but does not restore normal architecture and function. The process of regeneration is necessary to replace the scar and rebuild normal functioning tissue. Here, we address this problem in the context of heart disease, and discuss the origins and characteristics of cardiac fibroblasts, as well as the crucial role that they play in cardiac development and disease. We discuss the dual nature of cardiac fibroblasts, which can lead to scarring, pathological remodelling and functional deficit, but can also promote heart function in some contexts. Finally, we review current and proposed approaches whereby regeneration could be fostered by interventions that limit scar formation., (© 2016. Published by The Company of Biologists Ltd.)

The adult mammalian heart contains multiple cell types that work in unison under tightly regulated conditions to maintain homeostasis. Cardiac fibroblasts are a significant and unique population of non-muscle cells in the heart that have recently gained substantial interest in the cardiac biology community. To better understand this renaissance cell, it is essential to systematically survey what has been known in the literature about the cellular and molecular processes involved. We have built CARFMAP (http://visionet.erc.monash.edu.au/CARFMAP), an interactive cardiac fibroblast pathway map derived from the biomedical literature using a software-assisted manual data collection approach. CARFMAP is an information-rich interactive tool that enables cardiac biologists to explore the large body of literature in various creative ways. There is surprisingly little overlap between the cardiac fibroblast pathway map, a foreskin fibroblast pathway map, and a whole mouse organism signalling pathway map from the REACTOME database. Among the use cases of CARFMAP is a common task in our cardiac biology laboratory of identifying new genes that are (1) relevant to cardiac literature, and (2) differentially regulated in high-throughput assays. From the expression profiles of mouse cardiac and tail fibroblasts, we employed CARFMAP to characterise cardiac fibroblast pathways. Using CARFMAP in conjunction with transcriptomic data, we generated a stringent list of six genes that would not have been singled out using bioinformatics analyses alone. Experimental validation showed that five genes (Mmp3, Il6, Edn1, Pdgfc and Fgf10) are differentially regulated in the cardiac fibroblast. CARFMAP is a powerful tool for systems analyses of cardiac fibroblasts, facilitating systems-level cardiovascular research.

Background: Existing de novo software platforms have largely overlooked a valuable resource, the expertise of the intended biologist users. Typical data representations such as long gene lists, or highly dense and overlapping transcription factor networks often hinder biologists from relating these results to their expertise., Results: VISIONET, a streamlined visualisation tool built from experimental needs, enables biologists to transform large and dense overlapping transcription factor networks into sparse human-readable graphs via numerically filtering. The VISIONET interface allows users without a computing background to interactively explore and filter their data, and empowers them to apply their specialist knowledge on far more complex and substantial data sets than is currently possible. Applying VISIONET to the Tbx20-Gata4 transcription factor network led to the discovery and validation of Aldh1a2, an essential developmental gene associated with various important cardiac disorders, as a healthy adult cardiac fibroblast gene co-regulated by cardiogenic transcription factors Gata4 and Tbx20., Conclusions: We demonstrate with experimental validations the utility of VISIONET for expertise-driven gene discovery that opens new experimental directions that would not otherwise have been identified.

Understanding the complexity of cardiac physiology requires system-level studies of multiple cardiac cell types. Frequently, however, the end result of published research lacks the detail of the collaborative and integrative experimental design process, and the underlying conceptual framework. We review the recent progress in systems modelling and omics analysis of the heterocellular heart environment through complementary forward and inverse approaches, illustrating these conceptual and experimental frameworks with case studies from our own research program. The forward approach begins by collecting curated information from the niche cardiac biology literature, and connecting the dots to form mechanistic network models that generate testable system-level predictions. The inverse approach starts from the vast pool of public omics data in recent cardiac biological research, and applies bioinformatics analysis to produce novel candidates for further investigation. We also discuss the possibility of combining these two approaches into a hybrid framework, together with the benefits and challenges. These interdisciplinary research frameworks illustrate the interplay between computational models, omics analysis, and wet lab experiments, which holds the key to making real progress in improving human cardiac wellbeing., (Copyright © 2014 Elsevier Ltd. All rights reserved.)

Heart failure is one of the leading causes of death worldwide [1-4]. Current therapeutic strategies are inefficient and cannot cure this chronic and debilitating condition [5]. Ultimately, heart transplants are required for patient survival, but donor organs are scarce in availability and only prolong the life-span of patients for a limited time. Fibrosis is one of the main pathological features of heart failure [6,7], caused by inappropriate stimulation of fibroblasts and excessive extracellular matrix production. Therefore, an in-depth understanding of the cardiac fibroblast is essential to underpin effective therapeutic treatments for heart failure [5]. Fibroblasts in general have been an underappreciated cell type, regarded as relatively inert and providing only basic functionality; they are usually referred to as the  'biological glue' of all tissues in the body. However, more recent literature suggests that they actively participate in organ homeostasis and disease [7,8]. We have recently uncovered a unique molecular identity for fibroblasts isolated from the heart [9], expressing a set of cardiogenic transcription factors that have been previously associated with cardiomyocyte ontogenesis. This signature suggests that cardiac fibroblasts may be ideal for use in stem cell replacement therapies, as they may retain the memory of where they derive from embryologically. Our data also revealed that about 90% of fibroblasts from both tail and heart origins share a cell surface signature that has previously been described for mesenchymal stem cells (MSCs), raising the possibility that fibroblasts and MSCs may in fact be the same cell type. Thus, our findings carry profound implications for the field of regenerative medicine. Here, we describe detailed methodology and quality controls related to the gene expression profiling of cardiac fibroblasts, deposited at the Gene Expression Omnibus (GEO) under the accession number GSE50531. We also provide the R code to easily reproduce the data quantification and analysis processes.