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

Cis-regulatory evolution of the recently expanded Ly49 gene family

Comparative genomics has revealed the rapid expansion of multiple gene families involved in immunity. Members within each gene family often evolved distinct roles in immunity. However, less is known about the evolution of their epigenome and cis-regulation. Here we systematically profile the epigenome of the recently expanded murine Ly49 gene family that mainly encode either inhibitory or activating surface receptors on natural killer cells. We identify a set of cis-regulatory elements (CREs) for activating Ly49 genes. In addition, we show that in mice, inhibitory and activating Ly49 genes are regulated by two separate sets of proximal CREs, likely resulting from lineage-specific losses of CRE activity. Furthermore, we find that some Ly49 genes are cross-regulated by the CREs of other Ly49 genes, suggesting that the Ly49 family has begun to evolve a concerted cis-regulatory mechanism. Collectively, we demonstrate the different modes of cis-regulatory evolution for a rapidly expanding gene family.

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.

A single-cell genomic strategy for alternative transcript start sites identification

Alternative transcription start sites (TSSs) usage plays a critical role in gene transcription regulation in mammals. However, precisely identifying alternative TSSs remains challenging at the genome-wide level. We report a single-cell genomic technology for alternative TSSs annotation and cell heterogeneity detection. In the method, we utilize Fluidigm C1 system to capture individual cells of interest, SMARTer cDNA synthesis kit to recover full-length cDNAs, then dual priming oligonucleotide system to specifically enrich TSSs for genomic analysis. We apply this method to a genome-wide study of alternative TSSs identification in two different IFN-β stimulated mouse embryonic fibroblasts (MEFs). The data clearly discriminate two IFN-β stimulated MEFs. Moreover, our results indicate 81% expressed genes in these two cell types containing multiple TSSs, which is much higher than previous predictions based on Cap-Analysis Gene Expression (CAGE) (58%) or empirical determination (54%) in various cell types. This indicates that alternative TSSs are more pervasive than expected and implies our strategy could position them at an unprecedented sensitivity. It would be helpful for elucidating their biological insights in future.

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.

Epigenomic analysis reveals prevalent contribution of transposable elements to cis-regulatory elements, tissue-specific expression, and alternative promoters in zebrafish

Transposable elements (TEs) encode regulatory elements that impact gene expression in multiple species, yet a comprehensive analysis of zebrafish TEs in the context of gene regulation is lacking. Here, we systematically investigate the epigenomic and transcriptomic landscape of TEs across 11 adult zebrafish tissues using multidimensional sequencing data. We find that TEs contribute substantially to a diverse array of regulatory elements in the zebrafish genome and that 37% of TEs are positioned in active regulatory states in adult zebrafish tissues. We identify TE subfamilies enriched in highly specific regulatory elements among different tissues. We use transcript assembly to discover TE-derived transcriptional units expressed across tissues. Finally, we show that novel TE-derived promoters can initiate tissue-specific transcription of alternate gene isoforms. This work provides a comprehensive profile of TE activity across normal zebrafish tissues, shedding light on mechanisms underlying the regulation of gene expression in this widely used model organism.

Genome-Wide Profiling of Transcription Initiation with STRIPE-seq

Transcription start site (TSS) usage is a critical factor in the regulation of gene expression. A number of methods for global TSS mapping have been developed, but barriers of expense, technical difficulty, time, and/or cost have limited their broader adoption. To address these issues, we developed Survey of TRanscription Initiation at Promoter Elements with high-throughput sequencing (STRIPE-seq). Requiring only three enzymatic steps with intervening bead cleanups, a STRIPE-seq library can be prepared from as little as 50 ng total RNA in ~5 h at a cost of ~$12 (US). In addition to profiling TSS usage, STRIPE-seq provides information on transcript levels that can be used for differential expression analysis. Thanks to its simplicity and low cost, we envision that STRIPE-seq could be employed by any molecular biology laboratory interested in profiling transcription initiation.

Long Non-coding RNAs: Mechanisms, Experimental, and Computational Approaches in Identification, Characterization, and Their Biomarker Potential in Cancer

