Identification of Conserved Splicing Motifs in Mutually Exclusive Exons of 15 Insect Species

Background: During alternative splicing, the inclusion of an exon in the final mRNA molecule is determined by nuclear proteins that bind cis-regulatory sequences in a target pre-mRNA molecule. A recent study suggested that the regulatory codes of individual RNA-binding proteins may be nearly immutable between very diverse species such as mammals and insects. The model system Drosophila melanogaster therefore presents an excellent opportunity for the study of alternative splicing due to the availability of quality EST annotations in FlyBase. Methods: In this paper, we describe an in silico analysis pipeline to extract putative exonic splicing regulatory sequences from a multiple alignment of 15 species of insects. Our method, ESTs-to-ESRs (E2E), uses graph analysis of EST splicing graphs to identify mutually exclusive (ME) exons and combines phylogenetic measures, a sliding window approach along the multiple alignment and the Welch’s t statistic to extract conserved ESR motifs. Results: The most frequent 100% conserved word of length 5 bp in different insect exons was “ATGGA”. We identified 799 statistically significant “spike” hexamers, 218 motifs with either a left or right FDR corrected spike magnitude p-value < 0.05 and 83 with both left and right uncorrected p < 0.01. 11 genes were identified with highly significant motifs in one ME exon but not in the other, suggesting regulation of ME exon splicing through these highly conserved hexamers. The majority of these genes have been shown to have regulated spatiotemporal expression. 10 elements were found to match three mammalian splicing regulator databases. A putative ESR motif, GATGCAG, was identified in the ME-13b but not in the ME-13a of Drosophila N-Cadherin, a gene that has been shown to have a distinct spatiotemporal expression pattern of spliced isoforms in a recent study. Conclusions: Analysis of phylogenetic relationships and variability of sequence conservation as implemented in the E2E spikes method may lead to improved identification of ESRs. We found that approximately half of the putative ESRs in common between insects and mammals have a high statistical support (p < 0.01). Several Drosophila genes with spatiotemporal expression patterns were identified to contain putative ESRs located in one exon of the ME exon pairs but not in the other.

Functional characterization of the sciarid BhC4-1 core promoter in transgenic Drosophila

Background: Core promoters are cis-regulatory modules to which bind the basal transcriptional machinery and which participate in the regulation of transcription initiation. Although core promoters have not been extensively investigated through functional assays in a chromosomal context, the available data suggested that the response of a given core promoter might vary depending on the promoter context. Previous studies suggest that a (-57/+40) fragment constitutes the core promoter of the BhC4-1 gene which is located in DNA puff C4 of the sciarid fly Bradysia hygida. Here we tested this (-57/+40) fragment in distinct regulatory contexts in order to verify if promoter context affects its core promoter activity. Results: Consistent with the activity of a core promoter, we showed that in the absence of upstream regulatory sequences the (-57/+40) fragment drives low levels of reporter gene mRNA expression throughout development in transgenic Drosophila. By assaying the (-57/+40) fragment in two distinct regulatory contexts, either downstream of the previously characterized Fbp1 enhancer or downstream of the UAS element, we showed that the BhC4-1 core promoter drives regulated transcription in both the germline and in various tissues throughout development. Furthermore, the use of the BhC4-1 core promoter in a UAS construct significantly reduced salivary gland ectopic expression in third instar larvae, which was previously described to occur in the context of the GAL4/UAS system. Conclusions: Our results from functional analysis in transgenic Drosophila show that the BhC4-1 core promoter drives gene expression regardless of the promoter context that was assayed. New insights into the functioning of the GAL4/UAS system in Drosophila were obtained, indicating that the presence of the SV40 sequence in the 3' UTR of a UAS construct does not preclude expression in the germline. Furthermore, our analysis indicated that ectopic salivary gland expression in the GAL4/UAS system does not depend only on sequences present in the GAL4 construct, but can also be affected by the core promoter sequences in the UAS construct. In this context, we propose that the sciarid BhC4-1 core promoter constitutes a valuable core promoter which can be employed in functional assays in insects.

Site-directed mutagenesis of the ATM promoter: Consequences for response to proliferation and ionizing radiation

