Subpopulation behaviors in lactose metabolism by Streptococcus mutans.

When Streptococcus mutans is transferred from a preferred carbohydrate (glucose or fructose) to lactose, initiation of growth can take several hours, and substantial amounts of glucose are released during growth. Here, S. mutans strains UA159 and GS-5 were examined for stochastic behaviors in transcription of the lac operon. Using a gfp reporter fusion, we demonstrated that induction of the lac operon occurs in only a fraction of the population, with prior exposure to carbohydrate source and strain influencing the magniture of the sub-population response. Lower glucokinase activity in GS-5 was associated with release of substantially more glucose than UA159 and significantly lower lac expression. Mutants unable to use lactose grew on lactose as the sole carbohydrate when strains with an intact lac operon were also present in the cultures, indicative of the potential for population cheating. Utilizing a set of engineered obligate cheating and non-cheating strains, we confirmed that cheating can sustain a heterogeneous population. Futher, obligate cheaters of GS-5 competed well with the non-cheaters and showed a high degree of competitive fitness in a human-derived consortium biofilm model. The results show that bet-hedging behaviors in carbohydrate metabolism may substantially influence the composition and pathogenic potential of oral biofilms.

Relatively slow stochastic gene-state switching in the presence of positive feedback significantly broadens the region of bimodality through stabilizing the uninduced phenotypic state.

Within an isogenic population, even in the same extracellular environment, individual cells can exhibit various phenotypic states. The exact role of stochastic gene-state switching regulating the transition among these phenotypic states in a single cell is not fully understood, especially in the presence of positive feedback. Recent high-precision single-cell measurements showed that, at least in bacteria, switching in gene states is slow relative to the typical rates of active transcription and translation. Hence using the lac operon as an archetype, in such a region of operon-state switching, we present a fluctuating-rate model for this classical gene regulation module, incorporating the more realistic operon-state switching mechanism that was recently elucidated. We found that the positive feedback mechanism induces bistability (referred to as deterministic bistability), and that the parameter range for its occurrence is significantly broadened by stochastic operon-state switching. We further show that in the absence of positive feedback, operon-state switching must be extremely slow to trigger bistability by itself. However, in the presence of positive feedback, which stabilizes the induced state, the relatively slow operon-state switching kinetics within the physiological region are sufficient to stabilize the uninduced state, together generating a broadened parameter region of bistability (referred to as stochastic bistability). We illustrate the opposite phenotype-transition rate dependence upon the operon-state switching rates in the two types of bistability, with the aid of a recently proposed rate formula for fluctuating-rate models. The rate formula also predicts a maximal transition rate in the intermediate region of operon-state switching, which is validated by numerical simulations in our model. Overall, our findings suggest a biological function of transcriptional "variations" among genetically identical cells, for the emergence of bistability and transition between phenotypic states.

Folding LacZ in the periplasm of Escherichia coli.

Targeted, translational LacZ fusions provided the initial support for the signal sequence hypothesis in prokaryotes and allowed for selection of the mutations that identified the Sec translocon. Many of these selections relied on the fact that expression of targeted, translational lacZ fusions like malE-lacZ and lamB-lacZ42-1 causes lethal toxicity as folded LacZ jams the translocation pore. However, there is another class of targeted LacZ fusions that do not jam the translocon. These targeted, nonjamming fusions also show toxic phenotypes that may be useful for selecting mutations in genes involved in posttranslocational protein folding and targeting; however, they have not been investigated to the same extent as their jamming counterparts. In fact, it is still unclear whether LacZ can be fully translocated in these fusions. It may be that they simply partition into the inner membrane where they can no longer participate in folding or assembly. In the present study, we systematically characterize the nonjamming fusions and determine their ultimate localization. We report that LacZ can be fully translocated into the periplasm, where it is toxic. We show that this toxicity is likely due to LacZ misfolding and that, in the absence of the periplasmic disulfide bond catalyst DsbA, LacZ folds in the periplasm. Using the novel phenotype of periplasmic ß-galactosidase activity, we show that the periplasmic chaperone FkpA contributes to LacZ folding in this nonnative compartment. We propose that targeted, nonjamming LacZ fusions may be used to further study folding and targeting in the periplasm of Escherichia coli.

Nonlinear fitness landscape of a molecular pathway.

