Common genetic variation drives molecular heterogeneity in human iPSCs.

Technology utilizing human induced pluripotent stem cells (iPS cells) has enormous potential to provide improved cellular models of human disease. However, variable genetic and phenotypic characterization of many existing iPS cell lines limits their potential use for research and therapy. Here we describe the systematic generation, genotyping and phenotyping of 711 iPS cell lines derived from 301 healthy individuals by the Human Induced Pluripotent Stem Cells Initiative. Our study outlines the major sources of genetic and phenotypic variation in iPS cells and establishes their suitability as models of complex human traits and cancer. Through genome-wide profiling we find that 5-46% of the variation in different iPS cell phenotypes, including differentiation capacity and cellular morphology, arises from differences between individuals. Additionally, we assess the phenotypic consequences of genomic copy-number alterations that are repeatedly observed in iPS cells. In addition, we present a comprehensive map of common regulatory variants affecting the transcriptome of human pluripotent cells.

Corrigendum: Common genetic variation drives molecular heterogeneity in human iPSCs.

This corrects the article DOI: 10.1038/nature22403.

A high-content platform to characterise human induced pluripotent stem cell lines.

Induced pluripotent stem cells (iPSCs) provide invaluable opportunities for future cell therapies as well as for studying human development, modelling diseases and discovering therapeutics. In order to realise the potential of iPSCs, it is crucial to comprehensively characterise cells generated from large cohorts of healthy and diseased individuals. The human iPSC initiative (HipSci) is assessing a large panel of cell lines to define cell phenotypes, dissect inter- and intra-line and donor variability and identify its key determinant components. Here we report the establishment of a high-content platform for phenotypic analysis of human iPSC lines. In the described assay, cells are dissociated and seeded as single cells onto 96-well plates coated with fibronectin at three different concentrations. This method allows assessment of cell number, proliferation, morphology and intercellular adhesion. Altogether, our strategy delivers robust quantification of phenotypic diversity within complex cell populations facilitating future identification of the genetic, biological and technical determinants of variance. Approaches such as the one described can be used to benchmark iPSCs from multiple donors and create novel platforms that can readily be tailored for disease modelling and drug discovery.

Investigation of 15 of the top candidate genes for late-onset Alzheimer's disease.

The 12 genome-wide association studies (GWAS) published to-date for late-onset Alzheimer's disease (LOAD) have identified over 40 candidate LOAD risk modifiers, in addition to apolipoprotein (APOE) ε4. A few of these novel LOAD candidate genes, namely BIN1, CLU, CR1, EXOC3L2 and PICALM, have shown consistent replication, and are thus credible LOAD susceptibility genes. To evaluate other promising LOAD candidate genes, we have added data from our large, case-control series (n=5,043) to meta-analyses of all published follow-up case-control association studies for six LOAD candidate genes that have shown significant association across multiple studies (TNK1, GAB2, LOC651924, GWA_14q32.13, PGBD1 and GALP) and for an additional nine previously suggested candidate genes. Meta-analyses remained significant at three loci after addition of our data: GAB2 (OR=0.78, p=0.007), LOC651924 (OR=0.91, p=0.01) and TNK1 (OR=0.92, p=0.02). Breslow-Day tests revealed significant heterogeneity between studies for GAB2 (p<0.0001) and GWA_14q32.13 (p=0.006). We have also provided suggestive evidence that PGBD1 (p=0.04) and EBF3 (p=0.03) are associated with age-at-onset of LOAD. Finally, we tested for interactions between these 15 genes, APOE ε4 and the five novel LOAD genes BIN1, CLU, CR1, EXOC3L2 and PICALM but none were significant after correction for multiple testing. Overall, this large, independent follow-up study for 15 of the top LOAD candidate genes provides support for GAB2 and LOC651924 (6q24.1) as risk modifiers of LOAD and novel associations between PGBD1 and EBF3 with age-at-onset.

Differential Expression of Insulin-Like Growth Factor 1 and Wnt Family Member 4 Correlates With Functional Heterogeneity of Human Dermal Fibroblasts.

