CRISPR Screens to Discover Functional Noncoding Elements

Cross-Species Extrapolation of Biological Data to Guide the Environmental Safety Assessment of Pharmaceuticals-The State of the Art and Future Priorities

The extrapolation of biological data across species is a key aspect of biomedical research and drug development. In this context, comparative biology considerations are applied with the goal of understanding human disease and guiding the development of effective and safe medicines. However, the widespread occurrence of pharmaceuticals in the environment and the need to assess the risk posed to wildlife have prompted a renewed interest in the extrapolation of pharmacological and toxicological data across the entire tree of life. To address this challenge, a biological "read-across" approach, based on the use of mammalian data to inform toxicity predictions in wildlife species, has been proposed as an effective way to streamline the environmental safety assessment of pharmaceuticals. Yet, how effective has this approach been, and are we any closer to being able to accurately predict environmental risk based on known human risk? We discuss the main theoretical and experimental advancements achieved in the last 10 years of research in this field. We propose that a better understanding of the functional conservation of drug targets across species and of the quantitative relationship between target modulation and adverse effects should be considered as future research priorities. This pharmacodynamic focus should be complemented with the application of higher-throughput experimental and computational approaches to accelerate the prediction of internal exposure dynamics. The translation of comparative (eco)toxicology research into real-world applications, however, relies on the (limited) availability of experts with the skill set needed to navigate the complexity of the problem; hence, we also call for synergistic multistakeholder efforts to support and strengthen comparative toxicology research and education at a global level. Environ Toxicol Chem 2024;43:513-525. © 2023 The Authors. Environmental Toxicology and Chemistry published by Wiley Periodicals LLC on behalf of SETAC.

The immune system of prokaryotes: potential applications and implications for gene editing

Gene therapy has revolutionized the treatment of genetic diseases. Spearheading this revolution are sophisticated genome editing methods such as TALENs, ZFNs, and CRISPR-Cas, which trace their origins back to prokaryotic immune systems. Prokaryotes have developed various antiviral defense systems to combat viral attacks and the invasion of genetic elements. The comprehension of these defense mechanisms has paved the way for the development of indispensable tools in molecular biology. Among them, restriction endonuclease originates from the innate immune system of bacteria. The CRISPR-Cas system, a widely applied genome editing technology, is derived from the prokaryotic adaptive immune system. Single-base editing is a precise editing tool based on CRISPR-Cas system that involves deamination of target base. It is worth noting that prokaryotes possess deamination enzymes as part of their defense arsenal over foreign genetic material. Furthermore, prokaryotic Argonauts (pAgo) proteins, also function in anti-phage defense, play an important role in complementing the CRISPR-Cas system by addressing certain limitations it may have. Recent studies have also shed light on the significance of Retron, a reverse transcription transposon previously showed potential in genome editing, has also come to light in the realm of prokaryotic immunity. These noteworthy findings highlight the importance of studying prokaryotic immune system for advancing genome editing techniques. Here, both the origin of prokaryotic immunity underlying aforementioned genome editing tools, and potential applications of deaminase, pAgo protein and reverse transcriptase in genome editing among prokaryotes were introduced, thus emphasizing the fundamental mechanism and significance of prokaryotic immunity.

Advancements in CRISPR screens for the development of cancer immunotherapy strategies

CRISPR screen technology enables systematic and scalable interrogation of gene function by using the CRISPR-Cas9 system to perturb gene expression. In the field of cancer immunotherapy, this technology has empowered the discovery of genes, biomarkers, and pathways that regulate tumor development and progression, immune reactivity, and the effectiveness of immunotherapeutic interventions. By conducting large-scale genetic screens, researchers have successfully identified novel targets to impede tumor growth, enhance anti-tumor immune responses, and surmount immunosuppression within the tumor microenvironment (TME). Here, we present an overview of CRISPR screens conducted in tumor cells for the purpose of identifying novel therapeutic targets. We also explore the application of CRISPR screens in immune cells to propel the advancement of cell-based therapies, encompassing T cells, natural killer cells, dendritic cells, and macrophages. Furthermore, we outline the crucial components necessary for the successful implementation of immune-specific CRISPR screens and explore potential directions for future research.

CRISPR/dCas9 Tools: Epigenetic Mechanism and Application in Gene Transcriptional Regulation

