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.

Death and resurrection of the human IRGM gene

Immunity-related GTPases (IRG) play an important role in defense against intracellular pathogens. One member of this gene family in humans, IRGM, has been recently implicated as a risk factor for Crohn's disease. We analyzed the detailed structure of this gene family among primates and showed that most of the IRG gene cluster was deleted early in primate evolution, after the divergence of the anthropoids from prosimians ( about 50 million years ago). Comparative sequence analysis of New World and Old World monkey species shows that the single-copy IRGM gene became pseudogenized as a result of an Alu retrotransposition event in the anthropoid common ancestor that disrupted the open reading frame (ORF). We find that the ORF was reestablished as a part of a polymorphic stop codon in the common ancestor of humans and great apes. Expression analysis suggests that this change occurred in conjunction with the insertion of an endogenous retrovirus, which altered the transcription initiation, splicing, and expression profile of IRGM. These data argue that the gene became pseudogenized and was then resurrected through a series of complex structural events and suggest remarkable functional plasticity where alleles experience diverse evolutionary pressures over time. Such dynamism in structure and evolution may be critical for a gene family locked in an arms race with an ever-changing repertoire of intracellular parasites.

Pro-inflammatory gene expression and neurotoxic effects of activated microglia are attenuated by absence of CCAAT/enhancer binding protein beta

Background. Microglia and astrocytes respond to homeostatic disturbances with profound changes of gene expression. This response, known as glial activation or neuroinflammation, can be detrimental to the surrounding tissue. The transcription factor CCAAT/enhancer binding protein ß (C/EBPß) is an important regulator of gene expression in inflammation but little is known about its involvement in glial activation. To explore the functional role of C/EBPß in glial activation we have analyzed pro-inflammatory gene expression and neurotoxicity in murine wild type and C/EBPß-null glial cultures. Methods. Due to fertility and mortality problems associated with the C/EBPß-null genotype we developed a protocol to prepare mixed glial cultures from cerebral cortex of a single mouse embryo with high yield. Wild-type and C/EBPß-null glial cultures were compared in terms of total cell density by Hoechst-33258 staining; microglial content by CD11b immunocytochemistry; astroglial content by GFAP western blot; gene expression by quantitative real-time PCR, western blot, immunocytochemistry and Griess reaction; and microglial neurotoxicity by estimating MAP2 content in neuronal/microglial cocultures. C/EBPß DNA binding activity was evaluated by electrophoretic mobility shift assay and quantitative chromatin immunoprecipitation. Results. C/EBPß mRNA and protein levels, as well as DNA binding, were increased in glial cultures by treatment with lipopolysaccharide (LPS) or LPS + interferon ¿ (IFN¿). Quantitative chromatin immunoprecipitation showed binding of C/EBPß to pro-inflammatory gene promoters in glial activation in a stimulus- and gene-dependent manner. In agreement with these results, LPS and LPS+IFN¿ induced different transcriptional patterns between pro-inflammatory cytokines and NO synthase-2 genes. Furthermore, the expressions of IL-1ß and NO synthase-2, and consequent NO production, were reduced in the absence of C/EBPß. In addition, neurotoxicity elicited by LPS+IFN¿-treated microglia co-cultured with neurons was completely abolished by the absence of C/EBPß in microglia.

Nutritional and insulin regulation of fatty acid synthetase and leptin gene expression through ADD1/SREBP1.

The ability to regulate specific genes of energy metabolism in response to fasting and feeding is an important adaptation allowing survival of intermittent food supplies. However, little is known about transcription factors involved in such responses in higher organisms. We show here that gene expression in adipose tissue for adipocyte determination differentiation dependent factor (ADD) 1/sterol regulatory element binding protein (SREBP) 1, a basic-helix-loop-helix protein that has a dual DNA-binding specificity, is reduced dramatically upon fasting and elevated upon refeeding; this parallels closely the regulation of two adipose cell genes that are crucial in energy homeostasis, fatty acid synthetase (FAS) and leptin. This elevation of ADD1/SREBP1, leptin, and FAS that is induced by feeding in vivo is mimicked by exposure of cultured adipocytes to insulin, the classic hormone of the fed state. We also show that the promoters for both leptin and FAS are transactivated by ADD1/SREBP1. A mutation in the basic domain of ADD1/SREBP1 that allows E-box binding but destroys sterol regulatory element-1 binding prevents leptin gene transactivation but has no effect on the increase in FAS promoter function. Molecular dissection of the FAS promoter shows that most if not all of this action of ADD1/SREBP1 is through an E-box motif at -64 to -59, contained with a sequence identified previously as the major insulin response element of this gene. These results indicate that ADD1/SREBP1 is a key transcription factor linking changes in nutritional status and insulin levels to the expression of certain genes that regulate systemic energy metabolism.

