Disease gene characterization through large-scale co-expression analysis

Loss of Extreme Long-Range Enhancers in Human Neural Crest Drives a Craniofacial Disorder

Non-coding mutations at the far end of a large gene desert surrounding the SOX9 gene result in a human craniofacial disorder called Pierre Robin sequence (PRS). Leveraging a human stem cell differentiation model, we identify two clusters of enhancers within the PRS-associated region that regulate SOX9 expression during a restricted window of facial progenitor development at distances up to 1.45 Mb. Enhancers within the 1.45 Mb cluster exhibit highly synergistic activity that is dependent on the Coordinator motif. Using mouse models, we demonstrate that PRS phenotypic specificity arises from the convergence of two mechanisms: confinement of Sox9 dosage perturbation to developing facial structures through context-specific enhancer activity and heightened sensitivity of the lower jaw to Sox9 expression reduction. Overall, we characterize the longest-range human enhancers involved in congenital malformations, directly demonstrate that PRS is an enhanceropathy, and illustrate how small changes in gene expression can lead to morphological variation.

Comparative transcriptome analysis of the main beam and brow tine of sika deer antler provides insights into the molecular control of rapid antler growth

Background

Deer antlers have become a valuable model for biomedical research due to the capacities of regeneration and rapid growth. However, the molecular mechanism of rapid antler growth remains to be elucidated. The aim of the present study was to compare and explore the molecular control exerted by the main beam and brow tine during rapid antler growth.

Methods

The main beams and brow tines of sika deer antlers were collected from Chinese sika deer (Cervus nippon) at the rapid growth stage. Comparative transcriptome analysis was conducted using RNA-Seq technology. Differential expression was assessed using the DEGseq package. Functional Gene Ontology (GO) enrichment analysis was accomplished using a rigorous algorithm according to the GO Term Finder tool, and KEGG (Kyoto Encyclopedia of Genes and Genomes) pathway enrichment analysis was accomplished with the R function phyper, followed by the hypergeometric test and Bonferroni correction. Quantitative real-time polymerase chain reaction (qRT-PCR) was carried out to verify the RNA levels for differentially expressed mRNAs.

Results

The expression levels of 16 differentially expressed genes (DEGs) involved in chondrogenesis and cartilage development were identified as significantly upregulated in the main beams, including transcription factor SOX-9 (Sox9), collagen alpha-1(II) chain (Col2a1), aggrecan core protein (Acan), etc. However, the expression levels of 17 DEGs involved in endochondral ossification and bone formation were identified as significantly upregulated in the brow tines, including collagen alpha-1(X) chain (Col10a1), osteopontin (Spp1) and bone sialoprotein 2 (Ibsp), etc.

Conclusion

These results suggest that the antler main beam has stronger growth capacity involved in chondrogenesis and cartilage development compared to the brow tine during rapid antler growth, which is mainly achieved through regulation of Sox9 and its target genes, whereas the antler brow tine has stronger capacities of endochondral bone formation and resorption compared to the main beam during rapid antler growth, which is mainly achieved through the genes involved in regulating osteoblast and osteoclast activities. Thus, the current research has deeply expanded our understanding of the intrinsic molecular regulation displayed by the main beam and brow tine during rapid antler growth.

Meningioma transcription factors link cell lineage with systemic metabolic cues

Background

Tumor cells recapitulate cell-lineage transcriptional programs that are characteristic of normal tissues from which they arise. It is unclear why such lineage programs are fatefully maintained in tumors and if they contribute to cell proliferation and viability.

Methods

Here, we used the most common brain tumor, meningioma, which is strongly associated with female sex and high body mass index (BMI), as a model system to address these questions. We screened expression profiling data to identify the transcription factor (TF) genes which are highly enriched in meningioma, and characterized the expression pattern of those TFs and downstream genes in clinical meningioma samples as well as normal brain tissues. Meningioma patient-derived cell lines (PDCLs) were used for further validation and characterization.

Results

We identified 8 TFs highly enriched in meningioma. Expression of these TFs, which included sine oculis homeobox 1 (SIX1), readily distinguished meningiomas from other primary brain tumors and was maintained in PDCLs and even in pulmonary meningothelial nodules. In meningioma PDCLs, SIX1 and its coactivator eyes absent 2 (EYA2) supported the expression of the leptin receptor (LEPR), the cell-surface receptor for leptin (LEP), the adipose-specific hormone that is high in women and in individuals with high BMI. Notably, these transcriptional regulatory factors, LEPR and LEP, both contributed to support meningioma PDCLs proliferation and survival, elucidating a survival dependency on both a core transcriptional program and a metabolic cell-surface receptor.

