Symmetry in biology: from genetic code to stochastic gene regulation

Mathematical models, as instruments for understanding the workings of nature, are a traditional tool of physics, but they also play an ever increasing role in biology - in the description of fundamental processes as well as that of complex systems. In this review, the authors discuss two examples of the application of group theoretical methods, which constitute the mathematical discipline for a quantitative description of the idea of symmetry, to genetics. The first one appears, in the form of a pseudo-orthogonal (Lorentz like) symmetry, in the stochastic modelling of what may be regarded as the simplest possible example of a genetic network and, hopefully, a building block for more complicated ones: a single self-interacting or externally regulated gene with only two possible states: ` on` and ` off`. The second is the algebraic approach to the evolution of the genetic code, according to which the current code results from a dynamical symmetry breaking process, starting out from an initial state of complete symmetry and ending in the presently observed final state of low symmetry. In both cases, symmetry plays a decisive role: in the first, it is a characteristic feature of the dynamics of the gene switch and its decay to equilibrium, whereas in the second, it provides the guidelines for the evolution of the coding rules.

Gene Expression Noise in Spatial Patterning: hunchback Promoter Structure Affects Noise Amplitude and Distribution in Drosophila Segmentation

Positional information in developing embryos is specified by spatial gradients of transcriptional regulators. One of the classic systems for studying this is the activation of the hunchback (hb) gene in early fruit fly (Drosophila) segmentation by the maternally-derived gradient of the Bicoid (Bcd) protein. Gene regulation is subject to intrinsic noise which can produce variable expression. This variability must be constrained in the highly reproducible and coordinated events of development. We identify means by which noise is controlled during gene expression by characterizing the dependence of hb mRNA and protein output noise on hb promoter structure and transcriptional dynamics. We use a stochastic model of the hb promoter in which the number and strength of Bcd and Hb (self-regulatory) binding sites can be varied. Model parameters are fit to data from WT embryos, the self-regulation mutant hb(14F), and lacZ reporter constructs using different portions of the hb promoter. We have corroborated model noise predictions experimentally. The results indicate that WT (self-regulatory) Hb output noise is predominantly dependent on the transcription and translation dynamics of its own expression, rather than on Bcd fluctuations. The constructs and mutant, which lack self-regulation, indicate that the multiple Bcd binding sites in the hb promoter (and their strengths) also play a role in buffering noise. The model is robust to the variation in Bcd binding site number across a number of fly species. This study identifies particular ways in which promoter structure and regulatory dynamics reduce hb output noise. Insofar as many of these are common features of genes (e. g. multiple regulatory sites, cooperativity, self-feedback), the current results contribute to the general understanding of the reproducibility and determinacy of spatial patterning in early development.

Divergence Involving Global Regulatory Gene Mutations in an Escherichia coli Population Evolving under Phosphate Limitation

Many of the important changes in evolution are regulatory in nature. Sequenced bacterial genomes point to flexibility in regulatory circuits but we do not know how regulation is remodeled in evolving bacteria. Here, we study the regulatory changes that emerge in populations evolving under controlled conditions during experimental evolution of Escherichia coli in a phosphate-limited chemostat culture. Genomes were sequenced from five clones with different combinations of phenotypic properties that coexisted in a population after 37 days. Each of the distinct isolates contained a different mutation in 1 of 3 highly pleiotropic regulatory genes (hfq, spoT, or rpoS). The mutations resulted in dissimilar proteomic changes, consistent with the documented effects of hfq, spoT, and rpoS mutations. The different mutations do share a common benefit, however, in that the mutations each redirect cellular resources away from stress responses that are redundant in a constant selection environment. The hfq mutation lowers several individual stress responses as well the small RNA-dependent activation of rpoS translation and hence general stress resistance. The spoT mutation reduces ppGpp levels, decreasing the stringent response as well as rpoS expression. The mutations in and upstream of rpoS resulted in partial or complete loss of general stress resistance. Our observations suggest that the degeneracy at the core of bacterial stress regulation provides alternative solutions to a common evolutionary challenge. These results can explain phenotypic divergence in a constant environment and also how evolutionary jumps and adaptive radiations involve altered gene regulation.

MicroRNA transcriptome profiles during swine skeletal muscle development