dKDM5/LID regulates H3K4me3 dynamics at the transcription-start site (TSS) of actively transcribed developmental genes

H3K4me3 is a histone modification that accumulates at the transcription-start site (TSS) of active genes and is known to be important for transcription activation. The way in which H3K4me3 is regulated at TSS and the actual molecular basis of its contribution to transcription remain largely unanswered. To address these questions, we have analyzed the contribution of dKDM5/LID, the main H3K4me3 demethylase in Drosophila, to the regulation of the pattern of H3K4me3. ChIP-seq results show that, at developmental genes, dKDM5/LID localizes at TSS and regulates H3K4me3. dKDM5/LID target genes are highly transcribed and enriched in active RNApol II and H3K36me3, suggesting a positive contribution to transcription. Expression-profiling show that, though weakly, dKDM5/LID target genes are significantly downregulated upon dKDM5/LID depletion. Furthermore, dKDM5/LID depletion results in decreased RNApol II occupancy, particularly by the promoter-proximal Pol lloser5 form. Our results also show that ASH2, an evolutionarily conserved factor that locates at TSS and is required for H3K4me3, binds and positively regulates dKDM5/LID target genes. However, dKDM5/LID and ASH2 do not bind simultaneously and recognize different chromatin states, enriched in H3K4me3 and not, respectively. These results indicate that, at developmental genes, dKDM5/LID and ASH2 coordinately regulate H3K4me3 at TSS and that this dynamic regulation contributes to transcription.

Thermal evolution of gene expression profiles in Drosophila subobscura

Background: Despite its pervasiveness, the genetic basis of adaptation resulting in variation directly or indirectly related to temperature (climatic) gradients is poorly understood. By using 3-fold replicated laboratory thermal stocks covering much of the physiologically tolerable temperature range for the temperate (i.e., cold tolerant) species Drosophila subobscura we have assessed whole-genome transcriptional responses after three years of thermal adaptation, when the populations had already diverged for inversion frequencies, pre-adult life history components, and morphological traits. Total mRNA from each population was compared to a reference pool mRNA in a standard, highly replicated two-colour competitive hybridization experiment using cDNA microarrays.Results: A total of 306 (6.6%) cDNA clones were identified as 'differentially expressed' (following a false discovery rate correction) after contrasting the two furthest apart thermal selection regimes (i.e., 13°C vs . 22°C), also including four previously reported candidate genes for thermotolerance in Drosophila (Hsp26, Hsp68, Fst, and Treh). On the other hand, correlated patterns of gene expression were similar in cold- and warm-adapted populations. Analysis of functional categories defined by the Gene Ontology project point to an overrepresentation of genes involved in carbohydrate metabolism, nucleic acids metabolism and regulation of transcription among other categories. Although the location of differently expressed genes was approximately at random with respect to chromosomes, a physical mapping of 88 probes to the polytene chromosomes of D. subobscura has shown that a larger than expected number mapped inside inverted chromosomal segments.Conclusion: Our data suggest that a sizeable number of genes appear to be involved in thermal adaptation in Drosophila, with a substantial fraction implicated in metabolism. This apparently illustrates the formidable challenge to understanding the adaptive evolution of complex trait variation. Furthermore, some clustering of genes within inverted chromosomal sections was detected. Disentangling the effects of inversions will be obviously required in any future approach if we want to identify the relevant candidate genes.

IFN-gamma-dependent transcription of MHC class II IA is impaired in macrophages from aged mice

To determine the effect of aging on IFN-gamma-induced MHC class II antigen expression, we produced bone marrow¿derived macrophages in vitro. In these conditions, we analyzed the effect of aging on the genomic expression of macrophages without the influence of other cell types that may be affected by aging. Although macrophages from young and aged mice showed an identical degree of differentiation, after incubation with IFN-gamma, the expression at the cell surface of the IA complex and the levels of IAbeta protein and mRNA were lower in aged macrophages. Moreover, the transcription of the IAbeta gene was impaired in aged macrophages. The amount of transcription factors that bound to the W and X, but not to the Y, boxes of the IAbeta promoter gene was lower in aged macrophages. Similar levels of CIITA mRNA were found after IFN-gamma treatment of both young and aged macrophages. This shows that neither the initial cascade that starts after the interaction of IFN-gamma with the receptor nor the second signals involved in the expression of CIITA are impaired in aged macrophages. These data indicate that aging is associated with low levels of MHC class II gene induction by IFN-gamma because of impaired transcription.

Functional interaction between peroxisome proliferator-activated receptors-alpha and Mef-2C on human carnitine palmitoyltransferase 1beta (CPT1beta) gene activation

