Prevalent and distinct spliceosomal 3′-end processing mechanisms for fungal telomerase RNA

Telomerase RNA (TER) is an essential component of the telomerase ribonucleoprotein complex. The mechanism for TER 3′-end processing is highly divergent among different organisms. Here we report a unique spliceosome-mediated TER 3′-end cleavage mechanism in Neurospora crassa which is distinct from that found specifically in the fission yeast Schizosaccharomyces pombe. While the S. pombe TER intron contains the canonical 5′-splice site GUAUGU, the N. crassa TER intron contains a non-canonical 5′-splice site AUAAGU that alone prevents the second step of splicing and promotes spliceosomal cleavage. The unique N. crassa TER 5′-splice site sequence is evolutionarily conserved in TERs from Pezizomycotina and early branching Taphrinomycotina species. This suggests that the widespread and basal N. crassa-type spliceosomal cleavage mechanism is more ancestral than the S. pombe-type. The discovery of a prevalent, yet distinct, spliceosomal cleavage mechanism throughout diverse fungal clades furthers our understanding of TER evolution and non-coding RNA processing.


INTRODUCTION
The telomerase enzyme is responsible for the addition of telomeric DNA repeats onto the ends of eukaryotic chromosomes to maintain genome stability and cellular replicative capacity 1 . The core telomerase ribonucleoprotein (RNP) enzyme is composed of the catalytic telomerase reverse transcriptase (TERT) and telomerase RNA (TER) which contains a short region as the template for telomeric DNA repeat synthesis 2 . TERs are highly divergent across eukaryotic clades, differing immensely in nucleotide sequence, length and biogenesis pathway 3,4 . The diversity in TER biogenesis is evidenced by the different RNA polymerase (pol) machineries employed for transcription. While TER is transcribed by RNA pol II in most species, it is transcribed by RNA pol III in ciliates 5,6,7,8 . Furthermore, TERs from evolutionary separate groups of species associate with entirely unique groups of accessory proteins-the dyskerin complex and TCAB1 in vertebrates 9, 10 ; Est1, Est3 and Sm-ring in budding yeasts 11,12,13 ; the dyskerin complex and POT1 in the plant Arabidopsis 14,15 ; and p50 and p65 in ciliates 16 . These separate sets of TER accessory proteins are essential for localization, RNP assembly and TER 3′-end processing 17 .
The identification of the Schizosaccharomycetes pombe TER (SpoTER1) revealed a spliceosome-mediated cleavage mechanism for RNA 3′-end processing 18,19,20 . The nascent SpoTER1 transcript harbors a downstream intron followed by a second exon and a 3′poly(A) tail. The mature SpoTER1 3′-end is generated by a spliceosomal cleavage at the 5′splice site (5′-SS), removing the intron, second exon and poly(A) tail. This is in contrast to splicing-whereby the intervening intron sequence is excised and the neighboring exons are joined 21 . The SpoTER1 intron itself uncouples these two tightly associated splicing steps by inhibiting the second transesterification step at the 3′-splice site (3′-SS). The second transesterification step is abated by a strong branch point sequence (BPS), an extended distance from the BPS to the 3′-SS and a weak polypyrimidine (Py) in the SpoTER1 intron 19 .
Here we report in N. crassa a spliceosome-mediated mechanism for TER 3′-end processing that is distinct from the S. prombe-type. The N. crassa-type TER spliceosome-mediated 3′end cleavage relies exclusively on a unique 5′-SS AUAAGU, which alone is necessary and sufficient for inhibiting the second transesterification step in splicing. This is the first biologically functional instance of spliceosome-mediated cleavage relying exclusively on the 5′-SS. Additionally, the N. crassa-type TER intron is pervasive across Ascomycota, found universally in TERs from filamentous fungi and early branching fission yeasts. This preponderance and conservation among distinct, evolutionary distant and basal fungal species suggests this mechanism is more ancestral than the S. pombe-type. The discovery of this unique spliceosome-mediated cleavage mechanism provides new insights into TER biogenesis and evolution.

