Vertebrate Embryogenesis Research Articles

Feedback regulation of retinaldehyde reductase DHRS3, a critical determinant of retinoic acid homeostasis.

Retinoic acid is crucial for vertebrate embryogenesis, influencing anterior-posterior patterning and organogenesis through its interaction with nuclear hormone receptors comprising heterodimers of retinoic acid receptors (RARα, β, or γ) and retinoid X receptors (RXRα, β, or γ). Tissue retinoic acid levels are tightly regulated since both its excess and deficiency are deleterious. Dehydrogenase/reductase 3 (DHRS3) plays a critical role in this regulation by converting retinaldehyde to retinol, preventing excessive retinoic acid formation. Mutations in DHRS3 can result in embryonic lethality and congenital defects. This study shows that mouse Dhrs3 expression is responsive to vitamin A status and is directly regulated by the RAR/RXR complex through cis-regulatory elements. This highlights a negative feedback mechanism that ensures retinoic acid homeostasis.

The dynamics and functional impact of tRNA repertoires during early embryogenesis in zebrafish.

Embryogenesis entails dramatic shifts in mRNA translation and turnover that reprogram gene expression during cellular proliferation and differentiation. Codon identity modulates mRNA stability during early vertebrate embryogenesis, but how the composition of tRNA pools is matched to translational demand is unknown. By quantitative profiling of tRNA repertoires in zebrafish embryos during the maternal-to-zygotic transition, we show that zygotic tRNA repertoires are established after the onset of gastrulation, succeeding the major wave of zygotic mRNA transcription. Maternal and zygotic tRNA pools are distinct, but their reprogramming does not result in a better match to the codon content of the zygotic transcriptome. Instead, we find that an increase in global translation at gastrulation sensitizes decoding rates to tRNA supply, thus destabilizing maternal mRNAs enriched in slowly translated codons. Translational activation and zygotic tRNA expression temporally coincide with an increase of TORC1 activity at gastrulation, which phosphorylates and inactivates the RNA polymerase III repressor Maf1a/b. Our data indicate that a switch in global translation, rather than tRNA reprogramming, determines the onset of codon-dependent maternal mRNA decay during zebrafish embryogenesis.

A population of Vasa2 and Piwi1 expressing cells generates germ cells and neurons in a sea anemone

Germline segregation, essential for protecting germ cells against mutations, occurs during early embryogenesis in vertebrates, insects and nematodes. Highly regenerative animals (e.g., cnidarians), however, retain stem cells with both germinal and somatic potentials throughout adulthood, but their biology and evolution remain poorly understood. Among cnidarians (e.g., sea anemones, jellyfish), stem cells are only known in few hydrozoans (e.g., Hydra). Here, we identify and characterize a rare, multipotent population of stem and/or progenitor cells expressing the conserved germline and multipotency proteins Vasa2 and Piwi1 in the sea anemone Nematostella vectensis. Using piwi1 and vasa2 transgenic reporter lines, we reveal that the Vasa2+/Piwi1+ cell population generates not only gametes, but also a diversity of proliferative somatic cells, including neural progenitors, in juveniles and adults. Our work has uncovered a multipotent population of Vasa2+/Piwi1+ stem/progenitor cells that forms the cellular basis to understand body plasticity and regenerative capacities in sea anemones and corals.

Expansion of a neural crest gene signature following ectopic MYCN expression in sympathoadrenal lineage cells in vivo.

Neural crest cells (NCC) are multipotent migratory stem cells that originate from the neural tube during early vertebrate embryogenesis. NCCs give rise to a variety of cell types within the developing organism, including neurons and glia of the sympathetic nervous system. It has been suggested that failure in correct NCC differentiation leads to several diseases, including neuroblastoma (NB). During normal NCC development, MYCN is transiently expressed to promote NCC migration, and its downregulation precedes neuronal differentiation. Overexpression of MYCN has been linked to high-risk and aggressive NB progression. For this reason, understanding the effect overexpression of this oncogene has on the development of NCC-derived sympathoadrenal progenitors (SAP), which later give rise to sympathetic nerves, will help elucidate the developmental mechanisms that may prime the onset of NB. Here, we found that overexpressing human EGFP-MYCN within SAP lineage cells in zebrafish led to the transient formation of an abnormal SAP population, which displayed expanded and elevated expression of NCC markers while paradoxically also co-expressing SAP and neuronal differentiation markers. The aberrant NCC signature was corroborated with in vivo time-lapse confocal imaging in zebrafish larvae, which revealed transient expansion of sox10 reporter expression in MYCN overexpressing SAPs during the early stages of SAP development. In these aberrant MYCN overexpressing SAP cells, we also found evidence of dampened BMP signaling activity, indicating that BMP signaling disruption occurs following elevated MYCN expression. Furthermore, we discovered that pharmacological inhibition of BMP signaling was sufficient to create an aberrant NCC gene signature in SAP cells, phenocopying MYCN overexpression. Together, our results suggest that MYCN overexpression in SAPs disrupts their differentiation by eliciting abnormal NCC gene expression programs, and dampening BMP signaling response, having developmental implications for the priming of NB in vivo.

