Sites of transcription initiation drive mRNA isoform selection.

The generation of distinct messenger RNA isoforms through alternative RNA processing modulates the expression and function of genes, often in a cell-type-specific manner. Here, we assess the regulatory relationships between transcription initiation, alternative splicing, and 3' end site selection. Applying long-read sequencing to accurately represent even the longest transcripts from end to end, we quantify mRNA isoforms in Drosophila tissues, including the transcriptionally complex nervous system. We find that in Drosophila heads, as well as in human cerebral organoids, 3' end site choice is globally influenced by the site of transcription initiation (TSS). "Dominant promoters," characterized by specific epigenetic signatures including p300/CBP binding, impose a transcriptional constraint to define splice and polyadenylation variants. In vivo deletion or overexpression of dominant promoters as well as p300/CBP loss disrupted the 3' end expression landscape. Our study demonstrates the crucial impact of TSS choice on the regulation of transcript diversity and tissue identity.

Interplay between the transcription preinitiation complex and the +1 nucleosome.

Eukaryotic transcription starts with the assembly of a preinitiation complex (PIC) on core promoters. Flanking this region is the +1 nucleosome, the first nucleosome downstream of the core promoter. While this nucleosome is rich in epigenetic marks and plays a key role in transcription regulation, how the +1 nucleosome interacts with the transcription machinery has been a long-standing question. Here, we summarize recent structural and functional studies of the +1 nucleosome in complex with the PIC. We specifically focus on how differently organized promoter-nucleosome templates affect the assembly of the PIC and PIC-Mediator on chromatin and result in distinct transcription initiation.

PD-1 Imposes Qualitative Control of Cellular Transcriptomes in Response to T Cell Activation.

Targeted blockade of programmed cell death 1 (PD-1), an immune-checkpoint receptor that inhibits T cell activation, provides clinical benefits in various cancers. However, how PD-1 modulates gene expression in T cells remains enigmatic. Here we investigated how PD-1 affects transcriptome changes induced by T cell receptor (TCR) activation. Intriguingly, we identified a huge variance in PD-1 sensitivity among TCR-inducible genes. When we quantified the half maximal effective concentration (EC50) as the relationship between change in gene expression and TCR signal strength, we found that genes associated with survival and proliferation were efficiently expressed upon TCR activation and resistant to PD-1-mediated inhibition. Conversely, genes encoding cytokines and effector molecules were expressed less efficiently and sensitive to PD-1-mediated inhibition. We further demonstrated that transcription factor binding motifs and CpG frequency in the promoter region affect EC50 and thus the PD-1 sensitivity of genes. Our findings explain how PD-1, dependent on the TCR signal strength, calibrates cellular transcriptomes to shape functional properties of T cell populations.

Transcription factor IID parks and drives preinitiation complexes at sharp or broad promoters.

Core promoters are sites where transcriptional regulatory inputs of a gene are integrated to direct the assembly of the preinitiation complex (PIC) and RNA polymerase II (Pol II) transcription output. Until now, core promoter functions have been investigated by distinct methods, including Pol II transcription initiation site mappings and structural characterization of PICs on distinct promoters. Here, we bring together these previously unconnected observations and hypothesize how, on metazoan TATA promoters, the precisely structured building up of transcription factor (TF) IID-based PICs results in sharp transcription start site (TSS) selection; or, in contrast, how the less strictly controlled positioning of the TATA-less promoter DNA relative to TFIID-core PIC components results in alternative broad TSS selections by Pol II.

CapZyme-Seq Comprehensively Defines Promoter-Sequence Determinants for RNA 5' Capping with NAD.

Nucleoside-containing metabolites such as NAD+ can be incorporated as 5' caps on RNA by serving as non-canonical initiating nucleotides (NCINs) for transcription initiation by RNA polymerase (RNAP). Here, we report CapZyme-seq, a high-throughput-sequencing method that employs NCIN-decapping enzymes NudC and Rai1 to detect and quantify NCIN-capped RNA. By combining CapZyme-seq with multiplexed transcriptomics, we determine efficiencies of NAD+ capping by Escherichia coli RNAP for ∼16,000 promoter sequences. The results define preferred transcription start site (TSS) positions for NAD+ capping and define a consensus promoter sequence for NAD+ capping: HRRASWW (TSS underlined). By applying CapZyme-seq to E. coli total cellular RNA, we establish that sequence determinants for NCIN capping in vivo match the NAD+-capping consensus defined in vitro, and we identify and quantify NCIN-capped small RNAs (sRNAs). Our findings define the promoter-sequence determinants for NCIN capping with NAD+ and provide a general method for analysis of NCIN capping in vitro and in vivo.

