Integrating population genetics, stem cell biology and cellular genomics to study complex human diseases.

Human pluripotent stem (hPS) cells can, in theory, be differentiated into any cell type, making them a powerful in vitro model for human biology. Recent technological advances have facilitated large-scale hPS cell studies that allow investigation of the genetic regulation of molecular phenotypes and their contribution to high-order phenotypes such as human disease. Integrating hPS cells with single-cell sequencing makes identifying context-dependent genetic effects during cell development or upon experimental manipulation possible.

Molecular profiling of human substantia nigra identifies diverse neuron types associated with vulnerability in Parkinson's disease.

Parkinson's disease (PD) is characterized pathologically by the loss of dopaminergic (DA) neurons in the substantia nigra (SN). Whether cell types beyond DA neurons in the SN show vulnerability in PD remains unclear. Through transcriptomic profiling of 315,867 high-quality single nuclei in the SN from individuals with and without PD, we identified cell clusters representing various neuron types, glia, endothelial cells, pericytes, fibroblasts, and T cells and investigated cell type-dependent alterations in gene expression in PD.

Novel human pluripotent stem cell-derived hypothalamus organoids demonstrate cellular diversity.

The hypothalamus is a region of the brain that plays an important role in regulating body functions and behaviors. There is a growing interest in human pluripotent stem cells (hPSCs) for modeling diseases that affect the hypothalamus. Here, we established an hPSC-derived hypothalamus organoid differentiation protocol to model the cellular diversity of this brain region.

A village in a dish model system for population-scale hiPSC studies.

The mechanisms by which DNA alleles contribute to disease risk, drug response, and other human phenotypes are highly context-specific, varying across cell types and different conditions. Human induced pluripotent stem cells are uniquely suited to study these context-dependent effects but cell lines from hundreds or thousands of individuals are required. Village cultures, where multiple induced pluripotent stem lines are cultured and differentiated in a single dish, provide an elegant solution for scaling induced pluripotent stem experiments to the necessary sample sizes required for population-scale studies.

Pitfalls and opportunities for applying latent variables in single-cell eQTL analyses.

Using latent variables in gene expression data can help correct unobserved confounders and increase statistical power for expression quantitative trait Loci (eQTL) detection. The probabilistic estimation of expression residuals (PEER) and principal component analysis (PCA) are widely used methods that can remove unwanted variation and improve eQTL discovery power in bulk RNA-seq analysis. However, their performance has not been evaluated extensively in single-cell eQTL analysis, especially for different cell types.

Shared genetic susceptibility between trigger finger and carpal tunnel syndrome: a genome-wide association study.

Background: Trigger finger and carpal tunnel syndrome are the two most common non-traumatic connective tissue disorders of the hand. Both of these conditions frequently co-occur, often in patients with rheumatoid arthritis. However, this phenotypic association is poorly understood.

Transcriptomic and proteomic retinal pigment epithelium signatures of age-related macular degeneration.

There are currently no treatments for geographic atrophy, the advanced form of age-related macular degeneration. Hence, innovative studies are needed to model this condition and prevent or delay its progression. Induced pluripotent stem cells generated from patients with geographic atrophy and healthy individuals were differentiated to retinal pigment epithelium.

Single cell eQTL analysis identifies cell type-specific genetic control of gene expression in fibroblasts and reprogrammed induced pluripotent stem cells.

Background: The discovery that somatic cells can be reprogrammed to induced pluripotent stem cells (iPSCs) has provided a foundation for in vitro human disease modelling, drug development and population genetics studies. Gene expression plays a critical role in complex disease risk and therapeutic response. However, while the genetic background of reprogrammed cell lines has been shown to strongly influence gene expression, the effect has not been evaluated at the level of individual cells which would provide significant resolution.

Prediction of short-term antidepressant response using probabilistic graphical models with replication across multiple drugs and treatment settings.

Heterogeneity in the clinical presentation of major depressive disorder and response to antidepressants limits clinicians' ability to accurately predict a specific patient's eventual response to therapy. Validated depressive symptom profiles may be an important tool for identifying poor outcomes early in the course of treatment. To derive these symptom profiles, we first examined data from 947 depressed subjects treated with selective serotonin reuptake inhibitors (SSRIs) to delineate the heterogeneity of antidepressant response using probabilistic graphical models (PGMs).

Acylcarnitine metabolomic profiles inform clinically-defined major depressive phenotypes.

Background: Acylcarnitines have important functions in mitochondrial energetics and β-oxidation, and have been implicated to play a significant role in metabolic functions of the brain. This retrospective study examined whether plasma acylcarnitine profiles can help biochemically distinguish the three phenotypic subtypes of major depressive disorder (MDD): core depression (CD+), anxious depression (ANX+), and neurovegetative symptoms of melancholia (NVSM+).

Methods: Depressed outpatients (n = 240) from the Mayo Clinic Pharmacogenomics Research Network were treated with citalopram or escitalopram for eight weeks.

Dual Roles for the TSPYL Family in Mediating Serotonin Transport and the Metabolism of Selective Serotonin Reuptake Inhibitors in Patients with Major Depressive Disorder.