CRISPR/Cas9-mediated cleavage of DNA, which depends on the endonuclease activity of Cas9, has been widely used for gene editing due to its excellent programmability and specificity. However, the changes to the DNA sequence that are mediated by CRISPR/Cas9 affect the structures and stability of the genome, which may affect the accuracy of results. Mutations in the RuvC and HNH regions of the Cas9 protein lead to the inactivation of Cas9 into dCas9 with no endonuclease activity. Despite the loss of endonuclease activity, dCas9 can still bind the DNA strand using guide RNA. Recently, proteins with active/inhibitory effects have been linked to the end of the dCas9 protein to form fusion proteins with transcriptional active/inhibitory effects, named CRISPRa and CRISPRi, respectively. These CRISPR tools mediate the transcription activity of protein-coding and non-coding genes by regulating the chromosomal modification states of target gene promoters, enhancers, and other functional elements. Here, we highlight the epigenetic mechanisms and applications of the common CRISPR/dCas9 tools, by which we hope to provide a reference for future related gene regulation, gene function, high-throughput target gene screening, and disease treatment.

Discovery of target genes and pathways at GWAS loci by pooled single-cell CRISPR screens

Most variants associated with complex traits and diseases identified by genome-wide association studies (GWAS) map to noncoding regions of the genome with unknown effects. Using ancestrally diverse, biobank-scale GWAS data, massively parallel CRISPR screens, and single-cell transcriptomic and proteomic sequencing, we discovered 124 cis-target genes of 91 noncoding blood trait GWAS loci. Using precise variant insertion through base editing, we connected specific variants with gene expression changes. We also identified trans-effect networks of noncoding loci when cis target genes encoded transcription factors or microRNAs. Networks were themselves enriched for GWAS variants and demonstrated polygenic contributions to complex traits. This platform enables massively parallel characterization of the target genes and mechanisms of human noncoding variants in both cis and trans.

Biallelic and gene-wide genomic substitution for endogenous intron and retroelement mutagenesis in human cells

Functional annotation of the vast noncoding landscape of the diploid human genome still remains a major challenge of genomic research. An efficient, scarless, biallelic, and gene-wide mutagenesis approach is needed for direct investigation of the functional significance of endogenous long introns in gene regulation. Here we establish a genome substitution platform, the Universal Knock-in System or UKiS, that meets these requirements. For proof of concept, we first used UKiS on the longest intron of TP53 in the pseudo-diploid cell line HCT116. Complete deletion of the intron, its substitution with mouse and zebrafish syntenic introns, and specific removal of retrotransposon-derived elements (retroelements) were all efficiently and accurately achieved in both alleles, revealing a suppressive role of intronic Alu elements in TP53 expression. We also used UKiS for TP53 intron deletion in human induced pluripotent stem cells without losing their stemness. Furthermore, UKiS enabled biallelic removal of all introns from three human gene loci of ~100 kb and longer to demonstrate that intron requirements for transcriptional activities vary among genes. UKiS is a standard platform with which to pursue the design of noncoding regions for genome writing in human cells.

Enhancer transcription detected in the nascent transcriptomic landscape of bread wheat

The precise spatiotemporal gene expression is orchestrated by enhancers that lack general sequence features and thus are difficult to be computationally identified. By nascent RNA sequencing combined with epigenome profiling, we detect active transcription of enhancers from the complex bread wheat genome. We find that genes associated with transcriptional enhancers are expressed at significantly higher levels, and enhancer RNA is more precise and robust in predicting enhancer activity compared to chromatin features. We demonstrate that sub-genome-biased enhancer transcription could drive sub-genome-biased gene expression. This study highlights enhancer transcription as a hallmark in regulating gene expression in wheat.

CRISPR-Based Approaches for the High-Throughput Characterization of Long Non-Coding RNAs

Over the last decade, tens of thousands of new long non-coding RNAs (lncRNAs) have been identified in the human genome. Nevertheless, except for a handful of genes, the genetic characteristics and functions of most of these lncRNAs remain elusive; this is partially due to their relatively low expression, high tissue specificity, and low conservation across species. A major limitation for determining the function of lncRNAs was the lack of methodologies suitable for studying these genes. The recent development of CRISPR/Cas9 technology has opened unprecedented opportunities to uncover the genetic and functional characteristics of the non-coding genome via targeted and high-throughput approaches. Specific CRISPR/Cas9-based approaches were developed to target lncRNA loci. Some of these approaches involve modifying the sequence, but others were developed to study lncRNAs by inducing transcriptional and epigenetic changes. The discovery of other programable Cas proteins broaden our possibilities to target RNA molecules with greater precision and accuracy. These approaches allow for the knock-down and characterization of lncRNAs. Here, we review how various CRISPR-based strategies have been used to characterize lncRNAs with important functions in different biological contexts and how these approaches can be further utilized to improve our understanding of the non-coding genome.

