Experimental analysis and computer prediction of CTF/NFI transcription factor DNA binding sites.

Accurate prediction of transcription factor binding sites is needed to unravel the function and regulation of genes discovered in genome sequencing projects. To evaluate current computer prediction tools, we have begun a systematic study of the sequence-specific DNA-binding of a transcription factor belonging to the CTF/NFI family. Using a systematic collection of rationally designed oligonucleotides combined with an in vitro DNA binding assay, we found that the sequence specificity of this protein cannot be represented by a simple consensus sequence or weight matrix. For instance, CTF/NFI uses a flexible DNA binding mode that allows for variations of the binding site length. From the experimental data, we derived a novel prediction method using a generalised profile as a binding site predictor. Experimental evaluation of the generalised profile indicated that it accurately predicts the binding affinity of the transcription factor to natural or synthetic DNA sequences. Furthermore, the in vitro measured binding affinities of a subset of oligonucleotides were found to correlate with their transcriptional activities in transfected cells. The combined computational-experimental approach exemplified in this work thus resulted in an accurate prediction method for CTF/NFI binding sites potentially functioning as regulatory regions in vivo.

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 ( approximately 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.

Complete nucleotide sequence of Birmingham IncP alpha plasmids. Compilation and comparative analysis.

The IncP alpha promiscuous plasmid (R18, R68, RK2, RP1 and RP4) comprises 60,099 bp of nucleotide sequence, encoding at least 74 genes. About 40 kb of the genome, designated the IncP core and including all essential replication and transfer functions, can be aligned with equivalent sequences in the IncP beta plasmid R751. The compiled IncP alpha sequence revealed several previously unidentified reading frames that are potential genes. IncP alpha plasmids carry genetic information very efficiently: the coding sequences of the genes are closely packed but rarely overlap, and occupy almost 86% of the genome's nucleotide sequence. All of the 74 genes should be expressed, although there is as yet experimental evidence for expression of only 60 of them. Six examples of tandem-in-frame initiation sites specifying two gene products each are known. Two overlapping gene arrangements occupy different reading frames of the same region. Intergenic regions include most of the 25 promoters; transcripts are usually polycistronic. Translation of most of the open reading frames seems to be initiated independently, each from its own ribosomal binding and initiation site, although, a few cases of coupled translation have been reported. The most frequently used initiation codon is AUG but translation for a few open reading frames begins at GUG or UUG. The most common stop-codon is UGA followed by UAA and then UAG. Regulatory circuits are complex and largely dependent on two components of the central control operon. KorA and KorB are transcriptional repressors controlling at least seven operons. KorA and KorB act synergistically in several cases by recognizing and binding to conserved nucleotide sequences. Twelve KorB binding sites were found around the IncP alpha sequence and these are conserved in R751 (IncP beta) with respect to both sequence and location. Replication of IncP alpha plasmids requires oriV and the plasmid-encoded initiator protein TrfA in combination with the host-encoded replication machinery. Conjugative plasmid transfer depends on two separate regions occupying about half of the genome. The primary segregational stability system designated Par/Mrs consists of a putative site-specific recombinase, a possible partitioning apparatus and a post-segregational lethality mechanism, all encoded in two divergent operons. Proteins related to the products of F sop and P1 par partitioning genes are separately encoded in the central control operon.

Structure-function analysis of the human TFIIB-related factor II protein reveals an essential role for the C-terminal domain in RNA polymerase III transcription.

The transcription factors TFIIB, Brf1, and Brf2 share related N-terminal zinc ribbon and core domains. TFIIB bridges RNA polymerase II (Pol II) with the promoter-bound preinitiation complex, whereas Brf1 and Brf2 are involved, as part of activities also containing TBP and Bdp1 and referred to here as Brf1-TFIIIB and Brf2-TFIIIB, in the recruitment of Pol III. Brf1-TFIIIB recruits Pol III to type 1 and 2 promoters and Brf2-TFIIIB to type 3 promoters such as the human U6 promoter. Brf1 and Brf2 both have a C-terminal extension absent in TFIIB, but their C-terminal extensions are unrelated. In yeast Brf1, the C-terminal extension interacts with the TBP/TATA box complex and contributes to the recruitment of Bdp1. Here we have tested truncated Brf2, as well as Brf2/TFIIB chimeric proteins for U6 transcription and for assembly of U6 preinitiation complexes. Our results characterize functions of various human Brf2 domains and reveal that the C-terminal domain is required for efficient association of the protein with U6 promoter-bound TBP and SNAP(c), a type 3 promoter-specific transcription factor, and for efficient recruitment of Bdp1. This in turn suggests that the C-terminal extensions in Brf1 and Brf2 are crucial to specific recruitment of Pol III over Pol II.

Cross-species analysis of the replication complex of Old World arenaviruses reveals two nucleoprotein sites involved in L protein function.

Lassa virus (LASV) causing hemorrhagic Lassa fever in West Africa, Mopeia virus (MOPV) from East Africa, and lymphocytic choriomeningitis virus (LCMV) are the main representatives of the Old World arenaviruses. Little is known about how the components of the arenavirus replication machinery, i.e., the genome, nucleoprotein (NP), and L protein, interact. In addition, it is unknown whether these components can function across species boundaries. We established minireplicon systems for MOPV and LCMV in analogy to the existing LASV system and exchanged the components among the three systems. The functional and physical integrity of the resulting complexes was tested by reporter gene assay, Northern blotting, and coimmunoprecipitation studies. The minigenomes, NPs, and L proteins of LASV and MOPV could be exchanged without loss of function. LASV and MOPV L protein was also active in conjunction with LCMV NP, while the LCMV L protein required homologous NP for activity. Analysis of LASV/LCMV NP chimeras identified a single LCMV-specific NP residue (Ile-53) and the C terminus of NP (residues 340 to 558) as being essential for LCMV L protein function. The defect of LASV and MOPV NP in supporting transcriptional activity of LCMV L protein was not caused by a defect in physical NP-L protein interaction. In conclusion, components of the replication complex of Old World arenaviruses have the potential to functionally and physically interact across species boundaries. Residue 53 and the C-terminal domain of NP are important for function of L protein during genome replication and transcription.

Multicentric carpotarsal osteolysis is caused by mutations clustering in the amino-terminal transcriptional activation domain of MAFB.

Multicentric carpotarsal osteolysis (MCTO) is a rare skeletal dysplasia characterized by aggressive osteolysis, particularly affecting the carpal and tarsal bones, and is frequently associated with progressive renal failure. Using exome capture and next-generation sequencing in five unrelated simplex cases of MCTO, we identified previously unreported missense mutations clustering within a 51 base pair region of the single exon of MAFB, validated by Sanger sequencing. A further six unrelated simplex cases with MCTO were also heterozygous for previously unreported mutations within this same region, as were affected members of two families with autosomal-dominant MCTO. MAFB encodes a transcription factor that negatively regulates RANKL-induced osteoclastogenesis and is essential for normal renal development. Identification of this gene paves the way for development of novel therapeutic approaches for this crippling disease and provides insight into normal bone and kidney development.

Glucocorticoid-induced tumor necrosis factor receptor is a p21Cip1/WAF1 transcriptional target conferring resistance of keratinocytes to UV light-induced apoptosis.

Glucocorticoid-induced tumor necrosis factor receptor (GITR) is a member of the tumor necrosis factor receptor superfamily, is expressed in T lymphocytes, and exerts an anti-apoptotic function in these cells. We reported that GITR is also highly expressed in the skin, specifically in keratinocytes, and that it is under negative transcriptional control of p21(Cip1/WAF1), independently from the cell cycle. Although GITR expression is higher in p21-deficient keratinocytes and skin, it is down-modulated with differentiation and in response to UVB. The combined analysis of keratinocytes with increased GITR expression versus normal keratinocytes and skin of mice with a disruption of the GITR gene indicates that this protein protects keratinocytes from UVB-induced apoptosis both in vitro and in vivo.

Dissection of local and systemic transcriptional responses to phosphate starvation in Arabidopsis.

Phosphate is a crucial and often limiting nutrient for plant growth. To obtain inorganic phosphate (P(i) ), which is very insoluble, and is heterogeneously distributed in the soil, plants have evolved a complex network of morphological and biochemical processes. These processes are controlled by a regulatory system triggered by P(i) concentration, not only present in the medium (external P(i) ), but also inside plant cells (internal P(i) ). A 'split-root' assay was performed to mimic a heterogeneous environment, after which a transcriptomic analysis identified groups of genes either locally or systemically regulated by P(i) starvation at the transcriptional level. These groups revealed coordinated regulations for various functions associated with P(i) starvation (including P(i) uptake, P(i) recovery, lipid metabolism, and metal uptake), and distinct roles for members in gene families. Genetic tools and physiological analyses revealed that genes that are locally regulated appear to be modulated mostly by root development independently of the internal P(i) content. By contrast, internal P(i) was essential to promote the activation of systemic regulation. Reducing the flow of P(i) had no effect on the systemic response, suggesting that a secondary signal, independent of P(i) , could be involved in the response. Furthermore, our results display a direct role for the transcription factor PHR1, as genes systemically controlled by low P(i) have promoters enriched with P1BS motif (PHR1-binding sequences). These data detail various regulatory systems regarding P(i) starvation responses (systemic versus local, and internal versus external P(i) ), and provide tools to analyze and classify the effects of P(i) starvation on plant physiology.

Different transcriptional control of metabolism and extracellular matrix in visceral and subcutaneous fat of obese and rimonabant treated mice.

BACKGROUND: The visceral (VAT) and subcutaneous (SCAT) adipose tissues play different roles in physiology and obesity. The molecular mechanisms underlying their expansion in obesity and following body weight reduction are poorly defined. METHODOLOGY: C57Bl/6 mice fed a high fat diet (HFD) for 6 months developed low, medium, or high body weight as compared to normal chow fed mice. Mice from each groups were then treated with the cannabinoid receptor 1 antagonist rimonabant or vehicle for 24 days to normalize their body weight. Transcriptomic data for visceral and subcutaneous adipose tissues from each group of mice were obtained and analyzed to identify: i) genes regulated by HFD irrespective of body weight, ii) genes whose expression correlated with body weight, iii) the biological processes activated in each tissue using gene set enrichment analysis (GSEA), iv) the transcriptional programs affected by rimonabant. PRINCIPAL FINDINGS: In VAT, "metabolic" genes encoding enzymes for lipid and steroid biosynthesis and glucose catabolism were down-regulated irrespective of body weight whereas "structure" genes controlling cell architecture and tissue remodeling had expression levels correlated with body weight. In SCAT, the identified "metabolic" and "structure" genes were mostly different from those identified in VAT and were regulated irrespective of body weight. GSEA indicated active adipogenesis in both tissues but a more prominent involvement of tissue stroma in VAT than in SCAT. Rimonabant treatment normalized most gene expression but further reduced oxidative phosphorylation gene expression in SCAT but not in VAT. CONCLUSION: VAT and SCAT show strikingly different gene expression programs in response to high fat diet and rimonabant treatment. Our results may lead to identification of therapeutic targets acting on specific fat depots to control obesity.

