Somatic mutations alter the differentiation outcomes of iPSC-derived neurons.

The use of induced pluripotent stem cells (iPSC) as models for development and human disease has enabled the study of otherwise inaccessible tissues. A remaining challenge in developing reliable models is our limited understanding of the factors driving irregular differentiation of iPSCs, particularly the impact of acquired somatic mutations. We leveraged data from a pooled dopaminergic neuron differentiation experiment of 238 iPSC lines profiled with single-cell RNA and whole-exome sequencing to study how somatic mutations affect differentiation outcomes. We found that deleterious somatic mutations in key developmental genes, notably the BCOR gene, are strongly associated with failure in dopaminergic neuron differentiation and a larger proliferation rate in culture. We further identified broad differences in cell type composition between incorrectly and successfully differentiating lines, as well as significant changes in gene expression contributing to the inhibition of neurogenesis. Our work calls for caution in interpreting differentiation-related phenotypes in disease-modeling experiments.

Histone lysine methyltransferase-related neurodevelopmental disorders: current knowledge and saRNA future therapies.

Neurodevelopmental disorders encompass a group of debilitating diseases presenting with motor and cognitive dysfunction, with variable age of onset and disease severity. Advances in genetic diagnostic tools have facilitated the identification of several monogenic chromatin remodeling diseases that cause Neurodevelopmental disorders. Chromatin remodelers play a key role in the neuro-epigenetic landscape and regulation of brain development; it is therefore not surprising that mutations, leading to loss of protein function, result in aberrant neurodevelopment. Heterozygous, usually de novo mutations in histone lysine methyltransferases have been described in patients leading to haploinsufficiency, dysregulated protein levels and impaired protein function. Studies in animal models and patient-derived cell lines, have highlighted the role of histone lysine methyltransferases in the regulation of cell self-renewal, cell fate specification and apoptosis. To date, in depth studies of histone lysine methyltransferases in oncology have provided strong evidence of histone lysine methyltransferase dysregulation as a determinant of cancer progression and drug resistance. As a result, histone lysine methyltransferases have become an important therapeutic target for the treatment of different cancer forms. Despite recent advances, we still lack knowledge about the role of histone lysine methyltransferases in neuronal development. This has hampered both the study and development of precision therapies for histone lysine methyltransferases-related Neurodevelopmental disorders. In this review, we will discuss the current knowledge of the role of histone lysine methyltransferases in neuronal development and disease progression. We will also discuss how RNA-based technologies using small-activating RNAs could potentially provide a novel therapeutic approach for the future treatment of histone lysine methyltransferase haploinsufficiency in these Neurodevelopmental disorders, and how they could be first tested in state-of-the-art patient-derived neuronal models.

Highly accurate quantification of allelic gene expression for population and disease genetics.

Analysis of allele-specific gene expression (ASE) is a powerful approach for studying gene regulation, particularly when sample sizes are small, such as for rare diseases, or when studying the effects of rare genetic variation. However, detection of ASE events relies on accurate alignment of RNA sequencing reads, where challenges still remain, particularly for reads containing genetic variants or those that align to many different genomic locations. We have developed the Personalised ASE Caller (PAC), a tool that combines multiple steps to improve the quantification of allelic reads, including personalized (i.e., diploid) read alignment with improved allocation of multimapping reads. Using simulated RNA sequencing data, we show that PAC outperforms standard alignment approaches for ASE detection, reducing the number of sites with incorrect biases (>10%) by ∼80% and increasing the number of sites that can be reliably quantified by ∼3%. Applying PAC to real RNA sequencing data from 670 whole-blood samples, we show that genetic regulatory signatures inferred from ASE data more closely match those from population-based methods that are less prone to alignment biases. Finally, we use PAC to characterize cell type-specific ASE events that would be missed by standard alignment approaches, and in doing so identify disease relevant genes that may modulate their effects through the regulation of gene expression. PAC can be applied to the vast quantity of existing RNA sequencing data sets to better understand a wide array of fundamental biological and disease processes.

CRISPR activation enables high-fidelity reprogramming into human pluripotent stem cells.

