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.

Differential contribution of transcription factors to Arabidopsis thaliana defense against Spodoptera littoralis

In response to insect herbivory, Arabidopsis plants activate the synthesis of the phytohor- mone jasmonate-isoleucine, which binds to a complex consisting of the receptor COI1 and JAZ repressors. Upon proteasome-mediated JAZ degradation, basic helix-loop-helix tran- scription factors (TFs) MYC2, MYC3, and MYC4 become activated and this results in the expression of defense genes. Although the jasmonate (JA) pathway is known to be essen- tial for the massive transcriptional reprogramming that follows herbivory, there is however little information on other TFs that are required for defense against herbivores and whether they contribute significantly to JA-dependent defense gene expression. By transcriptome profiling, we identified 41TFs that were induced in response to herbivory by the generalist Spodoptera littoralis. Among them, nine genes, including WRKY18, WRKY40, ANAC019, ANAC055, ZAT10, ZAT12, AZF2, ERF13, and RRTF1, were found to play a significant role in resistance to S. littoralis herbivory. Compared to the triple mutant myc234 that is as sensitive as coi1-1 to herbivory, knockout lines of these nine TFs were only partially more sensitive to S. littoralis but, however, some displayed distinct gene expression changes at the whole-genome level. Data thus reveal that MYC2, MYC3, and MYC4 are master regu- lators of Arabidopsis resistance to a generalist herbivore and identify new genes involved in insect defense.

Transcription, tagging and coding of bilingual corpora via LIDES

L'objectiu d'aquest informe és presentar l'aplicació d'una sèrie de propostes sobre transcripció, etiquetatge i codificació a dos corpus: el corpus bilingüe LC (La Canonja (Català-Espanyol)) i el corpus trilingüe CSCD (Code-switching as Communicative Design (Català-Espanyol-Anglès)). Aquestes propostes, que constitueixen l'aportació de l'equip IULA-LIPPS (Language Interaction in Plurilingual and Plurilectal Speakers) al manual de codificació del sistema LIDES (Language Interaction Database Exchange System), adoptat pel grup europeu LIPPS, poden ser útils per transcriure, etiquetar i codificar dades provinents de llengües tipològicament properes i distants.

CNOT3 is a modifier of PRPF31 mutations in retinitis pigmentosa with incomplete penetrance.

Heterozygous mutations in the PRPF31 gene cause autosomal dominant retinitis pigmentosa (adRP), a hereditary disorder leading to progressive blindness. In some cases, such mutations display incomplete penetrance, implying that certain carriers develop retinal degeneration while others have no symptoms at all. Asymptomatic carriers are protected from the disease by a higher than average expression of the PRPF31 allele that is not mutated, mainly through the action of an unknown modifier gene mapping to chromosome 19q13.4. We investigated a large family with adRP segregating an 11-bp deletion in PRPF31. The analysis of cell lines derived from asymptomatic and affected individuals revealed that the expression of only one gene among a number of candidates within the 19q13.4 interval significantly correlated with that of PRPF31, both at the mRNA and protein levels, and according to an inverse relationship. This gene was CNOT3, encoding a subunit of the Ccr4-not transcription complex. In cultured cells, siRNA-mediated silencing of CNOT3 provoked an increase in PRPF31 expression, confirming a repressive nature of CNOT3 on PRPF31. Furthermore, chromatin immunoprecipitation revealed that CNOT3 directly binds to a specific PRPF31 promoter sequence, while next-generation sequencing of the CNOT3 genomic region indicated that its variable expression is associated with a common intronic SNP. In conclusion, we identify CNOT3 as the main modifier gene determining penetrance of PRPF31 mutations, via a mechanism of transcriptional repression. In asymptomatic carriers CNOT3 is expressed at low levels, allowing higher amounts of wild-type PRPF31 transcripts to be produced and preventing manifestation of retinal degeneration.

Altered expression of the transcription factor Mef2c during retinal degeneration in Rpe65-/- mice.