An operon of three transcriptional regulators controls horizontal gene transfer of the integrative and conjugative element ICEclc in Pseudomonas knackmussii B13.

The integrative and conjugative element ICEclc is a mobile genetic element in Pseudomonas knackmussii B13, and an experimental model for a widely distributed group of elements in Proteobacteria. ICEclc is transferred from specialized transfer competent cells, which arise at a frequency of 3-5% in a population at stationary phase. Very little is known about the different factors that control the transfer frequency of this ICE family. Here we report the discovery of a three-gene operon encoded by ICEclc, which exerts global control on transfer initiation. The operon consists of three consecutive regulatory genes, encoding a TetR-type repressor MfsR, a MarR-type regulator and a LysR-type activator TciR. We show that MfsR autoregulates expression of the operon, whereas TciR is a global activator of ICEclc gene expression, but no clear role was yet found for MarR. Deletion of mfsR increases expression of tciR and marR, causing the proportion of transfer competent cells to reach almost 100% and transfer frequencies to approach 1 per donor. mfsR deletion also caused a two orders of magnitude loss in population viability, individual cell growth arrest and loss of ICEclc. This indicates that autoregulation is an important feature maintaining ICE transfer but avoiding fitness loss. Bioinformatic analysis showed that the mfsR-marR-tciR operon is unique for ICEclc and a few highly related ICE, whereas tciR orthologues occur more widely in a large variety of suspected ICE among Proteobacteria.

