A conserved BDNF, glutamate- and GABA-enriched gene module related to human depression identified by coexpression meta-analysis and DNA variant genome-wide association studies

Large scale gene expression (transcriptome) analysis and genome-wide association studies (GWAS) for single nucleotide polymorphisms have generated a considerable amount of gene- and disease-related information, but heterogeneity and various sources of noise have limited the discovery of disease mechanisms. As systematic dataset integration is becoming essential, we developed methods and performed meta-clustering of gene coexpression links in 11 transcriptome studies from postmortem brains of human subjects with major depressive disorder (MDD) and non-psychiatric control subjects. We next sought enrichment in the top 50 meta-analyzed coexpression modules for genes otherwise identified by GWAS for various sets of disorders. One coexpression module of 88 genes was consistently and significantly associated with GWAS for MDD, other neuropsychiatric disorders and brain functions, and for medical illnesses with elevated clinical risk of depression, but not for other diseases. In support of the superior discriminative power of this novel approach, we observed no significant enrichment for GWAS-related genes in coexpression modules extracted from single studies or in meta-modules using gene expression data from non-psychiatric control subjects. Genes in the identified module encode proteins implicated in neuronal signaling and structure, including glutamate metabotropic receptors (GRM1, GRM7), GABA receptors (GABRA2, GABRA4), and neurotrophic and development-related proteins [BDNF, reelin (RELN), Ephrin receptors (EPHA3, EPHA5)]. These results are consistent with the current understanding of molecular mechanisms of MDD and provide a set of putative interacting molecular partners, potentially reflecting components of a functional module across cells and biological pathways that are synchronously recruited in MDD, other brain disorders and MDD-related illnesses. Collectively, this study demonstrates the importance of integrating transcriptome data, gene coexpression modules and GWAS results for providing novel and complementary approaches to investigate the molecular pathology of MDD and other complex brain disorders. © 2014 Chang et al

A computational method for genotype calling in family-based sequencing data

Background: As sequencing technologies can help researchers detect common and rare variants across the human genome in many individuals, it is known that jointly calling genotypes across multiple individuals based on linkage disequilibrium (LD) can facilitate the analysis of low to modest coverage sequence data. However, genotype-calling methods for family-based sequence data, particularly for complex families beyond parent-offspring trios, are still lacking. Results: In this study, first, we proposed an algorithm that considers both linkage disequilibrium (LD) patterns and familial transmission in nuclear and multi-generational families while retaining the computational efficiency. Second, we extended our method to incorporate external reference panels to analyze family-based sequence data with a small sample size. In simulation studies, we show that modeling multiple offspring can dramatically increase genotype calling accuracy and reduce phasing and Mendelian errors, especially at low to modest coverage. In addition, we show that using external panels can greatly facilitate genotype calling of sequencing data with a small number of individuals. We applied our method to a whole genome sequencing study of 1339 individuals at ~10X coverage from the Minnesota Center for Twin and Family Research. Conclusions: The aggregated results show that our methods significantly outperform existing ones that ignore family constraints or LD information. We anticipate that our method will be useful for many ongoing family-based sequencing projects. We have implemented our methods efficiently in a C++ program FamLDCaller, which is available from http://www.pitt.edu/~wec47/famldcaller.html

Meta-analysis methods for combining multiple expression profiles: Comparisons, statistical characterization and an application guideline

Background: As high-throughput genomic technologies become accurate and affordable, an increasing number of data sets have been accumulated in the public domain and genomic information integration and meta-analysis have become routine in biomedical research. In this paper, we focus on microarray meta-analysis, where multiple microarray studies with relevant biological hypotheses are combined in order to improve candidate marker detection. Many methods have been developed and applied in the literature, but their performance and properties have only been minimally investigated. There is currently no clear conclusion or guideline as to the proper choice of a meta-analysis method given an application; the decision essentially requires both statistical and biological considerations.Results: We performed 12 microarray meta-analysis methods for combining multiple simulated expression profiles, and such methods can be categorized for different hypothesis setting purposes: (1) HSA: DE genes with non-zero effect sizes in all studies, (2) HSB: DE genes with non-zero effect sizes in one or more studies and (3) HSr: DE gene with non-zero effect in "majority"of studies. We then performed a comprehensive comparative analysis through six large-scale real applications using four quantitative statistical evaluation criteria: detection capability, biological association, stability and robustness. We elucidated hypothesis settings behind the methods and further apply multi-dimensional scaling (MDS) and an entropy measure to characterize the meta-analysis methods and data structure, respectively.Conclusions: The aggregated results from the simulation study categorized the 12 methods into three hypothesis settings (HSA, HSB, and HSr). Evaluation in real data and results from MDS and entropy analyses provided an insightful and practical guideline to the choice of the most suitable method in a given application. All source files for simulation and real data are available on the author's publication website. © 2013 Chang et al.; licensee BioMed Central Ltd