Conclusions

These findings provide one rationale for why lineage TF expression is maintained in meningioma and for the epidemiological association of female sex and obesity with meningioma risk.

Gene multifunctionality scoring using gene ontology

Multifunctional genes are important genes because of their essential roles in human cells. Studying and analyzing multifunctional genes can help understand disease mechanisms and drug discovery. We propose a computational method for scoring gene multifunctionality based on functional annotations of the target gene from the Gene Ontology. The method is based on identifying pairs of GO annotations that represent semantically different biological functions and any gene annotated with two annotations from one pair is considered multifunctional. The proposed method can be employed to identify multifunctional genes in the entire human genome using solely the GO annotations. We evaluated the proposed method in scoring multifunctionality of all human genes using four criteria: gene-disease associations; protein-protein interactions; gene studies with PubMed publications; and published known multifunctional gene sets. The evaluation results confirm the validity and reliability of the proposed method for identifying multifunctional human genes. The results across all four evaluation criteria were statistically significant in determining multifunctionality. For example, the method confirmed that multifunctional genes tend to be associated with diseases more than other genes, with significance [Formula: see text]. Moreover, consistent with all previous studies, proteins encoded by multifunctional genes, based on our method, are involved in protein-protein interactions significantly more ([Formula: see text]) than other proteins.

Lack of GAS2L2 Causes PCD by Impairing Cilia Orientation and Mucociliary Clearance

Primary ciliary dyskinesia (PCD) is a genetic disorder in which impaired ciliary function leads to chronic airway disease. Exome sequencing of a PCD subject identified an apparent homozygous frameshift variant, c.887_890delTAAG (p.Val296Glyfs^∗13), in exon 5; this frameshift introduces a stop codon in amino acid 308 of the growth arrest-specific protein 2-like 2 (GAS2L2). Further genetic screening of unrelated PCD subjects identified a second proband with a compound heterozygous variant carrying the identical frameshift variant and a large deletion (c.867_^∗343+1207del; p.?) starting in exon 5. Both individuals had clinical features of PCD but normal ciliary axoneme structure. In this research, using human nasal cells, mouse models, and X.laevis embryos, we show that GAS2L2 is abundant at the apical surface of ciliated cells, where it localizes with basal bodies, basal feet, rootlets, and actin filaments. Cultured GAS2L2-deficient nasal epithelial cells from one of the affected individuals showed defects in ciliary orientation and had an asynchronous and hyperkinetic (GAS2L2-deficient = 19.8 Hz versus control = 15.8 Hz) ciliary-beat pattern. These results were recapitulated in Gas2l2^-/- mouse tracheal epithelial cell (mTEC) cultures and in X. laevis embryos treated with Gas2l2 morpholinos. In mice, the absence of Gas2l2 caused neonatal death, and the conditional deletion of Gas2l2 impaired mucociliary clearance (MCC) and led to mucus accumulation. These results show that a pathogenic variant in GAS2L2 causes a genetic defect in ciliary orientation and impairs MCC and results in PCD.

Genes uniquely expressed in human growth plate chondrocytes uncover a distinct regulatory network

Background

Chondrogenesis is the earliest stage of skeletal development and is a highly dynamic process, integrating the activities and functions of transcription factors, cell signaling molecules and extracellular matrix proteins. The molecular mechanisms underlying chondrogenesis have been extensively studied and multiple key regulators of this process have been identified. However, a genome-wide overview of the gene regulatory network in chondrogenesis has not been achieved.

Results

In this study, employing RNA sequencing, we identified 332 protein coding genes and 34 long non-coding RNA (lncRNA) genes that are highly selectively expressed in human fetal growth plate chondrocytes. Among the protein coding genes, 32 genes were associated with 62 distinct human skeletal disorders and 153 genes were associated with skeletal defects in knockout mice, confirming their essential roles in skeletal formation. These gene products formed a comprehensive physical interaction network and participated in multiple cellular processes regulating skeletal development. The data also revealed 34 transcription factors and 11,334 distal enhancers that were uniquely active in chondrocytes, functioning as transcriptional regulators for the cartilage-selective genes.

