Analysis of candidate genes for morphine preference quantitative trait locus Mop2

Compared to DBA/2J (D2), C57BL/6J (B6) inbred mice exhibit strong morphine preference when tested using a two-bottle choice drinking paradigm. A morphine preference quantitative trait locus (QTL), Mop2, was originally mapped to proximal chromosome (Chr) 10 using a B6xD2 F2 intercross population, confirmed with reciprocal congenic strains and fine mapped with recombinant congenic strains. These efforts identified a ∼ 10-Million base pair (Mbp) interval, underlying Mop2, containing 35 genes. To further reduce the interval, mice from the D2.B6-Mop2-P1 congenic strain were backcrossed to parental D2 mice and two new recombinant strains of interest were generated: D2.B6-Mop2-P1.pD.dB and D2.B6-Mop2-P1.pD.dD. Results obtained from testing these strains in the two-bottle choice drinking paradigm suggest that the gene(s) responsible for the Mop2 QTL is one or more of 22 remaining within the newly defined interval (∼ 7.6 Mbp) which includes Oprm1 and several other genes related to opioid pharmacology. Real-time qRT-PCR analysis of Oprm1 and opioid-related genes Rgs17, Ppp1r14c, Vip, and Iyd revealed both between-strain and within-strain expression differences in comparisons of saline- and morphine-treated B6 and D2 mice. Analysis of Rgs17 protein levels also revealed both between-strain and within-strain differences in comparisons of saline- and morphine-treated B6 and D2 mice. Results suggest that the Mop2 QTL represents the combined influence of multiple genetic variants on morphine preference in these two strains. Relative contributions of each variant remain to be determined.

Quantitative trait loci for electrical seizure threshold mapped in C57BLKS/J and C57BL/10SnJ mice

We mapped the quantitative trait loci (QTL) that contribute to the robust difference in maximal electroshock seizure threshold (MEST) between C57BLKS/J (BKS) and C57BL10S/J (B10S) mice. BKS, B10S, BKS × B10S F1 and BKS × B10S F2 intercross mice were tested for MEST at 8-9 weeks of age. Results of F2 testing showed that, in this cross, MEST is a continuously distributed trait determined by polygenic inheritance. Mice from the extremes of the trait distribution were genotyped using microarray technology. MEST correlated significantly with body weight and sex; however, because of the high correlation between these factors, the QTL mapping was conditioned on sex alone. A sequential series of statistical analyses was used to map QTLs including single-point, multipoint and multilocus methods. Two QTLs reached genome-wide levels of significance based upon an empirically determined permutation threshold: chromosome 6 (LOD = 6.0 at ∼69 cM) and chromosome 8 (LOD = 5.7 at ∼27 cM). Two additional QTLs were retained in a multilocus regression model: chromosome 3 (LOD = 2.1 at ∼68 cM) and chromosome 5 (LOD = 2.7 at ∼73 cM). Together the four QTLs explain one third of the total phenotypic variance in the mapping population. Lack of overlap between the major MEST QTLs mapped here in BKS and B10S mice and those mapped previously in C57BL/6J and DBA/2J mice (strains that are closely related to BKS and B10S) suggest that BKS and B10S represent a new polygenic mouse model for investigating susceptibility to seizures.

Confirmation of multiple seizure susceptibility QTLs on chromosome 15 in C57BL/6J and DBA/2J inbred mice

To confirm seizure susceptibility (SZS) quantitative trait loci (QTLs) on chromosome (chr) 15 identified previously using C57BL/6J (B6) and DBA/2J (D2) mice and to refine their genomic map position, we studied a set of three congenic strains in which overlapping segments of chr 15 from D2 were transferred onto the B6 background. We measured thresholds for generalized electroshock seizure (GEST) and maximal electroshock seizure (MEST) in congenic strains and B6-like littermates and also tested their responses to kainic acid (KA) and pentylenetetrazol (PTZ). Results document that MEST is significantly lower in strains 15M and 15D, which harbor medial and distal (telomeric) segments of chr 15 (respectively) from D2, compared with strain 15P, which harbors the proximal (acromeric) segment of chr 15 from D2, and with control littermates. Congenic strains 15P and 15M exhibited greater KA SZS compared with strain 15D and B6-like controls. All congenic strains were similar to controls with regard to PTZ SZS. Taken together, results suggest there are multiple SZS QTLs on chr 15 and that two QTLs harbor gene variants that affect MEST and KA SZS independently. The MEST QTL is refined to a 19 Mb region flanked by rs13482630 and D15Mit159. This interval contains 350 genes, 183 of which reside in areas where the polymorphism rate between B6 and D2 is high. The KA QTL interval spans a 65 Mb region flanked by markers D15Mit13 and rs31271969. It harbors 83 genes in highly polymorphic areas, 310 genes in all. Complete dissection of these loci will lead to identification of genetic variants that influence SZS in mice and provide a better understanding of seizure biology.

