Effects of advanced paternal age on trajectories of social behavior in offspring.

Our study is the first investigation of the effects of advanced paternal age (APA) on the developmental trajectory of social behavior in rodent offspring. Given the strong epidemiological association between APA and sexually dimorphic neurodevelopmental disorders that are characterized by abnormalities in social behavior (autism, schizophrenia), we assessed sociability in male and female inbred mice (C57BL/6J) across postnatal development (N = 104) in relation to paternal age. We found differences in early social behavior in both male and female offspring of older breeders, with differences in this social domain persisting into adulthood in males only. We showed that these social deficits were not present in the fathers of these offspring, confirming a de novo origin of an altered social trajectory in the offspring generation. Our results, highly novel in rodent research, support the epidemiological observations in humans and provide evidence for a causal link between APA, age-related changes in the paternal sperm DNA and neurodevelopmental disorders in their offspring.

Identification of genes and gene pathways associated with major depressive disorder by integrative brain analysis of rat and human prefrontal cortex transcriptomes.

Despite moderate heritability estimates, progress in uncovering the molecular substrate underpinning major depressive disorder (MDD) has been slow. In this study, we used prefrontal cortex (PFC) gene expression from a genetic rat model of MDD to inform probe set prioritization in PFC in a human post-mortem study to uncover genes and gene pathways associated with MDD. Gene expression differences between Flinders sensitive (FSL) and Flinders resistant (FRL) rat lines were statistically evaluated using the RankProd, non-parametric algorithm. Top ranking probe sets in the rat study were subsequently used to prioritize orthologous selection in a human PFC in a case-control post-mortem study on MDD from the Stanley Brain Consortium. Candidate genes in the human post-mortem study were then tested against a matched control sample using the RankProd method. A total of 1767 probe sets were differentially expressed in the PFC between FSL and FRL rat lines at (q⩽0.001). A total of 898 orthologous probe sets was found on Affymetrix's HG-U95A chip used in the human study. Correcting for the number of multiple, non-independent tests, 20 probe sets were found to be significantly dysregulated between human cases and controls at q⩽0.05. These probe sets tagged the expression profile of 18 human genes (11 upregulated and seven downregulated). Using an integrative rat-human study, a number of convergent genes that may have a role in pathogenesis of MDD were uncovered. Eighty percent of these genes were functionally associated with a key stress response signalling cascade, involving NF-κB (nuclear factor kappa-light-chain-enhancer of activated B cells), AP-1 (activator protein 1) and ERK/MAPK, which has been systematically associated with MDD, neuroplasticity and neurogenesis.

Differential methylation of the TRPA1 promoter in pain sensitivity.

Chronic pain is a global public health problem, but the underlying molecular mechanisms are not fully understood. Here we examine genome-wide DNA methylation, first in 50 identical twins discordant for heat pain sensitivity and then in 50 further unrelated individuals. Whole-blood DNA methylation was characterized at 5.2 million loci by MeDIP sequencing and assessed longitudinally to identify differentially methylated regions associated with high or low pain sensitivity (pain DMRs). Nine meta-analysis pain DMRs show robust evidence for association (false discovery rate 5%) with the strongest signal in the pain gene TRPA1 (P=1.2 × 10(-13)). Several pain DMRs show longitudinal stability consistent with susceptibility effects, have similar methylation levels in the brain and altered expression in the skin. Our approach identifies epigenetic changes in both novel and established candidate genes that provide molecular insights into pain and may generalize to other complex traits.

Methylomic analysis of monozygotic twins discordant for autism spectrum disorder and related behavioural traits.

Autism spectrum disorder (ASD) defines a group of common, complex neurodevelopmental disorders. Although the aetiology of ASD has a strong genetic component, there is considerable monozygotic (MZ) twin discordance indicating a role for non-genetic factors. Because MZ twins share an identical DNA sequence, disease-discordant MZ twin pairs provide an ideal model for examining the contribution of environmentally driven epigenetic factors in disease. We performed a genome-wide analysis of DNA methylation in a sample of 50 MZ twin pairs (100 individuals) sampled from a representative population cohort that included twins discordant and concordant for ASD, ASD-associated traits and no autistic phenotype. Within-twin and between-group analyses identified numerous differentially methylated regions associated with ASD. In addition, we report significant correlations between DNA methylation and quantitatively measured autistic trait scores across our sample cohort. This study represents the first systematic epigenomic analyses of MZ twins discordant for ASD and implicates a role for altered DNA methylation in autism.

DNA methylation in interleukin-11 predicts clinical response to antidepressants in GENDEP.

