Information analysis of human splice site mutations

Functional assessment of somatic STK11 variants identified in primary human non-small cell lung cancers

Serine/Threonine Kinase 11 (STK11) encodes an important tumor suppressor that is frequently mutated in lung adenocarcinoma. Clinical studies have shown that mutations in STK11 resulting in loss of function correlate with resistance to anti-PD-1 monoclonal antibody therapy in KRAS-driven non-small cell lung cancer (NSCLC), but the molecular mechanisms responsible remain unclear. Despite this uncertainty, STK11 functional status is emerging as a reliable biomarker for predicting non-response to anti-PD-1 therapy in NSCLC patients. The clinical utility of this biomarker ultimately depends upon accurate classification of STK11 variants. For nonsense variants occurring early in the STK11 coding region, this assessment is straightforward. However, rigorously demonstrating the functional impact of missense variants remains an unmet challenge. Here we present data characterizing four STK11 splice-site variants by analyzing tumor mRNA, and 28 STK11 missense variants using an in vitro kinase assay combined with a cell-based p53-dependent luciferase reporter assay. The variants we report were identified in primary human NSCLC biopsies in collaboration with the University of Vermont Genomic Medicine group. Additionally, we compare our experimental results with data from 22 in silico predictive algorithms. Our work highlights the power, utility and necessity of functional variant assessment and will aid STK11 variant curation, provide a platform to assess novel STK11 variants and help guide anti-PD-1 therapy utilization in KRAS-driven NSCLCs.

Interpretable prioritization of splice variants in diagnostic next-generation sequencing

A critical challenge in genetic diagnostics is the computational assessment of candidate splice variants, specifically the interpretation of nucleotide changes located outside of the highly conserved dinucleotide sequences at the 5' and 3' ends of introns. To address this gap, we developed the Super Quick Information-content Random-forest Learning of Splice variants (SQUIRLS) algorithm. SQUIRLS generates a small set of interpretable features for machine learning by calculating the information-content of wild-type and variant sequences of canonical and cryptic splice sites, assessing changes in candidate splicing regulatory sequences, and incorporating characteristics of the sequence such as exon length, disruptions of the AG exclusion zone, and conservation. We curated a comprehensive collection of disease-associated splice-altering variants at positions outside of the highly conserved AG/GT dinucleotides at the termini of introns. SQUIRLS trains two random-forest classifiers for the donor and for the acceptor and combines their outputs by logistic regression to yield a final score. We show that SQUIRLS transcends previous state-of-the-art accuracy in classifying splice variants as assessed by rank analysis in simulated exomes, and is significantly faster than competing methods. SQUIRLS provides tabular output files for incorporation into diagnostic pipelines for exome and genome analysis, as well as visualizations that contextualize predicted effects of variants on splicing to make it easier to interpret splice variants in diagnostic settings.

Mutations primarily alter the inclusion of alternatively spliced exons

Genetic analyses and systematic mutagenesis have revealed that synonymous, non-synonymous and intronic mutations frequently alter the inclusion levels of alternatively spliced exons, consistent with the concept that altered splicing might be a common mechanism by which mutations cause disease. However, most exons expressed in any cell are highly-included in mature mRNAs. Here, by performing deep mutagenesis of highly-included exons and by analysing the association between genome sequence variation and exon inclusion across the transcriptome, we report that mutations only very rarely alter the inclusion of highly-included exons. This is true for both exonic and intronic mutations as well as for perturbations in trans. Therefore, mutations that affect splicing are not evenly distributed across primary transcripts but are focussed in and around alternatively spliced exons with intermediate inclusion levels. These results provide a resource for prioritising synonymous and other variants as disease-causing mutations.

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.

Functional characterization of biallelic RTTN variants identified in an infant with microcephaly, simplified gyral pattern, pontocerebellar hypoplasia, and seizures

BACKGROUND

Biallelic deleterious variants in RTTN, which encodes rotatin, are associated with primary microcephaly, polymicrogyria, seizures, intellectual disability, and primordial dwarfism in human infants.

METHODS AND RESULTS

We performed exome sequencing of an infant with primary microcephaly, pontocerebellar hypoplasia, and intractable seizures and his healthy, unrelated parents. We cultured the infant's fibroblasts to determine primary ciliary phenotype.

RESULTS

We identified biallelic variants in RTTN in the affected infant: a novel missense variant and a rare, intronic variant that results in aberrant transcript splicing. Cultured fibroblasts from the infant demonstrated reduced length and number of primary cilia.

CONCLUSION

Biallelic variants in RTTN cause primary microcephaly in infants. Functional characterization of primary cilia length and number can be used to determine pathogenicity of RTTN variants.

Discovery and validation of information theory-based transcription factor and cofactor binding site motifs

