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.

Integrating Computational and Biological Hemodynamic Approaches to Improve Modeling of Atherosclerotic Arteries.

Atherosclerosis is the primary cause of cardiovascular disease, resulting in mortality, elevated healthcare costs, diminished productivity, and reduced quality of life for individuals and their communities. This is exacerbated by the limited understanding of its underlying causes and limitations in current therapeutic interventions, highlighting the need for sophisticated models of atherosclerosis. This review critically evaluates the computational and biological models of atherosclerosis, focusing on the study of hemodynamics in atherosclerotic coronary arteries. Computational models account for the geometrical complexities and hemodynamics of the blood vessels and stenoses, but they fail to capture the complex biological processes involved in atherosclerosis. Different in vitro and in vivo biological models can capture aspects of the biological complexity of healthy and stenosed vessels, but rarely mimic the human anatomy and physiological hemodynamics, and require significantly more time, cost, and resources. Therefore, emerging strategies are examined that integrate computational and biological models, and the potential of advances in imaging, biofabrication, and machine learning is explored in developing more effective models of atherosclerosis.

Deep sequencing of proteotoxicity modifier genes uncovers a Presenilin-2/beta-amyloid-actin genetic risk module shared among alpha-synucleinopathies.

Whether neurodegenerative diseases linked to misfolding of the same protein share genetic risk drivers or whether different protein-aggregation pathologies in neurodegeneration are mechanistically related remains uncertain. Conventional genetic analyses are underpowered to address these questions. Through careful selection of patients based on protein aggregation phenotype (rather than clinical diagnosis) we can increase statistical power to detect associated variants in a targeted set of genes that modify proteotoxicities. Genetic modifiers of alpha-synuclein (ɑS) and beta-amyloid (Aß) cytotoxicity in yeast are enriched in risk factors for Parkinson's disease (PD) and Alzheimer's disease (AD), respectively. Here, along with known AD/PD risk genes, we deeply sequenced exomes of 430 ɑS/Aß modifier genes in patients across alpha-synucleinopathies (PD, Lewy body dementia and multiple system atrophy). Beyond known PD genes GBA1 and LRRK2, rare variants AD genes (CD33, CR1 and PSEN2) and Aß toxicity modifiers involved in RhoA/actin cytoskeleton regulation (ARGHEF1, ARHGEF28, MICAL3, PASK, PKN2, PSEN2) were shared risk factors across synucleinopathies. Actin pathology occurred in iPSC synucleinopathy models and RhoA downregulation exacerbated ɑS pathology. Even in sporadic PD, the expression of these genes was altered across CNS cell types. Genome-wide CRISPR screens revealed the essentiality of PSEN2 in both human cortical and dopaminergic neurons, and PSEN2 mutation carriers exhibited diffuse brainstem and cortical synucleinopathy independent of AD pathology. PSEN2 contributes to a common-risk signal in PD GWAS and regulates ɑS expression in neurons. Our results identify convergent mechanisms across synucleinopathies, some shared with AD.

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.

Mimicry of embryonic circulation enhances the hoxa hemogenic niche and human blood development.

Precursors of the adult hematopoietic system arise from the aorta-gonad-mesonephros (AGM) region shortly after the embryonic circulation is established. Here, we develop a microfluidic culture system to mimic the primitive embryonic circulation and address the hypothesis that circulatory flow and shear stress enhance embryonic blood development. Embryonic (HOXA+) hematopoiesis was derived from human pluripotent stem cells and induced from mesoderm by small-molecule manipulation of TGF-ß and WNT signaling (SB/CHIR). Microfluidic and orbital culture promoted the formation of proliferative CD34+RUNX1C-GFP+SOX17-mCHERRY+ precursor cells that were released into the artificial circulation from SOX17+ arterial-like structures. Single-cell transcriptomic analysis delineated extra-embryonic (yolk sac) and HOXA+ embryonic blood differentiation pathways. SB/CHIR and circulatory flow enhance hematopoiesis by the formation of proliferative HOXA+RUNX1C+CD34+ precursor cells that differentiate into monocyte/macrophage, granulocyte, erythrocyte, and megakaryocyte progenitors.

Hif-1a suppresses ROS-induced proliferation of cardiac fibroblasts following myocardial infarction.

