Human glucocorticoid receptor isoform beta: recent understanding of its potential implications in physiology and pathophysiology.

The human glucocorticoid receptor (GR) gene expresses two splicing isoforms alpha and beta through alternative use of specific exons 9alpha and 9beta. In contrast to the classic receptor GRalpha, which mediates most of the known actions of glucocorticoids, the functions of GRbeta have been largely unexplored. Owing to newly developed methods, for example microarrays and the jellyfish fluorescence proteins, we and others have recently revealed novel functions of GRbeta. Indeed, this enigmatic GR isoform influences positively and negatively the transcriptional activity of large subsets of genes, most of which are not responsive to glucocorticoids, in addition to its well-known dominant negative effect against GRalpha-mediated transcriptional activity. A recent report suggested that the "ligand-binding domain" of GRbeta is active, forming a functional ligand-binding pocket associated with the synthetic compound RU 486. In this review, we discuss the functions of GRbeta, its mechanisms of action, and its pathologic implications.

Glucocorticoid receptor (GR) beta has intrinsic, GRalpha-independent transcriptional activity.

The human glucocorticoid receptor (GR) gene produces C-terminal GRbeta and GRalpha isoforms through alternative use of specific exons 9beta and alpha, respectively. We explored the transcriptional activity of GRbeta on endogenous genes by developing HeLa cells stably expressing EGFP-GRbeta or EGFP. Microarray analyses revealed that GRbeta had intrinsic gene-specific transcriptional activity, regulating mRNA expression of a large number of genes negatively or positively. Majority of GRbeta-responsive genes was distinct from those modulated by GRalpha, while GRbeta and GRalpha mutually modulated each other's transcriptional activity in a subpopulation of genes. We did not observe in HCT116 cells nuclear translocation of GRbeta and activation of this receptor by RU 486, a synthetic steroid previously reported to bind GRbeta and to induce nuclear translocation. Our results indicate that GRbeta has intrinsic, GRalpha-independent, gene-specific transcriptional activity, in addition to its previously reported dominant negative effect on GRalpha-induced transactivation of GRE-driven promoters.

Molecular mechanism underlying differential apoptosis between human melanoma cell lines UACC903 and UACC903(+6) revealed by mitochondria-focused cDNA microarrays.

Human malignant melanoma cell line UACC903 is resistant to apoptosis while chromosome 6-mediated suppressed cell line UACC903(+6) is sensitive. Here, we describe identification of differential molecular pathways underlying this difference. Using our recently developed mitochondria-focused cDNA microarrays, we identified 154 differentially expressed genes including proapoptotic (BAK1 [6p21.3], BCAP31, BNIP1, CASP3, CASP6, FAS, FDX1, FDXR, TNFSF10 and VDAC1) and antiapoptotic (BCL2L1, CLN3 and MCL1) genes. Expression of these pro- and anti-apoptotic genes was higher in UACC903(+6) than in UACC903 before UV treatment and was altered after UV treatment. qRT-PCR and Western blots validated microarray results. Our bioinformatic analysis mapped these genes to differential molecular pathways that predict resistance and sensitivity of UACC903 and UACC903(+6) to apoptosis respectively. The pathways were functionally confirmed by the FAS ligand-induced cell death and by siRNA knockdown of BAK1 protein. These results demonstrated the differential molecular pathways underlying survival and apoptosis of UACC903 and UACC903(+6) cell lines.

Mitochondria as key components of the stress response.

The exquisitely orchestrated adaptive response to stressors that challenge the homeostasis of the cell and organism involves important changes in mitochondrial function. A complex signaling network enables mitochondria to sense internal milieu or environmental changes and to adjust their bioenergetic, thermogenic, oxidative and/or apoptotic responses accordingly, aiming at re-establishment of homeostasis. Mitochondrial dysfunction is increasingly recognized as a key component in both acute and chronic allostatic states, although the extent of its role in the pathogenesis of such conditions remains controversial. Genetic and environmental factors that determine mitochondrial function might contribute to the significant variation of the stress response. Understanding the often reciprocal interplay between stress mediators and mitochondrial function is likely to help identify potential therapeutic targets for many stress and mitochondria-related pathologies.

Third-generation human mitochondria-focused cDNA microarray and its bioinformatic tools for analysis of gene expression.

