AT-AC pre-mRNA splicing mechanisms and conservation of minor introns in voltage-gated ion channel genes

Widespread cleavage of A-to-I hyperediting substrates

A-to-I RNA editing is the conversion of adenosine to inosine in double-stranded cellular and viral RNAs. Recently, abundant hyperediting of human transcripts, affecting thousands of genes, has been reported. Most of these editing sites are confined to intramolecular hairpin double-stranded RNA (dsRNA) structures formed by pairing of neighboring, reversely oriented, primate-specific Alu repeats. The biological implication of this extensive modification is still a mystery. A number of studies have shown that heavily edited transcripts are often retained in the nucleus. A recent study found that the edited region in transcripts of the mouse Slc7a2 gene is post-transcriptionally cleaved upon stress, enabling the release of the mRNA to the cytoplasm, followed by its translation. Here, we aim to test whether this scenario might be relevant for many other hyperedited Alu targets. Bioinformatics analysis of publicly available mRNA and expressed sequence tag data provides evidence showing that neighboring, reversely oriented, Alu elements are often cleaved at both ends of the region harboring the inverted repeats followed by rejoining of the two parts of the transcript on both sides of the inverted repeats, resulting in almost inosine-free mRNA products. Deleted segments vary among transcripts of the same gene and are not flanked by the canonical splicing signal sequences. The tissue distribution of these events seems to correlate with known A-to-I editing patterns, suggesting that it depends on the dsRNA structure being edited. Results are experimentally verified by polymerase chain reaction and cloning data. A database of 566 human and 107 mouse putative cleavage loci is supplied.

Oxalate decarboxylase of the white-rot fungus Dichomitus squalens demonstrates a novel enzyme primary structure and non-induced expression on wood and in liquid cultures

Oxalate decarboxylase (ODC) catalyses the conversion of oxalic acid to formic acid and CO2 in bacteria and fungi. In wood-decaying fungi the enzyme has been linked to the regulation of intra- and extracellular quantities of oxalic acid, which is one of the key components in biological decomposition of wood. ODC enzymes are biotechnologically interesting for their potential in diagnostics, agriculture and environmental applications, e.g. removal of oxalic acid from industrial wastewaters. We identified a novel ODC in mycelial extracts of two wild-type isolates of Dichomitus squalens, and cloned the corresponding Ds-odc gene. The primary structure of the Ds-ODC protein contains two conserved Mn-binding cupin motifs, but at the N-terminus, a unique, approximately 60 aa alanine-serine-rich region is found. Real-time quantitative RT-PCR analysis confirmed gene expression when the fungus was cultivated on wood and in liquid medium. However, addition of oxalic acid in liquid cultures caused no increase in transcript amounts, thereby indicating a constitutive rather than inducible expression of Ds-odc. The detected stimulation of ODC activity by oxalic acid is more likely due to enzyme activation than to transcriptional upregulation of the Ds-odc gene. Our results support involvement of ODC in primary rather than secondary metabolism in fungi.

Identification and expression analysis of two zebrafish E2F5 genes during oogenesis and development

During the process of screening genes that are preferentially expressed in zebrafish germ cells, we isolated an E2F family member zE2F5b, which is a truncated isoform of zebrafish E2F5a. The E2F transcription factors are thought to play an essential role in cell cycle progression. Reverse transcription polymerase chain reaction and whole mount in situ hybridization analyses have demonstrated that both zE2F5a and zE2F5b mRNAs are predominantly expressed in stage I and II oocytes. In early development, the two zE2F5 transcripts are expressed by many parts of the embryo and are then restricted to the brain area. Interestingly, zE2F5b mRNA is detected later in the lens and optic tectum while zE2F5a mRNA is absent. Our data suggests that E2F5 may play important roles in zebrafish brain, eye and ovarian development.

The conserved 3' end domain of U6atac snRNA can direct U6 snRNA to the minor spliceosome

U6 and U6atac snRNAs play analogous critical roles in the major U2-dependent and minor U12-dependent spliceosomes, respectively. Previous results have shown that most of the functional cores of these two snRNAs are either highly similar in sequence or functionally interchangeable. Thus, a mechanism must exist to restrict each snRNA to its own spliceosome. Here we show that a chimeric U6 snRNA containing the unique and highly conserved 3' end domain of U6atac snRNA is able to function in vivo in U12-dependent spliceosomal splicing. Function of this chimera required the coexpression of a modified U4atac snRNA; U4 snRNA could not substitute. Partial deletions of this element in vivo, as well as in vitro antisense experiments, showed that the 3' end domain of U6atac snRNA is necessary to direct the U4atac/U6atac.U5 tri-snRNP to the forming U12-dependent spliceosome. In vitro experiments also uncovered a role for U4atac snRNA in this targeting.

Solution structure of the U11-48K CHHC zinc-finger domain that specifically binds the 5' splice site of U12-type introns

The formation of stable 18S U11/U12 di-snRNPs before their association with the pre-mRNA is a characteristic feature of the minor spliceosome. During the spliceosomal assembly, the 18S snRNP binds cooperatively to the introns' 5' splice and branch point site. The molecular basis for this recognition is still unknown. Here, we report the solution structure of the U11-48K CHHC Zn finger, a domain unique to the minor spliceosome. The CHHC Zn-finger structure revealed an unexpected similarity to the TFIIIA domains, with distinct features originating from the type and separation of the zinc-coordinating residues. We show that this domain specifically binds the 5' splice site sequence of U12-type introns when base paired to U11 snRNA in vitro and hence may contribute to the U12 intron recognition. We propose a model in which the U11-48K Zn finger stabilizes U11-5' splice site base pairing and thus plays an important role during the minor spliceosome assembly.

Porcine Fc gammaRIII isoforms are generated by alternative splicing

The Fc gammaRIII plays an essential role in antibody-mediating effector functions of immune cells. Here, we report that transcripts encoding porcine Fc gammaRIII isoforms are generated by alternative splicing. Fc gammaRIII a.1 is expressed on the cell surface while Fc gammaRIII a.2 is secreted from the transfected cells due to a partial deletion of the transmembrane domain. Interestingly, a putative soluble Fc gammaRIII (sCD16) is detected in circulating plasma. Both Fc gammaRIII a.2 and sCD16 exhibit a similar molecular mass (approximately 35 kDa) and contain the extracellular D2 domains that are structurally intact. Despite the D2 domain deletion, Fc gammaRIII a.3 associates with Fc gammaRIII a.1. Hence, these results suggest one possibility that three Fc gammaRIII isoforms differentially modulate Fc gammaRIII-mediated immune responses in the porcine system. Furthermore, we demonstrate that a cytolytic triggering G7 monoclonal antibody recognizes the D2 domain that is responsible for interacting with the immune complexes and subsequent activations of porcine innate immune cells. This result suggests that the D2 domain is a target region to develop therapeutic antibodies that regulate the FcR-mediated immune responses.

Rearrangement of actin cytoskeleton mediates invasion of Lotus japonicus roots by Mesorhizobium loti