Global transcriptional programs in peripheral nerve endoneurium and DRG are resistant to the onset of type 1 diabetic neuropathy in Ins2 mice.

While the morphological and electrophysiological changes underlying diabetic peripheral neuropathy (DPN) are relatively well described, the involved molecular mechanisms remain poorly understood. In this study, we investigated whether phenotypic changes associated with early DPN are correlated with transcriptional alterations in the neuronal (dorsal root ganglia [DRG]) or the glial (endoneurium) compartments of the peripheral nerve. We used Ins2(Akita/+) mice to study transcriptional changes underlying the onset of DPN in type 1 diabetes mellitus (DM). Weight, blood glucose and motor nerve conduction velocity (MNCV) were measured in Ins2(Akita/+) and control mice during the first three months of life in order to determine the onset of DPN. Based on this phenotypic characterization, we performed gene expression profiling using sciatic nerve endoneurium and DRG isolated from pre-symptomatic and early symptomatic Ins2(Akita/+) mice and sex-matched littermate controls. Our phenotypic analysis of Ins2(Akita/+) mice revealed that DPN, as measured by reduced MNCV, is detectable in affected animals already one week after the onset of hyperglycemia. Surprisingly, the onset of DPN was not associated with any major persistent changes in gene expression profiles in either sciatic nerve endoneurium or DRG. Our data thus demonstrated that the transcriptional programs in both endoneurial and neuronal compartments of the peripheral nerve are relatively resistant to the onset of hyperglycemia and hypoinsulinemia suggesting that either minor transcriptional alterations or changes on the proteomic level are responsible for the functional deficits associated with the onset of DPN in type 1 DM.

