Analyzing the functional effects of DNA variants with gene editing

Continual advancements in genomics have led to an ever-widening disparity between the rate of discovery of genetic variants and our current understanding of their functions and potential roles in disease. Systematic methods for phenotyping DNA variants are required to effectively translate genomics data into improved outcomes for patients with genetic diseases. To make the biggest impact, these approaches must be scalable and accurate, faithfully reflect disease biology, and define complex disease mechanisms. We compare current methods to analyze the function of variants in their endogenous DNA context using genome editing strategies, such as saturation genome editing, base editing and prime editing. We discuss how these technologies can be linked to high-content readouts to gain deep mechanistic insights into variant effects. Finally, we highlight key challenges that need to be addressed to bridge the genotype to phenotype gap, and ultimately improve the diagnosis and treatment of genetic diseases.

Dynein and dynactin move long-range but are delivered separately to the axon tip

Axonal transport is essential for neuronal survival. This is driven by microtubule motors including dynein, which transports cargo from the axon tip back to the cell body. This function requires its cofactor dynactin and regulators LIS1 and NDEL1. Due to difficulties imaging dynein at a single-molecule level, it is unclear how this motor and its regulators coordinate transport along the length of the axon. Here, we use a neuron-inducible human stem cell line (NGN2-OPTi-OX) to endogenously tag dynein components and visualize them at a near-single molecule regime. In the retrograde direction, we find that dynein and dynactin can move the entire length of the axon (>500 µm). Furthermore, LIS1 and NDEL1 also undergo long-distance movement, despite being mainly implicated with the initiation of dynein transport. Intriguingly, in the anterograde direction, dynein/LIS1 moves faster than dynactin/NDEL1, consistent with transport on different cargos. Therefore, neurons ensure efficient transport by holding dynein/dynactin on cargos over long distances but keeping them separate until required.

scSNV-seq: high-throughput phenotyping of single nucleotide variants by coupled single-cell genotyping and transcriptomics

CRISPR screens with single-cell transcriptomic readouts are a valuable tool to understand the effect of genetic perturbations including single nucleotide variants (SNVs) associated with diseases. Interpretation of these data is currently limited as genotypes cannot be accurately inferred from guide RNA identity alone. scSNV-seq overcomes this limitation by coupling single-cell genotyping and transcriptomics of the same cells enabling accurate and high-throughput screening of SNVs. Analysis of variants across the JAK1 gene with scSNV-seq demonstrates the importance of determining the precise genetic perturbation and accurately classifies clinically observed missense variants into three functional categories: benign, loss of function, and separation of function.

Mitochondrial and Endoplasmic Reticulum Stress Trigger Triglyceride Accumulation in Models of Parkinson's Disease Independent of Mutations in MAPT

The metabolic basis of Parkinson's disease pathology is poorly understood. However, the involvement of mitochondrial and endoplasmic reticulum stress in dopamine neurons in disease aetiology is well established. We looked at the effect of rotenone- and tunicamycin-induced mitochondrial and ER stress on the metabolism of wild type and microtubule-associated protein tau mutant dopamine neurons. Dopamine neurons derived from human isolated iPSCs were subjected to mitochondrial and ER stress using RT and TM, respectively. Comprehensive metabolite profiles were generated using a split phase extraction analysed by reversed phase lipidomics whilst the aqueous phase was measured using HILIC metabolomics. Mitochondrial and ER stress were both shown to cause significant dysregulation of metabolism with RT-induced stress producing a larger shift in the metabolic profile of both wild type and MAPT neurons. Detailed analysis showed that accumulation of triglycerides was a significant driver of metabolic dysregulation in response to both stresses in both genotypes. Whilst the consequence is similar, the mechanisms by which triglyceride accumulation occurs in dopamine neurons in response to mitochondrial and ER stress are very different. Thus, improving our understanding of how these mechanisms drive the observed triglyceride accumulation can potentially open up new therapeutic avenues.

Single-cell transcriptomics defines an improved, validated monoculture protocol for differentiation of human iPSC to microglia