Conclusions

Our findings revealed a complex gene regulatory network controlling skeletal development whereby transcription factors, enhancers and lncRNAs participate in chondrogenesis by transcriptional regulation of key genes. Additionally, the cartilage-selective genes represent candidate genes for unsolved human skeletal disorders.

Electronic supplementary material

The online version of this article (10.1186/s12864-017-4378-y) contains supplementary material, which is available to authorized users.

Cortical gene expression: prognostic value for seizure outcome following temporal lobectomy and amygdalohippocampectomy

Whole genome analyses were performed to test the hypothesis that temporal cortical gene expression differs between epilepsy patients rendered seizure-free versus non-seizure-free following anterior temporal lobectomy with amygdalohippocampectomy (ATL/AH). Twenty four patients underwent ATL/AH to treat medically intractable seizures of temporal lobe origin (mean age 35.5 years, mean follow-up 42.2 months); they were then dichotomized into seizure-free and non-seizure-free groups. Tissue RNA was isolated from the lateral temporal cortex and gene expression analysis was performed. Whole genome data were analyzed for prognostic value for seizure-free outcome following ATL/AH by logistic regression. Genes that could distinguish seizure outcome groups were identified based on providing an accuracy of >0.90 judging by area under the receiver operating characteristic curve, AUC, with a P value of the slope coefficient of <0.05. Four genes and seven RNA probes were with prognostic value for post-operative seizure-free outcome. Gene expression associated with seizure-free outcome included relative down-regulation of zinc finger protein 852 (ZNF852), CUB domain-containing protein 2 (CDCP2), proline-rich transmembrane protein 1 (PRRT1), hypothetical LOC440200 (FLJ41170), RNA probe 8047763, RNA probe 8126238, RNA probe 8113489, RNA probe 8092883, RNA probe 7935228, RNA probe 806293, and RNA probe 8104131. This study describes the predictive value of temporal cortical gene expression for seizure-free outcome after ATL/AH. Four genes and seven RNA probes were found to predict post-operative seizure-free outcome. Future prospective investigation of these genes and probes in human brain tissue and blood could establish new biomarkers predictive of seizure outcome following ATL/AH.

HSP47 and FKBP65 cooperate in the synthesis of type I procollagen

Osteogenesis imperfecta (OI) is a genetic disorder that results in low bone mineral density and brittle bones. Most cases result from dominant mutations in the type I procollagen genes, but mutations in a growing number of genes have been identified that produce autosomal recessive forms of the disease. Among these include mutations in the genes SERPINH1 and FKBP10, which encode the type I procollagen chaperones HSP47 and FKBP65, respectively, and predominantly produce a moderately severe form of OI. Little is known about the biochemical consequences of the mutations and how they produce OI. We have identified a new OI mutation in SERPINH1 that results in destabilization and mislocalization of HSP47 and secondarily has similar effects on FKBP65. We found evidence that HSP47 and FKBP65 act cooperatively during posttranslational maturation of type I procollagen and that FKBP65 and HSP47 but fail to properly interact in mutant HSP47 cells. These results thus reveal a common cellular pathway in cases of OI caused by HSP47 and FKBP65 deficiency.

Orchestration of neurodevelopmental programs by RBFOX1: implications for autism spectrum disorder

Neurodevelopmental and neuropsychiatric disorders result from complex interactions between critical genetic factors and as-yet-unknown environmental components. To gain clinical insight, it is critical to develop a comprehensive understanding of these genetic components. RBFOX1, an RNA splicing factor, regulates expression of large genetic networks during early neuronal development, and haploinsufficiency causes severe neurodevelopmental phenotypes including autism spectrum disorder (ASD), intellectual disability, and epilepsy. Genomic testing in individuals and large patient cohorts has identified phenotypically similar cases possessing copy number variations in RBFOX1, implicating the gene as an important cause of neurodevelopmental disease. However, a significant proportion of the observed structural variation is inherited from phenotypically normal individuals, raising questions regarding overall pathogenicity of variation at the RBFOX1 locus. In this chapter, we discuss the molecular, cellular, and clinical evidence supporting the role of RBFOX1 in neurodevelopment and present a comprehensive model for the contribution of structural variation in RBFOX1 to ASD.

Harnessing gene expression networks to prioritize candidate epileptic encephalopathy genes

