Using information content and base frequencies to distinguish mutations from genetic polymorphisms in splice junction recognition sites

A proposed molecular mechanism for pathogenesis of severe RNA-viral pulmonary infections

Background: Certain riboviruses can cause severe pulmonary complications leading to death in some infected patients. We propose that DNA damage induced-apoptosis accelerates viral release, triggered by depletion of host RNA binding proteins (RBPs) from nuclear RNA bound to replicating viral sequences. Methods: Information theory-based analysis of interactions between RBPs and individual sequences in the Severe Acute Respiratory Syndrome CoronaVirus 2 (SARS-CoV-2), Influenza A (H3N1), HIV-1, and Dengue genomes identifies strong RBP binding sites in these viral genomes. Replication and expression of viral sequences is expected to increasingly sequester RBPs - SRSF1 and RNPS1. Ordinarily, RBPs bound to nascent host transcripts prevents their annealing to complementary DNA. Their depletion induces destabilizing R-loops. Chromosomal breakage occurs when an excess of unresolved R-loops collide with incoming replication forks, overwhelming the DNA repair machinery. We estimated stoichiometry of inhibition of RBPs in host nuclear RNA by counting competing binding sites in replicating viral genomes and host RNA. Results: Host RBP binding sites are frequent and conserved among different strains of RNA viral genomes. Similar binding motifs of SRSF1 and RNPS1 explain why DNA damage resulting from SRSF1 depletion is complemented by expression of RNPS1. Clustering of strong RBP binding sites coincides with the distribution of RNA-DNA hybridization sites across the genome. SARS-CoV-2 replication is estimated to require 32.5-41.8 hours to effectively compete for binding of an equal proportion of SRSF1 binding sites in host encoded nuclear RNAs. Significant changes in expression of transcripts encoding DNA repair and apoptotic proteins were found in an analysis of influenza A and Dengue-infected cells in some individuals. Conclusions: R-loop-induced apoptosis indirectly resulting from viral replication could release significant quantities of membrane-associated virions into neighboring alveoli. These could infect adjacent pneumocytes and other tissues, rapidly compromising lung function, causing multiorgan system failure and other described symptoms.

A proposed molecular mechanism for pathogenesis of severe RNA-viral pulmonary infections

Background: Certain riboviruses can cause severe pulmonary complications leading to death in some infected patients. We propose that DNA damage induced-apoptosis accelerates viral release, triggered by depletion of host RNA binding proteins (RBPs) from nuclear RNA bound to replicating viral sequences. Methods: Information theory-based analysis of interactions between RBPs and individual sequences in the Severe Acute Respiratory Syndrome CoronaVirus 2 (SARS-CoV-2), Influenza A (H3N2), HIV-1, and Dengue genomes identifies strong RBP binding sites in these viral genomes. Replication and expression of viral sequences is expected to increasingly sequester RBPs - SRSF1 and RNPS1. Ordinarily, RBPs bound to nascent host transcripts prevents their annealing to complementary DNA. Their depletion induces destabilizing R-loops. Chromosomal breakage occurs when an excess of unresolved R-loops collide with incoming replication forks, overwhelming the DNA repair machinery. We estimated stoichiometry of inhibition of RBPs in host nuclear RNA by counting competing binding sites in replicating viral genomes and host RNA. Results: Host RBP binding sites are frequent and conserved among different strains of RNA viral genomes. Similar binding motifs of SRSF1 and RNPS1 explain why DNA damage resulting from SRSF1 depletion is complemented by expression of RNPS1. Clustering of strong RBP binding sites coincides with the distribution of RNA-DNA hybridization sites across the genome. SARS-CoV-2 replication is estimated to require 32.5-41.8 hours to effectively compete for binding of an equal proportion of SRSF1 binding sites in host encoded nuclear RNAs. Significant changes in expression of transcripts encoding DNA repair and apoptotic proteins were found in an analysis of influenza A and Dengue-infected cells in some individuals. Conclusions: R-loop-induced apoptosis indirectly resulting from viral replication could release significant quantities of membrane-associated virions into neighboring alveoli. These could infect adjacent pneumocytes and other tissues, rapidly compromising lung function, causing multiorgan system failure and other described symptoms.

Expression Changes Confirm Genomic Variants Predicted to Result in Allele-Specific, Alternative mRNA Splicing

