Migration and differentiation of muscle stem cells are coupled by RhoA signalling during regeneration.

Skeletal muscle is highly regenerative and is mediated by a population of migratory adult muscle stem cells (muSCs). Effective muscle regeneration requires a spatio-temporally regulated response of the muSC population to generate sufficient muscle progenitor cells that then differentiate at the appropriate time. The relationship between muSC migration and cell fate is poorly understood and it is not clear how forces experienced by migrating cells affect cell behaviour. We have used zebrafish to understand the relationship between muSC cell adhesion, behaviour and fate in vivo. Imaging of pax7-expressing muSCs as they respond to focal injuries in trunk muscle reveals that they migrate by protrusive-based means. By carefully characterizing their behaviour in response to injury we find that they employ an adhesion-dependent mode of migration that is regulated by the RhoA kinase ROCK. Impaired ROCK activity results in reduced expression of cell cycle genes and increased differentiation in regenerating muscle. This correlates with changes to focal adhesion dynamics and migration, revealing that ROCK inhibition alters the interaction of muSCs to their local environment. We propose that muSC migration and differentiation are coupled processes that respond to changes in force from the environment mediated by RhoA signalling.

Endocardium-to-coronary artery differentiation during heart development and regeneration involves sequential roles of Bmp2 and Cxcl12/Cxcr4.

Endocardial cells lining the heart lumen are coronary vessel progenitors during embryogenesis. Re-igniting this developmental process in adults could regenerate blood vessels lost during cardiac injury, but this requires additional knowledge of molecular mechanisms. Here, we use mouse genetics and scRNA-seq to identify regulators of endocardial angiogenesis and precisely assess the role of CXCL12/CXCR4 signaling. Time-specific lineage tracing demonstrated that endocardial cells differentiated into coronary endothelial cells primarily at mid-gestation. A new mouse line reporting CXCR4 activity-along with cell-specific gene deletions-demonstrated it was specifically required for artery morphogenesis rather than angiogenesis. Integrating scRNA-seq data of endocardial-derived coronary vessels from mid- and late-gestation identified a Bmp2-expressing transitioning population specific to mid-gestation. Bmp2 stimulated endocardial angiogenesis in vitro and in injured neonatal mouse hearts. Our data demonstrate how understanding the molecular mechanisms underlying endocardial angiogenesis can identify new potential therapeutic targets promoting revascularization of the injured heart.

Quantitative lineage analysis identifies a hepato-pancreato-biliary progenitor niche.

Studies based on single cells have revealed vast cellular heterogeneity in stem cell and progenitor compartments, suggesting continuous differentiation trajectories with intermixing of cells at various states of lineage commitment and notable degrees of plasticity during organogenesis1-5. The hepato-pancreato-biliary organ system relies on a small endoderm progenitor compartment that gives rise to a variety of different adult tissues, including the liver, pancreas, gall bladder and extra-hepatic bile ducts6,7. Experimental manipulation of various developmental signals in the mouse embryo has underscored important cellular plasticity in this embryonic territory6. This is reflected in the existence of human genetic syndromes as well as congenital malformations featuring multi-organ phenotypes in liver, pancreas and gall bladder6. Nevertheless, the precise lineage hierarchy and succession of events leading to the segregation of an endoderm progenitor compartment into hepatic, biliary and pancreatic structures have not yet been established. Here we combine computational modelling approaches with genetic lineage tracing to accurately reconstruct the hepato-pancreato-biliary lineage tree. We show that a multipotent progenitor subpopulation persists in the pancreato-biliary organ rudiment, contributing cells not only to the pancreas and gall bladder but also to the liver. Moreover, using single-cell RNA sequencing and functional experiments we define a specialized niche that supports this subpopulation in a multipotent state for an extended time during development. Together these findings indicate sustained plasticity underlying hepato-pancreato-biliary development that might also explain the rapid expansion of the liver while attenuating pancreato-biliary growth.

Genomic landscape and clonal architecture of mouse oral squamous cell carcinomas dictate tumour ecology.

To establish whether 4-nitroquinoline N-oxide-induced carcinogenesis mirrors the heterogeneity of human oral squamous cell carcinoma (OSCC), we have performed genomic analysis of mouse tongue lesions. The mutational signatures of human and mouse OSCC overlap extensively. Mutational burden is higher in moderate dysplasias and invasive SCCs than in hyperplasias and mild dysplasias, although mutations in p53, Notch1 and Fat1 occur in early lesions. Laminin-α3 mutations are associated with tumour invasiveness and Notch1 mutant tumours have an increased immune infiltrate. Computational modelling of clonal dynamics indicates that high genetic heterogeneity may be a feature of those mild dysplasias that are likely to progress to more aggressive tumours. These studies provide a foundation for exploring OSCC evolution, heterogeneity and progression.

