Linking promoters to functional transcripts in small samples with nanoCAGE and CAGEscan

Mechanobiological stimulation in organ-on-a-chip systems reduces hepatic drug metabolic capacity in favor of regenerative specialization

Hepatic physiology depends on the liver's complex structural composition which among others, provides high oxygen supply rates, locally differential oxygen tension, endothelial paracrine signaling, as well as residual hemodynamic shear stress to resident hepatocytes. While functional improvements were shown by implementing these factors into hepatic culture systems, direct cause-effect relationships are often not well characterized-obfuscating their individual contribution in more complex microphysiological systems. By comparing increasingly complex hepatic in vitro culture systems that gradually implement these parameters, we investigate the influence of the cellular microenvironment to overall hepatic functionality in pharmacological applications. Here, hepatocytes were modulated in terms of oxygen tension and supplementation, endothelial coculture, and exposure to fluid shear stress delineated from oxygen influx. Results from transcriptomic and metabolomic evaluation indicate that particularly oxygen supply rates are critical to enhance cellular functionality-with cellular drug metabolism remaining comparable to physiological conditions after prolonged static culture. Endothelial signaling was found to be a major contributor to differential phenotype formation known as metabolic zonation, indicated by WNT pathway activity. Lastly, oxygen-delineated shear stress was identified to direct cellular fate towards increased hepatic plasticity and regenerative phenotypes at the cost of drug metabolic functionality - in line with regenerative effects observed in vivo. With these results, we provide a systematic evaluation of critical parameters and their impact in hepatic systems. Given their adherence to physiological effects in vivo, this highlights the importance of their implementation in biomimetic devices, such as organ-on-a-chip systems. Considering recent advances in basic liver biology, direct translation of physiological structures into in vitro models is a promising strategy to expand the capabilities of pharmacological models.

Dynamic, IPSC-derived hepatic tissue tri-culture system for the evaluation of liver physiologyin vitro

Availability of hepatic tissue for the investigation of metabolic processes is severely limited. While primary hepatocytes or animal models are widely used in pharmacological applications, a change in methodology towards more sustainable and ethical assays is highly desirable. Stem cell derived hepatic cells are generally regarded as a viable alternative for the above model systems, if current limitations in functionality and maturation can be overcome. By combining microfluidic organ-on-a-chip technology with individually differentiated, multicellular hepatic tissue fractions, we aim to improve overall functionality of hepatocyte-like cells, as well as evaluate cellular composition and interactions with non-parenchymal cell populations towards the formation of mature liver tissue. Utilizing a multi-omic approach, we show the improved maturation profiles of hepatocyte-like cells maintained in a dynamic microenvironment compared to standard tissue culture setups without continuous perfusion. In order to evaluate the resulting tissue, we employ single cell sequencing to distinguish formed subpopulations and spatial localization. While cellular input was strictly defined based on established differentiation protocols of parenchyma, endothelial and stellate cell fractions, resulting hepatic tissue was shown to comprise a complex mixture of epithelial and non-parenchymal fractions with specific local enrichment of phenotypes along the microchannel. Following this approach, we show the importance of passive, paracrine developmental processes in tissue formation. Using such complex tissue models is a crucial first step to develop stem cell-derivedin vitrosystems that can compare functionally with currently used pharmacological and toxicological applications.

Nano-CUT&Tag for multimodal chromatin profiling at single-cell resolution

The ability to comprehensively analyze the chromatin state with single-cell resolution is crucial for understanding gene regulatory principles in heterogenous tissues or during development. Recently, we developed a nanobody-based single-cell CUT&Tag (nano-CT) protocol to simultaneously profile three epigenetic modalities-two histone marks and open chromatin state-from the same single cell. Nano-CT implements a new set of secondary nanobody-Tn5 fusion proteins to direct barcoded tagmentation by Tn5 transposase to genomic targets labeled by primary antibodies raised in different species. Such nanobody-Tn5 fusion proteins are currently not commercially available, and their in-house production and purification can be completed in 3-4 d by following our detailed protocol. The single-cell indexing in nano-CT is performed on a commercially available platform, making it widely accessible to the community. In comparison to other multimodal methods, nano-CT stands out in data complexity, low sample requirements and the flexibility to choose two of the three modalities. In addition, nano-CT works efficiently with fresh brain samples, generating multimodal epigenomic profiles for thousands of brain cells at single-cell resolution. The nano-CT protocol can be completed in just 3 d by users with basic skills in standard molecular biology and bioinformatics, although previous experience with single-cell assay for transposase-accessible chromatin using sequencing (scATAC-seq) is beneficial for more in-depth data analysis. As a multimodal assay, nano-CT holds immense potential to reveal interactions of various chromatin modalities, to explore epigenetic heterogeneity and to increase our understanding of the role and interplay that chromatin dynamics has in cellular development.

Multiplex generation and single cell analysis of structural variants in a mammalian genome

The functional consequences of structural variants (SVs) in mammalian genomes are challenging to study. This is due to several factors, including: 1) their numerical paucity relative to other forms of standing genetic variation such as single nucleotide variants (SNVs) and short insertions or deletions (indels); 2) the fact that a single SV can involve and potentially impact the function of more than one gene and/or cis regulatory element; and 3) the relative immaturity of methods to generate and map SVs, either randomly or in targeted fashion, in in vitro or in vivo model systems. Towards addressing these challenges, we developed Genome-Shuffle-seq, a straightforward method that enables the multiplex generation and mapping of several major forms of SVs (deletions, inversions, translocations) throughout a mammalian genome. Genome-Shuffle-seq is based on the integration of "shuffle cassettes" to the genome, wherein each shuffle cassette contains components that facilitate its site-specific recombination (SSR) with other integrated shuffle cassettes (via Cre-loxP), its mapping to a specific genomic location (via T7-mediated in vitro transcription or IVT), and its identification in single-cell RNA-seq (scRNA-seq) data (via T7-mediated in situ transcription or IST). In this proof-of-concept, we apply Genome-Shuffle-seq to induce and map thousands of genomic SVs in mouse embryonic stem cells (mESCs) in a single experiment. Induced SVs are rapidly depleted from the cellular population over time, possibly due to Cre-mediated toxicity and/or negative selection on the rearrangements themselves. Leveraging T7 IST of barcodes whose positions are already mapped, we further demonstrate that we can efficiently genotype which SVs are present in association with each of many single cell transcriptomes in scRNA-seq data. Finally, preliminary evidence suggests our method may be a powerful means of generating extrachromosomal circular DNAs (ecDNAs). Looking forward, we anticipate that Genome-Shuffle-seq may be broadly useful for the systematic exploration of the functional consequences of SVs on gene expression, the chromatin landscape, and 3D nuclear architecture. We further anticipate potential uses for in vitro modeling of ecDNAs, as well as in paving the path to a minimal mammalian genome.

