Gene expression profiling of U12-type spliceosome mutant Drosophila reveals widespread changes in metabolic pathways

Background

The U12-type spliceosome is responsible for the removal of a subset of introns from eukaryotic mRNAs. U12-type introns are spliced less efficiently than normal U2-type introns, which suggests a rate-limiting role in gene expression. The Drosophila genome contains about 20 U12-type introns, many of them in essential genes, and the U12-type spliceosome has previously been shown to be essential in the fly.

Methodology/Principal Findings

We have used a Drosophila line with a P-element insertion in U6atac snRNA, an essential component of the U12-type spliceosome, to investigate the impact of U12-type introns on gene expression at the organismal level during fly development. This line exhibits progressive accumulation of unspliced U12-type introns during larval development and the death of larvae at the third instar stage. Surprisingly, microarray and RT-PCR analyses revealed that most genes containing U12-type introns showed only mild perturbations in the splicing of U12-type introns. In contrast, we detected widespread downstream effects on genes that do not contain U12-type introns, with genes related to various metabolic pathways constituting the largest group.

Conclusions/Significance

U12-type intron-containing genes exhibited variable gene-specific responses to the splicing defect, with some genes showing up- or downregulation, while most did not change significantly. The observed residual U12-type splicing activity could be explained with the mutant U6atac allele having a low level of catalytic activity. Detailed analysis of all genes suggested that a defect in the splicing of the U12-type intron of the mitochondrial prohibitin gene may be the primary cause of the various downstream effects detected in the microarray analysis.     

Evolutionary dynamics of U12-type spliceosomal introns

Background  

Many multicellular eukaryotes have two types of spliceosomes for the removal of introns from messenger RNA precursors. The major (U2) spliceosome processes the vast majority of introns, referred to as U2-type introns, while the minor (U12) spliceosome removes a small fraction (less than 0.5%) of introns, referred to as U12-type introns. U12-type introns have distinct sequence elements and usually occur together in genes with U2-type introns. A phylogenetic distribution of U12-type introns shows that the minor splicing pathway appeared very early in eukaryotic evolution and has been lost repeatedly.     

Alternative splicing and bioinformatic analysis of human U12-type introns

U12-type introns exist, albeit rarely, in a variety of multicellular organisms. Splicing of U12 intron-containing precursor mRNAs takes place in the U12-type spliceosome that is distinct from the major U2-type spliceosome. Due to incompatibility of these two spliceosomes, alternative splicing involving a U12-type intron may give rise to a relatively complicated impact on gene expression. We studied alternative U12-type intron splicing in an attempt to gain more mechanistic insights. First, we characterized mutually exclusive exon selection of the human JNK2 gene, which involves an unusual intron possessing the U12-type 5′ splice site and the U2-type 3′ splice site. We demonstrated that the long and evolutionary conserved polypyrimidine tract of this hybrid intron provides important signals for inclusion of its downstream alternative exon. In addition, we examined the effects of single nucleotide polymorphisms in the human WDFY1 U12-type intron on pre-mRNA splicing. These results provide mechanistic implications on splice-site selection of U12-type intron splicing. We finally discuss the potential effects of splicing of a U12-type intron with genetic defects or within a set of genes encoding RNA processing factors on global gene expression.     

A computational scan for U12-dependent introns in the human genome sequence

U12-dependent introns are found in small numbers in most eukaryotic genomes, but their scarcity makes accurate characterisation of their properties challenging. A computational search for U12-dependent introns was performed using the draft version of the human genome sequence. Human expressed sequences confirmed 404 U12-dependent introns within the human genome, a 6-fold increase over the total number of non-redundant U12-dependent introns previously identified in all genomes. Although most of these introns had AT-AC or GT-AG terminal dinucleotides, small numbers of introns with a surprising diversity of termini were found, suggesting that many of the non-canonical introns found in the human genome may be variants of U12-dependent introns and, thus, spliced by the minor spliceosome. Comparisons with U2-dependent introns revealed that the U12-dependent intron set lacks the ‘short intron’ peak characteristic of U2-dependent introns. Analysis of this U12-dependent intron set confirmed reports of a biased distribution of U12-dependent introns in the genome and allowed the identification of several alternative splicing events as well as a surprising number of apparent splicing errors. This new larger reference set of U12-dependent introns will serve as a resource for future studies of both the properties and evolution of the U12 spliceosome.     

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.     