We previously reported that testis-specific Y-encoded-like protein (TSPYLs) are transcription regulators for CYP3A4, CYP2C9, and CYP2C19. Here, we observed dual roles for TSPYLs in mediating serotonin transport and the metabolism of selective serotonin reuptake inhibitors (SSRIs) in patients with major depressive disorder (MDD). The widely prescribed SSRIs, citalopram, and escitalopram are metabolized mainly by CYP2C19.

Single Nucleotide Polymorphisms at a Distance from Aryl Hydrocarbon Receptor (AHR) Binding Sites Influence AHR Ligand-Dependent Gene Expression.

Greater than 90% of significant genome-wide association study (GWAS) single-nucleotide polymorphisms (SNPs) are in noncoding regions of the genome, but only 25.6% are known expression quantitative trait loci (eQTLs). Therefore, the function of many significant GWAS SNPs remains unclear.

Metabolomic signature of exposure and response to citalopram/escitalopram in depressed outpatients.

Metabolomics provides valuable tools for the study of drug effects, unraveling the mechanism of action and variation in response due to treatment. In this study we used electrochemistry-based targeted metabolomics to gain insights into the mechanisms of action of escitalopram/citalopram focusing on a set of 31 metabolites from neurotransmitter-related pathways. Overall, 290 unipolar patients with major depressive disorder were profiled at baseline, after 4 and 8 weeks of drug treatment.

Pharmacogenomics-Driven Prediction of Antidepressant Treatment Outcomes: A Machine-Learning Approach With Multi-trial Replication.

We set out to determine whether machine learning-based algorithms that included functionally validated pharmacogenomic biomarkers joined with clinical measures could predict selective serotonin reuptake inhibitor (SSRI) remission/response in patients with major depressive disorder (MDD). We studied 1,030 white outpatients with MDD treated with citalopram/escitalopram in the Mayo Clinic Pharmacogenomics Research Network Antidepressant Medication Pharmacogenomic Study (PGRN-AMPS; n = 398), Sequenced Treatment Alternatives to Relieve Depression (STAR*D; n = 467), and International SSRI Pharmacogenomics Consortium (ISPC; n = 165) trials. A genomewide association study for PGRN-AMPS plasma metabolites associated with SSRI response (serotonin) and baseline MDD severity (kynurenine) identified single nucleotide polymorphisms (SNPs) in DEFB1, ERICH3, AHR, and TSPAN5 that we tested as predictors.

Pharmacogenomic Next-Generation DNA Sequencing: Lessons from the Identification and Functional Characterization of Variants of Unknown Significance in and .

and are highly polymorphic pharmacogenes; however, clinically actionable genetic variability in drug metabolism due to these genes has been limited to a few common alleles. The identification and functional characterization of less-common open reading frame sequence variation might help to individualize therapy with drugs that are substrates for the enzymes encoded by these genes. The present study identified seven uncharacterized variants each in and using next-generation sequence data for 1013 subjects, and functionally characterized the encoded proteins.

The Role of the Aryl Hydrocarbon Receptor (AHR) in Immune and Inflammatory Diseases.

The aryl hydrocarbon receptor (AHR) is a nuclear receptor that modulates the response to environmental stimuli. It was recognized historically for its role in toxicology but, in recent decades, it has been increasingly recognized as an important modulator of disease-especially for its role in modulating immune and inflammatory responses. AHR has been implicated in many diseases that are driven by immune/inflammatory processes, including major depressive disorder, multiple sclerosis, rheumatoid arthritis, asthma, and allergic responses, among others.

Augmentation of Physician Assessments with Multi-Omics Enhances Predictability of Drug Response: A Case Study of Major Depressive Disorder.

This work proposes a "" workflow to sequentially augment physician assessments of patients' symptoms and their socio-demographic measures with heterogeneous biological measures to accurately predict treatment outcomes using machine learning. Across many psychiatric illnesses, ranging from major depressive disorder to schizophrenia, symptom severity assessments are subjective and do not include biological measures, making predictability in eventual treatment outcomes a challenge. Using data from the Mayo Clinic PGRN-AMPS SSRI trial as a case study, this work demonstrates a significant improvement in the prediction accuracy for antidepressant treatment outcomes in patients with major depressive disorder from 35% to 80% individualized by patient, compared to using only a physician's assessment as the predictors.

Beta-defensin 1, aryl hydrocarbon receptor and plasma kynurenine in major depressive disorder: metabolomics-informed genomics.

Major depressive disorder (MDD) is a heterogeneous disease. Efforts to identify biomarkers for sub-classifying MDD and antidepressant therapy by genome-wide association studies (GWAS) alone have generally yielded disappointing results. We applied a metabolomics-informed genomic research strategy to study the contribution of genetic variation to MDD pathophysiology by assaying 31 metabolites, including compounds from the tryptophan, tyrosine, and purine pathways, in plasma samples from 290 MDD patients.

Pharmacometabolomics informs Pharmacogenomics.

Introduction: The initial decades of the 21 century have witnessed striking technical advances that have made it possible to detect, identify and quantitatively measure large numbers of plasma or tissue metabolites. In parallel, similar advances have taken place in our ability to sequence DNA and RNA. Those advances have moved us beyond studies of single metabolites and single genetic polymorphisms to the study of hundreds or thousands of metabolites and millions of genomic variants in a single cell or subject.