Although ATM, the protein defective in ataxia-telangiectasia (A-T), is activated primarily by radiation, there is also evidence that expression of the protein can be regulated by both radiation and growth factors. Computer analysis of the ATM promoter proximal 700-bp sequence reveals a number of potentially important cis-regulatory sequences. Using nucleotide substitutions to delete putative functional elements in the promoter of ATM, we examined the importance of some of these sites for both the basal and the radiation-induced activity of the promoter. In lymphoblastoid cells, most of the mutations in transcription factor consensus sequences [Sp1(1), Sp1(2), Cre, Ets, Xre, gammaIre(2), a modified AP1 site (Fse), and GCF] reduced basal activity to various extents, whereas others [gammaIre(1), NF1, Myb] left basal activity unaffected. In human skin fibroblasts, results were generally the same, but the basal activity varied up to 8-fold in these and other cell lines. Radiation activated the promoter approximately 2.5-fold in serum-starved lymphoblastoid cells, reaching a maximum by 3 hr, and all mutated elements equally blocked this activation. Reduction in Sp1 and AP1 DNA binding activity by serum starvation was rapidly reversed by exposure of cells to radiation. This reduction was not evident in A-T cells, and the response to radiation was less marked. Data provided for interaction between ATM and Sp1 by protein binding and co-immunoprecipitation could explain the altered regulation of Sp1 in A-T cells. The data described here provide additional evidence that basal and radiation-induced regulation of the ATM promoter is under multifactorial control. (C) 2003 Wiley-Liss, Inc.

Identification and analysis of functional elements in 1% of the human genome by the ENCODE pilot project.

We report the generation and analysis of functional data from multiple, diverse experiments performed on a targeted 1% of the human genome as part of the pilot phase of the ENCODE Project. These data have been further integrated and augmented by a number of evolutionary and computational analyses. Together, our results advance the collective knowledge about human genome function in several major areas. First, our studies provide convincing evidence that the genome is pervasively transcribed, such that the majority of its bases can be found in primary transcripts, including non-protein-coding transcripts, and those that extensively overlap one another. Second, systematic examination of transcriptional regulation has yielded new understanding about transcription start sites, including their relationship to specific regulatory sequences and features of chromatin accessibility and histone modification. Third, a more sophisticated view of chromatin structure has emerged, including its inter-relationship with DNA replication and transcriptional regulation. Finally, integration of these new sources of information, in particular with respect to mammalian evolution based on inter- and intra-species sequence comparisons, has yielded new mechanistic and evolutionary insights concerning the functional landscape of the human genome. Together, these studies are defining a path for pursuit of a more comprehensive characterization of human genome function.

Cloning and expression of rat pancreatic beta-cell malonyl-CoA decarboxylase.

To gain insight into the function and regulation of malonyl-CoA decarboxylase (MCD) we have cloned rat MCD cDNA from a differentiated insulin-secreting pancreatic beta-cell-line cDNA library. The full-length cDNA sequence shows 69% identity with the cDNA cloned previously from the goose uropygial gland, and predicts a 492 amino acid protein of 54.7 kDa. The open reading frame contains an N-terminal mitochondrial targeting sequence and the C-terminal part of the enzyme ends with a peroxisomal (Ser-Lys-Leu) targeting motif. Since the sequence does not reveal hydrophobic domains, MCD is most likely expressed in the mitochondrial matrix and inside the peroxisomes. A second methionine residue, located 3' of the mitochondrial presequence, might be the first amino acid of a putative cytosolic MCD, since the nucleotide sequence around it fits fairly well with a consensus Kozak site for translation initiation. However, primer extension detects the presence of only one transcript initiating upstream of the first ATG, indicating that the major, if not exclusive, transcript expressed in the pancreatic beta-cell encodes MCD with its mitochondrial presequence. The sequence also shows multiple possible sites of phosphorylation by casein kinase II and protein kinase C. mRNA tissue-distribution analysis indicates a transcript of 2.2 kb, and that the MCD gene is expressed over a wide range of rat tissues. The distribution of the enzyme shows a broad range of activities from very low in the brain to elevated in the liver and heart. The results provide the foundations for further studies of the role of MCD in lipid metabolism and metabolic signalling in various tissues.

The ENCODE (ENCyclopedia Of DNA Elements) Project.