Genome-wide transcriptional responses to shade: Linking shade avoidance and auxin

Plants have the ability to use the composition of incident light as a cue to adapt development and growth to their environment. Arabidopsis thaliana as well as many crops are best adapted to sunny habitats. When subjected to shade, these plants exhibit a variety of physiological responses collectively called shade avoidance syndrome (SAS). It includes increased growth of hypocotyl and petioles, decreased growth rate of cotyledons and reduced branching and crop yield. These responses are mainly mediated by phytochrome photoreceptors, which exist either in an active, far-red light (FR) absorbing or an inactive, red light (R) absorbing isoform. In direct sunlight, the R to FR light (R/FR) ratio is high and converts the phytochromes into their physiologically active state. The phytochromes interact with downstream transcription factors such as PHYTOCHROME INTERACTING FACTOR (PIF), which are subsequently degraded. Light filtered through a canopy is strongly depleted in R, which result in a low R/FR ratio and renders the phytochromes inactive. Protein levels of downstream transcription factors are stabilized, which initiates the expression of shade-induced genes such as HFR1, PIL1 or ATHB-2. In my thesis, I investigated transcriptional responses mediated by the SAS in whole Arabidopsis seedlings. Using microarray and chromatin immunoprecipitation data, we identified genome-wide PIF4 and PIF5 dependent shade regulated gene as well as putative direct target genes of PIF5. This revealed evidence for a direct regulatory link between phytochrome signaling and the growth promoting phytohormone auxin (IAA) at the level of biosynthesis, transport and signaling. Subsequently, it was shown, that free-IAA levels are upregulated in response to shade. It is assumed that shade-induced auxin production takes predominantly place in cotyledons of seedlings. This implies, that IAA is subsequently transported basipetally to the hypocotyl and enhances elongation growth. The importance of auxin transport for growth responses has been established by chemical and genetic approaches. To gain a better understanding of spatio-temporal transcriptional regulation of shade-induce auxin, I generated in a second project, an organ specific high throughput data focusing on cotyledon and hypocotyl of young Arabidopsis seedlings. Interestingly, both organs show an opposite growth regulation by shade. I first investigated the spatio-transcriptional regulation of auxin re- sponsive gene, in order to determine how broad gene expression pattern can be explained by the hypothesized movement of auxin from cotyledons to hypocotyls in shade. The analysis suggests, that several genes are indeed regulated according to our prediction and others are regulated in a more complex manner. In addition, analysis of gene families of auxin biosynthetic and transport components, lead to the identification of essential family members for shade-induced growth re- sponses, which were subsequently experimentally confirmed. Finally, the analysis of expression pattern identified several candidate genes, which possibly explain aspects of the opposite growth response of the different organs.

