Chien-Wei Lin

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.

Hypermethylation of BDNF and SST Genes in the Orbital Frontal Cortex of Older Individuals: A Putative Mechanism for Declining Gene Expression with Age

Expression of brain-derived neurotrophic factor (BDNF) and somatostatin (SST) mRNAs in the brain decreases progressively and robustly with age, and lower BDNF and SST expression in the brain has been observed in many brain disorders. BDNF is known to regulate SST expression; however, the mechanisms underlying decreased expression of both genes are not understood. DNA methylation (DNAm) is an attractive candidate mechanism. To investigate the contribution of DNAm to the age-related decline in BDNF and SST expression, the Illumina Infinium HumanMethylation450 Beadchip Array was used to quantify DNAm of BDNF (26 CpG loci) and SST (9 CpG loci) in the orbital frontal cortices of postmortem brains from 22 younger (age <42 years) and 22 older individuals (age >60 years) with known age-dependent BDNF and SST expression differences. Relative to the younger individuals, 10 of the 26 CpG loci in BDNF and 8 of the 9 CpG loci in SST were significantly hypermethylated in the older individuals. DNAm in BDNF exons/promoters I, II, and IV negatively correlated with BDNF expression (r=−0.37, p<0.05; r=−0.40, p<0.05; r=−0.24, p=0.07), and DNAm in SST 5′ UTR and first exon/intron negatively correlated with SST expression (r=−0.48, p<0.01; r=−0.63, p<0.001), respectively. An expanded set of BDNF- and GABA-related genes exhibited similar age-related changes in DNAm and correlation with gene expression. These results suggest that DNAm may be a proximal mechanism for decreased expression of BDNF, SST, and other BDNF- and GABA-related genes with brain aging and, by extension, for brain disorders in which their expression is decreased.

SMART: Statistical Metabolomics Analysis—An R Tool

Metabolomics data provide unprecedented opportunities to decipher metabolic mechanisms by analyzing hundreds to thousands of metabolites. Data quality concerns and complex batch effects in metabolomics must be appropriately addressed through statistical analysis. This study developed an integrated analysis tool for metabolomics studies to streamline the complete analysis flow from initial data preprocessing to downstream association analysis. We developed Statistical Metabolomics Analysis-An R Tool (SMART), which can analyze input files with different formats, visually represent various types of data features, implement peak alignment and annotation, conduct quality control for samples and peaks, explore batch effects, and perform association analysis. A pharmacometabolomics study of antihypertensive medication was conducted and data were analyzed using SMART. Neuromedin N was identified as a metabolite significantly associated with angiotensin-converting-enzyme inhibitors in our metabolome-wide association analysis (p = 1.56 × 10(-4) in an analysis of covariance (ANCOVA) with an adjustment for unknown latent groups and p = 1.02 × 10(-4) in an ANCOVA with an adjustment for hidden substructures). This endogenous neuropeptide is highly related to neurotensin and neuromedin U, which are involved in blood pressure regulation and smooth muscle contraction. The SMART software, a user guide, and example data can be downloaded from http://www.stat.sinica.edu.tw/hsinchou/metabolomics/SMART.htm .

Predisposition to treatment response in major depressive episode: A peripheral blood gene coexpression network analysis

Antidepressant efficacy is insufficient, unpredictable and poorly understood in major depressive episode (MDE). Gene expression studies allow for the identification of significantly dysregulated genes but can limit the exploration of biological pathways. In the present study, we proposed a gene coexpression analysis to investigate biological pathways associated with treatment response predisposition and their regulation by microRNAs (miRNAs) in peripheral blood samples of MDE and healthy control subjects. We used a discovery cohort that included 34 MDE patients that were given 12-week treatment with citalopram and 33 healthy controls. Two replication cohorts with similar design were also analyzed. Expression-based gene network was built to define clusters of highly correlated sets of genes, called modules. Association between each modules first principal component of the expression data and clinical improvement was tested in the three cohorts. We conducted gene ontology analysis and miRNA prediction based on the module gene list. Nine of the 59 modules from the gene coexpression network were associated with clinical improvement. The association was partially replicated in other cohorts. Gene ontology analysis demonstrated that 4 modules were associated with cytokine production, acute inflammatory response or IL-8 functions. Finally, we found 414 miRNAs that may regulate one or several modules associated with clinical improvement. By contrast, only 12 miRNAs were predicted to specifically regulate modules unrelated to clinical improvement. Our gene coexpression analysis underlines the importance of inflammation-related pathways and the involvement of a large miRNA program as biological processes predisposing associated with antidepressant response.

