Inhibition of DNA polymerase reaction by pyrimidine nucleotide analogues lacking the 2-keto group

To investigate the influence of the pyrimidine 2-keto group on selection of nucleotides for incorporation into DNA by polymerases, we have prepared two C nucleoside triphosphates that are analogues of dCTP and dTTP, namely 2-amino-5-(2'-deoxy-beta-d-ribofuranosyl)pyridine-5'-triphosphate (d*CTP) and 5-(2'-deoxy- beta-d-ribofuranosyl)-3-methyl-2-pyridone-5'-triphosphate (d*TTP) respectively. Both proved strongly inhibitory to PCR catalysed by Taq polymerase; d*TTP rather more so than d*CTP. In primer extension experiments conducted with either Taq polymerase or the Klenow fragment of Escherichia coli DNA polymerase I, both nucleotides failed to substitute for their natural pyrimidine counterparts. Neither derivative was incorporated as a chain terminator. Their capacity to inhibit DNA polymerase activity may well result from incompatibility with the correctly folded form of the polymerase enzyme needed to stabilize the transition state and catalyse phosphodiester bond formation.

tRNA-derived fragments target the small ribosomal subunit to fine-tune translation

Post-transcriptional cleavage of RNA molecules to generate smaller fragments is a widespread mechanism that enlarges the structural and functional complexity of cellular RNomes. In particular, fragments deriving from both precursor and mature tRNAs represent one of the rapidly growing classes of post-transcriptional RNA pieces. Importantly, these tRNA-derived fragments (tRFs) possess distinct expression patterns, abundance, cellular localizations, or biological roles compared with their parental tRNA molecules [1]. Here we present evidence that tRFs from the archaeon Haloferax volcanii directly bind to ribosomes. In a previous genomic screen for ribosome-associated small RNAs we have identified a 26 residue long fragment originating from the 5’ part of valine tRNA (Val-tRF) to be by far the most abundant tRF in H. volcanii [2]. The Val-tRF is processed in a stress- dependent manner and was found to primarily target the small ribosomal subunit in vitro and in vivo. Translational activity was markedly reduced in the presence of Val-tRF, while control RNA fragments of similar length did not show inhibition of protein biosynthesis. Crosslinking experiments and subsequent primer extension analyses revealed the Val-tRF interaction site to surround the mRNA path in the 30S subunit. In support of this, binding experiments demonstrated that Val-tRF does compete with mRNAs for ribosome binding. Therefore this tRF represents a ribosome-bound non-protein-coding RNA (ncRNA) capable of regulating gene expression in H. volcanii under environmental stress conditions probably by fine-tuning the rate of protein production [1].

tRNA-derived fragments target the small ribosomal subunit to fine-tune translation

Post-transcriptional cleavage of RNA molecules to generate smaller fragments is a widespread mechanism that enlarges the structural and functional complexity of cellular RNomes. In particular, fragments deriving from both precursor and mature tRNAs represent one of the rapidly growing classes of post-transcriptional RNA pieces. Importantly, these tRNA-derived fragments (tRFs) possess distinct expression patterns, abundance, cellular localizations, or biological roles compared with their parental tRNA molecules [1]. Here we present evidence that tRFs from the archaeon Haloferax volcanii directly bind to ribosomes. In a previous genomic screen for ribosome-associated small RNAs we have identified a 26 residue long fragment originating from the 5’ part of valine tRNA (Val-tRF) to be by far the most abundant tRF in H. volcanii [2]. The Val-tRF is processed in a stress- dependent manner and was found to primarily target the small ribosomal subunit in vitro and in vivo. Translational activity was markedly reduced in the presence of Val-tRF, while control RNA fragments of similar length did not show inhibition of protein biosynthesis. Crosslinking experiments and subsequent primer extension analyses revealed the Val-tRF interaction site to surround the mRNA path in the 30S subunit. In support of this, binding experiments demonstrated that Val-tRF does compete with mRNAs for ribosome binding. Therefore this tRF represents a ribosome-associated non-protein-coding RNA (rancRNA) capable of regulating gene expression in H. volcanii under environmental stress conditions probably by fine-tuning the rate of protein production [3].

tRNA-derived fragments target the small ribosomal subunit to fine-tune translation

Post-transcriptional cleavage of RNA molecules to generate smaller fragments is a widespread mechanism that enlarges the structural and functional complexity of cellular RNomes. In particular, fragments deriving from both precursor and mature tRNAs represent one of the rapidly growing classes of post-transcriptional RNA pieces. Importantly, these tRNA-derived fragments (tRFs) possess distinct expression patterns, abundance, cellular localizations, or biological roles compared with their parental tRNA molecules [1]. Here we present evidence that tRFs from the halophilic archaeon Haloferax volcanii directly bind to ribosomes. In a previous genomic screen for ribosome-associated small RNAs we have identified a 26 residue long fragment originating from the 5’ part of valine tRNA (Val-tRF) to be by far the most abundant tRF in H. volcanii [2]. The Val-tRF is processed in a stress-dependent manner and was found to primarily target the small ribosomal subunit in vitro and in vivo. Translational activity was markedly reduced in the presence of Val-tRF, while control RNA fragments of similar length did not show inhibition of protein biosynthesis. Crosslinking experiments and subsequent primer extension analyses revealed the Val-tRF interaction site to surround the mRNA path in the 30S subunit. In support of this, binding experiments demonstrated that Val-tRF does compete with mRNAs for ribosome binding. Therefore this tRF represents a ribosome-associated non-coding RNA (rancRNA) capable of regulating gene expression in H. volcanii under environmental stress conditions probably by fine-tuning the rate of protein production [3].

STRUCTURAL ANALYSIS OF THE TRANSCRIPTION PRODUCTS OF A TEMPERATURE-SENSITIVE PHENOTYPIC MUTANT OF MOLONEY SARCOMA VIRUS--MUSV TS110

Cells infected with a temperature sensitive phenotypic mutant of Moloney sarcoma virus (MuSVts110) exhibit a transformed phenotype at 33('(DEGREES)) and synthesize two virus specific proteins, p85('gag-mos), a gag-mos fusion protein and p58('gag), a truncated gag precursor protein (the gag gene codes for viral structural proteins and mos is the MuSV transforming gene). At 39('(DEGREES)) only p58('gag) is synthesized and the morphology of the cells is similar to uninfected NRK parental cells. Two MuSVts110 specific RNAs are made in MuSVts110-infected cells, one of 4.0 kb in length, the other of 3.5 kb. Previous work indicated that each of these RNAs arose by a single central deletion of parental MuSV genetic material, and that p58('gag) was made by the 4.0 kb RNA and p85('gag-mos) from the 3.5 kb RNA. The objective of my dissertation research was to map precisely the deletion boundaries of both of the MuSVts110 RNAs, and to determine the proper reading frame across both deletion borders. This work succeeded in arriving at the following conclusions: (a) Using S-1 nuclease analysis and primer extension sequencing, it was found that the 4.0 kb MuSVts110 RNA arose by a 1488 base deletion of 5.2 kb parental MuSV genomic RNA. This deletion resulted in an out of frame fusion of the gag and mos genes that resulted in the formation of a "stop" codon which causes termination of translation just beyond the c-terminus of the gag region. Thus, this RNA can only be translated into the truncated gag protein p58('gag). (b) S-1 analysis of RNA from cells cultivated at different temperatures demonstrated that the 4.0 kb RNA was synthesized at all temperatures but that synthesis of the 3.5 kb RNA was temperature sensitive. These observations supported the data derived from blot hybridization experiments the interpretation of which argued for the existence of a single provirus in MuSVts110 infected cells, and hence only a single primary transcript (the 4.0 kb RNA). (c) Analyses similar to those described in (a) above showed that the 3.5 kb RNA was derived from the 4.0 kb MuSVts110 RNA by a further deletion of 431 bases, fusing the gag and mos genes into a continuous reading frame capable of directing synthesis of the p85('gag-mos) protein. These sequence data and the presence of only one MuSVts110-specific provirus, indicate that a splice mechanism is employed to generate the 3.5 kb RNA since the gag and mos genes are observed to be fused in frame in this RNA. . . . (Author's abstract exceeds stipulated maximum length. Discontinued here with permission of author.) UMI ^

Cloning, characterization and regulation of the soluble guanylyl cyclase alpha1 and beta1 genes

Cell signaling by nitric oxide (NO) through soluble guanylyl cyclase (sGC) and cGMP production regulates physiological responses such as smooth muscle relaxation, neurotransmission, and cell growth and differentiation. Although the NO receptor, sGC, has been studied extensively at the protein level, information on regulation of the sGC genes remains elusive. In order to understand the molecular mechanisms involved at the level of gene expression, cDNA and genomic fragments of the murine sGCα1 subunit gene were obtained through library screenings. Using the acquired clones, the sGCα 1 gene structure was determined following primer extension, 3 ′RACE and intron/exon boundary analyses. The basal activity of several 5′-flanking regions (putative promoter regions) for both the α1 and β1 sGC subunits were determined following their transfection into mouse N1E-115 neuroblastoma and rat RENE1Δ14 uterine epithelial cells using a luciferase reporter plasmid. Using the sGC sequences, real-time RT-PCR assays were designed to measure mRNA levels of the sGC α1 and β1 genes in rat, mouse and human. Subsequent studies found that uterine sGC mRNA and protein levels decreased rapidly in response to 17β-estradiol (estrogen) in an in vivo rat model. As early as 1 hour following treatment, mRNA levels of both sGC mRNAs decreased, and reached their lowest level of expression after 3 hours. This in vivo response was completely blocked by the pure estrogen receptor antagonist, ICI 182,780, was not seen in several other tissues examined, did not occur in response to other steroid hormones, and was due to a post-transcriptional mechanism. Additional studies ex vivo and in various cell culture models suggested that the estrogen-mediated decreased sGC mRNA expression did not require signals from other tissues, but may require cell communication or paracrine factors between different cell types within the uterus. Using chemical inhibitors and molecular targeting in other related studies, it was revealed that c-Jun-N-terminal kinase (JNK) signaling was responsible for decreased sGC mRNA expression in rat PC12 and RFL-6 cells, two models previously determined to exhibit rapid decreased sGC mRNA expression in response to different stimuli. To further investigate the post-transcriptional gene regulation, the full length sGCα1 3′-untranslated region (3′UTR) was cloned from rat uterine tissue and ligated downstream of the rabbit β-globin gene and expressed as a chimeric mRNA in the rat PC12 and RFL-6 cell models. Expression studies with the chimeric mRNA showed that the sGCα 1 3′UTR was not sufficient to mediate the post-transcriptional regulation of its mRNA by JNK or cAMP signaling in PC12 and RFL-6 cells. This study has provided numerous valuable tools for future studies involving the molecular regulation of the sGC genes. Importantly, the present results identified a novel paradigm and a previously unknown signaling pathway for sGC mRNA regulation that could potentially be exploited to treat diseases such as uterine cancers, neuronal disorders, hypertension or various inflammatory conditions. ^

Regulation and transcript analysis of the ceroperon responsible for quorum sensing in Rhodobacter sphaeroides 2.4.1

Rhodobacter sphaeroides 2.4.1 is a Gram negative facultative photoheterotrophic bacterium that has been shown to have an N-acyl homoserine lactone-based quorum sensing system called cer for c&barbelow;ommunity e&barbelow;scape r&barbelow;esponse. The cer ORFs are cerR, the transcriptional regulator, cerI, the autoinducer synthase and cerA , whose function is unknown. The autoinducer molecule, 7,8- cis-N-(tetradecenoyl) homoserine lactone, has been characterized. The objective of this study was to identify an environmental stimulus that influences the regulation of cerRAI and, to characterize transcription of the cer operon. ^ A cerR::lacZ transcriptional fusion was made and β-Galactosidase assays were performed in R. sphaeroides 2.4.1 strains, wild type, AP3 (CerI−) and AP4 (CerR−). The cerR::lacZ β-Galactosidase assays were used as an initial survey of the mode of regulation of the Cer system. A cerA::lacZ translational fusion was created and was used to show that cerA can be translated. The presence of 7,8-cis-N-(tetradecenoyl) homoserine lactone was detected from R. sphaeroides strains wild type and AP4 (CerR−) using a lasR::lacZ translational fusion autoinducer bioassay. The cerR::lacZ transcriptional fusion in R. sphaeroides 2.4.1 wild type was tested under different environmental stimuli, such as various carbon sources, oxygen tensions, light intensities and culture media to determine if they influence transcription of the cer ORFs. Although lacZ assay data implicated high light intensity at 100 W/m2 to stimulate cer transcription, quantitative Northern RNA data of the cerR transcript showed that low light intensity at 3 W/m2 is at least one environmental stimulus that induces cer transcription. This finding was supported by DNA microarray analysis. Northern analysis of the cerRAI transcript provided evidence that the cer ORFs are co-transcribed, and that the cer operon contains two additional genes. Bioinformatics was used to identify genes that may be regulated by the Cer system by identifying putative lux box homologue sequences in the presumed promoter region of these genes. Genes that were identified were fliQ, celB and calsymin, all implicated in interacting with plants. Primer extension was used to help localize cis-elements in the promoter region. The cerR::lacZ transcriptional fusion was monitored in a subset of different global DNA binding transcriptional regulator mutant strains of R. sphaeroides 2.4.1. Those regulators involved in maintaining an anaerobic photosynthetic lifestyle appeared to have an effect. Collectively, the data imply that R. sphaeroides 2.4.1 activates the Cer system when grown anaerobic photosynthetically at low light intensity, 3 W/m2, and it may be involved in an interaction with plants. ^

An extracytoplasmic function sigma factor operon regulates Myxococcus xanthus developmental gene expression

Myxococcus xanthus is a Gram-negative soil bacterium that undergoes multicellular development when high-density cells are starved on a solid surface. Expression of the 4445 gene, predicted to encode a periplasmic protein, commences 1.5 h after the initiation of development and requires starvation and high density conditions. Addition of crude or boiled supernatant from starving high-density cells restored 4445 expression to starving low-density cells. Addition of L-threonine or L-isoleucine to starving low-density cells also restored 4445 expression, indicating that the high-density signaling activity present in the supernatant might be composed of extracellular amino acids or small peptides. To investigate the circuitry integrating these starvation and high-density signals, the cis- and trans-acting elements controlling 4445 expression were identified. The 4445 transcription start site was determined by primer extension analysis to be 58 by upstream of the predicted translation start site. The promoter region contained a consensus sequence characteristic of e&barbelow;xtrac&barbelow;ytoplasmic f&barbelow;unction (ECF) sigma factor-dependent promoters, suggesting that 4445 expression might be regulated by an ECF sigma factor-dependent pathway, which are known to respond to envelope stresses. The small size of the minimum regulatory region, identified by 5′-end deletion analysis as being only 66 by upstream of the transcription start site, suggests that RNA polymerase could be the sole direct regulator of 4445 expression. To identify trans-acting negative regulators of 4445 expression, a strain containing a 4445-lacZ was mutagenized using the Himar1-tet transposon. The four transposon insertions characterized mapped to an operon encoding a putative ECF sigma factor, ecfA; an anti-sigma factor, reaA; and a negative regulator, reaB. The reaA and the reaB mutants expressed 4445 during growth and development at levels almost 100-fold higher than wild type, indicating that these genes encode negative regulators. The ecfA mutant expressed 4445-lacZ at basal levels, indicating that ecfA is a positive regulator. High Mg2+ concentrations over-stimulated this ecfA pathway possibly due to the depletion of exopolysaccharides and assembled type IV pili. These data indicate that the ecfA operon encodes a new regulatory stress pathway that integrates and transduces starvation and cell density cues during early development and is also responsive to cell-surface alterations.^

Expression regulation of AP -2 gamma gene

AP-2γ is a member of the AP-2 transcription factor family, is highly enriched in the trophoblast cell lineage, and is essential for placenta development. In an effort to identify factors regulating AP-2γ gene expression we isolated and characterized the promoter and 5′ flanking region of the mouse and human AP-2γ genes. The transcription start site of the mouse AP-2γ gene was mapped by primer extension and 5′ RACE. Transient gene transfer studies showed that basal promoter activity resides within a highly conserved ∼200 by DNA sequence located immediately upstream of the transcription start site. The conserved region is highly GC-rich and lacks typical TATA or CCAAT boxes. Multiple potential Sp and AP-2 binding sites are clustered within this region. Electrophoretic mobility shift assays demonstrated that Sp1 and Sp3 bind to three sites in the promoter region of the mouse AP-2γ gene. Combined mutation of the three putative Sp sites reduced promoter activity by 80% in trophoblast and non-trophoblast cells, demonstrating the functional importance of these sites in AP-2γ gene expression. ^ Mutational analysis of the 5′-flanking region revealed a 117-bp positive regulatory region of the mouse AP-2γ gene located between −5700 and −5583 upstream of the transcription start site. This 117-bp positive regulatory element provided approximately 7-fold enhancement of reporter gene expression in cultured trophoblast cells. A C/EBP-Sp1 transcription factor-binding module is located in this DNA sequence. Electrophoretic mobility shift assays demonstrated that transcription factors Sp1, Sp3 and C/EBP bind to the enhancer element. Mutation of each protein-binding site reduced the enhanced expression significantly. Mutagenesis assays showed that two other protein-binding sites also contribute to the enhancer activity. In summary, we have shown that Sp1 and Sp3 bind to cis-regulatory elements located in the promoter region and contribute to basal promoter activity. We have identified a 117-bp positive regulatory element of AP-2γ gene, and we have shown that Sp and C/EBP proteins bind to the cis -regulatory elements and contribute to the enhanced gene expression. ^

Fluorescence energy transfer detection as a homogeneous DNA diagnostic method

A homogeneous DNA diagnostic assay based on template-directed primer extension detected by fluorescence resonance energy transfer, named template-directed dye-terminator incorporation (TDI) assay, has been developed for mutation detection and high throughput genome analysis. Here, we report the successful application of the TDI assay to detect mutations in the cystic fibrosis transmembrane conductance regulator (CFTR) gene, the human leukocyte antigen H (HLA-H) gene, and the receptor tyrosin kinase (RET) protooncogene that are associated with cystic fibrosis, hemochromatosis, and multiple endocrine neoplasia, type 2, respectively. Starting with total human DNA, the samples are amplified by the PCR followed by enzymatic degradation of excess primers and deoxyribonucleoside triphosphates before the primer extension reaction is performed. All these standardized steps are performed in the same tube, and the fluorescence changes are monitored in real time, making it a useful clinical DNA diagnostic method.

Construction of a Z-DNA-specific restriction endonuclease

Novel restriction enzymes can be created by fusing the nuclease domain of FokI endonuclease with defined DNA binding domains. Recently, we have characterized a domain (Zα) from the N-terminal region of human double-stranded RNA adenosine deaminase (hADAR1), which binds the Z-conformation with high specificity. Here we report creation of a conformation-specific endonuclease, Zα nuclease, which is a chimera of Zα and FokI nuclease. Purified Zα nuclease cleaves negatively supercoiled plasmids only when they contain a Z-DNA forming insert, such as (dC-dG)13. The precise location of the cleavage sites was determined by primer extension. Cutting has been mapped to the edge of the B-Z junction, suggesting that Zα nuclease binds within the Z-DNA insert, but cleaves in the nearby B-DNA, by using a mechanism similar to type IIs restriction enzymes. These data show that Zα binds Z-DNA in an environment similar to that in a cell. Zα nuclease, a structure-specific restriction enzyme, may be a useful tool for further study of the biological role of Z-DNA.

A nucleolar protein related to ribosomal protein L7 is required for an early step in large ribosomal subunit biogenesis

The Saccharomyces cerevisiae Rlp7 protein has extensive identity and similarity to the large ribosomal subunit L7 proteins and shares an RNA-binding domain with them. Rlp7p is not a ribosomal protein; however, it is encoded by an essential gene and therefore must perform a function essential for cell growth. In this report, we show that Rlp7p is a nucleolar protein that plays a critical role in processing of precursors to the large ribosomal subunit RNAs. Pulse–chase labeling experiments with Rlp7p-depleted cells reveal that neither 5.8SS, 5.8SL, nor 25S is produced, indicating that both the major and minor processing pathways are affected. Analysis of processing intermediates by primer extension indicates that Rlp7p-depleted cells accumulate the 27SA3 precursor RNA, which is normally the major substrate (85%) used to produce the 5.8S and 25S rRNAs, and the ratio of 27SBL to 27SBS precursors changes from approximately 1:8 to 8:1 (depleted cells). Because 27SA3 is the direct precursor to 27SBS, we conclude that Rlp7p is specifically required for the 5′ to 3′ exonucleolytic trimming of the 27SA3 into the 27SBS precursor. As it is essential for processing in both the major and minor pathways, we propose that Rlp7p may act as a specificity factor that binds precursor rRNAs and tethers the enzymes that carry out the early 5′ to 3′ exonucleolytic reactions that generate the mature rRNAs. Rlp7p may also be required for the endonucleolytic cleavage in internal transcribed spacer 2 that separates the 5.8S rRNA from the 25S rRNA.

Characterization and mutational analysis of yeast Dbp8p, a putative RNA helicase involved in ribosome biogenesis

RNA helicases of the DEAD box family are involved in almost all cellular processes involving RNA molecules. Here we describe functional characterization of the yeast RNA helicase Dbp8p (YHR169w). Our results show that Dbp8p is an essential nucleolar protein required for biogenesis of the small ribosomal subunit. In vivo depletion of Dbp8p resulted in a ribosomal subunit imbalance due to a deficit in 40S ribosomal subunits. Subsequent analyses of pre-rRNA processing by pulse–chase labeling, northern hybridization and primer extension revealed that the early steps of cleavage of the 35S precursor at sites A1 and A2 are inhibited and delayed at site A0. Synthesis of 18S rRNA, the RNA moiety of the 40S subunit, is thereby blocked in the absence of Dbp8p. The involvement of Dbp8p as a bona fide RNA helicase in ribosome biogenesis is strongly supported by the loss of Dbp8p in vivo function obtained by site-directed mutagenesis of some conserved motifs carrying the enzymatic properties of the protein family.

Identification of 10 novel snoRNA gene clusters from Arabidopsis thaliana

Ten novel small nucleolar RNA (snoRNA) gene clusters, consisting of two or three snoRNA genes, respectively, were identified from Arabidopsis thaliana. Twelve of the 25 snoRNA genes in these clusters are homologous to those of yeast and mammals according to the conserved antisense sequences that guide 2′-O-ribose methylation of rRNA. The remaining 13 snoRNA genes, including two 5.8S rRNA methylation guides, are new genes identified from A.thaliana. Interestingly, seven methylated nucleotides, predicted by novel snoRNAs Z41a–Z46, are methylated neither in yeast nor in vertebrates. Using primer extension at low dNTP concentration the six methylation sites were determined as expected. These snoRNAs were recognized as specific guides for 2′-O-ribose methylation of plant rRNAs. Z42, however, did not guide the expected methylation of 25S rRNA in our assay. Thus, its function remains to be elucidated. The intergenic spacers of the gene clusters are rich in uridine (up to 40%) and most of them range in size from 35 to 100 nt. Lack of a conserved promoter element in each spacer and the determination of polycistronic transcription from a cluster by RT–PCR assay suggest that the snoRNAs encoded in the clusters are transcribed as a polycistron under an upstream promoter, and individual snoRNAs are released after processing of the precursor. Numerous snoRNA gene clusters identified from A.thaliana and other organisms suggest that the snoRNA gene cluster is an ancient gene organization existing abundantly in plants.

Endogenous Msx1 antisense transcript: In vivo and in vitro evidences, structure, and potential involvement in skeleton development in mammals

Msx1 is a key factor for the development of tooth and craniofacial skeleton and has been proposed to play a pivotal role in terminal cell differentiation. In this paper, we demonstrated the presence of an endogenous Msx1 antisense RNA (Msx1-AS RNA) in mice, rats, and humans. In situ analysis revealed that this RNA is expressed only in differentiated dental and bone cells with an inverse correlation with Msx1 protein. These in vivo data and overexpression of Msx1 sense and AS RNA in an odontoblastic cell line (MO6-G3) showed that the balance between the levels of the two Msx1 RNAs is related to the expression of Msx1 protein. To analyze the impact of this balance in the Msx-Dlx homeoprotein pathway, we analyzed the effect of Msx1, Msx2, and Dlx5 overexpression on proteins involved in skeletal differentiation. We showed that the Msx1-AS RNA is involved in crosstalk between the Msx-Dlx pathways because its expression was abolished by Dlx5. Msx1 was shown to down-regulate a master gene of skeletal cells differentiation, Cbfa1. All these data strongly suggest that the ratio between Msx1 sense and antisense RNAs is a very important factor in the control of skeletal terminal differentiation. Finally, the initiation site for Msx1-AS RNA transcription was located by primer extension in both mouse and human in an identical region, including a consensus TATA box, suggesting an evolutionary conservation of the AS RNA-mediated regulation of Msx1 gene expression.