To facilitate profiling mitochondrial transcriptomes, we developed a third-generation human mitochondria-focused cDNA microarray (hMitChip3) and its bioinformatic tools. hMitChip3 consists of the 37 mitochondrial DNA-encoded genes, 1098 nuclear DNA-encoded and mitochondria-related genes, and 225 controls, each in triplicate. The bioinformatic tools included data analysis procedures and customized database for interpretation of results. The database associated 645 molecular functions with 946 hMitChip3 genes, 612 biological processes with 930 genes, 172 cellular components with 869 genes, 107 biological chemistry pathways with 476 genes, 23 reactome events with 227 genes, 320 genetic disorders with 237 genes, and 87 drugs targets with 55 genes. To test these tools, hMitChip3 was used to compare expression profiles between human melanoma cell lines UACC903 (rapidly dividing) and UACC903(+6) (slowly dividing). Our results demonstrated internal gene-set consistency (correlation R > or = 0.980 +/- 0.005) and interexperimental reproducibility (R > or = 0.931 +/- 0.013). Expression patterns of 16 genes, involved in DNA, RNA, or protein biosyntheses in mitochondrial and other organelles, were consistent with the proliferation rates of both cell lines, and the patterns of 6 tested genes were verified by quantitative reverse transcription PCR (RT-PCR). Thus, hMitChip3 and its bioinformatics software provide an integrated tool for profiling mitochondria-focused gene expression.

Signature patterns revealed by microarray analyses of mice infected with influenza virus A and Streptococcus pneumoniae.

We used cDNA microarrays to identify differentially expressed genes in mice in response to infections with influenza virus A/PR/8/34 (H1N1) and Streptococcus pneumoniae. Expression microarray analysis showed up-regulation and down-regulation of many genes involved in the defense, inflammatory response and intracellular signaling pathways including chemokine, apoptosis, MAPK, Notch, Jak-STAT, T-cell receptor and complement and coagulation cascades. We have revealed signature patterns of gene expression in mice infected with two different classes of pathogens: influenza virus A and S. pneumoniae. Quantitative real-time RT-PCR results confirmed microarray results for most of the genes tested. These studies document clear differences in gene expression profiles between mice infected with influenza virus A and S. pneumoniae. Identification of genes that are differentially expressed after respiratory infections can provide insights into the mechanisms by which the host interacts with different pathogens, useful information about stage of diseases and selection of suitable targets for early diagnosis and treatments. The advantage of this novel approach is that the detection of pathogens is based on the differences in host gene expression profiles in response to different pathogens instead of detecting pathogens directly.

Monoamine oxidase-A is a major target gene for glucocorticoids in human skeletal muscle cells.

Skeletal myopathy is a common complication of endogenous and exogenous glucocorticoid excess, yet its pathogenetic mechanisms remain unclear. There is accumulating evidence that mitochondrial dysfunction and oxidative stress are involved in this process. To explore the glucocorticoid-induced transcriptional adaptations that may affect mitochondrial function in skeletal muscle, we studied gene expression profiles in dexamethasone-treated primary human skeletal myocytes using a cDNA microarray, which contains 501 mitochondria-related genes. We found that monoamine oxidase A (MAO-A) was the most significantly up-regulated gene. MAO-A is the primary enzyme metabolizing catecholamines and dietary amines, and its role in skeletal muscle remains largely unexplored. Dexamethasone induced dose- and time-dependent increases of MAO-A gene and protein expression, while its effects on MAO-B were minimal. Both the glucocorticoid receptor (GR) and the Sp1 transcription factor were required for dexamethasone-induced MAO-A mRNA expression, as blockade of the GR with RU 486 or ablation of Sp1 binding with mithramycin abrogated MAO-A mRNA induction. The observed dexamethasone effect was biologically functional, as this steroid significantly increased MAO-mediated hydrogen peroxide production. We suggest that MAO-A-mediated oxidative stress can lead to cell damage, representing a novel pathogenetic mechanism for glucocorticoid-induced myopathy and a potential target for therapeutic intervention.

Identification of differentially expressed genes in mouse spermatogenesis.

Complementary DNA microarray and quantitative polymerase chain reaction were used as tools for discovering genes that are differentially expressed in the mouse under normal physiological conditions at distinctive stages of male germ cell development, that is, type A spermatogonia, pachytene spermatocytes, and round spermatids. By using this strategy, we identified a set of genes exhibiting differential expression patterns in spermatogenesis, suggesting that specific functions of the encoded products occurred during the developmental process. Among them were several genes previously not known to be active in testis, which signified undiscovered functional roles of these genes during spermatogenesis. Many of the genes identified were not previously characterized. This study highlights new targets for manipulation to unravel the molecular mechanism of spermatogenesis.