Genes are regulated because their expression involves a fitness cost to the organism. The production of proteins by transcription and translation is a well-known cost factor, but the enzymatic activity of the proteins produced can also reduce fitness, depending on the internal state and the environment of the cell. Here, we map the fitness costs of a key metabolic network, the lactose utilization pathway in Escherichia coli. We measure the growth of several regulatory lac operon mutants in different environments inducing expression of the lac genes. We find a strikingly nonlinear fitness landscape, which depends on the production rate and on the activity rate of the lac proteins. A simple fitness model of the lac pathway, based on elementary biophysical processes, predicts the growth rate of all observed strains. The nonlinearity of fitness is explained by a feedback loop: production and activity of the lac proteins reduce growth, but growth also affects the density of these molecules. This nonlinearity has important consequences for molecular function and evolution. It generates a cliff in the fitness landscape, beyond which populations cannot maintain growth. In viable populations, there is an expression barrier of the lac genes, which cannot be exceeded in any stationary growth process. Furthermore, the nonlinearity determines how the fitness of operon mutants depends on the inducer environment. We argue that fitness nonlinearities, expression barriers, and gene-environment interactions are generic features of fitness landscapes for metabolic pathways, and we discuss their implications for the evolution of regulation.

Integration of H-2Z1, a somatosensory cortex-expressed transgene, interferes with the expression of the Satb1 and Tbc1d5 flanking genes and affects the differentiation of a subset of cortical interneurons.

H-2Z1 is an enhancer trap transgenic mouse line in which the lacZ reporter delineates the somatosensory area of the cerebral cortex where it is expressed in a subset of layer IV neurons. In the search of somatosensory specific genes or regulatory sequences, we mapped the H-2Z1 transgene insertion site to chromosome 17, 100 and 460 kb away from Tbc1d5 and Satb1 flanking genes. We show here that insertion of the H-2Z1 transgene results in three distinct outcomes. First, a genetic background-sensitive expression of lacZ in several brain and body structures. While four genes in a 1 Mb region around the insertion are expressed in the barrel cortex, H-2Z1 expression resembles more that of its two direct neighbors. Moreover, H-2Z1 closely reports most of the body and brain expression sites of the Satb1 chromatin remodeling gene including tooth buds, thymic epithelium, pontine nuclei, fastigial cerebellar nuclei, and cerebral cortex. Second, the H-2Z1 transgene causes insertional mutagenesis of Tbc1d5 and Satb1, leading to a strong decrease in their expressions. Finally, insertion of H-2Z1 affects the differentiation of a subset of cortical GABAergic interneurons, a possible consequence of downregulation of Satb1 expression. Thus, the H-2Z1 "somatosensory" transgene is inserted in the regulatory landscape of two genes highly expressed in the developing somatosensory cortex and reports for a subdomain of their expression profiles. Together, our data suggest that regulation of H-2Z1 expression results from local and remote genetic interactions.

Abundant LacZ activity in the absence of Cre expression in the normal and inflamed synovium of adult Col2a1-Cre; ROSA26RLacZ reporter mice.

Recent analyses of Col2a1-Cre; ROSA26R reporter mice showed that synovial fibroblasts in 7-day mice were LacZ positive, due to a history of Col2a1-Cre expression conferred by their origin in the interzone of the developing joint. We have examined LacZ staining in adult Col2a1-Cre(+/0); ROSA26R(LacZ) mice, with and without inflammatory arthritis, and found that synovial fibroblasts in normal and inflamed synovium are LacZ positive, but Cre negative. Our results suggest that Cre-mediated recombination in joint interzone cells during development endure in adult synovial cells despite the absence of ongoing Cre expression. These findings have important implications and applications for the study of synovial inflammation in models of experimental arthritis.

Nrf2 degron-fused reporter system: a new tool for specific evaluation of Nrf2 inducers.

The transcription factor Nrf2 is degraded through the proteasome pathway, but is stabilized in response to oxidative and electrophilic stresses, and activates cytoprotective enzyme genes through binding to the antioxidant/electrophile response element (ARE/EpRE). Nrf2 inducers thus have considerable potential as therapeutic drugs. Although ARE-driven reporters are commonly employed to validate Nrf2 inducers, these reporters are relatively nonspecific. We have generated a new reporter, Nrf2d-LacZ, which may prove to be a better tool for validation of Nrf2 inducers. We made the Nrf2d-LacZ reporter by fusing the N-terminus of Nrf2 harboring a Neh2 degron to ß-galactosidase (LacZ), and compared its activity in immortalized mouse embryo fibroblasts (MEFs) with conventional ARE-luciferase (ARE-Luc) reporters in 293T cells, and in MEFs. Nrf2d-LacZ was degraded in unstressed conditions, but stabilized upon exposure to stresses. LacZ activity was induced by electrophiles in a dose-dependent manner, and the induction was detected much more rapidly compared with ARE-Luc. Nrf2d-LacZ was activated not only by electrophiles but also by a variety of other Nrf2 inducing stresses. Although ARE-Luc was activated by 12-O-Tetradecanoylphorbol 13-acetate in an Nrf2-independent manner, Nrf2d-LacZ was not activated by TPA, thus emphasizing the specificity of the Nrf2d-LacZ reporter system for validation of Nrf2 inducers.

