A comparative analysis of stably expressed genes across diverse angiosperms exposes flexibility in underlying promoter architecture

Promoters regulate both the amplitude and pattern of gene expression-key factors needed for optimization of many synthetic biology applications. Previous work in Arabidopsis found that promoters that contain a TATA-box element tend to be expressed only under specific conditions or in particular tissues, while promoters that lack any known promoter elements, thus designated as Coreless, tend to be expressed more uniformly. To test whether this trend represents a conserved promoter design rule, we identified stably expressed genes across multiple angiosperm species using publicly available RNA-seq data. Comparisons between core promoter architectures and gene expression stability revealed differences in core promoter usage in monocots and eudicots. Furthermore, when tracing the evolution of a given promoter across species, we found that core promoter type was not a strong predictor of expression pattern. Our analysis suggests that core promoter types are correlative rather than causative in promoter expression patterns and highlights the challenges in finding or building constitutive promoters that will work across diverse plant species.

Dicarboxylic Amino Acid Permease 7219 Regulates Fruiting Body Type of Auricularia heimuer

Auricularia heimuer is a widely cultivated jelly mushroom. The fruiting bodies are categorized into cluster and chrysanthemum types. With changing consumer demands and the need to reduce bio-waste, the demand for clustered fruiting bodies is increasing. Therefore, gene mining for fruiting body types is a matter of urgency. We determined that the A. heimuer locus for fruiting body type was located at one end of the genetic linkage map. The locus was localized between the markers D23860 and D389 by increasing the density of the genetic linkage map. BlastN alignment showed that the marker SCL-18 was also located between D23860 and D389, and a total of 25 coding genes were annotated within this interval. Through parental transcriptome analysis and qRT-PCR verification, the locus g7219 was identified as the gene controlling the fruiting body type. A single-nucleotide substitution in the TATA box of g7219 was detected between the parents. By PCR amplification of the promoter region of g7219, the TATA-box sequences of the cluster- and chrysanthemum-type strains were found to be CATAAAA and TATAAAA, respectively. This study provides a foundation for the breeding of fruiting body types and strain improvement of A. heimuer.

RNA Pol III promoters-key players in precisely targeted plant genome editing

The clustered regularly interspaced short palindrome repeat (CRISPR)/CRISPR-associated protein Cas) system is a powerful and highly precise gene-editing tool in basic and applied research for crop improvement programs. CRISPR/Cas tool is being extensively used in plants to improve crop yield, quality, and nutritional value and make them tolerant to environmental stresses. CRISPR/Cas system consists of a Cas protein with DNA endonuclease activity and one CRISPR RNA transcript that is processed to form one or several short guide RNAs that direct Cas9 to the target DNA sequence. The expression levels of Cas proteins and gRNAs significantly influence the editing efficiency of CRISPR/Cas-mediated genome editing. This review focuses on insights into RNA Pol III promoters and their types that govern the expression levels of sgRNA in the CRISPR/Cas system. We discussed Pol III promoters structural and functional characteristics and their comparison with Pol II promoters. Further, the use of synthetic promoters to increase the targeting efficiency and overcome the structural, functional, and expressional limitations of RNA Pol III promoters has been discussed. Our review reports various studies that illustrate the use of endogenous U6/U3 promoters for improving editing efficiency in plants and the applicative approach of species-specific RNA pol III promoters for genome editing in model crops like Arabidopsis and tobacco, cereals, legumes, oilseed, and horticultural crops. We further highlight the significance of optimizing these species-specific promoters' systematic identification and validation for crop improvement and biotic and abiotic stress tolerance through CRISPR/Cas mediated genome editing.

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.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-022-08829-6.

Genome-Wide Prediction of Transcription Start Sites in Conifers

The identification of promoters is an essential step in the genome annotation process, providing a framework for gene regulatory networks and their role in transcription regulation. Despite considerable advances in the high-throughput determination of transcription start sites (TSSs) and transcription factor binding sites (TFBSs), experimental methods are still time-consuming and expensive. Instead, several computational approaches have been developed to provide fast and reliable means for predicting the location of TSSs and regulatory motifs on a genome-wide scale. Numerous studies have been carried out on the regulatory elements of mammalian genomes, but plant promoters, especially in gymnosperms, have been left out of the limelight and, therefore, have been poorly investigated. The aim of this study was to enhance and expand the existing genome annotations using computational approaches for genome-wide prediction of TSSs in the four conifer species: loblolly pine, white spruce, Norway spruce, and Siberian larch. Our pipeline will be useful for TSS predictions in other genomes, especially for draft assemblies, where reliable TSS predictions are not usually available. We also explored some of the features of the nucleotide composition of the predicted promoters and compared the GC properties of conifer genes with model monocot and dicot plants. Here, we demonstrate that even incomplete genome assemblies and partial annotations can be a reliable starting point for TSS annotation. The results of the TSS prediction in four conifer species have been deposited in the Persephone genome browser, which allows smooth visualization and is optimized for large data sets. This work provides the initial basis for future experimental validation and the study of the regulatory regions to understand gene regulation in gymnosperms.