The ENCyclopedia Of DNA Elements (ENCODE) Project aims to identify all functional elements in the human genome sequence. The pilot phase of the Project is focused on a specified 30 megabases (approximately 1%) of the human genome sequence and is organized as an international consortium of computational and laboratory-based scientists working to develop and apply high-throughput approaches for detecting all sequence elements that confer biological function. The results of this pilot phase will guide future efforts to analyze the entire human genome.

Regulation of Placental Leptin Expression by Cyclic Adenosine 5 '-Monophosphate Involves Cross Talk between Protein Kinase A and Mitogen-Activated Protein Kinase Signaling Pathways

Leptin, a 16-kDa protein mainly produced by adipose tissue, has been involved in the control of energy balance through its hypothalamic receptor. However, pleiotropic effects of leptin have been identified in reproduction and pregnancy, particularly in placenta, where it was found to be expressed. In the current study, we examined the effect of cAMP in the regulation of leptin expression in trophoblastic cells. We found that dibutyryl cAMP [(Bu)(2)cAMP], a cAMP analog, showed an inducing effect on endogenous leptin expression in BeWo and JEG-3 cell lines when analyzed by Western blot analysis and quantitative RT-PCR. Maximal effect was achieved at 100 microM. Leptin promoter activity was also stimulated, evaluated by transient transfection with a reporter plasmid construction. Similar results were obtained with human term placental explants, thus indicating physiological relevance. Because cAMP usually exerts its actions through activation of protein kinase A (PKA) signaling, this pathway was analyzed. We found that cAMP response element-binding protein (CREB) phosphorylation was significantly increased with (Bu)(2)cAMP treatment. Furthermore, cotransfection with the catalytic subunit of PKA and/or the transcription factor CREB caused a significant stimulation on leptin promoter activity. On the other hand, the cotransfection with a dominant negative mutant of the regulatory subunit of PKA inhibited leptin promoter activity. We determined that cAMP effect could be blocked by pharmacologic inhibition of PKA or adenylyl ciclase in BeWo cells and in human placental explants. Thereafter, we decided to investigate the involvement of the MAPK/ERK signaling pathway in the cAMP effect on leptin induction. We found that 50 microm PD98059, a MAPK kinase inhibitor, partially blocked leptin induction by cAMP, measured both by Western blot analysis and reporter transient transfection assay. Moreover, ERK 1/2 phosphorylation was significantly increased with (Bu)(2)cAMP treatment, and this effect was dose dependent. Finally, we observed that 50 microm PD98059 inhibited cAMP-dependent phosphorylation of CREB in placental explants. In summary, we provide some evidence suggesting that cAMP induces leptin expression in placental cells and that this effect seems to be mediated by a cross talk between PKA and MAPK signaling pathways.

Up-regulation of placental leptin by human chorionic gonadotropin.

Leptin, the 16,000 molecular weight protein product of the obese gene, was originally considered as an adipocyte-derived signaling molecule for the central control of metabolism. However, leptin has been suggested to be involved in other functions during pregnancy, particularly in placenta, in which it was found to be expressed. In the present work, we have found that recombinant human chorionic gonadotropin (hCG) added to BeWo choriocarcinoma cell line showed a stimulatory effect on endogenous leptin expression, when analyzed by Western blot. This effect was time and dose dependent. Maximal effect was achieved at hCG 100 IU/ml. Moreover, hCG treatment enhanced leptin promoter activity up to 12.9 times, evaluated by transient transfection with a plasmid construction containing different promoter regions and the reporter gene luciferase. This effect was dose dependent and evidenced with all the promoter regions analyzed, regardless of length. Similar results were obtained with placental explants, thus indicating physiological relevance. Because hCG signal transduction usually involves cAMP signaling, this pathway was analyzed. Contrarily, we found that dibutyryl cAMP counteracted hCG effect on leptin expression. Furthermore, cotransfection with the catalytic subunit of PKA and/or the transcription factor cAMP response element binding protein repressed leptin expression. Thereafter we determined that hCG effect could be partially blocked by pharmacologic inhibition of MAPK pathway with 50 microM PD98059 but not by the inhibition of the phosphatidylinositol 3-kinase pathway with 0.1 microm wortmannin. Moreover, hCG treatment promoted MAPK kinase and ERK1/ERK2 phosphorylation in placental cells. Finally, cotransfection with a dominant-negative mutant of MAPK blocked the hCG-mediated activation of leptin expression. In conclusion, we provide some evidence suggesting that hCG induces leptin expression in trophoblastic cells probably involving the MAPK signal transduction pathway.

