Bronchial epithelium epithelial-mesenchymal plasticity forms aberrant basaloid-like cells in vitro.

Although epithelial-mesenchymal transition (EMT) is a common feature of fibrotic lung disease, its role in fibrogenesis is controversial. Recently, aberrant basaloid cells were identified in fibrotic lung tissue as a novel epithelial cell type displaying a partial EMT phenotype. The developmental origin of these cells remains unknown. To elucidate the role of EMT in the development of aberrant basaloid cells from the bronchial epithelium, we mapped EMT-induced transcriptional changes at the population and single-cell levels. Human bronchial epithelial cells grown as submerged or air-liquid interface (ALI) cultures with or without EMT induction were analyzed by bulk and single-cell RNA-Sequencing. Comparison of submerged and ALI cultures revealed differential expression of 8,247 protein coding (PC) and 1,621 long noncoding RNA (lncRNA) genes and revealed epithelial cell-type-specific lncRNAs. Similarly, EMT induction in ALI cultures resulted in robust transcriptional reprogramming of 6,020 PC and 907 lncRNA genes. Although there was no evidence for fibroblast/myofibroblast conversion following EMT induction, cells displayed a partial EMT gene signature and an aberrant basaloid-like cell phenotype. The substantial transcriptional differences between submerged and ALI cultures highlight that care must be taken when interpreting data from submerged cultures. This work supports that lung epithelial EMT does not generate fibroblasts/myofibroblasts and confirms ALI cultures provide a physiologically relevant system to study aberrant basaloid-like cells and mechanisms of EMT. We provide a catalog of PC and lncRNA genes and an interactive browser (https://bronc-epi-in-vitro.cells.ucsc.edu/) of single-cell RNA-Seq data for further exploration of potential roles in the lung epithelium in health and lung disease.

Effects of extended pharmacological disruption of zebrafish embryonic heart biomechanical environment on cardiac function, morphology, and gene expression.

BACKGROUND: Biomechanical stimuli are known to be important to cardiac development, but the mechanisms are not fully understood. Here, we pharmacologically disrupted the biomechanical environment of wild-type zebrafish embryonic hearts for an extended duration and investigated the consequent effects on cardiac function, morphological development, and gene expression. RESULTS: Myocardial contractility was significantly diminished or abolished in zebrafish embryonic hearts treated for 72 hours from 2 dpf with 2,3-butanedione monoxime (BDM). Image-based flow simulations showed that flow wall shear stresses were abolished or significantly reduced with high oscillatory shear indices. At 5 dpf, after removal of BDM, treated embryonic hearts were maldeveloped, having disrupted cardiac looping, smaller ventricles, and poor cardiac function (lower ejected flow, bulboventricular regurgitation, lower contractility, and slower heart rate). RNA sequencing of cardiomyocytes of treated hearts revealed 922 significantly up-regulated genes and 1,698 significantly down-regulated genes. RNA analysis and subsequent qPCR and histology validation suggested that biomechanical disruption led to an up-regulation of inflammatory and apoptotic genes and down-regulation of ECM remodeling and ECM-receptor interaction genes. Biomechanics disruption also prevented the formation of ventricular trabeculation along with notch1 and erbb4a down-regulation. CONCLUSIONS: Extended disruption of biomechanical stimuli caused maldevelopment, and potential genes responsible for this are identified.

A pragmatic randomised controlled trial of preferred intensity exercise in depressed adult women in the United Kingdom: secondary analysis of individual variability of depression.

BACKGROUND: This study is a secondary analysis of the trial by Callaghan et al. (2011), which reported higher antidepressant effects for preferred intensity (n = 19) vs. prescribed intensity (n = 19) exercise of three sessions/week over four weeks in depressed women. In particular, the present study sought to examine whether greater clinically significant individual change/recovery was observed in the preferred compared to the prescribed exercise group. METHODS: The reliable change index and the Ccutoff score criteria described by Jacobson and Truax (1991) were employed to determine clinical significance. These criteria examined if individual change in depression scores from pre- to post-intervention in the preferred intensity group were statistically significant beyond the standard error of difference derived from the active comparator prescribed group, and subsequently within a normal population range. Patients fulfilling the first or both criteria were classified as improved or recovered, respectively. RESULTS: Post-intervention depression scores of six patients in the preferred intensity exercise group (32%) demonstrated statistically reliable improvement (p < 0.05) and recovery. Half of this subgroup started as moderately depressed. No patient demonstrated a reliable deterioration in depression. Due to a small sample size, it was impossible to determine whether patients on psychiatric medication or medication-free patients were equally benefited from preferred intensity exercise. Thirteen patients in the preferred intensity group (68%) displayed non-statistically significant change in post-intervention depression scores (p > 0.05), although eight of them showed a non-significant improvement in post-intervention depression scores and three could not technically show an improvement in depression due to floor effects (baseline depression within normal range). CONCLUSIONS: Preferred intensity exercise of three sessions/week over four weeks led almost a third of the patients to record scores consistent with recovery from depression. Health professionals may consider that short-term preferred intensity exercise provides clinically significant antidepressant effects comparing favourably to exercise on prescription.

Redefinition of the human mast cell transcriptome by deep-CAGE sequencing.

Mast cells (MCs) mature exclusively in peripheral tissues, hampering research into their developmental and functional programs. Here, we employed deep cap analysis of gene expression on skin-derived MCs to generate the most comprehensive view of the human MC transcriptome ever reported. An advantage is that MCs were embedded in the FANTOM5 project, giving the opportunity to contrast their molecular signature against a multitude of human samples. We demonstrate that MCs possess a unique and surprising transcriptional landscape, combining hematopoietic genes with those exclusively active in MCs and genes not previously reported as expressed by MCs (several of them markers of unrelated tissues). We also found functional bone morphogenetic protein receptors transducing activatory signals in MCs. Conversely, several immune-related genes frequently studied in MCs were not expressed or were weakly expressed. Comparing MCs ex vivo with cultured counterparts revealed profound changes in the MC transcriptome in in vitro surroundings. We also determined the promoter usage of MC-expressed genes and identified associated motifs active in the lineage. Befitting their uniqueness, MCs had no close relative in the hematopoietic network (also only distantly related with basophils). This rich data set reveals that our knowledge of human MCs is still limited, but with this resource, novel functional programs of MCs may soon be discovered.

deltaGseg: macrostate estimation via molecular dynamics simulations and multiscale time series analysis.

SUMMARY: Binding free energy calculations obtained through molecular dynamics simulations reflect intermolecular interaction states through a series of independent snapshots. Typically, the free energies of multiple simulated series (each with slightly different starting conditions) need to be estimated. Previous approaches carry out this task by moving averages at certain decorrelation times, assuming that the system comes from a single conformation description of binding events. Here, we discuss a more general approach that uses statistical modeling, wavelets denoising and hierarchical clustering to estimate the significance of multiple statistically distinct subpopulations, reflecting potential macrostates of the system. We present the deltaGseg R package that performs macrostate estimation from multiple replicated series and allows molecular biologists/chemists to gain physical insight into the molecular details that are not easily accessible by experimental techniques. AVAILABILITY: deltaGseg is a Bioconductor R package available at http://bioconductor.org/packages/release/bioc/html/deltaGseg.html.

Computationally defined and in vitro validated putative genomic safe harbour loci for transgene expression in human cells.

Selection of the target site is an inherent question for any project aiming for directed transgene integration. Genomic safe harbour (GSH) loci have been proposed as safe sites in the human genome for transgene integration. Although several sites have been characterised for transgene integration in the literature, most of these do not meet criteria set out for a GSH and the limited set that do have not been characterised extensively. Here, we conducted a computational analysis using publicly available data to identify 25 unique putative GSH loci that reside in active chromosomal compartments. We validated stable transgene expression and minimal disruption of the native transcriptome in three GSH sites in vitro using human embryonic stem cells (hESCs) and their differentiated progeny. Furthermore, for easy targeted transgene expression, we have engineered constitutive landing pad expression constructs into the three validated GSH in hESCs.

Myocardial Biomechanics and the Consequent Differentially Expressed Genes of the Left Atrial Ligation Chick Embryonic Model of Hypoplastic Left Heart Syndrome.

Left atrial ligation (LAL) of the chick embryonic heart is a model of the hypoplastic left heart syndrome (HLHS) where a purely mechanical intervention without genetic or pharmacological manipulation is employed to initiate cardiac malformation. It is thus a key model for understanding the biomechanical origins of HLHS. However, its myocardial mechanics and subsequent gene expressions are not well-understood. We performed finite element (FE) modeling and single-cell RNA sequencing to address this. 4D high-frequency ultrasound imaging of chick embryonic hearts at HH25 (ED 4.5) were obtained for both LAL and control. Motion tracking was performed to quantify strains. Image-based FE modeling was conducted, using the direction of the smallest strain eigenvector as the orientations of contractions, the Guccione active tension model and a Fung-type transversely isotropic passive stiffness model that was determined via micro-pipette aspiration. Single-cell RNA sequencing of left ventricle (LV) heart tissues was performed for normal and LAL embryos at HH30 (ED 6.5) and differentially expressed genes (DEG) were identified.After LAL, LV thickness increased by 33%, strains in the myofiber direction increased by 42%, while stresses in the myofiber direction decreased by 50%. These were likely related to the reduction in ventricular preload and underloading of the LV due to LAL. RNA-seq data revealed potentially related DEG in myocytes, including mechano-sensing genes (Cadherins, NOTCH1, etc.), myosin contractility genes (MLCK, MLCP, etc.), calcium signaling genes (PI3K, PMCA, etc.), and genes related to fibrosis and fibroelastosis (TGF-ß, BMP, etc.). We elucidated the changes to the myocardial biomechanics brought by LAL and the corresponding changes to myocyte gene expressions. These data may be useful in identifying the mechanobiological pathways of HLHS.

C1 CAGE detects transcription start sites and enhancer activity at single-cell resolution.

Single-cell transcriptomic profiling is a powerful tool to explore cellular heterogeneity. However, most of these methods focus on the 3'-end of polyadenylated transcripts and provide only a partial view of the transcriptome. We introduce C1 CAGE, a method for the detection of transcript 5'-ends with an original sample multiplexing strategy in the C1TM microfluidic system. We first quantifiy the performance of C1 CAGE and find it as accurate and sensitive as other methods in the C1 system. We then use it to profile promoter and enhancer activities in the cellular response to TGF-ß of lung cancer cells and discover subpopulations of cells differing in their response. We also describe enhancer RNA dynamics revealing transcriptional bursts in subsets of cells with transcripts arising from either strand in a mutually exclusive manner, validated using single molecule fluorescence in situ hybridization.

Sense-antisense gene-pairs in breast cancer and associated pathological pathways.

More than 30% of human protein-coding genes form hereditary complex genome architectures composed of sense-antisense (SA) gene pairs (SAGPs) transcribing their RNAs from both strands of a given locus. Such architectures represent important novel components of genome complexity contributing to gene expression deregulation in cancer cells. Therefore, the architectures might be involved in cancer pathways and, in turn, be used for novel drug targets discovery. However, the global roles of SAGPs in cancer pathways has not been studied. Here we investigated SAGPs associated with breast cancer (BC)-related pathways using systems biology, prognostic survival and experimental methods. Gene expression analysis identified 73 BC-relevant SAGPs that are highly correlated in BC. Survival modelling and metadata analysis of the 1161 BC patients allowed us to develop a novel patient prognostic grouping method selecting the 12 survival-significant SAGPs. The qRT-PCR-validated 12-SAGP prognostic signature reproducibly stratified BC patients into low- and high-risk prognostic subgroups. The 1381 SAGP-defined differentially expressed genes common across three studied cohorts were identified. The functional enrichment analysis of these genes revealed the GABPA gene network, including BC-relevant SAGPs, specific gene sets involved in cell cycle, spliceosomal and proteasomal pathways. The co-regulatory function of GABPA in BC cells was supported using siRNA knockdown studies. Thus, we demonstrated SAGPs as the synergistically functional genome architectures interconnected with cancer-related pathways and associated with BC patient clinical outcomes. Taken together, SAGPs represent an important component of genome complexity which can be used to identify novel aspects of coordinated pathological gene networks in cancers.

Retinoic acid potentiates inflammatory cytokines in human mast cells: identification of mast cells as prominent constituents of the skin retinoid network.

Retinoic acid (RA), the active vitamin-A-metabolite, has well-established functions in skin homeostasis and in the immune system. Skin mast cells (MCs) combine traits of both structures, being of hematopoietic origin, but functional in the skin environment. It remains largely unknown whether mature MCs are targeted by the retinoid network. Here, we demonstrate that human skin MCs display substantial susceptibility to RA by which they are instructed to increase pro-inflammatory mediators (IL-1ß, IL-8, TNF-α) but not histamine release. The effects are observed at physiological RA levels, in different microenvironments, and are largely donor-independent. RA susceptibility is owed to the cells' abundant expression of RARA, the receptor mediating MC cytokine responses. Unexpectedly, bioinformatics calculations on the FANTOM5 expression atlas revealed general enrichment of retinoid network components in MCs against other skin cells, and MCs rapidly upregulated RA responsive genes. In conclusion, MCs are important yet hitherto overlooked retinoid targets in the skin.

Mast cell transcriptome elucidation: what are the implications for allergic disease in the clinic and where do we go next?

Mast cells (MCs) are principal effector cells of type-I-allergic reactions but still poorly defined in humans. The consortium Functional Annotation of the Mammalian Genome 5 has created a map of body-wide transcriptome data for a multitude of human cell types, including MCs. MCs were found to have a surprising transcriptional landscape expressing a range of genes not (or barely active) elsewhere in the body. Whereas several MC specific genes have no annotated function, others belong to networks defining specific MC traits and functions, such as granule architecture, IgER signaling, exocytosis and mediator production. Several of these genes are so highly enriched in MCs (versus all other cells), that they appear potentially specific targets for therapeutic interventions in diseases in which MCs are actively involved. We present some interesting candidates, highlight the uniqueness of MCs and discuss their role in allergy and itch sensation based on these renewed insights.