Blockade of Autocrine CCL5 Responses Inhibits Zika Virus Persistence and Spread in Human Brain Microvascular Endothelial Cells

Zika virus (ZIKV) is a neurovirulent flavivirus that uniquely causes fetal microcephaly, is sexually transmitted, and persists in patients for up to 6 months. ZIKV persistently infects human brain microvascular endothelial cells (hBMECs) that form the blood-brain barrier (BBB) and enables viral spread to neuronal compartments. We found that CCL5, a chemokine with prosurvival effects on immune cells, was highly secreted by ZIKV-infected hBMECs. Although roles for CCL5 in endothelial cell (EC) survival remain unknown, the presence of the CCL5 receptors CCR3 and CCR5 on ECs suggested that CCL5 could promote ZIKV persistence in hBMECs. We found that exogenous CCL5 induced extracellular signal-regulated kinase 1/2 (ERK1/2) phosphorylation in hBMECs and that ERK1/2 cell survival signaling was similarly activated by ZIKV infection. Neutralizing antibodies to CCL5, CCR3, or CCR5 inhibited persistent ZIKV infection of hBMECs. While knockout (KO) of CCL5 failed to prevent ZIKV infection of hBMECs, at 3 days postinfection (dpi), we observed a >90% reduction in ZIKV-infected CCL5-KO hBMECs and a multilog reduction in ZIKV titers. In contrast, the addition of CCL5 to CCL5-KO hBMECs dose-dependently rescued ZIKV persistence in hBMECs. Inhibiting CCL5 responses using CCR3 (UCB35625) and CCR5 (maraviroc) receptor antagonists reduced the number of ZIKV-infected hBMECs and ZIKV titers (50% inhibitory concentrations [IC₅₀s] of 2.5 to 12 μM), without cytotoxicity (50% cytotoxic concentration [CC₅₀] of >80 μM). These findings demonstrate that ZIKV-induced CCL5 directs autocrine CCR3/CCR5 activation of ERK1/2 survival responses that are required for ZIKV to persistently infect hBMECs. Our results establish roles for CCL5 in ZIKV persistence and suggest the potential for CCL5 receptor antagonists to therapeutically inhibit ZIKV spread and neurovirulence.

Computational prediction of CRISPR-impaired non-coding regulatory regions

Genome-wide CRISPR screens are becoming more widespread and allow the simultaneous interrogation of thousands of genomic regions. Although recent progress has been made in the analysis of CRISPR screens, it is still an open problem how to interpret CRISPR mutations in non-coding regions of the genome. Most of the tools concentrate on the interpretation of mutations introduced in gene coding regions. We introduce a computational pipeline that uses epigenomic information about regulatory elements for the interpretation of CRISPR mutations in non-coding regions. We illustrate our analysis protocol on the analysis of a genome-wide CRISPR screen in hTERT-RPE1 cells and reveal novel regulatory elements that mediate chemoresistance against doxorubicin in these cells. We infer links to established and to novel chemoresistance genes. Our analysis protocol is general and can be applied on any cell type and with different CRISPR enzymes.

Enhancer redundancy in development and disease

Shadow enhancers are seemingly redundant transcriptional cis-regulatory elements that regulate the same gene and drive overlapping expression patterns. Recent studies have shown that shadow enhancers are remarkably abundant and control most developmental gene expression in both invertebrates and vertebrates, including mammals. Shadow enhancers might provide an important mechanism for buffering gene expression against mutations in non-coding regulatory regions of genes implicated in human disease. Technological advances in genome editing and live imaging have shed light on how shadow enhancers establish precise gene expression patterns and confer phenotypic robustness. Shadow enhancers can interact in complex ways and may also help to drive the formation of transcriptional hubs within the nucleus. Despite their apparent redundancy, the prevalence and evolutionary conservation of shadow enhancers underscore their key role in emerging metazoan gene regulatory networks.

Integrative identification of the pathogenic role of a novel G6PD missense mutation c.697G>C

Background