Sex chromosome complement regulates expression of mood-related genes

Background: Studies on major depressive and anxiety disorders suggest dysfunctions in brain corticolimbic circuits, including altered gamma-aminobutyric acid (GABA) and modulatory (serotonin and dopamine) neurotransmission. Interestingly, sexual dimorphisms in GABA, serotonin, and dopamine systems are also reported. Understanding the mechanisms behind these sexual dimorphisms may help unravel the biological bases of the heightened female vulnerability to mood disorders. Here, we investigate the contribution of sex-related factors (sex chromosome complement, developmental gonadal sex, or adult circulating hormones) to frontal cortex expression of selected GABA-, serotonin-, and dopamine-related genes. Methods: As gonadal sex is determined by sex chromosome complement, the role of sex chromosomes cannot be investigated individually in humans. Therefore, we used the Four Core Genotypes (FCG) mouse model, in which sex chromosome complement and gonadal sex are artificially decoupled, to examine the expression of 13 GABA-related genes, 6 serotonin- and dopamine-related genes, and 8 associated signal transduction genes under chronic stress conditions. Results were analyzed by three-way ANOVA (sex chromosome complement × gonadal sex × circulating testosterone). A global perspective of gene expression changes was provided by heatmap representation and gene co-expression networks to identify patterns of transcriptional activities related to each main factor. Results: We show that under chronic stress conditions, sex chromosome complement influenced GABA/serotonin/dopamine- related gene expression in the frontal cortex, with XY mice consistently having lower gene expression compared to XX mice. Gonadal sex and circulating testosterone exhibited less pronounced, more complex, and variable control over gene expression. Across factors, male conditions were associated with a tightly co-expressed set of signal transduction genes. Conclusions: Under chronic stress conditions, sex-related factors differentially influence expression of genes linked to mood regulation in the frontal cortex. The main factor influencing expression of GABA-, serotonin-, and dopamine-related genes was sex chromosome complement, with an unexpected pro-disease effect in XY mice relative to XX mice. This effect was partially opposed by gonadal sex and circulating testosterone, although all three factors influenced signal transduction pathways in males. Since GABA, serotonin, and dopamine changes are also observed in other psychiatric and neurodegenerative disorders, these findings have broader implications for the understanding of sexual dimorphism in adult psychopathology. © 2013 Seney et al.; licensee BioMed Central Ltd

Discovering monotonic stemness marker genes from time-series stem cell microarray data

© 2015 Wang et al.; licensee BioMed Central Ltd. Background: Identification of genes with ascending or descending monotonic expression patterns over time or stages of stem cells is an important issue in time-series microarray data analysis. We propose a method named Monotonic Feature Selector (MFSelector) based on a concept of total discriminating error (DEtotal) to identify monotonic genes. MFSelector considers various time stages in stage order (i.e., Stage One vs. other stages, Stages One and Two vs. remaining stages and so on) and computes DEtotal of each gene. MFSelector can successfully identify genes with monotonic characteristics.Results: We have demonstrated the effectiveness of MFSelector on two synthetic data sets and two stem cell differentiation data sets: embryonic stem cell neurogenesis (ESCN) and embryonic stem cell vasculogenesis (ESCV) data sets. We have also performed extensive quantitative comparisons of the three monotonic gene selection approaches. Some of the monotonic marker genes such as OCT4, NANOG, BLBP, discovered from the ESCN dataset exhibit consistent behavior with that reported in other studies. The role of monotonic genes found by MFSelector in either stemness or differentiation is validated using information obtained from Gene Ontology analysis and other literature. We justify and demonstrate that descending genes are involved in the proliferation or self-renewal activity of stem cells, while ascending genes are involved in differentiation of stem cells into variant cell lineages.Conclusions: We have developed a novel system, easy to use even with no pre-existing knowledge, to identify gene sets with monotonic expression patterns in multi-stage as well as in time-series genomics matrices. The case studies on ESCN and ESCV have helped to get a better understanding of stemness and differentiation. The novel monotonic marker genes discovered from a data set are found to exhibit consistent behavior in another independent data set, demonstrating the utility of the proposed method. The MFSelector R function and data sets can be downloaded from: http://microarray.ym.edu.tw/tools/MFSelector/