Splice isoform structure and abundance can be affected by either noncoding or masquerading coding variants that alter the structure or abundance of transcripts. When these variants are common in the population, these nonconstitutive transcripts are sufficiently frequent so as to resemble naturally occurring, alternative mRNA splicing. Prediction of the effects of such variants has been shown to be accurate using information theory-based methods. Single nucleotide polymorphisms (SNPs) predicted to significantly alter natural and/or cryptic splice site strength were shown to affect gene expression. Splicing changes for known SNP genotypes were confirmed in HapMap lymphoblastoid cell lines with gene expression microarrays and custom designed q-RT-PCR or TaqMan assays. The majority of these SNPs (15 of 22) as well as an independent set of 24 variants were then subjected to RNAseq analysis using the ValidSpliceMut web beacon (http://validsplicemut.cytognomix.com), which is based on data from the Cancer Genome Atlas and International Cancer Genome Consortium. SNPs from different genes analyzed with gene expression microarray and q-RT-PCR exhibited significant changes in affected splice site use. Thirteen SNPs directly affected exon inclusion and 10 altered cryptic site use. Homozygous SNP genotypes resulting in stronger splice sites exhibited higher levels of processed mRNA than alleles associated with weaker sites. Four SNPs exhibited variable expression among individuals with the same genotypes, masking statistically significant expression differences between alleles. Genome-wide information theory and expression analyses (RNAseq) in tumor exomes and genomes confirmed splicing effects for 7 of the HapMap SNP and 14 SNPs identified from tumor genomes. q-RT-PCR resolved rare splice isoforms with read abundance too low for statistical significance in ValidSpliceMut. Nevertheless, the web-beacon provides evidence of unanticipated splicing outcomes, for example, intron retention due to compromised recognition of constitutive splice sites. Thus, ValidSpliceMut and q-RT-PCR represent complementary resources for identification of allele-specific, alternative splicing.

Pan-cancer repository of validated natural and cryptic mRNA splicing mutations

We present a major public resource of mRNA splicing mutations validated according to multiple lines of evidence of abnormal gene expression. Likely mutations present in all tumor types reported in the Cancer Genome Atlas (TCGA) and the International Cancer Genome Consortium (ICGC) were identified based on the comparative strengths of splice sites in tumor versus normal genomes, and then validated by respectively comparing counts of splice junction spanning and abundance of transcript reads in RNA-Seq data from matched tissues and tumors lacking these mutations. The comprehensive resource features 341,486 of these validated mutations, the majority of which (69.9%) are not present in the Single Nucleotide Polymorphism Database (dbSNP 150). There are 131,347 unique mutations which weaken or abolish natural splice sites, and 222,071 mutations which strengthen cryptic splice sites (11,932 affect both simultaneously). 28,812 novel or rare flagged variants (with <1% population frequency in dbSNP) were observed in multiple tumor tissue types. An algorithm was developed to classify variants into splicing molecular phenotypes that integrates germline heterozygosity, degree of information change and impact on expression. The classification thresholds were calibrated against the ClinVar clinical database phenotypic assignments. Variants are partitioned into allele-specific alternative splicing, likely aberrant and aberrant splicing phenotypes. Single variants or chromosome ranges can be queried using a Global Alliance for Genomics and Health (GA4GH)-compliant, web-based Beacon "Validated Splicing Mutations" either separately or in aggregate alongside other Beacons through the public Beacon Network, as well as through our website. The website provides additional information, such as a visual representation of supporting RNAseq results, gene expression in the corresponding normal tissues, and splicing molecular phenotypes.

Pan-cancer repository of validated natural and cryptic mRNA splicing mutations

We present a major public resource of mRNA splicing mutations validated according to multiple lines of evidence of abnormal gene expression. Likely mutations present in all tumor types reported in the Cancer Genome Atlas (TCGA) and the International Cancer Genome Consortium (ICGC) were identified based on the comparative strengths of splice sites in tumor versus normal genomes, and then validated by respectively comparing counts of splice junction spanning and abundance of transcript reads in RNA-Seq data from matched tissues and tumors lacking these mutations. The comprehensive resource features 341,486 of these validated mutations, the majority of which (69.9%) are not present in the Single Nucleotide Polymorphism Database (dbSNP 150). There are 131,347 unique mutations which weaken or abolish natural splice sites, and 222,071 mutations which strengthen cryptic splice sites (11,932 affect both simultaneously). 28,812 novel or rare flagged variants (with <1% population frequency in dbSNP) were observed in multiple tumor tissue types. An algorithm was developed to classify variants into splicing molecular phenotypes that integrates germline heterozygosity, degree of information change and impact on expression. The classification thresholds were calibrated against the ClinVar clinical database phenotypic assignments. Variants are partitioned into allele-specific alternative splicing, likely aberrant and aberrant splicing phenotypes. Single variants or chromosome ranges can be queried using a Global Alliance for Genomics and Health (GA4GH)-compliant, web-based Beacon "Validated Splicing Mutations" either separately or in aggregate alongside other Beacons through the public Beacon Network, as well as through our website. The website provides additional information, such as a visual representation of supporting RNAseq results, gene expression in the corresponding normal tissues, and splicing molecular phenotypes.

A heuristic model for computational prediction of human branch point sequence

BACKGROUND

Pre-mRNA splicing is the removal of introns from precursor mRNAs (pre-mRNAs) and the concurrent ligation of the flanking exons to generate mature mRNA. This process is catalyzed by the spliceosome, where the splicing factor 1 (SF1) specifically recognizes the seven-nucleotide branch point sequence (BPS) and the U2 snRNP later displaces the SF1 and binds to the BPS. In mammals, the degeneracy of BPS motifs together with the lack of a large set of experimentally verified BPSs complicates the task of BPS prediction in silico.

RESULTS

In this paper, we develop a simple and yet efficient heuristic model for human BPS prediction based on a novel scoring scheme, which quantifies the splicing strength of putative BPSs. The candidate BPS is restricted exclusively within a defined BPS search region to avoid the influences of other elements in the intron and therefore the prediction accuracy is improved. Moreover, using two types of relative frequencies for human BPS prediction, we demonstrate our model outperformed other current implementations on experimentally verified human introns.

CONCLUSION

We propose that the binding energy contributes to the molecular recognition involved in human pre-mRNA splicing. In addition, a genome-wide human BPS prediction is carried out. The characteristics of predicted BPSs are in accordance with experimentally verified human BPSs, and branch site positions relative to the 3'ss and the 5'end of the shortened AGEZ are consistent with the results of published papers. Meanwhile, a webserver for BPS predictor is freely available at http://biocomputer.bio.cuhk.edu.hk/BPS .

Prioritizing Variants in Complete Hereditary Breast and Ovarian Cancer Genes in Patients Lacking Known BRCA Mutations

BRCA1 and BRCA2 testing for hereditary breast and ovarian cancer (HBOC) does not identify all pathogenic variants. Sequencing of 20 complete genes in HBOC patients with uninformative test results (N = 287), including noncoding and flanking sequences of ATM, BARD1, BRCA1, BRCA2, CDH1, CHEK2, EPCAM, MLH1, MRE11A, MSH2, MSH6, MUTYH, NBN, PALB2, PMS2, PTEN, RAD51B, STK11, TP53, and XRCC2, identified 38,372 unique variants. We apply information theory (IT) to predict and prioritize noncoding variants of uncertain significance in regulatory, coding, and intronic regions based on changes in binding sites in these genes. Besides mRNA splicing, IT provides a common framework to evaluate potential affinity changes in transcription factor (TFBSs), splicing regulatory (SRBSs), and RNA-binding protein (RBBSs) binding sites following mutation. We prioritized variants affecting the strengths of 10 splice sites (four natural, six cryptic), 148 SRBS, 36 TFBS, and 31 RBBS. Three variants were also prioritized based on their predicted effects on mRNA secondary (2°) structure and 17 for pseudoexon activation. Additionally, four frameshift, two in-frame deletions, and five stop-gain mutations were identified. When combined with pedigree information, complete gene sequence analysis can focus attention on a limited set of variants in a wide spectrum of functional mutation types for downstream functional and co-segregation analysis.

Interpretation of mRNA splicing mutations in genetic disease: review of the literature and guidelines for information-theoretical analysis

The interpretation of genomic variants has become one of the paramount challenges in the post-genome sequencing era. In this review we summarize nearly 20 years of research on the applications of information theory (IT) to interpret coding and non-coding mutations that alter mRNA splicing in rare and common diseases. We compile and summarize the spectrum of published variants analyzed by IT, to provide a broad perspective of the distribution of deleterious natural and cryptic splice site variants detected, as well as those affecting splicing regulatory sequences. Results for natural splice site mutations can be interrogated dynamically with Splicing Mutation Calculator, a companion software program that computes changes in information content for any splice site substitution, linked to corresponding publications containing these mutations. The accuracy of IT-based analysis was assessed in the context of experimentally validated mutations. Because splice site information quantifies binding affinity, IT-based analyses can discern the differences between variants that account for the observed reduced (leaky) versus abolished mRNA splicing. We extend this principle by comparing predicted mutations in natural, cryptic, and regulatory splice sites with observed deleterious phenotypic and benign effects. Our analysis of 1727 variants revealed a number of general principles useful for ensuring portability of these analyses and accurate input and interpretation of mutations. We offer guidelines for optimal use of IT software for interpretation of mRNA splicing mutations.

Interpretation of mRNA splicing mutations in genetic disease: review of the literature and guidelines for information-theoretical analysis

The interpretation of genomic variants has become one of the paramount challenges in the post-genome sequencing era. In this review we summarize nearly 20 years of research on the applications of information theory (IT) to interpret coding and non-coding mutations that alter mRNA splicing in rare and common diseases. We compile and summarize the spectrum of published variants analyzed by IT, to provide a broad perspective of the distribution of deleterious natural and cryptic splice site variants detected, as well as those affecting splicing regulatory sequences. Results for natural splice site mutations can be interrogated dynamically with Splicing Mutation Calculator, a companion software program that computes changes in information content for any splice site substitution, linked to corresponding publications containing these mutations. The accuracy of IT-based analysis was assessed in the context of experimentally validated mutations. Because splice site information quantifies binding affinity, IT-based analyses can discern the differences between variants that account for the observed reduced (leaky) versus abolished mRNA splicing. We extend this principle by comparing predicted mutations in natural, cryptic, and regulatory splice sites with observed deleterious phenotypic and benign effects. Our analysis of 1727 variants revealed a number of general principles useful for ensuring portability of these analyses and accurate input and interpretation of mutations. We offer guidelines for optimal use of IT software for interpretation of mRNA splicing mutations.

Exploration of molecular factors impairing superoxide dismutase isoforms activity in human senile cataractous lenses

PURPOSE

To explore different molecular factors impairing the activities of superoxide dismutase (SOD) isoforms in senile cataractous lenses.

METHODS

Enzyme activity of SOD isoforms, levels of their corresponding cofactors copper (Cu), manganese (Mn), zinc (Zn), and expression of mRNA transcripts and proteins were determined in the lenses of human subjects with and without cataract. DNA from lens epithelium (LE) and peripheral blood was isolated. Polymerase chain reaction-single strand conformation polymorphism (PCR-SSCP) followed by sequencing was carried out to screen somatic mutations. The impact of intronic insertion/deletion (INDEL) variations on the splicing process and on the resultant transcript was evaluated. Genotyping of IVS4+42delG polymorphism of SOD1 gene was done by PCR-restriction fragment length polymorphism (RFLP).

RESULTS

A significant decrease in Cu/Zn- and Mn-SOD activity (P < 0.001) and in Cu/Zn-SOD transcript (P < 0.001) and its protein (P < 0.05) were found in cataractous lenses. No significant change in the level of copper (P = 0.36) and an increase in the level of manganese (P = 0.01) and zinc (P = 0.02) were observed in cataractous lenses. A significant positive correlation between the level of Cu/Zn-SOD activity and the levels of Cu (P = 0.003) and Zn (P = 0.005) was found in the cataractous lenses. DNA sequencing revealed three intronic INDEL variations in exon4 of SOD1 gene. Splice-junction analysis showed the potential of IVS4+42delG in creating a new cryptic acceptor site. If it is involved in alternate splicing, it could result in generation of SOD1 mRNA transcripts lacking exon4 region. Transcript analysis revealed the presence of complete SOD1 mRNA transcripts. Genotyping revealed the presence of IVS4+42delG polymorphism in all subjects.

CONCLUSIONS

The decrease in the activity of SOD1 isoform in cataractous lenses was associated with the decreased level of mRNA transcripts and their protein expression and was not associated with either modulation in the level of enzyme cofactors or with INDEL variations.

Pick one, but be quick: 5' splice sites and the problems of too many choices

Splice site selection is fundamental to pre-mRNA splicing and the expansion of genomic coding potential. 5' Splice sites (5'ss) are the critical elements at the 5' end of introns and are extremely diverse, as thousands of different sequences act as bona fide 5'ss in the human transcriptome. Most 5'ss are recognized by base-pairing with the 5' end of the U1 small nuclear RNA (snRNA). Here we review the history of research on 5'ss selection, highlighting the difficulties of establishing how base-pairing strength determines splicing outcomes. We also discuss recent work demonstrating that U1 snRNA:5'ss helices can accommodate noncanonical registers such as bulged duplexes. In addition, we describe the mechanisms by which other snRNAs, regulatory proteins, splicing enhancers, and the relative positions of alternative 5'ss contribute to selection. Moreover, we discuss mechanisms by which the recognition of numerous candidate 5'ss might lead to selection of a single 5'ss and propose that protein complexes propagate along the exon, thereby changing its physical behavior so as to affect 5'ss selection.

