Abstract
Histone 3 Lysine 9 (H3K9) methylation is known to be associated with pericentric heterochromatin and important in genomic stability. In this study, we show that trimethylation at H3K9 (H3K9me3) is enriched in an adult neural stem cell niche- the subventricular zone (SVZ) on the walls of the lateral ventricle in both rodent and non-human primate baboon brain. Previous studies have shown that there is significant correlation between baboon and human regarding genomic similarity and brain structure, suggesting that findings in baboon are relevant to human. To understand the function of H3K9me3 in this adult neurogenic niche, we performed genome-wide analyses using ChIP-Seq (chromatin immunoprecipitation and deep-sequencing) and RNA-Seq for in vivo SVZ cells purified from baboon brain. Through integrated analyses of ChIP-Seq and RNA-Seq, we found that H3K9me3-enriched genes associated with cellular maintenance, post-transcriptional and translational modifications, signaling pathways, and DNA replication are expressed, while genes involved in axon/neuron, hepatic stellate cell, or immune-response activation are not expressed. As neurogenesis progresses in the adult SVZ, cell fate restriction is essential to direct proper lineage commitment. Our findings highlight that H3K9me3 repression in undifferentiated SVZ cells is engaged in the maintenance of cell type integrity, implicating a role for H3K9me3 as an epigenetic mechanism to control cell fate transition within this adult germinal niche.
Introduction
Chromatin is functionally classified as euchromatin and heterochromatin, which are crucial for epigenetic controls of gene expression. Underlying the specialized chromatin structure around centromere and telomere, H3K9me3 was identified to be heterochromatin-enriched histone code to silence gene expression and prevent chromosomal instability (Czvitkovich et al., ; Lachner et al., ; Peters et al., , ; Black et al., ). For instance, a previous study has shown that loss of H3K9 methylation in Drosophila causes DNA damage in heterochromatin and mitotic defect (Peng and Karpen, ). In mice, loss of H3K9me2/me3 causes the disruption of heterochromatin and increases telomere length (Peters et al., ; Garcia-Cao et al., ; Benetti et al., ). H3K9 methylation is also involved in pluripotency of embryonic stem cells (ESCs) and multipotency of neural precursor cells (NPCs), in which the pluripotent genes (e.g., Nanog, Oct4) and non-neural genes (e.g., GATA4, NODAL) gain H3K9me3 that lead to long-term repression during differentiation of human ESCs into NPCs (Golebiewska et al., ; Hirabayashi and Gotoh, ). Thus, H3K9me3 plays a repressive role in numerous neuronal and non-neuronal genes (Roopra et al., ; Schaefer et al., ) in addition to its known function in genome stability. Among all lysine methyltransferases, the KMT1 family composed of G9a/GLP and Suv39h1/h2 are characterized to be essential for H3K9me1/2 and H3K9me3 modifications, respectively (Czvitkovich et al., ; Lachner et al., ; Peters et al., ; Black et al., ). Studies from KMT1 knock-out mice demonstrated that loss of H3K9 methylation contributes to behavioral abnormalities and cognitive impairment (Schaefer et al., ) in addition to its protective role in genome stability. In this work, we found that H3K9me3 is enriched in the subventricular zone (SVZ), where adult neurogenesis occurs.
The SVZ is the largest neural stem cell niche, which harbors stem/progenitor cells for adult neurogenesis. The SVZ contains slowly dividing neural stem cells (NSCs) with astrocyte-like morphology. In the rodent model, NSCs give rise to transit-amplifying cells, which subsequently give rise to immature neuroblasts. These neuroblasts migrate through the rostral migratory stream (RMS) and generate interneurons in the olfactory bulb (Alvarez-Buylla and Lim, ; Ihrie and Alvarez-Buylla, ). Numerous studies have demonstrated that extracellular signals such as growth factors or morphogens have significant effects on either self-renewal of NSCs or lineage commitment (Doetsch et al., ; Alvarez-Buylla and Lim, ; Zheng et al., ; Jackson et al., ; Ihrie et al., ). Additionally, the intracellular effectors of adult neurogenesis include cell-cycle inhibitors (p16/INK4A and p21) (Molofsky et al., ), transcription factors (Doetsch et al., ; Shi et al., ; Kohwi et al., ; Roybon et al., ; Qu et al., ; Ihrie and Alvarez-Buylla, ), and epigenetic mechanisms (Ming and Song, ). One such epigenetic mechanism includes histone modifications (Ming and Song, ). In this study, we demonstrate that H3K9me3 has a distinct distribution pattern within cell populations in the rodent and non-human primate SVZ. Yet, the molecular targets of H3K9me3 in this adult germinal niche remain unknown. Therefore, we developed a technique to purify subpopulations of SVZ cells from baboon brain (Papio anubis) (Sandstrom et al., ) for genome-wide analysis by using ChIP-Seq. In addition to genes involved in cell cycle and proliferation, we found that H3K9me3 enriched for genes functioning in axon and neuron projection, cellular maintenance/organization, cell signaling, and post-translational acetylation as well. A further integrated ChIP-Seq and RNA-Seq analysis revealed that 35% of H3K9me3-enriched genes are silenced, many of which are known to function in neuronal-, hepatic-, and immunological-cell type activation. In light of previous studies showing that H3K9me3 is a chromatin repressive mark, we anticipate that H3K9me3 is critical for the maintenance of cell identity within this adult neurogenic niche through its repressive function to protect against improper lineage differentiation within the SVZ.
Results
H3K9me3 is expressed in germinal zones within the adult brain
Using co-immunostaining with antibodies specific for H3K9me3 (a pericentric chromatin staining pattern, Figure 1B) and cell type specific markers in the SVZ (Figure 1I) for the gross anatomy of 8-week (P56) old adult mouse brain (Figure 1A), we found that H3K9me3-positive cells are co-localized with GFAP which labels quiescent and active NSCs (Figure 1C), Vimentin positive active NSC (Figure 1D), and PSA-NCAM positive neuroblast populations as well (Figure 1E). As Mash1 is commonly used in the mouse to denote cells as “transient amplifying cells,” we applied co-immunostaining of Mash1 and H3K9me3 in the mouse SVZ, and found colocalization of H3K9me3 and Mash1 (Figure 1F). We then performed a 2 h EdU administration in mice to label quickly dividing transit-amplifying cells and neuroblasts, and found colocalization of EdU and H3K9me3 (Figures 1G,H). These results show that H3K9me3 is present in undifferentiated SVZ cells within rodent brain.
Figure 1
The comparisons of brain volume and structure across primate species and human have shown significant correlation between baboon and human (Kochunov et al., ; Rogers et al., ), we therefore examined the H3K9me3 distribution pattern in the baboon SVZ (Figure 2A). It is worth noticing that GFAP-positive astrocytic ribbon lies in this niche and also extends toward the lateral ventricle (Figures 2B,C), resembling the previous finding in adult human brain (Sanai et al., ). Thus, the overall architecture of baboon SVZ geographically represents the SVZ in adult human brain (Sanai et al., , ), corroborating the extent to which findings in baboon SVZ are relevant to human SVZ. We found that H3K9me3 is associated with GFAP- and Vimentin-positive NSCs and also PSA-NCAM-positive neuroblasts in the baboon SVZ (Figures 2D–G), suggesting that H3K9me3 has function in these cell populations.
Figure 2
To further quantify the percentages of co-localization between H3K9me3 and SVZ subpopulations, we carried out flow cytometry analysis for dissociated SVZ cells after micro-dissection of SVZ from baboon brain. Flow cytometry analysis reveals approximately 45 and 20% of GFAP- and Doublecortin (DCX)-positive populations contain H3K9me3, respectively (Figure 3). Consistent with immunostaining results, the great majority of Vimentin- and PSA-NCAM-positive cells (~95%) are colocalized with H3K9me3 (Figure 3). This quantification confirms that H3K9me3 is enriched in undifferentiated SVZ cells, while the abundance of H3K9me3 varies across different cell populations within the SVZ.
Figure 3
H3K9me3 epigenetic landscape in the SVZ cells
While results from our study highlight the distinct localization pattern of H3K9me3 in the SVZ of rodent and non-human primate, the molecular targets of H3K9me3 in this germinal niche remain to be determined. Currently, most of the genome-wide associated studies in neural stem/progenitor cells are conducted in culture to obtain sufficient materials for analyses. However, the in vitro cell culture system cannot fully recapitulate in vivo epigenetic landscape since the metabolites yielded from cultured condition including acetyl and methyl donors can alter the status of histone acetylation and methylation (Black et al., ). To ascertain characteristics of SVZ cells as they exist in vivo for identification of the genomic loci that carry the H3K9me3 modification, we developed a technique to purify SVZ cells directly from the baboon brain within a short post-mortem interval (<20 min). We utilized conjugated Dynabeads with antibodies against SVZ cell type-specific markers to purify dissociated cells following dissection of the SVZ (Figure 4A) (Sandstrom et al., ). This technical innovation preserves the nature of distinct SVZ cell types and is ideal for genome-wide analysis to uncover H3K9me3 enriched loci in the SVZ through ChIP-Seq with antibodies specific for H3K9me3. Additionally, this approach ensures that H3K9me3 positive cells from the adjacent striatum are excluded from the genome-wide analysis of the SVZ cells. DNA obtained from each ChIP pull-down was sequenced to high depth (200 million tags; 36 bases) by using Illumina HiSeq2000 sequencer (Figure 4A). We did additional runs of ChIP-Seq as replicates with different vendor's antibody H3K9me3 and independent sample preparation. This approach yielded a confident list of H3K9me3 enriched loci across independent sets of deep-Seq. Because the baboon gene annotation is not currently available, we alternatively compiled, processed, and aligned the sequence reads to the Jan. 2006 rhesus macaque (Macaca mulatta) draft assembly, Mmul_051212 (Gibbs et al., ) (http://www.ncbi.nlm.nih.gov/genome/assembly/237568/) and the UCSC genome browser version: rheMac2 (http://hgdownload.cse.ucsc.edu/downloads.html#rhesus). Of note, there is only 2% difference at the genomic level between baboon and rhesus macaque, thus, we were able to identify 863 unique H3K9me3-enriched genes from independent sets of ChIP-Seq (FDR = 0.05) by using the MACS2 and closest-features program from the BEDOPS tool set (Neph et al., ) (Figure 4B; Supplemental Table 1i). Gene ontology with functional annotation analysis reveals that molecular functions associated with these genes include acetylation, axon/neuron projection, and protein targeting/import (Figure 4C; Supplemental Table 1ii) with connection to neurological disorders. While imprinted genes are known to be involved in broad aspects of biology including brain function, imprinting dysregulation has been associated with several neurodevelopmental and neurological disorders. To elucidate whether there are imprinted genes among H3K9me3 targets in the SVZ cells, we undertook an overlap comparison between the lists of identified imprinted genes (Luedi et al., ) and H3K9me3 targets identified from our ChIP-Seq analysis. We found a total of 11 genes that are imprinted and enriched with H3K9me3. These genes have known functions in CNS fate commitment, glial differentiation, neural projection, and neurite outgrowth, suggesting a role of H3K9me3 in collaboration with imprinting mechanism to maintain the populations of undifferentiated SVZ cells. Among the H3K9me3-enriched genes (n = 863), the top biological networks predicted by Ingenuity Pathways Analysis (IPA) are involved in (1) cell morphology and cellular assembly; (2) post-transcriptional and post-translational modifications; (3) protein synthesis; (4) cell cycle; and (5) cellular growth and proliferation (Figure 5; Supplemental Table 1iv). Additionally, the top biological pathways by IPA prediction include protein ubiquitination as well as signaling pathways involving AKT, BAX, c-JUN, MDM2, p300, P53, PP2A, and PTEN (Figure 6; Supplemental Figures 1–3).
Figure 4
Figure 5
Figure 6
To further characterize the functional consequence of H3K9me3 enrichment on target genes, we performed RNA-Seq for purified baboon SVZ cells and applied integrated analysis between ChIP-Seq and RNA-Seq. We identified 562 genes enriched with H3K9me3 are detectable by RNA-Seq, which are associated with cellular assembly/maintenance/organization, post-transcriptional and translational modifications, signaling pathways, and DNA replication. We also found that 301 of H3K9me3-enriched genes involved in axon/neuron or hepatic stellate cell activation and immunological response are not detectable by RNA-Seq (p-value = 2.2e-16) (Figure 4D; Supplemental Table 1iii). Given that H3K9me3 can act as short-term and long-term repression, we anticipate 35% (301/863) of H3K9me3-enriched genes to be under long-term repression. Lastly, H3K9me3-enriched genes in the acetylation category from our GO analysis (Figure 4C) were previously characterized with roles in histone acetylation that regulate gene expression or in protein acetylation that regulates protein stability, localization, and interactions with other molecules. Numerous H3K9me3-enriched genes associated with cellular assembly, maintenance, organization, and signaling are known to be regulated by histone acetylation (Cho and Cavalli, ). In addition to histone acetylation, a set of H3K9me3-enriched genes with connection to neurological disorders linked to basal ganglia malfunction are known to have post-translational acetylation on the protein as well (Lopez-Atalaya et al., ; Valor et al., ,). To further explore the relationship between H3K9me3 enrichment and the set of 162 H3K9me3-enriched genes in the acetylation category, we carried out an integrated analysis and found that 138 of H3K9me3-enriched genes in acetylation category are detectable by RNA-Seq, whereas 24 genes are not detectable by RNA-Seq (Supplemental Figure 4). IPA analysis reveals that the 24 undetectable genes are involved in protein degradation and endocrine system development, while the 138 detectable genes are associated with signaling, nervous system development and function, as well as cellular assembly, maintenance, morphology, and organization. Since histone acetylation is highly associated with activation of gene expression and protein acetylation is critical for stability and localization of protein, we anticipate that the interplay between H3K9me3 and acetylation on common sets of genes may modulate the balance of long- and short-term repression.
Discussion
Negative regulation of transcription is important in establishing and maintaining cell-type specific gene expression patterns. One such negative regulation can be achieved by H3K9me3 epigenetic repression. During development, many neuronal genes are subject to repression outside the nervous system to maintain neuronal specificity. For instance, neuronal genes in terminally differentiated fibroblasts are silenced. In this regard, the chromatin repressive mark H3K9me3 was identified to participate in either short- or long-term repression. Although epigenetic mechanisms controlling cell fate specification have been intensively studied in the developing embryonic central nervous system (Lim et al., ; Hirabayashi and Gotoh, ), whether the specification or maintenance of multipotency of neural progenitors in the adult brain relies on similar epigenetic regulations as in developing brains remains unknown. In this study, we demonstrate that H3K9me3 is associated with NSCs and is persistently enriched in the neuroblast population in both the adult rodent and baboon SVZ. Given the structural and genomic correlations between baboon and human, findings in baboon are more relevant to human with regards to the molecular targets of H3K9me3 in this germinal niche.
To uncover the molecular targets of H3K9me3 in adult baboon SVZ, we developed a technique to overcome the spatial complexity of the SVZ and to purify undifferentiated SVZ cells for genome-wide analysis. Our approach will be of considerable interest to those applying the genome-wide cutting edge techniques for cell lineage study. We identified that H3K9me3 is enriched for genes involved in the network of cell cycle and the signaling pathways of PTEN, MDM2, and AKT. These findings implicate that epigenetic regulation through H3K9me3 is associated with cellular growth and proliferation. Intriguingly, significant sets of H3K9me3-enriched genes with identified functions in axon/neuron projections and hepatic or immunological activity are not expressed. As mature axon/neuron and non-neuronal genes should be silenced in this germinal niche, we reason that the repression through H3K9me3 is a mechanism to secure the identity of undifferentiated SVZ cells. In summary, as neurogenesis proceeds, numerous genes must go through active, poised, and repressed states to coordinate lineage commitment. The extent of neurogenesis can be inferred from integrated regulations through signaling pathways and different genetic/epigenetic mechanisms. Our findings suggest that H3K9me3 regulates adult neurogenesis, at least in part, by suppressing a subset of genes to maintain SVZ niche properties and to tightly regulate lineage specification in order to coordinate proper timing for adult neurogenesis.
Materials and methods
All animal experiments were approved by the guidelines of the Institutional Animal Care and Use Committee of the University of Texas at San Antonio (UTSA) and Texas Biomedical Research Institute/Southwest National Primate Research Center (SNPRC) at San Antonio.
Immunostaining and confocal imaging
Mouse
Mice were trans-cardially perfused and fixed using 1X PBS(−) and 4% paraformaldehyde (PFA), respectively. Brains were cryoprotected in 30% sucrose prior to OCT embedding. For co-localization staining, 12 μm frozen sections were processed for immunostaining with antibodies against H3K9me3 (Upstate #07-422; 1:500), Glial Fibrillary Acidic Protein (GFAP)—clone GA5 (Millipore #MAB3402; 1:500), Vimentin (Sigma V2258; 1:500), Polysialic Acid-NCAM (PSANCAM)—clone 2-2b (Millipore #MAB5324; 1:500), Mash1(Abcam #ab38556; 1:500), and EdU 5-ethynyl-2′-deoxyuridine (Life Technology, #C10337 Click-iT® EdU).
Baboon
Coronal slices of fresh baboon forebrain were taken to obtain the SVZ and adjacent brain regions, which were subsequently, fixed in 4% PFA overnight and then cryoprotected in 30% sucrose before OCT embedding. For co-localization staining, 60 μm floating sections were processed for immunostaining with the antibodies listed above.
Vectashield with DAPI (Vector Laboratories Ltd # H-1200) was used for mounting medium and nuclear counter stain. Secondary antibodies AlexaFIuor 488 (Molecular Probes, 1:1000) and AlexaFluor 594 (Molecular Probes 1:1000) were utilized for fluorescent labeling. Mouse SVZ images were acquired under a Zeiss510 confocal microscope (40X and 100X oil immersion objectives). Baboon SVZ images were taken under a Zeiss710 two-photon confocal microscope (40X and 100X oil immersion objectives). Z-stacks were projected to single plane using Zen 2012 (Carl Zeiss Microscopy; black edition). All images were analyzed using ImageJ (NIH; version 1.47).
Cell types purification for ChIP-Seq analysis
Antibody against SVZ cell type markers, such as GFAP, Vimentin, PSA-NCAM, or Doublecortin was manually conjugated to Dynabeads (Dynabeads®-Protein A, Life Technology). We then used the Dynabeads-conjugated antibody to purify cells immediately dissociated from SVZ microdissection. Briefly, cells from fresh dissected baboon SVZ were immediately dissociated with Accutase, subsequently equilibrated in binding buffer containing phosphate-buffered saline (PBS), saponin, 1X protease inhibitor cocktail (Roche), and subjected to Dynabeads-conjugated antibody purification. The purified cells were crosslinked in 1.1% formaldehyde before chromatin shearing by Diagenode Bioruptor. The resulting sheared chromatin fragments in a size range between 200–500 base pairs were then incubated with H3K9me3 antibody-conjugated Protein A Dynabeads overnight. Additional runs of ChIP-Seq were performed as replicates. Each replicate included independent sample preparation and a different vendor H3K9me3 antibody to ensure no bias was introduced through a specific antibody (Millipore/Upstate #07-442; Active Motif #39162; Life Technology Dynabeads protein A).
For normalization, the aliquot of sheared chromatin fragments were incubated with antibody against total histone 3- conjugated Protein A Dynabeads (unmodified H3 antibody, Millipore #05-499; 1:1000). Subsequently, enriched chromatin fragments were eluted, de-crosslinked purified for library preparation (Illumina Library Kit), and sequenced with 200 million tags through Illumina HiSeq2000 sequencer. 56 ng of H3K9me3 ChIP-DNA were applied for library preparation and 5pM of libraries were loaded into Illumina sequencer. The resulting 121,743,940 pass filter reads were aligned for peak calls.
Sequence alignment and peak calling
The Rhesus macaque (rheMac2) gene annotation derived from the NCBI RefSeq project (http://nar.oxfordjournals.org/content/33/suppl_1/D501.full) was constructed and is maintained at the UCSC genome browser. Alignments were generated by using the Bowtie alignment program version 0.12.7., with a maximum of 2 mis-matches allowed in the mapping reads. Aligned read enrichments were detected using the MACS2 peak finding program. Peak calls were generated by normalizing to unmodified H3 DNA, and using an FDR of 0.05. ChIP-Sequencing raw data and processed peak calls have been deposited to GEO under the accession ID GSE59074.
Go, network, and pathway analysis
Duplicate gene references were removed prior to GO, network, and pathway analyses. DAVID Functional Annotation Tool (DAVID Bioinformatics Resources 6.7, NIAID/NIH) was utilized to perform GO analysis and significance cutoff was set at a p-value of <0.05. Network and canonical pathway analyses were performed using Ingenuity Pathway Analysis (IPA) (Ingenuity® Systems, Redwood City, CA, USA). The top 5 networks and pathways presented in this paper were determined by IPA (p-value < 0.003 and Fischer's Test Score > 34, respectively).
RNA-Seq analysis
Total RNA was extracted from purified baboon SVZ cells using TRIzol reagent and sequencing libraries were generated with Illumina RNA-Seq library kit. Paired-end RNA-deepSeq (76 base pair; >300 million tag reads; 269,081,636 mapped reads) were aligned to hg19. DESeq was used to normalize raw read counts; and Cufflink reports read counts and estimated FPKM (fragments per kilobase of exon per million fragments mapped; http://cufflinks.cbcb.umd.edu/faq.html#fpkm). Genes with expression values >1 FPKM were considered for subsequent analyses. RNA-Sequencing data have been deposited to GEO under the accession ID GSE58527.
Flow cytometry
Disassociated SVZ cells were incubated on ice for 10 min with anti-mouse CD16/CD32 (clone 2.4G2; BD Pharmingen) to block FcRs. Cells were then incubated on ice for 1 h with primary antibodies (anti-H3K9me3, GFAP, Vimentin, PSA-NCAM, and Doublecortin), and subsequently labeled with Alexa Fluor 488- and PE-conjugated secondary antibodies (BD Biosciences). Primary antibody resources: H3K9me3 (Millipore/Upstate #07-442; 1:500); Doublecortin (Millipore clone2G5 #MABN707; 1:250); Polysialic Acid-NCAM, clone 2-2b (Millipore #MAB5324, Lot# 1966892; 1:250); Glial Fibrillary Acidic Protein, clone GA5 (Millipore #MAB3402, Lot#1993774; 1:250); Vimentin (Sigma V2258; 1:250); For analysis, controls including the positive controls stained with each antibody separately, isotype controls, and the unstained cells were used for gate compensation. Flow data was acquired on a LSR-II flow cytometer (BD Biosciences) configured with an argon 488 laser with a 505 LP dichroic and 525/50 filter to detect Alexa fluor 488 and a green 510 laser with a 735 LP dichroic and a 575/26 filter to detect PE. Compensation and data analysis was performed using FlowJo software (Tree Star, Inc, Ashland, OR).
Conflict of interest statement
The Guest Associate Editor Michael T. Chin declares that, despite being affiliated to the same Institution as the author Richard S. Sandstrom, the review process was handled objectively. The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Statements
Acknowledgments
We thank Texas Biomedical Research Institute/Southwest National Primate Research Center (SNPRC) at San Antonio for baboon tissues. Authors are grateful for technical support from Drs. Clyde Phelix, Jera Pecotte, and Eric Haugen. Sequencing was performed at the High Throughput Genomics Center in the Department of Genome Sciences, University of Washington, Seattle, WA. This work received support from Computational System Biology Core and Imaging Core, funded by the National Institute on Minority Health and Health Disparities (G12MD007591) from the National Institutes of Health. Yufeng Wang is supported by NIH grants GM100806 and AI080579. This project was supported in part by the SPORE Grant 5 P50 CA097257-10 to Mitchel S. Berger and was supported by the 5P30GM092334 from the National Institutes of Health and TRAC award (Tenure-track Research Award Competition) to Chin-Hsing Annie Lin.
Conflict of interest
The Guest Associate Editor Michael T. Chin declares that, despite being affiliated to the same Institution as the author Richard S. Sandstrom, the review process was handled objectively. The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Supplementary material
The Supplementary Material for this article can be found online at: http://www.frontiersin.org/journal/10.3389/fgene.2014.00252/abstract
Supplemental Figure 1Protein Ubiquitination Pathway. One of top four canonical pathways predicted by Ingenuity Pathway Analysis for H3K9me3 enriched genes in baboon SVZ cells. The shaded focus genes (red highlight) were enriched with H3K9me3 identified by ChIP-Seq analysis. The node shape reflects the role of each element in the pathway and the direction and arrowhead shapes of each edge represent different types of interactions (see key at bottom-left panel).
Supplemental Figure 2P53 signaling pathway. One of top four canonical pathways predicted by Ingenuity Pathway Analysis for H3K9me3 enriched genes in baboon SVZ cells. The shaded focus genes (red highlight) were enriched with H3K9me3 identified by ChIP-Seq analysis. The node shape reflects the role of each element in the pathway and the direction and arrowhead shapes of each edge represent different types of interactions (see key in Supplemental Figure 1).
Supplemental Figure 3Hepatic stellate cell activation pathway. One of top four canonical pathways predicted by Ingenuity Pathway Analysis for H3K9me3 enriched genes in baboon SVZ cells. The shaded focus genes (red highlight) were enriched with H3K9me3 identified by ChIP-Seq analysis. The node shape reflects the role of each element in the pathway and the direction and arrowhead shapes of each edge represent different types of interactions (see key in Supplemental Figure 1).
Supplemental Figure 4The integrated analysis between RNA-Seq and the set of H3K9me3-enriched genes in the acetylation category from our GO analysis. 138 of H3K9me3-enriched genes in acetylation category are detectable by RNA-Seq, whereas 24 genes are not detectable by RNA-Seq. For characterization of these detectable or undetectable H3K9me3-enriched genes in baboon SVZ cells, IPA prediction reveals the top canonical networks and pathways as shown in the box.
Supplemental Table 1(i) Summary of Enriched Loci of H3K9me3; (ii) Gene Ontology by DAVID for H3K9me3-enriched genes in Baboon SVZ Cells; (iii) Integrated analyses for ChIP-Seq and RNA-Seq; (iv) A subset of H3K9me3-enriched genes was characterized as imprinting genes.
References
1
Alvarez-BuyllaA.LimD. A. (2004). For the long run: maintaining germinal niches in the adult brain. Neuron41, 683–686. 10.1016/S0896-6273(04)00111-4
2
BenettiR.Garcia-CaoM.BlascoM. A. (2007). Telomere length regulates the epigenetic status of mammalian telomeres and subtelomeres. Nat. Genet. 39, 243–250. 10.1038/ng1952
3
BlackJ. C.Van RechemC.WhetstineJ. R. (2012). Histone lysine methylation dynamics: establishment, regulation, and biological impact. Mol. Cell48, 491–507. 10.1016/j.molcel.2012.11.006
4
ChoY.CavalliV. (2014). HDAC signaling in neuronal development and axon regeneration. Curr. Opin. Neurobiol. 27C, 118–126. 10.1016/j.conb.2014.03.008
5
CzvitkovichS.SauerS.PetersA. H.DeinerE.WolfA.LaibleG.et al. (2001). Over-expression of the SUV39H1 histone methyltransferase induces altered proliferation and differentiation in transgenic mice. Mech. Dev. 107, 141–153. 10.1016/S0925-4773(01)00464-6
6
DoetschF.PetreanuL.CailleI.Garcia-VerdugoJ. M.Alvarez-BuyllaA. (2002). EGF converts transit-amplifying neurogenic precursors in the adult brain into multipotent stem cells. Neuron36, 1021–1034. 10.1016/S0896-6273(02)01133-9
7
Garcia-CaoM.O'SullivanR.PetersA. H.JenuweinT.BlascoM. A. (2004). Epigenetic regulation of telomere length in mammalian cells by the Suv39h1 and Suv39h2 histone methyltransferases. Nat. Genet. 36, 94–99. 10.1038/ng1278
8
GibbsR. A.RogersJ.KatzeM. G.BumgarnerR.WeinstockG. M.MardisE. R.et al. (2007). Evolutionary and biomedical insights from the rhesus macaque genome. Science316, 222–234. 10.1126/science.1139247
9
GolebiewskaA.AtkinsonS. P.LakoM.ArmstrongL. (2009). Epigenetic landscaping during hESC differentiation to neural cells. Stem Cells27, 1298–1308. 10.1002/stem.59
10
HirabayashiY.GotohY. (2010). Epigenetic control of neural precursor cell fate during development. Nat. rev. Neurosci. 11, 377–388. 10.1038/nrn2810
11
IhrieR. A.Alvarez-BuyllaA. (2011). Lake-front property: a unique germinal niche by the lateral ventricles of the adult brain. Neuron70, 674–686. 10.1016/j.neuron.2011.05.004
12
IhrieR. A.ShahJ. K.HarwellC. C.LevineJ. H.GuintoC. D.LezametaM.et al. (2011). Persistent sonic hedgehog signaling in adult brain determines neural stem cell positional identity. Neuron71, 250–262. 10.1016/j.neuron.2011.05.018
13
JacksonE. L.Garcia-VerdugoJ. M.Gil-PerotinS.RoyM.Quinones-HinojosaA.VandenBergS.et al. (2006). PDGFR alpha-positive B cells are neural stem cells in the adult SVZ that form glioma-like growths in response to increased PDGF signaling. Neuron51, 187–199. 10.1016/j.neuron.2006.06.012
14
KochunovP.GlahnD. C.FoxP. T.LancasterJ. L.SaleemK.ShelledyW.et al. (2010). Genetics of primary cerebral gyrification: heritability of length, depth and area of primary sulci in an extended pedigree of Papio baboons. Neuroimage53, 1126–1134. 10.1016/j.neuroimage.2009.12.045
15
KohwiM.OsumiN.RubensteinJ. L.Alvarez-BuyllaA. (2005). Pax6 is required for making specific subpopulations of granule and periglomerular neurons in the olfactory bulb. J. Neurosci. 25, 6997–7003. 10.1523/JNEUROSCI.1435-05.2005
16
LachnerM.O'CarrollD.ReaS.MechtlerK.JenuweinT. (2001). Methylation of histone H3 lysine 9 creates a binding site for HP1 proteins. Nature410, 116–120. 10.1038/35065132
17
LimD. A.HuangY. C.SwigutT.MirickA. L.Garcia-VerdugoJ. M.WysockaJ.et al. (2009). Chromatin remodelling factor Mll1 is essential for neurogenesis from postnatal neural stem cells. Nature458, 529–533. 10.1038/nature07726
18
Lopez-AtalayaJ. P.ItoS.ValorL. M.BenitoE.BarcoA. (2013). Genomic targets, and histone acetylation and gene expression profiling of neural HDAC inhibition. Nucleic Acids Res. 41, 8072–8084. 10.1093/nar/gkt590
19
LuediP. P.HarteminkA. J.JirtleR. L. (2005). Genome-wide prediction of imprinted murine genes. Genome Res. 15, 875–884. 10.1101/gr.3303505
20
MingG. L.SongH. (2011). Adult neurogenesis in the mammalian brain: significant answers and significant questions. Neuron70, 687–702. 10.1016/j.neuron.2011.05.001
21
MolofskyA. V.SlutskyS. G.JosephN. M.HeS.PardalR.KrishnamurthyJ.et al. (2006). Increasing p16INK4a expression decreases forebrain progenitors and neurogenesis during ageing. Nature443, 448–452. 10.1038/nature05091
22
NephS.KuehnM. S.ReynoldsA. P.HaugenE.ThurmanR. E.JohnsonA. K.et al. (2012). BEDOPS: high-performance genomic feature operations. Bioinformatics28, 1919–1920. 10.1093/bioinformatics/bts277
23
PengJ. C.KarpenG. H. (2009). Heterochromatic genome stability requires regulators of histone H3 K9 methylation. PLoS Genet. 5:e1000435. 10.1371/journal.pgen.1000435
24
PetersA. H.MermoudJ. E.O'CarrollD.PaganiM.SchweizerD.BrockdorffN.et al. (2002). Histone H3 lysine 9 methylation is an epigenetic imprint of facultative heterochromatin. Nat. Genet. 30, 77–80. 10.1038/ng789
25
PetersA. H.O'CarrollD.ScherthanH.MechtlerK.SauerS.SchoferC.et al. (2001). Loss of the Suv39h histone methyltransferases impairs mammalian heterochromatin and genome stability. Cell107, 323–337. 10.1016/S0092-8674(01)00542-6
26
QuQ.SunG.LiW.YangS.YeP.ZhaoC.et al. (2010). Orphan nuclear receptor TLX activates Wnt/beta-catenin signalling to stimulate neural stem cell proliferation and self-renewal. Nat. Cell Biol. 12, 31–40; sup. pp 1–9. 10.1038/ncb2001
27
RogersJ.KochunovP.ZillesK.ShelledyW.LancasterJ.ThompsonP.et al. (2010). On the genetic architecture of cortical folding and brain volume in primates. Neuroimage53, 1103–1108. 10.1016/j.neuroimage.2010.02.020
28
RoopraA.QaziR.SchoenikeB.DaleyT. J.MorrisonJ. F. (2004). Localized domains of G9a-mediated histone methylation are required for silencing of neuronal genes. Mol. Cell14, 727–738. 10.1016/j.molcel.2004.05.026
29
RoybonL.HjaltT.StottS.GuillemotF.LiJ. Y.BrundinP. (2009). Neurogenin2 directs granule neuroblast production and amplification while NeuroD1 specifies neuronal fate during hippocampal neurogenesis. PLoS ONE4:e4779. 10.1371/journal.pone.0004779
30
SanaiN.NguyenT.IhrieR. A.MirzadehZ.TsaiH. H.WongM.et al. (2011). Corridors of migrating neurons in the human brain and their decline during infancy. Nature478, 382–386. 10.1038/nature10487
31
SanaiN.TramontinA. D.Quinones-HinojosaA.BarbaroN. M.GuptaN.KunwarS.et al. (2004). Unique astrocyte ribbon in adult human brain contains neural stem cells but lacks chain migration. Nature427, 740–744. 10.1038/nature02301
32
SandstromR. S.ForetM. R.GrowD. A.HaugenE.RhodesC. T.CardonaA. E.et al. (2014). Epigenetic regulation by chromatin activation mark H3K4me3 in primate progenitor cells within adult neurogenic niche. Sci. Rep. 4:5371. 10.1038/srep05371
33
SchaeferA.SampathS. C.IntratorA.MinA.GertlerT. S.SurmeierD. J.et al. (2009). Control of cognition and adaptive behavior by the GLP/G9a epigenetic suppressor complex. Neuron64, 678–691. 10.1016/j.neuron.2009.11.019
34
ShiY.Chichung LieD.TaupinP.NakashimaK.RayJ.YuR. T.et al. (2004). Expression and function of orphan nuclear receptor TLX in adult neural stem cells. Nature427, 78–83. 10.1038/nature02211
35
ValorL. M.GuirettiD.Lopez-AtalayaJ. P.BarcoA. (2013a). Genomic landscape of transcriptional and epigenetic dysregulation in early onset polyglutamine disease. J. Neurosci. 33, 10471–10482. 10.1523/JNEUROSCI.0670-13.2013
36
ValorL. M.VioscaJ.Lopez-AtalayaJ. P.BarcoA. (2013b). Lysine acetyltransferases CBP and p300 as therapeutic targets in cognitive and neurodegenerative disorders. Curr. Pharm. Des. 19, 5051–5064. 10.1523/JNEUROSCI.0670-13.2013
37
ZhengW.NowakowskiR. S.VaccarinoF. M. (2004). Fibroblast growth factor 2 is required for maintaining the neural stem cell pool in the mouse brain subventricular zone. Dev. Neurosci. 26, 181–196. 10.1159/000082136
Summary
Keywords
epigenetics, histone methylation, pericentric chromatin, Papio anubis, SVZ, neurogenesis
Citation
Foret MR, Sandstrom RS, Rhodes CT, Wang Y, Berger MS and Lin C-HA (2014) Molecular targets of chromatin repressive mark H3K9me3 in primate progenitor cells within adult neurogenic niches. Front. Genet. 5:252. doi: 10.3389/fgene.2014.00252
Received
21 February 2014
Accepted
10 July 2014
Published
30 July 2014
Volume
5 - 2014
Edited by
Michael T. Chin, University of Washington, USA
Reviewed by
Francois P. Guillemot, National Institute of Medical Research, UK; Ian C. G. Weaver, Dalhousie University, Canada
Copyright
© 2014 Foret, Sandstrom, Rhodes, Wang, Berger and Lin.
This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) or licensor are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Chin-Hsing Annie Lin, Department of Biology, University of Texas at San Antonio, One UTSA Circle, San Antonio, TX 78249, USA e-mail: annie.lin@utsa.edu
†These authors have contributed equally to this work.
This article was submitted to Epigenomics and Epigenetics, a section of the journal Frontiers in Genetics.
Disclaimer
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.