We apply a novel gene expression network analysis to a cohort of 182 recently reported candidate Epileptic Encephalopathy genes to identify those most likely to be true Epileptic Encephalopathy genes. These candidate genes were identified as having single variants of likely pathogenic significance discovered in a large-scale massively parallel sequencing study. Candidate Epileptic Encephalopathy genes were prioritized according to their co-expression with 29 known Epileptic Encephalopathy genes. We utilized developing brain and adult brain gene expression data from the Allen Human Brain Atlas (AHBA) and compared this to data from Celsius: a large, heterogeneous gene expression data warehouse. We show replicable prioritization results using these three independent gene expression resources, two of which are brain-specific, with small sample size, and the third derived from a heterogeneous collection of tissues with large sample size. Of the nineteen genes that we predicted with the highest likelihood to be true Epileptic Encephalopathy genes, two (GNAO1 and GRIN2B) have recently been independently reported and confirmed. We compare our results to those produced by an established in silico prioritization approach called Endeavour, and finally present gene expression networks for the known and candidate Epileptic Encephalopathy genes. This highlights sub-networks of gene expression, particularly in the network derived from the adult AHBA gene expression dataset. These networks give clues to the likely biological interactions between Epileptic Encephalopathy genes, potentially highlighting underlying mechanisms and avenues for therapeutic targets.

Mutations in SPAG1 cause primary ciliary dyskinesia associated with defective outer and inner dynein arms

Primary ciliary dyskinesia (PCD) is a genetically heterogeneous, autosomal-recessive disorder, characterized by oto-sino-pulmonary disease and situs abnormalities. PCD-causing mutations have been identified in 20 genes, but collectively they account for only ∼65% of all PCDs. To identify mutations in additional genes that cause PCD, we performed exome sequencing on three unrelated probands with ciliary outer and inner dynein arm (ODA+IDA) defects. Mutations in SPAG1 were identified in one family with three affected siblings. Further screening of SPAG1 in 98 unrelated affected individuals (62 with ODA+IDA defects, 35 with ODA defects, 1 without available ciliary ultrastructure) revealed biallelic loss-of-function mutations in 11 additional individuals (including one sib-pair). All 14 affected individuals with SPAG1 mutations had a characteristic PCD phenotype, including 8 with situs abnormalities. Additionally, all individuals with mutations who had defined ciliary ultrastructure had ODA+IDA defects. SPAG1 was present in human airway epithelial cell lysates but was not present in isolated axonemes, and immunofluorescence staining showed an absence of ODA and IDA proteins in cilia from an affected individual, thus indicating that SPAG1 probably plays a role in the cytoplasmic assembly and/or trafficking of the axonemal dynein arms. Zebrafish morpholino studies of spag1 produced cilia-related phenotypes previously reported for PCD-causing mutations in genes encoding cytoplasmic proteins. Together, these results demonstrate that mutations in SPAG1 cause PCD with ciliary ODA+IDA defects and that exome sequencing is useful to identify genetic causes of heterogeneous recessive disorders.

Exome sequencing identifies mutations in CCDC114 as a cause of primary ciliary dyskinesia

Primary ciliary dyskinesia (PCD) is a genetically heterogeneous, autosomal-recessive disorder, characterized by oto-sino-pulmonary disease and situs abnormalities. PCD-causing mutations have been identified in 14 genes, but they collectively account for only ~60% of all PCD. To identify mutations that cause PCD, we performed exome sequencing on six unrelated probands with ciliary outer dynein arm (ODA) defects. Mutations in CCDC114, an ortholog of the Chlamydomonas reinhardtii motility gene DCC2, were identified in a family with two affected siblings. Sanger sequencing of 67 additional individuals with PCD with ODA defects from 58 families revealed CCDC114 mutations in 4 individuals in 3 families. All 6 individuals with CCDC114 mutations had characteristic oto-sino-pulmonary disease, but none had situs abnormalities. In the remaining 5 individuals with PCD who underwent exome sequencing, we identified mutations in two genes (DNAI2, DNAH5) known to cause PCD, including an Ashkenazi Jewish founder mutation in DNAI2. These results revealed that mutations in CCDC114 are a cause of ciliary dysmotility and PCD and further demonstrate the utility of exome sequencing to identify genetic causes in heterogeneous recessive disorders.

