Recurrent CNVs and SNVs at the NPHP1 locus contribute pathogenic alleles to Bardet-Biedl syndrome

Homozygosity for a recurrent 290 kb deletion of NPHP1 is the most frequent cause of isolated nephronophthisis (NPHP) in humans. A deletion of the same genomic interval has also been detected in individuals with Joubert syndrome (JBTS), and in the mouse, Nphp1 interacts genetically with Ahi1, a known JBTS locus. Given these observations, we investigated the contribution of NPHP1 in Bardet-Biedl syndrome (BBS), a ciliopathy of intermediate severity. By using a combination of array-comparative genomic hybridization, TaqMan copy number assays, and sequencing, we studied 200 families affected by BBS. We report a homozygous NPHP1 deletion CNV in a family with classical BBS that is transmitted with autosomal-recessive inheritance. Further, we identified heterozygous NPHP1 deletions in two more unrelated persons with BBS who bear primary mutations at another BBS locus. In parallel, we identified five families harboring an SNV in NPHP1 resulting in a conserved missense change, c.14G>T (p.Arg5Leu), that is enriched in our Hispanic pedigrees; in each case, affected individuals carried additional bona fide pathogenic alleles in another BBS gene. In vivo functional modeling in zebrafish embryos demonstrated that c.14G>T is a loss-of-function variant, and suppression of nphp1 in concert with each of the primary BBS loci found in our NPHP1-positive pedigrees exacerbated the severity of the phenotype. These results suggest that NPHP1 mutations are probably rare primary causes of BBS that contribute to the mutational burden of the disorder.

Draft genome sequences and description of Lactobacillus rhamnosus strains L31, L34, and L35

Lactobacillus rhamnosus is a facultative, lactic acid bacterium in the phylum Firmicutes. Lactobacillus spp. are generally considered beneficial, and specific strains of L. rhamnosus are validated probiotics. We describe the draft genomes of three L. rhamnosus strains (L31, L34, and L35) isolated from the feces of Thai breastfed infants, which exhibit anti-inflammatory properties in vitro. The three genomes range between 2.8 - 2.9 Mb, and contain approximately 2,700 protein coding genes.

Abstract IA21: Developing a new functional classification of lung cancer based on tumor acquired vulnerabilities

We have been comprehensively screening for “vulnerabilities” that have been acquired during the multi-step pathogenesis of lung cancer cells but are not present in normal lung epithelial cells to identify genetic and chemical perturbations that will selectively kill lung cancer. We think many of these have occurred to allow the lung cancers to undergo/tolerate “oncogene addiction.” We tested a sub-panel of 12-15 non-small cell lung cancer (NSCLC) lines that covers the known molecular spectra of lung cancer with genome wide siRNA and large scale chemical library (~250,000 compounds) and natural products in vitro screens to identify “hits” that will kill (suppress the growth of) lung cancer cells but not normal human bronchial epithelial cells and that also only kill a subset of lung cancer cells providing two types of specificity. “Hits” from these broad screens are then tested (including detailed drug concentration curves) across a large panel of lung cancer lines (~100) representing a variety of lung cancer histologic and molecular oncogenotypes. Other versions of these screens include the intensive use of “mini-libraries” each containing 50 – 150 gene targets by siRNAs or shRNAs, or ~200 defined drugs to explore pathways in detail in tests of over 70 NSCLCs. Examples include: nuclear receptors and their co-regulators (120 genes); cancer stem cell pathways (50 genes); chromatin remodelers (75 genes) and identified lung cancer mutated driver oncogenes (175 genes). In addition to the in vitro tests, we have developed in vivo (xenograft) tests where shRNA mini-libraries are introduced into tumor cells at high representation which are grown as xenografts, analyzed by NexGen sequencing and shRNAs identified that drop out or are retained in xenografts compared to in vitro grown cells to identify vulnerabilities that are only detected in the in vivo situation. All of the data are then related to the large legacy molecular datasets associated with the lung cancer lines (including whole exome sequence analyses and genome wide mRNA, copy number variation, methylation, miR expression data and proteomics data). In addition, detailed chemical and pharmacokinetic analyses for favorable drug properties and subsequent chemical modifications also occur for the chemical compounds to progress those towards potential clinical studies. The results of all of these analyses have identified ~300 new chemical compounds and ~300 genetic hits all of which show selectivity for lung cancer over normal lung cells and selectivity for subtypes of lung cancer. The chemical and genes hits are being compared to the tumor molecular information and integrated in turn through a “connectivity map” type of approach – to identify drugs and gene hits involving the same pathways. The molecular correlates of the tumor lines are related to similar molecular changes in patient derived xenografts and patient tumor specimens to provide a connection of the molecular subtype-selective vulnerabilities (“enrollment biomarkers”) between the preclinical response phenotypes and patient tumor specimens. From these data we find lung cancers can be classified into groups (“clades”) that represent functional vulnerabilities to the gene and chemical compound hits and these in turn can be related to molecular abnormalities in tumors. One example of this is our detailed analyses a matched lung adenocarcinoma/normal lung epithelial cell model derived from the same patient which identified three distinct target/response-indicator pairings that are represented a significant frequencies (6-16%) in the lung adenocarcinoma population (Kim et al. Cell 155:552, 2013). These include three totally novel lung cancer selective targeted therapies: NLRP3 mutation/inflammasome activation-dependent FLIP addiction; co-occurring KRAS and LKB1 mutation-driven COPI addiction; and selective sensitivity to a synthetic indolotriazine that is specified by a seven-gene expression signature. Our panel of “hits” provide the opportunity to identify all potential therapeutic targets for lung cancer, while the molecular correlates will allow “personalization” of these new therapies going forward in preclinical and clinical translation. (Supported by NCI SPORE P50CA70907, NCI CTD2N, CPRIT, UTSW CCSG P30CA142543)