An improved method for the highly specific detection of transcription start sites

Precise detection of the transcriptional start site (TSS) is a key for characterizing transcriptional regulation of genes and for annotation of newly sequenced genomes. Here, we describe the development of an improved method, designated 'TSS-seq2.' This method is an iterative improvement of TSS-seq, a previously published enzymatic cap-structure conversion method to detect TSSs in base sequences. By modifying the original procedure, including by introducing split ligation at the key cap-selection step, the yield and the accuracy of the reaction has been substantially improved. For example, TSS-seq2 can be conducted using as little as 5 ng of total RNA with an overall accuracy of 96%; this yield a less-biased and more precise detection of TSS. We then applied TSS-seq2 for TSS analysis of four plant species that had not yet been analyzed by any previous TSS method.

Using long-read CAGE sequencing to profile cryptic-promoter-derived transcripts and their contribution to the immunopeptidome

Recent studies have shown that the noncoding genome can produce unannotated proteins as antigens that induce immune response. One major source of this activity is the aberrant epigenetic reactivation of transposable elements (TEs). In tumors, TEs often provide cryptic or alternate promoters, which can generate transcripts that encode tumor-specific unannotated proteins. Thus, TE-derived transcripts (TE transcripts) have the potential to produce tumor-specific, but recurrent, antigens shared among many tumors. Identification of TE-derived tumor antigens holds the promise to improve cancer immunotherapy approaches; however, current genomics and computational tools are not optimized for their detection. Here we combined CAGE technology with full-length long-read transcriptome sequencing (long-read CAGE, or LRCAGE) and developed a suite of computational tools to significantly improve immunopeptidome detection by incorporating TE and other tumor transcripts into the proteome database. By applying our methods to human lung cancer cell line H1299 data, we show that long-read technology significantly improves mapping of promoters with low mappability scores and that LRCAGE guarantees accurate construction of uncharacterized 5' transcript structure. Augmenting a reference proteome database with newly characterized transcripts enabled us to detect noncanonical antigens from HLA-pulldown LC-MS/MS data. Lastly, we show that epigenetic treatment increased the number of noncanonical antigens, particularly those encoded by TE transcripts, which might expand the pool of targetable antigens for cancers with low mutational burden.

Deep annotation of long noncoding RNAs by assembling RNA-seq and small RNA-seq data

Long noncoding RNAs (lncRNAs) are increasingly being recognized as modulators in various biological processes. However, due to their low expression, their systematic characterization is difficult to determine. Here, we performed transcript annotation by a newly developed computational pipeline, termed RNA-seq and small RNA-seq combined strategy (RSCS), in a wide variety of cellular contexts. Thousands of high-confidence potential novel transcripts were identified by the RSCS, and the reliability of the transcriptome was verified by analysis of transcript structure, base composition, and sequence complexity. Evidenced by the length comparison, the frequency of the core promoter and the polyadenylation signal motifs, and the locations of transcription start and end sites, the transcripts appear to be full length. Furthermore, taking advantage of our strategy, we identified a large number of endogenous retrovirus-associated lncRNAs, and a novel endogenous retrovirus-lncRNA that was functionally involved in control of Yap1 expression and essential for early embryogenesis was identified. In summary, the RSCS can generate a more complete and precise transcriptome, and our findings greatly expanded the transcriptome annotation for the mammalian community.

Implications of tumour heterogeneity on cancer evolution and therapy resistance: lessons from breast cancer

Tumour heterogeneity is pervasive amongst many cancers and leads to disease progression, and therapy resistance. In this review, using breast cancer as an exemplar, we focus on the recent advances in understanding the interplay between tumour cells and their microenvironment using single cell sequencing and digital spatial profiling technologies. Further, we discuss the utility of lineage tracing methodologies in pre-clinical models of breast cancer, and how these are being used to unravel new therapeutic vulnerabilities and reveal biomarkers of breast cancer progression. © 2023 The Pathological Society of Great Britain and Ireland.

ZBTB12 is a molecular barrier to dedifferentiation in human pluripotent stem cells

Development is generally viewed as one-way traffic of cell state transition from primitive to developmentally advanced states. However, molecular mechanisms that ensure the unidirectional transition of cell fates remain largely unknown. Through exact transcription start site mapping, we report an evolutionarily conserved BTB domain-containing zinc finger protein, ZBTB12, as a molecular barrier for dedifferentiation of human pluripotent stem cells (hPSCs). Single-cell RNA sequencing reveals that ZBTB12 is essential for three germ layer differentiation by blocking hPSC dedifferentiation. Mechanistically, ZBTB12 fine-tunes the expression of human endogenous retrovirus H (HERVH), a primate-specific retrotransposon, and targets specific transcripts that utilize HERVH as a regulatory element. In particular, the downregulation of HERVH-overlapping long non-coding RNAs (lncRNAs) by ZBTB12 is necessary for a successful exit from a pluripotent state and lineage derivation. Overall, we identify ZBTB12 as a molecular barrier that safeguards the unidirectional transition of metastable stem cell fates toward developmentally advanced states.

Functional evaluation of LTR-derived lncRNAs in porcine oocytes and zygotes with RNA-seq and small RNA-seq