Dysfunction of CYC2g is responsible for the evolutionary shift from radiate to disciform flowerheads in the Chrysanthemum group (Asteraceae: Anthemideae)

Evolutionary shifts among radiate, disciform and discoid flowerheads have occurred repeatedly in a number of major lineages across the Asteraceae phylogeny; such transitions may also appear within evolutionarily young groups. Although several studies have demonstrated that CYC2 genes partake in regulating floral morphogenesis in Asteraceae, the evolution of capitulum forms within a recently diverging lineage has remained poorly understood. Here, we study the molecular regulation of the shift from a radiate to a disciform capitulum within the Chrysanthemum group. This is a recently radiating group mainly comprising two genera, Chrysanthemum and Ajania, that are phylogenetically intermingled but distinct in flowerhead morphology: Chrysanthemum spp. with radiate capitula and Ajania spp. with disciform capitula. We found that the morphogenesis of zygomorphy in the marginal floret in Ajania was disrupted soon after floral primordium emergence; CYC2g, one of the CYC2 copies that was expressed prominently in the ray floret of Chrysanthemum was not expressed in flowerheads of Ajania. Weakening the expression of ClCYC2g in Chrysanthemum lavandulifolium led to the gradual transition of a ray flower toward the disc-like form. Molecular evolutionary analyses indicated that the disciform capitulum might have evolved only once, approximately 8 Mya, arising from dysfunction of the CYC2g orthologs. A 20-nt deletion, including a putative TATA-box of the Ajania-type CYC2g promoter, appeared to inhibit the expression of the gene. Considering the divergent habitats of Chrysanthemum and Ajania, we propose that the shift from radiate to disciform capitulum must have been related to changes in pollination strategies under selective pressure.

The highly diverse TATA box-binding proteins among protists: A review

Transcription is the first step of gene expression regulation and is a fundamental mechanism for establishing the viability and development of a cell. The TATA box-binding protein (TBP) interaction with a TATA box in a promoter is one of the best studied mechanisms in transcription initiation. TBP is a transcription factor that is highly conserved from archaea to humans and is essential for the transcription initiated by each of the three RNA polymerases. In addition, the discovery of TBP-related factor 1 (TRF1) and other factors related to TBP shed light on the variability among transcription initiation complexes, thus demonstrating that the compositions of these complexes are, in fact, more complicated than originally believed. Despite these facts, the majority of studies on transcription have been performed on animal, plant and fungal cells, which serve as canonical models, and information regarding protist cells is relatively scarce. The aim of this work is to review the diversity of the TBPs that have been documented in protists and describe some of the specific features that differentiate them from their counterparts in higher eukaryotes.

Natural Selection Equally Supports the Human Tendencies in Subordination and Domination: A Genome-Wide Study With in silico Confirmation and in vivo Validation in Mice