Conventional reprogramming methods rely on the ectopic expression of transcription factors to reprogram somatic cells into induced pluripotent stem cells (iPSCs). The forced expression of transcription factors may lead to off-target gene activation and heterogeneous reprogramming, resulting in the emergence of alternative cell types and aberrant iPSCs. Activation of endogenous pluripotency factors by CRISPR activation (CRISPRa) can reduce this heterogeneity. Here, we describe a high-efficiency reprogramming of human somatic cells into iPSCs using optimized CRISPRa. Efficient reprogramming was dependent on the additional targeting of the embryo genome activation-enriched Alu-motif and the miR-302/367 locus. Single-cell transcriptome analysis revealed that the optimized CRISPRa reprogrammed cells more directly and specifically into the pluripotent state when compared to the conventional reprogramming method. These findings support the use of CRISPRa for high-quality pluripotent reprogramming of human cells.

Differentiation of human induced pluripotent stem cells into cortical neural stem cells.

Efficient and effective methods for converting human induced pluripotent stem cells into differentiated derivatives are critical for performing robust, large-scale studies of development and disease modelling, and for providing a source of cells for regenerative medicine. Here, we describe a 14-day neural differentiation protocol which allows for the scalable, simultaneous differentiation of multiple iPSC lines into cortical neural stem cells We currently employ this protocol to differentiate and compare sets of engineered iPSC lines carrying loss of function alleles in developmental disorder associated genes, alongside isogenic wildtype controls. Using RNA sequencing (RNA-Seq), we can examine the changes in gene expression brought about by each disease gene knockout, to determine its impact on neural development and explore mechanisms of disease. The 10-day Neural Induction period uses the well established dual-SMAD inhibition approach combined with Wnt/ß-Catenin inhibition to selectively induce formation of cortical NSCs. This is followed by a 4-day Neural Maintenance period facilitating NSC expansion and rosette formation, and NSC cryopreservation. We also describe methods for thawing and passaging the cryopreserved NSCs, which are useful in confirming their viability for further culture. Routine implementation of immunocytochemistry Quality Control confirms the presence of PAX6-positive and/or FOXG1-positive NSCs and the absence of OCT4-positive iPSCs after differentiation. RNA-Seq, flow cytometry, immunocytochemistry (ICC) and RT-qPCR provide additional confirmation of robust presence of NSC markers in the differentiated cells. The broader utility and application of our protocol is demonstrated by the successful differentiation of wildtype iPSC lines from five additional independent donors. This paper thereby describes an efficient method for the production of large numbers of high purity cortical NSCs, which are widely applicable for downstream research into developmental mechanisms, further differentiation into postmitotic cortical neurons, or other applications such as large-scale drug screening experiments.

Identification of rare and common regulatory variants in pluripotent cells using population-scale transcriptomics.

Induced pluripotent stem cells (iPSCs) are an established cellular system to study the impact of genetic variants in derived cell types and developmental contexts. However, in their pluripotent state, the disease impact of genetic variants is less well known. Here, we integrate data from 1,367 human iPSC lines to comprehensively map common and rare regulatory variants in human pluripotent cells. Using this population-scale resource, we report hundreds of new colocalization events for human traits specific to iPSCs, and find increased power to identify rare regulatory variants compared with somatic tissues. Finally, we demonstrate how iPSCs enable the identification of causal genes for rare diseases.

Population-scale proteome variation in human induced pluripotent stem cells.

Human disease phenotypes are driven primarily by alterations in protein expression and/or function. To date, relatively little is known about the variability of the human proteome in populations and how this relates to variability in mRNA expression and to disease loci. Here, we present the first comprehensive proteomic analysis of human induced pluripotent stem cells (iPSC), a key cell type for disease modelling, analysing 202 iPSC lines derived from 151 donors, with integrated transcriptome and genomic sequence data from the same lines. We characterised the major genetic and non-genetic determinants of proteome variation across iPSC lines and assessed key regulatory mechanisms affecting variation in protein abundance. We identified 654 protein quantitative trait loci (pQTLs) in iPSCs, including disease-linked variants in protein-coding sequences and variants with trans regulatory effects. These include pQTL linked to GWAS variants that cannot be detected at the mRNA level, highlighting the utility of dissecting pQTL at peptide level resolution.

Genomic basis for RNA alterations in cancer.

