Andrew J. Knights

The genome of the social amoeba Dictyostelium discoideum

The social amoebae are exceptional in their ability to alternate between unicellular and multicellular forms. Here we describe the genome of the best-studied member of this group, Dictyostelium discoideum. The gene-dense chromosomes of this organism encode approximately 12,500 predicted proteins, a high proportion of which have long, repetitive amino acid tracts. There are many genes for polyketide synthases and ABC transporters, suggesting an extensive secondary metabolism for producing and exporting small molecules. The genome is rich in complex repeats, one class of which is clustered and may serve as centromeres. Partial copies of the extrachromosomal ribosomal DNA (rDNA) element are found at the ends of each chromosome, suggesting a novel telomere structure and the use of a common mechanism to maintain both the rDNA and chromosomal termini. A proteome-based phylogeny shows that the amoebozoa diverged from the animal–fungal lineage after the plant–animal split, but Dictyostelium seems to have retained more of the diversity of the ancestral genome than have plants, animals or fungi.

Sequence of Plasmodium falciparum chromosomes 1, 3–9 and 13

Since the sequencing of the first two chromosomes of the malaria parasite, Plasmodium falciparum, there has been a concerted effort to sequence and assemble the entire genome of this organism. Here we report the sequence of chromosomes 1, 3–9 and 13 of P. falciparum clone 3D7—these chromosomes account for approximately 55% of the total genome. We describe the methods used to map, sequence and annotate these chromosomes. By comparing our assemblies with the optical map, we indicate the completeness of the resulting sequence. During annotation, we assign Gene Ontology terms to the predicted gene products, and observe clustering of some malaria-specific terms to specific chromosomes. We identify a highly conserved sequence element found in the intergenic region of internal var genes that is not associated with their telomeric counterparts.

Fine-mapping cellular QTLs with RASQUAL and ATAC-seq

When cellular traits are measured using high-throughput DNA sequencing, quantitative trait loci (QTLs) manifest as fragment count differences between individuals and allelic differences within individuals. We present RASQUAL (Robust Allele-Specific Quantitation and Quality Control), a new statistical approach for association mapping that models genetic effects and accounts for biases in sequencing data using a single, probabilistic framework. RASQUAL substantially improves fine-mapping accuracy and sensitivity relative to existing methods in RNA-seq, DNase-seq and ChIP-seq data. We illustrate how RASQUAL can be used to maximize association detection by generating the first map of chromatin accessibility QTLs (caQTLs) in a European population using ATAC-seq. Despite a modest sample size, we identified 2,707 independent caQTLs (at a false discovery rate of 10%) and demonstrated how RASQUAL and ATAC-seq can provide powerful information for fine-mapping gene-regulatory variants and for linking distal regulatory elements with gene promoters. Our results highlight how combining between-individual and allele-specific genetic signals improves the functional interpretation of noncoding variation.

Molecular and functional variation in iPSC-derived sensory neurons

Induced pluripotent stem cells (iPSCs), and cells derived from them, have become key tools for modeling biological processes, particularly in cell types that are difficult to obtain from living donors. Here we present a map of regulatory variants in iPSC-derived neurons, based on 123 differentiations of iPSCs to a sensory neuronal fate. Gene expression was more variable across cultures than in primary dorsal root ganglion, particularly for genes related to nervous system development. Using single-cell RNA-sequencing, we found that the number of neuronal versus contaminating cells was influenced by iPSC culture conditions before differentiation. Despite high differentiation-induced variability, our allele-specific method detected thousands of quantitative trait loci (QTLs) that influenced gene expression, chromatin accessibility, and RNA splicing. On the basis of these detected QTLs, we estimate that recall-by-genotype studies that use iPSC-derived cells will require cells from at least 20–80 individuals to detect the effects of regulatory variants with moderately large effect sizes.This study identifies regulatory variants in sensory neurons derived from induced pluripotent stem cells. Despite differentiation-induced variability, an allele-specific method allowed detection of loci influencing gene expression, chromatin accessibility and RNA splicing.

Shared genetic effects on chromatin and gene expression reveal widespread enhancer priming in immune response

Noncoding regulatory variants play an important role in the genetics of complex traits. Although quantitative trait locus (QTL) mapping is a powerful approach to identify these variants, many genetic effects may remain unobserved when cells are sampled in only one of a large number of possible environments. Using a novel induced pluripotent stem cell-derived system, we mapped QTLs regulating chromatin accessibility and gene expression in macrophages in four conditions mimicking the interplay between interferon-gamma response and Salmonella infection. We found that approximately 50% of condition-specific effects on gene expression altered chromatin accessibility prior to stimulation. Furthermore, 6% of the chromatin accessibility QTLs regulated multiple neighbouring regions and these interactions were modulated by stimulation, occasionally producing condition-specific changes in gene expression. Profiling additional states also doubled the number of expression QTLs that could be confidently colocalised with disease associations. Thus, a substantial fraction of disease-associated variants may affect ‘primed’ regulatory elements in naive cells.Noncoding regulatory variants are often highly context-specific, modulating gene expression in a small subset of possible cellular states. Although these genetic effects are likely to play important roles in disease, the molecular mechanisms underlying context-specificity are not well understood. Here, we identify shared quantitative trait loci (QTLs) for chromatin accessibility and gene expression (eQTLs) and show that a large fraction (~60%) of eQTLs that appear following macrophage immune stimulation alter chromatin accessibility in unstimulated cells, suggesting they perturb enhancer priming. We show that such variants are likely to influence the binding of cell type specific transcription factors (TFs), such as PU.1, which then indirectly alter the binding of stimulus-specific TFs, such as NF-κB or STAT2. Our results imply that, although chromatin accessibility assays are powerful for fine mapping causal noncoding variants, detecting their downstream impact on gene expression will be challenging, requiring profiling of large numbers of stimulated cellular states and timepoints.