Purpose. To investigate the role of the myocyte enhancer factor 2 (Mef2) transcription factor family in retinal diseases, Mef2c expression was assessed during retinal degeneration in the Rpe65(-/-) mouse model of Leber's congenital amaurosis (LCA). Mef2c-dependent expression of photoreceptor-specific genes was further addressed. Methods. Expression of Mef2 members was analyzed by oligonucleotide microarray, quantitative PCR (qPCR) and in situ hybridization. Mef2c-dependent transcriptional activity was assayed by luciferase assay in HEK293T cells. Results. Mef2c was the only Mef2 member markedly downregulated during retinal degeneration in Rpe65(-/-) mice. Mef2c mRNA level was decreased by more than 2 fold at 2 and 4 months and by 3.5 fold at 6 months in retinas of Rpe65(-/-) mice. Downregulation of Mef2c at the protein level was confirmed in Rpe65(-/-) retinas. The decrease in Mef2c mRNA levels in the developing Rpe65(-/-) retinas, from post-natal day (P)13 onward, was concomitant with the decreased expression of the rod-specific transcription factors Nrl and Nr2e3. Nrl was further shown to drive Mef2c transcriptional activity, supporting a physiological role for Mef2c in the retina. In addition, Mef2c appeared to act as a transcriptional repressor of its own expression, as well as those of the retina-specific retinal G-protein coupled receptor (Rgr), rhodopsin and M-opsin genes. Conclusions. These findings highlight the early altered regulation of the rod-specific transcriptional network in Rpe65-related disease. They further indicate that Mef2c may act as a novel transcription factor involved in the development and the maintenance of photoreceptor cells.

Histone deacetylase inhibitors repress macrophage migration inhibitory factor (MIF) expression by targeting MIF gene transcription through a local chromatin deacetylation.

The cytokine macrophage migration inhibitory factor plays a central role in inflammation, cell proliferation and tumorigenesis. Moreover, macrophage migration inhibitory factor levels correlate with tumor aggressiveness and metastatic potential. Histone deacetylase inhibitors are potent antitumor agents recently introduced in the clinic. Therefore, we hypothesized that macrophage migration inhibitory factor would represent a target of histone deacetylase inhibitors. Confirming our hypothesis, we report that histone deacetylase inhibitors of various chemical classes strongly inhibited macrophage migration inhibitory factor expression in a broad range of cell lines, in primary cells and in vivo. Nuclear run on, transient transfection with macrophage migration inhibitory factor promoter reporter constructs and transduction with macrophage migration inhibitory factor expressing adenovirus demonstrated that trichostatin A (a prototypical histone deacetylase inhibitor) inhibited endogenous, but not episomal, MIF gene transcription. Interestingly, trichostatin A induced a local and specific deacetylation of macrophage migration inhibitory factor promoter-associated H3 and H4 histones which did not affect chromatin accessibility but was associated with an impaired recruitment of RNA polymerase II and Sp1 and CREB transcription factors required for basal MIF gene transcription. Altogether, this study describes a new molecular mechanism by which histone deacetylase inhibitors inhibit MIF gene expression, and suggests that macrophage migration inhibitory factor inhibition by histone deacetylase inhibitors may contribute to the antitumorigenic effects of histone deacetylase inhibitors.

Phosphorylation regulates FOXC2-mediated transcription in lymphatic endothelial cells.

One of the key mechanisms linking cell signaling and control of gene expression is reversible phosphorylation of transcription factors. FOXC2 is a forkhead transcription factor that is mutated in the human vascular disease lymphedema-distichiasis and plays an essential role in lymphatic vascular development. However, the mechanisms regulating FOXC2 transcriptional activity are not well understood. We report here that FOXC2 is phosphorylated on eight evolutionarily conserved proline-directed serine/threonine residues. Loss of phosphorylation at these sites triggers substantial changes in the FOXC2 transcriptional program. Through genome-wide location analysis in lymphatic endothelial cells, we demonstrate that the changes are due to selective inhibition of FOXC2 recruitment to chromatin. The extent of the inhibition varied between individual binding sites, suggesting a novel rheostat-like mechanism by which expression of specific genes can be differentially regulated by FOXC2 phosphorylation. Furthermore, unlike the wild-type protein, the phosphorylation-deficient mutant of FOXC2 failed to induce vascular remodeling in vivo. Collectively, our results point to the pivotal role of phosphorylation in the regulation of FOXC2-mediated transcription in lymphatic endothelial cells and underscore the importance of FOXC2 phosphorylation in vascular development.

