Software | Open | Published:
Development of a versatile enrichment analysis tool reveals associations between the maternal brain and mental health disorders, including autism
BMC Neurosciencevolume 14, Article number: 147 (2013)
A recent study of lateral septum (LS) suggested a large number of autism-related genes with altered expression in the postpartum state. However, formally testing the findings for enrichment of autism-associated genes proved to be problematic with existing software. Many gene-disease association databases have been curated which are not currently incorporated in popular, full-featured enrichment tools, and the use of custom gene lists in these programs can be difficult to perform and interpret. As a simple alternative, we have developed the Modular Single-set Enrichment Test (MSET), a minimal tool that enables one to easily evaluate expression data for enrichment of any conceivable gene list of interest.
The MSET approach was validated by testing several publicly available expression data sets for expected enrichment in areas of autism, attention deficit hyperactivity disorder (ADHD), and arthritis. Using nine independent, unique autism gene lists extracted from association databases and two recent publications, a striking consensus of enrichment was detected within gene expression changes in LS of postpartum mice. A network of 160 autism-related genes was identified, representing developmental processes such as synaptic plasticity, neuronal morphogenesis, and differentiation. Additionally, maternal LS displayed enrichment for genes associated with bipolar disorder, schizophrenia, ADHD, and depression.
The transition to motherhood includes the most fundamental social bonding event in mammals and features naturally occurring changes in sociability. Some individuals with autism, schizophrenia, or other mental health disorders exhibit impaired social traits. Genes involved in these deficits may also contribute to elevated sociability in the maternal brain. To date, this is the first study to show a significant, quantitative link between the maternal brain and mental health disorders using large scale gene expression data. Thus, the postpartum brain may provide a novel and promising platform for understanding the complex genetics of improved sociability that may have direct relevance for multiple psychiatric illnesses. This study also provides an important new tool that fills a critical analysis gap and makes evaluation of enrichment using any database of interest possible with an emphasis on ease of use and methodological transparency.
Large scale, genome wide expression studies (such as exon microarrays, promoter microarrays, ChIP-on-CHIPS, and next gen sequencing) have allowed researchers to move from traditional approaches that focus on analyzing one or a few genes to new methods that generate high volumes of simultaneously collected expression data for tens of thousands of genes. This represents a powerful opportunity to identify alterations in large scale biological processes relevant to the study being performed, but it has also created a highly complex data environment with a unique set of challenges. Over the last decade, new software for analyzing enrichment of functionally related gene groups in large scale gene expression data has been developed and implemented successfully. There is an abundance of programs currently available that perform enrichment analyses in a variety of ways, each with its own advantages, disadvantages, and suitable applications [1–6].
Our laboratory recently carried out microarray experiments to explore gene expression changes in the lateral septum (LS) of outbred mice associated with the transition from a virgin to lactating, maternal state. The LS is part of an interconnected network of brain regions known to be critically important in social and maternal behavior [7, 8]. Several enrichment analyses were performed, including the Broad Institute’s Gene Set Enrichment Analysis (GSEA) and the NIH’s Database for Annotation, Visualization and Integrated Discovery (DAVID) [1, 3, 9, 10]. These programs proved to be fruitful in profiling large scale alterations in functionally related networks of genes in the postpartum LS, and enrichment was discovered for pathways related to ion channel activity, developmental processes, cyclic nucleotide metabolism, nucleosome components, and the Ras superfamily of small GTPases .
We began to notice that numerous genes among the most significant expression results in maternal LS had strong links to autism. This was intriguing because both autism and maternity involve behavioral and emotional alterations. Autism is a pervasive class of spectrum disorders typically characterized by deficits in social interaction, impaired communication abilities, and patterns of repetitive behavior [12, 13]. Autistic individuals can have difficulty empathizing with others, developing spoken language, and maintaining relationships [14, 15]. In comparison, the transition from a virgin to postpartum state represents a profound transformation in which an animal that was previously concerned mostly with its own survival develops heightened aspects of sociability and becomes critically focused on the wellbeing of its offspring . Because the establishment of the maternal phenotype is a natural process, the postpartum, outbred mouse provides a useful model for exploring the genetics of social behavior and its dysregulation . Oxytocin and vasopressin signaling in rodent maternal care has recently been investigated for its translational relevance to autism spectrum disorders in humans [18–22]. Moving from a single molecule focus to large scale genetic comparisons between the maternal brain and autism may provide broader insight to the neurological basis of sociability and contribute to our understanding of mental health disorders.
We recognized that it would be valuable to formally test for enrichment of autism-related genes in the maternal LS, but we quickly encountered several difficulties. Most available enrichment analysis tools operate solely on functional and pathway association databases, such as Gene Ontology (GO) terms or Kyoto Encyclopedia of Genes and Genomes (KEGG) pathways, and do not include curated groups of disease-linked genes . The handful of programs that do support disease ontology make use of only few predominant disease databases, such as the Online Mendelian Inheritance in Man (OMIM) and the Genetic Association Database (GAD) [23–25]. These databases are detailed and comprehensive, but because enrichment analysis is critically sensitive to the gene lists used as input, it is detrimental to be restricted to a single source of disease associations. There is a vibrant landscape of researchers and organizations that build and maintain numerous genetic disease association databases, but these resources are not linked to any formal means of enrichment analysis (Additional file 1: Table S1) [26–29]. In addition to comprehensive collections of gene-disease associations that encompass many diseases, it is not uncommon for a high profile disorder to engender a community of foundations, initiatives, and projects dedicated to furthering genetic research for that disease alone. For autism in particular, such notable groups include AutismKB by Peking University’s Center for Bioinformatics , AutDB by the nonprofit organization Mindspec , and the Autism Genetic Database (AGD) . Furthermore, ongoing research continuously generates new sets of candidate autism genes through a variety of novel techniques. For example, one recent publication used a novel Genome-wide Association Study method that reduces statistical noise and minimizes false positives to build a list of autism candidate genes , while another group used a functional profile of 84 autism-linked genes to screen the human genome for a larger set of predictive autism susceptibility genes .
To make use of this deep and varied field of autism genetics in enrichment analysis, it was clear that researchers must be allowed to create and use custom gene lists in a statistically robust hypothesis testing method. While some popular enrichment packages permit the inclusion of user-generated gene lists, the practicality of this feature often suffers from overly complex, opaque computations and internal significance corrections, a loss of large scale expression data fidelity through forced species-to-species conversions, and a heavy reliance on gene ranking methods and myriad configuration options that can easily lead to vastly different results that may be spurious or open to misinterpretation. As a simple alternative that focuses on the use of custom gene lists, we developed the Modular Single-set Enrichment Test (MSET) – a minimal randomization script that allows researchers to create and utilize any list of disease-associated genes as modules for enrichment testing and gene identification. MSET is designed to be accessible and easy for researchers of any background to use; it runs in R (free to download), requires only a single command line to execute, and subsequently interacts with the user via a simple window and text prompt interface. The lightweight architecture of the program confers independence from platform specific or species specific annotation systems and introduces minimal additional handling of expression data. Such an approach trades the exploratory aspect of full-featured enrichment software that searches for and scores enrichment across numerous biological pathways for a more focused hypothesis testing analysis for enrichment in diseases of interest. This allows investigators to use the most current database for a given disorder, customize their own, or repeat the analysis using numerous databases to generate a comprehensive meta-analysis and obtain high confidence in their results. Another important aspect of MSET is that it quickly highlights genes of interest that can be examined in follow up studies. As a result of its versatility, modular power, and ease of use, MSET represents an improvement over existing tools for performing enrichment analysis with independently assembled disease-associated gene sets.
Using the MSET approach, we discovered a striking consensus of enrichment for expression changes of autism-related genes across nine independent databases in the maternal LS compared to virgin. We found additional evidence of enrichment in the maternal LS for several other mental health diseases, including bipolar disorder, depression, and schizophrenia. This article will serve to present the biological findings from our data, as well as to introduce and demonstrate the novel enrichment tool used to generate them.
MSET is a script written in R that calls upon text files saved within the same folder as the script itself. Outlined simply, the MSET script is a randomization test that calculates the probability of randomly generating a set of microarray results that includes as many disease-associated genes as does a set of results from an actual experiment. It requires two input files: 1) a gene list of interest from a disease association database (autism-associated genes, schizophrenia-associated genes, etc.) for which enrichment will be evaluated, and 2) a full list of all microarray gene results ordered by significance. After choosing a threshold for how many of the top genes will be examined for enrichment, MSET saves them as a list of significant expression results. This list is separate from the full microarray results, which is also known as the microarray background. MSET then generates a specified number of simulated results by sampling randomly and without replacement from the microarray background. It counts the number of genes in each set of simulated results that appear in the gene list of interest (matches to database), and the p-value is calculated as the proportion of simulated results which contained at least as many matches to database as the actual significant expression results.
MSET analysis is performed at the gene level, but microarray platforms can include probes that target intergenic areas and lack a gene ID annotation. Prior to permutation testing, these non-annotated reads are removed from the significant microarray results. Because some genes in microarrays may be assigned multiple probe sets, duplicate gene IDs are also removed from the top selected significant results list, such that each gene is represented only once. The same treatments are applied to the gene list of interest database to ensure completeness and uniqueness. For the microarray background, targets without annotations to genes are removed, but duplicate targets for a single gene are retained. This is part of a strategy to preserve true-to-platform probabilities when sampling targets from the background to build simulated results.
After the user specifies the number of simulated results to be generated from the microarray background (e.g., 10,000), the script begins sampling without replacement. For each simulated set of results, MSET actually builds a list that is twice the length of the true significant results, removes duplicate genes from it, and then discards the excess genes such that the simulated results list is then the same length as the actual results. In this way, genes that are represented by many probes have a greater chance to be included in the simulated results, corresponding to their “weight” in the microarray platform. This allows for simulated results that are true to the particular microarray platform.
MSET output consists of two elements; a graph of the probability distribution of matches to database in simulated results, and a text readout with a variety of information about the input gene lists, the simulated results, and enrichment p-value. It also displays the individual disease-associated genes from the database of interest that appeared in the actual microarray results.
All expression data used in this study were taken from published microarray studies that have been uploaded to NCBI’s Gene Expression Omnibus (GEO, http://0-www.ncbi.nlm.nih.gov.brum.beds.ac.uk/geo/) or an institutional website. Details of experimental design are provided in the original publications. Postpartum LS expression data were prepared in our laboratory with the Probe Logarithmic Intensity Error (PLIER) algorithm. For data accessed via the GEO, gene expression changes between experimental groups were ranked by the built in GEO2R tool. While a researcher can decide which probe set algorithm is optimal for their particular microarray experiment, differences in the current study are small enough that both methods yield approximately similarly ordered results.
The data sets supporting the results of this article are available via NCBI’s Gene Expression Omnibus (GEO) with the following links and accession numbers.
Nine independent online databases were selected as sources for autism-associated genes (Additional file 1: Table S1). Of these, four are general disease association databases, three are autism-specific aggregations, and two are recent publications that generated candidate autism gene lists. The autism gene lists extracted from the selected databases exhibited varying sizes and degrees of similarity as assessed by pairwise comparisons, but are generally non-redundant. Figure 1 illustrates the relatedness of the autism databases and shows how they fit into the modular workflow of the enrichment analysis.
A set of 809 genes displayed significant changes in expression in the postpartum LS relative to virgin (FDR-adjusted p < 0.25) as measured by microarray analysis. Within these results, enrichment for all nine autism gene lists was detected (Figure 1, p < 0.05). Collectively, 160 autism-associated genes were identified in the significant maternal LS expression results (Additional file 2: Table S2). To highlight the genes in this group that have the most evidence for an association to autism, 36 differentially expressed genes in the postpartum LS which appeared in three or more autism databases are shown in Table 1. A functional profile of all 160 autism-linked genes in the postpartum LS results generated by NIH’s DAVID functional annotation clustering tool revealed striking themes of differentiation/development, synaptic plasticity, and neuronal signaling (Table 2). Using gene lists extracted from the general disease association databases, enrichment was additionally discovered in the maternal LS for genes related to bipolar disorder, schizophrenia, ADHD, and depression. This is shown in Figure 2, which is a comprehensive heat map that summarizes enrichment results for all data sets analyzed in this study.
As a demonstration of MSET’s applicability, enrichment of autism gene lists was evaluated in a set of microarray results from the developing neocortex of T-box brain gene 1 (Tbr1) null mice, a putative model for autism genetics. In the 809 most significant results, these expression data exhibited compelling enrichment in seven out of the nine autism gene lists (P < 0.05), with subtler enrichment in two gene lists that bordered on or failed to reach significance (Figure 3). In contrast to the maternal LS results, the Tbr1 data showed little to no enrichment in genes associated with other mental health disorders (Figure 2).
Furthermore, enrichment analysis for the full complement of mental health disorders was conducted within a set of 809 significantly altered genes in microarray results from mice treated with methylphenidate, a common treatment for ADHD. These expression data exhibited some enrichment for autism gene lists (two out of nine, p < 0.05), but notably showed high enrichment specifically for ADHD gene lists (Figure 2).
As a test for specificity of MSET analysis, enrichment for autism and arthritis-associated gene lists was compared within two sets of expression data; the 809 significant maternal LS microarray results and 809 microarray results from a serum induced murine arthritis model. Enrichment of arthritis-associated genes was detected in the arthritis model expression data, but not in the maternal LS results (p < 0.05). Additionally, enrichment of autism gene lists was not detected within the arthritis model results (Figure 4). The arthritis model results failed to show enrichment for any of the mental health disorders (Figure 2).
In this study, we have introduced the Modular Single-set Enrichment Test (MSET), a newly developed tool designed for assessing enrichment of disease/disorder-associated gene sets within microarray results. In addition to demonstrating the capabilities and limitations of this novel software, we used it to discover a strong link between the maternal brain and autism, as well as several other mental health disorders. We subsequently identified a network of candidate genes that may influence sociability in mothers and revealed the functional character of this network to be primarily related to developmental and neuronal signaling processes.
Enrichment in postpartum LS for autism, and other mental health disorders
The compelling enrichment of autism-associated genes found in expression changes of the postpartum LS (Figure 1) is, to our knowledge, the first demonstration of a genetic link between the maternal brain and pathways involved in autism. The discovery of enrichment in postpartum LS for several mood/social disorders (Figure 2) suggests that the phenotypic consequence of LS gene changes in the transition to motherhood possesses a significant behavioral and emotional component. Because the mother-infant relationship is the first and foremost social bond formed in mammals, it has been suggested that the genetic and neural networks underlying sociability in this ancestral event might serve as an evolutionary template from which sociability in other contexts is derived . Our data indirectly support this concept on a large-scale genetic level. While autism rates are higher in males , it could be the conserved use of the same core genes for sociability that provides the connection between autism and the maternal brain. Furthermore, the severity of autism symptoms is often described as spectral, rather than binary. It is therefore plausible that subtle dysregulation of genes which are naturally modulated in the control of sociability, such as in the transition to motherhood, would more likely contribute to this observed phenomenon than would rarer gain-of-function or loss-of-function mutations.
Table 1 presents 36 of the 160 autism-associated genes that MSET identified in the postpartum LS, as well as the number of autism databases in which they are featured. This is presumably a reflection of the strength of their association with autism based on past studies, with consensus genes having the most widely recognized evidence. However, it is not a perfect indicator because it only counts positive association discoveries, and does not consider the existence of any potential contradictory evidence or disagreement. For example, Foxp2, a forkhead/winged helix (FOX) transcription factor, is found in seven of the nine autism databases used (Table 1). It is located in a region of chromosome 7q that has been linked to autism in the past , and mutations in Foxp2 cause speech and language acquisition pathologies in humans . However, more recent evidence suggests that the language deficits are more directly related to a developmental impairment of motor brain regions, rather than to social behavior, and several recent reports conclude that Foxp2 does not contribute to autism susceptibility [39–41]. Even if Foxp2 were omitted from significant postpartum LS expression results, the observed enrichment would be highly significant. This illustrates the important point that, although assessing the degree of enrichment using MSET is robust and largely resistant to single gene false positives in upstream databases, caution must be exercised when interpreting the biological importance of individual genes identified by MSET in the testing procedure. Another advantage of MSET is that the user can manually annotate any file, remove genes that are considered to be inappropriate, or even create novel gene lists for testing.
Table 1 includes several autism-linked genes that were identified in our original microarray analysis as particularly interesting based on their biological function and relevance to emotional state and behavior. These include the GABAA receptor subunits α4 and δ, four potassium channel subunits (Kcnd2, Kcnd3, Kcnh7, and Kcnj4), dopamine receptors Drd1a and Drd2, the kappa opioid receptor Oprk1, fatty acid binding protein 7 (Fabp7), and suppressor of cytokine signaling 2 (Socs2). The biology of these genes is discussed in greater detail in our original report .
NIH DAVID’s functional annotation clustering tool was used to generate a functional profile of the 160 autism-associated genes found to be differentially expressed in the postpartum LS (Table 2). The most highly enriched pathways were primarily developmental, involving processes such as synaptic plasticity, neuronal morphogenesis/differentiation, and cell motility. Several clusters related to synaptic transmission also showed high levels of enrichment. Because these biological processes have now been implicated in both autism and the maternal LS, it is likely that aspects of sociability modulated in both phenomena are influenced by structural changes in the brain, including axonal/dendritic growth, and even neurogenesis. This possibility is supported by a body of literature which has revealed that diverse regions of the adult brain contain multipotent stem cells capable of generating new neurons [42–49], and it has been shown that maternal behavior is associated with the stimulation of neurogenesis in the subventricular zone .
In addition to autism, MSET analysis revealed that significant postpartum LS expression results exhibit enrichment for bipolar disorder (BPD), schizophrenia, ADHD, and depression-associated genes (Figure 2). These gene lists were extracted from the four general disease association databases that were also used in the autism enrichment analysis (Additional file 1: Table S1). Links were particularly strong for both BPD and schizophrenia. BPD and depression links are of interest because rates of depression increase in the postpartum state, with postpartum depression affecting 1-10% of mothers . Thus, some of the normal changes that occur in the maternal brain likely lead to a vulnerability of key depression type pathways. Positive associations have been consistently found for an elevated risk of BPD in women after childbirth , which is considered to be part of a suite of diagnosable “postpartum psychoses”. Also among this class of diseases is schizophrenia, which, in addition to its well-known cognitive dysfunction, is also characterized by emotional deficits . Recent studies highlight that a subset of genes contribute to multiple mental health disorders [54–57], so it is not completely surprising that a behavioral transformation as fundamental as the transition to motherhood might have links to multiple disorders. To ensure that this multitude of positive enrichment was not due to an artifact in MSET analysis, we tested the postpartum LS expression results for enrichment of arthritis-associated genes (Figure 4), which proved to be absent. The MSET tool has been used successfully in our laboratory to detect enrichment of mental health-related gene sets in other areas within the maternal brain, such as the medial preoptic area (unpublished observations). While there were similarities in enrichment across regions, there were also differences in enrichment patterns and in the individual genes which accounted for enrichment. This indicates that there may be common, global expression changes in the maternal brain, but also that each region has its own genetic “signature”. Future work will characterize the genetic profile of the maternal brain more comprehensively.
Enrichment analysis in expression data from a murine model of induced arthritis
To validate and demonstrate the applicability of MSET, we performed a series of analyses on expression data taken from several independently conducted microarray experiments. These data are publicly available through NCBI’s Gene Expression Omnibus (http://0-www.ncbi.nlm.nih.gov.brum.beds.ac.uk/geo/) or through institutional hosting.
To test for expected specificity of MSET analysis, we assessed enrichment for the full range of disease-associated gene sets within microarray expression data from a murine arthritis model study in which an arthritic state was induced via the transfer of serum from a knockout mouse into a wild type animal . Complementary findings to the postpartum LS results were observed, in which the arthritis model data showed enrichment specifically for arthritis-associated gene sets, but not for autism (Figure 4) or any other mental health disorders (Figure 2). Collectively, these results demonstrate that the enrichment analysis performed in the present study is reliable and specific. Specificity may not be expected in every application, as different models and experimental treatments used in microarray studies can affect broad or unanticipated gene pathways.
Enrichment of autism-associated genes in Tbr1 null transgenic mice
To showcase the broader applicability of MSET, we performed enrichment analysis for the full range of disease-associated gene lists in a set of expression data collected from murine T-box brain gene 1 (Tbr1) null developing neocortex . The Tbr1 null animal was chosen because Tbr1 is a developmentally related transcription factor that binds, among other targets, the promoter of a gene called autism susceptibility candidate 2 (Auts2), named for its implication in autism susceptibility in the frontal cortex [59, 60]. The Tbr1 null neocortex was observed to be enriched specifically for autism-associated gene sets (Figure 3), and not for any other mental health disorders included in our analysis (Figure 2). These findings suggest that, although inviable shortly after birth, the Tbr1 knockout animal may provide a valuable model for the study of autism-related biology. The Tbr1 null expression data also showed enrichment in two out of the four arthritis-associated gene sets. While this is not particularly strong enrichment, the observed variability could be due to broader physiological changes across numerous systems (possibly including the immune response) that must undoubtedly be affected by the fatal null mutation.
MSET enrichment analysis in expression data from methylphenidate treated mice
In addition to using MSET to analyze enrichment in expression results from animals that have undergone a natural change (mothers) and transgenic animals (Tbr1 null), we also tested its capabilities in a set of expression data from mice that were subjected to a pharmacological treatment. In the study, mice were treated with chronic (90 days) exposure to methylphenidate, commonly used to treat ADHD, and microarray analysis was performed on microdissected substantia nigra pars compacta (SNpc) . In our enrichment analysis of these data, we observed a subtle degree of enrichment for autism-related gene sets (in three out of nine lists), but found that a consensus of enrichment was only detected for ADHD-related gene lists, and not for any other mental health disorder or arthritis (Figure 2). This shows that MSET can be effectively utilized with sensitivity in microarray data collected from a variety of different experimental protocols and treatments, providing a promising new strategy for exploring the genetics underlying mental health disorders from numerous, complementary angles.
Considerations and limitations of MSET analysis
MSET allows for powerful research possibilities, but there are numerous considerations that must be made regarding its appropriate application and the input parameters used. MSET utilizes a fairly simple gene randomization testing procedure to determine if members of a disease-associated gene set are overrepresented within significant microarray results compared to what would be expected by chance. This is in contrast to programs like GSEA, in which the coincident distribution of gene set members is characterized within a ranked list of microarray results using a running-sum statistic and correlated to phenotype with individual sample expression values . Accordingly, MSET calls for only one simple input file of summarized microarray gene results (in addition to disease-associated gene sets of interest), and does not require expression values, chip annotations, or phenotype/trait files. Some web applications exist for performing overrepresentation analysis (such as GOHyperGAll in Bioconductor, InnateDB, and GenMAPP-CS in the GO-Elite program), but they include problematic gene ID conversions, species limitations, a strict dependence on GO terms and existing ontologies, and inflexibility in generating custom gene sets. MSET represents an advancement in versatility and ease of use over the existing landscape of tools for testing enrichment of independently curated disease-associated gene sets.
While MSET is theoretically capable of testing for enrichment of genes linked to functional pathways, using other, more full-featured programs for this purpose is recommended. It has been proposed that gene independence is a safe assumption for enrichment analysis . Others have countered that gene-gene interaction can inflate p-values and generate false positives in functional enrichment. Due to the relatively simple nature of the randomization algorithm, independence is assumed in MSET analysis for disease enrichment. This is a safe assumption because disease gene sets are heterogeneous groups curated by phenotypic associations, rather than functional relatedness. However, it would be more conservative to use other programs that account for potential gene-gene interaction in functional enrichment analyses. GSEA, for instance, preserves gene-gene interaction by permuting labels of whole samples, rather than at the individual gene level .
Because enrichment analysis is highly sensitive to the input gene lists used , this will be a focus for much of the discussion regarding the performance of MSET. It can be seen in Figure 1 that the nine autism databases used in our analysis vary in both their identity and size. They also differ in the methods used to produce candidate autism gene lists; therefore, one should be aware that some databases may be more robustly assembled than others, and confidence in MSET results relies critically on confidence in upstream database quality. MSET’s reliability is bolstered by its capacity to test enrichment for multiple gene lists associated with the same disease; this feature minimizes the effects of weak associations on enrichment testing significance.
There is generally a balance between specificity of enrichment and the accuracy of its detection, which is related to gene list length. This is particularly relevant when extracting gene lists from general disease association databases (such as the DISEASES database, GAD, HuGE Phenopedia, and Malacards), which compile positive associations broadly across many diseases. Smaller gene lists may be more specific to their associated disease, but the MSET suffers from a decrease in the accuracy of hypothesis testing as the average number of matches found in simulated results becomes small. This can be seen in the probability density curves for database matches generated with the AutismKB and Malacards autism gene lists in Figure 1. Their “spikey” appearance reflects the highly discrete nature of distributions with a very small range. Consequently, chance variation in the number of matches in the microarray results being analyzed, even by a single gene, represents a disproportionately large jump in p-value from peak to peak. The smoother distributions generated from larger databases provide a much greater resolution for hypothesis testing; however, larger gene lists may be less biologically specific to the associated disease, and extremely large gene lists can result in false positive enrichment results. While there is assuredly some “true” degree of genetic overlap underlying various diseases, there is probably an additional level of similarity across seemingly unrelated conditions introduced artificially through the methodology of association studies and their aggregation. For example, one might expect genes featured in centrally important signaling pathways to show positive associations with many diseases and experimental conditions in microarray studies, leading to false positive results for enrichment of extremely large gene lists. Specificity can be further complicated by the detailed nature of disease association labels in comprehensive databases. For instance, the DISEASES database has separate gene lists for arthritis, psoriatic arthritis, osteoarthritis, rheumatoid arthritis, septic arthritis, and more. The primary capacity of MSET to overcome these factors is rooted in its ability to be repeated modularly to generate a meta-analysis. This allows for isolated enrichment findings to be interpreted within the context of larger patterns. Also, because a deeper and more refined body of resources exists for autism genetics than for the other disorders featured in this study, we have relatively greater confidence in the downstream autism enrichment results. As ongoing research adds to our genetic understanding of various diseases, the MSET tool is in an ideal position to allow researchers to swiftly adapt and make use of updated knowledge bases in the future.
The tradeoff between specificity and accuracy also applies to the significant microarray results in which enrichment testing is performed. For postpartum LS expression data, we used an FDR-adjusted p-value of 0.25 as a significance threshold, which produced 809 genes from the microarray background. Researchers using other model organisms or biological systems may want to use different criteria for statistical significance. Other microarray studies may not yield a high number of significant gene changes by FDR-adjusted p-values. In these situations, a less stringent significance threshold may be applied to make use of a larger number of results, but the greater inclusivity and incidence of false discovery may render them somewhat less biologically meaningful. The subsequent enrichment analysis must therefore be taken with an accordingly critical interpretation.
MSET is designed to allow the user to conduct the most appropriate examination possible for enrichment of one or more disorders in a particular set of expression data. Due to the necessarily customizable nature of the input parameters that make for a quality assessment in one set of microarray results, it is difficult to objectively compare enrichment across numerous expression results. In the current study, we have done so by standardizing both the number of significant expression results selected from the background and the databases used. It cannot be assumed that the 809 most significant genes from one study are as meaningful or specific as those from another study, but the comprehensive and identically repeated analysis performed here is a valuable preliminary comparison. Collectively, the analyses undertaken in the current study provide a promising indication that the MSET method can be a valuable and informative approach to large scale genetic questions.
In this article, we present the novel enrichment tool MSET as a lightweight alternative for performing enrichment analysis with custom-built, disease-associated gene lists. The value of MSET is twofold; firstly in its ease of use for researchers of any background due to its methodological simplicity and transparency, and secondly in the statistical hypothesis testing power offered by the ability to test enrichment within expression data modularly with numerous, independent gene-association databases. Using the MSET approach, we discovered that the maternal LS is highly enriched for autism-associated genes, as well as for genes linked to related mental health disorders, thus supporting a role for the maternal brain as a valuable translational model to the genetics of social and emotional disorders. By matching genes within significant expression results to autism genetic databases, we constructed a network of candidate genes that may regulate sociability in new mothers as well as in contexts of psychiatric illness. We further revealed that this network is primarily related to developmental processes, indicating that structural changes in LS likely underlie the modulation of social interaction and emotional reactivity in the postpartum state.
Availability and requirements
The Modular Single-set Enrichment Test (MSET) is available for download as Additional file 3, and is hosted online at https://sourceforge.net/projects/mset2013. It is written in R and requires the R command console to run, which can be downloaded from http://www.r-project.org/. It is licensed under the open source Apache 2.0 license. Refer to the MSET manual contained within the MSET folder for instructions on installation and use.
Multiple single-set enrichment test
Attention deficit hyperactivity disorder
Gene set enrichment analysis
Database for annotation, visualization and integrated discovery
Kyoto encyclopedia of genes and genomes
Online mendelian inheritance in man
Genetic association database
Autism genetic database
Gene expression omnibus.
Da Huang W, Sherman BT, Lempicki RA: Bioinformatics enrichment tools: paths toward the comprehensive functional analysis of large gene lists. Nucleic Acids Res. 2009, 37 (1): 1-13. 10.1093/nar/gkn923.
Curtis RK, Oresic M, Vidal-Puig A: Pathways to the analysis of microarray data. Trends Biotechnol. 2005, 23 (8): 429-435. 10.1016/j.tibtech.2005.05.011.
Da Huang W, Sherman BT, Lempicki RA: Systematic and integrative analysis of large gene lists using DAVID bioinformatics resources. Nat Protoc. 2009, 4 (1): 44-57.
Ackermann M, Strimmer K: A general modular framework for gene set enrichment analysis. BMC Bioinformatics. 2009, 10: 47-10.1186/1471-2105-10-47.
Smyth GK: Linear models and empirical bayes methods for assessing differential expression in microarray experiments. Stat Appl Genet Mol Biol. 2004, 3: Article3.
Kim SY, Volsky DJ: PAGE: parametric analysis of gene set enrichment. BMC Bioinformatics. 2005, 6: 144-10.1186/1471-2105-6-144.
Sheehan TP, Chambers RA, Russell DS: Regulation of affect by the lateral septum: implications for neuropsychiatry. Brain Res Brain Res Rev. 2004, 46 (1): 71-117. 10.1016/j.brainresrev.2004.04.009.
Olazabal DE, Pereira M, Agrati D, Ferreira A, Fleming AS, Gonzalez-Mariscal G, Levy F, Lucion AB, Morrell JI, Numan M, et al: Flexibility and adaptation of the neural substrate that supports maternal behavior in mammals. Neurosci Biobehav Rev. 2013, 37 (8): 1857-1892.
Subramanian A, Tamayo P, Mootha VK, Mukherjee S, Ebert BL, Gillette MA, Paulovich A, Pomeroy SL, Golub TR, Lander ES, et al: Gene set enrichment analysis: a knowledge-based approach for interpreting genome-wide expression profiles. Proc Natl Acad Sci USA. 2005, 102 (43): 15545-15550. 10.1073/pnas.0506580102.
Mootha VK, Lindgren CM, Eriksson KF, Subramanian A, Sihag S, Lehar J, Puigserver P, Carlsson E, Ridderstrale M, Laurila E, et al: PGC-1alpha-responsive genes involved in oxidative phosphorylation are coordinately downregulated in human diabetes. Nat Genet. 2003, 34 (3): 267-273. 10.1038/ng1180.
Eisinger BE, Zhao C, Driessen TM, Saul MC, Gammie SC: Large scale expression changes of genes related to neuronal signaling and developmental processes found in lateral septum of postpartum outbred mice. PLoS One. 2013, 8 (5): e63824-10.1371/journal.pone.0063824.
Mittal VA, Walker EF: Diagnostic and statistical manual of mental disorders. Psychiatry Res. 2011, 189 (1): 158-159. 10.1016/j.psychres.2011.06.006.
Lord C, Risi S, Lambrecht L, Cook EH, Leventhal BL, DiLavore PC, Pickles A, Rutter M: The autism diagnostic observation schedule-generic: a standard measure of social and communication deficits associated with the spectrum of autism. J Autism Dev Disord. 2000, 30 (3): 205-223. 10.1023/A:1005592401947.
Noens IL, Van Berckelaer-Onnes IA: Captured by details: sense-making, language and communication in autism. J Commun Disord. 2005, 38 (2): 123-141. 10.1016/j.jcomdis.2004.06.002.
Landry SH, Loveland KA: The effect of social context on the functional communication skills of autistic children. J Autism Dev Disord. 1989, 19 (2): 283-299. 10.1007/BF02211847.
Kuroda KO, Tachikawa K, Yoshida S, Tsuneoka Y, Numan M: Neuromolecular basis of parental behavior in laboratory mice and rats: with special emphasis on technical issues of using mouse genetics. Prog Neuropsychopharmacol Biol Psychiatry. 2011, 35 (5): 1205-1231. 10.1016/j.pnpbp.2011.02.008.
Kumar RC: The maternal brain as a model for investigating mental illness. Prog Brain Res. 2001, 133: 333-338.
Insel TR: The challenge of translation in social neuroscience: a review of oxytocin, vasopressin, and affiliative behavior. Neuron. 2010, 65 (6): 768-779. 10.1016/j.neuron.2010.03.005.
Neumann ID: Brain mechanisms underlying emotional alterations in the peripartum period in rats. Depress Anxiety. 2003, 17 (3): 111-121. 10.1002/da.10070.
Veenema AH, Neumann ID: Central vasopressin and oxytocin release: regulation of complex social behaviours. Prog Brain Res. 2008, 170: 261-276.
Neumann ID: Brain oxytocin: a key regulator of emotional and social behaviours in both females and males. J Neuroendocrinol. 2008, 20 (6): 858-865. 10.1111/j.1365-2826.2008.01726.x.
Lukas M, Neumann ID: Oxytocin and vasopressin in rodent behaviors related to social dysfunctions in autism spectrum disorders. Behav Brain Res. 2013, 251: 85-94.
Becker KG, Barnes KC, Bright TJ, Wang SA: The genetic association database. Nat Genet. 2004, 36 (5): 431-432. 10.1038/ng0504-431.
Zhang Y, De S, Garner JR, Smith K, Wang SA, Becker KG: Systematic analysis, comparison, and integration of disease based human genetic association data and mouse genetic phenotypic information. BMC Med Genomics. 2010, 3: 1-10.1186/1755-8794-3-1.
Hamosh A, Scott AF, Amberger JS, Bocchini CA, McKusick VA: Online Mendelian Inheritance in Man (OMIM), a knowledgebase of human genes and genetic disorders. Nucleic Acids Res. 2005, 33: D514-517.
Yu W, Clyne M, Khoury MJ, Gwinn M: Phenopedia and Genopedia: disease-centered and gene-centered views of the evolving knowledge of human genetic associations. Bioinformatics. 2010, 26 (1): 145-146. 10.1093/bioinformatics/btp618.
HuGE Navigator. http://www.hugenavigator.net/HuGENavigator/home.do.
Sune Frankild LJJ: DISEASES: Disease-gene associations mined from literature. http://diseases.jensenlab.org/Search.
Xu LM, Li JR, Huang Y, Zhao M, Tang X, Wei L: AutismKB: an evidence-based knowledgebase of autism genetics. Nucleic Acids Res. 2012, 40: D1016-1022. 10.1093/nar/gkr1145.
Basu SN, Kollu R, Banerjee-Basu S: AutDB: a gene reference resource for autism research. Nucleic Acids Res. 2009, 37: D832-836. 10.1093/nar/gkn835.
Matuszek G, Talebizadeh Z: Autism Genetic Database (AGD): a comprehensive database including autism susceptibility gene-CNVs integrated with known noncoding RNAs and fragile sites. BMC Med Genet. 2009, 10: 102-10.1186/1471-2350-10-102.
Hussman JP, Chung RH, Griswold AJ, Jaworski JM, Salyakina D, Ma D, Konidari I, Whitehead PL, Vance JM, Martin ER, et al: A noise-reduction GWAS analysis implicates altered regulation of neurite outgrowth and guidance in autism. Molecular Autism. 2011, 2 (1): 1-10.1186/2040-2392-2-1.
Kumar A, Swanwick CC, Johnson N, Menashe I, Basu SN, Bales ME, Banerjee-Basu S: A brain region-specific predictive gene map for autism derived by profiling a reference gene set. PLoS One. 2011, 6 (12): e28431-10.1371/journal.pone.0028431.
Broad KD, Curley JP, Keverne EB: Mother-infant bonding and the evolution of mammalian social relationships. Philos Trans R Soc Lond B Biol Sci. 2006, 361 (1476): 2199-2214. 10.1098/rstb.2006.1940.
Fombonne E: The epidemiology of autism: a review. Psychol Med. 1999, 29 (4): 769-786. 10.1017/S0033291799008508.
Bonora E, Lamb JA, Barnby G, Sykes N, Moberly T, Beyer KS, Klauck SM, Poustka F, Bacchelli E, Blasi F, et al: Mutation screening and association analysis of six candidate genes for autism on chromosome 7q. European Journal of Human Genetics : EJHG. 2005, 13 (2): 198-207. 10.1038/sj.ejhg.5201315.
Lai CS, Fisher SE, Hurst JA, Vargha-Khadem F, Monaco AP: A forkhead-domain gene is mutated in a severe speech and language disorder. Nature. 2001, 413 (6855): 519-523. 10.1038/35097076.
Newbury DF, Bonora E, Lamb JA, Fisher SE, Lai CS, Baird G, Jannoun L, Slonims V, Stott CM, Merricks MJ, et al: FOXP2 is not a major susceptibility gene for autism or specific language impairment. Am J Hum Genet. 2002, 70 (5): 1318-1327. 10.1086/339931.
Marui T, Koishi S, Funatogawa I, Yamamoto K, Matsumoto H, Hashimoto O, Nanba E, Kato C, Ishijima M, Watanabe K, et al: No association of FOXP2 and PTPRZ1 on 7q31 with autism from the Japanese population. Neurosci Res. 2005, 53 (1): 91-94. 10.1016/j.neures.2005.05.003.
Wassink TH, Piven J, Vieland VJ, Pietila J, Goedken RJ, Folstein SE, Sheffield VC: Evaluation of FOXP2 as an autism susceptibility gene. Am J Med Genet. 2002, 114 (5): 566-569. 10.1002/ajmg.10415.
Reynolds BA, Weiss S: Generation of neurons and astrocytes from isolated cells of the adult mammalian central nervous system. Science. 1992, 255 (5052): 1707-1710. 10.1126/science.1553558.
Richards LJ, Kilpatrick TJ, Bartlett PF: De novo generation of neuronal cells from the adult mouse brain. Proc Natl Acad Sci USA. 1992, 89 (18): 8591-8595. 10.1073/pnas.89.18.8591.
Gage FH, Coates PW, Palmer TD, Kuhn HG, Fisher LJ, Suhonen JO, Peterson DA, Suhr ST, Ray J: Survival and differentiation of adult neuronal progenitor cells transplanted to the adult brain. Proc Natl Acad Sci USA. 1995, 92 (25): 11879-11883. 10.1073/pnas.92.25.11879.
Palmer TD, Takahashi J, Gage FH: The adult rat hippocampus contains primordial neural stem cells. Mol Cell Neurosci. 1997, 8 (6): 389-404. 10.1006/mcne.1996.0595.
Palmer TD, Ray J, Gage FH: FGF-2-responsive neuronal progenitors reside in proliferative and quiescent regions of the adult rodent brain. Mol Cell Neurosci. 1995, 6 (5): 474-486. 10.1006/mcne.1995.1035.
Palmer TD, Markakis EA, Willhoite AR, Safar F, Gage FH: Fibroblast growth factor-2 activates a latent neurogenic program in neural stem cells from diverse regions of the adult CNS. The Journal of Neuroscience: the official Journal of the Society for Neuroscience. 1999, 19 (19): 8487-8497.
Shihabuddin LS, Horner PJ, Ray J, Gage FH: Adult spinal cord stem cells generate neurons after transplantation in the adult dentate gyrus. The Journal of Neuroscience: the official Journal of the Society for Neuroscience. 2000, 20 (23): 8727-8735.
Kondo T, Raff M: Oligodendrocyte precursor cells reprogrammed to become multipotential CNS stem cells. Science. 2000, 289 (5485): 1754-1757. 10.1126/science.289.5485.1754.
Furuta M, Bridges RS: Effects of maternal behavior induction and pup exposure on neurogenesis in adult, virgin female rats. Brain Res Bull. 2009, 80 (6): 408-413. 10.1016/j.brainresbull.2009.08.011.
Bloch M, Daly RC, Rubinow DR: Endocrine factors in the etiology of postpartum depression. Compr Psychiatry. 2003, 44 (3): 234-246. 10.1016/S0010-440X(03)00034-8.
Tsuchiya KJ, Byrne M, Mortensen PB: Risk factors in relation to an emergence of bipolar disorder: a systematic review. Bipolar Disord. 2003, 5 (4): 231-242. 10.1034/j.1399-5618.2003.00038.x.
Rogers R, Cavanaugh JL, Dolmetsch R: Schedule of Affective Disorders and Schizophrenia, a diagnostic interview in evaluations of insanity: an exploratory study. Psychol Rep. 1981, 49 (1): 135-138. 10.2466/pr0.19188.8.131.52.
Ronald A, Simonoff E, Kuntsi J, Asherson P, Plomin R: Evidence for overlapping genetic influences on autistic and ADHD behaviours in a community twin sample. J Child Psychol Psychiatry. 2008, 49 (5): 535-542. 10.1111/j.1469-7610.2007.01857.x.
Rommelse NN, Franke B, Geurts HM, Hartman CA, Buitelaar JK: Shared heritability of attention-deficit/hyperactivity disorder and autism spectrum disorder. Eur Child Adolesc Psychiatry. 2010, 19 (3): 281-295. 10.1007/s00787-010-0092-x.
Lichtenstein P, Carlstrom E, Rastam M, Gillberg C, Anckarsater H: The genetics of autism spectrum disorders and related neuropsychiatric disorders in childhood. Am J Psychiatry. 2010, 167 (11): 1357-1363. 10.1176/appi.ajp.2010.10020223.
Green EK, Grozeva D, Jones I, Jones L, Kirov G, Caesar S, Gordon-Smith K, Fraser C, Forty L, Russell E, et al: The bipolar disorder risk allele at CACNA1C also confers risk of recurrent major depression and of schizophrenia. Mol Psychiatry. 2010, 15 (10): 1016-1022. 10.1038/mp.2009.49.
Jacobs JP, Ortiz-Lopez A, Campbell JJ, Gerard CJ, Mathis D, Benoist C: Deficiency of CXCR2, but not other chemokine receptors, attenuates autoantibody-mediated arthritis in a murine model. Arthritis Rheum. 2010, 62 (7): 1921-1932.
Bedogni F, Hodge RD, Elsen GE, Nelson BR, Daza RA, Beyer RP, Bammler TK, Rubenstein JL, Hevner RF: Tbr1 regulates regional and laminar identity of postmitotic neurons in developing neocortex. Proc Natl Acad Sci USA. 2010, 107 (29): 13129-13134. 10.1073/pnas.1002285107.
Englund C, Fink A, Lau C, Pham D, Daza RA, Bulfone A, Kowalczyk T, Hevner RF: Pax6, Tbr2, and Tbr1 are expressed sequentially by radial glia, intermediate progenitor cells, and postmitotic neurons in developing neocortex. The Journal of Neuroscience: the Official Journal of the Society for Neuroscience. 2005, 25 (1): 247-251. 10.1523/JNEUROSCI.2899-04.2005.
Sadasivan S, Pond BB, Pani AK, Qu C, Jiao Y, Smeyne RJ: Methylphenidate exposure induces dopamine neuron loss and activation of microglia in the basal ganglia of mice. PLoS One. 2012, 7 (3): e33693-10.1371/journal.pone.0033693.
Irizarry RA, Wang C, Zhou Y, Speed TP: Gene set enrichment analysis made simple. Stat Methods Med Res. 2009, 18 (6): 565-575. 10.1177/0962280209351908.
This work was primarily funded by United States National Institutes of Health Grant R01 MH 085642 to S.G. Additional support was provided by the National Science Foundation Grant IOS-0921706 to S.G. The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript.
The authors declare that they have no competing interests.
BE, MS, TD, and SG conceived and developed the MSET script. BE carried out MSET validation on independent microarray data sets and performed MSET analysis on postpartum LS expression data. BE and SG drafted the manuscript. All authors read and approved the final manuscript.