HIV-1 usurps transcription start site heterogeneity of host RNA polymerase II to maximize replication fitness.

HIV-1 relies on host RNA polymeraseII (Pol II) to transcribe its genome and uses multiple transcription start sites (TSS), including three consecutive guanosines located near the U3-R junction, to generate transcripts containing three, two, and one guanosine at the 5' end, referred to as 3G, 2G, and 1G RNA, respectively. The 1G RNA is preferentially selected for packaging, indicating that these 99.9% identical RNAs exhibit functional differences and highlighting the importance of TSS selection. Here, we demonstrate that TSS selection is regulated by sequences between the CATA/TATA box and the beginning of R. Furthermore, we have generated two HIV-1 mutants with distinct 2-nucleotide modifications that predominantly express 3G RNA or 1G RNA. Both mutants can generate infectious viruses and undergo multiple rounds of replication in T cells. However, both mutants exhibit replication defects compared to the wild-type virus. The 3G-RNA-expressing mutant displays an RNA genome-packaging defect and delayed replication kinetics, whereas the 1G-RNA-expressing mutant exhibits reduced Gag expression and a replication fitness defect. Additionally, reversion of the latter mutant is frequently observed, consistent with sequence correction by plus-strand DNA transfer during reverse transcription. These findings demonstrate that HIV-1 maximizes its replication fitness by usurping the TSS heterogeneity of host RNA Pol II to generate unspliced RNAs with different specialized roles in viral replication. The three consecutive guanosines at the junction of U3 and R may also maintain HIV-1 genome integrity during reverse transcription. These studies reveal the intricate regulation of HIV-1 RNA and complex replication strategy.

The human initiator is a distinct and abundant element that is precisely positioned in focused core promoters.

DNA sequence signals in the core promoter, such as the initiator (Inr), direct transcription initiation by RNA polymerase II. Here we show that the human Inr has the consensus of BBCA+1BW at focused promoters in which transcription initiates at a single site or a narrow cluster of sites. The analysis of 7678 focused transcription start sites revealed 40% with a perfect match to the Inr and 16% with a single mismatch outside of the CA+1 core. TATA-like sequences are underrepresented in Inr promoters. This consensus is a key component of the DNA sequence rules that specify transcription initiation in humans.

Finding the start site: redefining the human initiator element.

Transcription by RNA polymerase II (Pol II) is dictated in part by core promoter elements, which are DNA sequences flanking the transcription start site (TSS) that help direct the proper initiation of transcription. Taking advantage of recent advances in genome-wide sequencing approaches, Vo ngoc and colleagues (pp. 6-11) identified transcripts with focused sites of initiation and found that many were transcribed from promoters containing a new consensus sequence for the human initiator (Inr) core promoter element.

Common occurrence of hotspots of single strand DNA breaks at transcriptional start sites.

BACKGROUND: We recently developed two high-resolution methods for genome-wide mapping of two prominent types of DNA damage, single-strand DNA breaks (SSBs) and abasic (AP) sites and found highly complex and non-random patterns of these lesions in mammalian genomes. One salient feature of SSB and AP sites was the existence of single-nucleotide hotspots for both lesions. RESULTS: In this work, we show that SSB hotspots are enriched in the immediate vicinity of transcriptional start sites (TSSs) in multiple normal mammalian tissues, however the magnitude of enrichment varies significantly with tissue type and appears to be limited to a subset of genes. SSB hotspots around TSSs are enriched on the template strand and associate with higher expression of the corresponding genes. Interestingly, SSB hotspots appear to be at least in part generated by the base-excision repair (BER) pathway from the AP sites. CONCLUSIONS: Our results highlight complex relationship between DNA damage and regulation of gene expression and suggest an exciting possibility that SSBs at TSSs might function as sensors of DNA damage to activate genes important for DNA damage response.

Characterization of gene expression in the kidney of renal tubular cell-specific NFAT5 knockout mice.