Background: MicroRNA (miR) are a class of small RNAs that regulate gene expression by inhibiting translation of protein encoding transcripts. To evaluate the role of miR in skeletal muscle of swine, global microRNA abundance was measured at specific developmental stages including proliferating satellite cells, three stages of fetal growth, day-old neonate, and the adult. Results: Twelve potential novel miR were detected that did not match previously reported sequences. In addition, a number of miR previously reported to be expressed in mammalian muscle were detected, having a variety of abundance patterns through muscle development. Muscle-specific miR-206 was nearly absent in proliferating satellite cells in culture, but was the highest abundant miR at other time points evaluated. In addition, miR-1 was moderately abundant throughout developmental stages with highest abundance in the adult. In contrast, miR-133 was moderately abundant in adult muscle and either not detectable or lowly abundant throughout fetal and neonate development. Changes in abundance of ubiquitously expressed miR were also observed. MiR-432 abundance was highest at the earliest stage of fetal development tested (60 day-old fetus) and decreased throughout development to the adult. Conversely, miR-24 and miR-27 exhibited greatest abundance in proliferating satellite cells and the adult, while abundance of miR-368, miR-376, and miR-423-5p was greatest in the neonate. Conclusion: These data present a complete set of transcriptome profiles to evaluate miR abundance at specific stages of skeletal muscle growth in swine. Identification of these miR provides an initial group of miR that may play a vital role in muscle development and growth.

ProbFAST: Probabilistic Functional Analysis System Tool

Background: The post-genomic era has brought new challenges regarding the understanding of the organization and function of the human genome. Many of these challenges are centered on the meaning of differential gene regulation under distinct biological conditions and can be performed by analyzing the Multiple Differential Expression (MDE) of genes associated with normal and abnormal biological processes. Currently MDE analyses are limited to usual methods of differential expression initially designed for paired analysis. Results: We proposed a web platform named ProbFAST for MDE analysis which uses Bayesian inference to identify key genes that are intuitively prioritized by means of probabilities. A simulated study revealed that our method gives a better performance when compared to other approaches and when applied to public expression data, we demonstrated its flexibility to obtain relevant genes biologically associated with normal and abnormal biological processes. Conclusions: ProbFAST is a free accessible web-based application that enables MDE analysis on a global scale. It offers an efficient methodological approach for MDE analysis of a set of genes that are turned on and off related to functional information during the evolution of a tumor or tissue differentiation. ProbFAST server can be accessed at http://gdm.fmrp.usp.br/probfast.

Differential expression of microRNAs in mouse pain models

Background: MicroRNAs (miRNAs) are short non-coding RNAs that inhibit translation of target genes by binding to their mRNAs. The expression of numerous brain-specific miRNAs with a high degree of temporal and spatial specificity suggests that miRNAs play an important role in gene regulation in health and disease. Here we investigate the time course gene expression profile of miR-1, -16, and -206 in mouse dorsal root ganglion (DRG), and spinal cord dorsal horn under inflammatory and neuropathic pain conditions as well as following acute noxious stimulation. Results: Quantitative real-time polymerase chain reaction analyses showed that the mature form of miR-1, -16 and -206, is expressed in DRG and the dorsal horn of the spinal cord. Moreover, CFA-induced inflammation significantly reduced miRs-1 and -16 expression in DRG whereas miR-206 was downregulated in a time dependent manner. Conversely, in the spinal dorsal horn all three miRNAs monitored were upregulated. After sciatic nerve partial ligation, miR-1 and -206 were downregulated in DRG with no change in the spinal dorsal horn. On the other hand, axotomy increases the relative expression of miR-1, -16, and 206 in a time-dependent fashion while in the dorsal horn there was a significant downregulation of miR-1. Acute noxious stimulation with capsaicin also increased the expression of miR-1 and -16 in DRG cells but, on the other hand, in the spinal dorsal horn only a high dose of capsaicin was able to downregulate miR-206 expression. Conclusions: Our results indicate that miRNAs may participate in the regulatory mechanisms of genes associated with the pathophysiology of chronic pain as well as the nociceptive processing following acute noxious stimulation. We found substantial evidence that miRNAs are differentially regulated in DRG and the dorsal horn of the spinal cord under different pain states. Therefore, miRNA expression in the nociceptive system shows not only temporal and spatial specificity but is also stimulus-dependent.