Long noncoding RNAs (lncRNAs) are increasingly being recognized as modulators of early embryonic development in mammals. However, they are seldom investigated in pigs. Here, to annotate full-length RNA transcripts, we performed annotation using a newly developed computational pipeline—an RNA-seq and small RNA-seq combined strategy—using our previously obtained RNA-seq and small RNA-seq data from porcine oocytes and zygotes. As evidenced by the length comparison, the frequency of the core promoter, and the polyadenylation signal motifs, the transcripts appear to be full-length. Furthermore, our strategy allowed the identification of a large number of endogenous retrovirus-associated lncRNAs (ERV-lncRNAs) and found that some of them were highly expressed in porcine zygotes, as compared to oocytes. Through the knockdown strategy, two ERV-lncRNAs (TCONS_00035465 and TCONS_00031520) were identified as playing potential roles in the early embryo development of pigs, laying a foundation for future research.

Complete Transcriptome Analysis by 5'-End Single-Cell RNA-Seq with Random Priming

Single-cell transcriptome analysis reveals heterogeneous cell types in complex tissues and leads to unexpected biological findings when compared to bulk populations. However most of the methods focus on the 3'-end of polyadenylated transcripts using droplet-based technology. To achieve complete transcriptome, we describe single-cell 5'-end transcriptome protocol with random primed-cDNA harvesting on the Fluidigm C1™ platform which can isolate and process up to 96 cells from a single run with custom library preparation. The method enables detection of Transcription Start Site (TSS) at the single-cell resolution yielding a more comprehensive overview of gene regulatory elements governing in the EpiSC-like cell (EpiLC) including non-polyadenylated RNA and enhancer RNA activities.

Dynamic regulation of N6,2'-O-dimethyladenosine (m6Am) in obesity

The prevalent m⁶Am mRNA cap modification was recently identified as a valid target for removal by the human obesity gene FTO along with the previously established m⁶A mRNA modification. However, the deposition and dynamics of m⁶Am in regulating obesity are unknown. Here, we investigate the liver m⁶A/m methylomes in mice fed on a high fat Western-diet and in ob/ob mice. We find that FTO levels are elevated in fat mice, and that genes which lost m⁶Am marking under obesity are overly downregulated, including the two fatty-acid-binding proteins FABP2, and FABP5. Furthermore, the cellular perturbation of FTO correspondingly affect protein levels of its targets. Notably, generally m⁶Am- but not m⁶A-methylated genes, are found to be highly enriched in metabolic processes. Finally, we deplete all m⁶A background via Mettl3 knockout, and unequivocally uncover the association of m⁶Am methylation with increased mRNA stability, translation efficiency, and higher protein expression. Together, these results strongly implicate a dynamic role for m⁶Am in obesity-related translation regulation.

Global approaches for profiling transcription initiation

Transcription start site (TSS) selection influences transcript stability and translation as well as protein sequence. Alternative TSS usage is pervasive in organismal development, is a major contributor to transcript isoform diversity in humans, and is frequently observed in human diseases including cancer. In this review, we discuss the breadth of techniques that have been used to globally profile TSSs and the resulting insights into gene regulation, as well as future prospects in this area of inquiry.

Cap analysis of gene expression (CAGE) and noncoding regulatory elements

Cap analysis of gene expression (CAGE) was developed to detect the 5' end of RNA. Trapping of the RNA 5'-cap structure enables the enrichment and selective sequencing of complete transcripts. Upscaled high-throughput versions of CAGE have enabled the genome-wide identification of transcription start sites, including transcriptionally active promoters and enhancers. CAGE sequencing can be exploited to draw comprehensive maps of active genomic regulatory elements in a cell type- and activation-specific manner. The cells of the immune system are among the best candidates to be analyzed in humans, since they are easily accessible. In this review, we discuss how CAGE data are instrumental for integrative analyses with quantitative trait loci and omics data, and their usefulness in the mechanistic interpretation of the effects of genetic variations over the entire human genome. Integrating CAGE data with the currently available omics information will contribute to better understanding of the genome-wide association study variants that lie outside of annotated genes, deepening our knowledge on human diseases, and enabling the targeted design of more specific therapeutic interventions.

A Comprehensive Toolbox to Analyze Enhancer-Promoter Functions

Knowledge in gene transcription and chromatin regulation has been intensely studied for decades, but thanks to next-generation sequencing (NGS) techniques there has been a major leap forward in the last few years. Historically, identification of specific enhancer elements has led to the identification of master transcription factors (TFs) in the 1990s. Genetic and biochemical experiments have identified the key regulators controlling RNA polymerase II (RNAPII) transcription and structurally analyses have elucidated detailed mechanisms. NGS and the development of chromatin immunoprecipitation (ChIP) have accelerated the gain of knowledge in the recent years. By now, we have a dazzling wealth of techniques that are currently used to put gene expression into a genome-wide context. This book is an attempt to assemble useful protocols for many researchers within and nearby research areas. In general, these innovative techniques focus on enhancer and promoter studies. The techniques should also be of interest for related fields such as DNA repair and replication.

De novo activated transcription of inserted foreign coding sequences is inheritable in the plant genome

The manner in which inserted foreign coding sequences become transcriptionally activated and fixed in the plant genome is poorly understood. To examine such processes of gene evolution, we performed an artificial evolutionary experiment in Arabidopsis thaliana. As a model of gene-birth events, we introduced a promoterless coding sequence of the firefly luciferase (LUC) gene and established 386 T2-generation transgenic lines. Among them, we determined the individual LUC insertion loci in 76 lines and found that one-third of them were transcribed de novo even in the intergenic or inherently unexpressed regions. In the transcribed lines, transcription-related chromatin marks were detected across the newly activated transcribed regions. These results agreed with our previous findings in A. thaliana cultured cells under a similar experimental scheme. A comparison of the results of the T2-plant and cultured cell experiments revealed that the de novo-activated transcription concomitant with local chromatin remodelling was inheritable. During one-generation inheritance, it seems likely that the transcription activities of the LUC inserts trapped by the endogenous genes/transcripts became stronger, while those of de novo transcription in the intergenic/untranscribed regions became weaker. These findings may offer a clue for the elucidation of the mechanism by which inserted foreign coding sequences become transcriptionally activated and fixed in the plant genome.

MIR-NATs repress MAPT translation and aid proteostasis in neurodegeneration