Infection thread-dependent invasion of legume roots by rhizobia leads to internalization of bacteria into the plant cells, which is one of the salient features of root nodule symbiosis. We found that two genes, Nap1 (for Nck-associated protein 1) and Pir1 (for 121F-specific p53 inducible RNA), involved in actin rearrangements were essential for infection thread formation and colonization of Lotus japonicus roots by its natural microsymbiont, Mesorhizobium loti. nap1 and pir1 mutants developed an excess of uncolonized nodule primordia, indicating that these two genes were not essential for the initiation of nodule organogenesis per se. However, both the formation and subsequent progression of infection threads into the root cortex were significantly impaired in these mutants. We demonstrate that these infection defects were due to disturbed actin cytoskeleton organization. Short root hairs of the mutants had mostly transverse or web-like actin filaments, while bundles of actin filaments in wild-type root hairs were predominantly longitudinal. Corroborating these observations, temporal and spatial differences in actin filament organization between wild-type and mutant root hairs were also observed after Nod factor treatment, while calcium influx and spiking appeared unperturbed. Together with various effects on plant growth and seed formation, the nap1 and pir1 alleles also conferred a characteristic distorted trichome phenotype, suggesting a more general role for Nap1 and Pir1 in processes establishing cell polarity or polar growth in L. japonicus.

A putative autonomous 20.5 kb-CACTA transposon insertion in an F3'H allele identifies a new CACTA transposon subfamily in Glycine max

BACKGROUND

The molecular organization of very few genetically defined CACTA transposon systems have been characterized thoroughly as those of Spm/En in maize, Tam1 of Antirrhinum majus Candystripe1 (Cs1) from Sorghum bicolor and CAC1 from Arabidopsis thaliana, for example. To date, only defective deletion derivatives of CACTA elements have been described for soybean, an economically important plant species whose genome sequence will be completed in 2008.

RESULTS

We identified a 20.5 kb insertion in a soybean flavonoid 3'-hydroxylase (F3'H) gene representing the t* allele (stable gray trichome color) whose origin traces to a single mutable chimeric plant displaying both tawny and gray trichomes. This 20.5 kb insertion has the molecular structure of a putative autonomous transposon of the CACTA family, designated Tgmt*. It encodes a large gene that was expressed in two sister isolines (T* and tm) of the stable gray line (t*) from which Tgmt* was isolated. RT-PCR derived cDNAs uncovered the structure of a large precursor mRNA as well as alternatively spliced transcripts reminiscent of the TNPA-mRNA generated by the En-1 element of maize but without sequence similarity to the maize TNPA. The larger mRNA encodes a transposase with a tnp2 and TNP1-transposase family domains. Because the two soybean lines expressing Tgmt* were derived from the same mutable chimeric plant that created the stable gray trichome t* allele line from which the element was isolated, Tgmt* has the potential to be an autonomous element that was rapidly inactivated in the stable gray trichome t* line. Comparison of Tgmt* to previously described Tgm elements demonstrated that two subtypes of CACTA transposon families exist in soybean based on divergence of their characteristic subterminal repeated motifs and their transposases. In addition, we report the sequence and annotation of a BAC clone containing the F3'H gene (T locus) which was interrupted by the novel Tgmt* element in the gray trichome allele t*.

CONCLUSION

The molecular characterization of a 20.5 kb insertion in the flavonoid 3'-hydroxylase (F3'H) gene of a soybean gray pubescence allele (t*) identified the structure of a CACTA transposon designated Tgmt*. Besides the terminal inverted repeats and subterminal repeated motifs,Tgmt* encoded a large gene with two putative functions that are required for excision and transposition of a CACTA element, a transposase and the DNA binding protein known to associate to the subterminal repeated motifs. The degree of dissimilarity between Tgmt* transposase and subterminal repeated motifs with those of previously characterized defective CACTA elements (Tgm1-7) were evidence of the existence of two subfamilies of CACTA transposons in soybean, an observation not previously reported in other plants. In addition, our analyses of a genetically active and potentially autonomous element sheds light on the complete structure of a soybean element that is useful for annotation of the repetitive fraction of the soybean genome sequence and may prove useful for transposon tagging or transposon display experiments in different genetic lines.

Noncanonical and canonical splice sites: a novel mutation at the rare noncanonical splice-donor cut site (IVS4+1A>G) of SEDL causes variable splicing isoforms in X-linked spondyloepiphyseal dysplasia tarda

X-linked spondyloepiphyseal dysplasia tarda can be caused by mutations in the SEDL gene. This study describes an interesting novel mutation (IVS4+1A>G) located exactly at the rare noncanonical AT-AC consensus splicing donor point of SEDL, which regained the canonical GT-AG consensus splicing junction in addition to several other rarer noncanonical splice patterns. The mutation activated several cryptic splice sites and generated the production of seven erroneous splicing isoforms, which we confirmed by sequencing of RT-PCR products and resequencing of cDNA clones. All the practical splice donors/acceptors were further assessed using FSPLICE 1.0 and SPL(M) Platforms to predict potential splice sites in genomic DNA. Subsequently, the expression levels of SEDL among the affected patients, carriers and controls were estimated using real-time quantitative PCR. Expression analyses showed that the expression levels of SEDL in both patients and carriers were decreased. Taken together, these results illustrated how disruption of the AT donor site in a rare AT-AC intron, leading to a canonical GT donor site, resulted in a multitude of aberrant transcripts, thus impairing exon definition. The unexpected splicing patterns resulting from the special mutation provide additional challenges and opportunities for understanding splicing mechanisms and specificity.

U12-dependent intron splicing in plants

U12-dependent (U12) introns have persisted in the genomes of plants since the ancestral divergence between plants and metazoans. These introns, which are rare, are found in a range of genes that include essential functions in DNA replication and RNA metabolism and are implicated in regulating the expression of their host genes. U12 introns are removed from pre-mRNAs by a U12 intron-specific spliceosome. Although this spliceosome shares many properties with the more abundant U2-dependent (U2) intron spliceosome, four of the five small nuclear RNAs (snRNAs) required for splicing are different and specific for the unique splicing of U12 introns. Evidence in plants so far indicates that splicing signals of plant U12 introns and their splicing machinery are similar to U12 intron splicing in other eukaryotes. In addition to the high conservation of splicing signals, plant U12 introns also retain unique characteristic features of plant U2 introns, such as UA-richness, which suggests a requirement for plant-specific components for both the U2 and U12 splicing reaction. This chapter compares U12 and U2 splicing and reviews what is known about plant U12 introns and their possible role in gene expression.

Novel mRNA isoforms of the sodium channels Na(v)1.2, Na(v)1.3 and Na(v)1.7 encode predicted two-domain, truncated proteins

The expression of voltage-gated sodium channels is regulated at multiple levels, and in this study we addressed the potential for alternative splicing of the Na(v)1.2, Na(v)1.3, Na(v)1.6 and Na(v)1.7 mRNAs. We isolated novel mRNA isoforms of Na(v)1.2 and Na(v)1.3 from adult mouse and rat dorsal root ganglia (DRG), Na(v)1.3 and Na(v)1.7 from adult mouse brain, and Na(v)1.7 from neonatal rat brain. These alternatively spliced isoforms introduce an additional exon (Na(v)1.2 exon 17A and topologically equivalent Na(v)1.7 exon 16A) or exon pair (Na(v)1.3 exons 17A and 17B) that contain an in-frame stop codon and result in predicted two-domain, truncated proteins. The mouse and rat orthologous exon sequences are highly conserved (94-100% identities), as are the paralogous Na(v)1.2 and Na(v)1.3 exons (93% identity in mouse) to which the Na(v)1.7 exon has only 60% identity. Previously, Na(v)1.3 mRNA has been shown to be upregulated in rat DRG following peripheral nerve injury, unlike the downregulation of all other sodium channel transcripts. Here we show that the expression of Na(v)1.3 mRNA containing exons 17A and 17B is unchanged in mouse following peripheral nerve injury (axotomy), whereas total Na(v)1.3 mRNA expression is upregulated by 33% (P=0.003), suggesting differential regulation of the alternatively spliced transcripts. The alternatively spliced rodent exon sequences are highly conserved in both the human and chicken genomes, with 77-89% and 72-76% identities to mouse, respectively. The widespread conservation of these sequences strongly suggests an additional level of regulation in the expression of these channels, that is also tissue-specific.