Data from ChIP-seq experiments can derive the genome-wide binding specificities of transcription factors (TFs) and other regulatory proteins. We analyzed 765 ENCODE ChIP-seq peak datasets of 207 human TFs with a novel motif discovery pipeline based on recursive, thresholded entropy minimization. This approach, while obviating the need to compensate for skewed nucleotide composition, distinguishes true binding motifs from noise, quantifies the strengths of individual binding sites based on computed affinity and detects adjacent cofactor binding sites that coordinate with the targets of primary, immunoprecipitated TFs. We obtained contiguous and bipartite information theory-based position weight matrices (iPWMs) for 93 sequence-specific TFs, discovered 23 cofactor motifs for 127 TFs and revealed six high-confidence novel motifs. The reliability and accuracy of these iPWMs were determined via four independent validation methods, including the detection of experimentally proven binding sites, explanation of effects of characterized SNPs, comparison with previously published motifs and statistical analyses. We also predict previously unreported TF coregulatory interactions (e.g. TF complexes). These iPWMs constitute a powerful tool for predicting the effects of sequence variants in known binding sites, performing mutation analysis on regulatory SNPs and predicting previously unrecognized binding sites and target genes.

Prevalence and spectrum of germline rare variants in BRCA1/2 and PALB2 among breast cancer cases in Sarawak, Malaysia

PURPOSE

To characterize the spectrum of germline mutations in BRCA1, BRCA2, and PALB2 in population-based unselected breast cancer cases in an Asian population.

METHODS

Germline DNA from 467 breast cancer patients in Sarawak General Hospital, Malaysia, where 93% of the breast cancer patients in Sarawak are treated, was sequenced for the entire coding region of BRCA1; BRCA2; PALB2; Exons 6, 7, and 8 of TP53; and Exons 7 and 8 of PTEN. Pathogenic variants included known pathogenic variants in ClinVar, loss of function variants, and variants that disrupt splice site.

RESULTS

We found 27 pathogenic variants (11 BRCA1, 10 BRCA2, 4 PALB2, and 2 TP53) in 34 patients, which gave a prevalence of germline mutations of 2.8, 3.23, and 0.86% for BRCA1, BRCA2, and PALB2, respectively. Compared to mutation non-carriers, BRCA1 mutation carriers were more likely to have an earlier age at onset, triple-negative subtype, and lower body mass index, whereas BRCA2 mutation carriers were more likely to have a positive family history. Mutation carrier cases had worse survival compared to non-carriers; however, the association was mostly driven by stage and tumor subtype. We also identified 19 variants of unknown significance, and some of them were predicted to alter splicing or transcription factor binding sites.

CONCLUSION

Our data provide insight into the genetics of breast cancer in this understudied group and suggest the need for modifying genetic testing guidelines for this population with a much younger age at diagnosis and more limited resources compared with Caucasian populations.

A unified analytic framework for prioritization of non-coding variants of uncertain significance in heritable breast and ovarian cancer

BACKGROUND

Sequencing of both healthy and disease singletons yields many novel and low frequency variants of uncertain significance (VUS). Complete gene and genome sequencing by next generation sequencing (NGS) significantly increases the number of VUS detected. While prior studies have emphasized protein coding variants, non-coding sequence variants have also been proven to significantly contribute to high penetrance disorders, such as hereditary breast and ovarian cancer (HBOC). We present a strategy for analyzing different functional classes of non-coding variants based on information theory (IT) and prioritizing patients with large intragenic deletions.

METHODS

We captured and enriched for coding and non-coding variants in genes known to harbor mutations that increase HBOC risk. Custom oligonucleotide baits spanning the complete coding, non-coding, and intergenic regions 10 kb up- and downstream of ATM, BRCA1, BRCA2, CDH1, CHEK2, PALB2, and TP53 were synthesized for solution hybridization enrichment. Unique and divergent repetitive sequences were sequenced in 102 high-risk, anonymized patients without identified mutations in BRCA1/2. Aside from protein coding and copy number changes, IT-based sequence analysis was used to identify and prioritize pathogenic non-coding variants that occurred within sequence elements predicted to be recognized by proteins or protein complexes involved in mRNA splicing, transcription, and untranslated region (UTR) binding and structure. This approach was supplemented by in silico and laboratory analysis of UTR structure.

RESULTS

15,311 unique variants were identified, of which 245 occurred in coding regions. With the unified IT-framework, 132 variants were identified and 87 functionally significant VUS were further prioritized. An intragenic 32.1 kb interval in BRCA2 that was likely hemizygous was detected in one patient. We also identified 4 stop-gain variants and 3 reading-frame altering exonic insertions/deletions (indels).

CONCLUSIONS

We have presented a strategy for complete gene sequence analysis followed by a unified framework for interpreting non-coding variants that may affect gene expression. This approach distills large numbers of variants detected by NGS to a limited set of variants prioritized as potential deleterious changes.

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.