Endogenous U2·U5·U6 snRNA complexes in S. pombe are intron lariat spliceosomes

Excision of introns from pre-mRNAs is mediated by the spliceosome, a multi-megadalton complex consisting of U1, U2, U4/U6, and U5 snRNPs plus scores of associated proteins. Spliceosome assembly and disassembly are highly dynamic processes involving multiple stable intermediates. In this study, we utilized a split TAP-tag approach for large-scale purification of an abundant endogenous U2·U5·U6 complex from Schizosaccharomyces pombe. RNAseq revealed this complex to largely contain excised introns, indicating that it is primarily ILS (intron lariat spliceosome) complexes. These endogenous ILS complexes are remarkably resistant to both high-salt and nuclease digestion. Mass spectrometry analysis identified 68, 45, and 43 proteins in low-salt-, high-salt-, and micrococcal nuclease-treated preps, respectively. The protein content of a S. pombe ILS complex strongly resembles that previously reported for human spliced product (P) and Saccharomyces cerevisiae ILS complexes assembled on single pre-mRNAs in vitro. However, the ATP-dependent RNA helicase Brr2 was either substoichiometric in low-salt preps or completely absent from high-salt and MNase preps. Because Brr2 facilitates spliceosome disassembly, its relative absence may explain why the ILS complex accumulates logarithmically growing cultures and the inability of S. pombe extracts to support in vitro splicing.     

Cwc21p promotes the second step conformation of the spliceosome and modulates 3′ splice site selection

Pre-mRNA splicing involves two transesterification steps catalyzed by the spliceosome. How RNA substrates are positioned in each step and the molecular rearrangements involved, remain obscure. Here, we show that mutations in PRP16, PRP8, SNU114 and the U5 snRNA that affect this process interact genetically with CWC21, that encodes the yeast orthologue of the human SR protein, SRm300/SRRM2. Our microarray analysis shows changes in 3′ splice site selection at elevated temperature in a subset of introns in cwc21Δ cells. Considering all the available data, we propose a role for Cwc21p positioning the 3′ splice site at the transition to the second step conformation of the spliceosome, mediated through its interactions with the U5 snRNP. This suggests a mechanism whereby SRm300/SRRM2, might influence splice site selection in human cells.     

Zinc-dependent cleavage in the catalytic core of the hammerhead ribozyme: evidence for a pH-dependent conformational change

We have characterized a novel Zn²⁺-catalyzed cleavage site between nucleotides C3 and U4 in the catalytic core of the hammerhead ribozyme. In contrast to previously described divalent metal-ion-dependent cleavage of RNA, U4 cleavage is only observed in the presence of Zn²⁺. This new cleavage site has an unusual pH dependence, in that U4 cleavage products are only observed above pH 7.9 and reach a maximum yield at about pH 8.5. These data, together with the fact that no metal ion-binding site is observed in proximity to the U4 cleavage site in either of the crystal structures, point toward a pH-dependent conformational change in the hammerhead ribozyme. We have described previously Zn²⁺-dependent cleavage between G8 and A9 in the hammerhead ribozyme and have discovered that U4 cleavage occurs only after A9 cleavage. To our knowledge, this is the first example of sequential cleavage events as a possible regulatory mechanism in ribozymes.     

Guard Cells Possess a Calcium-Dependent Protein Kinase That Phosphorylates the KAT1 Potassium Channel