Ribonucleoprotein assembly defects correlate with spinal muscular atrophy severity and preferentially affect a subset of spliceosomal snRNPs

Spinal muscular atrophy (SMA) is a motor neuron disease caused by reduced levels of the survival motor neuron (SMN) protein. SMN together with Gemins2-8 and unrip proteins form a macromolecular complex that functions in the assembly of small nuclear ribonucleoproteins (snRNPs) of both the major and the minor splicing pathways. It is not known whether the levels of spliceosomal snRNPs are decreased in SMA. Here we analyzed the consequence of SMN deficiency on snRNP metabolism in the spinal cord of mouse models of SMA with differing phenotypic severities. We demonstrate that the expression of a subset of Gemin proteins and snRNP assembly activity are dramatically reduced in the spinal cord of severe SMA mice. Comparative analysis of different tissues highlights a similar decrease in SMN levels and a strong impairment of snRNP assembly in tissues of severe SMA mice, although the defect appears smaller in kidney than in neural tissue. We further show that the extent of reduction in both Gemin proteins expression and snRNP assembly activity in the spinal cord of SMA mice correlates with disease severity. Remarkably, defective SMN complex function in snRNP assembly causes a significant decrease in the levels of a subset of snRNPs and preferentially affects the accumulation of U11 snRNP--a component of the minor spliceosome--in tissues of severe SMA mice. Thus, impairment of a ubiquitous function of SMN changes the snRNP profile of SMA tissues by unevenly altering the normal proportion of endogenous snRNPs. These findings are consistent with the hypothesis that SMN deficiency affects the splicing machinery and in particular the minor splicing pathway of a rare class of introns in SMA.

Characterization of Bruton's tyrosine kinase mutations in Mexican patients with X-linked agammaglobulinemia

X-linked agammaglobulinemia (XLA) is a humoral primary immunodeficiency in which affected patients have very low levels of peripheral B cells and a profound deficiency of all immunoglobulin isotypes. Mutations in the gene encoding for Bruton's tyrosine kinase (Btk) are responsible for most of the agammaglobulinemia. In this work, 14 Btk mutations responsible of causing XLA are described; eight of which are novel and six are mutations previously reported. Seven of the mutations were due to deletions and insertions of exons and introns, respectively, which suggest splicing defects. The others were missense mutations, five of which affect arginine residues and have been described, and two new which affect leucine and glutamine residues (L111P and E605G). Most of these mutations were located at the kinase domain of Btk and, less frequently, they were found in PH and SH2 domains. Protein expression was also affected since most of the patients did not express or express very low Btk.

Characterization and evolution of the novel gene family FAM90A in primates originated by multiple duplication and rearrangement events

Genomic plasticity of human chromosome 8p23.1 region is highly influenced by two groups of complex segmental duplications (SDs), termed REPD and REPP, that mediate different kinds of rearrangements. Part of the difficulty to explain the wide range of phenotypes associated with 8p23.1 rearrangements is that REPP and REPD are not yet well characterized, probably due to their polymorphic status. Here, we describe a novel primate-specific gene family, named FAM90A (family with sequence similarity 90), found within these SDs. According to the current human reference sequence assembly, the FAM90A family includes 24 members along 8p23.1 region plus a single member on chromosome 12p13.31, showing copy number variation (CNV) between individuals. These genes can be classified into subfamilies I and II, which differ in their upstream and 5'-untranslated region sequences, but both share the same open reading frame and are ubiquitously expressed. Sequence analysis and comparative fluorescence in situ hybridization studies showed that FAM90A subfamily II suffered a big expansion in the hominoid lineage, whereas subfamily I members were likely generated sometime around the divergence of orangutan and African great apes by a fusion process. In addition, the analysis of the Ka/Ks ratios provides evidence of functional constraint of some FAM90A genes in all species. The characterization of the FAM90A gene family contributes to a better understanding of the structural polymorphism of the human 8p23.1 region and constitutes a good example of how SDs, CNVs and rearrangements within themselves can promote the formation of new gene sequences with potential functional consequences.

Nucleases of the metallo-beta-lactamase family and their role in DNA and RNA metabolism

Proteins of the metallo-beta-lactamase family with either demonstrated or predicted nuclease activity have been identified in a number of organisms ranging from bacteria to humans and has been shown to be important constituents of cellular metabolism. Nucleases of this family are believed to utilize a zinc-dependent mechanism in catalysis and function as 5' to 3' exonucleases and or endonucleases in such processes as 3' end processing of RNA precursors, DNA repair, V(D)J recombination, and telomere maintenance. Examples of metallo-beta-lactamase nucleases include CPSF-73, a known component of the cleavage/polyadenylation machinery, which functions as the endonuclease in 3' end formation of both polyadenylated and histone mRNAs, and Artemis that opens DNA hairpins during V(D)J recombination. Mutations in two metallo-beta-lactamase nucleases have been implicated in human diseases: tRNase Z required for 3' processing of tRNA precursors has been linked to the familial form of prostate cancer, whereas inactivation of Artemis causes severe combined immunodeficiency (SCID). There is also a group of as yet uncharacterized proteins of this family in bacteria and archaea that based on sequence similarity to CPSF-73 are predicted to function as nucleases in RNA metabolism. This article reviews the cellular roles of nucleases of the metallo-beta-lactamase family and the recent advances in studying these proteins.

Spliceosomal small nuclear RNA genes in 11 insect genomes

The removal of introns from the primary transcripts of protein-coding genes is accomplished by the spliceosome, a large macromolecular complex of which small nuclear RNAs (snRNAs) are crucial components. Following the recent sequencing of the honeybee (Apis mellifera) genome, we used various computational methods, ranging from sequence similarity search to RNA secondary structure prediction, to search for putative snRNA genes (including their promoters) and to examine their pattern of conservation among 11 available insect genomes (A. mellifera, Tribolium castaneum, Bombyx mori, Anopheles gambiae, Aedes aegypti, and six Drosophila species). We identified candidates for all nine spliceosomal snRNA genes in all the analyzed genomes. All the species contain a similar number of snRNA genes, with the exception of A. aegypti, whose genome contains more U1, U2, and U5 genes, and A. mellifera, whose genome contains fewer U2 and U5 genes. We found that snRNA genes are generally more closely related to homologs within the same genus than to those in other genera. Promoter regions for all spliceosomal snRNA genes within each insect species share similar sequence motifs that are likely to correspond to the PSEA (proximal sequence element A), the binding site for snRNA activating protein complex, but these promoter elements vary in sequence among the five insect families surveyed here. In contrast to the other insect species investigated, Dipteran genomes are characterized by a rapid evolution (or loss) of components of the U12 spliceosome and a striking loss of U12-type introns.

EMMA: an efficient massive mapping algorithm using improved approximate mapping filtering

