Avoiding transcription factor competition at promoter level increases the chances of obtaining oscillations

Background: The ultimate goal of synthetic biology is the conception and construction of genetic circuits that are reliable with respect to their designed function (e.g. oscillators, switches). This task remains still to be attained due to the inherent synergy of the biological building blocks and to an insufficient feedback between experiments and mathematical models. Nevertheless, the progress in these directions has been substantial. Results: It has been emphasized in the literature that the architecture of a genetic oscillator must include positive (activating) and negative (inhibiting) genetic interactions in order to yield robust oscillations. Our results point out that the oscillatory capacity is not only affected by the interaction polarity but by how it is implemented at promoter level. For a chosen oscillator architecture, we show by means of numerical simulations that the existence or lack of competition between activator and inhibitor at promoter level affects the probability of producing oscillations and also leaves characteristic fingerprints on the associated period/amplitude features. Conclusions: In comparison with non-competitive binding at promoters, competition drastically reduces the region of the parameters space characterized by oscillatory solutions. Moreover, while competition leads to pulse-like oscillations with long-tail distribution in period and amplitude for various parameters or noisy conditions, the non-competitive scenario shows a characteristic frequency and confined amplitude values. Our study also situates the competition mechanism in the context of existing genetic oscillators, with emphasis on the Atkinson oscillator.

Pseudogenes in the ENCODE regions: consensus annotation, analysis of transcription, and evolution

Arising from either retrotransposition or genomic duplication of functional genes, pseudogenes are “genomic fossils” valuable for exploring the dynamics and evolution of genes and genomes. Pseudogene identification is an important problem in computational genomics, and is also critical for obtaining an accurate picture of a genome’s structure and function. However, no consensus computational scheme for defining and detecting pseudogenes has been developed thus far. As part of the ENCyclopedia Of DNA Elements (ENCODE) project, we have compared several distinct pseudogene annotation strategies and found that different approaches and parameters often resulted in rather distinct sets of pseudogenes. We subsequently developed a consensus approach for annotating pseudogenes (derived from protein coding genes) in the ENCODE regions, resulting in 201 pseudogenes, two-thirds of which originated from retrotransposition. A survey of orthologs for these pseudogenes in 28 vertebrate genomes showed that a significant fraction (∼80%) of the processed pseudogenes are primate-specific sequences, highlighting the increasing retrotransposition activity in primates. Analysis of sequence conservation and variation also demonstrated that most pseudogenes evolve neutrally, and processed pseudogenes appear to have lost their coding potential immediately or soon after their emergence. In order to explore the functional implication of pseudogene prevalence, we have extensively examined the transcriptional activity of the ENCODE pseudogenes. We performed systematic series of pseudogene-specific RACE analyses. These, together with complementary evidence derived from tiling microarrays and high throughput sequencing, demonstrated that at least a fifth of the 201 pseudogenes are transcribed in one or more cell lines or tissues.

The DART classification of unannotated transcription within the ENCODE regions: associating transcription with known and novel loci

For the ∼1% of the human genome in the ENCODE regions, only about half of the transcriptionally active regions (TARs) identified with tiling microarrays correspond to annotated exons. Here we categorize this large amount of “unannotated transcription.” We use a number of disparate features to classify the 6988 novel TARs—array expression profiles across cell lines and conditions, sequence composition, phylogenetic profiles (presence/absence of syntenic conservation across 17 species), and locations relative to genes. In the classification, we first filter out TARs with unusual sequence composition and those likely resulting from cross-hybridization. We then associate some of those remaining with proximal exons having correlated expression profiles. Finally, we cluster unclassified TARs into putative novel loci, based on similar expression and phylogenetic profiles. To encapsulate our classification, we construct a Database of Active Regions and Tools (DART.gersteinlab.org). DART has special facilities for rapidly handling and comparing many sets of TARs and their heterogeneous features, synchronizing across builds, and interfacing with other resources. Overall, we find that ∼14% of the novel TARs can be associated with known genes, while ∼21% can be clustered into ∼200 novel loci. We observe that TARs associated with genes are enriched in the potential to form structural RNAs and many novel TAR clusters are associated with nearby promoters. To benchmark our classification, we design a set of experiments for testing the connectivity of novel TARs. Overall, we find that 18 of the 46 connections tested validate by RT-PCR and four of five sequenced PCR products confirm connectivity unambiguously.

Prominent use of distal 5’ transcription start sites and discovery of a large number of additional exons in ENCODE regions

