Single-cell genomics and regulatory networks for 388 human brains.

Single-cell genomics is a powerful tool for studying heterogeneous tissues such as the brain. Yet little is understood about how genetic variants influence cell-level gene expression. Addressing this, we uniformly processed single-nuclei, multiomics datasets into a resource comprising >2.8 million nuclei from the prefrontal cortex across 388 individuals. For 28 cell types, we assessed population-level variation in expression and chromatin across gene families and drug targets. We identified >550,000 cell type-specific regulatory elements and >1.4 million single-cell expression quantitative trait loci, which we used to build cell-type regulatory and cell-to-cell communication networks. These networks manifest cellular changes in aging and neuropsychiatric disorders. We further constructed an integrative model accurately imputing single-cell expression and simulating perturbations; the model prioritized ~250 disease-risk genes and drug targets with associated cell types.

Single-cell genomics and regulatory networks for 388 human brains.

Single-cell genomics is a powerful tool for studying heterogeneous tissues such as the brain. Yet, little is understood about how genetic variants influence cell-level gene expression. Addressing this, we uniformly processed single-nuclei, multi-omics datasets into a resource comprising >2.8M nuclei from the prefrontal cortex across 388 individuals. For 28 cell types, we assessed population-level variation in expression and chromatin across gene families and drug targets. We identified >550K cell-type-specific regulatory elements and >1.4M single-cell expression-quantitative-trait loci, which we used to build cell-type regulatory and cell-to-cell communication networks. These networks manifest cellular changes in aging and neuropsychiatric disorders. We further constructed an integrative model accurately imputing single-cell expression and simulating perturbations; the model prioritized ~250 disease-risk genes and drug targets with associated cell types.

Validation of Enhancer Regions in Primary Human Neural Progenitor Cells using Capture STARR-seq.

Genome-wide association studies (GWAS) and expression analyses implicate noncoding regulatory regions as harboring risk factors for psychiatric disease, but functional characterization of these regions remains limited. We performed capture STARR-sequencing of over 78,000 candidate regions to identify active enhancers in primary human neural progenitor cells (phNPCs). We selected candidate regions by integrating data from NPCs, prefrontal cortex, developmental timepoints, and GWAS. Over 8,000 regions demonstrated enhancer activity in the phNPCs, and we linked these regions to over 2,200 predicted target genes. These genes are involved in neuronal and psychiatric disease-associated pathways, including dopaminergic synapse, axon guidance, and schizophrenia. We functionally validated a subset of these enhancers using mutation STARR-sequencing and CRISPR deletions, demonstrating the effects of genetic variation on enhancer activity and enhancer deletion on gene expression. Overall, we identified thousands of highly active enhancers and functionally validated a subset of these enhancers, improving our understanding of regulatory networks underlying brain function and disease.

Synthetic enhancer design by in silico compensatory evolution reveals flexibility and constraint in cis-regulation.

BACKGROUND: Models that incorporate specific chemical mechanisms have been successful in describing the activity of Drosophila developmental enhancers as a function of underlying transcription factor binding motifs. Despite this, the minimum set of mechanisms required to reconstruct an enhancer from its constituent parts is not known. Synthetic biology offers the potential to test the sufficiency of known mechanisms to describe the activity of enhancers, as well as to uncover constraints on the number, order, and spacing of motifs. RESULTS: Using a functional model and in silico compensatory evolution, we generated putative synthetic even-skipped stripe 2 enhancers with varying degrees of similarity to the natural enhancer. These elements represent the evolutionary trajectories of the natural stripe 2 enhancer towards two synthetic enhancers designed ab initio. In the first trajectory, spatially regulated expression was maintained, even after more than a third of binding sites were lost. In the second, sequences with high similarity to the natural element did not drive expression, but a highly diverged sequence about half the length of the minimal stripe 2 enhancer drove ten times greater expression. Additionally, homotypic clusters of Zelda or Stat92E motifs, but not Bicoid, drove expression in developing embryos. CONCLUSIONS: Here, we present a functional model of gene regulation to test the degree to which the known transcription factors and their interactions explain the activity of the Drosophila even-skipped stripe 2 enhancer. Initial success in the first trajectory showed that the gene regulation model explains much of the function of the stripe 2 enhancer. Cases where expression deviated from prediction indicates that undescribed factors likely act to modulate expression. We also showed that activation driven Bicoid and Hunchback is highly sensitive to spatial organization of binding motifs. In contrast, Zelda and Stat92E drive expression from simple homotypic clusters, suggesting that activation driven by these factors is less constrained. Collectively, the 40 sequences generated in this work provides a powerful training set for building future models of gene regulation.

Diverse patterns of genomic targeting by transcriptional regulators in Drosophila melanogaster.

Annotation of regulatory elements and identification of the transcription-related factors (TRFs) targeting these elements are key steps in understanding how cells interpret their genetic blueprint and their environment during development, and how that process goes awry in the case of disease. One goal of the modENCODE (model organism ENCyclopedia of DNA Elements) Project is to survey a diverse sampling of TRFs, both DNA-binding and non-DNA-binding factors, to provide a framework for the subsequent study of the mechanisms by which transcriptional regulators target the genome. Here we provide an updated map of the Drosophila melanogaster regulatory genome based on the location of 84 TRFs at various stages of development. This regulatory map reveals a variety of genomic targeting patterns, including factors with strong preferences toward proximal promoter binding, factors that target intergenic and intronic DNA, and factors with distinct chromatin state preferences. The data also highlight the stringency of the Polycomb regulatory network, and show association of the Trithorax-like (Trl) protein with hotspots of DNA binding throughout development. Furthermore, the data identify more than 5800 instances in which TRFs target DNA regions with demonstrated enhancer activity. Regions of high TRF co-occupancy are more likely to be associated with open enhancers used across cell types, while lower TRF occupancy regions are associated with complex enhancers that are also regulated at the epigenetic level. Together these data serve as a resource for the research community in the continued effort to dissect transcriptional regulatory mechanisms directing Drosophila development.