Efficient massive mapping algorithm (EMMA), an algorithm on efficiently mapping massive cDNAs onto genomic sequences, has recently been developed. The process of mapping massive cDNAs onto genomic sequences has been improved using more approximate mapping filtering based on an enhanced suffix array coupled with a pruned fast hash table, algorithms of block alignment extensions, and k-longest paths. When compared with the classical BLAT software in this field, the computing of EMMA ranges from two to forty-one times faster under similar prediction precisions.

A pollen-specific polygalacturonase from lily is related to major grass pollen allergens

A pollen-specific gene from lily (Lilium longiflorum Thunb. cv. Snow Queen), designated LLP-PG, was characterized. Southern blots of lily genomic DNA indicated that LLP-PG is a member of a small gene family. A thorough sequence analysis revealed that the LLP-PG gene is interrupted by two introns and encodes a protein of 413 amino acids, with a calculated molecular mass of 44 kDa, and a pI of 8.1. Evaluation of the hydropathy profile showed that the protein has a hydrophobic segment at the N-terminus, indicating the presence of a putative signal peptide. A sequence similarity search showed a significant homology of the encoded protein to pollen polygalacturonases (PGs) from various plant species and to an important group (group 13) of grass pollen allergens. The LLP-PG transcript is pollen-specific and it accumulates only at the latest stage during pollen development, in the mature pollen. In contrast to other "late genes" LLP-PG transcript can neither be induced by abscisic acid (ABA) nor by dehydration. Immunoblot analyses of pollen protein extracts from lily, timothy grass and tobacco with IgG antibodies directed against LLP-PG and against the timothy grass pollen allergen, Phl p 13, indicated that lily LLP-PG shares surface-exposed epitopes with pollen PGs from monocotyledonous and dicotyledonous plants. Enzyme-linked immunosorbent assay (ELISA) analyses and inhibition ELISA assays with patients' IgE demonstrated a very low IgE reactivity of lily rLLP-PG and a lack of cross-reactivity between rLLP-PG and the timothy grass pollen allergen, rPhl p 13. These data demonstrated that despite the significant sequence homology and the conserved surface-exposed epitopes LLP-PG represents a low-allergenic member of pollen PGs.

The PKA-CREB system encoded by the honeybee genome

The cAMP-dependent kinase (PKA) plays a crucial part in long-term memory formation in the honeybee (Apis mellifera). One of the putative substrates of the PKA activity is the cAMP response element binding protein (CREB), a transcription factor in the bZIP protein family. We searched the honeybee genome to characterize genes from the CREB/CREM and the PKA families. We identified two genes that encode regulatory subunits and three genes encode catalytic subunits of PKA. Eight genes code for bZIP proteins, but only one gene was found that encodes a member of the CREB/CREM family. The phylogenetic relationship of these genes was analysed with their Drosophila and human counterparts.

Comprehensive splice-site analysis using comparative genomics

We have collected over half a million splice sites from five species-Homo sapiens, Mus musculus, Drosophila melanogaster, Caenorhabditis elegans and Arabidopsis thaliana-and classified them into four subtypes: U2-type GT-AG and GC-AG and U12-type GT-AG and AT-AC. We have also found new examples of rare splice-site categories, such as U12-type introns without canonical borders, and U2-dependent AT-AC introns. The splice-site sequences and several tools to explore them are available on a public website (SpliceRack). For the U12-type introns, we find several features conserved across species, as well as a clustering of these introns on genes. Using the information content of the splice-site motifs, and the phylogenetic distance between them, we identify: (i) a higher degree of conservation in the exonic portion of the U2-type splice sites in more complex organisms; (ii) conservation of exonic nucleotides for U12-type splice sites; (iii) divergent evolution of C.elegans 3' splice sites (3'ss) and (iv) distinct evolutionary histories of 5' and 3'ss. Our study proves that the identification of broad patterns in naturally-occurring splice sites, through the analysis of genomic datasets, provides mechanistic and evolutionary insights into pre-mRNA splicing.

Cloning and characterization of rat importin 9: implication for its neuronal function

We describe the structure of the rat importin 9 gene, together with its transcripts and the encoded protein with its putative functional domains. The importin 9 gene contains 24 exons in a genomic region spanning >52,000 bp. It is transcribed into two mRNAs, generated by means of alternative polyadenylation site usage arranged in tandem. Both transcripts possess the same noncanonical polyadenylation signal (AGUAAA) in rat, this hexamer being conserved in all vertebrates examined. Additionally, intron 8 is bordered by AT-AC dinucleotides. Importin 9 is expressed throughout adult rat tissues, but the 114-kDa Importin 9 protein was detected only in the brain. The localization of the Importin 9 protein was examined by immunohistochemistry in both adult rat tissues and primary hippocampal cell cultures. The strongest labeling was detected in vivo in areas populated by neurons in high density and also in the dendritic processes emanating from these cells. This protein was clearly concentrated in the nuclei of these cells, although their cytoplasms too were heavily labeled. Strong cytoplasmic and very strong nuclear staining was found in a vast majority of the cells with neuronal morphology in vitro. Cultured cells with glial morphology generally exhibited a weaker cytoplasmic labeling. In these cells, the signal decorated the nuclear envelope without nuclear staining and gradually dwindled toward the cell periphery. These results hint at the cell- or tissue-type specific functions of this type of importin protein.

Free energy landscapes of RNA/RNA complexes: with applications to snRNA complexes in spliceosomes

We develop a statistical mechanical model for RNA/RNA complexes with both intramolecular and intermolecular interactions. As an application of the model, we compute the free energy landscapes, which give the full distribution for all the possible conformations, for U4/U6 and U2/U6 in major spliceosome and U4atac/U6atac and U12/U6atac in minor spliceosome. Different snRNA experiments found contrasting structures, our free energy landscape theory shows why these structures emerge and how they compete with each other. For yeast U2/U6, the model predicts that the two distinct experimental structures, the four-helix junction structure and the helix Ib-containing structure, can actually coexist and specifically compete with each other. In addition, the energy landscapes suggest possible mechanisms for the conformational switches in splicing. For instance, our calculation shows that coaxial stacking is essential for stabilizing the four-helix junction in yeast U2/U6. Therefore, inhibition of the coaxial stacking possibly by protein-binding may activate the conformational switch from the four-helix junction to the helix Ib-containing structure. Moreover, the change of the energy landscape shape gives information about the conformational changes. We find multiple (native-like and misfolded) intermediates formed through base-pairing rearrangements in snRNA complexes. For example, the unfolding of the U2/U6 undergoes a transition to a misfolded state which is functional, while in the unfolding of U12/U6atac, the functional helix Ib is found to be the last one to unfold and is thus the most stable structural component. Furthermore, the energy landscape gives the stabilities of all the possible (functional) intermediates and such information is directly related to splicing efficiency.

The Rice Mutant dwarf bamboo shoot 1: A Leaky Mutant of the NACK-type Kinesin-like Gene Can Initiate Organ Primordia but not Organ Development