There is increasing genetic evidence for the role of microglia in neurodegenerative diseases, including Alzheimer's, Parkinson's, and motor neuron disease. Therefore, there is a need to generate authentic in vitro models to study human microglial physiology. Various methods have been developed using human induced Pluripotent Stem Cells (iPSC) to generate microglia, however, systematic approaches to identify which media components are actually essential for functional microglia are mostly lacking. Here, we systematically assess medium components, coatings, and growth factors required for iPSC differentiation to microglia. Using single-cell RNA sequencing, qPCR, and functional assays, with validation across two labs, we have identified several medium components from previous protocols that are redundant and do not contribute to microglial identity. We provide an optimised, defined medium which produces both transcriptionally and functionally relevant microglia for modelling microglial physiology in neuroinflammation and for drug discovery.

Substantial somatic genomic variation and selection for BCOR mutations in human induced pluripotent stem cells

We explored human induced pluripotent stem cells (hiPSCs) derived from different tissues to gain insights into genomic integrity at single-nucleotide resolution. We used genome sequencing data from two large hiPSC repositories involving 696 hiPSCs and daughter subclones. We find ultraviolet light (UV)-related damage in ~72% of skin fibroblast-derived hiPSCs (F-hiPSCs), occasionally resulting in substantial mutagenesis (up to 15 mutations per megabase). We demonstrate remarkable genomic heterogeneity between independent F-hiPSC clones derived during the same round of reprogramming due to oligoclonal fibroblast populations. In contrast, blood-derived hiPSCs (B-hiPSCs) had fewer mutations and no UV damage but a high prevalence of acquired BCOR mutations (26.9% of lines). We reveal strong selection pressure for BCOR mutations in F-hiPSCs and B-hiPSCs and provide evidence that they arise in vitro. Directed differentiation of hiPSCs and RNA sequencing showed that BCOR mutations have functional consequences. Our work strongly suggests that detailed nucleotide-resolution characterization is essential before using hiPSCs.

Single-Cell Transcriptomics of Parkinson's Disease Human In Vitro Models Reveals Dopamine Neuron-Specific Stress Responses

The advent of induced pluripotent stem cell (iPSC)-derived neurons has revolutionized Parkinson's disease (PD) research, but single-cell transcriptomic analysis suggests unresolved cellular heterogeneity within these models. Here, we perform the largest single-cell transcriptomic study of human iPSC-derived dopaminergic neurons to elucidate gene expression dynamics in response to cytotoxic and genetic stressors. We identify multiple neuronal subtypes with transcriptionally distinct profiles and differential sensitivity to stress, highlighting cellular heterogeneity in dopamine in vitro models. We validate this disease model by showing robust expression of PD GWAS genes and overlap with postmortem adult substantia nigra neurons. Importantly, stress signatures are ameliorated using felodipine, an FDA-approved drug. Using isogenic SNCA-A53T mutants, we find perturbations in glycolysis, cholesterol metabolism, synaptic signaling, and ubiquitin-proteasomal degradation. Overall, our study reveals cell type-specific perturbations in human dopamine neurons, which will further our understanding of PD and have implications for cell replacement therapies.

The Ncoa7 locus regulates V-ATPase formation and function, neurodevelopment and behaviour

Members of the Tre2/Bub2/Cdc16 (TBC), lysin motif (LysM), domain catalytic (TLDc) protein family are associated with multiple neurodevelopmental disorders, although their exact roles in disease remain unclear. For example, nuclear receptor coactivator 7 (NCOA7) has been associated with autism, although almost nothing is known regarding the mode-of-action of this TLDc protein in the nervous system. Here we investigated the molecular function of NCOA7 in neurons and generated a novel mouse model to determine the consequences of deleting this locus in vivo. We show that NCOA7 interacts with the cytoplasmic domain of the vacuolar (V)-ATPase in the brain and demonstrate that this protein is required for normal assembly and activity of this critical proton pump. Neurons lacking Ncoa7 exhibit altered development alongside defective lysosomal formation and function; accordingly, Ncoa7 deletion animals exhibited abnormal neuronal patterning defects and a reduced expression of lysosomal markers. Furthermore, behavioural assessment revealed anxiety and social defects in mice lacking Ncoa7. In summary, we demonstrate that NCOA7 is an important V-ATPase regulatory protein in the brain, modulating lysosomal function, neuronal connectivity and behaviour; thus our study reveals a molecular mechanism controlling endolysosomal homeostasis that is essential for neurodevelopment.