Glucose-6-phosphate dehydrogenase (G6PD) deficiency is a hereditary disease caused by pathogenic mutations of G6PD. While most of the pathogenic variants of G6PD have been annotated, hemolysis of unknown etiology but analogous to that in G6PD deficiency persists, implying the existence of undocumented pathogenic variants. In our previous study, we reported four novel G6PD variants in China, for which the pathogenicity remains to be verified.

Methods

The variants were verified by exogenous expression in HEK-293 cells, and their functions were predicted by PolyPhen-2 and SIFT. The CRISPR/Cas9 system was exploited to edit the G6PD c.697G>C variant in HEK-293 cells and K562 cells. The expression of G6PD was detected by quantitative PCR (qPCR) and western blotting. The cell growth capacity was detected by the CCK-8 assay and crystal violet staining. The G6PD enzyme activity was reflected by the G6P/6PG ratio test. The apoptosis of cells was detected by Annexin V-APC/7-AAD staining. The secondary and crystallographic structures were denoted according to the literature and PyMOL software. The G6PD protein was purified from lysis of transformed Escherichia coli (E. coli) cell with Ni-charged Resin Column. The enzymatic activity was detected at different temperatures.

Results

The G6PD activity of exogenous G6PD c.697G>C in HEK-293 cells was significantly lower than that of wild type (WT) G6PD, a finding that was consistent with the observation in clinical samples. The functional predictions conducted by different algorithms indicated the damage role of the G6PD c.697G>C variant in its enzymatic activity. We recapitulated the G6PD c.697G>C variant both in HEK-293 cells and K562 cells by adapting the CRISPR/Cas9 strategy. Using distinct cell lines expressing the G6PD c.697G>C variant endogenously, we confirmed the deteriorative role of the G6PD c.697G>C variant in its enzymatic activity. Regarding the secondary and crystallographic structure, we found a mutated amino acid approaching the structural NADP⁺ binding site. Finally, we demonstrated the c.697G>C variant compromised the thermal stability of G6PD protein.

Conclusions

Our data delineated the pathogenic role of G6PD c.697G>C variant for G6PD deficiency, implying the wide usage of CRISPR/Cas9 for genetic disease research.

Towards a comprehensive catalogue of validated and target-linked human enhancers

The human gene catalogue is essentially complete, but we lack an equivalently vetted inventory of bona fide human enhancers. Hundreds of thousands of candidate enhancers have been nominated via biochemical annotations; however, only a handful of these have been validated and confidently linked to their target genes. Here we review emerging technologies for discovering, characterizing and validating human enhancers at scale. We furthermore propose a new framework for operationally defining enhancers that accommodates the heterogeneous and complementary results that are emerging from reporter assays, biochemical measurements and CRISPR screens.

Functional genetic analysis in a jawless vertebrate, the sea lamprey: insights into the developmental evolution of early vertebrates

Lampreys and hagfishes are the only surviving relicts of an ancient but ecologically dominant group of jawless fishes that evolved in the seas of the Cambrian era over half a billion years ago. Because of their phylogenetic position as the sister group to all other vertebrates (jawed vertebrates), comparisons of embryonic development between jawless and jawed vertebrates offers researchers in the field of evolutionary developmental biology the unique opportunity to address fundamental questions related to the nature of our earliest vertebrate ancestors. Here, we describe how genetic analysis of embryogenesis in the sea lamprey (Petromyzon marinus) has provided insight into the origin and evolution of developmental-genetic programs in vertebrates. We focus on recent work involving CRISPR/Cas9-mediated genome editing to study gene regulatory mechanisms involved in the development and evolution of neural crest cells and new cell types in the vertebrate nervous system, and transient transgenic assays that have been instrumental in dissecting the evolution of cis-regulatory control of gene expression in vertebrates. Finally, we discuss the broad potential for these functional genomic tools to address previously unanswerable questions related to the evolution of genomic regulatory mechanisms as well as issues related to invasive sea lamprey population control.

CNEr: A toolkit for exploring extreme noncoding conservation