Actin Depolymerizing Factor Destrin Regulates Cilia Development and Function during Vertebrate Embryogenesis.

The actin cytoskeleton plays fundamental roles in ciliogenesis and the actin depolymerizing factor destrin regulates actin dynamics by treadmilling actin filaments and increasing globular actin pools. However, the specific developmental roles of destrin in ciliogenesis have not been fully elucidated. Here, we investigated the function of destrin in ciliogenesis using Xenopus laevis and human retinal pigmented epithelial (hRPE1) cells. We discovered the loss of destrin increased the number of multiciliated cells in the Xenopus epithelium and impeded cilia motility. Additionally, destrin depletion remarkably reduced the length of primary cilia in the Xenopus neural tube and hRPE1 cells by affecting actin dynamics. Immunofluorescence using markers of ciliary components indicated that destrin controls the directionality and polarity of basal bodies and axonemal elongation by modulating actin dynamics, independent of basal body docking. In conclusion, destrin plays a significant role during vertebrate ciliogenesis regulating both primary and multicilia development. Our data suggest new insights for understanding the roles of actin dynamics in cilia development.

A dual ribosomal system in the zebrafish soma and germline.

Protein synthesis during vertebrate embryogenesis is driven by ribosomes of two distinct origins: maternal ribosomes synthesized during oogenesis and stored in the egg, and somatic ribosomes, produced by the developing embryo after zygotic genome activation (ZGA). In zebrafish, these two ribosome types are expressed from different genomic loci and also differ in their ribosomal RNA (rRNA) sequence. To characterize this dual ribosome system further, we examined the expression patterns of maternal and somatic rRNAs during embryogenesis and in adult tissues. We found that maternal rRNAs are not only expressed during oogenesis but are continuously produced in the zebrafish germline. Proteomic analyses of maternal and somatic ribosomes unveiled differences in core ribosomal protein composition. Most nucleotide differences between maternal and somatic rRNAs are located in the flexible, structurally not resolved expansion segments. Our in vivo data demonstrated that both maternal and somatic ribosomes can be translationally active in the embryo. Using transgenically tagged maternal or somatic ribosome subunits, we experimentally confirm the presence of hybrid 80S ribosomes composed of 40S and 60S subunits from both origins and demonstrate the preferential in vivo association of maternal ribosomes with germline-specific transcripts. Our study identifies a distinct type of ribosomes in the zebrafish germline and thus presents a foundation for future explorations into possible regulatory mechanisms and functional roles of heterogeneous ribosomes.

Dissecting the regulatory logic of specification and differentiation during vertebrate embryogenesis.

The interplay between transcription factors and chromatin accessibility regulates cell type diversification during vertebrate embryogenesis. To systematically decipher the gene regulatory logic guiding this process, we generated a single-cell multi-omics atlas of RNA expression and chromatin accessibility during early zebrafish embryogenesis. We developed a deep learning model to predict chromatin accessibility based on DNA sequence and found that a small number of transcription factors underlie cell-type-specific chromatin landscapes. While Nanog is well-established in promoting pluripotency, we discovered a new function in priming the enhancer accessibility of mesendodermal genes. In addition to the classical stepwise mode of differentiation, we describe instant differentiation, where pluripotent cells skip intermediate fate transitions and terminally differentiate. Reconstruction of gene regulatory interactions reveals that this process is driven by a shallow network in which maternally deposited regulators activate a small set of transcription factors that co-regulate hundreds of differentiation genes. Notably, misexpression of these transcription factors in pluripotent cells is sufficient to ectopically activate their targets. This study provides a rich resource for analyzing embryonic gene regulation and reveals the regulatory logic of instant differentiation.

