- Research article
- Open Access
- Open Peer Review
Exploring the functional role of the CHRM2 gene in human cognition: results from a dense genotyping and brain expression study
BMC Medical Geneticsvolume 8, Article number: 66 (2007)
The CHRM2 gene, located on the long arm of chromosome 7 (7q31-35), is involved in neuronal excitability, synaptic plasticity and feedback regulation of acetylcholine release, and has been implicated in higher cognitive processing. The aim of this study is the identification of functional (non)coding variants underlying cognitive phenotypic variation.
We previously reported an association between polymorphisms in the 5'UTR regions of the CHRM2 gene and intelligence.. However, no functional variants within this area have currently been identified. In order to identify the relevant functional variant(s), we conducted a denser coverage of SNPs, using two independent Dutch cohorts, consisting of a children's sample (N = 371 ss; mean age 12.4) and an adult sample (N= 391 ss; mean age 37.6). For all individuals standardized intelligence measures were available. Subsequently, we investigated genotype-dependent CHRM2 gene expression levels in the brain, to explore putative enhancer/inhibition activity exerted by variants within the muscarinic acetylcholinergic receptor.
Using a test of within-family association two of the previously reported variants – rs2061174, and rs324650 – were again strongly associated with intelligence (P < 0.01). A new SNP (rs2350780) showed a trend towards significance. SNP rs324650, is located within a short interspersed repeat (SINE). Although the function of short interspersed repeats remains contentious, recent research revealed potential functionality of SINE repeats in a gene-regulatory context. Gene-expression levels in post-mortem brain material, however were not dependent on rs324650 genotype.
Using a denser coverage of SNPs in the CHRM2 gene, we confirmed the 5'UTR regions to be most interesting in the context of intelligence, and ruled out other regions of this gene. Although no correlation between genomic variants and gene expression was found, it would be interesting to examine allele-specific effects on CHRM2 transcripts expression in much more detail, for example in relation to transcripts specific halve-life and their relation to LTP and memory.
Identifying genes for variation in the range of normal intelligence could provide important clues to the genetic etiology of disturbed cognition in e.g. autism, reading disorder, and ADHD. Since the earliest 90's several groups have focussed on the identification – and subsequent replication – of common genetic polymorphisms underlying normal variation in cognitive abilities [1–5]. Among a handful of candidate genes that have been investigated in relation to normal cognitive variation as summarized in Posthuma & De Geus 2006 , the muscarinic 2 cholinergic receptor gene (CHRM2) has been consistently found to be associated with cognitive ability, and currently is the best replicated gene associated with general intelligence. A population-based association study conducted by Comings et al. (2003)  reported an association between a 3'UTR variant of the cholinergic muscarinic receptor 2 (CHRM2) gene explaining 1% of the variance in scores on full-scale IQ (FSIQ), and years of education. Suggestive evidence for linkage with performance IQ was found at 7q31-36, in the vicinity of the CHRM2 gene in a genome scan for intelligence based on 329 Australian families and 100 Dutch families, totalling 625 sib-pairs . We subsequently reported association between genetic variants within the CHRM2 gene and intelligence quotient (IQ) using two independent Dutch cohorts . This finding was then replicated by Dick and colleagues . All three association studies (Comings et al., 2003; Gosso et al., 2006; Dick et al., 2007) report significant association with IQ and non coding regions within in the CHRM2 gene (rs81919992 located in the 3' untranslated region (UTR) , and rs2061174 , and rs324650  in introns 4 and 5, respectively).
The CHRM2 gene belongs to the superfamily of G-protein-coupled receptors (GPCRs). Muscarinic acetylcholine receptors (M1-M5) activate a multitude of signaling pathways important for modulating neuronal excitability, synaptic plasticity and feedback regulation of acetylcholine (ACh) release [10, 11]. Combined behavioral and pharmacological animal studies involving M2 antagonists have shown the importance of cholinergic receptor activity for acquisition and retrieval of several learning tasks [12–16].
Despite its confirmed putative role in cognitive processes, further evidence for genetic regulatory variants on the CHRM2 gene have been difficult to assess, mainly due to its complex transcriptional expression patterns. Three different CHRM2 promoters have been reported based on work performed on different human cell lines . In combination with alternative splicing patterns this results in, at least, 6 different mRNA transcripts encoding for the same receptor protein (isoforms A till F)[17, 18]. Promoter activity for the CHRM2 gene was postulated to be tissue specific. The first promoter located upstream of exon 1, was preferentially used in cardiac cells (isoforms A and B); promoter 2 on intron 1 alternatively expressed on brain (isoforms C and D); and a third promoter located on intro 2 non-tissue specific (isoforms E and F). Independently, Zhou and coworkers  reported a fourth putative promoter region on intron 5, but this last result has not been independently confirmed yet . Although CHRM2 promoter usage is believed to be tissue specific, a single protein receptor is encoded. The functional significance of these transcripts is still unknown.
To fine-map the CHRM2 gene and to detect its functional role in cognitive ability, we genotyped a dense set of tag-SNPs within and flanking the CHRM2 gene in a sample of 762 Dutch individuals from 358 twin families belonging to two different age cohorts (mean ages 12.4 and 37.6). A family based genetic association test was used, which allows evaluating evidence for association free from spurious effects of population stratification [20–22]. In addition, gene expression assays were performed on brain controls to determine whether a significant correlation exists between the associated SNPs and CHRM2 gene expression levels.
All young and adult twins and their siblings were part of two larger cognitive studies and were recruited from the Netherlands Twin Registry [23, 24]. We have shown previously that the adult participants are representative of the Dutch population with respect to intelligence . Informed consent was obtained from the participants (adult cohort) or from their parents if they were under 18 (young cohort). The study was approved by the institutional review board of the VU University Medical Center. None of the individuals tested suffered from severe physical or mental handicaps, as assessed through surveys sent out to participants or their parents every two years.
The young cohort consisted of 177 twin pairs born between 1990 and 1992, and 55 siblings [6, 26], of which 371 were available for genotyping. Mean age of the genotyped twins was 12.4 (SD = 0.9) years of age and the siblings were between 8 and 15 years old at the time of testing. There were 35 monozygotic male twin pairs (MZM), 28 dizygotic male twin pairs (DZM), 48 monozygotic female twin pairs (MZF), 23 dizygotic female twin pairs (DZF), 26 dizygotic opposite-sex twin pairs (DOS), 24 male siblings and 24 female siblings, and 3 subjects form incomplete twin pairs (1 male, 2 females). Participation in this study included a voluntary agreement to provide buccal swabs for DNA extraction.
This sample is similar to the sample used in our initial analyses, except for twenty individuals that were deleted from analyses in the current sample due to additional genotyping and a more stringent threshold of genotyping failure per individual.
A total of 793 family members from 317 extended twin families participated in the adult cognition study . Participation in this study did not automatically include DNA collection, however, part of the sample (276 subjects) returned to the lab to provide blood samples. The sample characteristics have been described elsewhere . One hundred fifteen additional individuals provided buccal swabs via our biobanking project  for DNA extraction. Mean age of the total genotyped sample was 36.2 years (SD = 12.6). There were 25 monozygotic male twin pairs (MZM), 15 dizygotic male twin pairs (DZM), 1 DZM triplet, 20 monozygotic female twin pairs (MZF), 28 dizygotic female twin pairs (DZF) and 23 dizygotic opposite-sex twin pairs (DOS), 29 female siblings and 28 male siblings, and 109 subjects from incomplete twin pairs (41 males, 68 females).
In the young cohort, cognitive ability was assessed with the Dutch adaptation of the WISC-R , and consisted of four verbal subtests (similarities, vocabulary, arithmetic, and digit span) and two performance subtests (block design, and object assembly).
In the adult cohort, the Dutch adaptation of the WAISIII-R , assessed IQ and consisted of four verbal subtests (VIQ: information, similarities, vocabulary, and arithmetic) and four performance subtests (PIQ: picture completion, block design, matrix reasoning, and digit-symbol substitution). The correlation between verbal IQ and performance IQ is usually around 0.50 (0.53 in our data), implying that only 25% of the variance in PIQ and VIQ is shared. Thus, a substantial part of the variance in these two measures is non-overlapping, and theoretically they are expected to capture different aspects of cognitive ability. We therefore included VIQ and PIQ as measures of the two different aspects of intelligence as well as Full scale IQ (FSIQ) as a general measure of intelligence. In both cohorts, VIQ, PIQ and FSIQ were normally distributed, (see Table 1).
For both cohorts IQ scores standardized for the effects of age and sex were calculated. These were then z-transformed within cohorts to allow easy comparison across cohorts and across different tests.
DNA collection and isolation
Buccal swabs were collected from 371 children; DNA in adults was collected from blood samples in 391 adults. The DNA isolation from buccal swabs was performed using a cloroform/isopropanol extraction [31, 32]. DNA was extracted from blood samples using the salting out protocol described elsewhere . Zygosity was assessed using 11 highly polymorphic microsatellite markers (Heterozygosity > 0.80). Genotyping was performed blind to familial status and phenotypic data.
DNA and RNA extraction from tissue homogenates
Control brains from 50 individuals, 23 males with a mean age of 70.3 years (SD = 9.38), and 27 females with a mean age of 73.3 years (SD = 10.50) were obtained at autopsy from The Netherlands Brain Bank (NBB) . This material comes mainly from the superior and inferior parietal lobe. DNA isolation from 0.20 gram of frozen tissue was performed using the Puregene™ Kit (Gentra Systems, USA) according to standard protocol and doubled volume of all reagents per tissue weight. To verify DNA isolation, products were run on a 1% agarose gel.
Total RNA was isolated from 0.10 gram of frozen brain tissue with RNA-Bee™ following the manufacturer's recommendations (Isotex Diagnostics, Inc., USA). RNA was purified using the Qiagen RNeasy Mini kit (Qiagen Benelux B.V., The Netherlands) and verified on a 2% agarose gel. Five μg RNA was used to make cDNA using 200 U of Superscript™ III Reverse Transcriptase (Invitrogen, The Netherlands) in First Strand Buffer (Invitrogen, The Netherlands), 3.4 * 10-2 μg/μl random hexamer oligo, 3.4 * 10-2 μg/μl poly d(T) 12–18, 1.3 mM dNTPs, 1.1 μM DTT (Invitrogen, The Netherlands), 10 U RNaseOUT™ Ribonuclease Inhibitor Recombinant (Invitrogen, The Netherlands) and incubated two hours at 50°C. Subsequently, 20 U RNase H (Invitrogen, The Netherlands) was added and incubated 30 minutes at 37°C. Products were run on a 1% agarose gel to examine the quality.
Single nucleotide polymorphisms (SNPs) were selected using the information available from the International HapMap Project. SNP selection was based on a randomly selected population with northern and western European ancestry by the Centre d'Etude du polymorphisme Humain (CEPH) . The Minor Allele Frequency MAF had to be > 0.05 in order to exclude rare homozygous genotypes. Forty-two SNPs within the CHRM2 gene were thus selected from the CEPH population using Haploview version 3.32 (NCBI build 36.1).
SNP genotyping was performed using the SNPlex® assay platform. The SNPlex assay was conducted following the manufacturer's recommendations (Applied Biosystems, Foster city, CA, USA). All pre-PCR steps were performed on a cooled block. Reactions were carried out in Gene Amp 9700 Thermocycler (Applied Biosystems, Foster city, CA, USA). Data was analyzed using Genemapper v3.7 (Applied Biosystems, Foster city, CA, USA).
CHRM2 transcripts at brain level
Three different primer combinations were used to investigate the presence of CHRM2 transcript variants in normal brain controls. Forward primers FA&BGAGGCATCCAGGTCTCCAT, FC&DCGCAGCTCTCGCCA-GAGCCTT, and FE&FAAAGGACTCCTCGCTCCTTC were used in combination with a unique reverse primer RA-FCCCGATAATGGTCACCAAAC in order to tag isoforms A till F. PCR was performed at 94°C for 30 sec, 55°C for 30 sec, and 72°C for 1:30 min, for 40 cycles, followed by a 7 min extension at 72°C. To verify primers specificity PCR products were run on a 2% agarose gel.
Gene expression assay
RT-PCR was performed using specific primers encompassing the untranslated exon 5 (the last untranslated exon), which is present in all mRNA transcripts, and the coding sequence (CDS) of the CHRM2 gene; F-GAAACCAGCGACAGGTTTAAATG, R-GCTATTGTTAGAGGA-GTTTGTTGAGTTATTC. PCR was carried out at 94°C for 1 min, 64°C for 1 min, and 72°C for 1 min, for 40 cycles, followed by a 10 min extension at 72°C. Optimization of primer concentration and cDNA input was performed and dissociation curves for the selected primers obtained. Two housekeeping genes – β-actin and HPRT – were used as internal controls. RT-PCR reactions were performed twice independently, each time in duplicate.
Allele frequencies of all SNPs were estimated in both the children and adult cohorts using Haploview  in which a Hardy-Weinberg test is implemented, based on an exact calculation of the probability of observing a certain number of heterozygotes conditional on the number of copies of the minor SNP allele.
Genetic association tests were conducted using the program QTDT which implements the orthogonal model proposed by Abecasis et al., 2000  (see also Fulker et al., 1999; Posthuma et al., 2004 [21, 22]). This model allows one to decompose the genotypic effect into orthogonal between- (βb) and within- (βw) family components, and also models the residual sib-correlation as a function of polygenic or environmental factors. MZ twins can be included and are modelled as such, by adding zygosity status to the datafile. They are not informative to the within family component (unless they are paired with non-twin siblings), but are informative for the between family component. The between-family association component is sensitive to population admixture, whereas the within-family component is significant only in the presence of LD due to close linkage. The models used in QTDT take into account additive allelic between- and within family effects.
It is worth noting that, if population stratification acts to create a false association, the test for association using the within family component is still valid. More importantly, if population stratification acts to hide a genuine association, the test for association using the within family component has more power to detect this association than a population based association test. A significance level α of 0.01 was chosen.
Genotyping success rate was 95.36 (SD = 3.80) among both cohorts. Six tag-SNPs, (rs6957496, rs1424569, rs10488600, rs17494540, rs324582, and rs11773032), although with high genotyping rate, deviated from HWE (P < 0.05) despite a high genotype call rate. One tag-SNP, rs11773032 showed no variation in our population and was thus deleted from further analysis. LD parameters D' and r2 were obtained for all successfully genotyped SNPs. LD blocks were generated applying the algorithm defined by Gabriel et al., 2002  in which confidence bounds on D' are generated if 95% of the information shows "strong LD". By default, this method ignores markers with MAF < 0.05 (see Figure 1 and Table 2).
Two 5'UTR SNPs, previously reported, showed the strongest association with IQ, rs2061174 (intron 4) in the adult cohort and rs324650 (intron 5) in the young cohort  (see Figure 2). Within-family genetic effects were reflected in an increased IQ of 6.89 (PIQ) points for those individuals carrying the "A" allele of rs2061174 within the adult cohort. individuals in the young cohort bearing the "T" allele of rs324650 showed an increment of 5.30 IQ (VIQ) points (see Tables 3, 4 and 5). Interestingly, the most significant variant in the young cohort, rs324650, is part of a short interspersed repeat (SINE), namely a MIRb (mammalian-wide interspersed repeat) repeat of 160 bp long. The derived "T" allele contained in this repeat seems to be human-specific. In addition this MIRb repeat is also present in non-human primate linages – rhesus (macaca mulatta) and chimpanzee (pan troglodytes) – but not in other mammalian linages. Such an allele-specific effect may reflect that the variant is in LD with the causal allele, or that the "T" allele is directly modifying binding-properties of transcription starting sites (TSS) .
CHRM2 transcripts expression at brain level and correlations with CHRM2 tag-SNPs
Previous studies have shown that of the six known isoforms of CHRM2 only C and D are expressed in the brain [17, 18]. In contrast to this, we observed all six CHRM2 transcripts isoforms in brain material(data not shown).
After normalizing raw gene expression data to expression level of the housekeeping genes, no correlation between gene expression and CHRM2 gene genotypes for SNPs rs2061174, rs324640 or rs324650 was observed (data not shown).
Converging evidence from previous studies [7–9] has pointed to a role of the CHRM2 gene in intelligence. None of these studies, however, have identified the functional polymorphism explaining its role at a molecular level. The present study aimed to zoom in on the functional variants, by fine-mapping the most significant areas within this gene and also investigating differential brain expression as a function of different genotypes on the SNPs most strongly related to intelligence.
A total of 42 SNPs within the CHRM2 gene were genotyped in a young and adult cohort. Association analysis was conducted separately in both age cohorts to detect possible age dependent gene effects. Associations were found in different regions of the gene for each age cohort. Our current analyses showed that the same SNPs that were associated previously with intelligence, were again most significant, whereas a new SNP (rs2350780) showed a trend towards significance. Because of the dense coverage of SNPs used in this study, this confirms the importance of intron 4 and intron 5 regions, but rules out association with SNPs located elsewhere in the gene.
Four new SNPs in the intron 3 region, (rs2350780, rs1364409, rs7782965, and 1378646) showed association with PIQ in the adult cohort. These SNPs are in high LD (r2 between 0.58 – 0.72) between the most significant SNPs. SNP rs2350780 and rs2061174 were also found to be associated with intelligence by Dick and co-workers . These intronic SNPs are located 68 kb apart in introns 3 and 4, respectively. In our cohort, LD between these two variants is 0.58.
We found the most significant associations with PIQ in adults (rs2061174, χ2 = 9.14; P = 0.003) and with VIQ in children (rs323650, χ2 = 9.50; P = 0.002). Because only part of the variance in PIQ and VIQ is shared, and these results might reflect brain maturation processes and age-related genetic effects. Alternatively, the results could point to, and potentially explain, the genetic overlap between PIQ and VIQ, in which common genetic variants do not only interact modulating hippocampal neurotransmitter activity, but also and even more interesting from the epigenetic point of view, they might modulate priming and dendritic outgrowth underlying synaptic plasticity during embryogenesis  and at a post-natal stage , reflecting phenotypic variation at different IQ domains across the lifespan.
From a developmental perspective, brain maturation can be considered the most complex and dynamic lifelong process taking place in humans. Neuronal plasticity patterns (e.g. dendritic "pruning", synapse elimination, myelination) have been shown to vary significantly across life and among diverse brain structures (for a review see Toga et al., 2006 ). Variation in cognitive phenotypes may be the result of diverse allele-dependent effects that, although small in effect size, may contribute to cognitive phenotypes outcomes across life.
In situ hybridization experiments on mammals (e.g. mice)  have been of great utility to aid specific localization and interpretation of gene expression patterns. However, the localization of CHRM2 receptors transcripts has been conducted using probe sequences that did not distinguish between alternatively spliced transcripts. Our gene expression analyses showed that, in contrast to previously reported findings [17, 18], all six currently known transcripts (isoforms A till F) of the CHRM2 gene were present in brain tissue.
Our genotype-dependent CHRM2 expression, did not reveal functional significance of any of the SNPs that were significantly related to intelligence. However, one should keep in mind that at this point we were only able to study material from superior and inferior parietal lobe and further studies on other brain regions might give different results. Furthermore it would be of interest to examine allele-specific effects on CHRM2 transcripts expression in much more detail, for example in relation to transcripts specific halve-life and their relation to LTP and memory.
Although brain expression analysis did not reveal differential expression of CHRM2 transcripts, our study further zooms in on the CHRM2 gene, clearly confirming two regions of most importance to intelligence within introns 4 and 5. These regions are poorly conserved regions among relatively distant species, although they are conserved among primate species. Interestingly, the variant associated in the young cohort (rs324650) is located within a SINE repeat (MIRb). SINE repeats belongs to a wide family of transposable elements, which constitute the largest class of interspersed repeats that are found in our genome (12%) together with long interspersed repeats (LINE) an long terminal repeats (LTRs) . SINE repeats transpose through a RNA intermediate (reverse transcription process). All eukaryotic genomes contain mobile elements (retrosposable elements), although the proportion and activity of the classes of elements varies widely between genomes . The CHRM2 gene, like its G-protein receptor counterparts, shares the interestingly feature – at least form a functional perspective – of being an intronless protein , which is also observed among dopamine receptors , widely studied in relation to attention deficits.
Recent research has revealed a potential functionality of retroposons in a gene-regulatory context [38, 47–50]. It has been postulated that retroposon insertion processes may favour the generation of intronless proteins (for a review see Flavell 1995 and Brosius 2003 [51, 52]). If this hypothesis holds, the resulting intronless proteins are expected to contain exons among their 5'UTR region. Not surprisingly, among G-proteins with intronless open reading frames (ORFs), about 18% have been reported to contain untranslated exons on their 5'UTR [46, 53].
Multiple promoters and transcripts have been reported for the CHRM2 gene suggesting that the associated regions we identified harbour functional elements involved in regulation of transcription and/or alternative splicing [17–19]. Further investigation involving functional assays and non-coding polymorphisms may aid the search and subsequent identification of regulatory variants underlying normal cognitive variation.
Cherny SS, Cardon LR: General cognitive ability. Nature and nurture during middle childhood. Edited by: de Fries J, Plomin R, Fulker D. 1994, Oxford , Blackwell Publishers, 46-56.
Plomin R, Turic DM, Hill L, Turic DE, Stephens M, Williams J, Owen MJ, O'Donovan MC: A functional polymorphism in the succinate-semialdehyde dehydrogenase (aldehyde dehydrogenase 5 family, member A1) gene is associated with cognitive ability. Mol Psychiatry. 2004, 9 (6): 582-586. 10.1038/sj.mp.4001441.
Buyske S, Bates ME, Gharani N, Matise TC, Tischfield JA, Manowitz P: Cognitive traits link to human chromosomal regions. Behav Genet. 2006, 36 (1): 65-76. 10.1007/s10519-005-9008-9.
Posthuma D, Luciano M, Geus EJ, Wright MJ, Slagboom PE, Montgomery GW, Boomsma DI, Martin NG: A Genomewide Scan for Intelligence Identifies Quantitative Trait Loci on 2q and 6p. Am J Hum Genet. 2005, 77 (2): 318-326. 10.1086/432647.
Dick DM, Aliev F, Bierut L, Goate A, Rice J, Hinrichs A, Bertelsen S, Wang JC, Dunn G, Kuperman S, Schuckit M, Nurnberger J, Porjesz B, Beglieter H, Kramer J, Hesselbrock V: Linkage Analyses of IQ in the Collaborative Study on the Genetics of Alcoholism (COGA) Sample. Behav Genet. 2006, 36 (1): 77-86. 10.1007/s10519-005-9009-8.
Polderman TJC, Gosso MF, Posthuma D, Van Beijsterveldt TC, Heutink P, Verhulst FC, Boomsma DI: A longitudinal twin study on IQ, executive functioning, and attention problems during childhood and early adolescence. Acta Neurol Belg. 2006, 106 (4): 191-207.
Comings DE, Wu S, Rostamkhani M, McGue M, Lacono WG, Cheng LS, MacMurray JP: Role of the cholinergic muscarinic 2 receptor (CHRM2) gene in cognition. Mol Psychiatry. 2003, 8 (1): 10-11. 10.1038/sj.mp.4001095.
Gosso MF, van Belzen MJ, de Geus EJC, Polderman JC, Heutink P, Boomsma DI, Posthuma D: Association between the CHRM2 gene and intelligence in a sample of 304 Dutch families. Genes, Brain and Behavior. 2006, Blackwell Publishing, 5: 577-584. 10.1111/j.1601-183X.2006.00211.x.
Dick DM, Aliev F, Kramer J, Wang JC, Hinrichs A, Bertelsen S, Kuperman S, Schuckit M, Nurnberger J, Edenberg HJ, Porjesz B, Begleiter H, Hesselbrock V, Goate A, Bierut L: Association of CHRM2 with IQ: Converging Evidence for a Gene Influencing Intelligence. Behav Genet. 2007, 37 (2): 265-272. 10.1007/s10519-006-9131-2.
Levey AI, Kitt CA, Simonds WF, Price DL, Brann MR: Identification and localization of muscarinic acetylcholine receptor proteins in brain with subtype-specific antibodies. J Neurosci. 1991, 11 (10): 3218-3226.
Volpicelli LA, Levey AI: Muscarinic acetylcholine receptor subtypes in cerebral cortex and hippocampus. Prog Brain Res. 2004, 145: 59-66.
Orsetti M, Casamenti F, Pepeu G: Enhanced acetylcholine release in the hippocampus and cortex during acquisition of an operant behavior. Brain Res. 1996, 724 (1): 89-96. 10.1016/0006-8993(96)00292-2.
Carey GJ, Billard W, Binch H, Cohen-Williams M, Crosby G, Grzelak M, Guzik H, Kozlowski JA, Lowe DB, Pond AJ, Tedesco RP, Watkins RW, Coffin VL: SCH 57790, a selective muscarinic M(2) receptor antagonist, releases acetylcholine and produces cognitive enhancement in laboratory animals. Eur J Pharmacol. 2001, 431 (2): 189-200. 10.1016/S0014-2999(01)01440-6.
Accili D, Fishburn CS, Drago J, Steiner H, Lachowicz JE, Park BH, Gauda EB, Lee EJ, Cool MH, Sibley DR, Gerfen CR, Westphal H, Fuchs S: A targeted mutation of the D3 dopamine receptor gene is associated with hyperactivity in mice. Proc Natl Acad Sci U S A. 1996, 93 (5): 1945-1949. 10.1073/pnas.93.5.1945.
Seeger T, Fedorova I, Zheng F, Miyakawa T, Koustova E, Gomeza J, Basile AS, Alzheimer C, Wess J: M2 muscarinic acetylcholine receptor knock-out mice show deficits in behavioral flexibility, working memory, and hippocampal plasticity. J Neurosci. 2004, 24 (45): 10117-10127. 10.1523/JNEUROSCI.3581-04.2004.
Gautam D, Duttaroy A, Cui Y, Han S, Deng C, Seeger T, Alzheimer C, Wess J: M1-M3 muscarinic acetylcholine receptor-deficient mice: novel phenotypes. J Mol Neur. 2006, 30: 157-160. 10.1385/JMN:30:1:157.
Fenech AG, Billington CK, Swan C, Richards S, Hunter T, Ebejer MJ, Felice AE, Ellul-Micallef R, Hall IP: Novel polymorphisms influencing transcription of the human CHRM2 gene in airway smooth muscle. Am J Respir Cell Mol Biol. 2004, 30 (5): 678-686. 10.1165/rcmb.2003-0011OC.
Krejci A, Bruce AW, Dolezal V, Tucek S, Buckley NJ: Multiple promoters drive tissue-specific expression of the human M muscarinic acetylcholine receptor gene. J Neurochem. 2004, 91 (1): 88-98. 10.1111/j.1471-4159.2004.02694.x.
Zhou C, Fryer AD, Jacoby DB: Structure of the human M(2) muscarinic acetylcholine receptor gene and its promoter. Gene. 2001, 271 (1): 87-92. 10.1016/S0378-1119(01)00494-2.
Abecasis GR, Cardon LR, Cookson WO: A general test of association for quantitative traits in nuclear families. Am J Hum Genet. 2000, 66 (1): 279-292. 10.1086/302698.
Fulker DW, Cherny SS, Sham PC, Hewitt JK: Combined linkage and association sib-pair analysis for quantitative traits. Am J Hum Genet. 1999, 64 (1): 259-267. 10.1086/302193.
Posthuma D, de Geus EJ, Boomsma DI, Neale MC: Combined linkage and association tests in mx. Behav Genet. 2004, 34 (2): 179-196. 10.1023/B:BEGE.0000013732.19486.74.
Boomsma DI, de Geus EJ, Vink JM, Stubbe JH, Distel MA, Hottenga JJ, Posthuma D, van Beijsterveldt TC, Hudziak JJ, Bartels M, Willemsen G: Netherlands Twin Register: from twins to twin families. Twin Res Hum Genet. 2006, 9 (6): 849-857. 10.1375/twin.9.6.849.
Bartels M, Rietveld MJ, Van Baal GC, Boomsma DI: Genetic and environmental influences on the development of intelligence. Behav Genet. 2002, 32 (4): 237-249. 10.1023/A:1019772628912.
Posthuma D, de Geus EJ, Boomsma DI: Perceptual speed and IQ are associated through common genetic factors. Behav Genet. 2001, 31 (6): 593-602. 10.1023/A:1013349512683.
Polderman TJC, Stins JF, Posthuma D, Gosso MF, Verhulst FC, Boomsma DI: The phenotypic and genotypic relation between working memory speed and capacity. Intelligence. 2006, 34: 549-560. 10.1016/j.intell.2006.03.010.
Gosso MF, van Belzen M, de Geus EJ, Polderman JC, Heutink P, Boomsma DI, Posthuma D: Association between the CHRM2 gene and intelligence in a sample of 304 Dutch families. Genes Brain Behav. 2006, 5 (8): 577-584. 10.1111/j.1601-183X.2006.00211.x.
Hoekstra C, Meijer P, Kluft C, Heutink P, Smit G, de Geus E, Smit JH, van Bruggen A, Montgomery GW, Boomsma DI: Genetics of dizygotic twinning: a feasibility study for a biobank. Twin Res. 2004, 7 (6): 556-563. 10.1375/1369052042663751.
Wechsler D: Wechsler Intelligence Scale for Children-Revised (WISC-R). 1986, Lisse , Swets & Zeitlinger
Wechsler D: WAIS-III Wechsler Adult Intelligence Scale . 1997, San Antonio, Texas , Psychological Corporation
Meulenbelt I, Droog S, Trommelen GJ, Boomsma DI, Slagboom PE: High-yield noninvasive human genomic DNA isolation method for genetic studies in geographically dispersed families and populations. Am J Hum Genet. 1995, 57 (5): 1252-1254.
Min JL, Lakenberg N, Bakker-Verweij M, Suchiman E, Boomsma DI, Slagboom PE, Meulenbelt I: High microsatellite and SNP genotyping success rates established in a large number of genomic DNA samples extracted from mouth swabs and genotypes. Twin Res Hum Genet. 2006, 9: 501-506. 10.1375/twin.9.4.501.
Miller SA, Dykes DD, Polesky HF: A simple salting out procedure for extracting DNA from human nucleated cells. Nucleic Acids Res. 1988, 16 (3): 1215-10.1093/nar/16.3.1215.
Nederland Brain Bank. [http://www.brainbank.nl]
Centre d’Etude du polymorphisme Humain. [http://www.hapmap.org/thehapmap.html.en]
Gabriel SB, Schaffner SF, Nguyen H, Moore JM, Roy J, Blumenstiel B, Higgins J, DeFelice M, Lochner A, Faggart M, Liu-Cordero SN, Rotimi C, Adeyemo A, Cooper R, Ward R, Lander ES, Daly MJ, Altshuler D: The structure of haplotype blocks in the human genome. Science. 2002, 296 (5576): 2225-2229. 10.1126/science.1069424.
Mei L, Yan JC: [Advances in subclassifications of muscarinic acetylcholine receptors]. Sheng Li Ke Xue Jin Zhan. 1985, 16 (3): 248-252.
Niculescu MD, Craciunescu CN, Zeisel SH: Dietary choline deficiency alters global and gene-specific DNA methylation in the developing hippocampus of mouse fetal brains. Faseb J. 2006, 20 (1): 43-49. 10.1096/fj.05-4707com.
Olivera S, Henley JM, Rodriguez-Ithurralde D: AMPA receptor potentiation by acetylcholinesterase is age-dependently upregulated at synaptogenesis sites of the rat brain. Int J Dev Neurosci. 2003, 21 (1): 49-61. 10.1016/S0736-5748(02)00083-7.
Toga AW, Thompson PM, Sowell ER: Mapping brain maturation. Trends Neurosci. 2006, 29 (3): 148-159. 10.1016/j.tins.2006.01.007.
Lein ES, Hawrylycz MJ, Ao N, Ayres M, Bensinger A, Bernard A, Boe AF, Boguski MS, Brockway KS, Byrnes EJ, Chen L, Chen L, Chen TM, Chi Chin M, Chong J, Crook BE, Czaplinska A, Dang CN, Datta S, Dee NR, Desaki AL, Desta T, Diep E, Dolbeare TA, Donelan MJ, Dong HW, Dougherty JG, Duncan BJ, Ebbert AJ, Eichele G, Estin LK, Faber C, Facer BA, Fields R, Fischer SR, Fliss TP, Frensley C, Gates SN, Glattfelder KJ, Halverson KR, Hart MR, Hohmann JG, Howell MP, Jeung DP, Johnson RA, Karr PT, Kawal R, Kidney JM, Knapik RH, Kuan CL, Lake JH, Laramee AR, Larsen KD, Lau C, Lemon TA, Liang AJ, Liu Y, Luong LT, Michaels J, Morgan JJ, Morgan RJ, Mortrud MT, Mosqueda NF, Ng LL, Ng R, Orta GJ, Overly CC, Pak TH, Parry SE, Pathak SD, Pearson OC, Puchalski RB, Riley ZL, Rockett HR, Rowland SA, Royall JJ, Ruiz MJ, Sarno NR, Schaffnit K, Shapovalova NV, Sivisay T, Slaughterbeck CR, Smith SC, Smith KA, Smith BI, Sodt AJ, Stewart NN, Stumpf KR, Sunkin SM, Sutram M, Tam A, Teemer CD, Thaller C, Thompson CL, Varnam LR, Visel A, Whitlock RM, Wohnoutka PE, Wolkey CK, Wong VY, Wood M, Yaylaoglu MB, Young RC, Youngstrom BL, Feng Yuan X, Zhang B, Zwingman TA, Jones AR: Genome-wide atlas of gene expression in the adult mouse brain. Nature. 2007, 445 (7124): 168-10.1038/nature05453.
Smit AF: Interspersed repeats and other mementos of transposable elements in mammalian genomes. Curr Opin Genet Dev. 1999, 9 (6): 657-663. 10.1016/S0959-437X(99)00031-3.
Kidwell MG, Lisch DR: Transposable elements and host genome evolution. Trends In Ecology And Evolution. 2000, 15 (3): 95-99. 10.1016/S0169-5347(99)01817-0.
Fredriksson R, Lagerstrom MC, Lundin LG, Schioth HB: The G-Protein-Coupled Receptors in the Human Genome Form Five Main Families. Phylogenetic Analysis, Paralogon Groups, and Fingerprints. Mol Pharmacol. 2003, 63 (6): 1256-1272. 10.1124/mol.63.6.1256.
O'Dowd BF: Structures of dopamine receptors. J Neurochem. 1993, 60 (3): 804-816. 10.1111/j.1471-4159.1993.tb03224.x.
Bejerano G, Lowe CB, Ahituv N, King B, Siepel A, Salama SR, Rubin EM, Kent WJ, Haussler D: A distal enhancer and an ultraconserved exon are derived from a novel retroposon. Nature. 2006, 441 (7089): 87-90. 10.1038/nature04696.
Han JS, Boeke JD: LINE-1 retrotransposons: modulators of quantity and quality of mammalian gene expression?. Bioessays. 2005, 27 (8): 775-784. 10.1002/bies.20257.
Hellmann-Blumberg U, Hintz MF, Gatewood JM, Schmid CW: Developmental differences in methylation of human Alu repeats. Mol Cell Biol. 1993, 13 (8): 4523-4530.
Muratani K, Hada T, Yamamoto Y, Kaneko T, Shigeto Y, Ohue T, Furuyama J, Higashino K: Inactivation of the cholinesterase gene by Alu insertion: possible mechanism for human gene transposition. Proc Natl Acad Sci U S A. 1991, 88 (24): 11315-11319. 10.1073/pnas.88.24.11315.
Flavell AJ: Retroelements, reverse transcriptase and evolution. Comp Biochem Physiol B Biochem Mol Biol. 1995, 110 (1): 3-15. 10.1016/0305-0491(94)00122-B.
Brosius J: The contribution of RNAs and retroposition to evolutionary novelties. Genetica. 2003, 118 (2-3): 99-116. 10.1023/A:1024141306559.
Rouquier S, Taviaux S, Trask BJ, Brand-Arpon V, van den Engh G, Demaille J, Giorgi D: Distribution of olfactory receptor genes in the human genome. Nat Genet. 1998, 18 (3): 243-250. 10.1038/ng0398-243.
Lee CJ, Irizarry K: Alternative splicing in the nervous system: an emerging source of diversity and regulation. Biological Psychiatry. 2003, 54 (8): 771-10.1016/S0006-3223(03)00375-5.
The pre-publication history for this paper can be accessed here:http://www.biomedcentral.com/1471-2350/8/66/prepub
Supported by the Universitair Stimulerings Fonds (grant number 96/22), the Human Frontiers of Science Program (grant number rg0154/1998-B), the Netherlands Organization for Scientific Research (NWO) grants 904-57-94 and NWO/SPI 56-464-14192. DP is supported by GenomEUtwin grant (EU/QLRT-2001-01254) and by NWO/MaGW Vernieuwingsimpuls 016-065-318. This study was supported by the Centre for Medical Systems Biology (CMSB), a centre of excellence approved by the Netherlands Genomics Initiative/Netherlands Organization for Scientific Research (NWO). We thank all the twin families participating in this study as well as the Dutch Brain Bank for samples, Saskia van Mil and David Sondervan for technical support and Dina Ruano for valuable comments.
The author(s) declare that they have no competing interests.
MFG conducted the SNP selection and genotyping. MFG and DP performed the statistical analyses. DNA was provided by DIB. Phenotypic data was provided by DIB, EJC, TJC and DP. MFG drafted the manuscript under DP and PH supervision. DP and PH supervised the study. All authors read and approved the final manuscript.