That plant dwarfism is caused by hormonal defects related to gibberellin and brassinosteroid has been well documented. Other contributing elements, however, have not been elucidated. Here, we report on one of the most severe dwarf mutants of rice, dwarf bamboo shoot 1 (dbs1). Most mutant plants died within 1 month after sowing, but a few (5.2%) survived and grew. Vacuolation enlarged cells in the leaf primordia and seminal root before abortion, which disrupted the organized cell files in these organs. Relative to the severe defects in shoot and root growth, the overall structure of the dbs1 embryo was almost normal. Similarly, initiation and organogenesis of the leaf primordia at the shoot apical meristem and those of the lateral root primordia at the root elongation zone occurred normally. These observations suggest that DBS1 is involved in the growth and development of organs but not in organ initiation or organogenesis. Positional cloning of DBS1 revealed that it encoded a NACK-type kinesin-like protein (OsNACK), homologous to the essential components of a mitogen-activated protein kinase cascade during plant cytokinesis. A BLAST search indicated that DBS1 was the only gene encoding the OsNACK-type protein in the rice genome, and the dbs1 mutant produced only small amounts of the translatable DBS1 mRNA. Thus, we conclude that the dbs1 mutation causes a severe defect in DBS1 function but does not completely shut it down. We discuss the leaky phenotype of dbs1 under the restricted functioning of OsNACK.

Splicing of a rare class of introns by the U12-dependent spliceosome

Pre-mRNA splicing is catalyzed by two unique spliceosomes, designated U2- or U12-dependent. In contrast to the well-characterized U2-dependent spliceosome, much remains to be learned about the less abundant U12-type spliceosome. This review focuses on recent advances in elucidating the structure and function of the minor U12-dependent spliceosome. Interesting similarities and differences between the U12- and U2-dependent spliceosomes are also highlighted.

Identification and analysis of U5 snRNA variants in Drosophila

Distinct isoforms of spliceosomal RNAs may be involved in regulating pre-messenger RNA splicing in eukaryotic cells. During a large-scale effort to identify small noncoding RNAs in Drosophila, we isolated a U5 snRNA-like molecule containing a 5' segment identical to that of the canonical (major) U5 snRNA but with a variant Sm binding site and a distinct 3' hairpin sequence. Based on this finding, another six similar U5 snRNA-like sequences were identified within the Drosophila genome by sequence similarity to the invariant loop in the 5' half of U5. Interestingly, although all of these variants are expressed in vivo, each shows a distinct temporal expression profile during Drosophila development, and one is expressed primarily in fly heads. The presence of these U5 snRNA variants within RNP particles suggests their role in splicing and implies a possible connection to regulation of developmental and tissue-specific gene expression.

Identification and characterization of multiple isoforms of a murine and human tumor suppressor, patched, having distinct first exons

Mutations in mouse and human patched (PTCH) genes are associated with birth defects and cancer. PTCH, a 12-pass transmembrane protein, is a receptor for Sonic hedgehog (Shh) signaling proteins. Shh proteins activate transcription of target genes, including PTCH, via GLI transcription factors. Here we identified seven and five isoforms of human and mouse PTCH mRNA, respectively, which are generated by the complex alternative use of five exons as the first exon (exons 1a to 1e in the 5'-to-3' order). Although expression profiles of these isoforms were highly variable among human tissues, three of them, PTCHa, PTCHb, and PTCHd, were predominantly expressed in most tissues, PTCHd being most ubiquitous. In contrast, PTCHb was always predominant and reached a maximum at E10.5 during mouse development. These three mRNA isoforms encode three PTCH proteins with distinct N-termini, PTCH(L), PTCH(M), and PTCH(S). The expression of these three isoforms was regulated by GLI transcription factors, and at least two functional GLI-binding sequences were identified, one in exon 1a and the other between exon 1a and exon 1b. PTCH(L) and PTCH(M) were equally active in terms of suppressing GLI-mediated transcription and inducing apoptosis. PTCH(S) protein (encoded by PTCHd), lacking the first transmembrane domain, was more unstable than the other two, resulting in a reduced activity. This study may shed light on the mechanism whereby a single PTCH gene plays a role in both tumor cell growth and embryonic development.

GObar: a gene ontology based analysis and visualization tool for gene sets

BACKGROUND

Microarray experiments, as well as other genomic analyses, often result in large gene sets containing up to several hundred genes. The biological significance of such sets of genes is, usually, not readily apparent. Identification of the functions of the genes in the set can help highlight features of interest. The Gene Ontology Consortium 1 has annotated genes in several model organisms using a controlled vocabulary of terms and placed the terms on a Gene Ontology (GO), which comprises three disjoint hierarchies for Molecular functions, Biological processes and Cellular locations. The annotations can be used to identify functions that are enriched in the set, but this analysis can be misleading since the underlying distribution of genes among various functions is not uniform. For example, a large number of genes in a set might be kinases just because the genome contains many kinases.

RESULTS

We use the Gene Ontology hierarchy and the annotations to pick significant functions and pathways by comparing the distribution of functions in a given gene list against the distribution of all the genes in the genome, using the hypergeometric distribution to assign probabilities. GObar is a web-based visualizer that implements this algorithm. The public website for GObar 2 can analyse gene lists from the yeast (S. cervisiae), fly (D. Melanogaster), mouse (M. musculus) and human (H. sapiens) genomes. It also allows visualization of the GO tree, as well as placement of a single gene on the GO hierarchy. We analyse a gene list from a genomic study of pre-mRNA splicing to demonstrate the utility of GObar.

CONCLUSION

GObar is freely available as a web-based tool at http://katahdin.cshl.org:9331/GO2 and can help analyze and visualize gene lists from genomic analyses.

Different splicing defects lead to differential effects downstream of the lipid and protein phosphatase activities of PTEN

PTEN, encoding a dual phosphatase tumor suppressor, is mutated in 85 and 65% of individuals with Cowden syndrome (CS) and Bannayan-Riley-Ruvalcaba syndrome (BRRS), respectively. Approximately 23 germline mutations in putative splice sites have been published, but resulting downstream outcome data are limited. We determined splicing defects in PTEN in 40 germline PTEN mutation positive cases and 33 mutation negative cases with classic CS, BRRS and CS- or BRRS-like features. Altered splicing was observed in 4/40 mutation positive probands and 2/33 mutation negative probands. We then sought to characterize the transcriptional and biochemical outcomes of the five distinct splice-site mutations, which led to the skipping of exon 3, 4 or 6. Two mutation negative BRRS patients also showed exon 3 skipping, and later, genomic sequencing revealed a mutation deep in intron 2. The splice-site mutations leading to the deletions of exon 3, 4 or 6 resulted in reduced dual phosphatase activities of PTEN. Deletion of exon 4 was associated with severely reduced lipid phosphatase activity, whereas exon 3 skipping resulted in markedly reduced protein phosphatase activity. In addition, exon 3 deleted transcript and protein were stable and localized to the nucleus more efficiently than the wild-type PTEN. In contrast, exon 4 skipping resulted in unstable transcripts and severely truncated unstable PTEN protein lacking its phosphatase domain. We have not only described for the first time, the effect of a deep intronic/branch-site mutation on exon skipping in PTEN but also found that different splice-site mutations resulting in the deletion of different exons lead to distinct outcomes.

Evolutionary conservation of minor U12-type spliceosome between plants and humans