INSIGHT: A population-scale COVID-19 testing strategy combining point-of-care diagnosis with centralized high-throughput sequencing

We present INSIGHT [isothermal NASBA (nucleic acid sequence-based amplification) sequencing-based high-throughput test], a two-stage coronavirus disease 2019 testing strategy, using a barcoded isothermal NASBA reaction. It combines point-of-care diagnosis with next-generation sequencing, aiming to achieve population-scale testing. Stage 1 allows a quick decentralized readout for early isolation of presymptomatic or asymptomatic patients. It gives results within 1 to 2 hours, using either fluorescence detection or a lateral flow readout, while simultaneously incorporating sample-specific barcodes. The same reaction products from potentially hundreds of thousands of samples can then be pooled and used in a highly multiplexed sequencing-based assay in stage 2. This second stage confirms the near-patient testing results and facilitates centralized data collection. The 95% limit of detection is <50 copies of viral RNA per reaction. INSIGHT is suitable for further development into a rapid home-based, point-of-care assay and is potentially scalable to the population level.

Screening for functional transcriptional and splicing regulatory variants with GenIE

Genome-wide association studies (GWAS) have identified numerous genetic loci underlying human diseases, but a fundamental challenge remains to accurately identify the underlying causal genes and variants. Here, we describe an arrayed CRISPR screening method, Genome engineering-based Interrogation of Enhancers (GenIE), which assesses the effects of defined alleles on transcription or splicing when introduced in their endogenous genomic locations. We use this sensitive assay to validate the activity of transcriptional enhancers and splice regulatory elements in human induced pluripotent stem cells (hiPSCs), and develop a software package (rgenie) to analyse the data. We screen the 99% credible set of Alzheimer's disease (AD) GWAS variants identified at the clusterin (CLU) locus to identify a subset of likely causal variants, and employ GenIE to understand the impact of specific mutations on splicing efficiency. We thus establish GenIE as an efficient tool to rapidly screen for the role of transcribed variants on gene expression.

The MITF-SOX10 regulated long non-coding RNA DIRC3 is a melanoma tumour suppressor

The MITF and SOX10 transcription factors regulate the expression of genes important for melanoma proliferation, invasion and metastasis. Despite growing evidence of the contribution of long noncoding RNAs (lncRNAs) in cancer, including melanoma, their functions within MITF-SOX10 transcriptional programmes remain poorly investigated. Here we identify 245 candidate melanoma associated lncRNAs whose loci are co-occupied by MITF-SOX10 and that are enriched at active enhancer-like regions. Our work suggests that one of these, Disrupted In Renal Carcinoma 3 (DIRC3), may be a clinically important MITF-SOX10 regulated tumour suppressor. DIRC3 depletion in human melanoma cells leads to increased anchorage-independent growth, a hallmark of malignant transformation, whilst melanoma patients classified by low DIRC3 expression have decreased survival. DIRC3 is a nuclear lncRNA that activates expression of its neighbouring IGFBP5 tumour suppressor through modulating chromatin structure and suppressing SOX10 binding to putative regulatory elements within the DIRC3 locus. In turn, DIRC3 dependent regulation of IGFBP5 impacts the expression of genes involved in cancer associated processes and is needed for DIRC3 control of anchorage-independent growth. Our work indicates that lncRNA components of MITF-SOX10 networks are an important new class of melanoma regulators and candidate therapeutic targets that can act not only as downstream mediators of MITF-SOX10 function but as feedback regulators of MITF-SOX10 activity.

Editing the Genome of Human Induced Pluripotent Stem Cells Using CRISPR/Cas9 Ribonucleoprotein Complexes