HSF1 drives a transcriptional program distinct from heat shock to support highly malignant human cancers

Heat-Shock Factor 1 (HSF1), master regulator of the heat-shock response, facilitates malignant transformation, cancer cell survival, and proliferation in model systems. The common assumption is that these effects are mediated through regulation of heat-shock protein (HSP) expression. However, the transcriptional network that HSF1 coordinates directly in malignancy and its relationship to the heat-shock response have never been defined. By comparing cells with high and low malignant potential alongside their nontransformed counterparts, we identify an HSF1-regulated transcriptional program specific to highly malignant cells and distinct from heat shock. Cancer-specific genes in this program support oncogenic processes: cell-cycle regulation, signaling, metabolism, adhesion and translation. HSP genes are integral to this program, however, many are uniquely regulated in malignancy. This HSF1 cancer program is active in breast, colon and lung tumors isolated directly from human patients and is strongly associated with metastasis and death. Thus, HSF1 rewires the transcriptome in tumorigenesis, with prognostic and therapeutic implications.

Exome sequencing identifies CCDC8 mutations in 3-M syndrome, suggesting that CCDC8 contributes in a pathway with CUL7 and OBSL1 to control human growth

3-M syndrome, a primordial growth disorder, is associated with mutations in CUL7 and OBSL1. Exome sequencing now identifies mutations in CCDC8 as a cause of 3-M syndrome. CCDC8 is a widely expressed gene that is transcriptionally associated to CUL7 and OBSL1, and coimmunoprecipitation indicates a physical interaction between CCDC8 and OBSL1 but not CUL7. We propose that CUL7, OBSL1, and CCDC8 are members of a pathway controlling mammalian growth.

BioGraph: unsupervised biomedical knowledge discovery via automated hypothesis generation

We present BioGraph, a data integration and data mining platform for the exploration and discovery of biomedical information. The platform offers prioritizations of putative disease genes, supported by functional hypotheses. We show that BioGraph can retrospectively confirm recently discovered disease genes and identify potential susceptibility genes, outperforming existing technologies, without requiring prior domain knowledge. Additionally, BioGraph allows for generic biomedical applications beyond gene discovery. BioGraph is accessible at http://www.biograph.be.

BioGraph: unsupervised biomedical knowledge discovery via automated hypothesis generation

High-throughput technologies for DNA sequencing and for analyses of transcriptomes, proteomes and metabolomes have provided the foundations for deciphering the structure, variation and function of the human genome and relating them to health and disease states. The increased efficiency of DNA sequencing opens up the possibility of analyzing a large number of individual genomes and transcriptomes, and complete reference proteomes and metabolomes are within reach using powerful analytical techniques based on chromatography, mass spectrometry and nuclear magnetic resonance. Computational and mathematical tools have enabled the development of systems approaches for deciphering the functional and regulatory networks underlying the behavior of complex biological systems. Further conceptual and methodological developments of these tools are needed for the integration of various data types across the multiple levels of organization and time frames that are characteristic of human development, physiology and disease. Medical genomics has attempted to overcome the initial limitations of genome-wide association studies and has identified a limited number of susceptibility loci for many complex and common diseases. Iterative systems approaches are starting to provide deeper insights into the mechanisms of human diseases, and to facilitate the development of better diagnostic and prognostic biomarkers for cancer and many other diseases. Systems approaches will transform the way drugs are developed through academy-industry partnerships that will target multiple components of networks and pathways perturbed in diseases. They will enable medicine to become predictive, personalized, preventive and participatory, and, in the process, concepts and methods from Western and oriental cultures can be combined. We recommend that systems medicine should be developed through an international network of systems biology and medicine centers dedicated to inter-disciplinary training and education, to help reduce the gap in healthcare between developed and developing countries.

Genome-wide association and linkage identify modifier loci of lung disease severity in cystic fibrosis at 11p13 and 20q13.2

A combined genome-wide association and linkage study was used to identify loci causing variation in cystic fibrosis lung disease severity. We identified a significant association (P = 3.34 × 10(-8)) near EHF and APIP (chr11p13) in p.Phe508del homozygotes (n = 1,978). The association replicated in p.Phe508del homozygotes (P = 0.006) from a separate family based study (n = 557), with P = 1.49 × 10(-9) for the three-study joint meta-analysis. Linkage analysis of 486 sibling pairs from the family based study identified a significant quantitative trait locus on chromosome 20q13.2 (log(10) odds = 5.03). Our findings provide insight into the causes of variation in lung disease severity in cystic fibrosis and suggest new therapeutic targets for this life-limiting disorder.