Whole-Genome Sequences of Two Clinical Isolates of Extensively Drug-Resistant Mycobacterium tuberculosis from Zunyi, China.

Before 2013, 92 countries reported extensively drug-resistant Mycobacterium tuberculosis cases to the WHO. Here, we announce the genome sequences of two clinical isolates of extensively drug-resistant tuberculosis (XDR-TB) from Zunyi, China. The genome sequences are composed of 4,411,507 bp and 4,411,515 bp with 2,210 and 2,071 variants, respectively, when compared to the H37Rv genome.

Mitochondrial Gene Expression Profiles and Metabolic Pathways in the Amygdala Associated with Exaggerated Fear in an Animal Model of PTSD.

The metabolic mechanisms underlying the development of exaggerated fear in post-traumatic stress disorder (PTSD) are not well defined. In the present study, alteration in the expression of genes associated with mitochondrial function in the amygdala of an animal model of PTSD was determined. Amygdala tissue samples were excised from 10 non-stressed control rats and 10 stressed rats, 14 days post-stress treatment. Total RNA was isolated, cDNA was synthesized, and gene expression levels were determined using a cDNA microarray. During the development of the exaggerated fear associated with PTSD, 48 genes were found to be significantly upregulated and 37 were significantly downregulated in the amygdala complex based on stringent criteria (p < 0.01). Ingenuity pathway analysis revealed up- or downregulation in the amygdala complex of four signaling networks - one associated with inflammatory and apoptotic pathways, one with immune mediators and metabolism, one with transcriptional factors, and one with chromatin remodeling. Thus, informatics of a neuronal gene array allowed us to determine the expression profile of mitochondrial genes in the amygdala complex of an animal model of PTSD. The result is a further understanding of the metabolic and neuronal signaling mechanisms associated with delayed and exaggerated fear.

Mitochondrial targeted neuron focused genes in hippocampus of rats with traumatic brain injury.

CONTEXT: Mild traumatic brain injury (mTBI) represents a major health problem in civilian populations as well as among the military service members due to (1) lack of effective treatments, and (2) our incomplete understanding about the progression of secondary cell injury cascades resulting in neuronal cell death due to deficient cellular energy metabolism and damaged mitochondria. AIMS: The aim of this study was to identify and delineate the mitochondrial targeted genes responsible for altered brain energy metabolism in the injured brain. SETTINGS AND DESIGN: Rats were either grouped into naïve controls or received lateral fluid percussion brain injury (2-2.5 atm) and followed up for 7 days. MATERIALS AND METHODS: Rats were either grouped into naïve controls or received lateral fluid percussion brain injury (2-2.5 atm) and followed for 7 days. The severity of brain injury was evaluated by the neurological severity scale-revised (NSS-R) at 3 and 5 days post TBI and immunohistochemical analyses at 7 days post TBI. The expression profiles of mitochondrial-targeted genes across the hippocampus from TBI and naïe rats were also examined by oligo-DNA microarrays. RESULTS: NSS-R scores of TBI rats (5.4 ± 0.5) in comparison to naïe rats (3.9 ± 0.5) and H and E staining of brain sections suggested a mild brain injury. Bioinformatics and systems biology analyses showed 31 dysregulated genes, 10 affected canonical molecular pathways including a number of genes involved in mitochondrial enzymes for oxidative phosphorylation, mitogen-activated protein Kinase (MAP), peroxisome proliferator-activated protein (PPAP), apoptosis signaling, and genes responsible for long-term potentiation of Alzheimer's and Parkinson's diseases. CONCLUSIONS: Our results suggest that dysregulated mitochondrial-focused genes in injured brains may have a clinical utility for the development of future therapeutic strategies aimed at the treatment of TBI.

Biomarkers in an animal model for revealing neural, hematologic, and behavioral correlates of PTSD.