Although human dermis contains distinct fibroblast subpopulations, the functional heterogeneity of fibroblast lines from different donors is under-appreciated. We identified one commercially sourced fibroblast line (c64a) that failed to express α-smooth muscle actin (α-SMA), a marker linked to fibroblast contractility, even when treated with transforming growth factor-ß1 (TGF-ß1). Gene expression profiling identified insulin-like growth factor 1 (IGF1) as being expressed more highly, and Asporin (ASPN) and Wnt family member 4 (WNT4) expressed at lower levels, in c64a fibroblasts compared to three fibroblast lines that had been generated in-house, independent of TGF-ß1 treatment. TGF-ß1 increased expression of C-X-C motif chemokine ligand 1 (CXCL1) in c64a cells to a greater extent than in the other lines. The c64a gene expression profile did not correspond to any dermal fibroblast subpopulation identified by single-cell RNAseq of freshly isolated human skin cells. In skin reconstitution assays, c64a fibroblasts did not support epidermal stratification as effectively as other lines tested. In fibroblast lines generated in-house, shRNA-mediated knockdown of IGF1 increased α-SMA expression without affecting epidermal stratification. Conversely, WNT4 knockdown had no consistent effect on α-SMA expression, but increased the ability of fibroblasts to support epidermal stratification. Thus, by comparing the properties of different lines of cultured dermal fibroblasts, we have identified IGF1 and WNT4 as candidate mediators of two distinct dermal functions: myofibroblast formation and epidermal maintenance.

Identifying Extrinsic versus Intrinsic Drivers of Variation in Cell Behavior in Human iPSC Lines from Healthy Donors.

Large cohorts of human induced pluripotent stem cells (iPSCs) from healthy donors are a potentially powerful tool for investigating the relationship between genetic variants and cellular behavior. Here, we integrate high content imaging of cell shape, proliferation, and other phenotypes with gene expression and DNA sequence datasets from over 100 human iPSC lines. By applying a dimensionality reduction approach, Probabilistic Estimation of Expression Residuals (PEER), we extracted factors that captured the effects of intrinsic (genetic concordance between different cell lines from the same donor) and extrinsic (cell responses to different fibronectin concentrations) conditions. We identify genes that correlate in expression with intrinsic and extrinsic PEER factors and associate outlier cell behavior with genes containing rare deleterious non-synonymous SNVs. Our study, thus, establishes a strategy for examining the genetic basis of inter-individual variability in cell behavior.

Dermal Blimp1 Acts Downstream of Epidermal TGFß and Wnt/ß-Catenin to Regulate Hair Follicle Formation and Growth.

B-lymphocyte-induced maturation protein 1 (Blimp1) is a transcriptional repressor that regulates cell growth and differentiation in multiple tissues, including skin. Although in the epidermis Blimp1 is important for keratinocyte and sebocyte differentiation, its role in dermal fibroblasts is unclear. Here we show that Blimp1 is dynamically regulated in dermal papilla cells during hair follicle (HF) morphogenesis and the postnatal hair cycle, preceding dermal Wnt/ß-catenin activation. Blimp1 ablation in E12.5 mouse dermal fibroblasts delayed HF morphogenesis and growth and prevented new HF formation after wounding. By combining targeted quantitative PCR screens with bioinformatic analysis and experimental validation we demonstrated that Blimp1 is both a target and a mediator of key dermal papilla inductive signaling pathways including transforming growth factor-ß and Wnt/ß-catenin. Epidermal overexpression of stabilized ß-catenin was able to override the HF defects in Blimp1 mutant mice, underlining the close reciprocal relationship between the dermal papilla and adjacent HF epithelial cells. Overall, our study reveals the functional role of Blimp1 in promoting the dermal papilla inductive signaling cascade that initiates HF growth.

Erratum to: Genetically-controlled Vesicle-Associated Membrane Protein 1 expression may contribute to Alzheimer's pathophysiology and susceptibility.

Following publication of this work, we noticed that we inadvertently failed to include Dr Ferenc Deák in the author list. The author list has now been corrected and the amended authors' contributions section has been modified accordingly below.

Genetically-controlled Vesicle-Associated Membrane Protein 1 expression may contribute to Alzheimer's pathophysiology and susceptibility.