We report that cardiac fibroblasts (CFs) and mesenchymal progenitors are more hypoxic than other cardiac interstitial populations, express more hypoxia-inducible factor 1α (HIF-1α), and exhibit increased glycolytic metabolism. CF-specific deletion of Hif-1a resulted in decreased HIF-1 target gene expression and increased mesenchymal progenitors in uninjured hearts and increased CF activation without proliferation following sham injury, as demonstrated using single-cell RNA sequencing (scRNA-seq). After myocardial infarction (MI), however, there was ∼50% increased CF proliferation and excessive scarring and contractile dysfunction, a scenario replicated in 3D engineered cardiac microtissues. CF proliferation was associated with higher reactive oxygen species (ROS) as occurred also in wild-type mice treated with the mitochondrial ROS generator MitoParaquat (MitoPQ). The mitochondrial-targeted antioxidant MitoTEMPO rescued Hif-1a mutant phenotypes. Thus, HIF-1α in CFs provides a critical braking mechanism against excessive post-ischemic CF activation and proliferation through regulation of mitochondrial ROS. CFs are potential cellular targets for designer antioxidant therapies in cardiovascular disease.

Retinal ganglion cell-specific genetic regulation in primary open-angle glaucoma.

To assess the transcriptomic profile of disease-specific cell populations, fibroblasts from patients with primary open-angle glaucoma (POAG) were reprogrammed into induced pluripotent stem cells (iPSCs) before being differentiated into retinal organoids and compared with those from healthy individuals. We performed single-cell RNA sequencing of a total of 247,520 cells and identified cluster-specific molecular signatures. Comparing the gene expression profile between cases and controls, we identified novel genetic associations for this blinding disease. Expression quantitative trait mapping identified a total of 4,443 significant loci across all cell types, 312 of which are specific to the retinal ganglion cell subpopulations, which ultimately degenerate in POAG. Transcriptome-wide association analysis identified genes at loci previously associated with POAG, and analysis, conditional on disease status, implicated 97 statistically significant retinal ganglion cell-specific expression quantitative trait loci. This work highlights the power of large-scale iPSC studies to uncover context-specific profiles for a genetically complex disease.

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. Integrating transcriptional profiles of 127,659 retinal pigment epithelium cells generated from 43 individuals with geographic atrophy and 36 controls with genotype data, we identify 445 expression quantitative trait loci in cis that are asssociated with disease status and specific to retinal pigment epithelium subpopulations. Transcriptomics and proteomics approaches identify molecular pathways significantly upregulated in geographic atrophy, including in mitochondrial functions, metabolic pathways and extracellular cellular matrix reorganization. Five significant protein quantitative trait loci that regulate protein expression in the retinal pigment epithelium and in geographic atrophy are identified - two of which share variants with cis- expression quantitative trait loci, including proteins involved in mitochondrial biology and neurodegeneration. Investigation of mitochondrial metabolism confirms mitochondrial dysfunction as a core constitutive difference of the retinal pigment epithelium from patients with geographic atrophy. This study uncovers important differences in retinal pigment epithelium homeostasis associated with geographic atrophy.

FACS Enrichment of Total Interstitial Cells and Fibroblasts from Adult Mouse Ventricles.

Besides cardiomyocytes, the heart contains numerous interstitial cell types, including cardiac fibroblasts, endothelial cells, immune (myeloid and lymphoid) cells, and mural cells (pericytes and vascular smooth muscle cells), which play key roles in heart repair, regeneration, and disease. We recently published a comprehensive map of cardiac stromal cell heterogeneity and flux in healthy and infarcted hearts using single-cell RNA sequencing (scRNA-seq) ( Farbehi et al., 2019 ). Here, we describe the FACS (Fluorescent Activated Cell Sorting)-based method used in that study for isolation of two cardiac cell fractions from adult mouse ventricles: the total interstitial cell population (TIP; non-cardiomyocytes) and enriched (Pdgfra-GFP+) cardiac fibroblasts.

Role of PDGF-A/B Ligands in Cardiac Repair After Myocardial Infarction.