Long non-coding RNAs are diverse class of non-coding RNA molecules >200 base pairs of length having various functions like gene regulation, dosage compensation, epigenetic regulation. Dysregulation and genomic variations of several lncRNAs have been implicated in several diseases. Their tissue and developmental specific expression are contributing factors for them to be viable indicators of physiological states of the cells. Here we present an comprehensive review the molecular mechanisms and functions, state of the art experimental and computational pipelines and challenges involved in the identification and functional annotation of lncRNAs and their prospects as biomarkers. We also illustrate the application of co-expression networks on the TCGA-LIHC dataset for putative functional predictions of lncRNAs having a therapeutic potential in Hepatocellular carcinoma (HCC).

Investigation of the hepatic development in the coculture of hiPSCs-derived LSECs and HLCs in a fluidic microenvironment

Interactions between the different liver cell types are critical to the maintenance or induction of their function in vitro. In this work, human-induced Pluripotent Stem Cells (hiPSCs)-derived Liver Sinusoidal Endothelial Cells (LSECs) and Hepatocytes-Like Cells (HLCs) were cultured and matured in a microfluidic environment. Both cell populations were differentiated in Petri dishes, detached, and inoculated in microfluidic biochips. In cocultures of both cell types, the tissue has exhibited a higher production of albumin (3.19 vs 5.31 μg/mL/106 cells in monocultures and cocultures) as well as a higher inducibility CYP450 over monocultures of HLCs. Tubular-like structures composed of LSECs and positive for the endothelial marker PECAM1, as well as a tissue more largely expressing Stabilin-2 were detected in cocultures only. In contrast, monocultures exhibited no network and less specific endothelial markers. The transcriptomic analysis did not reveal a marked difference between the profiles of both culture conditions. Nevertheless, the analysis allowed us to highlight different upstream regulators in cocultures (SP1, EBF1, and GATA3) and monocultures (PML, MECP2, and NRF1). In cocultures, the multi-omics dataset after 14 days of maturation in biochips has shown the activation of signaling related to hepatic maturation, angiogenesis, and tissue repair. In this condition, inflammatory signaling was also found to be reduced when compared to monocultures as illustrated by the activation of NFKB and by the detection of several cytokines involved in tissue injury in the latter. Finally, the extracted biological processes were discussed regarding the future development of a new generation of human in vitro hepatic models.

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.

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.

Embryonic tissue differentiation is characterized by transitions in cell cycle dynamic-associated core promoter regulation

The core-promoter, a stretch of DNA surrounding the transcription start site (TSS), is a major integration-point for regulatory-signals controlling gene-transcription. Cellular differentiation is marked by divergence in transcriptional repertoire and cell-cycling behaviour between cells of different fates. The role promoter-associated gene-regulatory-networks play in development-associated transitions in cell-cycle-dynamics is poorly understood. This study demonstrates in a vertebrate embryo, how core-promoter variations define transcriptional output in cells transitioning from a proliferative to cell-lineage specifying phenotype. Assessment of cell proliferation across zebrafish embryo segmentation, using the FUCCI transgenic cell-cycle-phase marker, revealed a spatial and lineage-specific separation in cell-cycling behaviour. To investigate the role differential promoter usage plays in this process, cap-analysis-of-gene-expression (CAGE) was performed on cells segregated by cycling dynamics. This analysis revealed a dramatic increase in tissue-specific gene expression, concurrent with slowed cycling behaviour. We revealed a distinct sharpening in TSS utilization in genes upregulated in slowly cycling, differentiating tissues, associated with enhanced utilization of the TATA-box, in addition to Sp1 binding-sites. In contrast, genes upregulated in rapidly cycling cells carry broad distribution of TSS utilization, coupled with enrichment for the CCAAT-box. These promoter features appear to correspond to cell-cycle-dynamic rather than tissue/cell-lineage origin. Moreover, we observed genes with cell-cycle-dynamic-associated transitioning in TSS distribution and differential utilization of alternative promoters. These results demonstrate the regulatory role of core-promoters in cell-cycle-dependent transcription regulation, during embryo-development.

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.

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.