Citation Format: John D. Minna, Adi Gazdar, Alexander Augustyn, Rebecca Britt, Ryan Carstens, Patrick Dospoy, Boning Gao, Luc Girard, Suzie Hight, Kenneth Huffman, Jill Larsen, Michael Peyton, Chunli Shao, David Mangelsdorf, Rolf Brekken, Ralph Deberardinis, Pei-Hsuan Chen, Carmen Behrens, Lauren Byers, John Heymach, Jack Roth, Ignacio Wistuba, Yang Xie, Caleb Davis, David Wheeler, Richard Gibbs, Edward Marcotte, Joseph Ready, Deepak Nijhawan, Noelle Williams, Steven McKnight, Bruce Posner, John MacMillan, Michael Roth, Michael White. Developing a new functional classification of lung cancer based on tumor acquired vulnerabilities. [abstract]. In: Proceedings of the AACR-IASLC Joint Conference on Molecular Origins of Lung Cancer; 2014 Jan 6-9; San Diego, CA. Philadelphia (PA): AACR; Clin Cancer Res 2014;20(2Suppl):Abstract nr IA21.

The somatic genomic landscape of glioblastoma

We describe the landscape of somatic genomic alterations based on multidimensional and comprehensive characterization of more than 500 glioblastoma tumors (GBMs). We identify several novel mutated genes as well as complex rearrangements of signature receptors, including EGFR and PDGFRA. TERT promoter mutations are shown to correlate with elevated mRNA expression, supporting a role in telomerase reactivation. Correlative analyses confirm that the survival advantage of the proneural subtype is conferred by the G-CIMP phenotype, and MGMT DNA methylation may be a predictive biomarker for treatment response only in classical subtype GBM. Integrative analysis of genomic and proteomic profiles challenges the notion of therapeutic inhibition of a pathway as an alternative to inhibition of the target itself. These data will facilitate the discovery of therapeutic and diagnostic target candidates, the validation of research and clinical observations and the generation of unanticipated hypotheses that can advance our molecular understanding of this lethal cancer.

Integrative annotation of variants from 1092 humans: application to cancer genomics

Interpreting variants, especially noncoding ones, in the increasing number of personal genomes is challenging. We used patterns of polymorphisms in functionally annotated regions in 1092 humans to identify deleterious variants; then we experimentally validated candidates. We analyzed both coding and noncoding regions, with the former corroborating the latter. We found regions particularly sensitive to mutations ("ultrasensitive") and variants that are disruptive because of mechanistic effects on transcription-factor binding (that is, "motif-breakers"). We also found variants in regions with higher network centrality tend to be deleterious. Insertions and deletions followed a similar pattern to single-nucleotide variants, with some notable exceptions (e.g., certain deletions and enhancers). On the basis of these patterns, we developed a computational tool (FunSeq), whose application to ~90 cancer genomes reveals nearly a hundred candidate noncoding drivers.