Genome editing using the CRISPR/Cas9 system has rapidly established itself as an essential tool in the genetic manipulation of many organisms, including human cell lines. Its application to human induced pluripotent stem cells (hiPSCs) allows for the generation of isogenic cell pairs that differ in a single genetic lesion, and therefore the identification and characterization of causal genetic variants. We describe a simple, effective approach to perform delicate manipulations of the genome of hiPSCs through delivery of Cas9 RNPs along with ssDNA oligonucleotide repair templates that can generate mutations in up to 98% of single cell clones and introduce single nucleotide changes at an efficiency of up to 40%. We describe our use of a T7 endonuclease assay to identify active guide RNAs, and a high-throughput sequencing genotyping strategy that allows the identification of correctly edited clones. We also present our experiences of generating single nucleotide changes at 15 sites, which show considerable variability between both guides and target sites in the efficiency at which such changes can be introduced.

Predicting the mutations generated by repair of Cas9-induced double-strand breaks

The DNA mutation produced by cellular repair of a CRISPR-Cas9-generated double-strand break determines its phenotypic effect. It is known that the mutational outcomes are not random, but depend on DNA sequence at the targeted location. Here we systematically study the influence of flanking DNA sequence on repair outcome by measuring the edits generated by >40,000 guide RNAs (gRNAs) in synthetic constructs. We performed the experiments in a range of genetic backgrounds and using alternative CRISPR-Cas9 reagents. In total, we gathered data for >109 mutational outcomes. The majority of reproducible mutations are insertions of a single base, short deletions or longer microhomology-mediated deletions. Each gRNA has an individual cell-line-dependent bias toward particular outcomes. We uncover sequence determinants of the mutations produced and use these to derive a predictor of Cas9 editing outcomes. Improved understanding of sequence repair will allow better design of gene editing experiments.

Alpha-synuclein induces the unfolded protein response in Parkinson's disease SNCA triplication iPSC-derived neurons

The recent generation of induced pluripotent stem cells (iPSCs) from a patient with Parkinson's disease (PD) resulting from triplication of the α-synuclein (SNCA) gene locus allows unprecedented opportunities to explore its contribution to the molecular pathogenesis of PD. We used the double-nicking CRISPR/Cas9 system to conduct site-specific mutagenesis of SNCA in these cells, generating an isogenic iPSC line with normalized SNCA gene dosage. Comparative gene expression analysis of neuronal derivatives from these iPSCs revealed an ER stress phenotype, marked by induction of the IRE1α/XBP1 axis of the unfolded protein response (UPR) and culminating in terminal UPR activation. Neuropathological analysis of post-mortem brain tissue demonstrated that pIRE1α is expressed in PD brains within neurons containing elevated levels of α-synuclein or Lewy bodies. Having used this pair of isogenic iPSCs to define this phenotype, these cells can be further applied in UPR-targeted drug discovery towards the development of disease-modifying therapeutics.

In situ functional dissection of RNA cis-regulatory elements by multiplex CRISPR-Cas9 genome engineering

RNA regulatory elements (RREs) are an important yet relatively under-explored facet of gene regulation. Deciphering the prevalence and functional impact of this post-transcriptional control layer requires technologies for disrupting RREs without perturbing cellular homeostasis. Here we describe genome-engineering based evaluation of RNA regulatory element activity (GenERA), a clustered regularly interspaced short palindromic repeats (CRISPR)-Cas9 platform for in situ high-content functional analysis of RREs. We use GenERA to survey the entire regulatory landscape of a 3′UTR, and apply it in a multiplex fashion to analyse combinatorial interactions between sets of miRNA response elements (MREs), providing strong evidence for cooperative activity. We also employ this technology to probe the functionality of an entire MRE network under cellular homeostasis, and show that high-resolution analysis of the GenERA dataset can be used to extract functional features of MREs. This study provides a genome editing-based multiplex strategy for direct functional interrogation of RNA cis-regulatory elements in a native cellular environment.

RNA regulatory elements (RREs) are important post-transcriptional control features but studying them requires disrupting their activity without disturbing cellular homeostasis. Here the authors present GenERA, a CRISPR-Cas9 screening platform of in situ analysis of native RREs.

Editing the genome of hiPSC with CRISPR/Cas9: disease models

