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

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.

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.

A novel splice-site mutation in angiotensin I-converting enzyme (ACE) gene, c.3691+1G>A (IVS25+1G>A), causes a dramatic increase in circulating ACE through deletion of the transmembrane anchor

BACKGROUND

Angiotensin-converting enzyme (ACE) (EC 4.15.1) metabolizes many biologically active peptides and plays a key role in blood pressure regulation and vascular remodeling. Elevated ACE levels are associated with different cardiovascular and respiratory diseases.

METHODS AND RESULTS

Two Belgian families with a 8-16-fold increase in blood ACE level were incidentally identified. A novel heterozygous splice site mutation of intron 25 - IVS25+1G>A (c.3691+1G>A) - cosegregating with elevated plasma ACE was identified in both pedigrees. Messenger RNA analysis revealed that the mutation led to the retention of intron 25 and Premature Termination Codon generation. Subjects harboring the mutation were mostly normotensive, had no left ventricular hypertrophy or cardiovascular disease. The levels of renin-angiotensin-aldosterone system components in the mutated cases and wild-type controls were similar, both at baseline and after 50 mg captopril. Compared with non-affected members, quantification of ACE surface expression and shedding using flow cytometry assay of dendritic cells derived from peripheral blood monocytes of affected members, demonstrated a 50% decrease and 3-fold increase, respectively. Together with a dramatic increase in circulating ACE levels, these findings argue in favor of deletion of transmembrane anchor, leading to direct secretion of ACE out of cells.

CONCLUSIONS

We describe a novel mutation of the ACE gene associated with a major familial elevation of circulating ACE, without evidence of activation of the renin-angiotensin system, target organ damage or cardiovascular complications. These data are consistent with the hypothesis that membrane-bound ACE, rather than circulating ACE, is responsible for Angiotensin II generation and its cardiovascular consequences.

Mutations in the LDL receptor gene in four Chinese homozygous familial hypercholesterolemia phenotype patients

BACKGROUND AND AIMS

Familial hypercholesterolemia (FH) is an autosomal dominant disorder of lipoprotein metabolism caused by mutations in the low-density lipoprotein receptor (LDL-R) gene, leading to elevated levels of cholesterol and an increased risk of coronary heart disease. In this article, from four homozygous FH phenotype probands we identified disease causing mutations and analyzed the relationship between genotype and phenotype.

METHODS AND RESULTS

DNA sequencing identified five LDL-R point mutations in four unrelated families. We found a novel homozygous mutation (C210R), a homozygous mutation at W462X, a compound heterozygous mutation of C122Y and T383I, and a G>A intron 3 splice site homozygous mutation. The functional alteration caused by the novel C210R mutation was confirmed by FACS analysis. Four probands have high low-density lipoprotein cholesterol (LDL-C) levels, ranging from 14.65 to 27.66 mmol/L. Their heterozygous parents had relatively low levels. B-mode ultrasound supplemented by Doppler was used to examine aortic/mitral valve structural alterations and carotid intima-media thickness (ITM) in all probands. The ITM values were between 1.2 and 2.3mm, much higher than the normal value of <0.8mm.

CONCLUSION

Our data demonstrated that all the probands were associated with severe hypercholesterolemia, thick carotid IMT and a low CFVR (coronary flow velocity reserve) value. The novel mutation (C120Y) is a disease causing mutation.

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.

Exonization of AluYa5 in the human ACE gene requires mutations in both 3' and 5' splice sites and is facilitated by a conserved splicing enhancer

Ancient Alu elements have been shown to be included in mature transcripts by point mutations that improve their 5′ or 3′ splice sites. We have examined requirements for exonization of a younger, disease-associated AluYa5 in intron 16 of the human ACE gene. A single G>C transversion in position −3 of the new Alu exon was insufficient for Alu exonization and a significant inclusion in mRNA was only observed when improving several potential splice donor sites in the presence of 3′ CAG. Since complete Alu exonization was not achieved by optimizing traditional splicing signals, including the branch site, we tested whether auxiliary elements in AluYa5 were required for constitutive inclusion. Exonization was promoted by a SELEX-predicted heptamer in Alu consensus sequence 222–228 and point mutations in highly conserved nucleotides of this heptamer decreased Alu inclusion. In addition, we show that Alu exonization was facilitated by a subset of serine/arginine-rich (SR) proteins through activation of the optimized 3′ splice site. Finally, the haplotype- and allele-specific ACE minigenes generated similar splicing patterns in both ACE-expressing and non-expressing cells, suggesting that previously reported allelic association with plasma ACE activity and cardiovascular disease is not attributable to differential splicing of introns 16 and 17.

Splice-site contribution in alternative splicing ofPLP1 andDM20: molecular studies in oligodendrocytes

Mutations in the proteolipid protein 1 (PLP1) gene cause the X-linked dysmyelinating diseases Pelizaeus-Merzbacher disease (PMD) and spastic paraplegia 2 (SPG2). We examined the severity of the following mutations that were suspected of affecting levels of PLP1 and DM20 RNA, the alternatively spliced products of PLP1: c.453G>A, c.453G>T, c.453G>C, c.453+2T>C, c.453+4A>G, c.347C>A, and c.453+28_+46del (the old nomenclature did not include the methionine codon: G450A, G450T, G450C, IVS3+2T>C, IVS3+4A>G, C344A, and IVS3+28-+46del). These mutations were evaluated by information theory-based analysis and compared with mRNA expression of the alternatively spliced products. The results are discussed relative to the clinical severity of disease. We conclude that the observed PLP1 and DM20 splicing patterns correlated well with predictions of information theory-based analysis, and that the relative strength of the PLP1 and DM20 donor splice sites plays an important role in PLP1 alternative splicing.

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.