Founder Effect of a c.828+3A&gt;T Splice Site Mutation in Peripherin 2 (PRPH2) Causing Autosomal Dominant Retinal Dystrophies

IMPORTANCE

Screening for splice site mutation c.828+3A>T in the peripherin 2 (PRPH2) gene should be a high priority in families with highly variable retinal dystrophies. The correction of missplicing is a potential therapeutic target.

OBJECTIVE

To determine the prevalence, genetic origin, and molecular mechanism of a donor c.828+3A>T mutation in the PRPH2 (peripherin 2, retinal degeneration slow) gene in individuals with retinal dystrophies.

DESIGN, SETTING, AND PARTICIPANTS

Case-control study that took place at the University of Texas Health Science Center, the University of Iowa, and the Retina Foundation of the Southwest, from January 1, 1987, to August 1, 2014, including affected individuals from 200 families with a diagnosis of autosomal dominant retinitis pigmentosa, 35 families with unspecified macular dystrophies, and 116 families with pattern dystrophy. Participants were screened for the c.828+3A>T mutation by restriction-enzyme digest, single-strand conformational polymorphism screening, or bidirectional sequencing. Haplotypes of polymorphic markers flanking the PRPH2 locus and sequence variants within the gene were determined by denaturing gel electrophoresis or automated capillary-based cycle sequencing. The effect of the splice site mutation on the PRPH2 transcript was analyzed using NetGene2, a splice prediction program and by the reverse transcription polymerase chain reaction of illegitimate transcripts from peripheral white blood cells.

MAIN OUTCOMES AND MEASURES

Results of testing for splice site mutation, haplotypes, and alternate transcripts.

RESULTS

The PRPH2 mutation was found in 97 individuals of 19 independently ascertained families with a clinical diagnosis of retinitis pigmentosa, macular dystrophy, and/or pattern dystrophy. All affected individuals also shared a rare haplotype of approximately 644 kilobase pairs containing the c.828+3A>T mutation, which extends from the short tandem repeat polymorphism D6S282 to c.1013G>A (rs434102, a single-nucleotide polymorphism) in exon 3 of PRPH2, suggesting this mutation is from a common ancestor and is a founder mutation. It has a prevalence of 2% in families diagnosed as having autosomal dominant retinitis pigmentosa and 10% in families with variable clinical diagnosis of pattern, macular, and retinal dystrophies. Individuals with the c.828+3A>T mutation expressed a PRPH2 transcript not found in control participants and that was consistent with abnormal splicing.

CONCLUSIONS AND RELEVANCE

The PRPH2 c.828+3A>T splice site mutation is a frequent cause of inherited retinal dystrophies and is owing to the founder effect. The likely cause of disease is the missplicing of the PRPH2 message that results in a truncated protein product. Identifying the genetic etiology assists in more accurate management and possible future therapeutic options.

Screening of Transforming Growth Factor Beta 3 and Jagged2 Genes in the Malay Population with Nonsyndromic Cleft Lip with or without Cleft Palate

Objective

To determine the prevalence of mutations in transforming growth factor beta 3 (TGFβ3) and Jagged2 genes and their association with nonsyndromic cleft lip with or without cleft palate (CL±P) patients.

Design

Cross-sectional study on nonsyndromic CL±P and noncleft patients.

Setting

Reconstructive clinic and outpatient dental clinic, Hospital Universiti Sains Malaysia.

Patients

Blood samples of 96 nonsyndromic CL±P and 96 noncleft subjects.

Main Outcome Measure

Prevalence and association of mutations in TGFβ3 and Jagged2 genes with nonsyndromic CL±P.

Results

Most of the nonsyndromic CL±P patients (53.1%) had left unilateral CLP. There were slightly more females (56.6%) compared with males. The prevalence of the mutations in the TGFβ3 gene was 17.7% (95% confidence interval [CI]: 9.5, 24.5) and in the Jagged2 gene was 12.5% (95% CI: 5.5, 18.5), which was higher compared with the noncleft group. For the TGFβ3 gene, there was no mutation in the coding region in either of the groups. All variants were single nucleotide polymorphisms located within the intronic flanking region. Two variants were identified (g.15812T>G and g.15966A>G) in both nonsyndromic CL±P and noncleft patients. However, the association was not significant ( P > .05). Three variants (g.19779C>T, g.19547G>A, and g.19712C>T) were identified in the Jagged2 gene among nonsyndromic CL±P and noncleft patients. Only g.19712C>T showed a significant association with nonsyndromic CL±P patients ( P = .039).

Conclusion

g.19712C>T might play a crucial role in the development of cleft lip and palate. To the best of our knowledge, this is the first report of the mutation found within intron 13 of the Jagged2 gene among nonsyndromic CL±P Malay patients.

