Divergent sensory and immune gene evolution in sea turtles with contrasting demographic and life histories.

Sea turtles represent an ancient lineage of marine vertebrates that evolved from terrestrial ancestors over 100 Mya. The genomic basis of the unique physiological and ecological traits enabling these species to thrive in diverse marine habitats remains largely unknown. Additionally, many populations have drastically declined due to anthropogenic activities over the past two centuries, and their recovery is a high global conservation priority. We generated and analyzed high-quality reference genomes for the leatherback (Dermochelys coriacea) and green (Chelonia mydas) turtles, representing the two extant sea turtle families. These genomes are highly syntenic and homologous, but localized regions of noncollinearity were associated with higher copy numbers of immune, zinc-finger, and olfactory receptor (OR) genes in green turtles, with ORs related to waterborne odorants greatly expanded in green turtles. Our findings suggest that divergent evolution of these key gene families may underlie immunological and sensory adaptations assisting navigation, occupancy of neritic versus pelagic environments, and diet specialization. Reduced collinearity was especially prevalent in microchromosomes, with greater gene content, heterozygosity, and genetic distances between species, supporting their critical role in vertebrate evolutionary adaptation. Finally, diversity and demographic histories starkly contrasted between species, indicating that leatherback turtles have had a low yet stable effective population size, exhibit extremely low diversity compared with other reptiles, and harbor a higher genetic load compared with green turtles, reinforcing concern over their persistence under future climate scenarios. These genomes provide invaluable resources for advancing our understanding of evolution and conservation best practices in an imperiled vertebrate lineage.

An evolutionary view of vanillylamine synthase pAMT, a key enzyme of capsaicinoid biosynthesis pathway in chili pepper.

Pungent capsaicinoid is synthesized only in chili pepper (Capsicum spp.). The production of vanillylamine from vanillin is a unique reaction in the capsaicinoid biosynthesis pathway. Although putative aminotransferase (pAMT) has been isolated as the vanillylamine synthase gene, it is unclear how Capsicum acquired pAMT. Here, we present a phylogenetic overview of pAMT and its homologs. The Capsicum genome contained 5 homologs, including pAMT, CaGABA-T1, CaGABA-T3, and two pseudogenes. Phylogenetic analysis indicated that pAMT is a member of the Solanaceae cytoplasmic GABA-Ts. Comparative genome analysis found that multiple copies of GABA-T exist in a specific Solanaceae genomic region, and the cytoplasmic GABA-Ts other than pAMT are located in the region. The cytoplasmic GABA-T was phylogenetically close to pseudo-GABA-T harboring a plastid transit peptide (pseudo-GABA-T3). This suggested that Solanaceae cytoplasmic GABA-Ts occurred via duplication of a chloroplastic GABA-T ancestor and subsequent loss of the plastid transit signal. The cytoplasmic GABA-T may have been translocated from the specific Solanaceae genomic region during Capsicum divergence, resulting in the current pAMT locus. A recombinant protein assay demonstrated that pAMT had higher vanillylamine synthase activity than those of other plant GABA-Ts. pAMT was expressed exclusively in the placental septum of mature green fruit, whereas tomato orthologs SlGABA-T2/4 exhibit a ubiquitous expression pattern in plants. These findings suggested that both the increased catalytic efficiency and transcriptional changes in pAMT may have contributed to establish vanillylamine synthesis in the capsaicinoid biosynthesis pathway. This study provides insights into the establishment of pungency in the evolution of chili peppers.

Out of the testis: biological impacts of new genes.

The study of newly evolved genes has long fascinated biologists, but large-scale studies of their expression dynamics and molecular function have provided conflicting interpretations of their biological impact. In this issue of Genes & Development, Kondo and colleagues (pp. 1841-1846) use extensive transcriptomic resources and current CRISPR/Cas9 technology to re-examine the functional impact of newly evolved genes in Drosophila and find evidence of their biological impact on male reproduction.

New genes often acquire male-specific functions but rarely become essential in Drosophila.