Identification of biomarkers representing the evolution of the pathophysiology of Post Traumatic Stress Disorder (PTSD) is vitally important, not only for objective diagnosis but also for the evaluation of therapeutic efficacy and resilience to trauma. Ongoing research is directed at identifying molecular biomarkers for PTSD, including traumatic stress induced proteins, transcriptomes, genomic variances and genetic modulators, using biologic samples from subjects' blood, saliva, urine, and postmortem brain tissues. However, the correlation of these biomarker molecules in peripheral or postmortem samples to altered brain functions associated with psychiatric symptoms in PTSD remains unresolved. Here, we present an animal model of PTSD in which both peripheral blood and central brain biomarkers, as well as behavioral phenotype, can be collected and measured, thus providing the needed correlation of the central biomarkers of PTSD, which are mechanistic and pathognomonic but cannot be collected from people, with the peripheral biomarkers and behavioral phenotypes, which can. Our animal model of PTSD employs restraint and tail shocks repeated for three continuous days - the inescapable tail-shock model (ITS) in rats. This ITS model mimics the pathophysiology of PTSD (17, 7, 4, 10). We and others have verified that the ITS model induces behavioral and neurobiological alterations similar to those found in PTSD subjects (17, 7, 10, 9). Specifically, these stressed rats exhibit (1) a delayed and exaggerated startle response appearing several days after stressor cessation, which given the compressed time scale of the rat's life compared to a humans, corresponds to the one to three months delay of symptoms in PTSD patients (DSM-IV-TR PTSD Criterian D/E (13)), (2) enhanced plasma corticosterone (CORT) for several days, indicating compromise of the hypothalamopituitary axis (HPA), and (3) retarded body weight gain after stressor cessation, indicating dysfunction of metabolic regulation. The experimental paradigms employed for this model are: (1) a learned helplessness paradigm in the rat assayed by measurement of acoustic startle response (ASR) and a charting of body mass; (2) microdissection of the rat brain into regions and nuclei; (3) enzyme-linked immunosorbent assay (ELISA) for blood levels of CORT; (4) a gene expression microarray plus related bioinformatics tools (18). This microarray, dubbed rMNChip, focuses on mitochondrial and mitochondria-related nuclear genes in the rat so as to specifically address the neuronal bioenergetics hypothesized to be involved in PTSD.

Rat mitochondrion-neuron focused microarray (rMNChip) and bioinformatics tools for rapid identification of differential pathways in brain tissues.

Mitochondrial function is of particular importance in brain because of its high demand for energy (ATP) and efficient removal of reactive oxygen species (ROS). We developed rat mitochondrion-neuron focused microarray (rMNChip) and integrated bioinformatics tools for rapid identification of differential pathways in brain tissues. rMNChip contains 1,500 genes involved in mitochondrial functions, stress response, circadian rhythms and signal transduction. The bioinformatics tool includes an algorithm for computing of differentially expressed genes, and a database for straightforward and intuitive interpretation for microarray results. Our application of these tools to RNA samples derived from rat frontal cortex (FC), hippocampus (HC) and hypothalamus (HT) led to the identification of differentially-expressed signal-transduction-bioenergenesis and neurotransmitter-synthesis pathways with a dominant number of genes (FC/HC = 55/6; FC/HT = 55/4) having significantly (p<0.05, FDR<10.70%) higher (≥1.25 fold) RNA levels in the frontal cortex than the others, strongly suggesting active generation of ATP and neurotransmitters and efficient removal of ROS. Thus, these tools for rapid and efficient identification of differential pathways in brain regions will greatly facilitate our systems-biological study and understanding of molecular mechanisms underlying complex and multifactorial neurodegenerative diseases.

Automatic Spot Identification for High Throughput Microarray Analysis.

High throughput microarray analysis has great potential in scientific research, disease diagnosis, and drug discovery. A major hurdle toward high throughput microarray analysis is the time and effort needed to accurately locate gene spots in microarray images. An automatic microarray image processor will allow accurate and efficient determination of spot locations and sizes so that gene expression information can be reliably extracted in a high throughput manner. Current microarray image processing tools require intensive manual operations in addition to the input of grid parameters to correctly and accurately identify gene spots. This work developed a method, herein called auto-spot, to automate the spot identification process. Through a series of correlation and convolution operations, as well as pixel manipulations, this method makes spot identification an automatic and accurate process. Testing with real microarray images has demonstrated that this method is capable of automatically extracting subgrids from microarray images and determining spot locations and sizes within each subgrid, regardless of variations in array patterns and background noises. With this method, we are one step closer to the goal of high throughput microarray analysis.

