Helena Kilpinen

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.

Coordinated Effects of Sequence Variation on DNA Binding, Chromatin Structure, and Transcription

DNA Differences The extent to which genetic variation affects an individuals phenotype has been difficult to predict because the majority of variation lies outside the coding regions of genes. Now, three studies examine the extent to which genetic variation affects the chromatin of individuals with diverse ancestry and genetic variation (see the Perspective by Furey and Sethupathy). Kasowski et al. (p. 750, published online 17 October) examined how genetic variation affects differences in chromatin states and their correlation to histone modifications, as well as more general DNA binding factors. Kilpinen et al. (p. 744, published online 17 October) document how genetic variation is linked to allelic specificity in transcription factor binding, histone modifications, and transcription. McVicker et al. (p. 747, published online 17 October) identified how quantitative trait loci affect histone modifications in Yoruban individuals and established which specific transcription factors affect such modifications. Human genetic variation results in coordinated allelic variation across molecular phenotypes. [Also see Perspective by Furey and Sethupathy] 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.

Association of DISC1 with autism and Asperger syndrome

The DISC1 gene at 1q42 has generated considerable interest in various psychiatric diseases, since a balanced translocation interrupting the gene was found to cosegregate with schizophrenia and related mental illnesses in a large Scottish pedigree. To date, linkage and association findings to this locus have been replicated in several study samples ascertained for psychotic disorders. However, the biological function of DISC1 in neuronal development would suggest a potential role for this gene also in other, early onset neuropsychiatric disorders. Here we have addressed the allelic diversity of the DISC1, DISC2 and TRAX genes, clustered in 1q42, in Finnish families ascertained for infantile autism (97 families, naffected=138) and Asperger syndrome (29 families, naffected=143). We established association between autism and a DISC1 intragenic microsatellite (D1S2709; P=0.004). In addition, evidence for association to Asperger syndrome was observed with an intragenic single nucleotide polymorphism (SNP) of DISC1 (rs1322784; P=0.0058), as well as with a three-SNP haplotype (P=0.0013) overlapping the HEP3 haplotype, that was previously observed to associate with schizophrenia in Finnish families. The strongest associations were obtained with broad diagnostic categories for both disorders and with affected males only, in agreement with the previous sex-dependent effects reported for DISC1. These results would further support the involvement of DISC1 gene also in the etiopathogenesis of early onset neuropsychiatric disorders.

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.

How next-generation sequencing is transforming complex disease genetics

Progress in understanding the genetics of human disease is closely tied to technological developments in DNA sequencing. Recently, next-generation technology has transformed the scale of sequencing; compared to the methods used in the Human Genome Project, modern sequencers are 50000-fold faster. Complex disease genetics presents an immediate opportunity to use this technology to move from approaches using only partial information (linkage and genome-wide association studies, GWAS) to complete analysis of the relationship between genomic variation and phenotype. We first describe sequence-based improvements to existing study designs, followed by prioritization of both samples and genomic regions to be sequenced, and then address the ultimate goal of analyzing thousands of whole-genome sequences. Finally, we discuss how the same technology will also fundamentally change the way we understand the biological mechanisms underlying disease associations discovered through sequencing.

Linkage and linkage disequilibrium scan for autism loci in an extended pedigree from Finland

Population isolates, such as Finland, have proved beneficial in mapping rare causative genetic variants due to a limited number of founders resulting in reduced genetic heterogeneity and extensive linkage disequilibrium (LD). We have here used this special opportunity to identify rare alleles in autism by genealogically tracing 20 autism families into one extended pedigree with verified genealogical links reaching back to the 17th century. In this unique pedigree, we performed a dense microsatellite marker genome-wide scan of linkage and LD and followed initial findings with extensive fine-mapping. We identified a putative autism susceptibility locus at 19p13.3 and obtained further evidence for previously identified loci at 1q23 and 15q11-q13. Most promising candidate genes were TLE2 and TLE6 clustered at 19p13 and ATP1A2 at 1q23.

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.

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.This study identifies regulatory variants in sensory neurons derived from induced pluripotent stem cells. Despite differentiation-induced variability, an allele-specific method allowed detection of loci influencing gene expression, chromatin accessibility and RNA splicing.

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.

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