Twenty Years of Delila and Molecular Information Theory: The Altenberg-Austin Workshop in Theoretical Biology Biological Information, Beyond Metaphor: Causality, Explanation, and Unification Altenberg, Austria, 11-14 July 2002

A brief personal history is given about how information theory can be applied to binding sites of genetic control molecules on nucleic acids. The primary example used is ribosome binding sites in Escherichia coli. Once the sites are aligned, the information needed to describe the sites can be computed using Claude Shannon's method. This is displayed by a computer graphic called a sequence logo. The logo represents an average binding site, and the mathematics easily allows one to determine the components of this average. That is, given a set of binding sites, the information for individual binding sites can also be computed. One can go further and predict the information of sites that are not in the original data set. Information theory also allows one to model the flexibility of ribosome binding sites, and this led us to a simple model for ribosome translational initiation in which the molecular components fit together only when the ribosome is at a good ribosome binding site. Since information theory is general, the same mathematics applies to human splice junctions, where we can predict the effect of sequence changes that cause human genetic diseases and cancer. The second example given is the Pribnow 'box' which, when viewed by the information theory method, reveals a mechanism for initiation of both transcription and DNA replication. Replication, transcription, splicing, and translation into protein represent the central dogma, so these examples show how molecular information theory is contributing to our knowledge of basic biology.

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.

PMD patient mutations reveal a long-distance intronic interaction that regulates PLP1/DM20 alternative splicing

Alternative splicing of the proteolipid protein 1 gene (PLP1) produces two forms, PLP1 and DM20, due to alternative use of 5' splice sites with the same acceptor site in intron 3. The PLP1 form predominates in central nervous system RNA. Mutations that reduce the ratio of PLP1 to DM20, whether mutant or normal protein is formed, result in the X-linked leukodystrophy Pelizaeus-Merzbacher disease (PMD). We investigated the ability of sequences throughout PLP1 intron 3 to regulate alternative splicing using a splicing minigene construct transfected into the oligodendrocyte cell line, Oli-neu. Our data reveal that the alternative splice of PLP1 is regulated by a long-distance interaction between two highly conserved elements that are separated by 581 bases within the 1071-base intron 3. Further, our data suggest that a base-pairing secondary structure forms between these two elements, and we demonstrate that mutations of either element designed to destabilize the secondary structure decreased the PLP1/DM20 ratio, while swap mutations designed to restore the structure brought the PLP1/DM20 ratio to near normal levels. Sequence analysis of intron 3 in families with clinical symptoms of PMD who did not have coding-region mutations revealed mutations that segregated with disease in three families. We showed that these patient mutations, which potentially destabilize the secondary structure, also reduced the PLP1/DM20 ratio. This is the first report of patient mutations causing disease by disruption of a long-distance intronic interaction controlling alternative splicing. This finding has important implications for molecular diagnostics of PMD.

Genetic testing in Tunisian families with heritable retinoblastoma using a low cost approach permits accurate risk prediction in relatives and reveals incomplete penetrance in adults

Heritable retinoblastoma is caused by oncogenic mutations in the RB1 tumor suppressor gene. Identification of these mutations in patients is important for genetic counseling and clinical management of relatives at risk. In order to lower analytical efforts, we designed a stepwise mutation detection strategy that was adapted to the spectrum of oncogenic RB1 gene mutations. We applied this strategy on 20 unrelated patients with familial and/or de novo bilateral retinoblastoma from Tunisia. In 19 (95%) patients, we detected oncogenic mutations including base substitutions, small length mutations, and large deletions. Further analyses on the origin of the mutations showed mutational mosaicism in one unilaterally affected father of a bilateral proband and incomplete penetrance in two mothers. In a large family with several retinoblastoma patients, the mutation identified in the index patient was also detected in several non-penetrant relatives. RNA analyses showed that this mutation results in an in-frame loss of exon 9. In summary, our strategy can serve as a model for RB1 mutation identification with high analytical sensitivity. Our results point out that genetic testing is needed to reveal or exclude incomplete penetrance specifically in parents of patients with sporadic disease.

A rare sequence variant in intron 1 of THAP1 is associated with primary dystonia