Relatively little is known about the in vivo functions of newly emerging genes, especially in metazoans. Although prior RNAi studies reported prevalent lethality among young gene knockdowns, our phylogenomic analyses reveal that young Drosophila genes are frequently restricted to the nonessential male reproductive system. We performed large-scale CRISPR/Cas9 mutagenesis of "conserved, essential" and "young, RNAi-lethal" genes and broadly confirmed the lethality of the former but the viability of the latter. Nevertheless, certain young gene mutants exhibit defective spermatogenesis and/or male sterility. Moreover, we detected widespread signatures of positive selection on young male-biased genes. Thus, young genes have a preferential impact on male reproductive system function.

Genome-wide identification and evolution of the tubulin gene family in Camelina sativa.

BACKGROUND: Tubulins play crucial roles in numerous fundamental processes of plant development. In flowering plants, tubulins are grouped into α-, ß- and γ-subfamilies, while α- and ß-tubulins possess a large isotype diversity and gene number variations among different species. This circumstance leads to insufficient recognition of orthologous isotypes and significantly complicates extrapolation of obtained experimental results, and brings difficulties for the identification of particular tubulin isotype function. The aim of this research is to identify and characterize tubulins of an emerging biofuel crop Camelina sativa. RESULTS: We report comprehensive identification and characterization of tubulin gene family in C. sativa, including analyses of exon-intron organization, duplicated genes comparison, proper isotype designation, phylogenetic analysis, and expression patterns in different tissues. 17 α-, 34 ß- and 6 γ-tubulin genes were identified and assigned to a particular isotype. Recognition of orthologous tubulin isotypes was cross-referred, involving data of phylogeny, synteny analyses and genes allocation on reconstructed genomic blocks of Ancestral Crucifer Karyotype. An investigation of expression patterns of tubulin homeologs revealed the predominant role of N6 (A) and N7 (B) subgenomes in tubulin expression at various developmental stages, contrarily to general the dominance of transcripts of H7 (C) subgenome. CONCLUSIONS: For the first time a complete set of tubulin gene family members was identified and characterized for allohexaploid C. sativa species. The study demonstrates the comprehensive approach of precise inferring gene orthology. The applied technique allowed not only identifying C. sativa tubulin orthologs in model Arabidopsis species and tracking tubulin gene evolution, but also uncovered that A. thaliana is missing orthologs for several particular isotypes of α- and ß-tubulins.

The role of Na+-coupled bicarbonate transporters (NCBT) in health and disease.

