- Research article
- Open Access
Transcriptome analysis of human brain tissue identifies reduced expression of complement complex C1Q Genes in Rett syndrome
BMC Genomics volume 17, Article number: 427 (2016)
MECP2, the gene mutated in the majority of Rett syndrome cases, is a transcriptional regulator that can activate or repress transcription. Although the transcription regulatory function of MECP2 has been known for over a decade, it remains unclear how transcriptional dysregulation leads to the neurodevelopmental disorder. Notably, little convergence was previously observed between the genes abnormally expressed in the brain of Rett syndrome mouse models and those identified in human studies.
Here we carried out a comprehensive transcriptome analysis of human brain tissue from Rett syndrome brain using both RNA-seq and microarrays.
We identified over two hundred differentially expressed genes, and identified the complement C1Q complex genes (C1QA, C1QB and C1QC) as a point of convergence between gene expression changes in human and mouse Rett syndrome brain.
The results of our study support a role for alterations in the expression level of C1Q complex genes in RTT pathogenesis.
Rett syndrome (RTT) is an X-linked neurodevelopmental disorder primarily affecting girls at a frequency of 1/10,000 live female births. The core manifestations of RTT include intellectual disability, intractable seizures, spasticity and stereotypic hand movements . More than 95 % of classic RTT cases are caused by sporadic mutations in the gene encoding methyl-CpG binding protein 2 (MECP2) .
MECP2 encodes a nuclear protein, which belongs to the methyl-CpG binding protein family and comprises several distinct domains (Fig. 1): a methyl-binding domain (MBD), intervening domain (ID), transcriptional repression domain (TRD), and C-terminal domain (CTD) (Hendrich and Bird, 1998). MECP2 binds with high affinity to mCG dinucleotides, and it has been recently shown to also bind mCA dinucleotides . Although present in most somatic cells, MECP2 expression is most abundant in the brain, with a high ratio of neuronal to glial expression levels .
How MECP2 dysfunction leads to the RTT phenotype remains unclear, despite extensive work on multiple mouse models that partially recapitulate the neurological abnormalities of RTT . Gene expression studies of MECP2 transgenic and knockout mice have shown that MECP2 both activates and represses transcription [5, 6] in the mouse brain. While initial studies pointed out to a transcriptional repressor role of MECP2 [7, 8], involving recruitment of NCoR, HDAC3 and Sin3a, more recent data have uncovered an important role of MECP2 as a transcriptional activator . Remarkably, MECP2’s abundance in neuronal nuclei is similar to that of the histone octamer, and thus MECP2 binds globally across the methylated genome .
Given the transcriptional regulatory role of MECP2, the underlying mechanism of the RTT phenotype likely results from the dysregulation of MECP2 target genes. Therefore, identifying genes dysregulated in human RTT brain could provide important clues into the mechanism of the disease. However, limited data is available on transcriptome changes in human RTT brain. Previous postmortem brain studies using microarray platforms either failed to identify genes that passed statistical significance criteria after false discovery rate correction or lacked age-matched controls [9–11], and thus little consensus was observed between distinct human studies or between human and mouse data.
Here we carried out the first RNA-seq analysis of brain tissue from RTT patients, and identified several hundreds of differentially expressed genes. We also found a significant overlap between genes downregulated in Rett brain samples and genes activated by MECP2 in mouse models, which included genes involved in the complement cascade.
Results and discussion
To investigate genome-wide transcriptome changes in human RTT brain, we obtained postmortem tissue from four RTT cases and four age-, sex-, and ethnicity- matched controls (Methods). For each individual brain, we obtained tissue samples from both frontal and temporal cortex (Fig. 1). After quality control assessment of RNA samples two RTT and two control samples were eliminated from further analyses due to low RNA quality (see Methods). The remaining samples, include three Rett cases with the following MECP2 mutations (Fig. 1a): R1 carries a splice site mutation at the intron3-exon4 junction (c.378-2 A > G), while R2 and R3 carry exon4 mutations: c.763 C > T (p. 255 R > X; transcriptional repression domain) and c.451 G > T (p.151 D > Y; C-terminal domain) respectively. Despite the distinct mutations, the clinical picture was very similar for the three RTT cases. R4852 was diagnosed with RTT with secondary generalised epilepsy and kyphoscoliosis; R4516 was diagnosed with RTT and epilepsy with complex partial seizures; R1815 was diagnosed with RTT, complex partial seizures and scoliosis.
For microarray analyses, each RNA sample was analyzed separately (Fig. 1b, Methods), while for RNA-seq, due to limited amounts of RNA, frontal and temporal samples from the same individual were pooled as shown in Fig. 1b. We obtained an average of 62 million sequencing reads per sample, a high sequencing depth that allowed us to detect lowly expressed transcripts.
Since dissected tissue consists of a mixture of neuronal and glial cells, the results of gene expression analyses could potentially be skewed by differences in cell-type composition across samples. This issue, although important in all transcriptome studies, is particularly important in studies with a small number of samples. However, the potentially confounding effect of cellular composition is rarely addressed due to the difficulty of obtaining good quality RNA from postmortem tissues after laser-microdissection or cell sorting. Thus we used an in silico decomposition method  to estimate the proportion of neurons and astrocytes in each dissected tissue sample, based on the microarray expression profiles and neuronal and astrocyte cell markers  (Methods). We found that although the proportion of neurons showed only minor variations across samples (Fig. 2a, Additional file 1: Table S1, range: 43 % to 60 %, p = 0.06, two-tailed t-test), the differences in cell type composition affected the clustering of expression profiles (Fig. 2b). The RTT samples that were more similar to controls in terms of cell composition (R1815) clustered with controls, while the rest of RTT samples formed a distinct group (Fig. 2b). To address this issue, we used a normalization method proposed for eliminating unwanted variation, such as experimental co-variates, from transcriptome data (Remove Unwanted Variation (RUV) normalization [14, 15], Methods). RUV estimates a parameter of unwanted variation in the data, i.e. variation not related to the case–control differences, and uses this parameter as covariate in the generalized linear model applied for differential expression (DE) analysis. After applying RUV normalization to both the RNA-seq and microarray data, we found that sample clustering no longer depended on cell-type composition (Fig. 2).
We next analyzed the RNA-seq data to identify genes differentially expressed between RTT samples and controls (Methods). We found that 244 genes (including 14 non-coding RNAs) showed significant expression changes between RTT samples and controls (Fig. 3, FDR < 0.05, Methods), of which 151 were down-regulated and 93 were up-regulated in Rett brain samples.
The use of two independent transcriptome analysis methods (RNA-seq and microarrays), with distinct chemistries, allowed us to assess the technical validation rate of the results on a genome-wide scale, rather than assessing a limited number of genes by a low-throughput method (such as qRT-PCR). Of the 244 differentially expressed genes, 174 were detected on the microarray platform, and 73 % of these were significantly differentially expressed based on the microarray data (FDR < 0.05, see Methods), with the same directionality of expression changes as identified in the RNA-seq data (Additional file 2: Table S2), indicating an appropriate technical validation rate [16, 17].
Given the limited availability of human brain tissue from Rett Synrome patients, we sought to use data from RTT mouse models to assess our results in an independent dataset. To this end, we compared our DE genes with those observed in RTT knock-out and transgenic mice . Ben-Shachar et al. identified genes with consistent expression changes in the cerebellum of MECP2 transgenic and knock-out mice (i.e. significant changes in both mouse models, with opposite directionality). Of the 270 genes identified as activated or repressed by MECP2 in mouse, and passing the detection threshold in our data, 13 genes were also significantly differentially expressed in human RTT samples (Table 1), with the same directionality as in the knock-out mice. These results indicated a significant overlap between the human and mouse differentially expressed genes (p = 1.29E-5, hypergeometric test; see Methods), despite the assessment of distinct brain regions. Remarkably, all three genes encoding subunits of the complement C1Q complex (C1QA, C1QB, C1QC) were downregulated in human RTT samples and in MECP2 knock-out mice, and upregulated in MECP2 transgenic mice, strongly suggesting that the expression of these genes is regulated by MECP2.
Gene ontology and pathway enrichment analyses (Methods) of down-regulated DE genes showed a significant overrepresentation of genes implicated in immune responses, and the complement cascade in particular (Fig. 4a), while up-regulated genes showed no significant functional enrichment. Notably, in addition to the C1Q complex genes, several other genes that belong to the REACTOME complement cascade pathway (C3, TGFBR2, CX3CR1 and TYROBP) also showed reduced expression in RTT brain (Fig. 4b).
Among the top 10 downregulated genes, two were non-coding RNAs specifically expressed in human brain and testis: RP11-178 F10.3 and RP11-122 F24.1 (-Additional file 3: Figure S1). Their functions are yet uncharacterized, and our data highlights them as relevant for further investigation.
Since a recent study reported that MECP2 primarily represses long genes in the human brain, we assessed the dependency of gene expression changes (log2 fold change) on gene length. We did not observe a significant correlation between gene expression changes and gene length (r2 = -0.02, Additional file 4: Figure S2). Although the top three longest genes were upregulated in RTT brain (Additional file 4: Figure S2), it is difficult to conclude that this observation reflects a general phenomenon of length-dependent transcriptional repression.
We were next interested to know how genes downregulated in RTT brains varied in expression in the normal human brain during fetal and postnatal development. Since the onset of RTT occurs around two years of age, after an initial period of normal development, we hypothesized that MECP2 transcriptional targets might play a particularly important role during this developmental time period.
To this end, we used transcriptome data from a recent study  of gene expression variation in over 200 human prefrontal cortex samples, including 71 samples from fetal and postnatal (<10 years) brains. We assessed the correlation of gene expression with age for all genes detected in our RNA-seq experiment and present in the Colantuoni et al.  dataset. As previously reported, the fetal to postnatal transition is associated with major gene expression changes, as most genes undergo either an increase or a decrease in expression during this developmental period. Interestingly, the genes differentially expressed in RTT brain showed primarily an increase in expression levels at the fetal to post-natal transition, a significant shift from the bi-modal distribution of all genes in the dataset (p < 0.05, Kolmogorov Smirnov test). By stratifying the genes into up- and down-regulated, we found that the increase in gene expression at the fetal-to-postnatal transition primarily reflected the behavior of downregulated genes (Fig. 5). Thus genes downregulated in RTT brain increase in expression in normal postnatal brain, suggesting that they may play a functional role during this developmental period.
Our study is the first RNA-seq analysis of post-mortem brain samples from RTT cases and age-, gender- and ethnicity-matched controls. By careful consideration of experimental co-variates such as cellular composition of brain tissue samples, and effective normalization of the data, we were able to detect several hundreds of differentially expressed genes, and observe significant convergence with gene expression changes previously reported in RTT mouse models.
Our data points out a transcriptional deficit of genes encoding components of the C1Q complement complex, as well as several other genes implicated in the complement pathway: C3, TGFBR2, CX3CR1 and TYROBP.
C1Q complex genes are expressed in the brain primarily in microglia, at low levels in resting state, and at higher levels following microglial activation . Although their role has been primarily characterized in neuronflammation, complement proteins play diverse roles in the brain, including a role in synaptic pruning . Interestingly, C1Q genes are expressed in neurons in response to TGFβ receptor (TGFBR2) activation during synaptic pruning . The expression of C1Q genes in response to TGFβ release from astrocytes leads to activation of the complement cascade, and microglial-dependent synaptic pruning .
An early transcriptome study of human RTT showed upregulation of glial transcripts and downregulation of neuron-specific mRNAs in post-mortem RTT brain . Whether those results were influenced by the tissue composition of brain samples remains unclear, but our study highlights the importance of considering tissue composition as a covariate in transcriptome analyses of brain tissue samples. Two other studies have previously investigated gene expression changes in RTT genome-wide. Gibson et al  compared gene expression in frontal and occipital cortex between six RTT cases and six controls. Due to significant age differences between the Rett and control group, 4–11 years old and 43–52 years old respectively, this study focused on inter-regional differences within each group. Our study on the other hand, was not adequately powered to assess gene expression differences between frontal and temporal cortex for RTT cases and controls separately, and thus larger datasets will be required to further evaluate brain-region specific gene expression differences in RTT.
Our gene expression data is consistent with one of three possible scenarios (a) a reduction of the total number of microglia in RTT, (b) normal number of resting microglia but reduced microglial activation, (c) reduced expression of C1Q genes in neurons in RTT brain.
It has been previously observed that microglia are depleted in the later developmental stages of RTT mouse brain, after an initial over-activation during early postnatal development . This observation is consistent with our data, showing reduced expression of C1Q complex genes in human RTT brain during adolescence. Whether microglial activation also occurs in the human RTT brain during early postnatal development, is yet unclear. However, reduced dendritic spine densities, which could result from excessive synaptic pruning during early postnatal development, is a feature of RTT documented both in mouse and human brains .
Previous mouse studies have highlighted a potential role of astrocytes and microglia in RTT [22, 24, 25]. Restoring MECP2 expression in astrocytes partially restored the neurological phenotypes in a RTT mouse model . In addition, human iPSC-derived astrocytes carrying RTT mutations adversely affect the morphology and function of co-cultured wild-type neurons . Moreover, it has been proposed that repopulating the brain with wild-type microglia by irradiation followed by bone-marrow transplantation could arrest the pathology in a mouse model of RTT . However, the latter result could not be replicated in independent studies , indicating that there is yet insufficient evidence to support the use of bone marrow transplantation in RTT patients. Taken together, the results of our study, the first genome-wide assessment of transcriptional changes in human RTT brain, support the notion that alterations in the expression level of C1Q complex genes may contribute to RTT pathogenesis.
Brain tissue samples
Post-mortem brain tissue from four female RTT patients and four age-, sex- and ethnicity- matched controls (Additional file 5: Table S3a) were obtained from the NICHD Brain and Tissue Bank for Developmental Disorders (http://medschool.umaryland.edu/btbank/). All work with these samples was completed in accordance with UNSW ethics requirements for work with human samples (project approval #HC13110). For each brain, approximately 500 mg of frozen tissue was obtained from frontal cortex (BA9 or PFC) and temporal cortex (BA21, BA22, or BA41/42). MECP2 mutation information was obtained from the NICHD brain bank and confirmed by Sanger sequencing (Additional file 5: Table S3a).
Total RNA was extracted from around 100 mg of brain tissue using the Qiagen miRNEasy kit following the manufacturer’s protocol, including an on-column DNase I treatment. An additional DNase I (New England Biolabs) treatment was carried out in solution, in order to eliminate any contaminating genomic DNA. Total RNA was eluted in 40 μl RNase-free water and stored at −80 °C. The concentration of RNA and double-stranded DNA (dsDNA) were measured using the Qubit® 2.0 Fluorometer (Life Technologies), and the RNA integrity number (RIN) of each sample was assessed using the Agilent 2100 Bioanalyzer (Agilent Technologies). The RNA and dsDNA concentrations and RIN values are listed in Additional file 5: Table S3b. Only samples with RIN ≥ 5.0 were used in downstream microarray and qRT-PCR analysis. Four RNA samples were excluded from further analysis due to low RNA quality: both RNA samples from the RTT case R4882 (R4882-F and R4882-T), aswell as the control samples C1078-T and C1541-F.
Six control RNA samples and six RTT RNA samples (Fig. 1) were analysed on Illumina HumanHT-12 v4 Expression BeadChip arrays. 500 ng of total RNA was processed according to the manufacturer’s protocol and hybridised at the UNSW Ramaciotti Centre for Genomics. All samples were run on the same chip in order to avoid batch effects. Microarray expression data was analysed using the R software (http://www.r-project.org), R package ruv13 and the Bioconductor package Lumi . Raw expression data was log2 transformed and normalized by quantile normalization. Microarray quality control criteria included high inter-array Pearson correlation coefficients, low variance of mean inter-array correlation and probe detection P values <0.05 in at least 30 % of the samples. Further, the ruv package was used to estimate the unwanted variation in the data. As recommended for this method , we obtained negative control genes by eliminating the top 5000 genes differentially expressed between RTT samples and controls in a first pass differential expression analysis. We applied the RUV-2 function with default parameters and k = 2, to obtain w coefficients, which were then used as covariates in the differential expression (DE) analysis. The choice k = 2 is the minimum k that effectively eliminates the effect of cellular composition on sample clustering. DE was carried out using limma , followed by Benjamini and Hochberg (BH) correction for multiple comparisons. Genes with FDR-corrected p-value < 0.05 were considered differentially expressed.
Strand-specific RNA-seq data was generated from six RNA samples (Fig. 1): a pool of equal volumes of frontal and temporal cortex RNA for each RTT brain (R1815-pool, R4516-pool, R4852-pool) and two of the control brains (C1078-pool, C1571-pool), as well as temporal cortex RNA from a third control sample (C1541-T). 1–5 μg of total RNA was depleted of ribosomal RNA using the Epicentre Ribo-zero kit, according to the manufacturer’s protocol. Library preparation using the Illumina TruSeq Stranded kit (http://www.illumina.com/products/truseq_stranded_total_rna_library_prep_kit.html) and sequencing on an Illumina HiSeq 2500 sequencer were carried out at the UNSW Ramaciotti Centre for Genomics. Libraries were barcoded and sequenced in one lane, thus avoiding lane effects, to obtain 100 bp paired-end reads. An average of 60 million reads were obtained for each sample (Additional file 6: Table S4). Sequencing reads were trimmed to eliminate adaptor sequences using Trimmomatic  aligned to the Human Genome (hg19) using the TopHat spliced aligner  with default parameters. Quality-control assessment of RNA-seq data was carried out using the RNAseqQC software . Gene-level expression quantification was carried out using the Rsubread featureCount function with the following parameters: useMetaFeatures = TRUE, allowMultiOverlap = FALSE, minMQS = 10, requireBothEndsMapped = TRUE, countChimericFragments = FALSE . Genes with expression level >0.5 FPKM in at least 3 of the 6 samples were retained for further analysis.
TMM (Trimmed Mean of M-values) normalised data was analysed for a first pass differential expression analysis using edgeR , in order to define negative control genes as described above for microarrays. The unwanted variation coefficients (w) were then estimated using the RUVg function  with default parameters and k = 1, and further used as covariate in the differential expression analysis. Differential expression analysis was carried out using edgeR with BH correction for multiple comparisons. Genes with FDR-corrected p-value < 0.05 were considered differentially expressed. The differential expression analysis results for RNA-seq and microarray data are provided in Additional file 7: Table S5 and Additional file 8: Table S6 respectively.
Tissue composition analysis
The cellular composition of brain tissue samples was assessed in silico, using the ssKL algorithm in CellMix R package , based on quantile normalised (but not log2 transformed) microarray expression data. Cell-type specific marker gene lists for neurons and astrocytes were obtained from Cahoy et al. . The estimated proportions of astrocytes and neurons in each sample are listed in Additional file 1: Table S1.
Gene ontology and pathway enrichment analysis was carried out using the GOseq package , with gene-length adjustment for RNA-seq data. Gene ontology terms with BH corrected p-value < 0.05 were considered significant. Pathway enrichment analysis was carried out using GSEA (http://www.broadinstitute.org/gsea/index.jsp) and the Reactome pathway database (http://www.reactome.org/). The gene ontology and pathway enrichment results are listed in Additional file 9: Table S7 and Additional file 10: Table S8.
Comparison with mouse data
The overlap between the 244 genes identified as differentially expressed in the present study (i.e. 1.3 % of the entire dataset) and genes differentially expressed in mouse models of RTT  was carried out using a hypergeometric test implemented in the phyper function in R. 270 genes had been identified as differentially expressed in MECP2 knockout and transgenic mice, and were quantified in our dataset, and 13 of those (i.e. 4.8 %) were also identified as dysregulated in human RTT brain, a statistically significant overlap (p = 1.29E-5).
BA, brain area; BH, Benjamini and Hochberg; CTD, C-terminal domain; DE, differentially expressed; ID, intervening domain; MBD, methyl-binding domain; PFC, prefrontal cortex; RTT, Rett syndrome; TRD, transcriptional repression domain.
Samaco R, Neul J. Complexities of Rett syndrome and MeCP2. The Journal of Neuroscience. 2011;31(22):7951–9.
Amir R, Van den Veyver I, Wan M, Tran C, Francke U, Zoghbi H. Rett syndrome is caused by mutations in X-linked MECP2, encoding methyl-CpG-binding protein 2. Nature Genetics. 1999;23(2):185–8.
Gabel HW, Kinde B, Stroud H, Gilbert CS, Harmin DA, Kastan NR, Hemberg M, Ebert DH, Greenberg ME. Disruption of DNA-methylation-dependent long gene repression in Rett syndrome. Nature. 2015;522(7554):89–93.
Skene PJ, Illingworth RS, Webb S, Kerr AR, James KD, Turner DJ, Andrews R, Bird AP. Neuronal MeCP2 is expressed at near histone-octamer levels and globally alters the chromatin state. Mol Cell. 2010;37(4):457–68.
Chahrour M, Jung S, Shaw C, Zhou X, Wong S, Qin J, Zoghbi H. MeCP2, a key contributor to neurological disease, activates and represses transcription. Science (New York, NY). 2008;320(5880):1224–9.
Ben-Shachar S, Chahrour M, Thaller C, Shaw CA, Zoghbi HY. Mouse models of MeCP2 disorders share gene expression changes in the cerebellum and hypothalamus. Hum Mol Genet. 2009;18(13):2431–42.
Nan X, Ng HH, Johnson CA, Laherty CD, Turner BM, Eisenman RN, Bird A. Transcriptional repression by the methyl-CpG-binding protein MeCP2 involves a histone deacetylase complex. Nature. 1998;393(6683):386–9.
Kimura H, Shiota K. Methyl-CpG-binding protein, MeCP2, is a target molecule for maintenance DNA methyltransferase, Dnmt1. J Biol Chem. 2003;278(7):4806–12.
Colantuoni C, Jeon O, Hyder K, Chenchik A, Khimani A, Narayanan V, Hoffman E, Kaufmann W, Naidu S, Pevsner J. Gene expression profiling in postmortem Rett Syndrome brain: differential gene expression and patient classification. Neurobiology of Disease. 2001;8(5):847–65.
Deng V, Matagne V, Banine F, Frerking M, Ohliger P, Budden S, Pevsner J, Dissen GA, Sherman LS, Ojeda SR. FXYD1 is an MeCP2 target gene overexpressed in the brains of Rett syndrome patients and Mecp2-null mice. Hum Mol Genet. 2007;16(6):640–50.
Gibson JH, Slobedman B, K NH, Williamson SL, Minchenko D, El-Osta A, Stern JL, Christodoulou J. Downstream targets of methyl CpG binding protein 2 and their abnormal expression in the frontal cortex of the human Rett syndrome brain. BMC Neurosci. 2010;11:53.
Gaujoux R, Seoighe C. Cell Mix: a comprehensive toolbox for gene expression deconvolution. Bioinformatics. 2013;29(17):2211–2.
Cahoy JD, Emery B, Kaushal A, Foo LC, Zamanian JL, Christopherson KS, Xing Y, Lubischer JL, Krieg PA, Krupenko SA, et al. A transcriptome database for astrocytes, neurons, and oligodendrocytes: a new resource for understanding brain development and function. The Journal of Neuroscience. 2008;28(1):264–78.
Gagnon-Bartsch JA, Speed TP. Using control genes to correct for unwanted variation in microarray data. Biostatistics. 2012;13(3):539–52.
Risso D, Ngai J, Speed TP, Dudoit S. Normalization of RNA-seq data using factor analysis of control genes or samples. Nat Biotechnol. 2014;32(9):896–902.
Wang C, Gong B, Bushel PR, Thierry-Mieg J, Thierry-Mieg D, Xu J, Fang H, Hong H, Shen J, Su Z, et al. The concordance between RNA-seq and microarray data depends on chemical treatment and transcript abundance. Nat Biotechnol. 2014;32(9):926–32.
Fang Z, Cui X. Design and validation issues in RNA-seq experiments. Briefings in Bioinformatics. 2011;12(3):280–7.
Colantuoni C, Lipska BK, Ye T, Hyde TM, Tao R, Leek JT, Colantuoni EA, Elkahloun AG, Herman MM, Weinberger DR, et al. Temporal dynamics and genetic control of transcription in the human prefrontal cortex. Nature. 2011;478(7370):519–23.
Kouser L, Madhukaran SP, Shastri A, Saraon A, Ferluga J, Al-Mozaini M, Kishore U. Emerging and Novel Functions of Complement Protein C1q. Front Immunol. 2015;6:317.
Veerhuis R, Nielsen HM, Tenner AJ. Complement in the brain. Mol Immunol. 2011;48(14):1592–603.
Bialas AR, Stevens B. TGF-beta signaling regulates neuronal C1q expression and developmental synaptic refinement. Nature Neuroscience. 2013;16(12):1773–82.
Cronk JC, Derecki NC, Ji E, Xu Y, Lampano AE, Smirnov I, Baker W, Norris GT, Marin I, Coddington N, et al. Methyl-CpG Binding Protein 2 Regulates Microglia and Macrophage Gene Expression in Response to Inflammatory Stimuli. Immunity. 2015;42(4):679–91.
Xu X, Miller EC, Pozzo-Miller L. Dendritic spine dysgenesis in Rett syndrome. Front Neuroanat. 2014;8:97.
Williams EC, Zhong X, Mohamed A, Li R, Liu Y, Dong Q, Ananiev GE, Mok JC, Lin BR, Lu J, et al. Mutant astrocytes differentiated from Rett syndrome patients-specific iPSCs have adverse effects on wild-type neurons. Hum Mol Genet. 2014;23(11):2968–80.
Lioy DT, Garg SK, Monaghan CE, Raber J, Foust KD, Kaspar BK, Hirrlinger PG, Kirchhoff F, Bissonnette JM, Ballas N, et al. A role for glia in the progression of Rett’s syndrome. Nature. 2011;475(7357):497–500.
Derecki NC, Cronk JC, Lu Z, Xu E, Abbott SB, Guyenet PG, Kipnis J. Wild-type microglia arrest pathology in a mouse model of Rett syndrome. Nature. 2012;484(7392):105–9.
Wang J, Wegener JE, Huang TW, Sripathy S, De Jesus-Cortes H, Xu P, Tran S, Knobbe W, Leko V, Britt J, et al. Wild-type microglia do not reverse pathology in mouse models of Rett syndrome. Nature. 2015;521(7552):E1–4.
Du P, Kibbe WA, Lin SM. lumi: a pipeline for processing Illumina microarray. Bioinformatics. 2008;24(13):1547–8.
Ritchie ME, Phipson B, Wu D, Hu Y, Law CW, Shi W, Smyth GK. limma powers differential expression analyses for RNA-sequencing and microarray studies. Nucleic Acids Res. 2015;43(7):e47.
Bolger AM, Lohse M, Usadel B. Trimmomatic: a flexible trimmer for Illumina sequence data. Bioinformatics. 2014;30(15):2114–20.
Trapnell C, Pachter L, Salzberg SL. TopHat: discovering splice junctions with RNA-Seq. Bioinformatics. 2009;25(9):1105–11.
DeLuca DS, Levin JZ, Sivachenko A, Fennell T, Nazaire MD, Williams C, Reich M, Winckler W, Getz G. RNA-SeQC: RNA-seq metrics for quality control and process optimization. Bioinformatics. 2012;28(11):1530–2.
Liao Y, Smyth GK, Shi W. featureCounts: an efficient general purpose program for assigning sequence reads to genomic features. Bioinformatics. 2014;30(7):923–30.
Robinson MD, McCarthy DJ, Smyth GK. edgeR: a Bioconductor package for differential expression analysis of digital gene expression data. Bioinformatics. 2010;26(1):139–40.
Young MD, Wakefield MJ, Smyth GK, Oshlack A. Gene ontology analysis for RNA-seq: accounting for selection bias. Genome Biol. 2010;11(2):R14.
This work was supported by a Ramaciotti establishment grant, an NHMRC project grant (APP1062510), and an Autism CRC Innovation grant to IV.
Availability of supporting data
RNA-seq and Microarray data are deposited in SRA (PRJNA302685) and GEO (GSE75303) respectively.
LN and ZF processed the brain and RNA samples. PL, IV, HA, AAA,TGA analyzed the data. RJE supervised data analysis. IV designed the study and supervised all aspects of the project. All authors have read and approved the manuscript.
The authors have no competing interests to declare.
All work with human brain samples was completed in accordance with UNSW ethics requirements for work with human samples (project approval #HC13110).
Tissue Composition. (XLS 6 kb)
List of Differentially Expressed Genes (RNA-seq). (XLS 50 kb)
Expression levels of RP11-178F10.3 and RP11-122F24.1 across human tissues. Boxplots display expression levels based on RNA-seq data from the GTEx Consortium. The plots were generated using the GTEx potal online tool: http://www.gtexportal.org. (JPG 280 kb)
Dependence between gene expression changes and gene length in Rett syndrome brain. Scatterplot displays gene length on the X-axis and fold-changes in Rett syndrome samples relative to control samples, as measured by microarrays, on the Y-axis, for all genes in the dataset. Genes identified as significantly differentially expressed are shown in red, while all other genes are shown in black. (JPG 96 kb)
Brain tissue and RNA samples. (XLS 10 kb)
RNA-seq data overview. (XLS 7 kb)
Differential Expression Results, RNA-seq. (XLS 2488 kb)
Differential Expression Results, Microarrays. (XLS 2356 kb)
Gene Ontology Enrichment Results for Downregulated Genes. (XLS 68 kb)
Pathway Enrichment Results for Downregulated Genes. (XLS 12 kb)
About this article
Cite this article
Lin, P., Nicholls, L., Assareh, H. et al. Transcriptome analysis of human brain tissue identifies reduced expression of complement complex C1Q Genes in Rett syndrome. BMC Genomics 17, 427 (2016). https://doi.org/10.1186/s12864-016-2746-7
- Rett Syndrome
- Transcriptome profiling