Combined inadequacies of multiple B vitamins amplify colonic Wnt signaling and promote intestinal tumorigenesis in BAT-LacZxApc1638N mice.

The Wnt pathway is a pivotal signaling cascade in colorectal carcinogenesis. The purpose of this work is to determine whether depletion of folate and other metabolically related B vitamins induces in vivo activation of intestinal Wnt signaling and whether this occurs in parallel with increased tumorigenesis. A hybrid mouse was created by crossing a Wnt-reporter animal (BAT-LacZ) with a model of colorectal cancer (Apc1638N). A mild depletion of folate and vitamins B2, B6, and B12 was induced over 16 wk, and the control animals in each instance were pair fed a diet containing the basal requirement of these nutrients. The multiplicity of macroscopic tumors and aberrant crypt foci both increased by ~50% in the hybrid mice fed the depletion diet (P<0.05). A 4-fold elevation in Wnt signaling was produced by the depletion diet (P<0.05) and was accompanied by significant changes in the expression of a number of Wnt-related genes in a pattern consistent with its activation. Proliferation and apoptosis of the colonic mucosa both changed in a protransformational direction (P<0.05). In summary, mild depletion of multiple B vitamins produces in vivo activation of colonic Wnt signaling, implicating it as a key pathway by which B-vitamin inadequacies enhance intestinal tumorigenesis.

Reciprocal expression of connexin 40 and 45 during phenotypical changes in renin-secreting cells.

Gap junctional coupling of renin-producing cells is of major functional importance for the control of renin synthesis and release. This study was designed to determine the relevance of the vascular gap junction protein connexin 45 (Cx45) for the control of renin expression and secretion. By crossbreeding mice which drive Cre recombinase under the control of the endogenous renin promoter with mice harboring floxed Cx45 gene alleles, we generated viable mice with a deletion of Cx45 in the renin cell lineage. These mice were normotensive, and renin cells in their kidneys were normal with regard to localization and number. Sodium deficiency induced typical recruitment of renin-producing cells along afferent arterioles, whereas sodium overload resulted in a decrease in the number of cells expressing renin. Regulation of renin secretion by perfusion pressure, catecholamines, and angiotensin II from isolated kidneys of mice with renin cell-specific deletion of Cx45 was normal. Analyzing Cx45 promoter activity in cells of the preglomerular arteriolar tree by using mice driving the reporter gene LacZ under the control of the Cx45 promoter revealed strong staining in smooth muscle cells of the media, whereas renin-expressing cells were almost devoid of LacZ staining. Conversely, renin-producing cells, but not vascular smooth muscle cells expressed the gap junction protein Cx40. These findings suggest that Cx45 plays no major functional role in renin-producing cells, probably because the expression of Cx45 is downregulated in these cells. Since renin-producing cells in the adult kidney can reversibly transform into vascular smooth muscle cells and vice versa, our findings on connexin expression indicate that these phenotype switches are paralleled by characteristic reciprocal changes in the transcriptional activity of Cx40 and Cx45 genes.

Spatial p21 expression profile in the mid-term mouse embryo.

Throughout development cells make the decision to proliferate, arrest or die. Control of this process is essential for normal development, with unrestrained cell proliferation and cell death underlying the origin and progression of disease. The cell-cycle is tightly regulated by a number of factors including the cyclin-dependent kinase inhibitor 1A (Cdkn1a), termed p21 (or Cip1 or WAF1). p21 acts as a negative regulator of cell-cycle progression by binding and inhibiting complexes formed between the cyclin-dependent kinases and their catalytic partners the cyclins. In this report we identify the temporal spatial expression profile of p21 in the developing mid-term mouse embryo using a p21-LacZ reporter mouse line. Expression of p21 was restricted to specific regions with a correspondence to both areas of terminal differentiation and active remodelling. A complex temporal and spatial relationship between p21 expression and regions of apoptosis was evident. A protective role with regard to apoptosis for p21 is proposed.