Pre- and peri-implantation Zika virus infection impairs fetal development by targeting trophectoderm cells

Zika virus (ZIKV) infection results in an increased risk of spontaneous abortion and poor intrauterine growth although the underlying mechanisms remain undetermined. Little is known about the impact of ZIKV infection during the earliest stages of pregnancy, at pre- and peri-implantation, because most current ZIKV pregnancy studies have focused on post-implantation stages. Here, we demonstrate that trophectoderm cells of pre-implantation human and mouse embryos can be infected with ZIKV, and propagate virus causing neural progenitor cell death. These findings are corroborated by the dose-dependent nature of ZIKV susceptibility of hESC-derived trophectoderm cells. Single blastocyst RNA-seq reveals key transcriptional changes upon ZIKV infection, including nervous system development, prior to commitment to the neural lineage. The pregnancy rate of mice is >50% lower in pre-implantation infection than infection at E4.5, demonstrating that pre-implantation ZIKV infection leads to miscarriage. Cumulatively, these data elucidate a previously unappreciated association of pre- and peri-implantation ZIKV infection and microcephaly.

Broad Heterochromatic Domains Open in Gonocyte Development Prior to De Novo DNA Methylation

Facultative heterochromatin forms and reorganizes in response to external stimuli. However, how the initial establishment of such a chromatin state is regulated in cell-cycle-arrested cells remains unexplored. Mouse gonocytes are arrested male germ cells, at which stage the genome-wide DNA methylome forms. Here, we discovered transiently accessible heterochromatin domains of several megabases in size in gonocytes and named them differentially accessible domains (DADs). Open DADs formed in gene desert and gene cluster regions, primarily at transposons, with the reprogramming of histone marks, suggesting DADs as facultative heterochromatin. De novo DNA methylation took place with two waves in gonocytes: the first region specific and the second genome-wide. DADs were resistant to the first wave and their opening preceded the second wave. In addition, the higher-order chromosome architecture was reorganized with less defined chromosome compartments in gonocytes. These findings suggest that multiple layers of chromatin reprogramming facilitate de novo DNA methylation.

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.

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.

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.

Systematic analysis of transcription start sites in avian development

Cap Analysis of Gene Expression (CAGE) in combination with single-molecule sequencing technology allows precision mapping of transcription start sites (TSSs) and genome-wide capture of promoter activities in differentiated and steady state cell populations. Much less is known about whether TSS profiling can characterize diverse and non-steady state cell populations, such as the approximately 400 transitory and heterogeneous cell types that arise during ontogeny of vertebrate animals. To gain such insight, we used the chick model and performed CAGE-based TSS analysis on embryonic samples covering the full 3-week developmental period. In total, 31,863 robust TSS peaks (>1 tag per million [TPM]) were mapped to the latest chicken genome assembly, of which 34% to 46% were active in any given developmental stage. ZENBU, a web-based, open-source platform, was used for interactive data exploration. TSSs of genes critical for lineage differentiation could be precisely mapped and their activities tracked throughout development, suggesting that non-steady state and heterogeneous cell populations are amenable to CAGE-based transcriptional analysis. Our study also uncovered a large set of extremely stable housekeeping TSSs and many novel stage-specific ones. We furthermore demonstrated that TSS mapping could expedite motif-based promoter analysis for regulatory modules associated with stage-specific and housekeeping genes. Finally, using Brachyury as an example, we provide evidence that precise TSS mapping in combination with Clustered Regularly Interspaced Short Palindromic Repeat (CRISPR)-on technology enables us, for the first time, to efficiently target endogenous avian genes for transcriptional activation. Taken together, our results represent the first report of genome-wide TSS mapping in birds and the first systematic developmental TSS analysis in any amniote species (birds and mammals). By facilitating promoter-based molecular analysis and genetic manipulation, our work also underscores the value of avian models in unravelling the complex regulatory mechanism of cell lineage specification during amniote development.