Whole-genome sequence-based analysis of high-density lipoprotein cholesterol

We describe initial steps for interrogating whole genome sequence (WGS) data to characterize the genetic architecture of a complex trait, such as high density lipoprotein cholesterol (HDL-C). We estimate that common variation contributes more to HDL-C heritability than rare variation, and screening for Mendelian dyslipidemia variants identified individuals with extreme HDL-C. WGS analyses highlight the value of regulatory and non-protein coding regions of the genome in addition to protein coding regions.

Abstract 1114: Polymerase epsilon (POLE) mutations and mutator phenotypes in colorectal and endometrial tumors

Colorectal rectal and endometrial cancers are divided into microsatellite instable (MSI) and microsatellite stable (MSS) types. MSS tumors have chromosome instability with mutation rates of 1-10/Mb. MSI patients have a better prognosis than MSS patients and are treated less aggressively. These tumors exhibit microsatellite instability (MSI), and the CpG island methylator phenotype (CIMP), an increased mutation rate (10-100/Mb, hypermutated). Here we demonstrate a novel class of tumor in these two diseases exhibiting an ultramutator phenotype with mutation rates exceeding 100/ Mb, MSS and chromosome stable. The ultra-high mutation rates appear to be caused by recurrent mutations in the exonuclease domain of DNA polymerase epsilon (POLE).

Matched tumor and normal whole exomes for 512 colorectal and 248 endometrial cancers (TCGA data control center) were evaluated for mutation frequency and MSI status. Among the colorectal tumors, 70/512(14%) were hyper mutated, endometroid tumors had 65/248(26%) hypermutated tumors. Within the hypermutated tumors, we identified tumors with extremely high mutation rates (100 mutations/Mb) termed “ultramutated”. 14 (3%) colorectal and 17 (7%) endometrial were ultramutated. 100% of the colorectal and endometrial ultramutated contain an exonuclease domain mutation in POLE with three recurrently mutated positions: P286R/H/S and V411L and S459F, these mutations are found exclusively in ultramutated samples. The functional relevance of these mutations has been demonstrated by previous mutational analysis of the exonuclease domain in bacteria phage (T4) yeast,and mice. Disruption of residues in the active site leads to high error rates, a mutator phenotype and tumor formation. The defined exonuclease active site aligns closely with the recurrent mutation found in the ultramutated tumors.

The ultramuted tumors have a distinct mutation spectrum, occur in context and exhibit strand biases. POLE exonuclease domain mutation tumors have high rates of C to A and T to G mutations with low rates of C to G and T to A. Changes of C to A occur predominately in context of 3’ and 5’ flanking T bases. This mutation spectrum and context of flanking bases suggests replicative strand bias. We found that mutations on the leading strand, exhibited a preference for C to A over G to T (60:40). Suggesting transcription-coupled repair (TCR) plays a role in ultramutation.

Our results demonstrate that ultra mutated tumors are driven by recurrent mutations in the exonuclease domain of DNA polymerase epsilon. The distinct mutation spectrum, strand bias and mutation context all suggest TCR is involved in ultramutated tumors. We describe a previously unknown role of POLE in colorectal tumors as well as identification of a previously undescribed hot spot for cancer mutations.

Citation Format: Eve Shinbrot, Nils Weinhold, Nikolaus Schultz, Lawrence A. Donehower, Jennifer Drummond, Kyle Chang, Richard Gibbs, Chris Sander, David A. Wheeler. Polymerase epsilon (POLE) mutations and mutator phenotypes in colorectal and endometrial tumors. [abstract]. In: Proceedings of the 104th Annual Meeting of the American Association for Cancer Research; 2013 Apr 6-10; Washington, DC. Philadelphia (PA): AACR; Cancer Res 2013;73(8 Suppl):Abstract nr 1114. doi:10.1158/1538-7445.AM2013-1114

Abstract 5066: Genome sequence and epigenome map of bladder cancer development from occult field effects