The human genome expresses thousands of natural antisense transcripts (NAT) that can regulate epigenetic state, transcription, RNA stability or translation of their overlapping genes1,2. Here we describe MAPT-AS1, a brain-enriched NAT that is conserved in primates and contains an embedded mammalian-wide interspersed repeat (MIR), which represses tau translation by competing for ribosomal RNA pairing with the MAPT mRNA internal ribosome entry site3. MAPT encodes tau, a neuronal intrinsically disordered protein (IDP) that stabilizes axonal microtubules. Hyperphosphorylated, aggregation-prone tau forms the hallmark inclusions of tauopathies4. Mutations in MAPT cause familial frontotemporal dementia, and common variations forming the MAPT H1 haplotype are a significant risk factor in many tauopathies5 and Parkinson's disease. Notably, expression of MAPT-AS1 or minimal essential sequences from MAPT-AS1 (including MIR) reduces-whereas silencing MAPT-AS1 expression increases-neuronal tau levels, and correlate with tau pathology in human brain. Moreover, we identified many additional NATs with embedded MIRs (MIR-NATs), which are overrepresented at coding genes linked to neurodegeneration and/or encoding IDPs, and confirmed MIR-NAT-mediated translational control of one such gene, PLCG1. These results demonstrate a key role for MAPT-AS1 in tauopathies and reveal a potentially broad contribution of MIR-NATs to the tightly controlled translation of IDPs6, with particular relevance for proteostasis in neurodegeneration.

Approaches to Identify and Characterise the Post-Transcriptional Roles of lncRNAs in Cancer

It is becoming increasingly evident that the non-coding genome and transcriptome exert great influence over their coding counterparts through complex molecular interactions. Among non-coding RNAs (ncRNA), long non-coding RNAs (lncRNAs) in particular present increased potential to participate in dysregulation of post-transcriptional processes through both RNA and protein interactions. Since such processes can play key roles in contributing to cancer progression, it is desirable to continue expanding the search for lncRNAs impacting cancer through post-transcriptional mechanisms. The sheer diversity of mechanisms requires diverse resources and methods that have been developed and refined over the past decade. We provide an overview of computational resources as well as proven low-to-high throughput techniques to enable identification and characterisation of lncRNAs in their complex interactive contexts. As more cancer research strategies evolve to explore the non-coding genome and transcriptome, we anticipate this will provide a valuable primer and perspective of how these technologies have matured and will continue to evolve to assist researchers in elucidating post-transcriptional roles of lncRNAs in cancer.

Multi-omics analysis of hiPSCs-derived HLCs matured on-chip revealed patterns typical of liver regeneration

Maturation of human-induced pluripotent stem cells (hiPSCs)-derived hepatocytes-like cells (HLCs) toward a complete hepatocyte phenotype remains a challenge as primitiveness patterns are still commonly observed. In this study, we propose a modified differentiation protocol for those cells which includes a prematuration in Petri dishes and a maturation in microfluidic biochip. For the first time, a large range of biomolecular families has been extracted from the same sample to combine transcriptomic, proteomic, and metabolomic analysis. After integration, these datasets revealed specific molecular patterns and highlighted the hepatic regeneration profile in biochips. Overall, biochips exhibited processes of cell proliferation and inflammation (via TGFB1) coupled with anti-fibrotic signaling (via angiotensin 1-7, ATR-2, and MASR). Moreover, cultures in this condition displayed physiological lipid-carbohydrate homeostasis (notably via PPAR, cholesterol metabolism, and bile synthesis) coupled with cell respiration through advanced oxidative phosphorylation (through the overexpression of proteins from the third and fourth complex). The results presented provide an original overview of the complex mechanisms involved in liver regeneration using an advanced in vitro organ-on-chip technology.

Low Quantity Single Strand CAGE (LQ-ssCAGE) Maps Regulatory Enhancers and Promoters

The Cap Analysis of Gene Expression (CAGE) is a powerful method to identify Transcription Start Sites (TSSs) of capped RNAs while simultaneously measuring transcripts expression level. CAGE allows mapping at single nucleotide resolution at all active promoters and enhancers. Large CAGE datasets have been produced over the years from individual laboratories and consortia, including the Encyclopedia of DNA Elements (ENCODE) and Functional Annotation of the Mammalian Genome (FANTOM) consortia. These datasets constitute open resource for TSS annotations and gene expression analysis. Here, we provide an experimental protocol for the most recent CAGE method called Low Quantity (LQ) single strand (ss) CAGE "LQ-ssCAGE", which enables cost-effective profiling of low quantity RNA samples. LQ-ssCAGE is especially useful for samples derived from cells cultured in small volumes, cellular compartments such as nuclear RNAs or for samples from developmental stages. We demonstrate the reproducibility and effectiveness of the method by constructing 240 LQ-ssCAGE libraries from 50 ng of THP-1 cell extracted RNAs and discover lowly expressed novel enhancer and promoter-derived lncRNAs.

Use of Cap Analysis Gene Expression to detect human papillomavirus promoter activity patterns at different disease stages

Transcription of human papillomavirus (HPV) genes proceeds unidirectionally from multiple promoters. Direct profiling of transcription start sites (TSSs) by Cap Analysis Gene Expression (CAGE) is a powerful strategy for examining individual HPV promoter activity. The objective of this study was to evaluate alterations of viral promoter activity during infection using CAGE technology. We used CAGE-based sequencing of 46 primary cervical samples, and quantitatively evaluated TSS patterns in the HPV transcriptome at a single-nucleotide resolution. TSS patterns were classified into two types: early promoter-dominant type (Type A) and late promoter-dominant type (Type B). The Type B pattern was more frequently found in CIN1 and CIN2 lesions than in CIN3 and cancer samples. We detected transcriptomes from multiple HPV types in five samples. Interestingly, in each sample, the TSS patterns of both HPV types were the same. The viral gene expression pattern was determined by the differentiation status of the epithelial cells, regardless of HPV type. We performed unbiased analyses of TSSs across the HPV genome in clinical samples. Visualising TSS pattern dynamics, including TSS shifts, provides new insights into how HPV infection status relates to disease state.

Simple and efficient profiling of transcription initiation and transcript levels with STRIPE-seq