Although coding variants in THAP1 have been causally associated with primary dystonia, the contribution of noncoding variants remains uncertain. Herein, we examine a previously identified Intron 1 variant (c.71+9C>A, rs200209986). Among 1672 subjects with mainly adult-onset primary dystonia, 12 harbored the variant in contrast to 1/1574 controls (P < 0.01). Dystonia classification included cervical dystonia (N = 3), laryngeal dystonia (adductor subtype, N = 3), jaw-opening oromandibular dystonia (N = 1), blepharospasm (N = 2), and unclassified (N = 3). Age of dystonia onset ranged from 25 to 69 years (mean = 54 years). In comparison to controls with no identified THAP1 sequence variants, the c.71+9C>A variant was associated with an elevated ratio of Isoform 1 (NM_018105) to Isoform 2 (NM_199003) in leukocytes. In silico and minigene analyses indicated that c.71+9C>A alters THAP1 splicing. Lymphoblastoid cells harboring the c.71+9C>A variant showed extensive apoptosis with relatively fewer cells in the G2 phase of the cell cycle. Differentially expressed genes from lymphoblastoid cells revealed that the c.71+9C>A variant exerts effects on DNA synthesis, cell growth and proliferation, cell survival, and cytotoxicity. In aggregate, these data indicate that THAP1 c.71+9C>A is a risk factor for adult-onset primary dystonia.

Validation of predicted mRNA splicing mutations using high-throughput transcriptome data

Interpretation of variants present in complete genomes or exomes reveals numerous sequence changes, only a fraction of which are likely to be pathogenic. Mutations have been traditionally inferred from allele frequencies and inheritance patterns in such data. Variants predicted to alter mRNA splicing can be validated by manual inspection of transcriptome sequencing data, however this approach is intractable for large datasets. These abnormal mRNA splicing patterns are characterized by reads demonstrating either exon skipping, cryptic splice site use, and high levels of intron inclusion, or combinations of these properties. We present, Veridical, an in silico method for the automatic validation of DNA sequencing variants that alter mRNA splicing. Veridical performs statistically valid comparisons of the normalized read counts of abnormal RNA species in mutant versus non-mutant tissues. This leverages large numbers of control samples to corroborate the consequences of predicted splicing variants in complete genomes and exomes.

Validation of predicted mRNA splicing mutations using high-throughput transcriptome data

Interpretation of variants present in complete genomes or exomes reveals numerous sequence changes, only a fraction of which are likely to be pathogenic. Mutations have been traditionally inferred from allele frequencies and inheritance patterns in such data. Variants predicted to alter mRNA splicing can be validated by manual inspection of transcriptome sequencing data, however this approach is intractable for large datasets. These abnormal mRNA splicing patterns are characterized by reads demonstrating either exon skipping, cryptic splice site use, and high levels of intron inclusion, or combinations of these properties. We present, Veridical, an in silico method for the automatic validation of DNA sequencing variants that alter mRNA splicing. Veridical performs statistically valid comparisons of the normalized read counts of abnormal RNA species in mutant versus non-mutant tissues. This leverages large numbers of control samples to corroborate the consequences of predicted splicing variants in complete genomes and exomes.

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.

Global remodeling of nucleosome positions in C. elegans

Background

Eukaryotic chromatin architecture is affected by intrinsic histone-DNA sequence preferences, steric exclusion between nucleosome particles, formation of higher-order structures, and in vivo activity of chromatin remodeling enzymes.

Results

To disentangle sequence-dependent nucleosome positioning from the other factors, we have created two high-throughput maps of nucleosomes assembled in vitro on genomic DNA from the nematode worm Caenorhabditis elegans. A comparison of in vitro nucleosome positions with those observed in a mixed-stage, mixed-tissue population of C. elegans cells reveals that in vivo sequence preferences are modified on the genomic scale. Indeed, G/C dinucleotides are predicted to be most favorable for nucleosome formation in vitro but not in vivo. Nucleosome sequence read coverage in vivo is distinctly lower in chromosome arms than in central regions; the observed changes in apparent nucleosome sequence specificity, likely due to genome-wide chromatin remodeler activity, contribute to the formation of these megabase-scale chromatin domains. We also observe that the majority of well-positioned in vivo nucleosomes do not occupy thermodynamically favorable sequences observed in vitro. Finally, we find that exons are intrinsically more amenable to nucleosome formation compared to introns. Nucleosome occupancy of introns and exons consistently increases with G/C content in vitro but not in vivo, in agreement with our observation that G/C dinucleotide enrichment does not strongly promote in vivo nucleosome formation.

Conclusions

Our findings highlight the importance of both sequence specificity and active nucleosome repositioning in creating large-scale chromatin domains, and the antagonistic roles of intrinsic sequence preferences and chromatin remodelers in C. elegans.