This report presents systematic empirical annotation of transcript products from 399 annotated protein-coding loci across the 1% of the human genome targeted by the Encyclopedia of DNA elements (ENCODE) pilot project using a combination of 5' rapid amplification of cDNA ends (RACE) and high-density resolution tiling arrays. We identified previously unannotated and often tissue- or cell-line-specific transcribed fragments (RACEfrags), both 5' distal to the annotated 5' terminus and internal to the annotated gene bounds for the vast majority (81.5%) of the tested genes. Half of the distal RACEfrags span large segments of genomic sequences away from the main portion of the coding transcript and often overlap with the upstream-annotated gene(s). Notably, at least 20% of the resultant novel transcripts have changes in their open reading frames (ORFs), most of them fusing ORFs of adjacent transcripts. A significant fraction of distal RACEfrags show expression levels comparable to those of known exons of the same locus, suggesting that they are not part of very minority splice forms. These results have significant implications concerning (1) our current understanding of the architecture of protein-coding genes; (2) our views on locations of regulatory regions in the genome; and (3) the interpretation of sequence polymorphisms mapping to regions hitherto considered to be "noncoding," ultimately relating to the identification of disease-related sequence alterations.

Tandem chimeric transcription as a means to increase protein diversity in the human genome

The “one-gene, one-protein” rule, coined by Beadle and Tatum, has been fundamental to molecular biology. The rule implies that the genetic complexity of an organism depends essentially on its gene number. The discovery, however, that alternative gene splicing and transcription are widespread phenomena dramatically altered our understanding of the genetic complexity of higher eukaryotic organisms; in these, a limited number of genes may potentially encode a much larger number of proteins. Here we investigate yet another phenomenon that may contribute to generate additional protein diversity. Indeed, by relying on both computational and experimental analysis, we estimate that at least 4%–5% of the tandem gene pairs in the human genome can be eventually transcribed into a single RNA sequence encoding a putative chimeric protein. While the functional significance of most of these chimeric transcripts remains to be determined, we provide strong evidence that this phenomenon does not correspond to mere technical artifacts and that it is a common mechanism with the potential of generating hundreds of additional proteins in the human genome.

Transcription factor map alignment of promoter regions

We address the problem of comparing and characterizing the promoter regions of genes with similar expression patterns. This remains a challenging problem in sequence analysis, because often the promoter regions of co-expressed genes do not show discernible sequence conservation. In our approach, thus, we have not directly compared the nucleotide sequence of promoters. Instead, we have obtained predictions of transcription factor binding sites, annotated the predicted sites with the labels of the corresponding binding factors, and aligned the resulting sequences of labels—to which we refer here as transcription factor maps (TF-maps). To obtain the global pairwise alignment of two TF-maps, we have adapted an algorithm initially developed to align restriction enzyme maps. We have optimized the parameters of the algorithm in a small, but well-curated, collection of human–mouse orthologous gene pairs. Results in this dataset, as well as in an independent much larger dataset from the CISRED database, indicate that TF-map alignments are able to uncover conserved regulatory elements, which cannot be detected by the typical sequence alignments.

Role of transcription factor PROX1 in colon cancer metastasis