Splicing of rare, U12-type or AT-AC introns is mediated by a distinct spliceosome that assembles from U11, U12, U4atac, U6atac, and U5 snRNPs. Although in human cells the protein composition of minor and major snRNPs is similar, differences, particularly in U11 and U12 snRNPs, have been recently described. We have identified an Arabidopsis U11 snRNP-specific 35K protein as an interacting partner of an RS-domain-containing cyclophilin. By using a transient expression system in Arabidopsis protoplasts, we show that the 35K protein incorporates into snRNP. Oligo affinity selection and glycerol gradient centrifugation revealed that the Arabidopsis 35K protein is present in monomeric U11 snRNP and in U11/U12-di snRNP. The interaction of the 35K protein with Arabidopsis SR proteins together with its strong sequence similarity to U1-70K suggests that its function in splicing of minor introns is analogous to that of U1-70K. Analysis of Arabidopsis and Oryza sativa genome sequences revealed that all U11/U12-di-snRNP-specific proteins are conserved in dicot and monocot plants. In addition, we have identified an Arabidopsis gene encoding the homolog of U4atac snRNA and a second Arabidopsis gene encoding U6atac snRNA. Secondary structure predictions indicate that the Arabidopsis U4atac is able to form dimeric complexes with both Arabidopsis U6atac snRNAs. As revealed by RNaseA/T1 protection assay, the U4atac snRNA gene is expressed as an ~160-nt RNA, whereas the second U6atac snRNA gene seems to be a pseudogene. Taken together, our data indicate that recognition and splicing of minor, AT-AC introns in plants is highly similar to that in humans.

The mouse vesicular inhibitory amino acid transporter gene: Expression during embryogenesis, analysis of its core promoter in neural stem cells and a reconsideration of its alternate splicing

The vesicular inhibitory amino acid transporter, VIAAT (also known as vesicular GABA transporter VGAT) transports GABA or glycine into synaptic vesicles. To initiate an analysis of the expression and regulation of VIAAT during neurogenesis we have cloned and characterized the mouse Viaat gene. We find that the mouse Viaat coding sequence is encoded by two exons spanning 5.3 kb. A survey of expression by whole mount in situ hybridization of mouse embryos indicates that Viaat is activated early in neuron differentiation and is expressed widely within the developing CNS; however, we did not detect expression in the superficial non-neural structures that express the GABA synthase Gad1. Analysis of the Viaat promoter indicates that a minimal promoter region containing a CG rich sequence is sufficient for efficient expression in neural stem and precursor cells. Our analysis of the Viaat sequence and splicing does not support the existence of two Viaat isoforms as previously proposed [Ebihara et al., Brain Res. Mol Brain Res. 110 (2003), 126-139]. Instead, the alternative isoform Viaat-a appears to be due to PCR artifacts that have occurred independently in multiple labs.

Comparative genomics and diversifying selection of the clustered vertebrate protocadherin genes

To explain the mechanism for specifying diverse neuronal connections in the brain, Sperry proposed that individual cells carry chemoaffinity tags on their surfaces. The enormous complexity of these connections requires a tremendous diversity of cell-surface proteins. A large number of neural transmembrane protocadherin (Pcdh) proteins is encoded by three closely linked human and mouse gene clusters (alpha, beta, and gamma). To gain insight into Pcdh evolution, I performed comprehensive comparative cDNA and genomic DNA analyses for the three clusters in the chimpanzee, rat, and zebrafish genomes. I found that there are species-specific duplications in vertebrate Pcdh genes and that additional diversity is generated through alternative splicing within the zebrafish "variable" and "constant" regions. Moreover, different codons (sites) in the mammalian Pcdh ectodomains (ECs) are under diversifying selection, with some under diversity-enhancing positive Darwinian selection and others, including calcium-binding sites, under strong purifying selection. Interestingly, almost all positively selected codon positions are located on the surface of ECs 2 and 3. These diversified residues likely play an important role in combinatorial interactions of Pcdh proteins, which could provide the staggering diversity required for neuronal connections in the brain. These results also suggest that adaptive selection is an additional evolutionary factor for increasing Pcdh diversity.

An LKB1 AT-AC intron mutation causes Peutz-Jeghers syndrome via splicing at noncanonical cryptic splice sites

Peutz-Jeghers syndrome (PJS) is an autosomal dominant disorder associated with gastrointestinal polyposis and an increased cancer risk. PJS is caused by germline mutations in the tumor suppressor gene LKB1. One such mutation, IVS2+1A>G, alters the second intron 5' splice site, which has sequence features of a U12-type AT-AC intron. We report that in patients, LKB1 RNA splicing occurs from the mutated 5' splice site to several cryptic, noncanonical 3' splice sites immediately adjacent to the normal 3' splice site. In vitro splicing analysis demonstrates that this aberrant splicing is mediated by the U12-dependent spliceosome. The results indicate that the minor spliceosome can use a variety of 3' splice site sequences to pair to a given 5' splice site, albeit with tight constraints for maintaining the 3' splice site position. The unusual splicing defect associated with this PJS-causing mutation uncovers differences in splice-site recognition between the major and minor pre-mRNA splicing pathways.

Regulation of alternative RNA splicing by exon definition and exon sequences in viral and mammalian gene expression

Intron removal from a pre-mRNA by RNA splicing was once thought to be controlled mainly by intron splicing signals. However, viral and other eukaryotic RNA exon sequences have recently been found to regulate RNA splicing, polyadenylation, export, and nonsense-mediated RNA decay in addition to their coding function. Regulation of alternative RNA splicing by exon sequences is largely attributable to the presence of two major cis-acting elements in the regulated exons, the exonic splicing enhancer (ESE) and the suppressor or silencer (ESS). Two types of ESEs have been verified from more than 50 genes or exons: purine-rich ESEs, which are the more common, and non-purine-rich ESEs. In contrast, the sequences of ESSs identified in approximately 20 genes or exons are highly diverse and show little similarity to each other. Through interactions with cellular splicing factors, an ESE or ESS determines whether or not a regulated splice site, usually an upstream 3' splice site, will be used for RNA splicing. However, how these elements function precisely in selecting a regulated splice site is only partially understood. The balance between positive and negative regulation of splice site selection likely depends on the cis-element's identity and changes in cellular splicing factors under physiological or pathological conditions.

Analysis of the Set of GABAA Receptor Genes in the Human Genome

The genes of the ionotropic gamma-aminobutyric acid receptor (GABR) subunits have shown an unusual chromosomal clustering, but only now can this be fully specified by analyses of the human genome. We have characterized the genes encoding the 18 known human GABR subunits, plus one now located here, for their precise locations, sizes, and exon/intron structures. Clusters of 17 of the 19, distributed between five chromosomes, are specified in detail, and their possible significance is considered. By applying search algorithms designed to recognize sequences of all known GABR-type subunits in species from man down to nematodes, we found no new GABR subunit is detectable in the human genome. However, the sequence of the human orthologue of the rat GABR rho3 receptor subunit was uncovered by these algorithms, and its gene could be analyzed. Consistent with those search results, orthologues of the beta4 and gamma4 subunits from the chicken, not cloned from mammals, were not detectable in the human genome by specific searches for them. The relationships are consistent with the mammalian subunit being derived from the beta line and epsilon from the gamma line, with mammalian loss of beta4 and gamma4. In their structures the human GABR genes show a basic pattern of nine coding exons, with six different genomic mechanisms for the alternative splicing found in various subunits. Additional noncoding exons occur for certain subunits, which can be regulatory. A dicysteine loop and its exon show remarkable constancy between all GABR subunits and species, of deduced functional significance.

Determinants of plant U12-dependent intron splicing efficiency