We proposed the following heuristic decision-making rule: "IF {an excess of a protein relating to the nervous system is an experimentally known physiological marker of low pain sensitivity, fast postinjury recovery, or aggressive, risk/novelty-seeking, anesthetic-like, or similar agonistic-intolerant behavior} AND IF {a single nucleotide polymorphism (SNP) causes overexpression of the gene encoding this protein} THEN {this SNP can be a SNP marker of the tendency in dominance} WHILE {underexpression corresponds to subordination} AND vice versa." Using this decision-making rule, we analyzed 231 human genes of neuropeptidergic, non-neuropeptidergic, and neurotrophinergic systems that encode neurotrophic and growth factors, interleukins, neurotransmitters, receptors, transporters, and enzymes. These proteins are known as key factors of human social behavior. We analyzed all the 5,052 SNPs within the 70 bp promoter region upstream of the position where the protein-coding transcript starts, which were retrieved from databases Ensembl and dbSNP using our previously created public Web service SNP_TATA_Comparator (http://beehive.bionet.nsc.ru/cgi-bin/mgs/tatascan/start.pl). This definition of the promoter region includes all TATA-binding protein (TBP)-binding sites. A total of 556 and 552 candidate SNP markers contributing to the dominance and the subordination, respectively, were uncovered. On this basis, we determined that 231 human genes under study are subject to natural selection against underexpression (significance p < 0.0005), which equally supports the human tendencies in domination and subordination such as the norm of a reaction (plasticity) of the human social hierarchy. These findings explain vertical transmission of domination and subordination traits previously observed in rodent models. Thus, the results of this study equally support both sides of the century-old unsettled scientific debate on whether both aggressiveness and the social hierarchy among humans are inherited (as suggested by Freud and Lorenz) or are due to non-genetic social education, when the children are influenced by older individuals across generations (as proposed by Berkowitz and Fromm).

Modeling the Function of TATA Box Binding Protein in Transcriptional Changes Induced by HIV-1 Tat in Innate Immune Cells and the Effect of Methamphetamine Exposure

Innate immune cells are targets of HIV-1 infection in the Central Nervous System (CNS), generating neurological deficits. Infected individuals with substance use disorders as co-morbidities, are more likely to have aggravated neurological disorders, higher CNS viral load and inflammation. Methamphetamine (Meth) is an addictive stimulant drug, commonly among HIV+ individuals. The molecular basis of HIV direct effects and its interactions with Meth in host response, at the gene promoter level, are not well understood. The main HIV-1 peptide acting on transcription is the transactivator of transcription (Tat), which promotes replication by recruiting a Tata-box binding protein (TBP) to the virus long-terminal repeat (LTR). We tested the hypothesis that Tat can stimulate host gene expression through its ability to increase TBP, and thus promoting its binding to promoters that bear Tata-box binding motifs. Genes with Tata-box domains are mainly inducible, early response, and involved in inflammation, regulation and metabolism, relevant in HIV pathogenesis. We also tested whether Tat and Meth interact to trigger the expression of Tata-box bearing genes. The THP1 macrophage cell line is a well characterized innate immune cell system for studying signal transduction in inflammation. These cells are responsive to Tat, as well as to Meth, by recruiting RNA Polymerase (RNA Pol) to inflammatory gene promoters, within 15 min of stimulation (1). THP-1 cells, including their genetically engineered derivatives, represent valuable tools for investigating monocyte structure and function in both health and disease, as a consistent system (2). When differentiated, they mimic several aspects of the response of macrophages, and innate immune cells that are the main HIV-1 targets within the Central Nervous System (CNS). THP1 cells have been used to characterize the impact of Meth and resulting neurotransmitters on HIV entry (1), mimicking the CNS micro-environment. Integrative consensus sequence analysis in genes with enriched RNA Pol, revealed that TBP was a major transcription factor in Tat stimulation, while the co-incubation with Meth shifted usage to a distinct and diversified pattern. For validating these findings, we engineered a THP1 clone to be deficient in the expression of all major TBP splice variants, and tested its response to Tat stimulation, in the presence or absence of Meth. Transcriptional patterns in TBP-sufficient and deficient clones confirmed TBP as a dominant transcription factor in Tat stimulation, capable of inducing genes with no constitutive expression. However, in the presence of Meth, TBP was no longer necessary to activate the same genes, suggesting promoter plasticity. These findings demonstrate TBP as mechanism of host-response activation by HIV-1 Tat, and suggest that promoter plasticity is a challenge imposed by co-morbid factors such as stimulant drug addiction. This may be one mechanism responsible for limited efficacy of therapeutic approaches in HIV+ Meth abusers.

Core cis-element variation confers subgenome-biased expression of a transcription factor that functions in cotton fiber elongation

Cotton cultivars have evolved to produce extensive, long, seed-born fibers important for the textile industry, but we know little about the molecular mechanism underlying spinnable fiber formation. Here, we report how PACLOBUTRAZOL RESISTANCE 1 (PRE1) in cotton, which encodes a basic helix-loop-helix (bHLH) transcription factor, is a target gene of spinnable fiber evolution. Differential expression of homoeologous genes in polyploids is thought to be important to plant adaptation and novel phenotypes. PRE1 expression is specific to cotton fiber cells, upregulated during their rapid elongation stage and A-homoeologous biased in allotetraploid cultivars. Transgenic studies demonstrated that PRE1 is a positive regulator of fiber elongation. We determined that the natural variation of the canonical TATA-box, a regulatory element commonly found in many eukaryotic core promoters, is necessary for subgenome-biased PRE1 expression, representing a mechanism underlying the selection of homoeologous genes. Thus, variations in the promoter of the cell elongation regulator gene PRE1 have contributed to spinnable fiber formation in cotton. Overexpression of GhPRE1 in transgenic cotton yields longer fibers with improved quality parameters, indicating that this bHLH gene is useful for improving cotton fiber quality.

Distinct core promoter codes drive transcription initiation at key developmental transitions in a marine chordate

Background

Development is largely driven by transitions between transcriptional programs. The initiation of transcription at appropriate sites in the genome is a key component of this and yet few rules governing selection are known. Here, we used cap analysis of gene expression (CAGE) to generate bp-resolution maps of transcription start sites (TSSs) across the genome of Oikopleura dioica, a member of the closest living relatives to vertebrates.

Results

Our TSS maps revealed promoter features in common with vertebrates, as well as striking differences, and uncovered key roles for core promoter elements in the regulation of development. During spermatogenesis there is a genome-wide shift in mode of transcription initiation characterized by a novel core promoter element. This element was associated with > 70% of male-specific transcription, including the use of cryptic internal promoters within operons. In many cases this led to the exclusion of trans-splice sites, revealing a novel mechanism for regulating which mRNAs receive the spliced leader. Binding of the cell cycle regulator, E2F1, is enriched at the TSS of maternal genes in endocycling nurse nuclei. In addition, maternal promoters lack the TATA-like element found in zebrafish and have broad, rather than sharp, architectures with ordered nucleosomes. Promoters of ribosomal protein genes lack the highly conserved TCT initiator. We also report an association between DNA methylation on transcribed gene bodies and the TATA-box.

Conclusions

Our results reveal that distinct functional promoter classes and overlapping promoter codes are present in protochordates like in vertebrates, but show extraordinary lineage-specific innovations. Furthermore, we uncover a genome-wide, developmental stage-specific shift in the mode of TSS selection. Our results provide a rich resource for the study of promoter structure and evolution in Metazoa.

Electronic supplementary material

The online version of this article (10.1186/s12864-018-4504-5) contains supplementary material, which is available to authorized users.

Effect of temperature on the structure and hydration layer of TATA-box DNA: A molecular dynamics simulation study

DNA within the living cells experiences a diverse range of temperature, ranging from freezing condition to hot spring water. How the structure, the mechanical properties of DNA, and the solvation dynamics around DNA changes with the temperature is important to understand the functionality of DNA under those acute temperature conditions. In that notion, we have carried out molecular dynamics simulations of a DNA oligomer, containing TATA-box sequence for three different temperatures (250K, 300K and 350K). We observed that the structure of the DNA, in terms of backbone torsion angles, sugar pucker, base pair parameters, and base pair step parameters, did not show any unusual properties within the studied range of temperatures, but significant structural alteration was noticed between BI and BII forms at higher temperature. As expected, the flexibility of the DNA, in terms of the torsional rigidity and the bending rigidity is highly temperature dependent, confirming that flexibility increases with increase in temperature. Additionally, the groove widths of the studied DNA showed temperature sensitivity, specifically, the major groove width decreases and the minor groove width increases, respectively, with the increase in temperature. We observed that at higher temperature, water around both the major and the minor groove of the DNA is less structured. However, the water dynamics around the minor groove of the DNA is more restricted as compared to the water around the major groove throughout the studied range of temperatures, without any anomalous behavior.

A single nucleotide change in a core promoter is involved in the progressive overexpression of the duplicated CYP9M10 haplotype lineage in Culex quinquefasciatus

Although the importance of cis-acting mutations on detoxification enzyme genes for insecticide resistance is widely accepted, only a few of them have been determined as concrete mutations present in genomic DNA till date. The overexpression of a cytochrome P450 gene, CYP9M10, is associated with pyrethroid resistance in the southern house mosquito Culex quinquefasciatus. The haplotypes of CYP9M10 exhibiting overexpression (resistant haplotypes) belong to one specific phylogenetic lineage that shares high nucleotide sequence homology and the same insertion of a transposable element. Among the resistant haplotypes, allelic progression involving an additional cis-acting mutation and gene duplication evolved a CYP9M10 haplotype associated with extremely high transcription and strong pyrethroid resistance. Here we show that a single nucleotide substitution G-27A, which is located near the transcription start site of CYP9M10, is involved in the progression of the duplicated haplotype lineage. The deletion of a 7-bp AT-rich sequence that includes nucleotide -27 inhibited the initiation of transcription from the original transcriptional initiation site. The mutation was suspected to reside within a core promoter, TATA-box, of CYP9M10.

Molecular organization of the 5S rDNA gene type II in elasmobranchs

The 5S rDNA gene is a non-coding RNA that can be found in 2 copies (type I and type II) in bony and cartilaginous fish. Previous studies have pointed out that type II gene is a paralog derived from type I. We analyzed the molecular organization of 5S rDNA type II in elasmobranchs. Although the structure of the 5S rDNA is supposed to be highly conserved, our results show that the secondary structure in this group possesses some variability and is different than the consensus secondary structure. One of these differences in Selachii is an internal loop at nucleotides 7 and 112. These mutations observed in the transcribed region suggest an independent origin of the gene among Batoids and Selachii. All promoters were highly conserved with the exception of BoxA, possibly due to its affinity to polymerase III. This latter enzyme recognizes a dT4 sequence as stop signal, however in Rajiformes this signal was doubled in length to dT8. This could be an adaptation toward a higher efficiency in the termination process. Our results suggest that there is no TATA box in elasmobranchs in the NTS region. We also provide some evidence suggesting that the complexity of the microsatellites present in the NTS region play an important role in the 5S rRNA gene since it is significantly correlated with the length of the NTS.

Distinct role of core promoter architecture in regulation of light-mediated responses in plant genes

In the present study, we selected four distinct classes of light-regulated promoters. The light-regulated promoters can be distinctly grouped into either TATA-box-containing or TATA-less (initiator-containing) promoters. Further, using either native promoters or their swapped versions of core promoter elements, we established that TATA-box and Inr (Initiator) elements have distinct mechanisms which are involved in light-mediated regulation, and these elements are not swappable. We identified that mutations in either functional TATA-box or Inr elements lead to the formation of nucleosomal structure. The nucleotide diversity in either the TATA-box or Inr element in Arabidopsis ecotypes proposes that the nucleotide variation in core promoters can alter the gene expression. We show that motif overrepresentation in light-activated promoters encompasses different specific regulatory motifs present downstream of TSS (transcription start site), and this might serve as a key factor in regulating light promoters which are parallel with these elements. Finally, we conclude that the TATA-box or Inr element does not act in isolation, but our results clearly suggests the probable involvement of other distinct core promoter elements in concurrence with the TATA-box or Inr element to impart selectivity to light-mediated transcription.

Investigation of dmyc Promoter and Regulatory Regions

Products of the myc gene family integrate extracellular signals by modulating a wide range of their targets involved in cellular biogenesis and metabolism; the purpose of this integration is to regulate cell death, proliferation, and differentiation. However, understanding the regulation of myc at the transcription level remains a challenge. We performed rapid amplification of dmyc cDNA ends (5' RACE) and mapped the transcription start site at P1 promoter, 18 base pairs upstream of the start of the known EST GM01143 and within the 5' UTR. Our data show that the first TATA box, previously computationally predicted, is utilized to generate dmyc full length mRNA. The largest transcript contains all three exons, generated after the removal of the introns by constitutively regulated splicing events. Further investigation of Downstream Promoter Element (DPE) was achieved by studying lacZ reporter activity; investigation revealed that this element and its upstream cluster of binding sites are required for the dmyc intron 2 activity. These findings may provide valuable tools for further analysis of dmyc cis-elements.

Promoter Structure of the RNA Polymerase II Large Subunit Gene in Caenorhabditis elegans and C. briggsae

The 5'-end of the Caenorhabditis elegans ama-1 gene transcript, which encodes the largest subunit of RNA polymerase II, was cloned. Sequencing revealed that the message is trans-spliced. To characterize the Ce-ama-1 promoter, DNA sequence spanning 3 kb upstream from the initiation codon was determined. Typical elements, such as TATA and Spl sites, were absent. The homologue of ama-1 in C. briggsae, Cb-ama-1, was isolated and its 5' flanking sequence compared with that of Ce-ama-1, revealing only limited similarity, although both sequences included a potential initiator-class transcriptional regulator and phased repeats of an ATC motif. The latter elements are postulated to facilitate DNA bending and may play a role in transcription regulation.