Autosomal recessive retinitis pigmentosa due to ABCA4 mutations: clinical, pathologic, and molecular characterization

PURPOSE

Autosomal recessive retinitis pigmentosa (ARRP) is a genetically heterogeneous condition characterized by progressive loss of retinal photoreceptor cells. In order to gain new insights into the pathogenesis of ARRP, we evaluated the morphological, biochemical, and gene expression changes in eyes from a human donor with ARRP due to mutations in the ABCA4 gene.

METHODS

Eyes were obtained postmortem from a donor with end-stage retinitis pigmentosa. The coding sequences of the RDS, RHO, and ABCA4 genes were screened for disease-causing mutations. Morphological changes in different regions of the retina were examined histologically, and levels of lipofuscin-associated bisretinoids were measured. Gene expression was examined in retinal/choroidal tissue using microarray analysis, and all parameters were compared to those in unaffected control donors.

RESULTS

Genetic analysis of the donor's DNA identified two mutations in the ABCA4 gene, IVS14+1G > C and Phe1440del1 cT, each on a separate allele. Morphological evaluation revealed complete loss of the outer nuclear layer, remodeling of the inner retina, loss of retinal vasculature, and regional neovascularization. The retinal pigment epithelium and choriocapillaris exhibited regional preservation. Microarray analysis revealed loss of photoreceptor cell-associated transcripts, with preservation of multiple genes expressed specifically in inner retinal neurons.

