RNA-sequencing of the brain transcriptome implicates dysregulation of neuroplasticity, circadian rhythms and GTPase binding in bipolar disorder.

RNA-sequencing (RNA-seq) is a powerful technique to investigate the complexity of gene expression in the human brain. We used RNA-seq to survey the brain transcriptome in high-quality postmortem dorsolateral prefrontal cortex from 11 individuals diagnosed with bipolar disorder (BD) and from 11 age- and gender-matched controls. Deep sequencing was performed, with over 350 million reads per specimen. At a false discovery rate of <5%, we detected five differentially expressed (DE) genes and 12 DE transcripts, most of which have not been previously implicated in BD. Among these, Prominin 1/CD133 and ATP-binding cassette-sub-family G-member2 (ABCG2) have important roles in neuroplasticity. We also show for the first time differential expression of long noncoding RNAs (lncRNAs) in BD. DE transcripts include those of serine/arginine-rich splicing factor 5 (SRSF5) and regulatory factor X4 (RFX4), which along with lncRNAs have a role in mammalian circadian rhythms. The DE genes were significantly enriched for several Gene Ontology categories. Of these, genes involved with GTPase binding were also enriched for BD-associated SNPs from previous genome-wide association studies, suggesting that differential expression of these genes is not simply a consequence of BD or its treatment. Many of these findings were replicated by microarray in an independent sample of 60 cases and controls. These results highlight common pathways for inherited and non-inherited influences on disease risk that may constitute good targets for novel therapies.

Upregulation of α-synuclein during localized radiation therapy signals the association of cancer-related fatigue with the activation of inflammatory and neuroprotective pathways.

PURPOSE: Neuroinflammatory mechanisms are associated with fatigue in neurodegenerative conditions such as Parkinson's. The symptoms in Parkinson's including fatigue are thought to be related to α-synuclein overexpression. This study investigated genomic correlates of fatigue experienced by men with prostate cancer receiving external beam radiation therapy (EBRT). PATIENTS AND METHODS: Sixteen men with non-metastatic prostate cancer who were scheduled to receive EBRT were enrolled. Fatigue scores and blood were obtained at baseline (prior to EBRT, D0); one hour following initiation of EBRT (D1), day 7 (D7), day 14 (D14), midpoint (days 19-21, D21), completion (days 38-42, D42), and four weeks post-EBRT (days 68-72, D72). Gene expression profiling using microarray analysis was performed from peripheral blood and confirmatory qPCR and protein (ELISA) analyses verified the microarray results. Correlations between fatigue and gene/protein expressions were determined using a mixed model approach. RESULTS: Microarray data showed significant, differential expression of 463 probesets following EBRT. SNCA had a 2.95-fold change at D21 from baseline. SNCA expression was confirmed by qPCR (p<0.001) and ELISA (p<0.001) over time during EBRT. Fatigue scores were significantly correlated with SNCA gene expression on D14 (r=0.55, p<0.05) and plasma α-synuclein concentrations on D42 of EBRT (r=0.54, p=0.04). CONCLUSION: Fatigue experienced during EBRT may be mediated by α-synuclein overexpression. Alpha-synuclein may serve as a useful biomarker to understand the mechanisms and pathways related to the development of fatigue in this population.

Cooperative and redundant signaling of leukotriene B4 and leukotriene D4 in human monocytes.