Increasing evidence suggests that changes in cytosolic Ca²⁺ levels and phosphorylation play important roles in the regulation of stomatal aperture and as ion transporters of guard cells. However, protein kinases responsible for Ca²⁺ signaling in guard cells remain to be identified. Using biochemical approaches, we have identified a Ca²⁺-dependent protein kinase with a calmodulin-like domain (CDPK) in guard cell protoplasts of Vicia faba. Both autophosphorylation and catalytic activity of CDPK are Ca²⁺ dependent. CDPK exhibits a Ca²⁺-induced electrophoretic mobility shift and its Ca²⁺-dependent catalytic activity can be inhibited by the calmodulin antagonists trifluoperazine and N-(6-aminohexyl)-5-chloro-1-naphthalenesulfonamide. Antibodies to soybean CDPKα cross-react with CDPK. Micromolar Ca²⁺ concentrations stimulate phosphorylation of several proteins from guard cells; cyclosporin A, a specific inhibitor of the Ca²⁺-dependent protein phosphatase calcineurin enhances the Ca²⁺-dependent phosphorylation of several soluble proteins. CDPK from guard cells phosphorylates the K⁺ channel KAT1 protein in a Ca²⁺-dependent manner. These results suggest that CDPK may be an important component of Ca²⁺ signaling in guard cells.     

Ca2+ Regulates Reactive Oxygen Species Production and pH during Mechanosensing in Arabidopsis Roots

Mechanical stimulation of plants triggers a cytoplasmic Ca²⁺ increase that is thought to link the touch stimulus to appropriate growth responses. We found that in roots of Arabidopsis thaliana, external and endogenously generated mechanical forces consistently trigger rapid and transient increases in cytosolic Ca²⁺ and that the signatures of these Ca²⁺ transients are stimulus specific. Mechanical stimulation likewise elicited an apoplastic alkalinization and cytoplasmic acidification as well as apoplastic reactive oxygen species (ROS) production. These responses showed the same kinetics as mechanically induced Ca²⁺ transients and could be elicited in the absence of a mechanical stimulus by artificially increasing Ca²⁺ concentrations. Both pH changes and ROS production were inhibited by pretreatment with a Ca²⁺ channel blocker, which also inhibited mechanically induced elevations in cytosolic Ca²⁺. In trichoblasts of the Arabidopsis root hair defective2 mutant, which lacks a functional NADPH oxidase RBOH C, touch stimulation still triggered pH changes but not the local increase in ROS production seen in wild-type plants. Thus, mechanical stimulation likely elicits Ca²⁺-dependent activation of RBOH C, resulting in ROS production to the cell wall. This ROS production appears to be coordinated with intra- and extracellular pH changes through the same mechanically induced cytosolic Ca²⁺ transient.     

The Crystal Structure of PPIL1 Bound to Cyclosporine A Suggests a Binding Mode for a Linear Epitope of the SKIP Protein

Background

The removal of introns from pre-mRNA is carried out by a large macromolecular machine called the spliceosome. The peptidyl-prolyl cis/trans isomerase PPIL1 is a component of the human spliceosome and binds to the spliceosomal SKIP protein via a binding site distinct from its active site.

Principal Findings

Here, we have studied the PPIL1 protein and its interaction with SKIP biochemically and by X-ray crystallography. A minimal linear binding epitope derived from the SKIP protein could be determined using a peptide array. A 36-residue region of SKIP centred on an eight-residue epitope suffices to bind PPIL1 in pull-down experiments. The crystal structure of PPIL1 in complex with the inhibitor cyclosporine A (CsA) was obtained at a resolution of 1.15 Å and exhibited two bound Cd²⁺ ions that enabled SAD phasing. PPIL1 residues that have previously been implicated in binding of SKIP are involved in the coordination of Cd²⁺ ions in the present crystal structure. Employing the present crystal structure, the determined minimal binding epitope and previously published NMR data [1], a molecular docking study was performed. In the docked model of the PPIL1·SKIP interaction, a proline residue of SKIP is buried in a hydrophobic pocket of PPIL1. This hydrophobic contact is encircled by several hydrogen bonds between the SKIP peptide and PPIL1.

Conclusion

We characterized a short, linear epitope of SKIP that is sufficient to bind the PPIL1 protein. Our data indicate that this SKIP peptide could function in recruiting PPIL1 into the core of the spliceosome. We present a molecular model for the binding mode of SKIP to PPIL1 which emphasizes the versatility of cyclophilin-type PPIases to engage in additional interactions with other proteins apart from active site contacts despite their limited surface area.     