Gfi1-Cre knock-in mouse line: A tool for inner ear hair cell-specific gene deletion.

Gfi1 encodes a zinc-finger transcription factor essential for the development and maintenance of haematopoiesis and the inner ear. In mouse inner ear, Gfi1 expression is confined to hair cells during development and in adulthood. To construct a genetic tool for inner ear hair cell-specific gene deletion, we generated a Gfi1-Cre mouse line by knocking-in Cre coding sequences into the Gfi1 locus and inactivating the endogenous Gfi1. The specificity and efficiency of Gfi1-Cre recombinase-mediated recombination in the developing inner ear was revealed through the expression of the conditional R26R-lacZ reporter gene. The onset of lacZ expression in the Gfi1(Cre/+) inner ear was first detected at E13.5 in the vestibule and at E15.5 in the cochlea, coinciding with the generation of hair cells. Throughout inner ear development, lacZ expression was detected only in hair cells. Thus, Gfi1-Cre knock-in mouse line provides a useful tool for gene manipulations specifically in inner ear hair cells.

Description of distributed features of the nestin-containing cells in brains of adult mice: a potential source of neural precursor cells.

The distribution of neural precursor cells (NPCs) in adult mice brain has so far not been described. Therefore, we investigated the distribution of NPCs by analyzing the nestin-containing cells (NCCs) in distinct brain regions of adult nestin second-intron enhancer-controlled LacZ reporter transgenic mice through LacZ staining. Results showed that NCCs existed in various regions of adult mouse brain. In cerebellum, the greatest number of NCCs existed in cortex of the simple lobule, followed by cortex of the cerebellar lobule. In olfactory bulb, NCCs were most numerous in the granular cell layer, followed by the mitral cell layer and the internal plexiform, glomerular, and external plexiform layers. In brain nuclei (nu), NCCs were most numerous in the marginal nu, followed by the brainstem and diencephalon nu. NCCs in sensory nu of brainstem were more numerous than in motor nu, and NCCs in the dorsal of sensory nu were more numerous than in the ventral part. In brain ventricle systems, NCCs were largely distributed in the center of and external to the lateral ventricle, the inferior part of the third ventricle, the dorsal and inferior parts of the fourth ventricle, and the gray matter around the cerebral aqueduct. NCCs in the left vs. right brain were not significantly different. These data collectively indicate that NCCs were extensively distributed in the cerebellum and olfactory bulb, the partial nu of the marginal system, the partial brain nu adjacent to brain ventricle systems, the subependymal zone, and the cerebral cortex around the marginal lobe and were a potential source of NPCs.

Low levels of citrin (SLC25A13) expression in adult mouse brain restricted to neuronal clusters.

The mitochondrial aspartate-glutamate carriers (AGC) aralar (SLC25A12) and citrin (SLC25A13) are components of the malate aspartate shuttle (MAS), a major intracellular pathway to transfer reducing equivalents from NADH to the mitochondrial matrix. Aralar is the main AGC isoform present in the adult brain, and it is expressed mainly in neurons. To search for the other AGC isoform, citrin, in brain glial cells, we used a citrin knockout mouse in which the lacZ gene was inserted into the citrin locus as reporter gene. In agreement with the low citrin levels known to be present in the adult mouse brain, beta-galactosidase expression was very low. Surprisingly, unlike the case with astroglial cultures that express citrin, no beta-galactosidase was found in brain glial cells. It was confined to neuronal cells within discrete neuronal clusters. Double-immunolabelling experiments showed that beta-galactosidase colocalized not with glial cell markers but with the pan-neuronal marker NeuN. The deep cerebellar nuclei and a few midbrain nuclei (reticular tegmental pontine nuclei; magnocellular red nuclei) were the regions where beta-galactosidase expression was highest, and it was up-regulated in fasted mice, as was also the case for liver beta-galactosidase. The results support the notion that glial cells have much lower AGC levels and MAS activity than neurons.

Math1 expression redefines the rhombic lip derivatives and reveals novel lineages within the brainstem and cerebellum.