Transcript alterations often result from somatic changes in cancer genomes1. Various forms of RNA alterations have been described in cancer, including overexpression2, altered splicing3 and gene fusions4; however, it is difficult to attribute these to underlying genomic changes owing to heterogeneity among patients and tumour types, and the relatively small cohorts of patients for whom samples have been analysed by both transcriptome and whole-genome sequencing. Here we present, to our knowledge, the most comprehensive catalogue of cancer-associated gene alterations to date, obtained by characterizing tumour transcriptomes from 1,188 donors of the Pan-Cancer Analysis of Whole Genomes (PCAWG) Consortium of the International Cancer Genome Consortium (ICGC) and The Cancer Genome Atlas (TCGA)5. Using matched whole-genome sequencing data, we associated several categories of RNA alterations with germline and somatic DNA alterations, and identified probable genetic mechanisms. Somatic copy-number alterations were the major drivers of variations in total gene and allele-specific expression. We identified 649 associations of somatic single-nucleotide variants with gene expression in cis, of which 68.4% involved associations with flanking non-coding regions of the gene. We found 1,900 splicing alterations associated with somatic mutations, including the formation of exons within introns in proximity to Alu elements. In addition, 82% of gene fusions were associated with structural variants, including 75 of a new class, termed 'bridged' fusions, in which a third genomic location bridges two genes. We observed transcriptomic alteration signatures that differ between cancer types and have associations with variations in DNA mutational signatures. This compendium of RNA alterations in the genomic context provides a rich resource for identifying genes and mechanisms that are functionally implicated in cancer.

Publisher Correction: Molecular and functional variation in iPSC-derived sensory neurons.

An amendment to this paper has been published and can be accessed via a link at the top of the paper.

Molecular and functional variation in iPSC-derived sensory neurons.

Induced pluripotent stem cells (iPSCs), and cells derived from them, have become key tools for modeling biological processes, particularly in cell types that are difficult to obtain from living donors. Here we present a map of regulatory variants in iPSC-derived neurons, based on 123 differentiations of iPSCs to a sensory neuronal fate. Gene expression was more variable across cultures than in primary dorsal root ganglion, particularly for genes related to nervous system development. Using single-cell RNA-sequencing, we found that the number of neuronal versus contaminating cells was influenced by iPSC culture conditions before differentiation. Despite high differentiation-induced variability, our allele-specific method detected thousands of quantitative trait loci (QTLs) that influenced gene expression, chromatin accessibility, and RNA splicing. On the basis of these detected QTLs, we estimate that recall-by-genotype studies that use iPSC-derived cells will require cells from at least 20-80 individuals to detect the effects of regulatory variants with moderately large effect sizes.

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

This corrects the article DOI: 10.1038/nature22403.

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.

Mondo-Mlx Mediates Organismal Sugar Sensing through the Gli-Similar Transcription Factor Sugarbabe.

The ChREBP/Mondo-Mlx transcription factors are activated by sugars and are essential for sugar tolerance. They promote the conversion of sugars to lipids, but beyond this, their physiological roles are insufficiently understood. Here, we demonstrate that in an organism-wide setting in Drosophila, Mondo-Mlx controls the majority of sugar-regulated genes involved in nutrient digestion and transport as well as carbohydrate, amino acid, and lipid metabolism. Furthermore, human orthologs of the Mondo-Mlx targets display enrichment among gene variants associated with high circulating triglycerides. In addition to direct regulation of metabolic genes, Mondo-Mlx maintains metabolic homeostasis through downstream effectors, including the Activin ligand Dawdle and the Gli-similar transcription factor Sugarbabe. Sugarbabe controls a subset of Mondo-Mlx-dependent processes, including de novo lipogenesis and fatty acid desaturation. In sum, Mondo-Mlx is a master regulator of other sugar-responsive pathways essential for adaptation to a high-sugar diet.

Population Variation and Genetic Control of Modular Chromatin Architecture in Humans.

Chromatin state variation at gene regulatory elements is abundant across individuals, yet we understand little about the genetic basis of this variability. Here, we profiled several histone modifications, the transcription factor (TF) PU.1, RNA polymerase II, and gene expression in lymphoblastoid cell lines from 47 whole-genome sequenced individuals. We observed that distinct cis-regulatory elements exhibit coordinated chromatin variation across individuals in the form of variable chromatin modules (VCMs) at sub-Mb scale. VCMs were associated with thousands of genes and preferentially cluster within chromosomal contact domains. We mapped strong proximal and weak, yet more ubiquitous, distal-acting chromatin quantitative trait loci (cQTL) that frequently explain this variation. cQTLs were associated with molecular activity at clusters of cis-regulatory elements and mapped preferentially within TF-bound regions. We propose that local, sequence-independent chromatin variation emerges as a result of genetic perturbations in cooperative interactions between cis-regulatory elements that are located within the same genomic domain.

Synaptic, transcriptional and chromatin genes disrupted in autism.