Accurate mapping of transcription start sites (TSSs) is key for understanding transcriptional regulation. However, current protocols for genome-wide TSS profiling are laborious and/or expensive. We present Survey of TRanscription Initiation at Promoter Elements with high-throughput sequencing (STRIPE-seq), a simple, rapid, and cost-effective protocol for sequencing capped RNA 5' ends from as little as 50 ng total RNA. Including depletion of uncapped RNA and reaction cleanups, a STRIPE-seq library can be constructed in about 5 h. We show application of STRIPE-seq to TSS profiling in yeast and human cells and show that it can also be effectively used for quantification of transcript levels and analysis of differential gene expression. In conjunction with our ready-to-use computational workflows, STRIPE-seq is a straightforward, efficient means by which to probe the landscape of transcriptional initiation.

Characterization of liver zonation-like transcriptomic patterns in HLCs derived from hiPSCs in a microfluidic biochip environment

The liver zonation is an important phenomenon characterized by a gradient of several functions along the liver acinus. However, this gradient remains difficult to reproduce in in-vitro conditions, making the obtention of an in-vitro method to recapitulate the liver zonation a challenging issue. In this study, we evaluated the spatial evolution of the transcriptome profile of human induced pluripotent stem cells (hiPSCs) differentiated toward hepatocytes-like cells (HLCs) phenotype in a microfluidic biochip environment. Cells collected at the inlet of the biochip, where the oxygen concentration is higher, were identified by the expression of genes involved in metabolic pathways related to cellular reorganization and cell proliferation. Cells collected in the middle and at the outlet of the biochips, where oxygen concentrations are lower, were characterized by the upregulation of genes involved in cellular detoxification processes (CYP450), PPAR signaling or arginine biosynthesis. A subset of 16 transcription factors (TFs) was extracted and identified as upstream regulators to HNF1A and PPARA. These TFs are also known as regulators to target genes engaged in the Wnt/βcatenin pathway, in the TGFβ/BMP/SMAD signaling, in the transition between epithelial mesenchymal transition (EMT) and mesenchymal epithelial transition (MET), in the homeostasis of lipids, bile acids and carbohydrates homeostasis, in drug metabolism, in the estrogen processing and in the oxidative stress response. Overall, the analysis allowed to confirm a partial zonation-like pattern in hiPSCs-derived HLCs in the microfluidic biochip environment. These results provide important insights into the reproduction of liver zonation in-vitro for a better understanding of the phenomenon.

Analysis of hiPSCs differentiation toward hepatocyte-like cells upon extended exposition to oncostatin

The capability to produce and maintain functional human adult hepatocytes remains one of the major challenges for the use of in-vitro models toward liver cell therapy and industrial drug-screening applications. Among the suggested strategies to solve this issue, the use of human-induced pluripotent stem cells (hiPSCs), differentiated toward hepatocyte-like cells (HLCs) is promising. In this work, we propose a 31-day long protocol, that includes a final 14-day long phase of oncostatin treatment, as opposed to a 7-day treatment which led to the formation of a hepatic tissue functional for CYP1A2, CYP2B6, CYP2C8, CYP2D6, and CYP3A4. The production of albumin, as well as bile acid metabolism and transport, were also detected. Transcriptome profile comparisons and liver transcription factors (TFs) motif dynamics revealed increased expression of typical hepatic markers such as HNF1A and of important metabolic markers like PPARA. The performed analysis has allowed for the extraction of potential targets and pathways which would allow enhanced hepatic maturation in-vitro. From this investigation, NRF1 and SP3 appeared as transcription factors of importance. Complex epithelial-mesenchymal transition (EMT) and mesenchymal-epithelial transition (MET) patterns were also observed during the differentiation process. Moreover, whole transcriptome analysis highlighted a response typical of the one observed in liver regeneration and hepatocyte proliferation. While a complete maturation of hepatocytes was yet to be obtained, the results presented in this work provide new insights into the process of liver development and highlight potential targets aimed to improve in-vitro liver regeneration.

Machine-driven parameter screen of biochemical reactions

The development of complex methods in molecular biology is a laborious, costly, iterative and often intuition-bound process where optima are sought in a multidimensional parameter space through step-by-step optimizations. The difficulty of miniaturizing reactions under the microliter volumes usually handled in multiwell plates by robots, plus the cost of the experiments, limit the number of parameters and the dynamic ranges that can be explored. Nevertheless, because of non-linearities of the response of biochemical systems to their reagent concentrations, broad dynamic ranges are necessary. Here we use a high-performance nanoliter handling platform and computer generation of liquid transfer programs to explore in quadruplicates 648 combinations of 4 parameters of a biochemical reaction, the reverse-transcription, which lead us to uncover non-linear responses, parameter interactions and novel mechanistic insights. With the increased availability of computer-driven laboratory platforms for biotechnology, our results demonstrate the feasibility and advantage of methods development based on reproducible, computer-aided exhaustive characterization of biochemical systems.

A human minisatellite hosts an alternative transcription start site for NPRL3 driving its expression in a repeat number-dependent manner

Minisatellites, also called variable number of tandem repeats (VNTRs), are a class of repetitive elements that may affect gene expression at multiple levels and have been correlated to disease. Their identification and role as expression quantitative trait loci (eQTL) have been limited by their absence in comparative genomic hybridization and single nucleotide polymorphisms arrays. By taking advantage of cap analysis of gene expression (CAGE), we describe a new example of a minisatellite hosting a transcription start site (TSS) which expression is dependent on the repeat number. It is located in the third intron of the gene nitrogen permease regulator like protein 3 (NPRL3). NPRL3 is a component of the GAP activity toward rags 1 protein complex that inhibits mammalian target of rapamycin complex 1 (mTORC1) activity and it is found mutated in familial focal cortical dysplasia and familial focal epilepsy. CAGE tags represent an alternative TSS identifying TAGNPRL3 messenger RNAs (mRNAs). TAGNPRL3 is expressed in red blood cells both at mRNA and protein levels, it interacts with its protein partner NPRL2 and its overexpression inhibits cell proliferation. This study provides an example of a minisatellite that is both a TSS and an eQTL as well as identifies a new VNTR that may modify mTORC1 activity.