BACKGROUND: Leukotriene B(4) (LTB(4)) and cysteinyl leukotrienes (cysLTs) are important immune mediators, often found concomitantly at sites of inflammation. Although some of the leukotriene-mediated actions are distinctive (e.g., bronchial constriction for cysLTs), many activities such as leukocyte recruitment to tissues and amplification of inflammatory responses are shared by both classes of leukotrienes. OBJECTIVE: We used human monocytes to characterize leukotriene-specific signaling, gene expression signatures, and functions and to identify interactions between LTB(4)- and cysLTs-induced pathways. METHODS: Responsiveness to leukotrienes was assessed using oligonucleotide microarrays, real-time PCR, calcium mobilization, kinase activation, and chemotaxis assays. RESULTS: Human monocytes were found to express mRNA for high- and low-affinity LTB(4) receptors, BLT(1) and BLT(2), but signal predominantly through BLT(1) in response to LTB(4) stimulation as shown using selective agonists, inhibitors, and gene knock down experiments. LTB(4) acting through BLT(1) coupled to G-protein α inhibitory subunit activated calcium signaling, p44/42 mitogen-activated protein kinase, gene expression, and chemotaxis. Twenty-seven genes, including immediate early genes (IEG), transcription factors, cytokines, and membrane receptors were significantly up-regulated by LTB(4). LTB(4) and LTD(4) had similar effects on signaling, gene expression, and chemotaxis indicating redundant cell activation pathways but costimulation with both lipid mediators was additive for many monocyte functions. CONCLUSION: Leukotriene B(4) and LTD(4) display both redundant and cooperative effects on intracellular signaling, gene expression, and chemotaxis in human monocytes. These findings suggest that therapies targeting either leukotriene alone may be less effective than approaches directed at both.

Transcriptome analysis and molecular signature of human retinal pigment epithelium.