Factors affecting splicing of plant U12-dependent introns have been examined by extensive mutational analyses in an in vivo tobacco (Nicotiana tabacum) protoplast system using introns from three different Arabidopsis thaliana genes: CBP20, GSH2, and LD. The results provide evidence that splicing efficiency of plant U12 introns depends on a combination of factors, including UA content, exon bridging interactions between the U12 intron and flanking U2-dependent introns, and exon splicing enhancer sequences (ESEs). Unexpectedly, all three plant U12 introns required an adenosine at the upstream purine position in the branchpoint consensus UCCUURAUY. The exon upstream of the LD U12 intron is a major determinant of its higher level of splicing efficiency and potentially contains two ESE regions. These results suggest that in plants, U12 introns represent a level at which expression of their host genes can be regulated.

Regulation of alternative RNA splicing by exon definition and exon sequences in viral and mammalian gene expression

Intron removal from a pre-mRNA by RNA splicing was once thought to be controlled mainly by intron splicing signals. However, viral and other eukaryotic RNA exon sequences have recently been found to regulate RNA splicing, polyadenylation, export, and nonsense-mediated RNA decay in addition to their coding function. Regulation of alternative RNA splicing by exon sequences is largely attributable to the presence of two major cis-acting elements in the regulated exons, the exonic splicing enhancer (ESE) and the suppressor or silencer (ESS). Two types of ESEs have been verified from more than 50 genes or exons: purine-rich ESEs, which are the more common, and non-purine-rich ESEs. In contrast, the sequences of ESSs identified in approximately 20 genes or exons are highly diverse and show little similarity to each other. Through interactions with cellular splicing factors, an ESE or ESS determines whether or not a regulated splice site, usually an upstream 3' splice site, will be used for RNA splicing. However, how these elements function precisely in selecting a regulated splice site is only partially understood. The balance between positive and negative regulation of splice site selection likely depends on the cis-element's identity and changes in cellular splicing factors under physiological or pathological conditions.

Recycling of the U12-type spliceosome requires p110, a component of the U6atac snRNP

U12-dependent introns are spliced by the so-called minor spliceosome, requiring the U11, U12, and U4atac/U6atac snRNPs in addition to the U5 snRNP. We have recently identified U6-p110 (SART3) as a novel human recycling factor that is related to the yeast splicing factor Prp24. U6-p110 transiently associates with the U6 and U4/U6 snRNPs during the spliceosome cycle, regenerating functional U4/U6 snRNPs from singular U4 and U6 snRNPs. Here we investigated the involvement of U6-p110 in recycling of the U4atac/U6atac snRNP. In contrast to the major U6 and U4/U6 snRNPs, p110 is primarily associated with the U6atac snRNP but is almost undetectable in the U4atac/U6atac snRNP. Since p110 does not occur in U5 snRNA-containing complexes, it appears to be transiently associated with U6atac during the cycle of the minor spliceosome. The p110 binding site was mapped to U6 nucleotides 38 to 57 and U6atac nucleotides 10 to 30, which are highly conserved between these two functionally related snRNAs. With a U12-dependent in vitro splicing system, we demonstrate that p110 is required for recycling of the U4atac/U6atac snRNP.

U4 small nuclear RNA can function in both the major and minor spliceosomes

U4 small nuclear RNA (snRNA) and U6 snRNA form a base-paired di-snRNP complex that is essential for pre-mRNA splicing of the major class of metazoan nuclear introns. The functionally analogous but highly diverged U4atac and U6atac snRNAs form a similar complex that is involved in splicing of the minor class of introns. Previous results with mutants of U6atac in which a substructure was replaced by the analogous structure from U6 snRNA suggested that wild-type U4 snRNA might be able to interact productively with the mutant U6atac snRNA. Here we show that a mutant U4 snRNA designed to base pair with a mutant U6atac snRNA can activate U12-dependent splicing when coexpressed in an in vivo genetic suppression assay. This genetic interaction could also be demonstrated in an in vitro crosslinking assay. These results show that a U4/U6atac di-snRNP can correctly splice a U12-dependent intron and suggest that the specificity for spliceosome type resides in the U6 and U6atac small nuclear ribonucleoproteins. Further experiments suggest that expression of a mutant U4 snRNA that can bind to wild-type U6atac snRNA alters the specificity of some splice sites, providing an evolutionary rationale for maintaining two U4-like snRNAs.

The impact of splice isoforms on voltage-gated calcium channel alpha1 subunits

Semi-conserved exon boundaries in members of the CACNA1 gene family result in recurring pre-mRNA splicing patterns. The resulting variations in the encoded pore-forming subunit of the voltage-gated calcium channel affect functionally significant regions, such as the vicinity of the voltage-sensing S4 segments or the intracellular loops that are important for protein interaction. In addition to generating functional diversity, RNA splicing regulates the quantitative expression of other splice isoforms of the same gene by producing transcripts with premature stop codons which encode two-domain or three-domain channels. An overview of some of the known splice isoforms of the alpha(1) calcium channel subunits and their significance is given.

Identification, characterization and molecular phylogeny of U12-dependent introns in the Arabidopsis thaliana genome

U12-dependent introns are spliced by the minor U12-type spliceosome and occur in a variety of eukaryotic organisms, including Arabidopsis. In this study, a set of putative U12-dependent introns was compiled from a large collection of cDNA/EST- confirmed introns in the Arabidopsis thaliana genome by means of high-throughput bioinformatic analysis combined with manual scrutiny. A total of 165 U12-type introns were identified based upon stringent criteria. This number of sequences well exceeds the total number of U12-type introns previously reported for plants and allows a more thorough statistical analysis of U12-type signals. Of particular note is the discovery that the distance between the branch site adenosine and the acceptor site ranges from 10 to 39 nt, significantly longer than the previously postulated limit of 21 bp. Further analysis indicates that, in addition to the spacing constraint, the sequence context of the potential acceptor site may have an important role in 3' splice site selection. Several alternative splicing events involving U12-type introns were also captured in this study, providing evidence that U12-dependent acceptor sites can also be recognized by the U2-type spliceosome. Furthermore, phylogenetic analysis suggests that both U12-type AT-AC and U12-type GT-AG introns occurred in Na+/H+ antiporters in a progenitor of animals and plants.

Refined annotation of the Arabidopsis genome by complete expressed sequence tag mapping

Expressed sequence tags (ESTs) currently encompass more entries in the public databases than any other form of sequence data. Thus, EST data sets provide a vast resource for gene identification and expression profiling. We have mapped the complete set of 176,915 publicly available Arabidopsis EST sequences onto the Arabidopsis genome using GeneSeqer, a spliced alignment program incorporating sequence similarity and splice site scoring. About 96% of the available ESTs could be properly aligned with a genomic locus, with the remaining ESTs deriving from organelle genomes and non-Arabidopsis sources or displaying insufficient sequence quality for alignment. The mapping provides verified sets of EST clusters for evaluation of EST clustering programs. Analysis of the spliced alignments suggests corrections to current gene structure annotation and provides examples of alternative and non-canonical pre-mRNA splicing. All results of this study were parsed into a database and are accessible via a flexible Web interface at http://www.plantgdb.org/AtGDB/.

Sequence, gene structure, and expression pattern of CTNNBL1, a minor-class intron-containing gene--evidence for a role in apoptosis