Reverse transcription quantitative real-time polymerase chain reaction reference genes in the spared nerve injury model of neuropathic pain: validation and literature search.

BACKGROUND: The reverse transcription quantitative real-time polymerase chain reaction (RT-qPCR) is a widely used, highly sensitive laboratory technique to rapidly and easily detect, identify and quantify gene expression. Reliable RT-qPCR data necessitates accurate normalization with validated control genes (reference genes) whose expression is constant in all studied conditions. This stability has to be demonstrated.We performed a literature search for studies using quantitative or semi-quantitative PCR in the rat spared nerve injury (SNI) model of neuropathic pain to verify whether any reference genes had previously been validated. We then analyzed the stability over time of 7 commonly used reference genes in the nervous system - specifically in the spinal cord dorsal horn and the dorsal root ganglion (DRG). These were: Actin beta (Actb), Glyceraldehyde-3-phosphate dehydrogenase (GAPDH), ribosomal proteins 18S (18S), L13a (RPL13a) and L29 (RPL29), hypoxanthine phosphoribosyltransferase 1 (HPRT1) and hydroxymethylbilane synthase (HMBS). We compared the candidate genes and established a stability ranking using the geNorm algorithm. Finally, we assessed the number of reference genes necessary for accurate normalization in this neuropathic pain model. RESULTS: We found GAPDH, HMBS, Actb, HPRT1 and 18S cited as reference genes in literature on studies using the SNI model. Only HPRT1 and 18S had been once previously demonstrated as stable in RT-qPCR arrays. All the genes tested in this study, using the geNorm algorithm, presented gene stability values (M-value) acceptable enough for them to qualify as potential reference genes in both DRG and spinal cord. Using the coefficient of variation, 18S failed the 50% cut-off with a value of 61% in the DRG. The two most stable genes in the dorsal horn were RPL29 and RPL13a; in the DRG they were HPRT1 and Actb. Using a 0.15 cut-off for pairwise variations we found that any pair of stable reference gene was sufficient for the normalization process. CONCLUSIONS: In the rat SNI model, we validated and ranked Actb, RPL29, RPL13a, HMBS, GAPDH, HPRT1 and 18S as good reference genes in the spinal cord. In the DRG, 18S did not fulfill stability criteria. The combination of any two stable reference genes was sufficient to provide an accurate normalization.

Functional interactions between the estrogen receptor and the transcription activator Sp1 regulate the estrogen-dependent transcriptional activity of the vitellogenin A1 io promoter.

Two distinct, TATA box-containing promoters regulate the transcriptional activity of the Xenopus vitellogenin A1 gene. These two promoters are of different strength and are separated by 1.8 kilobase pairs of untranslated sequence. Estrogen receptor (ER) and its ligand, 17beta-estradiol, induce the activity of both promoters. The estrogen response elements (EREs) are located proximal to the downstream i promoter while no ERE-like sequences have been identified in the vicinity of the upstream io promoter. We show here, that transcriptional activity of the upstream io promoter is Sp1-dependent. Moreover, we demonstrate that estrogen inducibility of the io promoter results from functional interactions between the io bound Sp1 and the ER bound at the proximity of i. Functional interactions between Sp1 and ER do not require the presence of a TATA box for transcriptional activation, as is demonstrated using the acyl-CoA oxidase promoter. The relative positions that ER and Sp1 occupy with respect to the initiation site determines whether these two transcription activators can synergize for transcription initiation.