We are pursuing a unique strategy to identify early events in bladder cancer development by the genome-wide analysis of somatic changes in tissues from resected bladders and characterized pathologically as urothelial carcinoma. We then extend this analysis into precursor intraurothelial lesions and normal tissue. Here we report on the construction of a high resolution whole-organ histologic and genetic map of bladder cancer development that tracks its progression from in situ precursor conditions to invasive clinically aggressive cancer using total genome and exome sequencing and complemented with HumanOmni 2.5-8 v1 SNP-based Illumina genotyping, combined with CpG island methylation and cDNA microarray platforms derived from 12 whole-organ cystectomy specimens with invasive urothelial carcinoma. This approach provided a high resolution integrated map of structural genome, epigenome, and transcriptome alterations that parallel the development of bladder cancer from early field effects to invasive disease. By analyzing genome-wide patterns of structural genomic, epigenomic, and transcriptome changes, we found three major waves of alterations associated with growth advantage of successive clones reflecting a stepwise conversion of normal urothelial cells into cancer cells. This approach is significant for the identification of specific chromosomal regions as well as epigenetically and transcriptionally modified gene ontologies involved in early clonal expansion of intraurothelial neoplasia and their clonal evolution to carcinoma in situ progressing to invasive cancer. These studies suggest involvement of specific gene networks in early clonal expansion of intraurothelial neoplasia antecedent to microscopically recognizable precursor lesions such as dysplasia and carcinoma in situ. This map represents a prerequisite for specific studies on the mechanisms of incipient phases of bladder cancer development and should have a significant impact on the understanding of human carcinogenesis.

Citation Format: {Authors}. {Abstract title} [abstract]. In: Proceedings of the 103rd Annual Meeting of the American Association for Cancer Research; 2012 Mar 31-Apr 4; Chicago, IL. Philadelphia (PA): AACR; Cancer Res 2012;72(8 Suppl):Abstract nr 5066. doi:1538-7445.AM2012-5066

The midterm experience of tapered stent grafts in the endovascular management of iliac artery aneurysms with unfavorable anatomy

We report our experience and the midterm results of a modern technique for endovascular management of isolated iliac artery aneurysms (IAAs) with unfavorable neck anatomy, which involves the inversion of an iliac leg of a Zenith stent graft. Patients who underwent endovascular IAA repair from 2002 to 2010 were reviewed. A total of 12 patients, with a mean age of 77.6 years, underwent endovascular repair of 13 IAAs. Mean size of the aneurysms was 54.6 mm (range 34-133 mm). Mean proximal neck diameter was 18 mm (range 15-22 mm). In 7 patients, the length of the proximal neck was <15 mm (10-14 mm). Only 1 patient developed thrombosis of the stent graft immediately after the operation. Patients were followed up for a mean of 31.5 months (range 18-72 months). Our midterm results demonstrate the durability of this technique in the management of iliac aneurysms with unfavorable anatomy.

SSTR5 P335L monoclonal antibody differentiates pancreatic neuroendocrine neuroplasms with different SSTR5 genotypes

BACKGROUND

Somatostatin receptor type 5 (SSTR5) P335L is a hypofunctional, single nucleotide polymorphism of SSTR5 with implications in the diagnostics and therapy of pancreatic neuroendocrine neoplasms. The purpose of this study is to determine whether a SSTR5 P335L-specific monoclonal antibody could sufficiently differentiate pancreatic neuroendocrine neoplasms (PNENs) with different SSTR5 genotypes.

METHODS

Cellular proliferation rate, SSTR5 mRNA level, and SSTR5 protein level were measured by performing MTS assay, a quantitative reverse transcription polymerase chain reaction study, Western blot analysis, and immunohistochemistry, respectively. SSTR5 genotype was determined with the TaqMan SNP Genotyping assay (Applied Biosystems, Foster City, CA).

RESULTS

We found that the SSTR5 analogue RPL-1980 inhibited cellular proliferation of CAPAN-1 cells more than that of PANC-1 cells. Only PANC-1 (TT) cells, but not CAPAN-1 (CC) cells expressed SSTR5 P335L. In 29 white patients with PNENs, 38% had a TT genotype for SSTR5 P335L, 24% had a CC genotype for WT SSTR5, and 38% hada CT genotype for both SSTR5 P335L and WT SSTR5. Immunohistochemistry using SSTR5 P335L monoclonal antibody detected immunostaining signals only from the neuroendocrine specimens with TT and CT genotypes, but not those with CC genotypes.

CONCLUSION