Single Cell Transcriptomics-Informed Induced Pluripotent Stem Cells Differentiation to Tenogenic Lineage.

During vertebrate embryogenesis, axial tendons develop from the paraxial mesoderm and differentiate through specific developmental stages to reach the syndetome stage. While the main roles of signaling pathways in the earlier stages of the differentiation have been well established, pathway nuances in syndetome specification from the sclerotome stage have yet to be explored. Here, we show stepwise differentiation of human iPSCs to the syndetome stage using chemically defined media and small molecules that were modified based on single cell RNA-sequencing and pathway analysis. We identified a significant population of branching off-target cells differentiating towards a neural phenotype overexpressing Wnt. Further transcriptomics post-addition of a WNT inhibitor at the somite stage and onwards revealed not only total removal of the neural off-target cells, but also increased syndetome induction efficiency. Fine-tuning tendon differentiation in vitro is essential to address the current challenges in developing a successful cell-based tendon therapy.

ULI-ssDRIP-seq revealed R-loop dynamics during vertebrate early embryogenesis

R-loop, a chromatin structure containing one RNA:DNA hybrid and one unpaired single-stranded DNA, plays multiple biological roles. However, due to technical limitations, the landscapes and potential functions of R-loops during embryogenesis remain elusive. Here, we developed a quantitative and high-resolution ultra-low input R-loop profiling method, named ULI-ssDRIP-seq, which can map global R-loops with as few as 1000 cells. By using ULI-ssDRIP-seq, we reveal the R-loop dynamics in the zebrafish from gametes to early embryos. In oocytes, the R-loop level is relatively low in most regions of the nuclear genome, except maternal-inherited rDNA and mitochondrial genome. The correlation between R-loop and CG methylation dynamics during early development is relatively weak. Furthermore, either up- or down-regulation of global R-loops by knockdown or overexpression of RNase H1 causes a delay of embryonic development with dramatic expression changes in zygotic and maternal genes. This study provides comprehensive R-loop landscapes during early vertebrate embryogenesis and demonstrates the implication of R-loops in embryonic development.

Systemic and strict regulation of the glutathione redox state in mitochondria and cytosol is needed for zebrafish ontogeny

BackgroundRedox control seems to be indispensable for proper embryonic development. The ratio between glutathione (GSH) and its oxidized disulfide (GSSG) is the most abundant cellular redox circuit. MethodsWe used zebrafish harboring the glutaredoxin 1-redox sensitive green fluorescent protein (Grx1-roGFP) probe either in mitochondria or cytosol to test the hypothesis that the GSH:GSSG ratio is strictly regulated through zebrafish embryogenesis to sustain the different developmental processes of the embryo. ResultsFollowing the GSSG:GSH ratio as a proxy for the GSH-dependent reduction potential (EhGSH) revealed increasing mitochondrial and cytosolic EhGSH during cleavage and gastrulation. During organogenesis, cytosolic EhGSH decreased, while that of mitochondria remained high. The similarity between EhGSH in brain and muscle suggests a central regulation. Modulation of GSH metabolism had only modest effects on the GSSG:GSH ratios of newly hatched larvae. However, inhibition of GSH reductase directly after fertilization led to dead embryos already 10 h later. Exposure to the emerging environmental pollutant Perfluorooctane Sulfonate (PFOS) disturbed the apparent regulated EhGSH as well. ConclusionsMitochondrial and cytosolic GSSG:GSH ratios are almost identical in different organs during zebrafish development indicating that the EhGSH might follow H2O2 levels and rather indirectly affect specific enzymatic activities needed for proper embryogenesis. General significanceOur data confirm that vertebrate embryogenesis depends on strictly regulated redox homeostasis. Disturbance of the GSSG:GSH circuit, e.g. induced by environmental pollution, leads to malformation and death.

Dissection of N-deacetylase and N-sulfotransferase activities of NDST1 and their effects on Wnt8 distribution and signaling in Xenopus embryos.