Neurospora crassa TER precursor harbors an intron
Neurospora crassa TER (NcrTER) was recently identified as a 2,049 nt RNA 22 . The 3′-end of the mature NcrTER was determined by rapid amplification of cDNA ends (RACE) from an RNA sample treated with poly(A) polymerase (Fig. 1a, lane 1). However, 3′-RACE with the RNA untreated with poly(A) polymerase generated two larger products (Fig. 1a, lane 2). The smaller of these two products was not readily visible by ethidium bromide staining, while discernable for gel extraction. Cloning and subsequent sequencing of the two PCR products revealed distinct NcrTER variants of approximately 400 and 470 bp in length. The 470 bp product harbored a 74 nt intron sequence, while the 400 bp product lacked this intron and had the two exons joined (Fig. 1b). The intron contained the conserved features of a canonical intron: 5′-SS, BPS and 3′-SS (Fig. 1c). The presence of an intron in NcrTER transcripts suggested that NcrTER undergoes spliceosome-mediated cleavage for 3′-end processing, similar to that discovered in the fission yeast SpoTER1 18 . Interestingly, the NcrTER intron 5′-SS AUAAGU is non-canonical, a stark contrast to the canonical 5′-SS GUAUGU sequence found in the SpoTER1 and mRNA introns 18 . We then determined the relative abundance of the three NcrTER variants by performing northern blot analysis of total RNA with probes targeting either the first or second exon. The 2,049 nt NcrTER was readily detected by northern blot, while the precursor and spliced variants were below the threshold for detection (Fig. 1d). This result suggests that the mature 2,049 nt NcrTER variant has far greater abundance. This is consistent with our initial identification of only the 2,049 nt NcrTER species in the purified telomerase holoenzyme 22 .

NcrTER-type intron is prevalent in Ascomycota fungal TERs
Analysis of 72 additional Pezizomycotina filamentous fungal TERs that we previously identified revealed the absolute conservation of the NcrTER-type intron in all species examined (Fig. 2a). All known Pezizomycotina TER introns contain the unique noncanonical 5′-SS AUAAGU sequence, the BPS sequence RCURAY that includes an invariant adenine residue as a putative branch point residue and the putative 3′-SS YAG ( Supplementary Fig. 1). The 3′-SS was experimentally verified by 3′-RACE for N. crassa, Aspergillus nidulans and Mycosphaerella graminicola ( Supplementary Fig. 1). The NcrTER 5′-SS AUAAGU is extremely conserved in filamentous fungal TERs and distinct from the canonical, yet less conserved 5′-SS GURNGU found in yeast TER introns from the fission yeast S. pombe as well as species from the budding yeast Candida and Hansenula genera 18,23,24 .
We then investigated whether the S. pombe-type or the N. crassa-type TER intron is more ancestral. S. pombe belongs to subphylum Taphrinomycotina that together with the other two subphyla Pezizomycotina and Saccharomycotina constitute the Ascomycota phylum. In addition to the TER that we previously cloned from a basal Taphrinomycotina species Saitoella complicata 22 , we identified two additional Taphrinomycotina TERs from Schizosaccharomyces octosporus and Schizosaccharomyces japonicus by a bioinformatics approach (see Methods). The 3′-ends and intron sequences of these newly identified TERs were experimentally determined by RACE. Surprisingly, there was a clear divide in the intron 5′-SS sequence among Taphrinomycotina TERs. The basal species S. complicate and S. japonicus harbored the N. crassa-type 5′-SS AUAAGU, while the later evolving species S. octosporus contained the S. pombe-type 5′-SS GURNGU (Fig. 2b). Thus, the more prevalent N. crassa-type intron appears to be more basal and thus ancestral than the S. pombe-type intron (Fig. 2c).