Nuclear factor of activated T cells 5 (NFAT5; also called TonEBP/OREBP) is a transcription factor that is activated by hypertonicity and induces osmoprotective genes to protect cells against hypertonic conditions. In the kidney, renal tubular NFAT5 is known to be involved in the urine concentration mechanism. Previous studies have suggested that NFAT5 modulates the immune system and exerts various effects on organ damage, depending on organ and disease states. Pathophysiological roles of NFAT5 in renal tubular cells, however, still remain obscure. We conducted comprehensive analysis by performing transcription start site (TSS) sequencing on the kidney of inducible and renal tubular cell-specific NFAT5 knockout (KO) mice. Mice were subjected to unilateral ureteral obstruction to examine the relevance of renal tubular NFAT5 in renal fibrosis. TSS sequencing analysis identified 722 downregulated TSSs and 1,360 upregulated TSSs, which were differentially regulated ≤-1.0 and ≥1.0 in log2 fold, respectively. Those TSSs were annotated to 532 downregulated genes and 944 upregulated genes, respectively. Motif analysis showed that sequences that possibly bind to NFAT5 were enriched in TSSs of downregulated genes. Gene Ontology analysis with the upregulated genes suggested disorder of innate and adaptive immune systems in the kidney. Unilateral ureteral obstruction significantly exacerbated renal fibrosis in the renal medulla in KO mice compared with wild-type mice, accompanied by enhanced activation of immune responses. In conclusion, NFAT5 in renal tubules could have pathophysiological roles in renal fibrosis through modulating innate and adaptive immune systems in the kidney.NEW & NOTEWORTHY TSS-Seq analysis of the kidney from renal tubular cell-specific NFAT5 KO mice uncovered novel genes that are possibly regulated by NFAT5 in the kidney under physiological conditions. The study further implied disorders of innate and adaptive immune systems in NFAT5 KO mice, thereby exacerbating renal fibrosis at pathological states. Our results may implicate the involvement of renal tubular NFAT5 in the progression of renal fibrosis. Further studies would be worthwhile for the development of novel therapy to treat chronic kidney disease.

Identification of high-confidence human poly(A) RNA isoform scaffolds using nanopore sequencing.

Nanopore sequencing devices read individual RNA strands directly. This facilitates identification of exon linkages and nucleotide modifications; however, using conventional direct RNA nanopore sequencing, the 5' and 3' ends of poly(A) RNA cannot be identified unambiguously. This is due in part to RNA degradation in vivo and in vitro that can obscure transcription start and end sites. In this study, we aimed to identify individual full-length human RNA isoforms among ∼4 million nanopore poly(A)-selected RNA reads. First, to identify RNA strands bearing 5' m7G caps, we exchanged the biological cap for a modified cap attached to a 45-nt oligomer. This oligomer adaptation method improved 5' end sequencing and ensured correct identification of the 5' m7G capped ends. Second, among these 5'-capped nanopore reads, we screened for features consistent with a 3' polyadenylation site. Combining these two steps, we identified 294,107 individual high-confidence full-length RNA scaffolds from human GM12878 cells, most of which (257,721) aligned to protein-coding genes. Of these, 4876 scaffolds indicated unannotated isoforms that were often internal to longer, previously identified RNA isoforms. Orthogonal data for m7G caps and open chromatin, such as CAGE and DNase-HS seq, confirmed the validity of these high-confidence RNA scaffolds.

Selective packaging of HIV-1 RNA genome is guided by the stability of 5' untranslated region polyA stem.

To generate infectious virus, HIV-1 must package two copies of its full-length RNA into particles. HIV-1 transcription initiates from multiple, neighboring sites, generating RNA species that only differ by a few nucleotides at the 5' end, including those with one (1G) or three (3G) 5' guanosines. Strikingly, 1G RNA is preferentially packaged into virions over 3G RNA. We investigated how HIV-1 distinguishes between these nearly identical RNAs using in-gel chemical probing combined with recently developed computational tools for determining RNA conformational ensembles, as well as cell-based assays to quantify the efficiency of RNA packaging into viral particles. We found that 1G and 3G RNAs fold into distinct structural ensembles. The 1G RNA, but not the 3G RNA, primarily adopts conformations with an intact polyA stem, exposed dimerization initiation site, and multiple, unpaired guanosines known to mediate Gag binding. Furthermore, we identified mutants that exhibited altered genome selectivity and packaged 3G RNA efficiently. In these mutants, both 1G and 3G RNAs fold into similar conformational ensembles, such that they can no longer be distinguished. Our findings demonstrate that polyA stem stability guides RNA-packaging selectivity. These studies also uncover the mechanism by which HIV-1 selects its genome for packaging: 1G RNA is preferentially packaged because it exposes structural elements that promote RNA dimerization and Gag binding.