Wnt is a family of secreted signaling proteins involved in the regulation of cellular processes, including maintenance of stem cells, carcinogenesis, and cell differentiation. In the context of early vertebrate embryogenesis, graded distribution of Wnt proteins has been thought to regulate positional information along the antero-posterior axis. However, understanding of the molecular basis for Wnt spatial distribution remains poor. Modified states of heparan sulfate (HS) proteoglycans are essential for Wnt8 localization, because depletion of N-deacetylase/N-sulfotransferase 1 (NDST1), a modification enzyme of HS chains, decreases Wnt8 levels and NDST1 overexpression increases Wnt8 levels on the cell surface. Since overexpression of NDST1 increases both deacetylation and N-sulfation of HS chains, it is not clear which function of NDST1 is actually involved in Wnt8 localization. In the present study, we generated an NDST1 mutant that specifically increases deacetylation, but not N-sulfation, of HS chains in Xenopus embryos. Unlike wild-type NDST1, this mutant did not increase Wnt8 accumulation on the cell surface, but it reduced canonical Wnt signaling, as determined with the TOP-Flash reporter assay. These results suggest that N-sulfation of HS chains is responsible for localization of Wnt8 and Wnt8 signaling, whereas deacetylation has an inhibitory effect on canonical Wnt signaling. Consistently, overexpression of wild-type NDST1, but not the mutant, resulted in small eyes in Xenopus embryos. Thus, our NDST1 mutant enables us to dissect the regulation of Wnt8 localization and signaling by HS proteoglycans by specifically manipulating the enzymatic activities of NDST1.

Live Imaging Transverse Sections of Zebrafish Embryo Explants.

Vertebrate embryogenesis is a highly dynamic process involving coordinated cell and tissue movements that generate the final embryonic body plan. Many of these movements are difficult to image at high resolution because they occur deep within the embryo along the midline, causing light scattering and requiring longer working distances. Here, we present an explant-based method to image transverse cross sections of living zebrafish embryos. This method allows for the capture of all cell movements at high-resolution throughout the embryonic trunk, including hard-to-image deep tissues. This technique offers an alternative to expensive or computationally difficult microscopy methods. Key features • Generates intact zebrafish explants with minimal tissue disturbance. • Allows for live imaging of deep tissues normally obscured by common confocal microscopy techniques. • Immobilizes tissues for extended periods required for time-lapse imaging. • Utilizes readily available reagents and tools, which can minimize the time and cost of the procedure.

Single-nucleus chromatin landscapes during zebrafish early embryogenesis

Vertebrate embryogenesis is a remarkable process, during which numerous cell types of different lineages arise within a short time frame. An overwhelming challenge to understand this process is the lack of dynamic chromatin accessibility information to correlate cis-regulatory elements (CREs) and gene expression within the hierarchy of cell fate decisions. Here, we employed single-nucleus ATAC-seq to generate a chromatin accessibility dataset on the first day of zebrafish embryogenesis, including 3.3 hpf, 5.25 hpf, 6 hpf, 10 hpf, 12 hpf, 18 hpf and 24 hpf, obtained 51,620 high-quality nuclei and 23 clusters. Furthermore, by integrating snATAC-seq data with single-cell RNA-seq data, we described the dynamics of chromatin accessibility and gene expression across developmental time points, which validates the accuracy of the chromatin landscape data. Together, our data could serve as a fundamental resource for revealing the epigenetic regulatory mechanisms of zebrafish embryogenesis.

Transmembrane protein 150b attenuates BMP signaling in the Xenopus organizer.

The vertebrate organizer is a specified embryonic tissue that regulates dorsoventral patterning and axis formation. Although numerous cellular signaling pathways have been identified as regulators of the organizer's dynamic functions, the process remains incompletely understood, and as-yet unknown pathways remain to be explored for sophisticated mechanistic understanding of the vertebrate organizer. To identify new potential key factors of the organizer, we performed complementary DNA(cDNA) microarray screening using organizer-mimicking Xenopus laevis tissue. This analysis yielded a list of prospective organizer genes, and we determined the role of six-transmembrane domain containing transmembrane protein 150b (Tmem150b) in organizer function. Tmem150b was expressed in the organizer region and induced by Activin/Nodal signaling. In X. laevis, Tmem150b knockdown resulted in head defects and a shortened body axis. Moreover, Tmem150b negatively regulated bone morphogenetic protein (BMP) signaling, likely via physical interaction with activin receptor-like kinase 2 (ALK2). These findings demonstrated that Tmem150b functions as a novel membrane regulatory factor of BMP signaling with antagonistic effects, contributing to the understanding of regulatory molecular mechanisms of organizer axis function. Investigation of additional candidate genes identified in the cDNA microarray analysis could further delineate the genetic networks of the organizer during vertebrate embryogenesis.