Genetic determinants of arterial calcification associated with atherosclerosis

Increasing research interest has focused on arterial calcification in the setting of atherosclerosis. Many features of atherosclerosis-related calcification provide useful clinical information. For example, calcium mineral deposits frequently form in atherosclerotic plaque, and intimal arterial calcification can be used as a surrogate marker for atherosclerosis; also, calcium deposits are readily and noninvasively quantified, which is useful because greater amounts of coronary calcification predict a higher risk of myocardial infarction and death. Several mechanisms leading to calcification associated with atherosclerosis have been proposed; however, no direct testing of proposed mechanisms has yet been reported. Studies in genetically altered animals and in humans have shed light on potential genetic determinants, which in turn could form the basis for a more comprehensive understanding of the factors affecting calcification within plaque and the associated pathobiologic implications. We review proposed molecular and cellular mechanisms of atherosclerosis-associated arterial calcification, summarize genetic influences, and suggest areas in which further investigation is needed. Understanding the molecular and genetic determinants of specific structural plaque components such as calcification can provide a solid foundation for the development of novel therapeutic approaches to favorably alter plaque structure and minimize vulnerability to arterial rupture.

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.

Atypical 5' splice sites cause CFTR exon 9 to be vulnerable to skipping

The molecular basis of the skipping of constitutive exons in many messenger RNAs is not fully understood. A well-studied example is exon 9 of the human cystic fibrosis transmembrane conductance regulator gene (CFTR), in which an abbreviated polypyrimidine tract between the branch point A and the 3' splice site is associated with increased exon skipping and disease. However, many exons, both in CFTR and in other genes and have short polypyrimidine tracts in their 3' splice sites, yet they are not skipped. Inspection of the 5' splice sites immediately up- and downstream of exon 9 revealed deviations from consensus sequence, so we hypothesized that this exon may be inherently vulnerable to skipping. To test this idea, we constructed a CFTR minigene and replicated exon 9 skipping associated with the length of the polypyrimidine tract upstream of exon 9. We then mutated the flanking 5' splice sites and determined the effect on exon skipping. Conversion of the upstream 5' splice site to consensus by replacing a pyrimidine at position +3 with a purine resulted in increased exon skipping. In contrast, conversion of the downstream 5' splice site to consensus by insertion of an adenine at position +4 resulted in a substantial reduction in exon 9 skipping, regardless of whether the upstream 5' splice site was consensus or not. These results suggested that the native downstream 5' splice site plays an important role in CFTR exon 9 skipping, a hypothesis that was supported by data from sheep and mouse genomes. Although CFTR exon 9 in sheep is preceded by a long polypyrimidine tract (Y(14)), it skips exon 9 in vivo and has a nonconsensus downstream 5' splice site identical to that in humans. On the other hand, CFTR exon 9 in mice is preceded by a short polypyrimidine tract (Y(5)) but is not skipped in vivo. Its downstream 5' splice site differs from that in humans by a 2-nt insertion, which, when introduced into the human CFTR minigene, abolished exon 9 skipping. Taken together, these observations place renewed emphasis on deviations at 5' splice sites in nucleotides other than the invariant GT, particularly when such changes are found in conjunction with other altered splicing sequences, such as a shortened polypyrimidine tract. Thus, careful inspection of entire 5' splice sites may identify constitutive exons that are vulnerable to skipping.

Characterization of the effects of mutations in the putative branchpoint sequence of intron 4 on the splicing within the human lecithin:cholesterol acyltransferase gene

We have previously identified a point mutation (intervening sequence (IVS) 4: T --> C) in the branchpoint consensus sequence of intron 4 of the lecithin:cholesterol acyltransferase (LCAT) gene in patients with fish-eye disease. To investigate the possible mechanisms responsible for the defective splicing, we made a series of mutations in the branchpoint sequence and expressed these mutants in HEK-293 cells followed by the analysis of pre-mRNA splicing using reverse transcriptase-polymerase chain reaction as well as LCAT activity assay. The results reveal that 1) the mutation of the branchpoint adenosine to any other nucleotide completely abolishes splicing; 2) the insertion of a normal branch site into the intronic sequence of the natural (IVS4-22c) or the branchpoint (IVS4-20t) mutant completely restores splicing; 3) the natural mutation can be partially rescued by making a single nucleotide change (G --> A) within the branchpoint consensus sequence; and 4) other single base changes, particularly around the branchpoint adenosine residue, significantly decrease the efficiency of splicing and thus enzyme activity. Surprisingly, the nucleotide transversion at the last position of the branchpoint sequence (i.e. IVS4-25a or -25g) results in a 2.7-fold increase in splicing efficiency. Therefore, these observations clearly establish the functional significance of the branchpoint sequence of intron 4 for the splicing of the human LCAT mRNA precursors.