Cell Tracking Profiler - a user-driven analysis framework for evaluating 4D live-cell imaging data.

Accurate measurements of cell morphology and behaviour are fundamentally important for understanding how disease, molecules and drugs affect cell function in vivo Here, by using muscle stem cell (muSC) responses to injury in zebrafish as our biological paradigm, we established a 'ground truth' for muSC behaviour. This revealed that segmentation and tracking algorithms from commonly used programs are error-prone, leading us to develop a fast semi-automated image analysis pipeline that allows user-defined parameters for segmentation and correction of cell tracking. Cell Tracking Profiler (CTP) is a package that runs two existing programs, HK Means and Phagosight within the Icy image analysis suite, to enable user-managed cell tracking from 3D time-lapse datasets to provide measures of cell shape and movement. We demonstrate how CTP can be used to reveal changes to cell behaviour of muSCs in response to manipulation of the cell cytoskeleton by small-molecule inhibitors. CTP and the associated tools we have developed for analysis of outputs thus provide a powerful framework for analysing complex cell behaviour in vivo from 4D datasets that are not amenable to straightforward analysis.

Single cell analysis of human foetal liver captures the transcriptional profile of hepatobiliary hybrid progenitors.

The liver parenchyma is composed of hepatocytes and bile duct epithelial cells (BECs). Controversy exists regarding the cellular origin of human liver parenchymal tissue generation during embryonic development, homeostasis or repair. Here we report the existence of a hepatobiliary hybrid progenitor (HHyP) population in human foetal liver using single-cell RNA sequencing. HHyPs are anatomically restricted to the ductal plate of foetal liver and maintain a transcriptional profile distinct from foetal hepatocytes, mature hepatocytes and mature BECs. In addition, molecular heterogeneity within the EpCAM+ population of freshly isolated foetal and adult human liver identifies diverse gene expression signatures of hepatic and biliary lineage potential. Finally, we FACS isolate foetal HHyPs and confirm their hybrid progenitor phenotype in vivo. Our study suggests that hepatobiliary progenitor cells previously identified in mice also exist in humans, and can be distinguished from other parenchymal populations, including mature BECs, by distinct gene expression profiles.

Identifying Extrinsic versus Intrinsic Drivers of Variation in Cell Behavior in Human iPSC Lines from Healthy Donors.

Large cohorts of human induced pluripotent stem cells (iPSCs) from healthy donors are a potentially powerful tool for investigating the relationship between genetic variants and cellular behavior. Here, we integrate high content imaging of cell shape, proliferation, and other phenotypes with gene expression and DNA sequence datasets from over 100 human iPSC lines. By applying a dimensionality reduction approach, Probabilistic Estimation of Expression Residuals (PEER), we extracted factors that captured the effects of intrinsic (genetic concordance between different cell lines from the same donor) and extrinsic (cell responses to different fibronectin concentrations) conditions. We identify genes that correlate in expression with intrinsic and extrinsic PEER factors and associate outlier cell behavior with genes containing rare deleterious non-synonymous SNVs. Our study, thus, establishes a strategy for examining the genetic basis of inter-individual variability in cell behavior.

pATsi: Paralogs and Singleton Genes from Arabidopsis thaliana.

Arabidopsis thaliana is widely accepted as a model species in plant biology. Its genome, due to its small size and diploidy, was the first to be sequenced among plants, making this species also a reference for plant comparative genomics. Nevertheless, the evolutionary mechanisms that shaped the Arabidopsis genome are still controversial. Indeed, duplications, translocations, inversions, and gene loss events that contributed to the current organization are difficult to be traced. A reliable identification of paralogs and single-copy genes is essential to understand these mechanisms. Therefore, we implemented a dedicated pipeline to identify paralog genes and classify single-copy genes into opportune categories. PATsi, a web-accessible database, was organized to allow the straightforward access to the paralogs organized into networks and to the classification of single-copy genes. This permits to efficiently explore the gene collection of Arabidopsis for evolutionary investigations and comparative genomics.

Intergenic Alu exonisation facilitates the evolution of tissue-specific transcript ends.