RNA interference-mediated knockdown of CD49e (alpha 5 integrin chain) in human thymic epithelial cells modulates the expression of multiple genes and decreases thymocyte adhesion

Background: The thymus is a central lymphoid organ, in which bone marrow-derived T cell precursors undergo a complex process of maturation. Developing thymocytes interact with thymic microenvironment in a defined spatial order. A component of thymic microenvironment, the thymic epithelial cells, is crucial for the maturation of T-lymphocytes through cell-cell contact, cell matrix interactions and secretory of cytokines/chemokines. There is evidence that extracellular matrix molecules play a fundamental role in guiding differentiating thymocytes in both cortical and medullary regions of the thymic lobules. The interaction between the integrin alpha 5 beta 1 (CD49e/CD29; VLA-5) and fibronectin is relevant for thymocyte adhesion and migration within the thymic tissue. Our previous results have shown that adhesion of thymocytes to cultured TEC line is enhanced in the presence of fibronectin, and can be blocked with anti-VLA-5 antibody. Results: Herein, we studied the role of CD49e expressed by the human thymic epithelium. For this purpose we knocked down the CD49e by means of RNA interference. This procedure resulted in the modulation of more than 100 genes, some of them coding for other proteins also involved in adhesion of thymocytes; others related to signaling pathways triggered after integrin activation, or even involved in the control of F-actin stress fiber formation. Functionally, we demonstrated that disruption of VLA-5 in human TEC by CD49e-siRNA-induced gene knockdown decreased the ability of TEC to promote thymocyte adhesion. Such a decrease comprised all CD4/CD8-defined thymocyte subsets. Conclusion: Conceptually, our findings unravel the complexity of gene regulation, as regards key genes involved in the heterocellular cell adhesion between developing thymocytes and the major component of the thymic microenvironment, an interaction that is a mandatory event for proper intrathymic T cell differentiation.

Cell homeostasis in a Leishmania major mutant overexpressing the spliced leader RNA is maintained by an increased proteolytic activity

Although several stage-specific genes have been identified in Leishmania, the molecular mechanisms governing developmental gene regulation in this organism are still not well understood. We have previously reported an attenuation of virulence in Leishmania major and L braziliensis carrying extra-copies of the spliced leader RNA gene. Here, we surveyed the major differences in proteome and transcript expression profiles between the spliced leader RNA overexpressor and control lines using two-dimensional gel electrophoresis and differential display reverse transcription PCR, respectively. Thirty-nine genes related to stress response, cytoskeleton, proteolysis, cell cycle control and proliferation, energy generation, gene transcription, RNA processing and post-transcriptional regulation have abnormal patterns of expression in the spliced leader RNA overexpressor line. The evaluation of proteolytic pathways in the mutant revealed a selective increase of cysteine protease activity and an exacerbated ubiquitin-labeled protein population. Polysome profile analysis and measurement of cellular protein aggregates showed that protein translation in the spliced leader RNA overexpressor line is increased when compared to the control line. We found that L major promastigotes maintain homeostasis in culture when challenged with a metabolic imbalance generated by spliced leader RNA surplus through modulation of intracellular proteolysis. However, this might interfere with a fine-tuned gene expression control necessary for the amastigote multiplication in the mammalian host. (c) 2010 Elsevier Ltd. All rights reserved.

Myosin Va phosphorylated on ser(1650) is found in nuclear speckles and redistributes to nucleoli upon inhibition of transcription

Nuclear actin and nuclear myosins have been implicated in the regulation of geneexpression in vertebrate cells. Myosin V is a class of actin-based motor proteins involved in cytoplasmic vesicle transport and anchorage, spindle-pole alignment and mRNA translocation. In this study, myosin-Va, phosphorylated on a conserved serine in the tail domain (phospho-ser(1650) MVa), was localized to subnuclear compartments. A monoclonal antibody, 9E6, raised against a peptide corresponding to phosphoserine(1650) and flanking regions of the murine myosin Va sequence, was immunoreactive to myosin Va heavy chain in cellular and nuclear extracts of HeLa cells, PC12 cells and B16-F10 melanocytes. Immunofluorescence microscopy with this antibody revealed discrete irregular spots within the nucleoplasm that colocalized with SC35, a splicing factor that earmarks nuclear speckles. Phospho-ser(1650) MVa was not detected in other nuclear compartments, such as condensed chromatin, Cajal bodies, gems and perinucleolar caps. Although nucleoli also were not labeled by 9E6 under normal conditions, inhibition of transcription in HeLa cells by actinomycin D caused the redistribution of phospho-ser(1650) MVa to nucleoli, as well as separating a fraction of phosphoser(1650) MVa from SC35 into near-neighboring particles. These observations indicate a novel role for myosin Va in nuclear compartmentalization and offer a new lead towards the understanding of actomyosin-based gene regulation.