Determinant in Human Immunodeficiency Virus Type 1 for Efficient Replication under Cytokine-Induced CD4+ T-Helper 1 (Th1)- and Th2-Type Conditions

Cytokines are potent stimuli for CD4⁺-T-cell differentiation. Among them, interleukin-12 (IL-12) and IL-4 induce naive CD4⁺ T cells to become T-helper 1 (Th1) or Th2 cells, respectively. In this study we found that macrophage-tropic human immunodeficiency virus type 1 (HIV-1) strains replicated more efficiently in IL-12-induced Th1-type cultures derived from normal CD4⁺ T cells than did T-cell-line-tropic (T-tropic) strains. In contrast, T-tropic strains preferentially infected IL-4-induced Th2-type cultures derived from the same donor CD4⁺ T cells. Additional studies using chimeric viruses demonstrated that the V3 region of HIV-1 gp120 was the principal determinant for efficiency of replication. Cell fusion analysis showed that cells expressing envelope protein from a T-tropic strain effectively fused with IL-4-induced Th2-type culture cells. Flow cytometric analysis showed that the level of CCR5 expression was higher on IL-12-induced Th1-type culture cells, whereas CXCR4 was highly expressed on IL-4-induced Th2-type culture cells, although a low level of CXCR4 expression was observed on IL-12-induced Th1-type culture cells. These results indicate that HIV-1 isolates exhibit differences in the ability to infect CD4⁺-T-cell subsets such as Th1 or Th2 cells and that this difference may partly correlate with the expression of particular chemokine receptors on these cells. The findings suggest that immunological conditions are one of the factors responsible for inducing selection of HIV-1 strains.     

The Arabidopsis U12-type spliceosomal protein U11/U12-31K is involved in U12 intron splicing via RNA chaperone activity and affects plant development

U12 introns are removed from precursor-mRNA by a U12 intron-specific spliceosome that contains U11 and U12 small nuclear ribonucleoproteins. Although several proteins unique to the U12-type spliceosome have been identified, the manner by which they affect U12-dependent intron splicing as well as plant growth and development remain largely unknown. Here, we assessed the role of U11/U12-31K, a U12-type spliceosomal protein in Arabidopsis thaliana. T-DNA-tagged homozygote lines for U11/U12-31K could not be obtained, and heterozygote mutants were defective for seed maturation, indicating that U11/U12-31K is essential for the normal development of Arabidopsis. Knockdown of U11/U12-31K by artificial microRNA caused a defect in proper U12 intron splicing, resulting in abnormal stem growth and development of Arabidopsis. This defect in proper splicing was not restricted to specific U12-type introns, but most U12 intron splicing was influenced by U11/U12-31K. The stunted inflorescence stem growth was recovered by exogenously applied gibberellic acid (GA), but not by cytokinin, auxin, or brassinosteroid. GA metabolism-related genes were highly downregulated in U11/U12-31K knockdown plants. Importantly, U11/U12-31K was determined to harbor RNA chaperone activity. We propose that U11/U12-31K is an RNA chaperone that is indispensible for proper U12 intron splicing and for normal growth and development of plants.     

Human GC-AG alternative intron isoforms with weak donor sites show enhanced consensus at acceptor exon positions

It has been previously observed that the intrinsically weak variant GC donor sites, in order to be recognized by the U2-type spliceosome, possess strong consensus sequences maximized for base pair formation with U1 and U5/U6 snRNAs. However, variability in signal strength is a fundamental mechanism for splice site selection in alternative splicing. Here we report human alternative GC-AG introns (for the first time from any species), and show that while constitutive GC-AG introns do possess strong signals at their donor sites, a large subset of alternative GC-AG introns possess weak consensus sequences at their donor sites. Surprisingly, this subset of alternative isoforms shows strong consensus at acceptor exon positions 1 and 2. The improved consensus at the acceptor exon can facilitate a strong interaction with U5 snRNA, which tethers the two exons for ligation during the second step of splicing. Further, these isoforms nearly always possess alternative acceptor sites and exhibit particularly weak polypyrimidine tracts characteristic of AG-dependent introns. The acceptor exon nucleotides are part of the consensus required for the U2AF³⁵-mediated recognition of AG in such introns. Such improved consensus at acceptor exons is not found in either normal or alternative GT-AG introns having weak donor sites or weak polypyrimidine tracts. The changes probably reflect mechanisms that allow GC-AG alternative intron isoforms to cope with two conflicting requirements, namely an apparent need for differential splice strength to direct the choice of alternative sites and a need for improved donor signals to compensate for the central mismatch base pair (C-A) in the RNA duplex of U1 snRNA and the pre-mRNA. The other important findings include (i) one in every twenty alternative introns is a GC-AG intron, and (ii) three of every five observed GC-AG introns are alternative isoforms.     