Retinoic acid signaling pathway perturbation impacts mesodermal-tissue development in the zebrafish embryo: Biomarker candidate identification using transcriptomics.

The zebrafish embryo (ZE) model provides a developmental model well conserved throughout vertebrate embryogenesis, with relevance for early human embryo development. It was employed to search for gene expression biomarkers of compound-induced disruption of mesodermal development. We were particularly interested in the expression of genes related to the retinoic acid signaling pathway (RA-SP), as a major morphogenetic regulating mechanism. We exposed ZE to teratogenic concentrations of valproic acid (VPA) and all-trans retinoic acid (ATRA), using folic acid (FA) as a non-teratogenic control compound shortly after fertilization for 4 h, and performed gene expression analysis by RNA sequencing. We identified 248 genes specifically regulated by both teratogens but not by FA. Further analysis of this gene set revealed 54 GO-terms related to the development of mesodermal tissues, distributed along the paraxial, intermediate, and lateral plate sections of the mesoderm. Gene expression regulation was specific to tissues and was observed for somites, striated muscle, bone, kidney, circulatory system, and blood. Stitch analysis revealed 47 regulated genes related to the RA-SP, which were differentially expressed in the various mesodermal tissues. These genes provide potential molecular biomarkers of mesodermal tissue and organ (mal)formation in the early vertebrate embryo.

Histone deacetylase 1 maintains lineage integrity through histone acetylome refinement during early embryogenesis.

Histone acetylation is a pivotal epigenetic modification that controls chromatin structure and regulates gene expression. It plays an essential role in modulating zygotic transcription and cell lineage specification of developing embryos. While the outcomes of many inductive signals have been described to require enzymatic activities of histone acetyltransferases and deacetylases (HDACs), the mechanisms by which HDACs confine the utilization of the zygotic genome remain to be elucidated. Here, we show that histone deacetylase 1 (Hdac1) progressively binds to the zygotic genome from mid blastula and onward. The recruitment of Hdac1 to the genome at blastula is instructed maternally. Cis-regulatory modules (CRMs) bound by Hdac1 possess epigenetic signatures underlying distinct functions. We highlight a dual function model of Hdac1 where Hdac1 not only represses gene expression by sustaining a histone hypoacetylation state on inactive chromatin, but also maintains gene expression through participating in dynamic histone acetylation-deacetylation cycles on active chromatin. As a result, Hdac1 maintains differential histone acetylation states of bound CRMs between different germ layers and reinforces the transcriptional program underlying cell lineage identities, both in time and space. Taken together, our study reveals a comprehensive role for Hdac1 during early vertebrate embryogenesis.

Rgp1 contributes to craniofacial cartilage development and Rab8a-mediated collagen II secretion.

Rgp1 was previously identified as a component of a guanine nucleotide exchange factor (GEF) complex to activate Rab6a-mediated trafficking events in and around the Golgi. While the role of Rgp1 in protein trafficking has been examined in vitro and in yeast, the role of Rgp1 during vertebrate embryogenesis and protein trafficking in vivo is unknown. Using genetic, CRISPR-induced zebrafish mutants for Rgp1 loss-of-function, we found that Rgp1 is required for craniofacial cartilage development. Within live rgp1-/- craniofacial chondrocytes, we observed altered movements of Rab6a+ vesicular compartments, consistent with a conserved mechanism described in vitro. Using transmission electron microscopy (TEM) and immunofluorescence analyses, we show that Rgp1 plays a role in the secretion of collagen II, the most abundant protein in cartilage. Our overexpression experiments revealed that Rab8a is a part of the post-Golgi collagen II trafficking pathway. Following loss of Rgp1, chondrocytes activate an Arf4b-mediated stress response and subsequently respond with nuclear DNA fragmentation and cell death. We propose that an Rgp1-regulated Rab6a-Rab8a pathway directs secretion of ECM cargoes such as collagen II, a pathway that may also be utilized in other tissues where coordinated trafficking and secretion of collagens and other large cargoes is required for normal development and tissue function.