We have identified and characterized a cDNA designated CTNNBL1 (catenin (cadherin-associated protein), beta-like 1) coding for a protein of 563 amino acids having predicted structural homology to beta-catenin and other armadillo (arm) family proteins. CTNNBL1 is expressed in multiple human tissues, and its sequence is conserved across widely divergent species. The human CTNNBL1 gene on chromosome 20q11.2 contains 16 exons spanning > 178 kb. Intron 4 is a minor-class intron bearing AT at the 5' splice site and AC at the 3' splice site. An acidic domain, as well as a putative bipartite nuclear localization signal, a nuclear export signal, a leucine-isoleucine zipper, and phosphorylation motifs are present in the protein sequence. Transient expression of CTNNBL1 in CHO cells results in localization to the nucleus and apoptosis. The rate of cell death was higher when cells were transfected with a carboxy-terminal fragment of CTNNBL1, suggesting that the apoptosis-inducing activity is a function of this region.

Domains of human U4atac snRNA required for U12-dependent splicing in vivo

U4atac snRNA forms a base-paired complex with U6atac snRNA. Both snRNAs are required for the splicing of the minor U12-dependent class of eukaryotic nuclear introns. We have developed a new genetic suppression assay to investigate the in vivo roles of several regions of U4atac snRNA in U12-dependent splicing. We show that both the stem I and stem II regions, which have been proposed to pair with U6atac snRNA, are required for in vivo splicing. Splicing activity also requires U4atac sequences in the 5' stem-loop element that bind a 15.5 kDa protein that also binds to a similar region of U4 snRNA. In contrast, mutations in the region immediately following the stem I interaction region, as well as a deletion of the distal portion of the 3' stem-loop element, were active for splicing. Complete deletion of the 3' stem-loop element abolished in vivo splicing function as did a mutation of the Sm protein binding site. These results show that the in vivo sequence requirements of U4atac snRNA are similar to those described previously for U4 snRNA using in vitro assays and provide experimental support for models of the U4atac/U6atac snRNA interaction.

Structure of the sodium channel gene SCN11A: evidence for intron-to-exon conversion model and implications for gene evolution

Exon/intron boundaries in the regions encoding the trans-membrane segments of voltage-gated Na channel genes are conserved, supporting their proposed evolution from a single domain channel, while the exons encoding the cytoplasmic loops are less conserved with their evolutionary heritage being less defined. SCN11A encodes the tetrodotoxin-resistant (TTX-R) sodium channel Nav1.9a/NaN, which is preferentially expressed in nociceptive primary sensory neurons of dorsal root ganglia (DRG) and trigeminal ganglia. SCN11A is localized to human chromosome 3 (3p21-24) close to the other TTX-R sodium channel genes SCN5A and SCN10A. An alternative transcript, Nav1.9b, has been detected in rat DRG and trigeminal ganglion. Nav1.9b is predicted to produce a truncated protein due to a frame-shift, which is introduced by the new sequence of exon 23c (E23c). In human and mouse SCN11A, divergent splicing signals prevent utilization of E23c. Unlike exons 5A/N in genes encoding TTX-sensitive sodium channels, which appear to have resulted from exon duplication, E23c might have evolved from the conversion of an intronic sequence. Although a functional role for Nav1.9b has yet to be established, intron-to-exon conversion may represent a mechanism for ion channels to acquire novel features.

An aminophospholipid translocase associated with body fat and type 2 diabetes phenotypes

OBJECTIVE

We have shown that a region on proximal mouse chromosome 7, near the pink-eyed (p) dilution locus, contains an ATPase (pfatp), a putative aminophospholipid translocase. Studies have suggested that this gene is a prime candidate for modulating body fat or involved in lipid metabolism in mouse and humans. Toward further analyses, our objective was to generate the complete genomic structures of mouse and human genes.

RESEARCH METHODS AND PROCEDURES

The genomic structure of mouse pfatp was deduced by comparing the full-length cDNA sequence with the genomic sequence derived from a mouse BAC. The human ortholog was identified from the National Center for Biotechnology Information database. Full-length cDNA was generated, and the corresponding genomic structure was deduced from the Human Genome Database.

RESULTS

Murine pfatp, and its human ortholog, PFATP, belong to class V of the third subfamily of P-type ATPases. The gene organization is strikingly similar in both organisms and all exon-intron junctions are conserved. A putative promoter region of PFATP contains a strong CpG island. The 5' untranslated regions of the two cDNAs have potential binding sites for multiple transcription factors, including Sp1, USF, AP1, and AP2, involved in adipogenesis and adipocyte metabolism.

DISCUSSION

We report the generation of the complete genomic structure of a novel aminophospholipid translocase in mice and humans. Because the exact biological role and the subsequent relevance of these ATPases to obesity and diabetes are unknown, these data help to delineate the role of these genes in lipid/adipocyte metabolism.

Branchpoint selection in the splicing of U12-dependent introns in vitro

In metazoans, splicing of introns from pre-mRNAs can occur by two pathways: the major U2-dependent or the minor U12-dependent pathways. Whereas the U2-dependent pathway has been well characterized, much about the U12-dependent pathway remains to be discovered. Most of the information regarding U12-type introns has come from in vitro studies of a very few known introns of this class. To expand our understanding of U12-type splicing, especially to test the hypothesis that the simple base-pairing mechanism between the intron and U12 snRNA defines the branchpoint of U12-dependent introns, additional in vitro splicing substrates were created from three putative U12-type introns: the third intron of the Xenopus RPL1 a gene (XRP), the sixth intron of the Xenopus TFIIS.oA gene (XTF), and the first intron of the human Sm E gene (SME). In vitro splicing in HeLa nuclear extract confirmed U12-dependent splicing of each of these introns. Surprisingly, branchpoint mapping of the XRP splicing intermediate shows use of the upstream rather than the downstream of two consecutive adenosines within the branchpoint sequence (BPS), contrary to the prediction based on alignment with the sixth intron of human P120, a U12-dependent intron whose branch site was previously determined. Also, in the SME intron, the position of the branchpoint A residue within the region base paired with U12 differs from that in P120 and XTF. Analysis of these three additional introns therefore rules out simple models for branchpoint selection by the U12-type spliceosome.

Characterization of the gene EPAC2: structure, chromosomal localization, tissue expression, and identification of the liver-specific isoform

The liver-specific protein cAMP-GEFII (also known as Epac2) belongs to a family of cyclic adenosine monophosphate (cAMP) binding proteins having guanine nucleotide exchange factor (GEF) activity (the cAMP-GEF family). Here we clone the gene EPAC2, encoding cAMP-GEFII, from a human liver cDNA library. Human EPAC2 has at least 31 exons and is mapped to human chromosome 2q31. Analyses by primer extension, reverse transcriptase-polymerase chain reaction, and in situ hybridization revealed the presence of three transcription start sites of liver-specific Epac2: two major sites located in exon 10 and a minor site in intron 9. The same translation start site is used in all three transcripts. Liver-specific cAMP-GEFII protein, which lacks the first cAMP-binding domain and the Dishevelled/Egl-10/Pleckstrin domain, was detected at 79 kDa by immunoblot analysis, confirming the presence of the short form of cAMP-GEFII in the liver. Liver-specific cAMP-GEFII also has GEF activity toward Rap1. These results demonstrate the presence of liver-specific cAMP-GEFII. Together with the previous finding that cAMP-GEFII is responsible for cAMP-dependent exocytosis in secretory cells, our study suggests that cAMP-GEFII may have a distinct role in liver.