AtCCX3 is an Arabidopsis endomembrane H+ -dependent K+ transporter

The Arabidopsis (Arabidopsis thaliana) cation calcium exchangers (CCXs) were recently identified as a subfamily of cation transporters; however, no plant CCXs have been functionally characterized. Here, we show that Arabidopsis AtCCX3 (At3g14070) and AtCCX4 (At1g54115) can suppress yeast mutants defective in Na⁺, K⁺, and Mn²⁺ transport. We also report high-capacity uptake of ⁸⁶Rb⁺ in tonoplast-enriched vesicles from yeast expressing AtCCX3. Cation competition studies showed inhibition of ⁸⁶Rb⁺ uptake in AtCCX3 cells by excess Na⁺, K⁺, and Mn²⁺. Functional epitope-tagged AtCCX3 fusion proteins were localized to endomembranes in plants and yeast. In Arabidopsis, AtCCX3 is primarily expressed in flowers, while AtCCX4 is expressed throughout the plant. Quantitative polymerase chain reaction showed that expression of AtCCX3 increased in plants treated with NaCl, KCl, and MnCl₂. Insertional mutant lines of AtCCX3 and AtCCX4 displayed no apparent growth defects; however, overexpression of AtCCX3 caused increased Na⁺ accumulation and increased ⁸⁶Rb⁺ transport. Uptake of ⁸⁶Rb⁺ increased in tonoplast-enriched membranes isolated from Arabidopsis lines expressing CCX3 driven by the cauliflower mosaic virus 35S promoter. Overexpression of AtCCX3 in tobacco (Nicotiana tabacum) produced lesions in the leaves, stunted growth, and resulted in the accumulation of higher levels of numerous cations. In summary, these findings suggest that AtCCX3 is an endomembrane-localized H⁺-dependent K⁺ transporter with apparent Na⁺ and Mn²⁺ transport properties distinct from those of previously characterized plant transporters.     

Degeneration of a homing endonuclease and its target sequence in a wild yeast strain

Mobile introns and inteins self-propagate by ‘homing’, a gene conversion process initiated by site-specific homing endonucleases. The VMA intein, which encodes the PI-SceI endonuclease in Saccharomyces cerevisiae, is present in several different yeast strains. Surprisingly, a wild wine yeast (DH1-1A) contains not only the intein⁺ allele, but also an inteinless allele that has not undergone gene conversion. To elucidate how these two alleles co-exist, we characterized the endonuclease encoded by the DH1-1A intein⁺ allele and the target site in the intein^– allele. Sequence analysis reveals seven mutations in the 31 bp recognition sequence, none of which occurs at positions that are individually critical for activity. However, binding and cleavage of the sequence by PI-SceI is reduced 10-fold compared to the S.cerevisiae target. The PI-SceI analog encoded by the DH1-1A intein⁺ allele contains 11 mutations at residues in the endonuclease and protein splicing domains. None affects protein splicing, but one, a R417Q substitution, accounts for most of the decrease in DNA cleavage and DNA binding activity of the DH1-1A protein. Loss of activity in the DH1-1A endonuclease and target site provides one explanation for co-existence of the intein⁺ and intein^– alleles.     

Shewanella spp. Genomic Evolution for a Cold Marine Lifestyle and In-Situ Explosive Biodegradation