Cpeb1b-mediated cytoplasmic polyadenylation of shha mRNA modulates zebrafish definitive hematopoiesis

During vertebrate embryogenesis, hematopoietic stem and progenitor cell (HSPC) production through endothelial-to-hematopoietic transition requires suitable developmental signals, but how these signals are accurately regulated remains incompletely understood. Cytoplasmic polyadenylation, which is one of the posttranscriptional regulations, plays a crucial role in RNA metabolism. Here, we report that Cpeb1b-mediated cytoplasmic polyadenylation is important for HSPC specification by translational control of Hedgehog (Hh) signaling during zebrafish early development. Cpeb1b is highly expressed in notochord and its deficiency results in defective HSPC production. Mechanistically, Cpeb1b regulates hemogenic endothelium specification by the Hedgehog-Vegf-Notch axis. We demonstrate that the cytoplasmic polyadenylation element motif-dependent interaction between Cpeb1b and shha messenger RNA (mRNA) in the liquid-like condensates, which are induced by Pabpc1b phase separation, is required for cytoplasmic polyadenylation of shha mRNA. Intriguingly, the cytoplasmic polyadenylation regulates translation but not stability of shha mRNA, which further enhances the Shha protein level and Hh signal transduction. Taken together, our findings uncover the role of Cpeb1b-mediated cytoplasmic polyadenylation in HSPC development and provide insights into how posttranscriptional regulation can direct developmental signals with high fidelity to translate them into cell fate transition.

Somitic mesoderm morphogenesis is necessary for neural tube closure during Xenopus development.

Neural tube closure is a fundamental process during vertebrate embryogenesis, which leads to the formation of the central nervous system. Defective neural tube closure leads to neural tube defects which are some of the most common human birth defects. While the intrinsic morphogenetic events shaping the neuroepithelium have been studied extensively, how tissues mechanically coupled with the neural plate influence neural tube closure remains poorly understood. Here, using Xenopus laevis embryos, live imaging in combination with loss of function experiments and morphometric analysis of fixed samples we explore the reciprocal mechanical communication between the neural plate and the somitic mesoderm and its impact on tissue morphogenesis. We show that although somitic mesoderm convergent extension occurs independently from neural plate morphogenesis neural tube closure depends on somitic mesoderm morphogenesis. Specifically, impaired somitic mesoderm remodelling results in defective apical constriction within the neuroepithelium and failure of neural tube closure. Last, our data reveal that mild abnormalities in somitic mesoderm and neural plate morphogenesis have a synergistic effect during neurulation, leading to severe neural tube closure defects. Overall, our data reveal that defective morphogenesis of tissues mechanically coupled with the neural plate can not only drastically exacerbate mild neural tube defects that may arise from abnormalities within the neural tissue but can also elicit neural tube defects even when the neural plate is itself free of inherent defects.

Nano3P-seq: transcriptome-wide analysis of gene expression and tail dynamics using end-capture nanopore cDNA sequencing

RNA polyadenylation plays a central role in RNA maturation, fate, and stability. In response to developmental cues, polyA tail lengths can vary, affecting the translation efficiency and stability of mRNAs. Here we develop Nanopore 3′ end-capture sequencing (Nano3P-seq), a method that relies on nanopore cDNA sequencing to simultaneously quantify RNA abundance, tail composition, and tail length dynamics at per-read resolution. By employing a template-switching-based sequencing protocol, Nano3P-seq can sequence RNA molecule from its 3′ end, regardless of its polyadenylation status, without the need for PCR amplification or ligation of RNA adapters. We demonstrate that Nano3P-seq provides quantitative estimates of RNA abundance and tail lengths, and captures a wide diversity of RNA biotypes. We find that, in addition to mRNA and long non-coding RNA, polyA tails can be identified in 16S mitochondrial ribosomal RNA in both mouse and zebrafish models. Moreover, we show that mRNA tail lengths are dynamically regulated during vertebrate embryogenesis at an isoform-specific level, correlating with mRNA decay. Finally, we demonstrate the ability of Nano3P-seq in capturing non-A bases within polyA tails of various lengths, and reveal their distribution during vertebrate embryogenesis. Overall, Nano3P-seq is a simple and robust method for accurately estimating transcript levels, tail lengths, and tail composition heterogeneity in individual reads, with minimal library preparation biases, both in the coding and non-coding transcriptome.