A SSTR5 P335L monoclonal antibody that specifically recognizes SSTR5 P335L but not WT SSTR5 could differentiate PNENs with different SSTR5 genotypes, thereby providing a potential tool for the clinical diagnosis of PNEN.

Mutations in SERPINF1 cause osteogenesis imperfecta type VI

Osteogenesis imperfecta (OI) is a spectrum of genetic disorders characterized by bone fragility. It is caused by dominant mutations affecting the synthesis and/or structure of type I procollagen or by recessively inherited mutations in genes responsible for the posttranslational processing/trafficking of type I procollagen. Recessive OI type VI is unique among OI types in that it is characterized by an increased amount of unmineralized osteoid, thereby suggesting a distinct disease mechanism. In a large consanguineous family with OI type VI, we performed homozygosity mapping and next-generation sequencing of the candidate gene region to isolate and identify the causative gene. We describe loss of function mutations in serpin peptidase inhibitor, clade F, member 1 (SERPINF1) in two affected members of this family and in an additional unrelated patient with OI type VI. SERPINF1 encodes pigment epithelium-derived factor. Hence, loss of pigment epithelium-derived factor function constitutes a novel mechanism for OI and shows its involvement in bone mineralization.

Gastrointestinal microbiome signatures of pediatric patients with irritable bowel syndrome

BACKGROUND & AIMS

The intestinal microbiomes of healthy children and pediatric patients with irritable bowel syndrome (IBS) are not well defined. Studies in adults have indicated that the gastrointestinal microbiota could be involved in IBS.

METHODS

We analyzed 71 samples from 22 children with IBS (pediatric Rome III criteria) and 22 healthy children, ages 7-12 years, by 16S ribosomal RNA gene sequencing, with an average of 54,287 reads/stool sample (average 454 read length = 503 bases). Data were analyzed using phylogenetic-based clustering (Unifrac), or an operational taxonomic unit (OTU) approach using a supervised machine learning tool (randomForest). Most samples were also hybridized to a microarray that can detect 8741 bacterial taxa (16S rRNA PhyloChip).

RESULTS

Microbiomes associated with pediatric IBS were characterized by a significantly greater percentage of the class γ-proteobacteria (0.07% vs 0.89% of total bacteria, respectively; P < .05); 1 prominent component of this group was Haemophilus parainfluenzae. Differences highlighted by 454 sequencing were confirmed by high-resolution PhyloChip analysis. Using supervised learning techniques, we were able to classify different subtypes of IBS with a success rate of 98.5%, using limited sets of discriminant bacterial species. A novel Ruminococcus-like microbe was associated with IBS, indicating the potential utility of microbe discovery for gastrointestinal disorders. A greater frequency of pain correlated with an increased abundance of several bacterial taxa from the genus Alistipes.

CONCLUSIONS

Using 16S metagenomics by PhyloChip DNA hybridization and deep 454 pyrosequencing, we associated specific microbiome signatures with pediatric IBS. These findings indicate the important association between gastrointestinal microbes and IBS in children; these approaches might be used in diagnosis of functional bowel disorders in pediatric patients.

The functional spectrum of low-frequency coding variation

Background

Rare coding variants constitute an important class of human genetic variation, but are underrepresented in current databases that are based on small population samples. Recent studies show that variants altering amino acid sequence and protein function are enriched at low variant allele frequency, 2 to 5%, but because of insufficient sample size it is not clear if the same trend holds for rare variants below 1% allele frequency.

Results

The 1000 Genomes Exon Pilot Project has collected deep-coverage exon-capture data in roughly 1,000 human genes, for nearly 700 samples. Although medical whole-exome projects are currently afoot, this is still the deepest reported sampling of a large number of human genes with next-generation technologies. According to the goals of the 1000 Genomes Project, we created effective informatics pipelines to process and analyze the data, and discovered 12,758 exonic SNPs, 70% of them novel, and 74% below 1% allele frequency in the seven population samples we examined. Our analysis confirms that coding variants below 1% allele frequency show increased population-specificity and are enriched for functional variants.

Conclusions

This study represents a large step toward detecting and interpreting low frequency coding variation, clearly lays out technical steps for effective analysis of DNA capture data, and articulates functional and population properties of this important class of genetic variation.

Exome sequencing of ion channel genes reveals complex profiles confounding personal risk assessment in epilepsy