Missing value imputation in high-dimensional phenomic data: Imputable or not, and how?

Background: In modern biomedical research of complex diseases, a large number of demographic and clinical variables, herein called phenomic data, are often collected and missing values (MVs) are inevitable in the data collection process. Since many downstream statistical and bioinformatics methods require complete data matrix, imputation is a common and practical solution. In high-throughput experiments such as microarray experiments, continuous intensities are measured and many mature missing value imputation methods have been developed and widely applied. Numerous methods for missing data imputation of microarray data have been developed. Large phenomic data, however, contain continuous, nominal, binary and ordinal data types, which void application of most methods. Though several methods have been developed in the past few years, not a single complete guideline is proposed with respect to phenomic missing data imputation. Results: In this paper, we investigated existing imputation methods for phenomic data, proposed a self-training selection (STS) scheme to select the best imputation method and provide a practical guideline for general applications. We introduced a novel concept of "imputability measure" (IM) to identify missing values that are fundamentally inadequate to impute. In addition, we also developed four variations of K-nearest-neighbor (KNN) methods and compared with two existing methods, multivariate imputation by chained equations (MICE) and missForest. The four variations are imputation by variables (KNN-V), by subjects (KNN-S), their weighted hybrid (KNN-H) and an adaptively weighted hybrid (KNN-A). We performed simulations and applied different imputation methods and the STS scheme to three lung disease phenomic datasets to evaluate the methods. An R package "phenomeImpute" is made publicly available. Conclusions: Simulations and applications to real datasets showed that MICE often did not perform well, KNN-A, KNN-H and random forest were among the top performers although no method universally performed the best. Imputation of missing values with low imputability measures increased imputation errors greatly and could potentially deteriorate downstream analyses. The STS scheme was accurate in selecting the optimal method by evaluating methods in a second layer of missingness simulation. All source files for the simulation and the real data analyses are available on the author's publication website

Patient-provider communication regarding drug costsin Medicare Part D beneficiaries with diabetes: a TRIAD Study

<p>Abstract</p> <p>Background</p> <p>Little is known about drug cost communications of Medicare Part D beneficiaries with chronic conditions such as diabetes. The purpose of this study is to assess Medicare Part D beneficiaries with diabetes' levels of communication with physicians regarding prescription drug costs; the perceived importance of these communications; levels of prescription drug switching due to cost; and self-reported cost-related medication non-adherence.</p> <p>Methods</p> <p>Data were obtained from a cross-sectional survey (58% response rate) of 1,458 Medicare beneficiaries with diabetes who entered the coverage gap in 2006; adjusted percentages of patients with communication issues were obtained from multivariate regression analyses adjusting for patient demographics and clinical characteristics.</p> <p>Results</p> <p>Fewer than half of patients reported discussing the cost of medications with their physicians, while over 75% reported that such communications were important. Forty-eight percent reported their physician had switched to a less expensive medication due to costs. Minorities, females, and older adults had significantly lower levels of communication with their physicians regarding drug costs than white, male, and younger patients respectively. Patients with < $25 K annual household income were more likely than higher income patients to have talked about prescription drug costs with doctors, and to report cost-related non-adherence (27% vs. 17%, p < .001).</p> <p>Conclusions</p> <p>Medicare Part D beneficiaries with diabetes who entered the coverage gap have low levels of communication with physicians about drug costs, despite the high perceived importance of such communication. Understanding patient and plan-level characteristics differences in communication and use of cost-cutting strategies can inform interventions to help patients manage prescription drug costs.</p

Phytoplankton dynamics in relation to seasonal variability and upwelling and relaxation patterns at the mouth of Ria de Aveiro (West Iberian Margin) over a four-year period