Prognostic impact of MGMT promoter methylation and MGMT and CD133 expression in colorectal adenocarcinoma.

BACKGROUND New biomarkers are needed for the prognosis of advanced colorectal cancer, which remains incurable by conventional treatments. O6-methylguanine DNA methyltransferase (MGMT) methylation and protein expression have been related to colorectal cancer treatment failure and tumor progression. Moreover, the presence in these tumors of cancer stem cells, which are characterized by CD133 expression, has been associated with chemoresistance, radioresistance, metastasis, and local recurrence. The objective of this study was to determine the prognostic value of CD133 and MGMT and their possible interaction in colorectal cancer patients. METHODS MGMT and CD133 expression was analyzed by immunohistochemistry in 123 paraffin-embedded colorectal adenocarcinoma samples, obtaining the percentage staining and intensity. MGMT promoter methylation status was obtained by using bisulfite modification and methylation-specific PCR (MSP). These values were correlated with clinical data, including overall survival (OS), disease-free survival (DFS), tumor stage, and differentiation grade. RESULTS Low MGMT expression intensity was significantly correlated with shorter OS and was a prognostic factor independently of treatment and histopathological variables. High percentage of CD133 expression was significantly correlated with shorter DFS but was not an independent factor. Patients with low-intensity MGMT expression and ≥50% CD133 expression had the poorest DFS and OS outcomes. CONCLUSIONS Our results support the hypothesis that MGMT expression may be an OS biomarker as useful as tumor stage or differentiation grade and that CD133 expression may be a predictive biomarker of DFS. Thus, MGMT and CD133 may both be useful for determining the prognosis of colorectal cancer patients and to identify those requiring more aggressive adjuvant therapies. Future studies will be necessary to determine its clinical utility.

The N-terminal DNA-binding 'zinc finger' of the oestrogen and glucocorticoid receptors determines target gene specificity.

Steroid hormone receptors activate specific gene transcription by binding as hormone-receptor complexes to short DNA enhancer-like elements termed hormone response elements (HREs). We have shown previously that a highly conserved 66 amino acid region of the oestrogen (ER) and glucocorticoid (GR) receptors, which corresponds to part of the receptor DNA binding domain (region C) is responsible for determining the specificity of target gene activation. This region contains two sub-regions (CI and CII) analogous to the 'zinc-fingers' of the transcription factor TFIIIA. We show here that CI and CII appear to be separate domains both involved in DNA binding. Furthermore, using chimaeric ERs in which either the first (N-terminal) (CI) or second (CII) 'zinc finger' region has been exchanged with that of the GR, indicates that it is the first 'zinc finger' which largely determines target gene specificity. We suggest that receptor recognition of the HRE is analogous to that of the helix-turn-helix DNA binding motif in that the receptor binds to DNA as a dimer with the first 'zinc finger' lying in the major groove recognizing one half of the palindromic HRE, and that protein-DNA interaction is stabilized through non-specific DNA binding and dimer interactions contributed by the second 'zinc finger'.

Point Mutations in the Monocarboxylate Transporter SLC16A12 Lead to Juvenile and Age-Related Cataract

Purpose: Previously we reported on a premature termination mutation in SLC16A12 that leads to dominant juvenile cataract and renal glucosuria. To assess the mutation rate and genotype-phenotype correlations of SLC16A12 in juvenile or age-related forms of cataract, we performed a mutation screen in cataract patients. Methods: Clinical data of approximately 660 patients were collected, genomic DNA was isolated and analyzed. Exons 3 to 8 including flanking intron sequences of SLC16A12 were PCR amplified and DNA sequence was determined. Selected mutations were tested by cell culture assays, in silico analysis and RT-PCR. Results: We found sequence alterations at a rate of approximately 1/75 patients. None of them was found in 360 control alleles. Alterations affect splice site and regulatory region but most mutations caused an amino acid substitution. The majority of the coding region mutations maps to trans-membrane domains. One mutation located to the 5'UTR. It affects translational efficiency of SLC16A12. In addition, we identified a cataract-predisposing SNP in the non-coding region that causes allele-specific splicing of the 5'UTR region. Conclusions: Altered translational efficiency of the solute carrier SLC16A12 and its allele-specific splicing strongly support a model of challenged homeostasis to cause various forms of cataract. In addition, the pathogenic property of the here reported sequence alterations is supported by the lack of known sequence variations within the coding region of SLC16A12. Due to the relatively high mutation rate, we suggest to include SLC16A12 in diagnostic cataract screening. Generally, our data recommend the assessment of regulatory sequences for diagnostic purposes.