CONCLUSIONS

The persistence of transcripts expressed by inner retinal neurons suggests that despite significant plasticity that occurs during retinal degeneration, bipolar cells and ganglion cells remain at least partially differentiated. Findings from this study suggest that some forms of therapy currently under investigation may have benefit even in advanced retinal degeneration.

A brief review of molecular information theory

The idea that we could build molecular communications systems can be advanced by investigating how actual molecules from living organisms function. Information theory provides tools for such an investigation. This review describes how we can compute the average information in the DNA binding sites of any genetic control protein and how this can be extended to analyze its individual sites. A formula equivalent to Claude Shannon's channel capacity can be applied to molecular systems and used to compute the efficiency of protein binding. This efficiency is often 70% and a brief explanation for that is given. The results imply that biological systems have evolved to function at channel capacity, which means that we should be able to build molecular communications that are just as robust as our macroscopic ones.

Validation of Fanconi anemia complementation Group A assignment using molecular analysis

PURPOSE

Fanconi anemia is a genetically heterogeneous chromosomal breakage disorder exhibiting a high degree of clinical variability. Clinical diagnoses are confirmed by testing patient cells for increased sensitivity to crosslinking agents. Fanconi anemia complementation group assignment, essential for efficient molecular diagnosis of the disease, had not been validated for clinical application before this study. The purpose of this study was (1) confirmation of the accuracy of Fanconi anemia complementation group assignment to Group A (FANCA) and (2) development of a rapid mutation detection strategy that ensures the efficient capture of all FANCA mutations.

METHODS

Using fibroblasts from 29 patients, diagnosis of Fanconi anemia and assignment to complementation Group A was made through breakage analysis studies. FANCA coding and flanking sequences were analyzed using denaturing high pressure liquid chromatography, sequencing, and multiplex ligation-dependent probe amplification. Patients in which two mutations were not identified were analyzed by cDNA sequencing. Patients with no mutations were sequenced for mutations in FANCC, G, E, and F.

RESULTS

Of the 56 putative mutant alleles studied, 89% had an identifiable FANCA pathogenic mutation. Eight unique novel mutations were identified.

CONCLUSION

Complementation assignment to Group A was validated in a clinical laboratory setting using our FANCA rapid molecular testing strategy.

Can donor splice site recognition occur without the involvement of U1 snRNP?

Many disease-causing mutations affecting donor splice site recognition are reported in the literature. One of the more frequently observed nucleotide changes causing aberrant splicing are due to mutations in the donor splice site which lower the strength of base pairing with U1 snRNA (small nuclear RNA). However, recent data have highlighted the possibility of a recognition mechanism for weak donor splice sites that are at least partially U1-independent. This is important as most of the donor splice site prediction programs currently in use are based on the U1 snRNA 5′-splice site base pairing and would not pick this up. We review these mechanisms and how an up-to-date donor splice site mutation repertoire indicates the heterogeneity of the molecular mechanism. We suggest that, in clinical molecular genetics, it is important to evaluate sequence variants for aberrant splicing even in those cases where the variant is not thought to alter the U1 snRNA interaction.