Transcriptome profiling of hiPSC-derived LSECs with nanoCAGE

Liver Sinusoidal Endothelial Cells (LSECs) are an important component of the liver as they compose the microvasculature which allows the supply of oxygen, blood, and nutrients. However, maintenance of these cells in vitro remains challenging as they tend to rapidly lose some of their characteristics such as fenestration or as their immortalized counterparts present poor characteristics. In this work, human induced pluripotent stem cells (hiPSCs) have been differentiated toward an LSEC phenotype. After differentiation, the RNA quantification allowed demonstration of high expression of specific vascular markers (CD31, CD144, and STAB2). Immunostaining performed on the cells was found to be positive for both Stabilin-1 and Stabilin-2. Whole transcriptome analysis performed with the nanoCAGE method further confirmed the overall vascular commitment of the cells. The gene expression profile revealed the upregulation of the APLN, LYVE1, VWF, ESAM and ANGPT2 genes while VEGFA appeared to be downregulated. Analysis of promoter motif activities highlighted several transcription factors (TFs) of interest in LSECs (IRF2, ERG, MEIS2, SPI1, IRF7, WRNIP1, HIC2, NFIX_NFIB, BATF, and PATZ1). Based on this investigation, we compiled the regulatory network involving the relevant TFs, their target genes as well as their related signaling pathways. The proposed hiPSC-derived LSEC model and its regulatory network were then confirmed by comparing the experimental data to primary human LSEC reference datasets. Thus, the presented model appears as a promising tool to generate more complex in vitro liver multi-cellular tissues.

Cap Analysis of Gene Expression (CAGE): A Quantitative and Genome-Wide Assay of Transcription Start Sites

Cap analysis of gene expression (CAGE) is an approach to identify and monitor the activity (transcription initiation frequency) of transcription start sites (TSSs) at single base-pair resolution across the genome. It has been effectively used to identify active promoter and enhancer regions in cancer cells, with potential utility to identify key factors to immunotherapy. Here, we overview a series of CAGE protocols and describe detailed experimental steps of the latest protocol based on the Illumina sequencing platform; both experimental steps (see Subheadings 3.1-3.11) and computational processing steps (see Subheadings 3.12-3.20) are described.

MAPCap allows high-resolution detection and differential expression analysis of transcription start sites

The position, shape and number of transcription start sites (TSS) are critical determinants of gene regulation. Most methods developed to detect TSSs and study promoter usage are, however, of limited use in studies that demand quantification of expression changes between two or more groups. In this study, we combine high-resolution detection of transcription start sites and differential expression analysis using a simplified TSS quantification protocol, MAPCap (Multiplexed Affinity Purification of Capped RNA) along with the software icetea. Applying MAPCap on developing Drosophila melanogaster embryos and larvae, we detected stage and sex-specific promoter and enhancer activity and quantify the effect of mutants of maleless (MLE) helicase at X-chromosomal promoters. We observe that MLE mutation leads to a median 1.9 fold drop in expression of X-chromosome promoters and affects the expression of several TSSs with a sexually dimorphic expression on autosomes. Our results provide quantitative insights into promoter activity during dosage compensation.

The position, shape and number of transcription start sites (TSS) regulate gene expression. Here authors present MAPCap, a method for high-resolution detection and differential expression analysis of TSS, and apply MAPCap to early fly development, detecting stage and sex-specific promoter and enhancer activity.

Using RAMPAGE to Identify and Annotate Promoters in Insect Genomes