BACKGROUND: Alzheimer's disease is a neurodegenerative disorder in which extracellular deposition of ß-amyloid (Aß) oligomers causes synaptic injury resulting in early memory loss, altered homeostasis, accumulation of hyperphosphorylated tau and cell death. Since proteins in the SNAP (Soluble N-ethylmaleimide-sensitive factor Attachment Protein) REceptors (SNARE) complex are essential for neuronal Aß release at pre-synaptic terminals, we hypothesized that genetically controlled SNARE expression could alter neuronal Aß release at the synapse and hence play an early role in Alzheimer's pathophysiology. RESULTS: Here we report 5 polymorphisms in Vesicle-Associated Membrane Protein 1 (VAMP1), a gene encoding a member of the SNARE complex, associated with bidirectionally altered cerebellar VAMP1 transcript levels (all p<0.05). At the functional level, we demonstrated that control of VAMP1 expression by heterogeneous knockdown in mice resulted in up to 74% reduction in neuronal Aß exocytosis (p<0.001). We performed a case-control association study of the 5 VAMP1 expression regulating polymorphisms in 4,667 Alzheimer's disease patients and 6,175 controls to determine their contribution to Alzheimer's disease risk. We found that polymorphisms associated with increased brain VAMP1 transcript levels conferred higher risk for Alzheimer's disease than those associated with lower VAMP1 transcript levels (p=0.03). Moreover, we also report a modest protective association for a common VAMP1 polymorphism with Alzheimer's disease risk (OR=0.88, p=0.03). This polymorphism was associated with decreased VAMP1 transcript levels (p=0.02) and was functionally active in a dual luciferase reporter gene assay (p<0.01). CONCLUSIONS: Genetically regulated VAMP1 expression in the brain may modify both Alzheimer's disease risk and may contribute to Alzheimer's pathophysiology.

Human pluripotent stem cells on artificial microenvironments: a high content perspective.

Self-renewing stem cell populations are increasingly considered as resources for cell therapy and tools for drug discovery. Human pluripotent stem (hPS) cells in particular offer a virtually unlimited reservoir of homogeneous cells and can be differentiated toward diverse lineages. Many diseases show impairment in self-renewal or differentiation, abnormal lineage choice or other aberrant cell behavior in response to chemical or physical cues. To investigate these responses, there is a growing interest in the development of specific assays using hPS cells, artificial microenvironments and high content analysis. Several hurdles need to be overcome that can be grouped into three areas: (i) availability of robust, homogeneous, and consistent cell populations as a starting point; (ii) appropriate understanding and use of chemical and physical microenvironments; (iii) development of assays that dissect the complexity of cell populations in tissues while mirroring specific aspects of their behavior. Here we review recent progress in the culture of hPS cells and we detail the importance of the environment surrounding the cells with a focus on synthetic material and suitable high content analysis approaches. The technologies described, if properly combined, have the potential to create a paradigm shift in the way diseases are modeled and drug discovery is performed.

Linking protective GAB2 variants, increased cortical GAB2 expression and decreased Alzheimer's disease pathology.

GRB-associated binding protein 2 (GAB2) represents a compelling genome-wide association signal for late-onset Alzheimer's disease (LOAD) with reported odds ratios (ORs) ranging from 0.75-0.85. We tested eight GAB2 variants in four North American Caucasian case-control series (2,316 LOAD, 2,538 controls) for association with LOAD. Meta-analyses revealed ORs ranging from (0.61-1.20) with no significant association (all p>0.32). Four variants were hetergeneous across the populations (all p<0.02) due to a potentially inflated effect size (OR = 0.61-0.66) only observed in the smallest series (702 LOAD, 209 controls). Despite the lack of association in our series, the previously reported protective association for GAB2 remained after meta-analyses of our data with all available previously published series (11,952-22,253 samples; OR = 0.82-0.88; all p<0.04). Using a freely available database of lymphoblastoid cell lines we found that protective GAB2 variants were associated with increased GAB2 expression (p = 9.5×10(-7)-9.3×10(-6)). We next measured GAB2 mRNA levels in 249 brains and found that decreased neurofibrillary tangle (r = -0.34, p = 0.0006) and senile plaque counts (r = -0.32, p = 0.001) were both good predictors of increased GAB2 mRNA levels albeit that sex (r = -0.28, p = 0.005) may have been a contributing factor. In summary, we hypothesise that GAB2 variants that are protective against LOAD in some populations may act functionally to increase GAB2 mRNA levels (in lymphoblastoid cells) and that increased GAB2 mRNA levels are associated with significantly decreased LOAD pathology. These findings support the hypothesis that Gab2 may protect neurons against LOAD but due to significant population heterogeneity, it is still unclear whether this protection is detectable at the genetic level.