Search

Despite their widespread use in research, there has not yet been a systematic genomic analysis of human embryonic stem cell (hESC) lines at a single-nucleotide resolution. We therefore performed whole-genome sequencing (WGS) of 143 hESC lines and annotated their single-nucleotide and structural genetic variants. We found that while a substantial fraction of hESC lines contained large deleterious structural variants, finer-scale structural and single-nucleotide variants (SNVs) that are ascertainable only through WGS analyses were present in hESC genomes and human blood-derived genomes at similar frequencies. Moreover, WGS allowed us to identify SNVs associated with cancer and other diseases that could alter cellular phenotypes and compromise the safety of hESC-derived cellular products transplanted into humans. As a resource to enable reproducible hESC research and safer translation, we provide a user-friendly WGS data portal and a data-driven scheme for cell line maintenance and selection.

SUMMARYVariation in the human genome contributes to abundant diversity in human traits and vulnerabilities, but the underlying molecular and cellular mechanisms are not yet known, and will need scalable approaches to accelerate their recognition. Here, we advanced and applied an experimental platform that analyzes genetic, molecular, and phenotypic heterogeneity across cells from very many human donors cultured in a single, shared in vitro environment, with algorithms (Dropulation and Census-seq) for assigning phenotypes to individual donors. We used natural genetic variation and synthetic (CRISPR-Cas9) genetic perturbations to analyze the vulnerability of neural progenitor cells to infection with Zika virus. These analyses identified a common variant in the antiviral IFITM3 gene that regulated IFITM3 expression and explained most inter-individual variation in NPCs’ susceptibility to Zika virus infectivity. These and other approaches could provide scalable ways to recognize the impact of genes and genetic variation on cellular phenotypes.HIGHLIGHTSMeasuring cellular phenotypes in iPSCs and hPSC-derived NPCs from many donorsEffects of donor sex, cell source, genetic and other variables on hPSC RNA expressionNatural genetic variation and synthetic perturbation screens both identify IFITM3 in NPC susceptibility to Zika virusA common genetic variant in IFITM3 explains most inter-individual variation in NPC susceptibility to Zika virus

SummaryTens of thousands of genetic variants shape human phenotypes, mostly by unknown cellular mechanisms. Here we describe Census-seq, a way to measure cellular phenotypes in cells from many people simultaneously. Analogous to pooled CRISPR screens but for natural variation, Census-seq associates cellular phenotypes to donors’ genotypes by quantifying the presence of each donor’s DNA in cell “villages” before and after sorting or selection for cellular traits of interest. Census-seq enables population-scale cell-biological phenotyping with low cost and high internal control. We demonstrate Census-seq through investigation of genetic effects on the SMN protein whose deficiency underlies spinal muscular atrophy (SMA). Census-seq quantified and mapped effects of many common alleles on SMN protein levels and response to SMN-targeted therapeutics, including a common, cryptic non-responder allele. We provide tools enabling population-scale cell experiments and explain how Census-seq can be used to map genetic effects on diverse cell phenotypes.Abstract FigureHighlightsCensus-seq reveals how inherited genetic variation affects cell phenotypesGenetic analysis of cellular traits in cell villages of >100 donorsCharacterizing human alleles that shape SMN protein expression and drug responsesDevelopment of protocols and software to enable cellular population genetics

The identification of convergent phenotypes in different models of psychiatric illness highlights robust phenotypes that are more likely to be implicated in disease pathophysiology. Here, we utilize human iPSCs harboring distinct mutations in DISC1 that have been found in families with major mental illness. One mutation was engineered to mimic the consequences on DISC1 protein of a balanced translocation linked to mental illness in a Scottish pedigree; the other mutation was identified in an American pedigree with a high incidence of mental illness. Directed differentiation of these iPSCs using NGN2 expression shows rapid conversion to a homogenous population of mature excitatory neurons. Both DISC1 mutations result in reduced DISC1 protein expression, and show subtle effects on certain presynaptic proteins. In addition, RNA sequencing and qPCR showed decreased expression of UNC5D, DPP10, PCDHA6, and ZNF506 in neurons with both DISC1 mutations. Longitudinal analysis of neurite outgrowth revealed decreased neurite outgrowth in neurons with each DISC1 mutation, which was mimicked by UNC5D knockdown and rescued by transient upregulation of endogenous UNC5D. This study shows a narrow range of convergent phenotypes of two mutations found in families with major mental illness, and implicates dysregulated netrin signaling in DISC1 biology.

While a mutation in C9ORF72 is the most common genetic contributor to amyotrophic lateral sclerosis (ALS), much remains to be learned concerning the function of the protein normally encoded at this locus. To elaborate further on functions for C9ORF72, we used quantitative mass spectrometry-based proteomics to identify interacting proteins in motor neurons and found that its long isoform complexes with and stabilizes SMCR8, which further enables interaction with WDR41. To study the organismal and cellular functions for this tripartite complex, we generated Smcr8 loss-of-function mutant mice and found that they developed phenotypes also observed in C9orf72 loss-of-function animals, including autoimmunity. Along with a loss of tolerance for many nervous system autoantigens, we found increased lysosomal exocytosis in Smcr8 mutant macrophages. In addition to elevated surface Lamp1 (lysosome-associated membrane protein 1) expression, we also observed enhanced secretion of lysosomal components—phenotypes that we subsequently observed in C9orf72 loss-of-function macrophages. Overall, our findings demonstrate that C9ORF72 and SMCR8 have interdependent functions in suppressing autoimmunity as well as negatively regulating lysosomal exocytosis—processes of potential importance to ALS.