Sequence read data has been deposited into Sequence Read Archive (http://www.ncbi.nlm.nih.gov/sra; accession number SRA050182). Additional data, software and computational predictions are available on the Nucleosome Explorer website (http://nucleosome.rutgers.edu).

Interpretation, stratification and evidence for sequence variants affecting mRNA splicing in complete human genome sequences

Information theory-based methods have been shown to be sensitive and specific for predicting and quantifying the effects of non-coding mutations in Mendelian diseases. We present the Shannon pipeline software for genome-scale mutation analysis and provide evidence that the software predicts variants affecting mRNA splicing. Individual information contents (in bits) of reference and variant splice sites are compared and significant differences are annotated and prioritized. The software has been implemented for CLC-Bio Genomics platform. Annotation indicates the context of novel mutations as well as common and rare SNPs with splicing effects. Potential natural and cryptic mRNA splicing variants are identified, and null mutations are distinguished from leaky mutations. Mutations and rare SNPs were predicted in genomes of three cancer cell lines (U2OS, U251 and A431), which were supported by expression analyses. After filtering, tractable numbers of potentially deleterious variants are predicted by the software, suitable for further laboratory investigation. In these cell lines, novel functional variants comprised 6–17 inactivating mutations, 1–5 leaky mutations and 6–13 cryptic splicing mutations. Predicted effects were validated by RNA-seq analysis of the three aforementioned cancer cell lines, and expression microarray analysis of SNPs in HapMap cell lines.

Prediction of mutant mRNA splice isoforms by information theory-based exon definition

Mutations that affect mRNA splicing often produce multiple mRNA isoforms, resulting in complex molecular phenotypes. Definition of an exon and its inclusion in mature mRNA relies on joint recognition of both acceptor and donor splice sites. This study predicts cryptic and exon-skipping isoforms in mRNA produced by splicing mutations from the combined information contents (R(i), which measures binding-site strength, in bits) and distribution of the splice sites defining these exons. The total information content of an exon (R(i),total) is the sum of the R(i) values of its acceptor and donor splice sites, adjusted for the self-information of the distance separating these sites, that is, the gap surprisal. Differences between total information contents of an exon (ΔR(i,total)) are predictive of the relative abundance of these exons in distinct processed mRNAs. Constraints on splice site and exon selection are used to eliminate nonconforming and poorly expressed isoforms. Molecular phenotypes are computed by the Automated Splice Site and Exon Definition Analysis (http://splice.uwo.ca) server. Predictions of splicing mutations were highly concordant (85.2%; n = 61) with published expression data. In silico exon definition analysis will contribute to streamlining assessment of abnormal and normal splice isoforms resulting from mutations.

A cryptic BAP1 splice mutation in a family with uveal and cutaneous melanoma, and paraganglioma

Inactivating germ line BRCA1-associated protein-1 (BAP1) mutations have recently been reported in families with uveal or cutaneous malignant melanoma (UMM, CMM), mesothelioma, and meningioma. Although apparently predisposing to a wide range of tumors, the exact tumor spectrum associated with germ line BAP1 mutations has yet to be established. Here, we report a novel germ line BAP1 splice mutation, c.1708C>G (p.Leu570fs*40), in a multiple-case Danish UMM family with a spectrum of other tumors. Whole-exome sequencing identified an apparent missense mutation of BAP1 in UMM, CMM, as well as paraganglioma, breast cancer, and suspected mesothelioma cases in the family. Bioinformatic analysis and splicing assays demonstrated that this mutation creates a strong cryptic splice donor, resulting in aberrant splicing and a truncating frameshift of the BAP1 transcript. Somatic loss of the wild-type allele was also confirmed in the UMM and paraganglioma tumors. Our findings further support BAP1 as a melanoma susceptibility gene and extend the potential predisposition spectrum to paraganglioma.

Promoter variants in the MSMB gene associated with prostate cancer regulate MSMB/NCOA4 fusion transcripts

Beta-microseminoprotein (MSP)/MSMB is an immunoglobulin superfamily protein synthesized by prostate epithelial cells and secreted into seminal plasma. Variants in the promoter of the MSMB gene have been associated with the risk of prostate cancer (PCa) in several independent genome-wide association studies. Both MSMB and an adjacent gene, NCOA4, are subjected to transcriptional control via androgen response elements. The gene product of NCOA4 interacts directly with the androgen receptor as a co-activator to enhance AR transcriptional activity. Here, we provide evidence for the expression of full-length MSMB-NCOA4 fusion transcripts regulated by the MSMB promoter. The predominant MSMB-NCOA4 transcript arises by fusion of the 5'UTR and exons 1-2 of the MSMB pre-mRNA, with exons 2-10 of the NCOA4 pre-mRNA, producing a stable fusion protein, comprising the essential domains of NCOA4. Analysis of the splice sites of this transcript shows an unusually strong splice acceptor at NCOA4 exon 2 and the presence of Alu repeats flanking the exons potentially involved in the splicing event. Transfection experiments using deletion clones of the promoter coupled with luciferase reporter assays define a core MSMB promoter element located between -27 and -236 of the gene, and a negative regulatory element immediately upstream of the start codon. Computational network analysis reveals that the MSMB gene is functionally connected to NCOA4 and the androgen receptor signaling pathway. The data provide an example of how GWAS-associated variants may have multiple genetic and epigenetic effects.

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.

Identification of activated cryptic 5' splice sites using structure profiles and odds measure

The activation of cryptic 5′ splice sites (5′ SSs) is often related to human hereditary diseases. The DNA-based mutation screening strategies are commonly used to recognize the cryptic 5′ SSs, because features of the local DNA sequence can influence the choice of cryptic 5′ SSs. To improve the identification of the cryptic 5′ SSs, we developed a structure-based method, named SPO (structure profiles and odds measure), which combines two parameters, the structural feature derived from hydroxyl radical cleavage pattern and odds measure, to assess the likelihood of a cryptic 5′ SS activation in competing with its paired authentic 5′ SS. Compared to the current tools for identifying activated cryptic 5′ SSs, the SPO algorithm achieves higher prediction accuracy than the other methods, including MaxEnt, MDD, Markov model, weight matrix model, Shapiro and Senapathy matrix, R_i and ΔG. In addition, the predicted ΔSPO scores from the SPO algorithm exhibited a greater degree of correlation with the strength of cryptic 5′ SS activation than that measured from the other seven methods. In conclusion, the SPO algorithm provides an optimal identification of cryptic 5′ SSs, can be applied in designing mutagenesis experiments for various splicing events and may be helpful to investigate the relationship between structural variants and human hereditary diseases.

VAX1 mutation associated with microphthalmia, corpus callosum agenesis, and orofacial clefting: the first description of a VAX1 phenotype in humans

Vax1 and Vax2 have been implicated in eye development and the closure of the choroid fissure in mice and zebrafish. We sequenced the coding exons of VAX1 and VAX2 in 70 patients with anophthalmia/microphthalmia (A/M). In VAX1, we observed homozygosity for two successive nucleotide substitutions c.453G>A and c.454C>A, predicting p.Arg152Ser, in a proband of Egyptian origin with microphthalmia, small optic nerves, cleft lip/palate, and corpus callosum agenesis. This mutation affects an invariant residue in the homeodomain of VAX1 and was absent from 96 Egyptian controls. It is likely that the mutation results in a loss of function, as the mutation results in a phenotype similar to the Vax1 homozygous null mouse. We did not identify any mutations in VAX2. This is the first description of a phenotype associated with a VAX1 mutation in humans and establishes VAX1 as a new causative gene for A/M.

Disruption of PTPRO causes childhood-onset nephrotic syndrome

Idiopathic nephrotic syndrome (INS) is a genetically heterogeneous group of disorders characterized by proteinuria, hypoalbuminemia, and edema. Because it typically results in end-stage kidney disease, the steroid-resistant subtype (SRNS) of INS is especially important when it occurs in children. The present study included 29 affected and 22 normal individuals from 17 SRNS families; genome-wide analysis was performed with Affymetrix 250K SNP arrays followed by homozygosity mapping. A large homozygous stretch on chromosomal region 12p12 was identified in one consanguineous family with two affected siblings. Direct sequencing of protein tyrosine phosphatase receptor type O (PTPRO; also known as glomerular epithelial protein-1 [GLEPP1]) showed homozygous c.2627+1G>T donor splice-site mutation. This mutation causes skipping of the evolutionarily conserved exon 16 (p.Glu854_Trp876del) at the RNA level. Immunohistochemistry with GLEPP1 antibody showed a similar staining pattern in the podocytes of the diseased and control kidney tissues. We used a highly polymorphic intragenic DNA marker-D12S1303-to search for homozygosity in 120 Turkish and 13 non-Turkish individuals in the PodoNet registry. This analysis yielded 17 candidate families, and a distinct homozygous c.2745+1G>A donor splice-site mutation in PTPRO was further identified via DNA sequencing in a second Turkish family. This mutation causes skipping of exon 19, and this introduces a premature stop codon at the very beginning of exon 20 (p.Asn888Lysfs*3) and causes degradation of mRNA via nonsense-mediated decay. Immunohistochemical analysis showed complete absence of immunoreactive PTPRO. Ultrastructural alterations, such as diffuse foot process fusion and extensive microvillus transformation of podocytes, were observed via electron microscopy in both families. The present study introduces mutations in PTPRO as another cause of autosomal-recessive nephrotic syndrome.

Effects of splicing mutations on NF2-transcripts: transcript analysis and information theoretic predictions

This study examined the effects of 22 putative splicing mutations in the NF2 gene by means of transcript analysis and information theory based prediction. Fourteen mutations were within the dinucleotide acceptor and donor regions, often referred to as (AG/GT) sequences. Six were outside these dinucleotide regions but within the more broadly defined splicing regions used in the information theory based model. Two others were in introns and outside the broadly defined regions. Transcript analysis revealed exon skipping or activation of one or more cryptic splicing sites for 17 mutations. No alterations were found for the two intronic mutations and for three mutations in the broadly defined splicing regions. Concordance and partial concordance between the calculated predictions and the results of transcript analysis were found for 14 and 6 mutations, respectively. For two mutations, the predicted alteration was not found in the transcripts. Our results demonstrate that the effects of splicing mutations in NF2 are often complex and that information theory based analysis is helpful in elucidating the consequences of these mutations.

Manitoba-oculo-tricho-anal (MOTA) syndrome is caused by mutations in FREM1

BACKGROUND

Manitoba-oculo-tricho-anal (MOTA) syndrome is a rare condition defined by eyelid colobomas, cryptophthalmos and anophthalmia/microphthalmia, an aberrant hairline, a bifid or broad nasal tip, and gastrointestinal anomalies such as omphalocele and anal stenosis. Autosomal recessive inheritance had been assumed because of consanguinity in the Oji-Cre population of Manitoba and reports of affected siblings, but no locus or cytogenetic aberration had previously been described.

METHODS AND RESULTS

This study shows that MOTA syndrome is caused by mutations in FREM1, a gene previously mutated in bifid nose, renal agenesis, and anorectal malformations (BNAR) syndrome. MOTA syndrome and BNAR syndrome can therefore be considered as part of a phenotypic spectrum that is similar to, but distinct from and less severe than, Fraser syndrome. Re-examination of Frem1(bat/bat) mutant mice found new evidence that Frem1 is involved in anal and craniofacial development, with anal prolapse, eyelid colobomas, telecanthus, a shortened snout and reduced philtral height present in the mutant mice, similar to the human phenotype in MOTA syndrome.

CONCLUSIONS

The milder phenotypes associated with FREM1 deficiency in humans (MOTA syndrome and BNAR syndrome) compared to that resulting from FRAS1 and FREM2 loss of function (Fraser syndrome) are also consistent with the less severe phenotypes resulting from Frem1 loss of function in mice. Together, Fraser, BNAR and MOTA syndromes constitute a clinically overlapping group of FRAS-FREM complex diseases.

Identification of a novel splicing mutation in the ARSA gene in a patient with late-infantile form of metachromatic leukodystrophy

Metachromatic leukodystrophy (MLD; MIM 250100), a severe neurodegenerative disorder inherited as an autosomal recessive trait, is caused by mutations in the arylsulfatase A (ARSA) gene. Although several germ line ARSA mutations have been identified in patients with MLD of various ethnic backgrounds elsewhere in the world, no genetically confirmed cases of MLD have been reported in Korea. Recently, we identified a mutation in the ARSA gene of a Korean male with MLD. A male infant with late-infantile form of MLD had been admitted to our hospital for further examination. His neuromuscular symptoms, which included inability to walk at the age of 12 months, gradually worsened, even after allograft bone marrow transplantation; he died at the age of 9 yr. His elder brother had also been diagnosed with MLD. To confirm the presence of a genetic abnormality, all the coding exons of the ARSA gene and the flanking introns were amplified by PCR. A molecular analysis of the ARSA gene revealed both a novel heterozygous splicing mutation (c.1101+1G>T) in intron 6 and a heterozygous missense mutation in exon 2 (c.296G>A; Gly99Asp). The patient's elder brother who had MLD is believed to have had the same mutation, which may be correlated with a rapidly deteriorating clinical course. This study identified a novel mutation in the ARSA gene, related to a late-infantile form of MLD with a lethal clinical course and suggested that molecular diagnosis of patients may be useful in early diagnosis and for deciding intervention measures for their family members.

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.

Ontogenomic study of the relationship between number of gene splice variants and GO categorization

MOTIVATION

Splice variation plays important roles in evolution and cancer. Different splice variants of a gene may be characteristic of particular cellular processes, subcellular locations or organs. Although several genomic projects have identified splice variants, there have been no large-scale computational studies of the relationship between number of splice variants and biological function. The Gene Ontology (GO) and tools for leveraging GO, such as GoMiner, now make such a study feasible.

RESULTS

We partitioned genes into two groups: those with numbers of splice variants <or=b and >b (b=1,..., 10). Then we used GoMiner to determine whether any GO categories are enriched in genes with particular numbers of splice variants. Since there was no a priori 'appropriate' partition boundary, we studied those 'robust' categories whose enrichment did not depend on the selection of a particular partition boundary. Furthermore, because the distribution of splice variant number was a snapshot taken at a particular point in time, we confirmed that those observations were stable across successive builds of GenBank. A small number of categories were found for genes in the lower partitions. A larger number of categories were found for genes in the higher partitions. Those categories were largely associated with cell death and signal transduction. Apoptotic genes tended to have a large repertoire of splice variants, and genes with splice variants exhibited a distinctive 'apoptotic island' in clustered image maps (CIMs).

AVAILABILITY

Supplementary tables and figures are available at URL http://discover.nci.nih.gov/OG/supplementaryMaterials.html. The Safari browser appears to perform better than Firefox for these particular items.