Alternative transcripts and evidence of imprinting of GNAL on 18p11.2

Genetic studies implicating the region of human chromosome 18p11.2 in susceptibility to bipolar disorder and schizophrenia have observed parent-of-origin effects that may be explained by genomic imprinting. We have identified a transcriptional variant of the GNAL gene in this region, employing an alternative first exon that is 5' to the originally identified start site. This alternative GNAL transcript encodes a longer functional variant of the stimulatory G-protein alpha subunit, Golf. The isoforms of Golf display different expression patterns in the CNS and functionally couple to the dopamine D1 receptor when heterologously expressed in Sf9 cells. In addition, there are CpG islands in the vicinity of both first exons that are differentially methylated, a hallmark of genomic imprinting. These results suggest that GNAL, and possibly other genes in the region, is subject to epigenetic regulation and strengthen the case for a susceptibility gene in this region.

Association between variation in the human KCNJ10 potassium ion channel gene and seizure susceptibility

PURPOSE

Our research program uses genetic linkage and association analysis to identify human seizure sensitivity and resistance alleles. Quantitative trait loci mapping in mice led to identification of genetic variation in the potassium ion channel gene Kcnj10, implicating it as a putative seizure susceptibility gene. The purpose of this work was to translate these animal model data to a human genetic association study.

METHODS