Undetectable and Decreased Expression of KIAA1949 (Phostensin) Encoded on Chromosome 6p21.33 in Human Breast Cancers Revealed by Transcriptome Analysis.

Cytogenetic aberration and loss of heterozygosity (LOH) are documented on chromosome 6 in many cancers and the introduction of a neo-tagged chromosome 6 into breast cancer cell lines mediates suppression of tumorigenicity. In this study, we described the identification of KIAA1949 (phostensin) as a putative tumor suppressor gene. Our microarray analysis screened 25,985 cDNAs between a tumorigenic and metastatic breast cancer cell line MDA-MB-231 and the chromosome 6-mediated suppressed, non-tumorigenic and non-metastatic derivative cell line MDA/H6, resulting in the identification of 651 differentially expressed genes. Using customized microarrays containing these 651 cDNAs and 117 controls, we identified 200 frequently dysregulated genes in 10 breast cancer cell lines and 5 tumor tissues using MDA/H6 as reference. Our bioinformatics analysis revealed that chromosome 6 encodes 25 of these 200 genes, with 4 downregulation and 21 upergulation. Northern analysis validated microarray results and was used to detect the size and number of RNA species of 2 downregulated (KIAA1949, PTK7) and 2 upregulated (SFRS3, HMGN3) genes in the cell lines. Particularly, the KIAA1949 gene at 6p21.33, a band region with the frequent cytogenetic aberration and LOH encodes 4 poly(A)-RNA species (4.6-, 4-, 3- and 1.5-kb) in multiple normal and breast cancer samples. Microarray analysis revealed KIAA1949 downregulation in 86.7% (n=15) of breast cancer cell lines and tumor tissues. Northern analysis demonstrated undetectable and decreased expression of KIAA1949 in 100% (n=10) of breast cancer cell lines. Taken together, these results strongly suggest KIAA1949 is a novel candidate breast cancer suppressor gene.

Increased invasiveness and aggressiveness in breast epithelia with cytoplasmic p63 expression.

Our previous studies revealed that pregnancy associated breast cancer (PABC) had significantly reduced nuclear p63 expression in myoepithelia, while intense cytoplasmic p63 expression in associated epithelia. Our current study assessed these epithelia using immunohistochemistry with a panel of aggressiveness and invasiveness related markers and comparative genomic hybridization (array-CGH) with over 30,000 DNA probes. These epithelia showed several unique alterations, including (1) immunohistochemical and morphological resemblance to invasive cancer, (2) significant gain in copy numbers of DNA coding genes for morphogenesis, angiogenesis, and metastasis, and (3) significant loss in copy numbers of DNA coding genes for tumor suppressors, cell adhesion, and macromolecular complex assembly or intra-cellular trafficking. Detected array-CGH alterations correlated well with in vivo expression of a number of corresponding proteins tested. These findings suggest that aberrant sub-cellular localization of p63 expression in normal or hyperplastic appearing epithelial cells may significant contribute to increased invasiveness and aggressiveness of these cells.

Mitochondrial gene profiling: translational perspectives.

The last decade has witnessed the development of multiple microarray platforms designed to study, in a comprehensive fashion, the expression and sequence of both mitochondrial and nuclear genes that encode mitochondrial proteins. Mitochondrial dysfunction has been implicated in a number of severe medical conditions including cancer, metabolic diseases (i.e., cardiovascular, diabetes and obesity) and neurodegenerative disorders and it is responsible for the adverse effects of numerous drugs. Profiling of the genetic and genomic status of mitochondria with focused microarrays offers the promise of rapidly and robustly identifying novel biomarkers for early disease diagnoses and prognoses, predicting of drug safety, liability, and selecting and stratifying of patients in clinical trials.

Brx mediates the response of lymphocytes to osmotic stress through the activation of NFAT5.

Extracellular hyperosmolarity, or osmotic stress, generally caused by differences in salt and macromolecule concentrations across the plasma membrane, occurs in lymphoid organs and at inflammatory sites. The response of immune cells to osmotic stress is regulated by nuclear factor of activated T cells 5 (NFAT5), a transcription factor that induces the expression of hyperosmolarity-responsive genes and stimulates cytokine production. We report that the guanine nucleotide exchange factor (GEF) Brx [also known as protein kinase A-anchoring protein 13 (AKAP13)] is essential for the expression of nfat5 in response to osmotic stress, thus transmitting the extracellular hyperosmolarity signal and enabling differentiation of splenic B cells and production of immunoglobulin. This process required the activity of p38 mitogen-activated protein kinase (MAPK) and NFAT5 and involved a physical interaction between Brx and c-Jun N-terminal kinase (JNK)-interacting protein 4 (JIP4), a scaffold molecule specific to activation of the p38 MAPK cascade. Our results indicate that Brx integrates the responses of immune cells to osmotic stress and inflammation by elevating intracellular osmolarity and stimulating the production of cytokines.

