Biomedical subjects

Christopher B Burge

Publications and source records attributed to Christopher B Burge.

At least 19 recordsLinked to original sources

PubMed2007

Weak definition of IKBKAP exon 20 leads to aberrant splicing in familial dysautonomia.

Splicing mutations that lead to devastating genetic diseases are often located in nonconserved or weakly conserved sequences that normally do not affect splicing. Thus, the underlying reason for the splicing defect is not immediately obvious. An example of this phenomenon is observed in the neurodevelopmental disease familial dysautonomia (FD), which is caused by a single-base change in the 5' splice site (5'ss) of intron 20 in the IKBKAP gene (c.2204+6T>C). This mutation, which is in the sixth position of the intron and results in exon 20 skipping, has no phenotype in many other introns. To determine why the position 6 mutation causes aberrant splicing only in certain cases, we first used an in silico approach to identify potential sequences involved in exon 20 skipping. Computational analyses of the exon 20 5'ss itself predicted that this nine-nucleotide splicing signal, even when it contains the T>C mutation, is not sufficiently weak to explain the FD phenotype. However, the computational analysis predicted that both the upstream 3' splice site (3'ss) and exon 20 contain weak splicing signals, indicating that the FD 5'ss, together with the surrounding splicing signals, are not adequate for defining exon 20. These in silico predictions were corroborated using IKBKAP minigenes in a new rapid and simple in vitro coupled RNA polymerase (RNAP) II transcription/splicing assay. Finally, the weak splicing signals that flank the T>C mutation were validated as the underlying cause of familial dysautonomia in vivo using transient transfection assays. Together, our study demonstrates the general utility of combining in silico data with an in vitro RNAP II transcription/splicing system for rapidly identifying critical sequences that underlie the numerous splicing diseases caused by otherwise silent mutations.

Alternative Splicing↗

PubMed2006

Inference of splicing regulatory activities by sequence neighborhood analysis.

Sequence-specific recognition of nucleic-acid motifs is critical to many cellular processes. We have developed a new and general method called Neighborhood Inference (NI) that predicts sequences with activity in regulating a biochemical process based on the local density of known sites in sequence space. Applied to the problem of RNA splicing regulation, NI was used to predict hundreds of new exonic splicing enhancer (ESE) and silencer (ESS) hexanucleotides from known human ESEs and ESSs. These predictions were supported by cross-validation analysis, by analysis of published splicing regulatory activity data, by sequence-conservation analysis, and by measurement of the splicing regulatory activity of 24 novel predicted ESEs, ESSs, and neutral sequences using an in vivo splicing reporter assay. These results demonstrate the ability of NI to accurately predict splicing regulatory activity and show that the scope of exonic splicing regulatory elements is substantially larger than previously anticipated. Analysis of orthologous exons in four mammals showed that the NI score of ESEs, a measure of function, is much more highly conserved above background than ESE primary sequence. This observation indicates a high degree of selection for ESE activity in mammalian exons, with surprisingly frequent interchangeability between ESE sequences.

Binding Sites↗

PubMed2006

Detection of broadly expressed neuronal genes in C. elegans.

The genes that are expressed in most or all types of neurons define generic neuronal features and provide a window into the developmental origin and function of the nervous system. Few such genes (sometimes referred to as pan-neuronal or broadly expressed neuronal genes) have been defined to date and the mechanisms controlling their regulation are not well understood. As a first step in investigating their regulation, we used a computational approach to detect sequences overrepresented in their promoter elements. We identified a ten-nucleotide cis-regulatory motif shared by many broadly expressed neuronal genes and demonstrated that it is involved in control of neuronal expression. Our results further suggest that global and cell-type-specific controls likely act in concert to establish pan-neuronal gene expression. Using the newly discovered motif and genome-level gene expression data, we identified a set of 234 candidate broadly expressed genes. The known involvement of many of these genes in neurogenesis and physiology of the nervous system supports the utility of this set for future targeted analyses.

Animals↗

PubMed2006

General and specific functions of exonic splicing silencers in splicing control.

Correct splice site recognition is critical in pre-mRNA splicing. We find that almost all of a diverse panel of exonic splicing silencer (ESS) elements alter splice site choice when placed between competing sites, consistently inhibiting use of intron-proximal 5' and 3' splice sites. Supporting a general role for ESSs in splice site definition, we found that ESSs are both abundant and highly conserved between alternative splice site pairs and that mutation of ESSs located between natural alternative splice site pairs consistently shifted splicing toward the intron-proximal site. Some exonic splicing enhancers (ESEs) promoted use of intron-proximal 5' splice sites, and tethering of hnRNP A1 and SF2/ASF proteins between competing splice sites mimicked the effects of ESS and ESE elements, respectively. Further, we observed that specific subsets of ESSs had distinct effects on a multifunctional intron retention reporter and that one of these subsets is likely preferred for regulation of endogenous intron retention events. Together, our findings provide a comprehensive picture of the functions of ESSs in the control of diverse types of splicing decisions.

3' Untranslated Regions↗

PubMed2006

Untangling influences of hydrophobicity on protein sequences and structures.

We perform a statistical analysis of solvent accessibility and hydrophobicity profiles of a representative set of proteins. The joint probability distribution is well fitted to a multivariable Gaussian, which takes a relatively simple form when expressed in terms of the Fourier transforms of the profiles. This allows us to quantify the asymmetric manner by which these profiles influence each other. For example, the alpha-helix periodicity in sequence hydrophobicity is dictated by the solvent accessibility of structures, and not vice versa, possibly indicating the faster evolution of sequences compared to structures. The decorrelated hydrophobicity and solvent accessibility profiles show distinct behaviors at long periods, where sequence hydrophobicity fluctuates less, while solvent accessibility fluctuates more than average. The correlations between the two profiles can be interpreted as the Boltzmann weight of the solvation energy at room temperature, consistent with earlier observations.

Amino Acid Sequence↗

PubMed2006

HOLLYWOOD: a comparative relational database of alternative splicing.

RNA splicing is an essential step in gene expression, and is often variable, giving rise to multiple alternatively spliced mRNA and protein isoforms from a single gene locus. The design of effective databases to support experimental and computational investigations of alternative splicing (AS) is a significant challenge. In an effort to integrate accurate exon and splice site annotation with current knowledge about splicing regulatory elements and predicted AS events, and to link information about the splicing of orthologous genes in different species, we have developed the Hollywood system. This database was built upon genomic annotation of splicing patterns of known genes derived from spliced alignment of complementary DNAs (cDNAs) and expressed sequence tags, and links features such as splice site sequence and strength, exonic splicing enhancers and silencers, conserved and non-conserved patterns of splicing, and cDNA library information for inferred alternative exons. Hollywood was implemented as a relational database and currently contains comprehensive information for human and mouse. It is accompanied by a web query tool that allows searches for sets of exons with specific splicing characteristics or splicing regulatory element composition, or gives a graphical or sequence-level summary of splicing patterns for a specific gene. A streamlined graphical representation of gene splicing patterns is provided, and these patterns can alternatively be layered onto existing information in the UCSC Genome Browser. The database is accessible at http://hollywood.mit.edu.

Alternative Splicing↗

PubMed2005

The widespread impact of mammalian MicroRNAs on mRNA repression and evolution.

Thousands of mammalian messenger RNAs are under selective pressure to maintain 7-nucleotide sites matching microRNAs (miRNAs). We found that these conserved targets are often highly expressed at developmental stages before miRNA expression and that their levels tend to fall as the miRNA that targets them begins to accumulate. Nonconserved sites, which outnumber the conserved sites 10 to 1, also mediate repression. As a consequence, genes preferentially expressed at the same time and place as a miRNA have evolved to selectively avoid sites matching the miRNA. This phenomenon of selective avoidance extends to thousands of genes and enables spatial and temporal specificities of miRNAs to be revealed by finding tissues and developmental stages in which messages with corresponding sites are expressed at lower levels.

Animals↗

PubMed2005

A transcriptional response to replication status mediated by the conserved bacterial replication protein DnaA.

Organisms respond to perturbations in DNA replication. We characterized the global transcriptional response to inhibition of DNA replication in Bacillus subtilis. We focused on changes that were independent of the known recA-dependent global DNA damage (SOS) response. We found that overlapping sets of genes are affected by perturbations in replication elongation or initiation and that this transcriptional response serves to inhibit cell division and maintain cell viability. Approximately 20 of the operons (>50 genes) affected have potential DnaA-binding sites and are probably regulated directly by DnaA, the highly conserved replication initiation protein and transcription factor. Many of these genes have homologues and recognizable DnaA-binding sites in other bacteria, indicating that a DnaA-mediated response, elicited by changes in DNA replication status, may be conserved.

Bacillus subtilis↗

PubMed2005

Recognition of unknown conserved alternatively spliced exons.

The split structure of most mammalian protein-coding genes allows for the potential to produce multiple different mRNA and protein isoforms from a single gene locus through the process of alternative splicing (AS). We propose a computational approach called UNCOVER based on a pair hidden Markov model to discover conserved coding exonic sequences subject to AS that have so far gone undetected. Applying UNCOVER to orthologous introns of known human and mouse genes predicts skipped exons or retained introns present in both species, while discriminating them from conserved noncoding sequences. The accuracy of the model is evaluated on a curated set of genes with known conserved AS events. The prediction of skipped exons in the approximately 1% of the human genome represented by the ENCODE regions leads to more than 50 new exon candidates. Five novel predicted AS exons were validated by RT-PCR and sequencing analysis of 15 introns with strong UNCOVER predictions and lacking EST evidence. These results imply that a considerable number of conserved exonic sequences and associated isoforms are still completely missing from the current annotation of known genes. UNCOVER also identifies a small number of candidates for conserved intron retention.

Journal Article↗

PubMed2005

A combinatorial code for splicing silencing: UAGG and GGGG motifs.

Alternative pre-mRNA splicing is widely used to regulate gene expression by tuning the levels of tissue-specific mRNA isoforms. Few regulatory mechanisms are understood at the level of combinatorial control despite numerous sequences, distinct from splice sites, that have been shown to play roles in splicing enhancement or silencing. Here we use molecular approaches to identify a ternary combination of exonic UAGG and 5'-splice-site-proximal GGGG motifs that functions cooperatively to silence the brain-region-specific CI cassette exon (exon 19) of the glutamate NMDA R1 receptor (GRIN1) transcript. Disruption of three components of the motif pattern converted the CI cassette into a constitutive exon, while predominant skipping was conferred when the same components were introduced, de novo, into a heterologous constitutive exon. Predominant exon silencing was directed by the motif pattern in the presence of six competing exonic splicing enhancers, and this effect was retained after systematically repositioning the two exonic UAGGs within the CI cassette. In this system, hnRNP A1 was shown to mediate silencing while hnRNP H antagonized silencing. Genome-wide computational analysis combined with RT-PCR testing showed that a class of skipped human and mouse exons can be identified by searches that preserve the sequence and spatial configuration of the UAGG and GGGG motifs. This analysis suggests that the multi-component silencing code may play an important role in the tissue-specific regulation of the CI cassette exon, and that it may serve more generally as a molecular language to allow for intricate adjustments and the coordination of splicing patterns from different genes.

Alternative Splicing↗

PubMed2005

Linking C5 deficiency to an exonic splicing enhancer mutation.

As an important component of the innate immune system, complement provides the initial response to prevent infections by pathogenic microorganisms. Patients with dysfunction of C5 display a propensity for severe recurrent infections. In this study, we present a patient with C5 deficiency demonstrated by immunochemical and functional analyses. Direct sequencing of all C5 exons displayed no mutation of obvious functional significance, except for an A to G transition in exon 10 predicting an exchange from lysine to arginine. This sequence alteration was present in only one allele of family members with a reduced serum C5 concentration and in both alleles of the patient with almost complete C5 deficiency, suggesting that this alteration may be producing the phenotype. Recent findings indicate that distinct nucleotide sequences, termed exonic splicing enhancers (ESEs), influence the splicing process. cDNA from all family members harboring the mutated allele showed skipping of exon 10, which resulted in a premature STOP codon, explaining the lack of C5 in the propositus. Sequence analysis of the mutated region revealed the substitution to be located within an ESE, as predicted by the RESCUE-ESE program. The altered ESE sequence is located close to the 5' splicing site and also lowers the predicted strength of the splice site itself. This apparently inconsequential sequence alteration represents a noncanonical splicing mutation altering an ESE. Our finding sheds a new light on the role of putative silent/conservative mutations in disease-associated genes.

Alternative Splicing↗

PubMed2005

Identification and analysis of alternative splicing events conserved in human and mouse.

Alternative pre-mRNA splicing affects a majority of human genes and plays important roles in development and disease. Alternative splicing (AS) events conserved since the divergence of human and mouse are likely of primary biological importance, but relatively few of such events are known. Here we describe sequence features that distinguish exons subject to evolutionarily conserved AS, which we call alternative conserved exons (ACEs), from other orthologous human/mouse exons and integrate these features into an exon classification algorithm, acescan. Genome-wide analysis of annotated orthologous human-mouse exon pairs identified approximately 2,000 predicted ACEs. Alternative splicing was verified in both human and mouse tissues by using an RT-PCR-sequencing protocol for 21 of 30 (70%) predicted ACEs tested, supporting the validity of a majority of acescan predictions. By contrast, AS was observed in mouse tissues for only 2 of 15 (13%) tested exons that had EST or cDNA evidence of AS in human but were not predicted ACEs, and AS was never observed for 11 negative control exons in human or mouse tissues. Predicted ACEs were much more likely to preserve the reading frame and less likely to disrupt protein domains than other AS events and were enriched in genes expressed in the brain and in genes involved in transcriptional regulation, RNA processing, and development. Our results also imply that the vast majority of AS events represented in the human EST database are not conserved in mouse.

Alternative Splicing↗

PubMed2005

Conserved seed pairing, often flanked by adenosines, indicates that thousands of human genes are microRNA targets.

We predict regulatory targets of vertebrate microRNAs (miRNAs) by identifying mRNAs with conserved complementarity to the seed (nucleotides 2-7) of the miRNA. An overrepresentation of conserved adenosines flanking the seed complementary sites in mRNAs indicates that primary sequence determinants can supplement base pairing to specify miRNA target recognition. In a four-genome analysis of 3' UTRs, approximately 13,000 regulatory relationships were detected above the estimate of false-positive predictions, thereby implicating as miRNA targets more than 5300 human genes, which represented 30% of our gene set. Targeting was also detected in open reading frames. In sum, well over one third of human genes appear to be conserved miRNA targets.

3' Untranslated Regions↗

PubMed2004

Systematic identification and analysis of exonic splicing silencers.

Exonic splicing silencers (ESSs) are cis-regulatory elements that inhibit the use of adjacent splice sites, often contributing to alternative splicing (AS). To systematically identify ESSs, an in vivo splicing reporter system was developed to screen a library of random decanucleotides. The screen yielded 141 ESS decamers, 133 of which were unique. The silencer activity of over a dozen of these sequences was also confirmed in a heterologous exon/intron context and in a second cell type. Of the unique ESS decamers, most could be clustered into groups to yield seven putative ESS motifs, some resembling known motifs bound by hnRNPs H and A1. Potential roles of ESSs in constitutive splicing were explored using an algorithm, ExonScan, which simulates splicing based on known or putative splicing-related motifs. ExonScan and related bioinformatic analyses suggest that these ESS motifs play important roles in suppression of pseudoexons, in splice site definition, and in AS.

Algorithms↗

PubMed2004

Patterns of intron gain and loss in fungi.

Little is known about the patterns of intron gain and loss or the relative contributions of these two processes to gene evolution. To investigate the dynamics of intron evolution, we analyzed orthologous genes from four filamentous fungal genomes and determined the pattern of intron conservation. We developed a probabilistic model to estimate the most likely rates of intron gain and loss giving rise to these observed conservation patterns. Our data reveal the surprising importance of intron gain. Between about 150 and 250 gains and between 150 and 350 losses were inferred in each lineage. We discuss one gene in particular (encoding 1-phosphoribosyl-5-pyrophosphate synthetase) that displays an unusually high rate of intron gain in multiple lineages. It has been recognized that introns are biased towards the 5' ends of genes in intron-poor genomes but are evenly distributed in intron-rich genomes. Current models attribute this bias to 3' intron loss through a poly-adenosine-primed reverse transcription mechanism. Contrary to standard models, we find no increased frequency of intron loss toward the 3' ends of genes. Thus, recent intron dynamics do not support a model whereby 5' intron positional bias is generated solely by 3'-biased intron loss.

Adenosine↗

PubMed2004

Variation in sequence and organization of splicing regulatory elements in vertebrate genes.

Although core mechanisms and machinery of premRNA splicing are conserved from yeast to human, the details of intron recognition often differ, even between closely related organisms. For example, genes from the pufferfish Fugu rubripes generally contain one or more introns that are not properly spliced in mouse cells. Exploiting available genome sequence data, a battery of sequence analysis techniques was used to reach several conclusions about the organization and evolution of splicing regulatory elements in vertebrate genes. The classical splice site and putative branch site signals are completely conserved across the vertebrates studied (human, mouse, pufferfish, and zebrafish), and exonic splicing enhancers also appear broadly conserved in vertebrates. However, another class of splicing regulatory elements, the intronic splicing enhancers, appears to differ substantially between mammals and fish, with G triples (GGG) very abundant in mammalian introns but comparatively rare in fish. Conversely, short repeats of AC and GT are predicted to function as intronic splicing enhancers in fish but are not enriched in mammalian introns. Consistent with this pattern, exonic splicing enhancer-binding SR proteins are highly conserved across all vertebrates, whereas heterogeneous nuclear ribonucleoproteins, which bind many intronic sequences, vary in domain structure and even presence/absence between mammals and fish. Exploiting differences in intronic sequence composition, a statistical model was developed to predict the splicing phenotype of Fugu introns in mammalian systems and was used to engineer the spliceability of a Fugu intron in human cells by insertion of specific sequences, thereby rescuing splicing in human cells.

Animals↗

PubMed2004

Variation in alternative splicing across human tissues.

BACKGROUND: Alternative pre-mRNA splicing (AS) is widely used by higher eukaryotes to generate different protein isoforms in specific cell or tissue types. To compare AS events across human tissues, we analyzed the splicing patterns of genomically aligned expressed sequence tags (ESTs) derived from libraries of cDNAs from different tissues. RESULTS: Controlling for differences in EST coverage among tissues, we found that the brain and testis had the highest levels of exon skipping. The most pronounced differences between tissues were seen for the frequencies of alternative 3' splice site and alternative 5' splice site usage, which were about 50 to 100% higher in the liver than in any other human tissue studied. Quantifying differences in splice junction usage, the brain, pancreas, liver and the peripheral nervous system had the most distinctive patterns of AS. Analysis of available microarray expression data showed that the liver had the most divergent pattern of expression of serine-arginine protein and heterogeneous ribonucleoprotein genes compared to the other human tissues studied, possibly contributing to the unusually high frequency of alternative splice site usage seen in liver. Sequence motifs enriched in alternative exons in genes expressed in the brain, testis and liver suggest specific splicing factors that may be important in AS regulation in these tissues. CONCLUSIONS: This study distinguishes the human brain, testis and liver as having unusually high levels of AS, highlights differences in the types of AS occurring commonly in different tissues, and identifies candidate cis-regulatory elements and trans-acting factors likely to have important roles in tissue-specific AS in human cells.

Alternative Splicing↗

PubMed2004

Single nucleotide polymorphism-based validation of exonic splicing enhancers.

Because deleterious alleles arising from mutation are filtered by natural selection, mutations that create such alleles will be underrepresented in the set of common genetic variation existing in a population at any given time. Here, we describe an approach based on this idea called VERIFY (variant elimination reinforces functionality), which can be used to assess the extent of natural selection acting on an oligonucleotide motif or set of motifs predicted to have biological activity. As an application of this approach, we analyzed a set of 238 hexanucleotides previously predicted to have exonic splicing enhancer (ESE) activity in human exons using the relative enhancer and silencer classification by unanimous enrichment (RESCUE)-ESE method. Aligning the single nucleotide polymorphisms (SNPs) from the public human SNP database to the chimpanzee genome allowed inference of the direction of the mutations that created present-day SNPs. Analyzing the set of SNPs that overlap RESCUE-ESE hexamers, we conclude that nearly one-fifth of the mutations that disrupt predicted ESEs have been eliminated by natural selection (odds ratio = 0.82 +/- 0.05). This selection is strongest for the predicted ESEs that are located near splice sites. Our results demonstrate a novel approach for quantifying the extent of natural selection acting on candidate functional motifs and also suggest certain features of mutations/SNPs, such as proximity to the splice site and disruption or alteration of predicted ESEs, that should be useful in identifying variants that might cause a biological phenotype.

Alleles↗