Human branch point consensus sequence is yUnAy

Yeast carries a strictly conserved branch point sequence (BPS) of UACUAAC, whereas the human BPS is degenerative and is less well characterized. The human consensus BPS has never been extensively explored in vitro to date. Here, we sequenced 367 clones of lariat RT-PCR products arising from 52 introns of 20 human housekeeping genes. Among the 367 clones, a misincorporated nucleotide at the branch point was observed in 181 clones, for which we can precisely pinpoint the branch point. The branch points were comprised of 92.3% A, 3.3% C, 1.7% G and 2.8% U. Our analysis revealed that the human consensus BPS is simply yUnAy, where the underlined is the branch point at position zero and the lowercase pyrimidines ('y') are not as well conserved as the uppercase U and A. We found that the branch points are located 21-34 nucleotides upstream of the 3' end of an intron in 83% clones. We also found that the polypyrimidine tract spans 4-24 nucleotides downstream of the branch point. Our analysis demonstrates that the human BPSs are more degenerative than we have expected and that the human BPSs are likely to be recognized in combination with the polypyrimidine tract and/or the other splicing cis-elements.

In vitro and in silico analysis reveals an efficient algorithm to predict the splicing consequences of mutations at the 5' splice sites

We have found that two previously reported exonic mutations in the PINK1 and PARK7 genes affect pre-mRNA splicing. To develop an algorithm to predict underestimated splicing consequences of exonic mutations at the 5' splice site, we constructed and analyzed 31 minigenes carrying exonic splicing mutations and their derivatives. We also examined 189,249 U2-dependent 5' splice sites of the entire human genome and found that a new variable, the SD-Score, which represents a common logarithm of the frequency of a specific 5' splice site, efficiently predicts the splicing consequences of these minigenes. We also employed the information contents (R(i)) to improve the prediction accuracy. We validated our algorithm by analyzing 32 additional minigenes as well as 179 previously reported splicing mutations. The SD-Score algorithm predicted aberrant splicings in 198 of 204 sites (sensitivity = 97.1%) and normal splicings in 36 of 38 sites (specificity = 94.7%). Simulation of all possible exonic mutations at positions -3, -2 and -1 of the 189 249 sites predicts that 37.8, 88.8 and 96.8% of these mutations would affect pre-mRNA splicing, respectively. We propose that the SD-Score algorithm is a practical tool to predict splicing consequences of mutations affecting the 5' splice site.

Comparative analysis of splice site regions by information content

We have applied concepts from information theory for a comparative analysis of donor (gt) and acceptor (ag) splice site regions in the genes of five different organisms by calculating their mutual information content (relative entropy) over a selected block of nucleotides. A similar pattern that the information content decreases as the block size increases was observed for both regions in all the organisms studied. This result suggests that the information required for splicing might be contained in the consensus of ~6–8 nt at both regions. We assume from our study that even though the nucleotides are showing some degrees of conservation in the flanking regions of the splice sites, certain level of variability is still tolerated, which leads the splicing process to occur normally even if the extent of base pairing is not fully satisfied. We also suggest that this variability can be compensated by recognizing different splice sites with different spliceosomal factors.

Measurement of biological information with applications from genes to landscapes

Biological diversity is quantified for reasons ranging from primer design, to bioprospecting, and community ecology. As a common index for all levels, we suggest Shannon's (S)H, already used in information theory and biodiversity of ecological communities. Since Lewontin's first use of this index to describe human genetic variation, it has been used for variation of viruses, splice-junctions, and informativeness of pedigrees. However, until now there has been no theory to predict expected values of this index under given genetic and demographic conditions. We present a new null theory for (S)H at the genetic level, and show that this index has advantages including (i) independence of measures at each hierarchical level of organization; (ii) robust estimation of genetic exchange over a wide range of conditions; (iii) ability to incorporate information on population size; and (iv) explicit relationship to standard statistical tests. Utilization of this index in conjunction with other existing indices offers powerful insights into genetic processes. Our genetic theory is also extendible to the ecological community level, and thus can aid the comparison and integration of diversity at the genetic and community levels, including the need for measures of community diversity that incorporate the genetic differentiation between species.

Predicting severity of haemophilia A and B splicing mutations by information analysis

Bleeding symptoms and clotting activity vary among mutations that alter mRNA splicing of either the factor VIII or factor IX genes. We analyzed splicing mutations in both genes for changes in individual information (R(i), in bits) involving both donor or acceptor sites. Mutations with low or negative R(i) values (<2.4 bits) or significant changes in R(i) (DeltaR(i) > or = 7 bits) exhibited either reduced protein activity, increased clotting time and bleeding frequency and were predictive of severe disease. Thus, information analysis of splicing mutations may be useful in predicting phenotypes in hemophilia.

Automated splicing mutation analysis by information theory