Muscle-type carnitine palmitoyltransferase 1 (CPT1β) is considered to be the gene that controls fatty acid mitochondrial β-oxidation. A functional peroxisome proliferator-activated receptor (PPAR) responsive element (PPRE) and a myocite-specific (MEF2) site that binds MEF2A and MEF2C in the promoter of this gene had been previously identified. We investigated the roles of the PPRE and the MEF2 binding sites and the potential interaction between PPARα and MEF2C regulating the CPT1β gene promoter. Mutation analysis indicated that the MEF2 site contributed to the activation of the CPT1β promoter by PPAR in C2C12 cells. The reporter construct containing the PPRE and the MEF2C site was synergistically activated by co-expression of PPAR, retinoid X receptor (RXR) and MEF2C in non-muscle cells. Moreover, protein-binding assays demonstrated that MEF2C and PPAR specifically bound to one another in vitro. Also for the synergistic activation of the CPT1β gene promoter by MEF2C and PPARα-RXRα, a precise arrangement of its binding sites was essential.

Hog1 bypasses stress-mediated down-regulation of transcription by RNA polymerase II redistribution and chromatin remodeling

Cells are subjected to dramatic changes of gene expression upon environmental changes. Stresscauses a general down-regulation of gene expression together with the induction of a set of stress-responsivegenes. The p38-related stress-activated protein kinase Hog1 is an important regulator of transcription uponosmostress in yeast. Genome-wide localization studies of RNA polymerase II (RNA Pol II) and Hog1 showed that stress induced major changes in RNA Pol II localization, with a shift toward stress-responsive genes relative to housekeeping genes. RNA Pol II relocalization required Hog1, which was also localized to stress-responsive loci. In addition to RNA Pol II-bound genes, Hog1 also localized to RNA polymerase III-bound genes, pointing to a wider role for Hog1 in transcriptional control than initially expected. Interestingly, an increasing association of Hog1 with stressresponsive genes was strongly correlated with chromatin remodeling and increased gene expression. Remarkably, MNase-Seq analysis showed that although chromatin structure was not significantly altered at a genome-wide level in response to stress, there was pronounced chromatin remodeling for those genes that displayed Hog1 association. Hog1 serves to bypass the general down-regulation of gene expression that occurs in response to osmostress, and does so both by targeting RNA Pol II machinery and by inducing chromatin remodeling at stressresponsive loci.

IFN-gamma-dependent transcription of MHC class II IA is impaired in macrophages from aged mice

To determine the effect of aging on IFN-gamma-induced MHC class II antigen expression, we produced bone marrow¿derived macrophages in vitro. In these conditions, we analyzed the effect of aging on the genomic expression of macrophages without the influence of other cell types that may be affected by aging. Although macrophages from young and aged mice showed an identical degree of differentiation, after incubation with IFN-gamma, the expression at the cell surface of the IA complex and the levels of IAbeta protein and mRNA were lower in aged macrophages. Moreover, the transcription of the IAbeta gene was impaired in aged macrophages. The amount of transcription factors that bound to the W and X, but not to the Y, boxes of the IAbeta promoter gene was lower in aged macrophages. Similar levels of CIITA mRNA were found after IFN-gamma treatment of both young and aged macrophages. This shows that neither the initial cascade that starts after the interaction of IFN-gamma with the receptor nor the second signals involved in the expression of CIITA are impaired in aged macrophages. These data indicate that aging is associated with low levels of MHC class II gene induction by IFN-gamma because of impaired transcription.

dKDM5/LID regulates H3K4me3 dynamics at the transcription-start site (TSS) of actively transcribed developmental genes

H3K4me3 is a histone modification that accumulates at the transcription-start site (TSS) of active genes and is known to be important for transcription activation. The way in which H3K4me3 is regulated at TSS and the actual molecular basis of its contribution to transcription remain largely unanswered. To address these questions, we have analyzed the contribution of dKDM5/LID, the main H3K4me3 demethylase in Drosophila, to the regulation of the pattern of H3K4me3. ChIP-seq results show that, at developmental genes, dKDM5/LID localizes at TSS and regulates H3K4me3. dKDM5/LID target genes are highly transcribed and enriched in active RNApol II and H3K36me3, suggesting a positive contribution to transcription. Expression-profiling show that, though weakly, dKDM5/LID target genes are significantly downregulated upon dKDM5/LID depletion. Furthermore, dKDM5/LID depletion results in decreased RNApol II occupancy, particularly by the promoter-proximal Pol lloser5 form. Our results also show that ASH2, an evolutionarily conserved factor that locates at TSS and is required for H3K4me3, binds and positively regulates dKDM5/LID target genes. However, dKDM5/LID and ASH2 do not bind simultaneously and recognize different chromatin states, enriched in H3K4me3 and not, respectively. These results indicate that, at developmental genes, dKDM5/LID and ASH2 coordinately regulate H3K4me3 at TSS and that this dynamic regulation contributes to transcription.