Shewanella halifaxensis and Shewanella sediminis were among a few aquatic γ-proteobacteria that were psychrophiles and the first anaerobic bacteria that degraded hexahydro-1,3,5-trinitro-1,3,5-triazine (RDX). Although many mesophilic or psychrophilic strains of Shewanella and γ-proteobacteria were sequenced for their genomes, the genomic evolution pathways for temperature adaptation were poorly understood. On the other hand, the genes responsible for anaerobic RDX mineralization pathways remain unknown. To determine the unique genomic properties of bacteria responsible for both cold-adaptation and RDX degradation, the genomes of S. halifaxensis and S. sediminis were sequenced and compared with 108 other γ-proteobacteria including Shewanella that differ in temperature and Na⁺ requirements, as well as RDX degradation capability. Results showed that for coping with marine environments their genomes had extensively exchanged with deep sea bacterial genomes. Many genes for Na⁺-dependent nutrient transporters were recruited to use the high Na⁺ content as an energy source. For coping with low temperatures, these two strains as well as other psychrophilic strains of Shewanella and γ-proteobacteria were found to decrease their genome G+C content and proteome alanine, proline and arginine content (p-value <0.01) to increase protein structural flexibility. Compared to poorer RDX-degrading strains, S. halifaxensis and S. sediminis have more number of genes for cytochromes and other enzymes related to RDX metabolic pathways. Experimentally, one cytochrome was found induced in S. halifaxensis by RDX when the chemical was the sole terminal electron acceptor. The isolated protein degraded RDX by mono-denitration and was identified as a multiheme 52 kDa cytochrome using a proteomic approach. The present analyses provided the first insight into divergent genomic evolution of bacterial strains for adaptation to the specific cold marine conditions and to the degradation of the pollutant RDX. The present study also provided the first evidence for the involvement of a specific c-type cytochrome in anaerobic RDX metabolism.     

Computational screen for spliceosomal RNA genes aids in defining the phylogenetic distribution of major and minor spliceosomal components

The RNA molecules of the spliceosome are critical for specificity and catalysis during splicing of eukaryotic pre-mRNA. In order to examine the evolution and phylogenetic distribution of these RNAs, we analyzed 149 eukaryotic genomes representing a broad range of phylogenetic groups. RNAs were predicted using high-sensitivity local alignment methods and profile HMMs in combination with covariance models. The results provide the most comprehensive view so far of the phylogenetic distribution of spliceosomal RNAs. RNAs were predicted in many phylogenetic groups where these RNA were not previously reported. Examples are RNAs of the major (U2-type) spliceosome in all fungal lineages, in lower metazoa and many protozoa. We also identified the minor (U12-type) spliceosomal U11 and U6atac RNAs in Acanthamoeba castellanii, where U12 spliceosomal RNA as well as minor introns were reported recently. In addition, minor-spliceosome-specific RNAs were identified in a number of phylogenetic groups where previously such RNAs were not observed, including the nematode Trichinella spiralis, the slime mold Physarum polycephalum and the fungal lineages Zygomycota and Chytridiomycota. The detailed map of the distribution of the U12-type RNA genes supports an early origin of the minor spliceosome and points to a number of occasions during evolution where it was lost.     

U12-type spliceosomal introns of Insecta

Most of eukaryotic genes are interrupted by introns that need to be removed from pre-mRNAs before they can perform their function. This is done by complex machinery called spliceosome. Many eukaryotes possess two separate spliceosomal systems that process separate sets of introns. The major (U2) spliceosome removes majority of introns, while minute fraction of intron repertoire is processed by the minor (U12) spliceosome. These two populations of introns are called U2-type and U12-type, respectively. The latter fall into two subtypes based on the terminal dinucleotides. The minor spliceosomal system has been lost independently in some lineages, while in some others few U12-type introns persist. We investigated twenty insect genomes in order to better understand the evolutionary dynamics of U12-type introns. Our work confirms dramatic drop of U12-type introns in Diptera, leaving these genomes just with a handful cases. This is mostly the result of intron deletion, but in a number of dipteral cases, minor type introns were switched to a major type, as well. Insect genes that harbor U12-type introns belong to several functional categories among which proteins binding ions and nucleic acids are enriched and these few categories are also overrepresented among these genes that preserved minor type introns in Diptera.