The advent of human-induced pluripotent stem cell (hiPSC) technology has provided a unique opportunity to establish cellular models of disease from individual patients, and to study the effects of the underlying genetic aberrations upon multiple different cell types, many of which would not normally be accessible. Combining this with recent advances in genome editing techniques such as the clustered regularly interspaced short palindromic repeat (CRISPR) system has provided an ability to repair putative causative alleles in patient lines, or introduce disease alleles into a healthy “WT” cell line. This has enabled analysis of isogenic cell pairs that differ in a single genetic change, which allows a thorough assessment of the molecular and cellular phenotypes that result from this abnormality. Importantly, this establishes the true causative lesion, which is often impossible to ascertain from human genetic studies alone. These isogenic cell lines can be used not only to understand the cellular consequences of disease mutations, but also to perform high throughput genetic and pharmacological screens to both understand the underlying pathological mechanisms and to develop novel therapeutic agents to prevent or treat such diseases. In the future, optimising and developing such genetic manipulation technologies may facilitate the provision of cellular or molecular gene therapies, to intervene and ultimately cure many debilitating genetic disorders.

Most microRNAs in the single-cell alga Chlamydomonas reinhardtii are produced by Dicer-like 3-mediated cleavage of introns and untranslated regions of coding RNAs

We describe here a forward genetic screen to investigate the biogenesis, mode of action, and biological function of miRNA-mediated RNA silencing in the model algal species, Chlamydomonas reinhardtii. Among the mutants from this screen, there were three at Dicer-like 3 that failed to produce both miRNAs and siRNAs and others affecting diverse post-biogenesis stages of miRNA-mediated silencing. The DCL3-dependent siRNAs fell into several classes including transposon- and repeat-derived siRNAs as in higher plants. The DCL3-dependent miRNAs differ from those of higher plants, however, in that many of them are derived from mRNAs or from the introns of pre-mRNAs. Transcriptome analysis of the wild-type and dcl3 mutant strains revealed a further difference from higher plants in that the sRNAs are rarely negative switches of mRNA accumulation. The few transcripts that were more abundant in dcl3 mutant strains than in wild-type cells were not due to sRNA-targeted RNA degradation but to direct DCL3 cleavage of miRNA and siRNA precursor structures embedded in the untranslated (and translated) regions of the mRNAs. Our analysis reveals that the miRNA-mediated RNA silencing in C. reinhardtii differs from that of higher plants and informs about the evolution and function of this pathway in eukaryotes.

A genome-wide CRISPR library for high-throughput genetic screening in Drosophila cells

The simplicity of the CRISPR/Cas9 system of genome engineering has opened up the possibility of performing genome-wide targeted mutagenesis in cell lines, enabling screening for cellular phenotypes resulting from genetic aberrations. Drosophila cells have proven to be highly effective in identifying genes involved in cellular processes through similar screens using partial knockdown by RNAi. This is in part due to the lower degree of redundancy between genes in this organism, whilst still maintaining highly conserved gene networks and orthologs of many human disease-causing genes. The ability of CRISPR to generate genetic loss of function mutations not only increases the magnitude of any effect over currently employed RNAi techniques, but allows analysis over longer periods of time which can be critical for certain phenotypes. In this study, we have designed and built a genome-wide CRISPR library covering 13,501 genes, among which 8989 genes are targeted by three or more independent single guide RNAs (sgRNAs). Moreover, we describe strategies to monitor the population of guide RNAs by high throughput sequencing (HTS). We hope that this library will provide an invaluable resource for the community to screen loss of function mutations for cellular phenotypes, and as a source of guide RNA designs for future studies.

Understanding functional miRNA-target interactions in vivo by site-specific genome engineering

MicroRNA (miRNA) target recognition is largely dictated by short 'seed' sequences, and single miRNAs therefore have the potential to regulate a large number of genes. Understanding the contribution of specific miRNA-target interactions to the regulation of biological processes in vivo remains challenging. Here we use transcription activator-like effector nuclease (TALEN) and clustered regularly interspaced short palindromic repeat (CRISPR)/Cas9 technologies to interrogate the functional relevance of predicted miRNA response elements (MREs) to post-transcriptional silencing in zebrafish and Drosophila. We also demonstrate an effective strategy that uses CRISPR-mediated homology-directed repair with short oligonucleotide donors for the assessment of MRE activity in human cells. These methods facilitate analysis of the direct phenotypic consequences resulting from blocking specific miRNA-MRE interactions at any point during development.