Human pluripotent stem cells (hPS cells) can self-renew indefinitely, making them an attractive source for regenerative therapies. This expansion potential has been linked with the acquisition of large copy number variants that provide mutated cells with a growth advantage in culture. The nature, extent and functional effects of other acquired genome sequence mutations in cultured hPS cells are not known. Here we sequence the protein-coding genes (exomes) of 140 independent human embryonic stem cell (hES cell) lines, including 26 lines prepared for potential clinical use. We then apply computational strategies for identifying mutations present in a subset of cells in each hES cell line. Although such mosaic mutations were generally rare, we identified five unrelated hES cell lines that carried six mutations in the TP53 gene that encodes the tumour suppressor P53. The TP53 mutations we observed are dominant negative and are the mutations most commonly seen in human cancers. We found that the TP53 mutant allelic fraction increased with passage number under standard culture conditions, suggesting that the P53 mutations confer selective advantage. We then mined published RNA sequencing data from 117 hPS cell lines, and observed another nine TP53 mutations, all resulting in coding changes in the DNA-binding domain of P53. In three lines, the allelic fraction exceeded 50%, suggesting additional selective advantage resulting from the loss of heterozygosity at the TP53 locus. As the acquisition and expansion of cancer-associated mutations in hPS cells may go unnoticed during most applications, we suggest that careful genetic characterization of hPS cells and their differentiated derivatives be carried out before clinical use.

To discover novel genes underlying amyotrophic lateral sclerosis (ALS), we aggregated exomes from 3,864 cases and 7,839 ancestry-matched controls. We observed a significant excess of rare protein-truncating variants among ALS cases, and these variants were concentrated in constrained genes. Through gene level analyses, we replicated known ALS genes including SOD1, NEK1 and FUS. We also observed multiple distinct protein-truncating variants in a highly constrained gene, DNAJC7. The signal in DNAJC7 exceeded genome-wide significance, and immunoblotting assays showed depletion of DNAJC7 protein in fibroblasts in a patient with ALS carrying the p.Arg156Ter variant. DNAJC7 encodes a member of the heat-shock protein family, HSP40, which, along with HSP70 proteins, facilitates protein homeostasis, including folding of newly synthesized polypeptides and clearance of degraded proteins. When these processes are not regulated, misfolding and accumulation of aberrant proteins can occur and lead to protein aggregation, which is a pathological hallmark of neurodegeneration. Our results highlight DNAJC7 as a novel gene for ALS.

SUMMARYTo discover novel genetic risk factors underlying amyotrophic lateral sclerosis (ALS), we aggregated exomes from 3,864 cases and 7,839 ancestry matched controls. We observed a significant excess of ultra-rare and rare protein-truncating variants (PTV) among ALS cases, which was primarily concentrated in constrained genes; however, a significant enrichment in PTVs does persist in the remaining exome. Through gene level analyses, known ALS genes, SOD1, NEK1, and FUS, were the most strongly associated with disease status. We also observed suggestive statistical evidence for multiple novel genes including DNAJC7, which is a highly constrained gene and a member of the heat shock protein family (HSP40). HSP40 proteins, along with HSP70 proteins, facilitate protein homeostasis, such as folding of newly synthesized polypeptides, and clearance of degraded proteins. When these processes are not regulated, misfolding and accumulation of degraded proteins can occur leading to aberrant protein aggregation, one of the pathological hallmarks of neurodegeneration.

SummaryGenomic DNA replicates in a choreographed temporal order that impacts the distribution of mutations along the genome. We show here that DNA replication timing is shaped by genetic polymorphisms that act in cis upon megabase-scale DNA segments. In genome sequences from proliferating cells, read depth along chromosomes reflected DNA replication activity in those cells. We used this relationship to analyze variation in replication timing among 161 individuals sequenced by the 1000 Genomes Project. Genome-wide association of replication timing with genetic variation identified 16 loci at which inherited alleles associate with replication timing. We call these “replication timing quantitative trait loci” (rtQTLs). rtQTLs involved the differential use of replication origins, exhibited allele-specific effects on replication timing, and associated with gene expression variation at megabase scales. Our results show replication timing to be shaped by genetic polymorphism and identify a means by which inherited polymorphism regulates the mutability of nearby sequences.

Genetic drift is a mechanism of population divergence that is important in the evolution of plants and animals but is thought to be rare in free-living microorganisms because of their typically large population sizes and unrestricted means of dispersal. We used both phylogenetic and insertion sequence (IS) element analyses in hyperthermophilic archaea of the genus Pyrococcus to test the hypothesis that genetic drift played an important role in the diversification of these microorganisms. Multilocus sequence typing of a collection of 36 isolates of Pyrococcus, from different hydrothermal systems in the Pacific Ocean and the Mediterranean Sea, revealed that Pyrococcus populations from different geographic locations are genetically differentiated. Analysis of IS elements in these isolates exposed their presence in all individuals of only one geographically isolated lineage, that of Vulcano Island in the Mediterranean Sea. Detailed sequence analysis of six selected IS elements in the Vulcano population showed that these elements cause deleterious genomic alterations, including inactivation of gene function. The high frequency of IS elements in the sampled population together with their observed harmful effects in the genome of Pyrococcus provide molecular evidence that the Vulcano Island population of Pyrococcus is geographically isolated and that those genetic mobile elements have been brought up to high frequency by genetic drift. Thus, genetic drift resulting from physical isolation should be considered as a factor influencing differentiation in prokaryotes.