Information theory-based software tools have been useful in interpreting noncoding sequence variation within functional sequence elements such as splice sites. Individual information analysis detects activated cryptic splice sites and associated splicing regulatory sites and is capable of distinguishing null from partially functional alleles. We present a server (https://splice.cmh.edu) designed to analyze splicing mutations in binding sites in either human genes, genome-mapped mRNAs, user-defined sequences, or dbSNP entries. Standard HUGO-approved gene symbols and HGVS-approved systematic mutation nomenclature (or dbSNP format) are entered via a web portal. After verifying the accuracy of input variant(s), the surrounding interval is retrieved from the human genome or user-supplied reference sequence. The server then computes the information contents (Ri) of all potential constitutive and/or regulatory splice sites in both the reference and variant sequences. Changes in information content are color-coded, tabulated, and visualized as sequence walkers, which display the binding sites with the reference sequence. The software was validated by analyzing approximately 1,300 mutations from Human Mutation as well as eight mapped SNPs from dbSNP designated as splice site variants. All of the splicing mutations and variants affected splice site strength or activated cryptic splice sites. The server also detected several missense mutations that were unexpectedly predicted to have concomitant effects on splicing or appeared to activate cryptic splicing.

Screening for inter-individual splicing differences in human GSTM4 and the discovery of a single nucleotide substitution related to the tandem skipping of two exons

The glutathione S-transferase Mu class (GSTM) genes encode phase II metabolism enzymes that are involved in the detoxification of various carcinogens and drugs. Some genetic polymorphisms in GSTM genes are related to disease phenotypes and drug-metabolism differences in the population. Polymorphisms that alter gene-splicing patterns are functionally very important because they often lead to the insertion or deletion of many amino acids. To identify inter-individual differences in the splicing pattern of the GSTM4 gene, we used reverse transcriptase polymerase chain reaction (RT-PCR) to screen cDNA from 96 human liver samples. We discovered a novel splice variant of GSTM4 that resulted from tandem skipping of exons 4 and 5. This exon-skipping event is associated with a mutation at the splice acceptor site in intron 4.

Development and refinement of pregnane X receptor (PXR) DNA binding site model using information theory: insights into PXR-mediated gene regulation

The pregnane X receptor (PXR) acts as a receptor to induce gene expression in response to structurally diverse xenobiotics through binding as a heterodimer with the 9-cis retinoic acid receptor (RXR) to enhancers in target gene promoters. We identified and estimated the affinities of novel PXR/RXR binding sites in regulated genes and additional genomic targets of PXR with an information theory-based model of the PXR/RXR binding site. Our initial PXR/RXR model, the result of the alignment of 15 previously characterized binding sites, was used to scan the promoters of known PXR target genes. Sites from these genes, with information contents of >8 bits bound by PXR/RXR in vitro, were used to revise the information weight matrix; this procedure was repeated by screening for progressively weaker binding sites. After three iterations of refinement, the model was based on 48 validated PXR/RXR binding sites and has an average information content (Rsequence) of 14.43 +/- 3.21 bits. A scan of the human genome predicted novel PXR/RXR binding sites in the promoters of UGT1A3 (19.78 bits at -8040 and 16.37 bits at -6930) and UGT1A6 (12.74 bits at -9216), both of which were identified previously as targets for PXR. These sites were subsequently demonstrated to specifically bind PXR/RXR in competition electrophoretic mobility shift assays. A strong PXR site was also predicted upstream of the CASP10 gene (18.69 bits at -7872) and was validated by binding studies and reporter assays as a PXR responsive element. This suggests that the PXR-mediated response extends beyond genes involved in drug biotransformation and transport.

Computational comparative analyses of alternative splicing regulation using full-length cDNA of various eukaryotes

We previously reported a computational approach to infer alternative splicing patterns from Mus musculus full-length cDNA clones and microarray data. Although we predicted a large number of unreported splice variants, the general mechanisms regulating alternative splicing were yet unknown. In the present study, we compared alternative exons and constitutive exons in terms of splice-site strength and frequency of potential regulatory sequences. These regulatory features were further compared among five different species: Homo sapiens, M. musculus, Arabidopsis thaliana, Oryza sativa, and Drosophila melanogaster. Solid statistical validations of our comparative analyses indicated that alternative exons have (1) weaker splice sites and (2) more potential regulatory sequences than constitutive exons. Based on our observations, we propose a combinatorial model of alternative splicing mechanisms, which suggests that alternative exons contain weak splice sites regulated alternatively by potential regulatory sequences on the exons.

Polymorphisms in the Human Monomethylarsonic Acid (MMAV) Reductase/hGSTO1 Gene and Changes in Urinary Arsenic Profiles

Large interindividual variability in urinary arsenic profiles, following chronic inorganic arsenic exposure, is well-known in humans. To understand this variability, we studied the relationship between polymorphisms in the gene for human monomethylarsonic acid (MMA(V)) reductase/hGSTO1 and the urinary arsenic profiles of individuals chronically exposed to arsenic in their drinking water. To ensure that we did not overlook rare polymorphisms, not included in the public databases, we amplified and sequenced all six exons of the gene and their flanking regions, using DNA isolated from peripheral blood samples of 75 subjects, living in the vicinity of Torreon, Mexico. Four groups, based on the levels of arsenic (9-100 microg/L) in their drinking water, were studied. We identified six novel polymorphisms and two reported previously. The novel polymorphisms were a three base pair deletion (delGGC) in the first intron; a G > C transversion, leading to a serine-to-cysteine substitution at amino acid 86; a G > T transversion and a A > T transversion in intron 5; a G > A transition resulting in glutamate-to-lysine substitution in amino acid 208; and a C > T transition producing an alanine-to-valine substitution in amino acid 236. Two subjects displayed significant differences in patterns of urinary arsenic; they had increased levels of urinary inorganic arsenic and reduced levels of methylated urinary arsenic species as compared to the rest of the study population. These two subjects had the same unique polymorphisms in hGSTO1 in that they were heterozygous for E155del and Glu208Lys. The identified SNPs may be one of the reasons for the large interindividual variability in the response of humans to chronic inorganic arsenic exposure. The findings suggest the need for further studies to identify unambiguously specific polymorphisms that may account for interindividual variability in the human response to chronic inorganic arsenic exposure.

Clinical heterogeneity and molecular findings in five Polish patients with glycerol kinase deficiency: investigation of two splice site mutations with computerized splice junction analysis and Xp21 gene-specific mRNA analysis

Five cases of glycerol kinase deficiency are presented with clinical, biochemical, and genetic results. Two had the glycerol kinase deficiency as part of an Xp21 contiguous gene deletion syndrome-complex form-and three had an isolated form of the enzyme deficiency. In these we found two splice site mutations (IVS1+4A>G, IVS9-1G>T) and one insertion (1393_1394insG). In patients with the complex form, a deletion of the DAX1, GK genes and the distal part of the DMD gene was found. A computerized study was performed to predict the effects of the splice site mutations. It showed that the IVS9-1G>T mutation substantially altered and removed the wild-type site and enhanced a cryptic site seven nucleotides downstream, and that the IVS1+4A>G diminished the strength of the wild-type donor site from strong to leaky. To verify these predictions, we developed an RT-PCR system with gene-specific primers that exclusively amplifies the Xp21 glycerol kinase gene transcript. Identification of individuals at risk is motivated by a need to avoid delay in a correct diagnosis. For reliable identification of heterozygotes for isolated glycerol kinase deficiency, knowledge of the specific mutation in the proband is required. This is easily obtained with the RT-PCR analyses developed in this study.

Information theory-based analysis of CYP2C19, CYP2D6 and CYP3A5 splicing mutations

Several mutations are known or suspected to affect mRNA splicing of CYP2C19, CYP2D6 and CYP3A5 genes; however, little experimental evidence exists to support these conclusions. The present study applies mathematical models that measure changes in information content of splice sites in these genes to demonstrate the relationship between the predicted phenotypes of these variants to the corresponding genotypes. Based on information analysis, the CYP2C19*2 variant activates a new cryptic site 40 nucleotides downstream of the natural splice site. CYP2C19*7 abolishes splicing at the exon 5 donor site. The CYP2D6*4 allele similarly inactivates splicing at the acceptor site of exon 4 and activates a new cryptic site one nucleotide downstream of the natural acceptor. CYP2D6*11 inactivates the acceptor site of exon 2. The CYP3A5*3 allele activates a new cryptic site 236 nucleotides upstream of the exon 4 natural acceptor site. CYP3A5*5 inactivates the exon 5 donor site and CYP3A5*6 strengthens a site upstream of the natural donor site, resulting in skipping of exon 7. Other previously described missense and nonsense mutations at terminal codons of exons in these genes affected splicing. CYP2D6*8 and CYP2D6*14 both decrease the strength of the exon 3 donor site, producing transcripts lacking this exon. The results of information analysis are consistent with the poor metabolizer phenotypes observed in patients with these mutations, and illustrate the potential value of these mathematical models to quantitatively evaluate the functional consequences of new mutations suspected of altering mRNA splicing.

Newfoundland rod-cone dystrophy, an early-onset retinal dystrophy, is caused by splice-junction mutations in RLBP1

Some isolated populations exhibit an increased prevalence of rare recessive diseases. The island of Newfoundland is a characteristic geographic isolate, settled by a small number of families primarily during the late 1700s and early 1800s. During our studies of this population, we identified a group of families exhibiting a retinal dystrophy reminiscent of retinitis punctata albescens but with a substantially lower age at onset and more-rapid and distinctive progression, a disorder that we termed "Newfoundland rod-cone dystrophy" (NFRCD). The size of one of these families was sufficient to allow us to perform a genomewide screen to map the NFRCD locus. We detected significant linkage to markers on the long arm of chromosome 15, in a region encompassing RLBP1, the gene encoding the cellular retinaldehyde-binding protein. Previously, mutations in RLBP1 have been associated with other retinal dystrophies, leading us to hypothesize that RLBP1 mutations might also cause NFRCD. To test this hypothesis, we sequenced all coding exons and splice junctions of RLBP1. We detected two sequence alterations, each of which is likely to be pathogenic, since each segregates with the disease and is predicted to interfere with mRNA splicing. In contrast to some previously reported RLBP1 mutations, which yield a protein that may retain some residual activity, each NFRCD mutation is likely to give rise to a null allele. This difference may account for the severe phenotype in these families and exemplifies the molecular continuum that underlies clinically distinct but genetically related entities.

Information analysis of human splice site mutations

Splice site nucleotide substitutions can be analyzed by comparing the individual information contents (Ri, bits) of the normal and variant splice junction sequences [Rogan and Schneider, 1995]. In the present study, we related splicing abnormalities to changes in Ri values of 111 previously reported splice site substitutions in 41 different genes. Mutant donor and acceptor sites have significantly less information than their normal counterparts. With one possible exception, primary mutant sites with <2.4 bits were not spliced. Sites with Ri values > or = 2.4 bits but less than the corresponding natural site usually decreased, but did not abolish splicing. Substitutions that produced small changes in Ri probably do not impair splicing and are often polymorphisms. The Ri values of activated cryptic sites were generally comparable to or greater than those of the corresponding natural splice sites. Information analysis revealed preexisting cryptic splice junctions that are used instead of the mutated natural site. Other cryptic sites were created or strengthened by sequence changes that simultaneously altered the natural site. Comparison between normal and mutant splice site Ri values distinguishes substitutions that impair splicing from those which do not, distinguishes null alleles from those that are partially functional, and detects activated cryptic splice sites.

Splice-site mutations in atherosclerosis candidate genes: relating individual information to phenotype

BACKGROUND

Nucleotide variants in several genes for lipid and methionine metabolism influence the risk of premature atherosclerosis. Ten percent of single nucleotide substitutions in these genes involve mRNA splice sites. The effects of some of these changes on splicing and on phenotypic severity are not inherently obvious.

METHODS AND RESULTS

Using an information theory-based model, we measured the individual information content (R(i), in bits) of splice sites adjacent to 289 mutations (including 31 splice-site mutations) in the atherosclerosis candidate genes APOAII, APOB, APOCII, APOE, CBS, CETP, LCAT, LIPA, LDLR, and LPL. The predictions of information analysis were then corroborated by published mRNA analyses. The R(i) values of mutant sites were consistent with either complete (n=17) or partial (n=8) inactivation of these sites. Seven mutations were predicted to activate cryptic splice sites. Predicted inactive mutant sites were associated with either "average" or "severe" dyslipidemia and commensurate reductions in protein levels or activity, whereas mutations expected to exhibit residual splicing had average or "mild" effects on lipid and protein expression.

CONCLUSIONS

Information analysis of splice-junction variants in atherosclerosis candidate genes distinguishes inactive from leaky splice sites and identifies activated cryptic sites. Predicted changes in splicing were related to phenotypic severity.

Mutations that alter RNA splicing of the human HPRT gene: a review of the spectrum

The human HPRT gene contains spans approximately 42,000 base pairs in genomic DNA, has a mRNA of approximately 900 bases and a protein coding sequence of 657 bases (initiation codon AUG to termination codon UAA). This coding sequence is distributed into 9 exons ranging from 18 (exon 5) to 184 (exon 3) base pairs. Intron sizes range from 170 (intron 7) to 13,075 (intron 1) base pairs. In a database of human HPRT mutations, 277 of 2224 (12.5%) mutations result in alterations in splicing of the mRNA as analyzed by both reverse transcriptase mediated production of a cDNA followed by PCR amplification and cDNA sequencing and by genomic DNA PCR amplification and sequencing. Mutations have been found in all eight 5' (donor) and 3' (acceptor) splice sequences. Mutations in the 5' splice sequences of introns 1 and 5 result in intron inclusion in the cDNA due to the use of cryptic donor splice sequences within the introns; mutations in the other six 5' sites result in simple exon exclusion. Mutations in the 3' splice sequences of introns 1, 3, 7 and 8 result in partial exon exclusion due to the use of cryptic acceptor splice sequences within the exons; mutations in the other four 3' sites result in simple exon exclusion. A base substitution in exon 3 (209G-->T) creates a new 5' (donor) splice site which results in the exclusion of 110 bases of exon 3 from the cDNA. Two base substitutions in intron 8 (IVS8-16G-->A and IVS8-3T-->G) result in the inclusion of intron 8 sequences in the cDNA due to the creation of new 3' (acceptor) splice sites. Base substitution within exons 1, 3, 4, 6 and 8 also result in splice alterations in cDNA. Those in exons 1 and 6 are at the 3' end of the exon and may directly affect splicing. Those within exons 3 and 4 may be the result of the creation of nonsense codons, while those in exon 8 cannot be explained by this mechanism. Lastly, many mutations that affect splicing of the HPRT mRNA have pleiotropic effects in that multiple cDNA products are found.

Sequence walkers: a graphical method to display how binding proteins interact with DNA or RNA sequences

A graphical method is presented for displaying how binding proteins and other macromolecules interact with individual bases of nucleotide sequences. Characters representing the sequence are either oriented normally and placed above a line indicating favorable contact, or upside-down and placed below the line indicating unfavorable contact. The positive or negative height of each letter shows the contribution of that base to the average sequence conservation of the binding site, as represented by a sequence logo. These sequence 'walkers' can be stepped along raw sequence data to visually search for binding sites. Many walkers, for the same or different proteins, can be simultaneously placed next to a sequence to create a quantitative map of a complex genetic region. One can alter the sequence to quantitatively engineer binding sites. Database anomalies can be visualized by placing a walker at the recorded positions of a binding molecule and by comparing this to locations found by scanning the nearby sequences. The sequence can also be altered to predict whether a change is a polymorphism or a mutation for the recognizer being modeled.

Detection and characterization of a novel splice mutation in the LDL receptor intron 12 resulting in two different mutant mRNA variants

Using a simple, standardized denaturing gradient gel electrophoresis (DGGE) based mutation screening technique, a novel G-to-A mutation in the last base of the intron 12 splice acceptor site of the LDL receptor gene was found in 2 Danish families with familial hypercholesterolemia (FH). The mutation is shown to result in 2 mRNA splice variants, both leading to truncated LDLR proteins, containing only the first 594 of the normal 839 amino acids. In one of the FH-families harbouring the mutation, a striking difference in the clinical picture amongst biochemically diagnosed FH patients was clarified when genetic analysis showed that 2 hypercholesterolemic family members, who despite advanced age had no atherosclerotic disease, had not inherited the family LDLR mutation. DGGE analyses of the LDLR exons, LDLR promoter, and apolipoprotein B codon 3456-3553 as well as Southern blotting of the LDLR gene were without signs of other mutations in the non-atherosclerotic hypercholesterolemics of the family. Availability of the clinically applicable mutation screening assay for FH may thus aid in defining reasons for phenotypic differences in FH families and potentially supply information allowing a more differentiated therapeutic approach to individual members of FH families.