Initiation and progression of most colorectal cancers (CRCs) are driven by hyper-activation of the canonical Wnt/ß-catenin/TCF signaling pathway. However, a basal level of activation of this pathway is necessary for intestinal cell homeostasis; thus only CRC-specific effectors of this pathway could be exploited as potential clinical targets. PROX1 is an evolutionary conserved transcription factor with multiple roles in several tissues in embryogenesis, and increasing relevance in cancer. PROX1 is a colon cancer-specific Wnt target in the intestine, thus it might represent a therapeutic target. The role of PROX1 in promoting the transition from early to highly-dysplastic adenoma was previously described [1], Importantly, tumor metastasis is a leading cause of cancer-related mortality. Frequently, micrometastases are already present in patients at the time of diagnosis, therefore better understanding of the mechanisms regulating growth of macrometastatic lesions is important for the development of novel treatment approaches. In this study we showed that PROX1 is expressed in colon cancer stem cell and promotes the outgrowth of metastatic lesions. Firstly, we analyzed the expression of PROX1 in advanced CRCs and their metastases. We found that PROX1 over-expression is a feature of microsatellite stable tumors (~85% of microsatellite stable (MSS) CRCs), which generally have worse prognosis in comparison to microsatellite unstable CRCs. Analysis of primary CRCs and corresponding metastatic lesions showed that PROX1 expression is conserved, or increased in metastases. Further bioinformatics analysis of tumor and metastases gene expression profiles showed that PROX1 is co- expressed with stem cell and progenitor markers. Moreover, in inducible ApcmLgr5-EGFP-lres-CreERT2 model, Prox1+ cells marked a sub-population of Lgr5+ stem cells and subsequent transient amplifying cell population. Orthotopic model of CRC and lung colonization assays in mice demonstrated that PROX1 promotes tumor cell outgrowth in metastatic lesions, while it has no effect on primary tumor growth, invasion, and survival in circulation or cell extravasation. In vitro, PROX1 expressing tumor cells demonstrated strongly increased capacity to form spheroids, and increased survival and proliferation under hypoxic or nutrient-deprivation conditions. By monitoring cellular respiration under these conditions, we found that PROX1 expressing cells exhibit a better metabolic adaptation to changes in fuel source. Autophagy inhibitors, prevented growth both in vitro and in vivo of PROX1 expressing cells. Importantly, conditional inactivation of PROX1 after the establishment of metastases prevented further growth of macroscopic lesions resulting in stable disease. In summary, we identified a novel mechanism underlying the ability of metastatic colon cancer stem and progenitor cells to survive and grow in target organs through metabolic adaptation. Our results establish PROX1 as a key factor of CRC metastatic disease where it promotes survival of metastatic colon cancer stem-like cells, through their metabolic adaptation in sub-optimal microenvironments - L'initiation et la progression de la plupart des cancers colorectaux (CRC) sont entraînées par une hyper-activation de la voie métabolique Wnt/ß- caténine/TCF. Toutefois, un niveau d'activation minimal de Wnt est nécessaire pour l'homéostasie des cellules intestinales ; ainsi seuls des effecteurs spécifiques du CRC- de cette voie pourraient être exploités comme des cibles cliniques potentielles. PROX1 est un facteur de transcription évolutif conservé avec de multiples rôles dans plusieurs tissus durant l'embryogenèse et une pertinence croissante dans le cancer. PROX1 est une cible Wnt spécifique dans le cancer de l'intestin, donc il pourrait représenter une cible thérapeutique. Le rôle de PROX1 durant l'évolution de la maladie d'un stade précoce jusqu'à l'adénome hautement dysplasique a été décrit précédemment. Surtout, la métastase des tumeurs est une cause majeure de mortalité liée au cancer. Souvent, les micro-métastases sont déjà présentes chez les patients au moment du diagnostic, c'est pourquoi une meilleure compréhension des mécanismes régulant la croissance des lésions macrométastatiques est importante pour le développement de nouvelles approches thérapeutiques. Dans cette étude, nous avons prouvé que PROX1 est exprimé dans les cellules souches du cancer du côlon et favorise l'apparition de lésions métastatiques. Nous avons d'abord analysé l'expression de PROX1 dans des CRC avancés ainsi que dans leurs métastases. Nous avons constaté que la surexpression de PROX1 est une caractéristique des tumeurs stables microsatellites (~85% du MSS CRC), qui ont généralement un pronostic défavorable par rapport aux microsatellites CRC instables. L'analyse des CRC primaires et de leurs métastases liées a montré que l'expression de PROX1 est conservée, voire augmentée dans les métastases. A l'aide d'une base de données de tumeurs et métastases, nous avons observé une co- régulation de PROX1 entre cellules souches et marqueurs de progéniteurs mais pas avec des cellules différenciées. De plus, en utilisant un modèle Apcm Lgr5-EGFP-IRES-CreERT2 inductible, les cellules Prox1+ ont marqué une sous-population de cellules LGR& capable de produire une lignée. Un modèle orthotopique de cancer colorectal et des essais de colonisation du poumon chez la souris ont démontré que PROX1 favorise l'excroissance des cellules tumorales dans les lésions métastatiques, alors qu'il n'a aucun effet sur la croissance tumorale primaire, l'invasion ou une extravasation des cellules. In vitro, les cellules tumorales exprimant PROX1 ont démontré une forte augmentation de leur capacité à former des sphéroïdes, ainsi qu'une augmentation de la survie et de la prolifération dans des conditions hypoxiques ou lors de privation de nutriments. En contrôlant la respiration cellulaire dans ces conditions, nous avons constaté que les cellules exprimant PROX1 présentent une meilleure adaptation métabolique à l'évolution des sources de carburant. Des inhibiteurs de l'autophagie, suggérant une approche thérapeutique potentielle, ont tué à la fois in vitro et in vivo les cellules exprimant PROX1. Surtout, l'inactivation conditionnelle de PROX1 après l'apparition de métastases a empêché la croissance des lésions macroscopiques résultant en une maladie stable. En résumé, nous avons identifié un nouveau mécanisme mettant en évidence la capacité des cellules souches du cancer du côlon métastatique à survivre et à se développer dans les organes cibles grâce à l'adaptation métabolique. Nos résultats définissent PROX1 comme un facteur clé du cancer colorectal métastatique en favorisant la survie des cellules souches métastatiques apparentées au cancer du colon grâce à leur adaptation métabolique aux microenvironnements défavorables.