Transcriptional profiling reveals genes in the human pathogen Trichophyton rubrum that are expressed in response to pH signaling

Trichophyton rubrum is a dermatophyte that infects human skin and nails. Its growth on keratin as its carbon source shifts the ambient pH from acidic to alkaline, which may be an efficient strategy for its successful infection and maintenance in the host. In this study, we used suppression subtractive hybridization to identify genes preferentially expressed in T rubrum incubated at either pH 5.0 or pH 8.0. The functional grouping of the 341 overexpressed unigenes indicated proteins putatively involved in diverse cellular processes, such as membrane remodeling, cellular transport, metabolism, cellular protection, fungal pathogenesis, gene regulation, interaction with the environment, and iron uptake. Although the basic metabolic machinery identified under both growth conditions seems to be functionally similar, distinct genes are upregulated at acidic or alkaline pHs. We also isolated a large number of genes of unknown function, probably unique to T rubrum or dermatophytes. Interestingly, the transcriptional profiling of several genes in a pacC mutant suggests that, in T rubrum, the transcription factor PacC has a diversity of metabolic functions, in response to either acidic or alkaline ambient pH. (C) 2009 Elsevier Ltd. All rights reserved.

Fur controls iron homeostasis and oxidative stress defense in the oligotrophic alpha-proteobacterium Caulobacter crescentus

In most bacteria, the ferric uptake regulator (Fur) is a global regulator that controls iron homeostasis and other cellular processes, such as oxidative stress defense. In this work, we apply a combination of bioinformatics, in vitro and in vivo assays to identify the Caulobacter crescentus Fur regulon. A C. crescentus fur deletion mutant showed a slow growth phenotype, and was hypersensitive to H(2)O(2) and organic peroxide. Using a position weight matrix approach, several predicted Fur-binding sites were detected in the genome of C. crescentus, located in regulatory regions of genes not only involved in iron uptake and usage but also in other functions. Selected Fur-binding sites were validated using electrophoretic mobility shift assay and DNAse I footprinting analysis. Gene expression assays revealed that genes involved in iron uptake were repressed by iron-Fur and induced under conditions of iron limitation, whereas genes encoding iron-using proteins were activated by Fur under conditions of iron sufficiency. Furthermore, several genes that are regulated via small RNAs in other bacteria were found to be directly regulated by Fur in C. crescentus. In conclusion, Fur functions as an activator and as a repressor, integrating iron metabolism and oxidative stress response in C. crescentus.

Pseudomonas aeruginosa PA14 cupD transcription is activated by the RcsB response regulator, but repressed by its putative cognate sensor RcsC

The opportunistic pathogen Pseudomonas aeruginosa PA14 possesses four fimbrial cup clusters, which may confer the ability to adapt to different environments. cupD lies in the pathogenicity island PAPI-1 next to genes coding for a putative phosphorelay system composed of the hybrid histidine kinase RcsC and the response regulator RcsB. The main focus of this work was the regulation of cupD at the mRNA level. It was found that the HN-S-like protein MvaT does not exert a strong influence on cupD transcript levels, as it does for cupA. cupD transcription is higher in cultures grown at 28 degrees C, which agrees with a cupD mutant presenting attenuated virulence only in a plant model, but not in a mouse model of infection. Whereas an rcsC in-frame deletion mutant presented higher levels of cupD mRNA, rcsB deletion had the opposite effect. Accordingly, overexpression of RcsB increased the levels of cupD transcription, and promoted biofilm formation and the appearance of fimbriae. A single transcription start site was determined for cupD and transcription from this site was induced by RcsB. A motif similar to the enterobacterial RcsB/RcsA-binding site was detected adjacent to the -35 region, suggesting that this could be the RcsB-binding site. Comparison of P. aeruginosa and Escherichia coli Rcs may provide insights into how similar systems can be used by different bacteria to control gene expression and to adapt to various environmental conditions.