Transcriptional differences in interleukin-11 (IL11) after antidepressant treatment have been found to correspond to clinical response in major depressive disorder (MDD) patients. Expression differences were partly mediated by a single-nucleotide polymorphism (rs1126757), identified as a predictor of antidepressant response as part of a genome-wide association study. Here we attempt to identify whether DNA methylation, another baseline factor known to affect transcription factor binding, might also predict antidepressant response, using samples collected from the Genome-based Therapeutic Drugs for Depression project (GENDEP). DNA samples from 113 MDD individuals from the GENDEP project, who were treated with either escitalopram (n=80) or nortriptyline (n=33) for 12 weeks, were randomly selected. Percentage change in Montgomery-Åsberg Depression Rating Scale scores between baseline and week 12 were utilized as our measure of antidepressant response. The Sequenom EpiTYPER platform was used to assess DNA methylation across the only CpG island located in the IL11 gene. Regression analyses were then used to explore the relationship between CpG unit methylation and antidepressant response. We identified a CpG unit predictor of general antidepressant response, a drug by CpG unit interaction predictor of response, and a CpG unit by rs1126757 interaction predictor of antidepressant response. The current study is the first to investigate the potential utility of pharmaco-epigenetic biomarkers for the prediction of antidepressant response. Our results suggest that DNA methylation in IL11 might be useful in identifying those patients likely to respond to antidepressants, and if so, the best drug suited to each individual.

Maternal separation is associated with strain-specific responses to stress and epigenetic alterations to Nr3c1, Avp, and Nr4a1 in mouse.

Stressful events early in life have been widely linked to behavioral phenotypes and have been implicated in the development of psychiatric disorders. Using a maternal separation paradigm, we investigated phenotypic and epigenetic changes following early life stress in two inbred strains of mice, C57BL/6J and DBA/2J. We found an increase in the corticosterone response to stress in male, C57BL/6J mice that had undergone maternal separation compared to controls. In addition, early life stress induced a number of mild but significant behavioral changes, many of which were sex and strain dependent. Following maternal separation anxiety was decreased in males but increased in DBA/2J females, DBA/2J males displayed reduced exploration of a novel object, and baseline activity was altered in males of both strains. Finally, we examined DNA methylation levels in the hippocampus across promoter regions of Nr3c1, Avp, and Nr4a1, and found altered levels at several CpG sites in maternally separated male mice compared to controls. This study contributes to a growing body of recent literature suggesting that epigenetic changes may mediate the impact of early life stress on behavior. In particular, we establish that the phenotypic and epigenetic responses to an adverse environment differ as a function of genetic background.