Two types of human malignant melanoma cell lines revealed by expression patterns of mitochondrial and survival-apoptosis genes: implications for malignant melanoma therapy.

Human malignant melanoma has poor prognosis because of resistance to apoptosis and therapy. We describe identification of the expression profile of 1,037 mitochondria-focused genes and 84 survival-apoptosis genes in 21 malignant melanoma cell lines and 3 normal melanocyte controls using recently developed hMitChip3 cDNA microarrays. Unsupervised hierarchical clustering analysis of 1,037 informative genes, and 84 survival-apoptosis genes, classified these malignant melanoma cell lines into type A (n = 12) and type B (n = 9). Three hundred fifty-five of 1,037 (34.2%) genes displayed significant (P ≤ 0.030; false discovery rate ≤ 3.68%) differences (± ≥ 2.0-fold) in average expression, with 197 genes higher and 158 genes lower in type A than in type B. Of 84 genes with known survival-apoptosis functions, 38 (45.2%) displayed the significant (P < 0.001; false discovery rate < 0.15%) difference. Antiapoptotic (BCL2, BCL2A1, PPARD, and RAF1), antioxidant (MT3, PRDX5, PRDX3, GPX4, GLRX2, and GSR), and proapoptotic (BAD, BNIP1, APAF1, BNIP3L, CASP7, CYCS, CASP1, and VDAC1) genes expressed at higher levels in type A than in type B, whereas the different set of antiapoptotic (PSEN1, PPP2CA, API5, PPP2R1B, PPP2R1A, and FIS1), antioxidant (HSPD1, GSS, SOD1, ATOX1, and CAT), and proapoptotic (ENDOG, BAK1, CASP2, CASP4, PDCD5, HTRA2, SEPT4, TNFSF10, and PRODH) genes expressed at lower levels in type A than in type B. Microarray data were validated by quantitative reverse transcription-PCR. These results showed the presence of two types of malignant melanoma, each with a specific set of dysregulated survival-apoptosis genes, which may prove useful for development of new molecular targets for therapeutic intervention and novel diagnostic biomarkers for treatment and prognosis of malignant melanoma.

Dysregulated mitochondrial genes and networks with drug targets in postmortem brain of patients with posttraumatic stress disorder (PTSD) revealed by human mitochondria-focused cDNA microarrays.

Posttraumatic stress disorder (PTSD) is associated with decreased activity in the dorsolateral prefrontal cortex (DLPFC), the brain region that regulates working memory and preparation and selection of fear responses. We investigated gene expression profiles in DLPFC Brodmann area (BA) 46 of postmortem patients with (n=6) and without PTSD (n=6) using human mitochondria-focused cDNA microarrays. Our study revealed PTSD-specific expression fingerprints of 800 informative mitochondria-focused genes across all of these 12 BA46 samples, and 119 (+/->1.25, p<0.05) and 42 (+/->1.60, p<0.05) dysregulated genes between the PTSD and control samples. Quantitative RT-PCR validated the microarray results. These fingerprints can essentially distinguish the PTSD DLPFC BA46 brains from controls. Of the 119 dysregulated genes (+/-> or =125%, p<0.05), the highest percentages were associated with mitochondrial dysfunction (4.8%, p=6.61 x 10(-6)), oxidative phosphorylation (3.8%, p=9.04 x 10(-4)), cell survival-apoptosis (25.2%, p<0.05) and neurological diseases (23.5%, p<0.05). Fifty (50) dysregulated genes were present in the molecular networks that are known to be involved in neuronal function-survival and contain 7 targets for neuropsychiatric drugs. Thirty (30) of the dysregulated genes are associated with a number of neuropsychiatric disorders. Our results indicate mitochondrial dysfunction in the PTSD DLPFC BA46 and provide the expression fingerprints that may ultimately serve as biomarkers for PTSD diagnosis and the drugs and molecular targets that may prove useful for development of remedies for prevention and treatment of PTSD.