51
|
Kember RL, Fernandes C, Tunbridge EM, Liu L, Payá-Cano JL, Parsons MJ, Schalkwyk LC. A B2 SINE insertion in the Comt1 gene (Comt1(B2i)) results in an overexpressing, behavior modifying allele present in classical inbred mouse strains. GENES BRAIN AND BEHAVIOR 2011; 9:925-32. [PMID: 20618449 PMCID: PMC3037563 DOI: 10.1111/j.1601-183x.2010.00614.x] [Citation(s) in RCA: 17] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
Catechol-O-methyltransferase (COMT) is a key enzyme for dopamine catabolism and COMT is a candidate gene for human psychiatric disorders. In mouse it is located on chromosome 16 in a large genomic region of extremely low variation among the classical inbred strains, with no confirmed single nucleotide polymorphisms (SNPs) between strains C57BL/6J and DBA/2J within a 600-kB window. We found a B2 SINE in the 3′ untranslated region (UTR) of Comt1 which is present in C57BL/6J (Comt1B2i) and other strains including 129 (multiple sublines), but is not found in DBA/2J (Comt1+) and many other strains including wild-derived Mus domesticus, M. musculus, M. molossinus, M.castaneus and M. spretus. Comt1B2i is absent in strains closely related to C57BL/6, such as C57L and C57BR, indicating that it was polymorphic in the cross that gave rise to these strains. The strain distribution of Comt1B2i indicates a likely origin of the allele in the parental Lathrop stock. A stringent association test, using 670 highly outbred mice (Boulder Heterogeneous Stock), indicates that this insertion allele may be responsible for a difference in behavior related to exploration. Gene expression differences at the mRNA and enzyme activity level (1.7-fold relative to wild type) indicate a mechanism for this behavioral effect. Taken together, these findings show that Comt1B2i (a B2 SINE insertion) results in a relatively modest difference in Comt1 expression and enzyme activity (comparable to the human Val-Met polymorphism) which has a demonstrable behavioral phenotype across a variety of outbred genetic backgrounds.
Collapse
|
52
|
Schosser A, Pirlo K, Gaysina D, Cohen-Woods S, Schalkwyk LC, Elkin A, Korszun A, Gunasinghe C, Gray J, Jones L, Meaburn E, Farmer AE, Craig IW, McGuffin P. Utility of the pooling approach as applied to whole genome association scans with high-density Affymetrix microarrays. BMC Res Notes 2010; 3:274. [PMID: 21040578 PMCID: PMC2984392 DOI: 10.1186/1756-0500-3-274] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/14/2010] [Accepted: 11/01/2010] [Indexed: 11/20/2022] Open
Abstract
Background We report an attempt to extend the previously successful approach of combining SNP (single nucleotide polymorphism) microarrays and DNA pooling (SNP-MaP) employing high-density microarrays. Whereas earlier studies employed a range of Affymetrix SNP microarrays comprising from 10 K to 500 K SNPs, this most recent investigation used the 6.0 chip which displays 906,600 SNP probes and 946,000 probes for the interrogation of CNVs (copy number variations). The genotyping assay using the Affymetrix SNP 6.0 array is highly demanding on sample quality due to the small feature size, low redundancy, and lack of mismatch probes. Findings In the first study published so far using this microarray on pooled DNA, we found that pooled cheek swab DNA could not accurately predict real allele frequencies of the samples that comprised the pools. In contrast, the allele frequency estimates using blood DNA pools were reasonable, although inferior compared to those obtained with previously employed Affymetrix microarrays. However, it might be possible to improve performance by developing improved analysis methods. Conclusions Despite the decreasing costs of genome-wide individual genotyping, the pooling approach may have applications in very large-scale case-control association studies. In such cases, our study suggests that high-quality DNA preparations and lower density platforms should be preferred.
Collapse
|
53
|
Meaburn EL, Schalkwyk LC, Mill J. Allele-specific methylation in the human genome: implications for genetic studies of complex disease. Epigenetics 2010; 5:578-82. [PMID: 20716955 DOI: 10.4161/epi.5.7.12960] [Citation(s) in RCA: 101] [Impact Index Per Article: 7.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/19/2022] Open
Abstract
Across the genome, outside of a small number of known imprinted genes and regions subject to X-inactivation in females, DNA methylation at CpG dinucleotides is often assumed to be complementary across both alleles in a diploid cell. However, recent findings suggest the reality is more complex, with the discovery that allele-specific methylation (ASM) is a common feature across the genome. A key observation is that the majority of ASM is associated with genetic variation in cis, although a noticeable proportion is also non-cis in nature and mediated, for example, by parental origin. ASM appears to be both quantitative, characterized by subtle skewing of DNA methylation between alleles, and heterogeneous, varying across tissues and between individuals. These findings have important implications for complex disease genetics; whilst cis-mediated ASM provides a functional consequence for non-coding genetic variation, heterogeneous and quantitative ASM complicates the identification of disease-associated loci. We propose that non-cis ASM could contribute toward the 'missing heritability' of complex diseases, rendering certain loci hemizygous and masking the direct association between genotype and phenotype. We suggest that the interpretation of results from genomewide association studies can be improved by the incorporation of epi-allelic information, and that in order to fully understand the extent and consequence of ASM in the human genome, a comprehensive sequencing-based analysis of allelic methylation patterns across tissues and individuals is required.
Collapse
|
54
|
Davies MN, Meaburn EL, Schalkwyk LC. Gene set enrichment; a problem of pathways. Brief Funct Genomics 2010; 9:385-90. [PMID: 20861160 DOI: 10.1093/bfgp/elq021] [Citation(s) in RCA: 5] [Impact Index Per Article: 0.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/19/2022] Open
Abstract
Gene Set Enrichment (GSE) is a computational technique which determines whether a priori defined set of genes show statistically significant differential expression between two phenotypes. Currently, the gene sets used for GSE are derived from annotation or pathway databases, which often contain computationally based and unrepresentative data. Here, we propose a novel approach for the generation of comprehensive and biologically derived gene sets, deriving sets through the application of machine learning techniques to gene expression data. These gene sets can be produced for specific tissues, developmental stages or environments. They provide a powerful and functionally meaningful way in which to mine genomewide association and next generation sequencing data in order to identify disease-associated variants and pathways.
Collapse
|
55
|
Docherty SJ, Davis OSP, Kovas Y, Meaburn EL, Dale PS, Petrill SA, Schalkwyk LC, Plomin R. A genome-wide association study identifies multiple loci associated with mathematics ability and disability. GENES, BRAIN, AND BEHAVIOR 2010; 9:234-47. [PMID: 20039944 PMCID: PMC2855870 DOI: 10.1111/j.1601-183x.2009.00553.x] [Citation(s) in RCA: 59] [Impact Index Per Article: 4.2] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Subscribe] [Scholar Register] [Received: 07/01/2009] [Revised: 09/08/2009] [Accepted: 11/02/2009] [Indexed: 12/01/2022]
Abstract
Numeracy is as important as literacy and exhibits a similar frequency of disability. Although its etiology is relatively poorly understood, quantitative genetic research has demonstrated mathematical ability to be moderately heritable. In this first genome-wide association study (GWAS) of mathematical ability and disability, 10 out of 43 single nucleotide polymorphism (SNP) associations nominated from two high- vs. low-ability (n = 600 10-year-olds each) scans of pooled DNA were validated (P < 0.05) in an individually genotyped sample of (*)2356 individuals spanning the entire distribution of mathematical ability, as assessed by teacher reports and online tests. Although the effects are of the modest sizes now expected for complex traits and require further replication, interesting candidate genes are implicated such as NRCAM which encodes a neuronal cell adhesion molecule. When combined into a set, the 10 SNPs account for 2.9% (F = 56.85; df = 1 and 1881; P = 7.277e-14) of the phenotypic variance. The association is linear across the distribution consistent with a quantitative trait locus (QTL) hypothesis; the third of children in our sample who harbour 10 or more of the 20 risk alleles identified are nearly twice as likely (OR = 1.96; df = 1; P = 3.696e-07) to be in the lowest performing 15% of the distribution. Our results correspond with those of quantitative genetic research in indicating that mathematical ability and disability are influenced by many genes generating small effects across the entire spectrum of ability, implying that more highly powered studies will be needed to detect and replicate these QTL associations.
Collapse
|
56
|
Schalkwyk LC, Meaburn EL, Smith R, Dempster EL, Jeffries AR, Davies MN, Plomin R, Mill J. Allelic skewing of DNA methylation is widespread across the genome. Am J Hum Genet 2010; 86:196-212. [PMID: 20159110 DOI: 10.1016/j.ajhg.2010.01.014] [Citation(s) in RCA: 195] [Impact Index Per Article: 13.9] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/13/2009] [Revised: 12/21/2009] [Accepted: 01/18/2010] [Indexed: 12/28/2022] Open
Abstract
DNA methylation is assumed to be complementary on both alleles across the genome, although there are exceptions, notably in regions subject to genomic imprinting. We present a genome-wide survey of the degree of allelic skewing of DNA methylation with the aim of identifying previously unreported differentially methylated regions (DMRs) associated primarily with genomic imprinting or DNA sequence variation acting in cis. We used SNP microarrays to quantitatively assess allele-specific DNA methylation (ASM) in amplicons covering 7.6% of the human genome following cleavage with a cocktail of methylation-sensitive restriction enzymes (MSREs). Selected findings were verified using bisulfite-mapping and gene-expression analyses, subsequently tested in a second tissue from the same individuals, and replicated in DNA obtained from 30 parent-child trios. Our approach detected clear examples of ASM in the vicinity of known imprinted loci, highlighting the validity of the method. In total, 2,704 (1.5%) of our 183,605 informative and stringently filtered SNPs demonstrate an average relative allele score (RAS) change > or =0.10 following MSRE digestion. In agreement with previous reports, the majority of ASM ( approximately 90%) appears to be cis in nature, and several examples of tissue-specific ASM were identified. Our data show that ASM is a widespread phenomenon, with >35,000 such sites potentially occurring across the genome, and that a spectrum of ASM is likely, with heterogeneity between individuals and across tissues. These findings impact our understanding about the origin of individual phenotypic differences and have implications for genetic studies of complex disease.
Collapse
|
57
|
Smith RG, Kember RL, Mill J, Fernandes C, Schalkwyk LC, Buxbaum JD, Reichenberg A. Advancing paternal age is associated with deficits in social and exploratory behaviors in the offspring: a mouse model. PLoS One 2009; 4:e8456. [PMID: 20041141 PMCID: PMC2794376 DOI: 10.1371/journal.pone.0008456] [Citation(s) in RCA: 70] [Impact Index Per Article: 4.7] [Reference Citation Analysis] [Abstract] [MESH Headings] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/28/2009] [Accepted: 12/02/2009] [Indexed: 01/16/2023] Open
Abstract
Background Accumulating evidence from epidemiological research has demonstrated an association between advanced paternal age and risk for several psychiatric disorders including autism, schizophrenia and early-onset bipolar disorder. In order to establish causality, this study used an animal model to investigate the effects of advanced paternal age on behavioural deficits in the offspring. Methods C57BL/6J offspring (n = 12 per group) were bred from fathers of two different ages, 2 months (young) and 10 months (old), and mothers aged 2 months (n = 6 breeding pairs per group). Social and exploratory behaviors were examined in the offspring. Principal Findings The offspring of older fathers were found to engage in significantly less social (p = 0.02) and exploratory (p = 0.02) behaviors than the offspring of younger fathers. There were no significant differences in measures of motor activity. Conclusions Given the well-controlled nature of this study, this provides the strongest evidence for deleterious effects of advancing paternal age on social and exploratory behavior. De-novo chromosomal changes and/or inherited epigenetic changes are the most plausible explanatory factors.
Collapse
|
58
|
Meaburn EL, Fernandes C, Craig IW, Plomin R, Schalkwyk LC. Assessing individual differences in genome-wide gene expression in human whole blood: reliability over four hours and stability over 10 months. Twin Res Hum Genet 2009; 12:372-80. [PMID: 19653838 DOI: 10.1375/twin.12.4.372] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/30/2023]
Abstract
Studying the causes and correlates of natural variation in gene expression in healthy populations assumes that individual differences in gene expression can be reliably and stably assessed across time. However, this is yet to be established. We examined 4-hour test-retest reliability and 10 month test-retest stability of individual differences in gene expression in ten 12-year-old children. Blood was collected on four occasions: 10 a.m. and 2 p.m. on Day 1 and 10 months later at 10 a.m. and 2 p.m. Total RNA was hybridized to Affymetrix-U133 plus 2.0 arrays. For each probeset, the correlation across individuals between 10 a.m. and 2 p.m. on Day 1 estimates test-retest reliability. We identified 3,414 variable and abundantly expressed probesets whose 4-hour test-retest reliability exceeded .70, a conventionally accepted level of reliability, which we had 80% power to detect. Of the 3,414 reliable probesets, 1,752 were also significantly reliable 10 months later. We assessed the long-term stability of individual differences in gene expression by correlating the average expression level for each probe-set across the two 4-hour assessments on Day 1 with the average level of each probe-set across the two 4-hour assessments 10 months later. 1,291 (73.7%) of the 1,752 probe-sets that reliably detected individual differences across 4 hours on two occasions, 10 months apart, also stably detected individual differences across 10 months. Heritability, as estimated from the MZ twin intraclass correlations, is twice as high for the 1,752 reliable probesets versus all present probesets on the array (0.68 vs 0.34), and is even higher (0.76) for the 1,291 reliable probesets that are also stable across 10 months. The 1,291 probesets that reliably detect individual differences from a single peripheral blood collection and stably detect individual differences over 10 months are promising targets for research on the causes (e.g., eQTLs) and correlates (e.g., psychopathology) of individual differences in gene expression.
Collapse
|
59
|
Ronald A, Butcher LM, Docherty SJ, Davis OSP, Schalkwyk LC, Craig IW, Plomin R. Response to comment by Stuart Macgregor. Behav Genet 2009. [DOI: 10.1007/s10519-009-9319-3] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
|
60
|
Ronald A, Butcher LM, Docherty S, Davis OSP, Schalkwyk LC, Craig IW, Plomin R. A genome-wide association study of social and non-social autistic-like traits in the general population using pooled DNA, 500 K SNP microarrays and both community and diagnosed autism replication samples. Behav Genet 2009; 40:31-45. [PMID: 20012890 PMCID: PMC2797846 DOI: 10.1007/s10519-009-9308-6] [Citation(s) in RCA: 37] [Impact Index Per Article: 2.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/21/2009] [Accepted: 10/14/2009] [Indexed: 10/28/2022]
Abstract
Two separate genome-wide association studies were conducted to identify single nucleotide polymorphisms (SNPs) associated with social and nonsocial autistic-like traits. We predicted that we would find SNPs associated with social and non-social autistic-like traits and that different SNPs would be associated with social and nonsocial. In Stage 1, each study screened for allele frequency differences in approximately 430,000 autosomal SNPs using pooled DNA on microarrays in high-scoring versus low-scoring boys from a general population sample (N = approximately 400/group). In Stage 2, 22 and 20 SNPs in the social and non-social studies, respectively, were tested for QTL association by individually genotyping an independent community sample of 1,400 boys. One SNP (rs11894053) was nominally associated (P < .05, uncorrected for multiple testing) with social autistic-like traits. When the sample was increased by adding females, 2 additional SNPs were nominally significant (P < .05). These 3 SNPs, however, showed no significant association in transmission disequilibrium analyses of diagnosed ASD families.
Collapse
|
61
|
Kõks S, Velthut A, Sarapik A, Altmäe S, Reinmaa E, Schalkwyk LC, Fernandes C, Lad HV, Soomets U, Jaakma U, Salumets A. The differential transcriptome and ontology profiles of floating and cumulus granulosa cells in stimulated human antral follicles. Mol Hum Reprod 2009; 16:229-40. [PMID: 19933312 DOI: 10.1093/molehr/gap103] [Citation(s) in RCA: 46] [Impact Index Per Article: 3.1] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/13/2022] Open
Abstract
Communication between various ovarian cell types is a prerequisite for folliculogenesis and ovulation. In antral follicles granulosa cells divide into two distinct populations of mural and cumulus granulosa cells (CGC), enveloping the antrum and surrounding the oocyte, respectively. Both cell types, with the mural compartment in excess, contribute to the floating granulosa cell (FGC) population in the follicular fluid. The aim of this study was to compare the transcriptomes of FGC and CGC in stimulated antral follicles obtained from 19 women undergoing IVF-ICSI procedure. FGC were obtained from follicular fluid during the follicle puncture procedure and CGC were acquired after oocyte denudation for micromanipulation. Gene expression analysis was conducted using the genome-wide Affymetrix transcriptome array. The expression profile of the two granulosa cell populations varied significantly. Out of 28 869 analysed transcripts 4480 were differentially expressed (q-value < 10(-4)) and 489 showed > or =2-fold difference in the expression level with 222 genes up-regulated in FGC and 267 in CGC. The transcriptome of FGC showed higher expression of genes involved in immune response, hematological system function and organismal injury, although CGC had genes involved in protein degradation and nervous system function up-regulated. Cell-to-cell signalling and interaction pathways were noted in both cell populations. Furthermore, numerous novel transcripts that have not been previously described in follicular physiology were identified. In conclusion, our results provide a solid basis for future studies in follicular biology that will help to identify molecular markers for oocyte and embryo viability in IVF.
Collapse
|
62
|
Overall RW, Kempermann G, Peirce J, Lu L, Goldowitz D, Gage FH, Goodwin S, Smit AB, Airey DC, Rosen GD, Schalkwyk LC, Sutter TR, Nowakowski RS, Whatley S, Williams RW. Genetics of the hippocampal transcriptome in mouse: a systematic survey and online neurogenomics resource. Front Neurosci 2009; 3:55. [PMID: 20582282 PMCID: PMC2858614 DOI: 10.3389/neuro.15.003.2009] [Citation(s) in RCA: 58] [Impact Index Per Article: 3.9] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/01/2009] [Accepted: 10/26/2009] [Indexed: 11/13/2022] Open
Abstract
Differences in gene expression in the CNS influence behavior and disease susceptibility. To systematically explore the role of normal variation in expression on hippocampal structure and function, we generated an online microarray database for a diverse panel of strains of mice, including most common inbred strains and numerous recombinant inbred lines (www.genenetwork.org). Using this resource, coexpression networks for families of genes can be generated rapidly to test causal models related to function. The data set is optimized for quantitative trait locus (QTL) mapping and was used to identify over 5500 QTLs that modulate mRNA levels. We describe a wide variety of analyses and novel synthetic approaches that take advantage of this resource, and demonstrate how both the data and associated tools can be applied to the study of gene regulation in the hippocampus and relations to structure and function.
Collapse
|
63
|
Davies MN, Lawn S, Whatley S, Fernandes C, Williams RW, Schalkwyk LC. To What Extent is Blood a Reasonable Surrogate for Brain in Gene Expression Studies: Estimation from Mouse Hippocampus and Spleen. Front Neurosci 2009; 3:54. [PMID: 20582281 PMCID: PMC2858613 DOI: 10.3389/neuro.15.002.2009] [Citation(s) in RCA: 10] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/12/2009] [Accepted: 09/17/2009] [Indexed: 01/08/2023] Open
Abstract
Microarrays are designed to measure genome-wide differences in gene expression. In cases where a tissue is not accessible for analysis (e.g. human brain), it is of interest to determine whether a second, accessible tissue could be used as a surrogate for transcription profiling. Surrogacy has applications in the study of behavioural and neurodegenerative disorders. Comparison between hippocampus and spleen mRNA obtained from a mouse recombinant inbred panel indicates a high degree of correlation between the tissues for genes that display a high heritability of expression level. This correlation is not limited to apparent expression differences caused by sequence polymorphisms in the target sequences and includes both cis and trans genetic effects. A tissue such as blood could therefore give surrogate information on expression in brain for a subset of genes, in particular those co-expressed between the two tissues, which have heritably varying expression.
Collapse
|
64
|
Kas MJH, Gelegen C, Schalkwyk LC, Collier DA. Interspecies comparisons of functional genetic variations and their implications in neuropsychiatry. Am J Med Genet B Neuropsychiatr Genet 2009; 150B:309-17. [PMID: 18561257 DOI: 10.1002/ajmg.b.30815] [Citation(s) in RCA: 19] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Submit a Manuscript] [Subscribe] [Scholar Register] [Indexed: 01/30/2023]
Abstract
Animal studies are important for the identification and functional characterization of the biological substrates underlying complex psychiatric disorders. However, novel insights into the relationship between the genome and behavior are needed for the development of fully valid translational models. Based on the notion that in different species, the same genes may independently give rise to alleles with similar functional and phenotypic effects, either under similar selection or through similar genomic mechanisms, we propose the use of genetic validity as a tool for identifying analogous pathology between animals and human neuropsychiatric disorders. Furthermore, the identification of copy number variants which disrupt entire genes, reinforces the notion that transgenic animals, such as knockouts or knock-ins, may provide unexpectedly valid disease models for psychiatric traits. To illustrate interspecies comparison of genetic variations in relation to neurobehavioral traits, examples are provided for the BDNF, COMT, and DISC1 genes in mouse and man. We propose that alignment of individual genetic variations with endophenotypes obtained from mice and across categories of neuropsychiatric disorders will provide an important step in translational research.
Collapse
|
65
|
Kõks S, Soomets U, Paya-Cano JL, Fernandes C, Luuk H, Plaas M, Terasmaa A, Tillmann V, Noormets K, Vasar E, Schalkwyk LC. Wfs1 gene deletion causes growth retardation in mice and interferes with the growth hormone pathway. Physiol Genomics 2009; 37:249-59. [PMID: 19293327 DOI: 10.1152/physiolgenomics.90407.2008] [Citation(s) in RCA: 36] [Impact Index Per Article: 2.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/22/2022] Open
Abstract
The aim of present study was to describe changes in gene expression in the temporal lobe of mice induced by deletion of the Wfs1 gene. Temporal lobes samples were analyzed using Affymetrix Mouse Genome 420 2 GeneChips and expression profiles were functionally annotated with GSEA and Ingenuity Pathway Analysis. We found that Wfs1 mutant mice are significantly smaller (20.9 +/- 1.6 g) than their wild-type counterparts (31.0 +/- 0.6 g, P < 0.0001). This difference existed in 129S6 and C57B6 backgrounds. Interestingly, microarray analysis identified upregulation of growth hormone (GH) transcripts and functional analysis revealed activation of GH pathways. In line with microarray data, the level of IGF-1 in the plasma of Wfs1 mutant mice was significantly increased (P < 0.05). Thus, Wfs1 deletion induces growth retardation, whereas the GH pathway is activated. To test the interaction between the Wfs1 deletion and genomic background, mutant mice were backcrossed to two different genetic backgrounds. In line with previous studies, an interaction between a gene knockout and genetic background was found in gene expression profiles in the congenic region. However, genetic background did not alter the effect of the Wfs1 mutation on either body weight or GH pathway activation. Further studies are needed to describe biochemical and molecular changes of the growth hormone axis as well as in other hormones to clarify their role in growth retardation in the Wfs1 mutant mice.
Collapse
|
66
|
Davis OSP, Plomin R, Schalkwyk LC. The SNPMaP package for R: a framework for genome-wide association using DNA pooling on microarrays. Bioinformatics 2008; 25:281-3. [PMID: 19008252 PMCID: PMC2639010 DOI: 10.1093/bioinformatics/btn587] [Citation(s) in RCA: 24] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/11/2023] Open
Abstract
Summary: Large-scale genome-wide association (GWA) studies using thousands of high-density SNP microarrays are becoming an essential tool in the search for loci related to heritable variation in many phenotypes. However, the cost of GWA remains beyond the reach of many researchers. Fortunately, the majority of statistical power can still be obtained by estimating allele frequencies from DNA pools, reducing the cost to that of tens, rather than thousands of arrays. We present a set of software tools for processing SNPMaP (SNP microarrays and pooling) data from CEL files to Relative Allele Scores in the rich R statistical computing environment. Availability: The SNPMaP package is available from http://cran.r-project.org/ under the GNU General Public License version 3 or later. Contact:snpmap@iop.kcl.ac.uk Supplementary information: Additional resources and test datasets are available at http://sgdp.iop.kcl.ac.uk/snpmap/
Collapse
|
67
|
Parsons MJ, Grimm CH, Paya-Cano JL, Sugden K, Nietfeld W, Lehrach H, Schalkwyk LC. Using hippocampal microRNA expression differences between mouse inbred strains to characterise miRNA function. Mamm Genome 2008; 19:552-60. [PMID: 18594913 PMCID: PMC2755775 DOI: 10.1007/s00335-008-9116-y] [Citation(s) in RCA: 32] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/06/2008] [Accepted: 05/17/2008] [Indexed: 11/13/2022]
Abstract
Micro-RNAs (miRNAs) are short, single-stranded, noncoding RNAs that are involved in the regulation of protein-coding genes at the level of messenger RNA (mRNA). They are involved in the regulation of numerous traits, including developmental timing, apoptosis, immune function, and neuronal development. To better understand how the expression of the miRNAs themselves is regulated, we looked for miRNA expression differences among four mouse inbred strains, A/J, BALB/cJ, C57BL/6J, and DBA/2J, in one tissue, the hippocampus. A total of 166 miRNA RT-PCR assays were used to screen RNA pools for each strain. Twenty miRNA species that were markedly different between strains were further investigated using eight individual samples per strain, and 11 miRNAs showed significant differences across strains (p < 0.05). This is the first observation of miRNA expression differences across inbred mice strains. We conducted an in silico correlation analysis of the expression of these differentially expressed miRNAs with phenotype data and mRNA expression to better characterise the effects of these miRNAs on both phenotype and the regulation of mRNA expression. This approach has allowed us to nominate miRNAs that have potential roles in anxiety, exploration, and learning and memory.
Collapse
|
68
|
Meaburn EL, Harlaar N, Craig IW, Schalkwyk LC, Plomin R. Quantitative trait locus association scan of early reading disability and ability using pooled DNA and 100K SNP microarrays in a sample of 5760 children. Mol Psychiatry 2008; 13:729-40. [PMID: 17684495 DOI: 10.1038/sj.mp.4002063] [Citation(s) in RCA: 88] [Impact Index Per Article: 5.5] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Submit a Manuscript] [Subscribe] [Scholar Register] [Indexed: 11/09/2022]
Abstract
Quantitative genetic research suggests that reading disability is the quantitative extreme of the same genetic and environmental factors responsible for normal variation in reading ability. This finding warrants a quantitative trait locus (QTL) strategy that compares low versus high extremes of the normal distribution of reading in the search for QTLs associated with variation throughout the distribution. A low reading ability group (N=755) and a high reading group (N=747) were selected from a representative UK sample of 7-year-olds assessed on two measures of reading that we have shown to be highly heritable and highly genetically correlated. The low and high reading ability groups were each divided into 10 independent DNA pools and the 20 pools were assayed on 100 K single nucleotide polymorphism (SNP) microarrays to screen for the largest allele frequency differences between the low and high reading ability groups. Seventy five of these nominated SNPs were individually genotyped in an independent sample of low (N=452) and high (N=452) reading ability children selected from a second sample of 4258 7-year-olds. Nine of the seventy-five SNPs were nominally significant (P<0.05) in the predicted direction. These 9 SNPs and 14 other SNPs showing low versus high allele frequency differences in the predicted direction were genotyped in the rest of the second sample to test the QTL hypothesis. Ten SNPs yielded nominally significant linear associations in the expected direction across the distribution of reading ability. However, none of these SNP associations accounted for more than 0.5% of the variance of reading ability, despite 99% power to detect them. We conclude that QTL effect sizes, even for highly heritable common disorders and quantitative traits such as early reading disability and ability, might be much smaller than previously considered.
Collapse
|
69
|
Caspi A, Williams B, Kim-Cohen J, Craig IW, Milne BJ, Poulton R, Schalkwyk LC, Taylor A, Werts H, Moffitt TE. Moderation of breastfeeding effects on the IQ by genetic variation in fatty acid metabolism. Proc Natl Acad Sci U S A 2007; 104:18860-5. [PMID: 17984066 PMCID: PMC2141867 DOI: 10.1073/pnas.0704292104] [Citation(s) in RCA: 266] [Impact Index Per Article: 15.6] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/14/2007] [Indexed: 11/18/2022] Open
Abstract
Children's intellectual development is influenced by both genetic inheritance and environmental experiences. Breastfeeding is one of the earliest such postnatal experiences. Breastfed children attain higher IQ scores than children not fed breast milk, presumably because of the fatty acids uniquely available in breast milk. Here we show that the association between breastfeeding and IQ is moderated by a genetic variant in FADS2, a gene involved in the genetic control of fatty acid pathways. We confirmed this gene-environment interaction in two birth cohorts, and we ruled out alternative explanations of the finding involving gene-exposure correlation, intrauterine growth, social class, and maternal cognitive ability, as well as maternal genotype effects on breastfeeding and breast milk. The finding shows that environmental exposures can be used to uncover novel candidate genes in complex phenotypes. It also shows that genes may work via the environment to shape the IQ, helping to close the nature versus nurture debate.
Collapse
|
70
|
Lad HV, Liu L, Payá-Cano JL, Fernandes C, Schalkwyk LC. Quantitative traits for the tail suspension test: automation, optimization, and BXD RI mapping. Mamm Genome 2007; 18:482-91. [PMID: 17619104 PMCID: PMC1998875 DOI: 10.1007/s00335-007-9029-1] [Citation(s) in RCA: 19] [Impact Index Per Article: 1.1] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/12/2007] [Accepted: 04/17/2007] [Indexed: 01/08/2023]
Abstract
Immobility in the tail suspension test (TST) is considered a model of despair in a stressful situation, and acute treatment with antidepressants reduces immobility. Inbred strains of mouse exhibit widely differing baseline levels of immobility in the TST and several quantitative trait loci (QTLs) have been nominated. The labor of manual scoring and various scoring criteria make obtaining robust data and comparisons across different laboratories problematic. Several studies have validated strain gauge and video analysis methods by comparison with manual scoring. We set out to find objective criteria for automated scoring parameters that maximize the biological information obtained, using a video tracking system on tapes of tail suspension tests of 24 lines of the BXD recombinant inbred panel and the progenitor strains C57BL/6J and DBA/2J. The maximum genetic effect size is captured using the highest time resolution and a low mobility threshold. Dissecting the trait further by comparing genetic association of multiple measures reveals good evidence for loci involved in immobility on chromosomes 4 and 15. These are best seen when using a high threshold for immobility, despite the overall better heritability at the lower threshold. A second trial of the test has greater duration of immobility and a completely different genetic profile. Frequency of mobility is also an independent phenotype, with a distal chromosome 1 locus.
Collapse
|
71
|
Docherty SJ, Butcher LM, Schalkwyk LC, Plomin R. Applicability of DNA pools on 500 K SNP microarrays for cost-effective initial screens in genomewide association studies. BMC Genomics 2007; 8:214. [PMID: 17610740 PMCID: PMC1925094 DOI: 10.1186/1471-2164-8-214] [Citation(s) in RCA: 38] [Impact Index Per Article: 2.2] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/31/2007] [Accepted: 07/04/2007] [Indexed: 01/02/2023] Open
Abstract
Background Genetic influences underpinning complex traits are thought to involve multiple quantitative trait loci (QTLs) of small effect size. Detection of such QTL associations requires systematic screening of large numbers of DNA markers within large sample populations. Using pooled DNA on SNP microarrays to screen for allelic frequency differences between groups such as cases and controls (called SNP Microarray and Pooling, or SNP-MaP) has been validated as an efficient solution on both 10 k and 100 k platforms. We demonstrate that this approach can be effectively applied to the truly genomewide Affymetrix GeneChip® Mapping 500 K Array. Results In comparisons between five independent DNA pools (N ~200 per pool) on separate Affymetrix GeneChip® Mapping 500 K Array sets, we show that, for SNPs with minor allele frequencies > 0.05, the reliability of the rank order of estimated allele frequencies, assessed as the average correlation between allele frequency estimates across the DNA pools, was 0.948 (average mean difference across the five pools = 0.069). Similarly, validity of the SNP-MaP approach was demonstrated by a rank-order correlation of 0.937 (average mean difference = 0.095) between the average DNA pool allele frequency estimates and the allele frequencies of an independent (CEPH) sample of 60 unrelated individually genotyped subjects. Conclusion We conclude that SNP-MaP can be extended for use on the Affymetrix GeneChip® Mapping 500 K Array, providing a cost-effective, reliable and valid initial screen of 500 K SNP microarrays in genomewide association scans.
Collapse
|
72
|
Kas MJH, Fernandes C, Schalkwyk LC, Collier DA. Genetics of behavioural domains across the neuropsychiatric spectrum; of mice and men. Mol Psychiatry 2007; 12:324-30. [PMID: 17389901 DOI: 10.1038/sj.mp.4001979] [Citation(s) in RCA: 90] [Impact Index Per Article: 5.3] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Submit a Manuscript] [Subscribe] [Scholar Register] [Indexed: 11/09/2022]
Abstract
Family and twin studies have revealed that genetic factors play a major role in psychiatric disorders, however, attempts to find susceptibility genes for these complex disorders have been largely unsuccessful. Therefore, new research strategies are required to tackle the complex interactions of genes, developmental, and environmental events. Here, we will address a behavioural domain concept that focuses on the genetics of behavioural domains relevant to both animal behaviour and across human psychiatric disorders. We believe that interspecies trait genetics rather than complex syndrome genetics will optimize genotype-phenotype relationships for psychiatric disorders and facilitate the identification of biological substrates underlying these disorders.
Collapse
|
73
|
Abstract
Microarrays are revolutionizing genetics by making it possible to genotype hundreds of thousands of DNA markers and to assess the expression (RNA transcripts) of all of the genes in the genome. Microarrays are slides the size of a postage stamp that contain millions of DNA sequences to which single-stranded DNA or RNA can hybridize. This miniaturization requires little DNA or RNA and makes the method fast and inexpensive; multiple assays of each target make the method highly accurate. DNA microarrays with hundreds of thousands of DNA markers have made it possible to conduct systematic scans of the entire genome to identify genetic associations with complex disorders or dimensions likely to be influenced by many genes of small effect size. RNA microarrays can provide snapshots of gene expression across all of the genes in the genome at any time in any tissue, which has far-reaching applications such as structural and functional 'genetic neuroimaging' and providing a biological basis for understanding environmental influence.
Collapse
|
74
|
Schalkwyk LC, Fernandes C, Nash MW, Kurrikoff K, Vasar E, Kõks S. Interpretation of knockout experiments: the congenic footprint. GENES BRAIN AND BEHAVIOR 2007; 6:299-303. [PMID: 17331107 PMCID: PMC1890814 DOI: 10.1111/j.1601-183x.2007.00304.x] [Citation(s) in RCA: 40] [Impact Index Per Article: 2.4] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Subscribe] [Scholar Register] [Indexed: 12/04/2022]
Abstract
In gene targeting experiments, the importance of genetic background is now widely appreciated, and knockout alleles are routinely backcrossed onto a standard inbred background. This produces a congenic strain with a substantial segment of embryonic stem (ES)-cell-derived chromosome still flanking the knockout allele, a phenomenon often neglected in knockout studies. In cholecystokynin 2 (Cckbr) knockout mice backcrossed with C57BL/6, we have found a clear ‘congenic footprint’ of expression differences in at least 10 genes across 40 Mb sequence flanking the Cckbr locus, each of which is potentially responsible for aspects of the ‘knockout’ phenotype. The expression differences are overwhelmingly in the knockout-low direction, which may point to a general phenomenon of background dependence. This finding emphasizes the need for caution in using gene knockouts to attribute phenotypic effects to genes. This is especially the case when the gene is of unknown function or the phenotype is unexpected, and is a particular concern for large-scale knockout and phenotypic screening programmes. However, the impact of genetic background should not be simply viewed as a potential confound, but as a unique opportunity to study the broader responses of a system to a specific (genetic) perturbation.
Collapse
|
75
|
Fernandes C, Hoyle E, Dempster E, Schalkwyk LC, Collier DA. Performance deficit of α7 nicotinic receptor knockout mice in a delayed matching-to-place task suggests a mild impairment of working/episodic-like memory. GENES BRAIN AND BEHAVIOR 2006; 5:433-40. [PMID: 16923147 DOI: 10.1111/j.1601-183x.2005.00176.x] [Citation(s) in RCA: 74] [Impact Index Per Article: 4.1] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 10/25/2022]
Abstract
Patients with schizophrenia exhibit deficits in a range of cognitive functions, particularly working and episodic memory, which are thought to be core features of the disorder. Memory dysfunction in schizophrenia is familial and thus a promising endophenotype for genetic studies. Both human and animal studies suggest a role for the neural nicotinic acid receptor family in cognition and specifically the alpha7-receptor subunit in schizophrenia and its endophenotypes. Consequently, we tested mice lacking the alpha7 subunit of the neural nicotinic receptor (B6.129S7-Chrna7(tm1Bay)/J) in the delayed matching-to-place (DMP) task of the Morris water maze, a measure of working/episodic memory akin to human episodic memory. We report that a minor impairment in alpha7 knockout mice was observed in the DMP task, with knockout mice taking longer to find the hidden platform than their wildtype controls. This suggests a role for the alpha7 subunit in working/episodic memory and a potential role for the alpha7 neural nicotinic receptor gene (CHRNA7) in schizophrenia and its endophenotypes.
Collapse
|