Loss of CHSY1, a secreted FRINGE enzyme, causes syndromic brachydactyly in humans via increased NOTCH signaling

We delineated a syndromic recessive preaxial brachydactyly with partial duplication of proximal phalanges to 16.8 Mb over 4 chromosomes. High-throughput sequencing of all 177 candidate genes detected a truncating frameshift mutation in the gene CHSY1 encoding a chondroitin synthase with a Fringe domain. CHSY1 was secreted from patients' fibroblasts and was required for synthesis of chondroitin sulfate moieties. Noticeably, its absence triggered massive production of JAG1 and subsequent NOTCH activation, which could only be reversed with a wild-type but not a Fringe catalytically dead CHSY1 construct. In vitro, depletion of CHSY1 by RNAi knockdown resulted in enhanced osteogenesis in fetal osteoblasts and remarkable upregulation of JAG2 in glioblastoma cells. In vivo, chsy1 knockdown in zebrafish embryos partially phenocopied the human disorder; it increased NOTCH output and impaired skeletal, pectoral-fin, and retinal development. We conclude that CHSY1 is a secreted FRINGE enzyme required for adjustment of NOTCH signaling throughout human and fish embryogenesis and particularly during limb patterning.

Fibrochondrogenesis results from mutations in the COL11A1 type XI collagen gene

Fibrochondrogenesis is a severe, autosomal-recessive, short-limbed skeletal dysplasia. In a single case of fibrochondrogenesis, whole-genome SNP genotyping identified unknown ancestral consanguinity by detecting three autozygous regions. Because of the predominantly skeletal nature of the phenotype, the 389 genes localized to the autozygous intervals were prioritized for mutation analysis by correlation of their expression with known cartilage-selective genes via the UCLA Gene Expression Tool, UGET. The gene encoding the α1 chain of type XI collagen (COL11A1) was the only cartilage-selective gene among the three candidate intervals. Sequence analysis of COL11A1 in two genetically independent fibrochondrogenesis cases demonstrated that each was a compound heterozygote for a loss-of-function mutation on one allele and a mutation predicting substitution for a conserved triple-helical glycine residue on the other. The parents who were carriers of missense mutations had myopia. Early-onset hearing loss was noted in both parents who carried a loss-of-function allele, suggesting COL11A1 as a locus for mild, dominantly inherited hearing loss. These findings identify COL11A1 as a locus for fibrochondrogenesis and indicate that there might be phenotypic manifestations among carriers.

Altered gene synchrony suggests a combined hormone-mediated dysregulated state in major depression

Coordinated gene transcript levels across tissues (denoted “gene synchrony”) reflect converging influences of genetic, biochemical and environmental factors; hence they are informative of the biological state of an individual. So could brain gene synchrony also integrate the multiple factors engaged in neuropsychiatric disorders and reveal underlying pathologies? Using bootstrapped Pearson correlation for transcript levels for the same genes across distinct brain areas, we report robust gene transcript synchrony between the amygdala and cingulate cortex in the human postmortem brain of normal control subjects (n = 14; Control/Permutated data, p<0.000001). Coordinated expression was confirmed across distinct prefrontal cortex areas in a separate cohort (n = 19 subjects) and affected different gene sets, potentially reflecting regional network- and function-dependent transcriptional programs. Genewise regional transcript coordination was independent of age-related changes and array technical parameters. Robust shifts in amygdala-cingulate gene synchrony were observed in subjects with major depressive disorder (MDD, denoted here “depression”) (n = 14; MDD/Permutated data, p<0.000001), significantly affecting between 100 and 250 individual genes (10–30% false discovery rate). Biological networks and signal transduction pathways corresponding to the identified gene set suggested putative dysregulated functions for several hormone-type factors previously implicated in depression (insulin, interleukin-1, thyroid hormone, estradiol and glucocorticoids; p<0.01 for association with depression-related networks). In summary, we showed that coordinated gene expression across brain areas may represent a novel molecular probe for brain structure/function that is sensitive to disease condition, suggesting the presence of a distinct and integrated hormone-mediated corticolimbic homeostatic, although maladaptive and pathological, state in major depression.