Pancreatic-specific expression of the glucose transporter type 2 gene: identification of cis-elements and islet-specific trans-acting factors.

A defect in glucose sensing of the pancreatic beta-cells has been observed in several animal models of type II diabetes and has been correlated with a reduced gene expression of the glucose transporter type 2 (Glut2). In a transgenic mouse model, expression of Glut2 antisense RNA in pancreatic beta-cells has recently been shown to be associated with an impaired glucose-induced insulin secretion and the development of diabetes. To identify factors that may be involved in the specific decrease of Glut2 in the beta-cells of the diabetic animal, an attempt was made to localize the cis-elements and trans-acting factors involved in the control of Glut2 expression in the endocrine pancreas. It was demonstrated by transient transfection studies that only 338 base pairs (bp) of the murine Glut2 proximal promoter are needed for reporter gene expression in pancreatic islet-derived cell lines, whereas no activity was detected in nonpancreatic cells. Three cis-elements, GTI, GTII, and GTIII, have been identified by DNAse I footprinting and gel retardation experiments within these 338 bp. GTI and GTIII bind distinct but ubiquitously expressed trans-acting factors. On the other hand, nuclear proteins specifically expressed in pancreatic cell lines interact with GTII, and their relative abundance correlates with endogenous Glut2 expression. These GTII-binding factors correspond to nuclear proteins of 180 and 90 kilodaltons as defined by Southwestern analysis. The 180-kilodalton factor is present in pancreatic beta-cell lines but not in an alpha-cell line. Mutation of the GTI or GTIII cis-elements decreases transcriptional activity directed by the 338-bp promoter, whereas mutation of GTII increases gene transcription. Thus negative and positive regulatory sequences are identified within the proximal 338 bp of the GLUT2 promoter and may participate in the islet-specific expression of the gene by binding beta-cell specific trans-acting factors.

Differentiated evolutionary rates in alternative exons and the implications for splicing regulation

Background: Alternatively spliced exons play an important role in the diversification of gene function in most metazoans and are highly regulated by conserved motifs in exons and introns. Two contradicting properties have been associated to evolutionary conserved alternative exons: higher sequence conservation and higher rate of non-synonymous substitutions, relative to constitutive exons. In order to clarify this issue, we have performed an analysis of the evolution of alternative and constitutive exons, using a large set of protein coding exons conserved between human and mouse and taking into account the conservation of the transcript exonic structure. Further, we have also defined a measure of the variation of the arrangement of exonic splicing enhancers (ESE-conservation score) to study the evolution of splicing regulatory sequences. We have used this measure to correlate the changes in the arrangement of ESEs with the divergence of exon and intron sequences. Results: We find evidence for a relation between the lack of conservation of the exonic structure and the weakening of the sequence evolutionary constraints in alternative and constitutive exons. Exons in transcripts with non-conserved exonic structures have higher synonymous (dS) and non-synonymous (dN) substitution rates than exons in conserved structures. Moreover, alternative exons in transcripts with non-conserved exonic structure are the least constrained in sequence evolution, and at high EST-inclusion levels they are found to be very similar to constitutive exons, whereas alternative exons in transcripts with conserved exonic structure have a dS significantly lower than average at all EST-inclusion levels. We also find higher conservation in the arrangement of ESEs in constitutive exons compared to alternative ones. Additionally, the sequence conservation at flanking introns remains constant for constitutive exons at all ESE-conservation values, but increases for alternative exons at high ESE-conservation values. Conclusion: We conclude that most of the differences in dN observed between alternative and constitutive exons can be explained by the conservation of the transcript exonic structure. Low dS values are more characteristic of alternative exons with conserved exonic structure, but not of those with non-conserved exonic structure. Additionally, constitutive exons are characterized by a higher conservation in the arrangement of ESEs, and alternative exons with an ESE-conservation similar to that of constitutive exons are characterized by a conservation of the flanking intron sequences higher than average, indicating the presence of more intronic regulatory signals.

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.