Ion channel mutations are an important cause of rare Mendelian disorders affecting brain, heart, and other tissues. We performed parallel exome sequencing of 237 channel genes in a well-characterized human sample, comparing variant profiles of unaffected individuals to those with the most common neuronal excitability disorder, sporadic idiopathic epilepsy. Rare missense variation in known Mendelian disease genes is prevalent in both groups at similar complexity, revealing that even deleterious ion channel mutations confer uncertain risk to an individual depending on the other variants with which they are combined. Our findings indicate that variant discovery via large scale sequencing efforts is only a first step in illuminating the complex allelic architecture underlying personal disease risk. We propose that in silico modeling of channel variation in realistic cell and network models will be crucial to future strategies assessing mutation profile pathogenicity and drug response in individuals with a broad spectrum of excitability disorders.

Characterization of single-nucleotide variation in Indian-origin rhesus macaques (Macaca mulatta)

BACKGROUND

Rhesus macaques are the most widely utilized nonhuman primate model in biomedical research. Previous efforts have validated fewer than 900 single nucleotide polymorphisms (SNPs) in this species, which limits opportunities for genetic studies related to health and disease. Extensive information about SNPs and other genetic variation in rhesus macaques would facilitate valuable genetic analyses, as well as provide markers for genome-wide linkage analysis and the genetic management of captive breeding colonies.

RESULTS

We used the available rhesus macaque draft genome sequence, new sequence data from unrelated individuals and existing published sequence data to create a genome-wide SNP resource for Indian-origin rhesus monkeys. The original reference animal and two additional Indian-origin individuals were resequenced to low coverage using SOLiD™ sequencing. We then used three strategies to validate SNPs: comparison of potential SNPs found in the same individual using two different sequencing chemistries, and comparison of potential SNPs in different individuals identified with either the same or different sequencing chemistries. Our approach validated approximately 3 million SNPs distributed across the genome. Preliminary analysis of SNP annotations suggests that a substantial number of these macaque SNPs may have functional effects. More than 700 non-synonymous SNPs were scored by Polyphen-2 as either possibly or probably damaging to protein function and these variants now constitute potential models for studying functional genetic variation relevant to human physiology and disease.

CONCLUSIONS

Resequencing of a small number of animals identified greater than 3 million SNPs. This provides a significant new information resource for rhesus macaques, an important research animal. The data also suggests that overall genetic variation is high in this species. We identified many potentially damaging non-synonymous coding SNPs, providing new opportunities to identify rhesus models for human disease.

Scalloped thoracic stent-graft for treatment of aortic arch aneurysms with unfavourable landing zones

Endovascular treatments are limited in cases of thoracic aortic aneurysms extending up or proximal to the origin of the left subclavian artery (LSCA). In such cases, the LSCA is usually either occluded or revascularised. We report our first experience of four patients who underwent thoracic aneursym treatment with new custom-made grafts with a scallop in situ for the LSCA. The graft is tailor made per case, and a re-enforced scallop is positioned proximally allowing for the stent to be deployed beyond the origin of the LSCA; the origin of the LSCA remains patent, thus negating the need for revascularisation of the head and neck vessels on the left. The stent contains markers for identifying the scallop and are located along the midline to ensure correct alignment. All of the patients who underwent this procedure had technical success with flow through the LSCA both immediately after stent deployment and on follow-up imaging. This new stent has further expanded endovascular treatment options for patients with thoracic aneurysms extending up to and beyond the LSCA, which can play a part in improving outcome and decreasing mortality rates because surgery for revascularization will not be needed.

The Drosophila melanogaster transcriptome by paired-end RNA sequencing

RNA-seq was used to generate an extensive map of the Drosophila melanogaster transcriptome by broad sampling of 10 developmental stages. In total, 142.2 million uniquely mapped 64-100-bp paired-end reads were generated on the Illumina GA II yielding 356× sequencing coverage. More than 95% of FlyBase genes and 90% of splicing junctions were observed. Modifications to 30% of FlyBase gene models were made by extension of untranslated regions, inclusion of novel exons, and identification of novel splicing events. A total of 319 novel transcripts were identified, representing a 2% increase over the current annotation. Alternate splicing was observed in 31% of D. melanogaster genes, a 38% increase over previous estimations, but significantly less than that observed in higher organisms. Much of this splicing is subtle such as tandem alternate splice sites.