Considerations when investigating lncRNA function in vivo

Although a small number of the vast array of animal long non-coding RNAs (lncRNAs) have known effects on cellular processes examined in vitro, the extent of their contributions to normal cell processes throughout development, differentiation and disease for the most part remains less clear. Phenotypes arising from deletion of an entire genomic locus cannot be unequivocally attributed either to the loss of the lncRNA per se or to the associated loss of other overlapping DNA regulatory elements. The distinction between cis- or trans-effects is also often problematic. We discuss the advantages and challenges associated with the current techniques for studying the in vivo function of lncRNAs in the light of different models of lncRNA molecular mechanism, and reflect on the design of experiments to mutate lncRNA loci. These considerations should assist in the further investigation of these transcriptional products of the genome.

DOI:http://dx.doi.org/10.7554/eLife.03058.001

Mutagenesis and homologous recombination in Drosophila cell lines using CRISPR/Cas9

We have applied the CRISPR/Cas9 system to Drosophila S2 cells to generate targeted genetic mutations in more than 85% of alleles. By targeting a constitutive exon of the AGO1 gene, we demonstrate homozygous mutation in up to 82% of cells, thereby allowing the study of genetic knockouts in a Drosophila cell line for the first time. We have shown that homologous gene targeting is possible at 1–4% efficiency using this system, allowing for the construction of defined insertions and deletions. We demonstrate that a 1 kb homology arm length is optimal for integration by homologous gene targeting, and demonstrate its efficacy by tagging the endogenous AGO1 protein. This technology enables controlled genetic manipulation in Drosophila cell lines, and its simplicity offers the opportunity to study cellular phenotypes genome-wide.

CRISPR/Cas9 and genome editing in Drosophila

Recent advances in our ability to design DNA binding factors with specificity for desired sequences have resulted in a revolution in genetic engineering, enabling directed changes to the genome to be made relatively easily. Traditional techniques for generating genetic mutations in most organisms have relied on selection from large pools of randomly induced mutations for those of particular interest, or time-consuming gene targeting by homologous recombination. Drosophila melanogaster has always been at the forefront of genetic analysis, and application of these new genome editing techniques to this organism will revolutionise our approach to performing analysis of gene function in the future. We discuss the recent techniques that apply the CRISPR/Cas9 system to Drosophila, highlight potential uses for this technology and speculate upon the future of genome engineering in this model organism.

Identification and properties of 1,119 candidate lincRNA loci in the Drosophila melanogaster genome

The functional repertoire of long intergenic noncoding RNA (lincRNA) molecules has begun to be elucidated in mammals. Determining the biological relevance and potential gene regulatory mechanisms of these enigmatic molecules would be expedited in a more tractable model organism, such as Drosophila melanogaster. To this end, we defined a set of 1,119 putative lincRNA genes in D. melanogaster using modENCODE whole transcriptome (RNA-seq) data. A large majority (1.1 of 1.3 Mb; 85%) of these bases were not previously reported by modENCODE as being transcribed. Significant selective constraint on the sequences of these loci predicts that virtually all have sustained functionality across the Drosophila clade. We observe biases in lincRNA genomic locations and expression profiles that are consistent with some of these lincRNAs being involved in the regulation of neighboring protein-coding genes with developmental functions. We identify lincRNAs that may be important in the developing nervous system and in male-specific organs, such as the testes. LincRNA loci were also identified whose positions, relative to nearby protein-coding loci, are equivalent between D. melanogaster and mouse. This study predicts that the genomes of not only vertebrates, such as mammals, but also an invertebrate (fruit fly) harbor large numbers of lincRNA loci. Our findings now permit exploitation of Drosophila genetics for the investigation of lincRNA mechanisms, including lincRNAs with potential functional analogues in mammals.