Systematic mutagenesis of TFIIH subunit p52/Tfb2 identifies residues required for XPB/Ssl2 subunit function and genetic interactions with TFB6.

TFIIH is an evolutionarily conserved complex that plays central roles in both RNA polymerase II (pol II) transcription and DNA repair. As an integral component of the pol II preinitiation complex, TFIIH regulates pol II enzyme activity in numerous ways. The TFIIH subunit XPB/Ssl2 is an ATP-dependent DNA translocase that stimulates promoter opening prior to transcription initiation. Crosslinking-mass spectrometry and cryo-EM results have shown a conserved interaction network involving XPB/Ssl2 and the C-terminal Hub region of the TFIIH p52/Tfb2 subunit, but the functional significance of specific residues is unclear. Here, we systematically mutagenized the HubA region of Tfb2 and screened for growth phenotypes in a TFB6 deletion background in Saccharomyces cerevisiae. We identified six lethal and 12 conditional mutants. Slow growth phenotypes of all but three conditional mutants were relieved in the presence of TFB6, thus identifying a functional interaction between Tfb2 HubA mutants and Tfb6, a protein that dissociates Ssl2 from TFIIH. Our biochemical analysis of Tfb2 mutants with severe growth phenotypes revealed defects in Ssl2 association, with similar results in human cells. Further characterization of these tfb2 mutant cells revealed defects in GAL gene induction, and reduced occupancy of TFIIH and pol II at GAL gene promoters, suggesting that functionally competent TFIIH is required for proper pol II recruitment to preinitiation complexes in vivo. Consistent with recent structural models of TFIIH, our results identify key residues in the p52/Tfb2 HubA domain that are required for stable incorporation of XPB/Ssl2 into TFIIH and for pol II transcription.

Simple and accurate transcriptional start site identification using Smar2C2 and examination of conserved promoter features.

The precise and accurate identification and quantification of transcriptional start sites (TSSs) is key to understanding the control of transcription. The core promoter consists of the TSS and proximal non-coding sequences, which are critical in transcriptional regulation. Therefore, the accurate identification of TSSs is important for understanding the molecular regulation of transcription. Existing protocols for TSS identification are challenging and expensive, leaving high-quality data available for a small subset of organisms. This sparsity of data impairs study of TSS usage across tissues or in an evolutionary context. To address these shortcomings, we developed Smart-Seq2 Rolling Circle to Concatemeric Consensus (Smar2C2), which identifies and quantifies TSSs and transcription termination sites. Smar2C2 incorporates unique molecular identifiers that allowed for the identification of as many as 70 million sites, with no known upper limit. We have also generated TSS data sets from as little as 40 pg of total RNA, which was the smallest input tested. In this study, we used Smar2C2 to identify TSSs in Glycine max (soybean), Oryza sativa (rice), Sorghum bicolor (sorghum), Triticum aestivum (wheat) and Zea mays (maize) across multiple tissues. This wide panel of plant TSSs facilitated the identification of evolutionarily conserved features, such as novel patterns in the dinucleotides that compose the initiator element (Inr), that correlated with promoter expression levels across all species examined. We also discovered sequence variations in known promoter motifs that are positioned reliably close to the TSS, such as differences in the TATA box and in the Inr that may prove significant to our understanding and control of transcription initiation. Smar2C2 allows for the easy study of these critical sequences, providing a tool to facilitate discovery.

TSSFinder-fast and accurate ab initio prediction of the core promoter in eukaryotic genomes.