The 3' untranslated regions (3' UTRs) of transcripts serve as important hubs for posttranscriptional gene expression regulation. Here, we find that the exonisation of intergenic Alu elements introduced new terminal exons and polyadenylation sites during human genome evolution. While Alu exonisation from introns has been described previously, we shed light on a novel mechanism to create alternative 3' UTRs, thereby opening opportunities for differential posttranscriptional regulation. On the mechanistic level, we show that intergenic Alu exonisation can compete both with alternative splicing and polyadenylation in the upstream gene. Notably, the Alu-derived isoforms are often expressed in a tissue-specific manner, and the Alu-derived 3' UTRs can alter mRNA stability. In summary, we demonstrate that intergenic elements can affect processing of preceding genes, and elucidate how intergenic Alu exonisation can contribute to tissue-specific posttranscriptional regulation by expanding the repertoire of 3' UTRs.

hiCLIP reveals the in vivo atlas of mRNA secondary structures recognized by Staufen 1.

The structure of messenger RNA is important for post-transcriptional regulation, mainly because it affects binding of trans-acting factors. However, little is known about the in vivo structure of full-length mRNAs. Here we present hiCLIP, a biochemical technique for transcriptome-wide identification of RNA secondary structures interacting with RNA-binding proteins (RBPs). Using this technique to investigate RNA structures bound by Staufen 1 (STAU1) in human cells, we uncover a dominance of intra-molecular RNA duplexes, a depletion of duplexes from coding regions of highly translated mRNAs, an unexpected prevalence of long-range duplexes in 3' untranslated regions (UTRs), and a decreased incidence of single nucleotide polymorphisms in duplex-forming regions. We also discover a duplex spanning 858 nucleotides in the 3' UTR of the X-box binding protein 1 (XBP1) mRNA that regulates its cytoplasmic splicing and stability. Our study reveals the fundamental role of mRNA secondary structures in gene expression and introduces hiCLIP as a widely applicable method for discovering new, especially long-range, RNA duplexes.

Exploiting a Reference Genome in Terms of Duplications: The Network of Paralogs and Single Copy Genes in Arabidopsis thaliana.

Arabidopsis thaliana became the model organism for plant studies because of its small diploid genome, rapid lifecycle and short adult size. Its genome was the first among plants to be sequenced, becoming the reference in plant genomics. However, the Arabidopsis genome is characterized by an inherently complex organization, since it has undergone ancient whole genome duplications, followed by gene reduction, diploidization events and extended rearrangements, which relocated and split up the retained portions. These events, together with probable chromosome reductions, dramatically increased the genome complexity, limiting its role as a reference. The identification of paralogs and single copy genes within a highly duplicated genome is a prerequisite to understand its organization and evolution and to improve its exploitation in comparative genomics. This is still controversial, even in the widely studied Arabidopsis genome. This is also due to the lack of a reference bioinformatics pipeline that could exhaustively identify paralogs and singleton genes. We describe here a complete computational strategy to detect both duplicated and single copy genes in a genome, discussing all the methodological issues that may strongly affect the results, their quality and their reliability. This approach was used to analyze the organization of Arabidopsis nuclear protein coding genes, and besides classifying computationally defined paralogs into networks and single copy genes into different classes, it unraveled further intriguing aspects concerning the genome annotation and the gene relationships in this reference plant species. Since our results may be useful for comparative genomics and genome functional analyses, we organized a dedicated web interface to make them accessible to the scientific community.

Quantitative expression profiling of highly degraded RNA from formalin-fixed, paraffin-embedded breast tumor biopsies by oligonucleotide microarrays.

Microarray-based gene expression profiling is well suited for parallel quantitative analysis of large numbers of RNAs, but its application to cancer biopsies, particularly formalin-fixed, paraffin-embedded (FFPE) archived tissues, is limited by the poor quality of the RNA recovered. This represents a serious drawback, as FFPE tumor tissue banks are available with clinical and prognostic annotations, which could be exploited for molecular profiling studies, provided that reliable analytical technologies are found. We applied and evaluated here a microarray-based cDNA-mediated annealing, selection, extension and ligation (DASL) assay for analysis of 502 mRNAs in highly degraded total RNA extracted from cultured cells or FFPE breast cancer (MT) biopsies. The study included quantitative and qualitative comparison of data obtained by analysis of the same RNAs with genome-wide oligonucleotide microarrays vs DASL arrays and, by DASL, before and after extensive in vitro RNA fragmentation. The DASL-based expression profiling assay applied to RNA extracted from MCF-7 cells, before or after 24 h stimulation with a mitogenic dose of 17beta-estradiol, consistently allowed to detect hormone-induced gene expression changes following extensive RNA degradation in vitro. Comparable results where obtained with tumor RNA extracted from FFPE MT biopsies (6 to 19 years old). The method proved itself sensitive, reproducible and accurate, when compared to results obtained by microarray analysis of RNA extracted from snap-frozen tissue of the same tumor.