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. Here we discuss how the intersection of stem cell biology, population genetics and cellular genomics can help resolve the functional consequences of human genetic variation. We examine the critical challenges of integrating these fields and approaches to scaling them cost-effectively and practically. We highlight two areas of human biology that can particularly benefit from population-scale hPS cell studies, elucidating mechanisms underlying complex disease risk loci and evaluating relationships between common genetic variation and pharmacotherapeutic phenotypes.

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. Using an hPSC line with a tyrosine hydroxylase (TH)-TdTomato reporter for dopaminergic neurons (DNs) and other TH-expressing cells, we interrogated DN-specific pathways and functions in electrophysiologically active hypothalamus organoids. Single-cell RNA sequencing (scRNA-seq) revealed diverse neuronal and non-neuronal cell types in mature hypothalamus organoids. We identified several molecularly distinct hypothalamic DN subtypes that demonstrated different developmental maturities. Our in vitro 3D hypothalamus differentiation protocol can be used to study the development of this critical brain structure and can be applied to disease modeling to generate novel therapeutic approaches for disorders centered around the hypothalamus.

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. Here, we show the utility of village models, demonstrating how cells can be assigned to an induced pluripotent stem line using single-cell sequencing and illustrating that the genetic, epigenetic or induced pluripotent stem line-specific effects explain a large percentage of gene expression variation for many genes. We demonstrate that village methods can effectively detect induced pluripotent stem line-specific effects, including sensitive dynamics of cell states.

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. We have identified a novel type of eQTL for which SNPs distant from ligand-activated transcription factor (TF) binding sites can alter target gene expression in a SNP genotype-by-ligand-dependent fashion that we refer to as pharmacogenomic eQTLs (PGx-eQTLs)-loci that may have important pharmacotherapeutic implications. In the present study, we integrated chromatin immunoprecipitation-seq with RNA-seq and SNP genotype data for a panel of lymphoblastoid cell lines to identify 10 novel cis PGx-eQTLs dependent on the ligand-activated TF aryl hydrocarbon receptor (AHR)-a critical environmental sensor for xenobiotic (drug) and immune response. Those 10 cis PGx-eQTLs were eQTLs only after AHR ligand treatment, even though the SNPs did not create/destroy an AHR response element-the DNA sequence motif recognized and bound by AHR. Additional functional studies in multiple cell lines demonstrated that some cis PGx-eQTLs are functional in multiple cell types, whereas others displayed SNP-by-ligand-dependent effects in just one cell type. Furthermore, four of those cis PGx-eQTLs had previously been associated with clinical phenotypes, indicating that those loci might have the potential to inform clinical decisions. Therefore, SNPs across the genome that are distant from TF binding sites for ligand-activated TFs might function as PGx-eQTLs and, as a result, might have important clinical implications for interindividual variation in drug response. SIGNIFICANCE STATEMENT: More than 90% of single-nucleotide polymorphisms (SNPs) that are associated with clinical phenotypes are located in noncoding regions of the genome. However, the mechanisms of action of many of those SNPs have not been elucidated, and drugs may unmask functional expression quantitative trail loci (eQTLs). In the current study, we used drugs that bind to the ligand-activated transcription factor aryl hydrocarbon receptor (AHR) and identified SNPs that were associated with interindividual variation in gene expression following drug exposure-termed pharmacogenomic (PGx)-eQTLs. Possibly of greater significance, those PGx-eQTL SNPs were outside of AHR binding sites, indicating that they do not interrupt AHR DNA recognition. PGx-eQTLs such as those described in this work may have crucial implications for interindividual variation in drug.

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