Heat shock response in yeast involves changes in both transcription rates and mRNA stabilities

We have analyzed the heat stress response in the yeast Saccharomyces cerevisiae by determining mRNA levels and transcription rates for the whole transcriptome after a shift from 25uC to 37uC. Using an established mathematical algorithm, theoretical mRNA decay rates have also been calculated from the experimental data. We have verified the mathematical predictions for selected genes by determining their mRNA decay rates at different times during heat stress response using the regulatable tetO promoter. This study indicates that the yeast response to heat shock is not only due to changes in transcription rates, but also to changes in the mRNA stabilities. mRNA stability is affected in 62% of the yeast genes and it is particularly important in shaping the mRNA profile of the genes belonging to the environmental stress response. In most cases, changes in transcription rates and mRNA stabilities are homodirectional for both parameters, although some interesting cases of antagonist behavior are found. The statistical analysis of gene targets and sequence motifs within the clusters of genes with similar behaviors shows that both transcriptional and post-transcriptional regulons apparently contribute to the general heat stress response by means of transcriptional factors and RNA binding proteins.

A functional variant in the Stearoyl-CoA desaturase gene promoter enhances fatty acid desaturation in pork

There is growing public concern about reducing saturated fat intake. Stearoyl-CoA desaturase (SCD) is the lipogenic enzyme responsible for the biosynthesis of oleic acid (18:1) by desaturating stearic acid (18:0). Here we describe a total of 18 mutations in the promoter and 3′ non-coding region of the pig SCD gene and provide evidence that allele T at AY487830:g.2228T>C in the promoter region enhances fat desaturation (the ratio 18:1/18:0 in muscle increases from 3.78 to 4.43 in opposite homozygotes) without affecting fat content (18:0+18:1, intramuscular fat content, and backfat thickness). No mutations that could affect the functionality of the protein were found in the coding region. First, we proved in a purebred Duroc line that the C-T-A haplotype of the 3 single nucleotide polymorphisms (SNPs) (g.2108C>T; g.2228T>C; g.2281A>G) of the promoter region was additively associated to enhanced 18:1/18:0 both in muscle and subcutaneous fat, but not in liver. We show that this association was consistent over a 10-year period of overlapping generations and, in line with these results, that the C-T-A haplotype displayed greater SCD mRNA expression in muscle. The effect of this haplotype was validated both internally, by comparing opposite homozygote siblings, and externally, by using experimental Duroc-based crossbreds. Second, the g.2281A>G and the g.2108C>T SNPs were excluded as causative mutations using new and previously published data, restricting the causality to g.2228T>C SNP, the last source of genetic variation within the haplotype. This mutation is positioned in the core sequence of several putative transcription factor binding sites, so that there are several plausible mechanisms by which allele T enhances 18:1/18:0 and, consequently, the proportion of monounsaturated to saturated fat.

The Saccharomyces cerevisiae Hot1p regulated gene YHR087W (HGI1) has a role in translation upon high glucose concentration stress.

Background While growing in natural environments yeasts can be affected by osmotic stress provoked by high glucose concentrations. The response to this adverse condition requires the HOG pathway and involves transcriptional and posttranscriptional mechanisms initiated by the phosphorylation of this protein, its translocation to the nucleus and activation of transcription factors. One of the genes induced to respond to this injury is YHR087W. It encodes for a protein structurally similar to the N-terminal region of human SBDS whose expression is also induced under other forms of stress and whose deletion determines growth defects at high glucose concentrations. Results In this work we show that YHR087W expression is regulated by several transcription factors depending on the particular stress condition, and Hot1p is particularly relevant for the induction at high glucose concentrations. In this situation, Hot1p, together to Sko1p, binds to YHR087W promoter in a Hog1p-dependent manner. Several evidences obtained indicate Yhr087wp"s role in translation. Firstly, and according to TAP purification experiments, it interacts with proteins involved in translation initiation. Besides, its deletion mutant shows growth defects in the presence of translation inhibitors and displays a slightly slower translation recovery after applying high glucose stress than the wild type strain. Analyses of the association of mRNAs to polysome fractions reveals a lower translation in the mutant strain of the mRNAs corresponding to genes GPD1, HSP78 and HSP104. Conclusions The data demonstrates that expression of Yhr087wp under high glucose concentration is controlled by Hot1p and Sko1p transcription factors, which bind to its promoter. Yhr087wp has a role in translation, maybe in the control of the synthesis of several stress response proteins, which could explain the lower levels of some of these proteins found in previous proteomic analyses and the growth defects of the deletion strain. Keywords: Saccharomyces cerevisiae; High glucose osmotic stress; Gene YHR087W; Gene expression; Translation; Hot1p; Hog1p; Polysomes