The genetic architecture of autism spectrum disorder involves the interplay of common and rare variants and their impact on hundreds of genes. Using exome sequencing, here we show that analysis of rare coding variation in 3,871 autism cases and 9,937 ancestry-matched or parental controls implicates 22 autosomal genes at a false discovery rate (FDR) < 0.05, plus a set of 107 autosomal genes strongly enriched for those likely to affect risk (FDR < 0.30). These 107 genes, which show unusual evolutionary constraint against mutations, incur de novo loss-of-function mutations in over 5% of autistic subjects. Many of the genes implicated encode proteins for synaptic formation, transcriptional regulation and chromatin-remodelling pathways. These include voltage-gated ion channels regulating the propagation of action potentials, pacemaking and excitability-transcription coupling, as well as histone-modifying enzymes and chromatin remodellers-most prominently those that mediate post-translational lysine methylation/demethylation modifications of histones.

Identification and removal of low-complexity sites in allele-specific analysis of ChIP-seq data.

MOTIVATION: High-throughput sequencing technologies enable the genome-wide analysis of the impact of genetic variation on molecular phenotypes at unprecedented resolution. However, although powerful, these technologies can also introduce unexpected artifacts. RESULTS: We investigated the impact of library amplification bias on the identification of allele-specific (AS) molecular events from high-throughput sequencing data derived from chromatin immunoprecipitation assays (ChIP-seq). Putative AS DNA binding activity for RNA polymerase II was determined using ChIP-seq data derived from lymphoblastoid cell lines of two parent-daughter trios. We found that, at high-sequencing depth, many significant AS binding sites suffered from an amplification bias, as evidenced by a larger number of clonal reads representing one of the two alleles. To alleviate this bias, we devised an amplification bias detection strategy, which filters out sites with low read complexity and sites featuring a significant excess of clonal reads. This method will be useful for AS analyses involving ChIP-seq and other functional sequencing assays. AVAILABILITY: The R package abs filter for library clonality simulations and detection of amplification-biased sites is available from http://updepla1srv1.epfl.ch/waszaks/absfilter

Allele-specific regulation of DISC1 expression by miR-135b-5p.

Disrupted-in-schizophrenia-1 (DISC1) gene has been established as a risk factor for various neuropsychiatric phenotypes. Both coding and regulatory variants in DISC1 have been identified and associated with these phenotypes in genetic studies. MicroRNAs (miRNAs) are important regulators of protein coding genes. Since the miRNA-mRNA target recognition mechanism is vulnerable to disruption by DNA polymorphisms, we investigated whether polymorphisms in the DISC1 3'UTR affect binding of miRNAs and lead to allele-specific regulation of DISC1. We identified four predicted polymorphic miRNA target sites in the DISC1 3'UTR, and demonstrated that miR-135b-5p regulates the level of DISC1 mRNA. Moreover, DISC1 regulation by miR-135b-5p is allele specific: miR-135b-5p only binds to the major allele (A) of rs11122396, not to the minor allele (G). Thus, the G allele may be functionally related to the DISC1-associated phenotypes by abolishing regulation by miR-135b-5p, leading to elevated DISC1 levels.

Coordinated effects of sequence variation on DNA binding, chromatin structure, and transcription.

DNA sequence variation has been associated with quantitative changes in molecular phenotypes such as gene expression, but its impact on chromatin states is poorly characterized. To understand the interplay between chromatin and genetic control of gene regulation, we quantified allelic variability in transcription factor binding, histone modifications, and gene expression within humans. We found abundant allelic specificity in chromatin and extensive local, short-range, and long-range allelic coordination among the studied molecular phenotypes. We observed genetic influence on most of these phenotypes, with histone modifications exhibiting strong context-dependent behavior. Our results implicate transcription factors as primary mediators of sequence-specific regulation of gene expression programs, with histone modifications frequently reflecting the primary regulatory event.

Transcriptome and genome sequencing uncovers functional variation in humans.

Genome sequencing projects are discovering millions of genetic variants in humans, and interpretation of their functional effects is essential for understanding the genetic basis of variation in human traits. Here we report sequencing and deep analysis of messenger RNA and microRNA from lymphoblastoid cell lines of 462 individuals from the 1000 Genomes Project--the first uniformly processed high-throughput RNA-sequencing data from multiple human populations with high-quality genome sequences. We discover extremely widespread genetic variation affecting the regulation of most genes, with transcript structure and expression level variation being equally common but genetically largely independent. Our characterization of causal regulatory variation sheds light on the cellular mechanisms of regulatory and loss-of-function variation, and allows us to infer putative causal variants for dozens of disease-associated loci. Altogether, this study provides a deep understanding of the cellular mechanisms of transcriptome variation and of the landscape of functional variants in the human genome.