Oct-2 forms a complex with Oct-1 on the iNOS promoter and represses transcription by interfering with recruitment of RNA PolII by Oct-1

Oct-1 (POU2f1) and Oct-2 (POU2f2) are members of the POU family of transcription factors. They recognize the same DNA sequence but fulfil distinct functions: Oct-1 is ubiquitous and regulates a variety of genes while Oct-2 is restricted to B-cells and neurones. Here we examine the interplay and regulatory mechanisms of these factors to control the inducible nitric oxide synthase (iNOS, NOS2). Using two breast cancer cell lines as a comparative model, we found that MCF-7 express iNOS upon cytokine stimulation while MDA-MB-231 do not. Oct-1 is present in both cell lines but MDA-MB-231also express high levels of Oct-2. Manipulation of Oct-2 expression in these cell lines demonstrates that it is directly responsible for the repression of iNOS in MDA-MB-231. In MCF-7 cells Oct-1 binds the iNOS promoter, recruits RNA PolII and triggers initiation of transcription. In MDA-MB-231 cells, both Oct-1 and Oct-2 bind the iNOS promoter, forming a higher-order complex which fails to recruit RNA PolII, and as a consequence iNOS transcription does not proceed. Unravelling the mechanisms of transcription factor activity is paramount to the understanding of gene expression patterns that determine cell behaviour.

Shared genetic effects on chromatin and gene expression indicate a role for enhancer priming in immune response

Regulatory variants are often context specific, modulating gene expression in a subset of possible cellular states. Although these genetic effects can play important roles in disease, the molecular mechanisms underlying context specificity are poorly understood. Here, we identified shared quantitative trait loci (QTLs) for chromatin accessibility and gene expression in human macrophages exposed to IFNγ, Salmonella and IFNγ plus Salmonella. We observed that ~60% of stimulus-specific expression QTLs with a detectable effect on chromatin altered the chromatin accessibility in naive cells, thus suggesting that they perturb enhancer priming. Such variants probably influence binding of cell-type-specific transcription factors, such as PU.1, which can then indirectly alter the binding of stimulus-specific transcription factors, such as NF-κB or STAT2. Thus, although chromatin accessibility assays are powerful for fine-mapping causal regulatory variants, detecting their downstream effects on gene expression will be challenging, requiring profiling of large numbers of stimulated cellular states and time points.Analysis of chromatin accessibility and expression quantitative trait loci in stimulated or naïve macrophages identifies loci that constitutively alter chromatin but affect expression only after stimulation, thus indicating an effect on enhancer priming.

UTX-mediated enhancer and chromatin remodeling suppresses myeloid leukemogenesis through noncatalytic inverse regulation of ETS and GATA programs

The histone H3 Lys27-specific demethylase UTX (or KDM6A) is targeted by loss-of-function mutations in multiple cancers. Here, we demonstrate that UTX suppresses myeloid leukemogenesis through noncatalytic functions, a property shared with its catalytically inactive Y-chromosome paralog, UTY (or KDM6C). In keeping with this, we demonstrate concomitant loss/mutation of KDM6A (UTX) and UTY in multiple human cancers. Mechanistically, global genomic profiling showed only minor changes in H3K27me3 but significant and bidirectional alterations in H3K27ac and chromatin accessibility; a predominant loss of H3K4me1 modifications; alterations in ETS and GATA-factor binding; and altered gene expression after Utx loss. By integrating proteomic and genomic analyses, we link these changes to UTX regulation of ATP-dependent chromatin remodeling, coordination of the COMPASS complex and enhanced pioneering activity of ETS factors during evolution to AML. Collectively, our findings identify a dual role for UTX in suppressing acute myeloid leukemia via repression of oncogenic ETS and upregulation of tumor-suppressive GATA programs.This study shows that UTX (KDM6A) suppresses myeloid leukemogenesis through noncatalytic functions. UTX loss leads to alterations in H3K27ac, H3K4me1 and chromatin accessibility, and in gene-regulatory programs mediated by ETS and GATA transcription factors.

High resolution genetic mapping of causal regulatory interactions in the human genome

Physical interaction of distal regulatory elements in three-dimensional space poses a significant challenge for studies of common disease, because noncoding risk variants may be substantial distances from the genes they regulate. Experimental methods to capture these interactions, such as chromosome conformation capture (CCC), usually cannot assign causal direction of effect between regulatory elements, an important component of disease fine-mapping. Here, we developed a statistical model that uses Mendelian Randomisation within a Bayesian hierarchical model framework, and applied it to a novel ATAC-seq data from 100 individuals mapping over 15,000 putatively causal interactions between distal regions of open chromatin. Strikingly, the majority (>60%) of interactions we detected were over distances of <20Kb, a range where CCC-based methods perform poorly. Because we can infer the direction of causal interactions, the model also significantly improves our ability to fine-map: when we applied it to an eQTL data set we reduced the number of variants in the 90% credible set size by half. We experimentally validate one of our associations using CRISPR engineering of the BLK/FAM167A locus, which is associated with risk for a range of autoimmune diseases and show that the causal variant is likely to be a non-coding insertion within a CTCF binding motif. Our study suggests that many regulatory variants will be challenging to map to their cognate genes using CCC-based techniques, but association genetics of chromatin state can provide a powerful complement to these approaches.

Corrigendum: Fine-mapping cellular QTLs with RASQUAL and ATAC-seq

Nat. Genet. 48, 206–213 (2016); published online 14 December 2015; corrected after print 8 February 2016 In the version of this article initially published, the accession code for the ATAC-seq data was omitted. These data have been deposited in the European Nucleotide Archive under accession ERP011141.