VSNL1 Co-Expression Networks in Aging Include Calcium Signaling, Synaptic Plasticity, and Alzheimer’s Disease Pathways

The visinin-like 1 (VSNL1) gene encodes visinin-like protein 1, a peripheral biomarker for Alzheimer disease (AD). Little is known, however, about normal VSNL1 expression in brain and the biologic networks in which it participates. Frontal cortex gray matter obtained from 209 subjects without neurodegenerative or psychiatric illness, ranging in age from 16 to 91, was processed on Affymetrix GeneChip 1.1 ST and Human SNP Array 6.0. VSNL1 expression was unaffected by age and sex, and not significantly associated with SNPs in cis or trans. VSNL1 was significantly co-expressed with genes in pathways for calcium signaling, AD, long-term potentiation, long-term depression, and trafficking of AMPA receptors. The association with AD was driven, in part, by correlation with amyloid precursor protein (APP) expression. These findings provide an unbiased link between VSNL1 and molecular mechanisms of AD, including pathways implicated in synaptic pathology in AD. Whether APP may drive increased VSNL1 expression, VSNL1 drives increased APP expression, or both are downstream of common pathogenic regulators will need to be evaluated in model systems.

MetaKTSP: A Meta-Analytic Top Scoring Pair Method for Robust Cross-Study Validation of Omics Prediction Analysis