Conserved Noncoding Elements (CNEs) are elements exhibiting extreme noncoding conservation in Metazoan genomes. They cluster around developmental genes and act as long-range enhancers, yet nothing that we know about their function explains the observed conservation levels. Clusters of CNEs coincide with topologically associating domains (TADs), indicating ancient origins and stability of TAD locations. This has suggested further hypotheses about the still elusive origin of CNEs, and has provided a comparative genomics-based method of estimating the position of TADs around developmentally regulated genes in genomes where chromatin conformation capture data is missing. To enable researchers in gene regulation and chromatin biology to start deciphering this phenomenon, we developed CNEr, a R/Bioconductor toolkit for large-scale identification of CNEs and for studying their genomic properties. We apply CNEr to two novel genome comparisons—fruit fly vs tsetse fly, and two sea urchin genomes—and report novel insights gained from their analysis. We also show how to reveal interesting characteristics of CNEs by coupling CNEr with existing Bioconductor packages. CNEr is available at Bioconductor (https://bioconductor.org/packages/CNEr/) and maintained at github (https://github.com/ge11232002/CNEr).

Discovery of gene regulatory elements through a new bioinformatics analysis of haploid genetic screens

The systematic identification of regulatory elements that control gene expression remains a challenge. Genetic screens that use untargeted mutagenesis have the potential to identify protein-coding genes, non-coding RNAs and regulatory elements, but their analysis has mainly focused on identifying the former two. To identify regulatory elements, we conducted a new bioinformatics analysis of insertional mutagenesis screens interrogating WNT signaling in haploid human cells. We searched for specific patterns of retroviral gene trap integrations (used as mutagens in haploid screens) in short genomic intervals overlapping with introns and regions upstream of genes. We uncovered atypical patterns of gene trap insertions that were not predicted to disrupt coding sequences, but caused changes in the expression of two key regulators of WNT signaling, suggesting the presence of cis-regulatory elements. Our methodology extends the scope of haploid genetic screens by enabling the identification of regulatory elements that control gene expression.

Whole Genome Sequencing-Based Discovery of Structural Variants in Glioblastoma

Next-generation DNA sequencing (NGS) technologies are currently being applied in both research and clinical settings for the understanding and management of disease. The goal is to use high-throughput sequencing to identify specific variants that drive tumorigenesis within each individual's tumor genomic profile. The significance of copy number and structural variants in glioblastoma makes it essential to broaden the search beyond oncogenic single nucleotide variants toward whole genome profiles of genetic aberrations that may contribute to disease progression. The heterogeneity of glioblastoma and its variability of cancer driver mutations necessitate a more robust examination of a patient's tumor genome. Here, we present patient whole genome sequencing (WGS) information to identify oncogenic structural variants that may contribute to glioblastoma pathogenesis. We provide WGS protocols and bioinformatics approaches to identify copy number and structural variations in 41 glioblastoma patient samples. We present how WGS can identify structural diversity within glioblastoma samples. We specifically show how to apply current bioinformatics tools to detect EGFR variants and other structural aberrations from DNA whole genome sequencing and how to validate those variants within the laboratory. These comprehensive WGS protocols can provide additional information directing more precise therapeutic options in the treatment of glioblastoma.

CRISPR-Cas9: A New Addition to the Drug Metabolism and Disposition Tool Box

Clustered regularly interspaced short palindromic repeats (CRISPR)-CRISPR associated protein 9 (Cas9), i.e., CRISPR-Cas9, has been extensively used as a gene-editing technology during recent years. Unlike earlier technologies for gene editing or gene knockdown, such as zinc finger nucleases and RNA interference, CRISPR-Cas9 is comparably easy to use, affordable, and versatile. Recently, CRISPR-Cas9 has been applied in studies of drug absorption, distribution, metabolism, and excretion (ADME) and for ADME model generation. To date, about 50 papers have been published describing in vitro or in vivo CRISPR-Cas9 gene editing of ADME and ADME-related genes. Twenty of these papers describe gene editing of clinically relevant genes, such as ATP-binding cassette drug transporters and cytochrome P450 drug-metabolizing enzymes. With CRISPR-Cas9, the ADME tool box has been substantially expanded. This new technology allows us to develop better and more predictive in vitro and in vivo ADME models and map previously underexplored ADME genes and gene families. In this mini-review, we give an overview of the CRISPR-Cas9 technology and summarize recent applications of CRISPR-Cas9 within the ADME field. We also speculate about future applications of CRISPR-Cas9 in ADME research.

The accessible chromatin landscape during conversion of human embryonic stem cells to trophoblast by bone morphogenetic protein 4

Human embryonic stem cells (hESCs) exposed to the growth factor bone morphogenetic protein 4 (BMP4) in the absence of FGF2 have been used as a model to study the development of placental development. However, little is known about the cis-regulatory mechanisms underlying this important process. In this study, we used the public available chromatin accessibility data of hESC H1 cells and BMP4-induced trophoblast (TB) cell lines to identify DNase I hypersensitive sites (DHSs) in the two cell lines, as well as the transcription factor (TF) binding sites within the DHSs. By comparing read profiles in H1 and TB, we identified 17 472 TB-specific DHSs. The TB-specific DHSs are enriched in terms of "blood vessel" and "trophectoderm," consisting of TF motifs family: Leucine Zipper, Helix-Loop-Helix, GATA, and ETS. To validate differential expression of the TFs binding to these motifs, we analyzed public available RNA-seq and microarray data in the same context. Finally, by integrating the protein-protein interaction data, we constructed a TF network for placenta development and identified top 20 key TFs through centrality analysis in the network. Our results indicate BMP4-induced TB system provided an invaluable model for the study of TB development and highlighted novel candidate genes in placenta development in human.

Advancing Stem Cell Models of Alpha-Synuclein Gene Regulation in Neurodegenerative Disease

Alpha-synuclein (non A4 component of amyloid precursor, SNCA, NM_000345.3) plays a central role in the pathogenesis of Parkinson's disease (PD) and related Lewy body disorders such as Parkinson's disease dementia, Lewy body dementia, and multiple system atrophy. Since its discovery as a disease-causing gene in 1997, alpha-synuclein has been a central point of scientific interest both at the protein and gene level. Mutations, including copy number variants, missense mutations, short structural variants, and single nucleotide polymorphisms, can be causative for PD and affect conformational changes of the protein, can contribute to changes in expression of alpha-synuclein and its isoforms, and can influence regulation of temporal as well as spatial levels of alpha-synuclein in different tissues and cell types. A lot of progress has been made to understand both the physiological transcriptional and epigenetic regulation of the alpha-synuclein gene and whether changes in transcriptional regulation could lead to disease and neurodegeneration in PD and related alpha-synucleinopathies. Although the histopathological changes in these neurodegenerative disorders are similar, the temporal and spatial presentation and progression distinguishes them which could be in part due to changes or disruption of transcriptional regulation of alpha-synuclein. In this review, we describe different genetic alterations that contribute to PD and neurodegenerative conditions and review aspects of transcriptional regulation of the alpha-synuclein gene in the context of the development of PD. New technologies, advanced gene engineering and stem cell modeling, are on the horizon to shed further light on a better understanding of gene regulatory processes and exploit them for therapeutic developments.

A High-Resolution Genome-Wide CRISPR/Cas9 Viability Screen Reveals Structural Features and Contextual Diversity of the Human Cell-Essential Proteome

To interrogate genes essential for cell growth, proliferation and survival in human cells, we carried out a genome-wide clustered regularly interspaced short palindromic repeat (CRISPR)/Cas9 screen in a B-cell lymphoma line using a custom extended-knockout (EKO) library of 278,754 single-guide RNAs (sgRNAs) that targeted 19,084 RefSeq genes, 20,852 alternatively spliced exons, and 3,872 hypothetical genes. A new statistical analysis tool called robust analytics and normalization for knockout screens (RANKS) identified 2,280 essential genes, 234 of which were unique. Individual essential genes were validated experimentally and linked to ribosome biogenesis and stress responses. Essential genes exhibited a bimodal distribution across 10 different cell lines, consistent with a continuous variation in essentiality as a function of cell type. Genes essential in more lines had more severe fitness defects and encoded the evolutionarily conserved structural cores of protein complexes, whereas genes essential in fewer lines formed context-specific modules and encoded subunits at the periphery of essential complexes. The essentiality of individual protein residues across the proteome correlated with evolutionary conservation, structural burial, modular domains, and protein interaction interfaces. Many alternatively spliced exons in essential genes were dispensable and were enriched for disordered regions. Fitness defects were observed for 44 newly evolved hypothetical reading frames. These results illuminate the contextual nature and evolution of essential gene functions in human cells.

Quantifying the regulatory effect size of cis-acting genetic variation using allelic fold change

Mapping cis-acting expression quantitative trait loci (cis-eQTL) has become a popular approach for characterizing proximal genetic regulatory variants. In this paper, we describe and characterize log allelic fold change (aFC), the magnitude of expression change associated with a given genetic variant, as a biologically interpretable unit for quantifying the effect size of cis-eQTLs and a mathematically convenient approach for systematic modeling of cis-regulation. This measure is mathematically independent from expression level and allele frequency, additive, applicable to multiallelic variants, and generalizable to multiple independent variants. We provide efficient tools and guidelines for estimating aFC from both eQTL and allelic expression data sets and apply it to Genotype Tissue Expression (GTEx) data. We show that aFC estimates independently derived from eQTL and allelic expression data are highly consistent, and identify technical and biological correlates of eQTL effect size. We generalize aFC to analyze genes with two eQTLs in GTEx and show that in nearly all cases the two eQTLs act independently in regulating gene expression. In summary, aFC is a solid measure of cis-regulatory effect size that allows quantitative interpretation of cellular regulatory events from population data, and it is a valuable approach for investigating novel aspects of eQTL data sets.

Using Zebrafish to Test the Genetic Basis of Human Craniofacial Diseases

Genome-wide association studies (GWASs) opened an innovative and productive avenue to investigate the molecular basis of human craniofacial disease. However, GWASs identify candidate genes only; they do not prove that any particular one is the functional villain underlying disease or just an unlucky genomic bystander. Genetic manipulation of animal models is the best approach to reveal which genetic loci identified from human GWASs are functionally related to specific diseases. The purpose of this review is to discuss the potential of zebrafish to resolve which candidate genetic loci are mechanistic drivers of craniofacial diseases. Many anatomic, embryonic, and genetic features of craniofacial development are conserved among zebrafish and mammals, making zebrafish a good model of craniofacial diseases. Also, the ability to manipulate gene function in zebrafish was greatly expanded over the past 20 y, enabling systems such as Gateway Tol2 and CRISPR-Cas9 to test gain- and loss-of-function alleles identified from human GWASs in coding and noncoding regions of DNA. With the optimization of genetic editing methods, large numbers of candidate genes can be efficiently interrogated. Finding the functional villains that underlie diseases will permit new treatments and prevention strategies and will increase understanding of how gene pathways operate during normal development.

Functional interrogation of non-coding DNA through CRISPR genome editing

Methodologies to interrogate non-coding regions have lagged behind coding regions despite comprising the vast majority of the genome. However, the rapid evolution of clustered regularly interspaced short palindromic repeats (CRISPR)-based genome editing has provided a multitude of novel techniques for laboratory investigation including significant contributions to the toolbox for studying non-coding DNA. CRISPR-mediated loss-of-function strategies rely on direct disruption of the underlying sequence or repression of transcription without modifying the targeted DNA sequence. CRISPR-mediated gain-of-function approaches similarly benefit from methods to alter the targeted sequence through integration of customized sequence into the genome as well as methods to activate transcription. Here we review CRISPR-based loss- and gain-of-function techniques for the interrogation of non-coding DNA.

CRISPR/Cas9-mediated deletion of lncRNA Gm26878 in the distant Foxf1 enhancer region

Recent genome editing techniques, including CRISPR mutagenesis screens, offer unparalleled opportunities to study the regulatory non-coding genomic regions, enhancers, promoters, and functional non-coding RNAs. Heterozygous point mutations in FOXF1 and genomic deletion copy-number variants at chromosomal region 16q24.1 involving FOXF1 or its regulatory region mapping ~300 kb upstream of FOXF1 and leaving it intact have been identified in the vast majority of patients with a lethal neonatal lung disease, alveolar capillary dysplasia with misalignment of pulmonary veins (ACDMPV). Homozygous Foxf1 ^-/- mice have been shown to die by embryonic day 8.5 because of defects in the development of extraembryonic and lateral mesoderm-derived tissues, whereas heterozygous Foxf1 ^+/- mice exhibit features resembling ACDMPV. We have previously defined a human lung-specific enhancer region encoding two long non-coding RNAs, LINC01081 and LINC01082, expressed in the lungs. To investigate the biological significance of lncRNAs in the Foxf1 enhancer region, we have generated a CRISPR/Cas9-mediated ~2.4 kb deletion involving the entire lncRNA-encoding gene Gm26878, located in the mouse region syntenic with the human Foxf1 upstream enhancer. Very recently, this mouse genomic region has been shown to function as a Foxf1 enhancer. Our results indicate that homozygous loss of Gm26878 is neonatal lethal with low penetrance. No changes in Foxf1 expression were observed, suggesting that the regulation of Foxf1 expression differs between mouse and human.