The rhombic lip (RL) is an embryonic proliferative neuroepithelium that generates several groups of hindbrain neurons. However, the precise boundaries and derivatives of the RL have never been genetically identified. We use beta-galactosidase expressed from the Math1 locus in Math1-heterozygous and Math1-null mice to track RL-derived cells and to evaluate their developmental requirements for Math1. We uncover a Math1-dependent rostral rhombic-lip migratory stream (RLS) that generates some neurons of the parabrachial, lateral lemniscal, and deep cerebellar nuclei, in addition to cerebellar granule neurons. A more caudal Math1-dependent cochlear extramural stream (CES) generates the ventral cochlear nucleus and cochlear granule neurons. Similarly, mossy-fiber precerebellar nuclei require Math1, whereas the inferior olive and locus coeruleus do not. We propose that Math1 expression delimits the extent of the rhombic lip and is required for the generation of the hindbrain superficial migratory streams, all of which contribute neurons to the proprioceptive/vestibular/auditory sensory network.

A comprehensive negative regulatory program controlled by Brn3b to ensure ganglion cell specification from multipotential retinal precursors.

The retinal ganglion cells (RGCs) are the sole output neurons in the retina that form the optic nerve and convey light signals detected by photoreceptors to the higher visual system. Their degeneration and damage caused by glaucoma and injury can lead to blindness. During retinogenesis, RGCs are specified from a population of multipotential precursors capable of generating RGC, amacrine, horizontal, and cone cells. How the RGC fate is selected from these multiple neuron fates is unknown at present. Here we show that the previously unsuspected POU domain transcription factor Brn3b (brain-specific homeobox/POU domain protein 3b) plays such a critical role. Loss of Brn3b function in mice leads to misspecification of early RGC precursors as late-born RGC, amacrine, and horizontal cells, whereas misexpressed Brn3b suppresses non-RGC cell fates but promotes the RGC fate. Microarray profiling and other molecular analyses reveal that, in RGC precursors, Brn3b normally represses the expression of a network of retinogenic factor genes involved in fate commitment and differentiation of late-born RGC, amacrine, horizontal, and cone cells. Our data suggest that Brn3b specifies the RGC fate from multipotential precursors not only by promoting RGC differentiation but also by suppressing non-RGC differentiation programs as a safeguard mechanism.

Clock regulatory elements control cyclic expression of Lunatic fringe during somitogenesis.

Somitogenesis requires a segmentation clock and Notch signaling. Lunatic fringe (Lfng) expression in the presomitic mesoderm (PSM) cycles in the posterior PSM, is refined in the segmenting somite to the rostral compartment, and is required for segmentation. We identify distinct cis-acting regulatory elements for each aspect of Lfng expression. Fringe clock element 1 (FCE1) represents a conserved 110 bp region that is necessary to direct cyclic Lfng RNA expression in the posterior PSM. Mutational analysis of E boxes within FCE1 indicates a potential interplay of positive and negative transcriptional regulation by cyclically expressed bHLH proteins. A separable Lfng regulatory region directs expression to the prospective rostral aspect of the condensing somite. These independent Lfng regulatory cassettes advance a molecular framework for deciphering somite segmentation.

Separable cis-regulatory elements that contribute to tissue- and site-specific alpha 2(XI) collagen gene expression in the embryonic mouse cartilage.

Type XI collagen is a structural component of the cartilage extracellular matrix and plays an important role in skeletal morphogenesis. As a step toward defining the molecular mechanisms responsible for the regulation of type XI collagen expression, we characterized the promoter region of the mouse alpha 2(XI) collagen gene (Coll1a2). We also generated transgenic mice harboring various fragments of the promoter and the first intron of Coll1a2 linked to the Escherichia coli beta-galactosidase gene to identify the cis-acting elements responsible for tissue- and site-specific expression during development. Cloning and sequence analysis of the 5' flanking region of Coll1a2 showed that the putative 3' end of the retinoid X receptor beta gene was located 742 bp upstream of the Coll1a2 start site. This suggested that the promoter region of Coll1a2 was localized within this 742-bp sequence, which contained multiple consensus regulatory elements. Examination of the transgenic mice revealed that the longest DNA construct (containing the entire promoter and first intron sequences) directed lacZ expression in the notochord as well as in the primordial cartilage throughout the body, with the pattern of expression mimicking that of endogenous Coll1a2 transcripts. On the other hand, deletion of the upstream approximately 290 bp resulted in the elimination of lacZ expression in the primordial cartilage of the carpals, tarsals, and vertebral bodies, whereas lacZ expression in the notochord and in the other primordial cartilage elsewhere was not affected. Deletion of the first intron sequence also resulted in the loss of lacZ expression in the primordial cartilage of the carpals, tarsals, and vertebral bodies, as well as in the notochord. These results demonstrate that the upstream 742-bp and first intron segments of the mouse Coll1a2 gene contain the necessary information to confer high level tissue-specific expression in mouse embryos. In addition, our observations suggest the presence of site-specific cis-acting elements that control Coll11a2 gene expression in different cartilaginous components of the skeleton.