MOTIVATION Supervised machine learning is widely applied to transcriptomic data to predict disease diagnosis, prognosis or survival. Robust and interpretable classifiers with high accuracy are usually favored for their clinical and translational potential. The top scoring pair (TSP) algorithm is an example that applies a simple rank-based algorithm to identify rank-altered gene pairs for classifier construction. Although many classification methods perform well in cross-validation of single expression profile, the performance usually greatly reduces in cross-study validation (i.e. the prediction model is established in the training study and applied to an independent test study) for all machine learning methods, including TSP. The failure of cross-study validation has largely diminished the potential translational and clinical values of the models. The purpose of this article is to develop a meta-analytic top scoring pair (MetaKTSP) framework that combines multiple transcriptomic studies and generates a robust prediction model applicable to independent test studies. RESULTS We proposed two frameworks, by averaging TSP scores or by combining P-values from individual studies, to select the top gene pairs for model construction. We applied the proposed methods in simulated data sets and three large-scale real applications in breast cancer, idiopathic pulmonary fibrosis and pan-cancer methylation. The result showed superior performance of cross-study validation accuracy and biomarker selection for the new meta-analytic framework. In conclusion, combining multiple omics data sets in the public domain increases robustness and accuracy of the classification model that will ultimately improve disease understanding and clinical treatment decisions to benefit patients. AVAILABILITY AND IMPLEMENTATION An R package MetaKTSP is available online. (http://tsenglab.biostat.pitt.edu/software.htm). CONTACT ctseng@pitt.edu SUPPLEMENTARY INFORMATION Supplementary data are available at Bioinformatics online.

Enhanced Molecular Aging in Late-Life Depression: the Senescent-Associated Secretory Phenotype

OBJECTIVE This study aims to investigate whether a systemic molecular pattern associated with aging (senescent-associated secretory phenotype [SASP]) is elevated in adults with late-life depression (LLD), compared with never-depressed elderly comparison participants. DESIGN Cross-sectional study. PARTICIPANTS We included 111 older adults (80 with LLD and 31 comparison participants) in this study. MEASUREMENT A panel of 22 SASP-related proteins was extracted from a previous multiplex protein panel performed in these participants. We conducted a principal component analysis to create the SASP index based on individual weights of each of protein. RESULTS Participants with LLD showed a significantly increased SASP index compared with comparison participants, after controlling for age, depressive symptoms, medical comorbidity (CIRS-G) scores, sex, and cognitive performance (F(1,98) = 7.3, p = 0.008). Correlation analyses revealed that the SASP index was positively correlated with age (r = 0.2, p = 0.03) and CIRS score (r = 0.27, p = 0.005), and negatively correlated with information processing speed (r = -0.34, p = 0.001), executive function (r = -0.27, p = 0.004) and global cognitive performance (r = -0.28, p = 0.007). CONCLUSIONS To the best of our knowledge, this is the first study to show that a set of proteins (i.e., SASP index) primarily associated with cellular aging is abnormally regulated and elevated in LLD. These results suggest that individuals with LLD display enhanced aging-related molecular patterns that are associated with higher medical comorbidity and worse cognitive function. Finally, we provide a set of proteins that can serve as potential therapeutic targets and biomarkers to monitor the effects of therapeutic or preventative interventions in LLD.

Integrative analysis of single nucleotide polymorphisms and gene expression efficiently distinguishes samples from closely related ethnic populations

BackgroundAncestry informative markers (AIMs) are a type of genetic marker that is informative for tracing the ancestral ethnicity of individuals. Application of AIMs has gained substantial attention in population genetics, forensic sciences, and medical genetics. Single nucleotide polymorphisms (SNPs), the materials of AIMs, are useful for classifying individuals from distinct continental origins but cannot discriminate individuals with subtle genetic differences from closely related ancestral lineages. Proof-of-principle studies have shown that gene expression (GE) also is a heritable human variation that exhibits differential intensity distributions among ethnic groups. GE supplies ethnic information supplemental to SNPs; this motivated us to integrate SNP and GE markers to construct AIM panels with a reduced number of required markers and provide high accuracy in ancestry inference. Few studies in the literature have considered GE in this aspect, and none have integrated SNP and GE markers to aid classification of samples from closely related ethnic populations.ResultsWe integrated a forward variable selection procedure into flexible discriminant analysis to identify key SNP and/or GE markers with the highest cross-validation prediction accuracy. By analyzing genome-wide SNP and/or GE markers in 210 independent samples from four ethnic groups in the HapMap II Project, we found that average testing accuracies for a majority of classification analyses were quite high, except for SNP-only analyses that were performed to discern study samples containing individuals from two close Asian populations. The average testing accuracies ranged from 0.53 to 0.79 for SNP-only analyses and increased to around 0.90 when GE markers were integrated together with SNP markers for the classification of samples from closely related Asian populations. Compared to GE-only analyses, integrative analyses of SNP and GE markers showed comparable testing accuracies and a reduced number of selected markers in AIM panels.ConclusionsIntegrative analysis of SNP and GE markers provides high-accuracy and/or cost-effective classification results for assigning samples from closely related or distantly related ancestral lineages to their original ancestral populations. User-friendly BIASLESS (B iomarkers I dentification a nd S amples S ubdivision) software was developed as an efficient tool for selecting key SNP and/or GE markers and then building models for sample subdivision. BIASLESS was programmed in R and R-GUI and is available online at http://www.stat.sinica.edu.tw/hsinchou/genetics/prediction/BIASLESS.htm.

An Analytical Approach to Optimizing the Utility of ESP Games

In this paper, we propose an analytical model for computing the utility of ESP games, i.e., the throughput rate of appropriate labels for given puzzles. The model targets generalized games, where the number of players, the consensus threshold, and the stopping condition are variable. Via extensive simulations, we show that our model can accurately predict the stopping condition that will yield the optimal utility of an ESP game under a specific setting. A service provider can therefore utilize the model to ensure that the hosted ESP games produce high-quality labels efficiently, given that the number of players willing to invest time and effort in the game is limited.

SAQC: SNP array quality control.

BackgroundGenome-wide single-nucleotide polymorphism (SNP) arrays containing hundreds of thousands of SNPs from the human genome have proven useful for studying important human genome questions. Data quality of SNP arrays plays a key role in the accuracy and precision of downstream data analyses. However, good indices for assessing data quality of SNP arrays have not yet been developed.ResultsWe developed new quality indices to measure the quality of SNP arrays and/or DNA samples and investigated their statistical properties. The indices quantify a departure of estimated individual-level allele frequencies (AFs) from expected frequencies via standardized distances. The proposed quality indices followed lognormal distributions in several large genomic studies that we empirically evaluated. AF reference data and quality index reference data for different SNP array platforms were established based on samples from various reference populations. Furthermore, a confidence interval method based on the underlying empirical distributions of quality indices was developed to identify poor-quality SNP arrays and/or DNA samples. Analyses of authentic biological data and simulated data show that this new method is sensitive and specific for the detection of poor-quality SNP arrays and/or DNA samples.ConclusionsThis study introduces new quality indices, establishes references for AFs and quality indices, and develops a detection method for poor-quality SNP arrays and/or DNA samples. We have developed a new computer program that utilizes these methods called SNP Array Quality Control (SAQC). SAQC software is written in R and R-GUI and was developed as a user-friendly tool for the visualization and evaluation of data quality of genome-wide SNP arrays. The program is available online (http://www.stat.sinica.edu.tw/hsinchou/genetics/quality/SAQC.htm).