Cellular and organism survival depends upon the regulation of pH, which is regulated by highly specialized cell membrane transporters, the solute carriers (SLC) (For a comprehensive list of the solute carrier family members, see: https://www.bioparadigms.org/slc/ ). The SLC4 family of bicarbonate (HCO3-) transporters consists of ten members, sorted by their coupling to either sodium (NBCe1, NBCe2, NBCn1, NBCn2, NDCBE), chloride (AE1, AE2, AE3), or borate (BTR1). The ionic coupling of SLC4A9 (AE4) remains controversial. These SLC4 bicarbonate transporters may be controlled by cellular ionic gradients, cellular membrane voltage, and signaling molecules to maintain critical cellular and systemic pH (acid-base) balance. There are profound consequences when blood pH deviates even a small amount outside the normal range (7.35-7.45). Chiefly, Na+-coupled bicarbonate transporters (NCBT) control intracellular pH in nearly every living cell, maintaining the biological pH required for life. Additionally, NCBTs have important roles to regulate cell volume and maintain salt balance as well as absorption and secretion of acid-base equivalents. Due to their varied tissue expression, NCBTs have roles in pathophysiology, which become apparent in physiologic responses when their expression is reduced or genetically deleted. Variations in physiological pH are seen in a wide variety of conditions, from canonically acid-base related conditions to pathologies not necessarily associated with acid-base dysfunction such as cancer, glaucoma, or various neurological diseases. The membranous location of the SLC4 transporters as well as recent advances in discovering their structural biology makes them accessible and attractive as a druggable target in a disease context. The role of sodium-coupled bicarbonate transporters in such a large array of conditions illustrates the potential of treating a wide range of disease states by modifying function of these transporters, whether that be through inhibition or enhancement.

Genome-wide characterization of COMT family and regulatory role of CsCOMT19 in melatonin synthesis in Camellia sinensis.

BACKGROUND: Caffeic acid O-methyltransferase (COMT) is a key enzyme that regulates melatonin synthesis and is involved in regulating the growth, development, and response to abiotic stress in plants. Tea plant is a popular beverage consumed worldwide, has been used for centuries for its medicinal properties, including its ability to reduce inflammation, improve digestion, and boost immune function. By analyzing genetic variation within the COMT family, while helping tea plants resist adversity, it is also possible to gain a deeper understanding of how different tea varieties produce and metabolize catechins, then be used to develop new tea cultivars with desired flavor profiles and health benefits. RESULTS: In this study, a total of 25 CsCOMT genes were identified based on the high-quality tea (Camellia sinensis) plant genome database. Phylogenetic tree analysis of CsCOMTs with COMTs from other species showed that COMTs divided into four subfamilies (Class I, II, III, IV), and CsCOMTs was distributed in Class I, Class II, Class III. CsCOMTs not only undergoes large-scale gene recombination in pairs internally in tea plant, but also shares 2 and 7 collinear genes with Arabidopsis thaliana and poplar (Populus trichocarpa), respectively. The promoter region of CsCOMTs was found to be rich in cis-acting elements associated with plant growth and stress response. By analyzing the previously transcriptome data, it was found that some members of CsCOMT family exhibited significant tissue-specific expression and differential expression under different stress treatments. Subsequently, we selected six CsCOMTs to further validated their expression levels in different tissues organ using qRT-PCR. In addition, we silenced the CsCOMT19 through virus-induced gene silencing (VIGS) method and found that CsCOMT19 positively regulates the synthesis of melatonin in tea plant. CONCLUSION: These results will contribute to the understanding the functions of CsCOMT gene family and provide valuable information for further research on the role of CsCOMT genes in regulating tea plant growth, development, and response to abiotic stress.

Differentiation in detoxification gene complements, including neofunctionalization of duplicated cytochrome P450 genes, between lineages of cotton bollworm, Helicoverpa armigera.

Here we investigate the evolutionary dynamics of five enzyme superfamilies (CYPs, GSTs, UGTs, CCEs and ABCs) involved in detoxification in Helicoverpa armigera. The reference assembly for an African isolate of the major lineages, H. a. armigera, has 373 genes in the five superfamilies. Most of its CYPs, GSTs, UGTs and CCEs and a few of its ABCs occur in blocks and most of the clustered genes are in subfamilies specifically implicated in detoxification. Most of the genes have orthologues in the reference genome for the Oceania lineage, H. a. conferta. However, clustered orthologues and subfamilies specifically implicated in detoxification show greater sequence divergence and less constraint on non-synonymous differences between the two assemblies than do other members of the five superfamilies. Two duplicated CYPs, which were found in the H. a. armigera but not H. a. conferta reference genome, were also missing in 16 Chinese populations spanning two different lineages of H. a. armigera. The enzyme produced by one of these duplicates has higher activity against esfenvalerate than a previously described chimeric CYP mutant conferring pyrethroid resistance. Various transposable elements were found in the introns of most detoxification genes, generating diverse gene structures. Extensive resequencing data for the Chinese H. a. armigera and H. a. conferta lineages also revealed complex copy number polymorphisms in 17 CCE001s in a cluster also implicated in pyrethroid metabolism, with substantial haplotype differences between all three lineages. Our results suggest that cotton bollworm has a versatile complement of detoxification genes which are evolving in diverse ways across its range.

Range-limited Heaps' law for functional DNA words in the human genome.

Heaps' or Herdan-Heaps' law is a linguistic law describing the relationship between the vocabulary/dictionary size (type) and word counts (token) to be a power-law function. Its existence in genomes with certain definition of DNA words is unclear partly because the dictionary size in genome could be much smaller than that in a human language. We define a DNA word as a coding region in a genome that codes for a protein domain. Using human chromosomes and chromosome arms as individual samples, we establish the existence of Heaps' law in the human genome within limited range. Our definition of words in a genomic or proteomic context is different from other definitions such as over-represented k-mers which are much shorter in length. Although an approximate power-law distribution of protein domain sizes due to gene duplication and the related Zipf's law is well known, their translation to the Heaps' law in DNA words is not automatic. Several other animal genomes are shown herein also to exhibit range-limited Heaps' law with our definition of DNA words, though with various exponents. When tokens were randomly sampled and sample sizes reach to the maximum level, a deviation from the Heaps' law was observed, but a quadratic regression in log-log type-token plot fits the data perfectly. Investigation of type-token plot and its regression coefficients could provide an alternative narrative of reusage and redundancy of protein domains as well as creation of new protein domains from a linguistic perspective.

Genome-Wide Identification and Expression Analysis of the STAT Family in Reeve's Turtle (Mauremys reevesii).

The Stat (signal transducer and activator of transcription) gene family plays a vital role in regulating immunity and the processes of cellular proliferation, differentiation, and apoptosis across diverse organisms. Although the functions of Stat genes in immunity have been extensively documented in many mammals, limited data are available for reptiles. We used phylogenetic analysis to identify eight putative members of the Stat family (Stat1-1, Stat1-2, Stat2, Stat3, Stat4, Stat5b, Stat6-1, and Stat6-2) within the genome of M. reevesii, a freshwater turtle found in East Asia. Sequence analysis showed that the Stat genes contain four conserved structural domains protein interaction domain, coiled-coil domain, DNA-binding domain, and Src homology domain 2. In addition, Stat1, Stat2, and Stat6 contain TAZ2bind, Apolipo_F, and TALPID3 structural domains. The mRNA levels of Stat genes were upregulated in spleen tissues at 4, 8, 12, and 16 h after administration of lipopolysaccharide, a potent activator of the immune system. Stat5b expression at 12-h LPS post-injection exhibited the most substantial difference from the control. The expression of Stat5b in spleen tissue cellular was verified by immunofluorescence. These results suggest that Stat5b plays a role in the immune response of M. reevesii and may prove to be as a positive marker of an immune response in future studies.

Prevalent bee venom genes evolved before the aculeate stinger and eusociality.

BACKGROUND: Venoms, which have evolved numerous times in animals, are ideal models of convergent trait evolution. However, detailed genomic studies of toxin-encoding genes exist for only a few animal groups. The hyper-diverse hymenopteran insects are the most speciose venomous clade, but investigation of the origin of their venom genes has been largely neglected. RESULTS: Utilizing a combination of genomic and proteo-transcriptomic data, we investigated the origin of 11 toxin genes in 29 published and 3 new hymenopteran genomes and compiled an up-to-date list of prevalent bee venom proteins. Observed patterns indicate that bee venom genes predominantly originate through single gene co-option with gene duplication contributing to subsequent diversification. CONCLUSIONS: Most Hymenoptera venom genes are shared by all members of the clade and only melittin and the new venom protein family anthophilin1 appear unique to the bee lineage. Most venom proteins thus predate the mega-radiation of hymenopterans and the evolution of the aculeate stinger.

Evolution of a Human-Specific De Novo Open Reading Frame and Its Linked Transcriptional Silencer.

In the human genome, two short open reading frames (ORFs) separated by a transcriptional silencer and a small intervening sequence stem from the gene SMIM45. The two ORFs show different translational characteristics, and they also show divergent patterns of evolutionary development. The studies presented here describe the evolution of the components of SMIM45. One ORF consists of an ultra-conserved 68 amino acid (aa) sequence, whose origins can be traced beyond the evolutionary age of divergence of the elephant shark, ~462 MYA. The silencer also has ancient origins, but it has a complex and divergent pattern of evolutionary formation, as it overlaps both at the 68 aa ORF and the intervening sequence. The other ORF consists of 107 aa. It develops during primate evolution but is found to originate de novo from an ancestral non-coding genomic region with root origins within the Afrothere clade of placental mammals, whose evolutionary age of divergence is ~99 MYA. The formation of the complete 107 aa ORF during primate evolution is outlined, whereby sequence development is found to occur through biased mutations, with disruptive random mutations that also occur but lead to a dead-end. The 107 aa ORF is of particular significance, as there is evidence to suggest it is a protein that may function in human brain development. Its evolutionary formation presents a view of a human-specific ORF and its linked silencer that were predetermined in non-primate ancestral species. The genomic position of the silencer offers interesting possibilities for the regulation of transcription of the 107 aa ORF. A hypothesis is presented with respect to possible spatiotemporal expression of the 107 aa ORF in embryonic tissues.

Why old duplicated genes are not thrown away in (paleo)polyploids? Example from the petC gene in Brassica napus.

Duplication of genes at different time period, through recurrent and frequent polyploidization events, have played a major role in plant evolution, adaptation and diversification. Interestingly, some of the ancestral duplicated genes (referred as paleologs), have been maintained for millions of years, and there is still a poor knowledge of the reasons of their retention, especially when testing the phenotypic effect of individual copies by using functional genetic approaches. To fill this gap, we performed functional genetic (CRISPR-Cas9), physiological, transcriptomic and evolutionary studies to finely investigate this open question, taking the example of the petC gene (involved in cytochrome b6/f and thus impacting photosynthesis) that is present in four paleologous copies in the oilseed crop Brassica napus. RNA-Seq and selective pressure analyses suggested that all paleologous copies conserved the same function and that they were all highly transcribed. Thereafter, the Knock Out (K.O.) of one, several or all petC copies highlighted that all paleologous copies have to be K.O. to suppress the gene function. In addition, we could determine that phenotypic effects in single and double mutants could only be deciphered in high light conditions. Interestingly, we did not detect any significant differences between single mutants K.O. for either the A03 or A09 copy (despite being differentially transcribed), or even between mutants for a single or two petC copies. Altogether, this work revealed that petC paleologs have retained their ancestral function and that the retention of these copies is explained by their compensatory role, especially in optimal environmental conditions.

Sirtuin Evolution at the Dawn of Animal Life.

Sirtuins are a family of proteins that protect against cellular injury and aging; understanding their evolution should reveal fundamental mechanisms governing longevity. "Early-branching" animals such as sea sponges and jellyfish have been understudied in previous analyses of sirtuin diversity. These organisms not only hold important positions at the base of the evolutionary tree, but also have unique aging dynamics that defy convention, such as quasi-immortality and high regenerative capacity. In this study, we survey the evolution of sirtuin proteins in animals, with a focus on the oldest living lineages. We describe previously unrecognized expansions of "Class IV" and "Class I" sirtuins around the origin of animals, raising the number of sirtuin families in the last common ancestor to at least nine. Most of these undescribed sirtuins have been lost in vertebrates and other bilaterian animals. Our work also clarifies the evolution of PNC1 and NAMPT enzymes that carry out the rate-limiting step in sirtuin-related NAD+ biosynthesis. The genes for PNC1 and NAMPT enzymes were both present in the first animals, with the genes being lost a minimum of 11 and 13 times, respectively, over the course of animal evolution. We propose that species with these ancestral gene repertoires are ideal model organisms for studying the genetic regulation of animal longevity and will provide clues to increasing longevity in humans.

Phylogenomics insights into gene evolution, rapid species diversification, and morphological innovation of the apple tribe (Maleae, Rosaceae).

Maleae is one of the most widespread tribes of Rosaceae and includes several important fruit crops and ornamental plants. We used nuclear genes from 62 transcriptomes/genomes, including 26 newly generated transcriptomes, to reconstruct a well-supported phylogeny and study the evolution of fruit and leaf morphology and the possible effect of whole genome duplication (WGD). Our phylogeny recovered 11 well-supported clades and supported the monophyly of most genera (except Malus, Sorbus, and Pourthiaea) with at least two sampled species. A WGD was located to the most recent common ancestor (MRCA) of Maleae and dated to c. 54 million years ago (Ma) near the Early Eocene Climatic Optimum, supporting Gillenieae (x = 9) being a parental lineage of Maleae (x = 17) and including duplicate regulatory genes related to the origin of the fleshy pome fruit. Whole genome duplication-derived paralogs that are retained in specific lineages but lost in others are predicted to function in development, metabolism, and other processes. An upshift of diversification and innovations of fruit and leaf morphologies occurred at the MRCA of the Malinae subtribe, coinciding with the Eocene-Oligocene transition (c. 34 Ma), following a lag from the time of the WGD event. Our results provide new insights into the Maleae phylogeny, its rapid diversification, and morphological and molecular evolution.

Identification and characterization of sperm motility-initiating substance-2 gene in internally fertilizing Cynops species.

Sperm motility-initiating substance (SMIS) is an oviductal protein critical for internal fertilization in urodeles. It contributes to the establishment of various reproductive modes in amphibians and is thus a unique research model for the gene evolution of gamete-recognizing ligands that have diversified among animal species. In this study, a paralogous SMIS gene, smis2, was identified via the RNA sequencing of the oviduct of the newt, Cynops pyrrhogaster. The base sequence of the smis2 gene was homologous (˃90%) to that of the original smis gene (smis1), and deduced amino acid sequences of both genes conserved six cysteine residues essential for the cysteine knot motif. Furthermore, smis2 complementary DNA was identified in the oviduct of Cynops ensicauda, and the base substitution patterns also suggested that the smis gene was duplicated in the Salamandridae. Nonsynonymous/synonymous substitution ratios of smis1 and smis2 genes were 0.79 and 2.6, respectively, suggesting that smis2 gene evolution was independently driven by positive selection. Amino acid substitutions were concentrated in the cysteine knot motif of SMIS2. The smis2 gene was expressed in some organs in addition to the oviduct; in contrast, SMIS1 was only expressed in the oviduct. The SMIS2 protein was suggested to be produced and secreted at least in the oviduct and redundantly act in sperm. These results suggest that smis1 plays the original role in the oviduct, whereas smis2 may undergo neofunctionalization, which rarely occurs in gene evolution.

Gene Characterization of Nocturnin Paralogues in Goldfish: Full Coding Sequences, Structure, Phylogeny and Tissue Expression.

The aim of this work is the full characterization of all the nocturnin (noc) paralogues expressed in a teleost, the goldfish. An in silico analysis of the evolutive origin of noc in Osteichthyes is performed, including the splicing variants and new paralogues appearing after teleostean 3R genomic duplication and the cyprinine 4Rc. After sequencing the full-length mRNA of goldfish, we obtained two isoforms for noc-a (noc-aa and noc-ab) with two splice variants (I and II), and only one for noc-b (noc-bb) with two transcripts (II and III). Using the splicing variant II, the prediction of the secondary and tertiary structures renders a well-conserved 3D distribution of four α-helices and nine ß-sheets in the three noc isoforms. A synteny analysis based on the localization of noc genes in the patrilineal or matrilineal subgenomes and a phylogenetic tree of protein sequences were accomplished to stablish a classification and a long-lasting nomenclature of noc in goldfish, and valid to be extrapolated to allotetraploid Cyprininae. Finally, both goldfish and zebrafish showed a broad tissue expression of all the noc paralogues. Moreover, the enriched expression of specific paralogues in some tissues argues in favour of neo- or subfunctionalization.

[Structure and Evolution of the AqE Gene in Insects].

The AqE gene encodes a sulfolactate dehydrogenase-like enzyme of the LDH2/MDG2 oxidoreductase family. The gene is found in bacteria and fungi, as well as in animals and plants whose lifestyles are associated with aquatic environments. The AqE gene is present in arthropods and, in particular, insects that are predominantly terrestrial. The distribution and structure of AqE was studied in insects in order to trace its evolutionary fate. The AqE gene was found to be absent from certain insect orders and suborders, being apparently lost. AqE duplication or multiplication was observed in some orders. AqE was found to vary both in length and intron-exon structure, from intronless to multi-intron. An ancient nature was demonstrated for AqE multiplication in insects, while younger duplications were also detected. It was assumed that a new function might be acquired by the gene with the formation of paralogs.

Mitogenomes Reveal Alternative Initiation Codons and Lineage-Specific Gene Order Conservation in Echinoderms.

The mitochondrial genetic code is much more varied than the standard genetic code. The invertebrate mitochondrial code, for instance, comprises six initiation codons, including five alternative start codons. However, only two initiation codons are known in the echinoderm and flatworm mitochondrial code, the canonical ATG and alternative GTG. Here, we analyzed 23 Asteroidea mitogenomes, including ten newly sequenced species and unambiguously identified at least two other start codons, ATT and ATC, both of which also initiate translation of mitochondrial genes in other invertebrates. These findings underscore the diversity of the genetic code and expand upon the suite of initiation codons among echinoderms to avoid erroneous annotations. Our analyses have also uncovered the remarkable conservation of gene order among asteroids, echinoids, and holothuroids, with only an interchange between two gene positions in asteroids over ∼500 Ma of echinoderm evolution.

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.