In vivo localization of DNA sequences and visualization of large-scale chromatin organization using lac operator/repressor recognition.

We report a new method for in situ localization of DNA sequences that allows excellent preservation of nuclear and chromosomal ultrastructure and direct, in vivo observations. 256 direct repeats of the lac operator were added to vector constructs used for transfection and served as a tag for labeling by lac repressor. This system was first characterized by visualization of chromosome homogeneously staining regions (HSRs) produced by gene amplification using a dihydrofolate reductase (DHFR) expression vector with methotrexate selection. Using electron microscopy, most HSRs showed approximately 100-nm fibers, as described previously for the bulk, large-scale chromatin organization in these cells, and by light microscopy, distinct, large-scale chromatin fibers could be traced in vivo up to 5 microns in length. Subsequent experiments demonstrated the potential for more general applications of this labeling technology. Single and multiple copies of the integrated vector could be detected in living CHO cells before gene amplification, and detection of a single 256 lac operator repeat and its stability during mitosis was demonstrated by its targeted insertion into budding yeast cells by homologous recombination. In both CHO cells and yeast, use of the green fluorescent protein-lac repressor protein allowed extended, in vivo observations of the operator-tagged chromosomal DNA. Future applications of this technology should facilitate structural, functional, and genetic analysis of chromatin organization, chromosome dynamics, and nuclear architecture.

Preferential MyoD homodimer formation demonstrated by a general method of dominant negative mutation employing fusion with a lysosomal protease.

We report on a general strategy for engineering dominant negative mutations that, in principle, requires neither extensive structural or functional knowledge of the targeted protein. The approach consists of fusing the lysosomal protease cathepsin B (CB) to a subunit of a multimeric protein. The CB fusion polypeptide can proteolytically digest the multimer and/or detour the multimer from its usual subcellular destination to the lysosome. We first demonstrate the general validity of the approach with CB fusion to E. coli lacZ, encoding tetrameric beta-galactosidase. Cotransfection of NIH 3T3 cells with a vector expressing a CB-lacZ fusion inhibits the beta-galactosidase activity produced by transfection of lacZ alone. We infer that the dominant negative inhibition results from both direct proteolysis of the beta-galactosidase tetramer by the fusion subunit and detour of the tetramer to the lysosome. In a specific application of this strategy, we have fused CB to the dimeric bHLH skeletal muscle transcription factor MyoD. The CB-MyoD fusion protein localizes to the cytoplasm, presumably the lysosome, demonstrating the dominance of lysosomal localization to nuclear localization. The CB-MyoD fusion appears to divert homodimerizing native MyoD from its usual nuclear destination, consequently inhibiting MyoD-mediated transactivation and in vitro differentiation of C2C12 myoblasts. Surprisingly, the CB-MyoD fusion fails to interact with the bHLH heterodimerization partners, E12 and E47, suggesting preferential MyoD homodimer formation, at least in the prenuclear cellular compartments.

Stability and Hopf bifurcation analysis for a Lac operon model with nonlinear degradation rate and time delay.

In this paper, we construct a discrete time delay Lac operon model with nonlinear degradation rate for mRNA, resulting from the interaction among several identical mRNA pieces. By taking a discrete time delay as bifurcation parameter, we investigate the nonlinear dynamical behaviour arising from the model, using mathematical tools such as stability and bifurcation theory. Firstly, we discuss the existence and uniqueness of the equilibrium for this system and investigate the effect of discrete delay on its dynamical behaviour. Absence or limited delay causes the system to have a stable equilibrium, which changes into a Hopf point producing oscillations if time delay is increased. These sustained oscillation are shown to be present only if the nonlinear degradation rate for mRNA satisfies specific conditions. The direction of the Hopf bifurcation giving rise to such oscillations is also determined, via the use of the so-called multiple time scales technique. Finally, numerical simulations are shown to validate and expand the theoretical analysis. Overall, our findings suggest that the degree of nonlinearity of the model can be used as a control parameter for the stabilisation of the system.