Promoter annotation is an important task in the analysis of a genome. One of the main challenges for this task is locating the border between the promoter region and the transcribing region of the gene, the transcription start site (TSS). The TSS is the reference point to delimit the DNA sequence responsible for the assembly of the transcribing complex. As the same gene can have more than one TSS, so to delimit the promoter region, it is important to locate the closest TSS to the site of the beginning of the translation. This paper presents TSSFinder, a new software for the prediction of the TSS signal of eukaryotic genes that is significantly more accurate than other available software. We currently are the only application to offer pre-trained models for six different eukaryotic organisms: Arabidopsis thaliana, Drosophila melanogaster, Gallus gallus, Homo sapiens, Oryza sativa and Saccharomyces cerevisiae. Additionally, our software can be easily customized for specific organisms using only 125 DNA sequences with a validated TSS signal and corresponding genomic locations as a training set. TSSFinder is a valuable new tool for the annotation of genomes. TSSFinder source code and docker container can be downloaded from http://tssfinder.github.io. Alternatively, TSSFinder is also available as a web service at http://sucest-fun.org/wsapp/tssfinder/.

The single-cell landscape of alternative transcription start sites of diabetic retina.

More than half of mammalian protein-coding genes have multiple transcription start sites. Alternative transcription start site (TSS) modulate mRNA stability, localization, and translation efficiency on post-transcription level, and even generate novel protein isoforms. However, differential TSS usage among cell types in healthy and diabetic retina remains poorly characterized. In this study, by using 5'-tag-based single-cell RNA sequencing, we identified cell type-specific alternative TSS events and key transcription factors for each of retinal cell types. We observed that lengthening of 5'- UTRs in retinal cell types are enriched for multiple RNA binding protein binding sites, including splicing regulators Rbfox1/2/3 and Nova1. Furthermore, by comparing TSS expression between healthy and diabetic retina, we identified elevated apoptosis signal in Müller glia and microglia, which can be served as a putative early sign of diabetic retinopathy. By measuring 5'UTR isoforms in retinal single-cell dataset, our work provides a comprehensive panorama of alternative TSS and its potential consequence related to post-transcriptional regulation. We anticipate our assay can not only provide insights into cellular heterogeneity driven by transcriptional initiation, but also open up the perspectives for identification of novel diagnostic indexes for diabetic retinopathy.

Alternative transcription start sites contribute to acute-stress-induced transcriptome response in human skeletal muscle.

BACKGROUND: More than half of human protein-coding genes have an alternative transcription start site (TSS). We aimed to investigate the contribution of alternative TSSs to the acute-stress-induced transcriptome response in human tissue (skeletal muscle) using the cap analysis of gene expression approach. TSSs were examined at baseline and during recovery after acute stress (a cycling exercise). RESULTS: We identified 44,680 CAGE TSS clusters (including 3764 first defined) belonging to 12,268 genes and annotated for the first time 290 TSSs belonging to 163 genes. The transcriptome dynamically changes during the first hours after acute stress; the change in the expression of 10% of genes was associated with the activation of alternative TSSs, indicating differential TSSs usage. The majority of the alternative TSSs do not increase proteome complexity suggesting that the function of thousands of alternative TSSs is associated with the fine regulation of mRNA isoform expression from a gene due to the transcription factor-specific activation of various alternative TSSs. We identified individual muscle promoter regions for each TSS using muscle open chromatin data (ATAC-seq and DNase-seq). Then, using the positional weight matrix approach we predicted time course activation of "classic" transcription factors involved in response of skeletal muscle to contractile activity, as well as diversity of less/un-investigated factors. CONCLUSIONS: Transcriptome response induced by acute stress related to activation of the alternative TSSs indicates that differential TSSs usage is an essential mechanism of fine regulation of gene response to stress stimulus. A comprehensive resource of accurate TSSs and individual promoter regions for each TSS in muscle was created. This resource together with the positional weight matrix approach can be used to accurate prediction of TFs in any gene(s) of interest involved in the response to various stimuli, interventions or pathological conditions in human skeletal muscle.

Transcription start site mapping and small RNA profiling of Leptospira biflexa serovar Patoc.