CYP2C9 and CYP2C19 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 CYP2C9 and CYP2C19 using next-generation sequence data for 1013 subjects, and functionally characterized the encoded proteins. Constructs were created and transiently expressed in COS-1 cells for the assay of protein concentration and enzyme activities using fluorometric substrates and liquid chromatography- tandem mass spectrometry with tolbutamide (CYP2C9) and (S)-mephenytoin (CYP2C19) as prototypic substrates. The results were compared with the SIFT, Polyphen, and Provean functional prediction software programs. Cytochrome P450 oxidoreductase (CPR) activities were also determined. Positive correlations were observed between protein content and fluorometric enzyme activity for variants of CYP2C9 (P < 0.05) and CYP2C19 (P < 0.0005). However, CYP2C9 709G>C and CYP2C19 65A>G activities were much lower than predicted based on protein content. Substrate intrinsic clearance values for CYP2C9 218C>T, 343A>C, and CYP2C19 337G>A, 518C>T, 556C>T, and 557G>A were less than 25% of wild-type allozymes. CPR activity levels were similar for all variants. In summary, sequencing of CYP2C9 and CYP2C19 in 1013 subjects identified low-frequency variants that had not previously been functionally characterized. In silico predictions were not always consistent with functional assay results. These observations emphasize the need for high-throughput methods for pharmacogene variant mutagenesis and functional characterization.

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. The mechanisms by which AHR has been suggested to impact immune/inflammatory diseases include targeted gene expression and altered immune differentiation. It has been suggested that single nucleotide polymorphisms (SNPs) that are near AHR-regulated genes may contribute to AHR-dependent disease mechanisms/pathways. Further, we have found that SNPs that are outside of nuclear receptor binding sites (i.e., outside of AHR response elements (AHREs)) may contribute to AHR-dependent gene regulation in a SNP- and ligand-dependent manner. This review will discuss the evidence and mechanisms of AHR contributions to immune/inflammatory diseases and will consider the possibility that SNPs that are outside of AHR binding sites might contribute to AHR ligand-dependent inter-individual variation in disease pathophysiology and response to pharmacotherapeutics.

Treatment of Major Depressive Disorder in Pediatric Populations.

Major depressive disorder (MDD) is a severe illness that afflicts about 16.6% of people over their lifetime. MDD is highly correlated with suicidality, and often first presents in adolescence. Unfortunately, many pediatric patients suffering from MDD go undiagnosed, and current evidence-based treatment options in the U.S. are limited to psychotherapy and two selective serotonin reuptake inhibitors approved by the United States Food and Drug Administration. Molecular mechanisms have been shown to play a role in MDD pathogenesis, progression, and response to medication, yet few studies have explored the role of these pathways in pediatric MDD. In this review, we outline the gravity and importance of MDD in pediatric patients, some challenges in diagnosis and treatment, current treatments available for pediatric patients, and research to investigate differences between pediatric and adult MDD. We hope that this review will provide an outline of the current understanding and treatment of MDD in pediatric patients, and provide thoughtful insights for future work that could advance our understanding of MDD in pediatric populations, and also identify new therapeutic strategies.

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. Associations of metabolite concentrations with depressive symptoms were determined, followed by GWAS for selected metabolites and functional validation studies of the genes identified. Kynurenine (KYN), the baseline plasma metabolite that was most highly associated with depressive symptoms, was negatively correlated with severity of those symptoms. GWAS for baseline plasma KYN concentrations identified SNPs across the beta-defensin 1 (DEFB1) and aryl hydrocarbon receptor (AHR) genes that were cis-expression quantitative trait loci (eQTLs) for DEFB1 and AHR mRNA expression, respectively. Furthermore, the DEFB1 locus was associated with severity of MDD symptoms in a larger cohort of 803 MDD patients. Functional studies demonstrated that DEFB1 could neutralize lipopolysaccharide-stimulated expression of KYN-biosynthesizing enzymes in monocytic cells, resulting in altered KYN concentrations in the culture media. In addition, we demonstrated that AHR was involved in regulating the expression of enzymes in the KYN pathway and altered KYN biosynthesis in cell lines of hepatocyte and astrocyte origin. In conclusion, these studies identified SNPs that were cis-eQTLs for DEFB1 and AHR and, which were associated with variation in plasma KYN concentrations that were related to severity of MDD symptoms.