A B2 SINE insertion in the Comt1 gene (Comt1(B2i)) results in an overexpressing, behavior modifying allele present in classical inbred mouse strains.

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 (Comt1(B2i)) 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. Comt1(B2i) 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 Comt1(B2i) 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 Comt1(B2i) (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.

The differential transcriptome and ontology profiles of floating and cumulus granulosa cells in stimulated human antral follicles.

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.

A genome-wide association study identifies multiple loci associated with mathematics ability and disability.

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.

Wfs1 gene deletion causes growth retardation in mice and interferes with the growth hormone pathway.

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.

Quantitative trait locus association scan of early reading disability and ability using pooled DNA and 100K SNP microarrays in a sample of 5760 children.

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.

Interpretation of knockout experiments: the congenic footprint.

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.

Genetics of behavioural domains across the neuropsychiatric spectrum; of mice and men.

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.

Performance deficit of alpha7 nicotinic receptor knockout mice in a delayed matching-to-place task suggests a mild impairment of working/episodic-like memory.

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.

Association analysis of mild mental impairment using DNA pooling to screen 432 brain-expressed single-nucleotide polymorphisms.

We hypothesize that mild mental impairment (MMI) represents the low extreme of the same quantitative trait loci (QTLs) that operate throughout the distribution of intelligence. To detect QTLs of small effect size, we employed a direct association strategy by genotyping 432 presumably functional nonsynonymous single-nucleotide polymorphisms (nsSNPs) identified from public databases on DNA pools of 288 cases and 1025 controls. In total, 288 MMI cases were identified by in-home administration of McCarthy Scales of Children's Abilities to 836 twin pairs selected from a community sample of more than 14 000 children previously screened for nonverbal cognitive delay using parentally administered tests. Controls were selected from the community sample representing the full range of nonverbal intelligence. SNPs showing at least 7% allele frequency differences between case and control DNA pools were tested for their association with the full range of nonverbal intelligence using five DNA subpools, each representing quintiles of the normal quantitative trait scores from the 1025 controls. SNPs showing linear associations in the expected direction across quintiles using pooled DNA were individually genotyped for the 288 cases and 1025 controls and analyzed using standard statistical methods. One SNP (rs1136141) in HSPA8 met these criteria, yielding a significant (P=0.036) allelic frequency difference between cases and controls for individual genotyping and a significant (P=0.013) correlation within the control group that accounts for 0.5% of the variance. The present SNP strategy combined with DNA pooling and large samples represents a step towards identifying QTLs of small effect size associated with complex traits in the postgenomic era when all functional polymorphisms will be known.

Home-cage activity in heterogeneous stock (HS) mice as a model of baseline activity.

Behavioral genetic work in humans indicates that clinical hyperactivity is best viewed as the extreme end of activity levels in the population. However, current animal models of hyperactivity are not studied as quantitative traits as they are either knockout models or inbred strains. Furthermore, these animal models generally demonstrate elevated locomotion in novel environments, but not in their home-cages. This is the opposite of the symptoms seen in the human condition where childhood hyperactivity is generally more pronounced in constant, unstimulating situations. In this study we filmed an outbred population of 44 heterogeneous stock (HS) mice under red light during their active phase, to assess the reliability of individual differences in home-cage behavior and extract an index of home-cage activity (HCA) level. We then compared this measure to locomotor behavior in a novel environment--the open-field. Reliable individual differences in home-cage behaviors such as running, swinging on bars, and burrowing were found, and principal component factor analysis yielded a general activity factor, which accounted for 32% of the variance and correlated 0.90 with a subjective impression of activity level. The correlation between HCA and locomotor activity in the open-field was 0.23, which was non-significant. However, the association with HCA level appeared to increase over the five minutes of the open-field, presumably as the mice habituated. Furthermore, although mice displaying particularly high and low HCA were indistinguishable early in the open-field task, they became significantly differentiated over time. We conclude that home-cage behaviors and the open-field, after habituation, display good face and construct validity, and may provide a good model of baseline activity for quantitative trait loci (QTL) discovery and functional genomics in the HS mice.

Advanced integrated mouse YAC map including BAC framework.

Functional characterization of the mouse genome requires the availability of a comprehensive physical map to obtain molecular access to chromosomal regions of interest. Positional cloning remains a crucial way of linking phenotype with particular genes. A key step and frequent stumbling block in positional cloning is making a contig of a genetically defined candidate region. The most efficient first step is isolating YAC (Yeast Artificial Chromosome) clones. A robust, detailed YAC contig map is thus an important tool. Employing Interspersed Repetitive Sequence (IRS)-PCR genomics, we have generated an advanced second-generation YAC contig map of the mouse genome that doubles both the depth of clones and the density of markers available. In addition to the primarily YAC-based map, we located 1942 BAC (Bacterial Artificial Chromosome) clones. This allows us to present for the first time a dense framework of BACs spanning the genome of the mouse, which, for instance, can serve as a nucleus for genomic sequencing. Four large-insert mouse YAC libraries from three different strains are included in our data, and our analysis incorporates the data of Hunter et al. and Nusbaum et al. There is a total of 20,205 markers on the final map, 12,033 from our own data, and a total of 56,093 YACs, of which 44,401 are positive for more than one marker.

High-throughput scanning of the rat genome using interspersed repetitive sequence-PCR markers.

We report the establishment of a hybridization-based marker system for the rat genome based on the PCR amplification of interspersed repetitive sequences (IRS). Overall, 351 IRS markers were mapped within the rat genome. The IRS marker panel consists of 210 nonpolymorphic and 141 polymorphic markers that were screened for presence/absence polymorphism patterns in 38 different rat strains and substrains that are commonly used in biomedical research. The IRS marker panel was demonstrated to be useful for rapid genome screening in experimental rat crosses and high-throughput characterization of large-insert genomic library clones. Information on corresponding YAC clones is made available for this IRS marker set distributed over the whole rat genome. The two existing rat radiation hybrid maps were integrated by placing the IRS markers in both maps. The genetic and physical mapping data presented provide substantial information for ongoing positional cloning projects in the rat.

ETn insertion in the mouse Adcy1 gene: transcriptional and phylogenetic analyses.

Early retrotransposons (ETn) are murine transposable elements, bearing some structural similarity to integrated proviruses, and can be insertional mutagens. We have recently identified the causative mutation of the barrelless (Adcy1brl) phenotype as an integration of a 5.7-kb ETn in an intron of the adenylyl cyclase type I (Adcy1) gene. In the present study, Northern blot analysis shows that the ETn insertion results in loss of the normal Adcy1 transcript, a finding consistent with the loss-of-function Adcy1brl mutation, and generation of shorter transcripts. These aberrant transcripts are the products of abnormal RNA splicing and termination owing to the inserted sequence, and transcription initiation within the 3' long terminal repeat (LTR) of the ETn. The DNA sequences of the LTRs were compared in phylogenetic analyses with LTRs from 22 other ETn-related sequences. Three distinct families of ETn sequences can be identified on the basis of their LTRs. The ETn found in Adcy1brl is a member of a family that includes all classified ETn elements known to have recently transposed. Further, of the four known solitary (solo) LTRs, we have identified two that show evidence of recombination between LTRs from different ETn families.