Research resource: the dynamic transcriptional profile of sertoli cells during the progression of spermatogenesis.

Sertoli cells (SCs), the only somatic cells within seminiferous tubules, associate intimately with developing germ cells. They not only provide physical and nutritional support but also secrete factors essential to the complex developmental processes of germ cell proliferation and differentiation. The SC transcriptome must therefore adapt rapidly during the different stages of spermatogenesis. We report comprehensive genome-wide expression profiles of pure populations of SCs isolated at 5 distinct stages of the first wave of mouse spermatogenesis, using RNA sequencing technology. We were able to reconstruct about 13 901 high-confidence, nonredundant coding and noncoding transcripts, characterized by complex alternative splicing patterns with more than 45% comprising novel isoforms of known genes. Interestingly, roughly one-fifth (2939) of these genes exhibited a dynamic expression profile reflecting the evolving role of SCs during the progression of spermatogenesis, with stage-specific expression of genes involved in biological processes such as cell cycle regulation, metabolism and energy production, retinoic acid synthesis, and blood-testis barrier biogenesis. Finally, regulatory network analysis identified the transcription factors endothelial PAS domain-containing protein 1 (EPAS1/Hif2α), aryl hydrocarbon receptor nuclear translocator (ARNT/Hif1β), and signal transducer and activator of transcription 1 (STAT1) as potential master regulators driving the SC transcriptional program. Our results highlight the plastic transcriptional landscape of SCs during the progression of spermatogenesis and provide valuable resources to better understand SC function and spermatogenesis and its related disorders, such as male infertility.

Caveolin-1-mediated post-transcriptional regulation of inducible nitric oxide synthase in human colon carcinoma cells.

Reactive oxygen species are now widely recognized as important players contributing both to cell homeostasis and the development of disease. In this respect nitric oxide (NO) is no exception. The discussion here will center on regulation of the inducible form of nitric oxide synthase (iNOS) for two reasons. First, only iNOS produces micromolar NO concentrations, amounts that are high by comparison with the picomolar to nanomolar concentrations resulting from Ca2(+)-controlled NO production by endothelial eNOS or neuronal nNOS. Second, iNOS is not constitutively expressed in cells and regulation of this isoenzyme, in contrast to endothelial eNOS or neuronal nNOS, is widely considered to occur at the transcriptional level only. In particular, we were interested in the possibility that caveolin-1, a protein that functions as a tumor suppressor in colon carcinoma cells (Bender et al., 2002; this issue), might regulate iNOS activity. Our results provide evidence for the existence of a post-transcriptional mechanism controlling iNOS protein levels that involves caveolin-1-dependent sequestration of iNOS within a detergent-insoluble compartment. Interestingly, despite the high degree of conservation of the caveolin-1 scaffolding domain binding motif within all NOS enzymes, the interaction detected between caveolin-1 and iNOS in vitro is crucially dependent on presence of a caveolin-1 sequence element immediately adjacent to the scaffolding domain. A model is presented summarizing the salient aspects of these results. These observations are important in the context of tumor biology, since down-regulation of caveolin-1 is predicted to promote uncontrolled iNOS activity, genotoxic damage and thereby facilitate tumor development in humans.