NcrTER-type 5′-SS promotes spliceosomal cleavage
The 5′-SS of the NcrTER and canonical mRNA introns differ with respect to their first and forth residues (Fig. 3a). Newman et al. demonstrated in yeast that an introduced G-to-A substitution at the first residue in an mRNA intron suppresses the second step, yet permits the first step of splicing 25 . We thus hypothesized that the 5′-SS is sufficient for NcrTER 3′-end processing. To test this, we constructed a minimal exogenous expression cassette of the NcrTER gene. This minimal NcrTER expression cassette spanned a 3′-portion (residues 1888-2049) of exon 1, the intron (residues 2050-2123), exon 2 (residues 2124-2200), and the 500 bp downstream flanking genomic sequence that contains a polyadenylation signal ( Fig. 3b and 3c). This expression cassette also contained an eGFP gene as an internal control for RNA expression and normalization. NcrTER and eGFP mRNA were transcribed independently under separate ccg-1 (clock-controlled gene-1) promoters and were fused to the ccg-1 5′-UTR sequence, generating ccg-1::NcrTER and ccg-1::eGFP, respectively. The ccg-1 promoter is a light induced, strong RNA pol II promoter commonly used for ectopic gene expression in the N. crassa model system 26 . Northern blot analysis was performed to determine the steady state level of the NcrTER exon 1 fragment and eGFP mRNA with probes targeting the fused ccg-1 5′-UTR sequence (Fig. 3b). The ccg-1::NcrTER construct generated a stable RNA product detectable by northern blot analysis (Fig. 3d, lane 3). This RNA product had the expected 3′-end corresponding to the mature NcrTER 3′-end, verified by 3′-RACE. Thus, the exogenous overexpressed minimal NcrTER 3′-fragment was processed by spliceosome-mediated cleavage, similar to the full-length endogenous NcrTER. RT-PCR was performed to detect the precursor and spliced RNA transcripts, which were below the detection limit of northern blot ( Fig. 1d and 2d). While the unspliced precursor and spliced NcrTER 3′-fragments were both detected, the precursor appeared to be much more abundant (Fig. 3e, lane 3). This is consistent with the RT-PCR results of the endogenous NcrTER precursor and spliced product (Fig. 1a, lane 2). Thus, our ccg-1::NcrTER expression cassette provides a minimal system for the functional analysis of spliceosome-mediated NcrTER 3′-end processing.
Mutations were introduced to the NcrTER intron in 5′-SS, BPS and 3′-SS to identify the element(s) necessary for NcrTER 3′-end processing (Fig. 3c). The 5′-SS_1 mutant with the first residue mutated from A to G reduced the level of mature NcrTER exon-1 RNA to 8 % of wild-type ( We then examined if the 5′-SS_1 mutant-which has the complete splicing restored-could be converted for S. pombe-type spliceosomal cleavage by extending the distance between the BPS and the 3′-SS ( Supplementary Fig. 3). We generated a double mutant 5′-SS_1a that, in addition to the 5′-SS_1 mutation, disrupted the 3′-SS (Fig. 3c). A triple mutant 5′-SS_1b further disrupted a downstream cryptic 3′-SS (Fig. 3c). Both 5′-SS_1a and 5′-SS_1b mutants effectively restored spliceosomal cleavage and accumulated mature RNA at levels comparable to the wild-type NcrTER intron (Fig. 3d, lanes 5 and 6), while reducing the splicing efficiency (Fig. 3e, lanes 5 and 6). RT-PCR and sequencing of the spliced RNA products from 5′-SS_1a and 5′-SS_1b mutants showed that, when the preferred 3′-SS was mutated, two different cryptic 3′-SSs were used instead with low efficiencies ( Supplementary Fig. 2). Thus, disruption of the N. crassa-type spliceosomal cleavage in the 5′-SS_1 mutant was rescued by extending the distance between the BPS and the 3′-SS to employ S. pombe-type spliceosomal cleavage.
To verify that the NcrTER 3′-end processing is indeed spliceosome-mediated and not through non-specific endo-or exo-nucleolytic RNA degradation, we generated two mutations, BPS_1 and No_intron. Mutant BPS_1 abolished splicing by mutating the essential adenosine residue of the BPS and generated exclusively precursor RNA, while mutant No_intron produced only the spliced RNA product by deleting the entire intron. As expected, BPS_1 and No_intron mutants failed to accumulate mature RNA, indicating that spliceosome-mediated processing is required to generate the mature 3′-end of NcrTER ( The NcrTER intron boundary 5′-AU/AG-3′ is chimeric with a 3′-SS similar to the canonical GU/AG mRNA introns and a 5′-SS similar to rare AU/AC mRNA introns 27,28,29 . We investigated whether the NcrTER intron can also be converted to an AU/AC intron for splicing. Mutations were introduced to convert the NcrTER intron 3′-SS UAG to UAC. Mutant 3′-SS_1 had the proximal 3′-SS UAG changed to UAC, while mutant 3′-SS_1a had an additional mutation converting a downstream cryptic 3′-SS UAG to UAC (Fig. 3c). Mutant 3′-SS_1 showed a slightly decreased level of mature RNA at 80% of the wild-type (Fig. 3f, lane 2) and effectively increased splicing efficiency via the introduced 3′-SS UAC (Fig. 3g, lane 2), confirmed by RT-PCR and sequencing of the spliced products ( Supplementary Fig. 2). Mutant 3′-SS_1a had a further decrease in spliceosomal cleavage to 60% of the wild-type (Fig. 3f, lane 3) and also a greater increase in splicing efficiency (Fig.  3g, lane 3). Thus, the NcrTER intron can be readily converted to an AU/AC splicing intron by mutating the proximal 3′-SS UAG to UAC.

5′-SS AUAAGU converts an intron for spliceosomal cleavage
The NcrTER intron can be converted to favor splicing by a simple A-to-G substitution at the first adenine reside of the 5′-SS (Fig. 3). To further confirm that the NcrTER 5′-SS AUAAGU is the only key determinate for spliceosomal cleavage, we performed a reciprocal experiment to convert an mRNA intron to the NcrTER-type intron for spliceosomal cleavage. We modified our ccg1::NcrTER expression cassette by replacing the NcrTER intron with an mRNA intron from the ccg-1 gene (Fig. 4a).
Mutations were introduced to the ccg-1 mRNA at 5′-SS, BPS and 3′-SS to identify the element(s) necessary for conversion to NcrTER-type spliceosomal cleavage (Fig. 4b). The replacement of the NcrTER intron in our minimal expression cassette with the ccg-1 mRNA intron resulted in severely reduced spliceosomal cleavage (Fig. 4c, lane 1) and dramatically increased splicing efficiency (Fig. 4d, lane 1). In mutant 5′-SS_3, replacing the ccg-1 5′-SS with the NcrTER 5′-SS, by two G-to-A substitutions at the first and fourth position of the 5′-SS, significantly increased the accumulation of mature NcrTER exon-1 3′-fragment (Fig. 4c, lane 2) and reduced splicing efficiency (Fig. 4d, lane 2). Thus, the NcrTER 5′-SS AUAAGU alone appears to be sufficient to convert an mRNA intron for spliceosomal cleavage. Interestingly, sequencing the residual spliced products of the 5′-SS_3 mutant showed that a new 3′-SS CAC was utilized, similar to the AU/AC mRNA splicing introns (Supplemental Fig. S4). In the presence of the 5′-SS AUAAGU, there is a shift from the 3′-SS YAG employed by GU/AG mRNA introns to YAC employed by AU/AC mRNA introns. To minimize the 3′-SS YAC, we generated two double mutants, 5′-SS_3a and _3b, that changed the 3′-SS CAC employed by mutant 5′-SS_3 to CGC or CAU, respectively (Fig. 4b). As expected, preventing AU/AC mRNA splicing for each of these mutants resulted in far greater accumulation of mature NcrTER than the mutant 5′-SS_3 (Fig. 4c, lanes 3 and 4) and similarly reduced splicing efficiency (Fig. 4d, lanes 3 and 4). The residual spliced products from the mutants 5′-SS_3a and 5′-SS_3b arose from utilizing the 3′-SS CAG or UAG, respectively (Supplementary Fig. 4). Thus, the 5′-SS_3a and 5′-SS_3b mutants employed the same 3′-SS YAG as the wild-type NcrTER intron.
We next investigated whether a single substitution at the first residue of the ccg-1 5′-SS would convert the ccg-1 mRNA intron to the NcrTER intron. In contrast to mutant 5′-SS_3 that had both the first and fourth residues substituted and accumulated mature NcrTER, mutant 5′-SS_4 with only the first residue substituted failed to accumulate mature NcrTER above the level of the wild-type ccg-1 intron (Fig. 4c, lane 5). The double mutations 5′-SS_4a and 5′-SS_4b eliminated the 3′-SS CAC to prevent AU/AC splicing. The northern blot showed that mutant 5′-SS_4a did not increase mature RNA level (Fig. 4c, lane 6). However, mutant 5′-SS_4b had only a slight increase in mature NcrTER accumulation, which is relatively insignificant compared to the 5′-SS_3 mutant (Fig. 4c, lane 7). Thus, the ccg-1 mRNA intron requires a 5′-SS identical to the NcrTER-type intron for spliceosomal cleavage. The requirement of G-to-A substitutions for the first and fourth residues in the ccg-1 mRNA intron may relay on the first adenine substitution disrupting snRNA U1 interaction and the fourth position substitution compensating for this disruption.
Since the ccg-1 mRNA intron can be easily converted into the NcrTER-type intron, we examined whether the ccg-1 mRNA intron can be converted for the S. pombe-type spliceosomal cleavage by extending the distance between the BPS and 3′-SS CAG. Mutant 3′-SS_2 had the proximal 3′-SS eliminated, while mutant Linker_1 introduced a 34-nt insertion between the BPS and the 3′-SS CAG (Fig. 4b). The northern blot and RT-PCR results showed that both mutants had a significantly increased level of mature NcrTER exon-1 RNA (Fig. 4c, lanes 8 and 9) and decreased splicing efficiency (Fig. 4d, lanes 8 and  9). Thus, the ccg-1 mRNA splicing intron can be readily converted to either the NcrTERtype or SpoTER1-type intron for spliceosomal cleavage. We also examined the importance of the BPS in promoting spliceosomal cleavage. Mutant BPS_3 converted the ccg-1 mRNA intron BPS to NcrTER-type intron BPS, while mutant BPS_4 converted it to the strong BPS of spoTER1 intron (Supplementary Fig. 3). These BPS mutants did not significantly increase the accumulation of mature NcrTER (Fig. 4e) or decrease the splicing efficiency (Fig. 4f) compared with the wild-type ccg-1 mRNA intron.

DISCUSSION
Telomerase biogenesis requires a 3′-end maturation step for TER, which is remarkably divergent across eukaryotic taxa and utilizes a wide array of distinct mechanisms. Within the phylum Ascomycota, certain budding yeast TERs employ specific transcription termination pathways and nucleolytic cleavage for 3′-end maturation 5,30 , while the fission yeast S. pombe and a few select budding yeast species employ spliceosomal cleavage 19 . Herein we have identified a second, yet vastly more widespread, spliceosome-mediated TER 3′-end processing mechanism initially in N. crassa and virtually ubiquitous throughout Ascomycete fungi. The N. crassa-type spliceosomal cleavage is mechanistically distinct from the S. pombe-type and is dependent on a unique 5′-SS AUAAGU sequence ( Fig. 2a and Supplementary Fig. 1). While critical for the S. pombe-type mechanism, the BPS and the distance to the 3′-SS are seemingly unimportant for the N. crassa-type. The prevalence of the N. crassa-type TER intron throughout all known Pezizomycotina fungal TERs and within the early diverging Taphrinomycotina fission yeast suggests that the N. crassa-type spliceosomal cleavage mechanism is ancestral and the S. pombe-type is later evolved.
We propose that the highly prevalent and basal N. crassa-type TER intron is ancestral and the S. pombe-type intron is later evolved by an A-to-G substitution of the 5′-SS together with strengthening the BPS, extendeding distance to the 3′-SS and weakening the Py tract. It would seem less likely that multiple TER-intron fusion events occurred during the short evolution time within each of the three Ascomycete subphyla. Instead, it is more plausible that two transition events converted the N. crassa-type to the S. pombe-type intron independently along the fission and budding yeasts lineages. The fission yeast transition event seems to have occurred adjacent to the emergence of the basal S. japonicas species which employs the ancestral N. crassa-type intron (Fig. 5). Similarly, select budding yeast species from the genera Candida and Hansenula harbor the S. pombe-type intron for TER 3′-end maturation (Figs. 2c and 5), while later-branching budding yeast species appear to have lost the spliceosomal cleavage intron entirely. These Saccharomyces species with intron-less TERs instead employ the Nrd1/Nab3-dependent transcription termination pathway for generating non-polyadenylated TER 30 . Two independent transition events from a more ancestral N. crassa-type TER intron to the S. pombe-type or the loss of TER intron explain the scattered presence of the S. pombe-type intron within Ascomycete fungal TERs and the uniform presence of the N. crassa-type intron.
While there is strong support for the ancestry of the N. crassa-type TER intron, there are two conceivable possibilities for the origin of the progenitor fungal TER intron. The progenitor fungal TER intron probably arose from a single gene fusion event between an intron-less fungal TER gene and either a common GU/AG or a rare AU/AC mRNA intron (Fig. 5). This newly intron-containing TER gene then evolved into the N. crassa-type and later into the S. pombe-type TER intron. Deriving the N. crassa-type TER intron from a GU/AG intron requires alterations in the 5′-SS, while adoption from an AU/AC intron necessitates a C-to-G substitution of the last residue in the 3′-SS. Our mutagenesis studies demonstrate that the NcrTER spliceosomal cleaving intron can be readily converted to either a canonical GU/AG or a rare AU/AC intron-permitting efficient splicing-with only single or double mutations necessary (Figs. 3 and 4). Previous studies in yeast have shown that an introduced G-to-A substitution of the first residue in the 5′-SS of a common GU/AG mRNA intron strongly inhibits the second step in splicing 25,31 . Arising from an AU/AC intron would negate a G-to-A transition substitution and instead require a 3′-SS mutation for efficient spliceosomal cleavage. Given the abundance of GU/AG mRNA intron in fungal genomes 32 , we favor the possibility that the N. crassa-type intron evolved from a GU/AG mRNA intron. While all three NcrTER variants reconstituted equal telomerase activity in vitro ( Supplementary Fig. 5), it is evident that a spliceosomal cleaving TER intron is evolutionarily selected for. Spliceosomal cleaving TER intron removes the poly(A) tail from TER, which presumably prevents un-desired cytoplasmic localization or poly(A) tailmediated degradation 33 . Therefore, spliceosomal cleavage may promote nuclear retention and accumulation of functional TER. In S. pombe, disruption of TER spliceosome-mediated 3′-end processing results in reduced accumulation of SpoTER1 and telomere shortening 18 . The identification of additional fungal TERs from more basal species is necessary to determine the evolutionary origins of fungal TER introns.
The identification of the widespread N. crassa-type spliceosomal cleavage mechanism within fungal TERs provides insights into the evolution of TER 3′-end processing. Spliceosomal cleavage was initially viewed as a molecular idiosyncrasy, unique to a single RNA from a single species. Our findings demonstrate that spliceosomal cleavage is the dominant mechanism for Ascomycete fungal TER 3′-end processing, found in the vast majority of species with known TERs. It has yet to be determined whether this mechanism is conserved outside of the Ascomycota phylum, due to the complete lack of TERs identified from the Basidiomycota phylum. Moreover, it is unlikely that this RNA processing mechanism is exclusive to TER and is potentially employed by other non-coding RNAs for 3′-end maturation. This work expands our understanding of the diversity of spliceosomal cleavage mechanisms, the origins and evolution of TER maturation, as well as our knowledge of non-coding RNA processing and biogenesis.

Identification of Schizosaccharomyces TER sequences
Two protein coding genes (sub1, SPAC16A10.03c and rho4, dad1) flanking the SpoTER1 gene served as queries searching against S. japonicus (strain yFS275) and S. octosporus (strain yFS286) genomic sequencing data downloaded from the Broad Institute Database 35 . The intervening sequences between these putative protein coding genes were searched using the program Infernal 36 with a covariance model derived from the sequence alignment of Ascomycete TER three-way junctions 22 and yielded a unique hit for each species. The 5′-and 3′-ends of TERs from S. japonicus and S. octosporus were determined by RACE using the FirstChoice RLM-RACE kit (Life technologies). These cloned TER cDNAs sequences were deposited to GenBank with accession numbers KJ740640 (S. japonicas) and KJ740641 (S. octosporus).

Generation of N. crassa recombinant strains
The expression cassette containing the eGFP gene and the NcrTER 3′-portion that includes exon-1 (residues 1888-2049), the 74-residue intron (residues 2050-2123), exon 2 (residues 2124-2200), and the 500 bp downstream flanking genomic sequence that contains a polyadenylation signal was constructed by overlapping PCR and cloned into the pCCG::NcrTER vector 22,37 . The eGFP gene and the NcrTER 3′-portion were each fused immediately downstream of the ccg-1 5′-UTR and expressed under separate ccg-1 promoters. Intron mutants were generated by site-directed mutagenesis of this parental construct. For N. crassa transformation, the preserved conidia of N. crassa strain FGSC 9720 NC1 (mat A his-3 Δmus-52::bar + ) were inoculated into 25 ml Vogel's minimal solid media in a 500 ml flask. After growing at 30 °C for three days in total darkness, the mycelia were exposed to continuous light for 7 days. Conidia were harvested with 30 ml 1 M sorbitol and filtered through two layers of cheese cloth. Conidia were collected by centrifugation at 3,000 × g for 3 min at 4°C and washed 3 times with 1 M sorbitol. After the final wash the conidia were suspended in 1 M sorbitol to OD 600 of 0.5-1.0. Conidia were transformed with 1 μg of Stu I-linearized wild-type or mutant plasmid DNA by electroporation using a Gene Pulser Xcell system (Bio-Rad) set to 1500 volts and 25 μF. Selection of his + transformants was performed in Vogel's minimum medium supplemented with 2% sorbose in place of sucrose. Following growth at 30°C for three days, the colonies were inoculated onto Vogel's slants. Screening of homokaryon strains was performed by Southern blot analysis.

Southern blot analysis
Genomic DNA was isolated from 20 mg N. crassa mycelia using the Wizard Genomic DNA purification kit (Promega). The mycelia were lysed in 700 μl nuclei lysis solution in a 1.5 ml screw top tube containing 300 mg of 0.5-mm zirconia/silica beads (BioSpec Products) using a MiniBeadBeater-16 (BioSpec Products) at 3,450 rpm for 1 min. The homogenates were incubated at 65°C for 15 min and treated with 1 μl RNase A (50 mg/ml) at 37°C for 30 min. Following RNase treatment, 220 μl protein precipitation buffer (Promega) was added, followed by vortexing vigorously for 20 seconds and incubation on ice for 5 min. After centrifugation at 16,000 × g for 5 min, genomic DNA was precipitated by isopropanol and pelleted by centrifugation at 16,000 × g for 1 min. The DNA pellet was washed once with 75% ethanol, air dried and suspended in 100 μl TE buffer (10 mM Tris-HCl pH 8.0 and 1 mM EDTA pH 8.0). Ten micrograms of genomic DNAs were digested with 10 U BspEI (New England Biolabs) in a 100 μl reaction at 37°C overnight. Digested genomic DNA was phenol/chloroform extracted, ethanol precipitated and resolved on a 0.8% agarose gel. The agarose gel was denatured by immersion in a solution of 1.5 M NaCl and 0.5 M NaOH for 20 min twice, followed by neutralization in a solution of 1.5 M NaCl and 0.5 M Tris-HCl pH 7.5 for 15 min twice. The gel was dried in a gel dryer and pre-hybridized in 20 ml hybridization solution (5X SSC, 5X Denhardt's reagent, 0.1% SDS and 20 mM NaH 2 PO 4 , pH 7.0) at 58°C for 1 hour. In-gel hybridization was performed at 58°C overnight in 15 ml fresh hybridization solution with a 32 P end-labeled oligonucleotide probe targeting the his-3 locus. The gel was then washed at 58 °C for 20 min three times with a solution of 3X SSC and 0.1% SDS and three times with a solution of 2X SSC and 0.1% SDS. The gel was exposed to a phosphorstorage screen and imaged on a FX-Pro phosphorimager (Bio-Rad).

Northern blot analysis
Total RNA was isolated from N. crassa mycelia grown in Vogel's minimum medium at 30 °C with continuous light exposure for expressing ccg1::eGFP and ccg1::NcrTER fusion genes under the light-inducible ccg-1 promoter. Twenty milligrams of mycelia were homogenized in 1 ml TRI-reagent (Molecular Research Center, Inc.) at room temperature using 300 mg 0.5-mm zirconia/silica beads (BioSpec Products) and a MiniBeadBeater-16 (BioSpec Products) at 3,450 rpm for 1 min. Homogenate was supplemented with 0.1 ml 1bromo-3-chloropropane, vigorously mixed, and centrifuged at 4°C for phase separation. Residual phenol in the aqueous phase was removed by chloroform extraction followed by ethanol precipitation. For RNase H cleavage of ccg1::eGFP mRNA to remove the heterogeneous poly(A) tail, 16 μg total RNA was annealed with 600 pmol of DNA oligo (5′-AAGTCGTGCTGCTTCATGTG-3′) by incubation at 65 °C for 5 min and slowly cooled to room temperature. The RNA/DNA hybrid was digested with 5 U of RNase H (New England Biolabs) at 37°C for 30 min followed by ethanol-precipitation. RNase H treated-total RNA was resolved on a 4% polyacrylamide/8 M urea gel, electrotransferred and UV crosslinked to a Hybond-XL membrane (GE Healthcare). For the full length TER, untreated total RNA was resolved on a 1.5% agarose/formaldehyde denaturing gel, capillary transferred overnight and UV crosslinked to a Hybond-XL membrane (GE Healthcare). The membrane was pre-hybridized in ULTRAhyb hybridization buffer (Life technologies) at 65°C for 30 min and hybridized overnight with a 32 P radio-labeled riboprobe targeting the ccg-1 5′-UTR fused to the eGFP mRNA and the NcrTER 3′-fragment. The blot was washed at 65°C twice in a solution of 1X SSC and 0.2% SDS for 10 min, and twice in a solution of 0.2X SSC and 0.1% SDS for 30 min. The blot was exposed to a phosphorstorage screen and imaged on a FX-Pro phosphorimager (Bio-Rad).

RT-PCR
Total RNA was treated with RNase-free DNase I (Zymo Research Corporation) to remove any genomic DNA contamination, followed by phenol/chloroform extraction and ethanol precipitation. For cDNA synthesis, 2 μg DNase I-treated total RNA was annealed with 100 pmol oligo d(T) 25 DNA primer in the presence of 10 pmol each dNTP by incubation at 65 °C for 3 min and quickly cooled on ice. A 20 μl reverse-transcription reaction was performed in 1X Protoscript II Reverse Transcriptase Reaction buffer (50 mM Tris-HCl pH 8.3, 75 mM KCl and 3 mM MgCl 2 ), 10 mM DTT and 200 U ProtoScript II Reverse Transcriptase (New England Biolabs) at 48°C for 30 min. The reaction was terminated by incubation at 80°C for 5 min. PCR was performed in a 25 μl reaction with 1 μl of ½ X diluted reverse-transcription reaction in 1X Q5 Reaction buffer (25 mM TAPS-HCl pH 9.3, 50 mM KCl, 2 mM MgCl 2 and 1 mM β-mercaptoethanol), 0.2 mM each dNTP, 0.5 U of Q5 DNA Polymerase, 0.2 μM forward primer (5′-CACAACACTCAAACCACTTTCAC-3′) and 0.2 μM reverse (5′-TGTTCGACACCGGTCACACT-3′) primer. PCR products were resolved on a 3% agarose gel. The DNA bands corresponding to the spliced products were gel extracted using Wizard SV gel and PCR clean-up system (Promega) and sequenced.

Supplementary Material
Refer to Web version on PubMed Central for supplementary material.   The AU/AG intron is predominate and basal (thick blue line) among fungal TER introns. The GU/AG or AU/AC (U1/U2-type) introns are the putative predecessor of the fungal TER AU/AG intron with only a single point mutation (red) in either intron necessary to convert into this intron type (dashed red lines). From this putative ancestral fungal TER AU/AG intron, the TER GU/AG type intron appears to have arisen twice independently (purple). A subgroup of budding yeasts has lost TER spliceosomal cleavage.