Leptospirosis is an emerging zoonotic disease caused by bacterial species of the genus Leptospira. However, the regulatory mechanisms and pathways underlying the adaptation of pathogenic and non-pathogenic Leptospira spp. in different environmental conditions remain elusive. Leptospira biflexa is a non-pathogenic species of Leptospira that lives exclusively in a natural environment. It is an ideal model not only for exploring molecular mechanisms underlying the environmental survival of Leptospira species but also for identifying virulence factors unique to Leptospira's pathogenic species. In this study, we aim to establish the transcription start site (TSS) landscape and the small RNA (sRNA) profile of L. biflexa serovar Patoc grown to exponential and stationary phases via differential RNA-seq (dRNA-seq) and small RNA-seq (sRNA-seq) analyses, respectively. Our dRNA-seq analysis uncovered a total of 2726 TSSs, which are also used to identify other elements, e.g., promoter and untranslated regions (UTRs). Besides, our sRNA-seq analysis revealed a total of 603 sRNA candidates, comprising 16 promoter-associated sRNAs, 184 5'UTR-derived sRNAs, 230 true intergenic sRNAs, 136 5'UTR-antisense sRNAs, and 130 open reading frame (ORF)-antisense sRNAs. In summary, these findings reflect the transcriptional complexity of L. biflexa serovar Patoc under different growth conditions and help to facilitate our understanding of regulatory networks in L. biflexa. To the best of our knowledge, this is the first study reporting the TSS landscape of L. biflexa. The TSS and sRNA landscapes of L. biflexa can also be compared with its pathogenic counterparts, e.g., L. borgpetersenii and L. interrogans, to identify features contributing to their environmental survival and virulence.

iPromoter-Seqvec: identifying promoters using bidirectional long short-term memory and sequence-embedded features.

BACKGROUND: Promoters, non-coding DNA sequences located at upstream regions of the transcription start site of genes/gene clusters, are essential regulatory elements for the initiation and regulation of transcriptional processes. Furthermore, identifying promoters in DNA sequences and genomes significantly contributes to discovering entire structures of genes of interest. Therefore, exploration of promoter regions is one of the most imperative topics in molecular genetics and biology. Besides experimental techniques, computational methods have been developed to predict promoters. In this study, we propose iPromoter-Seqvec - an efficient computational model to predict TATA and non-TATA promoters in human and mouse genomes using bidirectional long short-term memory neural networks in combination with sequence-embedded features extracted from input sequences. The promoter and non-promoter sequences were retrieved from the Eukaryotic Promoter database and then were refined to create four benchmark datasets. RESULTS: The area under the receiver operating characteristic curve (AUCROC) and the area under the precision-recall curve (AUCPR) were used as two key metrics to evaluate model performance. Results on independent test sets showed that iPromoter-Seqvec outperformed other state-of-the-art methods with AUCROC values ranging from 0.85 to 0.99 and AUCPR values ranging from 0.86 to 0.99. Models predicting TATA promoters in both species had slightly higher predictive power compared to those predicting non-TATA promoters. With a novel idea of constructing artificial non-promoter sequences based on promoter sequences, our models were able to learn highly specific characteristics discriminating promoters from non-promoters to improve predictive efficiency. CONCLUSIONS: iPromoter-Seqvec is a stable and robust model for predicting both TATA and non-TATA promoters in human and mouse genomes. Our proposed method was also deployed as an online web server with a user-friendly interface to support research communities. Links to our source codes and web server are available at https://github.com/mldlproject/2022-iPromoter-Seqvec .

Kozak Sequence Acts as a Negative Regulator for De Novo Transcription Initiation of Newborn Coding Sequences in the Plant Genome.

The manner in which newborn coding sequences and their transcriptional competency emerge during the process of gene evolution remains unclear. Here, we experimentally simulated eukaryotic gene origination processes by mimicking horizontal gene transfer events in the plant genome. We mapped the precise position of the transcription start sites (TSSs) of hundreds of newly introduced promoterless firefly luciferase (LUC) coding sequences in the genome of Arabidopsis thaliana cultured cells. The systematic characterization of the LUC-TSSs revealed that 80% of them occurred under the influence of endogenous promoters, while the remainder underwent de novo activation in the intergenic regions, starting from pyrimidine-purine dinucleotides. These de novo TSSs obeyed unexpected rules; they predominantly occurred ∼100 bp upstream of the LUC inserts and did not overlap with Kozak-containing putative open reading frames (ORFs). These features were the output of the immediate responses to the sequence insertions, rather than a bias in the screening of the LUC gene function. Regarding the wild-type genic TSSs, they appeared to have evolved to lack any ORFs in their vicinities. Therefore, the repulsion by the de novo TSSs of Kozak-containing ORFs described above might be the first selection gate for the occurrence and evolution of TSSs in the plant genome. Based on these results, we characterized the de novo type of TSS identified in the plant genome and discuss its significance in genome evolution.