Human Molecular Genetics Advance Access originally published online on March 24, 2006
Human Molecular Genetics 2006 15(9):1497-1512; doi:10.1093/hmg/ddl068
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
A profile of alternative RNA splicing and transcript variation of CACNA1H, a human T-channel gene candidate for idiopathic generalized epilepsies
1Department of Physiology and 2Department of Neuroscience, Johns Hopkins University School of Medicine and 3Department of Medicine, University of Chicago School of Medicine, 725 N. Wolfe Street, Baltimore, MD 21205, USA
* To whom correspondence should be addressed. Tel: +1 4106141505; Fax: +1 4106148331; Email: wsagnew{at}jhmi.edu
Received January 13, 2006; Accepted March 16, 2006
| ABSTRACT |
|---|
|
|
|---|
Highly alternatively spliced genes may provide complex targets for disease mutations. Structural changes created by missense mutations may differentially affect the activity of alternative gene products, whereas missense, silent and non-coding mutations may alter developmental regulation of splice variant expression. CACNA1H is a human gene encoding Cav3.2 low-voltage-activated, T-type calcium channels associated with bursting behavior in neurons and has been linked to more than 30 mutations apparently predisposing to childhood absence epilepsy (CAE) and other idiopathic generalized epilepsies (IGEs). Biophysical properties, including the effects of missense mutations, have been evaluated previously for a single splice form of Cav3.2 expressed in transformed cell lines. We here show that CACNA1H is alternatively spliced at 1214 sites, capable of generating both functional and non-functional transcripts. Variable cytoplasmic and extracellular protein domains point to likely differences in gating behavior, sensitivity to neuromodulation and interactions with extracellular matrix. Biophysical profiles of selected physiological Cav3.2 forms reveal variations in kinetics and steady-state gating parameters, most likely to affect membrane firing. These were comparable to or larger than changes reported for previously studied mutations. Missense CAE and IGE mutations were clustered near segments associated with anomalous splicing. Missense and silent mutations were found to destroy, create or change the regulatory specificity of predicted exonic splicing enhancer sequences that may control splicing regulation. We discuss a paradigm for CACNA1H expression of Cav3.2 subunits, which may influence future basic and clinical studies.
| INTRODUCTION |
|---|
|
|
|---|
Elucidating the relationship of mutations to genetic disease requires explicit knowledge of the paradigm for gene expression. This may be complicated for highly alternatively spliced genes capable of expressing any of many different transcripts in different cells or in the same cell at different stages of development (1
T-type, low-voltage-activated calcium channels appear to play significant and possibly obligatory roles in idiopathic generalized epilepsies (IGEs) (5
,6
). Recently, a series of 33 mutations located in the T-channel subtype, CACNA1H, have been linked to IGEs, including childhood absence epilepsy (CAE) (7
11
). It remains unclear whether these are genuine disease-causing or predisposing mutations or simply rare polymorphisms. Contributing to this uncertainty is the fact that although 16 of the mutations are missense changes that alter channel structure and biophysics (8
,9
,11
), the remainder are synonymous polymorphisms found in the ORF or are mutations in non-coding intron sequences. Current models of Cav3.2 expression do not suggest how these silent mutations may contribute to disease.
Nearly all published reports, and all clinical studies of CACNA1H, have focussed on a single form of Cav3.2 subunit. This is perhaps surprising that insofar, as the close paralog, CACNA1G, has been shown to express a transcriptome of more than 30 variant Cav3.1 transcripts that are expressed in complex, developmentally regulated fashion during human brain development (12
14
). The present studies suggest a similar complexity for CACNA1H expression, which will influence the analysis of genetic errors.
IGEs include CAE, juvenile absence epilepsy, juvenile myoclonic epilepsy and myoclonic astatic epilepsy (15
). A common aspect of the IGEs is the generation of synchronous 36 Hz spike-wave-discharges (SWDs) that spread over both cerebral hemispheres. Recent investigations suggest that seizures are initiated in the neocortex, where CACNA1H is expressed, but rapidly engage the thalamus in patterns of anomalous firing that leads to synchronized discharges and generalized seizure (5
,16
,17
). Many lines of evidence suggest a role for T-type channels. SWDs are sensitive to ethosuximide, an anti-absence drug that reduces T-currents in thalamic relay neurons (18
20
) and also acts on non-inactivating Na+ currents and Ca2+-activated K currents (6
,21
). Transgenic mice lacking one subtype of T-type channel, CACNA1G, are resistant to pharmacological induction of seizures and possess thalamic neurons that lack normal SWD firing patterns (22
). Knockout mice lacking the P/Q type channels exhibit spontaneous SWDs: these do not occur when CACNA1G is simultaneously inactivated genetically (23
). Quantitative increases in T-currents were observed in the coloboma mouse model for absence epilepsy (24
,25
) and in the genetic absence epilepsy rat from Strasbourg (GAERS) (26
). Subsequently, in situ hybridization studies of GAERS brain tissues detected increased transcript abundance for CACNA1G and CACNA1H genes in thalamocortical neurons and neurons of the reticular nucleus (27
).
Three genes encode T-type channels, CACNA1G, CACNA1H and CACNA1I, expressing Cav3.1 (
1G), Cav3.2 (
1H) and Cav3.3 (
1I) subunits, respectively (reviewed in 28
,29
). All three types are expressed in scores of specialized regions of mammalian brain, including neocortex (30
). T-channels activate at negative membrane potentials and inactivate at potentials below than those that activate L-type channels: they generally exhibit slow gating kinetics, small single channel conductance, low permeation selectivity for Ba2+ over Ca2+, block by Ni2+ and characteristic sensitivity to non-dihydropyridine drugs, including mebifradil. Although capable of regenerative behavior, they typically mediate long-lasting depolarizing plateaus (window currents) that trigger bursts of action potentials by fast-gating sodium and potassium channels (6
). T-channels control a wide diversity of electrical events, including slow oscillations of neuronal firing threshold, low threshold spikes and rhythmic burst firing of thalamic relay neurons, pacemaker firing of cardiac myocytes, bursting responses of neurons to EPSPs, paradoxical rebound burst firing in response to IPSPs, synaptic coupling in dendrites, neurosecretion in exocrine cells, and electrical events associated with fertilization in reproductive cells. In addition to being abundant in the CNS (30
), CACNA1G and CACNA1H are found in various peripheral tissues (e.g. liver, pancreas, heart, kidney, adrenal cortical cells and vascular smooth muscle).
Cav3.2 cDNA from CACNA1H has been cloned from human heart (31
), carcinoma cell line (32
) and rat brain (33
) and characterized in respect to tissue distribution (30
): Cav3.2 expressed in heterologous cell lines has been extensively analyzed by patch-clamp recording (34
36
).
Cav3.2 channels are sensitive to second messenger neuromodulation by inhibitory G-protein subunits and calcium calmodulin-dependent protein kinase IIc (CamKIIc) (37
39
). Cav3.2 channels are the targets of both anti-convulsant and anti-hypertensive drugs (36
,40
). In nearly 30 studies, CACNA1H has been represented by a single Cav3.2 splice variant, 545 in our nomenclature. A second form (513), lacking exon 26, has been described in specialized human reproductive cells (41
) and rat brain (33
): no functional consequences of Cav3.2 splicing have previously been described.
We here find that CACNA1H is alternatively spliced at 1214 sites within the ORF, with the potential to generate more than 4000 alternative mRNA sequences. Cloning provides evidence for a family of physiological Cav3.2 transcript variations. The most previously studied variant, 545, comprised less than half of transcripts encountered in fetal human brain. We confirm that eight forms of Cav3.2 subunits expressed canonical T-currents in stably transformed mammalian cell lines. Six forms studied systematically show significant differences in gating parameters that are likely sufficient to affect membrane firing.
Missense mutations that affect gating may act differently in splice forms that, themselves, exhibit intrinsically variable biophysical properties. In the studies of Vitko et al., (11
), 11 of 12 CAE missense mutations were shown to change kinetic and steady-state gating parameters of Cav3.2, sufficiently to alter firing patterns in a computational model of thalamic relay neurons. Gating changes were similarly induced by a common polymorphism (R788C) that occurs in non-afflicted controls. A variant that combined R788C with the CAE polymorphism G773D displayed gating behavior different from either SNP alone, highlighting the influence of the peptide background on the functional manifestation of a mutation.
Both missense and silent mutations have the potential to alter splicing regulation (42
). We note that many CAE/IGE mutations cluster in or near two domains whose deletions predict truncated reading frames. We discuss three examples of missense or silent mutations that induced little or no change in channel gating, but destroy, create or change the predicted regulatory specificity of strong candidate exonic splicing enhancer (ESE) sequences that could regulate splicing biases at anomalous splice site boundaries.
| RESULTS |
|---|
|
|
|---|
CACNA1H transcripts are alternatively spliced at multiple sites in the ORF, each possessing viable donor and acceptor splice sequences (43
CACNA1H organization
CACNA1H encompasses >70 kb of human chromosome 16p13.3 (cf. GenBank entries AL031712
[GenBank]
, AL031703
[GenBank]
and AL 031715, compared with cloned cDNA channel sequences AF073931
[GenBank]
and AF070604
[GenBank]
). The ORF includes 36 exons, illustrated in Figure 1A, beginning at nucleotide 249 in exon 2 and ending at nucleotide 7310 in exon 36. In this report, exons are in bold face plain text, whereas splice variants are bold face italics. Variants are named according to a binary nomenclature that encodes variable domain content, described in the legend to Table 1.
|
|
Sites of variable splicing
Variable sites were identified by short amplicon exon scanning (48
The 12 alternative splice sites in the Cav3.2 ORF are shown in Figure 1A and B (see Supplementary Material, Table S4 for sequences). Splice variations included commonly recognized patterns of insertion/deletion (44
) (e.g. deletion of cassette exons, mutually exclusive cassette exon splicing, use of alternative donor or alternative acceptor sites within exon sequences and retained introns). Two novel splicing patterns were seen, in which new sequences were introduced by upstream frame-shift deletions coupled to concerted downstream intron-segment insertions that returned the normal reading frame.
Variable exons 14, 17 and 26 are canonical cassette exons. Use of alternative donor or acceptor sequences within exons leads to excision of partial exon segments 9B, 25A, 25C and 35A and two partial deletions in exon 31 (
31A and
31A'); of these, segment 25C was not known to comprise coding sequence for CACNA1H and, like 25C in Cav3.1, constitutes a retained intron. Other nominal intron segments, i23, i27B and i32B, may also be alternatively retained. Retention of segments i3A and i3C, suggested in published sources (TIGR human EST database NP 336223), was not confirmed in fetal human brain but is provisionally represented in Figure 1A. Of the variations detected, only the deletion of exon 26 had been previously observed in human reproductive cells (41
): in rat brain,
26 was the only configuration recovered by McRory et al. (33
). As illustrated, after the convention of Nidal Ginard et al. (49
) (Fig. 1C), re-assortment of 24 alternative configurations at 12 sites generates 4096 combinatorial transcripts. Future confirmation of i3a and i3c splicing could expand this 4-fold.
An explosion diagram of all possible Cav3.2 protein coding sequences is depicted in Figure 1D. Many variable sites lie in cytoplasmic regions, including the N- and C-terminus and interdomain loops, ID12, ID23 and ID34. One extracellular variation was contiguous with the pedestal that bears multiple consensus sequences for N-linked glycosylation, adjacent to the pore-forming P-loop of D3. When retained, intron i23 extends the extracellular pedestal from 30 to 92 amino acids. A BLAST search of the i23 predicted amino acid sequence reveals
40% amino acid identity with Type I and Type VII collagens.
A striking parallel in the alternative splicing pattern of segments 25C and 26 in ID34 of CACNA1H and paralog CACNA1G is shown in Figure 2. Although flanking segments 25AB and 27 are well conserved between these genes, 25C and 26 are of different length and bear no detectable homology at the nucleotide or amino acid. In both genes, these segments are spliced in either mutually exclusive fashion or are both deleted (
forms). Insertion of 25C or 26 in Cav3.2 transcripts creates strong consensus sites for protein kinase A phosphorylation at the splice junction by codon splitting. A strong protein kinase A site is also created for segment 25C, but not 26, in Cav3.1 transcripts by the same mechanism. Although we did not detect the double insertion 25C26 in our preliminary survey of Cav3.2 transcripts, such forms are found for Cav3.1 in fetal human brain and glioma cell lines (12
,50
) and in mouse EST databases (M. Emerick and W. Agnew, unpublished data).
|
A cytoplasmic domain encoded by exon 17 encompasses 127 amino acids in ID23: this segment possesses a quartet of consensus phosphorylation sites highly conserved in all three T-channel paralogs. This segment contributes half of the ID23 mass (14 kDa) and has been shown by Barrett and coworkers (39
2 and ß2
4. This segment also contains serine 1198, the site of phosphorylation by CamKIIc, which strongly modulates gating (37
17 variants may lack neuromodulation by both second messenger pathways. C-terminal segment 35A, encoding 11 amino acids, can be deleted through the use of an alternative internal acceptor sequence. By comparison, Cav3.1 C-termini exhibit alternative splicing of segments 34, 35 and 38B: 38B exhibits strong interactions with multiple other domains to affect gating. This raises the possibility that 35A might similarly modulate the behavior of Cav3.2.
Nominally deleterious variants, including
25A,
31A and
31A', are found in both CACNA1H and CACNA1G. Consensus splicing sequences for exon 31A, including a long, canonical polypyrimidine element, are well conserved. In Cav3.1, 31A deletion is predicted to create a downstream stop codon and PCT. This is likely to elicit nonsense-mediated decay (NMD), consistent with the low observed frequencies of
31A forms. These variants, nevertheless, showed strong evidence of developmental regulation, suggesting that they may provide a means of discontinuing Cav3.1 expression. Here, in Cav3.2,
31A forms were identified (discussed subsequently), in which the predicted downstream PCT was averted by concerted insertion of a small intron segment i32B, resulting in the insertion of a new coding segment in the region D4S5 (Fig. 1D, orange segments).
Splicing of 25A is conserved in both genes: although the specific acceptor splice sites differ slightly in sequence, they lie in the same location and flanking sequences are identical. In Cav3.1, 25A splicing is in frame and is strongly developmentally regulated. In Cav3.2, 25A leads to a frame shift. As for
31A,
25A forms were identified, in which downstream PCT was averted by insertion of a partial intron segment, i27B, which returns the normal reading frame, inserting a new sequence in the region D3S6ID34. These forms were detected by short amplicon PCR of library cDNAs; one example of both variations occurred in the survey of longer cDNAs from an independent cDNA library.
9B variations were of the order of >10% of PCR products from each of the tissues examined (cf. Fig. 3). Splicing at this site involves conventional U2 boundaries. Segment 9B encompasses 26 amino acids of cytoplasmic ID12, beginning at amino acid 462 and ending at 489: this deletion changes the reading frame and is predicted to elicit downstream PCT in exon 10 at amino acid 695. Interestingly, a cryptic in-frame translation initiation site lies at nucleotide 2571 (amino acid 775), 57 nt (19 amino acids) before S1 of D2, resulting in a second ORF encompassing the last homology domains and C-terminus. This splicing pattern, which does not create a chain termination codon immediately 5' to a normal splice site, may escape the NMD pathways for degrading invalid mRNA (51
,52
), perhaps contributing to the relative abundance of
9B transcripts.
|
Dolphin and coworkers (53
9B splicing, for several genes including CACNA1H. This suggests that activation of 9B splicing could terminate Cav3.2 expression and prevent synthesis from normal transcripts in the expressing cell. Mutations influencing
9B formation could be deleterious (cf. Discussion). There is some possibility that both peptide segments could be synthesized and assemble to form functional bi-molecular
1 subunits (54
Multiple tissues exhibit Cav3.2 splice variations
Cav3.2 alternative splicing was not limited to fetal brain. Figure 3 illustrates short amplicon PCR of representative sites in human fetal and adult brain, heart, kidney and liver cDNA libraries. Splicing variants at all identified sites were evident, varying in relative abundance with tissue source, especially noticeable for 
, 25C and 26 configurations of ID34, 35A and
9B.
17 forms were reproducibly detected at low levels in both adult and fetal brain tissues (data not shown).
Evidence for naturally occurring splice variants
To search directly for physiological Cav3.2 variants, we scored variable sites in fifty-three 3.3 kb Cav3.2 cDNAs encompassing exons 1936, inclusive.
Figure 4 illustrates a gallery of 12 splicing patterns detected among 53 Cav3.2 cDNAs. Relative abundance and splice site compositions are given in Table 1. Of these clones, 30% possessed the 3'-configurations of variant 545: as these clones lack 5'-coding exons 218, containing five additional sites of variation, this may overestimate 545 abundance. Of the remaining 70% of transcripts, variants 513, 512 and 521 constituted 24.5, 13.2 and 9.4% of the population, respectively; 517, 522, 544 and 577 were each 3.8%, with the remaining forms (640, 672, 688 and 800) represented by a single clone (1.9%). Full-length cDNAs in Figure 4B were constructed from the isolated variants with 545 3'-configurations. Of the 12 variants detected, six were functional when expressed in HEK293/FLP cells with the 545 configuration of exons 218 (544, 545, 512, 513, 577 and 800), as described subsequently: synthetic variant 576 was functional, as were constructs lacking exon 17. PCR and preliminary cloning data suggested that variants with 25C replacing 26 are more abundant in adult tissues, a pattern also seen for Cav3.1 (12
).
|
This inventory of Cav3.2 variants does not include variations in coding exons 218. Full-length cloning of this gene is hindered by the presence of GC-rich regions in the 5'-ORF, which will necessitate stringent methods to overcome: in addition, transcript abundance in other tissues, including adult brain, is one to two orders of magnitude lower than fetal tissues. Our previous survey of nearly 1600 CACNA1G transcripts suggests that marked differences in forms and relative abundance of Cav3.1 transcripts could occur in different tissues, such as adult and fetal brain. A similar large-scale survey of full-length cDNAs may thus be justified to determine the complete inventory and the regional and developmental specificity of Cav3.2 variant expression.
Alternative splicing alters Cav3.2 biophysics
As many as 32 in-frame variants of Cav3.2 could arise from splicing of segments 17, i23, 25C, 26 and 35A. No biophysical studies have been reported on the effects of splicing at any of these sites. To examine a systematic subset of Cav3.2 variants involving segments 25C, 26 and 35A, analogous to parallel studies of Cav3.1, we prepared six full-length cDNAs for expression in mammalian cell lines. These were selected because of their structural parallels with variations of Cav3.1 that are prominently regulated during development and that induce diverse changes in Cav3.1 gating (12
).
These were 545 (canonical Cav3.2), 544, 512, 513 and 577. Variant 576, not yet found as a natural form, was also prepared. These were expressed in stable lines of HEK/FLP cells (Materials and Methods); all of these clones expressed currents well in this cell background. Data for configurations 25C, 26 or 
are often presented in pairs, +/35A.
Figure 5A and B illustrates that these variable domains are not obligatory for channel activity and show the properties of variant 512, lacking 25C, 26 and 35A. Segments 17 and i23 and possible inserts i3A and i3B are also not required (data not shown): that is, the core subunit, lacking these cytoplasmic and extracellular domains, exhibits canonical T-channel behavior. Activation was evaluated by simple conductance transforms of peak currentvoltage relationships, as described in Materials and Methods. All six variants activated similarly, i.e. the half-points of peak conductance only varied between 39 and 44 mV (Table 2). Activation was not affected by segment 35A. When data were combined on the basis of the exon composition in the ID34 linker (25C/26/
), there was a small although significant, dependence on the linker (P=0.01, Students t-test). This is evident in Figure 5D.
|
|
The functional window supported by each variant was evaluated by measuring functional availability after 200 ms at various conditioning potentials. Typical data for 512 are shown in Figure 5C, and summary data for all of the variants are shown in Figure 5D and E. For Cav3.2, activation typically begins near 70 mV, with conductance midpoints between 45 and 48 mV and peak currents at 30 mV. Differences in activation were modest, varying over a 5 mV range, compared with an 11 mV range for Cav3.1. The midpoints of availability (200 ms) varied over a 10 mV range and were strongly, linearly correlated with the midpoint for activation (Fig. 5D). As found for Cav3.1, variants that gated at the more negative potentials had smaller functional windows than those gating at positive voltages (Fig. 5E). Functional window magnitude (200 ms) varied from a minimum of 22 to >30% and the voltage for peak window currents varied over a 9 mV range. Variants possessing exon 26 fell at the most negative potentials,

forms at the most positive potentials and 25C forms in between; this is opposite to the finding with the paralog gene CACNA1G, in which windows for forms expressing 25C were displaced
11 mV from forms lacking inserts, with 26 forms falling in between. This suggests that a broad shift with brain maturation toward Cav3.1 forms gating at more negative potentials may not occur for Cav3.2 variants.
In contrast to the small, although orderly, differences in midpoints and windows, there were larger differences in the time courses of the variants. Similar to all T-type channel data previously reported, all variants displayed little voltage dependence in the decay of current at positive potentials. Figure 6A illustrates a representative voltage-clamp time course for wild-type variant 545. Currents are relatively slow to develop: at negative potentials, peak currents occurred between 50 and 100 ms, whereas at positive potentials, current peaked at
6 ms. Figure 6B and C illustrates voltage dependence of
m and
h for all six variants; the data are well described by HodgkinHuxley formalism. The voltage ranges depicted lie between the approximate midpoint for activation (45 mV) and a potential that elicits maximum conductance (10 mV).
h approached a constant value at more positive potentials, as expected for a voltage-independent inactivation process. There was no significant difference in time constants between variants that contained and was missing exon segment 35A; therefore, for the determination of the potential roles of segments 25C and 26, data +/35A were pooled. Exon 26 produced the fastest currents (544 and 545), both for
h and
m. Conversely, in 
forms, both
m and
h were slowest (512 and 513). Interestingly, however, when 25C was present, the activation time constant
m was similar to when 26 was present, but current decay,
h, was slow, similar to when neither 25C or 26 was present.
|
None of the variants inactivated appreciably from closed states over the first 200 ms, i.e. there was virtually no reduction in peak current magnitude when cells were depolarized to voltages negative to those at which currents were detectable. As a result, the midpoint of voltage-dependent availability after 200 ms was only modestly more negative than the conductance relationship, by between 9 and 12 mV, depending on the variant. Figure 7A shows this for one of the variants (544: +26C/35A). Thus, there was a broad voltage range over which window currents were available after 200 ms (almost 30 mV for most variants) (see Table 2 for data summary).
|
Cav3.2 variants show marked difference in availability at 200 and 5000 ms. As illustrated in Figure 7A, the variant 544 (+26/35A) comes to activation equilibrium within 200 ms: inactivation, however, is incomplete in this time frame. At 5000 ms, the availability midpoint has shifted from 55 to 76 mV, with an increase in slope factor. This most likely results from slow steps of activation at negative potentials that, absent closed state inactivation, delay gating equilibrium. A small true window current persists after 5000 ms, (<12% of Imax) centered
60 mV.
These properties vary dramatically with splicing, as illustrated in Figure 7B. Variant 512, lacking 25C, 26 and 35A, exhibits a much smaller shift in the midpoint voltage for availability after 5000 ms (from 46 to 58 mV); thus, the availability midpoint for 512 at long times is 18 mV more positive than that for 544 (+26/35A). A pronounced anomaly in the availability curve, a shoulder of
20% of the maximum current, is present, with a midpoint voltage of 100 mV. This profile results in a much larger (610-fold) persistent window current (
12% of Imax at 5000 ms), with midpoint voltage 50 mV, 10 mV more positive than for 544. This rise could also be detected at 200 ms for availability measured between 80 and 110 mV: currents increased
10% over this range, small enough not to be detected easily in Boltzman fits to the 200 ms availability data. It is likely that the window was slightly underestimated as a result. Similar behavior was observed in five cells expressing 512: for other variants lacking 26, it was much less common (Fig. 7).
Figure 5D shows that the 200 ms window was smallest in magnitude, and most negative, when 26 was present, and largest and most positive in 
forms. Given that there were no statistical differences in conductance between the variants, it must be the case that recovery from inactivation was most rapid in 
forms and slowest when exons 26 and 35A were present, with other variants likely in between. This trend was generally observed. Figure 7C illustrates the time course of recovery from inactivation by the six tested variants. When segment 26 is present, recovery satisfies a mono-exponential time course with time constant of
1 s; when 25C is present, a similar mono-exponential recovery is well fit with a time constant from 0.9 to 1 s. When both segments are deleted (
), we found a nominally 3-fold overall faster recovery, following a bi-exponential time course with one time constant (3040% of current) of 0.1 s and a second (6070%) varying between 0.6 and 0.7 s. Variations of 35A had minimal impact on these profiles. Interestingly, those variants that recovered from inactivation most rapidly (
forms) also had the longest time constants for inactivation from the open state, taken as the time constant of decay at positive potentials. This suggested that 25C and 26 stabilized the inactivated configurations.
| DISCUSSION |
|---|
|
|
|---|
The present studies indicate that CACNA1H is extensively alternatively spliced and generates a family of variant transcripts. The emerging paradigm for CACNA1H expression can provide an informative context for considering disease mutations (7
The number and pattern of CACNA1H alternative splice sites are similar to that previously reported for CACNA1G. A schematic for regulated splice variant expression, adapted from that proposed for CACNA1G (12
), may be considered. In this model, splicing at each site in the precursor mRNA is modulated by trans-acting splicing regulatory factors, typically serine/arginine (SR)-rich proteins that bind to cis-ESE sequences, which typically reside near or overlap with 5'- and 3'-boundaries of each alternative segment. Several SR factors may act competitively or synergistically to influence the bias toward retaining or skipping the alternative segment (3
). To select a particular transcript, the ensemble of variable sites must be regulated in concert. In this model, transcription factors may promote the activation of gene transcription and may simultaneously activate a cascade of SR factors that select the variant for a particular cell lineage. Members of the various classes of SR factors may themselves be subject to regulated alternative RNA splicing (3
).
In this model, missense mutations may alter the function of expressed channels: such effects may be permuted among the various alternative subunits and differentially disturb disparate electrophysiological phenomena (e.g. synaptic coupling, firing threshold oscillations, SWDs and rebound burst firing). Missense, silent and even non-coding mutations that produce little or no direct changes in channel function may nevertheless perturb ESE regulatory sites and thus alter normal patterns of splice variant expression. Both of these phenomena are suggested by the distribution and properties of mutations that have been implicated in CAE and IGE.
First, we here find that Cav3.2 expression is diversified; only one of these variations has previously been detected. The principal form previously studied, variant 545, comprised a minority of the forms detected among 53 cDNAs from fetal human brain. Moreover, splicing of 25C, 26 and 35A shifted the voltage-domain for gating, the kinetics of activation, inactivation and recovery from inactivation and the magnitude and voltage midpoints for functional window currents (e.g. persistent currents in the 200 ms time frame). The magnitude and voltage-dependence of steady-state window currents (>5 s) were strongly affected: when accompanied by deletion of 35A in the C-terminus, deletion of both 25C and 26 resulted in a 612-fold increase in peak steady-state window current and a 12 mV shift in the voltage midpoint. This was accompanied by a relatively enormous +18 mV change in the availability midpoint. Among the six forms described here, we observed a continuum of inactivation and availability midpoints (cf. Fig. 5D). Splicing variations that involve segments 17, i23, 25C, 26 and 35A, as well as possibly i3A, i3B and 9B (discussed subsequently), could generate 32 or more variants, which may differ in regard to gating, localization or sensitivity to neuromodulation. Thus, both differential modification of function by missense mutations and altered patterns of expression of otherwise physiological variants by mutations of regulatory sequences could underlie pathological manifestations of CAE or IGEs more generally.
Although ion channel mutations unequivocally account for some monogenic diseases, the disease significance of mutations identified in genetic surveys of CAE and IGE patient cohorts remains to be fully established. The distributions and properties of some of the CAE and IGE associated mutations may be considered in relation to the emerging picture of CACNA1H alternative splicing. Thirty-three such mutations have been identified: Chen et al. (7
) reported 29 mutations, including 12 missense polymorphisms, among a cohort of 118 Chinese CAE patients; none of these was found in a normal population of 230. Eleven of these were subsequently found to alter gating parameters of human variant 545 sufficiently to change firing patterns in a computational model of thalamic relay neurons (11
). In addition, Heron et al. (10
) identified four missense mutations in a cohort of 192 patients displaying a variety of IGEs; none of the four mutations occurred in a set of 92 healthy controls. The association with disease, however, was ambiguous: no variant segregated with a single syndrome; in two instances, unaffected siblings were encountered which harbored the mutation; in two cases, siblings without the mutations nonetheless exhibited the disorder. The loose association with disease could in part reflect the polygenic nature of IGEs (6
,56
). Functional analyses of three of the four Heron mutations identified two, in which small, but statistically significant, gain of function changes were identified, also in the background of human variant 545 (9
). The fourth mutation produced a frame shift in ID12, generating a truncated, non-functional product.
Figure 9 examines the distributions of several of the 33 reported mutations. Interestingly, 11 of 16 misssense mutations lie in or flank segment 9B. Specifically, all four of the IGE associated mutations of Heron et al. (10
) as well as two CAE missense mutations (and one CAE silent mutation) lie within 9B. One CAE missense mutation (C456S) lies in 9A near the splicing junction, and four are found in exon 10 near the 9B boundary. This distribution suggests that this region may be a particular hotspot for disease-linked variations. Thirteen of these misssense mutations produce modest to significant changes in channel gating (9
,11
). Liu et al. (42
) have demonstrated that nonsense, missense and silent mutations can give rise to anomalous patterns of exon skipping, in addition to direct effects of altered peptide sequence on protein function.
|
We have noted earlier that 9B deletions are predicted to generate a truncated peptide comprising the N-terminus, D1 and a portion of cytoplasmic ID12. In well-studied systems, such as the Drosophila sex determination genes sex lethal and transformer, regulated splicing events that introduce stop codons and truncated proteins serve to turn off gene expression (57
9B transcripts could lead to significant increases in Cav3.2-mediated T-currents (if 9B deletion were diminished) or to decreases (if 9B deletion were enhanced).
As discussed by Mount (43
), it is not possible to reliably identify physiological splice sites by any strategy other than direct examination of cDNAs. Splicing can involve any number of nucleotides, including 0 (e.g. resplicing) and both conventional and non-canonical spliced donor or acceptor sites. He concludes that physiological patterns of alternative splicing can only be determined by direct examination of cDNAs. Rapid advances have been made, however, in predictive algorithms for locating candidate ESE sequences (46
,47
) specific for different classes of SR factors. In at least three instances, CACNA1H mutations were noted in previous studies to elicit minimal or no changes in Cav3.2 variant 545 function (C456S, D1463N and A1765A). We examined these mutations in relationship to candidate ESE sequences using the web-based programs RESCUE-ESE (46
) and ESE-Finder (47
).
Figure 9B indicates that six candidate ESE hexamer consensus sequences are located in the middle of exon 9, within a few base pairs of the 9B splice site. As noted by Fairbrother et al. (46
), strong ESE activity is often associated with clusters of three to six sites that neighbor or overlap the consensus splice site. Using the alternative criteria of Cartegni et al. (47
), we identified three strong SC35 specific sites, displaced slightly 5' from the clusters of six sites. This program distinguishes between four classes of SR factors. We note that for mutation C456S, a single base change (1455T
A) eliminates one of the overlapping SC35 sites and converts the specificity of the remaining site to SRp55. These changes in ESE specificity are perhaps as likely to be pathogenic as the minimal changes noted for the channel biophysics. According to Vitko et al. (11
) and Khrosvani et al. (9
), this mutation produced minimal changes in Cav3.2 channel gating in the human and rat forms of the 545 variant. If 9B splicing is subject to SR factor regulation, perhaps to silence gene expression, mutation C456S, which lies in 9A, just 5' to the 9B splice site, could cause 9B skipping or, alternatively, interfere with normal regulation of 9B splicing.
In Figure 9C, we make a similar analysis of the silent mutation A1765A, associated with the 31A/A' splice junctions: this mutation does not change the ORF and thus cannot elicit a direct change in the biophysics of Cav3.2. However, this segment is associated with two strong candidate ESE sequences for SR factors SF2/ASF and SRp55. The G
A change results in a slight increase in the predicted quantitative strength of the SF2/ASF splice site and introduces a third very strong overlapping site specific for SRp40. These changes could alter the sensitivity of 31A and/or 31A' splicing biases to expressed patterns of splicing regulatory factors. Insofar as these deletions appear likely to prevent normal channel expression, changes in regulation could result in either loss or gain of function by activating or silencing 31A/31A' splicing.
In Figure 9D, we illustrate candidate splice sites near the i23/23 splice junction. By the criteria of Fairbrother et al., we identify a cluster of three overlapping hexamer candidate ESE sequences, together with two additional flanking sequences. These are consistent with similar findings by the criteria of Carneci et al., which identify three clusters of overlapping sites of 14 candidate sequences, selective for four SR factors, SF2/ASF, SC35, SRp55 and SRp40. The most complex of these clusters involves eight ESEs of the four types. Such clustering may potentiate competitive or synergistic interactions in regulating splicing at this boundary. The mutation D1463N was reported by Vitko et al. (11
) to cause no detectable changes in Cav3.2 gating parameters in the 545 variant background. As illustrated, the G
A base conversion has several effects on the third ESE cluster. Three ESE sequences are eliminated altogether, and the strength of two remaining SR35 sites are altered. We have encountered cDNAs in which intron i23 is retained, creating an additional extracellular loop with significant homologies to collagen. These ESE sequence alterations at the i23/23 boundary could alter retention of i23 or the constitutive exon 23 itself, depending on SR factors active in the expressing cell.
Heterologous expression of functional splice variants provides a useful screen for altered channel function, as used here. However, transfected cDNAs are not typically subject to RNA splicing. Development of model systems involving mini-gene expression in cell lines with defined SR factor profiles will be required to reveal whether mutations noted here perturb splicing and may thereby provide a mechanism for disease potentiation (42
).
Collectively, the observations described here indicate that CACNA1H is an extensively alternatively spliced gene, with 12 or more sites of transcript variation within the ORF. A preliminary survey indicates that transcript variation is extensive. Eight structural variants, including six found naturally in a survey of fetal brain cDNAs, gave rise to canonical T-type currents in transformed mammalian cell lines. Variations in splicing induce changes in parameters likely to be crucial to T-channel function. Spliced domains include conservative variations that preserve membrane topology and channel function and anomalous variations that likely generate truncated products. Such truncated forms have been encountered in every calcium channel
1 subunit gene so far examined and may play a physiological role in regulating gene expression. Interestingly, one of these anomalous variable domains, 9B, encompasses or lies adjacent to 11 of 16 missense mutations, as well as one silent and one frame-shift polymorphism, reported to predispose for CAE or IGE. The model proposed for CACNA1H expression suggests that function-perturbing missense mutations could generate permutations of changes among variants likely to be expressed in different cells. Both missense and silent polymorphisms could alter the regulation of splicing, perturbing the normal timing and sites of Cav3.2 splice variant expression.
Taken together, these studies begin to outline a more complex context for the examination of disease mutations for CACNA1H than has been previously considered. Variant proteins could provide differentiated targets for therapeutic agents. Further studies may thus be warranted to identify the full physiological inventory of CACNA1H transcripts and determine the regional specificity of their expression with brain development.
| MATERIALS AND METHODS |
|---|
|
|
|---|
Exon scanning
Primer pairs were designed to generate fragments spanning a single exon or a set of tandem exons. Constitutive exons yield a single PCR product, whereas alternatively spliced segments typically produce multiple bands that may be resolved electrophoretically. Products were isolated, subcloned and DNA sequenced. Variants were compared with genomic data to confirm consensus splice donor and acceptor sequence (44
Cav3.2 library construction
ORF sequences containing exon 19 to the 3'-UTR of CACNA1H message were amplified from human Marathon-Ready fetal brain cDNA library (Clonetech) with Elongase amplification system (Life Technologies). Primer pairs were designed: anti-sense primer: 5'-GCCCTACAAGCCCCAGTGGT-3'; sense primer: 5'-CAGGGTCCATGCAGGGTTCT-3'. Wax beads (AmpliWax PCR Gem 100, PerkinElmer Applied Biosystems) were employed, following the manufacturer instructions, to perform a hot start to increase the specificity of the reaction. The reaction volume was 50 µl, containing 2 mM magnesium (10x Elongase buffer), 200 µM each dNTP (Life Technologies), 200 nM each primer, 2 µl cDNA and 1.4 µl Elongase enzyme mix (1 U/ml). The thermal program was: 94.5°C for 45 s, 94.5°C for 30 s/65°C for 30 s/68°C for 4 min for 24 cycles and a final extension at 68°C for 6 min. Ten microliter of each PCR product was run on an agarose gel adjacent to a DNA ladder to assess specificity and to verify that the product concentration did not exceed a 100 pM in order to keep template switching below a desired level (59
). The
3.4 kb PCR product was excised after electrophoresis on 0.8% agarose, purified from the gel slice on a Qiaquick column (Qiagen) and cloned into pCR-XL-TOPO with the TOPO-XL PCR cloning kit (Invitrogen). White colonies were picked into 200 µl LB with 50 µg/ml kanamycin and grown overnight at 37°C in 96-well tissue culture dishes. Aliquots were lysed for 20 min at 100°C in 1 mM TrisHCl and 0.1 mM Na EDTA, pH 8.0(1/10xTE). The lysates, for PCR screening reactions, were store frozen at 20°C. The remaining bacterial cultures were stored at 80°C in LB with 20% glycerol.
Splicing variant screening
Alternative exon combinations for full-length clones were detected by short-amplicon PCR screening with primer pairs surrounding sites of alternative splicing. The reaction volume was 25 µl, containing 5 µl of bacterial lysate as template and other reaction components as for exon scanning earlier (for primer pairs, see Supplementary Material, Table S2). Splice variants were identified by electrophoresis of PCR products on 4% argarose gel after band validation by sequencing. Clones with
31A and
31A' required direct sequencing.
Splicing in different tissues
Tissue comparisons of splice site configurations were made by amplification of Cav3.2 from Clontech Marathon-Ready human fetal brain, adult brain, heart, kidney and liver cDNA libraries (nos 639302, 639300, 639304, 639305 and 639303) using the primer pairs as Supplementary Material, Table S1. Amplification conditions were the same as earlier. PCR products were resolved on 1 or 4% argarose gel.
Electrophysiology
Cav3.2 splice variant cDNAs were subcloned into pcDNA5/pFRT/Zeocin, and stable cell lines in HEK293/FLP cells were created (Invitrogen, Carlsbad, CA, USA) for each of the variants. Cell lines were maintained with Dulbecco's eagle medium (GibcoInvitrogen) in 100 mM Corning culture dishes with 10% fetal bovine serum, 1% penicillinstreptomycin and 100 µg/ml Zeocin for selection. Media were changed every 48 h and cells were passaged weekly. Whole-cell voltage clamp was performed on trypsinized cells (0.25% TrypsinEDTA, Invitrogen) 36 days after plating. Whole-cell voltage clamp recordings were made using glass pipettes made from alumina silicate capillary tube (Fisherbrand) pulled to resistances of 1.02.0 M with a Flaming/Brown micropipette puller (Model 97, Sutter Instruments, Novato, CA, USA). The pipette solution (24
) contained (in mM) 115 CsCl, 1 tetrabutylammonium chloride, 1 MgCl, 5 Mg-ATP, 1 Li-GTP, 20 HEPES, pH.7.2 (adjusted with CsOH), and 11 BAPTA, with 0.9 mM CaCl2 to set the free Ca2+ at 27 nM. HEK293 cells were superfused with a solution containing the following (in mM): 127 TEA-Cl, 2 CaCl2, 0.5 MgCl2, 10 HEPES, 5 D-glucose and 32 sucrose, pH 7.4 (adjusted with CsOH). Recordings were made at room temperature (2224°C) using pCLAMP8 and an Axopatch 1D amplifier with a DigiData 1322A Series Interface digital-to-analog converter. Cells were held at 110 mV, and for currentvoltage relationship, determinations stepped to various potentials for 200 ms once each 5 s. Steady-state inactivation protocols used 200 ms prepulse steps and a 0 mV test step, with 8 s between trials at a holding potential of 120 or 110 mV. Data were low-pass filtered at 5 kHz and sampled at 20 kHz during recording and digitally filtered at 2 kHz during analysis. Specifics of other protocols are given in Results.
Data were analyzed using locally written programs in Matlab (Natick, MA, USA) and Origin 7 (Northampton, MA, USA). Data were capacity corrected with an average of 16 sub-threshold 20 MV steps and leak corrected using linear prediction from leak currents between 100 and 80 mV or between the current at the holding potential and zero. Cells with greater than 300 pA leak at 100 mV, activation slope factors steeper than 4 mV, or with Gmax<5 nS were excluded from analysis. Conductance and availability curves were constructed from normalized data, which were then fit using a Boltzmann equation, i.e. fraction=1/(1+exp(V1/2Vm)/dx), in order to estimate V1/2, the midpoint of the relationship, and dx, the slope factor, i.e. the number of millivolts over which the relationship changed e-fold. HodgkinHuxley formalism was used to estimate the time constants of activation (
m) and inactivation (
h) (I=m3h) from capacity corrected data, digitally filtered at 2 kHz excluding the first millisecond after the step to avoid filtering artifiacts. Statistical tests for differences used t-tests carried out in Origin or with, when appropriate, Boneferroni correction to minimize false positives from repeated measures in SAS (Cary, NC, USA). All data are reported as mean±SEM.
| SUPPLEMENTARY MATERIAL |
|---|
|
|
|---|
Supplementary Material is available at HMG Online.
| ACKNOWLEDGEMENTS |
|---|
This work was supported by NIH award HL52307 and a Holden Targeted Research Award from the Epilepsy Foundation to W.S.A. and NIH awards R01 HL065680 and RO1 HL62603 to D.A.H. We wish to thank Constance Mlecko and Timothy Millar for excellent technical assistance, Aaron Yamada-Hanff and Jason Yamada-Hanff for assistance with figures and Mark Emerick for critical discussions.
Conflict of Interest statement. The authors declare that they have no conflicts of interest or competing financial interests.
|
| REFERENCES |
|---|
|
|
|---|
- Black, D.L. and Grabowski, P.J. (2003) Alternative pre-mRNA splicing and neuronal function. Prog. Mol. Subcell. Biol., 31, 187216.[Medline]
- Grabowski, P.J. and Black, D.L. (2001) Alternative RNA splicing in the nervous system. Prog. Neurobiol., 65, 289308.[CrossRef][Web of Science][Medline]
-
Faustino, N.A. and Cooper, T.A. (2003) Pre-mRNA splicing and human disease. Genes Dev., 17, 419437.
[Free Full Text] -
Modrek, B., Resch, A., Grasso, C. and Lee, C. (2001) Genome-wide detection of alternative splicing in expressed sequences of human genes. Nucleic Acids Res., 29, 28502859.
[Abstract/Free Full Text] - Crunelli, V. and Leresche, N. (2002) Childhood absence epilepsy: genes, channels, neurons and networks. Nat. Rev. Neurosci., 3, 371382.[CrossRef][Web of Science][Medline]
- Crunelli, V. and Leresche, N. (2002) Block of thalamic T-type Ca(2+) channels by ethosuximide is not the whole story. Epilepsy Curr., 2, 5356.[CrossRef][Medline]
- Chen, Y., Lu, J., Pan, H., Zhang, Y., Wu, H., Xu, K., Liu, X., Jiang, Y., Bao, X., Yao, Z. et al. (2003) Association between genetic variation of CACNA1H and childhood epilepsy. Ann. Neurol., 54, 239243.[CrossRef][Web of Science][Medline]
-
Khosravani, H., Altier, C., Simms, B., Hamming, K.S., Snutch, T.P., Mezeyova, J., McRory, J.E. and Zamponi, G.W. (2004) Gating effects of mutations in the Cav3.2 T-type calcium channel associated with childhood absence epilepsy. J. Biol. Chem., 279, 96819684.
[Abstract/Free Full Text] - Khosravani, H., Bladen, C., Parker, D.B., Snutch, T.P., McRory, J.E. and Zamponi, G.W. (2005) Effects of Cav3.2 channel mutations linked to idiopathic generalized epilepsy. Ann. Neurol., 57, 745749.[CrossRef][Web of Science][Medline]
- Heron, S.E., Phillips, H.A., Mulley, J.C., Mazarib, A., Neufeld, M.Y., Berkovic, S.F. and Scheffer, I.E. (2004) Genetic variation of CACNA1H in idiopathic generalized epilepsy. Ann. Neurol., 55, 595596.[CrossRef][Web of Science][Medline]
-
Vitko, I., Chen, Y., Arias, J.M., Shen, Y., Wu, X.R. and Perez-Reyes, E. (2005) Functional characterization and neuronal modeling of the effects of childhood absence epilepsy variants of CACNA1H, a T-type calcium channel. J. Neurosci., 25, 48444855.
[Abstract/Free Full Text] - Emerick, M.C., Stein, R., Kunze, R., McNulty, M., Mitmann, S., Hanck, D.L. and Agnew, W.S. (2006) Profiling an array of Cav3.1 variants from the human T-type calcium channel CACNA1G: alternative structures, developmental expression and biophysical variations. Proteins Struct. Funct. Bioinform., in press.
-
Monteil, A., Chemin, J., Bourinet, E., Mennessier, G., Lory, P. and Nargeot, J. (2000) Molecular and functional properties of the human alpha(1G) subunit that forms T-type calcium channels. J. Biol. Chem., 275, 60906100.
[Abstract/Free Full Text] - Chemin, J., Monteil, A., Bourinet, E., Nargeot, J. and Lory, P. (2001) Alternatively spliced alpha(1G) (Ca(V)3.1) intracellular loops promote specific T-type Ca(2+) channel gating properties. Biophys. J., 80, 12381250.[Medline]
- Goetz, C. (1999) Textbook of Clinical Neurology. Saunders, Inc., Philadelphia, PA.
- Steriade, M. and Contreras, D. (1995) Relations between cortical and thalamic cellular events during transition from sleep patterns to paroxysmal activity. J. Neurosci., 15, 623642.[Abstract]
-
Meeren, H.K., Pijn, J.P., Van Luijtelaar, E.L., Coenen, A.M. and Lopes da Silva, F.H. (2002) Cortical focus drives widespread corticothalamic networks during spontaneous absence seizures in rats. J. Neurosci., 22, 14801495.
[Abstract/Free Full Text] - Coulter, D.A., Huguenard, J.R. and Prince, D.A. (1989) Characterization of ethosuximide reduction of low-threshold calcium current in thalamic neurons. Ann. Neurol., 25, 582593.[CrossRef][Web of Science][Medline]
- Coulter, D.A., Huguenard, J.R. and Prince, D.A. (1990) Differential effects of petit mal anticonvulsants and convulsants on thalamic neurones: calcium current reduction. Br. J. Pharmacol., 100, 800806.[Web of Science][Medline]
- Huguenard, J.R. and Prince, D.A. (1994) Intrathalamic rhythmicity studied in vitro: nominal T-current modulation causes robust antioscillatory effects. J. Neurosci., 14, 54855502.[Abstract]
- Manning, J.P., Richards, D.A., Leresche, N., Crunelli, V. and Bowery, N.G. (2004) Cortical-area specific block of genetically determined absence seizures by ethosuximide. Neuroscience, 123, 59.[CrossRef][Web of Science][Medline]
- Kim, D., Song, I., Keum, S., Lee, T., Jeong, M.J., Kim, S.S., McEnery, M.W. and Shin, H.S. (2001) Lack of the burst firing of thalamocortical relay neurons and resistance to absence seizures in mice lacking alpha(1G) T-type Ca(2+) channels. Neuron, 31, 3545.[CrossRef][Web of Science][Medline]
-
Song, I., Kim, D., Choi, S., Sun, M., Kim, Y. and Shin, H.S. (2004) Role of the alpha1G T-type calcium channel in spontaneous absence seizures in mutant mice. J. Neurosci., 24, 52495257.
[Abstract/Free Full Text] -
Zhang, Y., Vilaythong, A.P., Yoshor, D. and Noebels, J.L. (2004) Elevated thalamic low-voltage-activated currents precede the onset of absence epilepsy in the SNAP25-deficient mouse mutant coloboma. J. Neurosci., 24, 52395248.
[Abstract/Free Full Text] -
Zhang, Y., Mori, M., Burgess, D.L. and Noebels, J.L. (2002) Mutations in high-voltage-activated calcium channel genes stimulate low-voltage-activated currents in mouse thalamic relay neurons. J. Neurosci., 22, 63626371.
[Abstract/Free Full Text] - Tsakiridou, E., Bertollini, L., de Curtis, M., Avanzini, G. and Pape, H.C. (1995) Selective increase in T-type calcium conductance of reticular thalamic neurons in a rat model of absence epilepsy. J. Neurosci., 15, 31103117.[Abstract]
- Talley, E.M., Solorzano, G., Depaulis, A., Perez-Reyes, E. and Bayliss, D.A. (2000) Low-voltage-activated calcium channel subunit expression in a genetic model of absence epilepsy in the rat. Brain Res. Mol. Brain Res., 75, 159165.[Medline]
-
Perez-Reyes, E. (2003) Molecular physiology of low-voltage-activated t-type calcium channels. Physiol. Rev., 83, 117161.
[Abstract/Free Full Text] - Yunker, A.M. and McEnery, M.W. (2003) Low-voltage-activated (T-Type) calcium channels in review. J. Bioenerg. Biomembr., 35, 533575.[CrossRef][Web of Science][Medline]
-
Talley, E.M., Cribbs, L.L., Lee, J.H., Daud, A., Perez-Reyes, E. and Bayliss, D.A. (1999) Differential distribution of three members of a gene family encoding low voltage-activated (T-type) calcium channels. J. Neurosci., 19, 18951911.
[Abstract/Free Full Text] -
Cribbs, L.L., Lee, J.H., Yang, J., Satin, J., Zhang, Y., Daud, A., Barclay, J., Williamson, M.P., Fox, M., Rees, M. et al. (1998) Cloning and characterization of alpha1H from human heart, a member of the T-type Ca2+ channel gene family. Circ. Res., 83, 103109.
[Abstract/Free Full Text] - Williams, M.E., Washburn, M.S., Hans, M., Urrutia, A., Brust, P.F., Prodanovich, P., Harpold, M.M. and Stauderman, K.A. (1999) Structure and functional characterization of a novel human low-voltage activated calcium channel. J. Neurochem., 72, 791799.[CrossRef][Web of Science][Medline]
-
McRory, J.E., Santi, C.M., Hamming, K.S., Mezeyova, J., Sutton, K.G., Baillie, D.L., Stea, A. and Snutch, T.P. (2001) Molecular and functional characterization of a family of rat brain T-type calcium channels. J. Biol. Chem., 276, 39994011.
[Abstract/Free Full Text] -
Chemin, J., Monteil, A., Perez-Reyes, E., Bourinet, E., Nargeot, J. and Lory, P. (2002) Specific contribution of human T-type calcium channel isotypes (alpha(1G), alpha(1H) and alpha(1I)) to neuronal excitability. J. Physiol., 540, 314.
[Abstract/Free Full Text] - Park, J.Y., Jeong, S.W., Perez-Reyes, E. and Lee, J.H. (2003) Modulation of Ca(v)3.2 T-type Ca2+ channels by protein kinase C. FEBS Lett., 547, 3742.[CrossRef][Web of Science][Medline]
- Xia, M., Imredy, J.P., Santarelli, V.P., Liang, H.A., Condra, C.L., Bennett, P., Koblan, K.S. and Connolly, T.M. (2003) Generation and characterization of a cell line with inducible expression of Ca(v)3.2 (T-type) channels. Assay Drug Dev. Technol., 1, 637645.[CrossRef][Web of Science][Medline]
-
Welsby, P.J., Wang, H., Wolfe, J.T., Colbran, R.J., Johnson, M.L. and Barrett, P.Q. (2003) A mechanism for the direct regulation of T-type calcium channels by Ca2+/calmodulin-dependent kinase II. J. Neurosci., 23, 1011610121.
[Abstract/Free Full Text] -
Wolfe, J.T., Wang, H., Perez-Reyes, E. and Barrett, P.Q. (2002) Stimulation of recombinant Ca(v)3.2, T-type, Ca(2+) channel currents by CaMKIIgamma(C). J. Physiol., 538, 343355.
[Abstract/Free Full Text] - Wolfe, J.T., Wang, H., Howard, J., Garrison, J.C. and Barrett, P.Q. (2003) T-type calcium channel regulation by specific G-protein betagamma subunits. Nature, 424, 209213.[CrossRef][Medline]
- Heady, T.N., Gomora, J.C., Macdonald, T.L. and Perez-Reyes, E. (2001) Molecular pharmacology of T-type Ca2+ channels. Jpn. J. Pharmacol., 85, 339350.[CrossRef][Medline]
-
Jagannathan, S., Punt, E.L., Gu, Y., Arnoult, C., Sakkas, D., Barratt, C.L. and Publicover, S.J. (2002) Identification and localization of T-type voltage-operated calcium channel subunits in human male germ cells. Expression of multiple isoforms. J. Biol. Chem., 277, 84498456.
[Abstract/Free Full Text] - Liu, H.X., Cartegni, L., Zhang, M.Q. and Krainer, A.R. (2001) A mechanism for exon skipping caused by nonsense or missense mutations in BRCA1 and other genes. Nat. Genet., 27, 5558.[Web of Science][Medline]
- Mount, S.M. (2000) Genomic sequence, splicing, and gene annotation. Am. J. Hum. Genet., 67, 788792.[CrossRef][Web of Science][Medline]
-
Mount, S.M. (1982) A catalogue of splice junction sequences. Nucleic Acids Res., 10, 459472.
[Abstract/Free Full Text] - Mount, S.M. (1996) AT-AC introns: an ATtACk on dogma. Science, 271, 16901692.[CrossRef][Web of Science][Medline]
-
Fairbrother, W.G., Yeo, G.W., Yeh, R., Goldstein, P., Mawson, M., Sharp, P.A. and Burge, C.B. (2004) RESCUE-ESE identifies candidate exonic splicing enhancers in vertebrate exons. Nucleic Acids Res., 32, W187W190.
[Abstract/Free Full Text] -
Cartegni, L., Wang, J., Zhu, Z., Zhang, M.Q. and Krainer, A.R. (2003) ESEfinder: a web resource to identify exonic splicing enhancers. Nucleic Acids Res., 31, 35683571.
[Abstract/Free Full Text] - Mittman, S., Guo, J. and Agnew, W.S. (1999) Structure and alternative splicing of the gene encoding alpha1G, a human brain T calcium channel alpha1 subunit. Neurosci. Lett., 274, 143146.[CrossRef][Web of Science][Medline]
- Nadal-Ginard, B., Smith, C.W., Patton, J.G. and Breitbart, R.E. (1991) Alternative splicing is an efficient mechanism for the generation of protein diversity: contractile protein genes as a model system. Adv. Enzyme Regul., 31, 261286.[Web of Science][Medline]
- Latour, I., Gee, C.E., Robitaille, R. and Lacaille, J.C. (2001) Differential mechanisms of Ca2+ responses in glial cells evoked by exogenous and endogenous glutamate in rat hippocampus. Hippocampus, 11, 132145.[CrossRef][Web of Science][Medline]
- Hillman, R.T., Green, R.E. and Brenner, S.E. (2004) An unappreciated role for RNA surveillance. Genome Biol., 5, R8.[CrossRef][Medline]
- Mendell, J.T., Sharifi, N.A., Meyers, J.L., Martinez-Murillo, F. and Dietz, H.C. (2004) Nonsense surveillance regulates expression of diverse classes of mammalian transcripts and mutes genomic noise. Nat. Genet., 36, 10731078.[CrossRef][Web of Science][Medline]
-
Page, K.M., Heblich, F., Davies, A., Butcher, A.J., Leroy, J., Bertaso, F., Pratt, W.S. and Dolphin, A.C. (2004) Dominant-negative calcium channel suppression by truncated constructs involves a kinase implicated in the unfolded protein response. J. Neurosci., 24, 54005409.
[Abstract/Free Full Text] -
Raghib, A., Bertaso, F., Davies, A., Page, K.M., Meir, A., Bogdanov, Y. and Dolphin, A.C. (2001) Dominant-negative synthesis suppression of voltage-gated calcium channel Cav2.2 induced by truncated constructs. J. Neurosci., 21, 84958504.
[Abstract/Free Full Text] - Stuhmer, W., Conti, F., Suzuki, H., Wang, X.D., Noda, M., Yahagi, N., Kubo, H. and Numa, S. (1989) Structural parts involved in activation and inactivation of the sodium channel. Nature, 339, 597603.[CrossRef][Medline]
- Gargus, J.J. (2003) Unraveling monogenic channelopathies and their implications for complex polygenic disease. Am. J. Hum. Genet., 72, 785803.[CrossRef][Web of Science][Medline]
- Lopez, A.J. (1998) Alternative splicing of pre-mRNA: developmental consequences and mechanisms of regulation. Annu. Rev. Genet., 32, 279305.[CrossRef][Web of Science][Medline]
- Eng, L., Coutinho, G., Nahas, S., Yeo, G., Tanouye, R., Babaei, M., Dork, T., Burge, C. and Gatti, R.A. (2004) Nonclassical splicing mutations in the coding and noncoding regions of the ATM Gene: maximum entropy estimates of splice junction strengths. Hum. Mutat., 23, 6776.[CrossRef][Web of Science][Medline]
- Regan, M.R., Emerick, M.C. and Agnew, W.S. (2000) Full-length single-gene cDNA libraries: applications in splice variant analysis. Anal. Biochem., 286, 265276.[CrossRef][Web of Science][Medline]
-
Carninci, P., Kasukawa, T., Katayama, S., Gough, J., Frith, M.C., Maeda, N., Oyama, R., Ravasi, T., Lenhard, B., Wells, C. et al. (2005) The transcriptional landscape of the mammalian genome. Science, 309, 15591563.
[Abstract/Free Full Text]
This article has been cited by other articles:
![]() |
K. L. Powell, S. M. Cain, C. Ng, S. Sirdesai, L. S. David, M. Kyi, E. Garcia, J. R. Tyson, C. A. Reid, M. Bahlo, et al. A Cav3.2 T-Type Calcium Channel Point Mutation Has Splice-Variant-Specific Effects on Function and Segregates with Seizure Expression in a Polygenic Rat Model of Absence Epilepsy J. Neurosci., January 14, 2009; 29(2): 371 - 380. [Abstract] [Full Text] [PDF] |
||||
![]() |
W. B. Barbazuk, Y. Fu, and K. M. McGinnis Genome-wide analyses of alternative splicing in plants: Opportunities and challenges Genome Res., September 1, 2008; 18(9): 1381 - 1392. [Abstract] [Full Text] [PDF] |
||||
![]() |
J. Chemin, A. Mezghrani, I. Bidaud, S. Dupasquier, F. Marger, C. Barrere, J. Nargeot, and P. Lory Temperature-dependent Modulation of CaV3 T-type Calcium Channels by Protein Kinases C and A in Mammalian Cells J. Biol. Chem., November 9, 2007; 282(45): 32710 - 32718. [Abstract] [Full Text] [PDF] |
||||
![]() |
M. M. Jagodic, S. Pathirathna, M. T. Nelson, S. Mancuso, P. M. Joksovic, E. R. Rosenberg, D. A. Bayliss, V. Jevtovic-Todorovic, and S. M. Todorovic Cell-Specific Alterations of T-Type Calcium Current in Painful Diabetic Neuropathy Enhance Excitability of Sensory Neurons J. Neurosci., March 21, 2007; 27(12): 3305 - 3316. [Abstract] [Full Text] [PDF] |
||||
![]() |
J. Chemin, J. Nargeot, and P. Lory Chemical Determinants Involved in Anandamide-induced Inhibition of T-type Calcium Channels J. Biol. Chem., January 26, 2007; 282(4): 2314 - 2323. [Abstract] [Full Text] [PDF] |
||||
![]() |
I. Splawski, D. S. Yoo, S. C. Stotz, A. Cherry, D. E. Clapham, and M. T. Keating CACNA1H Mutations in Autism Spectrum Disorders J. Biol. Chem., August 4, 2006; 281(31): 22085 - 22091. [Abstract] [Full Text] [PDF] |
||||
![]() |
H. Khosravani and G. W. Zamponi Voltage-gated calcium channels and idiopathic generalized epilepsies. Physiol Rev, July 1, 2006; 86(3): 941 - 966. [Abstract] [Full Text] [PDF] |
||||
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||









ß or 