We used single stranded conformation polymorphism (SSCP) electrophoresis, DNA sequencing and database searching (NCBI) to identify variation in the human KCNJ10 gene. Restriction fragment length polymorphism (RFLP) analysis, SSCP and Pyrosequencing were used to genotype a single nucleotide polymorphism (SNP, dbSNP rs#1130183) in KCNJ10 in epilepsy patients (n = 407) and unrelated controls (n = 284). The epilepsy group was comprised of patients with refractory mesial temporal lobe epilepsy (n = 153), childhood absence (n = 84), juvenile myoclonic (n = 111) and idiopathic generalized epilepsy not otherwise specified (IGE-NOS, n = 59) and all were of European ancestry.

RESULTS

SNP rs#1130183 (C > T) alters amino acid 271 (of 379) from an arginine to a cysteine (R271C). The C allele (Arg) is common with conversion to the T allele (Cys) occurring twice as often in controls compared to epilepsy patients. Contingency analysis documented a statistically significant association between seizure resistance and allele frequency, Mantel-Haenszel chi square = 5.65, d.f. = 1, P = 0.017, odds ratio 0.52, 95% CI 0.33-0.82.

CONCLUSION

The T allele of SNP rs#1130183 is associated with seizure resistance when common forms of focal and generalized epilepsy are analyzed as a group. These data suggest that this missense variation in KCNJ10 (or a nearby variation) is related to general seizure susceptibility in humans.

Quantitative genetic study of maximal electroshock seizure threshold in mice: evidence for a major seizure susceptibility locus on distal chromosome 1

We conducted a quantitative trait locus (QTL) mapping study to dissect the multifactorial nature of maximal electroshock seizure threshold (MEST) in C57BL/6 (B6) and DBA/2 (D2) mice. MEST determination involved a standard paradigm in which 8- to 12-week-old mice received one shock per day with a daily incremental increase in electrical current until a maximal seizure (tonic hindlimb extension) was induced. Mean MEST values in parental strains were separated by over five standard deviation units, with D2 mice showing lower values than B6 mice. The distribution of MEST values in B6xD2 F2 intercrossed mice spanned the entire phenotypic range defined by parental strains. Statistical mapping yielded significant evidence for QTLs on chromosomes 1, 2, 5, and 15, which together explained over 60% of the phenotypic variance in the model. The chromosome 1 QTL represents a locus of major effect, accounting for about one-third of the genetic variance. Experiments involving a congenic strain (B6.D2-Mtv7(a)/Ty) enabled more precise mapping of the chromosome 1 QTL and indicate that it lies in the genetic interval between markers D1Mit145 and D1Mit17. These results support the hypothesis that the distal portion of chromosome 1 harbors a gene(s) that has a fundamental role in regulating seizure susceptibility.

Lack of association between an interleukin 1 beta (IL-1beta) gene variation and refractory temporal lobe epilepsy

PURPOSE

We attempted to confirm recent findings of Kanemoto et al. that demonstrated a positive association (p < 0.017) between a polymorphism in the promoter region of the interleukin 1-beta (IL-1beta) gene and the clinical phenotype of temporal lobe epilepsy with hippocampal sclerosis (TLE+HS).

METHODS

We determined the frequency of this polymorphism in a group of 61 TLE+HS patients of European ancestry and compared it with that found in 119 ethnically matched control subjects.

RESULTS

Analysis of genotype and allele frequencies showed no statistically significant difference in the distribution of the polymorphism between the two groups (p = 0.10).

CONCLUSIONS

These data suggest that this IL-1beta promoter polymorphism does not act as a strong susceptibility factor for TLE+HS in a population of individuals of European ancestry.

Human G(olf) gene polymorphisms and vulnerability to bipolar disorder

Two intronic polymorphisms of the human alpha subunit of the olfactory G-protein (G(olf)) are described. They were detected with single-stranded conformational polymorphism (SSCP) methods and confirmed by sequencing both strands. These single base pair (bp) substitutions occur in introns 3 (an A/G at 35 bp 3' from the exon 3/intron 3 5' splice site) and 10 (an T/G at 7 bp 5' from the 3' splice site). Both polymorphisms are relatively common, with minor allele frequencies of 31% (intron 3) and 16% (intron 10). The intron 3 variant shows no linkage disequilibrium with an intron 5 (CA)n microsatellite located approximately 50 kb 3' from the intron 3 variant, among a small group of German individuals with schizophrenia. The intron 3 variant is interesting because it may create an 'in-frame' cryptic splice site which, if activated, would add 12 residues to exon 3. The intron 10 variant is interesting because a purine is substituted for a pyrimidine in the 'polypyrimidine' tract of the 3' splice site, a single base substitution of the type which has been associated with aberrant splicing in the androgen receptor gene.

Lack of association between temporal lobe epilepsy and a novel polymorphism in the alpha 2 subunit gene (ATP1A2) of the sodium potassium transporting ATPase

Genetic linkage studies in rodents and humans have identified specific chromosomal regions harboring seizure susceptibility genes. We have identified a novel polymorphism in the human alpha 2 subunit gene (ATP1A2) of the sodium potassium transporting ATPase (NaK-pump), a candidate gene for human temporal lobe epilepsy (TLE) based on its chromosomal location and function in ion homeostasis. The polymorphism consists of a four base pair insertion 12 base pairs upstream of the start of exon 2. We performed an association study between this polymorphism and TLE. Our study did not find a significant difference in the frequency of this polymorphism between TLE patients and controls, indicating that this variation is not a major susceptibility factor. However, since the number of patients studied so far is small and the functional consequence of the polymorphism is unknown, the variation may yet be found to play a minor role in increased risk for seizure susceptibility. In contrast to the findings in TLE patients and controls, we did find a significant difference in the frequency of the variation between African Americans and persons of European descent. This finding demonstrates the potential effect of population stratification on studies of this type and supports the growing use of parental and familial samples for controls in association studies. Further study of this polymorphism is warranted as it may be involved in other disease processes for which there are known ethnic-specific susceptibilities. Am. J. Med. Genet. (Neuropsychiatr. Genet.) 96:79-83, 2000.

Hypoxic repression of lactate dehydrogenase-B in retina

In order to elucidate mechanisms controlling hypoxia induced gene repression in the retina, we studied expression of the lactate dehydrogenase gene LDH-B. RT/PCR was used to isolate cDNA fragments of the LDH-A and LDH-B genes from both rat and chick. Northern analysis was used to measure LDH-B mRNA expression and in situ hybridization to localize LDH-B mRNA during development of the avascular chick retina. Primary retinal cultures of rat and chick were subjected to hypoxia and reperfusion paradigms, and levels of LDH-A and LDH-B mRNA expression were measured by Northern and semi-quantitative RT/PCR analyses. Northern analysis demonstrated that LDH-B mRNA is developmentally regulated in the chick retina according to a time course that correlates with oxygen availability. In developing chick retina in situ hybridization localized LDH-B mRNA to aerobic regions of the retina, coincident with the distribution of active mitochondria. Northern and semi-quantitative RT/PCR analyses demonstrated that LDH-B mRNA levels decrease in primary rat and chick retinal cells in response to hypoxia. LDH-B mRNA increased to control levels in retinal cells exposed to hypoxia then reperfused with oxygen, while the opposite was true for LDH-A, an hypoxia inducible gene. These data demonstrate that LDH-B gene repression after hypoxia and reactivation after oxygen reperfusion occurs in both vascular and avascular retinal cells. Oxygen regulated expression of LDH-B is opposite and complementary to that of LDH-A and provides a system to elucidate the mechanisms underlying hypoxic gene repression and reactivation after reperfusion. Mechanisms controlling gene repression and reactivation based on oxygen availability are likely to be critically involved in ischemic damage and neovascular retinopathy.

Mapping loci for pentylenetetrazol-induced seizure susceptibility in mice

DBA/2J (D2) and C57BL/6J (B6) mice exhibit differential sensitivity to seizures induced by various chemical and physical methods, with D2 mice being relatively sensitive and B6 mice relatively resistant. We conducted studies in mature D2, B6, F1, and F2 intercross mice to investigate behavioral seizure responses to pentylenetetrazol (PTZ) and to map the location of genes that influence this trait. Mice were injected with PTZ and observed for 45 min. Seizure parameters included latencies to focal clonus, generalized clonus, and maximal seizure. Latencies were used to calculate a seizure score that was used for quantitative mapping. F2 mice (n = 511) exhibited a wide range of latencies with two-thirds of the group expressing maximal seizure. Complementary statistical analyses identified loci on proximal (near D1Mit11) and distal chromosome 1 (near D1Mit17) as having the strongest and most significant effects in this model. Another locus of significant effect was detected on chromosome 5 (near D5Mit398). Suggestive evidence for additional PTZ seizure-related loci was detected on chromosomes 3, 4, and 6. Of the seizure-related loci identified in this study, those on chromosomes 1 (distal), 4, and 5 map close to loci previously identified in a similar F2 population tested with kainic acid. Results document that the complex genetic influences controlling seizure response in B6 and D2 mice are partially independent of the nature of the chemoconvulsant stimulus with a locus on distal chromosome 1 being of fundamental importance.

Genotyping microsatellite polymorphisms by agarose gel electrophoresis with ethidium bromide staining: application to quantitative trait loci analysis of seizure susceptibility in mice

Agarose gel electrophoresis with ethidium bromide staining (AGE/EBS) is an efficient and reliable method for analyzing microsatellite polymorphisms. We report the use of AGE/EBS for analyzing DNA microsatellite polymorphisms in a preliminary quantitative trait loci (QTL) study of seizure susceptibility in which a candidate gene strategy was used to direct initial mapping efforts. F2 intercross progeny, derived from seizure-sensitive DBA/2J (D2) and seizure-resistant C57BL/6J (B6) inbred strains of mice, were tested for their sensitivity to the seizure-inducing effect of pentylenetetrazol (PTZ), a gamma-aminobutyric acid (GABA) receptor antagonist. A semi-automated method is described, in which DNA microsatellites were amplified by polymerase chain reaction (PCR) to yield products of 100-200 base pair (bp) in length. Alleles were separated on 3-6% MetaPhor agarose gels, stained with ethidium bromide, and visualized by ultraviolet (UV) illumination. Univariate analysis of genotype and phenotype data provides evidence for a seizure-related QTL on chromosome 5, near genes coding for the GABAA receptor subunits alpha 5 and gamma 3. Interestingly, this suggestive QTL derives from the more resistant B6 strain, but it nonetheless provides impetus for the characterization of possible strain differences in these two candidate genes. Overall, these results demonstrate that AGE/EBS can be useful for rapid screening of genomic regions of special interest in QTL mapping studies.

Genetic influences on electrical seizure threshold

C57BL/6J (B6) and DBA/2J (D2) mice have been characterized previously as seizure-resistant and seizure-sensitive, respectively, a distinction based primarily upon a differential response to the convulsant effects of various drugs. In the present study, electroconvulsive shock (ECS) was used to assess maximal electroshock threshold (MET) in B6, D2 and hybrid mice. Results revealed that D2 mice have a significantly lower MET compared to B6 mice. There was also a significant gender effect for B6 and F2 mice with females exhibiting a lower MET compared to males. METs for F1 and F2 intercross mice were intermediate between the two parental strains. The difference in variance between F2 and F1 generation mice indicated that about three-quarters of the total variance is due to genetic influence. Taken together, results of this study suggest that the large difference in MET between B6 and D2 mice is a highly heritable trait which may yield to genetic dissection through use of quantitative trait locus mapping.

Cloning of murine CDK9/PITALRE and its tissue-specific expression in development

The cdc2-family of serine/threonine kinases and their binding partners recently were implicated in developmental roles. We previously cloned a cdc2-related kinase, cdk9/PITALRE, that is able to phosphorylate the retinoblastoma protein in vitro. We describe here the cloning and the characterization of the mouse homolog of cdk9/PITALRE. The murine cDNA is 98% identical with humans and is expressed at high levels in brain and kidney tissues. The kinase activity and protein expression of cdk9/PITALRE were highest in terminally differentiated tissues such as the muscle and brain. In situ immunohistology and immunofluorescence detected cdk9/PITALRE protein not only within terminally differentiated cells such as muscle and neuronal cells, but also in proliferating cells. C2C12 and P19 cells induced to differentiate along muscle and neural lineages peaked in cdk9/PITALRE kinase activity at the end of differentiation. These results suggest that, among other roles, cdk9/PITALRE plays a role not unlike cdk5 in the differentiation of certain cell types.

The effect of ecdysis on DNA of the hepatopancreas and green gland of the Florida spiny lobster (Panulirus argus)

The DNA content per gram of wet hepatopancreas and green gland of Panulirus argus was found to vary with molt stage. In hepatopancreas of premolt and postmolt spiny lobsters, the DNA content was, respectively, 5-fold and 3-fold higher than in intermolt animals. Green gland DNA content showed a trend of varying in the same manner as hepatopancreas DNA. The DNA content per gram of wet green gland from intermolt lobsters averaged about 50 times higher than that of hepatopancreas. The changes in DNA content were not related to changes in wet tissue weights during the molting cycle. The ratio of hepatopancreas or green gland wet weights to total body weight remained constant. A ribosomal DNA probe hybridized to intermolt hepatopancreas DNA samples to a much greater extent than to premolt or postmolt samples, suggesting that rDNA sequences are under-represented in the amplified DNA which occurs prior to ecdysis.

Transient expression of RSVCAT in transgenic zebrafish made by electroporation

We report the fish use of an exponential decay electroporation system to introduce foreign DNA into fertilized zebrafish embryos. The plasmid RSVCAT (Rous sarcoma viral promoter (RSV) upstream from the chloramphenicol acetyltransferase gene (CAT)) was linearized and introduced into fertile zebrafish embryos by electroporation no later than the four-cell stage. Conditions for the procedure were empirically derived, and 68% of the treated animals survived through hatching to at least 6 days after fertilization and well beyond. Dot-blot analysis on DNA extracted from individual hatching fry demonstrated that 65% of the animals tested carried the foreign construct. Enzyme assays on the soluble proteins of treated animals were positive for chloramphenicol acetyltransferase activity. These data demonstrate that the foreign construct was being transiently expressed in the developing tissues of the embryo. The simplicity of this technique will greatly enhance the ability to analyze gene promoter regulation in vivo in transgenic zebrafish. The ability of the electroporated DNA to integrate into the host genome and to generate stable lines of transgenic fish is discussed.

Molecular analyses of carbonic anhydrase-II expression and regulation in the developing chicken lens

The expression of carbonic anhydrase-II (CA-II) in the developing chicken lens was examined and compared with that in the retina of the chicken embryo. CA-II expression was measured by immunohistochemistry and radioimmunoassay during development, and CA-II mRNA was quantified by Northern blot and densitometric scanning and localized by in situ hybridization. A functional promoter of the chicken CA-II gene was identified by transfection of primary embryonic chicken lens epithelial cells and analyzed in deletion mutants. The results establish that CA-II makes up about 0.1% of the total soluble protein of the embryonic chicken lens, an amount insufficient to make it a candidate for an enzyme crystallin in this species. Lens fiber differentiation coincided with a loss of CA-II mRNA and protein; by contrast, CA-II persisted in the epithelial cells of the embryonic and mature lens. This and previous studies showed that CA-II amounts to as much as 3% of the protein of the embryonic chicken retina and follows a different developmental time course of expression; like the lens, CA-II decreases until day 10 in the embryonic retina, but, unlike the lens, it increases thereafter and plateaus at hatching. Progressive deletions of the 5' flanking regions (from position -1314 to +32) of the CA-II gene fused to the bacterial chloramphenicol acetyltransferase (CAT) reporter gene resulted in a gradual loss of promoter activity, consistent with an additive effect of putative cis-regulatory elements found in many crystallin genes. These experiments provide the foundation for a molecular analysis of the developmental and differential regulation of the CA-II gene in lens and retina.

Changes in distribution of mitochondria in the developing chick retina

The distribution of mitochondria in the developing chick retina was examined by enzyme histochemistry, immunohistochemistry, and electron microscopy. Two distinct phenomena were observed: (1) progressive segregation of mitochondria in specific locations in the developing tissue; and (2) progressive loss of mitochondrial activity from the inner retina as it matures. Densitometric scans of stained tissue sections were used to quantitate the relative amounts of mitochondrial activity in the retinal layers. Mitochondria were localized to tissue regions by transmission electron microscopy, enzyme histochemistry for the inner membrane bound mitochondrial enzymes succinic dehydrogenase and cytochrome oxidase, a combined histochemical and ultrastructural method for cytochrome oxidase, and immunolocalization of the mitochondrial matrix enzyme glutamate dehydrogenase. Seven-day embryonic chick retina has a high number of mitochondrial structures and a high level of activity. As development proceeds, the mitochondria organize into layers within the tissue. However, the relative activity of mitochondria in much of the inner retina decreases. In the post-hatch retina, 50% of the mitochondrial activity is found in less than 10% of the tissue area, in the inner segments of the photoreceptor cells.

Changes in expression and distribution of lactate dehydrogenase isoenzymes in the developing chick retina

We have investigated the expression and distribution of lactate dehydrogenase (LDH) isoenzymes in the developing and adult chicken retina. Non-denaturing polyacrylamide gel electrophoresis was used to follow the appearance and expression of the five main LDH isoenzymes in tissue homogenates. Immunohistochemistry was used to define the distribution of the aerobic heart type LDH and the anaerobic muscle type LDH in paraffin sections of embryonic and adult chick retina. The electrophoretic results show that the expression of the anaerobic isoenzyme increases in the retina as development proceeds. Immunoreactivity against the anaerobic isoenzyme localizes to the inner plexiform layer and other regions of the inner chick retina. The aerobic isoenzyme immunoreactivity is localized to the inner segments of the photoreceptor cells as well as the ganglion cell layer. Both antibodies bind weakly in other layers of the retina with the notable exception that the anaerobic protein does not localize in the photoreceptor cell inner segments. These results provide further evidence for the anaerobic nature of the adult inner retina. As retinal development proceeds, the expression of the anaerobic isoenzyme in the retina increases.