Retinal pigment epithelium (RPE) is a polarized cell layer critical for photoreceptor function and survival. The unique physiology and relationship to the photoreceptors make the RPE a critical determinant of human vision. Therefore, we performed a global expression profiling of native and cultured human fetal and adult RPE and determined a set of highly expressed 'signature' genes by comparing the observed RPE gene profiles to the Novartis expression database (SymAtlas: http://wombat.gnf.org/index.html) of 78 tissues. Using stringent selection criteria of at least 10-fold higher expression in three distinct preparations, we identified 154 RPE signature genes, which were validated by qRT-PCR analysis in RPE and in an independent set of 11 tissues. Several of the highly expressed signature genes encode proteins involved in visual cycle, melanogenesis and cell adhesion and Gene ontology analysis enabled the assignment of RPE signature genes to epithelial channels and transporters (ClCN4, BEST1, SLCA20) or matrix remodeling (TIMP3, COL8A2). Fifteen RPE signature genes were associated with known ophthalmic diseases, and 25 others were mapped to regions of disease loci. An evaluation of the RPE signature genes in a recently completed AMD genomewide association (GWA) data set revealed that TIMP3, GRAMD3, PITPNA and CHRNA3 signature genes may have potential roles in AMD pathogenesis and deserve further examination. We propose that RPE signature genes are excellent candidates for retinal diseases and for physiological investigations (e.g. dopachrome tautomerase in melanogenesis). The RPE signature gene set should allow the validation of RPE-like cells derived from human embryonic or induced pluripotent stem cells for cell-based therapies of degenerative retinal diseases.

Distinct gene expression profiles in norepinephrine- and epinephrine-producing hereditary and sporadic pheochromocytomas: activation of hypoxia-driven angiogenic pathways in von Hippel-Lindau syndrome.

Pheochromocytomas in von Hippel-Lindau (VHL) syndrome produce exclusively norepinephrine, whereas those in multiple endocrine neoplasia type 2 (MEN 2) produce epinephrine. This study examined the pathways activated in VHL-associated pheochromocytomas by comparing gene expression profiles in VHL and MEN 2 tumors in relationship to profiles in sporadic norepinephrine- and epinephrine-producing tumors. Larger and more distinct differences in gene expression among hereditary than sporadic tumors indicated the importance of the underlying mutation to gene expression profiles. Many of the genes over-expressed in VHL compared with MEN 2 tumors were clearly linked to the hypoxia-driven angiogenic pathways that are activated in VHL-associated tumorigenesis. Such genes included those for the glucose transporter, vascular endothelial growth factor (VEGF), placental growth factor, angiopoietin 2, tie-1, VEGF receptor 2 and its coreceptor, neuropilin-1. Other up-regulated genes, such as connective tissue growth factor, cysteine-rich 61, matrix metalloproteinase 1, vascular endothelial cadherin, tenascin C, stanniocalcin 1, and cyclooxygenases 1 and 2 are known to be involved in VEGF-regulated angiogenesis. Shared differences in expression of subsets of genes in norepinephrine- versus epinephrine-producing hereditary and sporadic pheochromocytomas indicated other differences in gene expression that may underlie the biochemical phenotype. Over-expression of the hypoxia-inducible transcription factor, HIF-2alpha, in norepinephrine-predominant sporadic and VHL tumors compared with epinephrine-producing tumors indicates that expression of this gene depends on the noradrenergic biochemical phenotype. The findings fit with the known expression of HIF-2alpha in norepinephrine-producing cells of the sympathetic nervous system and might explain both the development and noradrenergic biochemical phenotype of pheochromocytomas in VHL syndrome.

Proteomic profiling of the cancer microenvironment by antibody arrays.

Critical changes in protein expression that enable tumors to initiate and progress originate in the local tissue microenvironment, and there are increasing indications that these microenvironmental alterations in protein expression play critical roles in shaping and directing this process. As a model to better understand how patterns of protein expression shape the tissue microenvironment, we analyzed protein expression in tissue derived from squamous cell carcinoma of the oral cavity through an antibody microarray approach for high-throughput proteomic analysis. Utilizing laser capture microdissection to procure total protein from specific microscopic cellular populations, we demonstrate that quantitative, and potentially qualitative, differences in expression patterns of multiple proteins within epithelial cells reproducibly correlate with oral cavity tumor progression. Furthermore, differential expression of multiple proteins was also found in stromal cells surrounding and adjacent to regions of diseased epithelium that directly correlated with tumor progression of the epithelium. Most of the proteins identified in both cell types are involved in signal transduction pathways, thus we hypothesize that extensive molecular communication involving complex cellular signaling between epithelium and stroma play a key role in driving oral cavity cancer progression.

Augmentation in expression of activation-induced genes differentiates memory from naive CD4+ T cells and is a molecular mechanism for enhanced cellular response of memory CD4+ T cells.

In an attempt to understand the molecular basis for the immunological memory response, we have used cDNA microarrays to measure gene expression of human memory and naive CD4+ T cells at rest and after activation. Our analysis of 54,768 cDNA clones provides the first glimpse into gene expression patterns of memory and naive CD4+ T cells at the genome-scale and reveals several novel findings. First, memory and naive CD4+ T cells expressed similar numbers of genes at rest and after activation. Second, we have identified 14 cDNA clones that expressed higher levels of transcripts in memory cells than in naive cells. Third, we have identified 135 (130 known genes and 5 expressed sequence tags) up-regulated and 68 (42 known genes and 26 expressed sequence tags) down-regulated cDNA clones in memory CD4+ T after in vitro stimulation with anti-CD3 plus anti-CD28. Interestingly, the increase in mRNA levels of up-regulated genes was greater in memory than in naive CD4+ T cells after in vitro stimulation and was higher with anti-CD3 plus anti-CD28 than with anti-CD3 alone in both memory and naive CD4+ T cells. Finally, the changes in expression of actin and cytokine genes identified by cDNA microarrays were confirmed by Northern and protein analyses. Together, we have identified approximately 200 cDNA clones whose expression levels changed after activation and suggest that the level of expression of up-regulated genes is a molecular mechanism that differentiates the response of memory from naive CD4+ T cells.

Sequence databases and microarrays as tools for identifying prostate cancer biomarkers.

Identification, acquisition, and assessment of molecular markers that could be adopted as surrogate endpoints for evaluating a response to prostate cancer intervention strategies is highly desirable. Recent advances in the fields of genomics and biotechnology have dramatically increased the quantity and accessibility of molecular information that is relevant to the study of prostate carcinogenesis. One major advance involves the construction of comprehensive databases that archive gene sequences and gene expression data. This information is in a format suitable for virtual queries designed to distinguish the molecular differences between normal and cancer cells. A second major advance uses robotic tools to construct microarrays comprising thousands of distinct genes expressed in prostate tissues. Such arrays offer a powerful approach for monitoring the expression of thousands of genes simultaneously and provide access for techniques designed to assess patterns or "fingerprints" of gene expression that may ultimately be used as signatures of response to therapeutic intervention.

Doctor lounge talk.

The discussions that take place in the doctor's lounge are essential to the well-being of the preceptor. Doctor lounge talk should be filtered for students with their level of development in mind.

Development of a prostate cDNA microarray and statistical gene expression analysis package.

A cDNA microarray comprising 5184 different cDNAs spotted onto nylon membrane filters was developed for prostate gene expression studies. The clones used for arraying were identified by cluster analysis of > 35 000 prostate cDNA library-derived expressed sequence tags (ESTs) present in the dbEST database maintained by the National Center for Biotechnology Information. Total RNA from two cell lines, prostate line 8.4 and melanoma line UACC903, was used to make radiolabeled probe for filter hybridizations. The absolute intensity of each individual cDNA spot was determined by phosphorimager scanning and evaluated by a bioinformatics package developed specifically for analysis of cDNA microarray experimentation. Results indicated 89% of the genes showed intensity levels above background in prostate cells compared with only 28% in melanoma cells. Replicate probe preparations yielded results with correlation values ranging from r = 0.90 to 0.93 and coefficient of variation ranging from 16 to 28%. Findings indicate that among others, the keratin 5 and vimentin genes were differentially expressed between these two divergent cell lines. Follow-up northern blot analysis verified these two expression changes, thereby demonstrating the reliability of this system. We report the development of a cDNA microarray system that is sensitive and reliable, demonstrates a low degree of variability, and is capable of determining verifiable gene expression differences between two distinct human cell lines. This system will prove useful for differential gene expression analysis in prostate-derived cells and tissue.

The diagnosis of low-grade peripheral B-cell neoplasms in bone marrow trephines.

The aim of this study was to establish how effective is the use of immunohistochemistry on formalin-fixed bone marrow in diagnosing low-grade B-cell neoplasms. We investigated a series of 41 consecutive patients with bone marrow involvement for whom no other diagnostic tissues were available. The sections were stained with the following antibodies: CD3, CD20, CD79a, CD5, CD10, CD23, anti-cyclin D1 and kappa and lambda light chains. Antigen retrieval was performed using either a microwave oven or a pressure cooker. Labelling was performed with an avidin-biotin-peroxidase labelling system. A final diagnosis was reached in 37 out of 41 cases (90%): B-chronic lymphocytic leukaemia (15 cases), follicular lymphoma (10 cases), mantle-cell lymphoma (eight cases) and lymphoplasmacytoid lymphoma/immunocytoma (four cases). In the remaining four cases, a generic diagnosis of low-grade B-cell neoplasm was made. The immunophenotyping of formalin-fixed marrow is a useful technique for diagnosing most of the low-grade B-cell neoplasms.

Pronounced and sustained central hypernoradrenergic function in major depression with melancholic features: relation to hypercortisolism and corticotropin-releasing hormone.

Both stress-system activation and melancholic depression are characterized by fear, constricted affect, stereotyped thinking, and similar changes in autonomic and neuroendocrine function. Because norepinephrine (NE) and corticotropin-releasing hormone (CRH) can produce these physiological and behavioral changes, we measured the cerebrospinal fluid (CSF) levels each hour for 30 consecutive hours in controls and in patients with melancholic depression. Plasma adrenocorticotropic hormone (ACTH) and cortisol levels were obtained every 30 min. Depressed patients had significantly higher CSF NE and plasma cortisol levels that were increased around the clock. Diurnal variations in CSF NE and plasma cortisol levels were virtually superimposable and positively correlated with each other in both patients and controls. Despite their hypercortisolism, depressed patients had normal levels of plasma ACTH and CSF CRH. However, plasma ACTH and CSF CRH levels in depressed patients were inappropriately high, considering the degree of their hypercortisolism. In contrast to the significant negative correlation between plasma cortisol and CSF CRH levels seen in controls, patients with depression showed no statistical relationship between these parameters. These data indicate that persistent stress-system dysfunction in melancholic depression is independent of the conscious stress of the disorder. These data also suggest mutually reinforcing bidirectional links between a central hypernoradrenergic state and the hyperfunctioning of specific central CRH pathways that each are driven and sustained by hypercortisolism. We postulate that alpha-noradrenergic blockade, CRH antagonists, and treatment with antiglucocorticoids may act at different loci, alone or in combination, in the treatment of major depression with melancholic features.

Analysis of gene expression in mutiple sclerosis lesions using cDNA microarrays.

In multiple sclerosis (MS) patients, a coordinated attack of the immune system against the primary constituents of oligodendrocytes and/or the myelin sheath of oligodendrocytes results in the formation of lesions in the brain and spinal cord. Thus far, however, a limited number of genes that potentially contribute to lesion pathology have been identified. Using cDNA microarray technology, we have performed experiments on MS tissue monitoring the expression pattern of over 5,000 genes and compared the gene expression profile of normal white matter with that found in acute lesions from the brain of a single MS patient. Sixty-two differentially expressed genes were identified, including the Duffy chemokine receptor, interferon regulatory factor-2, and tumor necrosis factor alpha receptor-2 among others. Thus, cDNA microarray technology represents a powerful new tool for the identification of genes not previously associated with the MS disease process.

Comparing protein sequence-based and predicted secondary structure-based methods for identification of remote homologs.

We have compared a novel sequence-structure matching technique, FORESST, for detecting remote homologs to three existing sequence based methods, including local amino acid sequence similarity by BLASTP, hidden Markov models (HMMs) of sequences of protein families using SAM, HMMs based on sequence motifs identified using meta-MEME. FORESST compares predicted secondary structures to a library of structural families of proteins, using HMMs. Altogether 45 proteins from nine structural families in the database CATH were used in a cross-validated test of the fold assignment accuracy of each method. Local sequence similarity of a query sequence to a protein family is measured by the highest segment pair (HSP) score. Each of the HMM-based approaches (FORESST, MEME, amino acid sequence-based HMM) yielded log-odds score for the query sequence. In order to make a fair comparison among these methods, the scores for each method were converted to Z-scores in a uniform way by comparing the raw scores of a query protein with the corresponding scores for a set of unrelated proteins. Z-Scores were analyzed as a function of the maximum pairwise sequence identity (MPSID) of the query sequence to sequences used in training the model. For MPSID above 20%, the Z-scores increase linearly with MPSID for the sequence-based methods but remain roughly constant for FORESST. Below 15%, average Z-scores are close to zero for the sequence-based methods, whereas the FORESST method yielded average Z-scores of 1.8 and 1.1, using observed and predicted secondary structures, respectively. This demonstrates the advantage of the sequence-structure method for detecting remote homologs.

Recognizing the pleckstrin homology domain fold in mammalian phospholipase D using hidden Markov models.

Phospholipase D was first described in plant tissue but has recently been shown to occur in mammalian cells where it is activated by cell surface receptors. Its mode of activation by receptors in unclear. Biochemical studies suggest that it may occur downstream of other effector proteins and that small GTP-dependent regulatory proteins may be involved. The sequence in a non-designated region of mammalian phospholipase D1 and 2 shows similarity to a structural domain that is present in signalling proteins that are regulated by protein kinases or heterotrimeric G-proteins. Mammalian phospholipase D has structural similarities with other lipid signalling phospholipases and thus may be regulated by receptors in an analogous fashion.

FORESST: fold recognition from secondary structure predictions of proteins.

MOTIVATION: A method for recognizing the three-dimensional fold from the protein amino acid sequence based on a combination of hidden Markov models (HMMs) and secondary structure prediction was recently developed for proteins in the Mainly-Alpha structural class. Here, this methodology is extended to Mainly-Beta and Alpha-Beta class proteins. Compared to other fold recognition methods based on HMMs, this approach is novel in that only secondary structure information is used. Each HMM is trained from known secondary structure sequences of proteins having a similar fold. Secondary structure prediction is performed for the amino acid sequence of a query protein. The predicted fold of a query protein is the fold described by the model fitting the predicted sequence the best. RESULTS: After model cross-validation, the success rate on 44 test proteins covering the three structural classes was found to be 59%. On seven fold predictions performed prior to the publication of experimental structure, the success rate was 71%. In conclusion, this approach manages to capture important information about the fold of a protein embedded in the length and arrangement of the predicted helices, strands and coils along the polypeptide chain. When a more extensive library of HMMs representing the universe of known structural families is available (work in progress), the program will allow rapid screening of genomic databases and sequence annotation when fold similarity is not detectable from the amino acid sequence. AVAILABILITY: FORESST web server at http://absalpha.dcrt.nih.gov:8008/ for the library of HMMs of structural families used in this paper. FORESST web server at http://www.tigr.org/ for a more extensive library of HMMs (work in progress). CONTACT: valedf@tigr.org; munson@helix.nih.gov; garnier@helix.nih.gov

The impact of early clinical training in medical education: a multi-institutional assessment.

With funding from The Robert Wood Johnson Foundation's Generalist Physician Initiative, Dartmouth Medical School (DMS), New York Medical College (NYMC), and Virginia Commonwealth University School of Medicine (VCU-SOM) adopted early community-based training models for longitudinal clinical experiences. These schools developed different evaluation strategies to assess these models. This paper describes each program, the method used to evaluate an aspect of the program, lessons learned about early clinical teaching and learning, and challenges encountered. Each program used cross-sectional evaluation, and the analysis methods included descriptive statistics, chi-square, t-tests, analysis of variance, and generalized linear models. Dartmouth determined that the type of preceptor does not greatly influence the development of clinical skills, although case-specific differences were discovered. NYMC learned that students taught clinical skills in community-based settings performed as well as or better than their peers who received early patient experience on hospital wards. Virginia Commonwealth discovered that community experiences contributed positively to students' education, critical thinking, and problem-solving skills. Students value early clinical experiences and make important achievements in clinical skills and knowledge development, although logistic challenges exist in conducting these courses. Evaluations are critical to ensure competency, and faculty development must be linked to the evaluation process.

Statistical significance of hierarchical multi-body potentials based on Delaunay tessellation and their application in sequence-structure alignment.

Statistical potentials based on pairwise interactions between C alpha atoms are commonly used in protein threading/fold-recognition attempts. Inclusion of higher order interaction is a possible means of improving the specificity of these potentials. Delaunay tessellation of the C alpha-atom representation of protein structure has been suggested as a means of defining multi-body interactions. A large number of parameters are required to define all four-body interactions of 20 amino acid types (20(4) = 160,000). Assuming that residue order within a four-body contact is irrelevant reduces this to a manageable 8,855 parameters, using a nonredundant dataset of 608 protein structures. Three lines of evidence support the significance and utility of the four-body potential for sequence-structure matching. First, compared to the four-body model, all lower-order interaction models (three-body, two-body, one-body) are found statistically inadequate to explain the frequency distribution of residue contacts. Second, coherent patterns of interaction are seen in a graphic presentation of the four-body potential. Many patterns have plausible biophysical explanations and are consistent across sets of residues sharing certain properties (e.g., size, hydrophobicity, or charge). Third, the utility of the multi-body potential is tested on a test set of 12 same-length pairs of proteins of known structure for two protocols: Sequence-recognizes-structure, where a query sequence is threaded (without gap) through the native and a non-native structure; and structure-recognizes-sequence, where a query structure is threaded by its native and another non-native sequence. Using cross-validated training, protein sequences correctly recognized their native structure in all 24 cases. Conversely, structures recognized the native sequence in 23 of 24 cases. Further, the score differences between correct and decoy structures increased significantly using the three- or four-body potential compared to potentials of lower order.