Application of Transcription Start Site (TSS) profiling technologies, coupled with large-scale next-generation sequencing (NGS) has yielded valuable insights into the location, structure, and activity of promoters across diverse metazoan model systems. In insects, TSS profiling has been used to characterize the promoter architecture of Drosophila melanogaster (Hoskins et al., Genome Res 21(2):182-192, 2011) and subsequently was employed to reveal widespread transposon-driven alternative promoter usage in the fruit fly (Batut et al., Genome Res 23:169-180, 2012).In this chapter we discuss the computational analysis of the experimental data derived from one of TSS profiling methods, RAMPAGE (RNA Annotation and Mapping of Promoters for Analysis of Gene Expression) that can be used for the precise, quantitative identification of promoters in insect genomes. We demonstrate this using the software tools GoRAMPAGE (Brendel and Raborn, GoRAMPAGE-A workflow for promoter detection by 5'-read mapping. https://github.com/BrendelGroup/GoRAMPAGE , 2016) and TSRchitect (Raborn and Brendel, TSRchitect: promoter identification from large-scale TSS profiling data. R Bioconductor package version 1.8.0 [Online]. Available: http://bioconductor.org/packages/release/bioc/html/TSRchitect.html , 2017), providing detailed instructions with the aim of taking the user from raw reads to processed results.

Alternative polyadenylation coordinates embryonic development, sexual dimorphism and longitudinal growth in Xenopus tropicalis

RNA alternative polyadenylation contributes to the complexity of information transfer from genome to phenome, thus amplifying gene function. Here, we report the first X. tropicalis resource with 127,914 alternative polyadenylation (APA) sites derived from embryos and adults. Overall, APA networks play central roles in coordinating the maternal-zygotic transition (MZT) in embryos, sexual dimorphism in adults and longitudinal growth from embryos to adults. APA sites coordinate reprogramming in embryos before the MZT, but developmental events after the MZT due to zygotic genome activation. The APA transcriptomes of young adults are more variable than growing adults and male frog APA transcriptomes are more divergent than females. The APA profiles of young females were similar to embryos before the MZT. Enriched pathways in developing embryos were distinct across the MZT and noticeably segregated from adults. Briefly, our results suggest that the minimal functional units in genomes are alternative transcripts as opposed to genes.

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 C1^TM 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.

Single-cell transcriptomic profiling allows the exploration of cellular heterogeneity but commonly focuses on the 3′-end of the transcript. Here the authors introduce C1 CAGE, which detects the 5′ transcript end in a multiplexed microfluidic system.

Analysis of the transcription factors and their regulatory roles during a step-by-step differentiation of induced pluripotent stem cells into hepatocyte-like cells

Human induced pluripotent stem cells have been investigated through a sequential in vitro step-by-step differentiation into hepatocyte-like cells using nanoCAGE, an original method for promoters, transcription factors, and transcriptome analysis.

SLIC-CAGE: high-resolution transcription start site mapping using nanogram-levels of total RNA

Cap analysis of gene expression (CAGE) is a methodology for genome-wide quantitative mapping of mRNA 5′ ends to precisely capture transcription start sites at a single nucleotide resolution. In combination with high-throughput sequencing, CAGE has revolutionized our understanding of the rules of transcription initiation, led to discovery of new core promoter sequence features, and discovered transcription initiation at enhancers genome-wide. The biggest limitation of CAGE is that even the most recently improved version (nAnT-iCAGE) still requires large amounts of total cellular RNA (5 µg), preventing its application to scarce biological samples such as those from early embryonic development or rare cell types. Here, we present SLIC-CAGE, a Super-Low Input Carrier-CAGE approach to capture 5′ ends of RNA polymerase II transcripts from as little as 5–10 ng of total RNA. This dramatic increase in sensitivity is achieved by specially designed, selectively degradable carrier RNA. We demonstrate the ability of SLIC-CAGE to generate data for genome-wide promoterome with 1000-fold less material than required by existing CAGE methods, by generating a complex, high-quality library from mouse embryonic day 11.5 primordial germ cells.

m1A within cytoplasmic mRNAs at single nucleotide resolution: a reconciled transcriptome-wide map

Following synthesis, RNA can be modified with over 100 chemically distinct modifications. Recently, two studies-one by our group-developed conceptually similar approaches to globally map N1-methyladenosine (m¹A) at single nucleotide resolution. Surprisingly, the studies diverged quite substantially in their estimates of the abundance, whereabouts, and stoichiometry of m¹A within internal sites in cytosolic mRNAs: One study reported it to be a very rare modification, present at very low stoichiometries, and invariably catalyzed by TRMT6/61A. The other found it to be present at >470 sites, often at high levels, and suggested that the vast majority were highly unlikely to be TRMT6/61A substrates. Here we reanalyze the data from the latter study, and demonstrate that the vast majority of the detected sites originate from duplications, misannotations, mismapping, SNPs, sequencing errors, and a set of sites from the very first transcribed base that appear to originate from nontemplated incorporations by reverse transcriptase. Only 53 of the sites detected in the latter study likely reflect bona-fide internal modifications of cytoplasmically encoded mRNA molecules, nearly all of which are likely TRMT6/TRMT61A substrates and typically modified at low to undetectable levels. The experimental data sets from both studies thus consistently demonstrate that within cytosolic mRNAs, m¹A is a rare internal modification where it is typically catalyzed at very low stoichiometries via a single complex. Our findings offer a clear and consistent view on the abundance and whereabouts of m¹A, and lay out directions for future studies.

m1A within cytoplasmic mRNAs at single nucleotide resolution: a reconciled transcriptome-wide map

Following synthesis, RNA can be modified with over 100 chemically distinct modifications. Recently, two studies-one by our group-developed conceptually similar approaches to globally map N1-methyladenosine (m¹A) at single nucleotide resolution. Surprisingly, the studies diverged quite substantially in their estimates of the abundance, whereabouts, and stoichiometry of m¹A within internal sites in cytosolic mRNAs: One study reported it to be a very rare modification, present at very low stoichiometries, and invariably catalyzed by TRMT6/61A. The other found it to be present at >470 sites, often at high levels, and suggested that the vast majority were highly unlikely to be TRMT6/61A substrates. Here we reanalyze the data from the latter study, and demonstrate that the vast majority of the detected sites originate from duplications, misannotations, mismapping, SNPs, sequencing errors, and a set of sites from the very first transcribed base that appear to originate from nontemplated incorporations by reverse transcriptase. Only 53 of the sites detected in the latter study likely reflect bona-fide internal modifications of cytoplasmically encoded mRNA molecules, nearly all of which are likely TRMT6/TRMT61A substrates and typically modified at low to undetectable levels. The experimental data sets from both studies thus consistently demonstrate that within cytosolic mRNAs, m¹A is a rare internal modification where it is typically catalyzed at very low stoichiometries via a single complex. Our findings offer a clear and consistent view on the abundance and whereabouts of m¹A, and lay out directions for future studies.

Deep Cap Analysis of Gene Expression (CAGE): Genome-Wide Identification of Promoters, Quantification of Their Activity, and Transcriptional Network Inference

Among the most significant findings of the post-genomic era, the discovery of pervasive transcription of mammalian genomes has tremendously modified our understanding of the genome output seen as RNA molecules. The increased focus on non-protein-coding genomic regions together with concomitant technological innovations has led to rapid discovery of numerous noncoding transcripts (ncRNAs). Biological relevance and functional roles of the vast majority of these ncRNAs remain largely unknown.The cap analysis of gene expression (CAGE) technology allows accurate transcript detection and quantification without relying on preexisting transcript models. In combination with complementary data sets, generated using other technologies, it has been shown as an efficient approach for exploring transcriptome complexity.Here, we describe the use of CAGE for the identification of novel noncoding transcripts in mammalian cells providing detailed information for basic data processing and advanced bioinformatics analyses.

NanoCAGE: A Method for the Analysis of Coding and Noncoding 5'-Capped Transcriptomes

Transcripts in all eukaryotes are characterized by the 5'-end specific cap structure in mRNAs. Cap Analysis Gene Expression or CAGE makes use of these caps to specifically obtain cDNA fragments from the 5'-end of RNA and sequences those at high throughput for transcript identification and genome-wide mapping of transcription start sites for coding and noncoding genes. Here, we provide an improved version of our nanoCAGE protocol that has been developed for preparing CAGE libraries from as little as 50 ng of total RNA within three standard working days. Key steps in library preparation have been improved over our previously published protocol to obtain libraries having a good 5'-end selection and a more equal size distribution for higher sequencing efficiency on Illumina MiSeq and HiSeq sequencers. We recommend nanoCAGE as the method of choice for transcriptome profiling projects even from limited amounts of RNA, and as the best approach for genome-wide mapping of transcription start sites within promoter regions.

Technological Developments in lncRNA Biology

It is estimated that more than 90% of the mammalian genome is transcribed as non-coding RNAs. Recent evidences have established that these non-coding transcripts are not junk or just transcriptional noise, but they do serve important biological purpose. One of the rapidly expanding fields of this class of transcripts is the regulatory lncRNAs, which had been a major challenge in terms of their molecular functions and mechanisms of action. The emergence of high-throughput technologies and the development in various conventional approaches have led to the expansion of the lncRNA world. The combination of multidisciplinary approaches has proven to be essential to unravel the complexity of their regulatory networks and helped establish the importance of their existence. Here, we review the current methodologies available for discovering and investigating functions of long non-coding RNAs (lncRNAs) and focus on the powerful technological advancement available to specifically address their functional importance.

Defining the 5΄ and 3΄ landscape of the Drosophila transcriptome with Exo-seq and RNaseH-seq

Cells regulate biological responses in part through changes in transcription start sites (TSS) or cleavage and polyadenylation sites (PAS). To fully understand gene regulatory networks, it is therefore critical to accurately annotate cell type-specific TSS and PAS. Here we present a simple and straightforward approach for genome-wide annotation of 5΄- and 3΄-RNA ends. Our approach reliably discerns bona fide PAS from false PAS that arise due to internal poly(A) tracts, a common problem with current PAS annotation methods. We applied our methodology to study the impact of temperature on the Drosophila melanogaster head transcriptome. We found hundreds of previously unidentified TSS and PAS which revealed two interesting phenomena: first, genes with multiple PASs tend to harbor a motif near the most proximal PAS, which likely represents a new cleavage and polyadenylation signal. Second, motif analysis of promoters of genes affected by temperature suggested that boundary element association factor of 32 kDa (BEAF-32) and DREF mediates a transcriptional program at warm temperatures, a result we validated in a fly line where beaf-32 is downregulated. These results demonstrate the utility of a high-throughput platform for complete experimental and computational analysis of mRNA-ends to improve gene annotation.

Untranslated Parts of Genes Interpreted: Making Heads or Tails of High-Throughput Transcriptomic Data via Computational Methods: Computational methods to discover and quantify isoforms with alternative untranslated regions

In this review we highlight the importance of defining the untranslated parts of transcripts, and present a number of computational approaches for the discovery and quantification of alternative transcription start and poly-adenylation events in high-throughput transcriptomic data. The fate of eukaryotic transcripts is closely linked to their untranslated regions, which are determined by the position at which transcription starts and ends at a genomic locus. Although the extent of alternative transcription starts and alternative poly-adenylation sites has been revealed by sequencing methods focused on the ends of transcripts, the application of these methods is not yet widely adopted by the community. We suggest that computational methods applied to standard high-throughput technologies are a useful, albeit less accurate, alternative to the expertise-demanding 5' and 3' sequencing and they are the only option for analysing legacy transcriptomic data. We review these methods here, focusing on technical challenges and arguing for the need to include better normalization of the data and more appropriate statistical models of the expected variation in the signal.

Linking FANTOM5 CAGE peaks to annotations with CAGEscan

The FANTOM5 expression atlas is a quantitative measurement of the activity of nearly 200,000 promoter regions across nearly 2,000 different human primary cells, tissue types and cell lines. Generation of this atlas was made possible by the use of CAGE, an experimental approach to localise transcription start sites at single-nucleotide resolution by sequencing the 5′ ends of capped RNAs after their conversion to cDNAs. While 50% of CAGE-defined promoter regions could be confidently associated to adjacent transcriptional units, nearly 100,000 promoter regions remained gene-orphan. To address this, we used the CAGEscan method, in which random-primed 5′-cDNAs are paired-end sequenced. Pairs starting in the same region are assembled in transcript models called CAGEscan clusters. Here, we present the production and quality control of CAGEscan libraries from 56 FANTOM5 RNA sources, which enhances the FANTOM5 expression atlas by providing experimental evidence associating core promoter regions with their cognate transcripts.

Context-specific effects of genetic variants associated with autoimmune disease

Autoimmune diseases such as rheumatoid arthritis and coeliac disease are typical examples of complex genetic diseases caused by a combination of genetic and non-genetic risk factors. Insight into the genetic risk factors (single nucleotide polymorphisms (SNPs)) has increased since genome-wide association studies (GWAS) became possible in 2007 and, for individual diseases, SNPs can now explain some 15-50% of genetic risk. GWAS have also shown that some 50% of the genetic risk factors for individual autoimmune diseases overlap between different diseases. Thus, shared risk factors may converge to pathways that, when perturbed by genetic variation, predispose to autoimmunity in general. This raises the question of what determines disease specificity, and suggests that identical risk factors may have different effects in various autoimmune diseases. Addressing this question requires translation of genetic risk factors to causal genes and then to molecular and cellular pathways. Since >90% of the genetic risk factors are found in the non-coding part of the genome (i.e. outside the exons of protein-coding genes) and can have an impact on gene regulation, there is an urgent need to better understand the non-coding part of the genome. Here, we will outline the methods being used to unravel the gene regulatory networks perturbed in autoimmune diseases and the importance of doing this in the relevant cell types. We will highlight findings in coeliac disease, which manifests in the small intestine, to demonstrate how cell type and disease context can impact on the consequences of genetic risk factors.

Acute doses of caffeine shift nervous system cell expression profiles toward promotion of neuronal projection growth

Caffeine is a widely consumed psychoactive substance, but little is known about the effects of caffeine stimulation on global gene expression changes in neurons. Here, we conducted gene expression profiling of human neuroepithelial stem cell-derived neurons, stimulated with normal consumption levels of caffeine (3 μM and 10 μM), over a period of 9 h. We found dosage-dependent activation of immediate early genes after 1 h. Neuronal projection development processes were up-regulated and negative regulation of axon extension processes were down-regulated at 3 h. In addition, genes involved in extracellular matrix organization, response for wound healing, and regulation of immune system processes were down-regulated by caffeine at 3 h. This study identified novel genes within the neuronal projection guidance pathways that respond to acute caffeine stimulation and suggests potential mechanisms for the effects of caffeine on neuronal cells.