Platelet-derived growth factors (PDGFs) are powerful inducers of cellular mitosis, migration, angiogenesis, and matrix modulation that play pivotal roles in the development, homeostasis, and healing of cardiac tissues. PDGFs are key signaling molecules and important drug targets in the treatment of cardiovascular disease as multiple researchers have shown that delivery of recombinant PDGF ligands during or after myocardial infarction can reduce mortality and improve cardiac function in both rodents and porcine models. The mechanism involved cannot be easily elucidated due to the complexity of PDGF regulatory activities, crosstalk with other protein tyrosine kinase activators, and diversity of the pathological milieu. This review outlines the possible roles of PDGF ligands A and B in the healing of cardiac tissues including reduced cell death, improved vascularization, and improved extracellular matrix remodeling to improve cardiac architecture and function after acute myocardial injury. This review may highlight the use of recombinant PDGF-A and PDGF-B as a potential therapeutic modality in the treatment of cardiac injury.

Single-cell transcriptomics reveals involution mimicry during the specification of the basal breast cancer subtype.

Basal breast cancer is associated with younger age, early relapse, and a high mortality rate. Here, we use unbiased droplet-based single-cell RNA sequencing (RNA-seq) to elucidate the cellular basis of tumor progression during the specification of the basal breast cancer subtype from the luminal progenitor population in the MMTV-PyMT (mouse mammary tumor virus-polyoma middle tumor-antigen) mammary tumor model. We find that basal-like cancer cells resemble the alveolar lineage that is specified upon pregnancy and encompass the acquisition of an aberrant post-lactation developmental program of involution that triggers remodeling of the tumor microenvironment and metastatic dissemination. This involution mimicry is characterized by a highly interactive multicellular network, with involution cancer-associated fibroblasts playing a pivotal role in extracellular matrix remodeling and immunosuppression. Our results may partially explain the increased risk and poor prognosis of breast cancer associated with childbirth.

Droplet-based single cell RNAseq tools: a practical guide.

Droplet based scRNA-seq systems such as Drop-seq, inDrop and Chromium 10X have been the catalyst for the wide adoption of high-throughput scRNA-seq technologies in the research laboratory. In order to understand the capabilities of these systems to deeply interrogate biology; here we provide a practical guide through all the steps involved in a typical scRNA-seq experiment. Through comparing and contrasting these three main droplet based systems (and their derivatives), we provide an overview of all critical considerations in obtaining high quality and biologically relevant data. We also discuss the limitations of these systems and how they fit into the emerging field of Genomic Cytometry.

Single-cell expression profiling reveals dynamic flux of cardiac stromal, vascular and immune cells in health and injury.

Besides cardiomyocytes (CM), the heart contains numerous interstitial cell types which play key roles in heart repair, regeneration and disease, including fibroblast, vascular and immune cells. However, a comprehensive understanding of this interactive cell community is lacking. We performed single-cell RNA-sequencing of the total non-CM fraction and enriched (Pdgfra-GFP+) fibroblast lineage cells from murine hearts at days 3 and 7 post-sham or myocardial infarction (MI) surgery. Clustering of >30,000 single cells identified >30 populations representing nine cell lineages, including a previously undescribed fibroblast lineage trajectory present in both sham and MI hearts leading to a uniquely activated cell state defined in part by a strong anti-WNT transcriptome signature. We also uncovered novel myofibroblast subtypes expressing either pro-fibrotic or anti-fibrotic signatures. Our data highlight non-linear dynamics in myeloid and fibroblast lineages after cardiac injury, and provide an entry point for deeper analysis of cardiac homeostasis, inflammation, fibrosis, repair and regeneration.

A rapid co-culture stamping device for studying intercellular communication.

Regulation of tissue development and repair depends on communication between neighbouring cells. Recent advances in cell micro-contact printing and microfluidics have facilitated the in-vitro study of homotypic and heterotypic cell-cell interaction. Nonetheless, these techniques are still complicated to perform and as a result, are seldom used by biologists. We report here development of a temporarily sealed microfluidic stamping device which utilizes a novel valve design for patterning two adherent cell lines with well-defined interlacing configurations to study cell-cell interactions. We demonstrate post-stamping cell viability of >95%, the stamping of multiple adherent cell types, and the ability to control the seeded cell density. We also show viability, proliferation and migration of cultured cells, enabling analysis of co-culture boundary conditions on cell fate. We also developed an in-vitro model of endothelial and cardiac stem cell interactions, which are thought to regulate coronary repair after myocardial injury. The stamp is fabricated using microfabrication techniques, is operated with a lab pipettor and uses very low reagent volumes of 20 µl with cell injection efficiency of >70%. This easy-to-use device provides a general strategy for micro-patterning of multiple cell types and will be important for studying cell-cell interactions in a multitude of applications.