AbstractmRNA structure is important for post-transcriptional regulation, largely because it affects binding of trans-acting factors 1 . However, little is known about the in vivo structure of full-length mRNAs. Here we present hiCLIP, a high-throughput technique to identify RNA secondary structures interacting with RNA-binding proteins (RBPs) in vivo. Using this technique to investigate RNA structures bound by Staufen 1 (STAU1), we uncover a dominance of intra-molecular RNA duplexes, a depletion of duplexes from coding regions of highly translated mRNAs, an unforeseen prevalence of long-range duplexes in 3′ untranslated regions (UTRs), and a decreased incidence of SNPs in duplex-forming regions. We also discover a duplex spanning 858nts 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 regulation and introduces hiCLIP as a widely applicable method for discovering novel, especially long-range, RNA duplexes.

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.

Summary 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.