A locus in Drosophila sechellia affecting tolerance of a host plant toxin.

Many insects feed on only one or a few types of host. These host specialists often evolve a preference for chemical cues emanating from their host and develop mechanisms for circumventing their host's defenses. Adaptations like these are central to evolutionary biology, yet our understanding of their genetics remains incomplete. Drosophila sechellia, an emerging model for the genetics of host specialization, is an island endemic that has adapted to chemical toxins present in the fruit of its host plant, Morinda citrifolia. Its sibling species, D. simulans, and many other Drosophila species do not tolerate these toxins and avoid the fruit. Earlier work found a region with a strong effect on tolerance to the major toxin, octanoic acid, on chromosome arm 3R. Using a novel assay, we narrowed this region to a small span near the centromere containing 18 genes, including three odorant binding proteins. It has been hypothesized that the evolution of host specialization is facilitated by genetic linkage between alleles contributing to host preference and alleles contributing to host usage, such as tolerance to secondary compounds. We tested this hypothesis by measuring the effect of this tolerance locus on host preference behavior. Our data were inconsistent with the linkage hypothesis, as flies bearing this tolerance region showed no increase in preference for media containing M. citrifolia toxins, which D. sechellia prefers. Thus, in contrast to some models for host preference, preference and tolerance are not tightly linked at this locus nor is increased tolerance per se sufficient to change preference. Our data are consistent with the previously proposed model that the evolution of D. sechellia as a M. citrifolia specialist occurred through a stepwise loss of aversion and gain of tolerance to M. citrifolia's toxins.

Transcriptional regulation of the heterocyst patterning gene patA from Anabaena sp. strain PCC 7120.

The filamentous cyanobacterium Anabaena sp. strain PCC 7120 forms a periodic pattern of nitrogen-fixing heterocysts when grown in the absence of combined nitrogen. PatA is necessary for proper patterning of heterocysts along filaments. In this study, apparent transcriptional start points (tsps) were identified at nucleotides -305, -614, and -645 relative to the translational start site (-305, -614, and -645 tsps). Transcriptional reporter fusions were used to show that transcription from the -305 tsp was induced in all cells of filaments in response to nitrogen deprivation, required hetR for induction, and increased in a patA mutant. Transcription from -614/-645 tsp reporter fusions was spatially regulated and occurred primarily in cells that would become heterocysts. Complementation of a patA mutant strain by alleles encoding substitutions in, or deletion of, the putative phosphoacceptor C-terminal domain indicates that the PATAN domain can function independently of the C-terminal domain of PatA. Localization of a ring of PatA-GFP at sites of cell division, as well as the formation of enlarged cells with altered cell morphology when patA was overexpressed, suggests that PatA may participate in cell division.

Proper control of genetic background with precise allele substitution: a comment on Coyne and Elwyn.

Coyne and Elwyn report that, using Drosophila stocks provided by us, they were unable to replicate our experiments measuring the effects of desaturase-2 on stress tolerance. In this note, we provide evidence that these authors did not properly control for the differences in genetic background between the lines. Their experiments are thus not meaningful replications of our previous ones. We discuss ways of circumventing this problem in studies involving induced mutations in single genes.

Adaptive loss of an old duplicated gene during incipient speciation.

To probe the role of natural selection in species origin, we performed a DNA polymorphism survey of the Drosophila melanogaster desaturase2 (ds2) locus. ds2 is responsible for a cuticular hydrocarbon difference between two behaviorally isolated races--Zimbabwe (Z) and Cosmopolitan (M). The ds2 allele prevalent in the Z populations is functional, while the allele from the M populations harbors a 16-bp deletion upstream of the gene which knocks out its expression. We find a signature of positive selection in the ds2 promoter, but not in the control gene, sas. This signature appears to be confined to the derived M population. We also find that the selection has been recent because the gene retains a signature of a selective sweep evidenced by the departure of Fay and Wu's H test from neutral expectation. We also find that ds2, as well as its duplicate pair ds1, has been maintained in the Drosophila genus for at least 40 Myr without any sign of adaptive change. Taken together with previous molecular genetic evidence, our results suggest that ds2 is one of the genes responsible for adaptive divergence of the Z and M races of D. melanogaster.

Ecological adaptation during incipient speciation revealed by precise gene replacement.

To understand the role of adaptation in speciation, one must characterize the ecologically relevant phenotypic effects of naturally occurring alleles at loci potentially causing reproductive isolation. The desaturase2 gene of Drosophila melanogaster is such a locus. Two geographically differentiated ds2 alleles underlie a pheromonal difference between the Zimbabwe and Cosmopolitan races. We used a site-directed gene replacement technique to introduce an allele of ds2 from the Zimbabwe population into Cosmopolitan flies. We show that the Cosmopolitan allele confers resistance to cold as well as susceptibility to starvation when the entire genetic background is otherwise identical. We conclude that ecological adaptation likely accompanies sexual isolation between the two behavioral races of D. melanogaster.