From June 2004 to December 2007, samples were weekly collected at a fixed station located at the mouth of Ria de Aveiro (West Iberian Margin). We examined the seasonal and inter-annual fluctuations in composition and community structure of the phytoplankton in relation to the main environmental drivers and assessed the influence of the oceano-graphic regime, namely changes in frequency and intensity of upwelling events, over the dynamics of the phytoplankton assemblage. The samples were consistently handled and a final subset of 136 OTUs (taxa with relative abundance > 0.01%) was subsequently submitted to various multivariate analyses. The phytoplankton assemblage showed significant changes at all temporal scales but with an overriding importance of seasonality over longer-(inter-annual) or shorter-term fluctuations (upwelling-related). Sea-surface temperature, salinity and maximum upwelling index were retrieved as the main driver of seasonal change. Seasonal signal was most evident in the fluctuations of chlorophyll a concentration and in the high turnover from the winter to spring phytoplankton assemblage. The seasonal cycle of production and succession was disturbed by upwelling events known to disrupt thermal stratification and induce changes in the phytoplankton assemblage. Our results indicate that both the frequency and intensity of physical forcing were important drivers of such variability, but the outcome in terms of species composition was highly dependent on the available local pool of species and the timing of those events in relation to the seasonal cycle. We conclude that duration, frequency and intensity of upwelling events, which vary seasonally and inter-annually, are paramount for maintaining long-term phytoplankton diversity likely by allowing unstable coexistence and incorporating species turnover at different scales. Our results contribute to the understanding of the complex mechanisms of coastal phytoplankton dynamics in relation to changing physical forcing which is fundamental to improve predictability of future prospects under climate change.Portuguese Foundation for Science and Technology (FCT, Portugal) [SFRH/BPD/ 94562/2013]; FEDER funds; national funds; CESAM [UID/AMB/50017]; FCT/MEC through national funds; FEDERinfo:eu-repo/semantics/publishedVersio

Difference-based clustering of short time-course microarray data with replicates

<p>Abstract</p> <p>Background</p> <p>There are some limitations associated with conventional clustering methods for short time-course gene expression data. The current algorithms require prior domain knowledge and do not incorporate information from replicates. Moreover, the results are not always easy to interpret biologically.</p> <p>Results</p> <p>We propose a novel algorithm for identifying a subset of genes sharing a significant temporal expression pattern when replicates are used. Our algorithm requires no prior knowledge, instead relying on an observed statistic which is based on the first and second order differences between adjacent time-points. Here, a pattern is predefined as the sequence of symbols indicating direction and the rate of change between time-points, and each gene is assigned to a cluster whose members share a similar pattern. We evaluated the performance of our algorithm to those of K-means, Self-Organizing Map and the Short Time-series Expression Miner methods.</p> <p>Conclusions</p> <p>Assessments using simulated and real data show that our method outperformed aforementioned algorithms. Our approach is an appropriate solution for clustering short time-course microarray data with replicates.</p

Understanding Sensory Nerve Mechanotransduction through Localized Elastomeric Matrix Control

BACKGROUND: While neural systems are known to respond to chemical and electrical stimulation, the effect of mechanics on these highly sensitive cells is still not well understood. The ability to examine the effects of mechanics on these cells is limited by existing approaches, although their overall response is intimately tied to cell-matrix interactions. Here, we offer a novel method, which we used to investigate stretch-activated mechanotransduction on nerve terminals of sensory neurons through an elastomeric interface. METHODOLOGY/PRINCIPAL FINDINGS: To apply mechanical force on neurites, we cultured dorsal root ganglion neurons on an elastic substrate, polydimethylsiloxane (PDMS), coated with extracellular matrices (ECM). We then implemented a controlled indentation scheme using a glass pipette to mechanically stimulate individual neurites that were adjacent to the pipette. We used whole-cell patch clamping to record the stretch-activated action potentials on the soma of the single neurites to determine the mechanotransduction-based response. When we imposed specific mechanical force through the ECM, we noted a significant neuronal action potential response. Furthermore, because the mechanotransduction cascade is known to be directly affected by the cytoskeleton, we investigated the cell structure and its effects. When we disrupted microtubules and actin filaments with nocodozale or cytochalasin-D, respectively, the mechanically induced action potential was abrogated. In contrast, when using blockers of channels such as TRP, ASIC, and stretch-activated channels while mechanically stimulating the cells, we observed almost no change in action potential signalling when compared with mechanical activation of unmodified cells. CONCLUSIONS/SIGNIFICANCE: These results suggest that sensory nerve terminals have a specific mechanosensitive response that is related to cell architecture