Down-regulation of the candidate tumor suppressor gene PAR-4 is associated with poor prognosis in breast cancer

Substantial experimental evidence indicates that PAWR gene (PKC apoptosis WT1 regulator; also named PAR-4, prostate apoptosis response-4) is a central player in cancer cell survival and a potential target for cancer-selective targeted therapeutics. However, little is known about the role of PAR-4 in breast cancer. We investigated the possible role of PAR-4 expression in breast cancer. IHC results on tissue microarrays containing 1,161 primary breast tumor samples showed that 57% (571/995) of analyzable cases were negative for PAR-4 nuclear staining. Down-regulation of nuclear PAR-4 protein expression predicted a poor prognosis for breast cancer patients (OS; P=0.041, log-rank test). PAR-4 down-regulation also correlates with poor survival in the group of patients with luminal A subtype breast cancer (P=0.028). Additionally, in this large series of breast cancer patients, we show that ERBB2/HER2, EGFR and pAKT protein expression are significantly associated with shorter disease-free survival and overall survival, but the prognosis was even worse for HER2-positive, EGFR-positive or pAKT-positive breast cancer patients with tumors negative for nuclear PAR-4 expression. Furthermore, using three-dimensional (3D) cell culture we provide preliminary results showing that PAR-4 is highly expressed in the MCF10A cells inside the acini structure, suggesting that PAR-4 might have a role in the lumen acini formation. Taken together, our results provide, for the first time, evidence that PAR-4 may have a role in the process of the mammary eland morphogenesis and its functional inactivation is associated with tumor aggressive phenotype and might represent an additional prognostic and predictive marker for breast cancer.

Gene expression down-regulation in CD90(+) prostate tumor-associated stromal cells involves potential organ-specific genes

Background: The prostate stroma is a key mediator of epithelial differentiation and development, and potentially plays a role in the initiation and progression of prostate cancer. The tumor-associated stroma is marked by increased expression of CD90/THYI. Isolation and characterization of these stromal cells could provide valuable insight into the biology of the tumor microenvironment. Methods: Prostate CD90(+) stromal fibromuscular cells from tumor specimens were isolated by cell-sorting and analyzed by DNA microarray. Dataset analysis was used to compare gene expression between histologically normal and tumor-associated stromal cells. For comparison, stromal cells were also isolated and analyzed from the urinary bladder. Results: The tumor-associated stromal cells were found to have decreased expression of genes involved in smooth muscle differentiation, and those detected in prostate but not bladder. Other differential expression between the stromal cell types included that of the CXC-chemokine genes. Conclusion: CD90(+) prostate tumor-associated stromal cells differed from their normal counterpart in expression of multiple genes, some of which are potentially involved in organ development.

Global gene expression under nitrogen starvation in Xylella fastidiosa: contribution of the sigma(54) regulon

Background: Xylella fastidiosa, a Gram-negative fastidious bacterium, grows in the xylem of several plants causing diseases such as citrus variegated chlorosis. As the xylem sap contains low concentrations of amino acids and other compounds, X. fastidiosa needs to cope with nitrogen limitation in its natural habitat. Results: In this work, we performed a whole-genome microarray analysis of the X. fastidiosa nitrogen starvation response. A time course experiment (2, 8 and 12 hours) of cultures grown in defined medium under nitrogen starvation revealed many differentially expressed genes, such as those related to transport, nitrogen assimilation, amino acid biosynthesis, transcriptional regulation, and many genes encoding hypothetical proteins. In addition, a decrease in the expression levels of many genes involved in carbon metabolism and energy generation pathways was also observed. Comparison of gene expression profiles between the wild type strain and the rpoN null mutant allowed the identification of genes directly or indirectly induced by nitrogen starvation in a sigma(54)-dependent manner. A more complete picture of the sigma(54) regulon was achieved by combining the transcriptome data with an in silico search for potential sigma(54)-dependent promoters, using a position weight matrix approach. One of these sigma(54)-predicted binding sites, located upstream of the glnA gene (encoding glutamine synthetase), was validated by primer extension assays, confirming that this gene has a sigma(54)-dependent promoter. Conclusions: Together, these results show that nitrogen starvation causes intense changes in the X. fastidiosa transcriptome and some of these differentially expressed genes belong to the sigma(54) regulon.