Publications

The ENCODE Consortium’s efforts to annotate noncoding cis-regulatory elements (CREs) have advanced our understanding of gene regulatory landscapes. Pooled, noncoding CRISPR screens offer a systematic approach to investigate cis-regulatory mechanisms. The ENCODE4 Functional Characterization Centers conducted 108 screens in human cell lines, comprising >540,000 perturbations across 24.85 megabases of the genome. Using 332 functionally confirmed CRE–gene links in K562 cells, we established guidelines for screening endogenous noncoding elements with CRISPR interference (CRISPRi), including accurate detection of CREs that exhibit variable, often low, transcriptional effects. Benchmarking five screen analysis tools, we find that CASA produces the most conservative CRE calls and is robust to artifacts of low-specificity single guide RNAs. We uncover a subtle DNA strand bias for CRISPRi in transcribed regions with implications for screen design and analysis. Together, we provide an accessible data resource, predesigned single guide RNAs for targeting 3,275,697 ENCODE SCREEN candidate CREs with CRISPRi and screening guidelines to accelerate functional characterization of the noncoding genome.

Conserved genomic sequences disrupted in humans may underlie uniquely human phenotypic traits. We identified and characterized 10,032 human-specific conserved deletions (hCONDELs). These short (average 2.56 base pairs) deletions are enriched for human brain functions across genetic, epigenomic, and transcriptomic datasets. Using massively parallel reporter assays in six cell types, we discovered 800 hCONDELs conferring significant differences in regulatory activity, half of which enhance rather than disrupt regulatory function. We highlight several hCONDELs with putative human-specific effects on brain development, including HDAC5, CPEB4, and PPP2CA. Reverting an hCONDEL to the ancestral sequence alters the expression of LOXL2 and developmental genes involved in myelination and synaptic function. Our data provide a rich resource to investigate the evolutionary mechanisms driving new traits in humans and other species.

The Fuller’s rose weevil Naupactus cervinus has become a globally invasive pest of several ornamental plants and fruit trees. This weevil has caused severe economic losses, and eggs laid on fruit are a quarantine barrier for several countries’ exports of fruit to markets in East Asia. Previous studies revealed that a genotype with high colonization ability successfully expanded throughout most continents, even in areas of inadequate environmental conditions, where this insect performs unexpectedly well, and that parthenogenesis may have helped to expand its geographic range. Pre-existing variation might have been eroded either by natural selection, leading to fixation of this variant able to cope with different environmental conditions to those in the native range, or by repeated bottlenecks during the process of invasion. To better understand the ecology of this invasive pest, we enlarged the sampling in areas of introduction, such as the southern United States, and surveyed genetic variation through mitochondrial and nuclear sequences in 13 localities across three states. Our results indicate that the invader genotype, already identified, has also colonized the continental United States (US), supporting the hypothesis of a general-purpose genotype capable of coping with adverse conditions and enlarging its geographical range. Parthenogenesis, and its associated lack of recombination, may help in maintaining a general-purpose genotype that facilitates the colonization of distant, unsuitable areas. However, demographic advantages linked to parthenogenesis as the sole mode of reproduction are also possible.

Food resource access can mediate establishment success in invasive species, and generalist herbivorous insects are thought to rely on mechanisms of transcriptional plasticity to respond to dietary variation. While asexually reproducing invasives typically have low genetic variation, the twofold reproductive capacity of asexual organisms is a marked advantage for colonization. We studied host-related transcriptional acclimation in parthenogenetic, invasive, and polyphagous weevils: Naupactus cervinus and N. leucoloma. We analyzed patterns of gene expression in three gene categories that can mediate weevil-host plant interactions through identification of suitable host plants, short-term acclimation to host plant defenses, and long-term adaptation to host plant defenses and their pathogens. This approach employed comparative transcriptomic methods to investigate differentially expressed host detection, detoxification, immune defense genes, and pathway-level gene set enrichment. Our results show that weevil gene expression responses can be host plant-specific, and that elements of that response can be maintained in the offspring. Some host plant groups, such as legumes, appear to be more taxing as they elicit a complex gene expression response which is both strong in intensity and specific in identity. However, the weevil response to taxing host plants shares many differentially expressed genes with other stressful situations, such as host plant cultivation conditions and transition to novel host, suggesting that there is an evolutionarily favorable shared gene expression regime for responding to different types of stressful situations. Modulating gene expression in the absence of other avenues for phenotypic adaptation may be an important mechanism of successful colonization for these introduced insects.

Effective interpretation of genome function and genetic variation requires a shift from epigenetic mapping of cis-regulatory elements (CREs) to characterization of endogenous function. We developed hybridization chain reaction fluorescence in situ hybridization coupled with flow cytometry (HCR–FlowFISH), a broadly applicable approach to characterize CRISPR-perturbed CREs via accurate quantification of native transcripts, alongside CRISPR activity screen analysis (CASA), a hierarchical Bayesian model to quantify CRE activity. Across >325,000 perturbations, we provide evidence that CREs can regulate multiple genes, skip over the nearest gene and display activating and/or silencing effects. At the cholesterol-level-associated FADS locus, we combine endogenous screens with reporter assays to exhaustively characterize multiple genome-wide association signals, functionally nominate causal variants and, importantly, identify their target genes.